The Role of mRNA Translation in the Regulation of Ribosome Trafficking to the Endoplasmic Reticulum

The goal of this research is to identify the trafficking patterns that direct ribosomes to the endoplasmic reticulum (ER). It is widely believed that the SRP pathway is the only mechanism that cells use to localize mRNA and ribosomes to the ER, but this has been found not to be a sufficient explanation for the patterns of RNA localization in cells, namely that non-signal sequence-containing mRNA are translated on the ER and that ribosomes retain their membrane association after translation termination. First, a summary of the history of the field is presented to provide context for the key, unanswered questions in the field. Then, experiments employing [32Pi] pulse-chase labeling of HeLa cells over a time course to follow nascent ribosome trafficking are presented. The purpose of the cell labeling was to track rRNA processing and assembly into nascent ribosomes, followed by their export into the cytoplasm and recruitment into active polysomes. A detergent-based cell fractionation procedure was also utilized to separate the cytosol and ER compartments in order to observe ribosomes on their path as they exit the nucleus and either localize to the ER or cytosolic cellular compartment. Through this method, it was seen that ribosomes appear in both compartments at the same time, suggesting a mechanism may be occurring in addition to SRP-dependent ribosome trafficking. This research provides an understanding toward a mechanism that is not currently known, but will one day more fully explain the patterns of ribosomal localization.

The DLK-1 kinase promotes mRNA stability and local translation in C. elegans synapses and axon regeneration.

Growth cone guidance and synaptic plasticity involve dynamic local changes in proteins at axons and dendrites. The Dual-Leucine zipper Kinase MAPKKK (DLK) has been previously implicated in synaptogenesis and axon outgrowth in C. elegans and other animals. Here we show that in C. elegans DLK-1 regulates not only proper synapse formation and axon morphology but also axon regeneration by influencing mRNA stability. DLK-1 kinase signals via a MAPKAP kinase, MAK-2, to stabilize the mRNA encoding CEBP-1, a bZip protein related to CCAAT/enhancer-binding proteins, via its 3'UTR. Inappropriate upregulation of cebp-1 in adult neurons disrupts synapses and axon morphology. CEBP-1 and the DLK-1 pathway are essential for axon regeneration after laser axotomy in adult neurons, and axotomy induces translation of CEBP-1 in axons. Our findings identify the DLK-1 pathway as a regulator of mRNA stability in synapse formation and maintenance and also in adult axon regeneration.

Characterization of a Full-Length TTP Family Member Association with RNA Sequence Elements

Post-transcriptional regulation of cytoplasmic mRNAs is an efficient mechanism of regulating the amounts of active protein within a eukaryotic cell. RNA sequence elements located in the untranslated regions of mRNAs can influence transcript degradation or translation through associations with RNA-binding proteins. Tristetraprolin (TTP) is the best known member of a family of CCCH zinc finger proteins that targets adenosine-uridine rich element (ARE) binding sites in the 3’ untranslated regions (UTRs) of mRNAs, promoting transcript deadenylation through the recruitment of deadenylases. More specifically, TTP has been shown to bind AREs located in the 3’-UTRs of transcripts with known roles in the inflammatory response. The mRNA-binding region of the protein is the highly conserved CCCH tandem zinc finger (TZF) domain. The synthetic TTP TZF domain has been shown to bind with high affinity to the 13-mer sequence of UUUUAUUUAUUUU. However, the binding affinities of full-length TTP family members to the same sequence and its variants are unknown. Furthermore, the distance needed between two overlapping or neighboring UUAUUUAUU 9-mers for tandem binding events of a full-length TTP family member to a target transcript has not been explored. To address these questions, we recombinantly expressed and purified the full-length C. albicans TTP family member Zfs1. Using full-length Zfs1, tagged at the N-terminus with maltose binding protein (MBP), we determined the binding affinities of the protein to the optimal TTP binding sequence, UUAUUUAUU. Fluorescence anisotropy experiments determined that the binding affinities of MBP-Zfs1 to non-canonical AREs were influenced by ionic buffer strength, suggesting that transcript selectivity may be affected by intracellular conditions. Furthermore, electrophoretic mobility shift assays (EMSAs) revealed that separation of two core AUUUA sequences by two uridines is sufficient for tandem binding of MBP-Zfs1. Finally, we found evidence for tandem binding of MBP-Zfs1 to a 27-base RNA oligonucleotide containing only a single ARE-binding site, and showed that this was concentration and RNA length dependent; this phenomenon had not been seen previously. These data suggest that the association of the TTP TZF domain and the TZF domains of other species, to ARE-binding sites is highly conserved. Domains outside of the TZF domain may mediate transcript selectivity in changing cellular conditions, and promote protein-RNA interactions not associated with the ARE-binding TZF domain.

In summary, the evidence presented here suggests that Zfs1-mediated decay of mRNA targets may require additional interactions, in addition to ARE-TZF domain associations, to promote transcript destabilization and degradation. These studies further our understanding of post-transcriptional steps in gene regulation.

RNA Localization and Translational Regulation on the Endoplasmic Reticulum

mRNA localization is emerging as a critical cellular mechanism for the spatiotemporal regulation of protein expression and serves important roles in oogenesis, embryogenesis, cell fate specification, and synapse formation. Signal sequence-encoding mRNAs are localized to the endoplasmic reticulum (ER) membrane by either of two mechanisms, a canonical mechanism of translation on ER-bound ribosomes (signal recognition particle pathway), or a poorly understood direct ER anchoring mechanism. In this study, we identify that the ER integral membrane proteins function as RNA-binding proteins and play important roles in the direct mRNA anchoring to the ER. We report that one of the ER integral membrane RNA-binding protein, AEG-1 (astrocyte elevated gene-1), functions in the direct ER anchoring and translational regulation of mRNAs encoding endomembrane transmembrane proteins. HITS-CLIP and PAR-CLIP analyses of the AEG-1 mRNA interactome of human hepatocellular carcinoma cells revealed a high enrichment for mRNAs encoding endomembrane organelle proteins, most notably encoding transmembrane proteins. AEG-1 binding sites were highly enriched in the coding sequence and displayed a signature cluster enrichment downstream of encoded transmembrane domains. In overexpression and knockdown models, AEG-1 expression markedly regulates translational efficiency and protein functions of two of its bound transcripts, MDR1 and NPC1. This study reveals a molecular mechanism for the selective localization of mRNAs to the ER and identifies a novel post-transcriptional gene regulation function for AEG-1 in membrane protein expression.

Visualizing Rare Watson-Crick-Like Tautomeric and Anionic Mismatches in DNA and RNA

The central dogma of molecular biology relies on the correct Watson-Crick (WC) geometry of canonical deoxyribonucleic acid (DNA) dG•dC and dA•dT base pairs to replicate and transcribe genetic information with speed and an astonishing level of fidelity. In addition, the Watson-Crick geometry of canonical ribonucleic acid (RNA) rG•rC and rA•rU base pairs is highly conserved to ensure that proteins are translated with high fidelity. However, numerous other potential nucleobase tautomeric and ionic configurations are possible that can give rise to entirely new pairing modes between the nucleotide bases. Very early on, James Watson and Francis Crick recognized their importance and in 1953 postulated that if bases adopted one of their less energetically disfavored tautomeric forms (and later ionic forms) during replication it could lead to the formation of a mismatch with a Watson-Crick-like geometry and could give rise to “natural mutations.”

Since this time numerous studies have provided evidence in support of this hypothesis and have expanded upon it; computational studies have addressed the energetic feasibilities of different nucleobases’ tautomeric and ionic forms in siico; crystallographic studies have trapped different mismatches with WC-like geometries in polymerase or ribosome active sites. However, no direct evidence has been given for (i) the direct existence of these WC-like mismatches in canonical DNA duplex, RNA duplexes, or non-coding RNAs; (ii) which, if any, tautomeric or ionic form stabilizes the WC-like geometry. This thesis utilizes nuclear magnetic resonance (NMR) spectroscopy and rotating frame relaxation dispersion (R1ρ RD) in combination with density functional theory (DFT), biochemical assays, and targeted chemical perturbations to show that (i) dG•dT mismatches in DNA duplexes, as well as rG•rU mismatches RNA duplexes and non-coding RNAs, transiently adopt a WC-like geometry that is stabilized by (ii) an interconnected network of rapidly interconverting rare tautomers and anionic bases. These results support Watson and Crick’s tautomer hypothesis, but additionally support subsequent hypotheses invoking anionic mismatches and ultimately tie them together. This dissertation shows that a common mismatch can adopt a Watson-Crick-like geometry globally, in both DNA and RNA, and whose geometry is stabilized by a kinetically linked network of rare tautomeric and anionic bases. The studies herein also provide compelling evidence for their involvement in spontaneous replication and translation errors.

Microbial inactivation of Pseudomonas putida and Pichia pastoris using gene silencing.

Antisense deoxyoligonucleotide (ASO) gene silencing was investigated as a potential disinfection tool for industrial and drinking water treatment application. ASOs bind with their reverse complementary mRNA transcripts thereby blocking protein translation. While ASO silencing has mainly been studied in medicine, it may be useful for modulating gene expression and inactivating microorganisms in environmental applications. In this proof of concept work, gene targets were sh ble (zeocin resistance) and todE (catechol-2,3-dioxygenase) in Pichia pastoris and npt (kanamycin resistance) in Pseudomonas putida. A maximum 0.5-fold decrease in P. pastoris cell numbers was obtained following a 120 min incubation with single-stranded DNA (ssDNA) concentrations ranging from 0.2 to 200 nM as compared to the no ssDNA control. In P. putida, a maximum 5.2-fold decrease was obtained after 90 min with 400 nM ssDNA. While the silencing efficiencies varied for the 25 targets tested, these results suggest that protein activity as well as microbial growth can be altered using ASO gene silencing-based tools. If successful, this technology has the potential to eliminate some of the environmental and health issues associated with the use of strong chemical biocides. However, prior to its dissemination, more research is needed to increase silencing efficiency and develop effective delivery methods.

SplicerAV: a tool for mining microarray expression data for changes in RNA processing.

BACKGROUND: Over the past two decades more than fifty thousand unique clinical and biological samples have been assayed using the Affymetrix HG-U133 and HG-U95 GeneChip microarray platforms. This substantial repository has been used extensively to characterize changes in gene expression between biological samples, but has not been previously mined en masse for changes in mRNA processing. We explored the possibility of using HG-U133 microarray data to identify changes in alternative mRNA processing in several available archival datasets. RESULTS: Data from these and other gene expression microarrays can now be mined for changes in transcript isoform abundance using a program described here, SplicerAV. Using in vivo and in vitro breast cancer microarray datasets, SplicerAV was able to perform both gene and isoform specific expression profiling within the same microarray dataset. Our reanalysis of Affymetrix U133 plus 2.0 data generated by in vitro over-expression of HRAS, E2F3, beta-catenin (CTNNB1), SRC, and MYC identified several hundred oncogene-induced mRNA isoform changes, one of which recognized a previously unknown mechanism of EGFR family activation. Using clinical data, SplicerAV predicted 241 isoform changes between low and high grade breast tumors; with changes enriched among genes coding for guanyl-nucleotide exchange factors, metalloprotease inhibitors, and mRNA processing factors. Isoform changes in 15 genes were associated with aggressive cancer across the three breast cancer datasets. CONCLUSIONS: Using SplicerAV, we identified several hundred previously uncharacterized isoform changes induced by in vitro oncogene over-expression and revealed a previously unknown mechanism of EGFR activation in human mammary epithelial cells. We analyzed Affymetrix GeneChip data from over 400 human breast tumors in three independent studies, making this the largest clinical dataset analyzed for en masse changes in alternative mRNA processing. The capacity to detect RNA isoform changes in archival microarray data using SplicerAV allowed us to carry out the first analysis of isoform specific mRNA changes directly associated with cancer survival.

permGPU: Using graphics processing units in RNA microarray association studies.

BACKGROUND: Many analyses of microarray association studies involve permutation, bootstrap resampling and cross-validation, that are ideally formulated as embarrassingly parallel computing problems. Given that these analyses are computationally intensive, scalable approaches that can take advantage of multi-core processor systems need to be developed. RESULTS: We have developed a CUDA based implementation, permGPU, that employs graphics processing units in microarray association studies. We illustrate the performance and applicability of permGPU within the context of permutation resampling for a number of test statistics. An extensive simulation study demonstrates a dramatic increase in performance when using permGPU on an NVIDIA GTX 280 card compared to an optimized C/C++ solution running on a conventional Linux server. CONCLUSIONS: permGPU is available as an open-source stand-alone application and as an extension package for the R statistical environment. It provides a dramatic increase in performance for permutation resampling analysis in the context of microarray association studies. The current version offers six test statistics for carrying out permutation resampling analyses for binary, quantitative and censored time-to-event traits.

Yeast screens identify the RNA polymerase II CTD and SPT5 as relevant targets of BRCA1 interaction.

BRCA1 has been implicated in numerous DNA repair pathways that maintain genome integrity, however the function responsible for its tumor suppressor activity in breast cancer remains obscure. To identify the most highly conserved of the many BRCA1 functions, we screened the evolutionarily distant eukaryote Saccharomyces cerevisiae for mutants that suppressed the G1 checkpoint arrest and lethality induced following heterologous BRCA1 expression. A genome-wide screen in the diploid deletion collection combined with a screen of ionizing radiation sensitive gene deletions identified mutants that permit growth in the presence of BRCA1. These genes delineate a metabolic mRNA pathway that temporally links transcription elongation (SPT4, SPT5, CTK1, DEF1) to nucleopore-mediated mRNA export (ASM4, MLP1, MLP2, NUP2, NUP53, NUP120, NUP133, NUP170, NUP188, POM34) and cytoplasmic mRNA decay at P-bodies (CCR4, DHH1). Strikingly, BRCA1 interacted with the phosphorylated RNA polymerase II (RNAPII) carboxy terminal domain (P-CTD), phosphorylated in the pattern specified by the CTDK-I kinase, to induce DEF1-dependent cleavage and accumulation of a RNAPII fragment containing the P-CTD. Significantly, breast cancer associated BRCT domain defects in BRCA1 that suppressed P-CTD cleavage and lethality in yeast also suppressed the physical interaction of BRCA1 with human SPT5 in breast epithelial cells, thus confirming SPT5 as a relevant target of BRCA1 interaction. Furthermore, enhanced P-CTD cleavage was observed in both yeast and human breast cells following UV-irradiation indicating a conserved eukaryotic damage response. Moreover, P-CTD cleavage in breast epithelial cells was BRCA1-dependent since damage-induced P-CTD cleavage was only observed in the mutant BRCA1 cell line HCC1937 following ectopic expression of wild type BRCA1. Finally, BRCA1, SPT5 and hyperphosphorylated RPB1 form a complex that was rapidly degraded following MMS treatment in wild type but not BRCA1 mutant breast cells. These results extend the mechanistic links between BRCA1 and transcriptional consequences in response to DNA damage and suggest an important role for RNAPII P-CTD cleavage in BRCA1-mediated cancer suppression.

MicroRNA antagonism of the picornaviral life cycle: alternative mechanisms of interference.

In addition to modulating the function and stability of cellular mRNAs, microRNAs can profoundly affect the life cycles of viruses bearing sequence complementary targets, a finding recently exploited to ameliorate toxicities of vaccines and oncolytic viruses. To elucidate the mechanisms underlying microRNA-mediated antiviral activity, we modified the 3' untranslated region (3'UTR) of Coxsackievirus A21 to incorporate targets with varying degrees of homology to endogenous microRNAs. We show that microRNAs can interrupt the picornavirus life-cycle at multiple levels, including catalytic degradation of the viral RNA genome, suppression of cap-independent mRNA translation, and interference with genome encapsidation. In addition, we have examined the extent to which endogenous microRNAs can suppress viral replication in vivo and how viruses can overcome this inhibition by microRNA saturation in mouse cancer models.

Use of 16S ribosomal RNA gene analyses to characterize the bacterial signature associated with poor oral health in West Virginia.

BACKGROUND: West Virginia has the worst oral health in the United States, but the reasons for this are unclear. This pilot study explored the etiology of this disparity using culture-independent analyses to identify bacterial species associated with oral disease. METHODS: Bacteria in subgingival plaque samples from twelve participants in two independent West Virginia dental-related studies were characterized using 16S rRNA gene sequencing and Human Oral Microbe Identification Microarray (HOMIM) analysis. Unifrac analysis was used to characterize phylogenetic differences between bacterial communities obtained from plaque of participants with low or high oral disease, which was further evaluated using clustering and Principal Coordinate Analysis. RESULTS: Statistically different bacterial signatures (P<0.001) were identified in subgingival plaque of individuals with low or high oral disease in West Virginia based on 16S rRNA gene sequencing. Low disease contained a high frequency of Veillonella and Streptococcus, with a moderate number of Capnocytophaga. High disease exhibited substantially increased bacterial diversity and included a large proportion of Clostridiales cluster bacteria (Selenomonas, Eubacterium, Dialister). Phylogenetic trees constructed using 16S rRNA gene sequencing revealed that Clostridiales were repeated colonizers in plaque associated with high oral disease, providing evidence that the oral environment is somehow influencing the bacterial signature linked to disease. CONCLUSIONS: Culture-independent analyses identified an atypical bacterial signature associated with high oral disease in West Virginians and provided evidence that the oral environment influenced this signature. Both findings provide insight into the etiology of the oral disparity in West Virginia.

Desensitization, internalization, and signaling functions of beta-arrestins demonstrated by RNA interference.

Beta-arrestins bind to activated G protein-coupled receptor kinase-phosphorylated receptors, which leads to their desensitization with respect to G proteins, internalization via clathrin-coated pits, and signaling via a growing list of "scaffolded" pathways. To facilitate the discovery of novel adaptor and signaling roles of beta-arrestins, we have developed and validated a generally applicable interfering RNA approach for selectively suppressing beta-arrestins 1 or 2 expression by up to 95%. Beta-arrestin depletion in HEK293 cells leads to enhanced cAMP generation in response to beta(2)-adrenergic receptor stimulation, markedly reduced beta(2)-adrenergic receptor and angiotensin II receptor internalization and impaired activation of the MAP kinases ERK 1 and 2 by angiotensin II. This approach should allow discovery of novel signaling and regulatory roles for the beta-arrestins in many seven-membrane-spanning receptor systems.

Analysis of the genome and transcriptome of Cryptococcus neoformans var. grubii reveals complex RNA expression and microevolution leading to virulence attenuation.

Cryptococcus neoformans is a pathogenic basidiomycetous yeast responsible for more than 600,000 deaths each year. It occurs as two serotypes (A and D) representing two varieties (i.e. grubii and neoformans, respectively). Here, we sequenced the genome and performed an RNA-Seq-based analysis of the C. neoformans var. grubii transcriptome structure. We determined the chromosomal locations, analyzed the sequence/structural features of the centromeres, and identified origins of replication. The genome was annotated based on automated and manual curation. More than 40,000 introns populating more than 99% of the expressed genes were identified. Although most of these introns are located in the coding DNA sequences (CDS), over 2,000 introns in the untranslated regions (UTRs) were also identified. Poly(A)-containing reads were employed to locate the polyadenylation sites of more than 80% of the genes. Examination of the sequences around these sites revealed a new poly(A)-site-associated motif (AUGHAH). In addition, 1,197 miscRNAs were identified. These miscRNAs can be spliced and/or polyadenylated, but do not appear to have obvious coding capacities. Finally, this genome sequence enabled a comparative analysis of strain H99 variants obtained after laboratory passage. The spectrum of mutations identified provides insights into the genetics underlying the micro-evolution of a laboratory strain, and identifies mutations involved in stress responses, mating efficiency, and virulence.

A dimensionless number for understanding the evolutionary dynamics of antigenically variable RNA viruses.

Antigenically variable RNA viruses are significant contributors to the burden of infectious disease worldwide. One reason for their ubiquity is their ability to escape herd immunity through rapid antigenic evolution and thereby to reinfect previously infected hosts. However, the ways in which these viruses evolve antigenically are highly diverse. Some have only limited diversity in the long-run, with every emergence of a new antigenic variant coupled with a replacement of the older variant. Other viruses rapidly accumulate antigenic diversity over time. Others still exhibit dynamics that can be considered evolutionary intermediates between these two extremes. Here, we present a theoretical framework that aims to understand these differences in evolutionary patterns by considering a virus's epidemiological dynamics in a given host population. Our framework, based on a dimensionless number, probabilistically anticipates patterns of viral antigenic diversification and thereby quantifies a virus's evolutionary potential. It is therefore similar in spirit to the basic reproduction number, the well-known dimensionless number which quantifies a pathogen's reproductive potential. We further outline how our theoretical framework can be applied to empirical viral systems, using influenza A/H3N2 as a case study. We end with predictions of our framework and work that remains to be done to further integrate viral evolutionary dynamics with disease ecology.

Nutritional control of mRNA isoform expression during developmental arrest and recovery in C. elegans.

Nutrient availability profoundly influences gene expression. Many animal genes encode multiple transcript isoforms, yet the effect of nutrient availability on transcript isoform expression has not been studied in genome-wide fashion. When Caenorhabditis elegans larvae hatch without food, they arrest development in the first larval stage (L1 arrest). Starved larvae can survive L1 arrest for weeks, but growth and post-embryonic development are rapidly initiated in response to feeding. We used RNA-seq to characterize the transcriptome during L1 arrest and over time after feeding. Twenty-seven percent of detectable protein-coding genes were differentially expressed during recovery from L1 arrest, with the majority of changes initiating within the first hour, demonstrating widespread, acute effects of nutrient availability on gene expression. We used two independent approaches to track expression of individual exons and mRNA isoforms, and we connected changes in expression to functional consequences by mining a variety of databases. These two approaches identified an overlapping set of genes with alternative isoform expression, and they converged on common functional patterns. Genes affecting mRNA splicing and translation are regulated by alternative isoform expression, revealing post-transcriptional consequences of nutrient availability on gene regulation. We also found that phosphorylation sites are often alternatively expressed, revealing a common mode by which alternative isoform expression modifies protein function and signal transduction. Our results detail rich changes in C. elegans gene expression as larvae initiate growth and post-embryonic development, and they provide an excellent resource for ongoing investigation of transcriptional regulation and developmental physiology.