Search

Age-related hearing loss (ARHL) is a major health concern among the elderly population. It is hoped that increasing our understanding of its underlying pathophysiological processes will lead to the development of novel therapies. Recent genome-wide association studies (GWASs) discovered a few dozen genetic variants in association with elevated risk for ARHL. Integrated analysis of GWAS results and transcriptomics data is a powerful approach for elucidating specific cell types that are involved in disease pathogenesis. Intriguingly, recent studies that applied such bioinformatics approaches to ARHL resulted in disagreeing findings as for the key cell types that are most strongly linked to the genetic pathogenesis of ARHL. These conflicting studies pointed either to cochlear sensory epithelial or to stria vascularis cells as the cell types most prominently involved in the genetic basis of ARHL. Seeking to resolve this discrepancy, we integrated the analysis of four ARHL GWAS datasets with four independent inner-ear single-cell RNA-sequencing datasets. Our analysis clearly points to the cochlear sensory epithelial cells as the key cells for the genetic predisposition to ARHL. We also explain the limitation of the bioinformatics analysis performed by previous studies that led to missing the enrichment for ARHL GWAS signal in sensory epithelial cells. Collectively, we show that cochlear epithelial cells, not stria vascularis cells, are the main inner-ear cells related to the genetic pathogenesis of ARHL., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2024 American Society of Human Genetics. All rights reserved.)

Motivation: Polygenic risk scores (PRSs) predict individuals' genetic risk of developing complex diseases. They summarize the effect of many variants discovered in genome-wide association studies (GWASs). However, to date, large GWASs exist primarily for the European population and the quality of PRS prediction declines when applied to other ethnicities. Genetic profiling of individuals in the discovery set (on which the GWAS was performed) and target set (on which the PRS is applied) is typically done by SNP arrays that genotype a fraction of common SNPs. Therefore, a key step in GWAS analysis and PRS calculation is imputing untyped SNPs using a panel of fully sequenced individuals. The imputation results depend on the ethnic composition of the imputation panel. Imputing genotypes with a panel of individuals of the same ethnicity as the genotyped individuals typically improves imputation accuracy. However, there has been no systematic investigation into the influence of the ethnic composition of imputation panels on the accuracy of PRS predictions when applied to ethnic groups that differ from the population used in the GWAS., Results: We estimated the effect of imputation of the target set on prediction accuracy of PRS when the discovery and the target sets come from different ethnic groups. We analyzed binary phenotypes on ethnically distinct sets from the UK Biobank and other resources. We generated ethnically homogenous panels, imputed the target sets, and generated PRSs. Then, we assessed the prediction accuracy obtained from each imputation panel. Our analysis indicates that using an imputation panel matched to the ethnicity of the target population yields only a marginal improvement and only under specific conditions., Availability and Implementation: The source code used for executing the analyses is this paper is available at https://github.com/Shamir-Lab/PRS-imputation-panels., (© The Author(s) 2024. Published by Oxford University Press.)

Introduction: Mutations in microRNA-96 (miR-96), a microRNA expressed within the hair cells (HCs) of the inner ear, result in progressive hearing loss in both mouse models and humans. In this study, we present the first HC-specific RNA-sequencing (RNA-seq) dataset from newborn Mir96 Dmdo heterozygous, homozygous mutant, and wildtype mice., Methods: Bulk RNA-seq was performed on HCs of newborn Mir96 Dmdo heterozygous, homozygous mutant, and wildtype mice. Differentially expressed gene analysis was conducted on Mir96 Dmdo homozygous mutant HCs compared to wildtype littermate controls, followed by GO term and protein-protein interaction analysis on these differentially expressed genes., Results: We identify 215 upregulated and 428 downregulated genes in the HCs of the Mir96 Dmdo homozygous mutant mice compared to their wildtype littermate controls. Many of the significantly downregulated genes in Mir96 Dmdo homozygous mutant HCs have established roles in HC development and/or known roles in deafness including Myo15a, Myo7a, Ush1c, Gfi1 , and Ptprq and have enrichment in gene ontology (GO) terms with biological functions such as sensory perception of sound. Interestingly, upregulated genes in Mir96 Dmdo homozygous mutants, including possible miR-96 direct targets, show higher wildtype expression in supporting cells compared to HCs., Conclusion: Our data further support a role for miR-96 in HC development, possibly as a repressor of supporting cell transcriptional programs in HCs. The HC-specific Mir96 Dmdo RNA-seq data set generated from this manuscript are now publicly available in a dedicated profile in the gene expression analysis resource (gEAR-https://umgear.org/p?l=miR96)., Competing Interests: Conflict of interest The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Skin pigmentation is paused after sun exposure; however, the mechanism behind this pausing is unknown. In this study, we found that the UVB-induced DNA repair system, led by the ataxia telangiectasia mutated (ATM) protein kinase, represses MITF transcriptional activity of pigmentation genes while placing MITF in DNA repair mode, thus directly inhibiting pigment production. Phosphoproteomics analysis revealed ATM to be the most significantly enriched pathway among all UVB-induced DNA repair systems. ATM inhibition in mouse or human skin, either genetically or chemically, induces pigmentation. Upon UVB exposure, MITF transcriptional activation is blocked owing to ATM-dependent phosphorylation of MITF on S414, which modifies MITF activity and interactome toward DNA repair, including binding to TRIM28 and RBBP4. Accordingly, MITF genome occupancy is enriched in sites of high DNA damage that are likely repaired. This suggests that ATM harnesses the pigmentation key activator for the necessary rapid, efficient DNA repair, thus optimizing the chances of the cell surviving. Data are available from ProteomeXchange with the identifier PXD041121., (Copyright © 2023 The Authors. Published by Elsevier Inc. All rights reserved.)

Background: Polygenic risk score (PRS), calculated based on genome-wide association studies (GWASs), can improve breast cancer (BC) risk assessment. To date, most BC GWASs have been performed in individuals of European (EUR) ancestry, and the generalisation of EUR-based PRS to other populations is a major challenge. In this study, we examined the performance of EUR-based BC PRS models in Ashkenazi Jewish (AJ) women., Methods: We generated PRSs based on data on EUR women from the Breast Cancer Association Consortium (BCAC). We tested the performance of the PRSs in a cohort of 2161 AJ women from Israel (1437 cases and 724 controls) from BCAC (BCAC cohort from Israel (BCAC-IL)). In addition, we tested the performance of these EUR-based BC PRSs, as well as the established 313-SNP EUR BC PRS, in an independent cohort of 181 AJ women from Hadassah Medical Center (HMC) in Israel., Results: In the BCAC-IL cohort, the highest OR per 1 SD was 1.56 (±0.09). The OR for AJ women at the top 10% of the PRS distribution compared with the middle quintile was 2.10 (±0.24). In the HMC cohort, the OR per 1 SD of the EUR-based PRS that performed best in the BCAC-IL cohort was 1.58±0.27. The OR per 1 SD of the commonly used 313-SNP BC PRS was 1.64 (±0.28)., Conclusions: Extant EUR GWAS data can be used for generating PRSs that identify AJ women with markedly elevated risk of BC and therefore hold promise for improving BC risk assessment in AJ women., Competing Interests: Competing interests: BCAC conflict of interest: MWB conducts research funded by Amgen, Novartis and Pfizer. PAF conducts research funded by Amgen, Novartis and Pfizer. He received Honoraria from Roche, Novartis and Pfizer. JV is one ofthe inventors of diagnosis and treatment of ERCC3-mutant cancer. AWK has a research funding for his institution from Myriad Genetics for an unrelated project (funding dates 2017–2019). UM has research collaborations with Mercy BioAnalytics, RNA Guardian, Dana Farber and iLOF (Intelligent Lab on Fiber). RAM is a consultant for Pharmavite., (© Author(s) (or their employer(s)) 2023. Re-use permitted under CC BY. Published by BMJ.)

Hypothesis: Metformin treatment will protect mice from noise-induced hearing loss (NIHL)., Background: We recently identified metformin as the top-ranking, Food and Drug Administration-approved drug to counter inner ear molecular changes induced by permanent threshold shift-inducing noise. This study is designed to functionally test metformin as a potential otoprotective drug against NIHL., Methods: Male and female B6CBAF1/J mice were obtained at 7 to 8 weeks of age. A cohort of the females underwent ovariectomy to simulate menopause and eliminate the effect of ovarian-derived estrogens. At 10 weeks of age, mice underwent a permanent threshold shift-inducing noise exposure (102.5 or 105 dB SPL, 8-16 kHz, 2 h). Auditory brainstem response (ABR) thresholds were obtained at baseline, 24 h after noise exposure, and 1 week after noise exposure. Mice were administered metformin (200 mg/kg/d) or a saline control in their drinking water after the baseline ABR and for the remainder of the study. After the 1-week ABR, mice were euthanized and cochlear tissue was analyzed., Results: Metformin treatment reduced the 1-week ABR threshold shift at 16 kHz ( p < 0.01; d = 1.20) and 24 kHz ( p < 0.01; d = 1.15) as well as outer hair cell loss in the 32-45.5 kHz range ( p < 0.0001; d = 2.37) in male mice. In contrast, metformin treatment did not prevent hearing loss or outer hair cell loss in the intact or ovariectomized female mice., Conclusions: Metformin exhibits sex-dependent efficacy as a therapeutic for NIHL. These data compel continued investigation into metformin's protective effects and demonstrate the importance of evaluating the therapeutic efficacy of drugs in subjects of both sexes., Competing Interests: The authors disclose no conflicts of interest., (Copyright © 2023 Written work prepared by employees of the Federal Government as part of their official duties is, under the U.S. Copyright Act, a “work of the United States Government” for which copyright protection under Title 17 of the United States Code is not available. As such, copyright does not extend to the contributions of employees of the Federal Government.)

During embryonic development, tissue-specific transcription factors and chromatin remodelers function together to ensure gradual, coordinated differentiation of multiple lineages. Here, we define this regulatory interplay in the developing retinal pigmented epithelium (RPE), a neuroectodermal lineage essential for the development, function and maintenance of the adjacent retina. We present a high-resolution spatial transcriptomic atlas of the developing mouse RPE and the adjacent ocular mesenchyme obtained by geographical position sequencing (Geo-seq) of a single developmental stage of the eye that encompasses young and more mature ocular progenitors. These transcriptomic data, available online, reveal the key transcription factors and their gene regulatory networks during RPE and ocular mesenchyme differentiation. Moreover, conditional inactivation followed by Geo-seq revealed that this differentiation program is dependent on the activity of SWI/SNF complexes, shown here to control the expression and activity of RPE transcription factors and, at the same time, inhibit neural progenitor and cell proliferation genes. The findings reveal the roles of the SWI/SNF complexes in controlling the intersection between RPE and neural cell fates and the coupling of cell-cycle exit and differentiation., Competing Interests: Competing interests The authors declare no competing or financial interests., (© 2023. Published by The Company of Biologists Ltd.)

The genetic disorder, ataxia-telangiectasia (A-T), is caused by loss of the homeostatic protein kinase, ATM, and combines genome instability, tissue degeneration, cancer predisposition, and premature aging. Primary fibroblasts from A-T patients exhibit premature senescence when grown at ambient oxygen concentration (21%). Here, we show that reducing oxygen concentration to a physiological level range (3%) dramatically extends the proliferative lifespan of human A-T skin fibroblasts. However, they still undergo senescence earlier than control cells grown under the same conditions and exhibit high genome instability. Comparative RNA-seq analysis of A-T and control fibroblasts cultured at 3% oxygen followed by cluster analysis of differentially expressed genes and functional enrichment analysis, revealed distinct transcriptional dynamics in A-T fibroblasts senescing in physiological oxygen concentration. While some transcriptional patterns were similar to those observed during replicative senescence of control cells, others were unique to the senescing A-T cells. We observed in them a robust activation of interferon-stimulated genes, with undetected expression the interferon genes themselves. This finding suggests an activation of a non-canonical cGAS-STING-mediated pathway, which presumably responds to cytosolic DNA emanating from extranuclear micronuclei detected in these cells. Senescing A-T fibroblasts also exhibited a marked, intriguely complex alteration in the expression of genes associated with extracellular matrix (ECM) remodeling. Notably, many of the induced ECM genes encode senescence-associated secretory phenotype (SASP) factors known for their paracrine pro-fibrotic effects. Our data provide a molecular dimension to the segmental premature aging observed in A-T patients and its associated symptoms, which develop as the patients advance in age., (© 2023 The Authors. Aging Cell published by Anatomical Society and John Wiley & Sons Ltd.)

Immunotherapy has revolutionized the treatment of advanced melanoma. Because the pathways mediating resistance to immunotherapy are largely unknown, we conducted transcriptome profiling of preimmunotherapy tumor biopsies from patients with melanoma that received PD-1 blockade or adoptive cell therapy with tumor-infiltrating lymphocytes. We identified two melanoma-intrinsic, mutually exclusive gene programs, which were controlled by IFNγ and MYC, and the association with immunotherapy outcome. MYC-overexpressing melanoma cells exhibited lower IFNγ responsiveness, which was linked with JAK2 downregulation. Luciferase activity assays, under the control of JAK2 promoter, demonstrated reduced activity in MYC-overexpressing cells, which was partly reversible upon mutagenesis of a MYC E-box binding site in the JAK2 promoter. Moreover, silencing of MYC or its cofactor MAX with siRNA increased JAK2 expression and IFNγ responsiveness of melanomas, while concomitantly enhancing the effector functions of T cells coincubated with MYC-overexpressing cells. Thus, we propose that MYC plays a pivotal role in immunotherapy resistance through downregulation of JAK2., (©2023 American Association for Cancer Research.)

Nutrient deprivation (starvation) induced by fasting and hypercaloric regimens are stress factors that can influence cell and tissue homeostasis in mammals. One of the key cellular responses to changes in nutrient availability is the cell survival pathway autophagy. While there has been much research into the protein networks regulating autophagy, less is known about the gene expression networks involved in this fundamental process. Here, we applied a network algorithm designed to analyse omics datasets, to identify sub-networks that are enriched for induced genes in response to starvation. This enabled us to identify two prominent active modules, one composed of key stress-induced transcription factors, including members of the Jun, Fos and ATF families, and the other comprising autophagosome sub-network genes, including ULK1. The results were validated in the brain, liver and muscle of fasting mice. Moreover, differential expression analysis of autophagy genes in the brain, liver and muscle of high-fat diet-exposed mice showed significant suppression of GABARAPL1 in the liver. Finally, our data provide a resource that may facilitate the future identification of regulators of autophagy.

Tissue-specific transcription factors (TFs) control the transcriptome through an association with noncoding regulatory regions (cistromes). Identifying the combination of TFs that dictate specific cell fate, their specific cistromes and examining their involvement in complex human traits remain a major challenge. Here, we focus on the retinal pigmented epithelium (RPE), an essential lineage for retinal development and function and the primary tissue affected in age-related macular degeneration (AMD), a leading cause of blindness. By combining mechanistic findings in stem-cell-derived human RPE, in vivo functional studies in mice and global transcriptomic and proteomic analyses, we revealed that the key developmental TFs LHX2 and OTX2 function together in transcriptional module containing LDB1 and SWI/SNF (BAF) to regulate the RPE transcriptome. Importantly, the intersection between the identified LHX2-OTX2 cistrome with published expression quantitative trait loci, ATAC-seq data from human RPE, and AMD genome-wide association study (GWAS) data, followed by functional validation using a reporter assay, revealed a causal genetic variant that affects AMD risk by altering TRPM1 expression in the RPE through modulation of LHX2 transcriptional activity on its promoter. Taken together, the reported cistrome of LHX2 and OTX2, the identified downstream genes and interacting co-factors reveal the RPE transcription module and uncover a causal regulatory risk single-nucleotide polymorphism (SNP) in the multifactorial common blinding disease AMD., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2023 Cohen-Gulkar et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

The peripheral and central auditory subsystems together form a complex sensory network that allows an organism to hear. The genetic programs of the two subsystems must therefore be tightly coordinated during development. Yet, their interactions and common expression pathways have never been systematically explored. MicroRNAs (miRNAs) are short non-coding RNAs that regulate gene expression and are essential for normal development of the auditory system. We performed mRNA and small-RNA sequencing of organs from both auditory subsystems at three critical developmental timepoints (E16, P0, P16) to obtain a comprehensive and unbiased insight of their expression profiles. Our analysis reveals common and organ-specific expression patterns for differentially regulated mRNAs and miRNAs, which could be clustered with a particular selection of functions such as inner ear development, Wnt signalling, K+ transport, and axon guidance, based on gene ontology. Bioinformatics detected enrichment of predicted targets of specific miRNAs in the clusters and predicted regulatory interactions by monitoring opposite trends of expression of miRNAs and their targets. This approach identified six miRNAs as strong regulatory candidates for both subsystems. Among them was miR-96, an established critical factor for proper development in both subsystems, demonstrating the strength of our approach. We suggest that other miRNAs identified by this analysis are also common effectors of proper hearing acquirement. This first combined comprehensive analysis of the developmental program of the peripheral and central auditory systems provides important data and bioinformatics insights into the shared genetic program of the two sensory subsystems and their regulation by miRNAs.

The use of alternative promoters, splicing, and cleavage and polyadenylation (APA) generates mRNA isoforms that expand the diversity and complexity of the transcriptome. Here, we uncovered thousands of previously undescribed 5' uncapped and polyadenylated transcripts (5' UPTs). We show that these transcripts resist exonucleases due to a highly structured RNA and N6-methyladenosine modification at their 5' termini. 5' UPTs appear downstream of APA sites within their host genes and are induced upon APA activation. Strong enrichment in polysomal RNA fractions indicates 5' UPT translational potential. Indeed, APA promotes downstream translation initiation, non-canonical protein output, and consistent changes to peptide presentation at the cell surface. Lastly, we demonstrate the biological importance of 5' UPTs using Bcl2, a prominent anti-apoptotic gene whose entire coding sequence is a 5' UPT generated from 5' UTR-embedded APA sites. Thus, APA is not only accountable for terminating transcripts, but also for generating downstream uncapped RNAs with translation potential and biological impact., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2022 The Authors. Published by Elsevier Inc. All rights reserved.)

Spatiotemporal gene expression patterns are governed to a large extent by the activity of enhancer elements, which engage in physical contacts with their target genes. Identification of enhancer-promoter (EP) links that are functional only in a specific subset of cell types is a key challenge in understanding gene regulation. We introduce CT-FOCS (cell type FOCS), a statistical inference method that uses linear mixed effect models to infer EP links that show marked activity only in a single or a small subset of cell types out of a large panel of probed cell types. Analyzing 808 samples from FANTOM5, covering 472 cell lines, primary cells and tissues, CT-FOCS inferred such EP links more accurately than recent state-of-the-art methods. Furthermore, we show that strictly cell type-specific EP links are very uncommon in the human genome., (© The Author(s) 2022. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Hearing loss is one of the top contributors to years lived with disability and is a risk factor for dementia. Molecular evidence on the cellular origins of hearing loss in humans is growing. Here, we performed a genome-wide association meta-analysis of clinically diagnosed and self-reported hearing impairment on 723,266 individuals and identified 48 significant loci, 10 of which are novel. A large proportion of associations comprised missense variants, half of which lie within known familial hearing loss loci. We used single-cell RNA-sequencing data from mouse cochlea and brain and mapped common-variant genomic results to spindle, root, and basal cells from the stria vascularis, a structure in the cochlea necessary for normal hearing. Our findings indicate the importance of the stria vascularis in the mechanism of hearing impairment, providing future paths for developing targets for therapeutic intervention in hearing loss., Competing Interests: Declaration of interests The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. C.R.C. is supported by the UK National Institute for Health Research (NIHR) Biomedical Research Center but the views expressed herein are his own and do not represent those of NIHR nor the UK Department of Health and Social Care. Y.Z. and V.M.L. are co-founders and shareholders of PersoMedix AB. In addition, V.M.L. is CEO and shareholder of HepaPredict AB and discloses support by the Robert Bosch Foundation, Merck KGaA, and Eli Lilly and Company. J.H.L. is co-founder and share-holder of Oscellaria AB. H.L. receives support from a consulting contract between Data Tecnica International and the National Institute on Aging (NIA), National Institutes of Health (NIH). M.A.N. received a competitive contract awarded to Data Tecnica International LLC by the National Institutes of Health to support open science research, and he also currently serves on the scientific advisory board for Clover Therapeutics and is an advisor to Neuron23 Inc as a data science fellow. P.F.S. is consultant and shareholder of Neumora., (Copyright © 2022 The Authors. Published by Elsevier Inc. All rights reserved.)

Motivation: Active module identification (AMI) is an essential step in many omics analyses. Such algorithms receive a gene network and a gene activity profile as input and report subnetworks that show significant over-representation of accrued activity signal ('active modules'). Such modules can point out key molecular processes in the analyzed biological conditions., Results: We recently introduced a novel AMI algorithm called DOMINO and demonstrated that it detects active modules that capture biological signals with markedly improved rate of empirical validation. Here, we provide an online server that executes DOMINO, making it more accessible and user-friendly. To help the interpretation of solutions, the server provides GO enrichment analysis, module visualizations and accessible output formats for customized downstream analysis. It also enables running DOMINO with various gene identifiers of different organisms., Availability and Implementation: The server is available at http://domino.cs.tau.ac.il. Its codebase is available at https://github.com/Shamir-Lab., (© The Author(s) 2022. Published by Oxford University Press.)

To date, genome-wide association studies have identified thousands of statistically-significant associations between genetic variants, and phenotypes related to a myriad of traits and diseases. A key goal for human-genetics research is to translate these associations into functional mechanisms. Popular gene-set analysis tools, like MAGMA, map variants to genes they might affect, and then integrate genome-wide association study data (that is, variant-level associations for a phenotype) to score genes for association with a phenotype. Gene scores are subsequently used in competitive gene-set analyses to identify biological processes that are enriched for phenotype association. By default, variants are mapped to genes in their proximity. However, many variants that affect phenotypes are thought to act at regulatory elements, which can be hundreds of kilobases away from their target genes. Thus, we explored the idea of augmenting a proximity-based mapping scheme with publicly-available datasets of regulatory interactions. We used MAGMA to analyze genome-wide association study data for ten different phenotypes, and evaluated the effects of augmentation by comparing numbers, and identities, of genes and gene sets detected as statistically significant between mappings. We detected several pitfalls and confounders of such "augmented analyses", and introduced ways to control for them. Using these controls, we demonstrated that augmentation with datasets of regulatory interactions only occasionally strengthened the enrichment for phenotype association amongst (biologically-relevant) gene sets for different phenotypes. Still, in such cases, genes and regulatory elements responsible for the improvement could be pinpointed. For instance, using brain regulatory-interactions for augmentation, we were able to implicate two acetylcholine receptor subunits involved in post-synaptic chemical transmission, namely CHRNB2 and CHRNE, in schizophrenia. Collectively, our study presents a critical approach for integrating regulatory interactions into gene-set analyses for genome-wide association study data, by introducing various controls to distinguish genuine results from spurious discoveries., Competing Interests: The authors have declared that no competing interests exist.

How the splicing machinery defines exons or introns as the spliced unit has remained a puzzle for 30 years. Here, we demonstrate that peripheral and central regions of the nucleus harbor genes with two distinct exon-intron GC content architectures that differ in the splicing outcome. Genes with low GC content exons, flanked by long introns with lower GC content, are localized in the periphery, and the exons are defined as the spliced unit. Alternative splicing of these genes results in exon skipping. In contrast, the nuclear center contains genes with a high GC content in the exons and short flanking introns. Most splicing of these genes occurs via intron definition, and aberrant splicing leads to intron retention. We demonstrate that the nuclear periphery and center generate different environments for the regulation of alternative splicing and that two sets of splicing factors form discrete regulatory subnetworks for the two gene architectures. Our study connects 3D genome organization and splicing, thus demonstrating that exon and intron definition modes of splicing occur in different nuclear regions., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2022 Elsevier Inc. All rights reserved.)

Cochlear supporting cells (SCs) are glia-like cells critical for hearing function. In the neonatal cochlea, the greater epithelial ridge (GER) is a mitotically quiescent and transient organ, which has been shown to nonmitotically regenerate SCs. Here, we ablated Lgr5+ SCs using Lgr5-DTR mice and found mitotic regeneration of SCs by GER cells in vivo. With lineage tracing, we show that the GER houses progenitor cells that robustly divide and migrate into the organ of Corti to replenish ablated SCs. Regenerated SCs display coordinated calcium transients, markers of the SC subtype inner phalangeal cells, and survive in the mature cochlea. Via RiboTag, RNA-sequencing, and gene clustering algorithms, we reveal 11 distinct gene clusters comprising markers of the quiescent and damaged GER, and damage-responsive genes driving cell migration and mitotic regeneration. Together, our study characterizes GER cells as mitotic progenitors with regenerative potential and unveils their quiescent and damaged translatomes., Competing Interests: The authors have declared that no competing interests exist.

Noise-induced hearing loss (NIHL) results from a complex interplay of damage to the sensory cells of the inner ear, dysfunction of its lateral wall, axonal retraction of type 1C spiral ganglion neurons, and activation of the immune response. We use RiboTag and single-cell RNA sequencing to survey the cell-type-specific molecular landscape of the mouse inner ear before and after noise trauma. We identify induction of the transcription factors STAT3 and IRF7 and immune-related genes across all cell-types. Yet, cell-type-specific transcriptomic changes dominate the response. The ATF3/ATF4 stress-response pathway is robustly induced in the type 1A noise-resilient neurons, potassium transport genes are downregulated in the lateral wall, mRNA metabolism genes are downregulated in outer hair cells, and deafness-associated genes are downregulated in most cell types. This transcriptomic resource is available via the Gene Expression Analysis Resource (gEAR; https://umgear.org/NIHL) and provides a blueprint for the rational development of drugs to prevent and treat NIHL., Competing Interests: Declaration of interests K.S.S., J.B.S., G.P., A.T.P., and J.B. are employees of Decibel Therapeutics. The data presented in this manuscript are registered for pending U.S. Provisional Patent Application, number: 63/151,249, title: “System and Methods for Cell Type-Specific Atlas for the Inner Ear Transcriptional Response to Acoustic Trauma,” and UMB docket number: RH-2021-073., (Copyright © 2021 The Authors. Published by Elsevier Inc. All rights reserved.)

The Wnt signaling pathways play fundamental roles during both development and adult homeostasis. Aberrant activation of the canonical Wnt signal transduction pathway is involved in many diseases including cancer, and is especially implicated in the development and progression of colorectal cancer. Although extensively studied, new genes, mechanisms and regulatory modulators involved in Wnt signaling activation or silencing are still being discovered. Here we applied a genome-scale CRISPR-Cas9 knockout (KO) screen based on Wnt signaling induced cell survival to reveal new inhibitors of the oncogenic, canonical Wnt pathway. We have identified several potential Wnt signaling inhibitors and have characterized the effects of the initiation factor DExH-box protein 29 (DHX29) on the Wnt cascade. We show that KO of DHX29 activates the Wnt pathway leading to upregulation of the Wnt target gene cyclin-D1, while overexpression of DHX29 inhibits the pathway. Together, our data indicate that DHX29 may function as a new canonical Wnt signaling tumor suppressor and demonstrates that this screening approach can be used as a strategy for rapid identification of novel Wnt signaling modulators., (© 2021. The Author(s).)

The auditory system is a complex sensory network with an orchestrated multilayer regulatory programme governing its development and maintenance. Accumulating evidence has implicated long non-coding RNAs (lncRNAs) as important regulators in numerous systems, as well as in pathological pathways. However, their function in the auditory system has yet to be explored. Using a set of specific criteria, we selected four lncRNAs expressed in the mouse cochlea, which are conserved in the human transcriptome and are relevant for inner ear function. Bioinformatic characterization demonstrated a lack of coding potential and an absence of evolutionary conservation that represent properties commonly shared by their class members. RNAscope®  analysis of the spatial and temporal expression profiles revealed specific localization to inner ear cells. Sub-cellular localization analysis presented a distinct pattern for each lncRNA and mouse tissue expression evaluation displayed a large variability in terms of level and location. Our findings establish the expression of specific lncRNAs in different cell types of the auditory system and present a potential pathway by which the lncRNA Gas5 acts in the inner ear. Studying lncRNAs and deciphering their functions may deepen our knowledge of inner ear physiology and morphology and may reveal the basis of as yet unresolved genetic hearing loss-related pathologies. Moreover, our experimental design may be employed as a reference for studying other inner ear-related lncRNAs, as well as lncRNAs expressed in other sensory systems.

Motivation: While promoter methylation is associated with reinforcing fundamental tissue identities, the methylation status of distant enhancers was shown by genome-wide association studies to be a powerful determinant of cell-state and cancer. With recent availability of long reads that report on the methylation status of enhancer-promoter pairs on the same molecule, we hypothesized that probing these pairs on the single-molecule level may serve the basis for detection of rare cancerous transformations in a given cell population. We explore various analysis approaches for deconvolving cell-type mixtures based on their genome-wide enhancer-promoter methylation profiles., Results: To evaluate our hypothesis we examine long-read optical methylome data for the GM12878 cell line and myoblast cell lines from two donors. We identified over 100 000 enhancer-promoter pairs that co-exist on at least 30 individual DNA molecules. We developed a detailed methodology for mixture deconvolution and applied it to estimate the proportional cell compositions in synthetic mixtures. Analysis of promoter methylation, as well as enhancer-promoter pairwise methylation, resulted in very accurate estimates. In addition, we show that pairwise methylation analysis can be generalized from deconvolving different cell types to subtle scenarios where one wishes to resolve different cell populations of the same cell-type., Availability and Implementation: The code used in this work to analyze single-molecule Bionano Genomics optical maps is available via the GitHub repository https://github.com/ebensteinLab/Single&#95;molecule&#95;methylation&#95;in&#95;EP., (© The Author(s) 2021. Published by Oxford University Press.)

Genome-wide association studies (GWAS) have identified numerous common genetic variants associated with complex human traits and diseases. However, the translation of GWAS discoveries into biological and clinical insights is highly challenging. In this study, we present a novel bioinformatics approach for enhancing the functional interpretation of GWAS signals, based on their integration with single-cell (sc)RNA-seq datasets that examine developmental processes. Our approach performs three tasks: (1) Identification of links between cell differentiation trajectories and traits; (2) Elucidation of biological processes and molecular pathways that underlie such trajectory-trait links; and (3) Prioritization of target genes that carry the links between trajectories, pathways and traits. We applied our method to a set of 11 traits of various pathologies, and 12 scRNA-seq datasets of diverse developmental processes, and it readily detected well-established biological connections, including those between the maturation of cortical inhibitory interneurons and schizophrenia, hepatocytes and cholesterol levels, and pancreatic beta-islet cells and type-2 diabetes. For each of these associations, our method pinpointed top candidate genes that are strongly associated with both the kinetics of the differentiation trajectory and the disease's genetic risk. By the identification of trajectory-disease links, molecular pathways that underlie them and prioritizing candidate risk genes, our method improves the understanding of the etiology of complex diseases, and thus holds promise for enhancing rational drug development that is aimed at targeting specific biological processes that mediate the genetic predisposition to diseases., Competing Interests: The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (© 2021 The Authors.)

The auditory system comprises the auditory periphery, engaged in sound transduction and the central auditory system, implicated in auditory information processing and perception. Recently, evidence mounted that the mammalian peripheral and central auditory systems share a number of genes critical for proper development and function. This bears implication for auditory rehabilitation and evolution of the auditory system. To analyze to which extent microRNAs (miRNAs) belong to genes shared between both systems, we characterize the expression pattern of 12 cochlea-abundant miRNAs in the central auditory system. Quantitative real-time PCR (qRT-PCR) demonstrated expression of all 12 genes in the cochlea, the auditory hindbrain and the non-auditory prefrontal cortex (PFC) at embryonic stage (E)16 and postnatal stages (P)0 and P30. Eleven of them showed differences in expression between tissues and nine between the developmental time points. Hierarchical cluster analysis revealed that the temporal expression pattern in the auditory hindbrain was more similar to the PFC than to the cochlea. Spatiotemporal expression analysis by RNA in situ hybridization demonstrated widespread expression throughout the cochlear nucleus complex (CNC) and the superior olivary complex (SOC) during postnatal development. Altogether, our data indicate that miRNAs represent a relevant class of genetic factors functioning across the auditory system. Given the importance of gene regulatory network (GRN) components for development, physiology and evolution, the 12 miRNAs provide promising entry points to gain insights into their molecular underpinnings in the auditory system.

Background: Frequent activation of the co-transcriptional factor YAP is observed in a large number of solid tumors. Activated YAP associates with enhancer loci via TEAD4-DNA-binding protein and stimulates cancer aggressiveness. Although thousands of YAP/TEAD4 binding-sites are annotated, their functional importance is unknown. Here, we aim at further identification of enhancer elements that are required for YAP functions., Results: We first apply genome-wide ChIP profiling of YAP to systematically identify enhancers that are bound by YAP/TEAD4. Next, we implement a genetic approach to uncover functions of YAP/TEAD4-associated enhancers, demonstrate its robustness, and use it to reveal a network of enhancers required for YAP-mediated proliferation. We focus on Enhancer TRAM2 , as its target gene TRAM2 shows the strongest expression-correlation with YAP activity in nearly all tumor types. Interestingly, TRAM2 phenocopies the YAP-induced cell proliferation, migration, and invasion phenotypes and correlates with poor patient survival. Mechanistically, we identify FSTL-1 as a major direct client of TRAM2 that is involved in these phenotypes. Thus, TRAM2 is a key novel mediator of YAP-induced oncogenic proliferation and cellular invasiveness., Conclusions: YAP is a transcription co-factor that binds to thousands of enhancer loci and stimulates tumor aggressiveness. Using unbiased functional approaches, we dissect YAP enhancer network and characterize TRAM2 as a novel mediator of cellular proliferation, migration, and invasion. Our findings elucidate how YAP induces cancer aggressiveness and may assist diagnosis of cancer metastasis.

Algorithms for active module identification (AMI) are central to analysis of omics data. Such algorithms receive a gene network and nodes' activity scores as input and report subnetworks that show significant over-representation of accrued activity signal ("active modules"), thus representing biological processes that presumably play key roles in the analyzed conditions. Here, we systematically evaluated six popular AMI methods on gene expression and GWAS data. We observed that GO terms enriched in modules detected on the real data were often also enriched on modules found on randomly permuted data. This indicated that AMI methods frequently report modules that are not specific to the biological context measured by the analyzed omics dataset. To tackle this bias, we designed a permutation-based method that empirically evaluates GO terms reported by AMI methods. We used the method to fashion five novel AMI performance criteria. Last, we developed DOMINO, a novel AMI algorithm, that outperformed the other six algorithms in extensive testing on GE and GWAS data. Software is available at https://github.com/Shamir-Lab., (© 2021 The Authors. Published under the terms of the CC BY 4.0 license.)

Despite the known importance of the transcription factors ATOH1, POU4F3 and GFI1 in hair cell development and regeneration, their downstream transcriptional cascades in the inner ear remain largely unknown. Here, we have used Gfi1 cre ;RiboTag mice to evaluate changes to the hair cell translatome in the absence of GFI1. We identify a systematic downregulation of hair cell differentiation genes, concomitant with robust upregulation of neuronal genes in the GFI1-deficient hair cells. This includes increased expression of neuronal-associated transcription factors (e.g. Pou4f1 ) as well as transcription factors that serve dual roles in hair cell and neuronal development (e.g. Neurod1 , Atoh1 and Insm1 ). We further show that the upregulated genes are consistent with the NEUROD1 regulon and are normally expressed in hair cells prior to GFI1 onset. Additionally, minimal overlap of differentially expressed genes in auditory and vestibular hair cells suggests that GFI1 serves different roles in these systems. From these data, we propose a dual mechanism for GFI1 in promoting hair cell development, consisting of repression of neuronal-associated genes as well as activation of hair cell-specific genes required for normal functional maturation., Competing Interests: Competing interestsThe authors declare no competing or financial interests., (© 2020. Published by The Company of Biologists Ltd.)

Alternative polyadenylation (APA) is emerging as a widespread regulatory layer since the majority of human protein-coding genes contain several polyadenylation (p(A)) sites in their 3'UTRs. By generating isoforms with different 3'UTR length, APA potentially affects mRNA stability, translation efficiency, nuclear export, and cellular localization. Polyadenylation sites are regulated by adjacent RNA cis-regulatory elements, the principals among them are the polyadenylation signal (PAS) AAUAAA and its main variant AUUAAA, typically located ~20-nt upstream of the p(A) site. Mutations in PAS and other auxiliary poly(A) cis-elements in the 3'UTR of several genes have been shown to cause human Mendelian diseases, and to date, only a few common SNPs that regulate APA were associated with complex diseases. Here, we systematically searched for SNPs that affect gene expression and human traits by modulation of 3'UTR APA. First, focusing on the variants most likely to exert the strongest effect, we identified 2,305 SNPs that interrupt the canonical PAS or its main variant. Implementing pA-QTL tests using GTEx RNA-seq data, we identified 330 PAS SNPs (called PAS pA-QTLs) that were significantly associated with the usage of their p(A) site. As expected, PAS-interrupting alleles were mostly linked with decreased cleavage at their p(A) site and the consequential 3'UTR lengthening. However, interestingly, in ~10% of the cases, the PAS-interrupting allele was associated with increased usage of an upstream p(A) site and 3'UTR shortening. As an indication of the functional effects of these PAS pA-QTLs on gene expression and complex human traits, we observed for few dozens of them marked colocalization with eQTL and/or GWAS signals. The PAS-interrupting alleles linked with 3'UTR lengthening were also strongly associated with decreased gene expression, indicating that shorter isoforms generated by APA are generally more stable than longer ones. Last, we carried out an extended, genome-wide analysis of 3'UTR variants and detected thousands of additional pA-QTLs having weaker effects compared to the PAS pA-QTLs., Competing Interests: The authors have declared that no competing interests exist.

Gene expression is regulated by the rates of synthesis and degradation of mRNAs, but how these processes are coordinated is poorly understood. Here, we show that reduced transcription dynamics of specific genes leads to enhanced m 6 A deposition, preferential activity of the CCR4-Not complex, shortened poly(A) tails, and reduced stability of the respective mRNAs. These effects are also exerted by internal ribosome entry site (IRES) elements, which we found to be transcriptional pause sites. However, when transcription dynamics, and subsequently poly(A) tails, are globally altered, cells buffer mRNA levels by adjusting the expression of mRNA degradation machinery. Stress-provoked global impediment of transcription elongation leads to a dramatic inhibition of the mRNA degradation machinery and massive mRNA stabilization. Accordingly, globally enhanced transcription, such as following B cell activation or glucose stimulation, has the opposite effects. This study uncovers two molecular pathways that maintain balanced gene expression in mammalian cells by linking transcription to mRNA stability., Competing Interests: Declaration of Interests The authors declare no competing interests., (Copyright © 2020 Elsevier Inc. All rights reserved.)

Data normalization is a critical step in RNA sequencing (RNA-seq) analysis, aiming to remove systematic effects from the data to ensure that technical biases have minimal impact on the results. Analyzing numerous RNA-seq datasets, we detected a prevalent sample-specific length effect that leads to a strong association between gene length and fold-change estimates between samples. This stochastic sample-specific effect is not corrected by common normalization methods, including reads per kilobase of transcript length per million reads (RPKM), Trimmed Mean of M values (TMM), relative log expression (RLE), and quantile and upper-quartile normalization. Importantly, we demonstrate that this bias causes recurrent false positive calls by gene-set enrichment analysis (GSEA) methods, thereby leading to frequent functional misinterpretation of the data. Gene sets characterized by markedly short genes (e.g., ribosomal protein genes) or long genes (e.g., extracellular matrix genes) are particularly prone to such false calls. This sample-specific length bias is effectively removed by the conditional quantile normalization (cqn) and EDASeq methods, which allow the integration of gene length as a sample-specific covariate. Consequently, using these normalization methods led to substantial reduction in GSEA false results while retaining true ones. In addition, we found that application of gene-set tests that take into account gene-gene correlations attenuates false positive rates caused by the length bias, but statistical power is reduced as well. Our results advocate the inspection and correction of sample-specific length biases as default steps in RNA-seq analysis pipelines and reiterate the need to account for intergene correlations when performing gene-set enrichment tests to lessen false interpretation of transcriptomic data., Competing Interests: The authors have declared that no competing interests exist.

Alternative polyadenylation (APA) is emerging as an important layer of gene regulation because the majority of mammalian protein-coding genes contain multiple polyadenylation (pA) sites in their 3' UTR. By alteration of 3' UTR length, APA can considerably affect post-transcriptional gene regulation. Yet, our understanding of APA remains rudimentary. Novel single-cell RNA sequencing (scRNA-seq) techniques allow molecular characterization of different cell types to an unprecedented degree. Notably, the most popular scRNA-seq protocols specifically sequence the 3' end of transcripts. Building on this property, we implemented a method for analysing patterns of APA regulation from such data. Analyzing multiple datasets from diverse tissues, we identified widespread modulation of APA in different cell types resulting in global 3' UTR shortening/lengthening and enhanced cleavage at intronic pA sites. Our results provide a proof-of-concept demonstration that the huge volume of scRNA-seq data that accumulates in the public domain offers a unique resource for the exploration of APA based on a very broad collection of cell types and biological conditions., (© The Author(s) 2019. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Estrogen receptor α (ERα) is an enhancer activating transcription factor, a key driver of breast cancer and a main target for cancer therapy. ERα-mediated gene regulation requires proper chromatin-conformation to facilitate interactions between ERα-bound enhancers and their target promoters. A major determinant of chromatin structure is the CCCTC-binding factor (CTCF), that dimerizes and together with cohesin stabilizes chromatin loops and forms the boundaries of topologically associated domains. However, whether CTCF-binding elements (CBEs) are essential for ERα-driven cell proliferation is unknown. To address this question in a global manner, we implemented a CRISPR-based functional genetic screen targeting CBEs located in the vicinity of ERα-bound enhancers. We identified four functional CBEs and demonstrated the role of one of them in inducing chromatin conformation changes in favor of activation of PREX1, a key ERα target gene in breast cancer. Indeed, high PREX1 expression is a bona-fide marker of ERα-dependency in cell lines, and is associated with good outcome after anti-hormonal treatment. Altogether, our data show that distinct CTCF-mediated chromatin structures are required for ERα- driven breast cancer cell proliferation., (© The Author(s) 2019. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Genome-wide analysis of cellular transcriptomes using RNA-seq or expression arrays is a major mainstay of current biological and biomedical research. EXPANDER (EXPression ANalyzer and DisplayER) is a comprehensive software package for analysis of expression data, with built-in support for 18 different organisms. It is designed as a "one-stop shop" platform for transcriptomic analysis, allowing for execution of all analysis steps starting with gene expression data matrix. Analyses offered include low-level preprocessing and normalization, differential expression analysis, clustering, bi-clustering, supervised grouping, high-level functional and pathway enrichment tests, and networks and motif analyses. A variety of options is offered for each step, using established algorithms, including many developed and published by our laboratory. EXPANDER has been continuously developed since 2003, having to date over 18,000 downloads and 540 citations. One of the innovations in the recent version is support for combined analysis of gene expression and ChIP-seq data to enhance the inference of transcriptional networks and their functional interpretation. EXPANDER implements cutting-edge algorithms and makes them accessible to users through user-friendly interface and intuitive visualizations. It is freely available to users at http://acgt.cs.tau.ac.il/expander/., (Copyright © 2019 The Authors. Published by Elsevier Ltd.. All rights reserved.)

Neuroadaptations in the brain reward system caused by excessive alcohol intake, lead to drinking escalation and alcohol use disorder phenotypes. Activity-dependent neuroprotective protein (ADNP) is crucial for brain development, and is implicated in neural plasticity in adulthood. Here, we discovered that alcohol exposure regulates Adnp expression in the mesolimbic system, and that Adnp keeps alcohol drinking in moderation, in a sex-dependent manner. Specifically, Sub-chronic alcohol treatment (2.5 g/kg/day for 7 days) increased Adnp mRNA levels in the dorsal hippocampus in both sexes, and in the nucleus accumbens of female mice, 24 h after the last alcohol injection. Long-term voluntary consumption of excessive alcohol quantities (~10-15 g/kg/24 h, 5 weeks) increased Adnp mRNA in the hippocampus of male mice immediately after an alcohol-drinking session, but the level returned to baseline after 24 h of withdrawal. In contrast, excessive alcohol consumption in females led to long-lasting reduction in hippocampal Adnp expression. We further tested the regulatory role of Adnp in alcohol consumption, using the Adnp haploinsufficient mouse model. We found that Adnp haploinsufficient female mice showed higher alcohol consumption and preference, compared to Adnp intact females, whereas no genotype difference was observed in males. Importantly, daily intranasal administration of the ADNP-snippet drug candidate NAP normalized alcohol consumption in Adnp haploinsufficient females. Finally, female Adnp haploinsufficient mice showed a sharp increase in alcohol intake after abstinence, suggesting that Adnp protects against relapse in females. The current data suggest that ADNP is a potential novel biomarker and negative regulator of alcohol-drinking behaviors.

Breast cancer is the most prevalent type of malignancy in women with ∼1.7 million new cases diagnosed annually, of which the majority express ERα (ESR1), a ligand-dependent transcription factor. Genome-wide chromatin binding maps suggest that ERα may control the expression of thousands of genes, posing a great challenge in identifying functional targets. Recently, we developed a CRISPR-Cas9 functional genetic screening approach to identify enhancers required for ERα-positive breast cancer cell proliferation. We validated several candidates, including CUTE, a putative ERα-responsive enhancer located in the first intron of CUEDC1 (CUE-domain containing protein). Here, we show that CUTE controls CUEDC1 expression, and that this interaction is essential for ERα-mediated cell proliferation. Moreover, ectopic expression of CUEDC1, but not a CUE-domain mutant, rescues the defects in CUTE activity. Finally, CUEDC1 expression correlates positively with ERα in breast cancer. Thus, CUEDC1 is a functional target gene of ERα and is required for breast cancer cell proliferation., (Copyright © 2018 The Netherlands cancer Institute. Published by Elsevier B.V. All rights reserved.)

The sensory cells that are responsible for hearing include the cochlear inner hair cells (IHCs) and outer hair cells (OHCs), with the OHCs being necessary for sound sensitivity and tuning 1 . Both cell types are thought to arise from common progenitors; however, our understanding of the factors that control the fate of IHCs and OHCs remains limited. Here we identify Ikzf2 (which encodes Helios) as an essential transcription factor in mice that is required for OHC functional maturation and hearing. Helios is expressed in postnatal mouse OHCs, and in the cello mouse model a point mutation in Ikzf2 causes early-onset sensorineural hearing loss. Ikzf2 cello/cello OHCs have greatly reduced prestin-dependent electromotile activity, a hallmark of OHC functional maturation, and show reduced levels of crucial OHC-expressed genes such as Slc26a5 (which encodes prestin) and Ocm. Moreover, we show that ectopic expression of Ikzf2 in IHCs: induces the expression of OHC-specific genes; reduces the expression of canonical IHC genes; and confers electromotility to IHCs, demonstrating that Ikzf2 can partially shift the IHC transcriptome towards an OHC-like identity.

Background: Our appreciation of the critical role of the genome's 3D organization in gene regulation is steadily increasing. Recent 3C-based deep sequencing techniques elucidated a hierarchy of structures that underlie the spatial organization of the genome in the nucleus. At the top of this hierarchical organization are chromosomal territories and the megabase-scale A/B compartments that correlate with transcriptional activity within cells. Below them are the relatively cell-type-invariant topologically associated domains (TADs), characterized by high frequency of physical contacts between loci within the same TAD, and are assumed to function as regulatory units. Within TADs, chromatin loops bring enhancers and target promoters to close spatial proximity. Yet, we still have only rudimentary understanding how differences in chromatin organization between different cell types affect cell-type-specific gene expression programs that are executed under basal and challenged conditions., Results: Here, we carried out a large-scale meta-analysis that integrated Hi-C data from thirteen different cell lines and dozens of ChIP-seq and RNA-seq datasets measured on these cells, either under basal conditions or after treatment. Pairwise comparisons between cell lines demonstrate a strong association between modulation of A/B compartmentalization, differential gene expression and transcription factor (TF) binding events. Furthermore, integrating the analysis of transcriptomes of different cell lines in response to various challenges, we show that A/B compartmentalization of cells under basal conditions significantly correlates not only with gene expression programs and TF binding profiles that are active under the basal condition but also with those induced in response to treatment. Yet, in pairwise comparisons between different cell lines, we find that a large portion of differential TF binding and gene induction events occur in genomic loci assigned to A compartment in both cell types, underscoring the role of additional critical factors in determining cell-type-specific transcriptional programs., Conclusions: Our results further indicate the role of dynamic genome organization in regulation of differential gene expression between different cell types and the impact of intra-TAD enhancer-promoter interactions that are established under basal conditions on both the basal and treatment-induced gene expression programs.

Background: Functional characterization of non-coding elements in the human genome is a major genomic challenge and the maturation of genome-editing technologies is revolutionizing our ability to achieve this task. Oncogene-induced senescence, a cellular state of irreversible proliferation arrest that is enforced following excessive oncogenic activity, is a major barrier against cancer transformation; therefore, bypassing oncogene-induced senescence is a critical step in tumorigenesis. Here, we aim at further identification of enhancer elements that are required for the establishment of this state., Results: We first apply genome-wide profiling of enhancer-RNAs (eRNAs) to systematically identify enhancers that are activated upon oncogenic stress. DNA motif analysis of these enhancers indicates AP-1 as a major regulator of the transcriptional program induced by oncogene-induced senescence. We thus constructed a CRISPR-Cas9 sgRNA library designed to target senescence-induced enhancers that are putatively regulated by AP-1 and used it in a functional screen. We identify a critical enhancer that we name Enh AP1-OIS1 and validate that mutating the AP-1 binding site within this element results in oncogene-induced senescence bypass. Furthermore, we identify FOXF1 as the gene regulated by this enhancer and demonstrate that FOXF1 mediates Enh AP1-OIS1 effect on the senescence phenotype., Conclusions: Our study elucidates a novel cascade mediated by AP-1 and FOXF1 that regulates oncogene-induced senescence and further demonstrates the power of CRISPR-based functional genomic screens in deciphering the function of non-coding regulatory elements in the genome.

The synchronized differentiation of neuronal and vascular tissues is crucial for normal organ development and function, although there is limited information about the mechanisms regulating the coordinated development of these tissues. The choroid vasculature of the eye serves as the main blood supply to the metabolically active photoreceptors, and develops together with the retinal pigmented epithelium (RPE). Here, we describe a novel regulatory relationship between the RPE transcription factors Pax6 and Sox9 that controls the timing of RPE differentiation and the adjacent choroid maturation. We used a novel machine learning algorithm tool to analyze high resolution imaging of the choroid in Pax6 and Sox9 conditional mutant mice. Additional unbiased transcriptomic analyses in mutant mice and RPE cells generated from human embryonic stem cells, as well as chromatin immunoprecipitation and high-throughput analyses, revealed secreted factors that are regulated by Pax6 and Sox9. These factors might be involved in choroid development and in the pathogenesis of the common blinding disease: age-related macular degeneration (AMD)., Competing Interests: Competing interestsThe authors declare no competing or financial interests., (© 2018. Published by The Company of Biologists Ltd.)

Oncogene-induced senescence (OIS), provoked in response to oncogenic activation, is considered an important tumor suppressor mechanism. Long non-coding RNAs (lncRNAs) are transcripts longer than 200 nt without a protein-coding capacity. Functional studies showed that deregulated lncRNA expression promote tumorigenesis and metastasis and that lncRNAs may exhibit tumor-suppressive and oncogenic function. Here, we first identified lncRNAs that were differentially expressed between senescent and non-senescent human fibroblast cells. Using RNA interference, we performed a loss-function screen targeting the differentially expressed lncRNAs, and identified lncRNA-OIS1 (lncRNA#32, AC008063.3 or ENSG00000233397) as a lncRNA required for OIS. Knockdown of lncRNA-OIS1 triggered bypass of senescence, higher proliferation rate, lower abundance of the cell-cycle inhibitor CDKN1A and high expression of cell-cycle-associated genes. Subcellular inspection of lncRNA-OIS1 indicated nuclear and cytosolic localization in both normal culture conditions as well as following oncogene induction. Interestingly, silencing lncRNA-OIS1 diminished the senescent-associated induction of a nearby gene (Dipeptidyl Peptidase 4, DPP4) with established role in tumor suppression. Intriguingly, similar to lncRNA-OIS1, silencing DPP4 caused senescence bypass, and ectopic expression of DPP4 in lncRNA-OIS1 knockdown cells restored the senescent phenotype. Thus, our data indicate that lncRNA-OIS1 links oncogenic induction and senescence with the activation of the tumor suppressor DPP4.