MetaMap: an atlas of metatranscriptomic reads in human disease-related RNA-seq data.

Background: With the advent of the age of big data in bioinformatics, large volumes of data and high-performance computing power enable researchers to perform re-analyses of publicly available datasets at an unprecedented scale. Ever more studies imply the microbiome in both normal human physiology and a wide range of diseases. RNA sequencing technology (RNA-seq) is commonly used to infer global eukaryotic gene expression patterns under defined conditions, including human disease-related contexts; however, its generic nature also enables the detection of microbial and viral transcripts. Findings: We developed a bioinformatic pipeline to screen existing human RNA-seq datasets for the presence of microbial and viral reads by re-inspecting the non-human-mapping read fraction. We validated this approach by recapitulating outcomes from six independent, controlled infection experiments of cell line models and compared them with an alternative metatranscriptomic mapping strategy. We then applied the pipeline to close to 150 terabytes of publicly available raw RNA-seq data from more than 17,000 samples from more than 400 studies relevant to human disease using state-of-the-art high-performance computing systems. The resulting data from this large-scale re-analysis are made available in the presented MetaMap resource. Conclusions: Our results demonstrate that common human RNA-seq data, including those archived in public repositories, might contain valuable information to correlate microbial and viral detection patterns with diverse diseases. The presented MetaMap database thus provides a rich resource for hypothesis generation toward the role of the microbiome in human disease. Additionally, codes to process new datasets and perform statistical analyses are made available.

An age-dependent model to analyse the evolutionary stability of bacterial quorum sensing.

Bacterial communication is enabled through the collective release and sensing of signalling molecules in a process called quorum sensing. Cooperative processes can easily be destabilized by the appearance of cheaters, who contribute little or nothing at all to the production of common goods. This especially applies for planktonic cultures. In this study, we analyse the dynamics of bacterial quorum sensing and its evolutionary stability under two levels of cooperation, namely signal and enzyme production. The model accounts for mutation rates and switches between planktonic and biofilm state of growth. We present a mathematical approach to model these dynamics using age-dependent colony models. We explore the conditions under which cooperation is stable and find that spatial structuring can lead to long-term scenarios such as coexistence or bistability, depending on the non-linear combination of different parameters like death rates and production costs.

Short term uptake and transport process for metformin in roots of Phragmites australis and Typha latifolia.

Metformin (MET) as an emerging contaminant has been detected in surface water and wastewater in numerous countries, due to insufficient retention in classical waste water treatment plants. In order to characterize the uptake of the compound during phytotreatment of waste water, a short term Pitman chamber experiment was carried out to assess the characteristics of MET uptake and transport by roots. Three different concentrations (0.5, 1.0 and 2.0 mmol L(-)(1)) were applied to cattail (Typha latifolia) and reed (Phragmites australis) roots which were used to investigate the uptake mechanism because they are frequently utilized in phytoremediation. In addition, quinidine was used as an inhibitor to assess the role of organic cation transporters (OCTs) in the uptake of MET by T. latifolia. The transport process of MET is different from carbamazepine (CBZ) and caffeine (CFN). In both T. latifolia and P. australis, the uptake processes were independent of initial concentrations. Quinidine, a known inhibitor of organic cation transporters, can significantly affect MET uptake by T. latifolia roots with inhibition ratios of 70-74%. Uptake into the root could be characterized by a linear model with R(2) values in the range of 0.881-0.999. Overall, the present study provides evidence that MET is taken up by plant roots and has the potential for subsequent translocation. OCTs could be one of the important pathways for MET uptake into the plant.

Model selection for microbial nutrient uptake using a cost-benefit approach.

We consider the uptake of various carbon sources by microorganisms based on four fundamental assumptions: (1) the uptake of nutrient follows a saturation characteristics (2) substrate processing has a benefit but comes at costs of maintaining the process chain (3) substrate uptake is controlled and (4) evolution optimized the control of substrate uptake. These assumptions result in relatively simple mathematical models. In case of two substrates, our main finding is the following: Depending on the overall topology of the metabolic pathway, three different behavioral patterns can be identified. (1) both substrates are consumed at a time, (2) one substrate is preferred and represses the uptake of the other (catabolite repression), or (3) a cell feeds exclusively on one or the other substrate, possibly leading to a population that splits in two sub-populations, each of them specialized on one substrate only. Batch-culture and retentostat data of toluene, benzoate, and acetate uptake by Geobacter metallireducens are used to demonstrate that the model structure is suited for a quantitative description of uptake dynamics.

Bet-hedging in stochastically switching environments.

We investigate the evolution of bet-hedging in a population that experiences a stochastically switching environment by means of adaptive dynamics. The aim is to extend known results to the situation at hand, and to deepen the understanding of the range of validity of these results. We find three different types of evolutionarily stable strategies (ESSs) depending on the frequency at which the environment changes: for a rapid change, a monomorphic phenotype adapted to the mean environment; for an intermediate range, a bimorphic bet-hedging phenotype; for slowly changing environments, a monomorphic phenotype adapted to the current environment. While the last result is only obtained by means of heuristic arguments and simulations, the first two results are based on the analysis of Lyapunov exponents for stochastically switching systems.

Stochastic modeling of Pseudomonas syringae growth in the phyllosphere.

Pseudomonas syringae is a gram-negative bacterium which lives on leaf surfaces. Its growth has been described using epifluorescence microscopy and image analysis; it was found to be growing in aggregates of a wide range of sizes. We develop a stochastic model to describe aggregate distribution and determine the mechanisms generating experimental observations. We found that a logistic birth-death model with migration (time-homogeneous Markov process) provides the best description of the observed data. We discuss how to analyze the joint distribution of the numbers of aggregates of different sizes at a given time and explore how to account for new aggregates being created, that is, the joint distribution of the family size statistics conditional on the total number of aggregates. We compute the first two moments. Through simulations we examine how the model's parameters affect the aggregate size distribution and successfully explain the quantitative experimental data available. Aggregation formation is thought to be the first step towards pathogenic behavior of this bacterium; understanding aggregate size distribution would prove useful to understand the switch from epiphytic to pathogenic behavior.

Oscillations of Hes7 caused by negative autoregulation and ubiquitination.

Loss of Hes7 function leads to irregular somite formation demonstrating that Hes7 is a crucial component of the segmentation clock during somitogenesis. Experiments revealed that not only the repressor functionality but also the half-life of the protein is crucial for oscillatory expression of Hes7 and regular somite formation. Numerical integration of a delay equation system supported this finding. However, in a recent paper it was shown that the number of binding sites is also decisive for damped or undamped oscillations. It was shown that for more than one binding site the Hill coefficient increases. This leads to a completely different behavior. The oscillations are undamped and thus the mathematical model can no longer explain the results observed in the experiments. In this paper we propose a more sophisticated model for the Hes7 oscillator. Since Hes7 is degraded by the ubiquitin-proteasome pathway we include Michaelis-Menten kinetics for the ubiquitination of Hes7. We identify the Michaelis-Menten constant as an additional model parameter for oscillatory behavior. By increasing the Michaelis-Menten constant we found damped oscillations even if the Hill coefficient is increased.

Diphtheria (D) and tetanus (T) antibody values in children with acute lymphoblastic leukaemia (ALL) after treatment according to Co-ALL 05/92.

BACKGROUND: Children and adolescents after acute lymphoblastic leukemia are at risk for a prolonged period of immunodeficiency. Normally within 6 to 9 months after the end of maintenance treatment an adequate immune recovery is present. Factors such as immunity against specific antigens prior to disease (applied baseline vaccination), intensity of treatment and age can play a role in the appearance of antibodies in serum. Diphtheria (D) and Tetanus (T) antibodies are known to appear within 3 to 6 months after end of treatment as a sign of immune recovery and the reinstatement of immunological memory. A number of different questions are of interest: What differences are seen in the antibodies to D and T in children of different ages after treatment with a standardized protocol? What is the influence of post-treatment revaccination with Diphtheria/Tetanus (D/T) and treatment group on the production of D/T antibodies? PATIENTS AND METHODS: Out of 142 children and adolescents until the age of 16, treated according to the Co-ALL 05/92 protocol, 59 patients were eligible for evaluation: 31 Low-Risk (LR)- and 28 High-Risk (HR) patients. Antibodies against Diphtheria (D) and Tetanus (T) were measured 3-12 months after the end of treatment and after revaccination in case of low antibody levels against D and/or T. In patients without adequate response after repeated revaccination the cellular immunity was examined with a skin test. RESULTS: After the end of treatment, children in the low-risk (LR)-group showed more frequently adequate antibody titres against D and T than children of the high-risk (HR)-group. Antibodies against T were present in 50% of all patients. After revaccination antibodies against T were found in nearly all patients whereas for D this is only the case in some children. Patients without sufficient antibody levels mainly showed an adequate cellular immunity. CONCLUSION: In children and adolescents with ALL after therapy antibody levels of D and T are dependent on treatment intensity. Revaccination leads to an adequate immunological answer against T in most patients , which is not the case for the diphtheria vaccination. Prospective multicenter trials starting together with the ALL-treatment should be able to gain more information about the behavior of antibody levels and the risk of infection from vaccine-preventable disease in immunocompromised patients and thus lead to standardized vaccination guidelines such as immunization with conjugate vaccines already during maintenance treatment.