A novel proposal of a simplified bacterial gene set and the neo-construction of a general minimized metabolic network.

A minimal gene set (MGS) is critical for the assembly of a minimal artificial cell. We have developed a proposal of simplifying bacterial gene set to approximate a bacterial MGS by the following procedure. First, we base our simplified bacterial gene set (SBGS) on experimentally determined essential genes to ensure that the genes included in the SBGS are critical. Second, we introduced a half-retaining strategy to extract persistent essential genes to ensure stability. Third, we constructed a viable metabolic network to supplement SBGS. The proposed SBGS includes 327 genes and required 431 reactions. This report describes an SBGS that preserves both self-replication and self-maintenance systems. In the minimized metabolic network, we identified five novel hub metabolites and confirmed 20 known hubs. Highly essential genes were found to distribute the connecting metabolites into more reactions. Based on our SBGS, we expanded the pool of targets for designing broad-spectrum antibacterial drugs to reduce pathogen resistance. We also suggested a rough semi-de novo strategy to synthesize an artificial cell, with potential applications in industry.

A minimal gene set (MGS) is defined as the smallest possible gene set necessary and sufficient to maintain a living organism1. The MGS represents the infrastructure of a minimal cell and could be taken as a foundation for synthetic biology to create strains with desirable traits23. MGS research is biologically significant because of the following reasons: (i) it can further our understanding of the origin and evolution of life - for example, it can aid in determining the last universal common ancestor (LUCA)14; (ii) a pathogen’s MGS can guide the development of antibiotics35; (iii) reducing metabolic pathways and downsizing the genome could result in useful industrial strains6; and, most importantly, (iv) it allows the construction of a minimal genome that can be used as a basis for synthetic biology3789. Numerous researchers have contributed to the study of determining MGS10. Mycoplasma genitalium, a free-living bacterium with the smallest gene repertoire among the organisms sequenced to date, is an ideal species for MGS research and synthetic biology11 and has become the first genome to be subjected to genome-scale gene essentiality screening12. Mushegian and Koonin pioneered identification of an MGS by cross-species comparison, and they defined the first MGS, which contains ~256 conserved genes shared by the Gram-negative bacterium Haemophilus influenzae and the Gram-positive bacterium M. genitalium13. Later, Gil and colleagues reported a core MGS of 206 genes, after performing a comprehensive study based on the comparative genomics analysis of all available reduced genomes and experimental essentiality studies published at that time14. A minimal metabolism chart was also proposed. The minimal metabolic machinery is that which is indispensable to sustain life. Many studies devoted to determining the minimal (or core) metabolic network have been published. Burgard and colleagues approximated the core metabolic reactions of Escherichia coli under different culture conditions15, and Pollack and colleagues determined a list of minimum enzymatic reactions by reviewing the metabolic activities of various Mycoplasma species16. Castellanos et al. modeled purine and pyrimidine metabolic pathways in a minimal E. coli cell17. Recently, Barve et al. identified 133 “absolutely superessential” reactions in the E. coli metabolic network18, and Yang et al. defined a core metabolic proteome in an E. coli model19. Gil and colleagues explored the stoichiometric consistency and some architectural properties of the minimal metabolic network proposed by them in 2004 on rich environment20. Recently, Gil et al. updated their MGS by adding some genes to improve cell viability and new genes for RNA processing and metabolism2122.

Gil et al. stated that the research of MGS should consider gene essentiality22. In recent years, the increasing availability of sequenced genomes and experimentally determined essential genes have prompted an update of the MGS. Thus, this project aimed to develop a novel proposal of simplifying bacterial gene set to approximate a bacterial MGS by combining metabolic network construction23 with comparative genomics based on known essential genes. In the proposal, a simplified bacterial gene set (SBGS) that preserves both the self-reproduction and self-maintenance systems is determined. We believe that this work provides valuable information for drug design and is a useful reference for building a novel bacterial chassis in synthetic biology.

Results and Discussions

Obtaining the initial SBGS by comparative genomics and neo-construction of the metabolic network

In this work, we obtained the initial SBGS by comparative genomics and supplemented it by neo-construction of a bacterial approximately minimal metabolic network. Persistent essential genes (PEGs) can be considered the backbone genes for all bacterial organisms. PEGs can be considered the backbone genes for all bacterial organisms. As described in the methods, we gathered 611 PEGs from the CEG database (Cluster of Essential Genes)24 (Supplementary File S1). Among them, 598 (97.88%) and 508 (83.14%) genes were determined to be functionally similar to the genes in E. coli MG 1655 and Bacillus subtilis 168, respectively. In PEGs, 11 genes were annotated as “hypothetical proteins”. All cluster sizes of these hypothetical proteins were three or four. Six genes with cluster size = 3 were annotated as “putative function”. Genes with larger cluster size were annotated as defined function. In total, 594 (97.22%, 611-11-6) genes with definite functions were re-annotated and distributed in 22 subsystems (Supplementary Table S1).

In contrast to previous studies wherein homologous genes are required to be conserved in all reference species, we developed a new strategy named half-retaining, which requires the homologous genes to be present in more than one half of the referenced species to compile a highly persistent essential gene set (HPEGS). The HPEGS) with 248 essential genes that persist in more than half (cluster size ≥ 15/2 ≈ 8) of the species in the CEG database was obtained (Supplementary Table S2).

We used the bootstrap method to demonstrate the stability of half-retaining strategy. The results are shown in Supplementary Tables S3 and S4. Regarding the gene identity, the overlap of the genes between HPEGS14 and HPEGS ranged from 87% to 95.9% (mean = 91.7% and variance = 0.0009%). The absence of high-level clades led to the lowest overlap, as seen in group 2 (i.e., deleting B. subtilis 168 led to no firmicutes phylum in the reference data). The gene numbers in each HPEGS were compared, and the overlapping rates were higher than the gene content, which ranged from 88.5% to 96.2% (mean = 92.4% and variance = 0.0006%). Thus, the SBGS is basically invariable, and may only exhibit slight destabilization when the reference set is changed. On the contrary, previous proposals included in the MGS only those genes conserved in all the reference species and it approaches a null set when the number of reference species increases gradually. As an example, an MGS of 256 genes was obtained in previous work25, but the number of genes was drastically reduced to 63 when 100 genomes were compared and was reduced to zero when using 1000 genomes25. Our loose definition has the advantage that if the reference species are randomly (with no bias and covering most major lineages) selected, the gene set remains stable regardless of the number of reference species used.

Metabolism is essential for organisms to sustain life. To artificially synthesize a cell, the minimal metabolism of a bacterial cell should be considered. Therefore, we constructed an approximate MMN to define the core genes involved in metabolism in the MGS. The top-down approach in synthetic biology is frequently used to downsize the function of the object26. Researchers can identify the essential components and functions of a cell based on its minimal function.

Therefore, on the next step to construct SBGS, we submitted the 594 genes included in the PEG set into Model_SEED to define a bacterial general metabolic network (GMN1). The GMN1 contains 555 reactions and 324 genes. By adding 84 essential genes (participating in 120 reactions), a viable bacterial metabolic model (GMN2) was obtained that encompassed 408 genes and 675 reactions.

Ultimately, 251 genes were removed from GMN2, resulting in a minimized bacterial metabolic network (MMN) with 157 genes that were involved in 431 reactions (Supplementary Fig. S1, Supplementary Tables S5 and 6 and Supplementary File S2). We only deleted genes with a cluster size <8, as well as their associated reactions, when we minimized the network. Hence, several redundant reactions associated with genes of cluster size >8 remain in our proposed bacterial MMN, and the bacterial MMN is approximately minimal. The redundancy of these reactions could increase the robustness of the network.

Consequently, by neo-constructing the approximate bacterial MMN, we obtained a gene set of 157 genes. Among these, 91 overlap with the collection of HPEGS. Hence, a SBGS including 314 (248 + 157 − 91) genes was obtained by the union of HPEGS and genes associated with the approximate bacterial MMN (Supplementary Table S2).

We analyzed the gene proportions involved in the metabolism in 40 species; the resulting distribution is shown in Fig. 1a. We found that the value of the SBGS is close to the proportion of those species with small genomes in the reference genome list used.

Figure 1

Comparison of the genes in our MGS and others.

(a). The distribution of genes involved in the metabolism of each species and in the minimal gene set. (b) A Venn diagram for the three MGSs showing that our MGS contains 91% of the genes (128 of 141) existing in both Koonin et al.’s and Gil et al.’s MGSs, as well as our 107 newly identified genes (underlined in Supplementary Table S2).

In the SBGS, 66 genes with a cluster size <8 were added by constructing an approximate bacterial MMN, which corresponded to two extreme cases. The first is that the included genes have cluster sizes close to 8 (e.g., 7, 6 and 5; 38%). The second is that the included genes have cluster sizes close to 0 (e.g., cluster size = 0, 1, 2; 46%). The first case might be caused by the limitations of the techniques used to determine essential genes in difficult situations27 in which some factually essential genes might be overlooked in a few organisms. Genes functioning in a synergetic network may explain the second case. For example, three genes related to L-glutamine transport, glnH, glnP and glnQ, only have cluster sizes of 0 or 1. However, they all cooperate with the gene glnS, which has a cluster size of 8 and participates in the “GLNabcpp” reaction (Supplementary Table S5). The inactivation of glnH, glnP and glnQ may not interrupt the “GLNabcpp” reaction (thus resulting in a cluster size far less than 8); however, if the four genes glnH, glnP, glnQ and glnS are inactivated together, the “GLNabcpp” reaction is interrupted. Because the “GLNabcpp” reaction is of critical importance and has a super-essentiality = 1, the glnS gene is regarded as essential and the other three as nonessential. This reason leads to the existence of 18 nonessential genes in the bacterial MGS. However, they are involved in essential reactions of the approximate bacterial MMN. Kemmeren and Holstege confirmed that observing the effect of deleting a single gene might not precisely reveal all essential genes28 because the genes perform functions within a synergistic network (or function), and one defective gene may be substituted for by others performing similar functions.

After comparing our SBGS with previous MGS studies, we identified 169 and 166 genes that overlapped with sets reported by Mushegian and Koonin13. and Gil et al.14, respectively, as shown in Fig. 1b. Remarkably, our SBGS includes 91% (128) of the genes that exist in both of previous sets. The remaining 13 genes that exist in both Mushegian and Koonin’s and Gil et al.’s MGSs were lost in our SBGS (Supplementary Table S7). In order to ensure obtaining a SBGS as complete as possible, the genes existing in both Mushegian and Koonin’s and Gil et al.’s MGSs are considered as reliably essential in MGS. Hence, we added the 13 genes into our SBGS. As a result, our SBGS includes 314 + 13 = 327 genes. Moreover, we identified 107 new genes in our SBGS, including 62 genes identified by the half-retaining strategy and 45 genes identified using the approximate MMN construction method. Although our set might still be incomplete to perform all essential cell functions, the other two sets must be more affected by this issue because they were designed mainly based on ubiquitous genes. Henry and coworkers remarks that Gil et al.’s MGS omitted important genes involved in physiological processes, such as dnaA (replication initiation factor) and dnaC (a loading factor for helicase dnaB)29. The work of Henry et al. found that both Mushegian and Koonin’s and Gil et al.’s sets omitted 40 genes involved in RNA processing, metabolism and translation29. In contrast, our half-retaining strategy loosens the request of retaining in the final set and hence could contain more genes encoding essential functions. Supplementing it with those genes needed for a viable metabolic network makes our set more reliable although it may not be minimal.

We used COG (Clusters of Orthologous Groups)30 to categorize the functions of genes in different genes sets (Table 1). Table 1 demonstrates that the HPEGS covers 17 COG categories. Genes related to metabolism cover 17 COG categories, and two additional COG categories were added to the MGS (i.e., signal transduction mechanisms and defense mechanisms). After supplementing the genes by the neo-construction of an approximate bacterial MMN, and similarly to the Koonin’s proposal, our MGS includes 19 categories of COGs categories in total.

Table 1

The number of genes in each COG categories of different gene sets.

COG category	Number of gene in gene sets
	HPEGS	MMN	SBGS	Koonin	Gil
Energy production and conversion (C)	12	13	16	17	9
Cell division and chromosome partitioning (D)	5	0	5	2	1
Amino acid transport and metabolism (E)	5	23	24	12	5
Nucleotide transport and metabolism (F)	13	15	14	17	15
Carbohydrate transport and metabolism (G)	11	17	18	15	17
Coenzyme metabolism (H)	14	19	19	8	12
Lipid metabolism (I)	25	28	29	4	7
Translation, ribosomal structure and biogenesis (J)	98	8	98	93	98
Transcription (K)	6	1	7	11	8
DNA replication, recombination and repair (L)	19	1	19	24	16
Cell envelope biogenesis, outer membrane (M)	20	18	25	8	2
Posttranslational modification, protein turnover, chaperones (O)	8	8	10	14	14
Inorganic ion transport and metabolism (P)	1	7	7	9	2
Secondary metabolites biosynthesis, transport and catabolism (Q)	3	3	3	1	0
General function prediction only (R)	13	1	13	11	0
Function unknown (S)	1	1	1	3	0
Signal transduction mechanisms (T)	0	2	2	1	2
Intracellular trafficking, secretion, and vesicular transport (U)	9	0	9	5	0
Defense mechanisms (V)	0	5	5	2	0

In summary, we have compiled a SBGS with 327 genes covering 19 COG categories sufficient to perform essential cellular functions by combining comparative genomics and metabolic network neo-construction. SBGS is a theoretical model and does not correspond to any existing species. The size of SBGS is consistent with the lower limit of theoretical gene numbers, as Koonin stated that the number of genes in the MGS would likely be in the range of 300–350 in nutrient-rich medium4.

Nevertheless, because MGS depends on the reference genomes used to extract it, there may not be a unique solution for a bacterial MGS. Our proposed SBGS provides an alternative reference for this issue.

Comparisons with previous works and assigning essentialities to reactions in the approximate bacterial MMN

Supplementary Table S6 lists the number of reactions in each subsystem for our MMN. We compared them with the hypothetical minimal metabolism of Gil et al., and found that all reactions in Gil et al. exist in our MMN. Due to the half-retaining strategy, our approximate MMN has some redundant reactions that, as above stated, will increase the robustness of the network.

Barve et al. pioneered a method to rank reaction essentiality and proposed a concept of “super-essentiality” to estimate the importance of all reactions in the E. coli metabolic network18. The super-essentiality ranges from 0 to 1. One reaction with super-essentiality = 0 means that it is non-essential in the network. 133 absolutely super-essential reactions were identified in this manner. We selected the common reactions between their work and our approximate bacterial MMN model to analyze the essentiality of the reactions.

As a result (Supplementary Fig. S2, Supplementary Table S5), 261 reactions (60.56%) were categorized as super-essential. Among these, 213 reactions (81.61%) had super-essentiality values larger than zero, and 101 reactions had a value of one (absolutely super-essential reactions). The latter group included most of the absolutely super-essential reactions reported in Barve’s work, demonstrating the consistency of reaction essentiality between our approximate bacterial MMN model and Barve’s previous work. However, 32 absolutely super-essential reactions were absent from our approximate bacterial MMN. The difference between ours and Barve’s work might be caused by the different nutrient compositions of the medium. Barve et al. identified super-essentialities based on different nutritional requirements, whereas we only used one type of medium (D-glucopyranose) when we performed Flux Balance Analysis (FBA).

Determining key metabolites and genes by analyzing the topology of the bacterial approximately MMN and the application on drug targets development

We identified 25 key metabolites according to the node connectivity of metabolites in the approximate bacterial MMN (Table 2). Jeong et al. found that the hub metabolites were similar among all organisms31. In fact, most of the hub metabolites identified are “current metabolites”, as defined by Ma and colleagues32. Among the 25 key metabolites identified in this work, 14 are “current metabolites” suggested by both Ma et al.32 and Jeong et al.31, whereas six are consistent with Ma et al.32 only. The other five key metabolites are newly identified in this work.

Table 2

Key metabolites in the network (out-degrees and in-degrees)a.

Metabolite name	Node connectivity	Metabolite description	Remarks
cpd00067	215	H	Current metabolite
cpd00001	122	H2O	Current metabolite
cpd11493	79	ACP	Current metabolite
cpd00002	73	ATP	Current metabolite
cpd00011	56	CO2	Current metabolite
cpd00008	54	ADP	Current metabolite
cpd00006	52	NADP	Current metabolite
cpd00005	49	NADPH	Current metabolite
cpd00012	45	PPi	Current metabolite
cpd00003	42	NAD	Current metabolite
cpd11492	40	Malonyl-acyl-carrierprotein	Ma
cpd00010	39	CoA	Ma
cpd00009	38	Phosphate	++
cpd00004	37	NADH	Current metabolite
cpd00080	33	Glycerol-3-phosphate	Ma
cpd00046	30	CMP	Current metabolite
cpd00052	19	CTP	Current metabolite
cpd00018	19	AMP	Current metabolite
cpd00033	16	Glycine	++
cpd00054	15	L-Serine	++
cpd00035	11	L-Alanine	++
cpd10516	11	fe3	++
cpd00023	10	L-Glutamate	Ma
cpd00020	9	Pyruvate	Ma
cpd00061	9	Phosphoenolpyruvate	Ma

aThe key metabolites in bold and italics marked with “++” are newly found in this work. The metabolites marked as “Current metabolite” are consistent with the work of Ma et al. (Ma & Zeng, 2003) and of Jeong et al. (Jeong et al.31). The metabolites marked as “Ma” are consistent with the work of Ma (Ma & Zeng,32).

Spearman correlation analysis showed that the node connectivity was significantly positively correlated with SBGS essentiality (p < 0.01, rho = 0.211; Table 3, see Methods). Both the average in-degree and out-degree of the metabolites were significantly negatively correlated with the bacterial MGS essentiality of a gene (p < 0.01, Table 3, see Methods). This result suggested that genes with a higher SBGS essentiality involve a lower average number of metabolites in a reaction.

Table 3

Spearman correlation between MGS essentiality and node connectivitya.

Object	rho	p-value
Essentiality ~ node connectivity	0.211	0.008198
Essentiality ~ average in-degree	−0.274	0.000512
Essentiality ~ average out-degree	−0.281	0.000369
Essentiality ~ average degree	−0.278	0.000424

aMGS essentiality is represented by cluster size; node connectivity is represented by the number of reactions associated with a gene; average in-degree is represented by the number of reactants divided by that of reactions associated with a gene; average out-degree is represented by the number of products divided by that of reactions associated with a gene; average degree is represented by the number of metabolites divided by that of reactions associated with a gene.

Looking back in history of drug discovery, we found that the highly essential genes are often used as effective drug targets. For example, Haydon et al. synthesized an antibacterial drug based on ftsZ, a cell division protein that is present in SBGS with a cluster size = 1433. Recently, Ravishankar et al. identified a target for anti-tubercular, topA (cluster size = 10)34. Tharinjaroen et al. found a novel target, lepB gene (cluster size = 9), for M. tuberculosis and M. bovis35. Another example is the fabG gene (cluster size = 13), which is associated with 38 reactions and acts as a 3-ketoacyl-(acyl-carrier-protein) reductase. According to our results, the fabG gene may be a good candidate target for developing antibacterial drugs. To confirm this idea, we searched the DrugBank database and found that it had already been used as a drug target36. These successful examples indicate that highly essential genes in SBGS have been used in previous antibacterial drug discovery programs. Meanwhile, we find the highly essential genes often are hub nodes in approximate bacterial MMN.

Currently, the rapid emergence of multidrug-resistant pathogens has led to the ineffectiveness of conventional antibiotics for combating super bacteria3738. Therefore, it is essential to develop new drugs to combat these pathogens. However, new drugs aimed at old targets may also meet with resistance. New drugs directed at new targets will be the most effective choice37. Considering that some targets are obsoleted prematurely, we must increase the pool of antibacterial drug targets39. We searched the DrugBank database for all genes in SBGS and found that 143 of 327 genes had been tagged as drug targets (Supplementary Table S8-1). This leaves 184 highly essential genes that have not yet been targeted by extant drugs (Supplementary Table S8-2). Among these, 102 genes do not have significant similarity with any human genes (BlastP E-value >10−3). These genes meet the most crucial criteria for broad spectrum antibacterial drug target selection: (1) a highly conserved function in a range of pathogens; (2) essentiality of the gene for the pathogens; and (3) no highly conserved function in humans40. Therefore, we suggest that these genes should be considered as targets in the development of new broad spectrum antibacterial drugs to expand the pool of targets for drug design.

The application of SBGS on synthesis of an artificial cell

Foley and Shuler noted that engineers are interested in synthetic biology to develop a self-replicating biological system41. Pohorille et al.42 and Rasmussen et al.43 proposed that a human-made system could be considered “living” if three criteria are met: self-maintenance (metabolism), self-reproduction, and the capacity for Darwinian evolution. The approximate bacterial MMN informed the creation of our SBGS. Indeed, our SBGS includes abundant genes related to DNA replication, translation, transcription and posttranslational modification, protein turnover, and chaperones. Thus, it preserves both self-reproduction and self-maintenance systems.

Previously, scientists have synthesized artificial cells and chromosomes in the wet laboratory444546 and have built computational simulation models894147484950. Shuler’s group created the first mathematical model50 and recently developed a minimal cell model that can be tested by chemically simulating the behavior of a whole cell9. Scientists at the JCVI have tried to define a minimal Mycoplasma genome by gene deletion techniques1251 and the chemically synthesized and assembled new a Mycoplasma genomes have been successfully introduced into a cell454651. There are two strategies to synthesize artificial cells at present, top-down (genome downsizing) and bottom-up (de novo synthesis)5253. Furthermore, Our SBGS could provide the third strategy for the synthesis of an applicable cell by the following procedure (Fig. 2). We suggest the semi-de novo synthesis of a cell, starting from M. genitalium. First, 101 genes of SBGS not present in M. genitalium would be integrated into its genome one by one (Supplementary Table S9). We would reference the gene order of the other species when transferring a new gene into the genome. After integrating one specific gene, we would verify the phenotype of the cell. If the cell does not thrive, we would abort the insertion of this gene into the genome. Thus, we would obtain an extensional genome with 475 + 101 = 576 genes. Second, we could knockout the 249 genes of M. genitalium that are absent in our MGS one by one from the extensional genome (Supplementary Table S10). Likewise, after knocking out a gene from the genome, we would determine its survival status. If the cell could not live normally, the gene would be retained. After completion of these steps, we would obtain an artificial cell with 327 genes that could live normally and could be regarded as a general bacterial chassis (Supplementary Table S11). Finally, we would supplement other genes into the chassis according to specific applications. Although it could provide a general reference of this issue, our proposal of semi-artificial bacterium might be vague in current form because the newly defined SBGS is just one additional theoretical proposal. A Mycoplasma has a quite special cell envelope and it may be not feasible to implement the proposal to simulate general metabolism that involves gram-positive and gram-negative bacteria. The minimal genome should not only include the MGS, but also contains the noncoding region such as noncoding RNA, UTR and gene control regions. For example, Serrano et al. have attempted to define the essential small ORFs and ncRNAs of a minimal cell54. Additionally, the gene order is highly variable among bacteria, and the chromosome architecture needs also to be taken into account. After all these points being addressed by us or the community in the near future, perhaps it will bring a novel bacterial species with great industrial applications.

Figure 2

Our design for the semi-de novo synthesis of an artificial cell based on our MGS.

(A) involves transferring the genes in the MGS but not in M. genitalium to the genome. (B) is knocking out the genes of M. genitalium that are absent from our MGS in the genome obtained in (A). (C) involves supplementing with genes required for specific applications.

Methods

Aiming at improving previous attempts to define a MGS, we proceeded through three steps (Fig. 3). First, we started from experimentally determined essential genes. Second, we developed a new strategy named half-retaining to identify essential genes conserved among over half of the reference species. Third, we supplemented our initial MGS of conserved essential genes by the neo-construction of a viable general metabolic network and subsequently downsizing it to an approximate minimal network.

Figure 3

The procedure of this work.

Data sources

The candidate essential genes were obtained from the CEG database24. The current version of CEG covers essential genes from 15 species of bacteria, which are listed in Supplementary Table S12. All clusters in the CEG correspond to essential genes for at least one reference species. Each cluster has a size value, indicating the number of reference species in which the corresponding gene is essential. Based on the cluster size, users can easily determine whether an essential gene is conserved in multiple species or is species-specific. Therefore, cluster size was used as a metric of conserved gene essentiality in the reference species and also of general gene essentiality in the SBGS. For example, the pgsA gene, with a cluster size of 15, is consistently essential in all 15 species. This result suggested that the pgsA gene is highly essential in the SBGS. The cluster size was used as a paramount index to determine the SBGS. To make the reconstructed model viable, we extracted a portion of the metabolic reactions from the iJR904 model55 and the iAF1260 model of E. coli56 (the two best annotated metabolic models) to fill the gaps in our reconstructed metabolic network.

The half-retaining strategy

In contrast to previous studies wherein homologous genes are required to be conserved in all reference species, we developed a new strategy named half-retaining, which requires the homologous genes to be present in more than one half of the referenced species to compile a highly persistent essential gene set (HPEGS). In this work, a gene is considered a persistent essential gene (PEG) if it is essential and shared by more than three reference species. The PEG is evolutionarily conserved and serves as the basis for the SBGS for bacterial life. Thus, we determined an initial SBGS via comparative genomics, using a half-retaining strategy to compile a HPEGS.

We used the bootstrap method to demonstrate the advantage of this new approach. For 15 organisms, one was picked out each time, and the remaining 14 species were used as reference species. Thus, we generated 15 groups of reference species. For each group, we used the half-retaining strategy to obtain an HPEGS14 (a highly conserved universal gene set based on the 14 retained species). Subsequently, the new HPEGS14 was compared with the HPEGS obtained based on all 15 reference species.

The construction of an approximate MMN

To construct a minimal metabolic network, the PEGs were re-annotated using the RAST tool57 of SEED58. The SEED annotator57 is based on the subsystems but not on sequence similarity. To estimate the reliability of the annotated genes, we compared them with E. coli MG1655 (Gram−) and B. subtitlis 168 (Gram+), which are two well-characterized model organisms, with the SEED tool. To downsize the genome and obtain an MGS with functional metabolic ability, we followed the procedure shown in Fig. 3.

First, GMN1 was neo-constructed based on the re-annotated PEG identified by the RAST tool of SEED. However, the GMN1 had gaps and was a dead network, lacking effective flux according to the FBA. The gaps were filled using the metabolic reactions extracted from the two metabolic models JR904 and iAF1260 of E. coli, which may be the most complete metabolic models thus far, until the viable (i.e., with effective flux in the biomass reaction) universal metabolic network GMN2 was obtained. To minimize the gene number, we deleted the disabled genes and disconnected reactions (i.e., those genes that could be knocked out and leave a viable network) from GMN2.

To refine the metabolic network according to the half-retaining strategy, genes with a cluster size ≥ 8 were reserved and regarded as the skeleton genes of the network, whereas those genes with a cluster size < 8 in GMN2 were deleted one by one in order of size. The temporary network was submitted to the MetaNetX website59 for flux analysis. We used growth medium with D-glucopyranose as the sole carbon source and rich in hydrogenphosphate, ammonium, water, proton, oxygen, and carbon dioxide for flux analysis. After single gene knockout from GMN2, if the flux is zero, the corresponding gene was retained in the network; otherwise it was dropped. New temporary networks were submitted to Model_Seed to regenerate a new biomass function to adapt to the new network.

Calculation of topological properties of the defined metabolic network

To identify the key metabolites in the MMN, we used the total number of reactions affected by a metabolite as its connectivity. For example, 8 reactions used metabolite M as substrate or product; hence, the connectivity of M is 8. Attacking at the hub nodes could paralyze the network. Thus, the key metabolites play important roles in the network.

In addition to the analysis of key metabolites, key genes were also analyzed. We proposed a criterion that a gene is more essential in the MGS if it is essential in more individual genomes. Based on this criterion, we used “cluster size” in CEG to indicate the essentiality of a gene in the MGS. The number of reactions associated with a gene was defined as the node connectivity of the gene. Furthermore, to determine whether the average metabolite number of multiple reactions of a gene relates to its MGS essentiality, we analyzed the average in-degree and out-degree of the relevant reactions. For example, if a gene is associated with 5 reactions and its cluster size is 9. These 5 reactions have 25 substrates as well as 35 products totally. So its MGS essentiality is 9, its connectivity is 5, the average in-degree is 25/5 = 5 and average out-degree is 35/5 = 7.

Additional Information

How to cite this article: Ye, Y.-N. et al. A novel proposal of a simplified bacterial gene set and the neo-construction of a general minimized metabolic network. Sci. Rep. 6, 35082; doi: 10.1038/srep35082 (2016).