Non-Darwinian evolution

The Creativity of Natural Selection? Part II: The Synthesis and Since

This is the second of a two-part essay on the history of debates concerning the creativity of natural selection, from Darwin through the evolutionary synthesis and up to the present. In the first part, I focussed on the mid-late nineteenth century to the early twentieth, with special emphasis on early Darwinism and its critics, the self-styled "mutationists." The second part focuses on the evolutionary synthesis and some of its critics, especially the "neutralists" and "neo-mutationists." Like Stephen Gould, I consider the creativity of natural selection to be a key component of what has traditionally counted as "Darwinism." I argue that the creativity of natural selection is best understood in terms of (1) selection initiating evolutionary change, and (2) selection directing evolutionary change, for example by creating the variation that it subsequently acts upon. I consider the respects in which both of these claims sound non-Darwinian, even though they have long been understood by supporters and critics alike to be virtually constitutive of Darwinism.

Speciation Post Synthesis: 1960-2000

Speciation-the origin of new species-has been one of the most active areas of research in evolutionary biology, both during, and since the Modern Synthesis. While the Modern Synthesis certainly shaped research on speciation in significant ways, providing a core framework, and set of categories and methods to work with, the history of work on speciation since the mid-twentieth century is a history of divergence and diversification. This piece traces this divergence, through both theoretical advances, and empirical insights into how different lineages, with different genetics and ecological conditions, are shaped by very different modes of diversification.

A-to-I RNA editing contributes to the persistence of predicted damaging mutations in populations

Adenosine-to-inosine (A-to-I) RNA editing is a very common co-/posttranscriptional modification that can lead to A-to-G changes at the RNA level and compensate for G-to-A genomic changes to a certain extent. It has been shown that each healthy individual can carry dozens of missense variants predicted to be severely deleterious. Why strongly detrimental variants are preserved in a population and not eliminated by negative natural selection remains mostly unclear. Here, we ask if RNA editing correlates with the burden of deleterious A/G polymorphisms in a population. Integrating genome and transcriptome sequencing data from 447 human lymphoblastoid cell lines, we show that nonsynonymous editing activities (prevalence/level) are negatively correlated with the deleteriousness of A-to-G genomic changes and positively correlated with that of G-to-A genomic changes within the population. We find a significantly negative correlation between nonsynonymous editing activities and allele frequency of A within the population. This negative editing-allele frequency correlation is particularly strong when editing sites are located in highly important genes/loci. Examinations of deleterious missense variants from the 1000 Genomes Project further show a significantly higher proportion of rare missense mutations for G-to-A changes than for other types of changes. The proportion for G-to-A changes increases with increasing deleterious effects of the changes. Moreover, the deleteriousness of G-to-A changes is significantly positively correlated with the percentage of editing enzyme binding motifs at the variants. Overall, we show that nonsynonymous editing is associated with the increased burden of G-to-A missense mutations in healthy individuals, expanding RNA editing in pathogenomics studies.

Many roads lead to Rome: Neutral phenotypes in microorganisms

John Bonner pointed out that microorganisms differ in several ways, some of which may reflect neutral phenotypic evolution. For making his case, Bonner referred to interspecies differences and morphological traits. Here we consider intraspecies differences and physiological traits. As a case-study, we examine the production of an extracellular cyclic 3 ' ,5 ' monophosphate phosphodiesterase in the cellular slime mold Dictyostelium discoideum. Temporal profiles of phosphodiesterase activity differ significantly between wild-type strains. From that we argue that the inference drawn initially from studies on a single wild-type, namely that the profile displayed by it pointed to an adaptive role, was mistaken. We generalize the conclusion to suggest that physiological differences exhibited by microorganisms of the same species may, but need not, reflect adaptations to different environments. Rather, the differences could be related to the fact that microorganisms live in groups whose composition can vary between homogeneous (clonal) and heterogeneous (polyclonal). More than one physiological profile is consistent with the normal development of the group in a given environment; the alternatives are neutral. When studying microbial physiology and behavior, it is expected that the observations are made on a clonal population; genetic (and so phenotypic) heterogeneity is carefully guarded against. As the example from D. discoideum shows, an unintended consequence of overlooking phenotypic heterogeneity is that one can fall into the trap of accepting a seemingly plausible, but possibly erroneous, adaptive explanation for a "normal" wild-type phenotype.

The Isolation of New Pore-Forming Toxins from the Sea Anemone Actinia fragacea Provides Insights into the Mechanisms of Actinoporin Evolution

Random mutations and selective pressure drive protein adaptation to the changing demands of the environment. As a consequence, nature favors the evolution of protein diversity. A group of proteins subject to exceptional environmental stress and known for their widespread diversity are the pore-forming hemolytic proteins from sea anemones, known as actinoporins. In this study, we identified and isolated new isoforms of actinoporins from the sea anemone Actinia fragacea (fragaceatoxins). We characterized their hemolytic activity, examined their stability and structure, and performed a comparative analysis of their primary sequence. Sequence alignment reveals that most of the variability among actinoporins is associated with non-functional residues. The differences in the thermal behavior among fragaceatoxins suggest that these variability sites contribute to changes in protein stability. In addition, the protein-protein interaction region showed a very high degree of identity (92%) within fragaceatoxins, but only 25% among all actinoporins examined, suggesting some degree of specificity at the species level. Our findings support the mechanism of evolutionary adaptation in actinoporins and reflect common pathways conducive to protein variability.

Understanding the Genetic Code

The universal triple-nucleotide genetic code is often viewed as a given, randomly selected through evolution. However, as summarized in this article, many observations and deductions within structural and thermodynamic frameworks help to explain the forces that must have shaped the code during the early evolution of life on Earth.

Genetic Diversity and Gene Family Expansions in Members of the Genus Entamoeba

Amoebiasis is the third-most common cause of mortality worldwide from a parasitic disease. Although the primary etiological agent of amoebiasis is the obligate human parasite Entamoeba histolytica, other members of the genus Entamoeba can infect humans and may be pathogenic. Here, we present the first annotated reference genome for Entamoeba moshkovskii, a species that has been associated with human infections, and compare the genomes of E. moshkovskii, E. histolytica, the human commensal Entamoeba dispar, and the nonhuman pathogen Entamoeba invadens. Gene clustering and phylogenetic analyses show differences in expansion and contraction of families of proteins associated with host or bacterial interactions. They intimate the importance to parasitic Entamoeba species of surface-bound proteins involved in adhesion to extracellular membranes, such as the Gal/GalNAc lectin and members of the BspA and Ariel1 families. Furthermore, E. dispar is the only one of the four species to lack a functional copy of the key virulence factor cysteine protease CP-A5, whereas the gene's presence in E. moshkovskii is consistent with the species' potentially pathogenic nature. Entamoeba moshkovskii was found to be more diverse than E. histolytica across all sequence classes. The former is ∼200 times more diverse than latter, with the four E. moshkovskii strains tested having a most recent common ancestor nearly 500 times more ancient than the tested E. histolytica strains. A four-haplotype test indicates that these E. moshkovskii strains are not the same species and should be regarded as a species complex.

Protein painting, an optimized MS-based technique, reveals functionally relevant interfaces of the PD-1/PD-L1 complex and the YAP2/ZO-1 complex

Protein-protein interactions lie at the heart of many biological processes and therefore represent promising drug targets. Despite this opportunity, identification of protein-protein interfaces remains challenging. We have previously developed a method that relies on coating protein surfaces with small-molecule dyes to discriminate between solvent-accessible protein surfaces and hidden interface regions. Dye-bound, solvent-accessible protein regions resist trypsin digestion, whereas hidden interface regions are revealed by denaturation and sequenced by MS. The small-molecule dyes bind promiscuously and with high affinity, but their binding mechanism is unknown. Here, we report on the optimization of a novel dye probe used in protein painting, Fast Blue B + naphthionic acid, and show that its affinity for proteins strongly depends on hydrophobic moieties that we call here "hydrophobic clamps." We demonstrate the utility of this probe by sequencing the protein-protein interaction regions between the Hippo pathway protein Yes-associated protein 2 (YAP2) and tight junction protein 1 (TJP1 or ZO-1), uncovering interactions via the known binding domain as well as ZO-1's MAGUK domain and YAP's N-terminal proline-rich domain. Additionally, we demonstrate how residues predicted by protein painting are present exclusively in the complex interface and how these residues may guide the development of peptide inhibitors using a case study of programmed cell death protein 1 (PD-1) and programmed cell death 1 ligand 1 (PD-L1). Inhibitors designed around the PD-1/PD-L1 interface regions identified via protein painting effectively disrupted complex formation, with the most potent inhibitor having an IC₅₀ of 5 μm.

Early in life effects and heredity: reconciling neo-Darwinism with neo-Lamarckism under the banner of the inclusive evolutionary synthesis

Recent discoveries show that early in life effects often have long-lasting influences, sometimes even spanning several generations. Such intergenerational effects of early life events appear not easily reconcilable with strict genetic inheritance. However, an integrative evolutionary medicine of early life effects needs a sound view of inheritance in development and evolution. Here, we show how to articulate the gene-centred and non-gene-centred visions of inheritance. We first recall the coexistence of two gene concepts in scientific discussions, a statistical one (focused on patterns of parent-offspring resemblance, and implicitly including non-DNA-sequence-based resemblance), and a molecular one (based on the DNA sequence). We then show how all the different mechanisms of inheritance recently discovered can be integrated into an inclusive theory of evolution where different mechanisms would enable adaptation to changing environments at different timescales. One surprising consequence of this integrative vision of inheritance is that early in life effects start much earlier than fertilization. This article is part of the theme issue 'Developing differences: early-life effects and evolutionary medicine'.

Neutral Theory and Phenotypic Evolution

Although the neutral theory of molecular evolution was proposed to explain DNA and protein sequence evolution, in principle it could also explain phenotypic evolution. Nevertheless, overall, phenotypes should be less likely than genotypes to evolve neutrally. I propose that, when phenotypic traits are stratified according to a hierarchy of biological organization, the fraction of evolutionary changes in phenotype that are adaptive rises with the phenotypic level considered. Consistently, molecular traits are frequently found to evolve neutrally whereas a large, random set of organismal traits were recently reported to vary largely adaptively. Many more studies of unbiased samples of phenotypic traits are needed to test the general validity of this hypothesis.

The Neutral Theory in Light of Natural Selection

In this perspective, we evaluate the explanatory power of the neutral theory of molecular evolution, 50 years after its introduction by Kimura. We argue that the neutral theory was supported by unreliable theoretical and empirical evidence from the beginning, and that in light of modern, genome-scale data, we can firmly reject its universality. The ubiquity of adaptive variation both within and between species means that a more comprehensive theory of molecular evolution must be sought.

Concise Synthesis of Potassium Acyltrifluoroborates from Aldehydes through Copper(I)-Catalyzed Borylation/Oxidation

Potassium acyltrifluoroborates (KATs) were prepared through copper(I)-catalyzed borylation of aldehydes and subsequent oxidation. This synthetic route is characterized by the wide range of aldehydes accessible, favorable step economy, mild reaction conditions, and tolerance of various functional groups, and it enables the facile generation of a range of KATs, for example, bearing halide, sulfide, acetal, or ester moieties. Moreover, this method was applied to the three-step synthesis of various α-amino acid analogues that bear a KAT moiety on the C-terminus by using naturally occurring amino acids as the starting material.

Thermophilic Adaptation in Prokaryotes Is Constrained by Metabolic Costs of Proteostasis

Prokaryotes evolved to thrive in an extremely diverse set of habitats, and their proteomes bear signatures of environmental conditions. Although correlations between amino acid usage and environmental temperature are well-documented, understanding of the mechanisms of thermal adaptation remains incomplete. Here, we couple the energetic costs of protein folding and protein homeostasis to build a microscopic model explaining both the overall amino acid composition and its temperature trends. Low biosynthesis costs lead to low diversity of physical interactions between amino acid residues, which in turn makes proteins less stable and drives up chaperone activity to maintain appropriate levels of folded, functional proteins. Assuming that the cost of chaperone activity is proportional to the fraction of unfolded client proteins, we simulated thermal adaptation of model proteins subject to minimization of the total cost of amino acid synthesis and chaperone activity. For the first time, we predicted both the proteome-average amino acid abundances and their temperature trends simultaneously, and found strong correlations between model predictions and 402 genomes of bacteria and archaea. The energetic constraint on protein evolution is more apparent in highly expressed proteins, selected by codon adaptation index. We found that in bacteria, highly expressed proteins are similar in composition to thermophilic ones, whereas in archaea no correlation between predicted expression level and thermostability was observed. At the same time, thermal adaptations of highly expressed proteins in bacteria and archaea are nearly identical, suggesting that universal energetic constraints prevail over the phylogenetic differences between these domains of life.

Model selection may not be a mandatory step for phylogeny reconstruction

Determining the most suitable model for phylogeny reconstruction constitutes a fundamental step in numerous evolutionary studies. Over the years, various criteria for model selection have been proposed, leading to debate over which criterion is preferable. However, the necessity of this procedure has not been questioned to date. Here, we demonstrate that although incongruency regarding the selected model is frequent over empirical and simulated data, all criteria lead to very similar inferences. When topologies and ancestral sequence reconstruction are the desired output, choosing one criterion over another is not crucial. Moreover, skipping model selection and using instead the most parameter-rich model, GTR+I+G, leads to similar inferences, thus rendering this time-consuming step nonessential, at least under current strategies of model selection.

Characterization of the mitochondrial genome of an ancient amphipod Halice sp. MT-2017 (Pardaliscidae) from 10,908 m in the Mariana Trench

Small amphipods (Halice sp. MT-2017) with body length <1 cm were collected from the Challenger Deep (~10,920 m below sea level). The divergence time of their lineage was approximately 109 Mya, making this group ancient compared to others under study. The mitochondrial genome of Halice sp. shared the usual gene components of metazoans, comprising 13 protein coding genes (PCGs), 22 transfer RNAs (tRNAs), and 2 ribosomal RNAs (rRNAs). The arrangement of these genes, however, differed greatly from that of other amphipods. Of the 15 genes that were rearranged with respect to the pancrustacean gene pattern, 12 genes (2 PCGs, 2 rRNAs, and 8 tRNAs) were both translocated and strand-reversed. In contrast, the mitochondrial genomes in other amphipods never show so many reordered genes, and in most instances, only tRNAs were involved in strand-reversion-coupled translocation. Other characteristics, including reversed strand nucleotide composition bias, relatively higher composition of non-polar amino acids, and lower evolutionary rate, were also identified. Interestingly, the latter two features were shared with another hadal amphipod, Hirondellea gigas, suggesting their possible associations with the adaptation to deep-sea extreme habitats. Overall, our data provided a useful resource for future studies on the evolutionary and adaptive mechanisms of hadal faunas.

Social evolution under demographic stochasticity

How social traits such as altruism and spite evolve remains an open question in evolutionary biology. One factor thought to be potentially important is demographic stochasticity. Here we provide a general theoretical analysis of the role of demographic stochasticity in social evolution. We show that the evolutionary impact of stochasticity depends on how the social action alters the recipient’s life cycle. If the action alters the recipient’s death rate, then demographic stochasticity always favours altruism and disfavours spite. On the other hand, if the action alters the recipient’s birth rate, then stochasticity can either favour or disfavour both altruism and spite depending on the ratio of the rate of population turnover to the population size. Finally, we also show that this ratio is critical to determining if demographic stochasticity can reverse the direction of selection upon social traits. Our analysis thus provides a general understanding of the role of demographic stochasticity in social evolution.

Explaining the evolution of social traits such as altruism and spite remains a key outstanding problem in evolutionary biology. Here we develop a simple theory for the effect of demographic stochasticity (random variation in an individual’s birth and death rates) on the evolution of social traits. Our results provide a clear set of predictions: whether a social trait is favoured or disfavoured is determined by how the social action alters the recipient’s life cycle. If the social action alters the recipient’s death rate, then altruism is favoured and spite disfavoured. If instead the social action alters the recipient’s birth rate, then both altruism and spite can be either favoured or disfavoured—the precise outcome depends upon the ratio of the population turnover rate to the population size.

Multivariate morphometrics, quantitative genetics, and neutral theory: Developing a "modern synthesis" for primate evolutionary morphology

Anthropologists are increasingly turning to explicit model-bound evolutionary approaches for understanding the morphological diversification of humans and other primate lineages. Such evolutionary morphological analyses rely on three interconnected conceptual frameworks; multivariate morphometrics for quantifying similarity and differences among taxa, quantitative genetics for modeling the inheritance and evolution of morphology, and neutral theory for assessing the likelihood that taxon diversification is due to stochastic processes such as genetic drift. Importantly, neutral theory provides a framework for testing more parsimonious explanations for observed morphological differences before considering more complex adaptive scenarios. However, the consistency with which these concepts are applied varies considerably, which mirrors some of the theoretical obstacles faced during the "modern synthesis" of classical population genetics in the early 20th century. Here, each framework is reviewed and some potential stumbling blocks to the full conceptual integration of multivariate morphometrics, quantitative genetics, and neutral theory are considered.

Genetic Code Error Minimization as a Non-Adaptive But Beneficial Trait

In a recent Letter, Di Giulio questions the use of the term 'neutral' when describing the process by which error minimization may have arisen as a side-product of genetic code expansion, resulting from the addition of similar amino acids to similar codons (Di Giulio, in J Mol Evol 86(9):593-597, 2018). However, I point out that in this scenario error minimization is non-adaptive, and so 'neutral' is an appropriate term to describe its imperviousness to direct selection. Error minimization is a form of mutational robustness, and so commonly viewed as beneficial. This in turn implies that not all beneficial traits may be adaptations generated by direct selection for that trait.

The importance of the Neutral Theory in 1968 and 50 years on: A response to Kern and Hahn 2018

A recent article reassessing the Neutral Theory of Molecular Evolution claims that it is no longer as important as is widely believed. The authors argue that "the neutral theory was supported by unreliable theoretical and empirical evidence from the beginning, and that in light of modern, genome-scale data, we can firmly reject its universality." Claiming that "the neutral theory has been overwhelmingly rejected," they propose instead that natural selection is the major force shaping both between-species divergence and within-species variation. Although this is probably a minority view, it is important to evaluate such claims carefully in the context of current knowledge, as inaccuracies can sometimes morph into an accepted narrative for those not familiar with the underlying science. We here critically examine and ultimately reject Kern and Hahn's arguments and assessment, and instead propose that it is now abundantly clear that the foundational ideas presented five decades ago by Kimura and Ohta are indeed correct.

pH-Induced Folding of the Caspase-Cleaved Par-4 Tumor Suppressor: Evidence of Structure Outside of the Coiled Coil Domain

Prostate apoptosis response-4 (Par-4) is a 38 kDa largely intrinsically disordered tumor suppressor protein that functions in cancer cell apoptosis. Par-4 down-regulation is often observed in cancer while up-regulation is characteristic of neurodegenerative conditions such as Alzheimer’s disease. Cleavage of Par-4 by caspase-3 activates tumor suppression via formation of an approximately 25 kDa fragment (cl-Par-4) that enters the nucleus and inhibits Bcl-2 and NF-ƙB, which function in pro-survival pathways. Here, we have investigated the structure of cl-Par-4 using biophysical techniques including circular dichroism (CD) spectroscopy, dynamic light scattering (DLS), and intrinsic tyrosine fluorescence. The results demonstrate pH-dependent folding of cl-Par-4, with high disorder and aggregation at neutral pH, but a largely folded, non-aggregated conformation at acidic pH.

Bioinformatics and Translation Elongation

Codon usage depends on mutation bias, tRNA-mediated selection, and the need for high efficiency and accuracy in translation. One codon in a synonymous codon family is often strongly over-used, especially in highly expressed genes, which often leads to a high dN/dS ratio because dS is very small. Many different codon usage indices have been proposed to measure codon usage and codon adaptation. Sense codon could be misread by release factors and stop codons misread by tRNAs, which also contribute to codon usage in rare cases. This chapter outlines the conceptual framework on codon evolution, illustrates codon-specific and gene-specific codon usage indices, and presents their applications. A new index for codon adaptation that accounts for background mutation bias (Index of Translation Elongation) is presented and contrasted with codon adaptation index (CAI) which does not consider background mutation bias. They are used to re-analyze data from a recent paper claiming that translation elongation efficiency matters little in protein production. The reanalysis disproves the claim.

Invoking adaptation to decipher the genetic legacy of past climate change

Persistence of natural populations during periods of climate change is likely to depend on migration (range shifts) or adaptation. These responses were traditionally considered discrete processes and conceptually divided into the realms of ecology and evolution. In a milestone paper, Davis and Shaw (2001) Science 292:673 argued that the interplay of adaptation and migration was central to biotic responses to Quaternary climate, but since then there has been no synthesis of efforts made to set up this research program. Here we review some of the salient findings from molecular genetic studies assessing ecological and evolutionary responses to Quaternary climate change. These studies have revolutionized our understanding of population processes associated with past species migration. However, knowledge remains limited about the role of natural selection for local adaptation of populations to Quaternary environmental fluctuations and associated range shifts, and for the footprints this might have left on extant populations. Next-generation sequencing technologies, high-resolution paleoclimate analyses, and advances in population genetic theory offer an unprecedented opportunity to test hypotheses about adaptation through time. Recent population genomics studies have greatly improved our understanding of the role of contemporary adaptation to local environments in shaping spatial patterns of genetic diversity across modern-day landscapes. Advances in this burgeoning field provide important conceptual and methodological bases to decipher the historical role of natural selection and assess adaptation to past environmental variation. We suggest that a process called "temporal conditional neutrality" has taken place: some alleles favored in glacial environments become selectively neutral in modern-day conditions, whereas some alleles that had been neutral during glacial periods become under selection in modern environments. Building on this view, we present a new integrative framework for addressing the interplay of demographic and adaptive evolutionary responses to Quaternary climate dynamics, the research agenda initially envisioned by Davis and Shaw (2001) Science 292:673.

The Experimental Study of Bacterial Evolution and Its Implications for the Modern Synthesis of Evolutionary Biology

Since the 1940s, microbiologists, biochemists and population geneticists have experimented with the genetic mechanisms of microorganisms in order to investigate evolutionary processes. These evolutionary studies of bacteria and other microorganisms gained some recognition from the standard-bearers of the modern synthesis of evolutionary biology, especially Theodosius Dobzhansky and Ledyard Stebbins. A further period of post-synthesis bacterial evolutionary research occurred between the 1950s and 1980s. These experimental analyses focused on the evolution of population and genetic structure, the adaptive gain of new functions, and the evolutionary consequences of competition dynamics. This large body of research aimed to make evolutionary theory testable and predictive, by giving it mechanistic underpinnings. Although evolutionary microbiologists promoted bacterial experiments as methodologically advantageous and a source of general insight into evolution, they also acknowledged the biological differences of bacteria. My historical overview concludes with reflections on what bacterial evolutionary research achieved in this period, and its implications for the still-developing modern synthesis.

Neutral Theory, Disease Mutations, and Personal Exomes

Genetic differences between species and within populations are two sides of the same coin under the neutral theory of molecular evolution. This theory posits that a vast majority of evolutionary substitutions, which appear as differences between species, are (nearly) neutral, that is, these substitutions are permitted without a significantly adverse impact on a species' survival. We refer to them as evolutionarily permissible (ePerm) variation. Evolutionary permissibility of any possible variant can be inferred from multispecies sequence alignments by applying sophisticated statistical methods to the evolutionary tree of species. Here, we explore the evolutionary permissibility of amino acid variants associated with genetic diseases and those observed in personal exomes. Consistent with the predictions of the neutral theory, disease associated amino acid variants are rarely ePerm, much more biochemically radical, and found predominantly at more conserved positions than their non-disease counterparts. Only 10% of amino acid mutations are ePerm, but these variants rise to become two-thirds of all substitutions in the human lineage (a 6-fold enrichment). In contrast, only a minority of the variants in a personal exome are ePerm, a seemingly counterintuitive pattern that results from a combination of mutational and evolutionary processes that are, in fact, broadly consistent with the neutral theory. Evolutionarily forbidden variants outnumber detrimental variants in individual exomes and may play an underappreciated role in protecting against disease. We discuss these observations and conclude that the long-term evolutionary history of species can illuminate functional biomedical properties of variation present in personal exomes.

Cyclic and Lasso Peptides: Sequence Determination, Topology Analysis, and Rotaxane Formation

A broadly applicable chemical cleavage methodology to facilitate MS/MS sequencing was developed for macrocyclic and lasso peptides, which hold promise as exciting new therapeutics. Existing methods such as Edman degradation, CNBr cleavage, and enzymatic digestion are either limited in scope or completely fail in cleavage of constrained nonribosomal peptides. Importantly, the new method was utilized for synthesizing a unique peptide-based rotaxane (both cyclic and threaded) from the lasso peptide, benenodin-1 ΔC5.

Proteomic and Metabolomic Characterization of a Mammalian Cellular Transition from Quiescence to Proliferation

There exist similarities and differences in metabolism and physiology between normal proliferative cells and tumor cells. Once a cell enters the cell cycle, metabolic machinery is engaged to facilitate various processes. The kinetics and regulation of these metabolic changes have not been properly evaluated. To correlate the orchestration of these processes with the cell cycle, we analyzed the transition from quiescence to proliferation of a non-malignant murine pro-B lymphocyte cell line in response to IL-3. Using multiplex mass-spectrometry-based proteomics, we show that the transition to proliferation shares features generally attributed to cancer cells: upregulation of glycolysis, lipid metabolism, amino-acid synthesis, and nucleotide synthesis and downregulation of oxidative phosphorylation and the urea cycle. Furthermore, metabolomic profiling of this transition reveals similarities to cancer-related metabolic pathways. In particular, we find that methionine is consumed at a higher rate than that of other essential amino acids, with a potential link to maintenance of the epigenome.

Real-Time Intrinsic Fluorescence Visualization and Sizing of Proteins and Protein Complexes in Microfluidic Devices

Optical detection has become a convenient and scalable approach to read out information from microfluidic systems. For the study of many key biomolecules, however, including peptides and proteins, which have low fluorescence emission efficiencies at visible wavelengths, this approach typically requires labeling of the species of interest with extrinsic fluorophores to enhance the optical signal obtained - a process which can be time-consuming, requires purification steps, and has the propensity to perturb the behavior of the systems under study due to interactions between the labels and the analyte molecules. As such, the exploitation of the intrinsic fluorescence of protein molecules in the UV range of the electromagnetic spectrum is an attractive path to allow the study of unlabeled proteins. However, direct visualization using 280 nm excitation in microfluidic devices has to date commonly required the use of coherent sources with frequency multipliers and devices fabricated out of materials that are incompatible with soft lithography techniques. Here, we have developed a simple, robust, and cost-effective 280 nm LED platform that allows real-time visualization of intrinsic fluorescence from both unlabeled proteins and protein complexes in polydimethylsiloxane microfluidic channels fabricated through soft lithography. Using this platform, we demonstrate intrinsic fluorescence visualization of proteins at nanomolar concentrations on chip and combine visualization with micron-scale diffusional sizing to measure the hydrodynamic radii of individual proteins and protein complexes under their native conditions in solution in a label-free manner.

Standard Codon Substitution Models Overestimate Purifying Selection for Nonstationary Data

Estimation of natural selection on protein-coding sequences is a key comparative genomics approach for de novo prediction of lineage-specific adaptations. Selective pressure is measured on a per-gene basis by comparing the rate of nonsynonymous substitutions to the rate of synonymous substitutions. All published codon substitution models have been time-reversible and thus assume that sequence composition does not change over time. We previously demonstrated that if time-reversible DNA substitution models are applied in the presence of changing sequence composition, the number of substitutions is systematically biased towards overestimation. We extend these findings to the case of codon substitution models and further demonstrate that the ratio of nonsynonymous to synonymous rates of substitution tends to be underestimated over three data sets of mammals, vertebrates, and insects. Our basis for comparison is a nonstationary codon substitution model that allows sequence composition to change. Goodness-of-fit results demonstrate that our new model tends to fit the data better. Direct measurement of nonstationarity shows that bias in estimates of natural selection and genetic distance increases with the degree of violation of the stationarity assumption. Additionally, inferences drawn under time-reversible models are systematically affected by compositional divergence. As genomic sequences accumulate at an accelerating rate, the importance of accurate de novo estimation of natural selection increases. Our results establish that our new model provides a more robust perspective on this fundamental quantity.

Neutral Evolution

Neutral evolution is the default process of genomic changes. This is because our world is finite, and the randomness, indispensable for neutral evolution, is important when we consider the history of a finite world. The random nature of DNA propagation is discussed using branching process, coalescent process, Markov process, and diffusion process. Expected evolutionary patterns under neutrality are then discussed on fixation probability, rate of evolution, and amount of DNA variation kept in population. We then discuss various features of neutral evolution starting from evolutionary rates, synonymous and nonsynonymous substitutions, junk DNA, and pseudogenes.

Key issues review: evolution on rugged adaptive landscapes

Adaptive landscapes represent a mapping between genotype and fitness. Rugged adaptive landscapes contain two or more adaptive peaks: allele combinations with higher fitness than any of their neighbors in the genetic space. How do populations evolve on such rugged landscapes? Evolutionary biologists have struggled with this question since it was first introduced in the 1930s by Sewall Wright. Discoveries in the fields of genetics and biochemistry inspired various mathematical models of adaptive landscapes. The development of landscape models led to numerous theoretical studies analyzing evolution on rugged landscapes under different biological conditions. The large body of theoretical work suggests that adaptive landscapes are major determinants of the progress and outcome of evolutionary processes. Recent technological advances in molecular biology and microbiology allow experimenters to measure adaptive values of large sets of allele combinations and construct empirical adaptive landscapes for the first time. Such empirical landscapes have already been generated in bacteria, yeast, viruses, and fungi, and are contributing to new insights about evolution on adaptive landscapes. In this Key Issues Review we will: (i) introduce the concept of adaptive landscapes; (ii) review the major theoretical studies of evolution on rugged landscapes; (iii) review some of the recently obtained empirical adaptive landscapes; (iv) discuss recent mathematical and statistical analyses motivated by empirical adaptive landscapes, as well as provide the reader with instructions and source code to implement simulations of evolution on adaptive landscapes; and (v) discuss possible future directions for this exciting field.

Synonymous Codon Usage Controls Various Molecular Aspects

Synonymous sites are generally considered to be functionally neutral. However, there are recent contradictory findings suggesting that synonymous alleles might have functional roles in various molecular aspects. For instance, a recent study demonstrated that synonymous single nucleotide polymorphisms have a similar effect size as nonsynonymous single nucleotide polymorphisms in human disease association studies. Researchers have recognized synonymous codon usage bias (SCUB) in the genomes of almost all species and have investigated whether SCUB is due to random nucleotide compositional bias or to natural selection of any functional exposure generated by synonymous mutations. One of the most prominent observations on the non-neutrality of synonymous codons is the correlation between SCUB and levels of gene expression, such that highly expressed genes tend to have a higher preference toward so-called optimal codons than lowly expressed genes. In relation, it is known that amounts of cognate tRNAs that bind to optimal codons are significantly higher than the amounts of cognate tRNAs that bind to non-optimal codons in genomes. In the present paper, we review various functions that synonymous codons might have other than regulating expression levels.

Genome variation and conserved regulation identify genomic regions responsible for strain specific phenotypes in rat

Background

The genomes of laboratory rat strains are characterised by a mosaic haplotype structure caused by their unique breeding history. These mosaic haplotypes have been recently mapped by extensive sequencing of key strains. Comparison of genomic variation between two closely related rat strains with different phenotypes has been proposed as an effective strategy for the discovery of candidate strain-specific regions involved in phenotypic differences. We developed a method to prioritise strain-specific haplotypes by integrating genomic variation and genomic regulatory data predicted to be involved in specific phenotypes. Specifically, we aimed to identify genomic regions associated with Metabolic Syndrome (MetS), a disorder of energy utilization and storage affecting several organ systems.

Results

We compared two Lyon rat strains, Lyon Hypertensive (LH) which is susceptible to MetS, and Lyon Low pressure (LL), which is susceptible to obesity as an intermediate MetS phenotype, with a third strain (Lyon Normotensive, LN) that is resistant to both MetS and obesity. Applying a novel metric, we ranked the identified strain-specific haplotypes using evolutionary conservation of the occupancy three liver-specific transcription factors (HNF4A, CEBPA, and FOXA1) in five rodents including rat. Consideration of regulatory information effectively identified regions with liver-associated genes and rat orthologues of human GWAS variants related to obesity and metabolic traits. We attempted to find possible causative variants and compared them with the candidate genes proposed by previous studies. In strain-specific regions with conserved regulation, we found a significant enrichment for published evidence to obesity—one of the metabolic symptoms shown by the Lyon strains—amongst the genes assigned to promoters with strain-specific variation.

Conclusions

Our results show that the use of functional regulatory conservation is a potentially effective approach to select strain-specific genomic regions associated with phenotypic differences among Lyon rats and could be extended to other systems.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-4351-9) contains supplementary material, which is available to authorized users.

Mineral surface chemistry control for origin of prebiotic peptides

Some seventy years ago, John Desmond Bernal proposed a role for clays in the origin of life. While much research has since been dedicated to the study of silicate clays, layered double hydroxides, believed to be common on the early Earth, have received only limited attention. Here we examine the role that layered hydroxides could have played in prebiotic peptide formation. We demonstrate how these minerals can concentrate, align and act as adsorption templates for amino acids, and during wetting—drying cycles, promote peptide bond formation. This enables us to propose a testable mechanism for the growth of peptides at layered double hydroxide interfaces in an early Earth environment. Our results provide insights into the potential role of mineral surfaces in mimicking aspects of biochemical reaction pathways.

Clay is thought to have played a part in the origin of life. Here, the authors show that layered double hydroxides, a type of clay little studied despite its presumed prevalence on the early Earth, can facilitate the formation of small proteins.

Most m6A RNA Modifications in Protein-Coding Regions Are Evolutionarily Unconserved and Likely Nonfunctional

Methylation of the adenosine base at the nitrogen-6 position (m6A) is the most prevalent internal posttranscriptional modification of mRNAs in many eukaryotes. Despite the rapid progress in the transcriptome-wide mapping of m6As, identification of proteins responsible for writing, reading, and erasing m6As, and elucidation of m6A functions in splicing, RNA stability, translation, and other processes, it is unknown whether most observed m6A modifications are functional. To address this question, we respectively analyze the evolutionary conservation of yeast and human m6As in protein-coding regions. Relative to comparable unmethylated As, m6As are overall no more conserved in yeasts and only slightly more conserved in mammals. Furthermore, yeast m6As and comparable unmethylated As have no significant difference in single nucleotide polymorphism (SNP) density or SNP site frequency spectrum. The same is true in human. The methylation status of a gene, not necessarily the specific sites methylated in the gene, is subject to purifying selection for no more than ∼20% of m6A-modified genes. These observations suggest that most m6A modifications in protein-coding regions are nonfunctional and nonadaptive, probably resulting from off-target activities of m6A methyltransferases. In addition, our reanalysis invalidates the recent claim of positive selection for newly acquired m6A modifications in human evolution. Regarding the small number of evolutionarily conserved m6As, evidence suggests that a large proportion of them are likely functional; they should be prioritized in future functional characterizations of m6As. Together, these findings have important implications for understanding the biological significance of m6A and other posttranscriptional modifications.

Stochasticity in cultural evolution: a revolution yet to happen

Over the last 40 years or so, there has been an explosion of cultural evolution research in anthropology and archaeology. In each discipline, cultural evolutionists investigate how interactions between individuals translate into group level patterns, with the aim of explaining the diachronic dynamics and diversity of cultural traits. However, while much attention has been given to deterministic processes (e.g. cultural transmission biases), we contend that current evolutionary accounts of cultural change are limited because they do not adopt a systematic stochastic approach (i.e. accounting for the role of chance). First, we show that, in contrast with the intense debates in ecology and population genetics, the importance of stochasticity in evolutionary processes has generated little discussion in the sciences of cultural evolution to date. Second, we speculate on the reasons, both ideological and methodological, why that should be so. Third, we highlight the inadequacy of genetically-inspired stochastic models in the context of cultural evolution modelling, and ask which fundamental stochastic processes might be more relevant to take up. We conclude that the field of cultural evolution would benefit from a stochastic revolution. For that to occur, stochastic models ought to be developed specifically for cultural data and not through a copy-pasting of neutral models from population genetics or ecology.

Molecular Evolution of Herpes Simplex Virus 2 Complete Genomes: Comparison between Primary and Recurrent Infections

Herpes simplex virus 1 (HSV-1) and HSV-2 are large, double-stranded DNA viruses that cause lifelong persistent infections characterized by periods of quiescence and recurrent disease. How HSV evolves within an infected individual experiencing multiple episodes of recurrent disease over time is not known. We determined the genome sequences of viruses isolated from two subjects in the Herpevac Trial for Women who experienced primary HSV-2 genital disease and compared them with sequences of viruses isolated from the subsequent fifth or sixth episode of recurrent disease in the same individuals. Each of the HSV-2 genome sequences was initially obtained using next-generation sequencing and completed with Sanger sequencing. Polymorphisms over the entire genomes were mapped, and amino acid variants resulting from nonsynonymous changes were analyzed based on the secondary and tertiary structures of a previously crystallized protein. A phylogenetic reconstruction was used to assess relationships among the four HSV-2 samples, other North American sequences, and reference sequences. Little genetic drift was detected in viruses shed by the same subjects following repeated reactivation events, suggesting strong selective pressure on the viral genome to maintain sequence fidelity during reactivations from its latent state within an individual host. Our results also demonstrate that some primary HSV-2 isolates from North America more closely resemble the HG52 laboratory strain from Scotland than the low-passage-number clinical isolate SD90e from South Africa or laboratory strain 333. Thus, one of the sequences reported here would be a logical choice as a reference strain for inclusion in future studies of North American HSV-2 isolates.IMPORTANCE The extent to which the HSV-2 genome evolves during multiple episodes of reactivation from its latent state within an infected individual is not known. We used next-generation sequencing techniques to determine whole-genome sequences of four viral samples from two subjects in the Herpevac Trial. The sequence of each subject's well-documented primary isolate was compared with the sequence of the isolate from their fifth or sixth episode of recurrent disease. Only 19 genetic polymorphisms unique to the primary or recurrent isolate were identified, 10 in subject A and 9 in subject B. These observations indicate remarkable genetic conservation between primary and recurrent episodes of HSV-2 infection and imply that strong selection pressures exist to maintain the fidelity of the viral genome during repeated reactivations from its latent state. The genome conservation observed also has implications for the potential success of a therapeutic vaccine.

Functionality and the Evolution of Marginal Stability in Proteins: Inferences from Lattice Simulations

It has been known for some time that many proteins are marginally stable. This has inspired several explanations. Having noted that the functionality of many enzymes is correlated with subunit motion, flexibility, or general disorder, some have suggested that marginally stable proteins should have an evolutionary advantage over proteins of differing stability. Others have suggested that stability and functionality are contradictory qualities, and that selection for both criteria results in marginally stable proteins, optimised to satisfy the competing design pressures. While these explanations are plausible, recent research simulating the evolution of model proteins has shown that selection for stability, ignoring any aspects of functionality, can result in marginally stable proteins because of the underlying makeup of protein sequence-space. We extend this research by simulating the evolution of proteins, using a computational protein model that equates functionality with binding and catalysis. In the model, marginal stability is not required for ligand-binding functionality and we observe no competing design pressures. The resulting proteins are marginally stable, again demonstrating that neutral evolution is sufficient for explaining marginal stability in observed proteins.

Evolutionary biology today and the call for an extended synthesis

Evolutionary theory has been extended almost continually since the evolutionary synthesis (ES), but except for the much greater importance afforded genetic drift, the principal tenets of the ES have been strongly supported. Adaptations are attributable to the sorting of genetic variation by natural selection, which remains the only known cause of increase in fitness. Mutations are not adaptively directed, but as principal authors of the ES recognized, the material (structural) bases of biochemistry and development affect the variety of phenotypic variations that arise by mutation and recombination. Against this historical background, I analyse major propositions in the movement for an 'extended evolutionary synthesis'. 'Niche construction' is a new label for a wide variety of well-known phenomena, many of which have been extensively studied, but (as with every topic in evolutionary biology) some aspects may have been understudied. There is no reason to consider it a neglected 'process' of evolution. The proposition that phenotypic plasticity may engender new adaptive phenotypes that are later genetically assimilated or accommodated is theoretically plausible; it may be most likely when the new phenotype is not truly novel, but is instead a slight extension of a reaction norm already shaped by natural selection in similar environments. However, evolution in new environments often compensates for maladaptive plastic phenotypic responses. The union of population genetic theory with mechanistic understanding of developmental processes enables more complete understanding by joining ultimate and proximate causation; but the latter does not replace or invalidate the former. Newly discovered molecular phenomena have been easily accommodated in the past by elaborating orthodox evolutionary theory, and it appears that the same holds today for phenomena such as epigenetic inheritance. In several of these areas, empirical evidence is needed to evaluate enthusiastic speculation. Evolutionary theory will continue to be extended, but there is no sign that it requires emendation.

Why we don't want another "Synthesis"

High-level debates in evolutionary biology often treat the Modern Synthesis as a framework of population genetics, or as an intellectual lineage with a changing distribution of beliefs. Unfortunately, these flexible notions, used to negotiate decades of innovations, are now thoroughly detached from their historical roots in the original Modern Synthesis (OMS), a falsifiable scientific theory. The OMS held that evolution can be adequately understood as a process of smooth adaptive change by shifting the frequencies of small-effect alleles at many loci simultaneously, without the direct involvement of new mutations. This shifting gene frequencies theory was designed to support a Darwinian view in which the course of evolution is governed by selection, and to exclude a mutation-driven view in which the timing and character of evolutionary change may reflect the timing and character of events of mutation. The OMS is not the foundation of current thinking, but a special case of a broader conception that includes (among other things) a mutation-driven view introduced by biochemists in the 1960s, and now widely invoked. This innovation is evident in mathematical models relating the rate of evolution directly to the rate of mutation, which emerged in 1969, and now represent a major branch of theory with many applications. In evo-devo, mutationist thinking is reflected by a concern for the "arrival of the fittest". Though evolutionary biology is not governed by any master theory, and incorporates views excluded from the OMS, the recognition of these changes has been hindered by woolly conceptions of theories, and by historical accounts, common in the evolutionary literature, that misrepresent the disputes that defined the OMS.

REVIEWERS

This article was reviewed by W. Ford Doolittle, Eugene Koonin and J. Peter Gogarten.