scholarly journals Identification of hypoxia-specific biomarkers in salmonids using RNA-sequencing and validation using high-throughput qPCR

2020 ◽  
Author(s):  
Arash Akbarzadeh ◽  
Aimee Lee S. Houde ◽  
Ben J.G. Sutherland ◽  
Oliver P. Günther ◽  
Kristina M. Miller

AbstractIdentifying early gene expression responses to hypoxia (i.e., low dissolved oxygen) as a tool to assess the degree of exposure to this stressor is crucial for salmonids, because they are increasingly exposed to hypoxic stress due to anthropogenic habitat change, e.g., global warming, excessive nutrient loading, and persistent algal blooms. Our goal was to discover and validate gill gene expression biomarkers specific to the hypoxia response in salmonids across multi-stressor conditions. Gill tissue was collected from 24 freshwater juvenile Chinook salmon (Oncorhynchus tshawytscha), held in normoxia [dissolved oxygen (DO) > 8 mg L−1] and hypoxia (DO = 4□5 mg L−1) in 10 and 18°C temperatures for up to six days. RNA-sequencing (RNA-seq) was then used to discover 240 differentially expressed genes between hypoxic and normoxic conditions, but not affected by temperature. The most significantly differentially expressed genes had functional roles in the cell cycle and suppression of cell proliferation associated with hypoxic conditions. The most significant genes (n = 30) were selected for real-time qPCR assay development. These assays demonstrated a strong correlation (r = 0.88; p < 0.001) between the expression values from RNA-seq and the fold changes from qPCR. Further, qPCR of the 30 candidate hypoxia biomarkers was applied to an additional 322 Chinook salmon exposed to hypoxic and normoxic conditions to reveal the top biomarkers to define hypoxic stress. Multivariate analyses revealed that smolt stage, water salinity, and morbidity status were relevant factors to consider with the expression of these genes in relation to hypoxic stress. These hypoxia candidate genes will be put into application screening Chinook salmon to determine the identity of stressors impacting the fish.

2020 ◽  
Vol 10 (9) ◽  
pp. 3321-3336
Author(s):  
Arash Akbarzadeh ◽  
Aimee Lee S Houde ◽  
Ben J G Sutherland ◽  
Oliver P Günther ◽  
Kristina M Miller

Abstract Identifying early gene expression responses to hypoxia (i.e., low dissolved oxygen) as a tool to assess the degree of exposure to this stressor is crucial for salmonids, because they are increasingly exposed to hypoxic stress due to anthropogenic habitat change, e.g., global warming, excessive nutrient loading, and persistent algal blooms. Our goal was to discover and validate gill gene expression biomarkers specific to the hypoxia response in salmonids across multi-stressor conditions. Gill tissue was collected from 24 freshwater juvenile Chinook salmon (Oncorhynchus tshawytscha), held in normoxia [dissolved oxygen (DO) &gt; 8 mg L-1] and hypoxia (DO = 4‒5 mg L-1) in 10 and 18° temperatures for up to six days. RNA-sequencing (RNA-seq) was then used to discover 240 differentially expressed genes between hypoxic and normoxic conditions, but not affected by temperature. The most significantly differentially expressed genes had functional roles in the cell cycle and suppression of cell proliferation associated with hypoxic conditions. The most significant genes (n = 30) were selected for real-time qPCR assay development. These assays demonstrated a strong correlation (r = 0.88; P &lt; 0.001) between the expression values from RNA-seq and the fold changes from qPCR. Further, qPCR of the 30 candidate hypoxia biomarkers was applied to an additional 322 Chinook salmon exposed to hypoxic and normoxic conditions to reveal the top biomarkers to define hypoxic stress. Multivariate analyses revealed that smolt stage, water salinity, and morbidity status were relevant factors to consider with the expression of these genes in relation to hypoxic stress. These hypoxia candidate genes will be put into application screening Chinook salmon to determine the identity of stressors impacting the fish.


2018 ◽  
Author(s):  
Adam McDermaid ◽  
Brandon Monier ◽  
Jing Zhao ◽  
Qin Ma

AbstractDifferential gene expression (DGE) is one of the most common applications of RNA-sequencing (RNA-seq) data. This process allows for the elucidation of differentially expressed genes (DEGs) across two or more conditions. Interpretation of the DGE results can be non-intuitive and time consuming due to the variety of formats based on the tool of choice and the numerous pieces of information provided in these results files. Here we present an R package, ViDGER (Visualization of Differential Gene Expression Results using R), which contains nine functions that generate information-rich visualizations for the interpretation of DGE results from three widely-used tools, Cuffdiff, DESeq2, and edgeR.


2019 ◽  
Vol 9 (2) ◽  
pp. 18 ◽  
Author(s):  
Isaac D. Raplee ◽  
Alexei V. Evsikov ◽  
Caralina Marín de Evsikova

The rapid expansion of transcriptomics and affordability of next-generation sequencing (NGS) technologies generate rocketing amounts of gene expression data across biology and medicine, including cancer research. Concomitantly, many bioinformatics tools were developed to streamline gene expression and quantification. We tested the concordance of NGS RNA sequencing (RNA-seq) analysis outcomes between two predominant programs for read alignment, HISAT2, and STAR, and two most popular programs for quantifying gene expression in NGS experiments, edgeR and DESeq2, using RNA-seq data from breast cancer progression series, which include histologically confirmed normal, early neoplasia, ductal carcinoma in situ and infiltrating ductal carcinoma samples microdissected from formalin fixed, paraffin embedded (FFPE) breast tissue blocks. We identified significant differences in aligners’ performance: HISAT2 was prone to misalign reads to retrogene genomic loci, STAR generated more precise alignments, especially for early neoplasia samples. edgeR and DESeq2 produced similar lists of differentially expressed genes, with edgeR producing more conservative, though shorter, lists of genes. Gene Ontology (GO) enrichment analysis revealed no skewness in significant GO terms identified among differentially expressed genes by edgeR versus DESeq2. As transcriptomics of FFPE samples becomes a vanguard of precision medicine, choice of bioinformatics tools becomes critical for clinical research. Our results indicate that STAR and edgeR are well-suited tools for differential gene expression analysis from FFPE samples.


Viruses ◽  
2021 ◽  
Vol 13 (2) ◽  
pp. 244 ◽  
Author(s):  
Antonio Victor Campos Coelho ◽  
Rossella Gratton ◽  
João Paulo Britto de Melo ◽  
José Leandro Andrade-Santos ◽  
Rafael Lima Guimarães ◽  
...  

HIV-1 infection elicits a complex dynamic of the expression various host genes. High throughput sequencing added an expressive amount of information regarding HIV-1 infections and pathogenesis. RNA sequencing (RNA-Seq) is currently the tool of choice to investigate gene expression in a several range of experimental setting. This study aims at performing a meta-analysis of RNA-Seq expression profiles in samples of HIV-1 infected CD4+ T cells compared to uninfected cells to assess consistently differentially expressed genes in the context of HIV-1 infection. We selected two studies (22 samples: 15 experimentally infected and 7 mock-infected). We found 208 differentially expressed genes in infected cells when compared to uninfected/mock-infected cells. This result had moderate overlap when compared to previous studies of HIV-1 infection transcriptomics, but we identified 64 genes already known to interact with HIV-1 according to the HIV-1 Human Interaction Database. A gene ontology (GO) analysis revealed enrichment of several pathways involved in immune response, cell adhesion, cell migration, inflammation, apoptosis, Wnt, Notch and ERK/MAPK signaling.


Blood ◽  
2012 ◽  
Vol 120 (21) ◽  
pp. 4582-4582
Author(s):  
Wei Liao ◽  
Gwen Jordaan ◽  
Artur Jaroszewicz ◽  
Matteo Pellegrini ◽  
Sanjai Sharma

Abstract Abstract 4582 High throughput sequencing of cellular mRNA provides a comprehensive analysis of the transcriptome. Besides identifying differentially expressed genes in different cell types, it also provides information of mRNA isoforms and splicing alterations. We have analyzed two CLL specimens and a normal peripheral blood B cells mRNA by this approach and performed data analysis to identify differentially expressed and spliced genes. The result showed CLLs specimens express approximately 40% more transcripts compared to normal B cells. The FPKM data (fragment per kilobase of exon per million) revealed a higher transcript expression on chromosome 12 in CLL#1 indicating the presence of trisomy 12, which was confirmed by fluorescent in-situ hybridization assay. With a two-fold change in FPKM as a cutoff and a p value cutoff of 0.05 as compared to the normal B cell control, 415 genes and 174 genes in CLL#1 and 676 and 235 genes in CLL#2 were up and downregulated or differentially expressed. In these two CLL specimens, 45% to 75% of differentially expressed genes are common to both the CLL specimens indicating that genetically disparate CLL specimens have a high percentage of a core set of genes that are potentially important for CLL biology. Selected differentially expressed genes with increased expression (selectin P ligand, SELPLG, and adhesion molecule interacts with CXADR antigen 1, AMICA) and decreased (Fos, Jun, CD69 and Rhob) expression based on the FPKM from RNA-sequencing data were also analyzed in additional CLL specimens by real time PCR analysis. The expression data from RNA-seq closely matches the fold-change in expression as measured by RT-PCR analysis and confirms the validity of the RNA-seq analysis. Interestingly, Fos was identified as one of the most downregulated gene in CLL. Using the Cufflinks and Cuffdiff software, the splicing patterns of genes in CLL specimens and normal B cells were analyzed. Approximately, 1100 to 1250 genes in the two CLL specimens were significantly differentially spliced as compared to normal B cells. In this analysis as well, there is a core set of 800 common genes which are differentially spliced in the two CLL specimens. The RNA-sequencing analysis accurately identifies differentially expressed novel genes and splicing variations that will help us understand the biology of CLL. Disclosures: No relevant conflicts of interest to declare.


2019 ◽  
Vol 9 (1) ◽  
Author(s):  
Kunzhe Dong ◽  
Shuang Chang ◽  
Qingmei Xie ◽  
Peng Zhao ◽  
Huanmin Zhang

Abstract Very virulent plus Marek’s disease (MD) virus (vv + MDV) induces tumors in relatively resistant lines of chickens and early mortality in highly susceptible lines of chickens. The vv + MDV also triggers a series of cellular responses in both types of chickens. We challenged birds sampled from a highly inbred chicken line (line 63) that is relatively resistant to MD and from another inbred line (line 72) that is highly susceptible to MD with a vv + MDV. RNA-sequencing analysis was performed with samples extracted from spleen tissues taken at 10-day and 21-day post infection (dpi). A total of 64 and 106 differentially expressed genes was identified in response to the vv + MDV challenge at latent phase in the resistant and susceptible lines of chickens, respectively. Direct comparisons between samples of the two lines identified 90 and 126 differentially expressed genes for control and MDV challenged groups, respectively. The differentially expressed gene profiles illustrated that intensive defense responses were significantly induced by vv + MDV at 10 dpi and 21 dpi but with slight changes in the resistant line. In contrast, vv + MDV induced a measurable suppression of gene expression associated with host defense at 10 dpi but followed by an apparent activation of the defense response at 21 dpi in the susceptible line of chickens. The observed difference in gene expression between the two genetic lines of chickens in response to MDV challenge during the latent phase provided a piece of indirect evidence that time points for MDV reactivation differ between the genetic lines of chickens with different levels of genetic resistance to MD. Early MDV reactivation might be necessary and potent to host defense system readiness for damage control of tumorigenesis and disease progression, which consequently results in measurable differences in phenotypic characteristics including early mortality (8 to 20 dpi) and tumor incidence between the resistant and susceptible lines of chickens. Combining differential gene expression patterns with reported GO function terms and quantitative trait loci, a total of 27 top genes was selected as highly promising candidate genes for genetic resistance to MD. These genes are functionally involved with virus process (F13A1 and HSP90AB1), immunity (ABCB1LB, RGS5, C10ORF58, OSF-2, MMP7, CXCL12, GAL1, GAL2, GAL7, HVCN1, PDE4D, IL4I1, PARP9, EOMES, MPEG1, PDK4, CCLI10, K60 and FST), and tumor suppression (ADAMTS2, LXN, ARRDC3, WNT7A, CLDN1 and HPGD). It is anticipated that these findings will facilitate advancement in the fundamental understanding on mechanisms of genetic resistance to MD. In addition, such advancement may also provide insights on tumor virus-induced tumorigenesis in general and help the research community recognize MD study may serve as a good model for oncology study involving tumor viruses.


2021 ◽  
Vol 39 (6_suppl) ◽  
pp. 478-478
Author(s):  
Zhichao Fu ◽  
Shenghua Liu ◽  
Jianfei Wang ◽  
Ning He ◽  
Yadong Yang ◽  
...  

478 Background: Bladder cancer is the ninth most common malignancy in the world, approximately 75% of patients are diagnosed with non-muscle invasive bladder cancer (NMIBC). Smoking has been established to be a carcinogenic risk factor of bladder cancer. Nevertheless, the detailed relationship between smoking and progression of NMIBC are poorly understood. In this study, we revealed high expressed genes in smoking patients were significantly related to tumor progression in NMIBC patients. Methods: A total of 54 NMIBC patients including 19 never smokers and 35 smokers (current smokers and previous smokers) were enrolled in this study.The gene expression profiles were obtained by RNA-seq and the differentially expressed genes between smoking and non-smoking patients were identified using DESeq2 .The further analysis of the association between genes expression and patient survival in NMIBC cohorts(Jakob et al., 2016)and IMvigor 210 cohorts(Jonathan et al., 2016)by Kaplan-Meier survival estimate. Results: We identified 46 differentially expressed genes (p<0.05) in smoking and non-somking NMIBC patients. IDO1 and KRT14 gene, which related to bladder cancer progression and poor prognosis, was identified significantly higher expressed in somking group compared with non-smoking and they have a logFC of 2.6,3.9 with FDR 1.83E-5,3.40E-5 respectively. The expression of other genes, including KRT6A, CASP14, SERPINA1, MYO3A and IL20RB, were significantly higher in smoking patients compared to non-somking. Notably, survival data analysis from 476 NMIBC cohorts showed that IL20RB had a significant relationship with poor PFS(p = 0.021) and in the Mvigor 210 Cohort including 310 advanced or metastatic urothelial carcinoma patients treated with atezolizumab, we found that the high expression of IL20RB was significantly related to poor OS(p = 0.002). Conclusions: We identified 14 genes related to tumor progression were significantly higher in smoking NMIBC patients than in non-smoking. Among these genes, the expression of IL20RB was related to the poor prognosis of NMIBC, and it may correlates with reduced clinical benefit of immunotherapeutic in patients with urothelial carcinoma.


PeerJ ◽  
2016 ◽  
Vol 4 ◽  
pp. e1991 ◽  
Author(s):  
Yanping Li ◽  
Shilin Tian ◽  
Xiaojun Yang ◽  
Xin Wang ◽  
Yuhai Guo ◽  
...  

Physcion and chrysophanol induce defense responses against powdery mildew in cucumbers. The combination of these two compounds has synergistic interaction against the disease. We performed RNA-seq on cucumber leaf samples treated with physcion and chrysophanol alone and with their combination. We generated 17.6 Gb of high-quality sequencing data (∼2 Gb per sample) and catalogued the expressions profiles of 12,293 annotated cucumber genes in each sample. We identified numerous differentially expressed genes that exhibited distinct expression patterns among the three treatments. The gene expression patterns of the Chr and Phy treatments were more similar to each other than to the Phy × Chr treatment. The Phy × Chr treatment induced the highest number of differentially expressed genes. This dramatic transcriptional change after Phy × Chr treatment leaves reflects that physcion combined with chrysophanol treatment was most closely associated with induction of disease resistance. The analysis showed that the combination treatment caused expression changes of numerous defense-related genes. These genes have known or potential roles in structural, chemical and signaling defense responses and were enriched in functional gene categories potentially responsible for cucumber resistance. These results clearly demonstrated that disease resistance in cucumber leaves was significantly influenced by the combined physcion and chrysophanol treatment. Thus, physcion and chrysophanol are appealing candidates for further investigation of the gene expression and associated regulatory mechanisms related to the defense response.


2020 ◽  
Vol 11 ◽  
Author(s):  
Xiaogang Cui ◽  
Shengli Zhang ◽  
Qin Zhang ◽  
Xiangyu Guo ◽  
Changxin Wu ◽  
...  

A total of 31 differentially expressed genes in the mammary glands were identified in our previous study using RNA sequencing (RNA-Seq), for lactating cows with extremely high and low milk protein and fat percentages. To determine the regulation of milk composition traits, we herein investigated the expression profiles of microRNA (miRNA) using small RNA sequencing based on the same samples as in the previous RNA-Seq experiment. A total of 497 known miRNAs (miRBase, release 22.1) and 49 novel miRNAs among the reads were identified. Among these miRNAs, 71 were found differentially expressed between the high and low groups (p &lt; 0.05, q &lt; 0.05). Furthermore, 21 of the differentially expressed genes reported in our previous RNA-Seq study were predicted as target genes for some of the 71 miRNAs. Gene ontology and KEGG pathway analyses showed that these targets were enriched for functions such as metabolism of protein and fat, and development of mammary gland, which indicating the critical role of these miRNAs in regulating the formation of milk protein and fat. With dual luciferase report assay, we further validated the regulatory role of 7 differentially expressed miRNAs through interaction with the specific sequences in 3′UTR of the targets. In conclusion, the current study investigated the complexity of the mammary gland transcriptome in dairy cattle using small RNA-seq. Comprehensive analysis of differential miRNAs expression and the data from previous study RNA-seq provided the opportunity to identify the key candidate genes for milk composition traits.


2018 ◽  
Vol 12 (1) ◽  
pp. 41-52 ◽  
Author(s):  
Bradford W. Lee ◽  
Virender B. Kumar ◽  
Pooja Biswas ◽  
Audrey C. Ko ◽  
Ramzi M. Alameddine ◽  
...  

Objective: This study utilized Next Generation Sequencing (NGS) to identify differentially expressed transcripts in orbital adipose tissue from patients with active Thyroid Eye Disease (TED) versus healthy controls. Method: This prospective, case-control study enrolled three patients with severe, active thyroid eye disease undergoing orbital decompression, and three healthy controls undergoing routine eyelid surgery with removal of orbital fat. RNA Sequencing (RNA-Seq) was performed on freshly obtained orbital adipose tissue from study patients to analyze the transcriptome. Bioinformatics analysis was performed to determine pathways and processes enriched for the differential expression profile. Quantitative Reverse Transcriptase-Polymerase Chain Reaction (qRT-PCR) was performed to validate the differential expression of selected genes identified by RNA-Seq. Results: RNA-Seq identified 328 differentially expressed genes associated with active thyroid eye disease, many of which were responsible for mediating inflammation, cytokine signaling, adipogenesis, IGF-1 signaling, and glycosaminoglycan binding. The IL-5 and chemokine signaling pathways were highly enriched, and very-low-density-lipoprotein receptor activity and statin medications were implicated as having a potential role in TED. Conclusion: This study is the first to use RNA-Seq technology to elucidate differential gene expression associated with active, severe TED. This study suggests a transcriptional basis for the role of statins in modulating differentially expressed genes that mediate the pathogenesis of thyroid eye disease. Furthermore, the identification of genes with altered levels of expression in active, severe TED may inform the molecular pathways central to this clinical phenotype and guide the development of novel therapeutic agents.


Sign in / Sign up

Export Citation Format

Share Document