Transcriptome-wide RNA Sequencing Research Articles

Abstract B01: Mechanisms of treatment resistance following Ras targeted therapy in acute myeloid leukemia.

Abstract Ras proto-oncogenes and numerous effectors of Ras signaling are recurrently mutated in acute myeloid leukemia (AML). While Ras represents an attractive therapeutic target, attempts at targeting oncogenic Ras directly have not been successful. To evaluate potential mechanisms of resistance to Ras targeted treatment in AML, our group leveraged a robust genetically-engineered mouse model in which AML is driven by expression of a tetracycline-repressible, constitutively active NRAS oncogene, tre-NRAS(V12), and a Mll-AF9 “knock-in” leukemogenic fusion gene (tre-NRAS(V12) & Mll-AF9 or tNM AML). The leukemia cells in this model are “addicted” to NRAS(V12), and mimicking Ras targeted treatment by doxycycline (Dox)-mediated suppression of the tre-NRAS(V12) oncogene expression results in rapid eradication of AML cells in vitro and in vivo. Furthermore, prolonged suppression of NRAS(V12) expression in NRAS(V12)-dependent (NRD) AML leads to the spontaneous development of NRAS(V12)-independent (NRI), “relapsed/refractory”, AML in about half of Dox treated mice. To elucidate the mechanism(s) that lead to the development of “relapsed/refractory” NRI AML following suppression of NRAS(V12) we performed transcriptome-wide next generation RNA-sequencing and comprehensive flow cytometric analysis of cancer signaling pathways comparing the parental NRD AML and two NRI AMLs. We confirmed that the NRAS(V12) oncogene is not aberrantly re-expressed. Furthermore, endogenous Ras gene expression is not up-regulated, and there is no evidence of reactivation of canonical Ras effector signaling pathways in either of the NRI AMLs by phospho-flow cytometry. Our preliminary analysis of cancer signaling pathways and transcriptome-wide RNA sequencing have identified candidate mediators of NRAS(V12)-independent AML growth and survival including the Myb proto-oncogene and anti-apoptotic Bcl2 that are enriched in both NRI AMLs relative to NRD AML at both the transcript and protein level. Furthermore, Myc protein is enriched in both NRI AMLs compared to NRD AML. We are currently performing a more comprehensive analysis of our next-generation RNA sequencing data to refine our list of candidate genes, and investigating the potential functional roles of Bcl2, Myb, and Myc in the development of NRI AML. In our initial functional studies, inhibition of Bcl2 activity suppresses NRI AML leukemic colony formation in vitro, and we are currently evaluating the ability of enforced expression of Bcl2, Myb, or Myc in NRD AML cells to render them NRAS(V12)-independent. To further investigate the translational potential of our findings, we are also evaluating the anti-leukemic effect of clinically relevant inhibitors of Bcl2 with inhibitors of canonical Ras effector pathways (RAF-MEK-ERK and PI3K-AKT-mTOR) in preclinical AML models. In this way we hope to gain a better understanding of mechanisms of treatment resistance to Ras targeted therapies, and thereby provide a foundation for the rational development of novel targeted treatment approaches for AML. Citation Format: Craig E. Eckfeldt, Robin DW Lee, Emily J. Pomeroy, Alpay N. Temiz, Susan K. Rathe, Jing Ma, Tanja A. Gruber, Ernesto Diaz-Flores, James R. Downing, Kevin M. Shannon, David A. Largaespada. Mechanisms of treatment resistance following Ras targeted therapy in acute myeloid leukemia. [abstract]. In: Proceedings of the AACR Special Conference on Hematologic Malignancies: Translating Discoveries to Novel Therapies; Sep 20-23, 2014; Philadelphia, PA. Philadelphia (PA): AACR; Clin Cancer Res 2015;21(17 Suppl):Abstract nr B01.

Exercise-induced differential changes in gene expression among arterioles of skeletal muscles of obese rats.

Using next-generation, transcriptome-wide RNA sequencing (RNA-Seq) technology we assessed the effects of exercise training on transcriptional profiles in skeletal muscle arterioles isolated from the soleus and gastrocnemius muscles of Otsuka Long Evans Tokushima Fatty (OLETF) rats that underwent an endurance exercise training program (EX; n = 13), interval sprint training program (SPRINT; n = 14), or remained sedentary (Sed; n = 12). We hypothesized that the greatest effects of exercise would be in the gastrocnemius arterioles. Results show that EX caused the largest number of changes in gene expression in the soleus and white gastrocnemius 2a arterioles with little to no changes in the feed arteries. In contrast, SPRINT caused substantial changes in gene expression in the feed arteries. IPA canonical pathway analysis revealed 18 pathways with significant changes in gene expression when analyzed across vessels and revealed that EX induces increased expression of the following genes in all arterioles examined: Shc1, desert hedgehog protein (Dhh), adenylate cyclase 4 (Adcy4), G protein binding protein, alpha (Gnat1), and Bcl2l1 and decreased expression of ubiquitin D (Ubd) and cAMP response element modulator (Crem). EX increased expression of endothelin converting enzyme (Ece1), Hsp90b, Fkbp5, and Cdcl4b in four of five arterioles. SPRINT had effects on expression of Crem, Dhh, Bcl2l1, and Ubd that were similar to EX. SPRINT also increased expression of Nfkbia, Hspa5, Tubb 2a and Tubb 2b, and Fkbp5 in all five arterioles and increased expression of Gnat1 in all but the soleus second-order arterioles. Many contractile and/or structural protein genes were increased by SPRINT in the gastrocnemius feed artery, but the same genes exhibited decreased expression in red gastrocnemius arterioles. We conclude that training-induced changes in arteriolar gene expression patterns differ by muscle fiber type composition and along the arteriolar tree.

Exercise training causes differential changes in gene expression in diaphragm arteries and 2A arterioles of obese rats.

We employed next-generation, transcriptome-wide RNA sequencing (RNA-Seq) technology to assess the effects of two different exercise training protocols on transcriptional profiles in diaphragm second-order arterioles (D2a) and in the diaphragm feed artery (DFA) from Otsuka Long Evans Tokushima Fatty (OLETF) rats. Arterioles were isolated from the diaphragm of OLETF rats that underwent an endurance exercise training program (EX; n = 13), interval sprint training program (SPRINT; n = 14), or remained sedentary (Sed; n = 12). Our hypothesis was that exercise training would have similar effects on gene expression in the diaphragm and soleus muscle arterioles because diaphragm blood flow increases during exercise to a similar extent as in soleus. Results reveal that several canonical pathways that were significantly altered by exercise in limb skeletal muscles were not among the pathways significantly changed in the diaphragm arterioles including actin cytoskeleton signaling, role of NFAT in regulation of immune response, protein kinase A signaling, and protein ubiquitination pathway. EX training altered the expression of a smaller number of genes than did SPRINT in the DFA but induced a larger number of genes with altered expression in the D2a than did SPRINT. In fact, FDR differential expression analysis (FDR, 10%) indicated that only two genes exhibited altered expression in D2a of SPRINT rats. Very few of the genes that exhibited altered expression in the DFA or D2a were also altered in limb muscle arterioles. Finally, results indicate that the 2a arterioles of soleus muscle (S2a) from endurance-trained animals and the DFA of SPRINT animals exhibited the largest number of genes with altered expression.

The C-terminal region of the RNA helicase CshA is required for the interaction with the degradosome and turnover of bulk RNA in the opportunistic pathogen Staphylococcus aureus

Staphylococcus aureus is a versatile opportunistic pathogen that adapts readily to a variety of different growth conditions. This adaptation requires a rapid regulation of gene expression including the control of mRNA abundance. The CshA DEAD-box RNA helicase was previously shown to be required for efficient turnover of the agr quorum sensing mRNA. Here we show by transcriptome-wide RNA sequencing and microarray analyses that CshA is required for the degradation of bulk mRNA. Moreover a subset of mRNAs is significantly stabilised in absence of CshA. Deletion of the C-terminal extension affects RNA turnover similar to the full deletion of the cshA gene. In accordance with RNA decay data, the C-terminal region of CshA is required for an RNA-independent interaction with components of the RNA degradation machinery. The C-terminal truncation of CshA reduces its ATPase activity and this reduction cannot be compensated at high RNA concentrations. Finally, the deletion of the C-terminal extension does affect growth at low temperatures, but to a significantly lesser degree than the full deletion, indicating that the core of the helicase can assume a partial function and opening the possibility that CshA is involved in different cellular processes.

EWS/FLI utilizes NKX2-2 to repress mesenchymal features of Ewing sarcoma.

In Ewing sarcoma, NKX2-2 is a critical activated target of the oncogenic transcription factor EWS/FLI that is required for transformation. However, its biological function in this malignancy is unknown. Here we provide evidence that NKX2-2 mediates the EWS/FLI-controlled block of mesenchymal features. Transcriptome-wide RNA sequencing revealed that NKX2-2 represses cell adhesion and extracellular matrix organization genes. NKX2-2-depleted cells form more focal adhesions and organized actin stress fibers, and spread over a wider area-hallmarks of mesenchymally derived cells. Furthermore, NKX2-2 represses the actin-stabilizing protein zyxin, suggesting that these morphological changes are attributable to zyxin de-repression. In addition, NKX2-2-knockdown cells display marked increases in migration and substrate adhesion. However, only part of the EWS/FLI phenotype is NKX2-2-dependent; consequently, NKX2-2 is insufficient to rescue EWS/FLI repression of mesenchymalization. Strikingly, we found that EWS/FLI-and NKX22-repressed genes are activated by ZEB2, which was previously shown to block Ewing sarcoma epithelialization. Together, these data support an emerging theme wherein Ewing sarcoma cells highly express transcription factors that maintain an undifferentiated state. Importantly, co-opting epithelial and mesenchymal traits by Ewing sarcoma cells may explain how the primary tumor grows rapidly while also "passively" metastasizing, without the need for transitions toward differentiated states, as in carcinomas.

Applied equine genetics.

Genome sequencing of the domestic horse and subsequent advancements in the field of equine genomics have led to an explosion in the development of tools for mapping traits and diseases and evaluating gene expression. The objective of this review is to discuss the current progress in the field of equine genomics, with specific emphasis on assembly and analysis of the reference sequence and subsequent sequencing of a Quarter Horse mare; the genomic tools currently available to researchers and their implications in genomic investigations in the horse; the genomics of Mendelian and non-Mendelian traits; the genomics of performance traits and considerations regarding genetic testing in the horse. The whole-genome sequencing of a Quarter Horse mare has provided additional variants within the equine genome that extend past single nucleotide polymorphisms to include insertions/deletions and copy number variants. Equine single nucleotide polymorphism arrays have allowed for the investigation of both simple and complex genetic traits while DNA microarrays have provided a tool for examining gene expression across various tissues and with certain disease conditions. Recently, next-generation sequencing has become more affordable and both whole-genome DNA sequencing and transcriptome-wide RNA sequencing are methodologies that are being applied to equine genomic research. Research in the field of equine genomics continues to expand rapidly as the cost of genotyping and sequencing decreases, resulting in a need for quality bioinformatics software and expertise to appropriately handle both the size and complexity of these data.

Transcriptome-wide RNA sequencing analysis of rat skeletal muscle feed arteries. I. Impact of obesity.

We employed next-generation RNA sequencing (RNA-Seq) technology to determine the influence of obesity on global gene expression in skeletal muscle feed arteries. Transcriptional profiles of the gastrocnemius and soleus muscle feed arteries (GFA and SFA, respectively) and aortic endothelial cell-enriched samples from obese Otsuka Long-Evans Tokushima Fatty (OLETF) and lean Long-Evans Tokushima Otsuka (LETO) rats were examined. Obesity produced 282 upregulated and 133 downregulated genes in SFA and 163 upregulated and 77 downregulated genes in GFA [false discovery rate (FDR) < 10%] with an overlap of 93 genes between the arteries. In LETO rats, there were 89 upregulated and 114 downregulated genes in the GFA compared with the SFA. There were 244 upregulated and 275 downregulated genes in OLETF rats (FDR < 10%) in the GFA compared with the SFA, with an overlap of 76 differentially expressed genes common to both LETO and OLETF rats in both the GFA and SFA. A total of 396 transcripts were found to be differentially expressed between LETO and OLETF in aortic endothelial cell-enriched samples. Overall, we found 1) the existence of heterogeneity in the transcriptional profile of the SFA and GFA within healthy LETO rats, 2) that this between-vessel heterogeneity was markedly exacerbated in the hyperphagic, obese OLETF rat, and 3) a greater number of genes whose expression was altered by obesity in the SFA compared with the GFA. Also, results indicate that in OLETF rats the GFA takes on a relatively more proatherogenic phenotype compared with the SFA.

Transcriptome-wide RNA sequencing analysis of rat skeletal muscle feed arteries. II. Impact of exercise training in obesity

We employed next-generation RNA sequencing (RNA-Seq) technology to determine the extent to which exercise training alters global gene expression in skeletal muscle feed arteries and aortic endothelial cells of obese Otsuka Long-Evans Tokushima Fatty (OLETF) rats. Transcriptional profiles of the soleus and gastrocnemius muscle feed arteries (SFA and GFA, respectively) and aortic endothelial cell-enriched samples from rats that underwent an endurance exercise training program (EndEx; n = 12) or a interval sprint training program (IST; n = 12) or remained sedentary (Sed; n = 12) were examined. In response to EndEx, there were 39 upregulated (e.g., MANF) and 20 downregulated (e.g., ALOX15) genes in SFA and 1 upregulated (i.e., Wisp2) and 1 downregulated (i.e., Crem) gene in GFA [false discovery rate (FDR) < 10%]. In response to IST, there were 305 upregulated (e.g., MANF, HSPA12B) and 324 downregulated genes in SFA and 101 upregulated and 66 downregulated genes in GFA, with an overlap of 32 genes between arteries. Furthermore, in aortic endothelial cells, there were 183 upregulated (e.g., eNOS, SOD-3) and 141 downregulated (e.g., ATF3, Clec1b, npy, leptin) genes with EndEx and 71 upregulated and 69 downregulated genes with IST, with an overlap of 35 between exercise programs. Expression of only two genes (Tubb2b and Slc9a3r2) was altered (i.e., increased) by exercise in all three arteries. The finding that both EndEx and IST produced greater transcriptional changes in the SFA compared with the GFA is intriguing when considering the fact that treadmill bouts of exercise are associated with greater relative increases in blood flow to the gastrocnemius muscle compared with the soleus muscle.

Transcriptome-wide sequencing reveals numerous APOBEC1 mRNA-editing targets in transcript 3′ UTRs

Apolipoprotein B-editing enzyme, catalytic polypeptide-1 (APOBEC1) is a cytidine deaminase, initially identified by its activity in converting a specific cytidine (C) to uridine (U) in apolipoprotein B (apoB) mRNA transcripts in the small intestine. Editing results in translation of a truncated apoB isoform with distinct functions in lipid transport. To address the possibility that APOBEC1 edits additional mRNAs, we developed a transcriptome-wide comparative RNA-Seq screen. We identified and validated 32 previously undescribed mRNA targets of APOBEC1 editing, all of which are located in AU-rich segments of transcript 3′ untranslated regions (3′ UTRs). Further analysis established several characteristic sequence features of editing targets, which were predictive for the identification of additional APOBEC1 substrates. The transcriptomics approach to RNA editing presented here dramatically expands the list of APOBEC1 mRNA editing targets and reveals a novel cellular mechanism for the modification of transcript 3′ UTRs.