Homorepeat variability within the human population.

Genetic variation within populations plays a crucial role in driving evolution. Unlike the average protein sequence, the evolution of homorepeats can be influenced by DNA replication slippage, when DNA polymerases either add or skip repeats of nucleotides. While there are some diseases known to be caused by abnormal changes in the length of amino acid homorepeats, naturally occurring variations in homorepeat length remain relatively unexplored. In our study, we examined the variation in amino acid homorepeat length of human individuals by analyzing 125 748 exomes, as well as 15 708 whole genomes. Our analyses revealed significant variability in homorepeat length across the human population, indicating that these motifs are prone to mutations at higher rates than non repeat sequences. We focused our study on glutamine homorepeats, also known as polyQ sequences, and found that shorter polyQ sequences tend to exhibit greater length variation, while longer ones primarily undergo deletions. Notably, polyQ sequencesthat are more conserved across primates tend to show less variation within the human population, indicating stronger selective pressure to maintain their length. Overall, our results demonstrate that there is large natural variation in the length of homorepeats within the human population, with no apparent impact on observable traits.

Correction: Whole genome variation in 27 Mexican indigenous populations, demographic and biomedical insights.

[This corrects the article DOI: 10.1371/journal.pone.0249773.].

Transcriptome Analysis Reveals Cr(VI) Adaptation Mechanisms in Klebsiella sp. Strain AqSCr.

Klebsiella sp. strain AqSCr, isolated from Cr(VI)-polluted groundwater, reduces Cr(VI) both aerobically and anaerobically and resists up 34 mM Cr(VI); this resistance is independent of the ChrA efflux transporter. In this study, we report the whole genome sequence and the transcriptional profile by RNA-Seq of strain AqSCr under Cr(VI)-adapted conditions and found 255 upregulated and 240 downregulated genes compared to controls without Cr(VI) supplementation. Genes differentially transcribed were mostly associated with oxidative stress response, DNA repair and replication, sulfur starvation response, envelope-osmotic stress response, fatty acid (FA) metabolism, ribosomal subunits, and energy metabolism. Among them, genes not previously associated with chromium resistance, for example, cybB, encoding a putative superoxide oxidase (SOO), gltA2, encoding an alternative citrate synthase, and des, encoding a FA desaturase, were upregulated. The sodA gene encoding a manganese superoxide dismutase was upregulated in the presence of Cr(VI), whereas sodB encoding an iron superoxide dismutase was downregulated. Cr(VI) resistance mechanisms in strain AqSCr seem to be orchestrated by the alternative sigma factors fecl, rpoE, and rpoS (all of them upregulated). Membrane lipid analysis of the Cr(IV)-adapted strain showed a lower proportion of unsaturated lipids with respect to the control, which we hypothesized could result from unsaturated lipid peroxidation followed by degradation, together with de novo synthesis mediated by the upregulated FA desaturase-encoding gene, des. This report helps to elucidate both Cr(VI) toxicity targets and global bacterial response to Cr(VI).

Whole genome variation in 27 Mexican indigenous populations, demographic and biomedical insights.

There has been limited study of Native American whole genome diversity to date, which impairs effective implementation of personalized medicine and a detailed description of its demographic history. Here we report high coverage whole genome sequencing of 76 unrelated individuals, from 27 indigenous groups across Mexico, with more than 97% average Native American ancestry. On average, each individual has 3.26 million Single Nucleotide Variants and short indels, that together comprise a catalog of 9,737,152 variants, 44,118 of which are novel. We report 497 common Single Nucleotide Variants (with allele frequency > 5%) mapped to drug responses and 316,577 in enhancer or promoter elements; interestingly we found some of these enhancer variants in PPARG, a nuclear receptor involved in highly prevalent health problems in Mexican population, such as obesity, diabetes, and insulin resistance. By detecting signals of positive selection we report 24 enriched key pathways under selection, most of them related to immune mechanisms. No missense variants in ACE2, the receptor responsible for the entry of the SARS CoV-2 virus, were found in any individual. Population genomics and phylogenetic analyses demonstrated stratification in a Northern-Central-Southern axis, with major substructure in the Central region. The Seri, a northern group with the most genetic divergence in our study, showed a distinctive genomic context with the most novel variants, and the most population specific genotypes. Genome-wide analysis showed that the average haplotype blocks are longer in Native Mexicans than in other world populations. With this dataset we describe previously undetected population level variation in Native Mexicans, helping to reduce the gap in genomic data representation of such groups.

High-throughput transcriptome sequencing and comparative analysis of Escherichia coli and Schizosaccharomyces pombe in respiratory and fermentative growth.

In spite of increased complexity in eukaryotes compared to prokaryotes, several basic metabolic and regulatory processes are conserved. Here we explored analogies in the eubacteria Escherichia coli and the unicellular fission yeast Schizosaccharomyces pombe transcriptomes under two carbon sources: 2% glucose; or a mix of 2% glycerol and 0.2% sodium acetate using the same growth media and growth phase. Overall, twelve RNA-seq libraries were constructed. A total of 593 and 860 genes were detected as differentially expressed for E. coli and S. pombe, respectively, with a log2 of the Fold Change ≥ 1 and False Discovery Rate ≤ 0.05. In aerobic glycolysis, most of the expressed genes were associated with cell proliferation in both organisms, including amino acid metabolism and glycolysis. In contrast in glycerol/acetate condition, genes related to flagellar assembly and membrane proteins were differentially expressed such as the general transcription factors fliA, flhD, flhC, and flagellum assembly genes were detected in E. coli, whereas in S. pombe genes for hexose transporters, integral membrane proteins, galactose metabolism, and ncRNAs related to cellular stress were overexpressed. In general, our study shows that a conserved "foraging behavior" response is observed in these eukaryotic and eubacterial organisms in gluconeogenic carbon sources.

Specialization of the Reiterated Copies of the Heterodimeric Integration Host Factor Genes in Geobacter sulfurreducens.

Integration host factor (IHF) is a widely distributed small heterodimeric protein member of the bacterial Nucleoid-Associated Proteins (NAPs), implicated in multiple DNA regulatory processes. IHF recognizes a specific DNA sequence and induces a large bend of the nucleic acid. IHF function has been mainly linked with the regulation of RpoN-dependent promoters, where IHF commonly recognizes a DNA sequence between the enhancer-binding region and the promoter, facilitating a close contact between the upstream bound activator and the promoter bound, RNA polymerase. In most proteobacteria, the genes encoding IHF subunits (ihfA and ihfB) are found in a single copy. However, in some Deltaproteobacteria, like Geobacter sulfurreducens, those genes are duplicated. To date, the functionality of IHF reiterated encoding genes is unknown. In this work, we achieved the functional characterization of the ihfA-1, ihfA-2, ihfB-1, and ihfB-2 from G. sulfurreducens. Unlike the ΔihfA-2 or ΔihfB-1 strains, single gene deletion in ihfA-1 or ihfB-2, provokes an impairment in fumarate and Fe(III) citrate reduction. Accordingly, sqRT-PCR experiments showed that ihfA-1 and ihfB-2 were expressed at higher levels than ihfA-2 and ihfB-1. In addition, RNA-Seq analysis of the ΔihfA-1 and ΔihfB-2 strains revealed a total of 89 and 122 differentially expressed genes, respectively. Furthermore, transcriptional changes in 25 genes were shared in both mutant strains. Among these genes, we confirmed the upregulation of the pilA-repressor, GSU1771, and downregulation of the triheme-cytochrome (pgcA) and the aconitate hydratase (acnA) genes by RT-qPCR. EMSA experiments also demonstrated the direct binding of IHF to the upstream promoter regions of GSU1771, pgcA and acnA. PilA changes in ΔihfA-1 and ΔihfB-2 strains were also verified by immunoblotting. Additionally, heme-staining of subcellular fractions in ΔihfA-1 and ΔihfB-2 strains revealed a remarkable deficit of c-type cytochromes. Overall, our data indicate that at least during fumarate and Fe(III) citrate reduction, the functional IHF regulator is likely assembled by the products of ihfA-1 and ihfB-2. Also, a role of IHF controlling expression of multiple genes (other than RpoN-dependent) affects G. sulfurreducens physiology and extracellular electron transfer.

Marine Archaeon Methanosarcina acetivorans Enhances Polyphosphate Metabolism Under Persistent Cadmium Stress.

Phosphate metabolism was studied to determine whether polyphosphate (polyP) pools play a role in the enhanced resistance against Cd2+ and metal-removal capacity of Cd2+-preadapted (CdPA) Methanosarcina acetivorans. Polyphosphate kinase (PPK), exopolyphosphatase (PPX) and phosphate transporter transcript levels and their activities increased in CdPA cells compared to control (Cnt) cells. K+ inhibited recombinant Ma-PPK and activated Ma-PPX, whereas divalent cations activated both enzymes. Metal-binding polyP and thiol-containing molecule contents, Cd2+-removal, and biofilm synthesis were significantly higher in CdPA cells >Cnt cells plus a single addition of Cd2+>Cnt cells. Also, CdPA cells showed a higher number of cadmium, sulfur, and phosphorus enriched-acidocalcisomes than control cells. Biochemical and physiological phenotype exhibited by CdPA cells returned to that of Cnt cells when cultured without Cd2+. Furthermore, no differences in the sequenced genomes upstream and downstream of the genes involved in Cd2+ resistance were found between CdPA and Cnt cells, suggesting phenotype loss rather than genome mutations induced by chronic Cd2+-exposure. Instead, a metabolic adaptation induced by Cd2+ stress was apparent. The dynamic ability of M. acetivorans to change its metabolism, depending on the environmental conditions, may be advantageous to remove cadmium in nature and biodigesters.

Demographic history and biologically relevant genetic variation of Native Mexicans inferred from whole-genome sequencing.

Understanding the genetic structure of Native American populations is important to clarify their diversity, demographic history, and to identify genetic factors relevant for biomedical traits. Here, we show a demographic history reconstruction from 12 Native American whole genomes belonging to six distinct ethnic groups representing the three main described genetic clusters of Mexico (Northern, Southern, and Maya). Effective population size estimates of all Native American groups remained below 2,000 individuals for up to 10,000 years ago. The proportion of missense variants predicted as damaging is higher for undescribed (~ 30%) than for previously reported variants (~ 15%). Several variants previously associated with biological traits are highly frequent in the Native American genomes. These findings suggest that the demographic and adaptive processes that occurred in these groups shaped their genetic architecture and could have implications in biological processes of the Native Americans and Mestizos of today.

The unphosphorylated form of the PilR two-component system regulates pilA gene expression in Geobacter sulfurreducens.

In Geobacter sulfurreducens, metal reduction and generation of bioelectricity require the participation of several elements, and among them, the type IV pili has an essential role. The pilus is composed of multiple PilA monomers. Expression of pilA gene depends mainly on the σ54 factor and the response regulator protein PilR. In this work, we characterized the role of the PilS-PilR two-component system in the regulation of the pilA gene expression. Experimental evidence indicates that PilS is autophosphorylated at the His-334 residue, which in turn is transferred to the conserved Asp-53 in PilR. Contrary to other PilS-PilR systems, substitution D53N in PilR resulted in higher activation of the pilA gene. By using a pilA::luxCDABE fusion with different promoter fragments and in vitro DNA-binding assays, we demonstrated the existence of multiple functional PilR binding sites. A regulatory model in which the non-phosphorylated PilR protein directs activation of pilA expression by binding to two sites in the promoter region of this gene is presented.

Acetate biostimulation as an effective treatment for cleaning up alkaline soil highly contaminated with Cr(VI).

Stimulation of microbial reduction of Cr(VI) to the less toxic and less soluble Cr(III) through electron donor addition has been regarded as a promising approach for the remediation of chromium-contaminated soil and groundwater sites. However, each site presents different challenges; local physicochemical characteristics and indigenous microbial communities influence the effectiveness of the biostimulation processes. Here, we show microcosm assays stimulation of microbial reduction of Cr(VI) in highly alkaline and saline soil samples from a long-term contaminated site in Guanajuato, Mexico. Acetate was effective promoting anaerobic microbial reduction of 15 mM of Cr(VI) in 25 days accompanied by an increase in pH from 9 to 10. Our analyses showed the presence of Halomonas, Herbaspirillum, Nesterenkonia/Arthrobacter, and Bacillus species in the soil sample collected. Moreover, from biostimulated soil samples, it was possible to isolate Halomonas spp. strains able to grow at 32 mM of Cr(VI). Additionally, we found that polluted groundwater has bacterial species different to those found in soil samples with the ability to resist and reduce chromate using acetate and yeast extract as electron donors.

Global transcriptional start site mapping in Geobacter sulfurreducens during growth with two different electron acceptors.

Geobacter sulfurreducens is an anaerobic soil bacterium that is involved in biogeochemical cycles of elements such as Fe and Mn. Although significant progress has been made in the understanding of the electron transfer processes in G. sulfurreducens, little is known about the regulatory mechanisms involved in their control. To expand the study of gene regulation in G. sulfurreducens, we carried out a genome-wide identification of transcription start sites (TSS) by 5'RACE and by deep RNA sequencing of primary mRNAs in two growth conditions. TSSs were identified along G. sulfurreducens genome and over 50% of them were located in the upstream region of the associated gene, and in some cases we detected genes with more than one TSS. Our global mapping of TSSs contributes with valuable information, which is needed for the study of transcript structure and transcription regulation signals and can ultimately contribute to the understanding of transcription initiation phenomena in G. sulfurreducens.

Structural properties of prokaryotic promoter regions correlate with functional features.

The structural properties of the DNA molecule are known to play a critical role in transcription. In this paper, the structural profiles of promoter regions were studied within the context of their diversity and their function for eleven prokaryotic species; Escherichia coli, Klebsiella pneumoniae, Salmonella Typhimurium, Pseudomonas auroginosa, Geobacter sulfurreducens Helicobacter pylori, Chlamydophila pneumoniae, Synechocystis sp., Synechoccocus elongates, Bacillus anthracis, and the archaea Sulfolobus solfataricus. The main anchor point for these promoter regions were transcription start sites identified through high-throughput experiments or collected within large curated databases. Prokaryotic promoter regions were found to be less stable and less flexible than the genomic mean across all studied species. However, direct comparison between species revealed differences in their structural profiles that can not solely be explained by the difference in genomic GC content. In addition, comparison with functional data revealed that there are patterns in the promoter structural profiles that can be linked to specific functional loci, such as sigma factor regulation or transcription factor binding. Interestingly, a novel structural element clearly visible near the transcription start site was found in genes associated with essential cellular functions and growth in several species. Our analyses reveals the great diversity in promoter structural profiles both between and within prokaryotic species. We observed relationships between structural diversity and functional features that are interesting prospects for further research to yet uncharacterized functional loci defined by DNA structural properties.

Exopolyphosphatase of Pseudomonas aeruginosa is essential for the production of virulence factors, and its expression is controlled by NtrC and PhoB acting at two interspaced promoters.

The exopolyphosphatase (Ppx) of Pseudomonas aeruginosa is encoded by the PA5241 gene (ppx). Ppx catalyses the hydrolysis of inorganic polyphosphates to orthophosphate (Pi). In the present work, we identified and characterized the promoter region of ppx and its regulation under environmental stress conditions. The role of Ppx in the production of several virulence factors was demonstrated through studies performed on a ppx null mutant. We found that ppx is under the control of two interspaced promoters, dually regulated by nitrogen and phosphate limitation. Under nitrogen-limiting conditions, its expression was controlled from a σ(54)-dependent promoter activated by the response regulator NtrC. However, under Pi limitation, the expression was controlled from a σ(70) promoter, activated by PhoB. Results obtained from the ppx null mutant demonstrated that Ppx is involved in the production of virulence factors associated with both acute infection (e.g. motility-promoting factors, blue/green pigment production, C6-C12 quorum-sensing homoserine lactones) and chronic infection (e.g. rhamnolipids, biofilm formation). Molecular and physiological approaches used in this study indicated that P. aeruginosa maintains consistently proper levels of Ppx regardless of environmental conditions. The precise control of ppx expression appeared to be essential for the survival of P. aeruginosa and the occurrence of either acute or chronic infection in the host.

Transcription factors in Escherichia coli prefer the holo conformation.

The transcriptional regulatory network of Escherichia coli K-12 is among the best studied gene networks of any living cell. Transcription factors bind to DNA either with their effector bound (holo conformation), or as a free protein (apo conformation) regulating transcription initiation. By using RegulonDB, the functional conformations (holo or apo) of transcription factors, and their mode of regulation (activator, repressor, or dual) were exhaustively analyzed. We report a striking discovery in the architecture of the regulatory network, finding a strong under-representation of the apo conformation (without allosteric metabolite) of transcription factors when binding to their DNA sites to activate transcription. This observation is supported at the level of individual regulatory interactions on promoters, even if we exclude the promoters regulated by global transcription factors, where three-quarters of the known promoters are regulated by a transcription factor in holo conformation. This genome-scale analysis enables us to ask what are the implications of these observations for the physiology and for our understanding of the ecology of E. coli. We discuss these ideas within the framework of the demand theory of gene regulation.

Evidence classification of high-throughput protocols and confidence integration in RegulonDB.

RegulonDB provides curated information on the transcriptional regulatory network of Escherichia coli and contains both experimental data and computationally predicted objects. To account for the heterogeneity of these data, we introduced in version 6.0, a two-tier rating system for the strength of evidence, classifying evidence as either 'weak' or 'strong' (Gama-Castro,S., Jimenez-Jacinto,V., Peralta-Gil,M. et al. RegulonDB (Version 6.0): gene regulation model of Escherichia Coli K-12 beyond transcription, active (experimental) annotated promoters and textpresso navigation. Nucleic Acids Res., 2008;36:D120-D124.). We now add to our classification scheme the classification of high-throughput evidence, including chromatin immunoprecipitation (ChIP) and RNA-seq technologies. To integrate these data into RegulonDB, we present two strategies for the evaluation of confidence, statistical validation and independent cross-validation. Statistical validation involves verification of ChIP data for transcription factor-binding sites, using tools for motif discovery and quality assessment of the discovered matrices. Independent cross-validation combines independent evidence with the intention to mutually exclude false positives. Both statistical validation and cross-validation allow to upgrade subsets of data that are supported by weak evidence to a higher confidence level. Likewise, cross-validation of strong confidence data extends our two-tier rating system to a three-tier system by introducing a third confidence score 'confirmed'. Database URL: http://regulondb.ccg.unam.mx/

RegulonDB v8.0: omics data sets, evolutionary conservation, regulatory phrases, cross-validated gold standards and more.

This article summarizes our progress with RegulonDB (http://regulondb.ccg.unam.mx/) during the past 2 years. We have kept up-to-date the knowledge from the published literature regarding transcriptional regulation in Escherichia coli K-12. We have maintained and expanded our curation efforts to improve the breadth and quality of the encoded experimental knowledge, and we have implemented criteria for the quality of our computational predictions. Regulatory phrases now provide high-level descriptions of regulatory regions. We expanded the assignment of quality to various sources of evidence, particularly for knowledge generated through high-throughput (HT) technology. Based on our analysis of most relevant methods, we defined rules for determining the quality of evidence when multiple independent sources support an entry. With this latest release of RegulonDB, we present a new highly reliable larger collection of transcription start sites, a result of our experimental HT genome-wide efforts. These improvements, together with several novel enhancements (the tracks display, uploading format and curational guidelines), address the challenges of incorporating HT-generated knowledge into RegulonDB. Information on the evolutionary conservation of regulatory elements is also available now. Altogether, RegulonDB version 8.0 is a much better home for integrating knowledge on gene regulation from the sources of information currently available.

Genetic changes during a laboratory adaptive evolution process that allowed fast growth in glucose to an Escherichia coli strain lacking the major glucose transport system.

BACKGROUND: Escherichia coli strains lacking the phosphoenolpyruvate: carbohydrate phosphotransferase system (PTS), which is the major bacterial component involved in glucose transport and its phosphorylation, accumulate high amounts of phosphoenolpyruvate that can be diverted to the synthesis of commercially relevant products. However, these strains grow slowly in glucose as sole carbon source due to its inefficient transport and metabolism. Strain PB12, with 400% increased growth rate, was isolated after a 120 hours adaptive laboratory evolution process for the selection of faster growing derivatives in glucose. Analysis of the genetic changes that occurred in the PB12 strain that lacks PTS will allow a better understanding of the basis of its growth adaptation and, therefore, in the design of improved metabolic engineering strategies for enhancing carbon diversion into the aromatic pathways. RESULTS: Whole genome analyses using two different sequencing methodologies: the Roche NimbleGen Inc. comparative genome sequencing technique, and high throughput sequencing with Illumina Inc. GAIIx, allowed the identification of the genetic changes that occurred in the PB12 strain. Both methods detected 23 non-synonymous and 22 synonymous point mutations. Several non-synonymous mutations mapped in regulatory genes (arcB, barA, rpoD, rna) and in other putative regulatory loci (yjjU, rssA and ypdA). In addition, a chromosomal deletion of 10,328 bp was detected that removed 12 genes, among them, the rppH, mutH and galR genes. Characterization of some of these mutated and deleted genes with their functions and possible functions, are presented. CONCLUSIONS: The deletion of the contiguous rppH, mutH and galR genes that occurred simultaneously, is apparently the main reason for the faster growth of the evolved PB12 strain. In support of this interpretation is the fact that inactivation of the rppH gene in the parental PB11 strain substantially increased its growth rate, very likely by increasing glycolytic mRNA genes stability. Furthermore, galR inactivation allowed glucose transport by GalP into the cell. The deletion of mutH in an already stressed strain that lacks PTS is apparently responsible for the very high mutation rate observed.

Evolution of a new function in an esterase: simple amino acid substitutions enable the activity present in the larger paralog, BioH.

Gene duplication and divergence are essential processes for the evolution of new activities. Divergence may be gradual, involving simple amino acid residue substitutions, or drastic, such that larger structural elements are inserted, deleted or rearranged. Vast protein sequence comparisons, supported by some experimental evidence, argue that large structural modifications have been necessary for certain catalytic activities to evolve. However, it is not clear whether these activities could not have been attained by gradual changes. Interestingly, catalytic promiscuity could play a fundamental evolutionary role: a preexistent secondary activity could be increased by simple amino acid residue substitutions that do not affect the enzyme's primary activity. The promiscuous profile of the enzyme may be modified gradually by genetic drift, making a pool of potentially useful activities that can be selected before duplication. In this work, we used random mutagenesis and in vivo selection to evolve the Pseudomonas aeruginosa PAO1 carboxylesterase PA3859, a small protein, to attain the function of BioH, a much larger paralog involved in biotin biosynthesis. BioH was chosen as a target activity because it provides a highly sensitive selection for evolved enzymatic activities by auxotrophy complementation. After only two cycles of directed evolution, mutants with the ability to efficiently complement biotin auxotrophy were selected. The in vivo and in vitro characterization showed that the activity of one of our mutant proteins was similar to that of the wild-type BioH enzyme. Our results demonstrate that it is possible to evolve enzymatic activities present in larger proteins by discrete amino acid substitutions.

Evolutionary walk between (ß/α)(8) barrels: catalytic migration from triosephosphate isomerase to thiamin phosphate synthase.

The functionally versatile (ß/α)(8) barrel scaffold was used to migrate triosephosphate isomerase (TPI) to thiamin phosphate synthase (TPS) activity, two enzymes that share the same fold but catalyze unrelated reactions through different mechanisms. The high sensitivity of the selection methodology was determinant to succeed in finding proteins with the desired activity. A combination of rational design and random mutagenesis was used to achieve the desired catalytic migration. One of the parallel directed evolution strategies followed resulted in TPI derivatives able to complement the thiamin phosphate auxotrophy phenotype of an Escherichia coli strain deleted of thiE, the gene that codes for TPS. Successive rounds of directed evolution resulted in better complementing TPI variants. Biochemical characterization of some of the evolved TPI clones demonstrated that the K(m) for the TPS substrates was similar to that of the native TPS; however and in agreement with the very slow complementation phenotype, the k(cat) was 4 orders of magnitude lower, indicating that substrate binding played a major role on selection. Interestingly, the crystal structure of the most proficient variant showed a slightly modified TPI active site occupied by a thiamin-phosphate-like molecule. Substitution of key residues in this region reduced TPS activity, strongly suggesting that this is also the catalytic site for the evolved TPS activity. The presence of the TPS reaction product at the active site explains the fast inactivation of the enzyme observed. In conclusion, by combining rational design, random mutagenesis and a very sensitive selection, it is possible to achieve enzymatic activity migration.

Choline catabolism, σ54 factor and NtrC are required for the full expression of the Pseudomonas aeruginosa phosphorylcholine phosphatase gene.

Choline favors the pathogenesis of Pseudomonas aeruginosa because hemolytic phospholipase C and phosphorylcholine phosphatase (PchP) are synthesized as a consequence of its catabolism. The experiments performed here resulted in the identification of the factors that regulate both the catabolism of choline and the gene coding for PchP. We have also identified and characterized the promoter of the pchP gene, its transcriptional organization and the factors that affect its expression. Deletion analyses reveal that the region between -188 and -68 contains all controlling elements necessary for pchP expression: a hypothetical -12/-24 promoter element, a consensus sequence for the integration host factor (-141/-133), and a palindromic sequence resembling a binding site for a potential enhancer binding protein (-190/-174). Our data also demonstrate that choline catabolism and NtrC (nitrogen regulatory protein) are necessary for the full expression of pchP and is partially dependent on σ(54) factor.