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.     

An Efficient Regression Type Estimator in Survey Sampling

For the estimation of population mean in simple random sampling, an efficient regression-type estimator is proposed which is more efficient than the conventional regression estimator and hence than mean per unit estimator, ratio and product estimators and many other estimators proposed by various authors. Some numerical examples are included for illustration.     

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.     

Some evidence of the enzymatic conversion of bovine suppressor phosphoseryl-tRNA to selenocysteyl-tRNA

In order to clarify the mechanisms of selenocysteine incorporation into glutathione peroxidase, some evidence to show the in vitro conversion of phosphoseryl-tRNA to selenocysteyl-tRNA is reported. [³H]Phosphoseryl-tRNA was incubated in a reaction mixture composed of SeO₂, glutathione and NADPH in the presence of selenium-transferase partially purified. Analyses of amino acids on the product tRNA showed that a part (4%) of [³H]phosphoseryl-tRNA was changed to [³H]selenocysteyl-tRNA. The conversion from seryl-tRNA^su or major seryl-tRNA_IGA was not found. Selenium-transferase was essential for the conversion. [³H]Selenocysteine, liberated from the tRNA, was modified with iodoacetic acid. The product was confirmed to be carboxymethyl-selenocysteine by two-dimensional TLC. Selenocysteyl-tRNA^su should be used to synthesize glutathione peroxidase by co-translational mechanisms.     

Correction of a bootstrap approach to testing for evolution along lines of least resistance

Testing for an association between the leading vectors of multivariate trait (co)variation within populations (the ‘line of least resistance’) and among populations is an important tool for exploring variational bias in evolution. In a recent study of stickleback fish populations, a bootstrap‐based test was introduced that takes into account estimation error in both vectors and hence improves the previously available bootstrap method. Because this test was implemented incorrectly, however, I here describe the correct test protocol and provide a reanalysis of the original data set. The application of this new test protocol should improve future investigations of evolution along lines of least resistance and other vector comparisons.     

Conditional Maximum Likelihood Estimation Following a Group Sequential Test

We consider estimation after a group sequential test. An estimator that is unbiased or has small bias may have substantial conditional bias (Troendle and Yu, 1999, Coburger and Wassmer, 2001). In this paper we derive the conditional maximum likelihood estimators of both the primary parameter and a secondary parameter, and investigate their properties within a conditional inference framework. The method applies to both the usual and adaptive group sequential test designs. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

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.     

Interaction of silent and replacement changes in eukaryotic coding sequences

Summary We examined the codon usages in wellconserved 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.     

Sampling bias in estimating Design II variance components with S1 families

Summary The use of several S₁ individuals to represent an S₀ individual permits the use of a Design II mating scheme for plants with only one pistillate flower per plant. Estimates of additive (V _A ) and dominance (V _D ) variance from this mating scheme will be biased upwards, when a small number (10) of individuals of each S₁ line are used. This bias can be computed, and the additive and dominance estimates can be corrected. Of particular interest is the observation that the additive genetic variance contributes to bias in estimates of V _D . When S₀ plants are non inbred and their selfedprogeny (S₁ lines) are used to represent them in developing families for use in the Design II, where m₁ is the number of individuals used to represent an S₁ line in developing half sib-families and m₂ is the number of individuals used to represent the S₁ line in making up full sib-families. For example, in a 3×3 Design II, with about 10 individuals used to represent each S₁ line in each cross, m₂ = 10 and m₁ = 30. When m₁ = m₂ = 1, and Joint contribution from Department of Agronomy, University of Nebraska 68583, and the S. S. Cameron Laboratory, Werribee, Victoria 3030, Australia. Published as paper No. 7395, Journal Series