The cells of the sensory epithelium, and not the stria vascularis, are the main cochlear cells related to the genetic pathogenesis of age-related hearing loss

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.

The predictive capacity of polygenic risk scores for disease risk is only moderately influenced by imputation panels tailored to the target population

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.

Evaluation of European-based polygenic risk score for breast cancer in Ashkenazi Jewish women in Israel

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.

Metformin Protects Against Noise-Induced Hearing Loss in Male Mice

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.

SWI/SNF complexes are required for retinal pigmented epithelium differentiation and for the inhibition of cell proliferation and neural differentiation programs

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.

Accelerated replicative senescence of ataxia-telangiectasia skin fibroblasts is retained at physiologic oxygen levels, with unique and common transcriptional patterns

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.

MYC Induces Immunotherapy and IFNγ Resistance Through Downregulation of JAK2

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.

Transcriptional profiling of the response to starvation and fattening reveals differential regulation of autophagy genes in mammals

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.

The LHX2-OTX2 transcriptional regulatory module controls retinal pigmented epithelium differentiation and underlies genetic risk for age-related macular degeneration

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.

Shared and organ-specific gene-expression programs during the development of the cochlea and the superior olivary complex

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.

Genome-wide association meta-analysis identifies 48 risk variants and highlights the role of the stria vascularis in hearing loss

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.

Incorporating regulatory interactions into gene-set analyses for GWAS data: A controlled analysis with the MAGMA tool

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.

Gene architecture directs splicing outcome in separate nuclear spatial regions

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.

The DOMINO web-server for active module identification analysis

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.

CT-FOCS: a novel method for inferring cell type-specific enhancer–promoter maps

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.

Lineage-tracing and translatomic analysis of damage-inducible mitotic cochlear progenitors identifies candidate genes regulating regeneration

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.

A cell-type-specific atlas of the inner ear transcriptional response to acoustic trauma

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.

Long reads capture simultaneous enhancer-promoter methylation status for cell-type deconvolution

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_molecule_methylation_in_EP.

Genetic mapping of developmental trajectories for complex traits and diseases

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.

Expression pattern of cochlear microRNAs in the mammalian auditory hindbrain

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.

A comprehensive enhancer screen identifies TRAM2 as a key and novel mediator of YAP oncogenesis

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 EnhancerTRAM2, 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.

DOMINO: a network-based active module identification algorithm with reduced rate of false calls

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.

Identification and characterization of key long non-coding RNAs in the mouse cochlea

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.

GFI1 functions to repress neuronal gene expression in the developing inner ear hair cells

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.

Systematic identification of functional SNPs interrupting 3'UTR polyadenylation signals

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.

mRNA molecules that encode for proteins end with a long stretch of adenosines, called poly(A) tail. The poly(A) tail contributes to the stability of the mRNA molecules, their translation to proteins and their import from the nucleus to the cytoplasm. The process of adding this tail to the mRNAs is called polyadenylation, and the termination site on the mRNAs at which the poly(A) tail is added is called the poly(A) site. In recent years it became evident that the vast majority of mRNAs of human genes contain several alternative poly(A) sites and their usage generates different mRNA isoforms that differ in their stability and translation efficiency. Therefore, alternative polyadenylation (APA) is emerging as a novel and important, yet underexplored, mechanism that regulate gene expression. The choice between alternative p(A) sites in an mRNA molecule is regulated by regulatory sequences located within a region in the mRNA called the 3’ untranslated region (3’UTR). A major challenge in present human genetics research is to understand how common genetic variants affect individuals’ health. In our study, we systematically identified dozens of genetic variants that affect the choice between alternative p(A) sites and demonstrated that by that, these variants influence the expression level of the target genes. Our results help to illuminate a novel mechanism by which genetic variants that are common in the population affect different traits including our risk for developing diseases.

Transcription Dynamics Regulate Poly(A) Tails and Expression of the RNA Degradation Machinery to Balance mRNA Levels

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⁶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.

Cell-type-specific analysis of alternative polyadenylation using single-cell transcriptomics data

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.

A CRISPR-Cas9 screen identifies essential CTCF anchor sites for estrogen receptor-driven breast cancer cell proliferation

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.

Recurrent functional misinterpretation of RNA-seq data caused by sample-specific gene length bias

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.

Analysis of numerous RNA-seq datasets reveals a recurrent sample-specific length bias that causes frequent false positive calls by gene-set enrichment analyses, leading to functional misinterpretation of the data. Its removal requires methods that allow the integration of gene length as sample-specific covariate.

LncRNA-OIS1 regulates DPP4 activation to modulate senescence induced by RAS

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.

Nuclear poly(A)-binding protein 1 is an ATM target and essential for DNA double-strand break repair

The DNA damage response (DDR) is an extensive signaling network that is robustly mobilized by DNA double-strand breaks (DSBs). The primary transducer of the DSB response is the protein kinase, ataxia-telangiectasia, mutated (ATM). Here, we establish nuclear poly(A)-binding protein 1 (PABPN1) as a novel target of ATM and a crucial player in the DSB response. PABPN1 usually functions in regulation of RNA processing and stability. We establish that PABPN1 is recruited to the DDR as a critical regulator of DSB repair. A portion of PABPN1 relocalizes to DSB sites and is phosphorylated on Ser95 in an ATM-dependent manner. PABPN1 depletion sensitizes cells to DSB-inducing agents and prolongs the DSB-induced G2/M cell-cycle arrest, and DSB repair is hampered by PABPN1 depletion or elimination of its phosphorylation site. PABPN1 is required for optimal DSB repair via both nonhomologous end-joining (NHEJ) and homologous recombination repair (HRR), and specifically is essential for efficient DNA-end resection, an initial, key step in HRR. Using mass spectrometry analysis, we capture DNA damage-induced interactions of phospho-PABPN1, including well-established DDR players as well as other RNA metabolizing proteins. Our results uncover a novel ATM-dependent axis in the rapidly growing interface between RNA metabolism and the DDR.

Activity-dependent neuroprotective protein (ADNP) is an alcohol-responsive gene and negative regulator of alcohol consumption in female mice

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.

Genomic meta-analysis of the interplay between 3D chromatin organization and gene expression programs under basal and stress conditions

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.

Electronic supplementary material

The online version of this article (10.1186/s13072-018-0220-2) contains supplementary material, which is available to authorized users.

Functional CRISPR screen identifies AP1-associated enhancer regulating FOXF1 to modulate oncogene-induced senescence

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 EnhAP1-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 EnhAP1-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.

Pax6 regulation of Sox9 in the mouse retinal pigmented epithelium controls its timely differentiation and choroid vasculature development

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).

FOCS: a novel method for analyzing enhancer and gene activity patterns infers an extensive enhancer-promoter map

Recent sequencing technologies enable joint quantification of promoters and their enhancer regions, allowing inference of enhancer–promoter links. We show that current enhancer–promoter inference methods produce a high rate of false positive links. We introduce FOCS, a new inference method, and by benchmarking against ChIA-PET, HiChIP, and eQTL data show that it results in lower false discovery rates and at the same time higher inference power. By applying FOCS to 2630 samples taken from ENCODE, Roadmap Epigenomics, FANTOM5, and a new compendium of GRO-seq samples, we provide extensive enhancer–promotor maps (http://acgt.cs.tau.ac.il/focs). We illustrate the usability of our maps for deriving biological hypotheses.

Proteomic analysis of polyribosomes identifies splicing factors as potential regulators of translation during mitosis

Precise regulation of mRNA translation is critical for proper cell division, but little is known about the factors that mediate it. To identify mRNA-binding proteins that regulate translation during mitosis, we analyzed the composition of polysomes from interphase and mitotic cells using unbiased quantitative mass-spectrometry (LC–MS/MS). We found that mitotic polysomes are enriched with a subset of proteins involved in RNA processing, including alternative splicing and RNA export. To demonstrate that these may indeed be regulators of translation, we focused on heterogeneous nuclear ribonucleoprotein C (hnRNP C) as a test case and confirmed that it is recruited to elongating ribosomes during mitosis. Then, using a combination of pulsed SILAC, metabolic labeling and ribosome profiling, we showed that knockdown of hnRNP C affects both global and transcript-specific translation rates and found that hnRNP C is specifically important for translation of mRNAs that encode ribosomal proteins and translation factors. Taken together, our results demonstrate how proteomic analysis of polysomes can provide insight into translation regulation under various cellular conditions of interest and suggest that hnRNP C facilitates production of translation machinery components during mitosis to provide daughter cells with the ability to efficiently synthesize proteins as they enter G1 phase.

Transcription Impacts the Efficiency of mRNA Translation via Co-transcriptional N6-adenosine Methylation

Transcription and translation are two main pillars of gene expression. Due to the different timings, spots of action, and mechanisms of regulation, these processes are mainly regarded as distinct and generally uncoupled, despite serving a common purpose. Here, we sought for a possible connection between transcription and translation. Employing an unbiased screen of multiple human promoters, we identified a positive effect of TATA box on translation and a general coupling between mRNA expression and translational efficiency. Using a CRISPR-Cas9-mediated approach, genome-wide analyses, and in vitro experiments, we show that the rate of transcription regulates the efficiency of translation. Furthermore, we demonstrate that m⁶A modification of mRNAs is co-transcriptional and depends upon the dynamics of the transcribing RNAPII. Suboptimal transcription rates lead to elevated m⁶A content, which may result in reduced translation. This study uncovers a general and widespread link between transcription and translation that is governed by epigenetic modification of mRNAs.

•

Transcription rates of mRNAs positively correlate with rates of their translation

•

Dynamics of RNA polymerase II impact the deposition of m⁶A on mRNAs

•

Suboptimal transcription enhances m⁶A modification of mRNAs

•

Excessive m⁶A modification is detrimental for the translation process

BRCA1185delAG tumors may acquire therapy resistance through expression of RING-less BRCA1

Heterozygous germline mutations in breast cancer 1 (BRCA1) strongly predispose women to breast cancer. BRCA1 plays an important role in DNA double-strand break (DSB) repair via homologous recombination (HR), which is important for tumor suppression. Although BRCA1-deficient cells are highly sensitive to treatment with DSB-inducing agents through their HR deficiency (HRD), BRCA1-associated tumors display heterogeneous responses to platinum drugs and poly(ADP-ribose) polymerase (PARP) inhibitors in clinical trials. It is unclear whether all pathogenic BRCA1 mutations have similar effects on the response to therapy. Here, we have investigated mammary tumorigenesis and therapy sensitivity in mice carrying the Brca1185stop and Brca15382stop alleles, which respectively mimic the 2 most common BRCA1 founder mutations, BRCA1185delAG and BRCA15382insC. Both the Brca1185stop and Brca15382stop mutations predisposed animals to mammary tumors, but Brca1185stop tumors responded markedly worse to HRD-targeted therapy than did Brca15382stop tumors. Mice expressing Brca1185stop mutations also developed therapy resistance more rapidly than did mice expressing Brca15382stop. We determined that both murine Brca1185stop tumors and human BRCA1185delAG breast cancer cells expressed a really interesting new gene domain-less (RING-less) BRCA1 protein that mediated resistance to HRD-targeted therapies. Together, these results suggest that expression of RING-less BRCA1 may serve as a marker to predict poor response to DSB-inducing therapy in human cancer patients.

3'UTR Shortening Potentiates MicroRNA-Based Repression of Pro-differentiation Genes in Proliferating Human Cells

Most mammalian genes often feature alternative polyadenylation (APA) sites and hence diverse 3’UTR lengths. Proliferating cells were reported to favor APA sites that result in shorter 3’UTRs. One consequence of such shortening is escape of mRNAs from targeting by microRNAs (miRNAs) whose binding sites are eliminated. Such a mechanism might provide proliferation-related genes with an expression gain during normal or cancerous proliferation. Notably, miRNA sites tend to be more active when located near both ends of the 3’UTR compared to those located more centrally. Accordingly, miRNA sites located near the center of the full 3’UTR might become more active upon 3'UTR shortening. To address this conjecture we performed 3' sequencing to determine the 3' ends of all human UTRs in several cell lines. Remarkably, we found that conserved miRNA binding sites are preferentially enriched immediately upstream to APA sites, and this enrichment is more prominent in pro-differentiation/anti-proliferative genes. Binding sites of the miR17-92 cluster, upregulated in rapidly proliferating cells, are particularly enriched just upstream to APA sites, presumably conferring stronger inhibitory activity upon shortening. Thus 3’UTR shortening appears not only to enable escape from inhibition of growth promoting genes but also to potentiate repression of anti-proliferative genes.

MicroRNAs (miRNAs) are regulators of gene expression. Typically they recognize a binding site in genes' sequences and exert a repressive effect. This scheme prescribes a regulatory network that determines which gene is regulated by which miRNA. Yet this is a static sequence-based scheme that might not support dynamic changes in network wiring. Can genes become subject to, or be released, from the regulation of a miRNA in a manner that depends on the physiological state of cells? Here we describe such a dynamic mechanism. It is established that miRNA regulation is often more effective when their binding sites reside near the end of the target mRNA or right after the coding sequence STOP codon. Thus, a site distant from the mRNA’s end might be latent as it will not bind efficiently its corresponding miRNA. Yet, in particular physiological states, e.g. cancer or rapidly proliferating cells, mRNA ends tend to become shortened. So far it was suggested that such shortening may serve to release proliferation-favoring genes from miRNA repression by eliminating their binding sites from the mRNA. We propose a mirror image, complementary mechanism that acts upon genes that need to be repressed during proliferation. Specifically, we propose that mRNA shortening can dynamically activate latent repressive miRNA binding sites by bringing the mRNA end close to them. We mapped the ends of all mRNAs in proliferating cells and found that cancer-enriched ends are strikingly positioned closely downstream to a high density of potentially latent binding sites, which are retained in the short mRNA but are now close to the new ends. This may enable such potentially latent miRNA sites to become dynamically activated upon proliferation. Remarkably, this mechanism targets preferentially pro-differentiation and antiproliferative genes, which are often repressed in cancer.

Genome-Wide Polyadenylation Maps Reveal Dynamic mRNA 3′-End Formation in the Failing Human Heart

Rationale:

Alternative cleavage and polyadenylation (APA) of mRNA represents a layer of gene regulation that to date has remained unexplored in the heart. This phenomenon may be relevant, as the positioning of the poly(A) tail in mRNAs influences the length of the 3′-untranslated region (UTR), a critical determinant of gene expression.

Objective:

To investigate whether the 3′UTR length is regulated by APA in the human heart and whether this changes in the failing heart.

Methods and Results:

We used 3′end RNA sequencing (e3′-Seq) to directly measure global patterns of APA in healthy and failing human heart specimens. By monitoring polyadenylation profiles in these hearts, we identified disease-specific APA signatures in numerous genes. Interestingly, many of the genes with shortened 3′UTRs in heart failure were enriched for functional groups such as RNA binding, whereas genes with longer 3′UTRs were enriched for cytoskeletal organization and actin binding. RNA sequencing in a larger series of human hearts revealed that these APA candidates are often differentially expressed in failing hearts, with an inverse correlation between 3′UTR length and the level of gene expression. Protein levels of the APA regulator, poly(A)-binding protein nuclear-1 were substantially downregulated in failing hearts.

Conclusions:

We provide genome-wide, high-resolution polyadenylation maps of the human heart and show that the 3′end formation of mRNA is dynamic in heart failure, suggesting that APA-mediated 3′UTR length modulation represents an additional layer of gene regulation in failing hearts.

Functional genetic screens for enhancer elements in the human genome using CRISPR-Cas9

Systematic identification of noncoding regulatory elements has, to date, mainly relied on large-scale reporter assays that do not reproduce endogenous conditions. We present two distinct CRISPR-Cas9 genetic screens to identify and characterize functional enhancers in their native context. Our strategy is to target Cas9 to transcription factor binding sites in enhancer regions. We identified several functional enhancer elements and characterized the role of two of them in mediating p53 (TP53) and ERα (ESR1) gene regulation. Moreover, we show that a genomic CRISPR-Cas9 tiling screen can precisely map functional domains within enhancer elements. Our approach expands the utility of CRISPR-Cas9 to elucidate the functions of the noncoding genome.

Myc coordinates transcription and translation to enhance transformation and suppress invasiveness

c‐Myc is one of the major human proto‐oncogenes and is often associated with tumor aggression and poor clinical outcome. Paradoxically, Myc was also reported as a suppressor of cell motility, invasiveness, and metastasis. Among the direct targets of Myc are many components of the protein synthesis machinery whose induction results in an overall increase in protein synthesis that empowers tumor cell growth. At present, it is largely unknown whether beyond the global enhancement of protein synthesis, Myc activation results in translation modulation of specific genes. Here, we measured Myc‐induced global changes in gene expression at the transcription, translation, and protein levels and uncovered extensive transcript‐specific regulation of protein translation. Particularly, we detected a broad coordination between regulation of transcription and translation upon modulation of Myc activity and showed the connection of these responses to mTOR signaling to enhance oncogenic transformation and to the TGFβ pathway to modulate cell migration and invasiveness. Our results elucidate novel facets of Myc‐induced cellular responses and provide a more comprehensive view of the consequences of its activation in cancer cells.

RFX transcription factors are essential for hearing in mice

Sensorineural hearing loss is a common and currently irreversible disorder, because mammalian hair cells (HCs) do not regenerate and current stem cell and gene delivery protocols result only in immature HC-like cells. Importantly, although the transcriptional regulators of embryonic HC development have been described, little is known about the postnatal regulators of maturating HCs. Here we apply a cell type-specific functional genomic analysis to the transcriptomes of auditory and vestibular sensory epithelia from early postnatal mice. We identify RFX transcription factors as essential and evolutionarily conserved regulators of the HC-specific transcriptomes, and detect Rfx1,2,3,5 and 7 in the developing HCs. To understand the role of RFX in hearing, we generate Rfx1/3 conditional knockout mice. We show that these mice are deaf secondary to rapid loss of initially well-formed outer HCs. These data identify an essential role for RFX in hearing and survival of the terminally differentiating outer HCs.

Parallel profiling of the transcriptome, cistrome, and epigenome in the cellular response to ionizing radiation

The DNA damage response (DDR) is a vast signaling network that is robustly activated by DNA double-strand breaks, the critical lesion induced by ionizing radiation (IR). Although much of this response operates at the protein level, a critical component of the network sustains many DDR branches by modulating the cellular transcriptome. Using deep sequencing, we delineated three layers in the transcriptional response to IR in human breast cancer cells: changes in the expression of genes encoding proteins or long noncoding RNAs, alterations in genomic binding by key transcription factors, and dynamics of epigenetic markers of active promoters and enhancers. We identified protein-coding and previously unidentified noncoding genes that were responsive to IR, and demonstrated that IR-induced transcriptional dynamics was mediated largely by the transcription factors p53 and nuclear factor κB (NF-κB) and was primarily dependent on the kinase ataxia-telangiectasia mutated (ATM). The resultant data set provides a rich resource for understanding a basic, underlying component of a critical cellular stress response.

Alternative cleavage and polyadenylation: extent, regulation and function

The 3' end of most protein-coding genes and long non-coding RNAs is cleaved and polyadenylated. Recent discoveries have revealed that a large proportion of these genes contains more than one polyadenylation site. Therefore, alternative polyadenylation (APA) is a widespread phenomenon, generating mRNAs with alternative 3' ends. APA contributes to the complexity of the transcriptome by generating isoforms that differ either in their coding sequence or in their 3' untranslated regions (UTRs), thereby potentially regulating the function, stability, localization and translation efficiency of target RNAs. Here, we review our current understanding of the polyadenylation process and the latest progress in the identification of APA events, mechanisms that regulate poly(A) site selection, and biological processes and diseases resulting from APA.

p53 induces transcriptional and translational programs to suppress cell proliferation and growth

BACKGROUND

Cell growth and proliferation are tightly connected to ensure that appropriately sized daughter cells are generated following mitosis. Energy stress blocks cell growth and proliferation, a critical response for survival under extreme conditions. Excessive oncogenic stress leads to p53 activation and the induction of senescence, an irreversible state of cell-cycle arrest and a critical component in the suppression of tumorigenesis. Nutrient-sensing and mitogenic cues converge on a major signaling node, which regulates the activity of the mTOR kinase. Although transcriptional responses to energy and oncogenic stresses have been examined by many gene-expression experiments, a global exploration of the modulation of mRNA translation in response to these conditions is lacking.

RESULTS

We combine RNA sequencing and ribosomal profiling analyses to systematically delineate modes of transcriptional and translational regulation induced in response to conditions of limited energy, oncogenic stress and cellular transformation. We detect a key role for mTOR and p53 in these distinct physiological states, and provide the first genome-wide demonstration that p53 activation results in mTOR inhibition and a consequent global repression of protein translation. We confirm the role of the direct p53 target genes Sestrin1 and Sestrin2 in this response, as part of the broad modulation of gene expression induced by p53 activation.

CONCLUSIONS

We delineate a bimodal tumor-suppressive regulatory program activated by p53, in which cell-cycle arrest is imposed mainly at the transcriptional level, whereas cell growth inhibition is enforced by global repression of the translation machinery.

Alternative cleavage and polyadenylation during colorectal cancer development

PURPOSE

Alternative cleavage and polyadenylation (APA) of mRNAs is a phenomenon that alters 3'-untranslated region length leading to altered posttranscriptional regulation of gene expression. Changing APA patterns have been shown to result in misregulation of genes involved in carcinogenesis; therefore, we hypothesized that altered APA contributes to progression of colorectal cancer, and that measurement of APA may lead to discovery of novel biomarkers.

EXPERIMENTAL DESIGN

We used next-generation sequencing to directly measure global patterns of APA changes during colorectal carcinoma progression in 15 human patient samples. Results were validated in a larger cohort of 50 patients, including 5 normal/carcinoma pairs from individuals.

RESULTS

We discovered numerous genes presenting progressive changes in APA. Genes undergoing untranslated region (3'UTR) shortening were enriched for functional groups such as cell-cycle and nucleic acid-binding and processing factors, and those undergoing 3'UTR lengthening or alternative 3'UTR usage were enriched for categories such as cell-cell adhesion and extracellular matrix. We found indications that APA changes result from differential processing of transcripts because of increased expression of cleavage and polyadenylation factors. Quantitative PCR analysis in a larger series of human patient samples, including matched pairs, confirmed APA changes in DMKN, PDXK, and PPIE genes.

CONCLUSIONS

Our results suggest that genes undergoing altered APA during human cancer progression may be useful novel biomarkers and potentially targeted for disease prevention and treatment. We propose that the strategy presented here may be broadly useful in discovery of novel biomarkers for other types of cancer and human disease.