Long read subcellular fractionation and sequencing reveals the translational fate of full-length mRNA isoforms during neuronal differentiation.

Alternative splicing (AS) alters the cis-regulatory landscape of mRNA isoforms leading to transcripts with distinct localization, stability and translational efficiency. To rigorously investigate mRNA isoform-specific ribosome association, we generated subcellular fractionation and sequencing (Frac-seq) libraries using both conventional short reads and long reads from human embryonic stem cells (ESC) and neural progenitor cells (NPC) derived from the same ESC. We performed de novo transcriptome assembly from high-confidence long reads from cytosolic, monosomal, light and heavy polyribosomal fractions and quantified their abundance using short reads from their respective subcellular fractions. Thousands of transcripts in each cell type exhibited association with particular subcellular fractions relative to the cytosol. Of the multi-isoform genes, 27% and 19% exhibited significant differential isoform sedimentation in ESC and NPC respectively. Alternative promoter usage and internal exon skipping accounted for the majority of differences between isoforms from the same gene. Random forest classifiers implicated coding sequence (CDS) and UTR lengths as important determinants of isoform-specific sedimentation profiles, and motif analyses reveal potential cell type-specific and subcellular fraction-associated RNA-binding protein signatures. Taken together our data demonstrate that alternative mRNA processing within the CDS and UTRs impacts the translational control of mRNA isoforms during stem cell differentiation, and highlights the utility of using a novel long-read sequencing-based method to study translational control.

Targeting IGF2BP3 enhances antileukemic effects of menin-MLL inhibition in MLL-AF4 leukemia.

ABSTRACT: RNA-binding proteins (RBPs) are emerging as a novel class of therapeutic targets in cancer, including in leukemia, given their important role in posttranscriptional gene regulation, and have the unexplored potential to be combined with existing therapies. The RBP insulin-like growth factor 2 messenger RNA-binding protein 3 (IGF2BP3) has been found to be a critical regulator of MLL-AF4 leukemogenesis and represents a promising therapeutic target. Here, we study the combined effects of targeting IGF2BP3 and menin-MLL interaction in MLL-AF4-driven leukemia in vitro and in vivo, using genetic inhibition with CRISPR-Cas9-mediated deletion of Igf2bp3 and pharmacologic inhibition of the menin-MLL interaction with multiple commercially available inhibitors. Depletion of Igf2bp3 sensitized MLL-AF4 leukemia to the effects of menin-MLL inhibition on cell growth and leukemic initiating cells in vitro. Mechanistically, we found that both Igf2bp3 depletion and menin-MLL inhibition led to increased differentiation in vitro and in vivo, seen in functional readouts and by gene expression analyses. IGF2BP3 knockdown had a greater effect on increasing survival and attenuating disease than pharmacologic menin-MLL inhibition with small molecule MI-503 alone and showed enhanced antileukemic effects in combination. Our work shows that IGF2BP3 is an oncogenic amplifier of MLL-AF4-mediated leukemogenesis and a potent therapeutic target, providing a paradigm for targeting leukemia at both the transcriptional and posttranscriptional level.

An intronic RNA element modulates Factor VIII exon-16 splicing.

Pathogenic variants in the human Factor VIII (F8) gene cause Hemophilia A (HA). Here, we investigated the impact of 97 HA-causing single-nucleotide variants on the splicing of 11 exons from F8. For the majority of F8 exons, splicing was insensitive to the presence of HA-causing variants. However, splicing of several exons, including exon-16, was impacted by variants predicted to alter exonic splicing regulatory sequences. Using exon-16 as a model, we investigated the structure-function relationship of HA-causing variants on splicing. Intriguingly, RNA chemical probing analyses revealed a three-way junction structure at the 3'-end of intron-15 (TWJ-3-15) capable of sequestering the polypyrimidine tract. We discovered antisense oligonucleotides (ASOs) targeting TWJ-3-15 partially rescue splicing-deficient exon-16 variants by increasing accessibility of the polypyrimidine tract. The apical stem loop region of TWJ-3-15 also contains two hnRNPA1-dependent intronic splicing silencers (ISSs). ASOs blocking these ISSs also partially rescued splicing. When used in combination, ASOs targeting both the ISSs and the region sequestering the polypyrimidine tract, fully rescue pre-mRNA splicing of multiple HA-linked variants of exon-16. Together, our data reveal a putative RNA structure that sensitizes F8 exon-16 to aberrant splicing.

Distinct oncogenic phenotypes in hematopoietic specific deletions of Trp53.

Loss of function in the tumor suppressor gene TP53 is the most common alteration seen in human cancer. In mice, P53 deletion in all cells leads predominantly to the development of T-cell lymphomas, followed by B-cell lymphomas, sarcomas and teratomas. In order to dissect the role of P53 in the hematopoietic system, we generated and analyzed two different mouse models deficient for P53. A pan-hematopoietic P53 deletion mouse was created using Vav1-Cre based deletion; and a B-cell-specific deletion mouse was created using a CD19-Cre based deletion. The Vav1-P53CKO mice predominantly developed T-cell malignancies in younger mice, and myeloid malignancies in older mice. In T-cell malignancies, there was accelerated thymic cell maturation with overexpression of Notch1 and its downstream effectors. CD19-P53CKO mice developed marginal zone expansion in the spleen, followed by marginal zone lymphoma, some of which progressed to diffuse large B-cell lymphomas. Interestingly, marginal zone and diffuse large B-cell lymphomas had a unique gene expression signature characterized by activation of the PI3K pathway, compared with wild type marginal zone or follicular cells of the spleen. This study demonstrates lineage specific P53 deletion leading to distinct phenotypes secondary to unique gene expression programs set in motion.

Isoform-specific translational control is evolutionarily conserved in primates.

Alternative splicing (AS) alters messenger RNA (mRNA) coding capacity, localization, stability, and translation. Here we use comparative transcriptomics to identify cis-acting elements coupling AS to translational control (AS-TC). We sequenced total cytosolic and polyribosome-associated mRNA from human, chimpanzee, and orangutan induced pluripotent stem cells (iPSCs), revealing thousands of transcripts with splicing differences between subcellular fractions. We found both conserved and species-specific polyribosome association patterns for orthologous splicing events. Intriguingly, alternative exons with similar polyribosome profiles between species have stronger sequence conservation than exons with lineage-specific ribosome association. These data suggest that sequence variation underlies differences in the polyribosome association. Accordingly, single nucleotide substitutions in luciferase reporters designed to model exons with divergent polyribosome profiles are sufficient to regulate translational efficiency. We used position specific weight matrices to interpret exons with species-specific polyribosome association profiles, finding that polymorphic sites frequently alter recognition motifs for trans-acting RNA binding proteins. Together, our results show that AS can regulate translation by remodeling the cis-regulatory landscape of mRNA isoforms.

Rescue of blood coagulation Factor VIII exon-16 mis-splicing by antisense oligonucleotides.

The human Factor VIII ( F8 ) protein is essential for the blood coagulation cascade and specific F8 mutations cause the rare bleeding disorder Hemophilia A (HA). Here, we investigated the impact of HA-causing single-nucleotide mutations on F8 pre-mRNA splicing. We found that 14/97 (∼14.4%) coding sequence mutations tested in our study induced exon skipping. Splicing patterns of 4/11 (∼36.4%) F8 exons tested were especially sensitive to the presence of common disease-causing mutations. RNA-chemical probing analyses revealed a three-way junction structure at the 3' end of intron 15 (TWJ-3-15). TWJ-3-15 sequesters the polypyrimidine tract, a key determinant of 3' splice site strength. Using exon-16 of the F8 gene as a model, we designed specific antisense oligonucleotides (ASOs) that target TWJ-3-15 and identified three that promote the splicing of F8 exon-16. Interaction of TWJ-3-15 with ASOs increases accessibility of the polypyrimidine tract and inhibits the binding of hnRNPA1-dependent splicing silencing factors. Moreover, ASOs targeting TWJ-3-15 rescue diverse splicing-sensitive HA-causing mutations, most of which are distal to the 3' splice site being impacted. The TWJ-3-15 structure and its effect on mRNA splicing provide a model for HA etiology in patients harboring specific F8 mutations and provide a framework for precision RNA-based HA therapies.

A potential role for SARS-CoV-2 small viral RNAs in targeting host microRNAs and modulating gene expression.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) causes coronavirus disease (COVID-19) in humans, with symptoms ranging from mild to severe, including fatality. The molecular mechanisms surrounding the effects of viral infection on the host RNA machinery remain poorly characterized. We used a comparative transcriptomics approach to investigate the effects of SARS-CoV-2 infection on the host mRNA and sRNA expression machinery in a human lung epithelial cell line (Calu-3) and an African green monkey kidney cell line (Vero-E6). Upon infection, we observed global changes in host gene expression and differential expression of dozens of host miRNAs, many with known links to viral infection and immune response. Additionally, we discovered an expanded landscape of more than a hundred SARS-CoV-2-derived small viral RNAs (svRNAs) predicted to interact with differentially expressed host mRNAs and miRNAs. svRNAs are derived from distinct regions of the viral genome and sequence signatures suggest they are produced by a non-canonical biogenesis pathway. 52 of the 67 svRNAs identified in Calu-3 cells are predicted to interact with differentially expressed miRNAs, with many svRNAs having multiple targets. Accordingly, we speculate that these svRNAs may play a role in SARS-CoV-2 propagation by modulating post-transcriptional gene regulation, and that methods for antagonizing them may have therapeutic value.

The RNA-binding protein IGF2BP3 is critical for MLL-AF4-mediated leukemogenesis.

Despite recent advances in therapeutic approaches, patients with MLL-rearranged leukemia still have poor outcomes. Here, we find that the RNA-binding protein IGF2BP3, which is overexpressed in MLL-translocated leukemia, strongly amplifies MLL-Af4-mediated leukemogenesis. Deletion of Igf2bp3 significantly increases the survival of mice with MLL-Af4-driven leukemia and greatly attenuates disease, with a minimal impact on baseline hematopoiesis. At the cellular level, MLL-Af4 leukemia-initiating cells require Igf2bp3 for their function in leukemogenesis. At the molecular level, IGF2BP3 regulates a complex posttranscriptional operon governing leukemia cell survival and proliferation. IGF2BP3-targeted mRNA transcripts include important MLL-Af4-induced genes, such as those in the Hoxa locus, and the Ras signaling pathway. Targeting of transcripts by IGF2BP3 regulates both steady-state mRNA levels and, unexpectedly, pre-mRNA splicing. Together, our findings show that IGF2BP3 represents an attractive therapeutic target in this disease, providing important insights into mechanisms of posttranscriptional regulation in leukemia.

A method for campus-wide SARS-CoV-2 surveillance at a large public university.

The systematic screening of asymptomatic and pre-symptomatic individuals is a powerful tool for controlling community transmission of infectious disease on college campuses. Faced with a paucity of testing in the beginning of the COVID-19 pandemic, many universities developed molecular diagnostic laboratories focused on SARS-CoV-2 diagnostic testing on campus and in their broader communities. We established the UC Santa Cruz Molecular Diagnostic Lab in early April 2020 and began testing clinical samples just five weeks later. Using a clinically-validated laboratory developed test (LDT) that avoided supply chain constraints, an automated sample pooling and processing workflow, and a custom laboratory information management system (LIMS), we expanded testing from a handful of clinical samples per day to thousands per day with the testing capacity to screen our entire campus population twice per week. In this report we describe the technical, logistical, and regulatory processes that enabled our pop-up lab to scale testing and reporting capacity to thousands of tests per day.

Alveolar progenitor differentiation and lactation depends on paracrine inhibition of notch via ROBO1/CTNNB1/JAG1.

In the mammary gland, how alveolar progenitor cells are recruited to fuel tissue growth with each estrus cycle and pregnancy remains poorly understood. Here, we identify a regulatory pathway that controls alveolar progenitor differentiation and lactation by governing Notch activation in mouse. Loss of Robo1 in the mammary gland epithelium activates Notch signaling, which expands the alveolar progenitor cell population at the expense of alveolar differentiation, resulting in compromised lactation. ROBO1 is expressed in both luminal and basal cells, but loss of Robo1 in basal cells results in the luminal differentiation defect. In the basal compartment, ROBO1 inhibits the expression of Notch ligand Jag1 by regulating ß-catenin (CTNNB1), which binds the Jag1 promoter. Together, our studies reveal how ROBO1/CTTNB1/JAG1 signaling in the basal compartment exerts paracrine control of Notch signaling in the luminal compartment to regulate alveolar differentiation during pregnancy.

A new SARS-CoV-2 lineage that shares mutations with known Variants of Concern is rejected by automated sequence repository quality control.

We report a SARS-CoV-2 lineage that shares N501Y, P681H, and other mutations with known variants of concern, such as B.1.1.7. This lineage, which we refer to as B.1.x (COG-UK sometimes references similar samples as B.1.324.1), is present in at least 20 states across the USA and in at least six countries. However, a large deletion causes the sequence to be automatically rejected from repositories, suggesting that the frequency of this new lineage is underestimated using public data. Recent dynamics based on 339 samples obtained in Santa Cruz County, CA, USA suggest that B.1.x may be increasing in frequency at a rate similar to that of B.1.1.7 in Southern California. At present the functional differences between this variant B.1.x and other circulating SARS-CoV-2 variants are unknown, and further studies on secondary attack rates, viral loads, immune evasion and/or disease severity are needed to determine if it poses a public health concern. Nonetheless, given what is known from well-studied circulating variants of concern, it seems unlikely that the lineage could pose larger concerns for human health than many already globally distributed lineages. Our work highlights a need for rapid turnaround time from sequence generation to submission and improved sequence quality control that removes submission bias. We identify promising paths toward this goal.

The role of exome sequencing in newborn screening for inborn errors of metabolism.

Public health newborn screening (NBS) programs provide population-scale ascertainment of rare, treatable conditions that require urgent intervention. Tandem mass spectrometry (MS/MS) is currently used to screen newborns for a panel of rare inborn errors of metabolism (IEMs)1-4. The NBSeq project evaluated whole-exome sequencing (WES) as an innovative methodology for NBS. We obtained archived residual dried blood spots and data for nearly all IEM cases from the 4.5 million infants born in California between mid-2005 and 2013 and from some infants who screened positive by MS/MS, but were unaffected upon follow-up testing. WES had an overall sensitivity of 88% and specificity of 98.4%, compared to 99.0% and 99.8%, respectively for MS/MS, although effectiveness varied among individual IEMs. Thus, WES alone was insufficiently sensitive or specific to be a primary screen for most NBS IEMs. However, as a secondary test for infants with abnormal MS/MS screens, WES could reduce false-positive results, facilitate timely case resolution and in some instances even suggest more appropriate or specific diagnosis than that initially obtained. This study represents the largest, to date, sequencing effort of an entire population of IEM-affected cases, allowing unbiased assessment of current capabilities of WES as a tool for population screening.

HNRNPA1 promotes recognition of splice site decoys by U2AF2 in vivo.

Alternative pre-mRNA splicing plays a major role in expanding the transcript output of human genes. This process is regulated, in part, by the interplay of trans-acting RNA binding proteins (RBPs) with myriad cis-regulatory elements scattered throughout pre-mRNAs. These molecular recognition events are critical for defining the protein-coding sequences (exons) within pre-mRNAs and directing spliceosome assembly on noncoding regions (introns). One of the earliest events in this process is recognition of the 3' splice site (3'ss) by U2 small nuclear RNA auxiliary factor 2 (U2AF2). Splicing regulators, such as the heterogeneous nuclear ribonucleoprotein A1 (HNRNPA1), influence spliceosome assembly both in vitro and in vivo, but their mechanisms of action remain poorly described on a global scale. HNRNPA1 also promotes proofreading of 3'ss sequences though a direct interaction with the U2AF heterodimer. To determine how HNRNPA1 regulates U2AF-RNA interactions in vivo, we analyzed U2AF2 RNA binding specificity using individual-nucleotide resolution crosslinking immunoprecipitation (iCLIP) in control and HNRNPA1 overexpression cells. We observed changes in the distribution of U2AF2 crosslinking sites relative to the 3'ss of alternative cassette exons but not constitutive exons upon HNRNPA1 overexpression. A subset of these events shows a concomitant increase of U2AF2 crosslinking at distal intronic regions, suggesting a shift of U2AF2 to "decoy" binding sites. Of the many noncanonical U2AF2 binding sites, Alu-derived RNA sequences represented one of the most abundant classes of HNRNPA1-dependent decoys. We propose that one way HNRNPA1 regulates exon definition is to modulate the interaction of U2AF2 with decoy or bona fide 3'ss.

The 3' end of the story: deciphering combinatorial interactions that control mRNA fate.

A new study investigates how microRNAs affect the binding of proteins to RNA.

Improved genome assembly of American alligator genome reveals conserved architecture of estrogen signaling.

The American alligator, Alligator mississippiensis, like all crocodilians, has temperature-dependent sex determination, in which the sex of an embryo is determined by the incubation temperature of the egg during a critical period of development. The lack of genetic differences between male and female alligators leaves open the question of how the genes responsible for sex determination and differentiation are regulated. Insight into this question comes from the fact that exposing an embryo incubated at male-producing temperature to estrogen causes it to develop ovaries. Because estrogen response elements are known to regulate genes over long distances, a contiguous genome assembly is crucial for predicting and understanding their impact. We present an improved assembly of the American alligator genome, scaffolded with in vitro proximity ligation (Chicago) data. We use this assembly to scaffold two other crocodilian genomes based on synteny. We perform RNA sequencing of tissues from American alligator embryos to find genes that are differentially expressed between embryos incubated at male- versus female-producing temperature. Finally, we use the improved contiguity of our assembly along with the current model of CTCF-mediated chromatin looping to predict regions of the genome likely to contain estrogen-responsive genes. We find that these regions are significantly enriched for genes with female-biased expression in developing gonads after the critical period during which sex is determined by incubation temperature. We thus conclude that estrogen signaling is a major driver of female-biased gene expression in the post-temperature sensitive period gonads.

Modulation of nonsense mediated decay by rapamycin.

Rapamycin is a naturally occurring macrolide whose target is at the core of nutrient and stress regulation in a wide range of species. Despite well-established roles as an inhibitor of cap-dependent mRNA translation, relatively little is known about its effects on other modes of RNA processing. Here, we characterize the landscape of rapamycin-induced post-transcriptional gene regulation. Transcriptome analysis of rapamycin-treated cells reveals genome-wide changes in alternative mRNA splicing and pronounced changes in NMD-sensitive isoforms. We demonstrate that despite well-documented attenuation of cap-dependent mRNA translation, rapamycin can augment NMD of certain transcripts. Rapamycin-treatment significantly reduces the levels of both endogenous and exogenous Premature Termination Codon (PTC)-containing mRNA isoforms and its effects are dose-, UPF1- and 4EBP-dependent. The PTC-containing SRSF6 transcript exhibits a shorter half-life upon rapamycin-treatment as compared to the non-PTC isoform. Rapamycin-treatment also causes depletion of PTC-containing mRNA isoforms from polyribosomes, underscoring the functional relationship between translation and NMD. Enhanced NMD activity also correlates with an enrichment of the nuclear Cap Binding Complex (CBC) in rapamycin-treated cells. Our data demonstrate that rapamycin modulates global RNA homeostasis by NMD.

IGF2BP3 Modulates the Interaction of Invasion-Associated Transcripts with RISC.

Insulin-like growth factor 2 mRNA binding protein 3 (IGF2BP3) expression correlates with malignancy, but its role(s) in pathogenesis remains enigmatic. We interrogated the IGF2BP3-RNA interaction network in pancreatic ductal adenocarcinoma (PDAC) cells. Using a combination of genome-wide approaches, we have identified 164 direct mRNA targets of IGF2BP3. These transcripts encode proteins enriched for functions such as cell migration, proliferation, and adhesion. Loss of IGF2BP3 reduced PDAC cell invasiveness and remodeled focal adhesion junctions. Individual nucleotide resolution crosslinking immunoprecipitation (iCLIP) revealed significant overlap of IGF2BP3 and microRNA (miRNA) binding sites. IGF2BP3 promotes association of the RNA-induced silencing complex (RISC) with specific transcripts. Our results show that IGF2BP3 influences a malignancy-associated RNA regulon by modulating miRNA-mRNA interactions.

RNA-binding protein IGF2BP3 targeting of oncogenic transcripts promotes hematopoietic progenitor proliferation.

Posttranscriptional control of gene expression is important for defining both normal and pathological cellular phenotypes. In vitro, RNA-binding proteins (RBPs) have recently been shown to play important roles in posttranscriptional regulation; however, the contribution of RBPs to cell specification is not well understood. Here, we determined that the RBP insulin-like growth factor 2 mRNA-binding protein 3 (IGF2BP3) is specifically overexpressed in mixed lineage leukemia-rearranged (MLL-rearranged) B-acute lymphoblastic leukemia (B-ALL), which constitutes a subtype of this malignancy associated with poor prognosis and high risk of relapse. IGF2BP3 was required for the survival of B-ALL cell lines, as knockdown led to decreased proliferation and increased apoptosis. Enforced expression of IGF2BP3 provided murine BM cells with a strong survival advantage, led to proliferation of hematopoietic stem and progenitor cells, and skewed hematopoietic development to the B cell/myeloid lineage. Cross-link immunoprecipitation and high throughput sequencing uncovered the IGF2BP3-regulated transcriptome, which includes oncogenes MYC and CDK6 as direct targets. IGF2BP3 regulated transcripts via targeting elements within 3' untranslated regions (3'UTR), and enforced IGF2BP3 expression in mice resulted in enhanced expression of Myc and Cdk6 in BM. Together, our data suggest that IGF2BP3-mediated targeting of oncogenic transcripts may represent a critical pathogenetic mechanism in MLL-rearranged B-ALL and support IGF2BP3 and its cognate RNA-binding partners as potential therapeutic targets in this disease.

Studying Isoform-Specific mRNA Recruitment to Polyribosomes with Frac-seq.

Gene expression profiling is widely used as a measure of the protein output of cells. However, it is becoming more evident that there are multiple layers of post-transcriptional gene regulation that greatly impact protein output (Battle et al., Science 347:664-667, 2014; Khan et al., Science 342:1100-1104, 2013; Vogel et al., Mol Syst Biol 6:400, 2010). Alternative splicing (AS) impacts the expression of protein coding genes in several ways. Firstly, AS increases exponentially the coding-capacity of genes generating multiple transcripts from the same genomic sequence. Secondly, alternatively spliced mRNAs are subjected differentially to RNA-degradation via pathways such as nonsense mediated decay (AS-NMD) or microRNAs (Shyu et al., EMBO J 27:471-481, 2008). And thirdly, cytoplasmic export from the nucleus and translation are regulated in an isoform-specific manner, adding an extra layer of regulation that impacts the protein output of the cell (Martin and Ephrussi, Cell 136:719-730, 2009; Sterne-Weiler et al., Genome Res 23:1615-1623, 2013). These data highlight the need of a method that allows analyzing both the nuclear events (AS) and the cytoplasmic fate (polyribosome-binding) of individual mRNA isoforms.In order to determine how alternative splicing determines the polyribosome association of mRNA isoforms we developed Frac-seq. Frac-seq combines subcellular fractionation and high throughput RNA sequencing (RNA-seq). Frac-seq gives a window onto the translational fate of specific alternatively spliced isoforms on a genome-wide scale. There is evidence of preferential translation of specific mRNA isoforms (Coldwell and Morley, Mol Cell Biol 26:8448-8460, 2006; Sanford et al., Genes Dev 18:755-768; Zhong et al., Mol Cell 35:1-10, 2009; Michlewski et al., Mol Cell 30:179-189, 2008); the advantage of Frac-seq is that it allows analyzing the binding of alternatively spliced isoforms to polyribosomes and comparing their relative abundance to the cytosolic fraction. Polyribosomes are resolved by sucrose gradient centrifugation of cytoplasmic extracts, subsequent reading and extraction. The total mRNA fraction is taken prior ultracentrifugation as a measure of all mRNAs present in the sample. Both populations of RNAs are then isolated using phenol-chloroform precipitation; polyadenylated RNAs are selected and converted into libraries and sequenced. Bioinformatics analysis is then performed to measure alternatively spliced isoforms; several tools can be used such as MISO, RSEM, or Cufflinks (Katz et al., Nat Methods 7:1009-1015, 2010; Li and Dewey, BMC Bioinformatics 12:323, 2011; Trapnell et al., Nat Protoc 7:562-578, 2012). Comparison of total mRNAs and polyribosome-bound mRNAs can be used as a measure of the polyribosome association of specific isoforms based on the presence/absence of specific alternative splicing events in each fraction. Frac-seq shows that not all isoforms from a gene are equally loaded into polyribosomes, that mRNA preferential loading does not always correlate to its expression in the cytoplasm and that the presence of specific events such as microRNA binding sites or Premature Termination Codons determine the loading of specific isoforms into polyribosomes.