Altered oncogene expression in cancer cells causes loss of redox homeostasis resulting in oxidative DNA damage, e.g. 8-oxoguanine (8-oxoG), repaired by base excision repair (BER). PARP1 coordinates BER and relies on the upstream 8-oxoguanine-DNA glycosylase (OGG1) to recognise and excise 8-oxoG. Here we hypothesize that OGG1 may represent an attractive target to exploit reactive oxygen species (ROS) elevation in cancer. Although OGG1 depletion is well tolerated in non-transformed cells, we report here that OGG1 depletion obstructs A3 T-cell lymphoblastic acute leukemia growth in vitro and in vivo, validating OGG1 as a potential anti-cancer target. In line with this hypothesis, we show that OGG1 inhibitors (OGG1i) target a wide range of cancer cells, with a favourable therapeutic index compared to non-transformed cells. Mechanistically, OGG1i and shRNA depletion cause S-phase DNA damage, replication stress and proliferation arrest or cell death, representing a novel mechanistic approach to target cancer. This study adds OGG1 to the list of BER factors, e.g. PARP1, as potential targets for cancer treatment.

YAP, the key protein effector of the Hippo pathway, is a transcriptional co-activator that controls the expression of cell cycle genes, promotes cell growth and proliferation and regulates organ size. YAP modulates gene transcription by binding to distal enhancers, but the mechanisms of gene regulation by YAP-bound enhancers remain poorly understood. Here we show that constitutive active YAP5SA leads to widespread changes in chromatin accessibility in untransformed MCF10A cells. Newly accessible regions include YAP-bound enhancers that mediate activation of cycle genes regulated by the Myb-MuvB (MMB) complex. By CRISPR-interference we identify a role for YAP-bound enhancers in phosphorylation of Pol II at Ser5 at MMB-regulated promoters, extending previously published studies that suggested YAP primarily regulates the pause-release step and transcriptional elongation. YAP5SA also leads to less accessible 'closed' chromatin regions, which are not directly YAP-bound but which contain binding motifs for the p53 family of transcription factors. Diminished accessibility at these regions is, at least in part, a consequence of reduced expression and chromatin-binding of the p53 family member ΔNp63 resulting in downregulation of ΔNp63-target genes and promoting YAP-mediated cell migration. In summary, our studies uncover changes in chromatin accessibility and activity that contribute to the oncogenic activities of YAP.

Efficient strategies for precise genome editing in human-induced pluripotent cells (hiPSCs) will enable sophisticated genome engineering for research and clinical purposes. The development of programmable sequence-specific nucleases such as Transcription Activator-Like Effectors Nucleases (TALENs) and Cas9-gRNA allows genetic modifications to be made more efficiently at targeted sites of interest. However, many opportunities remain to optimize these tools and to enlarge their spheres of application. We present several improvements: First, we developed functional re-coded TALEs (reTALEs), which not only enable simple one-pot TALE synthesis but also allow TALE-based applications to be performed using lentiviral vectors. We then compared genome-editing efficiencies in hiPSCs mediated by 15 pairs of reTALENs and Cas9-gRNA targeting CCR5 and optimized ssODN design in conjunction with both methods for introducing specific mutations. We found Cas9-gRNA achieved 7-8× higher non-homologous end joining efficiencies (3%) than reTALENs (0.4%) and moderately superior homology-directed repair efficiencies (1.0 versus 0.6%) when combined with ssODN donors in hiPSCs. Using the optimal design, we demonstrated a streamlined process to generated seamlessly genome corrected hiPSCs within 3 weeks.

Epitope tagging is a powerful and commonly used approach for studying the physical properties of proteins and their functions and localization in eukaryotic cells. In the case of Saccharomyces cerevisiae, it has been possible to exploit the high efficiency of homologous recombination to tag proteins by modifying their endogenous genes, making it possible to tag virtually every endogenous gene and perform genome-wide proteomics experiments. However, due to the relative inefficiency of homologous recombination in cultured human cells, epitope-tagging approaches have been limited to ectopically expressed transgenes, with the attendant limitations of their nonphysiological transcriptional regulation and levels of expression. To overcome this limitation, a modification and extension of adeno-associated virus-mediated human somatic cell gene targeting technology is described that makes it possible to simply and easily create an endogenous epitope tag in the same way that it is possible to knock out a gene. Using this approach, we have created and validated human cell lines with epitope-tagged alleles of two cancer-related genes in a variety of untransformed and transformed human cell lines. This straightforward approach makes it possible to study the physical and biological properties of endogenous proteins in human cells without the need for specialized antibodies for individual proteins of interest.

Alternative splicing is critical for development; however, its role in the specification of the three embryonic germ layers is poorly understood. By performing RNA-Seq on human embryonic stem cells (hESCs) and derived definitive endoderm, cardiac mesoderm, and ectoderm cell lineages, we detect distinct alternative splicing programs associated with each lineage. The most prominent splicing program differences are observed between definitive endoderm and cardiac mesoderm. Integrative multi-omics analyses link each program with lineage-enriched RNA binding protein regulators, and further suggest a widespread role for Quaking (QKI) in the specification of cardiac mesoderm. Remarkably, knockout of QKI disrupts the cardiac mesoderm-associated alternative splicing program and formation of myocytes. These changes arise in part through reduced expression of BIN1 splice variants linked to cardiac development. Mechanistically, we find that QKI represses inclusion of exon 7 in BIN1 pre-mRNA via an exonic ACUAA motif, and this is concomitant with intron removal and cleavage from chromatin. Collectively, our results uncover alternative splicing programs associated with the three germ lineages and demonstrate an important role for QKI in the formation of cardiac mesoderm.

Staufen1 (Stau1) is a ribonucleic acid (RNA)-binding protein involved in the post-transcriptional regulation of gene expression. Recent studies indicate that Stau1-bound messenger RNAs (mRNAs) mainly code for proteins involved in transcription and cell cycle control. Consistently, we report here that Stau1 abundance fluctuates through the cell cycle in HCT116 and U2OS cells: it is high from the S phase to the onset of mitosis and rapidly decreases as cells transit through mitosis. Stau1 down-regulation is mediated by the ubiquitin-proteasome system and the E3 ubiquitin ligase anaphase promoting complex/cyclosome (APC/C). Stau1 interacts with the APC/C co-activators Cdh1 and Cdc20 via its first 88 N-terminal amino acids. The importance of controlling Stau155 levels is underscored by the observation that its overexpression affects mitosis entry and impairs proliferation of transformed cells. Microarray analyses identified 275 Stau1(55)-bound mRNAs in prometaphase cells, an early mitotic step that just precedes Stau1 degradation. Interestingly, several of these mRNAs are more abundant in Stau155-containing complexes in cells arrested in prometaphase than in asynchronous cells. Our results point out for the first time to the possibility that Stau1 participates in a mechanism of post-transcriptional regulation of gene expression that is linked to cell cycle progression in cancer cells.

Normalization of single cell RNA-seq data remains a challenging task. The performance of different methods can vary greatly between datasets when unwanted factors and biology are associated. Most normalization methods also only remove the effects of unwanted variation for the cell embedding but not from gene-level data typically used for differential expression (DE) analysis to identify marker genes. We propose RUV-III-NB, a method that can be used to remove unwanted variation from both the cell embedding and gene-level counts. Using pseudo-replicates, RUV-III-NB explicitly takes into account potential association with biology when removing unwanted variation. The method can be used for both UMI or read counts and returns adjusted counts that can be used for downstream analyses such as clustering, DE and pseudotime analyses. Using published datasets with different technological platforms, kinds of biology and levels of association between biology and unwanted variation, we show that RUV-III-NB manages to remove library size and batch effects, strengthen biological signals, improve DE analyses, and lead to results exhibiting greater concordance with independent datasets of the same kind. The performance of RUV-III-NB is consistent and is not sensitive to the number of factors assumed to contribute to the unwanted variation.

Non-coding RNA regulatory elements are important for viral replication, making them promising targets for therapeutic intervention. However, regulatory RNA is challenging to detect and characterise using classical structure-function assays. Here, we present in cell Mutational Interference Mapping Experiment (in cell MIME) as a way to define RNA regulatory landscapes at single nucleotide resolution under native conditions. In cell MIME is based on (i) random mutation of an RNA target, (ii) expression of mutated RNA in cells, (iii) physical separation of RNA into functional and non-functional populations, and (iv) high-throughput sequencing to identify mutations affecting function. We used in cell MIME to define RNA elements within the 5' region of the HIV-1 genomic RNA (gRNA) that are important for viral replication in cells. We identified three distinct RNA motifs controlling intracellular gRNA production, and two distinct motifs required for gRNA packaging into virions. Our analysis reveals the 73AAUAAA78 polyadenylation motif within the 5' PolyA domain as a dual regulator of gRNA production and gRNA packaging, and demonstrates that a functional polyadenylation signal is required for viral packaging even though it negatively affects gRNA production.

Advances in transcriptomics have led to the discovery of a large number of long intergenic non-coding RNAs (lincRNAs), which are now recognized as important regulators of diverse cellular processes. Although originally thought to be non-coding, recent studies have revealed that many lincRNAs are bound by ribosomes, with a few lincRNAs even having ability to generate micropeptides. The question arises: how widespread the translation of lincRNAs may be and whether such translation is likely to be functional. To better understand biological relevance of lincRNA translation, we systematically characterized lincRNAs with ribosome occupancy by the expression, structural, sequence, evolutionary and functional features for eight human cell lines, revealed that lincRNAs with ribosome occupancy have remarkably distinctive properties compared with those without ribosome occupancy, indicating that translation has important biological implication in categorizing and annotating lincRNAs. Further analysis revealed lincRNAs exhibit remarkable cell-type specificity with differential translational repertoires and substantial discordance in functionality. Collectively, our analyses provide the first attempt to characterize global and cell-type specific properties of translation of lincRNAs in human cells, highlighting that translation of lincRNAs has clear molecular, evolutionary and functional implications. This study will facilitate better understanding of the diverse functions of lincRNAs.

Advanced sequencing technologies have generated a plethora of data for many chromatin marks in multiple tissues and cell types, yet there is lack of a generalized tool for optimal utility of those data. A major challenge is to quantitatively model the epigenetic dynamics across both the genome and many cell types for understanding their impacts on differential gene regulation and disease. We introduce IDEAS, an integrative and discriminative epigenome annotation system, for jointly characterizing epigenetic landscapes in many cell types and detecting differential regulatory regions. A key distinction between our method and existing state-of-the-art algorithms is that IDEAS integrates epigenomes of many cell types simultaneously in a way that preserves the position-dependent and cell type-specific information at fine scales, thereby greatly improving segmentation accuracy and producing comparable annotations across cell types.

E2F transcription factors (E2F1-8) are known to coordinately regulate the expression of a plethora of target genes, including those coding for microRNAs (miRNAs), to control cell cycle progression. Recent work has described the atypical E2F factor E2F7 as a transcriptional repressor of cell cycle-related protein-coding genes. However, the contribution of E2F7 to miRNA gene expression during the cell cycle has not been defined. We have performed a genome-wide RNA sequencing analysis to identify E2F7-regulated miRNAs and show that E2F7 plays as a major role in the negative regulation of a set of miRNAs that promote cellular proliferation. We provide mechanistic evidence for an interplay between E2F7 and the canonical E2F factors E2F1-3 in the regulation of multiple miRNAs. We show that miR-25, -26a, -27b, -92a and -7 expression is controlled at the transcriptional level by the antagonistic activity of E2F7 and E2F1-3. By contrast, let-7 miRNA expression is controlled indirectly through a novel E2F/c-MYC/LIN28B axis, whereby E2F7 and E2F1-3 modulate c-MYC and LIN28B levels to impact let-7 miRNA processing and maturation. Taken together, our data uncover a new regulatory network involving transcriptional and post-transcriptional mechanisms controlled by E2F7 to restrain cell cycle progression through repression of proliferation-promoting miRNAs.

The delivery of large DNA vectors (>100 000 bp) remains a limiting step in the engineering of mammalian cells and the development of human artificial chromosomes (HACs). Yeast is commonly used to assemble genetic constructs in the megabase size range, and has previously been used to transfer constructs directly into cultured cells. We improved this method to efficiently deliver large (1.1 Mb) synthetic yeast centromeric plasmids (YCps) to cultured cell lines at rates similar to that of 12 kb YCps. Synchronizing cells in mitosis improved the delivery efficiency by 10-fold and a statistical design of experiments approach was employed to boost the vector delivery rate by nearly 300-fold from 1/250 000 to 1/840 cells, and subsequently optimize the delivery process for multiple mammalian, avian, and insect cell lines. We adapted this method to rapidly deliver a 152 kb herpes simplex virus 1 genome cloned in yeast into mammalian cells to produce infectious virus.

N6-methyladenosine (m6A) is a reversible and dynamic RNA modification in eukaryotes. However, how cells establish cell-specific m6A methylomes is still poorly understood. Here, we developed a computational framework to systematically identify cell-specific trans regulators of m6A through integrating gene expressions, binding targets and binding motifs of large number of RNA binding proteins (RBPs) with a co-methylation network constructed using large-scale m6A methylomes across diverse cell states. We applied the framework and successfully identified 32 high-confidence m6A regulators that modulated the variable m6A sites away from stop codons in a cell-specific manner. To validate them, we knocked down three regulators respectively and found two of them (TRA2A and CAPRIN1) selectively promoted the methylations of the m6A sites co-localized with their binding targets on RNAs through physical interactions with the m6A writers. Knockdown of TRA2A increased the stabilities of the RNAs with TRA2A bound near the m6A sites and decreased the viability of cells. The successful identification of m6A regulators demonstrates a powerful and widely applicable strategy to elucidate the cell-specific m6A regulators. Additionally, our discovery of pervasive trans-acting regulating of m6A provides novel insights into the mechanisms by which spatial and temporal dynamics of m6A methylomes are established.

The three-dimensional (3D) structures of many biomacromolecules have been solved to reveal the functions of these molecules. However, these 3D structures have rarely been applied to constructing efficient molecular devices that function in living cells. Here, we demonstrate a 3D structure-based molecular design principle for constructing short hairpin RNA (shRNA)-mediated genetic information converters; these converters respond to specific proteins and trigger the desired gene expression by modulating the function of the RNA-processing enzyme Dicer. The inhibitory effect on Dicer cleavage against the shRNA designed to specifically bind to U1A spliceosomal protein was correlated with the degree of steric hindrance between Dicer and the shRNA-protein complex in vitro: The level of the hindrance was predicted based on the models. Moreover, the regulation of gene expression was achieved by using the shRNA converters designed to bind to the target U1A or nuclear factor-κB (NF-κB) p50 proteins expressed in human cells. The 3D molecular design approach is widely applicable for developing new devices in synthetic biology.

The molecular mechanism of human mitochondrial translation has yet to be fully described. We are particularly interested in understanding the process of translational termination and ribosome recycling in the mitochondrion. Several candidates have been implicated, for which subcellular localization and characterization have not been reported. Here, we show that the putative mitochondrial recycling factor, mtRRF, is indeed a mitochondrial protein. Expression of human mtRRF in fission yeast devoid of endogenous mitochondrial recycling factor suppresses the respiratory phenotype. Further, human mtRRF is able to associate with Escherichia coli ribosomes in vitro and can associate with mitoribosomes in vivo. Depletion of mtRRF in human cell lines is lethal, initially causing profound mitochondrial dysmorphism, aggregation of mitoribosomes, elevated mitochondrial superoxide production and eventual loss of OXPHOS complexes. Finally, mtRRF was shown to co-immunoprecipitate a large number of mitoribosomal proteins attached to other mitochondrial proteins, including putative members of the mitochondrial nucleoid.

Cells possess functional diversity hierarchically. However, most single-cell analyses neglect the nested structures while detecting and visualizing the functional diversity. Here, we incorporate cell hierarchy to study functional diversity at subpopulation, club (i.e., sub-subpopulation), and cell layers. Accordingly, we implement a package, SEAT, to construct cell hierarchies utilizing structure entropy by minimizing the global uncertainty in cell-cell graphs. With cell hierarchies, SEAT deciphers functional diversity in 36 datasets covering scRNA, scDNA, scATAC, and scRNA-scATAC multiome. First, SEAT finds optimal cell subpopulations with high clustering accuracy. It identifies cell types or fates from omics profiles and boosts accuracy from 0.34 to 1. Second, SEAT detects insightful functional diversity among cell clubs. The hierarchy of breast cancer cells reveals that the specific tumor cell club drives AREG-EGFT signaling. We identify a dense co-accessibility network of cis-regulatory elements specified by one cell club in GM12878. Third, the cell order from the hierarchy infers periodic pseudo-time of cells, improving accuracy from 0.79 to 0.89. Moreover, we incorporate cell hierarchy layers as prior knowledge to refine nonlinear dimension reduction, enabling us to visualize hierarchical cell layouts in low-dimensional space.

RAD51, a key factor in homology-directed repair (HDR), has long been considered an attractive target for cancer therapy, but few specific inhibitors have been found. A cell-penetrating, anti-DNA, lupus autoantibody, 3E10, was previously shown to inhibit HDR, sensitize tumors to radiation, and mediate synthetic lethal killing of BRCA2-deficient cancer cells, effects that were initially attributed to its affinity for DNA. However, as the molecular basis for its ability to inhibit DNA repair, we report that 3E10 directly binds to the N-terminus of RAD51, sequesters RAD51 in the cytoplasm, and impedes RAD51 binding to DNA. Further, we generate separation-of-function mutations in the complementarity-determining regions of 3E10 revealing that inhibition of HDR tracks with binding to RAD51 but not to DNA, whereas cell penetration is linked to DNA binding. The consequences of these mutations on putative 3E10 interactions with RAD51 and DNA are correlated with in silico molecular modeling. Taken together, the results identify 3E10 as a novel inhibitor of RAD51 by direct binding, accounting for its ability to suppress HDR and providing the molecular basis to guide pre-clinical development of 3E10 as an anti-cancer agent.

Regulated gene expression is essential for a proper progression through the cell cycle. The transcription factor NF-Y has a fundamental function in transcriptional regulation of cell cycle genes, particularly of G2/M genes. In order to investigate common and distinct functions of NF-Y subunits in cell cycle regulation, NF-YA, NF-YB and NF-YC have been silenced by shRNAs in HCT116 cells. NF-YA loss led to a delay in S-phase progression, DNA damage and apoptosis: we showed the activation of the replication checkpoint, through the recruitment of Δp53 and of the replication proteins PCNA and Mcm7 to chromatin. Differently, NF-YB depletion impaired cells from exiting G2/M, but did not interfere with S-phase progression. Gene expression analysis of NF-YA and NF-YB inactivated cells highlighted a common set of hit genes, as well as a plethora of uncommon genes, unveiling a different effect of NF-Y subunits loss on NF-Y binding to its target genes. Chromatin extracts and ChIP analysis showed that NF-YA depletion was more effective than NF-YB in hitting NF-Y recruitment to CCAAT-promoters. Our data suggest a critical role of NF-Y expression, highlighting that the lack of the single subunits are differently perceived by the cells, which activate diverse cell cycle blocks and signaling pathways.

ATAC-seq is a recently developed method to identify the areas of open chromatin in a cell. These regions usually correspond to active regulatory elements and their location profile is unique to a given cell type. When done at single-cell resolution, ATAC-seq provides an insight into the cell-to-cell variability that emerges from otherwise identical DNA sequences by identifying the variability in the genomic location of open chromatin sites in each of the cells. This paper presents Scasat (single-cell ATAC-seq analysis tool), a complete pipeline to process scATAC-seq data with simple steps. Scasat treats the data as binary and applies statistical methods that are especially suitable for binary data. The pipeline is developed in a Jupyter notebook environment that holds the executable code along with the necessary description and results. It is robust, flexible, interactive and easy to extend. Within Scasat we developed a novel differential accessibility analysis method based on information gain to identify the peaks that are unique to a cell. The results from Scasat showed that open chromatin locations corresponding to potential regulatory elements can account for cellular heterogeneity and can identify regulatory regions that separates cells from a complex population.

As molecular and cellular therapies advance in the clinic, the role of genetic regulation is becoming increasingly important for controlling therapeutic potency and safety. The emerging field of mammalian synthetic biology provides promising tools for the construction of regulatory platforms that can intervene with endogenous pathways and control cell behavior. Recent work has highlighted the development of synthetic biological systems that integrate sensing of molecular signals to regulated therapeutic function in various disease settings. However, the toxicity and limited dosing of currently available molecular inducers have largely inhibited translation to clinical settings. In this work, we developed synthetic microRNA-based genetic systems that are controlled by the pharmaceutical drug leucovorin, which is readily available and safe for prolonged administration in clinical settings. We designed microRNA switches to target endogenous cytokine receptor subunits (IL-2Rβ and γc) that mediate various signaling pathways in T cells. We demonstrate the function of these control systems by effectively regulating T cell proliferation with the drug input. Each control system produced unique functional responses, and combinatorial targeting of multiple receptor subunits exhibited greater repression of cell growth. This work highlights the potential use of drug-responsive genetic control systems to improve the management and safety of cellular therapeutics.