Contrast-FEL—A Test for Differences in Selective Pressures at Individual Sites among Clades and Sets of Branches

A number of evolutionary hypotheses can be tested by comparing selective pressures among sets of branches in a phylogenetic tree. When the question of interest is to identify specific sites within genes that may be evolving differently, a common approach is to perform separate analyses on subsets of sequences and compare parameter estimates in a post hoc fashion. This approach is statistically suboptimal and not always applicable. Here, we develop a simple extension of a popular fixed effects likelihood method in the context of codon-based evolutionary phylogenetic maximum likelihood testing, Contrast-FEL. It is suitable for identifying individual alignment sites where any among the K≥2 sets of branches in a phylogenetic tree have detectably different ω ratios, indicative of different selective regimes. Using extensive simulations, we show that Contrast-FEL delivers good power, exceeding 90% for sufficiently large differences, while maintaining tight control over false positive rates, when the model is correctly specified. We conclude by applying Contrast-FEL to data from five previously published studies spanning a diverse range of organisms and focusing on different evolutionary questions.     

Identification of the coding region for the putidaredoxin reductase gene from the plasmid of Pseudomonas putida

The first 12 NH₂-terminal amino acids of the Pseudomonas putida putidaredoxin reductase were shown to be Met-Asn-Ala-Asn-Asp-Asn-Val-Val-Ile-Val-Gly-Thr. Comparison of these data with the DNA sequence of the BamHI-HindIII 197-base fragment derived from the PstI 2.2-kb fragment obtained from the P. putida plasmid showed that the putidaredoxin reductase gene was downstream from the cytochrome P-450 gene and the intergenic region had the 24-nucleotide sequence TAAACACATGGGAGTGCGTGCTAA. The Shine-Dalgarno sequence GGAG was detected in this region. The initiating triplet for the reductase gene was GTG, which normally codes for valine, but in the initiating codon position codes for methionine. From the amino acid sequence and X-ray data comparisons with other flavoproteins, what appears to be the AMP binding region of the FAD can be recognized in the NH₂-terminal portion of the reductase involving residues 5–35.This article was presented during the proceedings of the International Conference on Macromolecular Structure and Function, held at the National Defence Medical College, Tokorozawa, Japan, December 1985.     

Nonuniform size distribution of nascent globin peptides, evidence for pause localization sites, and a cotranslational protein-folding model

Examination of nascent globin peptides accumulatingin vitro during globin synthesis in rabbit reticulocyte lysates was carried out. A view was supported that nonrandom distribution of codons with different usage frequencies in mRNA may determine the messenger's translation kinetics. Regions of reduced translation of - and -globin polypeptide chains were localized, and the cotranslational protein-folding model suggested previously was substantiated. An active conjunction of synthesis and folding of proteins was proposed as one of the main destinations of a translation nonuniformity.     

Sequence analysis of a Dictyostelium discoideum gene coding for an active dihydroorotate dehydrogenase in yeast

A Dictyostelium discoideum DNA fragment isolated on the basis of its ability to complement the ural mutation of yeast, codes for a dihydroorotate dehydrogenase activity. The complete nucleotide sequence of this 1898 bp fragment has been determined and reveals an open reading frame capable of coding for a 369 amino acid polypeptide of molecular mass 47.000. The gene shows preferential use of codons with weak pairing forces. Eleven codons, mainly those with a G in the third position, are absent. The flanking sequences are unusually rich in A + T (80%). Several direct and inverted repeats exist in the 5' flanking sequence.     

The araBAD operon of Salmonella typhimurium LT2. I. Nucleotide sequence of araB and primary structure of its product, ribulokinase

Hybrid plasmids containing the araBAD operon of Salmonella typhimurium LT2 were characterized by Southern blot and genetic analyses. The nucleotide sequence of araB was determined. The araB gene product, ribulokinase (EC 2.7.1.16), was purified and the results of amino acid composition analysis and partial amino acid sequence are in agreement with predictions from the DNA sequence. Ribulokinase is 569 amino acid residues long and has a calculated Mr of 61 793. Ribulokinase shares significant homology with xylulose kinase from Escherichia coli. Codon usage in the araB gene does not favor those codons which have intermediate codon-anticodon binding energy.     

Mitogenic response to epidermal growth factor (EGF) modulated by platelet-derived growth factor in cultured fibroblasts

The mitotic effects of epidermal growth factor (EGF) were investigated in two cultured fibroblast lines, BALB/c-3T3 and C3H 10T1/2 cells. EGF (30 ng/ml) added to quiescent 3T3 cells in medium containing either platelet-poor plasma or 10(-5) M insulin caused only minimal increases in the percentage of cells stimulated to initiate DNA synthesis. In contrast, EGF acted synergistically with either insulin or plasma to stimulate DNA synthesis in quiescent cultures of 10T1/2 cells, although the maximum effects of EGF were measured at concentrations several-fold greater than those found in either serum or plasma. In either 3T3 or 10T1/2 cells a transient preexposure to platelet-derived growth factor (PDGF) caused over a 10-fold increase in the sensitivity to the mitogenic effects of EGF. It is therefore possible that a primary action of PDGF is to increase the sensitivity of fibroblasts to EGF, independent of whether EGF alone is found to be mitogenic.     

The translational stop signal: Codon with a context, or extended factor recognition element?

Wide ranging studies of the readthrough of translational stop codons within the last 25 years have suggested that the stop codon might be only part of the molecular signature for recognition of the termination signal. Such studies do not distinguish between effects on suppression and effects on termination, and so we have used a number of different approaches to deduce whether the stop signal is a codon with a context or an extended factor recognition element. A data base of natural termination sites from a wide range of organisms (148 organisms, 40000 sequences) shows a very marked bias in the bases surrounding the stop codon in the genes for all organisms examined, with the most dramatic bias in the base following the codon (+4). The nature of this base determines the efficiency of the stop signal in vivo, and in Escherichia coli this is reinforced by overexpressing the stimulatory factor, release factor-3. Strong signals, defined by their high relative rates of selecting the decoding release factors, are enhanced whereas weak signals respond relatively poorly. Site-directed cross-linking from the +1, and bases up to +6 but not beyond make close contact with the bacterial release factor-2. The translational stop signal is deduced to be an extended factor recognition sequence with a core element, rather than simply a factor recognition triplet codon influenced by context.     

Codons AGA and AGG are read as glycine in ascidian mitochondria

Summary A 1.2-kb DNA fragment of the cytochrome oxidase subunit I (CO I) gene of mitochondria isolated from an ascidian,Halocynthia roretzi, was amplified by polymerase chain reaction (PCR) and sequenced. Codons AGA and AGG appeared in its reading frame, indicating that these are sense codons in this organelle. Sequence comparisons with the corresponding regions of other animal mitochondrial CO I genes suggest that codons AGA and AGG correspond to glycine in the ascidian mitochondrial genome, but not to serine as in most invertebrate genomes, nor to stops as in vertebrate genomes. The other codons are identical to those of vertebrate mitochondria.