Regulation of biosynthesis and transport of aromatic amino acids in low-GC Gram-positive bacteria

Understanding the transcriptomic response of Lactiplantibacillus pentosus LPG1 during Spanish-style green table olive fermentations

Lactiplantibacillus pentosus (Lbp. pentosus) is a species of lactic acid bacteria with a great relevance during the table olive fermentation process, with ability to form non-pathogenic biofilms on olive epidermis. The objective of this work is to deepen into the genetic mechanisms of adaptation of Lpb. pentosus LPG1 during Spanish-style green table olive fermentations, as well as to obtain a better understanding of the mechanisms of adherence of this species to the fruit surface. For this purpose, we have carried out a transcriptomic analysis of the differential gene expression of this bacterium during 60 days of fermentation in both brine and biofilms ecosystems. In brines, it was noticed that a total of 235 genes from Lpb. pentosus LPG1 were differentially expressed during course of fermentation and grouped into 9 clusters according to time-course analysis. Transport and metabolism of carbohydrates and amino acids, energy production, lactic acid and exopolysaccharide synthesis genes increased their expression in the planktonic cells during course of fermentation. On the other hand, expression of genes associated to stress response, bacteriocin synthesis and membrane protein decreased. A total of 127 genes showed significant differential expression between Lpb. pentosus LPG1 planktonic (brine) and sessile (biofilms) cells at the end of fermentation process (60 days). Among the 64 upregulated genes in biofilms, we found genes involved in adhesion (strA), exopolysaccharide production (ywqD, ywqE, and wbnH), cell shape and elongation (MreB), and well as prophage excision. Deeping into the genetic bases of beneficial biofilm formation by Lpb. pentosus strains with probiotic potential will help to turn this fermented vegetable into a carrier of beneficial microorganisms to the final consumers.

Meta-transcriptomic comparison of two sponge holobionts feeding on coral- and macroalgal-dissolved organic matter

Background

Sponge holobionts (i.e., the host and its associated microbiota) play a key role in the cycling of dissolved organic matter (DOM) in marine ecosystems. On coral reefs, an ecological shift from coral-dominated to algal-dominated ecosystems is currently occurring. Given that benthic corals and macroalgae release different types of DOM, in different abundances and with different bioavailability to sponge holobionts, it is important to understand how the metabolic activity of the host and associated microbiota change in response to the exposure to both DOM sources. Here, we look at the differential gene expression of two sponge holobionts 6 hours after feeding on naturally sourced coral- and macroalgal-DOM using RNA sequencing and meta-transcriptomic analysis.

Results

We found a slight, but significant differential gene expression in the comparison between the coral- and macroalgal-DOM treatments in both the high microbial abundance sponge Plakortis angulospiculatus and the low microbial abundance sponge Haliclona vansoesti. In the hosts, processes that regulate immune response, signal transduction, and metabolic pathways related to cell proliferation were elicited. In the associated microbiota carbohydrate metabolism was upregulated in both treatments, but coral-DOM induced further lipid and amino acids biosynthesis, while macroalgal-DOM caused a stress response. These differences could be driven by the presence of distinct organic macronutrients in the two DOM sources and of small pathogens or bacterial virulence factors in the macroalgal-DOM.

Conclusions

This work provides two new sponge meta-transcriptomes and a database of putative genes and genetic pathways that are involved in the differential processing of coral- versus macroalgal-DOM as food source to sponges with high and low abundances of associated microbes. These pathways include carbohydrate metabolism, signaling pathways, and immune responses. However, the differences in the meta-transcriptomic responses of the sponge holobionts after 6 hours of feeding on the two DOM sources were small. Longer-term responses to both DOM sources should be assessed to evaluate how the metabolism and the ecological function of sponges will be affected when reefs shift from coral towards algal dominance.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08893-y.

Streptococcus suis TrpX is part of a tryptophan uptake system, and its expression is regulated by a T-box regulatory element

Streptococcus suis, a common member of the porcine respiratory microbiota, can cause life-threatening diseases in pigs as well as humans. A previous study identified the gene trpX as conditionally essential for in vivo survival by intrathecal infection of pigs with a transposon library of S. suis strain 10. Here, we characterized trpX, encoding a putative tryptophan/tyrosine transport system substrate-binding protein, in more detail. We compared growth capacities of the isogenic trpX-deficient mutant derivative strain 10∆trpX with its parent. Growth experiments in chemically defined media (CDM) revealed that growth of 10∆trpX depended on tryptophan concentration, suggesting TrpX involvement in tryptophan uptake. We demonstrated that trpX is part of an operon structure and co-transcribed with two additional genes encoding a putative permease and ATPase, respectively. Bioinformatics analysis identified a putative tryptophan T-box riboswitch in the 5′ untranslated region of this operon. Finally, qRT-PCR and a reporter activation assay revealed trpX mRNA induction under tryptophan-limited conditions. In conclusion, our study showed that TrpX is part of a putative tryptophan ABC transporter system regulated by a T-box riboswitch probably functioning as a substrate-binding protein. Due to the tryptophan auxotrophy of S. suis, TrpX plays a crucial role for metabolic adaptation and growth during infection.

The scfCDE Operon Encodes a Predicted ABC Importer Required for Fitness and Virulence during Group A Streptococcus Invasive Infection

As a strict human pathogen, Streptococcus pyogenes (group A Streptococcus, or GAS) causes a wide range of infections, from superficial to life-threatening diseases, upon dissemination. Thus, it is necessary to gain a better understanding of how GAS successfully overcomes host-mediated challenges and infects various host niches. We previously identified subcutaneous fitness (scf) genes in the clinically relevant wild-type (WT) GAS M1T1 5448 strain that are critical for fitness during murine soft-tissue infection at both 24 h and 48 h postinfection. The uncharacterized locus scfCDE was transcribed as an operon and is predicted to encode an ABC importer for nutrient uptake (e.g., amino acids). Individual scfCDE deletion mutants grew comparably to WT 5448 in rich medium but exhibited reduced fitness during competitive growth in murine soft tissue and in nutrient-limiting chemically defined medium (CDM). A deletion of the permease gene scfD resulted in a monoculture growth defect in CDM that could be rescued by addition of excess peptides, suggesting a role as an amino acid importer. Interestingly, the ΔscfC substrate-binding and ΔscfD permease mutants, but not the ΔscfE ATPase mutant, were highly attenuated in murine soft tissue. Moreover, all three genes were required for GAS survival in human blood, indicating their impact is not limited to superficial infections. As such, scfCDE plays an integral role in enhancing GAS adaptation during localized infection as well as dissemination to deeper host environments. Since scfCDE is conserved throughout Firmicutes, this work may contribute to the development of therapeutic strategies against GAS and other Gram-positive pathogens.

Genome, transcriptome and fermentation analyses of Lactobacillus plantarum LY-78 provide new insights into the mechanism of phenyllactate biosynthesis in lactic acid bacteria

Phenyllactate (PLA) is found in a variety of fermented foods and is a promising antibacterial agent, drug and plastic synthetic precursor. Previous studies have shown that PLA is a product of Phe catabolism in lactic acid bacteria (LAB), and PLA biosynthesis is mainly related to lactate dehydrogenases (LDHs). Here, the genome, transcriptome and fermentation characteristics of PLA-producing Lactobacillus plantarum LY-78 were studied. The fermentation experiments demonstrated that L. plantarum LY-78 possesses the ability to synthesize PLA de novo. Secondly, the genome and transcriptome analyses revealed candidate pathways, operons and key genes for PLA biosynthesis in the strain. Finally, genome-wide transcriptome analysis revealed significant changes in the expression profile of strain LY-78 in the absence and presence of PPA. Overall, this work demonstrates for the first time that PLA can be a by-product of Phe anabolism in LAB, provides new insights and evidence for elucidating the mechanism of PLA biosynthesis in LAB, and may provide new candidate genes and research strategies for future PLA biosynthesis applications.

Enterococcal Genetics

The study of the genetics of enterococci has focused heavily on mobile genetic elements present in these organisms, the complex regulatory circuits used to control their mobility, and the antibiotic resistance genes they frequently carry. Recently, more focus has been placed on the regulation of genes involved in the virulence of the opportunistic pathogenic species Enterococcus faecalis and Enterococcus faecium. Little information is available concerning fundamental aspects of DNA replication, partition, and division; this article begins with a brief overview of what little is known about these issues, primarily by comparison with better-studied model organisms. A variety of transcriptional and posttranscriptional mechanisms of regulation of gene expression are then discussed, including a section on the genetics and regulation of vancomycin resistance in enterococci. The article then provides extensive coverage of the pheromone-responsive conjugation plasmids, including sections on regulation of the pheromone response, the conjugative apparatus, and replication and stable inheritance. The article then focuses on conjugative transposons, now referred to as integrated, conjugative elements, or ICEs, and concludes with several smaller sections covering emerging areas of interest concerning the enterococcal mobilome, including nonpheromone plasmids of particular interest, toxin-antitoxin systems, pathogenicity islands, bacteriophages, and genome defense.

Physicochemical Characterization of Biofluid Metabolites from Liquid Residual of Tuna Fish (Euthynnus affinis) throughout Refrigerated Storage Condition

The cold storage condition and use of chemical treatment to preserve the fish appearance sometimes cause difficulties to the consumers to estimate the freshness of fish in the market. However, during fish deterioration, some compound is released or formed due to microbial and biochemical process. Identification of released compound during fish spoilage is a crucial step in understanding the degree of spoilage. This study characterizes the physicochemical changes of metabolites biofluids from liquid residual of tuna fish (Euthynnus affinis) during refrigerated storage condition. Tuna fish were kept in ice at 0°C and stored in cold room (~4°C) for seven days in order to study the changes in fish freshness and loss of quality through the storage period. Liquid residual of fish was collected at 0, 1, 2, 3, 4, 5, 6, and 7 days of storage. LC-MS/MS analysis was carried out to determine the possible dominant compound which was later identified as creatine and phenylalanine. Quantification of creatine and phenylalanine using HPLC with UV detector found that creatine and phenylalanine increased significantly up to day 4 and day 5 upon storage time for creatine and phenylalanine, respectively (p<0.05). The liquid residual pH increased from 6.5 on day 0 to 7.5 on day 7 (p<0.05). Changes in chemical compounds were supported with physical analysis on gills colour of spoilage fish. L⁎ and a⁎ values decreased with storage time from 41.08 to 24.76 and 18.34 to 10.40, respectively, while b⁎ value increased from -3.80 to -0.46 (p<0.05). The finding of biofluid derived compounds was found as useful and alternative indicators of fish freshness in later study on the development of optical biosensor.

Deep sequencing of the oral microbiome reveals signatures of periodontal disease

The oral microbiome, the complex ecosystem of microbes inhabiting the human mouth, harbors several thousands of bacterial types. The proliferation of pathogenic bacteria within the mouth gives rise to periodontitis, an inflammatory disease known to also constitute a risk factor for cardiovascular disease. While much is known about individual species associated with pathogenesis, the system-level mechanisms underlying the transition from health to disease are still poorly understood. Through the sequencing of the 16S rRNA gene and of whole community DNA we provide a glimpse at the global genetic, metabolic, and ecological changes associated with periodontitis in 15 subgingival plaque samples, four from each of two periodontitis patients, and the remaining samples from three healthy individuals. We also demonstrate the power of whole-metagenome sequencing approaches in characterizing the genomes of key players in the oral microbiome, including an unculturable TM7 organism. We reveal the disease microbiome to be enriched in virulence factors, and adapted to a parasitic lifestyle that takes advantage of the disrupted host homeostasis. Furthermore, diseased samples share a common structure that was not found in completely healthy samples, suggesting that the disease state may occupy a narrow region within the space of possible configurations of the oral microbiome. Our pilot study demonstrates the power of high-throughput sequencing as a tool for understanding the role of the oral microbiome in periodontal disease. Despite a modest level of sequencing (~2 lanes Illumina 76 bp PE) and high human DNA contamination (up to ~90%) we were able to partially reconstruct several oral microbes and to preliminarily characterize some systems-level differences between the healthy and diseased oral microbiomes.

Fast, sensitive discovery of conserved genome-wide motifs

Regulatory sites that control gene expression are essential to the proper functioning of cells, and identifying them is critical for modeling regulatory networks. We have developed Magma (Multiple Aligner of Genomic Multiple Alignments), a software tool for multiple species, multiple gene motif discovery. Magma identifies putative regulatory sites that are conserved across multiple species and occur near multiple genes throughout a reference genome. Magma takes as input multiple alignments that can include gaps. It uses efficient clustering methods that make it about 70 times faster than PhyloNet, a previous program for this task, with slightly greater sensitivity. We ran Magma on all non-coding DNA conserved between Caenorhabditis elegans and five additional species, about 70 Mbp in total, in <4 h. We obtained 2,309 motifs with lengths of 6-20 bp, each occurring at least 10 times throughout the genome, which collectively covered about 566 kbp of the genomes, approximately 0.8% of the input. Predicted sites occurred in all types of non-coding sequence but were especially enriched in the promoter regions. Comparisons to several experimental datasets show that Magma motifs correspond to a variety of known regulatory motifs.

Biosynthesis of teichoic acids in Streptococcus pneumoniae and closely related species: lessons from genomes

The cell wall of Streptococcus pneumoniae contains an unusually complex wall teichoic acid (WTA), which has identical repeating units as the membrane-anchored lipoteichoic acid (LTA). Both polymers share a common cytoplasmic pathway of precursor synthesis, but several TA enzymes have remained elusive. Bioinformatic analysis of the genome of various pneumococcal strains, including choline-independent mutant strains, has allowed us to identify the missing TA genes. We present here the deduced complete pathways of WTA and LTA synthesis in S. pneumoniae and point to the variations occurring in different pneumococcal strains and in closely related species such as Streptococcus oralis and Streptococcus mitis.

Characterization of central carbon metabolism of Streptococcus pneumoniae by isotopologue profiling

The metabolism of Streptococcus pneumoniae was studied by isotopologue profiling after bacterial cultivation in chemically defined medium supplemented with [U-(13)C(6)]- or [1,2-(13)C(2)]glucose. GC/MS analysis of protein-derived amino acids showed lack of (13)C label in amino acids that were also essential for pneumococcal growth. Ala, Ser, Asp, and Thr displayed high (13)C enrichments, whereas Phe, Tyr, and Gly were only slightly labeled. The analysis of the labeling patterns showed formation of triose phosphate and pyruvate via the Embden-Meyerhof-Parnas pathway. The labeling patterns of Asp and Thr suggested formation of oxaloacetate exclusively via the phosphoenolpyruvate carboxylase reaction. Apparently, α-ketoglutarate was generated from unlabeled glutamate via the aspartate transaminase reaction. A fraction of Phe and Tyr obtained label via the chorismate route from erythrose 4-phosphate, generated via the pentose phosphate pathway, and phosphoenolpyruvate. Strikingly, the data revealed no significant flux from phosphoglycerate to Ser and Gly but showed formation of Ser via the reverse reaction, namely by hydroxymethylation of Gly. The essential Gly was acquired from the medium, and the biosynthesis pathway was confirmed in experiments using [U-(13)C(2)]glycine as a tracer. The hydroxymethyl group in Ser originated from formate, which was generated by the pyruvate formate-lyase. Highly similar isotopologue profiles were observed in corresponding experiments with pneumococcal mutants deficient in PavA, CodY, and glucose-6-phosphate dehydrogenase pointing to the robustness of the core metabolic network used by these facultative pathogenic bacteria. In conclusion, this study demonstrates the dual utilization of carbohydrates and amino acids under in vitro conditions and identifies the unconventional de novo biosynthesis of serine by pneumococci.

Roles of rel(Spn) in stringent response, global regulation and virulence of serotype 2 Streptococcus pneumoniae D39

RelA/SpoT homologue (RSH) proteins have (p)ppGpp synthetase and hydrolase activities that mediate major global responses to nutrient limitation and other stresses. RSH proteins are conserved in most bacteria and play diverse roles in bacterial pathogenesis. We report here that the RSH protein of Streptococcus pneumoniae, Rel(Spn), can be deleted and is the primary source of (p)ppGpp synthesis in virulent strain D39 under some conditions. A D39 Deltarel(Spn) mutant grew well in complex medium, but did not grow in chemically defined medium unless supplemented with the metals copper and manganese. Transcriptome analysis of D39 rel(+)(Spn) and Deltarel(Spn) strains treated with mupirocin revealed rel(Spn)-independent (translation stress), rel(Spn)-dependent (stringent response) and Deltarel(Spn)-dependent changes, suggesting that rel(Spn) and (p)ppGpp amount play wide-ranging homeostatic roles in pneumococcal physiology, besides adjusting macromolecular synthesis and transport in response to nutrient availability. Notably, the rel(Spn)-dependent response included significant upregulation of the ply operon encoding pneumolysin toxin, whereas the Deltarel(Spn)-dependent response affected expression linked to the VicRK and CiaRH two-component systems. Finally, a D39 Deltarel(Spn) mutant was severely attenuated and displayed a significantly altered course of disease progression in a mouse model of infection, which was restored to normal by an ectopic copy of rel(+)(Spn).

An in silico analysis of T-box regulated genes and T-box evolution in prokaryotes, with emphasis on prediction of substrate specificity of transporters

BACKGROUND

T-box anti-termination is an elegant and sensitive mechanism by which many bacteria maintain constant levels of amino acid-charged tRNAs. The amino acid specificity of the regulatory element is related to a so-called specifier codon and can in principle be used to guide the functional annotation of the genes controlled via the T-box anti-termination mechanism.

RESULTS

Hidden Markov Models were defined to search the T-box regulatory element and were applied to all completed prokaryotic genomes. The vast majority of the genes found downstream of the retrieved elements encoded functionalities related to transport and synthesis of amino acids and the charging of tRNA. This is completely in line with findings reported in literature and with the proposed biological role of the regulatory element. For several species, the functional annotation of a large number of genes encoding proteins involved in amino acid transport could be improved significantly on basis of the amino acid specificity of the identified T-boxes. In addition, these annotations could be extrapolated to a larger number of orthologous systems in other species. Analysis of T-box distribution confirmed that the element is restricted predominantly to species of the phylum Firmicutes. Furthermore, it appeared that the distribution was highly species specific and that in the case of amino acid transport some boxes seemed to "pop-up" only recently.

CONCLUSION

We have demonstrated that the identification of the molecular specificity of a regulatory element can be of great help in solving notoriously difficult annotation issues, e.g. by defining the substrate specificity of genes encoding amino acid transporters on basis of the amino acid specificity of the regulatory T-box. Furthermore, our analysis of the species-dependency of the occurrence of specific T-boxes indicated that these regulatory elements propagate in a semi-independent way from the genes that they control.

Characterizing riboswitch function: identification of Mg2+ binding site in T box antiterminator RNA

T box bacterial genes utilize a riboswitch mechanism to regulate gene expression at the transcriptional level. Complementary base pairing of the 5'-untranslated mRNA with uncharged cognate tRNA stabilizes formation of an antiterminator element and permits complete transcription. In the absence of tRNA, a mutually exclusive RNA terminator element forms and results in transcription termination. This regulatory mechanism requires divalent metal ions at the antitermination event. The structural effects of Mg(2+) binding to antiterminator model RNA were investigated to ascertain if this requirement is due to the presence of a specific metal ion binding site in the antiterminator. Spectroscopic analysis identified the presence of a hydrated, diffuse Mg(2+) binding site. The results indicate that the mechanistic requirement for divalent metal ions is not due to Mg(2+)-induced pre-formation of a functional antiterminator receptor; rather, Mg(2+) binds in a helical region of high phylogenetic sequence conservation adjacent to the tRNA binding site.

Comparative genomic analysis of T-box regulatory systems in bacteria

T-box antitermination is one of the main mechanisms of regulation of genes involved in amino acid metabolism in Gram-positive bacteria. T-box regulatory sites consist of conserved sequence and RNA secondary structure elements. Using a set of known T-box sites, we constructed the common pattern and used it to scan available bacterial genomes. New T-boxes were found in various Gram-positive bacteria, some Gram-negative bacteria (delta-proteobacteria), and some other bacterial groups (Deinococcales/Thermales, Chloroflexi, Dictyoglomi). The majority of T-box-regulated genes encode aminoacyl-tRNA synthetases. Two other groups of T-box-regulated genes are amino acid biosynthetic genes and transporters, as well as genes with unknown function. Analysis of candidate T-box sites resulted in new functional annotations. We assigned the amino acid specificity to a large number of candidate amino acid transporters and a possible function to amino acid biosynthesis genes. We then studied the evolution of the T-boxes. Analysis of the constructed phylogenetic trees demonstrated that in addition to the normal evolution consistent with the evolution of regulated genes, T-boxes may be duplicated, transferred to other genes, and change specificity. We observed several cases of recent T-box regulon expansion following the loss of a previously existing regulatory system, in particular, arginine regulon in Clostridium difficile and methionine regulon in Lactobacillaceae. Finally, we described a new structural class of T-boxes containing duplicated terminator-antiterminator elements and unusual reduced T-boxes regulating initiation of translation in the Actinobacteria.

Identification of neomycin B-binding site in T box antiterminator model RNA

The T box transcription antitermination mechanism regulates the expression of unique genes in many Gram-positive bacteria by responding, in a magnesium-dependent manner, to uncharged cognate tRNA base pairing with an antiterminator RNA element and other regions of the 5'-untranslated region. Model T box antiterminator RNA is known to bind aminoglycosides, ligands that typically bind RNA in divalent metal ion-binding sites. In this study, enzymatic footprinting and spectroscopic assays were used to identify and characterize the binding site of neomycin B to an antiterminator model RNA. Neomycin B binds the antiterminator bulge nucleotides in an electrostatic-dependent manner and displaces 3-4 monovalent cations, indicating that the antiterminator likely contains a divalent metal ion-binding site. Neomycin B facilitates rather than inhibits tRNA binding indicating that bulge-targeted inhibitors that bind the antiterminator via non-electrostatic interactions may be the more optimal candidates for antiterminator-targeted ligand design.

A model of evolution with constant selective pressure for regulatory DNA sites

BACKGROUND

Molecular evolution is usually described assuming a neutral or weakly non-neutral substitution model. Recently, new data have become available on evolution of sequence regions under a selective pressure, e.g. transcription factor binding sites. To reconstruct the evolutionary history of such sequences, one needs evolutionary models that take into account a substantial constant selective pressure.

RESULTS

We present a simple evolutionary model with a single preferred (consensus) nucleotide and the neutral substitution model adopted for all other nucleotides. This evolutionary model has a rate matrix in which all substitutions that do not involve the consensus nucleotide occur with the same rate. The model has two time scales for achieving a stationary distribution; in the general case only one of the two rate parameters can be evaluated from the stationary distribution. In the middle-time zone, a counterintuitive behavior was observed for some parameter values, with a probability of conservation for a non-consensus nucleotide greater than that for the consensus nucleotide. Such an effect can be observed only in the case of weak preference for the consensus nucleotide, when the probability to observe the consensus nucleotide in the stationary distribution is less than 1/2. If the substitution rate is represented as a product of mutation and fixation, only the fixation can be calculated from the stationary distribution. The exhibited conservation of non-consensus nucleotides does not take place if the elements of mutation matrix are identical, and can be related to the reduced mutation rate between the non-consensus nucleotides. This bias can have no effect on the stationary distribution of nucleotide frequencies calculated over the ensemble of multiple alignments, e.g. transcription factor binding sites upstream of different sets of co-regulated orthologous genes.

CONCLUSION

The derived model can be used as a null model when analyzing the evolution of orthologous transcription factor binding sites. In particular, our findings show that a nucleotide preferred at some position of a multiple alignment of binding sites for some transcription factor in the same genome is not necessarily the most conserved nucleotide in an alignment of orthologous sites from different species. However, this effect can take place only in the case of a mutation matrix whose elements are not identical.

Modeling classic attenuation regulation of gene expression in bacteria

A model is proposed primarily for the classical RNA attenuation regulation of gene expression through premature transcription termination. The model is based on the concept of the RNA secondary structure macrostate within the regulatory region between the ribosome and RNA-polymerase, on hypothetical equation describing deceleration of RNA-polymerase by a macrostate and on views of transcription and translation initiation and elongation, under different values of the four basic model parameters which were varied. A special effort was made to select adequate model parameters. We first discuss kinetics of RNA folding and define the concept of the macrostate as a specific parentheses structure used to construct a conventional set of hairpins. The originally developed software that realizes the proposed model offers functionality to fully model RNA secondary folding kinetics. Its performance is compared to that of a public server described in Ref. 1. We then describe the delay in RNA-polymerase shifting to the next base or its premature termination caused by an RNA secondary structure or, herefrom, a macrostate. In this description, essential concepts are the basic and excited states of the polymerase first introduced in Ref. 2: the polymerase shifting to the next base can occur only in the basic state, and its detachment from DNA strand - only in excited state. As to the authors' knowledge, such a model incorporating the above-mentioned attenuation characteristics is not published elsewhere. The model was implemented in an application with command line interface for running in batch mode in Windows and Linux environments, as well as a public web server.(3) The model was tested with a conventional Monte Carlo procedure. In these simulations, the estimate of correlation between the premature transcription termination probability p and concentration c of charged amino acyl-tRNA was obtained as function p(c) for many regulatory regions in many bacterial genomes, as well as for local mutations in these regions.

Evolutionary origins of the eukaryotic shikimate pathway: gene fusions, horizontal gene transfer, and endosymbiotic replacements

Currently the shikimate pathway is reported as a metabolic feature of prokaryotes, ascomycete fungi, apicomplexans, and plants. The plant shikimate pathway enzymes have similarities to prokaryote homologues and are largely active in chloroplasts, suggesting ancestry from the plastid progenitor genome. Toxoplasma gondii, which also possesses an alga-derived plastid organelle, encodes a shikimate pathway with similarities to ascomycete genes, including a five-enzyme pentafunctional arom. These data suggests that the shikimate pathway and the pentafunctional arom either had an ancient origin in the eukaryotes or was conveyed by eukaryote-to-eukaryote horizontal gene transfer (HGT). We expand sampling and analyses of the shikimate pathway genes to include the oomycetes, ciliates, diatoms, basidiomycetes, zygomycetes, and the green and red algae. Sequencing of cDNA from Tetrahymena thermophila confirmed the presence of a pentafused arom, as in fungi and T. gondii. Phylogenies and taxon distribution suggest that the arom gene fusion event may be an ancient eukaryotic innovation. Conversely, the Plantae lineage (represented here by both Viridaeplantae and the red algae) acquired different prokaryotic genes for all seven steps of the shikimate pathway. Two of the phylogenies suggest a derivation of the Plantae genes from the cyanobacterial plastid progenitor genome, but if the full Plantae pathway was originally of cyanobacterial origin, then the five other shikimate pathway genes were obtained from a minimum of two other eubacterial genomes. Thus, the phylogenies demonstrate both separate HGTs and shared derived HGTs within the Plantae clade either by primary HGT transfer or secondarily via the plastid progenitor genome. The shared derived characters support the holophyly of the Plantae lineage and a single ancestral primary plastid endosymbiosis. Our analyses also pinpoints a minimum of 50 gene/domain loss events, demonstrating that loss and replacement events have been an important process in eukaryote genome evolution.

Biochemical and genetic characterization of an early step in a novel pathway for the biosynthesis of aromatic amino acids and p-aminobenzoic acid in the archaeon Methanococcus maripaludis

Methanococcus maripaludis is a strictly anaerobic, methane-producing archaeon and facultative autotroph capable of biosynthesizing all the amino acids and vitamins required for growth. In this work, the novel 6-deoxy-5-ketofructose-1-phosphate (DKFP) pathway for the biosynthesis of aromatic amino acids (AroAAs) and p-aminobenzoic acid (PABA) was demonstrated in M. maripaludis. Moreover, PABA was shown to be derived from an early intermediate in AroAA biosynthesis and not from chorismate. Following metabolic labelling with [U-(13)C]-acetate, the expected enrichments for phenylalanine and arylamine derived from PABA were observed. DKFP pathway activity was reduced following growth with aryl acids, an alternative source of the AroAAs. Lastly, a deletion mutant of aroA', which encodes the first step in the DKFP pathway, required AroAAs and PABA for growth. Complementation of the mutants by an aroA' expression vector restored the wild-type phenotype. In contrast, a deletion of aroB', which encodes the second step in the DKFP pathway, did not require AroAAs or PABA for growth. Presumably, methanococci contain an alternative activity for this step. These results identify the initial reactions of a new pathway for the biosynthesis of PABA in methanococci.

Different levels of transcriptional regulation due to trophic constraints in the reduced genome of Buchnera aphidicola APS

Symbiotic associations involving intracellular microorganisms and animals are widespread, especially for species feeding on poor or unbalanced diets. Buchnera aphidicola, the obligate intracellular bacterium associated with most aphid species, provides its hosts with essential amino acids (EAAs), nutrients in short supply in the plant phloem sap. The Buchnera genome has undergone severe reductions during intracellular evolution. Genes for EAA biosynthesis are conserved, but most of the transcriptional regulatory elements are lost. This work addresses two main questions: is transcription in Buchnera (i) regulated and (ii) scaled to aphid EAA demand? Two microarray experiments were designed for profiling the gene expression in Buchnera. The first one was characterized by a specific depletion of tyrosine and phenylalanine in the aphid diet, and the second experiment combined a global diminution of EAAs in the aphid diet with a sucrose concentration increase to manipulate the aphid growth rate. Aphid biological performance and budget analysis (the balance between EAAs provided by the diet and those synthesized by Buchnera) were performed to quantify the nutritional demand from the aphids toward their symbiotic bacteria. Despite the absence of known regulatory elements, a significant transcriptional regulation was observed at different levels of organization in the Buchnera genome: between genes, within putative transcription units, and within specific metabolic pathways. However, unambiguous evidence for transcriptional changes underpinning the scaling of EAA biosynthesis to aphid demand was not obtained. The phenotypic relevance of the transcriptional response from the reduced genome of Buchnera is addressed.

Evolution of transcriptional regulatory networks in microbial genomes

Advances in sequencing and generating high-throughput expression data have created a situation in which it is possible to integrate comparative analysis with genome-wide studies of the structure and function of regulatory systems in model organisms. Recent studies have focused on topological properties and the evolution of regulatory networks. This problem can be addressed on several levels: evolution of binding sites upstream of orthologous or duplicated genes; co-evolution of transcription factors and the DNA motifs that they recognize; expansion, contraction and replacement of regulatory systems; the relationship between co-regulation and co-expression; and, finally, construction of evolutionary models that generate networks with realistic properties. This should eventually lead to the creation of a theory of regulatory evolution with a similar level of detail and understanding to the theory of molecular evolution of protein and DNA sequences.

More robust detection of motifs in coexpressed genes by using phylogenetic information

Background

Several motif detection algorithms have been developed to discover overrepresented motifs in sets of coexpressed genes. However, in a noisy gene list, the number of genes containing the motif versus the number lacking the motif might not be sufficiently high to allow detection by classical motif detection tools. To still recover motifs which are not significantly enriched but still present, we developed a procedure in which we use phylogenetic footprinting to first delineate all potential motifs in each gene. Then we mutually compare all detected motifs and identify the ones that are shared by at least a few genes in the data set as potential candidates.

Results

We applied our methodology to a compiled test data set containing known regulatory motifs and to two biological data sets derived from genome wide expression studies. By executing four consecutive steps of 1) identifying conserved regions in orthologous intergenic regions, 2) aligning these conserved regions, 3) clustering the conserved regions containing similar regulatory regions followed by extraction of the regulatory motifs and 4) screening the input intergenic sequences with detected regulatory motif models, our methodology proves to be a powerful tool for detecting regulatory motifs when a low signal to noise ratio is present in the input data set. Comparing our results with two other motif detection algorithms points out the robustness of our algorithm.

Conclusion

We developed an approach that can reliably identify multiple regulatory motifs lacking a high degree of overrepresentation in a set of coexpressed genes (motifs belonging to sparsely connected hubs in the regulatory network) by exploiting the advantages of using both coexpression and phylogenetic information.

Comparative analysis of RNA regulatory elements of amino acid metabolism genes in Actinobacteria

Background

Formation of alternative structures in mRNA in response to external stimuli, either direct or mediated by proteins or other RNAs, is a major mechanism of regulation of gene expression in bacteria. This mechanism has been studied in detail using experimental and computational approaches in proteobacteria and Firmicutes, but not in other groups of bacteria.

Results

Comparative analysis of amino acid biosynthesis operons in Actinobacteria resulted in identification of conserved regions upstream of several operons. Classical attenuators were predicted upstream of trp operons in Corynebacterium spp. and Streptomyces spp., and trpS and leuS genes in some Streptomyces spp. Candidate leader peptides with terminators were observed upstream of ilvB genes in Corynebacterium spp., Mycobacterium spp. and Streptomyces spp. Candidate leader peptides without obvious terminators were found upstream of cys operons in Mycobacterium spp. and several other species. A conserved pseudoknot (named LEU element) was identified upstream of leuA operons in most Actinobacteria. Finally, T-boxes likely involved in the regulation of translation initiation were observed upstream of ileS genes from several Actinobacteria.

Conclusion

The metabolism of tryptophan, cysteine and leucine in Actinobacteria seems to be regulated on the RNA level. In some cases the mechanism is classical attenuation, but in many cases some components of attenuators are missing. The most interesting case seems to be the leuA operon preceded by the LEU element that may fold into a conserved pseudoknot or an alternative structure. A LEU element has been observed in a transposase gene from Bifidobacterium longum, but it is not conserved in genes encoding closely related transposases despite a very high level of protein similarity. One possibility is that the regulatory region of the leuA has been co-opted from some element involved in transposition. Analysis of phylogenetic patterns allowed for identification of ML1624 of M. leprae and its orthologs as the candidate regulatory proteins that may bind to the LEU element. T-boxes upstream of the ileS genes are unusual, as their regulatory mechanism seems to be inhibition of translation initiation via a hairpin sequestering the Shine-Dalgarno box.

New insights into regulation of the tryptophan biosynthetic operon in Gram-positive bacteria

The tryptophan operon of Bacillus subtilis serves as an excellent model for investigating transcription regulation in Gram-positive bacteria. In this article, we extend this knowledge by analyzing the predicted regulatory regions in the trp operons of other fully sequenced Gram-positive bacteria. Interestingly, it appears that in eight of the organisms examined, transcription of the trp operon appears to be regulated by tandem T-box elements. These regulatory elements have recently been described in the trp operons of two bacterial species. Single T-box elements are commonly found in Gram-positive bacteria in operons encoding aminoacyl tRNA synthetases and proteins performing other functions. Different regulatory mechanisms appear to be associated with variations of trp gene organization within the trp operon.

Differences in LexA regulon structure among Proteobacteria through in vivo assisted comparative genomics

The LexA regulon encompasses an ensemble of genes involved in preserving cell viability under massive DNA damage and is present in most bacterial phyla. Up to date, however, the scope of this network had only been assessed in the Gamma Proteobacteria. Here, we report the structure of the LexA regulon in the Alpha Proteobacteria, using a combined approach that makes use of in vitro and in vivo techniques to assist and validate the comparative genomics in silico methodology. This leads to the first experimentally validated description of the LexA regulon in the Alpha Proteobacteria, and comparison of regulon core structures in both classes suggests that a least common multiple set of genes (recA, ssb, uvrA and ruvCAB) might be a defining property of the Proteobacteria LexA network.

Regulog analysis: detection of conserved regulatory networks across bacteria: application to Staphylococcus aureus

A transcriptional regulatory network encompasses sets of genes (regulons) whose expression states are directly altered in response to an activating signal, mediated by trans-acting regulatory proteins and cis-acting regulatory sequences. Enumeration of these network components is an essential step toward the creation of a framework for systems-based analysis of biological processes. Profile-based methods for the detection of cis-regulatory elements are often applied to predict regulon members, but they suffer from poor specificity. In this report we describe Regulogger, a novel computational method that uses comparative genomics to eliminate spurious members of predicted gene regulons. Regulogger produces regulogs, sets of coregulated genes for which the regulatory sequence has been conserved across multiple organisms. The quantitative method assigns a confidence score to each predicted regulog member on the basis of the degree of conservation of protein sequence and regulatory mechanisms. When applied to a reference collection of regulons from Escherichia coli, Regulogger increased the specificity of predictions up to 25-fold over methods that use cis-element detection in isolation. The enhanced specificity was observed across a wide range of biologically meaningful parameter combinations, indicating a robust and broad utility for the method. The power of computational pattern discovery methods coupled with Regulogger to unravel transcriptional networks was demonstrated in an analysis of the genome of Staphylococcus aureus. A total of 125 regulogs were found in this organism, including both well-defined functional groups and a subset with unknown functions.

Two biosynthetic pathways for aromatic amino acids in the archaeon Methanococcus maripaludis

Methanococcus maripaludis is a strictly anaerobic, methane-producing archaeon. Aromatic amino acids (AroAAs) are biosynthesized in this autotroph either by the de novo pathway, with chorismate as an intermediate, or by the incorporation of exogenous aryl acids via indolepyruvate oxidoreductase (IOR). In order to evaluate the roles of these pathways, the gene that encodes the third step in the de novo pathway, 3-dehydroquinate dehydratase (DHQ), was deleted. This mutant required all three AroAAs for growth, and no DHQ activity was detectible in cell extracts, compared to 6.0 +/- 0.2 mU mg(-1) in the wild-type extract. The growth requirement for the AroAAs could be fulfilled by the corresponding aryl acids phenylacetate, indoleacetate, and p-hydroxyphenylacetate. The specific incorporation of phenylacetate into phenylalanine by the IOR pathway was demonstrated in vivo by labeling with [1-(13)C]phenylacetate. M. maripaludis has two IOR homologs. A deletion mutant for one of these homologs contained 76, 74, and 42% lower activity for phenylpyruvate, p-hydoxyphenylpyruvate, and indolepyruvate oxidation, respectively, than the wild type. Growth of this mutant in minimal medium was inhibited by the aryl acids, but the AroAAs partially restored growth. Genetic complementation of the IOR mutant also restored much of the wild-type phenotype. Thus, aryl acids appear to regulate the expression or activity of the de novo pathway. The aryl acids did not significantly inhibit the activity of the biosynthetic enzymes chorismate mutase, prephenate dehydratase, and prephenate dehydrogenase in cell extracts, so the inhibition of growth was probably not due to an effect on these enzymes.

Comparative genomics of the methionine metabolism in Gram-positive bacteria: a variety of regulatory systems

Regulation of the methionine biosynthesis and transport genes in bacteria is rather diverse and involves two RNA-level regulatory systems and at least three DNA-level systems. In particular, the methionine metabolism in Gram-positive bacteria was known to be controlled by the S-box and T-box mechanisms, both acting on the level of premature termination of transcription. Using comparative analysis of genes, operons and regulatory elements, we described the methionine metabolic pathway and the methionine regulons in available genomes of Gram-positive bacteria. A large number of methionine-specific RNA elements were identified. S-boxes were shown to be widely distributed in Bacillales and Clostridia, whereas methionine-specific T-boxes occurred mostly in Lactobacillales. A candidate binding signal (MET-box) for a hypothetical methionine regulator, possibly MtaR, was identified in Streptococcaceae, the only family in the Bacillus/Clostridium group of Gram-positive bacteria having neither S-boxes, nor methionine-specific T-boxes. Positional analysis of methionine-specific regulatory sites complemented by genome context analysis lead to identification of new members of the methionine regulon, both enzymes and transporters, and reconstruction of the methionine metabolism in various bacterial genomes. In particular, we found candidate transporters for methionine (MetT) and methylthioribose (MtnABC), as well as new enzymes forming the S-adenosylmethionine recycling pathway. Methionine biosynthetic enzymes in various bacterial species are quite variable. In particular, Oceanobacillus iheyensis possibly uses a homolog of the betaine-homocysteine methyltransferase bhmT gene from vertebrates to substitute missing bacterial-type methionine synthases.