Influenza Viruses in Mice: Deep Sequencing Analysis of Serial Passage and Effects of Sialic Acid Structural Variation

Influenza A viruses have regularly jumped to new host species to cause epidemics or pandemics, an evolutionary process that involves variation in the viral traits necessary to overcome host barriers and facilitate transmission. Mice are not a natural host for influenza virus but are frequently used as models in studies of pathogenesis, often after multiple passages to achieve higher viral titers that result in clinical disease such as weight loss or death. Here, we examine the processes of influenza A virus infection and evolution in mice by comparing single nucleotide variations of a human H1N1 pandemic virus, a seasonal H3N2 virus, and an H3N2 canine influenza virus during experimental passage. We also compared replication and sequence variation in wild-type mice expressing N-glycolylneuraminic acid (Neu5Gc) with those seen in mice expressing only N-acetylneuraminic acid (Neu5Ac). Viruses derived from plasmids were propagated in MDCK cells and then passaged in mice up to four times. Full-genome deep sequencing of the plasmids, cultured viruses, and viruses from mice at various passages revealed only small numbers of mutational changes. The H3N2 canine influenza virus showed increases in frequency of sporadic mutations in the PB2, PA, and NA segments. The H1N1 pandemic virus grew well in mice, and while it exhibited the maintenance of some minority mutations, there was no clear evidence for adaptive evolution. The H3N2 seasonal virus did not establish in the mice. Finally, there were no clear sequence differences associated with the presence or absence of Neu5Gc.IMPORTANCE Mice are commonly used as a model to study the growth and virulence of influenza A viruses in mammals but are not a natural host and have distinct sialic acid receptor profiles compared to humans. Using experimental infections with different subtypes of influenza A virus derived from different hosts, we found that evolution of influenza A virus in mice did not necessarily proceed through the linear accumulation of host-adaptive mutations, that there was variation in the patterns of mutations detected in each repetition, and that the mutation dynamics depended on the virus examined. In addition, variation in the viral receptor, sialic acid, did not affect influenza virus evolution in this model. Overall, our results show that while mice provide a useful animal model for influenza virus pathology, host passage evolution will vary depending on the specific virus tested.

Comparative genomics can provide new insights into the evolutionary mechanisms and gene function in CAM plants

Crassulacean acid metabolism (CAM) photosynthesis is an important biological innovation enabling plant adaptation to hot and dry environments. CAM plants feature high water-use efficiency, with potential for sustainable crop production under water-limited conditions. A deep understanding of CAM-related gene function and molecular evolution of CAM plants is critical for exploiting the potential of engineering CAM into C3 crops to enhance crop production on semi-arid or marginal agricultural lands. With the newly emerging genomics resources for multiple CAM species, progress has been made in comparative genomics studies on the molecular basis and subsequently on the evolution of CAM. Here, recent advances in CAM comparative genomics research in constitutive and facultative CAM plants are reviewed, with a focus on the analyses of DNA/protein sequences and gene expression to provide new insights into the path and driving force of CAM evolution and to identify candidate genes involved in CAM-related biological processes. Potential applications of new computational and experimental technologies (e.g. CRISPR/Cas-mediated genome-editing technology) to the comparative and evolutionary genomics research on CAM plants are offered.

Bridging the Process-Pattern Divide to Understand the Origins and Early Stages of Adaptive Radiation: A Review of Approaches With Insights From Studies of Anolis Lizards

Understanding the origins and early stages of diversification is one of the most elusive tasks in adaptive radiation research. Classical approaches, which aim to infer past processes from present-day patterns of biological diversity, are fraught with difficulties and assumptions. An alternative approach has been to study young clades of relatively few species, which may represent the putative early stages of adaptive radiation. However, it is difficult to predict whether those groups will ever reach the ecological and morphological disparity observed in the sorts of clades usually referred to as adaptive radiations, thereby making their utility in informing the early stages of such radiations uncertain. Caribbean Anolis lizards are a textbook example of an adaptive radiation; anoles have diversified independently on each of the 4 islands in the Greater Antilles, producing replicated radiations of phenotypically diverse species. However, the underlying processes that drove these radiations occurred 30–65 million years ago and so are unobservable, rendering major questions about how these radiations came to be difficult to tackle. What did the ancestral species of the anole radiation look like? How did new species arise? What processes drove adaptive diversification? Here, we review what we have learned about the cryptic early stages of adaptive radiation from studies of Anolis lizards, and how these studies have attempted to bridge the process-pattern divide of adaptive radiation research. Despite decades of research, however, fundamental questions linking eco-evolutionary processes to macroevolutionary patterns in anoles remain difficult to answer.

What Happened to the Phycobilisome?

The phycobilisome (PBS) is the major light-harvesting complex of photosynthesis in cyanobacteria, red algae, and glaucophyte algae. In spite of the fact that it is very well structured to absorb light and transfer it efficiently to photosynthetic reaction centers, it has been completely lost in the green algae and plants. It is difficult to see how selection alone could account for such a major loss. An alternative scenario takes into account the role of chance, enabled by (contingent on) the evolution of an alternative antenna system early in the diversification of the three lineages from the first photosynthetic eukaryote.

Predicting evolutionary responses to interspecific interference in the wild

Many interspecifically territorial species interfere with each other reproductively, and in some cases, aggression towards heterospecifics may be an adaptive response to interspecific mate competition. This hypothesis was recently formalised in an agonistic character displacement (ACD) model which predicts that species should evolve to defend territories against heterospecific rivals above a threshold level of reproductive interference. To test this prediction, we parameterised the model with field estimates of reproductive interference for 32 sympatric damselfly populations and ran evolutionary simulations. Asymmetries in reproductive interference made the outcome inherently unpredictable in some cases, but 80% of the model's stable outcomes matched levels of heterospecific aggression in the field, significantly exceeding chance expectations. In addition to bolstering the evidence for ACD, this paper introduces a new, predictive approach to testing character displacement theory that, if applied to other systems, could help in resolving long-standing questions about the importance of character displacement processes in nature.

Introgression drives repeated evolution of winter coat color polymorphism in hares

Changing from summer-brown to winter-white pelage or plumage is a crucial adaptation to seasonal snow in more than 20 mammal and bird species. Many of these species maintain nonwhite winter morphs, locally adapted to less snowy conditions, which may have evolved independently. Mountain hares (Lepus timidus) from Fennoscandia were introduced into the Faroe Islands in 1855. While they were initially winter-white, within ∼65 y all Faroese hares became winter-gray, a morph that occurs in the source population at low frequency. The documented population history makes this a valuable model for understanding the genetic basis and evolution of the seasonal trait polymorphism. Through whole-genome scans of differentiation and single-nucleotide polymorphism (SNP) genotyping, we associated winter coat color polymorphism to the genomic region of the pigmentation gene Agouti, previously linked to introgression-driven winter coat color variation in the snowshoe hare (Lepus americanus). Lower Agouti expression in the skin of winter-gray individuals during the autumn molt suggests that regulatory changes may underlie the color polymorphism. Variation in the associated genomic region shows signatures of a selective sweep in the Faroese population, suggesting that positive selection drove the fixation of the variant after the introduction. Whole-genome analyses of several hare species revealed that the winter-gray variant originated through introgression from a noncolor changing species, in keeping with the history of ancient hybridization between the species. Our findings show the recurrent role of introgression in generating winter coat color variation by repeatedly recruiting the regulatory region of Agouti to modulate seasonal coat color change.

Convergent Evolution of Claw Shape in a Transcontinental Lizard Radiation

Species occupying similar selective environments often share similar phenotypes as the result of natural selection. Recent discoveries, however, have led to the understanding that phenotypes may also converge for other reasons than recurring selection. We argue that the vertebrate claw system constitutes a promising but understudied model system for testing the adaptive nature of phenotypic, functional, and genetic convergence. In this study, we combine basic morphometrics and advanced techniques in form analysis to examine claw shape divergence in a transcontinental lizard radiation (Lacertidae). We find substantial interspecific variation in claw morphology and phylogenetic comparative statistics reveal a strong correlation with structural habitat use: ground-dwelling species living in open areas are equipped with long, thick, weakly curved, slender-bodied claws, whereas climbing species carry high, short, strongly curved, full-bodied claws. Species occupying densely vegetated habitats tend to carry intermediately shaped claws. Evolutionary models suggest that claw shape evolves toward multiple adaptive peaks, with structural habitat use pulling species toward a specific selective optimum. Contrary to findings in several other vertebrate taxa, our analyses indicate that environmental pressures, not phylogenetic relatedness, drive convergent evolution of similarly shaped claws in lacertids. Overall, our study suggests that lacertids independently evolved similarly shaped claws as an adaptation to similar structural environments in order to cope with the specific locomotory challenges posed by the habitat. Future biomechanical studies that link form and function in combination with genomic and development research will prove valuable in better understanding the adaptive significance of claw shape divergence.

Public Microbial Resource Centers: Key Hubs for Findable, Accessible, Interoperable, and Reusable (FAIR) Microorganisms and Genetic Materials

In the context of open science, the availability of research materials is essential for knowledge accumulation and to maximize the impact of scientific research. In microbiology, microbial domain biological resource centers (mBRCs) have long-standing experience in preserving and distributing authenticated microbial strains and genetic materials (e.g., recombinant plasmids and DNA libraries) to support new discoveries and follow-on studies. These culture collections play a central role in the conservation of microbial biodiversity and have expertise in cultivation, characterization, and taxonomy of microorganisms. Information associated with preserved biological resources is recorded in databases and is accessible through online catalogues. Legal expertise developed by mBRCs guarantees end users the traceability and legality of the acquired material, notably with respect to the Nagoya Protocol. However, awareness of the advantages of depositing biological materials in professional repositories remains low, and the necessity of securing strains and genetic resources for future research must be emphasized. This review describes the unique position of mBRCs in microbiology and molecular biology through their history, evolving roles, expertise, services, challenges, and international collaborations. It also calls for an increased deposit of strains and genetic resources, a responsibility shared by scientists, funding agencies, and publishers. Journal policies requesting a deposit during submission of a manuscript represent one of the measures to make more biological materials available to the broader community, hence fully releasing their potential and improving openness and reproducibility in scientific research.

Historical contingency in the evolution of antibiotic resistance after decades of relaxed selection

Populations often encounter changed environments that remove selection for the maintenance of particular phenotypic traits. The resulting genetic decay of those traits under relaxed selection reduces an organism’s fitness in its prior environment. However, whether and how such decay alters the subsequent evolvability of a population upon restoration of selection for a previously diminished trait is not well understood. We addressed this question using Escherichia coli strains from the long-term evolution experiment (LTEE) that independently evolved for multiple decades in the absence of antibiotics. We first confirmed that these derived strains are typically more sensitive to various antibiotics than their common ancestor. We then subjected the ancestral and derived strains to various concentrations of these drugs to examine their potential to evolve increased resistance. We found that evolvability was idiosyncratic with respect to initial genotype; that is, the derived strains did not generally compensate for their greater susceptibility by “catching up” to the resistance level of the ancestor. Instead, the capacity to evolve increased resistance was constrained in some backgrounds, implying that evolvability depended upon prior mutations in a historically contingent fashion. We further subjected a time series of clones from one LTEE population to tetracycline and determined that an evolutionary constraint arose early in that population, corroborating the role of contingency. In summary, relaxed selection not only can drive populations to increased antibiotic susceptibility, but it can also affect the subsequent evolvability of antibiotic resistance in an unpredictable manner. This conclusion has potential implications for public health, and it underscores the need to consider the genetic context of pathogens when designing drug-treatment strategies.

Strains of Escherichia coli from the long-term evolution experiment (LTEE) have independently evolved for multiple decades and tens of thousands of generations in the absence of antibiotics; this study shows that evolution and diversification of a single bacterial species in an antibiotic-free environment can render the subsequent evolution of resistance unpredictable.

Evolution of Escherichia coli in different carbon environments for 2,000 generations

Cellular energetics is thought to have played a key role in dictating all major evolutionary transitions in the history of life on Earth. However, how exactly cellular energetics and metabolism come together to shape evolutionary paths is not well understood. In particular, when an organism is evolved in different energy environments, what are the phenomenological differences in the chosen evolutionary trajectories, is a question that is not well understood. In this context, starting from an Escherichia coli K-12 strain, we evolve the bacterium in five different carbon environments-glucose, arabinose, xylose, rhamnose and a mixture of these four sugars (in a predefined ratio) for approximately 2,000 generations. At the end of the adaptation period, we quantify and compare the growth dynamics of the strains in a variety of environments. The evolved strains show no specialized adaptation towards growth in the carbon medium in which they were evolved. Rather, in all environments, the evolved strains exhibited a reduced lag phase and an increased growth rate. Sequencing results reveal that these dynamical properties are not introduced via mutations in the precise loci associated with utilization of the sugar in which the bacterium evolved. These phenotypic changes are rather likely introduced via mutations elsewhere on the genome. Data from our experiments indicate that evolution in a defined environment does not alter hierarchy in mixed-sugar utilization in bacteria.

Estimating the predictability of cancer evolution

MOTIVATION

How predictable is the evolution of cancer? This fundamental question is of immense relevance for the diagnosis, prognosis and treatment of cancer. Evolutionary biologists have approached the question of predictability based on the underlying fitness landscape. However, empirical fitness landscapes of tumor cells are impossible to determine in vivo. Thus, in order to quantify the predictability of cancer evolution, alternative approaches are required that circumvent the need for fitness landscapes.

RESULTS

We developed a computational method based on conjunctive Bayesian networks (CBNs) to quantify the predictability of cancer evolution directly from mutational data, without the need for measuring or estimating fitness. Using simulated data derived from >200 different fitness landscapes, we show that our CBN-based notion of evolutionary predictability strongly correlates with the classical notion of predictability based on fitness landscapes under the strong selection weak mutation assumption. The statistical framework enables robust and scalable quantification of evolutionary predictability. We applied our approach to driver mutation data from the TCGA and the MSK-IMPACT clinical cohorts to systematically compare the predictability of 15 different cancer types. We found that cancer evolution is remarkably predictable as only a small fraction of evolutionary trajectories are feasible during cancer progression.

AVAILABILITY AND IMPLEMENTATION

https://github.com/cbg-ethz/predictability\_of\_cancer\_evolution.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Independent lineages in a common environment: the roles of determinism and contingency in shaping the migration timing of even- versus odd-year pink salmon over broad spatial and temporal scales

Studies of parallel evolution are seldom able to disentangle the influence of cryptic environmental variation from that of evolutionary history; whereas the unique life history of pink salmon (Oncorhynchus gorbuscha) presents an opportunity to do so. All pink salmon mature at age two and die after breeding. Hence, pink salmon bred in even years are completely reproductively isolated from those bred in odd years, even if the two lineages bred in same location. We used time series (mean = 7 years, maximum = 74 years) of paired even- and odd-year populations from 36 rivers spanning over 2000 km to explore parallelism in migration timing, a trait with a strong genetic basis. Migration timing was highly parallel, being determined almost entirely by local environmental differences among rivers. Interestingly, interannual changes in migration timing different somewhat between lineages. Overall, our findings indicate very strong determinism, with only a minor contribution of contingency.

What can be learned by scanning the genome for molecular convergence in wild populations?

Convergent evolution, where independent lineages evolve similar phenotypes in response to similar challenges, can provide valuable insight into how selection operates and the limitations it encounters. However, it has only recently become possible to explore how convergent evolution is reflected at the genomic level. The overlapping outlier approach (OOA), where genome scans of multiple independent lineages are used to find outliers that overlap and therefore identify convergently evolving loci, is becoming popular. Here, we present a quantitative analysis of 34 studies that used this approach across many sampling designs, taxa, and sampling intensities. We found that OOA studies with increased biological sampling power within replicates have increased likelihood of finding overlapping, "convergent" signals of adaptation between them. When identifying convergent loci as overlapping outliers, it is tempting to assume that any false-positive outliers derived from individual scans will fail to overlap across replicates, but this cannot be guaranteed. We highlight how population demographics and genomic context can contribute toward both true convergence and false positives in OOA studies. We finish with an exploration of emerging methods that couple genome scans with phenotype and environmental measures, leveraging added information from genome data to more directly test hypotheses of the likelihood of convergent evolution.

Quantifying the extent of morphological homoplasy: A phylogenetic analysis of 490 characters in Drosophila

Homoplasy is a fundamental phenomenon in evolutionary biology but an appraisal of its extent at the morphological level is still lacking. Here, we analyzed the evolution of 490 morphological characters conceptualized among 56 drosophilid species. We found that two thirds of morphological changes were homoplastic and that the level of homoplasy depended on the stage of development and the type of the organ, with the adult terminalia being the least homoplastic. In spite of its predominance at the character change level, homoplasy accounts for only ∼13% of between species similarities in pairwise comparisons. These results provide empirical insights on the limits of morphological changes and the frequency of recurrent evolution.

The Role of Developmental Integration and Historical Contingency in the Origin and Evolution of Cypriniform Trophic Novelties

While functional morphologists have long studied the evolution of anatomical structures, the origin of morphological novelties has received less attention. When such novelties first originate they must become incorporated into an integrated system to be rendered fully functional. Thus, developmental integration is key at the origin of morphological novelties. However, given enough evolutionary time such integration may be broken, allowing for a division of labor that is facilitated by subsequent decoupling of structures. Cypriniformes represent a diverse group of freshwater fishes characterized by several trophic novelties that include: kinethmoid-mediated premaxillary protrusion, a muscular palatal and post-lingual organ, hypertrophied lower pharyngeal jaws that masticate against the base of the neurocranium, novel pharyngeal musculature controlling movement of the hypertrophied lower pharyngeal jaws, and in a few species an incredibly complex epibranchial organ used to aggregate filtered phytoplankton. Here, we use the wealth of such trophic novelties in different cypriniform fishes to present case studies in which developmental integration allowed for the origin of morphological innovations. As proposed in case studies 1 and 2 trophic innovations may be associated with both morphological and lineage diversification. Alternatively, case studies 3 and 4 represent a situation where ecological niche was expanded but with no concomitant increase in species diversity.

UNVEILing connections between genotype, phenotype, and fitness in natural populations

Understanding the links between genetic variation and fitness in natural populations is a central goal of evolutionary genetics. This monumental task spans the fields of classical and molecular genetics, population genetics, biochemistry, physiology, developmental biology, and ecology. Advances to our molecular and developmental toolkits are facilitating integrative approaches across these traditionally separate fields, providing a more complete picture of the genotype-phenotype map in natural and non-model systems. Here, we summarize research presented at the first annual symposium of the UNVEIL Network, an NSF-funded collaboration between the University of Montana and the University of Nebraska, Lincoln, which took place from the 1st to the 3rd of June, 2018. We discuss how this body of work advances basic evolutionary science, what it implies for our ability to predict evolutionary change, and how it might inform novel conservation strategies.

Fragile DNA contributes to repeated evolution

Sequence features that affect DNA fragility might facilitate fast, repeated evolution by elevating mutation rates at genomic hotspots.

Experimental Evolution of Escherichia coli K-12 in the Presence of Proton Motive Force (PMF) Uncoupler Carbonyl Cyanide m-Chlorophenylhydrazone Selects for Mutations Affecting PMF-Driven Drug Efflux Pumps

Experimental evolution of Escherichia coli K-12 with benzoate, a partial uncoupler of the proton motive force (PMF), selects for mutations that decrease antibiotic resistance. We conducted experimental evolution in the presence of carbonyl cyanide m-chlorophenylhydrazone (CCCP), a strong uncoupler. Cultures were serially diluted daily 1:100 in LBK medium containing 20 to 150 µM CCCP buffered at pH 6.5 or at pH 8.0. After 1,000 generations, the populations tolerated up to 150 µM CCCP. Sequenced isolates had mutations in mprA (emrR), which downregulates the EmrAB-TolC pump that exports CCCP. A mprA::kanR deletion conferred growth at 60 μM CCCP, though not at the higher levels resisted by evolved strains (150 µM). Some mprA mutant strains also had point mutations affecting emrA, but deletion of emrA abolished the CCCP resistance. Thus, CCCP-evolved isolates contained additional adaptations. One isolate lacked emrA or mprA mutations but had mutations in cecR (ybiH), whose product upregulates drug pumps YbhG and YbhFSR, and in gadE, which upregulates the multidrug pump MdtEF. A cecR::kanR deletion conferred partial resistance to CCCP. Other multidrug efflux genes that had mutations included ybhR and acrAB The acrB isolate was sensitive to the AcrAB substrates chloramphenicol and tetracycline. Other mutant genes in CCCP-evolved strains include rng (RNase G) and cyaA (adenylate cyclase). Overall, experimental evolution revealed a CCCP-dependent fitness advantage for mutations increasing CCCP efflux via EmrA and for mutations that may deactivate proton-driven pumps for drugs not present (cecR, gadE, acrAB, and ybhR). These results are consistent with our previous report of drug sensitivity associated with evolved benzoate tolerance.IMPORTANCE The genetic responses of bacteria to depletion of proton motive force (PMF), and their effects on drug resistance, are poorly understood. PMF drives export of many antibiotics, but the energy cost may decrease fitness when antibiotics are absent. Our evolution experiment reveals genetic mechanisms of adaptation to the PMF uncoupler CCCP, including selection for increased CCCP efflux but also against the expression of PMF-driven pumps for drugs not present. The results have implications for our understanding of the gut microbiome, which experiences high levels of organic acids that decrease PMF.

Contingency in the convergent evolution of a regulatory network: Dosage compensation in Drosophila

The repeatability or predictability of evolution is a central question in evolutionary biology and most often addressed in experimental evolution studies. Here, we infer how genetically heterogeneous natural systems acquire the same molecular changes to address how genomic background affects adaptation in natural populations. In particular, we take advantage of independently formed neo-sex chromosomes in Drosophila species that have evolved dosage compensation by co-opting the dosage-compensation male-specific lethal (MSL) complex to study the mutational paths that have led to the acquisition of hundreds of novel binding sites for the MSL complex in different species. This complex recognizes a conserved 21-bp GA-rich sequence motif that is enriched on the X chromosome, and newly formed X chromosomes recruit the MSL complex by de novo acquisition of this binding motif. We identify recently formed sex chromosomes in the D. melanica and D. robusta species groups by genome sequencing and generate genomic occupancy maps of the MSL complex to infer the location of novel binding sites. We find that diverse mutational paths were utilized in each species to evolve hundreds of de novo binding motifs along the neo-X, including expansions of microsatellites and transposable element (TE) insertions. However, the propensity to utilize a particular mutational path differs between independently formed X chromosomes and appears to be contingent on genomic properties of that species, such as simple repeat or TE density. This establishes the "genomic environment" as an important determinant in predicting the outcome of evolutionary adaptations.

An Evolutionary Perspective of Neoplastic Diseases in the Universe

The existence of exoplanets orbiting low mass-stars is one of the most significant discoveries of our time. Especially intriguing to us is the possibility that Earth-sized exoplanets within a habitable zone might harbor life-forms that resemble our own RNA/DNA-based species. We further narrow this theoretical possibility with the following question: if alien life does indeed exist elsewhere, would extraterrestrial life be burdened with earthly diseases? Given that the chemistry of the universe is subject to specific rules, restraints, and predictable outcomes, we argue that cancer-signaling pathways might be programmed into the life cycle of habitable exoplanets. This hypothetical prediction is also based on evolutionary convergence, the repeated emergence of biological similarity that occurs when disparate life-forms adapt to comparable selection pressures. The possibility that mutations and nucleotide base rearrangements that drive cancer growth might be fixed in the chemical hardware of alien life provides us with the opportunity to wonder and consider the origins, evolution, and ubiquity of disease beyond Earth.

EvolvingSTEM: a microbial evolution-in-action curriculum that enhances learning of evolutionary biology and biotechnology

Evolution is a central, unifying theory for all of life science, yet the subject is poorly represented in most secondary-school biology courses, especially in the United States. One challenge to learning evolution is that it is taught as a conceptual, retrospective subject with few tangible outcomes for students. These typical passive learning strategies lead to student disengagement with the material and misunderstanding of evolutionary concepts. To promote greater investment and comprehension, we developed EvolvingSTEM, an inquiry-based laboratory curriculum that demonstrates concepts of natural selection, heredity, and ecological diversity through experimental evolution of a benign bacterium. Students transfer populations of Pseudomonas fluorescens growing on plastic beads, which selects for biofilm formation and mutants with new, conspicuous phenotypes. We introduced our curriculum to four introductory high school biology classes alongside their standard curriculum materials and found that students who learned evolution through EvolvingSTEM scored significantly better on a common assessment targeted to Next Generation Science Standards than students taught only the standard curriculum. This latter group subsequently achieved similar scores once they too completed our curriculum. Our work demonstrates that inquiry-based, hands-on experiences with evolving bacterial populations can greatly enhance student learning of evolutionary concepts.