Abstract
In this manuscript, we describe a case of type 2A von Willebrand disease (VWD) caused by the novel heterozygous G>A transition at nucleotide 3538, which should result in the putative, nonconservative substitution of G1180R. This mutation was reproduced by site-directed mutagenesis; however, the recombinant mutant protein was efficiently secreted from cells and assembled correctly into multimers. Because the substitution is located at the last nucleotide of exon 26, the patient's platelet von Willebrand factor (VWF) mRNA was analyzed and 3 transcripts were observed: the normal transcript without the 3538G>A transition, a transcript with the in-frame deletion of exon 26, and a transcript with the in-frame deletions of exons 23 and 26. These deletion VWF cDNA constructs were created and the resulting recombinant proteins were analyzed following transfection into COS-7 cells. Cotransfection results demonstrate that the exon-skipped transcripts led to intracellular retention, and the levels of VWF antigen (VWF:Ag) produced by these constructs were as follows: del23/26<del26< G1180R≤wild type. The homozygous exon-skipped transcripts show the presence of only the lowest molecular weight multimers. The G>A transition at nucleotide (nt) 3538 does not result in the expression of the G1180R missense mutation, but rather leads to exon skipping, which is the pathogenic basis of the patient's phenotype. This is the first report of a coding region mutation resulting in the skipping of 2 nonadjacent exons. (Blood. 2004;104:2739-2745)
Introduction
Von Willebrand disease (VWD) is the most common inherited bleeding disorder known in humans, with published prevalence rates between 0.01% and 1%.1,2 The clinical spectrum of the disease is broad, ranging from mild mucocutaneous bleeding observed in type 1 VWD (quantitative deficiency of von Willebrand factor [VWF]) to severe, spontaneous bleeding in type 3 VWD (complete absence of VWF). Type 2 VWD is characterized by the presence of qualitatively abnormal VWF and is a heterogeneous condition both phenotypically and genotypically. There are 4 distinct subtypes recognized.3 Type 2A and type 2B VWD are characterized by the loss of high-molecular-weight VWF multimers; however, these subtypes are best distinguished based on VWF interactions with platelets. This interaction is increased in type 2B VWD as is observed by the enhanced ristocetin-induced platelet agglutination, and significantly reduced in type 2A disease. Type 2M VWD refers to qualitatively abnormal variants with decreased platelet-dependent function not associated with the loss of high-molecular-weight multimers. Type 2N VWD demonstrates defective binding of VWF to factor VIII.
The accurate discrimination of subtypes is important, particularly when it affects patient management. This has become easier since our understanding of the condition has progressed. However, challenging cases remain, particularly when the diagnosis is based on phenotypic data alone. Coagulation-based laboratory tests vary over time within individuals, and an accurate diagnosis is sometimes difficult to achieve, even when repeated tests are performed. The VWF gene was cloned in 1985,4-7 enabling the identification of the molecular genetic basis in some of the subtypes. The ability to combine phenotypic and genotypic information has helped to improve the accuracy of diagnosis.
The inheritance pattern of type 2A VWD is usually autosomal dominant, and significant progress has been made in defining the molecular genetic pathology of this subtype. Missense mutations account for the majority of Type 2A cases, and most are located within the region of exon 28 that encodes the A2 domain. Type 2A mutations have also been described at the C-terminus of the protein and in the D2 and A1 domains. There have been 2 distinct pathogenic subgroups described within type 2A VWD8,9 ; group 1 is caused by defects in the synthesis and intracellular transport of multimers, and group 2 results from an increase in plasma proteolysis of VWF by the ADAMTS13 protease following secretion into blood. Type 2A mutations can result in defective dimerization or multimerization, and splicing defects found in the D3 domain leading to type 2A have previously been described in abstract form.10,11 VWF undergoes extensive posttranslational modifications by initially dimerizing through multiple intermolecular disulfide bonds between the carboxyl-terminal ends of the protein and then, once transported to the Golgi, forming interdimer disulfide bonds using cysteines found in the D3 domain. The resulting multimers range in size from 500 to 20 000 kDa, and it is the highest molecular weight multimers that are most functionally active in terms of VWF's adhesive properties.
VWF is the largest soluble plasma protein. The VWF gene encompasses 178 kb of genomic DNA and is comprised of 52 exons. Based on the number of exons that have to be spliced together, it is not surprising that several splicing mutations resulting in various types of VWD have previously been described (http://www.shef.ac.uk/vwf/).12 Splicing defects are most often the result of a mutation in the 5′ (donor) or 3′ (acceptor) splice sites located at the exon/intron boundary, however, changes in the last nucleotide of an exon have also been shown to cause exon skipping in a number of conditions.13-20 In this study we describe a patient with a novel VWF mutation, 3538G>A, which results in the in-frame exon skipping of exon 26, and exons 23 and 26. To our knowledge, this is the first description of skipping of 2 nonadjacent exons as a pathogenic mechanism in VWD or any other condition.
Patients, materials, and methods
Patients
The study family was submitted to the Canadian Type 1 VWD Study as a potential type 1 VWD family. Each individual was informed of the experimental nature of the study and gave his/her informed consent. The study was approved by the institutional review board at Queen's University, and at the source institution. To be eligible for the study, at least one member from the family had to have a personal history of excessive mucocutaneous bleeding, von Willebrand factor antigen (VWF:Ag), and ristocetin cofactor activity (VWF:RCo) between 0.05 and 0.50 U/mL measured on at least 2 different samples. Whole blood samples for DNA extraction were collected in 3.2% sodium citrate (at a ratio of 9:1 vol/vol) from the index case and available immediate family members, and a plasma sample was collected from the index case for repeat VWD phenotypic studies to be performed in a central laboratory.
Coagulation studies
Laboratory tests for VWF:Ag, VWF:RCo, and factor VIII coagulant activity (FVIII:C) were performed at the source clinic attended by the patient and were repeated on different plasma samples at the Clinical Hemostasis Laboratory at Kingston General Hospital (Canada). All available laboratory results were averaged. In the source clinic laboratories, VWF:Ag and VWF:RCo values were assayed using the locally validated methodology. In the central laboratory, VWF:Ag was measured using the IMMUBIND VWF enzyme-linked immunosorbent assay (ELISA) kit according to the manufacturer's instructions (American Diagnostica, Greenwich, CT). The VWF:RCo was measured by platelet aggregometry using freshly prepared, washed normal platelets, and the FVIII:C was measured using a one-stage assay. All measurements were made against a non-ABO-matched commercial reference plasma that was calibrated against the 91-666 or 97-586 World Health Organization (WHO) Plasma Standard. Ristocetin-induced platelet aggregation (RIPA) was performed on a Platelet Aggregometer (Chrono-log, Havertown, PA). VWF multimers were analyzed by electrophoresis using a 1.6% sodium dodecyl sulfate (SDS) agarose gel followed by electrotransfer to a nylon membrane. The multimers were visualized using the chemiluminescent visualization kit from Amersham Pharmacia Biotech (Baie D'Urfe, QC).
PCR amplification of genomic DNA and mutation identification
A blood sample was collected from the patient and each family member, and genomic DNA was isolated using a salt extraction method.21 DNA corresponding to the coding sequence of exons 18 to 52 of the VWF gene (encoding the mature VWF subunit) and at least 20 nucleotides of the adjacent intron were amplified by polymerase chain reaction (PCR). Primer sequences have been previously published.22,23 Additionally, genomic DNA was PCR amplified for introns 23, 25, and 26 and 300 bp of the 3′ end of intron 22 (primer sequences, Table 1). Using a DNA thermal cycler (Perkin Elmer Life Sciences, Shelton, CT), DNA was amplified for 35 cycles of 45 to 60 seconds at 94°C, 45 to 60 seconds at 53°C to 60°C, and 45 to 60 seconds at 72° C, and the amplified products were sequenced directly on an ABI model 373 automated sequencer (Cortec DNA Service Laboratories, Kingston, ON). All DNA sequences were compared with normal VWF DNA sequences with the assistance of Vector NTI Suite software (Infor-MAX, Bethesda, MD).
Primer . | Sequence (5′-3′) . | Location . |
---|---|---|
RT-PCR | ||
16F | CAGTGCCCCTGTTACTATG | 13/270 |
21R | CTCCCCGTCAAACAGCTCAAT | 18/319 |
20F | CCTTCGACGGGCTCAAATAC | 17/340 |
28R | GCGGTCGATCTTGCTGA | 24/653 |
28F | GCTCCCACGCCTACATC | 24/502 |
31R | AGTATCGCACAGCAAAGCC | 25/711 |
30F | GTGTCAGTGCTGCAGTATG | 25/275 |
36R | AAAGTCTTGGGGCTGAGCT | 29/223 |
35F | ACTCGGCACATCGTGACCT | 28/290 |
41R | AGGTGAAGTTGGGTCTGCA | 31/2594 |
40F | TGTGGCCTGTGTGAAGTAG | 31/616 |
52R | GAACATGCAGAGGACTGGC | 38/1052 |
Sequencing | ||
22F | ATGCTTGTTGGGGACAGT | int22nt(−343) |
22R | AGGTTGCTGCTGGTGAGGTCAT | 20/398 |
23F | TTCCCCTTTCCCCTGAGCCG | 20/259 |
23R | ACTCTGTGTCCATACCACCA | 20/865 |
25F | TGACCTGGAGGACGGCCACATT | 21/398 |
25R | TGGCTCAGGGTGCTGACACGTG | 22/292 |
26F | GGTGACGCCCATAGTCCCC | 22/61 |
26R | CCGGCCAGCCACCTCACACACT | 23/199 |
Mutagenesis | ||
mut26G/R-D | ACTGCCCTCCAAGGAAAATCCT | 22/333 |
mut26G/R-U | AGGATTTTCCTTGGAGGGCAGT | 23/135 |
Primer . | Sequence (5′-3′) . | Location . |
---|---|---|
RT-PCR | ||
16F | CAGTGCCCCTGTTACTATG | 13/270 |
21R | CTCCCCGTCAAACAGCTCAAT | 18/319 |
20F | CCTTCGACGGGCTCAAATAC | 17/340 |
28R | GCGGTCGATCTTGCTGA | 24/653 |
28F | GCTCCCACGCCTACATC | 24/502 |
31R | AGTATCGCACAGCAAAGCC | 25/711 |
30F | GTGTCAGTGCTGCAGTATG | 25/275 |
36R | AAAGTCTTGGGGCTGAGCT | 29/223 |
35F | ACTCGGCACATCGTGACCT | 28/290 |
41R | AGGTGAAGTTGGGTCTGCA | 31/2594 |
40F | TGTGGCCTGTGTGAAGTAG | 31/616 |
52R | GAACATGCAGAGGACTGGC | 38/1052 |
Sequencing | ||
22F | ATGCTTGTTGGGGACAGT | int22nt(−343) |
22R | AGGTTGCTGCTGGTGAGGTCAT | 20/398 |
23F | TTCCCCTTTCCCCTGAGCCG | 20/259 |
23R | ACTCTGTGTCCATACCACCA | 20/865 |
25F | TGACCTGGAGGACGGCCACATT | 21/398 |
25R | TGGCTCAGGGTGCTGACACGTG | 22/292 |
26F | GGTGACGCCCATAGTCCCC | 22/61 |
26R | CCGGCCAGCCACCTCACACACT | 23/199 |
Mutagenesis | ||
mut26G/R-D | ACTGCCCTCCAAGGAAAATCCT | 22/333 |
mut26G/R-U | AGGATTTTCCTTGGAGGGCAGT | 23/135 |
F indicates forward primer; R, reverse primer.
Primer sequences that are italicized and underlined were used to perform reverse-transcriptase reactions, underlined nucleotides are those that are different from the pseudogene, and bold nucleotides indicate base substitutions for mutagenesis reaction. Location indicates first nucleotide of primer with the numbering according to Mancuso et al,23 except for the location of 22F, which is shown as the distance from the start of exon 23.
Platelet RNA isolation and amplification
A blood sample was collected in EDTA (ethylenediaminetetraacetic acid) from the index case for VWF mRNA analysis. Platelet-rich plasma was separated by centrifugation at 800 rpm for 5 minutes at room temperature. Platelets were then pelleted by centrifugation at 1400 rpm for 10 minutes at room temperature. The platelet pellet was washed twice with phosphate-buffered saline (PBS) and the total RNA was extracted using the QIAamp RNA Blood Mini Kit (QIAGEN, Mississauga, ON). Reverse transcription of the RNA was performed according to the procedure provided by the manufacturer (Expand Reverse Transcriptase Kit; Roche, Laval, QC), and exons 20 to 28 (primer sequences, Table 1) of the VWF cDNA were PCR amplified (a region spanning the putative mutation site). The amplified product was electrophoresed through a 1% agarose gel. The bands of interest were excised from the agarose gel and purified using the QIAquick Gel Extraction Kit (QIAGEN), and the cDNA was sequenced directly on an ABI model 373 automated sequencer (University Core DNA & Protein Services DNA Sequencing Laboratory, University of Calgary, Calgary, AB). To identify rare transcripts, the patient's cDNA was ligated into the pCR 2.1 vector (Invitrogen, Carlsbad, CA), and 48 of the resulting colonies were analyzed by direct sequencing. The patient's VWF cDNA was also PCR amplified and sequenced (Mobix Lab, McMaster University, Hamilton, ON) in the following segments: exons 16 to 21, exons 28 to 31, exons 30 to 36, exons 35 to 41, and exons 40 to 52. A normal control platelet mRNA sample was treated in the same fashion.
Expression studies
The VWF expression vector pCIneoVWFES (kindly provided by Dr P. Kroner, Medical College of Wisconsin) consists of the full-length VWF cDNA inserted into the expression vector pCIneo (Invitrogen). A unique SalI restriction site was created by Dr Kroner in exon 27 at codon 3608 of the VWF gene. This SalI site in conjunction with a XbaI site at codon 357 was used to excise the 5′ end of the VWF gene, which was subsequently subcloned into the multicloning site of Bluescript to produce the vector BluescriptVWFXbaI-SaII. To introduce the G>A transition at nucleotide 3538, site-directed mutagenesis was performed on BluescriptVWFXbaI-SaII using the overlap-extension procedure of Ho et al24 (primers listed in Table 1), and the resulting fragment was subcloned into pCR 2.1 (Invitrogen). To create the constructs lacking exon 26 or exons 23/26, the patient's cDNA was amplified from exon 20 to exon 28 using the primers VWF-HS-ex20nt85-D and SalSDMex28-U, and the amplified products were subcloned into pCR 2.1; all pCR 2.1 inserts were sequenced to verify the integrity of the amplified VWF gene. The fragments were excised from pCR 2.1 with BamHI and SalI (located in VWF cDNA at 2717 and 3608, respectively) and ligated back into BluescriptVWFXbaI-SalI to produce the vectors BluescriptVWFG1180R, BluescriptVWFdel26, and BluescriptVWFdel23/26. Finally, the XbaI to SalI fragments were ligated back into the expression vector pCIneoVWFES to generate the cDNA expression vectors pCIneoG1180R, pCIneodel26, and pCIneodel23/26. All of the final constructs were sequenced to verify the integrity of the products.
COS-7 cells (catalog no. CRL-1651; American Type Culture Collection, Manassas, VA) were cultured in Dulbecco modified Eagle medium (DMEM) containing 2 mM l-glutamine, 100 U/mL penicillin, 100 μg/mL streptomycin, and 10% (vol/vol) fetal bovine serum (FBS) at 37°C and 5% CO2. Cells in the log phase of growth were used to seed tissue culture dishes such that the cells were approximately 50% confluent the following day. These cells were transfected using the calcium phosphate technique with a total of 20 μg DNA comprising 3.2 μg β galactosidase (βgal) reporter construct, 11.8 μg calf thymus DNA, and 5 μg of the mutant plasmid or 5 μg of the wild-type plasmid (pCIneoVWFES) to yield homozygous mutant or wild-type genotypes, respectively. The mutant and wild-type plasmids were also transfected in various ratios to yield a heterozygous genotype (wild-type-mutant ratios, 75:25, 50:50, 25:75). At 48 hours after transfection, media were collected and cells were lysed using a cell lysis solution supplied in the Galacto Light Plus reporter gene assay kit from Tropix (Bedford, MA). The transfection efficiency was determined for each experiment by measuring the βgal activity from the βgal reporter transcript (that was cotransfected with the pCIneo plasmid) using a Berthold Lumat LB 9501 luminometer (Fisher Scientific, Springfield, NJ) and the Galacto Light Plus reporter gene assay. The quantity of recombinant VWF:Ag (rVWF:Ag) present in the media and cell lysates was determined by ELISA using a Bio-Rad plate reader model 3550 (Mississauga, ON) and a polyclonal goat anti-human VWF antibody (Affinity Biologicals, Hamilton, ON), and the quantity of recombinant VWF protein measured in each experiment was normalized according to its transfection efficiency. The standard curve for the ELISA was generated using CryoCheck normal human reference plasma (lot no. 7070) from Precision Biologicals (Dart-mouth, NS). Recombinant VWF:Ag from the media was concentrated by dialysis for subsequent multimer analysis.
Results
Phenotype
The index patient is a 16-year-old female from the Philippines with a long-standing history of excessive mucocutaneous bleeding. She has had recurrent, significant epistaxis since a young age resulting in hospital admission on one occasion and many emergency room visits. The nosebleeds are often accompanied by welling of blood from her eyes. Shortly after menarche, she was admitted to hospital with severe menorrhagia and was treated with blood transfusions, DDAVP (1-desamino-8-D-arginine vasopressin), and cryoprecipitate. Her father has a history of mild mucocutaneous bleeding, although never requiring intervention. No other immediate family members were affected. An extended family history was not available.
The index patient's bleeding time was 8 minutes (normal range, 2.5-9.5 minutes), the partial thromboplastin time was prolonged at 43.2 seconds (normal range, 27.2-36.0 seconds), and the international normalized ratio was normal at 1.0. The mean VWF:Ag, VWF:RCo, and FVIII:C levels are 0.13 U/mL, 0.12 U/mL, and 0.15 U/mL, respectively, and the RIPA showed no response at 0.5 g/L ristocetin. There was also a loss of the high-molecular-weight VWF plasma multimers (Figure 1). There is no discrepancy between the VWF:Ag and the VWF:RCo in the pre-DDAVP results, which may reflect a lack of sensitivity for these assays (especially the VWF:RCo) at lower values. However, 30 minutes after DDAVP her plasma levels were as follows: VWF:Ag, 1.02 U/mL; VWF:RCo, 0.59 U/mL, and FVIII:C, 1.54 U/mL. It would have been interesting to evaluate the VWF multimers after DDAVP to see if they were consistent with a type 2A phenotype, however, this was not done. Based on this information, we reclassified the patient with type 2A VWD.
Sequence variation
Sequence analysis in the index case revealed a putative mutation at 3538G>A (Figure 2). Aside from common polymorphisms, no other sequence variations were observed in exons 18 to 52, introns 22, 23, 25, and 26, or in the flanking intronic sequence of any of the other introns. This change was also observed in the affected father, but was not present in the unaffected mother or brother. This change is located within the D3 domain of the mature VWF subunit, an area critical for proper multimerization.
RNA analysis for the 3538G> A sequence variation
The 3538G>A sequence variation is located at the last nucleotide of exon 26 and to examine its effect, RNA isolated from the platelets of the index case and a healthy control was reverse transcribed and PCR amplified. When exons 20 to 28 were amplified, 3 VWF cDNA products (984, 825, and 684 bp) were observed in the patient, but only 1 (984 bp) was identified in the healthy control (Figure 3). The 984-bp band observed in the healthy control and the patient revealed normal sequence; the 3538G>A transition was not observed in the patient's cDNA (Figure 4A). Sequence analysis of the 825-bp band from the patient revealed in-frame skipping of exon 26, plus a missense mutation P1127R caused by the last nucleotide of exon 25 creating a codon with the first 2 nucleotides of exon 27, instead of exon 26 (CCC to CGG) (Figure 4B-C). Sequence analysis of the patient-specific 684-bp band again showed in-frame skipping of exon 26 with P1127R, plus in-frame skipping of exon 23 with no additional amino acid changes (Figure 4D-E). There were no differences noted between the patient and healthy control for the reverse-transcription (RT)-PCR products derived from exons 16 to 21, 28 to 31, 30 to 36, 35 to 41, or 40 to 52 (Figure 3B-D). Therefore, the patient is making 3 VWF transcripts: 1 normal transcript, 1 with skipping of exon 26 plus P1127R, and 1 with skipping of both exons 23 and 26 plus P1127R. Because the 984-bp band does not appear to contain the 3538G>A transition, and the only other products identified are the exon-skipped products, it appears that the allele with the transition mutation does not produce appreciable amounts of the G1180R product; instead, the transition results primarily in exon skipping. Although we observed homozygosity for G at nucleotide position 3538, it is possible that there is a small proportion of transcripts containing the mutant A, undetectable by our sequencing methods. Approximately 50% of the cDNA products that were inserted into pCR 2.1 were missing exons 23 and 26, 33% were wild-type, and 17% were missing exon 26.
Expression and characterization of recombinant VWF
To determine the effect of exon skipping on VWF structure and function, the patient's mRNA was amplified and the 3 spliced products (1 normal and 2 exon-skipped transcripts) were introduced into the pCIneoVWFES expression vector containing the full-length VWF cDNA. The G1180R mutation was also examined to determine if this amino acid substitution could impair VWF synthesis and multimerization. The expression vectors pCIneoVWFES (wild type), pCIneoG1180R, pCIneodel26 (exon 26 skipped), and pCIneodel23/26 (exons 23 and 26 skipped) were transiently transfected into COS-7 cells alone or together at various ratios (wild type/mutant: 100:0, 75:25, 50:50, 25:75, 0:100) to generate wild-type, heterozygous, and homozygous mutant genotypes. To determine whether the mutant recombinant protein was retained within the cell or was efficiently secreted, VWF:Ag levels were assayed in cell lysates and in the conditioned media using an ELISA, and the data are presented in Figure 5. Secretion of rVWF protein from the heterozygous and homozygous G1180R transfections was marginally decreased by 14% and 19%, respectively, when compared with the secretion of the wild-type protein (not statistically significant, P = .36 for media, P = .10 for lysate). However, the exon-skipped transcripts resulted in a significant reduction of secreted rVWF protein. Secretion of protein missing exon 26 was reduced by 66% and 93% for heterozygous and homozygous mutant transfections, respectively. Secretion of rVWF protein missing exons 23 and 26 was even more disrupted, with secreted protein levels being reduced by 78% and essentially 100% for the heterozygous and homozygous mutant transfections, respectively. For all of the experiments, the decrease in recombinant protein found in the media corresponded to an increase in intracellular retention. These results illustrate the dominant-negative nature of these aberrantly spliced proteins. Furthermore, these results demonstrate that the mutant cDNA with both exons 23 and 26 spliced out causes more retention of the VWF protein compared with the construct missing only exon 26.
To determine the effect of the mutations on VWF multimers, conditioned media from COS-7 cells for the homozygous wild-type, homozygous mutant, and 50:50 heterozygous transfections were concentrated by dialysis and analyzed by SDS-agarose gel electrophoresis (Figure 6). All sizes of multimers were observed for the G1180R recombinant protein, suggesting that this mutation does not impair multimerization. A full range of multimers was also observed in the media from the heterozygous transfections (involving wild-type and exon-skipped cDNA constructs). Only the lowest molecular weight multimers are observed from recombinant protein missing exon 26 (homozygous) and the rVWF missing both exons 23/26 (homozygous), indicating that the missing cysteine residues are imperative in the formation of multimers.
Discussion
In this report, we present a novel single base-pair change at the last nucleotide of exon 26 that acts as the molecular basis for a young woman's type 2A VWD. We have shown that this mutation results in skipping of not only exon 26, but also of exon 23 in a population of mutant transcripts. To the best of our knowledge, this is the first demonstration of multiple, nonadjacent exon skipping as the pathogenic basis for an inherited disease.
Splicing is a complicated process that involves the precise identification and joining of exons and removal of introns. For this to happen, 5 small nuclear (sn) RNAs (U1, U2, and U4-6) and more than 60 polypeptides must correctly assemble into a large multicomponent ribonucleoprotein complex known as the spliceosome.25,26 The spliceosome recognizes conserved consensus sequences including the 5′ (donor site) and 3′ (acceptor site) splice sites that flank mammalian exons. In addition to the donor and acceptor sites, the last 2 nucleotides of an exon are also semiconserved; the second to last is an A 64% of the time, and the last is a G 73% of the time. It is conceivable that changes to the last nucleotides of an exon might lead to aberrant splicing due to the donor site not being recognized by the spliceosome complex. The mechanism by which the sequential pairing of splice sites is accomplished is not known, however the assembly of the spliceosome, while the RNA transcript is still being transcribed, is presumed to play a major role in the orderly pairing of appropriate splice sites.27,28 It is known that the splicing machinery associates with exons in a 5′ to 3′ order as they are synthesized, however, the covalent joining of exons does not necessarily proceed in a 5′ to 3′ direction due to inherent differences in rates of splicing catalysis between each exon pair.29
In order to understand the transcript that is missing exon 26, we must examine the possibility that a change in the last nucleotide of an exon could lead to aberrant splicing. There are many precedents in the literature for this, all based on the idea that mutations at this site would lead to inactivation of the 5′ donor splice site, and pairing of the donor site from one exon with the acceptor site of another, resulting in skipping of the adjacent exon. This has been described in a number of inherited diseases.13-20 The change in the last nucleotide of exon 26 is the only putative mutation identified in our patient, after sequencing of the coding regions of exons 18 to 52 with flanking introns, and the vast majority of the introns surrounding the skipped exons (intron 22 was not completely sequenced because of its large size and the very low possibility that changes so far from the mutant exon might lead to exon skipping). In the patient's cDNA, we observed only the G at nucleotide position 3538, although it is possible that there are small, unappreciable amounts of the mutant A allele present, which would produce the missense mutation G1180R. However, the expression studies on G1180R clearly show that it does not produce a type 2A VWD phenotype and is not likely contributing to the patient's phenotype. This phenomenon, of a putative missense mutation leading instead to a splicing mutation, has recently been described in a case of severe hypofibrinogenemia.30 Therefore, in our patient, in order to conclude that the single nucleotide change in the last nucleotide of exon 26 is the molecular basis of this patient's type 2A VWD, we must consider how this change leads not only to skipping of exon 26, but also to skipping of exon 23 in some transcripts.
There is strong evidence to support the theory that order and speed of splicing are important determinants of mRNA transcripts with exon skipping.31-34 Despite the fact that the spliceosome associates with the transcript as it emerges from the polymerase, the joining of exons does not occur in order from the 5′ end to the 3′ end of the transcript. Introns with strong splicing consensus sequences are often removed more efficiently than introns with weaker consensus sequences; thus, some exons located downstream will be joined together before slower splicing, upstream exons. Takahara et al31 describe a COL5A1 acceptor site mutation causing Ehlers-Danlos syndrome type 1, which results in the production of a number of transcripts, including one missing exons 5 and 6, one missing only exon 5, and other products formed by the activation of cryptic acceptor sites within exon 5. They, as well as other investigators,32,34 conclude that the order and speed of splicing can be responsible for different transcripts. Splicing mutations of fast-spliced introns often lead to exon skipping, however, splicing mutations of slower spliced introns often lead to cryptic splice site activation and intron inclusion.32,34 In this patient, the substitution is located in the last nucleotide of exon 26, leading to the pairing of the donor site from intron 25 with the 3′ acceptor site from intron 26. Because the donor and acceptor site of intron 26 have a strong consensus splicing sequence score, and no products that could result from the activation of a cryptic splice site were observed, we hypothesize that intron 26 is a fast-spliced intron. Possible cryptic donor sites were examined with the program NNSplice (http://www.fruitfly.org/seq_tools/splice.html), however, only the natural donor site had a high score.
In the Takahara et al paper,31 a transcript is described that is missing 2 adjacent exons; this is explainable by speed and order of splicing. However, it is difficult to apply this argument to the splicing of 2 nonadjacent exons, as in our patient's transcript with both exons 23 and 26 skipped; thus we must consider other potential mechanisms for this aberrant splicing. It is possible that the initial event of skipping exon 26 disrupts the subsequent upstream splicing of exon 23. Exonic sequences have been recently recognized as important controllers of alternative splicing. Exonic splicing enhancers (ESEs) act as binding sites for specific serine/arginine-rich (SR) proteins.35 These SR proteins recognize RNA motifs and act to enhance splicing by either directly recruiting the splicing machinery or by antagonizing the action of nearby silencer elements.36,37 Exonic splicing silencers (ESSs) are not as well understood but seem to work by interacting with negative regulators, often belonging to the heterogeneous nuclear ribonucleoprotein (hnRNP) family.38 In our patient, the change at the last nucleotide of exon 26, leading to the skipping of that exon, results in the destruction of an ESE for SRp40 (analysis done using the web-based program ESEfinder39 http://exon.cshl.edu/ESE/). It is possible that the correct splicing of exon 23 is dependent on that ESE and that its destruction leads to exon 23 skipping in some transcripts.
A second proposed model for exon 23 skipping is the interaction of G-rich and C-rich elements found in introns 22 and 23 that might be facilitated by the exon 26 mutation.40 This interaction could result in the formation of a hairpin structure that would cause looping out of exon 23 and result in the juxtaposition of the donor splice site from intron 22 with the acceptor splice site of intron 23. The secondary structure of VWF mRNA was examined using the program Mfold40 (http://www.bioinfo.rpi.edu/applications/mfold/old/rna/form1.cgi), and sequences from introns 22 and 23 were found to be complementary, making this proposal plausible.
The presence of 2 mutant transcripts and 1 wild-type transcript in our patient is key to understanding the dominant-negative effect of the mutant subunit on the correct multimerization of VWF. Skipping of exon 26 causes the in-frame deletion of 54 amino acids including 9 cysteine residues. It also causes the nonconservative amino acid change P1127R; however, the skipping of 2 critical exons likely contributes more significantly to the pathogenesis of the phenotype in this case. Skipping of exon 23 causes an in-frame deletion of 48 amino acids with the further loss of 3 cysteines. Exons 23 and 26 both encode regions of the D3 domain at the N-terminus of the mature VWF subunit, which plays a critical role in the multimerization of the VWF protein. In light of the location of the skipped exons, deletion of these 12 cysteines would be expected to have disastrous effects on the multimerization of the patient's VWF protein. The expression studies clearly show that the cDNA missing both exons 23 and 26 results in more intracellular retention than the exon 26-skipped construct, however both mutant constructs exhibit a dominant-negative effect on the VWF subunits translated from the normal cDNA, offering an explanation for the severity of the patient's phenotype. The probability that a dimer would contain a mutant subunit is 75%, assuming dimerization in the ER is random. Most of these mutant dimers will be retained within the cell, but some dimers containing a mutant subunit may be processed into multimers. One possible explanation for how this might lead to abnormal multimerization is that these mutant subunits may act as reducing ends during multimer formation; thus additional dimers can be added only to one end of the multimer not both. This would result in a reduction in the size of secreted VWF multimers, which is consistent with the patient's plasma multimer profile. An unexpected finding in our study was that all sizes of multimers were observed in the heterozygous transfection experiments using wild-type and exon-skipped recombinant cDNAs. This discrepancy between the multimer pattern of the heterozygous rVWF proteins and the patient's plasma multimers may well reflect the consequence of several limitations of the recombinant expression system, and is a phenomenon that has been documented in previous similar studies.41
Examination of the Intron Sequence Information System (ISIS)42 revealed 109 genes that demonstrated in-frame exon skipping, in one or more exons; however, none of these genes demonstrated skipping of nonadjacent exons. This is the first report that we are aware of in which a change in the last nucleotide of an exon is shown to cause skipping of 2 nonadjacent exons, and the molecular mechanism demonstrated in this case of type 2A VWD raises interesting questions about the mechanics of splicing in the VWF gene. Furthermore, it highlights the importance of examining putative mutations at the mRNA level.
Prepublished online as Blood First Edition Paper, June 29, 2004; DOI 10.1182/blood-2003-12-4286.
Supported in part by an operating grant from the Canadian Institutes of Health Research (MOP-42467). P.D.J. holds the Aventis-Behring-CHS-AHCDC Fellowship in Hemophilia. L.A.O. is supported by a Canadian Blood Services Graduate Fellowship. D.L. holds a Canada Research Chair in Molecular Hemostasis and is a Career Investigator of the Heart and Stroke Foundation of Ontario.
P.D.J. and L.A.O. contributed equally to this manuscript.
The publication costs of this article were defrayed in part by page charge payment. Therefore, and solely to indicate this fact, this article is hereby marked “advertisement” in accordance with 18 U.S.C. section 1734.
The authors wish to thank Morna Brown and the staff of the Molecular Hematology Laboratory at the Foothills Hospital, Calgary, for their invaluable contributions. We would also like to acknowledge the contribution of the Association of Hemophilia Clinic Directors of Canada to the Canadian Type 1 VWD Study.
This feature is available to Subscribers Only
Sign In or Create an Account Close Modal