Distribution of repetitive sequences on the leading and lagging strands of the Escherichia coli genome: comparative study of Long Direct Repeat (LDR) sequences.

In the present study, we developed a method for detecting sequences whose similarity to a target sequence is statistically significant and we examined the distribution of these sequences in the E. coli K-12 genome. Target sequences examined are as follows: (i) short repeat: Crossover hot-spot instigator (Chi) sequence, replication termination (Ter) sequence, and DnaA binding sequence (DnaA box); (ii) potential stem-loop structure repeats: palindromic unit (PU), boxC sequences, and intergenic repeat unit (IRU); (iii) potential RNA coding repeats: rRNAs, PAIR, TRIP, and QUAD; and (iv) potential protein coding repeats: insertion elements (ISs) and Long Direct Repeats (LDRs). We also examined the distribution of these sequences on leading and lagging strands. We obtained another four statistically significant LDR sequences with more than 187 bp matched to LDR-A near the LDR loci, suggesting that these regions might be used as high recombination hot spots for LDR. Adaptation of individual LDRs to E. coli genome is also discussed on the basis of codon usage.     

Genome sequence of Yersinia pestis KIM

We present the complete genome sequence of Yersinia pestis KIM, the etiologic agent of bubonic and pneumonic plague. The strain KIM, biovar Mediaevalis, is associated with the second pandemic, including the Black Death. The 4.6-Mb genome encodes 4,198 open reading frames (ORFs). The origin, terminus, and most genes encoding DNA replication proteins are similar to those of Escherichia coli K-12. The KIM genome sequence was compared with that of Y. pestis CO92, biovar Orientalis, revealing homologous sequences but a remarkable amount of genome rearrangement for strains so closely related. The differences appear to result from multiple inversions of genome segments at insertion sequences, in a manner consistent with present knowledge of replication and recombination. There are few differences attributable to horizontal transfer. The KIM and E. coli K-12 genome proteins were also compared, exposing surprising amounts of locally colinear "backbone," or synteny, that is not discernible at the nucleotide level. Nearly 54% of KIM ORFs are significantly similar to K-12 proteins, with conserved housekeeping functions. However, a number of E. coli pathways and transport systems and at least one global regulator were not found, reflecting differences in lifestyle between them. In KIM-specific islands, new genes encode candidate pathogenicity proteins, including iron transport systems, putative adhesins, toxins, and fimbriae.     

Identification of the Chi site of Haemophilus influenzae as several sequences related to the Escherichia coli Chi site

The Escherichia coli Chi site 5'-GCTGGTGG-3' modulates the activity of the powerful dsDNA exonuclease and helicase RecBCD. Genome sequence analyses revealed that Chi is frequent on the chromosome and oriented with respect to replication on the E . coli genome. Chi is also present much more frequently than predicted statistically for a random 8-mer sequence. Although it is assumed that Chi is ubiquitous, there is virtually no proof that its features are conserved in other microorganisms. We therefore identified and analysed the Chi sequence of an organism for which the full genome sequence was available, Haemophilus influenzae . The biological test we used is based on our finding that rolling circle plasmids provide a specific substrate for RecBCD analogues in different microorganisms. Unexpectedly, several related sequences, corresponding to 5'-GNTGGTGG-3' and 5'-G(G/C)TGGAGG-3', showed Chi activity. As in E . coli , the H . influenzae Chi sites are frequent on the genome, which is in keeping with the need for frequent Chi sites for dsDNA break repair of chromosomal DNA. Although statistically over-represented, this feature is less marked than that of the E . coli Chi site. In contrast to E . coli , the H . influenzae Chi motifs are only slightly oriented with respect to the replication strand. Thus, although Chi appears to have a highly conserved biological role in attenuating exonuclease activity, its sequence characteristics and statistical representation on the genome may differ according to the particular features of the host.     

Unexpected DNA context-dependence identifies a new determinant of Chi recombination hotspots

Homologous recombination occurs especially frequently near special chromosomal sites called hotspots. In Escherichia coli, Chi hotspots control RecBCD enzyme, a protein machine essential for the major pathway of DNA break-repair and recombination. RecBCD generates recombinogenic single-stranded DNA ends by unwinding DNA and cutting it a few nucleotides to the 3′ side of 5′ GCTGGTGG 3′, the sequence historically equated with Chi. To test if sequence context affects Chi activity, we deep-sequenced the products of a DNA library containing 10 random base-pairs on each side of the Chi sequence and cut by purified RecBCD. We found strongly enhanced cutting at Chi with certain preferred sequences, such as A or G at nucleotides 4–7, on the 3′ flank of the Chi octamer. These sequences also strongly increased Chi hotspot activity in E. coli cells. Our combined enzymatic and genetic results redefine the Chi hotspot sequence, implicate the nuclease domain in Chi recognition, indicate that nicking of one strand at Chi is RecBCD''s biologically important reaction in living cells, and enable more precise analysis of Chi''s role in recombination and genome evolution.     

Activation of Chi recombinational hotspots by RecBCD-like enzymes from enteric bacteria

Chi sites, 5'G-C-T-G-G-T-G-G-3', enhance homologous recombination in Escherichia coli and are activated by the RecBCD enzyme. To test the ability of Chi to be activated by analogous enzymes from other bacteria, we cloned recBCD-like genes from diverse bacteria into an E. coli recBCD deletion mutant. Clones from seven species of enteric bacteria conferred to this deletion mutant recombination proficiency, Chi hotspot activity in lambda Red- Gam- vegetative crosses, and RecBCD enzyme activities, including Chi-dependent DNA strand cleavage. Three clones from Pseudomonas aeruginosa and Ps. putida conferred recombination proficiency and ATP-dependent nuclease activity, but neither Chi hotspot activity nor Chi-dependent DNA cleavage. These results imply that Chi has been conserved as a recombination-promoting signal for RecBCD-like enzymes in enteric bacteria but not in more distantly related bacteria such as Pseudomonas spp. We discuss the possibility that other, presently unknown, nucleotide sequences serve the same function as Chi in Pseudomonas spp.     

Nucleotide sequence of the structural gene for tryptophanase of Escherichia coli K-12.

The tryptophanase structural gene, tnaA, of Escherichia coli K-12 was cloned and sequenced. The size, amino acid composition, and sequence of the protein predicted from the nucleotide sequence agree with protein structure data previously acquired by others for the tryptophanase of E. coli B. Physiological data indicated that the region controlling expression of tnaA was present in the cloned segment. Sequence data suggested that a second structural gene of unknown function was located distal to tnaA and may be in the same operon. The pattern of codon usage in tnaA was intermediate between codon usage in four of the ribosomal protein structural genes and the structural genes for three of the tryptophan biosynthetic proteins.     

crp genes of Shigella flexneri, Salmonella typhimurium, and Escherichia coli.

The complete nucleotide sequences of the Salmonella typhimurium LT2 and Shigella flexneri 2B crp genes were determined and compared with those of the Escherichia coli K-12 crp gene. The Shigella flexneri gene was almost like the E. coli crp gene, with only four silent base pair changes. The S. typhimurium and E. coli crp genes presented a higher degree of divergence in their nucleotide sequence with 77 changes, but the corresponding amino acid sequences presented only one amino acid difference. The nucleotide sequences of the crp genes diverged to the same extent as in the other genes, trp, ompA, metJ, and araC, which are structural or regulatory genes. An analysis of the amino acid divergence, however, revealed that the catabolite gene activator protein, the crp gene product, is the most conserved protein observed so far. Comparison of codon usage in S. typhimurium and E. coli for all genes sequenced in both organisms showed that their patterns were similar. Comparison of the regulatory regions of the S. typhimurium and E. coli crp genes showed that the most conserved sequences were those known to be essential for the expression of E. coli crp.     

The recA-recBCD dependent recombination pathways of Serratia marcescens and Proteus mirabilis in Escherichia coli: functions of hybrid enzymes and hybrid pathways.

The physical maps of cloned recBCD gene regions of Serratia marcescens and Proteus mirabilis were correlated to genes located in this region. The genes thyA, recC, recB, recD and argA were organized as in Escherichia coli. The 3 rec genes code for the 3 different subunits of the RecBCD enzyme and produced enzymes promoting recombination and repair of UV damage in E coli. The recBCD-dependent stimulation of recombination at specific nucleotide sequences called Chi (Chi-activation) was determined in lambda red-gam-crosses. Chi-activation by the different RecBCD enzymes decreased in the order E coli greater than S marcescens greater than P mirabilis. In E coli cloned subunits genes from S marcescens and P mirabilis led to the formation of functional hybrid enzymes consisting of subunits from 2 or even 3 species. The origin of the RecC subunit present in the hybrid enzymes affected the degree of Chi-activation. Further, changes in Chi-activation occurred when the RecD subunit in the enzyme from E coli was replaced by RecD proteins from S marcescens or P mirabilis. This suggested that the RecD subunit determines not only whether or not Chi-activation is possible but also to which extent it occurs. Finally we have reconstituted recombination pathways of S marcescens and P mirabilis by combining the cloned recA and recBCD genes from these species in E coli deleted for recA and recBCD. Both pathways can efficiently promote recombination and repair. Studies are summarized which showed that levels of repair and recombination promoted by the recA-recBCD genes are mostly higher when the recA and recBCD genes came from the same species than from 2 different species (hybrid RecBCD recombination pathway). The data are interpreted to provide evidence that in vivo the RecA protein co-operates with the RecBCD enzyme in recombination and repair of UV damage.     

Analysis of genome plasticity in pathogenic and commensal Escherichia coli isolates by use of DNA arrays

Genomes of prokaryotes differ significantly in size and DNA composition. Escherichia coli is considered a model organism to analyze the processes involved in bacterial genome evolution, as the species comprises numerous pathogenic and commensal variants. Pathogenic and nonpathogenic E. coli strains differ in the presence and absence of additional DNA elements contributing to specific virulence traits and also in the presence and absence of additional genetic information. To analyze the genetic diversity of pathogenic and commensal E. coli isolates, a whole-genome approach was applied. Using DNA arrays, the presence of all translatable open reading frames (ORFs) of nonpathogenic E. coli K-12 strain MG1655 was investigated in 26 E. coli isolates, including various extraintestinal and intestinal pathogenic E. coli isolates, 3 pathogenicity island deletion mutants, and commensal and laboratory strains. Additionally, the presence of virulence-associated genes of E. coli was determined using a DNA "pathoarray" developed in our laboratory. The frequency and distributional pattern of genomic variations vary widely in different E. coli strains. Up to 10% of the E. coli K-12-specific ORFs were not detectable in the genomes of the different strains. DNA sequences described for extraintestinal or intestinal pathogenic E. coli are more frequently detectable in isolates of the same origin than in other pathotypes. Several genes coding for virulence or fitness factors are also present in commensal E. coli isolates. Based on these results, the conserved E. coli core genome is estimated to consist of at least 3,100 translatable ORFs. The absence of K-12-specific ORFs was detectable in all chromosomal regions. These data demonstrate the great genome heterogeneity and genetic diversity among E. coli strains and underline the fact that both the acquisition and deletion of DNA elements are important processes involved in the evolution of prokaryotes.     

Conservation of translation initiation sites based on dinucleotide frequency and codon usage in Escherichia coli K-12 (W3110): non-random distribution of A/T-rich sequences immediately upstream of the translation initiation codon.

Dinucleotide frequencies are useful for characterizing consensus elements as a minimum unit of nucleotide sequence because the neighborhood relations of nucleotide sequences are reflected in dinucleotides. Using a consensus score based on dinucleotide frequencies and intra-species codon usage heterogeneity, denoted by the Z1 parameter, we report the relationship between nucleotide conservation at the translation initiation sites of genes in the Escherichia coli K-12 genome (W3110) and codon usage in its downstream genes. Significant positive correlations were obtained in three regions centered at -13, -4, and +7, which correspond to the Shine-Dalgarno element, the A + T element immediately upstream of the translation initiation site, and the downstream box, respectively.     

Codon usage in bacteria: correlation with gene expressivity

The nucleic acid sequence bank now contains over 600 protein coding genes of which 107 are from prokaryotic organisms. Codon frequencies in each new prokaryotic gene are given. Analysis of genetic code usage in the 83 sequenced genes of the Escherichia coli genome (chromosome, transposons and plasmids) is presented, taking into account new data on gene expressivity and regulation as well as iso-tRNA specificity and cellular concentration. The codon composition of each gene is summarized using two indexes: one is based on the differential usage of iso-tRNA species during gene translation, the other on choice between Cytosine and Uracil for third base. A strong relationship between codon composition and mRNA expressivity is confirmed, even for genes transcribed in the same operon. The influence of codon use of peptide elongation rate and protein yield is discussed. Finally, the evolutionary aspect of codon selection in mRNA sequences is studied.     

Chi Enhances Heteroduplex DNA Levels during Recombination

The major pathway of homologous recombination in Escherichia coli, the RecBCD pathway, is stimulated by Chi sites. To determine whether Chi enhances an early or late step in recombination, we measured formation of heteroduplex DNA (hDNA) in extracts of lambda-infected E. coli. Chi elevated hDNA levels in these extracts, supporting a role for Chi early (before hDNA formation) in recombination. RecA protein and RecBCD enzyme were both necessary for detection of hDNA, indicating that they, too, act early. Analysis of a panel of recBCD mutants indicated that Chi-nicking activity was needed for Chi's stimulation of hDNA formation. These results support a previously proposed model of recombination. Further results suggested that RecBCD enzyme has an additional role late in recombination.     

Di-codon usage for classification of genes

Genes are often classified into biologically related groups so that inferences on their functions can be made. This paper demonstrates that the di-codon usage is a useful feature for gene classification and gives better classification accuracy than the codon usage. Our experiments with different classifiers show that support vector machines performs better than other classifiers in classifying genes by using di-codon usage as features. The method is illustrated on 1841 HLA sequences which are classified into two major classes, HLA-I and HLA-II, and further classified into the subclasses of major classes. By using both codon and di-codon features, we show near perfect accuracies in the classification of HLA molecules into major classes and their sub-classes.     

Salmonella recD mutations increase recombination in a short sequence transduction assay.

We have identified recD mutants of Salmonella typhimurium by their ability to support growth of phage P22 abc (anti-RecBCD) mutants, whose growth is prevented by normal host RecBCD function. As in Escherichia coli, the recD gene of S. typhimurium lies between the recB and argA genes at min 61 of the genetic map. Plasmids carrying the Salmonella recBCD+ genes restore ATP-dependent exonuclease V activity to an E. coli recBCD deletion mutant. The new Salmonella recD mutations (placed on this plasmid) eliminate the exonuclease activity and enable the plasmid-bearing E. coli deletion mutant to support growth of phage T4 gene 2 mutants. The Salmonella recD mutations caused a 3- to 61-fold increase in the ability of a recipient strain to inherit (by transduction) a large inserted element (MudA prophage; 38 kb). In this cross, recombination events must occur in the short (3-kb) sequences that flank the element in the 44-kb transduced fragment. The effect of the recD mutation depends on the nature of the flanking sequences and is likely to be greatest when those sequences lack a Chi site. The recD mutation appears to minimize fragment degradation and/or cause RecBC-dependent recombination events to occur closer to the ends of the transduced fragment. The effect of a recipient recD mutation was eliminated if the donor P22 phage expressed its Abc (anti-RecBC) function. We hypothesize that in standard (high multiplicity of infection) P22-mediated transduction crosses, recombination is stimulated both by Chi sequences (when present in the transduced fragment) and by the phage-encoded Abc protein which inhibits the host RecBCD exonuclease.     

沙枣叶绿体全基因组序列及其使用密码子偏性分析

该研究以2株野生沙枣(Elaeagnus angustifolia Linn.)嫩枝经温室水培后的嫩叶为材料,采用CTAB法分别提取总DNA,并利用第二代测序技术进行总DNA从头测序,组装后得到2株沙枣叶绿体基因组全序列,并详细分析了其蛋白质编码基因密码子使用的偏好性及其原因,为沙枣叶绿体基因工程和分子系统进化等研究奠定基础。结果显示:(1)组装得到沙枣叶绿体基因组序列全长150 546 bp,由长度为81 113 bp的长单拷贝(LSC)区域和25 494 bp的短单拷贝(SSC)区域,以及1对分隔开它们的长18 445 bp的反向重复序列(IRS)组成;注释共得到132个基因,包括86个蛋白编码基因、38个tRNA基因和8个rRNA基因。(2)沙枣叶绿体基因组蛋白编码基因密码子的第三位碱基GC含量(GC_3)为28.47%,明显低于整个叶绿体基因组GC含量(37%),也低于第一位(GC_1)和第二位(GC_2)碱基的GC含量,说明密码子对AT碱基结尾有偏好性;其中, UCU、CCU、UGU、GCU、CUU、GAU、UCA和UAA为最优密码子。(3)同义密码子相对使用频率(RSCU)分析发现,影响密码子使用模式的因素并不单一,密码子的偏好性受到突变、选择及其他因素的共同影响,并且自然选择表达引起的序列差异比突变对密码子偏好性的影响要显著;中性绘图分析、有效密码子数(ENC-plot)分析和奇偶偏好性(PR2-plot)分析表明,沙枣叶绿体基因组使用密码子的偏性受选择的影响更大。(4)通过最大似然法、最大简约法和贝叶斯方法对胡颓子科6个物种和1个枣的叶绿体基因序列构建系统发育树,与它们使用密码子偏性聚类的结果一致,表明叶绿体基因组使用密码子偏性与物种的亲缘关系相关。     

杜仲基因密码子使用模式分析

为了解杜仲基因密码子使用模式,该文以杜仲基因组密码子为研究对象,运用CodonW软件对杜仲的320个蛋白编码基因进行同义密码子相对使用频率(RSCU)分析、ENC-GC3s关联分析编码基因的密码子ENC值、PR2-plot偏倚分析编码基因的密码子碱基使用频率,并运用CUSP软件与Codon Usage Database...     

DNA sequence analysis of the dnaK gene of Escherichia coli B and of two dnaK genes carrying the temperature-sensitive mutations dnaK7(Ts) and dnaK756(Ts).

The DNA sequence of the dnaK gene of Escherichia coli was analyzed. The nucleotide sequence of the wild-type dnaK gene of E. coli B differed from that of E. coli K-12 in 15 bp, none of which altered the amino acid sequence. Two temperature-sensitive dnaK mutations were examined by cloning and sequence analyses. Results showed that one dnaK mutation, dnaK7(Ts), was a one-base substitution of T for C at nucleotide position 448 in the open reading frame yielding an amber nonsense codon. The other mutation, dnaK756(Ts), consisted of base substitutions (A for G) at three nucleotide positions, 95, 1364, and 1403, in the open reading frame resulting in an aspartic acid codon in place of a glycine codon.     

The cobalamin (coenzyme B12) biosynthetic genes of Escherichia coli.

The enteric bacterium Escherichia coli synthesizes cobalamin (coenzyme B12) only when provided with the complex intermediate cobinamide. Three cobalamin biosynthetic genes have been cloned from Escherichia coli K-12, and their nucleotide sequences have been determined. The three genes form an operon (cob) under the control of several promoters and are induced by cobinamide, a precursor of cobalamin. The cob operon of E. coli comprises the cobU gene, encoding the bifunctional cobinamide kinase-guanylyltransferase; the cobS gene, encoding cobalamin synthetase; and the cobT gene, encoding dimethylbenzimidazole phosphoribosyltransferase. The physiological roles of these sequences were verified by the isolation of Tn10 insertion mutations in the cobS and cobT genes. All genes were named after their Salmonella typhimurium homologs and are located at the corresponding positions on the E. coli genetic map. Although the nucleotide sequences of the Salmonella cob genes and the E. coli cob genes are homologous, they are too divergent to have been derived from an operon present in their most recent common ancestor. On the basis of comparisons of G+C content, codon usage bias, dinucleotide frequencies, and patterns of synonymous and nonsynonymous substitutions, we conclude that the cob operon was introduced into the Salmonella genome from an exogenous source. The cob operon of E. coli may be related to cobalamin synthetic genes now found among non-Salmonella enteric bacteria.