Episodic evolution in the stomach lysozymes of ruminants

Summary By sequencing lysozymesc from deer and pig stomachs and comparing them to the known amino acid sequences of other lysozymesc, it was possible to examine the rate of sequence change during and after the period in which this enzyme acquired a new function. Evolutionary tree analysis suggests that the rate went up while lysozyme was being recruited to function as a digestive enzyme in the stomach of early ruminants. Later, presumably after lysozyme was well adapted for functioning in the new environment, which contains acid, pepsin, and fermentation products, the rate of amino acid replacement became subnormal.     

Molecular genetics and evolution of stomach and nonstomach lysozymes in the hoatzin

Multiple genes of the hoatzin encoding stomach lysozyme c and closely related members of this calcium-binding lysozyme c group were cloned from a genomic DNA library and sequenced. There are a minimum of five genes represented among these sequences that encode two distinct groups of protein sequences. One group of three genes corresponds to the stomach lysozyme amino acid sequences, and the remaining genes encode predicted proteins that are more basic in character and share several sequence identities with the pigeon egg-white lysozyme rather than with the hoatzin stomach lysozymes. Despite these structural similarities between some of the hoatzin gene products and the pigeon lysozyme, phylogenetic analyses indicate that all of the hoatzin sequences are closely related to one another. This is borne out by the relatively small genetic distances even in the intronic regions, which are not subject to the selective pressures operating on the coding regions of the stomach lysozymes. These results suggest that multiple gene duplication events have occurred during the evolution of hoatzin lysozymes.     

cDNA and amino acid sequences of rainbow trout (Oncorhynchus mykiss) lysozymes and their implications for the evolution of lysozyme and lactalbumin

Summary The complete 129-amino-acid sequences of two rainbow trout lysozymes (I and II) isolated from kidney were established using protein chemistry microtechniques. The two sequences differ only at position 86, I having aspartic acid and II having alanine. A cDNA clone coding for rainbow trout lysozyme was isolated from a cDNA library made from liver mRNA. Sequencing of the cloned cDNA insert, which was 1 kb in length, revealed a 432-bp open reading frame encoding an amino-terminal peptide of 15 amino acids and a mature enzyme of 129 amino acids identical in sequence to II. Forms I and II from kidney and liver were also analyzed using enzymatic amplification via PCR and direct sequencing; both organs contain mRNA encoding the two lysozymes. Evolutionary trees relating DNA sequences coding for lysozymesc and α-lactalbumins provide evidence that the gene duplication giving rise to conventional vertebrate lysozymesc and to lactalbumin preceded the divergence of fishes and tetrapods about 400 Myr ago. Evolutionary analysis also suggests that amino acid replacements may have accumulated more slowly on the lineage leading to fish lysozyme than on those leading to mammal and bird lysozymes.     

Amino acid sequences of ovomucoid third domains from 27 additional species of birds

Ovomucoids consist of a single polypeptide chain which is composed of three tandem Kazal domains. Each Kazal domain is an actual or putative protein inhibitor of serine proteinases. Ovomucoid third domains were already isolated and sequenced from 126 species of birds (Laskowskiet al., 1987, 1990). This paper adds 27 new species. A number of generalizations are made on the basis of sequences from 153 species. The residues that are in contact with the enzyme in enzyme-inhibitor complexes are strikingly hypervariable. While the primary specificity residue,P ₁, is the most variable; substitutions occur predominantly among aliphatic, hydrophobic residues. Consensus sequences for an avian ovomucoid third domain, for a b-type Kazal domain (i.e., a COOH terminal domain of multidomain inhibitors) and for a general Kazal domain are given. Finally, the individual new sequences are briefly discussed.     

Effectiveness of measures requiring and not requiring prior sequence alignment for estimating the dissimilarity of natural sequences

Summary Various measures of sequence dissimilarity have been evaluated by how well the additive least squares estimation of edges (branch lengths) of an unrooted evolutionary tree fit the observed pairwise dissimilarity measures and by how consistent the trees are for different data sets derived from the same set of sequences. This evaluation provided sensitive discrimination among dissimilarity measures and among possible trees. Dissimilarity measures not requiring prior sequence alignment did about as well as did the traditional mismatch counts requiring prior sequence alignment. Application of Jukes-Cantor correction to singlet mismatch counts worsened the results. Measures not requiring alignment had the advantage of being applicable to sequences too different to be critically alignable. Two different measures of pairwise dissimilarity not requiring alignment have been used: (1) multiplet distribution distance (MDD), the square of the Euclidean distance between vectors of the fractions of base singlets (or doublets, or triplets, or…) in the respective sequences, and (2) complements of long words (CLW), the count of bases not occurring in significantly long common words. MDD was applicable to sequences more different than was CLW (noncoding), but the latter often gave better results where both measures were available (coding). MDD results were improved by using longer multiplets and, if the sequences were coding, by using the larger amino acid and codon alphabets rather than the nucleotide alphabet. The additive least squares method could be used to provide a reasonable consensus of different trees for the same set of species (or related genes).     

Amino acid sequences of ovomucoid third domain from 25 additional species of birds

Ovomucoids were isolated from 25 avian species other than the 101 studied in Laskowskiet al. (1987,Biochemistry 26, 202–221). These were subjected to limited proteolysis with an appropriate enzyme, and connecting peptide extended ovomucoid third domains were isolated and sequenced to the end in a protein sequencer. Of the 25 new sequences, 13 duplicate ones were already known, and 12 are unique. Probably the most striking findings are a Pro¹⁴ Ser¹⁴ replacement in weka, an Ala¹⁴Thr¹⁵ replacement in Bulwer's pheasant, the discovery of two additional amino acid residues Ile¹⁸ and Gly¹⁸ at the P₁ reactive site position in Kalij pheasant and tawny frogmouth, respectively, and the first finding of a negative (Glu³⁴) rather than positive (Lys³⁴ or Arg³⁴) amino acid residue at the NH₂ terminus of the helix in caracara ovomucoid third domain. These results complete the determination of all the sequences of ovomucoid third domains in the four species genusGallus, in the five species genusSyrmaticus, and in the two species generaAix andPavo.     

Comparison of goose-type,chicken-type,and phage-type lysozymes illustrates the changes that occur in both amino acid sequence and three-dimensional structure during evolution

Summary The three-dimensional structure of goose-type lysozyme (GEWL), determined by x-ray crystallography and refined at high resolution, has similarities to the structures of hen (chicken) eggwhite lysozyme (HEWL) and bacteriophage T4 lysozyme (T4L). The nature of the structural correspondence suggests that all three classes of lysozyme diverged from a common evolutionary precursor, even though their amino acid sequences appear to be unrelated (Grütter et al. 1983).In this paper we make detailed comparisons of goose-type, chicken-type, and phage-type lysozymes. The lysozymes have undergone conformational changes at both the blobal and the local level. As in the globins, there are corresponding -helices that have rigid-body displacements relative to each other, but in some cases corresponding helices have increased or decreased in length, and in other cases there are helices in one structure that have no counterpart in another.Independent of the overall structural correspondence among the three lysozyme backbones is another, distinct correspondence between a set of three consecutive -helices in GEWL and three consecutive -helices in T4L. This structural correspondence could be due, in part, to a common energetically favorable contact between the first and the third helices.There are similarities in the active sites of the three lysozymes, but also one striking difference. Glu 73 (GEWL) spatially corresponds to Glu 35 (HEWL) and to Glu 11 (T4L). On the other hand, there are two aspartates in the GEWL active site, Asp 86 and Asp 97, neither of which corresponds exactly to Asp 52 (HEWL) or Asp 20 (T4L). (The discrepancy in the location of the carboxyl groups is about 10 Å for Asp 86 and 4 Å for Asp 97.) This lack of structural correspondence may reflect some differences in the mechanisms of action of three lysozymes. When the amino acid sequences of the three lysozyme types are aligned according to their structural correspondence, there is still no apparent relationship between the sequences except for possible weak matching in the vicinity of the active sites.     

Ancient origin of lactalbumin from lysozyme: Analysis of DNA and amino acid sequences

Summary Parsimony trees relating DNA sequences coding for lysozymesc and -lactalbumins suggest that the gene duplication that allowed lactalbumin to evolve from lysozyme preceded the divergence of mammals and birds. Comparisons of the amino acid sequences of additional lysozymes and lactalbumins are consistent with this view. When all base positions are considered, the probability that the duplication leading to the lactalbumin gene occurred after the start to mammalian evolution is estimated to be 0.05–0.10. Elimination of the phylogenetic noise generated by fast evolution and compositional bias at third positions of codons reduced this probability to 0.002–0.03. Thus the gene duplication may have long preceded the acquisition of lactalbumin function.     

Amino acid sequences of P1 protamines and the phylogeny of eutherian mammals: a cladistic study

Amino acid and cDNA sequences of eutherian PI protamines, known from publications of other authors, were compared by a cladistic method. Fish, toad and bird protamines were used for the pertinent “outgroup comparisons”, i.e. they provided relevant data for the comparative alignment of the sequences and for the recognition of evolutionary trends. In the sequence positions compared, each amino acid was individually assigned as a plesiomorphic or apomorphic character state (qualitative treatment). The resulting phylogenetic tree (Fig. 2) is only partially in accordance with common ideas on eutherian phylogeny. Disagreements refer to the branching points of Perissodactyla, Lagomorpha and Rodentia.     

Characterization of the cow stomach lysozyme genes: Repetitive DNA and concerted evolution

Cow stomach lysozyme genes have evolved in a mosaic pattern. The majority of the intronic and flanking sequences show an amount of sequence difference consistent with divergent evolution since duplication of the genes 40–50 million years ago. In contrast, exons 1, 2, and 4 and immediately adjacent intronic sequences differ little between genes and show evidence of recent concerted evolution. Exon 3 appears to be evolving divergently. The three characterized genes vary from 5.6 to 7.9 kilobases in length. Different distributions of repetitive DNA are found in each gene, which accounts for the majority of length differences between genes. The different distributions of repetitive DNA in each gene suggest the repetitive elements were inserted into each gene after the duplications that give rise to these three genes and provide additional support for divergent evolution for the majority of each gene. The observation that intronic and flanking sequences are evolving divergently suggests that the concerted evolution events involved in homogenizing the coding regions of lysozyme genes involve only one exon at a time. This model of concerted evolution would allow the shuffling of exon-sized pieces of information between genes, a phenomenon that may have aided in the early adaptive evolution of stomach lysozyme.Deceased July 21, 1991 Correspondence to: D.M. Irwin     

The alignment of sets of sequences and the construction of phyletic trees: An integrated method

Summary In this paper we argue that the alignment of sets of sequences and the construction of phyletic trees cannot be treated separately. The concept of good alignment is meaningless without reference to a phyletic tree, and the construction of phyletic trees presupposes alignment of the sequences.We propose an integrated method that generates both an alignment of a set of sequences and a phyletic tree. In this method a putative tree is used to align the sequences and the alignment obtained is used to adjust the tree; this process is iterated. As a demonstration we apply the method to the analysis of the evolution of 5S rRNA sequences in prokaryotes.     

Compositional statistics: An improvement of evolutionary parsimony and its application to deep branches in the tree of life

Summary We present compositional statistics, a new method of phylogenetic inference, which is an extension of evolutionary parsimony. Compositional statistics takes account of the base composition of the compared sequences by using nucleotide positions that evolutionary parsimony ignores. It shares with evolutionary parsimony the features of rate invariance and the fundamental distinction between transitions and transversions. Of the presently available methods of phylogenetic inference, compositional statistics is based on the fewest and mildest assumptions about the mode of DNA sequence evolution. It is therefore applicable to phylogenetic studies of the most distantly related organisms or molecules. This was illustrated by analyzing conservative positions in the DNA sequences of the large subunit of RNA polymerase from three archaebacterial groups, a eubacterium, a chloroplast, and the three eukaryotic polymerases. Internally consistent results, which are in accord with our knowledge of organelle origin and archaebacterial physiology, were achieved.