ABSTRACT
Steady-state levels of RNA transcripts are controlled by their rates of synthesis and degradation. Here we used nascent RNA Bru-seq and BruChase-seq to profile RNA dynamics across 16 human cell lines as part of ENCODE4 Deeply Profiled Cell Lines collection. We show that RNA turnover dynamics differ widely between transcripts of different genes and between different classes of RNA. Gene set enrichment analysis (GSEA) revealed that transcripts encoding proteins belonging to the same pathway often show similar turnover dynamics. Furthermore, transcript isoforms show distinct dynamics suggesting that RNA turnover is important in regulating mRNA isoform choice. Finally, splicing across newly made transcripts appears to be cooperative with either all or none type splicing. These data sets generated as part of ENCODE4 illustrate the intricate and coordinated regulation of RNA dynamics in controlling gene expression to allow for the precise coordination of cellular functions.
ABSTRACT
Arising as co-products of canonical gene expression, transcription-associated lincRNAs, such as promoter upstream transcripts (PROMPTs), enhancer RNAs (eRNAs), and readthrough (RT) transcripts, are often regarded as byproducts of transcription, although they may be important for the expression of nearby genes. We identified regions of nascent expression of these lincRNA in 16 human cell lines using Bru-seq techniques, and found distinctly regulated patterns of PROMPT, eRNA, and RT transcription using the diverse biochemical approaches in the ENCODE4 deeply profiled cell lines collection. Transcription of these lincRNAs was influenced by sequence-specific features and the local or 3D chromatin landscape. However, these sequence and chromatin features do not describe the full spectrum of lincRNA expression variability we identify, highlighting the complexity of their regulation. This may suggest that transcription-associated lincRNAs are not merely byproducts, but rather that the transcript itself, or the act of its transcription, is important for genomic function.
ABSTRACT
Chromosomal structural variation (SV) encompasses a heterogenous class of genetic variants that exerts strong influences on human health and disease. Despite their importance, many structural variants (SVs) have remained poorly characterized at even a basic level, a discrepancy predicated upon the technical limitations of prior genomic assays. However, recent advances in genomic technology can identify and localize SVs accurately, opening new questions regarding SV risk factors and their impacts in humans. Here, we first define and classify human SVs and their generative mechanisms, highlighting characteristics leveraged by various SV assays. We next examine the first-ever gapless assembly of the human genome and the technical process of assembling it, which required third-generation sequencing technologies to resolve structurally complex loci. The new portions of that "telomere-to-telomere" and subsequent pangenome assemblies highlight aspects of SV biology likely to develop in the near-term. We consider the strengths and limitations of the most promising new SV technologies and when they or longstanding approaches are best suited to meeting salient goals in the study of human SV in population-scale genomics research, clinical, and public health contexts. It is a watershed time in our understanding of human SV when new approaches are expected to fundamentally change genomic applications.
Subject(s)
Genomic Structural Variation , Genomics , Humans , Chromosome Aberrations , Genome, Human , Sequence Analysis, DNAABSTRACT
Error-corrected sequencing of genomic targets enriched by probe-based capture has become a standard approach for detecting single-nucleotide variants (SNVs) and small insertion/deletions (indels) present at very low variant allele frequencies. Less attention has been given to comparable strategies for rare structural variant (SV) junctions, where different error mechanisms must be addressed. Working from samples with known SV properties, we demonstrate that duplex sequencing (DuplexSeq), which demands confirmation of variants on both strands of a source DNA molecule, eliminates false SV junctions arising from chimeric PCR. DuplexSeq could not address frequent intermolecular ligation artifacts that arise during Y-adapter addition prior to strand denaturation without requiring multiple source molecules. In contrast, tagmentation libraries coupled with data filtering based on strand family size greatly reduced both artifact classes and enabled efficient and specific detection of single-molecule SV junctions. The throughput of SV capture sequencing (svCapture) and base-level accuracy of DuplexSeq provided detailed views of the microhomology profile and limited occurrence of de novo SNVs near the junctions of hundreds of newly created SVs, suggesting end joining as a possible formation mechanism. The open source svCapture pipeline enables rare SV detection as a routine addition to SNVs/indels in properly prepared capture sequencing libraries.
ABSTRACT
Retrotransposons are selfish genetic elements that encode an enzyme, reverse transcriptase (RT), which converts the element-encoded RNA into DNA prior to or during genomic integration. New studies provide compelling evidence that a bacterial group II intron-like RT has adapted enzymatic activities associated with RTs to function in host DNA repair.
Subject(s)
RNA-Directed DNA Polymerase , Retroelements , DNA Repair , DNA Transposable Elements/genetics , Introns , RNA , RNA-Directed DNA Polymerase/genetics , RNA-Directed DNA Polymerase/metabolismABSTRACT
The vacuole/lysosome plays essential roles in the growth and proliferation of many eukaryotic cells via the activation of target of rapamycin complex 1 (TORC1). Moreover, the yeast vacuole/lysosome is necessary for progression of the cell division cycle, in part via signaling through the TORC1 pathway. Here, we show that an essential cyclin-dependent kinase, Bur1, plays a critical role in cell cycle progression in cooperation with TORC1. A mutation in BUR1 combined with a defect in vacuole inheritance shows a synthetic growth defect. Importantly, the double mutant, as well as a bur1-267 mutant on its own, has a severe defect in cell cycle progression from G1 phase. In further support that BUR1 functions with TORC1, mutation of bur1 alone results in high sensitivity to rapamycin, a TORC1 inhibitor. Mechanistic insight for Bur1 function comes from the findings that Bur1 directly phosphorylates Sch9, a target of TORC1, and that both Bur1 and TORC1 are required for the activation of Sch9. Together, these discoveries suggest that multiple signals converge on Sch9 to promote cell cycle progression.
Subject(s)
Saccharomyces cerevisiae Proteins , Saccharomyces cerevisiae , Vacuoles , Cell Cycle/genetics , Cyclin-Dependent Kinases/genetics , Cyclin-Dependent Kinases/metabolism , Mechanistic Target of Rapamycin Complex 1/genetics , Mechanistic Target of Rapamycin Complex 1/metabolism , Saccharomyces cerevisiae/genetics , Saccharomyces cerevisiae/metabolism , Saccharomyces cerevisiae Proteins/genetics , Saccharomyces cerevisiae Proteins/metabolism , Transcription Factors , Vacuoles/metabolismABSTRACT
Resection of the 5'-terminated strand at DNA double-strand breaks (DSBs) is the critical regulated step in the transition to homologous recombination. Recent studies have described a multi-step model of DSB resection where endonucleolytic cleavage mediated by Mre11 and Sae2 leads to further degradation mediated by redundant pathways catalyzed by Exo1 and Sgs1/Dna2. These models have not been well tested at mitotic DSBs in vivo because most methods used to monitor resection cannot precisely map early cleavage events. Here we report resection monitoring with high-throughput sequencing using molecular identifiers, allowing exact counting of cleaved 5' ends at base resolution. Mutant strains, including exo1Δ, mre11-H125N and exo1Δ sgs1Δ, revealed a major Mre11-dependent cleavage position 60-70 bp from the DSB end whose exact position depended on local sequence. They further revealed an Exo1-dependent pause point approximately 200 bp from the DSB. Suppressing resection extension in exo1Δ sgs1Δ yeast exposed a footprint of regions where cleavage was restricted within 119 bp of the DSB. These results provide detailed in vivo views of prevailing models of DSB resection and extend them to show the combined influence of sequence specificity and access restrictions on Mre11 and Exo1 nucleases.
Subject(s)
DNA Breaks, Double-Stranded , Exodeoxyribonucleases/metabolism , Fungal Proteins/metabolism , MRE11 Homologue Protein/metabolism , Mitosis/genetics , Recombinational DNA Repair , Alleles , Base Sequence , DNA/chemistry , DNA End-Joining Repair , Exodeoxyribonucleases/genetics , Fungal Proteins/physiology , Gene Deletion , MRE11 Homologue Protein/physiology , RecQ Helicases/genetics , Saccharomycetales/enzymology , Saccharomycetales/geneticsABSTRACT
Impaired replication progression leads to de novo copy number variant (CNV) formation at common fragile sites (CFSs). We previously showed that these hotspots for genome instability reside in late-replicating domains associated with large transcribed genes and provided indirect evidence that transcription is a factor in their instability. Here, we compared aphidicolin (APH)-induced CNV and CFS frequency between wild-type and isogenic cells in which FHIT gene transcription was ablated by promoter deletion. Two promoter-deletion cell lines showed reduced or absent CNV formation and CFS expression at FHIT despite continued instability at the NLGN1 control locus. APH treatment led to critical replication delays that remained unresolved in G2/M in the body of many, but not all, large transcribed genes, an effect that was reversed at FHIT by the promoter deletion. Altering RNase H1 expression did not change CNV induction frequency and DRIP-seq showed a paucity of R-loop formation in the central regions of large genes, suggesting that R-loops are not the primary mediator of the transcription effect. These results demonstrate that large gene transcription is a determining factor in replication stress-induced genomic instability and support models that CNV hotspots mainly result from the transcription-dependent passage of unreplicated DNA into mitosis.
Subject(s)
Acid Anhydride Hydrolases/genetics , DNA Copy Number Variations , DNA Replication , Neoplasm Proteins/genetics , Transcription, Genetic , Acid Anhydride Hydrolases/biosynthesis , Animals , Aphidicolin/pharmacology , Cell Line , Chromosome Fragile Sites , Genetic Loci , Humans , Mice , Mutation , Neoplasm Proteins/biosynthesis , Promoter Regions, Genetic , R-Loop Structures , Ribonuclease H/metabolism , Stress, PhysiologicalABSTRACT
Pre-mRNA splicing is carried out by the spliceosome and involves splice site recognition, removal of introns, and ligation of exons. Components of the spliceosome have been shown to interact with the elongating RNA polymerase II (RNAPII) which is thought to allow splicing to occur concurrently with transcription. However, little is known about the regulation and efficiency of co-transcriptional splicing in human cells. In this study, we used Bru-seq and BruChase-seq to determine the co-transcriptional splicing efficiencies of 17,000 introns expressed across 6 human cell lines. We found that less than half of all introns across these 6 cell lines were co-transcriptionally spliced. Splicing efficiencies for individual introns showed variations across cell lines, suggesting that splicing may be regulated in a cell-type specific manner. Moreover, the splicing efficiency of introns varied within genes. The efficiency of co-transcriptional splicing did not correlate with gene length, intron position, splice site strengths, or the intron/neighboring exons GC content. However, we identified binding signals from multiple RNA binding proteins (RBPs) that correlated with splicing efficiency, including core spliceosomal machinery components-such as SF3B4, U2AF1 and U2AF2 showing higher binding signals in poorly spliced introns. In addition, multiple RBPs, such as BUD13, PUM1 and SND1, showed preferential binding in exons that flank introns with high splicing efficiencies. The nascent RNA splicing patterns presented here across multiple cell types add to our understanding of the complexity in RNA splicing, wherein RNA-binding proteins may play important roles in determining splicing outcomes in a cell type- and intron-specific manner.
ABSTRACT
MicroRNAs (miRNAs) are key contributors to gene regulatory networks. Because miRNAs are processed from RNA polymerase II transcripts, insight into miRNA regulation requires a comprehensive understanding of the regulation of primary miRNA transcripts. We used Bru-seq nascent RNA sequencing and hidden Markov model segmentation to map primary miRNA transcription units (TUs) across 32 human cell lines, allowing us to describe TUs encompassing 1443 miRNAs from miRBase and 438 from MirGeneDB. We identified TUs for 61 miRNAs with an unknown CAGE TSS signal for MirGeneDB miRNAs. Many primary transcripts containing miRNA sequences failed to generate mature miRNAs, suggesting that miRNA biosynthesis is under both transcriptional and post-transcriptional control. In addition to constitutive and cell-type specific TU expression regulated by differential promoter usage, miRNA synthesis can be regulated by transcription past polyadenylation sites (transcriptional read through) and promoter divergent transcription (PROMPTs). We identified 197 miRNA TUs with novel promoters, 97 with transcriptional read-throughs and 3 miRNA TUs that resemble PROMPTs in at least one cell line. The miRNA TU annotation data resource described here reveals a greater complexity in miRNA regulation than previously known and provides a framework for identifying cell-type specific differences in miRNA transcription in cancer and cell transition states.
ABSTRACT
DNA double-strand breaks (DSBs) must be rejoined properly to prevent the occurrence of serious genomic rearrangements associated with many human diseases. Non-homologous end joining (NHEJ) is a DSB repair mechanism known to protect genomic integrity that is also implicated in creating genomic translocations, inversions, deletions, and insertions. We recently investigated the impact of the pre-damage spatial proximity of DSB-bearing loci on the frequency of trans repair by NHEJ and surprisingly found no correlation between them. In this review, we consider various models that might account for these unexpected results. While DSB movement is necessary to explain our findings, many questions remain about the nature and timing of that motion.
Subject(s)
DNA Breaks, Double-Stranded , Translocation, Genetic , Animals , Chromosome Aberrations , DNA End-Joining Repair , Gene Rearrangement , Homologous Recombination , Humans , Models, GeneticABSTRACT
L1 retrotransposon-derived sequences comprise approximately 17% of the human genome. Darwinian selective pressures alter L1 genomic distributions during evolution, confounding the ability to determine initial L1 integration preferences. Here, we generated high-confidence datasets of greater than 88,000 engineered L1 insertions in human cell lines that act as proxies for cells that accommodate retrotransposition in vivo. Comparing these insertions to a null model, in which L1 endonuclease activity is the sole determinant dictating L1 integration preferences, demonstrated that L1 insertions are not significantly enriched in genes, transcribed regions, or open chromatin. By comparison, we provide compelling evidence that the L1 endonuclease disproportionately cleaves predominant lagging strand DNA replication templates, while lagging strand 3'-hydroxyl groups may prime endonuclease-independent L1 retrotransposition in a Fanconi anemia cell line. Thus, acquisition of an endonuclease domain, in conjunction with the ability to integrate into replicating DNA, allowed L1 to become an autonomous, interspersed retrotransposon.
Subject(s)
Long Interspersed Nucleotide Elements/genetics , Retroelements/genetics , Cell Line , Endonucleases/genetics , Endonucleases/metabolism , Genome, Human/genetics , Genome-Wide Association Study/methods , Genomics , HeLa Cells , Humans , Mutagenesis, Insertional/geneticsABSTRACT
DNA double-strand breaks (DSBs) are serious genomic insults that can lead to chromosomal rearrangements if repaired incorrectly. To gain insight into the nuclear mechanisms contributing to these rearrangements, we developed an assay in yeast to measure cis (same site) vs. trans (different site) repair for the majority process of precise nonhomologous end joining (NHEJ). In the assay, the HO endonuclease gene is placed between two HO cut sites such that HO expression is self-terminated upon induction. We further placed an additional cut site in various genomic loci such that NHEJ in trans led to expression of a LEU2 reporter gene. Consistent with prior reports, cis NHEJ was more efficient than trans NHEJ. However, unlike homologous recombination, where spatial distance between a single DSB and donor locus was previously shown to correlate with repair efficiency, trans NHEJ frequency remained essentially constant regardless of the position of the two DSB loci, even when they were on the same chromosome or when two trans repair events were put in competition. Repair of similar DSBs via single-strand annealing of short terminal direct repeats showed substantially higher repair efficiency and trans repair frequency, but still without a strong correlation of trans repair to genomic position. Our results support a model in which yeast cells mobilize, and perhaps compartmentalize, multiple DSBs in a manner that no longer reflects the predamage position of two broken loci.
Subject(s)
DNA Breaks, Double-Stranded , DNA End-Joining Repair/physiology , Gene Expression Regulation, Fungal/physiology , Genetic Loci/physiology , Genome, Fungal/physiology , Saccharomyces cerevisiae , 3-Isopropylmalate Dehydrogenase/biosynthesis , 3-Isopropylmalate Dehydrogenase/genetics , Saccharomyces cerevisiae/genetics , Saccharomyces cerevisiae/metabolism , Saccharomyces cerevisiae Proteins/biosynthesis , Saccharomyces cerevisiae Proteins/geneticsABSTRACT
Clathrin is a major coat protein involved in vesicle formation during endocytosis and transport in the endosomal/trans Golgi system. Clathrin is required for normal growth of yeast (Saccharomyces cerevisiae) and in some genetic backgrounds deletion of the clathrin heavy chain gene (CHC1) is lethal. Our lab defined a locus referred to as " s uppressor of c lathrin d eficiency" (SCD1). In the presence of the scd1-v allele ("v" - viable), yeast cells lacking clathrin heavy chain survive but grow slowly, are morphologically abnormal and have many membrane trafficking defects. In the presence of scd1-i ("i"- inviable), chc1∆ causes lethality. As a strategy to identify SCD1, we used pooled linkage analysis and whole genome sequencing. Here, we report that PAL2 (YHR097C) is the SCD1 locus. pal2∆ is synthetic lethal with chc1∆; whereas a deletion of its paralog, PAL1, is not synthetic lethal with clathrin deficiency. Like Pal1, Pal2 has two NPF motifs that are potential binding sites for EH domain proteins such as the early endocytic factor Ede1, and Pal2 associates with Ede1 Also, GFP-tagged Pal2p localizes to cortical patches containing other immobile phase endocytic coat factors. Overall, our data show that PAL2 is the SCD1 locus and the Pal2 protein has characteristics of an early factor involved in clathrin-mediated endocytosis.
Subject(s)
Clathrin , Endocytosis , Genetic Loci , Receptors, Cell Surface/metabolism , Saccharomyces cerevisiae Proteins/metabolism , Saccharomyces cerevisiae/metabolism , Receptors, Cell Surface/genetics , Saccharomyces cerevisiae/physiology , Saccharomyces cerevisiae Proteins/genetics , Whole Genome SequencingABSTRACT
Copy number variants (CNVs) are important in genome variation and genetic disease, with new mutations arising frequently in the germline and somatic cells. Replication stress caused by aphidicolin and hydroxyurea (HU) is a potent inducer of de novo CNVs in cultured mammalian cells. HU is used extensively for long-term management of sickle cell disease. Here, we examined the effects of HU treatment on germline CNVs in vivo in male mice to explore whether replication stress can act as a CNV mutagen in germline mitotic divisions as in cultured cells and whether this would support a concern for increased CNV mutations in offspring of men treated with HU. Several trials of HU administration were performed by oral gavage and subcutaneous pump, with CNVs characterized in C57BL/6 x C3H/HeJ hybrid mouse offspring by microarray and mate-pair sequencing. HU had a short half-life of ~14 min and a narrow dose window over which studies could be performed while maintaining fertility. Tissue histopathology and reticulocyte micronucleus assays verified that doses had a substantial tissue and genetic toxicity. CNVs were readily detected in offspring that originated in both paternal and maternal mouse strains, as de novo and inherited events. However, HU did not increase CNV formation above baseline levels. These results reveal a high rate of CNV mutagenesis in the mouse germline but do not support the hypothesis that HU would increase CNV formation during mammalian spermatogenesis, perhaps due to highly toxic effects on sperm development or experimental variables related to HU pharmacology in mice. Environ. Mol. Mutagen. 59:698-714, 2018. © 2018 Wiley Periodicals, Inc.
Subject(s)
DNA Copy Number Variations/genetics , DNA Replication/genetics , Germ Cells/drug effects , Hydroxyurea/toxicity , Spermatogenesis/drug effects , Spermatozoa/growth & development , Animals , DNA Copy Number Variations/drug effects , DNA Damage/drug effects , DNA Damage/genetics , DNA Replication/drug effects , Male , Mice , Mice, Inbred C3H , Mice, Inbred C57BLABSTRACT
Ever since initial suggestions that instability at common fragile sites (CFSs) could be responsible for chromosome rearrangements in cancers, CFSs and associated genes have been the subject of numerous studies, leading to questions and controversies about their role and importance in cancer. It is now clear that CFSs are not frequently involved in translocations or other cancer-associated recurrent gross chromosome rearrangements. However, recent studies have provided new insights into the mechanisms of CFS instability, their effect on genome instability, and their role in generating focal copy number alterations that affect the genomic landscape of many cancers.
Subject(s)
Chromosomal Instability , Chromosome Fragile Sites , DNA Copy Number Variations , Neoplasms/genetics , Oncogenes/genetics , Anaphase , Animals , Chromosome Breakage , DNA Breaks, Double-Stranded , DNA Replication , Gene Rearrangement , Humans , MetaphaseABSTRACT
In response to ionizing radiation (IR), cells activate a DNA damage response (DDR) pathway to re-program gene expression. Previous studies using total cellular RNA analyses have shown that the stress kinase ATM and the transcription factor p53 are integral components required for induction of IR-induced gene expression. These studies did not distinguish between changes in RNA synthesis and RNA turnover and did not address the role of enhancer elements in DDR-mediated transcriptional regulation. To determine the contribution of synthesis and degradation of RNA and monitor the activity of enhancer elements following exposure to IR, we used the recently developed Bru-seq, BruChase-seq and BruUV-seq techniques. Our results show that ATM and p53 regulate both RNA synthesis and stability as well as enhancer element activity following exposure to IR. Importantly, many genes in the p53-signaling pathway were coordinately up-regulated by both increased synthesis and RNA stability while down-regulated genes were suppressed either by reduced synthesis or stability. Our study is the first of its kind that independently assessed the effects of ionizing radiation on transcription and post-transcriptional regulation in normal human cells.
Subject(s)
Ataxia Telangiectasia Mutated Proteins/metabolism , Gene Expression Regulation/radiation effects , RNA Processing, Post-Transcriptional , Radiation, Ionizing , Transcription, Genetic , Tumor Suppressor Protein p53/metabolism , Cell Line , DNA Damage/radiation effects , Enhancer Elements, Genetic , Fibroblasts/drug effects , Fibroblasts/metabolism , Gene Expression Profiling , Humans , RNA Stability/radiation effects , Transcriptional ActivationABSTRACT
Human alkyladenine DNA glycosylase (AAG) initiates base excision repair (BER) to guard against mutations by excising alkylated and deaminated purines. Counterintuitively, increased expression of AAG has been implicated in increased rates of spontaneous mutation in microsatellite repeats. This microsatellite mutator phenotype is consistent with a model in which AAG excises bulged (unpaired) bases, altering repeat length. To directly test the role of base excision in AAG-induced mutagenesis, we conducted mutation accumulation experiments in yeast overexpressing different variants of AAG and detected mutations via high-depth genome resequencing. We also developed a new software tool, hp_caller, to perform accurate genotyping at homopolymeric repeat loci. Overexpression of wild-type AAG elevated indel mutations in homopolymeric sequences distributed throughout the genome. However, catalytically inactive variants (E125Q/E125A) caused equal or greater increases in frameshift mutations. These results disprove the hypothesis that base excision is the key step in mutagenesis by overexpressed wild-type AAG. Instead, our results provide additional support for the previously published model wherein overexpressed AAG interferes with the mismatch repair (MMR) pathway. In addition to the above results, we observed a dramatic mutator phenotype for N169S AAG, which has increased rates of excision of undamaged purines. This mutant caused a 10-fold increase in point mutations at G:C base pairs and a 50-fold increase in frameshifts in A:T homopolymers. These results demonstrate that it is necessary to consider the relative activities and abundance of many DNA replication and repair proteins when considering mutator phenotypes, as they are relevant to the development of cancer and its resistance to treatment.