Lactic acid bacteria as prime candidates for codon optimization

In species having a strong correlation of expressivity and codon bias it has been shown that heterologous expression can be optimized by changing codons of the introduced gene towards the set of codons that the host organism naturally uses in its highly expressed genes. Even though two lactic acid bacteria are fully sequenced, there are no reports on attempts of codon optimization in the literature. In this report it is demonstrated that codons used in highly expressed genes tend to differ from the codons in lowly expressed genes, and that there is a strong correlation of codon bias and empirical expressivity (codon adaptation index) in Lactococcus lactis and Lactobacillus plantarum. This strongly suggests that codon optimization strategies could be applied to expression systems with lactic acid bacteria as producer strains. A good example of a candidate for codon optimization is the mouse interleukin-2 gene, which in its natural form has an extremely low codon adaptation index for expression in Lc. lactis.     

Information gain for genetic parameter estimation with incorporation of marker data

Genetic marker data has been increasingly incorporated into segregation analysis, as combined segregation and linkage analysis has been performed more frequently. In this article, we study the extent of information gains with incorporation of marker data in segregation analysis, a topic that has not been investigated rigorously. Specifically, the current study is to investigate the influence of marker data on genetic model parameter estimation. A variance matrix criterion (as the inverse of the Fisher information matrix) and a relative entropy criterion (a measure of flatness of expected log-likelihood surface) are used to quantify the information gains. Our results indicate that substantial information gain can be achieved with the incorporation of marker data. The amount of variance reduction increases as the heterozygosity of the linked marker increases and as the trait gets closer to the linked marker(s). Incorporation of marker data in larger pedigrees also yields greater information gains based on both criteria. The effect of pedigree structure is also studied.     

Comprehensive sequence analysis of translation termination sites in various eukaryotes

Recent investigations into the translation termination sites of various organisms have revealed that not only stop codons but also sequences around stop codons have an effect on translation termination. To investigate the relationship between these sequence patterns and translation as well as its termination efficiency, we analysed the correlation between strength of consensus and translation efficiency, as predicted according to Codon Adaptation Index (CAI) value. We used RIKEN full-length mouse cDNA sequences and ten other eukaryotic UniGene datasets from NCBI for the analyses. First, we conducted sequence profile analyses following translation termination sites. We found base G and A at position +1 as a strong consensus for mouse cDNA. A similar consensus was found for other mammals, such as Homo sapiens, Rattus norvegicus and Bos taurus. However, some plants had different consensus sequences. We then analysed the correlation between the strength of consensus at each position and the codon biases of whole coding regions, using information content and CAI value. The results showed that in mouse cDNA, CAI value had a positive correlation with information content at positions +1. We also found that, for positions with strong consensus, the strength of the consensus is likely to have a positive correlation with CAI value in some other eukaryotes. Along with these observations, biological insights into the relationship between gene expression level, codon biases and consensus sequence around stop codons will be discussed.     

The Impact of Message Mutation on the Fitness of a Genetic Code

The Standard Genetic Code is organized such that similar codons encode similar amino acids. One explanation suggested that the Standard Code is the result of natural selection to reduce the fitness ``load' that derives from the mutation and mistranslation of protein-coding genes. We review the arguments against the mutational load-minimizing hypothesis and argue that they need to be reassessed. We review recent analyses of the organization of the Standard Code and conclude that under cautious interpretation they support the mutational load-minimizing hypothesis. We then present a deterministic asexual model with which we study the mode of selection for load minimization. In this model, individual fitness is determined by a protein phenotype resulting from the translation of a mutable set of protein-coding genes. We show that an equilibrium fitness may be associated with a population with the same genetic code and that genetic codes that assign similar codons to similar amino acids have a higher fitness. We also show that the number of mutant codons in each individual at equilibrium, which determines the strength of selection for load minimization, reflects a long-term evolutionary balance between mutations in messages and selection on proteins, rather than the number of mutations that occur in a single generation, as has been assumed by previous authors. We thereby establish that selection for mutational load minimization acts at the level of an individual in a single generation. We conclude with comments on the shortcomings and advantages of load minimization over other hypotheses for the origin of the Standard Code. Received: 4 April 2001 / Accepted: 22 October 2001     

Trends in Codon and Amino Acid Usage in Thermotoga maritima

The usage of synonymous codons and the frequencies of amino acids were investigated in the complete genome of the bacterium Thermotoga maritima using a multivariate statistical approach. The GC3 content of each gene was the most prominent source of variation of codon usage. Surprisingly the usage of UGU and UGC (synonymous triplets coding for Cys, the least frequent amino acid in this species) was detected as the second most prominent source of variation. However, this result is probably an artifact due to the very low frequency of Cys together with the nonbiased composition of this genome. The third trend was related to the preferential usage of a subset of codons among highly expressed genes, and these triplets are presumed to be translationally optimal. Concerning the amino acid usage, the hydropathy level of each protein (and therefore the frequency of charged residues) was the main trend, while the second factor was related to the frequency of usage of the smaller residues, suggesting that the cell economy strongly influences the architecture of the proteins. The third axis of the analysis discriminated the usage of Phe, Tyr, Trp (aromatic residues) plus Cys, Met, and His. These six residues have in common the property of being the preferential targets of reactive oxygen species, and therefore the anaerobic condition of T. maritima is an important factor for the amino acid frequencies. Finally, the Cys content of each protein was the fourth trend. Received: 22 June 2001 / Accepted: 1 October 2001     

Comparative Molecular Evolution of Primary (<Emphasis Type="BoldItalic">Buchnera</Emphasis>) and Secondary Symbionts of Aphids Based on Two Protein-Coding Genes

A+T content, phylogenetic relationships, codon usage, evolutionary rates, and ratio of synonymous versus non-synonymous substitutions have been studied in partial sequences of the atpD and aroQ/pheA genes of primary (Buchnera) and secondary symbionts of aphids and a set of selected non-symbiotic bacteria, belonging to the five subdivisions of the Proteobacteria. Compared to the homologous genes of the last group, both genes belonging to Buchnera behave in a similar way, showing a higher A+T content, forming a monophyletic group, a loss in codon bias, especially in third base position, an evolutionary acceleration and an increase in the number of non-synonymous substitutions, confirming previous results reported elsewhere for other genes. When available, these properties have been partly observed with the secondary symbionts, but with values that are intermediate between Buchnera and free living Proteobacteria. They show high A+T content, but not as high as Buchnera, a non-solved phylogenetic position between Buchnera, and the other γ-Proteobacteria, a loss in codon bias, again not as high as in Buchnera and a significant evolutionary acceleration in the case of the three atpD genes, but not when considering aroQ/pheA genes. These results give support to the hypothesis that they are symbionts at different stages of the symbiotic accommodation to the host.     

Evaluation of 17alpha-E-(trifluoromethylphenyl)vinyl estradiols as novel estrogen receptor ligands

As part of our program to develop novel ligands for the estrogen receptor, we synthesized the series of isomeric 17alpha-(trifluoromethyl)phenylvinyl estradiols using our solid-phase organic synthesis methodology. The compounds were evaluated for their relative binding affinity (RBA) using the ERalpha-LBD and in vivo potency using the immature rat uterotrophic growth assay. The ortho-isomer had the highest RBA values, 48-223, and the highest estrogenicity in vivo. The other isomers had significantly lower affinities and were weaker agonists in the uterotrophic assay. The results suggest that introduction of substituents at the 17alpha-position of estradiol is tolerated by the ER-LBD and permit agonist responses in the intact animal, however, the effect is sensitive to the position of groups on the phenyl ring. This study demonstrates that the 17alpha-position of estradiol is a reasonable site for modification but the position and physicochemical properties of such modifications may significantly affect the affinity and efficacy of the ligand.     

Cranial integration in Homo: singular warps analysis of the midsagittal plane in ontogeny and evolution

This study addresses some enduring issues of ontogenetic and evolutionary integration in the form of the hominid cranium. Our sample consists of 38 crania: 20 modern adult Homo sapiens, 14 sub-adult H. sapiens, and four archaic Homo. All specimens were CT-scanned except for two infant H. sapiens, who were imaged by MR instead. For each specimen 84 landmarks and semi-landmarks were located on the midsagittal plane and converted to Procrustes shape coordinates. Integration was quantified by the method of singular warps, a new geometric-statistical approach to visualizing correlations among regions. The two classic patterns of integration, evolutionary and ontogenetic, were jointly explored by comparing analyses of overlapping subsamples that span ranges of different hypothetical factors. Evolutionary integration is expressed in the subsample of 24 adult Homo, and ontogenetic integration in the subsample of 34 H. sapiens. In this data set, vault, cranial base, and face show striking and localized patterns of covariation over ontogeny, similar but not identical to the patterns seen over evolution. The principal differences between ontogeny and phylogeny pertain to the cranial base. There is also a component of cranial length to height ratio not reducible to either process. Our methodology allows a separation of these independent processes (and their impact on cranial shape) that conventional methods have not found.     

Insights on bias and information in group-level studies

Ecological and aggregate data studies are examples of group-level studies. Even though the link between the predictors and outcomes is not preserved in these studies, inference about individual-level exposure effects is often a goal. The disconnection between the level of inference and the level of analysis expands the array of potential biases that can invalidate the inference from group-level studies. While several sources of bias, specifically due to measurement error and confounding, may be more complex in group-level studies, two sources of bias, cross-level and model specification bias, are a direct consequence of the disconnection. With the goal of aligning inference from individual versus group-level studies, I discuss the interplay between exposure and study design. I specify the additional assumptions necessary for valid inference, specifically that the between- and within-group exposure effects are equal. Then cross-level inference is possible. However, all the information in the group-level analysis comes from between-group comparisons. Models where the group-level analysis provides even a small percentage of information about the within-group exposure effect are most susceptible to model specification bias. Model specification bias can be even more serious when the group-level model isn't derived from an individual-level model.     

Small-sample bias and corrections for conditional maximum-likelihood odds-ratio estimators

A number of small-sample corrections have been proposed for the conditional maximum-likelihood estimator of the odds ratio for matched pairs with a dichotomous exposure. I here contrast the rationale and performance of several corrections, specifically those that generalize easily to multiple conditional logistic regression. These corrections or Bayesian analyses with informative priors may serve as diagnostics for small-sample problems. Points are illustrated with a small exact performance comparison and with an example from a study of electrical wiring and childhood leukemia. The former comparison suggests that small-sample bias may be more prevalent than commonly realized.