A Dynamic Alternative Splicing Program Regulates Gene Expression In A Differentiation Stage-Specific Manner During Terminal Erythropoiesis

Blood ◽  
2013 ◽  
Vol 122 (21) ◽  
pp. 3413-3413 ◽  
Author(s):  
Harold Pimentel ◽  
Marilyn Parra ◽  
Jie Li ◽  
Sherry Gee ◽  
Dana Ghanem ◽  
...  

Abstract Spatio-temporal regulation of switches in alternative pre-mRNA splicing modulate exon usage, critically remodeling the transcriptome during development and differentiation of many tissues, while aberrant regulation of alternative splicing disrupts these processes and plays a role in numerous human diseases. Recently, the discovery of splicing factor mutations in myelodysplasia has increased interest in splicing regulation in hematology. Previously, a functionally critical erythroid splicing switch in protein 4.1 pre-mRNA has been reported, in which activation of alternative exon 16 splicing in late erythroblasts is required for assembly of a mechanically stable red cell membrane. To explore globally the landscape of important alternative splicing events in the erythroid lineage, we applied RNA-seq analysis to five highly FACS-purified populations of human erythroblasts, cultured from CD34+ cord blood progenitors, representing proerythroblasts, early and late basophilic erythroblasts, polychromatophilic erythroblasts, and orthochromatophilic erythroblasts. Alternative splicing events predicted by computational analysis were filtered to remove low expression genes and low frequency splicing events, to derive a list of >3000 ‘major’ alternative splicing events of potential importance in erythroid biology. Many of these were validated by inspection of RNA-seq reads mapped on the human genome, and/or by RT-PCR analysis. In this unique differentiation system we found an extensive and dynamic alternative splicing program enriched in genes that function in cell cycle regulation, organelle organization, chromatin structure and function, and RNA processing. For example, we identified alternative splicing events in ∼25 genes encoding chromatin modifying enzymes that methylate, demethylate, or acetylate specific lysine or arginine residues in histones; in transcription modulators such as ATRX and BCL11A that regulate normal globin gene expression; and in ∼50 RNA binding proteins with various roles in post-transcriptional gene regulation. Comparison of PSI (percent spliced in) values across the differentiation series revealed that dozens of alternative exons exhibit substantial switches in splicing efficiency during terminal erythropoiesis. The majority of splicing switches occur in late-stage polychromatophilic and orthochromatophilic erythroblasts, temporally correlated with changes in transcript abundance for many splicing factors and with substantial cell remodeling prior to enucleation. One of the biggest switches in late erythroblasts involves inclusion of a 35nt exon in the NDEL1 (nuclear distribution factor E-homolog-like1) gene, which alters C-terminal structure of a protein that functions in nuclear migration and nucleokinesis in nonerythroid cells and may have a role in erythroblast enucleation. Most of the regulated splicing events insert or delete sequences predicted to modulate protein structure and function in late erythroblasts. However, a subset of altered splicing events have a different effect on gene expression by introducing premature translation termination codons (PTCs), leading us to hypothesize that alternative splicing-coupled nonsense-mediated-decay (AS-NMD) contributes to stage-specific down-regulation of numerous erythroid transcripts. Consistent with such a model, most genes that up-regulate PTC exons in late erythroblasts exhibit reduction in overall expression levels, and inhibition of NMD increases the apparent expression of PTC isoforms. In contrast, genes that up-regulate coding exons are not preferentially down-regulated in late erythroblasts. We conclude that a dynamically regulated alternative splicing program in terminally differentiating erythroblasts plays a major post-transcriptional role in shaping gene expression as the cells transition from proliferation to differentiation, ensuring synthesis of the appropriate constellation of proteins as the cells prepare for enucleation and production of mature red cells. Disclosures: No relevant conflicts of interest to declare.

Blood ◽  
2014 ◽  
Vol 124 (21) ◽  
pp. 638-638 ◽  
Author(s):  
Naim Rashid ◽  
Stephane Minvielle ◽  
Florence Magrangeas ◽  
Mehmet Kemal Samur ◽  
Alice Clynen ◽  
...  

Abstract Alternative splicing is an important post-translational change that alters gene function. Misregulation of alternative splicing has been implicated in number of disease processes including cancer. Here we have analyzed alternative splicing in myeloma using high throughput RNA-seq. Our analytic pipeline for RNA-seq data used in this investigation not only provides information on expression levels for genes, but also provides information on the expression of known splice variants of genes (isoforms), and can identify novel exon level events across individuals (i.e. exon skipping events). We conducted a study of 328 newly-diagnosed patients with multiple myeloma treated homogeneously with novel agent combination containting lenalidomide, bortezomib and dexamethsone with or without high-dose melphalan followed by lenalidomide maintenance in the IFM/DFCI study. RNA isolated from purified CD138+ MM cells collected at the time of diagnosis and from 18 normal donor plasma cells were processed by RNA-seq (100 million paired end reads on Illumina HiSeq) and analyzed using a custom computational and statistical pipeline. Following read alignment to hg19, we utilized RSEM to quantify both gene-level and isoform-level expression of known ENSEMBL transcripts. We then implemented a novel testing approach based on compositional regression to discover genes that show significant isoform switching between the 328 MM samples and 18 Normal Plasma Cell (NPC) samples from healthy donors. Using various programs and their modifications, we also identified novel alternative splicing events, such as exon skipping and mutually exclusive exon usage, among others. Patient data for MM characteristics, cytogenetic and FISH as well as clinical survival outcomes were also analyzed and correlated with genomic data. We observed over 600 genes showing significant changes in relative isoform abundances (isoform switching) between MM and normal samples. A number of previously characterized genes including MYCL1 (adj. p = 0.0014) and CCND3 (adj. p = 0.0013), and MAP kinase-related genes (MAP3K8, MAPKAPK2, MAPKAPK3, MAP4K4) exhibited significant isoform switching compared to normal, in addition to some not well characterized genes. Genes showing the greatest magnitude of isoform switching include MEFV (adj. p = 2.7 x 10-5), showing a two fold change in the relative major isoform abundance compared to normal, and has been previously shown to have a role in lymphoid neoplasms. We applied hierarchical clustering to the isoforms showing significant changes in isoform-switching and identified 4 distinct clusters, which are currently being investigated for correlation with clinical subtypes of MM. Exon level analyses of alternative splicing events, such as exon skipping, are currently underway. Clinical data including MM characteristics, cytogenetics, FISH and survival outcomes was available for a subset of 265 patients. We found that 109 genes showed significant isoform switching between t(4;14) and non-t(4;14) patients, such as CD44 (adj. p =1.8 x 10-6) and WHSC1 (adj. p =5.1 x 10-28). Comparing del17p (28 in total) and non del17p patients, we found no significant splicing changes after multiple testing adjustment. Of these genes, only a subset (40%) were shown to be differentially expressed in terms of total gene expression, suggesting the importance of examining alternative splicing events in addition to total gene expression. With respect to treatment response, we compared the expression of gene isoforms between patients achieving complete response (CR) versus others and identified 38 isoforms associated with response to treatment (adj. p value < 0.05), with SEPT9, SLC2A5, and UBX6 having the strongest associations (adj. p-value < 3 x 10-4). Using a univariate cox regression model, 4 spliced isoforms relating to 3 genes were identified as having significant correlation with event-free survival (EFS) (FDR-adjusted cox p value < 0.05). We are in the process of now integrating the gene expression data with altered splicing data to develop an integrated survival model. In summary, this study highlights the significant frequency, biological and clinical importance of alternative splicing in MM and points to the need for evaluation of not only the expression level of genes but also post-translational modifications. The genes identified here are important targets for therapy as well as possible immune modulation. Disclosures Moreau: Celgene Corporation: Honoraria, Membership on an entity's Board of Directors or advisory committees.


Blood ◽  
2012 ◽  
Vol 120 (21) ◽  
pp. 978-978
Author(s):  
Harold Pimentel ◽  
Jie Li ◽  
Marilyn Parra ◽  
Dana Ghanem ◽  
Sherry Gee ◽  
...  

Abstract Abstract 978 Alternative pre-mRNA splicing plays a major role in development and differentiation by re-modeling the transcriptome to generate mRNAs that encode the biologically appropriate cell type-specific proteome. Earlier studies employing RT-PCR and exon microarrays demonstrated a small number of splicing changes during erythroid differentiation, one of which (in protein 4.1R) is critical for mechanical stability of the membrane skeleton. Here we report that the landscape of splicing changes executed during terminal erythropoiesis is far more extensive and highly dynamic, ultimately affecting the expression of many more proteins than previously recognized. Highly purified populations of FACS-sorted cells representing erythroblasts at distinct stages during terminal erythroid differentiation from proerythroblasts to orthochromatic erythroblasts of both human and mouse origin were used as the source of RNA for RNA-seq analysis. In total, hundreds of millions of sequence reads were obtained from three biological replicates for four (mouse) or five (human) cell populations, and reads were aligned to the Ensembl-annotated transcriptome using the Bowtie aligner. Transcript-level estimates were obtained using the streaming transcript abundance estimation tool, eXpress, expression of individual exons in “exon-inclusion” isoforms relative to total isoforms was represented as Ψ (psi), or percent spliced in, and statistical significance estimates adjusted for multiple comparisons by the Benjamini-Hochberg method. Thousands of alternative splicing events were predicted in genes with diverse functions in transcription, RNA processing, protein synthesis, membrane receptors, cytoskeletal structure, etc. Initial RT-PCR studies indicate that a high proportion of predicted alternative splicing events can be validated. Comparison of Ψ values across the differentiation series revealed that hundreds of alternative exons in erythroid transcripts exhibit substantial differences in splicing efficiency between proerythroblasts and orthrochromatic erythroblasts (ΔΨ>20%), suggesting that their splicing efficiency is regulated. Both increases and decreases in exon splicing efficiency were observed, indicating that multiple splicing regulatory pathways are active and that both splicing enhancer and splicing silencer factors are involved in the regulation. Interestingly, some of the splicing switches introduce premature translation termination codons, leading us to hypothesize that splicing-coupled nonsense mediated decay may down-regulate expression of a class of erythroid transcripts. To begin exploring mechanisms that regulate the late erythroid alternative splicing program, we used the RNA-seq data to derive differentiation stage-specific expression profiles of known splicing factors. Major changes in the expression profile of many splicing regulators were observed. hnRNPA1 was strongly down-regulated in late erythroblasts, in concert with up-regulation of the protein 4.1R splicing switch it has been shown to inhibit. Because many other exons are up-regulated with similar kinetics, hnRNPA1 may be a general inhibitor of alternative splicing in early erythroblasts. In contrast, RNA-seq data indicate that several other splicing factors including MBNL1, a known splicing factor in muscle and brain whose activity is disturbed in myotonic dystrophy, are substantially up-regulated in late erythropoiesis. We conclude that a highly dynamic alternative splicing program in terminally differentiating erythroblasts, in conjunction with the better studied transcriptional program, plays a major role in regulating gene expression to insure synthesis of the appropriate constellation of proteins both quantitatively and qualitatively as the cells are remodeled in preparation for production of mature red cells. Disclosures: No relevant conflicts of interest to declare.


Blood ◽  
2014 ◽  
Vol 124 (21) ◽  
pp. 449-449
Author(s):  
Harold Pimentel ◽  
Marilyn Parra ◽  
Sherry Gee ◽  
Narla Mohandas ◽  
Lior Pachter ◽  
...  

Abstract Erythroid RNAs, like their nonerythroid counterparts, are subject to post-transcriptional processing events that critically impact their coding capacity for the erythroid proteome. Previous studies have shown that differentiating human and mouse erythroblasts execute an extensive and dynamic alternative splicing program involving regulation of numerous alternative exons. Here we report that controlled excision of selected introns is also an important component of the erythroblast alternative splicing program. Intron retention (IR) patterns in differentiating human erythroblasts were determined via RNA-seq analysis of FACS-purified erythroblast populations. Comparison of IR among erythroblast populations and between erythroblasts and other hematopoietic cells suggests that regulation of IR occurs in a differentiation stage- and tissue-specific manner. For example, there was little overlap of intron retention events in erythroblasts with those reported in differentiating granulocytes. Moreover, the IR profile of proerythroblasts differed substantially from that in orthochromatic erythroblasts, with IR generally increasing in the more mature cells that are preparing for enucleation. IR in erythroblasts affected numerous genes functioning in RNA processing, iron homeostasis and heme biosynthesis, protein translation, and membrane properties. Mature erythroblasts exhibited retention of introns in several human disease genes including SF3B1, a splicing factor often mutated in myelodysplasia; TFR2, encoding transferrin receptor 2 that is mutated in a form of hemochromatosis; and FUS, an RNA binding protein implicated in ALS. Inspection of intronic RNA-seq reads in >60 genes with IR revealed that single or multiple introns can be retained within a transcript; however, other introns within the same genes, and indeed the great majority of introns in erythroblast-expressed genes, are efficiently spliced with minimal or no IR. Retained introns may be flanked by either constitutively or alternatively spliced exons, suggesting different regulatory mechanisms. Ongoing studies will explore whether IR in some transcripts might function to down-regulate gene expression by introduction of premature termination codons that induce nonsense-mediated decay, or alternatively, whether IR transcripts could represent a reserve of nearly-completed mRNAs that can be processed in response to appropriate physiological stimuli. In sum, these results suggest that a highly regulated IR program plays an important role in erythroid differentiation. Disclosures No relevant conflicts of interest to declare.


2021 ◽  
pp. 1-6
Author(s):  
Miriam C. Aziz ◽  
Patricia N. Schneider ◽  
Gemma L. Carvill

Developmental and epileptic encephalopathies (DEEs) describe a subset of neurodevelopmental disorders categorized by refractory epilepsy that is often associated with intellectual disability and autism spectrum disorder. The majority of DEEs are now known to have a genetic basis with de novo coding variants accounting for the majority of cases. More recently, a small number of individuals have been identified with intronic <i>SCN1A</i> variants that result in alternative splicing events that lead to ectopic inclusion of poison exons (PEs). PEs are short highly conserved exons that contain a premature truncation codon, and when spliced into the transcript, lead to premature truncation and subsequent degradation by nonsense-mediated decay. The reason for the inclusion/exclusion of these PEs is not entirely clear, but research suggests an autoregulatory role in gene expression and protein abundance. This is seen in proteins such as RNA-binding proteins and serine/arginine-rich proteins. Recent studies have focused on targeting these PEs as a method for therapeutic intervention. Targeting PEs using antisense oligonucleotides (ASOs) has shown to be effective in modulating alternative splicing events by decreasing the amount of transcripts harboring PEs, thus increasing the abundance of full-length transcripts and thereby the amount of protein in haploinsufficient genes implicated in DEE. In the age of personalized medicine, cellular and animal models of the genetic epilepsies have become essential in developing and testing novel precision therapeutics, including PE-targeting ASOs in a subset of DEEs.


2004 ◽  
Vol 18 (2) ◽  
pp. 167-183 ◽  
Author(s):  
Jianhua Zhang ◽  
Amy Moseley ◽  
Anil G. Jegga ◽  
Ashima Gupta ◽  
David P. Witte ◽  
...  

To understand the commitment of the genome to nervous system differentiation and function, we sought to compare nervous system gene expression to that of a wide variety of other tissues by gene expression database construction and mining. Gene expression profiles of 10 different adult nervous tissues were compared with that of 72 other tissues. Using ANOVA, we identified 1,361 genes whose expression was higher in the nervous system than other organs and, separately, 600 genes whose expression was at least threefold higher in one or more regions of the nervous system compared with their median expression across all organs. Of the 600 genes, 381 overlapped with the 1,361-gene list. Limited in situ gene expression analysis confirmed that identified genes did represent nervous system-enriched gene expression, and we therefore sought to evaluate the validity and significance of these top-ranked nervous system genes using known gene literature and gene ontology categorization criteria. Diverse functional categories were present in the 381 genes, including genes involved in intracellular signaling, cytoskeleton structure and function, enzymes, RNA metabolism and transcription, membrane proteins, as well as cell differentiation, death, proliferation, and division. We searched existing public sites and identified 110 known genes related to mental retardation, neurological disease, and neurodegeneration. Twenty-one of the 381 genes were within the 110-gene list, compared with a random expectation of 5. This suggests that the 381 genes provide a candidate set for further analyses in neurological and psychiatric disease studies and that as a field, we are as yet, far from a large-scale understanding of the genes that are critical for nervous system structure and function. Together, our data indicate the power of profiling an individual biologic system in a multisystem context to gain insight into the genomic basis of its structure and function.


2016 ◽  
Vol 14 (1) ◽  
pp. nrs.14002 ◽  
Author(s):  
Shailaja D. Divekar ◽  
Deanna M. Tiek ◽  
Aileen Fernandez ◽  
Rebecca B. Riggins

Estrogen-related receptors (ERRs) are founding members of the orphan nuclear receptor (ONR) subgroup of the nuclear receptor superfamily. Twenty-seven years of study have yet to identify cognate ligands for the ERRs, though they have firmly placed ERRα (ESRRA) and ERRγ (ESRRG) at the intersection of cellular metabolism and oncogenesis. The pace of discovery for novel functions of ERRβ (ESRRB), however, has until recently been somewhat slower than that of its family members. ERRβ has also been largely ignored in summaries and perspectives of the ONR literature. Here, we provide an overview of established and emerging knowledge of ERRβ in mouse, man, and other species, highlighting unique aspects of ERRβ biology that set it apart from the other two estrogen-related receptors, with a focus on the impact of alternative splicing on the structure and function of this receptor.


Blood ◽  
2019 ◽  
Vol 134 (Supplement_1) ◽  
pp. 812-812
Author(s):  
Mudit Chaand ◽  
Chris Fiore ◽  
Brian T Johnston ◽  
Diane H Moon ◽  
John P Carulli ◽  
...  

Human beta-like globin gene expression is developmentally regulated. Erythroblasts (EBs) derived from fetal tissues, such as umbilical cord blood (CB), primarily express gamma globin mRNA (HBG) and HbF, while EBs derived from adult tissues, such as bone marrow (BM), predominantly express beta globin mRNA (HBB) and adult hemoglobin. Human genetics has validated de-repression of HBG in adult EBs as a powerful therapeutic paradigm in diseases involving defective HBB, such as sickle cell anemia. To identify novel factors involved in the switch from HBG to HBB expression, and to better understand the global regulatory networks driving the fetal and adult cell states, we performed transcriptome profiling (RNA-seq) and chromatin accessibility profiling (ATAC-seq) on sorted EB cell populations from CB or BM. This approach improves upon previous studies that used unsorted cells (Huang J, Dev Cell 2016) or that did not measure chromatin accessibility (Yan H, Am J Hematol 2018). CD34+ cells from CB and BM were differentiated using a 3-phase in vitro culture system (Giarratana M, Blood 2011). Fluorescence-activated cell sorting and the cell surface markers CD36 and GYPA were used to isolate 7 discrete populations, with each sorting gate representing increasingly mature, stage-matched EBs from CB or BM (Fig 1A, B). RNA-seq analysis revealed expected expression patterns of the beta-like globins, with total levels increasing during erythroid maturation and primarily composed of HBB or HBG transcripts in BM or CB, respectively (Fig 1C). Erythroid maturation led to progressive increases in chromatin accessibility at the HBB promoter in BM populations. In CB-derived cells, erythroid maturation led to progressive increases in chromatin accessibility at the HBG promoters through the CD36+GYPA+ stage (Pops 1-5). Chromatin accessibility shifted from the HBG promoters to the HBB promoter during the final stages of differentiation (Pops 6-7), suggesting that HBG gene activation is transient in CB EBs (Fig 1D). Hierarchical clustering and principal component analysis of ATAC-seq data revealed that cell populations cluster based on differentiation stage rather than by BM or CB lineage, suggesting most molecular changes are stage-specific, not lineage-specific (Fig 2A, B). To identify transcription factors driving cell state, and potentially beta-like globin expression preference, we searched for DNA binding motifs within regions of differential chromatin accessibility and found NFI factor motifs enriched under peaks that were larger in BM relative to CB (Fig 2C). Transcription factor footprinting analysis showed that both flanking accessibility and footprint depth at NFI motifs were also increased in BM relative to CB (Fig 2D). Increased chromatin accessibility was observed at the NFIX promoter in BM relative to CB populations, and in HUDEP-2 relative to HUDEP-1 cell lines (Fig 2E). Furthermore, accessibility at the NFIX promoter correlated with elevated NFIX mRNA in BM and HUDEP-2 relative to CB and HUDEP-1, respectively. Together these data implicated NFIX in HbF repression, a finding consistent with previous genome-wide association and DNA methylation studies that suggested a possible role for NFIX in regulating beta-like globin gene expression (Fabrice D, Nat Genet 2016; Lessard S, Genome Med 2015). To directly test the hypothesis that NFIX represses HbF, short hairpin RNAs were used to knockdown (KD) NFIX in primary erythroblasts derived from human CD34+ BM cells (Fig 3A). NFIX KD led to a time-dependent induction of HBG mRNA, HbF, and F-cells comparable to KD of the known HbF repressor BCL11A (Fig 3B-D). A similar effect on HbF was observed in HUDEP-2 cells following NFIX KD (Fig 3E). Consistent with HbF induction, NFIX KD also increased chromatin accessibility and decreased DNA methylation at the HBG promoters in primary EBs (Fig 3F, G). NFIX KD led to a delay in erythroid differentiation as measured by CD36 and GYPA expression (Fig 3H). Despite this delay, by day 14 a high proportion of fully enucleated erythroblasts was observed, suggesting NFIX KD cells are capable of terminal differentiation (Fig 3H). Collectively, these data have enabled identification and validation of NFIX as a novel repressor of HbF, a finding that enhances the understanding of beta-like globin gene regulation and has potential implications in the development of therapeutics for sickle cell disease. Disclosures Chaand: Syros Pharmaceuticals: Employment, Equity Ownership. Fiore:Syros Pharmaceuticals: Employment, Equity Ownership. Johnston:Syros Pharmaceuticals: Employment, Equity Ownership. Moon:Syros Pharmaceuticals: Employment, Equity Ownership. Carulli:Syros Pharmaceuticals: Employment, Equity Ownership. Shearstone:Syros Pharmaceuticals: Employment, Equity Ownership.


2017 ◽  
Vol 63 (2) ◽  
pp. 89-99 ◽  
Author(s):  
Maria C. Davis ◽  
Christopher A. Kesthely ◽  
Emily A. Franklin ◽  
Shawn R. MacLellan

Transcription is the first and most heavily regulated step in gene expression. Sigma (σ) factors are general transcription factors that reversibly bind RNA polymerase (RNAP) and mediate transcription of all genes in bacteria. σ Factors play 3 major roles in the RNA synthesis initiation process: they (i) target RNAP holoenzyme to specific promoters, (ii) melt a region of double-stranded promoter DNA and stabilize it as a single-stranded open complex, and (iii) interact with other DNA-binding transcription factors to contribute complexity to gene expression regulation schemes. Recent structural studies have demonstrated that when σ factors bind promoter DNA, they capture 1 or more nucleotides that are flipped out of the helical DNA stack and this stabilizes the promoter open-complex intermediate that is required for the initiation of RNA synthesis. This review describes the structure and function of the σ70 family of σ proteins and the essential roles they play in the transcription process.


Sign in / Sign up

Export Citation Format

Share Document