Transcription profiles of the responses of chicken bursae of Fabricius to IBDV in different timing phases

Integrative RNA-seq and ChIP-seq analysis unveils metabolic regulation as a conserved antiviral mechanism of chicken p53

The tumor suppressor p53, primarily functioning as a transcription factor, has exhibited antiviral capabilities against various viruses in chickens, including infectious bursal disease virus (IBDV), avian leukosis virus subgroup J (ALV-J), and avian infectious laryngotracheitis virus (ILTV). Nevertheless, the existence of a universal antiviral mechanism employed by chicken p53 (chp53) against these viruses remains uncertain. This study conducted a comprehensive comparison of molecular networks involved in chp53's antiviral function against IBDV, ALV-J, and ILTV. This was achieved through an integrated analysis of ChIP-seq data, examining chp53's genome-wide chromatin occupancy, and RNA-seq data from chicken cells infected with these viruses. The consistent observation of chp53 target gene enrichment in metabolic pathways, confirmed via ChIP-qPCR, suggests a ubiquitous regulation of host cellular metabolism by chp53 across different viruses. Further genome binding motif conservation analysis and transcriptional co-factor prediction suggest conserved transcriptional regulation mechanism by which chp53 regulates host cellular metabolism during viral infection. These findings offer novel insights into the antiviral role of chp53 and propose that targeting the virus-host metabolic interaction through regulating p53 could serve as a universal strategy for antiviral therapies in chickens.IMPORTANCEThe current study conducted a comprehensive analysis, comparing molecular networks underlying chp53's antiviral role against infectious bursal disease virus (IBDV), avian leukosis virus subgroup J (ALV-J), and avian infectious laryngotracheitis virus (ILTV). This was achieved through a combined assessment of ChIP-seq and RNA-seq data obtained from infected chicken cells. Notably, enrichment of chp53 target genes in metabolic pathways was consistently observed across viral infections, indicating a universal role of chp53 in regulating cellular metabolism during diverse viral infections. These findings offer novel insights into the antiviral capabilities of chicken p53, laying a foundation for the potential development of broad-spectrum antiviral therapies in chickens.

RNA-seq reveals role of cell-cycle regulating genes in the pathogenicity of a field very virulent infectious bursal disease virus

Infectious bursal disease virus (IBDV) infection causes highly contagious and immunosuppressive disease in poultry. The thymus, serving as the primary organ for T cell maturation and differentiation, plays an important role in the pathogenicity of IBDV in the infected chickens. However, there are no reports on the molecular pathogenesis of IBDV in the thymus currently. The aim of the study was to elucidate the molecular mechanisms underlying the pathogenicity of a field very virulent (vv) IBDV strain NN1172 in the thymus of SPF chickens using integrative transcriptomic and proteomic analyses. Our results showed that a total of 4,972 Differentially expressed genes (DEGs) in the thymus of NN1172-infected chickens by transcriptomic analysis, with 2,796 up-regulated and 2,176 down-regulated. Meanwhile, the proteomic analysis identified 726 differentially expressed proteins (DEPs) in the infected thymus, with 289 up-regulated and 437 down-regulated. Overall, a total of 359 genes exhibited differentially expression at both mRNA and protein levels, with 134 consistently up-regulated and 198 genes consistently down-regulated, as confirmed through a comparison of the RNA-seq and the proteomic datasets. The gene ontology (GO) analysis unveiled the involvement of both DEGs and DEPs in diverse categories encompassing cellular components, biological processes, and molecular functions in the pathological changes in IBDV-infected thymus. The Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis revealed that the host mainly displayed severely disruption of cell survival/repair, proliferation and metabolism pathway, meanwhile, the infection triggers antiviral immune activation with a potential emphasis on the MDA5 pathway. Network inference analysis identified seven core hub genes, which include CDK1, TYMS, MCM5, KIF11, CCNB2, MAD2L1, and MCM4. These genes are all associated with cell-cycle regulating pathway and are likely key mediators in the pathogenesis induced by NN1172 infection in the thymus. This study discovered dominant pathways and genes which enhanced our understanding of the molecular mechanisms underlying IBDV pathogenesis in the thymus.

Infectious bursal disease virus replication is inhibited by avain T cell chemoattractant chemokine CCL19

Chemokine CCL19, together with its receptor CCR7, is one of the most important factors recruiting immune cells into target organ during virus infection. Our previous study has shown that CCL19 played a vital role in the process of T cell trafficking into bursae during bursal disease virus (IBDV) infection. In this study, we hypothesized that CCL19 could exert direct influences on IBDV replication other than recruiting immune cells. A eukaryotic expression vector of pEGFP-N1/CCL19 was successfully constructed and identified by PCR, double enzymes digestion, and sequencing. Different concentrations of pEGFP-N1/CCL19 plasmids were transfected into DF1 cells and CCL19 protein was highly expressed. Then, DF1 cells were infected with IBDV B87 strain post-transfection. Based on PCR and Western blot results, CCL19 could obviously decrease the gene levels of VP1 and VP2 and the protein levels of VP2 and VP3. When CCL19 was knocked down, the gene levels of VP1 and VP2 were significantly upregulated. Moreover, indirect immunostaining revealed that the IBDV content was largely decreased after CCL19 overexpression. Additionally, CCL19 inhibitory effects might rely on activation of the JNK signal pathway. Taken together, chemokine CCL19 directly blocks IBDV replication in DF1 cells, indicating that CCL19 could play crucial functions other than recruiting T cells during the pathogenesis of IBDV.

Transcriptome Sequencing of the Spleen Reveals Antiviral Response Genes in Chickens Infected with CAstV

Astrovirus infections pose a significant problem in the poultry industry, leading to multiple adverse effects such as a decreased egg production, breeding disorders, poor weight gain, and even increased mortality. The commonly observed chicken astrovirus (CAstV) was recently reported to be responsible for the "white chicks syndrome" associated with an increased embryo/chick mortality. CAstV-mediated pathogenesis in chickens occurs due to complex interactions between the infectious pathogen and the immune system. Many aspects of CAstV-chicken interactions remain unclear, and there is no information available regarding possible changes in gene expression in the chicken spleen in response to CAstV infection. We aim to investigate changes in gene expression triggered by CAstV infection. Ten 21-day-old SPF White Leghorn chickens were divided into two groups of five birds each. One group was inoculated with CAstV, and the other used as the negative control. At 4 days post infection, spleen samples were collected and immediately frozen at -70 °C for RNA isolation. We analyzed the isolated RNA, using RNA-seq to generate transcriptional profiles of the chickens' spleens and identify differentially expressed genes (DEGs). The RNA-seq findings were verified by quantitative reverse-transcription PCR (qRT-PCR). A total of 31,959 genes was identified in response to CAstV infection. Eventually, 45 DEGs (p-value < 0.05; log2 fold change > 1) were recognized in the spleen after CAstV infection (26 upregulated DEGs and 19 downregulated DEGs). qRT-PCR performed on four genes (IFIT5, OASL, RASD1, and DDX60) confirmed the RNA-seq results. The most differentially expressed genes encode putative IFN-induced CAstV restriction factors. Most DEGs were associated with the RIG-I-like signaling pathway or more generally with an innate antiviral response (upregulated: BLEC3, CMPK2, IFIT5, OASL, DDX60, and IFI6; downregulated: SPIK5, SELENOP, HSPA2, TMEM158, RASD1, and YWHAB). The study provides a global analysis of host transcriptional changes that occur during CAstV infection in vivo and proves that, in the spleen, CAstV infection in chickens predominantly affects the cell cycle and immune signaling.

Bacillus amyloliquefaciens TL Downregulates the Ileal Expression of Genes Involved in Immune Responses in Broiler Chickens to Improve Growth Performance

Bacillus amyloliquefaciens TL promotes broiler chicken performance by improving nutrient absorption and utilization and reducing intestinal inflammation. In this study, RNA-sequencing (RNA-seq)-based transcriptomes of ileal tissues collected from probiotic-fed and control broiler chickens were analyzed to elucidate the effects of the probiotic B. amyloliquefaciens TL, as a feed additive, on the gut immune function. In total, 475 genes were significantly differentially expressed between the ileum of probiotic-fed and control birds. The expression of genes encoding pyruvate kinase, prothymosin-α, and heat stress proteins was high in the ileum of probiotic-fed birds (FPKM > 500), but not in the control group. The gene ontology functional enrichment and pathway enrichment analyses revealed that the uniquely expressed genes in the control group were mostly involved in immune responses, whereas those in the probiotic group were involved in fibroblast growth factor receptor signaling pathways and positive regulation of cell proliferation. Bacillus amyloliquefaciens TL downregulated the expression of certain proinflammatory factors and affected the cytokine–cytokine receptor interaction pathway. Furthermore, B. amyloliquefaciens TL in broiler diets altered the expression of genes involved in immune functions in the ileum. Thus, it might contribute to improved broiler growth by regulating the immune system and reducing intestinal damage in broilers.

Transcriptome analysis of chicken intraepithelial lymphocyte natural killer cells infected with very virulent infectious bursal disease virus

The infectious bursal disease (IBD) is an acute immunosuppressive viral disease that significantly affects the economics of the poultry industry. The IBD virus (IBDV) was known to infect B lymphocytes and activate macrophage and T lymphocytes, but there are limited studies on the impact of IBDV infection on chicken intraepithelial lymphocyte natural killer (IEL-NK) cells. This study employed an mRNA sequencing approach to investigate the early regulation of gene expression patterns in chicken IEL-NK cells after infection with very virulent IBDV strain UPM0081. A total of 12,141 genes were expressed in uninfected chicken IEL-NK cells, and most of the genes with high expression were involved in the metabolic pathway, whereas most of the low expressed genes were involved in the cytokine-cytokine receptor pathway. A total of 1,266 genes were differentially expressed (DE) at 3 day-post-infection (dpi), and these DE genes were involved in inflammation, antiviral response and interferon stimulation. The innate immune response was activated as several genes involved in inflammation, antiviral response and recruitment of NK cells to the infected area were up-regulated. This is the first study to examine the whole transcriptome profile of chicken NK cells towards IBDV infection and provides better insight into the early immune response of chicken NK cells.

Genome-wide analysis of differentially expressed mRNAs, lncRNAs, and circRNAs in chicken bursae of Fabricius during infection with very virulent infectious bursal disease virus

Background

Infectious bursal disease virus (IBDV) causes acute, highly contagious, immunosuppressive, and lethal infectious disease in young chickens and mainly infects the bursa of Fabricius (BF). To investigate interactions between IBDV and its host, RNA sequencing was applied to analyze the responses of the differentially expressed transcriptional profiles of BF infected by very virulent IBDV (vvIBDV).

Results

In total, 317 upregulated and 94 downregulated mRNAs were found to be significantly differentially expressed in infected chickens, compared to controls. Long non-coding RNA (lncRNA) and circular RNA (circRNA) alterations were identified in IBDV-infected chickens, and significantly different expression was observed in 272 lncRNAs and 143 circRNAs. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analyses were performed to assess the functions of significantly dysregulated genes, which showed that the JAK-STAT signaling pathway, the NOD-like receptor signaling pathway, and apoptosis may be activated by IBDV infection. We predicted interactions between differentially expressed genes and produced lncRNA-mRNA and circRNA-miRNA-mRNA regulator network.

Conclusions

The present study identified the expression profiles of mRNAs, lncRNAs, and circRNAs during vvIBDV infection and provides new insights into the pathogenesis of IBDV and antiviral immunity of the host.

Integrating Genetic and Genomic Analyses of Combined Health Data Across Ecotypes to Improve Disease Resistance in Indigenous African Chickens

Poultry play an important role in the agriculture of many African countries. The majority of chickens in sub-Saharan Africa are indigenous, raised in villages under semi-scavenging conditions. Vaccinations and biosecurity measures rarely apply, and infectious diseases remain a major cause of mortality and reduced productivity. Genomic selection for disease resistance offers a potentially sustainable solution but this requires sufficient numbers of individual birds with genomic and phenotypic data, which is often a challenge to collect in the small populations of indigenous chicken ecotypes. The use of information across-ecotypes presents an attractive possibility to increase the relevant numbers and the accuracy of genomic selection. In this study, we performed a joint analysis of two distinct Ethiopian indigenous chicken ecotypes to investigate the genomic architecture of important health and productivity traits and explore the feasibility of conducting genomic selection across-ecotype. Phenotypic traits considered were antibody response to Infectious Bursal Disease (IBDV), Marek's Disease (MDV), Fowl Cholera (PM) and Fowl Typhoid (SG), resistance to Eimeria and cestode parasitism, and productivity [body weight and body condition score (BCS)]. Combined data from the two chicken ecotypes, Horro (n = 384) and Jarso (n = 376), were jointly analyzed for genetic parameter estimation, genome-wide association studies (GWAS), genomic breeding value (GEBVs) calculation, genomic predictions, whole-genome sequencing (WGS), and pathways analyses. Estimates of across-ecotype heritability were significant and moderate in magnitude (0.22-0.47) for all traits except for SG and BCS. GWAS identified several significant genomic associations with health and productivity traits. The WGS analysis revealed putative candidate genes and mutations for IBDV (TOLLIP, ANGPTL5, BCL9, THEMIS2), MDV (GRM7), SG (MAP3K21), Eimeria (TOM1L1) and cestodes (TNFAIP1, ATG9A, NOS2) parasitism, which warrant further investigation. Reliability of GEBVs increased compared to within-ecotype calculations but accuracy of genomic prediction did not, probably because the genetic distance between the two ecotypes offset the benefit from increased sample size. However, for some traits genomic prediction was only feasible in across-ecotype analysis. Our results generally underpin the potential of genomic selection to enhance health and productivity across-ecotypes. Future studies should establish the required minimum sample size and genetic similarity between ecotypes to ensure accurate joint genomic selection.

Identification of Reference Genes in Chicken Intraepithelial Lymphocyte Natural Killer Cells Infected with Very-virulent Infectious Bursal Disease Virus

Due to the limitations in the range of antibodies recognising avian viruses, quantitative real-time PCR (RT-qPCR) is still the most widely used method to evaluate the expression of immunologically related genes in avian viruses. The objective of this study was to identify suitable reference genes for mRNA expression analysis in chicken intraepithelial lymphocyte natural killer (IEL-NK) cells after infection with very-virulent infectious bursal disease virus (vvIBDV). Fifteen potential reference genes were selected based on the references available. The coefficient of variation percentage (CV%) and average count of these 15 genes were determined by NanoString technology for control and infected samples. The M and V values for shortlisted reference genes (ACTB, GAPDH, HMBS, HPRT1, SDHA, TUBB1 and YWHAZ) were calculated using geNorm and NormFinder. GAPDH, YWHAZ and HMBS were the most stably expressed genes. The expression levels of three innate immune response related target genes, CASP8, IL22 and TLR3, agreed in the NanoString and RNA sequencing (RNA-Seq) results using one or two reference genes for normalisation (not HMBS). In conclusion, GAPDH and YWHAZ could be used as reference genes for the normalisation of chicken IEL-NK cell gene responses to infection with vvIBDV.

Differential expression of pro-inflammatory and anti-inflammatory genes of layer chicken bursa after experimental infection with infectious bursal disease virus

Infectious bursal disease (IBD) is one of the most prevalent infectious diseases caused by IBD virus (IBDV), which results in bursal necrosis and immunosuppression that cause severe damage to the immune system in chickens. Cytokines are important mediators and regulators of both types of host responses. In the present study, layer chickens were artificially challenged with IBDV, and the differential expression of inflammatory genes was explored by using quantitative real-time PCR, which offered basic data for further study of IBDV pathogenesis. Data showed that after IBDV infection, the virus load in the bursa of Fabricius (BF) peaked at 96 h and then gradually decreased. Compared with those of the negative-infected group, the mRNA expression levels of pro-inflammatory cytokines (interleukin [IL]-1β, IL-6, IL-7, IL-8, tumor necrosis factor [TNF]-α, transforming growth factor [TGF]-β) and anti-inflammatory cytokine IL-10 in the infected group increased to varying degrees at 12 to 192 h, respectively. Furthermore, the IL-1β mRNA expression peaked at 48 h; the mRNA transcript levels of IL-6, IL-8, and IL-10 were the highest at 96 h; TNF-α mRNA expression peaked at 120 h; the IL-7 mRNA expression peaked at 144 h; and the TGF-β mRNA transcript level was the highest at 192 h. Taken together, these observations indicated that along with the change pattern of IBDV proliferation in BF, the mRNA expression of cytokines (IL-1β, IL-6, IL-7, IL-8, IL-10, TNF-α, TGF-β) obviously increased, and the kinetics of each of these cytokines was different. The kinetics of IL-6/IL-10 mRNA expression ratio was significantly positively correlated with that of the virus load. These results suggest that IBDV infection seriously interferes with the natural immune response mediated by inflammatory cytokines in chickens.

Bioinformatic identification and expression analysis of the chicken B cell lymphoma (BCL) gene

BACKGROUND

B cell lymphoma (BCL) families play an important role in apoptosis as a growth factor, cell death programming, cytokine expression and immune-related genes expression.

OBJECTIVES

In this study, to investigate the roles of BCLs, we performed genome-wide identification, expression and functional analyses of the BCL family in chicken.

METHODS

Chicken BCLs genes were identified and analyzed by using bioinformatics approach. Expression profiles and Hierarchical cluster analysis of the BCLs genes in different chicken tissues were obtained from the genome-wide RNA-seq in the GEO, and Cluster and Java Treeview, respectively.

RESULTS

A total of 16 BCLs genes were identified from the chicken genome, which could be further classified into five distinct groups in the phylogenetic tree. On the other hand, the interaction among BCLs proteins and between BCLs proteins with NF-κB subunits are limited, indicating that the remaining the functions of BCLs protein could be investigated in chicken. Moreover, KEGG pathway analysis indicated that BCL gene family was involved in regulation of apoptotic and immune response. Finally, BCL gene family was differentially expressed in chicken tissues, pathogen infection and growth stages of early chicken early embryo.

CONCLUSION

This study provides significant insights into the potential functions of BCLs in chicken, including the regulation of apoptosis, cell death and expression of immune-related genes.

CC chemokine ligand 19 might act as the main bursal T cell chemoattractant factor during IBDV infection

Infectious bursa disease virus (IBDV) pathogenesis is characterized by increased numbers of T cells and decreased numbers of B cells in the bursa. Currently, little is about the key factor that affects T migration into bursa. In humans, CC chemokine ligand 19 (CCL19) recruits monocytes and neutrophils and is usually involved in various inflammatory disorders. The aim of this study was to assess the roles of CCL19 in driving peripheral blood cells infiltration into bursa of Fabricius of chickens infected with IBDV. Bursal samples were collected from chickens of the infection group and the control group on day 1, 3, 5, and 7 post infection (dpi) with IBDV. The mRNA or protein levels of ccl19 and ccr7 genes in bursae were determined by real-time PCR and immunohistochemistry (IHC) methods. Moreover, an in vitro chemotaxis assay was performed to evaluate the chemotaxis ability of CCL19 and bursal total protein. The results have displayed that the mRNA levels of ccl19 were significantly increased on 1, 3, 5, and 7 dpi in the infection group. The highest value amounted to 73.4-fold of the control group. Also, the mRNA levels of CCR7, the receptor of CCL19, began to increase on 3 dpi and reached to the highest value of 206.3-fold on 5 dpi after IBDV infection. Then the gene expression of CCR7 in bursae of the infection group returned to the normal level. IHC results of CCL19 protein level accorded with the mRNA levels of CCL19, with the highest value on 5 dpi. Then, in vitro chemotaxis test demonstrated that the total bursal protein had the ability of recruiting peripheral white blood cells (PWBC) and the migration percentage was a little higher than that of the blank control with only basal medium (P < 0.05). Taken together, these data suggest that CCL19 acts as a chicken PWBC chemotactic factor and facilitate the infiltration of PWBC (especially T cells) into the bursae after IBDV infection.

An Ex Vivo Chicken Primary Bursal-cell Culture Model to Study Infectious Bursal Disease Virus Pathogenesis

Infectious bursal disease virus (IBDV) is a birnavirus of economic importance to the poultry industry. The virus infects B cells, causing morbidity, mortality, and immunosuppression in infected birds. In this study, we describe the isolation of chicken primary bursal cells from the bursa of Fabricius, the culture and infection of the cells with IBDV, and the quantification of viral replication. The addition of chicken CD40 ligand significantly increased cell proliferation fourfold over six days of culture and significantly enhanced cell viability. Two strains of IBDV, a cell-culture adapted strain, D78, and a very virulent strain, UK661, replicated well in the ex vivo cell cultures. This model will be of use in determining how cells respond to IBDV infection and will permit a reduction in the number of infected birds used in IBDV pathogenesis studies. The model can also be expanded to include other viruses and could be applied to different species of birds.

Bursal transcriptome profiling of different inbred chicken lines reveals key differentially expressed genes at 3 days post-infection with very virulent infectious bursal disease virus

Infectious bursal disease is a highly contagious disease in the poultry industry and causes immunosuppression in chickens. Genome-wide regulations of immune response genes of inbred chickens with different genetic backgrounds, following very virulent infectious bursal disease virus (vvIBDV) infection are poorly characterized. Therefore, this study aims to analyse the bursal tissue transcriptome of six inbred chicken lines 6, 7, 15, N, O and P following infection with vvIBDV strain UK661 using strand-specific next-generation sequencing, by highlighting important genes and pathways involved in the infected chicken during peak infection at 3 days post-infection. All infected chickens succumbed to the infection without major variations among the different lines. However, based on the viral loads and bursal lesion scoring, lines P and 6 can be considered as the most susceptible lines, while lines 15 and N were regarded as the least affected lines. Transcriptome profiling of the bursa identified 4588 genes to be differentially expressed, with 2985 upregulated and 1642 downregulated genes, in which these genes were commonly or uniquely detected in all or several infected lines. Genes that were upregulated are primarily pro-inflammatory cytokines, chemokines and IFN-related. Various genes that are associated with B-cell functions and genes related to apoptosis were downregulated, together with the genes involved in p53 signalling. In conclusion, bursal transcriptome profiles of different inbred lines showed differential expressions of pro-inflammatory cytokines and chemokines, Th1 cytokines, JAK-STAT signalling genes, MAPK signalling genes, and their related pathways following vvIBDV infection.

Differential gene expression in chicken primary B cells infected ex vivo with attenuated and very virulent strains of infectious bursal disease virus (IBDV)

Infectious bursal disease virus (IBDV) belongs to the family Birnaviridae and is economically important to the poultry industry worldwide. IBDV infects B cells in the bursa of Fabricius (BF), causing immunosuppression and morbidity in young chickens. In addition to strains that cause classical Gumboro disease, the so-called 'very virulent' (vv) strain, also in circulation, causes more severe disease and increased mortality. IBDV has traditionally been controlled through the use of live attenuated vaccines, with attenuation resulting from serial passage in non-lymphoid cells. However, the factors that contribute to the vv or attenuated phenotypes are poorly understood. In order to address this, we aimed to investigate host cell-IBDV interactions using a recently described chicken primary B-cell model, where chicken B cells are harvested from the BF and cultured ex vivo in the presence of chicken CD40L. We demonstrated that these cells could support the replication of IBDV when infected ex vivo in the laboratory. Furthermore, we evaluated the gene expression profiles of B cells infected with an attenuated strain (D78) and a very virulent strain (UK661) by microarray. We found that key genes involved in B-cell activation and signalling (TNFSF13B, CD72 and GRAP) were down-regulated following infection relative to mock, which we speculate could contribute to IBDV-mediated immunosuppression. Moreover, cells responded to infection by expressing antiviral type I IFNs and IFN-stimulated genes, but the induction was far less pronounced upon infection with UK661, which we speculate could contribute to its virulence.