Experimental Evaluation of a Direct Fitness Effect of the De Novo Evolved Mouse Gene Pldi.

De novo evolved genes emerge from random parts of noncoding sequences and have, therefore, no homologs from which a function could be inferred. While expression analysis and knockout experiments can provide insights into the function, they do not directly test whether the gene is beneficial for its carrier. Here, we have used a seminatural environment experiment to test the fitness of the previously identified de novo evolved mouse gene Pldi, which has been implicated to have a role in sperm differentiation. We used a knockout mouse strain for this gene and competed it against its parental wildtype strain for several generations of free reproduction. We found that the knockout (ko) allele frequency decreased consistently across three replicates of the experiment. Using an approximate Bayesian computation framework that simulated the data under a demographic scenario mimicking the experiment's demography, we could estimate a selection coefficient ranging between 0.21 and 0.61 for the wildtype allele compared to the ko allele in males, under various models. This implies a relatively strong selective advantage, which would fix the new gene in less than hundred generations after its emergence.

Living fast, dying young: Anthropogenic habitat modification influences the fitness and life history traits of a cooperative breeder.

Anthropogenic habitat modification can indirectly effect reproduction and survival in social species by changing the group structure and social interactions. We assessed the impact of habitat modification on the fitness and life history traits of a cooperative breeder, the Arabian babbler (Argya squamiceps). We collected spatial, reproductive and social data on 572 individuals belonging to 21 social groups over 6 years and combined it with remote sensing to characterize group territories in an arid landscape. In modified resource-rich habitats, groups bred more and had greater productivity, but individuals lived shorter lives than in natural habitats. Habitat modification favoured a faster pace-of-life with lower dispersal and dominance acquisition ages, which might be driven by higher mortality providing opportunities for the dominant breeding positions. Thus, habitat modification might indirectly impact fitness through changes in social structures. This study shows that trade-offs in novel anthropogenic opportunities might offset survival costs by increased productivity.

Genotypic-phenotypic landscape computation based on first principle and deep learning.

The relationship between genotype and fitness is fundamental to evolution, but quantitatively mapping genotypes to fitness has remained challenging. We propose the Phenotypic-Embedding theorem (P-E theorem) that bridges genotype-phenotype through an encoder-decoder deep learning framework. Inspired by this, we proposed a more general first principle for correlating genotype-phenotype, and the P-E theorem provides a computable basis for the application of first principle. As an application example of the P-E theorem, we developed the Co-attention based Transformer model to bridge Genotype and Fitness model, a Transformer-based pre-train foundation model with downstream supervised fine-tuning that can accurately simulate the neutral evolution of viruses and predict immune escape mutations. Accordingly, following the calculation path of the P-E theorem, we accurately obtained the basic reproduction number (${R}_0$) of SARS-CoV-2 from first principles, quantitatively linked immune escape to viral fitness and plotted the genotype-fitness landscape. The theoretical system we established provides a general and interpretable method to construct genotype-phenotype landscapes, providing a new paradigm for studying theoretical and computational biology.

Variation and plasticity in life-history traits and fitness of wild Arabidopsis thaliana populations are not related to their genotypic and ecological diversity.

BACKGROUND: Despite its implications for population dynamics and evolution, the relationship between genetic and phenotypic variation in wild populations remains unclear. Here, we estimated variation and plasticity in life-history traits and fitness of the annual plant Arabidopsis thaliana in two common garden experiments that differed in environmental conditions. We used up to 306 maternal inbred lines from six Iberian populations characterized by low and high genotypic (based on whole-genome sequences) and ecological (vegetation type) diversity. RESULTS: Low and high genotypic and ecological diversity was found in edge and core Iberian environments, respectively. Given that selection is expected to be stronger in edge environments and that ecological diversity may enhance both phenotypic variation and plasticity, we expected genotypic diversity to be positively associated with phenotypic variation and plasticity. However, maternal lines, irrespective of the genotypic and ecological diversity of their population of origin, exhibited a substantial amount of phenotypic variation and plasticity for all traits. Furthermore, all populations harbored maternal lines with canalization (robustness) or sensitivity in response to harsher environmental conditions in one of the two experiments. CONCLUSIONS: Overall, we conclude that the environmental attributes of each population probably determine their genotypic diversity, but all populations maintain substantial phenotypic variation and plasticity for all traits, which represents an asset to endure in changing environments.

Experimental evolution for the recovery of growth loss due to genome reduction.

As the genome encodes the information crucial for cell growth, a sizeable genomic deficiency often causes a significant decrease in growth fitness. Whether and how the decreased growth fitness caused by genome reduction could be compensated by evolution was investigated here. Experimental evolution with an Escherichia coli strain carrying a reduced genome was conducted in multiple lineages for approximately 1000 generations. The growth rate, which largely declined due to genome reduction, was considerably recovered, associated with the improved carrying capacity. Genome mutations accumulated during evolution were significantly varied across the evolutionary lineages and were randomly localized on the reduced genome. Transcriptome reorganization showed a common evolutionary direction and conserved the chromosomal periodicity, regardless of highly diversified gene categories, regulons, and pathways enriched in the differentially expressed genes. Genome mutations and transcriptome reorganization caused by evolution, which were found to be dissimilar to those caused by genome reduction, must have followed divergent mechanisms in individual evolutionary lineages. Gene network reconstruction successfully identified three gene modules functionally differentiated, which were responsible for the evolutionary changes of the reduced genome in growth fitness, genome mutation, and gene expression, respectively. The diversity in evolutionary approaches improved the growth fitness associated with the homeostatic transcriptome architecture as if the evolutionary compensation for genome reduction was like all roads leading to Rome.

Physiological Measurements and Transcriptomics Reveal the Fitness Costs of Monochamus saltuarius to Bursaphelenchus xylophilus.

The pine wood nematode (PWN) uses several Monochamus species as vehicles, through a temporary hitchhiking process known as phoresy, enabling it to access new host plant resources. Monochamus saltuarius acts as a new and major vector of the PWN in Northeastern China, showing lower PWN carrying capacity and a shorter transmission cycle compared to established vectors. The apparently altered symbiotic relationship offers an interesting area for researching the costs and adaptions involved in nematode-beetle, a specialized phoresy. We analyzed the response and fitness costs of M. saltuarius through physiological measurements and transcriptomics. The PWN exerted adverse repercussions on the growth and development of M. saltuarius. The PWN accelerated larval development into pupae, while beetle adults carrying the PWN exhibited an elevated abnormality rate and mortality, and reduced starvation resistance. During the pupal stage, the expression of growth-related genes, including ecdysone-inducible genes (E74EA), cuticle proteins, and chitin genes (CHTs), markedly increased. Meanwhile, the induced immune response, mainly by the IMD and Toll signaling pathways, could be a contributing factor to adult abnormality and mortality. Adult gonads and trachea exhibited enrichment in pathways related to fatty acid elongation, biosynthesis, and metabolism. FASN, ELOVL, and SCD possibly contributed to resistance against PWN. Our research indicated that phoretic interactions between vector beetles and PWN vary throughout the vector's lifespan, particularly before and after entry into the trachea. This study highlighted the fitness costs of immunity and metabolism on the vector beetle, indicating the adaptation mechanisms and evolutionary trade-offs to PWN.

Single-cell analysis reveals context-dependent, cell-level selection of mtDNA.

Heteroplasmy occurs when wild-type and mutant mitochondrial DNA (mtDNA) molecules co-exist in single cells1. Heteroplasmy levels change dynamically in development, disease and ageing2,3, but it is unclear whether these shifts are caused by selection or drift, and whether they occur at the level of cells or intracellularly. Here we investigate heteroplasmy dynamics in dividing cells by combining precise mtDNA base editing (DdCBE)4 with a new method, SCI-LITE (single-cell combinatorial indexing leveraged to interrogate targeted expression), which tracks single-cell heteroplasmy with ultra-high throughput. We engineered cells to have synonymous or nonsynonymous complex I mtDNA mutations and found that cell populations in standard culture conditions purge nonsynonymous mtDNA variants, whereas synonymous variants are maintained. This suggests that selection dominates over simple drift in shaping population heteroplasmy. We simultaneously tracked single-cell mtDNA heteroplasmy and ancestry, and found that, although the population heteroplasmy shifts, the heteroplasmy of individual cell lineages remains stable, arguing that selection acts at the level of cell fitness in dividing cells. Using these insights, we show that we can force cells to accumulate high levels of truncating complex I mtDNA heteroplasmy by placing them in environments where loss of biochemical complex I activity has been reported to benefit cell fitness. We conclude that in dividing cells, a given nonsynonymous mtDNA heteroplasmy can be harmful, neutral or even beneficial to cell fitness, but that the 'sign' of the effect is wholly dependent on the environment.

Fitness surfaces and local thermal adaptation in Drosophila along a latitudinal gradient.

Local adaptation is commonly cited to explain species distribution, but how fitness varies along continuous geographical gradients is not well understood. Here, we combine thermal biology and life-history theory to demonstrate that Drosophila populations along a 2500 km latitudinal cline are adapted to local conditions. We measured how heat tolerance and viability rate across eight populations varied with temperature in the laboratory and then simulated their expected cumulative Darwinian fitness employing high-resolution temperature data from their eight collection sites. Simulations indicate a trade-off between annual survival and cumulative viability, as both mortality and the recruitment of new flies are predicted to increase in warmer regions. Importantly, populations are locally adapted and exhibit the optimal combination of both traits to maximize fitness where they live. In conclusion, our method is able to reconstruct fitness surfaces employing empirical life-history estimates and reconstructs peaks representing locally adapted populations, allowing us to study geographic adaptation in silico.

Higher-order species interactions cause time-dependent niche and fitness differences: Experimental evidence in plant-feeding arthropods.

Species interact in different ways, including competition, facilitation and predation. These interactions can be non-linear or higher order and may depend on time or species densities. Although these higher-order interactions are virtually ubiquitous, they remain poorly understood, as they are challenging both theoretically and empirically. We propose to adapt niche and fitness differences from modern coexistence theory and apply them to species interactions over time. As such, they may not merely inform about coexistence, but provide a deeper understanding of how species interactions change. Here, we investigated how the exploitation of a biotic resource (plant) by phytophagous arthropods affects their interactions. We performed monoculture and competition experiments to fit a generalized additive mixed model to the empirical data, which allowed us to calculate niche and fitness differences. We found that species switch between different types of interactions over time, including intra- and interspecific facilitation, and strong and weak competition.

U-CAN-seq: A Universal Competition Assay by Nanopore Sequencing.

RNA viruses quickly evolve subtle genotypic changes that can have major impacts on viral fitness and host range, with potential consequences for human health. It is therefore important to understand the evolutionary fitness of novel viral variants relative to well-studied genotypes of epidemic viruses. Competition assays are an effective and rigorous system with which to assess the relative fitness of viral genotypes. However, it is challenging to quickly and cheaply distinguish and quantify fitness differences between very similar viral genotypes. Here, we describe a protocol for using reverse transcription PCR in combination with commercial nanopore sequencing services to perform competition assays on untagged RNA viruses. Our assay, called the Universal Competition Assay by Nanopore Sequencing (U-CAN-seq), is relatively cheap and highly sensitive. We used a well-studied N24A mutation in the chikungunya virus (CHIKV) nsp3 gene to confirm that we could detect a competitive disadvantage using U-CAN-seq. We also used this approach to show that mutations to the CHIKV 5' conserved sequence element that disrupt sequence but not structure did not affect the fitness of CHIKV. However, similar mutations to an adjacent CHIKV stem loop (SL3) did cause a fitness disadvantage compared to wild-type CHIKV, suggesting that structure-independent, primary sequence determinants in this loop play an important role in CHIKV biology. Our novel findings illustrate the utility of the U-CAN-seq competition assay.

Pre- and postnatal noise directly impairs avian development, with fitness consequences.

Noise pollution is expanding at an unprecedented rate and is increasingly associated with impaired reproduction and development across taxa. However, whether noise sound waves are intrinsically harmful for developing young-or merely disturb parents-and the fitness consequences of early exposure remain unknown. Here, by only manipulating the offspring, we show that sole exposure to noise in early life in zebra finches has fitness consequences and causes embryonic death during exposure. Exposure to pre- and postnatal traffic noise cumulatively impaired nestling growth and physiology and aggravated telomere shortening across life stages until adulthood. Consistent with a long-term somatic impact, early life noise exposure, especially prenatally, decreased individual offspring production throughout adulthood. Our findings suggest that the effects of noise pollution are more pervasive than previously realized.

SARS-CoV-2 variant of concern fitness and adaptation in primary human airway epithelia.

The severe acute respiratory syndrome coronavirus 2 pandemic is characterized by the emergence of novel variants of concern (VOCs) that replace ancestral strains. Here, we dissect the complex selective pressures by evaluating variant fitness and adaptation in human respiratory tissues. We evaluate viral properties and host responses to reconstruct forces behind D614G through Omicron (BA.1) emergence. We observe differential replication in airway epithelia, differences in cellular tropism, and virus-induced cytotoxicity. D614G accumulates the most mutations after infection, supporting zoonosis and adaptation to the human airway. We perform head-to-head competitions and observe the highest fitness for Gamma and Delta. Under these conditions, RNA recombination favors variants encoding the B.1.617.1 lineage 3' end. Based on viral growth kinetics, Alpha, Gamma, and Delta exhibit increased fitness compared to D614G. In contrast, the global success of Omicron likely derives from increased transmission and antigenic variation. Our data provide molecular evidence to support epidemiological observations of VOC emergence.

Antibiotic resistance begets more resistance: chromosomal resistance mutations mitigate fitness costs conferred by multi-resistant clinical plasmids.

Plasmids are the primary vectors of horizontal transfer of antibiotic resistance genes among bacteria. Previous studies have shown that the spread and maintenance of plasmids among bacterial populations depend on the genetic makeup of both the plasmid and the host bacterium. Antibiotic resistance can also be acquired through mutations in the bacterial chromosome, which not only confer resistance but also result in changes in bacterial physiology and typically a reduction in fitness. However, it is unclear whether chromosomal resistance mutations affect the interaction between plasmids and the host bacteria. To address this question, we introduced 13 clinical plasmids into a susceptible Escherichia coli strain and three different congenic mutants that were resistant to nitrofurantoin (ΔnfsAB), ciprofloxacin (gyrA, S83L), and streptomycin (rpsL, K42N) and determined how the plasmids affected the exponential growth rates of the host in glucose minimal media. We find that though plasmids confer costs on the susceptible strains, those costs are fully mitigated in the three resistant mutants. In several cases, this results in a competitive advantage of the resistant strains over the susceptible strain when both carry the same plasmid and are grown in the absence of antibiotics. Our results suggest that bacteria carrying chromosomal mutations for antibiotic resistance could be a better reservoir for resistance plasmids, thereby driving the evolution of multi-drug resistance.IMPORTANCEPlasmids have led to the rampant spread of antibiotic resistance genes globally. Plasmids often carry antibiotic resistance genes and other genes needed for its maintenance and spread, which typically confer a fitness cost on the host cell observed as a reduced growth rate. Resistance is also acquired via chromosomal mutations, and similar to plasmids they also reduce bacterial fitness. However, we do not know whether resistance mutations affect the bacterial ability to carry plasmids. Here, we introduced 13 multi-resistant clinical plasmids into a susceptible and three different resistant E. coli strains and found that most of these plasmids do confer fitness cost on susceptible cells, but these costs disappear in the resistant strains which often lead to fitness advantage for the resistant strains in the absence of antibiotic selection. Our results imply that already resistant bacteria are a more favorable reservoir for multi-resistant plasmids, promoting the ascendance of multi-resistant bacteria.

High temperature delays and low temperature accelerates evolution of a new protein phenotype.

Since the origin of life, temperatures on earth have fluctuated both on short and long time scales. How such changes affect the rate at which Darwinian evolution can bring forth new phenotypes remains unclear. On the one hand, high temperature may accelerate phenotypic evolution because it accelerates most biological processes. On the other hand, it may slow phenotypic evolution, because proteins are usually less stable at high temperatures and therefore less evolvable. Here, to test these hypotheses experimentally, we evolved a green fluorescent protein in E. coli towards the new phenotype of yellow fluorescence at different temperatures. Yellow fluorescence evolved most slowly at high temperature and most rapidly at low temperature, in contradiction to the first hypothesis. Using high-throughput population sequencing, protein engineering, and biochemical assays, we determined that this is due to the protein-destabilizing effect of neofunctionalizing mutations. Destabilization is highly detrimental at high temperature, where neofunctionalizing mutations cannot be tolerated. Their detrimental effects can be mitigated through excess stability at low temperature, leading to accelerated adaptive evolution. By modifying protein folding stability, temperature alters the accessibility of mutational paths towards high-fitness genotypes. Our observations have broad implications for our understanding of how temperature changes affect evolutionary adaptations and innovations.

Genetic architecture of inbreeding depression may explain its persistence in a population of wild red deer.

Inbreeding depression is of major concern in declining populations, but relatively little is known about its genetic architecture in wild populations, such as the degree to which it is composed of large or small effect loci and their distribution throughout the genome. Here, we combine fitness and genomic data from a wild population of red deer to investigate the genomic distribution of inbreeding effects. Based on the runs of homozygosity (ROH)-based inbreeding coefficient, FROH, we use chromosome-specific inbreeding coefficients (FROHChr) to explore whether the effect of inbreeding varies between chromosomes. Under the assumption that within an individual the probability of being identical-by-descent is equal across all chromosomes, we used a multi-membership model to estimate the deviation of FROHChr from the average inbreeding effect. This novel approach ensures effect sizes are not overestimated whilst maximising the power of our available dataset of >3000 individuals genotyped on >35,000 autosomal SNPs. We find that most chromosomes confer a minor reduction in fitness-related traits, which when these effects are summed, results in the observed inbreeding depression in birth weight, survival and lifetime breeding success. However, no chromosomes had a significant detrimental effect compared to the overall effect of inbreeding, indicating no major effect loci. We conclude that in this population, inbreeding depression is likely the result of multiple mildly or moderately deleterious mutations spread across all chromosomes, which are difficult to detect with statistical confidence. Such mutations will be inefficiently purged, which may explain the persistence of inbreeding depression in this population.

Metabolic disruptions and impaired reproductive fitness in wild-caught freshwater turtles (Emydura macquarii macquarii) exposed to elevated per- and polyfluoroalkyl substances (PFAS).

Per- and poly-fluoroalkyl substances (PFAS) pose a threat to organisms and ecosystems due to their persistent nature. Ecotoxicology endpoints used in regulatory guidelines may not reflect multiple, low-level but persistent stressors. This study examines the biological effects of PFAS on Eastern short-necked turtles in Queensland, Australia. In this study, blood samples were collected and analysed for PFAS, hormone levels, and functional omics endpoints. High levels of PFAS were found in turtles at the impacted site, with PFOS being the dominant constituent. The PFAS profiles of males and females differed, with males having higher PFAS concentrations. Hormone concentrations differed between impacted and reference sites in male turtles, with elevated testosterone and corticosterone indicative of stress. Further, energy utilisation, nucleotide synthesis, nitrogen metabolism, and amino acid synthesis were altered in both male and female turtles from PFAS-impacted sites. Both sexes show similar metabolic responses to environmental stressors from the PFAS-contaminated site, which may adversely affect their reproductive fitness. Purine metabolism, caffeine metabolism, and ferroptosis pathway changes in turtles can cause gout, cell death, and overall health problems. Further, the study showed that prolonged exposure to elevated PFAS levels in the wild could compromise turtle reproductive fitness by disrupting reproductive steroids and metabolic pathways.

Compensatory evolution in NusG improves fitness of drug-resistant M. tuberculosis.

Drug-resistant bacteria are emerging as a global threat, despite frequently being less fit than their drug-susceptible ancestors1-8. Here we sought to define the mechanisms that drive or buffer the fitness cost of rifampicin resistance (RifR) in the bacterial pathogen Mycobacterium tuberculosis (Mtb). Rifampicin inhibits RNA polymerase (RNAP) and is a cornerstone of modern short-course tuberculosis therapy9,10. However, RifR Mtb accounts for one-quarter of all deaths due to drug-resistant bacteria11,12. We took a comparative functional genomics approach to define processes that are differentially vulnerable to CRISPR interference (CRISPRi) inhibition in RifR Mtb. Among other hits, we found that the universally conserved transcription factor NusG is crucial for the fitness of RifR Mtb. In contrast to its role in Escherichia coli, Mtb NusG has an essential RNAP pro-pausing function mediated by distinct contacts with RNAP and the DNA13. We find this pro-pausing NusG-RNAP interface to be under positive selection in clinical RifR Mtb isolates. Mutations in the NusG-RNAP interface reduce pro-pausing activity and increase fitness of RifR Mtb. Collectively, these results define excessive RNAP pausing as a molecular mechanism that drives the fitness cost of RifR in Mtb, identify a new mechanism of compensation to overcome this cost, suggest rational approaches to exacerbate the fitness cost, and, more broadly, could inform new therapeutic approaches to develop drug combinations to slow the evolution of RifR in Mtb.

A quantitative genetic model of background selection in humans.

Across the human genome, there are large-scale fluctuations in genetic diversity caused by the indirect effects of selection. This "linked selection signal" reflects the impact of selection according to the physical placement of functional regions and recombination rates along chromosomes. Previous work has shown that purifying selection acting against the steady influx of new deleterious mutations at functional portions of the genome shapes patterns of genomic variation. To date, statistical efforts to estimate purifying selection parameters from linked selection models have relied on classic Background Selection theory, which is only applicable when new mutations are so deleterious that they cannot fix in the population. Here, we develop a statistical method based on a quantitative genetics view of linked selection, that models how polygenic additive fitness variance distributed along the genome increases the rate of stochastic allele frequency change. By jointly predicting the equilibrium fitness variance and substitution rate due to both strong and weakly deleterious mutations, we estimate the distribution of fitness effects (DFE) and mutation rate across three geographically distinct human samples. While our model can accommodate weaker selection, we find evidence of strong selection operating similarly across all human samples. Although our quantitative genetic model of linked selection fits better than previous models, substitution rates of the most constrained sites disagree with observed divergence levels. We find that a model incorporating selective interference better predicts observed divergence in conserved regions, but overall our results suggest uncertainty remains about the processes generating fitness variation in humans.

A fitness landscape instability governs the morphological diversity of tip-growing cells.

Cellular morphology affects many aspects of cellular and organismal physiology. This makes it challenging to dissect the evolutionary basis for specific morphologies since various cellular functions may exert competing selective pressures on this trait, and the influence of these pressures will depend on the specific mechanisms of morphogenesis. In this light, we combined experiment and theory to investigate the complex basis for morphological diversity among tip-growing cells from across the tree of life. We discovered that an instability in the widespread mechanism of "inflationary" tip growth leads directly to a bifurcation in the common fitness landscape of tip-growing cells, which imposes a strict global constraint on their morphologies. This result rationalizes the morphology of an enormous diversity of important fungal, plant, protistan, and bacterial systems. More broadly, our study elucidates the principle that strong evolutionary constraints on complex traits, like biological form, may emerge from emergent instabilities within developmental systems.

Dominance reversals: the resolution of genetic conflict and maintenance of genetic variation.

Beneficial reversals of dominance reduce the costs of genetic trade-offs and can enable selection to maintain genetic variation for fitness. Beneficial dominance reversals are characterized by the beneficial allele for a given context (e.g. habitat, developmental stage, trait or sex) being dominant in that context but recessive where deleterious. This context dependence at least partially mitigates the fitness consequence of heterozygotes carrying one non-beneficial allele for their context and can result in balancing selection that maintains alternative alleles. Dominance reversals are theoretically plausible and are supported by mounting empirical evidence. Here, we highlight the importance of beneficial dominance reversals as a mechanism for the mitigation of genetic conflict and review the theory and empirical evidence for them. We identify some areas in need of further research and development and outline three methods that could facilitate the identification of antagonistic genetic variation (dominance ordination, allele-specific expression and allele-specific ATAC-Seq (assay for transposase-accessible chromatin with sequencing)). There is ample scope for the development of new empirical methods as well as reanalysis of existing data through the lens of dominance reversals. A greater focus on this topic will expand our understanding of the mechanisms that resolve genetic conflict and whether they maintain genetic variation.