|Home | About | Journals | Submit | Contact Us | Français|
A massive hexanucleotide repeat expansion mutation (HREM) in C9ORF72 has recently been linked to amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD). Here we describe the frequency, origin and stability of this mutation in ALS+/−FTD from five European cohorts (total n=1347). Single-nucleotide polymorphisms defining the risk haplotype in linked kindreds were genotyped in cases (n=434) and controls (n=856). Haplotypes were analysed using PLINK and aged using DMLE+. In a London clinic cohort, the HREM was the most common mutation in familial ALS+/−FTD: C9ORF72 29/112 (26%), SOD1 27/112 (24%), TARDBP 1/112 (1%) and FUS 4/112 (4%) and detected in 13/216 (6%) of unselected sporadic ALS cases but was rare in controls (3/856, 0.3%). HREM prevalence was high for familial ALS+/−FTD throughout Europe: Belgium 19/22 (86%), Sweden 30/41 (73%), the Netherlands 10/27 (37%) and Italy 4/20 (20%). The HREM did not affect the age at onset or survival of ALS patients. Haplotype analysis identified a common founder in all 137 HREM carriers that arose around 6300 years ago. The haplotype from which the HREM arose is intrinsically unstable with an increased number of repeats (average 8, compared with 2 for controls, P<10−8). We conclude that the HREM has a single founder and is the most common mutation in familial and sporadic ALS in Europe.
Despite apparent differences in the clinical phenotype of amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD), evidence of an etiopathological link between these disorders is irrefutable. ALS due to motor neuron degeneration usually presents with focal weakness in a limb or mouth/throat muscles (bulbar) and spreads relentlessly causing widespread paralysis.1 FTD presents with changes in behaviour, personality and language due to degeneration of neurons in the frontal and temporal lobes.2 Both disorders can be familial and in a subset of these kindreds, individuals can present with either ALS or FTD, or features of both. In 2006, we reported linkage to a 11-Mb locus on chromosome 9p13.2–21.3 in Dutch and Scandanavian kindreds with autosomal-dominant ALS-FTD.3, 4 Linkage was subsequently confirmed in eight other dominant kindreds defining a minimal overlapping region of ~3.6Mb.5, 6 Genome-wide association studies in sporadic and familial ALS demonstrated highly significant association with single-nucleotide polymorphisms (SNPs) across a 170-Kb region at 9p21.2.7, 8, 9, 10, 11 A massive GGGGCC hexanucleotide repeat expansion mutation (HREM) has recently been identified within intron 1 of C9ORF72 as the pathogenic mutation responsible for familial and sporadic ALS and FTD in these cases.12, 13
Here we describe HREM mutation frequencies in ALS in five European populations. We have generated a detailed map of genetic variation across the locus that provides evidence of genomic instability, which on one occasion gave rise to a massive insertion, generating a single common founder for all the European HREM cases.
DNA was extracted from blood and post-mortem brain frontal cortex, using the standard procedures in patients diagnosed with ALS by the revised El Escorial criteria.14 All the ALS samples were of Northern European Caucasian origin and collected in specialist regional centres following informed consent. ALS cases were designated as familial if one or more first- or second-degree relatives developed ALS or FTD. A person was classified as having ALS+FTD if they presented with major cognitive or behavioural change at any stage during the course of their illness. Patients provided consent conforming to local and national ethics committee guidelines. In the London Clinic samples, mutations in SOD1, VCP, OPTN and UBQLN2 (all exons), TDP43 (exon 6) and FUS (exons 14 and 15) were screened and any positive samples were excluded from further analysis. Familial samples from the other European cohorts were screened and excluded for SOD1, FUS and TARDBP mutations.
DNA from 12 individuals carrying the disease haplotype and 4 without from our previously linked kindred,3 2 affected members from the previously published linked Scandanavian family,4 14 cases with suggested linkage to ch9p and 21 other individuals with familial ALS+/−FTD were processed for DNA capture using custom-designed overlapping probes (Roche Nimblegen, Madison, WI, USA) across the 3.6-Mb locus between D9S169 (27238617)5 and D9S251 (30819382).6 A total of 5μg of DNA was fragmented with a Bioruptor (Wolf Laboratories Ltd, York, UK) at 30s on/off bursts for 45mins to sizes of 200–300bp. End repair was followed by addition of adenine ends and ligation of adaptors (Illumina, Little Chesterford, Essex, UK) and peak sizes checked using a Bioanlayzer (Agilent Technologies, Wokingham, UK). Purification steps were conducted with SPRI beads according to the manufacturer's instructions (Beckman Coulter Genomics, High Wycombe, UK). Libraries were hybridized with 4.5μl of locus probes for 72h at 47°C, washed and bound to streptavidin beads, followed by PCR of 10 separate reactions per library. Individual reactions were pooled, cleaned using QIAquick (Qiagen, Crawley, UK), quantified by chip (DNA 1000, Agilent Technologies) and sequenced with 76 or 100bp paired-end reads on GAII and Hiseq Analyzers (Illumina).
Raw sequencing data were mapped to the human reference genome (Build hg19) using Novoalign (http://www.novocraft.com/) and processed using Picard tools v1.35 and the Genome Analysis Toolkit (GATK, version V1.1) to produce a ‘clean' BAM file.15 SNP and Indel calling was performed using the Unified Genotyper module in GATK in batch mode. The resulting Variant Call Format (VCF, version 4.0) file was annotated using Variant Filtration in GATK set as follows: QUAL<30.0||QD<5.0||HRun >10||SB>−5.00||DP<10 and cluster size 10. The VCF file was converted to ‘pedigree' format using vcftools v1.3.1 allowing us to phase all SNPs on the risk allele (http://vcftools.sourceforge.net).16 A PERL script was written to identify sequencing reads from the fastq files overlapping the HREM, and to count the numbers of repeats found within each.
A total of 82 SNPs spanning the locus that were shared among the affected individuals from the ch9p-linked families and highly represented in familial ALS/FTD cases were Kaspar genotyped (Kbiosciences Ltd, Hoddesdon, UK) in a cohort of 434 cases and 856 controls of European ancestry from Sweden, Belgium, England and Italy. In all, 16 cases from the locus-capture set were included to validate next-generation genotypes. Haplotypes were generated and frequencies were determined using PLINK v1.07,17 and phasing was corroborated using Snphap (D. Clayton; http://www-gene.cimr.cam.ac.uk/clayton/software). Hardy–Weinberg equilibrium was assessed by a χ2 test for quality control. DMLE+ v2.318 was used to estimate the age of the HREM via the expected relationship between HREM allele frequency, local linkage disequilibrium and population growth rate. For comparison, we also used a decay in linkage disequilibrium method (Equation 1),19 averaging the age estimates across all 82 SNPS. Between-SNP genetic distances were estimated using LDhat applied to HapMap Phase II (The International HapMap Consortium, 2007). The 82 SNPs span 110186 bases, with the HREM estimated to be 105131 bases from the telomeric SNP. We assumed that a random population sample of 39091 would be expected to yield 137 HREM-bearing individuals, based on our observed frequency of 3 in 856 controls, and we assumed that such a sample would form a fraction of 7.8 × 10−5 of all Europeans, based on a current population size of 500 million. We performed DMLE+ using a burn-in of 20000 iterations followed by runs of 100000 iterations, with population growth rates in Europe of 5%,20 with a lower limit of 2.5%21 and an upper limit of 8.5%,22 and with a 25-year inter-generation time, with lower and upper limits of 20 and 30 years respectively.
Exons 1A, 1B and the intervening intron were aligned from human, chimpanzee, gorilla, orangutan, mouse and rat reference genomes using ClustalW and manually edited in GeneDoc.
A total of 1347 ALS+/−FTD patients, including the 434 cases and 856 controls used in the haplotyping study, were screened for the GGGGCC HREM, using repeat primer PCR13 with a final concentration of 7% DMSO, 1 betaine, 0.17m of 7-deaza-2-deoxy GTP, 0.7–1.4μ of primer mix, 0.85m of MgCl2, 50% Applied Biosystems True Allele PCR Premix (Applied Biosystems, Warrington, UK) and 100ng of genomic DNA. Primers included a FAM-labelled reverse primer, one repeat-specific forward primer with an attached anchor sequence and the same anchor sequence as an independent forward primer. Cycling conditions were denaturation 95°C for 15mins and touchdown from 70 to 56°C with 3min extension. Fragment analysis was conducted on an ABI 3130 genetic Analyser and peaks visualized using Genemapper 4.0 (Applied Biosystems). Chromatograms were scored as mutant (sawtooth pattern) or wild type (<30 repeats). Direct sequencing of 48 cases and controls without the HREM was performed using Big Dye V1.1 chemistry and an ABI 3130 genetic analyser to validate repeat primer PCR genotypes, forward primer 5′-GGTTTAGGAGGTGTGTGTTTTTGT-3′, reverse primer 5′-CCAGCTTCGGTCAGAGAAAT-3′ and identical cycling conditions with two extra cycles at each stage of the touchdown protocol.
Unless otherwise stated, all calculations were performed in IBM SPSS v19 (SPSS Inc., Chicago, IL, USA) with two-sided significance tests. Independence of categorical variables was tested using the χ2 distribution. For small cell counts, the Fisher's exact test was used. Alleles of the highly polymorphic non-expanded hexanucleotide repeat were tested for association using Monte-Carlo simulation in the program CLUMP,23 which generates empirical P-values for observed χ2 tables, accounting for the multiple testing inherent in having multiple alleles at a locus. Age of onset and disease duration was tested for association with the HREM using Kaplan–Meier product limit estimate and the log rank test.
Mutations in the familial ALS+/−FTD cohort from the King's College Hospital, London clinic, were identified in 55% of all familial cases with the following frequency: C9ORF72 29/112 (26%), SOD1 27/112 (24%), FUS 4/112 (4%) and TARDBP 1/112 (1%). HREM mutations were also detected in 13/216 (6%) of unselected sporadic ALS cases from the same clinic. No mutations were identified in VCP, OPTN or UBQLN2.
Combining data from five European populations, (detailed individually in Table 1), the HREM in C9ORF72 was detected in 226/1347 (17%) of all ALS+/−FTD cases, in whom known ALS genes had been excluded, and 3/856 (0.3%) controls (Fisher's exact test P-value for allelic association=4.12 × 10−47; OR=57, 95% CI=17.7–224.6). The highest frequency was in familial ALS+FTD kindreds (48/67, 72%) but it was also prevalent in pure ALS kindreds (89/228, 39%), with the total familial frequency therefore being 46% (137/296, P-value 6.13 × 10−89; OR=244, 95% CI=74.4–974.3) In sporadic ALS+/−FTD, HREM frequencies across Europe were higher than for any other known gene at 87/1048 (8%) (P-value 1.1 × 10−19; OR 25.7, 95% CI=7.8–102). Given that sporadic ALS accounts for 95% of all cases, then sporadic ALS+/−FTD cases with the HREM outnumber familial by a ratio of 4:1. Frequencies of the HREM in familial ALS+/−FTD were high but showed considerable variation by country: 19/22 (86%) in Belgium, 30/41 (73%) in Sweden, 10/27 (37%) in the Netherlands, 73/185 (39%) in England and 4/20 (20%) in Italy.
Phenotypic data were available on 189 ALS cases with the HREM and 870 cases without HREM (Table 2). The male:female ratio in HREM-positive cases was 1.1:1 compared with non-HREM cases 1.8:1 (P=0.009), which is similar to population-based studies of familial and sporadic disease (Table 2). Patients with the expansion were more likely to present with cognitive/behavioural and bulbar symptoms than those without (P=0.02). Kaplan–Meier estimates showed no difference between the two groups in the age at onset (P=0.27) or disease duration (P=0.34) (Supplementary Figure 1).
Sequencing of DNA captured across the 3.6-MB locus in 53 individuals generated 1.2 billion reads with an average of 487-fold depth across the region. We identified 10604 SNPs that passed QC but no variants segregating with disease were identified in C9ORF11, MOB3B, IFNK, C9ORF72 or LINGO2 or other predicted genes within the locus. The largest number of GGGGCC repeats detected within intron 1 of C9ORF72 was 8, occurring at the end of a single read in one affected individual. The pathological HREM was not identified because it is so GC-rich and fails to amplify by PCR without a repeat-specific primer. Thus, it is not surprising that none of the variant-calling algorithms we used detected this polymorphism.
We were able to phase 82 informative SNPs within the linked kindreds that defined a shared haplotype across the locus. These were further genotyped in 433 cases and 856 controls and correlated with the presence of the HREM. Detailed inspection of the SNP haplotype in the 137 HREM-positive cases revealed that a full 82-SNP haplotype existed in the vast majority of cases (111/137, 81%, P=8.33 × 10−17) and 3 controls who were positive for the HREM. Despite significant recombination, alleles from the linked haplotype were always preserved in regions flanking one or other side of the HREM (80 SNPs telomeric and 2 SNPs centromeric), providing clear evidence of a single common founder in these European populations (Figure 1). The two SNPs centromeric to the HREM were rs11789520 and rs73440960, possessing allele frequencies in HREM-positive cases of 0.97% (P=0.00023) and 0.93 (P=0.00006), respectively, which demonstrate that the conserved haplotype straddles the expansion region.
Estimates of the age of the HREM using DMLE+ depend primarily on the growth rate of the population and generational interval, which are known to vary greatly over time (Table 3). Growth rates ranging from 2.5–8.5% and intergenerational intervals have been proposed between 20 and 30 years for most founder studies. If we take 5% as a conservative estimate of growth and 25 years as an intergenerational average, we estimate that the mutation arose around 251 generations ago, which equates roughly to 6300 years (see Table 3 for estimates based on a range of growth rates and intergenerational intervals). For comparison, we applied an alternative linkage disequilibrium method19 and estimated that the mutation arose 131 generations ago (3300 years ago assuming a 25-year intergeneration time). We acknowledge the limitations of these analytical tools but it is encouraging that these figures are not greatly disparate.
The human, chimpanzee and gorilla reference genomes contain three copies of the GGGGCC repeat (Figure 2a). Evidence from the NCBI Trace Archive database shows that chimpanzees can also have five or six repeats. No other species appear to contain the hexanucleotide motif, although orangutans possess a possible precursor sequence, 5′-GAGGCCGGGCCC-3′. Phylogenetic analysis of the human haplotypes reveals two main clades, one of which gave rise to the expansion mutation (Figure 2b).
The full background 82-SNP haplotype from which the HREM arose is present in almost all populations studied in the ‘1000 Genome' database of 1046 individuals (Figure 3). Its frequency in people of European ancestry averages at ~15.1%, which is nearly identical to the frequencies we derived from our analysis of controls of 262/1706 (14.9%) and our ALS+/−FTD cohort of 109/740 (14.8%). Repeat primer PCR genotyping does not give the number of repeats for the HREM allele, however, the longest number of repeats can be counted in cases without the expansion using fragment analysis. Sanger sequencing of 48 individuals showed a perfect correlation between the repeat number counted by fragment analysis and sequencing (Supplementary Figure 2). We measured the longest number of repeats in 1154 individuals and compared those with the background haplotype (r) to all other haplotypes (Figure 2c). The average number of repeats in those carrying haplotype (r) was 8 with a widespread of expanded alleles up to 26 (95% CI=4–13), whereas the most prevalent number of repeats in all the other haplotypes was only 2 (95% CI=1–7, P<10−8). This indicates that the background haplotype on which the expansion arose is intrinsically unstable, tending to generate longer repeats.
We have also identified that rs2492816 independently tagged a repeat number >2 for the risk allele (P<1.0E-13, Fisher's exact test) and a repeat number of 2 for the non-risk allele (P<1.0E-52, Fisher's exact test), which accounts for the apparent bimodality of the non-risk haplotype distribution.
We have demonstrated that the hexanucletotide repeat expansion mutation in C9ORF72 is the most common genetic cause of familial ALS+/−FTD across Europe, accounting for 20–86% of genetically undiagnosed familial cases, particularly where FTD and ALS co-segregate and in those presenting with bulbar or cognitive/behavioural symptoms. It is difficult to make robust conclusions about the origin of differences in frequency due to the small sample sizes for each country but they probably reflect the influence of a founder effect. In unselected patients from a London clinic with a family history of ALS (5–10% of all cases), a genetic diagnosis can now be confirmed in 55% of cases. C9ORF72 HREM was the most common mutation (26%) followed by SOD1 (24%), FUS (4%) and TARDBP (1%). The HREM is detected in 6% of sporadic ALS cases but is also present in the background population (0.3%). The penetrance of the HREM appears to be low, given that there is a common founder and the ratio of sporadic to familial HREM cases is 4:1. This is consistent with the incomplete penetrance reported in many linked kindreds. Further work is required to generate figures for age-related penetrance that can be used in genetic counselling and predictive gene testing.
Following exhaustive sequencing, we have confidently identified a haplotype that proves that all HREM carriers arose from a single common founder. We phased a 82-SNP haplotype within linked kindreds that is conserved in its entirety in the majority of all the HREM carriers and flanks at least in part all of the remaining HREM cases. The most economical explanation is that the expansion mutation arose on just one occasion in the European population, however, we cannot exclude the possibility that it arose on multiple occasions on the same background haplotype (r). Estimates of founder age depend heavily on estimates of population growth rates, with smaller rates leading to older estimates, and to a lesser extent on intergenerational interval. Historical evidence is that growth rates have varied greatly being much slower in the distant past than in the last century.24 Using averaged figures of a 5% growth rate and 25-year interval, we have estimated that the founder mutation arose around 6300 years ago (range 4400–8600 years).
A common founder was originally proposed for Finnish FALS cases based on a 42-SNP haplotype.8 A subsequent meta-analysis of genome-wide association study data from five European populations (Finnish, Irish, UK, US and Italian) reduced this to a common 20-SNP risk haplotype.25 In the original report of the HREM by the same group, however, only two-thirds of Finnish cases were reported to have a common haplotype, implying that the other third of their HREM cases may have different founders.13 By fine mapping across the locus in great detail, we have shown that all the HREM carriers (cases and controls) have conserved SNPs that flank one or other side of the HREM, confirming that all the carriers arose from a single founder haplotype (r). Given that the mean age at onset of our HREM cases is 60 years (see Supplementary Figure 1) and the penetrance is relatively low, we doubt that significant selection pressures would apply over past millennia where life expectancy was considerably lower. For these reasons, we would not expect the mutation to die out due to selective purification.
We have uncovered evidence that the GGGGCC repeat arose during primate evolution and is highly polymorphic but the biological significance of this is unknown. Haplotype frequencies in different ethnic populations from the 1000 Genomes database strongly suggest a European origin for the background 82-SNP haplotype (r) on which the HREM arose. The maximum number of repeats on either allele is much greater in those with the (r) haplotype than all other haplotypes combined. Nearly 50% of the individuals with non-(r) haplotypes have a maximum of two copies (295/641) compared with 5% (29/513) of individuals with the (r) haplotype, where the average is eight repeats and some individuals have many more. This difference in repeat number confirms initial observations based on a single SNP rs3849942 marker for the HREM risk haplotype.12 It is not clear why the (r) haplotype is prone to expansion but it is possible that 8–26 repeats, which are GC-rich, promote the formation of hairpin secondary loop structures that impair DNA replication. For instance, flap endonuclease 1 required for normal maturation of Okazaki fragments during replication fails to process flaps folded into aberrant hairpin structures and is thought to cause expansion at CAG repeats.26, 27 Alternatively, an independent de novo event may have occurred in a single person 6300 years ago, which affected the fidelity of DNA polymerase or a DNA mismatch repair enzyme, which in conjunction with the unstable repeat region resulted in the HREM.28, 29
The dominant pathology in 90% of ALS and tau-negative FTD inclusions contain the TAR DNA-binding protein (TDP-43) within the cytoplasm of neurons and glia.30 TDP-43 inclusions are also prominent in cases linked to chromosome 9p31 but HREM-specific pathology includes abundant cytoplasmic and intranuclear p62-positive inclusions in the hippocampus and cerebellum that are TDP-43-negative.32, 33 Precisely, how the HREM causes TDP-43 mislocalisation and neurodegeneration is not currently known. Evidence that the HREM reduces levels of C9ORF72 transcripts implicates a loss of function, however, probes detecting the HREM transcript identified RNA foci within the nuclei of neurons in the frontal cortex and spinal cord.12 In other dominant intronic repeat disorders, such as Myotonic Dystrophy (DM1), these foci have been shown to sequester RNA-binding proteins, which cause a range of deleterious changes in RNA processing.34, 35
We propose that all C9ORF72 HREM cases derive from a single common founder and are now the most common cause of familial and sporadic ALS in Western Europe. The GGGGCC repeat is highly polymorphic and particularly unstable in the context of a specific haplotype (r), but the massive pathogenic expansion may have arisen on just one occasion around 6300 years ago. Although gene testing will become widely available, further work is required to establish the disease risk for HREM carriers.
This work was supported by the NIHR Biomedical Research Centre for Mental Health at the South London and Maudsley NHS Foundation Trust and Institute of Psychiatry, Kings College London, and the SLAGEN consortium. We would like to thank Antonia Ratti, Cinzia Tiloca, Barbara Castellotti, Viviana Pensato, Stefania Corti, Roberto del Bo, Gianni Sorarù, Carla D'Ascenzo, Sandra D' Alfonso, Lucia Corrado, Cristina Cereda, Ceroni Mauro and Isabella Fogh for their help. This work was funded by the Medical Research Council, Motor Neuron disease Association (UK), American ALS Association and the Heaton-Ellis Trust. The research leading to these results has received funding from the European Community's Seventh Framework Programme (FP7/2007-2013) under the grant agreement number 259867.
The authors declare no conflict of interest.
Supplementary Information accompanies the paper on European Journal of Human Genetics website (http://www.nature.com/ejhg)