Novel PLA2G6 mutations associated with an exonic deletion due to non-allelic homologous recombination in a patient with infantile neuroaxonal dystrophy.

Novel PLA2G6 mutations associated with p.Asp283Asn and a unique intragenic deletion of exons 4 and 5 due to non-allelic homologous recombination were identified in a Japanese female patient with typical infantile neuroaxonal dystrophy. The patient showed progressive tetraplegia beginning at 9 months. An electroencephalogram showed a diffuse increase in fast waves, and brain magnetic resonance imaging showed progressive brain atrophy and T2 hypointensity in the globus pallidus.

Infantile neuroaxonal dystrophy (INAD; MIM #256600) is a rare autosomal-recessive neurodegenerative disorder involving both the central and peripheral nervous system. INAD is diagnosed when patients show several clinical features, including the onset of symptoms before 3 years of age, psychomotor regression, hypotonia, symmetrical pyramidal tract signs, a relentlessly progressive course leading to spastic tetraplegia, visual impairment with nystagmus due to optic atrophy, and dementia.[1,2] Brain pathological examinations demonstrated characteristic findings of axonal spheroids.[3] This disorder typically occurs because of mutations in the phospholipase A2 group VI gene (PLA2G6) located on the 22q13.1 region, which encodes a calcium-independent phospholipase A2 enzyme involved in phospholipid remodeling and catalyzes the hydrolysis of glycerophospholipids.[4] Recent molecular investigations identified nearly 100 PLA2G6 mutations not only in patients with INAD but also in patients with parkinsonism.[1,5-25] Thus, the disease concept has been expanded to include all of these disorders into a group called neurodegeneration with brain iron accumulation.[2,26] Currently, these broad clinical features due to PLA2G6 mutations are recognized as a phenotypic spectrum of PLA2G6-associated neurodegeneration (PLAN). Recently, we encountered a Japanese female patient with a typical form of INAD and identified compound heterozygous mutations.

At present, the patient is of 2 years and 8 months old. She is the only child of healthy non-consanguineous parents. Her family history was unremarkable. She was born after a 40-week-long uneventful gestation with a weight of 2,840 g without asphyxia. Her initial development was normal during early infancy: she attained head control at 3 months and sat without support at 7 months. Her development, however, slowed after she began to crawl at 9 months. Although she stood with a support at 14 months, her developmental deterioration gradually emerged at ~17 months. At 19 months, hypotonia was noted, particularly in the lower extremities. She lost the ability to sit and roll, and showed decreased voluntary movements at 23 months.

Neurological examination revealed generalized hypotonia, hyperreflexia in all extremities and bilateral Babinski reflexes. Electroencephalogram (EEG) showed a diffuse increase in fast waves without apparent rhythmic activity (Figures 1a and b). Brain magnetic resonance imaging (MRI) showed atrophy of the cerebellum (Figures 1d–g). Auditory brainstem responses were almost lacking and the cortical component of short-latency sensory-evoked potentials were not detected. Electromyogram showed signs of denervation. Other laboratory examinations revealed unremarkable results.

Figure 1

Clinical examinations of the present patient. Electroencephalogram (EEG) of the patient (a–c). Diffuse increase of low-amplitude 15–30 Hz fast waves are demonstrated during wakefulness at 23 months of age (a). At 2 years and 8 months, extremely fast waves during wakefulness (b) and multifocal spike-waves during sleep (c, arrowheads) are recorded. Brain magnetic resonance imaging (MRI) of the patient (d–i). Axial T1 (d) and T2 (e) images and a sagittal T1 (f) and a coronal T2 (g) images examined at 23 months. Mild enlargement of the lateral ventricles (d, e) and cerebellar atrophy (f, g) are demonstrated. These findings are remarkable in axial T1 (h) and T2 (i) images examined at 2 years and 8 months. Progressive T2 hypointensity in the globus pallidus is noted (i). ECG, electrocardiogram.

She gradually lost her visual ability to follow objects with emergence of strabismus and transient nystagmus at 2 years and 8 months. Her optic nerves showed initial signs of atrophy. Tetraplegia was progressive and the bilateral upper extremities were rigid. The EEG revealed a further increase in fast waves and multifocal spike-waves, particularly in the bilateral temporal regions (Figure 1c). She had no apparent seizures, although she often showed episodic stiffening and extension of the bilateral upper extremities associated with no ictal epileptic patterns in EEG. Brain MRI revealed a slight progression of diffuse cerebral atrophy and T2 hypointensity in the globus pallidus (Figures 1h and i). Based on these findings, she was clinically diagnosed with INAD.

Although PLA2G6 was listed as a candidate target of molecular diagnosis, the contribution of other genes could not be excluded. Thus, targeted re-sequencing was performed as the first step to screen single-nucleotide variants using TruSight One v1.0 sequencing panel (Illumina, San Diego, CA, USA), which includes 125,396 probes aimed to capture 11,946,514-bp targeted exon regions consisting of 4,813 genes associated with known clinical phenotypes. This study was performed in accordance with the Declaration of Helsinki Principles, and the ethics committee of Tokyo Women’s Medical University approved this study. Blood samples were obtained from the patient and her parents after receiving written informed consent thorough careful genetic counseling regarding the appropriate dealing of genetic information and possible incidental findings. Genomic DNA was extracted from blood samples and used for further examination. After constructing the sequence library using 50 ng of genomic DNA, the MiSeq next-generation sequencer (Illumina) was used to sequence 151-bp paired-end reads according to the manufacturer’s instructions. The extracted data were mapped to a reference genome (GRCh37/hg19) using BWA Enrichment v1.0 cloud software (Illumina).

The total aligned bases were 2.08-Gb with a mean coverage depth of 83.7 (the target coverage at 20× was 95.0%). A total of 7,919 variants were extracted. The extracted variants were annotated and filtered using Variant Studio software (Illumina) as described previously.[27] As a result, 20 variants were considered as candidates (Supplementary Table S1). Among them, a missense mutation, NM_003560.2(PLA2G6):c.847G>A (NP_003551.2(PLA2G6):p.Asp283Asn), was included (Figure 2a). Although the same variant has never been reported, this codon is located on the ankyrin repeat region (which is crucial for protein–protein interaction)[4,28] and was affected by a previously reported disease-causing mutation, p.Asp283Gly.[4] The prediction scores suggested pathogenesis, i.e., SIFT (0; deleterious) and PolyPhen-2 (0.999; probably damaging). From these findings, we considered this variant to be a disease-causing mutation. Subsequent Sanger sequencing confirmed that this variant was inherited from the patient’s father (Figure 2b). However, a maternally inherited variant was not detected through targeted re-sequencing.

Figure 2

Molecular analysis results. (a) Image view using Integrative Genomics Viewer (IGV) demonstrates an altered ‘T’ nucleotide in approximately half of the reads. The affected amino-acid ‘D’ is conserved among species. (b) Electropherograms of Sanger sequencing. The present patient and her father show heterozygous for c.847G>A. Thus, her father is an obligate carrier of this mutation. (c) Agarose gel electrophoresis of RT-PCR amplicons. Compared with the normal control (N), the patient (P) shows an aberrant short band (arrow). M4; phi-X HaeIII digest. (d) An electropherogram of Sanger sequencing for an aberrant short band demonstrates a skipping of exons 4 and 5. (e) Agarose gel electrophoresis of long PCR products. In addition to the normal bands with an expected size of 4,949 bp, the patient (P) and her mother (Mo) show aberrant short bands (arrow). Thus, the mother is an obligate carrier of this large deletion. (f) An electropherogram of Sanger sequence for the short band extracted from the agarose gel. A 12-bp homologous region is identified in the breakpoint. The genomic positions are indicated on the sequences. Because PLA2G6 is encoded in an anti-sense direction, the locations of genomic numbers are inverted. (g) The result of homology searching between two breakpoints. “Query” and “Sbjct” indicate the nucleotide sequences around the breakpoints in intron 3 and 5, respectively. An identity of 86% is calculated between them. The box indicates the homologous region in the breakpoint. Fa, patient’s father; M1, lambda/Hind III digest; N, normal control sample.

To identify an undiagnosed mutation in the homologous allele, PLA2G6 mRNA was analyzed by reverse transcription (RT)-PCR amplification. Total RNA was extracted from Epstein-Barr virus-transfected immortalized lymphocytes established from the patient and transcribed into complementary DNA. The RT-PCR product using the primer set (5ʹ-TGCTACCCTTCTATGAGAGC-3ʹ and 5ʹ-TGGTGCTGTTCACGTTGCAG-3ʹ) showed an aberrant shorter band in addition to a normal band (Figure 2c). DNA extracted from the short band was analyzed by Sanger sequencing and the result revealed deletions in exons 4 and 5 (Figure 2d), indicating a 372-bp in-flame deletion, c.426_797del372 (p.Ser142_Gly266delinsArg). Because no splicing mutations located at the exon–intron boundaries nearby this region were detected, we suspected exonic deletions of this region. Next, the genomic region encompassing exons 4 and 5 was amplified by PCR amplification using the primer set (5ʹ-ACCAGTTGGCCATCTTGTGC-3ʹ and 5ʹ-AAGGAGCACTGAAGCCATCG-3ʹ), which was designed to amplify a 4,949-bp PCR product, and an abnormal short band was identified (Figure 2e). DNA extracted from the short band was analyzed by Sanger sequencing and a genomic deletion of exons 4 and 5 was confirmed (Figure 2f). The size of aberrant short band was 640 bp, indicating a 4,309-bp genomic deletion. This deletion was also identified in the sample from the mother (Figure 2e). Thus, compound heterozygous mutations of PLA2G6 inherited from both parents were confirmed and a molecular diagnosis of PLAN was established in this patient.

A similar exonic deletion was previously reported by Tonelli et al. in exons 5 and 6.[11] The breakpoints of the reported deletion were analyzed and the homologous sequence was identified in the breakpoints. Based on their findings, they considered Alu-repeats as the genomic characteristic of this region. In this study, the exonic deletion involved exons 4 and 5 was identified, and the breakpoint sequence showed a 12-bp homologous sequence (Figure 2f). Using Genetyx software (Genetyx Corporation, Tokyo, Japan), we further analyzed the homology of the neighboring the 237-bp genomic regions encompassing the 12-bp homologous sequence, and 86% identity was confirmed (Figure 2g). This indicates that the mechanism of the identified exonic deletion is quite similar to that reported by Tonelli et al.[11] The genomic deletion between the similar genomic sequences suggests predisposition of the genomic rearrangements in this region.

Previously, nearly 100 PLA2G6 mutations were reported in patients with PLAN and/or parkinsonism (Supplementary Table S2). The identified mutations are distributed throughout the PLA2G6 complementary DNA sequence and there is no clear genotype–phenotype correlation. Particularly, Arg632Trp was identified in patients with both of INAD and parkinsonism.[4,10] Thus, phenotypic consequence of PLA2G6 mutations may be related to the combinations of the heterozygous mutations.