Search

Epistasis is caused by genetic interactions among mutations that affect fitness. To characterize properties and potential mechanisms of epistasis, we engineered eight double mutants that combined mutations from the rho and rpoB genes of Escherichia coli. The two genes encode essential functions for transcription, and the mutations in each gene were chosen because they were beneficial for adaptation to thermal stress (42.2 °C). The double mutants exhibited patterns of fitness epistasis that included diminishing returns epistasis at 42.2 °C, stronger diminishing returns between mutations with larger beneficial effects and both negative and positive (sign) epistasis across environments (20.0 °C and 37.0 °C). By assessing gene expression between single and double mutants, we detected hundreds of genes with gene expression epistasis. Previous work postulated that highly connected hub genes in coexpression networks have low epistasis, but we found the opposite: hub genes had high epistasis values in both coexpression and protein-protein interaction networks. We hypothesized that elevated epistasis in hub genes reflected that they were enriched for targets of Rho termination but that was not the case. Altogether, gene expression and coexpression analyses revealed that thermal adaptation occurred in modules, through modulation of ribonucleotide biosynthetic processes and ribosome assembly, the attenuation of expression in genes related to heat shock and stress responses, and with an overall trend toward restoring gene expression toward the unstressed state.

The relationship between genotype and phenotype, or the fitness landscape, is the foundation of genetic engineering and evolution. However, mapping fitness landscapes poses a major technical challenge due to the amount of quantifiable data that is required. Catalytic RNA is a special topic in the study of fitness landscapes due to its relatively small sequence space combined with its importance in synthetic biology. The combination of in vitro selection and high-throughput sequencing has recently provided empirical maps of both complete and local RNA fitness landscapes, but the astronomical size of sequence space limits purely experimental investigations. Next steps are likely to involve data-driven interpolation and extrapolation over sequence space using various machine learning techniques. We discuss recent progress in understanding RNA fitness landscapes, particularly with respect to protocells and machine representations of RNA. The confluence of technical advances may significantly impact synthetic biology in the near future.

Selection causes local adaptation across populations within species and simultaneously divergence between species. However, it is unclear if either the force of or the response to selection is similar across these scales. We show that natural selection drives divergence between closely related species in a pattern that is distinct from local adaptation within species. We use reciprocal transplant experiments across three species of Phlox wildflowers to characterize widespread adaptive divergence. Using provenance trials, we also find strong local adaptation between populations within a species. Comparing divergence and selection between these two scales of diversity we discover that one suite of traits predicts fitness differences between species and that an independent suite of traits predicts fitness variation within species. Selection drives divergence between species, contributing to speciation, while simultaneously favoring extensive diversity that is maintained across populations within a species. Our work demonstrates how the selection landscape is complex and multidimensional.

Many factors contribute to the ability of a microbial species to persist when encountering complexly contaminated environments including time of exposure, the nature and concentration of contaminants, availability of nutritional resources, and possession of a combination of appropriate molecular mechanisms needed for survival. Herein we sought to identify genes that are most important for survival of Gram-negative Enterobacteriaceae in contaminated groundwater environments containing high concentrations of nitrate and metals using the metal-tolerant Oak Ridge Reservation (ORR) isolate, Pantoea sp. MT58 (MT58). Survival fitness experiments in which a randomly barcoded transposon insertion (RB-TnSeq) library of MT58 was exposed directly to contaminated ORR groundwater samples from across a nitrate and mixed metal contamination plume were used to identify genes important for survival with increasing exposure times and concentrations of contaminants, and availability of a carbon source. Genes involved in controlling and using carbon, encoding transcriptional regulators, and related to Gram-negative outer membrane processes were among those found to be important for survival in contaminated ORR groundwater. A comparative genomics analysis of 75 Pantoea genus strains allowed us to further separate the survival determinants into core and non-core genes in the Pantoea pangenome, revealing insights into the survival of subsurface microorganisms during contaminant plume intrusion.

A major challenge in evolutionary biology is explaining how populations navigate rugged fitness landscapes without getting trapped on local optima. One idea illustrated by adaptive dynamics theory is that as populations adapt, their newly enhanced capacities to exploit resources alter fitness payoffs and restructure the landscape in ways that promote speciation by opening new adaptive pathways. While there have been indirect tests of this theory, to our knowledge none have measured how fitness landscapes deform during adaptation, or test whether these shifts promote diversification. Here, we achieve this by studying bacteriophage [Formula: see text], a virus that readily speciates into co-existing receptor specialists under controlled laboratory conditions. We use a high-throughput gene editing-phenotyping technology to measure [Formula: see text]'s fitness landscape in the presence of different evolved-[Formula: see text] competitors and find that the fitness effects of individual mutations, and their epistatic interactions, depend on the competitor. Using these empirical data, we simulate [Formula: see text]'s evolution on an unchanging landscape and one that recapitulates how the landscape deforms during evolution. [Formula: see text] heterogeneity only evolves in the shifting landscape regime. This study provides a test of adaptive dynamics, and, more broadly, shows how fitness landscapes dynamically change during adaptation, potentiating phenomena like speciation by opening new adaptive pathways.

The identification of catalytic RNAs is typically achieved through primarily experimental means. However, only a small fraction of sequence space can be analyzed even with high-throughput techniques. Methods to extrapolate from a limited data set to predict additional ribozyme sequences, particularly in a human-interpretable fashion, could be useful both for designing new functional RNAs and for generating greater understanding about a ribozyme fitness landscape. Using information theory, we express the effects of epistasis (i.e., deviations from additivity) on a ribozyme. This representation was incorporated into a simple model of the epistatic fitness landscape, which identified potentially exploitable combinations of mutations. We used this model to theoretically predict mutants of high activity for a self-aminoacylating ribozyme, identifying potentially active triple and quadruple mutants beyond the experimental data set of single and double mutants. The predictions were validated experimentally, with nine out of nine sequences being accurately predicted to have high activity. This set of sequences included mutants that form a previously unknown evolutionary bridge between two ribozyme families that share a common motif. Individual steps in the method could be examined, understood, and guided by a human, combining interpretability and performance in a simple model to predict ribozyme sequences by extrapolation.

Social behaviours are typically modelled using neighbour-modulated fitness, which focuses on individuals having their fitness altered by neighbours. However, these models are either interpreted using inclusive fitness, which focuses on individuals altering the fitness of neighbours, or not interpreted at all. This disconnect leads to interpretational mistakes and obscures the adaptive significance of behaviour. We bridge this gap by presenting a systematic methodology for constructing inclusive-fitness models. We find a behaviours inclusive-fitness effect by summing primary and secondary deviations in reproductive value. Primary deviations are the immediate result of a social interaction; for example, the cost and benefit of an altruistic act. Secondary deviations are compensatory effects that arise because the total reproductive value of the population is fixed; for example, the increased competition that follows an altruistic act. Compared to neighbour-modulated fitness methodologies, our approach is often simpler and reveals the models inclusive-fitness narrative clearly. We implement our methodology first in a homogeneous population, with supplementary examples of help under synergy, help in a viscous population and Creels paradox. We then implement our methodology in a class-structured population, where the advantages of our approach are most evident, with supplementary examples of altruism between age classes, and sex-ratio evolution.

The fitness effects of all possible mutations available to an organism largely shape the dynamics of evolutionary adaptation. Yet, whether and how this adaptive landscape changes over evolutionary times, especially upon ecological diversification and changes in community composition, remains poorly understood. We sought to fill this gap by analyzing a stable community of two closely related ecotypes ("L" and "S") shortly after they emerged within the E. coli Long-Term Evolution Experiment (LTEE). We engineered genome-wide barcoded transposon libraries to measure the invasion fitness effects of all possible gene knockouts in the coexisting strains as well as their ancestor, for many different, ecologically relevant conditions. We find consistent statistical patterns of fitness effect variation across both genetic background and community composition, despite the idiosyncratic behavior of individual knockouts. Additionally, fitness effects are correlated with evolutionary outcomes for a number of conditions, possibly revealing shifting patterns of adaptation. Together, our results reveal how ecological and epistatic effects combine to shape the adaptive landscape in a nascent ecological community.

Widespread insecticide resistance in African malaria vectors raises concerns over the potential to compromise malaria vector control interventions. Understanding the evolution of resistance mechanisms, and whether the selective disadvantages are large enough to be useful in resistance management or designing suitable control strategies is crucial. This study assessed whether insecticide resistance to pyrethroids has an effect on the gonotrophic cycle and reproductive potential of malaria vector Anopheles gambiae. Comparative tests were performed with pyrethroid-resistant and susceptible colonies of Anopheles gambiae colonized from the same geographical area, and the reference Kisumu strain was used as a control. Adult females aged 3 days old were given a blood meal and kept separately for individual egg-laying. The number of days taken to lay eggs post-blood-feeding was recorded to determine the length of the gonotrophic cycle. To measure adult longevity and reproduction potential, newly emerged males and females of equal numbers were aspirated into a cage and females allowed to blood feed daily. The number of eggs laid and the surviving mosquitoes were recorded daily to determine fecundity, net reproduction rate, intrinsic growth rate and adult longevity. Overall, the resistant females had a significantly longer (1.8 days) gonotrophic cycle than susceptible females (F2, 13 = 9. 836, P

Direct observation of evolution in response to natural environmental change can resolve fundamental questions about adaptation, including its pace, temporal dynamics, and underlying phenotypic and genomic architecture. We tracked the evolution of fitness-associated phenotypes and allele frequencies genome-wide in 10 replicate field populations of Drosophila melanogaster over 10 generations from summer to late fall. Adaptation was evident over each sampling interval (one to four generations), with exceptionally rapid phenotypic adaptation and large allele frequency shifts at many independent loci. The direction and basis of the adaptive response shifted repeatedly over time, consistent with the action of strong and rapidly fluctuating selection. Overall, we found clear phenotypic and genomic evidence of adaptive tracking occurring contemporaneously with environmental change, thus demonstrating the temporally dynamic nature of adaptation.

Estimating the complex relationship between fitness and genotype or phenotype (i.e. the adaptive landscape) is one of the central goals of evolutionary biology. However, adaptive walks connecting genotypes to organismal fitness, speciation, and novel ecological niches are still poorly understood and processes for surmounting fitness valleys remain controversial. One outstanding system for addressing these connections is a recent adaptive radiation of ecologically and morphologically novel pupfishes (a generalist, molluscivore, and scale-eater) endemic to San Salvador Island, Bahamas. We leveraged whole-genome sequencing of 139 hybrids from two independent field fitness experiments to identify the genomic basis of fitness, estimate genotypic fitness networks, and measure the accessibility of adaptive walks on the fitness landscape. We identified 132 single nucleotide polymorphisms (SNPs) that were significantly associated with fitness in field enclosures. Six out of the 13 regions most strongly associated with fitness contained differentially expressed genes and fixed SNPs between trophic specialists; one gene (mettl21e) was also misexpressed in lab-reared hybrids, suggesting a potential intrinsic genetic incompatibility. We then constructed genotypic fitness networks from adaptive alleles and show that scale-eating specialists are the most isolated of the three species on these networks. Intriguingly, introgressed and de novo variants reduced fitness landscape ruggedness as compared to standing variation, increasing the accessibility of genotypic fitness paths from generalist to specialists. Our results suggest that adaptive introgression and de novo mutations alter the shape of the fitness landscape, providing key connections in adaptive walks circumventing fitness valleys and triggering the evolution of novelty during adaptive radiation.

To date, only low pathogenic (LP) H5 and H7 avian influenza viruses (AIV) have been observed to naturally shift to a highly pathogenic (HP) phenotype after mutation of the monobasic hemagglutinin (HA) cleavage site (HACS) to polybasic motifs. The LPAIV monobasic HACS is activated by tissue-restricted trypsin-like enzymes, while the HPAIV polybasic HACS is activated by ubiquitous furin-like enzymes. However, glycosylation near the HACS can affect proteolytic activation and reduced virulence of some HPAIV in chickens. In 2012, a unique H4N2 virus with a polybasic HACS was isolated from quails but was LP in chickens. Whether glycosylation sites (GS) near the HACS hinder the evolution of HPAIV H4N2 remains unclear. Here, we analyzed the prevalence of potential GS in the N-terminus of HA1, 2NYT4 and 18NGT20, in all AIV sequences and studied their impact on H4N2 virus fitness. Although the two motifs are conserved, some non-H5/H7 subtypes lack one or both GS. Both sites were glycosylated in this H4N2 virus. Deglycosylation increased trypsin-independent replication in cell culture, cell-to-cell spread and syncytium formation at low-acidic pH, but negatively affected the thermostability and receptor-binding affinity. Alteration of 2NYT4 with or without 18NGT20 enabled systemic spread of the virus to different organs including the brain of chicken embryos. However, all intranasally inoculated chickens did not show clinical signs. Together, although the conserved GS near the HACS are important for HA stability and receptor binding, deglycosylation increased the H4N2 HA-activation, replication and tissue tropism suggesting a potential role for virus adaptation in poultry.

BackgroundAnti-malarial drug resistance may be limited by decreased fitness in resistant parasites. Important contributors to resistance are mutations in the Plasmodium falciparum putative drug transporter PfMDR1.MethodsImpacts on in vitro fitness of two common PfMDR1 polymorphisms, N86Y, which is associated with sensitivity to multiple drugs, and Y184F, which has no clear impact on drug sensitivity, were evaluated to study associations between resistance mediators and parasite fitness, measured as relative growth in competitive culture experiments. NF10 P. falciparum lines engineered to represent all PfMDR1 N86Y and Y184F haplotypes were co-cultured for 40 days, and the genetic make-up of the cultures was characterized every 4 days by pyrosequencing. The impacts of culture with anti-malarials on the growth of different haplotypes were also assessed. Lastly, the engineering of P. falciparum containing another common polymorphism, PfMDR1 D1246Y, was attempted.ResultsCo-culture results were as follows. With wild type (WT) Y184 fixed (N86/Y184 vs. 86Y/Y184), parasites WT and mutant at 86 were at equilibrium. With mutant 184 F fixed (N86/184F vs. 86Y/184F), mutants at 86 overgrew WT. With WT N86 fixed (N86/Y184 vs. N86/184F), WT at 184 overgrew mutants. With mutant 86Y fixed (86Y/Y184 vs. 86Y/184F), WT and mutant at 86 were at equilibrium. Parasites with the double WT were in equilibrium with the double mutant, but 86Y/Y184 overgrew N86/184F. Overall, WT N86/mutant 184F parasites were less fit than parasites with all other haplotypes. Parasites engineered for another mutation, PfMDR1 1246Y, were unstable in culture, with reversion to WT over time. Thus, the N86 WT is stable when accompanied by the Y184 WT, but incurs a fitness cost when accompanied by mutant 184F. Culturing in the presence of chloroquine favored 86Y mutant parasites and in the presence of lumefantrine favored N86 WT parasites; piperaquine had minimal impact.ConclusionsThese results are consistent with those for Ugandan field isolates, suggest reasons for varied haplotypes, and highlight the interplay between drug pressure and fitness that is guiding the evolution of resistance-mediating haplotypes in P. falciparum.

Influenza A (IAV) is a major human respiratory pathogen that causes illness, hospitalizations, and mortality annually worldwide. IAV is also a zoonotic pathogen with a multitude of hosts, allowing for interspecies transmission, reassortment events, and the emergence of novel pandemics, as was seen in 2009 with the emergence of a swine-origin H1N1 (pdmH1N1) virus into humans, causing the first influenza pandemic of the 21st century. While the 2009 pandemic was considered to have high morbidity and low mortality, studies have linked the pdmH1N1 virus and its gene segments to increased disease in humans and animal models. Genetic components of the pdmH1N1 virus currently circulate in the swine population, reassorting with endemic swine viruses that co-circulate and occasionally spillover into humans. This is evidenced by the regular detection of variant swine IAVs in humans associated with state fairs and other intersections of humans and swine. Defining genetic changes that support species adaptation, virulence, and cross-species transmission, as well as mutations that enhance or attenuate these features, will improve our understanding of influenza biology. It aids in surveillance and virus risk assessment and guides the establishment of counter measures for emerging viruses. Here, we review the current understanding of the determinants of specific IAV phenotypes, focusing on the fitness, transmission, and virulence determinants that have been identified in swine IAVs and/or in relation to the 2009 pdmH1N1 virus. [ABSTRACT FROM AUTHOR]

Historical contingency and diminishing returns epistasis have been typically studied for relatively divergent genotypes and/or over long evolutionary timescales. Here, we use Saccharomyces cerevisiae to study the extent of diminishing returns and the changes in the adaptive mutational spectra following a single first adaptive mutational step. We further evolve three clones that arose under identical conditions from a common ancestor. We follow their evolutionary dynamics by lineage tracking and determine adaptive outcomes using fitness assays and whole genome sequencing. We find that diminishing returns manifests as smaller fitness gains during the 2nd step of adaptation compared to the 1st step, mainly due to a compressed distribution of fitness effects. We also find that the beneficial mutational spectra for the 2nd adaptive step are contingent on the 1st step, as we see both shared and diverging adaptive strategies. Finally, we find that adaptive loss-of-function mutations, such as nonsense and frameshift mutations, are less common in the second step of adaptation than in the first step.

BackgroundDespite increasing documentation of insecticide resistance in malaria vectors against public health insecticides in sub-Saharan Africa, there is a paucity of information on the potential fitness costs of pyrethroid resistance in malaria vectors, which is important in improving the current resistant management strategies. This study aimed to assess the fitness cost effects of insecticide resistance on the development and survival of immature Anopheles gambiae from western Kenya.MethodsTwo-hour old, first instar larvae (L1) were introduced and raised in basins containing soil and rainwater in a semi-field set-up. Each day the number of surviving individuals per larval stage was counted and their stage of development were recorded until they emerged as adults. The larval life-history trait parameters measured include mean larval development time, daily survival and pupal emergence. Pyrethroid-resistant colony of An. gambiae sensu stricto and susceptible colony originating from the same site and with the same genetic background were used. Kisumu laboratory susceptible colony was used as a reference.ResultsThe resistant colony had a significantly longer larval development time through the developmental stages than the susceptible colony. The resistant colony took an average of 2 days longer to develop from first instar (L1) to fourth instar (L4) (8.8 ± 0.2 days) compared to the susceptible colony (6.6 ± 0.2 days). The development time from first instar to pupa formation was significantly longer by 3 days in the resistant colony (10.28 ± 0.3 days) than in susceptible colony (7.5 ± 0.2 days). The time from egg hatching to adult emergence was significantly longer for the resistant colony (12.1 ± 0.3 days) than the susceptible colony (9.6 ± 0.2 days). The pupation rate (80%; 95% (CI: 77.5-83.6) vs 83.5%; 95% (CI: 80.6-86.3)) and adult emergence rate (86.3% vs 92.8%) did not differ between the resistant and susceptible colonies, respectively. The sex ratio of the females to males for the resistant (1:1.2) and susceptible colonies (1:1.07) was significantly different.ConclusionThe study showed that pyrethroid resistance in An. gambiae had a fitness cost on their pre-imaginal development time and survival. Insecticide resistance delayed the development and reduced the survivorship of An. gambiae larvae. The study findings are important in understanding the fitness cost of insecticide resistance vectors that could contribute to shaping resistant management strategies.

Understanding when severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) emerged is critical to evaluating our current approach to monitoring novel zoonotic pathogens and understanding the failure of early containment and mitigation efforts for COVID-19. We used a coalescent framework to combine retrospective molecular clock inference with forward epidemiological simulations to determine how long SARS-CoV-2 could have circulated before the time of the most recent common ancestor of all sequenced SARS-CoV-2 genomes. Our results define the period between mid-October and mid-November 2019 as the plausible interval when the first case of SARS-CoV-2 emerged in Hubei province, China. By characterizing the likely dynamics of the virus before it was discovered, we show that more than two-thirds of SARS-CoV-2-like zoonotic events would be self-limited, dying out without igniting a pandemic. Our findings highlight the shortcomings of zoonosis surveillance approaches for detecting highly contagious pathogens with moderate mortality rates.

Domestic dogs have experienced population bottlenecks, recent inbreeding, and strong artificial selection. These processes have simplified the genetic architecture of complex traits, allowed deleterious variation to persist, and increased both identity-by-descent (IBD) segments and runs of homozygosity (ROH). As such, dogs provide an excellent model for examining how these evolutionary processes influence disease. We assembled a dataset containing 4,414 breed dogs, 327 village dogs, and 380 wolves genotyped at 117,288 markers and data for clinical and morphological phenotypes. Breed dogs have an enrichment of IBD and ROH, relative to both village dogs and wolves, and we use these patterns to show that breed dogs have experienced differing severities of bottlenecks in their recent past. We then found that ROH burden is associated with phenotypes in breed dogs, such as lymphoma. We next test the prediction that breeds with greater ROH have more disease alleles reported in the Online Mendelian Inheritance in Animals (OMIA). Surprisingly, the number of causal variants identified correlates with the popularity of that breed rather than the ROH or IBD burden, suggesting an ascertainment bias in OMIA. Lastly, we use the distribution of ROH across the genome to identify genes with depletions of ROH as potential hotspots for inbreeding depression and find multiple exons where ROH are never observed. Our results suggest that inbreeding has played a large role in shaping genetic and phenotypic variation in dogs and that future work on understudied breeds may reveal new disease-causing variation.

SARS-CoV-2 can mutate and evade immunity, with consequences for efficacy of emerging vaccines and antibody therapeutics. Here, we demonstrate that the immunodominant SARS-CoV-2 spike (S) receptor binding motif (RBM) is a highly variable region of S and provide epidemiological, clinical, and molecular characterization of a prevalent, sentinel RBM mutation, N439K. We demonstrate N439K S protein has enhanced binding affinity to the hACE2 receptor, and N439K viruses have similar in vitro replication fitness and cause infections with similar clinical outcomes as compared to wild type. We show the N439K mutation confers resistance against several neutralizing monoclonal antibodies, including one authorized for emergency use by the US Food and Drug Administration (FDA), and reduces the activity of some polyclonal sera from persons recovered from infection. Immune evasion mutations that maintain virulence and fitness such as N439K can emerge within SARS-CoV-2 S, highlighting the need for ongoing molecular surveillance to guide development and usage of vaccines and therapeutics.

Microbial interactions are expected to be major determinants of microbiome structure and function. Although fungi are found in diverse microbiomes, their interactions with bacteria remain largely uncharacterized. In this work, we characterize interactions in 16 different bacterial-fungal pairs, examining the impacts of 8 different fungi isolated from cheese rind microbiomes on 2 bacteria (Escherichia coli and a cheese-isolated Pseudomonas psychrophila). Using random barcode transposon-site sequencing with an analysis pipeline that allows statistical comparisons between different conditions, we observed that fungal partners caused widespread changes in the fitness of bacterial mutants compared to growth alone. We found that all fungal species modulated the availability of iron and biotin to bacterial species, which suggests that these may be conserved drivers of bacterial-fungal interactions. Species-specific interactions were also uncovered, a subset of which suggested fungal antibiotic production. Changes in both conserved and species-specific interactions resulted from the deletion of a global regulator of fungal specialized metabolite production. This work highlights the potential for broad impacts of fungi on bacterial species within microbiomes.

Streptococcus pneumoniae is an opportunistic human pathogen that causes invasive diseases, including pneumonia, with greater health risks upon influenza A virus (IAV) co-infection. To facilitate pathogenesis studies in vivo, we developed an inducible CRISPR interference system that enables genome-wide fitness testing in one sequencing step (CRISPRi-seq). We applied CRISPRi-seq to assess bottlenecks and identify pneumococcal genes important in a murine pneumonia model. A critical bottleneck occurs at 48 h with few bacteria causing systemic infection. This bottleneck is not present during IAV superinfection, facilitating identification of pneumococcal pathogenesis-related genes. Top in vivo essential genes included purA, encoding adenylsuccinate synthetase, and the cps operon required for capsule production. Surprisingly, CRISPRi-seq indicated no fitness-related role for pneumolysin during superinfection. Interestingly, although metK (encoding S-adenosylmethionine synthetase) was essential in vitro, it was dispensable in vivo. This highlights advantages of CRISPRi-seq over transposon-based genetic screens, as all genes, including essential genes, can be tested for pathogenesis potential.

The fitness landscape is a concept commonly used to describe evolution towards optimal phenotypes. It can be reduced to mechanistic detail using genome-scale models (GEMs) from systems biology. We use recently developed GEMs of Metabolism and protein Expression (ME-models) to study the distribution of Escherichia coli phenotypes on the rate-yield plane. We found that the measured phenotypes distribute non-uniformly to form a highly stratified fitness landscape. Systems analysis of the ME-model simulations suggest that this stratification results from discrete ATP generation strategies. Accordingly, we define "aero-types", a phenotypic trait that characterizes how a balanced proteome can achieve a given growth rate by modulating 1) the relative utilization of oxidative phosphorylation, glycolysis, and fermentation pathways; and 2) the differential employment of electron-transport-chain enzymes. This global, quantitative, and mechanistic systems biology interpretation of fitness landscape formed upon proteome allocation offers a fundamental understanding of bacterial physiology and evolution dynamics.

Dengue virus (DENV) cycles between mosquito and mammalian hosts. To examine how DENV populations adapt to these different host environments, we used serial passage in human and mosquito cell lines and estimated fitness effects for all single-nucleotide variants in these populations using ultra-deep sequencing. This allowed us to determine the contributions of beneficial and deleterious mutations to the collective fitness of the population. Our analysis revealed that the continuous influx of a large burden of deleterious mutations counterbalances the effect of rare, host-specific beneficial mutations to shape the path of adaptation. Beneficial mutations preferentially map to intrinsically disordered domains in the viral proteome and cluster to defined regions in the genome. These phenotypically redundant adaptive alleles may facilitate host-specific DENV adaptation. Importantly, the evolutionary constraints described in our simple system mirror trends observed across DENV and Zika strains, indicating it recapitulates key biophysical and biological constraints shaping long-term viral evolution.

Homing-based gene drives, engineered using CRISPR/Cas9, have been proposed to spread desirable genes throughout populations. However, invasion of such drives can be hindered by the accumulation of resistant alleles. To limit this obstacle, we engineer a confinable population modification home-and-rescue (HomeR) drive in Drosophila targeting an essential gene. In our experiments, resistant alleles that disrupt the target gene function were recessive lethal and therefore disadvantaged. We demonstrate that HomeR can achieve an increase in frequency in population cage experiments, but that fitness costs due to the Cas9 insertion limit drive efficacy. Finally, we conduct mathematical modeling comparing HomeR to contemporary gene drive architectures for population modification over wide ranges of fitness costs, transmission rates, and release regimens. HomeR could potentially be adapted to other species, as a means for safe, confinable, modification of wild populations.

The marine copepod, Tigriopus californicus, produces the red carotenoid pigment astaxanthin from yellow dietary precursors. This 'bioconversion' of yellow carotenoids to red is hypothesized to be linked to individual condition, possibly through shared metabolic pathways with mitochondrial oxidative phosphorylation. Experimental inter-population crosses of lab-reared T. californicus typically produces low-fitness hybrids is due in large part to the disruption of coadapted sets nuclear and mitochondrial genes within the parental populations. These hybrid incompatibilities can increase variability in life history traits and energy production among hybrid lines. Here, we tested if production of astaxanthin was compromised in hybrid copepods and if it was linked to mitochondrial metabolism and offspring development. We observed no clear mitonuclear dysfunction in hybrids fed a limited, carotenoid-deficient diet of nutritional yeast. However, when yellow carotenoids were restored to their diet, hybrid lines produced less astaxanthin than parental lines. We observed that lines fed a yeast diet produced less ATP and had slower offspring development compared to lines fed a more complete diet of algae, suggesting the yeast-only diet may have obscured effects of mitonuclear dysfunction. Astaxanthin production was not significantly associated with development among lines fed a yeast diet but was negatively related to development in early generation hybrids fed an algal diet. In lines fed yeast, astaxanthin was negatively related to ATP synthesis, but in lines fed algae, the relationship was reversed. Although the effects of the yeast diet may have obscured evidence of hybrid dysfunction, these results suggest that astaxanthin bioconversion may still be related to mitochondrial performance and reproductive success.

Natural selection is an important driver of genetic and phenotypic differentiation between species. For species in which potential gene flow is high but realized gene flow is low, adaptation via natural selection may be a particularly important force maintaining species. For a recent radiation of New World desert shrubs (Encelia: Asteraceae), we use fine-scale geographic sampling and population genomics to determine patterns of gene flow across two hybrid zones formed between two independent pairs of species with parapatric distributions. After finding evidence for extremely strong selection at both hybrid zones, we use a combination of field experiments, high-resolution imaging, and physiological measurements to determine the ecological basis for selection at one of the hybrid zones. Our results identify multiple ecological mechanisms of selection (drought, salinity, herbivory, and burial) that together are sufficient to maintain species boundaries despite high rates of hybridization. Given that multiple pairs of Encelia species hybridize at ecologically divergent parapatric boundaries, such mechanisms may maintain species boundaries throughout Encelia.

A recent hypothesis proposed that the total number of stem cell divisions in a tissue (TSCD model) determine its intrinsic cancer risk; however, a different model-the multistage model-has long been used to understand how cancer originates. Identifying the correct model has important implications for interpreting the frequency of cancers. Using worldwide cancer incidence data, we applied three tests to the TSCD model and an evolutionary multistage model of carcinogenesis (EMMC), a model in which cancer suppression is recognized as an evolving trait, with natural selection acting to suppress cancers causing a significant mean loss of Darwinian fitness. Each test supported the EMMC but contradicted the TSCD model. This outcome undermines results based on the TSCD model quantifying the relative importance of 'bad luck' (the random accumulation of somatic mutations) versus environmental and genetic factors in determining cancer incidence. Our testing supported the EMMC prediction that cancers of large rapidly dividing tissues predominate late in life. Another important prediction is that an indicator of recent oncogenic environmental change is an unusually high mean fitness loss due to cancer, rather than a high lifetime incidence. The evolutionary model also predicts that large and/or long-lived animals have evolved mechanisms of cancer suppression that may be of value in preventing or controlling human cancers.

Heteranthery, the presence of two or more anther types in the same flower, is taxonomically widespread among bee-pollinated angiosperms, yet has puzzled botanists since Darwin. We test two competing hypotheses for its evolution: the long-standing 'division of labour' hypothesis, which posits that some anthers are specialized as food rewards for bees whereas others are specialized for surreptitious pollination, and our new hypothesis that heteranthery is a way to gradually release pollen that maximizes pollen delivery. We examine the evolution of heteranthery and associated traits across the genus Clarkia (Onagraceae) and study plant-pollinator interactions in two heterantherous Clarkia species. Across species, heteranthery is associated with bee pollination, delayed dehiscence and colour crypsis of one anther whorl, and movement of that anther whorl upon dehiscence. Our mechanistic studies in heterantherous species show that bees notice, forage on and export pollen from each anther whorl when it is dehiscing, and that heteranthery promotes pollen export. We find no support for division of labour, but multifarious evidence that heteranthery is a mechanism for gradual pollen presentation that probably evolved through indirect male-male competition for siring success.

Acquired mutations are pervasive across normal tissues. However, understanding of the processes that drive transformation of certain clones to cancer is limited. Here we study this phenomenon in the context of clonal hematopoiesis (CH) and the development of therapy-related myeloid neoplasms (tMNs). We find that mutations are selected differentially based on exposures. Mutations in ASXL1 are enriched in current or former smokers, whereas cancer therapy with radiation, platinum and topoisomerase II inhibitors preferentially selects for mutations in DNA damage response genes (TP53, PPM1D, CHEK2). Sequential sampling provides definitive evidence that DNA damage response clones outcompete other clones when exposed to certain therapies. Among cases in which CH was previously detected, the CH mutation was present at tMN diagnosis. We identify the molecular characteristics of CH that increase risk of tMN. The increasing implementation of clinical sequencing at diagnosis provides an opportunity to identify patients at risk of tMN for prevention strategies.

Generalizing the effect of traits on performance across species may be achievable if traits explain variation in population fitness. However, testing relationships between traits and vital rates to infer effects on fitness can be misleading. Demographic trade-offs can generate variation in vital rates that yield equal population growth rates, thereby obscuring the net effect of traits on fitness. To address this problem, we describe a diversity of approaches to quantify intrinsic growth rates of plant populations, including experiments beyond range boundaries, density-dependent population models built from long-term demographic data, theoretical models, and methods that leverage widely available monitoring data. Linking plant traits directly to intrinsic growth rates is a fundamental step toward rigorous predictions of population dynamics and community assembly.

Mutant evolution in spatially structured systems is important for a range of biological systems, but aspects of it still require further elucidation. Adding to previous work, we provide a simple derivation of growth laws that characterize the number of mutants of different relative fitness in expanding populations in spatial models of different dimensionalities. These laws are universal and independent of "microscopic" modeling details. We further study the accumulation of mutants and find that, with advantageous and neutral mutants, more of them are present in spatially structured, compared to well-mixed colonies of the same size. The behavior of disadvantageous mutants is subtle: if they are disadvantageous through a reduction in division rates, the result is the same, and it is the opposite if the disadvantage is due to a death rate increase. Finally, we show that in all cases, the same results are observed in fragmented, nonspatial patch models. This suggests that the patterns observed are the consequence of population fragmentation, and not spatial restrictions per se We provide an intuitive explanation for the complex dependence of disadvantageous mutant evolution on spatial restriction, which relies on desynchronized dynamics in different locations/patches, and plays out differently depending on whether the disadvantage is due to a lower division rate or a higher death rate. Implications for specific biological systems, such as the evolution of drug-resistant cell mutants in cancer or bacterial biofilms, are discussed.

Adaptation via natural selection is an important driver of evolution, and repeatable adaptations of replicate populations, under conditions of a constant environment, have been extensively reported. However, isolated groups of populations in nature tend to harbor both genetic and physiological divergence due to multiple selective pressures that they have encountered. How this divergence affects adaptation of these populations to a new common environment remains unclear. To determine the impact of prior genetic and physiological divergence in shaping adaptive evolution to accommodate a new common environment, an experimental evolution study with the sulfate-reducing bacterium Desulfovibrio vulgaris Hildenborough (DvH) was conducted. Two groups of replicate populations with genetic and physiological divergence, derived from a previous evolution study, were propagated in an elevated-temperature environment for 1,000 generations. Ancestor populations without prior experimental evolution were also propagated in the same environment as a control. After 1,000 generations, all the populations had increased growth rates and all but one had greater fitness in the new environment than the ancestor population. Moreover, improvements in growth rate were moderately affected by the divergence in the starting populations, while changes in fitness were not significantly affected. The mutations acquired at the gene level in each group of populations were quite different, indicating that the observed phenotypic changes were achieved by evolutionary responses that differed between the groups. Overall, our work demonstrated that the initial differences in fitness between the starting populations were eliminated by adaptation and that phenotypic convergence was achieved by acquisition of mutations in different genes.IMPORTANCE Improving our understanding of how previous adaptation influences evolution has been a long-standing goal in evolutionary biology. Natural selection tends to drive populations to find similar adaptive solutions for the same selective conditions. However, variations in historical environments can lead to both physiological and genetic divergence that can make evolution unpredictable. Here, we assessed the influence of divergence on the evolution of a model sulfate-reducing bacterium, Desulfovibrio vulgaris Hildenborough, in response to elevated temperature and found a significant effect at the genetic but not the phenotypic level. Understanding how these influences drive evolution will allow us to better predict how bacteria will adapt to various ecological constraints.

Species-specific limits to lifespan (lifespan setpoint) determine the life expectancy of any given organism. Whether limiting lifespan provides an evolutionary benefit or is the result of an inevitable decline in fitness remains controversial. The identification of mutations extending lifespan suggests that aging is under genetic control, but the evolutionary driving forces limiting lifespan have not been defined. By examining the impact of lifespan on pathogen spread in a population, we propose that epidemics drive lifespan setpoints' evolution. Shorter lifespan limits infection spread and accelerates pathogen clearance when compared to populations with longer-lived individuals. Limiting longevity is particularly beneficial in the context of zoonotic transmissions, where pathogens must undergo adaptation to a new host. Strikingly, in populations exposed to pathogens, shorter-living variants outcompete individuals with longer lifespans. We submit that infection outbreaks can contribute to control the evolution of species' lifespan setpoints.

Ecological speciation occurs when reproductive isolation evolves as a byproduct of adaptive divergence between populations. Selection favouring gene regulatory divergence between species could result in transgressive levels of gene expression in F1 hybrids that may lower hybrid fitness. We combined 58 resequenced genomes with 124 transcriptomes to identify patterns of hybrid gene misexpression that may be driven by adaptive regulatory divergence within a young radiation of Cyprinodon pupfishes, which consists of a dietary generalist and two trophic specialists-a molluscivore and a scale-eater. We found more differential gene expression between closely related sympatric specialists than between allopatric generalist populations separated by 1,000 km. Intriguingly, 9.6% of genes that were differentially expressed between sympatric species were also misexpressed in F1 hybrids. A subset of these genes were in highly differentiated genomic regions and enriched for functions important for trophic specialization, including head, muscle and brain development. These regions also included genes that showed evidence of hard selective sweeps and were significantly associated with oral jaw length-the most rapidly diversifying skeletal trait in this radiation. Our results indicate that divergent ecological selection in sympatry can contribute to hybrid gene misexpression which may act as a reproductive barrier between nascent species.

Y chromosomes are typically viewed as genetic wastelands with few intact genes. Recent genomic analyses in Drosophila, however, show that gene gain is prominent on young Y chromosomes. Meiosis- and RNAi-related genes often coamplify on recently formed X and Y chromosomes, are testis-expressed, and produce antisense transcripts and short RNAs. RNAi pathways are also involved in suppressing sex ratio drive in Drosophila. These observations paint a dynamic picture of sex chromosome differentiation, suggesting that rapidly evolving genomic battles over segregation are rampant on young sex chromosomes and utilize RNAi to defend the genome against selfish elements that manipulate fair meiosis. Recurrent sex chromosome drive can have profound ecological, evolutionary, and cellular impacts and account for unique features of sex chromosomes.

Microbial interactions abound in natural ecosystems and shape community structure and function. Substantial attention has been given to cataloging mechanisms by which microbes interact, but there is a limited understanding of the genetic landscapes that promote or hinder microbial interactions. We previously developed a mutualistic coculture pairing Escherichia coli and Rhodopseudomonas palustris, wherein E. coli provides carbon to R. palustris in the form of glucose fermentation products and R. palustris fixes N2 gas and provides nitrogen to E. coli in the form of NH4+ The stable coexistence and reproducible trends exhibited by this coculture make it ideal for interrogating the genetic underpinnings of a cross-feeding mutualism. Here, we used random barcode transposon sequencing (RB-TnSeq) to conduct a genome-wide search for E. coli genes that influence fitness during cooperative growth with R. palustris RB-TnSeq revealed hundreds of genes that increased or decreased E. coli fitness in a mutualism-dependent manner. Some identified genes were involved in nitrogen sensing and assimilation, as expected given the coculture design. The other identified genes were involved in diverse cellular processes, including energy production and cell wall and membrane biogenesis. In addition, we discovered unexpected purine cross-feeding from R. palustris to E. coli, with coculture rescuing growth of an E. coli purine auxotroph. Our data provide insight into the genes and gene networks that can influence a cross-feeding mutualism and underscore that microbial interactions are not necessarily predictable a prioriIMPORTANCE Microbial communities impact life on Earth in profound ways, including driving global nutrient cycles and influencing human health and disease. These community functions depend on the interactions that resident microbes have with the environment and each other. Thus, identifying genes that influence these interactions will aid the management of natural communities and the use of microbial consortia as biotechnology. Here, we identified genes that influenced Escherichia coli fitness during cooperative growth with a mutualistic partner, Rhodopseudomonas palustris Although this mutualism centers on the bidirectional exchange of essential carbon and nitrogen, E. coli fitness was positively and negatively affected by genes involved in diverse cellular processes. Furthermore, we discovered an unexpected purine cross-feeding interaction. These results contribute knowledge on the genetic foundation of a microbial cross-feeding interaction and highlight that unanticipated interactions can occur even within engineered microbial communities.

BACKGROUND: Viridans group streptococci of the Streptococcus mitis-oralis subgroup are important endovascular pathogens. They can rapidly develop high-level and durable non-susceptibility to daptomycin both in vitro and in vivo upon exposure to daptomycin. Two consistent genetic adaptations associated with this phenotype (i.e., mutations in cdsA and pgsA) lead to the depletion of the phospholipids, phosphatidylglycerol and cardiolipin, from the bacterial membrane. Such alterations in phospholipid biosynthesis will modify carbon flow and change the bacterial metabolic status. To determine the metabolic differences between daptomycin-susceptible and non-susceptible bacteria, the physiology and metabolomes of S. mitis-oralis strains 351 (daptomycin-susceptible) and 351-D10 (daptomycin non-susceptible) were analyzed. S. mitis-oralis strain 351-D10 was made daptomycin non-susceptible through serial passage in the presence of daptomycin. RESULTS: Daptomycin non-susceptible S. mitis-oralis had significant alterations in glucose catabolism and a re-balancing of the redox status through amino acid biosynthesis relative to daptomycin susceptible S. mitis-oralis. These changes were accompanied by a reduced capacity to generate biomass, creating a fitness cost in exchange for daptomycin non-susceptibility. CONCLUSIONS: S. mitis-oralis metabolism is altered in daptomycin non-susceptible bacteria relative to the daptomycin susceptible parent strain. As demonstrated in Staphylococcus aureus, inhibiting the metabolic changes that facilitate the transition from a daptomycin susceptible state to a non-susceptible one, inhibits daptomycin non-susceptibility. By preventing these metabolic adaptations in S. mitis-oralis, it should be possible to deter the formation of daptomycin non-susceptibility.

Oxidative phosphorylation, the primary source of cellular energy in eukaryotes, requires gene products encoded in both the nuclear and mitochondrial genomes. As a result, functional integration between the genomes is essential for efficient adenosine triphosphate (ATP) generation. Although within populations this integration is presumably maintained by coevolution, the importance of mitonuclear coevolution in key biological processes such as speciation and mitochondrial disease has been questioned. In this study, we crossed populations of the intertidal copepod Tigriopus californicus to disrupt putatively coevolved mitonuclear genotypes in reciprocal F2 hybrids. We utilized interindividual variation in developmental rate among these hybrids as a proxy for fitness to assess the strength of selection imposed on the nuclear genome by alternate mitochondrial genotypes. Developmental rate varied among hybrid individuals, and in vitro ATP synthesis rates of mitochondria isolated from high-fitness hybrids were approximately two-fold greater than those of mitochondria isolated from low-fitness individuals. We then used Pool-seq to compare nuclear allele frequencies for high- or low-fitness hybrids. Significant biases for maternal alleles were detected on 5 (of 12) chromosomes in high-fitness individuals of both reciprocal crosses, whereas maternal biases were largely absent in low-fitness individuals. Therefore, the most fit hybrids were those with nuclear alleles that matched their mitochondrial genotype on these chromosomes, suggesting that mitonuclear effects underlie individual-level variation in developmental rate and that intergenomic compatibility is critical for high fitness. We conclude that mitonuclear interactions can have profound impacts on both physiological performance and the evolutionary trajectory of the nuclear genome.

Oxidative phosphorylation, the primary source of cellular energy in eukaryotes, requires gene products encoded in both the nuclear and mitochondrial genomes. As a result, functional integration between the genomes is essential for efficient adenosine triphosphate (ATP) generation. Although within populations this integration is presumably maintained by coevolution, the importance of mitonuclear coevolution in key biological processes such as speciation and mitochondrial disease has been questioned. In this study, we crossed populations of the intertidal copepod Tigriopus californicus to disrupt putatively coevolved mitonuclear genotypes in reciprocal F2 hybrids. We utilized interindividual variation in developmental rate among these hybrids as a proxy for fitness to assess the strength of selection imposed on the nuclear genome by alternate mitochondrial genotypes. Developmental rate varied among hybrid individuals, and in vitro ATP synthesis rates of mitochondria isolated from high-fitness hybrids were approximately two-fold greater than those of mitochondria isolated from low-fitness individuals. We then used Pool-seq to compare nuclear allele frequencies for high- or low-fitness hybrids. Significant biases for maternal alleles were detected on 5 (of 12) chromosomes in high-fitness individuals of both reciprocal crosses, whereas maternal biases were largely absent in low-fitness individuals. Therefore, the most fit hybrids were those with nuclear alleles that matched their mitochondrial genotype on these chromosomes, suggesting that mitonuclear effects underlie individual-level variation in developmental rate and that intergenomic compatibility is critical for high fitness. We conclude that mitonuclear interactions can have profound impacts on both physiological performance and the evolutionary trajectory of the nuclear genome.

Levels of gene expression underpin organismal phenotypes1,2, but the nature of selection that acts on gene expression and its role in adaptive evolution remain unknown1,2. Here we assayed gene expression in rice (Oryza sativa)3, and used phenotypic selection analysis to estimate the type and strength of selection on the levels of more than 15,000 transcripts4,5. Variation in most transcripts appears (nearly) neutral or under very weak stabilizing selection in wet paddy conditions (with median standardized selection differentials near zero), but selection is stronger under drought conditions. Overall, more transcripts are conditionally neutral (2.83%) than are antagonistically pleiotropic6 (0.04%), and transcripts that display lower levels of expression and stochastic noise7-9 and higher levels of plasticity9 are under stronger selection. Selection strength was further weakly negatively associated with levels of cis-regulation and network connectivity9. Our multivariate analysis suggests that selection acts on the expression of photosynthesis genes4,5, but that the efficacy of selection is genetically constrained under drought conditions10. Drought selected for earlier flowering11,12 and a higher expression of OsMADS18 (Os07g0605200), which encodes a MADS-box transcription factor and is a known regulator of early flowering13-marking this gene as a drought-escape gene11,12. The ability to estimate selection strengths provides insights into how selection can shape molecular traits at the core of gene action.

Our ability to predict how species will respond to human-induced rapid environmental change (HIREC) may depend upon our understanding of transgenerational plasticity (TGP), which occurs when environments experienced by previous generations influence phenotypes of subsequent generations. TGP evolved to help organisms cope with environmental stressors when parental environments are highly predictive of offspring environments. HIREC can alter conditions that favored TGP in historical environments by reducing parents' ability to detect environmental conditions, disrupting previous correlations between parental and offspring environments, and interfering with the transmission of parental cues to offspring. Because of the propensity to produce errors in these processes, TGP will likely generate negative fitness outcomes in response to HIREC, though beneficial fitness outcomes may occur in some cases.

Cancer progression is an evolutionary process. During this process, evolving cancer cell populations encounter restrictive ecological niches within the body, such as the primary tumor, circulatory system, and diverse metastatic sites. Efforts to prevent or delay cancer evolution-and progression-require a deep understanding of the underlying molecular evolutionary processes. Herein we discuss a suite of concepts and tools from evolutionary and ecological theory that can inform cancer biology in new and meaningful ways. We also highlight current challenges to applying these concepts, and propose ways in which incorporating these concepts could identify new therapeutic modes and vulnerabilities in cancer.

The extent to which sequence variation impacts plant fitness is poorly understood. High-resolution maps detailing the constraint acting on the genome, especially in regulatory sites, would be beneficial as functional annotation of noncoding sequences remains sparse. Here, we present a fitness consequence (fitCons) map for rice (Oryza sativa). We inferred fitCons scores (ρ) for 246 inferred genome classes derived from nine functional genomic and epigenomic datasets, including chromatin accessibility, messenger RNA/small RNA transcription, DNA methylation, histone modifications and engaged RNA polymerase activity. These were integrated with genome-wide polymorphism and divergence data from 1,477 rice accessions and 11 reference genome sequences in the Oryzeae. We found ρ to be multimodal, with ~9% of the rice genome falling into classes where more than half of the bases would probably have a fitness consequence if mutated. Around 2% of the rice genome showed evidence of weak negative selection, frequently at candidate regulatory sites, including a novel set of 1,000 potentially active enhancer elements. This fitCons map provides perspective on the evolutionary forces associated with genome diversity, aids in genome annotation and can guide crop breeding programs.

Consumers face numerous risks that can be minimized by incorporating different life-history strategies. How much and when a consumer adds to its energetic reserves or invests in reproduction are key behavioral and physiological adaptations that structure communities. Here we develop a theoretical framework that explicitly accounts for stochastic fluctuations of an individual consumers energetic reserves while foraging and reproducing on a landscape with resources that range from uniformly distributed to highly clustered. First, we show that the selection of alternative life histories depends on both the mean and variance of resource availability, where depleted and more stochastic environments promote investment in each reproductive event at the expense of future fitness as well as more investment per offspring. We then show that if resource variance scales with body size due to landscape clustering, consumers that forage for clustered foods are susceptible to strong Allee effects, increasing extinction risk. Finally, we show that the proposed relationship between resource distributions, consumer body size, and emergent demographic risk offers key ecological insights into the evolution of large-bodied grazing herbivores from small-bodied browsing ancestors.

Vaccination has essentially eradicated poliovirus. Yet, its mutation rate is higher than that of viruses like HIV, for which no effective vaccine exists. To investigate this, we infer a fitness model for the poliovirus viral protein 1 (vp1), which successfully predicts in vitro fitness measurements. This is achieved by first developing a probabilistic model for the prevalence of vp1 sequences that enables us to isolate and remove data that are subject to strong vaccine-derived biases. The intrinsic fitness constraints derived for vp1, a capsid protein subject to antibody responses, are compared with those of analogous HIV proteins. We find that vp1 evolution is subject to tighter constraints, limiting its ability to evade vaccine-induced immune responses. Our analysis also indicates that circulating poliovirus strains in unimmunized populations serve as a reservoir that can seed outbreaks in spatio-temporally localized sub-optimally immunized populations.

Drug-resistant mutations often have deleterious impacts on replication fitness, posing a fitness cost that can only be overcome by compensatory mutations. However, the role of fitness cost in the evolution of drug resistance has often been overlooked in clinical studies or in vitro selection experiments, as these observations only capture the outcome of drug selection. In this study, we systematically profile the fitness landscape of resistance-associated sites in HIV-1 protease using deep mutational scanning. We construct a mutant library covering combinations of mutations at 11 sites in HIV-1 protease, all of which are associated with resistance to protease inhibitors in clinic. Using deep sequencing, we quantify the fitness of thousands of HIV-1 protease mutants after multiple cycles of replication in human T cells. Although the majority of resistance-associated mutations have deleterious effects on viral replication, we find that epistasis among resistance-associated mutations is predominantly positive. Furthermore, our fitness data are consistent with genetic interactions inferred directly from HIV sequence data of patients. Fitness valleys formed by strong positive epistasis reduce the likelihood of reversal of drug resistance mutations. Overall, our results support the view that strong compensatory effects are involved in the emergence of clinically observed resistance mutations and provide insights to understanding fitness barriers in the evolution and reversion of drug resistance.

Consumers face numerous risks that can be minimized by incorporating different life-history strategies. How much and when a consumer adds to its energetic reserves or invests in reproduction are key behavioral and physiological adaptations that structure communities. Here we develop a theoretical framework that explicitly accounts for stochastic fluctuations of an individual consumer's energetic reserves while foraging and reproducing on a landscape with resources that range from uniformly distributed to highly clustered. First, we show that the selection of alternative life histories depends on both the mean and variance of resource availability, where depleted and more stochastic environments promote investment in each reproductive event at the expense of future fitness as well as more investment per offspring. We then show that if resource variance scales with body size due to landscape clustering, consumers that forage for clustered foods are susceptible to strong Allee effects, increasing extinction risk. Finally, we show that the proposed relationship between resource distributions, consumer body size, and emergent demographic risk offers key ecological insights into the evolution of large-bodied grazing herbivores from small-bodied browsing ancestors.

Life history theory examines how characteristics of organisms, such as age and size at maturity, may vary through natural selection as evolutionary responses that optimize fitness. Here we ask how predictions of age and size at maturity differ for the three classical fitness functions-intrinsic rate of natural increase r, net reproductive rate R0, and reproductive value Vx-for semelparous species. We show that different choices of fitness functions can lead to very different predictions of species behavior. In one's efforts to understand an organism's behavior and to develop effective conservation and management policies, the choice of fitness function matters. The central ingredient of our approach is the maturation reaction norm (MRN), which describes how optimal age and size at maturation vary with growth rate or mortality rate. We develop a practical geometric construction of MRNs that allows us to include different growth functions (linear growth and nonlinear von Bertalanffy growth in length) and develop two-dimensional MRNs useful for quantifying growth-mortality trade-offs. We relate our approach to Beverton-Holt life history invariants and to the Stearns-Koella categorization of MRNs. We conclude with a detailed discussion of life history parameters for Great Lakes Chinook Salmon and demonstrate that age and size at maturity are consistent with predictions using R0 (but not r or Vx) as the underlying fitness function.

Wolbachia is an endosymbiotic Alphaproteobacteria that can suppress insect-borne diseases through decreasing host virus transmission (population replacement) or through decreasing host population density (population suppression). We contrast natural Wolbachia infections in insect populations with Wolbachia transinfections in mosquitoes to gain insights into factors potentially affecting the long-term success of Wolbachia releases. Natural Wolbachia infections can spread rapidly, whereas the slow spread of transinfections is governed by deleterious effects on host fitness and demographic factors. Cytoplasmic incompatibility (CI) generated by Wolbachia is central to both population replacement and suppression programs, but CI in nature can be variable and evolve, as can Wolbachia fitness effects and virus blocking. Wolbachia spread is also influenced by environmental factors that decrease Wolbachia titer and reduce maternal Wolbachia transmission frequency. More information is needed on the interactions between Wolbachia and host nuclear/mitochondrial genomes, the interaction between invasion success and local ecological factors, and the long-term stability of Wolbachia-mediated virus blocking.