Identification and mapping of sigma-54 promoters in Chlamydia trachomatis

The first sigma(54) promoters in Chlamydia trachomatis L2 were mapped upstream of hypothetical proteins CT652.1 and CT683. Comparative genomics indicated that these sigma(54) promoters and potential upstream activation binding sites are conserved in orthologous C. trachomatis D, C. trachomatis mouse pneumonitis strain, and Chlamydia pneumoniae (CWL029 and AR39) genes.     

Amino Acid Biosynthetic Cost and Protein Conservation

Protein products of highly expressed genes tend to favor amino acids that have lower average biosynthetic costs (i.e., they exhibit metabolic efficiency). While this trend has been observed in several studies, the specific sites where cost-reducing substitutions accumulate have not been well characterized. Toward that end, weighted costs in conserved and variable positions were evaluated across a total of 9,119 homologous proteins in four mammalian orders (primate, carnivore, rodent, and artiodactyls), which together contain a total of 20,457,072 amino acids. Degree of conservation at homologous positions in these mammalian proteins and average-weighted cost across all positions within a single protein are significantly correlated. Dividing human genes into two classes (those with and those without CpG islands in their promoters) suggests that humans also preferentially utilize less costly amino acids in highly expressed genes. In contrast to the intuitive expectation that the relatively weak selective force associated with metabolic efficiency would be a selection pressure in complex multicellular organisms, the overall level of selective constraint within the variable regions of mammalian proteins allows the metabolic efficiency to derive a reduction of overall biosynthetic cost, particularly in genes with the highest levels of expression.     

Computational analysis of the polymorphic membrane protein superfamily of Chlamydia trachomatis and Chlamydia pneumoniae

Whole sequence genome analysis is invaluable in providing complete profiles of related proteins and gene families. The genome sequences of the obligate intracellular bacteria Chlamydia trachomatis and Chlamydia pneumoniae both encode proteins with similarity to several 90-kDa Chlamydia psittaci proteins. These proteins are members of a large superfamily, C. trachomatis with 9 members and C. pneumoniae with 21 members. All polymorphic membrane protein (Pmp) are heterogeneous, both in amino acid sequence and in predicted size. Most proteins have apparent signal peptide leader sequences and hence are predicted to be localized to the outer membrane. The unifying features of all proteins are the conserved amino acid motifs GGAI and FXXN repeated in the N-terminal half of each protein. In both genomes, the pmp genes are clustered at various locations on the chromosome. Phylogenetic analysis suggests six related families, each with at least one C. trachomatis and one C. pneumoniae orthologue. One of these families has seen prolific expansion in C. pneumoniae, resulting in 13 protein paralogues. The maintenance of orthologues from each species suggests specific functions for the proteins in chlamydial biology.     

Molecular analysis of two fructokinases involved in sucrose metabolism of enteric bacteria

Sucrose-positive derivatives of Escherichia coli K-12, containing the plasmid pUR400, and of Klebsiella pneumoniae hydrolyse intracellular sucrose 6-phosphate by means of an invertase into D-glucose 6-phosphate and free D-fructose. The latter is phosphorylated by an ATP-dependent fructokinase (gene scrK of an scr regulon) to D-fructose 6-phosphate. The lack of ScrK does not cause any visible phenotype in wild-type strains of both organisms. Using genes and enzymes normally involved in D-arabinitol metabolism from E. coli C and K. pneumoniae, derivatives of E. coli K-12 were constructed which allowed the identification of scrK mutations on conventional indicator plates. Cloning and sequencing of scrK from sucrose plasmid pUR400 and from the chromosome of K. pneumoniae revealed an open reading frame of 924 bp in both cases--the equivalent of a peptide containing 307 amino acid residues (Mr 39 and 34 kDa, respectively, on sodium dodecyl sulphate gels). The sequences showed overall identity among each other (69% identical residues) and to a kinase from Vibrio alginolyticus (57%) also involved in sucrose metabolism, lower overall identity (39%) to a D-ribose-kinase from E. coli, and local similarity to prokaryotic, and eukaryotic phosphofructokinases at the putative ATP-binding sites.     

Comparative analysis of amino acid usage and protein length distribution between alternatively and non-alternatively spliced genes across six eukaryotic genomes

Alternative splicing has been discovered in nearly all metazoan organisms as a mechanism to increase the diversity of gene products. However, the origin and evolution of alternatively spliced genes are still poorly understood. To understand the mechanisms for the evolution of alternatively spliced genes, it may be important to study the differences between alternatively and non-alternatively spliced genes. The aim of this research was to compare amino acid usage and protein length distribution between alternatively and non-alternatively spliced genes across six nearly complete eukaryotic genomes, including those of human (Homo sapiens), mouse (Mus musculus), rat (Rattus norvegicus), fruit fly (Drosophila melanogaster), Caenorhabditis elegans, and bovine (Bos taurus). Our results have suggested the following: (1) across the six species, alternatively and non-alternatively spliced genes have very similar tendency for amino acids usage for not only the overall scale but also those highly expressed genes, with all of the highly expressed genes having preferred amino acids including A, E, G, K, L, P, S, V, R, T, and D. (2) For not only the overall genes but also those highly expressed ones, the average length of the protein products of alternatively spliced genes is significantly greater than that of non-alternatively spliced ones. In contrast, distributions of protein lengths for the two groups of genes are very similar among all six species. Based on these results, we propose that alternatively spliced genes may have originated from non-alternatively spliced ones through events such as DNA mutations or gene fusion.     

Reinvestigating the codon and amino acid usage of S. cerevisiae genome: a new insight from protein secondary structure analysis

Biased usage of synonymous codons has been elucidated under the perspective of cellular tRNA abundance for quite a long time now. Taking advantage of publicly available gene expression data for Saccharomyces cerevisiae, a systematic analysis of the codon and amino acid usages in two different coding regions corresponding to the regular (helix and strand) as well as the irregular (coil) protein secondary structures, have been performed. Our analyses suggest that apart from tRNA abundance, mRNA folding stability is another major evolutionary force in shaping the codon and amino acid usage differences between the highly and lowly expressed genes in S. cerevisiae genome and surprisingly it depends on the coding regions corresponding to the secondary structures of the encoded proteins. This is obviously a new paradigm in understanding the codon usage in S. cerevisiae. Differential amino acid usage between highly and lowly expressed genes in the regions coding for the irregular protein secondary structure in S. cerevisiae is expounded by the stability of the mRNA folded structure. Irrespective of the protein secondary structural type, the highly expressed genes always tend to encode cheaper amino acids in order to reduce the overall biosynthetic cost of production of the corresponding protein. This study supports the hypothesis that the tRNA abundance is a consequence of and not a reason for the biased usage of amino acid between highly and lowly expressed genes.     

Do Amino Acid Biosynthetic Costs Constrain Protein Evolution in <Emphasis Type="Italic">Saccharomyces cerevisiae</Emphasis>?

Prokaryotic organisms preferentially utilize less energetically costly amino acids in highly expressed genes. Studies have shown that the proteome of Saccharomyces cerevisiae also exhibits this behavior, but only in broad terms. This study examines the question of metabolic efficiency as a proteome-shaping force at a finer scale, examining whether trends consistent with cost minimization as an evolutionary force are present independent of protein function and amino acid physicochemical property, and consistently with respect to amino acid biosynthetic costs. Inverse correlations between the average amino acid biosynthetic cost of the protein product and the levels of gene expression in S. cerevisiae are consistent with natural selection to minimize costs. There are, however, patterns of amino acid usage that raise questions about the strength (and possibly the universality) of this selective force in shaping S. cerevisiae’s proteome.     

Interaction between components of the type III secretion system of Chlamydiaceae

Members of the family Chlamydiaceae possess at least 13 genes, distributed throughout the chromosome, that are homologous with genes of known type III secretion systems (TTS). The aim of this study was to use putative TTS proteins of Chlamydophila pneumoniae, whose equivalents in other bacterial TTS function as chaperones, to identify interactions between chlamydial proteins. Using the BacterioMatch Two-Hybrid Vector system (Stratagene, La Jolla, Calif.), lcrH-2 and sycE, positions 1021 and 0325, respectively, from C. pneumoniae CM-1 were used as "bait" to identify target genes (positions 0324, 0705, 0708, 0808 to 0810, 1016 to 1020, and 1022) in close proximity on the chromosome. Interaction between the products of the lcrH-2 (1021) and lcrE (copN) (0324) genes was detected and confirmed by pull-down experiments and enzyme immunoassays using recombinant LcrH-2 and LcrE. As further confirmation of this interaction, the homologous genes from Chlamydia trachomatis, serovar E, and Chlamydophila psittaci, Texas turkey, were also cloned in the two-hybrid system to determine if LcrH-2 and LcrE would interact with their orthologs in other species. Consistent with their genetic relatedness, LcrH-2 from C. pneumoniae interacted with LcrE produced from the three species of Chlamydiaceae; LcrH-2 from C. psittaci reacted with LcrE from C. pneumoniae but not from C. trachomatis; and C. trachomatis LcrH-2 did not react with LcrE from the other two species. Deletions from the N and C termini of LcrE from C. pneumoniae identified the 50 C-terminal amino acids as essential for the interaction with LcrH-2. Thus, it appears that in the Chlamydiaceae TTS, LcrH-2 interacts with LcrE, and therefore it may serve as a chaperone for this protein.     

The major outer-membrane proteins of Chlamydia trachomatis serovars A and B: intra-serovar amino acid changes do not alter specificities of serovar- and C subspecies-reactive antibody-binding domains.

The major outer-membrane protein (MOMP) of Chlamydia trachomatis is a promising candidate antigen for chlamydial vaccine development. We have sequenced the MOMP genes for a serovar A and a serovar B isolate and have compared these new sequences with those already reported. Intra-serovar changes in the inferred amino acid sequences of the surface-exposed variable segments known to be responsible for binding of neutralizing antibody were observed. Nevertheless, epitope mapping with solid-phase peptides showed that these intra-serovar changes did not affect the binding of serovar- and subspecies-specific, potentially protective antibodies. Variable segment 1 of C. trachomatis serovar A contained two adjacent antibody-binding sites, one of which was C-subspecies specific while the other was serovar A specific. Therefore the subspecies binding site for C-complex organisms is in variable segment 1, whilst that for B-complex organisms is in variable segment 4. This work shows that MOMP sequences are relatively stable within the serovar categorization for isolates taken decades apart from different continents. Within a given serovar, however, limited interchange of functionally related amino acids may occur without impairing the binding of serovar-specific antibody.     

Bovine interferon alpha genes. Structure and expression

The bovine genome contains a gene family of interferon-alpha s (bIFN-alpha) that consists of at least five distinct members. Four of the bIFN-alpha genes isolated show a high degree of homology (97% in the nucleotide sequence and 93% in amino acid sequence). The overall homology in amino acid sequence of bIFN-alpha to human, murine, and rat IFNs-alpha is approximately 60%. Yet there are amino acid clusters (positions 28-41 and 118-146) which are highly conserved throughout the mammalian evolution and in which the overall homology can be as high as 86%. Within the C terminus conserved cluster there is a sequence containing 9 amino acids completely conserved in 16 mammalian IFNs-alpha and of these, 7 are also shared with a similar domain in some bacterial toxins, implying a common functional role for these domains. One of the genes, IFN-alpha C, was expressed in Escherichia coli. The purified bacterial IFN (specific activity, 2 X 10(8) units/mg) exhibited antiviral activity on bovine cells but no detectable activity was demonstrated on human and simian cells.     

Selection on Synthesis Cost Affects Interprotein Amino Acid Usage in All Three Domains of Life

Most investigations of the forces shaping protein evolution have focussed on protein function. However, cells are typically 50%–75% protein by dry weight, with protein expression levels distributed over five orders of magnitude. Cells may, therefore, be under considerable selection pressure to incorporate amino acids that are cheap to synthesize into proteins that are highly expressed. Such selection pressure has been demonstrated to alter amino acid usage in a few organisms, but whether “cost selection” is a general phenomenon remains unknown. One reason for this is that reliable protein expression level data is not available for most organisms. Accordingly, I have developed a new method for detecting cost selection. This method depends solely on interprotein gradients in amino acid usage. Applying it to an analysis of 43 whole genomes from all three domains of life, I show that selection on the synthesis cost of amino acids is a pervasive force in shaping the composition of proteins. Moreover, some amino acids have different price tags for different organisms—the cost of amino acids is changed for organisms living in hydrothermal vents compared with those living at the sea surface or for organisms that have difficulty acquiring elements such as nitrogen compared with those that do not—so I also investigated whether differences between organisms in amino acid usage might reflect differences in synthesis or acquisition costs. The results suggest that organisms evolve to alter amino acid usage in response to environmental conditions. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users. [Reviewing Editor: Hector Musto]     

A secondary structure motif predictive of protein localization to the chlamydial inclusion membrane

Chlamydiae are obligate intracellular pathogens that spend their entire growth phase sequestered in a membrane-bound vacuole called an inclusion. A set of chlamydial proteins, labelled Inc proteins, has been identified in the inclusion membrane (IM). The predicted IncA, IncB and IncC amino acid sequences share very limited similarity, but a common hydrophobicity motif is present within each Inc protein. In an effort to identify a relatively complete catalogue of Chlamydia trachomatis proteins present in the IM of infected cells, we have screened the genome for open reading frames encoding this structural motif. Hydropathy plot analysis was used to screen each translated open reading frame in the C. trachomatis genome database. Forty-six candidate IM proteins (C-lncs) that satisfied the criteria of containing a bilobed hydrophobic domain of at least 50 amino acids were identified. The genome of Chlamydia pneumoniae encodes a larger collection of C-lnc proteins, and only approximately half of the C-lncs are encoded within both genomes. In order to confirm the hydropathy plot screening method as a valid predictor of C-lncs, antisera and/or monoclonal antibodies were prepared against six of the C. trachomatis C-lncs. Immunofluorescence microscopy of C. trachomatis-infected cells probed with these antibodies showed that five out of six C-lncs are present in the chlamydial IM. Antisera were also produced against C. pneumoniae p186, a protein sharing identity with Chlamydia psittaci lncA and carrying a similar bilobed hydrophobic domain. These antisera labelled the inclusion membrane in C. pneumoniae infected cells, confirming that proteins sharing the unique secondary structural characteristic also localize to the inclusion membrane of C. pneumoniae. Sera from patients with high-titre antibodies to C. trachomatis were examined for reactivity with each tested C-lnc protein. Three out of six tested C-lncs were recognized by a majority of these patient sera. Collectively, these studies identify and characterize novel proteins localized to the chlamydial IM and demonstrate the existence of a potential secondary structural targeting motif for localization of chlamydial proteins to this unique intracellular environment.     

Metabolic and Translational Efficiency in Microbial Organisms

Metabolic efficiency, as a selective force shaping proteomes, has been shown to exist in Escherichia coli and Bacillus subtilis and in a small number of organisms with photoautotrophic and thermophilic lifestyles. Earlier attempts at larger-scale analyses have utilized proxies (such as molecular weight) for biosynthetic cost, and did not consider lifestyle or auxotrophy. This study extends the analysis to all currently sequenced microbial organisms that are amenable to these analyses while utilizing lifestyle specific amino acid biosynthesis pathways (where possible) to determine protein production costs and compensating for auxotrophy. The tendency for highly expressed proteins (with adherence to codon usage bias as a proxy for expressivity) to utilize less biosynthetically expensive amino acids is taken as evidence of cost selection. A comprehensive analysis of sequenced genomes to identify those that exhibit strong translational efficiency bias (389 out of 1,700 sequenced organisms) is also presented.     

Chlamydia species-dependent differences in the growth requirement for lysosomes

Genome reduction is a hallmark of obligate intracellular pathogens such as Chlamydia, where adaptation to intracellular growth has resulted in the elimination of genes encoding biosynthetic enzymes. Accordingly, chlamydiae rely heavily on the host cell for nutrients yet their specific source is unclear. Interestingly, chlamydiae grow within a pathogen-defined vacuole that is in close apposition to lysosomes. Metabolically-labeled uninfected host cell proteins were provided as an exogenous nutrient source to chlamydiae-infected cells, and uptake and subsequent labeling of chlamydiae suggested lysosomal degradation as a source of amino acids for the pathogen. Indeed, Bafilomycin A1 (BafA1), an inhibitor of the vacuolar H(+)/ATPase that blocks lysosomal acidification and functions, impairs the growth of C. trachomatis and C. pneumoniae, and these effects are especially profound in C. pneumoniae. BafA1 induced the marked accumulation of material within the lysosomal lumen, which was due to the inhibition of proteolytic activities, and this response inhibits chlamydiae rather than changes in lysosomal acidification per se, as cathepsin inhibitors also inhibit the growth of chlamydiae. Finally, the addition of cycloheximide, an inhibitor of eukaryotic protein synthesis, compromises the ability of lysosomal inhibitors to block chlamydial growth, suggesting chlamydiae directly access free amino acids in the host cytosol as a preferred source of these nutrients. Thus, chlamydiae co-opt the functions of lysosomes to acquire essential amino acids.     

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     

Gene expression level shapes the amino acid usages in Prochlorococcus marinus MED4

Prochlorococcus species are the first example of free-living bacteria with reduced genome. Codon and amino acid usages bias of Prochlorococcus marinus MED4 was investigated using all protein coding genes having length greater than or equal to 100 amino acids. Correspondence analysis on relative synonymous codon usage (RSCU) values shows that there is no such influence of translational selection in shaping the codon usage variation among the genes in this organism. However, amino acid usages were markedly different between the highly and lowly expressed genes in this organism and in particular, GC rich amino acids were found to occur significantly higher in highly expressed genes than the lowly expressed genes. Comparative analysis of the homologous genes of Synechococcus sp. WH8102 and Prochlorococcus marinus MED4 shows that amino acids conservation in highly expressed genes is significantly higher than lowly expressed genes. Based on our results we concluded that conservation of GC rich amino acids in the highly expressed genes to its ancestor is the major source of variation in amino acid usages in the organism.