BSim: an agent-based tool for modeling bacterial populations in systems and synthetic biology.

Open Access ArticleLarge-scale collective behaviors such as synchronization and coordination spontaneously arise in many bacterial populations. With systems biology attempting to understand these phenomena, and synthetic biology opening up the possibility of engineering them for our own benefit, there is growing interest in how bacterial populations are best modeled. Here we introduce BSim, a highly flexible agent-based computational tool for analyzing the relationships between single-cell dynamics and population level features. BSim includes reference implementations of many bacterial traits to enable the quick development of new models partially built from existing ones. Unlike existing modeling tools, BSim fully considers spatial aspects of a model allowing for the description of intricate micro-scale structures, enabling the modeling of bacterial behavior in more realistic three-dimensional, complex environments. The new opportunities that BSim opens are illustrated through several diverse examples covering: spatial multicellular computing, modeling complex environments, population dynamics of the lac operon, and the synchronization of genetic oscillators. BSim is open source software that is freely available from http://bsim-bccs.sf.net and distributed under the Open Source Initiative (OSI) recognized MIT license. Developer documentation and a wide range of example simulations are also available from the website. BSim requires Java version 1.6 or higher.Engineering and Physical Sciences Research Council (EPSRC)Biotechnology and Biological Sciences Research Council (BBSRC

Airway resistance at maximum inhalation as a marker of asthma and airway hyperresponsiveness

<p>Abstract</p> <p>Background</p> <p>Asthmatics exhibit reduced airway dilation at maximal inspiration, likely due to structural differences in airway walls and/or functional differences in airway smooth muscle, factors that may also increase airway responsiveness to bronchoconstricting stimuli. The goal of this study was to test the hypothesis that the minimal airway resistance achievable during a maximal inspiration (R_min) is abnormally elevated in subjects with airway hyperresponsiveness.</p> <p>Methods</p> <p>The R_minwas measured in 34 nonasthmatic and 35 asthmatic subjects using forced oscillations at 8 Hz. R_minand spirometric indices were measured before and after bronchodilation (albuterol) and bronchoconstriction (methacholine). A preliminary study of 84 healthy subjects first established height dependence of baseline R_minvalues.</p> <p>Results</p> <p>Asthmatics had a higher baseline R_min% predicted than nonasthmatic subjects (134 ± 33 vs. 109 ± 19 % predicted, p = 0.0004). Sensitivity-specificity analysis using receiver operating characteristic curves indicated that baseline R_minwas able to identify subjects with airway hyperresponsiveness (PC₂₀< 16 mg/mL) better than most spirometric indices (Area under curve = 0.85, 0.78, and 0.87 for R_min% predicted, FEV₁% predicted, and FEF_25-75% predicted, respectively). Also, 80% of the subjects with baseline R_min< 100% predicted did not have airway hyperresponsiveness while 100% of subjects with R_min> 145% predicted had hyperresponsive airways, regardless of clinical classification as asthmatic or nonasthmatic.</p> <p>Conclusions</p> <p>These findings suggest that baseline R_min, a measurement that is easier to perform than spirometry, performs as well as or better than standard spirometric indices in distinguishing subjects with airway hyperresponsiveness from those without hyperresponsive airways. The relationship of baseline R_minto asthma and airway hyperresponsiveness likely reflects a causal relation between conditions that stiffen airway walls and hyperresponsiveness. In conjunction with symptom history, R_mincould provide a clinically useful tool for assessing asthma and monitoring response to treatment.</p

Competition-based model of pheromone component ratio detection in the moth

For some moth species, especially those closely interrelated and sympatric, recognizing a specific pheromone component concentration ratio is essential for males to successfully locate conspecific females. We propose and determine the properties of a minimalist competition-based feed-forward neuronal model capable of detecting a certain ratio of pheromone components independently of overall concentration. This model represents an elementary recognition unit for the ratio of binary mixtures which we propose is entirely contained in the macroglomerular complex (MGC) of the male moth. A set of such units, along with projection neurons (PNs), can provide the input to higher brain centres. We found that (1) accuracy is mainly achieved by maintaining a certain ratio of connection strengths between olfactory receptor neurons (ORN) and local neurons (LN), much less by properties of the interconnections between the competing LNs proper. An exception to this rule is that it is beneficial if connections between generalist LNs (i.e. excited by either pheromone component) and specialist LNs (i.e. excited by one component only) have the same strength as the reciprocal specialist to generalist connections. (2) successful ratio recognition is achieved using latency-to-first-spike in the LN populations which, in contrast to expectations with a population rate code, leads to a broadening of responses for higher overall concentrations consistent with experimental observations. (3) when longer durations of the competition between LNs were observed it did not lead to higher recognition accuracy

A disk of dust and molecular gas around a high-mass protostar

The processes leading to the birth of low-mass stars such as our Sun have been well studied, but the formation of high-mass (> 8 x Sun's mass) stars has heretofore remained poorly understood. Recent observational studies suggest that high-mass stars may form in essentially the same way as low-mass stars, namely via an accretion process, instead of via merging of several low-mass (< 8 Msun) stars. However, there is as yet no conclusive evidence. Here, we report the discovery of a flattened disk-like structure observed at submillimeter wavelengths, centered on a massive 15 Msun protostar in the Cepheus-A region. The disk, with a radius of about 330 astronomical units (AU) and a mass of 1 to 8 Msun, is detected in dust continuum as well as in molecular line emission. Its perpendicular orientation to, and spatial coincidence with the central embedded powerful bipolar radio jet, provides the best evidence yet that massive stars form via disk accretion in direct analogy to the formation of low-mass stars

Differences in genotype and virulence among four multidrug-resistant <i>Streptococcus pneumoniae</i> isolates belonging to the PMEN1 clone

We report on the comparative genomics and characterization of the virulence phenotypes of four &lt;i&gt;S. pneumoniae&lt;/i&gt; strains that belong to the multidrug resistant clone PMEN1 (Spain&lt;sup&gt;23F&lt;/sup&gt; ST81). Strains SV35-T23 and SV36-T3 were recovered in 1996 from the nasopharynx of patients at an AIDS hospice in New York. Strain SV36-T3 expressed capsule type 3 which is unusual for this clone and represents the product of an in vivo capsular switch event. A third PMEN1 isolate - PN4595-T23 - was recovered in 1996 from the nasopharynx of a child attending day care in Portugal, and a fourth strain - ATCC700669 - was originally isolated from a patient with pneumococcal disease in Spain in 1984. We compared the genomes among four PMEN1 strains and 47 previously sequenced pneumococcal isolates for gene possession differences and allelic variations within core genes. In contrast to the 47 strains - representing a variety of clonal types - the four PMEN1 strains grouped closely together, demonstrating high genomic conservation within this lineage relative to the rest of the species. In the four PMEN1 strains allelic and gene possession differences were clustered into 18 genomic regions including the capsule, the blp bacteriocins, erythromycin resistance, the MM1-2008 prophage and multiple cell wall anchored proteins. In spite of their genomic similarity, the high resolution chinchilla model was able to detect variations in virulence properties of the PMEN1 strains highlighting how small genic or allelic variation can lead to significant changes in pathogenicity and making this set of strains ideal for the identification of novel virulence determinant

Association of the CpG Methylation Pattern of the Proximal Insulin Gene Promoter with Type 1 Diabetes

The insulin (INS) region is the second most important locus associated with Type 1 Diabetes (T1D). The study of the DNA methylation pattern of the 7 CpGs proximal to the TSS in the INS gene promoter revealed that T1D patients have a lower level of methylation of CpG -19, -135 and -234 (p = 2.10−16) and a higher methylation of CpG -180 than controls, while methylation was comparable for CpG -69, -102, -206. The magnitude of the hypomethylation relative to a control population was 8–15% of the corresponding levels in controls and was correlated in CpGs -19 and -135 (r = 0.77) and CpG -135 and -234 (r = 0.65). 70/485 (14%) of T1D patients had a simultaneous decrease in methylation of CpG -19, -135, -234 versus none in 317 controls. CpG methylation did not correlate with glycated hemoglobin or with T1D duration. The methylation of CpG -69, -102, -180, -206, but not CpG -19, -135, -234 was strongly influenced by the cis-genotype at rs689, a SNP known to show a strong association with T1D. We hypothesize that part of this genetic association could in fact be mediated at the statistical and functional level by the underlying changes in neighboring CpG methylation. Our observation of a CpG-specific, locus-specific methylation pattern, although it can provide an epigenetic biomarker of a multifactorial disease, does not indicate whether the reported epigenetic pattern preexists or follows the establishment of T1D. To explore the effect of chronic hyperglycemia on CpG methylation, we studied non obese patients with type 2 diabetes (T2D) who were found to have decreased CpG-19 methylation versus age-matched controls, similar to T1D (p = 2.10−6) but increased CpG-234 methylation (p = 5.10−8), the opposite of T1D. The causality and natural history of the different epigenetic changes associated with T1D or T2D remain to be determined

Study of Transcriptional Effects in Cis at the IFIH1 Locus

Background: The Thr allele at the non-synonymous single-nucleotide polymorphism (nsSNP) Thr946Ala in the IFIH1 gene confers risk for Type 1 diabetes (T1D). The SNP is embedded in a 236 kb linkage disequilibrium (LD) block that includes four genes: IFIH1, GCA, FAP and KCNH7. The absence of common nsSNPs in the other genes makes the IFIH1 SNP the strongest functional candidate, but it could be merely a marker of association, due to LD with a variant regulating expression levels of IFIH1 or neighboring genes. Methodology/Principal Findings: We investigated the effect of the T1D-associated variation on mRNA transcript expression of these genes. Heterozygous mRNA from lymphoblastoid cell lines (LCLs), pancreas and thymus was examined by allelic expression imbalance, to detect effects in cis on mRNA expression. Using single-nucleotide primer extension, we found no difference between mRNA transcripts in 9 LCLs, 6 pancreas and 13 thymus samples, suggesting that GCA and FAP are not involved. On the other hand, KCNH7 was not expressed at a detectable level in all tissues examined. Moreover, the association of the Thr946Ala SNP with T1D is not due to modulation of IFIH1 expression in organs involved in the disease, pointing to the IFIH1 nsSNP as the causal variant. Conclusions/Significance: The mechanism of the association of the nsSNP with T1D remains to be determined, but does not involve mRNA modulation. It becomes necessary to study differential function of the IFIH1 protein alleles at Thr946Al

Erasing Sensorimotor Memories via PKMζ Inhibition

Sensorimotor cortex has a role in procedural learning. Previous studies suggested that this learning is subserved by long-term potentiation (LTP), which is in turn maintained by the persistently active kinase, protein kinase Mzeta (PKMζ). Whereas the role of PKMζ in animal models of declarative knowledge is established, its effect on procedural knowledge is not well understood. Here we show that PKMζ inhibition, via injection of zeta inhibitory peptide (ZIP) into the rat sensorimotor cortex, disrupts sensorimotor memories for a skilled reaching task even after several weeks of training. The rate of relearning the task after the memory disruption by ZIP was indistinguishable from the rate of initial learning, suggesting no significant savings after the memory loss. These results indicate a shared molecular mechanism of storage for declarative and procedural forms of memory

Modeling the Population-Level Effects of Male Circumcision as an HIV-Preventive Measure: A Gendered Perspective

BACKGROUND: Evidence from biological, epidemiological, and controlled intervention studies has demonstrated that male circumcision (MC) protects males from HIV infection, and MC is now advocated as a public-health intervention against HIV. MC provides direct protection only to men, but is expected to provide indirect protection to women at risk of acquiring HIV from heterosexual transmission. How such indirect protection interacts with the possibility that MC campaigns will lead to behavior changes, however, is not yet well understood. Our objective here is to investigate the link between individual-level effects of MC campaigns and long-term population-level outcomes resulting from disease dynamics, looking at both genders separately, over a broad range of parameters. METHODS AND FINDINGS: We use simple mathematical models of heterosexual transmission to investigate the potential effects of a circumcision scale-up, combined with possible associated behavioral disinhibition. We examine patterns in expected long-term prevalence using a simple equilibrium model based on transmission factors, and validate our results with ODE-based simulations, focusing on the link between effects on females and those on males.We find that the long-term population-level effects on females and males are not strongly linked: there are many possible ways in which an intervention which reduces prevalence in males might nonetheless increase prevalence in females. CONCLUSIONS: Since an intervention that reduces long-term male prevalence could nonetheless increase long-term female prevalence, MC campaigns should explicitly consider both the short-term and long-term effects of MC interventions on females. Our findings strongly underline the importance of pairing MC programs with education, support programs and HIV testing and counseling, together with other prevention measures