Analytic approach to the evolutionary effects of genetic exchange

We present an approximate analytic study of our previously introduced model of evolution including the effects of genetic exchange. This model is motivated by the process of bacterial transformation. We solve for the velocity, the rate of increase of fitness, as a function of the fixed population size, $N$. We find the velocity increases with $\ln N$, eventually saturated at an $N$ which depends on the strength of the recombination process. The analytical treatment is seen to agree well with direct numerical simulations of our model equations

An Evolutionary Reduction Principle for Mutation Rates at Multiple Loci

A model of mutation rate evolution for multiple loci under arbitrary selection is analyzed. Results are obtained using techniques from Karlin (1982) that overcome the weak selection constraints needed for tractability in prior studies of multilocus event models. A multivariate form of the reduction principle is found: reduction results at individual loci combine topologically to produce a surface of mutation rate alterations that are neutral for a new modifier allele. New mutation rates survive if and only if they fall below this surface - a generalization of the hyperplane found by Zhivotovsky et al. (1994) for a multilocus recombination modifier. Increases in mutation rates at some loci may evolve if compensated for by decreases at other loci. The strength of selection on the modifier scales in proportion to the number of germline cell divisions, and increases with the number of loci affected. Loci that do not make a difference to marginal fitnesses at equilibrium are not subject to the reduction principle, and under fine tuning of mutation rates would be expected to have higher mutation rates than loci in mutation-selection balance. Other results include the nonexistence of 'viability analogous, Hardy-Weinberg' modifier polymorphisms under multiplicative mutation, and the sufficiency of average transmission rates to encapsulate the effect of modifier polymorphisms on the transmission of loci under selection. A conjecture is offered regarding situations, like recombination in the presence of mutation, that exhibit departures from the reduction principle. Constraints for tractability are: tight linkage of all loci, initial fixation at the modifier locus, and mutation distributions comprising transition probabilities of reversible Markov chains.Comment: v3: Final corrections. v2: Revised title, reworked and expanded introductory and discussion sections, added corollaries, new results on modifier polymorphisms, minor corrections. 49 pages, 64 reference

Does Mutation Rate Depend on Itself?

Recent evidence suggests that mutation rates are fitness-dependent, broadening our view of the impacts of mutation on the genetic health of populations

Structure and evolutionary history of a large family of NLR proteins in the zebrafish

Multicellular eukaryotes have evolved a range of mechanisms for immune recognition. A widespread family involved in innate immunity are the NACHT-domain and leucine-rich-repeat-containing (NLR) proteins. Mammals have small numbers of NLR proteins, whereas in some species, mostly those without adaptive immune systems, NLRs have expanded into very large families. We describe a family of nearly 400 NLR proteins encoded in the zebrafish genome. The proteins share a defining overall structure, which arose in fishes after a fusion of the core NLR domains with a B30.2 domain, but can be subdivided into four groups based on their NACHT domains. Gene conversion acting differentially on the NACHT and B30.2 domains has shaped the family and created the groups. Evidence of positive selection in the B30.2 domain indicates that this domain rather than the leucine-rich repeats acts as the pathogen recognition module. In an unusual chromosomal organization, the majority of the genes are located on one chromosome arm, interspersed with other large multigene families, including a new family encoding zinc-finger proteins. The NLR-B30.2 proteins represent a new family with diversity in the specific recognition module that is present in fishes in spite of the parallel existence of an adaptive immune system

Complex vaccination strategies prevent the emergence of vaccine resistance

Vaccination is the most effective tool to control infectious diseases. However, the evolution of vaccine resistance, exemplified by vaccine-resistance in SARS-CoV-2, remains a concern. Here, we model complex vaccination strategies against a pathogen with multiple epitopes - molecules targeted by the vaccine. We found that a vaccine targeting one epitope was ineffective in preventing vaccine escape. Vaccine resistance in highly infectious pathogens was prevented by the full-epitope vaccine, that is, one targeting all available epitopes, but only when the rate of pathogen evolution was low. Strikingly, a bet-hedging strategy of random administration of vaccines targeting different epitopes was the most effective in preventing vaccine resistance in pathogens with low rate of infection and high rate of evolution. Thus, complex vaccination strategies, when biologically feasible, may be preferable to the currently used single-vaccine approaches for long-term control of disease outbreaks, especially when applied to livestock with near 100% vaccination rates

Which solar EUV indices are best for reconstructing the solar EUV irradiance ?

The solar EUV irradiance is of key importance for space weather. Most of the time, however, surrogate quantities such as EUV indices have to be used by lack of continuous and spectrally resolved measurements of the irradiance. The ability of such proxies to reproduce the irradiance from different solar atmospheric layers is usually investigated by comparing patterns of temporal correlations. We consider instead a statistical approach. The TIMED/SEE experiment, which has been continuously operating since Feb. 2002, allows for the first time to compare in a statistical manner the EUV spectral irradiance to five EUV proxies: the sunspot number, the f10.7, Ca K, and Mg II indices, and the He I equivalent width. Using multivariate statistical methods such as multidimensional scaling, we represent in a single graph the measure of relatedness between these indices and various strong spectral lines. The ability of each index to reproduce the EUV irradiance is discussed; it is shown why so few lines can be effectively reconstructed from them. All indices exhibit comparable performance, apart from the sunspot number, which is the least appropriate. No single index can satisfactorily describe both the level of variability on time scales beyond 27 days, and relative changes of irradiance on shorter time scales.Comment: 6 figures, to appear in Adv. Space. Re

Analytical study of the effect of recombination on evolution via DNA shuffling

We investigate a multi-locus evolutionary model which is based on the DNA shuffling protocol widely applied in \textit{in vitro} directed evolution. This model incorporates selection, recombination and point mutations. The simplicity of the model allows us to obtain a full analytical treatment of both its dynamical and equilibrium properties, for the case of an infinite population. We also briefly discuss finite population size corrections

Small world effects in evolution

For asexual organisms point mutations correspond to local displacements in the genotypic space, while other genotypic rearrangements represent long-range jumps. We investigate the spreading properties of an initially homogeneous population in a flat fitness landscape, and the equilibrium properties on a smooth fitness landscape. We show that a small-world effect is present: even a small fraction of quenched long-range jumps makes the results indistinguishable from those obtained by assuming all mutations equiprobable. Moreover, we find that the equilibrium distribution is a Boltzmann one, in which the fitness plays the role of an energy, and mutations that of a temperature.Comment: 13 pages and 5 figures. New revised versio

Quasispecies Theory for Horizontal Gene Transfer and Recombination

We introduce a generalization of the parallel, or Crow-Kimura, and Eigen models of molecular evolution to represent the exchange of genetic information between individuals in a population. We study the effect of different schemes of genetic recombination on the steady-state mean fitness and distribution of individuals in the population, through an analytic field theoretic mapping. We investigate both horizontal gene transfer from a population and recombination between pairs of individuals. Somewhat surprisingly, these nonlinear generalizations of quasi-species theory to modern biology are analytically solvable. For two-parent recombination, we find two selected phases, one of which is spectrally rigid. We present exact analytical formulas for the equilibrium mean fitness of the population, in terms of a maximum principle, which are generally applicable to any permutation invariant replication rate function. For smooth fitness landscapes, we show that when positive epistatic interactions are present, recombination or horizontal gene transfer introduces a mild load against selection. Conversely, if the fitness landscape exhibits negative epistasis, horizontal gene transfer or recombination introduce an advantage by enhancing selection towards the fittest genotypes. These results prove that the mutational deterministic hypothesis holds for quasi-species models. For the discontinuous single sharp peak fitness landscape, we show that horizontal gene transfer has no effect on the fitness, while recombination decreases the fitness, for both the parallel and the Eigen models. We present numerical and analytical results as well as phase diagrams for the different cases.Comment: 54 pages; 8 figures; 12 tables; some typos corrected; to appear in Phys. Rev.

Twisting Flux Tubes as a cause of Micro-Flaring Activity

High-cadence optical observations of an H-alpha blue-wing bright point near solar AR NOAA 10794 are presented. The data were obtained with the Dunn Solar Telescope at the National Solar Observatory/Sacramento Peak using a newly developed camera system, the Rapid Dual Imager. Wavelet analysis is undertaken to search for intensity-related oscillatory signatures, and periodicities ranging from 15 to 370 s are found with significance levels exceeding 95%. During two separate microflaring events, oscillation sites surrounding the bright point are observed to twist. We relate the twisting of the oscillation sites to the twisting of physical flux tubes, thus giving rise to reconnection phenomena. We derive an average twist velocity of 8.1 km/s and detect a peak in the emitted flux between twist angles of 180 and 230 degrees.Comment: 8 pages, 10 figure