|
|
||||||||
1From the Departments of Human Genetics and 4Ophthalmology, Radboud University Nijmegen Medical Centre, Nijmegen, The Netherlands; the 2Nijmegen Centre for Molecular Life Sciences, Nijmegen, The Netherlands; 3Canisius-Wilhelmina Ziekenhuis, Nijmegen, The Netherlands; the 5Department of Clinical Genetics, Central Manchester and Manchester Childrens University Hospitals NHS Trust, Manchester, United Kingdom; and the 6Department of Ophthalmology, University Medical Centre Utrecht, Utrecht, The Netherlands.
| Abstract |
|---|
|
|
|---|
G) was recently reported that resulted in a 39-nucleotide exon 8 in-frame deletion. We investigated whether CSPG2/Versican was mutated in six Dutch families and one Chinese family with Wagner disease and in a family with ERVR. METHODS. In all families, extensive ophthalmic examinations, haplotype analysis of the 5q14.3 region, and sequence analysis of CSPG2/Versican were performed. The effects of splice site mutations were assessed by reverse transcriptionpolymerase chain reaction (RT-PCR) and real-time quantitative RT-PCR (QPCR).
RESULTS. Three novel intron 7 sequence variants (c.4004-5T
C, c.4004-5T
A, c.4004-1G
A) were identified in seven families. The c.4004-5T
C variant was identified in four families with Wagner disease and a family with ERVR. The families were shown to carry the same 5q14.3 haplotype, strongly suggesting that this is a common Dutch founder variant. All three changes segregated with the disease in the respective families and were absent in 250 healthy individuals. In patients with the c.4004-5T
A and c.4004-1G
A variants, RT-PCR analysis of CSPG2/Versican showed activation of a cryptic splice site resulting in a 39-nt exon 8 in-frame deletion in splice variant V0. QPCR revealed a highly significant (P < 0.0001) and consistent increase of the V2 (>38-fold) and V3 (>12-fold) splice variants in all patients with intron 7 nucleotide changes and in a Chinese Wagner disease family, in which the genetic defect remains to be found.
CONCLUSIONS. Wagner disease and ERVR are allelic disorders. Seven of the eight families exhibit a variant in intron 7 of CSPG2/Versican. The conspicuous clustering of sequence variants in the splice acceptor site of intron 7 and the consistent upregulation of the V2 and V3 isoforms strongly suggest that Wagner disease and ERVR may belong to a largely overlooked group of diseases that are caused by mRNA isoform balance shifts, representing a novel disease mechanism.
Linkage studies assigned the Wagner disease gene to a 2-cM interval at 5q14.3,4 5 but until recently, no definite causal mutations were identified. The critical region for a genetic defect underlying ERVR was found to overlap the critical region for Wagner disease, suggesting that these diseases are allelic.4 Recently, a heterozygous splice site change in CSPG2/Versican (c.4004-2A
G) was reported in a Japanese family with Wagner disease.6 This variant was shown to result in deletion of 13 amino acids from the beginning of the glycosaminoglycan (GAG)-ß domain of the mature protein.6
CSPG2/Versican has been an attractive candidate for Wagner disease and ERVR based on its genomic location and function. It encodes a large extracellular matrix proteoglycan (chondroitin sulfate proteoglycan type 2) that is present in virtually every human tissue, including different parts of the eye.7 CSPG2/Versican is believed to maintain the structure of the vitreous body in the human eye by keeping the collagen molecules apart.8 It has a tridomain structure, and the aminoterminal end binds to hyaluronan. The carboxyl-terminal domain has a C-type lectin domain adjacent to two epidermal growth factor domains and a complement regulatory region. The central area of CSPG2/Versican is encoded by two large exons (exons 7 and 8) and contains chondroitin sulfate attachment sites. The chondroitin sulfate chains of CSPG2/Versican are highly negatively charged and contribute to its antiadhesive properties. Alternative splicing of these two exons results in four different isoforms of the mature protein (see Fig. 6). V0 contains exons 7 and 8, V1 lacks exon 7 but contains exon 8, V2 contains exon 7 but lacks exon 8, and V3 lacks both these exons.7 Accordingly, the potential number of glycosoaminoglycan (GAG) attachment sites for each form are V0, 17 to 23; V1, 12 to 15; V2, 5 to 8; and V3, 0.7
In this study, we show that most Dutch families with Wagner disease and a family with ERVR share a common 5q14.3 founder haplotype and a CSPG2/Versican splice site variant. In two additional unrelated families with Wagner disease we found two other sequence variants in the same intron 7/exon 8 splice junction. We found a qualitative effect on CSPG2/Versican mRNA splice isoforms (i.e., a 39-nt in-frame truncation, in two families with Wagner disease). In addition, in all families, we found a significant increase of CSPG2/Versican mRNA isoforms V2 and V3 which leads us to propose a novel mechanism for disease pathogenesis.
| Methods |
|---|
|
|
|---|
|
Genotyping and Linkage Analysis
Genotyping at the Wagner disease/ERVR locus on chromosome 5 was performed by selecting 13 microsatellite markers from the Genome Database (www.gdb.org/ provided by RTI International, NC) covering a region of approximately 3 cM. For all the markers except D5S626, genotyping PCR was performed with primers with a fluorescent label and the results were analyzed (GeneMapper; Applied Biosystems, Inc. [ABI], IJssel, The Netherlands). For D5S626, 32P end-labeled primers were used. Haplotypes were constructed by Cyrillic 2.1 (Cyrillic Software, Wallingford, UK), and checked manually for confirmation of cosegregation of the marker alleles with the disease phenotype.
Sequence Analysis
The open reading frame (ORF) and the exonintron junctions of XRCC4 was analyzed for a proband of family W95-131, as described previously.5 The ORF and the exonintron junctions of the 15 exons of CSPG2/Versican were sequenced in a proband from each of the eight families (primer sequences available on request). Polymerase chain reaction (PCR) was performed according to the standard protocols.10 PCR products were purified (QiaQuick columns; Qiagen), either directly or after excision from the gel and were directly sequenced using dye termination chemistry (model 3700; ABI).
RT-PCR and Real-Time Quantitative RT-PCR Analysis
Total RNA was isolated from 10 mL of EDTA-anticoagulated venous blood within 4 hours after sampling (RNeasy minikit; Qiagen) according to the manufacturers protocol. cDNA was synthesized from 2.0 µg total RNA using random hexamers (GE Healthcare, Roosendaal, The Netherlands), and M-MLV reverse transcriptase (Invitrogen, Groningen, The Netherlands) in a total volume of 85.0 µL. After the synthesis, the cDNA was purified on PCR purification columns (QIAquick; Qiagen) and eluted in 30 µL sterilized water (milliQ; Millipore, Bedford, MA).
RT-PCR was performed using primer pairs that specifically amplify each of the four mRNA splice variants (see Fig. 6).6 In addition, we used a common primer pair that amplifies all mRNA variants, termed Vt for total expression (Supplementary Table S1, available online at http://www.iovs.org/cgi/content/full/47/8/3565/DC1).
For QPCR, 10x dilutions of the cDNA were made and 5.0 µL from the diluted stock was used for each PCR reaction. Primers were designed using Primer 3 software (http://frodo.wi.mit.edu/cgi-bin/primer3/primer3_www.cgi). Amplicon lengths were typically kept between 80 to 100 bp (Supplementary Table S1, http://www.iovs.org/cgi/content/full/47/8/3565/DC1). Quantitation was performed using relative quantification strategies (SYBR Green), according to the manufacturers protocol (Bio-Rad, Hercules, CA). PCR reactions were performed on a thermocycler (iCycler MyiQ Single-Color Real-Time Detection System; Bio-Rad). Initially, the primer conditions were optimized and a standard curve was determined for six different dilution points (10x320x) until for each of the primers PCR efficiencies close to 100% were obtained. For calculation purpose all the PCR efficiencies were assumed to be exactly 100%. For each sample triplicate PCR reactions were performed to ensure reproducibility. The differences in expression between patients and controls were calculated by using the 2
Ct method.11 12 By definition, Ct is the required number of amplification cycles to reach the threshold fluorescence level (automatically determined);
Ct implies the difference in the Ct between the reference gene (GUSB) and the gene of interest (CSPG2/Versican) and 
Ct represents the difference between two
Cts. GUSB was used as a reference gene, because of its stable expression in lymphocytes.13 14 
Ct is the difference between the mean
Ct of eight control subjects and the
Ct for individual patients. For statistical analysis of the QPCR data, the probabilities were calculated in the P-value calculator (www.graphpad.com/quickcalcs/PValue1.cfm/GraphPad, Inc., San Diego, CA) available in the public domain.
Splice-Site Score
The predicted effect of the identified intronic variants on splicing were estimated at the server of the Childrens Mercy Hospitals and Clinics at the University of Missouri-Kansas City (https://splice.cmh.edu).15
| Results |
|---|
|
|
|---|
|
|
120 individuals of the families W95-023, W95-038, W95-124, W95-131, W04-153, and W05-088, most of which are depicted in Figure 1 , employing 13 microsatellite markers. The phenotypes in these families fully cosegregated with marker alleles of the 5q14 region (data not shown) suggesting genetic homogeneity for Wagner disease. The haplotypes of patients carrying critical recombinations in 5q14 are depicted in Figure 3 . Recombinations observed in the affected individuals from families W95-023 and W95-131 restricted the region to approximately 850 kb with D5S626 and D5S107 as the flanking markers. All patients from four Wagner disease families (W95-023, W95-131, W04-153, W05-088) and from the ERVR family (W95-124) showed the same haplotype for the critical region, suggesting the presence of a founder allele in The Netherlands.
|
C change in CSPG2/Versican (Fig. 4) . The same variant was identified in the DNA of all patients from the three other pedigrees (W95-023, W04-153, and W05-088) belonging to the Wagner disease superfamily, and none of the unaffected individuals (Fig. 1) . In addition, we found a c.4004-5T
A change in the proband of family W95-038 and a c.4004-1G
A change in a proband from family W95-137 (Fig. 4) which also showed cosegregation with the disease in the respective families (Fig. 1) .
|
C change as observed in the other related families with Wagner disease (Fig. 1) . Sequence analysis of the intron 7exon 8 splice junction in 250 ethnically matched Dutch healthy individuals did not reveal any of these changes. These results strongly suggest the presence of a founder allele (c.4004-5T
C) in most of the Dutch patients with Wagner disease and in a family with ERVR, a clinically related phenotype. The causal variant in the Chinese family with Wagner disease has not yet been identified.
Effect of Intron 7 Variants on CSPG2/Versican Transcripts
The c.4004-1G
A variant affects one of two completely conserved nucleotides of the canonical 3' splice site of intron 7 and therefore is predicted to have a major effect on CSPG2/Versican splicing. The c.4004-5T
C and c.4004-5T
A variants are predicted to retain 64% and 14% of the original binding capacity, respectively (https://splice.cmh.edu). To investigate their effect on splicing, both qualitative and real-time quantitative RT-PCR analyses were subsequently performed.
The c.4004-1G
A and c.4004-5T
A nucleotide variants result in the activation of a cryptic downstream splice acceptor site of exon 8. V0 shows the wild-type and the 39-nt shortened cDNA fragments (Fig. 5 , lanes 1, 2, and 4), mimicking the previously described result.6 For V1, this effect is faintly visible. The c.4004-5T
C variant has only a small effect on the downstream splice acceptor (Fig. 5 , lane 3). RNA from both lymphocytes and fibroblasts show the same effect on V0 and V1 for c.4004-1G
A mutation (Fig. 5 , lanes 1 and 2). RT-PCR analysis of RNA from a patient of the Chinese family with Wagner disease (W05-282) did not show an activation of the downstream cryptic splice acceptor site (Fig. 5 , lane 5). For the ERVR family the RT-PCR result was similar to that of W95-131 (data not shown).
|
Ct and the x-fold changes (2
Ct), respectively, for the patients with Wagner disease and ERVR. The V2 isoform was upregulated between 38- and 478-fold (P < 0.0001), and the V3 isoform was upregulated between 12- and 52-fold (P < 0.0001). A 40-fold upregulation in patients indicates that the threshold fluorescence is reached approximately five amplification cycles earlier than that in control subjects. The patients from family W95-131 (Wagner disease) and W95-124 (ERVR), carrying the c.4004-5T
C change, showed comparable upregulation for V2 and V3 (Table 4) . Patients W95-137P1 and P2, carrying the c.4004-1G
A change, as well as the two Chinese patients with Wagner disease (W05-282P1 and -P2), showed the highest upregulation of V2 and V3. The V0 and V1 isoforms showed a tendency to be downregulated (Table 4) .
|
|
|
| Discussion |
|---|
|
|
|---|
C) in one of the probands we showed that this variant was present in all affected individuals of these four families with Wagner disease and the ERVR family. In addition, we found c.4004-5T
A and c.4004-1G
A variants in two unrelated families with Wagner disease.
Detailed RNA analysis of patients from these families, as well as from a Chinese Wagner disease family in which a sequence variant was not yet found, revealed two putatively important effects on CSPG2/Versican mRNA isoforms. First, we identified a 39-nt deletion in mRNA isoform V0 in patients carrying c.4004-1G
A (W95-137) and c.4004-5T
A (W95-038) variants. The same 39-nt deletion was found previously in a Japanese family with Wagner disease, carrying a c.4004-2A
G mutation.6 Because the deleted 13 amino acids of V0 does not harbor a GAG attachment site, this truncation is likely to have little or no effect on the function of the protein, although changes on secondary or tertiary structure cannot be ruled out. The c.4004-5T
C variant in the four Dutch families with Wagner disease and the ERVR family and the unknown mutation in the Chinese family were not associated with a clear-cut structural effect on the V0 isoform. However, we assume that all the identified intron 7 variants lead to a decrease of isoforms V0 and V1. Possibly, patients with the c.4004-5T
C variant harbor an additional, as yet unidentified, intron 7 sequence variant that contributes to the skipping of exon 8. The reduction in mRNA isoforms V0 and V1 is difficult to quantify, but is probably reflected in the second effect we observed on CSPG2/Versican mRNA, a significant increase in isoforms V2 (>38-fold) and V3 (>12-fold), which have five to eight or no GAG side chains, respectively.
Given the relatively low V2 and V3 splice isoform abundances in normal blood samples (Table 2) , even a small decrease in the V1 isoform is predicted to have a dramatic effect on the expression of the V2 and V3 isoforms (Table 3) . In view of the unknown absolute levels of CSPG2/Versican isoforms in different compartments of the human eye, it is difficult to predict the functional effect of the observed mRNA isoform changes. Recently, Zhao and Russell reported on the expression of CSPG2/Versican in human trabecular meshwork and ciliary muscle with V1, V0, V3, and V2 in a decreasing order of expression levels,16 which is similar to the expression levels we found in peripheral blood and skin fibroblasts. Of note, as described in the ocular phenotype section, the patients with ERVR clearly had a more severe eye phenotype than did patients with Wagner disease, which suggests that other genetic factors play a role. A more severe retinal phenotype in patients with ERVR is not reflected in the QPCR results of patients with ERVR (W95-124P1 and -P2) versus Wagner disease patients with the same intron 7 variant (W95-131P1 and -P2; Table 4 ), as no consistent difference in upregulation of the V2 and V3 isoforms was observed. We hypothesize that insufficient quantities of isoforms V0 and V1 may underlie both Wagner disease and ERVR. Whether subtle quantitative differences of V0 and/or V1 isoform can explain the difference in clinical severity of Wagner disease and ERVR remains to be investigated. The method used is not accurate enough to quantitate these subtle differences. Possibly, expression levels of the normal CSPG2/Versican gene also play a role.
Natural variation of CSPG2/Versican isoform expression may also explain the relatively high upregulation we identified in the Chinese family. If the normal expression of CSPG2/Versican isoforms V0 and V1 in the Asian population is higher than in the white population, or if isoforms V2 and V3 expression levels are lower in the Asian population than in the white population, disruption of V0 and V1 splicing would lead to higher relative upregulation of V2 and V3, as we observed. We propose that the causal variant resides in intron 7 or 8 causing similar balance shifts as found in the Dutch families with Wagner disease. The efforts to identify such variants are currently ongoing.
The absence of disease-causing mutations in the sizeable ORF (10,200 bp) of CSPG2/Versican in the families included in this study and in at least five other unrelated patients with Wagner disease (Black GCM, unpublished results, 1999), the conspicuous clustering of putative causal variants in the splice acceptor site of intron 7 found in this study and in the previously reported Japanese family,6 and the consistent V2 and V3 mRNA isoform upregulation lead us to believe that Wagner disease and ERVR are caused by an imbalance of the CSPG2/Versican isoforms, mediated by intronic variants. It is estimated that as much as 15% of point mutations causing human diseases affect the canonical splice sites,17 but this may be an underestimate. In most cases, the invariant AG or GT sequences of the 3' and 5' splice sites are mutated. The polypyrimidine stretch of 3' splice sites often harbor benign sequence variants. However, in a family with X-linked mental retardation, an 11C
T change of intron 1 of the ARHGEF6 gene resulted in skipping of the second exon.18 To the best of our knowledge, the only other inherited disorder that is exclusively caused by intronic mutations and exon skipping is autosomal dominant sensorineural deafness type 5.19
In conclusion, in this study the pathologic role of CSPG2/Versican mutations in Wagner disease has been significantly strengthened. We provide strong evidence that, in our cohort of patients, Wagner disease is associated with CSPG2/Versican intron 7 variants which result in a 39-nt truncation of isoform V0 in at least two of seven families. In addition, we show that Wagner disease and ERVR are allelic disorders and that they can be associated with the same intron7exon 8 splice site variant. The intron 7 variants consistently result in a dramatic upregulation of the V2 and V3 isoforms, most likely due to the suppression of the V0 and V1 mRNA isoforms (Fig. 6) . Both a V2/V3 overexpression model or a partial V0/V1 loss-of-function model could underlie Wagner disease and ERVR, though we favor the latter model. The pathogenic mechanism of the disease likely involves a reduction of chondroitin sulfate side chains of CSPG2/Versican. It can be predicted that other Wagner disease mutations can be found deep in intronic regions, most likely introns 7 (14,697 nt) or 8 (3,268 nt). Hence, Wagner disease and ERVR seem to belong to a largely overlooked group of diseases that are caused by a novel disease mechanism based on mRNA isoform balance shifts. Recent studies showed how splicing modulation may restore the function of the cystic fibrosis transmembrane conductance regulator.20 If such an approach becomes feasible in a therapeutic setting, Wagner disease may be amenable to this treatment.
|
| Acknowledgements |
|---|
| Footnotes |
|---|
Submitted for publication February 7, 2006; revised March 16, 2007; accepted May 24, 2006.
Disclosure: A. Mukhopadhyay, None; K. Nikopoulos, None; A. Maugeri, None; A.P.M. de Brouwer, None; C.E. van Nouhuys, None; C.J.F. Boon, None; R. Perveen, None; H.A.A. Zegers, None; D. Wittebol-Post, None; P.R. van den Biesen, None; S.D. van der Velde-Visser, None; H.G. Brunner, None; G.C.M. Black, None; C.B. Hoyng, None; F.P.M. Cremers, None
The publication costs of this article were defrayed in part by page charge payment. This article must therefore be marked "advertisement" in accordance with 18 U.S.C.
1734 solely to indicate this fact.
Corresponding author: Frans P. M. Cremers, Department of Human Genetics, Radboud University Nijmegen Medical Centre, PO Box 9101, 6500 HB Nijmegen, The Netherlands; f.cremers{at}antrg.umcn.nl.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
A. Ohno-Jinno, Z. Isogai, M. Yoneda, K. Kasai, O. Miyaishi, Y. Inoue, T. Kataoka, J.-S. Zhao, H. Li, M. Takeyama, et al. Versican and Fibrillin-1 Form a Major Hyaluronan-Binding Complex in the Ciliary Body Invest. Ophthalmol. Vis. Sci., July 1, 2008; 49(7): 2870 - 2877. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. P Meredith, A. J Richards, D. W Flanagan, J. D Scott, A. V Poulson, and M. P Snead Clinical characterisation and molecular analysis of Wagner syndrome Br. J. Ophthalmol., May 1, 2007; 91(5): 655 - 659. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. E. Keller, M. J. Kelley, and T. S. Acott Extracellular Matrix Gene Alternative Splicing by Trabecular Meshwork Cells in Response to Mechanical Stretching Invest. Ophthalmol. Vis. Sci., March 1, 2007; 48(3): 1164 - 1172. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |