Integrative Analysis Of RNA Sequencing Research Articles

Genome-wide RNA-seq, DNA methylation and small RNA-seq analysis unraveled complex gene regulatory networks in psoriasis pathogenesis

Psoriasis is a complex inflammatory skin disease characterized by reversible albeit relapsing red scaly plaques in the skin of a patient. In addition to the genetic predisposition, involvement of epigenetic and non-coding RNAs have also been liked with the disease. Nevertheless, any comprehensive study involving transcriptomic, small-RNA and DNA methylation at the genomic level from same patients is lacking. To investigate the complex regulation of molecular pathways in psoriasis, we carried out multi-omics integrative analysis of RNA-sequencing, small RNA-sequencing and DNA methylation profiling from the psoriatic and adjacent normal skin tissues. Our multi-omics analysis identified the genes and biological processes regulated either independently or in combination by DNA methylation and microRNAs. We identified miRNAs that specifically regulated keratinocyte hyper-proliferation, and cell cycle progression and checkpoint signaling in psoriasis. On contrary, DNA methylation was found to be more predominant in regulating immune and inflammatory responses, another causative factor in psoriasis pathogenesis. Many characteristic pathways in psoriasis e.g., Th17 cell differentiation and JAK-STAT signaling, were found to be regulated by both miRNAs and DNA methylation. We carried out functional characterization of a downregulated miRNA hsa-let-7c-5p, predicted to target upregulated genes in psoriasis involved in cell cycle processes, Th17 cell differentiation and JAK-STAT signaling pathways. Overexpression of hsa-let-7c-5p in keratinocytes caused the downregulation of its target genes, resulting in reduced cell proliferation and migration rates, demonstrating potential of miRNAs in regulating psoriasis pathogenesis. In conclusion, our findings identified distinct and shared gene-networks regulated by DNA methylation and miRNAs of a complex disease with reversible phenotype.

Novel T cell exhaustion gene signature to predict prognosis and immunotherapy response in thyroid carcinoma from integrated RNA-sequencing analysis.

Exhausted CD8+ T lymphocytes and tumor-associated macrophages play critical roles in determining cancer prognosis and the efficacy of immunotherapy. Our study revealed a negative correlation between exhausted CD8+ T lymphocytes and prognosis in thyroid carcinoma (THCA). Consensus clustering divided patients into two subgroups of exhaustion with different prognoses, as defined by marker genes of exhausted CD8+ T cells. Subsequently, we constructed an eight-gene prognostic signature, and developed a risk score named the exhaustion-related gene score (ERGS) to forecast both prognosis and immunotherapy response in THCA. Bulk RNA sequencing analysis revealed a higher prevalence of M2 macrophages, indicative of an immunosuppressive tumor microenvironment (TME), in the high-ERGS group. Single-cell RNA sequencing showed that SPP1+ macrophages and CD14+ monocytes infiltrations were positively associated with higher ERGS. Functionally, it was determined that SPP1+ macrophages exert an immunosuppressive role, while CD14+ monocytes were implicated in promoting tumor progression and angiogenesis. Analysis of cell-cell interactions between SPP1+ macrophages and T cells highlighted the activation of the SPP1-CD44 and MIF-CD74 axes, both of which could foster an immunosuppressive TME. Therapeutic strategies that target SPP1+ macrophages, CD14+ monocytes, and the SPP1-CD44 and MIF-CD74 axes may potentially improve the prognosis and amplify the immunotherapy response in THCA patients.

Abstract 4738: Overcoming therapy resistance in rhabdomyosarcoma by targeting the PI3K-MAPK pathway

Abstract Rhabdomyosarcoma (RMS) is a devastating pediatric cancer. There has not been any significant change in the treatment options for the last 4 decades, especially for relapsed or advanced disease. There is an urgent need to identify alternative treatment options to improve survival outcomes of RMS patients. RMS is comprised of two major subtypes as defined by the presence or absence of PAX3/PAX7-FOXO1 gene fusion. Fusion-negative RMS (FN RMS) is characterized by alterations in the RAS-PI3K-MAPK pathway in >90% of the cases. However, the role of the PI3K-MAPK pathway in therapy resistance and disease relapse remains poorly characterized. We have generated two FN-RMS lines, RD and SMS-CTR, with resistance to the standard-of-care chemotherapeutic agent, vincristine, through long-term culturing in incrementally increased concentrations of the drug. By RNA sequencing, vincristine-resistant RD and SMS-CTR cells demonstrated increased expression of genes in the RAS-PI3K-MAPK pathway compared to the parental lines. We hypothesize that increased activity of the PI3K-MAPK pathway contributes to therapy resistance in FN RMS. To test the hypothesis, we first showed that targeted disruption of PIK3CA by shRNA inhibited the cell growth of vincristine-resistant FN RMS cells. The same cells also showed increased sensitivity to the MEK inhibitor, trametinib, as well as the PIK3Cα inhibitor, alpelisib, compared to the parental cells. We then showed that zebrafish FN-RMS tumors treated with the combination of vincristine and alpelisib significantly inhibited tumor growth compared to the treatment with each agent alone. Through integrative analysis of RNA sequencing and ATAC-seq data, we identified MYOD1 and NFATC2 as among the top candidate upstream regulators. Targeted disruption of either MYOD1 or NFATC2 by CRISPR/Cas9 reduced cell growth of vincristine-resistant FN RMS cells and altered expression of the genes in the PI3K-MAPK pathway, implicating MYOD1 and NFATC2 in driving abnormal activity of this pathway. To assess the relevance of our findings in the clinical settings, we showed that selective components in the PI3K-MAPK pathways are up-regulated in 3 matched primary-recurrence/metastasis paraffinized tissue samples of human FN RMS by the NanoString GeoMx Digital Spatial Profiler (DSP) platform. Our findings indicate that targeting the PI3K-MAPK pathway is a promising targeted therapy option for overcoming therapy resistance in RMS. Citation Format: Yadong Wang, Eleanor Chen. Overcoming therapy resistance in rhabdomyosarcoma by targeting the PI3K-MAPK pathway [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 4738.

Integrative analysis of RNA-sequencing and microarray for the identification of adverse effects of UVB exposure on human skin.

Ultraviolet B (UVB) from sunlight represents a major environmental factor that causes toxic effects resulting in structural and functional cutaneous abnormalities in most living organisms. Although numerous studies have indicated the biological mechanisms linking UVB exposure and cutaneous manifestations, they have typically originated from a single study performed under limited conditions. We accessed all publicly accessible expression data of various skin cell types exposed to UVB, including skin biopsies, keratinocytes, and fibroblasts. We performed biological network analysis to identify the molecular mechanisms and identify genetic biomarkers. We interpreted the inflammatory response and carcinogenesis as major UVB-induced signaling alternations and identified three candidate biomarkers (IL1B, CCL2, and LIF). Moreover, we confirmed that these three biomarkers contribute to the survival probability of patients with cutaneous melanoma, the most aggressive and lethal form of skin cancer. Our findings will aid the understanding of UVB-induced cutaneous toxicity and the accompanying molecular mechanisms. In addition, the three candidate biomarkers that change molecular signals due to UVB exposure of skin might be related to the survival rate of patients with cutaneous melanoma.

Integrated RNA-Seq Analysis Uncovers the Potential Mechanism of the "Kidney Governing Bones" Theory of TCM.

Background As in philosophy of traditional Chinese medicine (TCM), the theory of “kidney governing bones” has been demonstrated by a series of scientific studies. Furthermore, many groups including ours have explored the molecular mechanisms related to bone development, growth, and regeneration using modern biology technologies, such as RNA sequencing (RNA-Seq) and isobaric tags for relative and absolute quantification (ITRAQ), and have demonstrated that the underlying molecular mechanisms were highly consistent with the “kidney governing bones” theory. Objective Kidney-yang deficiency (YD), as a pathological condition, has a passive effect on the skeleton growth; more specifically, it is a state of skeletal metabolic disorder. However, the exact molecular mechanisms related to the “kidney governing bones” theory under the control of multiple organs and systems are still unknown. Methods In this study, we performed RNA-Seq analysis to investigate and compare the gene expression patterns of six types of tissue (bone, cartilage, kidney, testicle, thyroid gland, and adrenal gland) from YD rats and normal rats and analyzed the interaction effects controlled by multiple functional genes and signaling pathways between those tissues. Results Our results showed that, in the state of YD, the functions of bone and cartilage were inhibited. Furthermore, multiple organs involving the reproductive, endocrine, and urinary systems were also investigated, and our results showed that YD could cause dysfunctions of these systems by downregulating multiple functional genes and signaling pathways that positively regulate the homeostasis of these tissues. Conclusion We ensure that “kidney governing bones” was not a simple change in a single gene but the changes in complex biological networks caused by functional changes in multiple genes. This also coincides with the holistic view of TCM, which holds that the human body itself is an organic whole and the functional activities of each organ coordinate with each other.

Integrative Analysis of Gene Expression Data by RNA Sequencing for Differential Diagnosis of Acute Leukemia: Potential Application of Machine Learning.

BCR-ABL1–positive acute leukemia can be classified into three disease categories: B-lymphoblastic leukemia (B-ALL), acute myeloid leukemia (AML), and mixed-phenotype acute leukemia (MPAL). We conducted an integrative analysis of RNA sequencing (RNA-seq) data obtained from 12 BCR-ABL1–positive B-ALL, AML, and MPAL samples to evaluate its diagnostic utility. RNA-seq facilitated the identification of all p190 BCR-ABL1 with accurate splicing sites and a new gene fusion involving MAP2K2. Most of the clinically significant mutations were also identified including single-nucleotide variations, insertions, and deletions. In addition, RNA-seq yielded differential gene expression profile according to the disease category. Therefore, we selected 368 genes differentially expressed between AML and B-ALL and developed two differential diagnosis models based on the gene expression data using 1) scoring algorithm and 2) machine learning. Both models showed an excellent diagnostic accuracy not only for our 12 BCR-ABL1–positive cases but also for 427 public gene expression datasets from acute leukemias regardless of specific genetic aberration. This is the first trial to develop models of differential diagnosis using RNA-seq, especially to evaluate the potential role of machine learning in identifying the disease category of acute leukemia. The integrative analysis of gene expression data by RNA-seq facilitates the accurate differential diagnosis of acute leukemia with successful detection of significant gene fusion and/or mutations, which warrants further investigation.

Perturbations in HNRNPH1 Splicing and Abundance Affect Global Splicing and Proliferation in Mantle Cell Lymphoma

Objectives Mantle cell lymphoma (MCL) is an uncommon B-cell non-Hodgkin lymphoma that is incurable with standard therapies. The genetic drivers of this cancer have not been firmly established and the features known to contribute to differences in clinical course remain limited. We previously discovered non-coding and silent mutations in HNRNPH1 that affect its splicing and contribute to poor outcomes for patients with MCL. We sought to extend our understanding of the mechanisms by which HNRNPH1 contributes to MCL pathology using a combination of in vitro models and integrative analysis of RNA sequencing from MCL tumors. Methods We previously sequenced ribosomal RNA-depleted RNA from 130 MCL tumors. Based on our earlier identification of mutations in HNRNPH1 and altered splicing of this gene, we performed differential splicing analyses using rMATS and leafcutter. We investigated the functional and phenotypic effect of deregulated hnRNP H1 protein through siRNA knockdown. Results Our previous work identified that splicing of HNRNPH1, and not total mRNA expression, correlated with protein abundance in MCL tumors. As a result, our analysis of alternative splicing focused on events associated with altered splicing of HNRNPH1. We identified 155 unique alternative splicing events (ΔPSI > 0.1, FDR < 0.1). Gene ontology analysis identified various aspects of RNA processing which are significantly enriched within this gene list, including mRNA splicing, transport, and metabolic process. This nominates HNRNPH1 as part of the complex network controlling alternative splicing within MCL. Available CLIP-seq in HeLa cells provides evidence for direct interactions between hnRNP H1 and transcripts identified by our analysis (e.g. RBM25, EIF4A1, HNRNPA2B1). Of the 155 events we identified, more than half involved retained introns. Generally, retained introns result in non-productive RNA species, which indicates that this program of intron retention in MCL is a mechanism by which protein abundance can be regulated by hnRNP H1. For all cases with available Mantle Cell Lymphoma International Prognostic Indicator (MIPI) classification, we determined the splicing ratio for HNRNPH1 and observed a general association between high MIPI scores and a lower ratio of non-productive HNRNPH1 transcripts. This suggested that the increased hnRNP H1 abundance we observed in HNRNPH1-mutant tumors contributes to increased proliferation of MCL cells. We verified this in vitro with siRNA knockdown of HNRNPH1 in HEK cells, which resulted in a significant decrease in cell proliferation. Conclusions We have described a pattern of alternative splicing in MCL that is associated with alterations in HNRNPH1 splicing and related protein abundance. The prevalence of retained introns suggests that hnRNP H1 regulates the abundance of protein-coding transcripts via alternative splicing coupled to nonsense-mediated decay. We continue to explore targets of hnRNP H1, a novel oncoprotein in MCL. Disclosures Morin: Celgene: Consultancy.

FRNAkenseq: a fully powered-by-CyVerse cloud integrated RNA-sequencing analysis tool.

BackgroundDecreasing costs make RNA sequencing technologies increasingly affordable for biologists. However, many researchers who can now afford sequencing lack access to resources necessary for downstream analysis. This means that even as algorithms to process RNA-Seq data improve, many biologists still struggle to manage the sheer volume of data produced by next generation sequencing (NGS) technologies. Scalable bioinformatics tools that exploit multiple platforms are needed to democratize bioinformatics resources in the sequencing era. This is essential for equipping many research groups in the life sciences with the tools to process the increasingly unwieldy datasets they produce.MethodsOne strategy to address this challenge is to develop a modern generation of sequence analysis tools capable of seamless data sharing and communication. Such tools will provide interoperability through offerings of interlinked resources. Systems of interlinked, scalable resources, which often incorporate cloud data storage, are broadly referred to as cyberinfrastructure. Cyberinfrastructure integrated tools will help researchers to robustly analyze large scale datasets by efficiently sharing data burdens across a distributed architecture. Additionally, interoperability will allow emerging tools to cross-adapt features of existing tools. It is important that these tools are designed to be easy to use for biologists.ResultsWe introduce fRNAkenseq, a powered-by-CyVerse RNA sequencing analysis tool that exhibits interoperability with other resources and meets the needs of biologists for comprehensive, easy to use RNA sequencing analysis. fRNAkenseq leverages a complex set of Application Programming Interfaces (APIs) associated with the NSF-funded cyberinfrastructure project, CyVerse, to execute FASTQ-to-differential expression RNA-Seq analyses. Integrating across bioinformatics platforms, fRNAkenseq also exploits cloud integration and cross-talk with another CyVerse associated tool, CoGe. fRNAkenseq offers novel features for the biologist such as more robust and comprehensive pipelines for enrichment than those currently available by default in a single tool, whether they are cloud-based or local installation. Importantly, cross-talk with CoGe allows fRNAkenseq users to execute RNA-Seq pipelines on an inventory of 47,000 archived genomes stored in CoGe or upload their own draft genome.

Integrative Analysis of Breast Cancer Cells Reveals an Epithelial-Mesenchymal Transition Role in Adaptation to Acidic Microenvironment.

Early ducts of breast tumors are unequivocally acidic. High rates of glycolysis combined with poor perfusion lead to a congestion of acidic metabolites in the tumor microenvironment, and pre-malignant cells must adapt to this acidosis to thrive. Adaptation to acidosis selects cancer cells that can thrive in harsh conditions and are capable of outgrowing the normal or non-adapted neighbors. This selection is usually accompanied by phenotypic change. Epithelial mesenchymal transition (EMT) is one of the most important switches correlated to malignant tumor cell phenotype and has been shown to be induced by tumor acidosis. New evidence shows that the EMT switch is not a binary system and occurs on a spectrum of transition states. During confirmation of the EMT phenotype, our results demonstrated a partial EMT phenotype in our acid-adapted cell population. Using RNA sequencing and network analysis we found 10 dysregulated network motifs in acid-adapted breast cancer cells playing a role in EMT. Our further integrative analysis of RNA sequencing and SILAC proteomics resulted in recognition of S100B and S100A6 proteins at both the RNA and protein level. Higher expression of S100B and S100A6 was validated in vitro by Immunocytochemistry. We further validated our finding both in vitro and in patients' samples by IHC analysis of Tissue Microarray (TMA). Correlation analysis of S100A6 and LAMP2b as marker of acidosis in each patient from Moffitt TMA approved the acid related role of S100A6 in breast cancer patients. Also, DCIS patients with higher expression of S100A6 showed lower survival compared to lower expression. We propose essential roles of acid adaptation in cancer cells EMT process through S100 proteins such as S100A6 that can be used as therapeutic strategy targeting both acid-adapted and malignant phenotypes.

Targeting Super-Enhancer-Driven Oncogenic Transcription by CDK7 Inhibition in Anaplastic Thyroid Carcinoma.

Background: Anaplastic thyroid carcinoma (ATC) is one of the most aggressive malignancies, with no effective treatment currently available. The molecular mechanisms of ATC carcinogenesis remain poorly understood. The objective of this study was to investigate the mechanisms and functions of super-enhancer (SE)-driven oncogenic transcriptional addiction in the progression of ATC and identify new drug targets for ATC treatments. Methods: High-throughput chemical screening was performed to identify new drugs inhibiting ATC cell growth. Cell viability assay, colony formation analysis, cell-cycle analysis, and animal study were used to examine the effects of drug treatments on ATC progression. Chromatin immunoprecipitation sequencing was conducted to establish a SE landscape of ATC. Integrative analysis of RNA sequencing, chromatin immunoprecipitation sequencing, and CRISPR/Cas9-mediated gene editing was used to identify THZ1 target genes. Drug combination analysis was performed to assess drug synergy. Patient samples were analyzed to evaluate candidate biomarkers of prognosis in ATC. Results: THZ1, a covalent inhibitor of cyclin-dependent kinase 7 (CDK7), was identified as a potent anti-ATC compound by high-throughput chemical screening. ATC cells, but not papillary thyroid carcinoma cells, are exceptionally sensitive to CDK7 inhibition. An integrative analysis of both gene expression profiles and SE features revealed that the SE-mediated oncogenic transcriptional amplification mediates the vulnerability of ATC cells to THZ1 treatment. Combining this integrative analysis with functional assays led to the discovery of a number of novel cancer genes of ATC, including PPP1R15A, SMG9, and KLF2. Inhibition of PPP1R15A with Guanabenz or Sephin1 greatly suppresses ATC growth. Significantly, the expression level of PPP1R15A is correlated with CDK7 expression in ATC tissue samples. Elevated expression of PPP1R15A and CDK7 are both associated with poor clinical prognosis in ATC patients. Importantly, CDK7 or PPP1R15A inhibition sensitizes ATC cells to conventional chemotherapy. Conclusions: Taken together, these findings demonstrate transcriptional addiction in ATC pathobiology and identify CDK7 and PPP1R15A as potential biomarkers and therapeutic targets for ATC.

Time Series Integrative Analysis of RNA Sequencing and MicroRNA Expression Data Reveals Key Biologic Wound Healing Pathways in Keloid-Prone Individuals

Time Series Integrative Analysis of RNA Sequencing and MicroRNA Expression Data Reveals Key Biologic Wound Healing Pathways in Keloid-Prone Individuals

RNA Binding Protein Ybx2 Regulates RNA Stability During Cold-Induced Brown Fat Activation.

Recent years have seen an upsurge of interest in brown adipose tissue (BAT) to combat the epidemic of obesity and diabetes. How its development and activation are regulated at the posttranscriptional level, however, has yet to be fully understood. RNA binding proteins (RBPs) lie in the center of posttranscriptional regulation. To systemically study the role of RBPs in BAT, we profiled >400 RBPs in different adipose depots and identified Y-box binding protein 2 (Ybx2) as a novel regulator in BAT activation. Knockdown of Ybx2 blocks brown adipogenesis, whereas its overexpression promotes BAT marker expression in brown and white adipocytes. Ybx2-knockout mice could form BAT but failed to express a full thermogenic program. Integrative analysis of RNA sequencing and RNA-immunoprecipitation study revealed a set of Ybx2's mRNA targets, including Pgc1α, that were destabilized by Ybx2 depletion during cold-induced activation. Thus, Ybx2 is a novel regulator that controls BAT activation by regulating mRNA stability.

Abstract 3459: SRC family kinase FYN promotes MET tyrosine kinase activation, epithelial to mesenchymal transition and metastasis in human prostate cancer

Abstract Introduction: FYN is a Src family kinase (SFK) member that has been shown to be upregulated in human prostate cancer tissues and cell lines (PCa). The majority of the studies on SRC kinases have focused on the classical c-SRC. We demonstrated previously that knockdown of FYN in PCa cell lines leads to reduced tumor growth in vitro and in vivo. We hypothesized that MET, an oncogenic tyrosine kinase receptor implicated in metastasis and therapeutic resistance of a variety of solid tumors is regulated by FYN expression in prostate cancer. Methods: We performed in vitro growth assays, migration and invasion assays and studied PCa cell response to HGF in wildtype, FYN-depleted and FYN-reconstituted cells. Additionally, we studied the role of FYN in promoting metastasis in vivo in xenograft mouse models. Summary: In this study, we demonstrate that knockdown of FYN leads to reduced migration, invasion and metastatic dissemination of PCa cells in vitro and in vivo. Conversely, ectopic expression of FYN promotes enhanced tumor growth and invasion. Hepatocyte growth factor (HGF) induced activation of MET is perturbed in FYN knockdown cells and overexpression of FYN leads to excessive MET activation suggesting a possible feedback loop between FYN and MET in prostate cancer. Consistent these findings, PC3 prostate cancer cells overexpressing FYN are resistant to treatment with a MET inhibitor upon ligand stimulation. More importantly, we show that FYN promotes epithelial to mesenchymal transition (EMT) of PCa cells and knockdown of FYN leads to reversal of EMT and alters radiation sensitivity of PCa cells. To further identify the molecular mechanism of FYN mediated regulation of tumor growth, we performed integrated RNA sequencing analysis of FYN knockdown cells and found MET and a network of molecular targets that appear to regulate FYN mediated tumor progression and metastasis. Conclusions: We demonstrate that FYN mediated MET activation promotes PCa invasion, EMT and metastasis. None of the SFK inhibitors currently in clinical trials are specific for SFKs alone. Therefore, it is imperative to develop and test novel inhibitors that selectively target FYN and other SFKs. Our findings have therapeutic implications since clinical trials are being performed for both SFK inhibitors and MET inhibitors in castration resistant PCa patients. Our studies reveal that FYN and MET interact in PCa cells and this interaction is clinically relevant given the success of XL-184 (a MET inhibitor) in phase III clinical trials in human prostate cancer. Clinical studies with translational endpoints are the need of the hour with agents that effectively inhibit FYN and MET in a variety of settings. Citation Format: Murali Gururajan, Margarit Sievert, Sheldon Mink, Chia-Yi (Gina) Chu, Matteo Morello, Jake Lichterman, Michael R. Freeman, Edwin M. Posadas. SRC family kinase FYN promotes MET tyrosine kinase activation, epithelial to mesenchymal transition and metastasis in human prostate cancer. [abstract]. In: Proceedings of the 105th Annual Meeting of the American Association for Cancer Research; 2014 Apr 5-9; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2014;74(19 Suppl):Abstract nr 3459. doi:10.1158/1538-7445.AM2014-3459