A, Brain magnetic resonance imaging studies showed microcephaly, a simplified gyral pattern, and normal pons and cerebellum at ages 4 years and 1 year in both affected siblings BAB3022 and BAB3280, respectively. B, Profound nonprogressive microcephaly was documented for all affected individuals. Upper panels show fronto-occipital circumference (FOC) measurements for both patients. Solid black lines indicate the mean FOC for age; dashed lines, ±2-SD FOC values for age; red and blue lines, corresponding FOC trend lines for each patient (red, BAB3280; blue, BAB3022); and dots, individual measurements at specific ages of the respective patients. Lower panels show standard deviations (SDs) below the mean for head circumference measurements in the 2 patients. C, Electron microscopy of a sural nerve biopsy specimen demonstrated the loss of large myelinated fibers, increased endoneural collagen, and rare degenerating myelin (original magnification ×1500). D, Hematoxylin-eosin staining of a muscle biopsy specimen showed neurogenic muscle atrophy with spinal muscular atrophy–like atrophy of large fascicles alternating with well-preserved nonatrophic fascicles (original magnification ×200). E, Family pedigree shows segregation of disease-causing VRK1 mutations, p.V236M and p.R89Q, with the neurological phenotype; patients BAB3280 and BAB3022 are compound heterozygotes for p.V236M and p.R89Q, their father (BAB3024) is a carrier for p.R89Q, their mother (BAB3023) is a carrier of p.V236M, and their unaffected sister (BAB3693) has no VRK1 mutant alleles.
A, In patient BAB5311, brain magnetic resonance imaging revealed microcephaly, a simplified gyral pattern, normal pons and cerebellar hemispheres, and vermian hypoplasia. B, He showed nonprogressive microcephaly, which was first noted by ultrasonography in utero at 30 weeks’ gestation and was recorded after birth. Upper panel shows fronto-occipital circumference (FOC) measurements. Solid black line without dots indicates the mean FOC for age; dashed lines, ±2-SD FOC values for age; line with dots, corresponding FOC trend line for the patient; and dots, individual measurements at specific ages of the patient. Lower panel shows standard deviations (SDs) below the mean for head circumference measurements. C, Family pedigree shows segregation of the nonsense p.R358X mutation; the affected proband (BAB5311) is homozygous, while both his parents (BAB5321 and BAB5320) and sister (BAB5322) are heterozygous carriers.
The data show that both unaffected parents of the patient are heterozygous, while the patient is homozygous for the common haplotype region in which the p.R358X mutation occurred. The region, which is approximately 500 673 base pairs, includes VRK1 only. In the upper plots, 1 pixel indicates 3 kilobases; in the lower plots, 1 pixel indicates 330 base pairs.
Gonzaga-Jauregui C, Lotze T, Jamal L, Penney S, Campbell IM, Pehlivan D, Hunter JV, Woodbury SL, Raymond G, Adesina AM, Jhangiani SN, Reid JG, Muzny DM, Boerwinkle E, Lupski JR, Gibbs RA, Wiszniewski W. Mutations in VRK1 Associated With Complex Motor and Sensory Axonal Neuropathy Plus Microcephaly. JAMA Neurol. 2013;70(12):1491-1498. doi:10.1001/jamaneurol.2013.4598
Patients with rare diseases and complex clinical presentations represent a challenge for clinical diagnostics. Genomic approaches are allowing the identification of novel variants in genes for very rare disorders, enabling a molecular diagnosis. Genomics is also revealing a phenotypic expansion whereby the full spectrum of clinical expression conveyed by mutant alleles at a locus can be better appreciated.
To elucidate the molecular cause of a complex neuropathy phenotype in 3 patients by applying genomic sequencing strategies.
Design, Setting, and Participants
Three affected individuals from 2 unrelated families presented with a complex neuropathy phenotype characterized by axonal sensorimotor neuropathy and microcephaly. They were recruited into the Centers for Mendelian Genomics research program to identify the molecular cause of their phenotype. Whole-genome, targeted whole-exome sequencing, and high-resolution single-nucleotide polymorphism arrays were performed in genetics clinics of tertiary care pediatric hospitals and biomedical research institutions.
Main Outcomes and Measures
Whole-genome and whole-exome sequencing identified the variants responsible for the patients’ clinical phenotype.
We identified compound heterozygous alleles in 2 affected siblings from 1 family and a homozygous nonsense variant in the third unrelated patient in the vaccinia-related kinase 1 gene (VRK1). In the latter subject, we found a common haplotype on which the nonsense mutation occurred and that segregates in the Ashkenazi Jewish population.
Conclusions and Relevance
We report the identification of disease-causing alleles in 3 children from 2 unrelated families with a previously uncharacterized complex axonal motor and sensory neuropathy accompanied by severe nonprogressive microcephaly and cerebral dysgenesis. Our data raise the question of whether VRK1 mutations disturb cell cycle progression and may result in apoptosis of cells in the nervous system. The application of unbiased genomic approaches allows the identification of potentially pathogenic mutations in unsuspected genes in highly genetically heterogeneous and uncharacterized neurological diseases.
Hereditary motor and sensory neuropathies (HMSNs) represent a group of slowly progressive neurological diseases caused by dysfunction of the peripheral nerves with secondary muscle wasting and weakness usually presenting as a distal symmetric polyneuropathy.1 Hereditary motor and sensory neuropathies are genetically heterogeneous, with more than 40 disease-causing genes identified to date. Patients most often present with isolated nerve disease, but on rare occasions neuropathy is accompanied by additional clinical signs and symptoms yielding a complex neuropathy. The development of genome-wide analysis methods has enabled the identification of novel HMSN-associated genes and has revealed that some patients with atypical peripheral neuropathy have mutations in known HMSN genes.
We studied 3 patients with a severe, rapidly progressing distal symmetric polyneuropathy who also had microcephaly. Using genomic sequencing approaches, we identified likely pathogenic variants in the 3 affected individuals from 2 different families in the vaccinia-related kinase 1 gene (VRK1; GenBank NM_003384).
This research project was approved by the Baylor College of Medicine Institutional Review Board. Written informed consent was obtained from the parents of the subjects.
Two sisters (patients BAB3022 and BAB3280), aged 10 and 2 years, respectively, and a 9-year-old boy (patient BAB5311) from an unrelated family presented with delayed motor development and microcephaly with normal cognition. Patient BAB3022 was born after an uncomplicated full-term pregnancy; delayed motor development and microcephaly were noted on the 18-month well-child care visit. She began rolling over at 12 months, sat at 24 months, and was never able to walk without support. She became wheelchair bound at age 4 years, and 2 years later she developed limited use of her arms and hands. The patient experienced worsening scoliosis, compromising her respiratory function and requiring surgical intervention at age 10 years. Due to worsening respiratory function and decreased oral intake, tracheostomy and gastrostomy placement were implemented. She attends regular school and her cognitive testing did not reveal any intellectual deficits. Physical evaluations revealed nonprogressive microcephaly, hypotonia, progressive muscle wasting, diffuse weakness progressing from distal to proximal muscle groups, and loss of deep tendon reflexes. The patient did not have hypertonia, abnormal movements, or nystagmus.
Extensive clinical evaluation included brain magnetic resonance imaging (MRI) at ages 1 and 4 years, showing microcephaly with a simplified gyral pattern but no evidence of other intracranial abnormalities (Figure 1A). Electrophysiological studies showed evidence of severe axonal motor and sensory peripheral neuropathy with severely decreased amplitudes; some sensory and motor nerves were nonresponsive (Table). Axonal etiology was further supported by histopathological studies of the sural nerve that showed mild to moderate reduction of the number of large myelinated fibers and a reduced number of unmyelinated axons (Figure 1C). Muscle biopsy showed marked variation in fiber size and shape and large-group atrophy with increased connective tissue that was suggestive of neurogenic myopathy (Figure 1D). The patient underwent extensive genetic laboratory evaluation, and findings were normal on metabolic, cytogenetic, and molecular studies. Despite numerous efforts, the cause of the neurological problems remained unexplained.
The affected sibling, patient BAB3280, presented with microcephaly in utero (Figure 1B). Her motor development was delayed; she was unable to sit or walk without support at age 20 months. Intellectual development was concordant with expected milestones. She manifested mild hypotonia but had normal muscle bulk, antigravity or greater power in all extremities, and normal deep tendon reflexes. No abnormal movements were appreciated during physical examination. Brain MRI showed microcephaly and a simplified gyral pattern but no other cranial abnormalities (Figure 1A). The electrophysiological studies found evidence for axonal motor and sensory neuropathy as in her older sibling (Table). Family history was negative for neurological disorders, and the parents denied consanguinity. The affected girls have another sister who at age 5 years is normocephalic with normal motor and intellectual development and normal findings on physical examination. The presence of clinical symptoms in 2 affected children suggested a genetic cause with a potential autosomal recessive mode of inheritance.
Patient BAB5311 was a fraternal twin, conceived by in vitro fertilization and born at 34.5 weeks of pregnancy. Prenatal surveillance ultrasonography revealed microcephaly at 30 weeks’ gestation and decreased fetal movements. The patient began to noticeably deviate from his developmental milestones at age 4 months. He began to sit at 8 months, walked at 18 months, and remained ambulatory until 6 years. Clinical evaluation at age 9 years revealed severe nonprogressive microcephaly (6 SDs below mean; Figure 2B), worsening hypotonia, muscle atrophy, tongue atrophy, decreased deep tendon reflexes, and preserved sensation and cognition. Medical history was also positive for sleeping problems, tremor, hypophonia, and dysarthric speech. Similar to patient BAB3022, the patient has had worsening scoliosis and required gastric feeding tube placement. Brain MRI showed microcephaly, simplified gyral pattern, normal pons and cerebellar hemispheres, and underdeveloped cerebellar vermis (Figure 2A). The electrophysiological studies revealed motor and sensory axonal neuropathy (Table). Neither the parents nor the twin sister have any evidence of neurological problems based on self-report of lack of neuropathy symptoms and objective clinical examination (Figure 2C).
The clinical evaluation was performed by several independent neurologists and medical geneticists. The MRI scans of all affected individuals were evaluated by a pediatric neuroradiologist. Electron microscopy of a sural nerve biopsy specimen and hematoxylin-eosin staining of a muscle biopsy specimen of patient BAB3022 were evaluated by a neuropathologist. Neurophysiological studies were carried out in the affected individuals; motor nerve conduction velocities and sensory nerve conduction velocities were determined.
Patient BAB3022 underwent whole-genome sequencing on the SOLiD sequencing platform (Life Technologies Corp) at approximately 30 times the average depth of coverage. Mapping to the reference human genome sequence and variant calling were performed using the SOLiD software tools and Corona Lite bioinformatic suite. Patient BAB3280 underwent whole-exome capture using the Baylor College of Medicine Human Genome Sequencing Center VCRome version 2.1 design (Roche NimbleGen, Inc) followed by sequencing on the HiSeq platform (Illumina, Inc) at approximately 100 times the average depth of coverage through the Centers for Mendelian Genomics research initiative. Sequence data were aligned and mapped to the human genome reference sequence (hg19) using the Mercury in-house bioinformatic pipeline.2 Variants were called using the ATLAS3 and SAMtools4 suites and annotated with an in-house–developed annotation pipeline that uses ANNOVAR5 and additional tools and databases. Patient BAB5311 underwent whole-exome sequencing in a clinical molecular laboratory. Raw sequence data (fastq files) were transferred to Baylor College of Medicine and analyzed using the bioinformatic approaches previously described for patient BAB3280.
Candidate variants in all the affected individuals were confirmed by Sanger sequencing after polymerase chain reaction amplification of the region containing the respective variants using specific primers. Segregation of the variants was tested in the 2 sets of parents and available unaffected siblings.
High-resolution single-nucleotide polymorphism arrays were performed for patients BAB5311, BAB5320, and BAB5321 using the Illumina HumanOmni2.5 BeadChip according to the manufacturer’s protocols. Analysis was performed using the genotyping application of the Illumina GenomeStudio software.
For sisters BAB3022 and BAB3280, analysis of the sequence data under a recessive model resulted in identification of 2 likely disease-causing variants in VRK1: chr14:96 391 443 (G>A) (c.G706A; p.V236M) and chr14:96 382 234 (G>A) (c.G266A; p.R89Q). No disease-associated variants were noted in any known HMSN-associated genes. Both variants were found in the affected sisters as compound heterozygous alleles that were inherited from nonconsanguineous heterozygous carrier parents (Figure 1E). The unaffected sister did not inherit either VRK1 mutation. Both variants are novel, not previously observed or reported in any database to our knowledge, and predicted to be damaging to protein function by independent bioinformatic algorithms.
Patient BAB5311 was found to have a homozygous nonsense mutation in VRK1: chr14:97 342 370 (C>T) (c.C1072T; p.R358X). This variant is a reported disease-causing mutation (rs137853063) as it had been previously identified in a patient with atypical pontocerebellar hypoplasia (PCH) type 1A (OMIM #607596).6 Currently, the frequency of this variant is thought to be approximately 0.007% in the general population and 0.0116% in the population of European descent.7 Both parents, who are self-reported to be of Ashkenazi Jewish ancestry but not related, and the unaffected twin sister of the patient were confirmed to be heterozygous carriers for this mutation. We performed high-resolution single-nucleotide polymorphism arrays on the affected patient and both parents to narrow a potential common haplotype region in which the mutation occurred. We confirmed by single-nucleotide polymorphism array data that the parents are not related, and we were able to narrow the common haplotype to a homozygous 500 673–base pair (g.chr14:97 072 050-97 572 723) genomic interval including only VRK1 and harboring the founder p.R358X mutation (Figure 3).
Whole-exome sequencing is a powerful research and diagnostic tool that allows identification of sequence variants in genetically and clinically heterogeneous conditions.8,9 Hereditary motor and sensory neuropathies represent a genetically heterogeneous group of disorders, with more than 40 disease-associated genes to date.1,10 All affected individuals in this study manifested an autosomal recessive axonal neuropathy associated with microcephaly and cerebral dysgenesis. The latter finding is inconsistent with mutations in known HMSN-associated genes. The patients did not meet clinical criteria for any known form of HMSN or other hereditary neurodevelopmental disorder, and extensive molecular investigations that included sequencing of known HMSN genes failed to delineate the molecular cause of the disease. Other possible disorders that may include neuropathy accompanied by microcephaly, such as some congenital disorders of glycosylation, were unlikely based on the lack of additional clinical features and normal findings on biochemical analyses.
Genomic studies identified many sequence variants in these 2 sisters (10 366 nonsynonymous variants in patient BAB3022 and 10 555 in patient BAB3280), of which they shared novel compound heterozygous mutations in 3 genes. However, only the compound heterozygous VRK1 mutations both could explain the observed phenotype and were consistent with mendelian expectations. It is a challenge to identify potential disease-causing variants in the myriad of variants produced by genomic sequencing approaches. Interpretation of candidate variants can be aided by prediction algorithms of the variant’s effect on the function and structure of the normal protein and the frequencies of rare and polymorphic variants in other individuals not affected by the observed phenotype. However, these data are most informative when applied in the context of a genetic model that considers the family history and overlapping clinical presentation between affected individuals.
The homozygous nonsense mutation in VRK1, identified in patient BAB5311, had been previously reported in only 1 consanguineous Israeli family with progressive microcephaly of prenatal onset, PCH, neurogenic muscle atrophy, and mild intellectual disability.6 Affected individuals in the family presented with early-onset ataxia and hyperreflexia followed by progressive muscle wasting and areflexia. The electrophysiological studies showed evidence of motor sensory neuropathy, but no further characterization was provided. Interestingly, this variant seems to have occurred as a founder variant in the Ashkenazi Jewish population, as both parents of our patient BAB5311 are of Ashkenazi Jewish ancestry and share the identical mutation and surrounding haplotype found in the family described by Renbaum et al.6
The PCHs are classified into categories PCH1 through PCH8 based on clinical findings observed on radiologic brain and spine imaging studies (ie, cerebellar and pontine hypoplasia), neurological signs and symptoms, and the natural history of the disease. Currently, classification of PCH does not consider underlying genetic causes such as mutations in genes like TSEN54, in which defects result in different PCH types.11 The original classification of PCH distinguished PCH1 from other forms of PCH based on early onset of respiratory distress, hypotonia, weakness, and central and peripheral motor dysfunction (the latter thought to occur due to anterior horn cell degeneration) and death occurring within the first few months of life.12- 14 Conversely, PCH2 through PCH8 present with progressive microcephaly and additional neurological findings including extrapyramidal dyskinesia, dystonia, and epilepsy, without motor or anterior horn cell involvement and degeneration.12,13
Concordant with the previous report of the family with mutations in VRK1, patients BAB3022, BAB3280, and BAB5311 share many clinical features with the described patients6 (microcephaly, peripheral neuropathy with secondary muscle atrophy), but their presentation is quite distinct owing to the lack of observed PCH on MRI, lack of central nervous system neurological symptoms (ataxia, hypertonia), and normal cognitive function. Additionally, although some features are shared by our cases and those that have been previously reported (microcephaly, hypotonia, motor neuron disease), the phenotype observed in our patients with mutations in VRK1 do not conform within the clinical definition of PCH1A as our patients do not have PCH and do have normal intellectual development.
The vaccinia-related kinase 1 gene encodes a serine/threonine kinase that is crucial for cell cycle progression and cell division and is proposed to be involved in nervous system development and maintenance. It is an early response gene that directly phosphorylates important transcription factors such as p53, c-Jun, and activating transcription factor 2 (ATF2).15 The gene is ubiquitously expressed in human tissues, with higher levels of expression in mitotically active cells. Null alleles of VRK1 result in embryonic lethality in Caenorhabditis elegans and Drosophila melanogaster, while hypomorphic alleles lead to sterility in mice due to the inability of germ cells to proliferate and differentiate.16- 18
The mutations observed in the 2 sisters of the first family, p.V236M and p.R89Q, localize to the kinase domain of the protein and are considered to be damaging by bioinformatic prediction algorithms. It is challenging to explain why VRK1 mutations result in pathology restricted to the nervous system. Nevertheless, mutations in many different genes involved in DNA replication, DNA damage response, and DNA single-strand break repair do result in neurological phenotypes. Moreover, mutations in genes whose proteins are important for DNA replication and cell cycle progression can result in syndromic microcephaly.
The VRK1 protein and p53 directly interact by phosphorylating each other in an autoregulatory loop.19 The VRK1 uniquely phosphorylates p53 in Thr18, stabilizing and activating it to interact with other transcription factors and preventing the binding of mouse double minute 2 homolog (MDM2).15 It has been proposed that phosphorylation of p53 in response to oxidative stress is connected to the kinase function of cleavage and polyadenylation factor I subunit 1 (CLP1).20 Disruption of the latter was shown recently to result in fatal deterioration of motor function due to progressive loss of motor neurons potentially resulting in abnormal apoptosis.19 We hypothesize that a similar mechanism of abnormal neural apoptosis may explain the neurological phenotype observed in patients with VRK1 mutations. Some evidence of this process was found on brain MRI studies for patient BAB5311 in which enlarged Virchow-Robin/perivascular spaces were observed in addition to progressive volume loss (Figure 2A). More experimental studies are required to better characterize the role of VRK1 in the nervous system and elucidate the biological consequences of mutations in human subjects. In conclusion, VRK1 may be a novel HMSN locus that can be associated with a complex peripheral neuropathy phenotype: autosomal recessive axonal motor sensory neuropathy and microcephaly.
Corresponding Author: James R. Lupski, MD, PhD, DSc, Department of Molecular and Human Genetics, Baylor College of Medicine, One Baylor Plaza, 604B, Houston, TX 77030 (email@example.com).
Accepted for Publication: July 31, 2013.
Published Online: October 14, 2013. doi:10.1001/jamaneurol.2013.4598.
Author Contributions: Drs Gonzaga-Jauregui and Lupski had full access to all of the data in the study and take responsibility for the integrity of the data and the accuracy of the data analysis.
Study concept and design: Lupski, Gibbs.
Acquisition of data: Gonzaga-Jauregui, Lotze, Jamal, Penney, Woodbury, Raymond, Adesina, Jhangiani, Muzny, Boerwinkle, Lupski, Wiszniewski.
Analysis and interpretation of data: Gonzaga-Jauregui, Campbell, Pehlivan, Hunter, Reid, Wiszniewski.
Drafting of the manuscript: Gonzaga-Jauregui, Lupski, Wiszniewski.
Critical revision of the manuscript for important intellectual content: All authors.
Statistical analysis: Gonzaga-Jauregui, Campbell.
Obtained funding: Lupski, Gibbs.
Administrative, technical, or material support: Jhangiani.
Study supervision: Muzny, Lupski.
Conflict of Interest Disclosures: Dr Hunter receives royalties for a chapter on pediatric neuroimaging in UpToDate. Dr Boerwinkle received honoraria from The Foundation of Rome for a project review. Dr Lupski is a consultant for Athena Diagnostics, 23andMe, and Ion Torrent Systems, Inc and holds multiple US and European patents for DNA diagnostics. Dr Gibbs was an owner of SeqWright and is an advisor to GE Healthcare/Clarient and the Allen Institute for Brain Science. The Medical Genetics Laboratories (http://www.bcm.edu/geneticlabs/) and Whole Genome Laboratory (http://www.bcm.edu/geneticlabs/index.cfm?PMID=21319) of the Department of Molecular and Human Genetics at Baylor College of Medicine derive revenue from molecular diagnostic testing. No other disclosures were reported.
Funding/Support: This work was supported by grant U54HG006542 from the National Human Genome Research Institute to the Baylor-Hopkins Center for Mendelian Genomics and grant U54HG003273 from the National Human Genome Research Institute to the Human Genome Sequencing Center, Baylor College of Medicine. Dr Hunter is supported by grants R01 NS21889, R03 NS069943, P01 NS056202, and PTO 74693P2 from the National Institutes of Health and by the Simons Foundation. Dr Boerwinkle is supported by grants 5 RC2 HL 102419-02, 5 N01HC 55016-40-0-1, 5 U01 HG 004803-04, 5 R01 HL 072810-09, 1 R01 HL 106034-02, 5 U01 HL 096917-02, 2 U01 GM 074492-08, 2 U54 HG003273-09, and 1 U54 HG006542-01 from the National Institutes of Health. Dr Lupski is supported by grants R01NS058529 from the National Institute of Neurological Disorders and Stroke and U54HG006542 from the National Human Genome Research Institute. Dr Gibbs is supported by the National Human Genome Research Institute and the National Cancer Institute. Dr Wiszniewski is supported by Career Development Award K23NS078056 from the National Institute of Neurological Disorders and Stroke and grant 2012/06/NZ2 from the Ministry of Science and Higher Education of Poland.
Role of the Sponsor: The funding organizations had no role in the design and conduct of the study; collection, management, analysis, and interpretation of the data; or preparation, review, or approval of the manuscript.
Additional Contributions: We thank the patients and families who participated in this study.