Conceived and designed the experiments: LM ADR EL TR AK. Performed the experiments: LM. Analyzed the data: LM ADR TR AK. Contributed reagents/materials/analysis tools: LM ADR EL. Wrote the paper: LM ADR EL TR AK.
Huntington's disease (HD) is one of several neurodegenerative disorders caused by expansion of CAG repeats in a coding gene. Somatic CAG expansion rates in HD vary between organs, and the greatest instability is observed in the brain, correlating with neuropathology. The fundamental mechanisms of somatic CAG repeat instability are poorly understood, but locally formed secondary DNA structures generated during replication and/or repair are believed to underlie triplet repeat expansion. Recent studies in HD mice have demonstrated that mismatch repair (MMR) and base excision repair (BER) proteins are expansion inducing components in brain tissues. This study was designed to simultaneously investigate the rates and modes of expansion in different tissues of HD R6/1 mice in order to further understand the expansion mechanisms
Huntington's disease (HD) is a genetically determined neurodegenerative disorder identified by the presence of a mutation for a long series of CAG repeats (>36 repeats) in the Huntingtin (
Huntington's disease (HD) is a genetically determined neurodegenerative disorder, the onset of which is known to depend upon the length of glutamine-encoding CAG-repeat sequences lying within the Huntingtin
Trinucleotide repeat (TNR) instability varies between organs in a variety of neurodegenerative disorders which are caused by expansion of CAG repeats in a coding gene, with the greatest instability observed in the brain
TNR sequences may form slipped strands during replication or repair, creating loops or hairpins, which protrude from the DNA duplex
In the earliest model for repeat expansion the DNA polymerase forms slip-outs on the nascent strand leading to small-scale repeat expansion in repetitive sequences
Nevertheless, the processes causing this expansion remain poorly understood, particularly in mammalian systems, although the formation of secondary DNA structures within the repeat sequence is thought to underlie the process
In order to understand the mechanisms underlying somatic CAG instability, 42 R6/1 HD exon 1 transgenic mice were sacrificed at either 10 or 21 weeks of age, whereupon tail, heart, lung, spleen, liver, cortex and striatum samples were taken for analysis of HD CAG repeat length. A tail biopsy at 3 weeks of age represents the reference level of CAG repeats present near birth in all tissues for each mouse
Representative examples of raw data from CAG-repeat sequences in tail (A), cortex (B) striatum (C) and liver (D) from individual male and female HD mice, aged 10 and 21 weeks are shown (blue) with the tail biopsy from the same 3-week old mouse overlaid (red). All traces demonstrate the increase in mean length of repeat sequences with time and the differing rates of expansion between tissue types. Of particular note is the strong periodicity shown in the older striatum samples. Size standard markers are shown for 118 (solid black arrow) and 138 (open white arrow) CAG repeats respectively.
Individual fragments from tail fit well to a normal distribution and are thus described by the mean (μ) and standard deviation (σ) of the curves fitted to raw data (
Size standard markers, are placed at 118 (solid black arrow) and 138 (open white arrow) CAG repeats. (A) A representative example of raw data from the 3-week biopsy (red) and 21-week sample (blue) are shown for one mouse. (B) Normal distributions are fitted to the data presented in (A), coloured as previously for 3-week (red) and 21-week (blue) samples. The resulting means (μ) and standard deviations (σ) are used to define the temporal change in repeat distribution within the sample, clearly demonstrating an increase in the mean number of repeats between 3 weeks (μ3) and 21 weeks of age (μ21). Likewise, broadening of the distribution is evident from the increase in standard deviation at 3 weeks (σ3) to that at 21 weeks (σ21). (C) Mean values (μ) for repeat lengths from the 3-week (red) and 21-week (blue) tails samples of 59 mice are compiled into a histogram, showing the systematic increase in repeat length with age. (D) Standard deviations of all 59 tail samples at 3-weeks (red) and 21-weeks (blue) are similarly compiled, with the histogram showing age dependent peak broadening. (E) A boxplot of expansions measured in the tails of 59 mice between 3 and 21 weeks of age shows a median expansion of 1.97 CAG repeats.
(A) Mean versus standard deviation points are plotted for all 103 mice in the data set, with the data divided into 3-week (red), 10-week (blue) and 21-week (green) age groups. To highlight the general trends in the data, the points representing the middle quartiles (by mean value) for each age group are shown with solid circles. The 3-week age group is shown with the weak positively-correlated trend line (red) implying a loose relation between standard deviation and mean value at a fixed age, which we assume to be caused by polymerase errors – during PCR of repeat sequences – which increase with sequence length. This trait is also present in 10-week and 21-week data, with parallel approximate trends shown (dashed) for emphasis. The standard deviations and means for each age group are, however, shown to increase systematically with age, along a trend-line (black arrow) that is completely separate from the PCR-dependent trend. This demonstrates the independence of the age-dependent increase in standard deviation and mean from the PCR induced variation. (B) Monte Carlo simulation of the proposed model mechanism for expansion is shown with a range of parameters for expansion and contraction probabilities
Having made these observations, it is necessary to consider them in the context of potential models for expansion, in order to fully investigate their implications and attempt to parameterize the processes involved. A continuous increase in both mean and standard deviation can be generally accounted for by multiple stochastic unitary (single CAG-repeats) extension and contraction events on the CAG tracts within the sample (
In contrast to fragment data from tail, heart, lung and spleen, raw data from 10- and 21-week cortex and striatum tissue show a peak retained at the 3-week repeat level alongside an age-dependent number of periodically spaced subsequent peaks (
(A) Fragment analysis curves from 3-week tail, 10- and 21-week striatum are shown sequentially, with fitted normal distributions (black) overlaid. The sum of all fits (red) demonstrates the curve-fit accuracy. Mean values (μa etc) for each fitted peak are shown. Mean μa coincides with the mean of the corresponding 3-week tail at all ages. This 21-week striatum sample is best-fit by six consecutive normal distributions. Periodicity is clear from the regularity of the intervals between consecutive means. (B) Multiple peaks fitted to another 21-week striatum dataset are shown, with the areas under each fitted distribution (Aa .. Ad), the mean values and the separations between the four peaks (S1 .. S3). The area of a peak represents the proportion of tissue containing the measured mean number of repeats, thus separation values represent a step-wise expansion from the previous mean. The age-dependent propagation of peaks with higher means, as seen in (A), is due to the stochastic insertion of short repeat sequences, of consistent length, into the CAG tracts of individual cells, which, over time, generates the observed periodicity. For analysis purposes, the sum of all areas (Aa,b,c,d) was rescaled to 10,000 (the approximate number of cells forming a sample), allowing the number and length of insertion events to be estimated. (C) A histogram of insertion lengths for all expansion events measured in 69 separate striatum samples is shown (blue). Both mean and median values of the distribution point to a dominant insertion length of seven CAG repeats. The separate contributions from 10-week (green) and 21-week (red) data are also shown. Inset figure shows a similar result for cortex, with the insertion length distributed at 7 repeats, despite the smaller number of insertion events observed in total.
Having shown data and analysis to define these two distinct modes of expansion, we place our findings into the context of existing literature, in order to develop reasonable hypothetical models for these two types of expansion.
The continuous expansion we observe shows a progressive increase in CAG tract length, which is comparable with the expansions observed in fibroblasts derived from an adult HD mouse
A hypothesis for periodic expansion is also presented (
Theory suggests that CAG flaps ranging from 4 to 16 triplet repeats in length can form thermodynamically unstable hairpin structures with an even number of repeats
It should be pointed out that it is possible that other stabilized structures, such as loop-outs or a stabilizing interaction with one of the many proteins and complexes that are in contact with the DNA may also serve to cause the observed periodicity. The majority of the potentially applicable models are covered in the literature
Cell proliferation in neurons and glial cells has been observed in the subependymal layer adjacent to the caudate nucleus in human HD brains
The liver may be particularly interesting in this regard, since this tissue potentially exhibits a mix of both types of expansion mechanisms. This could be attributed to different modes of expansion occurring in different cell types. Indeed, instability of the DM1 CTG•CAG repeat is known to occur in liver hepatocytes with high DNA ploidy
The question also arises as to why FEN1 did not influence CAG repeat expansion in the organs tested. One might expect a difference, since a recent
It is worth considering briefly why this periodicity has not been described before, since there are numerous potential reasons. One possibility could be that the mice used in the present study exhibit more instability due to environmental factors
So far, we have only studied periodicity in the R6/1 mouse model and without specific studies of other HD CAG mice models the generality of our data is unknown. Yet, the R6/1 transgenic mouse is a widely accepted and commonly used model for human HD that exhibits a progressive neurological phenotype that exhibits many of the features of HD
Importantly, CAG repeat expansion in human cortex is associated with an earlier age of disease onset, in addition to the role of the constitutional CAG repeat length
To summarize, we present two different mechanisms of somatic CAG repeat expansion; a continuous bi-directional expansion in tail, lung, heart and spleen tissues, and a dramatic periodic expansion centred around 7-repeat insertions in striatum and cortex. Further experiments are needed to determine whether the 7-repeat step-size is independent of CAG tract length and species and it remains to be shown whether these two models can explain the expansions observed in other brain tissues and organs, as well as in humans. Nevertheless, these results provide significant new insights into
All animal experiments have been approved by the local and national animal - and are carried out according to the regulation by FELASA (Federation of European Laboratory Animal Science Associations).
B6CBA-Tg(HDexon1)61pb/J mice of the R6/1 line
DNA from tail biopsies of 3-week old mice were used for genotyping. HD mice were genotyped with forward
CAG repeats were sized by PCR with primers 5′-FAM-
All raw data was processed through a masked Nelder-Mead simplex fitting method, optimising free parameters of standard deviation, mean and amplitude to fit consecutive normal distributions sufficient to account for ≥98% of the total area of the raw data set. In the case of tail data, only a single normal distribution was required. These optimised parameters were returned as the means (μ), and standard deviations (σ), which were used to define the TNR lengths present in each data set. CAG repeat tracts were flanked by sequences 86 bp in length as verified by sequencing. Thus, the mean number of CAG triplets (μt) present in a fragment analysis sample with a measured mean (μm) is defined by μt = (μm−86)/3. When analysing the periodicity present in striatum and cortex data, standard frequency analysis methods are not suitable, therefore our peak fitting method was used to fit consecutive normal distributions to the raw data (
Further examples of curve-fits to raw data from striatum. We present a further set of raw fragment analysis curves paired immediately below with the corresponding normal distribution curve-fits to samples from 21-week striatum samples, in order to confirm the prevalence of the periodicity seen in the data. Individual normal distribution fits are shown in magenta, with the sum of all fitted curves shown in cyan.
(1.49 MB PDF)
Click here for additional data file.
Further examples of curve-fits to raw data from cortex. We present a further set of raw fragment analysis curves paired immediately below with the corresponding normal distribution curve-fits to samples from 21-week cortex samples, in order to confirm the prevalence of the periodicity seen in the data. Individual normal distribution fits are shown in magenta, with the sum of all fitted curves shown in cyan.
(0.49 MB PDF)
Click here for additional data file.
Comparison of 6-week striatum with 3-week tail. In order to further confirm that expansion peaks in striatum tissue develop with time and that the 3-week tail sample is a legitimate representation of the condition of a given 3-week organ sample, we present the results from a single mouse which was sacrificed at 6 weeks of age. Here we show the results from the 3-week tail sample (A), to be compared with the striatum sample at 6-weeks with a mean of 122 repeats (B) and cortex with a main peak at 120 repeats and a small peak at 126 repeats (C). The striatum sample shows a mean slightly higher than the 3-week tail mean. However, the limitations of the analysis technique mean that the distribution could be similar to that shown for cortex, with two peaks which are bundled together by the curve-fitting algorithm. To maintain the unbiased nature of the analysis, it is necessary to admit the uncertainly and refrain from attempting guided curve-fitting. Irrespective of this, the small disparity between the 6-week striatum and 3-week tail mean values supports our stated assumption that the near-birth repeat level in all organs is well approximated by the level measured in a 3-week tail biopsy.
(0.24 MB PDF)
Click here for additional data file.
Proportions present in samples with mixed populations. Our ability to estimate the proportions present in samples with mixed populations of repeat lengths was tested, as shown below. Two separate samples containing 120 and 129 repeats were mixed in a variety of ratios, and then processed by fragment analysis. The spectra were analyzed as described in the paper and the fitted curves are shown below. The results show that the technique is more than adequate for the estimation of sample distributions, particularly when applied to an entire data set.
(0.27 MB PDF)
Click here for additional data file.
Striatum samples separated by repeat length after three consecutive 10× dilutions. To clarify that individual fragment analysis samples from striatum are composed of DNA, we decided to perform a limited set of serial dilutions on a single 21-week striatum sample and then proceed with fragment analysis, in order to see whether the original sample could be divided up into individual peaks, thereby further validating our multiple curve-fitting analysis method. A) Fragments from an original sample showing separable peaks (black below) were subjected to 3 successive 10× dilutions. This resulted in a group of samples in which there was a finite chance that fragments of only one, or a few, different lengths would exist, compared to the original sample. Several of these samples were then amplified by PCR and the peaks from three separate samples are shown here (blue, green and red), with clearly discernable means that correlate with the peaks shown in the original data. These results show that where individual length fragments are separated out by dilution, their means align well with the peaks that are visible in the original data, reinforcing the conclusion that striatum samples contain TNR tracts which have expanded periodically by multiples of ∼7 insertions. While a small pool PCR (spPCR) technique has been used to show that the fragment analysis curve resembles the distribution of individual fragment lengths measured in small pools (Gonitel et al. DNA instability in postmitotic neurons. ProcNatlAcadSci USA (2008) vol. 105 (9) pp. 3467–72
(1.14 MB PDF)
Click here for additional data file.
Simulation of periodic expansion with and without slipped-strand expansion. Here we present the results of three simulations of 7-repeat periodic expansion in a population of 20,000 cells, with varying levels of simultaneous slipped-strand expansion and contraction as observed in tail. In all cases the probability of a 7-repeat step, a unitary expansion and a unitary contraction are shown. In the first case (A), no slipped strand expansion and contraction are allowed, leaving clearly defined peaks throughout the data, similar to those observed in real striatum data. (B) is simulated with approximately a quarter of the level of slipped strand expansion and contraction, which renders the peaks indistinguishable. In (C) with the measured levels of expansion and contraction in tail, the distribution becomes utterly uninformative. This is the basis for our argument that slipped-strand and periodic expansions do not occur simultaneously in striatum.
(0.11 MB PDF)
Click here for additional data file.
Increased replication does not increase expansion rate in spleen. Proliferating cell nuclear antigen staining of frozen spleen sections from HD mice with Fen1 mutation (lymphoma) and without Fen1 mutation (normal spleen) are shown (A). Histograms of expansion levels in 21-week spleen samples presenting lymphomas from HD mice with Fen1 mutation (blue) and normal spleens from HD mice (red), show no significant difference (B). Increased replication in lymphoma tissues does not affect the rate of CAG repeat expansion. The continuous expansion in spleen therefore appears to be independent of replication.
(1.78 MB PDF)
Click here for additional data file.
Increased replication does not increase expansion rate in spleen. (A) A proposed model mechanism for the slight expansions measured in CAG tracts of tail tissue is shown. Since this expansion is the result of a multitude of small expansion and contraction events within individual cells, we present a model that accounts for both processes by alternate processing of migrated loop-outs formed as slipped-strand structures within the repeat sequence. Initiated by a bubble, loop-out, or cruciform structure, the loops on opposite strands may migrate apart, rather than resolving back into a duplex formation. Should a loop be processed as an error by either of the alternate mechanisms shown, a single contraction or expansion event can occur. Removal of a loop structure on either strand causes a single repeat contraction, while nicking and gap-filling on the strand opposite to a loop results in a single repeat expansion. A slight bias favouring expansion over contraction will result in the overall population expansion measured in tissue samples. (B) Explanation of simulations videos: All videos show the starting 3-week distribution (blue) and the daily progress of the distribution under the given expansion mechanism, up to 21-weeks (red). Additionally, for extra clarity, inset on the right is a matrix of 10,000 points coloured by the number of repeats in each individual cell, to represent the actual number of repeats in each cell. Meanwhile, inset on the left is a histogram of the probabilities for expansion and contraction for each CAG repeat insert length, which have been chosen for the simulation.
(2.11 MB PDF)
Click here for additional data file.
Hypothetical hairpin-based model for periodic expansion. There are several potential models that could explain a periodic expansion mechanism with an average expansion length of 7 repeats. All of these require some inherent stabilisation of a loop of DNA, whether it be self-hybridizing, or in coordination with a protein or complex which interacts with DNA. In the following model we propose a means by which a hairpin can be stabilized around an average of 7 repeats in length. The hairpin formation should be considered as a Markov chain of transitions between states, whereby the mean/most-stable state is gives a 7-repeat hairpin and other lengths are distributed around this mean. A proposed pathway for periodic expansion by the step-wise insertion of 7-repeat hairpins is shown. After a preliminary strand break the CAG triplets are displaced from the DNA duplex and initially form an unordered flap. As the flap length increases, the flap alternates between metastable folded and disordered states. Continued strand displacement increases the overall length of the flap up to and beyond, 6 repeats. At this point, depending upon the instantaneous state of flap folding, there are effectively two possible folding pathways. The first pathway is shown to the left (blue background), whereby the first six repeats of the flap fold stably into a hairpin, which is eventually recognised as an erroneous structure and correctly repaired. There is no resulting expansion. A finitely probable alternative pathway is shown to the right. A stable 6-repeat hairpin is formed from CAG triplets 2 to 7 of the flap, leaving a single CAG repeat on the 5-end of the hairpin. The overhanging CAG triplet may hybridize with an unpaired CTG on the complementary strand, temporarily stabilizing the hairpin on the duplex. The stabilizing flap facilitates gap-filling repair and ligation of the CAG loop to the duplex DNA. Subsequent repair of a nick or a lesion along the CTG strand causes the extra CAG repeats to be spliced into the sequence, generating a single 7-repeat expansion event. Inevitably there will be more or less stable - and consequently likely - flap conformations with lengths greater or less than 7 repeats, but a mean-length or most probable length of 7-repeats would produce the observed expansions. While this may not be the only means of inserting a longer stabilized repeat sequence into the original sequence, the essential length-dependent stability of loopout structures which may also function as the basis for expansion would be related to the stability of the hairpins illustrated here.
(0.54 MB PDF)
Click here for additional data file.
Replicate PCR, striatum samples from the same mouse. In order to further confirm the periodicity we observe in our striatum samples, we have performed replicate PCR on 9 samples from the same mouse, which we show below. As is clear from the figures, periodicity is a consistent feature of the samples, although the size and mean position of the peaks show some variation between replicate samples. The variability between peak heights could be attributable to sampling error during the PCR preparation or one-repeat expansion or contraction early in the PCR cycles as described in
(1.35 MB PDF)
Click here for additional data file.
Examples of data from Heart, Spleen and Lung, matching the pattern shown in tail data. In order to confirm our statement in the manuscript that Heart (A), Spleen (B) and Lung (C) data all show similar monomodal distributions to Tail data, we have included examples of curves from these three tissue types.
(0.81 MB PDF)
Click here for additional data file.
PCR of repeat sequences and the observation of periodicity. (A) Simulated error-prone PCR from a single molecule. Starting with a single 119 repeat template, and plotting all PCR products (coloured by length) on equal sized areas, we can see the development of length variability in the population shown. See
(1.97 MB PDF)
Click here for additional data file.
A deeper discussion of model parameters and estimations.
(0.03 MB DOC)
Click here for additional data file.
PCR of repeat sequences and the observation of periodicity.
(0.03 MB DOC)
Click here for additional data file.
CAG expansion/contraction 1.
(1.42 MB MOV)
Click here for additional data file.
CAG expansion/contraction 1 2 3.
(1.42 MB MOV)
Click here for additional data file.
CAG expansion/contraction 3.
(1.38 MB MOV)
Click here for additional data file.
CAG expansion 4 5 6 7 8 9 10.
(1.44 MB MOV)
Click here for additional data file.
CAG expansion/contraction 5.
(1.39 MB MOV)
Click here for additional data file.
CAG expansion/contraction 5 6 7 8 9.
(1.40 MB MOV)
Click here for additional data file.
CAG expansion/contraction equal probability 5 6 7 8 9 10.
(1.62 MB MOV)
Click here for additional data file.
CAG expansion/contraction 7.
(1.31 MB MOV)
Click here for additional data file.
The authors would like to thank Prof. Eric Rinvik at the University of Oslo for dissection of striatum and cortex. We also acknowledge Dr. Trine J. Meza for taking part in the initial stage of this project and Terese Kristoffersen, Ingvild Høkholm, Guro F. Lien, Anja Bjølgerud, Linda Ellevog, Ophélie Aussedat, Hege Wicksén, and Gaute Nesse for technical assistance.
The authors have declared that no competing interests exist.
This work was supported by the Norwegian Cancer Society and by the Research Council of Norway. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.