Transcriptome profile in bursa of Fabricius reveals potential mode for stress-influenced immune function in chicken stress model.

BACKGROUND: The molecular mechanisms underlying stress-influenced immune function of chicken (Gallus Gallus) are not clear. The stress models can be established effectively by feeding chickens corticosterone (CORT) hormone. The bursa of Fabricius is a unique central immune organ of birds. RNA-Seq technology was used to investigate differences in the expression profiles of immune-related genes and associated pathways in the bursa of Fabricius to clarify molecular mechanisms. The aim of this study was to broaden the understanding of the stress-influenced immune function in chickens. RESULTS: Differentially expressed genes (DEGs) in the bursa of Fabricius between experimental group (basal diet with added CORT 30 mg/kg; C_B group) and control group (basal diet; B_B group) were identified by using RNA-seq technology. In total, we found 1434 significant DEGs (SDEGs), which included 199 upregulated and 1235 downregulated genes in the C_B group compared with the B_B group. The immune system process GO term was the top significantly GO term, including MYD88, TLR4, IL15, VEGFA gene and so on. The cytokine-cytokine receptor interaction pathway and the Toll-like receptor signaling pathway were the key pathways affected by stress. The protein-protein interaction (PPI) analysis of the SDEGs showed that VEGFA, MyD88 and IL15 were hub genes and module analysis showed that MYD88, TLR4 and VEGFA play important roles in response to stress. CONCLUSION: This study showed that the VEGFA and ILs (such as IL15) via the cytokine-cytokine receptor interaction pathway, MYD88 and TLR4 via the Toll-like receptor signaling pathway may play important roles in the regulation of immune function under stress condition with CORT administration. The results of this study provide a reference for further studies of the molecular mechanisms of stress-influenced immune function.

Effects of immunosuppression-associated gga-miR-146a-5p on immune regulation in chicken macrophages by targeting the IRKA2 gene.

Stress-induced immunosuppression (SIIS) is one of the common problems in intensive poultry production, which brings enormous economic losses to the poultry industry. Accumulating evidence has shown that microRNAs (miRNAs) were important regulators of gene expression in the immune system. However, the miRNA-mediated molecular mechanisms underlying SIIS in chickens are still poorly understood. This study aimed to investigate the biological functions and regulatory mechanism of miRNAs in chicken SIIS. A stress-induced immunosuppression model was successfully established via daily injection of dexamethasone and analyzed miRNA expression in spleen. Seventy-four differentially expressed miRNAs (DEMs) was identified, and 229 target genes of the DEMs were predicted. Functional enrichment analysis the target genes revealed pathways related to immunity, such as MAPK signaling pathway and FoxO signaling pathway. The candidate miRNA, gga-miR-146a-5p, was found to be significantly downregulated in the Dex-induced chicken spleen, and we found that Dex stimulation significantly inhibited the expression of gga-miR-146a-5p in Chicken macrophages (HD11). Flow cytometry, 5-ethynyl-2'-deoxyuridine (EdU), cell counting kit-8 (CCK-8) and other assays indicated that gga-miR-146a-5p can promote the proliferation and inhibit apoptosis of HD11 cells. A dual-luciferase reporter assay suggested that the Interleukin 1 receptor associated kinase 2 (IRAK2) gene, which encoded a transcriptional factor, was a direct target of gga-miR-146a-5p, gga-miR-146a-5p suppressed the post-transcriptional activity of IRAK2. These findings not only improve our understanding of the specific functions of miRNAs in avian stress but also provide potential targets for genetic improvement of stress resistance in poultry.

Analysis of miRNA and mRNA reveals core interaction networks and pathways of dexamethasone-induced immunosuppression in chicken bursa of Fabricius.

Stress-induced immunosuppression is a serious problem affecting the production value of poultry, but its specific molecular mechanism has not yet been elucidated. We selected 7-day-old Gushi cocks as test animals and successfully established a stress-induced immunosuppression model by injecting 2.0 mg/kg (body weight) dexamethasone (Dex). We then constructed six cDNA libraries and two small RNA libraries of Bursa of Fabricius from the control group and the Dex group. RNA-seq results revealed 21,028 transcripts including 3920 novel transcripts; 500 miRNAs including 68 novel miRNAs were identified. Correlation analysis of miRNA, target genes and mRNA results indicated that the gga-miR-15 family, gga-miR-103-3p, gga-miR-456-3p, and gga-miR-27b-3p, as core differentially expressed miRNAs, may potentially regulate multiple genes which are involved in immune-related pathways; and that the core genes Suppressor of IKBKE 1 (SIKE1) and high mobility group AT-hook 2 (HMGA2) are associated with the miR-17 family (gga-miR-20a-5p, gga-miR-20b-5p, gga-miR-106-5p, and gga-miR-17-5p) and gga-let -7 family (gga-let-7b, gga-let-7i, gga-let-7c-5p, and gga-let-7f-5p). The interaction networks of mRNAs of significantly enrichment pathways and PPI (protein-protein interaction) networks showed that IL6, IL1B, IL8L1, CCL5, SOCS3, SOCS1, ITGB5, GSTA3, SQLE, FDFT1, FN1, IL18, IL10, MAPK11 and MAPK12 are network core nodes and that most of them are strongly associated with immune response. One of the candidate miRNAs, gga-miR-20b-5p, may play an important role in stress-induced immunosuppression. Luciferase assay and over-expression experiments suggested that gga-miR-20b-5p negatively regulated the expression of target gene SIKE1. These results provide better understanding of the mechanism of stress-induced immunosuppression in Gushi chicken bursa, and provide novel targets for subsequent research to improve poultry anti-stress capability.

Identification and expression analysis of MicroRNAs in chicken spleen in a corticosterone-induced stress model.

For investigating the effects of stress on the immune response of chickens, we established a corticosterone (CORT)-induced stress model by exogenous intake of CORT. Control group was fed with a basal diet and the stress model group was fed with a 30 mg/Kg CORT-treated diet in ad libitum conditions for 7 days. Then, we used RNA-seq technology to identify the expression pattern of miRNAs, target genes, and relevant pathways in chicken spleen. Results showed that 71 differentially expressed miRNAs (DEMs) were determined, 9 of which were significantly differentially expressed miRNAs (SDEMs), and 241 target genes of DEMs were predicted. GO annotation and KEGG pathway analysis were carried out to understand the role of the DEMs. Out of 287 significantly enriched GO terms, 37 were stress- or immune-related, such as response to light stimulus, detection of oxidative stress, and immune response in mucosal-associated lymphoid tissue. Out of 85 KEGG pathways, 8 were related to stress or immunity, such as cytokine-cytokine receptor interaction, JAK-STAT signaling pathway, and RLR signaling pathway. We then constructed the interaction networks between target genes from immune-related pathways and their DEMs. The analysis results suggested that some DEMs (gga-miR-17 family, gga-miR-15/16 family, gga-miR-2954 and gga-miR-34b-5p) and target genes (SIKE1, CX3CL1, IL11Ra, PIGR, and CDKN1A) were core miRNAs and genes. This study revealed the dynamic miRNA transcriptome, target genes and related pathways in chicken spleen under CORT-induced stress model, which provided a basis for studying the molecular mechanism of stress affecting immune function.

Identification of genes related to stress affecting thymus immune function in a chicken stress model using transcriptome analysis.

With the rapid development of the poultry breeding industry and highly intensive production management, the losses caused by stress responses are becoming increasingly serious. To screen candidate genes related to chicken stress and provide a basis for future research on the molecular mechanisms governing the effects of stress on chicken immune function, we successfully constructed a chicken stress model by exogenously introducing corticosterone (CORT). RNA-seq technology was used to identify and analyze the mRNA and enrichment pathways of the thymus in the stress model group and the control group. The results showed that there were 101 significantly differentially expressed genes (SDEGs) (Padj < 0.05, |log2fold changes| ≥ 1 and FPKM >1), of which 44 were upregulated genes, while 57 were downregulated genes. Gene Ontology (GO) enrichment analysis found that the terms related to immunity or stress mainly included antigen processing and presentation, positive regulation of T cell-mediated immunity, and immune effector process. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis showed that the main pathways related to immunity or stress were the PPAR signaling pathway, NOD-like receptor signaling pathway, and intestinal immune network for IgA production. Among the SDEGs, XCL1, HSPA8, DMB1 and BAG3 are strongly related to immunity or stress and may be important genes involved in regulating stress affecting the immune function of chickens. The above results provide a theoretical reference for subsequent research on the molecular regulatory mechanisms by which stress affects the immune function of poultry.

TMT-based quantitative proteomic analysis reveals the spleen regulatory network of dexamethasone-induced immune suppression in chicks.

Stress-induced immunosuppression is one of the most widespread problems in the poultry industry. Understanding the molecular regulatory mechanism of immunosuppression induced by stress in the chicken spleen would provide a scientific foundation for the prevention of stress reactions and antistress molecular breeding in poultry. To assess the protein expression profile of spleen tissue in a stress-included immunosuppression model, we performed a TMT-based proteomic analysis of chicken spleen tissue in a Dex-induced immunosuppression model (group C) and a control group (group A). We identified 590 differentially abundant proteins (DAPs) in chicken spleen tissue. These DAPs were significantly enriched in the following functional categories: ECM-receptor interaction, DNA replication, p53 signaling pathway, PI3K-Akt signaling pathway and NF-kappa B signaling pathway. Integrative analysis of the proteome and our previous transcriptome data revealed 62 DAPs showing correlations with the expression of their encoding mRNAs. Complementary proteome- and transcriptome-level analyses revealed a complex molecular network of stress-included immunosuppression. DPP4 and ALDH1A3 were the most significantly upregulated DAPs. GBP and OASL were identified as important nodes in the network related to stress-induced immunosuppression. The candidate genes identified in this study may be useful for the marker-based breeding of new chicken varieties with reduced stress levels. SIGNIFICANCE: This study provides a large amount of new information about the spleen proteome of the Dex-induced immunosuppression in chicks, as well as the correlation of transcriptome and proteome. Analysis of this resource has enabled us to examine mechanism of protein and transcript diversification, which expands the understanding of the complexity of the mechanism of stress-induced immunosuppression.

Use of transcriptomic analysis to identify microRNAs related to the effect of stress on thymus immune function in a chicken stress model.

In modern poultry production, stress-induced immunosuppression leads to serious economic losses and harm to animals, but the molecular mechanisms governing the effects of stress on the chicken thymus have not been elucidated. In this study, we successfully constructed a stress model of 7-day-old Gushi chickens by adding exogenous corticosterone (CORT) to their diet and determined the microRNA (miRNA) expression profile of thymus tissues using RNA-seq technology. The results identified 51 differentially expressed miRNAs (DEMs), including 30 upregulated miRNAs and 21 downregulated miRNAs. A total of 164 target genes of the DEMs were predicted based on bioinformatic analysis methods, and Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses of these target genes were performed. The results from the GO enrichment analysis of the target genes identified 349 significantly enriched terms, including terms associated with the stress response and immune function that are primarily involved in the negative regulation of phagocytosis, the response to stress and the cellular response to stimulus. The KEGG pathway analysis indicated that the enriched pathways related to immunity or stress included the MAPK signaling pathway, lysosomes, endocytosis, and the RIG-I-like receptor signaling pathway. Among these pathways, DEMs (such as gga-miR-2954, gga-miR-106-5p, and gga-miR-16-5p) and corresponding target genes (such as IL11Ra, SIKE1, and CX3CL1) might be strongly correlated with thymic immunity in chickens. The results of this study provide a reference for further research on the molecular regulatory mechanisms governing the effect of stress on the immune function of the chicken thymus.

Identification of genes related to dexamethasone-induced immunosuppression in chicken thymus using transcriptome analysis.

The molecular mechanism of stress-induced immunosuppression (SIS) in certain poultry immune organs is not completely clear. In this study, we constructed a stress immunosuppression model by selecting 180 healthy 7-day-old Gushi chickens and dividing them randomly into two groups: a D_T group and a B_T group. The D_T group was given dexamethasone, and the B_T group was given normal saline, according to the treatment method established and reported in our previous study. Thymus samples were subsequently taken from both groups. RNA-seq was used to sequence the transcriptomes of the thymus samples from both groups, and 1278 significant differentially expressed genes (DEGs) were obtained, of which 845 genes were up-regulated and 433 genes were down-regulated (padj<0.05, |FC| ≥ 2, FPKM>1). We identified immune-related gene ontology (GO) terms including immune system processes, immune system process regulation, and T cell activation. The results of KEGG (http: //www.kegg.jp) analysis showed that the DEGs are involved in a variety of immune-related pathways, such as cytokine-cytokine receptor interactions, Jak-STAT signaling pathways, and cell adhesion molecules (CAMs). The cytokine-cytokine receptor interaction pathway involves the DEGs CCR6, CCR5, CD40LG and FAS. The DEGs in the Jak-STAT signaling pathway were SPRY2, BCL2L1. These DEGS play an important role in cell apoptosis. CD40L, CD8, among other genes, are involved in the CAMs pathway. The results of this study add to existing data on the genomic study of stress affecting immune function, and provide a basis for further studies of the molecular mechanisms of stress-influenced immune function.

Identification of genes related to effects of stress on immune function in the spleen in a chicken stress model using transcriptome analysis.

Stress is a physiological manifestation of the body's defense against adverse effects of external environment, but the molecular regulatory mechanism of stress effects on immune function of poultry has not been fully clarified. In this study, 7-day-old Chinese local breed Gushi cocks were used as model animal, and the stress model was successfully constructed by adding corticosterone (CORT) 30 mg/kg basic diet for 7 days. The spleen transcriptomes of the control group (B_S group) and the stress model group (C_S group) was determined by high-throughput mRNA sequencing (RNA-Seq) technology, and a total of 269 significantly differentially expressed genes (SDEGs) were obtained (Padj < 0.05, |FC| ≥ 2 and FPKM > 1). Compared with B_S group, there were 140 significantly up-regulated genes and 129 significantly down-regulated genes in C_S group. The immune/stress-related Gene Ontology (GO) terms included positive regulation of T cell mediated immunity, chemokine-mediated signaling pathway, T cell mediated immunity and so on. The SDEGs such as IL8L1, HSPA8, HSPA2, RSAD2, CCR8L and DMB1 were involved in these GO terms. Kyoto encyclopedia of genes and genomes (KEGG) enrichment analysis showed that the SDEGs participated in many immune-related signaling pathways. The immune-related genes HSPA2, HSPA8, HSP90AA1, HSPH1 and HERPUD1 were enriched in Protein processing in endoplasmic reticulum pathway, IL8L1, CXCL13L2, CCR6, LEPR, CCR9 and CCR8L were enriched in Cytokine-cytokine receptor interaction pathway. The protein-protein interactions (PPI) analysis showed HSPA8, HSPA2 and IL8L1 as key core nodes had 7 interactions and may play important roles in the regulation of CORT-induced stress effects on immune function. The data onto this study enriched the genomic study of stress effects on immune function, and provided unique insights into the molecular mechanism of stress effects on immune function, and the genes identified in this study can be candidates for future research on stress response.

Transcriptomic Analysis of Spleen Revealed Mechanism of Dexamethasone-Induced Immune Suppression in Chicks.

Stress-induced immunosuppression is a common problem in the poultry industry, but the specific mechanism of its effect on the immune function of chicken has not been clarified. In this study, 7-day-old Gushi cocks were selected as subjects, and a stress-induced immunosuppression model was successfully established via daily injection of 2.0 mg/kg (body weight) dexamethasone. We characterized the spleen transcriptome in the control (B_S) and model (D_S) groups, and 515 significant differentially expressed genes (SDEGs) (Fragments Per Kilobase of transcript sequence per Millions base pairs sequenced (FPKM) > 1, adjusted p-value (padj) < 0.05 and Fold change (|FC|) ≥ 2) were identified. The cytokine-cytokine receptor interaction signaling pathway was identified as being highly activated during stress-induced immunosuppression, including the following SDEGs-CXCL13L2, CSF3R, CSF2RB, CCR9, CCR10, IL1R1, IL8L1, IL8L2, GHR, KIT, OSMR, TNFRSF13B, TNFSF13B, and TGFBR2L. At the same time, immune-related SDEGs including CCR9, CCR10, DMB1, TNFRSF13B, TNFRSF13C and TNFSF13B were significantly enriched in the intestinal immune network for the IgA production signaling pathway. The SDEG protein-protein interaction module analysis showed that CXCR5, CCR8L, CCR9, CCR10, IL8L2, IL8L1, TNFSF13B, TNFRSF13B and TNFRSF13C may play an important role in stress-induced immunosuppression. These findings provide a background for further research on stress-induced immunosuppression. Thus, we can better understand the molecular genetic mechanism of chicken stress-induced immunosuppression.