- Methodology article
- Open Access

# Integrated application of uniform design and least-squares support vector machines to transfection optimization

- Jin-Shui Pan†
^{1}, - Mei-Zhu Hong†
^{2}, - Qi-Feng Zhou
^{3}, - Jia-Yan Cai
^{1}, - Hua-Zhen Wang
^{3}, - Lin-Kai Luo
^{3}, - De-Qiang Yang
^{3}, - Jing Dong
^{1}, - Hua-Xiu Shi
^{1}and - Jian-Lin Ren
^{1}Email author

**9**:52

https://doi.org/10.1186/1472-6750-9-52

© Pan et al; licensee BioMed Central Ltd. 2009

**Received:**28 January 2009**Accepted:**31 May 2009**Published:**31 May 2009

## Abstract

### Background

Transfection in mammalian cells based on liposome presents great challenge for biological professionals. To protect themselves from exogenous insults, mammalian cells tend to manifest poor transfection efficiency. In order to gain high efficiency, we have to optimize several conditions of transfection, such as amount of liposome, amount of plasmid, and cell density at transfection. However, this process may be time-consuming and energy-consuming. Fortunately, several mathematical methods, developed in the past decades, may facilitate the resolution of this issue. This study investigates the possibility of optimizing transfection efficiency by using a method referred to as least-squares support vector machine, which requires only a few experiments and maintains fairly high accuracy.

### Results

A protocol consists of 15 experiments was performed according to the principle of uniform design. In this protocol, amount of liposome, amount of plasmid, and the number of seeded cells 24 h before transfection were set as independent variables and transfection efficiency was set as dependent variable. A model was deduced from independent variables and their respective dependent variable. Another protocol made up by 10 experiments was performed to test the accuracy of the model. The model manifested a high accuracy. Compared to traditional method, the integrated application of uniform design and least-squares support vector machine greatly reduced the number of required experiments. What's more, higher transfection efficiency was achieved.

### Conclusion

The integrated application of uniform design and least-squares support vector machine is a simple technique for obtaining high transfection efficiency. Using this novel method, the number of required experiments would be greatly cut down while higher efficiency would be gained. Least-squares support vector machine may be applicable to many other problems that need to be optimized.

## Keywords

- Support Vector Machine
- Mean Square Error
- Transfection Efficiency
- Little Square Support Vector Machine
- Cationic Liposome

## Background

Central to life functions, protein expression in normal and diseased states is essential for quantifying altered patterns of gene expression. This is especially true in the era that the sequencing of human genome has been finished. To gain insights into protein expression, we have to transfect cells with kinds of expression vectors, based on plasmid, viral vector, or transposon, etc. Transfection may be one of the commonest but indispensable procedures for cellular biology. However, in the process of evolution, eukaryotic cells tend to have low transfection efficiency in order to protect their genomes from exogenous insults. Transfection difficulty manifests itself, especially in the cotransfection of mammalian cells. Theoretically, if the transfection efficiency of single kind of plasmid is E, which ranges from 0 to 1, the efficiency of double and triple cotransfection may decline to E^{2} and E^{3}, respectively. Therefore, it is of great importance to improve efficiency.

In order to enhance the transfection, several kinds of strategies are developed, which are categorized into two types: viral gene delivery carriers and non-viral gene delivery carriers. In non-viral gene delivery carriers, cationic liposomes has the widest application. Cationic liposomes are positively charged liposomes which interact with the negatively charged DNA molecules to form a stable complex. Cationic liposomes consist of a positively charged lipid and a co-lipid. A variety of positively charged lipid formulations are commercially available and many other are under development. Lipofection, one of the most frequently cited cationic lipids, was first reported by Felgner in 1987 to deliver genes to cells in culture [1]. Lipofection has been used to deliver linear DNA, plasmid DNA, and RNA to a variety of cells in culture. Liposomes offer several advantages in delivering genes to cells. (1) Liposomes have the ability to combine both with negatively and positively charged molecules. (2) Liposomes offer a degree of protection to the DNA from degradative processes. (3) Liposomes carry large pieces of DNA, potentially as large as a chromosome. (4) Liposomes can be targeted to specific cells or tissues. In addition, liposomes overcome problems inherent with viral vectors – specific concerns over immunogenicity and replication competent virus contamination. Liposomes resulted in a highly adaptable and flexible system capable of gene delivery both in vitro and in vivo. Current limitations regarding in vivo application of liposomes revolve around the low transfection efficiencies and transient gene expression. Also, liposomes display a small degree of cellular toxicity and appear to be inhibited by serum components. The ability to overcome these problems should greatly facilitate their application to a variety of gene delivery mechanisms.

Several factors have significant effects on the transfection efficiency of cationic liposomes, such as vigor of the host cells, the amount of plasmid, the amount of transfection agent, and the density of cells. However, it is hard to control vigor of host cells which has not a quantitative index. The other three factors are controllable in transfection, which can be adjusted according to the host cells and transfection agents. However, the adjustment of these three factors is a time and energy-consuming work. For most researchers, they may spend two to three months on optimizing transfection. Fortunately, several mathematical methods offer promising avenues to the resolution of this issue.

There are several ways to perform computer experiments, such as Latin Hypercube Sampling (LHS) and Uniform Design (UD). LHS was brought up by three scholars in the North American [2]. Uniform design, abbreviated as UD, was first developed by Fang *et al* in nineteen eighty [3]. UD seeks design points that are scattered uniformly on the domain. It has been popular since 1980. The main advantages of UD may be generalized as the following: first, it has the ability to greatly reduce the number of experiments while not to alter the representativeness; second, it generates a regression model based on the results and it's able to predict at what independent variables the dependent variable may gain the maximum.

As a relatively new algorithm used for classification and regression, support vector machine (SVM) was developed in the 1990s [4, 5]. It is a desired method for estimation based on finite-sample and therefore is able to solve a lot of practical problems in case of limited samples. Their practical successes may be attributed to solid theoretical foundations based on Vapnik Chervonenkis theory, and to the minimization of structural risk [6]. In order to implement the SVM into our transfection optimization, the least squares support vector machines (LS-SVM) was used, which has a growing popularity for regression problems [7]. It can be argued that LS-SVM would yield better generalization for regression problems on finite samples [8].

## Results and discussion

^{5}of cells, 0.66

*μ*g of plasmid and 1.32

*μ*g of LipofectAMINE were used. And this was accord with the observed data (Table 2 and Figure 2). More than that, there was a high degree of coincidence between calculated transfection efficiency and the deduced date from the model (Figure 3). Thus, by virtue of UD and LS-SVM, only 15 experiments, which can be performed in two 24-well plates, are needed to get the optimal transfection conditions whereas more than 15

^{3}experiments are needed to attain the expected purpose by using traditional method. What's more, if more accurate conditions were demanded, the number of experiments would greatly exceed 15

^{3}.

Protocol of experiments for generating the model

Sequence Number | Number of Cells (× 105)* | Amount of plasmid ( | Amount of LipofectAMINE ( | Expression ratio of GFP % | Fitted Values % |
---|---|---|---|---|---|

A | 11.0 | 100.50 | 80.88 | 64.74 | 64.50 |

B | 21.1 | 50.25 | 40.56 | 45.58 | 45.72 |

C | 31.2 | 140.70 | 111.12 | 68.37 | 68.06 |

D | 41.3 | 20.10 | 131.28 | 19.60 | 20.25 |

E | 51.4 | 120.60 | 10.32 | 30.61 | 31.04 |

F | 61.5 | 70.35 | 60.72 | 56.16 | 56.09 |

G | 71.6 | 80.40 | 151.44 | 75.24 | 74.80 |

H | 81.7 | 40.20 | 101.04 | 48.78 | 48.85 |

I | 91.8 | 150.75 | 50.64 | 59.19 | 59.06 |

J | 101.9 | 30.15 | 20.40 | 38.23 | 38.51 |

K | 112.0 | 110.55 | 141.36 | 82.53 | 81.94 |

L | 122.1 | 10.05 | 70.80 | 13.49 | 14.26 |

M | 132.2 | 130.65 | 90.96 | 71.67 | 71.29 |

N | 142.3 | 90.45 | 30.48 | 51.96 | 51.97 |

O | 152.4 | 60.30 | 121.20 | 61.85 | 61.67 |

Protocol of experiments for testing accuracy of the model

Sequence Number | Number of Cells (× 10 | Amount of plasmid ( | Amount of LipofectAMINE ( | Expression ratio of GFP (%) | Fitted Values (%) |
---|---|---|---|---|---|

A | 1.50 | 0.30 | 0.60 | 57.25 | 50.07 |

B | 1.60 | 0.36 | 0.72 | 63.56 | 58.98 |

C | 1.70 | 0.42 | 0.84 | 69.73 | 67.57 |

D | 1.80 | 0.48 | 0.96 | 77.64 | 75.13 |

E | 1.90 | 0.54 | 1.08 | 83.47 | 80.97 |

F | 2.00 | 0.60 | 1.20 | 86.42 | 84.47 |

G | 2.10 | 0.66 | 1.32 | 92.32 | 85.12 |

H | 2.20 | 0.72 | 1.44 | 85.19 | 82.59 |

I | 2.30 | 0.78 | 1.56 | 81.29 | 76.73 |

J | 2.40 | 0.84 | 1.68 | 73.61 | 67.62 |

*γ*= 42, C = 1500), and then the LS-SVM model obtains a peak estimated performance. After the optimal parameters for model construction are known, the according model (final model) is validated by predicting the validation data and comparing these predictions with the real observations.

MSE of LS-SVM model based on training parameters tuning

C | ... | 1200 | 1300 | 1400 | 1500 | 1600 | 1700 | ... |
---|---|---|---|---|---|---|---|---|

| ||||||||

... | ... | ... | ... | ... | ... | ... | ... | ... |

32 | 6.4129 | 6.4572 | 6.4789 | 6.511 | 6.5446 | 6.5762 | ... | |

36 | 6.2583 | 6.2794 | 6.2943 | 6.3152 | 6.323 | 6.3628 | ... | |

40 | 6.2155 | 6.207 | 6.2024 | 6.2068 | 6.2047 | 6.2135 | ... | |

42 | 6.2271 | 6.2013 | 6.1901 |
| 6.1815 | 6.1821 | ||

48 | 6.4039 | 6.3349 | 6.2821 | 6.2379 | 6.2074 | 6.1829 | ||

50 | 6.504 | 6.4205 | 6.3529 | 6.2984 | 6.254 | 6.2153 | ||

60 | 7.2761 | 7.1095 | 6.9685 | 6.848 | 6.7444 | 6.6548 | ||

... | ... | ... | ... | ... | ... | ... | ... | ... |

Coincidence between the observed data and the predicted values

real observation | predicted value | Net difference | Error ratio | |
---|---|---|---|---|

N1 | 57.25 | 51.91 | 5.34 | 0.093289 |

N2 | 63.56 | 59.86 | 3.70 | 0.058202 |

N3 | 69.73 | 67.43 | 2.30 | 0.033020 |

N4 | 77.64 | 74.10 | 3.54 | 0.045576 |

N5 | 83.47 | 79.38 | 4.09 | 0.048986 |

N6 | 86.42 | 82.81 | 3.61 | 0.041720 |

N7 | 92.32 | 84.04 | 8.28 | 0.089702 |

N8 | 85.19 | 82.81 | 2.38 | 0.027934 |

N9 | 81.29 | 79.03 | 2.26 | 0.027834 |

N10 | 73.61 | 72.74 | 0.87 | 0.011860 |

Contribution analysis of independent variables

variable | average predicated value | average MSE | |||||||||
---|---|---|---|---|---|---|---|---|---|---|---|

N1 | N2 | N3 | N4 | N5 | N6 | N7 | N8 | N9 | N10 | ||

x1 | 49.97 | 57.23 | 64.16 | 70.35 | 75.38 | 78.91 | 80.62 | 80.30 | 77.82 | 73.18 | 47.1928 |

x2 | 45.28 | 47.28 | 49.62 | 52.87 | 57.59 | 64.22 | 73.10 | 84.41 | 98.16 | 114.19 | 504.2573 |

x3 | 53.97 | 60.00 | 64.44 | 67.32 | 68.73 | 68.79 | 67.64 | 65.44 | 62.37 | 58.61 | 225.4837 |

Owing to UD, the amount of test points required can be enormously reduced, especially when the experimental region has many factors and multiple levels, while the results that reflect the major characteristics of the experimental system are ensured. As an efficient fractional factorial design, UD has been widely applied in manufacturing, system engineering, pharmaceutics, and natural sciences in the past decades [10–12]. The UD was used in this research to describe factors that significantly influence transfection efficiency to obtain a smaller, more manageable set. To perform a computer experiment, in order to have a wide coverage of the entire design region with a limited number of runs, UD is a good recommendation.

The SVM is a machine learning technique with a strong theoretical foundation that has been used to improve classification accuracy in biological applications [13–19]. The SVM is a maximum margin classifier that can solve non-linear classification problems by learning an optimal separating hyperplane in a higher-dimensional feature space. By use of non-linear kernel functions such as a Gaussian kernel, complex and non-linear decision functions can be learned by the SVM. LS-SVM is a reformulation to standard SVM. It is closely related to regularization networks and Gaussian processes but it additionally emphasizes and exploits primal-dual interpretations from optimization theory. In our experiment, LS-SVM mapped the original input space into a high dimension feature space by a Gaussian kernel and then learns a smoothest hyperplane to fit the training data. From statistical learning theory, it can be expected that this hyperplane would have excellent generalization ability and has minor local extreme value. Together, UD has the ability to greatly reduce the number of experiments while not to alter the representativeness and LS-SVM would yield better generalization for regression problems on finite samples. Thereupon, the integrated application of UD and LS-SVM would have high prediction accuracy and would contribute to transfection optimization.

## Conclusion

This paper investigates the integrated application of UD and LS-SVM to transfection, for obtaining precise information on the optimal conditions. Based on our experiments, UD and LS-SVM appear to have high efficiency and perform well even when undergone experiments are extremely scarce. With the established model, we are able to gain the optimal transfection conditions and the highest transfection efficiency that can be reached. Thus, the required time and experiments to improve transfection efficiency can be greatly reduced while the achieved efficiency may even be higher than traditional methods. It seems that LS-SVM has higher accuracy in the prediction of optimal transfection conditions than it does in the prediction of highest transfection efficiency; nevertheless, we usually have higher stringency of the information on optimal transfection conditions. It should be pointed out that the vigor of host cells and the purity of plasmid have crucial effect on the transfection efficiency too. However, these factors are uncontrollable in most settings. Further interpretation of the results obtained from other host cell lines is required. These issues are part of our ongoing research.

## Methods

### Cell Culture

The 293FT cell line was maintained in DMEM supplemented with 100 mL/L fetal calf serum, 2 mmol/mL L-glutamine, 100 *μ*g/mL penicillin and 100 units/mL streptomycin. The cells were incubated in a humidified incubator at 37°C containing 50 mL/L CO_{2}. Cell viability was estimated by the trypan blue dye exclusion method. The 293FT cells were seeded into 24-well plates 24 h prior to transfection. Three wells of cells were transfected for every experiment. The cells were transfected using LipofectAMINE 2000 (abbreviated as LipofectAMINE) cationic liposome (Invitrogen, Carlsbad, California, USA) and the cells were harvested 36 h after transfection. Transfection efficiency was evaluated by calculating the ratio of cells that express green fluorescent protein (GFP) by using flow cytometer (COULTER EPICS XL, Beckman, USA). The experiments were performed in duplicate.

### Uniform Design

On the basis of orthogonal design, UD as a new experimental design method was proposed by Fang *et al* in 1980s. The characteristics of UD are taking no account of regular comparability, completely ensuring the uniformity, and distributing the test points in the experimental scope adequately and uniformly. UD finds good representative points uniformly scattered over the sample space for a much more efficient parameter search. It is one kind of space filling designs that can be used for computer techniques. Suppose there are s samples of interest over a domain *C*^{
S
}. The goal here is to choose a set of m points *p*_{
m
}= {*θ*_{1}, ..., *θ*_{
m
}} ⊂ *C*^{
S
}such that these points are uniformly scattered on *C*^{
S
}.

### Experimentations

In a protocol consists of 15 experiments, amount of liposome, plasmid, and the number of seeded cells were set as independent variables while transfection efficiency was set as dependent variable. Each independent variable had 15 levels. The ranges of independent variables were set according to the instruction of manufacturer. The protocol was performed according to the principle of UD (Table 1). Each transfection efficiency (dependent variable) was calculated by flow cytometer. The expression of GFP in each experiment was also observed by fluorescence microscope (Eclipse 80I, Nikon, Tokyo, Japan). A model was constructed by using LS-SVM. The respective fitted value to each measured transfection efficiency was also deduced from the established model. Another protocol consisting of 10 experiments was designed centering on the predicted optimal conditions at which the dependent variable would reach the maximum (Table 2). And the observed GFP expression was shown in Figure 1 and Figure 2. All the observed data in Table 1 and Table 2 were the mean values of three independent experiments.

### Development of the LS-SVM based models for prediction of transfection efficiency

**x**

_{ i },

*y*

_{ i }), (

*i*= 1, ...,

*n*), where d-dimensional input is

**x**∈

*R*

^{ d }and the output is

*y*∈

*R*. In SVM regression formulations, the input X is first mapped into a m-dimensional feature space using some fixed (nonlinear) mapping, and then a linear model is constructed in this feature space [8]. Using mathematical notation, the linear model (in the feature space)

*f*(

**x**,

*ω*) is given by Equation (1), where

*g*

_{ j }(

**x**),

*j*= 1, ...,

*m*denotes a set of nonlinear transformations, and b is the "bias" term.

*L*(

*y*,

*f*(

**x**,

*ω*)). SVM regression uses a new type of loss function called

*ε*-insensitive loss function proposed by Vapnik [6]:

*ω*||

^{2}. In addition, it introduces (non-negative) slack variables

*ξ*

_{ i },

*i*= 1, ...

*n*to measure the deviation of training samples outside the

*ε*-insensitive zone. Thus, SVM regression is formulated as minimization of the following function:

*γ*should reflect the distribution/range of x-values of the training data.

The LS-SVM based model building processes were carried out by using the software package named LS-SVMlab1.5 available at http://www.esat.kuleuven.ac.be/sista/lssvmlab/. This toolbox provides in-depth functionality of SVM. Functions include tuning, optimizing, validating, and training SVMs. Significantly, it provides a good visual representation of the trained LS-SVMlab. The meta-parameters of the LS-SVM model with a Gaussian kernel function are *γ* (the width of the Gaussian kernels) and C (regularization factor). The model construction process consisted consecutively of: i) selection of the inputted variables, ii) selection of the training parameters (C and *γ*), iii) construction of the model, iv) performance evaluation by validation data.

### Contribution analysis of independent variables

The contribution of specific variable was evaluated by using a method based on LS-SVM. Random alteration of the value of a specific variable has a similar effect to take out of this variable. However, the feature space keeps unchanged when it undergoes this process. There are only three variables, so that we can analyze them one by one. The analysis process was described as follows: i) select a variable, such as the density of cells; ii) random exchange the 15 training samples value of cell density and keep other two variables (the amount of plasmid and the amount of LipofectAMINE) unchanged; iii) train another LS-SVM using the same parameters (*γ* = 42, C = 1500); iv) test the trained LS-SVM on 10 testing samples, get the predicated value and MSE; v) repeat step ii to step iv for 10 times, get the average predicated value and average MSE; vi) repeat step i to step v on other two variables. The variable with the biggest MSE has the most important effect on the depend variable.

## Notes

## Declarations

### Acknowledgements

We thank professor Tian-Hui Hu and associate professor Hong-Bin Dai for discussion and comments on the manuscript. The authors wish to thank anonymous referees for comments/suggestions that helped in improving the manuscript.

This research was supported by National Natural Science Foundation of China (No. 30750013) and Key Science Research Project Natural Science Foundation of Xiamen (No. WKZ0501).

## Authors’ Affiliations

## References

- Felgner PL, Gadek TR, Holm M, Roman R, Chan HW, Wenz M, Northrop JP, Ringold GM, Danielsen M: Lipofection: a highly efficient, lipid-mediated DNA-transfection procedure. Proc Natl Acad Sci USA. 1987, 84 (21): 7413-7417. 10.1073/pnas.84.21.7413.View ArticleGoogle Scholar
- McKay MD, Beckman RJ, Conover WJ: A Comparison of Three Methods for Selecting Values of Input Variables in the Analysis of Output From a Computer Code. Technometrics. 1979, 21: 239-245. 10.2307/1268522.Google Scholar
- Fang KT: The uniform design: application of number-theoretic methods in experimental design. Acta Mathematicae Applagatae Sinica. 1980, 3: 363-372.Google Scholar
- Vapnik VN: The Nature of Statistical Learning Theory [M]. 2000, New York: Springer-VerlagView ArticleGoogle Scholar
- Burges CJC: A Tutorial on Support Vector Machines for Pattern Recognition. Data Mining and Knowledge Discovery. 1998, 2 (2): 121-167. 10.1023/A:1009715923555.View ArticleGoogle Scholar
- Vapnik VN: An overview of statistical learning theory. IEEE TRANSACTIONS ON NEURAL NETWORKS. 1999, 10 (5): 988-999. 10.1109/72.788640.View ArticleGoogle Scholar
- Collobert R, Bengio S: SVMTorch: support vector machines for large-scale regression problems. J Mach Learn Res. 2001, 1: 143-160. 10.1162/15324430152733142.Google Scholar
- Suykens JAK, Vandewalle J: Least Squares Support Vector Machine Classifiers. Neural Processing Letters. 1999, 9 (3): 293-300. 10.1023/A:1018628609742.View ArticleGoogle Scholar
- Stone M: Cross-validatory choice and assessment of statistical predictions. Journal of the Royal Statistical Society. 1974, 36 (1): 111-147.Google Scholar
- Ji YB, Alaerts G, Xu CJ, Hu YZ, Heyden Vander Y: Sequential uniform designs for fingerprints development of Ginkgo biloba extracts by capillary electrophoresis. J Chromatogr A. 2006, 1128 (1–2): 273-281. 10.1016/j.chroma.2006.06.053.View ArticleGoogle Scholar
- Liu D, Wang P, Li F, Li J: Application of uniform design in L-isoleucine fermentation. Chinese journal of biotechnology. 1991, 7 (3): 207-212.Google Scholar
- Xu CP, Yun JW: Optimization of submerged-culture conditions for mycelial growth and exo-biopolymer production by Auricularia polytricha (wood ears fungus) using the methods of uniform design and regression analysis. Biotechnology and applied biochemistry. 2003, 38 (Pt 2): 193-199. 10.1042/BA20030020.View ArticleGoogle Scholar
- Baten AK, Halgamuge SK, Chang BC: Fast splice site detection using information content and feature reduction. BMC bioinformatics. 2008, 9 (Suppl 12): S8-10.1186/1471-2105-9-S12-S8.View ArticleGoogle Scholar
- Cheng CW, Su EC, Hwang JK, Sung TY, Hsu WL: Predicting RNA-binding sites of proteins using support vector machines and evolutionary information. BMC bioinformatics. 2008, 9 (Suppl 12): S6-10.1186/1471-2105-9-S12-S6.View ArticleGoogle Scholar
- Guo Y, Yu L, Wen Z, Li M: Using support vector machine combined with auto covariance to predict protein-protein interactions from protein sequences. Nucleic acids research. 2008, 36 (9): 3025-3030. 10.1093/nar/gkn159.View ArticleGoogle Scholar
- Han P, Zhang X, Norton RS, Feng ZP: Large-scale prediction of long disordered regions in proteins using random forests. BMC bioinformatics. 2009, 10 (1): 8-10.1186/1471-2105-10-8.View ArticleGoogle Scholar
- Kong L, Zhang Y, Ye ZQ, Liu XQ, Zhao SQ, Wei L, Gao G: CPC: assess the protein-coding potential of transcripts using sequence features and support vector machine. Nucleic acids research. 2007, W345-349. 10.1093/nar/gkm391. 35 Web ServerGoogle Scholar
- Krause L, McHardy AC, Nattkemper TW, Puhler A, Stoye J, Meyer F: GISMO – gene identification using a support vector machine for ORF classification. Nucleic acids research. 2007, 35 (2): 540-549. 10.1093/nar/gkl1083.View ArticleGoogle Scholar
- Verikas A, Gelzinis A, Bacauskiene M, Uloza V, Kaseta M: Using the patient's questionnaire data to screen laryngeal disorders. Computers in biology and medicine Comput Biol Med. 2009 Feb;39(2):148-55. 2009, 39 (2): 148-155.Google Scholar
- Cherkassky V, Ma Y: Practical selection of SVM parameters and noise estimation for SVM regression. Neural Networks. 2004, 17 (1): 113-126. 10.1016/S0893-6080(03)00169-2.View ArticleGoogle Scholar

## Copyright

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.