The sequence flanking translational initiation site in protozoa.

Nucleotide sequences flanking the translational initiation site were compiled from protozoan nuclear genes and were compared in every protozoan group. The entire 5'-untranslated sequences were very rich in A- and T-residues, but poor in G- and C-residues in most protozoan genes except for the flagellated ones. The sequence AAAAATTTTTAAAATTTAAAATGANAT emerged as a consensus sequence flanking the initiation site in the major protozoan group, although the sequences upstream from -4 (four nucleotides upstream from the ATG codon) were divergent among ciliates, sarcodinians, and sporozoans. On the other hand, the consensus sequence for flagellates was revealed to be a simple feature. Only the nucleotide position -3 was occupied with a high frequency of A-residue, in other positions it appeared randomly. These facts suggest that the strong preference for A-residue at the position -3 is a universal feature in nuclear genes for all eukaryotes.     

G and T Nucleotide Contents Show Specie-Invariant Negative Correlation for All Three Codon Positions

Abstract

The nucleotide contents of the three codon positions show a number of statistical pairwise correlations, some of which are universal for all analysed genomes. Among the most prominent of these correlations are negative correlations between G and T contents found in genes of all species analysed. The pair A/C, which is complementary to G/T shows similar negative correlation in genes of most species. In the genes of several species including all mammalian genes studied, positive correlations between A and T contents, and G and C contents are found. Since these regularities are observed in all three codon positions they are connected with amino-acid content of proteins. Such correlations may origin from features of the mutation process or/and translation reading frame check. The well-known bias of the preference for G in the first codon position and its deficiency in the second is accompanied by opposite bias in T content. In the third codon position there is no general nucleotide preference, but its content is often biased with regard to GC content of the gene. G and T contents in this case are always shifted in the opposite directions Several ideas are drawn to explain this preference.     

Positive and negative regulatory DNA elements including a CCArGG box are involved in the cell type-specific expression of the human muscle dystrophin gene.

The muscle-specific promoter of the dystrophin gene is active in skeletal, cardiac, and smooth muscles and is specifically stimulated during differentiation of myoblasts into multinucleated myotubes. An 850-base pair (bp) DNA fragment upstream from the cap site is able to confer a partial muscle specificity to a reporter gene. The region between -850 and -140 bp includes nonspecific negative and positive regulatory sequences. A continuous stretch of 140 bp upstream from the cap site exhibits a striking conservation between rodents and human (93% homology) and still retains muscle preference of expression. It contains two putative binding sites for factors involved in regulation of other muscle-specific genes, a CCArGG box and an E box. This latter element, however, is unable to confer the ability to be transactivated by MyoD1 to the dystrophin promoter. The -140-bp promoter fragment exhibits antagonist effects contributed by one inhibiting sequence (nucleotide -140/-96), active in all cell types, and one activating region, from nucleotide -96 to the cap site, sufficient to confer a muscle preference of expression, in which the CCArGG box seems to play a major role.     

The complete nucleotide sequence of the Pseudomonas gene coding for carboxypeptidase G2

The complete nucleotide sequence of the Pseudomonas chromosomal gene coding for the enzyme carboxypeptidase G2 (CPG2) has been determined. The nucleotide sequence obtained has been confirmed by comparing the predicted amino acid sequence with that of randomly derived peptide fragments and by N-terminal sequencing of the purified protein. The gene has been shown to code for a 22 amino acid signal peptide at its N-terminus which closely resembles the signal peptides of other secreted proteins. An alternative 36 amino acid signal peptide which may function in Pseudomonas has also been identified. The codon utilisation of the gene is influenced by the high G + C (67.2%) content of the DNA and exhibits a 92.8% preference for codons ending in G or C. This unusual codon preference may contribute to the generally observed weak expression of Pseudomonas genes in Escherichia coli. A region of DNA upstream of the structural gene has also been sequenced and a ribosome binding site and two putative promoter sequences identified.     

Translation initiation AUG context varies with codon usage bias and gene length in Drosophila melanogaster

The relationship between the codon usage bias and the sequence context surrounding the AUG translation initiation codon was examined in 1100 Drosophila melanogaster mRNA sequences. The codon usage bias measured by the "codon adaptation index" (CAI), and the effectiveness of the AUG context for translation initiation assessed by the "AUG context adaptation index" (AUGCAI), showed a significant positive relationship (correlation coefficient: r = 0.34, p <0.0001), indicating that these two factors are evolutionally under a similar natural selection constraint at the translational level. The importance of each position of the AUG context in relation to codon usage bias was examined, and the preference for the nucleotide at the -13, -12, -11, -10, -7, -6, -5, -4, -3, -2, and -1 positions showed a significant positive correlation to the codon usage bias, suggesting the action of natural selection on these very specific positions of the Drosophila genome. The relationship between AUGCAI value and gene length was also examined, and a significant negative relationship was found (r = -0.15, p <0.0001), suggesting a general tendency of higher expressivity of shorter genes, and of lower expressivity of longer genes in D. melanogaster.     

A new method for estimating synonymous and nonsynonymous rates of nucleotide substitution considering the relative likelihood of nucleotide and codon changes

A new method is proposed for estimating the number of synonymous and nonsynonymous nucleotide substitutions between homologous genes. In this method, a nucleotide site is classified as nondegenerate, twofold degenerate, or fourfold degenerate, depending on how often nucleotide substitutions will result in amino acid replacement; nucleotide changes are classified as either transitional or transversional, and changes between codons are assumed to occur with different probabilities, which are determined by their relative frequencies among more than 3,000 changes in mammalian genes. The method is applied to a large number of mammalian genes. The rate of nonsynonymous substitution is extremely variable among genes; it ranges from 0.004 X 10(-9) (histone H4) to 2.80 X 10(-9) (interferon gamma), with a mean of 0.88 X 10(-9) substitutions per nonsynonymous site per year. The rate of synonymous substitution is also variable among genes; the highest rate is three to four times higher than the lowest one, with a mean of 4.7 X 10(-9) substitutions per synonymous site per year. The rate of nucleotide substitution is lowest at nondegenerate sites (the average being 0.94 X 10(-9), intermediate at twofold degenerate sites (2.26 X 10(-9)). and highest at fourfold degenerate sites (4.2 X 10(-9)). The implication of our results for the mechanisms of DNA evolution and that of the relative likelihood of codon interchanges in parsimonious phylogenetic reconstruction are discussed.     

An hour in the morning is worth two in the evening: association of morning component of morningness–eveningness with single nucleotide polymorphisms in circadian clock genes

Sub-constructs of morning–evening preference might be differentially related to polymorphisms in circadian clock genes. We previously reported significant association between a single nucleotide polymorphism in PER3 (rs2640909) and Morning but not Evening Lateness scale of the Sleep–Wake Pattern Assessment Questionnaire. To further explore such a scale-specific relationship, seven single nucleotide polymorphisms in five circadian clock genes were studied using exploratory and confirmatory samples (in total, n = 698). The association of rs2640909 with Morning Lateness scale was not replicated in the confirmatory sample but remained significant in the merged sample. Moreover, we found and confirmed an association of this scale with rs1159814 in RORα. The results provided further evidence for differential relationship of polymorphisms in circadian clock genes with morning and evening components of morning–evening preference. We also suggested possibility to take into account the pattern of geographic variation in allele frequency for prioritization of circadian clock polymorphisms in candidate gene studies.     

Comparison of the differential context-dependence of DNA deamination by APOBEC enzymes: correlation with mutation spectra in vivo

To investigate the extent to which in vivo mutation spectra might reflect the intrinsic specificities of active mutators, genetic and biochemical assays were used to analyse the DNA target specificities of cytidine deaminases of the APOBEC family. The results reveal the critical importance of nucleotides immediately 5' of the targeted C for the specificity of all three enzymes studied (AID, APOBEC1 and APOBEC3G). At position -1, APOBEC1 showed a marked preference for dT, AID for dA/dG and APOBEC3G a strong preference for dC. Furthermore, AID and APOBEC3G showed distinct dependence on the nucleotide at position -2 with dA/dT being favoured by AID and dC by APOBEC3G. Most if not all activity of the recombinant deaminases on free dC could be attributed to low-level contamination by host enzymes. The target preference of APOBEC3G supports it being a major but possibly not sole contributor to HIV hypermutation without making it a dominant contribution to general HIV sequence variation. The specificity of AID as deduced from the genetic assay (which relies on inactivation of sacB of Bacillus subtilis) agrees well with that deduced by Pham et al. using an in vitro assay although we postulate that major intrinsic mutational hotspots in immunoglobulin V genes in vivo might reflect favoured sites of AID action being generated by proximal DNA targets located on opposite DNA strands. The target specificity of AID also accords with the spectrum of mutations observed in B lymphoma-associated oncogenes. The possibility of deaminase involvement in non-lymphoid human tumours is hinted at by tissue-specific differences in the spectra of dC transitions in tumour-suppressor genes. Thus, the patterns of hypermutation in antibodies and retroviruses owe much to the intrinsic sequence preferences of the AID/APOBEC family of DNA deaminases: analogous biases might also contribute to the spectra of cancer-associated mutation.     

Maize Mu transposons are targeted to the 5' untranslated region of the gl8 gene and sequences flanking Mu target-site duplications exhibit nonrandom nucleotide composition throughout the genome

The widespread use of the maize Mutator (Mu) system to generate mutants exploits the preference of Mu transposons to insert into genic regions. However, little is known about the specificity of Mu insertions within genes. Analysis of 79 independently isolated Mu-induced alleles at the gl8 locus established that at least 75 contain Mu insertions. Analysis of the terminal inverted repeats (TIRs) of the inserted transposons defined three new Mu transposons: Mu10, Mu 11, and Mu12. A large percentage (>80%) of the insertions are located in the 5' untranslated region (UTR) of the gl8 gene. Ten positions within the 5' UTR experienced multiple independent Mu insertions. Analyses of the nucleotide composition of the 9-bp TSD and the sequences directly flanking the TSD reveals that the nucleotide composition of Mu insertion sites differs dramatically from that of random DNA. In particular, the frequencies at which C's and G's are observed at positions -2 and +2 (relative to the TSD) are substantially higher than expected. Insertion sites of 315 RescueMu insertions displayed the same nonrandom nucleotide composition observed for the gl8-Mu alleles. Hence, this study provides strong evidence for the involvement of sequences flanking the TSD in Mu insertion-site selection.     

A single nucleotide linked to a switch in metal ion reactivity preference in the HDV ribozymes

The two ribozymes of hepatitis delta virus (HDV) cleave faster in divalent metal ions than in monovalent cations, and a variety of divalent metal ions can act as catalysts in supporting these higher rates. Although the ribozymes are closely related in sequence and structure, they display a different metal ion preference; the genomic form cleaves moderately faster in Mg2+ than in Ca2+ while the reverse is true for the antigenomic ribozyme. This difference raises questions about understanding the catalytic role of the metal ion in the reaction. We found that the metal ion reactivity preference correlated with the identity of a single nucleotide 5' of the cleavage site (-1 position). It is a U in the genomic sequence and a C in the antigenomic sequence. With both ribozymes, the reactivity preference for Mg2+ and Ca2+ could be reversed with a change at this position (C to U or U to C). Moreover, with an A at position -1, there was a relative increase in cleavage rates in low concentrations of Mn2+ for both ribozymes. Metal ion reactivity preference was also linked to changes in pH, and the pH-rate profiles could be shifted with nucleotide changes at position -1. Together, the data provide biochemical evidence in support of an organized active site, as seen in the crystal structures, where at least one metal ion, an ionizable group, and the conformation of the phosphate backbone at the cleavage site interact in concert to promote cleavage.     

Translation of chloroplast-encoded mRNA: potential initiation and termination signals.

A survey of 196 protein-coding chloroplast DNA sequences demonstrated the preference for AUG and UAA codons for initiation and termination of translation, respectively. As in prokaryotes at every nucleotide position from -25 to +25 (AUG is +1 to +3) and for 25 nucleotides 5' and 3' to the termination codon an A or U is predominant, except for C at +5 and G at +22. A Shine-Dalgarno (SD) sequence (GGAGG or tri- or tetranucleotide variant) was found within 100 bp 5' to the AUG codon in 92% of the genes. In 40% of these cases, the location of the SD sequence was similar to that of the consensus for prokaryotes (-12 to -7 5' to AUG), presumed to be optimal for translation initiation. A SD sequence could not be located in 6% of the chloroplast sequences. We propose that mRNA secondary structures may be required for the relocation of a distal SD sequences to within the optimal region (-12 to -7) for initiation of translation. We further suggest that termination at UGA codons in chloroplast genes may occur by a mechanism, involving 16S rRNA secondary structure, which has been proposed for UGA termination in E. coli.     

Cleavage site determinants in the mammalian polyadenylation signal.

Using a series of position and nucleotide variants of the SV40 late polyadenylation signal we have demonstrated that three sequence elements determine the precise site of 3-end cleavage in mammalian pre-mRNAs: an upstream AAUAAA element, a down-stream U-rich element consisting of five nucleotides, at least four of which are uridine, and a nucleotide preference at the site of cleavage in the order A > U > C >> G. Cleavage occurs no closer than 11 bases, but no further than 23 bases from the AAUAAA element. The downstream U-rich element is usually located 10-30 bases from the cleavage site. The relative position of the AAUAAA and the U-rich elements define the approximate region within a 13 base domain in which cleavage will occur. The exact position of cleavage is then determined by the local nucleotide sequence in the order of preference noted above. This model accounts for nearly three quarters of polyadenylation signals surveyed and is consistent with previous experimental observations.     

Nucleotides downstream of start codons show marked non-randomness in Escherichia coli but not in Bacillus subtilis

This study aimed at measuring the nucleotide non-randomness in the region downstream of start codons in bacterial genes and to see if the non-randomness differs between biased and unbiased genes, in terms of the effective number of codons (Nc) and the codon adaptation index (CAI). In Escherichia coli, there was a marked elevation in nucleotide conservation for the genes having low Nc-values compared to the genes having high Nc-values, i.e the more biased genes showed a higher level of non-randomness. Likewise, the genes displaying high CAI-values showed stronger nucleotide conservation than the genes of low CAI-values. This elevated conservation is visible up to approximately 15-17 nucleotides downstream of the start codon, after which there is little difference. This indicates that there may be distinct selectional mechanisms acting upon the first 5-6 codons within genes in E. coli. In B. subtilis, these effects are less pronounced, if present at all. Furthermore, analyses of codons used in this region were not in support of the hypothesis that the elevation in nucleotide non-randomness is a question of selection for certain optimal codons.     

Homology Requirements for Unequal Crossing over in Humans

To gain insight into mechanisms of unequal homologous recombination in vivo, genes generated by homologous unequal crossovers in the human beta-globin gene cluster were examined by nucleotide sequencing and hybridization experiments. The naturally occurring genes studied included one delta-beta Lepore-Baltimore fusion gene, one delta-beta Lepore-Hollandia fusion gene, 12 delta-beta Lepore-Boston genes, one A gamma-beta fusion Kenya gene, one A gamma-G gamma fusion (the central gene of a triplication) and one G gamma-A gamma fusion. A comparison of the nucleotide sequences of three Lepore-Boston genes indicates that they were derived from at least two independent homologous but unequal crossover events, although the crossovers occurred within the same 58-bp region. Nine additional Lepore-Boston genes from individuals of various ethnic origins were shown, by hybridization to specific oligonucleotide probes, to have been generated by a crossover in the same region as the sequenced genes. Evidence for gene conversion accompanying a homologous unequal crossover event was found in only one case (although some of the single nucleotide differences observed in other genes in this study may be related to the crossover events in ways that we do not presently understand). Thus, as judged by this limited sample, concurrent gene conversions are not commonly associated with homologous but unequal exchange in humans in vivo. Classification of the recombinant chromosomes by their polymorphic restriction sites in the beta-globin gene cluster indicated that the Lepore-Boston genes are found in at least six different haplotype backgrounds. Therefore the total number of independent examples in this study is at least 6, and at most 12. We have shown that in at least six cases of genes that have arisen by homologous but unequal crossing over in vivo, each event occurred in a relatively extensive region of uninterrupted identity between the parental genes. This preference cannot be explained by a mechanism whereby crossovers occur at random within misaligned related but not identical genes. In general, crossovers occur in regions that are among the largest available stretches of identity for a particular pair of mismatched genes. Our data are in agreement with those of other types of studies of homologous recombination, and support the idea that sequence identity, rather than general homology, is a critical factor in homologous recombination.     

Nucleotide sequence of gene pfkB encoding the minor phosphofructokinase of Escherichia coli K-12

The nucleotide sequence of a 1.3-kb DNA fragment containing the entire pfkB gene which codes for Pfk-2 of Escherichia coli, a minor phosphofructokinase (Pfk) enzyme, is reported. The Pfk-2 protein subunit is encoded by 924 bp, has 308 amino acids and an Mr of 33 000. Like other weakly expressed E. coli genes the codon usage in the pfkB gene is random; there is no strong bias for the usage of major tRNA isoaccepting species, and the codon preference rules of Grosjean and Fiers [Gene, 18 (1982) 199-209] are followed. This is the first report of the complete gene sequence of a phosphofructokinase.     

Molecular evolution of mammalian lactate dehydrogenase-A genes and pseudogenes: association of a mouse processed pseudogene with a B1 repetitive sequence

A mouse genomic clone containing a lactate dehydrogenase-A (LDH-A) processed pseudogene and a B1 repetitive element was isolated, and a nucleotide sequence of approximately 3 kb was determined. The pseudogene and B1 element are flanked by perfect 13-bp repeats, and the B1 sequence starts at 14 nucleotides 3' to the presumptive polyadenylation signal of the pseudogene. The nucleotide sequences of the LDH-A genes and processed pseudogenes from mouse, rat, and human were compared, and a phylogenetic tree was constructed. The rate and pattern of nucleotide substitutions in the LDH-A pseudogenes are similar to previously reported results (Li et al. 1984). The average rate of nucleotide substitutions in the LDH-A pseudogenes is 4.3 X 10(- 9)/site/year. The substitutions of C----T and G----A are most frequent, and A----G substitutions are relatively high. The rate of synonymous substitutions in the LDH-A genes is 5.3 X 10(-9), which is not significantly higher than the average rate of 4.7 X 10(-9) for 35 mammalian genes. The rate of nonsynonymous substitutions in the LDH-A genes is 0.20 X 10(-9), which is considerably lower than the average rate of 0.88 X 10(-9) for 35 mammalian genes. Thus, the mammalian LDH-A gene appears to be highly conserved in evolution.      

Novel third-letter bias in Escherichia coli codons revealed by rigorous treatment of coding constraints

A novel bias in codon third-letter usage was found in Escherichia coli genes with low fractions of "optimal codons", by comparing intact sequences with control random sequences. Third-letter usage has been found to be biased according to preference in codon usage and to doublet preference from the following first letter. The present study examines third-letter usage in the context of the nucleotide sequence when these preferences are considered. In order to exclude any influence by these factors, the random sequences were generated such that the amino acid sequence, codon usage, and the doublet frequency in each gene were all preserved. Comparison of intact sequences with these randomly generated sequences reveals that third letters of codons show a strong preference for the purine/pyrimidine pattern of the next codons: purine (R) is preferred to pyrimidine (Y) at the third site when followed by an R-Y-R codon, and pyrimidine is preferred when followed by an R-R-Y, an R-Y-Y or a Y-R-Y codon. This bias is probably related to interactions of tRNA molecules in the ribosome.     

Codon usage in Giardia lamblia.

A codon usage table for the intestinal parasite Giardia lamblia was generated by analysis of the nucleotide sequences of eight genes comprising 3,135 codons. Codon usage revealed a biased use of synonymous codons with a preference for NNC codons (42.1%). The codon usage of G. lamblia more closely resembles that of the prokaryote Halobacterium halobium (correlation coefficient r = 0.73) rather than that of other eukaryotic protozoans, i.e. Trypanosoma brucei (r = 0.434) and Plasmodium falciparum (r = -0.31). These observations are consistent with the view that G. lamblia represents the first line of descent from the ancestral cells that first took on eukaryotic features.