Mutation hotspots in the p53 gene in tumors of different origin: correlation with evolutionary conservation and signs of positive selection

We present a classification analysis of the mutation spectra of the p53 gene and construct maps of hotspots for the germline (Li-Fraumein syndrome), different types of tumors and their derived cell lines. While spectra from solid tumors share common hotspots with the germline spectrum, they also contain unique sets of somatic hotspots that are not observed in the germline. All these hotspots correspond to amino acid replacements in the DNA-binding interface of p53. The mutation spectra of lymphomas and cell lines derived from lymphomas and lung cancers contained few hotspots compared to solid tumors. Thus, the distribution of hotspots in the p53 gene appears to depend on the tumor type and cell growth conditions; this specificity is missed by the bulk hotspot analysis. A negative correlation was detected between the amino acid replacement propensity in tumors and evolutionary variability: the hotspots are located in the positions that are highly conserved in p53 and its paralogs, p63 and p73. In all the mutation spectra, substitutions leading to amino acid replacements strongly dominate over silent substitutions, indicating that functional sites evolving under strong purifying selection are subject to intensive positive selection in p53-dependent tumors. These results are compatible with the gain-of-function concept of the role of p53 in tumorigenesis.     

Detecting compensatory covariation signals in protein evolution using reconstructed ancestral sequences

When protein sequences divergently evolve under functional constraints, some individual amino acid replacements that reverse the charge (e.g. Lys to Asp) may be compensated by a replacement at a second position that reverses the charge in the opposite direction (e.g. Glu to Arg). When these side-chains are near in space (proximal), such double replacements might be driven by natural selection, if either is selectively disadvantageous, but both together restore fully the ability of the protein to contribute to fitness (are together "neutral"). Accordingly, many have sought to identify pairs of positions in a protein sequence that suffer compensatory replacements, often as a way to identify positions near in space in the folded structure. A "charge compensatory signal" might manifest itself in two ways. First, proximal charge compensatory replacements may occur more frequently than predicted from the product of the probabilities of individual positions suffering charge reversing replacements independently. Conversely, charge compensatory pairs of changes may be observed to occur more frequently in proximal pairs of sites than in the average pair. Normally, charge compensatory covariation is detected by comparing the sequences of extant proteins at the "leaves" of phylogenetic trees. We show here that the charge compensatory signal is more evident when it is sought by examining individual branches in the tree between reconstructed ancestral sequences at nodes in the tree. Here, we find that the signal is especially strong when the positions pairs are in a single secondary structural unit (e.g. alpha helix or beta strand) that brings the side-chains suffering charge compensatory covariation near in space, and may be useful in secondary structure prediction. Also, "node-node" and "node-leaf" compensatory covariation may be useful to identify the better of two equally parsimonious trees, in a way that is independent of the mathematical formalism used to construct the tree itself. Further, compensatory covariation may provide a signal that indicates whether an episode of sequence evolution contains more or less divergence in functional behavior. Compensatory covariation analysis on reconstructed evolutionary trees may become a valuable tool to analyze genome sequences, and use these analyses to extract biomedically useful information from proteome databases.     

Primary structure of pronghorn pancreatic ribonuclease: Close relationship between giraffe and pronghorn

Summary Pancreatic ribonuclease from pronghorn (Antilocapra americana) was isolated and its amino acid sequence was determined from a tryptic digest of the performic acid-oxidized protein. Peptides were positioned by homology with other ribonucleases. Only peptides that differed in amino acid composition from the corresponding peptides of ox or goat ribonucleases were sequenced.In a most parsimonious tree of pancreatic ribonucleases, pronghorn and giraffe were placed together and these two were placed with the bovids, leaving the deer as a taxon separate from the other ruminants. The amino acid replacements that determine this tree topology are three rarely occurring replacements shared by pronghorn and giraffe. Notwithstanding their close phylogenetic relationship, both ribonucleases differ strongly in extent of glycosidation, net charge and antigenic properties.     

The molecular basis of adaptive evolution of squirrelfish rhodopsins

The wavelengths of maximal absorption (lambdamax) of the rhodopsins of nine squirrelfishes (N. sammara, N. argenteus, S. punctatissimum, S. microstoma, S. diadema, S. xantherythrum, S. spiniferum, N. aurolineatus, and S. tiere) and two soldierfishes (M. violacea and M. berndti) vary between 481 and 502 nm. Phylogenetic and mutagenesis analyses suggest that the common ancestor of these pigments had a lambdamax value of approximately 493 nm, and the contemporary lambdamax values were generated mostly by amino acid replacements E122M, F261Y, and A292S. The probability of observing all these amino acid replacements at specific branches of the phylogenetic tree is only 2.5 x 10(-9); it is highly unlikely that these changes have occurred by neutral evolution. Because of a close association between the lambdamax values of these pigments and the wavelengths of light available to the corresponding species, the excess number of amino acid changes at specific branches in the phylogenetic tree strongly suggests that the rhodopsins have undergone adaptive changes at various stages of the holocentrid evolution.     

Copper/Zinc superoxide dismutase: How likely is gene transfer from ponyfish toPhotobacterium leiognathi?

Summary The proposed transfer of the gene for Cu/Zn superoxide dismutase from the ponyfish to its symbiotic bacteriumPhotobacterium leiognathi has been evaluated by an extensive analysis of all available Cu/Zn superoxide dismutase sequences. By the use of four different computer programs, phylogenetic trees were constructed from the sequences of the superoxide dismutases of human, ox, pig, horse, swordfish, fruit fly, yeast, andNeurospora crassa to find out whether superoxide dismutase sequences can reliably be used for the reconstruction of genealogical relationships. All programs arrived at the same most parsimonious tree (one requiring 232 amino acid replacements), the topology of which conformed to established opinions about the phylogenetic relations among these eukaryotes, except that it placed humans closer to the artiodactyls ox and pig than it placed horses. This could be corrected at the cost of two amino acid replacements. The sequence ofP. leiognathi superoxide dismutase was then connected at all possible positions to the corrected eukaryotic tree. It was slighly more parsimonious to link the bacterial sequence to the root of the tree than to the fish branch: The former required 316 (or 317) amino acid replacements, versus 319 for the latter. This relative lack of discrimination between such distinct alternative topologies may be a general complication in the comparison of prokaryotic and eukaryotic proteins: Bacterial cytochrome c sequences also were found to be connected as parsimoniously to the root of the eukaryotic tree as to any terminal or ancestral branch. It was calculated that the rate of evolution of the bacterial superoxide dismutase gene, if transfer occurred 30 million years (Myr) ago, must have amounted to 487 amino acid replacements per 100 residues per 100 Myr. This is more than 5 times the highest rate observed in any protein (that found for fibrinopeptides), and even much higher than the maximum rate of protein evolution that can be deduced from the neutral mutation rate of unconstrained DNA. Also, no significant evidence that shared derived amino acid replacements are present in swordfish andP. leiognathi superoxide dismutase, as might be expected had gene transfer occurred, was found. On the basis of the available data it seems more reasonable to ascribe the isolated occurrence of Cu/Zn superoxide dismutase inP. leiognathi (as well as inCaulobacter crescentus) to irregular patterns of gene expression and inactivation in the course of divergent evolution than to undocumented processes of gene transfer from eukaryotes to prokaryotes.     

Genetic studies of the lac repressor. XII. Amino acid replacements in the DNA binding domain of the Escherichia coli lac repressor

Out of the first 62 residues of the lac repressor, 38 positions have been substituted by at least one amino acid exchange. The total number of replacements in this region is 131. Data from several studies are considered.     

Mutations and selection in the generation of class II histocompatibility antigen polymorphism.

A comparison of seven human DR and DC class II histocompatibility antigen beta-chain amino acid sequences indicates that the allelic variation is of comparable magnitude within the DR and DC beta-chain genes. Silent and replacement nucleotide substitutions in six DR and DC beta-chain sequences, as well as in seven murine class II sequences (three I-A beta and four I-A alpha alleles) were analyzed. The results suggest that the mutation rates are of a comparable magnitude in the nucleotide sequences encoding the first and second external domains of the class II molecules. Nevertheless, the allelic amino acid replacements are predominantly located in the first domains. We conclude that a conservative selective pressure acts on the second domains, whereas in many positions in the first domains replacement substitutions are selectively neutral or maybe even favoured. Thus, the difference between the first and second domains as regards the number of amino acid replacements is mainly due to selection.     

On the PAM matrix model of protein evolution

The internal consistency of the PAM matrix model of protein evolution is here investigated. The 1 PAM matrix has been constructed from amino acid replacements observed in closely related sequences. Such replacements are of two types, those that do not require an intermediate amino acid replacement and those that do. The second type of replacement must generally be produced by a repetition of the first. This allows data on the first type to be used in predicting data on the second type so that some elements of the 1 PAM matrix may be used to predict others. A discrepancy of more than two orders of magnitude is found between the predictions and the data when this is carried out. This is partly accounted for by an error in constructing the matrix. However, it also seems necessary that the basic model be modified. Several possibilities are considered. One of these is to incorporate a site-dependent spectrum of mutabilities associated with each amino acid.      

Patterns of nucleotide substitutions inferred from the phylogenies of the class I major histocompatibility complex genes

Summary Patterns of nucleotide substitutions in human major histocompatibility complex (MHC) class I genes were estimated by using phylogenetic trees of DNA sequences. The pattern is defined as a set of 12 parameters, each of which represents the relative frequency of substitutions from a particular nucleotide to another. The pattern at the antigen recognition sites (ARS) in functional MHC genes was remarkably different from that at the remaining coding region (non-ARS). In particular, the proportion of transitions among all the nucleotide substitutions (P _s) was extremely low at the third codon positions of ARS. In the HLA-A genes, P _s at the third codon positions was only 6% in ARS, whereas it was 69% in non-ARS. In HLA-B, the corresponding values were 30% in ARS and 80% in non-ARS, respectively. On the other hand, P _s in a class I pseudogene (HLA-H) was 57%, which was in good agreement with P _s in other pseudogenes. Because pseudogenes are selectively neutral, the pattern in pseudogenes is regarded as the pattern of spontaneous substitution mutations. In general, the pattern in functional genes that are subject to selective forces deviates from the pattern in pseudogenes. At the third codon positions in coding regions, transitions scarcely cause amino acid replacements, whereas about half of transversions do cause replacements. Accordingly, P _s at the third codon positions decreases if amino acid replacements are accelerated by natural selection but increases if amino acids are conserved by functional constraint. Our observations imply that the ARS region is subject to natural selection favoring amino acid replacements, whereas the non-ARS region is subject to functional constraint. Offprint requests to: T. Gojobori     

Determination of a specific region of the purine–cytosine permease involved in the recognition of its substrates

Three u.v.-induced mutants of the purine-cytosine permease gene of Saccharomyces cerevisiae, with altered apparent Michaelis constant of transport (Kmapp), were cloned and sequenced. One of the mutants had extensive nucleotide replacement, whereas the other two had a single mutation. To evaluate the contribution of the different amino acid replacements to the phenotype of the complex mutant, simpler mutants were created by site-directed mutagenesis. All the amino acid replacements found in the segment from amino acids 371 to 377 inclusive, contribute to the determination of the phenotype. According to the model postulated this segment lies on the cell surface. In particular, amino acids at position 374 and 377 modulate the affinity of the permease towards its substrates. In the wild-type, when asparagine is present at both of these positions, the lowest Kmapp values are found.     

A guest-host approach to oligodepsipeptide structure

In this paper we investigate the effect of main chain isosteric replacement of specific amino acid residues by α-hydroxy acids. As part of a long term program specifically protected heptaglutamates were prepared and their circular dichroism and nuclear magnetic resonance spectra in various solvents were examined. From these experiments conformational preferences were deduced. We have also prepared oligo-(γ-methyl-glutamates) replacing the amino acids at specific positions along the chain with S-lactic acid and have elucidated the effect of these main chain isosteric replacements on oligopeptide structure. Analogues of collagen also have been prepared with glycolic acid replacing specific glycine residues. We synthesized the model hexamers Ac-Ala-Gly-Pro-Ala-Gly-Pro-NHMe, Ac-Ala-Glc-Pro-Ala-Gly-Pro-NHMe, and Ac-Ala-Gly-Pro-Ala-Glc-Pro-NHMe in order to study their structural characteristics under various conditions. Preliminary nuclear magnetic resonance and circular dichroism results are presented.     

Protein sequence evidence for monophyly of the carnivore families Procyonidae and Mustelidae

The amino acid sequence of the eye lens protein alpha-crystallin A of the ring-tailed cat, Bassariscus astutus, has been determined. The sequence of the Bassariscus alpha A chain, which is 173 residues long, was compared with the previously determined set of 41 mammalian alpha A sequences. Among the investigated carnivores (dog, cat, sloth bear, American mink, gray seal, and California sea lion) the Bassariscus alpha A sequence exclusively shares two amino acid replacements with the alpha A chain of the mink, Mustela vison: 7 His----Gln and 61 Ile--- -Val. The Mustela and Bassariscus alpha A sequences differ at only three positions and have no replacements in common with any of the other investigated carnivore alpha A chains. Furthermore, the replacement 7 His----Gln has only been found in three-toed sloth, whereas 61 Ile----Val occurs scattered in three other taxa: pig, rhinoceros, and prosimians. It thus is most parsimonious to join Bassariscus and Mustela--and consequently their respective families, Procyonidae and Mustelidae--as sister groups in the phylogenetic tree of mammalian alpha A sequences.      

Bayesian analysis suggests that most amino acid replacements in Drosophila are driven by positive selection

One of the principal goals of population genetics is to understand the processes by which genetic variation within species (polymorphism) becomes converted into genetic differences between species (divergence). In this transformation, selective neutrality, near neutrality, and positive selection may each play a role, differing from one gene to the next. Synonymous nucleotide sites are often used as a uniform standard of comparison across genes on the grounds that synonymous sites are subject to relatively weak selective constraints and so may, to a first approximation, be regarded as neutral. Synonymous sites are also interdigitated with nonsynonymous sites and so are affected equally by genomic context and demographic factors. Hence a comparison of levels of polymorphism and divergence between synonymous sites and amino acid replacement sites in a gene is potentially informative about the magnitude of selective forces associated with amino acid replacements. We have analyzed 56 genes in which polymorphism data from D. simulans are compared with divergence from a reference strain of D. melanogaster. The framework of the analysis is Bayesian and assumes that the distribution of selective effects (Malthusian fitnesses) is Gaussian with a mean that differs for each gene. In such a model, the average scaled selection intensity (gamma = N(e)s) of amino acid replacements eligible to become polymorphic or fixed is -7.31, and the standard deviation of selective effects within each locus is 6.79 (assuming homoscedasticity across loci). For newly arising mutations of this type that occur in autosomal or X-linked genes, the average proportion of beneficial mutations is 19.7%. Among the amino acid polymorphisms in the sample, the expected average proportion of beneficial mutations is 47.7%, and among amino acid replacements that become fixed the average proportion of beneficial mutations is 94.3%. The average scaled selection intensity of fixed mutations is +5.1. The presence of positive selection is pervasive with the single exception of kl-5, a Y-linked fertility gene. We find no evidence that a significant fraction of fixed amino acid replacements is neutral or nearly neutral or that positive selection drives amino acid replacements at only a subset of the loci. These results are model dependent and we discuss possible modifications of the model that might allow more neutral and nearly neutral amino acid replacements to be fixed.     

Inferring the number of evolutionary events from DNA coding sequence differences

The estimation of the amount of evolutionary divergence that has taken place between two DNA coding sequences depends strongly on the degree of constraint on amino acid replacements. If amino acid replacements are relatively unconstrained, the individual nucleotide is the appropriate unit of analysis and the method of Tajima and Nei can be used. If amino acid replacements are constrained, however, this method is shown to be inapplicable. For sequences with strong amino acid constraints, a method is outlined analogous to the Tajima and Nei method using codons as the unit of analysis. Only synonymous substitutions are used. Codon usage data can be employed to estimate the necessary parameters of the calculation, or a priori models of substitution may be employed. Sequences with significant but intermediate constraints on amino acid replacements are, in principle, unanalyzable.      

Interaction of silent and replacement changes in eukaryotic coding sequences

We examined the codon usages in well-conserved and less-well-conserved regions of vertebrate protein genes and found them to be similar. Despite this similarity, there is a statistically significant decrease in codon bias in the less-well-conserved regions. Our analysis suggests that although those codon changes initially fixed under amino acid replacements tend to follow the overall codon usage pattern, they also reduce the bias in codon usage. This decrease in codon bias leads one to predict that the rate of change of synonymous codons should be greater in those regions that are less well conserved at the amino acid level than in the better-conserved regions. Our analysis supports this prediction. Furthermore, we demonstrate a significantly elevated rate of change of synonymous codons among the adjacent codons 5' to amino acid replacement positions. This provides further support for the idea that there are contextual constraints on the choice of synonymous codons in eukaryotes.     

Functional Nonequivalence of Drosophila Actin Isoforms

We show that different Drosophila actinisoforms are not interchangeable. We sequenced the sixgenes that encode conventional Drosophilaactins and found that they specify amino acidreplacements in 27 of 376 positions. To test the significance ofthese changes we used directed mutagenesis to introduce10 such conversions, independently, into the Act88Fflight muscle-specific actin gene. We challenged these variant actins to replace the nativeprotein by transforming germline chromosomes of aDrosophila strain lacking flight muscle actin.Only one of the 10 reproducibly perturbed myofibrillarfunction, demonstrating that most isoform-specific aminoacid replacements are of minor significance. In order toestablish the consequences of multiple amino acidreplacements, we substituted portions of theDrosophila Act88F actin gene with correspondingregions of genes encoding other isoforms. Only one offive constructs tested engendered normally functioningflight muscles, and the severity of myofibrillar defects correlated with the number of replacementswithin the chimeric genes. Finally, we completelyconverted the flight muscle actin-encoding gene to onespecifying a nonmuscle isoform, a change entailing atotal of 18 amino acid replacements. Transformationof flies with this construct resulted in disruption offlight muscle structure and function. We conclude thatactin isoform sequences are not equivalent and that effects of the amino acid replacements,while minor individually, collectively confer uniqueproperties.     

Porcine thyrotropin. The amino-acid sequence of the alpha and beta subunits.

The amino acid sequences of both the alpha and beta subunits of porcine thyrotropin have been studied. Bovine thyrotropin primary structure was taken as a model for ordering the tryptic peptides of porcine thyrotropin. The amino acid sequence of the alpha subunit is identical to that of porcine luteinizing hormone, while oligosaccharide side-chains differ in composition. The primary structure of the beta subunit differs from that of bovine thyrotropin by six amino acid replacements, in positions 22, 24, 26, 36, 62 and 69, and by the absence of a methionyl residue at the carboxy terminus. Chemical evolutions of thyrotropin and luteinizing hormone are compared.     

Homoplasy in genome-wide analysis of rare amino acid replacements: the molecular-evolutionary basis for Vavilov's law of homologous series

Background  

Rare genomic changes (RGCs) that are thought to comprise derived shared characters of individual clades are becoming an increasingly important class of markers in genome-wide phylogenetic studies. Recently, we proposed a new type of RGCs designated RGC_CAMs (after Conserved Amino acids-Multiple substitutions) that were inferred using genome-wide identification of amino acid replacements that were: i) located in unambiguously aligned regions of orthologous genes, ii) shared by two or more taxa in positions that contain a different, conserved amino acid in a much broader range of taxa, and iii) require two or three nucleotide substitutions. When applied to animal phylogeny, the RGC_CAM approach supported the coelomate clade that unites deuterostomes with arthropods as opposed to the ecdysozoan (molting animals) clade. However, a non-negligible level of homoplasy was detected.     

Amino acid replacements in yeast iso-1-cytochrome c. Comparison with the phylogenetic series and the tertiary structure of related cytochromes c

The structural and folding requirements of eukaryotic cytochromes c have been investigated by determining the appropriate DNA sequences of a collection of 46 independent cyc 1 missense mutations obtained in the yeast Saccharomyces cerevisiae and by deducing the corresponding amino acid replacements that abolish function of iso-1-cytochrome c. A total of 33 different replacements at 19 amino acid positions were uncovered in this and previous studies. Because all of these nonfunctional iso-1-cytochromes c are produced at far below the normal level and because a representative number are labile in vitro, most of the replacements appear to be affecting stability of the protein or heme attachment. By considering the tertiary structure of related cytochromes c, the loss of function of most of the mutant iso-1-cytochromes c could be attributed to either replacements of critical residues that directly interact with the heme group or to replacements that disrupt the proper folding of the protein. The replacements of residues interacting with the heme group include those required for covalent attachment (Cys-19 and Cys-22), ligand formation (His-23 and Met-85), and formation of the immediate heme environment (Leu-37, Tyr-53, Trp-64, and Leu-73). Proper folding of the protein is prevented by replacements of glycine residues at sites that cannot accommodate side chains (Gly-11 and Gly-34); by replacements of residues with proline, which limit the torsion angle (Leu-14 and His-38); and by replacements apparently unable to direct the local folding of the backbone into the proper conformation (Pro-35, Tyr-72, Asn-75, Pro-76, Lys-84, Leu-99, and Leu-103). Even though most of the missense mutations occurred at sites corresponding to evolutionarily invariant or conserved residues, a consideration of the replacements in functional revertants indicates that the requirement for residues evolutionarily preserved is less stringent than commonly assumed.