# Statistical model of natural stimuli predicts edge-like pooling of spatial frequency channels in V2

- Aapo Hyvärinen
^{1}Email author, - Michael Gutmann
^{1}and - Patrik O Hoyer
^{1}

**6**:12

https://doi.org/10.1186/1471-2202-6-12

© Hyvärinen et al; licensee BioMed Central Ltd. 2005

**Received: **12 October 2004

**Accepted: **16 February 2005

**Published: **16 February 2005

## Abstract

### Background

It has been shown that the classical receptive fields of simple and complex cells in the primary visual cortex emerge from the statistical properties of natural images by forcing the cell responses to be maximally sparse or independent. We investigate how to learn features beyond the primary visual cortex from the statistical properties of modelled complex-cell outputs. In previous work, we showed that a new model, non-negative sparse coding, led to the emergence of features which code for contours of a given spatial frequency band.

### Results

We applied ordinary independent component analysis to modelled outputs of complex cells that span different frequency bands. The analysis led to the emergence of features which pool spatially coherent across-frequency activity in the modelled primary visual cortex. Thus, the statistically optimal way of processing complex-cell outputs abandons separate frequency channels, while preserving and even enhancing orientation tuning and spatial localization. As a technical aside, we found that the non-negativity constraint is not necessary: ordinary independent component analysis produces essentially the same results as our previous work.

### Conclusion

We propose that the pooling that emerges allows the features to code for realistic low-level image features related to step edges. Further, the results prove the viability of statistical modelling of natural images as a framework that produces quantitative predictions of visual processing.

## Background

A number of models approach the computational modelling of primary visual cortex by using two processing stages. First, there is a linear filtering with filters that are bandpass, oriented, and spatially localized. In some models, the outputs of the linear filters are half-wave rectified, but this difference is inessential because a rectification is done in the second stage anyway. The second stage then consists of pooling together rectified outputs of the first stage, so that cells that have the same orientation and frequency, as well as similar spatial locations, are pooled together. This pooling is then essentially a summation of rectified outputs of filters of different phases. These two processing steps are assumed to roughly correspond to simple and complex cells in V1, respectively. While there is controversy of the validity of such models, see e.g. [1–3], this is probably the simplest and most succesful approach.

Recent research has seen a number of models that attempt to explain these processing stages based on statistical modelling of natural images (ecologically valid input). First, application of independent component analysis (ICA) [4] or sparse coding [5] shows that the statistically optimal linear features of natural images are very similar to those computed in simple cells in V1 [6–12]. Second, application of a variant of ICA in which some pooling is done in a second stage leads to processing that is similar to what is done in complex cells [13]. Thus, models based on natural image statistics have been able to succesfully reproduce the above-mentioned two stages, and many well-known observations on V1.

It would be most useful if we could use this modelling endeavour in a *predictive* manner, so that we would be able to predict properties of cells in the visual cortex, in cases where the properties have not yet been demonstrated experimentally. This would give testable, quantitative hypotheses that might lead to great advances, especially in the research in extrastriate areas such as V2, whose function is not well understood at this point.

Here, we attempt to accomplish such predictive modelling in order to predict properties of a third processing step, following the two described above. Previously, we have applied a modification of the ICA / sparse coding model on the outputs of modelled complex cells whose input consisted of natural images [14]. The modification consisted of assuming that the coefficients in the generative decomposition, as well as the values of the higher-order features, were all non-negative.

We extend our previous results in two ways. The complex cells in our previous work were all constrained to have the same frequency, which was done in order to reduce the computational load. Here, we first report a technical advance: it is not necessary to make the assumptions of nonnegativity as in [14]. Thus, we are able to use the conventional, computationally optimized ICA algorithms, in particular the FastICA algorithm [15]. We are then easily able to incorporate complex cells of different frequencies in the input without exceeding available computational resources. This enables us to study whether some kind of interaction between different frequencies emerges in the statistically optimal higher-order representation.

## Results

### Experiment 1: Using ordinary ICA with no constraints

As described in Methods, we input a large number of natural image patches into model complex cells that computed the sum of squares of outputs of two simple cells, one odd-symmetric and the other even-symmetric. Then, we performed independent component analysis of the complex cell outputs using the FastICA algorithm.

*f*

_{2}in Figure 1. The purpose was to show that the results in [14] can be replicated using ordinary ICA methods.

**a**

_{ i }which show the contribution of the third-stage feature of index

*i*on the activities of complex cells. A collection of the obtained basis vectors is shown in Figure 2 for the nonlinearity

*g*

_{1}(see Table 1), visualized in the same way as in [14], see Methods. We can see the same kind of emergence of collinear features as in [14]. That is, the higher-order features code for the simultaneous activation of complex cells that together form something similar to a straight line segment.

Nonlinearities *g* used in FastICA. The nonlinearities probe the non-Gaussianity of the estimated components in different ways.

| Classic measure of sparseness |

| More robust variant of |

| Skewness (asymmetry) |

| Robust variant of |

Those coefficients that are clearly different from zero have almost always the same sign in a single basis vector. Defining the sign as explained in Methods, this means that the coefficients are essentially non-negative. We thus see that the constraint of non-negativity of the basis vectors imposed in [14] has little impact on the results: even without this constraint, the system learns basis vectors which are mainly non-negative.

*g*

_{4}(robust skewness) seems to lead to the largest number of long contours.

### Experiment 2: Emergence of pooling over frequencies

In the second experiment, the complex-cell set was expanded to include cells of three different preferred frequencies. In total, there were now 432 complex cells. We performed ICA on the complex-cell outputs when their input consisted of natural images. Thus, we obtained 432 higher-order basis vectors (features) **a**_{
i
}with corresponding activities *s*_{
i
}.

*P*

_{ i }in Equation (4) of Methods for the learned basis vectors. The distribution of this measure for natural image input and white Gaussian noise input is shown in Figure 6. The figure shows that frequency pooling according to this measure was essentially nonexistent for white Gaussian noise input, but relatively strong for many basis vectors when the input consisted of natural images. To express this more quantitatively, we computed the 99% quantile for the white Gaussian noise input. Then, 59% of the basis vectors for natural image input had a pooling index

*P*

_{ i }that was larger than this quantile. (For the 95% quantile the proportion was 63%.) Thus, we can say that more than half of the higher-order basis vectors, when learned from natural images, have a pooling over frequencies that is significantly above chance level.

*P*

_{ i }above its 99% quantile for white Gaussian noise). These are plotted in Figure 7. Visual inspection shows that in this subset, all basis vectors exhibit pooling over frequencies that respects the orientation tuning and collinearity properties.

We also examined quantitatively whether the higher-order features are tuned to orientation. We investigated which complex cell has the maximum weight in **a**_{
i
}for each *i* in each frequency band. When the input consisted of natural images, in 86% of the cells the maximally weighted complex cells were found to be located at the hot-spot (*x*_{
i
}, *y*_{
i
})* (i.e., point of maximum activity, see Methods for exact definition) and tuned to the preferred orientation of the higher-order feature for *every* frequency *f*. This shows how the higher-order features are largely selective to a single orientation. When the input consisted of Gaussian white noise, only 34% of the cells were found to be orientation-selective according to this criterion.

## Discussion

### Frequency channels and edges

What is the functional meaning of the pooling we have found? We propose that this spatially coherent pooling of multiple frequencies leads to representation of an edge that is more realistic than the band-pass edges given by typical Gabor filters [16]. Presumably, this is largely due to the fact that natural images contain many sharp, step-like edges that are not contained in a single frequency band. Thus, representation of such edges is difficult unless information from different frequency bands is combined.

In terms of frequency channels, the model predicts that frequency channels should be pooled together after complex cell processing. Models based on frequency channels and related concepts have been most prominent in image coding literature in recent years, both in biological and computer vision circles. The utility of frequency channels in the initial processing stages is widely acknowledged, and it is not put into question by our results – in fact, the statistical modelling framework does show that using band-pass simple and complex cells is statistically optimal [6, 13]. However, the question of when the frequency channels should be pooled or otherwise combined has received little attention [17, 18]. Our results point out that a statistically optimal way is to pool them together right after the complex cell "stage", and this pooling should be done among cells of a given orientation which form a local, collinear configuration.

### Related work

Several investigators have looked at the connection between natural image statistics, Gestalt grouping rules, and local interactions in the visual cortex [14, 19–21]. However, few has considered the statistical relations between features of different frequencies so far. It should be noted that some related work on interactions of different frequencies does exist in the models of contrast gain control [22].

Compared to our own previous work [14], the main difference seems to be in the frequency tuning of the model complex cells. In [14], the complex cells were all constrained to have the same spatial frequency tuning – just as in Experiment 1 of the present paper. Therefore, it was impossible to obtain results related to frequency pooling. It seems that any differences in the results are not due to differences in the statistical analysis of the complex-cell outputs or the natural image data set used, because in Experiment 1 of the present paper, we essentially replicated the results of [14]. The statistical model for analyzing the outputs of complex cells was somewhat different in our earlier work: the components *s*_{
i
}and the coefficients *a*_{
ki
}were constrained to be non-negative, following proposals by [23, 24]. However, this constraint seems to be immaterial, because even without imposing the constraint, the coefficients turned out to be essentially non-negative (after defining the global sign as described in Methods).

Recent measurements from cat area 18 (somewhat analogous to V2) emphasize responses to "second-order" or "non-Fourier" stimuli, typically sine-wave gratings whose amplitudes are modulated [17, 25]. These results and the proposed models are related to our results and predictions, yet fundamentally different. In the model in [25], a higher-order cell pools outputs of complex cells in the same frequency band to find contours that are defined by texture-like cues instead of luminance. The same cell also receives direct input from simple cells of a different frequency, which enables the cell to combine luminance and second-order cues. This is in stark contrast to higher-order cells in our model, which pool outputs of complex cells of different frequencies. They can hardly find contours defined by second-order cues; instead they seem to be good for coding broad-band contours. Furthermore, in [17, 25], any collinearity of pooling seems to be absent. This naturally leads to the question: Why are our predictions so different from these results from area 18? We suspect this is because it is customary to think of visual processing in terms of division into frequency channels – "second-order" stimuli are just an extension of this conceptualization. Therefore, not much attempt has been made to find cells that break the division into frequency channels according to our prediction. On the other hand, one can presume that the cells found in area 18 in [17, 25] are different from our predictions because they use a different coding strategy from the one used in our model, perhaps related to the temporal aspects of natural image sequences [26, 27].

Another closely related line of work is by Zetzsche and coworkers [28, 29] who emphasize the importance of decomposing the image information to local phase and amplitude information. The local amplitude is basically given by complex-cell outputs, whereas the physiological coding of the local phases is not known. An important question for future work is how to incorporate phase information in the higher-order units. Some models by Zetzsche et al actually predict some kind of pooling over frequencies, but rather directly after the simple cell stage (see Fig. 16 in [29]).

### Towards predictive modelling

The present results are an instance of predictive modelling, where we attempt to predict properties of cells and cell assemblies that have not yet been observed in experiments. To be precise, the prediction is that in V2 (or some related area) there should be cells whose optimal stimulus is a broad-band edge that has no sidelobes while being relatively sharp, i.e. the optimal stimulus is closer to a step-edge than the band-pass edges that tend to be optimal for V1 simple and complex cells. The optimal stimulus should also be more elongated [30, 31] than what is usually observed in V1, while being highly selective for orientation.

Statistical models of natural images offer a framework that lends itself to predictive modelling of the visual cortex. First, they offer a framework where we often see emergence of new kinds of feature detectors – sometimes very different from what was expected when the model was formulated. Second, the framework is highly constrained and data-driven. The rigorous theory of statistical estimation makes it rather difficult to insert the theorist's subjective expectations in the model, and therefore the results are strongly determined by the data. Third, the framework is very constructive. From just a couple of simple theoretical specifications, e.g. non-Gaussianity, natural images lead to the emergence of complex phenomena.

We hope that the present work as well as future results in the same direction will serve as a basis for a new kind of synergy between theoretical and experimental neuroscience.

## Conclusion

We have shown that pooling over complex cells of different frequency preferences emerges when we model the statistical properties of natural images. This is accomplished by applying ordinary ICA on a set of modelled complex cells with multiple frequencies, and inputting natural images to the complex cells. The resulting independent components, as represented by the corresponding basis vectors, code for simultaneous activation of complex cells that have similar orientations, form a collinear configuration, and span multiple frequencies. Thus, statistical modelling of natural stimuli leads to an interesting hypothesis on the existence of a new kind of cells in the visual cortex.

## Methods

### Data and statistical analysis

The natural images were 1008 gray-scale images of size 1024 × 1536 pixels from van Hateren's database, available at http://hlab.phys.rug.nl/imlib/index.html (category "deblurred") [8]. We manually chose natural images in the narrower sense, i.e. only wildlife scenes. From the source images, 50,000 image patches of size 24 × 24 pixels were randomly extracted. The mean grey value of each image patch was subtracted and the pixel values were rescaled to unit variance. The resulting image patch will be denoted by *I*(*x*, *y*).

The complex-cell model was similar to our previous work [14]. The filter bank consisted of a number of complex cells arranged on a 6 × 6 grid. Complex-cell responses *x*_{
k
}to natural images were modelled with a classical energy model:

*W*were normalized to unit norm for visualization purposes. The actual normalization we used in the experiments consisted of standardizing the variances of the complex cell outputs so that they were equal to unity for natural image input. The number of complex cells totalled

*K*= 36 × 4 × 3 = 432. Note, however, that in Experiment 1 we only used a single frequency band.

Independent component analysis (ICA) was performed on the vector **x** = (*x*_{1},...,*x*_{
K
}) using the FastICA algorithm [15]. The orthogonalization approach was symmetric. Different nonlinearities *g* were used, see Table 1. Thus we learned (estimated) a linear decomposition of the form

or in vector form

where the vector **a**_{
i
}= (*a*_{1i},...,*a*_{
ki
}) gives a higher-order basis vector. The *s*_{
i
}define the values of the higher-order features in the third cortical processing stage.

Note that the signs of the basis vectors are not defined by the ICA model [4], i.e. the model does not distinguish between **a**_{
i
}and -**a**_{
i
}because any change in sign of the basis vector can be cancelled by changing the sign of *s*_{
i
}accordingly. Here, we defined the sign for each vector **a**_{
i
}so that the sign of the element with the maximal absolute value was positive.

To obtain a baseline with which to compare our results, and to show which part of the results is due to the statistical properties of natural images instead of some intrinsic properties of our filterbank and analysis methods, we did exactly the same kind of analysis for 24 × 24 image patches that consisted of white Gaussian noise, i.e. the gray-scale value in each pixel was randomly and independently drawn from a Gaussian distribution of zero mean and unit variance. The white Gaussian noise input provides a "chance level" for any quantities computed from the ICA results.

### Analysis of the ICA results

We visualized the resulting higher-order basis vectors **a**_{
i
}following [14] by plotting an ellipse at each centrepoint of complex cells. The orientation of the ellipse is the orientation of the complex cell *k*, and the brightness of the ellipse is proportional to the *a*_{
ki
}coefficient of the basis vector **a**_{
i
}, using a gray-scale coding of coefficient values. In Experiment 1, i.e. the case with a single frequency band, we used this method directly to visualize each higher-order basis vector in a single display. In Experiment 2, i.e. the multifrequency case, we visualized each frequency band separately.

In Experiment 2, we are interested in the frequency pooling of complex cells in different higher-order features. We quantified the pooling over frequencies using a simple measure defined as follows. Let us denote by *a*_{
i
}(*x*, *y*, *θ*, *f*_{
n
}) the coefficient in the higher-order basis vector **a**_{
i
}that corresponds to the complex cell with spatial location (*x*, *y*), orientation *θ* and preferred frequency *f*_{
n
}. We computed a quantity which is similar to the sums of correlations of the coefficients over the three frequency bands, but normalized in a slightly different way. This measure *P*_{
i
}was defined as follows:

where the normalization constant *C*_{
m
}is defined as

and likewise for *C*_{
n
}.

For further analysis of the estimated basis vectors, we defined the preferred orientation of a higher-order feature. First, let us define for a higher-order feature of index *i* the hot-spot (*x*_{
i
}, *y*_{
i
})* as the centre location (*x*, *y*) of complex cells where the higher-order component *s*_{
i
}generates the maximum amount of activity. That is, we sum the elements of **a**_{
i
}that correspond to a single spatial location, and choose the largest sum. This allows us to define the tuning to a given orientation of a higher-order feature *i* by summing over the elements of **a**_{
i
}that correspond to the spatial hotspot and a given orientation; the preferred orientation is the orientation for which this sum is maximized. We also computed the length of a higher-order feature as described in [14].

It is also possible to perform an image synthesis from a higher-order basis vector. However, the mapping from image to complex-cell outputs is not one-to-one. This means that the generation of the image is not uniquely defined given the activities of higher-order features alone. A unique definition can be achieved by constraining the phases of the complex cells. We assume that only odd-symmetric Gabor filters are active. Furthermore, we make the simplifying assumptions that the receptive fields *W* in simple cells are equal to the corresponding basis vectors, and that all the elements in the higher-order basis vector are non-negative (or small enough to be ignored). Then, the synthesized image
for higher-order basis vector **a**_{
i
}is given by

where the square root cancels the squaring operation in the computation of complex-cell responses, and *H* denotes the set of indices that correspond to complex cells of the preferred orientation at the hotspot. Negative values of *a*_{
ki
}were set to zero in this synthesis formula.

## Declarations

### Acknowledgements

A. H. was funded by the Academy of Finland, Academy Research Fellow position and project #48593. P.O.H. was funded by the Academy of Finland, project #204826. M.G. would like to thank Rodney Douglas for supporting this collaborative effort between HIIT and the Institute of Neuroinformatics, Zurich.

## Authors’ Affiliations

## References

- Chance FS, Nelson SB, Abbott LF: Complex cells as cortically amplified simple cells. Nature Neuroscience. 1999, 2 (3): 277-282. 10.1038/6381.View ArticlePubMedGoogle Scholar
- Mechler F, Ringach DL: On the classification of simple and complex cells. Vision Research. 2002, 42 (8): 1017-33. 10.1016/S0042-6989(02)00025-1.View ArticlePubMedGoogle Scholar
- Kagan I, Gur M, Snodderly D: Spatial organization of receptive fields of V1 neurons of alert monkeys: comparison with responses to gratings. J Neurophysiol. 2002, 88: 2557-2574.View ArticlePubMedGoogle Scholar
- Hyvärinen A, Karhunen J, Oja E: Independent Component Analysis. 2001, Wiley Interscience, [http://www.cis.hut.fi/projects/ica/book/]View ArticleGoogle Scholar
- Field D: What is the goal of sensory coding?. Neural Computation. 1994, 6: 559-601.View ArticleGoogle Scholar
- Olshausen BA, Field DJ: Emergence of simple-cell receptive field properties by learning a sparse code for natural images. Nature. 1996, 381: 607-609. 10.1038/381607a0.View ArticlePubMedGoogle Scholar
- Bell A, Sejnowski T: The 'Independent Components' of Natural Scenes are Edge Filters. Vision Research. 1997, 37: 3327-3338. 10.1016/S0042-6989(97)00121-1.PubMed CentralView ArticlePubMedGoogle Scholar
- van Hateren JH, van der Schaaf A: Independent component filters of natural images compared with simple cells in primary visual cortex. Proc Royal Society, Ser B. 1998, 265: 359-366. 10.1098/rspb.1998.0303.View ArticleGoogle Scholar
- van Hateren JH, Ruderman DL: Independent component analysis of natural image sequences yields spatiotem poral filters similar to simple cells in primary visual cortex. Proc Royal Society, Ser B. 1998, 265: 2315-2320. 10.1098/rspb.1998.0577.View ArticleGoogle Scholar
- Hoyer PO, Hyvärinen A: Independent Component Analysis Applied to Feature Extraction from Colour and Stereo Images. Network: Computation in Neural Systems. 2000, 11 (3): 191-210. 10.1088/0954-898X/11/3/302.View ArticleGoogle Scholar
- Caywood M, Willmore B, Tolhurst D: Independent Components of Color Natural Scenes Resemble V1 Neurons in Their Spatial and Color Tuning. J of Neurophysiology. 2004, 91: 2859-2873. 10.1152/jn.00775.2003.View ArticlePubMedGoogle Scholar
- Wachtler T, Lee TW, Sejnowski TJ: Chromatic structure of natural scenes. J Opt Soc Am A. 2001, 18: 65-77.View ArticleGoogle Scholar
- Hyvärinen A, Hoyer PO: Emergence of phase and shift invariant features by decomposition of natural images into independent feature subspaces. Neural Computation. 2000, 12 (7): 1705-1720. 10.1162/089976600300015312.View ArticlePubMedGoogle Scholar
- Hoyer PO, Hyvärinen A: A multi-layer sparse coding network learns contour coding from natural images. Vision Research. 2002, 42 (12): 1593-1605. 10.1016/S0042-6989(02)00017-2.View ArticlePubMedGoogle Scholar
- Hyvärinen A: Fast and Robust Fixed-Point Algorithms for Independent Component Analysis. IEEE Transactions on Neural Networks. 1999, 10 (3): 626-634. 10.1109/72.761722.View ArticlePubMedGoogle Scholar
- Griffin L, Lillholm M, Nielsen M: Natural image profiles are most likely to be step edges. Vision Research. 2004, 44 (4): 407-421. 10.1016/j.visres.2003.09.025.View ArticlePubMedGoogle Scholar
- Mareschal I, Baker CL: A cortical locus for the processing of contrast-defined contours. Nature Neuroscience. 1998, 1 (2): 150-154. 10.1038/401.View ArticlePubMedGoogle Scholar
- Olzak L, Wickens T: Discrimination of complex patterns: orientation information is integrated across spatial scale; spatial-frequency and contrast information are not. Perception. 1997, 26: 1101-1120.View ArticlePubMedGoogle Scholar
- Geisler WS, Perry JS, Super BJ, Gallogly DP: Edge co-occurence in natural images predicts contour grouping performance. Vision Research. 2001, 41: 711-724. 10.1016/S0042-6989(00)00277-7.View ArticlePubMedGoogle Scholar
- Sigman M, Cecchi GA, Gilbert CD, Magnasco MO: On a common circle: Natural scenes and Gestalt rules. Proceedings of the National Academy of Science, USA. 2001, 98: 1935-1940. 10.1073/pnas.031571498.View ArticleGoogle Scholar
- Elder J, Goldberg R: Ecological statistics of Gestalt laws for the perceptual organization of contours. Journal of Vision. 2002, 2 (4): 324-353.View ArticlePubMedGoogle Scholar
- Schwartz O, Simoncelli E: Natural signal statistics and sensory gain control. Nature Neuroscience. 2001, 4 (8): 819-825. 10.1038/90526.View ArticlePubMedGoogle Scholar
- Paatero P, Tapper U: Positive Matrix Factorization: A Non-negative Factor Model with Optimal Utilization of Error Estimates of Data Values. Environmetrics. 1994, 5: 111-126.View ArticleGoogle Scholar
- Lee DD, Seung HS: Learning the parts of objects by non-negative matrix factorization. Nature. 1999, 401: 788-791. 10.1038/44565.View ArticlePubMedGoogle Scholar
- Mareschal I, Baker C: Temporal and spatial response to second-order stimuli in cat area 18. J Neurophysiol. 1998, 80: 2811-2823.PubMedGoogle Scholar
- Hurri J, Hyvärinen A: Simple-Cell-Like Receptive Fields Maximize Temporal Coherence in Natural Video. Neural Computation. 2003, 15 (3): 663-691. 10.1162/089976603321192121.View ArticlePubMedGoogle Scholar
- Hyvärinen A, Hurri J, Väyrynen J: Bubbles: A unifying framework for low-level statistical properties of natural image sequences. J Opt Soc Am A Opt Image Sci Vis. 2003, 20 (7): 1237-1252.View ArticlePubMedGoogle Scholar
- Zetzsche C, Krieger G: Nonlinear neurons and high-order statistics: New approaches to human vision and electronic image processing. Human Vision and Electronic Imaging IV (Proc. SPIE vol. 3644). Edited by: Rogowitz B, Pappas T. 1999, SPIE, 2-33.View ArticleGoogle Scholar
- Zetzsche C, Röhrbein F: Nonlinear and extra-classical receptive field properties and the statistics of natural scenes. Network: Computation in Neural Systems. 2001, 12: 331-350. 10.1088/0954-898X/12/3/306.View ArticleGoogle Scholar
- Polat U, Tyler C: What pattern the eye sees best. Vision Research. 1999, 39 (5): 887-895. 10.1016/S0042-6989(98)00245-4.View ArticlePubMedGoogle Scholar
- Gilbert CD, Wiesel TN: Intrinsic connectivity and receptive field properties in visual cortex. Vision Research. 1985, 25 (3): 365-374. 10.1016/0042-6989(85)90061-6.View ArticlePubMedGoogle 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.