We present the results of microsecond molecular dynamics simulations carried out by the ABC group of laboratories on a set of B-DNA oligomers containing the 136 distinct tetranucleotide base sequences. We demonstrate that the resulting trajectories have extensively sampled the conformational space accessible to B-DNA at room temperature. We confirm that base sequence effects depend strongly not only on the specific base pair step, but also on the specific base pairs that flank each step. Beyond sequence effects on average helical parameters and conformational fluctuations, we also identify tetranucleotide sequences that oscillate between several distinct conformational substates. By analyzing the conformation of the phosphodiester backbones, it is possible to understand for which sequences these substates will arise, and what impact they will have on specific helical parameters.

Bayesian analysis of macroevolutionary mixtures (BAMM) is a statistical framework that uses reversible jump Markov chain Monte Carlo to infer complex macroevolutionary dynamics of diversification and phenotypic evolution on phylogenetic trees. A recent article by Moore et al. (MEA) reported a number of theoretical and practical concerns with BAMM. Major claims from MEA are that (i) BAMM’s likelihood function is incorrect, because it does not account for unobserved rate shifts; (ii) the posterior distribution on the number of rate shifts is overly sensitive to the prior; and (iii) diversification rate estimates from BAMM are unreliable. Here, we show that these and other conclusions from MEA are generally incorrect or unjustified. We first demonstrate that MEA’s numerical assessment of the BAMM likelihood is compromised by their use of an invalid likelihood function. We then show that “unobserved rate shifts” appear to be irrelevant for biologically plausible parameterizations of the diversification process. We find that the purportedly extreme prior sensitivity reported by MEA cannot be replicated with standard usage of BAMM v2.5, or with any other version when conventional Bayesian model selection is performed. Finally, we demonstrate that BAMM performs very well at estimating diversification rate variation across the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}${\sim}$\end{document}20% of simulated trees in MEA’s data set for which it is theoretically possible to infer rate shifts with confidence. Due to ascertainment bias, the remaining 80% of their purportedly variable-rate phylogenies are statistically indistinguishable from those produced by a constant-rate birth–death process and were thus poorly suited for the summary statistics used in their performance assessment. We demonstrate that inferences about diversification rates have been accurate and consistent across all major previous releases of the BAMM software. We recognize an acute need to address the theoretical foundations of rate-shift models for phylogenetic trees, and we expect BAMM and other modeling frameworks to improve in response to mathematical and computational innovations. However, we remain optimistic that that the imperfect tools currently available to comparative biologists have provided and will continue to provide important insights into the diversification of life on Earth.

Although DNA is frequently bent and supercoiled in the cell, much of the available information on DNA structure at the atomistic level is restricted to short linear sequences. We report atomistic molecular dynamics (MD) simulations of a series of DNA minicircles containing between 65 and 110 bp which we compare with a recent biochemical study of structural distortions in these tight DNA loops. We have observed a wealth of non-canonical DNA structures such as kinks, denaturation bubbles and wrinkled conformations that form in response to bending and torsional stress. The simulations show that bending alone is sufficient to induce the formation of kinks in circles containing only 65 bp, but we did not observe any defects in simulations of larger torsionally relaxed circles containing 110 bp over the same MD timescales. We also observed that under-winding in minicircles ranging in size from 65 to 110 bp leads to the formation of single stranded bubbles and wrinkles. These calculations are used to assess the ability of atomistic MD simulations to determine the structure of bent and supercoiled DNA.

Multiple myeloma (MM) is a plasma cell malignancy with a significant heritable basis. Genome-wide association studies have transformed our understanding of MM predisposition, but individual studies have had limited power to discover risk loci. Here we perform a meta-analysis of these GWAS, add a new GWAS and perform replication analyses resulting in 9,866 cases and 239,188 controls. We confirm all nine known risk loci and discover eight new loci at 6p22.3 (rs34229995, P=1.31 × 10−8), 6q21 (rs9372120, P=9.09 × 10−15), 7q36.1 (rs7781265, P=9.71 × 10−9), 8q24.21 (rs1948915, P=4.20 × 10−11), 9p21.3 (rs2811710, P=1.72 × 10−13), 10p12.1 (rs2790457, P=1.77 × 10−8), 16q23.1 (rs7193541, P=5.00 × 10−12) and 20q13.13 (rs6066835, P=1.36 × 10−13), which localize in or near to JARID2, ATG5, SMARCD3, CCAT1, CDKN2A, WAC, RFWD3 and PREX1. These findings provide additional support for a polygenic model of MM and insight into the biological basis of tumour development.

A Monte Carlo code applied to the cgDNA coarse-grain rigid-base model of B-form double-stranded DNA is used to predict a sequence-averaged persistence length of l = 53.5 nm in the sense of Flory, and of l = 160 bp or 53.5 nm in the sense of apparent tangent-tangent correlation decay. These estimates are slightly higher than the consensus experimental values of 150 bp or 50 nm, but we believe the agreement to be good given that the cgDNA model is itself parametrized from molecular dynamics simulations of short fragments of length 10-20 bp, with no explicit fit to persistence length. Our Monte Carlo simulations further predict that there can be substantial dependence of persistence lengths on the specific sequence [Formula: see text] of a fragment. We propose, and confirm the numerical accuracy of, a simple factorization that separates the part of the apparent tangent-tangent correlation decay [Formula: see text] attributable to intrinsic shape, from a part [Formula: see text] attributable purely to stiffness, i.e., a sequence-dependent version of what has been called sequence-averaged dynamic persistence length l̅ (=58.8 nm within the cgDNA model). For ensembles of both random and λ-phage fragments, the apparent persistence length [Formula: see text] has a standard deviation of 4 nm over sequence, whereas our dynamic persistence length [Formula: see text] has a standard deviation of only 1 nm. However, there are notable dynamic persistence length outliers, including poly(A) (exceptionally straight and stiff), poly(TA) (tightly coiled and exceptionally soft), and phased A-tract sequence motifs (exceptionally bent and stiff). The results of our numerical simulations agree reasonably well with both molecular dynamics simulation and diverse experimental data including minicircle cyclization rates and stereo cryo-electron microscopy images.

Convergence is widely regarded as compelling evidence for adaptation, often being portrayed as evidence that phenotypic outcomes are predictable from ecology, overriding contingencies of history. However, repeated outcomes may be very rare unless adaptive landscapes are simple, structured by strong ecological and functional constraints. One such constraint may be a limitation on body size because performance often scales with size, allowing species to adapt to challenging functions by modifying only size. When size is constrained, species might adapt by changing shape; convergent shapes may therefore be common when size is limiting and functions are challenging. We examine the roles of size and diet as determinants of jaw shape in Sciuridae. As expected, size and diet have significant interdependent effects on jaw shape and ecomorphological convergence is rare, typically involving demanding diets and limiting sizes. More surprising is morphological without ecological convergence, which is equally common between and within dietary classes. Those cases, like rare ecomorphological convergence, may be consequences of evolving on an adaptive landscape shaped by many-to-many relationships between ecology and function, many-to-one relationships between form and performance, and one-to-many relationships between functionally versatile morphologies and ecology. On complex adaptive landscapes, ecological selection can yield different outcomes.

Genome-wide association studies (GWAS) have transformed our understanding of susceptibility to multiple myeloma (MM), but much of the heritability remains unexplained. We report a new GWAS, a meta-analysis with previous GWAS and a replication series, totalling 9974 MM cases and 247,556 controls of European ancestry. Collectively, these data provide evidence for six new MM risk loci, bringing the total number to 23. Integration of information from gene expression, epigenetic profiling and in situ Hi-C data for the 23 risk loci implicate disruption of developmental transcriptional regulators as a basis of MM susceptibility, compatible with altered B-cell differentiation as a key mechanism. Dysregulation of autophagy/apoptosis and cell cycle signalling feature as recurrently perturbed pathways. Our findings provide further insight into the biological basis of MM.

The sudden environmental catastrophe in the wake of the end-Cretaceous asteroid impact had drastic effects that rippled through animal communities. To explore how these effects may have been exacerbated by prior ecological changes, we used a food-web model to simulate the effects of primary productivity disruptions, such as those predicted to result from an asteroid impact, on ten Campanian and seven Maastrichtian terrestrial localities in North America. Our analysis documents that a shift in trophic structure between Campanian and Maastrichtian communities in North America led Maastrichtian communities to experience more secondary extinction at lower levels of primary production shutdown and possess a lower collapse threshold than Campanian communities. Of particular note is the fact that changes in dinosaur richness had a negative impact on the robustness of Maastrichtian ecosystems against environmental perturbations. Therefore, earlier ecological restructuring may have exacerbated the impact and severity of the end-Cretaceous extinction, at least in North America.