Maximum-likelihood analysis of molecular adaptation in abalone sperm lysin reveals variable selective pressures among lineages and sites

Maximum-likelihood models of codon substitution were used to analyze sperm lysin genes of 25 abalone (HALIOTIS:) species to identify lineages and amino acid sites under diversifying selection. The models used the nonsynonymous/synonymous rate ratio (omega = d(N)/d(S)) as an indicator of selective pressure and allowed the ratio to vary among lineages or sites. Likelihood ratio tests suggested significant variation in selective pressure among lineages. The variable selective pressure provided an explanation for the previous observation that the omega ratio is >1 in comparisons of closely related species and <1 in comparisons of distantly related species. Computer simulations demonstrated that saturation of nonsynonymous substitutions and constraint on lysin structure were unlikely to account for the observed pattern. Lineages linking closely related sympatric species appeared to be under diversifying selection, while lineages separating distantly related species from different geographic locations were associated with low evolutionary rates. The selective pressure indicated by the omega ratio was found to vary greatly among amino acid sites in lysin. Sites under potential diversifying selection were identified. Ancestral lysins were inferred to trace the route of evolution at individual sites and to provide lysin sequences for future laboratory studies.     

Detecting the form of selection in the outer membrane protein C of Enterobacter aerogenes strains and Salmonella species

The types of selective pressure operating on the outer membrane protein C (ompC) of Enterobacter aerogenes strains, the causative agent for nosocomial infections, and Salmonella sp., the hazardous pathogen are investigated using the maximum likelihood-based codon substitution models. Although the rate of amino acid replacement to the silent substitution (omega) across the entire codon sites of ompC of E. aerogenes (omega=0.3194) and Salmonella sp. (omega=0.2047) indicate that the gene is subjected to purifying selection (i.e. omega<1), approximately 3.7% of ompC codon sites in E. aerogenes (omega=21.52) are under the influence of positive Darwinian selection (i.e. omega>1). Such contrast in the intensity of selective pressures in both pathogens could be associated with the differential response to the adverse environmental changes. In E. aerogenes, majority of the positively selected sites are located in the hypervariable cell-surface-exposed domains whereas the trans-membrane domains are functionally highly constrained.     

Accuracy and power of the likelihood ratio test in detecting adaptive molecular evolution

The selective pressure at the protein level is usually measured by the nonsynonymous/synonymous rate ratio (omega = dN/dS), with omega < 1, omega = 1, and omega > 1 indicating purifying (or negative) selection, neutral evolution, and diversifying (or positive) selection, respectively. The omega ratio is commonly calculated as an average over sites. As every functional protein has some amino acid sites under selective constraints, averaging rates across sites leads to low power to detect positive selection. Recently developed models of codon substitution allow the omega ratio to vary among sites and appear to be powerful in detecting positive selection in empirical data analysis. In this study, we used computer simulation to investigate the accuracy and power of the likelihood ratio test (LRT) in detecting positive selection at amino acid sites. The test compares two nested models: one that allows for sites under positive selection (with omega > 1), and another that does not, with the chi2 distribution used for significance testing. We found that use of the chi(2) distribution makes the test conservative, especially when the data contain very short and highly similar sequences. Nevertheless, the LRT is powerful. Although the power can be low with only 5 or 6 sequences in the data, it was nearly 100% in data sets of 17 sequences. Sequence length, sequence divergence, and the strength of positive selection also were found to affect the power of the LRT. The exact distribution assumed for the omega ratio over sites was found not to affect the effectiveness of the LRT.     

Molecular evolution of the ent-kaurenoic acid oxidase gene in Oryzeae

We surveyed the substitution patterns in the ent-kaurenoic acid oxidase (KAO) gene in 11 species of Oryzeae with an outgroup in the Ehrhartoidaea. The synonymous and non-synonymous substitution rates showed a high positive correlation with each other, but were negatively correlated with codon usage bias and GC content at third codon positions. The substitution rate was heterogenous among lineages. Likelihood-ratio tests showed that the non-synonymous/synonymous rate ratio changed significantly among lineages. Site-specific models provided no evidence for positive selection of particular amino acid sites in any codon of the KAO gene. This finding suggested that the significant rate heterogeneity among some lineages may have been caused by variability in the relaxation of the selective constraint among lineages or by neutral processes.     

Accuracy and power of bayes prediction of amino acid sites under positive selection

Bayes prediction quantifies uncertainty by assigning posterior probabilities. It was used to identify amino acids in a protein under recurrent diversifying selection indicated by higher nonsynonymous (d(N)) than synonymous (d(S)) substitution rates or by omega = d(N)/d(S) > 1. Parameters were estimated by maximum likelihood under a codon substitution model that assumed several classes of sites with different omega ratios. The Bayes theorem was used to calculate the posterior probabilities of each site falling into these site classes. Here, we evaluate the performance of Bayes prediction of amino acids under positive selection by computer simulation. We measured the accuracy by the proportion of predicted sites that were truly under selection and the power by the proportion of true positively selected sites that were predicted by the method. The accuracy was slightly better for longer sequences, whereas the power was largely unaffected by the increase in sequence length. Both accuracy and power were higher for medium or highly diverged sequences than for similar sequences. We found that accuracy and power were unacceptably low when data contained only a few highly similar sequences. However, sampling a large number of lineages improved the performance substantially. Even for very similar sequences, accuracy and power can be high if over 100 taxa are used in the analysis. We make the following recommendations: (1) prediction of positive selection sites is not feasible for a few closely related sequences; (2) using a large number of lineages is the best way to improve the accuracy and power of the prediction; and (3) multiple models of heterogeneous selective pressures among sites should be applied in real data analysis.     

Codon-substitution models to detect adaptive evolution that account for heterogeneous selective pressures among site classes.

The nonsynonymous to synonymous substitution rate ratio (omega = d(N)/d(S)) provides a sensitive measure of selective pressure at the protein level, with omega values <1, =1, and >1 indicating purifying selection, neutral evolution, and diversifying selection, respectively. Maximum likelihood models of codon substitution developed recently account for variable selective pressures among amino acid sites by employing a statistical distribution for the omega ratio among sites. Those models, called random-sites models, are suitable when we do not know a priori which sites are under what kind of selective pressure. Sometimes prior information (such as the tertiary structure of the protein) might be available to partition sites in the protein into different classes, which are expected to be under different selective pressures. It is then sensible to use such information in the model. In this paper, we implement maximum likelihood models for prepartitioned data sets, which account for the heterogeneity among site partitions by using different omega parameters for the partitions. The models, referred to as fixed-sites models, are also useful for combined analysis of multiple genes from the same set of species. We apply the models to data sets of the major histocompatibility complex (MHC) class I alleles from human populations and of the abalone sperm lysin genes. Structural information is used to partition sites in MHC into two classes: those in the antigen recognition site (ARS) and those outside. Positive selection is detected in the ARS by the fixed-sites models. Similarly, sites in lysin are classified into the buried and solvent-exposed classes according to the tertiary structure, and positive selection was detected at the solvent-exposed sites. The random-sites models identified a number of sites under positive selection in each data set, confirming and elaborating the results of the fixed-sites models. The analysis demonstrates the utility of the fixed-sites models, as well as the power of previous random-sites models, which do not use the prior information to partition sites.     

Human papillomavirus type 16 genetic variants: phylogeny and classification based on E6 and LCR

Naturally occurring genetic variants of human papillomavirus type 16 (HPV16) are common and have previously been classified into 4 major lineages; European-Asian (EAS), including the sublineages European (EUR) and Asian (As), African 1 (AFR1), African 2 (AFR2), and North-American/Asian-American (NA/AA). We aimed to improve the classification of HPV16 variant lineages by using a large resource of HPV16-positive cervical samples collected from geographically diverse populations in studies on HPV and/or cervical cancer undertaken by the International Agency for Research on Cancer. In total, we sequenced the entire E6 genes and long control regions (LCRs) of 953 HPV16 isolates from 27 different countries worldwide. Phylogenetic analyses confirmed previously described variant lineages and subclassifications. We characterized two new sublineages within each of the lineages AFR1 and AFR2 that are robustly classified using E6 and/or the LCR. We could differentiate previously identified AA1, AA2, and NA sublineages, although they could not be distinguished by E6 alone, requiring the LCR for correct phylogenetic classification. We thus provide a classification system for HPV16 genomes based on 13 and 32 phylogenetically distinguishing positions in E6 and the LCR, respectively, that distinguish nine HPV16 variant sublineages (EUR, As, AFR1a, AFR1b, AFR2a, AFR2b, NA, AA1, and AA2). Ninety-seven percent of all 953 samples fitted this classification perfectly. Other positions were frequently polymorphic within one or more lineages but did not define phylogenetic subgroups. Such a standardized classification of HPV16 variants is important for future epidemiological and biological studies of the carcinogenic potential of HPV16 variant lineages.     

Positive Darwinian selection at the pantophysin (Pan I) locus in marine gadid fishes

Maximum-likelihood models of codon substitution were used to test for positive Darwinian selection at the vesicle protein pantophysin in two allelic lineages segregating in the Atlantic cod Gadus morhua and in 18 related species of marine gadid fishes. Positive selection was detected in the two intravesicular loops of the integral membrane protein but not in four membrane-spanning regions or the 3' cytoplasmic tail. The proportion of positively selected sites (24.9%) and the mean nonsynonymous/synonymous rate ratio (omega = d(N)/d(S) = 5.35) were both greater in the first intravesicular (IV1) domain compared with the second intravesicular (IV2) domain (11.0% positively selected sites with mean omega = 3.76). Likelihood ratio tests comparing models that assume identical omega ratios along all branches of the phylogeny to those that allow omega ratios to vary among lineages were not significant for either the IV1 or IV2 domains, indicating that the selective pressures favoring amino acid replacements have operated consistently in both regions during the diversification of the group. Positive selection was observed in the IV1 domain in both G. morhua allelic lineages, and, although three of the four codons that differ between alleles were targets of positive selection in the broader group, no similar polymorphisms were detected in other taxa. The two G. morhua Pan I alleles appeared to have evolved before the speciation event separating it from its sister taxon, Theragra chalcogramma, and on the basis of a standard mtDNA clock are estimated to be at least 2 Myr old. Although the function of pantophysin remains unknown, the strong signal of positive selection at specific sites in the IV1 and IV2 domains may help clarify its role in cellular trafficking pathways.     

Likelihood Analysis of the Chalcone Synthase Genes Suggests the Role of Positive Selection in Morning Glories (<Emphasis Type="Italic">Ipomoea</Emphasis>)

Chalcone synthase (CHS) is a key enzyme in the biosynthesis of flavonoides, which are important for the pigmentation of flowers and act as attractants to pollinators. Genes encoding CHS constitute a multigene family in which the copy number varies among plant species and functional divergence appears to have occurred repeatedly. In morning glories (Ipomoea), five functional CHS genes (A–E) have been described. Phylogenetic analysis of the Ipomoea CHS gene family revealed that CHS A, B, and C experienced accelerated rates of amino acid substitution relative to CHS D and E. To examine whether the CHS genes of the morning glories underwent adaptive evolution, maximum-likelihood models of codon substitution were used to analyze the functional sequences in the Ipomoea CHS gene family. These models used the nonsynonymous/synonymous rate ratio ( = d_N/d_S) as an indicator of selective pressure and allowed the ratio to vary among lineages or sites. Likelihood ratio test suggested significant variation in selection pressure among amino acid sites, with a small proportion of them detected to be under positive selection along the branches ancestral to CHS A, B, and C. Positive Darwinian selection appears to have promoted the divergence of subfamily ABC and subfamily DE and is at least partially responsible for a rate increase following gene duplication.     

Evolutionary Dynamics of Variant Genomes of Human Papillomavirus Types 18, 45, and 97

Human papillomavirus type 18 (HPV18) and HPV45 account for approximately 20% of all cervix cancers. We show that HPV18, HPV45, and the recently discovered HPV97 comprise a clade sharing a most recent common ancestor within HPV α7 species. Variant lineages of these HPV types were classified by sequence analysis of the upstream regulatory region/E6 region among cervical samples from a population-based study in Costa Rica, and 27 representative genomes from each major variant lineage were sequenced. Nucleotide variation within HPV18 and HPV45 was 3.82% and 2.39%, respectively, and amino acid variation was 4.73% and 2.87%, respectively. Only 18 nucleotide variations, of which 10 were nonsynonymous, were identified among three HPV97 genomes. Full-genome comparisons revealed maximal diversity between HPV18 African and non-African variants (2.6% dissimilarity), whereas HPV18 Asian-American [E1 (AA)] and European (E2) variants were closely related (less than 0.5% dissimilarity); HPV45 genomes had a maximal difference of 1.6% nucleotides. Using a Bayesian Markov chain Monte Carlo (MCMC) method, the divergence times of HPV18, -45, and -97 from their most recent common ancestors indicated that HPV18 diverged approximately 7.7 million years (Myr) ago, whereas HPV45 and HPV97 split off around 5.7 Myr ago, in a period encompassing the divergence of the great ape species. Variants within the HPV18/45/97 lineages were estimated to have diverged from their common ancestors in the genus Homo within the last 1 Myr (<0.7 Myr). To investigate the molecular basis of HPV18, HPV45, and HPV97 evolution, regression models of codon substitution were used to identify lineages and amino acid sites under selective pressure. The E5 open reading frame (ORF) of HPV18 and the E4 ORFs of HPV18, HPV45, and HPV18/45/97 had nonsynonymous/synonymous substitution rate ratios (d_N/d_S) over 1 indicative of positive Darwinian selection. The L1 ORF of HPV18 genomes had an increased proportion of nonsynonymous substitutions (4.93%; average d_N/d_S ratio [M3] = 0.3356) compared to HPV45 (1.86%; M3 = 0.1268) and HPV16 (2.26%; M3 = 0.1330) L1 ORFs. In contrast, HPV18 and HPV16 genomes had similar amino acid substitution rates within the E1 ORF (2.89% and 3.24%, respectively), while HPV45 E1 was highly conserved (amino acid substitution rate was 0.77%). These data provide an evolutionary history of this medically important clade of HPVs and identify an unexpected divergence of the L1 gene of HPV18 that may have clinical implications for the long-term use of an L1-virus-like particle-based prophylactic vaccine.Papillomaviruses (PVs) are a large family of related viruses with circular double-stranded DNA genomes 8 kb in size. Some PV types cause epithelial hyperplasias ranging from benign exophytic warts to premalignant lesions that can progress to invasive cancer. Among the 61 currently designated alpha human PVs (HPVs), the majority have been isolated from the mucosal surface of the genital or oral region (8, 14). Of these, a select group have oncogenic potential and are associated with cervical cancer (11). Specifically, HPV type 16 (HPV16) and HPV18 have been identified in approximately two-thirds of cervical cancers, this tumor is the second most common cancer in women, and it is the principal cancer of women in developing countries (5, 24, 25, 30, 37).To date, studies of HPV18 variants have identified three lineages corresponding to the continental locations where the viral samples were obtained: European (E), Asian-American (AA), and African (Af) (29). The phylogeny of HPV18 variants is reflective of the migration patterns of Homo sapiens and suggests that HPV18 variant lineages might have diverged through genetic isolation at approximately the same time as Homo sapiens began establishing residence in different continental regions. Previous HPV18 intratypic phylogenetic analyses were limited to partial regions of the genome (3, 7, 29). Nevertheless, studies also suggest that HPV18 variants are associated with different levels of oncogenic potential and persistence and histological tumor types (1, 6, 35, 36, 46).HPV45 and HPV97 are the viral types most closely related to HPV18 and taken together form a clade and share a most recent common ancestor (MRCA). HPV97 is a recently described rare type (8, 17). HPV18 and HPV45 account for approximately 20% of all cervix cancers (25). Although HPV45 is a common type found in cervical cancer, its evolutionary history and sequence variability have not been extensively studied.In this report, 27 complete genomes representing the major variant lineages of HPV18, HPV45, and HPV97 were cloned and/or sequenced from clinical samples. Based on full genomes, the intratype/intertype evolutionary trees of HPV18, HPV45, and HPV97 were constructed. By examining the rate ratio of nonsynonymous (d_N) to synonymous (d_S) substitutions per site, diversifying selection acting on each of the eight protein-encoding regions of HPV18, HPV45, and HPV97 was evaluated. In addition, the times of divergence of HPV18/45/97 variants from their MRCA were investigated. These data provide an evolutionary history of this medically important clade of HPVs.     

Molecular Adaptation of rbcL in the Heterophyllous Aquatic Plant Potamogeton

Background

Heterophyllous aquatic plants show marked phenotypic plasticity. They adapt to environmental changes by producing different leaf types: submerged, floating and terrestrial leaves. By contrast, homophyllous plants produce only submerged leaves and grow entirely underwater. Heterophylly and submerged homophylly evolved under selective pressure modifying the species-specific optima for photosynthesis, but little is known about the evolutionary outcome of habit. Recent evolutionary analyses suggested that rbcL, a chloroplast gene that encodes a catalytic subunit of RuBisCO, evolves under positive selection in most land plant lineages. To examine the adaptive evolutionary process linked to heterophylly or homophylly, we analyzed positive selection in the rbcL sequences of ecologically diverse aquatic plants, Japanese Potamogeton.

Principal Findings

Phylogenetic and maximum likelihood analyses of codon substitution models indicated that Potamogeton rbcL has evolved under positive Darwinian selection. The positive selection has operated specifically in heterophyllous lineages but not in homophyllous ones in the branch-site models. This suggests that the selective pressure on this chloroplast gene was higher for heterophyllous lineages than for homophyllous lineages. The replacement of 12 amino acids occurred at structurally important sites in the quaternary structure of RbcL, two of which (residue 225 and 281) were identified as potentially under positive selection.

Conclusions/Significance

Our analysis did not show an exact relationship between the amino acid replacements and heterophylly or homophylly but revealed that lineage-specific positive selection acted on the Potamogeton rbcL. The contrasting ecological conditions between heterophyllous and homophyllous plants have imposed different selective pressures on the photosynthetic system. The increased amino acid replacement in RbcL may reflect the continuous fine-tuning of RuBisCO under varying ecological conditions.     

Differential selection of genes of cucumber mosaic virus subgroups

Cucumber mosaic virus (CMV) has an extremely broad plant-host range, a large number of vector species, and a wide geographical distribution. CMV is, therefore, a model by which to understand plant virus adaptation. The selective constraints exerted on the five proteins expressed from the CMV genome were evaluated by application of newly developed maximum-likelihood algorithms to analyze sequences available in data banks. The ratio between nonsynonymous and synonymous substitution rates (omega) was used to detect positive selection on particular codon sites. Amino acid sequences were conserved with omega ranging from 0.07 to 0.60 in different proteins. However, a small proportion of amino acids in proteins 1a, 2a, and 3b, the coat protein (CP), were positively selected (omega > 1). Moreover, the evolution of the CP in the three subgroups of CMV strains revealed different selection profiles along the sequence and significantly different speed of evolution at many positions. Constraints exerted by aphid transmission, rather than plant adaptation, seemed to be responsible for these patterns of evolution in the CP.     

Pattern of nucleotide substitution and divergence of prophenoloxidase in decapods

Despite the unprecedented development in identification and characterization of prophenoloxidase (proPO) in commercially important decapods, little is known about the evolutionary relationship, rate of amino acid replacement and differential selection pressures operating on proPO of different species of decapods. Here we report the evolutionary relationship among these nine decapod species based on proPO gene and types of selective pressures operating on proPO codon sites. Our analyses revealed that all the nine decapod species shared a common ancestor. The mean percentage sequence divergence at proPO gene was 34.4+/-0.6%. Pairwise estimates of nonsynonymous to synonymous ratio (omega) for Homarus americanus-H. gammarus is greater than one, therefore indicating adaptive evolution (functional diversification) of proPO in these two species. In contrast, strong purifying selection (omega<1) was observed in all other species pairs. However, phylogenetically closely related decapods revealed relatively higher omega value (omega=0.15+/-0.3) than the distantly related species pairs (omega=0.0075+/-0.005). These discrepancies could be due to higher fixation probability of beneficial mutation in closely related species. Maximum likelihood-based codon substitution analyses revealed a strong purifying selection operating on most of the codon sites, therefore suggesting proPO is functionally constrained (purifying selection). Codon substitution analyses have also revealed the evidence of strong purifying selection in haemocyanin subunits of decapods.     

Development of a DNA vaccine targeting human papillomavirus type 16 oncoprotein E6

Human papillomavirus (HPV), particularly type 16 (HPV-16), is present in more than 99% of cervical cancers. The HPV oncoproteins E6 and E7 are constantly expressed and therefore represent ideal targets for HPV vaccine development. We previously developed DNA vaccines encoding calreticulin (CRT) linked to HPV-16 E7 and generated potent E7-specific CD8(+) T-cell immune responses and antitumor effects against an E7-expressing tumor. Since vaccines targeting E6 also represent an important strategy for controlling HPV-associated lesions, we developed a DNA vaccine encoding CRT linked to E6 (CRT/E6). Our results indicated that the CRT/E6 DNA vaccine, but not a wild-type E6 DNA vaccine, generated significant E6-specific CD8(+) T-cell immune responses in vaccinated mice. Mapping of the immunodominant epitope of E6 revealed that an E6 peptide comprising amino acids (aa) 48 to 57 (E6 aa48-57), presented by H-2K(b), is the optimal peptide and that the region of E6 comprising aa 50 to 57 represents the minimal core sequence required for activating E6-specific CD8(+) T lymphocytes. We also demonstrated that E6 aa48-57 contains cytotoxic T-lymphocyte epitopes naturally presented by E6-expressing TC-1 cells. Vaccination with a CRT/E6 but not a CRT/mtE6 (lacking aa 50 to 57 of E6) DNA vaccine could protect vaccinated mice from challenge with E6-expressing TC-1 tumors. Thus, our data indicate that E6 aa48-57 contains the immunodominant epitope and that a CRT/E6 DNA vaccine may be useful for control of HPV infection and HPV-associated lesions.     

Characterization of human antibody-reactive epitopes encoded by human papillomavirus types 16 and 18.

We have previously reported that the most common human serum immunoglobulin G antibody reactivities to human papillomavirus type 16 and type 18 (HPV16 and HPV18)-encoded proteins are directed against the minor capsid proteins (HPV16 L2 and HPV18 L2) and to the E7 protein of HPV16 (S. A. Jenison, X.-P. Yu, J. M. Valentine, L. A. Koutsky, A. E. Christiansen, A. M. Beckmann, and D. A. Galloway, J. Infect. Dis. 162:60-69, 1990). In this study, the antibody-reactive segments of the HPV16 E7, HPV16 L2, and HPV18 L2 polypeptides were mapped by using nested sets of deleted recombinant proteins. A single major immunoreactive region was identified in the HPV16 E7 polypeptide between amino acids (aa) 21 and 34 (DLYCYE-QLNDSSEE). In contrast, three distinct immunoreactive regions of the HPV16 L2 polypeptide were present in the segment between aa149 and aa204, and three distinct immunoreactive regions of the HPV18 L2 polypeptide were present in the segment between aa110 and aa211. With the exception of one serum sample, serum immunoglobulin G antibodies which reacted with HPV16 L2 polypeptides or with HPV18 L2 polypeptides were not cross-reactive.     

The rate of adaptive evolution in enteric bacteria

Here we estimate the rate of adaptive substitution in a set of 410 genes that are present in 6 Escherichia coli and 6 Salmonella enterica genomes. We estimate that more than 50% of amino acid substitutions in this set of genes have been fixed by positive selection between the E. coli and S. enterica lineages. We also show that the proportion of adaptive substitutions is uncorrelated with the rate of amino acid substitution or gene function but that it may be correlated with levels of synonymous codon usage bias.     

Mapping and characterization of the interaction domains of human papillomavirus type 16 E1 and E2 proteins.

The papillomavirus E1 and E2 proteins are both necessary and sufficient in vivo for efficient origin-dependent viral DNA replication. The ability of E1 and E2 to complex with each other appears to be essential for efficient viral DNA replication. In this study, we used the yeast two-hybrid system and in vitro binding assays to map the domains of the human papillomavirus type 16 (HPV16) E1 and E2 proteins required for complex formation. The amino-terminal 190-amino-acid domain of HPV16 E2 was both required and sufficient for E1 binding. The carboxyl-terminal 229 amino acids of E 1 were essential for binding E2, and the amino-terminal 143 amino acids of HPV16 E1 were dispensable. Although the ability of the E1 minimal domain (amino acids [aa] 421 to 649) to interact with E2 was strong at 4 degrees C, it was significantly reduced at temperatures above 25 degrees C. A larger domain of E1 from aa 144 to 649 bound E2 efficiently at any temperature, suggesting that aa 144 to 420 of E1 may play a role in the HPV16 E1-E2 interaction at physiological temperatures.     

Adaptive Evolution in the GAF Domain of Phytochromes in Gymnosperms

The GAF domain of phytochrome is essential for photoconversion and signal transduction. In gymnosperms, it exists in all members of the phytochrome family that experience gene duplication. Maximum-likelihood models of codon substitution can provide a framework for constructing likelihood ratio tests of changes in selective pressure and make clear predictions about patterns of genetic change following gene duplication. In this study, 68 gymnosperm GAF sequences were analyzed to identify lineages and sites under positive selection. Our results indicate that (1) positive selection at a few sites (3.6%), rather than relaxation of selective constraints, has played a major role in the evolution of the gymnosperm GAF domain; (2) strong positive selective pressure tends to occur in the recent PHYP lineages of cogeneric species, but is absent in old lineages consisting of distantly related species; and (3) the selective pressure indicated by the ω ratio varies greatly among lineages and sites in the GAF domain.     

Codon-substitution models for detecting molecular adaptation at individual sites along specific lineages

The nonsynonymous (amino acid-altering) to synonymous (silent) substitution rate ratio (omega = d(N)/d(S)) provides a measure of natural selection at the protein level, with omega = 1, >1, and <1, indicating neutral evolution, purifying selection, and positive selection, respectively. Previous studies that used this measure to detect positive selection have often taken an approach of pairwise comparison, estimating substitution rates by averaging over all sites in the protein. As most amino acids in a functional protein are under structural and functional constraints and adaptive evolution probably affects only a few sites at a few time points, this approach of averaging rates over sites and over time has little power. Previously, we developed codon-based substitution models that allow the omega ratio to vary either among lineages or among sites. In this paper we extend previous models to allow the omega ratio to vary both among sites and among lineages and implement the new models in the likelihood framework. These models may be useful for identifying positive selection along prespecified lineages that affects only a few sites in the protein. We apply those branch-site models as well as previous branch- and site-specific models to three data sets: the lysozyme genes from primates, the tumor suppressor BRCA1 genes from primates, and the phytochrome (PHY) gene family in angiosperms. Positive selection is detected in the lysozyme and BRCA genes by both the new and the old models. However, only the new models detected positive selection acting on lineages after gene duplication in the PHY gene family. Additional tests on several data sets suggest that the new models may be useful in detecting positive selection after gene duplication in gene family evolution.     

Comparison of the complete sequence of the str operon in Salmonella typhimurium and Escherichia coli.

The nucleotide (nt) sequences of the str operon in Escherichia coli K-12 and Salmonella typhimurium LT2 were completed and compared at the nt and amino acid (aa) level. The order of conservation at the nt and aa level is rpsL greater than tufA greater than rpsG greater than f usA. A striking difference is that the rpsG-encoded ribosomal protein, S7, in E. coli K-12 is 23 aa longer than in S. typhimurium. The very low (0.18) codon adaptation index of this part of the E. coli K-12-encoding gene and the unusual stop codon (UGA) suggest that this is a relatively recent extension. A trend towards a higher G+C content in fusA (gene encoding elongation factor (EF)-G) and tufA (gene encoding EF-Tu) in S. typhimurium is noted. In fusA, nt substitutions at all three positions in a codon occur at a much higher frequency than expected from the number of nt substitutions in the gene, assuming they are random and independent events. An analysis of substitutions in this and other genes suggests that the triple substitutions in fusA, and some other genes, are the result of the sequential accumulation of individual mutations, probably driven by selection pressure for particular codons or aa.