A generally applicable validation scheme for the assessment of factors involved in reproducibility and quality of DNA-microarray data

Glutamate Dehydrogenase (GdhA) of Streptococcus pneumoniae Is Required for High Temperature Adaptation

During its progression from the nasopharynx to other sterile and nonsterile niches of its human host, Streptococcus pneumoniae must cope with changes in temperature. We hypothesized that the temperature adaptation is an important facet of pneumococcal survival in the host. Here, we evaluated the effect of temperature on pneumococcus and studied the role of glutamate dehydrogenase (GdhA) in thermal adaptation associated with virulence and survival. Microarray analysis revealed a significant transcriptional response to changes in temperature, affecting the expression of 252 genes in total at 34°C and 40°C relative to at 37°C. One of the differentially regulated genes was gdhA, which is upregulated at 40°C and downregulated at 34°C relative to 37°C. Deletion of gdhA attenuated the growth, cell size, biofilm formation, pH survival, and biosynthesis of proteins associated with virulence in a temperature-dependent manner. Moreover, deletion of gdhA stimulated formate production irrespective of temperature fluctuation. Finally, ΔgdhA grown at 40°C was less virulent than other temperatures or the wild type at the same temperature in a Galleria mellonella infection model, suggesting that GdhA is required for pneumococcal virulence at elevated temperature.

Transcriptional control of central carbon metabolic flux in Bifidobacteria by two functionally similar, yet distinct LacI-type regulators

Bifidobacteria resident in the gastrointestinal tract (GIT) are subject to constantly changing environmental conditions, which require rapid adjustments in gene expression. Here, we show that two predicted LacI-type transcription factors (TFs), designated AraQ and MalR1, are involved in regulating the central, carbohydrate-associated metabolic pathway (the so-called phosphoketolase pathway or bifid shunt) of the gut commensal Bifidobacterium breve UCC2003. These TFs appear to not only control transcription of genes involved in the bifid shunt and each other, but also seem to commonly and directly affect transcription of other TF-encoding genes, as well as genes related to uptake and metabolism of various carbohydrates. This complex and interactive network of AraQ/MalR1-mediated gene regulation provides previously unknown insights into the governance of carbon metabolism in bifidobacteria.

Metabolism of the predominant human milk oligosaccharide fucosyllactose by an infant gut commensal

A number of bifidobacterial species are found at a particularly high prevalence and abundance in faecal samples of healthy breastfed infants, a phenomenon that is believed to be, at least partially, due to the ability of bifidobacteria to metabolize Human Milk Oligosaccharides (HMOs). In the current study, we isolated a novel strain of Bifidobacterium kashiwanohense, named APCKJ1, from the faeces of a four-week old breastfed infant, based on the ability of the strain to utilise the HMO component fucosyllactose. We then determined the full genome sequence of this strain, and employed the generated data to analyze fucosyllactose metabolism in B. kashiwanohense APCKJ1. Transcriptomic and growth analyses, combined with metabolite analysis, in vitro hydrolysis assays and heterologous expression, allowed us to elucidate the pathway for fucosyllactose metabolism in B. kashiwanohense APCKJ1. Homologs of the key genes for this metabolic pathway were identified in particular in infant-derived members of the Bifdobacterium genus, revealing the apparent niche-specific nature of this pathway, and allowing a broad perspective on bifidobacterial fucosyllactose and L-fucose metabolism.

Comparative genome and methylome analysis reveals restriction/modification system diversity in the gut commensal Bifidobacterium breve

Bifidobacterium breve represents one of the most abundant bifidobacterial species in the gastro-intestinal tract of breast-fed infants, where their presence is believed to exert beneficial effects. In the present study whole genome sequencing, employing the PacBio Single Molecule, Real-Time (SMRT) sequencing platform, combined with comparative genome analysis allowed the most extensive genetic investigation of this taxon. Our findings demonstrate that genes encoding Restriction/Modification (R/M) systems constitute a substantial part of the B. breve variable gene content (or variome). Using the methylome data generated by SMRT sequencing, combined with targeted Illumina bisulfite sequencing (BS-seq) and comparative genome analysis, we were able to detect methylation recognition motifs and assign these to identified B. breve R/M systems, where in several cases such assignments were confirmed by restriction analysis. Furthermore, we show that R/M systems typically impose a very significant barrier to genetic accessibility of B. breve strains, and that cloning of a methyltransferase-encoding gene may overcome such a barrier, thus allowing future functional investigations of members of this species.

Reconstruction and inference of the Lactococcus lactis MG1363 gene co-expression network

Lactic acid bacteria are Gram-positive bacteria used throughout the world in many industrial applications for their acidification, flavor and texture formation attributes. One of the species, Lactococcus lactis, is employed for the production of fermented milk products like cheese, buttermilk and quark. It ferments lactose to lactic acid and, thus, helps improve the shelf life of the products. Many physiological and transcriptome studies have been performed in L. lactis in order to comprehend and improve its biotechnological assets. Using large amounts of transcriptome data to understand and predict the behavior of biological processes in bacterial or other cell types is a complex task. Gene networks enable predicting gene behavior and function in the context of transcriptionally linked processes. We reconstruct and present the gene co-expression network (GCN) for the most widely studied L. lactis strain, MG1363, using publicly available transcriptome data. Several methods exist to generate and judge the quality of GCNs. Different reconstruction methods lead to networks with varying structural properties, consequently altering gene clusters. We compared the structural properties of the MG1363 GCNs generated by five methods, namely Pearson correlation, Spearman correlation, GeneNet, Weighted Gene Co-expression Network Analysis (WGCNA), and Sparse PArtial Correlation Estimation (SPACE). Using SPACE, we generated an L. lactis MG1363 GCN and assessed its quality using modularity and structural and biological criteria. The L. lactis MG1363 GCN has structural properties similar to those of the gold-standard networks of Escherichia coli K-12 and Bacillus subtilis 168. We showcase that the network can be used to mine for genes with similar expression profiles that are also generally linked to the same biological process.

Adaption to glucose limitation is modulated by the pleotropic regulator CcpA, independent of selection pressure strength

Background

A central theme in (micro)biology is understanding the molecular basis of fitness i.e. which strategies are successful under which conditions; how do organisms implement such strategies at the molecular level; and which constraints shape the trade-offs between alternative strategies. Highly standardized microbial laboratory evolution experiments are ideally suited to approach these questions. For example, prolonged chemostats provide a constant environment in which the growth rate can be set, and the adaptive process of the organism to such environment can be subsequently characterized.

Results

We performed parallel laboratory evolution of Lactococcus lactis in chemostats varying the quantitative value of the selective pressure by imposing two different growth rates. A mutation in one specific amino acid residue of the global transcriptional regulator of carbon metabolism, CcpA, was selected in all of the evolution experiments performed. We subsequently showed that this mutation confers predictable fitness improvements at other glucose-limited growth rates as well. In silico protein structural analysis of wild type and evolved CcpA, as well as biochemical and phenotypic assays, provided the underpinning molecular mechanisms that resulted in the specific reprogramming favored in constant environments.

Conclusion

This study provides a comprehensive understanding of a case of microbial evolution and hints at the wide dynamic range that a single fitness-enhancing mutation may display. It demonstrates how the modulation of a pleiotropic regulator can be used by cells to improve one trait while simultaneously work around other limiting constraints, by fine-tuning the expression of a wide range of cellular processes.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1331-x) contains supplementary material, which is available to authorized users.

Staying alive: growth and survival of Bifidobacterium animalis subsp. animalis under in vitro and in vivo conditions

Members of the Bifidobacterium genus are widely used as probiotics in fermented milk products. Bifidobacterium animalis subsp. animalis CNCM I-4602 grows and survives poorly in reconstituted skimmed milk (RSM). Availing of genome and transcriptome information, this poor growth and survival phenotype in milk was substantially improved by the addition of certain compounds, such as yeast extract, uric acid, glutathione, cysteine, ferrous sulfate, and a combination of magnesium sulfate and manganese sulfate. Carbohydrate utilization of CNCM I-4602 was also investigated, allowing the identification of several carbohydrate utilization gene clusters, and highlighting this strain's inability to utilize lactose, unlike the type strain of this subspecies, B. animalis subsp. animalis ATCC25527 and the B. animalis subsp. lactis subspecies. In addition, the ability of B. animalis subsp. animalis CNCM I-4602 to colonize a murine model was investigated, which showed that this strain persists in the murine gut for a period of at least 4 weeks. Associated in vivo transcriptome analysis revealed that, among other genes, a gene cluster encoding a predicted type IVb tight adherence (Tad) pilus was upregulated, indicating that this extracellular structure plays a role in the colonization/adaptation of the murine gastrointestinal tract by this strain.

Bifidobacterium breve UCC2003 Employs Multiple Transcriptional Regulators To Control Metabolism of Particular Human Milk Oligosaccharides

Bifidobacterial carbohydrate metabolism has been studied in considerable detail for a variety of both plant- and human-derived glycans, particularly involving the bifidobacterial prototype strain Bifidobacterium breve UCC2003. We recently elucidated the metabolic pathways by which the human milk oligosaccharide (HMO) constituents lacto-N-tetraose (LNT), lacto-N-neotetraose (LNnT) and lacto-N-biose (LNB) are utilized by B. breve UCC2003. However, to date, no work has been carried out on the regulatory mechanisms that control the expression of the genetic loci involved in these HMO metabolic pathways. In this study, we describe the characterization of three transcriptional regulators and the corresponding operator and associated (inducible) promoter sequences, with the latter governing the transcription of the genetic elements involved in LN(n)T/LNB metabolism. The activity of these regulators is dependent on the release of specific monosaccharides, which are believed to act as allosteric effectors and which are derived from the corresponding HMOs targeted by the particular locus.IMPORTANCE Human milk oligosaccharides (HMOs) are a key factor in the development of the breastfed-infant microbiota. They function as prebiotics, selecting for a specific range of microbes, including a number of infant-associated species of bifidobacteria, which are thought to provide a range of health benefits to the infant host. While much research has been carried out on elucidating the mechanisms of HMO metabolism in infant-associated bifidobacteria, to date there is very little understanding of the transcriptional regulation of these pathways. This study reveals a multicomponent transcriptional regulation system that controls the recently identified pathways of HMO metabolism in the infant-associated Bifidobacterium breve prototype strain UCC2003. This not only provides insight into the regulatory mechanisms present in other infant-associated bifidobacteria but also provides an example of a network of sequential steps regulating microbial carbohydrate metabolism.

Characterizing a thermostable Cas9 for bacterial genome editing and silencing

CRISPR-Cas9-based genome engineering tools have revolutionized fundamental research and biotechnological exploitation of both eukaryotes and prokaryotes. However, the mesophilic nature of the established Cas9 systems does not allow for applications that require enhanced stability, including engineering at elevated temperatures. Here we identify and characterize ThermoCas9 from the thermophilic bacterium Geobacillus thermodenitrificans T12. We show that in vitro ThermoCas9 is active between 20 and 70 °C, has stringent PAM-preference at lower temperatures, tolerates fewer spacer-protospacer mismatches than SpCas9 and its activity at elevated temperatures depends on the sgRNA-structure. We develop ThermoCas9-based engineering tools for gene deletion and transcriptional silencing at 55 °C in Bacillus smithii and for gene deletion at 37 °C in Pseudomonas putida. Altogether, our findings provide fundamental insights into a thermophilic CRISPR-Cas family member and establish a Cas9-based bacterial genome editing and silencing tool with a broad temperature range.

Genome-Wide Search for Genes Required for Bifidobacterial Growth under Iron-Limitation

Bacteria evolved over millennia in the presence of the vital micronutrient iron. Iron is involved in numerous processes within the cell and is essential for nearly all living organisms. The importance of iron to the survival of bacteria is obvious from the large variety of mechanisms by which iron may be acquired from the environment. Random mutagenesis and global gene expression profiling led to the identification of a number of genes, which are essential for Bifidobacterium breve UCC2003 survival under iron-restrictive conditions. These genes encode, among others, Fe-S cluster-associated proteins, a possible ferric iron reductase, a number of cell wall-associated proteins, and various DNA replication and repair proteins. In addition, our study identified several presumed iron uptake systems which were shown to be essential for B. breve UCC2003 growth under conditions of either ferric and/or ferrous iron chelation. Of these, two gene clusters encoding putative iron-uptake systems, bfeUO and sifABCDE, were further characterised, indicating that sifABCDE is involved in ferrous iron transport, while the bfeUO-encoded transport system imports both ferrous and ferric iron. Transcription studies showed that bfeUO and sifABCDE constitute two separate transcriptional units that are induced upon dipyridyl-mediated iron limitation. In the anaerobic gastrointestinal environment ferrous iron is presumed to be of most relevance, though a mutation in the sifABCDE cluster does not affect B. breve UCC2003's ability to colonise the gut of a murine model.

CodY Regulates Thiol Peroxidase Expression as Part of the Pneumococcal Defense Mechanism against H2O2 Stress

Streptococcus pneumoniae is a facultative anaerobic pathogen. Although it maintains fermentative metabolism, during aerobic growth pneumococci produce high levels of H₂O₂, which can have adverse effects on cell viability and DNA, and influence pneumococcal interaction with its host. The pneumococcus is unusual in its dealing with toxic reactive oxygen species (ROS) in that it neither has catalase nor the global regulators of peroxide stress resistance. Previously, we identified pneumococcal thiol peroxidase (TpxD) as the key enzyme for enzymatic removal of H₂O₂, and showed that TpxD synthesis is up-regulated upon exposure to H₂O₂. This study aimed to reveal the mechanism controlling TpxD expression under H₂O₂ stress. We hypothesize that H₂O₂ activates a transcription factor which in turn up-regulates tpxD expression. Microarray analysis revealed a pneumococcal global transcriptional response to H₂O₂. Mutation of tpxD abolished H₂O₂-mediated response to high H₂O₂ levels, signifying the need for an active TpxD under oxidative stress conditions. Bioinformatic tools, applied to search for a transcription factor modulating tpxD expression, pointed toward CodY as a potential candidate. Indeed, a putative 15-bp consensus CodY binding site was found in the proximal region of tpxD-coding sequence. Binding of CodY to this site was confirmed by EMSA, and genetic engineering techniques demonstrated that this site is essential for TpxD up-regulation under H₂O₂ stress. Furthermore, tpxD expression was reduced in a ΔcodY mutant. These data indicate that CodY is an activator of tpxD expression, triggering its up-regulation under H₂O₂ stress. In addition we show that H₂O₂ specifically oxidizes the 2 CodY cysteines. This oxidation may trigger a conformational change in CodY, resulting in enhanced binding to DNA. A schematic model illustrating the contribution of TpxD and CodY to pneumococcal global transcriptional response to H₂O₂ is proposed.

Efficient Genome Editing of a Facultative Thermophile Using Mesophilic spCas9

Well-developed genetic tools for thermophilic microorganisms are scarce, despite their industrial and scientific relevance. Whereas highly efficient CRISPR/Cas9-based genome editing is on the rise in prokaryotes, it has never been employed in a thermophile. Here, we apply Streptococcus pyogenes Cas9 (spCas9)-based genome editing to a moderate thermophile, i.e., Bacillus smithii, including a gene deletion, gene knockout via insertion of premature stop codons, and gene insertion. We show that spCas9 is inactive in vivo above 42 °C, and we employ the wide temperature growth range of B. smithii as an induction system for spCas9 expression. Homologous recombination with plasmid-borne editing templates is performed at 45-55 °C, when spCas9 is inactive. Subsequent transfer to 37 °C allows for counterselection through production of active spCas9, which introduces lethal double-stranded DNA breaks to the nonedited cells. The developed method takes 4 days with 90, 100, and 20% efficiencies for gene deletion, knockout, and insertion, respectively. The major advantage of our system is the limited requirement for genetic parts: only one plasmid, one selectable marker, and a promoter are needed, and the promoter does not need to be inducible or well-characterized. Hence, it can be easily applied for genome editing purposes in both mesophilic and thermophilic nonmodel organisms with a limited genetic toolbox and ability to grow at, or tolerate, temperatures of 37 and at or above 42 °C.

Concordance analysis of microarray studies identifies representative gene expression changes in Parkinson's disease: a comparison of 33 human and animal studies

BACKGROUND

As the popularity of transcriptomic analysis has grown, the reported lack of concordance between different studies of the same condition has become a growing concern, raising questions as to the representativeness of different study types, such as non-human disease models or studies of surrogate tissues, to gene expression in the human condition.

METHODS

In a comparison of 33 microarray studies of Parkinson's disease, correlation and clustering analyses were used to determine the factors influencing concordance between studies, including agreement between different tissue types, different microarray platforms, and between neurotoxic and genetic disease models and human Parkinson's disease.

RESULTS

Concordance over all studies is low, with correlation of only 0.05 between differential gene expression signatures on average, but increases within human patients and studies of the same tissue type, rising to 0.38 for studies of human substantia nigra. Agreement of animal models, however, is dependent on model type. Studies of brain tissue from Parkinson's disease patients (specifically the substantia nigra) form a distinct group, showing patterns of differential gene expression noticeably different from that in non-brain tissues and animal models of Parkinson's disease; while comparison with other brain diseases (Alzheimer's disease and brain cancer) suggests that the mixed study types display a general signal of neurodegenerative disease. A meta-analysis of these 33 microarray studies demonstrates the greater ability of studies in humans and highly-affected tissues to identify genes previously known to be associated with Parkinson's disease.

CONCLUSIONS

The observed clustering and concordance results suggest the existence of a 'characteristic' signal of Parkinson's disease found in significantly affected human tissues in humans. These results help to account for the consistency (or lack thereof) so far observed in microarray studies of Parkinson's disease, and act as a guide to the selection of transcriptomic studies most representative of the underlying gene expression changes in the human disease.

Impact of aspirin on the transcriptome of Streptococcus pneumoniae D39

Aspirin or acetylsalicylic acid (ASA) is a medicine used to treat pain, fever, and inflammation. Here, we for the very first time reported the genome-wide transcriptional profiling of aspirin-regulated genes in Streptococcus pneumoniae in the presence of 5 mM aspirin in chemically-defined medium (CDM) using microarray analysis. Our results showed that expression of several genes was differentially expressed in the presence of aspirin. These genes were further grouped into COG (Clusters of Orthologous Groups) functional categories based on the putative functions of the corresponding proteins. Most of affected genes belong to COG category E (Amino acid transport and metabolism), G (Carbohydrate transport and metabolism), J (Translation, ribosomal structure and biogenesis), and I (Lipid transport and metabolism). Transcriptional profiling data of aspirin-regulated genes was deposited to Gene Expression Omnibus (GEO) database under accession number GSE94514.

Bifidobacterium breve UCC2003 metabolises the human milk oligosaccharides lacto-N-tetraose and lacto-N-neo-tetraose through overlapping, yet distinct pathways

In this study, we demonstrate that the prototype B. breve strain UCC2003 possesses specific metabolic pathways for the utilisation of lacto-N-tetraose (LNT) and lacto-N-neotetraose (LNnT), which represent the central moieties of Type I and Type II human milk oligosaccharides (HMOs), respectively. Using a combination of experimental approaches, the enzymatic machinery involved in the metabolism of LNT and LNnT was identified and characterised. Homologs of the key genetic loci involved in the utilisation of these HMO substrates were identified in B. breve, B. bifidum, B. longum subsp. infantis and B. longum subsp. longum using bioinformatic analyses, and were shown to be variably present among other members of the Bifidobacterium genus, with a distinct pattern of conservation among human-associated bifidobacterial species.

Glycosulfatase-Encoding Gene Cluster in Bifidobacterium breve UCC2003

Bifidobacteria constitute a specific group of commensal bacteria typically found in the gastrointestinal tract (GIT) of humans and other mammals. Bifidobacterium breve strains are numerically prevalent among the gut microbiota of many healthy breastfed infants. In the present study, we investigated glycosulfatase activity in a bacterial isolate from a nursling stool sample, B. breve UCC2003. Two putative sulfatases were identified on the genome of B. breve UCC2003. The sulfated monosaccharide N-acetylglucosamine-6-sulfate (GlcNAc-6-S) was shown to support the growth of B. breve UCC2003, while N-acetylglucosamine-3-sulfate, N-acetylgalactosamine-3-sulfate, and N-acetylgalactosamine-6-sulfate did not support appreciable growth. By using a combination of transcriptomic and functional genomic approaches, a gene cluster designated ats2 was shown to be specifically required for GlcNAc-6-S metabolism. Transcription of the ats2 cluster is regulated by a repressor open reading frame kinase (ROK) family transcriptional repressor. This study represents the first description of glycosulfatase activity within the Bifidobacterium genus.

IMPORTANCE

Bifidobacteria are saccharolytic organisms naturally found in the digestive tract of mammals and insects. Bifidobacterium breve strains utilize a variety of plant- and host-derived carbohydrates that allow them to be present as prominent members of the infant gut microbiota as well as being present in the gastrointestinal tract of adults. In this study, we introduce a previously unexplored area of carbohydrate metabolism in bifidobacteria, namely, the metabolism of sulfated carbohydrates. B. breve UCC2003 was shown to metabolize N-acetylglucosamine-6-sulfate (GlcNAc-6-S) through one of two sulfatase-encoding gene clusters identified on its genome. GlcNAc-6-S can be found in terminal or branched positions of mucin oligosaccharides, the glycoprotein component of the mucous layer that covers the digestive tract. The results of this study provide further evidence of the ability of this species to utilize mucin-derived sugars, a trait which may provide a competitive advantage in both the infant gut and adult gut.

Regulation of Cell Wall Plasticity by Nucleotide Metabolism in Lactococcus lactis

To ensure optimal cell growth and separation and to adapt to environmental parameters, bacteria have to maintain a balance between cell wall (CW) rigidity and flexibility. This can be achieved by a concerted action of peptidoglycan (PG) hydrolases and PG-synthesizing/modifying enzymes. In a search for new regulatory mechanisms responsible for the maintenance of this equilibrium in Lactococcus lactis, we isolated mutants that are resistant to the PG hydrolase lysozyme. We found that 14% of the causative mutations were mapped in the guaA gene, the product of which is involved in purine metabolism. Genetic and transcriptional analyses combined with PG structure determination of the guaA mutant enabled us to reveal the pivotal role of the pyrB gene in the regulation of CW rigidity. Our results indicate that conversion of l-aspartate (l-Asp) to N-carbamoyl-l-aspartate by PyrB may reduce the amount of l-Asp available for PG synthesis and thus cause the appearance of Asp/Asn-less stem peptides in PG. Such stem peptides do not form PG cross-bridges, resulting in a decrease in PG cross-linking and, consequently, reduced PG thickness and rigidity. We hypothesize that the concurrent utilization of l-Asp for pyrimidine and PG synthesis may be part of the regulatory scheme, ensuring CW flexibility during exponential growth and rigidity in stationary phase. The fact that l-Asp availability is dependent on nucleotide metabolism, which is tightly regulated in accordance with the growth rate, provides L. lactis cells the means to ensure optimal CW plasticity without the need to control the expression of PG synthesis genes.

Transcriptional Profile of Bacillus subtilis sigF-Mutant during Vegetative Growth

Sigma factor F is the first forespore specific transcription factor in Bacillus subtilis and controls genes required for the early stages of prespore development. The role of sigF is well studied under conditions that induce sporulation. Here, the impact of sigF disruption on the transcriptome of exponentially growing cultures is studied by micro-array analysis. Under these conditions that typically don’t induce sporulation, the transcriptome showed minor signs of sporulation initiation. The number of genes differentially expressed and the magnitude of expression were, as expected, quite small in comparison with sporulation conditions. The genes mildly down-regulated were mostly involved in anabolism and the genes mildly up-regulated, in particular fatty acid degradation genes, were mostly involved in catabolism. This is probably related to the arrest at sporulation stage II occurring in the sigF mutant, because continuation of growth from the formed disporic sporangia may require additional energy. The obtained knowledge is relevant for various experiments, such as industrial fermentation, prolonged experimental evolution or zero-growth studies, where sporulation is an undesirable trait that should be avoided, e.g by a sigF mutation.

Transcriptional and metabolic effects of glucose on Streptococcus pneumoniae sugar metabolism

Streptococcus pneumoniae is a strictly fermentative human pathogen that relies on carbohydrate metabolism to generate energy for growth. The nasopharynx colonized by the bacterium is poor in free sugars, but mucosa lining glycans can provide a source of sugar. In blood and inflamed tissues glucose is the prevailing sugar. As a result during progression from colonization to disease S. pneumoniae has to cope with a pronounced shift in carbohydrate nature and availability. Thus, we set out to assess the pneumococcal response to sugars found in glycans and the influence of glucose (Glc) on this response at the transcriptional, physiological, and metabolic levels. Galactose (Gal), N-acetylglucosamine (GlcNAc), and mannose (Man) affected the expression of 8 to 14% of the genes covering cellular functions including central carbon metabolism and virulence. The pattern of end-products as monitored by in vivo¹³C-NMR is in good agreement with the fermentation profiles during growth, while the pools of phosphorylated metabolites are consistent with the type of fermentation observed (homolactic vs. mixed) and regulation at the metabolic level. Furthermore, the accumulation of α-Gal6P and Man6P indicate metabolic bottlenecks in the metabolism of Gal and Man, respectively. Glc added to cells actively metabolizing other sugar(s) was readily consumed and elicited a metabolic shift toward a homolactic profile. The transcriptional response to Glc was large (over 5% of the genome). In central carbon metabolism (most represented category), Glc exerted mostly negative regulation. The smallest response to Glc was observed on a sugar mix, suggesting that exposure to varied sugars improves the fitness of S. pneumoniae. The expression of virulence factors was negatively controlled by Glc in a sugar-dependent manner. Overall, our results shed new light on the link between carbohydrate metabolism, adaptation to host niches and virulence.

Probing the regulatory effects of specific mutations in three major binding domains of the pleiotropic regulator CcpA of Bacillus subtilis

Carbon catabolite control is required for efficient use of available carbon sources to ensure rapid growth of bacteria. CcpA is a global regulator of carbon metabolism in Gram-positive bacteria like Bacillus subtilis. In this study the genome-wide gene regulation of a CcpA knockout and three specific CcpA mutants were studied by transcriptome analysis, to further elucidate the function of specific binding sites in CcpA. The following three amino acids were mutated to characterize their function: M17(R) which is involved in DNA binding, T62(H) which is important for the allosteric switch in CcpA upon HPr-Ser46-P binding, and R304(W) which is important for binding of the coeffectors HPr-Ser46-P and fructose-1,6-bisphosphate. The results confirm that CcpA was also involved in gene regulation in the absence of glucose. CcpA-M17R showed a small relief of Carbon Catabolite Control; the CcpA-M17R mutant regulates fewer genes than the CcpA-wt and the palindromicity of the cre site is less important for CcpA-M17R. CcpA-T62H was a stronger repressor than CcpA-wt and also acted as a strong repressor in the absence of glucose. CcpA-R304W was shown here to be less dependent on HPr-Ser46-P for its carbon catabolite control activities. The results presented here provide detailed information on alterations in gene regulation for each CcpA-mutant.

Transcriptome analysis of Streptococcus pneumoniae D39 in the presence of cobalt

Cobalt (Co(2 +)) is an important transition metal ion that plays a vital role in cellular physiology of bacteria. The role of Co(2 +) in the regulation of several genes/operons in Streptococcus pneumoniae has recently been reported [1]. The data described in this article relate to the genome-wide transcriptional profiling of Streptococcus pneumoniae D39, either in the presence or absence of 0.5 mM Co(2 +) in chemically defined medium (CDM) using DNA microarray analysis. Genes belonging to a broad range of cellular processes such as virulence, transport and efflux systems, stress response and surface attachment were differentially expressed in the presence of Co(2 +). We used transcriptional lacZ assays and electrophoretic mobility shift assays (EMSAs) to confirm our results [1]. The dataset is publicly available at the Gene Expression Omnibus (GEO) repository (http://www.ncbi.nlm.nih.gov/geo/) with accession number GSE57696.

A Duo of Potassium-Responsive Histidine Kinases Govern the Multicellular Destiny of Bacillus subtilis

Multicellular biofilm formation and surface motility are bacterial behaviors considered mutually exclusive. However, the basic decision to move over or stay attached to a surface is poorly understood. Here, we discover that in Bacillus subtilis, the key root biofilm-controlling transcription factor Spo0A~Pi (phosphorylated Spo0A) governs the flagellum-independent mechanism of social sliding motility. A Spo0A-deficient strain was totally unable to slide and colonize plant roots, evidencing the important role that sliding might play in natural settings. Microarray experiments plus subsequent genetic characterization showed that the machineries of sliding and biofilm formation share the same main components (i.e., surfactin, the hydrophobin BslA, exopolysaccharide, and de novo-formed fatty acids). Sliding proficiency was transduced by the Spo0A-phosphorelay histidine kinases KinB and KinC. We discovered that potassium, a previously known inhibitor of KinC-dependent biofilm formation, is the specific sliding-activating signal through a thus-far-unnoticed cytosolic domain of KinB, which resembles the selectivity filter sequence of potassium channels. The differential expression of the Spo0A~Pi reporter abrB gene and the different levels of the constitutively active form of Spo0A, Sad67, in Δspo0A cells grown in optimized media that simultaneously stimulate motile and sessile behaviors uncover the spatiotemporal response of KinB and KinC to potassium and the gradual increase in Spo0A~Pi that orchestrates the sequential activation of sliding, followed by sessile biofilm formation and finally sporulation in the same population. Overall, these results provide insights into how multicellular behaviors formerly believed to be antagonistic are coordinately activated in benefit of the bacterium and its interaction with the host.

IMPORTANCE

Alternation between motile and sessile behaviors is central to bacterial adaptation, survival, and colonization. However, how is the collective decision to move over or stay attached to a surface controlled? Here, we use the model plant-beneficial bacterium Bacillus subtilis to answer this question. Remarkably, we discover that sessile biofilm formation and social sliding motility share the same structural components and the Spo0A regulatory network via sensor kinases, KinB and KinC. Potassium, an inhibitor of KinC-dependent biofilm formation, triggers sliding via a potassium-perceiving cytosolic domain of KinB that resembles the selectivity filter of potassium channels. The spatiotemporal response of these kinases to variable potassium levels and the gradual increase in Spo0A~Pi levels that orchestrates the activation of sliding before biofilm formation shed light on how multicellular behaviors formerly believed to be antagonistic work together to benefit the population fitness.

The ParB-parS Chromosome Segregation System Modulates Competence Development in Streptococcus pneumoniae

ParB proteins bind centromere-like DNA sequences called parS sites and are involved in plasmid and chromosome segregation in bacteria. We previously showed that the opportunistic human pathogen Streptococcus pneumoniae contains four parS sequences located close to the origin of replication which are bound by ParB. Using chromatin immunoprecipitation (ChIP), we found here that ParB spreads out from one of these parS sites, parS(-1.6°), for more than 5 kb and occupies the nearby comCDE operon, which drives competence development. Competence allows S. pneumoniae to take up DNA from its environment, thereby mediating horizontal gene transfer, and is also employed as a general stress response. Mutating parS(-1.6°) or deleting parB resulted in transcriptional up-regulation of comCDE and ssbB (a gene belonging to the competence regulon), demonstrating that ParB acts as a repressor of competence. However, genome-wide transcription analysis showed that ParB is not a global transcriptional regulator. Different factors, such as the composition of the growth medium and antibiotic-induced stress, can trigger the sensitive switch driving competence. This work shows that the ParB-parS chromosome segregation machinery also influences this developmental process.

IMPORTANCE

Streptococcus pneumoniae (pneumococcus) is an important human pathogen responsible for more than a million deaths each year. Like all other organisms, S. pneumoniae must be able to segregate its chromosomes properly. Not only is understanding the molecular mechanisms underlying chromosome segregation in S. pneumoniae therefore of fundamental importance, but also, this knowledge might offer new leads for ways to target this pathogen. Here, we identified a link between the pneumococcal chromosome segregation system and the competence-developmental system. Competence allows S. pneumoniae to take up and integrate exogenous DNA in its chromosome. This process plays a crucial role in successful adaptation to--and escape from--host defenses, antibiotic treatments, and vaccination strategies. We show that the chromosome segregation protein ParB acts as a repressor of competence. To the best of our knowledge, this is the first example of a ParB protein controlling bacterial competence.

Spatio-temporal remodeling of functional membrane microdomains organizes the signaling networks of a bacterium

Lipid rafts are membrane microdomains specialized in the regulation of numerous cellular processes related to membrane organization, as diverse as signal transduction, protein sorting, membrane trafficking or pathogen invasion. It has been proposed that this functional diversity would require a heterogeneous population of raft domains with varying compositions. However, a mechanism for such diversification is not known. We recently discovered that bacterial membranes organize their signal transduction pathways in functional membrane microdomains (FMMs) that are structurally and functionally similar to the eukaryotic lipid rafts. In this report, we took advantage of the tractability of the prokaryotic model Bacillus subtilis to provide evidence for the coexistence of two distinct families of FMMs in bacterial membranes, displaying a distinctive distribution of proteins specialized in different biological processes. One family of microdomains harbors the scaffolding flotillin protein FloA that selectively tethers proteins specialized in regulating cell envelope turnover and primary metabolism. A second population of microdomains containing the two scaffolding flotillins, FloA and FloT, arises exclusively at later stages of cell growth and specializes in adaptation of cells to stationary phase. Importantly, the diversification of membrane microdomains does not occur arbitrarily. We discovered that bacterial cells control the spatio-temporal remodeling of microdomains by restricting the activation of FloT expression to stationary phase. This regulation ensures a sequential assembly of functionally specialized membrane microdomains to strategically organize signaling networks at the right time during the lifespan of a bacterium.

Cellular membranes organize proteins related to signal transduction, protein sorting and membrane trafficking into the so-called lipid rafts. It has been proposed that the functional diversity of lipid rafts would require a heterogeneous population of raft domains with varying compositions. However, a mechanism for such diversification is not known due in part to the complexity that entails the manipulation of eukaryotic cells. The recent discovery that bacteria organize many cellular processes in membrane microdomains (FMMs), functionally similar to the eukaryotic lipid rafts, prompted us to explore FMMs diversity in the bacterial model Bacillus subtilis. We show that diversification of FMMs occurs in cells and gives rise to functionally distinct microdomains, which compartmentalize distinct signal transduction pathways and regulate the expression of different genetic programs. We discovered that FMMs diversification does not occur randomly. Cells sequentially regulate the specialization of the FMMs during cell growth to ensure an effective and diverse activation of signaling processes.

Transcriptional profiling of UlaR-regulated genes in Streptococcus pneumoniae

The transcriptional regulator UlaR belongs to the family of PRD-containing transcriptional regulators, which are mostly involved in the regulation of carbohydrate metabolism. The role of the transcriptional regulator UlaR in Streptococcus pneumoniae has recently been described [1]. Here, we report detailed genome-wide transcriptional profiling of UlaR-regulated genes in S. pneumoniae D39 and its ∆ulaR derivative, either in the presence of 10 mM ascorbic acid in M17 medium using microarray analysis. 10 mM concentration of ascorbic acid was supplemented to the M17 medium because our lacZ-fusion studies indicated that UlaR acts as a transcriptional activator of its targets in the presence of ascorbic acid and the expression of the ula operon was maximal at a 10 mM ascorbic acid concentration [1]. All transcriptional profiling data of UlaR-regulated genes was deposited to Gene Expression Omnibus (GEO) database under accession number GSE61649.

Ascorbic acid-dependent gene expression in Streptococcus pneumoniae and the activator function of the transcriptional regulator UlaR2

In this study, we have explored the impact of ascorbic acid on the transcriptome of Streptococcus pneumoniae D39. The expression of several genes and operons, including the ula operon (which has been previously shown to be involved in ascorbic acid utilization), the AdcR regulon (which has been previously shown to be involved in zinc transport and virulence) and a PTS operon (which we denote here as ula2 operon) were altered in the presence of ascorbic acid. The ula2 operon consists of five genes, including the transcriptional activator ulaR2. Our β-galactosidase assay data and transcriptome comparison of the ulaR2 mutant with the wild-type demonstrated that the transcriptional activator UlaR2 in the presence of ascorbic acid activates the expression of the ula2 operon. We further predict a 16-bp regulatory site (5′-ATATTGTGCTCAAATA-3′) for UlaR2 in the Pula2. Furthermore, we have explored the effect of ascorbic acid on the expression of the AdcR regulon. Our ICP-MS analysis showed that addition of ascorbic acid to the medium causes zinc starvation in the cell which leads to the activation of the AdcR regulon.

Transcription of two adjacent carbohydrate utilization gene clusters in Bifidobacterium breve UCC2003 is controlled by LacI- and repressor open reading frame kinase (ROK)-type regulators

Members of the genus Bifidobacterium are commonly found in the gastrointestinal tracts of mammals, including humans, where their growth is presumed to be dependent on various diet- and/or host-derived carbohydrates. To understand transcriptional control of bifidobacterial carbohydrate metabolism, we investigated two genetic carbohydrate utilization clusters dedicated to the metabolism of raffinose-type sugars and melezitose. Transcriptomic and gene inactivation approaches revealed that the raffinose utilization system is positively regulated by an activator protein, designated RafR. The gene cluster associated with melezitose metabolism was shown to be subject to direct negative control by a LacI-type transcriptional regulator, designated MelR1, in addition to apparent indirect negative control by means of a second LacI-type regulator, MelR2. In silico analysis, DNA-protein interaction, and primer extension studies revealed the MelR1 and MelR2 operator sequences, each of which is positioned just upstream of or overlapping the correspondingly regulated promoter sequences. Similar analyses identified the RafR binding operator sequence located upstream of the rafB promoter. This study indicates that transcriptional control of gene clusters involved in carbohydrate metabolism in bifidobacteria is subject to conserved regulatory systems, representing either positive or negative control.

Physiological and cell morphology adaptation of Bacillus subtilis at near-zero specific growth rates: a transcriptome analysis

Nutrient scarcity is a common condition in nature, but the resulting extremely low growth rates (below 0.025 h(-1) ) are an unexplored research area in Bacillus subtilis. To understand microbial life in natural environments, studying the adaptation of B. subtilis to near-zero growth conditions is relevant. To this end, a chemostat modified for culturing an asporogenous B. subtilis sigF mutant strain at extremely low growth rates (also named a retentostat) was set up, and biomass accumulation, culture viability, metabolite production and cell morphology were analysed. During retentostat culturing, the specific growth rate decreased to a minimum of 0.00006 h(-1) , corresponding to a doubling time of 470 days. The energy distribution between growth and maintenance-related processes showed that a state of near-zero growth was reached. Remarkably, a filamentous cell morphology emerged, suggesting that cell separation is impaired under near-zero growth conditions. To evaluate the corresponding molecular adaptations to extremely low specific growth, transcriptome changes were analysed. These revealed that cellular responses to near-zero growth conditions share several similarities with those of cells during the stationary phase of batch growth. However, fundamental differences between these two non-growing states are apparent by their high viability and absence of stationary phase mutagenesis under near-zero growth conditions.

Cross-feeding by Bifidobacterium breve UCC2003 during co-cultivation with Bifidobacterium bifidum PRL2010 in a mucin-based medium

Background

Bifidobacteria constitute a specific group of commensal bacteria that commonly inhabit the mammalian gastrointestinal tract. Bifidobacterium breve UCC2003 was previously shown to utilize a variety of plant/diet/host-derived carbohydrates, including cellodextrin, starch and galactan, as well as the mucin and HMO-derived monosaccharide, sialic acid. In the current study, we investigated the ability of this strain to utilize parts of a host-derived source of carbohydrate, namely the mucin glycoprotein, when grown in co-culture with the mucin-degrading Bifidobacterium bifidum PRL2010.

Results

B. breve UCC2003 was shown to exhibit growth properties in a mucin-based medium, but only when grown in the presence of B. bifidum PRL2010, which is known to metabolize mucin. A combination of HPAEC-PAD and transcriptome analyses identified some of the possible monosaccharides and oligosaccharides which support this enhanced co-cultivation growth/viability phenotype.

Conclusion

This study describes the potential existence of a gut commensal relationship between two bifidobacterial species. We demonstrate the in vitro ability of B. breve UCC2003 to cross-feed on sugars released by the mucin-degrading activity of B. bifidum PRL2010, thus advancing our knowledge on the metabolic adaptability which allows the former strain to colonize the (infant) gut by its extensive metabolic abilities to (co-)utilize available carbohydrate sources.

Molecular characterization of three Lactobacillus delbrueckii subsp. bulgaricus phages

In this study, three phages infecting Lactobacillus delbrueckii subsp. bulgaricus, named Ld3, Ld17, and Ld25A, were isolated from whey samples obtained from various industrial fermentations. These phages were further characterized in a multifaceted approach: (i) biological and physical characterization through host range analysis and electron microscopy; (ii) genetic assessment through genome analysis; (iii) mass spectrometry analysis of the structural components of the phages; and (iv), for one phage, transcriptional analysis by Northern hybridization, reverse transcription-PCR, and primer extension. The three obtained phage genomes display high levels of sequence identity to each other and to genomes of the so-called group b L. delbrueckii phages c5, LL-Ku, and phiLdb, where some of the observed differences are believed to be responsible for host range variations.

Bacillus subtilis attachment to Aspergillus niger hyphae results in mutually altered metabolism

Interaction between microbes affects the growth, metabolism and differentiation of members of the microbial community. While direct and indirect competition, like antagonism and nutrient consumption have a negative effect on the interacting members of the population, microbes have also evolved in nature not only to fight, but in some cases to adapt to or support each other, while increasing the fitness of the community. The presence of bacteria and fungi in soil results in various interactions including mutualism. Bacilli attach to the plant root and form complex communities in the rhizosphere. Bacillus subtilis, when grown in the presence of Aspergillus niger, interacts similarly with the fungus, by attaching and growing on the hyphae. Based on data obtained in a dual transcriptome experiment, we suggest that both fungi and bacteria alter their metabolism during this interaction. Interestingly, the transcription of genes related to the antifungal and putative antibacterial defence mechanism of B. subtilis and A. niger, respectively, are decreased upon attachment of bacteria to the mycelia. Analysis of the culture supernatant suggests that surfactin production by B. subtilis was reduced when the bacterium was co-cultivated with the fungus. Our experiments provide new insights into the interaction between a bacterium and a fungus.

Lactobacillus reuteri 100-23 modulates urea hydrolysis in the murine stomach

Comparisons of in vivo (mouse stomach) and in vitro (laboratory culture) transcriptomes of Lactobacillus reuteri strain 100-23 were made by microarray analysis. These comparisons revealed the upregulation of genes associated with acid tolerance, including urease production, in the mouse stomach. Inactivation of the ureC gene reduced the acid tolerance of strain 100-23 in vitro, and the mutant was outcompeted by the wild type in the gut of ex-Lactobacillus-free mice. Urine analysis showed that stable isotope-labeled urea, administered by gavage, was metabolized to a greater extent in Lactobacillus-free mice than animals colonized by strain 100-23. This surprising observation was associated with higher levels of urease activity and fecal-type bacteria in the stomach digesta of Lactobacillus-free mice. Despite the modulation of urea hydrolysis in the stomach, recycling of urea nitrogen in the murine host was not affected since the essential amino acid isoleucine, labeled with a stable isotope, was detected in the livers of both Lactobacillus-free and 100-23-colonized animals. Therefore, our experiments reveal a new and unexpected impact of Lactobacillus colonization on urea hydrolysis in the murine gut.

Antibiotic-induced replication stress triggers bacterial competence by increasing gene dosage near the origin

Streptococcus pneumoniae (pneumococcus) kills nearly 1 million children annually, and the emergence of antibiotic-resistant strains poses a serious threat to human health. Because pneumococci can take up DNA from their environment by a process called competence, genes associated with antibiotic resistance can rapidly spread. Remarkably, competence is activated in response to several antibiotics. Here, we demonstrate that antibiotics targeting DNA replication cause an increase in the copy number of genes proximal to the origin of replication (oriC). As the genes required for competence initiation are located near oriC, competence is thereby activated. Transcriptome analyses show that antibiotics targeting DNA replication also upregulate origin-proximal gene expression in other bacteria. This mechanism is a direct, intrinsic consequence of replication fork stalling. Our data suggest that evolution has conserved the oriC-proximal location of important genes in bacteria to allow for a robust response to replication stress without the need for complex gene-regulatory pathways. PAPERCLIP:

Autoinducer-2 plays a crucial role in gut colonization and probiotic functionality of Bifidobacterium breve UCC2003

In the present study we show that luxS of Bifidobacterium breve UCC2003 is involved in the production of the interspecies signaling molecule autoinducer-2 (AI-2), and that this gene is essential for gastrointestinal colonization of a murine host, while it is also involved in providing protection against Salmonella infection in Caenorhabditis elegans. We demonstrate that a B. breve luxS-insertion mutant is significantly more susceptible to iron chelators than the WT strain and that this sensitivity can be partially reverted in the presence of the AI-2 precursor DPD. Furthermore, we show that several genes of an iron starvation-induced gene cluster, which are downregulated in the luxS-insertion mutant and which encodes a presumed iron-uptake system, are transcriptionally upregulated under in vivo conditions. Mutation of two genes of this cluster in B. breve UCC2003 renders the derived mutant strains sensitive to iron chelators while deficient in their ability to confer gut pathogen protection to Salmonella-infected nematodes. Since a functional luxS gene is present in all tested members of the genus Bifidobacterium, we conclude that bifidobacteria operate a LuxS-mediated system for gut colonization and pathogen protection that is correlated with iron acquisition.

Metabolism of sialic acid by Bifidobacterium breve UCC2003

Bifidobacteria constitute a specific group of commensal bacteria that inhabit the gastrointestinal tracts of humans and other mammals. Bifidobacterium breve UCC2003 has previously been shown to utilize several plant-derived carbohydrates that include cellodextrins, starch, and galactan. In the present study, we investigated the ability of this strain to utilize the mucin- and human milk oligosaccharide (HMO)-derived carbohydrate sialic acid. Using a combination of transcriptomic and functional genomic approaches, we identified a gene cluster dedicated to the uptake and metabolism of sialic acid. Furthermore, we demonstrate that B. breve UCC2003 can cross feed on sialic acid derived from the metabolism of 3'-sialyllactose, an abundant HMO, by another infant gut bifidobacterial strain, Bifidobacterium bifidum PRL2010.

Bacillus pumilus reveals a remarkably high resistance to hydrogen peroxide provoked oxidative stress

Bacillus pumilus is characterized by a higher oxidative stress resistance than other comparable industrially relevant Bacilli such as B. subtilis or B. licheniformis. In this study the response of B. pumilus to oxidative stress was investigated during a treatment with high concentrations of hydrogen peroxide at the proteome, transcriptome and metabolome level. Genes/proteins belonging to regulons, which are known to have important functions in the oxidative stress response of other organisms, were found to be upregulated, such as the Fur, Spx, SOS or CtsR regulon. Strikingly, parts of the fundamental PerR regulon responding to peroxide stress in B. subtilis are not encoded in the B. pumilus genome. Thus, B. pumilus misses the catalase KatA, the DNA-protection protein MrgA or the alkyl hydroperoxide reductase AhpCF. Data of this study suggests that the catalase KatX2 takes over the function of the missing KatA in the oxidative stress response of B. pumilus. The genome-wide expression analysis revealed an induction of bacillithiol (Cys-GlcN-malate, BSH) relevant genes. An analysis of the intracellular metabolites detected high intracellular levels of this protective metabolite, which indicates the importance of bacillithiol in the peroxide stress resistance of B. pumilus.

The YmdB phosphodiesterase is a global regulator of late adaptive responses in Bacillus subtilis

Bacillus subtilis mutants lacking ymdB are unable to form biofilms, exhibit a strong overexpression of the flagellin gene hag, and are deficient in SlrR, a SinR antagonist. Here, we report the functional and structural characterization of YmdB, and we find that YmdB is a phosphodiesterase with activity against 2',3'- and 3',5'-cyclic nucleotide monophosphates. The structure of YmdB reveals that the enzyme adopts a conserved phosphodiesterase fold with a binuclear metal center. Mutagenesis of a catalytically crucial residue demonstrates that the enzymatic activity of YmdB is essential for biofilm formation. The deletion of ymdB affects the expression of more than 800 genes; the levels of the σ(D)-dependent motility regulon and several sporulation genes are increased, and the levels of the SinR-repressed biofilm genes are decreased, confirming the role of YmdB in regulating late adaptive responses of B. subtilis.

Metabolism of four α-glycosidic linkage-containing oligosaccharides by Bifidobacterium breve UCC2003

Members of the genus Bifidobacterium are common inhabitants of the gastrointestinal tracts of humans and other mammals, where they ferment many diet-derived carbohydrates that cannot be digested by their hosts. To extend our understanding of bifidobacterial carbohydrate utilization, we investigated the molecular mechanisms by which 11 strains of Bifidobacterium breve metabolize four distinct α-glucose- and/or α-galactose-containing oligosaccharides, namely, raffinose, stachyose, melibiose, and melezitose. Here we demonstrate that all B. breve strains examined possess the ability to utilize raffinose, stachyose, and melibiose. However, the ability to metabolize melezitose was not common to all B. breve strains tested. Transcriptomic and functional genomic approaches identified a gene cluster dedicated to the metabolism of α-galactose-containing carbohydrates, while an adjacent gene cluster, dedicated to the metabolism of α-glucose-containing melezitose, was identified in strains that are able to use this carbohydrate.

Assessing the validity and reproducibility of genome-scale predictions

MOTIVATION

Validation and reproducibility of results is a central and pressing issue in genomics. Several recent embarrassing incidents involving the irreproducibility of high-profile studies have illustrated the importance of this issue and the need for rigorous methods for the assessment of reproducibility.

RESULTS

Here, we describe an existing statistical model that is very well suited to this problem. We explain its utility for assessing the reproducibility of validation experiments, and apply it to a genome-scale study of adenosine deaminase acting on RNA (ADAR)-mediated RNA editing in Drosophila. We also introduce a statistical method for planning validation experiments that will obtain the tightest reproducibility confidence limits, which, for a fixed total number of experiments, returns the optimal number of replicates for the study.

AVAILABILITY

Downloadable software and a web service for both the analysis of data from a reproducibility study and for the optimal design of these studies is provided at http://ccmbweb.ccv.brown.edu/reproducibility.html .

Metabolic and transcriptional analysis of acid stress in Lactococcus lactis, with a focus on the kinetics of lactic acid pools

The effect of pH on the glucose metabolism of non-growing cells of L. lactis MG1363 was studied by in vivo NMR in the range 4.8 to 6.5. Immediate pH effects on glucose transporters and/or enzyme activities were distinguished from transcriptional/translational effects by using cells grown at the optimal pH of 6.5 or pre-adjusted to low pH by growth at 5.1. In cells grown at pH 5.1, glucose metabolism proceeds at a rate 35% higher than in non-adjusted cells at the same pH. Besides the upregulation of stress-related genes (such as dnaK and groEL), cells adjusted to low pH overexpressed H(+)-ATPase subunits as well as glycolytic genes. At sub-optimal pHs, the total intracellular pool of lactic acid reached approximately 500 mM in cells grown at optimal pH and about 700 mM in cells grown at pH 5.1. These high levels, together with good pH homeostasis (internal pH always above 6), imply intracellular accumulation of the ionized form of lactic acid (lactate anion), and the concomitant export of the equivalent protons. The average number, n, of protons exported with each lactate anion was determined directly from the kinetics of accumulation of intra- and extracellular lactic acid as monitored online by (13)C-NMR. In cells non-adjusted to low pH, n varies between 2 and 1 during glucose consumption, suggesting an inhibitory effect of intracellular lactate on proton export. We confirmed that extracellular lactate did not affect the lactate: proton stoichiometry. In adjusted cells, n was lower and varied less, indicating a different mix of lactic acid exporters less affected by the high level of intracellular lactate. A qualitative model for pH effects and acid stress adaptation is proposed on the basis of these results.

Lytic infection of Lactococcus lactis by bacteriophages Tuc2009 and c2 triggers alternative transcriptional host responses

Here we present an entire temporal transcriptional profile of Lactococcus lactis subsp. cremoris UC509.9 undergoing lytic infection with two distinct bacteriophages, Tuc2009 and c2. Furthermore, corresponding high-resolution whole-phage genome tiling arrays of both bacteriophages were performed throughout lytic infection. Whole-genome microarrays performed at various time points postinfection demonstrated a rather modest impact on host transcription. The majority of changes in the host transcriptome occur during late infection stages; few changes in host gene transcription occur during the immediate and early infection stages. Alterations in the L. lactis UC509.9 transcriptome during lytic infection appear to be phage specific, with relatively few differentially transcribed genes shared between cells infected with Tuc2009 and those infected with c2. Despite the apparent lack of a coordinated general phage response, three themes common to both infections were noted: alternative transcription of genes involved in catabolic flux and energy production, differential transcription of genes involved in cell wall modification, and differential transcription of genes involved in the conversion of ribonucleotides to deoxyribonucleotides. The transcriptional profiles of both bacteriophages during lytic infection generally correlated with the findings of previous studies and allowed the confirmation of previously predicted promoter sequences. In addition, the host transcriptional response to lysogenization with Tuc2009 was monitored along with tiling array analysis of Tuc2009 in the lysogenic state. Analysis identified 44 host genes with altered transcription during lysogeny, 36 of which displayed levels of transcription significantly reduced from those for uninfected cells.

The transcriptional and gene regulatory network of Lactococcus lactis MG1363 during growth in milk

In the present study we examine the changes in the expression of genes of Lactococcus lactis subspecies cremoris MG1363 during growth in milk. To reveal which specific classes of genes (pathways, operons, regulons, COGs) are important, we performed a transcriptome time series experiment. Global analysis of gene expression over time showed that L. lactis adapted quickly to the environmental changes. Using upstream sequences of genes with correlated gene expression profiles, we uncovered a substantial number of putative DNA binding motifs that may be relevant for L. lactis fermentative growth in milk. All available novel and literature-derived data were integrated into network reconstruction building blocks, which were used to reconstruct and visualize the L. lactis gene regulatory network. This network enables easy mining in the chrono-transcriptomics data. A freely available website at http://milkts.molgenrug.nl gives full access to all transcriptome data, to the reconstructed network and to the individual network building blocks.

Functional analysis of the ComK protein of Bacillus coagulans

The genes for DNA uptake and recombination in Bacilli are commonly regulated by the transcriptional factor ComK. We have identified a ComK homologue in Bacillus coagulans, an industrial relevant organism that is recalcitrant for transformation. Introduction of B. coagulans comK gene under its own promoter region into Bacillus subtilis comK strain results in low transcriptional induction of the late competence gene comGA, but lacking bistable expression. The promoter regions of B. coagulans comK and the comGA genes are recognized in B. subtilis and expression from these promoters is activated by B. subtilis ComK. Purified ComK protein of B. coagulans showed DNA-binding ability in gel retardation assays with B. subtilis- and B. coagulans-derived probes. These experiments suggest that the function of B. coagulans ComK is similar to that of ComK of B. subtilis. When its own comK is overexpressed in B. coagulans the comGA gene expression increases 40-fold, while the expression of another late competence gene, comC is not elevated and no reproducible DNA-uptake could be observed under these conditions. Our results demonstrate that B. coagulans ComK can recognize several B.subtilis comK-responsive elements, and vice versa, but indicate that the activation of the transcription of complete sets of genes coding for a putative DNA uptake apparatus in B. coagulans might differ from that of B. subtilis.

DEAD-Box RNA helicases in Bacillus subtilis have multiple functions and act independently from each other

DEAD-box RNA helicases play important roles in remodeling RNA molecules and in facilitating a variety of RNA-protein interactions that are key to many essential cellular processes. In spite of the importance of RNA, our knowledge about RNA helicases is limited. In this study, we investigated the role of the four DEAD-box RNA helicases in the Gram-positive model organism Bacillus subtilis. A strain deleted of all RNA helicases is able to grow at 37°C but not at lower temperatures. The deletion of cshA, cshB, or yfmL in particular leads to cold-sensitive phenotypes. Moreover, these mutant strains exhibit unique defects in ribosome biogenesis, suggesting distinct functions for the individual enzymes in this process. Based on protein accumulation, severity of the cold-sensitive phenotype, and the interaction with components of the RNA degradosome, CshA is the major RNA helicase of B. subtilis. To unravel the functions of CshA in addition to ribosome biogenesis, we conducted microarray analysis and identified the ysbAB and frlBONMD mRNAs as targets that are strongly affected by the deletion of the cshA gene. Our findings suggest that the different helicases make distinct contributions to the physiology of B. subtilis. Ribosome biogenesis and RNA degradation are two of their major tasks in B. subtilis.

Genetic requirements for Moraxella catarrhalis growth under iron-limiting conditions

Iron sequestration by the human host is a first line defence against respiratory pathogens like Moraxella catarrhalis, which consequently experiences a period of iron starvation during colonization. We determined the genetic requirements for M. catarrhalis BBH18 growth during iron starvation using the high-throughput genome-wide screening technology genomic array footprinting (GAF). By subjecting a large random transposon mutant library to growth under iron-limiting conditions, mutants of the MCR_0996-rhlB-yggW operon, rnd, and MCR_0457 were negatively selected. Growth experiments using directed mutants confirmed the GAF phenotypes with ΔyggW (putative haem-shuttling protein) and ΔMCR_0457 (hypothetical protein) most severely attenuated during iron starvation, phenotypes which were restored upon genetic complementation of the deleted genes. Deletion of yggW resulted in similar attenuated phenotypes in three additional strains. Transcriptional profiles of ΔyggW and ΔMCR_0457 were highly altered with 393 and 192 differentially expressed genes respectively. In all five mutants, expression of nitrate reductase genes was increased and of nitrite reductase decreased, suggesting an impaired aerobic respiration. Alteration of iron metabolism may affect nasopharyngeal colonization as adherence of all mutants to respiratory tract epithelial cells was attenuated. In conclusion, we elucidated the genetic requirements for M. catarrhalis growth during iron starvation and characterized the roles of the identified genes in bacterial growth and host interaction.

Characterization of the ROK-family transcriptional regulator RokA of Streptococcus pneumoniae D39

The Gram-positive human pathogen Streptococcus pneumoniae possesses an unusually high number of gene clusters specific for carbohydrate utilization. This provides it with the ability to use a wide array of sugars, which may aid during infection and survival in different environmental conditions present in the host. In this study, the regulatory mechanism of transcription of a gene cluster, SPD0424-8, putatively encoding a cellobiose/lactose-specific phosphotransferase system is investigated. We demonstrate that this gene cluster is transcribed as one transcriptional unit directed by the promoter of the SPD0424 gene. Upstream of SPD0424, a gene was identified encoding a ROK-family transcriptional regulator (RokA: SPD0423). DNA microarray and transcriptional reporter analyses with a rokA mutant revealed that RokA acts as a transcriptional repressor of the SPD0424-8 operon. Furthermore, we identified a 25 bp AT-rich DNA operator site (5'-TATATTTAATTTATAAAAAATAAAA-3') in the promoter region of SPD0424, which was validated by promoter truncation studies, DNase I footprinting and electrophoretic mobility-shift assays. We tested a large range of different sugars for their effect on the expression of the SPD0424-8 operon, but only moderate variation in expression was observed in the conditions applied. Therefore, a co-factor for RokA-mediated transcriptional control could not be identified.

The putative lactococcal extracytoplasmic function anti-sigma factor llmg2447 determines resistance to the cell wall-active bacteriocin lcn972

Lactococcin 972 (Lcn972) is a cell wall-active bacteriocin that inhibits cell wall biosynthesis in Lactococcus lactis. In this work, the transcriptomes of the Lcn972-resistant (Lcn(r)) mutant L. lactis D1 and its parent strain were compared to identify factors involved in Lcn972 resistance. Upregulated genes included members of the cell envelope stress (CesSR) regulon, the penicillin-binding protein pbpX gene and gene llmg2447, which may encode a putative extracytoplasmic function (ECF) anti-sigma factor. The gene llmg2447 is located downstream of the nonfunctional ECF gene sigX(pseudo). Nisin-controlled expression of llmg2447 led to high Lcn972 resistance in L. lactis, with no cross-resistance to other cell wall-active antimicrobials. Upregulation of llmg2447 in L. lactis D1 (Lcn(r)) was linked to the integration of insertion element IS981 into the llmg2447 promoter region, replacing the native -35 box and activating the otherwise silent promoter P(2447). This is the first example of an orphan ECF anti-sigma factor involved in bacteriocin resistance. This new role in neutralizing cell wall-active compounds (e.g., Lcn972) could have evolved from a putative primary function of Llmg2447 in sensing cell envelope stress.

A two-component regulatory system controls autoregulated serpin expression in Bifidobacterium breve UCC2003

This work reports on the identification and molecular characterization of a two-component regulatory system (2CRS), encoded by serRK, which is believed to control the expression of the ser(2003) locus in Bifidobacterium breve UCC2003. The ser(2003) locus consists of two genes, Bbr_1319 (sagA) and Bbr_1320 (serU), which are predicted to encode a hypothetical membrane-associated protein and a serpin-like protein, respectively. The response regulator SerR was shown to bind to the promoter region of ser(2003), and the probable recognition sequence of SerR was determined by a combinatorial approach of in vitro site-directed mutagenesis coupled to transcriptional fusion and electrophoretic mobility shift assays (EMSAs). The importance of the serRK 2CRS in the response of B. breve to protease-mediated induction was confirmed by generating a B. breve serR insertion mutant, which was shown to exhibit altered ser(2003) transcriptional induction patterns compared to the parent strain, UCC2003. Interestingly, the analysis of a B. breve serU mutant revealed that the SerRK signaling pathway appears to include a SerU-dependent autoregulatory loop.

A conserved two-component signal transduction system controls the response to phosphate starvation in Bifidobacterium breve UCC2003

This work reports on the identification and molecular characterization of the two-component regulatory system (2CRS) PhoRP, which controls the response to inorganic phosphate (P(i)) starvation in Bifidobacterium breve UCC2003. The response regulator PhoP was shown to bind to the promoter region of pstSCAB, specifying a predicted P(i) transporter system, as well as that of phoU, which encodes a putative P(i)-responsive regulatory protein. This interaction is assumed to cause transcriptional modulation under conditions of P(i) limitation. Our data suggest that the phoRP genes are subject to positive autoregulation and, together with pstSCAB and presumably phoU, represent the complete regulon controlled by the phoRP-encoded 2CRS in B. breve UCC2003. Determination of the minimal PhoP binding region combined with bioinformatic analysis revealed the probable recognition sequence of PhoP, designated here as the PHO box, which together with phoRP is conserved among many high-GC-content Gram-positive bacteria. The importance of the phoRP 2CRS in the response of B. breve to P(i) starvation conditions was confirmed by analysis of a B. breve phoP insertion mutant which exhibited decreased growth under phosphate-limiting conditions compared to its parent strain UCC2003.