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

Effective Rheology of Bubbles Moving in a Capillary Tube

We calculate the average volumetric flux versus pressure drop of bubbles moving in a single capillary tube with varying diameter, finding a square-root relation from mapping the flow equations onto that of a driven overdamped pendulum. The calculation is based on a derivation of the equation of motion of a bubble train from considering the capillary forces and the entropy production associated with the viscous flow. We also calculate the configurational probability of the positions of the bubbles.Comment: 4 pages, 1 figur

Microsatellite discovery in an insular amphibian (Grandisonia alternans) with comments on cross-species utility and the accuracy of locus identification from unassembled Illumina data

The Seychelles archipelago is unique among isolated oceanic islands because it features an endemic radiation of caecilian amphibians (Gymnophiona). In order to develop population genetics resources for this system, we identified microsatellite loci using unassembled Illumina MiSeq data generated from a genomic library of Grandisonia alternans, a species that occurs on multiple islands in the archipelago. Applying a recently described method (PALFINDER) we identified 8001 microsatellite loci that were potentially informative for population genetics analyses. Of these markers, we screened 60 loci using five individuals, directly sequenced several amplicons to confirm their identity, and then used eight loci to score allele sizes in 64 G. alternans individuals originating from five islands. A number of these individuals were sampled using non-lethal methods, demonstrating the efficacy of non-destructive molecular sampling in amphibian research. Although two loci satisfied our criteria as diploid, neutrally evolving loci with the statistical power to detect population structure, our success in identifying reliable loci was very low. Additionally, we discovered some issues with primer redundancy and differences between Illumina and Sanger sequences that suggest some Illumina-inferred loci are invalid. We investigated cross-species utility for eight loci and found most could be successfully amplified, sequenced and aligned across other species and genera of caecilians from the Seychelles. Thus, our study in part supported the validity of using PALFINDER with unassembled reads for microsatellite discovery within and across species, but importantly identified major limitations to applying this approach to small datasets (ca. 1 million reads) and loci with small tandem repeat sizes

Genetic control of Eucalyptus urophylla and E. grandis resistance to canker caused by Chrysoporthe cubensis

Chrysophorte cubensis induced canker occurs in nearly all tropical and subtropical regions where eucalypts are planted, causing losses in both wood quality and volume productivity, especially so in the warmer and more humid regions of Brazil. The wide inter and intra-specific genetic variability of resistance to canker among Eucalyptus species facilitates the selection of resistant plants. In this study, we evaluated resistance to this pathogen in five Eucalyptus grandis (G) and 15 E. urophylla (U) trees, as well as in 495 individuals from 27 progenies derived from crosses between the trees. In the field, six-months-old test seedlings were inoculated with C. cubensis. Lesion length in the xylem and bark was measured eight months later. The results demonstrated that xylem lesions could preferentially be used for the selection of resistant clones. Eight trees (7 U and 1 G) were susceptible, and the remainder (8 U and 4 G) resistant. Individual narrow and broad sense heritability estimates were 17 and 81%, respectively, thereby suggesting that canker resistance is quantitative and highly dependent on dominance and epistasis

Phosphorylation of SOS1 on tyrosine 1196 promotes its RAC GEF activity and contributes to BCR-ABL leukemogenesis

Son of Sevenless 1 (SOS1) is a dual guanine nucleotide exchange factor (GEF) that activates the small GTPases RAC and RAS. Although the molecular mechanisms of RAS GEF catalysis have been unveiled, how SOS1 acquires RAC GEF activity and what is the physio-pathological relevance of this activity is much less understood. Here we show that SOS1 is tyrosine phosphorylated on Y1196 by ABL. Phosphorylation of Y1196 controls SOS1 inter-molecular interaction, is required to promote the exchange of nucleotides on RAC in vitro and for platelet-derived growth factor (PDGF) activation of RAC- and RAC-dependent actin remodeling and cell migration. SOS1 is also phosphorylated on Y1196 by BCR-ABL in chronic myelogenous leukemic cells. Importantly, in these cells, SOS1 is required for BCR-ABL-mediated activation of RAC, cell proliferation and transformation in vitro and in a xenograft mouse model. Finally, genetic removal of Sos1 in the bone marrow-derived cells (BMDCs) from Sos1fl/flmice and infected with BCR-ABL causes a significant delay in the onset of leukemogenesis once BMDCs are injected into recipient, lethally irradiated mice. Thus, SOS1 is required for full transformation and critically contribute to the leukemogenic potential of BCR-ABL

Motives and comprehension in a public goods game with induced emotions

This study analyses the sensitivity of public goods contributions through the lens of psychological motives. We report the results of a public goods experiment in which subjects were induced with the motives of care and anger through autobiographical recall. Subjects' preferences, beliefs, and perceptions under each motive are compared with those of subjects experiencing a neutral autobiographical recall control condition. We find, but only for those subjects with the highest comprehension of the game, that care elicits significantly higher contributions than anger, with the control treatment in between. This positive influence of the care motive on unconditional giving is accounted for partly by preferences for giving and partly by the beliefs concerning greater contributions by others. Anger also affects attention to own and other's payoffs (using mouse tracking) and perceptions of the game's incentive structure (cooperative or competitive)