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.     

Organization and structure of Volvox beta-tubulin genes

Summary Genomic clones encoding two Volvox -tubulin genes have been isolated and shown to represent the only two -tubulin genes in the genome. Restriction fragment length polymorphism analysis was used to demonstrate that the two genes are genetically linked. One of these genes was sequenced and the mRNA start site(s) determined by primer extension. A comparison of its sequence to those of the two -tubulin genes of Chlamydomonas revealed: (1) a high degree of conservation of the coding region, with the predicted amino acid sequence differing only in the C-terminal residue; (2) extensive sequence conservation in the 5 untranslated leader region and a 16 bp (putative regulatory) sequence in the promoter region; (3) the same number and location of introns, with a short region of homology in intron 1, but little significant homology in introns 2 and 3.     

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.     

Evolutionary aspects of trypanosomes: Analysis of genes

Summary The genes for four glycolytic enzymes ofTrypanosoma brucei have been analyzed. The proteins encoded by these genes show 38–57% identity with their counterparts in other organisms, whether pro- or eukaryotic. These data are consistent with a phylogenetic tree in which trypanosomes diverged very early from the main branch of the eukaryotic lineage. No definite conclusion can be drawn yet about the evolutionary origin of glycosomes, the microbodies of trypanosomes which contain most enzymes of the glycolytic pathway. A bias could be observed in the codon usage of the glycolytic genes and genes for other housekeeping proteins, indicating that trypanosomes may have selected a nucleotide sequence that enables efficient translation. However, the genes for variant surface glycoproteins (VSGs) do not show such a bias. This lack of preference for special codons is explained by the high evolutionary rate that could be observed for VSG genes.Presented at the FEBS Symposium on Genome Organization and Evolution, held in Crete, Greece, September 1–5, 1986     

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.