A human stem cell resource to decipher the biochemical and cellular basis of neurodevelopmental defects in Lowe syndrome

ABSTRACT Human brain development is a complex process where multiple cellular and developmental events are coordinated to generate normal structure and function. Alteration in any of these events can impact brain development, manifesting clinically as neurodevelopmental disorders. Human genetic disorders of lipid metabolism often present with features of altered brain function. Lowe syndrome (LS) is an X-linked recessive disease with features of altered brain function. LS results from mutations in OCRL1, which encodes a phosphoinositide 5-phosphatase enzyme. However, the cellular mechanisms by which loss of OCRL1 leads to brain defects remain unknown. Human brain development involves several cellular and developmental features not conserved in other species and understanding such mechanisms remains a challenge. Rodent models of LS have been generated but failed to recapitulate features of the human disease. Here we describe the generation of human stem cell lines from LS patients. Further, we present biochemical characterization of lipid metabolism in patient cell lines and demonstrate their use as a ‘disease-in-a-dish’ model for understanding the mechanism by which loss of OCRL1 leads to altered cellular and physiological brain development. This article has an associated First Person interview with the first author of the paper.


INTRODUCTION
Phosphoinositides are key regulators of the organization and function of eukaryotic cells (Schink et al., 2016). Of the seven species of phosphoinositides, phosphatidylinositol 4,5bisphosphate [PI(4,5)P 2 ] is the most abundant and is required to regulate multiple sub-cellular processes including membrane turnover, cytoskeletal function and the organization of the plasma membrane (Kolay et al., 2016). PI(4,5)P 2 exerts its control over cellular functions both through binding and allosteric regulation of protein activity and also through its ability to serve as a substrate for phospholipase C (PLC) and Class I PI3K signaling (Katan and Cockcroft, 2020). Therefore, the accurate regulation of PI(4,5)P 2 levels at cellular membranes is critical for normal function. PI(4,5)P 2 levels are regulated through enzymes that regulate its synthesis, the phosphatidylinositol 4-phosphate 5-kinases (PIP5K) (van den Bout and Divecha, 2009) and also by enzymes that catalyze its metabolism. In addition to PLC and Class I PI3K that utilize PI(4,5)P 2 to generate signaling molecules, lipid phosphatases that can dephosphorylate PI(4,5)P 2 have also been described. These include 4-phosphatase enzymes that generate PI5P, but their function in vivo remains unclear (Ungewickell et al., 2005). A large family of 5-phosphatases that can dephosphorylate PI(4,5)P 2 at position 5 to generate phosphatidylinositol 4-phosphate (PI4P) have been described (Ooms et al., 2009). This 5-phosphatase activity is encoded in all major eukaryotic genomes, with mammalian genomes encoding up to ten genes for this family of proteins; many of these gene-products have been linked to human diseases (Ramos et al., 2019). The significance of encoding a single enzyme activity through such a large gene family remains to be understood.
The oculocerebrorenal syndrome of Lowe (OCRL) gene encodes a 901 amino acid inositol polyphosphate 5-phosphatase enzyme, which is able to catalyze the removal of the 5′ phosphate from PI(4,5)P 2 to generate PI4P. OCRL was originally identified as the gene underlying the human inherited disease Lowe syndrome (LS) (Attree et al., 1992). In human patients with LS, sequencing studies have revealed a large diversity of mutations in OCRL, including deleterious missense and nonsense mutations in all of the major domains of the protein including the 5-phosphatase domain, PH, ASH and RhoGAP domain (Staiano et al., 2015). OCRL is widely expressed across many human tissues or organs and at all stages of life . The function of OCRL has been studied by overexpression or depletion in a number of common human cell lines and the protein has been reported to localize to and affect the function of many cellular organelles (Mehta et al., 2014).
LS is a rare (∼1/500,000 males), X-linked, recessive disorder characterized by the triad of congenital cataracts, intellectual or neurodevelopmental impairment and proximal renal tubular dysfunction (https://omim.org/entry/300535) (Bökenkamp and Ludwig, 2016;De Matteis et al., 2017). The brain phenotypes in LS include delayed and impaired cognitive milestones, hypotonia, febrile seizures and hyperechoic changes in the periventricular zone of the cerebral cortex. There are two enigmatic and unresolved observations in relation to the clinical presentation of LS: (i) Although LS is a monogenic disorder, there is substantial variability in the clinical presentation between individual patients, even in those individuals with mutations in OCRL with equivalent molecular consequences (e.g. truncating nonsense mutations prior to the start of the phosphatase domain). For example, while some patients present with severe neurodevelopmental phenotypes, others show relatively mild deficits in brain function. These disparities suggest that in addition to the mutation in OCRL, other changes in the genetic background of the individual may impact the clinical outcome of loss of OCRL. (ii) Although OCRL is widely expressed in human tissues, it remains a mystery as to why only three organs are affected in LS, namely the eye, brain and kidney. One possibility is that the phenotypic changes seen in human patients may arise due to the requirement of OCRL only in specific cell types of the affected organs. Therefore, the relevant cellular changes may only be seen when studying eye, brain or renal tissue during development. Although rodent models of OCRL were generated, they failed to show phenotypes that recapitulate the human disease (Jänne et al., 1998). A zebrafish model of OCRL depletion has been generated that recapitulates some aspects of the human phenotype but there remains a lack of models that allow the brain phenotype to be studied (Ramirez et al., 2012). A limited number of studies have been done on LS fibroblasts and renal biopsies, but there is presently no understanding of how loss of OCRL leads to neurodevelopmental phenotypes. Thus, there is a requirement for a model system in which the cellular and physiological changes in the brain during development can be studied.
One possible route to obtaining a suitable model system arises from the ability to use modern stem cell technology to generate human induced pluripotent stem cells (hiPSC) (Shi et al., 2017) from the somatic tissues of patients with LS. These hiPSC can then be differentiated into specific adult tissues and the trajectory of development along with the cellular and molecular changes in any particular patient derived line can be analyzed (Russo et al., 2015;Zeng et al., 2014). In the context of LS, a limited number of studies have reported individual hiPSC lines derived from patients (Barnes et al., 2018;Hsieh et al., 2018;Liu et al., 2021;Qian et al., 2021). In this study, we present the generation of hiPSC lines and neural derivatives from a family with a unique genetic structure and clinical features that should allow an understanding of the cellular basis of the neurodevelopmental phenotype in LS. We also present a biochemical analysis of phosphoinositide levels and an insight into the biochemical compensation for the loss of the 5-phosphatase activity of OCRL in LS cells.

Generation of hiPSCs from a family with LS
For this study, we selected a family whose genetic structure (Fig. 1Ai) is uniquely suited for the analysis of neurodevelopmental phenotypes in LS (Ahmed et al., 2021). Briefly, in this family, the patients with LS are the children of two female siblings, both of whom carry the identical mutation in OCRL. While all three children show the triad of eye, renal and brain phenotypes characteristic of LS, the brain phenotype of LSPH004 is much more severe than that of the identical twins LSPH002 and LSPH003. To understand the cellular and developmental mechanisms that underlie this neurodevelopmental defect, we generated hiPSC from each of these patients. Peripheral blood mononuclear cells (PBMC) were isolated from each patient, immortalized into lymphoblastoid cells lines (LCLs), which were then reprogrammed to generate hiPSC (Iyer et al., 2018) (Fig. 1Aii). The hiPSC lines so generated did not express the OCRL protein as determined by immunocytochemistry with an antibody to OCRL (Fig. 1B). The hiPSC lines showed expression of pluripotency markers SOX2 and SSEA4 as determined by immunocytochemistry (Fig. 1C) and SSEA4 and OCT4 as determined by single cell quantitative fluorescence activated cell sorting (FACS) analysis ( Fig. S1A-C). We differentiated these hiPSC into embryoid bodies and established their ability to differentiate into each of the three germ layers by detecting the expression of transcripts characteristic of each layer (SOX1, Nestin-ectoderm, Nodal-mesoderm and GATA4endoderm) (Fig. 1D). All hiPSC lines were determined to be of a normal karyotype ( Fig. 1E; Fig. S2D,E) and short tandem repeat analysis was used to determine and track the identity of each cell line (Table S1). These hiPSC lines offer a unique resource from which tissue specific differentiation, for example into brain tissue can be carried out; a comparison of the cellular and molecular differences between control and patient derived hiPSC lines during brain development can provide important insights into how loss of OCRL results in altered neurodevelopment.

Generation of neural stem cells (NSC) from hiPSC
During the development of the brain, a key step in the conversion of pluripotent, early embryonic stem cells into brain cells is the formation of NSC, which then both divide and differentiate to generate the different cell types of the brain. Thus, with the exception of microglia, NSC can be differentiated into all cell types in the brain. Since the neurodevelopmental phenotype was strongest in LSPH004, we generated NSC to understand the brain development phenotype of LS. As a control, we generated NSC from D149 hiPSC (Iyer et al., 2018) that was originally derived from an unaffected individual of similar Indian population background (Fig. 1Aii). The NSC so generated were characterized by confirming expression of established NSC protein markers such as Nestin, SOX1, SOX2, PAX6 and Musashi-1 along with the proliferation marker Ki-67 ( Fig. 2A). They were also confirmed to be karyotypically normal (Fig. S3A,B). Lastly, these NSC were differentiated into cultures of forebrain cortical neurons. The neurons so generated showed the characteristic morphology and molecular markers of neuronal development including MAP2, DCX and Synapsin-1 (Fig. 2B) as previously reported (Sharma et al., 2020). Western blot analysis of protein extracts from these lysates revealed the band corresponding to OCRL in wild-type neurons and absent in LSPH004 (Fig. 2C). To test the physiological status of these neurons, we monitored them for the presence of intracellular calcium transients [Ca 2+ ] i , a characteristic feature associated with neuronal development (Rosenberg and Spitzer, 2011). We found that in 30 days in vitro (DIV) cultures, robust [Ca 2+ ] i transients were observed in neurons derived from both D149 (Fig. 2Di) and LSPH004 (Fig. 2Dii). In addition, we also monitored the development of electrical activity in the differentiating neurons as a function of age in vitro using whole cell, patch clamp electrophysiology in both D149 and LSPH004. For example, in ten DIV cultures of D149, abortive action potentials (AP) were noted (Fig. 2Ei), by 20 DIV single AP were seen (Fig. 2Eii), by 30 DIV multiple AP were noted (Fig. 2Eiii) and by 40 DIV repetitive AP firing was observed (Fig. 2Eiv). Abortive and single action potentials are indicative of the incomplete physiological maturation of a neuron. In this study we recorded AP from ca. 130 individual cells; at early stages of differentiation (DIV 14-25) ca. 79% of cells showed an abortive AP, 19% showed single AP and no multiple or repetitive AP were noted. As the cultures increased in age, this pattern shifted to progressively more mature AP; by DIV 35-50, 31% cells showed multiple AP and 29% repetitive ones while the proportion of cells showing abortive or single AP reduced substantially (Fig. 2F). Thus, in this model system, we are able to monitor the progressive electrophysiological development of neurons.

Tools for controlled expression of proteins in NSC
In order to understand the cellular basis of the neurodevelopmental defect in LS, it is essential to monitor in real time, ongoing cellular and molecular processes using protein-based reporters such as those used to monitor phosphoinositide turnover at specific cellular membranes (Hammond and Balla, 2015). Likewise, in order to understand the function of specific domains of OCRL including the 5-phosphatase domain, it will be necessary to achieve carefully controlled reconstitution of OCRL variants in LS cells during development. Current methods of protein expression in neural cells, such as lentiviral mediated transduction of transgenes, result in variable levels of expression between cell lines and experiments. One strategy for the controlled expression of proteins in stem cells is integration of genetic constructs at specific locations in the genome that have been found to be suitable for expression of a transgene/ biosensor without any major adverse consequences due to insertion (Pei et al., 2015). Briefly, transcription activator-like effector nucleases (TALENs) are used to insert a Lox cassette at adenoassociated virus integration site 1 (AAVS1) [referred to as Genomic safe harbor sites (GSH)] locus on chromosome 19, thus generating a master cell line for subsequent insertion of transgenes at this location. Transgenes of interest, under suitable promoters, can be inserted into this safe harbor site by Lox-P/Cre-recombinase mediated cassette exchange (RMCE) (Fig. 3A). The advantage of this approach is that there is no rapid loss of transgene expression after transfection and no variation in copy number of transgenes between experiments or cell lines as might occur with transient transfection or lentiviral transduction of transgenes.
Using this approach, we inserted GSHs into D149 and LSPH004 NSCs. Insertion of GSHs was monitored by observing GFP expression and using a junction PCR (Fig. 3B) and the NSC line so generated continued to express characteristic protein markers (Fig. 3C). RMCE was used to insert a Lox cassette into GSHs of D149 (control) and LSPH004 ( patient) NSCs; this cassette included the cDNA for a protein probe for PI(4,5)P 2 , the PH domain of PLCδ fused to mCherry (mCherry:: PH-PLCδ) (Várnai and Balla, 1998). In this cassette, the mCherry::PH-PLCδ is expressed under a CAG promoter that is active in NSC. NSC in which the mCherry::PH-PLCδ transgene is recombined into the GSH show mCherry fluorescence, NSC in which only one of two copies of the GSH have undergone RMCE show both mCherry and GFP fluorescence and cells in which no recombination has taken place show only GFP fluorescence (Fig. 3D). We purified NSC with mCherry fluorescence only using FACS (Fig. 3E). Western blot analysis on these purified cells showed that they express a protein with a M r corresponding to that of mCherry::PH-PLCδ (Fig. 3F).
Impact of OCRL depletion on total cellular PI(4,5)P 2 Since OCRL is a PI(4,5)P 2 5-phosphatase, it is expected that cellular PI(4,5)P 2 levels might be elevated and PI4P levels reduced in LS patient cells. To test this, we extracted total lipids from hiPSC and NSC of both control and patient derived LSPH004 cells. PIP 2 levels were measured using liquid chromatography coupled with mass spectrometry (LCMS) from whole cell lysates (Sharma et al., 2019); although this method cannot distinguish between positional isomers of PIP 2 , the majority is expected to be PI(4,5)P 2 . In hiPSCs, a significant increase was seen in the total PIP (Fig. 4A) and PIP 2 (Fig. 4B) mass in LSPH004 compared to control (Fig. 4A,B). In experiments with NSCs, we compared D149 with LSPH004 and found no difference in the total PIP mass (Fig. 4C) or PIP 2 mass ( Fig. 4D) between these two lines. Thus, loss of OCRL results in a modest change in PI4P and PI(4,5)P 2 mass in LS patient derived hiPSC.
Elevated plasma membrane PI(4,5)P 2 levels in Lowe NSC Since we did not observe a significant change in the total mass of PI(4,5)P 2 in the LS patient NSC, we wondered whether the OCRL enzyme might control a relatively small but functional pool of PI(4,5)P 2 at a specific endomembrane in NSC and changes in this small pool might not be reflected in measurements of total PI(4,5)P 2 mass measurements. In cells, a key sub-cellular membrane where PI(4,5)P 2 is enriched is the plasma membrane. To measure plasma membrane PI(4,5)P 2 levels, we used the reporter lines expressing mCherry::PH-PLC-δ. Live cell imaging was performed on D149 and LSPH004 NSC expressing mCherry::PH-PLC-δ. We estimated PI(4,5)P 2 levels at the plasma membrane by imaging the biosensor expressing cells in a monolayer and calculating the plasma membrane to cytosolic (PM/Cyt) fluorescence ratio of the probe (Fig. 4E). Quantification of these data revealed that the PI(4,5)P 2 levels were significantly higher at the plasma membrane of the patient cell line LSPH004 compared to D149 (Fig. 4F). Thus, loss of OCRL alters the pool of PI(4,5)P 2 at the plasma membrane in NSC.

Compensatory mechanisms for loss of OCRL function
OCRL is part of a large family of lipid 5-phosphatases in the human genome (Ramos et al., 2019) and the function of this gene family in neurodevelopment has not been studied. Loss of OCRL function in LS patients might result in compensatory changes in the expression of other 5-phosphatase family members during brain development, thus leading to modest or no changes in PI4P and PI(4,5)P 2 levels.
To determine the expression pattern of these phosphatases during neural development, we performed qRT-PCR analysis for all the ten 5-phosphatases encoded in the human genome (Ramos et al., 2019) in hiPSC, NSC and 30 DIV neuronal cultures of D149. This analysis revealed an interesting and variable pattern of expression for each gene at these specific stages of neural differentiation in vitro. While some 5-phosphatases such as SYNJ1, SYNJ2, INPP5B and INPPL1 were expressed at similar levels across all three stages, INPP5D appeared downregulated during neuronal differentiation. By contrast, a set of 5-phosphatases including OCRL, INPP5F, INPP5K, INPP5E and INPP5J all showed upregulation during neuronal differentiation (Fig. 5A). We then compared expression of all ten phosphatases in each of the three developmental stages hiPSC, NSC and 30 DIV neurons between D149 and LSPH004. At the hiPSC stage, five of ten 5-phosphatases that we assayed were upregulated in the patient line LSPH004; these were INPP5D, INPP5E, INPP5F, INPP5K, INPPL1. By contrast, SYNJ1 and INPP5J were downregulated; as expected, OCRL was significantly downregulated in the patient hiPSC (Fig. 5B). In NSC, we observed five out of ten 5-phosphatases were upregulated in the patient line namely, INPP5J, INPP5E, INPP5F, SYNJ2, INPPL1 and the only 5-phosphatase downregulated was OCRL (Fig. 5C). We also differentiated the NSC into neurons and compared expression levels of the phosphatases at 30 DIV between D149 and LSPH004; this revealed that except for a modest downregulation of INPP5K, there were no compensatory changes in LSPH004 (Fig. 5D). Thus, loss of OCRL results in distinctive patterns of compensatory changes in 5phosphatase gene expression at various stages of neurodevelopment ( Fig. 5E).

Transcriptomic changes in LS cells
To characterize the gene expression changes resulting from OCRL loss of function that may lead to neurodevelopmental defects, we performed transcriptomic analysis from hiPSC and NSC of D149 and LSPH004. At the iPSC stage, corresponding to the earliest stages of human embryonic development, we found ca. 475 genes upregulated and ca. 400 genes downregulated (Fig. 6A). Analysis of this set of altered genes using Gene Ontology (GO) revealed a number of categories of highly enriched downregulated genes suggestive of altered brain development (GO terms: generation of neurons, neurogenesis, neuron differentiation, nervous system development) (Fig. 6C,D; Table S2). These included SLITRK1, DCX, MAP2, BRINP1, CEP290, PAX7, NCAM1, SEMA6D and CNTN2.
Likewise, a comparative analysis of D149 and LSPH004 NSC revealed more than 750 downregulated and 250 upregulated genes (Fig. 6B). GO analysis of this gene set also revealed strong enrichment of GO terms suggesting altered neural development with the most highly significant GO terms being neurogenesis and nervous system development and function ( Fig. 6E; Table S3). These include SEMA3F, PLXND1, SOX8, LRRC7, HAPLN2, NTRK2, NRG1 and SYT4. These findings suggest a role for OCRL in the development of the brain, starting with the earliest cellular stages of embryogenesis, that can be analyzed in this experimental system.

DISCUSSION
Although several sub-cellular functions have been described for OCRL (De Matteis et al., 2017;Mehta et al., 2014), the mechanism by which mutations in this gene result in human neurodevelopmental defects remains unknown. A particular challenge in understanding the brain phenotype in LS patients arises from the inability to obtain conventional biopsy samples from the human brain. This is in contrast to a limited number of studies where renal biopsies (De Leo et al., 2016) and skin samples (Wenk et al., 2003) have been used to address some specific aspects of the LS phenotype. In addition, since LS is neurodevelopmental in origin, understanding brain phenotypes present at birth requires the ability to study the cellular and physiological changes during brain development. To achieve this goal, it is necessary to have a model in which the development of human brain cells can be studied in vitro. A suitable model must recapitulate key aspects of brain development in vitro and also the key biochemical effects seen in LS patient cells.
We have identified a family with LS (Ahmed et al., 2021) and in this study, report the generation of hiPSC from patients in this family as a resource to analyze the cellular and molecular basis of brain development in LS. Since OCRL encodes a PI(4,5)P 2 phosphatase, it is expected that loss-of-function mutations in this gene should result in elevation of PI(4,5)P 2 levels and a drop in PI4P levels. Our measurements of the mass of PIP 2 and PIP in hiPSC cultures of the patient lines showed an elevation of total PIP 2 and PIP in LSPH004 cells. This parallels observations on cultured human fibroblast cell lines from LS patients that have been previously reported (Wenk et al., 2003). In this study, we also engineered the LSPH004 line to express a fluorescent reporter for plasma membrane PI(4,5)P 2 levels and found using this approach that the levels of this lipid at the plasma membrane of the patient line were elevated in comparison to control. Thus, the in vitro model system presented here recapitulates the key biochemical defect in phosphoinositide metabolism reported in LS patient tissues.
We differentiated these hiPSC lines using protocols based on the known principles of developmental neurobiology to generate 2D cultures of human neural tissue, primarily neurons in vitro. This developmental process in vitro recapitulates key aspects of brain development in vivo (Mertens et al., 2016). These hiPSC derived neural cultures show progressive increase of neural activity as a function of the age of the culture (DIV) (Fig. 2D,E). A comparison of neural activity development from LS patient derived cell lines reported here compared to those from control lines will provide an insight into the altered physiological development in brain tissue of LS patients. During metazoan development, gene expression is a key process that determines cell fate specification and differentiation. To understand gene expression changes that underlie the neurodevelopment phenotype in LS, we compared transcriptomes between control and LSPH004 at both the iPSC and NSC stages. Large differences in transcriptome were seen at both of these developmental stages (Fig. 6) revealing substantial changes in the expression of genes annotated as having roles related to nervous system development. We have presented an in-principle analysis of only LSPH004 in comparison to D149. However, it is to be noted that individual hiPSC lines for the same disease can have phenotypic differences, presumably due to reasons of genetic heterogeneity or other factors such as arising during line generation.  Table depicting the proportion of neurons that showed abortive, single, repetitive and multiple action potentials as a function of age. Columns depict the age of the neurons as DIV. Rows depict the type of action potential. Each number shown is the percentage of cells recorded at that age, which showed a specific type of action potential. Therefore, any conclusions on the phenotypes seen at the transcriptomic, cellular or physiological level should be validated carefully to link them to the mutation in the proposed gene under study (in this case OCRL). For example, phenotypes observed in LSPH004 in comparison to D149 should also be evaluated in additional ways that could include (i) comparing phenotypes of LSPH004 to that of LSPH002 and LSPH003 that carry the same mutation in OCRL (ii) comparing LSPH004 phenotypes to those in D149 hiPSC edited to introduce the OCRL mutation seen in the patient line (iii) reconstitution of LSPH004 with a wild-type copy of OCRL to rescue any phenotypes seen in this patient line. Phenotypes satisfying one or more of the above conditions will validate patient line phenotypes arising from OCRL loss of function. The use of the hiPSC lines and genome engineering tools described here are likely to be valuable in understanding the mechanism by which loss of OCRL leads to altered brain development. Altogether, the model system described here allows an analysis of the physiological development in LS patients as well as the biochemical and molecular correlates of this process.
In addition to neuronal differentiation, the LS patient derived hiPSCs can be differentiated into glia to study the contribution of OCRL function in these non-neuronal cell types to the altered brain development in LS patients. Hyperechoic periventricular lesions have been described in MRI scans of patients with LS with features suggestive of enhanced gliosis (Sener, 2004) and evidence suggestive of enhanced gliosis in a fish model of LS has been reported (Ramirez et al., 2012). However, its relevance in the brain tissue of human LS patients remains unknown and the hiPSC lines generated here will also allow this question to be studied. In addition, the use of co-cultures of one neural cell type (e.g. neurons) from LS patients with glial cells derived from controls and vice versa (Nadadhur et al., 2019) will allow an assessment of the cell autonomous and non-cell autonomous function of OCRL in supporting normal brain development. The developing brain is a complex 3D tissue and 2D cultures may not recapitulate all aspects of the cell-cell interactions or cortical layer formation that impact normal brain development. However, hiPSC can be used to generate 3D organoids of the developing brain (Di Lullo and Kriegstein, 2017); the use of the LS patient lines described in this study will also allow the role of OCRL in supporting the 3D architecture during brain development to be investigated. Finally, a recent study has proposed the use of the phosphoinositide 3-kinase inhibitor, alpelisib, for correcting renal defects in a mouse model of LS (Berquez et al., 2020). The availability of the hiPSC lines described in this study can help evaluate the effectiveness of such inhibitors in alleviating the phenotypes of LS in brain tissue.
It has been noted that although LS is a monogenic disorder, there is considerable variability in the brain phenotypes of individual patients. It has been proposed that such variability in clinical features between patients may arise from the impact of background mutations in the genomes of individual patients carrying functionally equivalent mutations in OCRL. However, the cellular and developmental basis of this has not been tested. The unique genetic structure and clinical features of the LS family we have studied (Ahmed et al., 2021) provides a unique opportunity to address this scientific question. The three patients in this study, LSPH002, LSPH003 and LSPH004 vary clinically in their neurodevelopmental phenotype with LSPH004 showing much severe brain phenotype than LSPH002 and LSPH003 who are identical twins. In this study, we have generated hiPSC from all three patients; these cell lines carry the identical mutation in OCRL in the background of the genome of the individual patient from whom they were derived. By differentiating these hiPSC into neural tissue and comparing LSPH004 with LSPH002 and LSPH003, it is likely that we will discover the cellular and developmental correlates of the variable brain phenotype. Since the whole genome sequence of each line has been determined, it is also possible to experimentally test the importance of specific variants in the patient genome by CRISPR/Cas9 genome editing of specific variants (Paquet et al., 2016) followed by phenotypic analysis.
Although we observed changes in the levels of PIP and PIP 2 in patient derived hiPSC from LS patient, several key aspects of altered phosphoinositide biochemistry in LS patient cells remain to be understood: (i) Changes were relatively modest and the changes in PIP and PIP 2 mass seen in hiPSC were no longer evident at the NSC stage. These observations indicate a plasticity in the control of PI(4,5)P 2 levels, perhaps by other enzymes that may also regulate PI(4,5)P 2 levels that needs to be understood. (ii) In addition to elevated PIP 2 levels in LS hiPSC, PIP levels were also elevated; this is unexpected though it has also previously been reported in fibroblasts from LS patients (Wenk et al., 2003). (iii) It is unclear if the elevated PIP 2 levels, PIP levels or both lead to cellular effects leading to patient phenotypes. In order to address these questions, the system described here offers many advantages. Our genome engineering approach using the GSH in these lines offers the opportunity for controlled expression of enzymes to modulate the levels of phosphoinositides acutely (Idevall-Hagren and De Camilli, 2014;Varnai et al., 2006) to test specific hypotheses related to the role of individual lipids in altered cellular and developmental phenotypes in developing neural cells. One likely mechanism underlying plasticity in the control of PI(4,5)P 2 levels in LS patient cells is changes in the expression of other genes encoding members of the 5-phosphatase family of enzymes. Indeed, RT-PCR analysis of the ten members of the 5-phosphatase in LSPH004 revealed that there were distinctive patterns of up-or downregulation in the expression of the 5-phosphatase gene family members in wild-type cells during development. Such patterns of 5-phosphatase expression and compensatory changes of these in LS patient cells during development may underlie the specific nature of the neurodevelopmental defects in LS patients. In summary, the  Plot showing total PIP 2 levels using LCMS in whole cell lipid extract from hiPSCs of Control (D149, n=9) and Patient lines (LSPH004, n=10). X-axis denoting samples and Y-axis represents normalized PIP 2 levels. Statistical test: Two-tailed unpaired t-test with Welch's correction. Whiskers at minimum and maximum values and a line at the median. (C) Box and Whisker Plot showing total PIP levels using LCMS in whole cell lipid extract from NSCs of Control (D149, n=12) and Patient line (LSPH004, n=12). X-axis denoting samples and Y-axis represents normalized PIP levels. Statistical test: Two-tailed unpaired t-test with Welch's correction. Whiskers at minimum and maximum values and a line at the median. (D) Total PIP 2 levels using LCMS in whole cell lipid extract from NSCs of Control (D149, n=12) and Patient line (LSPH004, n=12). X-axis denoting samples and Y-axis represents normalized PIP levels. Statistical test: Two-tailed unpaired t-test with Welch's correction. Whiskers at minimum and maximum values and a line at the median. (E) Representative confocal maximum z-projections of mCherry::PH-PLCδ expressing PIP 2 reporter lines used to estimate plasma membrane/cytosolic PIP 2 probe fluorescence ratio (PM/Cyt). Enlarged insert shows mCherry::PH-PLCδ biosensor localization to the plasma membrane in a single cell. (F) Dot plot denoting quantification of PIP 2 levels using mCherry::PH-PLCδ biosensor in Control (D149) and Patient (LSPH004) PIP 2 reporter neural stem cell lines. X-axis denotes samples, Y-axis represents plasma membrane/cytosolic fluorescence ratio (PM/Cyt) of mCherry::PH-PLCδ biosensor. Each dot represents PM/Cyt obtained from a cell. Statistical test: Two-tailed unpaired t-test with Welch's correction. ****P-value <0.0001. Error bars: s.e.m. hiPSC resources and their engineered derivatives described here offer powerful tools for understanding the regulation of PI(4,5)P 2 to PI4P balance in the developing nervous system by OCRL and the mechanism by which loss of this activity leads to neurodevelopmental defects.

Mycoplasma testing
Mycoplasma contamination was checked using spent media from subconfluent hiPSC and NSC dishes after 48 h in culture. MycoAlert™ (Lonza, #LT07-418) was used per the manufacturer's protocol.

Karyotyping
Overall chromosomal integrity of hiPSCs and NSCs was confirmed by karyotyping. For metaphase preparation, cells were arrested in log phase by treating with 0.1 µg/mL Colcemid™ (Gibco, #15212-012) treatment for 45 min at 37°C. Cells were harvested in fresh Carnoy's fixative (Methanol: Glacial Acetic Acid at 3:1) and G-banding karyotype analysis performed at a National Accreditation Board for Testing and Calibration Laboratories, India (NABL) accredited facility.

Generation of hiPSC lines
Blood was drawn from donors after informed consent and under aseptic conditions following IRB regulations. A peripheral blood mononuclear cell fraction was obtained by density gradient centrifugation and transformed into LCL by Epstein Barr virus (EBV) transformation using previously established protocol (Hui-Yuen et al., 2011). These LCLs were reprogrammed by electroporation of plasmids containing Yamanaka factors to generate hiPSCs as previously described (Iyer et al., 2018). Onfeeder hiPSC cultures were gradually transitioned to feeder-free conditions by weaning off from the standard HuES media to Essential 8™ medium (Gibco, #A1517001) and maintained on hESC-qualified Matrigel coated surface. hiPSCs were frozen at a density of 1×10 6 in 500 µl PSC cryomix (Gibco, #A26444-01). The plasmid footprint in feeder-free iPSC lines LSPH002 ( passage 7), LSPH003 ( passage 7) and LSPH004 ( passage 8) was tested using PCR analysis with relevant primer sets (Table S5), with 3% DMSO as a PCR additive. Plasmid pCXLE-hUL was used as a positive control while iPSC NIH5 was used as negative control. To ascertain differentiation potential, embryoid bodies (EB) were generated by transferring hiPSCs to a non-adherent surface in standard HuES media without bFGF for 48 h, and cDNA probed for primers specific for the three lineages.

Generation of NSCs
NSCs were generated as previously described (Mukherjee et al., 2019) with slight modifications. Briefly, hiPSCs were differentiated to form embryoid bodies in E6 medium (Gibco, #A1516401). Primary neural rosettes formed by this method were selected and manually passaged to obtain secondary and tertiary rosettes that were eventually triturated and plated as NSC monolayer in Neural Expansion Medium (NEM). As a measure to eliminate any non-NSC cells and obtain a reliable homogeneous NSC culture, the generated NSCs were subjected to CD133+ selection as previously described via FACS (Peh et al., 2009). A sub-confluent NSC culture maintained in NEM was enzymatically dissociated using Stempro Accutase (Gibco, #A11105-01) and following a wash in PBS, cells were immunolabelled using CD133-PE conjugated primary antibody (Abcam, #ab253271) at 10 µl antibody/million cells and incubated at room temperature for 30 min in dark. Cells were washed with PBS and resuspended in 1 mL sorting media [(1x DMEM/F12-without phenol red (Gibco, #21041025), 1% FBS (Gibco, #16000-044)+Penicillin-Streptomycin (Gibco, #15140-122)] keeping the concentration at 2 million cells per mL to obtain efficient sorting. The cells were sorted using FACS-Aria-III instrument (BD Biosciences) (Fig. S1D,E). Forward and side scatter parameters were adjusted so as to eliminate cell clumps and debris. Cells with highest fluorescence intensity were collected and plated on Matrigel at a concentration of 0.5 million cells per well of a 12-well tissue culture plate in NEM and expanded. Cells were periodically checked for bacterial or mycoplasma contamination. NSCs were further characterized by immunofluorescence and neuronal differentiation.

Reporter line generation
For generation of D149 and LSPH004 NSC safe harbor lines, cells growing in NEM without Penicillin-Streptomycin were harvested by enzymatic dissociation using Stempro Accutase for 5 min at 37°C. The cells were spun down at 1000 g for 3 min to remove Accutase post neutralization at room temperature. After a wash with PBS, cells were incubated in buffer-R (included in Invitrogen, #MPK1025) with 2 μg/μL each of plasmids: AAVS-1-TALEN-R (ZYP017), AAVS-1-TALEN-L (ZYP018) and ZYP037 (Pei et al., 2015) for 3 min at a concentration of 1 million cells per hit. The cell-plasmid mixture was electroplated using the Neon Electroporation System (Invitrogen, USA. #MPK1025) as per the manufacturer's protocol under the following standardized parameters (1100 volts, width-20, pulses-2) and plated in 1 well of a 12-well plate coated with Matrigel containing pre-warmed NEM without Penicillin-Streptomycin. A complete media change was performed 4 h post electroporation to remove dead cells and debris, which could potentially cause cytotoxicity in the culture if not removed. After 24 h, cells were checked under an epifluorescence microscope for expression of GFP and upon confirmation, the cells were subjected to Puromycin (Gibco, #A1113803) selection at a concentration of 0.4 μg/ml for 7-9 days with daily media change. After a week, large GFP positive colonies were visible in the plate, at this point the culture was subjected to FACS to remove any GFP-negative cell. The mixed cell culture was enzymatically dissociated using Stempro Accutase, washed once with PBS and resuspended in sorting media [(1x DMEM/F12-without phenol red (Gibco, #21041025), 1% FBS+Penicillin-Streptomycin] at a concentration of 2 million cells per ml. The cells were sorted using FACS-Aria-III (BD Biosciences). Forward and side scatter parameters were adjusted so as to eliminate cell clumps and debris. The gating parameter threshold for GFP-positive cells was set using non-electroporated control NSCs. The cells were collected in NEM and plated in pre-incubated Matrigel coated plate at a concentration of 0.5 million cells per well of a 12-well tissue culture plate in NEM and expanded for cryopreservation.
The Safe Harbor NSCs; D149-GSH and LSPH004-GSH were enzymatically dissociated using Stempro Accutase for 5 min at 37°C. The cells were spun down at 1000 g for 3 min at room temperature to remove Accutase. After a wash with PBS they were resuspended in Buffer-R (included in Invitrogen, #MPK1025) with 2 μg/μl each of plasmids: LoxP-CAG-mCherry::PH-PLCδ (ZYP070-PIP 2R ) and Cre-Recombinase (ZYP073) at a concentration of 1 million cells per hit of electroporation, which was performed using Neon Electroporation System (Invitrogen, USA. #MPK1025) as per the manufacturer's protocol under the following standardized parameters (1100 volts, width-20, pulses-2) and plated in 1 well of 12-well plate coated with Matrigel containing pre warmed NEM without Penicillin-Streptomycin. A complete media change was performed 4 h post electroporation to remove dead cells and debris which could potentially cause cytotoxicity in the culture if not removed. After 24 h, cells were checked under epifluorescence microscope for expression of mCherry and upon confirmation, the cells were subjected to G418 (Gibco, #10131035) selection at a concentration of 400 μg/ml for 15-20 days with daily media change. After about 2 weeks under selection, distinct mCherry positive colonies were observed, at this point the culture was subjected to FACS to remove GFP-positive cells. The cells were sorted at a concentration of 2 million cells per ml using FACS-Aria-III (BD Biosciences). Forward and side scatter parameters were adjusted so as to eliminate cell clumps and debris. The gating parameter threshold for mCherry-positive cells was set using GFP-positive and non-electroporated control NSCs. Post FACS cells were collected in NEM and plated in preincubated Matrigel coated plate at a concentration of 0.5 million cells per well of 12-well tissue culture plate in NEM and expanded to freeze down additional stock vials.

Isolation of RNA and cDNA synthesis
Total RNA was extracted from well characterized hiPSCs and NSCs of patient and control lines using TRIzol (Ambion, Life Technologies, #15596018) as per the manufacturer's protocol in six biological replicates and quantified using a Nanodrop 1000 spectrophotometer (Thermo Fisher Scientific). Following treatment with 1 U of DNase I (amplification grade, Thermo Fisher Scientific, #18068-015), 1 μg of the RNA from each replicate was used for cDNA synthesis in a reaction mixture containing 10 mM DTT and 40 U of RNase inhibitor (RNaseOUT, Thermo Fisher Scientific, #10777-019). The reaction mixture of 45.5 μl was incubated at 37°C for 30 min followed by heat inactivation at 70°C for 10 min, following which 200 U of Superscript II Reverse Transcriptase (Invitrogen, #18064-014) was added to the reaction volume along with 2.5 μM of random hexamers, and 0.5 mM of dNTPs making the final volume to 50 μl. The reaction was then incubated at 25°C for 10 min, followed by 42°C for 60 min and then heat inactivated at 70°C for 10 min on ProFlex PCR Systems (Life Technologies).

Real-Time Quantitative PCR
The primers used for qRT-PCR were designed using Primer-BLAST, NCBI (https://www.ncbi.nlm.nih.gov/tools/primer-blast/). With a condition of spanning exon-exon junction of their respective genes, their details are provided in supplementary Table 5. Real-Time qRT-PCR was performed in a volume of 10 μl with Power SYBR Green Master mix (Applied Biosystems, #4367659) on an Applied Biosystems ViiA7 system. It was performed with technical triplicates from the patient and control lines with primers for genes of interest and GAPDH was used as a house-keeping gene control. A no-reverse transcriptase control was also set up. The reaction was run under the following conditions: 50°C for 2 min, 95°C for 10 min, followed by 40 cycles of 95°C for 30 s (denaturation), 60°C for 30 s (annealing) and 72°C for 45 s (extension). The C t values obtained for individual genes were normalized to those of GAPDH from the same sample. The relative expression levels were calculated using ΔC t method, whereas the fold change between patient and control was calculated using ΔΔC t method.

Junction PCR
Junction PCR was performed using the genomic DNA of the edited cells, extracted using QIAamp DNA Mini Kit (Qiagen, #51304) using the manufacturer's protocol and quantified using a Nanodrop 1000 spectrophotometer (Thermo Fisher Scientific). Insertion of the ZYP037 donor template was confirmed using two pairs of primers; one of the primers in each pair annealed outside the region spanned by the homology arm, in the AAVS1 locus, while the other annealed within the inserted template, to avoid false detection of residual plasmid if any. 50 ng of genomic DNA was used as the template in 20 μl PCR reaction consisting of Phusion Pol (0.02 U/µl), 25 mM DNTPs, 5X HF-buffer and appropriate primers. The reaction was run as a touchdown PCR in two steps as follows: 98°C for 3 min (initial denaturation), step-1 [98°C for 20 s-denaturation, 66°C to 61°C (−0.5°C per cycle) for 30 s-annealing, 72°C for 2 min-extension] ×10 cycles, step-2 [98°C for 20 s-denaturation, 55°C for 30 s-annealing, 72°C for 2 min-extension] ×25 cycles, 72°C for 5 min (final extension). List of primers and antibodies used in this study is provided in Table S5.
Immunocytochemistry hiPSCs were characterized for pluripotency markers SOX2 and SSEA4 using the PSC immunocytochemistry kit (Invitrogen, #A24881) as per the manufacturer's instructions.
NSC cultures were fixed using 4% formaldehyde in phosphate-buffered saline (PBS) for 20 min, permeabilized using 0.1% Tx-100 for 5 min and incubated at room temperature for 1 h in a blocking solution of 5% BSA in PBS. Primary antibodies at respective dilutions were added and incubated overnight at 4°C in blocking solution, followed by incubation with secondary antibodies in blocking solution (Invitrogen) for 1 h.
Confocal images were recorded by collecting a range of z-stack using an Olympus FV 3000 confocal microscope. Epifluorescence images were captured using EVOS ® FL Cell Imaging System (Thermo Fisher Scientific). The image stack was merged using Z-project (maximum intensity projection) function using ImageJ (National Institute of Health, USA, http://imagej.nih.gov/ij).

Quantitative FACS analysis
hiPSCs were characterized quantitatively for pluripotency markers OCT4 and SSEA4 by flow cytometry. Briefly, the cells were detached using Stempro Accutase for 5 min at 37°C and washed with standard HuES media. The cells were then washed once with PBS and fixed with 4% PFA solution at a cell density of 1×10 6 cells per ml at room temperature for 15 min. The cells were then permeabilized with 0.1% v/v Triton X-100 in HBSS for 15 min at room temperature, washed once with HBSS and incubated with 1% BSA in HBSS for 30 min. Following this, the cells were incubated with primary antibodies in 1% BSA in combination at the following dilutions, Rabbit Anti-OCT4 (Thermo Fisher Scientific, #A24867) at 1:100 and Mouse IgG3 Anti-SSEA4 (Thermo Fisher Scientific, #A24866) at 1:100 and incubated for 3 h at room temperature with gentle mixing. The cells were then washed two to three times in HBSS Fig. 6. RNA-Seq analysis denoting DEGs and enriched GO biological process terms in patient iPSCs and NSCs. (A) Scatter plot showing log 2 fold change on Y-axis and mean of normalized expression counts of on X-axis from differential expression genes between control (D149) versus patient (LSPH004) hiPSC RNA-Seq. Red and blue dots indicate statistical significant up-and downregulated genes, respectively (−1.5>Log 2 FC>1.5, FDR≤0.05). (B) Scatter plot showing log 2 fold change on Y-axis and mean of normalized expression counts on X-axis from differential expression genes between control (D149) versus patient (LSPH004) NSC RNA-Seq. Red and blue dots indicate statistical significant up-and downregulated genes, respectively (−1.5>Log 2 FC>1.5, FDR≤ 0.05). (C-E) Dot-plot showing GO biological process term enrichment, Y-axis shows enriched GO terms whereas X-axis shows fold enrichment of each term in upregulated genes of hiPSC LSPH004 (C), downregulated genes of hiPSC LSPH004 (D) and downregulated genes of NSC LSPH004 (E). Colors indicate thelog 10 FDR from Fisher's exact test, and dots size is proportional to the number of differentially expressed genes (DEG) in the given pathway. and incubated with appropriate secondary antibodies, Alexa Fluor™ 594 donkey anti-rabbit Antibody (Thermo Fisher Scientific, #A24870) and Alexa Fluor™ 488 goat anti-mouse IgG3 (Thermo Fisher Scientific, #A24877) for 1 h at room temperature. Before use, the secondary antibodies were diluted at 1:250 in 1% BSA. The cells were washed twice and then analyzed using a BD FACSVerse cytometer (Fig. S1A-C).

Western blotting
NSCs and 30 DIV neurons were harvested using Stempro Accutase and pelleted at 1000 g for 5 min then washed three times with ice-cold PBS. The pelleted cells were homogenized in 1X RIPA lysis buffer containing freshly added phosphatase and protease inhibitor cocktail (Roche). To remove cellular debris, crude RIPA lysates were centrifuged at 13,000 rpm for 20 min at 4°C. The supernatant was transferred to a new tube and quantified with a Pierce BCA protein assay (Thermo Fisher Scientific, #23225). Thereafter, the samples were heated at 95°C with Laemmli loading buffer for 5 min and 20 ug protein was loaded onto Bolt™ 4 to 12%, Bis-Tris SDS gel (Invitrogen, #NW04120BOX). The proteins were then transferred onto a nitrocellulose membrane and incubated overnight at 4°C with indicated antibodies. The blots were then washed three times with Tris Buffer Saline containing 0.1% Tween-20 (0.1% TBS-T) and incubated with 1:10,000 concentration of appropriate HRP-conjugated secondary antibodies (Jackson Laboratories, Inc.) for 45 min. After three washes with 0.1% TBS-T, blots were developed using Clarity Western ECL substrate (Bio-Rad) on a GE ImageQuant LAS 4000 system.

Calcium imaging
Calcium imaging was performed in 30 DIV D149 and LSPH004 neurons, according to our previously published protocol with minor modifications (Sharma et al., 2020). Briefly, neurons were washed with Tyrode's buffer solution (5 mM KCl, 129 mM NaCl, 2 mM CaCl 2 , 1 mM MgCl 2 , 30 mM glucose and 25 mM HEPES, pH 7.4) for 10 min and later incubated with 4 uM fluo-4/AM (1 mM, Molecular probes, #F14201) and 0.002% pluronic F-127 (Sigma-Aldrich, #P2443) in the Tyrode's buffer solution in dark for 30-45 min at room temperature. Following dye loading, the cells were washed again with the buffer thrice, each wash for 5 min. Finally, cells were incubated for an additional 20 min at room temperature to facilitate deesterification. Ca 2+ imaging was performed for 8 min with a time interval of 1 s using CellSens Dimension software (Olympus, build 16686) at 20X objective of wide-field fluorescence microscope Olympus IX-83. A 4-min baseline measurement was recorded to visualize calcium transients, followed by the addition of TTX to abolish calcium transients for another 4 min. Calcium traces were obtained using the CellSens Dimensions software by drawing a region of interest (ROI) manually around each neuronal soma. To plot the calcium traces, the raw fluorescence intensity values from each neuron were normalized to the first fluorescence intensity signal of the baseline recording. GraphPad Prism 5.0 was used to plot calcium traces.

Lipid extraction
Cells in culture were washed with DMEM/F12 to remove debris and then harvested using Stempro Accutase over a few minutes at 37°C. The Accutase was then neutralized, and the cell suspension was then centrifuged at 1000 g (for iPSCs) or 2500 g (for NSCs) to pellet down the cells. The supernatant was discarded, and the cell pellet was resuspended in 1 ml of 1X PBS for a wash and transferred to a 2 ml low-bind tube. The tubes were centrifuged at the previously indicated speeds to pellet the cells again. The supernatant was discarded, and the cell pellet was resuspended in 340 µl of 1X PBS and divided into two aliquots of 170 µl each for subsequent processing.
To each tube, 750 µl of ice-cold quench mixture, followed by 15 µl of a pre-mixed ISD mixture containing 25 ng of 37:4 PIP, 25 ng of 37:4 PIP 2 and 50 ng of 31:1 PE was added. The tubes were vortexed for 2 min at about 1500 rpm. Thereafter, 725 µl of CHCl 3 and 170 µl of 2.4 M HCl was added. The tubes were again vortexed for 2 min and kept at room temperature on the bench for 5 min. Two separate phases can be seen with a whitish precipitate at this stage.
All tubes were then spun at 1500 rpm in a benchtop centrifuge for 3 min to clearly separate the phases. In fresh 2 mL low-bind tubes, one for each sample, 708 µl of LPWS was added and kept aside. Once the phases were observed to be separate in the original tubes upon centrifugation, 900 µl of the lower organic phase was pipetted out by piercing through the upper phase and the protein precipitate at the interface and added to the tubes containing the LPWS. The tubes were then vortexed for 2 min and spun at 1500 rpm in a benchtop centrifuge for 3 min to separate the phases. The lower phase was pipetted out to extent possible taking care not to aspirate any of the upper phase and collected into a fresh tube.

Lipid derivatization
The following steps were performed inside a chemical hood, while wearing appropriate respirator mask. 50 µl of TMS-Diazomethane was added to each sample and incubated on a shaker at 600 rpm for 10 min at room temperature. At the end of 10 min, TMS-Diazomethane in each tube was quenched using 10 µl of Glacial Acetic acid. The tubes were inverted a few times to complete the quenching and carefully snapped open once to let the Nitrogen released during quenching to escape. At this point, the samples were moved out of the hood. 500 µl of the upper phase of PDWS was added to each sample. Tubes were then vortexed for 2 min and spun at 1500 rpm in a benchtop centrifuge for 3 min to separate the phases. 400 µl of the upper phase was discarded and another 500 µl of upper phase of PDWS was added to each tube and the vortex and spin steps were repeated as done earlier to separate the phases. Finally, the entire upper phase was discarded from each sample. 45 µl MeOH and 5 µl H2O was added to each tube, mixed and spun down. All the samples were then dried in a SpeedVac at 500 rpm for around 2 h till only about 10-20 µl of solvent was remaining. 90 µl of MeOH was added to reconstitute the sample to a final volume of about 100-110 µl and taken for injection and analysis.
On the mass spectrometer, we first employed Neutral Loss Scans during pilot standardization experiments on biological samples and searched for parent ions that would lose neutral fragments corresponding to 490 a.m.u and 382 a.m.u indicative of PIP 2 and PIP species, respectively, and likewise 155 a.m.u for PE species as described in (Clark et al., 2011). Thereafter, we quantified PIP, PIP 2 and PE species in biological samples using the selective multiple reaction monitoring (MRM) method in the positive ion mode (Table S4). For each sample, PE levels measured were used to normalize for total cellular phospholipid content on the assumption that PE levels are not likely to be different between control and LS cell lines based on previous studies. The Sciex MultiQuant software was used to quantify the area under the peaks. For each run, the area under curve for each species of PIP, PIP 2 and PE was normalized to PIP, PIP 2 and PE internal standard peak, respectively, to account for differences in loading and ion response. Thereafter, the sum of normalized areas for all the species of PIP or PIP 2 was then divided by the sum of normalized areas for all the species of PE in each of the biological samples to account for differences in total phospholipid extracted across samples.
Live cell imaging of PI(4,5)P 2 probe 5-10,000 NSC of LSPH004-PIP 2R and the control line D149-PIP 2R were seeded on 35 mm glass bottom 15 mm bore confocal dishes (Biostar LifeTech LLP, #BDD011035) in NEM and cultured for 48 h to reach a uniform confluency. Prior to imaging, NEM was aspirated out of the dishes and cells were incubated with Hoechst (Invitrogen, USA, #H3570) at 5 μM final concentration in NEM for 10 min to stain the nucleus. The cells were then washed with PBS and NEM was added prior to imaging. Confocal images were recorded by collecting a z-stack using an Olympus FV 3000 confocal microscope. This was performed in three independent biological replicates.

Image analysis
The quantification of mCherry::PH-PLCδ probe was performed manually by generating the maximum z-projections of middle few planes of cells from confocal slices. Thereafter, line profiles were drawn across clearly identifiable plasma membrane regions and their adjacent cytosolic regions and the ratios of mean intensities (PM/Cyt) for these line profiles were calculated for each cell and used to generate statistics (Sharma et al., 2019).

NGS sequencing
Total RNA was extracted from well characterized hiPSCs and NSCs of patient and control lines using TRIzol (Ambion, Life Technologies, #5596018) as per the manufacturer's protocol in two biological replicates. The RNA was quantified using Qubit4 dsDNA HS Assay Kit (Thermo Fisher Scientific, #Q32854) and run on a Bio-analyzer chip (Agilent High Sensitivity DNA Chip, #5067-4626) to assess integrity. Post NEB Next Poly(A) mRNA Magnetic Isolation Module (#E74906), 150 ng of total RNA (RIN values >9) was used per sample for the library preparation using the NEB Next Ultra II Directional RNA Library Prep Kit for Illumina (New England Biolabs, #E7760 L). The libraries were then sequenced on Illumina HiSeq 2500 sequencing platform using 2×100 bp sequencing format.

Bioinformatics analysis
Illumina sequenced paired-end reads were obtained from sequencing as mentioned above. The quality of processed reads (adapter removal and trimming) was evaluated using FastQC (version: 0.11.9). The RNA-seq reads were then mapped onto to the human reference genome (hg38) using HISAT2 (version: 2.1.0), the resulting BAM files containing the aligned reads were provided to HTSeq (version: 0.12.4) to obtain gene-level read count table using the reference annotation file (GTF format). We further utilized the iDEP (0.91) (Ge et al., 2018) to transform the read counts data using the regularized log (rlog) transformation method that is originally implemented in the DESeq2 [https://bioconductor.org/packages/release/bioc/html/DESeq2. html]. Differential expression, Enrichment analysis and GO (Biological Process) were then conducted using DESeq2 and default settings within the iDEP (0.91). Genes with log 2 FC greater than +1.5 and lesser than −1.5 were considered as up-and downregulated, respectively, while adhering to P-value<0.05. Such filtered genes were used for further downstream analysis, which included gathering enriched GO (Biological Process) terms, for which the filtered gene list was provided to ShinyGO v0.66 (Ge et al., 2020) (http:// bioinformatics.sdstate.edu/go/) with default filtering parameters. To remove redundant GO terms, we evaluated the gene lists for statistical significance by Fisher's test using PANTHER (Mi et al., 2021) (http://pantherdb.org/). With an FDR <0.05 cut-off, we obtained the enrichment list which was further submitted to REViGO (Supek et al., 2011) (http://revigo.irb.hr/), which refines the number of redundant functional terms based on semantic similarity between the ontology terms. Using a tight filter of dispensability<0.05, the GO (Biological Process) terms were obtained.

Statistics and software
Two-tailed unpaired student's t-test was used to compare datasets of two. One-way ANOVA with post hoc Tukey's multiple pairwise comparison was used whenever the experiment consisted of more than two biological groups. All statistical analyses were performed on Graph Pad Prism (version. 8). Schematics were created with biorender.com.

Ethics approval
This work was carried out under the ethics approval provided by the Institutional Ethics Committee, St. John's Medical College & Hospital, Bangalore (IEC Study Ref. No. 28 / 2017) and the Institutional Ethics Committee, National Centre for Biological Sciences, Bangalore (NCBS/ IEC-8/002). Institutional Stem Cell Committee Approval-National Centre for Biological Sciences (01/ICSCR/IX-06.0 l.2020-RP2).