Noncanonical transcript forms in yeast and their regulation during environmental stress

A quartet of PIF bHLH factors provides a transcriptionally centered signaling hub that regulates seedling morphogenesis through differential expression-patterning of shared target genes in Arabidopsis

Dark-grown seedlings exhibit skotomorphogenic development. Genetic and molecular evidence indicates that a quartet of Arabidopsis Phytochrome (phy)-Interacting bHLH Factors (PIF1, 3, 4, and 5) are critically necessary to maintaining this developmental state and that light activation of phy induces a switch to photomorphogenic development by inducing rapid degradation of the PIFs. Here, using integrated ChIP–seq and RNA–seq analyses, we have identified genes that are direct targets of PIF3 transcriptional regulation, exerted by sequence-specific binding to G-box (CACGTG) or PBE-box (CACATG) motifs in the target promoters genome-wide. In addition, expression analysis of selected genes in this set, in all triple pif-mutant combinations, provides evidence that the PIF quartet members collaborate to generate an expression pattern that is the product of a mosaic of differential transcriptional responsiveness of individual genes to the different PIFs and of differential regulatory activity of individual PIFs toward the different genes. Together with prior evidence that all four PIFs can bind to G-boxes, the data suggest that this collective activity may be exerted via shared occupancy of binding sites in target promoters.

An important issue in understanding mechanisms of eukaryotic transcriptional regulation is how members of large transcription-factor families, with conserved DNA–binding domains (such as the 162-member Arabidopsis bHLH family), discriminate between target genes. However, the specific question of whether, and to what extent, closely related sub-family members, with potential overlapping functional redundancy (like the quartet of Phytochrome (phy)-Interacting bHLH transcription Factors (PIF1, 3, 4, and 5) studied here), share regulation of target genes through shared binding to promoter-localized consensus motifs does not appear to have been widely investigated. Here, using ChIP–seq analysis, we have identified genes that bind PIF3 to conserved, sequence-specific sites in their promoters; and, using RNA–seq, we have identified those genes displaying altered expression in various pif mutants. Integration of these data identifies those genes that are likely direct targets of transcriptional regulation by PIF3. Our data suggest that the PIF quartet members share directly in transcriptional activation of numerous target genes, potentially via redundant promoter occupancy, in a manner that varies quantitatively from gene to gene. This finding suggests that these PIFs function collectively as a signaling hub, selectively partitioning common upstream signals from light-activated phys at the transcriptional-network interface.

An efficient method for genome-wide polyadenylation site mapping and RNA quantification

The use of alternative poly(A) sites is common and affects the post-transcriptional fate of mRNA, including its stability, subcellular localization and translation. Here, we present a method to identify poly(A) sites in a genome-wide and strand-specific manner. This method, termed 3′T-fill, initially fills in the poly(A) stretch with unlabeled dTTPs, allowing sequencing to start directly after the poly(A) tail into the 3′-untranslated regions (UTR). Our comparative analysis demonstrates that it outperforms existing protocols in quality and throughput and accurately quantifies RNA levels as only one read is produced from each transcript. We use this method to characterize the diversity of polyadenylation in Saccharomyces cerevisiae, showing that alternative RNA molecules are present even in a genetically identical cell population. Finally, we observe that overlap of convergent 3′-UTRs is frequent but sharply limited by coding regions, suggesting factors that restrict compression of the yeast genome.

Decoupling epigenetic and genetic effects through systematic analysis of gene position

Classic "position-effect" experiments repositioned genes near telomeres to demonstrate that the epigenetic landscape can dramatically alter gene expression. Here, we show that systematic gene knockout collections provide an exceptional resource for interrogating position effects, not only near telomeres but at every genetic locus. Because a single reporter gene replaces each deleted gene, interrogating this reporter provides a sensitive probe into different chromatin environments while controlling for genetic context. Using this approach, we find that, whereas systematic replacement of yeast genes with the kanMX marker does not perturb the chromatin landscape, chromatin differences associated with gene position account for 35% of kanMX activity. We observe distinct chromatin influences, including a Set2/Rpd3-mediated antagonistic interaction between histone H3 lysine 36 trimethylation and the Rap1 transcriptional activation site in kanMX. This interaction explains why some yeast genes have been resistant to deletion and allows successful generation of these deletion strains through the use of a modified transformation procedure. These findings demonstrate that chromatin regulation is not governed by a uniform "histone code" but by specific interactions between chromatin and genetic factors.

Analysis of alternative cleavage and polyadenylation by 3' region extraction and deep sequencing

Alternative cleavage and polyadenylation (APA) leads to mRNA isoforms with different coding sequences (CDS) and/or 3′ untranslated regions (3′UTRs). Using 3′ Region Extraction And Deep Sequencing (3′READS), a method which addresses the internal priming and oligo(A) tail issues that commonly plague polyA site (pA) identification, we comprehensively mapped pAs in the mouse genome, thoroughly annotating 3′ ends of genes and revealing over five thousand pAs (~8% of total) flanked by A-rich sequences, which have hitherto been overlooked. About 79% of mRNA genes and 66% of long non-coding RNA (lncRNA) genes have APA; but these two gene types have distinct usage patterns for pAs in introns and upstream exons. Promoter-distal pAs become relatively more abundant during embryonic development and cell differentiation, a trend affecting pAs in both 3′-most exons and upstream regions. Upregulated isoforms generally have stronger pAs, suggesting global modulation of the 3′ end processing activity in development and differentiation.

Genetics and regulatory impact of alternative polyadenylation in human B-lymphoblastoid cells

Gene expression varies widely between individuals of a population, and regulatory change can underlie phenotypes of evolutionary and biomedical relevance. A key question in the field is how DNA sequence variants impact gene expression, with most mechanistic studies to date focused on the effects of genetic change on regulatory regions upstream of protein-coding sequence. By contrast, the role of RNA 3'-end processing in regulatory variation remains largely unknown, owing in part to the challenge of identifying functional elements in 3' untranslated regions. In this work, we conducted a genomic survey of transcript ends in lymphoblastoid cells from genetically distinct human individuals. Our analysis mapped the cis-regulatory architecture of 3' gene ends, finding that transcript end positions did not fall randomly in untranslated regions, but rather preferentially flanked the locations of 3' regulatory elements, including miRNA sites. The usage of these transcript length forms and motifs varied across human individuals, and polymorphisms in polyadenylation signals and other 3' motifs were significant predictors of expression levels of the genes in which they lay. Independent single-gene experiments confirmed the effects of polyadenylation variants on steady-state expression of their respective genes, and validated the regulatory function of 3' cis-regulatory sequence elements that mediated expression of these distinct RNA length forms. Focusing on the immune regulator IRF5, we established the effect of natural variation in RNA 3'-end processing on regulatory response to antigen stimulation. Our results underscore the importance of two mechanisms at play in the genetics of 3'-end variation: the usage of distinct 3'-end processing signals and the effects of 3' sequence elements that determine transcript fate. Our findings suggest that the strategy of integrating observed 3'-end positions with inferred 3' regulatory motifs will prove to be a critical tool in continued efforts to interpret human genome variation.

Incorporating RNA-seq data into the zebrafish Ensembl genebuild

Ensembl gene annotation provides a comprehensive catalog of transcripts aligned to the reference sequence. It relies on publicly available species-specific and orthologous transcripts plus their inferred protein sequence. The accuracy of gene models is improved by increasing the species-specific component that can be cost-effectively achieved using RNA-seq. Two zebrafish gene annotations are presented in Ensembl version 62 built on the Zv9 reference sequence. Firstly, RNA-seq data from five tissues and seven developmental stages were assembled into 25,748 gene models. A 3′-end capture and sequencing protocol was developed to predict the 3′ ends of transcripts, and 46.1% of the original models were subsequently refined. Secondly, a standard Ensembl genebuild, incorporating carefully filtered elements from the RNA-seq-only build, followed by a merge with the manually curated VEGA database, produced a comprehensive annotation of 26,152 genes represented by 51,569 transcripts. The RNA-seq-only and the Ensembl/VEGA genebuilds contribute contrasting elements to the final genebuild. The RNA-seq genebuild was used to adjust intron/exon boundaries of orthologous defined models, confirm their expression, and improve 3′ untranslated regions. Importantly, the inferred protein alignments within the Ensembl genebuild conferred proof of model contiguity for the RNA-seq models. The zebrafish gene annotation has been enhanced by the incorporation of RNA-seq data and the pipeline will be used for other organisms. Organisms with little species-specific cDNA data will generally benefit the most.

Meiosis-induced alterations in transcript architecture and noncoding RNA expression in S. cerevisiae

Changes in transcript architecture can have powerful effects on protein expression. Regulation of the transcriptome is often dramatically revealed during dynamic conditions such as development. To examine changes in transcript architecture we analyzed the expression and transcript boundaries of protein-coding and noncoding RNAs over the developmental process of meiosis in Saccharomyces cerevisiae. Custom-designed, high-resolution tiling arrays were used to define the time-resolved transcriptome of cells undergoing meiosis and sporulation. These arrays were specifically designed for the S. cerevisiae strain SK1 that sporulates with high efficiency and synchrony. In addition, new methods were created to define transcript boundaries and to identify dynamic changes in transcript expression and architecture over time. Of 8407 total segments, 699 (8.3%) were identified by our algorithm as regions containing potential transcript architecture changes. Our analyses reveal extensive changes to both the coding and noncoding transcriptome, including altered 5' ends, 3' ends, and splice sites. Additionally, 3910 (46.5%) unannotated expressed segments were identified. Interestingly, subsets of unannotated RNAs are located across from introns (anti-introns) or across from the junction between two genes (anti-intergenic junctions). Many of these unannotated RNAs are abundant and exhibit sporulation-specific changes in expression patterns. All work, including heat maps of the tiling array, annotation for the SK1 strain, and phastCONS conservation analysis, is available at http://groups.molbiosci.northwestern.edu/sontheimer/sk1meiosis.php. Our high-resolution transcriptome analyses reveal that coding and noncoding transcript architectures are exceptionally dynamic in S. cerevisiae and suggest a vast array of novel transcriptional and post-transcriptional control mechanisms that are activated upon meiosis and sporulation.

ePAT: a simple method to tag adenylated RNA to measure poly(A)-tail length and other 3' RACE applications

The addition of a poly(A)-tail to the 3' termini of RNA molecules influences stability, nuclear export, and efficiency of translation. In the cytoplasm, dynamic changes in the length of the poly(A)-tail have long been recognized as reflective of the switch between translational silence and activation. Thus, measurement of the poly(A)-tail associated with any given mRNA at steady-state can serve as a surrogate readout of its translation-state. Here, we describe a simple new method to 3'-tag adenylated RNA in total RNA samples using the intrinsic property of Escherichia coli DNA polymerase I to extend an RNA primer using a DNA template. This tag can serve as an anchor for cDNA synthesis and subsequent gene-specific PCR to assess poly(A)-tail length. We call this method extension Poly(A) Test (ePAT). The ePAT approach is as efficient as traditional Ligation-Mediated Poly(A) Test (LM-PAT) assays, avoids problems of internal priming associated with oligo-dT-based methods, and allows for the accurate analysis of both the poly(A)-tail length and alternate 3' UTR usage in 3' RACE applications.

A quantitative atlas of polyadenylation in five mammals

We developed PolyA-seq, a strand-specific and quantitative method for high-throughput sequencing of 3′ ends of polyadenylated transcripts, and used it to globally map polyadenylation (polyA) sites in 24 matched tissues in human, rhesus, dog, mouse, and rat. We show that PolyA-seq is as accurate as existing RNA sequencing (RNA-seq) approaches for digital gene expression (DGE), enabling simultaneous mapping of polyA sites and quantitative measurement of their usage. In human, we confirmed 158,533 known sites and discovered 280,857 novel sites (FDR < 2.5%). On average 10% of novel human sites were also detected in matched tissues in other species. Most novel sites represent uncharacterized alternative polyA events and extensions of known transcripts in human and mouse, but primarily delineate novel transcripts in the other three species. A total of 69.1% of known human genes that we detected have multiple polyA sites in their 3′UTRs, with 49.3% having three or more. We also detected polyadenylation of noncoding and antisense transcripts, including constitutive and tissue-specific primary microRNAs. The canonical polyA signal was strongly enriched and positionally conserved in all species. In general, usage of polyA sites is more similar within the same tissues across different species than within a species. These quantitative maps of polyA usage in evolutionarily and functionally related samples constitute a resource for understanding the regulatory mechanisms underlying alternative polyadenylation.

Dynamic changes in nucleosome occupancy are not predictive of gene expression dynamics but are linked to transcription and chromatin regulators

The response to stressful stimuli requires rapid, precise, and dynamic gene expression changes that must be coordinated across the genome. To gain insight into the temporal ordering of genome reorganization, we investigated dynamic relationships between changing nucleosome occupancy, transcription factor binding, and gene expression in Saccharomyces cerevisiae yeast responding to oxidative stress. We applied deep sequencing to nucleosomal DNA at six time points before and after hydrogen peroxide treatment and revealed many distinct dynamic patterns of nucleosome gain and loss. The timing of nucleosome repositioning was not predictive of the dynamics of downstream gene expression change but instead was linked to nucleosome position relative to transcription start sites and specific cis-regulatory elements. We measured genome-wide binding of the stress-activated transcription factor Msn2p over time and found that Msn2p binds different loci with different dynamics. Nucleosome eviction from Msn2p binding sites was common across the genome; however, we show that, contrary to expectation, nucleosome loss occurred after Msn2p binding and in fact required Msn2p. This negates the prevailing model that nucleosomes obscuring Msn2p sites regulate DNA access and must be lost before Msn2p can bind DNA. Together, these results highlight the complexities of stress-dependent chromatin changes and their effects on gene expression.

Genome-wide polyadenylation site mapping

Alternative polyadenylation site usage gives rise to variation in 3' ends of transcripts in diverse organisms ranging from yeast to human. Accurate mapping of polyadenylation sites of transcripts is of major biological importance, since the length of the 3'UTR can have a strong influence on transcript stability, localization, and translation. However, reads generated using total mRNA sequencing mostly lack the very 3' end of transcripts. Here, we present a method that allows simultaneous analysis of alternative 3' ends and transcriptome dynamics at high throughput. By using transcripts produced in vitro, the high precision of end mapping during the protocol can be controlled. This method is illustrated here for budding yeast. However, this method can be applied to any natural or artificially polyadenylated RNA.

Yeast pheromone receptor genes STE2 and STE3 are differently regulated at the transcription and polyadenylation level

The orderly expression of specific genes is the basis for cell differentiation. Saccharomyces cerevisiae has two haploid mating types, a and α cells, in which the mating-specific genes are differentially expressed. When a and α cells are committed to mate, their growth is arrested. Here we show that a cryptic polyadenylation site is present inside the coding region of the a-specific STE2 gene, encoding the receptor for the α-factor. The two cell types produce an incomplete STE2 transcript, but only a cells generate full-length STE2 mRNA. We eliminated the cryptic poly(A) signal, thereby allowing the production of a complete STE2 mRNA in α cells. We mutagenized α cells and isolated a mutant producing full-length STE2 mRNA. The mutation occurred in the ITC1 gene, whose product, together with the product of ISW2, is known to repress STE2 transcriptional initiation. We propose that the regulation of the yeast mating genes is achieved through a concerted mechanism involving transcriptional and posttranscriptional events. In particular, the early poly(A) site in STE2 could contribute to a complete shutoff of its expression in α cells, avoiding autocrine activation and growth arrest. Remarkably, no cryptic poly(A) sites are present in the a-factor receptor STE3 gene, indicating that S. cerevisiae has devised different strategies to regulate the two receptor genes. It is predictable that a correlation between the repression of a gene and the presence of a cryptic poly(A) site could also be found in other organisms, especially when expression of that gene may be harmful.

The evolution of gene expression regulatory networks in yeasts

Gene regulation is a major source of phenotypic diversity between and within species. This aspect of evolution has long been addressed from the sole point of view of the genome sequence. The incredible development of transcriptomics approaches now allows one to actually study the topology and the properties of regulatory networks on an evolutionary perspective. This new discipline is called comparative functional genomics or comparative transcriptomics. This article reviews some of the main advances made in this field, using yeast species, and especially the species sequenced in the frame of the Genolevures program, as a model.

Functional consequences of bidirectional promoters

Several studies have shown that promoters of protein-coding genes are origins of pervasive non-coding RNA transcription and can initiate transcription in both directions. However, only recently have researchers begun to elucidate the functional implications of this bidirectionality and non-coding RNA production. Increasing evidence indicates that non-coding transcription at promoters influences the expression of protein-coding genes, revealing a new layer of transcriptional regulation. This regulation acts at multiple levels, from modifying local chromatin to enabling regional signal spreading and more distal regulation. Moreover, the bidirectional activity of a promoter is regulated at multiple points during transcription, giving rise to diverse types of transcripts.

Comprehensive polyadenylation site maps in yeast and human reveal pervasive alternative polyadenylation

The emerging discoveries on the link between polyadenylation and disease states underline the need to fully characterize genome-wide polyadenylation states. Here, we report comprehensive maps of global polyadenylation events in human and yeast generated using refinements to the Direct RNA Sequencing technology. This direct approach provides a quantitative view of genome-wide polyadenylation states in a strand-specific manner and requires only attomole RNA quantities. The polyadenylation profiles revealed an abundance of unannotated polyadenylation sites, alternative polyadenylation patterns, and regulatory element-associated poly(A)(+) RNAs. We observed differences in sequence composition surrounding canonical and noncanonical human polyadenylation sites, suggesting novel noncoding RNA-specific polyadenylation mechanisms in humans. Furthermore, we observed the correlation level between sense and antisense transcripts to depend on gene expression levels, supporting the view that overlapping transcription from opposite strands may play a regulatory role. Our data provide a comprehensive view of the polyadenylation state and overlapping transcription.