A protein alignment scoring system sensitive at all evolutionary distances

Summary Protein sequence alignments generally are constructed with the aid of a substitution matrix that specifies a score for aligning each pair of amino acids. Assuming a simple random protein model, it can be shown that any such matrix, when used for evaluating variable-length local alignments, is implicitly a log-odds matrix, with a specific probability distribution for amino acid pairs to which it is uniquely tailored. Given a model of protein evolution from which such distributions may be derived, a substitution matrix adapted to detecting relationships at any chosen evolutionary distance can be constructed. Because in a database search it generally is not known a priori what evolutionary distances will characterize the similarities found, it is necessary to employ an appropriate range of matrices in order not to overlook potential homologies. This paper formalizes this concept by defining a scoring system that is sensitive at all detectable evolutionary distances. The statistical behavior of this scoring system is analyzed, and it is shown that for a typical protein database search, estimating the originally unknown evolutionary distance appropriate to each alignment costs slightly over two bits of information, or somewhat less than a factor of five in statistical significance. A much greater cost may be incurred, however, if only a single substitution matrix, corresponding to the wrong evolutionary distance, is employed.     

地衣芽孢杆菌2709碱性蛋白酶编码序列的PCR扩增、克隆及序列测定

在地衣芽孢杆菌NCIB 6816菌株碱性蛋白酶基因已知序列的基础上,通过设计合适的引物,利用PCR(Polymerase Chain Reaction)技术从地衣芽孢杆菌2709菌株的柒色体DNA中扩增了2709碱性蛋白酶的编码序列。对两个克隆的PCR片段的全序列分析结果显示,2709碱性蛋白酶的编码序列同相应的NCIB 6816序列相比有3％左右的碱基组成差异。由此推定的2709碱性蛋白酶的氨基酸序列肯定了2709碱性蛋白酶属典型的subtilisin Carlsberg类,同时还表明来源于不同地衣芽孢杆菌菌株的subtilisin Carlsberg存在着若干氨基酸组成上的差异。     

Amino acid sequence of the minor isomorph of the crustacean hyperglycemic hormone (CHH-II) of the Mexican crayfish Procambarus bouvieri (Ortmann): Presence of a d-amino acid

The primary structure of the neurohormone crustacean hyperglycemic hormone (CHH-II) was determined by means of enzymatic digestions, manual Edman degradation, and mass spectrometry. CHH-II is a 72 residue peptide (molecular mass 8388 Da), with six cysteines forming three disulfide bridges that connect residues 7–43, 23–39, and 26–52. The peptide has blocked N- and C-termini, and lacks tryptophan, histidine, and methionine. The CHH-I and CHH-II of Procambarus bouvieri have identical sequences and elicit levels of hyperglycemia that are not distinguishable. The difference between the two isomorphs consists in a posttranslational modification of a l-Phe in CHH-I to a d-Phe in CHH-II at the third position from the N-terminus.     

The size distribution of insertions and deletions in human and rodent pseudogenes suggests the logarithmic gap penalty for sequence alignment

The size distributions of deletions, insertions, and indels (i.e., insertions or deletions) were studied, using 78 human processed pseudogenes and other published data sets. The following results were obtained: (1) Deletions occur more frequently than do insertions in sequence evolution; none of the pseudogenes studied shows significantly more insertions than deletions. (2) Empirically, the size distributions of deletions, insertions, and indels can be described well by a power law, i.e., f _k = Ck ^–b , where f _k is the frequency of deletion, insertion, or indel with gap length k, b is the power parameter, and C is the normalization factor. (3) The estimates of b for deletions and insertions from the same data set are approximately equal to each other, indicating that the size distributions for deletions and insertions are approximately identical. (4) The variation in the estimates of b among various data sets is small, indicating that the effect of local structure exists but only plays a secondary role in the size distribution of deletions and insertions. (5) The linear gap penalty, which is most commonly used in sequence alignment, is not supported by our analysis; rather, the power law for the size distribution of indels suggests that an appropriate gap penalty is w _k = a + b ln k, where a is the gap creation cost and blnk is the gap extension cost. (6) The higher frequency of deletion over insertion suggests that the gap creation cost of insertion (a _i ) should be larger than that of deletion (a _d ); that is, a _i – a _d = In R, where R is the frequency ratio of deletions to insertions. Correspondence to: W.-H. Li     

Mutation Pattern Variation Among Regions of the Primate Genome

We sequenced three argininosuccinate-synthetase-processed pseudogenes (ΨAS-A1, ΨAS-A3, ΨAS-3) and their noncoding flanking sequences in human, orangutan, baboon, and colobus. Our data showed that these pseudogenes were incorporated into the genome of the Old World monkeys after the divergence of the Old World and New World monkey lineages. These pseudogene flanking regions show variable mutation rates and patterns. The variation in the G/C to A/T mutation rate (u) can account for the unequal GC contents at equilibrium: 34.9, 36.9, and 41.7% in the pseudogene ΨAS-A1, ΨAS-A3, and ΨAS-3 flanking regions, respectively. The A/T to G/C mutation rate (v) seems stable and the u/v ratios equal 1.9, 1.7, and 1.4 in the flanking regions of ΨAS-A1, ΨAS-A3, and ΨAS-3, respectively. These ``regional' variations of the mutation rate affect the evolution of the pseudogenes, too. The ratio u/v being greater than 1.0 in each case, the overall mutation rate in the GC-rich pseudogenes is, as expected, higher than in their GC-poor flanking regions. Moreover, a ``sequence effect' has been found. In the three cases examined u and v are higher (at least 20%) in the pseudogene than in its flanking region—i.e., the pseudogene appears as mutation ``hot' spots embedded in ``cold' regions. This observation could be partly linked to the fact that the pseudogene flanking regions are long-standing unconstrained DNA sequences, whereas the pseudogenes were relieved of selection on their coding functions only around 30–40 million years ago. We suspect that relatively more mutable sites maintained unchanged during the evolution of the argininosuccinate gene are able to change in the pseudogenes, such sites being eliminated or rare in the flanking regions which have been void of strong selective constraints over a much longer period. Our results shed light on (1) the multiplicity of factors that tune the spontaneous mutation rate and (2) the impact of the genomic position of a sequence on its evolution. Received: 10 February 1997 / Accepted: 21 April 1997     

Molecular Diversity and Evolution of bla TEM Genes Encoding β-Lactamases Resistant to Clavulanic Acid in Clinical E. coli

The molecular diversity of inhibitor-resistant TEM (IRT) enzymes was explored using a strategy which involved DNA amplification by polymerase chain reaction (PCR), analysis of restriction fragment length polymorphism (RFLP), and direct nucleotide sequencing. The study of plasmid-borne genes from 27 strains, resistant to amoxicillin and β-lactamase-inhibitor combinations, identified mutations resulting in amino acid change at positions 69, 244, 275, and 276 known to be associated with the IRT phenotype and a mutation at nucleotide position 162 in the promoter region. These mutations were found to lie on two different gene sequences, described here as ``TEM-1B like' and ``TEM-2 like' restriction linkage groups. Further analysis, of nucleotide sequences of promoter and coding regions of the β-lactamases, confirmed that a given mutation causing IRT phenotype could be associated with two different gene sequence frameworks and two different causal mutations could lie on identical gene sequence framework. These data argue in favor of convergent phenotypic evolution of IRT enzymes under the selective pressure imposed by the intensive clinical use of β-lactam–β-lactamase inhibitor combinations. Received: 18 March 1996 / Accepted: 15 July 1996     

Molecular Evolution and Phylogeny of Satellite RNA Associated with Bamboo Mosaic Potexvirus

Satellite RNA of bamboo mosaic potexvirus (satBaMV) is a linear RNA molecule which encodes a 20-kDa nonstructural protein. Sequences of seven different satBaMV isolates from bamboo hosts in three genera showed 0.7% to 7.5% base variation which spanned the whole RNA molecule. However, the putative 20-kDa open reading frame was all preserved in these isolates. The phylogenetic relationship based on the nucleotide sequence did not show particular grouping of satBaMV from the host in one genus; neither was the grouping of satBaMV evident by location of sampling. Putative secondary structures of the 3′ untranslated regions showed a basic pattern with conserved hexanucleotides (ACCUAA) and polyadenylation signal (AAUAAA) located in the loop regions. Although the satBaMV-encoded 20-kDa protein is a nonstructural protein, its predicted secondary structure contains eight-stranded β-sheets which may form ``jelly-roll' structure similar to that found in capsid protein encoded by satellite virus of panicum mosaic virus. Received: 26 June 1996 / Accepted: 9 September 1996