A code for transcription initiation in mammalian genomes

Widespread 3'UTR capped RNAs derive from G-rich regions in proximity to AGO2 binding sites

The 3' untranslated region (3'UTR) plays a crucial role in determining mRNA stability, localisation, translation and degradation. Cap analysis of gene expression (CAGE), a method for the detection of capped 5' ends of mRNAs, additionally reveals a large number of apparently 5' capped RNAs derived from locations within the body of the transcript, including 3'UTRs. Here, we provide direct evidence that these 3'UTR-derived RNAs are indeed capped and widespread in mammalian cells. By using a combination of AGO2 enhanced individual nucleotide resolution UV crosslinking and immunoprecipitation (eiCLIP) and CAGE following siRNA treatment, we find that these 3'UTR-derived RNAs likely originate from AGO2-binding sites, and most often occur at locations with G-rich motifs bound by the RNA-binding protein UPF1. High-resolution imaging and long-read sequencing analysis validate several 3'UTR-derived RNAs, showcase their variable abundance and show that they may not co-localise with the parental mRNAs. Taken together, we provide new insights into the origin and prevalence of 3'UTR-derived RNAs, show the utility of CAGE-seq for their genome-wide detection and provide a rich dataset for exploring new biology of a poorly understood new class of RNAs.

The DoGA consortium expression atlas of promoters and genes in 100 canine tissues

The dog, Canis lupus familiaris, is an important model for studying human diseases. Unlike many model organisms, the dog genome has a comparatively poor functional annotation, which hampers gene discovery for development, morphology, disease, and behavior. To fill this gap, we established a comprehensive tissue biobank for both the dog and wolf samples. The biobank consists of 5485 samples representing 132 tissues from 13 dogs, 12 dog embryos, and 24 wolves. In a subset of 100 tissues from nine dogs and 12 embryos, we characterized gene expression activity for each promoter, including alternative and novel, i.e., previously not annotated, promoter regions, using the 5' targeting RNA sequencing technology STRT2-seq. We identified over 100,000 promoter region candidates in the recent canine genome assembly, CanFam4, including over 45,000 highly reproducible sites with gene expression and respective tissue enrichment levels. We provide a promoter and gene expression atlas with interactive, open data resources, including a data coordination center and genome browser track hubs. We demonstrated the applicability of Dog Genome Annotation (DoGA) data and resources using multiple examples spanning canine embryonic development, morphology and behavior, and diseases across species.

High-quality chromosome scale genome assemblies of two important Sorghum inbred lines, Tx2783 and RTx436

Sorghum bicolor (L.) Moench is a significant grass crop globally, known for its genetic diversity. High quality genome sequences are needed to capture the diversity. We constructed high-quality, chromosome-level genome assemblies for two vital sorghum inbred lines, Tx2783 and RTx436. Through advanced single-molecule techniques, long-read sequencing and optical maps, we improved average sequence continuity 19-fold and 11-fold higher compared to existing Btx623 v3.0 reference genome and obtained 19 and 18 scaffolds (N50 of 25.6 and 14.4) for Tx2783 and RTx436, respectively. Our gene annotation efforts resulted in 29 612 protein-coding genes for the Tx2783 genome and 29 265 protein-coding genes for the RTx436 genome. Comparative analyses with 26 plant genomes which included 18 sorghum genomes and 8 outgroup species identified around 31 210 protein-coding gene families, with about 13 956 specific to sorghum. Using representative models from gene trees across the 18 sorghum genomes, a total of 72 579 pan-genes were identified, with 14% core, 60% softcore and 26% shell genes. We identified 99 genes in Tx2783 and 107 genes in RTx436 that showed functional enrichment specifically in binding and metabolic processes, as revealed by the GO enrichment Pearson Chi-Square test. We detected 36 potential large inversions in the comparison between the BTx623 Bionano map and the BTx623 v3.1 reference sequence. Strikingly, these inversions were notably absent when comparing Tx2783 or RTx436 with the BTx623 Bionano map. These inversion were mostly in the pericentromeric region which is known to have low complexity regions and harder to assemble and suggests the presence of potential artifacts in the public BTx623 reference assembly. Furthermore, in comparison to Tx2783, RTx436 exhibited 324 883 additional Single Nucleotide Polymorphisms (SNPs) and 16 506 more Insertions/Deletions (INDELs) when using BTx623 as the reference genome. We also characterized approximately 348 nucleotide-binding leucine-rich repeat (NLR) disease resistance genes in the two genomes. These high-quality genomes serve as valuable resources for discovering agronomic traits and structural variation studies.

Characterizing a complex CT-rich haplotype in intron 4 of SNCA using large-scale targeted amplicon long-read sequencing

Parkinson's disease (PD) is a common neurodegenerative disorder with a significant risk proportion driven by genetics. While much progress has been made, most of the heritability remains unknown. This is in-part because previous genetic studies have focused on the contribution of single nucleotide variants. More complex forms of variation, such as structural variants and tandem repeats, are already associated with several synucleinopathies. However, because more sophisticated sequencing methods are usually required to detect these regions, little is understood regarding their contribution to PD. One example is a polymorphic CT-rich region in intron 4 of the SNCA gene. This haplotype has been suggested to be associated with risk of Lewy Body (LB) pathology in Alzheimer's Disease and SNCA gene expression, but is yet to be investigated in PD. Here, we attempt to resolve this CT-rich haplotype and investigate its role in PD. We performed targeted PacBio HiFi sequencing of the region in 1375 PD cases and 959 controls. We replicate the previously reported associations and a novel association between two PD risk SNVs (rs356182 and rs5019538) and haplotype 4, the largest haplotype. Through quantitative trait locus analyzes we identify a significant haplotype 4 association with alternative CAGE transcriptional start site usage, not leading to significant differential SNCA gene expression in post-mortem frontal cortex brain tissue. Therefore, disease association in this locus might not be biologically driven by this CT-rich repeat region. Our data demonstrates the complexity of this SNCA region and highlights that further follow up functional studies are warranted.

High-density resolution of the Kaposi's sarcoma associated herpesvirus transcriptome identifies novel transcript isoforms generated by long-range transcription and alternative splicing

Kaposi's sarcoma-associated herpesvirus is the etiologic agent of Kaposi's sarcoma and two B-cell malignancies. Recent advancements in sequencing technologies have led to high resolution transcriptomes for several human herpesviruses that densely encode genes on both strands. However, for KSHV progress remained limited due to the overall low percentage of KSHV transcripts, even during lytic replication. To address this challenge, we have developed a target enrichment method to increase the KSHV-specific reads for both short- and long-read sequencing platforms. Furthermore, we combined this approach with the Transcriptome Resolution through Integration of Multi-platform Data (TRIMD) pipeline developed previously to annotate transcript structures. TRIMD first builds a scaffold based on long-read sequencing and validates each transcript feature with supporting evidence from Illumina RNA-Seq and deepCAGE sequencing data. Our stringent innovative approach identified 994 unique KSHV transcripts, thus providing the first high-density KSHV lytic transcriptome. We describe a plethora of novel coding and non-coding KSHV transcript isoforms with alternative untranslated regions, splice junctions and open-reading frames, thus providing deeper insights on gene expression regulation of KSHV. Interestingly, as described for Epstein-Barr virus, we identified transcription start sites that augment long-range transcription and may increase the number of latency-associated genes potentially expressed in KS tumors.

An atlas of transcribed enhancers across helper T cell diversity for decoding human diseases

Transcribed enhancer maps can reveal nuclear interactions underpinning each cell type and connect specific cell types to diseases. Using a 5' single-cell RNA sequencing approach, we defined transcription start sites of enhancer RNAs and other classes of coding and noncoding RNAs in human CD4+ T cells, revealing cellular heterogeneity and differentiation trajectories. Integration of these datasets with single-cell chromatin profiles showed that active enhancers with bidirectional RNA transcription are highly cell type-specific and that disease heritability is strongly enriched in these enhancers. The resulting cell type-resolved multimodal atlas of bidirectionally transcribed enhancers, which we linked with promoters using fine-scale chromatin contact maps, enabled us to systematically interpret genetic variants associated with a range of immune-mediated diseases.

TAD hierarchy restricts poised LTR activation and loss of TAD hierarchy promotes LTR co-option in cancer

Transposable elements (TEs) are abundant in the human genome, and they provide the sources for genetic and functional diversity. The regulation of TEs expression and their functional consequences in physiological conditions and cancer development remain to be fully elucidated. Previous studies suggested TEs are repressed by DNA methylation and chromatin modifications. The effect of 3D chromatin topology on TE regulation remains elusive. Here, by integrating transcriptome and 3D genome architecture studies, we showed that haploinsufficient loss of NIPBL selectively activates alternative promoters at the long terminal repeats (LTRs) of the TE subclasses. This activation occurs through the reorganization of topologically associating domain (TAD) hierarchical structures and recruitment of proximal enhancers. These observations indicate that TAD hierarchy restricts transcriptional activation of LTRs that already possess open chromatin features. In cancer, perturbation of the hierarchical chromatin topology can lead to co-option of LTRs as functional alternative promoters in a context-dependent manner and drive aberrant transcriptional activation of novel oncogenes and other divergent transcripts. These data uncovered a new layer of regulatory mechanism of TE expression beyond DNA and chromatin modification in human genome. They also posit the TAD hierarchy dysregulation as a novel mechanism for alternative promoter-mediated oncogene activation and transcriptional diversity in cancer, which may be exploited therapeutically.

Dissecting the cis-regulatory syntax of transcription initiation with deep learning

Despite extensive characterization of mammalian Pol II transcription, the DNA sequence determinants of transcription initiation at a third of human promoters and most enhancers remain poorly understood. Hence, we trained and interpreted a neural network called ProCapNet that accurately models base-resolution initiation profiles from PRO-cap experiments using local DNA sequence. ProCapNet learns sequence motifs with distinct effects on initiation rates and TSS positioning and uncovers context-specific cryptic initiator elements intertwined within other TF motifs. ProCapNet annotates predictive motifs in nearly all actively transcribed regulatory elements across multiple cell-lines, revealing a shared cis-regulatory logic across promoters and enhancers mediated by a highly epistatic sequence syntax of cooperative and competitive motif interactions. ProCapNet models of RAMPAGE profiles measuring steady-state RNA abundance at TSSs distill initiation signals on par with models trained directly on PRO-cap profiles. ProCapNet learns a largely cell-type-agnostic cis-regulatory code of initiation complementing sequence drivers of cell-type-specific chromatin state critical for accurate prediction of cell-type-specific transcription initiation.

How DNA encodes the start of transcription

A deep-learning model reveals the rules that define transcription initiation.

Cellular HSF1 expression is induced during HIV-1 infection by activation of its promoter mediated through the cooperative interaction of HSF1 and viral Nef protein

The Human Immunodeficiency Virus-1 (HIV-1) tends to activate cellular promoters driving expression of pro-viral genes by complex host-virus interactions for productive infection. We have previously demonstrated that expression of such a positive host factor HSF1 (heat shock factor 1) is elevated during HIV-1 infection; however, the mechanism remains to be elucidated. In the present study, we therefore examined whether HSF1 promoter is induced during HIV-1 infection leading to up-regulation of hsf1 gene expression. We mapped the putative transcription start site (TSS) predicted by Eukaryotic promoter database and deletion constructs of the predicted promoter region were tested through luciferase assay to identify the active promoter. The 347 bp upstream to 153 bp downstream region around the putative TSS displayed the highest activity and both Sp1 (stimulating protein 1) and HSF1 itself were identified to be important for its basal activation. Activity of HSF1 promoter was further stimulated during HIV-1 infection in CD4+ T cells, where interestingly the HSF1-site itself seems to play a major role. In addition, HIV-1 protein Nef (negative factor) was also observed to be responsible for the virus-mediated induction of hsf1 gene expression. Chromatin-immunoprecipitation assays further demonstrate that Nef and HSF1 are co-recruited to the HSF1-binding site and cooperatively act on this promoter. The interplay between host HSF1 and viral Nef on HSF1 promoter eventually leads to increase in HSF1 expression during HIV-1 infection. Understanding the mechanism of HSF1 up-regulation during HIV-1 infection might contribute to future antiviral strategies as HSF1 is a positive regulator of virus replication.

KSHV 3.0: a state-of-the-art annotation of the Kaposi's sarcoma-associated herpesvirus transcriptome using cross-platform sequencing

Kaposi's sarcoma-associated herpesvirus (KSHV) is a large, oncogenic DNA virus belonging to the gammaherpesvirus subfamily. KSHV has been extensively studied with various high-throughput RNA-sequencing approaches to map the transcription start and end sites, the splice junctions, and the translation initiation sites. Despite these efforts, the comprehensive annotation of the viral transcriptome remains incomplete. In the present study, we generated a long-read sequencing data set of the lytic and latent KSHV transcriptome using native RNA and direct cDNA-sequencing methods. This was supplemented with Cap Analysis of Gene Expression sequencing based on a short-read platform. We also utilized data sets from previous publications for our analysis. As a result of this combined approach, we have identified a number of novel viral transcripts and RNA isoforms and have either corroborated or improved the annotation of previously identified viral RNA molecules, thereby notably enhancing our comprehension of the transcriptomic architecture of the KSHV genome. We also evaluated the coding capability of transcripts previously thought to be non-coding by integrating our data on the viral transcripts with translatomic information from other publications.IMPORTANCEDeciphering the viral transcriptome of Kaposi's sarcoma-associated herpesvirus is of great importance because we can gain insight into the molecular mechanism of viral replication and pathogenesis, which can help develop potential targets for antiviral interventions. Specifically, the identification of substantial transcriptional overlaps by this work suggests the existence of a genome-wide interference between transcriptional machineries. This finding indicates the presence of a novel regulatory layer, potentially controlling the expression of viral genes.

Using long-read CAGE sequencing to profile cryptic-promoter-derived transcripts and their contribution to the immunopeptidome

Recent studies have shown that the noncoding genome can produce unannotated proteins as antigens that induce immune response. One major source of this activity is the aberrant epigenetic reactivation of transposable elements (TEs). In tumors, TEs often provide cryptic or alternate promoters, which can generate transcripts that encode tumor-specific unannotated proteins. Thus, TE-derived transcripts (TE transcripts) have the potential to produce tumor-specific, but recurrent, antigens shared among many tumors. Identification of TE-derived tumor antigens holds the promise to improve cancer immunotherapy approaches; however, current genomics and computational tools are not optimized for their detection. Here we combined CAGE technology with full-length long-read transcriptome sequencing (long-read CAGE, or LRCAGE) and developed a suite of computational tools to significantly improve immunopeptidome detection by incorporating TE and other tumor transcripts into the proteome database. By applying our methods to human lung cancer cell line H1299 data, we show that long-read technology significantly improves mapping of promoters with low mappability scores and that LRCAGE guarantees accurate construction of uncharacterized 5' transcript structure. Augmenting a reference proteome database with newly characterized transcripts enabled us to detect noncanonical antigens from HLA-pulldown LC-MS/MS data. Lastly, we show that epigenetic treatment increased the number of noncanonical antigens, particularly those encoded by TE transcripts, which might expand the pool of targetable antigens for cancers with low mutational burden.

Subcytoplasmic location of translation controls protein output

The cytoplasm is highly compartmentalized, but the extent and consequences of subcytoplasmic mRNA localization in non-polarized cells are largely unknown. We determined mRNA enrichment in TIS granules (TGs) and the rough endoplasmic reticulum (ER) through particle sorting and isolated cytosolic mRNAs by digitonin extraction. When focusing on genes that encode non-membrane proteins, we observed that 52% have transcripts enriched in specific compartments. Compartment enrichment correlates with a combinatorial code based on mRNA length, exon length, and 3' UTR-bound RNA-binding proteins. Compartment-biased mRNAs differ in the functional classes of their encoded proteins: TG-enriched mRNAs encode low-abundance proteins with strong enrichment of transcription factors, whereas ER-enriched mRNAs encode large and highly expressed proteins. Compartment localization is an important determinant of mRNA and protein abundance, which is supported by reporter experiments showing that redirecting cytosolic mRNAs to the ER increases their protein expression. In summary, the cytoplasm is functionally compartmentalized by local translation environments.

The regulatory landscape of 5' UTRs in translational control during zebrafish embryogenesis

The 5' UTRs of mRNAs are critical for translation regulation, but their in vivo regulatory features are poorly characterized. Here, we report the regulatory landscape of 5' UTRs during early zebrafish embryogenesis using a massively parallel reporter assay of 18,154 sequences coupled to polysome profiling. We found that the 5' UTR is sufficient to confer temporal dynamics to translation initiation, and identified 86 motifs enriched in 5' UTRs with distinct ribosome recruitment capabilities. A quantitative deep learning model, DaniO5P, revealed a combined role for 5' UTR length, translation initiation site context, upstream AUGs and sequence motifs on in vivo ribosome recruitment. DaniO5P predicts the activities of 5' UTR isoforms and indicates that modulating 5' UTR length and motif grammar contributes to translation initiation dynamics. This study provides a first quantitative model of 5' UTR-based translation regulation in early vertebrate development and lays the foundation for identifying the underlying molecular effectors.

Histone variant HTB4 delays leaf senescence by epigenetic control of Ib bHLH transcription factor-mediated iron homeostasis

Leaf senescence is an orderly process regulated by multiple internal factors and diverse environmental stresses including nutrient deficiency. Histone variants are involved in regulating plant growth and development. However, their functions and underlying regulatory mechanisms in leaf senescence remain largely unclear. Here, we found that H2B histone variant HTB4 functions as a negative regulator of leaf senescence. Loss of function of HTB4 led to early leaf senescence phenotypes that were rescued by functional complementation. RNA-seq analysis revealed that several Ib subgroup basic helix-loop-helix (bHLH) transcription factors (TFs) involved in iron (Fe) homeostasis, including bHLH038, bHLH039, bHLH100, and bHLH101, were suppressed in the htb4 mutant, thereby compromising the expressions of FERRIC REDUCTION OXIDASE 2 (FRO2) and IRON-REGULATED TRANSPORTER (IRT1), two important components of the Fe uptake machinery. Chromatin immunoprecipitation-quantitative polymerase chain reaction analysis revealed that HTB4 could bind to the promoter regions of Ib bHLH TFs and enhance their expression by promoting the enrichment of the active mark H3K4me3 near their transcriptional start sites. Moreover, overexpression of Ib bHLH TFs or IRT1 suppressed the premature senescence phenotype of the htb4 mutant. Our work established a signaling pathway, HTB4-bHLH TFs-FRO2/IRT1-Fe homeostasis, which regulates the onset and progression of leaf senescence.

KSHV 3.0: A State-of-the-Art Annotation of the Kaposi's Sarcoma-Associated Herpesvirus Transcriptome Using Cross-Platform Sequencing

Kaposi's sarcoma-associated herpesvirus (KSHV) is a large, oncogenic DNA virus belonging to the gammaherpesvirus subfamily. KSHV has been extensively studied with various high-throughput RNA-sequencing approaches to map the transcription start and end sites, the splice junctions, and the translation initiation sites. Despite these efforts, the comprehensive annotation of the viral transcriptome remains incomplete. In the present study, we generated a long-read sequencing dataset of the lytic and latent KSHV transcriptome using native RNA and direct cDNA sequencing methods. This was supplemented with CAGE sequencing based on a short-read platform. We also utilized datasets from previous publications for our analysis. As a result of this combined approach, we have identified a number of novel viral transcripts and RNA isoforms and have either corroborated or improved the annotation of previously identified viral RNA molecules, thereby notably enhancing our comprehension of the transcriptomic architecture of the KSHV genome. We also evaluated the coding capability of transcripts previously thought to be non-coding, by integrating our data on the viral transcripts with translatomic information from other publications.

Direct m6A recognition by IMP1 underlays an alternative model of target selection for non-canonical methyl-readers

m6A methylation provides an essential layer of regulation in organismal development, and is aberrant in a range of cancers and neuro-pathologies. The information encoded by m6A methylation is integrated into existing RNA regulatory networks by RNA binding proteins that recognise methylated sites, the m6A readers. m6A readers include a well-characterised class of dedicated proteins, the YTH proteins, as well as a broader group of multi-functional regulators where recognition of m6A is only partially understood. Molecular insight in this recognition is essential to build a mechanistic understanding of global m6A regulation. In this study, we show that the reader IMP1 recognises the m6A using a dedicated hydrophobic platform that assembles on the methyl moiety, creating a stable high-affinity interaction. This recognition is conserved across evolution and independent from the underlying sequence context but is layered upon the strong sequence specificity of IMP1 for GGAC RNA. This leads us to propose a concept for m6A regulation where methylation plays a context-dependent role in the recognition of selected IMP1 targets that is dependent on the cellular concentration of available IMP1, differing from that observed for the YTH proteins.

Comprehensive isoform-level analysis reveals the contribution of alternative isoforms to venom evolution and repertoire diversity

Animal venom systems have emerged as valuable models for investigating how novel polygenic phenotypes may arise from gene evolution by varying molecular mechanisms. However, a significant portion of venom genes produce alternative mRNA isoforms that have not been extensively characterized, hindering a comprehensive understanding of venom biology. In this study, we present a full-length isoform-level profiling workflow integrating multiple RNA sequencing technologies, allowing us to reconstruct a high-resolution transcriptome landscape of venom genes in the parasitoid wasp Pteromalus puparum Our findings demonstrate that more than half of the venom genes generate multiple isoforms within the venom gland. Through mass spectrometry analysis, we confirm that alternative splicing contributes to the diversity of venom proteins, acting as a mechanism for expanding the venom repertoire. Notably, we identified seven venom genes that exhibit distinct isoform usages between the venom gland and other tissues. Furthermore, evolutionary analyses of venom serpin3 and orcokinin further reveal that the co-option of an ancient isoform and a newly evolved isoform, respectively, contributes to venom recruitment, providing valuable insights into the genetic mechanisms driving venom evolution in parasitoid wasps. Overall, our study presents a comprehensive investigation of venom genes at the isoform level, significantly advancing our understanding of alternative isoforms in venom diversity and evolution and setting the stage for further in-depth research on venoms.

Improving the annotation of the cattle genome by annotating transcription start sites in a diverse set of tissues and populations using Cap Analysis Gene Expression sequencing

Understanding the genomic control of tissue-specific gene expression and regulation can help to inform the application of genomic technologies in farm animal breeding programs. The fine mapping of promoters [transcription start sites (TSS)] and enhancers (divergent amplifying segments of the genome local to TSS) in different populations of cattle across a wide diversity of tissues provides information to locate and understand the genomic drivers of breed- and tissue-specific characteristics. To this aim, we used Cap Analysis Gene Expression (CAGE) sequencing, of 24 different tissues from 3 populations of cattle, to define TSS and their coexpressed short-range enhancers (<1 kb) in the ARS-UCD1.2_Btau5.0.1Y reference genome (1000bulls run9) and analyzed tissue and population specificity of expressed promoters. We identified 51,295 TSS and 2,328 TSS-Enhancer regions shared across the 3 populations (dairy, beef-dairy cross, and Canadian Kinsella composite cattle from 2 individuals, 1 of each sex, per population). Cross-species comparative analysis of CAGE data from 7 other species, including sheep, revealed a set of TSS and TSS-Enhancers that were specific to cattle. The CAGE data set will be combined with other transcriptomic information for the same tissues to create a new high-resolution map of transcript diversity across tissues and populations in cattle for the BovReg project. Here we provide the CAGE data set and annotation tracks for TSS and TSS-Enhancers in the cattle genome. This new annotation information will improve our understanding of the drivers of gene expression and regulation in cattle and help to inform the application of genomic technologies in breeding programs.

nf-core/clipseq - a robust Nextflow pipeline for comprehensive CLIP data analysis

Crosslinking and immunoprecipitation (CLIP) technologies have become a central component of the molecular biologists' toolkit to study protein-RNA interactions and thus to uncover core principles of RNA biology. There has been a proliferation of CLIP-based experimental protocols, as well as computational tools, especially for peak-calling. Consequently, there is an urgent need for a well-documented bioinformatic pipeline that enshrines the principles of robustness, reproducibility, scalability, portability and flexibility while embracing the diversity of experimental and computational CLIP tools. To address this, we present nf-core/clipseq - a robust Nextflow pipeline for quality control and analysis of CLIP sequencing data. It is part of the international nf-core community effort to develop and curate a best-practice, gold-standard set of pipelines for data analysis. The standards enabled by Nextflow and nf-core, including workflow management, version control, continuous integration and containerisation ensure that these key needs are met. Furthermore, multiple tools are implemented ( e.g. for peak-calling), alongside visualisation of quality control metrics to empower the user to make their own informed decisions based on their data. nf-core/clipseq remains under active development, with plans to incorporate newly released tools to ensure that pipeline remains up-to-date and relevant for the community. Engagement with users and developers is encouraged through the nf-core GitHub repository and Slack channel to promote collaboration. It is available at https://nf-co.re/clipseq.

FIPRESCI: droplet microfluidics based combinatorial indexing for massive-scale 5'-end single-cell RNA sequencing

Single-cell RNA sequencing methods focusing on the 5'-end of transcripts can reveal promoter and enhancer activity and efficiently profile immune receptor repertoire. However, ultra-high-throughput 5'-end single-cell RNA sequencing methods have not been described. We introduce FIPRESCI, 5'-end single-cell combinatorial indexing RNA-Seq, enabling massive sample multiplexing and increasing the throughput of the droplet microfluidics system by over tenfold. We demonstrate FIPRESCI enables the generation of approximately 100,000 single-cell transcriptomes from E10.5 whole mouse embryos in a single-channel experiment, and simultaneous identification of subpopulation differences and T cell receptor signatures of peripheral blood T cells from 12 cancer patients.

Strand Asymmetries Across Genomic Processes

Across biological systems, a number of genomic processes, including transcription, replication, DNA repair, and transcription factor binding, display intrinsic directionalities. These directionalities are reflected in the asymmetric distribution of nucleotides, motifs, genes, transposon integration sites, and other functional elements across the two complementary strands. Strand asymmetries, including GC skews and mutational biases, have shaped the nucleotide composition of diverse organisms. The investigation of strand asymmetries often serves as a method to understand underlying biological mechanisms, including protein binding preferences, transcription factor interactions, retrotransposition, DNA damage and repair preferences, transcription-replication collisions, and mutagenesis mechanisms. Research into this subject also enables the identification of functional genomic sites, such as replication origins and transcription start sites. Improvements in our ability to detect and quantify DNA strand asymmetries will provide insights into diverse functionalities of the genome, the contribution of different mutational mechanisms in germline and somatic mutagenesis, and our knowledge of genome instability and evolution, which all have significant clinical implications in human disease, including cancer. In this review, we describe key developments that have been made across the field of genomic strand asymmetries, as well as the discovery of associated mechanisms.

Predictive model of transcriptional elongation control identifies trans regulatory factors from chromatin signatures

Promoter-proximal Polymerase II (Pol II) pausing is a key rate-limiting step for gene expression. DNA and RNA-binding trans-acting factors regulating the extent of pausing have been identified. However, we lack a quantitative model of how interactions of these factors determine pausing, therefore the relative importance of implicated factors is unknown. Moreover, previously unknown regulators might exist. Here we address this gap with a machine learning model that accurately predicts the extent of promoter-proximal Pol II pausing from large-scale genome and transcriptome binding maps and gene annotation and sequence composition features. We demonstrate high accuracy and generalizability of the model by validation on an independent cell line which reveals the model's cell line agnostic character. Model interpretation in light of prior knowledge about molecular functions of regulatory factors confirms the interconnection of pausing with other RNA processing steps. Harnessing underlying feature contributions, we assess the relative importance of each factor, quantify their predictive effects and systematically identify previously unknown regulators of pausing. We additionally identify 16 previously unknown 7SK ncRNA interacting RNA-binding proteins predictive of pausing. Our work provides a framework to further our understanding of the regulation of the critical early steps in transcriptional elongation.

Cisplatin resistance driver claspin is a target for immunotherapy in urothelial carcinoma

Bladder cancer is a major and fatal urological disease. Cisplatin is a key drug for the treatment of bladder cancer, especially in muscle-invasive cases. In most cases of bladder cancer, cisplatin is effective; however, resistance to cisplatin has a significant negative impact on prognosis. Thus, a treatment strategy for cisplatin-resistant bladder cancer is essential to improve the prognosis. In this study, we established a cisplatin-resistant (CR) bladder cancer cell line using an urothelial carcinoma cell lines (UM-UC-3 and J82). We screened for potential targets in CR cells and found that claspin (CLSPN) was overexpressed. CLSPN mRNA knockdown revealed that CLSPN had a role in cisplatin resistance in CR cells. In our previous study, we identified human leukocyte antigen (HLA)-A*02:01-restricted CLSPN peptide by HLA ligandome analysis. Thus, we generated a CLSPN peptide-specific cytotoxic T lymphocyte clone that recognized CR cells at a higher level than wild-type UM-UC-3 cells. These findings indicate that CLSPN is a driver of cisplatin resistance and CLSPN peptide-specific immunotherapy may be effective for cisplatin-resistant cases.

ATM suppresses c-Myc overexpression in the mammary epithelium in response to estrogen

ATM gene mutation carriers are predisposed to estrogen-receptor-positive breast cancer (BC). ATM prevents BC oncogenesis by activating p53 in every cell; however, much remains unknown about tissue-specific oncogenesis after ATM loss. Here, we report that ATM controls the early transcriptional response to estrogens. This response depends on topoisomerase II (TOP2), which generates TOP2-DNA double-strand break (DSB) complexes and rejoins the breaks. When TOP2-mediated ligation fails, ATM facilitates DSB repair. After estrogen exposure, TOP2-dependent DSBs arise at the c-MYC enhancer in human BC cells, and their defective repair changes the activation profile of enhancers and induces the overexpression of many genes, including the c-MYC oncogene. CRISPR/Cas9 cleavage at the enhancer also causes c-MYC overexpression, indicating that this DSB causes c-MYC overexpression. Estrogen treatment induced c-Myc protein overexpression in mammary epithelial cells of ATM-deficient mice. In conclusion, ATM suppresses the c-Myc-driven proliferative effects of estrogens, possibly explaining such tissue-specific oncogenesis.

Class II LitR serves as an effector of "short" LOV-type blue-light photoreceptor in Pseudomonas mendocina

PmlR2, a class II LitR/CarH family transcriptional regulator, and PmSB-LOV, a "short" LOV-type blue light photoreceptor, are adjacently encoded in Pseudomonas mendocina NBRC 14162. An effector protein for the "short" LOV-type photoreceptor in Pseudomonas has not yet been identified. Here, we show that PmlR2 is an effector protein of PmSB-LOV. Transcriptional analyses revealed that the expression of genes located near pmlR2 and its homolog gene, pmlR1, was induced in response to illumination. In vitro DNA-protein binding analyses showed that recombinant PmlR2 directly binds to the promoter region of light-inducible genes. Furthermore PmSB-LOV exhibited a typical LOV-type light-induced spectral change. Gel-filtration chromatography demonstrated that the illuminated PmSB-LOV was directly associated with PmlR2, whereas non-illuminated proteins did not interact. The inhibition of PmlR2 function following PmSB-LOV binding was verified by surface plasmon resonance: the DNA-binding ability of PmlR2 was specifically inhibited in the presence of blue light-illuminated-PmSB-LOV. An In vitro transcription assay showed a dose-dependent reduction in PmlR2 repressor activity in the presence of illuminated PmSB-LOV. Overall, evidence suggests that the DNA-binding activity of PmlR2 is inhibited by its direct association with blue light-activated PmSB-LOV, enabling transcription of light-inducible promoters by RNA polymerase.

DeepTSS: multi-branch convolutional neural network for transcription start site identification from CAGE data

BACKGROUND

The widespread usage of Cap Analysis of Gene Expression (CAGE) has led to numerous breakthroughs in understanding the transcription mechanisms. Recent evidence in the literature, however, suggests that CAGE suffers from transcriptional and technical noise. Regardless of the sample quality, there is a significant number of CAGE peaks that are not associated with transcription initiation events. This type of signal is typically attributed to technical noise and more frequently to random five-prime capping or transcription bioproducts. Thus, the need for computational methods emerges, that can accurately increase the signal-to-noise ratio in CAGE data, resulting in error-free transcription start site (TSS) annotation and quantification of regulatory region usage. In this study, we present DeepTSS, a novel computational method for processing CAGE samples, that combines genomic signal processing (GSP), structural DNA features, evolutionary conservation evidence and raw DNA sequence with Deep Learning (DL) to provide single-nucleotide TSS predictions with unprecedented levels of performance.

RESULTS

To evaluate DeepTSS, we utilized experimental data, protein-coding gene annotations and computationally-derived genome segmentations by chromatin states. DeepTSS was found to outperform existing algorithms on all benchmarks, achieving 98% precision and 96% sensitivity (accuracy 95.4%) on the protein-coding gene strategy, with 96.66% of its positive predictions overlapping active chromatin, 98.27% and 92.04% co-localized with at least one transcription factor and H3K4me3 peak.

CONCLUSIONS

CAGE is a key protocol in deciphering the language of transcription, however, as every experimental protocol, it suffers from biological and technical noise that can severely affect downstream analyses. DeepTSS is a novel DL-based method for effectively removing noisy CAGE signal. In contrast to existing software, DeepTSS does not require feature selection since the embedded convolutional layers can readily identify patterns and only utilize the important ones for the classification task. This study highlights the key role that DL can play in Molecular Biology, by removing the inherent flaws of experimental protocols, that form the backbone of contemporary research. Here, we show how DeepTSS can unleash the full potential of an already popular and mature method such as CAGE, and push the boundaries of coding and non-coding gene expression regulator research even further.

mRNA 5' terminal sequences drive 200-fold differences in expression through effects on synthesis, translation and decay

mRNA regulatory sequences control gene expression at multiple levels including translation initiation and mRNA decay. The 5' terminal sequences of mRNAs have unique regulatory potential because of their proximity to key post-transcriptional regulators. Here we have systematically probed the function of 5' terminal sequences in gene expression in human cells. Using a library of reporter mRNAs initiating with all possible 7-mer sequences at their 5' ends, we find an unexpected impact on transcription that underlies 200-fold differences in mRNA expression. Library sequences that promote high levels of transcription mirrored those found in native mRNAs and define two basic classes with similarities to classic Initiator (Inr) and TCT core promoter motifs. By comparing transcription, translation and decay rates, we identify sequences that are optimized for both efficient transcription and growth-regulated translation and stability, including variants of terminal oligopyrimidine (TOP) motifs. We further show that 5' sequences of endogenous mRNAs are enriched for multi-functional TCT/TOP hybrid sequences. Together, our results reveal how 5' sequences define two general classes of mRNAs with distinct growth-responsive profiles of expression across synthesis, translation and decay.

Promoter sequence and architecture determine expression variability and confer robustness to genetic variants

Genetic and environmental exposures cause variability in gene expression. Although most genes are affected in a population, their effect sizes vary greatly, indicating the existence of regulatory mechanisms that could amplify or attenuate expression variability. Here, we investigate the relationship between the sequence and transcription start site architectures of promoters and their expression variability across human individuals. We find that expression variability can be largely explained by a promoter's DNA sequence and its binding sites for specific transcription factors. We show that promoter expression variability reflects the biological process of a gene, demonstrating a selective trade-off between stability for metabolic genes and plasticity for responsive genes and those involved in signaling. Promoters with a rigid transcription start site architecture are more prone to have variable expression and to be associated with genetic variants with large effect sizes, while a flexible usage of transcription start sites within a promoter attenuates expression variability and limits genotypic effects. Our work provides insights into the variable nature of responsive genes and reveals a novel mechanism for supplying transcriptional and mutational robustness to essential genes through multiple transcription start site regions within a promoter.

SCAFE: a software suite for analysis of transcribed cis-regulatory elements in single cells

MOTIVATION

Cell type-specific activities of cis-regulatory elements (CRE) are central to understanding gene regulation and disease predisposition. Single-cell RNA 5'end sequencing (sc-end5-seq) captures the transcription start sites (TSS) which can be used as a proxy to measure the activity of transcribed CREs (tCREs). However, a substantial fraction of TSS identified from sc-end5-seq data may not be genuine due to various artifacts, hindering the use of sc-end5-seq for de novo discovery of tCREs.

RESULTS

We developed SCAFE-Single-Cell Analysis of Five-prime Ends-a software suite that processes sc-end5-seq data to de novo identify TSS clusters based on multiple logistic regression. It annotates tCREs based on the identified TSS clusters and generates a tCRE-by-cell count matrix for downstream analyses. The software suite consists of a set of flexible tools that could either be run independently or as pre-configured workflows.

AVAILABILITY AND IMPLEMENTATION

SCAFE is implemented in Perl and R. The source code and documentation are freely available for download under the MIT License from https://github.com/chung-lab/SCAFE. Docker images are available from https://hub.docker.com/r/cchon/scafe. The submitted software version and test data are archived at https://doi.org/10.5281/zenodo.7023163 and https://doi.org/10.5281/zenodo.7024060, respectively.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Deep sequencing of short capped RNAs reveals novel families of noncoding RNAs

In eukaryotes, capped RNAs include long transcripts such as messenger RNAs and long noncoding RNAs, as well as shorter transcripts such as spliceosomal RNAs, small nucleolar RNAs, and enhancer RNAs. Long capped transcripts can be profiled using cap analysis gene expression (CAGE) sequencing and other methods. Here, we describe a sequencing library preparation protocol for short capped RNAs, apply it to a differentiation time course of the human cell line THP-1, and systematically compare the landscape of short capped RNAs to that of long capped RNAs. Transcription initiation peaks associated with genes in the sense direction have a strong preference to produce either long or short capped RNAs, with one out of six peaks detected in the short capped RNA libraries only. Gene-associated short capped RNAs have highly specific 3' ends, typically overlapping splice sites. Enhancers also preferentially generate either short or long capped RNAs, with 10% of enhancers observed in the short capped RNA libraries only. Enhancers producing either short or long capped RNAs show enrichment for GWAS-associated disease SNPs. We conclude that deep sequencing of short capped RNAs reveals new families of noncoding RNAs and elucidates the diversity of transcripts generated at known and novel promoters and enhancers.

Isoform-resolved mRNA profiling of ribosome load defines interplay of HIF and mTOR dysregulation in kidney cancer

Hypoxia inducible factor (HIF) and mammalian target of rapamycin (mTOR) pathways orchestrate responses to oxygen and nutrient availability. These pathways are frequently dysregulated in cancer, but their interplay is poorly understood, in part because of difficulties in simultaneous measurement of global and mRNA-specific translation. Here, we describe a workflow for measurement of ribosome load of mRNAs resolved by their transcription start sites (TSSs). Its application to kidney cancer cells reveals extensive translational reprogramming by mTOR, strongly affecting many metabolic enzymes and pathways. By contrast, global effects of HIF on translation are limited, and we do not observe reported translational activation by HIF2A. In contrast, HIF-dependent alterations in TSS usage are associated with robust changes in translational efficiency in a subset of genes. Analyses of the interplay of HIF and mTOR reveal that specific classes of HIF1A and HIF2A transcriptional target gene manifest different sensitivity to mTOR, in a manner that supports combined use of HIF2A and mTOR inhibitors in treatment of kidney cancer.

Transcriptomic diversity in human medullary thymic epithelial cells

The induction of central T cell tolerance in the thymus depends on the presentation of peripheral self-epitopes by medullary thymic epithelial cells (mTECs). This promiscuous gene expression (pGE) drives mTEC transcriptomic diversity, with non-canonical transcript initiation, alternative splicing, and expression of endogenous retroelements (EREs) representing important but incompletely understood contributors. Here we map the expression of genome-wide transcripts in immature and mature human mTECs using high-throughput 5' cap and RNA sequencing. Both mTEC populations show high splicing entropy, potentially driven by the expression of peripheral splicing factors. During mTEC maturation, rates of global transcript mis-initiation increase and EREs enriched in long terminal repeat retrotransposons are up-regulated, the latter often found in proximity to differentially expressed genes. As a resource, we provide an interactive public interface for exploring mTEC transcriptomic diversity. Our findings therefore help construct a map of transcriptomic diversity in the healthy human thymus and may ultimately facilitate the identification of those epitopes which contribute to autoimmunity and immune recognition of tumor antigens.

Reactivity-dependent profiling of RNA 5-methylcytidine dioxygenases

Epitranscriptomic RNA modifications can regulate fundamental biological processes, but we lack approaches to map modification sites and probe writer enzymes. Here we present a chemoproteomic strategy to characterize RNA 5-methylcytidine (m⁵C) dioxygenase enzymes in their native context based upon metabolic labeling and activity-based crosslinking with 5-ethynylcytidine (5-EC). We profile m⁵C dioxygenases in human cells including ALKBH1 and TET2 and show that ALKBH1 is the major hm⁵C- and f⁵C-forming enzyme in RNA. Further, we map ALKBH1 modification sites transcriptome-wide using 5-EC-iCLIP and ARP-based sequencing to identify ALKBH1-dependent m⁵C oxidation in a variety of tRNAs and mRNAs and analyze ALKBH1 substrate specificity in vitro. We also apply targeted pyridine borane-mediated sequencing to measure f⁵C sites on select tRNA. Finally, we show that f⁵C at the wobble position of tRNA-Leu-CAA plays a role in decoding Leu codons under stress. Our work provides powerful chemical approaches for studying RNA m⁵C dioxygenases and mapping oxidative m⁵C modifications and reveals the existence of novel epitranscriptomic pathways for regulating RNA function.

Kleiner and co-workers profile RNA 5-methylcytidine (m⁵C) dioxygenase enzymes using an activity-based metabolic probing strategy. They reveal ALKBH1 as the major 5-formylcytidine (f⁵C) writer and characterize modification sites across mRNA and tRNA.

Sequence-based evaluation of promoter context for prediction of transcription start sites in Arabidopsis and rice

Genes are transcribed from transcription start sites (TSSs), and their position in a genome is strictly controlled to avoid mis-expression of undesired regions. In this study, we designed and developed a methodology for the evaluation of promoter context, which detects proximal promoter regions from − 200 to − 60 bp relative to a TSS, in Arabidopsis and rice genomes. The method positively evaluates spacer sequences and Regulatory Element Groups, but not core promoter elements like TATA boxes, and is able to predict the position of a TSS within a width of 200 bp. An important feature of the evaluation/prediction method is its independence of the core promoter elements, which was demonstrated by successful prediction of all the TATA, GA, and coreless types of promoters without notable differences in the accuracy of prediction. The positive relationship identified between the evaluation scores and gene expression levels suggests that this method is useful for the evaluation of promoter maturity.

A comparison of experimental assays and analytical methods for genome-wide identification of active enhancers

Mounting evidence supports the idea that transcriptional patterns serve as more specific identifiers of active enhancers than histone marks; however, the optimal strategy to identify active enhancers both experimentally and computationally has not been determined. Here, we compared 13 genome-wide RNA sequencing assays in K562 cells and showed that the nuclear run-on followed by cap-selection assay (GRO/PRO-cap) has advantages in eRNA detection and active enhancer identification. We also introduced a tool, Peak Identifier for Nascent Transcript Starts (PINTS), to identify active promoters and enhancers genome-wide and pinpoint the precise location of the 5′ transcription start sites. Finally, we compiled a comprehensive enhancer candidate compendium based on the detected eRNA TSSs available in 120 cell and tissue types that can be accessed at https://pints.yulab.org. With the knowledge of the best available assays and pipelines, this large-scale annotation of candidate enhancers will pave the way for selection and characterization of their functions in a time- and labor-efficient manner in the future.

Suppression of human trophoblast syncytialization by human cytomegalovirus infection

INTRODUCTION

Placental dysfunction triggers fetal growth restriction in congenital human cytomegalovirus (HCMV) infection. Studies suggest that HCMV infection interferes with the differentiation of human trophoblasts. However, the underlying mechanisms have not been clarified. This study investigated the impact of HCMV infection on gene transcriptomes in cytotrophoblasts (CTBs) associated with placental dysfunction.

METHODS

CTBs were isolated from human term placentas, and spontaneous syncytialization was observed in vitro. The transcriptome profiles were compared between CTB groups with and without HCMV infection by cap analysis gene expression sequencing. The effect of HCMV infection on trophoblast differentiation was evaluated by examining cell fusion status using immunocytochemical staining for desmoplakin and assessing the production of cell differentiation markers, including hCG, PlGF, and soluble Flt-1, using ELISA.

RESULTS

The expression of the genes categorized in the signaling pathways related to the cell cycle was significantly enhanced in CTBs with HCMV infection compared with uninfected CTBs. HCMV infection hindered the alteration of the gene expression profile associated with syncytialization. This suppressive effect under HCMV infection was concurrent with the reduction in hCG and PlGF secretion. Immunostaining for desmoplakin revealed that HCMV infection reduced the cell fusion of cultured CTBs. These findings imply that HCMV infection has a negative impact on syncytialization, which is indispensable for the maintenance of villous function.

DISCUSSION

HCMV infection interferes with gene expression profiles and functional differentiation of trophoblasts. Suppression of syncytialization may be a survival strategy for HCMV to expand infection and could be associated with placental dysfunction.

TSSr: an R package for comprehensive analyses of TSS sequencing data

Transcription initiation is regulated in a highly organized fashion to ensure proper cellular functions. Accurate identification of transcription start sites (TSSs) and quantitative characterization of transcription initiation activities are fundamental steps for studies of regulated transcriptions and core promoter structures. Several high-throughput techniques have been developed to sequence the very 5'end of RNA transcripts (TSS sequencing) on the genome scale. Bioinformatics tools are essential for processing, analysis, and visualization of TSS sequencing data. Here, we present TSSr, an R package that provides rich functions for mapping TSS and characterizations of structures and activities of core promoters based on all types of TSS sequencing data. Specifically, TSSr implements several newly developed algorithms for accurately identifying TSSs from mapped sequencing reads and inference of core promoters, which are a prerequisite for subsequent functional analyses of TSS data. Furthermore, TSSr also enables users to export various types of TSS data that can be visualized by genome browser for inspection of promoter activities in association with other genomic features, and to generate publication-ready TSS graphs. These user-friendly features could greatly facilitate studies of transcription initiation based on TSS sequencing data. The source code and detailed documentations of TSSr can be freely accessed at https://github.com/Linlab-slu/TSSr.

Activity-based RNA-modifying enzyme probing reveals DUS3L-mediated dihydrouridylation

Epitranscriptomic RNA modifications can regulate RNA activity; however, there remains a major gap in our understanding of the RNA chemistry present in biological systems. Here we develop RNA-mediated activity-based protein profiling (RNABPP), a chemoproteomic strategy that relies on metabolic RNA labeling, mRNA interactome capture and quantitative proteomics, to investigate RNA-modifying enzymes in human cells. RNABPP with 5-fluoropyrimidines allowed us to profile 5-methylcytidine (m⁵C) and 5-methyluridine (m⁵U) methyltransferases. Further, we uncover a new mechanism-based crosslink between 5-fluorouridine (5-FUrd)-modified RNA and the dihydrouridine synthase (DUS) homolog DUS3L. We investigate the mechanism of crosslinking and use quantitative nucleoside liquid chromatography-tandem mass spectrometry (LC-MS/MS) analysis and 5-FUrd-based crosslinking and immunoprecipitation (CLIP) sequencing to map DUS3L-dependent dihydrouridine (DHU) modifications across the transcriptome. Finally, we show that DUS3L-knockout (KO) cells have compromised protein translation rates and impaired cellular proliferation. Taken together, our work provides a general approach for profiling RNA-modifying enzyme activity in living cells and reveals new pathways for epitranscriptomic RNA regulation.

Global approaches for profiling transcription initiation

Transcription start site (TSS) selection influences transcript stability and translation as well as protein sequence. Alternative TSS usage is pervasive in organismal development, is a major contributor to transcript isoform diversity in humans, and is frequently observed in human diseases including cancer. In this review, we discuss the breadth of techniques that have been used to globally profile TSSs and the resulting insights into gene regulation, as well as future prospects in this area of inquiry.

TDP-43 condensation properties specify its RNA-binding and regulatory repertoire

Mutations causing amyotrophic lateral sclerosis (ALS) often affect the condensation properties of RNA-binding proteins (RBPs). However, the role of RBP condensation in the specificity and function of protein-RNA complexes remains unclear. We created a series of TDP-43 C-terminal domain (CTD) variants that exhibited a gradient of low to high condensation propensity, as observed in vitro and by nuclear mobility and foci formation. Notably, a capacity for condensation was required for efficient TDP-43 assembly on subsets of RNA-binding regions, which contain unusually long clusters of motifs of characteristic types and density. These "binding-region condensates" are promoted by homomeric CTD-driven interactions and required for efficient regulation of a subset of bound transcripts, including autoregulation of TDP-43 mRNA. We establish that RBP condensation can occur in a binding-region-specific manner to selectively modulate transcriptome-wide RNA regulation, which has implications for remodeling RNA networks in the context of signaling, disease, and evolution.

GRIK2 is a target for bladder cancer stem-like cell-targeting immunotherapy

Recent studies have revealed that treatment-resistant cancer stem-like cells (CSCs)/cancer-initiating cells (CICs) can be targeted by cytotoxic T lymphocytes (CTLs). CTLs recognize antigenic peptides derived from tumor-associated antigens; thus, the identification of tumor-associated antigens expressed by CSCs/CICs is essential. Human leucocyte antigen (HLA) ligandome analysis using mass spectrometry enables the analysis of naturally expressed antigenic peptides; however, HLA ligandome analysis requires a large number of cells and is challenging for CSCs/CICs. In this study, we established a novel bladder CSC/CIC model from a bladder cancer cell line (UM-UC-3 cells) using an ALDEFLUOR assay. CSCs/CICs were isolated as aldehyde dehydrogenase (ALDH)-high cells and several ALDH^high clone cells were established. ALDH^high clone cells were enriched with CSCs/CICs by sphere formation and tumorigenicity in immunodeficient mice. HLA ligandome analysis and cap analysis of gene expression using ALDH^high clone cells revealed a distinctive antigenic peptide repertoire in bladder CSCs/CICs, and we found that a glutamate receptor, ionotropic, kainite 2 (GRIK2)-derived antigenic peptide (LMYDAVHVV) was specifically expressed by CSCs/CICs. A GRIK2 peptide-specific CTL clone recognized GRIK2-overexpressing UM-UC-3 cells and ALDH^high clone cells, indicating that GRIK2 peptide can be a novel target for bladder CSC/CIC-targeting immunotherapy.

Computational identification and experimental characterization of preferred downstream positions in human core promoters

Metazoan core promoters, which direct the initiation of transcription by RNA polymerase II (Pol II), may contain short sequence motifs termed core promoter elements/motifs (e.g. the TATA box, initiator (Inr) and downstream core promoter element (DPE)), which recruit Pol II via the general transcription machinery. The DPE was discovered and extensively characterized in Drosophila, where it is strictly dependent on both the presence of an Inr and the precise spacing from it. Since the Drosophila DPE is recognized by the human transcription machinery, it is most likely that some human promoters contain a downstream element that is similar, though not necessarily identical, to the Drosophila DPE. However, only a couple of human promoters were shown to contain a functional DPE, and attempts to computationally detect human DPE-containing promoters have mostly been unsuccessful. Using a newly-designed motif discovery strategy based on Expectation-Maximization probabilistic partitioning algorithms, we discovered preferred downstream positions (PDP) in human promoters that resemble the Drosophila DPE. Available chromatin accessibility footprints revealed that Drosophila and human Inr+DPE promoter classes are not only highly structured, but also similar to each other, particularly in the proximal downstream region. Clustering of the corresponding sequence motifs using a neighbor-joining algorithm strongly suggests that canonical Inr+DPE promoters could be common to metazoan species. Using reporter assays we demonstrate the contribution of the identified downstream positions to the function of multiple human promoters. Furthermore, we show that alteration of the spacing between the Inr and PDP by two nucleotides results in reduced promoter activity, suggesting a spacing dependency of the newly discovered human PDP on the Inr. Taken together, our strategy identified novel functional downstream positions within human core promoters, supporting the existence of DPE-like motifs in human promoters.

Transcription of genes by the RNA polymerase II enzyme initiates at a genomic region termed the core promoter. The core promoter is a regulatory region that may contain diverse short DNA sequence motifs/elements that confer specific properties to it. Interestingly, core promoter motifs can be located both upstream and downstream of the transcription start site. Variable compositions of core promoter elements were identified. The initiator (Inr) motif and the downstream core promoter element (DPE) is a combination of elements that has been identified and extensively characterized in fruit flies. Although a few Inr+DPE -containing human promoters were identified, the presence of transcriptionally important downstream core promoter positions within human promoters has been a matter of controversy in the literature. Here, using a newly-designed motif discovery strategy, we discovered preferred downstream positions in human promoters that resemble fruit fly DPE. Clustering of the corresponding sequence motifs in eight additional species indicated that such promoters could be common to multicellular non-plant organisms. Importantly, functional characterization of the newly discovered preferred downstream positions supports the existence of Inr+DPE-containing promoters in human genes.

Single-cell alternative polyadenylation analysis delineates GABAergic neuron types

Background

Alternative polyadenylation (APA) is emerging as an important mechanism in the post-transcriptional regulation of gene expression across eukaryotic species. Recent studies have shown that APA plays key roles in biological processes, such as cell proliferation and differentiation. Single-cell RNA-seq technologies are widely used in gene expression heterogeneity studies; however, systematic studies of APA at the single-cell level are still lacking.

Results

Here, we described a novel computational framework, SAPAS, that utilizes 3′-tag-based scRNA-seq data to identify novel poly(A) sites and quantify APA at the single-cell level. Applying SAPAS to the scRNA-seq data of phenotype characterized GABAergic interneurons, we identified cell type-specific APA events for different GABAergic neuron types. Genes with cell type-specific APA events are enriched for synaptic architecture and communications. In further, we observed a strong enrichment of heritability for several psychiatric disorders and brain traits in altered 3′ UTRs and coding sequences of cell type-specific APA events. Finally, by exploring the modalities of APA, we discovered that the bimodal APA pattern of Pak3 could classify chandelier cells into different subpopulations that are from different laminar positions.

Conclusions

We established a method to characterize APA at the single-cell level. When applied to a scRNA-seq dataset of GABAergic interneurons, the single-cell APA analysis not only identified cell type-specific APA events but also revealed that the modality of APA could classify cell subpopulations. Thus, SAPAS will expand our understanding of cellular heterogeneity.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01076-3.

Redefining the PTEN promoter: Identification of novel upstream transcription start regions

Germline mutation of PTEN is causally observed in Cowden syndrome (CS) and is one of the most common, penetrant risk genes for autism spectrum disorder (ASD). However, the majority of individuals who present with CS-like clinical features are PTEN-mutation negative. Reassessment of PTEN promoter regulation may help explain abnormal PTEN dosage, as only the minimal promoter and coding regions are currently included in diagnostic PTEN mutation analysis. Therefore, we reanalyzed the architecture of the PTEN promoter using next-generation sequencing datasets. Specifically, run-on sequencing assays identified two additional transcription start regions (TSRs) at -2053 and - 1906 basepairs from the canonical start of PTEN, thus extending the PTEN 5'UTR and redefining the PTEN promoter. We show that these novel upstream TSRs are active in cancer cell lines, human cancer, and normal tissue. Further, these TSRs can produce novel PTEN transcripts due to the introduction of new splice donors at -2041, -1826, and - 1355, which may allow for splicing out of the PTEN 5'UTR or the first and second exon in upstream-initiated transcripts. Combining ENCODE ChIP-seq and pertinent literature, we also compile and analyze all transcription factors (TFs) binding at the redefined PTEN locus. Enrichment analyses suggest that TFs bind specifically to the upstream TSRs may be implicated in inflammatory processes. Together, these data redefine the architecture of the PTEN promoter, an important step toward a comprehensive model of PTEN transcription regulation, a basis for future investigations into the new promoters' role in disease pathogenesis.

A Transcription Start Site Map in Human Pancreatic Islets Reveals Functional Regulatory Signatures

Identifying the tissue-specific molecular signatures of active regulatory elements is critical to understand gene regulatory mechanisms. Here, we identify transcription start sites (TSS) using cap analysis of gene expression (CAGE) across 57 human pancreatic islet samples. We identify 9,954 reproducible CAGE tag clusters (TCs), ∼20% of which are islet specific and occur mostly distal to known gene TSS. We integrated islet CAGE data with histone modification and chromatin accessibility profiles to identify epigenomic signatures of transcription initiation. Using a massively parallel reporter assay, we validated the transcriptional enhancer activity for 2,279 of 3,378 (∼68%) tested islet CAGE elements (5% false discovery rate). TCs within accessible enhancers show higher enrichment to overlap type 2 diabetes genome-wide association study (GWAS) signals than existing islet annotations, which emphasizes the utility of mapping CAGE profiles in disease-relevant tissue. This work provides a high-resolution map of transcriptional initiation in human pancreatic islets with utility for dissecting active enhancers at GWAS loci.

Human White Adipose Tissue Displays Selective Insulin Resistance in the Obese State

Selective hepatic insulin resistance is a feature of obesity and type 2 diabetes. Whether similar mechanisms operate in white adipose tissue (WAT) of those with obesity and to what extent these are normalized by weight loss are unknown. We determined insulin sensitivity by hyperinsulinemic euglycemic clamp and insulin response in subcutaneous WAT by RNA sequencing in 23 women with obesity before and 2 years after bariatric surgery. To control for effects of surgery, women postsurgery were matched to never-obese women. Multidimensional analyses of 138 samples allowed us to classify the effects of insulin into three distinct expression responses: a common set was present in all three groups and included genes encoding several lipid/cholesterol biosynthesis enzymes; a set of obesity-attenuated genes linked to tissue remodeling and protein translation was selectively regulated in the two nonobese states; and several postobesity-enriched genes encoding proteins involved in, for example, one-carbon metabolism were only responsive to insulin in the women who had lost weight. Altogether, human WAT displays a selective insulin response in the obese state, where most genes are normalized by weight loss. This comprehensive atlas provides insights into the transcriptional effects of insulin in WAT and may identify targets to improve insulin action.

WWOX Loses the Ability to Regulate Oncogenic AP-2γ and Synergizes with Tumor Suppressor AP-2α in High-Grade Bladder Cancer

The cytogenic locus of the WWOX gene overlaps with the second most active fragile site, FRA16D, which is present at a higher frequency in bladder cancer (BLCA) patients with smoking habit, a known risk factor of this tumor. Recently, we demonstrated the relevance of the role of WWOX in grade 2 BLCA in collaboration with two AP-2 transcription factors whose molecular actions supported or opposed pro-cancerous events, suggesting a distinct character. As further research is needed on higher grades, the aim of the present study was to examine WWOX-AP-2 functionality in grade 3 and 4 BLCA using equivalent in vitro methodology with additional transcriptome profiling of cellular variants. WWOX and AP-2α demonstrated similar anti-cancer functionality in most biological processes with subtle differences in MMP-2/9 regulation; this contradicted that of AP-2γ, whose actions potentiated cancer progression. Simultaneous overexpression of WWOX and AP-2α/AP-2γ revealed that single discrepancies appear in WWOX-AP-2α collaboration but only at the highest BLCA grade; WWOX-AP-2α collaboration was considered anti-cancer. However, WWOX only appeared to have residual activity against oncogenic AP-2γ in grade 3 and 4: variants with either AP-2γ overexpression alone or combined WWOX and AP-2γ overexpression demonstrated similar pro-tumoral behavior. Transcriptome profiling with further gene ontology certified biological processes investigated in vitro and indicated groups of genes consisting of AP-2 targets and molecules worth investigation as biomarkers. In conclusion, tumor suppressor synergism between WWOX and AP-2α is unimpaired in high-grade BLCA compared to intermediate grade, yet the ability of WWOX to guide oncogenic AP-2γ is almost completely lost.

Discovery of widespread transcription initiation at microsatellites predictable by sequence-based deep neural network

Using the Cap Analysis of Gene Expression (CAGE) technology, the FANTOM5 consortium provided one of the most comprehensive maps of transcription start sites (TSSs) in several species. Strikingly, ~72% of them could not be assigned to a specific gene and initiate at unconventional regions, outside promoters or enhancers. Here, we probe these unassigned TSSs and show that, in all species studied, a significant fraction of CAGE peaks initiate at microsatellites, also called short tandem repeats (STRs). To confirm this transcription, we develop Cap Trap RNA-seq, a technology which combines cap trapping and long read MinION sequencing. We train sequence-based deep learning models able to predict CAGE signal at STRs with high accuracy. These models unveil the importance of STR surrounding sequences not only to distinguish STR classes, but also to predict the level of transcription initiation. Importantly, genetic variants linked to human diseases are preferentially found at STRs with high transcription initiation level, supporting the biological and clinical relevance of transcription initiation at STRs. Together, our results extend the repertoire of non-coding transcription associated with DNA tandem repeats and complexify STR polymorphism.

Characterization of the human TARDBP gene promoter

The expression of TDP-43, the main component of neuronal intracellular inclusions across a broad spectrum of ALS and FTD disorders, is developmentally regulated and studies in vivo have shown that TDP-43 overexpression can be toxic, even before observation of pathological aggregates. Starting from these observations, the regulation of its expression at transcriptional level might represent a further key element for the pathogenesis of neurodegenerative diseases. Therefore, we have characterized the human TARDBP promoter, in order to study the transcriptional mechanisms of expression. Mapping of cis-acting elements by luciferase assays in different cell outlined that the activity of the promoter seems to be higher in SH-SY5Y, Neuro2A, and HeLa than in HEK293. In addition, we tested effects of two SNPs found in the promoter region of ALS patients and observed no significant effect on transcription levels in all tested cell lines. Lastly, while TDP-43 overexpression did not affect significantly the activity of its promoter (suggesting that TDP-43 does not influence its own transcription), the presence of the 5'UTR sequence and of intron-1 splicing seem to impact positively on TDP-43 expression without affecting transcript stability. In conclusion, we have identified the region spanning nucleotides 451-230 upstream from the transcription start site as the minimal region with a significant transcription activity. These results lay an important foundation for exploring the regulation of the TARDBP gene transcription by exogenous and endogenous stimuli and the implication of transcriptional mechanisms in the pathogenesis of TDP-43 proteinopathies.