Development of a Fur-dependent and tightly regulated expression system in Escherichia colifor toxic protein synthesis
© Guan et al.; licensee BioMed Central Ltd. 2013
Received: 8 September 2012
Accepted: 8 March 2013
Published: 19 March 2013
There is a continuous demanding for tightly regulated prokaryotic expression systems, which allow functional synthesis of toxic proteins in Escherichia coli for bioscience or biotechnology application. However, most of the current promoter options either are tightly repressed only with low protein production levels, or produce substantial protein but lacking of the necessary repression to avoid mutations initiated by leaky expression in the absence of inducer. The aim of this study was to develop a tightly regulated, relatively high-efficient expression vector in E. coli based on the principle of iron uptake system.
By using GFP as reporter, PfhuA with the highest relative fluorescence units, but leaky expression, was screened from 23 iron-regulated promoter candidates. PfhuA was repressed by ferric uptake regulator (Fur)-Fe2+ complex binding to Fur box locating at the promoter sequence. Otherwise, PfhuA was activated without Fur-Fe2+ binding in the absence of iron. In order to improve the tightness of PfhuA regulation for toxic gene expression, Fur box in promoter sequence and fur expression were refined through five different approaches. Eventually, through substituting E. coli consensus Fur box for original one of PfhuA, the induction ratio of modified PfhuA (named PfhuA1) was improved from 3 to 101. Under the control of PfhuA1, strong toxic gene E was successfully expressed in high, middle, low copy-number vectors, and other two toxic proteins, Gef and MazF were functionally synthesized without E. coli death before induction.
The features of easy control, tight regulation and relatively high efficiency were combined in the newly engineered PfhuA1. Under this promoter, the toxic genes E, gef and mazF were functionally expressed in E. coli induced by iron chelator in a tightly controllable way. This study provides a tightly regulated expression system that might enable the stable cloning, and functional synthesis of toxic proteins for their function study, bacterial programmed cell death in biological containment system and bacterial vector vaccine development.
With the advent of the post-genomic era coming, the need is boosting to express a growing number of genes originating from different organisms . Unfortunately, many of these foreign genes severely interfere with the survival of Escherichia coli cells, which could lead to bacteria death or cause significant defects in bacteria growth. What’s more, following the development of biological technology, many genetically engineered E. coli have been constructed and developed for different purposes, such as bioremediation, biomedicine and bioenergy . However, their practical applications in the field are still restricted because engineered bacteria may cause new environmental contaminations. To minimize the potential risks, biological containment system was designed to monitor and restrict the distribution of engineered bacteria. Generally, containment systems are based on a ‘killer gene’ and a tight ‘regulatory circuit’ that controls expression of the killer gene in response to the presence or absence of environmental signals [3, 4]. In order to activate the killer gene expression only at expected condition, usually a specific environment, the promoter that controls bacteria programmed cell death needs to be easily controllable, tightly regulated and environment-responsive.
In general, there are five solutions to manage killer or highly toxic genes expression: manipulation of transcriptional and translational control elements of highly toxic genes, manipulation of the coding sequence, manipulation of the copy number, addition of stabilizing sequences, and empirical selection of E. coli strains . Among these solutions, manipulation of transcriptional control elements aims at blocking the leaky expression from the common used promoters, such as Plac, Ptrc, Ptac, PT7, PpL, PtetA or PlacUV5, which is critical in toxic gene expression [6, 7]. As a result, the manipulation strategies of transcriptional control elements have been described to reduce the possibility of a unexpected toxic event. The strategies include suppression of basal expression from leaky inducible promoters, suppression of read-through transcription from cryptic promoters, tight control of plasmid copy numbers and proteins production as inactive (but reversible) forms [5, 8–10].
Since the incomplete repression of promoter presents a major problem when cloning genes that encode lethal product to the bacterial host , it is obviously critical to tighten the control of gene expression by utilizing a tightly controllable promoter that permits E. coli normal growth until the very moment of highly toxic gene induction. Starting from this point, several tightly regulated expression circuits have been developed, such as Prha, hybrid Plac/ara-1, and PBAD-based vectors . The widely used PBAD derivative expression systems have been verified as useful solutions for their stringentness in toxic protein production in E. coli. Expression is induced to high levels on media containing L-arabinose, and tightly shuts off on media with glucose but without L-arabinose, which shows more stringent regulation of target gene expression than other expression systems. However, some problems still exist in this system, such as few vectors available and catabolism repression by glucose .
Usually, environment-responsive promoters are commonly engineered in inducible expression system, such as promoters sensing pH, temperature, oxygen concentration or iron availability . From iron-uptake systems that are evolved by bacteria growing in iron-limiting environment, many iron-related promoters have been reported [13–15]. The promoter from iron-uptake regulon is strongly repressed in iron rich conditions by Fur, but fully derepressed in absence of iron . Usually, Fur protein complexing with ferrous irons binds with high affinity to the 19-bp inverted repeat consensus sequence known as the Fur box (GATAATGAT [A/T] ATCATTATC) in the relevant promoter area, which controls transcription of iron-responsive genes in microorganisms . Fur inhibits transcription initiation by blocking the entry of RNA polymerase (RNAP) to the promoter. This ability is tightly dependent on the relative affinities of RNAP and Fur to binding sites in the DNA .
In this study, a novel tightly inducible expression system was developed as an alternative of the current ones. This system, designated pYPfhuA1, was capable of extremely tight regulation and allowed cloning of genes encoding highly toxic products within various copy-number plasmids. In detail, using E. coli Top10 as bacterial host, the strong iron-regulated promoter PfhuA was selected as the primary candidate. By modifying Fur box in promoter sequences and Fur repressor synthesis, the tightness of PfhuA was decreased to varying degrees. Thereinto, PfhuA1, which could be strict repressed by excess iron and efficiently induced by iron chelators, was the ideal promoter candidate for toxic gene expression in E. coli host.
Preliminary screening for iron-regulated promoters
To evaluate the regulation performance of PfhuA, the GFP synthesis was detected when Top10/ptPfhuAG growing in LB medium supplemented with repressor (40 μM FeSO4) or inducer (200 μM 2,2’-dipyridyl). As shown in Figure 1B, even with the repressor addition, an obvious leaky expression was detected under PfhuA transcription. Based on the data, the induction ratio that showed the tightness of promoter was calculated as 3.4, indicating that the tightness of PfhuA had to be improved for the toxic gene expression.
Modifications of PfhuAto improve its tightness
Comparison of different promoter modified strategies
Modified strategies and descriptions
Relative induction fold
E. coli Top10
200 ± 10
200 ± 9
260 ± 11
7812 ± 31
Original fhuA promoter
3292 ± 20
10116 ± 37
Changed the Fur box in −10 region into a conserved Fur box from E. coli
250 ± 12
5466 ± 26
Modified the Fur box in −10 region into an enhanced Fur box 
196 ± 5
1885 ± 17
Integrated a designed Fur box in −35 region
200 ± 9
662 ± 19
Inserted the SD-fur-TT circuit in the vector
230 ± 9
3464 ± 24
Inserted the Pfur-SDfur-TT circuit in the vector
195 ± 7
2389 ± 26
Performance evaluation of pPfhuA1as an inducible tightly regulated expression system
In order to see the time response of ptPfhuA1 by which it is effectively turned on, GFP yields along with E. coli growth were detected. GFP expression in the ptPfhuA1-bearing culture was continuously increased until 9 h post induction (Figure 3B). During the first 1 h induction, the relative fluorescence value increased by 3.4 folds, which meant the expression circuit had been switched on. Afterwards, the GFP expression increased very slowly from 1 to 5 h, and underwent a jump from 5 to 9 h. At last, the GFP expression achieved its peak at 9 h post induction, corresponding to the beginning of stationary phase in E. coli growth. Thus, the whole expression circle for PfhuA1 lasted for 9 h. In other side, the GFP synthesis always lagged behind the E. coli growth. This is desirable in toxic gene expression to reduce the damage of protein products to cells.
Toxic protein expression under the control of modified PfhuA1
With the successful application of pYPfhuA1 in the regulation of E gene expression, its broad availability was the next concern. Then, other two toxic proteins were synthesized using pYPfhuA1 in E. coli host. The gef encoding Gef of E. coli belongs to the hok killer gene family in Gram-negative bacteria and its expression kills the cell from the inside by interfering with a vital function in the cell membrane . The mazF from a toxin-antitoxin module mazEF specifies a stable toxin that cleaves mRNA at a specific site(s) which is responsible for programmed cell death in E. coli. Here patPfhuA1 was used, which is a middle-low-copy number vector with ~30 copies per cell . Colony forming units (CFUs) of the transformants were detected in the presence of inducer (Figure 4B). Comparing with the negative control Top10, all the recombinant strains were grown similarly in the presence of iron (data are not shown), and no leaky expression phenotype was appeared. However, the growth of Top10/patPfhuA1gef and Top10/patPfhuA1mazF was repressed when cultured with 2,2’-dipyridyl as the result of gef and mazF expression. At 0 h, 2,2’-dipyridyl was added. During the first 2 h after induction, all three strains performed to adapt the iron-limiting condition and maintained in a similar density of 109 CFU/ml. However, after 2 h, their growth behaviors were obviously different. Top10 strain grew fast to 2.4×109 CFU/ml in exponential way until 9 h post induction. However, the strains Top10/patPfhuA1gef and Top10/patPfhuA1mazF grew exponentially until 7 h post induction and then kept at 6×108 CFU/ml and 1×109 CFU/ml, respectively. Therefore, the growth repression shown by recombinant strains revealed that toxic genes gef and mazF were functionally expressed after induction under the control of PfhuA1. Furthermore, pYPfhuA1 could be applied in different toxic gene expressions.
One of the most notable merits for this system is its tight regulation by iron-limitation signal. Escolar et al. once pointed out that the actual sequences recognized by Fur consisted of a minimal array of three conserved 6 bp (GATAAT) units which could be extended laterally by discrete additions of repeats of the same unit, thereby giving rise to new sites to which the repressor binds with a range of affinities . It was assumed that the affinity would vary depending on both the number of repeats present on each operator and the conservation of their sequences, which would allow an entire range and hierarchy of transcription responses depending on small changes in the iron status in the cell . Based on this hypothesis, we designed three types of Fur boxes differing in sequence conservation and 6-bp unit repeat number. It was found that the more repeats of conserved 6-bp units or the more conserved Fur box, the stronger affinity between Fur and the promoter, which was following by the tighter promoter. This verified the hypothesis mentioned above. In another side, this strategy also provides a new way for modifying other inducible promoters. However, a problem in this strategy is that the transcription efficiency decreases along with the tightness improvement. This is in common, because most of the regulatory systems have an intrinsic limitation in the range of induction, and attempts to mutate promoters to reduce basal expression usually result in concomitant reduction of induced levels. To overcome this drawback, Royo et al. used the nasF attenuator and the NasR-dependent anti-termination system from Klebsiella oxytoca to construct a novel expression circuit which could conditionally prevent undesired transcription from any transcriptional initiation signal, while keeping induced levels intact [10, 25]. This point of view gives us a new reference for further optimization of expression vector constructed here.
Although iron is one of the most abundant elements on Earth, in aerobic environments it is predominantly found as ferric (hydro)oxides that are relatively insoluble at neutral pH, and thus, ionic ferric (Fe3+) concentrations are exceedingly low . Therefore, as an iron limitation responsive vector, pYPfhuA1 has great potential in developing ‘suicidal’ containment system for environmentally relevant application. Furthermore, in our previous research, iron-limiting condition was verified as an in vivo stimulating signal, and promoters derived from iron uptake system could be induced in vivo. Therefore, this means PfhuA1 could be used in the construction of in vivo inducible antigen synthesis system, which could alleviate the toxicity or metabolic burden of the host strain and improve the immunogenicity for bacterial vector vaccine [27–32]. Moreover, in vivo regulated and delayed attenuation strategy has been developed for maintaining the invasive abilities of bacterial vector to the greatest extent, such that the recombinant vaccine has the ability as virulent wild-type strain to reach effector lymphoid tissues before display of attenuation to preclude onset of any disease symptoms . pYPfhuA1 could also be applied in construction of in vivo delayed attenuation vaccine. In all, pYPfhuA1 could be induced by iron chelators in medium, in vitro aerobic environments and in vivo animal host, and has extensive application in biotechnology area.
In all, as the expression vector particularly developed for toxic protein synthesis, pYPfhuA1 owns its notable advantages, which include easy control, tight regulation, relatively high efficiency and multi-environments response potential. Consequently, pYPfhuA1 has great potential in functional study of toxic genes, construction of biological containment system, or bacterial vector vaccine development.
Bacterial strains and growth conditions
Plasmids used in this study
pUC18 derivative with rrnBT1T2 terminator, Apr
pUT derivative, with a reporter gene gfp, Apr
pUC18 derivative with rrnBT1T2 terminator, pBBR1 ori., Apr
pUC18 derivative with rrnBT1T2 terminator, p15A ori., Apr
pUC18 derivative with rrnBT1T2 terminator, pAT153 ori., Apr
pUT derivative containing PfhuAgfp TT, Apr
pUT derivative containing PfhuA1gfp TT, Apr
pUT derivative containing PfhuA2gfp TT, Apr
pUT derivative containing PfhuA3gfp TT, Apr
pUT derivative containing PfhuAgfp TT and SD-fur TT from E. coli, Apr
pUT derivative containing PfhuAgfp TT and Pfurfur TT from E. coli, Apr
pUT derivative containing PfhuA1E TT, Apr
pUTb derivative containing PfhuA1E TT, Apr
pUTa derivative containing PfhuA1E TT, Apr
pUTat derivative containing PfhuA1E TT, Apr
pUTat derivative containing PfhuA1gef TT, Apr
pUTat derivative containing PfhuA1mazF TT, Apr
A reporter plasmid pUTtG constructed previously was used as the promoter screening vector . The primers from the former work were applied to amplify the candidate promoters from bacterium chromosomes or plasmids . The amplified promoter products, possessing the native start codon, Shine-Dalgarno sequence, and −35 and −10 promoter elements plus additional upstream bases, were inserted into pUTtG, and the resultant plasmids were transformed into E. coli Top10 named Top10/ptPXG (X mean different promoters) for iron-regulated promoter screening.
PCR primers used in this study
General DNA procedure and online analysis tools
General DNA operations were carried out following the standard protocols. Automated DNA sequencing and primer synthesis were completed by Life Technologies (Shanghai, China). PfhuA characteristic regions were predicted by online tools: BPROM-bacterial promoter prediction program from Softberry (http://www.linux1.softberry.com/berry.phtml?topic=bprom%26group=programs%26subgroup=gfindb), BDGP-Promoter (http://www.fruitfly.org/seq_tools/promoter.html) and SCOPE-Suite for Computational identification Of Promoter Elements .
GFP synthesis detection
Overnight cell cultures were inoculated (1:100, v/v) into fresh LB medium containing appropriate antibiotics and cultured in a shaker at 200 rpm and 37°C. At middle log phase typically with an optical density at OD600 = 0.8-1.0, 2, 2’-dipyridyl was added to induce the expression of GFP. After 20 hours or defined time points of iron-limiting induction, 1 ml cell culture sample was taken, centrifuged at 12,000 g for 3 min, washed and resuspended in PBS (pH 7.2) to the same OD600 value (OD600 = 1.0). For each sample, 100 μl of cell suspension was added into a 96-well flat-bottom polystyrene plate (Costar, USA) and measured with a fluorescence plate reader (TECAN, GENios Pro, Austria). Excitation wavelength was set at 485 nm and emission was detected at 535 nm.
Toxic gene expression detection
The E. coli strains Top10/ptPfhuA1E, Top10/pbPfhuA1E, Top10/paPfhuA1E, Top10/patPfhuA1gef and Top10/patPfhuA1mazF were overnight grown at 37°C in LB medium supplemented with 40 μM FeSO4 to ensure tight repression of toxic gene. To induce toxic gene expression, the culture was diluted to an OD600 of 0.1 and cultured in a shaker at 200 rpm and 37°C. At early log phase (OD600 = 0.3-0.4), 2,2’-dipyridyl was added into culture to create iron-limiting condition. Cell samples were taken at different time points after induction to measure both the OD600 and the colony formation units (CFUs) to determine the bacteria growth in iron-limiting medium.
This work was supported by the joint project, National Natural Science Foundation of China- Austrian Science Fund (30811130545).
- Kimelman A, Levy A, Sberro H, Kidron S, Leavitt A, Amitai G, Yoder-Himes DR, Wurtzel O, Zhu Y, Rubin EM: A vast collection of microbial genes that are toxic to bacteria. Genome Res. 2012, 22 (9): 1-8.
- Singh JS, Abhilash P, Singh H, Singh RP, Singh D: Genetically engineered bacteria: An emerging tool for environmental remediation and future research perspectives. Gene. 2011, 480 (1–2): 1-9.View Article
- Atlas RM: Molecular methods for environmental monitoring and containment of genetically engineered microorganisms. Biodegradation. 1992, 3 (2–3): 137-146.View Article
- Molin S: Environmental potential of suicide genes. Curr Opin Biotechnol. 1993, 4 (3): 299-305. 10.1016/0958-1669(93)90099-I.View Article
- Saida F, Uzan M, Odaert B, Bontems F: Expression of highly toxic genes in E. coli: special strategies and genetic tools. Curr Protein Pept Sci. 2006, 7 (1): 47-56. 10.2174/138920306775474095.View Article
- Jana S, Deb J: Strategies for efficient production of heterologous proteins in Escherichia coli. Appl Microbiol Biotechnol. 2005, 67 (3): 289-298. 10.1007/s00253-004-1814-0.View Article
- Choi YJ, Morel L, Le Francois T, Bourque D, Bourget L, Groleau D, Massie B, Miguez CB: Novel, versatile, and tightly regulated expression system for Escherichia coli strains. Appl Environ Microbiol. 2010, 76 (15): 5058-5066. 10.1128/AEM.00413-10.View Article
- Giacalone MJ, Gentile AM, Lovitt BT, Berkley NL, Gunderson CW, Surber MW: Toxic protein expression in Escherichia coli using a rhamnose-based tightly regulated and tunable promoter system. Biotechniques. 2006, 40 (3): 355-363. 10.2144/000112112.View Article
- Glascock C, J Weickert M: Using chromosomal lacIQ1 to control expression of genes on high-copy-number plasmids in Escherichia coli. Gene. 1998, 223 (1): 221-231. 10.1016/S0378-1119(98)00240-6.View Article
- Royo JL, Manyani H, Cebolla A, Santero E: A new generation of vectors with increased induction ratios by overimposing a second regulatory level by attenuation. Nucleic Acids Res. 2005, 33 (19): e169-e169. 10.1093/nar/gni168.View Article
- Anthony LC, Suzuki H, Filutowicz M: Tightly regulated vectors for the cloning and expression of toxic genes. J Microbiol Methods. 2004, 58 (2): 243-250. 10.1016/j.mimet.2004.04.003.View Article
- Voigt CA: Genetic parts to program bacteria. Curr Opin Biotechnol. 2006, 17 (5): 548-557. 10.1016/j.copbio.2006.09.001.View Article
- Allen CE, Schmitt MP: HtaA is an iron-regulated hemin binding protein involved in the utilization of heme iron in Corynebacterium diphtheriae. J Bacteriol. 2009, 191 (8): 2638-2648. 10.1128/JB.01784-08.View Article
- Hoffelner H, Rieder G, Haas R: Helicobacter pylori vaccine development: Optimisation of strategies and importance of challenging strain and animal model. Int J Med Microbiol. 2008, 298 (1–2): 151-159.View Article
- Lim JM, Hong M, Kim S, Oh D, Kang H, Kwon O: Iron chelator-inducible expression system for Escherichia coli. J Microbiol Biotechnol. 2008, 18 (8): 1357-1363.
- Baichoo N, Helmann JD: Recognition of DNA by Fur: a reinterpretation of the Fur box consensus sequence. J Bacteriol. 2002, 184 (21): 5826-5832. 10.1128/JB.184.21.5826-5832.2002.View Article
- Escolar L, Pérez-Martín J, de Lorenzo V: Binding of the Fur (ferric uptake regulator) repressor of Escherichia coli to arrays of the GATAAT sequence. J Mol Biol. 1998, 283 (3): 537-547. 10.1006/jmbi.1998.2119.View Article
- Guan LY, Mu W, Champeimont J, Wang QY, HZh W, Xiao JF, Lubitz W, Zhang YX, Liu Q: Iron-regulated lysis of recombinant Escherichia coli in host releases protective antigen and confers biological containment. Infect Immun. 2011, 79 (7): 2608-2618. 10.1128/IAI.01219-10.View Article
- Rogers MB, Sexton JA, DeCastro GJ, Calderwood SB: Identification of an operon required for ferrichrome iron utilization in Vibrio cholerae. J Bacteriol. 2000, 182 (8): 2350-2353. 10.1128/JB.182.8.2350-2353.2000.View Article
- Mey AR, Wyckoff EE, Kanukurthy V, Fisher CR, Payne SM: Iron and fur regulation in Vibrio cholerae and the role of fur in virulence. Infect Immun. 2005, 73 (12): 8167-8178. 10.1128/IAI.73.12.8167-8178.2005.View Article
- Xiao Y, Liu Q, Chen H, Zhang Y: A stable plasmid system for heterologous antigen expression in attenuated Vibrio anguillarum. Vaccine. 2011, 29 (40): 6986-6993. 10.1016/j.vaccine.2011.07.036.View Article
- Witte A, Wanner G, Sulzner M, Lubitz W: Dynamics of PhiX174 protein E-mediated lysis of Escherichia coli. Arch Microbiol. 1992, 157 (4): 381-388. 10.1007/BF00248685.View Article
- Gerdes K, Poulsen L, Thisted T, Nielsen A, Martinussen J, Andreasen P: The hok killer gene family in gram-negative bacteria. New Biol. 1990, 2 (11): 946-956.
- Amitai S, Yassin Y, Engelberg-Kulka H: MazF-mediated cell death in Escherichia coli: a point of no return. J Bacteriol. 2004, 186 (24): 8295-8300. 10.1128/JB.186.24.8295-8300.2004.View Article
- Medina C, Camacho EM, Flores A, Mesa-Pereira B, Santero E: Improved Expression Systems for Regulated Expression in Salmonella Infecting Eukaryotic Cells. PLoS One. 2011, 6 (8): e23055-10.1371/journal.pone.0023055.View Article
- Stumm W, Morgan JJ: Aquatic chemistry. 1996, New York, USA: John Wiley & Sons, 3
- Mu W, Guan L, Yan Y, Liu Q, Zhang Y: A novel in vivo inducible expression system in Edwardsiella tarda for potential application in bacterial polyvalence vaccine. Fish Shellfish Immunol. 2011, 31 (6): 1097-1105. 10.1016/j.fsi.2011.09.015.View Article
- Loessner H, Endmann A, Leschner S, Bauer H, Zelmer A, Westphal K, Weiss S, zur Lage S: Improving live attenuated bacterial carriers for vaccination and therapy. Int J Med Microbiol. 2008, 298 (1–2): 21-26.View Article
- Wang S, Li Y, Shi H, Sun W, Roland KL, Curtiss R: Comparison of a regulated delayed antigen synthesis system with in vivo-inducible promoters for antigen delivery by live attenuated Salmonella vaccines. Infect Immun. 2011, 79 (2): 937-949. 10.1128/IAI.00445-10.View Article
- John M, Crean TI, Calderwood SB, Ryan ET: In vitro and in vivo analyses of constitutive and in vivo-induced promoters in attenuated vaccine and vector strains of Vibrio cholerae. Infect Immun. 2000, 68 (3): 1171-1175. 10.1128/IAI.68.3.1171-1175.2000.View Article
- Morin CE, Kaper JB: Use of stabilized luciferase‐expressing plasmids to examine in vivo‐induced promoters in the Vibrio cholerae vaccine strain CVD 103‐HgR. FEMS Immunol Med Mic. 2009, 57 (1): 69-79. 10.1111/j.1574-695X.2009.00580.x.View Article
- Husseiny MI, Hensel M: Evaluation of an intracellular-activated promoter for the generation of live Salmonella recombinant vaccines. Vaccine. 2005, 23 (20): 2580-2590. 10.1016/j.vaccine.2004.11.035.View Article
- Curtiss R, Xin W, Li Y, Kong W, Wanda S, Gunn B, Wang S: New technologies in using recombinant attenuated Salmonella vaccine vectors. Crit Rev Immunol. 2010, 30 (3): 255-270. 10.1615/CritRevImmunol.v30.i3.30.View Article
- Miller JH: Experiments in molecular genetics. 1972, Cold Spring Harbor, New York, USA: Cold Spring Harbor Laboratory
- Carlson JM, Chakravarty A, DeZiel CE: SCOPE: a web server for practical de novo motif discovery. Nucleic Acids Res. 2007, 35 (suppl 2): W259-W264.View Article
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.