Search

Background: Molecular subtyping of bladder cancer has revealed luminal tumors generally have a more favourable prognosis. However, some aggressive forms of variant histology, including micropapillary, are often classified luminal. In previous work, we found long non-coding RNA (lncRNA) expression profiles could identify a subgroup of luminal bladder tumors with less aggressive biology and better outcomes. Objective: In the present study, we aimed to investigate whether lncRNA expression profiles could identify high-grade T1 micropapillary bladder cancer with differential outcome. Design, setting, and participants: LncRNAs were quantified from RNA-seq data from a HGT1 bladder cancer cohort that was enriched for primary micropapillary cases (15/84). Unsupervised consensus clustering of variant lncRNAs identified a three-cluster solution, which was further characterised using a panel of micropapillary-associated biomarkers, molecular subtypes, gene signatures, and survival analysis. A single-sample genomic signature was trained using lasso-penalized logistic regression to classify micropapillary-like gene-expression, as characterised by lncRNA clustering. The genomic classifier (GC) was tested on luminal tumors derived from the TCGA cohort (N = 202). Outcome measurements and statistical analysis: Patient and tumor characteristics were compared between subgroups by using X2 tests and two-sided Wilcoxon rank-sum tests. Primary endpoints were overall, progression-free and high-grade recurrence-free survival, calculated as the date of high-grade T1 disease at TURBT till date of death from any cause, progression, or recurrence, respectively. Survival rates were estimated using weighted Kaplan–Meier (KM) curves. Results and limitations: Primary micropapillary HGT1 showed decreased FGFR3, SHH, and p53 pathway activity relative to tumors with conventional urothelial carcinoma. Many bladder cancer-associated lncRNAs were downregulated in micropapillary tumors, including UCA1, LINC00152, and MALAT1. Unsupervised consensus clustering resulted in a lncRNA cluster 1 (LC1) with worse prognosis that was enriched for primary micropapillary histology and the Luminal Unstable (LumU) molecular subtype. Interestingly, LC1 appeared to better identify aggressive HGT1 disease, compared to stratifying outcomes using primary histologic characteristics. A signature trained to identify LC1 cases showed good performance in the testing cohort, identifying seven cases with significantly worse survival (p < 0.001). Limitations include the retrospective nature of the study and the lack of a validation cohort. Conclusions: Using the lncRNA transcriptome we identified a subgroup of aggressive HGT1 bladder cancer that was enriched with micropapillary histology. These data suggest that lncRNAs can facilitate the identification of aggressive micropapillary-like tumors, potentially improving patient management.

The formation of new genes during evolution is an important motor of functional innovation, but the rate at which new genes originate and the likelihood that they persist over longer evolutionary periods are still poorly understood questions. Two important mechanisms by which new genes arise are gene duplication and de novo formation from a previously noncoding sequence. Does the mechanism of formation influence the evolutionary trajectories of the genes? Proteins arisen by gene duplication retain the sequence and structural properties of the parental protein, and thus they may be relatively stable. Instead, de novo originated proteins are often species specific and thought to be more evolutionary labile. Despite these differences, here we show that both types of genes share a number of similarities, including low sequence constraints in their initial evolutionary phases, high turnover rates at the species level, and comparable persistence rates in deeper branchers, in both yeast and flies. In addition, we show that putative de novo proteins have an excess of substitutions between charged amino acids compared with the neutral expectation, which is reflected in the rapid loss of their initial highly basic character. The study supports high evolutionary dynamics of different kinds of new genes at the species level, in sharp contrast with the stability observed at later stages. We acknowledge funding from Ministerio de Ciencia e Innovación Agencia Estatal de Investigación grant PGC2018-094091-B-I00 (cofunded by Fondo Europeo de Desarrollo Regional), as well as grants PID2021-122726NB-I00 and PID2021-122830OB-C43 funded by MCIN/AEI/10.13039/501100011033 and by “ERDF: A way of making Europe”, by the “European Union”. We also acknowledge funding from Generalitat de Catalunya, grant 2021SGR00042. The work was also funded by the European Union (ERC, NovoGenePop, project number 101052538).

The differential production of transcript isoforms from gene loci is a key cellular mechanism. Yet, its impact in protein production remains an open question. Here, we describe ORQAS (ORF quantification pipeline for alternative splicing), a pipeline for the translation quantification of individual transcript isoforms using ribosome-protected mRNA fragments (ribosome profiling). We find evidence of translation for 40–50% of the expressed isoforms in human and mouse, with 53% of the expressed genes having more than one translated isoform in human, and 33% in mouse. Differential splicing analysis revealed that about 40% of the splicing changes at RNA level are concordant with changes in translation. Furthermore, orthologous cassette exons between human and mouse preserve the directionality of the change, and are enriched in microexons in a comparison between glia and glioma. ORQAS leverages ribosome profiling to uncover a widespread and evolutionarily conserved impact of differential splicing on translation, particularly of microexon-containing isoforms., Genes express multiple mRNA isoforms through alternative processing. Here the authors analyze ribosome profiling data with ORQAS (ORF quantification pipeline for alternative splicing) and report that 40–50% of the expressed mRNA isoforms are translated.

The unicellular yeast Schizosaccharomyces pombe (fission yeast) retains many of the splicing features observed in humans and is thus an excellent model to study the basic mechanisms of splicing. Nearly half the genes contain introns, but the impact of alternative splicing in gene regulation and proteome diversification remains largely unexplored. Here we leverage Oxford Nanopore Technologies native RNA sequencing (dRNA), as well as ribosome profiling data, to uncover the full range of polyadenylated transcripts and translated open reading frames. We identify 332 alternative isoforms affecting the coding sequences of 262 different genes, 97 of which occur at frequencies higher than 20%, indicating that functional alternative splicing in S. pombe is more prevalent than previously suspected. Intron retention events make about 80% of the cases; these events may be involved in the regulation of gene expression and, in some cases, generate novel protein isoforms, as supported by ribosome profiling data in 18 of the intron retention isoforms. One example is the rpl22 gene, in which intron retention is associated with the translation of a protein of only 13 amino acids. We also find that lowly expressed transcripts tend to have longer poly(A) tails than highly expressed transcripts, highlighting an interdependence between poly(A) tail length and transcript expression level. Finally, we discover 214 novel transcripts that are not annotated, including 158 antisense transcripts, some of which also show translation evidence. The methodologies described in this work open new opportunities to study the regulation of splicing in a simple eukaryotic model. This work benefited from preliminary Nanopore RNA-seq data analyses performed by Bea Calvo and Audald Lloret-Villas, as well as discussions with Eduardo Eyras and Ivan de la Rubia. We acknowledge funding from Ministerio de Ciencia e Innovación (MCI), Agencia Estatal de Investigación (AEI) grant PGC2018–094091-B-I00, cofunded by Fondo Europeo de Desarrollo Regional (FEDER), and from Agència de Gestió d'Ajuts Universitaris i de Recerca (AGAUR), Generalitat de Catalunya, grant 2017SGR01020.

Ribosome profiling (Ribo-seq) has catalyzed a paradigm shift in our understanding of the translational ‘vocabulary’ of the human genome, discovering thousands of translated open reading frames (ORFs) within long non-coding RNAs and presumed untranslated regions of protein-coding genes. However, reference gene annotation projects have been circumspect in their incorporation of these ORFs due to uncertainties about their experimental reproducibility and physiological roles. Yet, it is indisputable that certain Ribo-seq ORFs make stable proteins, others mediate gene regulation, and many have medical implications. Ultimately, the absence of standardized ORF annotation has created a circular problem: while Ribo-seq ORFs remain unannotated by reference biological databases, this lack of characterisation will thwart research efforts examining their roles. Here, we outline the initial stages of a community-led effort supported by GENCODE / Ensembl, HGNC and UniProt to produce a consolidated catalog of human Ribo-seq ORFs.

Cells responds to diverse stimuli by changing the levels of specific effector proteins. These changes are usually examined using high throughput RNA sequencing data (RNA-Seq); transcriptional regulation is generally assumed to directly influence protein abundances. However, the correlation between RNA-Seq and proteomics data is in general quite limited owing to differences in protein stability and translational regulation. Here we perform RNA-Seq, ribosome profiling and proteomics analyses in baker's yeast cells grown in rich media and oxidative stress conditions to examine gene expression regulation at various levels. With the exception of a small set of genes involved in the maintenance of the redox state, which are regulated at the transcriptional level, modulation of protein expression is largely driven by changes in the relative ribosome density across conditions. The majority of shifts in mRNA abundance are compensated by changes in the opposite direction in the number of translating ribosomes and are predicted to result in no net change at the protein level. We also identify a subset of mRNAs which is likely to undergo specific translational repression during stress and which includes cell cycle control genes. The study suggests that post-transcriptional buffering of gene expression may be more common than previously anticipated. The work was funded by grants BFU2015–65235-P, BFU2015-68351-P and BFU2016-80039-R, from Ministerio de Economía e Innovación (Spanish Government) - FEDER (EU), and from grant PT17/0009/0014 from Instituto de Salud Carlos III – FEDER. We also received funding from the “Maria de Maeztu” Programme for Units of Excellence in R&D (MDM-2014-0370) and from Agència de Gestió d’Ajuts Universitaris i de Recerca Generalitat de Catalunya (AGAUR), grant number 2014SGR1121, 2014SGR0974, 2017SGR01020 and, predoctoral fellowship (FI) to W.B. We also acknowledge support from the EU Erasmus Programme to T.T.

The translatome can be defined as the sum of the RNA sequences that are translated into proteins in the cell by the ribosomal machinery. Until recently, it was generally assumed that the translatome was essentially restricted to evolutionary conserved proteins encoded by the set of annotated protein-coding genes. However, it has become increasingly clear that it also includes small regulatory open reading frames (ORFs), functional micropeptides, de novo proteins, and the pervasive translation of likely nonfunctional proteins. Many of these ORFs have been discovered thanks to the development of ribosome profiling, a technique to sequence ribosome-protected RNA fragments. To fully capture the diversity of translated ORFs, we propose a comprehensive classification that includes the new types of translated ORFs in addition to standard proteins.

Background A large fraction of genes contains upstream ORFs (uORFs) in the 5′ untranslated region (5’UTR). The translation of uORFs can inhibit the translation of the main coding sequence, for example by causing premature dissociation of the two ribosomal units or ribosome stalling. However, it is currently unknown if most uORFs are inhibitory or if this activity is restricted to specific cases. Here we interrogate ribosome profiling data from three different stress experiments in yeast to gain novel insights into this question. Results By comparing ribosome occupancies in different conditions and experiments we obtain strong evidence that, in comparison to primary coding sequences (CDS), which undergo translational arrest during stress, the translation of uORFs is mostly unaffected by changes in the environment. As a result, the relative abundance of uORF-encoded peptides increases during stress. In general, the changes in the translational efficiency of regions containing uORFs do not seem to affect downstream translation. The exception are uORFs found in a subset of genes that are significantly up-regulated at the level of translation during stress; these uORFs tend to be translated at lower levels in stress conditions than in optimal growth conditions, facilitating the translation of the CDS during stress. We find new examples of uORF-mediated regulation of translation, including the Gcn4 functional homologue fil1 and ubi4 genes in S. pombe. Conclusion We find evidence that the relative amount of uORF-encoded peptides increases during stress. The increased translation of uORFs is however uncoupled from the general CDS translational repression observed during stress. In a subset of genes that encode proteins that need to be rapidly synthesized upon stress uORFs act as translational switches.

Ligand-dependent corepressor (LCOR) mediates normal and malignant breast stem cell differentiation. Cancer stem cells (CSCs) generate phenotypic heterogeneity and drive therapy resistance, yet their role in immunotherapy is poorly understood. Here we show that immune-checkpoint blockade (ICB) therapy selects for LCORlow CSCs with reduced antigen processing/presentation machinery (APM) driving immune escape and ICB resistance in triple-negative breast cancer (TNBC). We unveil an unexpected function of LCOR as a master transcriptional activator of APM genes binding to IFN-stimulated response elements (ISREs) in an IFN signaling-independent manner. Through genetic modification of LCOR expression, we demonstrate its central role in modulation of tumor immunogenicity and ICB responsiveness. In TNBC, LCOR associates with ICB clinical response. Importantly, extracellular vesicle (EV) Lcor-messenger RNA therapy in combination with anti-PD-L1 overcame resistance and eradicated breast cancer metastasis in preclinical models. Collectively, these data support LCOR as a promising target for enhancement of ICB efficacy in TNBC, by boosting of tumor APM independently of IFN. This work was supported by the Instituto de Salud Carlos III-FSE (nos. MS17/00037 and PI18/00014) and the Cancer Research Institute, Clinic and Laboratory Integration Program (grant no. CRI2477 to T.C.-T). We also thank the AECC LAB (grant no. LABAE19007CELI) and the FERO foundation (to T.C.-T). This work was also supported by ISCIII (CIBERONC nos. CB16/12/00481, CB16/12/00241, PI18/00006 and PI21/00002), Generalitat de Catalunya (no. 2017 SGR 507) and the European Community through the Regional Development Funding Program to J. Albanell. Y.K. is supported by the Brewster Foundation, the Breast Cancer Research Foundation, the Susan G. Komen Foundation and the American Cancer Society. J. Arribas is funded by the Breast Cancer Research Foundation (no. BCRF-20-08), Instituto de Salud Carlos III (no. PI19/01181), Asociación Española Contra el Cáncer (no. GCAEC19017ARRI) and Fundación BBVA (no. CAIMI VHIO-FBBVA 2018-2021). M.M.A. is funded by MICINN and the Spanish Government (no. PGC2018-094091-B-I00). D.C. was funded by Instituto de Salud Carlos III (no. JR1800003). S.M. is funded by PERIS (no. SLT006/17/00040). The wotk of R.R.G. is funded by MICINN, the Spanish Government (no. PID2019-104948RB-I00) and the BBVA Foundation. S.B and J.P. acknowledge support from the Spanish Ministry of Science, the EMBL partnership and the CESO and CERCA Program.

Additional file 1: Table S1. Number of genes with higher number of reads in the 5’UTR vs CDS in stress than in normal conditions (S/N > 0) and the other way round (S/N

The aim of this study was to generate and analyze the atlas of the loggerhead turtle blood transcriptome by RNA-seq, as well as identify and characterize thioredoxin (Tnxs) and peroxiredoxin (Prdxs) antioxidant enzymes of the greatest interest in the control of peroxide levels and other biological functions. The transcriptome of loggerhead turtle was sequenced using the Illumina Hiseq 2000 platform and de novo assembly was performed using the Trinity pipeline. The assembly comprised 515,597 contigs with an N50 of 2,631 bp. Contigs were analyzed with CD-Hit obtaining 374,545 unigenes, of which 165,676 had ORFs encoding putative proteins longer than 100 amino acids. A total of 52,147 (31.5%) of these transcripts had significant homology matches in at least one of the five databases used. From the enrichment of GO terms, 180 proteins with antioxidant activity were identified, among these 28 Prdxs and 50 putative Tnxs. The putative proteins of loggerhead turtles encoded by the genes Prdx1, Prdx3, Prdx5, Prdx6, Txn and Txnip were predicted and characterized in silico. When comparing Prdxs and Txns of loggerhead turtle with homologous human proteins, they showed 18 (9%), 52 (18%) 94 (43%), 36 (16%), 35 (33%) and 74 (19%) amino acid mutations respectively. However, they showed high conservation in active sites and structural motifs (98%), with few specific modifications. Of these, Prdx1, Prdx3, Prdx5, Prdx6, Txn and Txnip presented 0, 25, 18, three, six and two deleterious changes. This study provides a high quality blood transcriptome and functional annotation of loggerhead sea turtles.

De novo gene origination has been recently established as an important mechanism for the formation of new genes. In organisms with a large genome, intergenic and intronic regions provide plenty of raw material for new transcriptional events to occur, but little is know about how de novo transcripts originate in more densely-packed genomes. Here, we identify 213 de novo originated transcripts in Saccharomyces cerevisiae using deep transcriptomics and genomic synteny information from multiple yeast species grown in two different conditions. We find that about half of the de novo transcripts are expressed from regions which already harbor other genes in the opposite orientation; these transcripts show similar expression changes in response to stress as their overlapping counterparts, and some appear to translate small proteins. Thus, a large fraction of de novo genes in yeast are likely to co-evolve with already existing genes., Genome-wide studies of de novo genes have tended to focus on genomic open reading frames (ORFs). Here, Blevins et al. use deep transcriptomics and synteny information to identify de novo transcripts in the yeast Saccharomyces cerevisiae, many of which are expressed from the alternative DNA strand.

Significance In plants, hormone signaling is, to some extent, controlled by transcriptional repressors that are degraded in a hormone-dependent manner. These repressor complexes consist of multiple proteins that assemble. MicroProteins are small single-domain proteins that are sequence-related to larger, multidomain proteins and act by disrupting protein complex formation. By forming protein complexes, or by sequestering signaling proteins, microProteins can act as effective modulators of protein activity. Our studies have discovered LITTLE NINJA, a NINJA-related microProtein that modulates JA signaling by attenuating the repression of JA-signaling. Transgenic plants ectopically expressing LITTLE NINJA are shorter and bushier than wild-type plants. Modulating LITTLE NINJA activity may provide strategies for future crop improvement and contribute to finding sustainable solutions for future agricultural production., MicroProteins are small, often single-domain proteins that are sequence-related to larger, often multidomain proteins. Here, we used a combination of comparative genomics and heterologous synthetic misexpression to isolate functional cereal microProtein regulators. Our approach identified LITTLE NINJA (LNJ), a microProtein that acts as a modulator of jasmonic acid (JA) signaling. Ectopic expression of LNJ in Arabidopsis resulted in stunted plants that resembled the decuple JAZ (jazD) mutant. In fact, comparing the transcriptomes of transgenic LNJ overexpressor plants and jazD revealed a large overlap of deregulated genes, suggesting that ectopic LNJ expression altered JA signaling. Transgenic Brachypodium plants with elevated LNJ expression levels showed deregulation of JA signaling as well and displayed reduced growth and enhanced production of side shoots (tiller). This tillering effect was transferable between grass species, and overexpression of LNJ in barley and rice caused similar traits. We used a clustered regularly interspaced short palindromic repeats (CRISPR) approach and created a LNJ-like protein in Arabidopsis by deleting parts of the coding sentence of the AFP2 gene that encodes a NINJA-domain protein. These afp2-crispr mutants were also stunted in size and resembled jazD. Thus, similar genome-engineering approaches can be exploited as a future tool to create LNJ proteins and produce cereals with altered architectures.

Hibernation is an adaptive strategy some mammals use to survive highly seasonal or unpredictable environments. We present the first investigation on the transcriptomics of hibernation in a natural population of primate hibernators: Crossley's dwarf lemurs (Cheirogaleus crossleyi). Using capture-mark-recapture techniques to track the same animals over a period of 7 months in Madagascar, we used RNA-seq to compare gene expression profiles in white adipose tissue (WAT) during three distinct physiological states. We focus on pathway analysis to assess the biological significance of transcriptional changes in dwarf lemur WAT and, by comparing and contrasting what is known in other model hibernating species, contribute to a broader understanding of genomic contributions of hibernation across Mammalia. The hibernation signature is characterized by a suppression of lipid biosynthesis, pyruvate metabolism and mitochondrial-associated functions, and an accumulation of transcripts encoding ribosomal components and iron-storage proteins. The data support a key role of pyruvate dehydrogenase kinase isoenzyme 4 (PDK4) in regulating the shift in fuel economy during periods of severe food deprivation. This pattern of PDK4 holds true across representative hibernating species from disparate mammalian groups, suggesting that the genetic underpinnings of hibernation may be ancestral to mammals.

High throughput RNA sequencing techniques have revealed that a large fraction of the genome is transcribed into long non-coding RNAs (lncRNAs). Unlike canonical protein-coding genes, lncRNAs do not contain long open reading frames (ORFs) and tend to be poorly conserved across species. However, many of them contain small ORFs (sORFs) that exhibit translation signatures according to ribosome profiling or proteomics data. These sORFs are a source of putative novel proteins; some of them may confer a selective advantage and be maintained over time, a process known as de novo gene birth. Here we review the mechanisms by which randomly occurring sORFs in lncRNAs can become new functional proteins.

The differential production of transcript isoforms from gene loci is a key cellular mechanism. Yet, its impact in protein production remains an open question. Here, we describe ORQAS (ORF quantification pipeline for alternative splicing), a new pipeline for the translation quantification of individual transcript isoforms using ribosome-protected mRNA fragments (Ribosome profiling). We found evidence of translation for 40-50% of the expressed transcript isoforms in human and mouse, with 53% of the expressed genes having more than one translated isoform in human, 33% in mouse. Differential analysis revealed that about 40% of the splicing changes at RNA level were concordant with changes in translation, with 21.7% of changes at RNA level and 17.8% at translation level conserved between human and mouse. Furthermore, orthologous cassette exons preserving the directionality of the change were conserved between human and mouse and enriched in microexons in a comparison between glia and glioma. ORQAS leverages ribosome profiling to uncover a widespread and evolutionary conserved impact of differential splicing on the translation of isoforms and, in particular, of microexon-containing isoforms. ORQAS is available at https://github.com/comprna/orqas

Evidence has accumulated that some genes originate directly from previously non-genic sequences, orde novo, rather than by the duplication or fusion of existing genes. However, howde novogenes emerge and eventually become functional is largely unknown. Here we perform the first study onde novogenes that uses transcriptomics data from eleven different yeast species, all grown identically in both rich media and in oxidative stress conditions. The genomes of these species are densely-packed with functional elements, leaving little room for the co-option of genomic sequences into new transcribed loci. Despite this, we find that at least 213 transcripts (~5%) have arisende novoin the past 20 million years of evolution of baker’s yeast-or approximately 10 new transcripts every million years. Nearly half of the total newly expressed sequences are generated from regions in which both DNA strands are used as templates for transcription, explaining the apparent contradiction between the limited ‘empty’ genomic space and high rate ofde novogene birth. In addition, we find that 40% of thesede novotranscripts are actively translated and that at least a fraction of the encoded proteins are likely to be under purifying selection. This study shows that even in very highly compact genomes,de novotranscripts are continuously generated and can give rise to new functional protein-coding genes.

Cells respond to changes in the environment by modifying the concentration of specific proteins. Paradoxically, the cellular response is usually examined by measuring variations in transcript abundance by high throughput RNA sequencing (RNA-Seq), instead of directly measuring protein concentrations. This happens because RNA-Seq-based methods provide better quantitative estimates, and more extensive gene coverage, than proteomics-based ones. However, variations in transcript abundance do not necessarily reflect changes in the corresponding protein abundance. How can we close this gap? Here we explore the use of ribosome profiling (Ribo-Seq) to perform differentially gene expression analysis in a relatively well-characterized system, oxidative stress in baker’s yeast. Ribo-Seq is an RNA sequencing method that specifically targets ribosome-protected RNA fragments, and thus is expected to provide a more accurate view of changes at the protein level than classical RNA-Seq. We show that gene quantification by Ribo-Seq is indeed more highly correlated with protein abundance, as measured from mass spectrometry data, than quantification by RNA-Seq. The analysis indicates that, whereas a subset of genes involved in oxidation-reduction processes is detected by both types of data, the majority of the genes that happen to be significant in the RNA-Seq-based analysis are not significant in the Ribo-Seq analysis, suggesting that they do not result in protein level changes. The results illustrate the advantages of Ribo-Seq to make inferences about changes in protein abundance in comparison with RNA-Seq.

SUMMARYThe mammalian transcriptome includes thousands of transcripts that do not correspond to annotated protein-coding genes. Although many of these transcripts show homology between human and mouse, only a small proportion of them have been functionally characterized. Here we use ribosome profiling data to identify translated open reading frames, as well as non-ribosomal protein-RNA interactions, in evolutionary conserved and non-conserved transcripts. We find that conserved regions are subject to significant evolutionary constraints and are enriched in translated open reading frames, as well as non-ribosomal protein-RNA interaction signatures, when compared to non-conserved regions. Translated ORFs can be divided in two classes, those encoding functional micropeptides and those that show no evidence of protein functionality. This study underscores the importance of combining evolutionary and biochemical measurements to advance in a more complete understanding of the transcriptome.

Accumulating evidence indicates that some protein-coding genes have originated de novo from previously non-coding genomic sequences. However, the processes underlying de novo gene birth are still enigmatic. In particular, the appearance of a new functional protein seems highly improbable unless there is already a pool of neutrally evolving peptides that are translated at significant levels and that can at some point acquire new functions. Here, we use deep ribosome-profiling sequencing data, together with proteomics and single nucleotide polymorphism information, to search for these peptides. We find hundreds of open reading frames that are translated and that show no evolutionary conservation or selective constraints. These data suggest that the translation of these neutrally evolving peptides may be facilitated by the chance occurrence of open reading frames with a favourable codon composition. We conclude that the pervasive translation of the transcriptome provides plenty of material for the evolution of new functional proteins. We are grateful for valuable discussions with many colleagues during this study. This work was funded by grants BFU2012-36820, BFU2015-65235-P and TIN2015-69175-C4-3-R from Ministerio de Economía e Innovación (Spanish Government) and co-funded by FEDER (EC). We also received funding from Agència de Gestió d’Ajuts Universitaris i de Recerca Generalitat de Catalunya (AGAUR), grant no. 2014SGR1121.

The birth of genes that encode new protein sequences is a major source of evolutionary innovation. However, we still understand relatively little about how these genes come into being and which functions they are selected for. To address these questions, we have obtained a large collection of mammalian-specific gene families that lack homologues in other eukaryotic groups. We have combined gene annotations and de novo transcript assemblies from 30 different mammalian species, obtaining ∼6,000 gene families. In general, the proteins in mammalian-specific gene families tend to be short and depleted in aromatic and negatively charged residues. Proteins which arose early in mammalian evolution include milk and skin polypeptides, immune response components, and proteins involved in reproduction. In contrast, the functions of proteins which have a more recent origin remain largely unknown, despite the fact that these proteins also have extensive proteomics support. We identify several previously described cases of genes originated de novo from noncoding genomic regions, supporting the idea that this mechanism frequently underlies the evolution of new protein-coding genes in mammals. Finally, we show that most young mammalian genes are preferentially expressed in testis, suggesting that sexual selection plays an important role in the emergence of new functional genes. The work was funded by grants BFU2012-36820 and BFU2015-65235-P from Ministerio de Economía e Innovación (Spanish Government) and co-funded by FEDER. We also received funding from Agència de Gestió d’Ajuts Universitaris i de Recerca Generatlitat de Catalunya (AGAUR), grant number 2014SGR1121

Recent years have witnessed the discovery of protein–coding genes which appear to have evolved de novo from previously non-coding sequences. This has changed the long-standing view that coding sequences can only evolve from other coding sequences. However, there are still many open questions regarding how new protein-coding sequences can arise from non-genic DNA. Two prerequisites for the birth of a new functional protein-coding gene are that the corresponding DNA fragment is transcribed and that it is also translated. Transcription is known to be pervasive in the genome, producing a large number of transcripts that do not correspond to conserved protein-coding genes, and which are usually annotated as long non-coding RNAs (lncRNA). Recently, sequencing of ribosome protected fragments (Ribo-Seq) has provided evidence that many of these transcripts actually translate small proteins. We have used mouse non-synonymous and synonymous variation data to estimate the strength of purifying selection acting on the translated open reading frames (ORFs). Whereas a subset of the lncRNAs are likely to actually be true protein-coding genes (and thus previously misclassified), the bulk of lncRNAs code for proteins which show variation patterns consistent with neutral evolution. We also show that the ORFs that have a more favorable, coding-like, sequence composition are more likely to be translated than other ORFs in lncRNAs. This study provides strong evidence that there is a large and ever-changing reservoir of lowly abundant proteins; some of these peptides may become useful and act as seeds for de novo gene evolution.

Recent years have witnessed the discovery of protein–coding genes which appear to have evolved de novo from previously non-coding sequences. This has changed the long-standing view that coding sequences can only evolve from other coding sequences. However, there are still many open questions regarding how new protein-coding sequences can arise from non-genic DNA. Two prerequisites for the birth of a new functional protein-coding gene are that the corresponding DNA fragment is transcribed and that it is also translated. Transcription is known to be pervasive in the genome, producing a large number of transcripts that do not correspond to conserved protein-coding genes, and which are usually annotated as long non-coding RNAs (lncRNA). Recently, sequencing of ribosome protected fragments (Ribo-Seq) has provided evidence that many of these transcripts actually translate small proteins. We have used mouse non-synonymous and synonymous variation data to estimate the strength of purifying selection acting on the translated open reading frames (ORFs). Whereas a subset of the lncRNAs are likely to actually be true protein-coding genes (and thus previously misclassified), the bulk of lncRNAs code for proteins which show variation patterns consistent with neutral evolution. We also show that the ORFs that have a more favorable, coding-like, sequence composition are more likely to be translated than other ORFs in lncRNAs. This study provides strong evidence that there is a large and ever-changing reservoir of lowly abundant proteins; some of these peptides may become useful and act as seeds for de novo gene evolution.

Phylostratigraphy is a computational framework for dating the emergence of sequences (usually genes) in a phylogeny. It has been extensively applied to make inferences on patterns of genome evolution, including patterns of disease gene evolution, ontogeny and de novo gene origination. Phylostratigraphy typically relies on BLAST searches along a species tree, but new simulation studies have raised concerns about the ability of BLAST to detect remote homologues and its impact on phylostratigraphic inferences. These simulations called into question some of our previously published work on patterns of gene emergence and evolution inferred from phylostratigraphy. Here, we re-assessed these simulations and found major problems including unrealistic parameter choices, irreproducibility, statistical flaws and partial representation of results. We found that, even with a possible overall BLAST false negative rate between 5-15%, the large majority (>74%) of sequences assigned to a recent evolutionary origin by phylostratigraphy is unaffected by technical concerns about BLAST. Where the results of the simulations did cast doubt on our previous findings, we repeated our analyses but now excluded all questionable sequences. The originally described patterns remained essentially unchanged. These new analyses strongly support our published inferences, including: genes that emerged after the origin of eukaryotes are more likely to be expressed in the ectoderm than in the endoderm or mesoderm in Drosophila, and the de novo emergence of protein-coding genes from non-genic sequences occurs through proto-gene intermediates in yeast. We conclude that BLAST is an appropriate and sufficiently sensitive tool in phylostratigraphic analysis.

In recent years, the study of the molecular processes involved in mammalian hibernation has shifted from investigating a few carefully selected candidate genes to large-scale analysis of differential gene expression. The availability of high-throughput data provides an unprecedented opportunity to ask whether phylogenetically distant species show similar mechanisms of genetic control, and how these relate to particular genes and pathways involved in the hibernation phenotype. In order to address these questions, we compare 11 datasets of differentially expressed (DE) genes from two ground squirrel species, one bat species, and the American black bear, as well as a list of genes extracted from the literature that previously have been correlated with the drastic physiological changes associated with hibernation. We identify several genes that are DE in different species, indicating either ancestral adaptations or evolutionary convergence. When we use a network approach to expand the original datasets of DE genes to large gene networks using available interactome data, a higher agreement between datasets is achieved. This indicates that the same key pathways are important for activating and maintaining the hibernation phenotype. Functional-term-enrichment analysis identifies several important metabolic and mitochondrial processes that are critical for hibernation, such as fatty acid beta-oxidation and mitochondrial transport. We do not detect any enrichment of positive selection signatures in the coding sequences of genes from the networks of hibernation-associated genes, supporting the hypothesis that the genetic processes shaping the hibernation phenotype are driven primarily by changes in gene regulation.

The mammalian transcriptome includes thousands of transcripts that do not correspond to annotated protein-coding genes and that are known as long non-coding RNAs (lncRNAs). A handful of lncRNAs have well-characterized regulatory functions but the biological significance of the majority of them is not well understood. LncRNAs that are conserved between mice and humans are likely to be enriched in functional sequences. Here, we investigate the presence of different types of ribosome profiling signatures in lncRNAs and how they relate to sequence conservation. We find that lncRNA-conserved regions contain three times more ORFs with translation evidence than non-conserved ones, and identify nine cases that display significant sequence constraints at the amino acid sequence level. The study also reveals that conserved regions in intergenic lncRNAs are significantly enriched in protein–RNA interaction signatures when compared to non-conserved ones; this includes sites in well-characterized lncRNAs, such as Cyrano, Malat1, Neat1 and Meg3, as well as in tens of lncRNAs of unknown function. This work illustrates how the analysis of ribosome profiling data coupled with evolutionary analysis provides new opportunities to explore the lncRNA functional landscape.

Gene duplication is widely regarded as a major mechanism modeling genome evolution and function. However, the mechanisms that drive the evolution of the two, initially redundant, gene copies are still ill defined. Many gene duplicates experience evolutionary rate acceleration, but the relative contribution of positive selection and random drift to the retention and subsequent evolution of gene duplicates, and for how long the molecular clock may be distorted by these processes, remains unclear. Focusing on rodent genes that duplicated before and after the mouse and rat split, we find significantly increased sequence divergence after duplication in only one of the copies, which in nearly all cases corresponds to the novel daughter copy, independent of the mechanism of duplication. We observe that the evolutionary rate of the accelerated copy, measured as the ratio of nonsynonymous to synonymous substitutions, is on average 5-fold higher in the period spanning 4-12 My after the duplication than it was before the duplication. This increase can be explained, at least in part, by the action of positive selection according to the results of the maximum likelihood-based branch-site test. Subsequently, the rate decelerates until purifying selection completely returns to preduplication levels. Reversion to the original rates has already been accomplished 40.5 My after the duplication event, corresponding to a genetic distance of about 0.28 synonymous substitutions per site. Differences in tissue gene expression patterns parallel those of substitution rates, reinforcing the role of neofunctionalization in explaining the evolution of young gene duplicates.

Cells express thousands of transcripts that show weak coding potential. Known as long non-coding RNAs (lncRNAs), they typically contain short open reading frames (ORFs) having no homology with known proteins. Recent studies show that a significant proportion of lncRNAs are translated, challenging the view that they are non-coding. These results are based on selective sequencing of ribosome-protected fragments, or ribosome profiling. The present study used ribosome profiling data from eight mouse tissues and cell types, combined with ~330,000 synonymous and non-synonymous single nucleotide variants, to dissect the patterns of purifying selection in proteins translated from lncRNAs. Using the three-nucleotide read periodicity that characterizes actively translated regions, we identified 832 mouse translated lncRNAs. Overall, they produced 1,489 different proteins, most of them smaller than 100 amino acids. Nearly half of the ORFs then showed sequence conservation in rat and/or human transcripts, and many of them are likely to encode functional micropeptides, including the recently discovered Myoregulin. For lncRNAs not conserved in rats or humans, the ORF codon usage bias distinguished between two classes, one with particularly high coding scores and evidence of purifying selection, consistent with the presence of lineage-specific functional proteins, and a second, larger, class of ORFs producing peptides with no significant purifying selection signatures. We obtained evidence that the translation of these lncRNAs depends on the chance occurrence of ORFs with a favorable codon composition. Some of these lncRNAs may be precursors of novel protein-coding genes, filling a gap in our current understanding of de novo gene birth.

Hibernation is a complex physiological response that some mammalian species employ to evade energetic demands. Previous work in mammalian hibernators suggests that hibernation is activated not by a set of genes unique to hibernators, but by differential expression of genes that are present in all mammals. This question of universal genetic mechanisms requires further investigation and can only be tested through additional investigations of phylogenetically dispersed species. To explore this question, we use RNA-Seq to investigate gene expression dynamics as they relate to the varying physiological states experienced throughout the year in a group of primate hibernators—Madagascar’s dwarf lemurs (genus Cheirogaleus). In a novel experimental approach, we use longitudinal sampling of biological tissues as a method for capturing gene expression profiles from the same individuals throughout their annual hibernation cycle. We identify 90 candidate genes that have variable expression patterns when comparing two active states (Active 1 and Active 2) with a torpor state. These include genes that are involved in metabolic pathways, feeding behavior, and circadian rhythms, as might be expected to correlate with seasonal physiological state changes. The identified genes appear to be critical for maintaining the health of an animal that undergoes prolonged periods of metabolic depression concurrent with the hibernation phenotype. By focusing on these differentially expressed genes in dwarf lemurs, we compare gene expression patterns in previously studied mammalian hibernators. Additionally, by employing evolutionary rate analysis, we find that hibernation-related genes do not evolve under positive selection in hibernating species relative to nonhibernators., This work was supported by grants from the Society of Integrative and Comparative Biology and Sigma Xi to S.L.F., a National Science Foundation Doctoral Research Improvement Grant (NSF-BCS 1455809 to A.D.Y. and S.L.F.), a grant from the Duke University Institute for Genome Science and Policy to A.D.Y., and funding from the Ministerio de Economía y Competitividad of the Spanish Government (BIO2009-08160 to J.L.V.C. and BFU2012-36820 to J.L.V.C. and M.M.A.); and Fundació ICREA to M.M.A.

We have investigated the previously uncharacterized human cytomegalovirus (HCMV) UL1 open reading frame (ORF), a member of the rapidly evolving HCMV RL11 family. UL1 is HCMV specific; the absence of UL1 in chimpanzee cytomegalovirus (CCMV) and sequence analysis studies suggest that UL1 may have originated by the duplication of an ancestor gene from the RL11-TRL cluster ( TRL11 , TRL12 , and TRL13 ). Sequence similarity searches against human immunoglobulin (Ig)-containing proteins revealed that HCMV pUL1 shows significant similarity to the cellular carcinoembryonic antigen-related (CEA) protein family N-terminal Ig domain, which is responsible for CEA ligand recognition. Northern blot analysis revealed that UL1 is transcribed during the late phase of the viral replication cycle in both fibroblast-adapted and endotheliotropic strains of HCMV. We characterized the protein encoded by hemagglutinin (HA)-tagged UL1 in the AD169-derived HB5 background. UL1 is expressed as a 224-amino-acid type I transmembrane glycoprotein which becomes detectable at 48 h postinfection. In infected human fibroblasts, pUL1 colocalized at the cytoplasmic site of virion assembly and secondary envelopment together with TGN-46, a marker for the trans -Golgi network, and viral structural proteins, including the envelope glycoprotein gB and the tegument phosphoprotein pp28. Furthermore, analyses of highly purified AD169 UL1-HA epitope-tagged virions revealed that pUL1 is a novel constituent of the HCMV envelope. Importantly, the deletion of UL1 in HCMV TB40/E resulted in reduced growth in a cell type-specific manner, suggesting that pUL1 may be implicated in regulating HCMV cell tropism.

Human cytomegalovirus (HCMV), the β-herpesvirus prototype, has evolved a wide spectrum of mechanisms to counteract host immunity. Among them, HCMV uses cellular captured genes encoding molecules capable of interfering with the original host function or of fulfilling new immunomodulatory tasks. Here, we report on UL7, a novel HCMV heavily glycosylated transmembrane protein, containing an Ig-like domain that exhibits remarkable amino acid similarity to CD229, a cell-surface molecule of the signalling lymphocyte-activation molecule (SLAM) family involved in leukocyte activation. The UL7 Ig-like domain, which is well-preserved in all HCMV strains, structurally resembles the SLAM-family N-terminal Ig-variable domain responsible for the homophilic and heterophilic interactions that trigger signalling. UL7 is transcribed with early-late kinetics during the lytic infectious cycle. Using a mAb generated against the viral protein, we show that it is constitutively shed, through its mucine-like stalk, from the cell-surface. Production of soluble UL7 is enhanced by PMA and reduced by a broad-spectrum metalloproteinase inhibitor. Although UL7 does not hold the ability to interact with CD229 or other SLAM-family members, it shares with them the capacity to mediate adhesion to leukocytes, specifically to monocyte-derived DCs. Furthermore, we demonstrate that UL7 expression attenuates the production of proinflammatory cytokines TNF, IL-8 and IL-6 in DCs and myeloid cell lines. Thus, the ability of UL7 to interfere with cellular proinflammatory responses may contribute to viral persistence. These results enhance our understanding of those HCMV-encoded molecules involved in sustaining the balance between HCMV and the host immune system.

The molecular clock hypothesis states that protein-coding genes evolve at an approximately constant rate. However, this is only expected to be true as long as the function and the tertiary structure of the molecule remain unaltered. An important implication of this statement is that significant deviations in the rate of evolution of a gene with respect to the species clock are likely to reflect functional and/or structural alterations. Here, we present a method to identify such deviations and apply it to a data set of 2,929 high-quality coding sequence alignments corresponding to one-to-one orthologous genes from six mammalian species--human, macaque, mouse, rat, cow, and dog. Deviated branches are defined as those that present significant alterations in both the rate of nonsynonymous substitutions (dN) and the selective pressure (dN/dS). Strikingly, we find that as many as 24.5% of the genes show branch-specific deviations in dN and dN/dS, though this is a relatively well-conserved set of genes. Around half of these genes show branch-specific acceleration of evolutionary rates. Positive selection (PS) tests based on divergence data only identify 17.7% of the accelerated branches. Failure to identify PS in accelerated branches with an excess of radical amino acid replacements suggests that these tests are conservative. Interestingly, genes with accelerated branches are significantly enriched in neural proteins, indicating that this type of protein might play a more important role than previously thought in species diversification, although they are generally not detected by PS tests. We discuss in detail several examples of genes that show lineage-specific evolutionary rate acceleration and are involved in synaptic transmission, chemosensory perception, and ubiquitination.

Notch has been linked to β-catenin-dependent tumorigenesis; however, the mechanisms leading to Notch activation and the contribution of the Notch pathway to colorectal cancer is not yet understood. By microarray analysis, we have identified a group of genes downstream of Wnt/β-catenin (down-regulated when blocking Wnt/β-catenin) that are directly regulated by Notch (repressed by γ-secretase inhibitors and up-regulated by active Notch1 in the absence of β-catenin signaling). We demonstrate that Notch is downstream of Wnt in colorectal cancer cells through β-catenin-mediated transcriptional activation of the Notch-ligand Jagged1. Consistently, expression of activated Notch1 partially reverts the effects of blocking Wnt/β-catenin pathway in tumors implanted s.c. in nude mice. Crossing APC Min /+ with Jagged1 +/Δ mice is sufficient to significantly reduce the size of the polyps arising in the APC mutant background indicating that Notch is an essential modulator of tumorigenesis induced by nuclear β-catenin. We show that this mechanism is operating in human tumors from Familial Adenomatous Polyposis patients. We conclude that Notch activation, accomplished by β-catenin-mediated up-regulation of Jagged1, is required for tumorigenesis in the intestine. The Notch-specific genetic signature is sufficient to block differentiation and promote vasculogenesis in tumors whereas proliferation depends on both pathways.

Genomes contain a large number of genes that do not have recognizable homologues in other species and that are likely to be involved in important species-specific adaptive processes. The origin of many such "orphan" genes remains unknown. Here we present the first systematic study of the characteristics and mechanisms of formation of primate-specific orphan genes. We determine that codon usage values for most orphan genes fall within the bulk of the codon usage distribution of bona fide human proteins, supporting their current protein-coding annotation. We also show that primate orphan genes display distinctive features in relation to genes of wider phylogenetic distribution: higher tissue specificity, more rapid evolution, and shorter peptide size. We estimate that around 24% are highly divergent members of mammalian protein families. Interestingly, around 53% of the orphan genes contain sequences derived from transposable elements (TEs) and are mostly located in primate-specific genomic regions. This indicates frequent recruitment of TEs as part of novel genes. Finally, we also obtain evidence that a small fraction of primate orphan genes, around 5.5%, might have originated de novo from mammalian noncoding genomic regions.

Single-amino-acid tandem repeats are very common in mammalian proteins but their function and evolution are still poorly understood. Here we investigate how the variability and prevalence of amino acid repeats are related to the evolutionary constraints operating on the proteins. We find a significant positive correlation between repeat size difference and protein nonsynonymous substitution rate in human and mouse orthologous genes. This association is observed for all the common amino acid repeat types and indicates that rapid diversification of repeat structures, involving both trinucleotide slippage and nucleotide substitutions, preferentially occurs in proteins subject to low selective constraints. However, strikingly, we also observe a significant negative correlation between the number of repeats in a protein and the gene nonsynonymous substitution rate, particularly for glutamine, glycine, and alanine repeats. This implies that proteins subject to strong selective constraints tend to contain an unexpectedly high number of repeats, which tend to be well conserved between the two species. This is consistent with a role for selection in the maintenance of a significant number of repeats. Analysis of the codon structure of the sequences encoding the repeats shows that codon purity is associated with high repeat size interspecific variability. Interestingly, polyalanine and polyglutamine repeats associated with disease show very distinctive features regarding the degree of repeat conservation and the protein sequence selective constraints.

Background The high density of tandem repeat sequences (satellites) in nematode genomes and the availability of genome sequences from several species in the group offer a unique opportunity to better understand the evolutionary dynamics and the functional role of these sequences. We take advantage of the previously developed SATFIND program to study the satellites in four Caenorhabditis species and investigate these questions. Methods The identification and comparison of satellites is carried out in three steps. First we find all the satellites present in each species with the SATFIND program. Each satellite is defined by its length, number of repeats, and repeat sequence. Only satellites with at least ten repeats are considered. In the second step we build satellite families with a newly developed alignment program. Satellite families are defined by a consensus sequence and the number of satellites in the family. Finally we compare the consensus sequence of satellite families in different species. Results We give a catalog of individual satellites in each species. We have also identified satellite families with a related sequence and compare them in different species. We analyze the turnover of satellites: they increased in size through duplications of fragments of 100-300 bases. It appears that in many cases they have undergone an explosive expansion. In C. elegans we have identified a subset of large satellites that have strong affinity for the centromere protein CENP-A. We have also compared our results with those obtained from other species, including one nematode and three mammals. Conclusions Most satellite families found in Caenorhabditis are species-specific; in particular those with long repeats. A subset of these satellites may facilitate the formation of kinetochores in mitosis. Other satellite families in C. elegans are either related to Helitron transposons or to meiotic pairing centers.

The high regulatory complexity of vertebrates has been related to two rounds of whole genome duplication (2R-WGD) that occurred before the divergence of the major vertebrate groups. Following these events, many developmental transcription factors (TFs) were retained in multiple copies and subsequently specialized in diverse functions, whereas others reverted to their singleton state. TFs are known to be generally rich in amino acid repeats or low-complexity regions (LCRs), such as polyalanine or polyglutamine runs, which can evolve rapidly and potentially influence the transcriptional activity of the protein. Here we test the hypothesis that LCRs have played a major role in the diversification of TF gene duplicates. We find that nearly half of the TF gene families originated during the 2R-WGD contains LCRs. The number of gene duplicates with LCRs is 155 out of 550 analyzed (28%), about twice as many as the number of single copy genes with LCRs (15 out of 115, 13%). In addition, duplicated TFs preferentially accumulate certain LCR types, the most prominent of which are alanine repeats. We experimentally test the role of alanine-rich LCRs in two different TF gene families, PHOX2A/PHOX2B and LHX2/LHX9. In both cases, the presence of the alanine-rich LCR in one of the copies (PHOX2B and LHX2) significantly increases the capacity of the TF to activate transcription. Taken together, the results provide strong evidence that LCRs are important driving forces of evolutionary change in duplicated genes. The work was supported by grants from the Ministerio de Innovacion y Tecnolog ıa (BIO2009-08160 to M.M.A.), Ministerio de Economıa y Competitividad (BFU2012-36820 co-funded by FEDER to M.M.A. and BFU2010-15347 and Centro de Excelencia Severo Ochoa 2013-2017-SEV-2012- 0208 to S.L.), the Secretariat of Universities and ResearchGovernment of Catalonia (2014SGR1121 to M.M.A. and 2014SGR674 to S.L.), Fundacio Javier Lamas Miralles (fellow- ship to N.R-T.), and Institucio Catalana de Recerca i Estudis Avanc¸ats (ICREA contract to M.M.A. and S.L.).