185 research outputs found
Trends to equilibrium for a class of relativistic diffusions
We address the question of the trends to equilibrium for a large class C of
relativistic diffusions. We show the existence of a spectral gap using the
Lyapounov method and deduce the exponential decay of the distance to
equilibrium in L2-norm and in total variation. A similar result was obtained
recently in arXiv:1009.5086 for a particular process of the class C.Comment: 10 page
Rphylopars: fast multivariate phylogenetic comparative methods for missing data and within-species variation
Over the past several years, phylogenetic comparative studies have increasingly approached trait evolution in a multivariate context, with a number of taxa that continues to rise dramatically. Recent methods for phylogenetic comparative studies have provided ways to incorporate measurement error and to address computational challenges. However, missing data remain a particularly common problem, in which data are unavailable for some but not all traits of interest for a given species (or individual), leaving researchers with the choice between omitting observations or utilizing imputation-based approaches.
Here, we introduce an r implementation of PhyloPars, a tool for phylogenetic imputation of missing data and estimation of trait covariance across species (phylogenetic covariance) and within species (phenotypic covariance). Rphylopars provides expanded capabilities over the original PhyloPars interface including a fast linear-time algorithm, thus allowing for extremely large data sets (which were previously computationally infeasible) to be analysed in seconds or minutes rather than hours.
In addition to providing fast and computationally efficient implementations, we introduce in Rphylopars methods to estimate macroevolutionary parameters under alternative evolutionary models (e.g. Early-Burst, multivariate Ornstein-Uhlenbeck).
By providing fast and computationally efficient methods with flexible options for various phylogenetic comparative approaches, Rphylopars expands the possibilities for researchers to analyse large and complex data with missing observations, within-species variation and deviations from Brownian motion
Dimension dependent hypercontractivity for Gaussian kernels
We derive sharp, local and dimension dependent hypercontractive bounds on the
Markov kernel of a large class of diffusion semigroups. Unlike the dimension
free ones, they capture refined properties of Markov kernels, such as trace
estimates. They imply classical bounds on the Ornstein-Uhlenbeck semigroup and
a dimensional and refined (transportation) Talagrand inequality when applied to
the Hamilton-Jacobi equation. Hypercontractive bounds on the Ornstein-Uhlenbeck
semigroup driven by a non-diffusive L\'evy semigroup are also investigated.
Curvature-dimension criteria are the main tool in the analysis.Comment: 24 page
Effective dynamics using conditional expectations
The question of coarse-graining is ubiquitous in molecular dynamics. In this
article, we are interested in deriving effective properties for the dynamics of
a coarse-grained variable , where describes the configuration of
the system in a high-dimensional space , and is a smooth function
with value in (typically a reaction coordinate). It is well known that,
given a Boltzmann-Gibbs distribution on , the equilibrium
properties on are completely determined by the free energy. On the
other hand, the question of the effective dynamics on is much more
difficult to address. Starting from an overdamped Langevin equation on , we propose an effective dynamics for using conditional
expectations. Using entropy methods, we give sufficient conditions for the time
marginals of the effective dynamics to be close to the original ones. We check
numerically on some toy examples that these sufficient conditions yield an
effective dynamics which accurately reproduces the residence times in the
potential energy wells. We also discuss the accuracy of the effective dynamics
in a pathwise sense, and the relevance of the free energy to build a
coarse-grained dynamics
A Malvaceae mystery: A mallow maelstrom of genome multiplications and maybe misleading methods?
Previous research suggests that Gossypium has undergone a 5- to 6-fold multiplication following its divergence from Theobroma. However, the number of events, or where they occurred in the Malvaceae phylogeny remains unknown. We analyzed transcriptomic and genomic data from representatives of eight of the nine Malvaceae subfamilies. Phylogenetic analysis of nuclear data placed Dombeya (Dombeyoideae) as sister to the rest of Malvadendrina clade, but the plastid DNA tree strongly supported Durio (Helicteroideae) in this position. Intraspecific Ks plots indicated that all sampled taxa, except Theobroma (Byttnerioideae), Corchorus (Grewioideae), and Dombeya (Dombeyoideae), have experienced whole genome multiplications (WGMs). Quartet analysis suggested WGMs were shared by Malvoideae-Bombacoideae and Sterculioideae-Tilioideae, but did not resolve whether these are shared with each other or Helicteroideae (Durio). Gene tree reconciliation and Bayesian concordance analysis suggested a complex history. Alternative hypotheses are suggested, each involving two independent autotetraploid and one allopolyploid event. They differ in that one entails an allopolyploid origin for the Durio lineage, whereas the other invokes an allopolyploid origin for Malvoideae-Bombacoideae. We highlight the need for more genomic information in the Malvaceae and improved methods to resolve complex evolutionary histories that may include allopolyploidy, incomplete lineage sorting, and variable rates of gene and genome evolution
The Probability of a Gene Tree Topology within a Phylogenetic Network with Applications to Hybridization Detection
Gene tree topologies have proven a powerful data source for various tasks, including species tree inference and species delimitation. Consequently, methods for computing probabilities of gene trees within species trees have been developed and widely used in probabilistic inference frameworks. All these methods assume an underlying multispecies coalescent model. However, when reticulate evolutionary events such as hybridization occur, these methods are inadequate, as they do not account for such events. Methods that account for both hybridization and deep coalescence in computing the probability of a gene tree topology currently exist for very limited cases. However, no such methods exist for general cases, owing primarily to the fact that it is currently unknown how to compute the probability of a gene tree topology within the branches of a phylogenetic network. Here we present a novel method for computing the probability of gene tree topologies on phylogenetic networks and demonstrate its application to the inference of hybridization in the presence of incomplete lineage sorting. We reanalyze a Saccharomyces species data set for which multiple analyses had converged on a species tree candidate. Using our method, though, we show that an evolutionary hypothesis involving hybridization in this group has better support than one of strict divergence. A similar reanalysis on a group of three Drosophila species shows that the data is consistent with hybridization. Further, using extensive simulation studies, we demonstrate the power of gene tree topologies at obtaining accurate estimates of branch lengths and hybridization probabilities of a given phylogenetic network. Finally, we discuss identifiability issues with detecting hybridization, particularly in cases that involve extinction or incomplete sampling of taxa
On Poincare and logarithmic Sobolev inequalities for a class of singular Gibbs measures
This note, mostly expository, is devoted to Poincar{\'e} and log-Sobolev
inequalities for a class of Boltzmann-Gibbs measures with singular interaction.
Such measures allow to model one-dimensional particles with confinement and
singular pair interaction. The functional inequalities come from convexity. We
prove and characterize optimality in the case of quadratic confinement via a
factorization of the measure. This optimality phenomenon holds for all beta
Hermite ensembles including the Gaussian unitary ensemble, a famous exactly
solvable model of random matrix theory. We further explore exact solvability by
reviewing the relation to Dyson-Ornstein-Uhlenbeck diffusion dynamics admitting
the Hermite-Lassalle orthogonal polynomials as a complete set of
eigenfunctions. We also discuss the consequence of the log-Sobolev inequality
in terms of concentration of measure for Lipschitz functions such as maxima and
linear statistics.Comment: Minor improvements. To appear in Geometric Aspects of Functional
Analysis -- Israel Seminar (GAFA) 2017-2019", Lecture Notes in Mathematics
225
Nitrogen fixation in a landrace of maize is supported by a mucilage-associated diazotrophic microbiota
© 2018 Van Deynze et al. http://creativecommons.org/licenses/by/4.0/. Plants are associated with a complex microbiota that contributes to nutrient acquisition, plant growth, and plant defense. Nitrogen-fixing microbial associations are efficient and well characterized in legumes but are limited in cereals, including maize. We studied an indigenous landrace of maize grown in nitrogen-depleted soils in the Sierra Mixe region of Oaxaca, Mexico. This landrace is characterized by the extensive development of aerial roots that secrete a carbohydrate-rich mucilage. Analysis of the mucilage microbiota indicated that it was enriched in taxa for which many known species are diazotrophic, was enriched for homologs of genes encoding nitrogenase subunits, and harbored active nitrogenase activity as assessed by acetylene reduction and 15 N 2 incorporation assays. Field experiments in Sierra Mixe using 15 N natural abundance or 15 N-enrichment assessments over 5 years indicated that atmospheric nitrogen fixation contributed 29%–82% of the nitrogen nutrition of Sierra Mixe maize
- …