Scalable inference of heterogeneous reaction kinetics from pooled single-cell recordings.

Mathematical methods combined with measurements of single-cell dynamics provide a means to reconstruct intracellular processes that are only partly or indirectly accessible experimentally. To obtain reliable reconstructions, the pooling of measurements from several cells of a clonal population is mandatory. However, cell-to-cell variability originating from diverse sources poses computational challenges for such process reconstruction. We introduce a scalable Bayesian inference framework that properly accounts for population heterogeneity. The method allows inference of inaccessible molecular states and kinetic parameters; computation of Bayes factors for model selection; and dissection of intrinsic, extrinsic and technical noise. We show how additional single-cell readouts such as morphological features can be included in the analysis. We use the method to reconstruct the expression dynamics of a gene under an inducible promoter in yeast from time-lapse microscopy data.

Health care provider choice in the case of patient-initiated contacts. An extended version of discrete choice of model demand

This paper analyzes the nature of health care provider choice inthe case of patient-initiated contacts, with special reference toa National Health Service setting, where monetary prices are zeroand general practitioners act as gatekeepers to publicly financedspecialized care. We focus our attention on the factors that mayexplain the continuously increasing use of hospital emergencyvisits as opposed to other provider alternatives. An extendedversion of a discrete choice model of demand for patient-initiatedcontacts is presented, allowing for individual and town residencesize differences in perceived quality (preferences) betweenalternative providers and including travel and waiting time asnon-monetary costs. Results of a nested multinomial logit model ofprovider choice are presented. Individual choice betweenalternatives considers, in a repeated nested structure, self-care,primary care, hospital and clinic emergency services. Welfareimplications and income effects are analyzed by computingcompensating variations, and by simulating the effects of userfees by levels of income. Results indicate that compensatingvariation per visit is higher than the direct marginal cost ofemergency visits, and consequently, emergency visits do not appearas an inefficient alternative even for non-urgent conditions.

Revenue management under general discrete choice model of consumer behavior

Customer choice behavior, such as 'buy-up' and 'buy-down', is an importantphe-nomenon in a wide range of industries. Yet there are few models ormethodologies available to exploit this phenomenon within yield managementsystems. We make some progress on filling this void. Specifically, wedevelop a model of yield management in which the buyers' behavior ismodeled explicitly using a multi-nomial logit model of demand. Thecontrol problem is to decide which subset of fare classes to offer ateach point in time. The set of open fare classes then affects the purchaseprobabilities for each class. We formulate a dynamic program todetermine the optimal control policy and show that it reduces to a dynamicnested allocation policy. Thus, the optimal choice-based policy caneasily be implemented in reservation systems that use nested allocationcontrols. We also develop an estimation procedure for our model based onthe expectation-maximization (EM) method that jointly estimates arrivalrates and choice model parameters when no-purchase outcomes areunobservable. Numerical results show that this combined optimization-estimation approach may significantly improve revenue performancerelative to traditional leg-based models that do not account for choicebehavior.

An enhanced concave program relaxation for choice network revenue management

The network choice revenue management problem models customers as choosing from an offer-set, andthe firm decides the best subset to offer at any given moment to maximize expected revenue. The resultingdynamic program for the firm is intractable and approximated by a deterministic linear programcalled the CDLP which has an exponential number of columns. However, under the choice-set paradigmwhen the segment consideration sets overlap, the CDLP is difficult to solve. Column generation has beenproposed but finding an entering column has been shown to be NP-hard. In this paper, starting with aconcave program formulation based on segment-level consideration sets called SDCP, we add a class ofconstraints called product constraints, that project onto subsets of intersections. In addition we proposea natural direct tightening of the SDCP called ?SDCP, and compare the performance of both methodson the benchmark data sets in the literature. Both the product constraints and the ?SDCP method arevery simple and easy to implement and are applicable to the case of overlapping segment considerationsets. In our computational testing on the benchmark data sets in the literature, SDCP with productconstraints achieves the CDLP value at a fraction of the CPU time taken by column generation and webelieve is a very promising approach for quickly approximating CDLP when segment consideration setsoverlap and the consideration sets themselves are relatively small.

Application of real-time PCR for total airborne bacterial assessment : comparison with epifluorescence microscopy and culture-dependent methods

Traditional culture-dependent methods to quantify and identify airborne microorganisms are limited by factors such as short-duration sampling times and inability to count nonculturableor non-viable bacteria. Consequently, the quantitative assessment of bioaerosols is often underestimated. Use of the real-time quantitative polymerase chain reaction (Q-PCR) to quantify bacteria in environmental samples presents an alternative method, which should overcome this problem. The aim of this study was to evaluate the performance of a real-time Q-PCR assay as a simple and reliable way to quantify the airborne bacterial load within poultry houses and sewage treatment plants, in comparison with epifluorescencemicroscopy and culture-dependent methods. The estimates of bacterial load that we obtained from real-time PCR and epifluorescence methods, are comparable, however, our analysis of sewage treatment plants indicate these methods give values 270-290 fold greater than those obtained by the ''impaction on nutrient agar'' method. The culture-dependent method of air impaction on nutrient agar was also inadequate in poultry houses, as was the impinger-culture method, which gave a bacterial load estimate 32-fold lower than obtained by Q-PCR. Real-time quantitative PCR thus proves to be a reliable, discerning, and simple method that could be used to estimate airborne bacterial load in a broad variety of other environments expected to carry high numbers of airborne bacteria. [Authors]

Impact of Matrix-Assisted Laser Desorption Ionization Time-of-Flight Mass Spectrometry on the Clinical Management of Patients With Gram-negative Bacteremia: A Prospective Observational Study.

Background. Early identification of pathogens from blood cultures using matrix-assisted laser desorption ionization time-of-flight (MALDI-TOF) mass spectrometry may optimize the choice of empirical antibiotic therapy in the setting of bloodstream infections. We aimed to assess the impact of this new technology on the use of antibiotic treatment in patients with gram-negative bacteremia. Methods. We conducted a prospective observational study from January to December 2010 to evaluate the sequential and separate impacts of Gram stain reporting and MALDI-TOF bacterial identification performed on blood culture pellets in patients with gram-negative bacteremia. The primary outcome was the impact of MALDI-TOF on empirical antibiotic choice. Results. Among 202 episodes of gram-negative bacteremia, Gram stain reporting had an impact in 42 cases (20.8%). MALDI-TOF identification led to a modification of empirical therapy in 71 of all 202 cases (35.1%), and in 16 of 27 cases (59.3%) of monomicrobial bacteremia caused by AmpC-producing Enterobacteriaceae. The most frequently observed impact was an early appropriate broadening of the antibiotic spectrum in 31 of 71 cases (43.7%). In total, 143 of 165 episodes (86.7%) of monomicrobial bacteremia were correctly identified at genus level by MALDI-TOF. Conclusions. In a low prevalence area for extended spectrum betalactamases (ESBL) and multiresistant gram-negative bacteria, MALDI-TOF performed on blood culture pellets had an impact on the clinical management of 35.1% of all gram-negative bacteremia cases, demonstrating a greater impact than Gram stain reporting. Thus, MALDI-TOF could become a vital second step beside Gram stain in guiding the empirical treatment of patients with bloodstream infection.

Issue of time in generally covariant theories and the Komar-Bergmann approach to observables in general relativity

Diffeomorphism-induced symmetry transformations and time evolution are distinct operations in generally covariant theories formulated in phase space. Time is not frozen. Diffeomorphism invariants are consequently not necessarily constants of the motion. Time-dependent invariants arise through the choice of an intrinsic time, or equivalently through the imposition of time-dependent gauge fixation conditions. One example of such a time-dependent gauge fixing is the Komar-Bergmann use of Weyl curvature scalars in general relativity. An analogous gauge fixing is also imposed for the relativistic free particle and the resulting complete set time-dependent invariants for this exactly solvable model are displayed. In contrast with the free particle case, we show that gauge invariants that are simultaneously constants of motion cannot exist in general relativity. They vary with intrinsic time.

Homeodomain binding motifs modulate the probability of odorant receptor gene choice in transgenic mice.

Odorant receptor (OR) genes constitute with 1200 members the largest gene family in the mouse genome. A mature olfactory sensory neuron (OSN) is thought to express just one OR gene, and from one allele. The cell bodies of OSNs that express a given OR gene display a mosaic pattern within a particular region of the main olfactory epithelium. The mechanisms and cis-acting DNA elements that regulate the expression of one OR gene per OSN - OR gene choice - remain poorly understood. Here, we describe a reporter assay to identify minimal promoters for OR genes in transgenic mice, which are produced by the conventional method of pronuclear injection of DNA. The promoter transgenes are devoid of an OR coding sequence, and instead drive expression of the axonal marker tau-β-galactosidase. For four mouse OR genes (M71, M72, MOR23, and P3) and one human OR gene (hM72), a mosaic, OSN-specific pattern of reporter expression can be obtained in transgenic mice with contiguous DNA segments of only ~300 bp that are centered around the transcription start site (TSS). The ~150bp region upstream of the TSS contains three conserved sequence motifs, including homeodomain (HD) binding sites. Such HD binding sites are also present in the H and P elements, DNA sequences that are known to strongly influence OR gene expression. When a 19mer encompassing a HD binding site from the P element is multimerized nine times and added upstream of a MOR23 minigene that contains the MOR23 coding region, we observe a dramatic increase in the number of transgene-expressing founders and lines and in the number of labeled OSNs. By contrast, a nine times multimerized 19mer with a mutant HD binding site does not have these effects. We hypothesize that HD binding sites in the H and P elements and in OR promoters modulate the probability of OR gene choice.

Detection of live and antibiotic-killed bacteria by quantitative real-time PCR of specific fragments of rRNA.

Assessing bacterial viability by molecular markers might help accelerate the measurement of antibiotic-induced killing. This study investigated whether rRNA could be suitable for this purpose. Cultures of penicillin-susceptible and penicillin-tolerant (Tol1 mutant) Streptococcus gordonii were exposed to mechanistically different penicillin and levofloxacin. Bacterial survival was assessed by viable counts and compared to quantitative real-time PCR amplification of either the 16S rRNA genes or the 16S rRNA, following reverse transcription. Penicillin-susceptible S. gordonii lost > or =4 log(10) CFU/ml of viability over 48 h of penicillin treatment. In comparison, the Tol1 mutant lost < or =1 log(10) CFU/ml. Amplification of a 427-bp fragment of 16S rRNA genes yielded amplicons that increased proportionally to viable counts during bacterial growth but did not decrease during drug-induced killing. In contrast, the same 427-bp fragment amplified from 16S rRNA paralleled both bacterial growth and drug-induced killing. It also differentiated between penicillin-induced killing of the parent and the Tol1 mutant (> or =4 log(10) CFU/ml and < or =1 log(10) CFU/ml, respectively) and detected killing by mechanistically unrelated levofloxacin. Since large fragments of polynucleotides might be degraded faster than smaller fragments, the experiments were repeated by amplifying a 119-bp region internal to the original 427-bp fragment. The amount of 119-bp amplicons increased proportionally to viability during growth but remained stable during drug treatment. Thus, 16S rRNA was a marker of antibiotic-induced killing, but the size of the amplified fragment was critical for differentiation between live and dead bacteria.

Front propagation in the A+B->2A reaction under subdiffusion

We consider an irreversible autocatalytic conversion reaction A+B->2A under subdiffusion described by continuous-time random walks. The reactants transformations take place independently of their motion and are described by constant rates. The analog of this reaction in the case of normal diffusion is described by the Fisher-Kolmogorov-Petrovskii-Piskunov equation leading to the existence of a nonzero minimal front propagation velocity, which is really attained by the front in its stable motion. We show that for subdiffusion, this minimal propagation velocity is zero, which suggests propagation failure.

State-to-state reaction probabilities within the quantum transition state framework

Rigorous quantum dynamics calculations of reaction rates and initial state-selected reaction probabilities of polyatomic reactions can be efficiently performed within the quantum transition state concept employing flux correlation functions and wave packet propagation utilizing the multi-configurational time-dependent Hartree approach. Here, analytical formulas and a numerical scheme extending this approach to the calculation of state-to-state reaction probabilities are presented. The formulas derived facilitate the use of three different dividing surfaces: two dividing surfaces located in the product and reactant asymptotic region facilitate full state resolution while a third dividing surface placed in the transition state region can be used to define an additional flux operator. The eigenstates of the corresponding thermal flux operator then correspond to vibrational states of the activated complex. Transforming these states to reactant and product coordinates and propagating them into the respective asymptotic region, the full scattering matrix can be obtained. To illustrate the new approach, test calculations study the D + H2(ν, j) → HD(ν′, j′) + H reaction for J = 0.

Development of a real-time PCR for the specific detection of Waddlia chondrophila in clinical samples.

Waddlia chondrophila is considered as an emerging human pathogen likely involved in miscarriage and lower respiratory tract infections. Given the low sensitivity of cell culture to recover such an obligate intracellular bacteria, molecular-based diagnostic approaches are warranted. We thus developed a real-time PCR that amplifies Waddlia chondrophila DNA. Specific primers and probe were selected to target the 16S rRNA gene. The PCR specifically amplified W. chondrophila but did not amplify other related-bacteria such as Parachlamydia acanthamoebae, Simkania negevensis and Chlamydia pneumoniae. The PCR exhibited a good intra-run and inter-run reproducibility and a sensitivity of less than ten copies of the positive control. This real-time PCR was then applied to 32 nasopharyngeal aspirates taken from children with bronchiolitis not due to respiratory syncytial virus (RSV). Three samples revealed to be Waddlia positive, suggesting a possible role of this Chlamydia-related bacteria in this setting.

Mating order-dependent female mate choice in the polygynandrous common lizard Lacerta vivipara.

Recent studies indicate that directional female mate choice and order-dependent female mate choice importantly contribute to non-random mating patterns. In species where females prefer larger sized males, disentangling different hypotheses leading to non-random mating patterns is especially difficult, given that male size usually correlates with behaviours that may lead to non-random mating (e.g. size-dependent emergence from hibernation, male fighting ability). Here we investigate female mate choice and order-dependent female mate choice in the polygynandrous common lizard (Lacerta vivipara). By sequentially presenting males in random order to females, we exclude non-random mating patterns potentially arising due to intra-sexual selection (e.g. male-male competition), trait-dependent encounter probabilities, trait-dependent conspicuousness, or trait-dependent emergence from hibernation. To test for order-dependent female mate choice we investigate whether the previous mating history affects female choice. We show that body size and body condition of the male with which a female mated for the first time were bigger and better, respectively, than the average body size and body condition of the rejected males. There was a negative correlation between body sizes of first and second copulating males. This indicates that female mate choice is dependent on the previous mating history and it shows that the female's choice criteria are non-static, i.e. non-directional. Our study therefore suggests that context-dependent female mate choice may not only arise due to genotype-environment interactions, but also due to other female mating strategies, i.e. order-dependent mate choice. Thus context-dependent female mate choice might be more frequent than previously thought.

Exposure to bioaerosols in poultry houses at different stages of fattening; use of real-time PCR for airborne bacterial quantification

Previous studies have demonstrated that poultry house workers are exposed to very high levels of organic dust and consequently have an increased prevalence of adverse respiratory symptoms. However, the influence of the age of broilers on bioaerosol concentrations has not been investigated. To evaluate the evolution of bioaerosol concentration during the fattening period, bioaerosol parameters (inhalable dust, endotoxin and bacteria) were measured in 12 poultry confinement buildings in Switzerland, at three different stages of the birds' growth; samples of air taken from within the breathing zones of individual poultry house employees as they caught the chickens ready to be transported for slaughter were also analysed. Quantitative polymerase chain reaction (Q-PCR) was used to assess the quantity of total airborne bacteria and total airborne Staphylococcus species. Bioaerosol levels increased significantly during the fattening period of the chickens. During the task of catching mature birds, the mean inhalable dust concentration for a worker was 26 +/- 1.9 mg m(-3) and endotoxin concentration was 6198 +/- 2.3 EU m(-3) air, >6-fold higher than the Swiss occupational recommended value (1000 EU m(-3)). The mean exposure level of bird catchers to total bacteria and Staphylococcus species measured by Q-PCR is also very high, respectively, reaching values of 53 (+/-2.6) x 10(7) cells m(-3) air and 62 (+/-1.9) x 10(6) m(-3) air. It was concluded that in the absence of wearing protective breathing apparatus, chicken catchers in Switzerland risk exposure beyond recommended limits for all measured bioaerosol parameters. Moreover, the use of Q-PCR to estimate total and specific numbers of airborne bacteria is a promising tool for evaluating any modifications intended to improve the safety of current working practices

Validation of real-time methylation-specific PCR to determine O6-methylguanine-DNA methyltransferase gene promoter methylation in glioma.

Epigenetic silencing of the DNA repair protein O(6)-methylguanine-DNA methyltransferase (MGMT) by promoter methylation predicts successful alkylating agent therapy, such as with temozolomide, in glioblastoma patients. Stratified therapy assignment of patients in prospective clinical trials according to tumor MGMT status requires a standardized diagnostic test, suitable for high-throughput analysis of small amounts of formalin-fixed, paraffin-embedded tumor tissue. A direct, real-time methylation-specific PCR (MSP) assay was developed to determine methylation status of the MGMT gene promoter. Assay specificity was obtained by selective amplification of methylated DNA sequences of sodium bisulfite-modified DNA. The copy number of the methylated MGMT promoter, normalized to the beta-actin gene, provides a quantitative test result. We analyzed 134 clinical glioma samples, comparing the new test with the previously validated nested gel-based MSP assay, which yields a binary readout. A cut-off value for the MGMT methylation status was suggested by fitting a bimodal normal mixture model to the real-time results, supporting the hypothesis that there are two distinct populations within the test samples. Comparison of the tests showed high concordance of the results (82/91 [90%]; Cohen's kappa = 0.80; 95% confidence interval, 0.82-0.95). The direct, real-time MSP assay was highly reproducible (Pearson correlation 0.996) and showed valid test results for 93% (125/134) of samples compared with 75% (94/125) for the nested, gel-based MSP assay. This high-throughput test provides an important pharmacogenomic tool for individualized management of alkylating agent chemotherapy.