scholarly journals RNA seq analyses of chicken reveals biological pathways involved in acclimation into different geographical locations

2020 ◽  
Vol 10 (1) ◽  
Author(s):  
Himansu Kumar ◽  
Hyojun Choo ◽  
Asankadyr U. Iskender ◽  
Krishnamoorthy Srikanth ◽  
Hana Kim ◽  
...  

Abstract Transcriptome expression reflects genetic response in diverse conditions. In this study, RNA sequencing was utilized to profile multiple tissues such as liver, breast, caecum, and gizzard of Korean commercial chicken raised in Korea and Kyrgyzstan. We analyzed ten samples per tissue from each location to identify candidate genes which are involved in the adaptation of Korean commercial chicken to Kyrgyzstan. At false discovery rate (FDR) < 0.05 and fold change (FC) > 2, we found 315, 196, 167 and 198 genes in liver, breast, cecum, and gizzard respectively as differentially expressed between the two locations. GO enrichment analysis showed that these genes were highly enriched for cellular and metabolic processes, catalytic activity, and biological regulations. Similarly, KEGG pathways analysis indicated metabolic, PPAR signaling, FoxO, glycolysis/gluconeogenesis, biosynthesis, MAPK signaling, CAMs, citrate cycles pathways were differentially enriched. Enriched genes like TSKU, VTG1, SGK, CDK2 etc. in these pathways might be involved in acclimation of organisms into diverse climatic conditions. The qRT-PCR result also corroborated the RNA-Seq findings with R2 of 0.76, 0.80, 0.81, and 0.93 for liver, breast, caecum, and gizzard respectively. Our findings can improve the understanding of environmental acclimation process in chicken.

Genes ◽  
2020 ◽  
Vol 11 (6) ◽  
pp. 631
Author(s):  
Zhicheng Sun ◽  
Fangrui Lou ◽  
Yuan Zhang ◽  
Na Song

Acanthogobius ommaturus is a euryhaline fish widely distributed in coastal, bay and estuarine areas, showing a strong tolerance to salinity. In order to understand the mechanism of adaptation to salinity stress, RNA-seq was used to compare the transcriptome responses of Acanthogobius ommaturus to the changes of salinity. Four salinity gradients, 0 psu, 15 psu (control), 30 psu and 45 psu were set to conduct the experiment. In total, 131,225 unigenes were obtained from the gill tissue of A. ommaturus using the Illumina HiSeq 2000 platform (San Diego, USA). Compared with the gene expression profile of the control group, 572 differentially expressed genes (DEGs) were screened, with 150 at 0 psu, 170 at 30 psu, and 252 at 45 psu. Additionally, among these DEGs, Gene Ontology (GO) analysis indicated that binding, metabolic processes and cellular processes were significantly enriched. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways analysis detected 3, 5 and 8 pathways related to signal transduction, metabolism, digestive and endocrine systems at 0 psu, 30 psu and 45 psu, respectively. Based on GO enrichment analysis and manual literature searches, the results of the present study indicated that A. ommaturus mainly responded to energy metabolism, ion transport and signal transduction to resist the damage caused by salinity stress. Eight DEGs were randomly selected for further validation by quantitative real-time PCR (qRT-PCR) and the results were consistent with the RNA-seq data.


Animals ◽  
2019 ◽  
Vol 9 (12) ◽  
pp. 1076 ◽  
Author(s):  
Himansu Kumar ◽  
Asankadyr U. Iskender ◽  
Krishnamoorthy Srikanth ◽  
Hana Kim ◽  
Asankadyr T. Zhunushov ◽  
...  

RNA sequencing was used to profile the liver transcriptome of a Korean commercial chicken (Hanhyup) at two different environments (Korea and Kyrgyzstan) to investigate their role during acclimatization into different climatic conditions. Ten samples from each location were analyzed to identify candidate genes that respond to environmental changes such as altitude, humidity, temperature, etc. Sequencing reads were preprocessed, aligned with the reference genome, assembled and expressions were estimated through bioinformatics approaches. At a false discovery rate (FDR) <0.05 and fold change (FC) ≥2, we found 315 genes were DE. Out of 315 DE genes, 174 and 141 were up- and down-regulated respectively in the Kyrgyz environment. Gene ontology (GO) enrichment analysis showed that the differentially expressed genes (DEGs) were associated with energy metabolism such as pyruvate and lactate metabolic processes, and glycerol catabolic process. Similarly, KEGG pathway analysis indicated pyruvate metabolism, glycolysis/gluconeogenesis, biosynthesis, citrate cycles were differentially enriched in the Kyrgyz environment. DEGs like TSKU, VTG1, SGK, CDK2, etc. in such pathways are highly involved in the adaptation of organisms into diverse climatic conditions. Our investigation may serve as a resource for the chicken industry, especially in exporting Hanhyup chicken from Korea to other countries.


2020 ◽  
Vol 48 (07) ◽  
pp. 1633-1650
Author(s):  
Xian Zhang ◽  
Xiaoxuan Zhao ◽  
Kaili Liu ◽  
Yuxuan Che ◽  
Xun Qiu ◽  
...  

Bufalin is an anticancer drug extract from traditional Chinese medicine. Several articles about bufalin have been published. However, the literature on bufalin has not yet been systematically studied. This study aimed to identify the study status and knowledge structures of bufalin and to summarize the antitumor mechanism. Data were retrieved and downloaded from the PubMed database. The softwares of BICOMB, gCLUTO, Ucinet 6.0, and NetDraw2.084 were used to analyze these publications. The bufalin related genes were recognized and tagged by ABNER software. Then these BF-related genes were performed by Gene Ontology (GO) enrichment analysis, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways analysis, and protein-protein interaction (PPI) network analysis. A total of 474 papers met the search criteria from 2000 to 2019. By biclustering clustering analysis, the 50 high-frequency main MeSH terms/subheadings were classified into 5 clusters. The clusters of drug therapy and the mechanism of bufalin were hotspot topics. A total of 50 genes were identified as BF-related genes. PPI network analysis showed that inducing apoptosis was the main effect of bufalin, and apoptosis-related gene Caspase 3 was the most reported by people. Bufalin could inhibit the proliferation, invasion, and metastasis of cancer cells through multiple signaling pathways, such as PI3K/AKT, Hedgehog, MAPK/JNK, Wnt/[Formula: see text]-catenin, TGF-[Formula: see text]/Smad, Integrin signaling pathway, and NF-KB signaling pathway via KEGG analysis. Through the quantitative analysis of bufalin literature, we revealed the research status and hot spots in this field and provided some guidance for further research.


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.


Animals ◽  
2020 ◽  
Vol 10 (2) ◽  
pp. 344 ◽  
Author(s):  
Bryan Irvine Lopez ◽  
Kier Gumangan Santiago ◽  
Donghui Lee ◽  
Seungmin Ha ◽  
Kangseok Seo

Immune response of 107 vaccinated Holstein cattle was initially obtained prior to the ELISA test. Five cattle with high and low bovine viral diarrhea virus (BVDV) type I antibody were identified as the final experimental animals. Blood samples from these animals were then utilized to determine significant differentially expressed genes (DEGs) using the RNA-seq transcriptome analysis and enrichment analysis. Our analysis identified 261 DEGs in cattle identified as experimental animals. Functional enrichment analysis in gene ontology (GO) annotations and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways revealed the DEGs potentially induced by the inactivated BVDV type I vaccine, and might be responsible for the host immune responses. Our findings suggested that inactivated vaccine induced upregulation of genes involved in different GO annotations, including antigen processing and presentation of peptide antigen (via MHC class I), immune response, and positive regulation of interferon-gamma production. The observed downregulation of other genes involved in immune response might be due to inhibition of toll-like receptors (TLRs) by the upregulation of the Bcl-3 gene. Meanwhile, the result of KEGG pathways revealed that the majority of DEGs were upregulated and enriched to different pathways, including cytokine-cytokine receptor interaction, platelet activation, extracellular matrix (ECM) receptor interaction, hematopoietic cell lineage, and ATP-binding cassette (ABC) transporters. These significant pathways supported our initial findings and are known to play a vital role in shaping adaptive immunity against BVDV type 1. In addition, type 1 diabetes mellitus pathways tended to be significantly enriched. Thus, further studies are needed to investigate the prevalence of type 1 diabetes mellitus in cattle vaccinated with inactivated and live BVDV vaccine.


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.


Blood ◽  
2019 ◽  
Vol 134 (Supplement_1) ◽  
pp. 160-160
Author(s):  
Ahmad Faisal Karim ◽  
Anthony R. Soltis ◽  
Nadia P Ewing ◽  
Clifton L. Dalgard ◽  
Matthew D. Wilkerson ◽  
...  

The formation of pathological anti-FVIII antibodies, referred to as "inhibitors", is the most serious complication of therapeutic FVIII infusions, affecting up to one third of severe Hemophilia A (HA) patients. Intensive FVIII therapy, i.e. "Immune Tolerance Induction" (ITI), enables ~2/3 of treated patients to achieve peripheral tolerance to FVIII. FVIII inhibitor formation is a classical T-cell dependent adaptive immune response. As such, it requires help from the innate immune system. However, the roles of innate immune cells and mechanisms of inhibitor development versus immune tolerance, achieved with or without ITI therapy, are not well understood. To address these questions, we carried out temporal transcriptomics profiling of FVIII-stimulated peripheral blood mononuclear cells (PBMCs) from HA subjects with and without a current or historic inhibitor using RNA-seq. PBMCs were isolated from 40 subjects in the following groups: (A) HA with an inhibitor that resolved either following ITI or spontaneously; (B) HA with a current inhibitor; (C) HA with no inhibitor history and (D) non-HA healthy controls. PBMCs were rested overnight and then stimulated with 5 nM FVIII, and total RNA was isolated 4, 16, 24 and 48 hours following stimulation. RNA from unstimulated cells at t = 4 hrs served as a negative control. Time-series differential expression analysis was performed with DESeq2 and genes with a log likelihood ratio test FDR &lt;0.05 and an absolute fold change &gt;1.25 at at least one stimulation time point compared to control were deemed significant. Subjects with a resolved past inhibitor (Group A) showed differential expression of only 15 genes. In contrast, subjects with a current inhibitor (Group B) showed differential expression of 56 genes. A clustering analysis divided the temporal trajectories of Group B genes into 3 distinct clusters. Twenty-three genes were up-regulated at 16 hr and 21 genes at 48 hr post-stimulation, respectively. Interestingly, gene ontology (GO) enrichment analysis of these genes revealed enrichments for innate immune modulators, including NLRP3, TLR8, IL32, CLEC10A and COLEC12.NLRP3 and TLR8 are associated with enhanced secretion of the pro-inflammatory cytokines IL-1beta and TNF-alpha, while IL32, which has several isoforms, has been associated with both inflammatory and regulatory immune processes. Expression levels of NLRP3, TLR8, CLEC10A and IL32 transcripts were validated by real time PCR, and changes in RNA transcript abundances correlated well with the RNA-seq results. IL-32 results were validated by both RT-qPCR on an aliquot of the original RNA sample and ELISA to measure the cytokine in supernatants at t=48 hrs. HA subjects with no inhibitor history (Group C) had 195 differentially expressed genes whose temporal profiles fell into 4 distinct clusters. GO enrichment analysis revealed biological processes related to epithelial cell proliferation, responses to toxic substances, and positive/negative regulation of cytokine secretion (TNF, NQO1, PMEPA1). The non-HA healthy control subjects (Group D) also showed cellular responses to ex vivo FVIII stimulation. A total of 63 differentially regulated genes fell into 4 distinct clusters. GO analysis identified expression patterns associated with leukocyte-mediated immunity, T-cell activation, and a hypoxia response. Overall, distinct transcriptional signatures were identified for each of the four groups, providing clues as to cellular mechanisms leading to or accompanying their disparate anti-FVIII antibody responses. We are currently characterizing PBMC immune cell subsets, e.g. macrophages and CD4+ T cells, to identify specific cell types responsible for the differentially regulated genes. Cellular responses of tolerized HA subjects and healthy non-HA controls were consistent with the known immunogenicity of FVIII, including persistence of FVIII-specific CD4+ T cells even in individuals with no measurable FVIII inhibitor. The inflammatory status of HA patients suffering from an ongoing inhibitor clearly includes up-regulation of innate immune modulators, some of which may act as ongoing danger signals that influence the responses to, and eventual outcome of, ITI therapy. Disclosures Pratt: Grifols, Inc: Research Funding; Bloodworks NW: Patents & Royalties: inventor on patents related to FVIII immunogenicity.


2020 ◽  
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 1,360 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.


Author(s):  
Er-Hu Chen ◽  
Jin-Yan Duan ◽  
Wei Song ◽  
Dian-Xuan Wang ◽  
Pei-An Tang

Abstract The rusty grain beetle, Cryptolestes ferrugineus (Stephens), is a serious pest of stored grain, which has developed high levels of resistance to phosphine. In this study, five geographically distant populations of C. ferrugineus had been collected in China, specifically in granaries where phosphine fumigant is used for pest control, and they showed a high resistance ratio up to 1,907 (LC50 = 21.0 mg/liter). Then, a reference transcriptome was constructed to use as a basis for investigating the molecular mechanisms of phosphine resistance in this species, which consisted of 47,006 unigenes with a mean length of 1,090. Subsequently, the RNA-Seq analysis of individuals from the most susceptible and resistant populations led to the identification of 54 genes that are differentially expressed. GO and KEGG analysis demonstrated that genes associated with mitochondrial and respiration functions were significantly enriched. Also, the ‘structural constituent of cuticle’ term was annotated in the GO enrichment analysis and further qRT-PCR confirmed that the expression levels of nine cuticular protein genes were significantly increased in the resistant population. In conclusion, we present here a transcriptome-wide overview of gene expression changes between resistant and susceptible populations of C. ferrugineus, and this in turn documents that mitochondria and cuticular protein genes may play together a crucial role in phosphine resistance. Further gene function analysis should enable the provision of advice to expedite resistance management decisions.


2020 ◽  
Author(s):  
Xin Li ◽  
Chenxin Wang ◽  
Xiaoqing Zhang ◽  
Jiali Liu ◽  
Yu Wang ◽  
...  

Abstract Objective: To reveal the molecular mechanism underlying the pathogenesis of HCM and find new effective therapeutic strategies using a systematic biological approach.Methods: The WGCNA algorithm was applied to building the co-expression network of HCM samples. A sample cluster analysis was performed using the hclust tool and a co-expression module was constructed. The WGCNA algorithm was used to study the interactive connection between co-expression modules and draw a heat map to show the strength of interactions between modules. The genetic information of the respective modules was mapped to the associated GO terms and KEGG pathways, and the Hub Genes with the highest connectivity in each module were identified. The Wilcoxon test was used to verify the expression level of hub genes between HCM and normal samples, and the "pROC" R package was used to verify the possibility of hub genes as biomarkers. Finally, the potential functions of hub genes were analyzed by GSEA software. Results: Seven co-expression modules were constructed using sample clustering analysis. GO and KEGG enrichment analysis judged that the turquoise module is an important module. The hub genes of each module are RPL35A for module Black, FH for module Blue, PREI3 for module Brown, CREB1 for module Green, LOC641848 for module Pink, MYH7 for module Turquoise and MYL6 for module Yellow. The results of the differential expression analysis indicate that MYH7 and FH are considered true hub genes. In addition, the ROC curves revealed their high diagnostic value as biomarkers for HCM. Finally, in the results of the GSEA analysis, MYH7 and FH highly expressed genes were enriched with the "proteasome" and a "PPAR signaling pathway," respectively.Conclusions: The MYH7 and FH genes may be the true hub genes of HCM. Their respective enriched pathways, namely the "proteasome" and the "PPAR signaling pathway," may play an important role in the development of HCM.


Sign in / Sign up

Export Citation Format

Share Document