“Genotype-first” approaches on a curious case of idiopathic progressive cognitive decline
© Shi et al.; licensee BioMed Central Ltd. 2014
Received: 10 July 2014
Accepted: 20 November 2014
Published: 3 December 2014
In developing countries, many cases with rare neurological diseases remain undiagnosed due to limited diagnostic experience. We encountered a case in China where two siblings both began to develop idiopathic progressive cognitive decline starting from age six, and were suspected to have an undiagnosed neurological disease.
Initial clinical assessments included review of medical history, comprehensive physical examination, genetic testing for metabolic diseases, blood tests and brain imaging. We performed exome sequencing with Agilent SureSelect exon capture and Illumina HiSeq2000 platform, followed by variant annotation and selection of rare, shared mutations that fit a recessive model of inheritance. To assess functional impacts of candidate variants, we performed extensive biochemical tests in blood and urine, and examined their possible roles by protein structure modeling.
Exome sequencing identified NAGLU as the most likely candidate gene with compound heterozygous mutations (chr17:40695717C > T and chr17:40693129A > G in hg19 coordinate), which were documented to be pathogenic. Sanger sequencing confirmed the recessive patterns of inheritance, leading to a genetic diagnosis of Sanfilippo syndrome (mucopolysaccharidosis IIIB). Biochemical tests confirmed the complete loss of activity of alpha-N-acetylglucosaminidase (encoded by NAGLU) in blood, as well as significantly elevated dermatan sulfate and heparan sulfate in urine. Structure modeling revealed the mechanism on how the two variants affect protein structural stability.
Successful diagnosis of a rare genetic disorder with an atypical phenotypic presentation confirmed that such “genotype-first” approaches can particularly succeed in areas of the world with insufficient medical genetics expertise and with cost-prohibitive in-depth phenotyping.
Next-generation sequencing (NGS) technologies have advanced many aspects of genomic sciences, and have accelerated the pace of discovery in biomedical research -. For example, NGS has been widely used in population genetics studies ,, metagenomics ,, agrigenomics ,, epigenetics  and gene expression studies ,. In particular, human genome and exome sequencing can now be used in dissecting the genetic basis of diseases and traits that have proven to be intractable to conventional gene-discovery strategies -. In addition, NGS empowers clinical diagnostics and other aspects of medical practice, such as genetic diagnosis, disease prognosis, therapeutic target identification, optimization of treatment regimens, prenatal testing and personalized disease-risk profiling.
Besides identifying new mutations in genes previously implicated in known diseases, genome/exome sequencing has also been used to characterize “idiopathic” or “mysterious” diseases, which present phenotypes that may be due to genetic causes. The vast majority of such diseases are indeed “known” diseases that cannot be easily diagnosed by conventional candidate gene approaches, or require complicated differential diagnosis. For example, a recent study from the Baylor Whole Genome Laboratory reported a success rate of 25 percent to provide a genetic diagnosis for 250 patients, who presented a range of phenotypes suggesting potential genetic causes and ~80% were children with neurologic phenotypes . An exome sequencing study on 100 patients with intellectual disability found that 53 of them carried de novo mutations, including 13 in known intellectual-disability genes, suggesting that exome sequencing is an effective diagnostic approach for diseases with extensive clinical heterogeneity . Additionally, several studies have reported the discovery of completely novel diseases through genome/exome sequencing. One of the early examples was Ogden syndrome, a previously unreported infantile lethal disorder, involving a mutation in NAA10 identified by us in 2011 using chromosome exon X capture and next-generation sequencing . A more recent example is a novel genetic disease which we refer to as “Bookman syndrome”, a pediatric onset disease with neuromuscular and cardiac involvement and with clinical features similar to Glycogen Storage Disease Type IV. Although exome sequencing failed to identify the causes for the disease due to technical reasons, we applied genome and transcriptome sequencing and identified a disease-contributory mutation in RBCK1 , which was further replicated by another group . Besides finding causal genes for genetic diseases, genome/exome sequencing may also aid in the treatment of rare diseases. One of the earliest high-profile examples is the Pulitzer Prize-winning story of Nicholas Volker, a young boy who at the age of 17 months experienced perineal fistulae coupled with leakage of stool via holes spanning from the intestine to the skin . Exome sequencing on this boy resulted in a diagnosis of an X-linked inhibitor of apoptosis deficiency. Based on this finding, an allogeneic hematopoietic progenitor cell transplant was performed on Volker. Altogether, these examples demonstrated the power of genome/exome sequencing in uncovering the genetic basis of idiopathic diseases and finding treatments for them, even when only a very small number of patient samples are available.
Among people with such diseases, those with idiopathic neuropathy present special challenges for genetic diagnosis ,. Due to the heterogeneous nature of clinical presentation, the difficulty with characterizing phenotypes precisely, and the locus heterogeneity for the same diseases, many patients with inherited neuropathy may not obtain a genetic diagnosis by candidate gene testing. However, these patients may be more likely to receive a positive diagnosis by exome sequencing interrogating all genes in a somewhat unbiased manner, and multiple recent studies have demonstrated successful examples ,. Frequently, after a candidate gene is identified, it is then found that an atypical phenotype accounts for earlier failed candidate gene approaches, so that genome/exome sequencing can indeed expand the clinical spectrum of previously reported or novel mutations in known disease genes, especially for patients with neuropathy . Perhaps more importantly, genome/exome sequencing can nominate important candidate genes or candidate diseases, which guides the selection of functional assays to confirm the diagnosis of the disease.
The two probands in our study are siblings from a family in Guangzhou, Guangdong Province in China. The healthy parents are non-consanguineous, and had one affected son and one affected daughter (Figure 1). No evidence of any similar genetic disease was reported in family members. The study to identify disease genes through exome sequencing was approved by Jinan University First Affiliated Hospital Institutional Review Board (reference number: 2013–004). An informed consent for participation of study and publication of genetic results as well as individual medical details was obtained from each family member (or their guardians) who participated in the study.
Individual II-2 is a 9 years old girl. She was born at term, also with normal birth parameters. She began to stand at 11 months, walk with aid at 13 months, and speak at 17 months. At age 5, she was just like other children of similar age, with the ability to dress and sing, and count by herself. Starting at 6 years of age, she began to show regression of developmental patterns: she could not dress by herself anymore, and could not express even a single sentence or count numbers. Clinical examination revealed a coarse face with low anterior and posterior hairlines, prominent frontal bossing, thick eyebrows, synophrys (unibrow), hypertelorism, and thick lips. Growth parameters were normal. Her clinical course was also severe, with progressive neurodegeneration, behavioral problems (including hyperactivity, impulsivity, obstinacy, anxious behaviors and autistic-like behaviors), and hearing loss. The diagnosis of severe intellectual disability was made, based on Wechsler Intelligence Scale examination. Measuring activities of daily living showed extreme disability. Brain MRI demonstrated cortical atrophy with enlargement of the subarachnoid spaces and ventricular dilatation (Figure 2). Brainstem evoked potentials showed moderate abnormalities. EEG recording showed abnormal sleep EEG, just like her brother’s manifestation.
We elected to perform exome sequencing on the probands, and to use parental samples for Sanger validation. Genomic DNA was extracted from peripheral blood for all family members by Qiagen DNeasy Blood & Tissue kit (Valencia, CA, USA). We obtained at least 5ug DNA from each sample, and all DNA samples passed quality control measures for exome sequencing.
Exome capture was carried out on 1 μg of genomic DNA, using the Agilent SureSelectXT Human All Exon 50 Mb kit (Agilent Technologies, Santa Clara, CA, USA), which targets 50 Mb of coding exonic sequences, which were annotated by the GENCODE project, CCDS, miRbaseV14 and Sanger, as well as 10 base pairs of flanking sequence of each exonic region and non-coding RNA. The captured and amplified library was then loaded onto the Illumina Hiseq2000 sequencer. We generated 100 bp paired-end reads in the sequencing run, according to the manufacturer’s protocols. Raw images were processed by Illumina Pipeline v1.3.4 for base-calling with default parameters. Illumina SCS and CASAVA software were used for raw data processing and FASTQ file generation. In total, we obtained ~70 million reads per sample.
Sequence alignment and variant calling
For bioinformatics analysis of the FASTQ data from exome sequencing, we used SeqMule (http://seqmule.usc.edu), which is an automated pipeline to execute multiple alignment algorithms and multiple variant calling software tools. Briefly, sequencing data was evaluated with FastQC (http://www.bioinformatics.babraham.ac.uk/projects/fastqc/). Short reads were aligned to reference genome (hg19) by BWA-MEM (version 0.7.4)  algorithm with default settings. Then we used three variant calling software tools, including the GATK version 3.1 , SAMtools version 0.1.19  and FreeBayes version 0.9.14  with default settings for exome sequencing. We selected the consensus variants using the “two of three algorithms” protocol, which specifies that variants that were detected by at least two algorithms will be present in the final consensus calling set.
Variant annotation and prioritization
We used the ANNOVAR software  for functional annotation of variants and the Tute web application (tutegenomics.com) for identifying candidate variants and genes for the clinical phenotype. We used a “single-case analysis” pipeline to identify a list of candidate genes with the following criteria: (1) identify variants causing splicing or protein coding changes, including stop loss and stop gain variants; (2) remove variants with minor allele frequency (MAF) greater than 1% in the PopFreqMax database from ANNOVAR package , which include maximum allele frequency from the 1000 Genomes Project April 2012 release (five ethnicity groups), the NHLBI-6500 Exomes (two ethnicity groups), and the Complete Genomics 46 genomes (CG46) database; (3) imposed a recessive mode of inheritance, with at least two deleterious mutations found in each candidate genes in the proband, and with the mutations shared by both probands; (4) The resulting variants and genes were sorted based on relevance to user-supplied phenotypes, as well as annotations from clinical databases such as ClinVar (version 20140303) and HGMD (version 2014_2) .
Validation by Sanger sequencing
Selected putative variants were examined among all family members using Sanger sequencing. Given the chromosomal position of variants, we designed PCR primers to amplify fragments harboring individual variants by Primer3 . The PCR primers were designed to encompass the candidate position, ensuring that common SNPs are not covered by the primers. The ABI 3730 XL sequencer was used for sequencing, and the resulting *.AB1 files were loaded into the ABI Sequence Scanner Software v1.0 for further analysis and genotype calling. All sequence traces were manually reviewed to ensure the reliability of the genotype calls.
The primers used for chr17 :40695717(C > T) mutation are CCCGCCTCTTCCCCAACTC (forward) and GGACGCCTCCAGCCCTCAA (reverse). The primers used for chr17 :40693129 (A > G) mutation are AAACCAGGAGCTGTAGAGAAGT (forward) and CTGCCTACCCCTACTGACATCT (reverse).
We searched the Protein Data Bank (PDB)  for structure models of NAGLU. Although the structure for human NAGLU was not available, a previous study solved the crystal structure, catalytic mechanism, and inhibition of CpGH89 from Clostridium perfringens, a close bacterial homolog of NAGLU . We downloaded the structure from X-ray diffraction with a resolution of 2.36 Å with binding of the ligand 2-Acetamindo-1,2-Dideoxynojirmycin (PDB identifier: 2VC9). We built a structure model for human NAGLU using I-TASSER , given high sequence identity. The template structure 2VC9 binds to 2-ACETAMIDO-1,2-DIDEOXYNOJIRMYCIN (PDB identifier: NOK). The position of NOK is decided by superpositioning our model to the template structure 2VC9 (shown in ball and stick). The 3D representation of the structure was generated by PYMOL (http://www.pymol.org/).
Biochemical test in blood/urine
Enzyme activity was measured on peripheral blood leukocytes using fluorogenic substrate obtained from Moscerdam Substrates (Rotterdam, Netherlands); peripheral blood leukocytes was collected and sonicated to protein homogenates. Protein homogenates were incubated with Moscerdam substrate. The parameters for enzyme activity was obtained by measuring the fluorescence intensity and comparing with standard fluorescence intensity and concentration. The assays were performed according to the protocol at Department of Endocrinology and Metabolism, Guangzhou Women and Children’s Medical Center -. Urinary GAG was measured using dimethylmethylene blue/Tris by spectrophotometry and corrected for urinary creatinine (Cr) content ,. The ratio of GAG/Cr (mg/mmol) was compared to age-matched normal controls (6 m-12 m, < 23.3; 12 m-36 m, <17.3; 3-5 years, <13.2; >5 years, < 9.7).
Assessment of clinical information
The two probands were brother and sister one year apart in age, and both of them had perfectly normal development patterns in the first five years of age. Starting from age six, the probands suddenly began to develop progressive speech and cognitive decline with behavioral difficulties, motor function decline and hearing loss: they gradually lost the ability to talk, count, or walk. Rough facial features began to develop including synophrys (unibrow), thick lips and enlarged mouth. Brain MRI demonstrated cortical atrophy with enlargement of the subarachnoid spaces and ventricular dilatation (Figure 2). An extended list of clinical features is given in Materials and Methods. In the course of the clinical workup, multiple attempts for genetic diagnosis were made, including tests of chromosomal abnormalities and sequencing on a panel for metabolic diseases. Multiple types of blood screening tests and physical examinations were performed on the pedigree, especially the probands, to exclude the common inherited metabolic diseases and nervous system abnormalities. Additional factors were considered to exclude the possibility of food poisoning, environmental toxification and other non-genetic reasons. However, none of these attempts showed positive results, and these two probands were eventually referred to as having “mental deterioration of unknown origin”.
To research the causes of the diseases in this family, we examined all available phenotype information on the probands. Both parents are healthy and are non-consanguineous. No evidence of other genetic disease was reported in any additional family members. We suspected that this is an extremely rare or novel genetic disease, possibly following a recessive pattern of inheritance. Exome sequencing appeared to be the next natural choice to identify candidate mutations possibly contributing to this idiopathic disease.
Exome sequencing identifies a prioritized list of candidate mutations
Bioinformatics predictions of deleteriousness on the two non-synonymous variants in NAGLU
Mutation Taster score
Mutation Taster prediction
Disease causing automatic
Mutation Assessor score
Mutation Assessor prediction
Sanger validation of the candidate variants
Functional impacts of the candidate variants by biochemical test
Results of biochemical tests on the activity of MPS-related enzymes
heparan sulfate sulfatase
galactosamine 6- sulfate Sulfatase
Results from urinary test of glycosaminoglycan related parameters
Urinary glycosaminoglycan (GAG)
Urinary glycosaminoglycan/Creatinine (GAG/Cr)
Functional impacts of the candidate variants by structural modeling
To further understand the functional impacts of the two variants, we attempted to examine structure models for NAGLU. A previous study solved the crystal structure, catalytic mechanism, and inhibition of CpGH89 from Clostridium perfringens, a close bacterial homolog of NAGLU . The structure was generated by X-ray diffraction with a resolution of 2.36 Å. We built a structure model for human NAGLU based on this template (PDB identifier: 2VC9) using I-TASSER , given the high sequence identity. The complete structure contains four domains (Figure 4C), including a N-terminal domain which is a putative family 32 carbohydrate-binding module (CBM) with the typical β-sandwich fold, a catalytic region comprises a small α/β domain, an elaborated (α/β)8 barrel, and an all α-helical domain that packs against the first three domains. We found that these two mutations were not directly within the active site (ligand-binding site) of NAGLU. Mutation Y309C is in the middle of the elaborated alpha/beta barrel (in blue) and R565W is located at the alpha-helical domain (in red) that packs against the binding site. Therefore, mutation Y309C may affect the structural stability of the alpha/beta barrel that constitutes the bottom of the active site. Mutation R565W may plays an important role in structural packing of the alpha-helical domain that contributes to the stability of the complete structure.
In this study, we present a case illustrating the power of NGS in clinical diagnosis of a rare disease in two siblings with somewhat atypical phenotypes, with a disease onset at ~ six years of age. Although this turned out to be a known disease with a known genetic cause, our experience in this study is nevertheless revealing in the context of a study conducted in China, where traditionally experienced medical geneticists are lacking. We suspected that had a more experienced medical geneticist reviewed the patient, biochemical tests on different types of MPS could have been ordered, and a diagnosis of MPS-IIIB may have been made without exome sequencing. However, in practice, numerous other tests were ordered and all turned out to be negative, to a point where the family was reluctant to pay for any additional tests. After exome sequencing, the mutation pointed us to a phenotype related to MPS, and therefore we performed biochemical tests in blood and urine to confirm this diagnosis. Therefore, exome sequencing in this particular context provided guidance for us to reach a final diagnosis, and this also illustrates the importance of supporting bioinformatics predictions with functional evidence.
Sanfilippo syndrome B (or MPS-IIIB) is an extremely rare disorder and the prevalence in China is unreported and unknown. Due to the rare nature of the disease, the incidence rate could vary significantly among different parts of the world. Using multiple ascertainment sources, one study obtained an incidence rate for Sanfilippo syndrome (form A + B + C) in western Australia for the period 1969 to 1996 of approximately 1 in 58,000 live births . Since there were 5 people with type B in the total of 11 cases, this suggested an incidence of MPS-IIIB of approximately 1 in 127,600 live births. In different studies, an incidence of 1:24,000 was reported for all subtypes of MPS III . A more recent study on the Australian population estimated an incidence rate of 1 in 200,000 . Northern Ireland tends to have much lower incidence of MPS, and a report showed the rate of 1 in 280,000 for all subtypes of MPS III . In the Netherland, the combined birth prevalence for all subtypes of MPS-III is 1.89 per 100,000 live births .
While we emphasize that the development of technologies has now enabled rapid and relatively cost-efficient identification of genetic causes in all parts of the world, the major bottleneck is still the identification of disease-contributory variants, which is heavily dependent on the bioinformatics analysis of raw data and functional interpretation of genetic variants. We used the SeqMule pipeline for variant calling, to take ‘2 out of 3’ consensus calls from three variant calling algorithms, in a way similar to what the 1000 Genomes Project has been implementing in their variant calling procedure. We also used ANNOVAR and the Tute web server for variant analysis focusing on missense, nonsense and splice variants. For the two missense variants in NAGLU, we used multiple prediction algorithms, including six from dbNSFP  (SIFT, PolyPhen, LRT, Mutation Taster, Mutation Assessor, FATHMM), two meta-scores built by us (MetaSVM, MetaLR), and the Tute score (tutegenomics.com), and found that these algorithms gave consistent predictions. However, the main reason we focused on NAGLU is that it is reported in clinical databases (ClinVar and HGMD) as containing mutations linked to disease. This fact underscores the importance of curated mutation databases in helping identify disease causal mutations from newly sequenced genomes or exomes.
In conclusion, we have presented an example in China where exome sequencing preceded and guided the selection of functional validations, to lead to a diagnosis of a known disease with a somewhat atypical presentation. With the ever decreasing cost of NGS and the rapid adoption of NGS in clinical settings, we expect that genome/exome sequencing will be increasingly useful in the diagnosis of rare genetic disorders with atypical phenotypic presentations, including in areas of the world where experienced medical geneticists are scarce.
The authors thank the probands and their family members for donating blood samples for scientific research, and for providing detailed phenotype information and medical records. We thank Dr. Feng Wang (Wuhan Institute of Technology and Kaigene Technology) for help with data analysis and archival. L.S. is supported by the Fundamental Research Funds for the Central Universities of China (21614331), the Medical Science Foundation of Guangdong Province (B2014213) and the National Natural Science Foundation of China (31400922). K.W. is supported by NIH/NHGRI grant number R01 HG006465. G.J.L is supported by funds from the Stanley Institute for Cognitive Genomics at Cold Spring Harbor Laboratory.
- Koboldt DC, Steinberg KM, Larson DE, Wilson RK, Mardis ER: The next-generation sequencing revolution and its impact on genomics. Cell. 2013, 155 (1): 27-38. 10.1016/j.cell.2013.09.006.View ArticlePubMedPubMed CentralGoogle Scholar
- Shendure J, Ji H: Next-generation DNA sequencing. Nat Biotechnol. 2008, 26 (10): 1135-1145. 10.1038/nbt1486.View ArticlePubMedGoogle Scholar
- Metzker ML: Sequencing technologies - the next generation. Nat Rev Genet. 2010, 11 (1): 31-46. 10.1038/nrg2626.View ArticlePubMedGoogle Scholar
- Durbin RM, Abecasis GR, Altshuler DL, Auton A, Brooks LD, Gibbs RA, Hurles ME, McVean GA: A map of human genome variation from population-scale sequencing. Nature. 2010, 467 (7319): 1061-1073. 10.1038/nature09534.View ArticlePubMedGoogle Scholar
- Pool JE, Hellmann I, Jensen JD, Nielsen R: Population genetic inference from genomic sequence variation. Genome Res. 2010, 20 (3): 291-300. 10.1101/gr.079509.108.View ArticlePubMedPubMed CentralGoogle Scholar
- Pallen MJ, Loman NJ, Penn CW: High-throughput sequencing and clinical microbiology: progress, opportunities and challenges. Curr Opin Microbiol. 2010, 13 (5): 625-631. 10.1016/j.mib.2010.08.003.View ArticlePubMedGoogle Scholar
- Schloss PD, Handelsman J: Metagenomics for studying unculturable microorganisms: cutting the Gordian knot. Genome Biol. 2005, 6 (8): 229-10.1186/gb-2005-6-8-229.View ArticlePubMedPubMed CentralGoogle Scholar
- Poland JA, Brown PJ, Sorrells ME, Jannink JL: Development of high-density genetic maps for barley and wheat using a novel two-enzyme genotyping-by-sequencing approach. PLoS One. 2012, 7 (2): e32253-10.1371/journal.pone.0032253.View ArticlePubMedPubMed CentralGoogle Scholar
- Goddard ME, Hayes BJ: Mapping genes for complex traits in domestic animals and their use in breeding programmes. Nat Rev Genet. 2009, 10 (6): 381-391. 10.1038/nrg2575.View ArticlePubMedGoogle Scholar
- Meaburn E, Schulz R: Next generation sequencing in epigenetics: insights and challenges. Semin Cell Dev Biol. 2012, 23 (2): 192-199. 10.1016/j.semcdb.2011.10.010.View ArticlePubMedGoogle Scholar
- Wang Z, Gerstein M, Snyder M: RNA-Seq: a revolutionary tool for transcriptomics. Nat Rev Genet. 2009, 10 (1): 57-63. 10.1038/nrg2484.View ArticlePubMedPubMed CentralGoogle Scholar
- Ozsolak F, Milos PM: RNA sequencing: advances, challenges and opportunities. Nat Rev Genet. 2011, 12 (2): 87-98. 10.1038/nrg2934.View ArticlePubMedGoogle Scholar
- Bamshad MJ, Ng SB, Bigham AW, Tabor HK, Emond MJ, Nickerson DA, Shendure J: Exome sequencing as a tool for Mendelian disease gene discovery. Nat Rev Genet. 2011, 12 (11): 745-755. 10.1038/nrg3031.View ArticlePubMedGoogle Scholar
- Rabbani B, Tekin M, Mahdieh N: The promise of whole-exome sequencing in medical genetics. J Hum Genet. 2014, 59 (1): 5-15. 10.1038/jhg.2013.114.View ArticlePubMedGoogle Scholar
- Ku CS, Naidoo N, Pawitan Y: Revisiting Mendelian disorders through exome sequencing. Hum Genet. 2011, 129 (4): 351-370. 10.1007/s00439-011-0964-2.View ArticlePubMedGoogle Scholar
- Gilissen C, Hoischen A, Brunner HG, Veltman JA: Unlocking Mendelian disease using exome sequencing. Genome Biol. 2011, 12 (9): 228-10.1186/gb-2011-12-9-228.View ArticlePubMedPubMed CentralGoogle Scholar
- Yang Y, Muzny DM, Reid JG, Bainbridge MN, Willis A, Ward PA, Braxton A, Beuten J, Xia F, Niu Z, Hardison M, Person R, Bekheirnia MR, Leduc MS, Kirby A, Pham P, Scull J, Wang M, Ding Y, Plon SE, Lupski JR, Beaudet AL, Gibbs RA, Eng CM: Clinical whole-exome sequencing for the diagnosis of mendelian disorders. N Engl J Med. 2013, 369 (16): 1502-1511. 10.1056/NEJMoa1306555.View ArticlePubMedPubMed CentralGoogle Scholar
- de Ligt J, Willemsen MH, van Bon BW, Kleefstra T, Yntema HG, Kroes T, Vulto-van Silfhout AT, Koolen DA, de Vries P, Gilissen C, del Rosario M, Hoischen A, Scheffer H, de Vries BB, Brunner HG, Veltman JA, Vissers LE: Diagnostic exome sequencing in persons with severe intellectual disability. N Engl J Med. 2012, 367 (20): 1921-1929. 10.1056/NEJMoa1206524.View ArticlePubMedGoogle Scholar
- Rope AF, Wang K, Evjenth R, Xing J, Johnston JJ, Swensen JJ, Johnson WE, Moore B, Huff CD, Bird LM, Carey JC, Opitz JM, Stevens CA, Jiang T, Schank C, Fain HD, Robison R, Dalley B, Chin S, South ST, Pysher TJ, Jorde LB, Hakonarson H, Lillehaug JR, Biesecker LG, Yandell M, Arnesen T, Lyon GJ: Using VAAST to identify an X-linked disorder resulting in lethality in male infants due to N-terminal acetyltransferase deficiency. Am J Hum Genet. 2011, 89 (1): 28-43. 10.1016/j.ajhg.2011.05.017.View ArticlePubMedPubMed CentralGoogle Scholar
- Wang K, Kim C, Bradfield J, Guo Y, Toskala E, Otieno FG, Hou C, Thomas K, Cardinale C, Lyon GJ, Golhar R, Hakonarson H: Whole-genome DNA/RNA sequencing identifies truncating mutations in RBCK1 in a novel Mendelian disease with neuromuscular and cardiac involvement. Genome Med. 2013, 5 (7): 67-10.1186/gm471.View ArticlePubMedPubMed CentralGoogle Scholar
- Nilsson J, Schoser B, Laforet P, Kalev O, Lindberg C, Romero NB, Dávila López M, Akman HO, Wahbi K, Iglseder S, Eggers C, Engel AG, Dimauro S, Oldfors A: Polyglucosan body myopathy caused by defective ubiquitin ligase RBCK1. Ann Neurol. 2013, 74 (6): 914-919. 10.1002/ana.23963.View ArticlePubMedGoogle Scholar
- Worthey EA, Mayer AN, Syverson GD, Helbling D, Bonacci BB, Decker B, Serpe JM, Dasu T, Tschannen MR, Veith RL, Basehore MJ, Broeckel U, Tomita-Mitchell A, Arca MJ, Casper JT, Margolis DA, Bick DP, Hessner MJ, Routes JM, Verbsky JW, Jacob HJ, Dimmock DP: Making a definitive diagnosis: successful clinical application of whole exome sequencing in a child with intractable inflammatory bowel disease. Genet Med. 2011, 13 (3): 255-262. 10.1097/GIM.0b013e3182088158.View ArticlePubMedGoogle Scholar
- Lyon GJ, O’Rawe J: Human genetics and clinical aspects of neurodevelopmental disorders. The Genetics of Neurodevelopmental Disorders. Edited by: Mitchell K. 2014, Wiley, Hoboken, NJ, USAGoogle Scholar
- Foo JN, Liu JJ, Tan EK: Whole-genome and whole-exome sequencing in neurological diseases. Nat Rev Neurol. 2012, 8 (9): 508-517. 10.1038/nrneurol.2012.148.View ArticlePubMedGoogle Scholar
- Zhao H, Race V, Matthijs G, De Jonghe P, Robberecht W, Lambrechts D, Van Damme P: Exome sequencing reveals HINT1 mutations as a cause of distal hereditary motor neuropathy. Eur J Hum Genet. 2014, 22 (6): 847-850. 10.1038/ejhg.2013.231.View ArticlePubMedGoogle Scholar
- Beetz C, Pieber TR, Hertel N, Schabhüttl M, Fischer C, Trajanoski S, Graf E, Keiner S, Kurth I, Wieland T, Varga RE, Timmerman V, Reilly MM, Strom TM, Auer-Grumbach M: Exome sequencing identifies a REEP1 mutation involved in distal hereditary motor neuropathy type V. Am J Hum Genet. 2012, 91 (1): 139-145. 10.1016/j.ajhg.2012.05.007.View ArticlePubMedPubMed CentralGoogle Scholar
- Klein CJ, Middha S, Duan X, Wu Y, Litchy WJ, Gu W, Dyck PJ, Gavrilova RH, Smith DI, Kocher JP, Dyck PJ: Application of whole exome sequencing in undiagnosed inherited polyneuropathies. J Neurol Neurosurg Psychiatry. 2014, 85 (11): 1265-1272. 10.1136/jnnp-2013-306740.View ArticlePubMedGoogle Scholar
- Stessman HA, Bernier R, Eichler EE: A genotype-first approach to defining the subtypes of a complex disease. Cell. 2014, 156 (5): 872-877. 10.1016/j.cell.2014.02.002.View ArticlePubMedPubMed CentralGoogle Scholar
- Lyon GJ, Wang K: Identifying disease mutations in genomic medicine settings: current challenges and how to accelerate progress. Genome Med. 2012, 4 (7): 58-10.1186/gm359.View ArticlePubMedPubMed CentralGoogle Scholar
- Li H: Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM. arXiv. 2013, 1303 (3997v2): 1-3.Google Scholar
- McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Kernytsky A, Garimella K, Altshuler D, Gabriel S, Daly M, DePristo MA: The genome analysis toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010, 20 (9): 1297-1303. 10.1101/gr.107524.110.View ArticlePubMedPubMed CentralGoogle Scholar
- Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R: The Sequence Alignment/Map format and SAMtools. Bioinformatics. 2009, 25 (16): 2078-2079. 10.1093/bioinformatics/btp352.View ArticlePubMedPubMed CentralGoogle Scholar
- Garrison E, Marth G: Haplotype-based variant detection from short-read sequencing. arXiv Preprint. 2012, arXiv: 1207-3907. arXivGoogle Scholar
- Wang K, Li M, Hakonarson H: ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Res. 2010, 38 (16): e164-10.1093/nar/gkq603.View ArticlePubMedPubMed CentralGoogle Scholar
- Landrum MJ, Lee JM, Riley GR, Jang W, Rubinstein WS, Church DM, Maglott DR: ClinVar: public archive of relationships among sequence variation and human phenotype. Nucleic Acids Res. 2014, 42 (Database issue): D980-D985. 10.1093/nar/gkt1113.View ArticlePubMedGoogle Scholar
- Untergasser A, Cutcutache I, Koressaar T, Ye J, Faircloth BC, Remm M, Rozen SG: Primer3–new capabilities and interfaces. Nucleic Acids Res. 2012, 40 (15): e115-10.1093/nar/gks596.View ArticlePubMedPubMed CentralGoogle Scholar
- Rose PW, Bi C, Bluhm WF, Christie CH, Dimitropoulos D, Dutta S, Green RK, Goodsell DS, Prlic A, Quesada M, Quinn GB, Ramos AG, Westbrook JD, Young J, Zardecki C, Berman HM, Bourne PE: The RCSB protein data bank: new resources for research and education. Nucleic Acids Res. 2013, 41 (Database issue): D475-D482. 10.1093/nar/gks1200.View ArticlePubMedGoogle Scholar
- Ficko-Blean E, Stubbs KA, Nemirovsky O, Vocadlo DJ, Boraston AB: Structural and mechanistic insight into the basis of mucopolysaccharidosis IIIB. Proc Natl Acad Sci U S A. 2008, 105 (18): 6560-6565. 10.1073/pnas.0711491105.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhang Y: I-TASSER server for protein 3D structure prediction. BMC Bioinform. 2008, 9: 40-10.1186/1471-2105-9-40.View ArticleGoogle Scholar
- Voznyi YV, Keulemans JL, van Diggelen OP: A fluorimetric enzyme assay for the diagnosis of MPS II (Hunter disease). J Inherit Metab Dis. 2001, 24 (6): 675-680. 10.1023/A:1012763026526.View ArticlePubMedGoogle Scholar
- Karpova EA, Voznyi Ya V, Keulemans JL, Hoogeveen AT, Winchester B, Tsvetkova IV, van Diggelen OP: A fluorimetric enzyme assay for the diagnosis of Sanfilippo disease type A (MPS IIIA). J Inherit Metab Dis. 1996, 19 (3): 278-285. 10.1007/BF01799255.View ArticlePubMedGoogle Scholar
- Marsh J, Fensom AH: 4-Methylumbelliferyl alpha-N-acetylglucosaminidase activity for diagnosis of Sanfilippo B disease. Clin Genet. 1985, 27 (3): 258-262. 10.1111/j.1399-0004.1985.tb00217.x.View ArticlePubMedGoogle Scholar
- van Diggelen OP, Zhao H, Kleijer WJ, Janse HC, Poorthuis BJ, van Pelt J, Kamerling JP, Galjaard H: A fluorimetric enzyme assay for the diagnosis of Morquio disease type A (MPS IV A). Clin Chim Acta. 1990, 187 (2): 131-139. 10.1016/0009-8981(90)90339-T.View ArticlePubMedGoogle Scholar
- Huang YL, Li SY, Zhao XY, Fan LP, Lin WC, Zhou ZH, Cheng J, Liu L: Enzymatic diagnosis and clinical characteristics of 52 children with mucopolysaccharidosis. Zhongguo Dang Dai Er Ke Za Zhi. 2012, 14 (7): 510-514.PubMedGoogle Scholar
- de Jong JG, Wevers RA, Liebrand-van Sambeek R: Measuring urinary glycosaminoglycans in the presence of protein: an improved screening procedure for mucopolysaccharidoses based on dimethylmethylene blue. Clin Chem. 1992, 38 (6): 803-807.PubMedGoogle Scholar
- Zhao XY, Huang YL, Li SY, Lin WC, Zhou ZH, Liu L: Quantitative measuring urinary glycosaminoglycan by dimethylmethylene-tris spectrophotometric method. Chin J Child Health Care. 2010, 11: 885-888.Google Scholar
- O'Rawe J, Jiang T, Sun G, Wu Y, Wang W, Hu J, Bodily P, Tian L, Hakonarson H, Johnson WE, Wei Z, Wang K, Lyon GJ: Low concordance of multiple variant-calling pipelines: practical implications for exome and genome sequencing. Genome Med. 2013, 5 (3): 28-10.1186/gm432.View ArticlePubMedPubMed CentralGoogle Scholar
- Liu X, Jian X, Boerwinkle E: dbNSFP: a lightweight database of human nonsynonymous SNPs and their functional predictions. Hum Mutat. 2011, 32 (8): 894-899. 10.1002/humu.21517.View ArticlePubMedPubMed CentralGoogle Scholar
- Schmidtchen A, Greenberg D, Zhao HG, Li HH, Huang Y, Tieu P, Zhao HZ, Cheng S, Zhao Z, Whitley CB, Di Natale P, Neufeld EF: NAGLU mutations underlying Sanfilippo syndrome type B. Am J Hum Genet. 1998, 62 (1): 64-69. 10.1086/301685.View ArticlePubMedPubMed CentralGoogle Scholar
- Andria G, Di Natale P, Del Giudice E, Strisciuglio P, Murino P: Sanfilippo B syndrome (MPS III B): mild and severe forms within the same sibship. Clin Genet. 1979, 15 (6): 500-504. 10.1111/j.1399-0004.1979.tb00832.x.View ArticlePubMedGoogle Scholar
- Chinen Y, Tohma T, Izumikawa Y, Uehara H, Ohta T: Sanfilippo type B syndrome: five patients with an R565P homozygous mutation in the alpha-N-acetylglucosaminidase gene from the Okinawa islands in Japan. J Hum Genet. 2005, 50 (7): 357-359. 10.1007/s10038-005-0258-4.View ArticlePubMedGoogle Scholar
- van de Kamp JJ, Niermeijer MF, von Figura K, Giesberts MA: Genetic heterogeneity and clinical variability in the Sanfilippo syndrome (types A, B, and C). Clin Genet. 1981, 20 (2): 152-160. 10.1111/j.1399-0004.1981.tb01821.x.View ArticlePubMedGoogle Scholar
- Stenson PD, Mort M, Ball EV, Shaw K, Phillips A, Cooper DN: The Human Gene Mutation Database: building a comprehensive mutation repository for clinical and molecular genetics, diagnostic testing and personalized genomic medicine. Hum Genet. 2014, 133 (1): 1-9. 10.1007/s00439-013-1358-4.View ArticlePubMedGoogle Scholar
- Weber B, Guo XH, Kleijer WJ, van de Kamp JJ, Poorthuis BJ, Hopwood JJ: Sanfilippo type B syndrome (mucopolysaccharidosis III B): allelic heterogeneity corresponds to the wide spectrum of clinical phenotypes. Eur J Hum Genet. 1999, 7 (1): 34-44. 10.1038/sj.ejhg.5200242.View ArticlePubMedGoogle Scholar
- Lee-Chen GJ, Lin SP, Lin SZ, Chuang CK, Hsiao KT, Huang CF, Lien WC: Identification and characterisation of mutations underlying Sanfilippo syndrome type B (mucopolysaccharidosis type IIIB). J Med Genet. 2002, 39 (2): E3-10.1136/jmg.39.2.e3.View ArticlePubMedPubMed CentralGoogle Scholar
- Nelson J, Crowhurst J, Carey B, Greed L: Incidence of the mucopolysaccharidoses in Western Australia. Am J Med Genet A. 2003, 123A (3): 310-313. 10.1002/ajmg.a.20314.View ArticlePubMedGoogle Scholar
- Meikle PJ, Hopwood JJ, Clague AE, Carey WF: Prevalence of lysosomal storage disorders. JAMA. 1999, 281 (3): 249-254. 10.1001/jama.281.3.249.View ArticlePubMedGoogle Scholar
- Nelson J: Incidence of the mucopolysaccharidoses in Northern Ireland. Hum Genet. 1997, 101 (3): 355-358. 10.1007/s004390050641.View ArticlePubMedGoogle Scholar
- Poorthuis BJ, Wevers RA, Kleijer WJ, Groener JE, de Jong JG, van Weely S, Niezen-Koning KE, van Diggelen OP: The frequency of lysosomal storage diseases in The Netherlands. Hum Genet. 1999, 105 (1–2): 151-156. 10.1007/s004399900075.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.