Nonsynonymous substitution in abalone sperm fertilization genes exceeds substitution in introns and mitochondrial DNA

Strong positive Darwinian selection acts on two sperm fertilization proteins, lysin and 18-kDa protein, from abalone (Haliotis). To understand the phylogenetic context for this dramatic molecular evolution, we obtained sequences of mitochondrial cytochrome c oxidase subunit I (mtCOI), and genomic sequences of lysin, 18-kDa, and a G protein subunit. Based on mtDNA differentiation, four north Pacific abalone species diverged within the past 2 million years (Myr), and remaining north Pacific species diverged over a period of 4–20 Myr. Between-species nonsynonymous differences in lysin and 18-kDa exons exceed nucleotide differences in introns by 3.5- to 24-fold. Remarkably, in some comparisons nonsynonymous substitutions in lysin and 18-kDa genes exceed synonymous substitutions in mtCOI. Lysin and 18-kDa intron/exon segments were sequenced from multiple red abalone individuals collected over a 1,200-km range. Only two nucleotide changes and two sites of slippage variation were detected in a total of >29,000 nucleotides surveyed. However, polymorphism in mtCOI and a G protein intron was found in this species. This finding suggests that positive selection swept one lysin allele and one 18-kDa allele to fixation. Similarities between mtCOI and lysin gene trees indicate that rapid adaptive evolution of lysin has occurred consistently through the history of the group. Comparisons with mtCOI molecular clock calibrations suggest that nonsynonymous substitutions accumulate 2–50 times faster in lysin and 18-kDa genes than in rapidly evolving mammalian genes.

Cloning and functional expression of a cDNA encoding a pheromone gland-specific acyl-CoA Δ11-desaturase of the cabbage looper moth, Trichoplusia ni

Desaturation of coenzyme-A esters of saturated fatty acids is a common feature of sex pheromone biosynthetic pathways in the Lepidoptera. The enzymes that catalyze this step share several biochemical properties with the ubiquitous acyl-CoA Δ9-desaturases of animals and fungi, suggesting a common ancestral origin. Unlike metabolic acyl-CoA Δ9-desaturases, pheromone desaturases have evolved unusual regio- and stereoselective activities that contribute to the remarkable diversity of chemical structures used as pheromones in this large taxonomic group. In this report, we describe the isolation of a cDNA encoding a pheromone gland desaturase from the cabbage looper moth, Trichoplusia ni, a species in which all unsaturated pheromone products are produced via a Δ11Z-desaturation mechanism. The largest ORF of the ≈1,250-bp cDNA encodes a 349-aa apoprotein (PDesat-Tn Δ11Z) with a predicted molecular mass of 40,240 Da. Its hydrophobicity profile is similar overall to those of rat and yeast Δ9-desaturases, suggesting conserved transmembrane topology. A 182-aa core domain delimited by conserved histidine-rich motifs implicated in iron-binding and catalysis has 72 and 58% similarity (including conservative substitutions) to acyl-CoA Δ9Z-desaturases of rat and yeast, respectively. Northern blot analysis revealed an ≈1,250-nt PDesat-Tn Δ11Z mRNA that is consistent with the spatial and temporal distribution of Δ11-desaturase enzyme activity. Genetic transformation of a desaturase-deficient strain of the yeast Saccharomyces cerevisiae with an expression plasmid encoding PDesat-Tn Δ11Z resulted in complementation of the strain’s fatty acid auxotrophy and the production of Δ11Z-unsaturated fatty acids.

An important 2′-OH group for an RNA–protein interaction

We have investigated the role of 2′-OH groups in the specific interaction between the acceptor stem of Escherichia coli tRNACys and cysteine-tRNA synthetase. This interaction provides for the high aminoacylation specificity observed for cysteine-tRNA synthetase. A synthetic RNA microhelix that recapitulates the sequence of the acceptor stem was used as a substrate and variants containing systematic replacement of the 2′-OH by 2′-deoxy or 2′-O-methyl groups were tested. Except for position U73, all substitutions had little effect on aminoacylation. Interestingly, the deoxy substitution at position U73 had no effect on aminoacylation, but the 2′-O-methyl substitution decreased aminoacylation by 10-fold and addition of the even bulkier 2′-O-propyl group decreased aminoacylation by another 2-fold. The lack of an effect by the deoxy substitution suggests that the hydrogen bonding potential of the 2′-OH at position U73 is unimportant for aminoacylation. The decrease in activity upon alkyl substitution suggests that the 2′-OH group instead provides a monitor of the steric environment during the RNA–synthetase interaction. The steric role was confirmed in the context of a reconstituted tRNA and is consistent with the observation that the U73 base is the single most important determinant for aminoacylation and therefore is a site that is likely to be in close contact with cysteine-tRNA synthetase. A steric role is supported by an NMR-based structural model of the acceptor stem, together with biochemical studies of a closely related microhelix. This role suggests that the U73 binding site for cysteine-tRNA synthetase is sterically optimized to accommodate a 2′-OH group in the backbone, but that the hydroxyl group itself is not involved in specific hydrogen bonding interactions.

Highly conserved features of DNA binding between two divergent members of the Myb family of transcription factors

Bas1p, a divergent yeast member of the Myb family of transcription factors, shares with the proteins of this family a highly conserved cysteine residue proposed to play a role in redox regulation. Substitutions of this residue in Bas1p (C153) allowed us to establish that, despite its very high conservation, it is not strictly required for Bas1p function: its substitution with a small hydrophobic residue led to a fully functional protein in vitro and in vivo. C153 was accessible to an alkylating agent in the free protein but was protected by prior exposure to DNA. The reactivity of cysteines in the first and third repeats was much lower than in the second repeat, suggesting a more accessible conformation of repeat 2. Proteolysis protection, fluorescence quenching and circular dichroism experiments further indicated that DNA binding induces structural changes making Bas1p less accessible to modifying agents. Altogether, our results strongly suggest that the second repeat of the DNA-binding domain of Bas1p behaves similarly to its Myb counterpart, i.e. a DNA-induced conformational change in the second repeat leads to formation of a full helix–turn–helix-related motif with the cysteine packed in the hydrophobic core of the repeat.

Induction and repression of DAN1 and the family of anaerobic mannoprotein genes in Saccharomyces cerevisiae occurs through a complex array of regulatory sites

The DAN/TIR mannoprotein genes of Saccharomyces cerevisiae (DAN1, DAN2, DAN3, DAN4, TIR1, TIR2, TIR3 and TIR4) are expressed in anaerobic cells while the predominant cell wall proteins Cwp1 and Cwp2 are down-regulated. Elements involved in activation and repression of the DAN/TIR genes were defined in this study, using the DAN1 promoter as a model. Nested deletions in a DAN1/lacZ reporter pinpointed regions carrying activation and repression elements. Inspection revealed two consensus sequences subsequently shown to be independent anaerobic response elements (AR1, consensus TCGTTYAG; AR2, consensus AAAAATTGTTGA). AR1 is found in all of the DAN/TIR promoters; AR2 is found in DAN1, DAN2 and DAN3. A 120 bp segment carrying two copies of AR1 preferentially activated transcription of lacZ under anaerobic conditions. A fusion of three synthetic copies of AR1 to MEL1 was also expressed anaerobically. Mutations in either AR1 site within the 120 bp segment caused a drastic loss of expression, indicating that both are necessary for activation and implying cooperativity between adjacent transcriptional activation complexes. A single AR2 site carried on a 46 bp fragment from the DAN1 promoter activated lacZ transcription under anaerobic conditions, as did a 26 bp synthetic AR2 fragment fused to MEL1. Nucleotide substitutions within the AR2 sequence eliminated the activity of the 46 bp segment. Ablation of the AR2 sequences in the full promoter caused a partial reduction of expression. The presence of the ATTGTT core (recognized by HMG proteins) in the AR2 sequence suggests that an HMG protein may activate through AR2. One region was implicated in aerobic repression of DAN1. It contains sites for the heme-induced Mot3 and Rox1 repressors.

Mutagenesis of the peptidyltransferase center of 23S rRNA: the invariant U2449 is dispensable

U2449 is one of many invariant residues in the central loop of domain V of 23S rRNA, a region that constitutes part of the peptidyltransferase center of the ribosome. In Escherichia coli, this U is post-transcriptionally modified to dihydrouridine (D) and is the only D modification found in E.coli rRNAs. To analyze the role of this base and its modification in ribosomal function, all three base substitutions were constructed on a plasmid copy of the rrnB operon and assayed for their ability to support cell growth in a strain of E.coli lacking chromosomal rrn operons. Both purine substitution mutations were not viable. However, growth and antibiotic sensitivity of cells expressing only the mutant D2449C rRNA was indistinguishable from wild type. We conclude that while a pyrimidine is required at position 2449 for proper ribosomal function, the D modification is dispensable.

The Catastrophe-promoting Activity of Ectopic Op18/Stathmin Is Required for Disruption of Mitotic Spindles But Not Interphase Microtubules

Oncoprotein18/stathmin (Op18) is a microtubule (MT) destabilizing protein that is inactivated during mitosis by phosphorylation at four Ser-residues. Op18 has at least two functions; the N-terminal region is required for catastrophe-promotion (i.e., transition from elongation to shortening), while the C-terminal region is required to inhibit MT-polymerization rate in vitro. We show here that a “pseudophosphorylation” derivative of Op18 (i.e., four Ser- to Glu-substitutions at phosphorylation sites) exhibits a selective loss of catastrophe-promoting activity. This is contrasted to authentic phosphorylation, which efficiently attenuates all activities except tubulin binding. In intact cells, overexpression of pseudophosphorylated Op18, which is not phosphorylated by endogenous kinases, is shown to destabilize interphase MTs but to leave spindle formation untouched. To test if the mitotic spindle is sensitive only to the catastrophe-promoting activity of Op18 and resistant to C-terminally associated activities, N- and C-terminal truncations with defined activity-profiles were employed. The cell-cycle phenotypes of nonphosphorylatable mutants (i.e., four Ser- to Ala-substitutions) of these truncation derivatives demonstrated that catastrophe promotion is required for interference with the mitotic spindle, while the C-terminally associated activities are sufficient to destabilize interphase MTs. These results demonstrate that specific Op18 derivatives with defined activity-profiles can be used as probes to distinguish interphase and mitotic MTs.

Reverse engineering the (β/α)8 barrel fold

The (β/α)8 barrel is the most commonly occurring fold among protein catalysts. To lay a groundwork for engineering novel barrel proteins, we investigated the amino acid sequence restrictions at 182 structural positions of the prototypical (β/α)8 barrel enzyme triosephosphate isomerase. Using combinatorial mutagenesis and functional selection, we find that turn sequences, α-helix capping and stop motifs, and residues that pack the interface between β-strands and α-helices are highly mutable. Conversely, any mutation of residues in the central core of the β-barrel, β-strand stop motifs, and a single buried salt bridge between amino acids R189 and D227 substantially reduces catalytic activity. Four positions are effectively immutable: conservative single substitutions at these four positions prevent the mutant protein from complementing a triosephosphate isomerase knockout in Escherichia coli. At 142 of the 182 positions, mutation to at least one amino acid of a seven-letter amino acid alphabet produces a triosephosphate isomerase with wild-type activity. Consequently, it seems likely that (β/α)8 barrel structures can be encoded with a subset of the 20 amino acids. Such simplification would greatly decrease the computational burden of (β/α)8 barrel design.

Insulin inhibits transcription of IRS-2 gene in rat liver through an insulin response element (IRE) that resembles IREs of other insulin-repressed genes

Recent data indicate that sustained elevations in plasma insulin suppress the mRNA for IRS-2, a component of the insulin signaling pathway in liver, and that this deficiency contributes to hepatic insulin resistance and inappropriate gluconeogenesis. Here, we use nuclear run-on assays to show that insulin inhibits transcription of the IRS-2 gene in the livers of intact rats. Insulin also inhibited transcription of a reporter gene driven by the human IRS-2 promoter that was transfected into freshly isolated rat hepatocytes. The human promoter contains a heptanucleotide sequence, TGTTTTG, that is identical to the insulin response element (IRE) identified previously in the promoters of insulin-repressed genes. Single base pair substitutions in this IRE decreased transcription of the IRS-2-driven reporter in the absence of insulin and abolished insulin-mediated repression. We conclude that insulin represses transcription of the IRS-2 gene by blocking the action of a positive factor that binds to the IRE. Sustained repression of IRS-2, as occurs in chronic hyperinsulinemia, contributes to hepatic insulin resistance and accelerates the development of the diabetic state.

Computational method to reduce the search space for directed protein evolution

We introduce a computational method to optimize the in vitro evolution of proteins. Simulating evolution with a simple model that statistically describes the fitness landscape, we find that beneficial mutations tend to occur at amino acid positions that are tolerant to substitutions, in the limit of small libraries and low mutation rates. We transform this observation into a design strategy by applying mean-field theory to a structure-based computational model to calculate each residue's structural tolerance. Thermostabilizing and activity-increasing mutations accumulated during the experimental directed evolution of subtilisin E and T4 lysozyme are strongly directed to sites identified by using this computational approach. This method can be used to predict positions where mutations are likely to lead to improvement of specific protein properties.

An Allele of the Ripening-Specific 1-Aminocyclopropane-1-Carboxylic Acid Synthase Gene (ACS1) in Apple Fruit with a Long Storage Life1

An allele of the 1-aminocyclopropane-1-carboxylic acid (ACC) synthase gene (Md-ACS1), the transcript and translated product of which have been identified in ripening apples (Malus domestica), was isolated from a genomic library of the apple cultivar, Golden Delicious. The predicted coding region of this allele (ACS1-2) showed that seven nucleotide substitutions in the corresponding region of ACS1-1 resulted in just one amino acid transition. A 162-bp sequence characterized as a short interspersed repetitive element retrotransposon was inserted in the 5′-flanking region of ACS1-2 corresponding to position −781 in ACS1-1. The XhoI site located near the 3′ end of the predicted coding region of ACS1-2 was absent from the reverse transcriptase-polymerase chain reaction product, revealing that exclusive transcription from ACS1-1 occurs during ripening of cv Golden Delicious fruit. DNA gel-blot and polymerase chain reaction analyses of genomic DNAs showed clearly that apple cultivars were either heterozygous for ACS1-1 and ACS1-2 or homozygous for each type. RNA gel-blot analysis of the ACS1-2 homozygous Fuji apple, which produces little ethylene and has a long storage life, demonstrated that the level of transcription from ACS1-2 during the ripening stage was very low.

Use of Ubiquitin Fusions to Augment Protein Expression in Transgenic Plants1

A major goal of plant biotechnology is the production of genetically engineered crops that express natural or foreign proteins at high levels. To enhance protein accumulation in transgenic plants, we developed a set of vectors that express proteins and peptides as C-terminal translational fusions with ubiquitin (UBQ). Studies of several proteins in tobacco (Nicotiana tabacum) showed that: (a) proteins can be readily expressed in plants as UBQ fusions; (b) by the action of endogenous UBQ-specific proteases (Ubps), these fusions are rapidly and precisely processed in vivo to release the fused protein moieties in free forms; (c) the synthesis of a protein as a UBQ fusion can significantly augment its accumulation; (d) proper processing and localization of a protein targeted to either the apoplast or the chloroplast is not affected by the N-terminal UBQ sequence; and (e) single amino acid substitutions surrounding the cleavage site can inhibit in vivo processing of the fusion by Ubps. Noncleavable UBQ fusions of β-glucuronidase became extensively modified, with additional UBQs in planta. Because multiubiquitinated proteins are the preferred substrates of the 26S proteasome, noncleavable fusions may be useful for decreasing protein half-life. Based on their ability to augment protein accumulation and the sequence specificity of Ubps, UBQ fusions offer a versatile way to express plant proteins.

Long term trends in the evolution of H(3) HA1 human influenza type A

We have studied the HA1 domain of 254 human influenza A(H3N2) virus genes for clues that might help identify characteristics of hemagglutinins (HAs) of circulating strains that are predictive of that strain’s epidemic potential. Our preliminary findings include the following. (i) The most parsimonious tree found requires 1,260 substitutions of which 712 are silent and 548 are replacement substitutions. (ii) The HA1 portion of the HA gene is evolving at a rate of 5.7 nucleotide substitutions/year or 5.7 × 10−3 substitutions/site per year. (iii) The replacement substitutions are distributed randomly across the three positions of the codon when allowance is made for the number of ways each codon can change the encoded amino acid. (iv) The replacement substitutions are not distributed randomly over the branches of the tree, there being 2.2 times more changes per tip branch than for non-tip branches. This result is independent of how the virus was amplified (egg grown or kidney cell grown) prior to sequencing or if sequencing was carried out directly on the original clinical specimen by PCR. (v) These excess changes on the tip branches are probably the result of a bias in the choice of strains to sequence and the detection of deleterious mutations that had not yet been removed by negative selection. (vi) There are six hypervariable codons accumulating replacement substitutions at an average rate that is 7.2 times that of the other varied codons. (vii) The number of variable codons in the trunk branches (the winners of the competitive race against the immune system) is 47 ± 5, significantly fewer than in the twigs (90 ± 7), which in turn is significantly fewer variable codons than in tip branches (175 ± 8). (viii) A minimum of one of every 12 branches has nodes at opposite ends representing viruses that reside on different continents. This is, however, no more than would be expected if one were to randomly reassign the continent of origin of the isolates. (ix) Of 99 codons with at least four mutations, 31 have ratios of non-silent to silent changes with probabilities less than 0.05 of occurring by chance, and 14 of those have probabilities <0.005. These observations strongly support positive Darwinian selection. We suggest that the small number of variable positions along the successful trunk lineage, together with knowledge of the codons that have shown positive selection, may provide clues that permit an improved prediction of which strains will cause epidemics and therefore should be used for vaccine production.

The high spontaneous mutation rate: Is it a health risk?*

The human mutation rate for base substitutions is much higher in males than in females and increases with paternal age. This effect is mainly, if not entirely, due to the large number of cell divisions in the male germ line. The mutation-rate increase is considerably greater than expected if the mutation rate were simply proportional to the number of cell divisions. In contrast, those mutations that are small deletions or rearrangements do not show the paternal age effect. The observed increase with the age of the father in the incidence of children with different dominant mutations is variable, presumably the result of different mixtures of base substitutions and deletions. In Drosophila, the rate of mutations causing minor deleterious effects is estimated to be about one new mutation per zygote. Because of a larger number of genes and a much larger amount of DNA, the human rate is presumably higher. Recently, the Drosophila data have been reanalyzed and the mutation-rate estimate questioned, but I believe that the totality of evidence supports the original conclusion. The most reasonable way in which a species can cope with a high mutation rate is by quasi-truncation selection, whereby a number of mutant genes are eliminated by one “genetic death.”

Dynamic evolution of plant mitochondrial genomes: Mobile genes and introns and highly variable mutation rates

We summarize our recent studies showing that angiosperm mitochondrial (mt) genomes have experienced remarkably high rates of gene loss and concomitant transfer to the nucleus and of intron acquisition by horizontal transfer. Moreover, we find substantial lineage-specific variation in rates of these structural mutations and also point mutations. These findings mostly arise from a Southern blot survey of gene and intron distribution in 281 diverse angiosperms. These blots reveal numerous losses of mt ribosomal protein genes but, with one exception, only rare loss of respiratory genes. Some lineages of angiosperms have kept all of their mt ribosomal protein genes whereas others have lost most of them. These many losses appear to reflect remarkably high (and variable) rates of functional transfer of mt ribosomal protein genes to the nucleus in angiosperms. The recent transfer of cox2 to the nucleus in legumes provides both an example of interorganellar gene transfer in action and a starting point for discussion of the roles of mechanistic and selective forces in determining the distribution of genetic labor between organellar and nuclear genomes. Plant mt genomes also acquire sequences by horizontal transfer. A striking example of this is a homing group I intron in the mt cox1 gene. This extraordinarily invasive mobile element has probably been acquired over 1,000 times separately during angiosperm evolution via a recent wave of cross-species horizontal transfers. Finally, whereas all previously examined angiosperm mtDNAs have low rates of synonymous substitutions, mtDNAs of two distantly related angiosperms have highly accelerated substitution rates.