mRNA Secondary Structure Research Articles

Crystal structure of the single-stranded RNA binding protein HutP from Geobacillus thermodenitrificans

RNA binding proteins control gene expression by the attenuation/antitermination mechanism. HutP is an RNA binding antitermination protein. It regulates the expression of hut operon when it binds with RNA by modulating the secondary structure of single-stranded hut mRNA. HutP necessitates the presence of l-histidine and divalent metal ion to bind with RNA. Herein, we report the crystal structures of ternary complex (HutP–l-histidine–Mg2+) and EDTA (0.5M) treated ternary complex (HutP–l-histidine–Mg2+), solved at 1.9Å and 2.5Å resolutions, respectively, from Geobacillus thermodenitrificans. The addition of 0.5M EDTA does not affect the overall metal-ion mediated ternary complex structure and however, the metal ions at the non-specific binding sites are chelated, as evidenced from the results of structural features.

The transcriptional activities and cellular localization of the human estrogen receptor alpha are affected by the synonymous Ala87 mutation

Until recently, synonymous mutations (which do not change amino acids) have been much neglected. Some evidence suggests that this kind of mutations could affect mRNA secondary structure or stability, translation kinetics and protein structure. To explore deeper the role of synonymous mutations, we studied their consequence on the functional activity of the estrogen receptor alpha (ERα). The ERα is a ligand-inducible transcription factor that orchestrates pleiotropic cellular effects, at both genomic and non-genomic levels in response to estrogens. In this work we analyzed in transient transfection experiments, the activity of ERα carrying the synonymous mutation Ala87, a polymorphism involving about 5–10% of the population. In comparison to the wild type receptor, our results show that ERαA87 mutation reduces the transactivation efficiency of ERα on an ERE reporter gene while its expression level remains similar. This mutation enhances 4-OHT-induced transactivation of ERα on an AP1 reporter gene. Finally, the mutation affects the subcellular localization of ERα in a cell type specific manner. It enhances the cytoplasmic location of ERα without significant changes in non-genomic effects of E2. The functional alteration of the ERαA87 determined in this work highlights the relevance of synonymous mutations for biomedical and pharmacological points of view.

Deciphering the rules by which dynamics of mRNA secondary structure affect translation efficiency in Saccharomyces cerevisiae

Messenger RNA (mRNA) secondary structure decreases the elongation rate, as ribosomes must unwind every structure they encounter during translation. Therefore, the strength of mRNA secondary structure is assumed to be reduced in highly translated mRNAs. However, previous studies in vitro reported a positive correlation between mRNA folding strength and protein abundance. The counterintuitive finding suggests that mRNA secondary structure affects translation efficiency in an undetermined manner. Here, we analyzed the folding behavior of mRNA during translation and its effect on translation efficiency. We simulated translation process based on a novel computational model, taking into account the interactions among ribosomes, codon usage and mRNA secondary structures. We showed that mRNA secondary structure shortens ribosomal distance through the dynamics of folding strength. Notably, when adjacent ribosomes are close, mRNA secondary structures between them disappear, and codon usage determines the elongation rate. More importantly, our results showed that the combined effect of mRNA secondary structure and codon usage in highly translated mRNAs causes a short ribosomal distance in structural regions, which in turn eliminates the structures during translation, leading to a high elongation rate. Together, these findings reveal how the dynamics of mRNA secondary structure coupling with codon usage affect translation efficiency.

ABCB1 haplotypes are associated with P-gp activity and affect a major molecular response in chronic myeloid leukemia patients treated with a standard dose of imatinib

Despite the high efficacy of imatinib mesylate (IM) treatment for chronic myeloid leukemia (CML) patients, some individuals develop resistance due to impaired bioavailability. It has been previously demonstrated that the haplotypes for ATP-binding cassette subfamily B member 1 (ABCB1)with c.1236C>T, c.3435C>T and c.2677G>T/A polymorphisms markedly affect the secondary structure of ABCB1 mRNA and its activity. These modifications may affect efflux transporter activity and response to treatment with IM. The aim of the present study was to investigate the influence of ABCB1 haplotypes on P-glycoprotein (P-gp) activity, IM plasma levels and IM response. In total, 28 chronic-phase CML patients treated with a standard dose of IM (400 mg/day) were studied. The patients were selected according to the haplotypes of ABCB1, with c.1236C>T, c.3435C>T and c.2677G>T polymorphisms, and were classified into two groups based on the presence of the mutated allele in each genotype for the three ABCB1 polymorphisms. In addition, expression of P-gp and breakpoint cluster region-abelson 1 (BCR-ABL1), ABCB1 and solute carrier family 22 member 1 (SLC22A1) mRNA were evaluated. The P-gp activity in the wild-type group was found to be higher than that in the mutated group (59.1 vs. 38.3%; P=0.001). Furthermore, the patients who did not achieve major molecular response (MMR) showed a higher rate of efflux mediated by P-gp when compared with individuals who achieved MMR (64.7 vs. 45.7%; P=0.001). All patients without MMR demonstrated effluxes of >60%. In addition, patients without MMR exhibited lower plasma concentrations of IM compared with those with MMR (0.51 vs. 1.42 μg/ml; P=0.001). Higher levels of SLC22A1 mRNA were observed in patients who achieved MMR and complete molecular response (P<0.05). In conclusion, the ABCB1 1236CT/3435CT/2677GT and 1236TT/3435TT/2677TT haplotypes are associated with reduced P-gp activity and MMR in chronic-phase CML patients treated with a standard dose of IM.

Cyanobacterial heat-shock response: role and regulation of molecular chaperones

Cyanobacteria constitute a morphologically diverse group of oxygenic photoautotrophic microbes which range from unicellular to multicellular, and non-nitrogen-fixing to nitrogen-fixing types. Sustained long-term exposure to changing environmental conditions, during their three billion years of evolution, has presumably led to their adaptation to diverse ecological niches. The ability to maintain protein conformational homeostasis (folding-misfolding-refolding or aggregation-degradation) by molecular chaperones holds the key to the stress adaptability of cyanobacteria. Although cyanobacteria possess several genes encoding DnaK and DnaJ family proteins, these are not the most abundant heat-shock proteins (Hsps), as is the case in other bacteria. Instead, the Hsp60 family of proteins, comprising two phylogenetically conserved proteins, and small Hsps are more abundant during heat stress. The contribution of the Hsp100 (ClpB) family of proteins and of small Hsps in the unicellular cyanobacteria (Synechocystis and Synechococcus) as well as that of Hsp60 proteins in the filamentous cyanobacteria (Anabaena) to thermotolerance has been elucidated. The regulation of chaperone genes by several cis-elements and trans-acting factors has also been well documented. Recent studies have demonstrated novel transcriptional and translational (mRNA secondary structure) regulatory mechanisms in unicellular cyanobacteria. This article provides an insight into the heat-shock response: its organization, and ecophysiological regulation and role of molecular chaperones, in unicellular and filamentous nitrogen-fixing cyanobacterial strains.

Molecular Dynamics Simulation of the Allosteric Regulation of eIF4A Protein from the Open to Closed State, Induced by ATP and RNA Substrates

BackgroundEukaryotic initiation factor 4A (eIF4A) plays a key role in the process of protein translation initiation by facilitating the melting of the 5′ proximal secondary structure of eukaryotic mRNA for ribosomal subunit attachment. It was experimentally postulated that the closed conformation of the eIF4A protein bound by the ATP and RNA substrates is coupled to RNA duplex unwinding to promote protein translation initiation, rather than an open conformation in the absence of ATP and RNA substrates. However, the allosteric process of eIF4A from the open to closed state induced by the ATP and RNA substrates are not yet fully understood.MethodologyIn the present work, we constructed a series of diplex and ternary models of the eIF4A protein bound by the ATP and RNA substrates to carry out molecular dynamics simulations, free energy calculations and conformation analysis and explore the allosteric properties of eIF4A.ResultsThe results showed that the eIF4A protein completes the conformational transition from the open to closed state via two allosteric processes of ATP binding followed by RNA and vice versa. Based on cooperative allosteric network analysis, the ATP binding to the eIF4A protein mainly caused the relative rotation of two domains, while the RNA binding caused the proximity of two domains via the migration of RNA bases in the presence of ATP. The cooperative binding of ATP and RNA for the eIF4A protein plays a key role in the allosteric transition.

Modified ‘one amino acid-one codon’ engineering of high GC content TaqII-coding gene from thermophilic Thermus aquaticus results in radical expression increase

BackgroundAn industrial approach to protein production demands maximization of cloned gene expression, balanced with the recombinant host’s viability. Expression of toxic genes from thermophiles poses particular difficulties due to high GC content, mRNA secondary structures, rare codon usage and impairing the host’s coding plasmid replication.TaqII belongs to a family of bifunctional enzymes, which are a fusion of the restriction endonuclease (REase) and methyltransferase (MTase) activities in a single polypeptide. The family contains thermostable REases with distinct specificities: TspGWI, TaqII, Tth111II/TthHB27I, TspDTI and TsoI and a few enzymes found in mesophiles. While not being isoschizomers, the enzymes exhibit amino acid (aa) sequence homologies, having molecular sizes of ~120 kDa share common modular architecture, resemble Type-I enzymes, cleave DNA 11/9 nt from the recognition sites, their activity is affected by S-adenosylmethionine (SAM).ResultsWe describe the taqIIRM gene design, cloning and expression of the prototype TaqII. The enzyme amount in natural hosts is extremely low. To improve expression of the taqIIRM gene in Escherichia coli (E. coli), we designed and cloned a fully synthetic, low GC content, low mRNA secondary structure taqIIRM, codon-optimized gene under a bacteriophage lambda (λ) P R promoter. Codon usage based on a modified ‘one amino acid–one codon’ strategy, weighted towards low GC content codons, resulted in approximately 10-fold higher expression of the synthetic gene. 718 codons of total 1105 were changed, comprising 65% of the taqIIRM gene. The reason for we choose a less effective strategy rather than a resulting in high expression yields ‘codon randomization’ strategy, was intentional, sub-optimal TaqII in vivo production, in order to decrease the high ‘toxicity’ of the REase-MTase protein.ConclusionsRecombinant wt and synthetic taqIIRM gene were cloned and expressed in E. coli. The modified ‘one amino acid–one codon’ method tuned for thermophile-coded genes was applied to obtain overexpression of the ‘toxic’ taqIIRM gene. The method appears suited for industrial production of thermostable ‘toxic’ enzymes in E. coli. This novel variant of the method biased toward increasing a gene’s AT content may provide economic benefits for industrial applications.

The mRNA Secondary Structure Impact on Translation of the KCNH2 Channel Protein

The gene KCNH2 encodes the K+ channel responsible for the repolarizing IKr current essential for normal cardiac physiology. Genetic, translational and post-translational perturbation of KCNH2 causes hereditary (Locus LQT2) and acquired Long QT Syndrome. Over 600 deleterious mutations have been identified as causes of LQT2, most leading to protein trafficking or assembly malfunction. Investigation has focused on the impact the changed amino acids have on the protein product, but little is known about mutational effects on mRNA function. KCNH2 is unusual in that its GC content is ∼66% and there is increased use of rare codons when compared to other human mRNA. Previous studies from our lab have shown that KCNH2 channel translation and trafficking efficiencies are dependent on “extra-coding” elements in the mRNA sequence. CM-KCNH2, a codon modified version of the KCNH2 mRNA, was made with reduced GC content (51%) and decreased use of rare codons, while maintaining the exact amino acid sequence. CM-KCNH2 protein product was a more efficiently trafficked protein, but was translated less efficiently when compared to native KCNH2 (NT-KCNH2). Whether codon-usage, GC content or mRNA 2° structure controls these processes is unknown. To address this, we employ a combined approach of ribosomal profiling to analyze precisely where translational slowing/pausing occurs, and RNA SHAPE analysis to investigate the secondary structure of the KCNH2 mRNA. We hypothesize that the indicated sites of ribosomal pausing, caused by secondary mRNA structure and rare codon placement, allow for correct secondary and tertiary structure formation of the KCNH2 protein, indicating that mRNA impacts not only the primary structure of a protein, but through ‘extra-coding’ features, impacts the folding efficiency and function of the protein product. Ongoing studies will have implications for evaluation and treatment of hereditary and acquired arrhythmias.

Influence of Downstream DNA/RNA Structure on Intersubunit Rotation of the Ribosome

The ribosome, composed of two subunits (30S and 50S), is a molecular machine that translates the genetic code inscribed in the mRNA into a protein sequence through repetitive cycles of aminoacyl-tRNA selection, peptide bond formation and translocation. However, the presence of mRNA secondary or tertiary structures change the rhythm of translation and in many instances may influence reading frame maintenance. We have shown with smFRET and SAXS that mRNA structures when present outside the entry tunnel of the ribosome induce a hyper-rotated state of the 30S subunit and render the L1 stalk primarily in an open configuration. In this study, we seek to investigate the distance dependence between downstream DNA oligos and the mRNA entry tunnel of the ribosome to determine how close structures need to be before they are sensed by the ribosome. A series of DNA-RNA hybrids were created to assemble with fluorescently labeled ribosomes and intersubunit rotation was followed by single molecule FRET.

Modeling Unpairing Costs for Fast Computation of the Net Binding Free Energy of an Oligo to an mRNA Target

View Large Image | View Hi-Res Image | Download PowerPoint SlideBinding of a small RNA (miRNA, snoRNA, siRNA) often regulates the expression of messenger RNA. The net free energy consists of the contribution of binding, which is well parameterized and fast to compute, offset by the penalty to remove pre-existing mRNA secondary structure, which is O(N3) to compute with previous methods. Here we introduce a base-composition-dependent model for the mean free energy to open the target which runs in O(N). Implemented in our BINDIGONET algorithm, this model for the net binding free energy of spliceosomal U1-snRNA to real and decoy donor splice sites enhances accuracy.

Finding the target sites of RNA‐binding proteins

RNA–protein interactions differ from DNA–protein interactions because of the central role of RNA secondary structure. Some RNA-binding domains (RBDs) recognize their target sites mainly by their shape and geometry and others are sequence-specific but are sensitive to secondary structure context. A number of small- and large-scale experimental approaches have been developed to measure RNAs associated in vitro and in vivo with RNA-binding proteins (RBPs). Generalizing outside of the experimental conditions tested by these assays requires computational motif finding. Often RBP motif finding is done by adapting DNA motif finding methods; but modeling secondary structure context leads to better recovery of RBP-binding preferences. Genome-wide assessment of mRNA secondary structure has recently become possible, but these data must be combined with computational predictions of secondary structure before they add value in predicting in vivo binding. There are two main approaches to incorporating structural information into motif models: supplementing primary sequence motif models with preferred secondary structure contexts (e.g., MEMERIS and RNAcontext) and directly modeling secondary structure recognized by the RBP using stochastic context-free grammars (e.g., CMfinder and RNApromo). The former better reconstruct known binding preferences for sequence-specific RBPs but are not suitable for modeling RBPs that recognize shape and geometry of RNAs. Future work in RBP motif finding should incorporate interactions between multiple RBDs and multiple RBPs in binding to RNA. WIREs RNA 2014, 5:111–130. doi: 10.1002/wrna.1201

Dynamic Recognition of the mRNA Cap by Saccharomyces cerevisiae eIF4E

Recognition of the mRNA 5' m⁷G(5')ppp(5')N cap is key to translation initiation for most eukaryotic mRNAs. The cap is bound by the eIF4F complex, consisting of a cap-binding protein (eIF4E), a "scaffold" protein (eIF4G), and an RNA helicase (eIF4A). As a central early step in initiation, regulation of eIF4F is crucial for cellular viability. Although the structure and function of eIF4E have been defined, a dynamic mechanistic picture of its activity at the molecular level in the eIF4F·mRNA complex is still unavailable. Here, using single-molecule fluorescence, we measured the effects of Saccharomyces cerevisiae eIF4F factors, mRNA secondary structure, and the poly(A)-binding protein Pab1p on eIF4E-mRNA binding dynamics. Our data provide an integrated picture of how eIF4G and mRNA structure modulate eIF4E-mRNA interaction, and uncover an eIF4G- and poly(A)-independent activity of poly(A)-binding protein that prolongs the eIF4E·mRNA complex lifetime.

Genetic Code-guided Protein Synthesis and Folding in Escherichia coli

Universal genetic codes are degenerated with 61 codons specifying 20 amino acids, thus creating synonymous codons for a single amino acid. Synonymous codons have been shown to affect protein properties in a given organism. To address this issue and explore how Escherichia coli selects its "codon-preferred" DNA template(s) for synthesis of proteins with required properties, we have designed synonymous codon libraries based on an antibody (scFv) sequence and carried out bacterial expression and screening for variants with altered properties. As a result, 342 codon variants have been identified, differing significantly in protein solubility and functionality while retaining the identical original amino acid sequence. The soluble expression level varied from completely insoluble aggregates to a soluble yield of ~2.5 mg/liter, whereas the antigen-binding activity changed from no binding at all to a binding affinity of > 10(-8) m. Not only does our work demonstrate the involvement of genetic codes in regulating protein synthesis and folding but it also provides a novel screening strategy for producing improved proteins without the need to substitute amino acids.

Novel mutations of KCNQ1 in Long QT syndrome

Autosomal recessive Long QT syndrome is characterized by prolonged QTc along with congenital bilateral deafness depends on mutations in K(+) channel genes. Afamily of a Long QT syndrome proband from India has been identified with novel indel variations. The molecular study of the proband revealed 4 novel indel variations in KCNQ1. In-silico analysis revealed the intronic variations has led to a change in the secondary structure of mRNA and splice site variations. The exonic variations leads to frameshift mutations. DNA analysis of the available family members revealed a carrier status. It is thus predicted that the variations may lead to a change in the position of the splicing enhancer/inhibitor in KCNQ1 leading to the formation of a truncated S2-S3 fragment of KCNQ1 transmembrane protein in cardiac cells as well as epithelial cells of inner ear leading to deafness and aberrant repolarization causing prolonged QTc.

P.6.6 Splice switching Antisense Oligonucleotides: A “kink” in the design?

Spinal Muscular Atrophy (SMA) is a neuromuscular disease caused by the lack of a functional SMN1 gene product, resulting in loss of motor neurons in the anterior horn. This leads to muscle wastage and early death. However, SMA phenotype and disease progression are modified by the presence of one or more copies of SMN2 a near identical gene that potentially encodes the same protein. While multiple copies of SMN2 can reduce disease severity, a base change in exon 7 alters splicing fidelity leading to exon 7 exclusion from the majority of transcripts and a non-functional protein. Antisense Oligonucleotides (AOs) to enhance recognition and retention of exon 7 in the SMN2 transcript have been explored as a therapy for SMA. Others have shown that a 20 base AO sequence targeted to ISS-N1 , a splice silencing domain within intron 7, can promote exon 7 recognition and retention in the mature SMN2 transcript. Five bases added to the 5′ end of this sequence further improves levels of exon 7 retention in SMA cells in vitro. However, contrary to expectation, we found that some mismatched AOs are as effective as the previously reported sequences. We tested a series of mismatch AO sequences for both Anti-ISS-N1 and Anti-ISS-N1+5. Different mismatches were evaluated, including the deletion of either a Uracil (U) or Adenine (A) base, and an U → A base change at various positions along the oligomer sequence. We show that AOs carrying the deletion of an Adenine base or the transversion from U → A to be more efficient, inducing 100% exon inclusion at an AO concentration of 10 nM. We speculate that the AOs are acting in a similar way to microRNAs, where selected mismatches between the mRNA target and the oligomer are biologically active. Perhaps the mismatch creates a “kink” or change in the mRNA secondary structure, ultimately affecting mRNA and potentially protein binding. Further studies aim to investigate this mechanism and its therapeutic potential for SMA and other disorders.

One antitoxin—two functions: SR4 controls toxin mRNA decay and translation

Type I toxin–antitoxin systems encoded on bacterial chromosomes became the focus of research during the past years. However, little is known in terms of structural requirements, kinetics of interaction with their targets and regulatory mechanisms of the antitoxin RNAs. Here, we present a combined in vitro and in vivo analysis of the bsrG/SR4 type I toxin–antitoxin system from Bacillus subtilis. The secondary structures of SR4 and bsrG mRNA and of the SR4/bsrG RNA complex were determined, apparent binding rate constants calculated and functional segments required for complex formation narrowed down. The initial contact between SR4 and its target was shown to involve the SR4 terminator loop and loop 3 of bsrG mRNA. Additionally, a contribution of the stem of SR4 stem-loop 3 to target binding was found. On SR4/bsrG complex formation, a 4 bp double-stranded region sequestering the bsrG Shine Dalgarno (SD) sequence was extended to 8 bp. Experimental evidence was obtained that this extended region caused translation inhibition of bsrG mRNA. Therefore, we conclude that SR4 does not only promote degradation of the toxin mRNA but also additionally inhibit its translation. This is the first case of a dual-acting antitoxin RNA.

Composability of regulatory sequences controlling transcription and translation in Escherichia coli

The inability to predict heterologous gene expression levels precisely hinders our ability to engineer biological systems. Using well-characterized regulatory elements offers a potential solution only if such elements behave predictably when combined. We synthesized 12,563 combinations of common promoters and ribosome binding sites and simultaneously measured DNA, RNA, and protein levels from the entire library. Using a simple model, we found that RNA and protein expression were within twofold of expected levels 80% and 64% of the time, respectively. The large dataset allowed quantitation of global effects, such as translation rate on mRNA stability and mRNA secondary structure on translation rate. However, the worst 5% of constructs deviated from prediction by 13-fold on average, which could hinder large-scale genetic engineering projects. The ease and scale this of approach indicates that rather than relying on prediction or standardization, we can screen synthetic libraries for desired behavior.

Deciphering the rules by which 5′-UTR sequences affect protein expression in yeast

The 5'-untranslated region (5'-UTR) of mRNAs contains elements that affect expression, yet the rules by which these regions exert their effect are poorly understood. Here, we studied the impact of 5'-UTR sequences on protein levels in yeast, by constructing a large-scale library of mutants that differ only in the 10 bp preceding the translational start site of a fluorescent reporter. Using a high-throughput sequencing strategy, we obtained highly accurate measurements of protein abundance for over 2,000 unique sequence variants. The resulting pool spanned an approximately sevenfold range of protein levels, demonstrating the powerful consequences of sequence manipulations of even 1-10 nucleotides immediately upstream of the start codon. We devised computational models that predicted over 70% of the measured expression variability in held-out sequence variants. Notably, a combined model of the most prominent features successfully explained protein abundance in an additional, independently constructed library, whose nucleotide composition differed greatly from the library used to parameterize the model. Our analysis reveals the dominant contribution of the start codon context at positions -3 to -1, mRNA secondary structure, and out-of-frame upstream AUGs (uAUGs) to phenotypic diversity, thereby advancing our understanding of how protein levels are modulated by 5'-UTR sequences, and paving the way toward predictably tuning protein expression through manipulations of 5'-UTRs.

Glucocorticoid Receptor Ligands ModulateCardiovirus encephalomyocarditis virusInternal Ribosome Entry Site Activity

The use of small molecules to modulate cellular processes is a powerful approach to investigate gene function as a complement to genetic approaches. The discovery and characterization of compounds that modulate translation initiation, the rate-limiting step of protein synthesis, is important both to provide tool compounds to explore this fundamental biological process and to further evaluate protein synthesis as a therapeutic target. While most messenger ribonucleic acids (mRNAs) recruit ribosomes via their 5' cap, some viral and cellular mRNAs initiate protein synthesis via an alternative "cap-independent" mechanism utilizing internal ribosome entry sites (IRES) elements, which are complex mRNA secondary structures, localized within the 5' nontranslated region of the mRNA upstream of the AUG start codon. This report describes the design of a functional, high throughput screen of small molecules miniaturized into a 1,536-well format and performed using the luciferase reporter gene under control of the viral Cardiovirus encephalomyocarditis virus (EMCV) IRES element to identify nontoxic compounds modulating translation initiated from the EMCV IRES. One activating compound, validated in a dose response manner, has previously been shown to bind the glucocorticoid receptor (GR). Subsequent testing of additional GR modulators further supported this as the possible mechanism of action. Detailed characterization of this compound activity supported the notion that this was due to an effect at the level of translation.

Genome-Wide Patterns of Codon Bias Are Shaped by Natural Selection in the Purple Sea Urchin,Strongylocentrotus purpuratus

Codon usage bias has been documented in a wide diversity of species, but the relative contributions of mutational bias and various forms of natural selection remain unclear. Here, we describe for the first time genome-wide patterns of codon bias at 4623 genes in the purple sea urchin, Strongylocentrotus purpuratus. Preferred codons were identified at 18 amino acids that exclusively used G or C at third positions, which contrasted with the strong AT bias of the genome (overall GC content is 36.9%). The GC content of third positions and coding regions exhibited significant correlations with the magnitude of codon bias. In contrast, the GC content of introns and flanking regions was indistinguishable from the genome-wide background, which suggested a limited contribution of mutational bias to synonymous codon usage. Five distinct clusters of genes were identified that had significantly different synonymous codon usage patterns. A significant correlation was observed between codon bias and mRNA expression supporting translational selection, but this relationship was driven by only one highly biased cluster that represented only 8.6% of all genes. In all five clusters preferred codons were evolutionarily conserved to a similar degree despite differences in their synonymous codon usage distributions and magnitude of codon bias. The third positions of preferred codons in two codon usage groups also paired significantly more often in stems than in loops of mRNA secondary structure predictions, which suggested that codon bias might also affect mRNA stability. Our results suggest that mutational bias has played a minor role in determining codon bias in S. purpuratus and that preferred codon usage may be heterogeneous across different genes and subject to different forms of natural selection.