A large fraction of protein-coding genes in metazoans undergo alternative pre-mRNA splicing in tissue- or cell-type-specific manners. Recent genome-wide approaches have identified many putative-binding sites for some of tissue-specific trans-acting splicing regulators. However, the mechanisms of splicing regulation in vivo remain largely unknown. To elucidate the modes of splicing regulation by the neuron-specific CELF family RNA-binding protein UNC-75 in Caenorhabditis elegans, we performed deep sequencing of poly(A)(+) RNAs from the unc-75(+)- and unc-75-mutant worms and identified more than 20 cassette and mutually exclusive exons repressed or activated by UNC-75. Motif searches revealed that (G/U)UGUUGUG stretches are enriched in the upstream and downstream introns of the UNC-75-repressed and -activated exons, respectively. Recombinant UNC-75 protein specifically binds to RNA fragments carrying the (G/U)UGUUGUG stretches in vitro. Bi-chromatic fluorescence alternative splicing reporters revealed that the UNC-75-target exons are regulated in tissue-specific and (G/U)UGUUGUG element-dependent manners in vivo. The unc-75 mutation affected the splicing reporter expression specifically in the nervous system. These results indicate that UNC-75 regulates alternative splicing of its target exons in neuron-specific and position-dependent manners through the (G/U)UGUUGUG elements in C. elegans. This study thus reveals the repertoire of target events for the CELF family in the living organism.

Alternative splicing of pre-mRNAs can regulate expression of protein-coding genes by generating unproductive mRNAs rapidly degraded by nonsense-mediated mRNA decay (NMD). Many of the genes directly regulated by alternative splicing coupled with NMD (AS-NMD) are related to RNA metabolism, but the repertoire of genes regulated by AS-NMD in vivo is to be determined. Here, we analyzed transcriptome data of wild-type and NMD-defective mutant strains of the nematode worm Caenorhabditis elegans and demonstrate that eight of the 82 cytoplasmic ribosomal protein (rp) genes generate unproductively spliced mRNAs. Knockdown of any of the eight rp genes exerted a dynamic and compensatory effect on alternative splicing of its own transcript and inverse effects on that of the other rp genes. A large subunit protein L10a, termed RPL-1 in nematodes, directly and specifically binds to an evolutionarily conserved 39-nt stretch termed L10ARE between the two alternative 5' splice sites in its own pre-mRNA to switch the splice site choice. Furthermore, L10ARE-mediated splicing autoregulation of the L10a-coding gene is conserved in vertebrates. These results indicate that L10a is an evolutionarily conserved splicing regulator and that homeostasis of a subset of the rp genes are regulated at the level of pre-mRNA splicing in vivo.

Combining our full-length cDNA method and the massively parallel sequencing technology, we developed a simple method to collect precise positional information of transcriptional start sites (TSSs) together with digital information of the gene-expression levels in a high throughput manner. We applied this method to observe gene-expression changes in a colon cancer cell line cultured in normoxic and hypoxic conditions. We generated more than 100 million 36-base TSS-tag sequences and revealed comprehensive features of hypoxia responsive alterations in the transcriptional landscape of the human genome. The features include presence of inducible 'hot regions' in 54 genomic regions, 220 novel hypoxia inducible promoters that may drive non-protein-coding transcripts, 191 hypoxia responsive alternative promoters and detailed views of 120 novel as well as known hypoxia responsive genes. We further analyzed hypoxic response of different cells using additional 60 million TSS-tags and found that the degree of the gene-expression changes were different among cell lines, possibly reflecting cellular robustness against hypoxia. The novel dynamic figure of the human gene transcriptome will deepen our understanding of the transcriptional program of the human genome as well as bringing new insights into the biology of cancer cells in hypoxia.

Gene overlap plays various regulatory functions on transcriptional and post-transcriptional levels. Most current studies focus on protein-coding genes overlapping with non-protein-coding counterparts, the so called natural antisense transcripts. Considerably less is known about the role of gene overlap in the case of two protein-coding genes. Here, we provide OverGeneDB, a database of human and mouse 5' end protein-coding overlapping genes. The database contains 582 human and 113 mouse gene pairs that are transcribed using overlapping promoters in at least one analyzed library. Gene pairs were identified based on the analysis of the transcription start site (TSS) coordinates in 73 human and 10 mouse organs, tissues and cell lines. Beside TSS data, resources for 26 human lung adenocarcinoma cell lines also contain RNA-Seq and ChIP-Seq data for seven histone modifications and RNA Polymerase II activity. The collected data revealed that the overlap region is rarely conserved between the studied species and tissues. In ∼50% of the overlapping genes, transcription started explicitly in the overlap regions. In the remaining half of overlapping genes, transcription was initiated both from overlapping and non-overlapping TSSs. OverGeneDB is accessible at http://overgenedb.amu.edu.pl.

In budding yeast, Set2 catalyzes di- and trimethylation of H3K36 (H3K36me2 and H3K36me3) via an interaction between its Set2-Rpb1 interaction (SRI) domain and C-terminal repeats of RNA polymerase II (Pol2) phosphorylated at Ser2 and Ser5 (CTD-S2,5-P). H3K36me2 is sufficient for recruitment of the Rpd3S histone deacetylase complex to repress cryptic transcription from transcribed regions. In fission yeast, Set2 is also responsible for H3K36 methylation, which represses a subset of RNAs including heterochromatic and subtelomeric RNAs, at least in part via recruitment of Clr6 complex II, a homolog of Rpd3S. Here, we show that CTD-S2P-dependent interaction of fission yeast Set2 with Pol2 via the SRI domain is required for formation of H3K36me3, but not H3K36me2. H3K36me3 silenced heterochromatic and subtelomeric transcripts mainly through post-transcriptional and transcriptional mechanisms, respectively, whereas H3K36me2 was not enough for silencing. Clr6 complex II appeared not to be responsible for heterochromatic silencing by H3K36me3. Our results demonstrate that H3K36 methylation has multiple outputs in fission yeast; these findings provide insights into the distinct roles of H3K36 methylation in metazoans, which have different enzymes for synthesis of H3K36me1/2 and H3K36me3.

To identify and characterize transcript structures ranging from transcriptional start sites (TSSs) to poly(A)-addition sites (PASs), we constructed and analyzed human TSS/PAS mate pair full-length cDNA libraries from 14 tissue types and four cell lines. The collected information enabled us to define TSS cluster (TSC) and PAS cluster (PAC) relationships for a total of 8530/9400 RefSeq genes, as well as 4251/5618 of their putative alternative promoters/terminators and 4619/4605 intervening transcripts, respectively. Analyses of the putative alternative TSCs and alternative PACs revealed that their selection appeared to be mostly independent, with rare exceptions. In those exceptional cases, pairs of transcript units rarely overlapped one another and were occasionally separated by Rad21/CTCF. We also identified a total of 172 similar cases in which TSCs and PACs spanned adjacent but distinct genes. In these cases, different transcripts may utilize different functional units of a particular gene or of adjacent genes. This approach was also useful for identifying fusion gene transcripts in cancerous cells. Furthermore, we could construct cDNA libraries in which 3'-end mate pairs were distributed randomly over the transcripts. These libraries were useful for assembling the internal structure of previously uncharacterized alternative promoter products, as well as intervening transcripts.

Long-read whole-genome sequencing analysis of DNA methylation would provide useful information on the chromosomal context of gene expression regulation. Here we describe the development of a method that improves the read length generated by using the bisulfite-sequencing-based approach. In this method, we combined recently developed enzymatic base conversion, where an unmethylated cytosine (C) should be converted to thymine (T), with nanopore sequencing. After methylation-sensitive base conversion, the sequencing library was constructed using long-range polymerase chain reaction. This type of analysis is possible using a minimum of 1 ng genomic DNA, and an N50 read length of 3.4-7.6 kb is achieved. To analyze the produced data, which contained a substantial number of base mismatches due to sequence conversion and an inaccurate base read of the nanopore sequencing, a new analytical pipeline was constructed. To demonstrate the performance of long-read methylation sequencing, breast cancer cell lines and clinical specimens were subjected to analysis, which revealed the chromosomal methylation context of key cancer-related genes, allele-specific methylated genes, and repetitive or deletion regions. This method should convert the intractable specimens for which the amount of available genomic DNA is limited to the tractable targets.