scholarly journals Comparative Analysis of the Liver Transcriptome among Cattle Breeds Using RNA-seq

2019 ◽  
Vol 6 (2) ◽  
pp. 36
Author(s):  
Chandra Pareek ◽  
Mateusz Sachajko ◽  
Jedrzej Jaskowski ◽  
Magdalena Herudzinska ◽  
Mariusz Skowronski ◽  
...  

Global gene expression in liver transcriptome varies among cattle breeds. The present investigation was aimed to identify the differentially expressed genes (DEGs), metabolic gene networks and metabolic pathways in bovine liver transcriptome of young bulls. In this study, we comparatively analyzed the bovine liver transcriptome of dairy (Polish Holstein Friesian (HF); n = 6), beef (Hereford; n = 6), and dual purpose (Polish-Red; n = 6) cattle breeds. This study identified 895, 338, and 571 significant (p < 0.01) differentially expressed (DE) gene-transcripts represented as 745, 265, and 498 hepatic DE genes through the Polish-Red versus Hereford, Polish-HF versus Hereford, and Polish-HF versus Polish-Red breeds comparisons, respectively. By combining all breeds comparisons, 75 hepatic DE genes (p < 0.01) were identified as commonly shared among all the three breed comparisons; 70, 160, and 38 hepatic DE genes were commonly shared between the following comparisons: (i) Polish-Red versus Hereford and Polish-HF versus Hereford; (ii) Polish-Red versus Hereford and Polish-HF versus Polish-Red; and (iii) Polish-HF versus Hereford and Polish-HF versus Polish-Red, respectively. A total of 440, 82, and 225 hepatic DE genes were uniquely observed for the Polish-Red versus Hereford, Polish-HF versus Hereford, and Polish-Red versus Polish-HF comparisons, respectively. Gene ontology (GO) analysis identified top-ranked enriched GO terms (p < 0.01) including 17, 16, and 31 functional groups and 151, 61, and 140 gene functions that were DE in all three breed liver transcriptome comparisons. Gene network analysis identified several potential metabolic pathways involved in glutamine family amino-acid, triglyceride synthesis, gluconeogenesis, p38MAPK cascade regulation, cholesterol biosynthesis (Polish-Red versus Hereford); IGF-receptor signaling, catecholamine transport, lipoprotein lipase, tyrosine kinase binding receptor (Polish-HF versus Hereford), and PGF-receptor binding, (Polish-HF versus Polish-Red). Validation results showed that the relative expression values were consistent to those obtained by RNA-seq, and significantly correlated between the quantitative reverse transcription PCR (RT-qPCR) and RNA-seq (Pearson’s r  > 0.90). Our results provide new insights on bovine liver gene expressions among dairy versus dual versus beef breeds by identifying the large numbers of DEGs markers submitted to NCBI gene expression omnibus (GEO) accession number GSE114233, which can serve as useful genetic tools to develop the gene assays for trait-associated studies as well as, to effectively implement in genomics selection (GS) cattle breeding programs in Poland.

2021 ◽  
Vol 80 (Suppl 1) ◽  
pp. 12.2-12
Author(s):  
I. Muller ◽  
M. Verhoeven ◽  
H. Gosselt ◽  
M. Lin ◽  
T. De Jong ◽  
...  

Background:Tocilizumab (TCZ) is a monoclonal antibody that binds to the interleukin 6 receptor (IL-6R), inhibiting IL-6R signal transduction to downstream inflammatory mediators. TCZ has shown to be effective as monotherapy in early rheumatoid arthritis (RA) patients (1). However, approximately one third of patients inadequately respond to therapy and the biological mechanisms underlying lack of efficacy for TCZ remain elusive (1). Here we report gene expression differences, in both whole blood and peripheral blood mononuclear cells (PBMC) RNA samples between early RA patients, categorized by clinical TCZ response (reaching DAS28 < 3.2 at 6 months). These findings could lead to identification of predictive biomarkers for TCZ response and improve RA treatment strategies.Objectives:To identify potential baseline gene expression markers for TCZ response in early RA patients using an RNA-sequencing approach.Methods:Two cohorts of RA patients were included and blood was collected at baseline, before initiating TCZ treatment (8 mg/kg every 4 weeks, intravenously). DAS28-ESR scores were calculated at baseline and clinical response to TCZ was defined as DAS28 < 3.2 at 6 months of treatment. In the first cohort (n=21 patients, previously treated with DMARDs), RNA-sequencing (RNA-seq) was performed on baseline whole blood PAXgene RNA (Illumina TruSeq mRNA Stranded) and differential gene expression (DGE) profiles were measured between responders (n=14) and non-responders (n=7). For external replication, in a second cohort (n=95 therapy-naïve patients receiving TCZ monotherapy), RNA-seq was conducted on baseline PBMC RNA (SMARTer Stranded Total RNA-Seq Kit, Takara Bio) from the 2-year, multicenter, double-blind, placebo-controlled, randomized U-Act-Early trial (ClinicalTrials.gov identifier: NCT01034137) and DGE was analyzed between 84 responders and 11 non-responders.Results:Whole blood DGE analysis showed two significantly higher expressed genes in TCZ non-responders (False Discovery Rate, FDR < 0.05): urotensin 2 (UTS2) and caveolin-1 (CAV1). Subsequent analysis of U-Act-Early PBMC DGE showed nine differentially expressed genes (FDR < 0.05) of which expression in clinical TCZ non-responders was significantly higher for eight genes (MTCOP12, ZNF774, UTS2, SLC4A1, FECH, IFIT1B, AHSP, and SPTB) and significantly lower for one gene (TND2P28M). Both analyses were corrected for baseline DAS28-ESR, age and gender. Expression of UTS2, with a proposed function in regulatory T-cells (2), was significantly higher in TCZ non-responders in both cohorts. Furthermore, gene ontology enrichment analysis revealed no distinct gene ontology or IL-6 related pathway(s) that were significantly different between TCZ-responders and non-responders.Conclusion:Several genes are differentially expressed at baseline between responders and non-responders to TCZ therapy at 6 months. Most notably, UTS2 expression is significantly higher in TCZ non-responders in both whole blood as well as PBMC cohorts. UTS2 could be a promising target for further analyses as a potential predictive biomarker for TCZ response in RA patients in combination with clinical parameters (3).References:[1]Bijlsma JWJ, Welsing PMJ, Woodworth TG, et al. Early rheumatoid arthritis treated with tocilizumab, methotrexate, or their combination (U-Act-Early): a multicentre, randomised, double-blind, double-dummy, strategy trial. Lancet. 2016;388(10042):343-55.[2]Bhairavabhotla R, Kim YC, Glass DD, et al. Transcriptome profiling of human FoxP3+ regulatory T cells. Human Immunology. 2016;77(2):201-13.[3]Gosselt HR, Verhoeven MMA, Bulatovic-Calasan M, et al. Complex machine-learning algorithms and multivariable logistic regression on par in the prediction of insufficient clinical response to methotrexate in rheumatoid arthritis. Journal of Personalized Medicine. 2021;11(1).Disclosure of Interests:None declared


2021 ◽  
Vol 15 (1) ◽  
Author(s):  
Weitong Cui ◽  
Huaru Xue ◽  
Lei Wei ◽  
Jinghua Jin ◽  
Xuewen Tian ◽  
...  

Abstract Background RNA sequencing (RNA-Seq) has been widely applied in oncology for monitoring transcriptome changes. However, the emerging problem that high variation of gene expression levels caused by tumor heterogeneity may affect the reproducibility of differential expression (DE) results has rarely been studied. Here, we investigated the reproducibility of DE results for any given number of biological replicates between 3 and 24 and explored why a great many differentially expressed genes (DEGs) were not reproducible. Results Our findings demonstrate that poor reproducibility of DE results exists not only for small sample sizes, but also for relatively large sample sizes. Quite a few of the DEGs detected are specific to the samples in use, rather than genuinely differentially expressed under different conditions. Poor reproducibility of DE results is mainly caused by high variation of gene expression levels for the same gene in different samples. Even though biological variation may account for much of the high variation of gene expression levels, the effect of outlier count data also needs to be treated seriously, as outlier data severely interfere with DE analysis. Conclusions High heterogeneity exists not only in tumor tissue samples of each cancer type studied, but also in normal samples. High heterogeneity leads to poor reproducibility of DEGs, undermining generalization of differential expression results. Therefore, it is necessary to use large sample sizes (at least 10 if possible) in RNA-Seq experimental designs to reduce the impact of biological variability and DE results should be interpreted cautiously unless soundly validated.


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.


2012 ◽  
Vol 111 (suppl_1) ◽  
Author(s):  
Emma L Robinson ◽  
Syed Haider ◽  
Hillary Hei ◽  
Richard T Lee ◽  
Roger S Foo

Heart failure comprises of clinically distinct inciting causes but a consistent pattern of change in myocardial gene expression supports the hypothesis that unifying biochemical mechanisms underlie disease progression. The recent RNA-seq revolution has enabled whole transcriptome profiling, using deep-sequencing technologies. Up to 70% of the genome is now known to be transcribed into RNA, a significant proportion of which is long non-coding RNAs (lncRNAs), defined as polyribonucleotides of ≥200 nucleotides. This project aims to discover whether the myocardium expression of lncRNAs changes in the failing heart. Paired end RNA-seq from a 300-400bp library of ‘stretched’ mouse myocyte total RNA was carried out to generate 76-mer sequence reads. Mechanically stretching myocytes with equibiaxial stretch apparatus mimics pathological hypertrophy in the heart. Transcripts were assembled and aligned to reference genome mm9 (UCSC), abundance determined and differential expression of novel transcripts and alternative splice variants were compared with that of control (non-stretched) mouse myocytes. Five novel transcripts have been identified in our RNA-seq that are differentially expressed in stretched myocytes compared with non-stretched. These are regions of the genome that are currently unannotated and potentially are transcribed into non-coding RNAs. Roles of known lncRNAs include control of gene expression, either by direct interaction with complementary regions of the genome or association with chromatin remodelling complexes which act on the epigenome.Changes in expression of genes which contribute to the deterioration of the failing heart could be due to the actions of these novel lncRNAs, immediately suggesting a target for new pharmaceuticals. Changes in the expression of these novel transcripts will be validated in a larger sample size of stretched myocytes vs non-stretched myocytes as well as in the hearts of transverse aortic constriction (TAC) mice vs Sham (surgical procedure without the aortic banding). In vivo investigations will then be carried out, using siLNA antisense technology to silence novel lncRNAs in mice.


Author(s):  
Crescenzio Gallo

The possible applications of modeling and simulation in the field of bioinformatics are very extensive, ranging from understanding basic metabolic paths to exploring genetic variability. Experimental results carried out with DNA microarrays allow researchers to measure expression levels for thousands of genes simultaneously, across different conditions and over time. A key step in the analysis of gene expression data is the detection of groups of genes that manifest similar expression patterns. In this chapter, the authors examine various methods for analyzing gene expression data, addressing the important topics of (1) selecting the most differentially expressed genes, (2) grouping them by means of their relationships, and (3) classifying samples based on gene expressions.


2019 ◽  
Vol 9 (1) ◽  
Author(s):  
Giulia Pintarelli ◽  
Sara Noci ◽  
Davide Maspero ◽  
Angela Pettinicchio ◽  
Matteo Dugo ◽  
...  

Abstract Alterations in the gene expression of organs in contact with the environment may signal exposure to toxins. To identify genes in lung tissue whose expression levels are altered by cigarette smoking, we compared the transcriptomes of lung tissue between 118 ever smokers and 58 never smokers. In all cases, the tissue studied was non-involved lung tissue obtained at lobectomy from patients with lung adenocarcinoma. Of the 17,097 genes analyzed, 357 were differentially expressed between ever smokers and never smokers (FDR < 0.05), including 290 genes that were up-regulated and 67 down-regulated in ever smokers. For 85 genes, the absolute value of the fold change was ≥2. The gene with the smallest FDR was MYO1A (FDR = 6.9 × 10−4) while the gene with the largest difference between groups was FGG (fold change = 31.60). Overall, 100 of the genes identified in this study (38.6%) had previously been found to associate with smoking in at least one of four previously reported datasets of non-involved lung tissue. Seven genes (KMO, CD1A, SPINK5, TREM2, CYBB, DNASE2B, FGG) were differentially expressed between ever and never smokers in all five datasets, with concordant higher expression in ever smokers. Smoking-induced up-regulation of six of these genes was also observed in a transcription dataset from lung tissue of non-cancer patients. Among the three most significant gene networks, two are involved in immunity and inflammation and one in cell death. Overall, this study shows that the lung parenchyma transcriptome of smokers has altered gene expression and that these alterations are reproducible in different series of smokers across countries. Moreover, this study identified a seven-gene panel that reflects lung tissue exposure to cigarette smoke.


Blood ◽  
2011 ◽  
Vol 118 (21) ◽  
pp. 2367-2367
Author(s):  
Mira Jeong ◽  
Deqiang Sun ◽  
Min Luo ◽  
Aysegul Ergen ◽  
Hongcang Gu ◽  
...  

Abstract Abstract 2367 Hematopoietic stem cell (HSC) Aging is a complex process linked to number of changes in gene expression and functional decline of self-renewal and differentiation potential. While epigenetic changes have been implicated in HSC aging, little direct evidence has been generated. DNA methylation is one of the major underlying mechanisms associated with the regulation of gene expression, but changes in DNA methylation patterns with HSC aging have not been characterized. We hypothesize that revealing the genome-wide DNA methylation and transcriptome signatures will lead to a greater understanding of HSC aging. Here, we report the first genome-scale study of epigenomic dynamics during normal mouse HSC aging. We isolated SP-KSL-CD150+ HSC populations from 4, 12, 24 month-old mouse bone marrow and carried out genome-wide reduced representative bisulfite sequencing (RRBS) and identified aging-associated differentially methylated CpGs. Three biological samples were sequenced from each aging group and we obtained 30–40 million high-quality reads with over 30X total coverage on ∼1.1M CpG sites which gives us adequate statistical power to infer methylation ratios. Bisulfite conversion rate of non-CpG cytosines was >99%. We analyzed a variety of genomic features to find that CpG island promoters, gene bodies, 5'UTRs, and 3'UTRs generally were associated with hypermethylation in aging HSCs. Overall, out of 1,777 differentially methylated CpGs, 92.8% showed age-related hypermethylation and 7.2% showed age-related hypomethylation. Gene ontology analyses have revealed that differentially methylated CpGs were significantly enriched near genes associated with alternative splicing, DNA binding, RNA-binding, transcription regulation, Wnt signaling and pathways in cancer. Most interestingly, over 579 splice variants were detected as candidates for age-related hypermethylation (86%) and hypomethylation (14%) including Dnmt3a, Runx1, Pbx1 and Cdkn2a. To quantify differentially expressed RNA-transcripts across the entire transcriptome, we performed RNA-seq and analyzed exon arrays. The Spearman's correlation between two different methods was good (r=0.80). From exon arrays, we identified 586 genes that were down regulated and 363 gene were up regulated with aging (p<0.001). Most interestingly, overall expression of DNA methyl transferases Dnmt1, Dnmt3a, Dnmt3b were down regulated with aging. We also found that Dnmt3a2, the short isoform of Dnmt3a, which lacks the N-terminal region of Dnmt3a and represents the major isoform in ES cells, is more expressed in young HSC. For the RNA-seq analysis, we focused first on annotated transcripts derived from cloned mRNAs and we found 307 genes were down regulated and 1015 gene were up regulated with aging (p<0.05). Secondly, we sought to identify differentially expressed isoforms and also novel transcribed regions (antisense and novel genes). To characterize the genes showing differential regulation, we analyzed their functional associations and observed that the highest scoring annotation cluster was enriched in genes associated with translation, the immune network and hematopoietic cell lineage. We expect that the results of these experiments will reveal the global effect of DNA methylation on transcript stability and the translational state of target genes. Our findings will lend insight into the molecular mechanisms responsible for the pathologic changes associated with aging in HSCs. Disclosures: No relevant conflicts of interest to declare.


2008 ◽  
Vol 35 (2) ◽  
pp. 182-190 ◽  
Author(s):  
Toshimori Kitami ◽  
Renee Rubio ◽  
William O'Brien ◽  
John Quackenbush ◽  
Joseph H. Nadeau

Dietary folate supplementation can dramatically reduce the severity and incidence of several common birth defects and adult diseases that are associated with anomalies in homocysteine and folate metabolism. The common polymorphisms that adversely affect these metabolic pathways do not fully account for the particular birth defects and adult diseases that occur in at-risk individuals. To test involvement of folate, homocysteine, and other pathways in disease pathogenesis and treatment response, we analyzed global and pathway-specific changes in gene expression and levels of selected metabolites after depletion and repletion of dietary folate in two genetically distinct inbred strains of mice. Compared with the C57BL/6J strain, A/J showed greater homeostatic response to folate perturbation by retaining a higher serum folate level and minimizing global gene expression changes. Remarkably, folate perturbation led to systematic strain-specific differences only in the expression profile of the cholesterol biosynthesis pathway and to changes in levels of serum and liver total cholesterol. By genetically increasing serum and liver total cholesterol levels in APOE-deficient mice, we modestly but significantly improved folate retention during folate depletion, suggesting that homeostasis among the homocysteine, folate and cholesterol metabolic pathways contributes to the beneficial effects of dietary folate supplementation.


2014 ◽  
Vol 42 (15) ◽  
pp. e121-e121 ◽  
Author(s):  
Hari Krishna Yalamanchili ◽  
Zhaoyuan Li ◽  
Panwen Wang ◽  
Maria P. Wong ◽  
Jianfeng Yao ◽  
...  

Abstract Conventionally, overall gene expressions from microarrays are used to infer gene networks, but it is challenging to account splicing isoforms. High-throughput RNA Sequencing has made splice variant profiling practical. However, its true merit in quantifying splicing isoforms and isoform-specific exon expressions is not well explored in inferring gene networks. This study demonstrates SpliceNet, a method to infer isoform-specific co-expression networks from exon-level RNA-Seq data, using large dimensional trace. It goes beyond differentially expressed genes and infers splicing isoform network changes between normal and diseased samples. It eases the sample size bottleneck; evaluations on simulated data and lung cancer-specific ERBB2 and MAPK signaling pathways, with varying number of samples, evince the merit in handling high exon to sample size ratio datasets. Inferred network rewiring of well established Bcl-x and EGFR centered networks from lung adenocarcinoma expression data is in good agreement with literature. Gene level evaluations demonstrate a substantial performance of SpliceNet over canonical correlation analysis, a method that is currently applied to exon level RNA-Seq data. SpliceNet can also be applied to exon array data. SpliceNet is distributed as an R package available at http://www.jjwanglab.org/SpliceNet.


2005 ◽  
Vol 20 (3) ◽  
pp. 224-232 ◽  
Author(s):  
Ioannis M. Stylianou ◽  
Michael Clinton ◽  
Peter D. Keightley ◽  
Clare Pritchard ◽  
Zuzzana Tymowska-Lalanne ◽  
...  

Obesity-related diseases are poised to become the primary cause of death in developed nations. While a number of monogenic causes of obesity have recently been identified, these are responsible for only a small proportion of human cases of obesity. Quantitative trait locus (QTL) studies using animal models have revealed hundreds of potential loci that affect obesity; however, few have been further analyzed beyond the original QTL scan. We previously mapped four QTL in an F2 between divergently selected Fat (F) and Lean (L) lines. A QTL of large effect on chromosome 15 ( Fob3) was subsequently mapped to a higher resolution into two smaller-effect QTL ( Fob3a and Fob3b) using crosses between the F-line and a congenic line containing L-line alleles at the Fob3 QTL region. Here we report the gene expression characterization of Fob3b. Microarray expression analysis using the NIA-NIH 15K cDNA array set containing 14,938 mouse ESTs was employed to identify candidate genes and pathways that are differentially expressed between the F-line and a congenic line containing only the Fob3b QTL ( Fob3b-line). Our study suggests squalene epoxidase (Sqle), a cholesterol biosynthesis enzyme, as a strong positional candidate gene for Fob3b. Several other cholesterol biosynthesis pathway genes unlinked to Fob3b were found to be differentially expressed, suggesting that a perturbation of this pathway could be in part responsible for the phenotypic difference between the F-line and Fob3b-line mice.


Sign in / Sign up

Export Citation Format

Share Document