- Research article
- Open Access
The genetic control of growth rate: a systems biology study in yeast
© Pir et al; licensee BioMed Central Ltd. 2012
- Received: 3 January 2012
- Accepted: 13 January 2012
- Published: 13 January 2012
Control of growth rate is mediated by tight regulation mechanisms in all free-living organisms since long-term survival depends on adaptation to diverse environmental conditions. The yeast, Saccharomyces cerevisiae, when growing under nutrient-limited conditions, controls its growth rate via both nutrient-specific and nutrient-independent gene sets. At slow growth rates, at least, it has been found that the expression of the genes that exert significant control over growth rate (high flux control or HFC genes) is not necessarily regulated by growth rate itself. It has not been determined whether the set of HFC genes is the same at all growth rates or whether it is the same in conditions of nutrient limitation or excess.
HFC genes were identified in competition experiments in which a population of hemizygous diploid yeast deletants were grown at, or close to, the maximum specific growth rate in either nutrient-limiting or nutrient-sufficient conditions. A hemizygous mutant is one in which one of any pair of homologous genes is deleted in a diploid, These HFC genes divided into two classes: a haploinsufficient (HI) set, where the hemizygous mutants grow slower than the wild type, and a haploproficient (HP) set, which comprises hemizygotes that grow faster than the wild type. The HI set was found to be enriched for genes involved in the processes of gene expression, while the HP set was enriched for genes concerned with the cell cycle and genome integrity.
A subset of growth-regulated genes have HFC characteristics when grown in conditions where there are few, or no, external constraints on the rate of growth that cells may attain. This subset is enriched for genes that participate in the processes of gene expression, itself (i.e. transcription and translation). The fact that haploproficiency is exhibited by mutants grown at the previously determined maximum rate implies that the control of growth rate in this simple eukaryote represents a trade-off between the selective advantages of rapid growth and the need to maintain the integrity of the genome.
- Relative Growth Rate
- Dilution Rate
- Gene Copy Number
- Dosage Compensation
- Maximum Specific Growth Rate
Appropriate control of the rate of cell growth is central to the long-term survival of species, particularly microorganisms. A fast growth rate is a competitive advantage when environmental conditions are favourable, while slow growth or even quiescence may allow survival under stress conditions such as rapid changes in the physicochemical environment, starvation, or exposure to toxins. Thus most organisms have evolved stringent controls over the rate of cell growth, and any disruption of these controls has severe consequences; an obvious example being the development of cancer as a result of uncontrolled cell division.
Our genome-scale approach to the investigation of growth rate control in the model eukaryote cell, S. cerevisiae, is based on the concepts of Metabolic Control Analysis (MCA; ). MCA defines the extent of control exerted over the flux through a pathway by its components. In our context, flux is the growth rate of cells, pathway components are the gene products, and control coefficients are the ratios of the fractional changes in growth rate to the fractional changes in the concentrations of gene products. Thus the term "control" has a special meaning in the context of MCA. Those cell components that exert control over flux are not necessarily regulatory molecules. Rather, they exert control because the flux through the pathway is sensitive to changes in the concentration (or activity) of the component. In this paper, we shall routinely use the term 'growth rate', rather than 'flux', since we have measured differences in growth rate consequent on reducing the copy number of individual genes from two to one in diploid yeast.
As part of our overall approach, we have designed experiments in two categories to elucidate the control coefficients of gene products for growth rate. In Category 1 experiments, we measured the change in the concentration of gene products as a response to a change in flux (;). We report no new Category 1 experiments in the present study, but we shall refer to our earlier experiments in this Category. In Category 2 experiments, the concentration of gene products is altered and any impact on growth rate measured ; all new data reported in this paper relates to Category 2 experiments.
In our earlier Category 1 experiments, we grew a reference yeast strain in carbon-, nitrogen-, phosphorus-, or sulphur-limited conditions in chemostat culture at dilution rates, D = 0.07, 0.1 and 0.2 h-1 (equivalent to doubling times ~ 10, 7, and 3.5 h, respectively), and collected samples at each of the different steady states. We analysed our samples for changes in the levels of mRNAs, proteins, and metabolites with respect to dilution rate and identified a set of growth-rate-regulated (GR) genes; i.e. a set of genes that significantly changed their expression levels in response to the change in growth rate, irrespective of the specific nutrient whose rate of supply determined the rate of growth . There is also a nutrient-regulated (NR) gene set, whose expression levels change according to the nutrient whose rate of supply is determining the growth rate; again, we have reported on these previously . The level of expression of NR genes may vary in either a growth-rate-dependent or growth-rate-independent manner.
In the set of Category 2 experiments that we reported previously, we established a growth-rate competition between yeast deletion mutants (each hemizygous for just one of the organisms 5,800 protein-encoding genes ) under carbon, nitrogen, or phosphorus limitation in chemostats at D = 0.1 h-1. We classified all the mutants according to the sign of their relative growth rate, the mutants with negative relative growth rate (i.e. whose proportion in the population fell significantly overtime) were classified as haploinsufficient (HI) and the mutants with positive relative growth rate (i.e. whose proportion in the population significantly increased with time) were classified as haploproficient (HP). The set of haploinsufficient and haploproficient genes together form the set of high flux control (HFC) genes, i.e. those genes for which a reduction of their copy number in diploid cells from 2 to 1 results in a significant change in growth rate (flux).
Our earlier results showed that, under all three nutrient limitations and at a dilution rate of 0.1 h-1, there was little overlap between HFC genes from Category 2 experiments and the GR genes defined in Category 1 experiments; i.e. the genes regulated by growth rate are not, themselves, regulators of growth. At the time , we noted that this result might, in the jargon of Systems Biology, represent a 'design rule' for the eukaryotic cell. However, since cells are the products of evolution rather than design, such rules may change according to the selection pressures to which an organism is exposed. Thus, while this rule about HFC genes holds for nutrient-limited environments, further studies are required to determine whether it has any greater generality. This paper reports these additional studies.
Since our previous set of Category 2 experiments was carried out at a low growth rate and under nutrient limitation, we decided to determine whether the lack of overlap between HFC and GR gene sets still held at high growth rates or under nutrient-unconstrained conditions. Chemostats tend to be unstable at dilution rates close to the cells' maximum specific growth rate (μmax) and so, to establish a steady state at μmax, we used turbidostats to allow us to monitor competition between the pool of hemizygous yeast deletants growing in a complex synthetic medium (FPM ; see Methods section for modifications). Turbidostats are continuous cultures in which the cells, rather than the experimenter, control the growth rate  since the nutrient supply is determined by the biomass concentration in the growth vessel. Accordingly, we established a turbidostat in which the biomass concentration was held at a value equivalent to that of a mid-exponential phase batch culture. The turbidostat is unconstrained by the supply of nutrients and so equilibrates at the maximum specific growth rate of the yeast strain used. However, if the culture is a pool of hemizygous mutants of S. cerevisiae, and if the different mutants in the pool can have different maximum growth rates, those mutants that can achieve a μmax greater than the population average will increase in the population over time, while those with a μmax less than the population average will decrease in the population. (Those readers unfamiliar with chemostats and turbidostats will find a fuller explanation in Additional File 1, Additional Text; also see reviews by Bull  and Pirt .)
Turbidostats, just like chemostats, represent a sensitive way of identifying haploinsufficient and haploproficient phenotypes, with the difference that (in a turbidostat) haploproficient mutants must be capable of growing at a rate greater than the previously recorded μmax. Thus, in this study, we have searched for sets of genes that exhibit either haploinsufficient or haploproficient phenotypes (these are both sub-sets of the class of HFC genes) in rapidly growing yeast cultures, either in nutrient-limited chemostats operated at D-values close to μmax or in the nutrient-unconstrained conditions of turbidostat culture. Further, we investigated subsets of HFC genes that have closely related cellular functions in order to elucidate the gene and protein features that mediate their high degree of control over growth rate. Finally, we propose an extended approach to growth rate (or flux) control, without limiting its definition to enzymes or haploinsufficiency.
Genes showing high growth-rate control during rapid growth in nutrient-unconstrained conditions
The mutants were then classified according to the sign of their relative growth rate. 796 mutants were identified with a positive slope (FDR < 0.05). These grew faster than the population average and are defined as haploproficient. 1932 mutants grew slower than the population average (FDR < 0.05) and are defined as haploinsufficient. In total, 2728 HFC genes were found for which the changes in the relative concentrations of their hemizygotes were not significant (see Additional File 1, Additional Text; Additional File 2, Figure S1; Additional File 3, Table S1).
GO terms enriched among HFC control (see Additional File 4, Table S2 for complete list of GO terms)
ENRICHED HAPLOINSUFFICIENT SET
Number of genes
RNA polymerase complex
nuclear DNA-directed RNA polymerase complex DNA-directed RNA polymerase II, core
DNA-directed RNA polymerase II, core complex
DNA-directed RNA polymerase ll, holoenzyme
DNA-directed RNA polymerase I complex
DNA-directed RNA polymerase III complex
telomere organization and biogenesis
ribosome biogenesis and assembly
proteasome regulatory particle
transcription regulator activity
ENRICHED HAPLOPROFICIENT SET
Number of genes
nucleobase, nucleoside, nucleotide and nucleic acid transmembrane transporter activity
G1/S transition of mitotic cell cycle
transcription factor TFIID complex
regulation of DNA metabolic process
trehalose metabolic process
Analysis of the smaller set of HP genes from the turbidostat selection showed that the over-represented functional categories included the cell cycle (e.g. spindle checkpoint genes) and macromolecule storage (e.g. genes encoding proteins involved in trehalose metabolism). The fraction of HFC genes in various organelles, selected protein complexes and cellular processes is shown in Additional File 5, Figure S2.
HFC at rapid growth in comparison to nutrient-limited chemostats at low growth rates
First, a yeast culture in synthetic complete medium (FPM) was grown in chemostat mode at D = 0.3 h-1. This means the same growth conditions were applied as in the turbidostat except that the cells were growing at a fixed dilution rate close to the maximum specific growth rate. The dilution rate in the turbidostats, and hence the maximum growth rate of the pool, stabilised at 0.32 h-1. However, 0.3 h-1 was used in the chemostat, because operating at 0.32 h-1 would be at the limit of the stability for a chemostat culture of the pool of mutants. We found that the set of HFC genes detected under these conditions was very similar to the set detected in turbidostat competitions (Additional File 3, Table S1). Thus, similar sets of HFC genes were found at similar growth rates despite the fact that one culture was nutrient-limited (the chemostat with FPM is leucine-limited), while the other was nutrient sufficient (turbidostat).
Next, we carried out competitions in ammonium-limited chemostats at higher dilution rates (D = 0.2 h-1 and 0.3 h-1). We will refer to these two conditions as N02 and N03, respectively, with a nitrogen-limited chemostat at D = 0.1 h-1 being termed N01; correspondingly cultures grown in FPM are termed FPMT, if run in turbidostatic mode, and FPM03, if run as a chemostat at D = 0.3 h-1. First of all, the results demonstrated that the HFC genes detected under nitrogen limitation vary considerably at different dilution rates. The correlation of FCC's (Figure 1B) from N01 and N02 is 0.54, indicating that most mutants behave similarly under the two conditions, however when dilution rate is increased further to 0.3, the similarity disappears. Competitions with different nutrient limitations (C01, N01 and P01) are highly correlated to each other at a slow growth rate, and cluster together on the dendrogram. The fast-growth rate experiments cluster together (FPMT, FPM03 and N03) and are highly correlated to each other. Hence the major difference between the competitions is the rate of growth achieved, rather than the nutrients being in plentiful or limiting supply.
Genetic control of growth rate is qualitatively, but not quantitatively, correlated with, growth rate regulation of gene expression
It should be noted that, although Category 1 experiments showed the transcript levels of genes encoding components of the mitochondrial ribosome to be up-regulated as a function of growth rate in chemostats, the competition experiments in turbidostats did not indicate that any growth-rate control was exerted by these components. Although mitochondrial activity has long been known to be an important factor that affects the growth rate in yeast (e.g. respiratory-deficient petite mutants are slow-growing ), our competition experiments do not suggest that genes encoding mitochondrial proteins exert any significant growth-rate control. Expression of the mitochondrial ribosomal proteins were up-regulated in C-limited chemostats at high dilution rates , indicating high levels of respiration were taking place at higher dilution rates under carbon limitation. The lack of growth-rate control by mitochondria in our experiments could be an outcome of the high glucose concentration in both the turbidostat and the nitrogen-limited chemostat at D = 0.3 h-1, where excess glucose would repress respiration and promote fermentation.
Range of flux control varies among genes within the same functional category
Most protein complexes are essential for viability, though not every component of a given complex is essential. By analogy, it can be hypothesized that most protein complexes exert growth-rate control when cells are growing rapidly, though not every gene that encodes a component of a particular has a significant flux control coefficient (FCC'). We have selected a subset of protein complexes and cellular processes to exemplify the range of FCC' values in turbidostats (Figure 3; Additional File 5, Figure S2). A highly haploinsufficient small protein complex like RNA polymerase I, can include a few HP genes with a small range of FCC's, while the genes encoding the rest of its subunits have a large range of negative growth-rate control coefficients. Moreover, the average FCC' values for members of different multiprotein complexes was found to vary in turbidostat competitions. For instance, genes encoding sub-units of the COPI vesicle coat and the exosome have high average positive FCC' values, while those encoding components of the RNA polymerase I complex and the chaperonin-containing T-complex have high average negative FCC's.
Genes with functions involved in the cell cycle and the maintenance of genome integrity are associated with haploproficiency
Genes involved in the cell cycle that determine the G1/S phase transition and progression through S phase were particularly enriched in our set of haploproficient genes. A subset of cell cycle related genes are expressed periodically as the cell cycle progresses  and regulated post-transcriptionally by a sequence of events  including 'just-in-time assembly' . A reduction in the copy number of some of these genes in the results in growth rate above μmax in turbidostats (Additional File 3, Table S1). This indicates that a reduction of controls that delay the cell cycle under normal conditions can, in fact, allow the cycle proceed more quickly in these HP hemizygotes.
Amongst those mutants that showed a large change in their cycle maps were deletants of CDC28 and two cyclin genes, CLB2 and CLN3, whose products promote the G2/M and G1/S transitions, respectively [16, 17]. Both CLN3 and CDC28 were HP at high growth rates, suggesting a link between altered progress through the cycle and faster growth. We have selected a subset of 9 genes that control the G2/M transition to investigate their effect on progress of cell cycle and phenotype in tetraploid deletion mutants (Alcasabas and co-workers, submitted). Other genes showing both haploproficiency and an altered cell cycle included DMA2, whose human ortholog RNF8 regulates the G2/M DNA damage checkpoint , and EPL1, a component of the NuA4 complex, which has also been linked to DNA repair . We will now discuss in more detail the set of haploproficient genes involved in the maintenance of genome integrity.
Genome integrity related functions are tightly linked to cell cycle as replication of genome requires timely functioning of various protein and protein complexes regulated by cell cycle. 207 of the 384 genes annotated to the GO term 'chromosome organisation' are HFC genes. Many of these genes, as well as others, are directly involved in the maintenance of genome integrity and also show a haploproficient phenotype. These genes and their functions (some 42 of them are listed in Additional File 6, Table S3). It is clear from the results of the turbidostat competition experiments that yeast is capable of growing at faster than maximum specific growth rate of the wild-type strain. However, the imperative of correctly replicating and segregating its chromosomes, repairing DNA damage, and generally maintaining its chromosome integrity has imposed restraints on the absolute rate of growth that can be sustained over evolutionary time.
Growth-Rate Control is a Complex Function of Protein Activity
The relationship between gene copy number and protein levels is an important factor that must be taken into consideration when investigating growth-rate control. Torres and co-workers  have reported that aneuploid yeast cells overexpress the genes on the extra chromosomes (93% of the genes on the extra chromosomes were found to be overexpressed by at least a factor of 1.3-fold compared to the wild-type haploid parent). Recently, it has also been reported that, for most genes tested, the average level of their protein product is proportional to gene copy number, which indicates that dosage compensation is a rare event .
Springer and co-workers  investigated 643 genes and 123 of them were found to have significant levels of dosage compensation when one copy was deleted from the genome, while 18 genes reduced their expression to less than half of the wild-type level, when one copy is deleted from the genome. We reasoned that more HI genes should be found in the set of genes that exhibit no dosage compensation and more HP or non-HFC genes in the set of genes that exhibit some dosage compensation. This is because a lack of dosage compensation is more likely to result in inadequate protein concentrations in the cell, and thus result in a reduction of fitness. Genes with dosage compensation would be expected to have less growth-rate control as the concentrations of their protein products in a hemizygote will be closer to those in the wild type. Furthermore, haploproficiency can be expected if dosage compensation leads to product levels that exceed those of the wild type. We have compared the dosage compensation values  (log2 transformed fold changes, 0 indicates no dosage compensation, 1 indicates dosage compensation that matches the wild-type level of the protein) to FCC' values from turbidostat competitions and the overlap of HFC genes to genes with or without significant dosage compensation.
It was further reported by Springer and co-workers  that dosage compensation remains the same for a protein across different growth conditions. If dosage compensation (or the lack of it) was a simple indication of high growth-rate control, then we would expect to find a similar set of HFC genes in different conditions, and flux control coefficients to be (negatively) correlated to dosage compensation, both of which we failed to find. Hence, it should be realised that the flux control exerted by a protein is a function of its biological activity (which is probably context dependent), rather than simply on its molecular abundance. None of the high-throughput studies (including our own) measures biological activity, and so it might be anticipated that FCC' values are a complex function of protein abundance.
Few other factors were shown to link protein abundance to growth-rate control. For example, Sopko and co-workers  found that ca. 15% of the genes in yeast are toxic when overexpressed; that is, their overexpression reduced the growth rate. This set of toxic genes was enriched in cell cycle genes, particularly in genes that are expressed periodically. This finding supports the hypothesis that a protein can be toxic even at its wild-type levels and thus the expression of its cognate gene has to be carefully controlled. Such toxicity may increase as the protein abundance increases. It follows that a reduction in the concentration of a protein that acts as an inhibitor of growth at its wild-type concentration can have a positive impact on growth rate. Yoshikawa and co-workers  reported that deletion and overexpression mutants of more than 400 genes have growth defects. (Note that very few deletion or over-expression mutants that have a haploproficient phenotype were found -probably because of the low sensitivity of measurements in batch culture in microtitre plates). Another factor that should be considered is dosage imbalance as a result of deviation from wild-type protein levels ([24–27]).
In order to explain our results with a dosage v. flux curve, we have to extend its applicability to all the cases we observe and we should expect to find a different curve for each protein. For example, the lack of correlation between dosage compensation and high flux control demonstrates that the threshold lies at a different level of activity for each gene or protein. Further, haploproficiency cannot be explained by a curve that converges to a saturation level, as the maximum growth rate of a heterozygous mutant can be no higher than the maximum growth rate of the wild type on such a curve. Therefore, the flux and the control exerted on it by a protein have to be the outcome of multiple factors, each of which can be hypothesized to be a non-linear function of gene (protein) dosage/activity. The activity v. flux curve shown in Figure 6 has to be coupled with additional functions that explain the contributions or limitations imposed on the growth rate by the protein - either due to its toxic effects and functional inhibition on growth, or to the imbalances that result from changes in its abundance.
The possibility of non-zero flux, even if the protein concentration is zero, and the possibility of achieving higher fluxes as the protein concentration is increased to higher levels than in the wild type. The former possibility applies to non-essential genes in a diploid organism, and the latter to proteins that increase the flux when they are over-expressed.
Functional inhibition: The function of many proteins is the inhibition or delay of processes that lead to cell growth. The timing and level of expression of genes encoding such proteins is usually fine-tuned by the cells to avoid unwanted inhibition of growth. Our results indicate that cell cycle checkpoint proteins may exert a significant inhibition on growth rate.
Hypothetical examples of such curves are given in Figures 7B-D. It is likely that genes with curves as shown in Figures 7C and 7D constitute a considerable fraction of the genome, rather than being exceptional cases. Revealing the superimposed fitness curve for each gene requires data on growth rates at various gene copy numbers (or protein concentrations), examples of which are starting to appear in the literature (Alcasabas and co-workers (submitted); ).
Haploinsufficiency is the situation where the reduction in the copy number of a gene in a diploid organism from two to one results in a significant loss of fitness. It is assumed that the reduction in copy number results in a corresponding reduction in the concentration of the gene's protein product, such that there is insufficient of that protein present to sustain the wild-type phenotype. The connection of haploinsufficiency to diseases has long been known [29–31] and various hypotheses have been proposed to explain the mechanisms of haploinsufficiency ([26, 25]). The kinetics of gene expression in connection to haploinsufficiency has also been studied . Haploinsufficiency was initially identified in gene-by-gene phenotype screenings in humans and other organisms. However, the availability of a genome-wide collection of deletion mutants, as well as tag arrays for quantifying the proportions of individual mutants in a population, has enabled genome-wide screening for phenotypes of S. cerevisiae deletion mutants under various conditions ([5, 33, 4]). Previously, sets of HFC genes have been reported for the following cell populations grown in batch culture: (i) homozygous deletion mutant pool growing in YPD ; (ii) both homozygous and heterozygous pools growing on fermentable and non-fermentable carbon sources ; and (iii) heterozygous and homozygous pools growing in YPD and minimal medium . Some of these pools have also been screened for haploinsufficiency when treated with drugs and sets of genes with drug sensitivity and resistance have been identified ([36–38]).
Haploproficiency under fast growth rate conditions has not been reported until recently  in genome-wide phenotype screenings of S. cerevisiae. However, it has been reported that it is possible to select for mutants with 5-50% increase in growth rate (as compared to wild type) under specific nutrient limitations . Some of these evolved strains had transporter gene amplifications, together with other mutations not associated to related physiology. It was later reported that Kluyveromyces marxianus, allowed to evolve in a mineral medium, can attain higher growth rates than the mother strain. Changes in the morphology of these cells were found to be linked to the increased maximum growth rate, probably because a larger surface area (relative to cell volume) of the elongated mutants allows transport processes to take place more efficiently. See the review by Cakar and co-workers for recent approaches in evolutionary engineering of strains .
Although data on the evolution of strains have provided a valuable resource for elucidating the mechanisms of growth rate control, systematic screens of single-gene deletion mutants provide genome-wide profiling of all known genes. While mutants with severe growth defects can be easily detected in such screenings, it has been difficult to detect single-gene deletion mutants which grow faster than the wild type, as the difference in growth rate can be small and below the threshold of detection. This hindered the identification of haploproficient hemizygotes in many studies. We have shown that competition at slow growth rates in nutrient-limited chemostats can be used to identify growth advantages as well as growth deficiencies at high sensitivity . Most of the genes only displayed haploproficiency under a particular nutrient limitation; for example, most genes encoding subunits of the 26S proteasome were HP under nitrogen limitation at D = 0.1 h-1. The reason for their haploproficiency under these conditions is unknown; however, it is likely to be a positive outcome of reducing the rate of protein turnover when there is competition for nitrogen resources.
In a recent study on Schizosaccharomyces pombe, 136 heterozygous deletion mutants in competition were shown to have faster growth rates than wild type when grown in batch on a rich medium. The same study also reported 183 haploproficient S. cerevisiae genes derived from a re-analysis of previously published data from a competition of heterozygous deletion pools under similar conditions . Screening of a library of hypomorphic alleles of essential genes also revealed a small set of faster-growing mutants with a reduced concentration of essential proteins . In this study, we have identified a large set of haploproficient as well as haploinsufficient hemizygotes using S. cerevisiae cultures growing at their maximum growth rates. Although these results are novel, it is not surprising that reducing the level of proteins that act as negative regulators of processes required for rapid growth, or that are toxic as monomers in the cytosol, can have a positive impact on growth rate. Although we might expect such growth-rate-limiting factors to be eliminated by evolution in order to attain yet higher growth rates and thus greater fitness, it is possible that these factors were optimized rather than eliminated, without compromising the basic function of the protein, which may be essential under sub-optimal conditions.
Our HI gene set includes 124 of the 186 HI genes reported by Deutschbauer and co-workers  (GR < 0.97) (our set had data on 166 of the 186), while 8 of the 166 are significantly HP in our set. Our HP gene set includes only 21 of the 181 HP genes reported by Deutschbauer and co-workers  (GR > 1.03), 178 of which exist in our data set. 46 of these genes are HI in our set, demonstrating that the results from the two experiments vary significantly with respect to the identification of the HP phenotype. There can be a number of reasons for the discrepancy between the two datasets. First, our classification criteria are based on goodness-of-fit rather than the percent change in growth rate - so mutants with small, but significant, changes are included in our HFC set, while mutants with large changes may have been excluded if their readings were noisy. Moreover, the experimental regimes used in the two studies are different since continuous turbidostats were used in this study as opposed to batch shake flasks. Finally, both experiments are expected to have some level of stochasticity in both the competition and sampling steps, and downstream processing of the samples may also introduce noise into the final results.
Features of genes and proteins are not globally correlated to flux control
Haploinsufficiency can be linked to a number of factors, and previous authors have found evidence for the HI phenotype of hemizygous deletants being the result of either the insufficient concentration of the protein encoded by the single-copy gene or of dosage imbalance . Noise in gene expression has also been proposed as a reason for haploinsufficiency ([32, 44]), while Veitia  suggested synergistic interactions between transcription factors (TF) as the cause of the haploinsufficiency of TF mutants. Accordingly, we compared our FCC' measurements to data from the literature to determine whether these, or any other factors, had an important impact on FCC' values in rapidly growing yeast.
Haploinsufficiency is often described as a direct consequence of inadequate concentrations of proteins , so it can be hypothesized that the proteins with higher abundance will have higher flux control when their gene copy number is reduced - particularly if dosage compensation in hemizygous mutants is rare, as reported by Springer and co-workers . However, we found no significant correlation between protein abundance  and high flux control at rapid growth. It was also hypothesized that proteins with noisy expression are more likely to fall below the threshold  that will cause irreversible growth defects. We have compared modified flux control coefficients to noise in mRNA expression (DV values from ). There was no genome-wide correlation between high growth-rate control and noise in mRNA expression, nor were the rates of mRNA and protein decay or half-lives [48–50] correlated to growth-rate control. In a recent study on haploinsufficiency in humans , genetic interaction network features were identified as the best predictors of haploinsufficiency. However, the correlation between the number of genetic interactions of genes and their FCC' value in our turbidostat study were insignificant (rho = -0.049 for positive interactions, -0.055 for negative interactions, -0.015 for synthetic lethal, rho = -0.065 for all physical and genetic interactions of yeast from BioGrid, ). We have also demonstrated that, although growth-regulated genes (from our Category 1 experiments) are, themselves, controllers of flux at the fastest growth rates, there is no simple relationship between the change in growth rate due to haploinsufficiency and the fold change in transcript level or protein concentration.
Haploinsufficiency has been linked to the stoichiometry of protein complexes as there is evidence that dosage imbalance has a profound effect on phenotype in some cases ([25, 29]). Previous competition experiments on homozygous and heterozygous deletion mutant pools grown in batch on YPD have also shown that genes encoding the components of protein complexes are enriched among the HI genes . Our HFC genes are also enriched in genes encoding components of protein complexes (hypergeometric distribution, p-value = 1.5e-7 for HI genes, p-value = 0.034 for HP genes; the genes assigned to GO Term 'protein complex' were used in the analysis (Saccharomyces Genome Database, SGD, http://www.yeastgenome.org). Moreover, HI genes and their protein products are involved in more genetic or protein-protein interactions than expected by chance (p-value = 8.7e-7 for HI as opposed to p-value = 0.22 for HP, data from Biogrid ). Both HI and HP genes are more likely to be essential than the average for the genome (hypergeometric distribution, p-value = 0.0051 for HI genes and p-value = 0.039 for HP genes; list of essential genes downloaded from SGD). These results demonstrate that protein complexes, or proteins that function in interaction with others, are more likely to be involved in high flux control when compared to proteins that act as singletons in the cell.
Our method has identified many more HI and HP genes than have been reported previously. Our results cover 124 of the 166 HI genes previously identified in YPD . The relatively poor overlap with previous data on HP genes indicates systematic differences between experimental design and data analysis regimes in the two different approaches.
We find that translation and transcription are the bottlenecks for fast growth when the nutrients are not limited. The enrichment of cell cycle genes in the set of genes that are toxic upon over-expression supports our hypothesis that reduction in protein levels increases growth rate as wild-type levels limit the growth rate by 'functional inhibition'. Accumulation of cells in the G1 phase may be an indication of its high flux control, although no direct link exists between the two. Genes involved in the maintenance of genome integrity also show a haploproficient phenotype under conditions of rapid growth. Thus the checks and balances that the yeast cell imposes on itself to ensure that its genome integrity is maintained and that the genome is faithfully replicated and segregated during the cell cycle limit the absolute rate of growth which is compatible with the organism's long-term survival.
Proteins may have different levels of flux control under different growth conditions; their flux control can even be in opposite directions depending on the nature of the competition. Hence, high flux control is context-dependent and is not a direct consequence of any single factor, such as protein abundance or growth-rate regulation of gene expression. It is more likely to be a complex function of protein activity, which can be hypothesised as a superimposition of functions set by multiple factors.
A 1 ml aliquot of the heterozygous deletion mutant (BY4743, MAT a /MATα his3Δ1/his3 Δ1 leu2Δ0/leu2Δ0 lys2Δ0/LYS2 MET15/met15Δ 0 ura3Δ 0/ura3 Δ0) pool (obtained from the Yeast Genome Deletion Project library and prepared as described in Additional File 1, Additional Methods) were grown in 100 ml of YPD media overnight at 30°C with shaking at 180 rpm. 10 ml of this preculture was used to inoculate 1L of the appropriate nitrogen-limited or FPM growth media in 2L fermentors. Cultures were grown overnight at 30°C with aeration and stirring at 750 rpm. After 24 hours the fermentors were switched to continuous growth with pH control (set point = 4.5) and either a fixed dilution rate (chemostat) or dilution rate determined by biomass (turbidostats). In turbidostat mode , the biomass corresponding to approximately half of the maximum reading from batch phase (corresponding to mid-exponential growth). 20 ml samples were taken from the preculture and overnight batch culture followed by at least four at approximately 24-hour intervals from the continuous cultures. Experiments were made with at least two biological replicates for each fermentation and two technical replicates for each sample. For further details on fermentations including instrumentation and sampling, see the Additional File 1, Additional Methods section.
DNA extraction, tag amplification, hybridization
Total DNA was isolated from each of the samples using a Promega Wizard Genomic DNA purification kit, according to the protocol provided together with the kit. Universal primers were then used to amplify the up-tag and down-tag sequences with separate PCR reactions as described previously . The PCR products were purified, quantified and hybridized on custom Affymetrix Tag3/4 arrays as described previously ([4, 54]).
Tag3 data were normalised and modelled as described previously . Tag4 data were normalized as described  and modelled as described previously . In both cases a relative growth rate and multiple testing corrected P-value (null hypothesis that relative growth rate is equal to 0) were obtained for each strain. The 182 deletion mutants listed in  were excluded from further analysis. The relative growth rates were used to calculate the Spearman correlation coefficients between experiments shown in Figure 1B. In Figure 1B the dendrogram showing the relationship between different experiments was constructed using the Euclidean distance between rows (experiments) and complete linkage hierarchical clustering. The functional analysis was made using Gene Ontology  as described in the Additional File 1, Additional Methods. All statistical methods were implemented either in R/Bioconductor  using RSRuby [http://rubyforge.org/projects/rsruby/] or Matlab [The Mathworks, Inc]. For further details on sample processing and data analysis, see the Additional File 1, Additional Methods section.
Derivation of modified flux control coefficients
where A is the scaling factor (initial copy number/D).
Note that the FCC' value of a gene depends on the local slope of a protein concentration vs fitness curve at the range of interest (Figures 6 and 7 give examples of such curves). For example, the FCC' of an essential gene can be very high in the range of 0 copies to 1 copy (change in flux in heterozygous deletion mutant with respect to the inviable null mutant), or the FCC' of a gene whose protein product is at a concentration higher than needed in the wild-type cell (i.e. a level beyond saturation) would be close to zero at proximity of the wild-type level. In this study, the relative growth rate (or FCC') gives the negative average slope of the protein concentration vs fitness curve between 1 copy and 2 copies of each gene, i.e. our HI genes with negative FCC's indicate a positive flux control by the gene when its copy number is increased from 1 to 2, and similarly our HP genes indicate a negative flux control. These principles apply both for gene copy number and for protein abundance, based on the assumption that diploid cells with only one copy of a gene always have a lower concentration of its protein product than the cells with two copies.
Cell cycle analysis
Cell cycle profiles were determined using the method described in Haase and co-workers . Briefly, cultures were grown to exponential phase before fixation in 70% v/v ethanol, followed by RNase A and protease treatment. Washed samples were then mixed with 1 μM Sytox Green and analysed for DNA content in a Beckman Coulter CyanADP flow cytometer. 20,000 events were captured for each sample with three biological replicates (separate cultures) for each strain. Two peaks were fit to the count histogram of DNA content using the Dean-Jett-Fox method  and the percentage of the total cells under the G1 peak was recorded for each sample.
We are grateful to Rachael Walker for assistance with the flow cytometry measurements, and thank Balázs Papp and Daniela Delneri for their comments on the manuscript. This work was supported by BBSRC grant BB/C505140/2 and a contract from the European Commission under the FP7 Collaborative Programme, UNICELLSYS (both to SGO).
- Kacser H, Burns JA: The molecular basis of dominance. Genetics 1981, 97: 639-666.Google Scholar
- Castrillo JI, Zeef LA, Hoyle DC, Zhang N, Hayes A, Gardner DC, Cornell MJ, Petty J, Hakes L, Wardleworth L, et al: Growth control of the eukaryote cell: a systems biology study in yeast. J Biol 2007, 6: 4. 10.1186/jbiol54View ArticleGoogle Scholar
- Gutteridge A, Pir P, Castrillo JI, Charles PD, Lilley KS, Oliver SG: Nutrient control of eukaryote cell growth: a systems biology study in yeast. BMC Biol 2010, 8: 68. 10.1186/1741-7007-8-68View ArticleGoogle Scholar
- Delneri D, Hoyle DC, Gkargkas K, Cross EJ, Rash B, Zeef L, Leong HS, Davey HM, Hayes A, Kell DB, et al: Identification and characterization of high-flux-control genes of yeast through competition analyses in continuous cultures. Nat Genet 2008, 40: 113-117. 10.1038/ng.2007.49View ArticleGoogle Scholar
- Giaever G, Chu AM, Ni L, Connelly C, Riles L, Veronneau S, Dow S, Lucau-Danila A, Anderson K, Andre B, et al: Functional profiling of the Saccharomyces cerevisiae genome. Nature 2002, 418: 387-391. 10.1038/nature00935View ArticleGoogle Scholar
- Allen J, Davey HM, Broadhurst D, Heald JK, Rowland JJ, Oliver SG, Kell DB: High-throughput classification of yeast mutants for functional genomics using metabolic footprinting. Nat Biotechnol 2003, 21: 692-696. 10.1038/nbt823View ArticleGoogle Scholar
- Pirt SJ: Principles of microbe and cell cultivation. Oxford: Blackwell Scientific; 1975.Google Scholar
- Bull AT: The renaissance of continuous culture in the post-genomics age. J Ind Microbiol Biotechnol 2010, 37: 993-1021. 10.1007/s10295-010-0816-4View ArticleGoogle Scholar
- Harris MA, Clark J, Ireland A, Lomax J, Ashburner M, Foulger R, Eilbeck K, Lewis S, Marshall B, Mungall C, et al: The Gene Ontology (GO) database and informatics resource. Nucleic Acids Res 2004, 32: D258-261. 10.1093/nar/gkh036View ArticleGoogle Scholar
- Sartor MA, Leikauf GD, Medvedovic M: LRpath: a logistic regression approach for identifying enriched biological groups in gene expression data. Bioinformatics 2009, 25: 211-217. 10.1093/bioinformatics/btn592View ArticleGoogle Scholar
- Ogur M, St John R: A differential and diagnostic plating method for population studies of respiration deficiency in yeast. J Bacteriol 1956, 72: 500-504.Google Scholar
- Spellman PT, Sherlock G, Zhang MQ, Iyer VR, Anders K, Eisen MB, Brown PO, Botstein D, Futcher B: Comprehensive identification of cell cycle-regulated genes of the yeast Saccharomyces cerevisiae by microarray hybridization. Mol Biol Cell 1998, 9: 3273-3297.View ArticleGoogle Scholar
- Chen KC, Calzone L, Csikasz-Nagy A, Cross FR, Novak B, Tyson JJ: Integrative analysis of cell cycle control in budding yeast. Mol Biol Cell 2004, 15: 3841-3862. 10.1091/mbc.E03-11-0794View ArticleGoogle Scholar
- de Lichtenberg U, Jensen LJ, Fausboll A, Jensen TS, Bork P, Brunak S: Comparison of computational methods for the identification of cell cycle-regulated genes. Bioinformatics 2005, 21: 1164-1171. 10.1093/bioinformatics/bti093View ArticleGoogle Scholar
- Brauer MJ, Huttenhower C, Airoldi EM, Rosenstein R, Matese JC, Gresham D, Boer VM, Troyanskaya OG, Botstein D: Coordination of growth rate, cell cycle, stress response, and metabolic activity in yeast. Mol Biol Cell 2008, 19: 352-367. 10.1091/mbc.E07-08-0779View ArticleGoogle Scholar
- Jorgensen P, Tyers M: How cells coordinate growth and division. Curr Biol 2004, 14: R1014-1027. 10.1016/j.cub.2004.11.027View ArticleGoogle Scholar
- Csikasz-Nagy A, Battogtokh D, Chen KC, Novak B, Tyson JJ: Analysis of a generic model of eukaryotic cell-cycle regulation. Biophys J 2006, 90: 4361-4379. 10.1529/biophysj.106.081240View ArticleGoogle Scholar
- Kolas NK, Chapman JR, Nakada S, Ylanko J, Chahwan R, Sweeney FD, Panier S, Mendez M, Wildenhain J, Thomson TM, et al: Orchestration of the DNA-damage response by the RNF8 ubiquitin ligase. Science 2007, 318: 1637-1640. 10.1126/science.1150034View ArticleGoogle Scholar
- Doyon Y, Selleck W, Lane WS, Tan S, Cote J: Structural and functional conservation of the NuA4 histone acetyltransferase complex from yeast to humans. Mol Cell Biol 2004, 24: 1884-1896. 10.1128/MCB.24.5.1884-1896.2004View ArticleGoogle Scholar
- Torres EM, Sokolsky T, Tucker CM, Chan LY, Boselli M, Dunham MJ, Amon A: Effects of aneuploidy on cellular physiology and cell division in haploid yeast. Science 2007, 317: 916-924. 10.1126/science.1142210View ArticleGoogle Scholar
- Springer M, Weissman JS, Kirschner MW: A general lack of compensation for gene dosage in yeast. Mol Syst Biol 6: 368. 10.1038/msb.2010.19Google Scholar
- Sopko R, Huang D, Preston N, Chua G, Papp B, Kafadar K, Snyder M, Oliver SG, Cyert M, Hughes TR, et al: Mapping pathways and phenotypes by systematic gene overexpression. Mol Cell 2006, 21: 319-330. 10.1016/j.molcel.2005.12.011View ArticleGoogle Scholar
- Yoshikawa K, Tanaka T, Ida Y, Furusawa C, Hirasawa T, Shimizu H: Comprehensive phenotypic analysis of single-gene deletion and overexpression strains of Saccharomyces cerevisiae. Yeast 2011, 28: 349-361. 10.1002/yea.1843View ArticleGoogle Scholar
- Katz W, Weinstein B, Solomon F: Regulation of tubulin levels and microtubule assembly in Saccharomyces cerevisiae: consequences of altered tubulin gene copy number. Mol Cell Biol 1990, 10: 5286-5294.View ArticleGoogle Scholar
- Papp B, Pal C, Hurst LD: Dosage sensitivity and the evolution of gene families in yeast. Nature 2003, 424: 194-197. 10.1038/nature01771View ArticleGoogle Scholar
- Veitia RA: Exploring the etiology of haploinsufficiency. Bioessays 2002, 24: 175-184. 10.1002/bies.10023View ArticleGoogle Scholar
- Kaizu K, Moriya H, Kitano H: Fragilities caused by dosage imbalance in regulation of the budding yeast cell cycle. PLoS Genet 2010, 6: e1000919. 10.1371/journal.pgen.1000919View ArticleGoogle Scholar
- Suk K, Choi J, Suzuki Y, Ozturk SB, Mellor JC, Wong KH, MacKay JL, Gregory RI, Roth FP: Reconstitution of human RNA interference in budding yeast. Nucleic Acids Res 2011, 39: e43. 10.1093/nar/gkq1321View ArticleGoogle Scholar
- Veitia RA: A generalized model of gene dosage and dominant negative effects in macromolecular complexes. FASEBJ 2010, 24: 994-1002. 10.1096/fj.09-146969View ArticleGoogle Scholar
- Seidman JG, Seidman C: Transcription factor haploinsufficiency: when half a loaf is not enough. J Clin Invest 2002, 109: 451-455.View ArticleGoogle Scholar
- Bosl WJ, Li R: The role of noise and positive feedback in the onset of autosomal dominant diseases. BMC Syst Biol 2010, 4: 93. 10.1186/1752-0509-4-93View ArticleGoogle Scholar
- Cook DL, Gerber AN, Tapscott SJ: Modeling stochastic gene expression: implications for haploinsufficiency. Proc Natl Acad Sci USA 1998, 95: 15641-15646. 10.1073/pnas.95.26.15641View ArticleGoogle Scholar
- Scherens B, Goffeau A: The uses of genome-wide yeast mutant collections. Genome Biol 2004, 5: 229. 10.1186/gb-2004-5-7-229View ArticleGoogle Scholar
- Steinmetz LM, Scharfe C, Deutschbauer AM, Mokranjac D, Herman ZS, Jones T, Chu AM, Giaever G, Prokisch H, Oefner PJ, Davis RW: Systematic screen for human disease genes in yeast. Nat Genet 2002, 31: 400-404.Google Scholar
- Deutschbauer AM, Jaramillo DF, Proctor M, Kumm J, Hillenmeyer ME, Davis RW, Nislow C, Giaever G: Mechanisms of haploinsufficiency revealed by genome-wide profiling in yeast. Genetics 2005, 169: 1915-1925. 10.1534/genetics.104.036871View ArticleGoogle Scholar
- Giaever G, Flaherty P, Kumm J, Proctor M, Nislow C, Jaramillo DF, Chu AM, Jordan MI, Arkin AP, Davis RW: Chemogenomic profiling: identifying the functional interactions of small molecules in yeast. Proc Natl Acad Sci USA 2004, 101: 793-798. 10.1073/pnas.0307490100View ArticleGoogle Scholar
- Hillenmeyer ME, Ericson E, Davis RW, Nislow C, Koller D, Giaever G: Systematic analysis of genome-wide fitness data in yeast reveals novel gene function and drug action. Genome Biol 2010, 11: R30. 10.1186/gb-2010-11-3-r30View ArticleGoogle Scholar
- Lanthaler K, Bilsland E, Dobson PD, Moss HJ, Pir P, Kell DB, Oliver SG: Genome-wide assessment of the carriers involved in the cellular uptake of drugs: a model system in yeast. BMC Biol 2011, 9: 70. 10.1186/1741-7007-9-70View ArticleGoogle Scholar
- Kim DU, Hayles J, Kim D, Wood V, Park HO, Won M, Yoo HS, Duhig T, Nam M, Palmer G, et al: Analysis of a genome-wide set of gene deletions in the fission yeast Schizosaccharomyces pombe. Nat Biotechnol 2010, 28: 617-623. 10.1038/nbt.1628View ArticleGoogle Scholar
- Gresham D, Desai MM, Tucker CM, Jenq HT, Pai DA, Ward A, DeSevo CG, Botstein D, Dunham MJ: The repertoire and dynamics of evolutionary adaptations to controlled nutrient-limited environments in yeast. PLoS Genet 2008, 4: e1000303. 10.1371/journal.pgen.1000303View ArticleGoogle Scholar
- Groeneveld P, Stouthamer AH, Westerhoff HV: Super life--how and why 'cell selection' leads to the fastest-growing eukaryote. FEBSJ 2009, 276: 254-270. 10.1111/j.1742-4658.2008.06778.xView ArticleGoogle Scholar
- Cakar ZP, Turanli-Yildiz B, Alkim C, Yilmaz U, Nielsen J: Evolutionary engineering of Saccharomyces cerevisiae for improved industrially important properties. FEMS Yeast Res 2011.Google Scholar
- Breslow DK, Cameron DM, Collins SR, Schuldiner M, Stewart-Ornstein J, Newman HW, Braun S, Madhani HD, Krogan NJ, Weissman JS: A comprehensive strategy enabling high-resolution functional analysis of the yeast genome. Nat Methods 2008, 5: 711-718. 10.1038/nmeth.1234View ArticleGoogle Scholar
- Batada NN, Hurst LD: Evolution of chromosome organization driven by selection for reduced gene expression noise. Nat Genet 2007, 39: 945-949. 10.1038/ng2071View ArticleGoogle Scholar
- Veitia RA: Exploring the molecular etiology of dominant-negative mutations. Plant Cell 2007, 19: 3843-3851. 10.1105/tpc.107.055053View ArticleGoogle Scholar
- Ghaemmaghami S, Huh WK, Bower K, Howson RW, Belle A, Dephoure N, O'Shea EK, Weissman JS: Global analysis of protein expression in yeast. Nature 2003, 425: 737-741. 10.1038/nature02046View ArticleGoogle Scholar
- Newman JR, Ghaemmaghami S, Ihmels J, Breslow DK, Noble M, DeRisi JL, Weissman JS: Single-cell proteomic analysis of S. cerevisiae reveals the architecture of biological noise. Nature 2006, 441: 840-846. 10.1038/nature04785View ArticleGoogle Scholar
- Belle A, Tanay A, Bitincka L, Shamir R, O'Shea EK: Quantification of protein half-lives in the budding yeast proteome. Proc Natl Acad Sci USA 2006, 103: 13004-13009. 10.1073/pnas.0605420103View ArticleGoogle Scholar
- Wang Y, Liu CL, Storey JD, Tibshirani RJ, Herschlag D, Brown PO: Precision and functional specificity in mRNA decay. Proc Natl Acad Sci USA 2002, 99: 5860-5865. 10.1073/pnas.092538799View ArticleGoogle Scholar
- Miller C, Schwalb B, Maier K, Schulz D, Dumcke S, Zacher B, Mayer A, Sydow J, Marcinowski L, Dolken L, et al: Dynamic transcriptome analysis measures rates of mRNA synthesis and decay in yeast. Mol Syst Biol 2011, 7: 458.View ArticleGoogle Scholar
- Huang N, Lee I, Marcotte EM, Hurles ME: Characterising and predicting haploinsufficiency in the human genome. PLoS Genet 2010, 6: e1001154. 10.1371/journal.pgen.1001154View ArticleGoogle Scholar
- Stark C, Breitkreutz BJ, Chatr-Aryamontri A, Boucher L, Oughtred R, Livstone MS, Nixon J, Van Auken K, Wang X, Shi X, et al: The BioGRID Interaction Database: 2011 update. Nucleic Acids Res 2011, 39: D698-704. 10.1093/nar/gkq1116View ArticleGoogle Scholar
- Davey HM, Davey CL, Woodward AM, Edmonds AN, Lee AW, Kell DB: Oscillatory, stochastic and chaotic growth rate fluctuations in permittistatically controlled yeast cultures. Biosystems 1996, 39: 43-61. 10.1016/0303-2647(95)01577-9View ArticleGoogle Scholar
- Pierce SE, Davis RW, Nislow C, Giaever G: Genome-wide analysis of barcoded Saccharomyces cerevisiae gene-deletion mutants in pooled cultures. Nat Protocols 2007, 2: 2958-2974. 10.1038/nprot.2007.427View ArticleGoogle Scholar
- Eason RG, Pourmand N, Tongprasit W, Herman ZS, Anthony K, Jejelowo O, Davis RW, Stolc V: Characterization of synthetic DNA bar codes in Saccharomyces cerevisiae gene-deletion strains. Proc Natl Acad Sci USA 2004, 101: 11046-11051. 10.1073/pnas.0403672101View ArticleGoogle Scholar
- The Gene Ontology project in 2008 Nucleic Acids Res 2008, 36: D440-444.Google Scholar
- Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, Ellis B, Gautier L, Ge Y, Gentry J, et al: Bioconductor: open software development for computational biology and bioinformatics. Genome Biol 2004, 5: R80. 10.1186/gb-2004-5-10-r80View ArticleGoogle Scholar
- Small, JR, Kacser H: Responses of metabolic systems to large changes in enzyme activities and effectors 2. The linear treatment of branched pathways and metabolite concentrations - assessment of the general non-linear case. Eur J Biochem 1993, 213: 625-640. 10.1111/j.1432-1033.1993.tb17802.xView ArticleGoogle Scholar
- Haase SB, Reed SI: Improved flow cytometric analysis of the budding yeast cell cycle. Cell Cycle 2002, 1: 132-136.View ArticleGoogle Scholar
- Fox MH: A model for the computer analysis of synchronous DNA distributions obtained by flow cytometry. Cytometry 1980, 1: 71-77. 10.1002/cyto.990010114View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.