Search

Valérie Brunner,1 Bernadette Maynadier,2 Laishun Chen,3 Louise Roques,2 Isabelle Hude,2 Sébastien Séguier,2 Laurence Barthe,1 Philippe Hermann11Pierre Fabre Médicament, Centre de R&D, Toulouse, 2Centre Experimental PreClinque, Campans, France; 3Forest Research Institute Inc., an affiliate of Actavis Inc., Jersey City, NJ, USAAbstract: Levomilnacipran is approved in the US for the treatment of major depressive disorder in adults. We characterized the metabolic profile of levomilnacipran in humans, monkeys, and rats after oral administration of [14C]-levomilnacipran. In vitro binding of levomilnacipran to human plasma proteins was also studied. Unchanged levomilnacipran was the major circulating compound after dosing in all species. Within 12 hours of dosing in humans, levomilnacipran accounted for 52.9% of total plasma radioactivity; the circulating metabolites N-desethyl levomilnacipran N-carbamoyl glucuronide, N-desethyl levomilnacipran, and levomilnacipran N-carbamoyl glucuronide accounted for 11.3%, 7.5%, and 5.6%, respectively. Similar results were seen in monkeys. N-Desethyl levomilnacipran and p-hydroxy levomilnacipran were the main circulating metabolites in rats. Mass balance results indicated that renal excretion was the major route of elimination with 58.4%, 35.5%, and 40.2% of total radioactivity being excreted as unchanged levomilnacipran in humans, monkeys, and rats, respectively. N-Desethyl levomilnacipran was detected in human, monkey, and rat urine (18.2%, 12.4%, and 7.9% of administered dose, respectively). Human and monkey urine contained measurable quantities of levomilnacipran glucuronide (3.8% and 4.1% of administered dose, respectively) and N-desethyl levomilnacipran glucuronide (3.2% and 2.3% of administered dose, respectively); these metabolites were not detected in rat urine. The metabolites p-hydroxy levomilnacipran and p-hydroxy levomilnacipran glucuronide were detected in human urine (≤1.2% of administered dose), and p-hydroxy levomilnacipran glucuronide was found in rat urine (4% of administered dose). None of the metabolites were pharmacologically active. Levomilnacipran was widely distributed with low plasma protein binding (22%).Keywords: FETZIMA, metabolites, mass balance, excretion, tissue distribution

We study parameter estimation for one-dimensional energy balance models with memory (EBMMs) given localized and noisy temperature measurements. Our results apply to a wide range of nonlinear, parabolic partial differential equations with integral memory terms. First, we show that a space-dependent parameter can be determined uniquely everywhere in the PDE's domain of definition D, using only temperature information in a small subdomain E ⊂D. This result is valid only when the data correspond to exact measurements of the temperature. We propose a method for estimating a model parameter of the EBMM using more realistic, error-contaminated temperature data derived, for example, from ice cores or marine-sediment cores. Our approach is based on a so-called mechanistic-statistical model that combines a deterministic EBMM with a statistical model of the observation process. Estimating a parameter in this setting is especially challenging, because the observation process induces a strong loss of information. Aside from the noise contained in past temperature measurements, an additional error is induced by the age-dating method, whose accuracy tends to decrease with a sample's remoteness in time. Using a Bayesian approach, we show that obtaining an accurate parameter estimate is still possible in certain cases. © 2014 The Author(s) Published by the Royal Society. All rights reserved.

The notion of inside dynamics of traveling waves has been introduced in the recent paper [14]. Assuming that a traveling wave u(t,x) = U(x − c t) is made of several components υi ≥ 0 (i ∈ I ⊂ N), the inside dynamics of the wave is then given by the spatio-temporal evolution of the densities of the components υi. For reaction-diffusion equations of the form ∂tu(t,x) = ∂xxu(t,x) + f(u(t,x)), where f is of monostable or bistable type, the results in [14] show that traveling waves can be classified into two main classes: pulled waves and pushed waves. Using the same framework, we study the pulled/pushed nature of the traveling wave solutions of delay equations ∂tu(t,x) = ∂xxu(t,x) + F(u(t −τ,x),u(t,x)) We begin with a review of the latest results on the existence of traveling wave solutions of such equations, for several classical reaction terms. Then, we give analytical and numerical results which describe the inside dynamics of these waves. From a point of view of population ecology, our study shows that the existence of a non-reproductive and motionless juvenile stage can slightly enhance the genetic diversity of a species colonizing an empty environment. [ABSTRACT FROM AUTHOR]

We propose a model to describe the adaptation of a phenotypically structured population in a H-patch environment connected by migration, with each patch associated with a different phenotypic optimum, and we perform a rigorous mathematical analysis of this model. We show that the large-time behaviour of the solution (persistence or extinction) depends on the sign of a principal eigenvalue, [Formula: see text], and we study the dependency of [Formula: see text] with respect to H. This analysis sheds new light on the effect of increasing the number of patches on the persistence of a population, which has implications in agroecology and for understanding zoonoses; in such cases we consider a pathogenic population and the patches correspond to different host species. The occurrence of a springboard effect, where the addition of a patch contributes to persistence, or on the contrary the emergence of a detrimental effect by increasing the number of patches on the persistence, depends in a rather complex way on the respective positions in the phenotypic space of the optimal phenotypes associated with each patch. From a mathematical point of view, an important part of the difficulty in dealing with [Formula: see text], compared to [Formula: see text] or [Formula: see text], comes from the lack of symmetry. Our results, which are based on a fixed point theorem, comparison principles, integral estimates, variational arguments, rearrangement techniques, and numerical simulations, provide a better understanding of these dependencies. In particular, we propose a precise characterisation of the situations where the addition of a third patch increases or decreases the chances of persistence, compared to a situation with only two patches., (© 2023. The Author(s), under exclusive licence to Springer-Verlag GmbH Germany, part of Springer Nature.)

The environmental factors affecting plant reproduction and effective dispersal, in particular biotic interactions, have a strong influence on plant expansion dynamics, but their demographic and genetic consequences remain an understudied body of theory. Here, we use a mathematical model in a one-dimensional space and on a single reproductive period to describe the joint effects of predispersal seed insect predators foraging strategy and plant reproduction strategy (masting) on the spatio-temporal dynamics of seed sources diversity in the colonisation front of expanding plant populations. We show that certain foraging strategies can result in a higher seed predation rate at the colonisation front compared to the core of the population, leading to an Allee effect. This effect promotes the contribution of seed sources from the core to the colonisation front, with long-distance dispersal further increasing this contribution. As a consequence, our study reveals a novel impact of the predispersal seed predation-induced Allee effect, which mitigates the erosion of diversity in expanding populations. We use rearrangement inequalities to show that masting has a buffering role: it mitigates this seed predation-induced Allee effect. This study shows that predispersal seed predation, plant reproductive strategies and seed dispersal patterns can be intermingled drivers of the diversity of seed sources in expanding plant populations, and opens new perspectives concerning the analysis of more complex models such as integro-difference or reaction-diffusion equations., (© 2023. The Author(s), under exclusive licence to Springer-Verlag GmbH Germany, part of Springer Nature.)

A profound transformation of agricultural production methods has become unavoidable due to the increase in the world's population, and environmental and climatic challenges. Agroecology is now recognized as a challenging model for agricultural systems, promoting their diversification and adaptation to environmental and socio-economic contexts, with consequences for the entire agri-food system and the development of rural and urban areas. Through a prospective exercise performed at a large interdisciplinary institute, INRAE, a research agenda for agroecology was built that filled a gap through its ambition and interdisciplinarity. It concerned six topics. For genetics, there is a need to study genetic aspects of complex systems (e.g., mixtures of genotypes) and to develop breeding methods for them. For landscapes, challenges lie in effects of heterogeneity at multiple scales, in multifunctionality and in the design of agroecological landscapes. Agricultural equipment and digital technologies show high potential for monitoring dynamics of agroecosystems. For modeling, challenges include approaches to complexity, consideration of spatial and temporal dimensions and representation of the cascade from cropping practices to ecosystem services. The agroecological transition of farms calls for modeling and observational approaches as well as for creating new design methods. Integration of agroecology into food systems raises the issues of product specificity, consumer behavior and organization of markets, standards and public policies. In addition, transversal priorities were identified: (i) generating sets of biological data, through research and participatory mechanisms, that are appropriate for designing agroecological systems and (ii) collecting and using coherent sets of data to enable assessment of vulnerability, resilience and risk in order to evaluate the performance of agroecological systems and to contribute to scaling up. The main lessons learned from this collective exercise can be useful for the entire scientific community engaged in research into agroecology., Competing Interests: Competing interestsThe authors declare no competing interests., (© The Author(s) 2022.)

Understanding how a diversity of plants in agroecosystems affects the adaptation of pathogens is a key issue in agroecology. We analyze PDE systems describing the dynamics of adaptation of two phenotypically structured populations, under the effects of mutation, selection and migration in a two-patch environment, each patch being associated with a different phenotypic optimum. We consider two types of growth functions that depend on the n-dimensional phenotypic trait: either local and linear or nonlocal nonlinear. In both cases, we obtain existence and uniqueness results as well as a characterization of the large-time behaviour of the solution (persistence or extinction) based on the sign of a principal eigenvalue. We show that migration between the two environments decreases the chances of persistence, with in some cases a 'lethal migration threshold' above which persistence is not possible. Comparison with stochastic individual-based simulations shows that the PDE approach accurately captures this threshold. Our results illustrate the importance of cultivar mixtures for disease prevention and control.

The management of public health and the preparedness for health emergencies partly rely on the collection and analysis of surveillance data, which become crucial in the context of an emergency such as the pandemic caused by COVID-19. For COVID-19, typically, numerous national and global initiatives have been set up from this perspective. Here, we propose to develop a shared vision of the country-level outbreaks during a pandemic, by enhancing, at the international scale, the foundations of the analysis of surveillance data and by adopting a unified and real-time approach to monitor and forecast the outbreak across time and across the world. This proposal, rolled out as a web platform, should contribute to strengthen epidemiological understanding, sanitary democracy as well as global and local responses to pandemics., Competing Interests: Authors declare no competing interests., (© 2020 The Author(s).)

Discrepancies in population structures, decision making, health systems and numerous other factors result in various COVID-19-mortality dynamics at country scale, and make the forecast of deaths in a country under focus challenging. However, mortality dynamics of countries that are ahead of time implicitly include these factors and can be used as real-life competing predicting models. We precisely propose such a data-driven approach implemented in a publicly available web app timely providing mortality curves comparisons and real-time short-term forecasts for about 100 countries. Here, the approach is applied to compare the mortality trajectories of second-line and front-line European countries facing the COVID-19 epidemic wave. Using data up to mid-April, we show that the second-line countries generally followed relatively mild mortality curves rather than fast and severe ones. Thus, the continuation, after mid-April, of the COVID-19 wave across Europe was likely to be mitigated and not as strong as it was in most of the front-line countries first impacted by the wave (this prediction is corroborated by posterior data)., Competing Interests: The authors have declared that no competing interests exist.

The COVID-19 epidemic was reported in the Hubei province in China in December 2019 and then spread around the world reaching the pandemic stage at the beginning of March 2020. Since then, several countries went into lockdown. Using a mechanistic-statistical formalism, we estimate the effect of the lockdown in France on the contact rate and the effective reproduction number R e of the COVID-19. We obtain a reduction by a factor 7 ( R e = 0.47, 95%-CI: 0.45-0.50), compared to the estimates carried out in France at the early stage of the epidemic. We also estimate the fraction of the population that would be infected by the beginning of May, at the official date at which the lockdown should be relaxed. We find a fraction of 3.7% (95%-CI: 3.0-4.8%) of the total French population, without taking into account the number of recovered individuals before April 1st, which is not known. This proportion is seemingly too low to reach herd immunity. Thus, even if the lockdown strongly mitigated the first epidemic wave, keeping a low value of R e is crucial to avoid an uncontrolled second wave (initiated with much more infectious cases than the first wave) and to hence avoid the saturation of hospital facilities., (Copyright © 2020 Roques, Klein, Papaïx, Sar and Soubeyrand.)

The number of screening tests carried out in France and the methodology used to target the patients tested do not allow for a direct computation of the actual number of cases and the infection fatality ratio (IFR). The main objective of this work is to estimate the actual number of people infected with COVID-19 and to deduce the IFR during the observation window in France. We develop a `mechanistic-statistical' approach coupling a SIR epidemiological model describing the unobserved epidemiological dynamics, a probabilistic model describing the data acquisition process and a statistical inference method. The actual number of infected cases in France is probably higher than the observations: we find here a factor ×8 (95%-CI: 5-12) which leads to an IFR in France of 0.5% (95%-CI: 0.3-0.8) based on hospital death counting data. Adjusting for the number of deaths in nursing homes, we obtain an IFR of 0.8% (95%-CI: 0.45-1.25). This IFR is consistent with previous findings in China (0.66%) and in the UK (0.9%) and lower than the value previously computed on the Diamond Princess cruse ship data (1.3%).

Populations may genetically adapt to severe stress that would otherwise cause their extirpation. Recent theoretical work, combining stochastic demography with Fisher's geometric model of adaptation, has shown how evolutionary rescue becomes unlikely beyond some critical intensity of stress. Increasing mutation rates may however allow adaptation to more intense stress, raising concerns about the effectiveness of treatments against pathogens. This previous work assumes that populations are rescued by the rise of a single resistance mutation. However, even in asexual organisms, rescue can also stem from the accumulation of multiple mutations in a single genome. Here, we extend previous work to study the rescue process in an asexual population where the mutation rate is sufficiently high so that such events may be common. We predict both the ultimate extinction probability of the population and the distribution of extinction times. We compare the accuracy of different approximations covering a large range of mutation rates. Moderate increase in mutation rates favors evolutionary rescue. However, larger increase leads to extinction by the accumulation of a large mutation load, a process called lethal mutagenesis. We discuss how these results could help design "evolution-proof" antipathogen treatments that even highly mutable strains could not overcome., (© 2019 The Author(s). Evolution © 2019 The Society for the Study of Evolution.)

Using a free boundary approach based on an analogy with ice melting models, we propose a deterministic PDE framework to describe the dynamics of fitness distributions in the presence of beneficial mutations with non-epistatic effects on fitness. Contrarily to most approaches based on deterministic models, our framework does not rely on an infinite population size assumption, and successfully captures the transient as well as the long time dynamics of fitness distributions. In particular, consistently with stochastic individual-based approaches or stochastic PDE approaches, it leads to a constant asymptotic rate of adaptation at large times, that most deterministic approaches failed to describe. We derive analytic formulas for the asymptotic rate of adaptation and the full asymptotic distribution of fitness. These formulas depend explicitly on the population size, and are shown to be accurate for a wide range of population sizes and mutation rates, compared to individual-based simulations. Although we were not able to derive an analytic description for the transient dynamics, numerical computations lead to accurate predictions and are computationally efficient compared to stochastic simulations. These computations show that the fitness distribution converges towards a travelling wave with constant speed, and whose profile can be computed analytically.

Various models describe asexual evolution by mutation, selection, and drift. Some focus directly on fitness, typically modeling drift but ignoring or simplifying both epistasis and the distribution of mutation effects (traveling wave models). Others follow the dynamics of quantitative traits determining fitness (Fisher's geometric model), imposing a complex but fixed form of mutation effects and epistasis, and often ignoring drift. In all cases, predictions are typically obtained in high or low mutation rate limits and for long-term stationary regimes, thus losing information on transient behaviors and the effect of initial conditions. Here, we connect fitness-based and trait-based models into a single framework, and seek explicit solutions even away from stationarity. The expected fitness distribution is followed over time via its cumulant generating function, using a deterministic approximation that neglects drift. In several cases, explicit trajectories for the full fitness distribution are obtained for arbitrary mutation rates and standing variance. For nonepistatic mutations, especially with beneficial mutations, this approximation fails over the long term but captures the early dynamics, thus complementing stationary stochastic predictions. The approximation also handles several diminishing returns epistasis models (e.g., with an optimal genotype); it can be applied at and away from equilibrium. General results arise at equilibrium, where fitness distributions display a "phase transition" with mutation rate. Beyond this phase transition, in Fisher's geometric model, the full trajectory of fitness and trait distributions takes a simple form; robust to the details of the mutant phenotype distribution. Analytical arguments are explored regarding why and when the deterministic approximation applies., (Copyright © 2016 by the Genetics Society of America.)