Noncanonical transcript forms in yeast and their regulation during environmental stress

Heat shock induces alternative polyadenylation through dynamic DNA methylation and chromatin looping

Alternative cleavage and polyadenylation (APA) is a gene regulatory mechanism used by cells under stress to upregulate proteostasis-promoting transcripts, but how cells achieve this remains poorly understood. Previously, we elucidated a DNA methylation-regulated APA mechanism, in which gene body DNA methylation enhances distal poly(A) isoform expression by blocking CTCF binding and chromatin loop formation at APA control regions. We hypothesized that DNA methylation-regulated APA is one mechanism cells employ to induce proteostasis-promoting poly(A) isoforms. At the DNAJB6 co-chaperone locus, acute heat shock resulted in binding of stress response transcription factors HSF1, ATF6, and YY1 at the APA control region and an increase in the expression of the proximal poly(A) isoform known to prevent protein aggregation. Furthermore, TET1 was recruited to rapidly demethylate DNA, facilitating CTCF binding and chromatin loop formation, thereby reinforcing preferential proximal poly(A) isoform expression. As cells recovered, the transcription factors vacated the APA control region, and DNMT1 was recruited to remethylate the region. This process resolved chromatin looping and reset the poly(A) isoform expression pattern. Our findings unveil an epigenetic mechanism enabling cells to dynamically modulate poly(A) isoforms in response to stress while shedding light on the interplay between DNA methylation, transcription factor binding, and chromatin looping.

Heat shock induces alternative polyadenylation through dynamic DNA methylation-regulated chromatin looping

Alternative cleavage and polyadenylation (APA) is a gene regulatory mechanism used by cells under stress to upregulate proteostasis-promoting transcripts, but how cells achieve this remains poorly understood. Previously, we elucidated a DNA methylation-regulated APA mechanism, in which gene body DNA methylation enhances distal poly(A) isoform expression by blocking CTCF binding and chromatin loop formation at APA control regions. We hypothesized that DNA methylation-regulated APA is one mechanism cells employ to induce proteostasis-promoting poly(A) isoforms. At the DNAJB6 co-chaperone gene locus, acute heat shock resulted in binding of stress response transcription factors HSF1, ATF6, and YY1 at the APA control region and an increase in the expression of the proximal poly(A) isoform known to prevent protein aggregation. Furthermore, TET1 was recruited to rapidly demethylate DNA, facilitating CTCF binding and chromatin loop formation, thereby reinforcing preferential proximal poly(A) isoform expression. As cells recovered, the transcription factors vacated the APA control region, and DNMT1 was recruited to remethylate the region. This process resolved chromatin looping and reset the poly(A) isoform expression pattern. Our findings unveil an epigenetic mechanism enabling cells to dynamically modulate poly(A) isoforms in response to stress while shedding light on the interplay between DNA methylation, transcription factors, and chromatin looping.

Chromatin Changes in Phytochrome Interacting Factor-Regulated Genes Parallel Their Rapid Transcriptional Response to Light

As sessile organisms, plants must adapt to a changing environment, sensing variations in resource availability and modifying their development in response. Light is one of the most important resources for plants, and its perception by sensory photoreceptors (e.g., phytochromes) and subsequent transduction into long-term transcriptional reprogramming have been well characterized. Chromatin changes have been shown to be involved in photomorphogenesis. However, the initial short-term transcriptional changes produced by light and what factors enable these rapid changes are not well studied. Here, we define rapidly light-responsive, Phytochrome Interacting Factor (PIF) direct-target genes (LRP-DTGs). We found that a majority of these genes also show rapid changes in Histone 3 Lysine-9 acetylation (H3K9ac) in response to the light signal. Detailed time-course analysis of transcript and chromatin changes showed that, for light-repressed genes, H3K9 deacetylation parallels light-triggered transcriptional repression, while for light-induced genes, H3K9 acetylation appeared to somewhat precede light-activated transcript accumulation. However, direct, real-time imaging of transcript elongation in the nucleus revealed that, in fact, transcriptional induction actually parallels H3K9 acetylation. Collectively, the data raise the possibility that light-induced transcriptional and chromatin-remodeling processes are mechanistically intertwined. Histone modifying proteins involved in long term light responses do not seem to have a role in this fast response, indicating that different factors might act at different stages of the light response. This work not only advances our understanding of plant responses to light, but also unveils a system in which rapid chromatin changes in reaction to an external signal can be studied under natural conditions.

Testing the Genomic Shock Hypothesis Using Transposable Element Expression in Yeast Hybrids

Transposable element (TE) insertions are a source of structural variation and can cause genetic instability and gene expression changes. A host can limit the spread of TEs with various repression mechanisms. Many examples of plant and animal interspecific hybrids show disrupted TE repression leading to TE propagation. Recent studies in yeast did not find any increase in transposition rate in hybrids. However, this does not rule out the possibility that the transcriptional or translational activity of TEs increases following hybridization because of a disruption of the host TE control mechanisms. Thus, whether total expression of a TE family is higher in hybrids than in their parental species remains to be examined. We leveraged publically available RNA-seq and ribosomal profiling data on yeast artificial hybrids of the Saccharomyces genus and performed differential expression analysis of their LTR retrotransposons (Ty elements). Our analyses of total mRNA levels show that Ty elements are generally not differentially expressed in hybrids, even when the hybrids are exposed to a low temperature stress condition. Overall, only 2/26 Ty families show significantly higher expression in the S. cerevisiae × S. uvarum hybrids while there are 3/26 showing significantly lower expression in the S. cerevisiae x S. paradoxus hybrids. Our analysis of ribosome profiling data of S. cerevisiae × S. paradoxus hybrids shows similar translation efficiency of Ty in both parents and hybrids, except for Ty1_cer showing higher translation efficiency. Overall, our results do not support the hypothesis that hybridization could act as a systematic trigger of TE expression in yeast and suggest that the impact of hybridization on TE activity is strain and TE specific.

PAS-seq 2: A fast and sensitive method for global profiling of polyadenylated RNAs

Alternative polyadenylation (APA) is a widespread phenomenon in eukaryotes that contributes to regulating gene expression and generating proteomic diversity. APA plays critical roles in development and its mis-regulation has been implicated in a wide variety of human diseases, including cancer. To study APA on the transcriptome-wide level, numerous deep sequencing methods that capture 3' end of mRNAs have been developed in the past decade, but they generally require a large amount of hands-on time and/or high RNA input. Here, we introduce PAS-seq 2, a fast and sensitive method for global and quantitative profiling of polyadenylated RNAs. Compared to our original PAS-seq, this method takes less time and requires much lower total RNA input due to improvement in the reverse transcription process. PAS-seq 2 can be applied to both APA and differential gene expression analyses.

Regulation of the Ysh1 endonuclease of the mRNA cleavage/polyadenylation complex by ubiquitin-mediated degradation

Mutation of the essential yeast protein Ipa1 has previously been demonstrated to cause defects in pre-mRNA 3' end processing and growth, but the mechanism underlying these defects was not clear. In this study, we show that the ipa1-1 mutation causes a striking depletion of Ysh1, the evolutionarily conserved endonuclease subunit of the 19-subunit mRNA Cleavage/Polyadenylation (C/P) complex, but does not decrease other C/P subunits. YSH1 overexpression rescues both the growth and 3' end processing defects of the ipa1-1 mutant. YSH1 mRNA level is unchanged in ipa1-1 cells, and proteasome inactivation prevents Ysh1 loss and causes accumulation of ubiquitinated Ysh1. Ysh1 ubiquitination is mediated by the Ubc4 ubiquitin-conjugating enzyme and Mpe1, which in addition to its function in C/P, is also a RING ubiquitin ligase. In summary, Ipa1 affects mRNA processing by controlling the availability of the C/P endonuclease and may represent a regulatory mechanism that could be rapidly deployed to facilitate reprogramming of cellular responses.

Ipa1 Is an RNA Polymerase II Elongation Factor that Facilitates Termination by Maintaining Levels of the Poly(A) Site Endonuclease Ysh1

The yeast protein Ipa1 was recently discovered to interact with the Ysh1 endonuclease of the prem-RNA cleavage and polyadenylation (C/P) machinery, and Ipa1 mutation impairs 3′end processing. We report that Ipa1 globally promotes proper transcription termination and poly(A) site selection, but with variable effects on genes depending upon the specific configurations of polyadenylation signals. Our findings suggest that the role of Ipa1 in termination is mediated through interaction with Ysh1, since Ipa1 mutation leads to decrease in Ysh1 and poor recruitment of the C/P complex to a transcribed gene. The Ipa1 association with transcriptionally active chromatin resembles that of elongation factors, and the mutant shows defective Pol II elongation kinetics in vivo. Ysh1 overexpression in the Ipa1 mutant rescues the termination defect, but not the mutant’s sensitivity to 6-azauracil, an indicator of defective elongation. Our findings support a model in which an Ipa1/Ysh1 complex helps coordinate transcription elongation and 3′ end processing.

The essential, uncharacterized Ipa1 protein was recently discovered to interact with the Ysh1 endonuclease of the pre-mRNA cleavage and polyadenylation machinery. Pearson et al. propose that the Ipa1/Ysh1 interaction provides the cell with a means to coordinate and regulate transcription elongation with 3′ end processing in accordance with the cell’s needs.

Accurate and Sensitive Quantitation of the Dynamic Heat Shock Proteome Using Tandem Mass Tags

Cells respond to environmental perturbations and insults through modulating protein abundance and function. However, the majority of studies have focused on changes in RNA abundance because quantitative transcriptomics has historically been more facile than quantitative proteomics. Modern Orbitrap mass spectrometers now provide sensitive and deep proteome coverage, allowing direct, global quantification of not only protein abundance but also post-translational modifications (PTMs) that regulate protein activity. We implemented and validated using the well-characterized heat shock response of budding yeast, a tandem mass tagging (TMT), triple-stage mass spectrometry (MS³) strategy to measure global changes in the proteome during the yeast heat shock response over nine time points. We report that basic-pH, ultra-high performance liquid chromatography (UPLC) fractionation of tryptic peptides yields superfractions of minimal redundancy, a crucial requirement for deep coverage and quantification by subsequent LC-MS³. We quantified 2275 proteins across three biological replicates and found that differential expression peaked near 90 min following heat shock (with 868 differentially expressed proteins at 5% false discovery rate). The sensitivity of the approach also allowed us to detect changes in the relative abundance of ubiquitination and phosphorylation PTMs over time. Remarkably, relative quantification of post-translationally modified peptides revealed striking evidence of regulation of the heat shock response by protein PTMs. These data demonstrate that the high precision of TMT-MS³ enables peptide-level quantification of samples, which can reveal important regulation of protein abundance and regulatory PTMs under various experimental conditions.

The polyadenylation factor FIP1 is important for plant development and root responses to abiotic stresses

Root development and its response to environmental changes is crucial for whole plant adaptation. These responses include changes in transcript levels. Here, we show that the alternative polyadenylation (APA) of mRNA is important for root development and responses. Mutations in FIP1, a component of polyadenylation machinery, affects plant development, cell division and elongation, and response to different abiotic stresses. Salt treatment increases the amount of poly(A) site usage within the coding region and 5' untranslated regions (5'-UTRs), and the lack of FIP1 activity reduces the poly(A) site usage within these non-canonical sites. Gene ontology analyses of transcripts displaying APA in response to salt show an enrichment in ABA signaling, and in the response to stresses such as salt or cadmium (Cd), among others. Root growth assays show that fip1-2 is more tolerant to salt but is hypersensitive to ABA or Cd. Our data indicate that FIP1-mediated alternative polyadenylation is important for plant development and stress responses.

Alternative cleavage and polyadenylation in health and disease

Most human genes have multiple sites at which RNA 3' end cleavage and polyadenylation can occur, enabling the expression of distinct transcript isoforms under different conditions. Novel methods to sequence RNA 3' ends have generated comprehensive catalogues of polyadenylation (poly(A)) sites; their analysis using innovative computational methods has revealed how poly(A) site choice is regulated by core RNA 3' end processing factors, such as cleavage factor I and cleavage and polyadenylation specificity factor, as well as by other RNA-binding proteins, particularly splicing factors. Here, we review the experimental and computational methods that have enabled the global mapping of mRNA and of long non-coding RNA 3' ends, quantification of the resulting isoforms and the discovery of regulators of alternative cleavage and polyadenylation (APA). We highlight the different types of APA-derived isoforms and their functional differences, and illustrate how APA contributes to human diseases, including cancer and haematological, immunological and neurological diseases.

Polyadenylation sites and their characteristics in the genome of channel catfish (Ictalurus punctatus) as revealed by using RNA-Seq data

Polyadenylation plays important roles in gene expression regulation in eukaryotes, which typically involves cleavage and poly(A) tail addition at the polyadenylation site (PAS) of the pre-mature mRNA. Many eukaryotic genes contain more than one PASs, termed as alternative polyadenylation (APA). As a crucial post-transcriptional regulation, polyadenylation affects various aspects of RNA metabolism such as mRNA stability, translocation, and translation. However, polyadenylation has been rarely studied in teleosts. Here we conducted polyadenylation analysis in channel catfish, a commercially important aquaculture species around the world. Using RNA-Seq data, we identified 20,320 PASs which were classified into 14,500 clusters by merging adjacent PASs. Most of the PASs were found in 3' UTRs, followed by intron regions based on the annotation of channel catfish reference genome. No apparent difference in PAS distribution was observed between the sense and antisense strand of the channel catfish genome. The sequence analysis of nucleotide composition and motif around PASs yielded a highly similar profile among various organisms, suggesting the conservation and importance of polyadenylation in evolution. Using APA genes with more than two PASs, gene ontology enrichment revealed genes particularly involved in RNA binding. Reactome pathway analysis showed the enrichment of the innate immune system, especially neutrophil degranulation.

Functional genomics for the oleaginous yeast Yarrowia lipolytica

Oleaginous yeasts are valuable systems for biosustainable production of hydrocarbon-based chemicals. Yarrowia lipolytica is one of the best characterized of these yeast with respect to genome annotation and flux analysis of metabolic processes. Nonetheless, progress is hampered by a dearth of genome-wide tools enabling functional genomics. In order to remedy this deficiency, we developed a library of Y. lipolytica insertion mutants via transposon mutagenesis. The Hermes DNA transposon was expressed to achieve saturation mutagenesis of the genome. Over 534,000 independent insertions were identified by next-generation sequencing. Poisson analysis of insertion density classified ~ 22% of genes as essential. As expected, most essential genes have homologs in Saccharomyces cerevisiae and Schizosaccharomyces pombe, and the majority of those are also essential. As an obligate aerobe, Y. lipolytica has significantly more respiration - related genes that are classified as essential than do S. cerevisiae and S. pombe. Contributions of non-essential genes to growth in glucose and glycerol carbon sources were assessed and used to evaluate two recent genome-scale models of Y. lipolytica metabolism. Fluorescence-activated cell sorting identified mutants in which lipid accumulation is increased. Our findings provide insights into biosynthetic pathways, compartmentalization of enzymes, and distinct functions of paralogs. This functional genomic analysis of the oleaginous yeast Y. lipolytica provides an important resource for modeling, bioengineering, and design of synthetic minimalized strains of respiratory yeasts.

Cellular stress alters 3'UTR landscape through alternative polyadenylation and isoform-specific degradation

Most eukaryotic genes express alternative polyadenylation (APA) isoforms with different 3'UTR lengths, production of which is influenced by cellular conditions. Here, we show that arsenic stress elicits global shortening of 3'UTRs through preferential usage of proximal polyadenylation sites during stress and enhanced degradation of long 3'UTR isoforms during recovery. We demonstrate that RNA-binding protein TIA1 preferentially interacts with alternative 3'UTR sequences through U-rich motifs, correlating with stress granule association and mRNA decay of long 3'UTR isoforms. By contrast, genes with shortened 3'UTRs due to stress-induced APA can evade mRNA clearance and maintain transcript abundance post stress. Furthermore, we show that stress causes distinct 3'UTR size changes in proliferating and differentiated cells, highlighting its context-specific impacts on the 3'UTR landscape. Together, our data reveal a global, 3'UTR-based mRNA stability control in stressed cells and indicate that APA can function as an adaptive mechanism to preserve mRNAs in response to stress.

Conserved non-AUG uORFs revealed by a novel regression analysis of ribosome profiling data

Upstream open reading frames (uORFs), located in transcript leaders (5' UTRs), are potent cis-acting regulators of translation and mRNA turnover. Recent genome-wide ribosome profiling studies suggest that thousands of uORFs initiate with non-AUG start codons. Although intriguing, these non-AUG uORF predictions have been made without statistical control or validation; thus, the importance of these elements remains to be demonstrated. To address this, we took a comparative genomics approach to study AUG and non-AUG uORFs. We mapped transcription leaders in multiple Saccharomyces yeast species and applied a novel machine learning algorithm (uORF-seqr) to ribosome profiling data to identify statistically significant uORFs. We found that AUG and non-AUG uORFs are both frequently found in Saccharomyces yeasts. Although most non-AUG uORFs are found in only one species, hundreds have either conserved sequence or position within Saccharomyces uORFs initiating with UUG are particularly common and are shared between species at rates similar to that of AUG uORFs. However, non-AUG uORFs are translated less efficiently than AUG-uORFs and are less subject to removal via alternative transcription initiation under normal growth conditions. These results suggest that a subset of non-AUG uORFs may play important roles in regulating gene expression.

mRNA detection in budding yeast with single fluorophores

Quantitative measurement of mRNA levels in single cells is necessary to understand phenotypic variability within an otherwise isogenic population of cells. Single-molecule mRNA Fluorescence In Situ Hybridization (FISH) has been established as the standard method for this purpose, but current protocols require a long region of mRNA to be targeted by multiple DNA probes. Here, we introduce a new single-probe FISH protocol termed sFISH for budding yeast, Saccharomyces cerevisiae using a single DNA probe labeled with a single fluorophore. In sFISH, we markedly improved probe specificity and signal-to-background ratio by using methanol fixation and inclined laser illumination. We show that sFISH reports mRNA changes that correspond to protein levels and gene copy number. Using this new FISH protocol, we can detect >50% of the total target mRNA. We also demonstrate the versatility of sFISH using FRET detection and mRNA isoform profiling as examples. Our FISH protocol with single-fluorophore sensitivity significantly reduces cost and time compared to the conventional FISH protocols and opens up new opportunities to investigate small changes in RNA at the single cell level.

A Genetic Tool to Track Protein Aggregates and Control Prion Inheritance

Protein aggregation is a hallmark of many diseases but also underlies a wide range of positive cellular functions. This phenomenon has been difficult to study because of a lack of quantitative and high-throughput cellular tools. Here, we develop a synthetic genetic tool to sense and control protein aggregation. We apply the technology to yeast prions, developing sensors to track their aggregation states and employing prion fusions to encode synthetic memories in yeast cells. Utilizing high-throughput screens, we identify prion-curing mutants and engineer "anti-prion drives" that reverse the non-Mendelian inheritance pattern of prions and eliminate them from yeast populations. We extend our technology to yeast RNA-binding proteins (RBPs) by tracking their propensity to aggregate, searching for co-occurring aggregates, and uncovering a group of coalescing RBPs through screens enabled by our platform. Our work establishes a quantitative, high-throughput, and generalizable technology to study and control diverse protein aggregation processes in cells.

From transcriptional complexity to cellular phenotypes: Lessons from yeast

Pervasive transcription has been reported in many eukaryotic organisms, revealing a highly interleaved transcriptome organization that involves thousands of coding and non-coding RNAs. However, to date, the biological impact of transcriptome complexity is still poorly understood. Here I will review how subtle variations of the transcriptome can lead to divergent cellular phenotypes by fine-tuning both its coding potential and regulation. I will discuss strategies that can be used to link molecular variations with divergent biological outcomes. Finally, I will explore the implication of transcriptional complexity for our understanding of gene expression in the context of cell-to-cell phenotypic variability. Copyright © 2017 John Wiley & Sons, Ltd.

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.

Acidophilic green algal genome provides insights into adaptation to an acidic environment

Some microalgae are adapted to extremely acidic environments in which toxic metals are present at high levels. However, little is known about how acidophilic algae evolved from their respective neutrophilic ancestors by adapting to particular acidic environments. To gain insights into this issue, we determined the draft genome sequence of the acidophilic green alga Chlamydomonas eustigma and performed comparative genome and transcriptome analyses between Ceustigma and its neutrophilic relative Chlamydomonas reinhardtii The results revealed the following features in Ceustigma that probably contributed to the adaptation to an acidic environment. Genes encoding heat-shock proteins and plasma membrane H⁺-ATPase are highly expressed in Ceustigma This species has also lost fermentation pathways that acidify the cytosol and has acquired an energy shuttle and buffering system and arsenic detoxification genes through horizontal gene transfer. Moreover, the arsenic detoxification genes have been multiplied in the genome. These features have also been found in other acidophilic green and red algae, suggesting the existence of common mechanisms in the adaptation to acidic environments.

Alternative polyadenylation in the regulation and dysregulation of gene expression

Transcriptional control shapes a cell's transcriptome composition, but it is RNA processing that refines its expression. The untranslated regions (UTRs) of mRNA are hotspots for regulatory control. Features in these can impact mRNA stability, localisation and translation. Here we describe how alternative cleavage and polyadenylation can change mRNA fate by changing the length of its 3'UTR.

Cleavage and polyadenylation: Ending the message expands gene regulation

Cleavage and polyadenylation (pA) is a fundamental step that is required for the maturation of primary protein encoding transcripts into functional mRNAs that can be exported from the nucleus and translated in the cytoplasm. 3'end processing is dependent on the assembly of a multiprotein processing complex on the pA signals that reside in the pre-mRNAs. Most eukaryotic genes have multiple pA signals, resulting in alternative cleavage and polyadenylation (APA), a widespread phenomenon that is important to establish cell state and cell type specific transcriptomes. Here, we review how pA sites are recognized and comprehensively summarize how APA is regulated and creates mRNA isoform profiles that are characteristic for cell types, tissues, cellular states and disease.

Noncanonical Alternative Polyadenylation Contributes to Gene Regulation in Response to Hypoxia

Stresses from various environmental challenges continually confront plants, and their responses are important for growth and survival. One molecular response to such challenges involves the alternative polyadenylation of mRNA. In plants, it is unclear how stress affects the production and fate of alternative mRNA isoforms. Using a genome-scale approach, we show that in Arabidopsis thaliana, hypoxia leads to increases in the number of mRNA isoforms with polyadenylated 3' ends that map to 5'-untranslated regions (UTRs), introns, and protein-coding regions. RNAs with 3' ends within protein-coding regions and introns were less stable than mRNAs that end at 3'-UTR poly(A) sites. Additionally, these RNA isoforms were underrepresented in polysomes isolated from control and hypoxic plants. By contrast, mRNA isoforms with 3' ends that lie within annotated 5'-UTRs were overrepresented in polysomes and were as stable as canonical mRNA isoforms. These results indicate that the generation of noncanonical mRNA isoforms is an important feature of the abiotic stress response. The finding that several noncanonical mRNA isoforms are relatively unstable suggests that the production of non-stop and intronic mRNA isoforms may represent a form of negative regulation in plants, providing a conceptual link with mechanisms that generate these isoforms (such as alternative polyadenylation) and RNA surveillance.

Genome-Wide Identification of Alternative Polyadenylation Events Using 3'T-Fill

Due to the increasing appreciation of the impact of alternative polyadenylation on cellular biology, our straightforward, scalable method is of interest to any researcher studying eukaryotic transcription. In addition to high quality gene expression measurements, it precisely maps poly(A) sites and thereby permits the distinction between differential 3'UTR isoforms. As sequencing through long homopolymer stretches is not possible on the Illumina platform, we developed a method that fills up the poly(A) stretch with dTTPs before the sequencing reaction starts.

The Antisense Transcriptome and the Human Brain

The transcriptome of a cell is made up of a varied array of RNA species, including protein-coding RNAs, long non-coding RNAs, short non-coding RNAs, and circular RNAs. The cellular transcriptome is dynamic and can change depending on environmental factors, disease state and cellular context. The human brain has perhaps the most diverse transcriptome profile that is enriched for many species of RNA, including antisense transcripts. Antisense transcripts are produced when both the plus and minus strand of the DNA helix are transcribed at a particular locus. This results in an RNA transcript that has a partial or complete overlap with an intronic or exonic region of the sense transcript. While antisense transcription is known to occur at some level in most organisms, this review focuses specifically on antisense transcription in the brain and how regulation of genes by antisense transcripts can contribute to functional aspects of the healthy and diseased brain. First, we discuss different techniques that can be used in the identification and quantification of antisense transcripts. This is followed by examples of antisense transcription and modes of regulatory function that have been identified in the brain.

PAT-seq: a method to study the integration of 3'-UTR dynamics with gene expression in the eukaryotic transcriptome

A major objective of systems biology is to quantitatively integrate multiple parameters from genome-wide measurements. To integrate gene expression with dynamics in poly(A) tail length and adenylation site, we developed a targeted next-generation sequencing approach, Poly(A)-Test RNA-sequencing. PAT-seq returns (i) digital gene expression, (ii) polyadenylation site/s, and (iii) the polyadenylation-state within and between eukaryotic transcriptomes. PAT-seq differs from previous 3' focused RNA-seq methods in that it depends strictly on 3' adenylation within total RNA samples and that the full-native poly(A) tail is included in the sequencing libraries. Here, total RNA samples from budding yeast cells were analyzed to identify the intersect between adenylation state and gene expression in response to loss of the major cytoplasmic deadenylase Ccr4. Furthermore, concordant changes to gene expression and adenylation-state were demonstrated in the classic Crabtree-Warburg metabolic shift. Because all polyadenylated RNA is interrogated by the approach, alternative adenylation sites, noncoding RNA and RNA-decay intermediates were also identified. Most important, the PAT-seq approach uses standard sequencing procedures, supports significant multiplexing, and thus replication and rigorous statistical analyses can for the first time be brought to the measure of 3'-UTR dynamics genome wide.

Mapping 3' mRNA isoforms on a genomic scale

Most eukaryotic genes are transcribed into mRNAs with alternative poly(A) sites. Emerging evidence suggests that mRNA isoforms with alternative poly(A) sites can perform critical regulatory functions in numerous biological processes. In recent years, a number of strategies utilizing high-throughput sequencing technologies have been developed to aid in the identification of genome-wide poly(A) sites. This unit describes a modified protocol for a recently published 3'READS (3' region extraction and deep sequencing) method that accurately identifies genome-wide poly(A) sites and that can be used to quantify the relative abundance of the resulting 3' mRNA isoforms. This approach minimizes nonspecific sequence reads due to internal priming and typically yields a high percentage of sequence reads that are ideally suited for accurate poly(A) identification.

Sequencing of first-strand cDNA library reveals full-length transcriptomes

Massively parallel strand-specific sequencing of RNA (ssRNA-seq) has emerged as a powerful tool for profiling complex transcriptomes. However, many current methods for ssRNA-seq suffer from the underrepresentation of both the 5' and 3' ends of RNAs, which can be attributed to second-strand cDNA synthesis. The 5' and 3' ends of RNA harbour crucial information for gene regulation; namely, transcription start sites (TSSs) and polyadenylation sites. Here we report a novel ssRNA-seq method that does not involve second-strand cDNA synthesis, as we Directly Ligate sequencing Adaptors to the First-strand cDNA (DLAF). This novel method with fewer enzymatic reactions results in a higher quality of the libraries than the conventional method. Sequencing of DLAF libraries followed by a novel analysis pipeline enables the profiling of both 5' ends and polyadenylation sites at near-base resolution. Therefore, DLAF offers the first genomics tool to obtain the 'full-length' transcriptome with a single library.

Alternative polyadenylation can regulate post-translational membrane localization

For many genomic loci, there are more than one potential cleavage and polyadenylation site, resulting in the generation of multiple distinct transcripts. When the proximal polyadenylation site is present within the coding region of the transcript, alternative polyadenylation can result in proteins with distinct amino acid sequences and potentially distinct functions. In most cases, the different possible polyadenylation sites are all present within the 3' untranslated regions (UTRs), and the amino acid sequence of the encoded proteins are not affected by polyadenylation site selection. In individual instances, the selection of the proximal versus distal polyadenylation site in the 3'UTR can dramatically affect transcript stability and translatability. In some instances, UTR alternative polyadenylation generates RNA isoforms that have distinct subcellular localization patterns, and that can regulate the location of the encoded protein in an RNA-guided manner. In a recent paper, the laboratory of Christine Mayr demonstrated that alternative polyadenylation of the transmembrane protein CD47 results in transcripts with the same localization pattern, but the encoded protein localizes to the endoplasmic reticulum when it is encoded by the transcript generated by using the proximal polyadenylation site in 3'UTR, and the identical protein localizes to the plasma membrane when the transcript is encoded by the distal polyadenylation site, also in the 3' UTR. Unlike previous studies, the mechanism of localization does not rely on differential trafficking of the mRNA and is instead, based on RNA-mediated recruitment of proteins to the cytoplasmic side of CD47 that support its plasma membrane localization. Other transmembrane proteins were discovered to be regulated similarly. The results demonstrate that the choice of polyadenylation site can affect protein localization and function, even when the sequence of the protein is unaffected. Further, the transcript encoding a protein can serve as a scaffold to recruit additional proteins that affect the protein's fate.

Using Klenow-mediated extension to measure poly(A)-tail length and position in the transcriptome

The poly(A)-tail that terminates most mRNA and many noncoding RNA is a convenient "hook" to isolate mRNA. However the length of this tail and its position within the primary RNA transcript can also hold diagnostic value for RNA metabolism. In general, mRNA with a long poly(A)-tail is well translated, whereas a short poly(A)-tail can indicate translational silencing. A short poly(A)-tail is also appended to RNA-decay intermediates via the TRAMP complex. A number of approaches have been developed to measure the length and position of the poly(A)-tail. Here, we describe a simple method to "tag" adenylated RNA using the native function of DNA polymerase I to extend an RNA primer on a DNA template in second-strand DNA synthesis. This function can be harnessed as a means to purify, visualize, and quantitate poly(A)-dynamics of individual RNA and the transcriptome en masse.

Methods for processing high-throughput RNA sequencing data

High-throughput sequencing (HTS) methods for analyzing RNA populations (RNA-Seq) are gaining rapid application to many experimental situations. The steps in an RNA-Seq experiment require thought and planning, especially because the expense in time and materials is currently higher and the protocols are far less routine than those used for other high-throughput methods, such as microarrays. As always, good experimental design will make analysis and interpretation easier. Having a clear biological question, an idea about the best way to do the experiment, and an understanding of the number of replicates needed will make the entire process more satisfying. Whether the goal is capturing transcriptome complexity from a tissue or identifying small fragments of RNA cross-linked to a protein of interest, conversion of the RNA to cDNA followed by direct sequencing using the latest methods is a developing practice, with new technical modifications and applications appearing every day. Even more rapid are the development and improvement of methods for analysis of the very large amounts of data that arrive at the end of an RNA-Seq experiment, making considerations regarding reproducibility, validation, visualization, and interpretation increasingly important. This introduction is designed to review and emphasize a pathway of analysis from experimental design through data presentation that is likely to be successful, with the recognition that better methods are right around the corner.

Transcriptional interference by antisense RNA is required for circadian clock function

Eukaryotic circadian oscillators consist of negative feedback loops that generate endogenous rhythmicities¹. Natural antisense RNAs are found in a wide range of eukaryotic organisms^2-5. Nevertheless, the physiological importance and mode of action of most antisense RNAs is not clear^6-9. frequency (frq) encodes a component of the Neurospora core circadian negative feedback loop which was thought to generate sustained rhythmicity¹⁰. Transcription of qrf, the long non-coding frq antisense RNA, is light induced, and its level oscillates in antiphase to frq sense RNA³. Here we show that qrf transcription is regulated by both light-dependent and -independent mechanisms. Light-dependent qrf transcription represses frq expression and regulates clock resetting. qrf expression in the dark, on the other hand, is required for circadian rhythmicity. frq transcription also inhibits qrf expression and surprisingly, drives the antiphasic rhythm of qrf transcripts. The mutual inhibition of frq and qrf transcription thus forms a double negative feedback loop that is interlocked with the core feedback loop. Genetic and mathematical modeling analyses indicate that such an arrangement is required for robust and sustained circadian rhythmicity. Moreover, our results suggest that antisense transcription inhibits sense expression by mediating chromatin modifications and premature transcription termination. Together, our results established antisense transcription as an essential feature in a circadian system and shed light on the importance and mechanism of antisense action.

Dom34 rescues ribosomes in 3' untranslated regions

Ribosomes that stall before completing peptide synthesis must be recycled and returned to the cytoplasmic pool. The protein Dom34 and cofactors Hbs1 and Rli1 can dissociate stalled ribosomes in vitro, but the identity of targets in the cell is unknown. Here, we extend ribosome profiling methodology to reveal a high-resolution molecular characterization of Dom34 function in vivo. Dom34 removes stalled ribosomes from truncated mRNAs, but, in contrast, does not generally dissociate ribosomes on coding sequences known to trigger stalling, such as polyproline. We also show that Dom34 targets arrested ribosomes near the ends of 3' UTRs. These ribosomes appear to gain access to the 3' UTR via a mechanism that does not require decoding of the mRNA. These results suggest that ribosomes frequently enter downstream noncoding regions and that Dom34 carries out the important task of rescuing them.

Altered gamma-aminobutyric acid type B receptor subunit 1 splicing in alcoholics

BACKGROUND

Chronic alcohol exposure can change splice variant expression. The gamma-aminobutyric acid type B (GABAB) receptor undergoes splicing and is an alcoholism treatment target, but there is little information about splicing changes in this receptor in alcoholics. We studied GABAB receptor subunit 1 (GABAB1) splicing in alcoholic postmortem brains.

METHODS

To maximize GABAB1 splice junction identification, we combined gene specific libraries with RNA-seq. Splice junctions and mapped reads were also found from intronic and intergenic regions. We compared GABAB1 splice junctions in prefrontal cortices from 14 alcoholic and 15 control subjects and introduced new strategies, reads per kilobase of splice junction model per million mapped reads and reads per kilobase of gene model per million mapped reads, for quantitating splice junction and gene expression.

RESULTS

Novel splice junction detection indicated that the GABAB1 gene is at least two times longer than the previously reported gene length. GABAB1 exon and intron expression data showed low expression at the 5' end exons and exon grouping. This indicated that there are short splicing variants in addition to GABAB receptor subunit GABAB1a, the longest known major transcript. We found that chronic alcohol altered exon/intron expression and splice junction levels. Decreased expression of the gamma-aminobutyric acid binding site, a transmembrane domain and a microRNA binding site may decrease normal GABAB1 transcript population and thereby decrease normal signal transduction in alcoholics.

CONCLUSIONS

We discovered novel, complex splicing of GABAB1 in human brain and showed that chronic alcohol produces additional splicing complexity.

Reconstructing and reprogramming the tumor-propagating potential of glioblastoma stem-like cells

Developmental fate decisions are dictated by master transcription factors (TFs) that interact with cis-regulatory elements to direct transcriptional programs. Certain malignant tumors may also depend on cellular hierarchies reminiscent of normal development but superimposed on underlying genetic aberrations. In glioblastoma (GBM), a subset of stem-like tumor-propagating cells (TPCs) appears to drive tumor progression and underlie therapeutic resistance yet remain poorly understood. Here, we identify a core set of neurodevelopmental TFs (POU3F2, SOX2, SALL2, and OLIG2) essential for GBM propagation. These TFs coordinately bind and activate TPC-specific regulatory elements and are sufficient to fully reprogram differentiated GBM cells to "induced" TPCs, recapitulating the epigenetic landscape and phenotype of native TPCs. We reconstruct a network model that highlights critical interactions and identifies candidate therapeutic targets for eliminating TPCs. Our study establishes the epigenetic basis of a developmental hierarchy in GBM, provides detailed insight into underlying gene regulatory programs, and suggests attendant therapeutic strategies. PAPERCLIP:

TAIL-seq: genome-wide determination of poly(A) tail length and 3' end modifications

Global investigation of the 3' extremity of mRNA (3'-terminome), despite its importance in gene regulation, has not been feasible due to technical challenges associated with homopolymeric sequences and relative paucity of mRNA. We here develop a method, TAIL-seq, to sequence the very end of mRNA molecules. TAIL-seq allows us to measure poly(A) tail length at the genomic scale. Median poly(A) length is 50-100 nt in HeLa and NIH 3T3 cells. Poly(A) length correlates with mRNA half-life, but not with translational efficiency. Surprisingly, we discover widespread uridylation and guanylation at the downstream of poly(A) tail. The U tails are generally attached to short poly(A) tails (<25 nt), while the G tails are found mainly on longer poly(A) tails (>40 nt), implicating their generic roles in mRNA stability control. TAIL-seq is a potent tool to dissect dynamic control of mRNA turnover and translational control, and to discover unforeseen features of RNA cleavage and tailing.

Alternative polyadenylation diversifies post-transcriptional regulation by selective RNA-protein interactions

Recent research has uncovered extensive variability in the boundaries of transcript isoforms, yet the functional consequences of this variation remain largely unexplored. Here, we systematically discriminate between the molecular phenotypes of overlapping coding and non‐coding transcriptional events from each genic locus using a novel genome‐wide, nucleotide‐resolution technique to quantify the half‐lives of 3′ transcript isoforms in yeast. Our results reveal widespread differences in stability among isoforms for hundreds of genes in a single condition, and that variation of even a single nucleotide in the 3′ untranslated region (UTR) can affect transcript stability. While previous instances of negative associations between 3′ UTR length and transcript stability have been reported, here, we find that shorter isoforms are not necessarily more stable. We demonstrate the role of RNA‐protein interactions in conditioning isoform‐specific stability, showing that PUF3 binds and destabilizes specific polyadenylation isoforms. Our findings indicate that although the functional elements of a gene are encoded in DNA sequence, the selective incorporation of these elements into RNA through transcript boundary variation allows a single gene to have diverse functional consequences.

Strand-Specific RNA-Seq Provides Greater Resolution of Transcriptome Profiling

RNA-Seq is a recently developed sequencing technology, that through the analysis of cDNA allows for unique insights into the transcriptome of a cell. The data generated by RNA-Seq provides information on gene expression, alternative splicing events and the presence of non-coding RNAs. It has been realised non-coding RNAs are more then just artefacts of erroneous transcription and play vital regulatory roles at the genomic, transcriptional and translational level. Transcription of the DNA sense strand produces antisense transcripts. This is known as antisense transcription and often results in the production of non-coding RNAs that are complementary to their associated sense transcripts. Antisense tran-scription has been identified in bacteria, fungi, protozoa, plants, invertebrates and mammals. It seems that antisense tran-scriptional ‘hot spots’ are located around nucleosome-free regions such as those associated with promoters, indicating that it is likely that antisense transcripts carry out important regulatory functions. This underlines the importance of identifying the presence and understanding the function of these antisense non-coding RNAs. The information concerning strand ori-gin is often lost during conventional RNA-Seq; capturing this information would substantially increase the worth of any RNA-Seq experiment. By manipulating the input cDNA during the template preparation stage it is possible to retain this vital information. This forms the basis of strand-specific RNA-Seq. With an ability to unlock immense portions of new in-formation surrounding the transcriptome, this cutting edge technology may hold the key to developing a greater under-standing of the transcriptome.

Alternative cleavage and polyadenylation: extent, regulation and function

The 3' end of most protein-coding genes and long non-coding RNAs is cleaved and polyadenylated. Recent discoveries have revealed that a large proportion of these genes contains more than one polyadenylation site. Therefore, alternative polyadenylation (APA) is a widespread phenomenon, generating mRNAs with alternative 3' ends. APA contributes to the complexity of the transcriptome by generating isoforms that differ either in their coding sequence or in their 3' untranslated regions (UTRs), thereby potentially regulating the function, stability, localization and translation efficiency of target RNAs. Here, we review our current understanding of the polyadenylation process and the latest progress in the identification of APA events, mechanisms that regulate poly(A) site selection, and biological processes and diseases resulting from APA.

Inferring evolutionary histories of pathway regulation from transcriptional profiling data

One of the outstanding challenges in comparative genomics is to interpret the evolutionary importance of regulatory variation between species. Rigorous molecular evolution-based methods to infer evidence for natural selection from expression data are at a premium in the field, and to date, phylogenetic approaches have not been well-suited to address the question in the small sets of taxa profiled in standard surveys of gene expression. We have developed a strategy to infer evolutionary histories from expression profiles by analyzing suites of genes of common function. In a manner conceptually similar to molecular evolution models in which the evolutionary rates of DNA sequence at multiple loci follow a gamma distribution, we modeled expression of the genes of an a priori-defined pathway with rates drawn from an inverse gamma distribution. We then developed a fitting strategy to infer the parameters of this distribution from expression measurements, and to identify gene groups whose expression patterns were consistent with evolutionary constraint or rapid evolution in particular species. Simulations confirmed the power and accuracy of our inference method. As an experimental testbed for our approach, we generated and analyzed transcriptional profiles of four Saccharomyces yeasts. The results revealed pathways with signatures of constrained and accelerated regulatory evolution in individual yeasts and across the phylogeny, highlighting the prevalence of pathway-level expression change during the divergence of yeast species. We anticipate that our pathway-based phylogenetic approach will be of broad utility in the search to understand the evolutionary relevance of regulatory change.

Comparative transcriptomic studies routinely identify thousands of genes differentially expressed between species. The central question in the field is whether and how such regulatory changes have been the product of natural selection. Can the signal of evolutionarily relevant expression divergence be detected amid the noise of changes resulting from genetic drift? Our work develops a theory of gene expression variation among a suite of genes that function together. We derive a formalism that relates empirical observations of expression of pathway genes in divergent species to the underlying strength of natural selection on expression output. We show that fitting this type of model to simulated data accurately recapitulates the parameters used to generate the simulation. We then make experimental measurements of gene expression in a panel of single-celled eukaryotic yeast species. To these data we apply our inference method, and identify pathways with striking evidence for accelerated or constrained regulatory evolution, in particular species and across the phylogeny. Our method provides a key advance over previous approaches in that it maximizes the power of rigorous molecular-evolution analysis of regulatory variation even when data are relatively sparse. As such, the theory and tools we have developed will likely find broad application in the field of comparative genomics.

Locus-specific editing of histone modifications at endogenous enhancers

Mammalian gene regulation is dependent on tissue-specific enhancers that can act across large distances to influence transcriptional activity. Mapping experiments have identified hundreds of thousands of putative enhancers whose functionality is supported by cell type-specific chromatin signatures and striking enrichments for disease-associated sequence variants. However, these studies did not address the in vivo functions of the putative elements or their chromatin states and did not determine which genes, if any, a given enhancer regulates. Here we present a strategy to investigate endogenous regulatory elements by selectively altering their chromatin state using programmable reagents. Transcription activator-like (TAL) effector repeat domains fused to the LSD1 histone demethylase efficiently remove enhancer-associated chromatin modifications from target loci, without affecting control regions. We find that inactivation of enhancer chromatin by these fusion proteins frequently causes downregulation of proximal genes, revealing enhancer target genes. Our study demonstrates the potential of epigenome editing tools to characterize an important class of functional genomic elements.

The genomic and transcriptomic landscape of a HeLa cell line

HeLa is the most widely used model cell line for studying human cellular and molecular biology. To date, no genomic reference for this cell line has been released, and experiments have relied on the human reference genome. Effective design and interpretation of molecular genetic studies performed using HeLa cells require accurate genomic information. Here we present a detailed genomic and transcriptomic characterization of a HeLa cell line. We performed DNA and RNA sequencing of a HeLa Kyoto cell line and analyzed its mutational portfolio and gene expression profile. Segmentation of the genome according to copy number revealed a remarkably high level of aneuploidy and numerous large structural variants at unprecedented resolution. Some of the extensive genomic rearrangements are indicative of catastrophic chromosome shattering, known as chromothripsis. Our analysis of the HeLa gene expression profile revealed that several pathways, including cell cycle and DNA repair, exhibit significantly different expression patterns from those in normal human tissues. Our results provide the first detailed account of genomic variants in the HeLa genome, yielding insight into their impact on gene expression and cellular function as well as their origins. This study underscores the importance of accounting for the strikingly aberrant characteristics of HeLa cells when designing and interpreting experiments, and has implications for the use of HeLa as a model of human biology.

Post-transcriptional regulation of iron homeostasis in Saccharomyces cerevisiae

Iron is an essential micronutrient for all eukaryotic organisms because it participates as a redox cofactor in a wide variety of biological processes. Recent studies in Saccharomyces cerevisiae have shown that in response to iron deficiency, an RNA-binding protein denoted Cth2 coordinates a global metabolic rearrangement that aims to optimize iron utilization. The Cth2 protein contains two Cx₈Cx₅Cx₃H tandem zinc fingers (TZFs) that specifically bind to adenosine/uridine-rich elements within the 3′ untranslated region of many mRNAs to promote their degradation. The Cth2 protein shuttles between the nucleus and the cytoplasm. Once inside the nucleus, Cth2 binds target mRNAs and stimulates alternative 3′ end processing. A Cth2/mRNA-containing complex is required for export to the cytoplasm, where the mRNA is degraded by the 5′ to 3′ degradation pathway. This post-transcriptional regulatory mechanism limits iron utilization in nonessential pathways and activates essential iron-dependent enzymes such as ribonucleotide reductase, which is required for DNA synthesis and repair. Recent findings indicate that the TZF-containing tristetraprolin protein also functions in modulating human iron homeostasis. Elevated iron concentrations can also be detrimental for cells. The Rnt1 RNase III exonuclease protects cells from excess iron by promoting the degradation of a subset of the Fe acquisition system when iron levels rise.

RNaseIII and T4 polynucleotide Kinase sequence biases and solutions during RNA-seq library construction

Background

RNA-seq is a next generation sequencing method with a wide range of applications including single nucleotide polymorphism (SNP) detection, splice junction identification, and gene expression level measurement. However, the RNA-seq sequence data can be biased during library constructions resulting in incorrect data for SNP, splice junction, and gene expression studies. Here, we developed new library preparation methods to limit such biases.

Results

A whole transcriptome library prepared for the SOLiD system displayed numerous read duplications (pile-ups) and gaps in known exons. The pile-ups and gaps of the whole transcriptome library caused a loss of SNP and splice junction information and reduced the quality of gene expression results. Further, we found clear sequence biases for both 5' and 3' end reads in the whole transcriptome library. To remove this bias, RNaseIII fragmentation was replaced with heat fragmentation. For adaptor ligation, T4 Polynucleotide Kinase (T4PNK) was used following heat fragmentation. However, its kinase and phosphatase activities introduced additional sequence biases. To minimize them, we used OptiKinase before T4PNK. Our study further revealed the specific target sequences of RNaseIII and T4PNK.

Conclusions

Our results suggest that the heat fragmentation removed the RNaseIII sequence bias and significantly reduced the pile-ups and gaps. OptiKinase minimized the T4PNK sequence biases and removed most of the remaining pile-ups and gaps, thus maximizing the quality of RNA-seq data.

Reviewers

This article was reviewed by Dr. A. Kolodziejczyk (nominated by Dr. Sarah Teichmann), Dr. Eugene Koonin, and Dr. Christoph Adami. For the full reviews, see the Reviewers' Comments section.

DNA damage induces targeted, genome-wide variation of poly(A) sites in budding yeast

Systemic response to DNA damage and other stresses is a complex process that includes changes in the regulation and activity of nearly all stages of gene expression. One gene regulatory mechanism used by eukaryotes is selection among alternative transcript isoforms that differ in polyadenylation [poly(A)] sites, resulting in changes either to the coding sequence or to portions of the 3' UTR that govern translation, stability, and localization. To determine the extent to which this means of regulation is used in response to DNA damage, we conducted a global analysis of poly(A) site usage in Saccharomyces cerevisiae after exposure to the UV mimetic, 4-nitroquinoline 1-oxide (4NQO). Two thousand thirty-one genes were found to have significant variation in poly(A) site distributions following 4NQO treatment, with a strong bias toward loss of short transcripts, including many with poly(A) sites located within the protein coding sequence (CDS). We further explored one possible mechanism that could contribute to the widespread differences in mRNA isoforms. The change in poly(A) site profile was associated with an inhibition of cleavage and polyadenylation in cell extract and a decrease in the levels of several key subunits in the mRNA 3'-end processing complex. Sequence analysis identified differences in the cis-acting elements that flank putatively suppressed and enhanced poly(A) sites, suggesting a mechanism that could discriminate between variable and constitutive poly(A) sites. Our analysis indicates that variation in mRNA length is an important part of the regulatory response to DNA damage.

Species and condition specific adaptation of the transcriptional landscapes in Candida albicans and Candida dubliniensis

BACKGROUND

Although Candida albicans and Candida dubliniensis are most closely related, both species behave significantly different with respect to morphogenesis and virulence. In order to gain further insight into the divergent routes for morphogenetic adaptation in both species, we investigated qualitative along with quantitative differences in the transcriptomes of both organisms by cDNA deep sequencing.

RESULTS

Following genome-associated assembly of sequence reads we were able to generate experimentally verified databases containing 6016 and 5972 genes for C. albicans and C. dubliniensis, respectively. About 95% of the transcriptionally active regions (TARs) contain open reading frames while the remaining TARs most likely represent non-coding RNAs. Comparison of our annotations with publically available gene models for C. albicans and C. dubliniensis confirmed approximately 95% of already predicted genes, but also revealed so far unknown novel TARs in both species. Qualitative cross-species analysis of these databases revealed in addition to 5802 orthologs also 399 and 49 species-specific protein coding genes for C. albicans and C. dubliniensis, respectively. Furthermore, quantitative transcriptional profiling using RNA-Seq revealed significant differences in the expression of orthologs across both species. We defined a core subset of 84 hyphal-specific genes required for both species, as well as a set of 42 genes that seem to be specifically induced during hyphal morphogenesis in C. albicans.

CONCLUSIONS

Species-specific adaptation in C. albicans and C. dubliniensis is governed by individual genetic repertoires but also by altered regulation of conserved orthologs on the transcriptional level.

Extensive transcript diversity and novel upstream open reading frame regulation in yeast

To understand the diversity of transcripts in yeast (Saccharomyces cerevisiae) we analyzed the transcriptional landscapes for cells grown under 18 different environmental conditions. Each sample was analyzed using RNA-sequencing, and a total of 670,446,084 uniquely mapped reads and 377,263 poly-adenylated end tags were produced. Consistent with previous studies, we find that the majority of yeast genes are expressed under one or more different conditions. By directly comparing the 5′ and 3′ ends of the transcribed regions, we find extensive differences in transcript ends across many conditions, especially those of stationary phase, growth in grape juice, and salt stimulation, suggesting differential choice of transcription start and stop sites is pervasive in yeast. Relative to the exponential growth condition (i.e., YPAD), transcripts differing at the 5′ ends and 3′ ends are predicted to differ in their annotated start codon in 21 genes and their annotated stop codon in 63 genes. Many (431) upstream open reading frames (uORFs) are found in alternate 5′ ends and are significantly enriched in transcripts produced during the salt response. Mutational analysis of five genes with uORFs revealed that two sets of uORFs increase the expression of a reporter construct, indicating a role in activation which had not been reported previously, whereas two other uORFs decreased expression. In addition, RNA binding protein motifs are statistically enriched for alternate ends under many conditions. Overall, these results demonstrate enormous diversity of transcript ends, and that this heterogeneity is regulated under different environmental conditions. Moreover, transcript end diversity has important biological implications for the regulation of gene expression. In addition, our data also serve as a valuable resource for the scientific community.

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.