Genome-wide annotation and quantitation of translation by ribosome profiling

Ribosome profiling reveals downregulation of UMP biosynthesis as the major early response to phage infection

Bacteria have evolved diverse defense mechanisms to counter bacteriophage attacks. Genetic programs activated upon infection characterize phage-host molecular interactions and ultimately determine the outcome of the infection. In this study, we applied ribosome profiling to monitor protein synthesis during the early stages of sk1 bacteriophage infection in Lactococcus cremoris. Our analysis revealed major changes in gene expression within 5 minutes of sk1 infection. Notably, we observed a specific and severe downregulation of several pyr operons which encode enzymes required for uridine monophosphate biosynthesis. Consistent with previous findings, this is likely an attempt of the host to starve the phage of nucleotides it requires for propagation. We also observed a gene expression response that we expect to benefit the phage. This included the upregulation of 40 ribosome proteins that likely increased the host's translational capacity, concurrent with a downregulation of genes that promote translational fidelity (lepA and raiA). In addition to the characterization of host-phage gene expression responses, the obtained ribosome profiling data enabled us to identify two putative recoding events as well as dozens of loci currently annotated as pseudogenes that are actively translated. Furthermore, our study elucidated alterations in the dynamics of the translation process, as indicated by time-dependent changes in the metagene profile, suggesting global shifts in translation rates upon infection. Additionally, we observed consistent modifications in the ribosome profiles of individual genes, which were apparent as early as 2 minutes post-infection. The study emphasizes our ability to capture rapid alterations of gene expression during phage infection through ribosome profiling.

IMPORTANCE

The ribosome profiling technology has provided invaluable insights for understanding cellular translation and eukaryotic viral infections. However, its potential for investigating host-phage interactions remains largely untapped. Here, we applied ribosome profiling to Lactococcus cremoris cultures infected with sk1, a major infectious agent in dairy fermentation processes. This revealed a profound downregulation of genes involved in pyrimidine nucleotide synthesis at an early stage of phage infection, suggesting an anti-phage program aimed at restricting nucleotide availability and, consequently, phage propagation. This is consistent with recent findings and contributes to our growing appreciation for the role of nucleotide limitation as an anti-viral strategy. In addition to capturing rapid alterations in gene expression levels, we identified translation occurring outside annotated regions, as well as signatures of non-standard translation mechanisms. The gene profiles revealed specific changes in ribosomal densities upon infection, reflecting alterations in the dynamics of the translation process.

It's a competitive business

A new in vitro system called Rec-Seq sheds light on how mRNA molecules compete for the machinery that translates their genetic sequence into proteins.

Balanced mitochondrial and cytosolic translatomes underlie the biogenesis of human respiratory complexes

BACKGROUND

Oxidative phosphorylation (OXPHOS) complexes consist of nuclear and mitochondrial DNA-encoded subunits. Their biogenesis requires cross-compartment gene regulation to mitigate the accumulation of disproportionate subunits. To determine how human cells coordinate mitochondrial and nuclear gene expression processes, we tailored ribosome profiling for the unique features of the human mitoribosome.

RESULTS

We resolve features of mitochondrial translation initiation and identify a small ORF in the 3' UTR of MT-ND5. Analysis of ribosome footprints in five cell types reveals that average mitochondrial synthesis levels correspond precisely to cytosolic levels across OXPHOS complexes, and these average rates reflect the relative abundances of the complexes. Balanced mitochondrial and cytosolic synthesis does not rely on rapid feedback between the two translation systems, and imbalance caused by mitochondrial translation deficiency is associated with the induction of proteotoxicity pathways.

CONCLUSIONS

Based on our findings, we propose that human OXPHOS complexes are synthesized proportionally to each other, with mitonuclear balance relying on the regulation of OXPHOS subunit translation across cellular compartments, which may represent a proteostasis vulnerability.

Using DMS-MaPseq to Uncover the Roles of DEAD-box Proteins in Ribosome Assembly

DMS (dimethylsulfate) is a time-tested chemical probe for nucleic acid secondary structure that has recently re-emerged as a powerful tool to study RNA structure and structural changes, by coupling it to high throughput sequencing techniques. This variant, termed DMS-MaPseq, allows for mapping of all RNAs in a cell at the same time. However, if an RNA adopts different structures, for example during the assembly of an RNA-protein complex, or as part of its functional cycle, then DMS-MaPseq cannot differentiate between these structures, and an ensemble average will be produced. This is especially challenging for long-lived RNAs, such as ribosomes, whose steady-state abundance far exceeds that of any assembly intermediates, rendering those inaccessible to DMS-MaPseq on total RNAs. These challenges can be overcome by purification of assembly intermediates stalled at specific assembly steps (or steps in the functional cycle), via a combination of affinity tags and mutants stalled at defined steps, and subsequent DMS probing of these intermediates. Interpretation of the differences in DMS accessibility is facilitated by additional structural information, e.g. from cryo-EM experiments, available for many functional RNAs. While this approach is generally useful for studying RNA folding or conformational changes within RNA-protein complexes, it can be particularly valuable for studying the role(s) of DEAD-box proteins, as these tend to lead to larger conformational rearrangements, often resulting from the release of an RNA-binding protein from a bound RNA. Here we provide an adaptation of the DMS-MaPseq protocol to study RNA conformational transitions during ribosome assembly, which addresses the challenges arising from the presence of many assembly intermediates, all at concentrations far below that of mature ribosomes. While this protocol was developed for the yeast S. cerevisiae, we anticipate that it should be readily transferable to other model organisms for which affinity purification has been established.

Optimization of Ribosome Footprinting Conditions for Ribo-Seq in Human and Drosophila melanogaster Tissue Culture Cells

Our understanding of mRNA translation and its regulation has been transformed by the development of ribosome profiling. This approach relies upon RNase footprinting of translating ribosomes in a precise manner to generate an accurate snapshot of ribosome positions with nucleotide resolution. Here we tested a variety of conditions, which contribute to the preciseness of ribosome footprinting and therefore the success of ribosome profiling. We found that NaCl concentration, RNaseI source, RNaseI amount, and temperature of footprinting all contributed to the quality of ribosome footprinting in human neuroblastoma SH-SY5Y cells. These ideal conditions for footprinting also improved footprint quality when used with Drosophila melanogaster S2 cells. Footprinting under the same conditions generated different footprints sizes and framing patterns in human and D. melanogaster cells. We also found that treatment of S2 cells with cycloheximide prior to footprinting impacted the distribution of footprints across ORFs, without affecting overall read length distribution and framing pattern, as previously found in other organisms. Together our results indicate that a variety of factors affect ribosome footprint quality and the nature of precise footprinting varies across species.

Translation complex stabilization on messenger RNA and footprint profiling to study the RNA responses and dynamics of protein biosynthesis in the cells

During protein biosynthesis, ribosomes bind to messenger (m)RNA, locate its protein-coding information, and translate the nucleotide triplets sequentially as codons into the corresponding sequence of amino acids, forming proteins. Non-coding mRNA features, such as 5' and 3' untranslated regions (UTRs), start sites or stop codons of different efficiency, stretches of slower or faster code and nascent polypeptide interactions can alter the translation rates transcript-wise. Most of the homeostatic and signal response pathways of the cells converge on individual mRNA control, as well as alter the global translation output. Among the multitude of approaches to study translational control, one of the most powerful is to infer the locations of translational complexes on mRNA based on the mRNA fragments protected by these complexes from endonucleolytic hydrolysis, or footprints. Translation complex profiling by high-throughput sequencing of the footprints allows to quantify the transcript-wise, as well as global, alterations of translation, and uncover the underlying control mechanisms by attributing footprint locations and sizes to different configurations of the translational complexes. The accuracy of all footprint profiling approaches critically depends on the fidelity of footprint generation and many methods have emerged to preserve certain or multiple configurations of the translational complexes, often in challenging biological material. In this review, a systematic summary of approaches to stabilize translational complexes on mRNA for footprinting is presented and major findings are discussed. Future directions of translation footprint profiling are outlined, focusing on the fidelity and accuracy of inference of the native in vivo translation complex distribution on mRNA.

Cytoplasmic long noncoding RNAs are differentially regulated and translated during human neuronal differentiation

The expression of long noncoding RNAs is highly enriched in the human nervous system. However, the function of neuronal lncRNAs in the cytoplasm and their potential translation remains poorly understood. Here we performed Poly-Ribo-Seq to understand the interaction of lncRNAs with the translation machinery and the functional consequences during neuronal differentiation of human SH-SY5Y cells. We discovered 237 cytoplasmic lncRNAs up-regulated during early neuronal differentiation, 58%-70% of which are associated with polysome translation complexes. Among these polysome-associated lncRNAs, we find 45 small ORFs to be actively translated, 17 specifically upon differentiation. Fifteen of 45 of the translated lncRNA-smORFs exhibit sequence conservation within Hominidea, suggesting they are under strong selective constraint in this clade. The profiling of publicly available data sets revealed that 8/45 of the translated lncRNAs are dynamically expressed during human brain development, and 22/45 are associated with cancers of the central nervous system. One translated lncRNA we discovered is LINC01116, which is induced upon differentiation and contains an 87 codon smORF exhibiting increased ribosome profiling signal upon differentiation. The resulting LINC01116 peptide localizes to neurites. Knockdown of LINC01116 results in a significant reduction of neurite length in differentiated cells, indicating it contributes to neuronal differentiation. Our findings indicate cytoplasmic lncRNAs interact with translation complexes, are a noncanonical source of novel peptides, and contribute to neuronal function and disease. Specifically, we demonstrate a novel functional role for LINC01116 during human neuronal differentiation.

Protein Stability Buffers the Cost of Translation Attenuation following eIF2α Phosphorylation

Phosphorylation of the translation initiation factor eIF2α is a rapid and vital response to many forms of stress, including protein-misfolding stress in the endoplasmic reticulum (ER stress). It is believed to cause a general reduction in protein synthesis while enabling translation of few transcripts. Such a reduction of protein synthesis comes with the threat of depleting essential proteins, a risk thought to be mitigated by its transient nature. Here, we find that translation attenuation is not uniform, with cytosolic and mitochondrial ribosomal subunits being prominently downregulated. Translation attenuation of these targets persists after translation recovery. Surprisingly, this occurs without a measurable decrease in ribosomal proteins. Explaining this conundrum, translation attenuation preferentially targets long-lived proteins, a finding not only demonstrated by ribosomal proteins but also observed at a global level. This shows that protein stability buffers the cost of translational attenuation, establishing an evolutionary principle of cellular robustness.

deltaTE: Detection of Translationally Regulated Genes by Integrative Analysis of Ribo-seq and RNA-seq Data

Ribosome profiling quantifies the genome‐wide ribosome occupancy of transcripts. With the integration of matched RNA sequencing data, the translation efficiency (TE) of genes can be calculated to reveal translational regulation. This layer of gene‐expression regulation is otherwise difficult to assess on a global scale and generally not well understood in the context of human disease. Current statistical methods to calculate differences in TE have low accuracy, cannot accommodate complex experimental designs or confounding factors, and do not categorize genes into buffered, intensified, or exclusively translationally regulated genes. This article outlines a method [referred to as deltaTE (ΔTE), standing for change in TE] to identify translationally regulated genes, which addresses the shortcomings of previous methods. In an extensive benchmarking analysis, ΔTE outperforms all methods tested. Furthermore, applying ΔTE on data from human primary cells allows detection of substantially more translationally regulated genes, providing a clearer understanding of translational regulation in pathogenic processes. In this article, we describe protocols for data preparation, normalization, analysis, and visualization, starting from raw sequencing files. © 2019 The Authors.

Basic Protocol: One‐step detection and classification of differential translation efficiency genes using DTEG.R

Alternate Protocol: Step‐wise detection and classification of differential translation efficiency genes using R

Support Protocol: Workflow from raw data to read counts

Active Tonic mTORC1 Signals Shape Baseline Translation in Naive T Cells

Naive CD4⁺ T cells are an example of dynamic cell homeostasis: T cells need to avoid autoreactivity while constantly seeing self-peptides, yet they must be primed to react to foreign antigens during infection. The instructive signals that balance this primed yet quiescent state are unknown. Interactions with self-peptides result in membrane-proximal, tonic signals in resting T cells. Here we reveal selective and robust tonic mTORC1 signals in CD4⁺ T cells that influence T cell fate decisions. We find that the Ras exchange factor Rasgrp1 is necessary to generate tonic mTORC1 signals. Genome-wide ribosome profiling of resting, primary CD4⁺ T cells uncovers a baseline translational landscape rich in mTOR targets linked to mitochondria, oxidative phosphorylation, and splicing. Aberrantly increased tonic mTORC1 signals from a Rasgrp1^Anaef allele result in immunopathology with spontaneous appearance of T peripheral helper cells, follicular helper T cells, and anti-nuclear antibodies that are preceded by subtle alterations in the translational landscape.

Myers et al. evaluate a mouse model of autoimmunity, Rasgrp1^Anaef. They find that T cells with the Rasgrp1^Anaef allele exhibit altered signaling from Rasgrp1 to the mTORC1 pathway in the basal state. They show that increased basal Rasgrp1^Anaef-mTORC1 signals lead to an altered translational landscape in T cells and immunopathology.

mTOR and other effector kinase signals that impact T cell function and activity

T cells play important roles in autoimmune diseases and cancer. Following the cloning of the T cell receptor (TCR), the race was on to map signaling proteins that contributed to T cell activation downstream of the TCR as well as co-stimulatory molecules such as CD28. We term this "canonical TCR signaling" here. More recently, it has been appreciated that T cells need to accommodate increased metabolic needs that stem from T cell activation in order to function properly. A central role herein has emerged for mechanistic/mammalian target of rapamycin (mTOR). In this review we briefly cover canonical TCR signaling to set the stage for discussion on mTOR signaling, mRNA translation, and metabolic adaptation in T cells. We also discuss the role of mTOR in follicular helper T cells, regulatory T cells, and other T cell subsets. Our lab recently uncovered that "tonic signals", which pass through proximal TCR signaling components, are robustly and selectively transduced to mTOR to promote baseline translation of various mRNA targets. We discuss insights on (tonic) mTOR signaling in the context of T cell function in autoimmune diseases such as lupus as well as in cancer immunotherapy through CAR-T cell or checkpoint blockade approaches.

mTORC1 signalling is not essential for the maintenance of muscle mass and function in adult sedentary mice

BACKGROUND

The balance between protein synthesis and degradation (proteostasis) is a determining factor for muscle size and function. Signalling via the mammalian target of rapamycin complex 1 (mTORC1) regulates proteostasis in skeletal muscle by affecting protein synthesis and autophagosomal protein degradation. Indeed, genetic inactivation of mTORC1 in developing and growing muscle causes atrophy resulting in a lethal myopathy. However, systemic dampening of mTORC1 signalling by its allosteric inhibitor rapamycin is beneficial at the organismal level and increases lifespan. Whether the beneficial effect of rapamycin comes at the expense of muscle mass and function is yet to be established.

METHODS

We conditionally ablated the gene coding for the mTORC1-essential component raptor in muscle fibres of adult mice [inducible raptor muscle-specific knockout (iRAmKO)]. We performed detailed phenotypic and biochemical analyses of iRAmKO mice and compared them with muscle-specific raptor knockout (RAmKO) mice, which lack raptor in developing muscle fibres. We also used polysome profiling and proteomics to assess protein translation and associated signalling in skeletal muscle of iRAmKO mice.

RESULTS

Analysis at different time points reveal that, as in RAmKO mice, the proportion of oxidative fibres decreases, but slow-type fibres increase in iRAmKO mice. Nevertheless, no significant decrease in body and muscle mass or muscle fibre area was detected up to 5 months post-raptor depletion. Similarly, ex vivo muscle force was not significantly reduced in iRAmKO mice. Despite stable muscle size and function, inducible raptor depletion significantly reduced the expression of key components of the translation machinery and overall translation rates.

CONCLUSIONS

Raptor depletion and hence complete inhibition of mTORC1 signalling in fully grown muscle leads to metabolic and morphological changes without inducing muscle atrophy even after 5 months. Together, our data indicate that maintenance of muscle size does not require mTORC1 signalling, suggesting that rapamycin treatment is unlikely to negatively affect muscle mass and function.

Accurate annotation of human protein-coding small open reading frames

Functional protein-coding small open reading frames (smORFs) are emerging as an important class of genes. However, the number of translated smORFs in the human genome is unclear because proteogenomic methods are not sensitive enough, and, as we show, Ribo-Seq strategies require additional measures to ensure comprehensive and accurate smORF annotation. Here, we integrate de novo transcriptome assembly and Ribo-Seq into an improved workflow that overcomes obstacles with previous methods to more confidently annotate thousands of smORFs. Evolutionary conservation analyses suggest that hundreds of smORF-encoded microproteins are likely functional. Additionally, many smORFs are regulated during fundamental biological processes, such as cell stress. Peptides derived from smORFs are also detectable on human leukocyte antigen complexes, revealing smORFs as a source of antigens. Thus, by including additional validation into our smORF annotation workflow, we accurately identify thousands of unannotated translated smORFs that will provide a rich pool of unexplored, functional human genes.

Polyamines Control eIF5A Hypusination, TFEB Translation, and Autophagy to Reverse B Cell Senescence

Failure to make adaptive immune responses is a hallmark of aging. Reduced B cell function leads to poor vaccination efficacy and a high prevalence of infections in the elderly. Here we show that reduced autophagy is a central molecular mechanism underlying immune senescence. Autophagy levels are specifically reduced in mature lymphocytes, leading to compromised memory B cell responses in old individuals. Spermidine, an endogenous polyamine metabolite, induces autophagy in vivo and rejuvenates memory B cell responses. Mechanistically, spermidine post-translationally modifies the translation factor eIF5A, which is essential for the synthesis of the autophagy transcription factor TFEB. Spermidine is depleted in the elderly, leading to reduced TFEB expression and autophagy. Spermidine supplementation restored this pathway and improved the responses of old human B cells. Taken together, our results reveal an unexpected autophagy regulatory mechanism mediated by eIF5A at the translational level, which can be harnessed to reverse immune senescence in humans.

Proteogenomic Annotation of Chinese Hamsters Reveals Extensive Novel Translation Events and Endogenous Retroviral Elements

A high-quality genome annotation greatly facilitates successful cell line engineering. Standard draft genome annotation pipelines are based largely on de novo gene prediction, homology, and RNA-Seq data. However, draft annotations can suffer from incorrect predictions of translated sequence, inaccurate splice isoforms, and missing genes. Here, we generated a draft annotation for the newly assembled Chinese hamster genome and used RNA-Seq, proteomics, and Ribo-Seq to experimentally annotate the genome. We identified 3529 new proteins compared to the hamster RefSeq protein annotation and 2256 novel translational events (e.g., alternative splices, mutations, and novel splices). Finally, we used this pipeline to identify the source of translated retroviruses contaminating recombinant products from Chinese hamster ovary (CHO) cell lines, including 119 type-C retroviruses, thus enabling future efforts to eliminate retroviruses to reduce the costs incurred with retroviral particle clearance. In summary, the improved annotation provides a more accurate resource for CHO cell line engineering, by facilitating the interpretation of omics data, defining of cellular pathways, and engineering of complex phenotypes.

RNA regulons are essential in intestinal homeostasis

Intestinal epithelial cells are among the most rapidly proliferating cell types in the human body. There are several different subtypes of epithelial cells, each with unique functional roles in responding to the ever-changing environment. The epithelium's ability for rapid and customized responses to environmental changes requires multitiered levels of gene regulation. An emerging paradigm in gastrointestinal epithelial cells is the regulation of functionally related mRNA families, or regulons, via RNA-binding proteins (RBPs). RBPs represent a rapid and efficient mechanism to regulate gene expression and cell function. In this review, we will provide an overview of intestinal epithelial RBPs and how they contribute specifically to intestinal epithelial stem cell dynamics. In addition, we will highlight key gaps in knowledge in the global understanding of RBPs in gastrointestinal physiology as an opportunity for future studies.

RibORF: Identifying Genome-Wide Translated Open Reading Frames Using Ribosome Profiling

Ribosome profiling identifies RNA fragments associated with translating ribosomes. The technology provides an opportunity to examine genome-wide translation events at single-nucleotide resolution and in an unbiased manner. Here I present a computational pipeline named RibORF to systematically identify translated open reading frames (ORFs), based on read distribution features representing active translation, including 3-nt periodicity and uniformness across codons. Analyses using the computational tool revealed pervasive translation in putative 'noncoding' regions, such as lncRNAs, pseudogenes, and 5'UTRs. The computational tool is useful for studying functional roles of non-canonical translation events in various biological processes. © 2018 by John Wiley & Sons, Inc.

Pimp My Ribosome: Ribosomal Protein Paralogs Specify Translational Control

The ability of cells to grow and divide, differentiate and function, and even senesce is dependent on the fine-tuning of both gene and protein expression. Protein concentration in the cell is regulated not only at the transcriptional and post-translational levels, but also at the level of translation. Ribosomes, the molecular machines behind translation, were once considered to be an invariant driving force behind protein expression. However, studies over the past decade paint a rather different picture; namely, that ribosomes constitute an additional layer of regulatory control that might define which subsets of mRNAs are translated, to what extent, and to what purpose. Recent works summarized herein directly implicate ribosome heterogeneity and, in particular, ribosomal protein (RP) paralog specificity in regulating mRNA translation and control of the cellular translatome.

Rfoot: Transcriptome-Scale Identification of RNA-Protein Complexes from Ribosome Profiling Data

Ribosome profiling was developed to identify genome-wide RNA fragments associated with translating ribosomes. However, no experimental procedure was developed to effectively purify ribosome-RNA complexes. Actually, ribosome profiling is a transcriptomic RNase footprinting assay, which can identify both ribosomal and non-ribosomal protein-RNA complexes. Many sequencing reads represent functional RNA footprints associated with non-ribosomal proteins. Here I present a computational pipeline named Rfoot to systematically identify genome-wide non-ribosomal RNA footprints, based on the highly localized read distribution feature. Analyses have revealed native functional protein-RNA complexes in lncRNAs, 3'UTRs of mRNAs, and all types of small noncoding RNAs. This computational tool is useful for discovering novel noncoding functions of RNAs. © 2018 by John Wiley & Sons, Inc.

Novel insights into global translational regulation through Pumilio family RNA-binding protein Puf3p revealed by ribosomal profiling

RNA binding proteins (RBPs) can regulate the stability, localization, and translation of their target mRNAs. Among them, Puf3p is a well-known Pumilio family RBP whose biology has been intensively studied. Nevertheless, the impact of Puf3p on the translational regulation of its downstream genes still remains to be investigated at the genome-wide level. In this study, we combined ribosome profiling and RNA-Seq in budding yeast (Saccharomyces cerevisiae) to investigate Puf3p's functions in translational regulation. Comparison of translational efficiency (TE) between wild-type and puf3Δ strains demonstrates extensive translational modulation in the absence of Puf3p (over 27% genes are affected at the genome level). Besides confirming its known role in regulating mitochondrial metabolism, our data demonstrate that Puf3p serves as a key post-transcriptional regulator of downstream RBPs by regulating their translational efficiencies, indicating a network of interactions among RBPs at the post-transcriptional level. Furthermore, Puf3p switches the balance of translational flux between mitochondrial and cytosolic ribosome biogenesis to adapt to changes in cellular metabolism. In summary, our results indicate that TE can be utilized as an informative index to interrogate the mechanism underlying RBP functions, and provide novel insights into Puf3p's mode-of-action.

The exon-intron gene structure upstream of the initiation codon predicts translation efficiency

Introns in mRNA leaders are common in complex eukaryotes, but often overlooked. These introns are spliced out before translation, leaving exon-exon junctions in the mRNA leaders (leader EEJs). Our multi-omic approach shows that the number of leader EEJs inversely correlates with the main protein translation, as does the number of upstream open reading frames (uORFs). Across the five species studied, the lowest levels of translation were observed for mRNAs with both leader EEJs and uORFs (29%). This class of mRNAs also have ribosome footprints on uORFs, with strong triplet periodicity indicating uORF translation. Furthermore, the positions of both leader EEJ and uORF are conserved between human and mouse. Thus, the uORF, in combination with leader EEJ predicts lower expression for nearly one-third of eukaryotic proteins.

The impact of ribosomal interference, codon usage, and exit tunnel interactions on translation elongation rate variation

Previous studies have shown that translation elongation is regulated by multiple factors, but the observed heterogeneity remains only partially explained. To dissect quantitatively the different determinants of elongation speed, we use probabilistic modeling to estimate initiation and local elongation rates from ribosome profiling data. This model-based approach allows us to quantify the extent of interference between ribosomes on the same transcript. We show that neither interference nor the distribution of slow codons is sufficient to explain the observed heterogeneity. Instead, we find that electrostatic interactions between the ribosomal exit tunnel and specific parts of the nascent polypeptide govern the elongation rate variation as the polypeptide makes its initial pass through the tunnel. Once the N-terminus has escaped the tunnel, the hydropathy of the nascent polypeptide within the ribosome plays a major role in modulating the speed. We show that our results are consistent with the biophysical properties of the tunnel.

The long non-coding RNA PARROT is an upstream regulator of c-Myc and affects proliferation and translation

Long non-coding RNAs are important regulators of gene expression and signaling pathways. The expression of long ncRNAs is dysregulated in cancer and other diseases. The identification and characterization of long ncRNAs is often challenging due to their low expression level and localization to chromatin. Here, we identify a functional long ncRNA, PARROT (Proliferation Associated RNA and Regulator Of Translation) transcribed by RNA polymerase II and expressed at a relatively high level in a number of cell lines. The PARROT long ncRNA is associated with proliferation in both transformed and normal cell lines. We characterize the long ncRNA PARROT as an upstream regulator of c-Myc affecting cellular proliferation and translation using RNA sequencing and mass spectrometry following depletion of the long ncRNA. PARROT is repressed during senescence of human mammary epithelial cells and overexpressed in some cancers, suggesting an important association with proliferation through regulation of c-Myc. With this study, we add to the knowledge of cytoplasmic functional long ncRNAs and extent the long ncRNA-Myc regulatory network in transformed and normal cells.

A Detailed Protocol for Subcellular RNA Sequencing (subRNA-seq)

In eukaryotic cells, RNAs at various maturation and processing levels are distributed across cellular compartments. The standard approach to determine transcript abundance and identity in vivo is RNA sequencing (RNA-seq). RNA-seq relies on RNA isolation from whole-cell lysates and thus mainly captures fully processed, stable, and more abundant cytoplasmic RNAs over nascent, unstable, and nuclear RNAs. Here, we provide a step-by-step protocol for subcellular RNA-seq (subRNA-seq). subRNA-seq allows the quantitative measurement of RNA polymerase II-generated RNAs from the chromatin, nucleoplasm, and cytoplasm of mammalian cells. This approach relies on cell fractionation prior to RNA isolation and sequencing library preparation. High-throughput sequencing of the subcellular RNAs can then be used to reveal the identity, abundance, and subcellular distribution of transcripts, thus providing insights into RNA processing and maturation. Deep sequencing of the chromatin-associated RNAs further offers the opportunity to study nascent RNAs. Subcellular RNA-seq libraries are obtained within 5 days. © 2017 by John Wiley & Sons, Inc.

Use of Antibiotics for Maintenance of Axenic Cultures of Amphidinium carterae for the Analysis of Translation

Most dinoflagellates in culture are bacterized, complicating the quantification of protein synthesis, as well as the analysis of its regulation. In bacterized cultures of Amphidinium carterae Hulbert, up to 80% of protein synthetic activity appears to be predominantly bacterial based on responses to inhibitors of protein synthesis. To circumvent this, axenic cultures of A. carterae were obtained and shown to respond to inhibitors of protein synthesis in a manner characteristic of eukaryotes. However, these responses changed with time in culture correlating with the reappearance of bacteria. Here we show that culture with kanamycin (50 μg/mL), carbenicillin (100 μg/mL), and streptomycin sulfate (50 μg/mL) (KCS), but not 100 units/mL of penicillin and streptomycin (PS), prevents the reappearance of bacteria and allows A. carterae protein synthesis to be quantified without the contribution of an associated bacterial community. We demonstrate that A. carterae can grow in the absence of a bacterial community. Furthermore, maintenance in KCS does not inhibit the growth of A. carterae cultures but slightly extends the growth phase and allows accumulation to somewhat higher saturation densities. We also show that cultures of A. carterae maintained in KCS respond to the eukaryotic protein synthesis inhibitors cycloheximide, emetine, and harringtonine. Establishment of these culture conditions will facilitate our ability to use polysome fractionation and ribosome profiling to study mRNA recruitment. Furthermore, this study shows that a simple and fast appraisal of the presence of a bacterial community in A. carterae cultures can be made by comparing responses to cycloheximide and chloramphenicol rather than depending on lengthier culture-based assessments.

Selective stalling of human translation through small-molecule engagement of the ribosome nascent chain

Proprotein convertase subtilisin/kexin type 9 (PCSK9) plays a key role in regulating the levels of plasma low-density lipoprotein cholesterol (LDL-C). Here, we demonstrate that the compound PF-06446846 inhibits translation of PCSK9 by inducing the ribosome to stall around codon 34, mediated by the sequence of the nascent chain within the exit tunnel. We further show that PF-06446846 reduces plasma PCSK9 and total cholesterol levels in rats following oral dosing. Using ribosome profiling, we demonstrate that PF-06446846 is highly selective for the inhibition of PCSK9 translation. The mechanism of action employed by PF-06446846 reveals a previously unexpected tunability of the human ribosome that allows small molecules to specifically block translation of individual transcripts.

Many disease-mediating proteins have proven difficult to target with traditional small-molecule pharmaceuticals. In this paper, we report that a small molecule, PF-06446846, directly inhibits translation of one such protein, proprotein convertase subtilisin/kexin type 9 (PCSK9), by acting on the translating human ribosome. PF-06446846 causes the translating ribosome to stall soon after translating the PCSK9 signal sequence. We further show that PF-06446846 activity is dependent on the amino acid sequence of the nascent chain inside the ribosome exit tunnel. In a rat safety study, we observe decreases in plasma PCSK9, total cholesterol, and low-density lipoprotein (LDL) cholesterol. Using mass spectrometry in cell culture and ribosome profiling, we demonstrate that despite acting on the ribosome, which synthesizes every protein in the cell, PF-06446846 displays a high level of selectivity for PCSK9. This unexpected potential for small molecules to selectively inhibit the human ribosome opens the possibility for future development of small molecules targeting disease-mediating proteins that were previously thought to be undruggable.

Kc167, a widely used Drosophila cell line, contains an active primary piRNA pathway

PIWI family proteins bind to small RNAs known as PIWI-interacting RNAs (piRNAs) and play essential roles in the germline by silencing transposons and by promoting germ cell specification and function. Here we report that the widely used Kc167 cell line, derived from Drosophila melanogaster embryos, expresses piRNAs that are loaded to Aub and Piwi. Kc167 piRNAs are produced by a canonical, primary piRNA biogenesis pathway, from phased processing of precursor transcripts by the Zuc endonuclease, Armi helicase, and dGasz mitochondrial scaffold protein. Kc167 piRNAs derive from cytoplasmic transcripts, notably tRNAs and mRNAs, and their abundance correlates with that of parent transcripts. The expression of Aub is robust in Kc167, that of Piwi is modest, while Ago3 is undetectable, explaining the lack of transposon-related piRNA amplification by the Aub-Ago3, ping-pong mechanism. We propose that the default state of the primary piRNA biogenesis machinery is random transcript sampling to allow generation of piRNAs from any transcript, including newly acquired retrotransposons. This state is unmasked in Kc167, likely because they do not express piRNA cluster transcripts in sufficient amounts and do not amplify transposon piRNAs. We use Kc167 to characterize an inactive isoform of Aub protein. Since most Kc167 piRNAs are genic, they can be mapped uniquely to the genome, facilitating computational analyses. Furthermore, because Kc167 is a widely used and well-characterized cell line that is easily amenable to experimental manipulations, we expect that it will serve as an excellent system to study piRNA biogenesis and piRNA-related factors.

Analysis of the translatome in solid tumors using polyribosome profiling/RNA-Seq

Gene expression involves multiple steps from the transcription of a mRNA in the nucleus to the production of the encoded protein in the cytoplasm. This final step occurs through a highly regulated process of mRNA translation on ribosomes that is required to maintain cell homeostasis. Alterations in the control of mRNA translation may lead to cell's transformation, a hallmark of cancer development. Indeed, recent advances indicated that increased translation of mRNAs encoding tumor-promoting proteins may be a key mechanism of tumor resistance in several cancers. Moreover, it was found that proteins whose encoding mRNAs are translated at higher efficiencies may be effective biomarkers. Evaluation of global changes in translation efficiency in human tumors has thus the potential of better understanding what can be used as biomarkers and therapeutic targets. Investigating changes in translation efficiency in human cancer cells has been made possible through the development and use of the polyribosome profiling combined with DNA microarray or deep RNA sequencing (RNA-Seq). While helpful, the use of cancer cell lines has many limitations and it is essential to define translational changes in human tumor samples in order to properly prioritize genes implicated in cancer phenotype. We present an optimized polyribosome RNA-Seq protocol suitable for quantitative analysis of mRNA translation that occurs in human tumor samples and murine xenografts. Applying this innovative approach to human tumors, which requires a complementary bioinformatics analysis, unlocks the potential to identify key mRNA which are preferentially translated in tumor tissue compared to benign tissue as well as translational changes which occur following treatment. These technical advances will be of interest to those researching all solid tumors, opening possibilities for understanding what may be therapeutic Achilles heels' or relevant biomarkers.

[An intimate look at the viral replication cycle through ribosome profiling]

Next Generation Sequencing (NGS) techniques have revolutionized most biomedical research fields over the past decade by allowing a broader vision on biological processes that occur at the molecular level. Among these, ribosome profiling or footprinting is a powerful tool to study mRNA translation in a transcriptome-wide manner. Ribosome profiling has been used to study the impact of translational control of gene expression under many different cellular conditions including viral infections. Indeed, translation is a critical step during the viral replication cycle in which the infected cell is embezzled to produce viral proteins. Ribosome profiling tools can provide new insights on viral translation by monitoring ribosome binding to viral and cellular RNAs with a high definition during the time course of an infection. Here, we describe the potential uses of ribosome profiling for the understanding of viral translational control and the impact of viral infection on host gene expression. We also discuss the main limitations and biases related to the technique that need to be taken into account for its use.

Mechanisms of protein balance in skeletal muscle

Increased global demand for adequate protein nutrition against a backdrop of climate change and concern for animal agriculture sustainability necessitates new and more efficient approaches to livestock growth and production. Anabolic growth is achieved when rates of new synthesis exceed turnover, producing a positive net protein balance. Conversely, deterioration or atrophy of lean mass is a consequence of a net negative protein balance. During early life and periods of growth, muscle mass is driven by increases in protein synthesis at the level of mRNA translation. Throughout life, muscle mass is further influenced by degradative processes such as autophagy and the ubiquitin proteasome pathway. Multiple signal transduction networks guide and coordinate these processes alongside quality control mechanisms to maintain protein homeostasis (proteostasis). Genetics, hormones, and environmental stimuli each influence proteostasis control, altering capacity and/or efficiency of muscle growth. An overview of recent findings and current methods to assess muscle protein balance and proteostasis is presented. Current efforts to identify novel control points have the potential through selective breeding design or development of hormetic strategies to better promote growth and health span during environmental stress.

Small Peptides as Newcomers in the Control of Drosophila Development

Throughout the last century, studies using the fruit fly have contributed to the discovery of many key genetic elements that control animal development. Recent work has shed light on an unexpectedly large number of RNAs that lack the classical hallmarks of protein-coding genes and are thus referred to as noncoding RNAs. However, there is mounting evidence that both mRNA and noncoding RNAs often contain small open reading frames (sORFs/smORFs), which can be translated into peptides. While genome-wide profiling supports a pervasive translation of these noncanonical sORF/smORF/SEP peptides, their functions remain poorly understood. Here, we review recent data obtained in Drosophila demonstrating the overlooked role of smORF peptides in the control of development and adult life. Focusing on a few smORF peptides whose functions have been elucidated recently, we discuss the importance of these newly identified regulatory molecules and how they act to regulate the building and function of the whole organism.

Exploring Ribosome Positioning on Translating Transcripts with Ribosome Profiling

Recent technological advances (e.g., microarrays and massively parallel sequencing) have facilitated genome-wide measurement of many aspects of gene regulation. Ribosome profiling is a high-throughput sequencing method used to measure gene expression at the level of translation. This is accomplished by quantifying both the number of translating ribosomes and their locations on mRNA transcripts. The inventors of this approach have published several methods papers detailing its implementation and addressing the basics of ribosome profiling data analysis. Here we describe our lab's procedure, which differs in some respects from those published previously. In addition, we describe a data analysis pipeline, Ribomap, for ribosome profiling data. Ribomap allocates sequence reads to alternative mRNA isoforms, normalizes sequencing bias along transcripts using RNA-seq data, and outputs count vectors of per-codon ribosome occupancy for each transcript.

A comparative genomics study on the effect of individual amino acids on ribosome stalling

BACKGROUND

During protein synthesis, the nascent peptide chain emerges from the ribosome through the ribosomal exit tunnel. Biochemical interactions between the nascent peptide and the tunnel may stall the ribosome movement and thus affect the expression level of the protein being synthesized. Earlier studies focused on one model organism (S. cerevisiae), have suggested that certain amino acid sequences may be responsible for ribosome stalling; however, the stalling effect at the individual amino acid level across many organisms has not yet been quantified.

RESULTS

By analyzing multiple ribosome profiling datasets from different organisms (including prokaryotes and eukaryotes), we report for the first time the organism-specific amino acids that significantly lead to ribosome stalling. We show that the identity of the stalling amino acids vary across the tree of life. In agreement with previous studies, we observed a remarkable stalling signal of proline and arginine in S. cerevisiae. In addition, our analysis supports the conjecture that the stalling effect of positively charged amino acids is not universal and that in certain conditions, negative charge may also induce ribosome stalling. Finally, we show that the beginning part of the tunnel tends to undergo more interactions with the translated amino acids than other positions along the tunnel.

CONCLUSIONS

The reported results support the conjecture that the ribosomal exit tunnel interacts with various amino acids and that the nature of these interactions varies among different organisms. Our findings should contribute towards better understanding of transcript and proteomic evolution and translation elongation regulation.

The use of duplex-specific nuclease in ribosome profiling and a user-friendly software package for Ribo-seq data analysis

Ribosome profiling is a technique that permits genome-wide, quantitative analysis of translation and has found broad application in recent years. Here we describe a modified profiling protocol and software package designed to benefit more broadly the translation community in terms of simplicity and utility. The protocol, applicable to diverse organisms, including organelles, is based largely on previously published profiling methodologies, but uses duplex-specific nuclease (DSN) as a convenient, species-independent way to reduce rRNA contamination. We show that DSN-based depletion compares favorably with other commonly used rRNA depletion strategies and introduces little bias. The profiling protocol typically produces high levels of triplet periodicity, facilitating the detection of coding sequences, including upstream, downstream, and overlapping open reading frames (ORFs) and an alternative ribosome conformation evident during termination of protein synthesis. In addition, we provide a software package that presents a set of methods for parsing ribosomal profiling data from multiple samples, aligning reads to coding sequences, inferring alternative ORFs, and plotting average and transcript-specific aspects of the data. Methods are also provided for extracting the data in a form suitable for differential analysis of translation and translational efficiency.

Allele-Selective Transcriptome Recruitment to Polysomes Primed for Translation: Protein-Coding and Noncoding RNAs, and RNA Isoforms

mRNA translation into proteins is highly regulated, but the role of mRNA isoforms, noncoding RNAs (ncRNAs), and genetic variants remains poorly understood. mRNA levels on polysomes have been shown to correlate well with expressed protein levels, pointing to polysomal loading as a critical factor. To study regulation and genetic factors of protein translation we measured levels and allelic ratios of mRNAs and ncRNAs (including microRNAs) in lymphoblast cell lines (LCL) and in polysomal fractions. We first used targeted assays to measure polysomal loading of mRNA alleles, confirming reported genetic effects on translation of OPRM1 and NAT1, and detecting no effect of rs1045642 (3435C>T) in ABCB1 (MDR1) on polysomal loading while supporting previous results showing increased mRNA turnover of the 3435T allele. Use of high-throughput sequencing of complete transcript profiles (RNA-Seq) in three LCLs revealed significant differences in polysomal loading of individual RNA classes and isoforms. Correlated polysomal distribution between protein-coding and non-coding RNAs suggests interactions between them. Allele-selective polysome recruitment revealed strong genetic influence for multiple RNAs, attributable either to differential expression of RNA isoforms or to differential loading onto polysomes, the latter defining a direct genetic effect on translation. Genes identified by different allelic RNA ratios between cytosol and polysomes were enriched with published expression quantitative trait loci (eQTLs) affecting RNA functions, and associations with clinical phenotypes. Polysomal RNA-Seq combined with allelic ratio analysis provides a powerful approach to study polysomal RNA recruitment and regulatory variants affecting protein translation.

Ribosome A and P sites revealed by length analysis of ribosome profiling data

The high-throughput sequencing of nuclease-protected mRNA fragments bound to ribosomes, a technique known as ribosome profiling, quantifies the relative frequencies with which different regions of transcripts are translated. This technique has revealed novel translation initiation sites with unprecedented scope and has furthered investigations into the connections between codon biases and translation rates. Yet the location of the codon being decoded in ribosome footprints is still unknown, and has been complicated by the recent observation of footprints with non-canonical lengths. Here we show how taking into account the variations in ribosome footprint lengths can reveal the ribosome aminoacyl (A) and peptidyl (P) site locations. These location assignments are in agreement with the proposed mechanisms for various ribosome pauses and further enhance the resolution of the profiling data. We also show that GC-rich motifs at the 5' ends of footprints are found in yeast, calling into question the anti-Shine-Dalgarno effect's role in ribosome pausing.

Ribosome profiling: a tool for quantitative evaluation of dynamics in mRNA translation

Translational regulation is important for plant growth, metabolism, and acclimation to environmental challenges. Ribosome profiling involves the nuclease digestion of mRNAs associated with ribosomes and mapping of the generated ribosome-protected footprints to transcripts. This is useful for investigation of translational regulation. Here we present a detailed method to generate, purify, and high-throughput-sequence ribosome footprints from Arabidopsis thaliana using two different isolation methods, namely, conventional differential centrifugation and the translating ribosome affinity purification (TRAP) technology. These methodologies provide researchers with an opportunity to quantitatively assess with high-resolution the translational activity of individual mRNAs by determination of the position and number of ribosomes in the corresponding mRNA. The results can provide insights into the translation of upstream open reading frames, alternatively spliced transcripts, short open reading frames, and other aspects of translation.

An optimized kit-free method for making strand-specific deep sequencing libraries from RNA fragments

Deep sequencing of strand-specific cDNA libraries is now a ubiquitous tool for identifying and quantifying RNAs in diverse sample types. The accuracy of conclusions drawn from these analyses depends on precise and quantitative conversion of the RNA sample into a DNA library suitable for sequencing. Here, we describe an optimized method of preparing strand-specific RNA deep sequencing libraries from small RNAs and variably sized RNA fragments obtained from ribonucleoprotein particle footprinting experiments or fragmentation of long RNAs. Our approach works across a wide range of input amounts (400 pg to 200 ng), is easy to follow and produces a library in 2–3 days at relatively low reagent cost, all while giving the user complete control over every step. Because all enzymatic reactions were optimized and driven to apparent completion, sequence diversity and species abundance in the input sample are well preserved.

The persistent contributions of RNA to eukaryotic gen(om)e architecture and cellular function

Currently, the best scenario for earliest forms of life is based on RNA molecules as they have the proven ability to catalyze enzymatic reactions and harbor genetic information. Evolutionary principles valid today become apparent in such models already. Furthermore, many features of eukaryotic genome architecture might have their origins in an RNA or RNA/protein (RNP) world, including the onset of a further transition, when DNA replaced RNA as the genetic bookkeeper of the cell. Chromosome maintenance, splicing, and regulatory function via RNA may be deeply rooted in the RNA/RNP worlds. Mostly in eukaryotes, conversion from RNA to DNA is still ongoing, which greatly impacts the plasticity of extant genomes. Raw material for novel genes encoding protein or RNA, or parts of genes including regulatory elements that selection can act on, continues to enter the evolutionary lottery.

Differential RNA-seq: the approach behind and the biological insight gained

RNA-sequencing has revolutionized the quantitative and qualitative analysis of transcriptomes in both prokaryotes and eukaryotes. It provides a generic approach for gene expression profiling, annotation of transcript boundaries and operons, as well as identifying novel transcripts including small noncoding RNA molecules and antisense RNAs. We recently developed a differential RNA-seq (dRNA-seq) method which in addition to the above, yields information as to whether a given RNA is a primary or processed transcript. Originally applied to describe the primary transcriptome of the gastric pathogen Helicobacter pylori, dRNA-seq has since provided global maps of transcriptional start sites in diverse species, informed new biology in the CRISPR-Cas9 system, advanced to a tool for comparative transcriptomics, and inspired simultaneous RNA-seq of pathogen and host.

A bicistronic MAVS transcript highlights a class of truncated variants in antiviral immunity

Bacterial and viral mRNAs are often polycistronic. Akin to alternative splicing, alternative translation of polycistronic messages is a mechanism to generate protein diversity and regulate gene function. Although a few examples exist, the use of polycistronic messages in mammalian cells is not widely appreciated. Here we report an example of alternative translation as a means of regulating innate immune signaling. MAVS, a regulator of antiviral innate immunity, is expressed from a bicistronic mRNA encoding a second protein, miniMAVS. This truncated variant interferes with interferon production induced by full-length MAVS, whereas both proteins positively regulate cell death. To identify other polycistronic messages, we carried out genome-wide ribosomal profiling and identified a class of antiviral truncated variants. This study therefore reveals the existence of a functionally important bicistronic antiviral mRNA and suggests a widespread role for polycistronic mRNAs in the innate immune system.

Library construction for next-generation sequencing: overviews and challenges

High-throughput sequencing, also known as next-generation sequencing (NGS), has revolutionized genomic research. In recent years, NGS technology has steadily improved, with costs dropping and the number and range of sequencing applications increasing exponentially. Here, we examine the critical role of sequencing library quality and consider important challenges when preparing NGS libraries from DNA and RNA sources. Factors such as the quantity and physical characteristics of the RNA or DNA source material as well as the desired application (i.e., genome sequencing, targeted sequencing, RNA-seq, ChIP-seq, RIP-seq, and methylation) are addressed in the context of preparing high quality sequencing libraries. In addition, the current methods for preparing NGS libraries from single cells are also discussed.

Comparing protein folding in vitro and in vivo: foldability meets the fitness challenge

In this review, we compare and contrast current knowledge about in vitro and in vivo protein folding. Major advances in understanding fundamental principles underlying protein folding in optimized in vitro conditions have yielded detailed physicochemical principles of folding landscapes for small, single domain proteins. In addition, there has been increased research focusing on the key features of protein folding in the cell that differentiate it from in vitro folding, such as co-translational folding, chaperone-facilitated folding, and folding in crowded conditions with many weak interactions. Yet these two research areas have not been bridged effectively in research carried out to date. This review points to gaps between the two that are ripe for future research. Moreover, we emphasize the biological selection pressures that impact protein folding in vivo and how fitness drives the evolution of protein sequences in ways that may place foldability in tension with other requirements on a given protein. We suggest that viewing the physicochemical process of protein folding through the lens of evolution will unveil new insights and pose novel challenges about in-cell folding landscapes.

Extensive localization of long noncoding RNAs to the cytosol and mono- and polyribosomal complexes

Background

Long noncoding RNAs (lncRNAs) form an abundant class of transcripts, but the function of the majority of them remains elusive. While it has been shown that some lncRNAs are bound by ribosomes, it has also been convincingly demonstrated that these transcripts do not code for proteins. To obtain a comprehensive understanding of the extent to which lncRNAs bind ribosomes, we performed systematic RNA sequencing on ribosome-associated RNA pools obtained through ribosomal fractionation and compared the RNA content with nuclear and (non-ribosome bound) cytosolic RNA pools.

Results

The RNA composition of the subcellular fractions differs significantly from each other, but lncRNAs are found in all locations. A subset of specific lncRNAs is enriched in the nucleus but surprisingly the majority is enriched in the cytosol and in ribosomal fractions. The ribosomal enriched lncRNAs include H19 and TUG1.

Conclusions

Most studies on lncRNAs have focused on the regulatory function of these transcripts in the nucleus. We demonstrate that only a minority of all lncRNAs are nuclear enriched. Our findings suggest that many lncRNAs may have a function in cytoplasmic processes, and in particular in ribosome complexes.

The Vertebrate Genome Annotation browser 10 years on

The Vertebrate Genome Annotation (VEGA) database (http://vega.sanger.ac.uk), initially designed as a community resource for browsing manual annotation of the human genome project, now contains five reference genomes (human, mouse, zebrafish, pig and rat). Its introduction pages have been redesigned to enable the user to easily navigate between whole genomes and smaller multi-species haplotypic regions of interest such as the major histocompatibility complex. The VEGA browser is unique in that annotation is updated via the Human And Vertebrate Analysis aNd Annotation (HAVANA) update track every 2 weeks, allowing single gene updates to be made publicly available to the research community quickly. The user can now access different haplotypic subregions more easily, such as those from the non-obese diabetic mouse, and display them in a more intuitive way using the comparative tools. We also highlight how the user can browse manually annotated updated patches from the Genome Reference Consortium (GRC).

RefSeq: an update on mammalian reference sequences

The National Center for Biotechnology Information (NCBI) Reference Sequence (RefSeq) database is a collection of annotated genomic, transcript and protein sequence records derived from data in public sequence archives and from computation, curation and collaboration (http://www.ncbi.nlm.nih.gov/refseq/). We report here on growth of the mammalian and human subsets, changes to NCBI’s eukaryotic annotation pipeline and modifications affecting transcript and protein records. Recent changes to NCBI’s eukaryotic genome annotation pipeline provide higher throughput, and the addition of RNAseq data to the pipeline results in a significant expansion of the number of transcripts and novel exons annotated on mammalian RefSeq genomes. Recent annotation changes include reporting supporting evidence for transcript records, modification of exon feature annotation and the addition of a structured report of gene and sequence attributes of biological interest. We also describe a revised protein annotation policy for alternatively spliced transcripts with more divergent predicted proteins and we summarize the current status of the RefSeqGene project.

Current status and new features of the Consensus Coding Sequence database

The Consensus Coding Sequence (CCDS) project (http://www.ncbi.nlm.nih.gov/CCDS/) is a collaborative effort to maintain a dataset of protein-coding regions that are identically annotated on the human and mouse reference genome assemblies by the National Center for Biotechnology Information (NCBI) and Ensembl genome annotation pipelines. Identical annotations that pass quality assurance tests are tracked with a stable identifier (CCDS ID). Members of the collaboration, who are from NCBI, the Wellcome Trust Sanger Institute and the University of California Santa Cruz, provide coordinated and continuous review of the dataset to ensure high-quality CCDS representations. We describe here the current status and recent growth in the CCDS dataset, as well as recent changes to the CCDS web and FTP sites. These changes include more explicit reporting about the NCBI and Ensembl annotation releases being compared, new search and display options, the addition of biologically descriptive information and our approach to representing genes for which support evidence is incomplete. We also present a summary of recent and future curation targets.