Selection in Complex Genetic Systems I. the Symmetric Equilibria of the Three-Locus Symmetric Viability Model

The symmetric equilibria of the three-locus symmetric viability model are determined and their stability analyzed. For tight linkage there may be four stable equilibria, each characterized by having one pair of complementary chromosomes in high frequencies, with all others low. For looser linkage the only stable symmetric equilibrium is that with complete linkage equilibrium. For intermediate recombination values both types of equilibria may be stable. A new class of equilibria with all pairwise linkage disequilibria zero, but with third order linkage disequilibrium, has been discovered. It may be stable for tight linkage.     

Effects of releasing maladapted individuals: a demographic-evolutionary model

A model of the joint dynamics of change in population size N and evolution in a quantitative trait z, as a result of a general form of density dependence, local stabilizing selection, and immigration of individuals deviating from the local optimum, is analyzed. For weak selection and migration, a reduction in total equilibrium population size below the initial level without immigration, K, is shown to occur if the immigrants deviates more than square root of 8 = 2.83 genetic standard deviations from the optimum and if the rate of migration m is sufficiently large relative to the strength of stabilizing selection s. For the Lotka-Volterra form of density dependence, two additional equilibria are shown to exist below K, provided that the strength of selection is large relative to the strength of density dependence. Reintroduction of an initially extinct population is possible if the immigrants are not too maladapted and if the genetic variance is sufficiently large. For a simplified version of the model corresponding to competition between similar species or different haplotypes, the equilibrium population size is always exactly at K if m < Ksz1(2) and is above K otherwise, which shows the importance of including recombination in the model.     

The multiple-locus symmetric fertility model

Selection due to variation in the fecundity among matings of genotypes with respect to many loci each with two alleles is studied. The fitness of a mating depends only on the genotypic distinction between homozygote and heterozygote at each locus in the two individuals, and differences among loci are allowed. This symmetric fertility model is therefore a generalization of the multiple-locus symmetric viability model. The phenomena seen in the two-locus symmetric fertility model generalize—e.g., the possibility of joint stability of equilibria with linkage equilibrium and with linkage disequilibrium, and the existence of different types of totally polymorphic equilibria with the gametic proportions in linkage equilibrium. The central equilibrium with genotypic frequencies in Hardy-Weinberg proportions and gametic frequencies in Robbins proportions exists for all symmetric fertility models. For some symmetric fertility regimes additional equilibria exist with gametic frequencies in linkage equilibrium and with genotypic frequencies in Hardy-Weinberg proportions at all except one locus. These equilibria may exist in the dioecious symmetric viability model, and then they will be locally stable. For free recombination the stable equilibria show linkage equilibrium, but several of these with different numbers of polymorphic loci may be stable simultaneously.     

Modeling the impact of immigration on the epidemiology of tuberculosis

This paper presents two new theoretical frameworks to investigate the impact of immigration on the transmission dynamics of tuberculosis. For the basic model, we present new analysis on the existence and stability of equilibria. Then, we use numerical simulations of the model to illustrate the behavior of the system. We apply the model to Canadian reported data on tuberculosis and observe a good agreement between the model prediction and the data. For the extended model, which incorporated the recruitment of the latent and infectious in immigrants to the basic model, we find that the usual threshold condition does not apply and a unique equilibrium exists for all parameter values. This indicates that the disease does not disappear and becomes endemic in host areas. This finding is also supported by numerical simulations with the extended model. Our study suggests that immigrants have a considerable influence on the overall transmission dynamics behavior of tuberculosis.     

The two-locus model of Gaussian stabilizing selection

We study the equilibrium structure of a well-known two-locus model in which two diallelic loci contribute additively to a quantitative trait that is under Gaussian stabilizing selection. The population is assumed to be infinitely large, randomly mating, and having discrete generations. The two loci may have arbitrary effects on the trait, the strength of selection and the recombination rate may also be arbitrary. We find that 16 different equilibrium patterns exist, having up to 11 equilibria; up to seven interior equilibria may coexist, and up to four interior equilibria, three in negative and one in positive linkage disequilibrium, may be simultaneously stable. Also, two monomorphic and two fully polymorphic equilibria may be simultaneously stable. Therefore, the result of evolution may be highly sensitive to perturbations in the initial conditions or in the underlying genetic parameters. For the special case of equal effects, global stability results are proved. In the general case, we rely in part on numerical computations. The results are compared with previous analyses of the special case of extremely strong selection, of an approximate model that assumes linkage equilibrium, and of the much simpler quadratic optimum model.     

Immigration can destabilize tri-trophic interactions: implications for conservation of top predators

Top predators often have large home ranges and thus are especially vulnerable to habitat loss and fragmentation. Increasing connectance among habitat patches is therefore a common conservation strategy, based in part on models showing that increased migration between subpopulations can reduce vulnerability arising from population isolation. Although three-dimensional models are appropriate for exploring consequences to top predators, the effects of immigration on tri-trophic interactions have rarely been considered. To explore the effects of immigration on the equilibrium abundances of top predators, we studied the effects of immigration in the three-dimensional Rosenzweig-MacArthur model. To investigate the stability of the top predator equilibrium, we used MATCONT to perform a bifurcation analysis. For some combinations of model parameters with low rates of top predator immigration, population trajectories spiral towards a stable focus. Holding other parameters constant, as immigration rate is increased, a supercritical Hopf bifurcation results in a stable limit cycle and thus top predator populations that cycle between high and low abundances. Furthermore, bistability arises as immigration of the intermediate predator is increased. In this case, top predators may exist at relatively low abundances while prey become extinct, or for other initial conditions, the relatively higher top predator abundance controls intermediate predators allowing for non-zero prey population abundance and increased diversity. Thus, our results reveal one of two outcomes when immigration is added to the model. First, over some range of top predator immigration rates, population abundance cycles between high and low values, making extinction from the trough of such cycles more likely than otherwise. Second, for relatively higher intermediate predator migration rates, top predators may exist at low values in a truncated system with impoverished diversity, again with extinction more likely.     

Competitive lipase-catalyzed ester hydrolysis and ammoniolysis in organic solvents; equilibrium model of a solid-liquid-vapor system.

Enzymatic ester hydrolysis and ammoniolysis were performed as competitive reactions in methyl isobutyl ketone without a separate aqueous phase. The reaction system contained solid ammonium bicarbonate, which dissolved as water, ammonia, and carbon dioxide. During the reaction an organic liquid phase, a vapor phase, and at least one solid phase are present. The overall equilibrium composition of this multiphase system is a complex function of the reaction equilibria and several phase equilibria. To gain a quantitative understanding of this system a mathematical model was developed and evaluated. The model is based on the mass balances for a closed batch system and straightforward relations for the reaction equilibria and the solubility equilibria of ammonium bicarbonate, the fatty acid ammonium salt, water, ammonia, and carbon dioxide. For butyl butyrate as a model ester and Candida antarctica lipase B as the biocatalyst this equilibrium model describes the experiments satisfactorily. The model predicts that high equilibrium yields of butyric acid can be achieved only in the absence of ammoniolysis or in the presence of a separate water phase. However, high yields of butyramide should be possible if the water concentration is fixed at a low level and a more suited source of ammonia is applied.     

A system dynamics model of island biogeography

The MacArthur-Wilson equilibrium theory of island biogeography has been one of the more influential concepts in modern biogeography and ecology. In this paper, we synthesize the theory and examine effects of different immigration/extinction rate-species diversity curves on original predictions from the theory by using the System Dynamics simulation modeling approach. Moreover, we develop a comprehensive and generic System Dynamics model to incorporate a variety of recent modifications and extensions of the theory, including area effect, distance effect, competition effect, habitat diversity effect, target effect, and rescue effect. Through computer simulation with STELLA, a more profound understanding of the theory of island biogeography can be gained. The System Dynamics modeling approach is especially appropriate for such a study because it maximizes the utilization of the ecological data by incorporating qualitative information so that a complex, imprecisely-defined ecological system can be studied quantitatively, effectively, and comprehensively. Our simulation results show that different monotonic rate-species diversity curves do not affect the essence of the theory of island biogeography, while the magnitude of equilibrium species diversity may be greatly affected. Non-monotonic rate-species diversity curves may result in potential multiple equilibria of species diversity. In addition, our model suggests that a non-monotonic relationship may exist between the equilibrium turnover rate and island area and between the equilibrium turnover rate and distance.     

Non-smooth plant disease models with economic thresholds

In order to control plant diseases and eventually maintain the number of infected plants below an economic threshold, a specific management strategy called the threshold policy is proposed, resulting in Filippov systems. These are a class of piecewise smooth systems of differential equations with a discontinuous right-hand side. The aim of this work is to investigate the global dynamic behavior including sliding dynamics of one Filippov plant disease model with cultural control strategy. We examine a Lotka–Volterra Filippov plant disease model with proportional planting rate, which is globally studied in terms of five types of equilibria. For one type of equilibrium, the global structure is discussed by the iterative equations for initial numbers of plants. For the other four types of equilibria, the bounded global attractor of each type is obtained by constructing appropriate Lyapunov functions. The ideas of constructing Lyapunov functions for Filippov systems, the methods of analyzing such systems and the main results presented here provide scientific support for completing control regimens on plant diseases in integrated disease management.     

The Worst-Case Weighted Multi-Objective Game with an Application to Supply Chain Competitions

In this paper, we propose a worst-case weighted approach to the multi-objective n-person non-zero sum game model where each player has more than one competing objective. Our “worst-case weighted multi-objective game” model supposes that each player has a set of weights to its objectives and wishes to minimize its maximum weighted sum objectives where the maximization is with respect to the set of weights. This new model gives rise to a new Pareto Nash equilibrium concept, which we call “robust-weighted Nash equilibrium”. We prove that the robust-weighted Nash equilibria are guaranteed to exist even when the weight sets are unbounded. For the worst-case weighted multi-objective game with the weight sets of players all given as polytope, we show that a robust-weighted Nash equilibrium can be obtained by solving a mathematical program with equilibrium constraints (MPEC). For an application, we illustrate the usefulness of the worst-case weighted multi-objective game to a supply chain risk management problem under demand uncertainty. By the comparison with the existed weighted approach, we show that our method is more robust and can be more efficiently used for the real-world applications.     

The role of screening and treatment in the transmission dynamics of HIV/AIDS and tuberculosis co-infection: a mathematical study

In this paper, we present a deterministic non-linear mathematical model for the transmission dynamics of HIV and TB co-infection and analyze it in the presence of screening and treatment. The equilibria of the model are computed and stability of these equilibria is discussed. The basic reproduction numbers corresponding to both HIV and TB are found and we show that the disease-free equilibrium is stable only when the basic reproduction numbers for both the diseases are less than one. When both the reproduction numbers are greater than one, the co-infection equilibrium point may exist. The co-infection equilibrium is found to be locally stable whenever it exists. The TB-only and HIV-only equilibria are locally asymptotically stable under some restriction on parameters. We present numerical simulation results to support the analytical findings. We observe that screening with proper counseling of HIV infectives results in a significant reduction of the number of individuals progressing to HIV. Additionally, the screening of TB reduces the infection prevalence of TB disease. The results reported in this paper clearly indicate that proper screening and counseling can check the spread of HIV and TB diseases and effective control strategies can be formulated around ‘screening with proper counseling’.     

Single-species metapopulation dynamics: concepts, models and observations

This paper outlines a conceptual and theoretical framework for single-species metapopulation dynamics based on the Levins model and its variants. The significance of the following factors to metapopulation dynamics are explored: evolutionary changes in colonization ability; habitat patch size and isolation; compensatory effects between colonization and extinction rates; the effect of immigration on local dynamics (the rescue effect); and heterogeneity among habitat patches. The rescue effect may lead to alternative stable equilibria in metapopulation dynamics. Heterogeneity among habitat patches may give rise to a bimodal equilibrium distribution of the fraction of patches occupied in an assemblage of species (the core-satellite distribution). A new model of incidence functions is described, which allows one to estimate species' colonization and extinction rates on islands colonized from mainland. Four distinct kinds of stochasticity affecting metapopulation dynamics are discussed with examples. The concluding section describes four possible scenarios of metapopulation extinction.     

若干具有非线性传染力的传染病模型的稳定性分析

讨论了具有常数迁入和非线性传染力的三类传染病模型,即SIRI模型,SIRI框架下的DS模型及SIR框架下的DI模型。给出了它们基本再生数R0的表达式,证明了R0≤1时无病平衡点是全局稳定的,同时证明了如果地方病平衡点存在,则必是全局稳定的结果（从而必唯一）对第一和第三个模型还给出了R0＞1时地方病平衡点的存在唯一性。     

Equilibria and dynamics of selection at a diallelic autosomal locus in the Nagylaki-Crow continuous model of a monoecious population

Let birth rates and death rates be constant, birth rates positive, fertilities additive, and each birth rate not larger than twice any other birth rate. Global convergence to equilibria is proved for the model in the title. There is at most one polymorphic equilibrium or there are a continuum of equilibria. The phase portraits are given. If there is a polymorphic equilibrium, then the largest negatively invariant set in the state space is a continuous curve connecting the two fixation equilibria. This curve coincides with the Hardy-Weinberg manifold exactly when the death rate is additive. Disregarding extinction, the polymorphic equilibria are the same for the continuous model as for the corresponding discrete model exactly when the death rate is additive.     

On some models of fertility selection

Additive, multiplicative and symmetric models of fertility controlled by one diallelic gene are studied. For the completely symmetric fertility system a complete equilibrium and local stability analysis is possible. Contrary to previous conjectures, asymmetric equilibria can be stable. Conditions are derived under which a multiplicative model can be regarded as equivalent to a symmetric fertility system.     

Linkage and the Maintenance of Heritable Variation by Mutation-Selection Balance

The role of linkage in influencing heritable variation maintained through a balance between mutation and stabilizing selection is investigated for two different models. In both cases one trait is considered and the interactions within and between loci are assumed to be additive. Contrary to most earlier investigations of this problem no a priori assumptions on the distribution of genotypic values are imposed. For a deterministic two-locus two-allele model with recombination and mutation, related to the symmetric viability model, a complete nonlinear analysis is performed. It is shown that, depending on the recombination rate, multiple stable equilibria may coexist. The equilibrium genetic and genic variances are calculated. For a polygenic trait in a finite population with a possible continuum of allelic effects a simulation study is performed. In both models the equilibrium genetic and genic variances are roughly equal to the house-of-cards prediction or its finite population counterpart as long as the recombination rate is not extremely low. However, negative linkage disequilibrium builds up. If the loci are very closely linked the equilibrium additive genetic variance is slightly lower than the house-of-cards prediction, but the genic variance is much higher. Depending on whether the parameters are in favor of the house-of-cards or the Gaussian approximation, different behavior of the genetic system occurs with respect to linkage.     

Structured models of metapopulation dynamics

I develop models of metapopulation dynamics that describe changes in the numbers of individuals within patches. These models are analogous to structured population models, with patches playing the role of individuals. Single species models which do not include the effect of immigration on local population dynamics of occupied patches typically lead to a unique equilibrium. The models can be used to study the distributions of numbers of individuals among patches, showing that both metapopulations with local outbreaks and metapopulations without outbreaks can occur in systems with no underlying environmental variability. Distributions of local population sizes (in occupied patches) can vary independently of the total population size, so both patterns of distributions of local population sizes are compatible with either rare or common species. Models which include the effect of immigration on local population dynamics can lead to two positive equilibria, one stable and one unstable, the latter representing a threshold between regional extinction and persistence.     

一个具暂时免疫且总人数可变的传染病动力学模型

建立了一个具常恢复率和接触率依赖于总人数的SIRS传染病动力学模型，讨论了系统平衡点的存在性和稳定性，对双线性传染率的特殊情形，给出了传染病平衡点的全局稳定性结论，推广和改进了已有的相应结果。     

A mathematical model for the spread of Strepotococcus pneumoniae with transmission dependent on serotype

We examine a mathematical model for the transmission of Streptococcus Pneumoniae amongst young children when the carriage transmission coefficient depends on the serotype. Carriage means pneumococcal colonization. There are two sequence types (STs) spreading in a population each of which can be expressed as one of two serotypes. We derive the differential equation model for the carriage spread and perform an equilibrium and global stability analysis on it. A key parameter is the effective reproduction number R (e). For R (e) ≤ 1, there is only the carriage-free equilibrium (CFE) and the carriage will die out whatever be the starting values. For R (e) > 1, unless the effective reproduction numbers of the two STs are equal, in addition to the CFE there are two carriage equilibria, one for each ST. If the ST with the largest effective reproduction number is initially present, then in the long-term the carriage will tend to the corresponding equilibrium.     

Conditions for the existence of stationary densities for some two-dimensional diffusion processes with applications in population biology

Conditions are derived that we conjecture are necessary and sufficient for the existence of stationary densities for a class of two-dimensional diffusion processes. The derivation of the conditions rests on the assumption that a two-dimensional stationary density (which can be viewed as a stable “internal equilibrium”) exists if and only if all “boundary equilibria” are unstable in the sense that small perturbations lead to moving away from the boundaries with high probability. For the models considered, the boundary equilibria are one-dimensional stationary densities and equilibrium points. To demonstrate the usefulness of the conditions, three random environment models are analyzed: a three-allele selection model, a two-species competition model, and a two-locus selection model. Several of the results obtained have been verified by alternate methods.