Skip to main content

Tissue-specific and repeat length-dependent somatic instability of the X-linked dystonia parkinsonism-associated CCCTCT repeat

A Correction to this article was published on 25 April 2022

This article has been updated

Abstract

X-linked dystonia-parkinsonism (XDP) is a progressive adult-onset neurodegenerative disorder caused by insertion of a SINE-VNTR-Alu (SVA) retrotransposon in the TAF1 gene. The SVA retrotransposon contains a CCCTCT hexameric repeat tract of variable length, whose length is inversely correlated with age at onset. This places XDP in a broader class of repeat expansion diseases, characterized by the instability of their causative repeat mutations. Here, we observe similar inverse correlations between CCCTCT repeat length with age at onset and age at death and no obvious correlation with disease duration. To gain insight into repeat instability in XDP we performed comprehensive quantitative analyses of somatic instability of the XDP CCCTCT repeat in blood and in seventeen brain regions from affected males. Our findings reveal repeat length-dependent and expansion-based instability of the XDP CCCTCT repeat, with greater levels of expansion in brain than in blood. The brain exhibits regional-specific patterns of instability that are broadly similar across individuals, with cerebellum exhibiting low instability and cortical regions exhibiting relatively high instability. The spectrum of somatic instability in the brain includes a high proportion of moderate repeat length changes of up to 5 repeats, as well as expansions of ~ 20- > 100 repeats and contractions of ~ 20–40 repeats at lower frequencies. Comparison with HTT CAG repeat instability in postmortem Huntington’s disease brains reveals similar brain region-specific profiles, indicating common trans-acting factors that contribute to the instability of both repeats. Analyses in XDP brains of expansion of a different SVA-associated CCCTCT located in the LIPG gene, and not known to be disease-associated, reveals repeat length-dependent expansion at overall lower levels relative to the XDP CCCTCT repeat, suggesting that expansion propensity may be modified by local chromatin structure. Together, the data support a role for repeat length-dependent somatic expansion in the process(es) driving the onset of XDP and prompt further investigation into repeat dynamics and the relationship to disease.

Introduction

X-linked dystonia parkinsonism (XDP, OMIM314250) is a progressive and fatal adult-onset neurodegenerative disease endemic to the island of Panay, Philippines [1, 2]. The clinical phenotype most commonly described consists of an initial presentation of focal dystonia that spreads to multiple body regions and combines with, or is replaced by, parkinsonism that predominates 10–15 years after onset [1, 3, 4]. The average age of symptom onset is 39–40 years, though the age at onset (AAO) can differ widely (12 to 64 years) [1, 2]. XDP principally affects males, with a frequency of 5.74 cases per 100,000 individuals in Panay, though female carriers are reported to have symptoms in a few cases [1, 5].

Limited neuropathological studies of post mortem XDP patient brain tissue have revealed changes to the neostriatum that include selective loss of medium-spiny neurons (MSNs) [6,7,8] as seen in Huntington’s disease (HD, OMIM 143,100) [9], with a selectivity towards the striosomal compartment [7] as well as evidence for pathological changes outside the neostriatum, including cortex, substantia nigra and cerebellum [10, 11]. Neuroimaging has demonstrated volume loss of the caudate, putamen and pallidum [12,13,14,15] as well as broad neuroanatomical changes across multiple brain regions, e.g. patterns of grey matter gain and loss in cerebellum and brainstem and extensive white matter abnormalities that include thinning of frontal and temporal cortices [10, 12, 15]. Interestingly, white matter changes were not apparent in the occipital cortex, potentially reflecting the lack of anatomical connection with the striatum [12], while increased connectivity was observed between the putamen and the paralimbic cortex (enriched for connections with striosomes), in particular the insula [14]. Overall, XDP pathology reflects prominent neuronal loss in the caudate and putamen as well as the involvement of multiple other brain regions and the connections between them.

Genetic linkage and refined mapping localized the causal locus of XDP to the X-chromosome [16,17,18], with recent work characterizing a thirteen-marker haplotype shared by all probands defining a minimal critical region of 219.7 kb with TATA-binding-protein (TBP)-associated factor-1 (TAF1) being the only gene within this region [19]. Among the thirteen disease-specific variants is a ~ 2.6 kb SINE-VNTR-Alu (SVA)-type retrotransposon [20] inserted in intron 32 of TAF1 [18]. XDP patient tissues and cell lines exhibit reduced TAF1 expression [18, 19, 21,22,23,24] as well as aberrant splicing that results in partial retention of intronic sequence proximal to the SVA insertion [19]. Reduced TAF1 expression, intron retention and aberrant splicing can be rescued by excision of the SVA [19, 22], suggesting that SVA-mediated TAF1 transcriptional dysregulation may contribute to disease pathogenesis.

The 5’ end of the SVA contains a hexameric CCCTCT repeat tract that varies in length from 30 to 55 repeats [4, 25]. Notably, repeat length is inversely correlated with AAO, as seen in other disorders caused by expanded microsatellite repeats [26], suggesting a critical role of CCCTCT repeat length in XDP pathogenesis. The length of the repeat was also associated with transcriptional activity in vitro [4] and its length inversely correlated with TAF1 expression in patient blood samples [25]. A common characteristic of repeat expansion disorders is the instability of the disease-associated repeat, both in the germline and in somatic tissues, where in the latter the repeat tends to expand in a length-dependent and tissue-specific manner [26,27,28]. In HD, genetic studies have provided strong evidence that somatic expansion of the HTT CAG repeat drives the rate of disease onset [29,30,31]. Studies of other repeat expansion diseases indicate that somatic expansion is a likely common mechanism driving pathogenesis [27, 32,33,34,35]. Significantly, a recent genome-wide association study (GWAS) for modifiers of XDP [36] identified genes (MSH3, PMS2) with known roles in in repeat instability [30, 37,38,39] that also modify HD [30, 40], indicating that a common mechanism at the level of repeat instability extends to XDP. The XDP CCCTCT repeat exhibits intergenerational instability, with repeat length tending to increase in transmissions from mothers and to decrease in transmissions from fathers [4, 25]. Patient cell lines show limited repeat instability [4, 25], while investigation of a small number of XDP individuals has provided evidence of somatic repeat expansion in post-mortem brain [25, 41].

Here, to gain a deeper understanding of somatic instability in XDP we have performed an extensive quantitative characterization of XDP CCCTCT repeat instability in blood, and in up to 17 brain regions from 41 XDP individuals. Our findings reveal repeat length- and tissue-dependent CCCTCT repeat expansion, suggesting that somatic expansion underlies the repeat length-dependent clinical onset of XDP.

Materials and methods

XDP patients and sample collection

Blood

Patients recruited for this study included individuals with XDP evaluated at Massachusetts General Hospital (MGH) (Boston, MA, USA), Jose R. Reyes Memorial Medical Center (JRRMMC) (Manila, Philippines), and regional clinics on the island of Panay (Panay, Philippines). All participants provided written informed consent, and the study was approved by local Institutional Review Boards (IRBs) at both MGH and JRRMMC. Patients enrolled were subjected to comprehensive neurological examinations and blood collection [42]. This study also included archival DNA specimens; collection methods and the clinical characterization of donor subjects who provided these specimens have been previously described [17]. Genomic DNA (gDNA) was extracted from blood using the Gentra Puregene kit (Qiagen). Enrolled patients were confirmed to be positive for the XDP mutation by PCR amplification for a known 48 bp deletion haplotype marker as previously described [4, 43]. Blood samples from 266 male XDP patients with known AAO were evaluated for correlation with repeat length. Somatic instability was analyzed in 164 blood samples, representing a subset of male XDP patients included in the cohort above.

Brain

Post-mortem brain tissue from XDP patients (n = 41; 40 with age at onset and death) was obtained in collaboration with the Collaborative Center for XDP (CCXDP), at MGH (Boston, MA, USA), Makati Medical Center (Makati City, Philippines), and the Sunshine Care Foundation (Panay, Philippines). Detailed descriptions of all methods related to donor consent, brain collection and tissue processing have been previously reported [44] and the use of XDP patient post-mortem brain tissue and all study procedures were approved by Institutional Review Boards at Makati Medical Center (Makati City, Philippines) and MGH (Boston, MA, USA). Genomic DNA was extracted from different brain regions using the DNeasy Blood and Tissue Kit (Qiagen), according to manufacturer’s instructions and with the following modifications: samples were incubated in buffer ATL and Proteinase K overnight at 56 °C; washes AW1 and AW2 were repeated; DNA was eluted in 100 μl of Qiagen Elution Buffer, preheated to 56°C, applied to the spin columns, and incubated at room temperature for 10 min before centrifugation. The sample was run through the spin column a second time before final centrifugation. The presence of the XDP mutation in each brain was confirmed as above.

Determination of XDP and LIPG CCCTCT repeat lengths and expansion indices

To determine the length of XDP and LIPG SVA CCCTCT repeats in blood and postmortem brain regions, we used fluorescent PCR-based assays, with the primers and conditions outlined in Additional File 1: Table S1. Both protocols used 125 ng of gDNA per reaction, in a 25 µl reaction volume with buffer and dNTPs provided with the PrimeSTAR GXL polymerase (Takara) according to the manufacturer’s protocol, and as previously described for the XDP repeat [4]. Following PCR, aliquots of each product were resolved via electrophoresis in agarose gels to confirm amplification of the SVA repeat sequence and then run on the ABI 3730 DNA sequencer (Applied Biosystems) with GeneScan 500 LIZ as internal size standard, and the data analyzed using GeneMapper v5 (Applied Biosystems) [4]. Repeat amplification resulted in a distribution of fragments separated by 6 bp and repeat size was defined as the tallest peak in this distribution. XDP repeat size was assigned relative to DNA standards of known repeat lengths (37, 44, 50) determined relative to a sequenced bacterial artificial chromosome with 35 repeats [19] and LIPG repeat size calculated based on fragment length (bp) and the number of repeats (39) in the reference (GRCh38) genome. To quantify XDP and LIPG CCCTCT repeat expansion, we generated an expansion index from the GeneMapper peak height data as previously described [45], using a 5% relative peak height threshold cut-off (i.e. excluding peaks whose height is less than 5% of the height of the modal allele). Because LIPG is autosomal, most individuals had two distinguishable allele lengths. In many cases, alleles were sufficiently separated to allow quantification of expansion peaks from each. In some individuals, when the alleles were too close, we only captured the expansion index from one allele.

Small pool-PCR Southern blot analyses

1 µg of gDNA was digested with HaeIII (37 °C for 12 h) and the enzyme subsequently inactivated at 80 °C for 20 min. Serial dilutions were made in water to a final concentration of 90 pg/µl and 1 µl (approx. 30 genome equivalents, g.e.) was used for PCR amplification using a non-FAM-labeled version of the XDP SVA hexamer primers with the small pool-PCR conditions outlined in Additional File 1: Table S1. For each sample, PCR amplifications of 36 replicates of 90 pg gDNA and 8 DNA-negative PCR controls were carried out in a 25 µl reaction volume with buffer and dNTPs provided with the PrimeSTAR GXL polymerase (Takara) according to the manufacturer’s protocol. 10 µl of each PCR product were run in 2% agarose gels alongside digoxigenin (DIG)-labeled size markers VII and VIII (Roche), for 16 h at 50 V then transferred to a positively charged nylon membrane (Roche) by common squash-blotting technique [25, 46]. The membrane was hybridized with 5 pmol/ml of a 5’ DIG-labeled (AGAGGG)10 probe (Sigma) in DIG Easy Hybridization Solution (Sigma) overnight at 45 °C and then washed twice each with 2 X SSC, 0.1% SDS at room temperature for 5 min, 0.1 X SSC, 0.1% SDS at 68 °C for 20 min, and 0.1 X SSC, 0.5% SDS at 68 °C for 20 min. DIG detection was carried out using the DIG Luminescent Detection system (Sigma) with CPSD substrate according to the manufacturer’s instructions.

Single molecule small pool-PCR sizing

1 µg gDNA was digested with HaeIII as above and the DNA serially diluted to a range of concentrations spanning 3 pg/ul to 300 pg/ul corresponding to approximately 0.5–100 diploid g.e/µl, respectively. For each sample, at least 10 PCR reactions of 1 µl DNA inputs were run for each dilution and resolved using the ABI 3730 DNA Sequencer. Poisson analysis was used to determine empirically for each sample the concentration of DNA that resulted in single molecule PCR amplification, i.e. the concentration that resulted in ~ 33% of all DNA input reactions having no product. We then ran, for each sample, at least three 96-well plates, each consisting of 72 replicates of the optimized single molecule amplifiable DNA amount, 18 DNA-negative PCR controls, the XDP repeat sizing standards (37, 44, 50) as above, and one empty well for machine control. PCR conditions for small pool-PCR were as described in Additional File 1: Table S1, and CCCTCT repeat size was determined as described above. Allele lengths between 330 and 560 bp (about 32–70 repeats) could be accurately determined based on the known repeat sizing controls. For PCR products falling outside of this range (49 out of a total of 1519 sized alleles) we estimated repeat length based on molecular weight. All peaks with heights >  = 150 were sized, and for each plate we verified that all of the no-DNA input wells were negative and that at least 1/3 of the DNA input wells were negative.

HD sample data

In this study we used HTT CAG repeat expansion data previously generated and reported from 8 postmortem brain regions from three HD individuals (HD1-3; CAG repeats 43/16, 44/17, 53/19) obtained from the New York Brain Bank under an approved MGH IRB protocol [28]. The data from a subset of eight tissues used in this study were chosen because they were identical or as close as possible to the XDP brain regions from our XDP cohort. Regions compared were: BA9 (HD and XDP), BA17 (HD) and occipital cortex (XDP), caudate accumbens and putamen (HD) and caudate (XDP), cerebellum (HD and XDP), cingulate gyrus (HD and XDP), globus pallidus putamen (HD) and putamen (XDP), hippocampal formation (HD) and hippocampus (XDP), subthalamic nucleus (HD and XDP), temporal pole (HD and XDP). For simplicity, in Fig. 5c we refer to all the regions according to the XDP labels. Somatic HTT CAG expansion indices were determined for this study using a 5% relative peak height threshold cut-off for comparison to the 5% threshold XDP CCCTCT expansion indices.

Statistical analysis

Data analysis and plots were generated using R/RStudio V.1.3. (https://cran.r-project.org/mirrors.html). Linear regression, stacked bars and scatter plots were generated using ggplot2 package (https://www.rdocumentation.org/packages/ggplot2/versions/3.3.5). Pearson or Spearman coefficients were determined using ggscatter package and used as appropriate when data distribution was Normal or not, respectively. The heatmap was generated from a scaled dataset using the heatmaply package followed by a clusterization method, based on Manhattan distance https://cran.r-project.org/web/packages/heatmaply/vignettes/heatmaply.html. Multiple pairwise comparison test was performed using Wilcoxon rank-sum test followed by Bonferroni Post Hoc method for P-value adjustment. X2 test was used to compare the numbers of events from single molecule SP-PCR data across brain tissues. P-value < 0.05 was considered significant.

Results

XDP CCCTCT repeat length inversely correlates with ages at onset and death but not with disease duration

We previously demonstrated in a cohort of 140 XDP males that CCCTCT repeat length in blood was inversely correlated with AAO [4]. This observation was subsequently confirmed in an independent cohort of 295 individuals [25]. Here, we have used an expanded dataset from our original sample of 140 comprising blood (n = 266) and brain (n = 40) DNA samples from clinically confirmed male XDP patients to examine further the relationship between CCCTCT repeat length and AAO, as well as age at death (AAD) and disease duration, defined as AAD minus AAO (Fig. 1). In these analyses, brain repeat length was determined in 40 postmortem samples with AAO (n = 39 in cerebellum and n = 1 in occipital cortex where cerebellum was not available). Both blood and brain tissue were available for 21 individuals; of these, blood and brain (cerebellar) repeat lengths were identical in 17 individuals and differed by one repeat in 4 individuals (17–012, 19–017, 19–021, 21–031; Additional file 1: Table S2). Mean (± SD) repeat lengths in blood and brain were 41.6 ± 3.9 (range:34–53) and 41.8 ± 4.6 (range:34–55), respectively. Mean (± SD) AAO of the blood and brain samples were 41.4 ± 8.3 (range:18–65) and 41.4 ± 8.7 (range:26–59) years, respectively. Blood repeat length inversely correlated with AAO and explained ~ 45% of the AAO variability (P = 7.7-36; Fig. 1a, red dots), consisted with previous studies [4, 25, 36]. A similar correlation was observed between brain repeat length and AAO, with repeat length explaining ~ 55% of the AAO variability (P = 4.7-08; Fig. 1a, blue dots). There was no difference in AAO-repeat length correlation between individuals exhibiting primarily dystonia at onset (N = 194) and those exhibiting primarily parkinsonism at onset (N = 43) (Additional file 2: Fig.S1), consistent with previous observations [25]. Both AAO and AAD were available for 68 individuals, 28 of whom had blood repeat sizing and 40 of whom had brain repeat sizing. As repeat length was largely identical between brain and blood for the individuals with both measures, we used a combined blood and brain dataset from these 68 individuals to examine relationships between repeat length and AAO, AAD or duration (Fig. 1. b-d). Mean (± SD) repeat length in these 68 individuals was 41.6 ± 4.4 (range:34–55), mean (± SD) AAO was 41.7 ± 4.4 (range:26–64), and mean (± SD) AAD was 50.7 ± 9.5 (range:30–69) years. Repeat length inversely correlated with AAO and AAD, explaining ~ 53% (P = 2.3-12) and ~ 42% (P = 2.5-09) of the AAO and AAD variability, respectively (Fig. 1b and c). In contrast, we found no significant correlation between repeat length and disease duration (AAD-AAO) (Fig. 1d). These data indicate that the length of the CCCTCT repeat is critical for process(es) driving XDP onset and death that ensues, though has no obvious effect or a weaker effect on duration.

Fig. 1
figure 1

Length of the CCCTCT repeat correlates with AAO and AAD in male XDP patients. a Inverse correlations between CCCTCT repeat length in blood (red dots and line, n = 266) and brain (blue dots and line, n = 40) with AAO. As this is a founder population, all the affected are related to each other to some extent. Repeating the regression analysis using relatedness as a covariate (where 126 “related” individuals are relative pairs as far removed as second cousins) still resulted in a highly significant association between repeat length and age at onset (p = 2e-16). b Inverse correlation between CCCTCT repeat length determined in a subgroup of blood and brain samples from deceased XDP patients (blood n = 28; brain n = 40) and AAO. c Inverse correlation between CCCTCT repeat length determined in a subgroup of blood and brain samples from deceased XDP patients (blood n = 28; brain n = 40) and AAD. d Length of the CCCTCT repeat determined in the subgroup of blood and brain samples from deceased XDP patients (blood n = 28; brain n = 40) is not correlated with disease duration (AAD-AAO). AAO, age at onset; AAD, age at death. Brain repeat lengths were determined in cerebellum (n = 39) or occipital cortex (n = 1). In a–d, blood (red dots) and brain samples (blue dots)

The XDP CCCTCT repeat exhibits tissue- and repeat length-dependent somatic expansion

The variation in repeat length between individuals reflects the instability of the CCCTCT repeat in germline transmissions [4, 25]. To gain insight into CCCTCT repeat instability in somatic tissues we have examined repeat length variation in blood (n = 164) and postmortem brain (n = 41) from affected males. In the brain, we analyzed between 1 and 17 brain regions in 41 individuals, including cerebellum only in 17 individuals and occipital cortex only in one individual (Additional File 1: Table S3). The XDP CCCTCT repeat was PCR-amplified using a previously established genotyping assay for repeat sizing [4]. PCR amplification of the repeat results in a distribution of fragment sizes, with repeat length determined as the modal allele in the distribution. Of the 23 postmortem samples in which multiple brain regions were analyzed, 4 (17–012, 17–17, 19–017 and 21–031) exhibited variation by one repeat unit (Additional File 2: Fig. S2) while in 19 individuals the modal repeat length was identical in all brain regions analyzed. Therefore, XDP CCCTCT repeat instability is not substantially reflected in differences in modal repeat length of the repeat-containing PCR amplicons.

We then analyzed XDP CCCTCT instability by quantifying an expansion index from repeat length distributions of GeneMapper outputs of the repeat-containing PCR products [45]. This relatively high throughput method is sensitive to subtle differences in repeat instability that are captured in the majority of alleles. Examples of GeneMapper traces from different tissues are shown in Additional File 2: Fig. S3. The peaks to the left of the modal allele are largely due to PCR slippage, and therefore we quantified only the expansion peaks to the right of the modal allele. These peaks are variable between tissues and are the result of somatic repeat length variation. Expansion indices in blood and brain regions are shown in Fig. 2a, ordered from left to right by the median expansion index per tissue. Very low levels of XDP CCCTCT expansion were detected in blood (median expansion index = 0.19, interquartile range [IQR] = 0.22). In contrast, all brain regions exhibited expansion indices that were significantly greater than those in blood (P < 0.05: Wilcoxon rank-sum tests with Bonferroni correction; Additional File 1: Table S4). Of the brain regions analyzed, cerebellum had the lowest expansion index (median expansion index = 0.77, IQR = 0.32), while occipital cortex exhibited the highest expansion index (median expansion index = 1.59, IQR = 0.7). Replicate PCR amplifications from the same DNA samples demonstrated that differences between brain regions are not due to technical variation (Additional File 2: Fig. S4). Statistically significant differences in expansion indices (P < 0.05: Wilcoxon rank sum tests with Bonferroni correction) were observed between some of the brain regions, most notably in comparisons with cerebellum or occipital cortex (Additional File 1: Table S4). Overall, there appeared to be a tendency towards higher expansion indices in cortical regions (cingulate gyrus, prefrontal cortex (BA9), parietal cortex, insula, temporal pole and occipital cortex) than subcortical areas (cerebellum, caudate, substantia nigra, inferior olivary nucleus, red nucleus, medial thalamus, hippocampus, putamen, lateral thalamus, deep cerebellar nuclei, sub-thalamic nucleus). Of the subcortical structures, there was no obvious distinction in expansion indices between forebrain (caudate, putamen, hippocampus, thalamus, subthalamic nucleus), midbrain (red nucleus) or hindbrain (deep cerebellar nuclei, inferior olivary nucleus) regions, with the exception of cerebellum (Fig. 2a). Due to the considerable variation in repeat expansion between individuals, we further evaluated tissue patterns of expansion by performing hierarchical clustering on a heatmap plot based on scaled expansion index values (Fig. 2b). The heatmap revealed similar patterns of brain region-specific expansion across individuals and distinguished two major clusters comprised of cortical and subcortical brain areas (Fig. 2b).

Fig. 2
figure 2

XDP CCCTCT repeat expansion index in blood and brain regions. a Distribution of expansion indices ranked by median values in blood and brain regions. Box-whisker plots show median ± interquartile range (IQR) and dots show values in individual patient samples b Heatmap of expansion indices values in different individuals (rows), scaled (z-score) across brain tissues (columns). To avoid poor normalization during scaling, brain regions with fewer than 12 measures and individuals with fewer than 6 tissue samples were excluded (Additional File 2: Fig. S4). Brain regions with no measurement are represented as black boxes. Patient IDs are shown with modal cerebellar repeat lengths in parentheses. Note that modal repeat lengths were the same across tissues for an individual with a few exceptions (Additional file 1: Fig.S2). c Linear regression analyses showing relationships between CCCTCT repeat length and expansion index in each tissue. The regression equations shown in bold font highlight those tissues (blood, cerebellum, subthalamic nuclei, cingulate gyrus, temporal pole, occipital cortex) showing a significant association of expansion index with repeat length. Grey shaded areas show 95% confidence interval. Blood (n = 164), Cereb = cerebellum (n = 40), Cd = caudate (n = 17), SbN = substantia nigra (n = 19), ION = inferior olivary nucleus (n = 9), RN = red nucleus (n = 11), Med.Th = medial thalamus (n = 20), Hip = hippocampus (n = 19), Pu = putamen (n = 19), Lat.Th = lateral thalamus (n = 20), DCN = deep cerebellar nuclei (n = 21), STh = subthalamic nucleus (n = 8), Cin.Gyrus = cingulate gyrus (n = 20), BA9 = frontal cortex Brodmann area 9 (n = 21), Parietal Cx = parietal cortex (n = 20), Insula = insular cortex (n = 19), TmP = temporal pole (n = 20), Oc. Cx = occipital cortex (n = 24).

As individuals differ in their repeat length, we investigated the extent to which repeat length might explain the variation in expansion index within any one tissue (Fig. 2c). Overall, the data showed positive correlations between expansion index and repeat length that were statistically significant in a subset of the tissues (blood, cerebellum, subthalamic nuclei, cingulate gyrus, temporal pole, occipital cortex). The proportion of the variation in expansion index explained by repeat length varied from 2% in blood to 45% in the red nucleus. The cerebellum exhibited the most significant correlation (P = 6.8 × 10–6), with repeat length explaining 37% of the expansion index variation. The various strengths of the associations with repeat length likely differ as a function of sample number, the magnitude of the instability, and the cell type heterogeneity in each tissue piece that is sampled. e.g. blood shows minimal repeat expansion, limiting the sensitivity to detect biological variation. In cerebellum, the relatively strong association with repeat length is likely contributed by both cell type homogeneity—99% of all cerebellar neurons are granule cells—and the greater number of cerebellar samples relative to the other brain regions.

Together, these data demonstrate greater somatic expansion of the XDP CCCTCT repeat in the brain than in blood as well as brain region-specific propensities for expansion that are similar across individuals. Significantly, we show that somatic CCCTCT expansion is dependent on repeat length, consistent with a contribution of somatic expansion to the onset of disease.

The XDP CCCTCT repeat exhibits large repeat length changes and expansion-biased instability in the brain

Analysis of repeat instability in fragment sizing data obtained from PCR-amplified “bulk” genomic DNA, as above, is limited by the lack of sensitivity to detect rare alleles and an upper limit for accurate fragment sizing of ~ 330–560 base pairs, equating to ~ 32–70 CCCTCT repeats. Further, while allele length distributions can be quantified in the PCR products, as with the expansion index metric, this may not accurately reflect the distribution of allele lengths present in genomic DNA due to contraction bias inherent to the PCR. Therefore, to investigate more fully the spectrum of repeat length mosaicism in XDP brains we employed two small pool-PCR (SP-PCR) approaches, providing the sensitivity to detect rare somatic events and to quantify allele size distributions. We analyzed a subset of the brain tissues, sampling across regions (occipital cortex, caudate, putamen, cerebellum) exhibiting a range of instabilities as determined from the GeneMapper-based analysis above, and across individuals with a range of repeat lengths (17–17: 54/55 repeats, 19–008: 41 repeats, 18–006: 35 repeats; Fig. 3, Table 1).

Fig. 3
figure 3

Southern blot images and estimated repeat lengths. a Representative Southern blot images for 17–17 occipital cortex (Oc. Cx) (lanes 1–10) and cerebellum (Cereb) (lanes 11–20) illustrating the varying degree of instability across brain regions. Each lane represents PCR amplification of ~ 30 g.e. M1 and M2 size markers are DIG VII and VIII, respectively and are shown with the corresponding base pair lengths. b Estimated CCCTCT repeat lengths based on distance migrated relative to the M1 and M2 markers. Repeat size data for each sample are obtained from 36 replicates (individual lanes), each with ~ 30 g.e. input DNA amount

Table 1 Summary of XDP CCCTCT repeat sizing and instability analyses

We first performed SP-PCR in conjunction with Southern blot detection, diluting the genomic DNA to approximately 30 genome equivalents (g.e) prior to PCR amplification of the CCCTCT repeat. Examples of the Southern blots are shown in Fig. 3a and a summary of the data is provided in Fig. 3b and Table 1, the latter indicating the approximate highest and lowest repeat lengths detectable for each sample. The approximate length ranges for the greatest density of signal on the Southern blots encompassed the repeat sizes determined by standard genotyping (Fig. 3b, Table 1). Notably, all samples showed distinct additional bands reflecting expansions or contractions, with a bias towards expansions. The largest alleles detected across all samples ranged from ~ 77 to ~ 149 repeats with increases in length relative to those determined by standard genotyping ranging from ~ 27 to > 100 units (Table 1). The smallest alleles detected ranged from ~ 22 to ~ 42 repeats, representing ~ 13–24 unit decreases relative to genotyped repeat lengths. The highest and lowest approximate repeat lengths detected were found in 18–006 occipital cortex (149 and 22, respectively) despite this sample having the shortest genotyped repeat and smallest expansion index (Table 1). Among the different tissues from individual 17–17 (54 repeats), occipital cortex exhibited the most instability, with repeats ranging from 36 to 129. Cerebellum was the most stable of these tissues, but nevertheless did show evidence for alleles ranging from 31 to as high as 82 repeats. Caudate and putamen exhibited degrees of mosaicism between those of occipital cortex and cerebellum. In occipital cortex from 19 to 008 (41 repeats) we detected a range of repeat lengths from 31 to 77. In general, qualitative patterns of instability observed on the Southern blots approximately parallel quantitative differences in expansion indices (Table 1, Fig. 3) but highlight the occurrence of rarer somatic events that are not detected in the bulk PCR-based analyses.

While the SP-PCR Southern blot analyses allow detection of large repeat length changes, input DNA amounts of multiple genomes do not allow for quantitative analyses of repeat length distributions in these samples as signals from individual amplification products are not necessarily distinguishable. To quantify repeat length distributions, we therefore performed SP-PCR of single input molecules. We targeted ~ 120–240 individual molecules per sample (Table 1) with the aim of capturing somatic events that occurred at a frequency of ~ 0.5–1%, and sized individual PCR products on the ABI sequencer to achieve single repeat resolution. It should be noted that fragment sizing of SP-PCR products has the same sizing limitations as bulk PCR and thus we were not able to assess the very large rare expansions that were seen on Southern blots. We examined the same brain samples as for the Southern blot-based analyses and extended the single molecule analyses to include putamen, caudate, cerebellum in addition to occipital cortex from 19–008 (Table 1, Fig. 4 and Additional File 1: Table S5). These data revealed a high proportion of alleles with lengths either expanded or contracted relative to the modal repeat length (Fig. 4b). Note that the modal repeat length in the single molecule input SP-PCR data was identical to the repeat length determined by standard genotyping of bulk genomic DNA with the exception of 17–17 cerebellum where SP-PCR modal allele was greater by one repeat (Table 1). Across these samples 65% to 84% (mean 74%) of alleles deviated from the modal allele length. The frequency of expansions ranged from 30 to 58% (mean 49%) while the frequency of contractions was lower overall, ranging from 16 to 45% (mean 26%) (Fig. 4a, Additional File 1: Table S5). The relative frequencies of contracted, modal and expanded alleles differed across the four brain regions of individual 17–17 (Chi2 = 33.30, df = 6, P < 0.0001) with a relatively high proportion of expansions in occipital cortex and a relatively low proportion of expansions in cerebellum. Relative frequencies of contracted, modal and expanded alleles were not significantly different between the four brain regions of individual 19–008 (Chi2 = 8.882, df = 6, P = 0.1803) but differed significantly between occipital cortices of the three individuals (Chi2 = 12.52, df = 4, P = 0.0139). The majority of the expanded alleles were 1–4 repeat units, with expansions of 5 or more repeats occurring in 2%-18% of alleles (mean 9%) and expansions of 20 or more repeats occurring in 0%-12% of alleles (mean 2%) (Additional File 1: Table S5). The majority of the contracted alleles were also in the range of 1–4 units, with contractions of 5 or more repeats in 0–11% of alleles (mean 3%) and contractions of 20 or more repeats in 0–3% of alleles (mean 0.8%) (Additional File 1: Table S5). Overall, the allele size distributions in the single molecule data capture both the tissue-specific and individual-specific differences in instability that are similarly reflected in the expansion index measure and SP-PCR Southern blot analyses.

Fig. 4
figure 4

XDP CCCTCT repeat length distributions in brain regions. Repeat lengths were determined by fragment sizing of amplicons obtained in single molecule input PCRs in four tissues (Pu = putamen, Cd = caudate, Cereb = cerebellum, Oc. Cx = occipital cortex) across three patients. a Percentages of expansions and contractions compared to the modal allele. b Histograms of repeat length frequencies. Data in a and b were derived from 121–243 single amplifiable molecules for each sample. Refer to Table 1 for summary data derived from these analyses

Features of XDP CCCTCT somatic expansion are shared among other microsatellite repeats

To gain additional insight into XDP CCCTCT repeat dynamics we were interested in exploring overlaps with other microsatellite repeats, in particular: 1) a different CCCTCT repeat, and 2) the unstable expanded HTT CAG repeat due to shared genetic and pathological features of HD and XDP. No other disease-causing CCCTCT repeats have been described to date, however CCCTCT repeats are common elements of SVA retrotransposons in the human genome [20]. To identify another CCCTCT repeat to study in comparison to the XDP repeat, we defined inclusion criteria as: 1) the repeat is similar in length to the XDP repeat (~ 35–50) and 2) the repeat-containing SVA is located in an intron and inserted in reverse orientation relative to the gene transcript, as it is for the XDP SVA. We thus identified a CCCTCT repeat of 39 units in the reference genome (hg19 chr18:47,105,372–47,105,605) within an SVA inserted in reverse orientation in intron 5 of the endothelial lipase G gene (LIPG), hereafter referred to as the LIPG CCCTCT repeat. We first PCR-amplified the LIPG CCCTCT repeat from a subset of XDP patient cerebellar DNAs. Repeat length varied from 39 to 71 (median = 53, IQR = 10), with two repeat lengths distinguishable in some individuals and only one in others (Additional File 1: Table S6). We then identified six individuals for analyses of LIPG CCCTCT repeat instability across brain regions (cerebellum, caudate, hippocampus, BA9, temporal pole and occipital cortex) that exhibited a range of XDP CCCTCT expansion levels. The six individuals were selected based both on tissue availability and having two LIPG CCCTCT repeat lengths sufficiently well-separated to allow quantification of an expansion index from each allele (Additional File 1: Table S6). Examples of GeneMapper outputs of LIPG CCCTCT repeat-containing PCR products are shown in Additional File 2: Fig.S5. Quantification of an expansion index across all the brain samples (Fig. 5a) revealed the lowest expansion index in cerebellum (median = 0.12, IQR = 0.06), and the highest expansion index in caudate (median = 0.95, IQR = 0.5), with significantly lower cerebellar expansion indices relative to other brain regions (P < 0.05: Wilcoxon rank-sum tests with Bonferroni correction, Additional File 1: Table S7). A comparison of LIPG and XDP CCCTCT expansion indices (Fig. 5a) revealed significantly lower values for cerebellum, hippocampus, BA9, temporal pole and occipital cortex brain regions despite the LIPG having longer repeats on average than the XDP repeat (P < 0.05: Wilcoxon rank-sum test, Additional File 1: Table S7). LIPG CCCTCT expansion indices also positively correlated with repeat length with the proportion of the variation in expansion index explained by repeat length varying from 54% in caudate to 69% in the temporal pole (Fig. 5b). It is worth noting that the variability in expansion index as a function of repeat length (R2) may be overestimated in these data due to the inclusion of two alleles from the same individual. Overall, despite the small sample size and lower absolute levels of expansion of the LIPG repeat compared to the XDP repeat, these data reveal that both repeats share properties of length-dependent expansion being relatively low in cerebellum.

Fig. 5
figure 5

Expansion of the CCCTCT LIPG repeat and HTT CAG repeat in comparison to the XDP CCCTCT repeat. a Distribution of expansion indices of LIPG and XDP CCCTCT repeats in XDP postmortem brain tissues. Box-whisker plots show median ± IQR and dots show values for individual alleles. XDP repeat: data are the same as in Fig. 2. LIPG repeat: Cereb = Cerebellum (n = 23 individuals, 40 alleles), Cd = Caudate (n = 5 patients, 10 alleles), Hip = Hippocampus (n = 6 patients, 10 alleles), BA9 = frontal cortex Brodmann area 9 (n = 6 patients, 11 alleles), TmP = Temporal pole (n = 6 patients, 12 alleles), Oc. Cx = Occipital cortex (n = 5 patients, 10 alleles). Note that some alleles that failed QC were excluded. b Linear regression analyses showing relationships between LIPG CCCCTC repeat length and expansion index in each brain region. The regression equations shown in bold font highlight those tissues (caudate, BA9, temporal pole and occipital cortex) showing a significant association of expansion index with repeat length. Grey shaded areas show 95% confidence interval. c Correlation of mean HTT CAG expansion index in three HD individuals (Materials and Methods) and mean XDP CCCTCT expansion indices. Refer to Additional File1: Table S3 and Table S6 for sample numbers for each brain region for XDP and LIPG respectively

We previously reported, using similar quantitative analyses, tissue-specific patterns of somatic expansion of the HTT CAG repeat in HD postmortem brains [28]. To compare tissue-specific instability of the XDP CCCTCT and HTT CAG repeats, we plotted mean expansion indices across all patient samples for nine brain regions (BA9, cerebellum, hippocampal formation, temporal pole, putamen, occipital cortex, subthalamic nuclei and caudate) that were shared across the HD study and this XDP study (see materials and methods). We found that XDP and HTT repeat expansion indices in XDP and HD patient brain tissues, respectively, were highly correlated (correlation coefficient r = 0.65, P = 0.0057, Fig. 5c), indicating shared tissue-specific expansion propensities of these two different disease-associated repeats. In contrast, and as indicated in Fig. 5a, the XDP and LIPG expansion indices are not correlated across the six tissues analyzed (correlation coefficient r = 0.14, P = 0.79).

Discussion

Previous studies have shown that the length of the XDP-associated CCCTCT repeat in blood is inversely correlated with AAO, accounting for ~ 50% of the AAO variance [4, 25]. There is also evidence for correlations between repeat length and other clinical disease measures [25]. The present study supports and extends these data; in our expanded blood dataset (N = 266), repeat length accounted for ~ 46% of the variance in AAO, and in as few as 40 individuals we detected a significant correlation between AAO and repeat length measured in brain DNA (R2 = 0.55). The different R2 values between the various studies and our cohorts [25, 47] may in part be explained by differences in the accuracy in determining AAO and warrants additional investigation. In addition, in a subset of individuals with known AAO and AAD we show for the first time that repeat length is inversely correlated with AAD, with a relationship paralleling that between repeat length and AAO. In contrast, we observed no significant correlation between repeat length and disease duration (the time between onset and death), an observation previously reported in HD [48]. However, this does not preclude a possible stronger effect of repeat length on duration that is counterbalanced by an effect of AAO [49] (i.e. longer repeat length resulting in a shorter duration, counterbalanced by longer repeat length resulting in earlier AAO and subsequent longer duration). As AAD was only available for 68 individuals in this study, additional patient data will be needed for further dissection of repeat length-dependent relationships with disease duration. Importantly, our data underscore the importance of CCCTCT repeat length in driving the rate of XDP, motivating the investigation of the instability of this repeat tract in somatic cells in patients.

To gain insight into the somatic instability of the XDP-associated CCCTCT repeat we have used multiple methodologies, including single molecule-based analyses, to probe the spectrum of repeat length mosaicism in blood and across seventeen brain regions from XDP patients. We demonstrate that the XDP CCCTCT repeat exhibits extensive somatic mosaicism, notably length-dependent and tissue-specific expansion that is measurable in the bulk of alleles, and the presence of rarer alleles in the brain that can be either substantially contracted or expanded on the order of ~ 10 s > 100 repeats relative to the repeat length determined using standard genotyping. Given the inverse correlation of CCCTCT repeat length with AAO these observations implicate somatic expansion as a driver of the rate of onset of XDP. Notably, a GWAS identified two genes, MSH3 and PMS2 as modifiers of the age of onset of XDP [36]. These genes are also modifiers of HD age at onset [30] and encode DNA mismatch repair proteins that modulate the somatic instability of disease-associated trinucleotide repeats, including the HTT CAG repeat [30, 37,38,39]. It is likely, therefore, that MSH3 and PMS2 modify XDP onset by altering the rate of somatic CCCTCT expansion.

We find greater levels of somatic expansion in all brain regions analyzed relative to levels in blood, supporting recent observations in two XDP patients [41]. Within the brain, we observe region-specific differences in the degree of repeat expansion that are reflected across the different XDP individuals, with cerebellum exhibiting the most stability and cortical structures tending to be the most unstable. Several other disease-associated microsatellite repeats are relatively stable in cerebellum [28, 50, 51]. Here, we show substantial correlation between brain region-specific levels of expansion of the XDP CCCTCT repeat and the HTT CAG repeat, as previously observed in a similar comparison between expansion of the HTT CAG repeat and of the ATXN1 CAG repeat underlying spinocerebellar ataxia type 1 (SCA1) [28]. These data provide support for common proteins (trans-acting factors) that modify tissue-specific levels of somatic expansion of both the XDP and HTT repeats, as well as other disease-associated repeats. We also found that the SVA-associated CCCTCT repeat within the LIPG gene, not known to be associated with any disease, exhibited repeat length-dependent expansion that was low in cerebellum compared to other brain regions analyzed. Interestingly, the LIPG CCCTCT repeat exhibited less instability than the XDP CCCTCT repeat in most of the brain regions analyzed, despite its relatively longer repeat lengths, pointing to potential modification of CCCTCT repeat instability in cis. It is plausible that local chromatin structure at the TAF1 SVA locus might predispose the CCCTCT repeat to expand, while at the LIPG SVA locus, expansion of the CCCTCT repeat is comparatively suppressed. In line with this idea, disease-associated short tandem repeats were found to be enriched at 3D chromatin boundaries; in contrast, matched non-disease-associated repeats did not exhibit such an enrichment [52]. Thus, insights into chromatin structural features at the XDP SVA locus relative to other non-disease associated SVAs may provide clues to the instability propensity of its CCCTCT repeat tract. We previously reported the high G-quadruplex-forming potential of the reverse orientation AGAGGG repeat in the XDP SVA sequence [4]; whether this plays a role in its repeat instability remains to be investigated. The LIPG gene is also expressed at low levels in brain tissues. Transcription has been proposed to play a role in promoting repeat instability [53], and therefore a low rate of transcription through the LIPG gene may contribute to the lower level of instability of the SVA-associated CCCTCT repeat within this gene.

Although the instability of the XDP repeat is biased towards expansions, SP-PCR revealed a substantial proportion of XDP repeat contractions. With the caveat of a small sample number of XDP brains analyzed by SP-PCR, the XDP repeat appears on the face of it to exhibit more somatic contraction than the HTT CAG repeat in HD postmortem brain in which similar SP-PCR studies have been conducted [29] Of note, the XDP repeat exhibits a greater frequency of contractions relative to expansions in transmission from males [4, 54] that contrasts with a strong bias in expansion frequency in male transmissions of the HTT CAG repeat [55] Therefore, despite similarities in tissue-specific expansion propensities of these two repeats, there appear to be differences in their dynamics that warrant further investigation.

Our analyses of repeat instability in different brain tissues do not immediately point to any clear correlation with brain regions implicated either through neuropathological or neuroimaging studies to be susceptible in XDP [10,11,12,13,14,15, 36]. e.g. neuropathological changes have been described in tissues that include caudate, putamen, cortex and cerebellum [10], yet these regions encompass both the lowest (cerebellum) and highest (cortex) levels of expansion. However, the association of repeat instability and cellular vulnerability is currently challenging due to: 1) the lack of cell type-specific resolution of repeat instability; 2) limited XDP neuropathology data; 3) neurodegeneration, notably of MSNs [6,7,8]. In HD, GWAS studies have provided support for a two-step model of pathogenesis that depends both on the rate of somatic CAG expansion and repeat length threshold(s) needed to trigger a toxic process(es) [30]. Both the rate of repeat expansion and toxicity-eliciting threshold may differ by cell type, and as both instability and toxicity components are needed for pathogenesis, high levels of expansion do not necessarily predict cellular vulnerability, e.g. this provides a logical explanation for high levels of HTT CAG expansion seen in the liver yet the absence of obvious liver pathology [28]. A two-step model provides a framework for other repeat expansion diseases, and similarly can explain why the striatum is not the primary target of pathogenesis in SCA1 despite high levels of CAG expansion in that tissue [28]. We propose that this model can also be applied to XDP, predicting that somatic expansion of the CCCTCT repeat in certain cell types will elicit a toxic process(es) ultimately culminating in clinical disease. A full understanding of XDP pathogenesis will therefore entail dissecting both instability and toxicity components in specific cell types. Further, there is evidence for altered brain connectivity in XDP [12, 14, 15, 56, 57], providing added complexity such that repeat expansion in one cell-type may trigger functional deficits at the level of a neuronal circuit. Notably our results provide evidence for a landscape of somatic events that include both repeat expansions and contractions, highlighting the importance of cell type-specific level resolution to understand relationships with disease processes. e.g. might repeat contraction in certain cell types confer protection? Currently, the nature of the toxic species is unclear, with reduced TAF1 levels and novel TAF1 isoforms being plausible candidates [18, 19, 21,22,23,24]. The inverse correlation of TAF1 mRNA levels with CCCTCT repeat length seen in blood [25] is consistent with a role of TAF1 levels in a pathological process triggered by CCCTCT repeat expansion. In the brain, a neuronal-specific TAF1 isoform (nTAF1), generated as a result of inclusion of microexon 34’, was reduced in caudate, nucleus accumbens and cortex from a single XDP individual [18]. However, a more recent study did not find evidence for altered splicing of microexon 34’ in additional XDP postmortem brains [58]. Further studies in postmortem brain, including cell type-specific analyses, will be needed to gain a more comprehensive understanding of the effects of the XDP mutation on TAF1 expression and splicing. Finally, while the identification of the MSH3 and PMS2 genes as XDP onset modifiers provides strong support for repeat expansion as the upstream driver of a toxic process(es), TAF1 itself has been implicated in promoting genome integrity [59,60,61]. Therefore, it is possible that altered TAF1 function in the disease process may further impact the DNA repair processes that underlie repeat instability, e.g. our data hint at differences in the instability of a non-disease-associated repeat (LIPG CCCTCT tract) in some tissues. Thus, genome-wide analyses of DNA instability/integrity in XDP patient brain would be of interest.

The prediction from our data is that that somatic CCCTCT repeat expansion contributes to length-dependent clinical measures, such as AAO. Of note, in the current dataset we find no difference in somatic expansion measured in blood between patients reporting symptom onset as either being predominantly dystonia or parkinsonism (Additional File 2: Fig.S1). This is consistent with the similar relationship between repeat length and AAO in these two patient subsets (Additional File 2: Fig.S1). Larger sample numbers will be needed to provide sufficient power for further tests of associations of repeat instability with clinical endpoints such as AAO. In addition, further studies will be needed to understand the relationship between repeat instability in blood and brain that will inform tests of association of instability with clinical measures.

Conclusions

These data demonstrate that the XDP CCCTCT repeat is unstable in somatic cells, exhibiting properties that are consistent with a role for somatic expansion in determining the timing of disease onset. Our data suggest further avenues of investigation aimed at understanding the dynamics of this repeat mutation and relationship to pathogenesis.

Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding authors on reasonable request. Requests for tissue specimens may be directed to xdp@partners.org.

Change history

Abbreviations

XDP:

X-linked dystonia-parkinsonism

SVA:

SINE-VNTR-Alu

HD:

Huntington’s disease

TAF1:

TATA-binding-protein (TBP)-associated factor-1

TBP:

TATA-binding-protein

LIPG:

Lipase G, Endothelial Type

MSNs:

Medium-spiny neurons

AAO:

Age at Onset

AAD:

Age at death

GWAS:

Genome-wide association

IRBs:

MGH: Massachusetts General Hospital

IRB:

Institutional Review Board

gDNA:

Genomic DNA

CCXDP:

Collaborative Center for XDP

DIG:

Digoxygenin

BA9:

Frontal cortex Brodmann area 9

SbN:

Substantia nigra

ION:

Inferior olivary nucleus

RN:

Red nucleus

DCN:

Deep cerebellar nuclei

STh:

Subthalamic nucleus

SP-PCR:

Small pool-PCR

HTT:

Huntingtin

MSH3:

MutS homolog 3: PMS2: PMS1 homolog2, Mistmatch repair system component

SCA1:

Spinocerebellar ataxia type 1

References

  1. Lee LV, Rivera C, Teleg RA, Dantes MB, Pasco PMD, Jamora RDG et al (2011) The unique phenomenology of sex-linked dystonia parkinsonism (XDP, DYT3, “Lubag”). Int J Neurosci 121:3–11

    Article  PubMed  Google Scholar 

  2. Lee LV, Maranon E, Demaisip C, Peralta O, Borres-Icasiano R, Arancillo J et al (2002) The natural history of sex-linked recessive dystonia parkinsonism of Panay, Philippines (XDP)*. Park Relat Disord 9:29–38

    Article  Google Scholar 

  3. Lee LV, Kupke KG, Caballar-Gonzaga F, Hebron-Ortiz M, Müller U (1991) The phenotype of the X-linked dystonia-parkinsonism syndrome: an assessment of 42 cases in the Philippines. Med (United States) 70:179–187

    CAS  Google Scholar 

  4. Bragg DC, Mangkalaphiban K, Vaine CA, Kulkarni NJ, Shin D, Yadav R et al (2017) Disease onset in X-linked dystonia-parkinsonism correlates with expansion of a hexameric repeat within an SVA retrotransposon in TAF1. Proc Natl Acad Sci U S A 114:E11020–E11028

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Domingo A, Westenberger A, Lee LV, Brænne I, Liu T, Vater I et al (2015) New insights into the genetics of X-linked dystonia-parkinsonism (XDP, DYT3). Eur J Hum Genet 23:1334–1340

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Waters CH, Faust PL, Powers J, Vinters H, Moskowitz C, Nygaard T et al (1993) Neuropathology of lubag (x-linked dystonia parkinsonism). Mov Disord 8:387–390

    Article  CAS  PubMed  Google Scholar 

  7. Goto S, Lee LV, Munoz EL, Tooyama I, Tamiya G, Makino S et al (2005) Functional anatomy of the basal ganglia in X-linked recessive dystonia-parkinsonism. Ann Neurol 58:7–17

    Article  PubMed  Google Scholar 

  8. Goto S, Kawarai T, Morigaki R, Okita S, Koizumi H, Nagahiro S et al (2013) Defects in the striatal neuropeptide y system in X-linked dystonia-parkinsonism. Brain 136:1555–1567

    Article  PubMed  Google Scholar 

  9. Vonsattel JP, Myers RH, Stevens TJ, Ferrante RJ, Bird ED, Richardson EP (1985) Neuropathological classification of huntington’s disease. J Neuropathol Exp Neurol 44:559–577

    Article  CAS  PubMed  Google Scholar 

  10. Arasaratnam CJ, Singh-Bains MK, Waldvogel HJ, Faull RLM (2021) Neuroimaging and neuropathology studies of X-linked dystonia parkinsonism. Neurobiol Dis 48:105186. https://doi.org/10.1016/j.nbd.2020.105186

    Article  CAS  Google Scholar 

  11. Petrozziello T, Mills AN, Vaine CA, Penney EB, Fernandez-Cerado C, Legarda GPA et al (2020) Neuroinflammation and histone H3 citrullination are increased in X-linked Dystonia Parkinsonism post-mortem prefrontal cortex. Neurobiol Dis. 144:05032. https://doi.org/10.1016/j.nbd.2020.105032

    Article  CAS  Google Scholar 

  12. Brüggemann N, Heldmann M, Klein C, Domingo A, Rasche D, Tronnier V et al (2016) Neuroanatomical changes extend beyond striatal atrophy in X-linked dystonia parkinsonism. Park Relat Disord 31:91–97

    Article  Google Scholar 

  13. Brüggemann N, Rosales RL, Waugh JL, Blood AJ, Domingo A, Heldmann M et al (2017) Striatal dysfunction in X-linked dystonia-parkinsonism is associated with disease progression. Eur J Neurol 24:680–686

    Article  PubMed  Google Scholar 

  14. Blood AJ, Waugh JL, Münte TF, Heldmann M, Domingo A, Klein C et al (2018) Increased insula-putamen connectivity in X-linked dystonia-parkinsonism. NeuroImage Clin 17:835–846

    Article  PubMed  Google Scholar 

  15. Hanssen H, Heldmann M, Prasuhn J, Tronnier V, Rasche D, Diesta CC et al (2018) Basal ganglia and cerebellar pathology in X-linked dystonia-parkinsonism. Brain 141:2995–3008

    Article  PubMed  Google Scholar 

  16. Németh AH, Nolte D, Dunne E, Niemann S, Kostrzewa M, Peters U et al (1999) Refined linkage disequilibrium and physical mapping of the gene locus for X-linked Dystonia-Parkinsonism (DYT3). Genomics 60:320–329

    Article  PubMed  Google Scholar 

  17. Nolte D, Niemann S, Müller U (2003) Specific sequence changes in multiple transcript system DYT3 are associated with X-linked dystonia parkinsonism. Proc Natl Acad Sci U S A 100:10347–10352

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Makino S, Kaji R, Ando S, Tomizawa M, Yasuno K, Goto S et al (2007) Reduced neuron-specific expression of the TAF1 gene is associated with X-linked dystonia-parkinsonism. Am J Hum Genet 80:393–406

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Aneichyk T, Hendriks WT, Yadav R, Shin D, Gao D, Vaine CA et al (2018) Dissecting the causal mechanism of X-linked dystonia-parkinsonism by integrating genome and transcriptome assembly. Cell 172:897-909.e21

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Hancks DC, Kazazian HH (2010) SVA retrotransposons: evolution and genetic instability. Semin Cancer Biol 20:234–245

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Ito N, Hendriks WT, Dhakal J, Vaine CA, Liu C, Shin D et al (2016) Decreased N-TAF1 expression in X-linked dystonia-parkinsonism patient-specific neural stem cells. DMM Dis Model Mech 9:451–462

    CAS  PubMed  Google Scholar 

  22. Rakovic A, Domingo A, Grütz K, Kulikovskaja L, Capetian P, Cowley SA et al (2018) Genome editing in induced pluripotent stem cells rescues TAF1 levels in X-linked dystonia-parkinsonism. Mov Disord 33:1108–1118

    Article  CAS  PubMed  Google Scholar 

  23. Domingo A, Amar D, Grütz K, Lee LV, Rosales R, Brüggemann N et al (2016) Evidence of TAF1 dysfunction in peripheral models of X-linked dystonia-parkinsonism. Cell Mol Life Sci 73:3205–3215

    Article  CAS  PubMed  Google Scholar 

  24. Al Ali J, Vaine CA, Shah S, Campion L, Hakoum A, Supnet ML et al (2021) TAF1 transcripts and neurofilament light chain as biomarkers for X-linked Dystonia-Parkinsonism. Mov Disord 36:206–215

    Article  CAS  PubMed  Google Scholar 

  25. Westenberger A, Reyes CJ, Saranza G, Dobricic V, Hanssen H, Domingo A et al (2019) A hexanucleotide repeat modifies expressivity of X-linked dystonia parkinsonism. Ann Neurol 85:812–822

    Article  CAS  PubMed  Google Scholar 

  26. Depienne C, Mandel JL (2021) 30 years of repeat expansion disorders: What have we learned and what are the remaining challenges? Am J Hum Genet 108:764–785

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Monckton DG (2021) The contribution of somatic expansion of the CAG repeat to symptomatic development in huntington’s disease: A historical perspective. J Huntingtons Dis 10:7–33

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Pinto RM, Arning L, Giordano JV, Razghandi P, Andrew MA, Gillis T et al (2020) Patterns of CAG repeat instability in the central nervous system and periphery in Huntington’s disease and in spinocerebellar ataxia type 1. Hum Mol Genet 29:2551–2567

    Article  CAS  Google Scholar 

  29. Swami M, Hendricks AE, Gillis T, Massood T, Mysore J, Myers RH et al (2009) Somatic expansion of the Huntington’s disease CAG repeat in the brain is associated with an earlier age of disease onset. Hum Mol Genet 18:3039–3047

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Lee JM, Correia K, Loupe J, Kim KH, Barker D, Hong EP et al (2019) CAG repeat not polyglutamine length determines timing of huntington’s disease onset. Cell 178:887-900.e14

    Article  CAS  Google Scholar 

  31. Hong EP, MacDonald ME, Wheeler VC, Jones L, Holmans P, Orth M et al (2021) Huntington’s disease pathogenesis: two sequential components. J Huntingtons Dis 10:35–51

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Cumming SA, Hamilton MJ, Robb Y, Gregory H, McWilliam C, Cooper A et al (2018) De novo repeat interruptions are associated with reduced somatic instability and mild or absent clinical features in myotonic dystrophy type 1. Eur J Hum Genet 26:1635–1647

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Zhao X, Kumari D, Miller CJ, Kim GY, Hayward B, Vitalo AG et al (2021) Modifiers of somatic repeat instability in mouse models of friedreich ataxia and the fragile x-related disorders: implications for the mechanism of somatic expansion in huntington’s disease. J Huntingtons Dis 10:149–163

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Bettencourt C, Hensman-Moss D, Flower M, Wiethoff S, Brice A, Goizet C et al (2016) DNA repair pathways underlie a common genetic mechanism modulating onset in polyglutamine diseases. Ann Neurol 79:983–990

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Morales F, Vásquez M, Corrales E, Vindas-Smith R, Santamaría-Ulloa C, Zhang B et al (2020) Longitudinal increases in somatic mosaicism of the expanded CTG repeat in myotonic dystrophy type 1 are associated with variation in age-at-onset. Hum Mol Genet 29:2496–2507

    Article  CAS  PubMed  Google Scholar 

  36. Laabs BH, Klein C, Pozojevic J, Domingo A, Brüggemann N, Grütz K et al (2021) Identifying genetic modifiers of age-associated penetrance in X-linked dystonia-parkinsonism. Nat Commun 12:3216. https://doi.org/10.1038/s41467-021-23491-4

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Ciosi M, Maxwell A, Cumming SA, Hensman Moss DJ, Alshammari AM, Flower MD et al (2019) A genetic association study of glutamine-encoding DNA sequence structures, somatic CAG expansion, and DNA repair gene variants, with Huntington disease clinical outcomes. EBioMedicine 48:568–580

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Gomes-Pereira M, Fortune MT, Ingram L, McAbney JP, Monckton DG (2004) Pms2 is a genetic enhancer of trinucleotide CAG-CTG repeat somatic mosaicism: Implications for the mechanism of triplet repeat expansion. Hum Mol Genet 13:1815–1825

    Article  CAS  PubMed  Google Scholar 

  39. Dragileva E, Hendricks A, Teed A, Gillis T, Lopez ET, Friedberg EC et al (2009) Intergenerational and striatal CAG repeat instability in Huntington’s disease knock-in mice involve different DNA repair genes. Neurobiol Dis 33:37–47

    Article  CAS  PubMed  Google Scholar 

  40. Moss DJH, Tabrizi SJ, Mead S, Lo K, Pardiñas AF, Holmans P et al (2017) Identification of genetic variants associated with Huntington’s disease progression: a genome-wide association study. Lancet Neurol 16:701–711

    Article  CAS  PubMed  Google Scholar 

  41. Reyes CJ, Laabs B-H, Schaake S, Lüth T, Ardicoglu R, Rakovic A et al (2021) Brain regional differences in hexanucleotide repeat length in X-linked dystonia-parkinsonism using nanopore sequencing. Neurol Genet. 7:e608

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  42. Albanese A, Del Sorbo F, Comella C, Jinnah HA, Mink JW, Post B et al (2013) Dystonia rating scales: critique and recommendations. Mov Disord 28:874–883

    Article  PubMed  PubMed Central  Google Scholar 

  43. Shiihashi G, Ito D, Yagi T, Nihei Y, Ebine T, Suzuki N (2016) Mislocated FUS is sufficient for gain-of-toxic-function amyotrophic lateral sclerosis phenotypes in mice. Brain 139:2380–2394

    Article  PubMed  Google Scholar 

  44. Fernandez-Cerado C, Legarda GP, Velasco-Andrada MS, Aguil A, Ganza-Bautista NG, Lagarde JBB et al (2021) Promise and challenges of dystonia brain banking: establishing a human tissue repository for studies of X-Linked Dystonia-Parkinsonism. J Neural Transm 128:575–587

    Article  CAS  PubMed  Google Scholar 

  45. Lee JM, Zhang J, Su AI, Walker JR, Wiltshire T, Kang K et al (2010) A novel approach to investigate tissue-specific trinucleotide repeat instability. BMC Syst Biol 4:29. https://doi.org/10.1186/1752-0509-4-29

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Gomes-Pereira M, Bidichandani SI, Monckton DG (2004) Analysis of unstable triplet repeats using small-pool polymerase chain reaction. Methods Mol Biol 277:61–76

    CAS  PubMed  Google Scholar 

  47. Bragg DC, Sharma N, Ozelius LJ (2019) X-linked Dystonia-Parkinsonism: recent advances. Curr Opin Neurol 32:604–609

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Keum JW, Shin A, Gillis T, Mysore JS, Abu Elneel K, Lucente D et al (2016) The HTT CAG-expansion mutation determines age at death but not disease duration in huntington disease. Am J Hum Genet 98:287–298

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Langbehn DR. Longer CAG repeat length is associated with shorter survival after disease onset in Huntington disease. Am J Hum Genet [Internet]. Cell Press; 2022 [cited 2022 Jan 19];109:172–9. Available from: https://linkinghub.elsevier.com/retrieve/pii/S0002929721004572

  50. Takano H, Onodera O, Takahashi H, Igarashi S, Yamada M, Oyake M et al (1996) Somatic mosaicism of expanded CAG repeats in brains of patients with dentatorubral-pallidoluysian atrophy: cellular population-dependent dynamics of mitotic instability. Am J Hum Genet 58:1212–1222

    CAS  PubMed  PubMed Central  Google Scholar 

  51. Cancel G, Gourfinkel-An I, Stevanin G, Didierjean O, Abbas N, Hirsch E et al (1998) Somatic mosaicism of the CAG repeat expansion in spinocerebellar ataxia type 3/Machado-Joseph disease. Hum Mutat 11:23–27

    Article  CAS  PubMed  Google Scholar 

  52. Sun JH, Zhou L, Emerson DJ, Phyo SA, Titus KR, Gong W et al (2018) Disease-associated short tandem repeats co-localize with chromatin domain boundaries. Cell 175:224-238.e15

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Lin Y, Dion V, Wilson JH (2006) Transcription promotes contraction of CAG repeat tracts in human cells. Nat Struct Mol Biol 13:179–180

    Article  CAS  PubMed  Google Scholar 

  54. Westenberger A, Rosales RL, Heinitz S, Freimann K, Lee LV, Jamora RD et al (2013) X-linked Dystonia-Parkinsonism manifesting in a female patient due to atypical turner syndrome. Mov Disord 28:675–678

    Article  PubMed  Google Scholar 

  55. Wheeler VC, Persichetti F, McNeil SM, Mysore JS, Mysore SS, MacDonald ME et al (2007) Factors associated with HD CAG repeat instability in Huntington disease. J Med Genet 44:695–701

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  56. Brüggemann N, Domingo A, Rasche D, Moll CKE, Rosales RL, Jamora RDG et al (2019) Association of pallidal neurostimulation and outcome predictors with X-linked dystonia parkinsonism. JAMA Neurol 76:211–216

    Article  PubMed  Google Scholar 

  57. Sprenger A, Hanssen H, Hagedorn I, Prasuhn J, Rosales RL, Jamora RDG et al (2019) Eye movement deficits in X-linked dystonia-parkinsonism are related to striatal degeneration. Park Relat Disord 61:170–178

    Article  Google Scholar 

  58. Capponi S, Stöffler N, Penney EB, Grütz K, Nizamuddin S, Vermunt MW et al (2021) Dissection of TAF1 neuronal splicing and implications for neurodegeneration in X-linked dystonia-parkinsonism. Brain Commun 3:fcab253. https://doi.org/10.1093/braincomms/fcab253

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  59. Kim JJ, Lee SY, Gong F, Battenhouse AM, Boutz DR, Bashyal A et al (2019) Systematic bromodomain protein screens identify homologous recombination and R-loop suppression pathways involved in genome integrity. Genes Dev 33:1751–1774

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  60. Peng H, Zhang S, Peng Y, Zhu S, Zhao X, Zhao X et al (2021) Yeast bromodomain factor 1 and its human homolog TAF1 play conserved roles in promoting homologous recombination. Adv Sci. 8:e2100753. https://doi.org/10.1002/advs.202100753

    Article  CAS  Google Scholar 

  61. Buchmann AM, Skaar JR, DeCaprio JA (2004) Activation of a DNA damage checkpoint response in a TAF1-defective cell line. Mol Cell Biol 24:5332–5339

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

We thank Marcy MacDonald, Jim Gusella and Jong-Min Lee for helpful input and discussion. We also thank XDP patients and their families for their invaluable contributions to this research.

Funding

This work is supported by funding from the CCXDP to VCW and LJO. RY is supported by the Massachusetts General Hospital Fund for Medical Discovery.

Author information

Authors and Affiliations

Authors

Contributions

Sample acquisition: EBP, MGM, CF-C, MSV-A, GPL, NGG-B, JBBL, PJA, TM-B, GA, MLS, JKD, CG, NS, ELM, MCA, CCED, DCB. Sample processing: EBP, MGM, LNC. Data acquisition: LNC, TG, AMM. Data analysis: LNC, AMM, RY, KC, TG, LJO, VCW. Data interpretation: LNC, AMM, LJO, VCW. Study conception: LJO, VCW. Design of the work: LNC, AMM, LJO, VCW. Manuscript drafting and revision: LNC, AMM, LJO, VCW. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Laurie J. Ozelius or Vanessa C. Wheeler.

Ethics declarations

Ethics for approval and consent to participate

All participants provided written informed consent, and the study was approved by Massachusetts General Hospital (Boston, MA, USA) and Jose R. Reyes Memorial Medical Center (Manila, Philippines) Institutional Review Boards (IRBs). Post-mortem brain tissue from XDP patients was obtained in collaboration with the Collaborative Center for XDP (CCXDP), at Massachusetts General Hospital (Boston, MA, USA), Makati Medical Center (Makati City, Philippines), and the Sunshine Care Foundation (Panay, Philippines). All procedures related to the collection, processing, and use of XDP patient post-mortem brain tissues were approved by IRBs at Makati Medical Center (Makati City, Philippines) and Massachusetts General Hospital (Boston, MA, USA).

Consent for publication

All authors consented to the publication of the manuscript.

Competing interests

V.C.W. is a scientific advisory board member of Triplet Therapeutics, Inc., a company developing new therapeutic approaches to address triplet repeat disorders such Huntington’s disease and Myotonic Dystrophy. Her financial interests in Triplet Therapeutics were reviewed and are managed by Massachusetts General Hospital and Mass General Brigham in accordance with their conflict of interest policies. She is a scientific advisory board member of LoQus23 Therapeutics, Ltd and has provided paid consulting services to Alnylam, Inc., Acadia Pharmaceuticals and Biogen, Inc. She has also received research support from Pfizer Inc. L.J.O. receives royalties from Athena Diagnostics.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

The original version of this article was revised: an error in the name of the second author in the author group was corrected.

Supplementary Information

Additional file 1

. Tables S1–S7.

Additional file 2

. Figures S1–S5.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Campion, L.N., Mejia Maza, A., Yadav, R. et al. Tissue-specific and repeat length-dependent somatic instability of the X-linked dystonia parkinsonism-associated CCCTCT repeat. acta neuropathol commun 10, 49 (2022). https://doi.org/10.1186/s40478-022-01349-0

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s40478-022-01349-0