PubMed HealthSearch

SEARCH · PubMed Health

Results for “Ribosome profiling”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 recordsLinked to original sources

RP-REP Ribosomal Profiling Reports: an open-source cloud-enabled framework for reproducible ribosomal profiling data processing, analysis, and result reporting.

Ribosomal profiling is an emerging experimental technology to measure protein synthesis by sequencing short mRNA fragments undergoing translation in ribosomes. Applied on the genome wide scale, this is a powerful tool to profile global protein synthesis within cell populations of interest. Such information can be utilized for biomarker discovery and detection of treatment-responsive genes. However, analysis of ribosomal profiling data requires careful preprocessing to reduce the impact of artifacts and dedicated statistical methods for visualizing and modeling the high-dimensional discrete read count data. Here we present Ribosomal Profiling Reports (RP-REP), a new open-source cloud-enabled software that allows users to execute start-to-end gene-level ribosomal profiling and RNA-Seq analysis on a pre-configured Amazon Virtual Machine Image (AMI) hosted on AWS or on the user's own Ubuntu Linux server. The software works with FASTQ files stored locally, on AWS S3, or at the Sequence Read Archive (SRA). RP-REP automatically executes a series of customizable steps including filtering of contaminant RNA, enrichment of true ribosomal footprints, reference alignment and gene translation quantification, gene body coverage, CRAM compression, reference alignment QC, data normalization, multivariate data visualization, identification of differentially translated genes, and generation of heatmaps, co-translated gene clusters, enriched pathways, and other custom visualizations. RP-REP provides functionality to contrast RNA-SEQ and ribosomal profiling results, and calculates translational efficiency per gene. The software outputs a PDF report and publication-ready table and figure files. As a use case, we provide RP-REP results for a dengue virus study that tested cytosol and endoplasmic reticulum cellular fractions of human Huh7 cells pre-infection and at 6 h, 12 h, 24 h, and 40 h post-infection. Case study results, Ubuntu installation scripts, and the most recent RP-REP source code are accessible at GitHub. The cloud-ready AMI is available at AWS (AMI ID: RPREP RSEQREP (Ribosome Profiling and RNA-Seq Reports) v2.1 (ami-00b92f52d763145d3)).

AMI

A custom library construction method for super-resolution ribosome profiling in Arabidopsis.

BACKGROUND: Ribosome profiling, also known as Ribo-seq, is a powerful technique to study genome-wide mRNA translation. It reveals the precise positions and quantification of ribosomes on mRNAs through deep sequencing of ribosome footprints. We previously optimized the resolution of this technique in plants. However, several key reagents in our original method have been discontinued, and thus, there is an urgent need to establish an alternative protocol. RESULTS: Here we describe a step-by-step protocol that combines our optimized ribosome footprinting in plants with available custom library construction methods established in yeast and bacteria. We tested this protocol in 7-day-old Arabidopsis seedlings and evaluated the quality of the sequencing data regarding ribosome footprint length, mapped genomic features, and the periodic properties corresponding to actively translating ribosomes through open resource bioinformatic tools. We successfully generated high-quality Ribo-seq data comparable with our original method. CONCLUSIONS: We established a custom library construction method for super-resolution Ribo-seq in Arabidopsis. The experimental protocol and bioinformatic pipeline should be readily applicable to other plant tissues and species.

3-nt periodicity

Age-related changes in ribosomal profiles and in vitro protein synthesis in skeletal muscle during fasting and subsequent refeeding of rats.

Age-dependent changes in ribosomal profiles and in in vitro protein synthesis in skeletal muscle (gastrocnemius muscle) were studied in rats subjected to fasting and subsequent refeeding. Male rats of 4, 7, and 21 weeks of age were used. Incorporation of 14C-leucine into protein in vitro using skeletal muscle ribosomes from fed rats decreased with the age of the animals. Sucrose density gradient analyses of ribosomal profiles revealed that the changes in protein synthesis were accompanied by progressive declines in the concentrations of polyribosomes and total RNA in skeletal muscle. After fasting for 1 or 2 days, protein synthesis in muscle of young rats dropped markedly and this change could be attributed to the fall in the concentrations of polyribosomes. Also, the concentration and total RNA in the muscles of the young rats decreased. During refeeding of rats after fasting for 2 days, the incorporation of 14C-amino acids into protein recovered more efficiently using ribosomes of skeletal muscle of young animals compared with using those of older animals. This rapid increase in the protein synthesizing activity of skeletal muscle may be related to the efficient recovery of the ribosomes toward heavier aggregates. The results indicate that young animals adapt more efficiently to fasting and refeeding than do older animals.

Aging

Localization of host poly(A)+ mRNA in the ribosome profile of mengovirus-infected Ehrlich ascites tumor cells.

The distribution of poly(A)+ mRNA among polysomes, monosomes, and ribosome-free supernatant fractions after mengovirus infection of Ehrlich ascites tumor (EAT) cells was investigated employing sucrose gradient centrifugation of their corresponding postnuclear supernatants. Poly(A)+ mRNA was isolated from sucrose gradient fractions and quantitated in a cell-free protein synthesizing system from uninfected EAT cells. It was also localized by annealing [3H]-poly(U) to the poly(A)-tracts of mRNA present in the sucrose gradient fractions. Both experiments revealed a gradual shift of host poly(A)+ mRNA from large to small polysomes and monosomes, respectively, with the time postinfection. The greatest part of host template RNA appears to remain ribosome-bound and only a fraction seems to be detached from the ribosomes in the course of mengovirus infection. At the end of the infectious cycle, 8 h postinfection, approximately 70% of the poly(A)+ mRNA detected in uninfected cells is still biologically active, but not translated in vivo, in agreement with data from the [3H]poly(U)hybridization experiment.

Animals

seq2ribo: structure-aware integration of machine learning and simulation to predict ribosome location profiles from RNA sequences.

MOTIVATION: Ribosome dynamics are vital in the process of protein expression. Current methods rely on ribosome profiling (Ribo-seq), RNA-seq profiles, and full genomic context. This restricts their use in de novo sequence design, like messenger RNA (mRNA) vaccines. Simulation-only approaches like the Totally Asymmetric Simple Exclusion Process (TASEP) oversimplify translation by focusing solely on codon elongation times. RESULTS: We present seq2ribo, a hybrid simulation and machine learning framework that predicts ribosome A-site locations using only an mRNA sequence as input. Our method first employs a novel structure-aware TASEP (sTASEP), which models translation using a comprehensive set of fitted parameters that include codon wait times and structural features, such as local angles, base-pairing, and discrete positional buckets. The ribosome locations generated by sTASEP are then processed by a polisher model, which learns to refine the simulated ribosome distributions. seq2ribo provides high-fidelity predictions of ribosome locations across diverse cell types (iPSC, HEK293, LCL, and RPE-1), significantly outperforming baselines. seq2ribo is the first method to achieve meaningful positional correlation with observed ribosome profiles from sequence alone, reaching transcript-level Pearson correlations up to 0.920 and within-transcript shape correlations up to 0.186, where all baselines yield near-zero values on these metrics. seq2ribo also reduces elementwise error by up to 37.7% relative to the sequence-only Translatomer baseline. By adding a task-specific head, seq2ribo achieves Pearson correlations up to 0.732 with experimental translation efficiency (TE) across several cell lines, and up to 0.903 with measured protein expression. By operating from sequence alone, seq2ribo provides a new tool for synthetic biology, enabling the rational design and optimization of mRNA sequences without the need for expression-level data or genomic context. AVAILABILITY: seq2ribo is available at https://github.com/Kingsford-Group/seq2ribo.

Machine Learning

seq2ribo: Structure-aware integration of machine learning and simulation to predict ribosome location profiles from RNA sequences.

MOTIVATION: Ribosome dynamics are vital in the process of protein expression. Current methods rely on ribosome profiling (Ribo-seq), RNA-seq profiles, and full genomic context. This restricts their use in de novo sequence design, like messenger RNA (mRNA) vaccines. Simulation-only approaches like the Totally Asymmetric Simple Exclusion Process (TASEP) oversimplify translation by focusing solely on codon elongation times. RESULTS: We present seq2ribo, a hybrid simulation and machine learning framework that predicts ribosome A-site locations using only an mRNA sequence as input. Our method first employs a novel structure-aware TASEP (sTASEP), which models translation using a comprehensive set of fitted parameters that include codon wait times and structural features, such as local angles, base-pairing, and discrete positional buckets. The ribosome locations generated by sTASEP are then processed by a polisher model, which learns to refine the simulated ribosome distributions. seq2ribo provides high-fidelity predictions of ribosome locations across diverse cell types (iPSC, HEK293, LCL, and RPE-1), significantly outperforming baselines. seq2ribo is the first method to achieve meaningful positional correlation with observed ribosome profiles from sequence alone, reaching transcript-level Pearson correlations up to 0.920 and within-transcript shape correlations up to 0.186, where all baselines yield near-zero values on these metrics. seq2ribo also reduces elementwise error by up to 37.7% relative to the sequence-only Translatomer baseline. By adding a task-specific head, seq2ribo achieves Pearson correlations up to 0.732 with experimental translation efficiency (TE) across several cell lines, and up to 0.903 with measured protein expression. By operating from sequence alone, seq2ribo provides a new tool for synthetic biology, enabling the rational design and optimization of mRNA sequences without the need for expression-level data or genomic context.

Journal Article

A novel regulation on the developmental checkpoint protein Sda that controls sporulation and biofilm formation in Bacillus subtilis.

UNLABELLED: Biofilm formation by Bacillus subtilis is triggered by an unusually simple environmental sensing mechanism. Certain serine codons, the four TCN codons (N for A, T, C, or G), in the gene for the biofilm repressor SinR caused lowered SinR translation and subsequent biofilm induction during transition from exponential to stationary growth. Global ribosome profiling showed that ribosomes pause when translating the four UCN (U for T on the mRNA) serine codons on mRNA, but not the two AGC/AGU serine codons. We proposed a serine codon hierarchy (AGC/AGT vs TCN) in that genes enriched in the TCN serine codons may experience reduced translation efficiency when serine is limited. In this study, we designed an algorithm to score all protein-coding genes in B. subtilis NCIB3610 based on the serine codon hierarchy. We generated a short list of 50 genes that could be subject to regulation by this novel mechanism. We further investigated one such gene from the list, sda, which encodes a developmental checkpoint protein regulating both sporulation and biofilm formation. We showed that synonymously switching the TCN serine codons to AGC in sda led to delayed biofilm formation and sporulation. This engineered strain also outgrew strains with other synonymously substituted sda alleles (TCN) in competition assays for biofilm formation and sporulation. Finally, we showed that the AGC serine codon substitutions in sda elevated the Sda protein levels. This serine codon hierarchy-based novel signaling mechanism could be exploited by bacteria in adapting to stationary phase and regulating important biological processes. IMPORTANCE: Genome-wide ribosome profiling in Bacillus subtilis shows that under serine limitation, ribosomes pause on the four TCN (N for A, C, G, and T), but not AGC/AGT serine codons, during translation at a global scale. This serine codon hierarchy (AGC/T vs TCN) differentially influences the translation efficiency of genes enriched in certain serine codons. In this study, we designed an algorithm to score all 4,000+ genes in the B. subtilis genome and generated a list of 50 genes that could be subject to this novel serine codon hierarchy-mediated regulation. We further investigated one such gene, sda, encoding a developmental checkpoint protein. We show that sda and cell developments controlled by Sda are also regulated by this novel mechanism.

Bacillus subtilis

Ribosome dynamics at the conserved PGP motif governs 2A peptide-bond-skipping efficiency.

Viral 2A oligopeptides drive an unusual ribosome recoding event in which peptide-bond formation fails at a conserved PG↓P motif, producing two discrete proteins without canonical termination. Despite decades of study, the molecular basis of 2A-mediated peptide-bond skipping remains poorly understood. Here, we combine quantitative 2A reporters with high-resolution ribosome profiling to interrogate ribosome dynamics at the core 2A sequences. We identify a pausing event at the terminal proline codon of the PGP motif that functions as a kinetic decision point: ribosome dwell time at this site inversely correlates with skipping efficiency. Increasing nascent chain flexibility by inserting linkers immediately upstream of the 2A sequence reduces ribosome occupancy at the terminal proline codon and enhances peptide-bond skipping. Strikingly, amino acid repeats positioned distally upstream also modulate 2A activity, indicating long-range coupling between nascent chain properties outside of the ribosome and the peptidyl transferase center inside the ribosome. In particular, hydrophobic residues potently suppress skipping, an effect that can be rescued by extending flexible segments within the peptide exit tunnel. Together, our findings support a model in which nascent chain features-beyond the core 2A motif-dynamically tune ribosomal recoding efficiency through co-translational feedback into the catalytic center.

Ribosomes

Ribo-ITP enables identification of translons from limited input samples.

In the last decade, an unexpectedly large number of translated regions (translons) have been discovered using ribosome profiling and proteomics. Translons can act as regulatory elements or encode functional micropeptides. However, identification of translons has been limited to cell lines or large organs due to high input requirements for conventional ribosome profiling and mass spectrometry. Here, we address this input limitation using Ribo-ITP on difficult-to-collect samples such as microdissected hippocampal tissues and single preimplantation embryos to identify thousands of translons. To test the translational capacity of the identified translons, we engineer a translon-dependent GFP reporter system and detect expression of translons initiating at ATG and near-cognate start codons in mouse embryonic stem cells (mESCs). We identify distinct expression patterns of translons using a comparative analysis of more than a thousand ribosome profiling datasets across a wide range of cell types. Further, using a machine learning model, we predict that specific upstream translons in synaptically enriched mRNAs regulate translation efficiency of the annotated coding region. Taken together, we present a proof-of-concept study to identify non-canonical translation events from low input samples which can be applied to cell and tissue types inaccessible to conventional methods.

Animals

Ribosomal distribution in a polyamine auxotroph of Escherichia coli.

The distribution of ribosomal particles has been studied in a polyamine-deficient mutant of Escherichia coli by sucrose gradient centrifugation analysis. Lysates from starved cells contained less 70S monomers and 30S subunits but more 50S particles than those prepared from bacteria supplemented with putrescine. The addition of the polyamine to putrescine-depleted cells induced a rapid change of the ribosomal profile. A similar effect could be obtained in vitro by equilibrium dialysis against a polyamine-containing solution. The ribosomal pattern obtained from starved bacteria was specific for polyamine deficiency. We conclude that the changes in ribosomal profiles upon restoration of putrescine levels in previously starved cells denote a shift of the equilibrium between 30S-50S couples and ribosomal subunits.

Escherichia coli

Translational control by RPL22L1-specific ribosomes enhances DNA repair and chemoresistance.

Ribosome heterogeneity has emerged as a regulatory layer in gene expression, yet its biological roles in cancers remain poorly characterized. Here, we identify RPL22L1, a paralog of the ribosomal protein RPL22, as a key modulator of DNA damage response (DDR) in colorectal cancer cells. DNA damage induces RPL22L1 upregulation and ribosomal incorporation, forming RPL22L1-specific ribosomes. Ribosome profiling reveals that RPL22L1-containing ribosomes preferentially translate mRNAs with highly structured 5' untranslated region (5'UTR). In particular, RPL22L1 enhances the translation of ATRX through a cap-independent mechanism. ATRX subsequently recruits DNA-PKcs to DNA damage sites, thereby enhancing the DNA repair capacity. RPL22L1 loss creates exploitable DDR vulnerabilities, sensitizing cancer cells to cisplatin and PARP inhibitors in vitro and in vivo. Collectively, these findings uncover a specialized ribosome-mediated translational program in DDR and highlight RPL22L1 as a potential therapeutic target in DDR-based cancer therapy.

DNA Repair

Developmental changes in the distribution of rat liver ribosomes in response to dietary manipulations.

Ribosomal profiles were prepared from liver postmitochondrial supernatant preparations of 2- to 21-day old and 60-day-old rats. Rats fed or starved for 3 or 15 to 16 hours were compared. Relative ribosomal distribution, calculated from planimetric analysis of the hepatic ribosomal profiles, was found to be approximately constant at all ages when fed animals were compared. In contrast, starvation for 15 hours resulted in profound changes in ribosomal distribution, i.e., a relative accumulation of oligosomes as compared with polysomes. This change in ribosomal distribution was most extensive in animals younger than 2 weeks of age, 60-day-old rats showing no effect. Refeeding of 6-day-old starved rats with diets containing 10% protein resulted in a prompt decrease in the planimetric yield of monsomes and disomes. Similarly, dietary protein was shown to be required to prevent accumulation of these two species. It is suggested that it is the response of the ribosomal cycle to starvation that changes with age rather than the maximum capacity for protein synthesis.

Age Factors

Comprehensive quantitative modeling of translation efficiency in a genome-reduced bacterium.

Translation efficiency has been mainly studied by ribosome profiling, which only provides an incomplete picture of translation kinetics. Here, we integrated the absolute quantifications of tRNAs, mRNAs, RNA half-lives, proteins, and protein half-lives with ribosome densities and derived the initiation and elongation rates for 475 genes (67% of all genes), 73 with high precision, in the bacterium Mycoplasma pneumoniae (Mpn). We found that, although the initiation rate varied over 160-fold among genes, most of the known factors had little impact on translation efficiency. Local codon elongation rates could not be fully explained by the adaptation to tRNA abundances, which varied over 100-fold among tRNA isoacceptors. We provide a comprehensive quantitative view of translation efficiency, which suggests the existence of unidentified mechanisms of translational regulation in Mpn.

RNA, Transfer

Effects of light on chloroplast translation in Marchantia polymorpha are similar to those in angiosperms and are not influenced by light-independent chlorophyll synthesis.

Translation of the chloroplast psbA mRNA in angiosperms is activated by photodamage of its gene product, the D1 subunit of photosystem II (PSII), providing nascent D1 for PSII repair. The involvement of chlorophyll in the regulatory mechanism has been suggested due to the regulatory roles of proteins proposed to mediate chlorophyll/D1 transactions and the fact that chlorophyll is synthesized only in the light in angiosperms. We used ribosome profiling and RNA-seq to address whether the effects of light on chloroplast translation are conserved in the liverwort Marchantia (Marchantia polymorpha), which synthesizes chlorophyll in both the dark and the light. As in angiosperms, ribosome occupancy on psbA mRNA decreased rapidly upon shifting plants to the dark and was rapidly restored upon a transfer back to the light, whereas ribosome occupancy on other chloroplast mRNAs changed very little. The results were similar in a Marchantia mutant unable to synthesize chlorophyll in the dark. Those results, in conjunction with pulse-labeling data, suggest that light elicits a plastome-wide activation of translation elongation and a specific increase in psbA translation initiation in Marchantia, as in angiosperms. These findings show that light regulates chloroplast translation similarly in vascular and non-vascular plants, and that constitutive chlorophyll synthesis does not affect light-regulated psbA translation initiation. Additionally, the translational outputs of chloroplast genes are similar in Marchantia and angiosperms but result from differing contributions of mRNA abundance and translational efficiencies. This adds to the evidence that chloroplast mRNA abundance and translational efficiencies co-evolve under selection to maintain protein outputs.

Chloroplasts

Beyond Canonical Neoantigens: Emerging Technologies for Identification of Noncanonical Antigens and Implications for Personalized Cancer Vaccines.

Over the past decade, advances in sequencing technologies and computational pipelines enabled the development of personalized cancer vaccines (PCVs). Current PCV strategies primarily target cancer neoantigens generated by non-synonymous DNA mutations, which can result in altered amino acid sequences capable of eliciting tumor-specific immune responses. More recently, a distinct class of tumor-specific antigens (TSA), termed noncanonical or cryptic antigens, has emerged as an additional source of immunogenic targets. Unlike canonical neoantigens, noncanonical antigens typically cannot be identified by tumor/normal whole-exome sequencing, as they do not arise from classical DNA mutations. Instead, they are often associated with less well recognized and/or aberrant processes in the pathways from DNA to human leukocyte antigen (HLA)-presented peptides. Examples include transposable elements, circular RNA, translation of alternative open reading frames and/or long non-coding RNA, among others. Emerging evidence suggests that noncanonical antigens represent a substantial portion of the tumor-specific immunopeptidome and, similar to canonical neoantigens, are absent during thymic selection and can evade central tolerance and elicit T cell responses. Technological advances have increasingly facilitated the identification of noncanonical antigens. Long-read RNA sequencing reveals noncanonical transcripts by improving transcriptome assembly, while ribosome profiling provides genome-wide maps of actively translated regions, facilitating the discovery of peptides from aberrant translation events. Specialized molecular approaches enable enrichment and sequencing of circular RNAs, and immunopeptidomics using mass spectrometry allows for direct characterization of HLA-presented peptides. Together, these technological advances have led to an increasing interest in prioritizing and targeting noncanonical antigens in the next generation of PCVs. This review provides an overview of the diverse origins of TSAs beyond classical neoantigens and discusses emerging approaches that may enable the integration of these antigens in future clinical trials.

circular RNA

IRES-like element-mediated translation of vsp1S4(-) suppresses BmCPV replication via RNAi antagonism.

Double-stranded RNA (dsRNA) viruses are thought to express proteins exclusively from their sense strand, while the antisense strand serves primarily as a replication template. Whether the antisense strand harbors hidden coding potential remains largely unexplored. Here, by integrating ribosome profiling and mass spectrometry, we identify a conserved 78-amino acid microprotein, vsp1S4(-), encoded by an antisense small open reading frame (sORFs) of the Bombyx mori cypovirus (BmCPV) genome. We demonstrate that vsp1S4(-) translation is driven by a previously unrecognized IRES-like element. Functional characterizations reveal that vsp1S4(-) localizes to the plasma membrane and acts as a negative regulator of viral replication. Mechanistically, vsp1S4(-) interacts directly with the viral RNAi suppressor NSP8, competitively disrupting the NSP8-AGO2 complex. This action restores the host's antiviral RNAi response, thereby limiting viral proliferation. Our findings challenge the conventional view of dsRNA virus coding capacity, unveil a novel viral immune evasion and replication control mechanism, and highlight antisense-encoded microproteins as potential targets for antiviral therapy.

Animals

PKD1 upstream open reading frames affect Polycystin-1 expression and polycystic kidney disease phenotypes.

Autosomal dominant polycystic kidney disease (ADPKD) accounts for 5%-10% of prevalent end-stage kidney failure (ESKD). ADPKD cysts result from a loss of sufficient functional expression of PKD1/Polycystin-1 (PC1) in approximately 80% of families. Kidney disease severity correlates with the extent to which PC1 dosage is reduced below a critical level, and evidence suggests therapeutic benefit from increasing PC1 expression in these conditions. Upstream open reading frame (uORF) translation can reduce translation of a protein's coding sequence. Ribosome profiling data and bioinformatic predictions suggested the presence of conserved PKD1 uORFs, so we sought to explore their biological role. We generated luciferase reporters and two humanized PKD1 5' UTR mouse models with or without single nucleotide edits removing uORF start codons (ΔuORF) to define active uORFs and test their impact on PC1 translation. PKD1 uORF start codons can robustly initiate translation, and ΔuORF conveys a 2-4 fold increase in PC1 protein expression and resultant prevention of kidney cysts in Dnajb11 as well as in Pkd1 missense models. PKD1 uORF1-blocking steric antisense oligonucleotides (ASOs) substantially increase PC1 expression in vitro. PKD1 uORFs play an important role in the low basal expression of WT PKD1, and their inhibition represents an opportunity to therapeutically increase PC1 translation in polycystic kidney and liver disease resulting from reduced dosage of PC1.

Animals

Hidden proteins encoded by non-canonical open reading frames: A review.

There is increasing evidence that translation is not limited to annotated protein-coding genes. Ribosome profiling sequencing, mass spectrometry-based proteomics, and immunopeptidomics have identified the productive translation of non-canonical open reading frames (ORFs). This suggests that the functional proteome includes not only conserved proteins but also proteins hidden in non-coding RNAs and de novo proteins. Some of these translated products are functional peptides, while others may be non-functional, potentially arising from evolutionary events. Several non-canonical ORF-encoded peptides have been found to regulate multiple physiological and pathological functions, particularly in cancer, immunity, and inflammation, indicating that they have potential as biomarkers and novel therapeutic targets. To better understand the diversity of functional peptides and translated non-canonical ORFs based on existing data, we summarize their classification according to transcriptional features and supporting evidence, including non-canonical ORFs located in ncRNAs and canonical mRNAs. This review provides a concise summary of the origin, discovery methods, and classification of non-canonical ORFs. It offers insights into the origins and functions of non-canonical ORF-encoded peptides from an evolutionary perspective, while also exploring the biological functions and regulatory mechanisms of these non-canonical ORF-encoded hidden proteins in tumorigenesis and progression.

Open Reading Frames