On the birth of limit cycles for non-smooth dynamical systems

The main objective of this work is to develop, via Brower degree theory and regularization theory, a variation of the classical averaging method for detecting limit cycles of certain piecewise continuous dynamical systems. In fact, overall results are presented to ensure the existence of limit cycles of such systems. These results may represent new insights in averaging, in particular its relation with non smooth dynamical systems theory. An application is presented in careful detail

Transcriptome‐based phylogeny of endemic Lake Baikal amphipod species flock: fast speciation accompanied by frequent episodes of positive selection

Endemic species flocks inhabiting ancient lakes, oceanic islands and other long‐lived isolated habitats are often interpreted as adaptive radiations. Yet molecular evidence for directional selection during species flocks radiation is scarce. Using partial transcriptomes of 64 species of Lake Baikal (Siberia, Russia) endemic amphipods and two nonendemic outgroups, we report a revised phylogeny of this species flock and analyse evidence for positive selection within the endemic lineages. We confirm two independent invasions of amphipods into Baikal and demonstrate that several morphological features of Baikal amphipods, such as body armour and reduction in appendages and sensory organs, evolved in several lineages in parallel. Radiation of Baikal amphipods has been characterized by short phylogenetic branches and frequent episodes of positive selection which tended to be more frequent in the early phase of the second invasion of amphipods into Baikal when the most intensive diversification occurred. Notably, signatures of positive selection are frequent in genes encoding mitochondrial membrane proteins with electron transfer chain and ATP synthesis functionality. In particular, subunits of both the membrane and substrate‐level ATP synthases show evidence of positive selection in the plankton species Macrohectopus branickii, possibly indicating adaptation to active plankton lifestyle and to survival under conditions of low temperature and high hydrostatic pressures known to affect membranes functioning. Other functional categories represented among genes likely to be under positive selection include Ca‐binding muscle‐related proteins, possibly indicating adaptation to Ca‐deficient low mineralization Baikal waters.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/136009/1/mec13927.pdfhttp://deepblue.lib.umich.edu/bitstream/2027.42/136009/2/mec13927_am.pd

Recombination rate and selection strength in HIV intra-patient evolution

The evolutionary dynamics of HIV during the chronic phase of infection is driven by the host immune response and by selective pressures exerted through drug treatment. To understand and model the evolution of HIV quantitatively, the parameters governing genetic diversification and the strength of selection need to be known. While mutation rates can be measured in single replication cycles, the relevant effective recombination rate depends on the probability of coinfection of a cell with more than one virus and can only be inferred from population data. However, most population genetic estimators for recombination rates assume absence of selection and are hence of limited applicability to HIV, since positive and purifying selection are important in HIV evolution. Here, we estimate the rate of recombination and the distribution of selection coefficients from time-resolved sequence data tracking the evolution of HIV within single patients. By examining temporal changes in the genetic composition of the population, we estimate the effective recombination to be r=1.4e-5 recombinations per site and generation. Furthermore, we provide evidence that selection coefficients of at least 15% of the observed non-synonymous polymorphisms exceed 0.8% per generation. These results provide a basis for a more detailed understanding of the evolution of HIV. A particularly interesting case is evolution in response to drug treatment, where recombination can facilitate the rapid acquisition of multiple resistance mutations. With the methods developed here, more precise and more detailed studies will be possible, as soon as data with higher time resolution and greater sample sizes is available.Comment: to appear in PLoS Computational Biolog

A mathematical framework for critical transitions: normal forms, variance and applications

Critical transitions occur in a wide variety of applications including mathematical biology, climate change, human physiology and economics. Therefore it is highly desirable to find early-warning signs. We show that it is possible to classify critical transitions by using bifurcation theory and normal forms in the singular limit. Based on this elementary classification, we analyze stochastic fluctuations and calculate scaling laws of the variance of stochastic sample paths near critical transitions for fast subsystem bifurcations up to codimension two. The theory is applied to several models: the Stommel-Cessi box model for the thermohaline circulation from geoscience, an epidemic-spreading model on an adaptive network, an activator-inhibitor switch from systems biology, a predator-prey system from ecology and to the Euler buckling problem from classical mechanics. For the Stommel-Cessi model we compare different detrending techniques to calculate early-warning signs. In the epidemics model we show that link densities could be better variables for prediction than population densities. The activator-inhibitor switch demonstrates effects in three time-scale systems and points out that excitable cells and molecular units have information for subthreshold prediction. In the predator-prey model explosive population growth near a codimension two bifurcation is investigated and we show that early-warnings from normal forms can be misleading in this context. In the biomechanical model we demonstrate that early-warning signs for buckling depend crucially on the control strategy near the instability which illustrates the effect of multiplicative noise.Comment: minor corrections to previous versio

Correlated Evolution of Nearby Residues in Drosophilid Proteins

Here we investigate the correlations between coding sequence substitutions as a function of their separation along the protein sequence. We consider both substitutions between the reference genomes of several Drosophilids as well as polymorphisms in a population sample of Zimbabwean Drosophila melanogaster. We find that amino acid substitutions are “clustered” along the protein sequence, that is, the frequency of additional substitutions is strongly enhanced within ≈10 residues of a first such substitution. No such clustering is observed for synonymous substitutions, supporting a “correlation length” associated with selection on proteins as the causative mechanism. Clustering is stronger between substitutions that arose in the same lineage than it is between substitutions that arose in different lineages. We consider several possible origins of clustering, concluding that epistasis (interactions between amino acids within a protein that affect function) and positional heterogeneity in the strength of purifying selection are primarily responsible. The role of epistasis is directly supported by the tendency of nearby substitutions that arose on the same lineage to preserve the total charge of the residues within the correlation length and by the preferential cosegregation of neighboring derived alleles in our population sample. We interpret the observed length scale of clustering as a statistical reflection of the functional locality (or modularity) of proteins: amino acids that are near each other on the protein backbone are more likely to contribute to, and collaborate toward, a common subfunction

Genetic Diversity in the SIR Model of Pathogen Evolution

We introduce a model for assessing the levels and patterns of genetic diversity in pathogen populations, whose epidemiology follows a susceptible-infected-recovered model (SIR). We model the population of pathogens as a metapopulation composed of subpopulations (infected hosts), where pathogens replicate and mutate. Hosts transmit pathogens to uninfected hosts. We show that the level of pathogen variation is well predicted by analytical expressions, such that pathogen neutral molecular variation is bounded by the level of infection and increases with the duration of infection. We then introduce selection in the model and study the invasion probability of a new pathogenic strain whose fitness (R0(1+s)) is higher than the fitness of the resident strain (R0). We show that this invasion probability is given by the relative increment in R0 of the new pathogen (s). By analyzing the patterns of genetic diversity in this framework, we identify the molecular signatures during the replacement and compare these with those observed in sequences of influenza A

The stability of ecosystems: a brief overview of the paradox of enrichment

In theory, enrichment of resource in a predator-prey model leads to destabilization of the system, thereby collapsing the trophic interaction, a phenomenon referred to as "the paradox of enrichment". After it was first proposed by Rosenzweig (1971), a number of subsequent studies were carried out on this dilemma over many decades. In this article, we review these theoretical and experimental works and give a brief overview of the proposed solutions to the paradox. The mechanisms that have been discussed are modifications of simple predator-prey models in the presence of prey that is inedible, invulnerable, unpalatable and toxic. Another class of mechanisms includes an incorporation of a ratio-dependent functional form, inducible defence of prey and density-dependent mortality of the predator. Moreover, we find a third set of explanations based on complex population dynamics including chaos in space and time. We conclude that, although any one of the various mechanisms proposed so far might potentially prevent destabilization of the predator-prey dynamics following enrichment, in nature different mechanisms may combine to cause stability, even when a system is enriched. The exact mechanisms, which may differ among systems, need to be disentangled through extensive field studies and laboratory experiments coupled with realistic theoretical models

The Inexorable Spread of a Newly Arisen Neo-Y Chromosome

A newly arisen Y-chromosome can become established in one part of a species range by genetic drift or through the effects of selection on sexually antagonistic alleles. However, it is difficult to explain why it should then spread throughout the species range after this initial episode. As it spreads into new populations, it will actually enter females. It would then be expected to perform poorly since it will have been shaped by the selective regime of the male-only environment from which it came. We address this problem using computer models of hybrid zone dynamics where a neo-XY chromosomal race meets the ancestral karyotype. Our models consider that the neo-Y was established by the fusion of an autosome with the ancestral X-chromosome (thereby creating the Y and the ‘fused X’). Our principal finding is that sexually antagonistic effects of the Y induce indirect selection in favour of the fused X-chromosomes, causing their spread. The Y-chromosome can then spread, protected behind the advancing shield of the fused X distribution. This mode of spread provides a robust explanation of how newly arisen Y-chromosomes can spread. A Y-chromosome would be expected to accumulate mutations that would cause it to be selected against when it is a rare newly arrived migrant. The Y can spread, nevertheless, because of the indirect selection induced by gene flow (which can only be observed in models comprising multiple populations). These results suggest a fundamental re-evaluation of sex-chromosome hybrid zones. The well-understood evolutionary events that initiate the Y-chromosome's degeneration will actually fuel its range expansion

Fitness Consequences of Advanced Ancestral Age over Three Generations in Humans

A rapid rise in age at parenthood in contemporary societies has increased interest in reports of higher prevalence of de novo mutations and health problems in individuals with older fathers, but the fitness consequences of such age effects over several generations remain untested. Here, we use extensive pedigree data on seven pre-industrial Finnish populations to show how the ages of ancestors for up to three generations are associated with fitness traits. Individuals whose fathers, grandfathers and great-grandfathers fathered their lineage on average under age 30 were ~13% more likely to survive to adulthood than those whose ancestors fathered their lineage at over 40 years. In addition, females had a lower probability of marriage if their male ancestors were older. These findings are consistent with an increase of the number of accumulated de novo mutations with male age, suggesting that deleterious mutations acquired from recent ancestors may be a substantial burden to fitness in humans. However, possible non-mutational explanations for the observed associations are also discussed