1
|
Fujita A, Patriota AG, Sato JR, Miyano S. The impact of measurement errors in the identification of regulatory networks. BMC Bioinformatics 2009; 10:412. [PMID: 20003382 PMCID: PMC2811120 DOI: 10.1186/1471-2105-10-412] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/16/2009] [Accepted: 12/13/2009] [Indexed: 11/21/2022] Open
Abstract
BACKGROUND There are several studies in the literature depicting measurement error in gene expression data and also, several others about regulatory network models. However, only a little fraction describes a combination of measurement error in mathematical regulatory networks and shows how to identify these networks under different rates of noise. RESULTS This article investigates the effects of measurement error on the estimation of the parameters in regulatory networks. Simulation studies indicate that, in both time series (dependent) and non-time series (independent) data, the measurement error strongly affects the estimated parameters of the regulatory network models, biasing them as predicted by the theory. Moreover, when testing the parameters of the regulatory network models, p-values computed by ignoring the measurement error are not reliable, since the rate of false positives are not controlled under the null hypothesis. In order to overcome these problems, we present an improved version of the Ordinary Least Square estimator in independent (regression models) and dependent (autoregressive models) data when the variables are subject to noises. Moreover, measurement error estimation procedures for microarrays are also described. Simulation results also show that both corrected methods perform better than the standard ones (i.e., ignoring measurement error). The proposed methodologies are illustrated using microarray data from lung cancer patients and mouse liver time series data. CONCLUSIONS Measurement error dangerously affects the identification of regulatory network models, thus, they must be reduced or taken into account in order to avoid erroneous conclusions. This could be one of the reasons for high biological false positive rates identified in actual regulatory network models.
Collapse
Affiliation(s)
- André Fujita
- Computational Science Research Program, RIKEN, 2-1 Hirosawa, Wako, Saitama, 351-0198, Japan
| | - Alexandre G Patriota
- Institute of Mathematics and Statistics, University of São Paulo, Rua do Matão, 1010 - São Paulo, 05508-090, Brazil
| | - João R Sato
- Center of Mathematics, Computation and Cognition, Universidade Federal do ABC, Rua Santa Adélia, 166 - Santo André, 09210-170, Brazil
| | - Satoru Miyano
- Computational Science Research Program, RIKEN, 2-1 Hirosawa, Wako, Saitama, 351-0198, Japan
- Human Genome Center, Institute of Medical Science, University of Tokyo, 4-6-1 Shirokanedai, Minato-ku, Tokyo, 108-8639, Japan
| |
Collapse
|
2
|
Kowalczyk M, Cocaign-Bousquet M, Loubiere P, Bardowski J. Identification and functional characterisation of cellobiose and lactose transport systems in Lactococcus lactis IL1403. Arch Microbiol 2007; 189:187-96. [PMID: 17909747 DOI: 10.1007/s00203-007-0308-8] [Citation(s) in RCA: 14] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/11/2007] [Revised: 08/02/2007] [Accepted: 09/13/2007] [Indexed: 11/24/2022]
Abstract
Physiological, biochemical and macroarray analyses of Lactococcus lactis IL1403 and its ccpA and bglR single and double mutants engaged in lactose and beta-glucosides catabolism were performed. The kinetic analysis indicated the presence of different transport systems for salicin and cellobiose. The control of salicin catabolism was found to be mediated by the transcriptional regulator BglR and the CcpA protein. The transcriptional analysis by macroarray technology of genes from the PEP:PTS regions showed that several genes, like ybhE, celB, ptcB and ptcA, were expressed at higher levels both in wild type cells exposed to cellobiose and in the ccpA mutant. We also demonstrated that in L. lactis IL1403 cultured on medium with cellobiose and lactose as carbon sources, after the first phase of cellobiose consumption and then co-metabolism of the two sugars, when cellobiose is exhausted the strain uses lactose as the only carbon source. These data could indicate that lactose and cellobiose are transported by a unique system-a PTS carrier induced by the presence of cellobiose, and negatively controlled by the CcpA regulator.
Collapse
Affiliation(s)
- Magdalena Kowalczyk
- Institute of Biochemistry and Biophysics, Polish Academy of Sciences, Pawinskiego 5a, 02-106 Warsaw, Poland
| | | | | | | |
Collapse
|
3
|
Abstract
The increasing number of genomic and post-genomic studies on Gram-positive organisms and especially on lactic acid bacteria brings a lot of information on sugar catabolism in these bacteria. Like for many other bacteria, glucose is the most preferred source of carbon and energy for Lactococcus lactis. Other carbon sources can induce their own utilization in the absence of well-metabolized sugar. These processes engage numbers of genes and undergo complex mechanisms of regulation. In this review, we discuss various biochemical and genetic control mechanisms involved in sugar catabolism, like regulation by repressors, activators, antiterminators or carbon catabolite repression control.
Collapse
Affiliation(s)
- Magdalena Kowalczyk
- Department of Microbial Biochemistry, Institute of Biochemistry and Biophysics PAS, Pawinskiego, Warszawa, Poland.
| | | |
Collapse
|
4
|
Mazzoli R, Pessione E, Giuffrida MG, Fattori P, Barello C, Giunta C, Lindley ND. Degradation of aromatic compounds by Acinetobacter radioresistens S13: growth characteristics on single substrates and mixtures. Arch Microbiol 2007; 188:55-68. [PMID: 17483933 DOI: 10.1007/s00203-007-0223-z] [Citation(s) in RCA: 52] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/04/2006] [Revised: 01/22/2007] [Accepted: 02/07/2007] [Indexed: 11/26/2022]
Abstract
Acinetobacter radioresistens S13 is able to grow on phenol or benzoate as the sole carbon and energy source: both these compounds are catabolized through the beta-ketoadipate pathway. Genes encoding the catabolic enzymes for degradation of aromatic compounds are localized on A. radioresistens S13 chromosome and organized in, at least, two distinct sets, one for benzoate degradation and another for phenol catabolism. In the present study, the growth and biodegradation kinetics for benzoate and phenol, and an easily metabolized substrate (acetate) were established. Benzoate was degraded slower and supports a less rapid and efficient growth than either acetate or phenol. A combined transcript-proteomic analysis of some of the major catabolic genes and their products nonetheless has shown that benzoate induces the expression of both benzoate and phenol catabolic operons. This result was confirmed by the fact that benzoate-acclimatized bacteria were rapidly able to degrade phenol too. Finally, the growth and biodegradation kinetics for different mixtures of acetate, benzoate and phenol were determined. Results indicate that a hierarchy of substrate utilization, benzoate > acetate > phenol, occurred: benzoate was the preferred substrate, despite its lower growth and biodegradation parameters. Hypotheses explaining these unusual metabolic features of A. radioresistens S13 are discussed.
Collapse
Affiliation(s)
- Roberto Mazzoli
- Dipartimento di Biologia Animale e dell'Uomo, Università di Torino, Via Accademia Albertina 13, 10123 Torino, Italy.
| | | | | | | | | | | | | |
Collapse
|
5
|
Redon E, Loubière P, Cocaign-Bousquet M. Role of mRNA stability during genome-wide adaptation of Lactococcus lactis to carbon starvation. J Biol Chem 2005; 280:36380-5. [PMID: 16131490 DOI: 10.1074/jbc.m506006200] [Citation(s) in RCA: 64] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022] Open
Abstract
The stability of mRNA was investigated for the first time at the genomic scale during carbon starvation adaptation of Lactococcus lactis IL1403. In exponential phase, mRNA half-lives were correlated positively to open reading frame length. A polypurine sequence, AGGAG, was identified as a putative 5'-stabilizer and inverted repeated sequences as a 3'-destabilizer. These original findings suggested that multiple pathways of mRNA degradation should coexist: internal cleavage, endonuclease cleavage initiated at the 5'-end, and exonuclease attack at the 3'-end. During carbon starvation adaptation, mRNA stability globally increased, but specific mechanisms allowing a wide range of stabilization factors between genes and differential kinetic evolution were involved. A formal method allowing the quantification of the relative influences of transcription and degradation on the mRNA pool control was developed and applied in L. lactis. Gene expression was mostly controlled by altered transcription prior to carbon source exhaustion, while the influence of mRNA stability increased during the starvation phase. This study highlighted that stability modulation in response to adverse growth conditions can govern gene regulation to the same extent as transcription in bacteria.
Collapse
Affiliation(s)
- Emma Redon
- Laboratoire Biotechnologie Bioprocédés, Unité Mixte de Recherche 5504 CNRS, Institut National des Sciences Appliquées, 135 Avenue de Rangueil, 31077 Toulouse Cedex 4, France
| | | | | |
Collapse
|
6
|
Kok J, Buist G, Zomer AL, van Hijum SA, Kuipers OP. Comparative and functional genomics of lactococci. FEMS Microbiol Rev 2005. [DOI: 10.1016/j.fmrre.2005.04.004] [Citation(s) in RCA: 23] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022] Open
|
7
|
Even S, Lindley ND, Cocaign-Bousquet M. Transcriptional, translational and metabolic regulation of glycolysis in Lactococcus lactis subsp. cremoris MG 1363 grown in continuous acidic cultures. MICROBIOLOGY (READING, ENGLAND) 2003; 149:1935-1944. [PMID: 12855744 DOI: 10.1099/mic.0.26146-0] [Citation(s) in RCA: 69] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
The physiological behaviour of Lactococcus lactis subsp. cremoris MG 1363 was characterized in continuous culture under various acidic conditions (pH 4.7-6.6). Biomass yield was diminished in cultures with low pH and the energy dedicated to maintenance increased due to organic acid inhibition and cytoplasmic acidification. Under such acidic conditions, the specific rate of glucose consumption by the bacterium increased, thereby enhancing energy supply. This acceleration of glycolysis was regulated by both an increase in the concentrations of glycolytic enzymes (hierarchical regulation) and the specific modulation of enzyme activities (metabolic regulation). However, when the inhibitory effect of intracellular pH on enzyme activity was taken into account in the model of regulation, metabolite regulation was shown to be the dominant factor controlling pathway flux. The changes in glycolytic enzyme concentrations were not correlated directly to modifications in transcript concentrations. Analyses of the relative contribution of the phenomena controlling enzyme synthesis indicated that translational regulation had a major influence compared to transcriptional regulation. An increase in the translation efficiency was accompanied by an important decrease of total cellular RNA concentrations, confirming that the translation apparatus of L. lactis was optimized under acid stress conditions.
Collapse
Affiliation(s)
- Sergine Even
- Laboratoire Biotechnologie-Bioprocédés, UMR 5504 INSA/CNRS, UMR 792 INSA/INRA, Institut National des Sciences Appliquées, 135 Avenue de Rangueil, 31077 Toulouse Cedex 4, France
| | - Nic D Lindley
- Laboratoire Biotechnologie-Bioprocédés, UMR 5504 INSA/CNRS, UMR 792 INSA/INRA, Institut National des Sciences Appliquées, 135 Avenue de Rangueil, 31077 Toulouse Cedex 4, France
| | - Muriel Cocaign-Bousquet
- Laboratoire Biotechnologie-Bioprocédés, UMR 5504 INSA/CNRS, UMR 792 INSA/INRA, Institut National des Sciences Appliquées, 135 Avenue de Rangueil, 31077 Toulouse Cedex 4, France
| |
Collapse
|
8
|
Bergès H, Lauber E, Liebe C, Batut J, Kahn D, de Bruijn FJ, Ampe F. Development of Sinorhizobium meliloti pilot macroarrays for transcriptome analysis. Appl Environ Microbiol 2003; 69:1214-9. [PMID: 12571049 PMCID: PMC143623 DOI: 10.1128/aem.69.2.1214-1219.2003] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022] Open
Abstract
In order to prepare for whole-genome expression analysis in Sinorhizobium meliloti, pilot DNA macroarrays were designed for 34 genes of known regulation. The experimental parameters assessed were the length of the PCR products, the influence of a tag at the 5' end of the primers, and the method of RNA labeling. Variance and principal-component analysis showed that the most important nonbiological parameter was the labeling method. The sizes of PCR products were also found to be important, whereas the influence of 5' tags was minimal. The variability between replicated spots on a membrane was found to be low. These experimental procedures were validated by analyzing the effects of microaerobic conditions on gene expression.
Collapse
Affiliation(s)
- Hélène Bergès
- Laboratoire de Biologie Moléculaire des Relations Plantes-Microorganismes, UMR215 CNRS-INRA, 31326 Castanet-Tolosan Cedex, France
| | | | | | | | | | | | | |
Collapse
|
9
|
Even S, Lindley ND, Loubière P, Cocaign-Bousquet M. Dynamic response of catabolic pathways to autoacidification in Lactococcus lactis: transcript profiling and stability in relation to metabolic and energetic constraints. Mol Microbiol 2002; 45:1143-52. [PMID: 12180931 DOI: 10.1046/j.1365-2958.2002.03086.x] [Citation(s) in RCA: 55] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022]
Abstract
The dynamic response of the central metabolic pathways to autoacidification (accumulation of organic acid fermentation products) in Lactococcus lactis was investigated in a global manner by integrating molecular data (cellular transcript concentrations, mRNA turnover) within physiological investigations of metabolic and energetic parameters. The decrease in pH associated with the accumulation of organic acids modified the physiological state of the cell considerably. Cytoplasmic acidification led to inhibition of enzyme activities and, consequently, to a diminished catabolic flux through glycolysis and a decreased rate of biochemical energy synthesis. This decrease in energy production together with the increased energy expenditure to counter cytoplasmic acidification led to energetic limitations for biomass synthesis. In these conditions, the specific growth rate decreased progressively, and growth ultimately stopped, although a diminished catabolic flux was maintained in the absence of growth. The cellular response to this phenomenon was to maintain significant levels of mRNA of catabolic genes, involving both continued transcription of the genes and also, in certain cases, an increase in transcript stability. Thus, translation was maintained, and intracellular concentration of certain enzymes increased, partially compensating for the inhibition of activity provoked by the diminished pH. When catabolic activity ceased after prolonged exposure to stress-induced stationary phase, endogenous RNA catabolism was observed.
Collapse
Affiliation(s)
- Sergine Even
- Laboratorie Biotechnologie-Bioprocédés, UMR 5504 INSA/CNRS and UMR 792 INSA/INRA, Centre de Bioingénierie Gilbert Durand, Institute National des Sciences Appliquées, Toulouse, France
| | | | | | | |
Collapse
|