The Role of Splicing Factors and Small Nuclear RNAS in Spliceosomal Formation

Protein coding genes are comprised of protein-coding exons and non-protein-coding introns. The process of splicing involves removal of the introns and joining of the exons to form a mature messenger RNA, which subsequently undergoes translation into polypeptide. The spliceosome is a large, RNA/protein assembly of five small nuclear RNAs as well as over 300 proteins, which catalyzes intron removal and exon ligation. The selection of specific exons for inclusion in the mature messenger RNA is spatio-temporally regulated and results in production of an enormous diversity of polypeptides from a single gene locus. This phenomenon, known as alternative splicing, is regulated, in part, by protein splicing factors, which target the spliceosome to exon/intron boundaries. The first part of my dissertation (Chapters II and III) focuses on the discovery and characterization of the 45 kilodalton FK506 binding protein (FKBP45), which I discovered in the silk moth, Bombyx mori, as a U1 small nuclear RNA binding protein. This protein family binds the immunosuppressants FK506 and rapamycin and contains peptidyl-prolyl cis-trans isomerase activity, which converts polypeptides from cis to trans about a proline residue. This is the first time that an FKBP has been identified in the spliceosome. The second section of my dissertation (Chapters IV, V, VI and VII) is an investigation of the potential role of small nuclear RNA sequence variants in the control of splicing. I identified 46 copies of small nuclear RNAs in the 6X whole genome shotgun of the Bombyx mori p50T strain. These variants may play a role in differential binding of specific proteins that mediate alternative splicing. Along these lines, further investigation of U2 snRNA sequence variants in Bombyx mori demonstrated that some U2 snRNAs preferentially assemble into high molecular weight spliceosomal complexes over others. Expression of snRNA variants may represent another mechanism by which the cell is able to fine tune the splicing process.

Mechanismen kleiner regulatorischer RNAs in Dictyostelium discoideum: Charakterisierung und Anwendung

Die RNA-Interferenz (RNAi) ist ein in Eukaryoten weit verbreiteter Mechanismus, der die transkriptionelle oder posttranskriptionelle Stilllegung von Genen beschreibt. Die Spezifität wird dabei durch die Sequenz einer kleinen, nicht-kodierenden RNA gewährleistet. Diese RNA leitet einen Effektorkomplex, dessen Zentrum immer von einem Argonautenprotein gebildet wird, üblicherweise zu einer komplementären mRNA. In der Folge kommt es zum Abbau des Transkripts oder zur Inhibierung der Translation. Aktuelle Veröffentlichungen konnten zudem das Aktivitätsprofil der Argonautenproteine beträchtlich erweitern: Im Zellkern vorkommende Argonautenproteine wurden beispielsweise mit Spleißvorgängen, der Promotorkontrolle von Genen und der DNA-Reparatur in Verbindung gebracht. In den letzten Jahren konnten weitreichende Kenntnisse bezüglich der Kontrolle einiger transposabler Elemente durch RNAi sowie der Biogenese kleiner regulatorischer RNAs in Dictyostelium discoideum und anderen Organismen gewonnen werden. Ein Fokus dieser Arbeit lag zunächst auf der Charakterisierung des Argonautenproteins AgnB und der Identifikation von Interaktionspartnern. Es konnte gezeigt werden, dass AgnB zumindest partiell im Zellkern der Amöbe lokalisiert und dort vermutlich regulatorische Aufgaben wahrnimmt. Gestützt wurde diese Annahme durch die massenspektrometrische und Western Blot basierte Detektion nukleärer Bindungspartner. Weiterführende Analysen konnten AgnB zudem als positiven Regulator für drei entwicklungsregulierte Gene beschreiben und so die Verbindung zum Prozess der RNA activation in der Amöbe herstellen. Identifizierte posttranslationale Modifikationen könnten diesbezüglich die Aktivität des Argonauten steuern. Mit Hilfe von Crosslink-RNA-Immunopräzipitation und anschließender Hochdurchsatz-sequenzierung oder der Northern Blot basierten Auswertung konnte eine Assoziation von AgnB und der Class I RNAs gezeigt werden. Diese Klasse umfasst nicht-kodierende RNAs mit einer Länge von etwa 42 bis 65 Nukleotiden und wurde bisher nicht als Substrat für die RNAi-Maschinerie in D. discoideum angesehen. Ein weiterer Teil dieser Arbeit beschäftigte sich mit dem Einfluss von AgnA und AgnB auf die Promotorbereiche des D. discoideum Retrotransposons DIRS-1. Im Verlauf der Untersuchungen konnte beobachtet werden, dass die Anordnung entgegengesetzt operierender DIRS-1 Promotor-sequenzen für die Stilllegung eines Reportergens ausreichend war. Darauf aufbauend konnte ein DIRS-1 basiertes knockdown System etabliert werden. Mit Hilfe dieses Systems konnten die Proteinmengen ausgewählter Zielgene so effektiv reduziert werden, dass die entsprechenden Stämme den Phänotyp des korrespondierenden knockout Stammes zeigten. Darüber hinaus war es möglich die Proteinreduktion zu modulieren, um so beispielsweise dosisabhängige Effekte zu registrieren. Vorangegangene Arbeiten zur Biogenese von micro (mi)RNAs konnten das RNA-bindende Protein RbdB als eine Hauptkomponente für die Prozessierung maturer miRNAs identifizieren. Der miRNA defiziente RbdB- Stamm wurde in dieser Arbeit zur Identifikation putativer miRNA Ziele genutzt. Dazu wurde sowohl das Transkriptom des D. discoideum Wildtyps als auch des rbdB knockout Stammes in hohem Durchsatz sequenziert, um so differentiell exprimierte Gene zu detektieren. Vielversprechende Kandidaten wurden mittels qRT-PCR verifiziert. Dabei wurde unter anderem ein putativer Transkriptionsfaktor als miRNA target identifiziert, der eine Vielzahl weiterer Gene regulieren könnte. Abschließend konnte in dieser Arbeit gezeigt werden, dass RbdB ebenfalls für die Generierung von small interfering (si)RNAs aus strukturierten Loci von Bedeutung ist. Dies weist auf mindestens zwei unterschiedliche Mechanismen zur siRNA Prozessierung in D. discoideum hin.

Investigation of the role of long non-coding RNAs in oncogene induced cellular senescence

Cellular senescence is a stable arrest of cell proliferation induced by several factors such as activated oncogenes, oxidative stress and shortening of telomeres. Senescence acts as a tumour suppression mechanism to halt the progression of cancer. However, senescence may also impact negatively upon tissue regeneration, thus contributing to the effects of ageing. The eukaryotic genome is controlled by various modes of transcriptional and translational regulation. Focus has therefore centred on the role of long non- coding RNAs (lncRNAs) in regulating the genome. Accordingly, understanding how lncRNAs function to regulate the senescent genome is integral to improving our knowledge and understanding of tumour suppression and ageing. Within this study, I set out to investigate the expression of lncRNAs’ expression within models of senescence. Through a custom expression array, I have shown that expression of multiple different lncRNAs is up-regulated and down regulated in IMR90 replicative senescent fibroblasts and oncogene-induced senescent melanocytes. LncRNA expression was determined to be specific to stable senescence-associated cell arrest and predominantly within the nucleus of senescent cells. In order to examine the function of lncRNA expression in senescence, I selected lncRNA transcript ENST0000430998 (lncRNA_98) to focus my investigations upon. LncRNA_98 was robustly upregulated within multiple models of senescence and efficiently depleted using anti-sense oligonucleotide technology. Characterisation and unbiased RNA-sequencing of lncRNA_98 deficient senescent cells highlighted a list of genes that are regulated by lncRNA_98 expression in senescent cells and may regulate aspects of the senescence program. Specifically, the formation of SAHF was impeded upon depletion of lncRNA_98 expression and levels of total pRB protein expression severely decreased. Validation and recapitulation of consequences of pRB depletion was confirmed through lncRNA_98 knock-out cells generated using CRISPR technology. Surprisingly, inhibition of ATM kinase functions permitted the restoration of pRB protein levels within lncRNA_98 deficient cells. I propose that lncRNA_98 antagonizes the ability of ATM kinase to downregulate pRB expression at a post-transcriptional level, thereby potentiating senescence. Furthermore, lncRNA expression was detected within fibroblasts of old individuals and visualised within senescent melanocytes in human benign nevi, a barrier to melanoma progression. Conversely, mining of 337 TCGA primary melanoma data sets highlighted that the lncRNA_98 gene and its expression was lost from a significant proportion of melanoma samples, consistent with lncRNA_98 having a tumour suppressor functions. The data presented in this study illustrates that lncRNA_98 expression has a regulatory role over pRB expression in senescence and may regulate aspects of tumourigenesis and ageing.

PAMP-triggered immunity against Pseudomonas syringae involves microRNA-mediated regulation of several uncharacterized R genes

Two main types of noncoding small RNA molecules have been found in plants: microRNAs (miRNAs) and small interfering RNAs (siRNAs). They differ in their biogenesis and mode of action, but share similar sizes (20-24 nt). Their precursors are processed by Dicer-Like RNase III (dcl) proteins present in Arabidopsis thaliana, and in their mature form can act as negative regulators of gene expression, being involved in a vast array of plant processes, including plant development, genomic integrity or response to stress. Small-RNA mediated regulation can occurs at transcriptional level (TGS) or at post-transcriptional level (PTGS). In recent years, the role of gene silencing in the regulation of expression of genes related to plant defence responses against bacterial pathogens is becoming clearer. Comparisons carried out in our lab between the expression profiles of different mutants affected in gene silencing, and plants challenged with Pseudomonas syringae pathovar tomato DC3000, led us to identify a set of uncharacterized R genes, belonging to the TIR-NBS-LRR gene family, differentially expressed in these conditions. Through the use of bioinformatics tools, we found a miRNA* of 22 nt putatively responsible for down-regulating expression of these R genes through the generation of siRNAs. We have also found that the corresponding pri-miRNA is down-regulated after PAMP-perception in a SA-dependent manner. We also demonstrate that plants with altered levels of miRNA* (knockdown lines or overexpression lines) exhibit altered PTI-associated phenotypes, suggesting a role for this miRNA* in this defence response against bacteria. In addition we identify one of the target genes as a negative regulator of defence response against Pseudomonas syringae.

Synthetic regulation of eukaryotic gene expression by noncoding RNA

Synthetic biological systems promise to combine the spectacular diversity of biological functionality with engineering principles to design new life to address many pressing needs. As these engineered systems advance in sophistication, there is ever-greater need for customizable, situation-specific expression of desired genes. However, existing gene control platforms are generally not modular, or do not display performance requirements required for robust phenotypic responses to input signals. This work expands the capabilities of eukaryotic gene control in two important directions.

For development of greater modularity, we extend the use of synthetic self-cleaving ribozyme switches to detect changes in input protein levels and convey that information into programmed gene expression in eukaryotic cells. We demonstrate both up- and down-regulation of levels of an output transgene by more than 4-fold in response to rising input protein levels, with maximal output gene expression approaching the highest levels observed in yeast. In vitro experiments demonstrate protein-dependent ribozyme activity modulation. We further demonstrate the platform in mammalian cells. Our switch devices do not depend on special input protein activity, and can be tailored to respond to any input protein to which a suitable RNA aptamer can be developed. This platform can potentially be employed to regulate the expression of any transgene or any endogenous gene by 3’ UTR replacement, allowing for more complex cell state-specific reprogramming.

We also address an important concern with ribozyme switches, and riboswitch performance in general, their dynamic range. While riboswitches have generally allowed for versatile and modular regulation, so far their dynamic ranges of output gene modulation have been modest, generally at most 10-fold. We address this shortcoming by developing a modular genetic amplifier for near-digital control of eukaryotic gene expression. We combine ribozyme switch-mediated regulation of a synthetic TF with TF-mediated regulation of an output gene. The amplifier platform allows for as much as 20-fold regulation of output gene expression in response to input signal, with maximal expression approaching the highest levels observed in yeast, yet being tunable to intermediate and lower expression levels. EC50 values are more than 4 times lower than in previously best-performing non-amplifier ribozyme switches. The system design retains the modular-input architecture of the ribozyme switch platform, and the near-digital dynamic ranges of TF-based gene control.

Together, these developments suggest great potential for the wide applicability of these platforms for better-performing eukaryotic gene regulation, and more sophisticated, customizable reprogramming of cellular activity.

Ghrelin gene-related peptides : multifunctional endocrine/autocrine modulators in health and disease

Ghrelin is a multi-functional peptide hormone which affects various processes including growth hormone and insulin release, appetite regulation, gut motility, metabolism and cancer cell proliferation. Ghrelin is produced in the stomach and in other normal and pathological cell types. It may act as an endocrine or autocrine/paracrine factor. The ghrelin gene encodes a precursor protein, preproghrelin, from which ghrelin and other potentially active peptides are derived by alternative mRNA splicing and/or proteolytic processing. The metabolic role of the peptide obestatin, derived from the preproghrelin C-terminal region, is controversial. However, it has direct effects on cancer cell proliferation. The regulation of ghrelin expression and the mechanisms through which the peptide products arise are unclear. We have recently re-examined the organisation of the ghrelin gene and identified several novel exons and transcripts. One transcript, which lacks the ghrelin-coding region of preproghrelin, contains the coding sequence of obestatin. Furthermore, we have identified an overlapping gene on the antisense strand of ghrelin, GHRLOS, which generates transcripts that may function as non-coding regulatory RNAs or code for novel, short bioactive peptides. The identification of these novel ghrelin-gene related transcripts and peptides raises critical questions regarding their physiological function and their role in obesity, diabetes and cancer.

Comparison of Nucleotide Sequences between Northern and Southern Philippine Isolates of Rice Grassy Stunt Virus Indicates Occurrence of Natural Genetic Reassortment

Rice grassy stunt virus is a member of the genus Tenuivirus, is persistently transmitted by a brown planthopper, and has occurred in rice plants in South, Southeast, and East Asia (similar to North and South America). We determined the complete nucleotide (nt) sequences of RNAs 1 (9760 nt), 2 (4069 nt), 3 (3127 nt), 4 (2909 nt), 5 (2704 nt), and 6 (2590 nt) of a southern Philippine isolate from South Cotabato and compared them with those of a northern Philippine isolate from Laguna (Toriyama et al., 1997, 1998). The numbers of nucleotides in the terminal untranslated regions and open reading frames were identical between the two isolates except for the 5′ untranslated region of the complementary strand of RNA 4. Overall nucleotide differences between the two isolates were only 0.08% in RNA 1, 0.58% in RNA 4, and 0.26% in RNA 5, whereas they were 2.19% in RNA 2, 8.38% in RNA 3, and 3.63% in RNA 6. In the intergenic regions, the two isolates differed by 9.12% in RNA 2, 11.6% in RNA 3, and 6.86% in RNA 6 with multiple consecutive nucleotide deletion/insertions, whereas they differed by only 0.78% in RNA 4 and 0.34% in RNA 5. The nucleotide variation in the intergenic region of RNA 6 within the South Cotabato isolate was only 0.33%. These differences in accumulation of mutations among individual RNA segments indicate that there was genetic reassortment in the two geographical isolates; RNAs 1, 4, and 5 of the two isolates came from a common ancestor, whereas RNAs 2, 3, and 6 were from two different ancestors.

Ghrelin axis genes, peptides and receptors : recent findings and future challenges

The ghrelin axis consists of the gene products of the ghrelin gene (GHRL), and their receptors, including the classical ghrelin receptor GHSR. While it is well-known that the ghrelin gene encodes the 28 amino acid ghrelin peptide hormone, it is now also clear that the locus encodes a range of other bioactive molecules, including novel peptides and non-coding RNAs. For many of these molecules, the physiological functions and cognate receptor(s) remain to be determined. Emerging research techniques, including proteogenomics, are likely to reveal further ghrelin axis-derived molecules. Studies of the role of ghrelin axis genes, peptides and receptors, therefore, promises to be a fruitful area of basic and clinical research in years to come.

Multiple evolutionary rate classes in animal genome evolution

The proportion of functional sequence in the human genome is currently a subject of debate. The most widely accepted figure is that approximately 5% is under purifying selection. In Drosophila, estimates are an order of magnitude higher, though this corresponds to a similar quantity of sequence. These estimates depend on the difference between the distribution of genomewide evolutionary rates and that observed in a subset of sequences presumed to be neutrally evolving. Motivated by the widening gap between these estimates and experimental evidence of genome function, especially in mammals, we developed a sensitive technique for evaluating such distributions and found that they are much more complex than previously apparent. We found strong evidence for at least nine well-resolved evolutionary rate classes in an alignment of four Drosophila species and at least seven classes in an alignment of four mammals, including human. We also identified at least three rate classes in human ancestral repeats. By positing that the largest of these ancestral repeat classes is neutrally evolving, we estimate that the proportion of nonneutrally evolving sequence is 30% of human ancestral repeats and 45% of the aligned portion of the genome. However, we also question whether any of the classes represent neutrally evolving sequences and argue that a plausible alternative is that they reflect variable structure-function constraints operating throughout the genomes of complex organisms.

Mutation analysis of the CDKN2A promoter in Australian melanoma families

Approximately 50% of all melanoma families worldwide show linkage to 9p21-22, but only about half of these have been shown to contain germ line CDKN2A mutations. It has been hypothesized that a proportion of these families carry mutations in the noncoding regions of CDKN2A. Several Canadian families have been reported to carry a mutation in the 5' UTR, at position -34 relative to the start site, which gives rise to a novel AUG translation initiation codon that markedly decreases translation from the wild-type AUG (Liu et al., 1999). Haplotype sharing in these Canadian families suggested that this mutation is of British origin. We sequenced 1,327 base pairs (bp) of CDKN2A, making up 1,116 bp of the 5' UTR and promoter, all of exon 1, and 61 bp of intron 1, in at least one melanoma case from 110 Australian families with three or more affected members known not to carry mutations within the p16 coding region. In addition, 431 bp upstream of the start codon was sequenced in an additional 253 affected probands from two-case melanoma families for which the CDKN2A mutation status was unknown. Several known polymorphisms at positions -33, -191, -493, and -735 were detected, in addition to four novel variants at positions 120, -252, -347, and -981 relative to the start codon. One of the probands from a two-case family was found to have the previously reported Q50R mutation. No family member was found to carry the mutation at position -34 or any other disease-associated mutation. For further investigation of noncoding CDKN2A mutations that may affect transcription, allele-specific expression analysis was carried out in 31 of the families with at least three affected members who showed either complete or "indeterminate" 9p haplotype sharing without CDKN2A exonic mutations. Reverse transcription polymerase chain reaction and automated sequencing showed expression of both CDKN2A alleles in all family members tested. The lack of CDKN2A promoter mutations and the absence of transcriptional silencing in the germ line of this cohort of families suggest that mutations in the promoter and 5' UTR play a very limited role in melanoma predisposition.

Evidence for time dependency of molecular rate estimates

Long-term changes in the genetic composition of a population occur by the fixation of new mutations, a process known as substitution. The rate at which mutations arise in a population and the rate at which they are fixed are expected to be equal under neutral conditions (Kimura, 1968). Between the appearance of a new mutation and its eventual fate of fixation or loss, there will be a period in which it exists as a transient polymorphism in the population (Kimura and Ohta, 1971). If the majority of mutations are deleterious (and nonlethal), the fixation probabilities of these transient polymorphisms are reduced and the mutation rate will exceed the substitution rate (Kimura, 1983). Consequently, different apparent rates may be observed on different time scales of the molecular evolutionary process (Penny, 2005; Penny and Holmes, 2001). The substitution rate of the mitochondrial protein-coding genes of birds and mammals has been traditionally recognized to be about 0.01 substitutions/site/million years (Myr) (Brown et al., 1979; Ho, 2007; Irwin et al., 1991; Shields and Wilson, 1987), with the noncoding D-loop evolving several times more quickly (e.g., Pesole et al., 1992; Quinn, 1992). Over the past decade, there has been mounting evidence that instantaneous mutation rates substantially exceed substitution rates, in a range of organisms (e.g., Denver et al., 2000; Howell et al., 2003; Lambert et al., 2002; Mao et al., 2006; Mumm et al., 1997; Parsons et al., 1997; Santos et al., 2005). The immediate reaction to the first of these findings was that the polymorphisms generated by the elevated mutation rate are short-lived, perhaps extending back only a few hundred years (Gibbons, 1998; Macaulay et al., 1997). That is, purifying selection was thought to remove these polymorphisms very rapidly.

Abrogation of contaminating RNA activity in HIV-1 Gag VLPs

Background: HIV-1 Gag virus like particles (VLPs) used as candidate vaccines are regarded as inert particles as they contain no replicative nucleic acid, although they do encapsidate cellular RNAs. During HIV-1 Gag VLP production in baculovirus-based expression systems, VLPs incorporate the baculovirus Gp64 envelope glycoprotein, which facilitates their entry into mammalian cells. This suggests that HIV-1 Gag VLPs produced using this system facilitate uptake and subsequent expression of encapsidated RNA in mammalian cells - an unfavourable characteristic for a vaccine. Methods. HIV-1 Gag VLPs encapsidating reporter chloramphenicol acetyl transferase (CAT) RNA, were made in insect cells using the baculovirus expression system. The presence of Gp64 on the VLPs was verified by western blotting and RT-PCR used to detect and quantitate encapsidated CAT RNA. VLP samples were heated to inactivate CAT RNA. Unheated and heated VLPs incubated with selected mammalian cell lines and cell lysates tested for the presence of CAT protein by ELISA. Mice were inoculated with heated and unheated VLPs using a DNA prime VLP boost regimen. Results: HIV-1 Gag VLPs produced had significantly high levels of Gp64 (∼1650 Gp64 molecules/VLP) on their surfaces. The amount of encapsidated CAT RNA/g Gag VLPs ranged between 0.1 to 7 ng. CAT protein was detected in 3 of the 4 mammalian cell lines incubated with VLPs. Incubation with heated VLPs resulted in BHK-21 and HeLa cell lysates showing reduced CAT protein levels compared with unheated VLPs and HEK-293 cells. Mice inoculated with a DNA prime VLP boost regimen developed Gag CD8 and CD4 T cell responses to GagCAT VLPs which also boosted a primary DNA response. Heating VLPs did not abrogate these immune responses but enhanced the Gag CD4 T cell responses by two-fold. Conclusions: Baculovirus-produced HIV-1 Gag VLPs encapsidating CAT RNA were taken up by selected mammalian cell lines. The presence of CAT protein indicates that encapsidated RNA was expressed in the mammalian cells. Heat-treatment of the VLPs altered the ability of protein to be expressed in some cell lines tested but did not affect the ability of the VLPs to stimulate an immune response when inoculated into mice. © 2011 Valley-Omar et al; licensee BioMed Central Ltd.

Single nucleotide polymorphism in hsa-mir-196a-2 and breast cancer risk : a case control study

microRNAs are small, non-coding RNAs that influence gene expression on a post-transcriptional level. They participate in diverse biological pathways and may act as either tumor suppressor genes or oncogenes. As they may have an effect on thousands of target mRNAs, single-nucleotide polymorphisms in microRNA genes might have major functional consequences, because the microRNA's properties and/or maturation may change. miR-196a has been reported to be aberrantly expressed in breast cancer tissue. Additionally, the SNP rs11614913 in hsa-mir-196a-2 has been found to be associated with breast cancer risk in some studies although not in others. This study evaluated the association between rs11614913 and breast cancer risk in a Caucasian case-control cohort in Queensland, Australia. Results do not support an association of the tested hsa-mir-196a-2 polymorphism with breast cancer susceptibility in this cohort. As there is a discrepancy between our results and previous findings, it is important to assess the role of rs11614913 in breast cancer by further larger studies investigating different ethnic groups.