Article Text

PDF

Wolfram syndrome: a clinical and molecular genetic analysis
  1. Philipp Eller,
  2. Bernhard Föger,
  3. Roland Gander,
  4. Teresa Sauper,
  5. Monika Lechleitner,
  6. Gerd Finkenstedt,
  7. Josef R Patsch
  1. Department of Medicine, Anichstrasse 35, A-6020 Innsbruck, Austria
  1. Dr Föger, bernhard.foeger{at}uibk.ac.at

Statistics from Altmetric.com

Editor—Wolfram syndrome (OMIM 222300) is a progressive neurodegenerative disorder characterised by the association of juvenile, non-autoimmune, insulin dependent diabetes mellitus and optic atrophy. The physician D J Wolfram, who reported four cases in 1938, is credited with the first description.1 With the identification of other clinical features, Wolfram syndrome was also referred to as DIDMOAD syndrome (diabetesinsipidus, diabetes mellitus,optic atrophy and deafness). The initial manifestation of Wolfram syndrome is typically, but not invariably, insulin deficient diabetes mellitus at a median age of 6 years, followed by optic atrophy at 11 years.2 In the second decade, many patients develop central diabetes insipidus and sensorineural deafness. Additional, but less frequent neurological and endocrinological abnormalities are atonic bladder, ataxia, myoclonus, peripheral neuropathy, hypogonadism, and a relatively high incidence of depression and psychotic behaviour.3-5 Death occurs between 25 and 49 years (median 30 years).

The prevalence of Wolfram syndrome was estimated at 1/100 000 in a North American population and 1/770 000 in the United Kingdom.2 6 Family studies indicate autosomal recessive inheritance.

Genetic linkage analysis has localised Wolfram syndrome to the short arm of chromosome 4 (4p16),7 8 and the primary genetic defect of the syndrome was attributed to theWFS1 gene, which consists of eight exons encoding a 890 amino acid polypeptide with an apparent molecular weight of 100 kDa.9 10 Hydrophobicity analysis suggested a transmembrane protein comprising three structural domains: a hydrophilic amino terminus separated from a hydrophilic carboxy-terminal tail by a hydrophobic region containing nine predicted transmembrane segments. Homologous genes in Mus musculus, Rattus norvegicus, andDrosophila melanogaster have been described (GenBank AF084482/AJ011971, AF136378; FlyBase FBgn0039003).

The majority of patients with Wolfram syndrome carry loss of function mutations in the WFS1gene.9-11 Defects in mitochondrial DNA were reported in some sporadic cases,12-14 but the mitochondrial genome was not found to be systematically involved in the disease.15-18 Moreover, a recent publication suggested a further nuclear gene locus (chromosome 4q22-24) associated with Wolfram syndrome.19

In the present study, we report three cases of Wolfram syndrome carrying novel mutations of the WFS1 gene and describe mRNA heterogeneity in the 5′ non-coding region of theWFS1 gene, which is caused by alternative splicing.

Methods

Blood samples were collected from available family members after informed consent was obtained. For genotyping, DNA was extracted from peripheral blood lymphocytes using QIAamp DNA Blood kit (Qiagen). The entire WFS1 coding sequence (GenBankAF084481/Y18064), consisting of exons 2-8, was amplified in the index patients with intronic primers, as described previously.9 10 Polymerase chain reaction (PCR) analysis was in a total volume of 20 μl containing 0.5 μmol/l of each forward and reverse primer, 0.5 U HotStarTaqpolymerase (Qiagen), 1.5 mmol/l MgCl, and 200 μmol/l dNTP mix. After an initial activation step at 95°C for 15 minutes, cycling was performed at 94°C for one minute, 57°C for one minute, 72°C for one minute (30 cycles), followed by 10 minutes at 72°C in a T3 Thermocycler (Biometra). Dye modified primers were used in conjunction with a commercial dideoxy sequencing kit (Amersham Pharmacia) to obtain sequence reaction products, which were electrophoresed in 4.3% acrylamide/7 mol/l urea Long Ranger gels (FMC Bioproducts) on a Licor 4000L fluorescent sequencer (MWG, Germany), and analysed by means of BaseImagIR version 4 software (Licor Biotechnologies).

The 1038ins(C) mutation was detected with a labelled PCR fragment on a denaturing polyacrylamide gel. As we could not directly sequence the PCR products in patient 2.1, a fragment of exon 8 containing the suspected mutations was subcloned into pCR2.1 vector (Invitrogen) as described by the manufacturer and at least four colonies sequenced for each allele. The 2315ins(T) mutation was also detected with PCR (primer 8-6, as previously described),10 digestion withBsaOI (Promega), and 2% agarose gel electrophoresis, with three bands in the wild type allele (165 bp, 128 bp, 64 bp), and two bands (229 bp, 128 bp) in the mutant. The 2638-2643del(GACTTC) mutation destroys aTaqI (Roche) cleavage site with an additional band of 212 bp in the mutant allele compared to the wild type (48 bp, 49 bp, 164 bp) site in a fragment of exon 8 (forward primer CCT CAA CTG CAT GGC CCA GC, reverse primer AGG TCT CTG CAG CCA CAG TC). The 1675G→A mutation eliminates aCfoI cleavage site (data not shown).20

Total RNA was extracted from peripheral lymphocytes using Rneasy kit (Qiagen). Poly A+ RNA of human pancreas, brain, heart, and lung was purchased from Clontech. First strand cDNA was synthesised using Superscript II (Gibco BRL) for reverse transcriptase-PCR (RT-PCR). cDNA was subcloned into pCR2.1 vector (Invitrogen) and at least three colonies sequenced for each splice variant. In order to estimate the relative relationship of the two isoforms to each other, we performed a competitive PCR with the fluorescently labelled forward primer (IRD800-CCC GAA CAA CTT TTC TGC CG) and reverse primer (AGG CTG TGG CAT TGA GTC GGG AAC G) (MWG). The PCR product was diluted to a final concentration of 5 pg/μl to guarantee an unsaturated signal-noise ratio. The splice variants of 171 and 167 bp were separated by polyacrylamide gel electrophoresis, and the signal intensity of the respective bands was measured by means of the Quantity One version 4.1.1 software (Biorad).

Statistical analyses were performed with SPSS version 8.0 software. Data are shown as means (SD) for the indicated number of experiments. The ratio of the splice variants in different tissues was compared by the one-way ANOVA analysis.

Patients

Patients were recruited at the local referral centre. Two unrelated families of Austrian descent were studied (table1).

Table 1

Clinical manifestations in patients with Wolfram syndrome

PATIENT 1.1

The patient was a 29 year old woman with unaffected parents, three healthy sibs, and one affected brother (patient 1.2). There was no reliable information in regard to consanguinity in the family. One grandmother had diabetes mellitus type 2 diagnosed at the age of 64.

At the age of 7, the patient developed insulin dependent diabetes mellitus (IDDM) after a measles infection. At the age of 11, bilateral optic atrophy developed with reduced colour vision and reduced visual acuity. There was no evidence of diabetic retinopathy. A visual field test showed central scotomas. At the same age, bilateral high tone sensorineural hearing loss greater than 50 decibels above 4000 Hz emerged. At the age of 20, desmopressin was started for central diabetes insipidus. The bone mineral density of the lumbar spine and proximal femur as measured by dual energy xray absorptiometry was normal. The neurogenic bladder (with overactive, hypocontractile, disinhibited detrusor function and low capacity of 180 ml) was initially treated with several intravesical electrical stimulations. The patient was managed conservatively with clean intermittent self-catheterisation and anticholinergic medication. As vesicoureteral reflux developed at the age of 28 years, an ileocaecal bladder augmentation using the Mainz Pouch I technique was performed to protect the upper urinary tract.21 22 The patient was prescribed a biphasic hormone combination (oestradiol valerate plus norgestrel) to regularise the menstrual cycle.

PATIENT 1.2

The 20 year old man presented with IDDM at the age of 6. Ophthalmological examination showed a progressive loss of visual acuity owing to bilateral optic atrophy (0.6 and 0.8 at the age of 7, 0.1, and 0.3 at the age of 13 for the right and left eye, respectively). Alcohol abuse precipitated several (hypoglycaemic) generalised seizures. Electroencephalographic recordings performed after the seizures were normal. GAD antibodies were normal (0.23 U/ml, reference values <6 U/ml), as determined at the age of 20.

PATIENT 2.1

The patient was a 22 year old woman, who developed IDDM at the age of 8 years. There was no family history of IDDM or consanguinity. Islet cell and glutamic acid decarboxylase (GAD-65) antibodies, as determined at the age of 21, were not detected. There was no macroalbuminuria and neurothesiometry and electromyography were normal. Interestingly, the patient underwent a cataract extraction of the right eye at the age of 2 even before the onset of diabetes mellitus, two operations for strabismus at the ages of 3 and 4, and a Nd:YAG posterior capsulotomy of the optic lens at the age of 18. She began to complain of progressive loss of visual acuity at the age of 16 and four years later bilateral atrophy of the optic nerves was detected. Neurological examination showed only low visual acuity, but cranial magnetic resonance imaging at 21 years showed atrophy of the median brainstem and the optic nerves. Additionally, there was a diffuse area of high signal on PD and T2 weighted images in the ventral pons, as recently described.23 The visual field test showed peripheral constriction, and ophthalmoscopic examination uncovered diabetic non-proliferative retinopathy with pale optic discs. A discrete nuclear cataract of the left myopic eye was present. There was no evidence of diabetes insipidus and the audiogram was normal.

Results

Sequences obtained from the family members were compared to theWFS1 sequence and to normal control chromosomes. Analysis of the complete coding sequence of theWFS1 gene in the index patients showed three novel mutations (table 2).

Table 2

Mutations and polymorphisms in the WFS1 gene

Fig 1 illustrates in index patient 1.1 homozygosity for a 1038ins(C) insertion of 1 bp at position 1038 in exon 8 predicting a frameshift starting at amino acid 346. The next out of frame stop codon at nucleotide positions 1185-1187 was predicted to be read through because of the homozygous polymorphism 1185T→C down to the stop codon of nucleotides 1622-1624. The putative protein would lack the transmembrane domains and the hydrophilic carboxy tail. Patient 1.2 was also homozygous for the 1038ins(C) mutation.

Figure 1

Homozygous 1038ins(C) mutation in patient 1.1. Sequence analysis of the WFS1 gene of patient 1.1. Genomic DNA was amplified, purified, and directly sequenced with fluorescently labelled primers. The region of exon 8 showing the homozygous insertion of 1 bp at position 1038 in patient 1.1 and the corresponding sequence from a normal control (wt) are presented.

As shown in fig 2, patient 2.1 was found to be a compound heterozygote for a (paternal) 2315ins(T) mutation resulting in an Ins772fs/ter776 frameshift and premature truncation of the predicted protein, and a (maternal) 2638-2643del(GACTTC) mutation resulting in an in frame deletion of aspartate and phenylalanine of codons 880-881. Both amino acids were reported to be conserved in the respectiveWFS1 sequences of Mus musculus and Rattus norvegicus(GenBank AF084482/AJ011971, AF136378).

Figure 2

Restriction fragment length polymorphism analysis of exon 8 in family 2. DNA fragments of exon 8 were amplified, digested with BsaOI and TaqI, respectively, and analysed by 2% agarose gel electrophoresis. Analysis of the family showed (A) that patient 2.1 and her father were heterozygous for the 2315ins(T) mutation, and (B) that patient 2.1 and her mother were heterozygous for the 2638-2643del(GACTTC) mutation. Lane M was loaded with DNA size marker pBR322/BsuRI and lane 1 with the undigested PCR product.

Sequence analysis also showed a number of polymorphisms in both families (table 2). Patient 1.1, for instance, was not only homozygous for the 1038ins(C) mutation, but also for the following polymorphisms: 684G→C, 997A→G, 1185T→C, 1500T→C, 1832A→G, 2433A→G, and 2565G→A. The rare coding variant 1675G→A was regarded recently as a risk allele for affective disorders.20 We found no evidence for psychiatric illness (as judged by history, examination, and hospital records) in patient 2.1 or in her unaffected father, both heterozygous for the 1675G→A variant. However, the functional relevance of this coding variant was difficult to clarify, as it was localised to the same allele as the 2315ins(C) mutation.

By screening cDNA of the WFS1 gene, we found mRNA heterogeneity in the 5′ untranslated region of theWFS1 gene. The two mRNA isoforms were isolated from human lymphocytes and designated type 1 and type 2, respectively. They differed from each other only in their Kozak sequence: type 2 was identical to type 1, except for being 4 bp shorter, corresponding to the four bases immediately upstream from the ATG translation initiation codon (fig 3A). Thus, the strong context of the putative ATG start codon of isoform type 1 was changed to a weak one of isoform type 2 with a cytidine instead of the original adenosine in nucleotide position –3 (that is, three nucleotides upstream from the ATG codon).24-26 We directly sequenced both exons 1 and 2 with the adjacent splice sites to exclude any genomic deletions. However, instead of reflecting a genomic mutation, the deletion of 4 bp in isoform type 2 was generated by the use of alternative 3′ splice sites (fig 3B). mRNA isoform type 1 corresponded to the previously published sequence (GenBank AF084481/Y18064), whereas isoform type 2 was a novel splice variant of the WFS1 gene. Each splice variant was found in patients with Wolfram syndrome as well as in all negative controls (n=9). To test the tissue distribution of the splicing variants we performed RT-PCR, competitive PCR with fluorescently labelled primers, and polyacrylamide gel electrophoresis. The ratio of the splice variants differed in brain (type 1/type 2 = 0.92 (0.34)) compared with pancreas (1.25 (0.41)), heart (1.29 (0.65)), and lung (1.10 (0.34)) by ANOVA analysis (means (SD), n=8, p<0.05). The occurrence of both splice variants in pancreas, brain, heart, and lung was confirmed by nested PCR, overdigestion with TurboNaeI (Promega), and 2% agarose gel electrophoresis. Thus, alternative splicing of theWFS1 heterogeneous nuclear RNA appeared to be ubiquitous.

Figure 3

Splice variants of the WFS1 gene. (A) cDNA was subcloned into the pCR2.1 vector and several colonies sequenced for each mRNA isoform. Two representative samples for both isoforms are shown. Isoform type 2 is identical to type 1, except for being 4 bp shorter. The deletion is localised to the Kozak sequence immediately upstream from the functional ATG start codon. (B) 5′ and 3′ splice site of intron 1 of the WFS1 gene. Two tandem repeats of the tetranucleotide sequence containing ag (caggcagg) at the transition from intron 1 to exon 2 lead to two alternative 3′ splice sites, which account for the different splice variants of the WFS1 gene. The consensus sequence for 3′ splice sites is shown. The exonic and intronic sequences are indicated in upper case and lower case letters, respectively. Nucleotide positions are given counting from the first base of the start codon. Nucleotides are written in the single letter code (A = adenosine, G = guanosine, R = A or G (puRine), C = cytidine, T = thymidine, Y = C or T (pYrimidine), N = any Nucleotide).

Discussion

In the work presented here, three novel mutations were identified in the WFS1 gene in patients with Wolfram syndrome. All three patients suffered from the typical association of insulin dependent, non-autoimmune diabetes mellitus and optic atrophy. Additional typical findings were hypoacusis, central diabetes insipidus, hypogonadism, and cataract.10 27 28 The sibs from family 1 both carried a homozygous 1038ins(C) mutation, predicting a severely truncated protein, whereas patient 2.1 appeared to be compound heterozygote for a 2315ins(T) and a 2638-2643del(GACTTC) mutation. All three mutations were localised to exon 8. The two insertions of 1 bp are typical loss of function mutations and confirm the notion that lack of the putative gene product causes the disease.9-11 The in frame deletion of 6 bp in patient 2.1 (del880-881DF), which is very close to the carboxy-terminus ofWFS1 pointed, together with other mutations (del882fs/ter937, del883fs/ter949, P885L), to the functional relevance of the cytoplasmic carboxy-terminus of theWFS1 gene product, which is hypothesised to interact with other, as yet unidentified, proteins.9 11The location and nature of WFS1 mutations might thus help in understanding the molecular basis of Wolfram syndrome and the role of the WFS1 protein in protecting neurones and islet β cells from premature cell death.29

Additionally, we present evidence for mRNA heterogeneity of the 5′ non-coding region of the WFS1 gene. mRNA heterogeneity is the result of alternative splicing of intron 1, which contains a redundant 3′ splice site. Competition of the tandemly repeated tetranucleotide CAGG (12-20 bp downstream from the putative branchpoint, fig 3) for the small nuclear ribonucleoprotein U2AF35 is in line with previous reports. Tandemly repeated AGs competed actively when placed 13-22 nucleotides downstream of the branchpoint in a yeast intron.30-32 When the branch to 3′ splice site distance was significantly longer (at least 35 nucleotides), a linear search mechanism in the 5′→3′ direction showed strong preference for the 5′ most AG.33 The two splice variants of the WFS1 gene differed from each other only in their Kozak consensus sequence, where a pyrimidine was substituted for the highly conserved purine at position –3. Neither the upstream leader sequence nor the open reading frame were affected. Thus, the final WFS1 protein of the two splice variants would be expected to be identical. Moreover, no major differences in the relative density of the two splice variants was evident in the tissues examined. To address the functional relevance, if any, of alternative splicing of WFS1, further molecular genetic and functional analysis will be required. Taken together, our data extend the mutational spectrum and provide, to our knowledge, the first evidence for natural variation in the translation regulatory region of the WFS1 gene.

Acknowledgments

We thank our patients for their participation in this study. This work was supported by a grant from the “Jubiläumsfond der Österreichischen Nationalbank” (grant 8539 to BF) and a grant from the FWF (grant 11693-MED to JP).

References

View Abstract

Request permissions

If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.