The feasibility of forecasting influenza epidemics in Cuba

A large influenza epidemic took place in Havana during the winter of 1988. The epidemiologic surveillance unit of the Pedro Kouri Institute of Tropical Medicine detected the begining of the epidemic wave. The Rvachev-Baroyan mathematical model of the geographic spread of an epidemic was used to forecast this epidemic under routine conditions of the public health system. The expected number of individuals who would attend outpatient services, because of influenza-like illness, was calculated and communicated to the health authorities within enough time to permit the introduction of available control measures. The approximate date of the epidemic peak, the daily expected number of individuals attending medical services, and the approximate time of the end of the epidemic wave were estimated. The prediction error was 12%. The model was sufficienty accurate to warrant its use as a pratical forecasting tool in the Cuban public health system.

As the world turns: short-term human spatial memory in egocentric and allocentric coordinates

We aimed to determine whether human subjects' reliance on different sources of spatial information encoded in different frames of reference (i.e., egocentric versus allocentric) affects their performance, decision time and memory capacity in a short-term spatial memory task performed in the real world. Subjects were asked to play the Memory game (a.k.a. the Concentration game) without an opponent, in four different conditions that controlled for the subjects' reliance on egocentric and/or allocentric frames of reference for the elaboration of a spatial representation of the image locations enabling maximal efficiency. We report experimental data from young adult men and women, and describe a mathematical model to estimate human short-term spatial memory capacity. We found that short-term spatial memory capacity was greatest when an egocentric spatial frame of reference enabled subjects to encode and remember the image locations. However, when egocentric information was not reliable, short-term spatial memory capacity was greater and decision time shorter when an allocentric representation of the image locations with respect to distant objects in the surrounding environment was available, as compared to when only a spatial representation encoding the relationships between the individual images, independent of the surrounding environment, was available. Our findings thus further demonstrate that changes in viewpoint produced by the movement of images placed in front of a stationary subject is not equivalent to the movement of the subject around stationary images. We discuss possible limitations of classical neuropsychological and virtual reality experiments of spatial memory, which typically restrict the sensory information normally available to human subjects in the real world.

Alteration of brain glycogen turnover in the conscious rat after 5h of prolonged wakefulness.

Although glycogen (Glyc) is the main carbohydrate storage component, the role of Glyc in the brain during prolonged wakefulness is not clear. The aim of this study was to determine brain Glyc concentration ([]) and turnover time (tau) in euglycemic conscious and undisturbed rats, compared to rats maintained awake for 5h. To measure the metabolism of [1-(13)C]-labeled Glc into Glyc, 23 rats received a [1-(13)C]-labeled Glc solution as drink (10% weight per volume in tap water) ad libitum as their sole source of exogenous carbon for a "labeling period" of either 5h (n=13), 24h (n=5) or 48 h (n=5). Six of the rats labeled for 5h were continuously maintained awake by acoustic, tactile and olfactory stimuli during the labeling period, which resulted in slightly elevated corticosterone levels. Brain [Glyc] measured biochemically after focused microwave fixation in the rats maintained awake (3.9+/-0.2 micromol/g, n=6) was not significantly different from that of the control group (4.0+/-0.1 micromol/g, n=7; t-test, P>0.5). To account for potential variations in plasma Glc isotopic enrichment (IE), Glyc IE was normalized by N-acetyl-aspartate (NAA) IE. A simple mathematical model was developed to derive brain Glyc turnover time as 5.3h with a fit error of 3.2h and NAA turnover time as 15.6h with a fit error of 6.5h, in the control rats. A faster tau(Glyc) (2.9h with a fit error of 1.2h) was estimated in the rats maintained awake for 5h. In conclusion, 5h of prolonged wakefulness mainly activates glycogen metabolism, but has minimal effect on brain [Glyc].

Why are slip lengths so large in carbon nanotubes?

The enhanced flow in carbon nanotubes is explained using a mathematical model that includes a depletion layer with reduced viscosity near the wall. In the limit of large tubes the model predicts no noticeable enhancement. For smaller tubes the model predicts enhancement that increases as the radius decreases. An analogy between the reduced viscosity and slip-length models shows that the term slip-length is misleading and that on surfaces which are smooth at the nanoscale it may be thought of as a length-scale associated with the size of the depletion region and viscosity ratio. The model therefore provides a physical interpretation of the classical Navier slip condition and explains why `slip-lengths' may be greater than the tube radius.

Evolution of decision-making and learning under fluctuating social environments

The capacity to learn to associate sensory perceptions with appropriate motor actions underlies the success of many animal species, from insects to humans. The evolutionary significance of learning has long been a subject of interest for evolutionary biologists who emphasize the bene¬fit yielded by learning under changing environmental conditions, where it is required to flexibly switch from one behavior to another. However, two unsolved questions are particularly impor¬tant for improving our knowledge of the evolutionary advantages provided by learning, and are addressed in the present work. First, because it is possible to learn the wrong behavior when a task is too complex, the learning rules and their underlying psychological characteristics that generate truly adaptive behavior must be identified with greater precision, and must be linked to the specific ecological problems faced by each species. A framework for predicting behavior from the definition of a learning rule is developed here. Learning rules capture cognitive features such as the tendency to explore, or the ability to infer rewards associated to unchosen actions. It is shown that these features interact in a non-intuitive way to generate adaptive behavior in social interactions where individuals affect each other's fitness. Such behavioral predictions are used in an evolutionary model to demonstrate that, surprisingly, simple trial-and-error learn¬ing is not always outcompeted by more computationally demanding inference-based learning, when population members interact in pairwise social interactions. A second question in the evolution of learning is its link with and relative advantage compared to other simpler forms of phenotypic plasticity. After providing a conceptual clarification on the distinction between genetically determined vs. learned responses to environmental stimuli, a new factor in the evo¬lution of learning is proposed: environmental complexity. A simple mathematical model shows that a measure of environmental complexity, the number of possible stimuli in one's environ¬ment, is critical for the evolution of learning. In conclusion, this work opens roads for modeling interactions between evolving species and their environment in order to predict how natural se¬lection shapes animals' cognitive abilities. - La capacité d'apprendre à associer des sensations perceptives à des actions motrices appropriées est sous-jacente au succès évolutif de nombreuses espèces, depuis les insectes jusqu'aux êtres hu¬mains. L'importance évolutive de l'apprentissage est depuis longtemps un sujet d'intérêt pour les biologistes de l'évolution, et ces derniers mettent l'accent sur le bénéfice de l'apprentissage lorsque les conditions environnementales sont changeantes, car dans ce cas il est nécessaire de passer de manière flexible d'un comportement à l'autre. Cependant, deux questions non résolues sont importantes afin d'améliorer notre savoir quant aux avantages évolutifs procurés par l'apprentissage. Premièrement, puisqu'il est possible d'apprendre un comportement incorrect quand une tâche est trop complexe, les règles d'apprentissage qui permettent d'atteindre un com¬portement réellement adaptatif doivent être identifiées avec une plus grande précision, et doivent être mises en relation avec les problèmes écologiques spécifiques rencontrés par chaque espèce. Un cadre théorique ayant pour but de prédire le comportement à partir de la définition d'une règle d'apprentissage est développé ici. Il est démontré que les caractéristiques cognitives, telles que la tendance à explorer ou la capacité d'inférer les récompenses liées à des actions non ex¬périmentées, interagissent de manière non-intuitive dans les interactions sociales pour produire des comportements adaptatifs. Ces prédictions comportementales sont utilisées dans un modèle évolutif afin de démontrer que, de manière surprenante, l'apprentissage simple par essai-et-erreur n'est pas toujours battu par l'apprentissage basé sur l'inférence qui est pourtant plus exigeant en puissance de calcul, lorsque les membres d'une population interagissent socialement par pair. Une deuxième question quant à l'évolution de l'apprentissage concerne son lien et son avantage relatif vis-à-vis d'autres formes plus simples de plasticité phénotypique. Après avoir clarifié la distinction entre réponses aux stimuli génétiquement déterminées ou apprises, un nouveau fac¬teur favorisant l'évolution de l'apprentissage est proposé : la complexité environnementale. Un modèle mathématique permet de montrer qu'une mesure de la complexité environnementale - le nombre de stimuli rencontrés dans l'environnement - a un rôle fondamental pour l'évolution de l'apprentissage. En conclusion, ce travail ouvre de nombreuses perspectives quant à la mo¬délisation des interactions entre les espèces en évolution et leur environnement, dans le but de comprendre comment la sélection naturelle façonne les capacités cognitives des animaux.

The stabilizing effects of the acquired immunity on the schistosomiasis transmission modeling - the sensitivity Analysis

A mathematical model is proposed to analyze the effects of acquired immunity on the transmission of schistosomiasis in the human host. From this model the prevalence curve dependent on four parameters can be obtained. These parameters were estimated fitting the data by the maximum likelihood method. The model showed a good retrieving capacity of real data from two endemic areas of schistosomiasis: Touros, Brazil (Schistosoma mansoni) and Misungwi, Tanzania (S. haematobium). Also, the average worm burden per person and the dispersion of parasite per person in the community can be obtained from the model. In this paper, the stabilizing effects of the acquired immunity assumption in the model are assessed in terms of the epidemiological variables as follows. Regarded to the prevalence curve, we calculate the confidence interval, and related to the average worm burden and the worm dispersion in the community, the sensitivity analysis (the range of the variation) of both variables with respect to their parameters is performed.

Differences in whole-body fat oxidation kinetics between cycling and running.

This study aimed to quantitatively describe and compare whole-body fat oxidation kinetics in cycling and running using a sinusoidal mathematical model (SIN). Thirteen moderately trained individuals (7 men and 6 women) performed two graded exercise tests, with 3-min stages and 1 km h(-1) (or 20 W) increment, on a treadmill and on a cycle ergometer. Fat oxidation rates were determined using indirect calorimetry and plotted as a function of exercise intensity. The SIN model, which includes three independent variables (dilatation, symmetry and translation) that account for main quantitative characteristics of kinetics, provided a mathematical description of fat oxidation kinetics and allowed for determination of the intensity (Fat(max)) that elicits maximal fat oxidation (MFO). While the mean fat oxidation kinetics in cycling formed a symmetric parabolic curve, the mean kinetics during running was characterized by a greater dilatation (i.e., widening of the curve, P < 0.001) and a rightward asymmetry (i.e., shift of the peak of the curve to higher intensities, P = 0.01). Fat(max) was significantly higher in running compared with cycling (P < 0.001), whereas MFO was not significantly different between modes of exercise (P = 0.36). This study showed that the whole-body fat oxidation kinetics during running was characterized by a greater dilatation and a rightward asymmetry compared with cycling. The greater dilatation may be mainly related to the larger muscle mass involved in running while the rightward asymmetry may be induced by the specific type of muscle contraction.

Uncertainties regarding dengue modeling in Rio de Janeiro, Brazil

Dengue fever is currently the most important arthropod-borne viral disease in Brazil. Mathematical modeling of disease dynamics is a very useful tool for the evaluation of control measures. To be used in decision-making, however, a mathematical model must be carefully parameterized and validated with epidemiological and entomological data. In this work, we developed a simple dengue model to answer three questions: (i) which parameters are worth pursuing in the field in order to develop a dengue transmission model for Brazilian cities; (ii) how vector density spatial heterogeneity influences control efforts; (iii) with a degree of uncertainty, what is the invasion potential of dengue virus type 4 (DEN-4) in Rio de Janeiro city. Our model consists of an expression for the basic reproductive number (R0) that incorporates vector density spatial heterogeneity. To deal with the uncertainty regarding parameter values, we parameterized the model using a priori probability density functions covering a range of plausible values for each parameter. Using the Latin Hypercube Sampling procedure, values for the parameters were generated. We conclude that, even in the presence of vector spatial heterogeneity, the two most important entomological parameters to be estimated in the field are the mortality rate and the extrinsic incubation period. The spatial heterogeneity of the vector population increases the risk of epidemics and makes the control strategies more complex. At last, we conclude that Rio de Janeiro is at risk of a DEN-4 invasion. Finally, we stress the point that epidemiologists, mathematicians, and entomologists need to interact more to find better approaches to the measuring and interpretation of the transmission dynamics of arthropod-borne diseases.

A Novel Pulse-Chase SILAC Strategy Measures Changes in Protein Decay and Synthesis Rates Induced by Perturbation of Proteostasis with an Hsp90 Inhibitor.

Standard proteomics methods allow the relative quantitation of levels of thousands of proteins in two or more samples. While such methods are invaluable for defining the variations in protein concentrations which follow the perturbation of a biological system, they do not offer information on the mechanisms underlying such changes. Expanding on previous work [1], we developed a pulse-chase (pc) variant of SILAC (stable isotope labeling by amino acids in cell culture). pcSILAC can quantitate in one experiment and for two conditions the relative levels of proteins newly synthesized in a given time as well as the relative levels of remaining preexisting proteins. We validated the method studying the drug-mediated inhibition of the Hsp90 molecular chaperone, which is known to lead to increased synthesis of stress response proteins as well as the increased decay of Hsp90 "clients". We showed that pcSILAC can give information on changes in global cellular proteostasis induced by treatment with the inhibitor, which are normally not captured by standard relative quantitation techniques. Furthermore, we have developed a mathematical model and computational framework that uses pcSILAC data to determine degradation constants kd and synthesis rates Vs for proteins in both control and drug-treated cells. The results show that Hsp90 inhibition induced a generalized slowdown of protein synthesis and an increase in protein decay. Treatment with the inhibitor also resulted in widespread protein-specific changes in relative synthesis rates, together with variations in protein decay rates. The latter were more restricted to individual proteins or protein families than the variations in synthesis. Our results establish pcSILAC as a viable workflow for the mechanistic dissection of changes in the proteome which follow perturbations. Data are available via ProteomeXchange with identifier PXD000538.

Radiation induced apoptosis and initial DNA damage are inversely related in locally advanced breast cancer patients.

Background. DNA-damage assays, quantifying the initial number of DNA double-strand breaks induced by radiation, have been proposed as a predictive test for radiation-induced toxicity. Determination of radiation-induced apoptosis in peripheral blood lymphocytes by flow cytometry analysis has also been proposed as an approach for predicting normal tissue responses following radiotherapy. The aim of the present study was to explore the association between initial DNA damage, estimated by the number of double-strand breaks induced by a given radiation dose, and the radio-induced apoptosis rates observed. Methods. Peripheral blood lymphocytes were taken from 26 consecutive patients with locally advanced breast carcinoma. Radiosensitivity of lymphocytes was quantified as the initial number of DNA double-strand breaks induced per Gy and per DNA unit (200 Mbp). Radio-induced apoptosis at 1, 2 and 8 Gy was measured by flow cytometry using annexin V/propidium iodide. Results. Radiation-induced apoptosis increased in order to radiation dose and data fitted to a semi logarithmic mathematical model. A positive correlation was found among radio-induced apoptosis values at different radiation doses: 1, 2 and 8 Gy (p < 0.0001 in all cases). Mean DSB/Gy/DNA unit obtained was 1.70 ± 0.83 (range 0.63-4.08; median, 1.46). A statistically significant inverse correlation was found between initial damage to DNA and radio-induced apoptosis at 1 Gy (p = 0.034). A trend toward 2 Gy (p = 0.057) and 8 Gy (p = 0.067) was observed after 24 hours of incubation. Conclusions. An inverse association was observed for the first time between these variables, both considered as predictive factors to radiation toxicity.

The demographic benefits of belligerence and bravery: defeated group repopulation or victorious group size expansion?

Intraspecific coalitional aggression between groups of individuals is a widespread trait in the animal world. It occurs in invertebrates and vertebrates, and is prevalent in humans. What are the conditions under which coalitional aggression evolves in natural populations? In this article, I develop a mathematical model delineating conditions where natural selection can favor the coevolution of belligerence and bravery between small-scale societies. Belligerence increases an actor's group probability of trying to conquer another group and bravery increase the actors's group probability of defeating an attacked group. The model takes into account two different types of demographic scenarios that may lead to the coevolution of belligerence and bravery. Under the first, the fitness benefits driving the coevolution of belligerence and bravery come through the repopulation of defeated groups by fission of victorious ones. Under the second demographic scenario, the fitness benefits come through a temporary increase in the local carrying capacity of victorious groups, after transfer of resources from defeated groups to victorious ones. The analysis of the model suggests that the selective pressures on belligerence and bravery are stronger when defeated groups can be repopulated by victorious ones. The analysis also suggests that, depending on the shape of the contest success function, costly bravery can evolve in groups of any size.

In vivo quantification of neuro-glial metabolism and glial glutamate concentration using 1H-[13C] MRS at 14.1T.

Astrocytes have recently become a major center of interest in neurochemistry with the discoveries on their major role in brain energy metabolism. An interesting way to probe this glial contribution is given by in vivo (13) C NMR spectroscopy coupled with the infusion labeled glial-specific substrate, such as acetate. In this study, we infused alpha-chloralose anesthetized rats with [2-(13) C]acetate and followed the dynamics of the fractional enrichment (FE) in the positions C4 and C3 of glutamate and glutamine with high sensitivity, using (1) H-[(13) C] magnetic resonance spectroscopy (MRS) at 14.1T. Applying a two-compartment mathematical model to the measured time courses yielded a glial tricarboxylic acid (TCA) cycle rate (Vg ) of 0.27 ± 0.02 μmol/g/min and a glutamatergic neurotransmission rate (VNT ) of 0.15 ± 0.01 μmol/g/min. Glial oxidative ATP metabolism thus accounts for 38% of total oxidative metabolism measured by NMR. Pyruvate carboxylase (VPC ) was 0.09 ± 0.01 μmol/g/min, corresponding to 37% of the glial glutamine synthesis rate. The glial and neuronal transmitochondrial fluxes (Vx (g) and Vx (n) ) were of the same order of magnitude as the respective TCA cycle fluxes. In addition, we estimated a glial glutamate pool size of 0.6 ± 0.1 μmol/g. The effect of spectral data quality on the fluxes estimates was analyzed by Monte Carlo simulations. In this (13) C-acetate labeling study, we propose a refined two-compartment analysis of brain energy metabolism based on (13) C turnover curves of acetate, glutamate and glutamine measured with state of the art in vivo dynamic MRS at high magnetic field in rats, enabling a deeper understanding of the specific role of glial cells in brain oxidative metabolism. In addition, the robustness of the metabolic fluxes determination relative to MRS data quality was carefully studied.

The risk of acquiring the new influenza A(H1N1) for Brazilian travelers to Chile, Argentina and the USA

We estimate the risk of acquiring the new influenza A(H1N1) for Brazilian travelers to Chile, Argentina and the USA. This is done by a mathematical model that quantifies the intensity of transmission of the new virus in those countries and the probability that one individual has of acquiring the influenza depending on the date of arrival and time spent in the area. The maximum estimated risk reached 7.5 cases per 10,000 visitors to Chile, 17 cases per 10,000 travelers to Argentina and 23 cases per 10,000 travelers to the USA. The estimated number of imported cases until 27 July is 57 ± 9 from Chile, 136 ± 27 from the USA and 301 ± 21 from Argentina, which are in accord with the official figures. Estimating the number of imported cases was particularly important for the moment of the disease introduction into this country, but it will certainly be important again as a tool to calculate the number of future imported cases from northern countries in our next inter-epidemic season, were imported cases can constitute again the majority of the new influenza burden to the Brazilian health services.

Projecting the effectiveness of RotaTeq® against rotavirus-related hospitalisations in Brazil

RotaTeq® (Merck & Company, Inc, Whitehouse Station, NJ, USA) is an oral pentavalent rotavirus vaccine (RV5) that has shown high and consistent efficacy in preventing rotavirus gastroenteritis (RGE) in randomised clinical trials previously conducted in industrialised countries with high medical care resources. To date, the efficacy and effectiveness data for RV5 are available in some Latin American countries, but not Brazil. In this analysis, we projected the effectiveness of RV5 in terms of the percentage reduction in RGE-related hospitalisations among children less than five years of age in four regions of Brazil, using a previously validated mathematical model. The model inputs included hospital-based rotavirus surveillance data from Goiânia, Porto Alegre, Salvador and São Paulo from 2005-2006, which provided the proportions of rotavirus attributable to serotypes G1, G2, G3, G4 and G9, and published rotavirus serotype-specific efficacy from the Rotavirus Efficacy and Safety Trial. The model projected an overall percentage reduction of 93% in RGE-related hospitalisations, with an estimated annual reduction in RGE-related hospitalisations between 42,991-77,383 in the four combined regions of Brazil. These results suggest that RV5 could substantially prevent RGE-related hospitalisations in Brazil.

Potential impact of a presumed increase in the biting activity of dengue-virus-infected Aedes aegypti (Diptera: Culicidae) females on virus transmission dynamics

Recently, we showed that infection with dengue virus increases the locomotor activity of Aedes aegypti females. We speculate that the observed increased locomotor activity could potentially increase the chances of finding a suitable host and, as a consequence, the relative biting rate of infected mosquitoes. We used a mathematical model to investigate the impact of the increased locomotor activity by assuming that this activity translated into an increased biting rate for infected mosquitoes. The results show that the increased biting rate resulted in dengue outbreaks with greater numbers of primary and secondary infections and with more severe biennial epidemics.