Next Article in Journal
Mutant CAG Repeats Effectively Targeted by RNA Interference in SCA7 Cells
Previous Article in Journal
Functions of the RNA Editing Enzyme ADAR1 and Their Relevance to Human Diseases
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Hybridization Capture-Based Next-Generation Sequencing to Evaluate Coding Sequence and Deep Intronic Mutations in the NF1 Gene

by
Karin Soares Cunha
1,2,3,*,
Nathalia Silva Oliveira
4,
Anna Karoline Fausto
4,
Carolina Cruz De Souza
5,
Audrey Gros
6,7,
Thomas Bandres
6,
Yamina Idrissi
7,
Jean-Philippe Merlio
6,7,
Rodrigo Soares De Moura Neto
8,
Rosane Silva
9,
Mauro Geller
10,11 and
David Cappellen
6,7
1
Graduate Program in Pathology, School of Medicine, Universidade Federal Fluminense, Niterói 24033-900, Brazil
2
Department of Pathology, School of Medicine, Universidade Federal Fluminense, Niterói 24033-900, Brazil
3
Neurofibromatosis National Center (Centro Nacional de Neurofibromatose), Rio de Janeiro 20011-330, Brazil
4
Anatomy Pathology Service, Hospital Universitário Antônio Pedro, Universidade Federal Fluminense, Niterói 24033-900, Brazil
5
School of Biomedicine, Universidade Federal Fluminense, Niterói 24210-130, Brazil
6
Service de Biologie des Tumeurs, Centre Hospitalier Universitaire de Bordeaux, Hôpital du Haut Lévêque, Pessac F-33604, France
7
Inserm (Institut National de la Santé et de la Recherche Médicale) U1053, Bordeaux Research in Translational Oncology (BaRITON) and University of Bordeaux, Bordeaux F-33076, France
8
Biology Institute, Universidade Federal do Rio de Janeiro, Rio de Janeiro 21941-599, Brazil
9
Carlos Chagas Filho Biophysics Institute, Universidade Federal do Rio de Janeiro, Rio de Janeiro 21941-599, Brazil
10
Department of Immunology and Microbiology, School of Medicine, Centro Universitário Serra dos Órgãos, Teresópolis 25964-004, Brazil
11
Martagão Gesteira Child Care and Pediatrics Institute, Universidade Federal do Rio de Janeiro, Rio de Janeiro 21941-912, Brazil
*
Author to whom correspondence should be addressed.
Genes 2016, 7(12), 133; https://doi.org/10.3390/genes7120133
Submission received: 23 October 2016 / Revised: 29 November 2016 / Accepted: 7 December 2016 / Published: 17 December 2016
(This article belongs to the Section Human Genomics and Genetic Diseases)

Abstract

:
Neurofibromatosis 1 (NF1) is one of the most common genetic disorders and is caused by mutations in the NF1 gene. NF1 gene mutational analysis presents a considerable challenge because of its large size, existence of highly homologous pseudogenes located throughout the human genome, absence of mutational hotspots, and diversity of mutations types, including deep intronic splicing mutations. We aimed to evaluate the use of hybridization capture-based next-generation sequencing to screen coding and noncoding NF1 regions. Hybridization capture-based next-generation sequencing, with genomic DNA as starting material, was used to sequence the whole NF1 gene (exons and introns) from 11 unrelated individuals and 1 relative, who all had NF1. All of them met the NF1 clinical diagnostic criteria. We showed a mutation detection rate of 91% (10 out of 11). We identified eight recurrent and two novel mutations, which were all confirmed by Sanger methodology. In the Sanger sequencing confirmation, we also included another three relatives with NF1. Splicing alterations accounted for 50% of the mutations. One of them was caused by a deep intronic mutation (c.1260 + 1604A > G). Frameshift truncation and missense mutations corresponded to 30% and 20% of the pathogenic variants, respectively. In conclusion, we show the use of a simple and fast approach to screen, at once, the entire NF1 gene (exons and introns) for different types of pathogenic variations, including the deep intronic splicing mutations.

1. Introduction

Neurofibromatosis 1 (NF1; MIM: 162200) is an autosomal dominant disorder with a complete penetrance and highly variable expressivity. It is one of the most common genetic disorders with an incidence of 1 case in 2000 births [1]. Neurofibromas, Lisch nodules, café-au-lait spots, as well as inguinal and axillary “freckles” develop in the majority of individuals with NF1. Other important features include learning disabilities, optic gliomas, specific bone lesions, and risk for the development of malignant neoplasms, mainly the malignant peripheral nerve sheath tumors (MPNSTs) [2].
NF1 is caused by a spectrum of mutations that affect the NF1 gene, which is a tumor suppressor gene located at chromosome 17q11.2. It spans 282 kb and comprises 60 exons (57 constitutive and 3 alternatively spliced exons), being one of the largest human genes [3]. The NF1 gene presents one of the highest frequencies of spontaneous mutations in the human genome. Therefore, half of individuals with NF1 present no family history of the disorder [4].
To date, according to the online Human Gene Mutation Database (data accessed on November 2016), more than 2600 different NF1 mutations dispersed throughout the gene have been reported and include missense/nonsense, splicing, small insertions or deletions (indels), and also gross chromosomal alterations. Although no true hotspots have been identified in the NF1 gene, recurrent mutations are not uncommon in individuals with NF1 [5].
Identification and characterization of NF1 gene mutations represent a considerable research and diagnostic challenge because of its large size, existence of a number of highly homologous pseudogenes located throughout the human genome (in chromosomes 2, 12, 14, 15, 18, 21, and 22), absence of mutational hotspots as well as the diversity of mutation types, including deep intronic splicing mutations [6,7]. The NF1 gene mutation detection rate in NF1 patients who meet the clinical diagnostic criteria is reported to be 50%–95% (reviewed by Upadhyaya [8]). The highest mutation detection rate reported so far (≈95%) is based on a multistep mutation analysis at the genomic and RNA level [9]. It includes a massive PCR genomic DNA (gDNA)/complementary DNA (cDNA) preparation and Sanger sequencing of multiple regions of the NF1 gene, which are associated with other methodologies to identify copy number alterations (high-performance liquid chromatography, multiplex ligand probe amplification (MLPA), microsatellite marker analysis, and fluorescence in situ hybridization (FISH)).
Therefore, current NF1 molecular diagnostic testing based on Sanger sequencing is labor-intensive, time-consuming (usually takes 3–4 weeks on average), and expensive [3,8]. Moreover, even using an exhaustive analysis, mutation screening tests do not achieve a sensitivity of 100%. Possibly many, if not all, cases of patients with classical NF1, with no NF1 mutation detected in the genetic tests, are related to mutations located in NF1 gene regions that are not commonly screened (e.g., the flanking 5′- and 3′-untranslated regions (UTRs)) [8]. It is also possible that the disease-causing pathogenic variations are related to deep intronic mutations [8], which are missed when traditional gDNA-based screening techniques, that sequence only the exonic regions and conserved splice sites, are applied.
In the last few years, DNA sequencing has been passing through a revolution with the use of next-generation sequencing (NGS). It permits a fast and simple sequencing, making possible the analysis of large genomes. With NGS, rather than just sequence the coding exons and conserved splice sites, it is possible to screen the entire sequence of a target large gene from a gDNA sample, including the intronic regions. Therefore, NGS seems to be a very useful tool to sequence the NF1 gene [10]. In this study, we aimed to evaluate the use of hybridization capture-based next-generation sequencing to screen not only coding, but also noncoding NF1 regions.

2. Materials and Methods

This study was approved by the ethical committee of Antônio Pedro University Hospital of Universidade Federal Fluminense, Niterói, Brazil (#121/11).

2.1. Patients

Some clinical features of the participants of this study are shown in Table 1. NGS samples were obtained from 12 individuals with NF1: 11 unrelated NF1 patients and one relative with NF1 (NF.85; son of the patient NF.84). For the Sanger methodology (see Section 2.6), we have also included another three relatives with NF1, as follows: NF.14 (cousin of patient NF.84), NF.87 (son of patient NF.91), and NF.27 (sister of patient NF.26). All the individuals had NF1 according to the National Institutes of Health (NIH) clinical criteria [11].
Extraction of DNA from leukocytes and purification were performed with QIAamp DNA Blood Mini kit (Qiagen, Valencia, CA, USA), according to the manufacturer’s instructions. Purity of DNA was assessed by spectrophotometry using NanoDrop 2000 (Thermo Fisher Scientific, Waltham, MA, USA). DNA yields were determined using a double stranded DNA (dsDNA) broad range (BR) kit on a Qubit® fluorometer (Thermo Fisher Scientific). Integrity of DNA was evaluated by agarose-gel electrophoresis.

2.2. Capture Probe Set Design

Capture probe set was designed and synthesized by Thermo Fisher Scientific for the Ion TargetSeq Custom Enrichment kit (Thermo Fisher Scientific) to target all the coding exons, UTRs, as well as the intronic regions of the NF1 gene. The designed nonoverlapping probe set covered 90.5% of the targeted region (chr17: 29421944-29704695).

2.3. Targeted Capture Library Construction and Next-Generation Sequencing

Libraries were constructed from 2–3 µg of DNA. DNA was sheared into 200 bp fragments on a Bioruptor® sonicator (Diagenode, Denville, NJ, USA), using two cycles of 15 min. Then, DNA fragments were end-repaired, nick repaired, and the barcodes and adaptors were ligated to fragmented DNA. After, we size-selected and amplified the library using the IonXpress Plus gDNA Fragment Library kit, Ion Xpress Barcode Adapters and E-Gel Size Select 2% agarose gels (Thermo Fisher Scientific), according to the manufacturer’s instructions. The amplified library was then hybridized to the NF1 capture probe set for two days (two hybridizations). Libraries were quantified in Qubit® fluorometer (Thermo Fisher Scientific) with dsDNA high sensitivity (HS) kit. The quantity and quality of the library were also verified with Agilent Bioanalyzer DNA High Sensitivity kit (Agilent, Richardson, TX, USA). Libraries were then pooled in equimolar concentrations and diluted to 26 pM, prior to clonal amplification onto Ion Sphere Particles using the Ion One Touch 2 System and the Ion OneTouch 200 Template kit v2 (Thermo Fisher Scientific), following the manufacturer’s instructions. Coated spheres were enriched on the Ion Torrent ES (Thermo Fisher Scientific) before loading onto Ion 316 or 318 sequencing chips (Thermo Fisher Scientific). NGS was performed on Ion Torrent Personal Genome Machine—Ion Torrent PGM™ (Thermo Fisher Scientific).

2.4. Bioinformatic Analysis of Next-Generation Sequencing Data

Primary analyses of NGS data were performed using Torrent Suite software on the Ion Torrent Browser (Thermo Fisher Scientific): base calling, quality score assignment, adapter trimming, and demultiplexing. Reads were generated using Ion Torrent platform pipeline software, which includes standard control of the read quality. Using the same software, reads were aligned against human whole genome (hg19). Only high-quality reads (Q ≥ 20) were aligned to the reference sequence using TMAP (Torrent Mapping Alignment Program).
Germline variants in the targeted regions were detected using the variant caller plug-in of the Ion Torrent Suite software. Variant caller was run with default germline low stringency settings. For each sample, all variants were exported through a vcf (Variant Call Format) file. Variants were filtered following the criteria of a minimum coverage of 30×, and a sample variant frequency between 30%–70%. For variant interpretation, Ingenuity® Variant Analysis™ software (Qiagen) and SNPnexus web server were used. Variants were also viewed manually on Integrative Genomics Viewer (IGV—version 2.3.40; Broad Institute, Cambridge, MA, USA) and on Alamut® Visual (version 2.6.1e; Interactive Biosoftware, Rouen, France) to eliminate strand bias and reduce false positive calls.

2.5. In Silico Predictions

To predict the potential effect of the mutation on the structure and function of the protein codified by the NF1 gene (neurofibromin), PolyPhen2 and SIFT were applied through Ingenuity® Variant Analysis™ software and SNPnexus web server. In silico splicing predictions were investigated with Alamut® Visual, which uses the following prediction methods: SpliceSiteFinder-like; MaxEntScan; NNSPLICE; GeneSplicer; Human Splicing Finder; Known Constitutive Signals; ESEFinder; RESCUE-ESE. Because Alamut® Visual integrates only Human Splicing Finder’s prediction capabilities on donor and acceptor splice sites, and does not investigate the mechanisms of the splicing abnormalities caused by mutations at the nonconsensus splice sites, we also used the online Human Splicing Finder (http://www.umd.be/HSF/).

2.6. Confirmatory Sanger Sequencing

All variants were confirmed by Sanger sequencing in an AB3130 Genetic Analyzer (Applied Biosystems, Carlsbad, CA, USA). Primers were designed using NCBI’s Primer BLAST tool (Table 2) and the sequences were analyzed with the Sequencing Analysis software (version 5.4; Applied Biosystems).

3. Results

3.1. Sequencing Statistics

The average number of bases with Q ≥ 20 across the three Ion Torrent PGMTM runs (12 samples) was 96.9 Mbp, constituting 83.3% of total output. Overall average depth was 177, ranging from 52 to 470. The alignment parameters of the reads for each sample and for each run are shown in Table S1. Average coverage depth for all NF1 exons and introns is shown in Table S2. With the exception of exon 1, which had an average coverage depth of 9.5×, all the other NF1 exons were highly covered, with an average coverage depth ranging from 68.5× to 280.4×. The average coverage depth of the intronic regions ranged from 46.1× to 228.4×. The distribution of reads into different chromosomes is shown in Table 3. About 84% of the reads were mapped to chromosome 17, which contains the NF1 gene.

3.2. Mutation Detection

To detect the pathogenic variants, we blindly inspected all mapped sequence reads from NGS from the 12 samples (11 unrelated patients). Following variant filtering, a total of 3554 variants were identified (≈238/patient; Table S1). With the criteria of minimum coverage of 30×, annotation and filtering of the variants resulted in correct interpretation of the pathogenic variation, as confirmed by Sanger sequencing (Figure 1 and Figure 2), except for one case (NF.19) where no mutation was identified with NGS. In this case, since exon 1 had low coverage, it was sequenced by Sanger method, but no mutation was found (primers are shown in Table 2).
Considering only the 11 unrelated NF1 patients, our approach showed a mutation detection rate of 91% (10 out of 11). We identified eight recurrent and two novel mutations. A summary of the mutations found in this study is shown in Table 4. As expected, the same mutations were found with NGS in the samples from patients NF.84 and NF.85, who were from the same family. Sanger sequencing showed the same mutations in patients NF.14 and NF.84, NF.87 and NF.91, as well as NF.27 and NF.26, who were also from the same family.
Point mutations accounted for 100% of the 10 NF1 pathogenic variants identified in this study. Splicing alterations (c.288 + 5G > A; c.1260 + 1604A > G; c.6483C > G; c.2850 + 1G > A; c.6757 − 2A > T) accounted for 50% (5 out of 10) of the mutations. All were recurrent, with the exception of the mutation c.288 + 5G > A; the splicing effect was predicted by Alamut® Visual. For the mutation c.6483C > G, Alamut® Visual predicted that the reading frame was interrupted by a premature stop codon. As demonstrated by Ars et al. [12], at the mRNA level, this is a nonsense mutation with a splicing effect, which leads to an exon skipping. Therefore, it was classified as a splicing mutation. One of the identified splicing alterations (c.1260 + 1604A > G) was caused by a deep intronic mutation located in intron 11 (1604 bp from exon 10).
Truncating mutations (c.185delT c.4949dupA; c.7867delG) corresponded to 30% (3 out of 10) of the pathogenic variants. The mutation c.185delT was already described by Forzan et al. [13]. This is a mutation that interrupts the reading frame prematurely; codon L62 is replaced by a STOP codon (p.L62*). The mutation c.4949dupA was previously described by Mattocks et al. [14] and van Minkelen et al. [15]. This duplication interrupts the reading frame prematurely: codon Y1650 is replaced by a STOP codon (p.Y1650*). The mutation c.7876delG is novel and, as predicted by Alamut® Visual, creates a frame shift starting at codon A2623. The new reading frame ends in a STOP codon 34 nucleotides downstream (p.A2623Qfs*35). Missense mutations (c.1139T > C, causing a p.L380P leucine to proline amino acid substitution; c.2540T > G, leading to a p.L847R leucine to arginine amino acid substitution) accounted for 20% (2 out of 10) of the pathogenic variants and were all recurrent.

4. Discussion

We developed an approach to analyze the entire NF1 gene, including the exonic and intronic regions by hybridization capture-based NGS. Although some previous studies had evaluated the use of NGS technology to investigate mutations in the NF1 gene [3,7,16,17,18,19,20,21], only exons and conserved splice sites were sequenced, which is the approach commonly used by other investigators based on Sanger methodology.
In our study, the choice of using a hybridization-based targeted selection to prepare the gDNA library, instead of using PCR-based enrichment to extract the NF1 region for sequencing, was based on its ability to capture a large target region in a single experiment, more conveniently, with less hands-on time [22]. We could design a probe set that covered 90.5% of the target bases of the NF1 gene. Most (99.3%) of the bases of the NF1 gene not covered by the designed probe set were located outside the exonic regions. Therefore, there was a small percentage of the exonic regions that was not covered by the capture probe set.
The average of the percentage of the target covered at least once and ≥30× was 97.7% and 96.0%, respectively, showing an excellent capture performance. Also, we were able to reach reads in the 150 bp range with total sequence throughput that guaranteed high coverage depth, which is important to achieve high sensitivity, preventing false-positive variant calls [23]. We identified 10 different pathogenic variations in 11 unrelated patients with clinical diagnosis of NF1, which represents a mutation detection rate of 91%. This high mutation rate is similar to the results of other studies that had used an exhaustive multistep methodology [9,24].
All exons and introns of the NF1 gene had a high depth of coverage with exception of exon 1. As it occurs in other genes, the first exon of the NF1 gene presents higher GC content (71%), compared to the other exons (42%) [7]. It is known that coverage may be reduced in regions with GC contents outside of an optimum range of 40%–60% [25,26], and low coverage in exon 1 of the NF1 gene had been previously reported in other studies with NGS, using both PCR-based enrichment and hybridization-based targeted selection of the NF1 gene for library preparation [7,19]. Regarding the use of PCR-based enrichment, high/low GC content reduces the efficiency of PCR amplification. With respect to the hybridization-based targeted selection, the problem is a reduced efficiency of capture probe hybridization to sequences with high/low GC content [26]. The latter is a property of the capture probes, and therefore with the advance of technologies, it is expected that, to some extent, this limitation is going to be compensated by a better probe set design [26], which would substantially improve the performance of our approach. Sequencing the exon 1 of the NF1 gene is also a challenge in cDNA-based protocols with Sanger methodology [24].
In the only case in which no mutation was found using NGS, we sequenced the coding region of exon 1 by Sanger method. Nevertheless, no pathogenic variation was found. It is possible that the disease-causing mutation was located outside exon 1, in regions of the NF1 gene not covered by the designed probe set or in regions with <30× of coverage depth. As suggested in the literature, in cases of patients with NF1 with no NF1 gene mutation detected, it is possible that the disease is caused by increased methylation of the NF1 promoter region, which would downregulate gene expression [6]. There is also the possibility that the mutation was a single or multiexon deletion or even a total deletion (microdeletion) of the NF1 gene, which are not easily detectable by our approach. Literature reports that ≈5% of NF1 mutations are microdeletions, which have been detected by FISH, MLPA analysis, and/or array comparative genomic hybridization (array-CGH) [2]. Moreover, ≈2.5% of all mutations found in the NF1 gene are small copy-number changes involving one to multiple exons detectable by long-range reverse transcription polymerase chain reaction (RT-PCR) and/or MLPA [2]. Although NGS has the potential to detect copy-number variations, its use for this aim has been limited by a lack of available and effective statistical approaches [29,30]. We believe that in the future, with improvement of bioinformatics tools, NGS will eliminate the necessity for using other methodologies to evaluate these types of mutations in the NF1 gene.
Another challenge of sequencing the NF1 gene is the high number of nonfunctional pseudogenes throughout the human genome. Because pseudogenes are not subject to selective pressure, their mutation rate is much higher than that of the functional genes from which they originate, and therefore their unintended sequencing can lead to false positive variant calls [31]. Target enrichment strategies and data analysis are particularly important to avoid interference of pseudogenes in NGS, allowing reliable detection of true pathogenic variations [31]. In our study, in order to avoid forcing alignment of possible amplified pseudogene regions on the NF1 sequence, reads were aligned against the whole human genome and only high-quality reads (Q ≥ 20) were aligned to the reference sequence. Sanger sequencing confirmed all the mutations identified by NGS, showing that there were no false positive variants from pseudogenes reads.
Naturally, with our approach, the number of total variants that were identified was high (average of 238 variants/patient) compared to other studies that have focused on the exonic regions (average of 30 variants/patient) [3]. Therefore, the use of in silico prediction tools was important to filter out the variants unlikely to be harmful, restricting the search to fewer candidate variants for more detailed manual inspection. Due to the filtering strategy we used, the size of the sequenced region was not an obstacle since only few variants were selected as final candidates.
Homopolymer errors are the dominant error type in Ion Torrent PGMTM data. By manually checking all possible pathogenic variants calls with the IGV software, we were able to identify some as homopolymer artifacts and thus they did not lead to false positive results.
Most (80%) of the mutations identified in this study were recurrent and two (20%) were new. Splicing defects accounted for 40% of the mutations in our sample series. In the literature, approximately 29% of all reported mutations in the NF1 gene cause aberrant splicing and over 57% of them reside outside the canonically conserved GT-AG dinucleotides at the splice donor and acceptor sites [32]. Moreover, around 10% of all splicing defect-causing mutations in the NF1 gene belong to the class of deep intronic mutations (located >500 bp away from exon boundaries), representing 1%–3% of all NF1 mutations [27,32,33,34]. The effect of deep intronic mutations on mRNA is the inclusion of a part of an intron due to the creation of a novel donor or acceptor site or an exonic splicing enhancer site, which, together with a nearby cryptic splice, defines a cryptic exon that is incorporated into the mature mRNA by the spliceosome [35].
Based on the high prevalence of splicing mutations in the NF1 gene, some authors propose the use of mRNA/cDNA for mutation detection, which also identifies deep intronic splicing mutations [32]. Over the past decades, virtually all reported deep intronic mutations in different target genes had been identified through mRNA analysis [10]. Nevertheless, mRNA-based assays present challenges, and precautions need to be taken to prevent problems inherent in this approach [6,10,24]. For example, short-term phytohemagglutinin-stimulated lymphocyte cultures are used to avoid illegitimate splicing that leads to multiple aberrant splice variants that impede the detection of mutations in an RNA-based approach. Moreover, before mRNA extraction, cultures should be treated with puromycin to prevent nonsense-mediated mRNA decay.
The use of gDNA is much easier in clinical settings in comparison to cDNA. Nevertheless, the approach that has been used in the literature to analyze the NF1 gene based on gDNA only sequences the exons and conserved splice sites, which span ≈5 kb, and therefore does not identify deep intronic mutations. Of course, it would be impracticable in clinical settings to use Sanger methodology to sequence the whole NF1 gene, including intronic and exonic regions, which spans ≈282 kb. One great advantage of the NGS technology associated with the in silico prediction tools is the capacity to analyze the entire intronic sequence of a target gene, making possible the identification of deep intronic mutations from a gDNA sample. Nowadays, there are a number of in silico tools that have been developed to assess the effect of DNA sequence variations on splicing, including deep intronic mutations. The online Human Splicing Finder correctly predicted the activation of an exonic cryptic acceptor site in one of our samples (c.1260 + 1604A > G). This deep intronic mutation is recurrent and was reported in three previous studies [15,24,27]. Its consequence had been previously demonstrated by Sabbagh et al. [27] who identified, by a combined cDNA/gDNA sequencing strategy, that this mutation causes a creation of an acceptor splice site, leading to an insertion of a cryptic exon (42 bp out of frame).
Currently, on a routine basis, there is a need to confirm the in silico-predicted splicing alterations [36]. Follow-up analysis of cDNA would be the simplest tool for this. Nevertheless, other techniques can also be used for this purpose using gDNA (i.e., in vitro splicing assay and minigene splicing assay [13,37]). We did not confirm the effect of any of the splicing mutations found in our sample. The majority of them (4/5) were recurrent and the splicing alteration effect had already been previously confirmed by other studies. One of the splicing mutations in our series is novel (c.288 + 5G > A). Mutation at this position with a possible splicing alteration effect (c.288 + 5G > C) had been previously reported [15]. The other novel pathogenic variant identified in our study (c.7867delG), as predicted by Alamut® Visual, was classified as a frameshift mutation with a truncating effect.
Developing a 100% sensitive genetic test for NF1 is necessary to help to establish the diagnosis of patients who do not fulfill the NIH clinical diagnostic criteria [11], which is often the case among young patients with sporadic disease and patients with atypical manifestations or unusual combinations of features [6,16]. Moreover, identification of the pathogenic NF1 mutation provides a presymptomatic/prenatal test in the offspring of sporadic patients [24]. In the research field, NGS will be extremely important to accelerate the rate of identification of mutations in the NF1 gene, therefore facilitating the study of the effect of specific mutations on structure and function of the gene product, and leading to a better understanding of the correlation between genotype and phenotype.
Limitations of this study include the small sample size and also the existence of small regions of the NF1 gene not covered by the designed probe set. Nevertheless, our results, together with the results of others that used NGS to sequence the NF1 gene [3,16,17,18,19,20], show that this technology is applicable to NF1 mutation screening. In this study, we show the results of a simple and fast approach, using hybridization capture-based NGS with gDNA as starting material, to screen, at once, the entire NF1 gene (exons and introns) for different types of pathogenic variations, including the deep intronic splicing mutations. More studies with young patients with sporadic disease and patients with atypical features would be important to demonstrate the clinical utility of our approach.

Supplementary Materials

The following are available online at www.mdpi.com/2073-4425/7/12/133/s1, Table S1: Alignment parameters of the reads, Table S2: Mean coverage of NF1 exons and introns among 12 patients.

Acknowledgments

We thank Ludwine Messiaen for her review and useful comments. We thank the individuals and families for their participation in the study. This research was supported by L’Oréal-UNESCO-Brazilian Academy of Sciences grant for Young Women Scientists (to K.S.C.), and by FAPERJ (Fundação de Amparo à Pesquisa do Estado do Rio de Janeiro, Brazil), under number E-26/110811/2013 (to K.S.C).

Author Contributions

K.S.C. conceived and designed the study. K.S.C., C.C.S. and M.G. obtained patient material. K.S.C., C.C.S., N.S.O., A.K.F., A.G., T.B., Y.I., R.S., J.-P.M. and D.C. performed experiments/analyzed and interpreted data. KSC drafted and wrote the manuscript. A.G., J.-P.M., R.S.M.N., R.S., M.G., and D.C. critically revised the manuscript. All authors read and approved the final manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Uusitalo, E.; Leppävirta, J.; Koffert, A.; Suominen, S.; Vahtera, J.; Vahlberg, T.; Pöyhönen, M.; Peltonen, J.; Peltonen, S. Incidence and mortality of neurofibromatosis: A total population study in Finland. J. Investig. Dermatol. 2015, 135, 904–906. [Google Scholar] [CrossRef] [PubMed]
  2. Cunha, K.S.G.; Geller, M. Advances in Neurofibromatosis Research; Nova Science Publishers Inc.: Hauppauge, NY, USA, 2011. [Google Scholar]
  3. Pasmant, E.; Parfait, B.; Luscan, A.; Goussard, P.; Briand-Suleau, A.; Laurendeau, I.; Fouveaut, C.; Leroy, C.; Montadert, A.; Wolkenstein, P.; et al. Neurofibromatosis type 1 molecular diagnosis: What can NGS do for you when you have a large gene with loss of function mutations? Eur. J. Hum. Genet. EJHG 2015, 23, 596–601. [Google Scholar] [CrossRef] [PubMed]
  4. Clementi, M.; Barbujani, G.; Turolla, L.; Tenconi, R. Neurofibromatosis-1: A maximum likelihood estimation of mutation rate. Hum. Genet. 1990, 84, 116–118. [Google Scholar] [CrossRef] [PubMed]
  5. De Luca, A.; Schirinzi, A.; Buccino, A.; Bottillo, I.; Sinibaldi, L.; Torrente, I.; Ciavarella, A.; Dottorini, T.; Porciello, R.; Giustini, S.; et al. Novel and recurrent mutations in the NF1 gene in Italian patients with neurofibromatosis type 1. Hum. Mutat. 2004, 23, 629. [Google Scholar] [CrossRef] [PubMed]
  6. Griffiths, S.; Thompson, P.; Frayling, I.; Upadhyaya, M. Molecular diagnosis of neurofibromatosis type 1: 2 years experience. Fam. Cancer 2007, 6, 21–34. [Google Scholar] [CrossRef] [PubMed]
  7. Uusitalo, E.; Hammais, A.; Palonen, E.; Brandt, A.; Mäkelä, V.-V.; Kallionpää, R.; Jouhilahti, E.-M.; Pöyhönen, M.; Soini, J.; Peltonen, J.; et al. Neurofibromatosis type 1 gene mutation analysis using sequence capture and high-throughput sequencing. Acta Derm. Venereol. 2014, 94, 663–666. [Google Scholar] [CrossRef] [PubMed]
  8. Upadhyaya, M. The Molecular Biology of Neurofibromatosis Type 1; Morgan & Claypool: Williston, VT, USA, 2014. [Google Scholar]
  9. Messiaen, L.M.; Callens, T.; Mortier, G.; Beysen, D.; Vandenbroucke, I.; Van Roy, N.; Speleman, F.; Paepe, A.D. 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, 541–555. [Google Scholar] [CrossRef]
  10. Flanagan, S.E.; Xie, W.; Caswell, R.; Damhuis, A.; Vianey-Saban, C.; Akcay, T.; Darendeliler, F.; Bas, F.; Guven, A.; Siklar, Z.; et al. Next-Generation Sequencing Reveals Deep Intronic Cryptic ABCC8 and HADH Splicing Founder Mutations Causing Hyperinsulinism by Pseudoexon Activation. Am. J. Hum. Genet. 2013, 92, 131–136. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. Stumpf, D.A.; Alksne, J.F.; Annegers, J.F.; Brown, S.S.; Conneally, P.M.; Housman, D.; Leppert, M.F.; Miller, J.P.; Moss, M.L.; Pileggi, A.J.; et al. Neurofibromatosis. Conference statement. National Institutes of Health Consensus Development Conference. Arch Neurol. 1988, 45, 575–578. [Google Scholar]
  12. Ars, E.; Kruyer, H.; Morell, M.; Pros, E.; Serra, E.; Ravella, A.; Estivill, X.; Lázaro, C. Recurrent mutations in the NF1 gene are common among neurofibromatosis type 1 patients. J. Med. Genet. 2003, 40, e82. [Google Scholar] [CrossRef] [PubMed]
  13. Forzan, M. Functional Splicing Assay mediante l’utilizzo di minigeni plasmidici nel gene NF1. Ph.D. Thesis, Università degli Studi di Padova, Padova, Italy, 2010. [Google Scholar]
  14. Mattocks, C.; Baralle, D.; Tarpey, P. Automated comparative sequence analysis identifies mutations in 89% of NF1 patients and confirms a mutation cluster in exons 11–17 distinct from the GAP related domain. Br. Med. J. 2004, 41, e48. [Google Scholar] [CrossRef]
  15. Van Minkelen, R.; van Bever, Y.; Kromosoeto, J.N.R.; Withagen-Hermans, C.J.; Nieuwlaat, A.; Halley, D.J.J.; van den Ouweland, A.M.W. A clinical and genetic overview of 18 years neurofibromatosis type 1 molecular diagnostics in the Netherlands. Clin. Genet. 2014, 85, 318–327. [Google Scholar] [CrossRef] [PubMed]
  16. Balla, B.; Árvai, K.; Horváth, P.; Tobiás, B.; Takács, I.; Nagy, Z.; Dank, M.; Fekete, G.; Kósa, J.P.; Lakatos, P. Fast and robust next-generation sequencing technique using ion torrent personal genome machine for the screening of neurofibromatosis type 1 (NF1) gene. J. Mol. Neurosci. MN 2014, 53, 204–210. [Google Scholar] [CrossRef] [PubMed]
  17. Chou, L.-S.; Liu, C.-S.J.; Boese, B.; Zhang, X.; Mao, R. DNA sequence capture and enrichment by microarray followed by next-generation sequencing for targeted resequencing: Neurofibromatosis type 1 gene as a model. Clin. Chem. 2010, 56, 62–72. [Google Scholar] [CrossRef] [PubMed]
  18. Hutter, S.; Piro, R.M.; Waszak, S.M.; Kehrer-Sawatzki, H.; Friedrich, R.E.; Lassaletta, A.; Witt, O.; Korbel, J.O.; Lichter, P.; Schuhmann, M.U.; et al. No correlation between NF1 mutation position and risk of optic pathway glioma in 77 unrelated NF1 patients. Hum. Genet. 2016, 135, 469–475. [Google Scholar] [CrossRef] [PubMed]
  19. Maruoka, R.; Takenouchi, T.; Torii, C.; Shimizu, A.; Misu, K.; Higasa, K.; Matsuda, F.; Ota, A.; Tanito, K.; Kuramochi, A.; et al. The use of next-generation sequencing in molecular diagnosis of neurofibromatosis type 1: A validation study. Genet. Test. Mol. Biomark. 2014, 18, 722–735. [Google Scholar] [CrossRef] [PubMed]
  20. Banerjee, S.; Dai, Y.; Liang, S.; Chen, H.; Wang, Y.; Tang, L.; Wu, J.; Huang, H. A novel mutation in NF1 is associated with diverse intra-familial phenotypic variation and astrocytoma in a Chinese family. J. Clin. Neurosci. Off. J. Neurosurg. Soc. Australas. 2016, 31, 182–184. [Google Scholar] [CrossRef] [PubMed]
  21. Calì, F.; Chiavetta, V.; Ruggeri, G.; Piccione, M.; Selicorni, A.; Palazzo, D.; Bonsignore, M.; Cereda, A.; Elia, M.; Failla, P.; et al. Mutation spectrum of NF1 gene in Italian patients with neurofibromatosis type 1 using Ion Torrent PGMTM platform. Eur. J. Med. Genet. 2016. [Google Scholar] [CrossRef] [PubMed]
  22. Mamanova, L.; Coffey, A.J.; Scott, C.E.; Kozarewa, I.; Turner, E.H.; Kumar, A.; Howard, E.; Shendure, J.; Turner, D.J. Target-enrichment strategies for next-generation sequencing. Nat. Methods 2010, 7, 111–118. [Google Scholar] [CrossRef] [PubMed]
  23. Galindo-González, L.; Pinzón-Latorre, D.; Bergen, E.A.; Jensen, D.C.; Deyholos, M.K. Ion Torrent sequencing as a tool for mutation discovery in the flax (Linum usitatissimum L.) genome. Plant Methods 2015, 11, 19. [Google Scholar] [CrossRef] [PubMed]
  24. Valero, M.C.; Martín, Y.; Hernández-Imaz, E.; Marina Hernández, A.; Meleán, G.; Valero, A.M.; Javier Rodríguez-Álvarez, F.; Tellería, D.; Hernández-Chico, C. A highly sensitive genetic protocol to detect NF1 mutations. J. Mol. Diagn. JMD 2011, 13, 113–122. [Google Scholar] [CrossRef] [PubMed]
  25. Summerer, D.; Schracke, N.; Wu, H.; Cheng, Y.; Bau, S.; Stähler, C.F.; Stähler, P.F.; Beier, M. Targeted high throughput sequencing of a cancer-related exome subset by specific sequence capture with a fully automated microarray platform. Genomics 2010, 95, 241–246. [Google Scholar] [CrossRef] [PubMed]
  26. Chilamakuri, C.S.R.; Lorenz, S.; Madoui, M.-A.; Vodák, D.; Sun, J.; Hovig, E.; Myklebost, O.; Meza-Zepeda, L.A. Performance comparison of four exome capture systems for deep sequencing. BMC Genom. 2014, 15. [Google Scholar] [CrossRef] [PubMed]
  27. Sabbagh, A.; Pasmant, E.; Imbard, A.; Luscan, A.; Soares, M.; Blanché, H.; Laurendeau, I.; Ferkal, S.; Vidaud, M.; Pinson, S.; et al. NF1 molecular characterization and neurofibromatosis type I genotype-phenotype correlation: The French experience. Hum. Mutat. 2013, 34, 1510–1518. [Google Scholar] [CrossRef] [PubMed]
  28. Ars, E.; Serra, E.; Garcia, J.; Kruyer, H.; Gaona, A.; Lazaro, C.; Estivill, X. Mutations affecting mRNA splicing are the most common molecular defects in patients with neurofibromatosis type 1. Hum. Mol. Genet. 2000, 9, 237. [Google Scholar] [CrossRef] [PubMed]
  29. Wang, H.; Nettleton, D.; Ying, K. Copy number variation detection using next generation sequencing read counts. BMC Bioinform. 2014, 15, 109. [Google Scholar] [CrossRef] [PubMed]
  30. Abou Tayoun, A.N.; Tunkey, C.D.; Pugh, T.J.; Ross, T.; Shah, M.; Lee, C.C.; Harkins, T.T.; Wells, W.A.; Tafe, L.J.; Amos, C.I.; et al. A comprehensive assay for CFTR mutational analysis using next-generation sequencing. Clin. Chem. 2013, 59, 1481–1488. [Google Scholar] [CrossRef] [PubMed]
  31. Poliseno, L. (Ed.) Pseudogenes; Methods in Molecular Biology Series; Springer: New York, NY, USA, 2014; Volume 1167.
  32. Messiaen, L.; Wimmer, K. Mutation analysis of the NF1 gene by cDNA-based sequencing of the coding region. In Advances in Neurofibromatosis Research; Cunha, K.S.G., Geller, M., Eds.; Nova Science Publishers Inc.: New York, NY, USA, 2012; pp. 89–108. [Google Scholar]
  33. Wimmer, K.; Roca, X.; Beiglböck, H.; Callens, T.; Etzler, J.; Rao, A.R.; Krainer, A.R.; Fonatsch, C.; Messiaen, L. Extensive in silico analysis of NF1 splicing defects uncovers determinants for splicing outcome upon 5′ splice-site disruption. Hum. Mutat. 2007, 28, 599–612. [Google Scholar] [CrossRef] [PubMed]
  34. Svaasand, E.K.; Engebretsen, L.F.; Ludvigsen, T.; Brechan, W.; Sjursen, W. A Novel Deep Intronic Mutation Introducing a Cryptic Exon Causing Neurofibromatosis Type 1 in a Family with Highly Variable Phenotypes: A Case Study. Hered. Genet. Curr. Res. 2015, 4, 3. [Google Scholar] [CrossRef]
  35. Castellanos, E.; Rosas, I.; Solanes, A.; Bielsa, I.; Lázaro, C.; Carrato, C.; Hostalot, C.; Prades, P.; Roca-Ribas, F.; Blanco, I.; et al. In vitro antisense therapeutics for a deep intronic mutation causing Neurofibromatosis type 2. Eur. J. Hum. Genet. 2013, 21, 769–773. [Google Scholar] [CrossRef] [PubMed]
  36. Coutinho, G.; Xie, J.; Du, L.; Brusco, A.; Krainer, A.R.; Gatti, R.A. Functional significance of a deep intronic mutation in the ATM gene and evidence for an alternative exon 28a. Hum. Mutat. 2005, 25, 118–124. [Google Scholar] [CrossRef] [PubMed]
  37. Jian, X.; Boerwinkle, E.; Liu, X. In silico tools for splicing defect prediction: A survey from the viewpoint of end users. Genet. Med. 2014, 16, 497–503. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Mutations found using next-generation sequencing (NGS) and illustrated in Integrative Genomics Viewer (IGV). Inserts: Sanger sequencing of the corresponding mutations: (A) Patient 1—c.7867delG; (B) Patient 2—c.1139T > C; (C) Patient 3—c.185delT; (D) Patient 4—c.288 + 5G > A; (E) Patient 5—c.2850 + 1G > A; (F) Patient 6—c.6757 − 2A > T.
Figure 1. Mutations found using next-generation sequencing (NGS) and illustrated in Integrative Genomics Viewer (IGV). Inserts: Sanger sequencing of the corresponding mutations: (A) Patient 1—c.7867delG; (B) Patient 2—c.1139T > C; (C) Patient 3—c.185delT; (D) Patient 4—c.288 + 5G > A; (E) Patient 5—c.2850 + 1G > A; (F) Patient 6—c.6757 − 2A > T.
Genes 07 00133 g001
Figure 2. Mutations found using NGS and illustrated in IGV. Inserts: Sanger sequencing of the corresponding mutations: (A) Patient 7—c.1260 + 1604A > G; (B) Patient 8—c.6483C > G; (C) Patient 9—c.2540T > G; (D) Patient 10—c.4949dupA.
Figure 2. Mutations found using NGS and illustrated in IGV. Inserts: Sanger sequencing of the corresponding mutations: (A) Patient 7—c.1260 + 1604A > G; (B) Patient 8—c.6483C > G; (C) Patient 9—c.2540T > G; (D) Patient 10—c.4949dupA.
Genes 07 00133 g002
Table 1. Summary of some clinical features of the participants.
Table 1. Summary of some clinical features of the participants.
Patient NumberGenderAgeFamily HistoryNumber of Skin NeurofibromasSuperficial Plexiform NeurofibromaNumber of Café-au-Lait MaculesInguinal/Axillary Freckling
NF.01Female46No200–499Yes3Yes
NF.14Female54Yes500–1000Yes20Yes
NF.19Female26No500–1000No35Yes
NF.26Male52Yes>1000No3Yes
NF.27Female55Yes>1000No4Yes
NF.42Male47No200–499No9Yes
NF.54Female32No200–499No12Yes
NF.63Female48Yes500–1000No12Yes
NF.78Female13Yes<10No44Yes
NF.83Male48No>1000Yes8Yes
NF.84Female45Yes200–499Yes28Yes
NF.85Male24Yes200–499Yes48Yes
NF.87Male34Yes200–499No6Yes
NF.91Male72No>1000No2Yes
NF.96Male41Yes100–199No17No
Table 2. Primers used for PCR amplification and sequencing of the neurofibromatosis (NF1) gene by the Sanger method.
Table 2. Primers used for PCR amplification and sequencing of the neurofibromatosis (NF1) gene by the Sanger method.
Patient NumberLocationSequences (5′–3′)TmAmplicon
NF.01Intron 11F: TAATGAGCCAGGGCATTGTACC66 °C385 pb
R: CTTTCACCAAGTACACTGAGGC66 °C
NF.19 *Exon 1F: CACAGACCCTCTCCTTGCCTCTTC71 °C243 pb
R: TACCTCCCCTCACCTACTCTGTCC68 °C
NF.26
NF.27
Exon 43F: TGCTGTTTGGCATTAGCAAAGT62 °C302 pb
R: TGTTACCAATAACACAGTCCATGC68 °C
NF.42Exon 37F: ATACCGGGCCTAGCAATCGC64 °C133 pb
R: TTGGTGTACTCCCTGACCCAGG70 °C
NF.54Exon 2F: AGCAGAACACACATACCAAAGTCAG72 °C158 pb
R: AATTCCCCAAAACACAGTAACCCAA72 °C
NF.63Exon 3/Intron 3F: GATGTGTGTTGATTGGTAGCAGA66 °C245 pb
R: GGACTGTCCTCTTGGTCCACA66 °C
NF.78Exon 54F: CTTGGCAGGCTACACTGGT60 °C158 pb
R: ACTTAAAGACAGGCACGAAGGT64 °C
NF.83Exon 21F: AAGAAATTTGACACTCGGCTGAT64 °C483 pb
R: TGCTGACAGGTGTATCTGCG62 °C
NF.85
NF.84
NF.14
Exon 10F: AGCTGGATTTTACTGCCATTTGTG68 °C233 pb
R: TAAAGTGTTGGTTGTTGTGAGGG66 °C
NF.91
NF.87
Exon 21/Intron 21F: CCTGCTCTGTATCCAATGCTAT64 °C133 pb
R: GCTTATTTCAAACAAGTCACTCT62 °C
NF.96Intron 45/Exon 46F: AGCTAGCTACCAAGATCACCA62 °C297 pb
R: ACACTGATACCCAAAATGAATGC64 °C
F: forward; R: reverse; Tm: primer’s melting temperature (2 × (AT) + 4 × (GC)); bp: size of the amplicons in base pairs. * In this case, no pathogenic variant was identified using next-generation sequencing (NGS). Since exon 1 showed low coverage, it was sequenced using Sanger method, but no mutation was found.
Table 3. Percentage of reads mapped to chromosome 17, which contains the NF1 gene, and to chromosomes that contain NF1 pseudogenes.
Table 3. Percentage of reads mapped to chromosome 17, which contains the NF1 gene, and to chromosomes that contain NF1 pseudogenes.
ChromosomePercentage (%)
Chr21.72
Chr120.84
Chr141.27
Chr152.57
Chr1783.94
Chr181.09
Chr211.06
Chr221.81
Other chromosomes5.71
Table 4. Positions of the mutations found in the study.
Table 4. Positions of the mutations found in the study.
Patient NumberExon/Intron (NG_009018.1)DNA Mutation gDNA Level (NG_009018.1)DNA Mutation cDNA Level (NM_000267.3)Predicted Protein (NP_000258.1)Mutation TypeVariant EffectReferences
NF.01IVS11g.29530107A > Gc.1260 + 1604A > Gp.N420_S421insLTT* (also noted as p.S421LfsX4)SubstitutionSplicing (deep intronic mutation)Valero et al. [24];
Sabbagh et al. [27];
van Minkelen et al. [15]
NF.26
NF.27
E43g.29664504C > Gc.6483C > Gp.Y2161*SubstitutionSplicing bArs et al. [28]
NF.42E37g.29653014dupc.4949dupAp.Y1650*DuplicationFrameshift truncationMattocks et al. [14];
van Minkelen et al. [15]
NF.54E2g.29483125delc.185delTp.L62*DeletionFrameshift truncationForzan [13]
NF.63IVS3g.29486116G > Ac.288 + 5G > Ap.(?)SubstitutionSplicingND
NF.78E54g.29684347delc.7867delGp.A2623Qfs*35DeletionFrameshift truncationND
NF.83E21g.29556173T > Gc.2540T > Gp.L847RSubstitutionMissensevan Minkelen et al. [15];
De Luca et al. [5];
Cali et al. [21]
NF.85
NF.84
NF.14
E10g.29528131T > Cc.1139T > Cp.L380PSubstitutionMissensevan Minkelen et al. [15]
NF.91
NF.87
IVS21g.29556484G > Ac.2850 + 1G > Ap.(?)SubstitutionSplicingArs et al. [28]
van Minkelen et al. [15]
NF.96IVS45g.29665720A > Tc.6757 − 2A > Tp.(?)SubstitutionSplicingPasmant et al. [3]
p.(?): Effect on protein level not known with certainty; ND: mutation not described until now; the mutations described by van Minkelen et al. [15] have been deposited in Leiden Open Variation Database (LOVD); a Nonsense mutation, as predicted by Alamut® Visual. It has a splicing alteration effect, causing an exon skipping, demonstrated at mRNA level by Ars et al. [28].

Share and Cite

MDPI and ACS Style

Cunha, K.S.; Oliveira, N.S.; Fausto, A.K.; De Souza, C.C.; Gros, A.; Bandres, T.; Idrissi, Y.; Merlio, J.-P.; De Moura Neto, R.S.; Silva, R.; et al. Hybridization Capture-Based Next-Generation Sequencing to Evaluate Coding Sequence and Deep Intronic Mutations in the NF1 Gene. Genes 2016, 7, 133. https://doi.org/10.3390/genes7120133

AMA Style

Cunha KS, Oliveira NS, Fausto AK, De Souza CC, Gros A, Bandres T, Idrissi Y, Merlio J-P, De Moura Neto RS, Silva R, et al. Hybridization Capture-Based Next-Generation Sequencing to Evaluate Coding Sequence and Deep Intronic Mutations in the NF1 Gene. Genes. 2016; 7(12):133. https://doi.org/10.3390/genes7120133

Chicago/Turabian Style

Cunha, Karin Soares, Nathalia Silva Oliveira, Anna Karoline Fausto, Carolina Cruz De Souza, Audrey Gros, Thomas Bandres, Yamina Idrissi, Jean-Philippe Merlio, Rodrigo Soares De Moura Neto, Rosane Silva, and et al. 2016. "Hybridization Capture-Based Next-Generation Sequencing to Evaluate Coding Sequence and Deep Intronic Mutations in the NF1 Gene" Genes 7, no. 12: 133. https://doi.org/10.3390/genes7120133

APA Style

Cunha, K. S., Oliveira, N. S., Fausto, A. K., De Souza, C. C., Gros, A., Bandres, T., Idrissi, Y., Merlio, J. -P., De Moura Neto, R. S., Silva, R., Geller, M., & Cappellen, D. (2016). Hybridization Capture-Based Next-Generation Sequencing to Evaluate Coding Sequence and Deep Intronic Mutations in the NF1 Gene. Genes, 7(12), 133. https://doi.org/10.3390/genes7120133

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop