The molecular evolution of pancreatic ribonuclease

Summary The primary structures of pancreatic ribonucleases from 26 species (18 artiodactyls, horse, whale, 5 rodents and turtle) are known. Several species contain identical ribonucleases (cow/bison; sheep/goat), other species show polymorphism (arabian camel) or the presence of two structural gene loci (guinea pig pancreas contains two ribonucleases that differ at 31 positions). 26 different sequences (including the ribonuclease from bovine seminal plasma which is paralogous to the pancreatic ribonucleases) were used to construct a most parsimonious tree. A second tree that most closely approximates current biological opinion requires 402 whereas the most parsimonious tree requires 389 nucleotide substitutions. The artiodactyl part of the most parsimonious tree conforms quite well with the biological one of this order, except for the position of the giraffe which is placed with the pronghorn. Other parts of the most parsimonious tree agree less with the biological tree, probably as a result of the occurrence of many parallel and back substitutions. Bovine seminal ribonuclease was found to be the result of a gene duplication which occurred before the divergence of the true ruminants, but after the divergence of this group from the cameloids.The evolutionary rate of ribonuclease was found to be 390, 3.0 and 11 nucleotide substitutions per 10⁹ yrs per ribonuclease gene, codon and covarion respectively. However, there is much variation in evolutionary rate in different taxa. Values ranging from about 100 (in the bovidae) to about 700 (in the rodents) nucleotide substitutions per 10⁹ yrs per gene were found.A method for counting parallel and back mutations is presented. The 389 nucleotide substitutions in the most parsimonious tree occur at 88 codon positions; 154 of them are the result of parallel and back mutations. Parallel evolution to a similar structure, including the presence of 2 sites with carbohydrate, was demonstrated in an extensive region at the surface of pig and guinea pig ribonuclease B. The presence of carbohydrate probably is important in a number of species. A correlation between the presence of heavily glycosidated ribonucleases and coecal digestion was observed. Hypothetical sequences of ancestral ungulate ribonucleases contain many recognition sites for carbohydrate attachment; this suggests that herbivores with coecal digestion might have preceded the true ruminants in mammalian evolution.     

Origin of the duplicated ribonuclease gene in guinea-pig: Comparison of the amino acid sequences with those of two close relatives: Capybara and cuis ribonuclease

The amino acid sequences of the pancreatic ribonuclease from capybara (Hydrochoerus hydrochaeris) and cuis (Galea musteloides) were determined. Both species belong to the same superfamily of the hystricomorph rodents as the guinea-pig. In guinea-pig pancreas two ribonucleases are present as a result of a recent gene duplication, but in capybara and cuis pancreas only one single ribonuclease has been found. A most parsimonious tree of ribonucleases indicates that the gene duplication leading to both guinea-pig ribonucleases occurred before the divergence of guinea-pig from capybara and cuis. This would mean that changes in expression of the ribonuclease genes have occurred in these taxa. Cuis and capybara ribonuclease have no Asn-X-Ser/Thr sequences and are carbohydrate-free proteins. Capybara ribonuclease has leucine at position 114, a position occupied by proline in the cis-configuration in bovine pancreatic ribonuclease.     

Secretory ribonuclease genes and pseudogenes in true ruminants

Mammalian pancreatic ribonucleases (RNase) form a family of extensively studied homologous proteins. Phylogenetic analyses, based on the primary structures of these enzymes, indicated that the presence of three homologous enzymes (pancreatic, seminal and brain ribonucleases) in the bovine species is due to gene duplication events, which occurred during the evolution of ancestral ruminants. In this paper the sequences are reported of the coding regions of the orthologues of the three bovine secretory ribonucleases in hog deer and roe deer, two deer species belonging to two different subfamilies of the family Cervidae. The sequences of the 3′ untranslated regions of the three different secretory RNase genes of these two deer species and giraffe are also presented. Comparison of these and previously determined sequences of ruminant ribonucleases showed that the brain-type enzymes of giraffe and these deer species exhibit variations in their C-terminal extensions. The seminal-type genes of giraffe, hog deer and roe deer show all the features of pseudogenes. Phylogenetic analyses, based on the complete coding regions and parts of the 3′ untranslated regions of the three different secretory ribonuclease genes of ox, sheep, giraffe and the two deer species, show that pancreatic, seminal- and brain-type RNases form three separate groups.     

Molecular evolution of mammalian ribonucleases 1

There have been many studies on the chemistry of mammalian pancreatic ribonucleases (ribonucleases 1), but the functional biology of this family of homologous proteins is still largely unknown. Many studies have been performed on the molecular evolution and properties of this enzyme from species belonging to a large number of mammalian taxa, including paralogous gene products resulting from recent gene duplications. Novel ribonuclease 1 sequences were determined for three rodent species (gundi, brush-tailed porcupine, and squirrel), rabbit, a fruit bat, elephant, and aardvark, and the new sequences were used for deriving most parsimonious networks of ribonucleases from different mammalian orders, including earlier determined nucleotide sequences and also a larger set of protein sequences. Weak support for interordinal relationships were obtained, except for an Afrotheria clade containing elephant and aardvark. Results of current analyses and also those obtained 20 years ago on amino acid sequences confirm conclusions derived recently from larger data sets of other molecules. Several examples of recent gene duplications in ribonucleases 1 are discussed, with respect to illustrate the concepts of orthology and paralogy. Previously evidence was presented for extensive parallelism between sequence regions with attached carbohydrate (about one quarter of the molecule) of unrelated species with cecal digestion (pig and guinea pig). These features are also present in the sequences of elephant and fruit bat, species with cecal digestion, but with a very low ribonuclease content in their pancreas.     

Phylogeny of ruminants secretory ribonuclease gene sequences of pronghorn (Antilocapra americana)

Phylogenetic analyses based on primary structures of mammalian ribonucleases, indicated that three homologous enzymes (pancreatic, seminal and brain ribonucleases) present in the bovine species are the results of gene duplication events, which occurred in the ancestor of the ruminants after divergence from other artiodactyls. In this paper sequences are presented of genes encoding pancreatic and brain-type ribonuclease genes of pronghorn (Antilocapra americana). The seminal-type ribonuclease gene could not be detected in this species, neither by PCR amplification nor by Southern blot analyses, indicating that it may be deleted completely in this species. Previously we demonstrated of a study of amino acid sequences of pancreatic ribonucleases of a large number of ruminants the monophyly of bovids and cervids, and that pronghorn groups with giraffe. Here we present phylogenetic analyses of nucleotide sequences of ribonucleases and other molecules from ruminant species and compare these with published data. Chevrotain (Tragulus) always groups with the other ruminants as separate taxon from the pecora or true ruminants. Within the pecora the relationships between Bovidae, Cervidae, Giraffidae, and pronghorn (Antilocapra) cannot be decided with certainty, although in the majority of analyses Antilocapra diverges first, separately or joined with giraffe. Broad taxon sampling and investigation of specific sequence features may be as important for reliable conclusions in phylogeny as the lengths of analyzed sequences.     

Primary structures of pancreatic ribonucleases from Bovidae: Impala,Thomson's gazelle,nilgai and water buffalo

The amino acid sequences of the pancreatic ribonucleases from impala, Thomson's gazelle, nilgai and two types of water buffalo were studied from several enzymic digests. Peptides were positioned by homology with other ribonucleases of known sequence. Only peptides that differ in amino acid composition from the corresponding peptides of ox and goat ribonuclease were sequenced.The primary structures of pancreatic ribonucleases from 11 species of the Bovidae family are known to date. The evolutionary rate of bovid ribonucleases is 2–3 times lower than the average rate observed in all mammals. A possible explanation is that the presence of a stable symbiotic system in the rumen of grass-eating ruminants has caused a slowing down of the evolutionary rate of pancreatic ribonuclease in this taxon. The subfamily of the Bovinae (five species) exhibits a slightly higher evolutionary rate; most replacements in the Bovinae occur near residues 34–36, the most commonly observed carbohydrate attachment site in ribonucleases.     

Sequences related to the ox pancreatic ribonuclease coding region in the genomic DNA of mammalian species

Mammalian pancreatic ribonucleases form a family of homologous proteins that has been extensively investigated. The primary structures of these enzymes were used to derive phylogenetic trees. These analyses indicate that the presence of three strictly homologous enzymes in the bovine species (the pancreatic, seminal, and cerebral ribonucleases) is due to gene duplication events which occurred during the evolution of ancestral ruminants.In this paper we present evidence that confirms this finding and that suggests an overall structural conservation of the putative ribonuclease genes in ruminant species.We could also demonstrate that the sequences related to ox ribonuclease coding regions present in genomic DNA of the giraffe species are the orthologues of the bovine genes encoding the three ribonucleases mentioned above.Correspondence to: A. Furia     

Steady-state enzyme kinetics of the pancreatic ribonucleases from five mannalian species.

The kinetic parameters Km, k+2 and k+2/Km of the pancreatic ribonucleases (EC 3.1.4.22) from cow, giraffe, horse, rat and lesser rorqual have been determined, using 2',3'-cyclic cytidine monophosphate and 2',3'-cuclic uridine monophosphate as substrates. No large differences were found between the activities of the five enzymes. The relative differences between the activities of the five enzymes are mainly due to differences in the rates of hydrolysis and not to differences in the affinities for the substrates.     

Amino acid sequence of human von Willebrand factor

The complete amino acid sequence of human von Willebrand factor (vWF) is presented. Most of the sequence was determined by analysis of the S-carboxymethylated protein. Some overlaps not provided by the protein sequence analysis were obtained from the sequence predicted by the nucleotide sequence of a cDNA clone [Sadler, J.E., Shelton-Inloes, B.B., Sorace, J., Harlan, M., Titani, K., & Davie, E.W. (1985) Proc. Natl. Acad. Sci. U.S.A. 82, 6391-6398]. The protein is composed of 2050 amino acid residues containing 12 Asn-linked and 10 Thr/Ser-linked oligosaccharide chains. One of the carbohydrate chains is linked to an Asn residue in the sequence Asn-Ser-Cys rather than the usual Asn-X-Ser/Thr sequence. The sequence of von Willebrand factor includes several regions bearing evidence of internal gene duplication of ancestral sequences. The protein also contains the tetrapeptide sequence Arg-Gly-Asp-Ser (at residues 1744-1747), which may be a cell attachment site, as in fibronectin. The amino- and carboxyl-terminal regions of the molecule contain clusters of half-cystinyl residues. The sequence is unique except for some homology to human complement factor B.     

Sequence differences between glycosylated and non-glycosylated Asn-X-Thr/Ser acceptor sites: implications for protein engineering

In N-glycosylated glycoproteins, carbohydrate is attached to Asn in the sequence Asn-X-Ser/Thr, where X denotes any amino acid. However, the presence of this consensus peptide does not always lead to glycosylation. We have compiled an extensive collection of glycosylated and non-glycosylated Asn-X-Thr/Ser sites and present a statistical study based on this data set. Our results indicate that non-glycosylated sites tend to be found more frequently towards the C termini of glycoproteins, and that proline residues in positions X and Y in the consensus Asn-X-Thr/Ser-Y strongly reduce the likelihood of N-linked glycosylation. Beyond this, there are no obvious local sequence features that seem to correlate with the absence or presence of N-linked glycosylation. These findings are discussed in terms of the prediction and engineering of glycosylation sites in secretory proteins.     

Natural human interferon-gamma. Complete amino acid sequence and determination of sites of glycosylation

Fresh human peripheral blood lymphocytes were induced with desacetylthymosin -alpha 1 and staphylococcal enterotoxin B. The induced gamma interferon (or IFN-gamma, immune interferon, type II interferon) was purified to homogeneity utilizing controlled-pore glass, concanavalin A-Sepharose, Bio-Gel P100, or Sephacryl S-200, and reversed phase high performance liquid chromatography. This procedure resulted in two active species with apparent Mr = 20,000 and 25,000 as judged by sodium dodecyl sulfate-polyacrylamide gel electrophoresis. Both species were found to have identical amino acid sequences with a pyroglutamate residue as NH2-terminus. In both cases six different COOH termini were found. They are, at least qualitatively, identical in both species. There are two possible Asn-X-Ser/Thr glycosylation sites. Both carry carbohydrates in the Mr = 25,000 species whereas in the Mr = 20,000 species only one site is glycosylated. This likely explains the difference in apparent molecular weight between the two species and the expected molecular weight based upon the amino acid sequence.     

Naturally occurring structural isomers in serum IgA1 o-glycosylation

IgA is the most abundantly produced antibody and plays an important role in the mucosal immune system. Human IgA is represented by two isotypes, IgA1 and IgA2. The major structural difference between these two subclasses is the presence of nine potential sites of O-glycosylation in the hinge region between the first and second constant region domains of the heavy chain. Thr(225), Thr(228), Ser(230), Ser(232) and Thr(236) have been identified as the predominant sites of O-glycan attachment. The range and distribution of O-glycan chains at each site within the context of adjacent sites in this clustered region create a complex heterogeneity of surface epitopes that is incompletely defined. We previously described the analysis of IgA1 O-glycan heterogeneity by use of high resolution LC-MS and electron capture dissociation tandem MS to unambiguously localize all amino acid attachment sites in IgA1 (Ale) myeloma protein. Here, we report the identification and elucidation of IgA1 O-glycopeptide structural isomers that occur based on amino acid position of the attached glycans (positional isomers) and the structure of the O-glycan chains at individual sites (glycan isomers). These isomers are present in a model IgA1 (Mce1) myeloma protein and occur naturally in normal human serum IgA1. Variable O-glycan chains attached to Ser(230), Thr(233) or Thr(236) produce the predominant positional isomers, including O-glycans composed of a single GalNAc residue. These findings represent the first definitive identification of structural isomeric IgA1 O-glycoforms, define the single-site heterogeneity for all O-glycan sites in a single sample, and have implications for defining epitopes based on clustered O-glycan variability.     

Complete sequence of the Rous sarcoma virus env gene: identification of structural and functional regions of its product.

The amino-terminal amino acid sequences of gp85 and gp37, the envelope glycoproteins of Rous sarcoma virus (RSV), were determined. Alignment of these sequences with the amino acid sequence predicted from the complete nucleotide sequence of the Prague strain of RSV, subgroup C (PR-C), has allowed us to delineate the env gene-coding region of this virus. The coding sequences for gp85 and gp37 have been placed in an open reading frame that extends from nucleotide 5045 to nucleotide 6862 and predict sizes of 341 amino acids (36,962 molecular weight) for gp85 and 198 amino acids (21,566 molecular weight) for gp37. Carbohydrate makes a significant contribution to the observed molecular weights of these polypeptides--the amino acid sequence contains 14 potential glycosylation sites (Asn-X-Ser/Thr) in gp85 and two in gp37. Experiments aimed at estimating the number of carbohydrate side chains yielded results consistent with most or all of these sites being occupied. Although an initiation codon is located early (codon 4) in the open reading frame, it is likely that splicing yields an mRNA on which translation initiates at the same AUG as that of the gag gene to produce a nascent polypeptide in which gp85 is preceded by a 62-amino-acid-long leader peptide. This leader contains the hydrophobic sequence (signal sequence) necessary for translocation across the endoplasmic reticulum and is completely removed from the env gene product during translation. The polyprotein precursor, Pr95env, is cleaved to gp85 and gp37 at the carboxyl side of the basic sequence:-Arg-Arg-Lys-Arg-. gp85 is attached through a disulphide linkage to gp37, and although the positions of the cysteines involved in this linkage are not known, the presence of a 27-amino-acid-long hydrophobic region at the carboxy-terminus of gp37 is consistent with its role as a membrane anchor for the viral glycoprotein complex. The location of host range variable regions with respect to the possible tertiary structure of the complex is discussed.     

Amino acid sequence of the nonsecretory ribonuclease of human urine

The amino acid sequence of a nonsecretory ribonuclease isolated from human urine was determined except for the identity of the residue at position 7. Sequence information indicates that the ribonucleases of human liver and spleen and an eosinophil-derived neurotoxin are identical or very closely related gene products. The sequence is identical at about 30% of the amino acid positions with those of all of the secreted mammalian ribonucleases for which information is available. Identical residues include active-site residues histidine-12, histidine-119, and lysine-41, other residues known to be important for substrate binding and catalytic activity, and all eight half-cystine residues common to these enzymes. Major differences include a deletion of six residues in the (so-called) S-peptide loop, insertions of two, and nine residues, respectively, in three other external loops of the molecule, and an addition of three residues at the amino terminus. The sequence shows the human nonsecretory ribonuclease to belong to the same ribonuclease superfamily as the mammalian secretory ribonucleases, turtle pancreatic ribonuclease, and human angiogenin. Sequence data suggest that a gene duplication occurred in an ancient vertebrate ancestor; one branch led to the nonsecretory ribonuclease, while the other branch led to a second duplication, with one line leading to the secretory ribonucleases (in mammals) and the second line leading to pancreatic ribonuclease in turtle and an angiogenic factor in mammals (human angiogenin). The nonsecretory ribonuclease has five short carbohydrate chains attached via asparagine residues at the surface of the molecule; these chains may have been shortened by exoglycosidase action.(ABSTRACT TRUNCATED AT 250 WORDS)     

Inclusion of Cetaceans Within the Order Artiodactyla Based on Phylogenetic Analysis of Pancreatic Ribonuclease Genes

Mammalian secretory ribonucleases (RNases 1) form a family of extensively studied homologous proteins that were already used for phylogenetic analyses at the protein sequence level previously. In this paper we report the determination of six ribonuclease gene sequences of Artiodactyla and two of Cetacea. These sequences have been used with ruminant homologues in phylogenetic analyses that supported a group including hippopotamus and toothed whales, a group of ruminant pancreatic and brain-type ribonucleases, and a group of tylopod sequences containing the Arabian camel pancreatic ribonuclease gene and Arabian and Bactrian camel and alpaca RNase 1 genes of unknown function. In all analyses the pig was the first diverging artiodactyl. This DNA-based tree is compatible to published trees derived from a number of other genes. The differences to those trees obtained with ribonuclease protein sequences can be explained by the influence of convergence of pancreatic RNases from hippopotamus, camel, and ruminants and by taking into account the information from third codon positions in the DNA-based analyses. The evolution of sequence features of ribonucleases such as the distribution of positively charged amino acids and of potential glycosylation sites is described with regard to increased double-stranded RNA cleavage that is observed in several cetacean and artiodactyl RNases which may have no role in ruminant or ruminant-like digestion. Received: 2 June 1998 / Accepted: 31 August 1998     

Evolution of glycoproteins as judged by the frequency of occurrence of the tripeptides Asn-X-Ser and Asn-X-Thr in proteins

Summary A statistical analysis of the data tabulated in the Atlas of Protein Sequence and Structure 1972 indicates that the observed frequency of occurrence of the tripeptides Asn-X-Ser and Asn-X-Thr is approximately one third of the expected in eukaryotic proteins, but in prokaryotic proteins the observation agrees closely with expectation. Thus the lowered frequency of these tripeptides found by Hunt and Dayhoff is restricted to eukaryotic proteins. Of all the Asn-X-Ser/Thr sequences examined, those which contain covalently attached carbohydrates are found only in the extracellular proteins of eukaryote. These observations are discussed in relation to the evolution of glycoproteins which seems to have occurred in the ancestor of eukaryotes after the divergence from prokaryotes.     

Primary structure of pronghorn pancreatic ribonuclease: Close relationship between giraffe and pronghorn

Summary Pancreatic ribonuclease from pronghorn (Antilocapra americana) was isolated and its amino acid sequence was determined from a tryptic digest of the performic acid-oxidized protein. Peptides were positioned by homology with other ribonucleases. Only peptides that differed in amino acid composition from the corresponding peptides of ox or goat ribonucleases were sequenced.In a most parsimonious tree of pancreatic ribonucleases, pronghorn and giraffe were placed together and these two were placed with the bovids, leaving the deer as a taxon separate from the other ruminants. The amino acid replacements that determine this tree topology are three rarely occurring replacements shared by pronghorn and giraffe. Notwithstanding their close phylogenetic relationship, both ribonucleases differ strongly in extent of glycosidation, net charge and antigenic properties.     

The primary structure of muskrat pancreatic ribonuclease.

Pancreatic ribonuclease from muskrat (Ondatra zibethica) was isolated and its amino acid sequence was determined from tryptic digests of the performic acid-oxidized and the reduced and aminoethylated enzyme. The peptides have been positioned in the sequence by homology with other ribonucleases. This could be done unambiguously for all peptides except Arg-Arg (tentative position 32-33) and Ser-Arg (tentative position 75-76). The amino acid sequences of the peptides were determined by the dansyl-Edman method, with the exception of residues 23-25 and 99-102, which were positioned by homology. The enzyme differs in 38 positions from the enzyme from rat and in 31-42 positions from other mammalian pancreatic ribonucleases, while rat ribonuclease differs at 44-52 positions from the other enzymes. These data point to a common ancestry of the enzymes from muskrat and rat and an increased evolution rate of rat ribonuclease after divergence of the ancestors of both species. Muskrat ribonuclease contains no carbohydrate, although the enzyme possesses a recognition site for carbohydrate attachment in the sequence Asn-Val-Thr (62-64).     

The amino-acid sequence of kangaroo pancreatic ribonuclease

Red kangaroo (Macropus rufus) ribonuclease was isolated from pancreatic tissue by affinity chromatography. The amino acid sequence was determined by automatic sequencing of overlapping large fragments and by analysis of shorter peptides obtained by digestion with a number of proteolytic enzymes. The polypeptide chain consists of 122 amino acid residues. Compared to other ribonucleases, the N-terminal residue and residue 114 are deleted. In other pancreatic ribonucleases position 114 is occupied by a cis proline residue in an external loop at the surface of the molecule. Other remarkable substitutions are the presence of a tyrosine residue at position 123 instead of a serine which forms a hydrogen bond with the pyrimidine ring of a nucleotide substrate, and a number of hydrophobichydrophilic interchanges in the sequence 51-55, which forms part of an alpha-helix in bovine ribonuclease and exhibits few substitutions in the placental mammals. Kangaroo ribonuclease contains no carbohydrate, although the enzyme possesses a recognition site for carbohydrate attachment in the sequence Asn-Val-Thr (62-64). The enzyme differs at about 35-40% of the positions from all other mammalian pancreatic ribonucleases sequenced to date, which is in agreement with the early divergence between the marsupials and the placental mammals. From fragmentary data a tentative sequence of red-necked wallaby (Macropus rufogriseus) pancreatic ribonuclease has been derived. Eight differences with the kangaroo sequence were found.