No Arabic abstract
This paper develops simplified mathematical models describing the mutation-selection balance for the asexual and sexual replication pathways in {it Saccharomyces cerevisiae}. We assume diploid genomes consisting of two chromosomes, and we assume that each chromosome is functional if and only if its base sequence is identical to some master sequence. The growth and replication of the yeast cells is modeled as a first-order process, with first-order growth rate constants that are determined by whether a given genome consists of zero, one, or two functional chromosomes. In the asexual pathway, we assume that a given diploid cell divides into two diploids. In the sexual pathway, we assume that a given diploid cell divides into two diploids, each of which then divide into two haploids. The resulting four haploids enter a haploid pool, where they grow and replicate until they meet another haploid with which to fuse. When the cost for sex is low, we find that the selective mating strategy leads to the highest mean fitness of the population, when compared to all of the other strategies. We also show that, at low to intermediate replication fidelities, sexual replication with random mating has a higher mean fitness than asexual replication, as long as the cost for sex is low. This is consistent with previous work suggesting that sexual replication is advantageous at high population densities, low replication rates, and intermediate replication fidelities. The results of this paper also suggest that {it S. cerevisiae} switches from asexual to sexual replication when stressed, because stressful growth conditions provide an opportunity for the yeast to clear out deleterious mutations from their genomes.
This paper studies the mutation-selection balance in three simplified replication models. The first model considers a population of organisms replicating via the production of asexual spores. The second model considers a sexually replicating population that produces identical gametes. The third model considers a sexually replicating population that produces distinct sperm and egg gametes. All models assume diploid organisms whose genomes consist of two chromosomes, each of which is taken to be functional if equal to some master sequence, and defective otherwise. In the asexual population, the asexual diploid spores develop directly into adult organisms. In the sexual populations, the haploid gametes enter a haploid pool, where they may fuse with other haploids. The resulting immature diploid organisms then proceed to develop into mature organisms. Based on an analysis of all three models, we find that, as organism size increases, a sexually replicating population can only outcompete an asexually replicating population if the adult organisms produce distinct sperm and egg gametes. A sexual replication strategy that is based on the production of large numbers of sperm cells to fertilize a small number of eggs is found to be necessary in order to maintain a sufficiently low cost for sex for the strategy to be selected for over a purely asexual strategy. We discuss the usefulness of this model in understanding the evolution and maintenance of sexual replication as the preferred replication strategy in complex, multicellular organisms.
Woolly mammoths (Mammuthus primigenius) populated Siberia, Beringia, and North America during the Pleistocene and early Holocene. Recent breakthroughs in ancient DNA sequencing have allowed for complete genome sequencing for two specimens of woolly mammoths (Palkopoulou et al. 2015). One mammoth specimen is from a mainland population ~45,000 years ago when mammoths were plentiful. The second, a 4300 yr old specimen, is derived from an isolated population on Wrangel island where mammoths subsisted with small effective population size more than 43-fold lower than previous populations. These extreme differences in effective population size offer a rare opportunity to test nearly neutral models of genome architecture evolution within a single species. Using these previously published mammoth sequences, we identify deletions, retrogenes, and non-functionalizing point mutations. In the Wrangel island mammoth, we identify a greater number of deletions, a larger proportion of deletions affecting gene sequences, a greater number of candidate retrogenes, and an increased number of premature stop codons. This accumulation of detrimental mutations is consistent with genomic meltdown in response to low effective population sizes in the dwindling mammoth population on Wrangel island. In addition, we observe high rates of loss of olfactory receptors and urinary proteins, either because these loci are non-essential or because they were favored by divergent selective pressures in island environments. Finally, at the locus of FOXQ1 we observe two independent loss-of-function mutations, which would confer a satin coat phenotype in this island woolly mammoth.
This paper develops a quasispecies model that incorporates the SOS response. We consider a unicellular, asexually replicating population of organisms, whose genomes consist of a single, double-stranded DNA molecule, i.e. one chromosome. We assume that repair of post-replication mismatched base-pairs occurs with probability $ lambda $, and that the SOS response is triggered when the total number of mismatched base-pairs exceeds $ l_S $. We further assume that the per-mismatch SOS elimination rate is characterized by a first-order rate constant $ kappa_{SOS} $. For a single fitness peak landscape where the master genome can sustain up to $ l $ mismatches and remain viable, this model is analytically solvable in the limit of infinite sequence length. The results, which are confirmed by stochastic simulations, indicate that the SOS response does indeed confer a fitness advantage to a population, provided that it is only activated when DNA damage is so extensive that a cell will die if it does not attempt to repair its DNA.
Exploring the genetic basis of heritable traits remains one of the central challenges in biomedical research. In simple cases, single polymorphic loci explain a significant fraction of the phenotype variability. However, many traits of interest appear to be subject to multifactorial control by groups of genetic loci instead. Accurate detection of such multivariate associations is nontrivial and often hindered by limited power. At the same time, confounding influences such as population structure cause spurious association signals that result in false positive findings if they are not accounted for in the model. Here, we propose LMM-Lasso, a mixed model that allows for both, multi-locus mapping and correction for confounding effects. Our approach is simple and free of tuning parameters, effectively controls for population structure and scales to genome-wide datasets. We show practical use in genome-wide association studies and linkage mapping through retrospective analyses. In data from Arabidopsis thaliana and mouse, our method is able to find a genetic cause for significantly greater fractions of phenotype variation in 91% of the phenotypes considered. At the same time, our model dissects this variability into components that result from individual SNP effects and population structure. In addition to this increase of genetic heritability, enrichment of known candidate genes suggests that the associations retrieved by LMM-Lasso are more likely to be genuine.
Increasing number in global COVID-19 cases demands for mathematical model to analyze the interaction between the virus dynamics and the response of innate and adaptive immunity. Here, based on the assumption of a weak and delayed response of the innate and adaptive immunity in SARS-CoV-2 infection, we constructed a mathematical model to describe the dynamic processes of immune system. Integrating theoretical results with clinical COVID-19 patients data, we classified the COVID-19 development processes into three typical modes of immune responses, correlated with the clinical classification of mild & moderate, severe and critical patients. We found that the immune efficacy (the ability of host to clear virus and kill infected cells) and the lymphocyte supply (the abundance and pool of naive T and B cell) play important roles in the dynamic process and determine the clinical outcome, especially for the severe and critical patients. Furthermore, we put forward possible treatment strategies for the three typical modes of immune response. We hope our results can help to understand the dynamical mechanism of the immune response against SARS-CoV-2 infection, and to be useful for the treatment strategies and vaccine design.