Skip to main content
  • Research article
  • Open access
  • Published:

MAPT rs242562 and GSK3B rs334558 are associated with Parkinson’s Disease in central China

Abstract

Background

Microtubule-associated protein tau (MAPT) is a neuronal protein involved in the pathogenesis of several neurodegenerative diseases including Parkinson’s Disease (PD). Glycogen synthase kinase 3 beta (GSK3B) catalyzes phosphorylation in multiple sites of tau protein. However, whether or not there is any association between the GSK3B gene alteration, MAPT haplotype and PD remains unexplored in Chinese population, especially in central Chinese population.

Results

Here, we aimed at studying the effect of MAPT rs242562 and GSK3B rs334558 on the risk of PD by performing a case-control association study in central China. Our data showed that all PD patients and controls were of MAPT H1/H1 diplotype in our study, thus confirming that the distribution of the MAPT H1 haplotype is common in China. GG genotype of MAPT rs242562 serves protection effect on PD risk in central Chinese population, while genotype of GSK3B rs334558 showed no difference between PD patients and controls.

Conclusions

We conclude that the MAPT rs242562 is associated with PD in central China in the background of MAPT H1/H1 diplotype. The GG genotype of rs242562 displays protection against PD in subgroup with GSK3B rs334558 T carrier.

Background

Parkinson’s disease (PD) is a common neurodegenerative disease in elders, characterized by the loss of dopaminergic neurons in substantial nigra and the formation of Lewy bodies [1]. PD affects approximately 1% of population aging over 60 years. Disease onset before 50 years of age defines as early-onset Parkinson's disease (EOPD) while that over 50 years of age, as late onset Parkinson’s disease (LOPD). The etiology and pathogenesis of PD, mainly considered to be related with aging, genetic and environmental factors, is incompletely understood. Recently, linkage analysis and genomewide association studies have shown that PARK1-PARK18, POLG, Glycogen Synthase kinase-3B (GSK3B), Microtubule-associated protein tau (MAPT) and other genes are associated with PD risk.

The MAPT gene is localized in a region of extended linkage disequilibrium (LD) on chromosome 17q21, containing 16 exons and encoding microtubule-associated proteins tau. Tau proteins, constituting a family of six isoforms ranging between 352-441 amino acids, are widely expressed in the nervous system. Excessive phosphorylation of tau proteins leads to abnormal intracellular aggregation, formatting of the double helix fiber filament, producing the nerve fiber tangles, and finally neurodegeneration [2]. MAPT gene, defines two extended haplotype, H1 and H2. The determination of MAPT haplotype H1/H2 is attributed to the presence of a 238 bp deletion between exons 9 and 10 of MAPT H2 haplotype [3]. H1 haplotype is reportedly related to PD susceptibility [4]. Concerning sub-haplotypes, there is a positive correlation between MAPT H1 haplotype and PD, including SNP rs242562 (A/G) and rs2435207 (G/A). The A-A sub-haplotype is associated with PD patients in a Norwegian population, while the G-A sub-haplotype for these two SNPs is associated with PD subjects in Greece [57]. However, at present, little is known about the correlation between MAPT gene polymorphism and PD patients in China, especially in central China area. Central China includes six widely agricultural provinces namely Hunan, Hubei, Jiangxi, Henan, Shanxi and Anhui.

GSK3B protein is widely expressed in all tissues, with particularly abundant levels in the brain [8]. The fundamental role of GSK3B in intracellular neuronal signaling systems is underpinned by its ability to phosphorylate several proteins that contribute to the structural characteristics and dynamics of neuronal cells. In addition, it has been proven that GSK3B catalyzes phosphorylation in multiple sites of tau protein [9], and it is an important pathogenic protein kinase for PD. Previous data indicate that the genetic alteration of GSK3B and its interaction with MAPT haplotypes are collectively related to PD morbidity rate in a Greek cohort [3, 10]. However, the association between the GSK3B gene alteration, MAPT haplotype and PD has not been previously explored in Chinese population.

In this study, our database contains 211 PD patients and 279 controls from central China. Our study examined the possible association between GSK3B promoter single nucleotide polymorphism (SNP) rs334558, MAPT haplotype, MAPT intron SNP rs242562 and PD in central China.

Methods

Subjects

We recruited 211 unrelated sporadic PD patients (average age at diagnosis: 57.23 ± 12.9 years; 137 males, 74 females) diagnosed according to the criteria of the UK Parkinson's Disease Society Brain Bank in Union hospital, Tongji Medical College, Huazhong University of Science and Technology (TJMC & HUST). The control group consisted of 279 individuals (average age: 53.72 ± 13.32 years; 178 males, 101 females). This study was approved by the Ethical Committee of TJMC & HUST. The ethics and written informed consent was obtained from all subjects.

Genotyping

Genomic DNA was extracted from venous blood using standard methods. Our primary objective was to select MAPT polymorphisms previously suggested to contribute to the risk of developing idiopathic PD in Chinese subjects. Therefore, we focused on the H1/H2 insertion/deletion polymorphism and rs242562 [11]. PCR products were generated with 50 ng DNA template in 12.5 μl 1 × taq PCR Master MIX (Bioteke Corporation, Beijing, China), 1 μl of 10 μmol/l each primer (Invitrogen, Carlsbad, CA, USA) and 12 μl Ultra-distilled water in a total volume of 25 μl. The PCR conditions used for the H1/H2 insertion/deletion polymorphism were: the H1/H2 insertion/deletion polymorphism forward primer GGAAGACGTTCTCACTGATCTG, reverse primer AGGAGTCTGGCTTCAGTCTCTC, initial denaturation at 94°C for 5 minutes followed by 33 cycles of 94°C for 30 seconds, 59°C for 30 seconds, 72°C for 30 seconds, and a final extension at 72°C for 10 minutes. The PCR conditions used for the rs242562 polymorphism were: rs242562 forward primer: CAGCCTTCCCTGTCCTTGATTC, rs242562 reverse primer; GCCTTCCCAACAGAGCAACC, initial denaturation at 94°C for 5 minutes, followed by 35 cycles of 94°C for 30 seconds, 58°C for 30 seconds and 72°C for 30 seconds, and a final extension at 72°C for 10 minutes. Digestion with XhoI restriction enzyme (New England Biolabs, Ipswich, MA, USA) at 37°C yielded a 385 bp band for the G-allele and a 287 bp and a 98 bp for the A-allele. The PCR conditions used for the rs334558 polymorphism were: rs334558 forward primer: GACGTCCGTGATTGGCTC, reverse primer: AGCCCAGAGCCCTGTCAG, initial denaturation at 94°C for 5 minutes, followed by 33 cycles of 94°C for 30 seconds, 62°C for 30s and 72°C for 15 seconds, with a final extension at 72°C for 10 minutes [12]. Digestion with AluI restriction enzyme (New England Biolabs) at 37°C yields a 344 bp band for the C allele and 220 and 124 bp for the T allele. Digestion products were resolved on a 2% agarose gel, stained in ethidium bromide solution and visualized with an ultraviolet light.

Statistical analysis

All statistical analysis was performed by SAS 9.2 (SAS Institute, Cary, NC). A log-odds ratio (OR) and its 95% confidence interval (CI) were estimated with logistic regression for allele wise, genotype-wise recessive/dominant genetic model of each polymorphism. The gene-gene interaction test was conducted by running two models: one with and the other without the interaction term (A-carrier * T-carrier). In addition, a trend test in logistic regression for allelic model was further employed for the effect of MAPT rs242562 and GSK3B rs334558 on risk for PD [13]. A probability (P) of less than 0.05 is considered statistically significant and P-values are two-tailed. Adjustment of P values for multiple testing uses Bonferroni correction. The power was calculated by Power V3.0 software (http://dceg.cancer.gov/tools/design/POWER) with the observed odds rates and minor allele frequencies for both the 2-df overall association test and the trend test [14].

Results

The allele and genotype distribution of MAPT rs242562 and GSK3B rs334558 in PD patients and controls were summarized in Table 1. The schematic of the MAPT region and associated polymorphisms including rs242562 and rs2435207 were shown in Figure 1. Genotype distributions of rs242562 and rs334558 followed Hardy-Weinberg equilibrium for PD patients and controls, no deviations from equilibrium were observed. All the PD patients and controls were of H1/H1 diplotype in our study.

Figure 1
figure 1

Schematic of the MAPT region and associated polymorphisms. Previous studies have identified a ~970 kb inversion polymorphism at chromosome 17q21.31, a region that contains MAPT and several other genes. Chromosomes with the inverted segment in different orientations represent two highly divergent MAPT haplotypes, H1 (direct orientation) and H2 (inverted orientation). No recombination has been identified between these two haplotypes over a region of ~1.5 Mb, but they have accumulated sequence variation independently. Polymorphisms including rs242562 and rs2435207 spanning MAPT exons 1 to 4 are specific to the H1 haplotype. Besides, there is a characteristic 238-bp deletion in the MAPT intron 9 of the H2 haplotype.

Table 1 Allele and Genotype frequencies of MAPT rs242562 and GSK3B rs334558 in PD cases and controls

For MAPT rs242562, we found positive results for genotype-wise dominant genetic model in involved subjects (Table 2, shown in bold) with a power of 74.5%. Genotype AG + AA is associated with the higher PD risk in overall involved subjects (AG + AA vs. GG, OR = 1.785, P = 0.016), compared to genotype GG, with A as the “putative” risk allele (Table 2). These data suggest that AG + AA genotype of rs242562 may confer risk effect on PD risk in central Chinese population. The ORs and 95% CIs for the case-control studies were not significant for allele-wise model and the recessive model in overall involved subjects with a lower power of 0.048-0.266 (Table 2). Additionally, the test for trend in logistic regression for allelic model also found no significant association with a power of 0.143 (Table 3), which is consistent with the results for allelic model in Table 2.

Table 2 Effect of MAPT rs242562 and GSK3B rs334558 on risk for PD
Table 3 The allelic-wise analysis with a test for trend for the effect of MAPT rs242562 and GSK3B rs334558 on risk for PD

For GSK3B rs334558, the allele-wise and genotype-wise ORs lack any statistical significance in overall involved subjects with an insufficient power of 0.029-0.064 (Table 2). The test for trend in logistic regression for allelic model with a power of 0.036 (Table 3) also found null association consistent with the allele-wise model in Table 2.

The gene-gene interaction test was conducted by running two models: with and without the interaction term (A-carrier * T-carrier), and likelihood ratio test indicated that interaction term was not statistically significant (P = 0.2140). The effect of GG and age on PD risk remained the same between the two models. We then studied the association between rs242562 and PD at different levels of the rs334558 genotype (Additional file 1: Table S1), and the association between rs334558 and PD at different levels of the rs242562 genotype (Additional file 2: Table S2). The results suggest that the GG genotype of rs242562 displays protection against PD in subgroup with T carrier, although the stratified analysis has a smaller sample size and a lower statistic power.

Discussion

In this study, we investigated the association of MAPT rs242562 and GSK3B rs334558 with PD in central China. The main findings of our study in central Chinese population include 1) All the PD patients and controls were in H1/H1 diplotype in our study, 2) GG genotype of MAPT rs242562 serves protection effect on PD risk in central Chinese population, 3) Genotype of GSK3B rs334558 showed no difference between PD patients and controls, and 4) GG genotype of rs242562 may display strong protective effect against PD risk in subgroup with GSK3B rs334558 T carrier.

Our data confirmed that the distribution of the MAPT H1 haplotype is common in China, which is in agreement with previous findings [15, 16]. Accumulating evidence showed that MAPT H1 haplotype-carrier group is a susceptible group of PD [16, 17]. This case-control study served as an association study based on population with susceptible haplotype H1 in central China. MAPT rs242562 regulates protein encoding despite it is in the intron of MAPT gene. The frequency of rs242562 GG genotype is higher in controls compared to that in cases, indicating it may play a protective role. These data are different from a previous study which showed that MAPT polymorphism rs242562 revealed no significant difference between PD patients and controls in Germany, Serbia and Greece [6, 18, 19]. The sample size, ethnic groups and environment may contribute to the different results.

GSK3B rs334558 is reportedly associated with PD, but the conclusions were inconsistent. A Greek study was the first to show that GSK3B rs334558 was related to PD, while CC served a protective effect and TT was overexpressed in PD [18]. Moreover, a study from Australia showed that the frequency of TT and H1/H1 diplotype in PD patients was significantly higher compared to control subjects [10]. However, the study about LOPD patients in India concluded that CC was a risk factor [20]. In the Spanish population, TT genotype and MAPT H1/H1 diplotype were associated with a decreased risk for PD [21]. In the present study, we did not find any significant difference in allele-wise and genotype-wise analysis for rs334558 between PD patients and controls. However, GG genotype of rs242562 may serve a strong protective effect on PD risk in subgroup with T carrier, although the stratified analysis has a small sample size and a low statistic power.

Sex may be involved in the association results of GSK3B and neurodegenerative disease [22]. In an Alzheimer's disease research, it was revealed that estrogen receptor and GSK3B could change the tau protein phosphorylation, indicating sex involved in pathogenesis of this neurodegenerative disease [22]. However, following sex stratification the sample size was smaller in our study, and whether hormonal or hormonal receptors change the association of GSK3B gene with PD needs to be further investigated.

Although MAPT rs242562 and GSK3B rs334558 are located in two different chromosomes, the encoded proteins interact with each other. Given the findings that H1 haplotype may act synergistically with variants in the MAPT and GSK3B genes in conferring risk for PD, gene-gene interactions will also be important to consider as they may provide critical insights into mechanisms of disease susceptibility. In this study, lower power was observed for rs334558 analysis. Thus, the null association should attribute to the insufficient power. Since the sample size of the genotype model including indicator variables of the two polymorphisms is considerably smaller, a larger sample size should be needed in order to further assess this association.

Conclusions

In summary, the MAPT rs242562 is associated with PD in central China, and GG genotype of rs242562 may provide protective effect against PD risk in subgroup with GSK3B rs334558 T carrier, although our results are limited by the sample size. Since MAPT gene and GSK3B gene confer genetic risk for PD, exploring how the gene-gene or gene-environment interactions contribute to dopaminergic neurodegeneration should be further considered. Such knowledge about the mechanisms could open up new windows for early diagnostic and therapeutic interventions in this important neurodegenerative disease.

Consent

Written informed consent was obtained from all patient and control for the association study. A copy of the written consent is available for review by the Series Editor of this journal.

Abbreviations

MAPT:

Microtubule-associated protein tau

PD:

Parkinson's disease

LOPD:

Late onset of Parkinson's disease

EOPD:

Early onset of Parkinson's disease

GSK3B:

Glycogen synthase kinase 3β

TJMC & HUST:

Tongji medical college, Huazhong University of Science

OR:

Odds ratios.

References

  1. Forno LS: Neuropathology of Parkinson's disease. J Neuropathol Exp Neurol. 1996, 55 (3): 259-272. 10.1097/00005072-199603000-00001.

    Article  CAS  PubMed  Google Scholar 

  2. Ballatore C, Lee VM, Trojanowski JQ: Tau-mediated neurodegeneration in Alzheimer's disease and related disorders. Nat Rev Neurosci. 2007, 8 (9): 663-672.

    Article  CAS  PubMed  Google Scholar 

  3. Kalinderi K, Fidani L, Katsarou Z, Clarimon J, Bostantjopoulou S, Kotsis A: GSK3beta polymorphisms, MAPT H1 haplotype and Parkinson's disease in a Greek cohort. Neurobiol Aging. 2011, 32 (3): 546-e541-545

    Article  PubMed  Google Scholar 

  4. Baker M, Litvan I, Houlden H, Adamson J, Dickson D, Perez-Tur J, Hardy J, Lynch T, Bigio E, Hutton M: Association of an extended haplotype in the tau gene with progressive supranuclear palsy. Hum Mol Genet. 1999, 8 (4): 711-715. 10.1093/hmg/8.4.711.

    Article  CAS  PubMed  Google Scholar 

  5. Fidani L, Kalinderi K, Bostantjopoulou S, Clarimon J, Goulas A, Katsarou Z, Hardy J, Kotsis A: Association of the Tau haplotype with Parkinson's disease in the Greek population. Mov Disord. 2006, 21 (7): 1036-1039. 10.1002/mds.20864.

    Article  PubMed  Google Scholar 

  6. Refenes N, Bolbrinker J, Tagaris G, Orlacchio A, Drakoulis N, Kreutz R: Role of the H1 haplotype of microtubule-associated protein tau (MAPT) gene in Greek patients with Parkinson's disease. BMC Neurol. 2009, 9: 26-10.1186/1471-2377-9-26.

    Article  PubMed Central  PubMed  Google Scholar 

  7. Skipper L, Wilkes K, Toft M, Baker M, Lincoln S, Hulihan M, Ross OA, Hutton M, Aasly J, Farrer M: Linkage disequilibrium and association of MAPT H1 in Parkinson disease. Am J Hum Genet. 2004, 75 (4): 669-677. 10.1086/424492.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  8. Embi N, Rylatt DB, Cohen P: Glycogen synthase kinase-3 from rabbit skeletal muscle, Separation from cyclic-AMP-dependent protein kinase and phosphorylase kinase. Eur J Biochem. 1980, 107 (2): 519-527.

    Article  CAS  PubMed  Google Scholar 

  9. Rankin CA, Sun Q, Gamblin TC: Tau phosphorylation by GSK-3beta promotes tangle-like filament morphology. Mol Neurodegener. 2007, 2: 12-10.1186/1750-1326-2-12.

    Article  PubMed Central  PubMed  Google Scholar 

  10. Kwok JB, Hallupp M, Loy CT, Chan DK, Woo J, Mellick GD, Buchanan DD, Silburn PA, Halliday GM, Schofield PR: GSK3B polymorphisms alter transcription and splicing in Parkinson's disease. Ann Neurol. 2005, 58 (6): 829-839. 10.1002/ana.20691.

    Article  CAS  PubMed  Google Scholar 

  11. Wider C, Vilarino-Guell C, Heckman MG, Jasinska-Myga B, Ortolaza-Soto AI, Diehl NN, Crook JE, Cobb SA, Bacon JA, Aasly JO, Gibson JM, Lynch T, Uitti RJ, Wszolek ZK, Farrer MJ, Ross OA: SNCA, MAPT, and GSK3B in Parkinson disease: a gene-gene interaction study. Eur J Neurol. 2011, 18 (6): 876-881. 10.1111/j.1468-1331.2010.03297.x.

    Article  CAS  PubMed  Google Scholar 

  12. Yuan Y, Tong Q, Zhou X, Zhang R, Qi Z, Zhang K: The association between glycogen synthase kinase 3 beta polymorphisms and Parkinson's disease susceptibility: a meta-analysis. Gene. 2013, 524 (2): 133-138. 10.1016/j.gene.2013.04.045.

    Article  CAS  PubMed  Google Scholar 

  13. Guedj M, Nuel G, Prum B: A note on allelic tests in case-control association studies. Ann Hum Genet. 2008, 72 (Pt 3): 407-409.

    Article  CAS  PubMed  Google Scholar 

  14. Faul F, Erdfelder E, Lang AG, Buchner A: G Power 3: a flexible statistical power analysis program for the social, behavioral, and biomedical sciences. Behav Res Methods. 2007, 39 (2): 175-191. 10.3758/BF03193146.

    Article  PubMed  Google Scholar 

  15. Evans W, Fung HC, Steele J, Eerola J, Tienari P, Pittman A, Silva R, Myers A, Vrieze FW, Singleton A, Hardy J: The tau H2 haplotype is almost exclusively Caucasian in origin. Neurosci Lett. 2004, 369 (3): 183-185. 10.1016/j.neulet.2004.05.119.

    Article  CAS  PubMed  Google Scholar 

  16. Zhang J, Song Y, Chen H, Fan D: The tau gene haplotype h1 confers a susceptibility to Parkinson's disease. Eur Neurol. 2005, 53 (1): 15-21. 10.1159/000082956.

    Article  CAS  PubMed  Google Scholar 

  17. Healy DG, Abou-Sleiman PM, Lees AJ, Casas JP, Quinn N, Bhatia K, Hingorani AD, Wood NW: Tau gene and Parkinson's disease: a case-control study and meta-analysis. J Neurol Neurosurg Psychiatry. 2004, 75 (7): 962-965. 10.1136/jnnp.2003.026203.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  18. Landrigan PJ, Sonawane B, Butler RN, Trasande L, Callan R, Droller D: Early environmental origins of neurodegenerative disease in later life. Environ Health Perspect. 2005, 113 (9): 1230-1233. 10.1289/ehp.7571.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  19. Winkler S, Konig IR, Lohmann-Hedrich K, Vieregge P, Kostic V, Klein C: Role of ethnicity on the association of MAPT H1 haplotypes and subhaplotypes in Parkinson's disease. Eur J Hum Genet. 2007, 15 (11): 1163-1168. 10.1038/sj.ejhg.5201901.

    Article  CAS  PubMed  Google Scholar 

  20. Das G, Misra AK, Das SK, Ray K, Ray J: Role of tau kinases (CDK5R1 and GSK3B) in Parkinson's disease: A study from India. Neurobiol Aging. 2012, 33 (7): 1485-e9-15.

    Article  PubMed  Google Scholar 

  21. Garcia-Gorostiaga I, Sanchez-Juan P, Mateo I, Rodriguez-Rodriguez E, Sanchez-Quintana C, del Olmo SC, Vazquez-Higuera JL, Berciano J, Combarros O, Infante J: Glycogen synthase kinase-3 beta and tau genes interact in Parkinson's and Alzheimer's diseases. Ann Neurol. 2009, 65 (6): 759-761. author reply 761-752

    Article  PubMed  Google Scholar 

  22. Goodenough S, Schleusner D, Pietrzik C, Skutella T, Behl C: Glycogen synthase kinase 3beta links neuroprotection by 17beta-estradiol to key Alzheimer processes. Neuroscience. 2005, 132 (3): 581-589. 10.1016/j.neuroscience.2004.12.029.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgments

This work was supported by grants 30870866, 81071021 and 31171211 from the National Natural Science Foundation of China (to TW), grant 81200983 from the National Natural Science Foundation of China (to NX), grant 81100958 from the National Natural Science Foundation of China (to ZTZ), grant 81301082 from the National Natural Science Foundation of China (to JSH), grant 2012B09 from China Medical Foundation (to NX) and grant 0203201343 from Hubei Molecular Imaging Key Laboratory (to NX). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Author information

Authors and Affiliations

Authors

Corresponding authors

Correspondence to Nian Xiong or Tao Wang.

Additional information

Competing interests

We confirm that we have read the Journal’s position on issues involved in ethical publication and affirm that this study is consistent with those guidelines. None of the authors have any conflict of interest to disclose.

Authors’ contributions

LY, JH, DZ, KG, NX, XL, JX, YZ, PW, ZL, JW, TW contributed to the conception and design. LY, DZ, KG, JX, YW, PW, JH, ZZ took care of the PCR studies. LY, JH, DZ, KG, LL, XL, JX, YZ, ZL, DW, JW, NX, TW analyzed and interoperated the data. LY, DZ, XL, LL, YZ, ZL,DW, JW, NX, TW coordinated all the experiments and helped to draft the manuscript. All authors read, revised and approved the final manuscript.

Lan Yu, Jinsha Huang, Desheng Zhai contributed equally to this work.

Electronic supplementary material

Authors’ original submitted files for images

Rights and permissions

Open Access  This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.

The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.

To view a copy of this licence, visit https://creativecommons.org/licenses/by/4.0/.

The Creative Commons Public Domain Dedication waiver (https://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Yu, L., Huang, J., Zhai, D. et al. MAPT rs242562 and GSK3B rs334558 are associated with Parkinson’s Disease in central China. BMC Neurosci 15, 54 (2014). https://0-doi-org.brum.beds.ac.uk/10.1186/1471-2202-15-54

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/1471-2202-15-54

Keywords