Simple, Fast and Accurate Implementation of the Diffusion Approximation Algorithm for Stochastic Ion Channels with Multiple States

The phenomena that emerge from the interaction of the stochastic opening and closing of ion channels (channel noise) with the non-linear neural dynamics are essential to our understanding of the operation of the nervous system. The effects that channel noise can have on neural dynamics are generally studied using numerical simulations of stochastic models. Algorithms based on discrete Markov Chains (MC) seem to be the most reliable and trustworthy, but even optimized algorithms come with a non-negligible computational cost. Diffusion Approximation (DA) methods use Stochastic Differential Equations (SDE) to approximate the behavior of a number of MCs, considerably speeding up simulation times. However, model comparisons have suggested that DA methods did not lead to the same results as in MC modeling in terms of channel noise statistics and effects on excitability. Recently, it was shown that the difference arose because MCs were modeled with coupled activation subunits, while the DA was modeled using uncoupled activation subunits. Implementations of DA with coupled subunits, in the context of a specific kinetic scheme, yielded similar results to MC. However, it remained unclear how to generalize these implementations to different kinetic schemes, or whether they were faster than MC algorithms. Additionally, a steady state approximation was used for the stochastic terms, which, as we show here, can introduce significant inaccuracies. We derived the SDE explicitly for any given ion channel kinetic scheme. The resulting generic equations were surprisingly simple and interpretable - allowing an easy and efficient DA implementation. The algorithm was tested in a voltage clamp simulation and in two different current clamp simulations, yielding the same results as MC modeling. Also, the simulation efficiency of this DA method demonstrated considerable superiority over MC methods.Comment: 32 text pages, 10 figures, 1 supplementary text + figur

Acceleration of generalized hypergeometric functions through precise remainder asymptotics

We express the asymptotics of the remainders of the partial sums {s_n} of the generalized hypergeometric function q+1_F_q through an inverse power series z^n n^l \sum_k c_k/n^k, where the exponent l and the asymptotic coefficients {c_k} may be recursively computed to any desired order from the hypergeometric parameters and argument. From this we derive a new series acceleration technique that can be applied to any such function, even with complex parameters and at the branch point z=1. For moderate parameters (up to approximately ten) a C implementation at fixed precision is very effective at computing these functions; for larger parameters an implementation in higher than machine precision would be needed. Even for larger parameters, however, our C implementation is able to correctly determine whether or not it has converged; and when it converges, its estimate of its error is accurate.Comment: 36 pages, 6 figures, LaTeX2e. Fixed sign error in Eq. (2.28), added several references, added comparison to other methods, and added discussion of recursion stabilit

IL1RL1 Gene Variants and Nasopharyngeal IL1RL-a Levels Are Associated with Severe RSV Bronchiolitis: A Multicenter Cohort Study

Targets for intervention are required for respiratory syncytial virus (RSV) bronchiolitis, a common disease during infancy for which no effective treatment exists. Clinical and genetic studies indicate that IL1RL1 plays an important role in the development and exacerbations of asthma. Human IL1RL1 encodes three isoforms, including soluble IL1RL1-a, that can influence IL33 signalling by modifying inflammatory responses to epithelial damage. We hypothesized that IL1RL1 gene variants and soluble IL1RL1-a are associated with severe RSV bronchiolitis.We studied the association between RSV and 3 selected IL1RL1 single-nucleotide polymorphisms rs1921622, rs11685480 or rs1420101 in 81 ventilated and 384 non-ventilated children under 1 year of age hospitalized with primary RSV bronchiolitis in comparison to 930 healthy controls. Severe RSV infection was defined by need for mechanical ventilation. Furthermore, we examined soluble IL1RL1-a concentration in nasopharyngeal aspirates from children hospitalized with primary RSV bronchiolitis. An association between SNP rs1921622 and disease severity was found at the allele and genotype level (p = 0.011 and p = 0.040, respectively). In hospitalized non-ventilated patients, RSV bronchiolitis was not associated with IL1RL1 genotypes. Median concentrations of soluble IL1RL1-a in nasopharyngeal aspirates were >20-fold higher in ventilated infants when compared to non-ventilated infants with RSV (median [and quartiles] 9,357 [936-15,528] pg/ml vs. 405 [112-1,193] pg/ml respectively; p<0.001).We found a genetic link between rs1921622 IL1RL1 polymorphism and disease severity in RSV bronchiolitis. The potential biological role of IL1RL1 in the pathogenesis of severe RSV bronchiolitis was further supported by high local concentrations of IL1RL1 in children with most severe disease. We speculate that IL1RL1a modifies epithelial damage mediated inflammatory responses during RSV bronchiolitis and thus may serve as a novel target for intervention to control disease severity

Do Larval Supply and Recruitment Vary among Chemosynthetic Environments of the Deep Sea?

BACKGROUND: The biological communities that inhabit chemosynthetic environments exist in an ephemeral and patchily distributed habitat with unique physicochemical properties that lead to high endemicity. Consequently, the maintenance and recovery from perturbation of the populations in these habitats is, arguably, mainly regulated by larval supply and recruitment. METHODOLOGY/PRINCIPAL FINDINGS: WE USE DATA FROM THE PUBLISHED SCIENTIFIC LITERATURE TO: (1) compare the magnitudes of and variability in larval supply and settlement and recruitment at hydrothermal vents, seeps, and whale, wood and kelp falls; (2) explore factors that affect these life history processes, when information is available; and (3) explore taxonomic affinities in the recruit assemblages of the different chemosynthetic habitats, using multivariate statistical techniques. Larval supply at vents can vary across segments by several orders of magnitude for gastropods; for bivalves, supply is similar at vents on different segments, and at cold seeps. The limited information on larval development suggests that dispersal potential may be highest for molluscs from cold seeps, intermediate for siboglinids at vents and lowest for the whale-bone siboglinid Osedax. Settlement is poorly studied and only at vents and seeps, but tends to be highest near an active source of emanating fluid in both habitats. Rate of recruitment at vents is more variable among studies within a segment than among segments. Across different chemosynthetic habitats, recruitment rate of bivalves is much more variable than that of gastropods and polychaetes. Total recruitment rate ranges only between 0.1 and 1 ind dm(-2) d(-1) across all chemosynthetic habitats, falling above rates in the non-reducing deep sea. The recruit assemblages at vents, seeps and kelp falls have lower taxonomic breadth, and include more families and genera that have many species more closely related to each other than those at whale and wood falls. Vents also have the most uneven taxonomic structure, with fewer recruits represented by higher taxonomic levels (phyla, orders, classes) compared to seeps and wood and kelp falls, whereas the opposite is true at whale falls. CONCLUSIONS/SIGNIFICANCE: Based on our evaluation of the literature, the patterns and regulatory factors of the early history processes in chemosynthetic environments in the deep sea remain poorly understood. More research focused on these early life history stages will allow us to make inferences about the ecological and biogeographic linkages among the reducing habitats in the deep sea

Lévy patterns in seabirds are multifaceted describing both spatial and temporal patterning

BACKGROUND: The flight patterns of albatrosses and shearwaters have become a touchstone for much of Lévy flight research, spawning an extensive field of enquiry. There is now compelling evidence that the flight patterns of these seabirds would have been appreciated by Paul Lévy, the mathematician after whom Lévy flights are named. Here we show that Lévy patterns (here taken to mean spatial or temporal patterns characterized by distributions with power-law tails) are, in fact, multifaceted in shearwaters being evident in both spatial and temporal patterns of activity. RESULTS: We tested for Lévy patterns in the at-sea behaviours of two species of shearwater breeding in the North Atlantic Ocean (Calonectris borealis) and the Mediterranean sea (C. diomedea) during their incubating and chick-provisioning periods. We found that distributions of flight durations, on/in water durations and inter-dive time-intervals have power-law tails and so bear the hallmarks of Lévy patterns. CONCLUSIONS: The occurrence of these statistical laws is remarkable given that bird behaviours are strongly shaped by an individual’s motivational state and by complex environmental interactions. Our observations could take Lévy patterns as models of animal behaviour to a new level by going beyond the characterisation of spatial movements to characterise how different behaviours are interwoven throughout daily animal life

The nuclear receptors of Biomphalaria glabrata and Lottia gigantea: Implications for developing new model organisms

© 2015 Kaur et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are creditedNuclear receptors (NRs) are transcription regulators involved in an array of diverse physiological functions including key roles in endocrine and metabolic function. The aim of this study was to identify nuclear receptors in the fully sequenced genome of the gastropod snail, Biomphalaria glabrata, intermediate host for Schistosoma mansoni and compare these to known vertebrate NRs, with a view to assessing the snail's potential as a invertebrate model organism for endocrine function, both as a prospective new test organism and to elucidate the fundamental genetic and mechanistic causes of disease. For comparative purposes, the genome of a second gastropod, the owl limpet, Lottia gigantea was also investigated for nuclear receptors. Thirty-nine and thirty-three putative NRs were identified from the B. glabrata and L. gigantea genomes respectively, based on the presence of a conserved DNA-binding domain and/or ligand-binding domain. Nuclear receptor transcript expression was confirmed and sequences were subjected to a comparative phylogenetic analysis, which demonstrated that these molluscs have representatives of all the major NR subfamilies (1-6). Many of the identified NRs are conserved between vertebrates and invertebrates, however differences exist, most notably, the absence of receptors of Group 3C, which includes some of the vertebrate endocrine hormone targets. The mollusc genomes also contain NR homologues that are present in insects and nematodes but not in vertebrates, such as Group 1J (HR48/DAF12/HR96). The identification of many shared receptors between humans and molluscs indicates the potential for molluscs as model organisms; however the absence of several steroid hormone receptors indicates snail endocrine systems are fundamentally different.The National Centre for the Replacement, Refinement and Reduction of Animals in Research, Grant Ref:G0900802 to CSJ, LRN, SJ & EJR [www.nc3rs.org.uk]

Occupational cancer in Britain: Statistical methodology

An approach using the attributable fraction (AF) has been developed to estimate the current burden of occupational cancer in Britain. The AF combines the relative risk (RR) associated with exposure with the proportion exposed. For each cancer–exposure pairing, the RR is selected from key epidemiological literature such as an industry, or population-based study, meta-analysis or review. The CARcinogen EXposure (CAREX) database provides point estimates for the number of workers exposed to a range of carcinogens; alternative sources are national surveys such as the Labour Force Survey and Census of Employment. The number of workers exposed are split between high and low exposure levels matched to appropriate RRs from the literature. The relevant period for cancer development during which exposure occurred is defined as the risk exposure period (REP). Estimation of the numbers ever exposed over the REP takes into account the changes in the number of people employed in primary and manufacturing industry and service sectors in Britain where appropriate, and adjustment is made for staff turnover over the period and for life expectancy. National estimates of the population ever of working age during the REP are used for the proportion denominator. Strategies have been developed to combine exposure AFs correctly while avoiding double counting and minimising bias. The AFs are applied to national cancer deaths and registrations to obtain occupation-attributable cancer numbers. The methods are adaptable for other diseases and other geographies, and are also adaptable to more sophisticated modelling if better exposure and dose–response data are available

Tropical-cyclone-driven erosion of the terrestrial biosphere from mountains

The transfer of organic carbon from the terrestrial biosphere to the oceans via erosion and riverine transport constitutes an important component of the global carbon cycle. More than one third of this organic carbon flux comes from sediment-laden rivers that drain the mountains in the western Pacific region. This region is prone to tropical cyclones, but their role in sourcing and transferring vegetation and soil is not well constrained. Here we measure particulate organic carbon load and composition in the LiWu River, Taiwan, during cyclone-triggered floods. We correct for fossil particulate organic carbon using radiocarbon, and find that the concentration of particulate organic carbon from vegetation and soils is positively correlated with water discharge. Floods have been shown to carry large amounts of clastic sediment. Non-fossil particulate organic carbon transported at the same time may be buried offshore under high rates of sediment accumulation. We estimate that on decadal timescales, 77–92% of non-fossil particulate organic carbon eroded from the LiWu catchment is transported during large, cyclone-induced floods. We suggest that tropical cyclones, which affect many forested mountains within the Intertropical Convergence Zone, may provide optimum conditions for the delivery and burial of non-fossil particulate organic carbon in the ocean. This carbon transfer is moderated by the frequency, intensity and duration of tropical cyclones