Epistasis not needed to explain low dN/dS

An important question in molecular evolution is whether an amino acid that occurs at a given position makes an independent contribution to fitness, or whether its effect depends on the state of other loci in the organism's genome, a phenomenon known as epistasis. In a recent letter to Nature, Breen et al. (2012) argued that epistasis must be "pervasive throughout protein evolution" because the observed ratio between the per-site rates of non-synonymous and synonymous substitutions (dN/dS) is much lower than would be expected in the absence of epistasis. However, when calculating the expected dN/dS ratio in the absence of epistasis, Breen et al. assumed that all amino acids observed in a protein alignment at any particular position have equal fitness. Here, we relax this unrealistic assumption and show that any dN/dS value can in principle be achieved at a site, without epistasis. Furthermore, for all nuclear and chloroplast genes in the Breen et al. dataset, we show that the observed dN/dS values and the observed patterns of amino acid diversity at each site are jointly consistent with a non-epistatic model of protein evolution.Comment: This manuscript is in response to "Epistasis as the primary factor in molecular evolution" by Breen et al. Nature 490, 535-538 (2012

Advancing Alternative Analysis: Integration of Decision Science.

Decision analysis-a systematic approach to solving complex problems-offers tools and frameworks to support decision making that are increasingly being applied to environmental challenges. Alternatives analysis is a method used in regulation and product design to identify, compare, and evaluate the safety and viability of potential substitutes for hazardous chemicals.Assess whether decision science may assist the alternatives analysis decision maker in comparing alternatives across a range of metrics.A workshop was convened that included representatives from government, academia, business, and civil society and included experts in toxicology, decision science, alternatives assessment, engineering, and law and policy. Participants were divided into two groups and prompted with targeted questions. Throughout the workshop, the groups periodically came together in plenary sessions to reflect on other groups' findings.We conclude the further incorporation of decision science into alternatives analysis would advance the ability of companies and regulators to select alternatives to harmful ingredients, and would also advance the science of decision analysis.We advance four recommendations: (1) engaging the systematic development and evaluation of decision approaches and tools; (2) using case studies to advance the integration of decision analysis into alternatives analysis; (3) supporting transdisciplinary research; and (4) supporting education and outreach efforts

First anatomical network analysis of fore- and hindlimb musculoskeletal modularity in bonobos, common chimpanzees, and humans

Studies of morphological integration and modularity, and of anatomical complexity in human evolution typically focus on skeletal tissues. Here we provide the first network analysis of the musculoskeletal anatomy of both the fore- and hindlimbs of the two species of chimpanzee and humans. Contra long-accepted ideas, network analysis reveals that the hindlimb displays a pattern opposite to that of the forelimb: Pan big toe is typically seen as more independently mobile, but humans are actually the ones that have a separate module exclusively related to its movements. Different fore- vs hindlimb patterns are also seen for anatomical network complexity (i.e., complexity in the arrangement of bones and muscles). For instance, the human hindlimb is as complex as that of chimpanzees but the human forelimb is less complex than in Pan. Importantly, in contrast to the analysis of morphological integration using morphometric approaches, network analyses do not support the prediction that forelimb and hindlimb are more dissimilar in species with functionally divergent limbs such as bipedal humans

smt: a Matlab structured matrices toolbox

We introduce the smt toolbox for Matlab. It implements optimized storage and fast arithmetics for circulant and Toeplitz matrices, and is intended to be transparent to the user and easily extensible. It also provides a set of test matrices, computation of circulant preconditioners, and two fast algorithms for Toeplitz linear systems.Comment: 19 pages, 1 figure, 1 typo corrected in the abstrac

Measuring Complex Morphological Traits with 3D Photogrammetry: A Case Study with Deer Antlers

The increasing availability of 3D-imaging technology provides new opportunities for measuring morphology. Photogrammetry enables easy 3D-data acquisition compared to conventional methods and here we assess its accuracy for measuring the size of deer antlers, a complex morphological structure. Using a proprietary photogrammetry software, we generated 3D images of antlers for 92 individuals from 29 species of cervids that vary widely in antler size and shape and used these to measure antler volume. By repeating the process, we found that the relative error averaged 8.5% of object size. Errors in converting arbitrary voxel units into real volumetric units accounted for 70% of the measurement variance and can therefore be reduced by replicating the conversion. We applied the method to clay models of known volume and found no indication of bias. The estimation was robust against variation in imaging device, distance and operator, but approximately 40 images per specimen were necessary to achieve good precision. We used the method to show that conventional measures of main-beam length are relatively poor estimators of antler volume. Using loose antlers of known weight, we also showed that the volume may be a relatively poor predictor of antler weight due to variation in bone density across species. We conclude that photogrammetry can be an efficient and accurate tool for measuring antlers, and likely many other complex morphological traits

Experimental evidence that livestock grazing intensity affects cyclic vole population regulation processes

Peer reviewedPublisher PD

How To Perform Meaningful Estimates of Genetic Effects

Although the genotype-phenotype map plays a central role both in Quantitative and Evolutionary Genetics, the formalization of a completely general and satisfactory model of genetic effects, particularly accounting for epistasis, remains a theoretical challenge. Here, we use a two-locus genetic system in simulated populations with epistasis to show the convenience of using a recently developed model, NOIA, to perform estimates of genetic effects and the decomposition of the genetic variance that are orthogonal even under deviations from the Hardy-Weinberg proportions. We develop the theory for how to use this model in interval mapping of quantitative trait loci using Halley-Knott regressions, and we analyze a real data set to illustrate the advantage of using this approach in practice. In this example, we show that departures from the Hardy-Weinberg proportions that are expected by sampling alone substantially alter the orthogonal estimates of genetic effects when other statistical models, like F2 or G2A, are used instead of NOIA. Finally, for the first time from real data, we provide estimates of functional genetic effects as sets of effects of natural allele substitutions in a particular genotype, which enriches the debate on the interpretation of genetic effects as implemented both in functional and in statistical models. We also discuss further implementations leading to a completely general genotype-phenotype map

Retrivability in The Danish National Hospital Registry of HIV and hepatitis B and C coinfection diagnoses of patients managed in HIV centers 1995–2004

<p>Abstract</p> <p>Background</p> <p>Hospital-based discharge registries are used increasingly for longitudinal epidemiological studies of HIV. We examined completeness of registration of HIV infections and of chronic hepatitis B (HBV) and hepatitis C (HCV) coinfections in the Danish National Hospital Registry (DNHR) covering all Danish hospitals.</p> <p>Methods</p> <p>The Danish HIV Cohort Study (DHCS) encompasses all HIV-infected patients treated in Danish HIV clinics since 1 January 1995. All 2,033 Danish patients in DHCS diagnosed with HIV-1 during the 10-year period from 1 January 1995 to 31 December 2004 were included in the current analysis. We used the DHCS as a reference to examine the completeness of HIV and of HBV and HCV coinfections recorded in DNHR. Cox regression analysis was used to estimate hazard ratios of time to diagnosis of HIV in DNHR compared to DHCS.</p> <p>Results</p> <p>Of the 2,033 HIV patients in DHCS, a total of 2,006 (99%) were registered with HIV in DNHR. Of these, 1,888 (93%) were registered in DNHR within one year of their first positive HIV test. A CD4 < 200 cells/μl, a viral load >= 100,000 copies/ml and being diagnosed after 1 January 2000, were associated with earlier registration in DNHR, both in crude and adjusted analyses. Thirty (23%) HIV patients registered with chronic HBV (n = 129) in DHCS and 126 (48%) of HIV patients with HCV (n = 264) in DHCS were registered with these diagnoses in the DNHR. Further 17 and 8 patients were registered with HBV and HCV respectively in DNHR, but not in DHCS. The positive predictive values of being registered with HBV and HCV in DHCS were thereby estimated to 0.88 and 0.97 and in DNHR to 0.32 and 0.54.</p> <p>Conclusion</p> <p>The study demonstrates that secondary data from national hospital databases may be reliable for identification of patients diagnosed with HIV infection. However, the predictive value of co-morbidity data may be low.</p

How predation and landscape fragmentation affect vole population dynamics

Background: Microtine species in Fennoscandia display a distinct north-south gradient from regular cycles to stable populations. The gradient has often been attributed to changes in the interactions between microtines and their predators. Although the spatial structure of the environment is known to influence predator-prey dynamics of a wide range of species, it has scarcely been considered in relation to the Fennoscandian gradient. Furthermore, the length of microtine breeding season also displays a north-south gradient. However, little consideration has been given to its role in shaping or generating population cycles. Because these factors covary along the gradient it is difficult to distinguish their effects experimentally in the field. The distinction is here attempted using realistic agent-based modelling. Methodology/Principal Findings: By using a spatially explicit computer simulation model based on behavioural and ecological data from the field vole (Microtus agrestis), we generated a number of repeated time series of vole densities whose mean population size and amplitude were measured. Subsequently, these time series were subjected to statistical autoregressive modelling, to investigate the effects on vole population dynamics of making predators more specialised, of altering the breeding season, and increasing the level of habitat fragmentation. We found that fragmentation as well as the presence of specialist predators are necessary for the occurrence of population cycles. Habitat fragmentation and predator assembly jointly determined cycle length and amplitude. Length of vole breeding season had little impact on the oscillations. Significance: There is good agreement between our results and the experimental work from Fennoscandia, but our results allow distinction of causation that is hard to unravel in field experiments. We hope our results will help understand the reasons for cycle gradients observed in other areas. Our results clearly demonstrate the importance of landscape fragmentation for population cycling and we recommend that the degree of fragmentation be more fully considered in future analyses of vole dynamics

The association between blood glucose and oxidized lipoprotein(a) in healthy young women

<p>Abstract</p> <p>Background</p> <p>Oxidized lipoproteins play important roles in the atherosclerotic processes. Oxidized lipoprotein(a) (oxLp(a)) may be more potent in atherosclerotic pathophysiology than native Lp(a), a cardiovascular disease-relevant lipoprotein. Increased blood glucose concentrations can induce oxidative modification of lipoproteins. The aim of this study was to investigate the association between circulating oxLp(a) and cardiometabolic variables including blood glucose in healthy volunteers within the normal range of blood glucose.</p> <p>Methods</p> <p>Several cardiometabolic variables and serum oxLp(a) (using an ELISA system) were measured among 70 healthy females (mean age, 22 years).</p> <p>Results</p> <p>Lp(a) and glucose were significantly and positively correlated with oxLp(a) in simple correlation test. Furthermore, a multiple linear regression analysis showed oxLp(a) to have a weakly, but significantly positive and independent correlation with only blood glucose (<it>β </it>= 0.269, <it>P </it>< 0.05).</p> <p>Conclusions</p> <p>These results suggest that increased glucose may enhance the oxidization of Lp(a) even at normal glucose levels.</p