Analysis of sequence patterns proximal to the stop codons in Oryza sativa

Abstract

A comprehensive analysis of sequence patterns around the stop codons was performed, by using more than 26,000 rice full-length cDNA sequences. Here it is shown that the bias was most outstanding at the position immediately before the stop codons (?1 codon), where the AAC codon was strongly preferred among ANC codons. Compared with other positions, the codon immediately after the stop codons (+1 codon) also displayed an apparent difference, and had a strong consensus for base A at the first, C at the second, and A at the third letters, respectively. Notably, the base biases at the positions directly downstream of the stop codons, such as the +4, +5 and +6 positions, were much stronger than other positions in the 3′-UTR region, suggesting that those base positions might act as an extended stop signal in the process of protein synthesis. Examination of the relationship between sequence pattern and gene expression level, assessed by CAI values and EST counting, revealed a tendency towards bigger base biases for highly expressed genes. It could be inferred that the translation stop signal is possibly involved in many sequence recognition elements other than the stop codons; highly expressed genes should hold strong sequence consensus around the stop codons for efficient translation termination.     

Correlation between sequence conservation of the 5' untranslated region and codon usage bias in Mus musculus genes

The codon adaptation index (CAI) values of all protein-coding sequences of the full-length cDNA libraries of Mus musculus were computed based on the RIKEN mouse full-length cDNA library. We have also computed the extent of consensus in flanking sequences of the initiator ATG codon based on the 'relative entropy' values of respective nucleotide positions (from -20 to +12 bp relative to the initiator ATG codon) for each group of genes classified by CAI values. With regard to the two nucleotides positions (-3 and +4) known to be highly conserved in Kozak's consensus sequence, a clear correlation between CAI values and relative entropy values was observed at position -3 but this was not significant at position +4, although a significant correlation was found at position -1 of the consensus sequence. Further, although no correlation was observed at any additional positions, relative entropy values were very high at positions -4, -6, and -8 in genes with high CAI values. These findings suggest that the extent of conservation in the flanking sequence of the initiator ATG codon including Kozak's consensus sequence was an important factor in modulation of the translation efficiency as well as synonymous codon usage bias particularly in highly expressed genes.     

Mutations in the K+ channel signature sequence.

Potassium channels share a highly conserved stretch of eight amino acids, a K+ channel signature sequence. The conserved sequence falls within the previously defined P-region of voltage-activated K+ channels. In this study we investigate the effect of mutations in the signature sequence of the Shaker channel on K+ selectivity determined under bi-ionic conditions. Nonconservative substitutions of two threonine residues and the tyrosine residue leave selectivity intact. In contrast, mutations at some positions render the channel nonselective among monovalent cations. These findings are consistent with a proposal that the signature sequence contributes to a selectivity filter. Furthermore, the results illustrate that the hydroxyl groups at the third and fourth positions, and the aromatic group at position seven, are not essential in determining K+ selectivity.     

Defining the sequence recognized with BmFTZ-F1, a sequence specific DNA binding factor in the silkworm, Bombyx mori, as revealed by direct sequencing of bound oligonucleotides and gel mobility shift competition analysis.

BmFTZ-F1 is a Bombyx mori homologue of FTZ-F1, a positive regulator of the fushi tarazu gene of Drosophila melanogaster. In order to determine the sequence recognized with this factor, we made three sets of oligonucleotide mixture which contain 4 possible nucleotides at different positions within the previously proposed 12-bp binding consensus sequence. Oligonucleotides which bound to purified BmFTZ-F1 were separated by a gel mobility shift procedure and a binding sequence was determined by direct sequencing through Maxam-Gilbert method. By this analysis, 7 positions showed clear sequence preference and 5 positions showed weak or no sequence preference. The importance of each nucleotide at each position was confirmed by a gel mobility shift competition analysis and results were presented as a quantitative difference in the binding affinity. From these analyses, we conclude that the best binding sequence of BmFTZ-F1 is 5'-PyCAAGGPyCPu-3'. This method may be useful for the determination of a binding sequence of other sequence specific DNA binding factor.     

A new suppressor of a lamB signal sequence mutation, prlZ1, maps to 69 minutes on the Escherichia coli chromosome.

Reversion analysis has been employed to isolate suppressors that restore export of a unique LamB signal sequence mutant. The mutation results in a substitution of Arg for Met at position 19, which prevents LamB export to the outer membrane and leads to a Dex- phenotype. Unlike other LamB signal sequence mutants utilized for reversion analysis, LamB19R becomes stably associated with the inner membrane in an export-specific manner. In this study, Dex+ revertants were selected and various suppressors were isolated. One of the extragenic suppressors, designated prlZ1, was chosen for further study. prlZ1 maps to 69 min on the Escherichia coli chromosome. The suppressor is dominant and SecB dependent. In addition to its effect on lamB19R, prlZ1 suppresses the export defect of signal sequence point mutations at positions 12, 15, and 16, as well as several point mutations in the maltose-binding protein signal sequence. prlZ1 does not suppress deletion mutations in either signal sequence. This pattern of suppression can be explained by interaction of a helical LamB signal sequence with the suppressor.     

Computational identification and sequence analysis of stop codon readthrough genes in Oryza sativa

Using an approach based on the Readthrough Candidate Extraction System (RCES), we extracted 111 candidates from 9620 gene sequences of rice. The results of homology search and sequence analysis demonstrated that these candidates included actual readthrough genes that would be important for further investigating the mechanism of translation termination regulated by readthrough event, and could also give some useful clues for functional genome annotation. Between the candidates and non-candidates of gene sequences in rice, there exist significant base biases at the positions surrounding the stop codons. These positions, especially both -1 and +4, are referred to as part of an extended stop signal. In candidates, G at position -1, and G or C at position +4 are much more favored than that in non-candidates. Both stop sequence patterns, GUAGC and GUGAG, might drive high readthrough efficiency in rice. Secondary structure analysis revealed that the -1 and +1 amino acids around the first stop codon of candidates have a strong bias toward arginine, particularly the +1 position (20.7%), which indicated that the amino acids at the readthrough region being frequently located in the hydrophilic region of beta-turn might be a determinant for efficient translation termination or not.     

Amino acid sequence of S-adenosyl-L-homocysteine hydrolase from Dictyostelium discoideum as deduced from the cDNA sequence

S-Adenosyl-L-homocysteine hydrolase has been cloned from a lambda gt11 cDNA library prepared from Dictyostelium discoideum that had been starved for 3 hours. The sequence of the cloned cDNA was determined and the deduced amino acid sequence was compared to the amino acid sequence of rat AdoHcy hydrolase. When the sequences from the two species were aligned, 74% of the amino acids were in identical positions. If conservative changes were taken into account the homology was 84%. Because differences have been reported in the binding characteristics of NAD+ to the D. discoideum and rat AdoHcy hydrolases, changes in the amino acids of the putative NAD+-binding site were of particular interest. Six changes were observed in this region but the changes appeared to be in regions that are not critical to the three dimensional folding of the NAD+-binding site.     

Extent of the DNA sequence required in integration of staphylococcal bacteriophage L54a.

We characterized the minimum length of the DNA sequence of the attachment sites involved in the integrative recombination of staphylococcal bacteriophage L54a. A DNA fragment carrying the functional viral attachment site (attP) or the bacterial attachment site (attB) was sequentially trimmed, recloned, and tested for integrative recombination in vivo. The size of the functional attP site was at least 228 base pairs (bp) but no more than 235 bp. The left endpoint of the attP site was located to between positions -142 and -140, whereas the right endpoint was located to between positions +86 and +93 with respect to the center of the core sequence. The attB site was located to within a 27-bp sequence, from position -15 to +12, which included the 18-bp core sequence.     

Molecular cloning, nucleotide sequence and expression of the tufB gene encoding elongation factor Tu from Thermus thermophilus HB8

The tufB gene encoding elongation factor Tu (EF-Tu) of Thermus thermophilus HB8 was cloned and expressed. Compared with the known tufA gene of T. thermophilus, nucleotide differences were found at 10 positions out of 1221 nucleotides, and amino acid substitutions were found at 4 positions out of 406 amino acids. The tufB product was 70.9% homologous to the corresponding sequence of the tufB product of E. coli. The G+C content of the third base of the codon in the tufB gene was 84.8% and G was especially preferred in this position.     

Comprehensive sequence analysis of translation termination sites in various eukaryotes

Recent investigations into the translation termination sites of various organisms have revealed that not only stop codons but also sequences around stop codons have an effect on translation termination. To investigate the relationship between these sequence patterns and translation as well as its termination efficiency, we analysed the correlation between strength of consensus and translation efficiency, as predicted according to Codon Adaptation Index (CAI) value. We used RIKEN full-length mouse cDNA sequences and ten other eukaryotic UniGene datasets from NCBI for the analyses. First, we conducted sequence profile analyses following translation termination sites. We found base G and A at position +1 as a strong consensus for mouse cDNA. A similar consensus was found for other mammals, such as Homo sapiens, Rattus norvegicus and Bos taurus. However, some plants had different consensus sequences. We then analysed the correlation between the strength of consensus at each position and the codon biases of whole coding regions, using information content and CAI value. The results showed that in mouse cDNA, CAI value had a positive correlation with information content at positions +1. We also found that, for positions with strong consensus, the strength of the consensus is likely to have a positive correlation with CAI value in some other eukaryotes. Along with these observations, biological insights into the relationship between gene expression level, codon biases and consensus sequence around stop codons will be discussed.     

The influence of flanking sequence on the O-glycosylation of threonine in vitro.

To investigate the influence of flanking amino acid sequence on the O-glycosylation of a single threonine residue in vitro, we have examined a series of 52 related peptides. The substrates were based upon a sequence from human von Willebrand factor which is known to be glycosylated in vivo (-6PHMAQVTVGPGL+5). Each residue of the parent peptide was substituted, in turn, with isoleucine, alanine, proline, glutamic acid, or arginine. Peptides were glycosylated using a UDP-GalNAc:polypeptide N-acetylgalactosaminyltransferase purified 15,000-fold from bovine colostrum by chromatography on DEAE-Sephacel, SP-Sephadex, Sephacryl S-300, Affi-Gel Blue, and 5-mercuri-UDP-GalNAc thiopropyl-Sepharose. Single amino acid changes in the sequences flanking the threonine could profoundly alter the glycosylation of the substrate peptides. Substitution of any amino acid tested at positions +3, -3, and -2 markedly decreased O-glycosylation, as did the presence of a charged residue at position -1. The substitution of amino acids at the other positions of the peptide substrate had little effect on the incorporation of GalNAc. Statistical analysis of sequences flanking known glycosylated threonine and serine residues suggests that they should be glycosylated with equal efficiency in the same sequence context (O'Connell et al., 1991). However, the bovine colostrum transferase failed to glycosylate a peptide derived from human erythropoietin which contains a serine that is glycosylated in vivo (-5PPDAASAAPLR+5). When a threonine was substituted for the serine in this peptide (-5PPDAATAAPLR+5), the substrate proved to be an excellent acceptor of GalNAc. These observations indicate that although flanking amino acid sequence is important for the O-glycosylation of specific hydroxyamino acids, discrete threonine- and serine-specific transferases may exist.     

Analysis of the polyadenylation consensus sequence context in the genes of nuclear encoded mitochondrial proteins

A compilation of the pre-mRNA ends of the genes of nuclear encoded mitochondrial proteins resulted in a consensus sequence of the type (T/A)NTTNNNNNTTTNAATAAA. Nucleotide positions +8, +13, +14, +16 and +17 downstream of the AATAAA sequence show also a predominance of nucleotide T. This consensus sequence suggests the importance of the immediate surroundings of the cannonical polyadenylation signal sequence AATAAA on the efficiency of the cleavage and polyadenylation of this specific group of pre-mRNAs.     

Effects on RNAi of the tight structure, sequence and position of the targeted region

RNA interference (RNAi) is a gene-silencing phenomenon that involves the double-stranded RNA-mediated cleavage of mRNA, and small interfering RNAs (siRNAs) can cause RNAi in mammalian cells. There have been many attempts to clarify the mechanism of RNAi, but information about the relationship between the sequence and structure, in particular, a tight structure, of the target RNA and the activities of siRNAs are limited. In the present study, we examined this relationship by introducing the TAR element, which adopts a very stable secondary structure, at different positions within target RNAs. Our results suggested that the activities of siRNAs were affected by the tight stem–loop structure of TAR. In contrast, the position of the target within the mRNA, the binding of the Tat protein to the TAR, and the location of the target within a translated or a noncoding region had only marginal effects on RNAi. When the target sequence was placed in two different orientations, only one orientation had a significant effect on the activities of siRNA, demonstrating that the presence of certain nucleotides at some specific positions was favorable for RNAi. Systematic analysis of 47 different sites within 47 plasmids under identical conditions indicated that it is the target sequence itself, rather than its location, that is the major determinant of siRNA activity.     

Calorimetric examination of high-affinity Src SH2 domain-tyrosyl phosphopeptide binding: dissection of the phosphopeptide sequence specificity and coupling energetics

SH2 domains are protein modules which interact with specific tyrosine phosphorylated sequences in target proteins. The SH2 domain of the Src kinase binds with high affinity to a tyrosine phosphorylated peptide containing the amino acids Glu, Glu, and Ile (EEI) at the positions +1, +2, and +3 C-terminal to the phosphotyrosine, respectively. To investigate the degree of selectivity of the Src SH2 domain for each amino acid of the EEI motif, the binding thermodynamics of a panel of substitutions at the +1 (Gln, Asp, Ala, Gly), +2 (Gln, Asp, Ala, Gly), and +3 (Leu, Val, Ala, Gly) positions were examined using titration microcalorimetry. It was revealed that the Src SH2 domain is insensitive (DeltaDeltaG degrees 相似文献