Analysis of coding variants in the betacellulin gene in type 2 diabetes and insulin secretion in African American subjects
© Elbein et al; licensee BioMed Central Ltd. 2006
Received: 16 May 2006
Accepted: 25 July 2006
Published: 25 July 2006
Betacellulin is a member of the epidermal growth factor family, expressed at the highest levels predominantly in the pancreas and thought to be involved in islet neogenesis and regeneration. Nonsynonymous coding variants were reported to be associated with type 2 diabetes in African American subjects. We tested the hypotheses that these previously identified variants were associated with type 2 diabetes in African Americans ascertained in Arkansas and that they altered insulin secretion in glucose tolerant African American subjects.
We typed three variants, exon1 Cys7Gly (C7G), exon 2 Leu44Phe (L44F), and exon 4 Leu124Met (L124M), in 188 control subjects and 364 subjects with type 2 diabetes. We tested for altered insulin secretion in 107 subjects who had undergone intravenous glucose tolerance tests to assess insulin sensitivity and insulin secretion.
No variant was associated with type 2 diabetes, and no variant altered insulin secretion or insulin sensitivity. However, an effect on lipids was observed for all 3 variants, and variant L124M was associated with obesity measures.
We were unable to confirm a role for nonsynonymous variants of betacellulin in the propensity to type 2 diabetes or to impaired insulin secretion.
Type 2 diabetes (T2DM) has a substantial genetic component, but identification of susceptibility genes has been limited by the large number of loci and small effect size. Defective insulin action is widely accepted as one element in the progression of metabolic syndrome to impaired glucose tolerance and T2DM, but increasing data support a role for early β-cell dysfunction in the pathogenesis of type 2 diabetes [1–3]. Impaired β-cell function predicts future diabetes , and work from our laboratory  and others [6, 7] suggest that the ability of pancreatic β-cell to compensate for prevailing insulin sensitivity is highly heritable. Mutations in β-cell transcription factors suggest that control of both β-cell mass and insulin secretion may play a role in the genetic susceptibility to β-cell failure and T2DM .
Betacellulin (BTC) is one of the factors potentially controlling β-cell growth. Although the BTC gene is located on chromosome 4q13-q21, which is not a replicated region of linkage to T2DM, BTC is nonetheless a strong candidate for T2DM. BTC, which was isolated from insulinoma cells, is highly expressed in pancreas and intestine . BTC belongs to the epidermal growth factor (EGF) family and appears to act through the EGF receptor, although other receptors including ErbB-4 and perhaps a specific BTC receptor have been proposed . Several models suggest that BTC can act to induce neogenesis of β-cells with resultant improvement in glucose homeostasis. Recombinant human BTC improved glucose tolerance and increased the number of islet-like cell clusters in alloxan-treated mice, suggesting increased islet neogenesis from ductal cells . Activin A and BTC together significantly reduced plasma glucose, increased pancreatic β-cell mass, and increased islet insulin content in streptozotocin treated rats . Rats treated with BTC alone after 90% pancreatectomy also experienced increased β-cell mass, increased islet insulin content, and improved glucose. Finally, BTC induced differentiation of the pancreatic exocrine cell line AR42J into insulin secreting cells [12, 13].
Silver et al. demonstrated that BTC was expressed in 9–24 week human fetal pancreas . They established that the human BTC gene contains 6 exons spanning at least 40 kb, of which the first 5 exons are translated. Three nonsynonymous coding variants were identified: Cys7Gly (C7G) in exon 1, Leu44Phe (L44F) in exon 2, and Leu124Met (L124M) in exon 4. The C7G variant showed a 13% lower minor allele frequency in African American cases (frequency 0.32; n = 185) than in controls (frequency 0.45; n = 149; p = 0.0004). Although no other variant was significantly associated with T2DM, and no variant was associated with T2DM in Caucasians, several haplotype combinations, including haplotypes comprising alleles at C7G and L44F, showed significant associations in African American subjects .
Recently, Silver and colleagues  reported the association of an intron 4 variant of BTC with decreased risk of type 1 diabetes in a case control study of 100 Caucasian cases and 282 Caucasian controls. This study was replicated in 113 informative trios, and provided additional evidence for BTC as a candidate for diabetes risk.
Based on these prior data and the strength of BTC as a candidate for the inherited β-cell defect that contributes to T2DM, we tested the role of these 3 coding variants in African American subjects ascertained in Arkansas in two studies: a case-control study of diabetic individuals with a family history of diabetes and normoglycemic control individuals, and a second study of glucose tolerant individuals who had undergone detailed assessment of insulin sensitivity and insulin secretion. We sought to test the hypothesis that coding defects of the BTC gene predispose to T2DM in African Americans and result in defective insulin secretion in nondiabetic individuals at risk.
Demographics of African American Case Control Population
Age at visit (years)
42.6 ± 13.1
54.8 ± 12.6*
Age of diagnosis
42.6 ± 11.9
Body Mass Index (kg/m 2 )
30.2 ± 7.2
32.7 ± 7.5*
Diabetic nephropathy (no nephropathy/nephropathy/unknown)
Characteristics of African American Subjects Who Underwent Intravenous Glucose Tolerance Testing
Mean ± standard deviation
37.2 ± 9.4
Body mass index (kg/m2)
30.5 ± 6.4
Waist circumference (cm)
97.6 ± 19.9
0.919 ± 0.108
Systolic blood pressure
106.4 ± 10.7
Diastolic blood pressure
63.1 ± 11.4
Percent fat (DEXA)
32.1 ± 11.2
SI (× 10-4 min-1/μU/ml)
4.95 ± 9.17
0.0194 ± 0.082
5013 ± 3996
2590 ± 2608
Fasting glucose (mmol/l)
4.66 ± 0.56
Fasting insulin (pmol/l)
54.3 ± 36.8
Genotyping of the three betacellulin polymorphisms
Conditions and the Primers Sequences
FORWARD PRIMER 5' to 3'
REVERSE PRIMER 5' to 3'
SEQUENCE PRIMER 5' to 3'
Insulin sensitivity (SI) was calculated from FSIGT glucose and insulin data using the MinMod Millenium software [17, 18]. Insulin secretion was measured as the acute 2 min – 10 min insulin response to the initial glucose bolus (AIRg), and disposition index (DI) was calculated as SI *AIRg, a measure of the β-cell response to insulin sensitivity .
Association of the three variants with T2DM were tested primarily for allelic association using either chi-squared or Fisher's Exact test, as implemented in either the 2by2 program  or the HaploView 3.2 program . Haplotypes were constructed from 2 and 3 locus combinations using both HaploView 3.2  and Phase 2.1 . Global significance was judged from the default 100 permutations. Impact of the three variants on quantitative traits in glucose tolerant individuals who had undergone FSIGT was tested using mixed effects models, in which ln-transformed BMI and age were covariates and protocol type (tolbutamide or insulin), gender, and genotype were fixed factors. Traits AIRg, DI, SI, triglycerides, and cholesterol were all ln-transformed to normality prior to analysis. Marginal means were tested for significance and converted back to the linear scale for presentation of results. Previous studies have suggested that although tolbutamide and insulin modified FSIGT are not interchangeable, the results are strongly correlated and thus the statistical method described here will be valid . Because fewer than 5 individuals were homozygous for the uncommon allele for both L124M and L44P variants, we only examined carrier status for the quantitative traits for both variants (for example, L124/L124 vs L124/M124 + M124/M124).
Association of Betacellulin Coding Variants with Type 2 Diabetes in African American Subjects
19 (Exon 1)
130 (Exon 2)
370 (Exon 4)
Genotypic Counts for Betacellulin Variants
Haplotype Analysis of Betacellulin Coding Variants
We next examined the role of the three BTC coding variants in 107 individuals who had normal glucose tolerance tests and who had undergone detailed assessment of insulin sensitivity (SI) and insulin secretion using the FSIGT. Only L124M altered the acute insulin response to glucose (AIRg), with carrier status increasing the response from 481 mg/dl among L124 homozygotes to 661 mg/dl in M124 carriers (p = 0.036; 95% CI 390 mg/dl – 593 mg/dl for L124 homozygotes, 531 mg/dl to 822 mg/dl for M124 carriers). No other variant altered the AIRg, the ability of the β-cell to compensate for insulin resistance (disposition index, DI = SI *AIRg), or in 40 individuals, the maximal insulin secretory response to arginine (AIRmax). Likewise, genotype did not have any main effects on insulin sensitivity, SI.
Marginal Means for Traits Associated with Betacellulin Polymorphisms in Non-diabetic Individuals
1.343 (1.131, 1.556)
1.025 (0.881, 1.169)
0.865 (0.653, 1.078)
HDL cholesterol mmol/l
1.372 (1.229, 1.515)
1.280 (1.183, 1.377)
1.536 (1.392, 1.679)
0.953 (0.853, 1.066)
1.174 (0.853, 1.066)
4.675 (4.400, 4.951)
4.2551 (3.967, 4.544)
Considerable data support heritable defects in β-cell function leading to increased diabetes susceptibility, but to date only the glutamine to lysine variant at position 23 in the β-cell potassium channel gene (KCNJ11 E23K) is widely accepted . That variant is uncommon in African Americans. Variants that explain the high prevalence of T2DM in African Americans are largely undiscovered. Betacellulin is a strong candidate to cause β-cell dysfunction and T2DM given its apparent role as a β-cell growth factor, the presence of three nonsynonymous coding variants, and the previous association of the C7G variant with T2DM in African Americans . Of the coding variants examined in this study, C7G in exon 1 is in the signal peptide. The minor allele, G7 in humans, is the wild type for mouse, rat, and cow BTC genes. Hence, this position is not conserved across species, and a functional role for this variant is uncertain. The other two positions are conserved across mouse, rat, and human, and might alter protein structure and function. L124M in the transmembrane domain is a conservative substitution, but L44F could alter tertiary protein structure. However, neither variant was associated with T2DM in the earlier study , nor do we find evidence for an association with T2DM. We did observe a modest association of homozygosity for L124 with reduced insulin secretion relative to M124 carrier status, but this association would not be significant after Bonferonni correction. Hence, our data do not support a role for BTC variants in either T2DM susceptibility or insulin secretion. BTC polymorphisms, either when considered individually or in haplotype combinations, were not associated with T2DM. Indeed, we found no trend to significance that would support the results of Silver et al . One limitation of our study is that our controls are relatively young, and despite our attempt to screen for family history and glucose tolerance, some controls may still develop T2DM in the future, which would reduce the power of our analysis. Furthermore, our analysis of insulin secretion was based on only 107 individuals who had undergone detailed phenotypic analysis; hence, power in this group was potentially limiting. Nonetheless, only the common homozygous state for L124M showed a trend to an association with altered insulin secretion.
The frequency of the G7 allele was similar in our control population to that reported by Silver et al.  (0.498 in the present study, 0.45 by Silver). However, we found a considerably higher frequency of this allele among African American diabetic subjects than was reported by Silver et al. (0.477 in the present study, 0.32 in the study of Silver et al.). The frequency of the L44F variant was similar in both cases and controls in our two studies, whereas L124M was somewhat more common in the present study in both cases (0.27 in the current study, 0.22 in Silver et al.) and controls (0.29 in the current study, 0.24 in Silver et al.). These differences likely represent chance fluctuations in allele frequencies given the size of the two studies. Based on the allele frequencies reported by Silver et al, we would have 85% power to find an association in our population, or based on the minor allele frequency of 0.32, 70% power to detect an OR of 1.6. However, using minor allele frequencies from our study and an OR of 1.5, our power was 62% for C7G, 65% for L124M, and only 33% for L44F. Hence, we could easily have missed effects at any of these variants of the size of PPARG  or KCNJ11  reported in previous studies in Caucasians (OR ≈ 1.2).
Alternatively, the differences may represent discrepant allele frequencies based on geographic differences in African American population structure . Notably, the data from Silver et al. showed marked differences in allele frequency between Caucasian and African American populations for both C7G and L44F variants, but not for the L124M variant. Hence, differences in population stratification between Arkansas and Maryland populations might explain the different results of our studies. We have examined over 40 SNPs in candidate genes not specifically selected for large differences between African American and Caucasian frequencies, and found no significant differences between cases and controls. More recently, we examined 25 markers that show marked differences between Caucasian and African American allele frequencies , and similarly found no significant association when tested by permutation (10,000 permutations; p > 0.08). Whereas we cannot completely exclude population stratification without a larger marker set and a larger sample size, these data make this an unlikely explanation for our lack of an association. Alternatively, differences in the ascertainment of our case and control populations in Arkansas and Baltimore might have resulted in different environmental or genetic factors that interact with BTC polymorphisms to increase the risk of T2DM.
We did find a significant effect of C7G on triglycerides and HDL cholesterol. Marginal means support a gene dosage effect. Furthermore, individuals homozygous for the common (C7) allele had the highest triglyceride levels, which is consistent with the association with T2DM found by Silver et al. We found a trend towards an association for L44F and triglycerides, with carrier status for the minor allele showing a trend towards higher triglyceride levels. We also found evidence for an association of reduced cholesterol in carriers of the L124 allele. These findings were not based on a priori hypotheses, and considering that we performed 8 analyses for 3 SNPs, none of these findings would reach statistical significance at a nominal p < 0.05 were a Bonferroni correction applied for 24 tests. Furthermore, a biological explanation for these findings is not obvious. BTC is expressed in the gut, but a role in lipid metabolism has not been reported. Hence, we cannot exclude a type II error as the explanation for these findings. This association will need confirmation in other sample sets given these limitations.
We were unable to find a role for coding variants of the BTC gene in either T2DM risk in African Americans or in altered β-cell function in glucose tolerant, African American individuals. This finding contrasts with the earlier report of an association of the C7G variant with T2DM in African American subjects, but is consistent with the apparently non-critical nature of this residue in the signal peptide and the lack of association observed in Caucasian individuals in the earlier study . Based on an allelic association using the case and control frequencies from Silver et al. and the sample size of the current study, we had 98% power to reject the null hypothesis of no difference in allele frequencies. Unidentified gene-gene or gene-environment interactions, differences in our ascertainment methods, or differences in African American populations in different geographical regions may account for our different findings.
- AIRG :
acute insulin response to glucose
body mass index
- SI :
insulin sensitivity index
type 2 diabetes
This work was supported by grants DK39311 and DK54636 from the National Institutes of Health/NIDDK, by the Research Service of the Department of Veterans Affairs, and by grant support for subject ascertainment from the American Diabetes Association. Subject ascertainment, DNA preparation, data management and statistical assistance were supported in part by grant M01RR14288 from National Institutes of Health/National Center for Research Resources to the General Clinical Research Center (GCRC) of the University of Arkansas for Medical Sciences, College of Medicine. We thank the GCRC nursing staff for the support of intravenous glucose tolerance tests, the GCRC laboratory for insulin assays and assistance with DNA extraction, Richard Harris for assistance with data management, Judith Johnson Cooper and Terri Hale for assistance with subject ascertainment.
- DeFronzo RA: Pathogenesis of type 2 diabetes: metabolic and molecular implications for identifying diabetes genes. Diabetes Reviews. 1997, 5: 177-270.Google Scholar
- Kahn SE: The relative contributions of insulin resistance and beta-cell dysfunction to the pathophysiology of Type 2 diabetes. Diabetologia. 2003, 46: 3-19. 10.1007/s00125-003-1190-9.View ArticlePubMedGoogle Scholar
- Ferrannini E, Mari A: Beta cell function and its relation to insulin action in humans: a critical appraisal. Diabetologia. 2004, 47: 943-956. 10.1007/s00125-004-1381-z.View ArticlePubMedGoogle Scholar
- Weyer C, Bogardus C, Mott DM, Pratley RE: The natural history of insulin secretory dysfunction and insulin resistance in the pathogenesis of type 2 diabetes mellitus. J Clin Invest. 1999, 104: 787-794.View ArticlePubMedPubMed CentralGoogle Scholar
- Elbein SC, Hasstedt SJ, Wegner K, Kahn SE: Heritability of pancreatic beta-cell function among nondiabetic members of Caucasian familial type 2 diabetic kindreds. J Clin Endocrinol Metab. 1999, 84: 1398-1403. 10.1210/jc.84.4.1398.PubMedGoogle Scholar
- Watanabe RM, Valle T, Hauser ER, Ghosh S, Eriksson J, Kohtamaki K, Ehnholm C, Tuomilehto J, Collins FS, Bergman RN, Boehnke M: Familiarity of quantitative metabolic traits in Finnish families with non-insulin-dependent diabetes mellitus. Finland-United States Investigation of NIDDM Genetics (FUSION) Study investigators. Hum Hered. 1999, 49: 159-168. 10.1159/000022865.View ArticlePubMedGoogle Scholar
- Poulsen P, Levin K, Petersen I, Christensen K, Beck-Nielsen H, Vaag A: Heritability of insulin secretion, peripheral and hepatic insulin action, and intracellular glucose partitioning in young and old danish twins. Diabetes. 2005, 54: 275-283.View ArticlePubMedGoogle Scholar
- Polonsky KS, Sturis J, Bell GI: Non-insulin-dependent diabetes mellitus – a genetically programmed failure of the beta cell to compensate for insulin resistance. N Engl J Med. 1996, 334: 777-783. 10.1056/NEJM199603213341207.View ArticlePubMedGoogle Scholar
- Yamamoto K, Miyagawa J, Waguri M, Sasada R, Igarashi K, Li M, Nammo T, Moriwaki M, Imagawa A, Yamagata K, Nakajima H, Namba M, Tochino Y, Hanafusa T, Matsuzawa Y: Recombinant human betacellulin promotes the neogenesis of beta-cells and ameliorates glucose intolerance in mice with diabetes induced by selective alloxan perfusion. Diabetes. 2000, 49: 2021-2027.View ArticlePubMedGoogle Scholar
- Ishiyama N, Kanzaki M, Seno M, Yamada H, Kobayashi I, Kojima I: Studies on the betacellulin receptor in pancreatic AR42J cells. Diabetologia. 1998, 41: 623-628. 10.1007/s001250050959.View ArticlePubMedGoogle Scholar
- Li L, Yi Z, Seno M, Kojima I: Activin A and betacellulin: effect on regeneration of pancreatic beta-cells in neonatal streptozotocin-treated rats. Diabetes. 2004, 53: 608-615.View ArticlePubMedGoogle Scholar
- Miyagawa J, Hanafusa O, Sasada R, Yamamoto K, Igarashi K, Yamamori K, Seno M, Tada H, Nammo T, Li M, Yamagata K, Nakajima H, Namba M, Kuwajima M, Matsuzawa Y: Immunohistochemical localization of betacellulin, a new member of the EGF family, in normal human pancreas and islet tumor cells. Endocr J. 1999, 46: 755-764.View ArticlePubMedGoogle Scholar
- Mashima H, Ohnishi H, Wakabayashi K, Mine T, Miyagawa J, Hanafusa T, Seno M, Yamada H, Kojima I: Betacellulin and activin A coordinately convert amylase-secreting pancreatic AR42J cells into insulin-secreting cells. J Clin Invest. 1996, 97: 1647-1654.View ArticlePubMedPubMed CentralGoogle Scholar
- Silver K, Tolea M, Wang J, Pollin TI, Yao F, Mitchell BD: The exon 1 Cys7Gly polymorphism within the betacellulin gene is associated with type 2 diabetes in African Americans. Diabetes. 2005, 54: 1179-1184.View ArticlePubMedGoogle Scholar
- Silver KD, Magnuson VL, Tolea M, Wang J, Hagopian WA, Mitchell BD: Association of a polymorphism in the betacellulin gene with type 1 diabetes mellitus in two populations. J Mol Med. 2006Google Scholar
- Ward WK, Bolgiano DC, McKnight B, Halter JB, Porte D: Diminished B cell secretory capacity in patients with noninsulin-dependent diabetes mellitus. J Clin Invest. 1984, 74: 1318-1328.View ArticlePubMedPubMed CentralGoogle Scholar
- Boston RC, Stefanovski D, Moate PJ, Sumner AE, Watanabe RM, Bergman RN: MINMOD Millennium: a computer program to calculate glucose effectiveness and insulin sensitivity from the frequently sampled intravenous glucose tolerance test. Diabetes Technol Ther. 2003, 5: 1003-1015. 10.1089/152091503322641060.View ArticlePubMedGoogle Scholar
- Pacini G, Bergman RN: MINMOD: a computer program to calculate insulin sensitivity and pancreatic responsivity from the frequently sampled intravenous glucose tolerance test. Comput Methods Programs Biomed. 1986, 23: 113-122. 10.1016/0169-2607(86)90106-9.View ArticlePubMedGoogle Scholar
- Kahn SE, Prigeon RL, McCulloch DK, Boyko EJ, Bergman RN, Schwartz MW, Neifing JL, Ward WK, Beard JC, Palmer JP, Porte D: Quantification of the relationship between insulin sensitivity and beta-cell function in human subjects: evidence for a hyperbolic function. Diabetes. 1993, 42: 1663-1672.View ArticlePubMedGoogle Scholar
- Ott J: Analysis of Human Genetic Linkage. 1999, Johns Hopkins University Press, Baltimore and London, ThirdGoogle Scholar
- Barrett JC, Fry B, Maller J, Daly MJ: Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics. 2005, 21: 263-265. 10.1093/bioinformatics/bth457.View ArticlePubMedGoogle Scholar
- Stephens M, Smith NJ, Donnelly P: A new statistical method for haplotype reconstruction from population data. Am J Hum Genet. 2001, 68: 978-989. 10.1086/319501.View ArticlePubMedPubMed CentralGoogle Scholar
- Saad MF, Steil GM, Kades WW, Ayad MF, Elsewafy WA, Boyadjian R, Jinagouda SD, Bergman RN: Differences between the tolbutamide-boosted and the insulin-modified minimal model protocols. Diabetes. 1997, 46: 1167-1171.View ArticlePubMedGoogle Scholar
- Gabriel SB, Schaffner SF, Nguyen H, Moore JM, Roy J, Blumenstiel B, Higgins J, DeFelice M, Lochner A, Faggart M, Liu-Cordero SN, Rotimi C, Adeyemo A, Cooper R, Ward R, Lander ES, Daly MJ, Altshuler D: The structure of haplotype blocks in the human genome. Science. 2002, 296: 2225-2229. 10.1126/science.1069424.View ArticlePubMedGoogle Scholar
- Gloyn AL, Weedon MN, Owen KR, Turner MJ, Knight BA, Hitman G, Walker M, Levy JC, Sampson M, Halford S, McCarthy MI, Hattersley AT, Frayling TM: Large-scale association studies of variants in genes encoding the pancreatic beta-cell KATP channel subunits Kir6.2 (KCNJ11) and SUR1 (ABCC8) confirm that the KCNJ11 E23K variant is associated with type 2 diabetes. Diabetes. 2003, 52: 568-572.View ArticlePubMedGoogle Scholar
- Altshuler D, Hirschhorn JN, Klannemark M, Lindgren CM, Vohl MC, Nemesh J, Lane CR, Schaffner SF, Bolk S, Brewer C, Tuomi T, Gaudet D, Hudson TJ, Daly M, Groop L, Lander ES: The common PPARgamma Pro12Ala polymorphism is associated with decreased risk of type 2 diabetes. Nat Genet. 2000, 26: 76-80. 10.1038/79839.View ArticlePubMedGoogle Scholar
- Florez JC, Burtt N, de Bakker PI, Almgren P, Tuomi T, Holmkvist J, Gaudet D, Hudson TJ, Schaffner SF, Daly MJ, Hirschhorn JN, Groop L, Altshuler D: Haplotype structure and genotype-phenotype correlations of the sulfonylurea receptor and the islet ATP-sensitive potassium channel gene region. Diabetes. 2004, 53: 1360-1368.View ArticlePubMedGoogle Scholar
- Tang H, Quertermous T, Rodriguez B, Kardia SL, Zhu X, Brown A, Pankow JS, Province MA, Hunt SC, Boerwinkle E, Schork NJ, Risch NJ: Genetic structure, self-identified race/ethnicity, and confounding in case-control association studies. Am J Hum Genet. 2005, 76: 268-275. 10.1086/427888.View ArticlePubMedGoogle Scholar
- Smith MW, Patterson N, Lautenberger JA, Truelove AL, McDonald GJ, Waliszewska A, Kessing BD, Malasky MJ, Scafe C, Le E, De Jager PL, Mignault AA, Yi Z, De The G, Essex M, Sankale JL, Moore JH, Poku K, Phair JP, Goedert JJ, Vlahov D, Williams SM, Tishkoff SA, Winkler CA, De La Vega FM, Woodage T, Sninsky JJ, Hafler DA, Altshuler D, Gilbert DA, O'Brien SJ, Reich D: A high-density admixture map for disease gene discovery in african americans. Am J Hum Genet. 2004, 74: 1001-1013. 10.1086/420856.View ArticlePubMedPubMed CentralGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2350/7/62/prepub
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.