Comprehensive map of ribosomal 2'-O-methylation and C/D box snoRNAs in Drosophila melanogaster.

During their maturation, ribosomal RNAs (rRNAs) are decorated by hundreds of chemical modifications that participate in proper folding of rRNA secondary structures and therefore in ribosomal function. Along with pseudouridine, methylation of the 2'-hydroxyl ribose moiety (Nm) is the most abundant modification of rRNAs. The majority of Nm modifications in eukaryotes are placed by Fibrillarin, a conserved methyltransferase belonging to a ribonucleoprotein complex guided by C/D box small nucleolar RNAs (C/D box snoRNAs). These modifications impact interactions between rRNAs, tRNAs and mRNAs, and some are known to fine tune translation rates and efficiency. In this study, we built the first comprehensive map of Nm sites in Drosophila melanogaster rRNAs using two complementary approaches (RiboMethSeq and Nanopore direct RNA sequencing) and identified their corresponding C/D box snoRNAs by whole-transcriptome sequencing. We de novo identified 61 Nm sites, from which 55 are supported by both sequencing methods, we validated the expression of 106 C/D box snoRNAs and we predicted new or alternative rRNA Nm targets for 31 of them. Comparison of methylation level upon different stresses show only slight but specific variations, indicating that this modification is relatively stable in D. melanogaster. This study paves the way to investigate the impact of snoRNA-mediated 2'-O-methylation on translation and proteostasis in a whole organism.

Translation is required for miRNA-dependent decay of endogenous transcripts.

Post-transcriptional repression of gene expression by miRNAs occurs through transcript destabilization or translation inhibition. mRNA decay is known to account for most miRNA-dependent repression. However, because transcript decay occurs co-translationally, whether target translation is a requirement for miRNA-dependent transcript destabilization remains unknown. To decouple these two molecular processes, we used cytosolic long noncoding RNAs (lncRNAs) as models for endogenous transcripts that are not translated. We show that, despite interacting with the miRNA-loaded RNA-induced silencing complex, the steady-state abundance and decay rates of these transcripts are minimally affected by miRNA loss. To further validate the apparent requirement of translation for miRNA-dependent decay, we fused two lncRNA candidates to the 3'-end of a protein-coding gene reporter and found this results in their miRNA-dependent destabilization. Further analysis revealed that the few natural lncRNAs whose levels are regulated by miRNAs in mESCs tend to associate with translating ribosomes, and possibly represent misannotated micropeptides, further substantiating the necessity of target translation for miRNA-dependent transcript decay. In summary, our analyses suggest that translation is required for miRNA-dependent transcript destabilization, and demonstrate that the levels of coding and noncoding transcripts are differently affected by miRNAs.

Analysis of Eukaryotic lincRNA Sequences Indicates Signatures of Hindered Translation Linked to Selection Pressure.

Long intergenic noncoding RNAs (lincRNAs) represent a large fraction of transcribed loci in eukaryotic genomes. Although classified as noncoding, most lincRNAs contain open reading frames (ORFs), and it remains unclear why cytoplasmic lincRNAs are not or very inefficiently translated. Here, we analyzed signatures of hindered translation in lincRNA sequences from five eukaryotes, covering a range of natural selection pressures. In fission yeast and Caenorhabditis elegans, that is, species under strong selection, we detected significantly shorter ORFs, a suboptimal sequence context around start codons for translation initiation, and trinucleotides ("codons") corresponding to less abundant tRNAs than for neutrally evolving control sequences, likely impeding translation elongation. For human, we detected signatures for cell-type-specific hindrance of lincRNA translation, in particular codons in abundant cytoplasmic lincRNAs corresponding to lower expressed tRNAs than control codons, in three out of five human cell lines. We verified that varying tRNA expression levels between cell lines are reflected in the amount of ribosomes bound to cytoplasmic lincRNAs in each cell line. We further propose that codons at ORF starts are particularly important for reducing ribosome-binding to cytoplasmic lincRNA ORFs. Altogether, our analyses indicate that in species under stronger selection lincRNAs evolved sequence features generally hindering translation and support cell-type-specific hindrance of translation efficiency in human lincRNAs. The sequence signatures we have identified may improve predicting peptide-coding and genuine noncoding lincRNAs in a cell type.

Transcriptome-wide sites of collided ribosomes reveal principles of translational pausing.

Translation initiation is the major regulatory step defining the rate of protein production from an mRNA. Meanwhile, the impact of nonuniform ribosomal elongation rates is largely unknown. Using a modified ribosome profiling protocol based on footprints from two closely packed ribosomes (disomes), we have mapped ribosomal collisions transcriptome-wide in mouse liver. We uncover that the stacking of an elongating onto a paused ribosome occurs frequently and scales with translation rate, trapping ∼10% of translating ribosomes in the disome state. A distinct class of pause sites is indicative of deterministic pausing signals. Pause site association with specific amino acids, peptide motifs, and nascent polypeptide structure is suggestive of programmed pausing as a widespread mechanism associated with protein folding. Evolutionary conservation at disome sites indicates functional relevance of translational pausing. Collectively, our disome profiling approach allows unique insights into gene regulation occurring at the step of translation elongation.

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.

EPD in 2020: enhanced data visualization and extension to ncRNA promoters.

The Eukaryotic Promoter Database (EPD), available online at https://epd.epfl.ch, provides accurate transcription start site (TSS) information for promoters of 15 model organisms plus corresponding functional genomics data that can be viewed in a genome browser, queried or analyzed via web interfaces, or exported in standard formats (FASTA, BED, CSV) for subsequent analysis with other tools. Recent work has focused on the improvement of the EPD promoter viewers, which use the UCSC Genome Browser as visualization platform. Thousands of high-resolution tracks for CAGE, ChIP-seq and similar data have been generated and organized into public track hubs. Customized, reproducible promoter views, combining EPD-supplied tracks with native UCSC Genome Browser tracks, can be accessed from the organism summary pages or from individual promoter entries. Moreover, thanks to recent improvements and stabilization of ncRNA gene catalogs, we were able to release promoter collections for certain classes of ncRNAs from human and mouse. Furthermore, we developed automatic computational protocols to assign orphan TSS peaks to downstream genes based on paired-end (RAMPAGE) TSS mapping data, which enabled us to add nearly 9000 new entries to the human promoter collection. Since our last article in this journal, EPD was extended to five more model organisms: rhesus monkey, rat, dog, chicken and Plasmodium falciparum.

Readthrough of stop codons under limiting ABCE1 concentration involves frameshifting and inhibits nonsense-mediated mRNA decay.

To gain insight into the mechanistic link between translation termination and nonsense-mediated mRNA decay (NMD), we depleted the ribosome recycling factor ABCE1 in human cells, resulting in an upregulation of NMD-sensitive mRNAs. Suppression of NMD on these mRNAs occurs prior to their SMG6-mediated endonucleolytic cleavage. ABCE1 depletion caused ribosome stalling at termination codons (TCs) and increased ribosome occupancy in 3' UTRs, implying enhanced TC readthrough. ABCE1 knockdown indeed increased the rate of readthrough and continuation of translation in different reading frames, providing a possible explanation for the observed NMD inhibition, since enhanced readthrough displaces NMD activating proteins from the 3' UTR. Our results indicate that stalling at TCs triggers ribosome collisions and activates ribosome quality control. Collectively, we show that improper translation termination can lead to readthrough of the TC, presumably due to ribosome collisions pushing the stalled ribosomes into the 3' UTR, where it can resume translation in-frame as well as out-of-frame.

MGA repository: a curated data resource for ChIP-seq and other genome annotated data.

The Mass Genome Annotation (MGA) repository is a resource designed to store published next generation sequencing data and other genome annotation data (such as gene start sites, SNPs, etc.) in a completely standardised format. Each sample has undergone local processing in order the meet the strict MGA format requirements. The original data source, the reformatting procedure and the biological characteristics of the samples are described in an accompanying documentation file manually edited by data curators. 10 model organisms are currently represented: Homo sapiens, Mus musculus, Danio rerio, Drosophila melanogaster, Apis mellifera, Caenorhabditis elegans, Arabidopsis thaliana, Zea mays, Saccharomyces cerevisiae and Schizosaccharomyces pombe. As of today, the resource contains over 24 000 samples. In conjunction with other tools developed by our group (the ChIP-Seq and SSA servers), it allows users to carry out a great variety of analysis task with MGA samples, such as making aggregation plots and heat maps for selected genomic regions, finding peak regions, generating custom tracks for visualizing genomic features in a UCSC genome browser window, or downloading chromatin data in a table format suitable for local processing with more advanced statistical analysis software such as R. Home page: http://ccg.vital-it.ch/mga/.

The eukaryotic promoter database in its 30th year: focus on non-vertebrate organisms.

We present an update of the Eukaryotic Promoter Database EPD (http://epd.vital-it.ch), more specifically on the EPDnew division, which contains comprehensive organisms-specific transcription start site (TSS) collections automatically derived from next generation sequencing (NGS) data. Thanks to the abundant release of new high-throughput transcript mapping data (CAGE, TSS-seq, GRO-cap) the database could be extended to plant and fungal species. We further report on the expansion of the mass genome annotation (MGA) repository containing promoter-relevant chromatin profiling data and on improvements for the EPD entry viewers. Finally, we present a new data access tool, ChIP-Extract, which enables computational biologists to extract diverse types of promoter-associated data in numerical table formats that are readily imported into statistical analysis platforms such as R.

Multilayered Organization of Jasmonate Signalling in the Regulation of Root Growth.

Physical damage can strongly affect plant growth, reducing the biomass of developing organs situated at a distance from wounds. These effects, previously studied in leaves, require the activation of jasmonate (JA) signalling. Using a novel assay involving repetitive cotyledon wounding in Arabidopsis seedlings, we uncovered a function of JA in suppressing cell division and elongation in roots. Regulatory JA signalling components were then manipulated to delineate their relative impacts on root growth. The new transcription factor mutant myc2-322B was isolated. In vitro transcription assays and whole-plant approaches revealed that myc2-322B is a dosage-dependent gain-of-function mutant that can amplify JA growth responses. Moreover, myc2-322B displayed extreme hypersensitivity to JA that totally suppressed root elongation. The mutation weakly reduced root growth in undamaged plants but, when the upstream negative regulator NINJA was genetically removed, myc2-322B powerfully repressed root growth through its effects on cell division and cell elongation. Furthermore, in a JA-deficient mutant background, ninja1 myc2-322B still repressed root elongation, indicating that it is possible to generate JA-responses in the absence of JA. We show that NINJA forms a broadly expressed regulatory layer that is required to inhibit JA signalling in the apex of roots grown under basal conditions. By contrast, MYC2, MYC3 and MYC4 displayed cell layer-specific localisations and MYC3 and MYC4 were expressed in mutually exclusive regions. In nature, growing roots are likely subjected to constant mechanical stress during soil penetration that could lead to JA production and subsequent detrimental effects on growth. Our data reveal how distinct negative regulatory layers, including both NINJA-dependent and -independent mechanisms, restrain JA responses to allow normal root growth. Mechanistic insights from this work underline the importance of mapping JA signalling components to specific cell types in order to understand and potentially engineer the growth reduction that follows physical damage.

Influence of Rotational Nucleosome Positioning on Transcription Start Site Selection in Animal Promoters.

The recruitment of RNA-Pol-II to the transcription start site (TSS) is an important step in gene regulation in all organisms. Core promoter elements (CPE) are conserved sequence motifs that guide Pol-II to the TSS by interacting with specific transcription factors (TFs). However, only a minority of animal promoters contains CPEs. It is still unknown how Pol-II selects the TSS in their absence. Here we present a comparative analysis of promoters' sequence composition and chromatin architecture in five eukaryotic model organisms, which shows the presence of common and unique DNA-encoded features used to organize chromatin. Analysis of Pol-II initiation patterns uncovers that, in the absence of certain CPEs, there is a strong correlation between the spread of initiation and the intensity of the 10 bp periodic signal in the nearest downstream nucleosome. Moreover, promoters' primary and secondary initiation sites show a characteristic 10 bp periodicity in the absence of CPEs. We also show that DNA natural variants in the region immediately downstream the TSS are able to affect both the nucleosome-DNA affinity and Pol-II initiation pattern. These findings support the notion that, in addition to CPEs mediated selection, sequence-induced nucleosome positioning could be a common and conserved mechanism of TSS selection in animals.

The Eukaryotic Promoter Database: expansion of EPDnew and new promoter analysis tools.

We present an update of EPDNew (http://epd.vital-it.ch), a recently introduced new part of the Eukaryotic Promoter Database (EPD) which has been described in more detail in a previous NAR Database Issue. EPD is an old database of experimentally characterized eukaryotic POL II promoters, which are conceptually defined as transcription initiation sites or regions. EPDnew is a collection of automatically compiled, organism-specific promoter lists complementing the old corpus of manually compiled promoter entries of EPD. This new part is exclusively derived from next generation sequencing data from high-throughput promoter mapping experiments. We report on the recent growth of EPDnew, its extension to additional model organisms and its improved integration with other bioinformatics resources developed by our group, in particular the Signal Search Analysis and ChIP-Seq web servers.

The ChIP-Seq tools and web server: a resource for analyzing ChIP-seq and other types of genomic data.

BACKGROUND: ChIP-seq and related high-throughput chromatin profilig assays generate ever increasing volumes of highly valuable biological data. To make sense out of it, biologists need versatile, efficient and user-friendly tools for access, visualization and itegrative analysis of such data. RESULTS: Here we present the ChIP-Seq command line tools and web server, implementing basic algorithms for ChIP-seq data analysis starting with a read alignment file. The tools are optimized for memory-efficiency and speed thus allowing for processing of large data volumes on inexpensive hardware. The web interface provides access to a large database of public data. The ChIP-Seq tools have a modular and interoperable design in that the output from one application can serve as input to another one. Complex and innovative tasks can thus be achieved by running several tools in a cascade. CONCLUSIONS: The various ChIP-Seq command line tools and web services either complement or compare favorably to related bioinformatics resources in terms of computational efficiency, ease of access to public data and interoperability with other web-based tools. The ChIP-Seq server is accessible at http://ccg.vital-it.ch/chipseq/ .

EPD and EPDnew, high-quality promoter resources in the next-generation sequencing era.

The Eukaryotic Promoter Database (EPD), available online at http://epd.vital-it.ch, is a collection of experimentally defined eukaryotic POL II promoters which has been maintained for more than 25 years. A promoter is represented by a single position in the genome, typically the major transcription start site (TSS). EPD primarily serves biologists interested in analysing the motif content, chromatin structure or DNA methylation status of co-regulated promoter subsets. Initially, promoter evidence came from TSS mapping experiments targeted at single genes and published in journal articles. Today, the TSS positions provided by EPD are inferred from next-generation sequencing data distributed in electronic form. Traditionally, EPD has been a high-quality database with low coverage. The focus of recent efforts has been to reach complete gene coverage for important model organisms. To this end, we introduced a new section called EPDnew, which is automatically assembled from multiple, carefully selected input datasets. As another novelty, we started to use chromatin signatures in addition to mRNA 5'tags to locate promoters of weekly expressed genes. Regarding user interfaces, we introduced a new promoter viewer which enables users to explore promoter-defining experimental evidence in a UCSC genome browser window.

RNA 3' processing functions of Arabidopsis FCA and FPA limit intergenic transcription.

The RNA-binding proteins FCA and FPA were identified based on their repression of the flowering time regulator FLC but have since been shown to have widespread roles in the Arabidopsis thaliana genome. Here, we use whole-genome tiling arrays to show that a wide spectrum of genes and transposable elements are misexpressed in the fca-9 fpa-7 (fcafpa) double mutant at two stages of seedling development. There was a significant bias for misregulated genomic segments mapping to the 3' region of genes. In addition, the double mutant misexpressed a large number of previously unannotated genomic segments corresponding to intergenic regions. We characterized a subset of these misexpressed unannotated segments and established that they resulted from extensive transcriptional read-through, use of downstream polyadenylation sites, and alternative splicing. In some cases, the transcriptional read-through significantly reduced expression of the associated genes. FCA/FPA-dependent changes in DNA methylation were found at several loci, supporting previous associations of FCA/FPA function with chromatin modifications. Our data suggest that FCA and FPA play important roles in the A. thaliana genome in RNA 3' processing and transcription termination, thus limiting intergenic transcription.

E2F transcription factor-1 modulates expression of glutamine metabolic genes in mouse embryonic fibroblasts and uterine sarcoma cells.

Metabolic reprogramming is considered as a hallmark of cancer and is clinically exploited as a novel target for therapy. The E2F transcription factor-1 (E2F1) regulates various cellular processes, including proliferative and metabolic pathways, and acts, depending on the cellular and molecular context, as an oncogene or tumor suppressor. The latter is evident by the observation that E2f1-knockout mice develop spontaneous tumors, including uterine sarcomas. This dual role warrants a detailed investigation of how E2F1 loss impacts metabolic pathways related to cancer progression. Our data indicate that E2F1 binds to the promoter of several glutamine metabolism-related genes. Interestingly, the expression of genes in the glutamine metabolic pathway were increased in mouse embryonic fibroblasts (MEFs) lacking E2F1. In addition, we confirm that E2f1-/- MEFs are more efficient in metabolizing glutamine and producing glutamine-derived precursors for proliferation. Mechanistically, we observe a co-occupancy of E2F1 and MYC on glutamine metabolic promoters, increased MYC binding after E2F1 depletion and that silencing of MYC decreased the expression of glutamine-related genes in E2f1-/- MEFs. Analyses of transcriptomic profiles in 29 different human cancers identified uterine sarcoma that showed a negative correlation between E2F1 and glutamine metabolic genes. CRISPR/Cas9 knockout of E2F1 in the uterine sarcoma cell line SK-UT-1 confirmed elevated glutamine metabolic gene expression, increased proliferation and increased MYC binding to glutamine-related promoters upon E2F1 loss. Together, our data suggest a crucial role of E2F1 in energy metabolism and metabolic adaptation in uterine sarcoma cells.

A conditional Smg6 mutant mouse model reveals circadian clock regulation through the nonsense-mediated mRNA decay pathway.

Nonsense-mediated messenger RNA (mRNA) decay (NMD) has been intensively studied as a surveillance pathway that degrades erroneous transcripts arising from mutations or RNA processing errors. While additional roles in physiological control of mRNA stability have emerged, possible functions in mammalian physiology in vivo remain unclear. Here, we created a conditional mouse allele that allows converting the NMD effector nuclease SMG6 from wild-type to nuclease domain-mutant protein. We find that NMD down-regulation affects the function of the circadian clock, a system known to require rapid mRNA turnover. Specifically, we uncover strong lengthening of free-running circadian periods for liver and fibroblast clocks and direct NMD regulation of Cry2 mRNA, encoding a key transcriptional repressor within the rhythm-generating feedback loop. Transcriptome-wide changes in daily mRNA accumulation patterns in the entrained liver, as well as an altered response to food entrainment, expand the known scope of NMD regulation in mammalian gene expression and physiology.

Molecular basis of translation termination at noncanonical stop codons in human mitochondria.

The genetic code that specifies the identity of amino acids incorporated into proteins during protein synthesis is almost universally conserved. Mitochondrial genomes feature deviations from the standard genetic code, including the reassignment of two arginine codons to stop codons. The protein required for translation termination at these noncanonical stop codons to release the newly synthesized polypeptides is not currently known. In this study, we used gene editing and ribosomal profiling in combination with cryo-electron microscopy to establish that mitochondrial release factor 1 (mtRF1) detects noncanonical stop codons in human mitochondria by a previously unknown mechanism of codon recognition. We discovered that binding of mtRF1 to the decoding center of the ribosome stabilizes a highly unusual conformation in the messenger RNA in which the ribosomal RNA participates in specific recognition of the noncanonical stop codons.

ROS-induced ribosome impairment underlies ZAKα-mediated metabolic decline in obesity and aging.

The ribotoxic stress response (RSR) is a signaling pathway in which the p38- and c-Jun N-terminal kinase (JNK)-activating mitogen-activated protein kinase kinase kinase (MAP3K) ZAKα senses stalling and/or collision of ribosomes. Here, we show that reactive oxygen species (ROS)-generating agents trigger ribosomal impairment and ZAKα activation. Conversely, zebrafish larvae deficient for ZAKα are protected from ROS-induced pathology. Livers of mice fed a ROS-generating diet exhibit ZAKα-activating changes in ribosomal elongation dynamics. Highlighting a role for the RSR in metabolic regulation, ZAK-knockout mice are protected from developing high-fat high-sugar (HFHS) diet-induced blood glucose intolerance and liver steatosis. Finally, ZAK ablation slows animals from developing the hallmarks of metabolic aging. Our work highlights ROS-induced ribosomal impairment as a physiological activation signal for ZAKα that underlies metabolic adaptation in obesity and aging.

Genetic and physiological analysis of Rht8 in bread wheat: an alternative source of semi-dwarfism with a reduced sensitivity to brassinosteroids.

Over the next decade, wheat grain production must increase to meet the demand of a fast growing human population. One strategy to meet this challenge is to raise wheat productivity by optimizing plant stature. The Reduced height 8 (Rht8) semi-dwarfing gene is one of the few, together with the Green Revolution genes, to reduce stature of wheat (Triticum aestivum L.), and improve lodging resistance, without compromising grain yield. Rht8 is widely used in dry environments such as Mediterranean countries where it increases plant adaptability. With recent climate change, its use could become increasingly important even in more northern latitudes. In the present study, the characterization of Rht8 was furthered. Morphological analyses show that the semi-dwarf phenotype of Rht8 lines is due to shorter internodal segments along the wheat culm, achieved through reduced cell elongation. Physiological experiments show that the reduced cell elongation is not due to defective gibberellin biosynthesis or signalling, but possibly to a reduced sensitivity to brassinosteroids. Using a fine-resolution mapping approach and screening 3104 F(2) individuals of a newly developed mapping population, the Rht8 genetic interval was reduced from 20.5 cM to 1.29 cM. Comparative genomics with model genomes confined the Rht8 syntenic intervals to 3.3 Mb of the short arm of rice chromosome 4, and to 2 Mb of Brachypodium distachyon chromosome 5. The very high resolution potential of the plant material generated is crucial for the eventual cloning of Rht8.