Molecular evolution of genes encoding ribonucleases in ruminant species

Evolution of Digestive Enzymes and RNASE1 Provides Insights into Dietary Switch of Cetaceans

Although cetaceans (whales, porpoises, and dolphins) have multi-chambered stomachs, feeding habits of modern cetaceans have dramatically changed from herbivorous to carnivorous. However, the genetic basis underlying this dietary switch remains unexplored. Here, we present the first systematic investigation of 10 digestive enzymes genes (i.e., CYP7A1, CTRC, LIPC, LIPF, PNLIP, PGC, PRSS1, SI, SLC5A1, and TMPRSS15) of representative cetaceans, and the evolutionary trajectory of RNASE1 in cetartiodactylans. Positive selections were detected with proteinases (i.e., CTRC, PRSS1, and TMPRSS15) and lipases (i.e., CYP7A1, LIPF, and PNLIP) suggesting that cetaceans have evolved an enhanced digestion capacity for proteins and lipids, the major nutritional components of their prey (fishes and invertebrates). In addition, it was found that RNASE1 gene duplicated after the cetartiodactylan speciation and two independent gene duplication events took place in Camelidae and Ruminantia. Positive selection was detected with RNASE1 of Camelidae and Bovidae, suggesting enhanced digestive efficiency in the ruminants. Remarkably, even though the ancestors of cetaceans were terrestrial artiodactyls that are herbivorous, modern cetaceans lost the pancreatic RNASE1 copy with digestive function, which is in accordance with the dietary change from herbivorous to carnivorous. In sum, this is the first study that provides new insights into the evolutionary mechanism of dietary switch in cetaceans.

Bovine brain ribonuclease is the functional homolog of human ribonuclease 1

Mounting evidence suggests that human pancreatic ribonuclease (RNase 1) plays important roles in vivo, ranging from regulating blood clotting and inflammation to directly counteracting tumorigenic cells. Understanding these putative roles has been pursued with continual comparisons of human RNase 1 to bovine RNase A, an enzyme that appears to function primarily in the ruminant gut. Our results imply a different physiology for human RNase 1. We demonstrate distinct functional differences between human RNase 1 and bovine RNase A. Moreover, we characterize another RNase 1 homolog, bovine brain ribonuclease, and find pronounced similarities between that enzyme and human RNase 1. We report that human RNase 1 and bovine brain ribonuclease share high catalytic activity against double-stranded RNA substrates, a rare quality among ribonucleases. Both human RNase 1 and bovine brain RNase are readily endocytosed by mammalian cells, aided by tight interactions with cell surface glycans. Finally, we show that both human RNase 1 and bovine brain RNase are secreted from endothelial cells in a regulated manner, implying a potential role in vascular homeostasis. Our results suggest that brain ribonuclease, not RNase A, is the true bovine homolog of human RNase 1, and provide fundamental insight into the ancestral roles and functional adaptations of RNase 1 in mammals.

Ancient expansion of the ribonuclease A superfamily revealed by genomic analysis of placental and marsupial mammals

Members of the ribonuclease (RNase) A superfamily participate in a diverse array of biological processes, including digestion, angiogenesis, innate immunity, and possibly male reproduction. The superfamily is vertebrate-specific, with 13-20 highly divergent members in primates and rodents, but only a few members in chicken and fish. This has led to the proposal that the superfamily started off from a progenitor with structural similarities to angiogenin and that the superfamily underwent a dramatic expansion during mammalian evolution. To date this evolutionary expansion and understand the functional diversification of the superfamily, we here determine its entire repertoire in the sequenced genomes of dog, cow, and opossum. We identified 7, 20, and 21 putatively functional RNase genes from these three species, respectively. Many of the identified genes are highly divergent from all previously known RNase genes, thus representing new lineages within the superfamily. Phylogenetic analysis indicates that the superfamily expansion predated the separation of placental and marsupial mammals and that differential gene loss and duplication occurred in different species, generating a great variation in gene number and content among extant mammals.

Pseudogenes: are they "junk" or functional DNA?

Pseudogenes have been defined as nonfunctional sequences of genomic DNA originally derived from functional genes. It is therefore assumed that all pseudogene mutations are selectively neutral and have equal probability to become fixed in the population. Rather, pseudogenes that have been suitably investigated often exhibit functional roles, such as gene expression, gene regulation, generation of genetic (antibody, antigenic, and other) diversity. Pseudogenes are involved in gene conversion or recombination with functional genes. Pseudogenes exhibit evolutionary conservation of gene sequence, reduced nucleotide variability, excess synonymous over nonsynonymous nucleotide polymorphism, and other features that are expected in genes or DNA sequences that have functional roles. We first review the Drosophila literature and then extend the discussion to the various functional features identified in the pseudogenes of other organisms. A pseudogene that has arisen by duplication or retroposition may, at first, not be subject to natural selection if the source gene remains functional. Mutant alleles that incorporate new functions may, nevertheless, be favored by natural selection and will have enhanced probability of becoming fixed in the population. We agree with the proposal that pseudogenes be considered as potogenes, i.e., DNA sequences with a potentiality for becoming new genes.

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.

Secretory ribonucleases in the primitive ruminant chevrotain (Tragulus javanicus)

Phylogenetic analyses of secretory ribonucleases or RNases 1 have shown that gene duplication events, giving rise to three paralogous genes (pancreatic, seminal and brain RNase), occurred during the evolution of ancestral ruminants. A higher number of paralogous sequences are present in chevrotain (Tragulus javanicus), the earliest diverged taxon within the ruminants. Two pancreatic RNase sequences were identified, one encoding the pancreatic enzyme, the other encoding a pseudogene. The identity of the pancreatic enzyme was confirmed by isolation of the protein and N-terminal sequence analysis. It is the most acidic pancreatic ribonuclease identified so far. Formation of the mature enzyme requires cleavage by signal peptidase of a peptide bond between two glutamic acid residues. The seminal-type RNase gene shows features of a pseudogene, like orthologous genes in other ruminants investigated with the exception of the bovine species. The brain-type RNase gene of chevrotain is expressed in brain tissue. A hybrid gene with a pancreatic-type N-terminal and a brain-type C-terminal sequence has been identified but nothing is known about its expression. Phylogenetic analysis of RNase 1 sequences of six ruminant, three other artiodactyl and two whale species support previous findings that two gene duplications occurred in a ruminant ancestor. Three distinct groups of pancreatic, seminal-type and brain-type RNases have been identified and within each group the chevrotain sequence it the first to diverge. In taxa with duplications of the RNase gene (ruminants and camels) the gene evolved at twice as fast than in taxa in which only one gene could be demonstrated; in ruminants there was an approximately fourfold increase directly after the duplications and then a slowing in evolutionary rate.

Ruminant brain ribonucleases: expression and evolution

Molecular evolutionary analyses of mammalian ribonucleases have shown that gene duplication events giving rise to three paralogous genes occurred in ruminant ancestors. One of these genes encodes a ribonuclease identified in bovine brain. A peculiar feature of this enzyme and orthologous sequences in other ruminants are C-terminal extensions consisting of 17-27 amino acid residues. Evidence was obtained by Western blot analysis for the presence of brain-type ribonucleases in brain tissue not only of ox, but also of sheep, roe deer and chevrotain (Tragulus javanicus), a member of the earliest diverged taxon of the ruminants. The C-terminal extension of brain-type ribonuclease from giraffe deviates much in sequence from orthologues in other ruminants, due to a change of reading frame. However, the gene encodes a functional enzyme, which could be expressed in heterologous systems. The messenger RNA of bovine brain ribonuclease is not only expressed at a high level in brain tissue but also in lactating mammary gland. The enzyme was isolated and identified from this latter tissue, but was not present in bovine milk, although pancreatic ribonucleases A and B could be isolated from both sources. This suggests different ways of secretion of the two enzyme types, possibly related to structural differences. The sequence of the brain-type RNase from chevrotain suggests that the C-terminal extensions of ruminant brain-type ribonucleases originate from deletions in the ancestral DNA (including a region with stop codons), followed by insertion of a 5-8-fold repeated hexanucleotide sequence, coding for a proline-rich polypeptide.

The melanocyte-stimulating hormone receptor (MC1-R) gene as a tool in evolutionary studies of artiodactyles

The complete coding region of the melanocyte-stimulating hormone receptor (MC1-R) gene was characterized in species belonging to the two families Bovidae and Cervidae; cattle (Bos taurus), sheep (Ovis aries), goat (Capra hircus), muskox (Ovibos moschatus), roe deer (Capreolus capreolus), reindeer (Rangifer tarandus), moose (Alces alces), red deer (Cervus elaphus) and fallow deer (Dama dama). This well conserved gene is a central regulator of mammalian coat colour. Examination of the interspecies variability revealed a 5.3-6.8% divergence between the Cervidae and Bovidae families, whereas the divergence within the families were 1.0-3.1% and 1.2-4.6%, respectively. Complete identity was found when two subspecies of reindeer, Eurasian tundra reindeer (R.t. tarandus) and Svalbard reindeer (R.t. platvrhynehus), were analyzed. An rooted phylogenetic tree based on Bovidae and Cervidae MC1-R DNA sequences was in complete agreement with current taxonomy, and was supported by bootstrapping analysis. Due to different frequencies of silent vs. replacement mutations, the amino acid based phylogenetic tree contains several dissimilarities when compared to the DNA based phylogenetic tree.

Evolution of oligomeric proteins. The unusual case of a dimeric ribonuclease

The model system made up of a monomeric and a dimeric ribonuclease of the pancreatic-type superfamily has recently attracted the attention of investigators interested in the evolution of oligomeric proteins. In this system, bovine pancreatic ribonuclease (RNase A) is the monomeric prototype, and bovine seminal ribonuclease (BS-RNase) is the dimeric counterpart. However, this evolutionary case is unusual, as BS-RNase is the only dimeric member of the whole large superfamily comprising more than 100 identified members from amphibia, aves, reptilia and mammalia. Furthermore, although the seminal-type RNase gene can be traced back to the divergence of the ruminants, it is expressed only in a single species (Bos taurus). These unusual findings are discussed, as well as previous hypotheses on the evolution of seminal RNase. Furthermore, a new 'minimalist' hypothesis is proposed, in line with basic principles of structural biology and molecular evolution.

Seminal-type ribonuclease genes in ruminants, sequence conservation without protein expression?

Bovine seminal ribonuclease (BS-RNase) is an interesting enzyme both for functional and structural reasons. The enzyme is the product of a gene duplication that occurred in an ancestral ruminant. It is possible to demonstrate the presence of seminal-type genes in all other investigated ruminant species, but they are not expressed and show features of pseudogenes. In this paper we report the determination of two pancreatic and one seminal-type ribonuclease gene sequences of swamp-type water buffalo (Bubalus bubalis). The two pancreatic sequences encode proteins with identical amino acid sequences as previously determined for the enzymes isolated from swamp-type and river-type water buffalo, respectively. The seminal-type sequence has no pseudogene features and codes for an enzyme with no unusual features compared with the active bovine enzyme, except for the replacement of one of the cysteines which takes part in the two intersubunit disulfide bridges. However, Western blotting demonstrates the presence of only small amounts of the pancreatic enzymes in water buffalo semen, suggesting that also in this species the seminal-type sequence is not expressed. But it is still possible that the gene is expressed somewhere else in the body or during development. Reconstruction of seminal-type ribonuclease sequences in ancestors of Bovinae and Bovidae indicates no serious abnormalities in the encoded proteins and leads us to the hypothesis that the ruminant seminal-type ribonuclease gene has not come to expression during most of its evolutionary history, but did not exhibit a high evolutionary rate that is generally observed in pseudogenes.

The differential pattern of tissue-specific expression of ruminant pancreatic type ribonucleases may help to understand the evolutionary history of their genes

Molecular evolutionary analyses of mammalian ribonucleases have shown that gene duplication events giving three paralogous genes occurred in ruminant ancestors. The enzymes of the bovine species encoded by these genes, isolated from pancreas, brain and seminal vesicles, present similar enzymological properties but distinct structural features. In other ruminant species, genomic sequences orthologous to the bovine genes of pancreas and brain ribonucleases encode active enzymes. In mammalian species other than ruminant artiodactyls, only one gene encoding ribonuclease of the pancreatic type is generally present. In this work, we describe a differential pattern of transcriptional expression of the pancreas and brain ribonuclease genes in the ox species and report transcription of the human ribonuclease gene in brain as well as in pancreas and in mammary gland. We also report the molecular cloning of the gene encoding the bovine seminal ribonuclease in which the structural organization already described for the two paralogous genes is conserved. The seminal RNAase is exclusively expressed in seminal vesicles of Bos taurus, whereas in other ruminant species, the orthologous sequence is a pseudogene. Previous studies from a number of research groups demonstrated that, unlike other mammalian ribonucleases, the seminal enzyme is a covalent dimer, and its unique quaternary structure correlates with special biological activities. The major determinant of dimer formation, i.e. the presence of two adjacent cysteine residues, is absent in the pseudogenes. We advance the hypothesis that the differentiation of distinct expression patterns could represent an important evolutionary determinant for the genes encoding pancreas and brain ribonucleases in ruminants, whereas the differentiation of a quaternary structure endowed with new biological functions could be the main determinant for the evolutionary success of the seminal gene in the bovine species.

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.

Ribonuclease k6: chromosomal mapping and divergent rates of evolution within the RNase A gene superfamily

We have localized the gene encoding human RNase k6 to within approximately 120 kb on the long (q) arm of chromosome 14 by HAPPY mapping. With this information, the relative positions of the six human RNase A ribonucleases that have been mapped to this locus can be inferred. To further our understanding of the individual lineages comprising the RNase A superfamily, we have isolated and characterized 10 novel genes orthologous to that encoding human RNase k6 from Great Ape, Old World, and New World monkey genomes. Each gene encodes a complete ORF with no less than 86% amino acid sequence identity to human RNase k6 with the eight cysteines and catalytic histidines (H15 and H123) and lysine (K38) typically observed among members of the RNase A superfamily. Interesting trends include an unusually low number of synonymous substitutions (Ks) observed among the New World monkey RNase k6 genes. When considering nonsilent mutations, RNase k6 is a relatively stable lineage, with a nonsynonymous substitution rate of 0.40 x 10(-9) nonsynonymous substitutions/nonsynonymous site/year (ns/ns/yr). These results stand in contrast to those determined for the primate orthologs of the two closely related ribonucleases, the eosinophil-derived neurotoxin (EDN) and eosinophil cationic protein (ECP), which have incorporated nonsilent mutations at very rapid rates (1.9 x 10(-9) and 2.0 x 10(-9) ns/ns/yr, respectively). The uneventful trends observed for RNase k6 serve to spotlight the unique nature of EDN and ECP and the unusual evolutionary constraints to which these two ribonuclease genes must be responding. [The sequence data described in this paper have been submitted to the GenBank data library under accession nos. AF037081-AF037090.]

Origin of dimeric structure in the ribonuclease superfamily

To enable application of postgenomic evolutionary approaches to understand the divergence of behavior and function in ribonucleases (RNases), the impact of divergent sequence on the divergence of tertiary and quaternary structure is analyzed in bovine pancreatic and seminal ribonucleases, which differ by 23 amino acids. In a crystal, seminal RNase is a homodimer joined by two "antiparallel" intersubunit disulfide bonds between Cys-31 from one subunit and Cys-32' from the other and having composite active sites arising from the "swap" of residues 1-20 from each subunit. Specialized Edman degradation techniques have completed the structural characterization of the dimer in solution, new cross-linking methods have been developed to assess the swap, and sequence determinants of quaternary structure have been explored by protein engineering using the reconstructed evolutionary history of the protein family as a guide. A single Cys at either position 32 (the first to be introduced during the divergent evolution of the family) or 31 converts monomeric RNase A into a dimer. Even with an additional Phe at position 31, another residue introduced early in the seminal lineage, swap is minimal. A hydrophobic contact formed by Leu-28, however, also introduced early in the seminal lineage, increases the amount of "antiparallel" connectivity of the two subunits and facilitates swapping of residues 1-20. Efficient swapping requires addition of a Pro at position 19, a residue also introduced early in the divergent evolution of the seminal RNase gene. Additional cysteines required for dimer formation are found to slow refolding of the protein through formation of incorrect disulfide bonds, suggesting a paradox in the biosynthesis of the protein. Further studies showed that the dimeric form of seminal RNase known in the crystal is not the only form in vivo, where a substantial amount of heterodimer is known. These data complete the acquisition of the background needed to understand the evolution of new structure, behavior, and function in the seminal RNase family of proteins.

Molecular cloning of four novel murine ribonuclease genes: unusual expansion within the ribonuclease A gene family

We have characterized four novel murine ribonuclease genes that, together with the murine eosinophil-associated ribonucleases 1 and 2, form a distinct and unusual cluster within the RNase A gene superfamily. Three of these genes (mR-3, mR-4, mR-5) include complete open reading frames, encoding ribonucleases with eight cysteines and appropriately spaced histidines (His11 and His124) and lysine (Lys35) that are characteristic of this enlarging protein family; the fourth sequence encodes a non-functional pseudogene (mR-6P). Although the amino acid sequence similarities among these murine ribonucleases varies from 60 to 94%, they form a unique cluster, as each sequence is found to be more closely related to another of this group than to either murine angiogenin or to murine pancreatic ribonuclease. Interestingly, the relationship between the six genes in this 'mR cluster' and the defined lineages of the RNase A gene family could not be determined by amino acid sequence homology, suggesting the possibility that there are one or more additional ribonuclease lineages that have yet to be defined. Although the nature of the evolutionary constraints promoting this unusual expansion and diversification remain unclear, the implications with respect to function are intriguing.

An angiogenic protein from bovine serum and milk--purification and primary structure of angiogenin-2

Bovine serum and milk contain a basic angiogenic protein that binds tightly to placental ribonuclease inhibitor. It was purified from both sources by ion-exchange and reversed-phase chromatographies. Its amino acid sequence revealed that it is a member of the ribonuclease superfamily. It contains 123 amino acids in a single polypeptide chain, is cross-linked by three disulfide bonds, is glycosylated at Asn33, and is 57% identical to bovine angiogenin. The amino-terminal and carboxyl-terminal residues are pyroglutamic acid and proline, respectively. The protein has ribonucleolytic activity that is similar to, but somewhat lower than, that of bovine angiogenin, i.e. very low relative to RNase. It is angiogenically potent on chicken chorioallantoic membrane, but less so than angiogenin. The sequence and activities demonstrate that this protein is a second, distinct, member of the angiogenin sub-family of pancreatic ribonucleases, and is referred to as angiogenin-2.

From housekeeper to microsurgeon: the diagnostic and therapeutic potential of ribonucleases

The RNA population in cells is controlled post-transcriptionally by ribonucleases (RNases) of varying specificity. Angiogenin, neurotoxins, and plant allergens are among many proteins with RNase activity or significant homology to known RNases. RNase activity in serum and cell extracts is elevated in a variety of cancers and infectious diseases. RNases are regulated by specific activators and inhibitors, including interferons. Many of these regulatory molecules are useful lead compounds for the design of drugs to control tumor angiogenesis, allergic reactions, and viral replication. One RNase (Onconase) and several RNase activators are now in clinical trials for cancer treatment or inhibition of chronic virus infections. Several others, alone or conjugated with specific cell binding molecules, are being developed for their antifungal, antiviral, and antitumor cell activity.

Molecular cloning and characterization of a novel human ribonuclease (RNase k6): increasing diversity in the enlarging ribonuclease gene family

The discovery of Ribonuclease k6 (RNase k6) was an unexpected result of our ongoing efforts to trace the evolutionary history of the ribonuclease gene family. The open reading frame of RNase k6, amplified from human genomic DNA, encodes a 150 amino acid polypeptide with eight cysteines and histidine and lysine residues corresponding to those found in the active site of the prototype, ribonuclease A. The single-copy gene encoding RNase k6 maps to human chromosome 14 and orthologous sequences were detected in both primate and non-primate mammalian species. A single mRNA transcript (1.5 kb) was detected in all human tissues tested, with lung representing the most abundant source. At the cellular level, transcripts encoding RNase k6 were detected in normal human monocytes and neutrophils (but not in eosinophils) suggesting a role for this ribonuclease in host defense. Of the five previously identified human ribonucleases of this group, RNase k6 is most closely related to eosinophil-derived neurotoxin (EDN), with 47% amino acid sequence identity; slight cross-reactivity between RNase k6 and EDN was observed on Western blots probed with polyclonal anti-EDN antiserum. The catalytic constants determined, Km = 5.0 microM and Kcat = 0.13 s-1, indicate that recombinant RNase k6 has approximately 40-fold less ribonuclease activity than recombinant EDN. The identification and characterization of RNase k6 has extended the ribonuclease gene family and suggests the possibility that there are others awaiting discovery.