Ten Simple Rules for Getting Help from Online Scientific Communities

The increasing complexity of research requires scientists to work at the intersection of multiple fields and to face problems for which their formal education has not prepared them. For example, biologists with no or little background in programming are now often using complex scripts to handle the results from their experiments; vice versa, programmers wishing to enter the world of bioinformatics must know about biochemistry, genetics, and other fields. In this context, communication tools such as mailing lists, web forums, and online communities acquire increasing importance. These tools permit scientists to quickly contact people skilled in a specialized field. A question posed properly to the right online scientific community can help in solving difficult problems, often faster than screening literature or writing to publication authors. The growth of active online scientific communities, such as those listed in Table S1, demonstrates how these tools are becoming an important source of support for an increasing number of researchers. Nevertheless, making proper use of these resources is not easy. Adhering to the social norms of World Wide Web communication—loosely termed “netiquette”—is both important and non-trivial. In this article, we take inspiration from our experience on Internet-shared scientific knowledge, and from similar documents such as “Asking the Questions the Smart Way” and “Getting Answers”, to provide guidelines and suggestions on how to use online communities to solve scientific problems

Predictability of evolutionary trajectories in fitness landscapes

Experimental studies on enzyme evolution show that only a small fraction of all possible mutation trajectories are accessible to evolution. However, these experiments deal with individual enzymes and explore a tiny part of the fitness landscape. We report an exhaustive analysis of fitness landscapes constructed with an off-lattice model of protein folding where fitness is equated with robustness to misfolding. This model mimics the essential features of the interactions between amino acids, is consistent with the key paradigms of protein folding and reproduces the universal distribution of evolutionary rates among orthologous proteins. We introduce mean path divergence as a quantitative measure of the degree to which the starting and ending points determine the path of evolution in fitness landscapes. Global measures of landscape roughness are good predictors of path divergence in all studied landscapes: the mean path divergence is greater in smooth landscapes than in rough ones. The model-derived and experimental landscapes are significantly smoother than random landscapes and resemble additive landscapes perturbed with moderate amounts of noise; thus, these landscapes are substantially robust to mutation. The model landscapes show a deficit of suboptimal peaks even compared with noisy additive landscapes with similar overall roughness. We suggest that smoothness and the substantial deficit of peaks in the fitness landscapes of protein evolution are fundamental consequences of the physics of protein folding.Comment: 14 pages, 7 figure

Middleborns disadvantaged? testing birth-order effects on fitness in pre-industrial finns

Parental investment is a limited resource for which offspring compete in order to increase their own survival and reproductive success. However, parents might be selected to influence the outcome of sibling competition through differential investment. While evidence for this is widespread in egg-laying species, whether or not this may also be the case in viviparous species is more difficult to determine. We use pre-industrial Finns as our model system and an equal investment model as our null hypothesis, which predicts that (all else being equal) middleborns should be disadvantaged through competition. We found no overall evidence to suggest that middleborns in a family are disadvantaged in terms of their survival, age at first reproduction or lifetime reproductive success. However, when considering birth-order only among same-sexed siblings, first-, middle-and lastborn sons significantly differed in the number of offspring they were able to rear to adulthood, although there was no similar effect among females. Middleborn sons appeared to produce significantly less offspring than first-or lastborn sons, but they did not significantly differ from lastborn sons in the number of offspring reared to adulthood. Our results thus show that taking sex differences into account is important when modelling birth-order effects. We found clear evidence of firstborn sons being advantaged over other sons in the family, and over firstborn daughters. Therefore, our results suggest that parents invest differentially in their offspring in order to both preferentially favour particular offspring or reduce offspring inequalities arising from sibling competition

Construction of a Mean Square Error Adaptive Euler--Maruyama Method with Applications in Multilevel Monte Carlo

A formal mean square error expansion (MSE) is derived for Euler--Maruyama numerical solutions of stochastic differential equations (SDE). The error expansion is used to construct a pathwise a posteriori adaptive time stepping Euler--Maruyama method for numerical solutions of SDE, and the resulting method is incorporated into a multilevel Monte Carlo (MLMC) method for weak approximations of SDE. This gives an efficient MSE adaptive MLMC method for handling a number of low-regularity approximation problems. In low-regularity numerical example problems, the developed adaptive MLMC method is shown to outperform the uniform time stepping MLMC method by orders of magnitude, producing output whose error with high probability is bounded by TOL>0 at the near-optimal MLMC cost rate O(TOL^{-2}log(TOL)^4).Comment: 43 pages, 12 figure

How Noisy Adaptation of Neurons Shapes Interspike Interval Histograms and Correlations

Channel noise is the dominant intrinsic noise source of neurons causing variability in the timing of action potentials and interspike intervals (ISI). Slow adaptation currents are observed in many cells and strongly shape response properties of neurons. These currents are mediated by finite populations of ionic channels and may thus carry a substantial noise component. Here we study the effect of such adaptation noise on the ISI statistics of an integrate-and-fire model neuron by means of analytical techniques and extensive numerical simulations. We contrast this stochastic adaptation with the commonly studied case of a fast fluctuating current noise and a deterministic adaptation current (corresponding to an infinite population of adaptation channels). We derive analytical approximations for the ISI density and ISI serial correlation coefficient for both cases. For fast fluctuations and deterministic adaptation, the ISI density is well approximated by an inverse Gaussian (IG) and the ISI correlations are negative. In marked contrast, for stochastic adaptation, the density is more peaked and has a heavier tail than an IG density and the serial correlations are positive. A numerical study of the mixed case where both fast fluctuations and adaptation channel noise are present reveals a smooth transition between the analytically tractable limiting cases. Our conclusions are furthermore supported by numerical simulations of a biophysically more realistic Hodgkin-Huxley type model. Our results could be used to infer the dominant source of noise in neurons from their ISI statistics

Mitochondrial Control Region and microsatellite analyses on harbour porpoise (Phocoena phocoena) unravel population differentiation in the Baltic Sea and adjacent waters

The population status of the harbour porpoise (Phocoena phocoena) in the Baltic area has been a continuous matter of debate. Here we present the by far most comprehensive genetic population structure assessment to date for this region, both with regard to geographic coverage and sample size: 497 porpoise samples from North Sea, Skagerrak, Kattegat, Belt Sea, and Inner Baltic Sea were sequenced at the mitochondrial Control Region and 305 of these specimens were typed at 15 polymorphic microsatellite loci. Samples were stratified according to sample type (stranding vs. by-caught), sex, and season (breeding vs. non-breeding season). Our data provide ample evidence for a population split between the Skagerrak and the Belt Sea, with a transition zone in the Kattegat area. Among other measures, this was particularly visible in significant frequency shifts of the most abundant mitochondrial haplotypes. A particular haplotype almost absent in the North Sea was the most abundant in Belt Sea and Inner Baltic Sea. Microsatellites yielded a similar pattern (i.e., turnover in occurrence of clusters identified by STRUCTURE). Moreover, a highly significant association between microsatellite assignment and unlinked mitochondrial haplotypes further indicates a split between North Sea and Baltic porpoises. For the Inner Baltic Sea, we consistently recovered a small, but significant separation from the Belt Sea population. Despite recent arguments that separation should exceed a predefined threshold before populations shall be managed separately, we argue in favour of precautionary acknowledging the Inner Baltic porpoises as a separate management unit, which should receive particular attention, as it is threatened by various factors, in particular local fishery measures. © Springer Science+Business Media B.V. 2009

Geographic determinants of reported human Campylobacter infections in Scotland

<p><b>Background:</b> Campylobacteriosis is the leading cause of bacterial gastroenteritis in most developed countries. People are exposed to infection from contaminated food and environmental sources. However, the translation of these exposures into infection in the human population remains incompletely understood. This relationship is further complicated by differences in the presentation of cases, their investigation, identification, and reporting; thus, the actual differences in risk must be considered alongside the artefactual differences.</p> <p><b>Methods:</b> Data on 33,967 confirmed Campylobacter infections in mainland Scotland between 2000 and 2006 (inclusive) that were spatially referenced to the postcode sector level were analysed. Risk factors including the Carstairs index of social deprivation, the easting and northing of the centroid of the postcode sector, measures of livestock density by species and population density were tested in univariate screening using a non-spatial generalised linear model. The NHS Health Board of the case was included as a random effect in this final model. Subsequently, a spatial generalised linear mixed model (GLMM) was constructed and age-stratified sensitivity analysis was conducted on this model.</p> <p><b>Results:</b> The spatial GLMM included the protective effects of the Carstairs index (relative risk (RR) = 0.965, 95% Confidence intervals (CIs) = 0.959, 0.971) and population density (RR = 0.945, 95% CIs = 0.916, 0.974. Following stratification by age group, population density had a significant protective effect (RR = 0.745, 95% CIs = 0.700, 0.792) for those under 15 but not for those aged 15 and older (RR = 0.982, 95% CIs = 0.951, 1.014). Once these predictors have been taken into account three NHS Health Boards remain at significantly greater risk (Grampian, Highland and Tayside) and two at significantly lower risk (Argyll and Ayrshire and Arran).</p> <p><b>Conclusions:</b> The less deprived and children living in rural areas are at the greatest risk of being reported as a case of Campylobacter infection. However, this analysis cannot differentiate between actual risk and heterogeneities in individual reporting behaviour; nevertheless this paper has demonstrated that it is possible to explain the pattern of reported Campylobacter infections using both social and environmental predictors.</p&gt

The Composition of the Cuticular and Internal Free Fatty Acids and Alcohols from Lucilia sericata Males and Females

GC, GC–MS, and HPLC–LLSD analyses were used to identify and quantify cuticular and internal lipids in males and females of the blow-fly (Lucilia sericata). Sixteen free fatty acids, seven alcohols and cholesterol were identified and quantitatively determined in the cuticular lipids of L. sericata. Cuticular fatty acids ranged from C6 to C20 and included unsaturated entities such as 16:1n-9, 18:1n-9, 20:4n-3 and 20:5n-3. Cuticular alcohols (only saturated and even-numbered) ranged from C12 to C20 in males and C10 to C22 in females. Only one sterol was found in the cuticular lipids of both males and females. 23 free fatty acids, five alcohols and cholesterol were identified in the internal lipids. Internal fatty acids were present in large amounts—7.4 mg/g (female) and 10.1 mg/g (male). Only traces of internal alcohols (from C14 to C26 in males, from C14 to C22 in females) were found in L. sericata. Large amounts of internal cholesterol were identified in L. sericata males and females (0.49 and 0.97 mg/g of the insect body, respectively)

Probing the urea dependence of residual structure in denatured human α-lactalbumin

Backbone 15N relaxation parameters and 15N–1HN residual dipolar couplings (RDCs) have been measured for a variant of human α-lactalbumin (α-LA) in 4, 6, 8 and 10 M urea. In the α-LA variant, the eight cysteine residues in the protein have been replaced by alanines (all-Ala α-LA). This protein is a partially folded molten globule at pH 2 and has been shown previously to unfold in a stepwise non-cooperative manner on the addition of urea. 15N R2 values in some regions of all-Ala α-LA show significant exchange broadening which is reduced as the urea concentration is increased. Experimental RDC data are compared with RDCs predicted from a statistical coil model and with bulkiness, average area buried upon folding and hydrophobicity profiles in order to identify regions of non-random structure. Residues in the regions corresponding to the B, D and C-terminal 310 helices in native α-LA show R2 values and RDC data consistent with some non-random structural propensities even at high urea concentrations. Indeed, for residues 101–106 the residual structure persists in 10 M urea and the RDC data suggest that this might include the formation of a turn-like structure. The data presented here allow a detailed characterization of the non-cooperative unfolding of all-Ala α-LA at higher concentrations of denaturant and complement previous studies which focused on structural features of the molten globule which is populated at lower concentrations of denaturant

Monte Carlo Analysis of Neck Linker Extension in Kinesin Molecular Motors

Kinesin stepping is thought to involve both concerted conformational changes and diffusive movement, but the relative roles played by these two processes are not clear. The neck linker docking model is widely accepted in the field, but the remainder of the step – diffusion of the tethered head to the next binding site – is often assumed to occur rapidly with little mechanical resistance. Here, we investigate the effect of tethering by the neck linker on the diffusive movement of the kinesin head, and focus on the predicted behavior of motors with naturally or artificially extended neck linker domains. The kinesin chemomechanical cycle was modeled using a discrete-state Markov chain to describe chemical transitions. Brownian dynamics were used to model the tethered diffusion of the free head, incorporating resistive forces from the neck linker and a position-dependent microtubule binding rate. The Brownian dynamics and chemomechanical cycle were coupled to model processive runs consisting of many 8 nm steps. Three mechanical models of the neck linker were investigated: Constant Stiffness (a simple spring), Increasing Stiffness (analogous to a Worm-Like Chain), and Reflecting (negligible stiffness up to a limiting contour length). Motor velocities and run lengths from simulated paths were compared to experimental results from Kinesin-1 and a mutant containing an extended neck linker domain. When tethered by an increasingly stiff spring, the head is predicted to spend an unrealistically short amount of time within the binding zone, and extending the neck is predicted to increase both the velocity and processivity, contrary to experiments. These results suggest that the Worm-Like Chain is not an adequate model for the flexible neck linker domain. The model can be reconciled with experimental data if the neck linker is either much more compliant or much stiffer than generally assumed, or if weak kinesin-microtubule interactions stabilize the diffusing head near its binding site