Early raise of BDNF in hippocampus suggests induction of posttranscriptional mechanisms by antidepressants
- Laura Musazzi†1,
- Annamaria Cattaneo†2,
- Daniela Tardito1,
- Alessandro Barbon3,
- Massimo Gennarelli2, 3,
- Sergio Barlati3,
- Giorgio Racagni1, 2 and
- Maurizio Popoli1Email author
© Musazzi et al; licensee BioMed Central Ltd. 2009
Received: 27 February 2009
Accepted: 13 May 2009
Published: 13 May 2009
The neurotrophin BDNF has been implicated in the regulation of neuroplasticity, gene expression, and synaptic function in the adult brain, as well as in the pathophysiology of neuropsychiatric disorders and the mechanism of action of antidepressants. Antidepressant treatments have been shown to increase the expression of BDNF mRNA, although the changes measured were found to be different depending on various factors. A few studies only have measured levels of BDNF protein after antidepressant treatments, and poor correlation was found between mRNA and protein changes. We studied the time course of expression of BDNF mRNA and protein during drug treatments, in order to elucidate the temporal profile of regulation of this effector and whether mRNA and protein levels correlate. Rat groups were treated for 1, 2 or 3 weeks with fluoxetine or reboxetine; in additional groups drug treatment was followed by a washout week (3+1). Total BDNF mRNA was measured by Real Time PCR, pro- and mature BDNF proteins were measured by Western blot.
We found that mature BDNF protein is induced more rapidly than mRNA, by both drugs in hippocampus (weeks 1–2) and by reboxetine in prefrontal/frontal cortex (week 1). The temporal profile of BDNF protein expression was largely inconsistent with that of mRNA, which followed the protein induction and reached a peak at week 3.
These results suggest that BDNF protein is rapidly elevated by antidepressant treatments by posttranscriptional mechanisms, and that induction of BDNF mRNA is a slower process.
Brain-Derived Neurotrophic Factor (BDNF) is an abundant neurotrophin regulating neuroplasticity, gene expression, synaptic function and cognition in the adult brain [1, 2], that has been implicated in the pathophysiology of various neuropsychiatric and neurodegenerative disorders as well as in the mechanism of action of antidepressant drugs [3–9]. It is generally assumed that antidepressant treatments increase the expression of BDNF, which, according to the neurotrophic hypothesis, represents an effector of changes in neuroplasticity and cellular resilience mediating the long-term therapeutic effect of antidepressants. In this framework BDNF has lately become, together with the activation of cAMP-response element binding protein (CREB), a sort of readout system in the study of antidepressant mechanisms [4, 5, 7, 8, 10–14].
However, several studies have shown that the changes in the expression of BDNF can be quite different depending on various factors, such as type of drug, dosage, route of administration, length of treatment (see ref.  for discussion). Moreover, only a few studies have measured levels of BDNF protein after antidepressant treatments [10, 15–19] and, when mRNA and protein levels have been measured at the same time, poor correlation was found between mRNA and protein changes [16, 17]. In addition, the rat BDNF gene has a complex structure with at least eight 5' exons that can be spliced to a single 3' exon containing the coding domain for the BDNF protein, generating 11 different transcripts according to the last nomenclature . This makes the effect of drug treatments on BDNF expression more complex to explain, because the changes in total BDNF can potentially be related to changes in different non-coding exons spliced to the coding exon (exon IX).
In this study, with regard to the action of antidepressants, we asked the following questions:
1. What is the time course of BDNF expression during antidepressant treatment? Is it consistent with the onset of therapeutic effect?
2. What is the expression level of BDNF after1 week-washout?
3. How does BDNF protein expression profile during treatment correlate with BDNF mRNA expression?
In order to assess the time-course of BDNF expression during long-term antidepressant treatments, we have treated rats with two different drugs endowed with complementary mechanisms: fluoxetine (FLX), a selective serotonin reuptake inhibitor (SSRI), and reboxetine (RBX), a selective norepinephrine reuptake inhibitor (NRI). The drug treatments were carried out for 1, 2 or 3 weeks and were followed by an additional washout week (3+1), that was added in order to study the fate of BDNF expression when antidepressant treatment is discontinued. In all these rat groups we assessed the expression of total BDNF at both mRNA and protein levels, measuring both the pro- and mature forms of BDNF [1, 21, 22]. We found that mature BDNF protein is induced in hippocampus more rapidly than mRNA, suggesting that antidepressants rapidly regulate BDNF at posttranscriptional level.
Distinct temporal profile of expression of total BDNF mRNA induced by fluoxetine and reboxetine
Expression profile of BDNF protein induced by fluoxetine and reboxetine consistent with early translation of mRNA
In P/FC (Fig. 2B–D), 2-Way ANOVA showed an effect of time (F4,53 = 7.69; p < 0.0001) and time/drug interaction (F4,53 = 4.07; p < 0.01) for pro-BDNF and an effect of time (F4,54 = 44.79; p < 0.0001) and drug (F1,54 = 12.53; p < 0.001) for mature BDNF. Indeed, both FLX and RBX up-regulated pro-BDNF, with the only exception of week 3 for FLX. Differently from the HPC, in P/FC the up-regulation of pro-BDNF by FLX sharply contrasted with the absence of changes in mature BDNF (Fig. 2B–D). As in HPC, RBX consistently up-regulated mature BDNF at all times.
Overall, the regulation of mRNA and protein was not consistent in HPC and P/FC for both drugs. No mature BDNF was induced by FLX in P/FC. The results showed that in HPC the level of mature BDNF protein reached a peak (weeks 1–2) before the actual peak of mRNA (week 3), suggesting early translation and processing of the protein during antidepressant treatment. The same was observed for RBX in P/FC.
The main findings of this work were that BDNF mRNA and protein levels during antidepressant treatments did not correlate, and that induction of mature protein preceded that of mRNA. The measurement of BDNF mRNA expression showed an overall similar temporal profile of induction with the two drugs. RBX acted somewhat faster, inducing BDNF mRNA already at week 1 in both areas, while FLX induced BDNF from week 2. The shape of this induction was also somewhat different; FLX progressively increased the expression up to week 3, while RBX induced BDNF in two distinct waves. After a precocious induction at week 1, the expression of BDNF was not further increased (P/FC) or was reduced (HPC) by RBX between week 1 and 2. A second, more pronounced, wave of BDNF induction by RBX was observed between week 2 and 3. Overall, RBX acted faster on BDNF transcription but FLX outsized RBX with regard to the maximum extent of mRNA induction. Interestingly, the washout week also differently affected the outcome of the two drug treatments in HPC, because it sharply reduced the effect of FLX in both areas and of RBX in P/FC, while in HPC of RBX-treated rats BDNF mRNA level was unchanged at week 3+1. In a companion study , we previously found that 1 washout week entirely abolished the activation of CREB induced by 3 weeks of RBX, but not of FLX, in line with the longer half-life of norFLX, the active FLX metabolite. This discrepancy between the effects of washout on activation of CREB vs. BDNF expression may underscore the complex regulation of the BDNF gene, that contains several different regulatory elements (other than CRE) at the level of promoters .
Surprisingly, the temporal profile of BDNF protein expression in HPC was largely inconsistent with that of mRNA expression (compare Fig. 1A with Fig. 2C). Although with FLX a significant increase of mRNA was not observed until week 2, in the same rats mature BDNF (14 kDa) was already significantly elevated at week 1 and reached maximal level at week 2, when mRNA level was still far below the maximal level attained at week 3. Instead, the level of pro-BDNF (32 kDa) in HPC was significantly elevated only at week 3. A possible explanation for this discrepancy is that, in the first several days of FLX treatment, BDNF mRNA already available is rapidly translated and the proBDNF produced is also rapidly processed to mature BDNF, suggesting that FLX has a more rapid effect on BDNF translation/processing than on the transcription of BDNF gene. Although more in line with the early induction of mRNA by RBX, the profile of induction of mature BDNF by this second drug was also inconsistent with the profile of mRNA induction, which reached maximal level only at week 3 (as for FLX). With RBX the level of mature BDNF was already maximal at week 1 and remained steadily elevated up to week 3, suggesting again a rapid effect of the drug on BDNF translation/processing, perhaps accelerated by the induction of new mRNA at week 1 by RBX.
In P/FC the effects of the two drugs on BDNF translation were different from each other. Maximal level of mature BDNF was induced by RBX already at week 1, as in the HPC, although the extent of increase was somewhat smaller than in HPC. Again, the temporal profile of the protein was inconsistent with that of mRNA, which reached maximal level only at week 3. By contrast, no mature BDNF was detected at any time point during FLX treatment, although moderate levels of proBDNF were attained, suggesting that in P/FC FLX was able to induce transcription of BDNF but little translation and no processing of the protein. Previous studies did not find induction of BDNF protein by FLX or other antidepressants in frontal cortex [15, 17]. This substantial lack of induction of BDNF protein by FLX in P/FC would suggest that downstream effects of this drug on synaptic plasticity are more powerful in HPC than in prefrontal/frontal areas of cerebral cortex.
The present results are in line with the several reports that previously showed up-regulation of BDNF mRNA induced by chronic antidepressant treatments, often increasing with the length of treatment [16, 25]. We also found here, as in two previous studies [16, 17], a lack of correlation between mRNA and protein, to the point that in HPC detectable levels of BDNF protein were raised faster than those of mRNA by the drug treatment. Overall, the differences in levels of BDNF protein we detected are somewhat bigger than in previous studies using the same drugs; this could be due to the fact that we measured here for the first time BDNF levels in a time-course treatment by Western analysis rather than by ELISA, a method that does not allow to distinguish between pro- and mature BDNF.
The present results would suggest that two antidepressants, FLX and RBX, have a particularly fast effect on posttranscriptional regulation of BDNF in HPC (RBX also in P/FC). It has already been suggested that BDNF synthesis could be regulated at posttranscriptional level , as found for other proteins . First, it is known that selected BDNF transcripts are located at dendrites, a mechanism that allows rapid regulation of translational control ; we recently found that in rats treated with FLX or RBX more BDNF mRNA is located at dendrites in CA3 hippocampal area (Tongiorgi E, Popoli M; unpublished). Second, it was recently suggested that microRNAs exert a robust effect on BDNF levels in mature human prefrontal cortex, participating in posttranscriptional fine tuning of BDNF expression . Third, it was recently shown that cleavage of pro- to mature BDNF is regulated by neuronal activity ; changing the processing of pro-BDNF by proteolytic enzymes could be an additional way whereby antidepressants rapidly increase the production of mature BDNF. This could explain why in HPC we found a faster rise of mature BDNF compared to pro-BDNF (compare Fig. 2A with Fig. 2C). Investigation of these mechanisms in subchronic and chronic treatments with antidepressants is warranted.
Finally, the finding that BDNF, a major readout effector, is rapidly elevated in HPC by antidepressant treatment may suggest that later downstream events are required to elicit systemic and behavioral effects of antidepressants. If BDNF production is faster than thus far envisaged, perhaps the correlates of the therapeutic effects should be searched in the actual cellular/molecular changes induced by the action of BDNF, such as the modifications in synaptic plasticity and the increment of neurogenesis. Further work is required to characterize these events.
The results of the present study show that BDNF protein is rapidly elevated by antidepressant treatments whereas the induction of BDNF mRNA is a slower process, suggesting that posttranscriptional mechanisms are involved in the mechanism of antidepressants.
Drug treatment and preparation of total extract
RNA isolation, cDNA synthesis and Real-time quantitative PCR
Total RNA was isolated with Trizol (Invitrogen, Carlsbad, CA). 2 μg of RNA were reverse transcribed using SuperScript™ II Reverse Transcriptase (Invitrogen) and 0.5 μg random hexamers (Promega, Madison, WI) in a 50 μl final volume. Real-time PCR was carried out using a LightCycler rapid thermal cycler System (Roche Diagnostics, Mannheim, Germany) and SYBR Green was used for the detection of double-strand DNA. The PCR reaction was set up into microcapillary tubes in a volume of 20 μl with 2 μl of cDNA and 2 μl of LightCycler DNA Master SYBR Green I (Roche Diagnostics). For BDNF quantitation, fwd AGCTGAGCGTGTGTGACAGT and rew ACCCATGGGATTACACTTGG primers were used, while for β-actin quantitation, used as internal control [32, 33], fwd GGGAAATCGTGCGTGACATT and rwd CGGATGTCAACGTCACACTT primers were used. The PCR program included an initial denaturation step during 30 s at 95°C, followed by 45 cycles of a 1 s melting step at 95°C, a 10 s annealing step at 58°C, and a 10 s elongation step at 72°C (all temperature transition rates were 20°C/s). At the end of each cycle, the fluorescence emitted by SYBR Green was measured. At the end of PCR reaction, samples were subjected to a temperature ramp (from 70°C to 95°C, 2°C/s) with continuous fluorescence monitoring for melting curve analysis. For each PCR product, a single narrow peak was obtained by melting curve analysis at the specific temperature. Each sample was assayed in duplicate and the analysis was performed with Light Cycler Relative Quantification Software. Samples containing no template were used as negative controls in each experiment.
Western blot Analysis
Western blot analysis was carried out as described previously , by incubating PVDF membranes, containing electrophoresed proteins from total homogenates, with polyclonal antibodies for BDNF 1:500 (Santa Cruz Biotechnologies, Santa Cruz, CA). A minimum of 2–3 replicates per rat were carried out in the different treatment groups (vehicle, FLX, RBX). Following incubation with peroxidase-coupled secondary antibodies, protein bands were detected by using ECL (GE Healthcare, Uppsala, Sweden). Standard curves were obtained by loading increasing amounts of samples on gels as described previously . All protein bands used were within the linear range of standard curves and normalized for β-actin level in the same membrane. Standardization and quantitation was as reported previously, except that Quantity One software (BioRad Laboratories, Hercules, CA) was used.
All data were analyzed by using two-way analysis of variance (2-Way ANOVA) for the variable treatment (FLX/RBX) and the variable time (Cnt, 1 week, 2 weeks, 3 weeks, 3+1 washout week). If the analysis of variance revealed significant group differences, post-hoc Bonferroni tests were carried out to elucidate the pattern of group differences. Significance was assumed at p < 0.05. Statistical analysis of the data was carried out using GraphPad Prism 4 (GraphPad Software, Inc., San Diego, CA).
This work was supported by Grants from the Ministry of University and Research, PRIN no. 2003053993 and 2005054953 (Italy).
- Lu B, Pang PT, Woo NH: The yin and yang of neurotrophin action. Nat Rev Neurosci. 2005, 6: 603-614. 10.1038/nrn1726.View ArticlePubMedGoogle Scholar
- Poo MM: Neurotrophins as synaptic modulators. Nat Rev Neurosci. 2001, 2: 24-32. 10.1038/35049004.View ArticlePubMedGoogle Scholar
- Castren E, Voikar V, Rantamaki T: Role of neurotrophic factors in depression. Curr Opin Pharmacol. 2007, 7: 18-21. 10.1016/j.coph.2006.08.009.View ArticlePubMedGoogle Scholar
- Duman RS, Monteggia LM: A neurotrophic model for stress-related mood disorders. Biol Psychiatry. 2006, 59: 1116-1127. 10.1016/j.biopsych.2006.02.013.View ArticlePubMedGoogle Scholar
- Martinowich K, Manji HK, Lu B: New insights into BDNF function in depression and anxiety. Nat Neurosci. 2007, 10: 1089-1093. 10.1038/nn1971.View ArticlePubMedGoogle Scholar
- Nestler EJ, Barrot M, DiLeone RJ, Eisch AJ, Gold SJ, Monteggia LM: Neurobiology of depression. Neuron. 2002, 34: 13-25. 10.1016/S0896-6273(02)00653-0.View ArticlePubMedGoogle Scholar
- Pittenger C, Duman RS: Stress, depression, and neuroplasticity: a convergence of mechanisms. Neuropsychopharmacol. 2008, 33: 88-109. 10.1038/sj.npp.1301574.View ArticleGoogle Scholar
- Russo-Neustadt AA, Chen MJ: Brain-derived neurotrophic factor and antidepressant activity. Curr Pharm Design. 2005, 11: 1495-1510. 10.2174/1381612053764788.View ArticleGoogle Scholar
- Tardito D, Perez J, Tiraboschi E, Musazzi L, Racagni G, Popoli M: Signaling pathways regulating gene expression, neuroplasticity and neurotrophic mechanisms in the action of antidepressants: a critical overview. Pharmacol Res. 2006, 58: 115-134. 10.1124/pr.58.1.7.Google Scholar
- Balu DT, Hoshaw BA, Malberg JE, Rosenzweig-Lipson S, Schechter LE, Lucki I: Differential regulation of central BDNF protein levels by antidepressant and non-antidepressant drug treatments. Brain Res. 2008, 1211: 37-43. 10.1016/j.brainres.2008.03.023.PubMed CentralView ArticlePubMedGoogle Scholar
- Donnici L, Tiraboschi E, Tardito D, Musazzi L, Racagni G, Popoli M: Time-dependent biphasic modulation of human BDNF by antidepressants in neuroblastoma cells. BMC Neurosci. 2008, 9: 61-10.1186/1471-2202-9-61.PubMed CentralView ArticlePubMedGoogle Scholar
- Khundakar AA, Zetterstrom TS: Biphasic change in BDNF gene expression following antidepressant drug treatment explained by differential transcript regulation. Brain Res. 2006, 1106: 12-20. 10.1016/j.brainres.2006.05.063.View ArticlePubMedGoogle Scholar
- Larsen MH, Hay-Schmidt A, Ronn LC, Mikkelsen JD: Temporal expression of brain-derived neurotrophic factor (BDNF) mRNA in the rat hippocampus after treatment with selective and mixed monoaminergic antidepressants. Euro J Pharmacol. 2008, 578: 114-122. 10.1016/j.ejphar.2007.08.050.View ArticleGoogle Scholar
- Malberg JE, Blendy JA: Antidepressant action: to the nucleus and beyond. Trends Pharmacol Sci. 2005, 26: 631-638. 10.1016/j.tips.2005.10.005.View ArticlePubMedGoogle Scholar
- Altar CA, Whitehead RE, Chen R, Wörtwein G, Madsen TM: Effects of electroconvulsive seizure and antidepressant drugs on brain-derived neurotrophic factor protein in rat brain. Biol Psychiatry. 2003, 54: 703-709. 10.1016/S0006-3223(03)00073-8.View ArticlePubMedGoogle Scholar
- De Foubert G, Carney SL, Robinson CS, Destexhe EJ, Tomlinson R, Hicks CA, Murray TK, Gaillard JP, Deville C, Xhenseval V, Thomas CE, O'Neill MJ, Zetterström TS: Fluoxetine-induced change in rat brain expression of brain-derived neurotrophic factor varies depending on length of treatment. Neurosci. 2004, 128: 597-604. 10.1016/j.neuroscience.2004.06.054.View ArticleGoogle Scholar
- Jacobsen JPR, Mork A: The effect of escitalopram, desipramine, electroconvulsive seizures and lithium on brain-derived neurotrophic factor mRNA and protein expression in the rat brain and the correlation to 5-HT and 5-HIAA levels. Brain Res. 2004, 1024: 183-192. 10.1016/j.brainres.2004.07.065.View ArticlePubMedGoogle Scholar
- Molteni R, Calabrese F, Cattaneo A, Mancini M, Gennarelli M, Racagni G, Riva MA: Acute Stress Responsiveness of the Neurotrophin BDNF in the Rat Hippocampus is Modulated by Chronic Treatment with the Antidepressant Duloxetine. Neuropsychopharmacol. 2009, 34: 1523-32. 10.1038/npp.2008.208.View ArticleGoogle Scholar
- Xu H, Richardson JS, Li XM: Dose-related effects of chronic antidepressants on neuroprotective proteins BDNF, bcl-2 and Cu/Zn-SOD in rat hippocampus. Neuropsychopharmacol. 2003, 28: 53-62. 10.1038/sj.npp.1300009.View ArticleGoogle Scholar
- Aid T, Kazantseva A, Piirsoo M, Palm K, Timmusk T: Mouse and rat BDNF gene structure and expression revisited. J Neurosci Res. 2007, 85: 525-535. 10.1002/jnr.21139.PubMed CentralView ArticlePubMedGoogle Scholar
- Mowla SJ, Farhadi HF, Pareek S, Atwal JK, Morris SJ, Seidah NG, Murphy RA: Biosynthesis and post-translational processing of the precursor to brain-derived neurotrophic factor. J Biol Chem. 2001, 276: 12660-12666. 10.1074/jbc.M008104200.View ArticlePubMedGoogle Scholar
- Nagappan G, Zaitsev E, Senatorov VV, Yang J, Hempstead BL, Lu B: Control of extracellular cleavage of ProBDNF by high frequency neuronal activity. Proc Natl Acad Sci USA. 2009, 106: 1267-1272. 10.1073/pnas.0807322106.PubMed CentralView ArticlePubMedGoogle Scholar
- Tardito D, Musazzi L, Tiraboschi E, Mallei A, Racagni G, Popoli M: Early induction of CREB activation and CREB-regulating signalling by antidepressants. Int J Neuropsychopharmacol. 2009, 1-15. 10.1017/S1461145709000376.Google Scholar
- West AE, Chen WG, Dalva MB, Dolmetsch RE, Kornhauser JM, Shaywitz AJ, Takasu MA, Tao X, Greenberg ME: Calcium regulation of neuronal gene expression. Proc Natl Acad Sci USA. 2001, 98: 11024-11031. 10.1073/pnas.191352298.PubMed CentralView ArticlePubMedGoogle Scholar
- Martínez-Turrillas R, Del Río J, Frechilla D: Sequential changes in BDNF mRNA expression and synaptic levels of AMPA receptor subunits in rat hippocampus after chronic antidepressant treatment. Neuropharmacol. 2005, 49: 1178-1188. 10.1016/j.neuropharm.2005.07.006.View ArticleGoogle Scholar
- Lipshitz HD, Smibert CA: Mechanisms of RNA localization and translational regulation. Curr Opin Genet Dev. 2000, 10: 476-488. 10.1016/S0959-437X(00)00116-7.View ArticlePubMedGoogle Scholar
- Tongiorgi E: Activity-dependent expression of brain-derived neurotrophic factor in dendrites: facts and open questions. Neurosci Res. 2008, 61: 335-346. 10.1016/j.neures.2008.04.013.View ArticlePubMedGoogle Scholar
- Mellios N, Huang HS, Grigorenko A, Rogaev E, Akbarian S: A set of differentially expressed miRNAs, including miR-30a-5p, act as post-transcriptional inhibitors of BDNF in prefrontal cortex. Hum Mol Genet. 2008, 17: 3030-3042. 10.1093/hmg/ddn201.PubMed CentralView ArticlePubMedGoogle Scholar
- Dow AL, Russell DS, Duman RS: Regulation of Activin mRNA and Smad2 phosphorylation by antidepressant treatment in the rat brain: effects in behavioral models. J Neurosci. 2005, 25: 4908-4916. 10.1523/JNEUROSCI.5155-04.2005.View ArticlePubMedGoogle Scholar
- Barbon A, Popoli M, La Via L, Moraschi S, Vallini I, Tardito D, Tiraboschi E, Musazzi L, Giambelli R, Gennarelli M, Racagni G, Barlati S: Regulation of editing and expression of glutamate alpha-amino-propionic-acid (AMPA)/kainate receptors by antidepressant drugs. Biol Psychiatry. 2006, 59: 713-720. 10.1016/j.biopsych.2005.10.018.View ArticlePubMedGoogle Scholar
- Glowinski J, Iversen LL: Regional studies of catecholamines in the rat brain. I. The disposition of [3H]norepinephrine, [3H]dopamine and [3H]dopa in various regions of the brain. J Neurochem. 1996, 13: 655-669. 10.1111/j.1471-4159.1966.tb09873.x.View ArticleGoogle Scholar
- Altieri M, Marini F, Arban R, Vitulli G, Jansson BO: Expression analysis of brain-derived neurotrophic factor (BDNF) mRNA isoforms after chronic and acute antidepressant treatment. Brain Res. 2004, 1000: 148-155. 10.1016/j.brainres.2003.12.028.View ArticlePubMedGoogle Scholar
- Alfonso J, Frick LR, Silberman DM, Palumbo ML, Genaro AM, Frasch AC: Regulation of hippocampal gene expression is conserved in two species subjected to different stressors and antidepressant treatments. Biol Psychiatry. 2006, 59: 244-251. 10.1016/j.biopsych.2005.06.036.View ArticlePubMedGoogle Scholar
- Tiraboschi E, Tardito D, Kasahara J, Moraschi S, Pruneri P, Gennarelli M, Racagni G, Popoli M: Selective phosphorylation of nuclear CREB by fluoxetine is linked to activation of CaM kinase IV and MAP kinase cascades. Neuropsychopharmacol. 2004, 29: 1831-1840. 10.1038/sj.npp.1300488.View 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.