PubMed Health⌕ Search

SEARCH · PubMed Health

Results for “Ka/Ks”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 recordsLinked to original sources
PubMed2005

Assessing the application of Ka/Ks ratio test to alternatively spliced exons.

SUMMARY: Recently, the Ka/Ks ratio test, which assesses the protein-coding potentials of genomic regions based on their non-synonymous to synonymous divergence rates, has been proposed and successfully used in genome annotations of eukaryotes. We systematically performed the Ka/Ks ratio test on 925 transcript-confirmed alternatively spliced exons in the human genome, which we describe in this manuscript. We found that 22.3% of evolutionarily conserved alternatively spliced exons cannot pass the Ka/Ks ratio test, compared with 9.8% for constitutive exons. The false negative rate was the highest (85.7%) for exons with low frequencies of transcript inclusion. Analyses of alternatively spliced exons supported by full-length mRNA sequences yielded similar results, and nearly half of exons involved in ancestral alternative splicing events could not pass this test. Our analysis suggests a future direction to incorporate comparative genomics-based alternative splicing predictions with the Ka/Ks ratio test in higher eukaryotes with extensive RNA alternative splicing.

Algorithms↗
PubMed2004

Average allozyme heterozygosity in vertebrates correlates with Ka/Ks measured in the human-mouse lineage.

It is well established that different allozyme proteins vary in heterozygosity in averages made over large numbers of species. For example, the enzyme 6-phosphogluconate dehydrogenase has a much higher average heterozygosity than glutamate dehydrogenase. Allozyme data alone provide insufficient power to determine the evolutionary cause of such a difference. Many studies have now been carried out on the DNA sequences coding for allozymes. These have identified diverse selective and nonselective causes of polymorphisms at individual loci. However the studies are mainly in a small number of model species; thus, it is difficult to identify from these DNA studies specific causes of global average heterozygosity differences among allozyme proteins. Here we demonstrate that estimates of average heterozygosity for 37 allozyme proteins in vertebrates correlate positively with Ka and Ka/Ks but not with Ks, measured in the human-mouse lineage. The values of Ka/Ks are less than 0.25, and Ka/Ks is negatively correlated with subunit number (quaternary structure), a measure of structural constraint. Proteins with lower levels of constraint have higher values of both Ka/Ks and heterozygosity. These results better support the hypothesis that differences in average allozyme diversity between proteins are more closely related to differences in the level of purifying selection than to differences in the underlying mutation rate or level of positive selection.

Animals↗
PubMed2006

Can RNA selection pressure distort the measurement of Ka/Ks?

Recently, an interesting question has emerged in the evolutionary interpretation of sequence substitution data as evidence of amino acid selection pressure. Specifically, the Ka/Ks metric was designed to measure selection pressure on amino acid substitutions, assuming that the synonymous substitution rate Ks reflects the neutral nucleotide substitution rate. However, there is increasing evidence for selection pressure at silent sites due to constraints of RNA splicing. Is Ka/Ks an appropriate metric for selection pressure on amino acid substitutions, in the presence of other selection pressures acting only at the RNA level (such as selection for exonic splicing enhancers)? Or can the resulting decreases in Ks from such selection pressures introduce bias into the Ka/Ks metric, so that it no longer gives an accurate measure of amino acid level selection pressure? In this review, we present both mathematical models and empirical evidence for these divergent points of view.

Amino Acid Substitution↗
PubMed2006

Distinguishing HIV-1 drug resistance, accessory, and viral fitness mutations using conditional selection pressure analysis of treated versus untreated patient samples.

BACKGROUND: HIV can evolve drug resistance rapidly in response to new drug treatments, often through a combination of multiple mutations 123. It would be useful to develop automated analyses of HIV sequence polymorphism that are able to predict drug resistance mutations, and to distinguish different types of functional roles among such mutations, for example, those that directly cause drug resistance, versus those that play an accessory role. Detecting functional interactions between mutations is essential for this classification. We have adapted a well-known measure of evolutionary selection pressure (Ka/Ks) and developed a conditional Ka/Ks approach to detect important interactions. RESULTS: We have applied this analysis to four independent HIV protease sequencing datasets: 50,000 clinical samples sequenced by Specialty Laboratories, Inc.; 1800 samples from patients treated with protease inhibitors; 2600 samples from untreated patients; 400 samples from untreated African patients. We have identified 428 mutation interactions in Specialty dataset with statistical significance and we were able to distinguish primary vs. accessory mutations for many well-studied examples. Amino acid interactions identified by conditional Ka/Ks matched 80 of 92 pair wise interactions found by a completely independent study of HIV protease (p-value for this match is significant: 10-70). Furthermore, Ka/Ks selection pressure results were highly reproducible among these independent datasets, both qualitatively and quantitatively, suggesting that they are detecting real drug-resistance and viral fitness mutations in the wild HIV-1 population. CONCLUSION: Conditional Ka/Ks analysis can detect mutation interactions and distinguish primary vs. accessory mutations in HIV-1. Ka/Ks analysis of treated vs. untreated patient data can distinguish drug-resistance vs. viral fitness mutations. Verification of these results would require longitudinal studies. The result provides a valuable resource for AIDS research and will be available for open access upon publication at http://www.bioinformatics.ucla.edu/HIV.

Journal Article↗
PubMed2006

The effects of multiple features of alternatively spliced exons on the K(A)/K(S) ratio test.

BACKGROUND: The evolution of alternatively spliced exons (ASEs) is of primary interest because these exons are suggested to be a major source of functional diversity of proteins. Many exon features have been suggested to affect the evolution of ASEs. However, previous studies have relied on the KA/KS ratio test without taking into consideration information sufficiency (i.e., exon length > 75 bp, cross-species divergence > 5%) of the studied exons, leading to potentially biased interpretations. Furthermore, which exon feature dominates the results of the KA/KS ratio test and whether multiple exon features have additive effects have remained unexplored. RESULTS: In this study, we collect two different datasets for analysis - the ASE dataset (which includes lineage-specific ASEs and conserved ASEs) and the ACE dataset (which includes only conserved ASEs). We first show that information sufficiency can significantly affect the interpretation of relationship between exons features and the KA/KS ratio test results. After discarding exons with insufficient information, we use a Boolean method to analyze the relationship between test results and four exon features (namely length, protein domain overlapping, inclusion level, and exonic splicing enhancer (ESE) frequency) for the ASE dataset. We demonstrate that length and protein domain overlapping are dominant factors, and they have similar impacts on test results of ASEs. In addition, despite the weak impacts of inclusion level and ESE motif frequency when considered individually, combination of these two factors still have minor additive effects on test results. However, the ACE dataset shows a slightly different result in that inclusion level has a marginally significant effect on test results. Lineage-specific ASEs may have contributed to the difference. Overall, in both ASEs and ACEs, protein domain overlapping is the most dominant exon feature while ESE frequency is the weakest one in affecting test results. CONCLUSION: The proposed method can easily find additive effects of individual or multiple factors on the KA/KS ratio test results of exons. Therefore, the system can analyze complex conditions in evolution where multiple features are involved. More factors can also be added into the system to extend the scope of evolutionary analysis of exons. In addition, our method may be useful when orthologous exons can not be found for the KA/KS ratio test.

Alternative Splicing↗
PubMed1998

Molecular evolution of an imprinted gene: repeatability of patterns of evolution within the mammalian insulin-like growth factor type II receptor.

The repeatability of patterns of variation in Ka/Ks and Ks is expected if such patterns are the result of deterministic forces. We have contrasted the molecular evolution of the mammalian insulin-like growth factor type II receptor (Igf2r) in the mouse-rat comparison with that in the human-cow comparison. In so doing, we investigate explanations for both the evolution of genomic imprinting and for Ks variation (and hence putatively for mutation rate evolution). Previous analysis of Igf2r, in the mouse-rat comparison, found Ka/Ks patterns that were suggested to be contrary to those expected under the conflict theory of imprinting. We find that Ka/Ks variation is repeatable and hence confirm these patterns. However, we also find that the molecular evolution of Igf2r signal sequences suggests that positive selection, and hence conflict, may be affecting this region. The variation in Ks across Igf2r is also repeatable. To the best of our knowledge this is the first demonstration of such repeatability. We consider three explanations for the variation in Ks across the gene: (1) that it is the result of mutational biases, (2) that it is the result of selection on the mutation rate, and (3) that it is the product of selection on codon usage. Explanations 2 and 3 predict a Ka-Ks correlation, which is not found. Explanation 3 also predicts a negative correlation between codon bias and Ks, which is also not found. However, in support of explanation 1 we do find that in rodents the rate of silent C --> T mutations at CpG sites does covary with Ks, suggesting that methylation-induced mutational patterns can explain some of the variation in Ks. We find evidence to suggest that this CpG effect is due to both variation in CpG density, and to variation in the frequency with which CpGs mutate. Interestingly, however, a GC4 analysis shows no covariance with Ks, suggesting that to eliminate methyl-associated effects CpG rates themselves must be analyzed. These results suggest that, in contrast to previous studies of intragenic variation, Ks patterns are not simply caused by the same forces responsible for Ka/Ks correlations.

Animals↗
PubMed2003

Asymmetric sequence divergence of duplicate genes.

Much like humans, gene duplicates may be created equal, but they do not stay that way for long. For four completely sequenced genomes we show that 20%-30% of duplicate gene pairs show asymmetric evolution in the amino acid sequence of their protein products. That is, one of the duplicates evolves much faster than the other. The greater this asymmetry, the greater the ratio Ka/Ks of amino acid substitutions (Ka) to silent substitutions (Ks) in a gene pair. This indicates that most asymmetric divergence may be caused by relaxed selective constraints on one of the duplicates. However, we also find some candidate duplicates where positive (directional) selection of beneficial mutations (Ka/Ks > 1) may play a role in asymmetric divergence. Our analysis rests on a codon-based model of molecular evolution that allows a test for asymmetric divergence in Ka. The method is also more sensitive in detecting positive selection (Ka/Ks > 1) than models relying only on pairwise gene comparisons.

Amino Acid Substitution↗
PubMed2003

Adaptive evolution of MRG, a neuron-specific gene family implicated in nociception.

The MRG gene family (also known as SNSR) belongs to the G-protein-coupled receptor (GPCR) superfamily, is expressed specifically in nociceptive neurons, and is implicated in the modulation of nociception. Here, we show that Ka/Ks (the ratio between nonsynonymous and synonymous substitution rates) displays distinct profiles along the coding regions of MRG, with peaks (Ka/Ks>1) corresponding to extracellular domains, and valleys (Ka/Ks<1) corresponding to transmembrane and cytoplasmic domains. The extracellular domains are also characterized by a significant excess of radical amino acid changes. Statistical analysis shows that positive selection is by far the most suitable model to account for the nucleotide substitution patterns in MRG. Together, these results demonstrate that the extracellular domains of the MRG receptor family, which presumably partake in ligand binding, have experienced strong positive selection. Such selection is likely directed at altering the sensitivity and/or selectivity of nociceptive neurons to aversive stimuli. Thus, our finding suggests pain perception as an aspect of the nervous system that may have experienced a surprising level of adaptive evolution.

Adaptation, Physiological↗
PubMed2005

Analyses of coding region polymorphisms in apical and basolateral human organic anion transporter (OAT) genes [OAT1 (NKT), OAT2, OAT3, OAT4, URAT (RST)].

BACKGROUND: Excretion by the kidney of a variety of organic anionic drugs and metabolites is mediated by a family of multispecific organic anion transporters (OAT genes) that are part of the SLC22 family of solute carriers. Different OATs localize to the apical (OAT2, OAT4, and RST/URAT) or basolateral (OAT1/NKT and OAT3) membranes of the renal proximal tubule; the net transport of organic anions from blood to urine is believed to require both apical and basolateral OATs. These genes are also thought to mediate transport of organic anionic drugs and metabolites (e.g., urate) across choroid plexus, retina, placenta, and possibly olfactory mucosa. The extent of functional redundancy among OATs remains uncertain, but closely related OAT genes are tightly linked in the genome. Hence, a better understanding of human variation in organic anionic drug excretion may be obtained by studying OAT genes in combination rather than individually. METHODS: We have analyzed single nucleotide polymorphisms (SNPs) in OAT1 (NKT), OAT2, OAT3, OAT4, and URAT1 (human homologue of RST) in an ethnically diverse sample of 96 individuals (192 haploid genomes). Ka/Ks analysis was also performed as well as haplotype reconstruction using the software program Arelquin. RESULTS: The data indicate that (1) nonsynonymous SNPs in OAT1 and OAT3 may not be frequent so it will be important to consider promoter region SNPs that regulate gene expression; (2) certain ethnic groups may have a high prevalence of nonsynonymous SNPs in particular OATs (e.g., OAT4 in Sub-Saharan Africans); (3) there are individuals who have nonsynonymous SNPs in apical and basolateral OATs; (4) nonsynonymous OAT4 SNPs may be more frequent, raising the possibility of altered maternofetal transport of drugs and metabolites; and (5) combinations of synonymous SNPs in OAT1 and OAT3 also occur in certain individuals. In addition, Ka/Ks analysis of human, chimp and rodent genes suggests that OAT4 is under accelerated selection pressure, perhaps reflecting specific human environmental exposures during evolution. In contrast, Ka/Ks analysis for URAT1 suggests decelerated selection pressure. Haplotype reconstruction also supports this view. CONCLUSION: Together, these data suggest that, in order to understand the effect of SNPs in genes of the SLC22 family on drug handling as well as excretion of metabolites like uric acid, it is important to consider the entire set of organic anion transporters. It will be particularly interesting to determine if individuals with nonsynonymous apical and basolateral SNPs have altered handling (and toxicity) of organic anionic drugs and metabolites. Certain OAT family members appear to be under greater evolutionary selection pressure.

Carrier Proteins↗
PubMed2006

A systematic search for positive selection in higher plants (Embryophytes).

BACKGROUND: Previously, a database characterizing examples of Embryophyte gene family lineages showing evidence of positive selection was reported. Of the gene family trees, 138 Embryophyte branches showed Ka/Ks>>1 and are candidates for functional adaptation. The database and these examples have now been studied in further detail to better understand the molecular basis for plant genome evolution. RESULTS: Neutral modeling showed an excess of positive and/or negative selection in the database over a neutral expectation centered on the mean Ka/Ks ratio. Out of 673 families with assigned structures, 490 have at least one branch with Ka/Ks >>1 in a region of the protein, enabling a picture of selective pressures delineated by protein structure. Most gene families allowed reconstruction back to the last common ancestor of flowering plants (Magnoliophytes) without saturation of 4- fold degenerate codon position. Positive selection occurred in a wide variety of gene families with different functions, including in the self incompatibility locus, in defense against pathogens, in embryogenesis, in cold acclimation, and in electrontransport. Structurally, selective pressures were similar between alpha-helices and beta- sheets, but were less negative and more variant on the surface and away from the hydrophobic core. CONCLUSION: Positive selection was detected statistically significantly in a small and nonrandom minority of gene families in a systematic analysis of embryophyte gene families. More sensitive methods increased the level of positive selection that was detected and presented a structural basis for the role of positive selection in plant genomes.

Adaptation, Physiological↗
PubMed2001

The Adaptive Evolution Database (TAED).

BACKGROUND: Developing an understanding of the molecular basis for the divergence of species lies at the heart of biology. The Adaptive Evolution Database (TAED) serves as a starting point to link events that occur at the same time in the evolutionary history (tree of life) of species, based upon coding sequence evolution analyzed with the Master Catalog. The Master Catalog is a collection of evolutionary models, including multiple sequence alignments, phylogenetic trees, and reconstructed ancestral sequences, for all independently evolving protein sequence modules encoded by genes in GenBank [1]. RESULTS: We have estimated from these models the ratio of nonsynonymous to synonymous nucleotide substitution (Ka/Ks), for each branch in their respective evolutionary trees of every subtree containing only chordata or only embryophyta proteins. Branches with high Ka/Ks values represent candidate episodes in the history of the family where the protein may have undergone positive selection, a phenomenon in molecular evolution where the mutant form of a gene must have conferred more fitness than the ancestral form. Such episodes are frequently associated with change in function. We have found that an unexpectedly large number of families (between 10 and 20% of those families examined) have at least one branch with a notably high Ka/Ks value (putative adaptive evolution). As a resource for biologists wishing to understand the interaction between protein sequences and the Darwinian processes that shape these sequences, we have collected these into The Adaptive Evolution Database (TAED). CONCLUSIONS: Placed in a phylogenetic perspective, candidate genes that are undergoing evolution at the same time in the same lineage can be viewed together. This framework based upon coding sequence evolution can be readily expanded to include other types of evolution. In its present form, TAED provides a resource for bioinformaticists interested in data mining and for experimental evolutionists seeking candidate examples of adaptive evolution for further experimental study.

Animals↗
PubMed2001

The adaptive evolution database (TAED).

BACKGROUND: The Master Catalog is a collection of evolutionary families, including multiple sequence alignments, phylogenetic trees and reconstructed ancestral sequences, for all protein-sequence modules encoded by genes in GenBank. It can therefore support large-scale genomic surveys, of which we present here The Adaptive Evolution Database (TAED). In TAED, potential examples of positive adaptation are identified by high values for the normalized ratio of nonsynonymous to synonymous nucleotide substitution rates (KA/KS values) on branches of an evolutionary tree between nodes representing reconstructed ancestral sequences. RESULTS: Evolutionary trees and reconstructed ancestral sequences were extracted from the Master Catalog for every subtree containing proteins from the Chordata only or the Embryophyta only. Branches with high KA/KS values were identified. These represent candidate episodes in the history of the protein family when the protein may have undergone positive selection, where the mutant form conferred more fitness than the ancestral form. Such episodes are frequently associated with change in function. An unexpectedly large number of families (between 10% and 20% of those families examined) were found to have at least one branch with high KA/KS values above arbitrarily chosen cut-offs (1 and 0.6). Most of these survived a robustness test and were collected into TAED. CONCLUSIONS: TAED is a raw resource for bioinformaticists interested in data mining and for experimental evolutionists seeking candidate examples of adaptive evolution for further experimental study. It can be expanded to include other evolutionary information (for example changes in gene regulation or splicing) placed in a phylogenetic perspective.

Adaptation, Physiological↗
PubMed2006

Multiple mechanisms promote the retained expression of gene duplicates in the tetraploid frog Xenopus laevis.

Gene duplication provides a window of opportunity for biological variants to persist under the protection of a co-expressed copy with similar or redundant function. Duplication catalyzes innovation (neofunctionalization), subfunction degeneration (subfunctionalization), and genetic buffering (redundancy), and the genetic survival of each paralog is triggered by mechanisms that add, compromise, or do not alter protein function. We tested the applicability of three types of mechanisms for promoting the retained expression of duplicated genes in 290 expressed paralogs of the tetraploid clawed frog, Xenopus laevis. Tests were based on explicit expectations concerning the ka/ks ratio, and the number and location of nonsynonymous substitutions after duplication. Functional constraints on the majority of paralogs are not significantly different from a singleton ortholog. However, we recover strong support that some of them have an asymmetric rate of nonsynonymous substitution: 6% match predictions of the neofunctionalization hypothesis in that (1) each paralog accumulated nonsynonymous substitutions at a significantly different rate and (2) the one that evolves faster has a higher ka/ks ratio than the other paralog and than a singleton ortholog. Fewer paralogs (3%) exhibit a complementary pattern of substitution at the protein level that is predicted by enhancement or degradation of different functional domains, and the remaining 13% have a higher average ka/ks ratio in both paralogs that is consistent with altered functional constraints, diversifying selection, or activity-reducing mutations after duplication. We estimate that these paralogs have been retained since they originated by genome duplication between 21 and 41 million years ago. Multiple mechanisms operate to promote the retained expression of duplicates in the same genome, in genes in the same functional class, over the same period of time following duplication, and sometimes in the same pair of paralogs. None of these paralogs are superfluous; degradation or enhancement of different protein subfunctions and neofunctionalization are plausible hypotheses for the retained expression of some of them. Evolution of most X. laevis paralogs, however, is consistent with retained expression via mechanisms that do not radically alter functional constraints, such as selection to preserve post-duplication stoichiometry or temporal, quantitative, or spatial subfunctionalization.

Amino Acid Substitution↗
PubMed2005

Molecular population genetics of accessory gland protein genes and testis-expressed genes in Drosophila mojavensis and D. arizonae.

Molecular population genetic investigation of Drosophila male reproductive genes has focused primarily on melanogaster subgroup accessory gland protein genes (Acp's). Consistent with observations from male reproductive genes of numerous taxa, Acp's evolve more rapidly than nonreproductive genes. However, within the Drosophila genus, large data sets from additional types of male reproductive genes and from different species groups are lacking. Here we report findings from a molecular population genetics analysis of male reproductive genes of the repleta group species, Drosophila arizonae and D. mojavensis. We find that Acp's have dramatically higher average pairwise Ka/Ks (0.93) than testis-enriched genes (0.19) and previously reported melanogaster subgroup Acp's (0.42). Overall, 10 of 19 Acp's have Ka/Ks > 1 either in nonpolarized analyses or in at least one lineage of polarized analyses. Of the nine Acp's for which outgroup data were available, average Ka/Ks was considerably higher in D. mojavensis (2.08) than in D. arizonae (0.87). Contrasts of polymorphism and divergence suggest that adaptive protein evolution at Acp's is more common in D. mojavensis than in D. arizonae.

Animals↗
PubMed2000

[Different evolution rule of SRY gene between hominoid and Old World monkey].

The complete SRY gene of P. nemaeus and M. arctiodes were amplified by using PCR and sequenced. Combined with the same regions of other primate animals, sequence analysis confirmed the conservation of HMGbox. In molecular phylogenetic tree constructed, we made the Ka/Ks value pairwise comparison of HMGbox flanking sequences both intra- and inter- Hominoid and the Old World monkey. To our interest, the pairwise comparison of species in Hominoid suggested a higher Ka/Ks ratio than that in the Old World monkey or between the Old World monkey and marmoset, which showed existent of different pattern of evolution. In addition, in Hominoid Ka/Ks ratio for HMGbox sequences was obviously on the increase. The rate of nonsynonymous substitutions per nonsynonymous site is faster than that of synonymous substitutions per synonymous site in Hominoid (driven by positive directional selection). Our analysis suggested that it might be resulted from functional change of this gene and the system error or the close relationship among Hominoid species. It is also supposed that the positive selection was functioned at the beginning of the divergence of Hominoid species.

Animals↗
PubMed2005

Worldwide haplotype diversity and coding sequence variation at human bitter taste receptor loci.

Bitter taste perception in humans is mediated by receptors encoded by 25 genes that together comprise the TAS2R (or T2R) gene family. The ability to identify the ligand(s) for each of these receptors is dependent on understanding allelic variation in TAS2R genes, which may have a significant effect on ligand recognition. To investigate the extent of coding variation among TAS2R alleles, we performed a comprehensive evaluation of sequence and haplotype variation in the human bitter taste receptor gene repertoire. We found that these genes exhibit substantial coding sequence diversity. In a worldwide population sample of 55 individuals, we found an average of 4.2 variant amino acid positions per gene. In aggregate, the 24 genes analyzed here, along with the phenylthiocarbamide (PTC) receptor gene analyzed previously, specify 151 different protein coding haplotypes. Analyses of the ratio of synonymous and nonsynonymous nucleotide substitutions using the Ka/Ks ratio revealed an excess of amino acid substitutions relative to most other genes examined to date (Ka/Ks = 0.94). In addition, comparisons with more than 1,500 other genes revealed that levels of diversity in the TAS2R genes were significantly greater than expected (pi = 0.11%; p < 0.01), as were levels of differentiation among continental populations (FST = 0.22; p < 0.05). These diversity patterns indicate that unusually high levels of allelic variation are found within TAS2R loci and that human populations differ appreciably with respect to TAS2R allele frequencies. Diversity in the TAS2R genes may be accounted for by natural selection, which may have favored alleles responsive to toxic, bitter compounds found in plants. These findings are consistent with the view that different alleles of the TAS2R genes encode receptors that recognize different ligands, and suggest that the haplotypes we have identified will be important in studies of receptor-ligand recognition.

Alleles↗
PubMed1991

Extreme variations in the ratios of non-synonymous to synonymous nucleotide substitution rates in signal peptide evolution.

Nucleotide sequences encoding signal peptides from the precursors of alpha-amylase/trypsin inhibitors from cereals are homologous to those corresponding to the precursors of thaumatin II and of plastocyanins. Non-synonymous (KA) and synonymous (KS) rates of nucleotide substitutions have been calculated for all possible binary combinations. Extreme variation in KA/KS ratios has been observed; from the 0.167 average found within the plastocyanin family to an average of 1.90 calculated for the inhibitors/thaumatin II transition. A similar calculation has been carried out for the signal peptide sequences of thionins, which are unrelated to those of the alpha-amylase/trypsin inhibitor family, and an average KA/KS of 0.12 has been obtained. This variation can be largely explained in terms of an empirical index of stability related to amino acid composition and seems to be independent of functional constraints.

Amino Acid Sequence↗