Abstract
More than half of the protein-coding genes in bacteria are organized in polycistronic operons composed of two or more genes. It remains under debate whether the operon organization maintains the stoichiometric expression of the genes within an operon. In this study, we performed a label-free data-independent acquisition hyper reaction monitoring mass-spectrometry (HRM-MS) experiment to quantify the Escherichia coli proteome in exponential phase and quantified 93.6% of the cytosolic proteins, covering 67.9% and 56.0% of the translating polycistronic operons in BW25113 and MG1655 strains, respectively. We found that the translational regulation contributes largely to the proteome complexity: the shorter operons tend to be more tightly controlled for stoichiometry than longer operons; the operons which mainly code for complexes is more tightly controlled for stoichiometry than the operons which mainly code for metabolic pathways. The gene interval (distance between adjacent genes in one operon) may serve as a regulatory factor for stoichiometry. The catalytic efficiency might be a driving force for differential expression of enzymes encoded in one operon. These results illustrated the multifaceted nature of the operon regulation: the operon unified transcriptional level and gene-specific translational level. This multi-level regulation benefits the host by optimizing the efficiency of the productivity of metabolic pathways and maintenance of different types of protein complexes.
Introduction
An operon is a cluster of genes transcribed in a single mRNA. This principle is conserved across bacterial and archaeal genomes, as well as mitochondria and chloroplast (; ; ). Operons are also found in virus and some lower eukaryotes, including yeasts, nematodes, and insects (; ; ; ; ). In a typical bacterial genome, more than half of the protein-coding genes are organized in multigene operons. A classical bacterial operon generates an mRNA strand with polycistronic structure containing multiple coding sequences and are translated together in the cytoplasm. These genes are often of related functions, for example, to build a protein complex or to participate in one metabolic pathway. Therefore, grouping related genes as operons under the control of a single promotor is often thought to simplify the regulation of gene expression for rapid adaptation to environmental changes.
An intuitive presumption of the operon organization is to maintain stoichiometry of the gene products. It was argued that co-regulation could be evolved by merging two independent genes in proximity together under the control of the same promoter, to reduce the control complexity (; ), measured protein synthesis rates by using ribosome profiling and implied that the synthesis rates quantitatively might reflect the stoichiometry of the protein complexes. Studies showed that an operon with one complex promoter might be better than two independent promoters; organization of genes in operons substantially reduces the shortfall in production of complex-forming individual proteins (; ). However, recent advances of omics techniques raised counter-arguments. A transcriptome-level study revealed that certain adjacent genes within one operon are not similarly transcribed in M. pneumoniae. In half of the polycistronic operons, genes exhibited a decaying expression according to its rank in the operon, which is termed “staircase-like decay behavior” (). Considering the widespread post-transcriptional regulations including translational control and protein turnover (), it is still under intensive debate whether this “staircase-behavior” influences the protein abundance (; ; ).
Theoretically, proteins in a complex should follow the stoichiometry, while the proteins involved in the same pathway may need differential expression controls (). For example, the enzymes in various amino acids synthesis pathways are regulated in single-input modules (SIMs). A series of such enzymes are successively expressed in one operon (; ). Meanwhile, the different catalytic kinetics of these enzymes determines that these enzymes should not be expressed at the same quantity (). These genes tend to duplicate to evolve a larger gene regulatory network (), indicating their regulation is less stringent, and an operon arrangement might be unnecessary. Therefore, a more detailed proteome-wide and quantitative investigation is necessary to discover the scope and impact of the operons in gene expression regulation.
A method capable to assess a quantification of the proteome should be used in this case. Although stable isotope labeling methods are more accurate than label-free mass-spectrometry (MS) methods (), the isotopes may affect the physiology of the bacteria (). The isotope labeling is more suitable for comparative quantification of multiple samples than estimating abundance of the proteins within one sample (). Therefore, label-free MS methods should be used. compared three label-free quantification methods (iBAQ, emPAI, and APEX) and found a staircase-like protein expression in most of the transcription units, and found high correlation abundances between some well-known complex subunits. In contrast, found only 5% “staircase behavior” for L. interrogans operons on the proteome level. These contradictory results reflected the cons of these label-free MS approaches: the technical variations and relatively low number of quantified proteins restricted the accurate and in-depth coverage of operon-controlled genes.
In this work, we set out to tackle these problems by employing a highly accurate label-free method, DIA (data-independent acquisition) (), to obtain quantification of the proteins constituting the Escherichia coli proteome with a high coverage and high accuracy. DIA is a MS-based proteomics method used in peptide quantification, in which all ions within a selected m/z range are fragmented and analyzed in a second stage of tandem mass spectrometry (). Although not suitable for discovery-based applications, DIA provides accurate peptide quantification without being limited to profiling predefined peptides of interest (; ). This allowed us to investigate the protein abundances within operons and thus to interrogate the possible stoichiometry in operons of different functions.
Materials and Methods
MS Sample Preparation
Escherichia coli K-12 sub-strains BW25113 and MG1655 were cultivated on glucose M9 minimal medium at 37°C in flasks to mid-exponential phase (OD600 = 0.6) and then harvested in 45 mL volume, immediately cooled in ice water, and then centrifuged at 10,000 × g for 5 min. The pellet was washed once with PBS, centrifuged at 10,000 × g for 5 min again. Pellet was re-suspended on ice with lysis buffer (5 M urea/2 M thiourea in 10 mM HEPES, pH 8.0), and were sonicated and centrifuged at 17,000 × g for 30 min in a table-top centrifuge to remove cell debris. Supernatant was collected, and protein concentrations were determined with a Bradford Protein Assay (Bio-Rad Protein Assay Dye Reagent Concentrate, Cat. #500-0006).
For proteome analysis, we employed in-solution protein digestion with a filter-aided sample preparation (FASP) method (). 1 mg of protein was subjected to reduction (8 M urea and 50 mM DTT at 37°C, 1 h), followed by alkylation with 100 mM iodoacetamide (IAA) in dark at room temperature for 30 min. The solution was transferred to the 30 kDa ultracentrifuge filters (Millipore). Proteins were washed with 8 M urea, and four sequential buffer changes were performed using 50 mM TEAB, respectively. Trypsin (Promega) was then added into the filter at a mass ratio of 1:20 for Proteins digested in 130 μL 50 mM TEAB at 37°C for 12 h. The released peptides were collected by centrifugation and dried with a cold-trap speed vacuum.
MS Experiments
One microgram of sample abovementioned peptides was analyzed on a C18 column (50 μm × 15 cm, 2 μm, Thermo Fisher) by using an EASY-nLC 1200 UHPLC connected to an Orbitrap Fusion Lumos mass spectrometer (Thermo Scientific). The peptides with the iRT-standard (1/10 by volume, Biognosys, HRM Calibration Kit: Ki-3003) were separated by a linear gradient from 6 to 30% ACN with 0.1% formic acid at 270 nL/min for 100–130 min and linearly increased to 90% ACN in 20 min. For the data-dependent acquisition (DDA), the source was operated at 2.0 kV. The DDA scheme included a full MS survey scan from m/z 400 to m/z 1500 at a resolution of 60,000 FWHM with AGC set to 4E5 (maximum injection time of 50 ms), followed by MS/MS scans at a resolution of 15,000 FWHM with AGC set to 5E4 (maximum injection time of 30 ms), data-dependent mode was set to top speed. Isolation window was 1.6. Dynamic exclusion was set to 90 s with a 10 ppm tolerance around the selected precursor. For the DIA hyper reaction monitoring (HRM-MS), individual tryptic peptide samples were mixed with the iRT-standard (1/10 by volume) and analyzed by the same method as DDA used. The method consisted of a full MS1 scan at a resolution of 60,000 FWHM from m/z 350 to m/z 1,200 with AGC set to 4E5 (maximum injection time of 30 ms) followed by 40 non-overlapping DIA windows acquired at a resolution of 30,000 FWHM with AGC set to 5E5 (maximum injection time of 50 ms), cycle time, 3.28 s. The MS/MS isolation windows were listed in Supplementary Table S1. For comparison, standard DDA MS experiment was performed as above. All MS raw data have been deposited in iProX with accession number: IPX0001095000 and ProteomeXchange with identifier PXD010126.
Spectral Library Generation
To generate the spectral library, three DDA measurements of the mixed samples were performed. Raw DDA datasets were searched against a combined database of the NCBI database of Escherichia coli str. K-12 (GCF_000005845.2_ASM584v2, 4140 entries) and the iRT standard peptides sequence using the Sequest HT (Proteome Discoverer v2.1) local server. Common contaminants in the database included trypsin and keratins. Precursor and product ion spectra were searched at an initial mass tolerance of 10 ppm and fragment mass tolerance 0.02 Da, respectively. Tryptic cleavage was selected, and up to two missed cleavages were allowed. Carbamidomethylation on cysteine (+57.021 Da) was set as a fixed modification, and oxidation (+15.995 Da) on methionine was assigned as a variable modification. A target-decoy-based strategy was applied to control peptide and protein false discovery rates (FDRs) at lower than 1%. Confident protein identifications should suit the following criteria: (1) protein level FDR ≤ 1%; (2) unique peptides ≥ 1 or 2; (3) peptide length ≥ 6 or 7 aa. The search result was exported in a pdResult file format containing the annotation of precursors and fragment ions and their exact retention times. The pdResult file was then imported into Spectronaut Pulsar 11 (Biognosys) to generate the spectral library used for HRM-MS data analysis, which yielded 14608 unique peptide sequences in 2041 protein groups with BW25113, and 8822 unique peptide sequences in 1607 protein groups with MG1655. A subset of identified peptides was used in library creation as modification parameter was set none. The generated spectral libraries were exported from Spectronaut as in Supplementary Table S2.
Protein Identification and Quantification
The DIA data were then analyzed with Spectronaut Pulsar 11 with the spectral library, which is a mass spectrometer vendor independent software for SWATH/DIA data analysis. Raw data were analyzed according to the user manual of the software. Default settings were setup for protein identification and peak area calculation. Raw data were converted into HTRMS files and imported to Spectronaut Pulsar 11 by choosing the matched database fasta file and spectral library, with the default settings of the Spectronaut Pulsar 11: (1) Calibration: calibration mode, automatic; iRT calibration strategy, non-liner iRT calibration. (2) Identification: decoy limit strategy, dynamic; decoy method, mutated; machine learning, per run; precursor q-value (peptide FDR) cutoff, 0.01; protein q-value (protein FDR) cutoff, 0.01; p-value estimator, kernel density estimator. (3) Workflow: default labeling type, label; profiling strategy, none; unify peptide peaks, false. (4) Quantification: interference correction, true; major(protein) grouping, by protein-group id; major group quantity, mean peptide quantity; minor (peptide) grouping, by stripped sequence; minor group quantity, mean precursor quantity; minor group top n, true; min, 1; max, 3; quantity MS-level, MS2; quantity type, area; data filtering, q-value; cross run normalization, true; row selection, q-value sparse; normalization strategy, local normalization. (5) Reporting: scoring histograms, true; pipeline report schema, protein quant; pipeline reporting unit, experiment. (6) Protein inference: protein inference workflow, automatic. (7) Data extraction: MS1 mass tolerance strategy, dynamic; correction factor, 1; MS2 mass tolerance strategy, dynamic; correction faction, 1. (8) Post analysis: differential abundance grouping, major group (quantification settings); smallest quantitative unit, precursor ion (summed fragment ions); use top n selection, false. (9) Retention time were used to assist identification. XIC extraction: XIC RT extraction window, dynamic; correction factor, 1. After peak extraction and area calculation were performed, the result was exported as the table format for further quantification analysis in Microsoft Excel. All MS raw data, Proteome Discoverer report (∗.msf file) and the constructed spectra library have been deposited in iProX with accession number: IPX0001095000 and ProteinXchange with identifier PXD010126. The relationship of submitted raw data are shown in Supplementary Data.
Protein abundances was calculated by using Spectronaut Pulsar protein pivot report, those proteins quantified by Spectronaut pulsar but not identified by Proteome Discoverer and not met the confident protein identifications criteria were removed. The abundances of identified proteins were calculated as follow procedure. Supposed that top 500 abundant proteins of E. coli can represent total protein copy numbers of a cell. Concentration of HRM-MS protein copies per cell was calculated based on the means of 500 most abundant protein quantities computed by other three label-free methods (APEX, iBAQ, PAI) downloaded from . HRM-MS intensity could be converted to protein copies per cells by coefficient k, which is defined by the following formula.
Bi is copy numbers of the gene in iBAQ dataset.
Pi is copy numbers of the gene in emPAI dataset.
Ai is copy numbers of the gene in APEX dataset.
Di is correspondence protein intensities of the gene quantified in HRM-MS dataset.
The amount of individual proteins was calculated as the product of conversion coefficient k to their intensity in the HRM-MS measured sample.
The calculations were performed by in-house generated python scripts. All scripts used in this study can be downloaded in the Supplementary Materials (Supplementary Scripts).
Coefficient of Variation of Protein Abundance in the Operons
Coefficient of variation (CV), which is defined as the ratio of the standard deviation to the arithmetic mean. Standard deviation is normalized by n − 1 by default (n is sample size). The CV of proteins within one operon is defined as the ratio of the standard deviation of protein quantities within this operon to the arithmetic mean of all protein quantities within this operon. For multi-gene operons (protein numbers ≥ 2), CV was calculate as follows:
where xi is the abundance of the i-th gene in this operon. To be noted, the CV calculation was only performed within one operon, not across the operons.
The CV of the protein half-life in the operons were calculated in the same way. Protein half-life time in the M9 minimal medium was from our previous work ().
Data Randomization
To compare with the real operon CV level if protein abundances in operons have stoichiometry control, we reshuffled the protein quantities detected in the polycistronic operons (846) randomly to each protein ID, the generated dataset was used as randomized negative control. Randomized protein quantities in “2-/3-/4-/≥5-protein” operons were extracted from this randomized negative control data.
Operon and Gene Ontology (GO) Analysis
Operon library of Escherichia coli str. K-12 were downloaded from the DOOR2 database () (NC_000913). Protein GI numbers were converted to proper identifiers by DAVID Gene Accession Conversion Tool (). The quantified proteins were integrated to the operon data. The PANTHER Version 13.0 (released 2017-11-121) () was used to perform the GO overrepresentation analysis with the significance threshold of 0.01, the quantified proteins of Escherichia coli in our work was selected as the background proteome, the Fisher’s Exact test was used to obtain p-values and ‘GO slim’ category were used. The protein subcellular localization data of E. coli was downloaded from EcoProDB ().
Complex and Pathway Classification
The operon contains more than or equal to two genes were called polycistronic operons. Among polycistronic operons, those ≥90% genes in operon encoded subunits of one protein complex is selected and classified as “Complex” group, others were classified as “Pathway” group.
Physical and Chemical Features of Proteins
The protein lengths in amino acids were obtained from the NCBI of Escherichia coli str. K-12 (GCF_000005845.2_ASM584v2, 4140 entries). Information of the hydrophobicity was calculated by Gravy Calculator2. In addition, the isoelectric point, protein length, instability and hydrophobicity distribution were calculated by using python 2.7 scripts and Biopython libraries.
Experimental Design and Statistical Rationale
To increases the precision of protein expression measurements of the entire E. coli proteome quantification, two biological replicates of BW25113 and MG1655 each were cultured in M9 minimal medium to mid-exponential phase and were harvest, then processed to HRM-MS analysis independently. iRT-standard (Biognosys, HRM Calibration Kit) was added to the peptides with 1/10 by volume. The peptide mixture of two biological replicates of each strain were used and performed LC-MS for three times for spectral library creation. Proteome Discoverer 2.1 and Spectronaut Pulsar 11 were used to generated spectral library, and Spectronaut Pulsar 11 was used to quantify the protein groups with little modified parameters. Kolmogorov–Smirnov test (KS-test) were used to compare the distributions between CVs at transcriptome, translatome and protein level, and Mann-Whitney U-test were used to compare the difference between complex and pathway operons.
mRNA Sequencing
E. coli strain BW25113 was cultivated on glucose M9 minimal medium at 37°C in flasks to mid-exponential phase (OD600 = 0.6) with 100 μg/mL chloramphenicol added 15 min before harvest, then the cells were centrifuged at 5,000 × g for 10 min at 4°C, followed by thrice washed with pre-chilled PBS. Cell pellet was then re-suspended in 6 mL pre-chilled sucrose-buffer solution [16 mM Tris (pH 8.1) supplemented with 0.5 M RNase-free sucrose, 50 mM KCl, 8.75 mM EDTA, 100 μg/mL chloramphenicol, 12.5 mg/mL lysozyme] and gently stirred for 5 min on ice. Then the cells were centrifuged 5,000 × g, 10 min. Total RNA was extracted by Trizol method, and mRNA-seq libraries were prepared using standard MGIEasyTM mRNA Library Prep Kit V2 following the manufacturer’s protocol. Sequencing was performed on a BGISEQ-500 sequencer for 50 cycles, single-ended mode. This dataset was deposited in the GEO database under the accession number GSM3489376, GSM3489377.
Analysis of Sequencing Data
The RNA-seq dataset of E. coli strain MG1655 were obtained from (GEO accession number: GSM1360030, GSM1360031, GSM1360042, GSM1360043) and (SRA accession number: SRR2016457). The datasets of strain BW25113 were generated as described above. For all datasets, adapters was trimmed from the reads. Reads were mapped to coding sequence of E. coli reference genome (GenBank: U00096) using FANSe3 algorithm () with the parameters -E3 -S10 –indel. Genes with at least 10 mapped reads were considered quantifiable (). The expression levels were estimated in rpkM.
Results
Near-Complete E. coli Cytosolic Proteome Quantification Using HRM-MS
To assess the quantification power and reproducibility of the HRM-MS method, we performed two biological replicates of E. coli total soluble proteins of strain BW25113 and MG1655. When using the previous identification criteria (at least one unique peptide, peptide length ≥6 amino acids) (), our HRM-MS results quantified 1951 and 1571 proteins in these two strains, respectively. The two biological replications identified almost identical proteins, demonstrating high robustness and reproducibility (Figure 1A). Two replicates quantified almost the same proteins: only a few proteins were quantified only in one replicate (Figure 1A). Under the stringent criteria as two unique peptides and at least seven amino acids peptide length. Even under the stringent criteria, we still quantified 1675 and 1252 proteins with a high reproducibility (Figure 1B,C and Supplementary Table S3). The number of quantified soluble proteins was almost doubled when compared to the previous results quantified by other methods (1021 proteins for APEX, 1183 for IBAQ and 1138 for emPAI) (). To rule out the difference of the instruments, we performed DDA MS experiments for the two strains in the same Orbitrap Fusion Lumos instrument. Proteins were identified under the stringent criteria and quantified using iBAQ method. The DDA MS quantified 1520 and 1482 proteins for BW25113 and MG1655 strains, respectively, comparable with the HRM-MS experiments. However, the correlation coefficients of the iBAQ quantification of two biological replicates were 0.932 and 0.939, respectively (Figure 1D,E), lower than the HRM-MS (R = 0.982 and 0.988 for the two strains, respectively). Each identified protein was covered by 19.20 and 15.45 peptides in average in two strains, covering 31.56% and 24.87% of the amino acid sequences, respectively (Figure 1F), which is higher than the typical peptide coverage of human proteome MS experiments (single search engine, up to ∼20% coverage) ().
FIGURE 1
We previously revealed 2922 genes which are being translated in the E. coli grown in the same condition using ribosome profiling (). Using HRM-MS method with stringent criteria, we quantified 55.2% of these translating genes in this work. We next analyzed the possible chemical and physical features of the translated but unquantified proteins in this work. The unquantified proteins are significantly more alkalic, less stable, shorter and more hydrophobic (Supplementary Figure S1). These are general factors that decreases the visibility of these proteins in shotgun MS experiments. Since our experiments were not optimized for membrane proteins, which are more prone to aggregate during the protein extraction, these proteins are expected to be less detected in the MS. Notably, we quantified 93.6% translating cytosolic proteins, showing a near-complete quantification of the soluble proteins.
Considering the advantage of the HRM-MS, all the subsequent analysis was based on the HRM-MS under the stringent criteria.
Operons Tend to Unify Gene Expression in General
The high coverage of proteome quantification allowed us to make an in-depth investigation of protein abundances in operons. Indeed, our quantification covered 67.9% and 56.0% of the translating polycistronic operons in BW25113 and MG1655, respectively. We calculated the coefficient of variation (CV, %) of proteins abundances within each operon. To compare with the real operon CV level, we random redistribution the protein quantities quantified in the experiment as randomized negative control. Smaller CV represent the unified expression level of the proteins within one operon. The mean of real CV was significantly smaller than that of randomized negative control data (Mann–Whitney U-test, p = 5.05 × 10−10 and p = 3.40 × 10−5 for BW25113 and MG1655, respectively) (Figure 2). The median CV of quantified proteins were also smaller than the randomized control (Kruskal–Wallis H-test, p = 1.44 × 10−8 and p = 1.18 × 10−6 for BW25113 and MG1655, respectively). As a positive control, the trend to unified expression is also valid in transcription (an operon is transcribed as an entire mRNA) (Supplementary Figure S2). This indicated that most polycistronic operons were co-regulated, and the stoichiometry balance of protein abundances within polycistronic operons exists in general, although with exceptions. These results were in general consistent with previous studies (; ).
FIGURE 2
Functional Enrichment of the Operon Stoichiometry Control
We then examined if the number of genes per operon could affect the stoichiometry of proteins encoded within operons in both BW25113 and MG1655 strains. We found all the medians of subgroups were lower than the randomized negative control (RD), which was in accordance with our abovementioned results (Figure 3A), indicating operon expression regulation exists in general. We continue to divide operons into two subgroups by their functions.
FIGURE 3
Next, each group was separated into two subgroups by the median CV of randomized data, the high CV subgroups (H2-H5, higher than median of randomized data) and the low CV subgroups (L2–L5, lower than median of randomized data). Gene ontology (GO) overrepresentation analysis was performed for each subgroup both in both strains (Figure 3B,C, see details in Supplementary Table S5). Most low CV subgroups (L2–L5) showed functional enrichment against the quantified proteins as background. The L2 group was overrepresented in almost all metabolic activities, and the L5 group was highly enriched in complexes and structural molecules in both E. coli strains. This provided a hint that the stringent stoichiometry control might be important for the efficient assembly of protein complexes. In contrast, there are little significant functional enrichment of GO terms enriched in H2–H5 subgroups in BW25113, and only the GO terms enriched in H2 and H5 subgroup in MG1655. These results indicated that the large CV of most of these operons might be caused by experimental error.
Since the operons in L2–L5 subgroups were enriched in metabolic pathways and complexes, we specifically divide these operons into “Complex” and “Pathway” groups. The group “Complex” are the operons encoding proteins for the same protein complex. The group “Pathway” are the operons encoding proteins involving in the same metabolic pathway. Similar to the Figure 3A, we performed randomization for each group multiple times for robustness. In almost all cases the 2-/3-/4-/≥5-protein “Complex” and “Pathway” operons exhibited lower CVs than the corresponding randomized data (Supplementary Figure S3), indicating that the stoichiometry is still maintained for the protein complexes and metabolic pathways in a certain extent, which is consistent with the traditional hypothesis.
However, there seems to be a trend that shorter operons (containing lower number of genes) possess lower CV within operon (Figure 3A), suggesting a length-dependent stoichiometry control. Since the mRNA of an entire operon is transcribed as a unit, the abovementioned phenomenon opened a question of the origin of the length-dependent stoichiometry control.
Length-Dependence of Stoichiometry Control Is Functional-Dependent
Linear regression analysis was performed to calculate the CVs at RNA and protein levels, within operon that encodes proteins forming complexes or involving in the metabolic pathways (Figure 4A and Supplementary Figure S4). Since the mRNA of one operon is transcribed as one unit, the CV within operon at RNA level is much lower than the randomized control, as expected (Figure 4A and Supplementary Figure S4). However, the “Pathway” operons did not show length dependence at protein level, while significant and positive correlation of protein CV within operon versus length were observed in “Complex” groups in both BW25113 and MG1655 (P-value of BW25113 = 0.0004, P-value of MG1655 = 0.022). In contrast, among the “Pathway” groups, those operons exhibited similar CV distribution regardless of their lengths (regression P > 0.05) (Figure 4A, “Protein” plots). This reflects necessity that the proteins operating in a pathway need to be more independently tuned and thus do not have to follow the stoichiometry. Those results indicating a length- and functional-dependence of operon stoichiometry control for larger operons. The RNA–protein correlation also echoed this trend (Figure 4B). In both strains, the Pearson R2 of the RNA–protein correlation in “Complex” group is considerably higher than in “Pathway” group (Figure 4C). These results indicated that there may be some inherent difference between “Pathway” and “Complex” type operons.
FIGURE 4
Significant difference (p < 0.01) of protein CV distribution were observed among the 2-protein operons against larger (4-/5-protein) operons only in “Complex” subgroups, but not observed in the “Pathway” operon subgroups in both BW25113 and MG1655 strains (Figure 4D, see details in Supplementary Table S4), consolidated our abovementioned observation that the shorter operons among “Complex” groups operons tend to be regulated more stringently.
Enzyme Activity Correlates to the Differential Translation of “Pathway” Genes
Figure 5A showed examples of “Complex” and “Pathway” operons, respectively. Genes in three operons showed almost same RNA abundance. However, the operon ID 3641 encoding ribosome proteins showed similar protein abundance, while the operon ID 3767 encoding enzymes in arginine synthesis pathway and the operon ID 3157 encoding enzymes in biotin synthesis pathway showed exaggerated difference in protein abundance. Both translation and degradation may affect the protein abundance. We found no significant difference of the “Complex” and “Pathway” operons in terms of the CV of protein half-life within operons (KS-test, P = 0.932) (Supplementary Figure S5). Therefore, the translational regulation should be the major factor leading to such differential protein expression in “Pathway” operons.
FIGURE 5
It is expected that the “Complex” proteins tend to be tightly stoichiometrically controlled to build functional complexes. In contrast, what benefit is related to the differential translation of the “Pathway” genes in one operon?
Since “Pathway” proteins comprise pathways, they sequentially catalyze conversion of a substrate to product via multiple and successive reaction steps. Therefore, we hypothesize that the bacteria require less high-efficiency enzymes to avoid energy waste. Supplementary Figure S6 showed the arginine synthesis pathway in E. coli, where ArgB and ArgC are two enzymes that relay. The Michalis constant (KM) of ArgB is more than 3 times higher than ArgC (Table 1), showing that the binding of ArgB to the substrate is weaker. Although lacking the measured value of ArgB kcat in E. coli, ArgB kcat value in yeast (56% homology in amino acid sequence) is approximately 1/3 of ArgC in E. coli, as a reference. These data indicated that the catalytic efficiency of ArgC is higher than ArgB. Therefore, ArgC is less needed in E. coli, which matched the proteome quantification (Figure 5A). No argH activity data is available in E. coli. The homologous enzyme in Anas platyrhynchos (75% homology in amino acid sequence) showed also lower catalytic efficiency than ArgC. Similar trend was also observed in the biotin synthesis pathway (Supplementary Figure S6). The enzyme BioF possess kcat value one order of magnitude higher than BioB, indicating higher conversion efficiency. Therefore, BioF is much less produced in E. coli. The enzyme BioC is extremely efficient with three orders of magnitude higher kcat value than BioF and much lower KM value than bioB and BioF. Therefore, BioC is orders of magnitude lower than the other two enzymes in E. coli.
Table 1
| Enzyme | kcat (s−1) | KM (mM) | Homology to E. coli protein (for non-E. coli proteins) | References |
|---|---|---|---|---|
| ArgB | 4.9 (S. cerevisiae) | 1.3 | 56% | ; |
| ArgC | 14 | 0.4 | ||
| ArgH | 4.9 (Anas platyrhynchos) | 0.4 (Anas platyrhynchos) | 75% | |
| BioB | 0.0039 | 0.002 | ; | |
| BioF | 0.05825 | 0.025 | ||
| BioC | 98.334 (Bacillus cereus) | 0.00108 (Bacillus cereus) | 62% | |
| GuaB | 13 | 0.061 | ||
| GuaA | 23 | 0.053 | ||
| AceB | 0.022 | |||
| AceA | 0.063 | |||
| BetA | 1.5 | |||
| BetB | 1.8 | |||
| UxuA | 4.79 | |||
| UxuB | 1 | |||
| DadA | 30 | |||
| DadX | 3.03 | |||
| ProB | 1.2 (Pseudomonas aeruginosa) | 74% | ||
| ProA | 300 | |||
| OtsA | 1 (Thermoplasma acidophilum) | 65% | ||
| OtsB | 0.61 | |||
| CysC | 0.0005 | |||
| CysD + CysN | 0.0045 (Thiobacillus denitrificans) | 67% (CysD) | ||
| PurD | 0.03 | |||
| PurH | 0.082 (Methanocaldococcus jannaschii) | 50% |
Kinetic parameters of the enzymes in one metabolic pathway.
Values are measured in E. coli unless specified. Non-E. coli data were used as a reference when no E. coli data is available.
To further validate this hypothesis, we manually went through all quantified “Pathway” operons and found other nine operons which fit the following criteria: (a) at least two quantified proteins in one operons, and they must exist in the same metabolic pathway; (b) their pathway must be branch-free, i.e., the substrate should be converted sequentially by these enzymes without introducing other rate-limiting metabolites as branches in pathway, at least in the range of the quantified enzymes (illustrations see Figure 5B); (c) enzyme activity parameters, e.g., kcat or KM, should be available, and the activity parameters of at least one enzyme should be available in E. coli. The RNA abundance, protein abundance and the KM values are illustrated in Figure 5C. All of the nine operons validated the correlation of the enzyme activity and the protein abundance without exception: less active enzymes (represented by higher KM values) were expressed in higher amount at protein level. To be noted, six out of nine operons showed inverse proportion of RNA and protein expression, suggesting that such regulation is mainly conducted at translation level.
Divergence of Stoichiometry Control Is Regulated at Translation Level by the Gene Intervals
Next, we investigated the factor that could determine the lower stringency of “Pathway” groups of operons. As both groups of operons contain a wide variety of genes, the major difference should lay on the gene structures. This is reflected by the gene intervals, defined as the distance from the stop codon of the first gene to the start codon of the next gene downstream within one operon (Figure 6A). A strongly significant difference on gene intervals was observed between two groups on their gene interval distributions (p = 1.314 × 10−5, KS-test, Figure 6B). The genes in “Complex” genes tend to be arranged very near to each other, while the “Pathway” genes tend to be located far from each other, reflected by the larger mean and median value of the gene interval.
FIGURE 6
Both a 70S ribosome and a 30S subunit cover about 40 nucleotides of the mRNA, roughly 18–20 nucleotides upstream and 16–18 nucleotides downstream of the P-site codon (). If a terminating 70S ribosome would have a downstream initiation site nearer than 40 nucleotides, a 70S termination event and a downstream 30S-binding initiation could not be independent events due to a steric clash between 70S and 30S. In this case, the 70S scanning initiation mode plays a major role for translation, meaning the first cistron is initiated by 30S and the downstream cistrons by 70S scanning to achieve a strict and precise 1:1 stoichiometric ratio (). This may be one of the factors explained that the “Complex” operons, with shorter gene intervals, are more stringently regulated for stoichiometry.
Discussion
Near-complete coverage of proteome identification and quantification is always a goal of proteomics research, as it reveals detailed global dynamics of important biological processes; for example, the dormancy and resuscitation of Mycobacterium tuberculosis (). The high proteomic resolution of quantification allows in-depth investigation of many scientific questions in debate for decades. In this study, we employed DIA based HRM-MS, a highly reproducible label-free quantification method in E. coli strains BW25113 and MG1655. Our datasets were better than the previously reported DDA-based E. coli proteome quantification datasets in terms of reproducibility. We therefore generated the most complete investigation on abundance of proteins encoded within operons in E. coli based on tryptic digestion up to date, which allows us to accurately evaluate the stoichiometry presumption of the operon organization. To further increase the sensitivity of identification, other complement approaches such as LysC digestion might be helpful (). Nevertheless, since we have already quantified 93.6% of the cytosolic proteins, specific methods dealing with the hydrophobic nature of membrane proteins should be employed to further expand the proteome coverage. The DIA based HRM-MS relies heavily on DDA spectral libraries. Therefore, they share the same shortcomings, e.g., the dependence on the physical and chemical properties of proteins.
We confirmed from our HRM-MS results that the operons coordinate the gene expression more stringent than the randomized control in general. In addition, we found a multifaceted nature of the operon regulation: operons are not created equal. The stringency is length-dependent and functional-dependent at protein level. Such multifaceted organization of operons revealed two-level control: the operon unified transcriptional level and gene-specific translational level, which benefits the host in different aspects.
Although the operon organization maintains in general the stoichiometry of the genes in the operons compared to fully randomized scenario, the operons for metabolic pathways are in general less strictly controlled for stoichiometry balance compared to those operons for protein complexes. Protein complexes needs stoichiometry to maintain their functions. Therefore, the operons encoding protein complexes are tightly regulated to ensure the equal expression, such as the ribosome protein operons (Figure 5A). In contrast, since the operons for metabolic pathways are not necessarily forming a complex for their functions, their distinct specific activities set their specific demand in quantity. All quantified operons which encode enzymes in branch-free pathways and with available enzyme activity data validated this hypothesis without exception (Figure 5C). Nevertheless, bacteria need to regulate related metabolic pathways in quick response to environmental stimuli. Therefore, organizing the proteins in the same pathway under the control of one promoter would minimize the regulatory complexity of the adjustments, leaving the delicate control of each individual gene to the translational level. The available data indicated that such translational regulation is quite common (Figure 5C).
Our data also showed that the shorter operons, whose products form complexes, tend to be more tightly controlled in stoichiometric expression (Figure 4). This could be understandable such as two-component protein complexes would be invalid if there were imbalanced expression of the components. For instance, the assembly efficiency decreased remarkably if two subunits of bacterial luciferase LuxA and LuxB were split at distant chromosomal sites (). Large operons tend to encode proteins for large complexes such as ribosomes (Figure 4). Such large complexes would sustain for relatively long time in cells to perform essential functions. For example, half-life of mammalian ribosomes can be as long as 300 h (). Many protein components of ribosomes are dissociable and interchangeable with unbound counterparts (). This allows the rapid renewal of the damaged proteins of the complex without degrading and re-synthesizing the entire complex, which is the most energy-efficient way to keep these valuable complexes in good condition. This requires the delicate expression control of these proteins within one operon to meet the actual demand.
In another aspect, the differential expression regulation within an operon is also important for bacteria. Previous studies proposed that such regulation happens via generating different transcripts from multiple promotors/terminators [e.g., Bacillus subtilis dnaK operon (), Vibrio vulnificus putAP operon (), Zymomonas mobilis gap operon (), E. coli glpEGR operon ()], or via differential degrading mRNAs [e.g., Acinetobacter calcoaceticus mop operon (), E. coli atp operon (; )]. However, these studies included only individual cases of specific operons. Taking advantage of deep coverage of E. coli cytosolic proteome, our data indicated that enzyme activity seems to be an additional driving force for the differential expression regulation within an operon. Highly efficient enzymes tended to be less produced than the other counterparts in the same pathway. In such cases, deviating from stoichiometry minimizes the energy waste and thus may provide survival advantage. This explained the fact that no significant length-dependent stoichiometry is observed in “Pathway” proteins. Our analysis was restricted by the very limited availability of the enzyme activity and kinetics data. Validation using more such data is necessary in the future.
In this study, we noticed that the gene intervals in operons may serve as a regulatory factor for stoichiometry. It is a general accepted notion that termination of bacterial protein synthesis is obligatorily followed by recycling step governed by the factors ribosomal recycling factor (RRF), EF-G, and IF3, where the ribosome dissociates into its subunits (). In contrast, a recently described 70S-scanning mode of initiation holds that after termination, the 70S ribosomes do not dissociate after termination step but rather scan along with the mRNA until reaching the initiation site of the downstream cistron of the same mRNA (). Binding of fMet-tRNA triggers 70S scanning, which occurs in the absence of energy-rich compounds (e.g., ATP, GTP) and seems to be driven by unidimensional diffusion (). Therefore, the 70S scanning initiation might be a mechanism to read out the “stoichiometry code” of closely located genes in operons. In addition, the rate of translation of an ORF is controlled by a number of other mRNA features. For example, the codon selection and the cognate tRNA concentration dictate the translational pausing, which is a strong determinant of co-translational folding for most proteins (; ; ). Shine–Dalgarno (SD) sequence accessibility and strength have been implicated in translational initiation (). Genome-wide mRNA secondary structure analysis indicated that ORF translation rate is correlated with its mRNA structure in bacteria (), but not in mammalian cells (). Although highly stable mRNA structures in direct proximity to the initiation codon diminish translation efficiency (), secondary structure hiding the SD sequence in front of the second cistron prevents 30S binding initiation, but the secondary structure can be resolved by scanning 70S ribosomes when the secondary structure has a comparable stability (ΔG ≥ −6 kcal/mol at 30°C) (; ). These molecular mechanisms should be universal beyond exponential growth condition, although further validation would be needed.
This two-level regulation mode involving transcription and translation would balance the regulation in different time-scale. As transcriptional control takes effects at least in half an hour, it is suitable for sustained alteration in gene expression and in pathway-scale. In contrast, rapid and fine adjustment can be only performed at the translational level, which takes effects in less than 1 min and occurs at the individual gene level (). This rapid responsiveness would be also ideal for real-time adjustment of the proteins needed in complexes and pathways. Translational regulation largely contributes to the proteome complexity and minimizes the energy waste on synthesizing unnecessary proteins. Thus, the delicate and differential translational regulation in bacteria maintains the functionality and efficiency of both macromolecular complexes and metabolic pathways, which is a desperate need of the survival and competence of bacteria.
Statements
Data availability statement
The datasets generated for this study can be found in GEO, GSM3489376, GSM3489377.
Author contributions
GZ, JZ, and CS: conceptualization. JZ and HZ: experiment and visualization. JZ, HZ, BQ, and RN: data analysis. GZ and Q-YH: supervision. JZ, HZ, BQ, RN, Q-YH, CS, and GZ: writing. GZ: funding acquisition.
Funding
This work was supported by the Ministry of Science and Technology of China, National Key Research and Development Program (2017YFA0505001/2017YFA0505101/2018YFC0 910201/2018YFC0910202), Guangdong Key R&D Program (2019B020226001), and the Distinguished Young Talent Award of National High-level Personnel Program of China.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2019.00473/full#supplementary-material
FIGURE S1Physical and chemical features of quantified and unquantified proteins in E. coli. Panels (A–D) are distributions of protein isoelectric point, instability, protein length, and hydrophobicity. Instability tests a protein for stability, which value above 40 means the protein has a shorter half-life. The Kolmogorov–Smirnov test (KS-test) was used to test the distribution between quantified and unquantified proteins. Significance different were found between quantified and unquantified proteins indicated different physical–chemical characterizations of them.
FIGURE S2CV within operons of measured data and randomized negative control at RNA level.
FIGURE S3Length- and functional-dependence of operon stoichiometry control compare to randomized data. (A–E) Distribution among “Complex,” “Complex_random,” “Pathway,” and “Pathway_random” groups in BW25113. (F–J) Distribution among “Complex,” “Complex_random,” “Pathway,” “Pathway_random” groups in MG1655. In most cases, the CV in “Complex” and “Pathway” groups were lower than corresponding “Complex_random” and “Pathway_random” groups in both strains, indicating real robust signals but not statistical artifacts in operons present in E. coli.
FIGURE S4Length dependence of the stoichiometry control at RNA and protein levels. Same figure in MG1655 strain, comparable to the Figure 4A.
FIGURE S5Distribution of protein half-life CV in “Complex” and “Pathway” operons.
FIGURE S6Pathways in Figure 5.
TABLE S1DIA MS/MS isolation windows table.
TABLE S2Spectral library used in HRM-MS analysis.
TABLE S3Protein abundances and properties used in this analysis.
TABLE S4P-value of Mann–Whitney U test between complex and pathway operons shown in Figure 3C.
TABLE S5P-value of gene ontology (GO) overrepresentation of the operons with lower and higher CV groups shown in Figure 3D,E.
SUPPLEMENTARY SCRIPTSThe python scripts used in this study.
SUPPLEMENTARY DATAThe relationship of submitted raw data.
References
1
ArikeL.ValgepeaK.PeilL.NahkuR.AdambergK.ViluR. (2012). Comparison and applications of label-free absolute proteome quantification methods on Escherichia coli.J. Proteom.755437–5448. 10.1016/j.jprot.2012.06.020
2
BartholomausA.FedyuninI.FeistP.SinC.ZhangG.VallerianiA.et al (2016). Bacteria differently regulate mRNA abundance to specifically respond to various stresses.Philos. Trans. A Math. Phys. Eng. Sci.374:20150069. 10.1098/rsta.2015.0069
3
Ben-ShaharY.NannapaneniK.CasavantT. L.ScheetzT. E.WelshM. J. (2007). Eukaryotic operon-like transcription of functionally related genes in Drosophila.Proc. Natl. Acad. Sci. U.S.A.104222–227. 10.1073/pnas.0609683104
4
BeyerD.SkripkinE.WadzackJ.NierhausK. H. (1994). How the ribosome moves along the mRNA during protein synthesis.J. Biol. Chem.26930713–30717.
5
BloomJ. S.KhanZ.KruglyakL.SinghM.CaudyA. A. (2009). Measuring differential gene expression by short read sequencing: quantitative comparison to 2-channel gene expression microarrays.BMC Genom.10:221. 10.1186/1471-2164-10-221
6
BlumenthalT. (2004). Operons in eukaryotes.Brief. Funct. Genomic Proteom.3199–211. 10.1093/bfgp/3.3.199
7
BurkhardtD. H.RouskinS.ZhangY.LiG. W.WeissmanJ. S.GrossC. A. (2017). Operon mRNAs are organized into ORF-centric structures that predict translation efficiency.eLife6:e22037. 10.7554/eLife.22037
8
ChakrabortyA. R.DavidsonA.HowellP. L. (1999). Mutational analysis of amino acid residues involved in argininosuccinate lyase activity in duck delta II crystallin.Biochemistry382435–2443. 10.1021/bi982150g
9
ChapmanJ. D.GoodlettD. R.MasselonC. D. (2014). Multiplexed and data-independent tandem mass spectrometry for global proteome profiling.Mass Spectrom. Rev.33452–470. 10.1002/mas.21400
10
ChengY. S.ShenY.RudolphJ.SternM.StubbeJ.FlanniganK. A.et al (1990). Glycinamide ribonucleotide synthetase from Escherichia coli: cloning, overproduction, sequencing, isolation, and characterization.Biochemistry29218–227. 10.1021/bi00453a030
11
de CimaS.Gil-OrtizF.CrabeelM.FitaI.RubioV. (2012). Insight on an arginine synthesis metabolon from the tetrameric structure of yeast acetylglutamate kinase.PLoS One7:e34734. 10.1371/journal.pone.0034734
12
de SmitM. H.van DuinJ. (1990). Secondary structure of the ribosome binding site determines translational efficiency: a quantitative analysis.Proc. Natl. Acad. Sci. U.S.A.877668–7672. 10.1073/pnas.87.19.7668
13
DoerrA. (2015). DIA mass spectrometry.Nat. Methods1235–35. 10.1038/nmeth.3234
14
EddyC. K.MejiaJ. P.ConwayT.IngramL. O. (1989). Differential expression of gap and pgk genes within the gap operon of Zymomonas mobilis.J. Bacteriol.1716549–6554. 10.1128/jb.171.12.6549-6554.1989
15
FarhL.HwangS. Y.SteinraufL.ChiangH. J.ShiuanD. (2001). Structure-function studies of Escherichia coli biotin synthase via a chemical modification and site-directed mutagenesis approach.J. Biochem.130627–635. 10.1093/oxfordjournals.jbchem.a003028
16
FranklinF. C.VenablesW. A. (1976). Biochemical, genetic, and regulatory studies of alanine catabolism in Escherichia coli K12.Mol. Gen. Genet.149229–237. 10.1007/bf00332894
17
GaoY.JiangY.LiuQ.WangR.LiuX.LiuB. (2014). Enzymatic and regulatory properties of the trehalose-6-phosphate synthase from the thermoacidophilic archaeon Thermoplasma acidophilum.Biochimie101215–220. 10.1016/j.biochi.2014.01.018
18
GayS. C.FribourghJ. L.DonohoueP. D.SegelI. H.FisherA. J. (2009). Kinetic properties of ATP sulfurylase and APS kinase from Thiobacillus denitrificans.Arch. Biochem. Biophys.489110–117. 10.1016/j.abb.2009.07.026
19
Gil-OrtizF.Ramon-MaiquesS.Fernandez-MurgaM. L.FitaI.RubioV. (2010). Two crystal structures of Escherichia coli N-acetyl-L-glutamate kinase demonstrate the cycling between open and closed conformations.J. Mol. Biol.399476–490. 10.1016/j.jmb.2010.04.025
20
GordonS. P.TsengE.SalamovA.ZhangJ.MengX.ZhaoZ.et al (2015). Widespread polycistronic transcripts in fungi revealed by single-molecule mRNA sequencing.PLoS One10:e0132628. 10.1371/journal.pone.0132628
21
GraupnerM.XuH.WhiteR. H. (2002). New class of IMP cyclohydrolases in Methanococcus jannaschii.J. Bacteriol.1841471–1473. 10.1128/jb.184.5.1471-1473.2002
22
GruezA.Roig-ZamboniV.GriselS.SalomoniA.ValenciaC.CampanacciV.et al (2004). Crystal structure and kinetics identify Escherichia coli YdcW gene product as a medium-chain aldehyde dehydrogenase.J. Mol. Biol.34329–41. 10.1016/j.jmb.2004.08.030
23
GuellM.van NoortV.YusE.ChenW. H.Leigh-BellJ.MichalodimitrakisK.et al (2009). Transcriptome complexity in a genome-reduced bacterium.Science3261268–1271. 10.1126/science.1176951
24
HaftR. J.KeatingD. H.SchwaeglerT.SchwalbachM. S.VinokurJ.TremaineM.et al (2014). Correcting direct effects of ethanol on translation and transcription machinery confers ethanol tolerance in bacteria.Proc. Natl. Acad. Sci. U.S.A.111E2576–E2585. 10.1073/pnas.1401853111
25
HickmanJ.AshwellG. (1960). Uronic acid metabolism in bacteria. II. Purification and properties of D-altronic acid and D-mannonic acid dehydrogenases in Escherichia coli.J. Biol. Chem.2351566–1570.
26
HirokawaG.DemeshkinaN.IwakuraN.KajiH.KajiA. (2006). The ribosome-recycling step: consensus or controversy?.Trends Biochem. Sci.31143–149. 10.1016/j.tibs.2006.01.007
27
HomuthG.MogkA.SchumannW. (1999). Post-transcriptional regulation of the Bacillus subtilis dnaK operon.Mol. Microbiol.321183–1197. 10.1046/j.1365-2958.1999.01428.x
28
HuangD. W.ShermanB. T.TanQ.KirJ.LiuD.BryantD.et al (2007). DAVID bioinformatics resources: expanded annotation database and novel algorithms to better extract biology from large gene lists.Nucleic Acids Res.35W169–W175. 10.1093/nar/gkm415
29
IberD. (2006). A quantitative study of the benefits of co-regulation using the spoIIA operon as an example.Mol. Syst. Biol.2:43. 10.1038/msb4100084
30
IshihamaY.SchmidtT.RappsilberJ.MannM.HartlF. U.KernerM. J.et al (2008). Protein abundance profiling of the Escherichia coli cytosol.BMC Genom.9:102. 10.1186/1471-2164-9-102
31
KerrK. M.HedstromL. (1997). The roles of conserved carboxylate residues in IMP dehydrogenase and identification of a transition state analog.Biochemistry3613365–13373. 10.1021/bi9714161
32
KrishnaR. V.LeisingerT. (1979). Biosynthesis of proline in Pseudomonas aeruginosa. Partial purification and characterization of gamma-glutamyl kinase.Biochem. J.181215–222. 10.1042/bj1810215
33
KuznetsovaE.ProudfootM.GonzalezC. F.BrownG.OmelchenkoM. V.BorozanI.et al (2006). Genome-wide analysis of substrate specificities of the Escherichia coli haloacid dehalogenase-like phosphatase family.J. Biol. Chem.28136149–36161. 10.1074/jbc.M605449200
34
LandfaldB.StromA. R. (1986). Choline-glycine betaine pathway confers a high level of osmotic tolerance in Escherichia coli.J. Bacteriol.165849–855. 10.1128/jb.165.3.849-855.1986
35
LawK. P.LimY. P. (2013). Recent advances in mass spectrometry: data independent analysis and hyper reaction monitoring.Expert Rev. Proteom.10551–566. 10.1586/14789450.2013.858022
36
LawrenceJ. G.RothJ. R. (1996). Selfish operons: horizontal transfer may drive the evolution of gene clusters.Genetics1431843–1860.
37
LeeJ. H.ParkN. Y.LeeM. H.ChoiS. H. (2003). Characterization of the Vibrio vulnificus putAP operon, encoding proline dehydrogenase and proline permease, and its differential expression in response to osmotic stress.J. Bacteriol.1853842–3852. 10.1128/jb.185.13.3842-3852.2003
38
LiG. W.BurkhardtD.GrossC.WeissmanJ. S. (2014). Quantifying absolute protein synthesis rates reveals principles underlying allocation of cellular resources.Cell157624–635. 10.1016/j.cell.2014.02.033
39
LianX.GuoJ.GuW.CuiY.ZhongJ.JinJ.et al (2016). Genome-wide and experimental resolution of relative translation elongation speed at individual gene level in human cells.PLoS Genet.12:e1005901. 10.1371/journal.pgen.1005901
40
LinS.CronanJ. E. (2012). The BioC O-methyltransferase catalyzes methyl esterification of malonyl-acyl carrier protein, an essential step in biotin synthesis.J. Biol. Chem.28737010–37020. 10.1074/jbc.M112.410290
41
LiuW.XiangL.ZhengT.JinJ.ZhangG. (2017). TranslatomeDB: a comprehensive database and cloud-based analysis platform for translatome sequencing data.Nucleic Acids Res.46D206–D212. 10.1093/nar/gkx1034
42
LohmanJ. R.OlsonA. C.RemingtonS. J. (2008). Atomic resolution structures of Escherichia coli and Bacillus anthracis malate synthase A: comparison with isoform G and implications for structure-based drug discovery.Protein Sci.171935–1945. 10.1110/ps.036269.108
43
MacKintoshC.NimmoH. G. (1988). Purification and regulatory properties of isocitrate lyase from Escherichia coli ML308.Biochem. J.25025–31. 10.1042/bj2500025
44
MaierT.SchmidtA.GuellM.KuhnerS.GavinA. C.AebersoldR.et al (2011). Quantification of mRNA and protein and integration with protein turnover in a bacterium.Mol. Syst. Biol.7:511. 10.1038/msb.2011.38
45
MaoX.MaQ.ZhouC.ChenX.ZhangH.YangJ.et al (2014). DOOR 2.0: presenting operons and their functions through dynamic and integrated views.Nucleic Acids Res.42D654–D659. 10.1093/nar/gkt1048
46
McCarthyJ. E. (1990). Post-transcriptional control in the polycistronic operon environment: studies of the atp operon of Escherichia coli.Mol. Microbiol.41233–1240. 10.1111/j.1365-2958.1990.tb00702.x
47
McCarthyJ. E.GerstelB.SurinB.WiedemannU.ZiemkeP. (1991). Differential gene expression from the Escherichia coli atp operon mediated by segmental differences in mRNA stability.Mol. Microbiol.52447–2458. 10.1111/j.1365-2958.1991.tb02090.x
48
McLoughlinS. Y.CopleyS. D. (2008). A compromise required by gene sharing enables survival: implications for evolution of new enzyme activities.Proc. Natl. Acad. Sci. U.S.A.10513497–13502. 10.1073/pnas.0804804105
49
MiH.HuangX.MuruganujanA.TangH.MillsC.KangD.et al (2017). PANTHER version 11: expanded annotation data from Gene Ontology and Reactome pathways, and data analysis tool enhancements.Nucleic Acids Res.45D183–D189. 10.1093/nar/gkw1138
50
NeilsonK. A.AliN. A.MuralidharanS.MirzaeiM.MarianiM.AssadourianG.et al (2011). Less label, more free: approaches in label-free quantitative mass spectrometry.Proteomics11535–553. 10.1002/pmic.201000553
51
NikolovE. N.DabevaM. D.NikolovT. K. (1983). Turnover of ribosomes in regenerating rat liver.Int. J. Biochem.151255–1260. 10.1016/0020-711x(83)90215-x
52
OliverJ. C.GudihalR.BurgnerJ. W.PedleyA. M.ZwierkoA. T.DavissonV. J.et al (2014). Conformational changes involving ammonia tunnel formation and allosteric control in GMP synthetase.Arch. Biochem. Biophys.54522–32. 10.1016/j.abb.2014.01.004
53
OsbournA. E.FieldB. (2009). Operons.Cell. Mol. Life Sci.663755–3775. 10.1007/s00018-009-0114-3
54
PiH.LeeL. W.LoS. J. (2009). New insights into polycistronic transcripts in eukaryotes.Chang Gung Med. J.32494–498.
55
PriceM. N.HuangK. H.AlmE. J.ArkinA. P. (2005). A novel method for accurate operon predictions in all sequenced prokaryotes.Nucleic Acids Res.33880–892. 10.1093/nar/gki232
56
PurvineS.EppelJ. T.YiE. C.GoodlettD. R. (2003). Shotgun collision-induced dissociation of peptides using a time of flight mass analyzer.Proteomics3847–850. 10.1002/pmic.200300362
57
QinB.YamamotoH.UedaT.VarshneyU.NierhausK. H. (2016). The termination phase in protein synthesis is not obligatorily followed by the RRF/EF-G-dependent recycling phase.J. Mol. Biol.4283577–3587. 10.1016/j.jmb.2016.05.019
58
QiuX.TaoY.ZhuY.YuanY.ZhangY.LiuH.et al (2012). Structural insights into decreased enzymatic activity induced by an insert sequence in mannonate dehydratase from Gram negative bacterium.J. Struct. Biol.180327–334. 10.1016/j.jsb.2012.06.013
59
SatishchandranC.MarkhamG. D. (2000). Mechanistic studies of Escherichia coli adenosine-5′-phosphosulfate kinase.Arch. Biochem. Biophys.378210–215. 10.1006/abbi.2000.1841
60
SchaferT.StraussD.PetfalskiE.TollerveyD.HurtE. (2003). The path from nucleolar 90S to cytoplasmic 40S pre-ribosomes.EMBO J.221370–1380. 10.1093/emboj/cdg121
61
SchirmerF.HillenW. (1998). The Acinetobacter calcoaceticus NCIB8250 mop operon mRNA is differentially degraded, resulting in a higher level of the 3’ CatA-encoding segment than of the 5′ phenolhydroxylase-encoding portion.Mol. Gen. Genet.257330–337. 10.1007/pl00008621
62
SchmidtA.BeckM.MalmstromJ.LamH.ClaassenM.CampbellD.et al (2011). Absolute quantification of microbial proteomes at different states by directed mass spectrometry.Mol. Syst. Biol.7:510. 10.1038/msb.2011.37
63
SchubertO. T.LudwigC.KogadeevaM.ZimmermannM.RosenbergerG.GengenbacherM.et al (2015). Absolute Proteome composition and dynamics during dormancy and resuscitation of Mycobacterium tuberculosis.Cell Host Microbe1896–108. 10.1016/j.chom.2015.06.001
64
SchwanhausserB.BusseD.LiN.DittmarG.SchuchhardtJ.WolfJ.et al (2011). Global quantification of mammalian gene expression control.Nature473337–342. 10.1038/nature10098
65
SeshasayeeA. S.BertoneP.FraserG. M.LuscombeN. M. (2006). Transcriptional regulatory networks in bacteria: from input signals to output responses.Curr. Opin. Microbiol.9511–519. 10.1016/j.mib.2006.08.007
66
ShiehY. W.MinguezP.BorkP.AuburgerJ. J.GuilbrideD. L.KramerG.et al (2015). Operon structure and cotranslational subunit association direct protein assembly in bacteria.Science350678–680. 10.1126/science.aac8171
67
SteitzJ. A.JakesK. (1975). How ribosomes select initiator regions in mRNA: base pair formation between the 3’ terminus of 16S rRNA and the mRNA during initiation of protein synthesis in Escherichia coli.Proc. Natl. Acad. Sci. U.S.A.724734–4738. 10.1073/pnas.72.12.4734
68
TaylorA. M.FarrarC. E.JarrettJ. T. (2008). 9-Mercaptodethiobiotin is formed as a competent catalytic intermediate by Escherichia coli biotin synthase.Biochemistry479309–9317. 10.1021/bi801035b
69
TeichmannS. A.BabuM. M. (2004). Gene regulatory network growth by duplication.Nat. Genet.36492–496. 10.1038/ng1340
70
TurbevilleT. D.ZhangJ.AdamsW. C.HunterG. A.FerreiraG. C. (2011). Functional asymmetry for the active sites of linked 5-aminolevulinate synthase and 8-amino-7-oxononanoate synthase.Arch. Biochem. Biophys.511107–117. 10.1016/j.abb.2011.05.002
71
WisniewskiJ. R.RakusD. (2014). Multi-enzyme digestion FASP and the ’Total Protein Approach’-based absolute quantification of the Escherichia coli proteome.J. Proteom.109322–331. 10.1016/j.jprot.2014.07.012
72
WolfY. I.RogozinI. B.KondrashovA. S.KooninE. V. (2001). Genome alignment, evolution of prokaryotic genome organization, and prediction of gene function using genomic context.Genome Res.11356–372. 10.1101/gr.161901
73
WuD.HuT.ZhangL.ChenJ.DuJ.DingJ.et al (2008). Residues Asp164 and Glu165 at the substrate entryway function potently in substrate orientation of alanine racemase from E. coli: enzymatic characterization with crystal structure analysis.Protein Sci.171066–1076. 10.1110/ps.083495908
74
XieX.ZubarevR. A. (2015). Isotopic resonance hypothesis: experimental verification by Escherichia coli growth measurements.Sci. Rep.5:9215. 10.1038/srep09215
75
YamamotoH.WittekD.GuptaR.QinB.UedaT.KrauseR.et al (2016). 70S-scanning initiation is a novel and frequent initiation mode of ribosomal translation in bacteria.Proc. Natl. Acad. Sci. U.S.A.113E1180–E1189. 10.1073/pnas.1524554113
76
YangB.LarsonT. J. (1998). Multiple promoters are responsible for transcription of the glpEGR operon of Escherichia coli K-12.Biochim. Biophys. Acta1396114–126. 10.1016/s0167-4781(97)00179-6
77
YunH.LeeJ. W.JeongJ.ChungJ.ParkJ. M.MyoungH. N.et al (2007). EcoProDB: the Escherichia coli protein database.Bioinformatics232501–2503. 10.1093/bioinformatics/btm351
78
ZaslaverA.MayoA. E.RosenbergR.BashkinP.SberroH.TsalyukM.et al (2004). Just-in-time transcription program in metabolic pathways.Nat. Genet.36486–491. 10.1038/ng1348
79
ZhangG.HubalewskaM.IgnatovaZ. (2009). Transient ribosomal attenuation coordinates protein synthesis and co-translational folding.Nat. Struct. Mol. Biol.16274–280. 10.1038/nsmb.1554
80
ZhangW.ChenX.YanZ.ChenY.CuiY.ChenB.et al (2017). Detergent-insoluble proteome analysis revealed aberrantly aggregated proteins in human preeclampsia placentas.J. Proteome Res.164468–4480. 10.1021/acs.jproteome.7b00352
81
ZhaoP.ZhongJ.LiuW.ZhaoJ.ZhangG. (2017). Protein-level integration strategy of multiengine MS spectra search results for higher confidence and sequence coverage.J. Proteome Res.164446–4454. 10.1021/acs.jproteome.7b00463
82
ZhengY.AntonB. P.RobertsR. J.KasifS. (2005). Phylogenetic detection of conserved gene clusters in microbial genomes.BMC Bioinform.6:243. 10.1186/1471-2105-6-243
83
ZhongJ.XiaoC.GuW.DuG.SunX.HeQ. Y.et al (2015). Transfer RNAs mediate the rapid adaptation of Escherichia coli to oxidative stress.PLoS Genet.11:e1005302. 10.1371/journal.pgen.1005302
Summary
Keywords
operons, proteome quantification, HRM-MS, multifaceted stoichiometry control, mass-spectrometry, DIA, translation
Citation
Zhao J, Zhang H, Qin B, Nikolay R, He Q-Y, Spahn CMT and Zhang G (2019) Multifaceted Stoichiometry Control of Bacterial Operons Revealed by Deep Proteome Quantification. Front. Genet. 10:473. doi: 10.3389/fgene.2019.00473
Received
24 February 2019
Accepted
01 May 2019
Published
24 May 2019
Volume
10 - 2019
Edited by
Michael Ibba, The Ohio State University, United States
Reviewed by
Jesse Rinehart, Yale University, United States; Michelle Gibbs, The Ohio State University, United States
Updates
Copyright
© 2019 Zhao, Zhang, Qin, Nikolay, He, Spahn and Zhang.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Christian M. T. Spahn, christian.spahn@charite.deGong Zhang, zhanggong-uni@qq.com
This article was submitted to RNA, a section of the journal Frontiers in Genetics
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.