Positive selection underlies repeated knockout of ORF8 in SARS-CoV-2 evolution.

Knockout of the ORF8 protein has repeatedly spread through the global viral population during SARS-CoV-2 evolution. Here we use both regional and global pathogen sequencing to explore the selection pressures underlying its loss. In Washington State, we identified transmission clusters with ORF8 knockout throughout SARS-CoV-2 evolution, not just on novel, high fitness viral backbones. Indeed, ORF8 is truncated more frequently and knockouts circulate for longer than for any other gene. Using a global phylogeny, we find evidence of positive selection to explain this phenomenon: nonsense mutations resulting in shortened protein products occur more frequently and are associated with faster clade growth rates than synonymous mutations in ORF8. Loss of ORF8 is also associated with reduced clinical severity, highlighting the diverse clinical impacts of SARS-CoV-2 evolution.

Hybrid speciation driven by multilocus introgression of ecological traits.

Hybridization allows adaptations to be shared among lineages and may trigger the evolution of new species1,2. However, convincing examples of homoploid hybrid speciation remain rare because it is challenging to demonstrate that hybridization was crucial in generating reproductive isolation3. Here we combine population genomic analysis with quantitative trait locus mapping of species-specific traits to examine a case of hybrid speciation in Heliconius butterflies. We show that Heliconius elevatus is a hybrid species that is sympatric with both parents and has persisted as an independently evolving lineage for at least 180,000 years. This is despite pervasive and ongoing gene flow with one parent, Heliconius pardalinus, which homogenizes 99% of their genomes. The remaining 1% introgressed from the other parent, Heliconius melpomene, and is scattered widely across the H. elevatus genome in islands of divergence from H. pardalinus. These islands contain multiple traits that are under disruptive selection, including colour pattern, wing shape, host plant preference, sex pheromones and mate choice. Collectively, these traits place H. elevatus on its own adaptive peak and permit coexistence with both parents. Our results show that speciation was driven by introgression of ecological traits, and that speciation with gene flow is possible with a multilocus genetic architecture.

Natural Selection Signatures in the Hondo and Ryukyu Japanese Subpopulations.

Natural selection signatures across Japanese subpopulations are under-explored. Here we conducted genome-wide selection scans with 622,926 single nucleotide polymorphisms for 20,366 Japanese individuals, who were recruited from the main-islands of Japanese Archipelago (Hondo) and the Ryukyu Archipelago (Ryukyu), representing two major Japanese subpopulations. The integrated haplotype score (iHS) analysis identified several signals in one or both subpopulations. We found a novel candidate locus at IKZF2, especially in Ryukyu. Significant signals were observed in the major histocompatibility complex region in both subpopulations. The lead variants differed and demonstrated substantial allele frequency differences between Hondo and Ryukyu. The lead variant in Hondo tags HLA-A*33:03-C*14:03-B*44:03-DRB1*13:02-DQB1*06:04-DPB1*04:01, a haplotype specific to Japanese and Korean. While in Ryukyu, the lead variant tags DRB1*15:01-DQB1*06:02, which had been recognized as a genetic risk factor for narcolepsy. In contrast, it is reported to confer protective effects against type 1 diabetes and human T lymphotropic virus type 1-associated myelopathy/tropical spastic paraparesis. The FastSMC analysis identified 8 loci potentially affected by selection within the past 20-150 generations, including 2 novel candidate loci. The analysis also showed differences in selection patterns of ALDH2 between Hondo and Ryukyu, a gene recognized to be specifically targeted by selection in East Asian. In summary, our study provided insights into the selection signatures within the Japanese and nominated potential sources of selection pressure.

High-resolution mapping reveals hotspots and sex-biased recombination in Populus trichocarpa.

Fine-scale meiotic recombination is fundamental to the outcome of natural and artificial selection. Here, dense genetic mapping and haplotype reconstruction were used to estimate recombination for a full factorial Populus trichocarpa cross of 7 males and 7 females. Genomes of the resulting 49 full-sib families (N = 829 offspring) were resequenced, and high-fidelity biallelic SNP/INDELs and pedigree information were used to ascertain allelic phase and impute progeny genotypes to recover gametic haplotypes. The 14 parental genetic maps contained 1,820 SNP/INDELs on average that covered 376.7 Mb of physical length across 19 chromosomes. Comparison of parental and progeny haplotypes allowed fine-scale demarcation of cross-over regions, where 38,846 cross-over events in 1,658 gametes were observed. Cross-over events were positively associated with gene density and negatively associated with GC content and long-terminal repeats. One of the most striking findings was higher rates of cross-overs in males in 8 out of 19 chromosomes. Regions with elevated male cross-over rates had lower gene density and GC content than windows showing no sex bias. High-resolution analysis identified 67 candidate cross-over hotspots spread throughout the genome. DNA sequence motifs enriched in these regions showed striking similarity to those of maize, Arabidopsis, and wheat. These findings, and recombination estimates, will be useful for ongoing efforts to accelerate domestication of this and other biomass feedstocks, as well as future studies investigating broader questions related to evolutionary history, perennial development, phenology, wood formation, vegetative propagation, and dioecy that cannot be studied using annual plant model systems.

Host adaptation of codon usage in SARS-CoV-2 from mammals indicates potential natural selection and viral fitness.

SARS-CoV-2 infection, which is the cause of the COVID-19 pandemic, has expanded across various animal hosts, and the virus can be transmitted particularly efficiently in minks. It is still not clear how SARS-CoV-2 is selected and evolves in its hosts, or how mutations affect viral fitness. In this report, sequences of SARS-CoV-2 isolated from human and animal hosts were analyzed, and the binding energy and capacity of the spike protein to bind human ACE2 and the mink receptor were compared. Codon adaptation index (CAI) analysis indicated the optimization of viral codons in some animals such as bats and minks, and a neutrality plot demonstrated that natural selection had a greater influence on some SARS-CoV-2 sequences than mutational pressure. Molecular dynamics simulation results showed that the mutations Y453F and N501T in mink SARS-CoV-2 could enhance the binding of the viral spike to the mink receptor, indicating the involvement of these mutations in natural selection and viral fitness. Receptor binding analysis revealed that the mink SARS-CoV-2 spike interacted more strongly with the mink receptor than the human receptor. Tracking the variations and codon bias of SARS-CoV-2 is helpful for understanding the fitness of the virus in virus transmission, pathogenesis, and immune evasion.

Pandemic-scale phylogenomics reveals the SARS-CoV-2 recombination landscape.

Accurate and timely detection of recombinant lineages is crucial for interpreting genetic variation, reconstructing epidemic spread, identifying selection and variants of interest, and accurately performing phylogenetic analyses1-4. During the SARS-CoV-2 pandemic, genomic data generation has exceeded the capacities of existing analysis platforms, thereby crippling real-time analysis of viral evolution5. Here, we use a new phylogenomic method to search a nearly comprehensive SARS-CoV-2 phylogeny for recombinant lineages. In a 1.6 million sample tree from May 2021, we identify 589 recombination events, which indicate that around 2.7% of sequenced SARS-CoV-2 genomes have detectable recombinant ancestry. Recombination breakpoints are inferred to occur disproportionately in the 3' portion of the genome that contains the spike protein. Our results highlight the need for timely analyses of recombination for pinpointing the emergence of recombinant lineages with the potential to increase transmissibility or virulence of the virus. We anticipate that this approach will empower comprehensive real-time tracking of viral recombination during the SARS-CoV-2 pandemic and beyond.

Toward a theory of evolution as multilevel learning.

We apply the theory of learning to physically renormalizable systems in an attempt to outline a theory of biological evolution, including the origin of life, as multilevel learning. We formulate seven fundamental principles of evolution that appear to be necessary and sufficient to render a universe observable and show that they entail the major features of biological evolution, including replication and natural selection. It is shown that these cornerstone phenomena of biology emerge from the fundamental features of learning dynamics such as the existence of a loss function, which is minimized during learning. We then sketch the theory of evolution using the mathematical framework of neural networks, which provides for detailed analysis of evolutionary phenomena. To demonstrate the potential of the proposed theoretical framework, we derive a generalized version of the Central Dogma of molecular biology by analyzing the flow of information during learning (back propagation) and predicting (forward propagation) the environment by evolving organisms. The more complex evolutionary phenomena, such as major transitions in evolution (in particular, the origin of life), have to be analyzed in the thermodynamic limit, which is described in detail in the paper by Vanchurin et al. [V. Vanchurin, Y. I. Wolf, E. V. Koonin, M. I. Katsnelson, Proc. Natl. Acad. Sci. U.S.A. 119, 10.1073/pnas.2120042119 (2022)].

Effect of clonal testing on the efficiency of genomic evaluation in forest tree breeding.

Through stochastic simulations, accuracies of breeding values and response to selection were assessed under traditional pedigree-(BLUP) and genomic-based evaluation methods (GBLUP) in forest tree breeding. The latter provides a methodological foundation for genomic selection. We evaluated the impact of clonal replication in progeny testing on the response to selection realized in seed orchards under variable marker density and target effective population sizes. We found that clonal replication in progeny trials boosted selection accuracy, thus providing additional genetic gains under BLUP. While a similar trend was observed for GBLUP, however, the added gains did not surpass those under BLUP. Therefore, breeding programs deploying extensive progeny testing with clonal propagation might not benefit from the deployment of genomic information. These findings could be helpful in the context of operational breeding programs.

Evolution enhances mutational robustness and suppresses the emergence of a new phenotype: A new computational approach for studying evolution.

The aim of this paper is two-fold. First, we propose a new computational method to investigate the particularities of evolution. Second, we apply this method to a model of gene regulatory networks (GRNs) and explore the evolution of mutational robustness and bistability. Living systems have developed their functions through evolutionary processes. To understand the particularities of this process theoretically, evolutionary simulation (ES) alone is insufficient because the outcomes of ES depend on evolutionary pathways. We need a reference system for comparison. An appropriate reference system for this purpose is an ensemble of the randomly sampled genotypes. However, generating high-fitness genotypes by simple random sampling is difficult because such genotypes are rare. In this study, we used the multicanonical Monte Carlo method developed in statistical physics to construct a reference ensemble of GRNs and compared it with the outcomes of ES. We obtained the following results. First, mutational robustness was significantly higher in ES than in the reference ensemble at the same fitness level. Second, the emergence of a new phenotype, bistability, was delayed in evolution. Third, the bistable group of GRNs contains many mutationally fragile GRNs compared with those in the non-bistable group. This suggests that the delayed emergence of bistability is a consequence of the mutation-selection mechanism.

The evolutionary maintenance of Lévy flight foraging.

Lévy flight is a type of random walk that characterizes the behaviour of many natural phenomena studied across a multiplicity of academic disciplines; within biology specifically, the behaviour of fish, birds, insects, mollusks, bacteria, plants, slime molds, t-cells, and human populations. The Lévy flight foraging hypothesis states that because Lévy flights can maximize an organism's search efficiency, natural selection should result in Lévy-like behaviour. Empirical and theoretical research has provided ample evidence of Lévy walks in both extinct and extant species, and its efficiency across models with a diversity of resource distributions. However, no model has addressed the maintenance of Lévy flight foraging through evolutionary processes, and existing models lack ecological breadth. We use numerical simulations, including lineage-based models of evolution with a distribution of move lengths as a variable and heritable trait, to test the Lévy flight foraging hypothesis. We include biological and ecological contexts such as population size, searching costs, lifespan, resource distribution, speed, and consider both energy accumulated at the end of a lifespan and averaged over a lifespan. We demonstrate that selection often results in Lévy-like behaviour, although conditional; smaller populations, longer searches, and low searching costs increase the fitness of Lévy-like behaviour relative to Brownian behaviour. Interestingly, our results also evidence a bet-hedging strategy; Lévy-like behaviour reduces fitness variance, thus maximizing geometric mean fitness over multiple generations.

Sequence locally, think globally: The Darwin Tree of Life Project.

The goals of the Earth Biogenome Project-to sequence the genomes of all eukaryotic life on earth-are as daunting as they are ambitious. The Darwin Tree of Life Project was founded to demonstrate the credibility of these goals and to deliver at-scale genome sequences of unprecedented quality for a biogeographic region: the archipelago of islands that constitute Britain and Ireland. The Darwin Tree of Life Project is a collaboration between biodiversity organizations (museums, botanical gardens, and biodiversity institutes) and genomics institutes. Together, we have built a workflow that collects specimens from the field, robustly identifies them, performs sequencing, generates high-quality, curated assemblies, and releases these openly for the global community to use to build future science and conservation efforts.

Emergence of methicillin resistance predates the clinical use of antibiotics.

The discovery of antibiotics more than 80 years ago has led to considerable improvements in human and animal health. Although antibiotic resistance in environmental bacteria is ancient, resistance in human pathogens is thought to be a modern phenomenon that is driven by the clinical use of antibiotics1. Here we show that particular lineages of methicillin-resistant Staphylococcus aureus-a notorious human pathogen-appeared in European hedgehogs in the pre-antibiotic era. Subsequently, these lineages spread within the local hedgehog populations and between hedgehogs and secondary hosts, including livestock and humans. We also demonstrate that the hedgehog dermatophyte Trichophyton erinacei produces two ß-lactam antibiotics that provide a natural selective environment in which methicillin-resistant S. aureus isolates have an advantage over susceptible isolates. Together, these results suggest that methicillin resistance emerged in the pre-antibiotic era as a co-evolutionary adaptation of S. aureus to the colonization of dermatophyte-infected hedgehogs. The evolution of clinically relevant antibiotic-resistance genes in wild animals and the connectivity of natural, agricultural and human ecosystems demonstrate that the use of a One Health approach is critical for our understanding and management of antibiotic resistance, which is one of the biggest threats to global health, food security and development.

Mutation bias reflects natural selection in Arabidopsis thaliana.

Since the first half of the twentieth century, evolutionary theory has been dominated by the idea that mutations occur randomly with respect to their consequences1. Here we test this assumption with large surveys of de novo mutations in the plant Arabidopsis thaliana. In contrast to expectations, we find that mutations occur less often in functionally constrained regions of the genome-mutation frequency is reduced by half inside gene bodies and by two-thirds in essential genes. With independent genomic mutation datasets, including from the largest Arabidopsis mutation accumulation experiment conducted to date, we demonstrate that epigenomic and physical features explain over 90% of variance in the genome-wide pattern of mutation bias surrounding genes. Observed mutation frequencies around genes in turn accurately predict patterns of genetic polymorphisms in natural Arabidopsis accessions (r = 0.96). That mutation bias is the primary force behind patterns of sequence evolution around genes in natural accessions is supported by analyses of allele frequencies. Finally, we find that genes subject to stronger purifying selection have a lower mutation rate. We conclude that epigenome-associated mutation bias2 reduces the occurrence of deleterious mutations in Arabidopsis, challenging the prevailing paradigm that mutation is a directionless force in evolution.

PopHumanVar: an interactive application for the functional characterization and prioritization of adaptive genomic variants in humans.

Adaptive challenges that humans faced as they expanded across the globe left specific molecular footprints that can be decoded in our today's genomes. Different sets of metrics are used to identify genomic regions that have undergone selection. However, there are fewer methods capable of pinpointing the allele ultimately responsible for this selection. Here, we present PopHumanVar, an interactive online application that is designed to facilitate the exploration and thorough analysis of candidate genomic regions by integrating both functional and population genomics data currently available. PopHumanVar generates useful summary reports of prioritized variants that are putatively causal of recent selective sweeps. It compiles data and graphically represents different layers of information, including natural selection statistics, as well as functional annotations and genealogical estimations of variant age, for biallelic single nucleotide variants (SNVs) of the 1000 Genomes Project phase 3. Specifically, PopHumanVar amasses SNV-based information from GEVA, SnpEFF, GWAS Catalog, ClinVar, RegulomeDB and DisGeNET databases, as well as accurate estimations of iHS, nSL and iSAFE statistics. Notably, PopHumanVar can successfully identify known causal variants of frequently reported candidate selection regions, including EDAR in East-Asians, ACKR1 (DARC) in Africans and LCT/MCM6 in Europeans. PopHumanVar is open and freely available at https://pophumanvar.uab.cat.

Domestication reshaped the genetic basis of inbreeding depression in a maize landrace compared to its wild relative, teosinte.

Inbreeding depression is the reduction in fitness and vigor resulting from mating of close relatives observed in many plant and animal species. The extent to which the genetic load of mutations contributing to inbreeding depression is due to large-effect mutations versus variants with very small individual effects is unknown and may be affected by population history. We compared the effects of outcrossing and self-fertilization on 18 traits in a landrace population of maize, which underwent a population bottleneck during domestication, and a neighboring population of its wild relative teosinte. Inbreeding depression was greater in maize than teosinte for 15 of 18 traits, congruent with the greater segregating genetic load in the maize population that we predicted from sequence data. Parental breeding values were highly consistent between outcross and selfed offspring, indicating that additive effects determine most of the genetic value even in the presence of strong inbreeding depression. We developed a novel linkage scan to identify quantitative trait loci (QTL) representing large-effect rare variants carried by only a single parent, which were more important in teosinte than maize. Teosinte also carried more putative juvenile-acting lethal variants identified by segregation distortion. These results suggest a mixture of mostly polygenic, small-effect partially recessive effects in linkage disequilibrium underlying inbreeding depression, with an additional contribution from rare larger-effect variants that was more important in teosinte but depleted in maize following the domestication bottleneck. Purging associated with the maize domestication bottleneck may have selected against some large effect variants, but polygenic load is harder to purge and overall segregating mutational burden increased in maize compared to teosinte.

Evaluation of Reference Genes for Transcriptional Profiling in Two Cockroach Models.

The German cockroach, Blattella germanica, and the American cockroach, Periplaneta americana are the most common and synanthropic household pests of interest to public health. While they have increasingly served as model systems in hemimetabolous insects for studying many biological issues, there is still a lack of stable reference gene evaluation for reliable quantitative real-time PCR (qPCR) outputs and functional genomics. Here, we evaluated the expression variation of common insect reference genes, including the historically used actin, across various tissues and developmental stages, and also under experimental treatment conditions in these two species by using three individual algorithms (geNorm, BestKeeper, and NormFinder) and a comprehensive program (RefFinder). RPL32 in B. germanica and EF1α in P. americana showed the overall lowest variation among all examined samples. Based on the stability rankings by RefFinder, the optimal but varied reference genes under specific conditions were selected for qPCR normalization. In addition, the combination of RPL32 and EF1α was recommended for all the tested tissues and stages in B. germanica, whereas the combination of multiple reference genes was unfavorable in P. americana. This study provides a condition-specific resource of reference gene selection for accurate gene expression profiling and facilitating functional genomics in these two important cockroaches.

The Impact of Omega-3 Fatty Acids on the Evolution of Acinetobacter baumannii Drug Resistance.

The bacterial pathogen Acinetobacter baumannii has emerged as an urgent threat to health care systems. The prevalence of multidrug resistance in this critical human pathogen is closely associated with difficulties in its eradication from the hospital environment and its recalcitrance to treatment during infection. The development of resistance in A. baumannii is in part due to substantial plasticity of its genome, facilitating spontaneous genomic evolution. Many studies have investigated selective pressures imposed by antibiotics on genomic evolution, but the influence of high-abundance bioactive molecules at the host-pathogen interface on mutation and rates of evolution is poorly understood. Here, we studied the roles of host fatty acids in the gain in resistance to common antibiotics. We defined the impact of the polyunsaturated fatty acids arachidonic acid and docosahexaenoic acid on the development of resistance to erythromycin in A. baumannii strain AB5075_UW using a microevolutionary approach. We employed whole-genome sequencing and various phenotypic analyses to characterize microbe-lipid-antibiotic interactions. Cells exposed to erythromycin in the presence of the fatty acids displayed significantly lower rates of development of resistance to erythromycin and, importantly, tetracycline. Subsequent analyses defined diverse means by which host fatty acids influence the mutation rates. This work has highlighted the critical need to consider the roles of host fatty acids in A. baumannii physiology and antimicrobial resistance. Collectively, we have identified a novel means to curb the development of resistance in this critical human pathogen. IMPORTANCE The global distribution of multidrug resistance in A. baumannii has necessitated seeking not only alternative therapeutic approaches but also the means to limit the development of resistance in clinical settings. Highly abundant host bioactive compounds, such as polyunsaturated fatty acids, are readily acquired by A. baumannii during infection and have been illustrated to impact the bacterium's membrane composition and antibiotic resistance. In this work, we show that in vitro supplementation with host polyunsaturated fatty acids reduces the rate at which A. baumannii gains resistance to erythromycin and tetracycline. Furthermore, we discover that the impact on resistance development is closely associated with the primary antimicrobial efflux systems of A. baumannii, which represent one of the major drivers of clinical resistance. Overall, this study emphasizes the potential of host macromolecules in novel approaches to circumvent the difficulties of multidrug resistance during A. baumannii treatment, with fatty acid supplements such as fish oil providing safe and cost-effective ways to enhance host tolerance to bacterial infections.

Evolution of Complexity. Molecular Aspects of Preassembly.

An extension of neo-Darwinism, termed preassembly, states that genetic material required for many complex traits, such as echolocation, was present long before emergence of the traits. Assembly of genes and gene segments had occurred over protracted time-periods within large libraries of non-coding genes. Epigenetic factors ultimately promoted transfers from noncoding to coding genes, leading to abrupt formation of the trait via de novo genes. This preassembly model explains many observations that to this present day still puzzle biologists: formation of super-complexity in the absence of multiple fossil precursors, as with bat echolocation and flowering plants; major genetic and physical alterations occurring in just a few thousand years, as with housecat evolution; lack of precursors preceding lush periods of species expansion, as in the Cambrian explosion; and evolution of costly traits that exceed their need during evolutionary times, as with human intelligence. What follows in this paper is a mechanism that is not meant to supplant neo-Darwinism; instead, preassembly aims to supplement current ideas when complexity issues leave them struggling.

COVID-19 pandemic and "survival of the fittest".

Occurrence and recurrence of COVID-19 cases have been observed globally. The complex relationship of host-pathogen and the environment plays a vital role in understanding the widespread recurrence of the SARS-CoV-2 among humans. Though the pathobiology of the disease is not completely understood, it is well established that COVID-19 poses a greater threat to individuals with co-morbidities and a weakened immune system. The article deals with the notion of innate immunity, natural selection, and the survival of the fittest during the COVID-19 outbreak. The article also attempts to introduce the concept of "lifestyle and cultural immunity" that needs to be addressed and incorporated at an early stage of childhood to boost up the human immune system. The communication further discusses the role of vaccination and micro-organisms pre-existing in the environment which are required to enhance the immunity of an individual.

Preimplantation genetic testing for BRCA gene mutation carriers: a cost effectiveness analysis.

BACKGROUND: Gynecologic oncologists should be aware of the option of conception through IVF/PGT-M for families with high BRCA related morbidity or mortality. Our objective was to investigate the cost-effectiveness of preimplantation genetic testing for selection and transfer of BRCA negative embryo in BRCA mutation carriers compared to natural conception. METHODS: Cost-effectiveness of two strategies, conception through IVF/PGT-M and BRCA negative embryo transfer versus natural conception with a 50% chance of BRCA positive newborn for BRCA mutation carriers was compared using a Markovian process decision analysis model. Costs of the two strategies were compared using quality adjusted life years (QALYs'). All costs were discounted at 3%. Incremental cost effectiveness ratio (ICER) compared to willingness to pay threshold was used for cost-effectiveness analysis. RESULTS: IVF/ PGT-M is cost-effective with an ICER of 150,219 new Israeli Shekels, per QALY gained (equivalent to 44,480 USD), at a 3% discount rate. CONCLUSIONS: IVF/ PGT-M and BRCA negative embryo transfer compared to natural conception among BRCA positive parents is cost effective and may be offered for selected couples with high BRCA mutation related morbidity or mortality. Our results could impact decisions regarding conception among BRCA positive couples and health care providers.