Glioma-specific alternative RNA splicing: A role for polypyrimidine tract binding protein-1

Alternative RNA splicing plays an integral role in cell fate determination and function, especially in the cells of the brain. Errors in RNA processing contribute to diseases such as cancer, where it leads to the production of oncogenic proteins or the loss of tumor suppressors. In silica mining suggests that hundreds of splice isoforms are misexpressed in the glial cell-derived glioma. However, there is little experimental evidence of the prevalence and contribution of these changes and whether they contribute to the formation and progression of this devastating malignancy. To determine the frequency of these aberrant events, global profiling of alternative RNA splice patterns in glioma and nontumor brain was conducted using an exon array. Most splicing changes were less than 5-fold in magnitude and 14 cassette exon events were validated, including 7 previously published events. To determine the possible causes of missplicing, the differential expression levels of splicing factors in these two tissues were also analyzed. Six RNA splicing factors had greater than 2-fold changes in expression. The highest differentially expressed factor was polypyrimidine tract binding protein-1 (PTB). Evaluation by immunohistochemistry determined that this factor was elevated in both early and late stages of glioma. Glial cell-specific PTB expression in the adult brain led me to examine the role of PTB in gliomagenesis. Downregulation of PTB slowed glioma cell proliferation and migration and enhanced cell adhesion to fibronectin and vitronectin. To determine whether PTB was affecting these processes through splicing, genome-wide exon expression levels were correlated with PTB levels. Surprisingly, previously reported PTB target transcripts were insensitive to changes in PTB levels in both patient samples and PTB-depleted glioma cells. Only one validated glioma-specific splice target, RTN4/Nogo, had a significant PTB-mediated splicing change. Downregulation of PTB enhanced inclusion of its alternative exon 3, which encodes an auxiliary domain within a neurite inhibitor protein. Overexpression of this splice isoform in glioma cells slowed proliferation in a manner similar to that observed in PTB knockdown cells. In summary, aberrant expression of splicing factors such as PTB in glioma may elicit changes in splicing patterns that enhance tumorigenesis. ^

Micro RNA expression alterations in patients with glioblastoma

Glioblastoma, also known as glioblastoma multiform or GBM, is the most common and most malignant primary brain tumor. The clinical history of patients with glioblastoma is short, usually less than 3 months in more than 50% of cases after diagnosis. Currently, the methods of glioblastoma treatment are chemotherapy, radiotherapy and surgery. Even with the more effective treatment options, patients with glioblastoma most likely have a median survival time of 10 to 12 months. It is necessary to seek other treatment methods, including gene-targeted treatment. The success of gene-targeted treatment depends critically on the knowledge of genes that may be the cause of, or contribute to disease. To establish a correlate between glioblastoma survival timeline and micro RNA expression alteration, a study of 91 glioblastoma patients was conducted at the University of Texas M. D. Anderson Cancer Center. These 91 glioblastoma patients were newly diagnosed from 2002 to 2007. Statistical analysis was conducted to test the association of miRNA expression alteration between long-term survival and short-term survival glioblastoma. The completion of this proposed study will provide a better understanding of the regulatory role of miRNA in glioblastoma progression.^

GENETIC INTERACTIONS OF FACTORS THAT REGULATE ALTERNATIVE RNA SPLICING IN THE MALE GERM LINE OF DROSOPHILA

Alternative RNA splicing is a critical process that contributes variety to protein functions, and further controls cell differentiation and normal development. Although it is known that most eukaryotic genes produce multiple transcripts in which splice site selection is regulated, how RNA binding proteins cooperate to activate and repress specific splice sites is still poorly understood. In addition how the regulation of alternative splicing affects germ cell development is also not well known. In this study, Drosophila Transformer 2 (Tra2) was used as a model to explore both the mechanism of its repressive function on its own pre-mRNA splicing, and the effect of the splicing regulation on spermatogenesis in testis. Half-pint (Hfp), a protein known as splicing activator, was identified in an S2 cell-based RNAi screen as a co-repressor that functions in combination with Tra2 in the splicing repression of the M1 intron. Its repressive splicing function is found to be sequence specific and is dependent on both the weak 3’ splice site and an intronic splicing silencer within the M1 intron. In addition we found that in vivo, two forms of Hfp are expressed in a cell type specific manner. These alternative forms differ at their amino terminus affecting the presence of a region with four RS dipeptides. Using assays in Drosophila S2 cells, we determined that the alternative N terminal domain is necessary in repression. This difference is probably due to differential localization of the two isoforms in the nucleus and cytoplasm. Our in vivo studies show that both Hfp and Tra2 are required for normal spermatogenesis and cooperate in repression of M1 splicing in spermatocytes. But interestingly, Tra2 and Hfp antagonize each other’s function in regulating germline specific alternative splicing of Taf1 (TBP associated factor 1). Genetic and cytological studies showed that mutants of Hfp and Taf1 both cause similar defects in meiosis and spermatogenesis. These results suggest Hfp regulates normal spermatogenesis partially through the regulation of taf1 splicing. These observations indicate that Hfp regulates tra2 and taf1 activity and play an important role in germ cell differentiation of male flies.

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.

Estrogen induced desensitization of c -{\it fos\/} messenger RNA expression {\it in vivo\/}

In this thesis, we investigated the regulation of the nuclear proto-oncogene, c-fos by estrogen in vivo. In the uterus, estrogen causes a rapid, dramatic and transient induction of c-fos mRNA and this occurs by transcriptional activation. We have discovered a previously unrecognized regulatory mechanism by which fos becomes desensitized to estrogen following the transient induction. We investigated three aspects of this desensitization: (1) the kinetics and general characteristics of the phenomenon; (2) the molecular mechanism of the desensitization; and (3) the relationship of desensitization to estrogen stimulated DNA synthesis. The desensitization occurs between 3-24 hours after initial hormonal stimulation and is reversible within 72 hours. The desensitization is not species specific, in that it occurs in both the rat and mouse. The desensitization also occurs in at least two estrogen responsive tissues, the uterus and vagina. The desensitization is not unique to c-fos, since both c-myc and c-jun show similar patterns of desensitization. However, the desensitization is not observed with creatine kinase B (CKB), indicating that not all estrogen inducible genes become desensitized. In the second general area, we determined the desensitization is at the transcriptional level. The desensitization is homologous, but not heterologous, since estrogen induction does not desensitize c-fos to other agents. Other studies show that the desensitization is not due to the lack of functional estrogen receptors. Taken together, these findings suggest that the desensitization occurs at the level of the estrogen responsive element. In the third major area, we demonstrated that the desensitization appears to be related to estrogen induced DNA synthesis. Support for this suggestion comes from the observation that short acting estrogens which induce fos, but not DNA synthesis, do not produce desensitization. ^

{\it cis\/} -acting determinants in the control of MuSVts110 RNA splicing

Like other simple retroviruses the murine sarcoma virus ts110 (MuSVts110) displays an inefficient mode of genome splicing. But, unlike the splicing phenotypic of other retroviruses, the splicing event effected upon the transcript of MuSVts110 is temperature sensitive. Previous work in this laboratory has established that the conditionally defective nature of MuSVts110 RNA splicing is mediated in cis by features in the viral transcript. Here we show that the 5$\sp\prime$ splice site of the MuSVts110 transcript acts as a point of control of the overall splicing efficiency at both permissive and nonpermissive temperatures for splicing. We strengthened and simultaneously weakened the nucleotide structure of the 5$\sp\prime$ splice site in an attempt to elucidate the differential effects each of the two known critical splicing components which interact with the 5$\sp\prime$ splice site have on the overall efficiency of intron excision. We found that a transversion of the sixth nucleotide, resulting in the formation of a near-consensus 5$\sp\prime$ splice site, dramatically increased the overall efficiency of MuSVts110 RNA splicing and abrogated the thermosensitive nature of this splicing event. Various secondary mutations within this original transversion mutant, designed to selectively decrease specific splicing component interactions, lead to recovery of inefficient and thermosensitive splicing. We have further shown that a sequence of 415 nucleotides lying in the downstream exon of the viral RNA and hypothesized to act as an element in the temperature-dependent inhibition of splicing displays a functional redundancy throughout its length; loss and/or replacement of any one sequence of 100 nucleotides within this sequence does not, with one exception detailed below, diminish the degree to which MuSVts110 RNA is inhibited to splice at the restrictive temperature. One specific deletion, though, fortuitously juxtaposed and activated cryptic consensus splicing signals for the excision of a cryptic intron within the downstream exon and markedly potentiated--across a newly defined cryptic exon--the splicing event effected upon the upstream, native intron. We have exploited this mutant of MuSVts110 to further an understanding of the process of exon definition and intron definition and show that the polypyrimidine tract and consensus 3$\sp\prime$ splice site, as well as the 5$\sp\prime$ splice site, within the intron at the 3$\sp\prime$ flank of the defined exon are required for the exon's definition; implying that definition of the downstream intron is required for the in vivo definition of the proximal, upstream exon. Finally; we have shown, through the construction of heterologous mutants of MuSVts110 employing a foreign 3$\sp\prime$ end-forming sequence, that efficiency of transcript splicing can be increased--to a degree which abrogates its thermosensitive nature--in direct proportion to increasing proximity of the 3$\sp\prime$ end-forming signal to the terminal 3$\sp\prime$ splice site. ^

Thermo -modulation of MuSVts110 splicing by inhibitory RNA sequences

Viral systems have contributed tremendously to the understanding of eukaryotic molecular biology. The proportional pattern of retroviral RNA expression offers many clues into the alternative splicing of cellular transcripts. The MuSVts110 virus presents an unusual expression system, where the mechanistic combination of RNA splicing and cellular transformation can be physiologically manipulated. Splicing of MuSVts110 pre-mRNA occurs inefficiently (30%-50%) at 33$\sp\circ$C or below and is subdued at 39$\sp\circ$C ($<$5%). Like most alternatively spliced cellular and retroviral transcripts, the MuSVts110 pre-mRNA contains cis-acting intron and exon sequences that attenuate splicing. These include a splicing inhibitory sequence at the 3$\prime$ end of the MuSVts110 v-mos exon, called the E2 Distal Element (E2DE), and a sub-optimal 3$\prime$ splice site. The E2DE directly inhibits MuSVts110 RNA splicing in a sequence-specific fashion at 39$\sp\circ$C but not at 28$\sp\circ$C, potentially through the association of cellular factors. Inefficient MuSVts110 splicing is pre-dominantly attributed to the utilization of multiple weak branchpoint sequences located between $-113$ and $-34$ nucleotides upstream of the 3$\prime$ splice site. The molecular control of MuSVts110 splicing, represented primarily by scattered multiple inefficient branchpoint sequences that are conditionally modulated by the E2DE at higher growth temperatures, is discussed. ^

Mechanism of mRNA stability control mediated by AU -rich element: Characterization of cis-elements and trans-factors

Regulation of cytoplasmic deadenylation, the first step in mRNA turnover, has direct impact on the fate of gene expression. AU-rich elements (AREs) found in the 3′ untranslated regions of many labile mRNAs are the most common RNA-destabilizing elements known in mammalian cells. Based on their sequence features and functional properties, AREs can be divided into three classes. Class I or class III ARE directs synchronous deadenylation, whereas class II ARE directs asynchronous deadenylation with the formation of poly(A)-intermediates. Through systematic mutagenesis study, we found that a cluster of five or six copies of AUUUA motifs forming various degrees of reiteration is the key feature dictating the choice between asynchronous versus synchronous deadenylation. A 20–30 nt AU-rich sequence immediately 5 ′ to this cluster of AUUUA motifs can greatly enhance its destabilizing ability and is an integral part of the AREs. These two features are the defining characteristics of class II AREs. ^ To better understand the decay mechanism of AREs, current methods have several limitations. Taking the advantage of tetracycline-regulated promoter, we developed a new transcriptional pulse strategy, Tet-system. By controlling the time and the amount of Tet addition, a pulse of RNA could be generated. Using this new system, we showed that AREs function in both growth- and density-arrested cells. The new strategy offers for the first time an opportunity to investigate control of mRNA deadenylation and decay kinetics in mammalian cells that exhibit physiologically relevant conditions. ^ As a member of heterogeneous nuclear RNA-binding protein, hnRNP D 0/AUF1 displays specific affinities for ARE sequences in vitro . But its in vivo function in ARE-mediated mRNA decay is unclear. AUF1/hnRNP D0 is composed of at least four isoforms derived by alternative RNA splicing. Each isoform exhibits different affinity for ARE sequence in vitro. Here, we examined in vivo effect of AUF1s/hnRNP D0s on degradation of ARE-containing mRNA. Our results showed that all four isoforms exhibit various RNA stabilizing effects in NIH3T3 cells, which are positively correlated with their binding affinities for ARE sequences. Further experiments indicated that AUF1/hnRNP D0 has a general role in modulating the stability of cytoplasmic mRNAs in mammalian cells. ^

Searching for the ideal DNA recognition code of covalent ligands

The combitiatorial approach restriction endonuclease protection selection and amplification REPSA was successfully used to determine ideal DNA interactions sites of covalent ligands. Unlike most other combinatorial methods, REPSA is based on inhibition of enzymatic cleavage by specific ligand-DNA complexes, which enables identification of binding sites of various ligands. However, the inherent nature of this technique posses a problem during selection of binding sites of covalent ligands. By modifying the technique according to the nature of the ligand, we demonstrate the flexibility of REPSA in identifying the preferred binding sites for monocovalent ligands, topoisomerase I and tallimustine, and the bicovalent ligand topoisomerase II. From among the preferred binding sites, we identified the consensus binding sequence of camptothecin induced topoisomerase I cleavage as ‘aGWT/Gc’, and tallimustine consensus sequences as ‘GTTCTA’ and ‘TTTTTTC’. We have shown for the first time that preferential binding of tallimustine occurs at sequences not previously reported. Furthermore, our data indicate that tallimustine is a novel DNA minor groove, guanine-specific alkylating agent. ^ Additionally, we have demonstrated in vivo that sequence-specific covalent DNA-binding small molecules have the ability to regulate transcription by inhibiting RNA polymerase II. Tallimustine, binding to its preferred sequences located in the 5′ untranslated region were an effective impediment for transcribing polymerase II. The ability of covalent binding small molecules to target predetermined DNA sequences located downstream of the promoter suggests a general approach for regulation of gene expression. ^

Regulations of SOX9 activity during chondrogenesis

Sox9 is a transcription factor required for chondrocyte differentiation and cartilage formation. In an effort to identify SOX9 interacting protein(s), we screened a chondrocyte cDNA library with a modified yeast two-hybrid method, Son of Sevenless (SOS) recruitment system (SRS). The catalytic subunit of cyclic AMP-dependent protein kinase A (PKA-Cα) and a new long form of c-Maf transcription factor (Lc-Maf) were found to interact specifically with SOX9. We showed here that two PKA phosphorylation consensus sites of SOX9 could be phosphorylated by PKA in vitro as well as in vivo. PKA phosphorylation of SOX9 increases its DNA binding and transcriptional activities on a Col2a1 chondrocyte-specific enhancer. Mutations of these two PKA phosphorylation sites markedly decreased the activation of SOX9 by PKA. ^ To test whether parathyroid hormone-related peptide (PTHrP) signaling results in SOX9 phosphorylation, we generated a phosphospecific antibody that specifically recognizes SOX9 that is phosphorylated at serine 181 (S 181) one of the two consensus PKA phosphorylation sites. Addition of PTHrP to COS7 cells cotransfected with SOX9 and PTH/PTHrP receptor strongly increased phosphorylation of SOX9 at S181; this phosphorylation was blocked by a PKA-specific inhibitor. In similar experiments we showed that PTHrP increased the activity of a SOX9-dependent Col2a1 enhancer. This increase in activity was abolished when a SOX9 mutant was used containing serine-to-alanine substitution in the two consensus PKA phosphorylation sites of SOX9. Using our phosphospecific SOX9 antibody we showed by immunohistochemistry of mouse embryos that Sox9 phosphorylated at S181 was localized almost exclusively in the pre-hypertrophic zone of the growth plate, an area corresponding to the major site of expression of PTH/PTHrP receptor. In contrast, no phosphorylation of Sox9 at S181 was detected in growth plates of PTH/PTHrP receptor null mutant mice. Sox9, regardless of phosphorylation state, was present in all chondrocytes of both genotypes except in hypertrophic chondrocytes. Thus, Sox9 is a target of PTHrP signaling and the PTHrP-dependent phosphorylation of SOX9 enhances its transcriptional activity. ^ In order to investigate the in vivo function of Sox9 phosphorylation by PKA, we are generating a mouse model of mutant Sox9 harboring point mutations in two PKA phosphorylation sites. Preliminary results indicated that heterozygous mice containing half amount of mutant Sox9 that can not be phosphorylated by PKA have normal skeletal phenotype and homozygous mice are being generated. ^ Lc-Maf encodes an extra ten amino acids at the carboxyl terminus of c-Maf and contains a completely different 3′ untranslated region. The interaction between SOX9 and Lc-Maf was further confirmed by co-immunoprecipitation and GST-pull down assays, which mapped the interacting domains of SOX9 to HMG DNA binding domain and that of Lc-Maf to basic leusine zipper motif. In situ hybridizations showed that RNA of Lc-Maf coexpressed with those of Sox9 and Col2a1 in areas of mesenchymal condensation during the early stages of mouse embryo development. A DNA binding site of Lc-Maf was identified at the 5′ part of a 48-bp Col2a1 enhancer element near the HMG binding site of SOX9. Lc-Maf and SOX9 synergistically activated a luciferase reporter plasmid containing a Col2al enhancer and increased the transcription of endogenous Col2a1 gene. In summary, Lc-Maf is the first identified SOX9-interating protein during chondrogenesis and may be an important activator of Col2a1 gene. ^

Cross -packing among retroviruses and characterization of the feline immunodeficiency virus (FIV) packaging signal: Implications for the development of FIV-based gene transfer systems

Feline immunodeficiency virus (FIV)-based gene transfer systems are being seriously considered for human gene therapy as an alternative to vectors based on primate lentiviruses, a genetically complex group of retroviruses capable of infecting non-dividing cells. The greater phylogenetic distance between the feline and primate lentiviruses is thought to reduce chances of the generation of recombinant viruses. However, safety of FIV-based vector systems has not been tested experimentally. Since primate lentiviruses such as human and simian immunodeficiency viruses (HIV/SIV) can cross-package each other's genomes, we tested this trait with respect to FIV. Unexpectedly, both feline and primate lentiviruses were reciprocally able to both cross-package and propagate each other's RNA genomes. This was largely due to the recognition of viral packaging signals by the heterologous proteins. However, a simple retrovirus such as Mason-Pfizer monkey virus (MPMV) was unable to package FIV RNA. Interestingly, FIV could package MPMV RNA, but not propagate it for further steps of replication. These findings suggest that upon co-infection of the same host, cross-packaging may allow distinct retroviruses to generate chimeric variants with unknown pathogenic potential. ^ In order to understand the packaging determinants in FIV, we conducted a detailed mutational analysis of the region thought to contain FIV packaging signal. We show that the first 90–120 nt of the 5′ untranslated region (UTR) and the first 90 nt of gag were simultaneously required for efficient FIV RNA packaging. These results suggest that the primary FIV packaging signal is multipartite and discontinuous, composed of two core elements separated by 150 nt of the 5 ′UTR. ^ The above studies are being used towards the development of safer FIV-based self-inactivating (SIN) vectors. These vectors are being designed to eliminate the ability of FIV transfer vector RNAs to be mobilized by primate lentiviral proteins that may be present in the target cells. Preliminary test of the first generation of these vectors has revealed that they are incapable of being propagated by feline proteins. The inability of FIV transfer vectors to express packageable vector RNA after integration should greatly increase the safety of FIV vectors for human gene therapy. ^

Seawater carbonate chemistry, clearance rates, respiration rates, condition index and cellular turnover (RNA: DNA) of Juvenile King Scallop, Pecten maximus in a laboratory experiment

The decline in ocean water pH and changes in carbonate saturation states through anthropogenically mediated increases in atmospheric CO2 levels may pose a hazard to marine organisms. This may be particularly acute for those species reliant on calcareous structures like shells and exoskeletons. This is of particular concern in the case of valuable commercially exploited species such as the king scallop, Pecten maximus. In this study we investigated the effects on oxygen consumption, clearance rates and cellular turnover in juvenile P. maximus following 3 months laboratory exposure to four pCO2 treatments (290, 380, 750 and 1140 µatm). None of the exposure levels were found to have significant effect on the clearance rates, respiration rates, condition index or cellular turnover (RNA: DNA) of individuals. While it is clear that some life stages of marine bivalves appear susceptible to future levels of ocean acidification, particularly under food limiting conditions, the results from this study suggest that where food is in abundance, bivalves like juvenile P. maximus may display a tolerance to limited changes in seawater chemistry.