Evolution of the cystatin B gene: implications for the origin of its variable dodecamer tandem repeat in humans

The human cystatin B gene contains a variable number of 12-bp tandem repeats in its promoter region, of which the common alleles contain two or three copies and unusual expansion causes progressive myoclonus epilepsy of the Unverricht-Lundborg type. We undertook a comprehensive analysis of the genomic sequence to address the evolutionary events of this variable repeat. By examination of a contiguous genome sequence spanning 5.0 kb and linkage analysis of detected polymorphic changes, we identified six major intragenic haplotypes in unrelated Japanese subjects. The number of normal repeats was closely correlated with these alleles, indicating that changes in the array should be comparatively rare events during human evolution. To examine the origin of the repeat array further, we also analyzed five primate genomes. Repetitive polymorphism was unlikely in hominoids, and the array originated with the dodecamer itself in the course of primate evolution. The variability conceivably developed after the separation to humans.     

Polymerase chain reaction-based analysis using deaminated DNA of dodecamer expansions in CSTB, associated with Unverricht-Lundborg myoclonus epilepsy

Progressive myoclonus epilepsy of the Unverricht-Lundborg type is an autosomal recessive disorder that is characterized clinically by myoclonic seizures and ataxia. The majority of affected individuals carry repeat expansions of a dodecamer in the promoter region of the cystatin B gene. The unusually high GC content of this tract is refractory to conventional polymerase chain reaction (PCR), and, as a result, a circumventive procedure involving the deamination of DNA with sodium bisulfite has been proposed. This study evaluates the effectiveness of this deamination modification for the detection of dodecamer repeat variants. An analysis of 258 healthy Japanese individuals revealed an allele with four copies of the dodecamer repeat with a frequency of 0.01, in addition to the more commonly observed two and three copy repeat alleles. Homozygous repeat expansions 600 and 680 base pairs in length were detected in the analyses of two affected individuals. For these cases, sequencing, along with an alternative PCR-stutter formation, revealed 41 and 48 copies, respectively, of the dodecamer repeat. The complete conversion of C to T was observed in the expanded tracts, indicating that no methylation occurred at the CpG sites. Based on these results, it was concluded that the use of deaminated DNA allows for a precise analysis of consecutive GC tracts.     

Progressive myoclonus epilepsy [EPM1] repeat d(CCCCGCCCCGCG)n forms folded hairpin structures at physiological pH.

The secondary structure of DNA has been shown to be an important component in the mechanism of expansion of the trinucleotide repeats that are associated with many neurodegenerative disorders. Recently, expansion of a dodecamer repeat, (CCCCGCCCCGCG)n upstream of cystatin B gene has been shown to be the most common mutation associated with Progressive Myoclonus Epilepsy (EPM1) of Unverricht-Lundborg type. We have investigated structure of oligonucleotides containing one, two and three copies of the EPM1 repeat sequences at physiological pH. CD spectra and anomalous faster gel electrophoretic mobilty indicates formation of intramolecularly folded structures that are formed independent of concentration. Hydroxylamine probing allowed us to identify the C residues that are involved in C.G base pairing. P1 nuclease studies elucidated the presence of unpaired regions in the folded back structures. UV melting studies show biphasic melting curves for the oligonucleotides containing two and three EPM1 repeats. Our data suggests multiple hairpin structures for two and three repeat containing oligonucleotides. In this paper we show that oligonucleotides containing EPM1 repeat adopt secondary structures that may facilitate strand slippage thereby causing the expansion.     

Dentatorubral-pallidoluysian atrophy (DRPLA): Close correlation of CAG repeat expansions with the wide spectrum of clinical presentations and prominent anticipation

Dentatorubral-pallidoluysian atrophy (DRPLA) is a rare autosomal dominant neurodegenerative disease characterized by various combinations of ataxia, choreoathetosis, myoclonus, epilepsy and dementia as well as various ages of onset. We have identified a specific unstable trinucleotide repeat expansion in a gene on the short arm of chromosome 12 as the pathogenic mutation for DRPLA. We investigated how the degree of the expansion of the CAG repeat affects the clinical manifestations of DRPLA. The sizes of the expanded alleles were well correlated with the ages of onset (r = −0.6955, P < 0.001). Patients with progressive myoclonus epilepsy (PME) phenotype had larger expansions (62–79 repeats) and earlier ages of onset (onset before age 20). Furthermore, most of the patients with PME phenotype inherited their expanded alleles from their affected fathers. On the other hand, patients with non-PME phenotype showed later ages of onset (onset after age 20) and smaller expansions (54–67 repeats). When ages of onset of each clinical symptom are compared with sizes of the CAG repeat, there is again a remarkably high correlation of the sizes of CAG repeat with each of the clinical symptoms. Thus the wide variation in clinical manifestations of DRPLA can now be clearly explained based on the degree of CAG repeat expansion, which strongly indicates that the expanded alleles are intimately involved in the neuronal degeneration in dentatofugal and pallidofugal systems.     

Identification of mutations in cystatin B, the gene responsible for the Unverricht-Lundborg type of progressive myoclonus epilepsy (EPM1).

Progressive myoclonus epilepsy (EPM1) is an autosomal recessive disorder, characterized by severe, stimulus-sensitive myoclonus and tonic-clonic seizures. The EPM1 locus was mapped to within 0.3 cM from PFKL in chromosome 21q22.3. The gene for the proteinase inhibitor cystatin B was recently localized in the EPM1 critical region, and mutations were identified in two EPM1 families. We have identified six nucleotide changes in the cystatin B gene of non-Finnish EPM1 families from northern Africa and Europe. The 426G-->C change in exon 1 results in a Gly4Arg substitution and is the first missense mutation described that is associated with EPM1. Molecular modeling predicts that this substitution severely affects the contact of cystatin B with papain. Mutations in the invariant AG dinucleotides of the acceptor sites of introns 1 and 2 probably result in abnormal splicing. A deletion of two nucleotides in exon 3 produces a frameshift and truncates the protein. Therefore, these four mutations are all predicted to impair the production of functional protein. These mutations were found in 7 of the 29 unrelated EPM1 patients analyzed, in homozygosity in 1, and in heterozygosity in the others. The remaining two sequence changes, 431G-->T and 2575A-->G, probably represent polymorphic variants. In addition, a tandem repeat in the 5' UTR (CCCCGCCCCGCG) is present two or three times in normal alleles. It is peculiar that in the majority of patients no mutations exist within the exons and splice sites of the cystatin B gene.     

Intramolecular i-motif structure at acidic pH for progressive myoclonus epilepsy (EPM1) repeat d(CCCCGCCCCGCG)n.

The most common mutation associated with Progressive Myoclonus Epilepsy (EPM1) of Unverricht-Lundberg type is the expansion of a dodecamer repeat, d(CCCCGCCCCGCG)n. We show that the C-rich strand of this repeat (2-3 copies) forms intercalated i-motif structure at acidic pH as judged by CD spectroscopy and anomalous gel electrophoretic mobility. The stability of the structure increases with the increase in the length of the repeat. Transient formation of stable, folded back structure like i-motif could play an important role in the mechanism of expansion of this repeat.     

Molecular dynamics studies of trinucleotide repeat DNA involved in neurodegenerative disorders.

Expansion of trinucleotide repeat DNA of the classes CAG-CTG, CGG-CCG and GAA-TTC are found to be associated with several neurodegenerative disorders. Different mechanisms have been attributed to the expansion of triplets, mainly involving the formation of alternate secondary structures by such repeats. This paper reports the molecular dynamics simulation of triplet repeat DNA sequences to study the basic structural features of DNA that are responsible for the formation of structures such as hairpins and slip-strand DNA leading to expansion. All the triplet repeat sequences studied were found to be more flexible compared to the control sequence unassociated with disease. Moreover, flexibility was found to be in the order CAG-CTG > CGG-CCG approximately GAA-TTC, the highly flexible CAG-CTG repeat being the most common cause of neurodegenerative disorders. In another simulation, a single G-C to T-A mutation at the 9th position of the CAG-CTG repeat exhibited a reduction in bending compared to the pure 15-mer CAG-CTG repeat. EPM1 dodecamer repeat associated with the pathogenesis of progressive myoclonus epilepsy was also simulated and showed flexible nature suggesting a similar expansion mechanism.     

Tetraplex formation by the progressive myoclonus epilepsy type-1 repeat: implications for instability in the repeat expansion diseases

The repeat expansion diseases are a group of genetic disorders resulting from an increase in size or expansion of a specific array of tandem repeats. It has been suggested that DNA secondary structures are responsible for this expansion. If this is so, we would expect that all unstable repeats should form such structures. We show here that the unstable repeat that causes progressive myoclonus epilepsy type-1 (EPM1), like the repeats associated with other diseases in this category, forms a variety of secondary structures. However, EPM1 is unique in that tetraplexes are the only structures likely to form in long unpaired repeat tracts under physiological conditions.     

CTG trinucleotide repeat "big jumps": large expansions, small mice

Trinucleotide repeat expansions are the genetic cause of numerous human diseases, including fragile X mental retardation, Huntington disease, and myotonic dystrophy type 1. Disease severity and age of onset are critically linked to expansion size. Previous mouse models of repeat instability have not recreated large intergenerational expansions ("big jumps"), observed when the repeat is transmitted from one generation to the next, and have never attained the very large tract lengths possible in humans. Here, we describe dramatic intergenerational CTG*CAG repeat expansions of several hundred repeats in a transgenic mouse model of myotonic dystrophy type 1, resulting in increasingly severe phenotypic and molecular abnormalities. Homozygous mice carrying over 700 trinucleotide repeats on both alleles display severely reduced body size and splicing abnormalities, notably in the central nervous system. Our findings demonstrate that large intergenerational trinucleotide repeat expansions can be recreated in mice, and endorse the use of transgenic mouse models to refine our understanding of triplet repeat expansion and the resulting pathogenesis.     

The relationship between trinucleotide (GAA) repeat length and clinical features in Friedreich ataxia.

Friedreich ataxia (FA) is associated with the expansion of a GAA trinucleotide repeat in the first intron of the X25 gene. We found both alleles expanded in 67 FA patients from 48 Italian families. Five patients from three families were compound heterozygotes with expansion on one allele and an isoleucine-->phenylalanine change at position 154 on the other one. We found neither expansions nor point mutations in three patients. The length of FA alleles ranged from 201 to 1,186 repeat units, with no overlap with the normal range, and showed a negatively skewed distribution with a peak between 800 and 1,000 repeats. The FA repeat showed meiotic instability with a median variation of 150 repeats. The lengths of both larger and smaller alleles in each patient inversely correlated with age at onset of the disorder. Smaller alleles showed the best correlation, accounting for approximately 50% of the variation of age at onset. Mean allele length was significantly higher in patients with diabetes and in those with cardiomyopathy.     

Cystatin B and its EPM1 mutants are polymeric and aggregate prone in vivo

Progressive myoclonus epilepsy type 1 (EPM1) is a neurodegenerative disease correlating with mutations of the cystatin B gene. Cystatin B is described as a monomeric protein with antiprotease function. This work shows that, in vivo, cystatin B has a polymeric structure, highly resistant to SDS, urea, boiling and sensitive to reducing agents and alkaline pH. Hydrogen peroxide increases the polymeric structure of the protein. Mass spectrometry analysis shows that the only component of the polymers is cystatin B. EPM1 mutants of cystatin B transfected in cultured cells are also polymeric. The banding pattern generated by a cysteine-minus mutant is different from that of the wild-type protein as it contains only monomers, dimers and some very high MW bands while misses components of MW intermediate between 25 and 250 kDa. Overexpression of wild-type or EPM1 mutants of cystatin B in neuroblastoma cells generates cytoplasmic aggregates. The cysteine-minus mutant is less prone to the formation of inclusion bodies. We conclude that cystatin B in vivo has a polymeric structure sensitive to the redox environment and that overexpression of the protein generates aggregates. This work describes a protein with a physiological role characterized by highly stable polymers prone to aggregate formation in vivo.     

Gonosomal mosaicism in myotonic dystrophy patients: involvement of mitotic events in (CTG)n repeat variation and selection against extreme expansion in sperm.

Myotonic dystrophy (DM) is caused by abnormal expansion of a polymorphic (CTG)n repeat, located in the DM protein kinase gene. We determined the (CTG)n repeat lengths in a broad range of tissue DNAs from patients with mild, classical, or congenital manifestation of DM. Differences in the repeat length were seen in somatic tissues from single DM individuals and twins. Repeats appeared to expand to a similar extent in tissues originating from the same embryonal origin. In most male patients carrying intermediate- or small-sized expansions in blood, the repeat lengths covered a markedly wider range in sperm. In contrast, male patients with large allele expansions in blood (> 700 CTGs) had similar or smaller repeats in sperm, when detectable. Sperm alleles with > 1,000 CTGs were not seen. We conclude that DM patients can be considered gonosomal mosaics, i.e., combined somatic and germ-line tissue mosaics. Most remarkably, we observed multiple cases where the length distributions of intermediate- or small-sized alleles in fathers'' sperm were significantly different from that in their offspring''s blood. Our combined findings indicate that intergenerational length changes in the unstable CTG repeat are most likely to occur during early embryonic mitotic divisions in both somatic and germ-line tissue formation. Both the initial CTG length, the overall number of cell divisions involved in tissue formation, and perhaps a specific selection process in spermatogenesis may influence the dynamics of this process. A model explaining mitotic instability and sex-dependent segregation phenomena in DM manifestation is discussed.     

Stabilizing effects of interruptions on trinucleotide repeat expansions in Saccharomyces cerevisiae

In most trinucleotide repeat (TNR) diseases, the primary factor determining the likelihood of expansions is the length of the TNR. In some diseases, however, stable alleles contain one to three base pair substitutions that interrupt the TNR tract. The unexpected stability of these alleles compared to the frequent expansions of perfect TNRs suggested that interruptions somehow block expansions and that expansions occur only upon loss of at least one interruption. The work in this study uses a yeast genetic assay to examine the mechanism of stabilization conferred by two interruptions of a 25-repeat tract. Expansion rates are reduced up to 90-fold compared to an uninterrupted allele. Stabilization is greatest when the interruption is replicated early on the lagging strand, relative to the rest of the TNR. Although expansions are infrequent, they are often polar, gaining new DNA within the largest available stretch of perfect repeats. Surprisingly, interruptions are always retained and sometimes even duplicated, suggesting that expansion in yeast cells can proceed without loss of the interruption. These findings support a stabilization model in which interruptions contribute in cis to reduce hairpin formation during TNR replication and thus inhibit expansion rates.     

PolyQ repeat expansions in ATXN2 associated with ALS are CAA interrupted repeats

Amyotrophic lateral sclerosis (ALS) is a devastating, rapidly progressive disease leading to paralysis and death. Recently, intermediate length polyglutamine (polyQ) repeats of 27-33 in ATAXIN-2 (ATXN2), encoding the ATXN2 protein, were found to increase risk for ALS. In ATXN2, polyQ expansions of ≥ 34, which are pure CAG repeat expansions, cause spinocerebellar ataxia type 2. However, similar length expansions that are interrupted with other codons, can present atypically with parkinsonism, suggesting that configuration of the repeat sequence plays an important role in disease manifestation in ATXN2 polyQ expansion diseases. Here we determined whether the expansions in ATXN2 associated with ALS were pure or interrupted CAG repeats, and defined single nucleotide polymorphisms (SNPs) rs695871 and rs695872 in exon 1 of the gene, to assess haplotype association. We found that the expanded repeat alleles of 40 ALS patients and 9 long-repeat length controls were all interrupted, bearing 1-3 CAA codons within the CAG repeat. 21/21 expanded ALS chromosomes with 3CAA interruptions arose from one haplotype (GT), while 18/19 expanded ALS chromosomes with <3CAA interruptions arose from a different haplotype (CC). Moreover, age of disease onset was significantly earlier in patients bearing 3 interruptions vs fewer, and was distinct between haplotypes. These results indicate that CAG repeat expansions in ATXN2 associated with ALS are uniformly interrupted repeats and that the nature of the repeat sequence and haplotype, as well as length of polyQ repeat, may play a role in the neurological effect conferred by expansions in ATXN2.     

Sex-dependent mechanisms for expansions and contractions of the CAG repeat on affected Huntington disease chromosomes.

A total of 254 affected parent-child pairs with Huntington disease (HD) and 440 parent-child pairs with CAG size in the normal range were assessed to determine the nature and frequency of intergenerational CAG changes in the HD gene. Intergenerational CAG changes are extremely rare (3/440 [0.68%]) on normal chromosomes. In contrast, on HD chromosomes, changes in CAG size occur in approximately 70% of meioses on HD chromosomes, with expansions accounting for 73% of these changes. These intergenerational CAG changes make a significant but minor contribution to changes in age at onset (r2 = .19). The size of the CAG repeat influenced larger intergenerational expansions (> 7 CAG repeats), but the likelihood of smaller expansions or contractions was not influenced by CAG size. Large expansions (> 7 CAG repeats) occur almost exclusively through paternal transmission (0.96%; P < 10(-7)), while offspring of affected mothers are more likely to show no change (P = .01) or contractions in CAG size (P = .002). This study demonstrates that sex of the transmitting parent is the major determinant for CAG intergenerational changes in the HD gene. Similar paternal sex effects are seen in the evolution of new mutations for HD from intermediate alleles and for large expansions on affected chromosomes. Affected mothers almost never transmit a significantly expanded CAG repeat, despite the fact that many have similar large-sized alleles, compared with affected fathers. The sex-dependent effects of major expansion and contractions of the CAG repeat in the HD gene implicate different effects of gametogenesis, in males versus females, on intergenerational CAG repeat stability.     

Phenotype correlation and intergenerational dynamics of the Friedreich ataxia GAA trinucleotide repeat.

The Friedreich ataxia (FA) mutation has recently been identified as an unstable trinucleotide GAA repeat present 7-22 times in the normal population but amplified as many as > 1,000 times in FA. Since it is an autosomal recessive disease, FA does not show typical features observed in other dynamic mutation disorders, such as genetic anticipation. We have analyzed the GAA repeat in 104 FA patients and 163 carrier relatives previously defined by linkage analysis. The GAA expansion was detected in all patients, most (94%) of them being homozygous for the mutation. We have demonstrated that clinical variability in FA is related to the size of the expanded alleles: milder forms of the disease-late-onset FA and FA with retained reflexes-are associated with shorter expansions, especially with the smaller of the two expanded alleles. Absence of cardiomyopathy is also associated with shorter alleles. Dynamics of the GAA repeat has been investigated in 212 parent-offspring pairs. Meiotic instability showed a sex bias: paternally transmitted alleles tend to decrease in a linear way that depends on the paternal expansion size, whereas maternal alleles can either increase or decrease. A different pattern of intergenerational variation was also observed, depending on the genetic status of the sib: patients had shorter expansions than were seen in heterozygous carriers. This finding has been interpreted as a postzygotic event. Finally, we have observed that the size of the expansion remains constant in the population through carriers.