Evolutionary honing in and mutational replacement: how long-term directed mutational responses to specific environmental pressures are possible.

Recent results have shown that the human malaria-resistant hemoglobin S mutation originates de novo more frequently in the gene and in the population where it is of adaptive significance, namely, in the hemoglobin subunit beta gene compared to the nonresistant but otherwise identical 20A[Formula: see text]T mutation in the hemoglobin subunit delta gene, and in sub-Saharan Africans, who have been subject to intense malarial pressure for many generations, compared to northern Europeans, who have not. This finding raises a fundamental challenge to the traditional notion of accidental mutation. Here, we address this finding with the replacement hypothesis, according to which preexisting genetic interactions can lead directly and mechanistically to mutations that simplify and replace them. Thus, an evolutionary process under selection can gradually hone in on interactions of importance for the currently evolving adaptations, from which large-effect mutations follow that are relevant to these adaptations. We exemplify this hypothesis using multiple types of mutation, including gene fusion mutations, gene duplication mutations, A[Formula: see text]G mutations in RNA-edited sites and transcription-associated mutations, and place it in the broader context of a system-level view of mutation origination called interaction-based evolution. Potential consequences include that similarity of mutation pressures may contribute to parallel evolution in genetically related species, that the evolution of genome organization may be driven by mutational mechanisms, that transposable element movements may also be explained by replacement, and that long-term directed mutational responses to specific environmental pressures are possible. Such mutational phenomena need to be further tested by future studies in natural and artificial settings.

Genes that are Used Together are More Likely to be Fused Together in Evolution by Mutational Mechanisms: A Bioinformatic Test of the Used-Fused Hypothesis.

Cases of parallel or recurrent gene fusions in evolution as well as in genetic disease and cancer are difficult to explain, because unlike point mutations, they can require the repetition of a similar configuration of multiple breakpoints rather than the repetition of a single point mutation. The used-together-fused-together hypothesis holds that genes that are used together repeatedly and persistently in a specific context are more likely to undergo fusion mutation in the course of evolution for mechanistic reasons. This hypothesis offers to explain gene fusion in both evolution and disease under one umbrella. Using bioinformatic data, we tested this hypothesis against alternatives, including that all gene pairs can fuse by random mutation, but among pairs thus fused, those that had interacted previously are more likely to be favored by selection. Results show that across multiple measures of gene interaction, human genes whose orthologs are fused in one or more species are more likely to interact with each other than random pairs of genes of the same genomic distance between pair members; that an overlap exists between genes that fused in the course of evolution in non-human species and genes that undergo fusion in human cancers; and that across six primate species studied, fusions predominate over fissions and exhibit substantial evolutionary parallelism. Together, these results support the used-together-fused-together hypothesis over its alternatives. Multiple implications are discussed, including the relevance of mutational mechanisms to the evolution of genome organization, to the distribution of fitness effects of mutation, to evolutionary parallelism and more. Supplementary Information: The online version contains supplementary material available at 10.1007/s11692-022-09579-9.

De novo mutation rates at the single-mutation resolution in a human HBB gene region associated with adaptation and genetic disease.

Although it is known that the mutation rate varies across the genome, previous estimates were based on averaging across various numbers of positions. Here, we describe a method to measure the origination rates of target mutations at target base positions and apply it to a 6-bp region in the human hemoglobin subunit beta (HBB) gene and to the identical, paralogous hemoglobin subunit delta (HBD) region in sperm cells from both African and European donors. The HBB region of interest (ROI) includes the site of the hemoglobin S (HbS) mutation, which protects against malaria, is common in Africa, and has served as a classic example of adaptation by random mutation and natural selection. We found a significant correspondence between de novo mutation rates and past observations of alleles in carriers, showing that mutation rates vary substantially in a mutation-specific manner that contributes to the site frequency spectrum. We also found that the overall point mutation rate is significantly higher in Africans than in Europeans in the HBB region studied. Finally, the rate of the 20A→T mutation, called the "HbS mutation" when it appears in HBB, is significantly higher than expected from the genome-wide average for this mutation type. Nine instances were observed in the African HBB ROI, where it is of adaptive significance, representing at least three independent originations; no instances were observed elsewhere. Further studies will be needed to examine mutation rates at the single-mutation resolution across these and other loci and organisms and to uncover the molecular mechanisms responsible.

Expanding the binding specificity for RNA recognition by a PUF domain.

The ability to design a protein to bind specifically to a target RNA enables numerous applications, with the modular architecture of the PUF domain lending itself to new RNA-binding specificities. For each repeat of the Pumilio-1 PUF domain, we generate a library that contains the 8,000 possible combinations of amino acid substitutions at residues critical for RNA contact. We carry out yeast three-hybrid selections with each library against the RNA recognition sequence for Pumilio-1, with any possible base present at the position recognized by the randomized repeat. We use sequencing to score the binding of each variant, identifying many variants with highly repeat-specific interactions. From these data, we generate an RNA binding code specific to each repeat and base. We use this code to design PUF domains against 16 RNAs, and find that some of these domains recognize RNAs with two, three or four changes from the wild type sequence.

Combining natural sequence variation with high throughput mutational data to reveal protein interaction sites.

Many protein interactions are conserved among organisms despite changes in the amino acid sequences that comprise their contact sites, a property that has been used to infer the location of these sites from protein homology. In an inter-species complementation experiment, a sequence present in a homologue is substituted into a protein and tested for its ability to support function. Therefore, substitutions that inhibit function can identify interaction sites that changed over evolution. However, most of the sequence differences within a protein family remain unexplored because of the small-scale nature of these complementation approaches. Here we use existing high throughput mutational data on the in vivo function of the RRM2 domain of the Saccharomyces cerevisiae poly(A)-binding protein, Pab1, to analyze its sites of interaction. Of 197 single amino acid differences in 52 Pab1 homologues, 17 reduce the function of Pab1 when substituted into the yeast protein. The majority of these deleterious mutations interfere with the binding of the RRM2 domain to eIF4G1 and eIF4G2, isoforms of a translation initiation factor. A large-scale mutational analysis of the RRM2 domain in a two-hybrid assay for eIF4G1 binding supports these findings and identifies peripheral residues that make a smaller contribution to eIF4G1 binding. Three single amino acid substitutions in yeast Pab1 corresponding to residues from the human orthologue are deleterious and eliminate binding to the yeast eIF4G isoforms. We create a triple mutant that carries these substitutions and other humanizing substitutions that collectively support a switch in binding specificity of RRM2 from the yeast eIF4G1 to its human orthologue. Finally, we map other deleterious substitutions in Pab1 to inter-domain (RRM2-RRM1) or protein-RNA (RRM2-poly(A)) interaction sites. Thus, the combined approach of large-scale mutational data and evolutionary conservation can be used to characterize interaction sites at single amino acid resolution.

Overexpression of eukaryotic initiation factor 5 rescues the translational defect of tpk1w in a manner that necessitates a novel phosphorylation site.

Cells respond to changes in their environment through mechanisms that often necessitate reprogramming of the translation machinery. The fastest and strongest of all tested responses is the translation inhibition observed following abrupt depletion of glucose from the media of yeast cells. The speed of the response suggests a post-translational modification of a key component of the translation machinery. This translation factor is as yet unknown. A cAMP-dependent protein kinase mutant yeast strain (tpk1(w)) that does not respond properly to glucose depletion and maintains translation was described previously. We hypothesized that the inability of tpk1(w) to arrest translation results from abnormal expression of key translation mediators. Genome-wide analysis of steady-state mRNA levels in tpk1(w) revealed underexpression of several candidates. Elevating the cellular levels of eukaryotic initiation factor (eIF) 5 by overexpression rescued the translational defect of tpk1(w). Restoring ribosomal dissociation by eIF5 necessitated an active GAP domain and multiple regions throughout this protein. Phosphoproteomics analysis of wild-type cells overexpressing eIF5 revealed increased phosphorylation in a novel site (Thr191) upon glucose depletion. Mutating this residue and introducing it into tpk1(w) abolished the ability of eIF5 to rescue the translational defect. Intriguingly, introducing this mutation into the wild-type strain did not hamper its translational response. We further show that Thr191 is phosphorylated in vitro by Casein Kinase II (CKII), and yeast cells with a mutated CKII have a reduced response to glucose depletion. These results implicate phosphorylation of eIF5 at Thr191 by CKII as one of the pathways for regulating translation upon glucose depletion.

Deep mutational scanning of an RRM domain of the Saccharomyces cerevisiae poly(A)-binding protein.

The RNA recognition motif (RRM) is the most common RNA-binding domain in eukaryotes. Differences in RRM sequences dictate, in part, both RNA and protein-binding specificities and affinities. We used a deep mutational scanning approach to study the sequence-function relationship of the RRM2 domain of the Saccharomyces cerevisiae poly(A)-binding protein (Pab1). By scoring the activity of more than 100,000 unique Pab1 variants, including 1246 with single amino acid substitutions, we delineated the mutational constraints on each residue. Clustering of residues with similar mutational patterns reveals three major classes, composed principally of RNA-binding residues, of hydrophobic core residues, and of the remaining residues. The first class also includes a highly conserved residue not involved in RNA binding, G150, which can be mutated to destabilize Pab1. A comparison of the mutational sensitivity of yeast Pab1 residues to their evolutionary conservation reveals that most residues tolerate more substitutions than are present in the natural sequences, although other residues that tolerate fewer substitutions may point to specialized functions in yeast. An analysis of ∼40,000 double mutants indicates a preference for a short distance between two mutations that display an epistatic interaction. As examples of interactions, the mutations N139T, N139S, and I157L suppress other mutations that interfere with RNA binding and protein stability. Overall, this study demonstrates that living cells can be subjected to a single assay to analyze hundreds of thousands of protein variants in parallel.

Genome-wide analysis of RNA extracted from isolated mitochondria.

Isolating mitochondria by subcellular fractionation is a well-established method for retrieving intact and functional mitochondria. This procedure has been used to identify proteins of the mitochondria and to explore import mechanisms. Using the same method, it was shown that mitochondria can be purified along with cytoplasmic ribosomes and nuclear-encoded mRNAs attached to the outer membrane. Combining this procedure with DNA microarray analysis allows for global identification of the mRNAs associated with mitochondria, and hence a better understanding of the underlying molecular mechanisms. In this chapter, we will describe a procedure for the isolation of mitochondria from yeast and RNA purification. We will then describe the process of labeling and hybridization to DNA microarrays, and comment on a few aspects of the data analysis.

Asc1 supports cell-wall integrity near bud sites by a Pkc1 independent mechanism.

BACKGROUND: The yeast ribosomal protein Asc1 is a WD-protein family member. Its mammalian ortholog, RACK1 was initially discovered as a receptor for activated protein C kinase (PKC) that functions to maintain the active conformation of PKC and to support its movement to target sites. In the budding yeast though, a connection between Asc1p and the PKC signaling pathway has never been reported. METHODOLOGY/PRINCIPAL FINDINGS: In the present study we found that asc1-deletion mutant (asc1Delta) presents some of the hallmarks of PKC signaling mutants. These include an increased sensitivity to staurosporine, a specific Pkc1p inhibitor, and susceptibility to cell-wall perturbing treatments such as hypotonic- and heat shock conditions and zymolase treatment. Microscopic analysis of asc1Delta cells revealed cell-wall invaginations near bud sites after exposure to hypotonic conditions, and the dynamic of cells' survival after this stress further supports the involvement of Asc1p in maintaining the cell-wall integrity during the mid-to late stages of bud formation. Genetic interactions between asc1 and pkc1 reveal synergistic sensitivities of a double-knock out mutant (asc1Delta/pkc1Delta) to cell-wall stress conditions, and high basal level of PKC signaling in asc1Delta. Furthermore, Asc1p has no effect on the cellular distribution or redistribution of Pkc1p at optimal or at cell-wall stress conditions. CONCLUSIONS/SIGNIFICANCE: Taken together, our data support the idea that unlike its mammalian orthologs, Asc1p acts remotely from Pkc1p, to regulate the integrity of the cell-wall. We speculate that its role is exerted through translation regulation of bud-site related mRNAs during cells' growth.

Tom20 mediates localization of mRNAs to mitochondria in a translation-dependent manner.

mRNAs encoding mitochondrial proteins are enriched in the vicinity of mitochondria, presumably to facilitate protein transport. A possible mechanism for enrichment may involve interaction of the translocase of the mitochondrial outer membrane (TOM) complex with the precursor protein while it is translated, thereby leading to association of polysomal mRNAs with mitochondria. To test this hypothesis, we isolated mitochondrial fractions from yeast cells lacking the major import receptor, Tom20, and compared their mRNA repertoire to that of wild-type cells by DNA microarrays. Most mRNAs encoding mitochondrial proteins were less associated with mitochondria, yet the extent of decrease varied among genes. Analysis of several mRNAs revealed that optimal association of Tom20 target mRNAs requires both translating ribosomes and features within the encoded mitochondrial targeting signal. Recently, Puf3p was implicated in the association of mRNAs with mitochondria through interaction with untranslated regions. We therefore constructed a tom20 Delta puf3 Delta double-knockout strain, which demonstrated growth defects under conditions where fully functional mitochondria are required. Mislocalization effects for few tested mRNAs appeared stronger in the double knockout than in the tom20 Delta strain. Taken together, our data reveal a large-scale mRNA association mode that involves interaction of Tom20p with the translated mitochondrial targeting sequence and may be assisted by Puf3p.

Exploring translation regulation by global analysis of ribosomal association.

Translation efficiency of an mRNA is related in most cases to its ribosomal association. This association can be readily measured through the separation of cellular complexes on sucrose gradients by velocity sedimentation, and identification of the sedimentation position of the mRNA in the gradient. Since ribosomes are the main driving force for mRNA sedimentation, sedimentation position is highly correlated with ribosomal association and thus translation efficiency. The advent of DNA microarrays allowed the determination of ribosomal association for many mRNAs in parallel through the combination of fractionation in a sucrose gradient followed by microarray analysis. This provided an enormous amount of novel information regarding translation control and regulation. Herein we provide a detailed protocol for performing such an analysis, indicating important points for consideration and discussing some of the advantages and limitations of this powerful approach.

Yeast translational response to high salinity: global analysis reveals regulation at multiple levels.

Genome-wide studies of steady-state mRNA levels revealed common principles underlying transcriptional changes in response to external stimuli. To uncover principles that govern other stages of the gene-expression response, we analyzed the translational response and its coordination with transcriptome changes following exposure to severe stress. Yeast cells were grown for 1 h in medium containing 1 M NaCl, which elicits a maximal but transient translation inhibition, and nonpolysomal or polysomal mRNA pools were subjected to DNA-microarray analyses. We observed a strong repression in polysomal association for most mRNAs, with no simple correlation with the changes in transcript levels. This led to an apparent accumulation of many mRNAs as a nontranslating pool, presumably waiting for recovery from the stress. However, some mRNAs demonstrated a correlated change in their polysomal association and their transcript levels (i.e., potentiation). This group was enriched with targets of the transcription factors Msn2/Msn4, and the translational induction of several tested mRNAs was diminished in an Msn2/Msn4 deletion strain. Genome-wide analysis of a strain lacking the high salinity response kinase Hog1p revealed that the group of translationally affected genes is significantly enriched with motifs that were shown to be associated with the ARE-binding protein Pub1. Since a relatively small number of genes was affected by Hog1p deletion, additional signaling pathways are likely to be involved in coordinating the translational response to severe salinity stress.

Genome-wide analysis of mRNA polysomal profiles with spotted DNA microarrays.

The sedimentation of an mRNA in sucrose gradients is highly affected by its ribosomal association. Sedimentation analysis has therefore become routine for studying changes in ribosomal association of mRNAs of interest. DNA microarray technology has been combined with sedimentation analysis to characterize changes in ribosomal association for thousands of mRNAs in parallel. Such analyses revealed mRNAs that are translationally regulated and have provided new insights into the translation process. In this chapter, we describe possible experimental designs for analyzing genome-wide changes in ribosomal association, and discuss some of their advantages and disadvantages. We then provide a detailed protocol for analysis of polysomal fractions using spotted DNA microarrays.

The RNA polymerase II subunit Rpb4p mediates decay of a specific class of mRNAs.

It is commonly appreciated that the mRNA level is determined by the balance between its synthetic and decay kinetics. Yet, little is known about coordination between these distinct processes. A major pathway of the eukaryotic mRNA decay initiates with shortening of the mRNA poly(A) tail (deadenylation), followed by removal of the mRNA 5' cap structure and its subsequent exonucleolytic degradation. Here we report that a subunit of RNA polymerase II, Rpb4p, is required for the decay of a class of mRNAs whose products are involved in protein synthesis. Cells lacking RPB4 are defective in the deadenylation and post-deadenylation steps of representatives of this class of mRNAs. Moreover, Rpb4p interacts with both the mRNP and with subunits of the mRNA decay complex Pat1/Lsm1-7 that enhances decapping. Consistently, a portion of Rpb4p is localized in P bodies, where mRNA decapping and degradation is executed, and mutations in RPB4 increase the number of P bodies per cell. We propose that Rpb4p has a dual function in mRNA decay. It promotes or enhances the deadenylation process of specific mRNAs and recruits Pat1/Lsm1-7 to these mRNAs, thus stimulating their decapping and further decay. In this way, Rpb4p might link the activity of the basal transcription apparatus with that of the mRNA decay machinery.

Single point mutations in the zinc finger motifs of the human immunodeficiency virus type 1 nucleocapsid alter RNA binding specificities of the gag protein and enhance packaging and infectivity.

A specific interaction between the nucleocapsid (NC) domain of the Gag polyprotein and the RNA encapsidation signal (Psi) is required for preferential incorporation of the retroviral genomic RNA into the assembled virion. Using the yeast three-hybrid system, we developed a genetic screen to detect human immunodeficiency virus type 1 (HIV-1) Gag mutants with altered RNA binding specificities. Specifically, we randomly mutated full-length HIV-1 Gag or its NC portion and screened the mutants for an increase in affinity for the Harvey murine sarcoma virus encapsidation signal. These screens identified several NC zinc finger mutants with altered RNA binding specificities. Furthermore, additional zinc finger mutants that also demonstrated this phenotype were made by site-directed mutagenesis. The majority of these mutants were able to produce normal virion-like particles; however, when tested in a single-cycle infection assay, some of the mutants demonstrated higher transduction efficiencies than that of wild-type Gag. In particular, the N17K mutant showed a seven- to ninefold increase in transduction, which correlated with enhanced vector RNA packaging. This mutant also packaged larger amounts of foreign RNA. Our results emphasize the importance of the NC zinc fingers, and not other Gag sequences, in achieving specificity in the genome encapsidation process. In addition, the described mutations may contribute to our understanding of HIV diversity resulting from recombination events between copackaged viral genomes and foreign RNA.

The conserved carboxy terminus of the capsid domain of human immunodeficiency virus type 1 gag protein is important for virion assembly and release.

The retroviral Gag precursor plays an important role in the assembly of virion particles. The capsid (CA) protein of the Gag molecule makes a major contribution to this process. In the crystal structure of the free CA protein of the human immunodeficiency virus type 1 (HIV-1), 11 residues of the C terminus were found to be unstructured, and to date no information exists on the structure of these residues in the context of the Gag precursor molecule. We performed phylogenetic analysis and demonstrated a high degree of conservation of these 11 amino acids. Deletion of this cluster or introduction of various point mutations into these residues resulted in significant impairment of particle infectivity. In this cluster, two putative structural regions were identified, residues that form a hinge region (353-VGGP-356) and those that contribute to an alpha-helix (357-GHKARVL-363). Overall, mutations in these regions resulted in inhibition of virion production, but mutations in the hinge region demonstrated the most significant reduction. Although all the Gag mutants appeared to have normal Gag-Gag and Gag-RNA interactions, the hinge mutants were characterized by abnormal formation of cytoplasmic Gag complexes. Gag proteins with mutations in the hinge region demonstrated normal membrane association but aberrant rod-like membrane structures. More detailed analysis of these structures in one of the mutants demonstrated abnormal trapped Gag assemblies. These data suggest that the conserved CA C terminus is important for HIV-1 virion assembly and release and define a putative target for drug design geared to inhibit the HIV-1 assembly process.