Recombination and selection trajectory of the monkeypox virus during its adaptation in the human population.

Monkeypox, caused by the monkeypox virus (MPXV), was historically confined to West and Central Africa but has now spread globally. Recombination and selection play crucial roles in the evolutionary adaptation of MPXV; however, the evolution of MPXV and its relationship with the recent, ground-breaking monkeypox epidemic remains poorly understood. To gain insights into the evolutionary dynamics of MPXV, comprehensive in silico recombination and selection analyses were conducted based on MPXV whole genome sequence data. Three types of recombination were identified: five ancestor-sharing interspecies recombination events, six specific interspecies recombination events and four intraspecies recombination events. The results highlight the prevalent occurrence of recombination in MPXV, with 73.3% occurring in variable regions of the genome. Selection analysis was performed from three dimensions: proteins around recombination regions, proteins from recombinant ancestors and MPXV branches, and whole-genome gene analysis. Results revealed 2 and 7 proteins under positive selection in the first two dimensions, respectively. These proteins are mainly involved in infection immunity, apoptosis regulation and viral virulence. Whole-genome analysis detected 25 genes under positive selection, mainly associated with immune response and viral regulation. Understanding their evolutionary patterns will help predict and prevent cross-species transmission, zoonotic outbreaks and potential human epidemics.

Genetic diversity accelerates canine distemper virus adaptation to ferrets.

RNA viruses adapt rapidly to new host environments by generating highly diverse genome sets, so-called "quasispecies." Minor genetic variants promote their rapid adaptation, allowing for the emergence of drug-resistance or immune-escape mutants. Understanding these adaptation processes is highly relevant to assessing the risk of cross-species transmission and the safety and efficacy of vaccines and antivirals. We hypothesized that genetic memory within a viral genome population facilitates rapid adaptation. To test this, we investigated the adaptation of the Morbillivirus canine distemper virus to ferrets and compared an attenuated, Vero cell-adapted virus isolate with its recombinant derivative over consecutive ferret passages. Although both viruses adapted to the new host, the reduced initial genetic diversity of the recombinant virus resulted in delayed disease onset. The non-recombinant virus gradually increased the frequencies of beneficial mutations already present at very low frequencies in the input virus. In contrast, the recombinant virus first evolved de novo mutations to compensate for the initial fitness impairments. Importantly, while both viruses evolved different sets of mutations, most mutations found in the adapted non-recombinant virus were identical to those found in a previous ferret adaptation experiment with the same isolate, indicating that mutations present at low frequency in the original virus stock serve as genetic memory. An arginine residue at position 519 in the carboxy terminus of the nucleoprotein shared by all adapted viruses was found to contribute to pathogenesis in ferrets. Our work illustrates the importance of genetic diversity for adaptation to new environments and identifies regions with functional relevance.IMPORTANCEWhen viruses encounter a new host, they can rapidly adapt to this host and cause disease. How these adaptation processes occur remains understudied. Morbilliviruses have high clinical and veterinary relevance and are attractive model systems to study these adaptation processes. The canine distemper virus is of particular interest, as it exhibits a broader host range than other morbilliviruses and frequently crosses species barriers. Here, we compared the adaptation of an attenuated virus and its recombinant derivative to that of ferrets. Pre-existing mutations present at low frequency allowed faster adaptation of the non-recombinant virus compared to the recombinant virus. We identified a common point mutation in the nucleoprotein that affected the pathogenesis of both viruses. Our study shows that genetic memory facilitates environmental adaptation and that erasing this genetic memory by genetic engineering results in delayed and different adaptation to new environments, providing an important safety aspect for the generation of live-attenuated vaccines.

Evolution across the adaptive landscape in a hyperdiverse beetle radiation.

The extraordinary diversification of beetles on Earth is a textbook example of adaptive evolution. Yet, the tempo and drivers of this super-radiation remain largely unclear. Here, we address this problem by investigating macroevolutionary dynamics in darkling beetles (Coleoptera: Tenebrionidae), one of the most ecomorphologically diverse beetle families (with over 30,000 species). Using multiple genomic datasets and analytical approaches, we resolve the long-standing inconsistency over deep relationships in the family. In conjunction with a landmark-based dataset of body shape morphology, we show that the evolutionary history of darkling beetles is marked by ancient rapid radiations, frequent ecological transitions, and rapid bursts of morphological diversification. On a global scale, our analyses uncovered a significant pulse of phenotypic diversification proximal to the Cretaceous-Palaeogene (K/Pg) mass extinction and convergence of body shape associated with recurrent ecological specializations. On a regional scale, two major Australasian radiations, the Adeliini and the Heleine clade, exhibited contrasting patterns of ecomorphological diversification, representing phylogenetic niche conservatism versus adaptive radiation. Our findings align with the Simpsonian model of adaptive evolution across the macroevolutionary landscape and highlight a significant role of ecological opportunity in driving the immense ecomorphological diversity in a hyperdiverse beetle group.

Genome sequencing of Coryphaenoides yaquinae reveals convergent and lineage-specific molecular evolution in deep-sea adaptation.

Abyssal (3501-6500 m) and hadal (>6500 m) fauna evolve under harsh abiotic stresses, characterized by high hydrostatic pressure, darkness and food shortage, providing unique opportunities to investigate mechanisms underlying environmental adaptation. Genomes of several hadal species have recently been reported. However, the genetic adaptation of deep sea species across a broad spectrum of ocean depths has yet to be thoroughly investigated, due to the challenges imposed by collecting the deep sea species. To elucidate the correlation between genetic innovation and vertical distribution, we generated a chromosome-level genome assembly of the macrourids Coryphaenoides yaquinae, which is widely distributed in the abyssal/hadal zone ranging from 3655 to 7259 m in depth. Genomic comparisons among shallow, abyssal and hadal-living species identified idiosyncratic and convergent genetic alterations underlying the extraordinary adaptations of deep-sea species including light perception, circadian regulation, hydrostatic pressure and hunger tolerance. The deep-sea fishes (Coryphaenoides Sp. and Pseudoliparis swirei) venturing into various ocean depths independently have undergone convergent amino acid substitutions in multiple proteins such as rhodopsin 1, pancreatic and duodenal homeobox 1 and melanocortin 4 receptor which are known or verified in zebrafish to be related with vision adaptation and energy expenditure. Convergent evolution events were also identified in heat shock protein 90 beta family member 1 and valosin-containing protein genes known to be related to hydrostatic pressure adaptation specifically in fishes found around the hadal range. The uncovering of the molecular convergence among the deep-sea species shed new light on the common genetic innovations required for deep-sea adaptation by the fishes.

Rapid Adaptation and Interspecific Introgression in the North American Crop Pest Helicoverpa zea.

Insect crop pests threaten global food security. This threat is amplified through the spread of nonnative species and through adaptation of native pests to control measures. Adaptations such as pesticide resistance can result from selection on variation within a population, or through gene flow from another population. We investigate these processes in an economically important noctuid crop pest, Helicoverpa zea, which has evolved resistance to a wide range of pesticides. Its sister species Helicoverpa armigera, first detected as an invasive species in Brazil in 2013, introduced the pyrethroid-resistance gene CYP337B3 to South American H. zea via adaptive introgression. To understand whether this could contribute to pesticide resistance in North America, we sequenced 237 H. zea genomes across 10 sample sites. We report H. armigera introgression into the North American H. zea population. Two individuals sampled in Texas in 2019 carry H. armigera haplotypes in a 4 Mbp region containing CYP337B3. Next, we identify signatures of selection in the panmictic population of nonadmixed H. zea, identifying a selective sweep at a second cytochrome P450 gene: CYP333B3. We estimate that its derived allele conferred a ∼5% fitness advantage and show that this estimate explains independently observed rare nonsynonymous CYP333B3 mutations approaching fixation over a ∼20-year period. We also detect putative signatures of selection at a kinesin gene associated with Bt resistance. Overall, we document two mechanisms of rapid adaptation: the introduction of fitness-enhancing alleles through interspecific introgression, and selection on intraspecific variation.

Conditionally Deleterious Mutation Load Accumulates in Genomic Islands of Local Adaptation but Can Be Purged with Sufficient Genotypic Redundancy.

AbstractLocal adaptation frequently evolves in patches or environments that are connected via migration. In these cases, genomic regions that are linked to a locally adapted locus experience reduced effective migration rates. Via individual-based simulations of a two-patch system, we show that this reduced effective migration results in the accumulation of conditionally deleterious mutations, but not universally deleterious mutations, adjacent to adaptive loci. When there is redundancy in the genetic basis of local adaptation (i.e., genotypic redundancy), turnover of locally adapted polymorphisms allows conditionally deleterious mutation load to be purged. The amount of mutational load that accumulates adjacent to locally adapted loci is dependent on redundancy, recombination rate, migration rate, population size, strength of selection, and the phenotypic effect size of adaptive alleles. Our results highlight the need to be cautious when interpreting patterns of local adaptation at the level of phenotype or fitness, as the genetic basis of local adaptation can be transient, and evolution may confer a degree of maladaptation to nonlocal environments.

Predictable and Divergent Change in the Multivariate P Matrix during Parallel Adaptation.

AbstractAdaptation to replicated environmental conditions can be remarkably predictable, suggesting that parallel evolution may be a common feature of adaptive radiation. An open question, however, is how phenotypic variation itself evolves during repeated adaptation. Here, we use a dataset of morphological measurements from 35 populations of threespine stickleback, consisting of 16 parapatric lake-stream pairs and three marine populations, to understand how phenotypic variation has evolved during transitions from marine to freshwater environments and during subsequent diversification across the lake-stream boundary. We find statistical support for divergent phenotypic covariance (P) across populations, with most diversification of P occurring among freshwater populations. Despite a close correspondence between within-population phenotypic variation and among-population divergence, we find that variation in P is unrelated to total variation in population means across the set of populations. For lake-stream pairs, we find that theoretical predictions for microevolutionary change can explain more than 30% of divergence in P matrices across the habitat boundary. Together, our results indicate that divergence in variance structure occurs primarily in dimensions of trait space with low phenotypic integration, correlated with disparate lake and stream environments. Our findings illustrate how conserved and divergent features of multivariate variation can underlie adaptive radiation.

Understanding evolutionary rescue and parallelism in response to environmental stress.

Evolutionary rescue, the process by which populations facing environmental stress avoid extinction through genetic adaptation, is a critical area of study in evolutionary biology. The order in which mutations arise and get established will be relevant to the population's rescue. This study investigates the degree of parallel evolution at the genotypic level between independent populations facing environmental stress and subject to different demographic regimes. Under density regulation, 2 regimes exist: In the first, the population can restore positive growth rates by adjusting its population size or through adaptive mutations, whereas in the second regime, the population is doomed to extinction unless a rescue mutation occurs. Analytical approximations for the likelihood of evolutionary rescue are obtained and contrasted with simulation results. We show that the initial level of maladaptation and the demographic regime significantly affect the level of parallelism. There is an evident transition between these 2 regimes. Whereas in the first regime, parallelism decreases with the level of maladaptation, it displays the opposite behavior in the rescue/extinction regime. These findings have important implications for understanding population persistence and the degree of parallelism in evolutionary responses as they integrate demographic effects and evolutionary processes.

Non-Poissonian Bursts in the Arrival of Phenotypic Variation Can Strongly Affect the Dynamics of Adaptation.

Modeling the rate at which adaptive phenotypes appear in a population is a key to predicting evolutionary processes. Given random mutations, should this rate be modeled by a simple Poisson process, or is a more complex dynamics needed? Here we use analytic calculations and simulations of evolving populations on explicit genotype-phenotype maps to show that the introduction of novel phenotypes can be "bursty" or overdispersed. In other words, a novel phenotype either appears multiple times in quick succession or not at all for many generations. These bursts are fundamentally caused by statistical fluctuations and other structure in the map from genotypes to phenotypes. Their strength depends on population parameters, being highest for "monomorphic" populations with low mutation rates. They can also be enhanced by additional inhomogeneities in the mapping from genotypes to phenotypes. We mainly investigate the effect of bursts using the well-studied genotype-phenotype map for RNA secondary structure, but find similar behavior in a lattice protein model and in Richard Dawkins's biomorphs model of morphological development. Bursts can profoundly affect adaptive dynamics. Most notably, they imply that fitness differences play a smaller role in determining which phenotype fixes than would be the case for a Poisson process without bursts.

Skin colour: A window into human phenotypic evolution and environmental adaptation.

As modern humans ventured out of Africa and dispersed around the world, they faced novel environmental challenges that led to geographic adaptations including skin colour. Over the long history of human evolution, skin colour has changed dramatically, showing tremendous diversity across different geographical regions, for example, the majority of individuals from the expansive lands of Africa have darker skin, whereas the majority of people from Eurasia exhibit lighter skin. What adaptations did lighter skin confer upon modern humans as they migrated from Africa to Eurasia? What genetic mechanisms underlie the diversity of skin colour observed in different populations? In recent years, scientists have gradually gained a deeper understanding of the interactions between pigmentation gene and skin colour through population-based genomic studies of different groups around the world, particularly in East Asia and Africa. In this review, we summarize our current understanding of 26 skin colour-related pigmentation genes and 48 SNPs that influence skin colour. Important pigmentation genes across three major populations are described in detail: MFSD12, SLC24A5, PDPK1 and DDB1/CYB561A3/TMEM138 influence skin colour in African populations; OCA2, KITLG, SLC24A2, GNPAT and PAH are key to the evolution of skin pigmentation in East Asian populations; and SLC24A5, SLC45A2, TYR, TYRP1, ASIP, MC1R and IRF4 significantly contribute to the lightening of skin colour in European populations. We summarized recent findings in genomic studies of skin colour in populations that implicate diverse geographic environments, local adaptation among populations, gene flow and multi-gene interactions as factors influencing skin colour diversity.

Epistasis, core-genome disharmony, and adaptation in recombining bacteria.

Recombination of short DNA fragments via horizontal gene transfer (HGT) can introduce beneficial alleles, create genomic disharmony through negative epistasis, and create adaptive gene combinations through positive epistasis. For non-core (accessory) genes, the negative epistatic cost is likely to be minimal because the incoming genes have not co-evolved with the recipient genome and are frequently observed as tightly linked cassettes with major effects. By contrast, interspecific recombination in the core genome is expected to be rare because disruptive allelic replacement is likely to introduce negative epistasis. Why then is homologous recombination common in the core of bacterial genomes? To understand this enigma, we take advantage of an exceptional model system, the common enteric pathogens Campylobacter jejuni and C. coli that are known for very high magnitude interspecies gene flow in the core genome. As expected, HGT does indeed disrupt co-adapted allele pairings, indirect evidence of negative epistasis. However, multiple HGT events enable recovery of the genome's co-adaption between introgressing alleles, even in core metabolism genes (e.g., formate dehydrogenase). These findings demonstrate that, even for complex traits, genetic coalitions can be decoupled, transferred, and independently reinstated in a new genetic background-facilitating transition between fitness peaks. In this example, the two-step recombinational process is associated with C. coli that are adapted to the agricultural niche.IMPORTANCEGenetic exchange among bacteria shapes the microbial world. From the acquisition of antimicrobial resistance genes to fundamental questions about the nature of bacterial species, this powerful evolutionary force has preoccupied scientists for decades. However, the mixing of genes between species rests on a paradox: 0n one hand, promoting adaptation by conferring novel functionality; on the other, potentially introducing disharmonious gene combinations (negative epistasis) that will be selected against. Taking an interdisciplinary approach to analyze natural populations of the enteric bacteria Campylobacter, an ideal example of long-range admixture, we demonstrate that genes can independently transfer across species boundaries and rejoin in functional networks in a recipient genome. The positive impact of two-gene interactions appears to be adaptive by expanding metabolic capacity and facilitating niche shifts through interspecific hybridization. This challenges conventional ideas and highlights the possibility of multiple-step evolution of multi-gene traits by interspecific introgression.

Limits to species' range: the tension between local and global adaptation.

We know that heritable variation is abundant, and that selection causes all but the smallest populations to rapidly shift beyond their original trait distribution. So then, what limits the range of a species? There are physical constraints and also population genetic limits to the effectiveness of selection, ultimately set by population size. Global adaptation, where the same genotype is favoured over the whole range, is most efficient when based on a multitude of weakly selected alleles and is effective even when local demes are small, provided that there is some gene flow. In contrast, local adaptation is sensitive to gene flow and may require alleles with substantial effect. How can populations combine the advantages of large effective size with the ability to specialise into local niches? To what extent does reproductive isolation help resolve this tension? I address these questions using eco-evolutionary models of polygenic adaptation, contrasting discrete demes with continuousspace.

Design Principles for Biological Adaptation: A Systems and Control-Theoretic Treatment.

Establishing a mapping between (from and to) the functionality of interest and the underlying network structure (design principles) remains a crucial step toward understanding and design of bio-systems. Perfect adaptation is one such crucial functionality that enables every living organism to regulate its essential activities in the presence of external disturbances. Previous approaches to deducing the design principles for adaptation have either relied on computationally burdensome brute-force methods or rule-based design strategies detecting only a subset of all possible adaptive network structures. This chapter outlines a scalable and generalizable method inspired by systems theory that unravels an exhaustive set of adaptation-capable structures. We first use the well-known performance parameters to characterize perfect adaptation. These performance parameters are then mapped back to a few parameters (poles, zeros, gain) characteristic of the underlying dynamical system constituted by the rate equations. Therefore, the performance parameters evaluated for the scenario of perfect adaptation can be expressed as a set of precise mathematical conditions involving the system parameters. Finally, we use algebraic graph theory to translate these abstract mathematical conditions to certain structural requirements for adaptation. The proposed algorithm does not assume any particular dynamics and is applicable to networks of any size. Moreover, the results offer a significant advancement in the realm of understanding and designing complex biochemical networks.

Probabilistic measures for biological adaptation and resilience.

This paper introduces an approach to quantifying ecological resilience in biological systems, particularly focusing on noisy systems responding to episodic disturbances with sudden adaptations. Incorporating concepts from nonequilibrium statistical mechanics, we propose a measure termed "ecological resilience through adaptation," specifically tailored to noisy, forced systems that undergo physiological adaptation in the face of stressful environmental changes. Randomness plays a key role, accounting for model uncertainty and the inherent variability in the dynamical response among components of biological systems. Our measure of resilience is rooted in the probabilistic description of states within these systems and is defined in terms of the dynamics of the ensemble average of a model-specific observable quantifying success or well-being. Our approach utilizes stochastic linear response theory to compute how the expected success of a system, originally in statistical equilibrium, dynamically changes in response to a environmental perturbation and a subsequent adaptation. The resulting mathematical derivations allow for the estimation of resilience in terms of ensemble averages of simulated or experimental data. Finally, through a simple but clear conceptual example, we illustrate how our resilience measure can be interpreted and compared to other existing frameworks in the literature. The methodology is general but inspired by applications in plant systems, with the potential for broader application to complex biological processes.

Fitness as the organismal performance measure guiding adaptive evolution.

A long-standing problem in evolutionary theory is to clarify in what sense (if any) natural selection cumulatively improves the design of organisms. Various concepts, such as fitness and inclusive fitness, have been proposed to resolve this problem. In addition, there have been attempts to replace the original problem with more tractable questions, such as whether a given gene or trait is favored by selection. Here, we ask what theoretical properties the concept fitness should possess to encapsulate the improvement criterion required to talk meaningfully about adaptive evolution. We argue that natural selection tends to shape phenotypes based on the causal properties of individuals and that this tendency is, therefore, best captured by a fitness concept that focuses on these properties. We highlight a fitness concept that meets this role under broad conditions but requires adjustments in our conceptual understanding of adaptive evolution. These adjustments combine elements of Dawkinsian gene selectionism and Egbert Leigh's "parliament of genes."

Linking morphology, performance, and habitat utilization: adaptation across biologically relevant 'levels' in tamarins.

BACKGROUND: Biological adaptation manifests itself at the interface of different biologically relevant 'levels', such as ecology, performance, and morphology. Integrated studies at this interface are scarce due to practical difficulties in study design. We present a multilevel analysis, in which we combine evidence from habitat utilization, leaping performance and limb bone morphology of four species of tamarins to elucidate correlations between these 'levels'. RESULTS: We conducted studies of leaping behavior in the field and in a naturalistic park and found significant differences in support use and leaping performance. Leontocebus nigrifrons leaps primarily on vertical, inflexible supports, with vertical body postures, and covers greater leaping distances on average. In contrast, Saguinus midas and S. imperator use vertical and horizontal supports for leaping with a relatively similar frequency. S. mystax is similar to S. midas and S. imperator in the use of supports, but covers greater leaping distances on average, which are nevertheless shorter than those of L. nigrifrons. We assumed these differences to be reflected in the locomotor morphology, too, and compared various morphological features of the long bones of the limbs. According to our performance and habitat utilization data, we expected the long bone morphology of L. nigrifrons to reflect the largest potential for joint torque generation and stress resistance, because we assume longer leaps on vertical supports to exert larger forces on the bones. For S. mystax, based on our performance data, we expected the potential for torque generation to be intermediate between L. nigrifrons and the other two Saguinus species. Surprisingly, we found S. midas and S. imperator having relatively more robust morphological structures as well as relatively larger muscle in-levers, and thus appearing better adapted to the stresses involved in leaping than the other two. CONCLUSION: This study demonstrates the complex ways in which behavioral and morphological 'levels' map onto each other, cautioning against oversimplification of ecological profiles when using large interspecific eco-morphological studies to make adaptive evolutionary inferences.

Convergence in carnivorous pitcher plants reveals a mechanism for composite trait evolution.

Composite traits involve multiple components that, only when combined, gain a new synergistic function. Thus, how they evolve remains a puzzle. We combined field experiments, microscopy, chemical analyses, and laser Doppler vibrometry with comparative phylogenetic analyses to show that two carnivorous Nepenthes pitcher plant species independently evolved similar adaptations in three distinct traits to acquire a new, composite trapping mechanism. Comparative analyses suggest that this new trait arose convergently through "spontaneous coincidence" of the required trait combination, rather than directional selection in the component traits. Our results indicate a plausible mechanism for composite trait evolution and highlight the importance of stochastic phenotypic variation as a facilitator of evolutionary novelty.

Dynamic Phytomeric Growth Contributes to Local Adaptation in Barley.

Vascular plants have segmented body axes with iterative nodes and internodes. Appropriate node initiation and internode elongation are fundamental to plant fitness and crop yield; however, how these events are spatiotemporally coordinated remains elusive. We show that in barley (Hordeum vulgare L.), selections during domestication have extended the apical meristematic phase to promote node initiation, but constrained subsequent internode elongation. In both vegetative and reproductive phases, internode elongation displays a dynamic proximal-distal gradient, and among subpopulations of domesticated barleys worldwide, node initiation and proximal internode elongation are associated with latitudinal and longitudinal gradients, respectively. Genetic and functional analyses suggest that, in addition to their converging roles in node initiation, flowering-time genes have been repurposed to specify the timing and duration of internode elongation. Our study provides an integrated view of barley node initiation and internode elongation and suggests that plant architecture should be recognized as a collection of dynamic phytomeric units in the context of crop adaptive evolution.

The evolutionary dynamics of local adaptations under genetic rescue is determined by mutational load and polygenicity.

Inbred populations often suffer from increased mutational load and reduced fitness due to lower efficacy of purifying selection in groups with small effective population sizes. Genetic rescue (GR) is a conservation tool that is studied and deployed with the aim of increasing the fitness of such inbred populations by assisted migration of individuals from closely related outbred populations. The success of GR depends on several factors-such as their demographic history and distribution of dominance effects of mutations-that may vary across populations. While we understand the impact of these factors on the dynamics of GR, their impact on local adaptations remains unclear. To this end, we conduct a population genetics simulation study to evaluate the impact of trait complexity (Mendelian vs. polygenic), dominance effects, and demographic history on the efficacy of GR. We find that the impact on local adaptations depends highly on the mutational load at the time of GR, which is in turn shaped dynamically by interactions between demographic history and dominance effects of deleterious variation. Over time local adaptations are generally restored post-GR, though in the short term they are often compromised in the process of purging deleterious variation. We also show that while local adaptations are almost always fully restored, the degree to which ancestral genetic variation affecting the trait is replaced by donor variation can vary drastically and is especially high for complex traits. Our results provide insights on the impact of GR on trait evolution and considerations for the practical implementation of GR.