Genome-wide analysis of poly(A) site selection in Schizosaccharomyces pombe

Transcriptional and post-transcriptional regulation of genes encoding secretory proteins in Aspergillus oryzae

Aspergillus oryzae, also known as the yellow koji mold, produces various hydrolytic enzymes that are widely used in different industries. Its high capacity to produce secretory proteins makes this filamentous fungus a suitable host for heterologous protein production. Amylolytic gene promoter is widely used to express heterologous genes in A. oryzae. The expression of this promoter is strictly regulated by several transcription factors, whose activation involves various factors. Furthermore, the expression levels of amylolytic and heterologous genes are post-transcriptionally regulated by mRNA degradation mechanisms in response to aberrant transcriptional termination or endoplasmic reticulum stress. This review discusses the transcriptional and post-transcriptional regulatory mechanisms controlling the expression of genes encoding secretory proteins in A. oryzae.

Pervasive mRNA uridylation in fission yeast is catalysed by both Cid1 and Cid16 terminal uridyltransferases

Messenger RNA uridylation is pervasive and conserved among eukaryotes, but the consequences of this modification for mRNA fate are still under debate. Utilising a simple model organism to study uridylation may facilitate efforts to understand the cellular function of this process. Here we demonstrate that uridylation can be detected using simple bioinformatics approach. We utilise it to unravel widespread transcript uridylation in fission yeast and demonstrate the contribution of both Cid1 and Cid16, the only two annotated terminal uridyltransferases (TUT-ases) in this yeast. To detect uridylation in transcriptome data, we used a RNA-sequencing (RNA-seq) library preparation protocol involving initial linker ligation to fragmented RNA-an approach borrowed from small RNA sequencing that was commonly used in older RNA-seq protocols. We next explored the data to detect uridylation marks. Our analysis show that uridylation in yeast is pervasive, similarly to the one in multicellular organisms. Importantly, our results confirm the role of the cytoplasmic uridyltransferase Cid1 as the primary uridylation catalyst. However, we also observed an auxiliary role of the second uridyltransferase, Cid16. Thus both fission yeast uridyltransferases are involved in mRNA uridylation. Intriguingly, we found no physiological phenotype of the single and double deletion mutants of cid1 and cid16 and only minimal impact of uridylation on steady-state mRNA levels. Our work establishes fission yeast as a potent model to study uridylation in a simple eukaryote, and we demonstrate that it is possible to detect uridylation marks in RNA-seq data without the need for specific methodologies.

Genetic screen for suppression of transcriptional interference identifies a gain-of-function mutation in Pol2 termination factor Seb1

The system of long noncoding RNA (lncRNA)-mediated transcriptional interference that represses fission yeast phosphate homoeostasis gene pho1 provides a sensitive readout of genetic influences on cotranscriptional 3'-processing and termination and a tool for discovery of regulators of this phase of the Pol2 transcription cycle. Here, we conducted a genetic screen for relief of transcriptional interference that unveiled a mechanism by which Pol2 termination is enhanced via a gain-of-function mutation, G476S, in the RNA-binding domain of an essential termination factor, Seb1. The genetic and physical evidence for gain-of-function is compelling: 1) seb1-G476S de-represses pho1 and tgp1, both of which are subject to lncRNA-mediated transcriptional interference; 2) seb1-G476S elicits precocious lncRNA transcription termination in response to lncRNA 5'-proximal poly(A) signals; 3) seb1-G476S derepression of pho1 is effaced by loss-of-function mutations in cleavage and polyadenylation factor (CPF) subunits and termination factor Rhn1; 4) synthetic lethality of seb1-G476S with pho1 derepressive mutants rpb1-CTD-S7A and aps1∆ is rescued by CPF/Rhn1 loss-of-function alleles; and 5) seb1-G476S elicits an upstream shift in poly(A) site preference in several messenger RNA genes. A crystal structure of the Seb1-G476S RNA-binding domain indicates potential for gain of contacts from Ser476 to RNA nucleobases. To our knowledge, this is a unique instance of a gain-of-function phenotype in a eukaryal transcription termination protein.

Nuclear mRNA maturation and mRNA export control: from trypanosomes to opisthokonts

The passage of mRNAs through the nuclear pores into the cytoplasm is essential in all eukaryotes. For regulation, mRNA export is tightly connected to the full machinery of nuclear mRNA processing, starting at transcription. Export competence of pre-mRNAs gradually increases by both transient and permanent interactions with multiple RNA processing and export factors. mRNA export is best understood in opisthokonts, with limited knowledge in plants and protozoa. Here, I review and compare nuclear mRNA processing and export between opisthokonts and Trypanosoma brucei. The parasite has many unusual features in nuclear mRNA processing, such as polycistronic transcription and trans-splicing. It lacks several nuclear complexes and nuclear-pore-associated proteins that in opisthokonts play major roles in mRNA export. As a consequence, trypanosome mRNA export control is not tight and export can even start co-transcriptionally. Whether trypanosomes regulate mRNA export at all, or whether leakage of immature mRNA to the cytoplasm is kept to a low level by a fast kinetics of mRNA processing remains to be investigated. mRNA export had to be present in the last common ancestor of eukaryotes. Trypanosomes are evolutionary very distant from opisthokonts and a comparison helps understanding the evolution of mRNA export.

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.

Common mechanism of transcription termination at coding and noncoding RNA genes in fission yeast

Termination of RNA polymerase II (RNAPII) transcription is a fundamental step of gene expression that is critical for determining the borders between genes. In budding yeast, termination at protein-coding genes is initiated by the cleavage/polyadenylation machinery, whereas termination of most noncoding RNA (ncRNA) genes occurs via the Nrd1-Nab3-Sen1 (NNS) pathway. Here, we find that NNS-like transcription termination is not conserved in fission yeast. Rather, genome-wide analyses show global recruitment of mRNA 3' end processing factors at the end of ncRNA genes, including snoRNAs and snRNAs, and that this recruitment coincides with high levels of Ser2 and Tyr1 phosphorylation on the RNAPII C-terminal domain. We also find that termination of mRNA and ncRNA transcription requires the conserved Ysh1/CPSF-73 and Dhp1/XRN2 nucleases, supporting widespread cleavage-dependent transcription termination in fission yeast. Our findings thus reveal that a common mode of transcription termination can produce functionally and structurally distinct types of polyadenylated and non-polyadenylated RNAs.

A compendium of conserved cleavage and polyadenylation events in mammalian genes

Cleavage and polyadenylation is essential for 3' end processing of almost all eukaryotic mRNAs. Recent studies have shown widespread alternative cleavage and polyadenylation (APA) events leading to mRNA isoforms with different 3' UTRs and/or coding sequences. Here, we present a compendium of conserved cleavage and polyadenylation sites (PASs) in mammalian genes, based on approximately 1.2 billion 3' end sequencing reads from more than 360 human, mouse, and rat samples. We show that ∼80% of mammalian mRNA genes contain at least one conserved PAS, and ∼50% have conserved APA events. PAS conservation generally reduces promiscuous 3' end processing, stabilizing gene expression levels across species. Conservation of APA correlates with gene age, gene expression features, and gene functions. Genes with certain functions, such as cell morphology, cell proliferation, and mRNA metabolism, are particularly enriched with conserved APA events. Whereas tissue-specific genes typically have a low APA rate, brain-specific genes tend to evolve APA. In addition, we show enrichment of mRNA destabilizing motifs in alternative 3' UTR sequences, leading to substantial differences in mRNA stability between 3' UTR isoforms. Using conserved PASs, we reveal sequence motifs surrounding APA sites and a preference of adenosine at the cleavage site. Furthermore, we show that mutations of U-rich motifs around the PAS often accompany APA profile differences between species. Analysis of lncRNA PASs indicates a mechanism of PAS fixation through evolution of A-rich motifs. Taken together, our results present a comprehensive view of PAS evolution in mammals, and a phylogenic perspective on APA functions.

A compendium of conserved cleavage and polyadenylation events in mammalian genes

Cleavage and polyadenylation is essential for 3' end processing of almost all eukaryotic mRNAs. Recent studies have shown widespread alternative cleavage and polyadenylation (APA) events leading to mRNA isoforms with different 3' UTRs and/or coding sequences. Here, we present a compendium of conserved cleavage and polyadenylation sites (PASs) in mammalian genes, based on approximately 1.2 billion 3' end sequencing reads from more than 360 human, mouse, and rat samples. We show that ∼80% of mammalian mRNA genes contain at least one conserved PAS, and ∼50% have conserved APA events. PAS conservation generally reduces promiscuous 3' end processing, stabilizing gene expression levels across species. Conservation of APA correlates with gene age, gene expression features, and gene functions. Genes with certain functions, such as cell morphology, cell proliferation, and mRNA metabolism, are particularly enriched with conserved APA events. Whereas tissue-specific genes typically have a low APA rate, brain-specific genes tend to evolve APA. In addition, we show enrichment of mRNA destabilizing motifs in alternative 3' UTR sequences, leading to substantial differences in mRNA stability between 3' UTR isoforms. Using conserved PASs, we reveal sequence motifs surrounding APA sites and a preference of adenosine at the cleavage site. Furthermore, we show that mutations of U-rich motifs around the PAS often accompany APA profile differences between species. Analysis of lncRNA PASs indicates a mechanism of PAS fixation through evolution of A-rich motifs. Taken together, our results present a comprehensive view of PAS evolution in mammals, and a phylogenic perspective on APA functions.

Long noncoding RNA repertoire and targeting by nuclear exosome, cytoplasmic exonuclease, and RNAi in fission yeast

Long noncoding RNAs (lncRNAs), which are longer than 200 nucleotides but often unstable, contribute a substantial and diverse portion to pervasive noncoding transcriptomes. Most lncRNAs are poorly annotated and understood, although several play important roles in gene regulation and diseases. Here we systematically uncover and analyze lncRNAs in Schizosaccharomyces pombe. Based on RNA-seq data from twelve RNA-processing mutants and nine physiological conditions, we identify 5775 novel lncRNAs, nearly 4× the previously annotated lncRNAs. The expression of most lncRNAs becomes strongly induced under the genetic and physiological perturbations, most notably during late meiosis. Most lncRNAs are cryptic and suppressed by three RNA-processing pathways: the nuclear exosome, cytoplasmic exonuclease, and RNAi. Double-mutant analyses reveal substantial coordination and redundancy among these pathways. We classify lncRNAs by their dominant pathway into cryptic unstable transcripts (CUTs), Xrn1-sensitive unstable transcripts (XUTs), and Dicer-sensitive unstable transcripts (DUTs). XUTs and DUTs are enriched for antisense lncRNAs, while CUTs are often bidirectional and actively translated. The cytoplasmic exonuclease, along with RNAi, dampens the expression of thousands of lncRNAs and mRNAs that become induced during meiosis. Antisense lncRNA expression mostly negatively correlates with sense mRNA expression in the physiological, but not the genetic conditions. Intergenic and bidirectional lncRNAs emerge from nucleosome-depleted regions, upstream of positioned nucleosomes. Our results highlight both similarities and differences to lncRNA regulation in budding yeast. This broad survey of the lncRNA repertoire and characteristics in S. pombe, and the interwoven regulatory pathways that target lncRNAs, provides a rich framework for their further functional analyses.

Tandem 3' UTR Patterns and Gene Expression Profiles of Marc-145 Cells During PRRSV Infection

Porcine reproductive and respiratory syndrome virus (PRRSV) causes substantial economic losses to the global pig industry. Alternative polyadenylation (APA) is a mechanism that diversifies gene expression, which is important for tumorigenesis, development, and cell differentiation. However, it is unclear whether APA plays a role in the course of PRRSV infection. To address this issue, in this study we carried out a whole-genome transcriptome analysis of PRRSV-infected Marc-145 African green monkey kidney cells and identified 185 APA switching genes and 393 differentially expressed genes (DEGs). Most of these genes were involved in cellular process, metabolism, and biological regulation, and there was some overlap between the two gene sets. DEGs were found to be more directly involved in the antiviral response than APA genes. These findings provide insight into the dynamics of host gene regulation during PRRSV infection and a basis for elucidating the pathogenesis of PRRSV.

Codon usage biases co-evolve with transcription termination machinery to suppress premature cleavage and polyadenylation

Codon usage biases are found in all genomes and influence protein expression levels. The codon usage effect on protein expression was thought to be mainly due to its impact on translation. Here, we show that transcription termination is an important driving force for codon usage bias in eukaryotes. Using Neurospora crassa as a model organism, we demonstrated that introduction of rare codons results in premature transcription termination (PTT) within open reading frames and abolishment of full-length mRNA. PTT is a wide-spread phenomenon in Neurospora, and there is a strong negative correlation between codon usage bias and PTT events. Rare codons lead to the formation of putative poly(A) signals and PTT. A similar role for codon usage bias was also observed in mouse cells. Together, these results suggest that codon usage biases co-evolve with the transcription termination machinery to suppress premature termination of transcription and thus allow for optimal gene expression.

PomBase: The Scientific Resource for Fission Yeast

The fission yeast Schizosaccharomyces pombe has become well established as a model species for studying conserved cell-level biological processes, especially the mechanics and regulation of cell division. PomBase integrates the S. pombe genome sequence with traditional genetic, molecular, and cell biological experimental data as well as the growing body of large datasets generated by emerging high-throughput methods. This chapter provides insight into the curation philosophy and data organization at PomBase, and provides a guide to using PomBase for infrequent visitors and anyone considering exploring S. pombe in their research.

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 RpII215^C4 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.

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.

Regulation by 3'-Untranslated Regions

3'-untranslated regions (3'-UTRs) are the noncoding parts of mRNAs. Compared to yeast, in humans, median 3'-UTR length has expanded approximately tenfold alongside an increased generation of alternative 3'-UTR isoforms. In contrast, the number of coding genes, as well as coding region length, has remained similar. This suggests an important role for 3'-UTRs in the biology of higher organisms. 3'-UTRs are best known to regulate diverse fates of mRNAs, including degradation, translation, and localization, but they can also function like long noncoding or small RNAs, as has been shown for whole 3'-UTRs as well as for cleaved fragments. Furthermore, 3'-UTRs determine the fate of proteins through the regulation of protein-protein interactions. They facilitate cotranslational protein complex formation, which establishes a role for 3'-UTRs as evolved eukaryotic operons. Whereas bacterial operons promote the interaction of two subunits, 3'-UTRs enable the formation of protein complexes with diverse compositions. All of these 3'-UTR functions are accomplished by effector proteins that are recruited by RNA-binding proteins that bind to 3'-UTR cis-elements. In summary, 3'-UTRs seem to be major players in gene regulation that enable local functions, compartmentalization, and cooperativity, which makes them important tools for the regulation of phenotypic diversity of higher organisms.

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.

Analysis of RNA Metabolism in Fission Yeast

Here we focus on the biogenesis and function of messenger RNA (mRNA) in fission yeast cells. Following a general introduction that also briefly touches on other classes of RNA, we provide an overview of methods used to analyze mRNAs throughout their life cycles.

The conserved protein Seb1 drives transcription termination by binding RNA polymerase II and nascent RNA

Termination of RNA polymerase II (Pol II) transcription is an important step in the transcription cycle, which involves the dislodgement of polymerase from DNA, leading to release of a functional transcript. Recent studies have identified the key players required for this process and showed that a common feature of these proteins is a conserved domain that interacts with the phosphorylated C-terminus of Pol II (CTD-interacting domain, CID). However, the mechanism by which transcription termination is achieved is not understood. Using genome-wide methods, here we show that the fission yeast CID-protein Seb1 is essential for termination of protein-coding and non-coding genes through interaction with S2-phosphorylated Pol II and nascent RNA. Furthermore, we present the crystal structures of the Seb1 CTD- and RNA-binding modules. Unexpectedly, the latter reveals an intertwined two-domain arrangement of a canonical RRM and second domain. These results provide important insights into the mechanism underlying eukaryotic transcription termination.

Genome-Wide Analysis and Functional Characterization of the Polyadenylation Site in Pigs Using RNAseq Data

Polyadenylation, a critical step in the production of mature mRNA for translation in most eukaryotes, involves cleavage and poly(A) tail addition at the 3′ end of mRNAs at the polyadenylation site (PAS). Sometimes, one gene can have more than one PAS, which can produce the alternative polyadenylation (APA) phenomenon and affect the stability, localization and translation of the mRNA. In this study, we discovered 28,363 PASs using pig RNAseq data, with 13,033 located in 7,403 genes. Among the genes, 41% were identified to have more than one PAS. PAS distribution analysis indicated that the PAS position was highly variable in genes. Additionally, the analysis of RNAseq data from the liver and testis showed a difference in their PAS number and usage. RT-PCR and qRT-PCR were performed to confirm our findings by detecting the expression of 3′UTR isoforms for five candidate genes. The analysis of RNAseq data under a different androstenone level and salmonella inoculation indicated that the functional usage of PAS might participate in the immune response and may be related to the androstenone level in pigs. This study provides new insights into pig PAS and facilitates further functional research of PAS.

Polyadenylation site selection: linking transcription and RNA processing via a conserved carboxy-terminal domain (CTD)-interacting protein

Despite the fact that the process of mRNA polyadenylation has been known for more than 40 years, a detailed understating of the mechanism underlying polyadenylation site selection is still far from complete. As 3' end processing is intimately associated with RNA polymerase II (RNAPII) transcription, factors that can successively interact with the transcription machinery and recognize cis-acting sequences on the nascent pre-mRNA would be well suited to contribute to poly(A) site selection. Studies using the fission yeast Schizosaccharomyces pombe have recently identified Seb1, a protein that shares homology with Saccharomyces cerevisiae Nrd1 and human SCAF4/8, and that is critical for poly(A) site selection. Seb1 binds to the C-terminal domain (CTD) of RNAPII via a conserved CTD-interaction domain and recognizes specific sequence motifs clustered downstream of the polyadenylation site on the uncleaved pre-mRNA. In this short review, we summarize insights into Seb1-dependent poly(A) site selection and discuss some unanswered questions regarding its molecular mechanism and conservation.

Determinants of RNA metabolism in the Schizosaccharomyces pombe genome

To decrypt the regulatory code of the genome, sequence elements must be defined that determine the kinetics of RNA metabolism and thus gene expression. Here, we attempt such decryption in an eukaryotic model organism, the fission yeast S. pombe. We first derive an improved genome annotation that redefines borders of 36% of expressed mRNAs and adds 487 non‐coding RNAs (ncRNAs). We then combine RNA labeling in vivo with mathematical modeling to obtain rates of RNA synthesis and degradation for 5,484 expressed RNAs and splicing rates for 4,958 introns. We identify functional sequence elements in DNA and RNA that control RNA metabolic rates and quantify the contributions of individual nucleotides to RNA synthesis, splicing, and degradation. Our approach reveals distinct kinetics of mRNA and ncRNA metabolism, separates antisense regulation by transcription interference from RNA interference, and provides a general tool for studying the regulatory code of genomes.

Genome-Wide Study of mRNA Isoform Half-Lives

In eukaryotes, RNA polymerase II-driven transcription and processing results in the formation of numerous mRNA 3' isoforms that for any given gene may differ from one another by as little as a single nucleotide. These 3' isoforms can vary in physical properties that may affect their function and stability. Here, we outline a systematic framework to measure individual mRNA 3' isoform half-lives on a genome-wide level in S. cerevisiae. Our approach utilizes the Anchor-Away system to sequester RNA polymerase II (Pol II) in the cytoplasm followed by direct single-molecule RNA sequencing to generate a highly detailed view of 3' isoform stability under most physiological conditions without many of the adverse effects associated with commonly used alternative approaches.

Whole transcriptome analysis with sequencing: methods, challenges and potential solutions

Whole transcriptome analysis plays an essential role in deciphering genome structure and function, identifying genetic networks underlying cellular, physiological, biochemical and biological systems and establishing molecular biomarkers that respond to diseases, pathogens and environmental challenges. Here, we review transcriptome analysis methods and technologies that have been used to conduct whole transcriptome shotgun sequencing or whole transcriptome tag/target sequencing analyses. We focus on how adaptors/linkers are added to both 5' and 3' ends of mRNA molecules for cloning or PCR amplification before sequencing. Challenges and potential solutions are also discussed. In brief, next generation sequencing platforms have accelerated releases of the large amounts of gene expression data. It is now time for the genome research community to assemble whole transcriptomes of all species and collect signature targets for each gene/transcript, and thus use known genes/transcripts to determine known transcriptomes directly in the near future.

The fission yeast MTREC complex targets CUTs and unspliced pre-mRNAs to the nuclear exosome

Cryptic unstable transcripts (CUTs) are rapidly degraded by the nuclear exosome. However, the mechanism by which they are recognized and targeted to the exosome is not fully understood. Here we report that the MTREC complex, which has recently been shown to promote degradation of meiotic mRNAs and regulatory ncRNAs, is also the major nuclear exosome targeting complex for CUTs and unspliced pre-mRNAs in Schizosaccharomyces pombe. The MTREC complex specifically binds to CUTs, meiotic mRNAs and unspliced pre-mRNA transcripts and targets these RNAs for degradation by the nuclear exosome, while the TRAMP complex has only a minor role in this process. The MTREC complex physically interacts with the nuclear exosome and with various RNA-binding and RNA-processing complexes, coupling RNA processing to the RNA degradation machinery. Our study reveals the central role of the evolutionarily conserved MTREC complex in RNA quality control, and in the recognition and elimination of CUTs.

APASdb: a database describing alternative poly(A) sites and selection of heterogeneous cleavage sites downstream of poly(A) signals

Increasing amounts of genes have been shown to utilize alternative polyadenylation (APA) 3'-processing sites depending on the cell and tissue type and/or physiological and pathological conditions at the time of processing, and the construction of genome-wide database regarding APA is urgently needed for better understanding poly(A) site selection and APA-directed gene expression regulation for a given biology. Here we present a web-accessible database, named APASdb (http://mosas.sysu.edu.cn/utr), which can visualize the precise map and usage quantification of different APA isoforms for all genes. The datasets are deeply profiled by the sequencing alternative polyadenylation sites (SAPAS) method capable of high-throughput sequencing 3'-ends of polyadenylated transcripts. Thus, APASdb details all the heterogeneous cleavage sites downstream of poly(A) signals, and maintains near complete coverage for APA sites, much better than the previous databases using conventional methods. Furthermore, APASdb provides the quantification of a given APA variant among transcripts with different APA sites by computing their corresponding normalized-reads, making our database more useful. In addition, APASdb supports URL-based retrieval, browsing and display of exon-intron structure, poly(A) signals, poly(A) sites location and usage reads, and 3'-untranslated regions (3'-UTRs). Currently, APASdb involves APA in various biological processes and diseases in human, mouse and zebrafish.

A rapid, simple, and inexpensive method for the preparation of strand-specific RNA-Seq libraries

High-throughput sequencing of short cDNA tags, or RNA-Seq, has become a staple of genome-wide gene expression studies in plants. RNA-Seq libraries necessarily contain tags that correspond to the mRNA-poly(A) junction, or polyadenylation site, and thus may be mined for data that can help study alternative polyadenylation. This report presents a simple, rapid, and inexpensive method for preparing strand-specific RNA-Seq libraries from varying quantities of total RNA.

PomBase 2015: updates to the fission yeast database

PomBase (http://www.pombase.org) is the model organism database for the fission yeast Schizosaccharomyces pombe. PomBase provides a central hub for the fission yeast community, supporting both exploratory and hypothesis-driven research. It provides users easy access to data ranging from the sequence level, to molecular and phenotypic annotations, through to the display of genome-wide high-throughput studies. Recent improvements to the site extend annotation specificity, improve usability and allow for monthly data updates. Both in-house curators and community researchers provide manually curated data to PomBase. The genome browser provides access to published high-throughput data sets and the genomes of three additional Schizosaccharomyces species (Schizosaccharomyces cryophilus, Schizosaccharomyces japonicus and Schizosaccharomyces octosporus).

Dicer promotes transcription termination at sites of replication stress to maintain genome stability

Nuclear RNAi is an important regulator of transcription and epigenetic modification, but the underlying mechanisms remain elusive. Using a genome-wide approach in the fission yeast S. pombe, we have found that Dcr1, but not other components of the canonical RNAi pathway, promotes the release of Pol II from the 3? end of highly transcribed genes, and, surprisingly, from antisense transcription of rRNA and tRNA genes, which are normally transcribed by Pol I and Pol III. These Dcr1-terminated loci correspond to sites of replication stress and DNA damage, likely resulting from transcription-replication collisions. At the rDNA loci, release of Pol II facilitates DNA replication and prevents homologous recombination, which would otherwise lead to loss of rDNA repeats especially during meiosis. Our results reveal a novel role for Dcr1-mediated transcription termination in genome maintenance and may account for widespread regulation of genome stability by nuclear RNAi in higher eukaryotes.

Motif types, motif locations and base composition patterns around the RNA polyadenylation site in microorganisms, plants and animals

Background

The polyadenylation of RNA is critical for gene functioning, but the conserved sequence motifs (often called signal or signature motifs), motif locations and abundances, and base composition patterns around mRNA polyadenylation [poly(A)] sites are still uncharacterized in most species. The evolutionary tendency for poly(A) site selection is still largely unknown.

Results

We analyzed the poly(A) site regions of 31 species or phyla. Different groups of species showed different poly(A) signal motifs: UUACUU at the poly(A) site in the parasite Trypanosoma cruzi; UGUAAC (approximately 13 bases upstream of the site) in the alga Chlamydomonas reinhardtii; UGUUUG (or UGUUUGUU) at mainly the fourth base downstream of the poly(A) site in the parasite Blastocystis hominis; and AAUAAA at approximately 16 bases and approximately 19 bases upstream of the poly(A) site in animals and plants, respectively. Polyadenylation signal motifs are usually several hundred times more abundant around poly(A) sites than in whole genomes. These predominant motifs usually had very specific locations, whether upstream of, at, or downstream of poly(A) sites, depending on the species or phylum. The poly(A) site was usually an adenosine (A) in all analyzed species except for B. hominis, and there was weak A predominance in C. reinhardtii. Fungi, animals, plants, and the protist Phytophthora infestans shared a general base abundance pattern (or base composition pattern) of “U-rich—A-rich—U-rich—Poly(A) site—U-rich regions”, or U-A-U-A-U for short, with some variation for each kingdom or subkingdom.

Conclusion

This study identified the poly(A) signal motifs, motif locations, and base composition patterns around mRNA poly(A) sites in protists, fungi, plants, and animals and provided insight into poly(A) site evolution.

Effects of codon optimization on the mRNA levels of heterologous genes in filamentous fungi

Filamentous fungi, particularly Aspergillus species, have recently attracted attention as host organisms for recombinant protein production. Because the secretory yields of heterologous proteins are generally low compared with those of homologous proteins or proteins from closely related fungal species, several strategies to produce substantial amounts of recombinant proteins have been conducted. Codon optimization is a powerful tool for improving the production levels of heterologous proteins. Although codon optimization is generally believed to improve the translation efficiency of heterologous genes without affecting their mRNA levels, several studies have indicated that codon optimization causes an increase in the steady-state mRNA levels of heterologous genes in filamentous fungi. However, the mechanism that determines the low mRNA levels when native heterologous genes are expressed was poorly understood. We recently showed that the transcripts of heterologous genes are polyadenylated prematurely within the coding region and that the heterologous gene transcripts can be stabilized significantly by codon optimization, which is probably attributable to the prevention of premature polyadenylation in Aspergillus oryzae. In this review, we describe the detailed mechanism of premature polyadenylation and the rapid degradation of mRNA transcripts derived from heterologous genes in filamentous fungi.