Joseph M. Mahaffy

The Marine Viromes of Four Oceanic Regions

Viruses are the most common biological entities in the marine environment. There has not been a global survey of these viruses, and consequently, it is not known what types of viruses are in Earths oceans or how they are distributed. Metagenomic analyses of 184 viral assemblages collected over a decade and representing 68 sites in four major oceanic regions showed that most of the viral sequences were not similar to those in the current databases. There was a distinct “marine-ness” quality to the viral assemblages. Global diversity was very high, presumably several hundred thousand of species, and regional richness varied on a North-South latitudinal gradient. The marine regions had different assemblages of viruses. Cyanophages and a newly discovered clade of single-stranded DNA phages dominated the Sargasso Sea sample, whereas prophage-like sequences were most common in the Arctic. However most viral species were found to be widespread. With a majority of shared species between oceanic regions, most of the differences between viral assemblages seemed to be explained by variation in the occurrence of the most common viral species and not by exclusion of different viral genomes. These results support the idea that viruses are widely dispersed and that local environmental conditions enrich for certain viral types through selective pressure.

Genomic analysis of uncultured marine viral communities.

Viruses are the most common biological entities in the oceans by an order of magnitude. However, very little is known about their diversity. Here we report a genomic analysis of two uncultured marine viral communities. Over 65% of the sequences were not significantly similar to previously reported sequences, suggesting that much of the diversity is previously uncharacterized. The most common significant hits among the known sequences were to viruses. The viral hits included sequences from all of the major families of dsDNA tailed phages, as well as some algal viruses. Several independent mathematical models based on the observed number of contigs predicted that the most abundant viral genome comprised 2–3% of the total population in both communities, which was estimated to contain between 374 and 7,114 viral types. Overall, diversity of the viral communities was extremely high. The results also showed that it would be possible to sequence the entire genome of an uncultured marine viral community.

Metagenomic Analyses of an Uncultured Viral Community from Human Feces

Here we present the first metagenomic analyses of an uncultured viral community from human feces, using partial shotgun sequencing. Most of the sequences were unrelated to anything previously reported. The recognizable viruses were mostly siphophages, and the community contained an estimated 1,200 viral genotypes.

Viral and microbial community dynamics in four aquatic environments.

The species composition and metabolic potential of microbial and viral communities are predictable and stable for most ecosystems. This apparent stability contradicts theoretical models as well as the viral–microbial dynamics observed in simple ecosystems, both of which show Kill-the-Winner behavior causing cycling of the dominant taxa. Microbial and viral metagenomes were obtained from four human-controlled aquatic environments at various time points separated by one day to >1 year. These environments were maintained within narrow geochemical bounds and had characteristic species composition and metabolic potentials at all time points. However, underlying this stability were rapid changes at the fine-grained level of viral genotypes and microbial strains. These results suggest a model wherein functionally redundant microbial and viral taxa are cycling at the level of viral genotypes and virus-sensitive microbial strains. Microbial taxa, viral taxa, and metabolic function persist over time in stable ecosystems and both communities fluctuate in a Kill-the-Winner manner at the level of viral genotypes and microbial strains.

Diversity and population structure of a near–shore marine–sediment viral community

Viruses, most of which are phage, are extremely abundant in marine sediments, yet almost nothing is known about their identity or diversity. We present the metagenomic analysis of an uncultured near–shore marine–sediment viral community. Three–quarters of the sequences in the sample were not related to anything previously reported. Among the sequences that could be identified, the majority belonged to double–stranded DNA phage. Temperate phage were more common than lytic phage, suggesting that lysogeny may be an important lifestyle for sediment viruses. Comparisons between the sediment sample and previously sequenced seawater viral communities showed that certain phage phylogenetic groups were abundant in all marine viral communities, while other phage groups were under–represented or absent. This ‘marineness’ suggests that marine phage are derived from a common set of ancestors. Several independent mathematical models, based on the distribution of overlapping shotgun sequence fragments from the library, were used to show that the diversity of the viral community was extremely high, with at least 104 viral genotypes per kilogram of sediment and a Shannon index greater than 9 nats. Based on these observations we propose that marine–sediment viral communities are one of the largest unexplored reservoirs of sequence space on the planet.

Viral diversity and dynamics in an infant gut

Metagenomic sequencing of DNA viruses from the feces of a healthy week-old infant revealed a viral community with extremely low diversity. The identifiable sequences were dominated by phages, which likely influence the diversity and abundance of co-occurring microbes. The most abundant fecal viral sequences did not originate from breast milk or formula, suggesting a non-dietary initial source of viruses. Certain sequences were stable in the infants gut over the first 3 months of life, but microarray experiments demonstrated that the overall viral community composition changed dramatically between 1 and 2 weeks of age.

Age-structured and two-delay models for erythropoiesis

An age-structured model is developed for erythropoiesis and is reduced to a system of threshold-type differential delay equations using the method of characteristics. Under certain assumptions, this model can be reduced to a system of delay differential equations with two delays. The parameters in the system are estimated from experimental data, and the model is simulated for a normal human subject following a loss of blood. The characteristic equation of the two-delay equation is analyzed and shown to exhibit Hopf bifurcations when the destruction rate of erythrocytes is increased. A numerical study for a rabbit with autoimmune hemolytic anemia is performed and compared with experimental data.

PHACCS, an online tool for estimating the structure and diversity of uncultured viral communities using metagenomic information.

BackgroundPhages, viruses that infect prokaryotes, are the most abundant microbes in the world. A major limitation to studying these viruses is the difficulty of cultivating the appropriate prokaryotic hosts. One way around this limitation is to directly clone and sequence shotgun libraries of uncultured viral communities (i.e., metagenomic analyses). PHACCS http://phage.sdsu.edu/phaccs, Phage Communities from Contig Spectrum, is an online bioinformatic tool to assess the biodiversity of uncultured viral communities. PHACCS uses the contig spectrum from shotgun DNA sequence assemblies to mathematically model the structure of viral communities and make predictions about diversity.ResultsPHACCS builds models of possible community structure using a modified Lander-Waterman algorithm to predict the underlying contig spectrum. PHACCS finds the most appropriate structure model by optimizing the model parameters until the predicted contig spectrum is as close as possible to the experimental one. This model is the basis for making estimates of uncultured viral community richness, evenness, diversity index and abundance of the most abundant genotype.ConclusionPHACCS analysis of four different environmental phage communities suggests that the power law is an important rank-abundance form to describe uncultured viral community structure. The estimates support the fact that the four phage communities were extremely diverse and that phage community biodiversity and structure may be correlated with that of their hosts.

A model for the initiation of replication in Escherichia coli

The role of the protein DnaA as the principal control of replication initiation is investigated by a mathematical model. Data showing that DnaA is growth rate regulated suggest that its concentration alone is not the only factor determining the timing of initiation. A mathematical model with stochastic and deterministic components is constructed from known experimental evidence and subdivides the total pool of DnaA protein into four forms. The active form, DnaA.ATP, can be bound to the origin of replication, oriC, where it is assumed that a critical level of these bound molecules is needed to initiate replication. The active form can also exist in a reserve pool bound to the chromosome or a free pool in the cytoplasm. Finally, a large inactive pool of DnaA protein completes the state variables and provides an explanation for how the DnaA.ATP form could be the principal controlling element in the timing of initiation. The fact that DnaA protein is an autorepressor is used to derive its synthesis rate. The model studies a single exponentially growing cell through a series of cell divisions. Computer simulations are performed, and the results compare favorably to data for different cell cycle times. The model shows synchrony of initiation events in agreement with experimental results.

Escherichia coli Cells with Increased Levels of DnaA and Deficient in Recombinational Repair Have Decreased Viability

The dnaA operon of Escherichia coli contains the genes dnaA, dnaN, and recF encoding DnaA, beta clamp of DNA polymerase III holoenzyme, and RecF. When the DnaA concentration is raised, an increase in the number of DNA replication initiation events but a reduction in replication fork velocity occurs. Because DnaA is autoregulated, these results might be due to the inhibition of dnaN and recF expression. To test this, we examined the effects of increasing the intracellular concentrations of DnaA, beta clamp, and RecF, together and separately, on initiation, the rate of fork movement, and cell viability. The increased expression of one or more of the dnaA operon proteins had detrimental effects on the cell, except in the case of RecF expression. A shorter C period was not observed with increased expression of the beta clamp; in fact, many chromosomes did not complete replication in runout experiments. Increased expression of DnaA alone resulted in stalled replication forks, filamentation, and a decrease in viability. When the three proteins of the dnaA operon were simultaneously overexpressed, highly filamentous cells were observed (>50 micro m) with extremely low viability and, in runout experiments, most chromosomes had not completed replication. The possibility that recombinational repair was responsible for the survival of cells overexpressing DnaA was tested by using mutants in different recombinational repair pathways. The absence of RecA, RecB, RecC, or the proteins in the RuvABC complex caused an additional approximately 100-fold drop in viability in cells with increased levels of DnaA, indicating a requirement for recombinational repair in these cells.