Convergent origin of repeats in genes coding for globular proteins. An analysis of the factors determining the presence of inverted and symmetrical repeats

The factors, determining the presence of inverted and symmetrical repeats in genes coding for globular proteins, have been analysed. An interesting property of genetical code has been revealed in the analysis of symmetrical repeats: the pairs of symmetrical codons corresponded to pairs of amino acids with mostly similar physical-chemical parameters. This property may explain the presence of symmetrical repeats and palindromes only in genes coding for beta-structural proteins-polypeptides, where amino acids with similar physical-chemical properties occupy symmetrical positions. A stochastic model of evolution of polynucleotide sequences has been used for analysis of inverted repeats. The modelling demonstrated that only limiting of sequences (uneven frequencies of used codons) is enough for arising of nonrandom inverted repeats in genes.     

Agmenellum quadruplicatum M.AquI, a novel modification methylase.

The complete type II modification methylase of Agmenellum quadruplicatum was cloned in Escherichia coli as an R.Sau3A fragment of approximately 4.5 kilobases. The coding sequence was contained in a stretch of 1,156 base pairs which was organized into two parallel, partly overlapping open reading frames of 248 and 139 codons. In vivo complementation experiments showed that the synthesis of both predicted peptides was required for full methylase activity. The amino acid sequences were considerably similar to regions of other deoxycytidylate methylases.     

The positive regulatory function of the 5''-proximal open reading frames in GCN4 mRNA can be mimicked by heterologous, short coding sequences.

Translational control of GCN4 expression in the yeast Saccharomyces cerevisiae is mediated by multiple AUG codons present in the leader of GCN4 mRNA, each of which initiates a short open reading frame of only two or three codons. Upstream AUG codons 3 and 4 are required to repress GCN4 expression in normal growth conditions; AUG codons 1 and 2 are needed to overcome this repression in amino acid starvation conditions. We show that the regulatory function of AUG codons 1 and 2 can be qualitatively mimicked by the AUG codons of two heterologous upstream open reading frames (URFs) containing the initiation regions of the yeast genes PGK and TRP1. These AUG codons inhibit GCN4 expression when present singly in the mRNA leader; however, they stimulate GCN4 expression in derepressing conditions when inserted upstream from AUG codons 3 and 4. This finding supports the idea that AUG codons 1 and 2 function in the control mechanism as translation initiation sites and further suggests that suppression of the inhibitory effects of AUG codons 3 and 4 is a general consequence of the translation of URF 1 and 2 sequences upstream. Several observations suggest that AUG codons 3 and 4 are efficient initiation sites; however, these sequences do not act as positive regulatory elements when placed upstream from URF 1. This result suggests that efficient translation is only one of the important properties of the 5' proximal URFs in GCN4 mRNA. We propose that a second property is the ability to permit reinitiation following termination of translation and that URF 1 is optimized for this regulatory function.     

Nucleotide sequence of human endogenous retrovirus genome related to the mouse mammary tumor virus genome.

We determined the complete nucleotide sequence of the human endogenous retrovirus genome HERV-K10 isolated as the sequence homologous to the Syrian hamster intracisternal A-particle (type A retrovirus) genome. HERV-K10 is 9,179 base pairs long with long terminal repeats of 968 base pairs at both ends; a sequence 290 base pairs long, however, was found to be deleted. It was concluded that a composite genome having the 290-base-pair fragment is the prototype HERV-K provirus gag (666 codons), protease (334 codons), pol (937 codons), and env (618 codons) genes. The size of the protease gene product of HERV-K is essentially the same as that of A- and D-type oncoviruses but nearly twice that of other retroviruses. A comparison of the deduced amino acid sequences encoded by the pol region showed HERV-K to be closely related to types A and D retroviruses and even more so to type B retrovirus. It was noted that the env gene product of HERV-K structurally resembles the mouse mammary tumor virus (type B retrovirus) env protein, and the possible expression of the HERV-K env gene in human breast cancer cells is discussed.     

Modified base compositions at degenerate positions of a mutagenic oligonucleotide enhance randomness in site-saturation mutagenesis.

Site-saturation mutagenesis, using degenerate oligonucleotide primers, is a frequently used method in introducing various mutations in a selected target codon. Oligonucleotides that are synthesized using equimolar concentrations of nucleoside phosphoramidites (dA, dC, dG, dT) in the positions to be saturated, result in a mutant population that is biased towards the original nucleotides. We found that this bias could be eliminated by modifying the concentrations of nucleoside phosphoramidites during the oligonucleotide synthesis. We synthesized eight degenerate oligonucleotides to saturate eight different codons, and sequenced a total of 344 mutagenized codons. In six of these eight oligonucleotides, we reduced to varying extents the concentrations of those nucleotides in the target positions that would form base pairs with the template. From the data, we analyzed the effects of different base compositions in the oligonucleotides when mutagenizing different codons, the influence of the positions of mismatches, and the significance of different non-Watson-Crick base pairs. Based on these results, we suggest levels to which different phosphoramidites should be reduced when synthesizing oligonucleotides for site-saturation mutagenesis.     

Analysis of the compositional biases in Plasmodium falciparum genome and proteome using Arabidopsis thaliana as a reference

Comparative genomic analysis of the malaria causative agent, Plasmodium falciparum, with other eukaryotes for which the complete genome is available, revealed that the genome from P. falciparum was more similar to the genome of a plant, Arabidopsis thaliana, than to other non-apicomplexan taxa. Plant-like sequences are thought to result from horizontal gene transfers after a secondary endosymbiosis involving an algal ancestor. The use of the A. thaliana genome and proteome as a reference gives an opportunity to refine our understanding of the extreme compositional bias in the P. falciparum genome that leads to a proteome-wide amino acid bias. A set of pairs of non-redundant protein homologues was selected owing to rigorous genome-wide sequence comparison methods. The introduction of A. thaliana as a reference was a mean to weight the magnitude of the protein evolutionary divergence in P. falciparum. The correlation of the amino acid proportions with evolutionary time supports the hypothesis that amino acids encoded by GC-rich codons are directionally substituted into amino acids encoded by AT-rich codons in the P. falciparum proteome. The long-term deviation of codons in malarial sequences appears as a possible consequence of a genome-wide tri-nucleotidic signature imprinting. Additionally, this study suggests possible working guidelines to improve the accuracy of P. falciparum sequence comparisons, for homology searches and phylogenetic studies.     

IS10 transposase mutations that specifically alter target site recognition.

IS10 inserts preferentially into particular hotspots. We describe here mutations of IS10 transposase, called 'ATS' that confer Altered Target Specificity. These mutations yield a general relaxation in target specificity but do not affect other aspects of transposition. Thus, the preference for specific nucleotide sequences at the target site can be cleanly separated from other steps of the transposition reaction. Eleven ATS mutations identified in a genetic screen occur at only two codons in transposase, one in each of two regions of the protein previously implicated in target site interactions (Patch I and Patch II). Genetic analysis suggests that mutations at the two ATS codons affect the same specific function of transposase, thus raising the possibility that Patch I and Patch II interact. For wild-type IS10, insertion specificity is determined in part by a specific 6 bp consensus sequence and in part by the immediately adjacent sequence context of the target DNA. The ATS mutations do not qualitatively alter the hierarchy with which base pairs are recognized in the consensus sequence; instead, sites selected by ATS transposase exhibit a reduction in the degree to which certain base pairs are preferred over others. Models for the basis of this phenotype are discussed.     

Rare codons in E. coli and S. typhimurium signal sequences

Codon usage has been examined in the signal sequences of 27 genes encoding proteins which possess leader peptides, and are inner-membrane located or exported. The results have been compared with codon usage in the corresponding coding sequences of most of the mature proteins. A bias is observed in the usage of rare codons for two of the three hydrophobic amino acids for which there are rare codons. Since hydrophobic residues are predominant in leader peptides, we suggest that a resulting concentration of rare codons in the signal sequence may play a role (or have played a role in the evolutionary past) in the secretion process by delaying translation.     

Recognizing shorter coding regions of human genes based on the statistics of stop codons.

With the quick progress of the Human Genome Project, a great amount of uncharacterized DNA sequences needs to be annotated copiously by better algorithms. Recognizing shorter coding sequences of human genes is one of the most important problems in gene recognition, which is not yet completely solved. This paper is devoted to solving the issue using a new method. The distributions of the three stop codons, i.e., TAA, TAG and TGA, in three phases along coding, noncoding, and intergenic sequences are studied in detail. Using the obtained distributions and other coding measures, a new algorithm for the recognition of shorter coding sequences of human genes is developed. The accuracy of the algorithm is tested based on a larger database of human genes. It is found that the average accuracy achieved is as high as 92.1% for the sequences with length of 192 base pairs, which is confirmed by sixfold cross-validation tests. It is hoped that by incorporating the present method with some existing algorithms, the accuracy for identifying human genes from unannotated sequences would be increased.     

Efficient expression of gene variants that harbour AGA codons next to the initiation codon

In an effort to improve the knowledge about the rules which direct the effect of the early ORF sequences on translation efficiency, we have analyzed the effect of pairs of the six arginine codons at the second and third positions on the expression of lacZ variants. Whereas the pairs of identical AGA or AGG codons were favorable for the gene expression, identical pairs of each of the four CGN codons were very inefficient. This result was unexpected because tandems of AGA or AGG codons located in more internal gene positions provoke deficient expression whilst internally located CGU and CGC are the most abundant and efficiently translated arginine codons. The mixed combinations of AGA and each of the CGN codons usually resulted in efficient rates of lacZ expression independently of the peptidyl-tRNA propensity to dissociate from the ribosome. Thus, the variant harboring the pair of AGA codons was expressed as efficiently as the variant carrying a pair of AAA codons in the same positions, a configuration reported as one of the most common and efficient for gene expression. We explain these results assuming that the presence of adenines in these early positions enhance gene expression. As expected, specific mRNA levels correlated with the intensity of lacZ expression for each variant. However, the induction of lacZ AGA AGA gene in pth cells accumulated peptidyl-tRNA^Arg4 as well as a short 5′-proximal lacZ mRNA fragment suggesting ribosome stalling due to depletion of aminoacylated-tRNA^Arg4.     

Determination of the amino acid requirements for a protein hinge in triosephosphate isomerase.

We have determined the sequence requirements for a protein hinge in triosephosphate isomerase. The codons encoding the hinge at the C-terminus of the active-site lid of triosephosphate isomerase were replaced with a genetic library of all possible 8,000 amino acid combinations. The most active of these 8,000 mutants were selected using in vivo complementation of a triosephosphate isomerase deficient strain of E. coli, DF502. Approximately 3% of the mutants complement DF502 with an activity that is above 70% of wild-type activity. The sequences of these hinge mutants reveal that the solutions to the hinge flexibility problem are varied. Moreover, these preferences are sequence dependent; that is, certain pairs occur frequently. They fall into six families of similar sequences. In addition to the hinge sequences expected on the basis of phylogenetic analysis, we selected three new families of 3-amino-acid hinges: X(A/S)(L/K/M), X(aromatic/beta-branched)(L/K), and XP(S/N). The absence of these hinge families in the more than 60 known species of triosephosphate isomerase suggests that during evolution, not all of sequence space is sampled, perhaps because there is no neutral mutation pathway to access the other families.     

Complete sequence of bluetongue virus L2 RNA that codes for the antigen recognized by neutralizing antibodies.

The complete sequence of the RNA which encodes the major outer-shell-neutralizing antigen (VP2) of bluetongue virus serotype 10 was determined from overlapping cDNA clones inserted into pBR322. The segment L2 RNA was 2,926 base pairs long (1.87 X 10(6) daltons) and had, in one strand, an open reading frame capable of coding for a protein that had a calculated size of 111,122 daltons (956 amino acids) and a +11.5 net charge. The coding strands of both the L2 gene and the group-specific L3 gene of bluetongue virus serotype 17 (M. Purdy, J. Petre, and P. Roy, J. Virol. 51:754-759, 1984) had common sequences of some six nucleotides at their 5' termini (namely, GUUAAA...) and eight nucleotides (namely, ...ACACUUAC) at their 3' termini. Both had short 5' noncoding regions with AUG codons at residues 20 to 22 (L2) and 18 to 20 (L3). The sequences flanking these AUG codons were similar (A/GCCAUGG). The 3' noncoding regions were longer (36 nucleotides for L2, 49 nucleotides for L3). The predicted amino acid sequence of the L2, compared with the similarly sized L3 gene product, was rich in cysteine residues and charged amino acids.     

Genetic analysis of the yopE region of Yersinia spp.: identification of a novel conserved locus, yerA, regulating yopE expression.

The yopE gene of Yersinia pseudotuberculosis was recently sequenced, and YopE was identified as an indispensable virulence determinant when tested in a mouse model (A. Forsberg and H. Wolf-Watz, Mol. Microbiol. 2:121-133, 1988). In the study described here, the DNA sequences of the yopE genes of Yersinia pestis EV76 and Yersinia enterocolitica 8081 were determined and compared with that of the Y. pseudotuberculosis gene. Only two codons were found to differ, both leading to amino acid replacements, when the gene from Y. pestis was compared. These two replacements were also present in the gene from Y. enterocolitica; in addition, 18 other codons were found to differ. Thirteen of these substitutions led to amino acid replacements. Downstream of the yopE gene, the plasmid partition locus par was found to be conserved in all three species. In Y. enterocolitica 8081, the sequence homology was interrupted by a putative insertion sequence element inserted between the yopE gene and the par region at a position only 5 base pairs downstream of the yopE stop codon. Upstream of the yopE gene, 620 base pairs were conserved in the three species. This region contained a 130-amino-acid-long open reading frame reading in the opposite direction to the yopE gene and expressed a 14-kilodalton protein in minicells. An insertion mutation in this region constructed in Y. pseudotuberculosis expressed significantly lower amounts of YopE protein in vitro than did the corresponding wild type. The expression level could be restored by transcomplementation. This new locus was designated yerA, for yopE-regulating gene A. The yerA mutant was avirulent when mice were challenged by oral infection.     

Unusually biased nucleotide sequences on sense strands of Flavobacterium sp. genes produce nonstop frames on the corresponding antisense strands.

From investigation of eight Flavobacterium sp. genes encoding enzyme proteins, it was found that six genes had nonstop frames (NSFs) on the antisense strands, and base sequences of the genes are mainly composed of repeating triplet sequence(s), 5'-GNC-3' (where G and C are guanine and cytosine, and N is either of the four bases), in the reading frames. Thus, we concluded that the biased nucleotide sequences on the sense strands produce NSFs on the corresponding antisense strands. Furthermore, from the precise alignments of both nucleotide and amino acid sequences of two related Flavobacterium sp. genes, nyIB and nyIB', it was found that base replacements might have occurred symmetrically in the codons. That is, transversions between G and C were observed at high frequencies at the first and third positions of codons, but not at the second positions. At the first position, AG base transitions were observed much more than similar CT transitions, whereas CT transitions were found at the third positions at a relatively high frequency. These suggest that symmetrical base replacements in codons might be the main contribution to evolution in Flavobacterium sp. genes.