- Research article
- Open Access
Genetic evaluation of Addison's disease in the Portuguese Water Dog
© Oberbauer et al; licensee BioMed Central Ltd. 2006
- Received: 19 January 2006
- Accepted: 02 May 2006
- Published: 02 May 2006
Addison's disease, also known as hypoadrenocorticism, has been reported in many individual dogs, although some breeds exhibit a greater incidence than the population as a whole. Addison's is presumed to be an autoimmune mediated hereditary defect but the mode of inheritance remains unclear. In particular, the heritability and mode of inheritance have not been defined for the Portuguese Water Dog although Addison's is known to be prevalent in the breed.
The analyses present clear evidence that establishes Addison's disease as an inherited disorder in the Portuguese Water Dog with an estimate of heritability of 0.49 (± 0.16); there were no differences in risk for disease across sexes (p > 0.49). Further, the complex segregation analysis provides suggestive evidence that Addison's disease in the Portuguese Water Dog is inherited under the control of a single, autosomal recessive locus.
The high heritability and mode of inheritance of Addison's disease in the Portuguese Water Dog should enable the detection of segregating markers in a genome-wide scan and the identification of a locus linked to Addison's. Though the confirmation of Addison's disease as an autosomal recessive disorder must wait until the gene is identified, breeders of these dogs may wish to keep the present findings in mind as they plan their breeding programs to select against producing affected dogs.
- Major Locus
- ACTH Stimulation
- High Density Region
- ACTH Stimulation Test
- Standard Poodle
Addison's disease, also referred to as chronic adrenal insufficiency or hypoadrenocorticism, has been reported in many individual dogs. Although mixed breeds present as the "breed" most frequently diagnosed with Addison's, certain pure-breds, including the Great Dane, Rottweiler, Standard Poodle, Portuguese Water Dog, West Highland White Terrier, and Wheaton Terrier, experience a disproportionate level of prevalence of the disease [1, 2]. The defining characteristic in Addison's disease is the failure of the adrenal cortex to produce sufficient quantities of steroid hormones. Addison's can be characterized as either primary or secondary depending upon the site of impairment, with primary being most common in dogs . The symptoms associated with Addison's are often diffuse reflecting the general metabolic influence of the hormones produced by the adrenal cortex. Symptoms include lethargy, vomiting, inappetence, depression, and weight loss, [2–4].
Addison's disease in humans has long been characterized as an inherited trait [5, 6] induced by autoimmune-mediated destruction of the adrenal cortex. Similar to the human, the majority of dogs present with primary Addison's disease in which the adrenal cortex exhibits severe deterioration mediated by an autoimmune reaction . An inherited susceptibility for Addison's has been proposed for more than 30 dog breeds . Yet even though canine Addison's is well documented to run in families [9–12], the mode of inheritance is unclear. One breed, the Standard Poodle, exhibits the presence of a single major autosomal locus affecting expression , while another breed, the Bearded Collie, has a less definitive mode of inheritance for Addison's . We undertook genetic analysis to characterize the inheritance of Addison's in a third breed having a low genetic relationship to the Standard Poodle  to clarify whether unrelated breeds demonstrate unique modes of inheritance for Addison's as has been demonstrated for other disorders such as blood clotting [14–16]. In particular, the heritability and mode of inheritance have not been defined for the Portuguese Water Dog although breeders recognize that Addison's is a concern in the breed (Canine Health Foundation, American Kennel Club 2001 survey). Knowledge of a genetic influence on the expression of Addison's in the Portuguese Water Dog would assist breeders in making breeding decisions that would decrease the incidence of Addison's in their breed. In this study, the heritability and the mode of inheritance, as determined by complex segregation analysis, are reported for the Portuguese Water Dog.
Estimate of heritability and gender contrast in a threshold model for Addison's disease.
Females – Males
Complex segregation analysis
Marginal posterior means, modes, standard deviations and limits to the 95% highest density regions (HDR) of model parameters for Addison's disease in Portuguese Water dogs in a Bayesian mixed-inheritance model with a completely recessive major locus.
Major Locus Variance
Additive Effect (a)
Dominance Deviation (d)
HDR 95% Low
HDR 95% High
Of the total dogs in this study, 590 dogs were inbred; their average inbreeding coefficient was 0.061. There were as many as seven overlapping generations included in the families under study and not all dogs within the study had equivalent number of generations evaluated. Nevertheless, with an eye toward making simple recommendations to breeders, the relationship between inbreeding and the risk for Addison's disease should be evaluated. While the data set does not permit computation of the level of inbreeding of Portuguese Water Dogs as a whole, the relationship between Addison's and inbreeding can be assessed. Classifying affected and unaffected dogs by their inbreeding coefficient can reveal a simple relationship between level of inbreeding and the presence of disease. In the present study, among dogs with an inbreeding coefficient between zero (0.0) and 0.05, 11 of 205 dogs (5.4%) were recorded as affected by Addison's. However, as the level of inbreeding increased, the risk of Addison's disease increased with it. Specifically, for Portuguese Water Dogs with an inbreeding coefficient between 0.05 and 0.10 there were 8 affected dogs among a total of 97 dogs (8.3%). For those dogs with an inbreeding coefficient between 0.10 and 0.15 there were 15 of 70 (21.4%) dogs affected by Addison's. Though a small sample, among dogs with an inbreeding coefficient greater than 0.15 there were 5 affected dogs among a total of 23 (21.7%). This result is consistent with the results of Table 2 in that inbreeding increases homozygosity, and a disease that is strongly influenced by a segregating recessive locus should see an increased incidence among a cluster of inbred animals.
The objective of this study was to assess the genetic contribution to the expression of Addison's disease in the Portuguese Water Dog. The results provide sufficient and convincing evidence that Addison's disease in the Portuguese Water Dog is an inherited disorder. The mode of inheritance, determined in the present data set, is consistent with a single, autosomal recessive locus. This finding should be helpful to breeders as they struggle to reduce the prevalence of this disorder in the breed. Unfortunately, given the relatively late onset of illness, the evaluation of dogs at risk for expressing Addison's disease, or identification of carriers, is not straightforward. However, the evidence for the contributions of a single locus (Table 2) suggests that a genome scan for the putative disease-causing allele is advisable. Identification of the allelic mutation would enable the identification of dogs at risk or carriers prior to their contribution to the gene pool.
Though not presented in Table 2 for the sake of brevity, several models were fit in the complex segregation analysis. Specifically, fitting these models are intended to reduce the probability of a false detection of a major locus ; models in which the assumption of Mendelian transmission of the alleles of the putative major locus are relaxed. The null model is Mendelian transmission, and this model should be a significant improvement for the likelihood over the no-major locus model if we are to draw the conclusion that a major locus is segregating in this population of dogs. "Improvement", when working in the Bayesian framework provided by iBay , is reliant upon the examination of the 95% HDR. In our examination of Addison's disease in the Portuguese Water dog, allowing for non-Mendelian transmission of the putative major alleles did not improve the fit: the 95% HDR for the transmission probabilities demonstrated considerable overlap with the Mendelian transmission values. For this reason we accept the conclusion of the simplest null model that of the major locus model with alleles moving from parent to offspring with the typically expected probabilities.
Though complex segregation analysis is fraught with challenging statistical issues, the estimate of the allele frequency for the putative disease allele in the base population is 0.49 (i.e., frequency of allele "A" = .51, with frequency of the disease-causing "B" allele = 1. – .51 = .49; see Table 2). Accordingly, for our suggested major locus model one would predict that 24% of dogs would be affected with Addison's disease in a randomly mating population. The observed incidence in the data set is considerably lower with 103 affected dogs out of 804 (103/804 = 0.128), suggesting that breeders may be selecting parents so as to reduce the risk of disease. Such a high expected prevalence of Addison's greatly increases the power of identifying linkage by facilitating the detection of segregating markers in a genome-wide scan. In an effort to simply assess whether the Portuguese Water Dog family in the present study offers the power to detect genetic linkage with the Addison's phenotype, we simulated a recessive locus perfectly linked with Addison's in a subset of the larger Portuguese Water Dog family and computed expected LOD scores in a linkage analysis. The predicted LOD score of a marker within 5 cM of the Addisonian locus in a small subset family of 142 dogs was 16.55 using the CRIMAP analysis program  giving support to the concept that the family under study offers the power to detect genetic linkage with the Addison's phenotype.
As in a similar analysis of the Standard Poodle , there is no evidence for a sex difference in risk of Addison's disease for Portuguese Water Dogs. This suggests that the literature reports of a greater prevalence among females may represent an artifact of the sample populations evaluated: dogs referred to teaching hospitals rather than dogs evaluated by all veterinarians [1, 12]. The results of the present study also are consistent with the published evidence of the mode of inheritance for Addison's disease in the Standard Poodle . In both the Portuguese Water Dog and the Standard Poodle, the data support the view that the disease is inherited as a single autosomal recessive trait. In contrast, although highly heritable, the mode of inheritance of Addison's disease in the Bearded Collie is less clear . This disparity in the mode of inheritance predictions possibly stems from the variability in the onset of the disorder and the classification of dogs as "unaffected" in the statistical analyses when in fact the dogs have not yet expressed the disease. The consistency of high heritability for the Addisonian phenotype across these different breeds supports the possibility that Addison's may have a common mode of inheritance for dogs in general. This concept is supported given the high incidence among mixed breeds  and the recent determination that large haplotype regions are shared across multiple dog breeds suggesting that genetic liability predisposition may likewise be shared .
Until early clinical detection of Addisonian status exists or genetic tests are available to identify dogs at risk for developing Addison's disease or those dogs that are carriers of the mutant allele, breeders of Portuguese Water Dogs should consider the genetic contribution to the expression of Addison's disease. Although the present study cannot confirm the mode of inheritance, the putative mode of inheritance as an autosomal recessive disorder suggests that breeders could actively select against generating affected dogs in their breeding programs.
Portuguese Water Dogs have been noted to have a disproportionate number of Addisonian cases [1, 2]. To clarify the heritability and predict the mode of inheritance of Addison's disease in the Portuguese Water Dog, a study was initiated in conjunction with the Portuguese Water Dog Club of America. Interested owners of Portuguese Water Dogs submitted survey questionnaires which asked for registered name, sire and dam, date of birth, sex, whether intact or altered, whether Addisonian or not, how diagnosis was made, age of diagnosis, and whether the dog received steroid treatment prior to diagnosis. Pedigree information was obtained from the American Kennel Club records.
A dog was designated as Addisonian if the diagnosis was made following an adrenocorticotropin hormone (ACTH) stimulation test administered by a veterinarian. The ACTH test measures the competence of the adrenal cortex by evaluating blood cortisol concentrations before and after the exogenous ACTH stimulation. In a clinically normal animal, the baseline blood cortisol prior to the ACTH stimulation would be 0.5 to 4.0 μg of cortisol/dl that would rise to 8.0 to 20 μg/dl following ACTH stimulation. Dogs failing to respond to the ACTH stimulation with elevated circulating cortisol are classified as Addisonian . Most owners were prompted to test their Portuguese Water Dogs for Addison's disease because of signs of lethargy, vomiting, or collapse. The ACTH stimulation criteria for diagnosis allows for an accurate and repeatable assessment of the disease. Dogs included in the present study were classified as either Addison-affected (based on the above criteria) or unaffected Portuguese Water Dogs.
Addisonian disease incidence data were collected on 804 Portuguese Water Dogs from the United States, though with pedigree information an additional 1,273 animals were included in the study; these additional animals did not have a recorded observation for Addison's disease. The total of 2,077 dogs included in the analysis is derived from one large family (with 2,051 dogs), three families of seven dogs, and 5 dogs of unassigned pedigree. The only additional variable included in our analysis was sex. No additional phenotypic information was recorded along with the diagnosis of Addison's disease (e.g., coat color). Specifically, there were 458 females (396 unaffected and 62 affected) and 346 males (305 unaffected and 41 affected) with known health status. Figure 1, representing a small subset of the dogs in the study, is provided to demonstrate transmission of the disorder. However, the parents of all affected Portuguese Water Dogs in this study trace back to common ancestors. The Portuguese Water Dog population in the United States was established from two separate Portuguese kennel lines (Algarbiorum and Alvalade) beginning in the late 1960s .
For the affected dogs included in the analyses, the mean age of Addisonian diagnosis was 46 months (median 36 months). The mean age of the dogs in the study designated unaffected was 48.5 months (median 37 months). Caution must be exercised to avoid designating a dog as being Addisonian when in fact the responsiveness of the adrenal has been blunted by exogenous administration of glucocorticoids. For 6.4% of the Addisonian dogs in the present study, owners reported that the dogs had received some sort of short term steroid treatment at some point prior to the Addisonian diagnosis. However, in the United States, the standard veterinary protocol for an ACTH stimulation test is to gradually withdraw exogenous corticosteroids prior to administration of the ACTH challenge (R. Nelson, personal communication). In addition, those survey questionnaires that included the ACTH stimulation test values, Addisonian dogs had pre and post-ACTH cortisol levels of less than 0.2 μg/dl which differ from the expected values of a clinically normal dog as stated earlier. Taken together, this suggests that the diagnosis of Addison's disease for the dogs in this study were Addisonian and not phenocopies. For the dogs with alteration status recorded, the vast majority of dogs (88%) were diagnosed after being reproductively altered.
Estimation of heritability
For the objective of estimating the heritability of Addison's disease in the Portuguese Water Dog, a threshold model for the liability to disease was used. This method assumes that a dog can be assigned to a specific disease class (unaffected/affected) when an underlying, unobservable risk (or liability) for disease exceeds a threshold of τ = 0. The distribution of the unobservable liability was assumed to be multivariate normal. The correlation in liability of two dogs i and j was modeled to be ρij = aij h2 + δij σe2 where ρij is the correlation in liability to disease between dogs i and j; aij is the additive relationship between dogs i and j; h2 is the narrow sense heritability of liability to disease; δij is the coefficient for the random environmental component for dogs i and j such that δij equals 1 if i equals j and zero otherwise; and σe2 = 1 - h2, with no loss of generality. The null hypothesis of no genetic contribution (i.e., h2 = 0) was tested with a likelihood ratio test, comparing the full model likelihood with that of the restricted model where h2 = 0. Heritability is expressed as the mean ± standard error of the mean. It should be noted that because the data represent owner submissions, the data were collected in a non-random fashion. Further, being a study of inheritance, the data set was constructed around probands. Such data require an adjustment for ascertainment bias. However, the mixed linear models utilized in this study accommodate nonrandomly sampled data  as long as the dogs added into the study to complete the pedigree associations can be considered a random sample of Portuguese Water Dogs. In addition, a test of the effect of sex on the liability of Addison's disease was also tested through the likelihood ratio test. Calculations were implemented through the computer program SOLAR [25, 26], making use of the binary trait analysis first described in Duggirala et al. .
Complex segregation analysis
The possibility that Addison's disease in Portuguese Water dogs is influenced by the action of a single segregating locus of large effect can also be examined. Complex segregation analysis, developed by Bonney , is intended to integrate Mendelian transmission genetics at a single locus with the patterns of covariance expected in polygenic inheritance. Lynch and Walsh  provide a more complete description of complex segregation analysis. Elston et al.  outlined the criteria that must be satisfied before acceptance of the single major locus model so as to reduce the risk of false positive declarations of a major locus model. Evaluation of the models necessary for complex segregation analysis was conducted with the Bayesian software package iBay http://www.lucjanss.com. The iBay software is an extension of MaGGic , rewritten to accommodate complex segregation analysis in binary traits for pedigrees that include inbreeding.
The software selected to conduct the complex segregation analysis is built upon a Bayesian foundation, making use of Monte Carlo Markov chains (MCMC) through a Gibbs sampler. Accordingly, point estimates of unknown parameters are not derived, but rather estimates of the posterior density for unknown parameters. The iBay  package was recently used to evaluate the contribution of a major locus to osteochodral diseases in pigs , where a more complete outline of the MCMC approach is detailed. The goal of this strategy was to simultaneously estimate the posterior density for a polygenic contribution to Addison's disease along with the contributions of a putative Mendelian locus. Specifically, for this mixed-inheritance model, the strategy allowed the evaluation of a polygenic variance component, the additive and dominance contributions of a single locus (the parameters -a, d, and a for the putative major locus genotypes AA, AB, and BB, respectively) and the frequency of allele A of the putative major locus (defined as "q"). Given our scoring of phenotypes, where affected is 1 and unaffected scored as zero, the "B" allele represent the putative disease-enhancing allele. Note also that the iBay software models the unobservable scale of this threshold trait such that the residual variance is fixed at 1.0 (i.e., σe2 = 1).
Creation of the Gibbs sample requires several key assumptions about the behavior of these unknown parameters. Though a variety of models can be considered, all are some variant of the following: sex as a fixed effect with a flat (i.e., uniform) prior density, the polygenic variance component with a flat prior density, as well as flat prior densities for the additive, dominance, and allele frequency parameters. A Gibbs sample of 9,000 was generated, beginning with the creation of 350,000 total samples, a "burn-in" of 50,000 and a sampling rate of every 100-th Gibbs value. This process was repeated two additional times, to create three replicate chains. From the 9,000 Gibbs samples, the mean, standard deviation, mode and the upper and lower limits of a 95% highest density region (HDR) was computed for each of the unknown parameters. HDR were computed according to Hyndman  with the public domain software hdrcde , a package within the R-program .
The authors acknowledge the dedication of Portuguese Water Dog owners and breeders who supported this research. This work was supported by the American Kennel Club Canine Health Foundation grants Number 2226 and Number 0225. The authors would like to extend thanks to Dr. L. Janss for making the iBay software available, error-free and ready to run.
- Peterson ME, Kintzer PP, Kass PH: Pretreatment clinical and laboratory findings in dogs with hypoadrenocorticism: 225 cases (1979–1993). J Am Vet Med Assoc. 1996, 208: 85-91.PubMedGoogle Scholar
- Feldman EC, Nelson RW: Canine and Feline Endocrinology and Reproduction Third edition. St. Louis, MO: W. B. Saunders, 1089: 2004.Google Scholar
- Willard MD, Schall WD, McCaw DE, Nachreiner RF: Canine hypoadrenocorticism: report of 37 cases and review of 39 previously reported cases. J Am Vet Med Assoc. 1982, 180: 59-62.PubMedGoogle Scholar
- Rakich N, Lorenz M: Clinical signs and laboratory abnormalities in 23 dogs with spontaneous hypoadrenocorticism. J Am Anim Hosp Assoc. 1984, 20: 647-649.Google Scholar
- Maclaren NK, Riley WJ: Inherited susceptibility to autoimmune Addison's disease is linked to human leukocyte antigens-DR3 and/or DR4, except when associated with type I autoimmune polyglandular syndrome. J Clin Endocrinol Metab. 1986, 62: 455-459.View ArticlePubMedGoogle Scholar
- Vaidya B, Pearce S: The emerging role of the CTLA-4 gene in autoimmune endocrinopathies. Eur J Endocrinol. 2004, 150: 619-626. 10.1530/eje.0.1500619.View ArticlePubMedGoogle Scholar
- Schaer MS, Riley WJ, Buergelt CD, D.J B, Senior DF, Burrows CF, Campbell GA: Autoimmunity and Addison's disease in the dog. J Am Anim Hosp Assoc. 1986, 22: 789-794.Google Scholar
- Padgett GA: Control of canine genetic diseases. New York, NY: Howell Book House: 1998.Google Scholar
- Burton S, DeLay J, Holmes A, Somerville C, Eye J, Shaw D, Wack O, Hanna P: Hypoadrenocorticism in young related Nova Scotia duck tolling retrievers. Can Vet J. 1997, 38: 231-234.PubMed CentralPubMedGoogle Scholar
- Famula TR, Belanger JM, Oberbauer AM: Heritability and complex segregation analysis of hypoadrenocorticism in the standard poodle. J Small Anim Pract. 2003, 44: 8-12. 10.1111/j.1748-5827.2003.tb00096.x.View ArticlePubMedGoogle Scholar
- Oberbauer AM, Benemann KS, Belanger JM, Wagner DR, Ward JH, Famula TR: Inheritance of hypoadrenocorticism in bearded collies. Am J Vet Res. 2002, 63: 643-647. 10.2460/ajvr.2002.63.643.View ArticlePubMedGoogle Scholar
- Smallwood LJ, Barsanti JA: Hypoadrenocorticism in a family of leonbergers. J Am Anim Hosp Assoc. 1995, 31: 301-305.View ArticlePubMedGoogle Scholar
- Parker HG, Kim LV, Sutter NB, Carlson S, Lorentzen TD, Malek TB, Johnson GS, DeFrance HB, Ostrander EA, Kruglyak L: Genetic structure of the purebred domestic dog. Science. 2004, 304: 1160-1164. 10.1126/science.1097406.View ArticlePubMedGoogle Scholar
- Fogh JM, Fogh IT: Inherited coagulation disorders. Vet Clin North Am Small Anim Pract. 1988, 18: 231-243.View ArticlePubMedGoogle Scholar
- Brooks M, Raymond S, Catalfamo J: Severe, recessive von Willebrand's disease in German Wirehaired Pointers. J Am Vet Med Assoc. 1996, 209: 926-929.PubMedGoogle Scholar
- Riehl J, Okura M, Mignot E, Nishino S: Inheritance of von Willebrand's disease in a colony of Doberman Pinschers. Am J Vet Res. 2000, 61: 115-120. 10.2460/ajvr.2000.61.115.View ArticlePubMedGoogle Scholar
- Elston RC, Nasmboodrii KK, Glueck CJ, Fallat R, Tsang R, Leuba V: Studies of the genetic transmission of hypercholesterolemia and hypertriglyceridemia in 195 member kindred. Annals of Human Genetics. 1975, 39: 67-87.View ArticlePubMedGoogle Scholar
- Romatowski J: Iatrogenic adrenocortical insufficiency in dogs. J Am Vet Med Assoc. 1990, 196: 1144-1146.PubMedGoogle Scholar
- Janss LLG: iBay. 1.0- A Short Manual. 2006Google Scholar
- Green P: CRI-MAP. version 2.4 edition. St. Louis, MO: Washington University School of Medicine: 1990Google Scholar
- Lindblad-Toh K, Wade CM, Mikkelsen TS, Karlsson EK, Jaffe DB, Kamal M, Clamp M, Chang JL, Kulbokas EJ, Zody MC, et al: Genome sequence, comparative analysis and haplotype structure of the domestic dog. Nature. 2005, 438: 803-819. 10.1038/nature04338.View ArticlePubMedGoogle Scholar
- Kaufman J: Diseases of the adrenal cortex of dogs and cats. Mod Vet Pract. 1984, 65: 513-516.PubMedGoogle Scholar
- Braund K, Miller DF: The complete Portuguese water dog1 st Edition. New York, N.Y.: Howell Book House: 1986.Google Scholar
- Henderson CR: Applications of Linear Models in Animal Breeding. Guelph, ON Canada: University of Guelph: 1984.Google Scholar
- Almasy L, Blangero J: Multipoint quantitative-trait linkage analysis in general pedigrees. Am J Hum Genet. 1998, 62: 1198-1211. 10.1086/301844.PubMed CentralView ArticlePubMedGoogle Scholar
- Blangero J, Almasy L, Goring H, Williams J, Dyer T, Peterson C: Sequential oligogenic linkage analysis routines (SOLAR). 3.0.4. 2005Google Scholar
- Duggirala R, Williams JT, Williams-Blangero S, Blangero J: A variance component approach to dichotomous trait linkage analysis using a threshold model. Genet Epidemiol. 1997, 14: 987-992. 10.1002/(SICI)1098-2272(1997)14:6<987::AID-GEPI71>3.0.CO;2-G.View ArticlePubMedGoogle Scholar
- Bonney GE: Regressive logistic models for familial disease and other binary traits. Biometrics. 1986, 42: 611-625.View ArticlePubMedGoogle Scholar
- Lynch M, Walsh B: Genetics and analysis of quantitative traits. Sunderland, MA: Sinauer Associates: 1998.Google Scholar
- Janss LLG: MAGGIC: A package of subroutines for genetic analyses with Gibbs sampling. 6th World Congress on Genetics Applied to Livestock Production; Armidale, Australia. Organizing Committee on 6th WCGLAP; 1998:459-460.Google Scholar
- Kadarmideen HN, Janss LL: Evidence of a major gene from Bayesian segregation analyses of liability to osteochondral diseases in pigs. Genetics. 2005, 171: 1195-1206. 10.1534/genetics.105.040956.PubMed CentralView ArticlePubMedGoogle Scholar
- Hyndman R: Computing and graphing highest density regions. The American Statistician. 1996, 50: 120-126. 10.2307/2684423.Google Scholar
- Hyndman R: hdrcde: Highest density regions and conditional density estimation. R package version 1.01. 2005Google Scholar
- Team RDC: R: A language and environment for statistical computing. Vienna, Austria:2005.Google 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.