Joshua B. Plotkin

Coding-sequence determinants of gene expression in Escherichia coli.

Synonymous mutations do not alter the encoded protein, but they can influence gene expression. To investigate how, we engineered a synthetic library of 154 genes that varied randomly at synonymous sites, but all encoded the same green fluorescent protein (GFP). When expressed in Escherichia coli, GFP protein levels varied 250-fold across the library. GFP messenger RNA (mRNA) levels, mRNA degradation patterns, and bacterial growth rates also varied, but codon bias did not correlate with gene expression. Rather, the stability of mRNA folding near the ribosomal binding site explained more than half the variation in protein levels. In our analysis, mRNA folding and associated rates of translation initiation play a predominant role in shaping expression levels of individual genes, whereas codon bias influences global translation efficiency and cellular fitness.

Synonymous but not the same: the causes and consequences of codon bias

Despite their name, synonymous mutations have significant consequences for cellular processes in all taxa. As a result, an understanding of codon bias is central to fields as diverse as molecular evolution and biotechnology. Although recent advances in sequencing and synthetic biology have helped to resolve longstanding questions about codon bias, they have also uncovered striking patterns that suggest new hypotheses about protein synthesis. Ongoing work to quantify the dynamics of initiation and elongation is as important for understanding natural synonymous variation as it is for designing transgenes in applied contexts.

The population genetics of dN/dS.

Evolutionary pressures on proteins are often quantified by the ratio of substitution rates at non-synonymous and synonymous sites. The dN/dS ratio was originally developed for application to distantly diverged sequences, the differences among which represent substitutions that have fixed along independent lineages. Nevertheless, the dN/dS measure is often applied to sequences sampled from a single population, the differences among which represent segregating polymorphisms. Here, we study the expected dN/dS ratio for samples drawn from a single population under selection, and we find that in this context, dN/dS is relatively insensitive to the selection coefficient. Moreover, the hallmark signature of positive selection over divergent lineages, dN/dS>1, is violated within a population. For population samples, the relationship between selection and dN/dS does not follow a monotonic function, and so it may be impossible to infer selection pressures from dN/dS. These results have significant implications for the interpretation of dN/dS measurements among population-genetic samples.

Seed Dispersal and Spatial Pattern in Tropical Trees

Theories of tropical tree diversity emphasize dispersal limitation as a potential mechanism for separating species in space and reducing competitive exclusion. We compared the dispersal morphologies, fruit sizes, and spatial distributions of 561 tree species within a fully mapped, 50-hectare plot of primary tropical forest in peninsular Malaysia. We demonstrate here that the extent and scale of conspecific spatial aggregation is correlated with the mode of seed dispersal. This relationship holds for saplings as well as for mature trees. Phylogenetically independent contrasts confirm that the relationship between dispersal and spatial pattern is significant even after controlling for common ancestry among species. We found the same qualitative results for a 50-hectare tropical forest plot in Panama. Our results provide broad empirical evidence for the importance of dispersal mode in establishing the long-term community structure of tropical forests.

Mutational robustness can facilitate adaptation

Robustness seems to be the opposite of evolvability. If phenotypes are robust against mutation, we might expect that a population will have difficulty adapting to an environmental change, as several studies have suggested. However, other studies contend that robust organisms are more adaptable. A quantitative understanding of the relationship between robustness and evolvability will help resolve these conflicting reports and will clarify outstanding problems in molecular and experimental evolution, evolutionary developmental biology and protein engineering. Here we demonstrate, using a general population genetics model, that mutational robustness can either impede or facilitate adaptation, depending on the population size, the mutation rate and the structure of the fitness landscape. In particular, neutral diversity in a robust population can accelerate adaptation as long as the number of phenotypes accessible to an individual by mutation is smaller than the total number of phenotypes in the fitness landscape. These results provide a quantitative resolution to a significant ambiguity in evolutionary theory.

Rate-Limiting Steps in Yeast Protein Translation

Summary Deep sequencing now provides detailed snapshots of ribosome occupancy on mRNAs. We leverage these data to parameterize a computational model of translation, keeping track of every ribosome, tRNA, and mRNA molecule in a yeast cell. We determine the parameter regimes in which fast initiation or high codon bias in a transgene increases protein yield and infer the initiation rates of endogenous Saccharomyces cerevisiae genes, which vary by several orders of magnitude and correlate with 5′ mRNA folding energies. Our model recapitulates the previously reported 5′-to-3′ ramp of decreasing ribosome densities, although our analysis shows that this ramp is caused by rapid initiation of short genes rather than slow codons at the start of transcripts. We conclude that protein production in healthy yeast cells is typically limited by the availability of free ribosomes, whereas protein production under periods of stress can sometimes be rescued by reducing initiation or elongation rates.

Hemagglutinin sequence clusters and the antigenic evolution of influenza A virus

Continual mutations to the hemagglutinin (HA) gene of influenza A virus generate novel antigenic strains that cause annual epidemics. Using a database of 560 viral RNA sequences, we study the structure and tempo of HA evolution over the past two decades. We detect a critical length scale, in amino acid space, at which HA sequences aggregate into clusters, or swarms. We investigate the spatio-temporal distribution of viral swarms and compare it to the time series of the influenza vaccines recommended by the World Health Organization. We introduce a method for predicting future dominant HA amino acid sequences and discuss its potential relevance to vaccine choice. We also investigate the relationship between cluster structure and the primary antibody-combining regions of the HA protein.

Reconciling molecular phylogenies with the fossil record

Historical patterns of species diversity inferred from phylogenies typically contradict the direct evidence found in the fossil record. According to the fossil record, species frequently go extinct, and many clades experience periods of dramatic diversity loss. However, most analyses of molecular phylogenies fail to identify any periods of declining diversity, and they typically infer low levels of extinction. This striking inconsistency between phylogenies and fossils limits our understanding of macroevolution, and it undermines our confidence in phylogenetic inference. Here, we show that realistic extinction rates and diversity trajectories can be inferred from molecular phylogenies. To make this inference, we derive an analytic expression for the likelihood of a phylogeny that accommodates scenarios of declining diversity, time-variable rates, and incomplete sampling; we show that this likelihood expression reliably detects periods of diversity loss using simulation. We then study the cetaceans (whales, dolphins, and porpoises), a group for which standard phylogenetic inferences are strikingly inconsistent with fossil data. When the cetacean phylogeny is considered as a whole, recently radiating clades, such as the Balaneopteridae, Delphinidae, Phocoenidae, and Ziphiidae, mask the signal of extinctions. However, when isolating these groups, we infer diversity dynamics that are consistent with the fossil record. These results reconcile molecular phylogenies with fossil data, and they suggest that most extant cetaceans arose from four recent radiations, with a few additional species arising from clades that have been in decline over the last ∼10 Myr.

Codon bias and frequency-dependent selection on the hemagglutinin epitopes of influenza A virus

Although the surface proteins of human influenza A virus evolve rapidly and continually produce antigenic variants, the internal viral genes acquire mutations very gradually. In this paper, we analyze the sequence evolution of three influenza A genes over the past two decades. We study codon usage as a discriminating signature of gene- and even residue-specific diversifying and purifying selection. Nonrandom codon choice can increase or decrease the effective local substitution rate. We demonstrate that the codons of hemagglutinin, particularly those in the antibody-combining regions, are significantly biased toward substitutional point mutations relative to the codons of other influenza virus genes. We discuss the evolutionary interpretation and implications of these biases for hemagglutinins antigenic evolution. We also introduce information-theoretic methods that use sequence data to detect regions of recent positive selection and potential protein conformational changes.

Reinforcement of pre-zygotic isolation and karyotype evolution in Agrodiaetus butterflies

The reinforcement model of evolution argues that natural selection enhances pre-zygotic isolation between divergent populations or species by selecting against unfit hybrids or costly interspecific matings. Reinforcement is distinguished from other models that consider the formation of reproductive isolation to be a by-product of divergent evolution. Although theory has shown that reinforcement is a possible mechanism that can lead to speciation, empirical evidence has been sufficiently scarce to raise doubts about the importance of reinforcement in nature. Agrodiaetus butterflies (Lepidoptera: Lycaenidae) exhibit unusual variability in chromosome number. Whereas their genitalia and other morphological characteristics are largely uniform, different species vary considerably in male wing colour, and provide a model system to study the role of reinforcement in speciation. Using comparative phylogenetic methods, we show that the sympatric distribution of 15 relatively young sister taxa of Agrodiaetus strongly correlates with differences in male wing colour, and that this pattern is most likely the result of reinforcement. We find little evidence supporting sympatric speciation: rather, in Agrodiaetus, karyotypic changes accumulate gradually in allopatry, prompting reinforcement when karyotypically divergent races come into contact.