Multi-locus approaches for the measurement of selection on correlated genetic loci

The study of ecological speciation is inherently linked to the study of selection. Methods for estimating phenotypic selection within a generation based on associations between trait values and fitness (e.g., survival) of individuals are established. These methods attempt to disentangle selection acting directly on a trait from indirect selection caused by correlations with other traits via multivariate statistical approaches (i.e., inference of selection gradients). The estimation of selection on genotypic or genomic variation could also benefit from disentangling direct and indirect selection on genetic loci. However, achieving this goal is difficult with genomic data because the number of potentially correlated genetic loci (p) is very large relative to the number of individuals sampled (n). In other words, the number of model parameters exceeds the number of observations (p ≫ n). We present simulations examining the utility of whole genome regression approaches (i.e., Bayesian sparse linear mixed models) for quantifying direct selection in cases where p ≫ n. Such models have been used for genome-wide association mapping and are common in artificial breeding. Our results show they hold promise for studies of natural selection in the wild, and thus of ecological speciation. But we also demonstrate important limitations to the approach and discuss study designs required for more robust inferences

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

Like trainer, like bot? Inheritance of bias in algorithmic content moderation

The internet has become a central medium through which `networked publics' express their opinions and engage in debate. Offensive comments and personal attacks can inhibit participation in these spaces. Automated content moderation aims to overcome this problem using machine learning classifiers trained on large corpora of texts manually annotated for offence. While such systems could help encourage more civil debate, they must navigate inherently normatively contestable boundaries, and are subject to the idiosyncratic norms of the human raters who provide the training data. An important objective for platforms implementing such measures might be to ensure that they are not unduly biased towards or against particular norms of offence. This paper provides some exploratory methods by which the normative biases of algorithmic content moderation systems can be measured, by way of a case study using an existing dataset of comments labelled for offence. We train classifiers on comments labelled by different demographic subsets (men and women) to understand how differences in conceptions of offence between these groups might affect the performance of the resulting models on various test sets. We conclude by discussing some of the ethical choices facing the implementers of algorithmic moderation systems, given various desired levels of diversity of viewpoints amongst discussion participants.Comment: 12 pages, 3 figures, 9th International Conference on Social Informatics (SocInfo 2017), Oxford, UK, 13--15 September 2017 (forthcoming in Springer Lecture Notes in Computer Science

Mobility promotes and jeopardizes biodiversity in rock-paper-scissors games

Biodiversity is essential to the viability of ecological systems. Species diversity in ecosystems is promoted by cyclic, non-hierarchical interactions among competing populations. Such non-transitive relations lead to an evolution with central features represented by the `rock-paper-scissors' game, where rock crushes scissors, scissors cut paper, and paper wraps rock. In combination with spatial dispersal of static populations, this type of competition results in the stable coexistence of all species and the long-term maintenance of biodiversity. However, population mobility is a central feature of real ecosystems: animals migrate, bacteria run and tumble. Here, we observe a critical influence of mobility on species diversity. When mobility exceeds a certain value, biodiversity is jeopardized and lost. In contrast, below this critical threshold all subpopulations coexist and an entanglement of travelling spiral waves forms in the course of temporal evolution. We establish that this phenomenon is robust, it does not depend on the details of cyclic competition or spatial environment. These findings have important implications for maintenance and evolution of ecological systems and are relevant for the formation and propagation of patterns in excitable media, such as chemical kinetics or epidemic outbreaks.Comment: Final submitted version; the printed version can be found at http://dx.doi.org/10.1038/nature06095 Supplementary movies are available at http://www.theorie.physik.uni-muenchen.de/lsfrey/images_content/movie1.AVI and http://www.theorie.physik.uni-muenchen.de/lsfrey/images_content/movie2.AV

The what and where of adding channel noise to the Hodgkin-Huxley equations

One of the most celebrated successes in computational biology is the Hodgkin-Huxley framework for modeling electrically active cells. This framework, expressed through a set of differential equations, synthesizes the impact of ionic currents on a cell's voltage -- and the highly nonlinear impact of that voltage back on the currents themselves -- into the rapid push and pull of the action potential. Latter studies confirmed that these cellular dynamics are orchestrated by individual ion channels, whose conformational changes regulate the conductance of each ionic current. Thus, kinetic equations familiar from physical chemistry are the natural setting for describing conductances; for small-to-moderate numbers of channels, these will predict fluctuations in conductances and stochasticity in the resulting action potentials. At first glance, the kinetic equations provide a far more complex (and higher-dimensional) description than the original Hodgkin-Huxley equations. This has prompted more than a decade of efforts to capture channel fluctuations with noise terms added to the Hodgkin-Huxley equations. Many of these approaches, while intuitively appealing, produce quantitative errors when compared to kinetic equations; others, as only very recently demonstrated, are both accurate and relatively simple. We review what works, what doesn't, and why, seeking to build a bridge to well-established results for the deterministic Hodgkin-Huxley equations. As such, we hope that this review will speed emerging studies of how channel noise modulates electrophysiological dynamics and function. We supply user-friendly Matlab simulation code of these stochastic versions of the Hodgkin-Huxley equations on the ModelDB website (accession number 138950) and http://www.amath.washington.edu/~etsb/tutorials.html.Comment: 14 pages, 3 figures, review articl

Consequences of converting graded to action potentials upon neural information coding and energy efficiency

Information is encoded in neural circuits using both graded and action potentials, converting between them within single neurons and successive processing layers. This conversion is accompanied by information loss and a drop in energy efficiency. We investigate the biophysical causes of this loss of information and efficiency by comparing spiking neuron models, containing stochastic voltage-gated Na+ and K+ channels, with generator potential and graded potential models lacking voltage-gated Na+ channels. We identify three causes of information loss in the generator potential that are the by-product of action potential generation: (1) the voltage-gated Na+ channels necessary for action potential generation increase intrinsic noise and (2) introduce non-linearities, and (3) the finite duration of the action potential creates a ‘footprint’ in the generator potential that obscures incoming signals. These three processes reduce information rates by ~50% in generator potentials, to ~3 times that of spike trains. Both generator potentials and graded potentials consume almost an order of magnitude less energy per second than spike trains. Because of the lower information rates of generator potentials they are substantially less energy efficient than graded potentials. However, both are an order of magnitude more efficient than spike trains due to the higher energy costs and low information content of spikes, emphasizing that there is a two-fold cost of converting analogue to digital; information loss and cost inflation

The Dynamics of Health and Return Migration

In the final article in a six-part PLoS Medicine; series on Migration & Health, Anita Davies and colleagues from the International Organization for Migration (IOM) discuss the specific health risks and policy needs associated with return migratio

Replication Pauses of the Wild-Type and Mutant Mitochondrial DNA Polymerase Gamma: A Simulation Study

The activity of polymerase γ is complicated, involving both correct and incorrect DNA polymerization events, exonuclease activity, and the disassociation of the polymerase:DNA complex. Pausing of pol-γ might increase the chance of deletion and depletion of mitochondrial DNA. We have developed a stochastic simulation of pol-γ that models its activities on the level of individual nucleotides for the replication of mtDNA. This method gives us insights into the pausing of two pol-γ variants: the A467T substitution that causes PEO and Alpers syndrome, and the exonuclease deficient pol-γ (exo−) in premature aging mouse models. To measure the pausing, we analyzed simulation results for the longest time for the polymerase to move forward one nucleotide along the DNA strand. Our model of the exo− polymerase had extremely long pauses, with a 30 to 300-fold increase in the time required for the longest single forward step compared to the wild-type, while the naturally occurring A467T variant showed at most a doubling in the length of the pauses compared to the wild-type. We identified the cause of these differences in the polymerase pausing time to be the number of disassociations occurring in each forward step of the polymerase

Medical student and academic staff perceptions of role models: an analytical cross-sectional study

BACKGROUND: This study explored the associations between the perceptions of students and the perceptions of academic staff about the characteristics of clinical lecturers at the Department of Internal Medicine at Kerman University of Medical Sciences (KUMS). It also assessed what characteristics constitute a 'role model' from the point of view of students and staff. METHODS: Staff and students were questioned about the characteristics of their colleagues and lecturers, respectively. They were asked about 15 characteristics under four headings: personality, teaching skill, group working and overall performance as a role model. Associations between lecturers' characteristics were explored using Pearson correlation and characteristics were allocated into groups by partition cluster method. In addition, predictors of being a valuable lecturer were assessed using logistic regression analysis. RESULTS: Based on staff responses, the strongest association observed was between honesty and being respectful (r = 0.93, p < 0.0001). Based on student responses, the strongest association observed was between being professional and honesty (r = 0.98, p < 0.0001). None of the correlations between student and staff perceptions were significant for any characteristic. Two groups were recognized among the characteristics. group one contained those characteristics which were related to the lecturer's activity; while the second group contained characteristics that were related to the personality or teaching performance of the lecturer. The predictors of lecturer as 'role model' (i.e., perceptions of students) consisted mostly of characteristics from the first group, while the predictors of a 'role model' by fellow academic staff consisted of characteristics that were in both groups. CONCLUSION: These findings showed considerable differences between the perceptions of students about their lecturers when compared with perceptions of staff about their colleagues. Students were more concerned with the personality of their lecturers, while staff also considered their ideas and behaviors. This suggests that a more comprehensive assessment of a lecturer's performance could be obtained by taking into account feedback from both students and colleagues

Mutations in protocadherin 15 and cadherin 23 affect tip links and mechanotransduction in mammalian sensory hair cells

Immunocytochemical studies have shown that protocadherin-15 (PCDH15) and cadherin-23 (CDH23) are associated with tip links, structures thought to gate the mechanotransducer channels of hair cells in the sensory epithelia of the inner ear. The present report describes functional and structural analyses of hair cells from Pcdh15av3J (av3J), Pcdh15av6J (av6J) and Cdh23v2J (v2J) mice. The av3J and v2J mice carry point mutations that are predicted to introduce premature stop codons in the transcripts for Pcdh15 and Cdh23, respectively, and av6J mice have an in-frame deletion predicted to remove most of the 9th cadherin ectodomain from PCDH15. Severe disruption of hair-bundle morphology is observed throughout the early-postnatal cochlea in av3J/av3J and v2J/v2J mice. In contrast, only mild-to-moderate bundle disruption is evident in the av6J/av6J mice. Hair cells from av3J/av3J mice are unaffected by aminoglycosides and fail to load with [3H]-gentamicin or FM1-43, compounds that permeate the hair cell's mechanotransducer channels. In contrast, hair cells from av6J/av6J mice load with both FM1-43 and [3H]-gentamicin, and are aminoglycoside sensitive. Transducer currents can be recorded from hair cells of all three mutants but are reduced in amplitude in all mutants and have abnormal directional sensitivity in the av3J/av3J and v2J/v2J mutants. Scanning electron microscopy of early postnatal cochlear hair cells reveals tip-link like links in av6J/av6J mice, substantially reduced numbers of links in the av3J/av3J mice and virtually none in the v2J/v2J mice. Analysis of mature vestibular hair bundles reveals an absence of tip links in the av3J/av3J and v2J/v2J mice and a reduction in av6J/av6J mice. These results therefore provide genetic evidence consistent with PCDH15 and CDH23 being part of the tip-link complex and necessary for normal mechanotransduction