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.

Gene expression and physiological changes of the Antarctic krill Euphausia superba under different hypoxia intensities

Antarctic krill (Euphausia superba) from South Georgia comprise one of the most northern and abundant krill stocks. South Georgia waters are undergoing rapid warming, as a result of climate change, which in turn could alter the oxygen concentration of the water. We investigated gene expression in Antarctic krill related to aerobic metabolism, antioxidant defence, and heat-shock response under severe (2.5% O2 saturation or 0.6 kPa) and threshold (20% O2 saturation or 4 kPa) hypoxia exposure compared to in situ levels (normoxic; 100% O2 saturation or 21 kPa). Biochemical metabolic and oxidative stress indicators complemented the genic expression analysis to detect in vivo signs of stress during the hypoxia treatments. Expression levels of the genes citrate synthase (CS), mitochondrial manganese superoxide dismutase (SODMn-m) and one heat-shock protein isoform (E) were higher in euphausiids incubated 6 h at 20% O2 saturation than in animals exposed to control (normoxic) conditions. All biochemical antioxidant defence parameters remained unchanged among treatments. Levels of lipid peroxidation were raised after 6 h of severe hypoxia. Overall, short-term exposure to hypoxia altered mitochondrial metabolic and antioxidant capacity, but did not induce anaerobic metabolism. Antarctic krill are swarming organisms and may experience short periods of hypoxia when present in dense swarms. A future, warmer Southern ocean, where oxygen saturation levels are decreased, may result in smaller, less dense swarms as they act to avoid greater levels of hypoxia.

The family of DOF transcription factors in Brachypodium distachyon: phylogenetic comparison with rice and barley DOFs and expression profiling

ABSTRACT: Transcription factors (TFs) are proteins that have played a central role both in evolution and in domestication, and are major regulators of development in living organisms. Plant genome sequences reveal that approximately 7% of all genes encode putative TFs. The DOF (DNA binding with One Finger) TF family has been associated with vital processes exclusive to higher plants and to their close ancestors (algae, mosses and ferns). These are seed maturation and germination, light-mediated regulation, phytohormone and plant responses to biotic and abiotic stresses, etc. In Hordeum vulgare and Oryza sativa, 26 and 30 different Dof genes, respectively, have been annotated. Brachypodium distachyon has been the first Pooideae grass to be sequenced and, due to its genomic, morphological and physiological characteristics, has emerged as the model system for temperate cereals, such as wheat and barley. RESULTS: Through searches in the B. distachyon genome, 27 Dof genes have been identified and a phylogenetic comparison with the Oryza sativa and the Hordeum vulgare DOFs has been performed. To explore the evolutionary relationship among these DOF proteins, a combined phylogenetic tree has been constructed with the Brachypodium DOFs and those from rice and barley. This phylogenetic analysis has classified the DOF proteins into four Major Cluster of Orthologous Groups (MCOGs). Using RT-qPCR analysis the expression profiles of the annotated BdDof genes across four organs (leaves, roots, spikes and seeds) has been investigated. These results have led to a classification of the BdDof genes into two groups, according to their expression levels. The genes highly or preferentially expressed in seeds have been subjected to a more detailed expression analysis (maturation, dry stage and germination). CONCLUSIONS: Comparison of the expression profiles of the Brachypodium Dof genes with the published functions of closely related DOF sequences from the cereal species considered here, deduced from the phylogenetic analysis, indicates that although the expression profile has been conserved in many of the putative orthologs, in some cases duplication followed by subsequent divergence may have occurred (neo-functionalization).

Genomic analysis of human and mouse TCL1 loci reveals a complex of tightly clustered genes

TCL1 and TCL1b genes on human chromosome 14q23.1 are activated in T cell leukemias by translocations and inversions at 14q32.1, juxtaposing them to regulatory elements of T cell receptor genes. In this report we present the cloning, mapping, and expression analysis of the human and murine TCL1/Tcl1 locus. In addition to TCL1 and TCL1b, the human locus contains two additional genes, TCL1-neighboring genes (TNG) 1 and 2, encoding proteins of 141 and 110 aa, respectively. Both genes show no homology to any known genes, but their expression profiles are very similar to those of TCL1 and TCL1b. TNG1 and TNG2 also are activated in T cell leukemias with rearrangements at 14q32.1. To aid in the development of a mouse model we also have characterized the murine Tcl1 locus and found five genes homologous to human TCL1b. Tcl1b1–Tcl1b5 proteins range from 117 to 123 aa and are 65–80% similar, but they show only a 30–40% similarity to human TCL1b. All five mouse Tcl1b and murine Tcl1 mRNAs are abundant in mouse oocytes and two-cell embryos but rare in various adult tissues and lymphoid cell lines. These data suggest a similar or complementary function of these proteins in early embryogenesis.

Development of an epithelium-specific expression cassette with human DNA regulatory elements for transgene expression in lung airways

The efficient expression of therapeutic genes in target cells or tissues is an important component of efficient and safe gene therapy. Utilizing regulatory elements from the human cytokeratin 18 (K18) gene, including 5′ genomic sequences and one of its introns, we have developed a novel expression cassette that can efficiently express reporter genes, as well as the human cystic fibrosis transmembrane conductance regulator (CFTR) gene, in cultured lung epithelial cells. CFTR transcripts expressed from the native K18 enhancer/promoter include two alternative splicing products, due to the activation of two cryptic splice sites in the CFTR coding region. Modification of the K18 intron and CFTR cDNA sequences eliminated the cryptic splice sites without changing the CFTR amino acid sequence, and led to enhanced CFTR mRNA and protein expression as well as biological function. Transgenic expression analysis in mice showed that the modified expression cassette can direct efficient and epithelium-specific expression of the Escherichia coli LacZ gene in the airways of fetal lungs, with no detectable expression in lung fibroblasts or endothelial cells. This is the first expression cassette which selectively directs lung transgene expression for CFTR gene therapy to airway epithelia.

Remodeling of Yeast Genome Expression in Response to Environmental Changes

We used genome-wide expression analysis to explore how gene expression in Saccharomyces cerevisiae is remodeled in response to various changes in extracellular environment, including changes in temperature, oxidation, nutrients, pH, and osmolarity. The results demonstrate that more than half of the genome is involved in various responses to environmental change and identify the global set of genes induced and repressed by each condition. These data implicate a substantial number of previously uncharacterized genes in these responses and reveal a signature common to environmental responses that involves ∼10% of yeast genes. The results of expression analysis with MSN2/MSN4 mutants support the model that the Msn2/Msn4 activators induce the common response to environmental change. These results provide a global description of the transcriptional response to environmental change and extend our understanding of the role of activators in effecting this response.

Involvement of the Octadecanoid Pathway and Protein Phosphorylation in Fungal Elicitor-Induced Expression of Terpenoid Indole Alkaloid Biosynthetic Genes in Catharanthus roseus

Two key genes in terpenoid indole alkaloid biosynthesis, Tdc and Str, encoding tryptophan decarboxylase and strictosidine synthase, respectively, are coordinately induced by fungal elicitors in suspension-cultured Catharanthus roseus cells. We have studied the roles of the jasmonate biosynthetic pathway and of protein phosphorylation in signal transduction initiated by a partially purified elicitor from yeast extract. In addition to activating Tdc and Str gene expression, the elicitor also induced the biosynthesis of jasmonic acid. The jasmonate precursor α-linolenic acid or methyl jasmonate (MeJA) itself induced Tdc and Str gene expression when added exogenously . Diethyldithiocarbamic acid, an inhibitor of jasmonate biosynthesis, blocked both the elicitor-induced formation of jasmonic acid and the activation of terpenoid indole alkaloid biosynthetic genes. The protein kinase inhibitor K-252a abolished both elicitor-induced jasmonate biosynthesis and MeJA-induced Tdc and Str gene expression. Analysis of the expression of Str promoter/gusA fusions in transgenic C. roseus cells showed that the elicitor and MeJA act at the transcriptional level. These results demonstrate that the jasmonate biosynthetic pathway is an integral part of the elicitor-triggered signal transduction pathway that results in the coordinate expression of the Tdc and Str genes and that protein kinases act both upstream and downstream of jasmonates.

Changes in global gene expression patterns during development and maturation of the rat kidney

We set out to define patterns of gene expression during kidney organogenesis by using high-density DNA array technology. Expression analysis of 8,740 rat genes revealed five discrete patterns or groups of gene expression during nephrogenesis. Group 1 consisted of genes with very high expression in the early embryonic kidney, many with roles in protein translation and DNA replication. Group 2 consisted of genes that peaked in midembryogenesis and contained many transcripts specifying proteins of the extracellular matrix. Many additional transcripts allied with groups 1 and 2 had known or proposed roles in kidney development and included LIM1, POD1, GFRA1, WT1, BCL2, Homeobox protein A11, timeless, pleiotrophin, HGF, HNF3, BMP4, TGF-α, TGF-β2, IGF-II, met, FGF7, BMP4, and ganglioside-GD3. Group 3 consisted of transcripts that peaked in the neonatal period and contained a number of retrotransposon RNAs. Group 4 contained genes that steadily increased in relative expression levels throughout development, including many genes involved in energy metabolism and transport. Group 5 consisted of genes with relatively low levels of expression throughout embryogenesis but with markedly higher levels in the adult kidney; this group included a heterogeneous mix of transporters, detoxification enzymes, and oxidative stress genes. The data suggest that the embryonic kidney is committed to cellular proliferation and morphogenesis early on, followed sequentially by extracellular matrix deposition and acquisition of markers of terminal differentiation. The neonatal burst of retrotransposon mRNA was unexpected and may play a role in a stress response associated with birth. Custom analytical tools were developed including “The Equalizer” and “eBlot,” which contain improved methods for data normalization, significance testing, and data mining.

Expression of biomineralization-related ion transport genes in Emiliania huxleyi

Biomineralization in the marine phytoplankton Emiliania huxleyi is a stringently controlled intracellular process. The molecular basis of coccolith production is still relatively unknown although its importance in global biogeochemical cycles and varying sensitivity to increased pCO2 levels has been well documented. This study looks into the role of several candidate Ca2+, H+ and inorganic carbon transport genes in E. huxleyi, using quantitative reverse transcriptase PCR. Differential gene expression analysis was investigated in two isogenic pairs of calcifying and non-calcifying strains of E. huxleyi and cultures grown at various Ca2+ concentrations to alter calcite production. We show that calcification correlated to the consistent upregulation of a putative HCO3- transporter belonging to the solute carrier 4 (SLC4) family, a Ca2+/H+ exchanger belonging to the CAX family of exchangers and a vacuolar H+-ATPase. We also show that the coccolith-associated protein, GPA is downregulated in calcifying cells. The data provide strong evidence that these genes play key roles in E. huxleyi biomineralization. Based on the gene expression data and the current literature a working model for biomineralization-related ion transport in coccolithophores is presented.

Probing the S100 protein family through genomic and functional analysis

The EF-hand superfamily of calcium binding proteins includes the S100, calcium binding protein, and troponin subfamilies. This study represents a genome, structure, and expression analysis of the S100 protein family, in mouse, human, and rat. We confirm the high level of conservation between mammalian sequences but show that four members, including S100A12, are present only in the human genome. We describe three new members of the S100 family in the three species and their locations within the S100 genomic clusters and propose a revised nomenclature and phylogenetic relationship between members of the EF-hand superfamily. Two of the three new genes were induced in bone-marrow-derived macrophages activated with bacterial lipopolysaccharide, suggesting a role in inflammation. Normal human and murine tissue distribution profiles indicate that some members of the family are expressed in a specific manner, whereas others are more ubiquitous. Structure-function analysis of the chemotactic properties of murine S100A8 and human S100A12, particularly within the active hinge domain, suggests that the human protein is the functional homolog of the murine protein. Strong similarities between the promoter regions of human S100A12 and murine S100A8 support this possibility. This study provides insights into the possible processes of evolution of the EF-hand protein superfamily. Evolution of the S100 proteins appears to have occurred in a modular fashion, also seen in other protein families such as the C2H2-type zinc-finger family. (C) 2004 Elsevier Inc. All rights reserved.

Correlating gene expression with larval competence and the effect of age and parentage on metamorphosis in the tropical abalone Haliotis asinina

The non-geniculate crustose coralline alga (CCA) Mastophora pacifica can induce the metamorphosis of competent Haliotis asinina (Vetigastropoda) larvae. The ability to respond to this natural cue varies considerably with larval age, with a higher proportion of older larvae (e.g. 90 h) able to metamorphose in response to M. pacifica than younger larvae (e.g. 66 h). Here we document the variation in time to acquisition of competence within a larval age class. For example, after 18 h of exposure to M. pacifica, approximately 15 and 36% of 84 and 90-h-old H. asinina larvae had initiated metamorphosis, respectively. This age-dependent response to M. pacifica is also observed when different aged larvae are exposed to CCA for varying periods. A higher proportion of older larvae require shorter periods of exposure to CCA than younger larvae in order to initiate metamorphosis. In this experiment, as in the previous, a small proportion of young larvae were able to respond to brief periods of CCA exposure, suggesting that they had developed the same state of competency as the majority of their older counterparts. Comparisons of the proportions of larvae undergoing metamorphosis between families reveals that parentage also has a significant (P < 0.05) affect on whether an individual will initiate metamorphosis at a given age. These familial differences are more pronounced when younger, largely pre-competent larvae (i.e. 66 h old) are exposed to M. pacifica, with proportions of larvae undergoing metamorphosis differing by as much as 10 fold between families. As these data suggest that variation in the rate of development of the competent state has a genetic basis, and as a first step towards identifying the molecular basis to this variation, we have identified numerous genes that are differentially expressed later in larval development using a differential display approach. Spatial expression analysis of these genes suggests that they may be directly involved in the acquisition of competence, or may play a functional role in the postlarva following metamorphosis.

Expression profiling in alcoholism research

This article represents the proceedings of a symposium at the 2004 International Society for Biomedical Research on Alcoholism in Mannheim, Germany, organized and co-chaired by Susan E. Bergeson and Wolfgang Sommer. The presentations and presenter were (1) Gene Expression in Brains of AlcoholPreferring and Non-Preferring Rats, by Howard J. Edenberg (2) Candidate Treatment Targets for Alcoholism: Leads from Functional Genomics Approaches, by Wolfgang Sommer (3) Microarray Analysis of Acute and Chronic Alcohol Response in Brain, by Susan E. Bergeson (4) On the Integration of QTL and Gene Expression Analysis, by Robert J. Hitzemann (5) Microarray and Proteomic Analysis of the Human Alcoholic Brain, by Peter R. Dodd.

A novel cis-acting element, ESP, contributes to high-level endosperm-specific expression in an oat globulin promoter

To examine the genetic controls of endosperm (ES) specificity, several cereal seed storage protein (SSP) promoters were isolated and studied using a transient expression analysis system. An oat globulin promoter (AsGlo1) capable of driving strong ES-specific expression in barley and wheat was identified. Progressive 5' deletions and cis element mutations demonstrated that the mechanism of specificity in the AsGlo1 promoter was distinct from that observed in glutelin and prolamin promoters. A novel interrupted palindromic sequence, ACATGTCAT-CATGT, was required for ES specificity and substantially contributed to expression strength of the AsGlo1 promoter. This sequence was termed the endosperm specificity palindrome (ESP) element. The GCN4 element, which has previously been shown to be required for ES specificity in cereal SSP promoters, had a quantitative role but was not required for tissue specificity. The 960-bp AsGlo1 promoter and a 251-bp deletion containing the ESP element also drove ES-specific expression in stably transformed barley. Reporter gene protein accumulated at very high levels (10% of total soluble protein) in ES tissues of plants transformed with an AsGlo1:GFP construct. Expression strength and tissue specificity were maintained over five transgenic generations. These attributes make the AsGlo1 promoter an ideal promoter for biotechnology applications. In conjunction with previous findings, our data demonstrate that there is more than one genetically distinct mechanism by which ES specificity can be achieved in cereal SSP promoters, and also suggest that there is redundancy between transcriptional and post-transcriptional tissue specificity mechanisms in cereal globulin genes.

Development of Gene Expression Markers of Acute Heat-Light Stress in Reef-Building Corals of the Genus Porites

Coral reefs are declining worldwide due to increased incidence of climate-induced coral bleaching, which will have widespread biodiversity and economic impacts. A simple method to measure the sub-bleaching level of heat-light stress experienced by corals would greatly inform reef management practices by making it possible to assess the distribution of bleaching risks among individual reef sites. Gene expression analysis based on quantitative PCR (qPCR) can be used as a diagnostic tool to determine coral condition in situ. We evaluated the expression of 13 candidate genes during heat-light stress in a common Caribbean coral Porites astreoides, and observed strong and consistent changes in gene expression in two independent experiments. Furthermore, we found that the apparent return to baseline expression levels during a recovery phase was rapid, despite visible signs of colony bleaching. We show that the response to acute heat-light stress in P. astreoides can be monitored by measuring the difference in expression of only two genes: Hsp16 and actin. We demonstrate that this assay discriminates between corals sampled from two field sites experiencing different temperatures. We also show that the assay is applicable to an Indo-Pacific congener, P. lobata, and therefore could potentially be used to diagnose acute heat-light stress on coral reefs worldwide.

Pax3 expression enhances PDGF-B-induced brainstem gliomagenesis and characterizes a subset of brainstem glioma.

High-grade Brainstem Glioma (BSG), also known as Diffuse Intrinsic Pontine Glioma (DIPG), is an incurable pediatric brain cancer. Increasing evidence supports the existence of regional differences in gliomagenesis such that BSG is considered a distinct disease from glioma of the cerebral cortex (CG). In an effort to elucidate unique characteristics of BSG, we conducted expression analysis of mouse PDGF-B-driven BSG and CG initiated in Nestin progenitor cells and identified a short list of expression changes specific to the brainstem gliomagenesis process, including abnormal upregulation of paired box 3 (Pax3). In the neonatal mouse brain, Pax3 expression marks a subset of brainstem progenitor cells, while it is absent from the cerebral cortex, mirroring its regional expression in glioma. Ectopic expression of Pax3 in normal brainstem progenitors in vitro shows that Pax3 inhibits apoptosis. Pax3-induced inhibition of apoptosis is p53-dependent, however, and in the absence of p53, Pax3 promotes proliferation of brainstem progenitors. In vivo, Pax3 enhances PDGF-B-driven gliomagenesis by shortening tumor latency and increasing tumor penetrance and grade, in a region-specific manner, while loss of Pax3 function extends survival of PDGF-B-driven;p53-deficient BSG-bearing mice by 33%. Importantly, Pax3 is regionally expressed in human glioma as well, with high PAX3 mRNA characterizing 40% of human BSG, revealing a subset of tumors that significantly associates with PDGFRA alterations, amplifications of cell cycle regulatory genes, and is exclusive of ACVR1 mutations. Collectively, these data suggest that regional Pax3 expression not only marks a novel subset of BSG but also contributes to PDGF-B-induced brainstem gliomagenesis.