Probabilistic Models to Describe the Dynamics of Migrating Microbial Communities

In all but the most sterile environments bacteria will reside in fluid being transported through conduits and some of these will attach and grow as biofilms on the conduit walls. The concentration and diversity of bacteria in the fluid at the point of delivery will be a mix of those when it entered the conduit and those that have become entrained into the flow due to seeding from biofilms. Examples include fluids through conduits such as drinking water pipe networks, endotracheal tubes, catheters and ventilation systems. Here we present two probabilistic models to describe changes in the composition of bulk fluid microbial communities as they are transported through a conduit whilst exposed to biofilm communities. The first (discrete) model simulates absolute numbers of individual cells, whereas the other (continuous) model simulates the relative abundance of taxa in the bulk fluid. The discrete model is founded on a birth-death process whereby the community changes one individual at a time and the numbers of cells in the system can vary. The continuous model is a stochastic differential equation derived from the discrete model and can also accommodate changes in the carrying capacity of the bulk fluid. These models provide a novel Lagrangian framework to investigate and predict the dynamics of migrating microbial communities. In this paper we compare the two models, discuss their merits, possible applications and present simulation results in the context of drinking water distribution systems. Our results provide novel insight into the effects of stochastic dynamics on the composition of non-stationary microbial communities that are exposed to biofilms and provides a new avenue for modelling microbial dynamics in systems where fluids are being transported

A metastable equilibrium model for the relative abundances of microbial phyla in a hot spring

Many studies link the compositions of microbial communities to their environments, but the energetics of organism-specific biomass synthesis as a function of geochemical variables has rarely been assessed. We describe a thermodynamic model that integrates geochemical and metagenomic data for biofilms sampled at five sites along a thermal and chemical gradient in the outflow channel of the hot spring known as ‘‘Bison Pool’’ in Yellowstone National Park. The relative abundances of major phyla in individual communities sampled along the outflow channel are modeled by computing metastable equilibrium among model proteins with amino acid compositions derived from metagenomic sequences. Geochemical conditions are represented by temperature and activities of basis species, including pH and oxidation-reduction potential quantified as the activity of dissolved hydrogen. By adjusting the activity of hydrogen, the model can be tuned to closely approximate the relative abundances of the phyla observed in the community profiles generated from BLAST assignments. The findings reveal an inverse relationship between the energy demand to form the proteins at equal thermodynamic activities and the abundance of phyla in the community.Although the metabolisms used by many members of these communities are driven by chemical disequilibria, the results support the possibility that higher-level patterns of chemotrophic microbial ecosystems are shaped by metastable equilibrium states that depend on both the composition of biomass and the environmental conditions

Robust estimation of microbial diversity in theory and in practice

Quantifying diversity is of central importance for the study of structure, function and evolution of microbial communities. The estimation of microbial diversity has received renewed attention with the advent of large-scale metagenomic studies. Here, we consider what the diversity observed in a sample tells us about the diversity of the community being sampled. First, we argue that one cannot reliably estimate the absolute and relative number of microbial species present in a community without making unsupported assumptions about species abundance distributions. The reason for this is that sample data do not contain information about the number of rare species in the tail of species abundance distributions. We illustrate the difficulty in comparing species richness estimates by applying Chao's estimator of species richness to a set of in silico communities: they are ranked incorrectly in the presence of large numbers of rare species. Next, we extend our analysis to a general family of diversity metrics ("Hill diversities"), and construct lower and upper estimates of diversity values consistent with the sample data. The theory generalizes Chao's estimator, which we retrieve as the lower estimate of species richness. We show that Shannon and Simpson diversity can be robustly estimated for the in silico communities. We analyze nine metagenomic data sets from a wide range of environments, and show that our findings are relevant for empirically-sampled communities. Hence, we recommend the use of Shannon and Simpson diversity rather than species richness in efforts to quantify and compare microbial diversity.Comment: To be published in The ISME Journal. Main text: 16 pages, 5 figures. Supplement: 16 pages, 4 figure

The treatment and outcomes of early-stage epithelial ovarian cancer: have we made any progress?

The objective of this study is to determine the progress and trends in the treatment and survival of women with early-stage (I–II) epithelial ovarian cancer. Data were obtained from the SEER database between 1988 and 2001. Kaplan–Meier and Cox regressions methods were employed for statistical analyses. Of the 8372 patients, the median age was 57 years (range: 12–99 years). A total of 6152 patients (73.4%) presented with stage I and 2220 (26.5%) with stage II disease. Over the periods 1988–1992, 1993–1997, and 1998–2001, 3-year disease-specific survivals increased from 86.1 to 87.2 to 88.8% (P=0.076). The number of patients that underwent lymphadenectomy has increased significantly from 26.2 to 38.7 to 54.2% over the study period (P<0.001). Of those patients who underwent staging procedures with lymphadenectomy, there was no improvement in survival over the three study periods (from 93.2 to 93.5 to 93.1%; P=0.978). On multivariate analysis, younger age, nonclear cell histology, earlier stage, lower grade, surgery, and lymphadenectomy were significant independent prognostic factors for improved survival. After adjusting for surgical staging with lymphadenectomy, the year of diagnosis was no longer an important prognostic factor. In conclusion, the use of lymphadenectomy during surgery for early-stage ovarian cancer has doubled over the last 14 years. The marginal improvement in survival demonstrated over time is potentially attributed to the increased use of staging procedures with lymphadenectomy

The potential to encode sex, age, and individual identity in the alarm calls of three species of Marmotinae

In addition to encoding referential information and information about the sender’s motivation, mammalian alarm calls may encode information about other attributes of the sender, providing the potential for recognition among kin, mates, and neighbors. Here, we examined 96 speckled ground squirrels (Spermophilus suslicus), 100 yellow ground squirrels (Spermophilus fulvus) and 85 yellow-bellied marmots (Marmota flaviventris) to determine whether their alarm calls differed between species in their ability to encode information about the caller’s sex, age, and identity. Alarm calls were elicited by approaching individually identified animals in live-traps. We assume this experimental design modeled a naturally occurring predatory event, when receivers should acquire information about attributes of a caller from a single bout of alarm calls. In each species, variation that allows identification of the caller’s identity was greater than variation allowing identification of age or sex. We discuss these results in relation to each species’ biology and sociality

Measurement of the inclusive and dijet cross-sections of b-jets in pp collisions at sqrt(s) = 7 TeV with the ATLAS detector

The inclusive and dijet production cross-sections have been measured for jets containing b-hadrons (b-jets) in proton-proton collisions at a centre-of-mass energy of sqrt(s) = 7 TeV, using the ATLAS detector at the LHC. The measurements use data corresponding to an integrated luminosity of 34 pb^-1. The b-jets are identified using either a lifetime-based method, where secondary decay vertices of b-hadrons in jets are reconstructed using information from the tracking detectors, or a muon-based method where the presence of a muon is used to identify semileptonic decays of b-hadrons inside jets. The inclusive b-jet cross-section is measured as a function of transverse momentum in the range 20 < pT < 400 GeV and rapidity in the range |y| < 2.1. The bbbar-dijet cross-section is measured as a function of the dijet invariant mass in the range 110 < m_jj < 760 GeV, the azimuthal angle difference between the two jets and the angular variable chi in two dijet mass regions. The results are compared with next-to-leading-order QCD predictions. Good agreement is observed between the measured cross-sections and the predictions obtained using POWHEG + Pythia. MC@NLO + Herwig shows good agreement with the measured bbbar-dijet cross-section. However, it does not reproduce the measured inclusive cross-section well, particularly for central b-jets with large transverse momenta.Comment: 10 pages plus author list (21 pages total), 8 figures, 1 table, final version published in European Physical Journal

Alternative Ii-independent antigen-processing pathway in leukemic blasts involves TAP-dependent peptide loading of HLA class II complexes

During HLA class II synthesis in antigen-presenting cells, the invariant chain (Ii) not only stabilizes HLA class II complexes in the endoplasmic reticulum, but also mediates their transport to specialized lysosomal antigen-loading compartments termed MIICs. This study explores an alternative HLA class II presentation pathway in leukemic blasts that involves proteasome and transporter associated with antigen processing (TAP)-dependent peptide loading. Although HLA-DR did associate with Ii, Ii silencing in the human class II-associated invariant chain peptide (CLIP)-negative KG-1 myeloid leukemic cell line did not affect total and plasma membrane expression levels of HLA-DR, as determined by western blotting and flow cytometry. Since HLA-DR expression does require peptide binding, we examined the role of endogenous antigen-processing machinery in HLA-DR presentation by CLIP− leukemic blasts. The suppression of proteasome and TAP function using various inhibitors resulted in decreased HLA-DR levels in both CLIP− KG-1 and ME-1 blasts. Simultaneous inhibition of TAP and Ii completely down-modulated the expression of HLA-DR, demonstrating that together these molecules form the key mediators of HLA class II antigen presentation in leukemic blasts. By the use of a proteasome- and TAP-dependent pathway for HLA class II antigen presentation, CLIP− leukemic blasts might be able to present a broad range of endogenous leukemia-associated peptides via HLA class II to activate leukemia-specific CD4+ T cells

Contrasting Microbial Community Assembly Hypotheses: A Reconciling Tale from the Río Tinto

The Río Tinto (RT) is distinguished from other acid mine drainage systems by its natural and ancient origins. Microbial life from all three domains flourishes in this ecosystem, but bacteria dominate metabolic processes that perpetuate environmental extremes. While the patchy geochemistry of the RT likely influences the dynamics of bacterial populations, demonstrating which environmental variables shape microbial diversity and unveiling the mechanisms underlying observed patterns, remain major challenges in microbial ecology whose answers rely upon detailed assessments of community structures coupled with fine-scale measurements of physico-chemical parameters.By using high-throughput environmental tag sequencing we achieved saturation of richness estimators for the first time in the RT. We found that environmental factors dictate the distribution of the most abundant taxa in this system, but stochastic niche differentiation processes, such as mutation and dispersal, also contribute to observed diversity patterns.We predict that studies providing clues to the evolutionary and ecological processes underlying microbial distributions will reconcile the ongoing debate between the Baas Becking vs. Hubbell community assembly hypotheses