The Mammalian 2′-5′ Oligoadenylate Synthetase Gene Family: Evidence for Concerted Evolution of Paralogous Oas1 Genes in Rodentia and Artiodactyla

Multiple 2′-5′ oligoadenylate (2-5A) synthetases are important components of innate immunity in mammals. Gene families encoding these proteins have previously been studied mainly in humans and mice. To reconstruct the evolution of this gene family in mammals, a search for additional 2-5A synthetase genes was performed in rat, cattle, pig, and dog. Twelve 2′-5′ oligoadenylate synthetase (Oas) genes were identified in the rat genome, including eight Oas1 genes, two Oas1 pseudogenes, single copies of Oas2 and Oas3, and two Oas-like genes, Oasl1 and Oasl2. Four OAS genes were detected in the pig genome and five OAS genes were found in both the cattle and dog genomes. An OAS3 gene was not found in either the cattle or the pig genome. While two tandemly duplicated OAS-like (OASL) genes were identified in the dog genome, only a single OASL orthologue was found in both the cattle and the pig genomes. The bovine and porcine OASL genes contain premature stop codons and encode truncated proteins, which lack the typical C-terminal double ubiquitin domains. The cDNA sequences of the rat, cattle, pig, and dog OAS genes were amplified, sequenced and compared with each other and with those in the human, mouse, horse, and chicken genomes. Evidence of concerted evolution of paralogous 2′-5′ oligoadenylate synthetase 1 genes was obtained in rodents (Rodentia) and even-toed ungulates (Artiodactyla). Calculations using the nonparametric Kolmogorov-Smirnov test suggested that the homogenization of paralogous OAS1 sequences was due to gene conversion rather than stabilizing selection. Electronic Supplementary Material Electronic Supplementary material is available for this article at and accessible for authorised users. Reviewing Editor: Dr. Martin Kreitman     

Origin and development of oligoadenylate synthetase immune system

Background

Oligoadenylate synthetases (OASs) are widely distributed in Metazoa including sponges, fish, reptiles, birds and mammals and show large variation, with one to twelve members in any given species. Upon double-stranded RNA (dsRNA) binding, avian and mammalian OASs generate the second messenger 2'-5'-linked oligoadenylate (2-5A), which activates ribonuclease L (RNaseL) and blocks viral replication. However, how Metazoa shape their OAS repertoires to keep evolutionary balance to virus infection is largely unknown. We performed comprehensive phylogenetic and functional analyses of OAS genes from evolutionarily lower to higher Metazoa to demonstrate how the OAS repertoires have developed anti-viral activity and diversified their functions.

Results

Ancient Metazoa harbor OAS genes, but lack both upstream and downstream genes of the OAS-related pathways, indicating that ancient OASs are not interferon-induced genes involved in the innate immune system. Compared to OASs of ancient Metazoa (i.e. sponge), the corresponding ones of higher Metazoa present an increasing number of basic residues on the OAS/dsRNA interaction interface. Such an increase of basic residues might improve their binding affinity to dsRNA. Moreover, mutations of functional residues in the active pocket might lead to the fact that higher Metazoan OASs lose the ability to produce 3'-5'-linked oligoadenylate (3-5A) and turn into specific 2-5A synthetases. In addition, we found that multiple rounds of gene duplication and domain coupling events occurred in the OAS family and mutations at functionally critical sites were observed in most new OAS members.

Conclusions

We propose a model for the expansion of OAS members and provide comprehensive evidence of subsequent neo-functionalization and sub-functionalization. Our observations lay the foundation for interrogating the evolutionary transition of ancient OAS genes to host defense genes and provide important information for exploring the unknown function of the OAS gene family.

     

Evolution of Chordate Actin Genes: Evidence from Genomic Organization and Amino Acid Sequences

The origin and evolutionary relationship of actin isoforms was investigated in chordates by isolating and characterizing two new ascidian cytoplasmic and muscle actin genes. The exon–intron organization and sequences of these genes were compared with those of other invertebrate and vertebrate actin genes. The gene HrCA1 encodes a cytoplasmic (nonmuscle)-type actin, whereas the MocuMA2 gene encodes an adult muscle-type actin. Our analysis of these genes showed that intron positions are conserved among the deuterostome actin genes. This suggests that actin gene families evolved from a single actin gene in the ancestral deuterostome. Sequence comparisons and molecular phylogenetic analyses also suggested a close relationship between the ascidian and vertebrate actin isoforms. It was also found that there are two distinct lineages of muscle actin isoforms in ascidians: the larval muscle and adult body-wall isoforms. The four muscle isoforms in vertebrates show a closer relationship to each other than to the ascidian muscle isoforms. Similarly, the two cytoplasmic isoforms in vertebrates show a closer relationship to each other than to the ascidian and echinoderm cytoplasmic isoforms. In contrast, the two types of ascidian muscle actin diverge from each other. The close relationship between the ascidian larval muscle actin and the vertebrate muscle isoforms was supported by both neighbor-joining and maximum parsimony analyses. These results suggest that the chordate ancestor had at least two muscle actin isoforms and that the vertebrate actin isoforms evolved after the separation of the vertebrates and urochordates. Received: 20 June 1996 / Accepted: 16 October 1996     

A new 5′ terminal murine GAPDH exon identified using 5′RACE LaNe

In this work, a ligation-independent, fully gene-specific, nested polymerase chain reaction (PCR) method for the elucidation of 5′ cDNA sequence is described and demonstrated for the first time. Two manifestations of the method, rapid amplification of cDNA ends (RACE) by lariat-dependent nested PCR 5′ (RACE LaNe), at least as simple to perform as conventional RACE, were successfully applied to the murine housekeeping genes phosphoglycerate kinase 1 (PGK1), β-actin (β-ACT), and glyceraldehyde-3-phosphate dehydrogenase (GAPDH) and the alpha thalassemia mental retardation Y homolog (ATRY) gene of the marsupial, Macropus eugenii. Significantly, a new murine GAPDH 5′ exon, separated by 365 kb of intronic sequence from previously annotated GAPDH sequence, was discovered using 5′RACE LaNe.     

Specificity of recognition sequence forEscherichia coli primase

Summary We have surveyed the frequency of each of 64 trinucleotide permutations at every nucleotide frame located from 1 to 15 nucleotides upstream of primer RNA-DNA transition sites mapped within a 1.5 kb region of the bacteriophage lambda genome and a 1.4 kb region of theEscherichia coli genome. We have demonstrated that in both systems initiation of DNA synthesis strongly correlates with a CAG sequence located 11 nucleotides upstream of the DNA start sites. Based on the examination of various reports of the priming reaction catalyzed byE. coli primase in vivo and in vitro, we propose that (i)E. coli primase itself recognizes a 3′GTC 5′ sequence on the template strand, (ii) DnaB helicase releases the specificity ofE. coli primase and, (iii) the consensus recognition sequence forE. coli primase associated with DnaB helicase is 3′PuPyPy 5′.     

Evolution of the 2′-5′-Oligoadenylate Synthetase Family in Eukaryotes and Bacteria

The 2′-5′-oligoadenylate synthetase (OAS) belongs to a nucleotidyl transferase family that includes poly(A) polymerases and CCA-adding enzymes. In mammals and birds, the OAS functions in the interferon system but it is also present in an active form in sponges, which are devoid of the interferon system. In view of these observations, we have pursued the idea that OAS genes could be present in other metazoans and in unicellular organisms as well. We have identified a number of OAS1 genes in annelids, mollusks, a cnidarian, chordates, and unicellular eukaryotes and also found a family of proteins in bacteria that contains the five OAS-specific motifs. This indicates a specific relationship to OAS. The wide distribution of the OAS genes has made it possible to suggest how the OAS1 gene could have evolved from a common ancestor to choanoflagellates and metazoans. Furthermore, we suggest that the OASL may have evolved from an ancestor of cartilaginous fishes, and that the OAS2 and the OAS3 genes evolved from a mammalian ancestor. OAS proteins function in the interferon system in mammals. This system is only found in jawed vertebrates. We therefore suggest that the original function of OAS may differ from its function in the interferon system, and that this original function of OAS is preserved even in OAS genes that code for proteins, which do not have 2′-5′-oligoadenylate synthetase activity.     

Diversity of sponge mitochondrial introns revealed by <Emphasis Type="Italic">cox 1</Emphasis>sequences of Tetillidae

Background  

Animal mitochondrial introns are rare. In sponges and cnidarians they have been found in the cox 1 gene of some spirophorid and homosclerophorid sponges, as well as in the cox 1 and nad 5 genes of some Hexacorallia. Their sporadic distribution has raised a debate as to whether these mobile elements have been vertically or horizontally transmitted among their hosts. The first sponge found to possess a mitochondrial intron was a spirophorid sponge from the Tetillidae family. To better understand the mode of transmission of mitochondrial introns in sponges, we studied cox 1 intron distribution among representatives of this family.     

Deuterostomic Actin Genes and the Definition of the Chordates: cDNA Cloning and Gene Organization for Cephalochordates and Hemichordates

The evolutionary relationship of muscle and nonmuscle actin isoforms in deuterostomia was studied by the isolation and characterization of two actin genes from the cephalochordate Branchiostoma lanceolatum and two from the hemichordate Saccoglossus kowalevskii The Branchiostoma genes specify a muscle and a nonmuscle actin type, respectively. Together with earlier results on muscle actins from vertebrates and urochordates, a N-terminal sequence signature is defined for chordate muscle actins. These diagnostic amino acid residues separate the chordates from the echinoderms and other metazoa. Although the two Saccoglossus actins characterized so far lack the diagnostic residues, in line with the presumptive phylogenetic position of hemichordates outside the chordates, a definitive conclusion can only be expected once the full complement of actin genes of Saccoglossus is established. Comparison of the intron patterns of the various deuterostomic actin genes shows that intron 330-3, which is present in all vertebrate genes, is conspicuously absent from nonvertebrate genes. The possible origin of this intron is discussed. Received: 4 July 1997 / Accepted: 29 August 1997     

Loss of Introns Along the Evolutionary Diversification Pathway of Snake Venom Disintegrins Evidenced by Sequence Analysis of Genomic DNA from <Emphasis Type="Italic">Macrovipera lebetina transmediterranea</Emphasis> and <Emphasis Type="Italic">Echis ocellatus</Emphasis>

Analysis of cDNAs from Macrovipera lebetina transmediterranea (Mlt) and Echis ocellatus (Eo) venom gland libraries encoding disintegrins argued strongly for a common ancestry of the messengers of short disintegrins and those for precursors of dimeric disintegrin chains. We now report the sequence analysis of disintegrin-coding genes from these two vipers. Genomic DNAs for dimeric disintegrin subunits Ml_G1 and Ml_G2 (Mlt) and Eo_D3 (Eo) contain single 1-kb introns exhibiting the 5′-GTAAG (donor)/3′-AG (acceptor) consensus intron splicing signature. On the other hand, the short RTS-disintegrins Ml_G3 (Mlt) and Eo_RTS (Eo) and the short RGD-disintegrin ocellatusin (Eo) are transcribed from intronless genomic DNA sequences, indicating that the evolutionary pathway leading to the emergence of short disintegrins involved the removal of all intronic sequences. The insertion position of the intron within Ml_G1, Ml_G2, and Eo_D3 is conserved in the genes for vertebrate ADAM (A disintegrin and metalloproteinase) protein disintegrin-like domains and within the gene for the medium-size snake disintegrins halystatins 2 and 3. However, a comparative analysis of currently available disintegrin(-like) genes outlines the view that a minimization of both the gene organization and the protein structure underlies the evolution of the snake venom disintegrin family. [Reviewing Editor: Dr. Bryan Fry] Amine Bazaa and Paula Juárez contributed equally to this work and may both be considered first authors.     

Structural Evidence for the Evolution of Pyrogenic Toxin Superantigens

Cells from metazoan organisms are eliminated in a variety of physiological and pathophysiological processes by apoptosis. In this report, we describe the cloning and characterization of molecules from the marine sponges Geodia cydonium and Suberites domuncula, whose domains show a high similarity to those that are found in molecules of the vertebrate Bcl-2 superfamily and of the death receptors. The Bcl-2 proteins contain up to four Bcl-2 homology regions (BH). Two Bcl-2-related molecules have been identified from sponges that are provided with two of those regions, BH1 and BH2, and are termed Bcl-2 homology proteins (BHP). The G. cydonium molecule, BHP1_GC, has a putative size of 28,164, while the related sequence from S. domuncula, BHP1_SD, has a M _r of 24,187. Phylogenetic analyses of the entire two sponge BHPs revealed a high similarity to members of the mammalian Bcl-2 superfamilies and to the Caenorhabditis elegans Ced-9. When the two domains, BH1 and BH2, are analyzed separately, again the highest similarity was found to the members of the Bcl-2 superfamily, but a clearly lower relationship to the C. elegans BH1 and BH2 domains in Ced-9. In unrooted phylogenetic trees the sponge BH1 and BH2 are grouped among the mammalian sequences and are only distantly related to the C. elegans BH domains. The analysis of the gene structure of the G. cydonium BHP showed that the single intron present is located within the BH2 domain at the same position as in C. elegans and rat Bcl-x_L. In addition, a sponge molecule comprising two death domains has been characterized from G. cydonium. The two death domains of the potential proapoptotic molecule GC_DD2, M _r 24,970, share a high similarity with the Fas-FADD/MORT1 domains. A death domain-containing molecule has not been identified in the C. elegans genome. The phylogenetic analysis revealed that the sponge domain originated from an ankyrin building block from which the mammalian Fas-FADD/MORT1 evolved. It is suggested that the apoptotic pathways that involve members of the Bcl-2 superfamily and of the death receptors are already present in the lowest metazoan phylum, the Porifera. Received: 27 July 1999 / Accepted: 28 December 1999     

The 172-kb genomic DNA region of the O. rufipogon yld1.1 locus: comparative sequence analysis with O. sativa ssp. japonica and O. sativa ssp. indica

Common wild rice (Oryza rufipogon) plays an important role by contributing to modern rice breeding. In this paper, we report the sequence and analysis of a 172-kb genomic DNA region of wild rice around the RM5 locus, which is associated with the yield QTL yld1.1. Comparative sequence analysis between orthologous RM5 regions from Oryza sativa ssp. japonica, O. sativa ssp. indica and O. rufipogon revealed a high level of conserved synteny in the content, homology, structure, orientation, and physical distance of all 14 predicted genes. Twelve of the putative genes were supported by matches to proteins with known function, whereas two were predicted by homology to rice and other plant expressed sequence tags or complementary DNAs. The remarkably high level of conservation found in coding, intronic and intergenic regions may indicate high evolutionary selection on the RM5 region. Although our analysis has not defined which gene(s) determine the yld1.1 phenotype, allelic variation and the insertion of transposable elements, among other nucleotide changes, represent potential variation responsible for the yield QTL. However, as suggested previously, two putative receptor-like protein kinase genes remain the key suspects for yld1.1. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Substrate specificity of a recombinant chicken <Emphasis Type="Italic">β</Emphasis>-carotene 15,15′-monooxygenase that converts <Emphasis Type="Italic">β</Emphasis>-carotene into retinal

The recombinant β-carotene 15,15′-monooxygenase from chicken liver was purified as a single 60 kDa band by His-Trap HP and Resource Q chromatography. It had a molecular mass of 240 kDa by gel filtration indicating the native form to be tetramer. The enzyme converted β-carotene under maximal conditions (pH 8.0 and 37°C) with a k _cat of 1.65 min⁻¹ and a K _m of 26 μM and its conversion yield of β-carotene to retinal was 120% (mol mol⁻¹). The enzyme displayed catalytic efficiency and conversion yield for β-carotene, β-cryptoxanthin, β-apo-8′-carotenal, β-apo-4′-carotenal, α-carotene and γ-carotene in decreasing order but not for zeaxanthin, lutein, β-apo-12′-carotenal and lycopene, suggesting that the presence of one unsubstituted β-ionone ring in a substrate with a molecular weight greater than C₃₀ seems to be essential for enzyme activity.     

Intron Conservation in a UV-Specific DNA Repair Gene Encoded by Chlorella Viruses

Large dsDNA-containing chlorella viruses encode a pyrimidine dimer-specific glycosylase (PDG) that initiates repair of UV-induced pyrimidine dimers. The PDG enzyme is a homologue of the bacteriophage T4-encoded endonuclease V. The pdg gene was cloned and sequenced from 42 chlorella viruses isolated over a 12-year period from diverse geographic regions. Surprisingly, the pdg gene from 15 of these 42 viruses contain a 98-nucleotide intron that is 100% conserved among the viruses and another 4 viruses contain an 81-nucleotide intron, in the same position, that is nearly 100% identical (one virus differed by one base). In contrast, the nucleotides in the pdg coding regions (exons) from the intron-containing viruses are 84 to 100% identical. The introns in the pdg gene have 5′-AG/GTATGT and 3′-TTGCAG/AA splice site sequences which are characteristic of nuclear-located, spliceosomal processed pre-mRNA introns. The 100% identity of the 98-nucleotide intron sequence in the 15 viruses and the near-perfect identity of an 81-nucleotide intron sequence in another 4 viruses imply strong selective pressure to maintain the DNA sequence of the intron when it is in the pdg gene. However, the ability of intron-plus and intron-minus viruses to repair UV-damaged DNA in the dark was nearly identical. These findings contradict the widely accepted dogma that intron sequences are more variable than exon sequences. Received: 13 May 1999 / Accepted: 20 August 1999     

Effects of dissolved oxygen levels on survival and growth in vitro of <Emphasis Type="Italic">Haliclona pigmentifera</Emphasis> (Demospongiae)

In vitro sponge cultures are considered as legitimate alternatives for utilizing marine sponges (Porifera) to yield bioactive molecules. Optimization of culture methodologies for enhancing sponge survival is in progress for the identification of the factors regulating sponge survival in vitro. Dissolved oxygen (DO) is an essential factor promoting sponge survival. However, the effects of variable DO levels on the in vitro survival responses of sponges are not fully understood. Hence, we have investigated the effects of variable DO levels on the survival of the marine sponge, Haliclona pigmentifera (Demospongiae), with no external nutritional supplementation in closed type incubator chambers. Our results indicate that, under hypoxic conditions (1.5-2.0 ppm DO), H. pigmentifera with intact ectodermal layers and subtle oscula show adherent growth for 42±3 days. Sponges with prominent oscula, foreign material, and damaged pinacoderm exhibit poor survival under similar conditions. Complete mortality occurs within 2 days under anoxia (<0.3 ppm DO), and survival for a few days has been observed at >4.0 ppm DO without adhesion. Cellular differences between the outer and inner zones and collagen-like extracellular matrix have been identified in adherent sponges. Based on the hypothesis that hypoxia-inducible factor1-α (HIF-1α) is a ubiquitous protein promoting hypoxic survival in animals, we have detected, by Western blot, a protein band corresponding to human HIF-1α-like protein from sponges exposed to hypoxia and to hypoxia-mimicking agents. We thus report, for the first time, adhesive growth and a protein band corresponding to human HIF-1α-like protein in sponges surviving hypoxia in vitro. This work was supported by the “Task Force Network Programme (CMM-004)” of the Council of Scientific and Industrial Research (CSIR), Government of India and by the Department of Ocean Development. CSIR is also acknowledged for providing a Senior Research Fellowship to V.G.G.     

Simple and efficient enzymatic transglycosylation of stevioside by β-cyclodextrin glucanotransferase from <Emphasis Type="Italic">Bacillus firmus</Emphasis>

Stevioside was subjected to 1,4-intermolecular transglycosylation using β-cyclodextrin glucanotransferase (β-CGtase) produced from an alkalophilic strain of Bacillus firmus. The reaction was carried out by traditional, ultrasound-assisted and microwave-assisted techniques. Reaction under microwave conditions was faster and was completed in 1 min yielding two 1,4 transglycosylated products, 4′-O-alpha-d-glycosyl stevioside (I) and 4′′-O-alpha-d-maltosyl stevioside (II) in 66% and 24%, respectively. The optimum transglycosylation occurred by using stevioside (1.24 mmol), β-CD (1.76 mmol) and β-CGtase (2 U/g) under microwave assisted reaction (MAR) in 5 ml sodium phosphate buffer (pH 7) at 50°C and 80 W power. MAR is therefore potentially a useful and economical method for faster transglycosylation of stevioside. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Organization of the genes encoding chalcone synthase in Pisum sativum

To analyze the regulation of defense-related genes by signal molecules produced by phytopathogens, we isolated genes that encode chalcone synthase (CHS) in Pisum sativum. We have obtained seven independent genomic clones that contain at least seven classes of CHS genes, identified by the hybridization analysis to CHS cDNA and by the restriction mapping analysis. Two of the genomic clones (clone 5 and 6) each contain two CHS genes in a tandem repeat. The nucleotide sequence analysis of CHS genomic clone 5 revealed that PsCHS1 and PsCHS2 were corresponding genes of the CHS cDNA clones, pCC6 and pCC2, respectively, as reported earlier. Both genes are interrupted by a single intron of 88 nucleotides with identical sequences, although exonic sequences and 5-flanking sequences are divergent. Nucleotide sequences of the introns in five other classes of CHS genes showed that three classes had an intron of 87 nt with a striking homology to each other, but that the intron of the other two classes of CHS genes showed heterogeneity both in size and nucleotide sequence. 5-upstream regions of PsCHS1 and PsCHS2 did not show sequence homology except the 31 bp identical sequence that contains the CCTACC motif resembling the box-1 sequence. Both PsCHS1 and PsCHS2 genes are shown to be induced by fungal elicitor by a primer extension analysis and a transient transformation analysis using pea protoplasts prepared from suspension cultured-cells.