Get 20M+ Full-Text Papers For Less Than $1.50/day. Start a 14-Day Trial for You and Your Team.

Learn More →

Asian Origin for the Worldwide-Spread Mutational Event in Machado-Joseph Disease

Asian Origin for the Worldwide-Spread Mutational Event in Machado-Joseph Disease Abstract Background Machado-Joseph disease is the most frequent dominant ataxia worldwide. Despite its frequency and presence in many populations, only 2 founder mutations have been suggested to explain its current geographic distribution. Objectives To trace back in history the main mutational events in Machado-Joseph disease, we aimed to assess ancestral haplotypes and population backgrounds, to date the mutations, and to trace the routes and time of introduction of the founder haplotypes in different populations. Design, Setting, and Participants We studied 264 families with Machado-Joseph disease from 20 different populations. Six intragenic single-nucleotide polymorphisms were used to determine ancestral mutational events; 4 flanking short tandem repeats were used to construct extended haplotypes and measure accumulation of genetic diversity over time within each lineage. Results The worldwide-spread lineage, TTACAC, had its highest diversity in the Japanese population, where we identified the ancestral short tandem repeat–based haplotype. Accumulated variability suggested a postneolithic mutation, about 5774 ± 1116 years old, with more recent introductions in North America, Germany, France, Portugal, and Brazil. As to the second mutational event, in the GTGGCA lineage, only 7 families (of 71 families) did not have Portuguese ancestry, although gene diversity was again smaller in Portuguese families (0.44) than in non-Portuguese families (0.93). Conclusions The worldwide-spread mutation may have first occurred in Asia and later been diffused throughout Europe, with a founder effect accounting for its high prevalence in Portugal; the other Machado-Joseph disease lineage is more recent, about 1416 ± 434 years old, and its dispersion may be explained mainly by recent Portuguese emigration. Machado-Joseph disease (MJD) is an autosomal dominant ataxia first described in families of Portuguese-Azorean extraction (Machado, Thomas, and Joseph families) in the United States.1-3 Since then, MJD has been reported in many ethnic backgrounds and is known to be the most common spinocerebellar ataxia worldwide. Its relative frequency among the spinocerebellar ataxias is higher in Portugal (49%),4 China (49%),5 Brazil (44%6; 92% in Rio Grande do Sul, Brazil7), the Netherlands (44%),8 Japan (43%),9 and Germany (42%)10; its relative frequency is lower in France (33%),11 the United States (21%),12 and Australia (12%)13; and it is rare in the United Kingdom (5%),14 India (3%),15 and Italy (1%).16 Machado-Joseph disease is a neurodegenerative disorder characterized by ataxia and limitation of eye movements, combined with more severe spasticity and dystonia in patients with earlier onset and with peripheral neuropathy in patients with later onset.17 The gene responsible, the ataxin 3 gene (ATXN3), includes a CAG repeat in exon 10. Although exceptions have been reported, normal alleles range from 12 to 44 CAG repeats, whereas patients with MJD most often have 1 expanded allele with 61 to 87 repeats.18 In our previous study, only 4 MJD intragenic haplotypes were found, but 2 (ACA or GGC) were present in 94% of all of the families.19 Among patients originating from the Azores, we found island-specific haplotypes: ACA in Flores and GGC in São Miguel; both haplotypes were found in mainland Portugal. The ACA haplotype was shared by 72% of all of the families worldwide, suggesting a founder mutation, although its origin and diffusion remained unresolved. The hypothesis was raised that the original event occurred in mainland Portugal, spread to the archipelago during its colonization, and from there spread to Brazil, the United States, and Canada, while the Portuguese seafare could explain its presence in Japan, China, and India. Nevertheless, the possibility that the mutational event associated with ACA occurred in other populations and was brought to the Azores and mainland Portugal could not be excluded. Here, we intend to clarify these questions, determining the origins, age, and spread of the 2 mutational events through more extensive haplotype analyses. Methods Families with mjd and controls We studied 264 families with MJD, from Portugal (104 families), Brazil (37 families), France (29 families), Japan (27 families), North America (23 families), Germany (14 families), India (4 families), Australia (4 families), Spain (3 families), the Netherlands (2 families), Norway (2 families), Yemen (2 families), Taiwan (2 families), England (1 family), the West Indies (1 family), Belgium (1 family), Algeria (1 family), Somalia (1 family), the Ivory Coast (1 family), Cambodia (1 family), and unknown origin (4 families). Family origin was defined as the family's most remote known ancestry in the affected line (families classified here as North American had no known previous origin; among the Australian families, 3 claimed to have English, Scottish, and Chinese ancestry, but this was rather uncertain). We genotyped healthy individuals as controls from Europe (277 phase-known chromosomes from Portugal), Asia (100 unrelated chromosomes from China), and Africa (26 phase-known and 28 unrelated chromosomes from Angola and Mozambique). Peripheral blood samples were collected after informed consent was obtained, and genomic DNA was extracted by standard procedures. Identification of mjd lineages To analyze independent mutational events, we extended haplotypes, including (CAG)n, the previously studied single-nucleotide polymorphisms (SNPs) A669TG/G669TG (rs1048755), C987GG/G987GG (rs12895357), and TAA1118/TAC1118 (rs7158733), and the 3 additional SNPs IVS6-30G>T (rs12590497), GTT527/GTC527 (rs1130166), and C1178/A1178 (rs3092822) (Figure 1). Haplotype analyses Extended haplotypes were based on 4 flanking short tandem repeats (STRs), (TAT)n, (CA)n, (AC)n, and (GT)n (Figure 1). Amplification and genotyping conditions were performed as previously described.20 PHASE software version 2.0 (http://www.stat.washington.edu/stephens/software.html) was used to reconstruct haplotypes from genotypic data when the phase could not be directly inferred by family structure or by allele-specific amplification. Allele frequencies and phase-known haplotypes were taken into account, but only haplotype pairs with a probability greater than 0.6 were used for further analyses. To estimate genetic distances among populations, we performed pairwise analyses with Arlequin software version 2.000 (http://anthro.unige.ch/software/arlequin/software), using the sum of square size difference (RST) as a measure of distance; RST is an analogue of FST suited for STR haplotype comparisons because it explicitly considers the number of single mutation steps between alleles. The resulting matrix of interpopulation values was plotted with multidimensional scaling using Statistica version 6 software (StatSoft, Tulsa, Oklahoma). Estimation of age of mutations We estimated the age of the 2 mutations from the variation accumulated in their ancestral haplotypes. This method was modified to include both recombination (c) and mutation (μ) rates in the generation of variation.21,22 The probability of change per generation (ε) was given by ε = 1 − [(1 − c)(1 − μ)], and the average of mutation and recombination events (λ) equals εt, where t is the number of generations. Results Mjd lineages The TTACAC lineage was shared by 19 of 20 populations, whereas GTGGCA was observed only in Portugal and 7 non-Portuguese families, which were from North America (4 families), Spain (1 family), the West Indies (1 family), and unknown origin (1 family). In the Portuguese families, both haplotypes were frequent but their distribution was heterogeneous: all of the 12 families from the northeastern mainland and virtually all of the families from Flores shared the TTACAC haplotype; in São Miguel, 10 families had the GTGGCA haplotype (also predominant in the central mainland) and 3 had the GTGCCA haplotype (Figure 2). Populations of origin and spread of mjd mutations For TTACAC, we found the highest molecular diversity among Japanese families, followed (in this order) by North American, French, German, Portuguese, and Brazilian families (Table 1). This suggests that the worldwide spread of the MJD lineage is more ancient in Japan than in any other populations studied. Further introductions would have occurred in North America, then throughout Europe and later in Brazil, according to lower diversity among these populations. A less clear picture emerged for the origin of the GTGGCA lineage. The vast majority of families shared a Portuguese origin, harboring 7 of 11 STR-based haplotypes (Figure 3); however, haplotypes in non-Portuguese families were more distant from each other, showing higher diversity (gene diversity ± SD, 0.93 ± 0.12) when compared with Portuguese haplotypes, which were separated by a single mutation or recombination and displayed much less diversity (gene diversity ± SD, 0.44 ± 0.08). Ancestral str-based haplotypes The most parsimonious relationships among flanking haplotypes are presented as phylogenetic networks (Figure 3 and Figure 4). The probability of mutation vs recombination was evaluated, considering the number of stepwise mutations required, observed intermediate haplotypes, and allele and haplotype frequency in controls (eFigure). To determine ancestral or founder haplotypes, we considered the following: (1) frequency in each population, (2) number of unchanged alleles, and (3) number of families with haplotypes just 1 step away from that putative ancestral haplotype. According to these criteria, the mutation in TTACAC seems to have occurred in the H1 (8-25-14-19) haplotype (Table 2). Introduction in North American, French, Portuguese, and Brazilian populations occurred through the H12 (11-21-14-15) haplotype, whereas the 1-step-away H13 (10-21-14-15) haplotype was introduced in Germany. In the United States, the founder haplotype H12 reflects the European contribution, whereas phylogenetically more distant haplotypes suggest another, more ancient introduction of the same lineage; whether this occurred directly from the ancestral background or another genetic background is to be determined. Although there were no MJD haplotypes common to North American and Japanese families, intermediate haplotypes may have been lost over time. Assuming the TTACAC lineage spread worldwide from East Asia, it is intriguing that we could not find MJD haplotypes of Japanese origin as founders among North American families or in higher frequency in European families. Other populations, however (not included in this analysis owing to low numbers of families), carried some haplotypes from the ancestral Japanese background, such as those from India (H7 haplotype), the Ivory Coast (H7 haplotype), and Yemen (H3 haplotype), in addition to new haplotypes, such as the H28 (11-21-12-19) and H29 (10-23-18-19) haplotypes in India, the H30 (10-20-14-15) and H31 (13-23-14-15) haplotypes in Taiwan, the H32 (10-22-14-15) haplotype in Cambodia, and the H33 (16-21-14-13) haplotype in Belgium. It is noteworthy that Indian families, closer to the Japanese families in the multidimensional scaling representation (Figure 5), shared a complete haplotype and many STR alleles with the putative ancestral population. As for GTGGCA, the H′1 (10-28-18-19) haplotype seems to be the ancestral STR background, as it is shared by 65% of the families, whereas the rare haplotypes H′8, H′9, H′10, and H′11 could have originated and migrated from mainland Portugal to Spain, North America, and the West Indies. On the other hand, haplotypes H′3, H′4, and H′5, exclusively found in the Azores, derived from the second most frequent haplotype H′2 (16-22-18-19), a recombinant of H′1. In this population, the finding of a third SNP-based haplotype (GTGCCA), explainable by a single (recurrent) mutation in C987GG/G987GG, is also noteworthy. Age of mutational events Assuming a generation time of 25 years and based on the most accurate recombination estimates (from family data), the age calculated for the worldwide-spread lineage TTACAC was 5774 ± 1116 years (Table 2). Introduction of this mutation in North America seems to have occurred much later (3396 ± 1182 years ago), in concordance with the genetic diversity found; nevertheless, these ignore the possibility of different introductions of mutations identical by descent in a given population. In Europe, TTACAC seems to be even more recent, with ages less than 2000 years in Germany and France and less than 1000 years in Portugal. As for GTGGCA, the estimated age regardless of its origin was 1416 ± 434 years. Comment This study focuses on the main MJD mutations spread worldwide, explaining the current disease distribution. Analysis of haplotype backgrounds of fast-evolving polymorphisms (STRs) within distinct lineages (identified by genotyping of SNPs, which are rather stable) allowed us to assess their origins and follow their spreading routes. It is expected that, in the absence of strong selective pressures, the ancestral population of a given lineage will show the largest diversity, as mutation and recombination would have more time to generate new STR haplotypes. This should be taken with caution owing to the large variance that stochastic events underlying the evolutionary process carry: a large sample size is needed to infer origin, and inferences based on small samples may not be reliable. Moreover, accumulation of variation is also proportional to the effective population size; in smaller populations, haplotypes will coalesce earlier and thus accumulate less variation. Ttacac lineage The high haplotype diversity among Japanese families as compared with all of the other families with the TTACAC lineage supports its Asian origin. Regarding the high frequency of this SNP-based haplotype among Asian controls as compared with the European and African control populations,20 it could be hypothesized that recurrent mutations have been contributing to this high diversity. However, contrary to other trinucleotide repeat disorders, there are no reports to our knowledge of de novo MJD expansions, and intermediate alleles are rarely found. This is in accordance with a model where very few mutational events account for the cases observed nowadays. The higher levels of diversity among French and German families as compared with Portuguese families do not support the hypothesis of an introduction of the TTACAC lineage in Europe exclusively from Portugal. In Europe, the presence of TTACAC expanded chromosomes is recent and indicates that this has not spread along with the Neolithic culture. Moreover, the low diversity and close phylogenetic relationships among haplotypes of Portuguese and Brazilian families suggest that historical links between these 2 populations are responsible for the presence of MJD in Brazil. In North America, separate introductions of the TTACAC lineage may have occurred directly from the ancestral Asian population, Europe, or both. Gtggca lineage As for GTGGCA, results are more controversial. Although it could be argued that haplotype variability was underestimated owing to sampling effects, this seems unlikely if we consider the high number of families studied and their rather distinct origins; thus, our findings might still be compatible with a Portuguese origin. Moreover, GTGGCA reaches its higher frequency in controls in Europe.20 Taken at face value, the lower diversity among Portuguese haplotypes would imply a non-Portuguese origin for the ancestral line of GTGGCA; nevertheless, most of these chromosomes came from small isolated populations such as the Azores, which could lead to an underestimation of its age. One could contemplate the hypothesis of a non-Portuguese origin with a subsequent expansion within Portugal and loss of the ancestral haplotype; this is, however, a less parsimonious explanation. Age and spread of ancestral mutational events The ancient TTACAC lineage may have occurred during the Bronze Age; the GTGGCA lineage, regardless of its population of origin, may have occurred less than 2000 years BP. Thus, both MJD mutations seem to be fairly recent as compared with other expanded repeat disorders. In Friedreich ataxia, the origin of expanded alleles was dated at about 25 000 years ago, following the population expansion of the Upper Paleolithic23; for type 2 myotonic dystrophy, the founding haplotype was calculated to be 5000 to 12 500 years old24; in Huntington disease, the presence of the ancestral haplotype only in the Land of Valencia (Spain) was estimated to range from 4700 to 10 000 years old.25 The age estimated for the expansion of TTACAC could explain its diffusion westward from East Asia through the trade routes established in central Asia known as the Silk Road; however, these results do not support admixture between Indian and Portuguese MJD populations and are against the hypothesis of MJD in India having been introduced from Portugal19 as its reverse.15 Historical links with Portugal, on the other hand, would be responsible for MJD among Brazilian families. If the GTGGCA lineage indeed has a Portuguese origin, we could speculate about the Azores vs mainland Portugal as its starting point. Taking into account only the haplotype diversity in the Azores together with a recurrent mutation in C987GG/G987GG (a rare event) in that population, the H′1 haplotype could be a recombinant from the ancestral H′2, whereas the high number of mainland Portuguese families with the H′1 haplotype would be explained by a founder effect. Nevertheless, phylogenetic relationships among haplotypes from different populations suggest that more diversity has been generated from H′1, supporting its mutational origin in mainland Portugal. This fits well with the fact that the Azores were inhabited when they were first discovered in 1439, subsequently being colonized mostly from mainland Portugal. Conclusions The scarcity of de novo expansions in MJD led us to question worldwide disease frequencies. Considering the clinical phenomenon of anticipation together with the rarity of intermediate alleles that could act as a reservoir for recurrent full mutations, the extinction of expanded MJD alleles would be expected; however, we must take into account that a theoretical progressive elimination of alleles from the human gene pool must be faced over an evolutionary timescale where other demographic events interfere as well. Moreover, distinct mutational mechanisms may apply for different repeat disorders; it would be interesting to explore cis- or trans-acting factors involved in repeat stability in MJD. Back to top Article Information Correspondence: Jorge Sequeiros, MD, PhD, Unidade de Investigação Genética e Epidemiológica em Doenças Neurológicas, Instituto de Biologia Molecular e Celular, Rua Campo Alegre 823, 4150-180 Porto, Portugal (jsequeir@ibmc.up.pt). Accepted for Publication: November 28, 2006. Author Contributions:Study concept and design: Martins, Coutinho, Amorim, and Sequeiros. Acquisition of data: Martins, Gaspar, Wong, Silveira, Nicholson, Brunt, Tranebjaerg, Stevanin, Hsieh, Soong, Loureiro, Dürr, Tsuji, Watanabe, Jardim, Giunti, Riess, Ranum, Brice, Rouleau, and Coutinho. Analysis and interpretation of data: Martins, Calafell, Silveira, Amorim, and Sequeiros. Drafting of the manuscript: Martins, Amorim, and Sequeiros. Critical revision of the manuscript for important intellectual content: Calafell, Gaspar, Wong, Silveira, Nicholson, Brunt, Tranebjaerg, Stevanin, Hsieh, Soong, Loureiro, Dürr, Tsuji, Watanabe, Jardim, Giunti, Riess, Ranum, Brice, Rouleau, and Coutinho. Statistical analysis: Martins, Calafell, and Amorim. Obtained funding: Amorim and Sequeiros. Study supervision: Calafell, Coutinho, Amorim, and Sequeiros. Financial Disclosure: None reported. Funding/Support: This work was supported in part by Fundação para a Ciência e Tecnologia, through a research grant from Programa Operacional Ciência e Inovação (Instituto de Patologia e Imunologia Molecular da Universidade do Porto), and through Financiamento Plurianual de Unidades (Unidade de Investigação Genética e Epidemiológica em Doenças Neurológicas, Instituto de Biologia Molecular e Celular). Dr Martins was a recipient of scholarship SFRH/BD/8880/2002 from Fundação para a Ciência e Tecnologia. Additional Contributions: Albertino Damasceno, MD, and Benilde Soares, MD, of the Eduardo Mondlane University, Maputo, provided the Mozambican control samples. Nicholas Wood, MD, PhD, provided collaboration. Paula Magalhães, MSc, Susana Carrilho, MSc, and Joana Cerqueira, MSc, provided efficient technical assistance at Centro de Genética Preditiva e Preventiva, Instituto de Biologia Molecular e Celular. We are grateful to the patients and relatives who agreed to collaborate in this project. References 1. Nakano KKDawson DMSpence A Machado disease: a hereditary ataxia in Portuguese emigrants to Massachusetts. Neurology 1972;22 (1) 49- 55PubMedGoogle ScholarCrossref 2. Woods BTSchaumburg HH Nigro-spino-dentatal degeneration with nuclear ophthalmoplegia: a unique and partially treatable clinico-pathological entity. J Neurol Sci 1972;17 (2) 149- 166PubMedGoogle ScholarCrossref 3. Rosenberg RNNyhan WLBay CShore P Autosomal dominant striato-nigral degeneration: a clinical, pathologic and biochemical study of a new genetic disorder. Neurology 1976;26 (8) 703- 714PubMedGoogle ScholarCrossref 4. Silveira IMiranda CGuimarães L et al. Trinucleotide repeats in 202 families with ataxia: a small expanded (CAG)n allele at the SCA17 locus. Arch Neurol 2002;59 (4) 623- 629PubMedGoogle ScholarCrossref 5. Jiang HTang BXu B et al. Frequency analysis of autosomal dominant spinocerebellar ataxias in Han population in the Chinese mainland and clinical and molecular characterization of spinocerebellar ataxia type 6 [in Chinese]. Zhonghua Yi Xue Yi Chuan Xue Za Zhi 2005;22 (1) 1- 4PubMedGoogle Scholar 6. Lopes-Cendes ITeive HGCalcagnotto ME et al. Frequency of the different mutations causing spinocerebellar ataxia (SCA1, SCA2, MJD/SCA3 and DRPLA) in a large group of Brazilian patients. Arq Neuropsiquiatr 1997;55(3B)519- 529PubMedGoogle ScholarCrossref 7. Jardim LBSilveira IPereira ML et al. A survey of spinocerebellar ataxia in South Brazil: 66 new cases with Machado-Joseph disease, SCA7, SCA8, or unidentified disease-causing mutations. J Neurol 2001;248 (10) 870- 876PubMedGoogle ScholarCrossref 8. van de Warrenburg BPSinke RJVerschuuren-Bemelmans CC et al. Spinocerebellar ataxias in the Netherlands: prevalence and age at onset variance analysis. Neurology 2002;58 (5) 702- 708PubMedGoogle ScholarCrossref 9. Takano HCancel GIkeuchi T et al. Close associations between prevalences of dominantly inherited spinocerebellar ataxias with CAG-repeat expansions and frequencies of large normal CAG alleles in Japanese and Caucasian populations. Am J Hum Genet 1998;63 (4) 1060- 1066PubMedGoogle ScholarCrossref 10. Schöls LAmoiridis GBüttner TPrzuntek HEpplen JTRiess O Autosomal dominant cerebellar ataxia: phenotypic differences in genetically defined subtypes? Ann Neurol 1997;42 (6) 924- 932PubMedGoogle ScholarCrossref 11. Stevanin GDürr ADavid G et al. Clinical and molecular features of spinocerebellar ataxia type 6. Neurology 1997;49 (5) 1243- 1246PubMedGoogle ScholarCrossref 12. Moseley MLBenzow KASchut LJ et al. Incidence of dominant spinocerebellar and Friedreich triplet repeats among 361 ataxia families. Neurology 1998;51 (6) 1666- 1671PubMedGoogle ScholarCrossref 13. Storey Edu Sart DShaw JH et al. Frequency of spinocerebellar ataxia types 1, 2, 3, 6, and 7 in Australian patients with spinocerebellar ataxia. Am J Med Genet 2000;95 (4) 351- 357PubMedGoogle ScholarCrossref 14. Giunti PSabbadini GSweeney MG et al. The role of the SCA2 trinucleotide repeat expansion in 89 autosomal dominant cerebellar ataxia families: frequency, clinical and genetic correlates. Brain 1998;121459- 467PubMedGoogle ScholarCrossref 15. Mittal USrivastava AKJain SJain SMukerji M Founder haplotype for Machado-Joseph disease in the Indian population: novel insights from history and polymorphism studies. Arch Neurol 2005;62 (4) 637- 640PubMedGoogle ScholarCrossref 16. Brusco AGellera CCagnoli C et al. Molecular genetics of hereditary spinocerebellar ataxia: mutation analysis of spinocerebellar ataxia genes and CAG/CTG repeat expansion detection in 225 Italian families. Arch Neurol 2004;61 (5) 727- 733PubMedGoogle ScholarCrossref 17. Coutinho PAndrade C Autosomal dominant system degeneration in Portuguese families of the Azores Islands: a new genetic disorder involving cerebellar, pyramidal, extrapyramidal and spinal cord motor functions. Neurology 1978;28 (7) 703- 709PubMedGoogle ScholarCrossref 18. Maciel PCosta MCFerro A et al. Improvement in the molecular diagnosis of Machado-Joseph disease. Arch Neurol 2001;58 (11) 1821- 1827PubMedGoogle ScholarCrossref 19. Gaspar CLopes-Cendes IHayes S et al. Ancestral origins of the Machado-Joseph disease mutation: a worldwide haplotype study. Am J Hum Genet 2001;68 (2) 523- 528PubMedGoogle ScholarCrossref 20. Martins SCalafell FWong VCNSequeiros JAmorim A A multi-step mutation mechanism drives the evolution of the CAG repeat at MJD/SCA3 locus. Eur J Hum Genet 2006;14 (8) 932- 940PubMedGoogle ScholarCrossref 21. Morral NBertranpetit JEstivill X et al. The origin of the major cystic fibrosis mutation (delta F508) in European populations. Nat Genet 1994;7 (2) 169- 175PubMedGoogle ScholarCrossref 22. Cobo AMSaenz APoza JJ et al. A common haplotype associated with the Basque 2362AG→TCATCT mutation in the muscular calpain-3 gene. Hum Biol 2004;76 (5) 731- 741PubMedGoogle ScholarCrossref 23. Labuda MLabuda DMiranda C et al. Unique origin and specific ethnic distribution of the Friedreich ataxia GAA expansion. Neurology 2000;54 (12) 2322- 2324PubMedGoogle ScholarCrossref 24. Bachinski LLUdd BMeola G et al. Confirmation of the type 2 myotonic dystrophy (CCTG)n expansion mutation in patients with proximal myotonic myopathy/proximal myotonic dystrophy of different European origins: a single shared haplotype indicates an ancestral founder effect. Am J Hum Genet 2003;73 (4) 835- 848PubMedGoogle ScholarCrossref 25. García-Planells JBurguera JASolis P et al. Ancient origin of the CAG expansion causing Huntington disease in a Spanish population. Hum Mutat 2005;25 (5) 453- 459PubMedGoogle ScholarCrossref http://www.deepdyve.com/assets/images/DeepDyve-Logo-lg.png Archives of Neurology American Medical Association

Loading next page...
 
/lp/american-medical-association/asian-origin-for-the-worldwide-spread-mutational-event-in-machado-l6LCVjS8Dm
Publisher
American Medical Association
Copyright
Copyright © 2007 American Medical Association. All Rights Reserved.
ISSN
0003-9942
DOI
10.1001/archneur.64.10.1502
pmid
17923634
Publisher site
See Article on Publisher Site

Abstract

Abstract Background Machado-Joseph disease is the most frequent dominant ataxia worldwide. Despite its frequency and presence in many populations, only 2 founder mutations have been suggested to explain its current geographic distribution. Objectives To trace back in history the main mutational events in Machado-Joseph disease, we aimed to assess ancestral haplotypes and population backgrounds, to date the mutations, and to trace the routes and time of introduction of the founder haplotypes in different populations. Design, Setting, and Participants We studied 264 families with Machado-Joseph disease from 20 different populations. Six intragenic single-nucleotide polymorphisms were used to determine ancestral mutational events; 4 flanking short tandem repeats were used to construct extended haplotypes and measure accumulation of genetic diversity over time within each lineage. Results The worldwide-spread lineage, TTACAC, had its highest diversity in the Japanese population, where we identified the ancestral short tandem repeat–based haplotype. Accumulated variability suggested a postneolithic mutation, about 5774 ± 1116 years old, with more recent introductions in North America, Germany, France, Portugal, and Brazil. As to the second mutational event, in the GTGGCA lineage, only 7 families (of 71 families) did not have Portuguese ancestry, although gene diversity was again smaller in Portuguese families (0.44) than in non-Portuguese families (0.93). Conclusions The worldwide-spread mutation may have first occurred in Asia and later been diffused throughout Europe, with a founder effect accounting for its high prevalence in Portugal; the other Machado-Joseph disease lineage is more recent, about 1416 ± 434 years old, and its dispersion may be explained mainly by recent Portuguese emigration. Machado-Joseph disease (MJD) is an autosomal dominant ataxia first described in families of Portuguese-Azorean extraction (Machado, Thomas, and Joseph families) in the United States.1-3 Since then, MJD has been reported in many ethnic backgrounds and is known to be the most common spinocerebellar ataxia worldwide. Its relative frequency among the spinocerebellar ataxias is higher in Portugal (49%),4 China (49%),5 Brazil (44%6; 92% in Rio Grande do Sul, Brazil7), the Netherlands (44%),8 Japan (43%),9 and Germany (42%)10; its relative frequency is lower in France (33%),11 the United States (21%),12 and Australia (12%)13; and it is rare in the United Kingdom (5%),14 India (3%),15 and Italy (1%).16 Machado-Joseph disease is a neurodegenerative disorder characterized by ataxia and limitation of eye movements, combined with more severe spasticity and dystonia in patients with earlier onset and with peripheral neuropathy in patients with later onset.17 The gene responsible, the ataxin 3 gene (ATXN3), includes a CAG repeat in exon 10. Although exceptions have been reported, normal alleles range from 12 to 44 CAG repeats, whereas patients with MJD most often have 1 expanded allele with 61 to 87 repeats.18 In our previous study, only 4 MJD intragenic haplotypes were found, but 2 (ACA or GGC) were present in 94% of all of the families.19 Among patients originating from the Azores, we found island-specific haplotypes: ACA in Flores and GGC in São Miguel; both haplotypes were found in mainland Portugal. The ACA haplotype was shared by 72% of all of the families worldwide, suggesting a founder mutation, although its origin and diffusion remained unresolved. The hypothesis was raised that the original event occurred in mainland Portugal, spread to the archipelago during its colonization, and from there spread to Brazil, the United States, and Canada, while the Portuguese seafare could explain its presence in Japan, China, and India. Nevertheless, the possibility that the mutational event associated with ACA occurred in other populations and was brought to the Azores and mainland Portugal could not be excluded. Here, we intend to clarify these questions, determining the origins, age, and spread of the 2 mutational events through more extensive haplotype analyses. Methods Families with mjd and controls We studied 264 families with MJD, from Portugal (104 families), Brazil (37 families), France (29 families), Japan (27 families), North America (23 families), Germany (14 families), India (4 families), Australia (4 families), Spain (3 families), the Netherlands (2 families), Norway (2 families), Yemen (2 families), Taiwan (2 families), England (1 family), the West Indies (1 family), Belgium (1 family), Algeria (1 family), Somalia (1 family), the Ivory Coast (1 family), Cambodia (1 family), and unknown origin (4 families). Family origin was defined as the family's most remote known ancestry in the affected line (families classified here as North American had no known previous origin; among the Australian families, 3 claimed to have English, Scottish, and Chinese ancestry, but this was rather uncertain). We genotyped healthy individuals as controls from Europe (277 phase-known chromosomes from Portugal), Asia (100 unrelated chromosomes from China), and Africa (26 phase-known and 28 unrelated chromosomes from Angola and Mozambique). Peripheral blood samples were collected after informed consent was obtained, and genomic DNA was extracted by standard procedures. Identification of mjd lineages To analyze independent mutational events, we extended haplotypes, including (CAG)n, the previously studied single-nucleotide polymorphisms (SNPs) A669TG/G669TG (rs1048755), C987GG/G987GG (rs12895357), and TAA1118/TAC1118 (rs7158733), and the 3 additional SNPs IVS6-30G>T (rs12590497), GTT527/GTC527 (rs1130166), and C1178/A1178 (rs3092822) (Figure 1). Haplotype analyses Extended haplotypes were based on 4 flanking short tandem repeats (STRs), (TAT)n, (CA)n, (AC)n, and (GT)n (Figure 1). Amplification and genotyping conditions were performed as previously described.20 PHASE software version 2.0 (http://www.stat.washington.edu/stephens/software.html) was used to reconstruct haplotypes from genotypic data when the phase could not be directly inferred by family structure or by allele-specific amplification. Allele frequencies and phase-known haplotypes were taken into account, but only haplotype pairs with a probability greater than 0.6 were used for further analyses. To estimate genetic distances among populations, we performed pairwise analyses with Arlequin software version 2.000 (http://anthro.unige.ch/software/arlequin/software), using the sum of square size difference (RST) as a measure of distance; RST is an analogue of FST suited for STR haplotype comparisons because it explicitly considers the number of single mutation steps between alleles. The resulting matrix of interpopulation values was plotted with multidimensional scaling using Statistica version 6 software (StatSoft, Tulsa, Oklahoma). Estimation of age of mutations We estimated the age of the 2 mutations from the variation accumulated in their ancestral haplotypes. This method was modified to include both recombination (c) and mutation (μ) rates in the generation of variation.21,22 The probability of change per generation (ε) was given by ε = 1 − [(1 − c)(1 − μ)], and the average of mutation and recombination events (λ) equals εt, where t is the number of generations. Results Mjd lineages The TTACAC lineage was shared by 19 of 20 populations, whereas GTGGCA was observed only in Portugal and 7 non-Portuguese families, which were from North America (4 families), Spain (1 family), the West Indies (1 family), and unknown origin (1 family). In the Portuguese families, both haplotypes were frequent but their distribution was heterogeneous: all of the 12 families from the northeastern mainland and virtually all of the families from Flores shared the TTACAC haplotype; in São Miguel, 10 families had the GTGGCA haplotype (also predominant in the central mainland) and 3 had the GTGCCA haplotype (Figure 2). Populations of origin and spread of mjd mutations For TTACAC, we found the highest molecular diversity among Japanese families, followed (in this order) by North American, French, German, Portuguese, and Brazilian families (Table 1). This suggests that the worldwide spread of the MJD lineage is more ancient in Japan than in any other populations studied. Further introductions would have occurred in North America, then throughout Europe and later in Brazil, according to lower diversity among these populations. A less clear picture emerged for the origin of the GTGGCA lineage. The vast majority of families shared a Portuguese origin, harboring 7 of 11 STR-based haplotypes (Figure 3); however, haplotypes in non-Portuguese families were more distant from each other, showing higher diversity (gene diversity ± SD, 0.93 ± 0.12) when compared with Portuguese haplotypes, which were separated by a single mutation or recombination and displayed much less diversity (gene diversity ± SD, 0.44 ± 0.08). Ancestral str-based haplotypes The most parsimonious relationships among flanking haplotypes are presented as phylogenetic networks (Figure 3 and Figure 4). The probability of mutation vs recombination was evaluated, considering the number of stepwise mutations required, observed intermediate haplotypes, and allele and haplotype frequency in controls (eFigure). To determine ancestral or founder haplotypes, we considered the following: (1) frequency in each population, (2) number of unchanged alleles, and (3) number of families with haplotypes just 1 step away from that putative ancestral haplotype. According to these criteria, the mutation in TTACAC seems to have occurred in the H1 (8-25-14-19) haplotype (Table 2). Introduction in North American, French, Portuguese, and Brazilian populations occurred through the H12 (11-21-14-15) haplotype, whereas the 1-step-away H13 (10-21-14-15) haplotype was introduced in Germany. In the United States, the founder haplotype H12 reflects the European contribution, whereas phylogenetically more distant haplotypes suggest another, more ancient introduction of the same lineage; whether this occurred directly from the ancestral background or another genetic background is to be determined. Although there were no MJD haplotypes common to North American and Japanese families, intermediate haplotypes may have been lost over time. Assuming the TTACAC lineage spread worldwide from East Asia, it is intriguing that we could not find MJD haplotypes of Japanese origin as founders among North American families or in higher frequency in European families. Other populations, however (not included in this analysis owing to low numbers of families), carried some haplotypes from the ancestral Japanese background, such as those from India (H7 haplotype), the Ivory Coast (H7 haplotype), and Yemen (H3 haplotype), in addition to new haplotypes, such as the H28 (11-21-12-19) and H29 (10-23-18-19) haplotypes in India, the H30 (10-20-14-15) and H31 (13-23-14-15) haplotypes in Taiwan, the H32 (10-22-14-15) haplotype in Cambodia, and the H33 (16-21-14-13) haplotype in Belgium. It is noteworthy that Indian families, closer to the Japanese families in the multidimensional scaling representation (Figure 5), shared a complete haplotype and many STR alleles with the putative ancestral population. As for GTGGCA, the H′1 (10-28-18-19) haplotype seems to be the ancestral STR background, as it is shared by 65% of the families, whereas the rare haplotypes H′8, H′9, H′10, and H′11 could have originated and migrated from mainland Portugal to Spain, North America, and the West Indies. On the other hand, haplotypes H′3, H′4, and H′5, exclusively found in the Azores, derived from the second most frequent haplotype H′2 (16-22-18-19), a recombinant of H′1. In this population, the finding of a third SNP-based haplotype (GTGCCA), explainable by a single (recurrent) mutation in C987GG/G987GG, is also noteworthy. Age of mutational events Assuming a generation time of 25 years and based on the most accurate recombination estimates (from family data), the age calculated for the worldwide-spread lineage TTACAC was 5774 ± 1116 years (Table 2). Introduction of this mutation in North America seems to have occurred much later (3396 ± 1182 years ago), in concordance with the genetic diversity found; nevertheless, these ignore the possibility of different introductions of mutations identical by descent in a given population. In Europe, TTACAC seems to be even more recent, with ages less than 2000 years in Germany and France and less than 1000 years in Portugal. As for GTGGCA, the estimated age regardless of its origin was 1416 ± 434 years. Comment This study focuses on the main MJD mutations spread worldwide, explaining the current disease distribution. Analysis of haplotype backgrounds of fast-evolving polymorphisms (STRs) within distinct lineages (identified by genotyping of SNPs, which are rather stable) allowed us to assess their origins and follow their spreading routes. It is expected that, in the absence of strong selective pressures, the ancestral population of a given lineage will show the largest diversity, as mutation and recombination would have more time to generate new STR haplotypes. This should be taken with caution owing to the large variance that stochastic events underlying the evolutionary process carry: a large sample size is needed to infer origin, and inferences based on small samples may not be reliable. Moreover, accumulation of variation is also proportional to the effective population size; in smaller populations, haplotypes will coalesce earlier and thus accumulate less variation. Ttacac lineage The high haplotype diversity among Japanese families as compared with all of the other families with the TTACAC lineage supports its Asian origin. Regarding the high frequency of this SNP-based haplotype among Asian controls as compared with the European and African control populations,20 it could be hypothesized that recurrent mutations have been contributing to this high diversity. However, contrary to other trinucleotide repeat disorders, there are no reports to our knowledge of de novo MJD expansions, and intermediate alleles are rarely found. This is in accordance with a model where very few mutational events account for the cases observed nowadays. The higher levels of diversity among French and German families as compared with Portuguese families do not support the hypothesis of an introduction of the TTACAC lineage in Europe exclusively from Portugal. In Europe, the presence of TTACAC expanded chromosomes is recent and indicates that this has not spread along with the Neolithic culture. Moreover, the low diversity and close phylogenetic relationships among haplotypes of Portuguese and Brazilian families suggest that historical links between these 2 populations are responsible for the presence of MJD in Brazil. In North America, separate introductions of the TTACAC lineage may have occurred directly from the ancestral Asian population, Europe, or both. Gtggca lineage As for GTGGCA, results are more controversial. Although it could be argued that haplotype variability was underestimated owing to sampling effects, this seems unlikely if we consider the high number of families studied and their rather distinct origins; thus, our findings might still be compatible with a Portuguese origin. Moreover, GTGGCA reaches its higher frequency in controls in Europe.20 Taken at face value, the lower diversity among Portuguese haplotypes would imply a non-Portuguese origin for the ancestral line of GTGGCA; nevertheless, most of these chromosomes came from small isolated populations such as the Azores, which could lead to an underestimation of its age. One could contemplate the hypothesis of a non-Portuguese origin with a subsequent expansion within Portugal and loss of the ancestral haplotype; this is, however, a less parsimonious explanation. Age and spread of ancestral mutational events The ancient TTACAC lineage may have occurred during the Bronze Age; the GTGGCA lineage, regardless of its population of origin, may have occurred less than 2000 years BP. Thus, both MJD mutations seem to be fairly recent as compared with other expanded repeat disorders. In Friedreich ataxia, the origin of expanded alleles was dated at about 25 000 years ago, following the population expansion of the Upper Paleolithic23; for type 2 myotonic dystrophy, the founding haplotype was calculated to be 5000 to 12 500 years old24; in Huntington disease, the presence of the ancestral haplotype only in the Land of Valencia (Spain) was estimated to range from 4700 to 10 000 years old.25 The age estimated for the expansion of TTACAC could explain its diffusion westward from East Asia through the trade routes established in central Asia known as the Silk Road; however, these results do not support admixture between Indian and Portuguese MJD populations and are against the hypothesis of MJD in India having been introduced from Portugal19 as its reverse.15 Historical links with Portugal, on the other hand, would be responsible for MJD among Brazilian families. If the GTGGCA lineage indeed has a Portuguese origin, we could speculate about the Azores vs mainland Portugal as its starting point. Taking into account only the haplotype diversity in the Azores together with a recurrent mutation in C987GG/G987GG (a rare event) in that population, the H′1 haplotype could be a recombinant from the ancestral H′2, whereas the high number of mainland Portuguese families with the H′1 haplotype would be explained by a founder effect. Nevertheless, phylogenetic relationships among haplotypes from different populations suggest that more diversity has been generated from H′1, supporting its mutational origin in mainland Portugal. This fits well with the fact that the Azores were inhabited when they were first discovered in 1439, subsequently being colonized mostly from mainland Portugal. Conclusions The scarcity of de novo expansions in MJD led us to question worldwide disease frequencies. Considering the clinical phenomenon of anticipation together with the rarity of intermediate alleles that could act as a reservoir for recurrent full mutations, the extinction of expanded MJD alleles would be expected; however, we must take into account that a theoretical progressive elimination of alleles from the human gene pool must be faced over an evolutionary timescale where other demographic events interfere as well. Moreover, distinct mutational mechanisms may apply for different repeat disorders; it would be interesting to explore cis- or trans-acting factors involved in repeat stability in MJD. Back to top Article Information Correspondence: Jorge Sequeiros, MD, PhD, Unidade de Investigação Genética e Epidemiológica em Doenças Neurológicas, Instituto de Biologia Molecular e Celular, Rua Campo Alegre 823, 4150-180 Porto, Portugal (jsequeir@ibmc.up.pt). Accepted for Publication: November 28, 2006. Author Contributions:Study concept and design: Martins, Coutinho, Amorim, and Sequeiros. Acquisition of data: Martins, Gaspar, Wong, Silveira, Nicholson, Brunt, Tranebjaerg, Stevanin, Hsieh, Soong, Loureiro, Dürr, Tsuji, Watanabe, Jardim, Giunti, Riess, Ranum, Brice, Rouleau, and Coutinho. Analysis and interpretation of data: Martins, Calafell, Silveira, Amorim, and Sequeiros. Drafting of the manuscript: Martins, Amorim, and Sequeiros. Critical revision of the manuscript for important intellectual content: Calafell, Gaspar, Wong, Silveira, Nicholson, Brunt, Tranebjaerg, Stevanin, Hsieh, Soong, Loureiro, Dürr, Tsuji, Watanabe, Jardim, Giunti, Riess, Ranum, Brice, Rouleau, and Coutinho. Statistical analysis: Martins, Calafell, and Amorim. Obtained funding: Amorim and Sequeiros. Study supervision: Calafell, Coutinho, Amorim, and Sequeiros. Financial Disclosure: None reported. Funding/Support: This work was supported in part by Fundação para a Ciência e Tecnologia, through a research grant from Programa Operacional Ciência e Inovação (Instituto de Patologia e Imunologia Molecular da Universidade do Porto), and through Financiamento Plurianual de Unidades (Unidade de Investigação Genética e Epidemiológica em Doenças Neurológicas, Instituto de Biologia Molecular e Celular). Dr Martins was a recipient of scholarship SFRH/BD/8880/2002 from Fundação para a Ciência e Tecnologia. Additional Contributions: Albertino Damasceno, MD, and Benilde Soares, MD, of the Eduardo Mondlane University, Maputo, provided the Mozambican control samples. Nicholas Wood, MD, PhD, provided collaboration. Paula Magalhães, MSc, Susana Carrilho, MSc, and Joana Cerqueira, MSc, provided efficient technical assistance at Centro de Genética Preditiva e Preventiva, Instituto de Biologia Molecular e Celular. We are grateful to the patients and relatives who agreed to collaborate in this project. References 1. Nakano KKDawson DMSpence A Machado disease: a hereditary ataxia in Portuguese emigrants to Massachusetts. Neurology 1972;22 (1) 49- 55PubMedGoogle ScholarCrossref 2. Woods BTSchaumburg HH Nigro-spino-dentatal degeneration with nuclear ophthalmoplegia: a unique and partially treatable clinico-pathological entity. J Neurol Sci 1972;17 (2) 149- 166PubMedGoogle ScholarCrossref 3. Rosenberg RNNyhan WLBay CShore P Autosomal dominant striato-nigral degeneration: a clinical, pathologic and biochemical study of a new genetic disorder. Neurology 1976;26 (8) 703- 714PubMedGoogle ScholarCrossref 4. Silveira IMiranda CGuimarães L et al. Trinucleotide repeats in 202 families with ataxia: a small expanded (CAG)n allele at the SCA17 locus. Arch Neurol 2002;59 (4) 623- 629PubMedGoogle ScholarCrossref 5. Jiang HTang BXu B et al. Frequency analysis of autosomal dominant spinocerebellar ataxias in Han population in the Chinese mainland and clinical and molecular characterization of spinocerebellar ataxia type 6 [in Chinese]. Zhonghua Yi Xue Yi Chuan Xue Za Zhi 2005;22 (1) 1- 4PubMedGoogle Scholar 6. Lopes-Cendes ITeive HGCalcagnotto ME et al. Frequency of the different mutations causing spinocerebellar ataxia (SCA1, SCA2, MJD/SCA3 and DRPLA) in a large group of Brazilian patients. Arq Neuropsiquiatr 1997;55(3B)519- 529PubMedGoogle ScholarCrossref 7. Jardim LBSilveira IPereira ML et al. A survey of spinocerebellar ataxia in South Brazil: 66 new cases with Machado-Joseph disease, SCA7, SCA8, or unidentified disease-causing mutations. J Neurol 2001;248 (10) 870- 876PubMedGoogle ScholarCrossref 8. van de Warrenburg BPSinke RJVerschuuren-Bemelmans CC et al. Spinocerebellar ataxias in the Netherlands: prevalence and age at onset variance analysis. Neurology 2002;58 (5) 702- 708PubMedGoogle ScholarCrossref 9. Takano HCancel GIkeuchi T et al. Close associations between prevalences of dominantly inherited spinocerebellar ataxias with CAG-repeat expansions and frequencies of large normal CAG alleles in Japanese and Caucasian populations. Am J Hum Genet 1998;63 (4) 1060- 1066PubMedGoogle ScholarCrossref 10. Schöls LAmoiridis GBüttner TPrzuntek HEpplen JTRiess O Autosomal dominant cerebellar ataxia: phenotypic differences in genetically defined subtypes? Ann Neurol 1997;42 (6) 924- 932PubMedGoogle ScholarCrossref 11. Stevanin GDürr ADavid G et al. Clinical and molecular features of spinocerebellar ataxia type 6. Neurology 1997;49 (5) 1243- 1246PubMedGoogle ScholarCrossref 12. Moseley MLBenzow KASchut LJ et al. Incidence of dominant spinocerebellar and Friedreich triplet repeats among 361 ataxia families. Neurology 1998;51 (6) 1666- 1671PubMedGoogle ScholarCrossref 13. Storey Edu Sart DShaw JH et al. Frequency of spinocerebellar ataxia types 1, 2, 3, 6, and 7 in Australian patients with spinocerebellar ataxia. Am J Med Genet 2000;95 (4) 351- 357PubMedGoogle ScholarCrossref 14. Giunti PSabbadini GSweeney MG et al. The role of the SCA2 trinucleotide repeat expansion in 89 autosomal dominant cerebellar ataxia families: frequency, clinical and genetic correlates. Brain 1998;121459- 467PubMedGoogle ScholarCrossref 15. Mittal USrivastava AKJain SJain SMukerji M Founder haplotype for Machado-Joseph disease in the Indian population: novel insights from history and polymorphism studies. Arch Neurol 2005;62 (4) 637- 640PubMedGoogle ScholarCrossref 16. Brusco AGellera CCagnoli C et al. Molecular genetics of hereditary spinocerebellar ataxia: mutation analysis of spinocerebellar ataxia genes and CAG/CTG repeat expansion detection in 225 Italian families. Arch Neurol 2004;61 (5) 727- 733PubMedGoogle ScholarCrossref 17. Coutinho PAndrade C Autosomal dominant system degeneration in Portuguese families of the Azores Islands: a new genetic disorder involving cerebellar, pyramidal, extrapyramidal and spinal cord motor functions. Neurology 1978;28 (7) 703- 709PubMedGoogle ScholarCrossref 18. Maciel PCosta MCFerro A et al. Improvement in the molecular diagnosis of Machado-Joseph disease. Arch Neurol 2001;58 (11) 1821- 1827PubMedGoogle ScholarCrossref 19. Gaspar CLopes-Cendes IHayes S et al. Ancestral origins of the Machado-Joseph disease mutation: a worldwide haplotype study. Am J Hum Genet 2001;68 (2) 523- 528PubMedGoogle ScholarCrossref 20. Martins SCalafell FWong VCNSequeiros JAmorim A A multi-step mutation mechanism drives the evolution of the CAG repeat at MJD/SCA3 locus. Eur J Hum Genet 2006;14 (8) 932- 940PubMedGoogle ScholarCrossref 21. Morral NBertranpetit JEstivill X et al. The origin of the major cystic fibrosis mutation (delta F508) in European populations. Nat Genet 1994;7 (2) 169- 175PubMedGoogle ScholarCrossref 22. Cobo AMSaenz APoza JJ et al. A common haplotype associated with the Basque 2362AG→TCATCT mutation in the muscular calpain-3 gene. Hum Biol 2004;76 (5) 731- 741PubMedGoogle ScholarCrossref 23. Labuda MLabuda DMiranda C et al. Unique origin and specific ethnic distribution of the Friedreich ataxia GAA expansion. Neurology 2000;54 (12) 2322- 2324PubMedGoogle ScholarCrossref 24. Bachinski LLUdd BMeola G et al. Confirmation of the type 2 myotonic dystrophy (CCTG)n expansion mutation in patients with proximal myotonic myopathy/proximal myotonic dystrophy of different European origins: a single shared haplotype indicates an ancestral founder effect. Am J Hum Genet 2003;73 (4) 835- 848PubMedGoogle ScholarCrossref 25. García-Planells JBurguera JASolis P et al. Ancient origin of the CAG expansion causing Huntington disease in a Spanish population. Hum Mutat 2005;25 (5) 453- 459PubMedGoogle ScholarCrossref

Journal

Archives of NeurologyAmerican Medical Association

Published: Oct 1, 2007

References

You’re reading a free preview. Subscribe to read the entire article.


DeepDyve is your
personal research library

It’s your single place to instantly
discover and read the research
that matters to you.

Enjoy affordable access to
over 18 million articles from more than
15,000 peer-reviewed journals.

All for just $49/month

Explore the DeepDyve Library

Search

Query the DeepDyve database, plus search all of PubMed and Google Scholar seamlessly

Organize

Save any article or search result from DeepDyve, PubMed, and Google Scholar... all in one place.

Access

Get unlimited, online access to over 18 million full-text articles from more than 15,000 scientific journals.

Your journals are on DeepDyve

Read from thousands of the leading scholarly journals from SpringerNature, Wiley-Blackwell, Oxford University Press and more.

All the latest content is available, no embargo periods.

See the journals in your area

DeepDyve

Freelancer

DeepDyve

Pro

Price

FREE

$49/month
$499/year

Save searches from
Google Scholar,
PubMed

Create folders to
organize your research

Export folders, citations

Read DeepDyve articles

Abstract access only

Unlimited access to over
18 million full-text articles

Print

20 pages / month