Systematic Characterization of p53-Regulated Long Noncoding RNAs across Human Cancers Reveals Remarkable Heterogeneity among Different Tumor Types.

The p53 tumor suppressor protein, a sequence-specific DNA binding transcription factor, regulates the expression of a large number of genes, in response to various forms of cellular stress. Although the protein coding target genes of p53 have been well studied, less is known about its role in regulating long noncoding genes and their functional relevance to cancer. Here we report the genome-wide identification of a large set (>1,000) of long noncoding RNAs (lncRNA), which are putative p53 targets in a colon cancer cell line and in human patient datasets from five different common types of cancer. These lncRNAs have not been annotated by other studies of normal unstressed systems. In the colon cancer cell line, a high proportion of these lncRNAs are uniquely induced by different chemotherapeutic agents that activate p53, whereas others are induced by more than one agent tested. Further, subsets of these lncRNAs independently predict overall and disease-free survival of patients across the five different common cancer types. Interestingly, both genetic alterations and patient survival associated with different lncRNAs are unique to each cancer tested, indicating extraordinary tissue-specific variability in the p53 noncoding response. The newly identified noncoding p53 target genes have allowed us to construct a classifier for tumor diagnosis and prognosis. IMPLICATIONS: Our results not only identify myriad p53-regulated long noncoding (lncRNA), they also reveal marked drug-induced, as well as tissue- and tumor-specific heterogeneity in these putative p53 targets and our findings have enabled the construction of robust classifiers for diagnosis and prognosis.

Time-Course Progression of Whole Transcriptome Expression Changes of Trigeminal Ganglia Compared to Dorsal Root Ganglia in Rats Exposed to Nerve Injury.

Mechanisms underlying neuropathic pain (NP) are complex with multiple genes, their interactions, environmental and epigenetic factors being implicated. Transcriptional changes in the trigeminal (TG) and dorsal root (DRG) ganglia have been implicated in the development and maintenance of NP. Despite efforts to unravel molecular mechanisms of NP, many remain unknown. Also, most of the studies focused on the spinal system. Although the spinal and trigeminal systems share some of the molecular mechanisms, differences exist. We used RNA-sequencing technology to identify differentially expressed genes (DEGs) in the TG and DRG at baseline and 3 time points following the infraorbital or sciatic nerve injuries, respectively. Pathway analysis and comparison analysis were performed to identify differentially expressed pathways. Additionally, upstream regulator effects were investigated in the two systems. DEG (differentially expressed genes) analyses identified 3,225 genes to be differentially expressed between TG and DRG in naïve animals, 1,828 genes 4 days post injury, 5,644 at day 8 and 9,777 DEGs at 21 days postinjury. A comparison of top enriched canonical pathways revealed that a number of signaling pathway was significantly inhibited in the TG and activated in the DRG at 21 days postinjury. Finally, CORT upstream regulator was predicted to be inhibited in the TG while expression levels of the CSF1 upstream regulator were significantly elevated in the DRG at 21 days postinjury. This study provides a basis for further in-depth studies investigating transcriptional changes, pathways, and upstream regulation in TG and DRG in rats exposed to peripheral nerve injuries. PERSPECTIVE: Although trigeminal and dorsal root ganglia are homologs of each other, they respond differently to nerve injury and therefore treatment. Activation/inhibition of number of biological pathways appear to be ganglion/system specific suggesting that different approaches might be required to successfully treat neuropathies induced by injuries in spinal and trigeminal systems.

SETX (senataxin), the helicase mutated in AOA2 and ALS4, functions in autophagy regulation.

SETX (senataxin) is an RNA/DNA helicase that has been implicated in transcriptional regulation and the DNA damage response through resolution of R-loop structures. Mutations in SETX result in either of two distinct neurodegenerative disorders. SETX dominant mutations result in a juvenile form of amyotrophic lateral sclerosis (ALS) called ALS4, whereas recessive mutations are responsible for ataxia called ataxia with oculomotor apraxia type 2 (AOA2). How mutations in the same protein can lead to different phenotypes is still unclear. To elucidate AOA2 disease mechanisms, we first examined gene expression changes following SETX depletion. We observed the effects on both transcription and RNA processing, but surprisingly observed decreased R-loop accumulation in SETX-depleted cells. Importantly, we discovered a strong connection between SETX and the macroautophagy/autophagy pathway, reflecting a direct effect on transcription of autophagy genes. We show that SETX depletion inhibits the progression of autophagy, leading to an accumulation of ubiquitinated proteins, decreased ability to clear protein aggregates, as well as mitochondrial defects. Analysis of AOA2 patient fibroblasts also revealed a perturbation of the autophagy pathway. Our work has thus identified a novel function for SETX in the regulation of autophagy, whose modulation may have a therapeutic impact for AOA2.Abbreviations: 3'READS: 3' region extraction and deep sequencing; ACTB: actin beta; ALS4: amyotrophic lateral sclerosis type 4; AOA2: ataxia with oculomotor apraxia type 2; APA: alternative polyadenylation; AS: alternative splicing; ATG7: autophagy-related 7; ATP6V0D2: ATPase H+ transporting V0 subunit D2; BAF: bafilomycin A1; BECN1: beclin 1; ChIP: chromatin IP; Chloro: chloroquine; CPT: camptothecin; DDR: DNA damage response; DNMT1: DNA methyltransferase 1; DRIP: DNA/RNA IP; DSBs: double strand breaks; EBs: embryoid bodies; FTD: frontotemporal dementia; GABARAP: GABA type A receptor-associated protein; GO: gene ontology; HR: homologous recombination; HTT: huntingtin; IF: immunofluorescence; IP: immunoprecipitation; iPSCs: induced pluripotent stem cells; KD: knockdown; MAP1LC3/LC3: microtubule associated protein 1 light chain 3; MN: motor neuron; MTORC1: mechanistic target of rapamycin kinase complex 1; PASS: PolyA Site Supporting; PFA: paraformaldehyde; RNAPII: RNA polymerase II; SCA: spinocerebellar ataxia; SETX: senataxin; SMA: spinal muscular atrophy; SMN1: survival of motor neuron 1, telomeric; SQSTM1/p62: sequestosome 1; TFEB: transcription factor EB; TSS: transcription start site; TTS: transcription termination site; ULK1: unc-51 like autophagy activating kinase 1; WB: western blot; WIPI2: WD repeat domain, phosphoinositide interacting 2; XRN2: 5'-3' exoribonuclease 2.

Widespread transcript shortening through alternative polyadenylation in secretory cell differentiation.

Most eukaryotic genes produce alternative polyadenylation (APA) isoforms. Here we report that, unlike previously characterized cell lineages, differentiation of syncytiotrophoblast (SCT), a cell type critical for hormone production and secretion during pregnancy, elicits widespread transcript shortening through APA in 3'UTRs and in introns. This global APA change is observed in multiple in vitro trophoblast differentiation models, and in single cells from placentas at different stages of pregnancy. Strikingly, the transcript shortening is unrelated to cell proliferation, a feature previously associated with APA control, but instead accompanies increased secretory functions. We show that 3'UTR shortening leads to transcripts with higher mRNA stability, which augments transcriptional activation, especially for genes involved in secretion. Moreover, this mechanism, named secretion-coupled APA (SCAP), is also executed in B cell differentiation to plasma cells. Together, our data indicate that SCAP tailors the transcriptome during formation of secretory cells, boosting their protein production and secretion capacity.

Rapidly Correcting Frameshift Mutations in the Mycobacterium tuberculosis orn Gene Produce Reversible Ethambutol Resistance and Small-Colony-Variant Morphology.

We have identified a previously unknown mechanism of reversible high-level ethambutol (EMB) resistance in Mycobacterium tuberculosis that is caused by a reversible frameshift mutation in the M. tuberculosisorn gene. A frameshift mutation in orn produces the small-colony-variant (SCV) phenotype, but this mutation does not change the MICs of any drug for wild-type M. tuberculosis However, the same orn mutation in a low-level EMB-resistant double embB-aftA mutant (MIC = 8 µg/ml) produces an SCV with an EMB MIC of 32 µg/ml. Reversible resistance is indistinguishable from a drug-persistent phenotype, because further culture of these orn-embB-aftA SCV mutants results in rapid reversion of the orn frameshifts, reestablishing the correct orn open reading frame, returning the culture to normal colony size, and reversing the EMB MIC back to that (8 µg/ml) of the parental strain. Transcriptomic analysis of orn-embB-aftA mutants compared to wild-type M. tuberculosis identified a 27-fold relative increase in the expression of embC, which is a cellular target for EMB. Expression of embC in orn-embB-aftA mutants was also increased 5-fold compared to that in the parental embB-aftA mutant, whereas large-colony orn frameshift revertants of the orn-embB-aftA mutant had levels of embC expression similar to that of the parental embB-aftA strain. Reversible frameshift mutants may contribute to a reversible form of microbiological drug resistance in human tuberculosis.

Phase variation in Mycobacterium tuberculosis glpK produces transiently heritable drug tolerance.

The length and complexity of tuberculosis (TB) therapy, as well as the propensity of Mycobacterium tuberculosis to develop drug resistance, are major barriers to global TB control efforts. M. tuberculosis is known to have the ability to enter into a drug-tolerant state, which may explain many of these impediments to TB treatment. We have identified a mechanism of genetically encoded but rapidly reversible drug tolerance in M. tuberculosis caused by transient frameshift mutations in a homopolymeric tract (HT) of 7 cytosines (7C) in the glpK gene. Inactivating frameshift mutations associated with the 7C HT in glpK produce small colonies that exhibit heritable multidrug increases in minimal inhibitory concentrations and decreases in drug-dependent killing; however, reversion back to a fully drug-susceptible large-colony phenotype occurs rapidly through the introduction of additional insertions or deletions in the same glpK HT region. These reversible frameshift mutations in the 7C HT of M. tuberculosis glpK occur in clinical isolates, accumulate in M. tuberculosis-infected mice with further accumulation during drug treatment, and exhibit a reversible transcriptional profile including induction of dosR and sigH and repression of kstR regulons, similar to that observed in other in vitro models of M. tuberculosis tolerance. These results suggest that GlpK phase variation may contribute to drug tolerance, treatment failure, and relapse in human TB. Drugs effective against phase-variant M. tuberculosis may hasten TB treatment and improve cure rates.

Transcriptome 3'end organization by PCF11 links alternative polyadenylation to formation and neuronal differentiation of neuroblastoma.

Diversification at the transcriptome 3'end is an important and evolutionarily conserved layer of gene regulation associated with differentiation and dedifferentiation processes. Here, we identify extensive transcriptome 3'end-alterations in neuroblastoma, a tumour entity with a paucity of recurrent somatic mutations and an unusually high frequency of spontaneous regression. Utilising extensive RNAi-screening we reveal the landscape and drivers of transcriptome 3'end-diversification, discovering PCF11 as critical regulator, directing alternative polyadenylation (APA) of hundreds of transcripts including a differentiation RNA-operon. PCF11 shapes inputs converging on WNT-signalling, and governs cell cycle, proliferation, apoptosis and neurodifferentiation. Postnatal PCF11 down-regulation induces a neurodifferentiation program, and low-level PCF11 in neuroblastoma associates with favourable outcome and spontaneous tumour regression. Our findings document a critical role for APA in tumorigenesis and describe a novel mechanism for cell fate reprogramming in neuroblastoma with potentially important clinical implications. We provide an interactive data repository of transcriptome-wide APA covering > 170 RNAis, and an APA-network map with regulatory hubs.

The C9ORF72 Gene, Implicated in Amyotrophic Lateral Sclerosis and Frontotemporal Dementia, Encodes a Protein That Functions in Control of Endothelin and Glutamate Signaling.

A GGGGCC repeat expansion in the C9ORF72 (C9) gene is the most common known cause of amyotrophic lateral sclerosis (ALS) and frontotemporal dementia. Several mechanisms have been proposed to account for its toxicity, including the possibility that reduced C9 protein levels contribute to disease. To investigate this possibility, we examined the effects of reduced C9 levels in several cell systems. We first showed that C9 knockdown (KD) in U87 glioblastoma cells results in striking morphological changes, including vacuolization and alterations in cell size. Unexpectedly, RNA analysis revealed changes in expression of many genes, including genes involved in endothelin (EDN) signaling and immune system pathways and multiple glutamate cycling genes (e.g., EAAT2), which were verified in several cell models, including astrocytes and brain samples from C9-positive patients. Consistent with deregulation of the glutamate cycling genes, elevated intracellular glutamate was detected in both KD cells and patient astrocytes. Importantly, levels of mRNAs encoding EDN1 and its receptors, known to be elevated in ALS, were sharply increased by C9 KD, likely resulting from an observed activation of NF-κB signaling and/or a possible role of a C9 isoform in gene control.

Identifying Cellular Nonsense-Mediated mRNA Decay (NMD) Targets: Immunoprecipitation of Phosphorylated UPF1 Followed by RNA Sequencing (p-UPF1 RIP-Seq).

Recent progress in the technology of transcriptome-wide high-throughput sequencing has revealed that nonsense-mediated mRNA decay (NMD) targets ~10% of physiologic transcripts for the purpose of tuning gene expression in response to various environmental conditions. Regardless of the eukaryote studied, NMD requires the ATP-dependent RNA helicase upframeshift 1 (UPF1). It was initially thought that cellular NMD targets could be defined by their binding to steady-state UPF1, which is largely hypophosphorylated. However, the propensity for steady-state UPF1 to bind RNA nonspecifically, coupled with regulated phosphorylation of UPF1 on an NMD target serving as the trigger for NMD, made it clear that it is phosphorylated UPF1 (p-UPF1), rather than steady-state UPF1, that can be used to distinguish cellular NMD targets from cellular RNAs that are not. Here, we describe the immunoprecipitation of p-UPF1 followed by RNA sequencing (p-UPF1 RIP-seq) as a transcriptome-wide approach to define physiologic NMD targets.

MPK1/SLT2 Links Multiple Stress Responses with Gene Expression in Budding Yeast by Phosphorylating Tyr1 of the RNAP II CTD.

The RNA polymerase II largest subunit C-terminal domain consists of repeated YSPTSPS heptapeptides. The role of tyrosine-1 (Tyr1) remains incompletely understood, as, for example, mutating all Tyr1 residues to Phe (Y1F) is lethal in vertebrates but a related mutant has only a mild phenotype in S. pombe. Here we show that Y1F substitution in budding yeast resulted in a strong slow-growth phenotype. The Y1F strain was also hypersensitive to several different cellular stresses that involve MAP kinase signaling. These phenotypes were all linked to transcriptional changes, and we also identified genetic and biochemical interactions between Tyr1 and both transcription initiation and termination factors. Further studies uncovered defects related to MAP kinase I (Slt2) pathways, and we provide evidence that Slt2 phosphorylates Tyr1 in vitro and in vivo. Our study has thus identified Slt2 as a Tyr1 kinase, and in doing so provided links between stress response activation and Tyr1 phosphorylation.

Regulation of gene expression by translation factor eIF5A: Hypusine-modified eIF5A enhances nonsense-mediated mRNA decay in human cells.

Nonsense-mediated mRNA decay (NMD) couples protein synthesis to mRNA turnover. It eliminates defective transcripts and controls the abundance of certain normal mRNAs. Our study establishes a connection between NMD and the translation factor eIF5A (eukaryotic initiation factor 5A) in human cells. eIF5A modulates the synthesis of groups of proteins (the eIF5A regulon), and undergoes a distinctive two-step post-translational modification (hypusination) catalyzed by deoxyhypusine synthase and deoxyhypusine hydroxylase. We show that expression of NMD-susceptible constructs was increased by depletion of the major eIF5A isoform, eIF5A1. NMD was also attenuated when hypusination was inhibited by RNA interference with either of the two eIF5A modifying enzymes, or by treatment with the drugs ciclopirox or deferiprone which inhibit deoxyhypusine hydroxylase. Transcriptome analysis by RNA-Seq identified human genes whose expression is coordinately regulated by eIF5A1, its modifying enzymes, and the pivotal NMD factor, Upf1. Transcripts encoding components of the translation system were highly represented, including some encoding ribosomal proteins controlled by alternative splicing coupled to NMD (AS-NMD). Our findings extend and strengthen the association of eIF5A with NMD, previously inferred in yeast, and show that hypusination is important for this function of human eIF5A. In addition, they advance drug-mediated NMD suppression as a therapeutic opportunity for nonsense-associated diseases. We propose that regulation of mRNA stability contributes to eIF5A's role in selective gene expression.

Comparative analysis of alternative polyadenylation in S. cerevisiae and S. pombe.

Alternative polyadenylation (APA) is a widespread mechanism that generates mRNA isoforms with distinct properties. Here we have systematically mapped and compared cleavage and polyadenylation sites (PASs) in two yeast species, S. cerevisiae and S. pombe Although >80% of the mRNA genes in each species were found to display APA, S. pombe showed greater 3' UTR size differences among APA isoforms than did S. cerevisiae PASs in different locations of gene are surrounded with distinct sequences in both species and are often associated with motifs involved in the Nrd1-Nab3-Sen1 termination pathway. In S. pombe, strong motifs surrounding distal PASs lead to higher abundances of long 3' UTR isoforms than short ones, a feature that is opposite in S. cerevisiae Differences in PAS placement between convergent genes lead to starkly different antisense transcript landscapes between budding and fission yeasts. In both species, short 3' UTR isoforms are more likely to be expressed when cells are growing in nutrient-rich media, although different gene groups are affected in each species. Significantly, 3' UTR shortening in S. pombe coordinates with up-regulation of expression for genes involved in translation during cell proliferation. Using S. pombe strains deficient for Pcf11 or Pab2, we show that reduced expression of 3'-end processing factors lengthens 3' UTR, with Pcf11 having a more potent effect than Pab2. Taken together, our data indicate that APA mechanisms in S. pombe and S. cerevisiae are largely different: S. pombe has many of the APA features of higher species, and Pab2 in S. pombe has a different role in APA regulation than its mammalian homolog, PABPN1.

Distinct regulation of alternative polyadenylation and gene expression by nuclear poly(A) polymerases.

Polyadenylation of nascent RNA by poly(A) polymerase (PAP) is important for 3' end maturation of almost all eukaryotic mRNAs. Most mammalian genes harbor multiple polyadenylation sites (PASs), leading to expression of alternative polyadenylation (APA) isoforms with distinct functions. How poly(A) polymerases may regulate PAS usage and hence gene expression is poorly understood. Here, we show that the nuclear canonical (PAPα and PAPγ) and non-canonical (Star-PAP) PAPs play diverse roles in PAS selection and gene expression. Deficiencies in the PAPs resulted in perturbations of gene expression, with Star-PAP impacting lowly expressed mRNAs and long-noncoding RNAs to the greatest extent. Importantly, different PASs of a gene are distinctly regulated by different PAPs, leading to widespread relative expression changes of APA isoforms. The location and surrounding sequence motifs of a PAS appear to differentiate its regulation by the PAPs. We show Star-PAP-specific PAS usage regulates the expression of the eukaryotic translation initiation factor EIF4A1, the tumor suppressor gene PTEN and the long non-coding RNA NEAT1. The Star-PAP-mediated APA of PTEN is essential for DNA damage-induced increase of PTEN protein levels. Together, our results reveal a PAS-guided and PAP-mediated paradigm for gene expression in response to cellular signaling cues.

Transcription elongation rate has a tissue-specific impact on alternative cleavage and polyadenylation in Drosophila melanogaster.

Alternative polyadenylation (APA) is a mechanism that generates multiple mRNA isoforms with different 3'UTRs and/or coding sequences from a single gene. Here, using 3' region extraction and deep sequencing (3'READS), we have systematically mapped cleavage and polyadenylation sites (PASs) in Drosophila melanogaster, expanding the total repertoire of PASs previously identified for the species, especially those located in A-rich genomic sequences. Cis-element analysis revealed distinct sequence motifs around fly PASs when compared to mammalian ones, including the greater enrichment of upstream UAUA elements and the less prominent presence of downstream UGUG elements. We found that over 75% of mRNA genes in Drosophila melanogaster undergo APA. The head tissue tends to use distal PASs when compared to the body, leading to preferential expression of APA isoforms with long 3'UTRs as well as with distal terminal exons. The distance between the APA sites and intron location of PAS are important parameters for APA difference between body and head, suggesting distinct PAS selection contexts. APA analysis of the RpII215C4 mutant strain, which harbors a mutant RNA polymerase II (RNAPII) with a slower elongation rate, revealed that a 50% decrease in transcriptional elongation rate leads to a mild trend of more usage of proximal, weaker PASs, both in 3'UTRs and in introns, consistent with the "first come, first served" model of APA regulation. However, this trend was not observed in the head, suggesting a different regulatory context in neuronal cells. Together, our data expand the PAS collection for Drosophila melanogaster and reveal a tissue-specific effect of APA regulation by RNAPII elongation rate.

An Mtr4/ZFC3H1 complex facilitates turnover of unstable nuclear RNAs to prevent their cytoplasmic transport and global translational repression.

Many long noncoding RNAs (lncRNAs) are unstable and rapidly degraded in the nucleus by the nuclear exosome. An exosome adaptor complex called NEXT (nuclear exosome targeting) functions to facilitate turnover of some of these lncRNAs. Here we show that knockdown of one NEXT subunit, Mtr4, but neither of the other two subunits, resulted in accumulation of two types of lncRNAs: prematurely terminated RNAs (ptRNAs) and upstream antisense RNAs (uaRNAs). This suggested a NEXT-independent Mtr4 function, and, consistent with this, we isolated a distinct complex containing Mtr4 and the zinc finger protein ZFC3H1. Strikingly, knockdown of either protein not only increased pt/uaRNA levels but also led to their accumulation in the cytoplasm. Furthermore, all pt/uaRNAs examined associated with active ribosomes, but, paradoxically, this correlated with a global reduction in heavy polysomes and overall repression of translation. Our findings highlight a critical role for Mtr4/ZFC3H1 in nuclear surveillance of naturally unstable lncRNAs to prevent their accumulation, transport to the cytoplasm, and resultant disruption of protein synthesis.

Normalization of TAM post-receptor signaling reveals a cell invasive signature for Axl tyrosine kinase.

BACKGROUND: Tyro3, Axl, and Mertk (TAMs) are a family of three conserved receptor tyrosine kinases that have pleiotropic roles in innate immunity and homeostasis and when overexpressed in cancer cells can drive tumorigenesis. METHODS: In the present study, we engineered EGFR/TAM chimeric receptors (EGFR/Tyro3, EGFR/Axl, and EGF/Mertk) with the goals to interrogate post-receptor functions of TAMs, and query whether TAMs have unique or overlapping post-receptor activation profiles. Stable expression of EGFR/TAMs in EGFR-deficient CHO cells afforded robust EGF inducible TAM receptor phosphorylation and activation of downstream signaling. RESULTS: Using a series of unbiased screening approaches, that include kinome-view analysis, phosphor-arrays, RNAseq/GSEA analysis, as well as cell biological and in vivo readouts, we provide evidence that each TAM has unique post-receptor signaling platforms and identify an intrinsic role for Axl that impinges on cell motility and invasion compared to Tyro3 and Mertk. CONCLUSION: These studies demonstrate that TAM show unique post-receptor signatures that impinge on distinct gene expression profiles and tumorigenic outcomes.

Intronic cleavage and polyadenylation regulates gene expression during DNA damage response through U1 snRNA.

The DNA damage response involves coordinated control of gene expression and DNA repair. Using deep sequencing, we found widespread changes of alternative cleavage and polyadenylation site usage on ultraviolet-treatment in mammalian cells. Alternative cleavage and polyadenylation regulation in the 3' untranslated region is substantial, leading to both shortening and lengthening of 3' untranslated regions of genes. Interestingly, a strong activation of intronic alternative cleavage and polyadenylation sites is detected, resulting in widespread expression of truncated transcripts. Intronic alternative cleavage and polyadenylation events are biased to the 5' end of genes and affect gene groups with important functions in DNA damage response and cancer. Moreover, intronic alternative cleavage and polyadenylation site activation during DNA damage response correlates with a decrease in U1 snRNA levels, and is reversible by U1 snRNA overexpression. Importantly, U1 snRNA overexpression mitigates ultraviolet-induced apoptosis. Together, these data reveal a significant gene regulatory scheme in DNA damage response where U1 snRNA impacts gene expression via the U1-alternative cleavage and polyadenylation axis.

The Nrd1-like protein Seb1 coordinates cotranscriptional 3' end processing and polyadenylation site selection.

Termination of RNA polymerase II (RNAPII) transcription is associated with RNA 3' end formation. For coding genes, termination is initiated by the cleavage/polyadenylation machinery. In contrast, a majority of noncoding transcription events in Saccharomyces cerevisiae does not rely on RNA cleavage for termination but instead terminates via a pathway that requires the Nrd1-Nab3-Sen1 (NNS) complex. Here we show that the Schizosaccharomyces pombe ortholog of Nrd1, Seb1, does not function in NNS-like termination but promotes polyadenylation site selection of coding and noncoding genes. We found that Seb1 associates with 3' end processing factors, is enriched at the 3' end of genes, and binds RNA motifs downstream from cleavage sites. Importantly, a deficiency in Seb1 resulted in widespread changes in 3' untranslated region (UTR) length as a consequence of increased alternative polyadenylation. Given that Seb1 levels affected the recruitment of conserved 3' end processing factors, our findings indicate that the conserved RNA-binding protein Seb1 cotranscriptionally controls alternative polyadenylation.

PAF Complex Plays Novel Subunit-Specific Roles in Alternative Cleavage and Polyadenylation.

The PAF complex (Paf1C) has been shown to regulate chromatin modifications, gene transcription, and RNA polymerase II (PolII) elongation. Here, we provide the first genome-wide profiles for the distribution of the entire complex in mammalian cells using chromatin immunoprecipitation and high throughput sequencing. We show that Paf1C is recruited not only to promoters and gene bodies, but also to regions downstream of cleavage/polyadenylation (pA) sites at 3' ends, a profile that sharply contrasted with the yeast complex. Remarkably, we identified novel, subunit-specific links between Paf1C and regulation of alternative cleavage and polyadenylation (APA) and upstream antisense transcription using RNAi coupled with deep sequencing of the 3' ends of transcripts. Moreover, we found that depletion of Paf1C subunits resulted in the accumulation of PolII over gene bodies, which coincided with APA. Depletion of specific Paf1C subunits led to global loss of histone H2B ubiquitylation, although there was little impact of Paf1C depletion on other histone modifications, including tri-methylation of histone H3 on lysines 4 and 36 (H3K4me3 and H3K36me3), previously associated with this complex. Our results provide surprising differences with yeast, while unifying observations that link Paf1C with PolII elongation and RNA processing, and indicate that Paf1C subunits could play roles in controlling transcript length through suppression of PolII accumulation at transcription start site (TSS)-proximal pA sites and regulating pA site choice in 3'UTRs.