Functions and intramolecular interactions of ribosomal RNA in decoding of termination codons

Two regions in the 3$\prime$ domain of 16S rRNA (the RNA of the small ribosomal subunit) have been implicated in decoding of termination codons. Using segment-directed PCR random mutagenesis, I isolated 33 translational suppressor mutations in the 3$\prime$ domain of 16S rRNA. Characterization of the mutations by both genetic and biochemical methods indicated that some of the mutations are defective in UGA-specific peptide chain termination and that others may be defective in peptide chain termination at all termination codons. The studies of the mutations at an internal loop in the non-conserved region of helix 44 also indicated that this structure, in a non-conserved region of 16S rRNA, is involved in both peptide chain termination and assembly of 16S rRNA.^ With a suppressible trpA UAG nonsense mutation, a spontaneously arising translational suppressor mutation was isolated in the rrnB operon cloned into a pBR322-derived plasmid. The mutation caused suppression of UAG at two codon positions in trpA but did not suppress UAA or UGA mutations at the same trpA positions. The specificity of the rRNA suppressor mutation suggests that it may cause a defect in UAG-specific peptide chain termination. The mutation is a single nucleotide deletion (G2484$\Delta$) in helix 89 of 23S rRNA (the large RNA of the large ribosomal subunit). The result indicates a functional interaction between two regions of 23S rRNA. Furthermore, it provides suggestive in vivo evidence for the involvement of the peptidyl-transferase center of 23S rRNA in peptide chain termination. The $\Delta$2484 and A1093/$\Delta$2484 (double) mutations were also observed to alter the decoding specificity of the suppressor tRNA lysT(U70), which has a mutation in its acceptor stem. That result suggests that there is an interaction between the stem-loop region of helix 89 of 23S rRNA and the acceptor stem of tRNA during decoding and that the interaction is important for the decoding specificity of tRNA.^ Using gene manipulation procedures, I have constructed a new expression vector to express and purify the cellular protein factors required for a recently developed, realistic in vitro termination assay. The gene for each protein was cloned into the newly constructed vector in such a way that expression yielded a protein with an N-terminal affinity tag, for specific, rapid purification. The amino terminus was engineered so that, after purification, the unwanted N-terminal tag can be completely removed from the protein by thrombin cleavage, yielding a natural amino acid sequence for each protein. I have cloned the genes for EF-G and all three release factors into this new expression vector and the genes for all the other protein factors into a pCAL-n expression vector. These constructs will allow our laboratory group to quickly and inexpensively purify all the protein factors needed for the new in vitro termination assay. (Abstract shortened by UMI.) ^

Involvement of a region of 23S ribosomal RNA of {\it E. coli\/} in decoding of termination codons

Involvement of E. coli 23S ribosomal RNA (rRNA) in decoding of termination codons was first indicated by the characterization of a 23S rRNA mutant that causes UGA-specific nonsense suppression. The work described here was begun to test the hypothesis that more 23S rRNA suppressors of specific nonsense mutations can be isolated and that they would occur non-randomly in the rRNA genes and be clustered in specific, functionally significant regions of rRNA.^ Approximately 2 kilobases of the gene for 23S rRNA were subjected to PCR random mutagenesis and the amplified products screened for suppression of nonsense mutations in trpA. All of the suppressor mutations obtained were located in a thirty-nucleotide part of the GTPase center, a conserved rRNA sequence and structure, and they and others made in that region by site-directed mutagenesis were shown to be UGA-specific in their suppression of termination codon mutations. These results proved the initial hypothesis and demonstrated that a group of nucleotides in this region are involved in decoding of the UGA termination codon. Further, it was shown that limitation of cellular availability or synthesis of L11, a ribosomal protein that binds to the GTPase center rRNA, resulted in suppression of termination codon mutations, suggesting the direct involvement of L11 in termination in vivo.^ Finally, in vivo analysis of certain site-specific mutations made in the GTPase center RNA demonstrated that (a) the G$\cdot$A base pair closing the hexanucleotide hairpin loop was not essential for normal termination, (b) the "U-turn" structure in the 1093 to 1098 hexaloop is critical for normal termination, (c) nucleotides A1095 and A1067, necessary for the binding to ribosomes of thiostrepton, an antibiotic that inhibits polypeptide release factor binding to ribosomes in vitro, are also necessary for normal peptide chain termination in vivo, and (d) involvement of this region of rRNA in termination is determined by some unique subset structure that includes particular nucleotides rather than merely by a general structural feature of the GTPase center.^ This work advances the understanding of peptide chain termination by demonstrating that the GTPase region of 23S rRNA participates in recognition of termination codons, through an associated ribosomal protein and specific conserved nucleotides and structural motifs in its RNA. ^

The sequence comparison index: A novel method for comparing proteins and proteomes

Historically morphological features were used as the primary means to classify organisms. However, the age of molecular genetics has allowed us to approach this field from the perspective of the organism's genetic code. Early work used highly conserved sequences, such as ribosomal RNA. The increasing number of complete genomes in the public data repositories provides the opportunity to look not only at a single gene, but at organisms' entire parts list. ^ Here the Sequence Comparison Index (SCI) and the Organism Comparison Index (OCI), algorithms and methods to compare proteins and proteomes, are presented. The complete proteomes of 104 sequenced organisms were compared. Over 280 million full Smith-Waterman alignments were performed on sequence pairs which had a reasonable expectation of being related. From these alignments a whole proteome phylogenetic tree was constructed. This method was also used to compare the small subunit (SSU) rRNA from each organism and a tree constructed from these results. The SSU rRNA tree by the SCI/OCI method looks very much like accepted SSU rRNA trees from sources such as the Ribosomal Database Project, thus validating the method. The SCI/OCI proteome tree showed a number of small but significant differences when compared to the SSU rRNA tree and proteome trees constructed by other methods. Horizontal gene transfer does not appear to affect the SCI/OCI trees until the transferred genes make up a large portion of the proteome. ^ As part of this work, the Database of Related Local Alignments (DaRLA) was created and contains over 81 million rows of sequence alignment information. DaRLA, while primarily used to build the whole proteome trees, can also be applied shared gene content analysis, gene order analysis, and creating individual protein trees. ^ Finally, the standard BLAST method for analyzing shared gene content was compared to the SCI method using 4 spirochetes. The SCI system performed flawlessly, finding all proteins from one organism against itself and finding all the ribosomal proteins between organisms. The BLAST system missed some proteins from its respective organism and failed to detect small ribosomal proteins between organisms. ^

Control of chromosomal rearrangements in {\it Escherichia coli\/}

A complete physical map of Escherichia coli K-12 strain MG1655 was constructed by digesting chromosomal DNA with the infrequently cutting restriction enzymes NotI, SfiI and XbaI and separating the fragments by pulsed field gel electrophoresis. The map was used to compare six K-12 strains of E. coli. Although several differences were noted and localized, the map of MG1655 was representative of all the K-12 strains tested. The maps were also used to analyze chromosomal rearrangements in the E. coli strain MG1655. The spontaneous and UV induced frequencies of tandem duplication formation were measured at several loci distributed around the chromosome. The spontaneous duplication frequency varied from 10$\sp{-5}$ to 10$\sp{-3}$ and increased at least ten-fold following mild UV irradiation treatment. Duplications of several regions of the chromosome, including the serA region and the metE region, were mapped using pulsed field gel electrophoresis. Duplications of serA were found to be large, ranging in size from 600 kb to 2100 kb. Several of the duplications isolated at serA were caused by ectopic recombination between IS5 elements and between IS186 elements. Duplications of the metE region, however, were almost exclusively the result of ectopic recombination between ribosomal RNA cistrons. Duplication frequencies were determined at both serA and metE in wild type and mismatch repair mutant strains (mutL, mutS, uvrD and recF). Even though all of the mismatch repair mutations increased duplication frequency of metE, the largest increases were observed in the mutL and mutS strains. Duplication frequency of serA was increased less dramatically by mutations in mismatch repair. Several duplications of metE isolated in a wild type and a mismatch repair mutant were mapped. The results showed that the same repeated sequences were used for duplication formation in the mismatch repair mutant as were used in the wild type strain. Several isolates showed evidence of multiple rearrangements indicating that mismatch repair may play a role in stabilizing the genome by controlling chromosomal rearrangement. ^

A NEW TUMOR SUPPRESSOR GENE CANDIDATE REGULATED BY THE NON-CODING RNA PCA3 IN HUMAN PROSTATE CANCER

Prostate cancer is the second leading cause of cancer-related death and the most common non-skin cancer in men in the USA. Considerable advancements in the practice of medicine have allowed a significant improvement in the diagnosis and treatment of this disease and, in recent years, both incidence and mortality rates have been slightly declining. However, it is still estimated that 1 man in 6 will be diagnosed with prostate cancer during his lifetime, and 1 man in 35 will die of the disease. In order to identify novel strategies and effective therapeutic approaches in the fight against prostate cancer, it is imperative to improve our understanding of its complex biology since many aspects of prostate cancer initiation and progression still remain elusive. The study of tumor biomarkers, due to their specific altered expression in tumor versus normal tissue, is a valid tool for elucidating key aspects of cancer biology, and may provide important insights into the molecular mechanisms underlining the tumorigenesis process of prostate cancer. PCA3, is considered the most specific prostate cancer biomarker, however its biological role, until now, remained unknown. PCA3 is a long non-coding RNA (ncRNA) expressed from chromosome 9q21 and its study led us to the discovery of a novel human gene, PC-TSGC, transcribed from the opposite strand and in an antisense orientation to PCA3. With the work presented in this thesis, we demonstrate that PCA3 exerts a negative regulatory role over PC-TSGC, and we propose PC-TSGC to be a new tumor suppressor gene that contrasts the transformation of prostate cells by inhibiting Rho-GTPases signaling pathways. Our findings provide a biological role for PCA3 in prostate cancer and suggest a new mechanism of tumor suppressor gene inactivation mediated by non-coding RNA. Also, the characterization of PCA3 and PC-TSGC led us to propose a new molecular pathway involving both genes in the transformation process of the prostate, thus providing a new piece of the jigsaw puzzle representing the complex biology of prostate cancer.

Regulation of mucin gene expression by CREB via a nonclassical retinoic acid signaling pathway.

Vitamin A and its metabolite retinoic acid (RA) are essential elements for normal lung development and the differentiation of lung epithelial cells. We previously showed that RA rapidly activated cyclic AMP response element-binding protein (CREB) in a nonclassical manner in normal human tracheobronchial epithelial (NHTBE) cells. In the present study, we further demonstrated that this nonclassical signaling of RA on the activation of CREB plays a critical role in regulating the expression of airway epithelial cell differentiation markers, the MUC2, MUC5AC, and MUC5B genes. We found that RA rapidly activates the protein kinase Calpha isozyme and transmits the activation signal to CREB via the Raf/MEK/extracellular signal-regulated kinase/p90 ribosomal S6 kinase (RSK) pathway. Activated RSK translocated from the cytoplasm to the nucleus, where it phosphorylates CREB. Activated CREB then binds to a cis-acting replication element motif on the promoter (at nucleotides [nt] -878 to -871) of the MUC5AC gene. The depletion of CREB using small interfering RNA abolished not only the RA-induced MUC5AC but also RA-induced MUC2 and MUC5B. Taken together, our findings demonstrate that CREB activation via this nonclassical RA signaling pathway may play an important role in regulating the expression of mucin genes and mediating the early biological effects of RA during normal mucous differentiation in NHTBE cells.

Clinical and microbiological aspects of linezolid resistance mediated by the cfr gene encoding a 23S rRNA methyltransferase.

The cfr (chloramphenicol-florfenicol resistance) gene encodes a 23S rRNA methyltransferase that confers resistance to linezolid. Detection of linezolid resistance was evaluated in the first cfr-carrying human hospital isolate of linezolid and methicillin-resistant Staphylococcus aureus (designated MRSA CM-05) by dilution and diffusion methods (including Etest). The presence of cfr was investigated in isolates of staphylococci colonizing the patient's household contacts and clinical isolates recovered from patients in the same unit where MRSA CM-05 was isolated. Additionally, 68 chloramphenicol-resistant Colombian MRSA isolates recovered from hospitals between 2001 and 2004 were screened for the presence of the cfr gene. In addition to erm(B), the erm(A) gene was also detected in CM-05. The isolate belonged to sequence type 5 and carried staphylococcal chromosomal cassette mec type I. We were unable to detect the cfr gene in any of the human staphylococci screened (either clinical or colonizing isolates). Agar and broth dilution methods detected linezolid resistance in CM-05. However, the Etest and disk diffusion methods failed to detect resistance after 24 h of incubation. Oxazolidinone resistance mediated by the cfr gene is rare, and acquisition by a human isolate appears to be a recent event in Colombia. The detection of cfr-mediated linezolid resistance might be compromised by the use of the disk diffusion or Etest method.

Evolution of sensory complexity recorded in a myxobacterial genome.

Myxobacteria are single-celled, but social, eubacterial predators. Upon starvation they build multicellular fruiting bodies using a developmental program that progressively changes the pattern of cell movement and the repertoire of genes expressed. Development terminates with spore differentiation and is coordinated by both diffusible and cell-bound signals. The growth and development of Myxococcus xanthus is regulated by the integration of multiple signals from outside the cells with physiological signals from within. A collection of M. xanthus cells behaves, in many respects, like a multicellular organism. For these reasons M. xanthus offers unparalleled access to a regulatory network that controls development and that organizes cell movement on surfaces. The genome of M. xanthus is large (9.14 Mb), considerably larger than the other sequenced delta-proteobacteria. We suggest that gene duplication and divergence were major contributors to genomic expansion from its progenitor. More than 1,500 duplications specific to the myxobacterial lineage were identified, representing >15% of the total genes. Genes were not duplicated at random; rather, genes for cell-cell signaling, small molecule sensing, and integrative transcription control were amplified selectively. Families of genes encoding the production of secondary metabolites are overrepresented in the genome but may have been received by horizontal gene transfer and are likely to be important for predation.

Coordinated changes in mRNA turnover, translation, and RNA processing bodies in bronchial epithelial cells following inflammatory stimulation.

Bronchial epithelial cells play a pivotal role in airway inflammation, but little is known about posttranscriptional regulation of mediator gene expression during the inflammatory response in these cells. Here, we show that activation of human bronchial epithelial BEAS-2B cells by proinflammatory cytokines interleukin-4 (IL-4) and tumor necrosis factor alpha (TNF-alpha) leads to an increase in the mRNA stability of the key chemokines monocyte chemotactic protein 1 and IL-8, an elevation of the global translation rate, an increase in the levels of several proteins critical for translation, and a reduction of microRNA-mediated translational repression. Moreover, using the BEAS-2B cell system and a mouse model, we found that RNA processing bodies (P bodies), cytoplasmic domains linked to storage and/or degradation of translationally silenced mRNAs, are significantly reduced in activated bronchial epithelial cells, suggesting a physiological role for P bodies in airway inflammation. Our study reveals an orchestrated change among posttranscriptional mechanisms, which help sustain high levels of inflammatory mediator production in bronchial epithelium during the pathogenesis of inflammatory airway diseases.

Acquisition of a natural resistance gene renders a clinical strain of methicillin-resistant Staphylococcus aureus resistant to the synthetic antibiotic linezolid.

Linezolid, which targets the ribosome, is a new synthetic antibiotic that is used for treatment of infections caused by Gram-positive pathogens. Clinical resistance to linezolid, so far, has been developing only slowly and has involved exclusively target site mutations. We have discovered that linezolid resistance in a methicillin-resistant Staphylococcus aureus hospital strain from Colombia is determined by the presence of the cfr gene whose product, Cfr methyltransferase, modifies adenosine at position 2503 in 23S rRNA in the large ribosomal subunit. The molecular model of the linezolid-ribosome complex reveals localization of A2503 within the drug binding site. The natural function of cfr likely involves protection against natural antibiotics whose site of action overlaps that of linezolid. In the chromosome of the clinical strain, cfr is linked to ermB, a gene responsible for dimethylation of A2058 in 23S rRNA. Coexpression of these two genes confers resistance to all the clinically relevant antibiotics that target the large ribosomal subunit. The association of the ermB/cfr operon with transposon and plasmid genetic elements indicates its possible mobile nature. This is the first example of clinical resistance to the synthetic drug linezolid which involves a natural resistance gene with the capability of disseminating among Gram-positive pathogenic strains.

Accessing novel developmental mechanisms in the mouse by gene trapping

Genetic analysis is a powerful method for analyzing the function of specific genes in development. I sought to identify novel genes in the mouse using a genetic analysis relying on the expression pattern and phenotype of mutated genes. To this end, I have conducted a gene trap screen using the vector $\rm SA\beta geo,$ a promoterless DNA construct that encodes a fusion protein with lacZ and neomycin resistance activities. Productive integration and expression of the $\beta$geo protein in embryonic stem (ES) cells requires integration into an active transcription unit. The endogenous regulatory elements direct reporter gene expression which reflects the expression of the endogenous gene. Of eight mouse lines generated from gene trap ES cell clones, four showed differential regulation of $\beta$geo activity during embryogenesis. These four were analyzed in more detail.^ Three of the lines RNA 1, RNA2 and RNA 3 had similar expression patterns, within subsets of cells in sites of embryonic hematopoiesis. Cloning of the trapped genes revealed that all three integrations had occurred within 45S rRNA precursor transcription units. These results imply that there exists in these cells some mechanism responsible for the efficient production of the $\beta$geo protein from an RNA polymerase I transcript that is not present in most of the cells in the embryo.^ The fourth line, GT-2, showed widespread, dynamic expression. Many of the sites of expression were important classic embryonic induction systems. Cloning of the sequences fused to the $5\sp\prime$ end of the $\beta$geo sequence revealed that the trapped gene contained significant sequence homology with a previously identified human sequence HumORF5. An open reading frame of this sequence is homologous to a group of eukaryotic proteins that are members of the RNA helicase superfamily I.^ Analysis of the gene trap lines suggests that potentially novel developmental mechanisms have been uncovered. In the case of RNA 1, 2 and 3, the differential production of ribosomal RNAs may be required for differentiation or function of the $\beta$geo positive hematopoietic cells. In the GT-2 line, a previously unsuspected temporal and spatial regulation of a putative RNA helicase implies a role for this activity during specific aspects of mouse development. ^

The role of specific regions of ribosomal RNAs in translation termination

The ribosome is a molecular machine that produces proteins in a cell. It consists of RNAs (rRNAs) and proteins. The rRNAs have been implicated in various aspects of protein biosynthesis supporting the idea that they function directly in translation. In this study the direct involvement of rRNA in translation termination was hypothesized and both genetic and biochemical strategies were designed to test this hypothesis. As a result, several regions of rRNAs from both ribosomal subunits were implicated in termination. More specifically, the mutation G1093A in an RNA of the large subunit (23S rRNA) and the mutation C1054A in the small subunit RNA (16S rRNA) of the Escherichia coli ribosome, were shown to affect the binding of the proteins that drive termination, RF1 and RF2. These mutations also caused defects in catalysis of peptidyl-tRNA hydrolysis, the last step of termination. Furthermore, the mutations affected the function of RF2 to a greater extent than that of RF1, a striking result considering the similarity of the RFs. The major defect in RF2 function was consistent with in vivo characteristics of the mutants and can be explained by the inability of the mutant rRNA sites to activate the hydrolytic center, that is the catalytic site for peptidyl-tRNA hydrolysis. Consistent with this explanation is the possibility of a direct interaction between the G1093-region (domain II of 23S rRNA) and the hydrolytic center (most likely domains IV–VI of 23S rRNA). To test that interaction hypothesis selections were performed for mutations in domains IV–VI that compensated for the growth defects caused by G1093A. Several compensatory mutations were isolated which not only restored growth in the presence of G1093A but also appeared to compensate for the termination defects caused by the G1093A. Therefore these results provided genetic evidence for an intramolecular interaction that might lead to peptidyl-tRNA hydrolysis. Finally, a new approach to the study of rRNA involvement in termination was designed. By screening a library of rRNA fragments, a fragment of the 23S rRNA (nt 74-136) was identified that caused readthrough of UGA. The antisense RNA fragment produced a similar effect. The data implicated the corresponding segment of intact 23S rRNA in termination. ^

Post-transcriptional Regulation of Mammalian Gene Expression in Non-coding Region of Target RNA

Tumor Suppressor Candidate 2 (TUSC2) is a novel tumor suppressor gene located in the human chromosome 3p21.3 region. TUSC2 mRNA transcripts could be detected on Northern blots in both normal lung and some lung cancer cell lines, but no endogenous TUSC2 protein could be detected in a majority of lung cancer cell lines. Mechanisms regulating TUSC2 protein expression and its inactivation in primary lung cancer cells are largely unknown. We investigated the role of the 5’- and 3’-untranslated regions (UTRs) of the TUSC2 gene in the regulation of TUSC2 protein expression. We found that two small upstream open-reading frames (uORFs) in the 5’UTR of TUSC2 could markedly inhibit the translational initiation of TUSC2 protein by interfering with the “scanning” of the ribosome initiation complexes. Site-specific stem-loop array reverse transcription-polymerase chain reaction (SLA-RT-PCR) verified several micoRNAs (miRNAs) targeted at 3’UTR and directed TUSC2 cleavage and degradation. In addition, we used the established let-7-targeted high mobility group A2 (Hmga2) mRNA as a model system to study the mechanism of regulation of target mRNA by miRNAs in mammalian cells under physiological conditions. There have been no evidence of direct link between mRNA downregulation and mRNA cleavages mediated by miRNAs. Here we showed that the endonucleolytic cleavages on mRNAs were initiated by mammalian miRNA in seed pairing style. Let-7 directed cleavage activities among the eight predicted potential target sites have varied efficiency, which are influenced by the positional and the structural contexts in the UTR. The 5’ cleaved RNA fragments were mostly oligouridylated at their 3’-termini and accumulated for delayed 5’–3’ degradation. RNA fragment oligouridylation played important roles in marking RNA fragments for delayed bulk degradation and in converting RNA degradation mode from 3’–5’ to 5’–3’ with cooperative efforts from both endonucleolytic and non-catalytic miRNA-induced silencing complex (miRISC). Our findings point to a mammalian miRNA-mediated mechanism for the regulation of mRNA that miRNA can decrease target mRNA through target mRNA cleavage and uridine addition

RNA-SEQUENCING APPLICATIONS: GENE EXPRESSION QUANTIFICATION AND METHYLATOR PHENOTYPE IDENTIFICATION

My dissertation focuses on two aspects of RNA sequencing technology. The first is the methodology for modeling the overdispersion inherent in RNA-seq data for differential expression analysis. This aspect is addressed in three sections. The second aspect is the application of RNA-seq data to identify the CpG island methylator phenotype (CIMP) by integrating datasets of mRNA expression level and DNA methylation status. Section 1: The cost of DNA sequencing has reduced dramatically in the past decade. Consequently, genomic research increasingly depends on sequencing technology. However it remains elusive how the sequencing capacity influences the accuracy of mRNA expression measurement. We observe that accuracy improves along with the increasing sequencing depth. To model the overdispersion, we use the beta-binomial distribution with a new parameter indicating the dependency between overdispersion and sequencing depth. Our modified beta-binomial model performs better than the binomial or the pure beta-binomial model with a lower false discovery rate. Section 2: Although a number of methods have been proposed in order to accurately analyze differential RNA expression on the gene level, modeling on the base pair level is required. Here, we find that the overdispersion rate decreases as the sequencing depth increases on the base pair level. Also, we propose four models and compare them with each other. As expected, our beta binomial model with a dynamic overdispersion rate is shown to be superior. Section 3: We investigate biases in RNA-seq by exploring the measurement of the external control, spike-in RNA. This study is based on two datasets with spike-in controls obtained from a recent study. We observe an undiscovered bias in the measurement of the spike-in transcripts that arises from the influence of the sample transcripts in RNA-seq. Also, we find that this influence is related to the local sequence of the random hexamer that is used in priming. We suggest a model of the inequality between samples and to correct this type of bias. Section 4: The expression of a gene can be turned off when its promoter is highly methylated. Several studies have reported that a clear threshold effect exists in gene silencing that is mediated by DNA methylation. It is reasonable to assume the thresholds are specific for each gene. It is also intriguing to investigate genes that are largely controlled by DNA methylation. These genes are called “L-shaped” genes. We develop a method to determine the DNA methylation threshold and identify a new CIMP of BRCA. In conclusion, we provide a detailed understanding of the relationship between the overdispersion rate and sequencing depth. And we reveal a new bias in RNA-seq and provide a detailed understanding of the relationship between this new bias and the local sequence. Also we develop a powerful method to dichotomize methylation status and consequently we identify a new CIMP of breast cancer with a distinct classification of molecular characteristics and clinical features.