Mutations in Diphosphoinositol-Pentakisphosphate Kinase PPIP5K2 are associated with hearing loss in human and mouse

Autosomal recessive nonsyndromic hearing loss is a genetically heterogeneous disorder. Here, we report a severe-to-profound sensorineural hearing loss locus, DFNB100 on chromosome 5q13.2-q23.2. Exome enrichment followed by massive parallel sequencing revealed a c.2510G>A transition variant in PPIP5K2 that segregated with DFNB100-associated hearing loss in two large apparently unrelated Pakistani families. PPIP5Ks enzymes interconvert 5-IP7 and IP8, two key members of the inositol pyrophosphate (PP-IP) cell-signaling family. Their actions at the interface of cell signaling and bioenergetic homeostasis can impact many biological processes. The c.2510G>A transition variant is predicted to substitute a highly invariant arginine residue with histidine (p.Arg837His) in the phosphatase domain of PPIP5K2. Biochemical studies revealed that the p.Arg837His variant reduces the phosphatase activity of PPIP5K2 and elevates its kinase activity. We found that in mouse inner ear, PPIP5K2 is expressed in the cochlear and vestibular sensory hair cells, supporting cells and spiral ganglion neurons. Mice homozygous for a targeted deletion of the Ppip5k2 phosphatase domain exhibit degeneration of cochlear outer hair cells and elevated hearing thresholds. Our demonstration that PPIP5K2 has a role in hearing in humans indicates that PP-IP signaling is important to hair cell maintenance and function within inner ear.


Introduction
Hearing loss (HL) is a heterogeneous neurosensory deficiency that occurs at all ages with varying severities, affecting 1 in 500 newborns and >360 million people worldwide [1]. Many forms of HL are inherited, and~400 syndromic forms occur with linked medical comorbidities. Genetically complex non-syndromic recessively inherited HL (NSRHL) comprises~75% of hereditary deafness [2]. RNA profiling studies revealed expression of 18,133 genes in the organ of Corti (OC) hair and supporting cells [3], underscoring the complexity of inner ear development and function. Elucidating the roles of each of these genes by conventional methods would require generating mouse or other animal models for each. By contrast, genetic studies of human families segregating HL has significantly helped in overcoming the barriers of sheer scale, functional redundancy, and non-essential roles of some genes. Discovering a variant associated with human deafness rules out functional redundancy at least in the auditory system. Genetic and functional studies of the deafness associated proteins have been pivotal in elucidating various molecular networks essential for hearing [4][5][6].
Our in vitro biochemical studies indicate that the p.Arg837His deafness-associated variant of PPIP5K2 has reduced phosphatase activity and increased kinase activity. As a consequence, the variant has the capacity to synthesize more of the IP8 signal in vivo. An elevated cellular production of IP8 can be recapitulated in mice harboring a homozygous deletion of the phosphatase domain of PPIP5K2, thereby unmasking higher kinase activity. These animals exhibited degeneration of cochlear outer hair cells and progressive hearing loss. To the best of our knowledge, this is the first study demonstrating the necessity of PP-IP metabolism for inner ear development and hearing function.

Families PKDF041 and PKDF751 have prelingual profound HL
Two large apparently unrelated families PKDF041 and PKDF751 segregating NSRHL were enrolled from Punjab province of Pakistan (Fig 1A and 1B). Affected individuals of both families have prelingual bilateral severe to profound sensorineural hearing loss (Fig 1C and 1D). According to family history no hearing was noted since birth in all of the affected individuals. Tympanometry of three of the affected [V:2 (19 yrs), V:3 (23 yrs), V:5 (26 yrs)], and one normal hearing individual (V:7; 22 yrs) of family PKDF751 revealed no abnormalities of the tympanic membrane or middle ear. However, these affected individuals (V:2, V:3, V:5) failed a transient evoked otoacoustic emission test, which is indicative of defective outer hair cell function. Romberg and tandem gait tests did not reveal any overt vestibular dysfunction in affected individuals of both families. Funduscopy examinations of affected individuals of PKDF041, V:5 (42 years), VI:3 (16 years), VI:4 (10 years) and VI:5 (20 years) revealed no evidence of retinitis pigmentosa. In addition, ERG a-and b-waves amplitudes for the affected individuals [V:5 (42 yrs), V:7 (30 yrs)] of PKDF041, were normal [13], further excluding a retinal degeneration phenotype.

Linkage analysis mapped a new deafness locus, DFNB100 on chromosome 5q
Linkage analysis of NSRHL segregating in family PKDF041 with STR markers across the genome, previously revealed linkage to DFNB49 locus on chromosome 5, with a maximum two-point lod score (Zmax) of 4.44 (θ = 0) for marker D5S2055 [13]. Subsequently, mutations in TRIC were identified as the cause of DFNB49-linked HL [14,15]. Full sequencing of TRIC in the genomic DNA from two affected individuals from the PKDF041 family along with a normal hearing sibling did not reveal any pathogenic variants, suggesting there is an additional gene in the PKDF041 linkage interval in which a pathogenic variant is associated with deafness.
Subsequently, genome wide linkage analysis in DNA samples from a second family, PKDF751, revealed a large region of homozygosity shared among the affected individuals on human chromosome 5q12.2-q23.3, which completely overlaps with the linkage interval defined by family PKDF041 (S1 Fig). Similar to family PKDF041, mutation in the protein coding exons, non-coding exons or in the splice junctions of TRIC were not detected in the affected individuals of family PKDF751. Therefore, the HUGO nomenclature committee assigned the DFNB100 designation for the locus defined by families PKDF041 and PKDF751 (S1 Fig). Besides TRIC, other candidate genes include OCLN, SHROOM1, GPR98, KCNN2 and SLC12A2 (S1 Fig) [16][17][18][19][20][21]. However, sequencing of all these genes did not reveal a potentially pathogenic variant among the affected individuals of families PKDF041 and PKDF751. Rather than continue hierarchical sequencing of candidate genes based on function or expression in the inner ear, we employed an exome sequencing approach to identify a mutant variant responsible for deafness at the DFNB100 locus.

A missense variant in PPIP5K2 is associated with DFNB100-deafness
To enrich and capture coding regions, we used genomic DNA samples from one affected individual and one unaffected control from family PKDF751, and performed massively parallel sequencing of the exome. Copy number variants (CNVs) analysis did not reveal any indels (!50bp) associated with deafness within DFNB100 linkage interval. Our filtering of the exome sequencing data revealed three genes with predicted pathogenic changes within DFNB100 linkage interval (S1 Table). Sanger sequencing revealed a single point mutation (c.2510G>A) in the PPIP5K2 gene (NCBI RefSeq NM_001276277) segregating with HL in family PKDF751 (Fig 2A), which is predicted to replace an evolutionarily-conserved arginine with histidine at residue 837 (NP_001263206; Fig 2B, S2 Table). Sanger sequencing of all coding and non-coding exons of PPIP5K2 in the DNA samples from family PKDF041 confirmed association of the same c.2510G>A variant with deafness. SNPs linked to PPIP5K2 were genotyped in affected individuals of the PKDF041 and PKDF751 families, and the flanking haplotype was consistent  Table). We did not detect this variant in 180 ethnically-matched Pakistani control samples, nor in small publicly-available databases (1000 Genome [22], or NHLBI_EP [23]). At present, the c.2510G>A variant of PPIP5K2 is at a low frequency (0.000146) in the ExAC database [24], and thus c.2510G>A is not a common polymorphism.

The PPIP5K2 R837H variant modifies the catalytic activities of the enzyme
Human PPIP5K2 encodes at least eleven alternatively spliced isoforms, all of which are affected by the p.(Arg837His) substitution identified in the two DFNB100 HL families ( Fig 2B). PPIP5K2 has three distinct functional regions: a kinase domain, a phosphatase domain and an intrinsically disordered region that may mediate protein-protein interactions ( Fig 2B). The Arg837His substitution is located in the phosphatase domain ( Fig 2C). Human PPIP5Ks is a mutually competitive kinase and phosphatase that interconverts PP-IPs ( Fig 2C). In silico analyses predict that the p.Arg837His residue is likely of structural and/or functional importance (S2 Table). Indeed, Arg837 is conserved throughout the animal kingdom (Fig 2D), and lies nine residues C-terminal to a conserved histidine residue that is catalytically-essential for IP8 phosphatase activity in a yeast orthologue of PPIP5K2 [25]. Thus, p.Arg837 is a candidate for binding the negatively charged PP-IP substrate. The substitution of histidine would be less effective in this role as its side chain is less polar at physiological pH [26]. Another possibility is that p.Arg837 acts in a structurally-stabilizing salt bridge, which could be disrupted by substitution with His [27].
To pursue these ideas, we first assayed directly the impact of the p.Arg837His variant upon the phosphatase activity of PPIP5K2 in vitro, using recombinant, FLAG-tagged WT and PPIP5K2 R837H which were expressed in HEK293 cells and then immuno-affinity purified ( Fig  2E). The IP8 phosphatase activity of the PPIP5K2 R837H protein was reduced by 21% ( Ã p< 0.05) as compared to the WT protein ( Fig 2E). PPIP5K2 interconverts 5-IP7 to IP8 through the actions of mutually competing phosphatase and kinase domains. Thus, we hypothesized based upon the in vitro kinase assays of 5-IP7 phosphorylation that PPIP5K2 R837H protein would accumulate more IP8 than WT PPIP5K2. That proposal presumes that a potentially higher level of accumulation of IP8 by the PPIP5K2 R837H protein is not negated by an increased rate of IP8 dephosphorylation. However, that self-correction did not occur. Assays that used PPIP5K2 R837H accumulated 60% higher levels of IP8 compared to assays that contained WT protein ( Fig 2F). The proportionately higher effect of the p.Arg837His variant upon the kinase activity may reflect an impact upon conformational coupling between the two domains [8], in addition to the reduction in phosphatase activity. Note that, in vitro, the phosphatase and kinase domains of PPIP5K2 can also be shown to interconvert IP6 and 1-IP7 [12]. However, these particular reactions proceed relatively slowly [12]; their general biological significance is in doubt, especially as levels of 1-IP7 in mammalian cells are virtually undetectable [11].
families. (B) The human PPIP5K2 gene has 33 exons. Alternate splicing of PPIP5K2 gives rise to many different transcripts that are predicted to affect the composition of the non-catalytic intrinsically disordered region (IDR). The non-coding 5' and 3' UTRs are gray rectangles while sequences encoding the kinase, phosphatase (pase) and IDR domains, and other coding regions of exons are denoted by orange, light green, dark green, and black boxes, respectively. The nucleotide and amino acid variants of PPIP5K2 are also shown. The blue bar under the phosphatase domain indicates the antigenic region of the commercially available polyclonal antibodies. (C) Schematic representation of PPIP5Ks biochemical function. Arrows indicate direction of the reaction. (D) PPIP5K2 amino acid conservation from eleven vertebrates and one non-vertebrate (Drosophila). Conserved amino acids are shaded in gray, while the p.Arg837 amino acid affected due to c.2510G>A variant, identified in our DFNB100 families, is noted with an arrow head. (E) The phosphatase activities of FLAG-tagged WT PPIP5K2 (21.5 ± 1.5 nmol/mg protein/min) and PPIP5K2 R837H (16.9 ± 2.1 nmol/mg protein/min) were determined as described in the Methods section ( Ã p<0.05, n = 5; mean ± SE). (F) The kinase activities of FLAG-tagged WT PPIP5K2 (0.9 ± 0.21 nmol/mg protein/min) and PPIP5K2 R837H (1.44 ± 0.29 nmol/mg protein/min) were determined as described in the Methods Section ( Ã p<0.03, n = 5; mean ± SE). https://doi.org/10.1371/journal.pgen.1007297.g002

Ppip5k2 is ubiquitously expressed in mice
Our analysis by real-time PCR showed Ppip5k2 and Ppip5k1 to be widely expressed in mouse tissues, including inner ear ( Fig 3A, S3A Fig). Transcriptome analysis indicates that Ppip5k2 is the major Ppip5k gene to be expressed in the mouse cochlear and vestibular sensory hair cells, supporting cells as well as in ganglion neurons, while levels of Ppip5k1 expression in these cells are up to 100-fold lower (gEAR, SHIELD) [28]. Similarly, PPIP5K2 immunoreactivity, assessed using an antibody directed against a region within the phosphatase domain (Fig 2B), was observed in all these cell types at various developmental stages (Fig 3). Moreover, all three layers of the stria vascularis (marginal, intermediate and basal cells) exhibited PPIP5K2 immunoreactivity (Fig 3B and 3C). The distribution of PPIP5K1 immunoreactivity within the inner ear, assessed using isoform-specific antibodies (S2 Fig

Disruption of murine Ppip5k2 promotes auditory deficit and hair cell degeneration
To determine the role of PPIP5K2 in the inner ear, we obtained a mouse from the Knockout Mouse Phenotyping Program (KOMP) that has a gene trap cassette with a LacZ reporter in the intron between exons 13 and 14 of Ppip5k2 (Fig 4A). The gene trap leads to the translation of a protein in which the N-terminal kinase domain is intact, but the adjoining phosphatase domain is truncated to just 84 residues (Fig 4A), which is insufficient to encode phosphatase activity [29,30]. To study the impact upon the enzymatic properties of the protein product of the targeted allele, we prepared a recombinant, C-terminally truncated human PPIP5K2 (PPIP5K2  ). This is equivalent to the murine protein that is expected to be expressed by the Ppip5k2 K^a llele. The shorter human protein converted 5-IP7 to IP8 at a 23-fold higher rate than that catalyzed by the kinase domain of the full-length human PPIP5K2 (Fig 4B). We therefore named the murine allele as Ppip5k2 K^. It is impractical to assay the impact of the Ppip5k2 K^a llele upon PP-IPs in vivo, due to PP-IPs sub-to low-micromolar levels, and the relative insensitivity of mass assays for these molecules. Instead, we used [ 3 H]-inositol radiolabeled PPIP5K -/-HEK293 cells [12], as a host for over-expression of either WT human PPIP5K2 or the truncated PPIP5K2 1-466 mutant. The latter supported the synthesis of 4-fold (n = 3) higher cellular IP8 levels than those sustained by the WT protein (Fig 4C), despite the expression level of truncated PPIP5K2 1-466 being much lower than that of WT PPIP5K2 (S4 Fig). Based on these findings, we conclude that the Ppip5k2 mutant mice likely expresses a "hyper kinase" activity. Note that the changes in cellular levels of IP7 and IP8 resulting from over-expression of PPIP5K2 R837H in PPIP5K -/-HEK293 cells [12] were similar to those observed after over-expression of WT PPIP5K2 ( Fig  4C). The catalytic differences between these particular enzymes are less prominent (Fig 2E and  2F), making it harder to detect their differential impact upon PP-IP turnover in intact cells. PP-IP levels will also be influenced by cell-type dependent differences in expression of other PP-IP kinases and phosphatases.
In Ppip5k2 K^/K^m ice, we observed a significant upregulation of Ppip5k1 message in several tissues at P150, including a nearly 3-fold elevation in the inner ear (S3F Fig). Any increased expression of Ppip5k1 can be predicted to further enhance the synthesis of IP8 [12]. However, the metabolic consequences of this effect are likely minor in the inner ear, where relative levels of Ppip5k1 are 100-fold less than those of Ppip5k2 (see above).
To characterize inner ear function, we measured auditory-evoked brainstem responses (ABR) in Ppip5k2 K^/ K^a nd WT mice at various developmental stages (S5 Fig). No statistically significant difference in hearing thresholds was observed in Ppip5k2 K^/ K^a s compare to WT at postnatal day 60 (P60), when tested for broad band clicks as well as pure tone burst of frequencies 8, 16, and 24 kHz (Fig 4D). However, Ppip5k2 K^/K^m ice had slightly elevated, but not statistically significant, thresholds at 32 kHz compared to WT mice (Fig 4D, top left). The difference progressively increased with age. At P90, the threshold at 32 kHz of Ppip5k2 K^/Km ice were statistically significantly (p<0.05) elevated as compare to WT (Fig 4D, top right). At  P120 and P150, significantly (p<0.001) higher threshold was observed at 24 kHz as well in Ppip5k2 K^/K^ ( Fig 4D, bottom), which suggest progression of hearing loss from base towards middle cochlear turn. The Ppip5k2 +/K^m ice at P150 also exhibited significantly higher threshold (p < 0.05) at 24 kHz, which suggest that reduced amount of WT protein is also not sufficient to maintain hearing function in older mice. As anticipated, considering the C57BL/6 genetic background, which is susceptible to age-related HL [31,32], at P120 and P150, the thresholds for 32 kHz were observed to be elevated in both genotypes. We also compared the ABR wave I amplitudes and latencies in control and Ppip5k2 K^m utant mice (Fig 5). Significant reduction in wave I amplitudes at higher frequency (24kHz) was observed in older homozygous mutant mice. Overall, ABR auditory data from mice are consistent with the conclusion that elevated PPIP5K2 kinase activity is a causal factor for inherited auditory defects, although the deficit in the Ppip5k2 K^/K^m ice is less severe than that observed in humans homozygous for the p.(Arg837His) substitution.
To determine whether the auditory deficit in Ppip5k2 K^/K^m ice was secondary to hair cell loss, spiral ganglion degeneration or stria vascularis atrophy, we analyzed whole mounts preparations as well as serial sections of the cochlear tissue at P150 (Fig 6). The cytoarchitecture and morphology of the cochlea in Ppip5k2 K^/K^m ice appear normal at P150. No obvious degeneration of inner (IHCs) and outer (OHCs) hair cells was observed in the apical and middle cochlear turn (Fig 6). However, consistent with high frequency hearing loss observed at P150, a greater degree of degeneration of OHCs in the basal coil was evident in Ppip5k2 +/Kâ nd Ppip5k2 K^/K^m ice (Fig 6). Although no difference was observed in the number of IHCs, a statistically significant (p<0.01) reduction in OHCs was observed in the basal turn (Fig 5B). No obvious indications of stria vascularis atrophy or spiral ganglion neuron degeneration were found in Ppip5k2 K^/ K^m ice at P150 (S6 Fig). Thus, the degeneration of OHCs is likely contributing to the progressive, high frequency hearing deficit in these Ppip5k2 K^/K^m ice. Overall, our data indicate that genetic perturbation of PP-IP turnover has consequences for cellular signaling that promotes an auditory deficit, in both humans and mice.

Discussion
This study points to the importance of PP-IP signaling and the role of the bifunctional PP-IP kinase/phosphatase, PPIP5K2, for normal mammalian inner ear function. This conclusion is based upon a novel single variant [c.2510G>A, p. (Arg837His)] in the PPIP5K2 gene that we associated with autosomal recessive, nonsyndromic, prelingual sensorineural deafness. Biochemical analyses revealed that, compared to WT human PPIP5K2, the PPIP5K2 R837H variant exhibited lower phosphatase activity and higher kinase activity, indicating it promotes increased metabolic flux from 5-IP7 to IP8 in vivo. The availability of a gene-trap murine model for elevated PPIP5K2 kinase activity is consistent with the association of p.(Arg837His) with an auditory deficit.
Humans who are homozygous for the c.2510G>A variant failed transient evoked otoacoustic emission tests, indicating an outer hair cell functional defect. Thus, it is significant that, in the Ppip5k2 K^/K^m ice, we found that hearing loss is associated with selective degeneration of cochlear outer hair cells. How might the latter phenomenon be promoted by elevated IP8 levels? One possibility is that hair cell loss is due to apoptosis, since previous work has shown that increased PP-IP synthesis is pro-apoptotic [33]. Another possibility emerges from recent work showing that loss of outer hair cells can result from bioenergetic imbalance [34]; the latter could also result from perturbation to PP-IP turnover and signaling [11,12]. Deficits in vesicle trafficking and endocytosis might affect the function of the inner hair cell ribbon synapse [35,36], and thus contribute to hearing loss observed in the Ppip5k2 K^m ice.
It is intriguing that the hearing deficit in the Ppip5k2 K^/K^m ice is less severe than that observed in humans homozygous for the p.(Arg837His) variant. Perhaps in mice there is greater functional redundancy of PPIP5K2 with other proteins. The particular susceptibility of humans to a single arginine-to-histidine substitution in PPIP5K2 may also reflect genetic or environmental factors. Elucidation of the causes of this dissimilarity may reveal molecular or cellular pathways for potential interventions to prevent HL in DFNB100 families. Future studies with mice in which the Ppip5k2 gene is targeted on a genetic background resistant to agerelated HL, may also help in deciphering the roles of PP-IPs in development, maintenance and function of sensory cells in the inner ear.
PPIP5Ks are pivotal enzymes for regulating PP-IP turnover. Their actions at the interface of cell signaling and bioenergetic homeostasis can impact many biological processes [11]. Yet surprisingly, as far as we are aware, our work provides the first description of any amino acid variant in either PPIP5K1 or PPIP5K2 that is both functionally-significant and associates with a human disorder. Two candidate single nucleotide polymorphisms, rs35671301 (p.Ser419Ala) and rs17155147 (p.Thr1267Met), in human PPIP5K2 were reported to be enriched in individuals with autism spectrum disorder [37, 38]. However, in those studies, there is no empirical or in silico evidence that either amino-acid replacement might alter affect protein function. Thus, PPIP5K2 is given new clinical significance by our observations.
In conclusion, the identification of the PPIP5K2 variant as a risk factor for deafness raises the possibility that there may be other variants in its coding sequence that exhibit a similar hearing deficit. Moreover, our study reinforces the concept that studies into PP-IP turnover and signaling are relevant to human health and well-being.

Family enrollment and clinical evaluation
This study was conducted under IRB-approved protocols, by IRB committees at University of Maryland School of Medicine (UMSOM), USA (HP-00059851), National Centre of Excellence in Molecular Biology, Lahore, Pakistan (FWA00001758), and the National Institutes Health, USA (Combined Neuroscience Blue Panel IRB; OH-93-N-016), and in accord with the Declaration of Helsinki for the release of clinical information, family history, and blood draw. Written informed consent was obtained from all adult participants and parents of minor subjects. Genomic DNA was extracted from blood samples. Clinical histories were obtained for all the individuals of families PKDF041 and PKDF751. A general physical examination was performed to evaluate general health. Pure tone audiometry tests for air and bone conduction were performed at frequencies from 250 to 8,000 Hz. Tympanometry and otoacoustic emission tests were performed to assess middle ear and cochlear function. Vestibular function was evaluated by testing tandem gait ability and by using the Romberg test.

Exome sequencing and segregation analysis
Exome sequencing on the genomic DNA of one affected individual of family PKDF751 was performed and analyzed as described previously [39].

Fig 5. Ppip5k2 K^/K^m ice exhibit late onset high-frequency decline in ABR wave I amplitudes.
Average pure tone-evoked ABR wave I amplitudes and latencies in Ppip5k2 +/+ , Ppip5k2 +/K^, and Ppip5k2 K^/K^a t different time points (P60, P90, P120, P150) were assessed (n = 7; mean ± SEM). Significantly reduced wave I amplitudes at 24KHz were detected at P120 in Ppip5k2 K^/K^m ice ( Ã p < 0.01). At P150 both the homozygous and heterozygous Ppip5k2 K^m ice had significantly reduced wave I amplitudes at 24KHz mice ( ÃÃ p < 0.05). These results suggest impaired late onset auditory neuronal function in PPIP5K2 mutant mice. https://doi.org/10.1371/journal.pgen.1007297.g005 We used the XHMM and CoNIFER methods [40, 41] to call CNV events in WES data. The segregation of the candidate variants in participating member of family PKDF751 was performed through Sanger sequencing. Primers were designed with Primer3 software, and used to amplify exons as well as flanking introns and untranslated regions (S4 Table).

Ppip5k2 gene-trap mouse model
All experiments were approved by the Animal Care and Use Committees at the University of Maryland, School of Medicine in accordance with the National Institutes of Health (NIH) Guide for the Care and Use of Laboratory Animals. The B6N(Cg)-Ppip5k2 tm1a(EUCOMM)Wtsi/J strain was generated by KOMP. In this study, we designated this strain as Ppip5k2 K^d ue to its intact kinase domain. The Ppip5k2 K^m ice had a gene-trap cassette downstream of exon 13, which is expected to cause premature truncation.
In some experiments, the wild type (WT) and mutant PPIP5K2 constructs were expressed in PPIP5K -/-HEK293 cells that had been radiolabeled by growing them in media containing [ 3 H]-inositol. Cellular levels of [ 3 H]IP7 and [ 3 H]IP8 were determined by HPLC [12].
For the purification of WT and mutant PPIP5K2 proteins, the HEK293 cells hosting the various constructs were harvested 16-20 hr after transfection and lysed in ice-cold buffer containing 10 mM HEPES, 130 mM NaCl, 1% Triton X-100, 10 mM NaF, 10 mM Na 2 HPO 4 , 10 mM Na pyrophosphate and protease inhibitor cocktail (Roche). All subsequent steps were performed on ice in an anaerobic chamber (Bactron CAT180). The FLAG-tagged PPIP5K2 proteins were immunopurified using FLAG M2 affinity gel (Sigma). Purified FLAG-PPIP5K2 proteins were analyzed by SDS-PAGE and stained with Coomassie Blue. Densitometry analysis (ImageJ) was used to calculate the amounts of PPIP5K2 loaded onto the gel, by reference to additional lanes loaded with known amounts of the human recombinant PPIP5K2 kinase domain.

Phosphatase and kinase assays
All enzyme assays were performed using physiologically relevant substrate concentrations at 37˚C; the extent of substrate conversion was up to 30%. IP8 phosphatase activities were Confocal images of the organ of Corti, labeled with antisera against MYO7A, a hair cell marker, at the apical, middle and basal turns of the cochlea in Ppip5k2 +/+ , Ppip5k2 +/K^, and Ppip5k2 K^/K^m ice at P150. In WT mice, the position of dead OHCs are marked with asterisks ( Ã ), observed only in the basal turn. In contrast, significantly higher number of OHCs were missing in Ppip5k2 +/K^, and Ppip5k2 K^/K^m ice, which correlates with the high-frequency HL observed in these mice. Scale bar: 20μm. (B) Inner hair cells and outer hair cells were quantified in Ppip5k2 +/+ , Ppip5k2 +/K^, and Ppip5k2 K^/K^m ice at P150. For quantification purposes, organ of Corti were isolated from three different mice of each genotype, and the hair cells were counted in 200μm stretches of apical, middle and basal coil regions. No significant difference was observed for the number of IHCs. Whereas in the case of OHC number (mean ± SEM), a statistically significant ( Ã p<0.05, ÃÃ p<0.01) decrease in OHCs were observed in the basal coil of Ppip5k2 +/K^, and Ppip5k2 K^/K^m ice, which coincide with the high-frequency HL observed at this age in both genotypes.

Confocal imaging
Inner ear tissues were harvested from temporal bones of the WT and Ppip5k2 K^m ice, fixed and processed for immunostaining as previously described

Assessment of hearing function
Briefly, mice were anesthetized with an IP injection of ketamine/xylazine mix at 80-100 mg/kg and 10-12.5 mg/kg, body weight respectively. ABR recordings were performed using a TDT RZ6/BioSigRZ system along with RA4PA preamplifier/digitizer. The system was calibrated between 4-40KHz frequency range from 90 to 5 db using the ACO Pacific type 7016 microphone. For simultaneous recording from both ears, two open filed speakers were placed 10 cm apart from each ear. Two electrodes were placed behind each ear. The reference was placed at the base of the skull, whereas the ground was placed at the base of the tail. Mice were tested at broad band clicks and pure tone frequencies of 8, 16, 24, and 32 kHz for ABRs (n = 7 per genotype at each time point). All the stimuli were tested from 90 dB SPL to10 dB SPL in 5 dB decrements, with a total of 512 responses averaged at each level. Results are displayed as a mean and standard error of the mean (SEM). When no ABR signals were observed at 90db, for calculation purposes, mice were assigned a threshold of 95db. Statistical significance was determined using a 2-way ANOVA. A Bonferroni correction was performed to adjust for multiple statistical tests.

RNA extraction and quantitative real-time PCR
The whole inner ear was isolated from P150 WT and Ppip5k2 K^/K^m ice using the RiboPure RNA isolation kit (Life Technologies, Grand Island, NY). cDNA was prepared using an oligo-dT primer and SMARTScribe Reverse Transcriptase enzymes (Clontech, Mountain View, CA). To determine the differential expression of Ppip5k1, SYBRGreen based real-time primers were designed using Integrated DNA Technologies online PrimeTime qPCR assay design tool. The real-time PCR assays were performed in triplicate using an ABI StepOnePlus Real-Time thermal cycler (ABI, Foster City, CA). CT values were normalized using Gapdh as an endogenous control, and fold changes of Ppip5k1 transcripts in different tissues were calculated using 2^-ΔΔCt standard formula. An expression with a 2-fold change and with a P-value less than 0.05 based on a Student's t-test analysis was considered significant.

Web resources
The URLs for data presented herein are as follows: