Effect of biogeographic history on population vulnerability in European amphibians.

The genetic diversity of populations, which contributes greatly to their adaptive potential, is negatively affected by anthropogenic habitat fragmentation and destruction. However, continental-scale losses of genetic diversity also resulted from the population expansions that followed the end of the last glaciation, an element that is rarely considered in a conservation context. We addressed this issue in a meta-analysis in which we compared the spatial patterns of vulnerability of 18 widespread European amphibians in light of phylogeographic histories (glacial refugia and postglacial routes) and anthropogenic disturbances. Conservation statuses significantly worsened with distances from refugia, particularly in the context of industrial agriculture; human population density also had a negative effect. These findings suggest that features associated with the loss of genetic diversity in post-glacial amphibian populations (such as enhanced fixation load or depressed adaptive potential) may increase their susceptibility to current threats (e.g., habitat fragmentation and pesticide use). We propose that the phylogeographic status of populations (i.e., refugial vs. post-glacial) should be considered in conservation assessments for regional and national red lists.

Coal deposition in carbonate rich shallow Lacustrine systems: The Calaf and Mequinenza sequences (Oligocene, eastern Ebro Basin, NE Spain)

Two main coal-bearing sequences developed during the Oligocene in the Tertiary Ebro Basin, the Calaf (early Oligocene) and Mequinenza (late Oligocene) coal basins. Coal deposition took place in shallow marsh-swamp-lake complexes which sometimes became closed and evolved under warm climatic conditions with fluctuating humidity. These shallow lacustrine systems are closely interrelated with the terminal parts of the distributive fluvial systems which spread from the tectonically active Ebro basin margins. Laterally extensive lignite-bearing sequences, including rather thin, lenticular autochthonous and/or hypautochthonous coal seams with high ash and sulphur contents, characterized coal deposition in the shallow lacustrine systems. Coal seam geometry, which makes them nearly subeconomic, resulted from the tectonic instability during basin margin evolution and the sometimes closed, arid conditions under which the lacustrine systems evolved. High ash and sulphur contents resulted from the inadequate isolation of peat forming environments from clastic influx and from the very low acidity and sometimes high sulphate contents of the lacustrine waters. Coal exploration in shallow lacustrine sequences similar to those described here must take into account that the spread of coal-forming environments and maxima of coal deposition are usually coincident with lake expansions and retraction or shifting of the terminal fluvial zones interrelated with the lacustrine areas.

Incompressibility of neutron-rich matter

The saturation properties of neutron-rich matter are investigated in a relativistic mean-field formalism using two accurately calibrated models: NL3 and FSUGold. The saturation properties density, binding energy per nucleon, and incompressibility coefficient are calculated as a function of the neutron-proton asymmetry α≡(N-Z)/A to all orders in α. Good agreement (at the 10% level or better) is found between these numerical calculations and analytic expansions that are given in terms of a handful of bulk parameters determined at saturation density. Using insights developed from the analytic approach and a general expression for the incompressibility coefficient of infinite neutron-rich matter, i.e., K0(α)=K0+Kτα2+ , we construct a hybrid model with values for K0 and Kτ as suggested by recent experimental findings. Whereas the hybrid model provides a better description of the measured distribution of isoscalar monopole strength in the Sn isotopes relative to both NL3 and FSUGold, it significantly underestimates the distribution of strength in 208Pb. Thus, we conclude that the incompressibility coefficient of neutron-rich matter remains an important open problem.

The role of intramolecular barriers on the glass transition of polymers: computer simulations versus mode coupling theory

We present computer simulations of a simple bead-spring model for polymer melts with intramolecular barriers. By systematically tuning the strength of the barriers, we investigate their role on the glass transition. Dynamic observables are analyzed within the framework of the mode coupling theory (MCT). Critical nonergodicity parameters, critical temperatures, and dynamic exponents are obtained from consistent fits of simulation data to MCT asymptotic laws. The so-obtained MCT λ-exponent increases from standard values for fully flexible chains to values close to the upper limit for stiff chains. In analogy with systems exhibiting higher-order MCT transitions, we suggest that the observed large λ-values arise form the interplay between two distinct mechanisms for dynamic arrest: general packing effects and polymer-specific intramolecular barriers. We compare simulation results with numerical solutions of the MCT equations for polymer systems, within the polymer reference interaction site model (PRISM) for static correlations. We verify that the approximations introduced by the PRISM are fulfilled by simulations, with the same quality for all the range of investigated barrier strength. The numerical solutions reproduce the qualitative trends of simulations for the dependence of the nonergodicity parameters and critical temperatures on the barrier strength. In particular, the increase in the barrier strength at fixed density increases the localization length and the critical temperature. However the qualitative agreement between theory and simulation breaks in the limit of stiff chains. We discuss the possible origin of this feature.

Interplay between surface anisotropy and dipolar interactions in an assembly of nanomagnets

We study the interplay between the effects of surface anisotropy and dipolar interactions in monodisperse assemblies of nanomagnets with oriented anisotropy. We derive asymptotic formulas for the assembly magnetization, taking into account temperature, applied field, core and surface anisotropy, and dipolar interparticle interactions. We find that the interplay between surface anisotropy and dipolar interactions is well described by the analytical expression of the assembly magnetization derived here: the overall sign of the product of the two parameters governing the surface and the dipolar contributions determines whether intrinsic and collective terms compete or have synergistic effects on the magnetization. This is illustrated by the magnetization curves of γ-Fe2O3 nanoparticle assemblies in the low concentration limit.

Robust Estimates of the Negative Binomial Model

We consider robust parametric procedures for univariate discrete distributions, focusing on the negative binomial model. The procedures are based on three steps: ?First, a very robust, but possibly inefficient, estimate of the model parameters is computed. ?Second, this initial model is used to identify outliers, which are then removed from the sample. ?Third, a corrected maximum likelihood estimator is computed with the remaining observations. The final estimate inherits the breakdown point (bdp) of the initial one and its efficiency can be significantly higher. Analogous procedures were proposed in [1], [2], [5] for the continuous case. A comparison of the asymptotic bias of various estimates under point contamination points out the minimum Neyman's chi-squared disparity estimate as a good choice for the initial step. Various minimum disparity estimators were explored by Lindsay [4], who showed that the minimum Neyman's chi-squared estimate has a 50% bdp under point contamination; in addition, it is asymptotically fully efficient at the model. However, the finite sample efficiency of this estimate under the uncontaminated negative binomial model is usually much lower than 100% and the bias can be strong. We show that its performance can then be greatly improved using the three step procedure outlined above. In addition, we compare the final estimate with the procedure described in

Spectral properties of the Laplacian of multiplex networks.

One of the more challenging tasks in the understanding of dynamical properties of models on top of complex networks is to capture the precise role of multiplex topologies. In a recent paper, Gómez et al. [ Phys. Rev. Lett. 110 028701 (2013)], some of the authors proposed a framework for the study of diffusion processes in such networks. Here, we extend the previous framework to deal with general configurations in several layers of networks and analyze the behavior of the spectrum of the Laplacian of the full multiplex. We derive an interesting decoupling of the problem that allow us to unravel the role played by the interconnections of the multiplex in the dynamical processes on top of them. Capitalizing on this decoupling we perform an asymptotic analysis that allow us to derive analytical expressions for the full spectrum of eigenvalues. This spectrum is used to gain insight into physical phenomena on top of multiplex, specifically, diffusion processes and synchronizability.

Empirical study of correlated survival times for recurrent events with proportional hazards margins and the effect of correlation and censoring.

Background: In longitudinal studies where subjects experience recurrent incidents over a period of time, such as respiratory infections, fever or diarrhea, statistical methods are required to take into account the within-subject correlation. Methods: For repeated events data with censored failure, the independent increment (AG), marginal (WLW) and conditional (PWP) models are three multiple failure models that generalize Cox"s proportional hazard model. In this paper, we revise the efficiency, accuracy and robustness of all three models under simulated scenarios with varying degrees of within-subject correlation, censoring levels, maximum number of possible recurrences and sample size. We also study the methods performance on a real dataset from a cohort study with bronchial obstruction. Results: We find substantial differences between methods and there is not an optimal method. AG and PWP seem to be preferable to WLW for low correlation levels but the situation reverts for high correlations. Conclusions: All methods are stable in front of censoring, worsen with increasing recurrence levels and share a bias problem which, among other consequences, makes asymptotic normal confidence intervals not fully reliable, although they are well developed theoretically.

The scope of Baker's law.

656 I. 657 II. 658 III. 660 IV. 661 V. 663 VI. 663 VII. 664 VIII. 664 665 References 665 SUMMARY: Baker's law refers to the tendency for species that establish on islands by long-distance dispersal to show an increased capacity for self-fertilization because of the advantage of self-compatibility when colonizing new habitat. Despite its intuitive appeal and broad empirical support, it has received substantial criticism over the years since it was proclaimed in the 1950s, not least because it seemed to be contradicted by the high frequency of dioecy on islands. Recent theoretical work has again questioned the generality and scope of Baker's law. Here, we attempt to discern where the idea is useful to apply and where it is not. We conclude that several of the perceived problems with Baker's law fall away when a narrower perspective is adopted on how it should be circumscribed. We emphasize that Baker's law should be read in terms of an enrichment of a capacity for uniparental reproduction in colonizing situations, rather than of high selfing rates. We suggest that Baker's law might be tested in four different contexts, which set the breadth of its scope: the colonization of oceanic islands, metapopulation dynamics with recurrent colonization, range expansions with recurrent colonization, and colonization through species invasions.

Progress in front propagation research

We review the progress in the field of front propagation in recent years. We survey many physical, biophysical and cross-disciplinary applications, including reduced-variable models of combustion flames, Reid's paradox of rapid forest range expansions, the European colonization of North America during the 19th century, the Neolithic transition in Europe from 13 000 to 5000 years ago, the description of subsistence boundaries, the formation of cultural boundaries, the spread of genetic mutations, theory and experiments on virus infections, models of cancer tumors, etc. Recent theoretical advances are unified in a single framework, encompassing very diverse systems such as those with biased random walks, distributed delays, sequential reaction and dispersion, cohabitation models, age structure and systems with several interacting species. Directions for future progress are outlined

Efficient total variation algorithm for fetal MRI reconstruction

Fetal MRI reconstruction aims at finding a high-resolution image given a small set of low-resolution images. It is usually modeled as an inverse problem where the regularization term plays a central role in the reconstruction quality. Literature has considered several regularization terms s.a. Dirichlet/Laplacian energy [1], Total Variation (TV)based energies [2,3] and more recently non-local means [4]. Although TV energies are quite attractive because of their ability in edge preservation, standard explicit steepest gradient techniques have been applied to optimize fetal-based TV energies. The main contribution of this work lies in the introduction of a well-posed TV algorithm from the point of view of convex optimization. Specifically, our proposed TV optimization algorithm for fetal reconstruction is optimal w.r.t. the asymptotic and iterative convergence speeds O(1/n(2)) and O(1/root epsilon), while existing techniques are in O(1/n) and O(1/epsilon). We apply our algorithm to (1) clinical newborn data, considered as ground truth, and (2) clinical fetal acquisitions. Our algorithm compares favorably with the literature in terms of speed and accuracy.

Mother-offspring competition promotes colonization success.

Colonization is the crucial process underlying range expansions, biological invasions, and metapopulation dynamics. Which individuals leave their natal population to colonize empty habitats is a crucial question and is presently unresolved. Dispersal is the first step in colonization. However, not all dispersing individuals are necessarily good colonizers. Indeed, in some species, the phenotype of dispersers differs depending on the selective pressures that induce dispersal. In particular, kin-based interactions, a factor driving social evolution, should induce different social response profiles in nondispersing and dispersing individuals. Kin competition (defined here as between the mother and offspring) has been proven to produce dispersers with a particular phenotype that may enhance their colonizing ability. By using the common lizard (Lacerta vivipara), we conducted a multipopulation experiment to study the effect of kin competition on dispersal and colonization success. We manipulated mother-offspring interactions, which are the most important component of kin competition in the studied species, at the family and population levels and measured the consequences on colonization success. We demonstrate that mother-offspring competition at the population level significantly influences colonization success. Increased competition at the population level enhanced the colonization rate of the largest juveniles as well as the growth and survival of the colonizers. Based on these results, we calculated that kin-induced colonization halves the extinction probability of a newly initiated population. Because interactions between relatives are likely to affect the ability of a species to track habitat modifications, kin-based dispersal should be considered in the study of invasion dynamics and metapopulation functioning.

Genomic insights into the Ixodes scapularis tick vector of Lyme disease.

Ticks transmit more pathogens to humans and animals than any other arthropod. We describe the 2.1 Gbp nuclear genome of the tick, Ixodes scapularis (Say), which vectors pathogens that cause Lyme disease, human granulocytic anaplasmosis, babesiosis and other diseases. The large genome reflects accumulation of repetitive DNA, new lineages of retro-transposons, and gene architecture patterns resembling ancient metazoans rather than pancrustaceans. Annotation of scaffolds representing ∼57% of the genome, reveals 20,486 protein-coding genes and expansions of gene families associated with tick-host interactions. We report insights from genome analyses into parasitic processes unique to ticks, including host 'questing', prolonged feeding, cuticle synthesis, blood meal concentration, novel methods of haemoglobin digestion, haem detoxification, vitellogenesis and prolonged off-host survival. We identify proteins associated with the agent of human granulocytic anaplasmosis, an emerging disease, and the encephalitis-causing Langat virus, and a population structure correlated to life-history traits and transmission of the Lyme disease agent.

Estimating extreme value cumulative distribution functions using bias-corrected kernel approaches

We propose a new kernel estimation of the cumulative distribution function based on transformation and on bias reducing techniques. We derive the optimal bandwidth that minimises the asymptotic integrated mean squared error. The simulation results show that our proposed kernel estimation improves alternative approaches when the variable has an extreme value distribution with heavy tail and the sample size is small.

Vacancy-assisted domain-growth in asymmetric binary alloys: A Monte Carlo Study

A Monte Carlo simulation study of the vacancy-assisted domain growth in asymmetric binary alloys is presented. The system is modeled using a three-state ABV Hamiltonian which includes an asymmetry term. Our simulated system is a stoichiometric two-dimensional binary alloy with a single vacancy which evolves according to the vacancy-atom exchange mechanism. We obtain that, compared to the symmetric case, the ordering process slows down dramatically. Concerning the asymptotic behavior it is algebraic and characterized by the Allen-Cahn growth exponent x51/2. The late stages of the evolution are preceded by a transient regime strongly affected by both the temperature and the degree of asymmetry of the alloy. The results are discussed and compared to those obtained for the symmetric case.