Skip Navigation

This Article
Right arrow Abstract Freely available
Right arrow Print PDF (140K) Freely available
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Right arrow Search for citing articles in:
ISI Web of Science (77)
Right arrowRequest Permissions
Right arrow Commercial Re-use Guidelines
for Open Access NAR Content
Google Scholar
Right arrow Articles by Petruska, J
Right arrow Articles by Goodman, M.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Petruska, J
Right arrow Articles by Goodman, M.
Social Bookmarking
 Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?

© 1996 Oxford University Press 1992-1999

Footnote

Stability of intrastrand hairpin structures formed by the CAG/CTG class of DNA triplet repeats associated with neurological diseases

Stability of intrastrand hairpin structures formed by the CAG/CTG class of DNA triplet repeats associated with neurological diseases John Petruska , Norman Arnheim and Myron F. Goodman*

Department of Biological Sciences, Hedco Molecular Biology Laboratories, Molecular Biology Program, University of Southern California, Los Angeles , CA 90089-1340, USA

Received March 26, 1996; Revised and Accepted April 19, 1996

ABSTRACT

Expansions of trinucleotide repeats in DNA, a novel source of mutations associated with human disease, may arise by DNA replication slippage initiated by hairpin folding of primer or template strands containing such repeats. To evaluate the stability of single-strand folding by repeating triplets of DNA bases, thermal melting profiles of (CAG) 10 , (CTG) 10 , (GAC) 10 and (GTC) 10 strands are determined at low and physiological salt concentrations, and measurements of melting temperature and enthalpy change are made in each case. Comparisons are made to strands with three times as many repeats, (CAG) 30 and (CTG) 30 . Evidence is presented for stable intrastrand folding by the CAG/CTG class of triplet repeats. Relative to the GAC/GTC class not associated with disease, the order of folding stability is found to be CTG > GAC [approx] CAG > GTC for 10 repeats. Surprisingly, the folds formed by 30 repeats of CTG or CAG have no higher melting temperature and are only 40% more stable in free energy than those formed by 10 repeats. This finding suggests that triplet expansions with higher repeat number may result from the formation of more folded structures with similar stability rather than fewer but longer folds of greater stability.

INTRODUCTION

At least nine human diseases, most of which are neurodegenerative, are known to result from expansions of trinucleotide repeats in genomic DNA ( 1 - 5 ). Of all possible classes of repeating triplets in DNA, three have been associated with disease; namely, the CAG/CTG class, the CGG/CCG class and, most recently, the GAA/TTC class ( 6 ). Expansions of the CAG/CTG class of triplet repeats, including permutations AGC/GCT and GCA/TGC, are associated with Huntington's (HD) and six other neurological diseases; those of the CGG/CCG class, with Fragile X syndromes, FRAXA and FRAXE.

Genetic instability in the number of tandemly repeating triplets is a major characteristic of these disorders (reviewed in 1 - 5 ). The number of triplet repeats in disease-associated genes tends to increase in germline transmission from parent to offspring. Usually, the larger the repeat number the higher the probability of expansion and the earlier and more severe the disease.

In the present paper, we examine the single-strand folding properties of the CAG/CTG of triplet repeats associated with neurological diseases. For comparison, we also examine structures formed by repeating triplets of the GAC/GTC class, not associated with any known disease. A favored molecular explanation for genetic expansions of tandemly repeating triplets is primer/template slippage in DNA replication, an idea first proposed to explain frameshift mutations in simple sequence repeats ( 7 ). The tendency for primer to slip on template during replication of triplet repeat regions is attributed to the ability of primer or template strands to fold into hairpin structures stabilized by intrastrand basepairing (8-13). By characterising heat denaturation profiles of disease-associated (CAG/CTG) and unassociated (GAC/GTC) triplet repeat sequences, we are able to evaluate thermal stabilities of triplet repeat folding that may be important for understanding the slippage process at the molecular level.

MATERIALS AND METHODS

Strand preparations

Synthetic 30mer DNA strands with 10 triplet repeats-(CAG) 10 , (CTG) 10 , (GAC) 10 and (GTC) 10 -were made using an Applied Biosystems 392 DNA/RNA synthesizer with [beta]-cyanoethyl phosphoramidites. Each strand was purified by electrophoresis on 12% polyacrylamide gel in 8 M urea, cutting out the desired band and extracting DNA from the band. Extracted samples at high strand concentration (absorbance A 260 = 5-10/cm) were dialyzed extensively against low ionic strength buffer (5 mM NaH 2 PO 4 , 5 mM Na 2 HPO 4 , 1 mM Na 4 EDTA, pH 7.0) and stored frozen at -70oC.

Synthetic 90mer DNA strands-(CTG) 30 , (CAG) 30 and random (C,A,G) 30 -were purchased from Operon Technologies, Inc. The purchased samples, at high strand concentration (absorbance A 260 [approx] 10/cm), were dialyzed against low ionic strength buffer and stored frozen as above.

The DNA strands used in our melting studies each had an OH group on the 5'-end as well as on the 3'-end.

Melting studies

Thermal denaturation of DNA samples is first examined in the low ionic strength buffer (19 mM Na + ) used to store frozen samples. At this low counterion concentration, melting of secondary structure occurs <80oC even for correctly paired duplexes such as (CAG) 30 @(CTG) 30 . We examine melting profiles obtained by measuring A 260 versus temperature T, while T is raised from 20 to 80oC at a constant rate, starting at 2oC/min. A slower rate of heating (0.5oC/min) is used to establish equilibrium conditions needed to evaluate standard enthalpy change ([Delta]Ho) from melting curve slope by a van't Hoff expression ( 14 ). When an ordinary cuvette is heated to 2oC/min, a teflon cover tightly wrapped with parafilm is sufficient to prevent significant evaporation; at slower rates of heating, a layer of mineral oil, >= 5 mm, on the sample is also needed to prevent excessive evaporation. To obtain corresponding curves at physiological ionic strength, 1/33 vol 5 M NaCl is added, raising [Na + ] to 167 mM. Unless stated otherwise, the total concentration of DNA bases in each case is ~0.1 mM, with A 260 [approx] 1 for 1 cm pathlength.

To reveal differences in secondary structure, we present melting curves obtained at 2oC/min, showing the A 260 ratio, A 260 (ToC)/ A 260 (80oC), with 80oC representing the melted state. The observed A 260 value at 80oC in each case is given in the figure legends. The curve of A 260 ratio versus T is presented for each strand and also for the correctly paired duplex formed by annealing equimolar amounts of the two complementary strands. To establish equilibrium, the following melting profiles were carried out at a slower heating rate (0.5oC/min): (CTG) 10 , (CAG) 10 , (CAG) 30 , (GAC) 10 , (GTC) 10 @(GAC) 10 and (CTG) 10 @(CAG) 10 in low salt (19 mM Na + ); (CTG) 10 was also examined at slow rate in high salt (167 mM Na + ). The slower heating caused essentially no change in the slopes for any of the curves >37oC obtained at 2oC/min, indicating that stably folded structures are evaluated close to equilibrium at both rates, and only these structures are included in Table 1 .

Evaluation of T m and [Delta]Ho from melting curves

A sigmoidal (S-shaped) melting profile represents a transition from ordered to disordered states. The temperature at the sigmoidal midpoint, where the curve becomes linear with steepest slope, is identified as T m , the melting temperature. At T = T m ordered and disordered structures are considered equally probable, so [Delta]G = 0 and [Delta]S = [Delta]H/T m are the corresponding free energy and entropy changes for the melting transition centered on T m . When temperature is raised slowly enough to assure equilibrium and the melting curve is properly normalized, the slope of the linear region centered on T m can be used to evaluate the standard enthalpy change [Delta]Ho in the transition, according to a `van't Hoff' formula derived by Marky and Breslauer ( 14 ). The normalized slope ([sigma]) around T m is measured after drawing linear baselines at low and high temperatures ( 14 ) and normalizing absorbance values between the baselines so that they range from 0 to 1, with 0.5 being the value at T m . The `van't Hoff' value of [Delta]Ho is obtained by the formula (14), [Delta]Ho = (2 + 2n)RT m 2 [sigma], where [sigma] is the normalized slope or first derivative (with respect to T) at T = T m , with R being the gas constant (1.987 cal/mol-K) and n the molecularity or number of strands in the ordered state (n = 1 for single-stranded state, 2 for double-stranded). To measure [sigma] as the first derivative, we use a linear least-squares fit to at least 10 data points in the linear region centered on T m.

End-labeling and electrophoretic analysis

To distinguish between single-strand structures and dimeric or multimeric associations, samples of each strand are labeled at the 5'-end with 32 P, using [[gamma]- 32 P]ATP and T 4 polynucleotide kinase, and subjected to 12% polyacrylamide gel electrophoresis (PAGE) in 0.5* TBE (45 mM Tris-borate pH 8.0, 1 mM EDTA). After frozen samples are thawed, end-labeling is carried out at low temperature (<20oC) to retain unstable associations promoted by freezing. Each labeled sample is divided into two parts, one kept on ice until loading on PAGE apparatus, the other kept on ice until ~30 min before loading, when it is preheated to 37oC for 30 min to deliberately denature unstable components. During PAGE, conducted in cold room, temperature is maintained at 4oC and current is kept low enough (5 mA at 25 V/cm) to avoid heat denaturation during electrophoresis.

RESULTS

The thermal stability of single-strand folding by triplet repeat sequences is examined at low (19 mM) and physiological (167 mM) salt concentrations, using melting curves obtained by measuring A 260 versus T at a constant heating rate. At a rate of 2oC/min, we observe in some cases low-temperaure (<37oC) sigmoidal features of unstable associations between strands promoted by freezing as well as higher temperature sigmoidal profiles of stably folded structure. At a slower rate (0.5oC/min) the low-temperature features become much less conspicuous, because sigmoidal slope is reduced as expected for non-equilibrium associations between two or more strands, but the higher-temperature profiles show very little change as expected for single strand folding close to equilibrium. Since slower heating, which requires an oil layer to retard evaporation, does not significantly alter the melting profiles >37oC, only melting curves obtained at 2oC/min (without oil) are presented.

The disease-associated triplet sequence, 5'-CAG-3', and its complement 5'-CTG-3', are each examined in two strand lengths, 30mer strands with 10 triplet repeats and 90mer strands with 30 repeats. Also, for comparison, the `sister' sequences having the same bases in reverse order, 5'-GAC-3' and 5'-GTC-3', which are not associated with any known disease, are each examined in 30mer strands with 10 repeats. The denaturation profiles of single-strand folded structures are also compared with those of their correctly paired double-helical structures formed by slowly melting and annealing equimolar mixtures of the two complementary strands.

Melting curves for CAG and CTG repeats

At low ionic strength (19 mM Na + ), where even normal DNA duplexes melt <80oC, sigmoidal melting curves are obtained (Fig. 1 a) for each of the 30mer strands, (CAG) 10 and (CTG) 10 , as well as for their correctly paired duplex, (CAG) 10 @(CTG) 10 . The corresponding curves for the 90mer counterparts, (CAG) 30, (CTG) 30 and (CAG) 30 @(CTG) 30 , are shown in Figure 1 b. Included for comparison is the approximately linear curve (Fig. 1 b, light dashed line) obtained for (C,A,G) 30, a `randomized' version of (CAG) 30 .


Figure 1 . Thermal melting profiles of secondary structure formed by CAG and CTG triplet repeats in synthetic 30mer and 90mer DNA strands at low ionic strength. Shown plotted against temperature T in each case is the A 260 ratio, defined as A 260 (T)/ A 260 (80oC) with 80oC representing the melted state. The midpoint (+) of a sigmoidal (S-shaped) region indicates the melting temperature T m of secondary structure. ( a ) Separate 30mer strands, (CAG) 10 and (CTG) 10 , and their annealed duplex, (CAG) 10 @(CTG) 10 . The observed value of A 260 (80oC) in each case is 1.32, 0.90 and 1.10 respectively. ( b ) Corresponding results for the 90mer strands, (CAG) 30 and (CTG) 30 , and their duplex, (CAG) 30 @(CTG) 30 ; with A 260 (80oC) = 1.24, 1.00 and 1.12 respectively. Also shown in (b) is a single-stranded `random' DNA, (C,A,G) 30 , similar in length and composition to (CAG) 30 but having a sequence of C, A, G bases chosen at random.

The sigmoidal character, indicative of a cooperative transition from ordered to disordered states, is more pronounced for the CTG repeats (Fig. 1 a and b, thick line) than for the complementary CAG repeats (Fig. 1 a and b, thin line). The corresponding melting temperature, T m , evaluated at the sigmoidal midpoint (+) is also higher. As seen in Figure 1 a, T m is 51oC for (CTG) 10 compared with 47oC for (CAG) 10 . The duplex (CAG) 10 @(CTG) 10 , formed by annealing the two complementary strands, exhibits a sharper transition at higher temperature, T m = 67oC (Fig. 1 a, dashed curve), as expected for cooperative melting of 30 base pairs (bp) in a normal DNA double-helix stabilized by nearest-neighbor base stacking ( 15 - 18 ).

For the (CTG) 10 strand alone, the steepness or slope of the sigmoidal region is almost a third that of the 30 bp duplex, suggesting that this strand has a secondary structure with almost a third as many stacked base pairs ( 14 , 18 ). The (CAG) 10 strand shows a considerably lower slope indicating looser or shorter structure with less favorable base stacking. For each single-strand structure, we find no change in T m over a 32-fold range of strand dilution, consistent with the melting of single-stranded intramolecular structure ( 14 ).

The (CTG) 10 melting curve (Fig. 1 a) also reveals a minor sigmoidal component with T m = 29oC in addition to the major (T m = 51oC) component. This small component is a non-equilibrium feature that appears only when concentrated strand solutions are frozen and gradually disappears as solutions are diluted (data not shown). Slower melting (at 0.5oC/min) makes this feature much less conspicuous by reducing its sigmoidal slope, but has practically no effect on the slope or T m of the major component. As verified later by electrophoretic studies, the minor component represents an unstable dimeric association between like strands, while the major component represents the stable intrastrand folding.

The longer strands, (CAG) 30 and (CTG) 30 , show a somewhat sharper melting profile (Fig. 1 b) consistent with a greater number of stacked bases ( 14 , 18 ). However, surprisingly each of these 90-base strands has no higher melting temperature than its 30-base counterpart (Fig. 1 a). The (CAG) 30 strand by itself has a T m value (46oC) 1oC lower than found for (CAG) 10 in Figure 1 a, while (CTG) 30 has the same value (51oC) as observed for the major single-strand component of (CTG) 10 . In contrast, the 90-bp duplex, (CAG) 30 @(CTG) 30 , has T m = 74oC (Fig. 1 b) or 7oC higher than the corresponding 30-bp duplex (Fig. 1 a), consistent with known T m increases with increasing duplex length ( 15 ).

Also, we see that a randomly selected sequence, (C,A,G) 30 , shows too little sigmoidal character for a T m assignment (Fig. 1 b, light dashed line). This serves as a control to confirm that the sigmoidal features of (CAG) 30 and (CAG) 10 single strands, while not as pronounced as those of (CTG) 30 and (CTG) 10 , are still much greater than exhibited by a nonrepetitive sequence of similar base compostion.

Melting curves for GAC and GTC repeats

Sigmoidal melting features are also observed for GAC and GTC repeats (Fig. 2 a), which have the same bases as CAG and CTG respectively, but in reverse order. The curves shown for (GAC) 10 and (GTC) 10 in Figure 2 a are obtained with strand concentrations similar to those in Figure 1 a, with nearly equivalent absorbance in the denatured (80oC) state (see legend to Fig. 2 a). The (GAC) 10 @(GTC) 10 duplex formed by annealing the two strands has a melting curve (Fig. 2 a, dashed line) almost identical to that observed for (CAG) 10 @(CTG) 10 in Figure 1 a. However, the single strands have very different profiles. The strand with GTC repeats shows only one sigmoidal transition (T m = 39oC), whereas that with GAC repeats shows two (34 and 47oC).


Figure 2 . Melting profiles observed for repeats of the `sister' triplets, GAC and GTC, in 30mer DNA strands at low ionic strength. ( a ) Separate (GAC) 10 and (GTC) 10 strands and (GAC) 10 @(GTC) 10 duplex; with A 260 (80oC) = 1.66, 1.34 and 1.50 respectively. ( b ) (GAC) 10 profile obtained after freezing, same as in (a), compared with profile obtained after 8-fold dilution, (GAC) 10 -dilute, or after melting and reannealing without freezing, (GAC) 10 -remelt.

The low-temperature transition of GAC repeats, like the minor one of CTG repeats (Fig. 1 a), is a non-equilibrium feature arising from strand-strand interactions promoted by freezing. This feature is observed after strand solutions are frozen and thawed (Fig. 2 b, solid line), but not if solutions are made sufficiently dilute (Fig. 2 b, dashed line) or if a remelt is done after heating and cooling without freezing (Fig. 2 b, dotted line). As in the case of CTG repeats, slower heating reduces the sigmoidal character of the low-temperature transition, but has much less effect on the slope or T m of the higher-temperature transition representing the unfolding of single-stranded structure.

Analysis by end-labeling and electrophoresis

By labeling strand 5'-ends with 32 P and using electrophoresis (12% PAGE) to separate structures by size, we confirm that the unstable structures promoted by freezing are mainly dimeric associations between like strands. However, because strands are diluted by end-labeling and by loading and subsequent electrophoresis on the gel, the amount of inter-strand association seen in the gel is likely to be less than in melting curves. Figure 3 shows PAGE results obtained with and without preheating (to 37oC for 30 min) to melt the unstable components of (CTG) 10 and (GAC) 10 having T m = 29 and 34oC respectively, seen in Figures 1 a and 2 a. The preheated samples are in lanes on the left (Fig. 3 ); the unheated samples are on the right, in the same order. By comparing each `preheated' lane on the left with the corresponding `unheated' lane on the right, and referring to the size of duplex DNA markers (on the extreme right), we can gauge the approximate size of structure removed by preheating.


Figure 3 . Electrophoretic examination of triplet repeat secondary structures. Samples of DNA used in melting studies at low ionic strength (Figs 1a and b and 2a) were 32 P-end-labeled and subjected to electrophoresis on 12% polyacrylamide gel with and without preheating to 37oC to remove unstable structures formed by freezing concentrated solution. The preheated samples are shown in the lanes on the left, the unheated samples in the corresponding set of lanes on the right. Shown in the lane between the two sets and also in the lane on the extreme right are double-stranded DNA markers of known size.

In the case of the 30-base strand (CTG) 10 , for example, a minor band about the size of a 30 bp duplex is present in the unheated (CTG) 10 lane but not in the heated lane, while the major band remains at ~15 bp in size in both lanes (Fig. 3 ). Thus, the minor component melting at 29oC (Fig. 1 a) is apparently a dimer, while the major stable component with T m = 51oC is a monomer of a single-strand structure resembling a 15 bp duplex in size. Similarly, the low-melting component of (GAC) 10 with T m = 34oC (Fig. 2 a), is seen as a dimer while the component with T m = 47oC is a monomer. For the other 30mer strands, (CAG) 10 and (GTC) 10 , showing only single components with melting temperatures of 47 and 39oC respectively, only one gel band is observed, as expected for single-strand hairpin folding in each case.

For folded 90mer strands, we observe a main electrophoretic band at ~45 bp consistent with single-strand folding (Fig. 3 ). The (CTG) 30 band migrates slightly faster than (CAG) 30, in agreement with the more rapid (CTG) 10 migration compared with (CAG) 10 , as expected since base T is smaller than A. Also a minor dimer band (~90 bp) is evident for (CTG) 30 but not for CAG 30 , consistent with our observations of dimerization by (CTG) 10 but not (CAG) 10 . Apparently, the dimer formed by 30 CTG repeats is more stable than that formed by 10 repeats because it withstands heating to 37oC (Fig. 3 ). Perhaps this interstrand association is almost as stable as the intrastrand fold so that only one sigmodial region is observed in the (CTG) 30 melting curve (Fig. 1 b).

Stabilization of folded structures by added salt

By adding NaCl to raise counterion concentration from low (19 mM) to physiological (167 mM) levels, we find that single-strand folding of triplet repeats is greatly stabilized, much like normal duplex DNA ( 15 - 18 ). For example, in the case of single strands with 10 or 30 repeats, increases in T m of 13-15oC are observed (Table 1 ), compared with 16-17oC for their correctly paired duplexes.

In Table 1 , along with T m are presented [Delta]Ho, the standard enthalpy change upon melting as deduced from melting curve slope, and the corresponding free energy change [Delta]Go evaluated at 37oC. To obtain [Delta]Go = [Delta]Ho - T[Delta]So, we use the entropy expression, [Delta]So = [Delta]Ho/T m , with T m and T in degrees K, and [Delta]Ho obtained from the normalized slope at T m . The result is [Delta]Go = [Delta]Ho(T m - T)/T m whose value at T = 37oC is shown in each case (Table 1 ). Since [Delta]Ho is positive (heat absorbed) and T m in degrees K is positive, the calculated [Delta]Go has the same sign as T m - 37oC and indicates the stability of folded structure at 37oC

Table 1 . Comparison of melting temperature T m , standard enthalpy change [Delta]H o , and corresponding free energy change [Delta]Go at 37oC required to melt secondary structure formed by DNA triplet repeats at low and physiological ionic strengths a
Repeat structure

Ionic strength, [Na + ]

0.02 M

0.17 M

T m

[Delta]Ho

[Delta]Go (37oC)

T m

[Delta]Ho

[Delta]Go (37oC)

oC

kcal/mol

kcal/mol

oC

kcal/mol

kcal/mol

Hairpin folding

(CTG) 30

51

72

3.1

66

76

6.5

(CAG) 30

46

48

1.4

60

50

3.5

(CTG) 10

51

52

2.2

66

55

4.7

(CAG) 10

47

36

1.1

60

38

2.6

(GAC) 10

47

39

1.2

62

40

3.0

(GTC) 10

39

48

0.3

53

50

2.4

Base-paired duplex

(CAG) 10 @(CTG) 10

67

190

16.7

83

220

28.3

(67)

(188)

(16.6)

(85)

(226)

(30.3)

(GAC) 10 @(GTC) 10

66

200

16.9

83

230

29.2

(74)

(200)

(21.3)

(87)

(230)

(31.9)

a The standard free energy change for melting, [Delta]Go = [Delta]Ho(T m - T)/T m , is evaluated at T = 37oC = 310 K, using T m and [Delta]Ho obtained from sigmoidal midpoint and normalized slope of melting curve as described in Materials and Methods. Also shown for each duplex are expected values (in brackets), based on nearest-neighbor doublet evaluations in normal double-helical DNA in 0.02 M and higher Na + concentrations ( 15 ).

As shown in Table 1 , the addition of 0.15 M NaCl, to raise [Na + ] from 0.02 to 0.17 M, greatly stabilizes the single-strand folding of triplet repeats, raising both T m and [Delta]Go by substantial amounts. We note that [Delta]Go, being equivalent to the product ([Delta]Ho/T m )(T m - T), tends to change in proportion to T m - T for strands of equal length, since [Delta]Ho/T m stays relatively unchanged. The folding of 10 or 30 CTG repeats has T m elevated by 15oC (51-66oC), the somewhat less stable folding of CAG repeats, by 13-14oC. Comparable elevations are also seen for 10 GAC and GTC repeats. These elevations in T m are almost as large as we find for the correctly paired duplexes, (CAG) 10 @(CTG) 10 and (GAC) 10 @(GTC) 10 , namely ~16-17oC, which agree with known salt-induced T m increases in normal DNA duplexes ( 15 - 18 ).

In contrast, the interstrand associations promoted by freezing show much less stabilization by salt. The dimeric association of (CTG) 10 shows only a 7oC increase in T m (29-36oC), while that of (GAC) 10 shows a 3oC decrease (34-31oC). Since their T m values remain <37oC, such associations are unstable under physiological conditions and therefore are not included in Table 1 .

DISCUSSION

Single-strand folding by the CAG/CTG class of triplet repeats has now been examined with chemical probes, NMR, electrophoresis, thermal melting analysis and theoretical calculations ( 10 - 13 , 19 , 20 ). All these studies, including our own, are consistent with the idea that such repeats form stable hairpin folds that may be obstacles causing slippage in DNA replication. Our results suggest, however, that thinking of these structures as `classical' hairpins with stem lengths proportional to repeat number ( 21 , 22 ) may be an oversimplification. The melting profiles of strands with 30 repeats provide an important new insight when compared with their 10 repeat counterparts. The T m values for (CTG) 30 and (CAG) 30 (Fig. 1 b) are no higher than those of (CTG) 10 and (CAG) 10 (Fig. 1 a) and their [Delta]Ho and resultant [Delta]Go for melting are only 40% higher (Table 1 ). Since T m stays the same, and [Delta]Ho (as indicated by sigmoidal slope) only increases by 40%, instead of by 200% as predicted in proportion to repeat number ( 10 , 12 ), it seems likely that 30 repeats tend to form more complex hairpin folds with stems not much longer than those formed by 10 repeats.

To our knowledge, this is the first time [Delta]Ho has been experimentally evaluated in addition to T m for single-strand folding of triplet repeats. The combination of [Delta]Ho and T m embodied in [Delta]Go provides a more reliable measure of folding stability than T m alone. The T m and [Delta]Go values for 10 repeats (Table 1 ) both indicate that folding stability in 19 mM Na + counterion decreases in the order, (CTG) 10 > (GAC) 10 [approx] (CAG) 10 > (GTC) 10 . A previous study of 15 repeats, using much less Na + (1 mM) and T m evaluated by electrophoretic mobility melting profiles ( 19 , 20 ), indicates (GTC) 15 [approx] (GAC) 15 > (CAG) 15 [approx] (GTC) 15 . The main difference is that we find T m = 47oC for (CAG) 10 , the same as for (CAG) 30 (Table 1 ), whereas the value reported for (CAG) 15 is only 38oC ( 20 ). Since T m and [Delta]Go are strongly salt-dependent, as shown in Table 1 , some disagreement may be expected when different salt concentrations are used by different methods.

We note that for CNG repeats (CNGCNG...), single-strand folding yields a hairpin stem stabilized by base pairing between GC doublets, known to have attractive base stacking interactions even in low salt ( 15 , 16 ). However, for GNC repeats (GNCGNC...), hairpin folding results in base pairing between CG doublets, which have less favorable stacking in low salt ( 15 , 16 ). As salt is added, CG stacking becomes more attractive while GC stacking stays the same, so that CG and GC doublets form equally stable stacked base pairs in ~1 M salt ( 15 ). At this time little is known about how the mispairs flanking base-paired doublets in hairpin stems affect doublet stacking interactions.

Hairpin stems of CTG or GTC repeats have pyrimidine mispairs (T opposite T) whereas those of CAG or GAC repeats have purine mispairs (A opposite A). As seen in Table 1 , the purine mispairs result in lower enthalpy changes for melting. At each salt concentration, both (CAG) 10 and (GAC) 10 show lower [Delta]Ho than (CTG) 10 and (GTC) 10 . This result is consistent with expectation, since purines are larger than pyrimidines and therefore more likely to interfere with normal stacking of GC and CG doublets. However, [Delta]Ho is not a good index of stability at physiological temperature. As pointed out in the Results section, the proper stability index ([Delta]Go at 37oC) correlates much better with T m than with [Delta]Ho. The reason is that stability is determined by the product of two terms one of which (T m - 37oC) changes much more than the other ([Delta]Ho/T m ), for equal numbers of repeats.

We see that (CTG) 10 folding is ~2 kcal/mol more stable than (GTC) 10 at low to physiological salt concentrations (Table 1 ). This is the biggest difference in stability we find for single-strands with 10 repeats. The difference is consistent with the idea that GC doublets form stronger stacked base pairs than CG doublets even when flanked by T mispairs. However, when flanked by A mispairs, the difference in stability between GC and CG doublets is no longer apparent. We find (CAG) 10 and (GAC) 10 have similar folding stability, intermediate between (CTG) 10 and (GTC) 10 . As seen in Table 1 , hairpin folds of triplet repeat sequences are much less stable than their correctly paired duplexes, their [Delta]Go values being about an order of magnitude smaller.

Upon increasing repeat number, we find (CTG) 30 and (CAG) 30 single strands have [Delta]Ho and [Delta]Go values only 1.4 times as large as their shorter counterparts with 10 repeats, and no higher T m (Table 1 ). This is perhaps the most surprising result, considering that strands with three times as many repeats might be expected to form a hairpins three times as long and therefore have [Delta]Ho and [Delta]Go about three times as large.

Models for hairpin folding by 10 triplet repeats

Shown in Figure 4 are some hairpin folding models that may explain our results for single strands with 10 repeats. To account for the high stability of (CTG) 10 folding (Fig. 1 a and Table 1 ), hairpin structure I is suggested for this strand (Fig. 4 ). This compact fold yields the longest possible duplex stem made with basepairs between GC doublets (G@C pair followed by C@G pair) and mispairs of T opposite T. Structure II for (CTG) 10 , with one less G@C bp and a short protruding `sticky end' end permitting weak dimerization, is also suggested (Fig. 4 ) to account for the minor component melting at 29oC (Fig. 1 a).


Figure 4 . Models of 30mer strand folding proposed to explain the melting profiles in Figures 1a and 2a and b and the electrophoretic results in Figure 3. Mispairs of T opposite T are shown stacked, but those of A opposite A are shown unstacked to explain differences in cooperativity observed in melting curves. Two structures (I and II) are indicated for (CTG) 10 and also for (GAC) 10 , to account for the observation of unstable interstrand associations promoted by freezing.

A compact hairpin fold similar to structure I for (CTG) 10 , but containing base-paired CG doublets instead of GC, is proposed for the (GTC) 10 strand in Figure 4 . As pointed out earlier, base pairs formed between CG doublets (C@G followed by G@C) are expected to have less favorable base stacking interactions (15,16) than those formed between GC doublets. Thus hairpins formed by GTC repeats should be less stable than those formed by CTG repeats, in agreement with Table 1 showing (GTC) 10 has lower T m and [Delta]G for melting than (CTG) 10 and also with electrophoretic mobility melting comparisons indicating (GTC) 15 is less stably folded than (CTG) 15 ( 19 ).

The corresponding triplet repeats containing A in place of T, (CAG) 10 and (GAC) 10 , show melting curves with lower slopes, implying that A mispairs interfere with stacking more than T mispairs. Previous studies of mispairing have shown that A opposite A destabilizes DNA double helices more than T opposite T ( 23 ). The bulkiness of As may cause them to remain extrahelical, keeping them exposed to water and thereby reducing cooperative or sigmoidal behavior (Figs 1 a and 2 a and b). Having mispaired As more exposed may explain why (CAG) 10 forms less stable hairpins than (CTG) 10 . Also, to explain our finding that (GAC) 10 forms interstrand dimers but (CAG) 10 does not, we suggest that the somewhat more stable (GAC) 10 folding, in addition to yielding hairpin I, may also yield a `slipped' structure II (Fig. 4 ), with a `sticky end' long enough to permit fairly stable dimerization with T m = 34oC (Fig. 2 a and b).

Comparison with theory

Theoretical calculations have been made assuming that larger numbers of CAG or CTG repeats form proportionally longer and more stable hairpins ( 10 , 12 ). According to these calculations, enthalpy and resultant free energy changes for hairpin melting should be about three times as high for 30 repeats as for 10 repeats, whereas we find them to be only 1.4 times as high (Table 1 ). Such calculations may not fully take into account the flexibility of structures containing mispairs and their potential to bend or kink into shorter domains of folding. The greater degrees of freedom available to structures with repeated mismatches may allow multiple short hairpins to form preferentially over longer ones.

The exact structure of mispairs in triplet repeat hairpins is not known, but there are some interesting theoretical possibilities. For example, in the case of CAG hairpins, with mispaired As flanked by normal C@G and G@C base pairs, a `triad' structure has been proposed ( 11 ), in which each A is unstacked and hydrogen-bonded to G on the opposite strand, while allowing G to remain hydrogen-bonded to C in the normal manner. The proposed bonding of A to G, if it occurs, apparently does not make CAG hairpins more stable than CTG hairpins which lack such bonding opportunities. On the contrary, CTG hairpins are significantly more stable, with both higher T m and larger [Delta]H o for melting than those formed by CAG repeats (Table 1 ).

Relation to slippage

The exact cause of a DNA slippage event, the number of repeats added or subtracted per event and the likelihood of such events occurring as DNA polymerase traverses a trinucleotide repeat region is unknown. Formation of a hairpin structure on the template strand might result in a deletion if polymerase can bypass the hairpin. Alternatively, hairpin formation on the elongating primer strand might result in addition of repeats. Paternal transmissions of disease alleles in pedigrees (reviewed in 27 - 30 ) and in single sperm and single molecule studies ( 27 - 30 ), reveal a preponderance of expansions suggesting that deletion caused by polymerase bypass of template hairpin structures is relatively rare. Data from single sperm typing studies on Huntington's disease germline mutations support the idea that blocks to polymerase lead to additions of a random number of repeats per slippage event ( 26 ). Alleles with more repeats have a higher germline mutation frequency and add a higher average number of repeats. Such alleles may form more complex secondary structures, e.g., several smaller hairpins instead of one long one, consistent with the physical-biochemical data presented here and lead to the addition of more repeats.

It has been suggested that the potential for slippage is greatest when the lagging strand is the template ( 27 - 30 ). The idea proposed is that when an Okazaki fragment is initiated within a repeated region, slippage of the elongating strand can occur at both its ends. Each slippage event in a triplet repeat region must involve the melting of an integral number of three basepair units in the primer/template duplex. The heat absorption required is >20 kcal/mol of base-paired CAG/CTG triplet at physiological ionic strength (15), in accordance with our finding of [Delta]Ho = 220 kcal/mol for the duplex with 10 base-paired triplets (Table 1 ). This energy requirement can only be partially offset by hairpin folding in the dissociated strand, since hairpins are much less stable than correct duplexes. However, the repeat region in the template lagging strand, rendered single stranded by the movement of the replication fork or by exonuclease activity in mismatch repair, may form hairpin structures spontaneously. The presence of stable hairpin folds in the template may contribute to primer slippage by acting as physical blocks to primer elongation by polymerase. Thus stalling of the polymerase at template folds might be a primary cause for slippage of the nascent strand leading to triplet additions.

ACKNOWLEDGEMENT

This work was supported by National Institutes of Health grants AG11398, GM21422 and GM36745.

REFERENCES

1 Ashley, C.T., Jr and Warren, S.T. (1995) Annu. Rev. Genet. 29, 703-728.

2 Sutherland, G.R. and Richards, R.J. (1995) Proc. Natl Acad. Sci. USA 92, 3636-3641. MEDLINE Abstract

3 Kunkel, T.A. (1993) Nature 365, 207-208. MEDLINE Abstract

4 Mandel, J.-L. (1994) Nature Genet. 7, 453-455. MEDLINE Abstract

5 Nelson, D.L. and Warren, S.T. (1993) Nature Genet. 4, 107-108. MEDLINE Abstract

6 Campuzano, V., Montermini, L., Molto, M.D., Pianese, L., Cossee, M., Cavalcanti, F., Monros, E., Rodius, F., Duclos, F., Monticelli, A., Zara, F., Canizares, J., Koutnikova, H., Bidichandani, S.I., Gellera, C., Filla, A., De Frutos, R., Palau, F., Patel, P.I., Di Donato, S., Mandel, J.-L., Cocozza, S., Koenig, M. and Pandolfo, M. (1996) Science 271, 1423-1427. MEDLINE Abstract

7 Streisinger, G., Okada, Y., Emrich, J., Newton, J., Tsugita, A., Terzaghi, E. and Inouye, M. (1966) Cold Spring Harbor Symp. Quant. Biol. 31, 77-84.

8 Sinden, R.R. and Wells, R.D. (1992) Curr. Opin. Biotechnol. 3, 612-622. MEDLINE Abstract

9 Fry, M. and Loeb, L.A. (1994) Proc. Natl Acad. Sci. USA 91, 4950-4954. MEDLINE Abstract

10 Mitas, M., Yi, A., Dill, J., Kamp, T.J., Chambers, E.J. and Haworth, I.S. (1995) Nucleic Acids Res. 23, 1050-1059. MEDLINE Abstract

11 Kuryavyi, V.V. and Jovin, T.M. (1995) Nature Genet. 9, 339-341. MEDLINE Abstract

12 Gacy, A.M., Goeliner, G., Juranic, N., Macura, S. and McMurray, C.T. (1995) Cell 81, 533-540. MEDLINE Abstract

13 Mitas, M., Yu, A., Dill, J. and Haworth, I.S. (1995) Biochemistry 34, 12803-12811. MEDLINE Abstract

14 Marky, L.A. and Breslauer, K.J. (1987) Biopolymers 26, 1601-1620. MEDLINE Abstract

15 Petruska, J. and Goodman, M.F. (1995) J. Biol. Chem. 270, 746-750. MEDLINE Abstract

16 Gotoh, O. and Tagashira, Y. (1981) Biopolymers 20, 1033-1042.

17 Delcourt, S.G. and Blake, R.D. (1991) J. Biol. Chem. 266, 15160-15169. MEDLINE Abstract

18 Breslauer, K.J., Frank, R., Blocker, H. and Marky, L.A. (1986) Proc. Natl Acad. Sci. USA 83, 3746-3750. MEDLINE Abstract

19 Yu, A., Dill, J., Wirth, S.S., Huang, G., Lee, V.H., Haworth, I.S. and Mitas, M. (1995) Nucleic Acids Res. 23, 2706-2714. MEDLINE Abstract

20 Yu, A., Dill, J. and Mitas, M. (1995) Nucleic Acids Res. 23, 4055-4057. MEDLINE Abstract

21 Antao, V.P. and Tinoco, I., Jr (1992) Nucleic Acids Res. 20, 819-824. MEDLINE Abstract

22 Erie, D.A., Suri, A.K., Breslauer, K.J., Jones, R.A. and Olson, W.K. (1993) Biochemistry 32, 436-454. MEDLINE Abstract

23 Aboul-ela, F., Koh, D., Tinoco, I.J. and Martin, F.H. (1985) Nucleic Acids Res. 13, 4811-4825.

24 Zhang, L., Fishbeck, K. and Arnheim, N. (1995) Hum. Mol. Genet. 4, 303-305. MEDLINE Abstract

25 Monckton, D.G., Wong, L.-J.C.A.T. and Casky, C.T. (1995) Hum. Mol. Genet. 4, 1-8. MEDLINE Abstract

26 Leeflang, E.P., Zhang, L., Tavare, S., Hubert, R., Srinidhi, J., MacDonald, M.E., Myers, R.M., de Young, M., Wexler, N.S., Gusella, J.F. and Arnheim, N. (1995) Hum. Mol. Genet. 4, 1519-1526. MEDLINE Abstract

27 Richards, R.I. and Sutherland, G.R. (1994) Nature Genet. 6, 114-116. MEDLINE Abstract

28 Kunst, C.B. and Warren, S.T. (1994) Cell 77, 853-861. MEDLINE Abstract

29 Eichler, E.E., Holden, J.J.A., Popovich, B.W., Reiss, A.L., Snow, K., hibodeau, S.N., Richards, S., Ward, P.A. and Nelson, D.L. (1994) Nature Genet. 8, 88-94. MEDLINE Abstract

30 Kang, S., Jaworski, A., Ohshima, K. and Wells, R.D. (1995) Nature Genet. 10, 213-218. MEDLINE Abstract


Return

* To whom correspondence should be addressed
Add to CiteULike CiteULike   Add to Connotea Connotea   Add to Del.icio.us Del.icio.us    What's this?


This article has been cited by other articles:


Home page
J. Biol. Chem.Home page
E. Delagoutte, G. M. Goellner, J. Guo, G. Baldacci, and C. T. McMurray
Single-stranded DNA-binding Protein in Vitro Eliminates the Orientation-dependent Impediment to Polymerase Passage on CAG/CTG Repeats
J. Biol. Chem., May 9, 2008; 283(19): 13341 - 13356.
[Abstract] [Full Text] [PDF]


Home page
Mol. Cell. Biol.Home page
Y. Lin and J. H. Wilson
Transcription-Induced CAG Repeat Contraction in Human Cells Is Mediated in Part by Transcription-Coupled Nucleotide Excision Repair
Mol. Cell. Biol., September 1, 2007; 27(17): 6209 - 6217.
[Abstract] [Full Text] [PDF]


Home page
Nucleic Acids ResHome page
E. Grabczyk, M. Mancuso, and M. C. Sammarco
A persistent RNA{middle dot}DNA hybrid formed by transcription of the Friedreich ataxia triplet repeat in live bacteria, and by T7 RNAP in vitro
Nucleic Acids Res., August 13, 2007; 35(16): 5351 - 5359.
[Abstract] [Full Text] [PDF]


Home page
GeneticsHome page
R. Zahra, J. K. Blackwood, J. Sales, and D. R. F. Leach
Proofreading and Secondary Structure Processing Determine the Orientation Dependence of CAG{middle dot}CTG Trinucleotide Repeat Instability in Escherichia coli
Genetics, May 1, 2007; 176(1): 27 - 41.
[Abstract] [Full Text] [PDF]


Home page
Microbiol. Mol. Biol. Rev.Home page
E. V. Mirkin and S. M. Mirkin
Replication Fork Stalling at Natural Impediments
Microbiol. Mol. Biol. Rev., March 1, 2007; 71(1): 13 - 35.
[Abstract] [Full Text] [PDF]


Home page
Mol Biol EvolHome page
P. Siwach, S. D. Pophaly, and S. Ganesh
Genomic and Evolutionary Insights into Genes Encoding Proteins with Single Amino Acid Repeats
Mol. Biol. Evol., July 1, 2006; 23(7): 1357 - 1369.
[Abstract] [Full Text] [PDF]


Home page
J. Biol. Chem.Home page
K. Walter, G. Warnecke, R. Bowater, W. Deppert, and E. Kim
Tumor Suppressor p53 Binds with High Affinity to CTG{middle dot}CAG Trinucleotide Repeats and Induces Topological Alterations in Mismatched Duplexes
J. Biol. Chem., December 30, 2005; 280(52): 42497 - 42507.
[Abstract] [Full Text] [PDF]


Home page
Nucleic Acids ResHome page
S. Amrane, B. Sacca, M. Mills, M. Chauhan, H. H. Klump, and J.-L. Mergny
Length-dependent energetics of (CTG)n and (CAG)n trinucleotide repeats
Nucleic Acids Res., July 21, 2005; 33(13): 4065 - 4077.
[Abstract] [Full Text] [PDF]


Home page
Nucleic Acids ResHome page
L. M. Chi and S. L. Lam
Structural roles of CTG repeats in slippage expansion during DNA replication
Nucleic Acids Res., March 14, 2005; 33(5): 1604 - 1617.
[Abstract] [Full Text] [PDF]


Home page
Nucleic Acids ResHome page
L. H. Mochmann and R. D. Wells
Transcription influences the types of deletion and expansion products in an orientation-dependent manner from GAC*GTC repeats
Nucleic Acids Res., August 18, 2004; 32(15): 4469 - 4479.
[Abstract] [Full Text] [PDF]


Home page
Nucleic Acids ResHome page
B. L. Heidenfelder and M. D. Topal
Effects of sequence on repeat expansion during DNA replication
Nucleic Acids Res., December 15, 2003; 31(24): 7159 - 7164.
[Abstract] [Full Text] [PDF]


Home page
Nucleic Acids ResHome page
M.-H. Hou, H. Robinson, Y.-G. Gao, and A. H.-J. Wang
Crystal structure of actinomycin D bound to the CTG triplet repeat sequences linked to neurological diseases
Nucleic Acids Res., November 15, 2002; 30(22): 4910 - 4917.
[Abstract] [Full Text] [PDF]


Home page
Proc. Natl. Acad. Sci. USAHome page
J. Volker, N. Makube, G. E. Plum, H. H. Klump, and K. J. Breslauer
Conformational energetics of stable and metastable states formed by DNA triplet repeat oligonucleotides: Implications for triplet expansion diseases
PNAS, November 12, 2002; 99(23): 14700 - 14705.
[Abstract] [Full Text] [PDF]


Home page
J. Biol. Chem.Home page
M. J. Hartenstine, M. F. Goodman, and J. Petruska
Weak Strand Displacement Activity Enables Human DNA Polymerase beta to Expand CAG/CTG Triplet Repeats at Strand Breaks
J. Biol. Chem., October 25, 2002; 277(44): 41379 - 41389.
[Abstract] [Full Text] [PDF]


Home page
Nucleic Acids ResHome page
C. E. Pearson, M. Tam, Y.-H. Wang, S. E. Montgomery, A. C. Dar, J. D. Cleary, and K. Nichol
Slipped-strand DNAs formed by long (CAG){middle dot}(CTG) repeats: slipped-out repeats and slip-out junctions
Nucleic Acids Res., October 15, 2002; 30(20): 4534 - 4547.
[Abstract] [Full Text] [PDF]


Home page
J. Biol. Chem.Home page
L. A. Henricksen, J. Veeraraghavan, D. R. Chafin, and R. A. Bambara
DNA Ligase I Competes with FEN1 to Expand Repetitive DNA Sequences in Vitro
J. Biol. Chem., June 14, 2002; 277(25): 22361 - 22369.
[Abstract] [Full Text] [PDF]


Home page
Nucleic Acids ResHome page
G. Raca, E. Yu. Siyanova, C. T. McMurray, and S. M. Mirkin
Expansion of the (CTG)n repeat in the 5'-UTR of a reporter gene impedes translation
Nucleic Acids Res., October 15, 2000; 28(20): 3943 - 3949.
[Abstract] [Full Text] [PDF]


Home page
Genome ResHome page
L. N. Bull, C. R. Pabón-Peña, and N. B. Freimer
Compound Microsatellite Repeats: Practical and Theoretical Features
Genome Res., September 1, 1999; 9(9): 830 - 838.
[Abstract] [Full Text]


Home page
Proc. Natl. Acad. Sci. USAHome page
H. Moore, P. W. Greenwell, C.-P. Liu, N. Arnheim, and T. D. Petes
Triplet repeats form secondary structures that escape DNA repair in yeast
PNAS, February 16, 1999; 96(4): 1504 - 1509.
[Abstract] [Full Text] [PDF]


Home page
J. Biol. Chem.Home page
R. R. Iyer and R. D. Wells
Expansion and Deletion of Triplet Repeat Sequences in Escherichia coli Occur on the Leading Strand of DNA Replication
J. Biol. Chem., February 5, 1999; 274(6): 3865 - 3877.
[Abstract] [Full Text] [PDF]


Home page
J. Biol. Chem.Home page
J. Petruska, M. J. Hartenstine, and M. F. Goodman
Analysis of Strand Slippage in DNA Polymerase Expansions of CAG/CTG Triplet Repeats Associated with Neurodegenerative Disease
J. Biol. Chem., February 27, 1998; 273(9): 5204 - 5210.
[Abstract] [Full Text] [PDF]


Home page
J. Biol. Chem.Home page
M. Napierala and W. J. Krzyzosiak
CUG Repeats Present in Myotonin Kinase RNA Form Metastable "Slippery" Hairpins
J. Biol. Chem., December 5, 1997; 272(49): 31079 - 31085.
[Abstract] [Full Text] [PDF]


Home page
J. Biol. Chem.Home page
K. Ohshima and R. D. Wells
Hairpin Formation during DNA Synthesis Primer Realignment in Vitro in Triplet Repeat Sequences from Human Hereditary Disease Genes
J. Biol. Chem., July 4, 1997; 272(27): 16798 - 16806.
[Abstract] [Full Text] [PDF]


Home page
Proc. Natl. Acad. Sci. USAHome page
R. Cox and S. M. Mirkin
Characteristic enrichment of DNA repeats in different genomes
PNAS, May 13, 1997; 94(10): 5237 - 5242.
[Abstract] [Full Text] [PDF]


Home page
J. Biol. Chem.Home page
M. J. Hartenstine, M. F. Goodman, and J. Petruska
Base Stacking and Even/Odd Behavior of Hairpin Loops in DNA Triplet Repeat Slippage and Expansion with DNA Polymerase
J. Biol. Chem., June 9, 2000; 275(24): 18382 - 18390.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Print PDF (140K) Freely available
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Right arrow Search for citing articles in:
ISI Web of Science (77)
Right arrowRequest Permissions
Right arrow Commercial Re-use Guidelines
for Open Access NAR Content
Google Scholar
Right arrow Articles by Petruska, J
Right arrow Articles by Goodman, M.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Petruska, J
Right arrow Articles by Goodman, M.
Social Bookmarking
 Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?