Conserved restriction sites within the ribosomal RNA genes of vertebrates

We have mapped the cleavage sites of four restriction enzymes which recognize six-base sequences within the nuclear ribosomal (rRNA) genes of twelve vertebrates, including several placental mammals (Homo sapiens, man; Bos taurus, cow; Equus caballus, horse; Sus scofra, pig; Ovis aries, sheep; Rattus rattus, rat), a marsupial (Didelphis marsupialis, opossum), a bird (Gallus domesticus, chicken), an amphibian (Xenopus laevis), a reptile (Alligator mississipiensis), a bony fish (Cynoscion nebulosus, sea trout), and a cartilagenous fish (Carcharhinus species, requiem shark). These animals represent a span of approx. 400 million years of evolutionary divergence. Our data identify restriction sites in the rRNA genes which are highly conserved among higher vertebrates and therefore are likely to be in functionally important regions. Additionally, the restriction enzyme sites identified will be useful in cloning and sequencing the rRNA genes in any vertebrate. Finally, the consistent size and conserved sequence homology suggests that these rRNA gene segments will be useful as internal controls in hybridization experiments involving other genomic regions in vertebrates.     

Erythrocyte acid phosphatase: species specificity in activity modulation by purine analogs

The previous studies of the interaction of purine analogs and human erythrocyte acid phosphatase isozymes were extended to include erythrocyte acid phosphatase from seven other species. Consistent responses, similar to the observations with the several genotypically different human isozymes, were observed. The isozyme from chimpanzee erythrocytes was similar to the human B-type isozyme while the baboon and cow isozymes were at the other extreme in responsiveness and were more divergent from the B-isozyme than was the human A-type isozyme. The ACP from rabbit, dog, sheep and rhesus erythrocytes exhibited intermediate levels of responsiveness but did differ from the human A-type isozyme. Additional studies indicated some differences between the responsiveness of the partially purified erythrocyte enzyme and the low molecular weight ACP from liver.     

Identification and molecular evolution of cow CENP-A gene family

The centromere protein A (CENP-A), a histone H3-like protein, provides an essential role for chromosomal segregation during mitosis and meiosis. In this study we identified ten new CENP-A-like genes (excluding the original CENP-A gene) in cow by searching its genome database, while all other examined mammals contained only a single copy of the CENP-A gene. Phylogenetic analysis shows that the expansion of these genes occurred before the divergence of cow and sheep but after the artiodactyls diverged from other mammals. Our analyses also indicate that multiple gene duplication and intron loss events have occurred during the formation of this gene family, and the high similarity in intron sequences between the new genes and the old one suggests very recent gene expansion events. Furthermore, evidence from the comparisons of the synonymous and nonsynonymous substitutions and the distributions of the mutation sites suggested that the CENP-A-L-1 gene, a potentially functional member of cow CENP-A gene family, evolved under positive selection. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Estradiol stabilizes estrogen receptor messenger ribonucleic acid in sheep endometrium via discrete sequence elements in its 3'-untranslated region

The preovulatory surge of estrogen up-regulates estrogen receptor-alpha (ER) gene expression in the uterus during the estrous/menstrual cycles of female mammals. Previously, we demonstrated that the 5-fold increase in ER mRNA levels in endometrium of ovariectomized ewes treated with a physiological dose of estradiol (E2) is entirely due to an increase in ER mRNA stability. Our current work confirms that the E2 effect is specific to ER mRNA. The sequence of ER mRNA, cloned from sheep endometrium, shows a high degree of conservation with those of other species, even in the 5'- and the very long 3'-untranslated regions. In a cell-free assay, ER mRNA demonstrates greater stability with endometrial extracts from E2-treated ewes compared with those from untreated ovariectomized ewes. The E2-enhanced stability of ER mRNA was ablated by prior treatment of the extracts with proteinase K, 70 C heat, and oxidizing and alkylating reagents, indicating that a protein is responsible for stabilization of the message. The 3'-untranslated region of ER mRNA contains discrete sequences required for E2-enhanced stability, four of which were identified by extensive deletion mutant analyses. Transfer of two of the four minimal E2-modulated stability sequences conferred E2-enhanced stability to a heterologous RNA. These minimal E2-modulated stability sequences contain a common 10-base, uridine-rich sequence that is predicted to reside in a loop structure. Throughout our studies, estrogen stabilization of ER mRNA in sheep endometrium resembled that of vitellogenin mRNA in frog liver, indicating conservation of this ancient mechanism for enhancing gene expression in response to estrogen.     

The MB1 family of repeats in clones from the genomes of mammals]

Members of MB1 family repeats are revealed in genomes of many mammals (cow, rabbit, opossum, horse, ...). The MB1 repeats from cow and rabbit genomes are mirror-reflected about the SINE families repeats from cow and rabbit genomes. The life time of MB1 repeats are no less than 100 million years. Classification of MB1 repeats from human genome using the information similarity was performed. This classification has revealed two subfamily MB1 repeats in human genome. Possible processes of creation of MB1 family repeats common for many mammals are discussed.     

Selection on amino acid substitutions in Arabidopsis

Studies of nucleotide diversity have found an excess of low-frequency amino acid polymorphisms segregating in Arabidopsis thaliana, suggesting a predominance of weak purifying selection acting on amino acid polymorphism in this inbreeding species. Here, we investigate levels of diversity and divergence at synonymous and nonsynonymous sites in 6 circumpolar populations of the outbreeding Arabidopsis lyrata and compare these results with A. thaliana, to test for differences in mutation and selection parameters across genes, populations, and species. We find that A. lyrata shows an excess of low-frequency nonsynonymous polymorphisms both within populations and species wide, consistent with weak purifying selection similar to the patterns observed in A. thaliana. Furthermore, nonsynonymous polymorphisms tend to be more restricted in their population distribution in A. lyrata, consistent with purifying selection preventing their geographic spread. Highly expressed genes show a reduced ratio of amino acid to synonymous change for both polymorphism and fixed differences, suggesting a general pattern of stronger purifying selection on high-expression proteins.     

Pattern of nucleotide substitution and divergence of prophenoloxidase in decapods

Despite the unprecedented development in identification and characterization of prophenoloxidase (proPO) in commercially important decapods, little is known about the evolutionary relationship, rate of amino acid replacement and differential selection pressures operating on proPO of different species of decapods. Here we report the evolutionary relationship among these nine decapod species based on proPO gene and types of selective pressures operating on proPO codon sites. Our analyses revealed that all the nine decapod species shared a common ancestor. The mean percentage sequence divergence at proPO gene was 34.4+/-0.6%. Pairwise estimates of nonsynonymous to synonymous ratio (omega) for Homarus americanus-H. gammarus is greater than one, therefore indicating adaptive evolution (functional diversification) of proPO in these two species. In contrast, strong purifying selection (omega<1) was observed in all other species pairs. However, phylogenetically closely related decapods revealed relatively higher omega value (omega=0.15+/-0.3) than the distantly related species pairs (omega=0.0075+/-0.005). These discrepancies could be due to higher fixation probability of beneficial mutation in closely related species. Maximum likelihood-based codon substitution analyses revealed a strong purifying selection operating on most of the codon sites, therefore suggesting proPO is functionally constrained (purifying selection). Codon substitution analyses have also revealed the evidence of strong purifying selection in haemocyanin subunits of decapods.     

Structure of the goat psi beta y beta-globin pseudogene. Analysis of goat pseudogene evolutionary patterns

The 12-member beta-globin gene locus of the goat contains three beta(adult)-type pseudogenes, one in each of three four-gene subsets of the locus. We have determined the complete nucleotide sequence of psi beta y, the pseudogene present in the most downstream four-gene subset, which also contains the functional fetal gene, beta F. psi beta y contains, throughout its length, numerous incapacitating mutations in common with the previously sequenced goat psi beta x and psi beta z pseudogenes consistent with the model that all were descended from a common pseudogene ancestor which became defective prior to the expansion of the beta-globin locus in the goat lineage. Evolutionary analysis of the psi beta y sequence in comparison to psi beta x and psi beta z provides evidence that nucleotide substitutions were fixed in a random manner within these pseudogenes with respect to polarity, coding versus non-coding regions, and replacement sites versus silent sites. However, substitutions appear to have accumulated asymmetrically between different pseudogenes in a manner that provides evidence for partial gene conversion. Moreover, the presence of deletions in goat psi beta y, which are also observed in the cow pseudogene psi 2, but not in the cow psi 1 pseudogene, indicate that goat psi beta y and cow psi 2 are orthologous but cow psi 1 actually arose prior to the goat/cow divergence. The authentic goat orthologue to cow psi 1 temporarily existed in the goat lineage but was deleted, probably prior to the divergence of goats and sheep.     

Testing the neutral theory of molecular evolution using genomic data: a comparison of the human and bovine transcriptome

Despite growing evidence of rapid evolution in protein coding genes, the contribution of positive selection to intra- and interspecific differences in protein coding regions of the genome is unclear. We attempted to see if genes coding for secreted proteins and genes with narrow expression, specifically those preferentially expressed in the mammary gland, have diverged at a faster rate between domestic cattle (Bos taurus) and humans (Homo sapiens) than other genes and whether positive selection is responsible. Using a large data set, we identified groups of genes based on secretion and expression patterns and compared them for the rate of nonsynonymous (dN) and synonymous (dS) substitutions per site and the number of radical (Dr) and conservative (Dc) amino acid substitutions. We found evidence of rapid evolution in genes with narrow expression, especially for those expressed in the liver and mammary gland and for genes coding for secreted proteins. We compared common human polymorphism data with human-cattle divergence and found that genes with high evolutionary rates in human-cattle divergence also had a large number of common human polymorphisms. This argues against positive selection causing rapid divergence in these groups of genes. In most cases dN/dS ratios were lower in human-cattle divergence than in common human polymorphism presumably due to differences in the effectiveness of purifying selection between long-term divergence and short-term polymorphism.     

Recombination and selection in the evolution of picornaviruses and other Mammalian positive-stranded RNA viruses

Picornaviridae are a large virus family causing widespread, often pathogenic infections in humans and other mammals. Picornaviruses are genetically and antigenically highly diverse, with evidence for complex evolutionary histories in which recombination plays a major part. To investigate the nature of recombination and selection processes underlying the evolution of serotypes within different picornavirus genera, large-scale analysis of recombination frequencies and sites, segregation by serotype within each genus, and sequence selection and composition was performed, and results were compared with those for other nonenveloped positive-stranded viruses (astroviruses and human noroviruses) and with flavivirus and alphavirus control groups. Enteroviruses, aphthoviruses, and teschoviruses showed phylogenetic segregation by serotype only in the structural region; lack of segregation elsewhere was attributable to extensive interserotype recombination. Nonsegregating viruses also showed several characteristic sequence divergence and composition differences between genome regions that were absent from segregating virus control groups, such as much greater amino acid sequence divergence in the structural region, markedly elevated ratios of nonsynonymous-to-synonymous substitutions, and differences in codon usage. These properties were shared with other picornavirus genera, such as the parechoviruses and erboviruses. The nonenveloped astroviruses and noroviruses similarly showed high frequencies of recombination, evidence for positive selection, and differential codon use in the capsid region, implying similar underlying evolutionary mechanisms and pressures driving serotype differentiation. This process was distinct from more-recent sequence evolution generating diversity within picornavirus serotypes, in which neutral or purifying selection was prominent. Overall, this study identifies common themes in the diversification process generating picornavirus serotypes that contribute to understanding of their evolution and pathogenicity.     

Discordant evolution of the adjacent antiretroviral genes TRIM22 and TRIM5 in mammals

TRIM5alpha provides a cytoplasmic block to retroviral infection, and orthologs encoded by some primates are active against HIV. Here, we present an evolutionary comparison of the TRIM5 gene to its closest human paralogs: TRIM22, TRIM34, and TRIM6. We show that TRIM5 and TRIM22 have a dynamic history of gene expansion and loss during the evolution of mammals. The cow genome contains an expanded cluster of TRIM5 genes and no TRIM22 gene, while the dog genome encodes TRIM22 but has lost TRIM5. In contrast, TRIM6 and TRIM34 have been strictly preserved as single gene orthologs in human, dog, and cow. A more focused analysis of primates reveals that, while TRIM6 and TRIM34 have evolved under purifying selection, TRIM22 has evolved under positive selection as was previously observed for TRIM5. Based on TRIM22 sequences obtained from 27 primate genomes, we find that the positive selection of TRIM22 has occurred episodically for approximately 23 million years, perhaps reflecting the changing pathogenic landscape. However, we find that the evolutionary episodes of positive selection that have acted on TRIM5 and TRIM22 are mutually exclusive, with generally only one of these genes being positively selected in any given primate lineage. We interpret this to mean that the positive selection of one gene has constrained the adaptive flexibility of its neighbor, probably due to genetic linkage. Finally, we find a striking congruence in the positions of amino acid residues found to be under positive selection in both TRIM5alpha and TRIM22, which in both proteins fall predominantly in the beta2-beta3 surface loop of the B30.2 domain. Astonishingly, this same loop is under positive selection in the multiple cow TRIM5 genes as well, indicating that this small structural loop may be a viral recognition motif spanning a hundred million years of mammalian evolution.     

Transposable elements as a source of genetic innovation: expression and evolution of a family of retrotransposon-derived neogenes in mammals

A family of functional neogenes called Mart, related to the gag gene of Sushi-like long terminal repeat retrotransposons from fish and amphibians, is present in the genome of human (11 genes) and other primates, as well as in mouse (11 genes), rat, dog (12 genes), cat, and cow. Mart genes have lost their capacity of retrotransposition through non-functionalizing rearrangements having principally affected long terminal repeats and pol open reading frame. Most Mart genes are located on the X chromosome in different mammals. Sequence database analysis suggested that Mart genes are present in opossum (marsupial), but absent from the genome of chicken. Hence, the Mart gene family might have been formed from Sushi-like retrotransposon(s) after the split of birds and mammals (310 myr ago), but before the divergence between placental mammals and marsupials (170 myr ago). RT-PCR analysis showed that at least six Mart genes are expressed during mouse embryonic development, with in situ hybridization analysis revealing rather ubiquitous expression patterns. Mart expression was also detected in adult mice, with some genes being expressed in all tissues tested, while others showed a much more restricted expression pattern. Although additional analysis will be required to establish the function of the retrotransposon-derived Mart neogenes, these observations support the evolutionary importance of retrotransposable elements as a source of genetic novelty.     

Ovine placental lactogen specifically binds to endometrial glands of the ovine uterus

A hormonal servomechanism has been proposed to regulate differentiation and function of the endometrial glandular epithelium (GE) in the ovine uterus during pregnancy. This mechanism involves sequential actions of estrogen, progesterone, ovine interferon tau (IFNtau), placental lactogen (oPL), and placental growth hormone (oGH). The biological actions of oPL in vitro are mediated by homodimerization of the prolactin receptor (oPRLR) and heterodimerization of the oPRLR and oGH receptor. The objectives of the study were to determine the effects of intrauterine oPL, oGH, and their combination on endometrial histoarchitecture and gene expression and to localize and characterize binding sites for oPL in the ovine uterus in vivo using an in situ ligand binding assay. Intrauterine infusion of oPL and/or oGH following IFNtau into ovariectomized ewes treated with progesterone daily differentially affected endometrial gland number and expression of uterine milk proteins and osteopontin. However, neither hormone affected PRLR, insulin-like growth factor (IGF)-I, or IGF-II mRNA levels in the endometrium. A chimeric protein of placental secretory alkaline phosphatase (SEAP) and oPL was used to identify and characterize binding sites for oPL in frozen sections of interplacentomal endometrium from pregnant ewes. Specific binding of SEAP-oPL was detected in the endometrial GE on Days 30, 60, 90, and 120 of pregnancy. In Day 90 endometrium, SEAP-oPL binding to the endometrial GE was displaced completely by oPL and prolactin (oPRL) but only partially by oGH. Binding experiments using the extracellular domain of the oPRLR also showed that iodinated oPL binding sites could be competed for by oPRL and oPL but not by oGH. Collectively, results indicate that oPL binds to receptors in the endometrial glands and that oPRL is more effective than oGH in competing for these binding sites. Thus, effects of oPL on the endometrial glands may be mediated by receptors for oPRL and oGH.     

Immunocytochemical localization and messenger ribonucleic acid levels of a progesterone-dependent endometrial secretory protein (cathepsin L) in the pregnant cat uterus.

The purpose of this study was to localize immunocytochemically a progesterone-dependent protein (PDP) and to determine PDP mRNA levels during the initial stage of the implantation period. Uterine tissue was collected from Day 0-18 postcoital animals. The tissue was processed for immunocytochemical localization of PDP, and the endometrial RNA was isolated and analyzed for PDP gene expression by slot-blot hybridization. PDP was detected immunocytochemically as early as Day 5 postcoitus in the epithelial cells of the deep uterine glands, and the intensity of immunostaining appeared to peak by Day 12 postcoitus. PDP was absent in the endometrium obtained from implantation sites after Day 16 postcoitus, but the synthesis of PDP was maintained in the endometrium obtained from nonimplantation sites. Immunogold electron microscopy demonstrated that PDP was present in electron-dense granules of the glandular epithelial cells. PDP mRNA was detectable in the endometrium at Day 5 postcoitus and peaked around Day 10 postcoitus. PDP mRNA was absent in the endometrium from implantation sites after Day 16 postcoitus, but was maintained in the endometrium from nonimplantation sites. In summary, the results of this study illustrate that PDP is synthesized within the epithelial cells of the deep uterine glands, packaged within membrane-bound secretory granules, and released into the uterine lumen. Also, the process of implantation alters the gene expression in a very localized way since PDP mRNA and PDP-positive granules were absent in the endometrial glands obtained from the implantation site within 1-2 days of the onset of implantation, whereas both PDP mRNA and PDP-positive granules were maintained in the endometrial glands from nonimplantation-site regions.(ABSTRACT TRUNCATED AT 250 WORDS)     

Diversifying and Purifying Selection in the Peptide Binding Region of <Emphasis Type="Italic">DRB</Emphasis> in Mammals

The class II genes of the major histocompatibility complex encode proteins which play a crucial role in antigen presentation. They are among the most polymorphic proteins known, and this polymorphism is thought to be the result of natural selection. To understand the selective pressure acting on the protein and to examine possible differences in the evolutionary dynamics among species, we apply maximum likelihood models of codon substitution to analyze the DRB genes of six mammalian species: human, chimpanzee, macaque, tamarin, dog, and cow. The models account for variable selective pressures across codons in the gene and have the power to detect amino acid residues under either positive or negative selection. Our analysis detected positive selection in the DRB genes in each of the six mammals examined. Comparison with structural data reveals that almost all amino acid residues inferred to be under positive selection in humans are in the peptide binding region (PBR) and are in contact with the antigen side chains, although residues outside of but close to the PBR are also detected. Strong purifying selection is also detected in the PBR, at sites which contact the antigen and at sites which may be involved in dimerization or T cell binding. The analysis demonstrates the utility of the random-sites analysis even when structural information is available. The different mammalian species are found to share many positively or negatively selected sites, suggesting that their functional roles have remained very similar in the different species, despite the different habitats and pathogens of the species.     

Endogenous inhibition of arachidonic acid metabolism in the endometrium of the sheep

Arachidonic acid is metabolised via the cyclo-oxygenase pathway to several biologically active metabolites. These metabolites control important reproductive functions like luteolysis of the corpus luteum. The metabolism of arachidonic acid was studied by the enzymatic conversion of [1-14C]-labelled arachidonic acid in sheep endometrial tissue. The inhibitory capacity of sheep endometrial tissue was measured by the enzymatic conversion of [1-14C]-arachidonic acid by sheep seminal vesicular gland microsomes. Endometrial microsomes converted arachidonic acid into different prostaglandins and monohydroxy acids but at a low rate. A factor(s) inhibiting both prostaglandin and monohydroxy acid synthesis was found in both the microsomal and cytosolic fractions of endometrial tissue. A very high inhibitory potency of prostaglandin and monohydroxy acid synthesis, calculated as IC50 values, was found in cytosolic fractions. For comparison IC50 values of indomethacin, mefenamic acid, carprofen and acetylsalicylic acid were also calculated in this in vitro system. These data indicate that both prostaglandin and monohydroxy acid synthesizing capacities and an inhibitory factor(s) are present in sheep endometrium and possibly regulate arachidonic acid metabolism in this tissue.     

Molecular phylogeny,long-term evolution,and functional divergence of flavin-containing monooxygenases

Flavin-containing monooxygenases (FMOs) metabolize xenobiotic compounds, many of which are clinically important, as well as endogenous substrates as part of a discrete physiological process. The FMO gene family is conserved and ancient with representatives present in all phyla so far examined. However, there is a lack of information regarding the long-term evolution and functional divergence of these proteins. This study represents the first attempt to characterize the long-term evolution followed by the members of this family. Our analysis shows that there is extensive silent divergence at the nucleotide level suggesting that this family has been subject to strong purifying selection at the protein level. Invertebrate FMOs have a polyphyletic origin. The functional divergence of FMOs 1–5 started before the split between amphibians and mammals. The vertebrate FMO5 is more ancestral than other four FMOs. Moreover, the existence of higher levels of codon bias was detected at the N-terminal ends, which can be ascribed to the critical role played by the FAD binding motif in this region. Finally, critical amino acid residues for FMO functional divergence (type I & II) after gene duplication were detected and characterized.     

The nonsynonymous/synonymous substitution rate ratio versus the radical/conservative replacement rate ratio in the evolution of mammalian genes

There are 2 ways to infer selection pressures in the evolution of protein-coding genes, the nonsynonymous and synonymous substitution rate ratio (K(A)/K(S)) and the radical and conservative amino acid replacement rate ratio (K(R)/K(C)). Because the K(R)/K(C) ratio depends on the definition of radical and conservative changes in the classification of amino acids, we develop an amino acid classification that maximizes the correlation between K(A)/K(S) and K(R)/K(C). An analysis of 3,375 orthologous gene groups among 5 mammalian species shows that our classification gives a significantly higher correlation coefficient between the 2 ratios than those of existing classifications. However, there are many orthologous gene groups with a low K(A)/K(S) but a high K(R)/K(C) ratio. Examining the functions of these genes, we found an overrepresentation of functional categories related to development. To determine if the overrepresentation is stage specific, we examined the expression patterns of these genes at different developmental stages of the mouse. Interestingly, these genes are highly expressed in the early middle stage of development (blastocyst to amnion). It is commonly thought that developmental genes tend to be conservative in evolution, but some molecular changes in developmental stages should have contributed to morphological divergence in adult mammals. Therefore, we propose that the relaxed pressures indicated by the K(R)/K(C) ratio but not by K(A)/K(S) in the early middle stage of development may be important for the morphological divergence of mammals at the adult stage, whereas purifying selection detected by K(A)/K(S) occurs in the early middle developmental stage.