- Research article
- Open Access
- Open Peer Review
Functional analysis of splicing mutations in exon 7 of NF1gene
BMC Medical Geneticsvolume 8, Article number: 4 (2007)
Neurofibromatosis type 1 is one of the most common autosomal dominant disorders, affecting about 1:3,500 individuals. NF1 exon 7 displays weakly defined exon-intron boundaries, and is particularly prone to missplicing.
In this study we investigated the expression of exon 7 transcripts using bioinformatic identification of splicing regulatory sequences, and functional minigene analysis of four sequence changes [c.910C>T (R304X), c.945G>A/c.946C>A (Q315Q/L316M), c.1005T>C (N335N)] identified in exon 7 of three different NF1 patients.
Our results detected the presence of three exonic splicing enhancers (ESEs) and one putative exonic splicing silencer (ESS) element. The wild type minigene assay resulted in three alternative isoforms, including a transcript lacking NF1 exon 7 (NF1ΔE7). Both the wild type and the mutated constructs shared NF1ΔE7 in addition to the complete messenger, but displayed a different ratio between the two transcripts. In the presence of R304X and Q315Q/L316M mutations, the relative proportion between the different isoforms is shifted toward the expression of NF1ΔE7, while in the presence of N335N variant, the NF1ΔE7 expression is abolished.
In conclusion, it appears mandatory to investigate the role of each nucleotide change within the NF1 coding sequence, since a significant proportion of NF1 exon 7 mutations affects pre-mRNA splicing, by disrupting exonic splicing motifs and modifying the delicate balance between aberrantly and correctly spliced transcripts.
Alternative splicing, the process by which exons are included or excluded in the mature mRNA, is an important mechanism whereby different transcripts are generated from the same gene unit. In fact, most human genes are transcribed in multiple alternative mRNAs, according to different regulatory programs, resulting in functionally different protein isoforms . In the best characterized models of vertebrate cell-specific alternative splicing, post-transcriptional regulation is tissue or developmental stage specific and may be mediated by intronic and exonic cis elements. These elements, which are important for the correct splice-site identification, can act by stimulating (exonic splicing enhancers, ESEs) or repressing (exonic splicing silencers, ESSs) the exon's splicing .
Neurofibromatosis type 1 (NF1, MIM#162200) is one of the most common autosomal dominant disorders, affecting about 1:3,500 individuals in all ethnic groups. The NF1 gene is approximately 280 kb in size and maps to chromosome 17q11.2 [3–5]. NF1 contains 60 exons, with an 11- to 13-kb transcript and an open reading frame coding for 2,818 amino acids . The disease is fully penetrant and the diagnosis of NF1 is based on the clinical criteria recommended by NIH Consensus Conference (Stumpf, et al., 1988), which include multiple cafè-au-lait spots, cutaneous or subcutaneous neurofibromas, plexiform neurofibromas, axillary or inguinal freckling, optic gliomas, and iris Lish nodules. Although NF1 mutations are distributed along the entire coding sequence, no genotype-phenotype correlation has been found so far , with the exception of the recurrent and atypical deletions underlying NF1 microdeletion syndrome .
The NF1 gene is ubiquitously expressed, and four normal in-frame NF1 splice isoforms are known, brain specific 9br isoform (30 bp) , exon 10a-2 isoform (45 bp) , exon 23a isoform (63 bp), which are expressed in all tissues at various levels , and the muscle specific exon 48a isoform (54 bp) . In addition, several other alternative transcripts have been described, including ex29-, ex30-, ex29/30- and the N-isoform [13, 14].
Mutation analysis has shown that approximately 50% of NF1 mutations result in splicing alterations [15–18]. In some cases, splicing mutations do not occur at the conserved AG/GT dinucleotides of the splice sites. For example, mutations leading to stop codons in exon 7 and 37 of NF1 gene have been reported to be involved in exon skipping [18–21]. In addition, mutational analysis of the NF1 gene disclosed several additional splice variants in which specific exons are skipped in fresh lymphocytes of unaffected persons, albeit typically at low levels [22, 23]. A number of studies have also reported that some of these transcripts are more abundant when RNA from aged blood or from blood kept at non-physiological temperatures is analyzed [15, 22, 24, 25].
The expression of an alternative transcript lacking exon 7 has been demonstrated [17, 26]. Indeed, NF1 exon 7 displays weakly defined exon-intron boundaries, and is particularly prone to aberrant splicing. In the present study, we have used in silico and in vitro analysis to evaluate the functional consequences on gene expression of four nucleotide variants detected in NF1 exon 7, including a nonsense mutation (R304X), a missense mutation (L316M), and two silent changes (Q315Q and N335N). Since both Q315Q and L316M mutations were together in cis in the same patients , our analysis aim was to understand their effect together and independently.
DNA mutation analysis
The nucleotide variants investigated in this study include a recurrent nonsense mutation [c.910C>T (R304X)] [22, 27–29] and a novel silent change [c.1005T>C (N335N)] identified using denaturing high performance liquid chromatography (dHPLC) followed by bidirectional sequencing, as well as a silent change and a missense mutation occurring together in cis in the same patient [c.945G>A/c.946C>A (Q315Q/L316M)], and previously reported by us . PCR conditions, amplicon length, and resolution temperatures for dHPLC analysis are reported elsewhere [28, 30]. The N335N silent change was found in a two generation NF1 family (family NF-01) carrying another NF1 gene mutation, a frameshift deletion (c.476delC) in exon 4a. The silent change N335N was found in the proband (II-2) and her child (III-1), both affected by NF1, and in the proband's father (I-1) presenting out of NF1 clinical signs only three cutaneous neurofibromas. Frameshift mutation c.476delC was detected in the proband and her child, but not in the proband's father. Both N335N and c.476delC were not found in 200 healthy subjects. Microsatellite analysis performed using 10 markers tightly linked to the NF1 locus (D17S1873, D17S841, D17S1863, D17S635, D17S1166, IVS-38, 3'NF1-1, D17S1800, 3'NF1-2, D17S798) showed that N335N and c.476delC co-segregated on the same chromosome 17. Pedigree microsatellite analysis details of family NF-01 are shown in Figure 1. The project was approved by the institutional review board and all participants provided informed consent.
Prediction of exonic splice enhancers
The NF1 exon 7 coding sequence including sequence changes Q315Q, L316M, R304X and N335N was analyzed using the online ESEfinder web interface (Release 1.1)  which performs searches for putative ESEs in query sequences using weight matrices corresponding to the motifs of four different human serine/arginine-residue proteins (SR proteins): SF2/ASF, SC35, SRp40, and SRp55. The score for a given sequence was considered to be potentially significant if greater than the default threshold value defined in the input page. These values were set as follows: 1.956 for SF2/ASF, 2.383 for SC35, 2.670 for SRp40, and 2.676 for the SRp55.
Minigene construction and expression
Minigene constructs were generated and transfected into cos-1 cells. Exon 7 of the NF1 gene and flanking intronic sequence were amplified using as template the DNA of a non-affected individual. Primers NF1-7F (5'-cacacactcgagAACAGCTTGTTTGGGAAGGA-3') and NF1-7R (5'-cacacaggatccGGCCCTAATTGCCACATTATT-3') were used to generate a fragment containing exon 7 as well as 258 bp of 5' and 248 bp of 3' of intronic flanking sequence. The forward primer contains a 5'-XhoI linker, whereas the reverse primer contains a 5'-BamHI linker. The PCR amplification reaction was performed on 100 ng genomic DNA in a standard 50 μL volume, containing 1× Optimase reaction buffer (Transgenomic, Crewe, UK), 1.5 mM MgSO4, 1 μM of each primer, 0.2 μM dNTPs, and 2.5 U Optimase TaqDNA polymerase (Transgenomic, Crewe, UK) in a 9700 thermal cycler (Applied Biosystem, FosterCity, CA). Thermal conditions were 35 cycles of 95°C for 30 seconds, 57°C for 30 seconds, and 72°C for 1 minute, preceded by 5 minutes at 95°C and followed by a final elongation step at 72°C for 10 minutes.
After confirmation of successful amplification through detection of the expected 680 bp band on the agarose gel, the products were digested with XhoI (Promega, Madison, WI, USA) and BamHI (Promega, Madison, WI, USA) restriction enzymes. The exon trapping expression vector pSPL3 (Invitrogen Corporation, Carlsbad, CA) contains a replicon and Apr marker for growth in E. coli, an SV40 segment for replication and transcription in cos-1 cells, HIV-1 tat splicing signals, and a multiple cloning site. The tat segment contains an intron, splice donor (SD) and splice acceptor (SA) sites, and flanking exon sequences. The insert was directly ligated between the SD and SA sites into the XhoI/BamHI restriction points. Ligation into pSPL3 was performed at room temperature for 30 min, using T4 DNA ligase (Invitrogen Corporation, Carlsbad, CA). E. coli DH5α competent cells (Invitrogen Corporation, Carlsbad, CA) were transformed with the plasmid constructs and plated overnight. The resulting clones were checked for fragment orientation and sequenced. Minigene constructs were isolated using a midiprep kit (Qiagen, Hilden, Germany). The resulting pSPL3-NF1-7-wt minigene construct is shown in Figure 2. The sequence changes were independently introduced in the pSPL3-NF1-7-wt by means of the QuickChange Site-Direct Mutagenesis Kit (Stratagene, La Jolla, CA), as instructed by the manufacturer. The changes made in pSPL3-NF1-7-wt are shown in Figure 2, and Table 1 lists the mutagenic primers used. All the mutants were sequenced to confirm that only the desired changes were introduced, and were then isolated with a miniprep kit (Qiagen, Hilden, Germany). The minigene constructs containing either the wild type sequence or a exon 7 variant were transfected into cos-1 cells by electroporation.
Approximately 2.5 × 106 cos-1 cells were trypsinized, washed with cold PBS 1× buffer (Mg2+ and Ca2+ free) and then resuspended in 800 μl of the same PBS solution. Two μg of normal or mutant minigene DNAs were added to the cell suspension in 0.4 cm gap sterile cuvettes. The cell-DNA mixture was then incubated on ice for 10 min and electroporated at 4°C using a Bio-Rad Gene Pulser II (Bio-Rad, Richmond, CA) at a setting of 300 Volts, 250 μF. Electroporated cells were then incubated on ice for 10 min, diluted 1:20 in complete medium and cultured in T-75 flasks at 37°C in a humidified atmosphere under 5% CO2 for at least 48 h. Total cellular RNA from the transfected cos-1 cells was purified by the acid guanidine-phenol-chloroform method and used for RT-PCR to confirm the splicing patterns. First strand cDNA was synthesized from 2 to 3 μg of total RNA by random-primed reverse transcription with Superscript II Reverse Transcriptase (Invitrogen Corporation, Carlsbad, CA). To evaluate the pattern of transcripts from the transfected minigenes, the following vector-specific primers were used for RT-PCR amplification: a forward primer SD6 (5'-TCTGAGTCACCTGGACAACC-3') and a reverse primer SA2 (5'-ATCTCAGTGGTATTTGTGAGC-3'). The PCR amplification reaction was performed as follows: in 50 μL volume, 2 μl of cDNA, 5 μl of Expand High Fidelity buffer 3 (Roche, Mannheim, Germany), 1 μM of each primer, 0.8 μM dNTPs, and 2.6 U Expand High Fidelity enzyme mix (Roche, Mannheim, Germany) in a 9700 (Applied Biosystem, FosterCity, CA, USA) thermal cycler. Thermal conditions were 30 cycles of 95°C for 30 seconds, 58°C for 30 seconds, and 68°C for 1 minute, preceded by 2 minutes at 95°C, and followed by a final elongation step at 68°C for 10 minutes. The PCR products were separated by electrophoresis on a 3% agarose gel and each band signal was quantified by Quantity One software (Bio-Rad, Richmond, CA). All transcripts were analyzed by sequencing.
The human NF1 exon\intron 7 sequence was derived from genomic NF1 sequence [GenBank: NM_000267].
Using bioinformatic identification of putative regulatory sequences and functional minigene analysis we have evaluated the effect on the NF1 pre-mRNA splicing process of four nucleotide variants detected in NF1 exon 7 of patients with a clinical diagnosis of NF1. These included a novel silent change (N335N), and three mutations (R304X, Q315Q and L316M) previously reported [22, 27–29]. All variants, except R304X, were predicted by ESEfinder to destroy or create an ESE element. The predictions for the mutated exonic sequences are summarized in Table 2. The R304X mutation does not add or abolish any ESE sequence, while the Q315Q and the L316M changes abolish respectively a SF2/ASF and a SC35 ESE motif. As a consequence of the N335N change, a new ESE motif for the SRp40 protein is introduced.
To confirm these predictions, we performed in vitro experiments that tested the splicing enhancement capacity of the wild type and mutant sequences comprising the predicted ESEs motifs. Figure 3 shows the amplified PCR spliced products of the minigene constructs, produced from primer SD6 and primer SA2 after the transfection. Each band was quantified by the Quantity One software (Bio-Rad, Richmond, CA). The wild type construct showed three splicing products in different proportions: a fragment of 435 bp (86%), another 403-bp fragment (2%), and a third of 261 bp (11%). After gel extraction, each PCR product was directly sequenced. The 435 bp PCR product included exon 7 and its intronic flanking sequences, the 403 bp fragment matched with the expected exon 7, but lacked the last 32 nucleotides, and the third band of 261 bp corresponded to a transcript lacking the entire sequence of exon 7 (NF1ΔE7). The construct containing both the Q315Q and the L316M changes showed only a single PCR product matching with NF1ΔE7. The constructs containing respectively the L316M, the Q315Q and the R304X changes disclosed the same three bands as the wild type, but in different proportions, with the NF1ΔE7 strongly represented, and corresponding respectively to 75%, 86% and 44% of the three fragments. The L316M, Q315Q and R304X constructs expressed mRNAs that included the 403-bp exon 7 product. The construct containing the N335N change showed only the 435 bp PCR product corresponding to exon 7 and its flanking intronic sequences.
In this study four NF1 exon 7 variants detected in three different patients with NF1 have been investigated using in silico prediction to determine their effect on ESE sequences and by in vitro studies to assess their functional significance on the splicing process. The mutant series included a recurrent nonsense mutation [c.910C>T (R304X)] [22, 29], a novel silent change [c.1005T>C (N335N)], and a silent change and a missense mutation [c.945G>A/c.946C>A (Q315Q/L316M)] occurring together in cis in the same patient . In silico analysis showed that the presence of the Q315Q and the L316M mutations each caused the loss of an ESE motif (SF2/ASF and SC35 protein respectively), suggesting a role for these nucleotide sequences in exon 7 retention in the NF1 message. Using the minigene assay, we demonstrated that the presence of the Q315Q and L316M mutations consistently shifted the proportion between the isoforms toward the expression of NF1ΔE7, enhancing exon 7 skipping and resulting in decreased levels of full length neurofibromin in vivo. Whereas the construct containing both Q315Q and L316M changes showed only the presence of NF1ΔE7 transcript, the constructs containing either the Q315Q or the L316M alone retained some wild type transcript. As predicted by ESEfinder software, Q315Q and L316M are located in two ESE motifs recognized by two different SR proteins. Each mutation individually affects only one ESE element and thus allows for the expression of a small amount of full length transcript. Conversely, the presence of both changes destroys two ESE elements and completely inhibits any residual inclusion of exon 7. According to the minigene results, RT-PCR analysis of the patient carrying the Q315Q and L316M mutations showed the presence of an additional transcript lacking NF1 exon 7, which accounts for 61% of the entire message . Based on the ESEfinder results, the R304X mutation was not found to substantially affect any ESE element. However, when this mutation was placed in a heterologous minigene context, it was able to enhance the skipping of exon 7. Similar results were obtained by others expressing the R304X into a different vector . Using RESCUE-ESE, a computational method alternative to ESEfinder, the same authors were able to predict that mutation R304X causes the loss of a cluster of two hexamer motifs in exon 7, but concomitantly creates two novel hexamer motifs. They speculated that the loss of the predicted ESE sites as a result of R304X has a dominant effect over the increase in the score of the new ESE sites. The results of another study suggested that R304X mutation substantially alters the predicted minimum free-energy structure of exon 7 , although this result has not been confirmed by others .
In our experiment, the search for possible SR protein binding motifs in NF1 exon 7 using the ESEfinder web interface  showed the presence of three ESE elements in the wild type sequence. This finding and the knowledge of the weak exon-intron boundaries of NF1 exon 7 [17, 26, 27] are consistent with the occurrence of exon 7 alternative splicing also seen in non-pathological conditions. Although exon 7 skipping without sequence alterations is mainly induced by stressful factors such as the amplification of RNA from aged blood [22, 23], a small amount of NF1ΔE7 isoform (<1%) was also observed from RNA extracted immediately following blood drawing [17, 26]. Consistently, the wild type minigene construct analysis showed the presence of three transcripts, one including exon 7, one lacking the entire exon 7 and a third lacking the last 32 nucleotides of exon 7, indicating the in vitro existence of an equilibrium between splicing products of the NF1 gene. Regarding the wild type construct, NF1ΔE7 is less expressed than the other transcripts. R304X in silico analysis did not produce relevant results, while in vitro analysis shows a reduced expression of full length transcripts, even if weaker than Q315Q and L316M mutations. The transcript lacking the last 32 nucleotides of exon 7 was expressed by all minigene constructs. To our knowledge, a messanger lacking the last 32 nucleotides of exon 7 has not been observed previously . Although this transcript could represent a NF1 isoform expressed at a very low level, we cannot exclude that this message is an artefact of the minigene assay. In fact, the proportion of exon-including and exon lacking spliced minigene transcripts does not always correspond precisely to the levels measured in patients . It should also be considered that in vivo a NF1 message lacking 32 nucleotides would result in an out-of-frame transcript, possibly degraded by nonsense-mediated mRNA decay pathway.
Whereas Q315Q, L316M and R304X mutations shift the equilibrium between NF1 transcripts towards the skipping of exon 7, minigene analysis of the N335N silent change demonstrated that this variant completely shifts the equilibrium between the NF1 messengers towards the full length transcript, being the only detectable transcript. It was previously reported that the sequence around nucleotide 1007G had a weak ESE activity and a mutation at this site caused only a mild splicing defect . It was speculated that this ESE sequence was part of a composite exonic regulatory element of splicing (CERES), defined as an exonic regulatory splicing element, having overlapping enhancer and silencer functions . Using ESEfinder analysis we were able to predict that the consequence of the N335N change is the introduction of a novel ESE motif for the SRp40 protein. However, ExonScan [34, 35], another algorithm which simulates splicing based on known or putative splicing-related motifs, gives a different prediction showing that the mutation (c.910T>C) N335N results in the loss of an ESS (ATCAAT) motif in exon 7, but concurrently creates another novel ESS motif (AACTGG) (data not shown). As previously hypothesized for the R304X mutation and also for the N335N silent change, the loss of a putative splicing-related motif has a dominant effect over the acquisition of a new splicing-related motif. Accordingly, in vitro analysis has proven that nucleotide 1005T could partially belong to an ESS, since a change destroying this sequence completely eliminates NF1ΔE7 expression. The presence of some ESE/ESS elements in NF1 exon 7 and the existence of NF1ΔE7 in healthy individuals suggest that the splicing of exon 7 is subject to fine regulation. Notably, NF1 exon 7 is in-frame, thus the resulting protein could theoretically retain some function, missing only the residues encoded by this exon. As reported in our previous studies even if NF1 exon 7 does not belong to the neurofibromin active site it is a mutation hot spot, suggesting that exon 7 is important for the pathogenesis of NF1 [27, 30, 36]. It could be speculated that the NF1ΔE7 isoform could have a specific role in cellular metabolism, or that its' expression could represent a mechanism to regulate intracellular neurofibromin levels. However, the N335N change, which eliminates the in vitro expression of NF1ΔE7, was detected in a subject who, before clinical evaluation of full blown NF1 in his daughter and grandson, did not fit with the NF1 diagnostic criteria, manifesting only three cutaneous neurofibromas. This incomplete NF1 phenotype could be explained by the presence of undetectable mosaicism for the c.476delC mutation which was absent in the germline DNA of this subject, but was found in cis with the inherited N335N change in the germline DNA of his daughter and grandson. Unfortunately, none of the neurofibromas could be verified by biopsy to prove this hypothesis. Alternatively, it could be possible that c.476delC frameshift mutation could have arisen de novo on the same allele of N335N during spermatogenesis. However, the contribution of the N335N variant to the NF1 clinical phenotype is difficult to infer. In fact, it is difficult to imagine for this silent substitution to have the same effect on exon 7 skipping regardless of the upstream c.476delC frameshift mutation, as this type of mutations are expected to trigger NMD of the corresponding allele. However, lack of this change in normal individuals appears to rule out its polymorphic nature. Conversely, a pathogenic role for the Q315Q and L316M is easier to infer, as they are de novo mutations belonging to contiguous codons . Moreover, as shown by the minigene analysis, the presence of each mutation individually is sufficient to enhance the expression of NF1ΔE7, thus suggesting that each of these mutations per se can result in NF1. In the present case they acted additively to cause the NF1 patient phenotype.
In the present study we observed that a significant proportion of neurofibromatosis associated mutations (missense, nonsense or silent changes) residing in NF1 exon 7 affect pre-mRNA splicing by disrupting exonic splicing motifs, shifting the transcripts balance towards aberrantly spliced transcripts and producing two types of transcripts: those with the original mutation, and those with the skipping . Therefore, it should be mandatory to investigate each nucleotide change occurring within the NF1 coding sequence for their potential pathogenetic role in NF1.
Matlin AJ, Clark F, Smith CW: Understanding alternative splicing: towards a cellular code. Nat Rev Mol Cell Biol. 2005, 6 (5): 386-398. 10.1038/nrm1645.
Cartegni L, Chew SL, Krainer AR: Listening to silence and understanding nonsense: exonic mutations that affect splicing. Nat Rev Genet. 2002, 3 (4): 285-298. 10.1038/nrg775.
Cawthon RM, Weiss R, Xu GF, Viskochil D, Culver M, Stevens J, Robertson M, Dunn D, Gesteland R, O'Connell P, et al: A major segment of the neurofibromatosis type 1 gene: cDNA sequence, genomic structure, and point mutations. Cell. 1990, 62 (1): 193-201. 10.1016/0092-8674(90)90253-B.
Viskochil D, Buchberg AM, Xu G, Cawthon RM, Stevens J, Wolff RK, Culver M, Carey JC, Copeland NG, Jenkins NA, et al: Deletions and a translocation interrupt a cloned gene at the neurofibromatosis type 1 locus. Cell. 1990, 62 (1): 187-192. 10.1016/0092-8674(90)90252-A.
Wallace MR, Marchuk DA, Andersen LB, Letcher R, Odeh HM, Saulino AM, Fountain JW, Brereton A, Nicholson J, Mitchell AL, et al: Type 1 neurofibromatosis gene: identification of a large transcript disrupted in three NF1 patients. Science. 1990, 249 (4965): 181-186. 10.1126/science.2134734.
Marchuk DA, Saulino AM, Tavakkol R, Swaroop M, Wallace MR, Andersen LB, Mitchell AL, Gutmann DH, Boguski M, Collins FS: cDNA cloning of the type 1 neurofibromatosis gene: complete sequence of the NF1 gene product. Genomics. 1991, 11 (4): 931-940. 10.1016/0888-7543(91)90017-9.
Castle B, Baser ME, Huson SM, Cooper DN, Upadhyaya M: Evaluation of genotype-phenotype correlations in neurofibromatosis type 1. J Med Genet. 2003, 40 (10): e109-10.1136/jmg.40.10.e109.
Venturin M, Guarnieri P, Natacci F, Stabile M, Tenconi R, Clementi M, Hernandez C, Thompson P, Upadhyaya M, Larizza L, Riva P: Mental retardation and cardiovascular malformations in NF1 microdeleted patients point to candidate genes in 17q11.2. J Med Genet. 2004, 41 (1): 35-41. 10.1136/jmg.2003.014761.
Danglot G, Regnier V, Fauvet D, Vassal G, Kujas M, Bernheim A: Neurofibromatosis 1 (NF1) mRNAs expressed in the central nervous system are differentially spliced in the 5' part of the gene. Hum Mol Genet. 1995, 4 (5): 915-920. 10.1093/hmg/4.5.915.
Kaufmann D, Muller R, Kenner O, Leistner W, Hein C, Vogel W, Bartelt B: The N-terminal splice product NF1-10a-2 of the NF1 gene codes for a transmembrane segment. Biochem Biophys Res Commun. 2002, 294 (2): 496-503. 10.1016/S0006-291X(02)00501-6.
Andersen LB, Ballester R, Marchuk DA, Chang E, Gutmann DH, Saulino AM, Camonis J, Wigler M, Collins FS: A conserved alternative splice in the von Recklinghausen neurofibromatosis (NF1) gene produces two neurofibromin isoforms, both of which have GTPase-activating protein activity. Mol Cell Biol. 1993, 13 (1): 487-495.
Gutmann DH, Geist RT, Rose K, Wright DE: Expression of two new protein isoforms of the neurofibromatosis type 1 gene product, neurofibromin, in muscle tissues. Dev Dyn. 1995, 202 (3): 302-311.
Park VM, Kenwright KA, Sturtevant DB, Pivnick EK: Alternative splicing of exons 29 and 30 in the neurofibromatosis type 1 gene. Hum Genet. 1998, 103 (4): 382-385. 10.1007/s004390050837.
Suzuki Y, Suzuki H, Kayama T, Yoshimoto T, Shibahara S: Brain tumors predominantly express the neurofibromatosis type 1 gene transcripts containing the 63 base insert in the region coding for GTPase activating protein-related domain. Biochem Biophys Res Commun. 1991, 181 (3): 955-961. 10.1016/0006-291X(91)92029-J.
Ars E, Serra E, de la Luna S, Estivill X, Lazaro C: Cold shock induces the insertion of a cryptic exon in the neurofibromatosis type 1 (NF1) mRNA. Nucleic Acids Res. 2000, 28 (6): 1307-1312. 10.1093/nar/28.6.1307.
Fahsold R, Hoffmeyer S, Mischung C, Gille C, Ehlers C, Kucukceylan N, Abdel-Nour M, Gewies A, Peters H, Kaufmann D, Buske A, Tinschert S, Nurnberg P: Minor lesion mutational spectrum of the entire NF1 gene does not explain its high mutability but points to a functional domain upstream of the GAP-related domain. Am J Hum Genet. 2000, 66 (3): 790-818. 10.1086/302809.
Vandenbroucke I, Vandesompele J, De Paepe A, Messiaen L: Quantification of NF1 transcripts reveals novel highly expressed splice variants. FEBS Lett. 2002, 522 (1-3): 71-76. 10.1016/S0014-5793(02)02887-9.
Hoffmeyer S, Nurnberg P, Ritter H, Fahsold R, Leistner W, Kaufmann D, Krone W: Nearby stop codons in exons of the neurofibromatosis type 1 gene are disparate splice effectors. Am J Hum Genet. 1998, 62 (2): 269-277. 10.1086/301715.
Vandenbroucke, Vandesompele J, Paepe AD, Messiaen L: Quantification of splice variants using real-time PCR. Nucleic Acids Res. 2001, 29 (13): E68-8. 10.1093/nar/29.13.e68.
Messiaen L, Callens T, De Paepe A, Craen M, Mortier G: Characterisation of two different nonsense mutations, C6792A and C6792G, causing skipping of exon 37 in the NF1 gene. Hum Genet. 1997, 101 (1): 75-80. 10.1007/s004390050590.
Baralle M, Skoko N, Knezevich A, De Conti L, Motti D, Bhuvanagiri M, Baralle D, Buratti E, Baralle FE: NF1 mRNA biogenesis: effect of the genomic milieu in splicing regulation of the NF1 exon 37 region. FEBS Lett. 2006, 580 (18): 4449-4456. 10.1016/j.febslet.2006.07.018.
Messiaen LM, Callens T, Mortier G, Beysen D, Vandenbroucke I, Van Roy N, Speleman F, Paepe AD: Exhaustive mutation analysis of the NF1 gene allows identification of 95% of mutations and reveals a high frequency of unusual splicing defects. Hum Mutat. 2000, 15 (6): 541-555. 10.1002/1098-1004(200006)15:6<541::AID-HUMU6>3.0.CO;2-N.
Thomson SA, Fishbein L, Wallace MR: NF1 mutations and molecular testing. J Child Neurol. 2002, 17 (8): 555-61; discussion 571-2, 646-51.
Ars E, Kruyer H, Morell M, Pros E, Serra E, Ravella A, Estivill X, Lazaro C: Recurrent mutations in the NF1 gene are common among neurofibromatosis type 1 patients. J Med Genet. 2003, 40 (6): e82-10.1136/jmg.40.6.e82.
Wimmer K, Eckart M, Rehder H, Fonatsch C: Illegitimate splicing of the NF1 gene in healthy individuals mimics mutation-induced splicing alterations in NF1 patients. Hum Genet. 2000, 106 (3): 311-313. 10.1007/s004390051043.
Wimmer K, Eckart M, Stadler PF, Rehder H, Fonatsch C: Three different premature stop codons lead to skipping of exon 7 in neurofibromatosis type I patients. Hum Mutat. 2000, 16 (1): 90-91. 10.1002/1098-1004(200007)16:1<90::AID-HUMU20>3.0.CO;2-J.
Colapietro P, Gervasini C, Natacci F, Rossi L, Riva P, Larizza L: NF1 exon 7 skipping and sequence alterations in exonic splice enhancers (ESEs) in a neurofibromatosis 1 patient. Hum Genet. 2003, 113 (6): 551-554. 10.1007/s00439-003-1009-2.
De Luca A, Schirinzi A, Buccino A, Bottillo I, Sinibaldi L, Torrente I, Ciavarella A, Dottorini T, Porciello R, Giustini S, Calvieri S, Dallapiccola B: Novel and recurrent mutations in the NF1 gene in Italian patients with neurofibromatosis type 1. Hum Mutat. 2004, 23 (6): 629-10.1002/humu.9245.
Pros E, Larriba S, Lopez E, Ravella A, Gili ML, Kruyer H, Valls J, Serra E, Lazaro C: NF1 mutation rather than individual genetic variability is the main determinant of the NF1-transcriptional profile of mutations affecting splicing. Hum Mutat. 2006, 27 (11): 1104-1114. 10.1002/humu.20396.
De Luca A, Buccino A, Gianni D, Mangino M, Giustini S, Richetta A, Divona L, Calvieri S, Mingarelli R, Dallapiccola B: NF1 gene analysis based on DHPLC. Hum Mutat. 2003, 21 (2): 171-172. 10.1002/humu.9111.
Cartegni L, Wang J, Zhu Z, Zhang MQ, Krainer AR: ESEfinder: A web resource to identify exonic splicing enhancers. Nucleic Acids Res. 2003, 31 (13): 3568-3571. 10.1093/nar/gkg616.
Zatkova A, Messiaen L, Vandenbroucke I, Wieser R, Fonatsch C, Krainer AR, Wimmer K: Disruption of exonic splicing enhancer elements is the principal cause of exon skipping associated with seven nonsense or missense alleles of NF1. Hum Mutat. 2004, 24 (6): 491-501. 10.1002/humu.20103.
Pagani F, Stuani C, Tzetis M, Kanavakis E, Efthymiadou A, Doudounakis S, Casals T, Baralle FE: New type of disease causing mutations: the example of the composite exonic regulatory elements of splicing in CFTR exon 12. Hum Mol Genet. 2003, 12 (10): 1111-1120. 10.1093/hmg/ddg131.
Wang Z, Rolish ME, Yeo G, Tung V, Mawson M, Burge CB: Systematic identification and analysis of exonic splicing silencers. Cell. 2004, 119 (6): 831-845. 10.1016/j.cell.2004.11.010.
Origone P, De Luca A, Bellini C, Buccino A, Mingarelli R, Costabel S, La Rosa C, Garre C, Coviello DA, Ajmar F, Dallapiccola B, Bonioli E: Ten novel mutations in the human neurofibromatosis type 1 (NF1) gene in Italian patients. Hum Mutat. 2002, 20 (1): 74-75. 10.1002/humu.9039.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2350/8/4/prepub
We thank patients and families enrolled in this study and physicians who referred these subjects. This work was supported in part by the Italian Ministry of Health, Ricerca Corrente 2005–2006.
The author(s) declare that they have no competing interests.
IB carried out the molecular genetic studies and drafted the manuscript. ADL conceived of the study and helped to draft the manuscript. AS and VG carried out the molecular genetic studies. IT collected the samples and helped to draft the manuscript. SC and CG collected the samples. LL collected the samples and revised the manuscript. AP participated in the design of the study and revised the manuscript. BD provided overall direction to the project and revised the manuscript as well as approved the final version of it. All authors read and approved the final manuscript.