Nucleotide sequence analysis of the phosphomannose isomerase gene (pmi) of Pseudomonas aeruginosa and comparison with the corresponding Escherichia coli gene manA

Phosphomannose isomerase (PMI) has been proposed to catalyze the first step of the alginic acid biosynthetic pathway in Pseudomonas aeruginosa. The nucleotide sequence of the P. aeruginosa pmi gene contained on a 2.0-kb BamHI-SstI DNA fragment has been determined. The gene was defined by the start and stop codons and by in vitro disruption of an open reading frame of 1440 bp corresponding to a polypeptide product with a predicted Mr of 52 860. This polypeptide displayed an apparent Mr of approx. 56 000 upon electrophoresis of a maxicell extract on sodium dodecyl sulfate-polyacrylamide gels. The codon utilization of the pmi gene was distinct in the wobble base preference and influenced by the high G + C content (66 mol%) of the P. aeruginosa DNA. Computer assisted matching analysis failed to demonstrate any significant homology at the nucleotide level between the P. aeruginosa pmi and Escherichia coli manA (pmi) genes. However, sequences homologous to the P. aeruginosa pmi gene were found in other Pseudomonas species, such as P. putida and P. mendocina, and in Azotobacter vinelandii, all capable of producing alginic acid.     

Real-time PCR detection of Pseudomonas aeruginosa in clinical and municipal wastewater and genotyping of the ciprofloxacin-resistant isolates

Real-time quantification of Pseudomonas aeruginosa was performed in various wastewater systems including clinical, municipal wastewaters and inflow from a wastewater treatment plant. The highest concentrations of P. aeruginosa-specific targets were detected in clinical wastewaters. Limitations of the detection system resulting from inhibition or cross-reaction were identified. Ciprofloxacin-resistant P. aeruginosa strains were isolated after specific enrichment from clinical and municipal wastewaters. In some cases they were also cultivated from effluent of a wastewater treatment plant, and from its downstream river water. A total of 119 isolates were phenotypically characterized as ciprofloxacin-resistant via antibiogram testing. Subsequently, the fluoroquinolone-resistance-mediating mutations in the genes gyrA codon positions 83 and 87, gyrB codon position 466 and parC codon positions 87 and 91 were determined by mini-sequencing. Ciprofloxacin resistance was mainly associated with mutations in gyrA codon position 83 and parC mutation in codon positions 87 or 91 of the bacterial gyrase and topoisomerase II genes. All ciprofloxacin-resistant P. aeruginosa strains were compared with genotypes from clinical data of fluoroquinolone-resistant P. aeruginosa infections. The results were in agreement with data from clinical analyses, with the exception that no gyrA 87 and no gyrB mutations were found in ciprofloxacin-resistant P. aeruginosa wastewater isolates.     

Gene expression and protein length influence codon usage and rates of sequence evolution in Populus tremula

Codon bias is generally thought to be determined by a balance between mutation, genetic drift, and natural selection on translational efficiency. However, natural selection on codon usage is considered to be a weak evolutionary force and selection on codon usage is expected to be strongest in species with large effective population sizes. In this paper, I study associations between codon usage, gene expression, and molecular evolution at synonymous and nonsynonymous sites in the long-lived, woody perennial plant Populus tremula (Salicaceae). Using expression data for 558 genes derived from expressed sequence tags (EST) libraries from 19 different tissues and developmental stages, I study how gene expression levels within single tissues as well as across tissues affect codon usage and rates sequence evolution at synonymous and nonsynonymous sites. I show that gene expression have direct effects on both codon usage and the level of selective constraint of proteins in P. tremula, although in different ways. Codon usage genes is primarily determined by how highly expressed a genes is, whereas rates of sequence evolution are primarily determined by how widely expressed genes are. In addition to the effects of gene expression, protein length appear to be an important factor influencing virtually all aspects of molecular evolution in P. tremula.     

Inferring gene function from evolutionary change in signatures of translation efficiency

Background

The genetic code is redundant, meaning that most amino acids can be encoded by more than one codon. Highly expressed genes tend to use optimal codons to increase the accuracy and speed of translation. Thus, codon usage biases provide a signature of the relative expression levels of genes, which can, uniquely, be quantified across the domains of life.

Results

Here we describe a general statistical framework to exploit this phenomenon and to systematically associate genes with environments and phenotypic traits through changes in codon adaptation. By inferring evolutionary signatures of translation efficiency in 911 bacterial and archaeal genomes while controlling for confounding effects of phylogeny and inter-correlated phenotypes, we linked 187 gene families to 24 diverse phenotypic traits. A series of experiments in Escherichia coli revealed that 13 of 15, 19 of 23, and 3 of 6 gene families with changes in codon adaptation in aerotolerant, thermophilic, or halophilic microbes. Respectively, confer specific resistance to, respectively, hydrogen peroxide, heat, and high salinity. Further, we demonstrate experimentally that changes in codon optimality alone are sufficient to enhance stress resistance. Finally, we present evidence that multiple genes with altered codon optimality in aerobes confer oxidative stress resistance by controlling the levels of iron and NAD(P)H.

Conclusions

Taken together, these results provide experimental evidence for a widespread connection between changes in translation efficiency and phenotypic adaptation. As the number of sequenced genomes increases, this novel genomic context method for linking genes to phenotypes based on sequence alone will become increasingly useful.     

The complete nucleotide sequence of phi CTX, a cytotoxin-converting phage of Pseudomonas aeruginosa: implications for phage evolution and horizontal gene transfer via bacteriophages

phi CTX is a cytotoxin-converting phage isolated from Pseudomonas aeruginosa. In this study, we determined the complete nucleotide sequence of the phi CTX phage genome. The precise genome size was 35,538 bp with 21 base 5'-extruding cohesive ends. Forty-seven open reading frames (ORFs) were identified on the phi CTX genome, including two previously identified genes, ctx and int. Among them, 15 gene products were identified in the phage particle by protein microsequencing. The most striking feature of the phi CTX genome was an extensive homology with the coliphage P2 and P2-related phages; more than half of the ORFs (25 ORFs) had marked homology to P2 genes with 28.9-65.8% identity. The gene arrangement on the genome was also highly conserved for the two phages, although the G + C content and codon usage of most phi CTX genes were similar to those of the host P. aeruginosa chromosome. In addition, phi CTX was found to share several common features with P2, including the morphology, non-inducibility, use of lipopolysaccharide core oligosaccharide as receptor and Ca(2+)-dependent receptor binding. These findings indicate that phi CTX is a P2-like phage well adapted to P. aeruginosa, and provide clear evidence of the intergeneric spread and evolution of bacteriophages. Furthermore, comparative analysis of genome structures of phi CTX, P2 and other P2 relatives revealed the presence of several hot-spots where foreign DNAs, including the cytotoxin gene, were inserted. They appear to be deeply concerned in the acquisition of various genes that are horizontally transferred by bacteriophage infection.     

Selection on the codon bias ofChlamydomonas reinhardtii chloroplast genes and the plantpsbA gene

Plant chloroplast genes have a codon use that reflects the genome compositional bias of a high A+T content with the single exception of the highly translatedpsbA gene which codes for the photosystem II D1 protein. The codon usage of plantpsbA corresponds more closely to the limited tRNA population of the chloroplast and is very similar to the codon use observed in the chloroplast genes of the green algaChlamydomonas reinhardtii. This pattern of codon use may be an adaptation for increased translation efficiency. A correspondence between codon use of plantpsbA andChlamydomonas chloroplast genes and the tRNAs coded by the chloroplast genome, however, is not observed in all synonymous codon groups. It is shown here that the degree of correspondence between codon use and tRNA population in different synonymous groups is correlated with the second codon position composition. Synonymous groups with an A or T at the second codon position have a high representation of codons for which a complementary tRNA is coded by the chloroplast genome. Those with a G or C at the second position have an increased representation of codons that bind a chloroplast tRNA by wobble. It is proposed that the difference between synonymous groups in terms of codon adaptation to the tRNA population in plantpsbA andChlamydomonas chloroplast genes may be the result of differences in second position composition.     

Glutathione S-transferase GSTM1, GSTT1 and p53 codon 72 polymorphisms in human tumor cells

The genes of the glutathione S-transferase (GST) family encode enzymes that appear to be critical in cellular protection against the cytotoxic effects, whereas p53 is a tumor suppressor gene. Despite a large number of studies on germline polymorphisms of GSTM1, GSTT1 and p53 genes, there have been very few reports on genotyping of these genes in human malignant tumor cells. In this study, we investigated GSTM1, GSTT1 and p53 codon 72 polymorphisms in a variety of human tumor cell lines originating from different organs to clarify tissue-specific polymorphic frequency of these genes in human solid tumors. The GSTM1 and GSTT1 genetic polymorphisms were evaluated using multiplex PCR techniques and PCR-RFLP analysis was conducted to identify p53 codon 72 genotypes. Gene expression of GSTM1 or GSTT1 was detected by RT-PCR in the cells with respective present genotype for each. Polymorphisms of p53 codon 72 detected by PCR-RFLP were also confirmed using SSCP and sequence analyses. GSTM1 and GSTT1 genotypes were various in 104 cell lines examined. Null GSTM1 genotype was dominant in small cell lung, kidney and ovarian carcinoma cells, whereas null GSTT1 genotype was dominant in cervical and endometrial carcinoma cells. GSTM1 and GSTT1 genotypes in ovarian carcinoma cells were quite similar to those in small cell lung carcinoma cells. Polymorphic frequency of p53 codon 72 was also various among the cells, however, the Pro allele was found in only 1 of 6 kidney, 14 cervical and 4 endometrial carcinoma cell lines. There was a significant difference in GSTM1 and p53 genotypes between 34 small cell and 24 non small cell lung carcinoma cells (P < 0.01). Combined study on the distribution of GSTM1, GSTT1 and p53 genotypes revealed that null GSTM1 genotype was associated with the Arg allele of p53 codon 72 in 58 lung carcinoma cells and null GSTT1 genotype was associated with the Pro/Pro homozygote in 104 tumor cell lines examined. This is the first study examining GSTM1, GSTT1 and p53 codon 72 polymorphisms in a variety of human solid tumor cells and suggesting that polymorphic frequency of these genes may be tissue- and organ-specific. The molecular interaction between GST gene defects and p53 codon 72 genotype in the development of human malignant tumors should be further investigated.     

The genetics and genomics of the drought response in Populus

The genetic nature of tree adaptation to drought stress was examined by utilizing variation in the drought response of a full-sib second generation (F(2)) mapping population from a cross between Populus trichocarpa (93-968) and P. deltoides Bart (ILL-129) and known to be highly divergent for a vast range of phenotypic traits. We combined phenotyping, quantitative trait loci (QTL) analysis and microarray experiments to demonstrate that 'genetical genomics' can be used to provide information on adaptation at the species level. The grandparents and F(2) population were subjected to soil drying, and contrasting responses to drought across genotypes, including leaf coloration, expansion and abscission, were observed, and QTL for these traits mapped. A subset of extreme genotypes exhibiting extreme sensitivity and insensitivity to drought on the basis of leaf abscission were defined, and microarray experiments conducted on these genotypes and the grandparent species. The extreme genotype groups induced a different set of genes: 215 and 125 genes differed in their expression response between groups in control and drought, respectively, suggesting species adaptation at the gene expression level. Co-location of differentially expressed genes with drought-specific and drought-responsive QTLs was examined, and these may represent candidate genes contributing to the variation in drought response.     

Codon Usage in Plastid Genes Is Correlated with Context, Position Within the Gene, and Amino Acid Content

Highly expressed plastid genes display codon adaptation, which is defined as a bias toward a set of codons which are complementary to abundant tRNAs. This type of adaptation is similar to what is observed in highly expressed Escherichia coli genes and is probably the result of selection to increase translation efficiency. In the current work, the codon adaptation of plastid genes is studied with regard to three specific features that have been observed in E. coli and which may influence translation efficiency. These features are (1) a relatively low codon adaptation at the 5′ end of highly expressed genes, (2) an influence of neighboring codons on codon usage at a particular site (codon context), and (3) a correlation between the level of codon adaptation of a gene and its amino acid content. All three features are found in plastid genes. First, highly expressed plastid genes have a noticeable decrease in codon adaptation over the first 10–20 codons. Second, for the twofold degenerate NNY codon groups, highly expressed genes have an overall bias toward the NNC codon, but this is not observed when the 3′ neighboring base is a G. At these sites highly expressed genes are biased toward NNT instead of NNC. Third, plastid genes that have higher codon adaptations also tend to have an increased usage of amino acids with a high G + C content at the first two codon positions and GNN codons in particular. The correlation between codon adaptation and amino acid content exists separately for both cytosolic and membrane proteins and is not related to any obvious functional property. It is suggested that at certain sites selection discriminates between nonsynonymous codons based on translational, not functional, differences, with the result that the amino acid sequence of highly expressed proteins is partially influenced by selection for increased translation efficiency. Received: 21 July 1999 / Accepted: 5 November 1999     

Molecular analysis of the folC gene of Pseudomonas aeruginosa

We have cloned the Pseudomonas aeruginosa folC gene coding for folylpolyglutamate synthetase-dihydrofolate synthetase, which was located between the trpF and purF loci, and determined the nucleotide sequence of the folC gene and its flanking region. The deduced amino acid sequence of P. aeruginosa FolC was highly homologous to that of Escherichia coli FolC. The cloned gene complemented E. coli folC mutations and was found to encode both folylpolyglutamate synthetase and dihydrofolate synthetase activities. The gene organization around the folC gene in P. aeruginosa was completely conserved with that in E. coli; the accD gene was located upstream of the folC gene, and dedD, cvpA and purF genes followed the folC gene in this order. The gene arrangement and the result of the promoter activity assay suggested that the P. aeruginosa accD and folC genes were co-transcribed.     

Characterization of the GO system of Pseudomonas aeruginosa

The mutT, mutM, and mutY genes of the GO system of the Pseudomonas aeruginosa PAO1 strain have been characterized by cloning, sequencing, and complementation analysis. The three genes, when cloned in a plasmid, were able to complement the high mutation frequency of the corresponding Escherichia coli deficient strains. Our results demonstrate that the putative mutT, mutM, and mutY gene products from P. aeruginosa are able to perform the expected activity. In addition, the sequence of the P. aeruginosa mutT gene strongly suggested that the product of this gene has a bifunctional activity in P. aeruginosa, being the C-terminal part 40% identical to a consensus sequence of thiamine monophosphate synthases. Our results also demonstrated that the N-terminal part of the protein is necessary and sufficient for the 8-oxodGTP hydrolase activity.     

A novel type II secretion system in Pseudomonas aeruginosa

The genome sequence of Pseudomonas aeruginosa strain PAO1 has been determined to facilitate postgenomic studies aimed at understanding the capacity of adaptation of this ubiquitous opportunistic pathogen. P. aeruginosa produces toxins and hydrolytic enzymes that are secreted via the type II secretory pathway using the Xcp machinery or 'secreton'. In this study, we characterized a novel gene cluster, called hxc for homologous to xcp. Characterization of an hxcR mutant, grown in phosphate-limiting medium, revealed the absence of a 40 kDa protein found in the culture supernatant of wild-type or xcp derivative mutant strains. The protein corresponded to the alkaline phosphatase L-AP, renamed LapA, which is secreted in an xcp-independent but hxc-dependent manner. Finally, we showed that expression of the hxc gene cluster is under phosphate regulation. This is the first report of the existence of two functional type II secretory pathways within the same organism, which could be related to the high adaptation potential of P. aeruginosa.     

用图论方法研究核酸序列的密码子使用与基因表达水平的关系

研究了Escherichiacoli（115个基因）和SacharomycesYeast（97个基因）核酸序列的密码子使用频率与基因表达水平的关系．将同义密码子按使用频率统计值分成三种特性的密码子：最适密码子（H）、非最适密码子（L）和稀有密码子（R）,对每一基因序列的编码区,算出它们各自出现的概率P（H）,P（L）和P（R）．以P（H）和P（R）为指标,用图论法聚类,发现每种生物的高低表达基因明显分开,基因表达水平被分为四级：甚高表达基因（VH）、高表达基因（H）、较低表达基因（LM）和低表达基因（LL）．每类基因的表达水平与实验结果保持了很好的相关性,与E．coli和Yeast的现有资料相比,符合很好．     

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.