Propensity for somatic expansion increases over the course of life in Huntington disease

  1. Radhia Kacher
  2. François-Xavier Lejeune
  3. Sandrine Noël
  4. Cécile Cazeneuve
  5. Alexis Brice
  6. Sandrine Humbert  Is a corresponding author
  7. Alexandra Durr  Is a corresponding author
  1. Sorbonne Université, Paris Brain Institute (ICM Institut du Cerveau), AP-HP, INSERM, CNRS, University Hospital Pitié-Salpêtrière, France
  2. Univ. Grenoble Alpes, INSERM, U 1216, Grenoble Institut Neurosciences, France
  3. Paris Brain Institute’s Data and Analysis Core, University Hospital Pitié-Salpêtrière, France
  4. Neurogenetics Laboratory, Department of Genetics, Assistance Publique–Hôpitaux de Paris, University Hospital Pitié-Salpêtrière, France
6 figures, 6 tables and 1 additional file

Figures

Figure 1 with 1 supplement
Somatic instability is negligible during gestation but increases with age.

(A) Left panel: Changes in expansion index (EI) for each Huntington disease (HD) individual, across two or sometimes three visits over multiple years (see Table 5). Different colors indicate different reference (CAG)n at diagnosis (main peak on PCR profile as shown in Figure 1—figure supplement 1). For each value, the disease status is indicated with an empty triangle (premanifest) or filled triangle (manifest; score >5 on the Unified HD Rating Scale total motor score [UHDRS-TMS]). Note that the measurement of the inherited allele remained stable across visits except for patient 44 (Table 2), who went from 54 to 55 (CAG)n at the second sampling, which was taken into account for the EI. Right panel: A closer look at the values clustered at the bottom of the axis in A ((CAG)n 39–46). The EI increases with progression to manifest state even for individuals with relatively small reference repeats. (B) Scatter plot and regression lines show the linear relation between EI and (CAG)n as observed in adult blood from the longitudinal cohort (left panel, n = 50 patients with at least two samples each; of these, 12 had a third sample; yellow: first sample, orange: second sample, and red triangles: third sample) and cortical tissue (right panel, n = 7 fetal brains, green triangles, and n = 15 adult brains, blue triangles). Pearson’s correlation coefficients and estimated regression slopes with p-values, indicated in the upper portion of each graph, reveal a positive linear relation between EI and reference CAG repeat length. (C) EI values from seven fetal samples according to the reference CAG repeat, ranging from 40 to 46 (at 13 weeks gestation). The instability indices of the cortical samples (green) overlap with those of the trophoblast samples (yellow); indices from carrier parents’ blood are in orange (ages 25 to 34 years). Two of the fetuses had the same (CAG)n of 46 and thus overlap on the graph. Left panel: Comparison of brain tissue instability from HD carrier fetuses at 13 weeks (green), to the corresponding trophoblasts sampled for prenatal testing (yellow) and the premanifest carrier parents’ blood (orange). Right panel: The percentage of mutant alleles bearing the different somatic expansions ascertained from the peak heights. Four graphs were plotted for the reference CAGs (40, 42, 45, and 46) determined on the fetal tissues. The parental blood samples show significant somatic expansions, whereas the trophoblast and the developing cortex show very little.

Figure 1—figure supplement 1
PCR profile and analysis of the CAG repeat length and instability.

(A) PCR profile of a post-mortem brain (from a 74-year-old Huntington disease [HD] individual) with a non-expanded allele of 17 CAG and an expanded allele of 40 CAG. The non-expanded allele consists of one main peak (with a bit of PCR stutter to the left), whereas the expanded allele shows several small peaks to the right of the main peak. The highest peak determines the reference CAG size; the additional peaks correspond to larger expansions and reflect CAG somatic mosaicism. The y axis registers the level of fluorescence for each fragment. (B) Method for calculating the expansion index (EI) and the percentage of mutant alleles. In this study, we considered peaks to represent somatic expansions only if they reached a threshold of at least 3% of the height of the main peak. (C) PCR profile of a HD post-mortem brain from a 8-year-old child carrying an expanded allele of 128 CAG and a non-expanded allele of 30 CAG. Top: Complete PCR profile with the normal allele peak. Bottom: Zoom on the expanded allele.

Somatic mosaicism increases in blood and post-mortem cortex over time.

Comparison of mosaicism in cortical tissue from Huntington disease (HD) carrier fetuses at 13 weeks (green), blood samples over time (gold, orange, and red for t1, t2, and t3, respectively) and adult post-mortem cortices (blue). We ascertained the ‘% mutant alleles’ (as in Figure 1C) from the peak heights from PCR profiles obtained on GeneMapper. Three reference CAG lengths (41, 43, or 45 CAG) were chosen from our cohort to illustrate the evolution of instability, and each graph represents one individual patient (from top to bottom): for the fetal samples, patients 2, 4, and 5 (Table 1); for the blood samples, the repeated measures from patients 4, 19, and 34 (Table 2); for the post-mortem samples, patients 2, 5, 10 (Table 3).

Figure 3 with 1 supplement
Somatic expansion correlates with reference CAG repeat length, with age at onset (AO) and age at death (AD).

Scatter plots showing (A) reference CAG repeat length, (B) expansion index (EI) estimated at age at onset (EI-AO) and (C) expansion rate (ER) with respect to AO (left panels) and AD (right panels). p-Values of the slopes and adjusted R-squared for the linear regressions of log(AO) and log(AD) are in the upper right corner of the plots. Red curves denote locally estimated scatterplot smoothing (LOESS) of the data with 95% confidence intervals shaded in gray. (A) Left panel: The reference CAG repeat length explains roughly half the variability in AO (R2adj = 47.6%, p=1.8e-8). HD individuals were classified as having onset earlier than expected (green), as expected (grey), or later than expected (red) according to their distance from the linear predictions given by the CAG repeat (see 'Materials and methods'). Right panel: The reference CAG repeat length explains 68% of the variability in AD (p=1.7e-04). (B) EI-AO explains 20.7% of the variability in AO (p=6.0e-04, left panel) and 49.7% of the variability in AD (p=2.9e-03, right panel). (C) ER explains 33% of the variability in AO (p=9.3e-06, left panel) and does not explain the variability in AD (R2adj < 0, p=0.348, right panel).

Figure 3—figure supplement 1
Pearson’s correlation analysis.

Correlations between age at onset (AO), age at death (AD) and reference CAG repeat length (CAG) with expansion rate (ER) and expansion index (EI) estimated at age at onset (EI-AO). On the lower left, plots show the relation between the different variables, with a regression line for each scatter plot. On the top right are Pearson’s correlation coefficients with p-values.

Figure 4 with 1 supplement
Residual expansion index (EI) and expansion rate (ER) correlate with disease status.

(A) Boxplots showing the distribution of the EI estimated at age at onset and ER, between the patients’ groups classified in Figure 3A as having earlier-than-expected, expected, or later-than-expected onset. p-Values for the Kruskal-Wallis test are at the top of each plot; thick horizontal lines indicate the median, diamonds indicate the mean. Earlier and later onset patients had similar ERs compared to the as-expected group (p=0.810 and 0.181). (B) Boxplots showing the distribution of the residual EI at the first visit and residual ER, between the patients’ groups classified as premanifest and manifest. p-Values for the Wilcoxon test are at the top of each plot.

Figure 4—figure supplement 1
Comparison of the relation between age at death, CAG repeat length, and expansion index (EI) as measured in blood samples from deceased longitudinal patients and in post-mortem HD brains.

(A) For the 14 patients from the longitudinal cohort who died during the course of the study, we calculated age at death (AD) relative to CAG repeat length (left) and expansion rate (ER) (right). Data shown with a LOESS curve in red and 95% confidence interval shaded in gray. p-Value of the slope and adjusted R-squared for the linear regressions of log(AD) on CAG repeat length, and residual AD on ER are provided inside the plots. Colors represent different CAG repeat lengths, as indicated by the legend. (B) For the 14 post-mortem brain sample subjects we compared reference (CAG)n to AD, and EI as measured in the frontal cortex to residual age at death. Data are shown with a LOESS curve in red and 95% confidence interval shaded in gray. p-Value of the slope and adjusted R-squared for the linear regressions of log(AD) on CAG repeat length and residual AD on EI are inside the plots. Colors represent different CAG repeat lengths, as indicated.

Evolution of somatic expansions in HD patient blood is a function of age, CAG repeat length, and the interaction between age and repeat length.

A linear mixed model was fitted to the longitudinal data using all blood samples collected. The fitted lines show that the predicted somatic expansion increases with age for all reference repeat sizes, most notably at greater reference repeat lengths. The sex of the patients is indicated by solid and dashed lines (male and female respectively). Each curve covers the same age interval observed in our cohort for a given repeat length.

Author response image 1

Tables

Table 1
Descriptive data on fetal tissues.
Fetal brain (GW13)Trophoblast (GW11)Parent blood
CAGExpansion indexCAGExpansion indexCAGExpansion index
1400.0441400.0407400.1255
2410.0494410.0495400.1073
3420.0434420.0473420.2383
4430.0489430.0482450.3434
5450.0533450.0516450.3692
6460.0523460.0571450.3692
7460.0599460.0589450.2390
Summary43.3 ± 2.40.0502 ± 0.00643.3 ± 2.40.0505 ± 0.00642.8 ± 2.50.2560 ± 0.1103
  1. GW13: gestational week 13. GW11: gestational week 11.

    Pearson's correlation between EI and (CAG)n for fetal brain: R = 0.8244, R2 = 0.6796, p=0.023.

Table 2
Data on the 50 patients in the longitudinal study, arranged according to reference CAG repeat length.


First sample


Second sample


Third sample
CAGMotoronsetAgeStatusUHDRSEIAgeStatusUHDRSEIAgeStatusUHDRSEIEREI-AOAO GroupAD
1394939P00.118750M70.13690.001660.1355Earlier
2405840P00.159158M170.19100.001770.1908Later
3408076M00.204887M360.24990.004090.2210Later91
4414446M00.212065M340.397069M750.42650.009460.1941Earlier77
5413728M00.292942M270.380645M420.42740.00740.2905Earlier
6423932P00.168846M280.217353M180.34460.007670.2079Earlier
7424232P00.196839P00.291751M250.35960.008260.2939Earlier
8424334P00.199147M310.33870.010740.2959As
expected
9424539P00.214552M220.23800.001810.2255As
expected
10424850M210.298259M760.359661M760.32640.003860.2943As
expected
66
11425051M180.194359M390.26060.008290.1859As
expected
12425050P40.407061M400.42390.001540.4072As
expected
13426062M450.377473M780,42490.004320.3685Later74
14425756M00.398170M700.49050.00660.3980Later
15425660M200.431168M450.45810.003380.4176Later75
16426163M360.456965M530.4563−0.00030.4576Later
17433729P00.180653M310.36410.007650.2418Earlier
18434543P00.315959M400.42750.006980.3299As
expected
19434537P00.287952M380.462854M430.50880.012490.3865As
expected
20434746P00.325556M320.42280.009730.3351As
expected
64
21434750M180.443058M570.51640.009170.4153As
expected
22444821P00.173129P00.206347M180.28660.004380.2862Later
23445639P00.200262M380.41360.009280.3581Later
24444432P00.217635P00.22140.001270.2327As
expected
56
25444425P00.216242M120.53790.018920.5377As
expected
26443731P00.373944M320.54770.013370.4541Earlier49
27444437P00.363550P00.46870.008090.4201As
expected
57
28443728P00.375143M350.51750.009490.4604Earlier
29445042P00.422954M380.51620.007780.4851Later
30442626M110.427936M220.51130.008340.4279Earlier
31444649MNA0.455060M860.58540.011850.4193As
expected
61
32443036P00.427250M490.56670.009970.3674Earlier54
33453720P00.355138M70.54070.010310.5303As
expected
34454035P00.452249M520.679451M550.70500.015960.5325As
expected
35454046MNA0.740858MNA1.04570.025410.5884As
expected
36464027P00.294848M430.61260.015140.4917As
expected
37463626P00.512239M580.70440.014780.6599As
expected
38463535P00.603145M910.76750.016450.6031As
expected
49
39464547M180.770155M590.94310.021630.7270Later
40473646M451.568757M922.67500.100570.5629As
expected
64
41483021P00.780038M571.63410.050241.2321Earlier39
42483233M80.963043M841.216253M921.80550.042130.8650As
expected
43493929P01.270942M841.7097421.72080.034191.6129Later
44492720P00.843927M631.35350.07281.3534Earlier
45493341M431.651243M531.978851M702.68550.099120.9159As
expected
46502531MNA2.727841MNA3.61780.0892.1936Earlier
47522520P00.847933M321.706039M692.00140.061561.1696Earlier
48522736MNA3.154236M493.3273NANAAs
expected
49532626M151.981934M223.07560.136721.982As
expected
50543620P01.944936M133.87450.12063.8745Later
44.6 ± 3.542.2±
10.8
37.8±
12.6
19 M/
31P
6.6±
12.9
0.620±
0.655
49.7 ±
12.3
46 M/
4P
41.1±
25.1
0.881±
0.929
51.3±
8.1
12M/0P53±
25.6
0.967 ± 0.8410.0236±
0.0332
0.625±
0.649
62.9±
13.6
  1. UHDRS: United Huntington's Disease Rating Scale/124; P: premanifest; M: manifest; ER : expansion rate; EI: expansion index; AD: age at death; AO: age at onset.

Table 3
Descriptive data on post-mortem brain donors.
SexCAGExpansion indexAge at death
1F401.450674
2F411.319784
3F423.244475
4M431.683455
5M432.600768
6F441.2883NA
7M443.255141
8F443.044059
9F443.708355
10F452.972259
11F461.519143
12F479.094243
13M482.850641
14F503.796158
15F508.591245
Summary11F/5M44.7 ± 3.03.3612 ± 2.390057.4 ± 13.9
Table 4
Descriptive data on the Linear mixed model analyzing the longitudinal data.

Linear mixed model fit by REML. t-tests use Satterthwaite's method ['lmerModLmerTest']. Formula: log(ei)~age + cag+age:cag +sex + (1 | id).

OutcomeFixed effectsRandom effects
EstimateStd. Errordft valuePr(>|t|)SignifGroupsNameVarianceStd.Dev.
log(ei)(Intercept)−0.6030.04848.31−12.51.03e-16***id(Intercept)0.0670.258
age0.0280.00170.3922.83.40–34***Residual0.0090.095
cag0.2760.01254.2124.01.63e-30***Number of obs: 112, groups: id, 50
sexMale0.0280.07847.030.360.719ns
age:cag0.0023.8e-477.225.692.21e-07***
  1. Signif: ***p<0.001, nsp >0.05.

Table 5
Descriptive data on longitudinal cohort.
t1t2t3
n505012
Sex, F/M (%F)31/19 (62%)31/19 (62%)7/5 (58.3%)
CAG44.6 ± 3.5
(n = 50, r = 39–54)
NANA
Expansion index0.620 ± 0.655
(n = 50, r = 0.119–3.154)
0.881 ± 0.929
(n = 50, r = 0.137–3.874)
0.967 ± 0.841
(n = 12, r = 0.287–2.686)
Age at sampling37.8 ± 12.6
(n = 50, r = 20–76)
49.7 ± 12.3
(n = 50, r = 27–87)
51.3 ± 8.1
(n = 12, r = 39–69)
Status M/P (%M)19/31 (38%)46/4 (92%)12/0 (100%)
Chorea onset yes/no (%yes)20/30 (40%)46/4 (92%)11/1 (91.7%)
UHDRS6.6 ± 12.9
(n = 46, r = 0–45)
41.1 ± 25.1
(n = 48, r = 0–92)
53 ± 25.6
(n = 11, r = 18–92)
t2 - t1t3 - t2
Delta age12 ± 4.9 (n = 50, r = 0–24)6.2 ± 5.2 (n = 12, r = 0–18)
Delta expansion0.16 ± 0.22
(n = 50, r = –0.01–1.37)
0.11 ± 0.12
(n = 12, r = –0.03–0.37)
Rate (delta expansion/delta age)0.0136 ± 0.016
(n = 49, r = –0.0036–0.0859)
0.0182 ± 0.0182
(n = 11, r = –0.0153–0.0461)
Expansion rate (ER)0.0236 ± 0.0332 (n = 49, r = −0.0003–0.1367)
Expansion index at onset (EI-AO)0.625 ± 0.649 (n = 49, r = 0.136–3.874)
Key resources table
Reagent type (species) or resourceDesignationSource or referenceIdentifiersAdditional
information
gene (human)HttNCBING_009378.1
sequence-based reagentHD-F2This paperPCR primersGGGAGACCGCCATGGCGACCCTGGA
sequence-based reagentHD-WR2-hexhttps://doi.org/10.1006/mcpr.1993.1034PCR primersHEXGGCGGTGGCGGCTGTTGCTGCTGCT
sequence-based reagentHD-WCAAM4-R-famThis paperPCR primers[6FAM]GGCGGTGGCGGCTGTTGCTGTTGAT
commercial assay or kitQIAamp Fast DNA Tissue KitQiagen51404
commercial assay or kitMaxwell RSC Blood DNA kitPromegaAS1400
commercial assay or kitTaq DNA PolymeraseQiagen201205
chemical compound, drugGenescan-400HD Rox dye size standardApplied Biosystems15829736
chemical compound, drugHi-Di FormamideApplied Biosystems15803570
software, algorithmGeneMapper software v5.0Applied BiosystemsRRID:SCR_014290
software, algorithmR version 3.6.1R Development Core TeamRRID:SCR_001905https://www.R-project.org/

Additional files

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Radhia Kacher
  2. François-Xavier Lejeune
  3. Sandrine Noël
  4. Cécile Cazeneuve
  5. Alexis Brice
  6. Sandrine Humbert
  7. Alexandra Durr
(2021)
Propensity for somatic expansion increases over the course of life in Huntington disease
eLife 10:e64674.
https://doi.org/10.7554/eLife.64674