scholarly journals GenomicSuperSignature: interpretation of RNA-seq experiments through robust, efficient comparison to public databases

2021 ◽  
Author(s):  
Sehyun Oh ◽  
Ludwig Geistlinger ◽  
Marcel Ramos ◽  
Jaclyn N. Taroni ◽  
Vincent J. Carey ◽  
...  

Millions of transcriptomic profiles have been deposited in public archives, yet remain underused for the interpretation of new experiments. We present a novel method for interpreting new transcriptomic datasets through near-instantaneous comparison to public archives without high-performance computing requirements. We performed Principal Component Analysis on 536 studies comprising 44,890 RNA sequencing profiles. Sufficiently similar loading vectors were aggregated to form Replicable Axes of Variation (RAV). RAVs were annotated with metadata of originating studies and by gene set enrichment analysis, forming a knowledge graph. Functionality to associate new datasets with RAVs, extract interpretable annotations, and provide intuitive visualization are implemented as the GenomicSuperSignature R/Bioconductor package. We demonstrated the efficient and coherent database searching, robustness to batch effects and heterogeneous training data, and transfer learning capacity of our method using TCGA and rare diseases datasets. GenomicSuperSignature will aid analyzing new gene expression data in the context of existing databases using minimal computing resources.

2021 ◽  
Author(s):  
Chengang Guo ◽  
Zhimin wei ◽  
Wei Lyu ◽  
Yanlou Geng

Abstract Quinoa saponins have complex, diverse and evident physiologic activities. However, the key regulatory genes for quinoa saponin metabolism are not yet well studied. The purpose of this study was to explore genes closely related to quinoa saponin metabolism. In this study, the significantly differentially expressed genes in yellow quinoa were firstly screened based on RNA-seq technology. Then, the key genes for saponin metabolism were selected by gene set enrichment analysis (GSEA) and principal component analysis (PCA) statistical methods. Finally, the specificity of the key genes was verified by hierarchical clustering. The results of differential analysis showed that 1654 differentially expressed genes were achieved after pseudogenes deletion. Therein, there were 142 long non-coding genes and 1512 protein-coding genes. Based on GSEA analysis, 116 key candidate genes were found to be significantly correlated with quinoa saponin metabolism. Through PCA dimension reduction analysis, 57 key genes were finally obtained. Hierarchical cluster analysis further demonstrated that these key genes can clearly separate the four groups of samples. The present results could provide references for the breeding of sweet quinoa and would be helpful for the rational utilization of quinoa saponins.


2021 ◽  
Vol 9 (3) ◽  
pp. e001610
Author(s):  
Incheol Seo ◽  
Hye Won Lee ◽  
Sang Jun Byun ◽  
Jee Young Park ◽  
Hyeonji Min ◽  
...  

BackgroundNeoadjuvant chemoradiation therapy (CRT) is a widely used preoperative treatment strategy for locally advanced rectal cancer (LARC). However, a few studies have evaluated the molecular changes caused by neoadjuvant CRT in these cancer tissues. Here, we aimed to investigate changes in immunotherapy-related immunogenic effects in response to preoperative CRT in LARC.MethodsWe analyzed 60 pairs of human LARC tissues before and after irradiation from three independent LARC cohorts, including a LARC patient RNA sequencing (RNA-seq) dataset from our cohort and GSE15781 and GSE94104 datasets.ResultsGene ontology analysis showed that preoperative CRT significantly enriched the immune response in LARC tissues. Moreover, gene set enrichment analysis revealed six significantly enriched Kyoto Encyclopedia of Genes and Genomes pathways associated with downregulated genes, including mismatch repair (MMR) genes, in LARC tissues after CRT in all three cohorts. Radiation also induced apoptosis and downregulated various MMR system-related genes in three colorectal cancer cells. One patient with LARC showed a change in microsatellite instability (MSI) status after CRT, as demonstrated by the loss of MMR protein and PCR for MSI. Moreover, CRT significantly increased tumor mutational burden in LARC tissues. CIBERSORT analysis revealed that the proportions of M2 macrophages and CD8 T cells were significantly increased after CRT in both the RNA-seq dataset and GSE94104. Notably, preoperative CRT increased various immune biomarker scores, such as the interferon-γ signature, the cytolytic activity and the immune signature.ConclusionsTaken together, our findings demonstrated that neoadjuvant CRT modulated the immune-related characteristics of LARC, suggesting that neoadjuvant CRT may enhance the responsiveness of LARC to immunotherapy.


2019 ◽  
Author(s):  
rui kong ◽  
Nan Wang ◽  
Wei Han ◽  
Yuejuan Zheng ◽  
Jie Lu

Abstract Background: In recent years, long non-coding RNAs (lncRNAs) are emerging as crucial regulators in the immunological process of liver hepatocellular carcinoma (LIHC). Increasing studies have found that some lncRNAs could be used as a diagnostic or therapeutic target for clinical management, but little research has investigated the role of immune-related lncRNA in tumor prognosis. In this study, we aimed to develop an immune lncRNA signature for the precise diagnosis and prognosis of liver hepatocellular carcinoma. Methods: Gene expression profiles of LIHC samples obtained from TCGA were screened for immune-related genes using two reference gene sets. The optimal immune-related lncRNA signature was built via correlational analysis, univariate and multivariate cox analysis. Then the Kaplan-Meier plot, ROC curve, clinical analysis, gene set enrichment analysis, and principal component analysis were carried out to evaluate the capability of immune lncRNA signature as a prognostic indicator. Results: Six long non-coding RNA MSC−AS1, AC009005.1, AL117336.3, AL031985.3, AL365203.2, AC099850.3 were identified via correlation analysis and cox regression analysis considering their interactions with immune genes. Next, tumor samples were separated into two risk groups by the signature with different clinical outcomes. Stratification analysis showed the prognostic ability of this signature acted as an independent factor. The AUC value of ROC curve was 0.779. The Kaplan-Meier method was used in survival analysis and results showed a statistical difference between the two risk groups. The predictive performance of this signature was validated by principal component analysis (PCA). Data from gene set enrichment analysis (GSEA) further unveiled several potential biological processes of these biomarkers may involve in. Conclusion: In summary, the study demonstrated the potential role of the six-lncRNA signature served as an independent prognostic factor for LIHC patients.


Genes ◽  
2020 ◽  
Vol 11 (4) ◽  
pp. 412
Author(s):  
Zhe Zhang ◽  
Zizengchen Wang ◽  
Yanna Dang ◽  
Jinyang Wang ◽  
Sakthidasan Jayaprakash ◽  
...  

The newly identified liver-enriched gene 1 (LEG1) encodes a protein with a characteristic domain of unknown function 781 (DUF781/LEG1), constituting a protein family with only one member in mammals. A functional study in zebrafish suggested that LEG1 genes are involved in liver development, while the platypus LEG1 homolog, Monotreme Lactation Protein (MLP), which is enriched in the mammary gland and milk, acts as an antibacterial substance. However, no functional studies on eutherian LEG1s have been published to date. Thus, we here report the first functional prediction study at the cellular level. As previously reported, eutherian LEG1s can be classified into three paralogous groups. Pigs have all three LEG1 genes (pLEG1s), while humans and mice have retained only LEG1a. Hence, pLEG1s might represent an ideal model for studying LEG1 gene functions. RNA-seq was performed by the overexpression of pLEG1s and platypus MLP in HepG2 cells. Enrichment analysis showed that pLEG1a and pLEG1b might exhibit little function in liver cells; however, pLEG1c is probably involved in the endoplasmic reticulum (ER) stress response and protein folding. Additionally, gene set enrichment analysis revealed that platypus MLP shows antibacterial activity, confirming the functional study in platypus. Therefore, our study showed from the transcriptomic perspective that mammalian LEG1s have different functions in liver cells due to the subfunctionalization of paralogous genes.


Neurosurgery ◽  
2019 ◽  
Vol 66 (Supplement_1) ◽  
Author(s):  
Robert Y North ◽  
Yan Li ◽  
Pradipta Ray ◽  
Laurence D Rhines ◽  
Claudio E Tatsui ◽  
...  

Abstract INTRODUCTION Women are at greater risk to suffer from many chronic pain conditions, more often report painful symptoms in epidemiological studies, and demonstrate greater pain sensitivity to experimentally measured pain responses. There is growing evidence from animal models for sex-specific biological differences in nociception, particularly involving primary afferent neurons, that may contribute to these differences. However, the details and extent of sex-specific differences associated with pain in human afferent neurons has not been previously investigated. METHODS Human dorsal root ganglia (DRG) and medical histories were obtained from patients undergoing spinal surgery that necessitated sacrifice of spinal nerve roots as part of standard of care. Clinical data for presence of painful radiculopathy was obtained through retrospective review of medical records or collected at study enrollment. RNA sequencing (RNA-seq) was performed on 21 DRG from 15 patients with variable presence of radicular pain reported in a corresponding dermatome. Differential expression analysis for male w/pain (MP) vs female w/pain (FP) samples was performed with thresholds for robustly expressed autosomal genes (TPM >3.0), fold change of 2.0 or higher, with false discovery rate (FDR) <0.05. RESULTS Comparison of the MP and FP cohorts yielded 575 differentially expressed genes with 426 upregulated in MP and 149 upregulated in FP. Gene set enrichment analysis demonstrated significant differences in genes related to inflammation and immune regulation (increased MAPK and BDNF signaling in MP, increased Rhodopsin-like GPCR in FP) and differing clusters of spinal cord injury-associated genes (TLR4, AIF1, OMG, C1QB increased in FP, EGR1, NR4A1, ZFP36, BTG2, MYC in MP). CONCLUSION Utilizing RNA-seq of human DRG innervating regions of pain, this study provides the first demonstration of sex-specific differences for the biology of pain within the dorsal root ganglion in humans and implicates the immune system as a critical influence in these differences.


Genes ◽  
2020 ◽  
Vol 11 (11) ◽  
pp. 1359
Author(s):  
Ho-Yeon Lee ◽  
Jae-Yoon Kim ◽  
Kyoung Hyoun Kim ◽  
Seongmun Jeong ◽  
Youngbum Cho ◽  
...  

Horses have been studied for exercise function rather than food production, unlike most livestock. Therefore, the role and characteristics of tissue landscapes are critically understudied, except for certain muscles used in exercise-related studies. In the present study, we compared RNA-Seq data from 18 Jeju horse skeletal muscles to identify differentially expressed genes (DEGs) between tissues that have similar functions and to characterize these differences. We identified DEGs between different muscles using pairwise differential expression (DE) analyses of tissue transcriptome expression data and classified the samples using the expression values of those genes. Each tissue was largely classified into two groups and their subgroups by k-means clustering, and the DEGs identified in comparison between each group were analyzed by functional/pathway level using gene set enrichment analysis and gene level, confirming the expression of significant genes. As a result of the analysis, the differences in metabolic properties like glycolysis, oxidative phosphorylation, and exercise adaptation of the groups were detected. The results demonstrated that the biochemical and anatomical features of a wide range of muscle tissues in horses could be determined through transcriptome expression analysis, and provided proof-of-concept data demonstrating that RNA-Seq analysis can be used to classify and study in-depth differences between tissues with similar properties.


2021 ◽  
Vol 11 (1) ◽  
Author(s):  
Ludmila Mudri Hul ◽  
Adriana Mércia Guaratini Ibelli ◽  
Igor Ricardo Savoldi ◽  
Débora Ester Petry Marcelino ◽  
Lana Teixeira Fernandes ◽  
...  

AbstractLocomotor problems are among one of the main concerns in the current poultry industry, causing major economic losses and affecting animal welfare. The most common bone anomalies in the femur are dyschondroplasia, femoral head separation (FHS), and bacterial chondronecrosis with osteomyelitis (BCO), also known as femoral head necrosis (FHN). The present study aimed to identify differentially expressed (DE) genes in the articular cartilage (AC) of normal and FHS-affected broilers by RNA-Seq analysis. In the transcriptome analysis, 12,169 genes were expressed in the femur AC. Of those, 107 genes were DE (FDR < 0.05) between normal and affected chickens, of which 9 were downregulated and 98 were upregulated in the affected broilers. In the gene-set enrichment analysis using the DE genes, 79 biological processes (BP) were identified and were grouped into 12 superclusters. The main BP found were involved in the response to biotic stimulus, gas transport, cellular activation, carbohydrate-derived catabolism, multi-organism regulation, immune system, muscle contraction, multi-organism process, cytolysis, leukocytes and cell adhesion. In this study, the first transcriptome analysis of the broilers femur articular cartilage was performed, and a set of candidate genes (AvBD1, AvBD2, ANK1, EPX, ADA, RHAG) that could trigger changes in the broiler´s femoral growth plate was identified. Moreover, these results could be helpful to better understand FHN in chickens and possibly in humans.


2019 ◽  
Vol 37 (7_suppl) ◽  
pp. 146-146
Author(s):  
William S. Chen ◽  
Rahul Raj Aggarwal ◽  
Li Zhang ◽  
Shuang Zhao ◽  
Tomasz M. Beer ◽  
...  

146 Background: Metastatic castration-resistant prostate cancer (mCRPC) is the lethal form of the disease. Several recent efforts have identified genomic alterations in mCRPC, but the clinical implications of these alterations have not been fully elucidated. We conducted a prospective cohort study (n = 101) using whole genome sequencing (WGS) to analyze the association between key driver gene alterations and overall survival. We also performed whole-transcriptome RNA sequencing (RNA-seq) analyses to identify potential mechanisms of enzalutamide resistance in mCRPC. Methods: Metastasis biopsies were obtained in 101 mCRPC patients as part of the multi-institutional West Coast Prostate Cancer Dream Team project. Samples underwent WGS and RNA-seq. The resulting mutation, copy number, and structural variant calls were integrated to determine functional copy number status of candidate genes for downstream clinical analyses. We performed univariate and multivariable analyses to assess the prognostic significance of candidate genomic events with respect to overall survival. To nominate and investigate genomic pathways associated with enzalutamide resistance, we performed expression-based gene set enrichment analysis followed by cross-sectional enrichment and survival analyses related to the top nominated pathway. Results: RB1 loss was associated with poor overall survival (median 14.1 vs. 42.0 months, p < 0.001). When we compared enzalutamide resistant versus naïve samples using gene set enrichment analysis, we identified the Wnt/beta-catenin pathway as the top differentially expressed pathway in enzalutamide-resistant patients. Furthermore, CTNNB1 (beta-catenin) activating mutations were exclusive to enzalutamide-resistant patients (p = 0.013) and predictive of poor overall survival (median 13.6 vs. 41.7 months, p < 0.001). Conclusions: Impaired survival in mCRPC patients is associated with RB1 loss, identified by integrated genomic analysis of CRPC metastasis biopsies. Among men with mCRPC that was enzalutamide-resistant, the Wnt/beta-catenin pathway is nominated as an important predictive (and potentially therapeutic) pathway.


2021 ◽  
Vol 11 ◽  
Author(s):  
Ming Gao ◽  
Xinzhuang Wang ◽  
Dayong Han ◽  
Enzhou Lu ◽  
Jian Zhang ◽  
...  

Glioblastoma multiforme (GBM) is the most aggressive primary tumor of the central nervous system. As biomedicine advances, the researcher has found the development of GBM is closely related to immunity. In this study, we evaluated the GBM tumor immunoreactivity and defined the Immune-High (IH) and Immune-Low (IL) immunophenotypes using transcriptome data from 144 tumors profiled by The Cancer Genome Atlas (TCGA) project based on the single-sample gene set enrichment analysis (ssGSEA) of five immune expression signatures (IFN-γ response, macrophages, lymphocyte infiltration, TGF-β response, and wound healing). Next, we identified six immunophenotype-related long non-coding RNA biomarkers (im-lncRNAs, USP30-AS1, HCP5, PSMB8-AS1, AL133264.2, LINC01684, and LINC01506) by employing a machine learning computational framework combining minimum redundancy maximum relevance algorithm (mRMR) and random forest model. Moreover, the expression level of identified im-lncRNAs was converted into an im-lncScore using the normalized principal component analysis. The im-lncScore showed a promising performance for distinguishing the GBM immunophenotypes with an area under the curve (AUC) of 0.928. Furthermore, the im-lncRNAs were also closely associated with the levels of tumor immune cell infiltration in GBM. In summary, the im-lncRNA signature had important clinical implications for tumor immunophenotyping and guiding immunotherapy in glioblastoma patients in future.


Blood ◽  
2018 ◽  
Vol 132 (Supplement 1) ◽  
pp. 41-41
Author(s):  
Krystle Nomie ◽  
Liang Zhang ◽  
Yixin Yao ◽  
Yang Liu ◽  
Shaojun Zhang ◽  
...  

Abstract Introduction Mantle cell lymphoma (MCL) is an incurable B-cell lymphoma subtype and constitutive activation of the B-cell receptor pathway is a hallmark of B-cell lymphomas. Bruton's tyrosine kinase (BTK) is a critical component of the B-cell receptor pathway, and ibrutinib, a first-in-class, once-daily, and oral covalent inhibitor of BTK, was developed to reduce/silence B-cell receptor pathway activity, leading to clinically remarkable anti-tumor activity. In our prior multiple-center Phase II clinical trial, the overall response rate in relapsed/refractory MCL patients was 68% (Wang et al., NEJM, 2013), surpassing the effectiveness of other therapies. Although ibrutinib is extremely efficacious in patients with relapsed/refractory MCL, the one-year overall survival rate of ibrutinib-exposed patients who relapse is only 22%. Methods Patient primary cells were isolated from MCL patients treated with ibrutinib either prior to treatment or at treatment discontinuation. Whole exome sequencing (WES) was performed to determine the mutational landscape of ibrutinib resistance. RNA-seq was employed to compare the gene expression profiles between ibrutinib-sensitive and -resistant patient samples. Gene set enrichment analysis was utilized to identify dysregulated molecular pathways associated with the resistant phenotype. The RNA-seq data were then validated with reverse phase protein array (RPPA) analysis of ibrutinib-sensitive and -resistant MCL cell lines. Metabolic assays including the measurement of mitochondria respiration rates with the Seahorse analyzer and reactive oxygen species (ROS) levels, targeted metabolomics, and ATP analysis. Functional studies targeting this molecular pathway were conducted, including in vitro cell viability and apoptosis assays, as well as in vivo efficacy studies in an ibrutinib-resistant MCL patient-derived xenograft mouse model. Results WES data analysis identified frequent inactivating somatic alterations in ATM, KMT2D, and TP53 in both the ibrutinib-sensitive and -resistant tumors. CDKN2A (5/7, 71%) was frequently deleted, and the deletion was only observed in the ibrutinib-resistant tumors (p = 0.010). The RNA-seq analysis identified a total of 63 protein-coding genes as the most differentially expressed genes (DEGs) between the ibrutinib-resistant and -sensitive groups, with a fold change of ≥ 2 or ≤ -2 and the false discovery rate (FDR q-value) ≤ 0.01. Among the DEGs, 26 genes were upregulated in ibrutinib-resistant tumors. In addition, gene set enrichment analysis (GSEA) revealed the marked upregulation of oncogenic pathways including c-MYC, mTOR (mTORC1), Wnt, and NF-ĸb signaling, followed by cell cycle, apoptosis, BCR signaling and DNA repair in the ibrutinib-resistant tumors. Notably, in addition to these oncogenic pathways, the metabolic pathways, including oxidative phosphorylation (OXPHOS), were significantly enriched in the ibrutinib-resistant tumors (normalized enrichment score > 3 and FDR q-value < 1e-5). In further support of this finding, metabolomics analysis and the measurement of ATP production and mitochondrial respiration indicated that the OXPHOS pathway is the predominant metabolic pathway employed by ibrutinib-resistant MCL cells. To determine the effects of targeting these pathways, OXPHOS was inhibited with a novel electron transport complex I inhibitor (IACS-010759, developed by MD Anderson Cancer Center) in both MCL cell lines and ibrutinib-resistant MCL patient-derived xenograft (PDX) models. Single agent IACS-010759 treatment at 10 mg/kg oral gavage for 5 consecutive days/week completely prevented tumor growth compared with the vehicle control as shown by measuring tumor volume (n = 5, p < 0.0001) and human β2M levels (n = 5; p < 0.0001) throughout treatment. No apparent toxicities were observed in the IACS-010759-treated MCL PDX mice. Conclusion This current study warrants the exploitation of active cancer metabolic pathways, especially OXPHOS, to improve the clinical outcomes of MCL and additional lymphoma, which is actively being investigated in a Phase I lymphoma clinical trial (NCT03291938). Disclosures Wang: AstraZeneca: Consultancy, Research Funding; Celgene: Honoraria, Membership on an entity's Board of Directors or advisory committees, Research Funding; Juno: Research Funding; Novartis: Research Funding; Janssen: Consultancy, Honoraria, Membership on an entity's Board of Directors or advisory committees, Research Funding; MoreHealth: Consultancy; Acerta Pharma: Honoraria, Research Funding; Kite Pharma: Research Funding; Pharmacyclics: Honoraria, Research Funding; Dava Oncology: Honoraria.


Sign in / Sign up

Export Citation Format

Share Document