scholarly journals Gene Expression Profile in Similar Tissues Using Transcriptome Sequencing Data of Whole-Body Horse Skeletal Muscle

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 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.


2021 ◽  
Vol 22 (1) ◽  
Author(s):  
Chunyan Li ◽  
Xiaoyun He ◽  
Zijun Zhang ◽  
Chunhuan Ren ◽  
Mingxing Chu

Abstract Background Long noncoding RNA (lncRNA) has been identified as important regulator in hypothalamic-pituitary-ovarian axis associated with sheep prolificacy. However, little is known of their expression pattern and potential roles in the pineal gland of sheep. Herein, RNA-Seq was used to detect transcriptome expression pattern in pineal gland between follicular phase (FP) and luteal phase (LP) in FecBBB (MM) and FecB++ (ww) STH sheep, respectively, and differentially expressed (DE) lncRNAs and mRNAs associated with reproduction were identified. Results Overall, 135 DE lncRNAs and 1360 DE mRNAs in pineal gland between MM and ww sheep were screened. Wherein, 39 DE lncRNAs and 764 DE mRNAs were identified (FP vs LP) in MM sheep, 96 DE lncRNAs and 596 DE mRNAs were identified (FP vs LP) in ww sheep. Moreover, GO and KEGG enrichment analysis indicated that the targets of DE lncRNAs and DE mRNAs were annotated to multiple biological processes such as phototransduction, circadian rhythm, melanogenesis, GSH metabolism and steroid biosynthesis, which directly or indirectly participate in hormone activities to affect sheep reproductive performance. Additionally, co-expression of lncRNAs-mRNAs and the network construction were performed based on correlation analysis, DE lncRNAs can modulate target genes involved in related pathways to affect sheep fecundity. Specifically, XLOC_466330, XLOC_532771, XLOC_028449 targeting RRM2B and GSTK1, XLOC_391199 targeting STMN1, XLOC_503926 targeting RAG2, XLOC_187711 targeting DLG4 were included. Conclusion All of these differential lncRNAs and mRNAs expression profiles in pineal gland provide a novel resource for elucidating regulatory mechanism underlying STH sheep prolificacy.


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.


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.


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.


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.


F1000Research ◽  
2017 ◽  
Vol 6 ◽  
pp. 2010 ◽  
Author(s):  
Monther Alhamdoosh ◽  
Charity W. Law ◽  
Luyi Tian ◽  
Julie M. Sheridan ◽  
Milica Ng ◽  
...  

Gene set enrichment analysis is a popular approach for prioritising the biological processes perturbed in genomic datasets. The Bioconductor project hosts over 80 software packages capable of gene set analysis. Most of these packages search for enriched signatures amongst differentially regulated genes to reveal higher level biological themes that may be missed when focusing only on evidence from individual genes. With so many different methods on offer, choosing the best algorithm and visualization approach can be challenging. The EGSEA package solves this problem by combining results from up to 12 prominent gene set testing algorithms to obtain a consensus ranking of biologically relevant results.This workflow demonstrates how EGSEA can extend limma-based differential expression analyses for RNA-seq and microarray data using experiments that profile 3 distinct cell populations important for studying the origins of breast cancer. Following data normalization and set-up of an appropriate linear model for differential expression analysis, EGSEA builds gene signature specific indexes that link a wide range of mouse or human gene set collections obtained from MSigDB, GeneSetDB and KEGG to the gene expression data being investigated. EGSEA is then configured and the ensemble enrichment analysis run, returning an object that can be queried using several S4 methods for ranking gene sets and visualizing results via heatmaps, KEGG pathway views, GO graphs, scatter plots and bar plots. Finally, an HTML report that combines these displays can fast-track the sharing of results with collaborators, and thus expedite downstream biological validation. EGSEA is simple to use and can be easily integrated with existing gene expression analysis pipelines for both human and mouse data.


2019 ◽  
Author(s):  
Alemu Takele Assefa ◽  
Jo Vandesompele ◽  
Olivier Thas

SummarySPsimSeq is a semi-parametric simulation method for bulk and single cell RNA sequencing data. It simulates data from a good estimate of the actual distribution of a given real RNA-seq dataset. In contrast to existing approaches that assume a particular data distribution, our method constructs an empirical distribution of gene expression data from a given source RNA-seq experiment to faithfully capture the data characteristics of real data. Importantly, our method can be used to simulate a wide range of scenarios, such as single or multiple biological groups, systematic variations (e.g. confounding batch effects), and different sample sizes. It can also be used to simulate different gene expression units resulting from different library preparation protocols, such as read counts or UMI counts.Availability and implementationThe R package and associated documentation is available from https://github.com/CenterForStatistics-UGent/SPsimSeq.Supplementary informationSupplementary data are available at bioRχiv online.


Sign in / Sign up

Export Citation Format

Share Document