Evidence that alternative transcriptional initiation is largely nonadaptive

mRNA Transcript Variants Expressed in Mammalian Cells

Gene expression or gene regulation studies often assume one gene expresses one mRNA. However, contrary to the conventional idea, a single gene in mammalian cells can express multiple transcript variants translated into several different proteins. The transcript variants are generated through transcription from alternative start sites and alternative post-transcriptional processing of the precursor mRNA (pre-mRNA). In addition, gene mutations and RNA editing further enhance the diversity of the transcript variants. The transcript variants can encode proteins with various domains, expanding the functional repertoire of a single gene. Some transcript variants may not encode proteins but function as non-coding RNAs and regulate gene expression. The expression level of the transcript variants may vary between cell types or within the same cells under different biological conditions. Transcript variants are characteristic of cell differentiation in a particular tissue, and the variants may play a key role in normal development and aging. Studies also reported that some transcript variants may have roles in disease pathogenesis. The biological significances urge studying the complexity of gene expression at the transcript level. This article updates the molecular basis of transcript variants in mammalian cells, including the formation mechanisms and potential roles in host biology. Gaining insight into the transcript variants will not only identify novel mechanisms of gene regulation but also unravel the role of the variants in health and disease.

Piwi regulates the usage of alternative transcription start sites in the Drosophila ovary

Alternative transcription initiation, which refers to the transcription of a gene from different transcription start sites (TSSs), is prevalent across metazoans and has important biological functions. Although transcriptional regulation has been extensively studied, the mechanism that selects one TSS over others within a gene remains elusive. Using the Cap Analysis of Gene Expression sequencing (CAGE-seq) method, we discovered that Piwi, an RNA-binding protein, regulates TSS usage in at least 87 genes. In piwi-deficient Drosophila ovaries, these genes displayed significantly altered TSS usage (ATU). The regulation of TSS usage occurred in both germline and somatic cells in ovaries, as well as in cultured ovarian somatic cells (OSCs). Correspondingly, RNA Polymerase II (Pol II) initiation and elongation at the TSSs of ATU genes were affected in germline-piwi-knockdown ovaries and piwi-knockdown OSCs. Furthermore, we identified a Facilitates Chromatin Transcription (FACT) complex component, Ssrp, that is essential for mRNA elongation, as a novel interactor of Piwi in the nucleus. Temporally controlled knockdown of ssrp affected TSS usage in ATU genes, whereas overexpression of ssrp partially rescued the TSS usage of ATU genes in piwi mutant ovaries. Thus, Piwi may interact with Ssrp to regulate TSS usage in Drosophila ovaries by affecting Pol II initiation and elongation.

Most m5C Modifications in Mammalian mRNAs are Nonadaptive

5-Methylation (m5C) on mRNA molecules is a prevalent internal posttranscriptional modification in eukaryotes. Although m5C modification has been reported to regulate some biological processes, whether most mRNA m5C modifications are functional is unknown. To address this question, we analyzed the genome-wide evolutionary characteristics of m5C modifications in protein-coding genes of humans and mice. Our analysis of RNA sequencing data from 13 tissues of both species revealed that (i) the occurrence of m5C modification is exceedingly low, (ii) the fraction of m5Cs decreases with the amount of Cs across genes or tissues, (iii) m5C modifications are mostly unshared between species, and (iv) m5C sites and motifs do not exhibit greater evolutionary conservation. Additionally, we estimate that a large fraction of the observed mRNA m5C modifications may be deleterious. Together, these observations suggest that most m5C modifications in mammalian mRNAs are nonadaptive, which has important implications for understanding the biological significance of m5C and other posttranscriptional modifications.

Re-appraising the evidence for the source, regulation and function of p53-family isoforms

The p53 family of proteins evolved from a common ancestor into three separate genes encoding proteins that act as transcription factors with distinct cellular roles. Isoforms of each member that lack specific regions or domains are suggested to result from alternative transcription start sites, alternative splicing or alternative translation initiation, and have the potential to exponentially increase the functional repertoire of each gene. However, evidence supporting the presence of individual protein variants at functional levels is often limited and is inferred by mRNA detection using highly sensitive amplification techniques. We provide a critical appraisal of the current evidence for the origins, expression, functions and regulation of p53-family isoforms. We conclude that despite the wealth of publications, several putative isoforms remain poorly established. Future research with improved technical approaches and the generation of isoform-specific protein detection reagents is required to establish the physiological relevance of p53-family isoforms in health and disease. In addition, our analyses suggest that p53-family variants evolved partly through convergent rather than divergent evolution from the ancestral gene.

A compendium of genetic variations associated with promoter usage across 49 human tissues

Promoters play a crucial role in regulating gene transcription. However, our understanding of how genetic variants influence alternative promoter selection is still incomplete. In this study, we implement a framework to identify genetic variants that affect the relative usage of alternative promoters, known as promoter usage quantitative trait loci (puQTLs). By constructing an atlas of human puQTLs across 49 different tissues from 838 individuals, we have identified approximately 76,856 independent loci associated with promoter usage, encompassing 602,009 genetic variants. Our study demonstrates that puQTLs represent a distinct type of molecular quantitative trait loci, effectively uncovering regulatory targets and patterns. Furthermore, puQTLs are regulating in a tissue-specific manner and are enriched with binding sites of epigenetic marks and transcription factors, especially those involved in chromatin architecture formation. Notably, we have also found that puQTLs colocalize with complex traits or diseases and contribute to their heritability. Collectively, our findings underscore the significant role of puQTLs in elucidating the molecular mechanisms underlying tissue development and complex diseases.

Evolution and Spatiotemporal Expression of ankha and ankhb in Zebrafish

Craniometaphyseal Dysplasia (CMD) is a rare skeletal disorder that can result from mutations in the ANKH gene. This gene encodes progressive anksylosis (ANK), which is responsible for transporting inorganic pyrophosphate (PPi) and ATP from the intracellular to the extracellular environment, where PPi inhibits bone mineralization. When ANK is dysfunctional, as in patients with CMD, the passage of PPi to the extracellular environment is reduced, leading to excess mineralization, particularly in bones of the skull. Zebrafish may serve as a promising model to study the mechanistic basis of CMD. Here, we provide a detailed analysis of the zebrafish Ankh paralogs, Ankha and Ankhb, in terms of their phylogenic relationship with ANK in other vertebrates as well as their spatiotemporal expression patterns during zebrafish development. We found that a closer evolutionary relationship exists between the zebrafish Ankhb protein and its human and other vertebrate counterparts, and stronger promoter activity was predicted for ankhb compared to ankha. Furthermore, we noted distinct temporal expression patterns, with ankha more prominently expressed in early development stages, and both paralogs also being expressed at larval growth stages. Whole-mount in situ hybridization was used to compare the spatial expression patterns of each paralog during bone development, and both showed strong expression in the craniofacial region as well as the notochord and somites. Given the substantial overlap in spatiotemporal expression but only subtle patterning differences, the exact roles of these genes remain speculative. In silico analyses predicted that Ankha and Ankhb have the same function in transporting PPi across the membrane. Nevertheless, this study lays the groundwork for functional analyses of each ankh paralog and highlights the potential of using zebrafish to find possible targeted therapies for CMD.

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

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

GTDrift: a resource for exploring the interplay between genetic drift, genomic and transcriptomic characteristics in eukaryotes

We present GTDrift, a comprehensive data resource that enables explorations of genomic and transcriptomic characteristics alongside proxies of the intensity of genetic drift in individual species. This resource encompasses data for 1506 eukaryotic species, including 1413 animals and 93 green plants, and is organized in three components. The first two components contain approximations of the effective population size, which serve as indicators of the extent of random genetic drift within each species. In the first component, we meticulously investigated public databases to assemble data on life history traits such as longevity, adult body length and body mass for a set of 979 species. The second component includes estimations of the ratio between the rate of non-synonymous substitutions and the rate of synonymous substitutions (dN/dS) in protein-coding sequences for 1324 species. This ratio provides an estimate of the efficiency of natural selection in purging deleterious substitutions. Additionally, we present polymorphism-derived N e estimates for 66 species. The third component encompasses various genomic and transcriptomic characteristics. With this component, we aim to facilitate comparative transcriptomics analyses across species, by providing easy-to-use processed data for more than 16 000 RNA-seq samples across 491 species. These data include intron-centered alternative splicing frequencies, gene expression levels and sequencing depth statistics for each species, obtained with a homogeneous analysis protocol. To enable cross-species comparisons, we provide orthology predictions for conserved single-copy genes based on BUSCO gene sets. To illustrate the possible uses of this database, we identify the most frequently used introns for each gene and we assess how the sequencing depth available for each species affects our power to identify major and minor splice variants.

Identification and Functional Characterization of Alternative Transcripts of LncRNA HNF1A-AS1 and Their Impacts on Cell Growth, Differentiation, Liver Diseases, and in Response to Drug Induction

The long non-coding RNA (lncRNA) hepatocyte nuclear factor-1 alpha (HNF1A) antisense RNA 1 (HNF1A-AS1) is an important lncRNA for liver growth, development, cell differentiation, and drug metabolism. Like many lncRNAs, HNF1A-AS1 has multiple annotated alternative transcripts in the human genome. Several fundamental biological questions are still not solved: (1) How many transcripts really exist in biological samples, such as liver samples and liver cell lines? (2) What are the expression patterns of different alternative HNF1A-AS1 transcripts at different conditions, including during cell growth and development, after exposure to xenobiotics (such as drugs), and in disease conditions, such as metabolic dysfunction-associated steatotic liver disease (MASLD), alcohol-associated liver disease (ALD) cirrhosis, and obesity? (3) Does the siRNA used in previous studies knock down one or multiple transcripts? (4) Do different transcripts have the same or different functions for gene regulation? The presented data confirm the existence of several annotated HNF1A-AS1 transcripts in liver samples and cell lines, but also identify some new transcripts, which are not annotated in the Ensembl genome database. Expression patterns of the identified HNF1A-AS1 transcripts are highly correlated with the cell differentiation of matured hepatocyte-like cells from human embryonic stem cells (hESC), growth and differentiation of HepaRG cells, in response to rifampicin induction, and in various liver disease conditions. The expression levels of the HNF1A-AS1 transcripts are also highly correlated to the expression of cytochrome P450 enzymes, such as CYP3A4, during HepaRG growth, differentiation, and in response to rifampicin induction.

Updating mRNA variants of the human RSK4 gene and their expression in different stressed situations

We determined RNA spectrum of the human RSK4 (hRSK4) gene (also called RPS6KA6) and identified 29 novel mRNA variants derived from alternative splicing, which, plus the NCBI-documented ones and the five we reported previously, totaled 50 hRSK4 RNAs that, by our bioinformatics analyses, encode 35 hRSK4 protein isoforms of 35-762 amino acids. Many of the mRNAs are bicistronic or tricistronic for hRSK4. The NCBI-normalized NM_014496.5 and the protein it encodes are designated herein as the Wt-1 mRNA and protein, respectively, whereas the NM_001330512.1 and the long protein it encodes are designated as the Wt-2 mRNA and protein, respectively. Many of the mRNA variants responded differently to different situations of stress, including serum starvation, a febrile temperature, treatment with ethanol or ethanol-extracted clove buds (an herbal medicine), whereas the same stressed situation often caused quite different alterations among different mRNA variants in different cell lines. Mosifloxacin, an antibiotics and also a functional inhibitor of hRSK4, could inhibit the expression of certain hRSK4 mRNA variants. The hRSK4 gene likely uses alternative splicing as a handy tool to adapt to different stressed situations, and the mRNA and protein multiplicities may partly explain the incongruous literature on its expression and comports.

Random genetic drift sets an upper limit on mRNA splicing accuracy in metazoans

Most eukaryotic genes undergo alternative splicing (AS), but the overall functional significance of this process remains a controversial issue. It has been noticed that the complexity of organisms (assayed by the number of distinct cell types) correlates positively with their genome-wide AS rate. This has been interpreted as evidence that AS plays an important role in adaptive evolution by increasing the functional repertoires of genomes. However, this observation also fits with a totally opposite interpretation: given that 'complex' organisms tend to have small effective population sizes (Ne), they are expected to be more affected by genetic drift, and hence more prone to accumulate deleterious mutations that decrease splicing accuracy. Thus, according to this 'drift barrier' theory, the elevated AS rate in complex organisms might simply result from a higher splicing error rate. To test this hypothesis, we analyzed 3496 transcriptome sequencing samples to quantify AS in 53 metazoan species spanning a wide range of Ne values. Our results show a negative correlation between Ne proxies and the genome-wide AS rates among species, consistent with the drift barrier hypothesis. This pattern is dominated by low abundance isoforms, which represent the vast majority of the splice variant repertoire. We show that these low abundance isoforms are depleted in functional AS events, and most likely correspond to errors. Conversely, the AS rate of abundant isoforms, which are relatively enriched in functional AS events, tends to be lower in more complex species. All these observations are consistent with the hypothesis that variation in AS rates across metazoans reflects the limits set by drift on the capacity of selection to prevent gene expression errors.

Multi-landmark alignment of genomic signals reveals conserved expression patterns across transcription start sites

The prevalent one-dimensional alignment of genomic signals to a reference landmark is a cornerstone of current methods to study transcription and its DNA-dependent processes but it is prone to mask potential relations among multiple DNA elements. We developed a systematic approach to align genomic signals to multiple locations simultaneously by expanding the dimensionality of the genomic-coordinate space. We analyzed transcription in human and uncovered a complex dependence on the relative position of neighboring transcription start sites (TSSs) that is consistently conserved among cell types. The dependence ranges from enhancement to suppression of transcription depending on the relative distances to the TSSs, their intragenic position, and the transcriptional activity of the gene. Our results reveal a conserved hierarchy of alternative TSS usage within a previously unrecognized level of genomic organization and provide a general methodology to analyze complex functional relationships among multiple types of DNA elements.

Genome-Wide Mapping of 5' Isoforms with 5'-Seq

The transcriptome is far more complex than previously assumed. Transcripts from the same gene can differ in terms of transcription start site, transcription end site, or pattern of splicing, and growing evidence supports the functional importance of these distinct transcript isoforms. Easily identifying these isoforms experimentally via library construction and high-throughput sequencing is crucial. Current library construction methods for identifying transcription start sites (5' transcript isoforms) involve large number of steps and (expensive) reagents, utilization of cDNA intermediates for adapter ligation, and are less suitable for studying low-abundance isoforms. Here, I describe a quick protocol for the generation of sequencing libraries to define capped 5' isoforms (5'-Seq) of various abundances in yeast and suggest a 5' isoform data analysis pipeline. The protocol relies on the utilization of a dephosphorylation-decapping method (oligo-capping) to generate a sequencing library from mRNA fragments and is a simplification of previously published 5' isoform protocols in terms of the handling steps, time, and cost. This method is exemplified using Saccharomyces cerevisiae mRNA, but it can be applied to various cellular conditions to study the effects of 5' transcript isoforms on transcriptional and/or translational regulation. © 2023 Wiley Periodicals LLC. Basic Protocol: Construction of a DNA sequencing library from capped 5' isoforms Support Protocol: Sequencing data analysis.

Analysis of regulatory sequences in exosomal DNA of NANOGP8

Exosomes participate in intercellular communication by transporting functionally active molecules. Such cargo from the original cells comprising proteins, micro-RNA, mRNA, single-stranded (ssDNA) and double-stranded DNA (dsDNA) molecules pleiotropically transforms the target cells. Although cancer cells secrete exosomes carrying a significant level of DNA capable of modulating oncogene expression in a recipient cell, the regulatory mechanism is unknown. We have previously reported that cancer cells produce exosomes containing NANOGP8 DNA. NANOGP8 is an oncogenic paralog of embryonic stem cell transcription factor NANOG and does not express in cells since it is a pseudogene. However, in this study, we evaluated NANOGP8 expression in glioblastoma multiforme (GBM) tissue from a surgically removed brain tumor of a patient. Significantly higher NANOGP8 transcription was observed in GBM cancer stem cells (CSCs) than in GBM cancer cells or neural stem cells (NSCs), despite identical sequences of NANOGP8-upstream genomic region in all the cell lines. This finding suggests that upstream genomic sequences of NANOGP8 may have environment-dependent promoter activity. We also found that the regulatory sequences upstream of exosomal NANOGP8 GBM DNA contain multiple core promoter elements, transcription factor binding sites, and segments of human viruses known for their oncogenic role. The exosomal sequence of NANOGP8-upstream GBM DNA is different from corresponding genomic sequences in CSCs, cancer cells, and NSCs as well as from the sequences reported by NCBI. These sequence dissimilarities suggest that exosomal NANOGP8 GBM DNA may not be a part of the genomic DNA. Exosomes possibly acquire this DNA from other sources where it is synthesized by an unknown mechanism. The significance of exosome-bestowed regulatory elements in the transcription of promoter-less retrogene such as NANOGP8 remains to be determined.

Gene product diversity: adaptive or not?

One gene does not equal one RNA or protein. The genomic revolution has revealed numerous different RNA and protein molecules that can be produced from one gene, such as circular RNAs generated by back-splicing, proteins with residues mismatching the genomic encoding because of RNA editing, and proteins extended in the C terminus via stop codon readthrough in translation. Are these diverse products results of exquisite gene regulations or imprecise biological processes? While there are cases where the gene product diversity appears beneficial, genome-scale patterns suggest that much of this diversity arises from nonadaptive, molecular errors. This finding has important implications for studying the functions of diverse gene products and for understanding the fundamental properties and evolution of cellular life.

The impact of Piscirickettsia salmonis infection on genome-wide DNA methylation profile in Atlantic Salmon

Salmon rickettsial septicaemia (SRS), caused by the bacteria Piscirickettsia salmonis (P. salmonis), is responsible for significant mortality in farmed Atlantic salmon in Chile. Currently there are no effective treatments or preventive measures for this disease, although genetic selection or genome engineering to increase salmon resistance to SRS are promising strategies. The accuracy and efficiency of these strategies are usually influenced by the available biological background knowledge of the disease. The aim of this study was to investigate DNA methylation changes in response to P. salmonis infection in the head kidney and liver tissue of Atlantic salmon, and the interaction between gene expression and DNA methylation in the same tissues. The head kidney and liver methylomes of 66 juvenile salmon were profiled using reduced representation bisulphite sequencing (RRBS), and compared between P. salmonis infected animals (3 and 9 days post infection) and uninfected controls, and between SRS resistant and susceptible fish. Methylation was correlated with matching RNA-Seq data from the same animals, revealing that methylation in the first exon leads to an important repression of gene expression. Head kidney methylation showed a clear response to the infection, associated with immunological processes such as actin cytoskeleton regulation, phagocytosis, endocytosis and pathogen associated pattern receptor signaling. Our results contribute to the growing understanding of the role of methylation in regulation of gene expression and response to infectious diseases and could inform the incorporation of epigenetic markers into genomic selection for disease resistant and the design of diagnostic epigenetic markers to better manage fish health in salmon aquaculture.

Building integrative functional maps of gene regulation

Every cell in the human body inherits a copy of the same genetic information. The three billion base pairs of DNA in the human genome, and the roughly 50 000 coding and non-coding genes they contain, must thus encode all the complexity of human development and cell and tissue type diversity. Differences in gene regulation, or the modulation of gene expression, enable individual cells to interpret the genome differently to carry out their specific functions. Here we discuss recent and ongoing efforts to build gene regulatory maps, which aim to characterize the regulatory roles of all sequences in a genome. Many researchers and consortia have identified such regulatory elements using functional assays and evolutionary analyses; we discuss the results, strengths and shortcomings of their approaches. We also discuss new techniques the field can leverage and emerging challenges it will face while striving to build gene regulatory maps of ever-increasing resolution and comprehensiveness.

Alternative Transcription Start Site Usage and Functional Implications in Pathogenic Fungi

Pathogenic fungi require delicate gene regulation mechanisms to adapt to diverse living environments and escape host immune systems. Recent advances in sequencing technology have exposed the complexity of the fungal genome, thus allowing the gradual disentanglement of multiple layers of gene expression control. Alternative transcription start site (aTSS) usage, previously reported to be prominent in mammals and to play important roles in physiopathology, is also present in fungi to fine-tune gene expression. Depending on the alteration in their sequences, RNA isoforms arising from aTSSs acquire different characteristics that significantly alter their stability and translational capacity as well as the properties and biologic functions of the resulting proteins. Disrupted control of aTSS usage has been reported to severely impair growth, virulence, and the infectious capacity of pathogenic fungi. Here, we discuss principle concepts, mechanisms, and the functional implication of aTSS usage in fungi.

Alternative splicing as a source of phenotypic diversity

A major goal of evolutionary genetics is to understand the genetic processes that give rise to phenotypic diversity in multicellular organisms. Alternative splicing generates multiple transcripts from a single gene, enriching the diversity of proteins and phenotypic traits. It is well established that alternative splicing contributes to key innovations over long evolutionary timescales, such as brain development in bilaterians. However, recent developments in long-read sequencing and the generation of high-quality genome assemblies for diverse organisms has facilitated comparisons of splicing profiles between closely related species, providing insights into how alternative splicing evolves over shorter timescales. Although most splicing variants are probably non-functional, alternative splicing is nonetheless emerging as a dynamic, evolutionarily labile process that can facilitate adaptation and contribute to species divergence.

PAMP-triggered genetic reprogramming involves widespread alternative transcription initiation and an immediate transcription factor wave

Immune responses triggered by pathogen-associated molecular patterns (PAMPs) are key to pathogen defense, but drivers and stabilizers of the growth-to-defense genetic reprogramming remain incompletely understood in plants. Here, we report a time-course study of the establishment of PAMP-triggered immunity (PTI) using cap analysis of gene expression. We show that around 15% of all transcription start sites (TSSs) rapidly induced during PTI define alternative transcription initiation events. From these, we identify clear examples of regulatory TSS change via alternative inclusion of target peptides or domains in encoded proteins, or of upstream open reading frames in mRNA leader sequences. We also find that 60% of PAMP response genes respond earlier than previously thought. In particular, a cluster of rapidly and transiently PAMP-induced genes is enriched in transcription factors (TFs) whose functions, previously associated with biological processes as diverse as abiotic stress adaptation and stem cell activity, appear to converge on growth restriction. Furthermore, examples of known potentiators of PTI, in one case under direct mitogen-activated protein kinase control, support the notion that the rapidly induced TFs could constitute direct links to PTI signaling pathways and drive gene expression changes underlying establishment of the immune state.

Genome-wide analysis of the CAD gene family reveals two bona fide CAD genes in oil palm

Cinnamyl alcohol dehydrogenase (CAD) is the key enzyme for lignin biosynthesis in plants. In this study, genome-wide analysis was performed to identify CAD genes in oil palm (Elaeis guineensis). Phylogenetic analysis was then conducted to select the bona fide EgCADs. The bona fide EgCAD genes and their respective 5' flanking regions were cloned and analysed. Their expression profiles were evaluated in various organs using RT-PCR. Seven EgCAD genes (EgCAD1-7) were identified and divided into four phylogenetic groups. EgCAD1 and EgCAD2 display high sequence similarities with other bona fide CADs and possess all the signature motifs of the bona fide CAD. They also display similar 3D protein structures. Gene expression analysis showed that EgCAD1 was expressed most abundantly in the root tissues, while EgCAD2 was expressed constitutively in all the tissues studied. EgCAD1 possesses only one transcription start site, while EgCAD2 has five. Interestingly, a TC microsatellite was found in the 5' flanking region of EgCAD2. The 5' flanking regions of EgCAD1 and EgCAD2 contain lignin-associated regulatory elements i.e. AC-elements, and other defence-related motifs, including W-box, GT-1 motif and CGTCA-motif. Altogether, these results imply that EgCAD1 and EgCAD2 are bona fide CAD involved in lignin biosynthesis during the normal development of oil palm and in response to stresses. Our findings shed some light on the roles of the bona fide CAD genes in oil palm and pave the way for manipulating lignin content in oil palm through a genetic approach.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-022-03208-0.

In-Depth Temporal Transcriptome Profiling of an Alphaherpesvirus Using Nanopore Sequencing

In this work, a long-read sequencing (LRS) technique based on the Oxford Nanopore Technology MinION platform was used for quantifying and kinetic characterization of the poly(A) fraction of bovine alphaherpesvirus type 1 (BoHV-1) lytic transcriptome across a 12-h infection period. Amplification-based LRS techniques frequently generate artefactual transcription reads and are biased towards the production of shorter amplicons. To avoid these undesired effects, we applied direct cDNA sequencing, an amplification-free technique. Here, we show that a single promoter can produce multiple transcription start sites whose distribution patterns differ among the viral genes but are similar in the same gene at different timepoints. Our investigations revealed that the circ gene is expressed with immediate–early (IE) kinetics by utilizing a special mechanism based on the use of the promoter of another IE gene (bicp4) for the transcriptional control. Furthermore, we detected an overlap between the initiation of DNA replication and the transcription from the bicp22 gene, which suggests an interaction between the two molecular machineries. This study developed a generally applicable LRS-based method for the time-course characterization of transcriptomes of any organism.

Unusual mammalian usage of TGA stop codons reveals that sequence conservation need not imply purifying selection

The assumption that conservation of sequence implies the action of purifying selection is central to diverse methodologies to infer functional importance. GC-biased gene conversion (gBGC), a meiotic mismatch repair bias strongly favouring GC over AT, can in principle mimic the action of selection, this being thought to be especially important in mammals. As mutation is GC→AT biased, to demonstrate that gBGC does indeed cause false signals requires evidence that an AT-rich residue is selectively optimal compared to its more GC-rich allele, while showing also that the GC-rich alternative is conserved. We propose that mammalian stop codon evolution provides a robust test case. Although in most taxa TAA is the optimal stop codon, TGA is both abundant and conserved in mammalian genomes. We show that this mammalian exceptionalism is well explained by gBGC mimicking purifying selection and that TAA is the selectively optimal codon. Supportive of gBGC, we observe (i) TGA usage trends are consistent at the focal stop codon and elsewhere (in UTR sequences); (ii) that higher TGA usage and higher TAA→TGA substitution rates are predicted by a high recombination rate; and (iii) across species the difference in TAA <-> TGA substitution rates between GC-rich and GC-poor genes is largest in genomes that possess higher between-gene GC variation. TAA optimality is supported both by enrichment in highly expressed genes and trends associated with effective population size. High TGA usage and high TAA→TGA rates in mammals are thus consistent with gBGC’s predicted ability to “drive” deleterious mutations and supports the hypothesis that sequence conservation need not be indicative of purifying selection. A general trend for GC-rich trinucleotides to reside at frequencies far above their mutational equilibrium in high recombining domains supports the generality of these results.

The BAM7 gene in Zea mays encodes a protein with similar structural and catalytic properties to Arabidopsis BAM2

Starch accumulates in the plastids of green plant tissues during the day to provide carbon for metabolism at night. Starch hydrolysis is catalyzed by members of the β-amylase (BAM) family, which in Arabidopsis thaliana (At) includes nine structurally and functionally diverse members. One of these enzymes, AtBAM2, is a plastid-localized enzyme that is unique among characterized β-amylases since it is tetrameric and exhibits sigmoidal kinetics. Sequence alignments show that the BAM domains of AtBAM7, a catalytically inactive, nuclear-localized transcription factor with an N-terminal DNA-binding domain, and AtBAM2 are more closely related to each other than they are to any other AtBAM. Since the BAM2 gene is found in more ancient lineages, it was hypothesized that the BAM7 gene evolved from BAM2. However, analysis of the genomes of 48 flowering plants revealed 12 species that appear to possess a BAM7 gene but lack a BAM2 gene. Upon closer inspection, these BAM7 proteins have a greater percent identity to AtBAM2 than to AtBAM7, and they share all of the AtBAM2 functional residues that BAM7 proteins normally lack. It is hypothesized that these genes may encode BAM2-like proteins although they are currently annotated as BAM7-like genes. To test this hypothesis, a cDNA for the short form of corn BAM7 (ZmBAM7-S) was designed for expression in Escherichia coli. Small-angle X-ray scattering data indicate that ZmBAM7-S has a tetrameric solution structure that is more similar to that of AtBAM2 than to that of AtBAM1. In addition, partially purified ZmBAM7-S is catalytically active and exhibits sigmoidal kinetics. Together, these data suggest that some BAM7 genes may encode a functional BAM2. Exploring and understanding the β-amylase gene structure could have an impact on the current annotation of genes.

MAGE genes encoding for embryonic development in cattle is mainly regulated by zinc finger transcription factor family and slightly by CpG Islands

Background

Melanoma Antigen Genes (MAGEs) are a family of genes that have piqued the interest of scientists for their unique expression pattern. The MAGE genes can be classified into type I MAGEs that expressed in testis and other reproductive tissues while type II MAGEs that have broad expression in many tissues. Several MAGE gene families are expressed in embryonic tissues in almost all eukaryotes, which is essential for embryo development mainly during germ cell differentiation. The aim of this study was to analyze the promoter regions and regulatory elements (transcription factors and CpG islands) of MAGE genes encoding for embryonic development in cattle.

Results

The in silico analysis revealed the highest promoter prediction scores (1.0) for TSS were obtained for two gene sequences (MAGE B4-like and MAGE-L2) while the lowest promoter prediction scores (0.8) was obtained for MAGE B17-like. It also revealed that the best common motif, motif IV, bear a resemblance with three TF families including Zinc-finger family, SMAD family and E2A related factors. From thirteen identified TFs candidates, majority of them (11/13) were clustered to Zinc-finger family serving as transcriptionally activator role whereas three (SP1, SP3 and Znf423) of them as activator or repressor in response to physiological and pathological stimuli. On the other hand we revealed slightly rich CpG islands in the gene body and promoter regions of MAGE genes encoding for embryonic development in cattle.

Conclusion

This in silico analysis of gene promoter regions and regulatory elements in MAGE genes could be useful for understanding regulatory networks and gene expression patterns during embryo development in bovine.

TSS-seq of Toxoplasma gondii sporozoites revealed a novel motif in stage-specific promoters

Toxoplasma gondii is one of the most common zoonotic protozoan parasites. It has three major infectious stages: rapidly multiplying tachyzoites (Tz), slowly replicating bradyzoites (Bz) and a resting/free-living stage, sporozoites (Sz). The regulatory mechanisms governing stage-specific gene expression are not fully understood. Few transcriptional start sites (TSS) are known for Sz. In this study, we obtained TSS of Sz using an oligo-capping method and RNA-seq analysis. We identified 1,043,503 TSS in the Sz transcriptome. These defined 38,973 TSS clusters, of which, 11,925 were expressed in Sz and 1535 TSS differentially expressed in Sz. Based on these data, we defined promoter regions and novel sporozoite stage-specific motifs using MEME. TGTANNTACA was distributed around -55 to -75 regions from each TSS. Interestingly, the same motif was reported in another apicomplexan, Plasmodium berghei, as a cis-element of female-specific gametocyte genes, implying the presence of common regulatory machinery. Further comparative analysis should better define the distribution and function of these elements in other members of this important parasitic phylum.

Mammalian splicing divergence is shaped by drift, buffering in trans, and a scaling law

This study globally investigates the allelic splicing pattern in multiple tissues of an F1 hybrid mouse and reveals the underlying driving forces shaping such tissue-dependent splicing divergence.

Alternative splicing is ubiquitous, but the mechanisms underlying its pattern of evolutionary divergence across mammalian tissues are still underexplored. Here, we investigated the cis-regulatory divergences and their relationship with tissue-dependent trans-regulation in multiple tissues of an F1 hybrid between two mouse species. Large splicing changes between tissues are highly conserved and likely reflect functional tissue-dependent regulation. In particular, micro-exons frequently exhibit this pattern with high inclusion levels in the brain. Cis-divergence of splicing appears to be largely non-adaptive. Although divergence is in general associated with higher densities of sequence variants in regulatory regions, events with high usage of the dominant isoform apparently tolerate more mutations, explaining why their exon sequences are highly conserved but their intronic splicing site flanking regions are not. Moreover, we demonstrate that non-adaptive mutations are often masked in tissues where accurate splicing likely is more important, and experimentally attribute such buffering effect to trans-regulatory splicing efficiency.

Identifying transcript 5' capped ends in Plasmodium falciparum

Background

The genome of the human malaria parasite Plasmodium falciparum is poorly annotated, in particular, the 5' capped ends of its mRNA transcripts. New approaches are needed to fully catalog P. falciparum transcripts for understanding gene function and regulation in this organism.

Methods

We developed a transcriptomic method based on next-generation sequencing of complementary DNA (cDNA) enriched for full-length fragments using eIF4E, a 5' cap-binding protein, and an unenriched control. DNA sequencing adapter was added after enrichment of full-length cDNA using two different ligation protocols. From the mapped sequence reads, enrichment scores were calculated for all transcribed nucleotides and used to calculate P-values of 5' capped nucleotide enrichment. Sensitivity and accuracy were increased by combining P-values from replicate experiments. Data were obtained for P. falciparum ring, trophozoite and schizont stages of intra-erythrocytic development.

Results

5' capped nucleotide signals were mapped to 17,961 non-overlapping P. falciparum genomic intervals. Analysis of the dominant 5' capped nucleotide in these genomic intervals revealed the presence of two groups with distinctive epigenetic features and sequence patterns. A total of 4,512 transcripts were annotated as 5' capped based on the correspondence of 5' end with 5' capped nucleotide annotated from full-length cDNA data.

Discussion

The presence of two groups of 5' capped nucleotides suggests that alternative mechanisms may exist for producing 5' capped transcript ends in P. falciparum. The 5' capped transcripts that are antisense, outside of, or partially overlapping coding regions may be important regulators of gene function in P. falciparum.

Mammalian circular RNAs result largely from splicing errors

Ubiquitous in eukaryotes, circular RNAs (circRNAs) comprise a large class of mostly non-coding RNAs produced by back-splicing. Although some circRNAs have demonstrated biochemical activities, whether most circRNAs are functional is unknown. Here, we test the hypothesis that circRNA production primarily results from splicing error and so is deleterious instead of beneficial. In support of the error hypothesis, our analysis of RNA sequencing data from 11 shared tissues of humans, macaques, and mice finds that (1) back-splicing is much rarer than linear-splicing, (2) the rate of back-splicing diminishes with the splicing amount, (3) the overall prevalence of back-splicing in a species declines with its effective population size, and (4) circRNAs are overall evolutionarily unconserved. We estimate that more than 97% of the observed circRNA production is deleterious. We identify a small number of functional circRNA candidates, and the genome-wide trend strongly suggests that circRNAs are largely non-functional products of splicing errors.

Evolution of tissue and developmental specificity of transcription start sites in Bos taurus indicus

To further the understanding of the evolution of transcriptional regulation, we profiled genome-wide transcriptional start sites (TSSs) in two sub-species, Bos taurus taurus and Bos taurus indicus, that diverged approximately 500,000 years ago. Evolutionary and developmental-stage differences in TSSs were detected across the sub-species, including translocation of dominant TSS and changes in TSS distribution. The 16% of all SNPs located in significant differentially used TSS clusters across sub-species had significant shifts in allele frequency (472 SNPs), indicating they may have been subject to selection. In spleen and muscle, a higher relative TSS expression was observed in Bos indicus than Bos taurus for all heat shock protein genes, which may be responsible for the tropical adaptation of Bos indicus.

Forutan et al. measure RNA expression and map transcription start sites (TSSs) between two sub-species of cattle in adult and fetal tissues and demonstrate translocation and changes in TSS distribution between the sub-species. This study provides insight into cattle adaptation and provides support of rapidly evolving TSS between two recently diverged subspecies.

Assessing the functional relevance of splice isoforms

Alternative splicing of messenger RNA can generate an array of mature transcripts, but it is not clear how many go on to produce functionally relevant protein isoforms. There is only limited evidence for alternative proteins in proteomics analyses and data from population genetic variation studies indicate that most alternative exons are evolving neutrally. Determining which transcripts produce biologically important isoforms is key to understanding isoform function and to interpreting the real impact of somatic mutations and germline variations. Here we have developed a method, TRIFID, to classify the functional importance of splice isoforms. TRIFID was trained on isoforms detected in large-scale proteomics analyses and distinguishes these biologically important splice isoforms with high confidence. Isoforms predicted as functionally important by the algorithm had measurable cross species conservation and significantly fewer broken functional domains. Additionally, exons that code for these functionally important protein isoforms are under purifying selection, while exons from low scoring transcripts largely appear to be evolving neutrally. TRIFID has been developed for the human genome, but it could in principle be applied to other well-annotated species. We believe that this method will generate valuable insights into the cellular importance of alternative splicing.

Readthrough Errors Purge Deleterious Cryptic Sequences, Facilitating the Birth of Coding Sequences

De novo protein-coding innovations sometimes emerge from ancestrally noncoding DNA, despite the expectation that translating random sequences is overwhelmingly likely to be deleterious. The "preadapting selection" hypothesis claims that emergence is facilitated by prior, low-level translation of noncoding sequences via molecular errors. It predicts that selection on polypeptides translated only in error is strong enough to matter and is strongest when erroneous expression is high. To test this hypothesis, we examined noncoding sequences located downstream of stop codons (i.e., those potentially translated by readthrough errors) in Saccharomyces cerevisiae genes. We identified a class of "fragile" proteins under strong selection to reduce readthrough, which are unlikely substrates for co-option. Among the remainder, sequences showing evidence of readthrough translation, as assessed by ribosome profiling, encoded C-terminal extensions with higher intrinsic structural disorder, supporting the preadapting selection hypothesis. The cryptic sequences beyond the stop codon, rather than spillover effects from the regular C-termini, are primarily responsible for the higher disorder. Results are robust to controlling for the fact that stronger selection also reduces the length of C-terminal extensions. These findings indicate that selection acts on 3' UTRs in Saccharomyces cerevisiae to purge potentially deleterious variants of cryptic polypeptides, acting more strongly in genes that experience more readthrough errors.

Mammalian Alternative Translation Initiation Is Mostly Nonadaptive

Alternative translation initiation (ATLI) refers to the existence of multiple translation initiation sites per gene and is a widespread phenomenon in eukaryotes. ATLI is commonly assumed to be advantageous through creating proteome diversity or regulating protein synthesis. We here propose an alternative hypothesis that ATLI arises primarily from nonadaptive initiation errors presumably due to the limited ability of ribosomes to distinguish sequence motifs truly signaling translation initiation from similar sequences. Our hypothesis, but not the adaptive hypothesis, predicts a series of global patterns of ATLI, all of which are confirmed at the genomic scale by quantitative translation initiation sequencing in multiple human and mouse cell lines and tissues. Similarly, although many codons differing from AUG by one nucleotide can serve as start codons, our analysis suggests that using non-AUG start codons is mostly disadvantageous. These and other findings strongly suggest that ATLI predominantly results from molecular error, requiring a major revision of our understanding of the precision and regulation of translation initiation.

Transcription initiation of distant core promoters in a large-sized genome of an insect

Background

Core promoters have a substantial influence on various steps of transcription, including initiation, elongation, termination, polyadenylation, and finally, translation. The characterization of core promoters is crucial for exploring the regulatory code of transcription initiation. However, the current understanding of insect core promoters is focused on those of Diptera (especially Drosophila) species with small genome sizes.

Results

Here, we present an analysis of the transcription start sites (TSSs) in the migratory locust, Locusta migratoria, which has a genome size of 6.5 Gb. The genomic differences, including lower precision of transcription initiation and fewer constraints on the distance from transcription factor binding sites or regulatory elements to TSSs, were revealed in locusts compared with Drosophila insects. Furthermore, we found a distinct bimodal log distribution of the distances from the start codons to the core promoters of locust genes. We found stricter constraints on the exon length of mRNA leaders and widespread expression activity of the distant core promoters in locusts compared with fruit flies. We further compared core promoters in seven arthropod species across a broad range of genome sizes to reinforce our results on the emergence of distant core promoters in large-sized genomes.

Conclusions

In summary, our results provide novel insights into the effects of genome size expansion on distant transcription initiation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01004-5.

Characterization of Bicistronic Transcription in Budding Yeast

Bicistronic transcripts (operon-like transcripts) have occasionally been reported in eukaryotes, including unicellular yeasts, plants, and humans, despite the fact that they lack trans-splice mechanisms. However, the characteristics of eukaryotic bicistronic transcripts are poorly understood, except for those in nematodes. Here, we describe the genomic, transcriptomic, and ribosome profiling features of bicistronic transcripts in unicellular yeasts. By comparing the expression level of bicistronic transcripts with their monocistronic equivalents, we identify two main categories of bicistronic transcripts: highly and lowly expressed. These two categories exhibit quite different features. First, highly expressed bicistronic transcripts have higher conservation within and between strains and shorter intergenic spacers with higher GC content and less stable secondary structure. Second, genes in highly expressed bicistronic transcripts have lower translation efficiency, with the second gene showing statistically significant lower translation efficiency than the first. Finally, the genes found in these highly expressed bicistronic transcripts tend to be younger, with more recent origins. Together, these results suggest that bicistronic transcripts in yeast are heterogeneous. We further propose that at least some highly expressed bicistronic transcripts appear to play a role in modulating monocistronic translation.IMPORTANCE Operons, where a single mRNA transcript encodes multiple adjacent proteins, are a widespread feature of bacteria and archaea. In contrast, the genes of eukaryotes are generally considered monocistronic. However, a number of studies have revealed the presence of bicistronic transcripts in eukaryotes, including humans. The basic features of these transcripts are largely unknown in eukaryotes, especially in organisms lacking trans-splice mechanisms. Our analyses characterize bicistronic transcripts in one such eukaryotic group, yeasts. We show that highly expressed bicistronic transcripts have unusual features compared to lowly expressed bicistronic transcripts, with several features influencing translational modulation.

In Vivo Validation of Alternative FDXR Transcripts in Human Blood in Response to Ionizing Radiation

Following cell stress such as ionising radiation (IR) exposure, multiple cellular pathways are activated. We recently demonstrated that ferredoxin reductase (FDXR) has a remarkable IR-induced transcriptional responsiveness in blood. Here, we provided a first comprehensive FDXR variant profile following DNA damage. First, specific quantitative real-time polymerase chain reaction (qPCR) primers were designed to establish dose-responses for eight curated FDXR variants, all up-regulated after IR in a dose-dependent manner. The potential role of gender on the expression of these variants was tested, and neither the variants response to IR nor the background level of expression was profoundly affected; moreover, in vitro induction of inflammation temporarily counteracted IR response early after exposure. Importantly, transcriptional up-regulation of these variants was further confirmed in vivo in blood of radiotherapy patients. Full-length nanopore sequencing was performed to identify other FDXR variants and revealed the high responsiveness of FDXR-201 and FDXR-208. Moreover, FDXR-218 and FDXR-219 showed no detectable endogenous expression, but a clear detection after IR. Overall, we characterised 14 FDXR transcript variants and identified for the first time their response to DNA damage in vivo. Future studies are required to unravel the function of these splicing variants, but they already represent a new class of radiation exposure biomarkers.

High Transcriptional Error Rates Vary as a Function of Gene Expression Level

Errors in gene transcription can be costly, and organisms have evolved to prevent their occurrence or mitigate their costs. The simplest interpretation of the drift barrier hypothesis suggests that species with larger population sizes would have lower transcriptional error rates. However, Escherichia coli seems to have a higher transcriptional error rate than species with lower effective population sizes, for example Saccharomyces cerevisiae. This could be explained if selection in E. coli were strong enough to maintain adaptations that mitigate the consequences of transcriptional errors through robustness, on a gene by gene basis, obviating the need for low transcriptional error rates and associated costs of global proofreading. Here, we note that if selection is powerful enough to evolve local robustness, selection should also be powerful enough to locally reduce error rates. We therefore predict that transcriptional error rates will be lower in highly abundant proteins on which selection is strongest. However, we only expect this result when error rates are high enough to significantly impact fitness. As expected, we find such a relationship between expression and transcriptional error rate for non-C→U errors in E. coli (especially G→A), but not in S. cerevisiae. We do not find this pattern for C→U changes in E. coli, presumably because most deamination events occurred during sample preparation, but do for C→U changes in S. cerevisiae, supporting the interpretation that C→U error rates estimated with an improved protocol, and which occur at rates comparable with E. coli non-C→U errors, are biological.

Universal promoter scanning by Pol II during transcription initiation in Saccharomyces cerevisiae

BACKGROUND

The majority of eukaryotic promoters utilize multiple transcription start sites (TSSs). How multiple TSSs are specified at individual promoters across eukaryotes is not understood for most species. In Saccharomyces cerevisiae, a pre-initiation complex (PIC) comprised of Pol II and conserved general transcription factors (GTFs) assembles and opens DNA upstream of TSSs. Evidence from model promoters indicates that the PIC scans from upstream to downstream to identify TSSs. Prior results suggest that TSS distributions at promoters where scanning occurs shift in a polar fashion upon alteration in Pol II catalytic activity or GTF function.

RESULTS

To determine the extent of promoter scanning across promoter classes in S. cerevisiae, we perturb Pol II catalytic activity and GTF function and analyze their effects on TSS usage genome-wide. We find that alterations to Pol II, TFIIB, or TFIIF function widely alter the initiation landscape consistent with promoter scanning operating at all yeast promoters, regardless of promoter class. Promoter architecture, however, can determine the extent of promoter sensitivity to altered Pol II activity in ways that are predicted by a scanning model.

CONCLUSIONS

Our observations coupled with previous data validate key predictions of the scanning model for Pol II initiation in yeast, which we term the shooting gallery. In this model, Pol II catalytic activity and the rate and processivity of Pol II scanning together with promoter sequence determine the distribution of TSSs and their usage.

The preponderance of nonsynonymous A-to-I RNA editing in coleoids is nonadaptive

A-to-I editing enzymatically converts the base adenosine (A) in RNA molecules to inosine (I), which is recognized as guanine (G) in translation. Exceptionally abundant A-to-I editing was recently discovered in the neural tissues of coleoids (octopuses, squids, and cuttlefishes), with a greater fraction of nonsynonymous sites than synonymous sites subject to high levels of editing. Although this phenomenon is thought to indicate widespread adaptive editing, its potential advantage is unknown. Here we propose an alternative, nonadaptive explanation. Specifically, increasing the cellular editing activity permits some otherwise harmful G-to-A nonsynonymous substitutions, because the As are edited to Is at sufficiently high levels. These high editing levels are constrained upon substitutions, resulting in the predominance of nonsynonymous editing at highly edited sites. Our evidence for this explanation suggests that the prevalent nonsynonymous editing in coleoids is generally nonadaptive, as in species with much lower editing activities.

The neural tissues of coleoids have a greater fraction of nonsynonymous sites than synonymous sites subject to high levels of A-to-I RNA editing, a pattern thought to indicate widespread adaptive editing. Here the authors propose and provide evidence for an alternative, nonadaptive explanation.

In eubacteria, unlike eukaryotes, there is no evidence for selection favouring fail-safe 3' additional stop codons

Errors throughout gene expression are likely deleterious, hence genomes are under selection to ameliorate their consequences. Additional stop codons (ASCs) are in-frame nonsense ‘codons’ downstream of the primary stop which may be read by translational machinery should the primary stop have been accidentally read through. Prior evidence in several eukaryotes suggests that ASCs are selected to prevent potentially-deleterious consequences of read-through. We extend this evidence showing that enrichment of ASCs is common but not universal for single cell eukaryotes. By contrast, there is limited evidence as to whether the same is true in other taxa. Here, we provide the first systematic test of the hypothesis that ASCs act as a fail-safe mechanism in eubacteria, a group with high read-through rates. Contra to the predictions of the hypothesis we find: there is paucity, not enrichment, of ASCs downstream; substitutions that degrade stops are more frequent in-frame than out-of-frame in 3’ sequence; highly expressed genes are no more likely to have ASCs than lowly expressed genes; usage of the leakiest primary stop (TGA) in highly expressed genes does not predict ASC enrichment even compared to usage of non-leaky stops (TAA) in lowly expressed genes, beyond downstream codon +1. Any effect at the codon immediately proximal to the primary stop can be accounted for by a preference for a T/U residue immediately following the stop, although if anything, TT- and TC- starting codons are preferred. We conclude that there is no compelling evidence for ASC selection in eubacteria. This presents an unusual case in which the same error could be solved by the same mechanism in eukaryotes and prokaryotes but is not. We discuss two possible explanations: that, owing to the absence of nonsense mediated decay, bacteria may solve read-through via gene truncation and in eukaryotes certain prion states cause raised read-through rates.

In all organisms, gene expression is error-prone. One such error, translational read-through, occurs where the primary stop codon of an expressed gene is missed by the translational machinery. Failure to terminate is likely to be costly, hence genomes are under selection to prevent this from happening. One proposed error-proofing strategy involves in-frame proximal additional stop codons (ASCs) which may act as a ‘fail-safe’ mechanism by providing another opportunity for translation to terminate. There is evidence for ASC enrichment in several eukaryotes. We extend this evidence showing it to be common but not universal in single celled eukaryotes. However, the situation in bacteria is poorly understood, despite bacteria having high read-through rates. Here, we test the fail-safe hypothesis within a broad range of bacteria. To our surprise, we find that not only are ASCs not enriched, but they may even be selected against. This provides evidence for an unusual circumstance where eukaryotes and prokaryotes could solve the same problem the same way but don’t. What are we to make of this? We suggest that if read-through is the problem, ASCs are not necessarily the expected solution. Owing to the absence of nonsense-mediated decay, a process that makes gene truncation in eukaryotes less viable, we propose bacteria may rescue a leaky stop by mutation that creates a new stop upstream. Alternatively, raised read-through rates in some particular conditions in eukaryotes might explain the difference.

Stop-codon read-through arises largely from molecular errors and is generally nonadaptive

Stop-codon read-through refers to the phenomenon that a ribosome goes past the stop codon and continues translating into the otherwise untranslated region (UTR) of a transcript. Recent ribosome-profiling experiments in eukaryotes uncovered widespread stop-codon read-through that also varies among tissues, prompting the adaptive hypothesis that stop-codon read-through is an important, regulated mechanism for generating proteome diversity. Here we propose and test a competing hypothesis that stop-codon read-through arises mostly from molecular errors and is largely nonadaptive. The error hypothesis makes distinct predictions about the probability of read-through, frequency of sequence motifs for read-through, and conservation of the read-through region, each of which is supported by genome-scale data from yeasts and fruit flies. Thus, except for the few cases with demonstrated functions, stop-codon read-through is generally nonadaptive. This finding, along with other molecular errors recently quantified, reveals a much less precise or orderly cellular life than is commonly thought.