Locus-specific transcription silencing at the FHIT gene suppresses replication stress-induced copy number variant formation and associated replication delay.

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.

Characterization of novel primary miRNA transcription units in human cells using Bru-seq nascent RNA sequencing.

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.

Fragile sites in cancer: more than meets the eye.

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.

Transcriptional and post-transcriptional regulation of the ionizing radiation response by ATM and p53.

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.

Mechanisms of glycosylase induced genomic instability.

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.

Identifying transcription start sites and active enhancer elements using BruUV-seq.

BruUV-seq utilizes UV light to introduce transcription-blocking DNA lesions randomly in the genome prior to bromouridine-labeling and deep sequencing of nascent RNA. By inhibiting transcription elongation, but not initiation, pre-treatment with UV light leads to a redistribution of transcription reads resulting in the enhancement of nascent RNA signal towards the 5'-end of genes promoting the identification of transcription start sites (TSSs). Furthermore, transcripts associated with arrested RNA polymerases are protected from 3'-5' degradation and thus, unstable transcripts such as putative enhancer RNA (eRNA) are dramatically increased. Validation of BruUV-seq against GRO-cap that identifies capped run-on transcripts showed that most BruUV-seq peaks overlapped with GRO-cap signal over both TSSs and enhancer elements. Finally, BruUV-seq identified putative enhancer elements induced by tumor necrosis factor (TNF) treatment concomitant with expression of nearby TNF-induced genes. Taken together, BruUV-seq is a powerful new approach for identifying TSSs and active enhancer elements genome-wide in intact cells.

Rate of elongation by RNA polymerase II is associated with specific gene features and epigenetic modifications.

The rate of transcription elongation plays an important role in the timing of expression of full-length transcripts as well as in the regulation of alternative splicing. In this study, we coupled Bru-seq technology with 5,6-dichlorobenzimidazole 1-ß-D-ribofuranoside (DRB) to estimate the elongation rates of over 2000 individual genes in human cells. This technique, BruDRB-seq, revealed gene-specific differences in elongation rates with a median rate of around 1.5 kb/min. We found that genes with rapid elongation rates showed higher densities of H3K79me2 and H4K20me1 histone marks compared to slower elongating genes. Furthermore, high elongation rates had a positive correlation with gene length, low complexity DNA sequence, and distance from the nearest active transcription unit. Features that negatively correlated with elongation rate included the density of exons, long terminal repeats, GC content of the gene, and DNA methylation density in the bodies of genes. Our results suggest that some static gene features influence transcription elongation rates and that cells may alter elongation rates by epigenetic regulation. The BruDRB-seq technique offers new opportunities to interrogate mechanisms of regulation of transcription elongation.

Use of Bru-Seq and BruChase-Seq for genome-wide assessment of the synthesis and stability of RNA.

Gene expression studies commonly examine total cellular RNA, which only provides information about its steady-state pool of RNA. It remains unclear whether differences in the steady-state reflects variable rates of transcription or RNA degradation. To specifically monitor RNA synthesis and degradation genome-wide, we developed Bru-Seq and BruChase-Seq. These assays are based on metabolic pulse-chase labeling of RNA using bromouridine (Bru). In Bru-Seq, recently labeled RNAs are sequenced to reveal spans of nascent transcription in the genome. In BruChase-Seq, cells are chased in uridine for different periods of time following Bru-labeling, allowing for the isolation of RNA populations of specific ages. Here we describe these methodologies in detail and highlight their usefulness in assessing RNA synthesis and stability as well as splicing kinetics with examples of specific genes from different human cell lines.

Copy number variants are produced in response to low-dose ionizing radiation in cultured cells.

Despite their importance to human genetic variation and disease, little is known about the molecular mechanisms and environmental risk factors that impact copy number variant (CNV) formation. While it is clear that replication stress can lead to de novo CNVs, for example, following treatment of cultured mammalian cells with aphidicolin (APH) and hydroxyurea (HU), the effect of different types of mutagens on CNV induction is unknown. Here we report that ionizing radiation (IR) in the range of 1.5-3.0 Gy effectively induces de novo CNV mutations in cultured normal human fibroblasts. These IR-induced CNVs are found throughout the genome, with the same hotspot regions seen after APH- and HU-induced replication stress. IR produces duplications at a higher frequency relative to deletions than do APH and HU. At most hotspots, these duplications are physically shifted from the regions typically deleted after APH or HU, suggesting different pathways involved in their formation. CNV breakpoint junctions from irradiated samples are characterized by microhomology, blunt ends, and insertions like those seen in spontaneous and APH/HU-induced CNVs and most nonrecurrent CNVs in vivo. The similarity to APH/HU-induced CNVs suggests that low-dose IR induces CNVs through a replication-dependent mechanism, as opposed to replication-independent repair of DSBs. Consistent with this mechanism, a lower yield of CNVs was observed when cells were held for 48 hr before replating after irradiation. These results predict that any environmental DNA damaging agent that impairs replication is capable of creating CNVs.

Genome-wide transcriptional effects of the anti-cancer agent camptothecin.

The anti-cancer drug camptothecin inhibits replication and transcription by trapping DNA topoisomerase I (Top1) covalently to DNA in a "cleavable complex". To examine the effects of camptothecin on RNA synthesis genome-wide we used Bru-Seq and show that camptothecin treatment primarily affected transcription elongation. We also observed that camptothecin increased RNA reads past transcription termination sites as well as at enhancer elements. Following removal of camptothecin, transcription spread as a wave from the 5'-end of genes with no recovery of transcription apparent from RNA polymerases stalled in the body of genes. As a result, camptothecin preferentially inhibited the expression of large genes such as proto-oncogenes, and anti-apoptotic genes while smaller ribosomal protein genes, pro-apoptotic genes and p53 target genes showed relative higher expression. Cockayne syndrome group B fibroblasts (CS-B), which are defective in transcription-coupled repair (TCR), showed an RNA synthesis recovery profile similar to normal fibroblasts suggesting that TCR is not involved in the repair of or RNA synthesis recovery from transcription-blocking Top1 lesions. These findings of the effects of camptothecin on transcription have important implications for its anti-cancer activities and may aid in the design of improved combinatorial treatments involving Top1 poisons.

Recurrent reciprocal RNA chimera involving YPEL5 and PPP1CB in chronic lymphocytic leukemia.

Chronic lymphocytic leukemia (CLL) is the most common form of leukemia in adults in the Western hemisphere. Tumor-specific chromosomal translocations, characteristic findings in several human malignancies that directly lead to malignant transformation, have not been identified in CLL. Using paired-end transcriptome sequencing, we identified recurrent and reciprocal RNA chimeras involving yippee like 5 (YPEL5) and serine/threonine-protein phosphatase PP1-beta-catalytic subunit (PPP1CB) in CLL. Two of seven index cases (28%) harbored the reciprocal RNA chimeras in our initial screening. Using quantitative real-time PCR (q real-time PCR), YPEL5/PPP1CB and PPP1CB/YPEL5 fusion transcripts were detected in 97 of 103 CLL samples (95%) but not in paired normal samples, benign lymphocytes, or various unrelated cancers. Whole-genome sequencing and Southern blotting demonstrated no evidence for a genomic fusion between YPEL5 and PPP1CB. YPEL5/PPP1CB chimera, when introduced into mammalian cells, expressed a truncated PPP1CB protein that demonstrated diminished phosphatase activity. PPP1CB silencing resulted in enhanced proliferation and colony formation of MEC1 and JVM3 cells, implying a role in the pathogenesis of mature B-cell leukemia. These studies uncover a potential role for recurrent RNA chimeras involving phosphatases in the pathogenesis of a common form of leukemia.

Breaking news on fragile sites in cancer.

Chromosome rearrangements in B lymphocytes can be initiated by AID-associated double strand breaks (DSBs), with others arising by unclear mechanisms. A recent study by Barlow and colleagues in Cell reports on genomic regions, termed early replicating fragile sites, that may explain many AID-independent DSBs and creates a compelling link between replication stress, transcription, and chromosome rearrangements.

Coordinated regulation of synthesis and stability of RNA during the acute TNF-induced proinflammatory response.

Steady-state gene expression is a coordination of synthesis and decay of RNA through epigenetic regulation, transcription factors, micro RNAs (miRNAs), and RNA-binding proteins. Here, we present bromouride labeling and sequencing (Bru-Seq) and bromouridine pulse-chase and sequencing (BruChase-Seq) to assess genome-wide changes to RNA synthesis and stability in human fibroblasts at homeostasis and after exposure to the proinflammatory tumor necrosis factor (TNF). The inflammatory response in human cells involves rapid and dramatic changes in gene expression, and the Bru-Seq and BruChase-Seq techniques revealed a coordinated and complex regulation of gene expression both at the transcriptional and posttranscriptional levels. The combinatory analysis of both RNA synthesis and stability using Bru-Seq and BruChase-Seq allows for a much deeper understanding of mechanisms of gene regulation than afforded by the analysis of steady-state total RNA and should be useful in many biological settings.

Harnessing genomics to identify environmental determinants of heritable disease.

Next-generation sequencing technologies can now be used to directly measure heritable de novo DNA sequence mutations in humans. However, these techniques have not been used to examine environmental factors that induce such mutations and their associated diseases. To address this issue, a working group on environmentally induced germline mutation analysis (ENIGMA) met in October 2011 to propose the necessary foundational studies, which include sequencing of parent-offspring trios from highly exposed human populations, and controlled dose-response experiments in animals. These studies will establish background levels of variability in germline mutation rates and identify environmental agents that influence these rates and heritable disease. Guidance for the types of exposures to examine come from rodent studies that have identified agents such as cancer chemotherapeutic drugs, ionizing radiation, cigarette smoke, and air pollution as germ-cell mutagens. Research is urgently needed to establish the health consequences of parental exposures on subsequent generations.

De novo CNV formation in mouse embryonic stem cells occurs in the absence of Xrcc4-dependent nonhomologous end joining.

Spontaneous copy number variant (CNV) mutations are an important factor in genomic structural variation, genomic disorders, and cancer. A major class of CNVs, termed nonrecurrent CNVs, is thought to arise by nonhomologous DNA repair mechanisms due to the presence of short microhomologies, blunt ends, or short insertions at junctions of normal and de novo pathogenic CNVs, features recapitulated in experimental systems in which CNVs are induced by exogenous replication stress. To test whether the canonical nonhomologous end joining (NHEJ) pathway of double-strand break (DSB) repair is involved in the formation of this class of CNVs, chromosome integrity was monitored in NHEJ-deficient Xrcc4(-/-) mouse embryonic stem (ES) cells following treatment with low doses of aphidicolin, a DNA replicative polymerase inhibitor. Mouse ES cells exhibited replication stress-induced CNV formation in the same manner as human fibroblasts, including the existence of syntenic hotspot regions, such as in the Auts2 and Wwox loci. The frequency and location of spontaneous and aphidicolin-induced CNV formation were not altered by loss of Xrcc4, as would be expected if canonical NHEJ were the predominant pathway of CNV formation. Moreover, de novo CNV junctions displayed a typical pattern of microhomology and blunt end use that did not change in the absence of Xrcc4. A number of complex CNVs were detected in both wild-type and Xrcc4(-/-) cells, including an example of a catastrophic, chromothripsis event. These results establish that nonrecurrent CNVs can be, and frequently are, formed by mechanisms other than Xrcc4-dependent NHEJ.

Hydroxyurea induces de novo copy number variants in human cells.

Copy number variants (CNVs) are widely distributed throughout the human genome, where they contribute to genetic variation and phenotypic diversity. Spontaneous CNVs are also a major cause of genetic and developmental disorders and arise frequently in cancer cells. As with all mutation classes, genetic and environmental factors almost certainly increase the risk for new and deleterious CNVs. However, despite the importance of CNVs, there is limited understanding of these precipitating risk factors and the mechanisms responsible for a large percentage of CNVs. Here we report that low doses of hydroxyurea, an inhibitor of ribonucleotide reductase and an important drug in the treatment of sickle cell disease and other diseases induces a high frequency of de novo CNVs in cultured human cells that resemble pathogenic and aphidicolin-induced CNVs in size and breakpoint structure. These CNVs are distributed throughout the genome, with some hotspots of de novo CNV formation. Sequencing revealed that CNV breakpoint junctions are characterized by short microhomologies, blunt ends, and short insertions. These data provide direct experimental support for models of replication-error origins of CNVs and suggest that any agent or condition that leads to replication stress has the potential to induce deleterious CNVs. In addition, they point to a need for further study of the genomic consequences of the therapeutic use of hydroxyurea.

Genetic interactions between HNT3/Aprataxin and RAD27/FEN1 suggest parallel pathways for 5' end processing during base excision repair.

Mutations in Aprataxin cause the neurodegenerative syndrome ataxia oculomotor apraxia type 1. Aprataxin catalyzes removal of adenosine monophosphate (AMP) from the 5' end of a DNA strand, which results from an aborted attempt to ligate a strand break containing a damaged end. To gain insight into which DNA lesions are substrates for Aprataxin action in vivo, we deleted the Saccharomyces cerevisiae HNT3 gene, which encodes the Aprataxin homolog, in combination with known DNA repair genes. While hnt3Delta single mutants were not sensitive to DNA damaging agents, loss of HNT3 caused synergistic sensitivity to H(2)O(2) in backgrounds that accumulate strand breaks with blocked termini, including apn1Delta apn2Delta tpp1Delta and ntg1Delta ntg2Delta ogg1Delta. Loss of HNT3 in rad27Delta cells, which are deficient in long-patch base excision repair (LP-BER), resulted in synergistic sensitivity to H(2)O(2) and MMS, indicating that Hnt3 and LP-BER provide parallel pathways for processing 5' AMPs. Loss of HNT3 also increased the sister chromatid exchange frequency. Surprisingly, HNT3 deletion partially rescued H(2)O(2) sensitivity in recombination-deficient rad51Delta and rad52Delta cells, suggesting that Hnt3 promotes formation of a repair intermediate that is resolved by recombination.

Recruitment and dissociation of nonhomologous end joining proteins at a DNA double-strand break in Saccharomyces cerevisiae.

Nonhomologous end joining (NHEJ) is an important DNA double-strand-break (DSB) repair pathway that requires three protein complexes in Saccharomyces cerevisiae: the Ku heterodimer (Yku70-Yku80), MRX (Mre11-Rad50-Xrs2), and DNA ligase IV (Dnl4-Lif1), as well as the ligase-associated protein Nej1. Here we use chromatin immunoprecipitation from yeast to dissect the recruitment and release of these protein complexes at HO-endonuclease-induced DSBs undergoing productive NHEJ. Results revealed that Ku and MRX assembled at a DSB independently and rapidly after DSB formation. Ligase IV appeared at the DSB later than Ku and MRX and in a strongly Ku-dependent manner. Ligase binding was extensive but slightly delayed in rad50 yeast. Ligase IV binding occurred independently of Nej1, but instead promoted loading of Nej1. Interestingly, dissociation of Ku and ligase from unrepaired DSBs depended on the presence of an intact MRX complex and ATP binding by Rad50, suggesting a possible role of MRX in terminating a NHEJ repair phase. This activity correlated with extended DSB resection, but limited degradation of DSB ends occurred even in MRX mutants with persistently bound Ku. These findings reveal the in vivo assembly of the NHEJ repair complex and shed light on the mechanisms controlling DSB repair pathway utilization.

Abrogation of the Chk1-Pds1 checkpoint leads to tolerance of persistent single-strand breaks in Saccharomyces cerevisiae.

In budding yeast, Apn1, Apn2, Tpp1, and Rad1/Rad10 are important enzymes in the removal of spontaneous DNA lesions. apn1 apn2 rad1 yeast are inviable due to accumulation of abasic sites and strand breaks with 3' blocking lesions. We found that tpp1 apn1 rad1 yeast exhibited slow growth but frequently gave rise to spontaneous slow growth suppressors that segregated as single-gene mutations. Using a candidate gene approach, we identified several tpp1 apn1 rad1 suppressors. Deleting uracil glycosylase suppressed both tpp1 apn1 rad1 and apn1 apn2 rad1 growth defects by reducing the abasic site burden. Mutants affecting the Chk1-Pds1 metaphase-anaphase checkpoint only suppressed tpp1 apn1 rad1 slow growth. In contrast, most S-phase checkpoint mutants were synthetically lethal in a tpp1 apn1 rad1 background. Epistasis analyses showed an additive effect between chk1 and ung1, indicating different mechanisms of suppression. Loss of Chk1 partially restored cell-growth parameters in tpp1 apn1 rad1 yeast, but at the same time exacerbated chromosome instability. We propose a model in which recombinational repair during S phase coupled with failure of the metaphase-anaphase checkpoint allows for tolerance of persistent single-strand breaks at the expense of genome stability.