Identification of putative 1p36 tumor suppressor genes involved in oligodendroglioma development

Oligodendrogliomas are primary neoplasms of the central nervous system (CNS). One of the most common and characteristic chromosomal abnormalities observed in oligodendroglioma is allelic loss of 1p (Reifenberger et al., 1994; Bello et al., 1995). Since 1p loss has been reported for both well-differentiated and anaplastic oligodendroglioma, it is believed to occur early in tumor development (Bello et al., 1995). This allelic loss also has clinical significance, for oligodendroglioma patients with 1p loss generally respond significantly better to combination chemotherapy and have longer average survival than do oligodendroglioma patients without 1p loss (Cairncross et al., 1998). To date, no genes on 1p have been implicated as essential to the development or treatment response of oligodendroglioma. In order to localize and/or identify a gene involved in oligodendroglioma development, I tested 170 oligodendrogliomas for deletions of 1p and tested 26 tumors for differential expression of genes in the region of 1p36. Evidence obtained from these methods implicated two genes, SHREW1 and the gene encoding DNA fragmentation factor beta (DFFB). The function for the SHREW1 locus is currently not well known, but preliminary data suggests that it a novel member of adherens junctions. The DFFB gene is an enhancer for apoptosis. Thus, both SHREW1 and DFFB may be candidates for an oligodendroglioma tumor suppressor. Mutational analysis of both genes did not uncover any mutations. Future studies will evaluate other mechanisms that may be responsible for inactivation of these genes in oligodendrogliomas. ^

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.

Identification of differentially expressed mRNA in prokaryotic organisms by customized amplification libraries (DECAL): The effect of isoniazid on gene expression in Mycobacterium tuberculosis

Understanding the effects of the external environment on bacterial gene expression can provide valuable insights into an array of cellular mechanisms including pathogenesis, drug resistance, and, in the case of Mycobacterium tuberculosis, latency. Because of the absence of poly(A)+ mRNA in prokaryotic organisms, studies of differential gene expression currently must be performed either with large amounts of total RNA or rely on amplification techniques that can alter the proportional representation of individual mRNA sequences. We have developed an approach to study differences in bacterial mRNA expression that enables amplification by the PCR of a complex mixture of cDNA sequences in a reproducible manner that obviates the confounding effects of selected highly expressed sequences, e.g., ribosomal RNA. Differential expression using customized amplification libraries (DECAL) uses a library of amplifiable genomic sequences to convert total cellular RNA into an amplified probe for gene expression screens. DECAL can detect 4-fold differences in the mRNA levels of rare sequences and can be performed on as little as 10 ng of total RNA. DECAL was used to investigate the in vitro effect of the antibiotic isoniazid on M. tuberculosis, and three previously uncharacterized isoniazid-induced genes, iniA, iniB, and iniC, were identified. The iniB gene has homology to cell wall proteins, and iniA contains a phosphopantetheine attachment site motif suggestive of an acyl carrier protein. The iniA gene is also induced by the antibiotic ethambutol, an agent that inhibits cell wall biosynthesis by a mechanism that is distinct from isoniazid. The DECAL method offers a powerful new tool for the study of differential gene expression.

Identification of genes differentially expressed by nerve growth factor- and neurotrophin-3-dependent sensory neurons

The neurotrophins nerve growth factor (NGF) and neurotrophin-3 (NT3) support the survival of subpopulations of primary sensory neurons with defined and distinct physiological characteristics. Only a few genes have been identified as being differentially expressed in these subpopulations, and not much is known about the nature of the molecules involved in the processing of sensory information in NGF-dependent nociceptive neurons or NT3-dependent proprioceptive neurons. We devised a simple dorsal root ganglion (DRG) explant culture system, allowing the selection of neuronal populations preferentially responsive to NGF or NT3. The reliability of this assay was first monitored by the differential expression of the NGF and NT3 receptors trkA and trkC, as well as that of neuropeptides and calcium-binding proteins. We then identified four differentially expressed sodium channels, two enriched in the NGF population and two others in the NT3 population. Finally, using an optimized RNA fingerprinting protocol, we identified 20 additional genes, all differentially expressed in DRG explants cultured with NGF or NT3. This approach thus allows the identification of large number of genes expressed in subpopulations of primary sensory neurons and opens the possibility of studying the molecular mechanisms of nociception and proprioception.

The Yeast Saccharomyces cerevisiae Contains Two Glutaredoxin Genes That Are Required for Protection against Reactive Oxygen Species

Glutaredoxins are small heat-stable proteins that act as glutathione-dependent disulfide oxidoreductases. Two genes, designated GRX1 and GRX2, which share 40–52% identity and 61–76% similarity with glutaredoxins from bacterial and mammalian species, were identified in the yeast Saccharomyces cerevisiae. Strains deleted for both GRX1 and GRX2 were viable but lacked heat-stable oxidoreductase activity using β-hydroxyethylene disulfide as a substrate. Surprisingly, despite the high degree of homology between Grx1 and Grx2 (64% identity), the grx1 mutant was unaffected in oxidoreductase activity, whereas the grx2 mutant displayed only 20% of the wild-type activity, indicating that Grx2 accounted for the majority of this activity in vivo. Expression analysis indicated that this difference in activity did not arise as a result of differential expression of GRX1 and GRX2. In addition, a grx1 mutant was sensitive to oxidative stress induced by the superoxide anion, whereas a strain that lacked GRX2 was sensitive to hydrogen peroxide. Sensitivity to oxidative stress was not attributable to altered glutathione metabolism or cellular redox state, which did not vary between these strains. The expression of both genes was similarly elevated under various stress conditions, including oxidative, osmotic, heat, and stationary phase growth. Thus, Grx1 and Grx2 function differently in the cell, and we suggest that glutaredoxins may act as one of the primary defenses against mixed disulfides formed following oxidative damage to proteins.

Lamprey Dlx genes and early vertebrate evolution

Gnathostome vertebrates have multiple members of the Dlx family of transcription factors that are expressed during the development of several tissues considered to be vertebrate synapomorphies, including the forebrain, cranial neural crest, placodes, and pharyngeal arches. The Dlx gene family thus presents an ideal system in which to examine the relationship between gene duplication and morphological innovation during vertebrate evolution. Toward this end, we have cloned Dlx genes from the lamprey Petromyzon marinus, an agnathan vertebrate that occupies a critical phylogenetic position between cephalochordates and gnathostomes. We have identified four Dlx genes in P. marinus, whose orthology with gnathostome Dlx genes provides a model for how this gene family evolved in the vertebrate lineage. Differential expression of these lamprey Dlx genes in the forebrain, cranial neural crest, pharyngeal arches, and sensory placodes of lamprey embryos provides insight into the developmental evolution of these structures as well as a model of regulatory evolution after Dlx gene duplication events.

The 5' coding region of Paramecium surface antigen genes controls mutually exclusive transcription.

Paramecium tetraurelia stock 51 can express at least 11 different types of surface antigens, yet only a single type is expressed on the surface of an individual cell at any one time. The differential expression of stock 51 type A and B surface antigen genes (51A and 51B) is regulated at the level of transcription. Previously, we reported that nucleotide sequences upstream of position -26 (relative to the start of translation) in the 51A and 51B surface antigen genes are necessary for transcriptional activity but are not sufficient to direct differential transcriptional control. In this report we demonstrate that at least some of the critical elements necessary for differential transcription of the 51A and 51B genes lie within the 5' coding region. A hybrid gene that contains 51B upstream sequences (-475 to +1) attached to the ATG start codon of 51A is not cotranscribed with the 51B gene. In contrast, further substitution with 51B sequences (-1647 to +885) allows the chimeric gene to be coexpressed with 51B. A different hybrid gene containing a substitution of 51B sequence from -26 to +885 in the 51A gene is also coexpressed with 51B, revealing that the critical elements within the coding region of 51B do not require 51B upstream sequences for their effect. Coinjection of the 51A gene with the chimeric gene that contains 51B up to +885 showed that the same sequences that allow coexpression with 51B prevent cotranscription with 51A. Together, these results demonstrate that a region downstream of the transcriptional start site between nucleotide positions +1 and +885 (relative to translational start) is necessary to control differential transcriptional activity.

Acute Synaptic Activity Causes Differential miRNA Expression in The Drosophila melanogaster Larval Central Nervous System

The primary goal of this thesis was to determine if spaced synaptic stimulation induced the differential expression of microRNAs (miRNAs) in the Drosophila melanogaster central nervous system (CNS). Prior to attaining this goal, we needed to identify and validate a spaced stimulation paradigm that could induce the formation of new synaptic growth at a model synapse, the larval neuromuscular junction (NMJ). Both Channelrhodopsin- and high potassium-based stimulation paradigms adapted from (Ataman, et al. 2008) were tested. Once validation of these paradigms was complete, we sought to characterize the miRNA expression profile of the larval CNS by miRNA array. Following attainment of these data, we used quantitative real-time PCR (RT-qPCR) to determine if acute synaptic stimulation caused the differential expression of neuronal miRNAs. We found that upon high potassium spaced training in a wild type (Canton S) genotype, 5 miRNAs showed significant differential expression when normalized to a validated reference gene, the U1 snRNA. Moreover, absolute quantification of our RT-qPCR study implicated one miRNA: miR-958 as being significantly regulated by activity. Investigation into potential targets for miR-958 revealed it to be a potential regular of Dlar, a protein tyrosine phosphatase implicated in synapse development. This investigation provides the foundation to directly test our underlying hypothesis that, following spaced training, differential expression of miRNAs alters the translation of proteins required to induce and maintain these structural changes at the synapse.

The Plasmodium serine-type SERA proteases display distinct expression patterns and non-essential in vivo roles during life cycle progression of the malaria parasite

Parasite proteases play key roles in several fundamental steps of the Plasmodium life cycle, including haemoglobin degradation, host cell invasion and parasite egress. Plasmodium exit from infected host cells appears to be mediated by a class of papain-like cysteine proteases called 'serine repeat antigens' (SERAs). A SERA subfamily, represented by Plasmodium falciparum SERA5, contains an atypical active site serine residue instead of a catalytic cysteine. Members of this SERAser subfamily are abundantly expressed in asexual blood stages, rendering them attractive drug and vaccine targets. In this study, we show by antibody localization and in vivo fluorescent tagging with the red fluorescent protein mCherry that the two P. berghei serine-type family members, PbSERA1 and PbSERA2, display differential expression towards the final stages of merozoite formation. Via targeted gene replacement, we generated single and double gene knockouts of the P. berghei SERAser genes. These loss-of-function lines progressed normally through the parasite life cycle, suggesting a specialized, non-vital role for serine-type SERAs in vivo. Parasites lacking PbSERAser showed increased expression of the cysteine-type PbSERA3. Compensatory mechanisms between distinct SERA subfamilies may thus explain the absence of phenotypical defect in SERAser disruptants, and challenge the suitability to develop potent antimalarial drugs based on specific inhibitors of Plasmodium serine-type SERAs.

Tissue-specific gene expression in soybean (Glycine max) detected by cDNA microarray analysis

We have constructed cDNA microarrays for soybean (Glycine max L. Merrill), containing approximately 4,100 Unigene ESTs derived from axenic roots, to evaluate their application and utility for functional genomics of organ differentiation in legumes. We assessed microarray technology by conducting studies to evaluate the accuracy of microarray data and have found them to be both reliable and reproducible in repeat hybridisations. Several ESTs showed high levels (>50 fold) of differential expression in either root or shoot tissue of soybean. A small number of physiologically interesting, and differentially expressed sequences found by microarray analysis were verified by both quantitative real-time RT-PCR and Northern blot analysis. There was a linear correlation (r(2) = 0.99, over 5 orders of magnitude) between microarray and quantitative real-time RT-PCR data. Microarray analysis of soybean has enormous potential not only for the discovery of new genes involved in tissue differentiation and function, but also to study the expression of previously characterised genes, gene networks and gene interactions in wild-type, mutant or transgenic; plants.

Identification of candidate genes at the corticoseptal boundary during development

Cortical midline glia are critical to the formation of the corpus callosum during development. The glial wedge is a Population of midline glia that is located at the corticoseptal boundary and expresses repulsive/growth-inhibitory molecules that guide callosal axons as they cross the midline. The glial wedge are the first cells within the cortex to express GFAP and thus may express molecules specific for glial maturation. The corticoseptal boundary is a genetically defined boundary between the cingulate cortex (dorsal telencephalon) and the septum (ventral telencephalon). The correct dorso-ventral position of this boundary is vital to the formation of both the glial wedge and the corpus callosum. Our aim was to identify genes expressed specifically within the glial wedge that might be involved in either glial differentiation, formation of the corticoseptal boundary or development of the corpus callosum. To identify such genes we have performed a differential display PCR screen comparing RNA isolated from the glial wedge with RNA isolated from control tissues such as the neocortex and septum, of embryonic day 17 mouse brains. Using 200 different combinations of primers, we identified and cloned 67 distinct gene fragments. In situ hybridization analysis confirmed the differential expression of many of the genes, and showed that clones G24F3, G39F8 and transcription factor LZIP have specific expression patterns in the telencephalon of embryonic and postnatal brains. An RNase Protection Assay (RPA) revealed that the expression of G39F8, G24173 and LZIP increase markedly in the telencephalon at E16 and continue to be expressed until at least PO, during the period when the corpus callosum is forming. (c) 2005 Elsevier B.V. All rights reserved.

Effects of lithium and valproic acid on gene expression and phenotypic markers in an NT2 neurosphere model of neural development

Mood stabilising drugs such as lithium (LiCl) and valproic acid (VPA) are the first line agents for treating conditions such as Bipolar disorder and Epilepsy. However, these drugs have potential developmental effects that are not fully understood. This study explores the use of a simple human neurosphere-based in vitro model to characterise the pharmacological and toxicological effects of LiCl and VPA using gene expression changes linked to phenotypic alterations in cells. Treatment with VPA and LiCl resulted in the differential expression of 331 and 164 genes respectively. In the subset of VPA targeted genes, 114 were downregulated whilst 217 genes were upregulated. In the subset of LiCl targeted genes, 73 were downregulated and 91 were upregulated. Gene ontology (GO) term enrichment analysis was used to highlight the most relevant GO terms associated with a given gene list following toxin exposure. In addition, in order to phenotypically anchor the gene expression data, changes in the heterogeneity of cell subtype populations and cell cycle phase were monitored using flow cytometry. Whilst LiCl exposure did not significantly alter the proportion of cells expressing markers for stem cells/undifferentiated cells (Oct4, SSEA4), neurons (Neurofilament M), astrocytes (GFAP) or cell cycle phase, the drug caused a 1.4-fold increase in total cell number. In contrast, exposure to VPA resulted in significant upregulation of Oct4, SSEA, Neurofilament M and GFAP with significant decreases in both G2/M phase cells and cell number. This neurosphere model might provide the basis of a human-based cellular approach for the regulatory exploration of developmental impact of potential toxic chemicals.

Human Synaptic Plasticity Gene Expression Profile and Dendritic Spine Density Changes in HIV-Infected Human CNS Cells: Role in HIV-Associated Neurocognitive Disorders (HAND)

HIV-associated neurocognitive disorders (HAND) is characterized by development of cognitive, behavioral and motor abnormalities, and occur in approximately 50% of HIV infected individuals. Our current understanding of HAND emanates mainly from HIV-1 subtype B (clade B), which is prevalent in USA and Western countries. However very little information is available on neuropathogenesis of HIV-1 subtype C (clade C) that exists in Sub-Saharan Africa and Asia. Therefore, studies to identify specific neuropathogenic mechanisms associated with HAND are worth pursuing to dissect the mechanisms underlying this modulation and to prevent HAND particularly in clade B infection. In this study, we have investigated 84 key human synaptic plasticity genes differential expression profile in clade B and clade C infected primary human astrocytes by using RT2 Profile PCR Array human Synaptic Plasticity kit. Among these, 31 and 21 synaptic genes were significantly (≥3 fold) down-regulated and 5 genes were significantly (≥3 fold) up-regulated in clade B and clade C infected cells, respectively compared to the uninfected control astrocytes. In flow-cytometry analysis, down-regulation of postsynaptic density and dendrite spine morphology regulatory proteins (ARC, NMDAR1 and GRM1) was confirmed in both clade B and C infected primary human astrocytes and SK-N-MC neuroblastoma cells. Further, spine density and dendrite morphology changes by confocal microscopic analysis indicates significantly decreased spine density, loss of spines and decreased dendrite diameter, total dendrite and spine area in clade B infected SK-N-MC neuroblastoma cells compared to uninfected and clade C infected cells. We have also observed that, in clade B infected astrocytes, induction of apoptosis was significantly higher than in the clade C infected astrocytes. In conclusion, this study suggests that down-regulation of synaptic plasticity genes, decreased dendritic spine density and induction of apoptosis in astrocytes may contribute to the severe neuropathogenesis in clade B infection.

Transcriptome-Wide Mapping of Pea Seed Ageing Reveals a Pivotal Role for Genes Related to Oxidative Stress and Programmed Cell Death

Understanding of seed ageing, which leads to viability loss during storage, is vital for ex situ plant conservation and agriculture alike. Yet the potential for regulation at the transcriptional level has not been fully investigated. Here, we studied the relationship between seed viability, gene expression and glutathione redox status during artificial ageing of pea (Pisum sativum) seeds. Transcriptome-wide analysis using microarrays was complemented with qRT-PCR analysis of selected genes and a multilevel analysis of the antioxidant glutathione. Partial degradation of DNA and RNA occurred from the onset of artificial ageing at 60% RH and 50 degrees C, and transcriptome profiling showed that the expression of genes associated with programmed cell death, oxidative stress and protein ubiquitination were altered prior to any sign of viability loss. After 25 days of ageing viability started to decline in conjunction with progressively oxidising cellular conditions, as indicated by a shift of the glutathione redox state towards more positive values (>-190 mV). The unravelling of the molecular basis of seed ageing revealed that transcriptome reprogramming is a key component of the ageing process, which influences the progression of programmed cell death and decline in antioxidant capacity that ultimately lead to seed viability loss.

Caracterização e comparação de genes expressos em ápices meristemáticos de cana-de-açúcar cultivados em SP e RN

In this work, we used sugarcane as a model due to its importance for sugar and ethanol production. Unlike the current plant models, sugarcane presents a complex genetics and an enormous allelic variation. Here, we report the analysis of SAGE libraries produced using the shoot apical meristem from contrasted genotypes by flowering induction (non-flowering vs. early-flowering varieties) grown under São Paulo state conditions. The expression pattern was analyzed using samples from São Paulo (SP) and Rio Grande do Norte (RN) states. These results showed that cDNAs identified by SAGE libraries had differential expression only in São Paulo state samples. Furthermore, the cDNA identified CYP (Citocrome P450) was chosen for in silico and genome characterization because it was found in SAGE libraries and subtractive libraries from samples from RN. Phylogenetic trees showed the relationship for these sequences. Furthermore, the qRT-PCR for CYP showed a potential role as flowering indutor for RN samples considering different isophorms. Considering the results present here, it can be consider that CYP gene may be used as molecular marker