Modeling transmission dynamics of lyme disease: Multiple vectors, seasonality, and vector mobility.

Lyme disease is the most prevalent tick-borne disease in the United States, which humans acquire from an infected tick of the genus Ixodes (primarily Ixodes scapularis). While previous studies have provided useful insights into various aspects of Lyme disease, the tick's host preference in the presence of multiple hosts has not been considered in the existing models. In this study, we develop a transmission dynamics model that includes the interactions between the primary vectors involved: blacklegged ticks (I. scapularis), white-footed mice (Peromyscus leucopus), and white-tailed deer (Odocoileus virginianus). Our model shows that the presence of multiple vectors may have a significant impact on the dynamics and spread of Lyme disease. Based on our model, we also calculate the basic reproduction number, R 0 , a threshold value that predicts whether a disease exists or dies out. Subsequent extensions of the model consider seasonality of the tick's feeding period and mobility of deer between counties. Our results suggest that a longer tick peak feeding period results in a higher infection prevalence. Moreover, while the deer mobility may not be a primary factor for short-term emergence of Lyme disease epidemics, in the long-run it can significantly contribute to local infectiousness in neighboring counties, which eventually reach the endemic steady state.

Normal and pathological dynamics of platelets in humans.

We develop a mathematical model of platelet, megakaryocyte, and thrombopoietin dynamics in humans. We show that there is a single stationary solution that can undergo a Hopf bifurcation, and use this information to investigate both normal and pathological platelet production, specifically cyclic thrombocytopenia. Carefully estimating model parameters from laboratory and clinical data, we then argue that a subset of parameters are involved in the genesis of cyclic thrombocytopenia based on clinical information. We provide model fits to the existing data for both platelet counts and thrombopoietin levels by changing four parameters that have physiological correlates. Our results indicate that the primary change in cyclic thrombocytopenia is an interference with, or destruction of, the thrombopoietin receptor with secondary changes in other processes, including immune-mediated destruction of platelets and megakaryocyte deficiency and failure in platelet production. This study contributes to the understanding of the origin of cyclic thrombocytopenia as well as extending the modeling of thrombopoiesis.

An expanded Notch-Delta model exhibiting long-range patterning and incorporating MicroRNA regulation.

Notch-Delta signaling is a fundamental cell-cell communication mechanism that governs the differentiation of many cell types. Most existing mathematical models of Notch-Delta signaling are based on a feedback loop between Notch and Delta leading to lateral inhibition of neighboring cells. These models result in a checkerboard spatial pattern whereby adjacent cells express opposing levels of Notch and Delta, leading to alternate cell fates. However, a growing body of biological evidence suggests that Notch-Delta signaling produces other patterns that are not checkerboard, and therefore a new model is needed. Here, we present an expanded Notch-Delta model that builds upon previous models, adding a local Notch activity gradient, which affects long-range patterning, and the activity of a regulatory microRNA. This model is motivated by our experiments in the ascidian Ciona intestinalis showing that the peripheral sensory neurons, whose specification is in part regulated by the coordinate activity of Notch-Delta signaling and the microRNA miR-124, exhibit a sparse spatial pattern whereby consecutive neurons may be spaced over a dozen cells apart. We perform rigorous stability and bifurcation analyses, and demonstrate that our model is able to accurately explain and reproduce the neuronal pattern in Ciona. Using Monte Carlo simulations of our model along with miR-124 transgene over-expression assays, we demonstrate that the activity of miR-124 can be incorporated into the Notch decay rate parameter of our model. Finally, we motivate the general applicability of our model to Notch-Delta signaling in other animals by providing evidence that microRNAs regulate Notch-Delta signaling in analogous cell types in other organisms, and by discussing evidence in other organisms of sparse spatial patterns in tissues where Notch-Delta signaling is active.

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.

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 infant's 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.

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 Earth's 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.

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

BACKGROUND: Phages, 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. RESULTS: PHACCS 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. CONCLUSION: PHACCS 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.

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 10(4) 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.

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.

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.

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.