Ignoring heterozygous sites biases phylogenomic estimates of divergence times: implications for the evolutionary history of Microtus voles

Phylogenetic reconstruction of the evolutionary history of closely related organisms may be difficult because of the presence of unsorted lineages and of a relatively high proportion of heterozygous sites that are usually not handled well by phylogenetic programs. Genomic data may provide enough fixed polymorphisms to resolve phylogenetic trees, but the diploid nature of sequence data remains analytically challenging. Here, we performed a phylogenomic reconstruction of the evolutionary history of the common vole (Microtus arvalis) with a focus on the influence of heterozygosity on the estimation of intraspecific divergence times. We used genome-wide sequence information from 15 voles distributed across the European range. We provide a novel approach to integrate heterozygous information in existing phylogenetic programs by repeated random haplotype sampling from sequences with multiple unphased heterozygous sites. We evaluated the impact of the use of full, partial, or no heterozygous information for tree reconstructions on divergence time estimates. All results consistently showed four deep and strongly supported evolutionary lineages in the vole data. These lineages undergoing divergence processes split only at the end or after the last glacial maximum based on calibration with radiocarbon-dated paleontological material. However, the incorporation of information from heterozygous sites had a significant impact on absolute and relative branch length estimations. Ignoring heterozygous information led to an overestimation of divergence times between the evolutionary lineages of M. arvalis. We conclude that the exclusion of heterozygous sites from evolutionary analyses may cause biased and misleading divergence time estimates in closely related taxa.

Models of hybridization during range expansions and their application to recent human evolution

Several lines of genetic, archeological and paleontological evidence suggest that anatomically modern humans (Homo sapiens) colonized the world in the last 60,000 years by a series of migrations originating from Africa (e.g. Liu et al., 2006; Handley et al., 2007; Prugnolle, Manica, and Balloux, 2005; Ramachandran et al. 2005; Li et al. 2008; Deshpande et al. 2009; Mellars, 2006a, b; Lahr and Foley, 1998; Gravel et al., 2011; Rasmussen et al., 2011). With the progress of ancient DNA analysis, it has been shown that archaic humans hybridized with modern humans outside Africa. Recent direct analyses of fossil nuclear DNA have revealed that 1–4 percent of the genome of Eurasian has been likely introgressed by Neanderthal genes (Green et al., 2010; Reich et al., 2010; Vernot and Akey, 2014; Sankararaman et al., 2014; Prufer et al., 2014; Wall et al., 2013), with Papua New Guineans and Australians showing even larger levels of admixture with Denisovans (Reich et al., 2010; Skoglund and Jakobsson, 2011; Reich et al., 2011; Rasmussen et al., 2011). It thus appears that the past history of our species has been more complex than previously anticipated (Alves et al., 2012), and that modern humans hybridized several times with local hominins during their expansion out of Africa, but the exact mode, time and location of these hybridizations remain to be clarifi ed (Ibid.; Wall et al., 2013). In this context, we review here a general model of admixture during range expansion, which lead to some predictions about expected patterns of introgression that are relevant to modern human evolution.

Consistency of Concave Regression with an Application to Current-Status Data

We consider the problem of nonparametric estimation of a concave regression function F. We show that the supremum distance between the least square s estimatorand F on a compact interval is typically of order(log(n)/n)2/5. This entails rates of convergence for the estimator’s derivative. Moreover, we discuss the impact of additional constraints on F such as monotonicity and pointwise bounds. Then we apply these results to the analysis of current status data, where the distribution function of the event times is assumed to be concave.

Reduction of cement leakage by sequential PMMA application in a vertebroplasty model

PURPOSE Leakage is the most common complication of percutaneous cement augmentation of the spine. The viscosity of the polymethylmethacrylate (PMMA) cement is strongly correlated with the likelihood of cement leakage. We hypothesized that cement leakage can be reduced by sequential cement injection in a vertebroplasty model. METHODS A standardized vertebral body substitute model, consisting of aluminum oxide foams coated by acrylic cement with a preformed leakage path, simulating a ventral vein, was developed. Three injection techniques of 6 ml PMMA were assessed: injection in one single step (all-in-one), injection of 1 ml at the first and 5 ml at the second step with 1 min latency in-between (two-step), and sequential injection of 0.5 ml with 1-min latency between the sequences (sequential). Standard PMMA vertebroplasty cement was used; each injection type was tested on ten vertebral body substitute models with two possible leakage paths per model. Leakage was assessed by radiographs using a zonal graduation: intraspongious = no leakage and extracortical = leakage. RESULTS The leakage rate was significantly lower in the "sequential" technique (2/20 leakages) followed by "two-step" (15/20) and "all-in-one" (20/20) techniques (p < 0.001). The RR for a cement leakage was 10.0 times higher in the "all-in-one" compared to the "sequential" group (95 % confidence intervals 2.7-37.2; p < 0.001). CONCLUSIONS The sequential cement injection is a simple approach to minimize the risk for leakage. Taking advantage of the temperature gradient between body and room temperature, it is possible to increase the cement viscosity inside the vertebra while keeping it low in the syringe. Using sequential injection of small cement volumes, further leakage paths are blocked before further injection of the low-viscosity cement.

European summer temperatures since Roman times

The spatial context is critical when assessing present-day climate anomalies, attributing them to potential forcings and making statements regarding their frequency and severity in a long-term perspective. Recent international initiatives have expanded the number of high-quality proxy-records and developed new statistical reconstruction methods. These advances allow more rigorous regional past temperature reconstructions and, in turn, the possibility of evaluating climate models on policy-relevant, spatio-temporal scales. Here we provide a new proxy-based, annually-resolved, spatial reconstruction of the European summer (June–August) temperature fields back to 755 CE based on Bayesian hierarchical modelling (BHM), together with estimates of the European mean temperature variation since 138 BCE based on BHM and composite-plus-scaling (CPS). Our reconstructions compare well with independent instrumental and proxy-based temperature estimates, but suggest a larger amplitude in summer temperature variability than previously reported. Both CPS and BHM reconstructions indicate that the mean 20th century European summer temperature was not significantly different from some earlier centuries, including the 1st, 2nd, 8th and 10th centuries CE. The 1st century (in BHM also the 10th century) may even have been slightly warmer than the 20th century, but the difference is not statistically significant. Comparing each 50 yr period with the 1951–2000 period reveals a similar pattern. Recent summers, however, have been unusually warm in the context of the last two millennia and there are no 30 yr periods in either reconstruction that exceed the mean average European summer temperature of the last 3 decades (1986–2015 CE). A comparison with an ensemble of climate model simulations suggests that the reconstructed European summer temperature variability over the period 850–2000 CE reflects changes in both internal variability and external forcing on multi-decadal time-scales. For pan-European temperatures we find slightly better agreement between the reconstruction and the model simulations with high-end estimates for total solar irradiance. Temperature differences between the medieval period, the recent period and the Little Ice Age are larger in the reconstructions than the simulations. This may indicate inflated variability of the reconstructions, a lack of sensitivity and processes to changes in external forcing on the simulated European climate and/or an underestimation of internal variability on centennial and longer time scales.