Publications

2026
Hamada, N. et al. Distinct pathophysiological mechanisms of CEP152 variants in microcephaly and brain abnormalities. EMBO Mol. Med. (2026).at <> Publisher's VersionAbstract
CEP152 is essential for centriole function and neurodevelopment, and pathogenic recessive variants in CEP152 cause primary microcephaly. We identified new compound heterozygous CEP152 variants, c.314 G > A,p.(W105*) and c.2689 A > T,p.(K897*), in a microcephalic patient and analyzed them alongside a homozygous variant c.95 A > C,p.(Q32P) associated with severe microcephaly with marked gyral simplification. In vitro assays revealed distinct effects: p.K897* prevented centrosomal localization, p.W105* led to protein degradation, and p.Q32P retained centrosomal targeting but disrupted binding to Polo-like kinase 4, a key centriole biogenesis kinase and CEP152 partner. In vivo, both Cep152W105*/K897* and Cep152Q32P/Q32P knock-in mice displayed microcephaly; notably, Cep152Q32P/Q32P mice also exhibited severe cortical defects during brain development. Cellular analyses revealed centrosome dysfunction, mitotic errors, and increased apoptosis, which were exacerbated in Cep152Q32P/Q32P brains. Morphological examination, including electron microscopy, further demonstrated structural abnormalities of the centrosomes and centrioles in Cep152Q32P/Q32P brains. Electrophysiological and gene expression analyses confirmed variant-specific neuronal impairments, which correlate with clinical severity. Collectively, these findings demonstrate that distinct CEP152 variants disrupt neurodevelopment through different mechanisms, thereby explaining the spectrum of microcephaly severity and associated phenotypes.
Dvir, E., Meshorer, E. & Shifman, S. Complex genotype–phenotype relationships in neurodevelopmental disorders. Trends in Genetics (2026). Publisher's VersionAbstract
With the advent of sequencing technologies in recent years, hundreds of high-confidence risk genes have been implicated in neurodevelopmental disorders (NDDs). However, individuals carrying pathogenic variants in the same gene frequently exhibit diverse clinical presentations, including varied symptoms and diagnoses. We propose that this heterogeneity arises from different interacting factors that modulate the phenotypic outcomes of pathogenic variants, including variant-level features, modifying variation across the genome, prenatal and early-life environmental exposures, and developmental noise. Resolving these factors requires integrative approaches that combine population-scale genetics and functional genomics with environmental monitoring and quantitative assessments of stochastic developmental variation. Advancing our understanding of these factors is critical to elucidating the etiology of NDDs and improving diagnostic and personalized therapeutic strategies.
Nishijo, T. et al. Kcnq2 R213 knock-in mice reveal variant- and region-specific mechanisms underlying self-limited familial neonatal-infantile epilepsy and early infantile developmental and epileptic encephalopathy. Acta Neuropathologica Communications 14, 76 (2026). Publisher's VersionAbstract
KCNQ2 variants cause a spectrum of neonatal epilepsies, ranging from self-limited familial neonatal-infantile epilepsy (SeLFNIE) to early infantile developmental and epileptic encephalopathy (EIDEE). Two distinct missense variants at the same residue, p.R213W and p.R213Q, are associated with SeLFNIE and EIDEE, respectively. This study aimed to elucidate the in vivo effects of these variants on brain development and neuronal excitability using two knock-in mouse models, Kcnq2(R213W/+) and Kcnq2(R213Q/+). We assessed survival, seizure susceptibility, histological and molecular phenotypes, and electrophysiological properties in cortical and hippocampal neurons, and conducted RNA sequencing analyses of cortical tissue to identify transcriptional alterations. Kcnq2(R213Q/+) mice exhibited tonic-clonic seizures, shortened lifespan, delayed cortical neuron migration, abnormal elongation of the axon initial segment in cortical neurons, and dentate gyrus-specific gliosis. In contrast, Kcnq2(R213W/+) mice showed a milder phenotype with transient seizures and largely preserved cortical function. RNA sequencing analyses revealed that p.R213Q selectively upregulated genes involved in endoplasmic reticulum stress and synaptic regulation, together with compensatory upregulation of potassium channel subunits. These findings demonstrate that the two Kcnq2 variants lead to distinct neurodevelopmental phenotypes, attributable not only to differential impairment of the M-current but also to aberrant cortical development and stress response pathways. In particular, p.R213Q induces sustained cortical hyperexcitability and axon initial segment abnormalities, whereas p.R213W results in the milder phenotype. The established knock-in models provide powerful tools for elucidating disease mechanisms of EIDEE and SeLFNIE, and developing targeted therapies for KCNQ2-related epilepsies.
Oweis, W. et al. ATF3-dependent formation of inclusion bodies in polyQ-expressing human iPSC-derived neurons confers cellular protection. Cell death and differentiation (2026). Publisher's VersionAbstract
Huntington's disease (HD) is an incurable, neurodegenerative disorder. While the causative mutation - CAG expansions within the coding region of the Huntingtin (HTT) gene - has been identified over 30 years ago, the pathological mechanisms underlying HD are still not clear. The abnormal CAG track encodes a polyglutamine (polyQ) expanded protein, which leads to HTT protein misfolding. These polyQ aggregates can form insoluble inclusion bodies (IBs); however, whether IBs are protective or detrimental remains debatable. Here we developed fluorescent iPSC-based human neuronal models for polyQ-related disorders. Comparing cell death in IB+ and IB- iPSC-derived neurons, growing side-by-side, we demonstrate that polyQ IBs have a significant protective effect. Remarkably, knocking out ATF3 prevented polyQ-IB formation and rendered the cells more vulnerable to induced stress. Taken together, our results reveal ATF3's role in polyQ IB formation in human NPCs, and demonstrate that polyQ IBs protect cells from stress-induced death.
Amelan, A. et al. CRISPR knockout screens reveal genes and pathways essential for neuronal differentiation and implicate PEDS1 in neurodevelopment. Nature Neuroscience 29, 592–603 (2026). Publisher's VersionAbstract
Neurodevelopmental disorders (NDDs) arise from disruptions in brain development, yet the underlying pathways remain incompletely understood. Here we demonstrate that genome-wide CRISPR knockout screens in mouse embryonic stem cells differentiating into neural lineages identify hundreds of essential genes, only a minority of which are currently implicated in NDDs. Dominant NDD genes were enriched for transcriptional regulators, whereas recessive NDD genes were predominantly involved in metabolic processes. Mouse models for eight genes (Eml1, Dusp26, Dynlrb2, Mta3, Peds1, Sgms1, Slitrk4 and Vamp3) revealed marked neuroanatomical abnormalities, including microcephaly in half of the cases. Focusing on PEDS1, a key enzyme in plasmalogen biosynthesis, we identified a bi-allelic variant in individuals with microcephaly, global developmental delay and congenital cataracts. In mice, Peds1 deficiency led to accelerated cell-cycle exit and impaired neuronal differentiation and migration. These pathways required for neural differentiation provide a genetic framework for discovering additional NDD genes.
2025
Rabany, O. et al. Dynamic rRNA Methylation Regulates Translation in the Hematopoietic System and is Essential for Stem Cell Fitness. Blood (2025).doi:10.1182/blood.2024028300 Publisher's VersionAbstract

Self-renewal and differentiation are at the basis of hematopoiesis. While it is known that tight regulation of translation is vital for hematopoietic stem cells' (HSCs) biology, the mechanisms underlying translation regulation across the hematopoietic system remain obscure. Here we reveal a novel mechanism of translation regulation in the hematopoietic hierarchy, which is mediated by ribosomal RNA (rRNA) methylation dynamics. Using ultra-low input ribosome-profiling, we characterized cell-type-specific translation capacity during erythroid differentiation. We found that translation efficiency changes progressively with differentiation and can distinguish between discrete cell populations as well as to define differentiation trajectories. To reveal the underlying mechanism, we performed comprehensive mapping of the most abundant rRNA modification - 2'-O-methyl (2'OMe). We found that, like translation efficiency, 2'OMe dynamics followed a distinct trajectory during erythroid differentiation.Genetic perturbation of individual 2'OMe sites demonstrated their distinct roles in modulating proliferation and differentiation. By combining CRISPR screening, molecular and functional analyses, we identified a specific methylation site, 28S-Gm4588, which is progressively lost during differentiation, as a key regulator of HSC self-renewal. We showed that low methylation at this site led to translational skewing, mediated mainly by codon frequency, which promoted differentiation. Functionally, HSCs with diminished 28S-Gm4588 methylation exhibited impaired self-renewal capacity ex-vivo, and loss of fitness in-vivo in bone marrow transplantations.Extending our findings beyond the hematopoietic system, we also found distinct dynamics of 2'OMe profiles during differentiation of non-hematopoietic stem cells. Our findings reveal rRNA methylation dynamics as a general mechanism for cell-type-specific translation, required for cell function and differentiation.

2024
Hamada, N., Nishijo, T., Iwamoto, I., Shifman, S. & Nagata, K.-I. Analyses of Conditional Knockout Mice for , a Gene Responsible for Neurodevelopmental Disorders in Excitatory and Inhibitory Neurons in the Brain. Cells 13, 540 (2024). Publisher's VersionAbstract

POGZ (Pogo transposable element derived with ZNF domain) is known to function as a regulator of gene expression. While variations in the gene have been associated with intellectual disabilities and developmental delays in humans, the exact pathophysiological mechanisms remain unclear. To shed light on this, we created two lines of conditional knockout mice for , one specific to excitatory neurons (Emx1-Pogz mice) and the other to inhibitory neurons (Gad2-Pogz mice) in the brain. Emx1-Pogz mice showed a decrease in body weight, similar to total knockout mice. Although the two lines did not display significant morphological abnormalities in the telencephalon, impaired POGZ function affected the electrophysiological properties of both excitatory and inhibitory neurons differently. These findings suggest that these mouse lines could be useful tools for clarifying the precise pathophysiological mechanisms of neurodevelopmental disorders associated with gene abnormalities.

Deshpande, D. et al. The evolution of computational research in a data-centric world. Cell 187, 4449-4457 (2024). Publisher's VersionAbstract

Computational data-centric research techniques play a prevalent and multi-disciplinary role in life science research. In the past, scientists in wet labs generated the data, and computational researchers focused on creating tools for the analysis of those data. Computational researchers are now becoming more independent and taking leadership roles within biomedical projects, leveraging the increased availability of public data. We are now able to generate vast amounts of data, and the challenge has shifted from data generation to data analysis. Here we discuss the pitfalls, challenges, and opportunities facing the field of data-centric research in biology. We discuss the evolving perception of computational data-driven research and its rise as an independent domain in biomedical research while also addressing the significant collaborative opportunities that arise from integrating computational research with experimental and translational biology. Additionally, we discuss the future of data-centric research and its applications across various areas of the biomedical field.

Moses, E. et al. The killifish germline regulates longevity and somatic repair in a sex-specific manner. Nature Aging 4, 791-813 (2024). Publisher's VersionAbstract

Classical evolutionary theories propose tradeoffs among reproduction, damage repair and lifespan. However, the specific role of the germline in shaping vertebrate aging remains largely unknown. In this study, we used the turquoise killifish (Nothobranchius furzeri) to genetically arrest germline development at discrete stages and examine how different modes of infertility impact life history. We first constructed a comprehensive single-cell gonadal atlas, providing cell-type-specific markers for downstream phenotypic analysis. We show here that germline depletion-but not arresting germline differentiation-enhances damage repair in female killifish. Conversely, germline-depleted males instead showed an extension in lifespan and rejuvenated metabolic functions. Through further transcriptomic analysis, we highlight enrichment of pro-longevity pathways and genes in germline-depleted male killifish and demonstrate functional conservation of how these factors may regulate longevity in germline-depleted Caenorhabditis elegans. Our results, therefore, demonstrate that different germline manipulation paradigms can yield pronounced sexually dimorphic phenotypes, implying alternative responses to classical evolutionary tradeoffs.

2023
Rosenski, J., Shifman, S. & Kaplan, T. Predicting gene knockout effects from expression data. BMC Med Genomics 16, 26 (2023). Publisher's VersionAbstract
BACKGROUND: The study of gene essentiality, which measures the importance of a gene for cell division and survival, is used for the identification of cancer drug targets and understanding of tissue-specific manifestation of genetic conditions. In this work, we analyze essentiality and gene expression data from over 900 cancer lines from the DepMap project to create predictive models of gene essentiality. METHODS: We developed machine learning algorithms to identify those genes whose essentiality levels are explained by the expression of a small set of "modifier genes". To identify these gene sets, we developed an ensemble of statistical tests capturing linear and non-linear dependencies. We trained several regression models predicting the essentiality of each target gene, and used an automated model selection procedure to identify the optimal model and hyperparameters. Overall, we examined linear models, gradient boosted trees, Gaussian process regression models, and deep learning networks. RESULTS: We identified nearly 3000 genes for which we accurately predict essentiality using gene expression data of a small set of modifier genes. We show that both in the number of genes we successfully make predictions for, as well as in the prediction accuracy, our model outperforms current state-of-the-art works. CONCLUSIONS: Our modeling framework avoids overfitting by identifying the small set of modifier genes, which are of clinical and genetic importance, and ignores the expression of noisy and irrelevant genes. Doing so improves the accuracy of essentiality prediction in various conditions and provides interpretable models. Overall, we present an accurate computational approach, as well as interpretable modeling of essentiality in a wide range of cellular conditions, thus contributing to a better understanding of the molecular mechanisms that govern tissue-specific effects of genetic disease and cancer.
Alkelai, A. et al. Genetic insights into childhood-onset schizophrenia: The yield of clinical exome sequencing. Schizophr Res 252, 138-145 (2023). Publisher's VersionAbstract
Childhood-onset schizophrenia (COS) is a rare form of schizophrenia with an onset prior to 13 years of age. Although genetic factors play a role in COS etiology, only a few causal variants have been reported to date. This study presents a diagnostic exome sequencing (ES) in 37 Israeli Jewish families with a proband diagnosed with COS. By implementing a trio/duo ES approach and applying a well-established diagnostic pipeline, we detected clinically significant variants in 7 probands (19 %). These single nucleotide variants and indels were mostly inherited. The implicated genes were ANKRD11, GRIA2, CHD2, CLCN3, CLTC, IGF1R and MICU1. In a secondary analysis that compared COS patients to 4721 healthy controls, we observed that patients had a significant enrichment of rare loss of function (LoF) variants in LoF intolerant genes associated with developmental diseases. Taken together, ES could be considered as a valuable tool in the genetic workup for COS patients.
2022
Shohat, S., Vol, E. & Shifman, S. Gene essentiality in cancer cell lines is modified by the sex chromosomes. Genome Res (2022).doi:10.1101/gr.276488.121 Publisher's VersionAbstract
Human sex differences arise from gonadal hormones and sex chromosomes. Studying the direct effects of sex chromosomes in humans is still challenging. Here we studied how the sex chromosomes can modulate gene expression and the outcome of mutations across the genome by exploiting the tendency of cancer cell lines to lose or gain sex chromosomes. We inferred the dosage of the sex chromosomes in 355 female and 408 male cancer cell lines and used it to dissect the contribution of the Y and X Chromosomes to sex-biased gene expression. Furthermore, based on genome-wide CRISPR screens, we identified genes whose essentiality is different between male and female cells depending on the sex chromosomes. The most significant genes were X-linked genes compensated by Y-linked paralogs. Our sex-based analysis identifies genes that, when mutated, can affect male and female cells differently and reinforces the role of the X and Y-Chromosomes in sex-specific cell function.
Herman, N., Kadener, S. & Shifman, S. The chromatin factor ROW cooperates with BEAF-32 in regulating long-range inducible genes. EMBO Rep e54720 (2022).doi:10.15252/embr.202254720 Publisher's VersionAbstract
Insulator proteins located at the boundaries of topological associated domains (TAD) are involved in higher-order chromatin organization and transcription regulation. However, it is still not clear how long-range contacts contribute to transcriptional regulation. Here, we show that relative-of-WOC (ROW) is essential for the long-range transcription regulation mediated by the boundary element-associated factor of 32kD (BEAF-32). We find that ROW physically interacts with heterochromatin proteins (HP1b and HP1c) and the insulator protein (BEAF-32). These proteins interact at TAD boundaries where ROW, through its AT-hook motifs, binds AT-rich sequences flanked by BEAF-32-binding sites and motifs. Knockdown of row downregulates genes that are long-range targets of BEAF-32 and bound indirectly by ROW (without binding motif). Analyses of high-throughput chromosome conformation capture (Hi-C) data reveal long-range interactions between promoters of housekeeping genes bound directly by ROW and promoters of developmental genes bound indirectly by ROW. Thus, our results show cooperation between BEAF-32 and the ROW complex, including HP1 proteins, to regulate the transcription of developmental and inducible genes through long-range interactions.
Dvir, E., Shohat, S., Flint, J. & Shifman, S. Identification of genetic mechanisms for tissue-specific genetic effects based on CRISPR screens. Genetics (2022).doi:10.1093/genetics/iyac134 Publisher's VersionAbstract
A major challenge in genetic studies of complex diseases is to determine how the action of risk genes is restricted to a tissue or cell type. Here we investigate tissue specificity of gene action using CRISPR screens from 786 cancer cell lines originating from 24 tissues. We find that the expression pattern of the gene across tissues explains only a minority of cases of tissue-specificity (9%), while gene amplification and the expression levels of paralogs account for 39.5% and 15.5%, respectively. Additionally, the transfer of small molecules to mutant cells explains tissue-specific gene action in blood. The tissue-specific genes we found are not specific just for human cancer cell lines: we found that the tissue-specific genes are intolerant to functional mutations in the human population and are associated with human diseases more than genes that are essential across all cell types. Our findings offer important insights into genetic mechanisms for tissue specificity of human diseases.
2021
Alkelai, A. et al. Expansion of the GRIA2 phenotypic representation: a novel de novo loss of function mutation in a case with childhood onset schizophrenia. J Hum Genet (2021).doi:10.1038/s10038-020-00846-1 Publisher's VersionAbstract
Childhood-onset schizophrenia (COS) is a rare form of schizophrenia with an onset before 13 years of age. There is rising evidence that genetic factors play a major role in COS etiology, yet, only a few single gene mutations have been discovered. Here we present a diagnostic whole-exome sequencing (WES) in an Israeli Jewish female with COS and additional neuropsychiatric conditions such as obsessive-compulsive disorder (OCD), anxiety, and aggressive behavior. Variant analysis revealed a de novo novel stop gained variant in GRIA2 gene (NM_000826.4: c.1522 G > T (p.Glu508Ter)). GRIA2 encodes for a subunit of the AMPA sensitive glutamate receptor (GluA2) that functions as ligand-gated ion channel in the central nervous system and plays an important role in excitatory synaptic transmission. GluA2 subunit mutations are known to cause variable neurodevelopmental phenotypes including intellectual disability, autism spectrum disorder, epilepsy, and OCD. Our findings support the potential diagnostic role of WES in COS, identify GRIA2 as possible cause to a broad psychiatric phenotype that includes COS as a major manifestation and expand the previously reported GRIA2 loss of function phenotypes.
Shohat, S., Amelan, A. & Shifman, S. Convergence and Divergence in the Genetics of Psychiatric Disorders From Pathways to Developmental Stages. Biol Psychiatry (2021).doi:10.1016/j.biopsych.2020.05.019 Publisher's VersionAbstract
In the past decade, the identification of susceptibility genes for psychiatric disorders has become routine, but understanding the biology underlying these discoveries has proven extremely difficult. The large number of potential risk genes and the genetic overlap between disorders are major obstacles for studying the etiology of these conditions. Systems biology approaches relying on gene ontologies, gene coexpression, and protein-protein interactions are used to identify convergence of the genes in relation to biological processes, cell types, brain areas, and developmental stages. Across psychiatric disorders, there is a clear enrichment for genes expressed in the brain and especially in the cortex, but a higher resolution is vastly dependent on sample size and statistical power. There is indication that susceptibility genes tend to be expressed in the brain during periods preceding the typical onset of the disorders. Thus, the role of genes in prenatal brain development is more pronounced for childhood-onset disorders, such as autism spectrum disorder and attention-deficit/hyperactivity disorder, but is much less so for bipolar disorder and depression. One of the most consistent findings across multiple disorders and classes of genetic variants is the role of genes intolerant to mutations in psychiatric disorders, yet this association is more pronounced for disorders with a clear neurodevelopmental component. Notwithstanding, a detailed understanding of the neurobiology of psychiatric disorders is still lacking. It is possible that it will only be revealed by studying the risk genes at the level of the development and function of neuronal networks and circuits.
2020
Winek, K. et al. Transfer RNA fragments replace microRNA regulators of the cholinergic poststroke immune blockade. Proc Natl Acad Sci U S A (2020).doi:10.1073/pnas.2013542117 Publisher's VersionAbstract
Stroke is a leading cause of death and disability. Recovery depends on a delicate balance between inflammatory responses and immune suppression, tipping the scale between brain protection and susceptibility to infection. Peripheral cholinergic blockade of immune reactions fine-tunes this immune response, but its molecular regulators are unknown. Here, we report a regulatory shift in small RNA types in patient blood sequenced 2 d after ischemic stroke, comprising massive decreases of microRNA levels and concomitant increases of transfer RNA fragments (tRFs) targeting cholinergic transcripts. Electrophoresis-based size-selection followed by qRT-PCR validated the top six up-regulated tRFs in a separate cohort of stroke patients, and independent datasets of small and long RNA sequencing pinpointed immune cell subsets pivotal to these responses, implicating CD14 monocytes in the cholinergic inflammatory reflex. In-depth small RNA targeting analyses revealed the most-perturbed pathways following stroke and implied a structural dichotomy between microRNA and tRF target sets. Furthermore, lipopolysaccharide stimulation of murine RAW 264.7 cells and human CD14 monocytes up-regulated the top six stroke-perturbed tRFs, and overexpression of stroke-inducible tRF-22-WE8SPOX52 using a single-stranded RNA mimic induced down-regulation of immune regulator Z-DNA binding protein 1. In summary, we identified a "changing of the guards" between small RNA types that may systemically affect homeostasis in poststroke immune responses, and pinpointed multiple affected pathways, which opens new venues for establishing therapeutics and biomarkers at the protein and RNA level.
Amir, N. et al. Value-complexity tradeoff explains mouse navigational learning. PLoS Comput Biol 16, e1008497 (2020). Publisher's VersionAbstract
We introduce a novel methodology for describing animal behavior as a tradeoff between value and complexity, using the Morris Water Maze navigation task as a concrete example. We develop a dynamical system model of the Water Maze navigation task, solve its optimal control under varying complexity constraints, and analyze the learning process in terms of the value and complexity of swimming trajectories. The value of a trajectory is related to its energetic cost and is correlated with swimming time. Complexity is a novel learning metric which measures how unlikely is a trajectory to be generated by a naive animal. Our model is analytically tractable, provides good fit to observed behavior and reveals that the learning process is characterized by early value optimization followed by complexity reduction. Furthermore, complexity sensitively characterizes behavioral differences between mouse strains.
Suliman-Lavie, R. et al. Pogz deficiency leads to transcription dysregulation and impaired cerebellar activity underlying autism-like behavior in mice. Nat Commun 11, 5836 (2020). Publisher's VersionAbstract
Several genes implicated in autism spectrum disorder (ASD) are chromatin regulators, including POGZ. The cellular and molecular mechanisms leading to ASD impaired social and cognitive behavior are unclear. Animal models are crucial for studying the effects of mutations on brain function and behavior as well as unveiling the underlying mechanisms. Here, we generate a brain specific conditional knockout mouse model deficient for Pogz, an ASD risk gene. We demonstrate that Pogz deficient mice show microcephaly, growth impairment, increased sociability, learning and motor deficits, mimicking several of the human symptoms. At the molecular level, luciferase reporter assay indicates that POGZ is a negative regulator of transcription. In accordance, in Pogz deficient mice we find a significant upregulation of gene expression, most notably in the cerebellum. Gene set enrichment analysis revealed that the transcriptional changes encompass genes and pathways disrupted in ASD, including neurogenesis and synaptic processes, underlying the observed behavioral phenotype in mice. Physiologically, Pogz deficiency is associated with a reduction in the firing frequency of simple and complex spikes and an increase in amplitude of the inhibitory synaptic input in cerebellar Purkinje cells. Our findings support a mechanism linking heterochromatin dysregulation to cerebellar circuit dysfunction and behavioral abnormalities in ASD.
Dinstein, I. et al. The National Autism Database of Israel: a Resource for Studying Autism Risk Factors, Biomarkers, Outcome Measures, and Treatment Efficacy. J Mol Neurosci 70, 1303-1312 (2020). Publisher's Version