The non-equilibrium nature of culinary evolution

Food is an essential part of civilization, with a scope that ranges from the biological to the economic and cultural levels. Here, we study the statistics of ingredients and recipes taken from Brazilian, British, French and Medieval cookery books. We find universal distributions with scale invariant behaviour. We propose a copy-mutate process to model culinary evolution that fits our empirical data very well. We find a cultural 'founder effect' produced by the non-equilibrium dynamics of the model. Both the invariant and idiosyncratic aspects of culture are accounted for by our model, which may have applications in other kinds of evolutionary processes.

Active Dendrites Enhance Neuronal Dynamic Range

Since the first experimental evidences of active conductances in dendrites, most neurons have been shown to exhibit dendritic excitability through the expression of a variety of voltage-gated ion channels. However, despite experimental and theoretical efforts undertaken in the past decades, the role of this excitability for some kind of dendritic computation has remained elusive. Here we show that, owing to very general properties of excitable media, the average output of a model of an active dendritic tree is a highly non-linear function of its afferent rate, attaining extremely large dynamic ranges (above 50 dB). Moreover, the model yields double-sigmoid response functions as experimentally observed in retinal ganglion cells. We claim that enhancement of dynamic range is the primary functional role of active dendritic conductances. We predict that neurons with larger dendritic trees should have larger dynamic range and that blocking of active conductances should lead to a decrease in dynamic range.

Loxosceles gaucho Venom-Induced Acute Kidney Injury - In Vivo and In Vitro Studies

Background: Accidents caused by Loxosceles spider may cause severe systemic reactions, including acute kidney injury (AKI). There are few experimental studies assessing Loxosceles venom effects on kidney function in vivo. Methodology/Principal Findings: In order to test Loxosceles gaucho venom (LV) nephrotoxicity and to assess some of the possible mechanisms of renal injury, rats were studied up to 60 minutes after LV 0.24 mg/kg or saline IV injection (control). LV caused a sharp and significant drop in glomerular filtration rate, renal blood flow and urinary output and increased renal vascular resistance, without changing blood pressure. Venom infusion increased significantly serum creatine kinase and aspartate aminotransferase. In the LV group renal histology analysis found acute epithelial tubular cells degenerative changes, presence of cell debris and detached epithelial cells in tubular lumen without glomerular or vascular changes. Immunohistochemistry disclosed renal deposition of myoglobin and hemoglobin. LV did not cause injury to a suspension of fresh proximal tubules isolated from rats. Conclusions/Significance: Loxosceles gaucho venom injection caused early AKI, which occurred without blood pressure variation. Changes in glomerular function occurred likely due to renal vasoconstriction and rhabdomyolysis. Direct nephrotoxicity could not be demonstrated in vitro. The development of a consistent model of Loxosceles venom-induced AKI and a better understanding of the mechanisms involved in the renal injury may allow more efficient ways to prevent or attenuate the systemic injury after Loxosceles bite.

Immunomodulatory effects of recombinant BCG expressing pertussis toxin on TNF-alpha and IL-10 in a bladder cancer model

Background: Since successful treatment of superficial bladder cancer with BCG requires proper induction of Th1 immunity, we have developed a rBCG-S1PT strain that induced a stronger cellular immune response than BCG. This preclinical study was designed to compare the modulatory effects of BCG and rBCG-S1PT on bladder TNF-alpha and IL-10 expression and to evaluate antitumour activity. Methods: For Experiment I, the MB49 bladder cancer cell line was used in C57BL/6 mice. Chemical cauterization of the bladder was performed to promote intravesical tumor implantation. Mice were treated by intravesical instillation with BCG, rBCG-S1PT or PBS once a week for four weeks. After 35 days the bladders were removed and weighed. TNF-

Comparison of the effects of carbon dioxide and helium pneumoperitoneum on renal function

Background and Purpose: Carbon dioxide pneumoperitoneum is associated with significant hypercarbia and acidosis. The aim of this study is to evaluate the effects of carbon dioxide and helium pneumoperitoneum on renal function. Materials and Methods: Thirty adult dogs were put randomly into one of three groups ( n = 10 animals each): group A - pneumoperitoneum not performed; group B - CO2 pneumoperitoneum; and group C - helium pneumoperitoneum. The groups were analyzed with consideration given to body weight, hematologic values, hemodynamic parameters ( heart rate, mean arterial pressure, central venous pressure, cardiac output, stroke volume, systemic vascular resistance, pulmonary vascular resistance, left cardiac work index, cardiac index, mean pulmonary artery pressure, and pulmonary capillary wedge pressure), and renal function ( plasma renin activity, urinary output, creatinine clearance, and sodium excretory fraction). Results: An accentuated decrease in urinary output was observed during pneumoperitoneum in groups B and C compared to the control group. In groups B and C, creatinine clearance declined significantly during pneumoperitoneum in comparison to group A, but after deflation a faster recovery of glomerular filtration was noticed for group C, and a significant increase in sodium excretory fraction was seen for group B. On the other hand, in comparison to the control group, group B had a significant increase in plasma renin activity, with late recovery of glomerular function. Conclusion: Helium ameliorates renal alterations when used for pneumoperitoneum, and it might be used for patients with compromised renal function who have to undergo laparoscopic surgery.

An empirical evaluation of imputation accuracy for association statistics reveals increased type-I error rates in genome-wide associations

Background: Genome wide association studies (GWAS) are becoming the approach of choice to identify genetic determinants of complex phenotypes and common diseases. The astonishing amount of generated data and the use of distinct genotyping platforms with variable genomic coverage are still analytical challenges. Imputation algorithms combine directly genotyped markers information with haplotypic structure for the population of interest for the inference of a badly genotyped or missing marker and are considered a near zero cost approach to allow the comparison and combination of data generated in different studies. Several reports stated that imputed markers have an overall acceptable accuracy but no published report has performed a pair wise comparison of imputed and empiric association statistics of a complete set of GWAS markers. Results: In this report we identified a total of 73 imputed markers that yielded a nominally statistically significant association at P < 10(-5) for type 2 Diabetes Mellitus and compared them with results obtained based on empirical allelic frequencies. Interestingly, despite their overall high correlation, association statistics based on imputed frequencies were discordant in 35 of the 73 (47%) associated markers, considerably inflating the type I error rate of imputed markers. We comprehensively tested several quality thresholds, the haplotypic structure underlying imputed markers and the use of flanking markers as predictors of inaccurate association statistics derived from imputed markers. Conclusions: Our results suggest that association statistics from imputed markers showing specific MAF (Minor Allele Frequencies) range, located in weak linkage disequilibrium blocks or strongly deviating from local patterns of association are prone to have inflated false positive association signals. The present study highlights the potential of imputation procedures and proposes simple procedures for selecting the best imputed markers for follow-up genotyping studies.

Relatedness and HLA-DRB1 typing may discriminate the magnitude of the genetic susceptibility to tuberculosis using a household contact model

Background Tuberculosis clusters in families may be due to increased household exposure, shared genetic factors, or both. Household contact studies are useful to control exposure because socioeconomic and environmental conditions are similar to all subjects, allowing the evaluation of the contribution of relatedness to disease development. Methods In this study, the familial aggregation of tuberculosis using relatedness and a specific inherited marker (HLA-DRB1) was evaluated. Fifty families, which had at least two cases of tuberculosis diagnosed within the past 5 years, were selected from a cohort of tuberculosis carried out in Recife, Brazil. The first case diagnosed was considered to be a primary case. The secondary attack rate of tuberculosis in household contacts was estimated according to the degree of relatedness. The relative risk of having tuberculosis based on the degree of relatedness household and the population attributable fraction to relatedness were also estimated. HLA-DRB1 typing and attributable etiologic/preventive fractions were calculated among sick and healthy household contacts. Results Compared to unrelated contacts, the relative risk for tuberculosis adjusted for age was 1.38 (95% CI 0.86 to 2.21). Relatedness contributed 23% to the development of tuberculosis at the population levels. The HLA-DRB1*04 allele group (OR = 2.44; p =0.0324; etiologic fraction =0.15) was overrepresented and the DRB1*15 allele group (OR=0.48; p=0.0488; protective fraction=0.19) was underrepresented among household contacts exhibiting tuberculosis. The presence of DRB1 shared alleles between primary cases and their contacts was a risk factor for tuberculosis (p=0.0281). Conclusion This household contact model together with the utilisation of two genetic variables permitted the evaluation of genetic factors contributing towards tuberculosis development.

Chaotic model and memory in single calcium-activated potassium channel kinetics

Ion channels are pores formed by proteins and responsible for carrying ion fluxes through cellular membranes. The ion channels can assume conformational states thereby controlling ion flow. Physically, the conformational transitions from one state to another are associated with energy barriers between them and are dependent on stimulus, such as, electrical field, ligands, second messengers, etc. Several models have been proposed to describe the kinetics of ion channels. The classical Markovian model assumes that a future transition is independent of the time that the ion channel stayed in a previous state. Others models as the fractal and the chaotic assume that the rate of transitions between the states depend on the time that the ionic channel stayed in a previous state. For the calcium activated potassium channels of Leydig cells the R/S Hurst analysis has indicated that the channels are long-term correlated with a Hurst coefficient H around 0.7, showing a persistent memory in this kinetic. Here, we applied the R/S analysis to the opening and closing dwell time series obtained from simulated data from a chaotic model proposed by L. Liebovitch and T. Toth [J. Theor. Biol. 148, 243 (1991)] and we show that this chaotic model or any model that treats the set of channel openings and closings as independent events is inadequate to describe the long-term correlation (memory) already described for the experimental data. (C) 2008 American Institute of Physics.

Tumor slices as a model to evaluate doxorubicin in vitro treatment and expression of trios of genes PRSS11, MTSS1, CLPTM1 and PRSS11, MTSS1, SMYD2 in canine mammary gland cancer

Background: In women with breast cancer submitted to neoadjuvant chemotherapy based in doxorubicin, tumor expression of groups of three genes (PRSS11, MTSS1, CLPTM1 and PRSS11, MTSS1, SMYD2) have classified them as responsive or resistant. We have investigated whether expression of these trios of genes could predict mammary carcinoma response in dogs and whether tumor slices, which maintain epithelial-mesenchymal interactions, could be used to evaluate drug response in vitro. Methods: Tumors from 38 dogs were sliced and cultured with or without doxorubicin 1 mu M for 24 h. Tumor cells were counted by two observers to establish a percentage variation in cell number, between slices. Based on these results, a reduction in cell number between treated and control samples >= 21.7%, arbitrarily classified samples, as drug responsive. Tumor expression of PRSS11, MTSS1, CLPTM1 and SMYD2, was evaluated by real time PCR. Relative expression results were then transformed to their natural logarithm values, which were spatially disposed according to the expression of trios of genes, comprising PRSS11, MTSS1, CLPTM1 and PRSS11, MTSS1, SMYD2. Fisher linear discrimination test was used to generate a separation plane between responsive and non-responsive tumors. Results: Culture of tumor slices for 24 h was feasible. Nine samples were considered responsive and 29 non-responsive to doxorubicin, considering the pre-established cut-off value of cell number reduction = 21.7%, between doxorubicin treated and control samples. Relative gene expression was evaluated and tumor samples were then spatially distributed according to the expression of the trios of genes: PRSS11, MTSS1, CLPTM1 and PRSS11, MTSS1, SMYD2. A separation plane was generated. However, no clear separation between responsive and non-responsive samples could be observed. Conclusion: Three-dimensional distribution of samples according to the expression of the trios of genes PRSS11, MTSS1, CLPTM1 and PRSS11, MTSS1, SMYD2 could not predict doxorubicin in vitro responsiveness. Short term culture of mammary gland cancer slices may be an interesting model to evaluate chemotherapy activity.

A New Mouse Model for Marfan Syndrome Presents Phenotypic Variability Associated with the Genetic Background and Overall Levels of Fbn1 Expression

Marfan syndrome is an autosomal dominant disease of connective tissue caused by mutations in the fibrillin-1 encoding gene FBN1. Patients present cardiovascular, ocular and skeletal manifestations, and although being fully penetrant, MFS is characterized by a wide clinical variability both within and between families. Here we describe a new mouse model of MFS that recapitulates the clinical heterogeneity of the syndrome in humans. Heterozygotes for the mutant Fbn1 allele mg Delta(loxPneo), carrying the same internal deletion of exons 19-24 as the mg Delta mouse model, present defective microfibrillar deposition, emphysema, deterioration of aortic wall and kyphosis. However, the onset of a clinical phenotypes is earlier in the 129/Sv than in C57BL/6 background, indicating the existence of genetic modifiers of MFS between these two mouse strains. In addition, we characterized a wide clinical variability within the 129/Sv congenic heterozygotes, suggesting involvement of epigenetic factors in disease severity. Finally, we show a strong negative correlation between overall levels of Fbn1 expression and the severity of the phenotypes, corroborating the suggested protective role of normal fibrillin-1 in MFS pathogenesis, and supporting the development of therapies based on increasing Fbn1 expression.

Simplified Mathematical Model for an Anaerobic Sequencing Batch Biofilm Reactor Treating Lipid-Rich Wastewater Subject to Rising Organic Loading Rates

This study proposes a simplified mathematical model to describe the processes occurring in an anaerobic sequencing batch biofilm reactor (ASBBR) treating lipid-rich wastewater. The reactor, subjected to rising organic loading rates, contained biomass immobilized cubic polyurethane foam matrices, and was operated at 32 degrees C +/- 2 degrees C, using 24-h batch cycles. In the adaptation period, the reactor was fed with synthetic substrate for 46 days and was operated without agitation. Whereas agitation was raised to 500 rpm, the organic loading rate (OLR) rose from 0.3 g chemical oxygen demand (COD) . L(-1) . day(-1) to 1.2 g COD . L(-1) . day(-1). The ASBBR was fed fat-rich wastewater (dairy wastewater), in an operation period lasting for 116 days, during which four operational conditions (OCs) were tested: 1.1 +/- 0.2 g COD . L(-1) . day(-1) (OC1), 4.5 +/- 0.4 g COD . L(-1) . day(-1) (OC2), 8.0 +/- 0.8 g COD . L(-1) . day(-1) (OC3), and 12.1 +/- 2.4 g COD . L(-1) . day(-1) (OC4). The bicarbonate alkalinity (BA)/COD supplementation ratio was 1:1 at OC1, 1:2 at OC2, and 1:3 at OC3 and OC4. Total COD removal efficiencies were higher than 90%, with a constant production of bicarbonate alkalinity, in all OCs tested. After the process reached stability, temporal profiles of substrate consumption were obtained. Based on these experimental data a simplified first-order model was fit, making possible the inference of kinetic parameters. A simplified mathematical model correlating soluble COD with volatile fatty acids (VFA) was also proposed, and through it the consumption rates of intermediate products as propionic and acetic acid were inferred. Results showed that the microbial consortium worked properly and high efficiencies were obtained, even with high initial substrate concentrations, which led to the accumulation of intermediate metabolites and caused low specific consumption rates.

The Coupled Aerosol and Tracer Transport model to the Brazilian developments on the Regional Atmospheric Modeling System (CATT-BRAMS) - Part 1: Model description and evaluation

We introduce the Coupled Aerosol and Tracer Transport model to the Brazilian developments on the Regional Atmospheric Modeling System (CATT-BRAMS). CATT-BRAMS is an on-line transport model fully consistent with the simulated atmospheric dynamics. Emission sources from biomass burning and urban-industrial-vehicular activities for trace gases and from biomass burning aerosol particles are obtained from several published datasets and remote sensing information. The tracer and aerosol mass concentration prognostics include the effects of sub-grid scale turbulence in the planetary boundary layer, convective transport by shallow and deep moist convection, wet and dry deposition, and plume rise associated with vegetation fires in addition to the grid scale transport. The radiation parameterization takes into account the interaction between the simulated biomass burning aerosol particles and short and long wave radiation. The atmospheric model BRAMS is based on the Regional Atmospheric Modeling System (RAMS), with several improvements associated with cumulus convection representation, soil moisture initialization and surface scheme tuned for the tropics, among others. In this paper the CATT-BRAMS model is used to simulate carbon monoxide and particulate material (PM(2.5)) surface fluxes and atmospheric transport during the 2002 LBA field campaigns, conducted during the transition from the dry to wet season in the southwest Amazon Basin. Model evaluation is addressed with comparisons between model results and near surface, radiosondes and airborne measurements performed during the field campaign, as well as remote sensing derived products. We show the matching of emissions strengths to observed carbon monoxide in the LBA campaign. A relatively good comparison to the MOPITT data, in spite of the fact that MOPITT a priori assumptions imply several difficulties, is also obtained.

Turbulent magnetic pumping in a Babcock-Leighton solar dynamo model

Context. The turbulent pumping effect corresponds to the transport of magnetic flux due to the presence of density and turbulence gradients in convectively unstable layers. In the induction equation it appears as an advective term and for this reason it is expected to be important in the solar and stellar dynamo processes. Aims. We explore the effects of turbulent pumping in a flux-dominated Babcock-Leighton solar dynamo model with a solar-like rotation law. Methods. As a first step, only vertical pumping has been considered through the inclusion of a radial diamagnetic term in the induction equation. In the second step, a latitudinal pumping term was included and then, a near-surface shear was included. Results. The results reveal the importance of the pumping mechanism in solving current limitations in mean field dynamo modeling, such as the storage of the magnetic flux and the latitudinal distribution of the sunspots. If a meridional flow is assumed to be present only in the upper part of the convective zone, it is the full turbulent pumping that regulates both the period of the solar cycle and the latitudinal distribution of the sunspot activity. In models that consider shear near the surface, a second shell of toroidal field is generated above r = 0.95 R(circle dot) at all latitudes. If the full pumping is also included, the polar toroidal fields are efficiently advected inwards, and the toroidal magnetic activity survives only at the observed latitudes near the equator. With regard to the parity of the magnetic field, only models that combine turbulent pumping with near-surface shear always converge to the dipolar parity. Conclusions. This result suggests that, under the Babcock-Leighton approach, the equartorward motion of the observed magnetic activity is governed by the latitudinal pumping of the toroidal magnetic field rather than by a large scale coherent meridional flow. Our results support the idea that the parity problem is related to the quadrupolar imprint of the meridional flow on the poloidal component of the magnetic field and the turbulent pumping positively contributes to wash out this imprint.

uvby-beta photometry of solar twins The solar colors, model atmospheres, and the T(eff) and metallicity scales

Aims. Solar colors have been determined on the uvby-beta photometric system to test absolute solar fluxes, to examine colors predicted by model atmospheres as a function of stellar parameters (T(eff), log g, [Fe/H]), and to probe zero-points of T(eff) and metallicity scales. Methods. New uvby-beta photometry is presented for 73 solar-twin candidates. Most stars of our sample have also been observed spectroscopically to obtain accurate stellar parameters. Using the stars that most closely resemble the Sun, and complementing our data with photometry available in the literature, the solar colors on the uvby-beta system have been inferred. Our solar colors are compared with synthetic solar colors computed from absolute solar spectra and from the latest Kurucz (ATLAS9) and MARCS model atmospheres. The zero-points of different T(eff) and metallicity scales are verified and corrections are proposed. Results. Our solar colors are (b - y)(circle dot) = 0.4105 +/- 0.0015, m(1,circle dot) = 0.2122 +/- 0.0018, c(1,circle dot) = 0.3319 +/- 0.0054, and beta(circle dot) = 2.5915 +/- 0.0024. The (b - y)(circle dot) and m(1,circle dot) colors obtained from absolute spectrophotometry of the Sun agree within 3-sigma with the solar colors derived here when the photometric zero-points are determined from either the STIS HST observations of Vega or an ATLAS9 Vega model, but the c(1,circle dot) and beta(circle dot) synthetic colors inferred from absolute solar spectra agree with our solar colors only when the zero-points based on the ATLAS9 model are adopted. The Kurucz solar model provides a better fit to our observations than the MARCS model. For photometric values computed from the Kurucz models, (b - y)(circle dot) and m(1,circle dot) are in excellent agreement with our solar colors independently of the adopted zero-points, but for c(1,circle dot) and beta circle dot agreement is found only when adopting the ATLAS9 zero-points. The c(1,circle dot) color computed from both the Kurucz and MARCS models is the most discrepant, probably revealing problems either with the models or observations in the u band. The T(eff) calibration of Alonso and collaborators has the poorest performance (similar to 140 K off), while the relation of Casagrande and collaborators is the most accurate (within 10 K). We confirm that the Ramirez & Melendez uvby metallicity calibration, recommended by Arnadottir and collaborators to obtain [Fe/H] in F, G, and K dwarfs, needs a small (similar to 10%) zero-point correction to place the stars and the Sun on the same metallicity scale. Finally, we confirm that the c(1) index in solar analogs has a strong metallicity sensitivity.

Cyclic variability of the circumstellar disk of the Be star zeta Tauri II. Testing the 2D global disk oscillation model

Context. About 2/3 of the Be stars present the so-called V/R variations, a phenomenon characterized by the quasi-cyclic variation in the ratio between the violet and red emission peaks of the HI emission lines. These variations are generally explained by global oscillations in the circumstellar disk forming a one-armed spiral density pattern that precesses around the star with a period of a few years. Aims. This paper presents self-consistent models of polarimetric, photometric, spectrophotometric, and interferometric observations of the classical Be star zeta Tauri. The primary goal is to conduct a critical quantitative test of the global oscillation scenario. Methods. Detailed three-dimensional, NLTE radiative transfer calculations were carried out using the radiative transfer code HDUST. The most up-to-date research on Be stars was used as input for the code in order to include a physically realistic description for the central star and the circumstellar disk. The model adopts a rotationally deformed, gravity darkened central star, surrounded by a disk whose unperturbed state is given by a steady-state viscous decretion disk model. It is further assumed that this disk is in vertical hydrostatic equilibrium. Results. By adopting a viscous decretion disk model for zeta Tauri and a rigorous solution of the radiative transfer, a very good fit of the time-average properties of the disk was obtained. This provides strong theoretical evidence that the viscous decretion disk model is the mechanism responsible for disk formation. The global oscillation model successfully fitted spatially resolved VLTI/AMBER observations and the temporal V/R variations in the H alpha and Br gamma lines. This result convincingly demonstrates that the oscillation pattern in the disk is a one-armed spiral. Possible model shortcomings, as well as suggestions for future improvements, are also discussed.