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

The long noncoding RNA loc107053557 acts as a gga-miR-3530-5p sponge to suppress the replication of vvIBDV through regulating STAT1 expression

Infectious bursal disease virus (IBDV) causes immunosuppression and high mortality in young chickens. Long non-coding RNAs (lncRNAs) and microRNAs (miRNAs) are important regulators during viral infection. However, detailed the regulatory mechanisms of lncRNA-miRNA-mRNA have not yet been described in IBDV infection. Here, we analysed the role of lncRNA53557/gga-miR-3530-5p/STAT1 axis in very virulent IBDV (vvIBDV) infection. Evidently upregulated expression of lncRNA53557 was observed in bursa of Fabricius and DT40 cells. Meanwhile, overexpression of lncRNA53557 promoted STAT1 expression and inhibited vvIBDV replication and vice versa, indicating that the upregulation of lncRNA53557 was part of the host antiviral defence. The subcellular fractionation assay confirmed that lncRNA53557 can be localized in the cytoplasm. Further, dual-luciferase reporter, RNA pulldown, FISH and RT-qPCR assays revealed that lncRNA53557 were directly bound to gga-miR-3530-5p and had a negative regulatory relationship between them. Subsequent mechanistic analysis showed that lncRNA53557 acted as a competing endogenous RNA (ceRNA) of gga-miR-3530-5p to relieve the repressive effect of gga-miR-3530-5p on its target STAT1, as well as Mx1, OASL, and ISG15, thereby suppressing vvIBDV replication. The study reveals that a network of enriched lncRNAs and lncRNA-associated ceRNA is involved in the regulation of IBDV infection, offering new insight into the mechanisms underlying IBDV-host interaction.

CircMYO1B/miR-155 pathway is a common mechanism of stress-induced immunosuppression affecting immune response to three vaccines in chicken

Stress-induced immunosuppression (SIIS) can weaken the immune response effect of poultry vaccination, and bring huge hidden dangers and economic losses to the poultry industry. However, the detailed molecular mechanisms are still not fully understood. Unveiling the common mechanism of SIIS affecting the immune response to different vaccines is critical for detecting and minimizing the losses caused by SIIS. This study used glucocorticoid dexamethasone (Dex) to simulate SIIS, and three classic avian vaccines (including avian influenza virus (AIV), Newcastle disease virus (NDV), and infectious bursal disease virus (IBDV)) were used to induce immune responses in chicken. Quantitative real-time PCR (qRT-PCR) revealed the expression characteristics and functions of circMYO1B and miR-155 in the processes of SIIS affecting the immune response to the aforementioned avian vaccines, as well as their targeted regulatory relationship. Subsequent bioinformatics analysis predicted FOS, one of the potential target genes of miR-155. The results showed that circMYO1B/miR-155 pathway served as a key common mechanism by which SIIS affected the immune response to the three vaccines. Both heart and proventriculus appeared to be the crucial tissues for this process, with five days post immunization (dpi) emerging as the primary time of interest. Moreover, mitogen-activated protein kinase (MAPK) signaling system played a key role in modulating the immune response subsequent to SIIS administration. Our findings provide new insights into the immune function of competitive endogenous RNA (ceRNA), which have important function in the detection and treatment of SIIS affecting vaccine immunity.

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.

Functional Involvement of circRNAs in the Innate Immune Responses to Viral Infection

Effective viral clearance requires fine-tuned immune responses to minimize undesirable inflammatory responses. Circular RNAs (circRNAs) are a class of non-coding RNAs that are abundant and highly stable, formed by backsplicing pre-mRNAs, and expressed ubiquitously in eukaryotic cells, emerging as critical regulators of a plethora of signaling pathways. Recent progress in high-throughput sequencing has enabled a better understanding of the physiological and pathophysiological functions of circRNAs, overcoming the obstacle of the sequence overlap between circRNAs and their linear cognate mRNAs. Some viruses also encode circRNAs implicated in viral replication or disease progression. There is increasing evidence that viral infections dysregulate circRNA expression and that the altered expression of circRNAs is critical in regulating viral infection and replication. circRNAs were shown to regulate gene expression via microRNA and protein sponging or via encoding small polypeptides. Recent studies have also highlighted the potential role of circRNAs as promising diagnostic and prognostic biomarkers, RNA vaccines and antiviral therapy candidates due to their higher stability and lower immunogenicity. This review presents an up-to-date summary of the mechanistic involvement of circRNAs in innate immunity against viral infections, the current understanding of their regulatory roles, and the suggested applications.

Circular RNA circAKIRIN2 participates in the process of stress-induced immunosuppression affecting immune response to infectious bursal disease virus vaccine in chicken

At present, stress-induced immunosuppression is still a hidden threat that leads to immunization failure and outbreaks of poultry diseases, and causes huge economic losses to the modern poultry industry. However, the molecular mechanisms of stress-induced immunosuppression affecting viral vaccine immunity are still poorly understood. Here, we identified circAKIRIN2 as a conserved circular transcript in chicken, and explored its expression patterns in different immune states by quantitative real-time PCR (qRT-PCR), then conducted bioinformatics analysis. The results showed that circAKIRIN2 actively participated in the process of stress-induced immunosuppression affecting the immune response to infectious bursal disease virus (IBDV) vaccine. The key time points for circAKIRIN2 involving in the process were 2 day post immunization (dpi), 5 dpi, and 28 dpi, especially at the acquired immune stage. The important tissues that responded to the process included the heart, liver, and lung, all of which changed significantly. In addition, circAKIRIN2 as a competing endogenous RNA (ceRNA) sponging zinc finger and BTB domain containing 20 (ZBTB20) was a potential molecular mechanism for regulating immune functions in the process. In conclusion, circAKIRIN2 is a key regulatory factor for stress-induced immunosuppression affecting the IBDV vaccine immune response, and this study can provide a new perspective for exploring the molecular regulatory mechanisms of stress-induced immunosuppression affecting immune response.

Integrated Transcriptome Analysis Reveals mRNA-miRNA Pathway Crosstalk in Roman Laying Hens' Immune Organs Induced by AFB1

Aflatoxin B1 (AFB1) is a widely distributed contaminant in moldy corn, rice, soybean, and oil crops. Many studies have revealed its adverse effects, such as carcinogenicity, immunotoxicity, and hepatotoxicity, on the health of humans and animals. To investigate the immunotoxic effects on chicken immune organs induced by AFB1, we integrated RNA and small-RNA sequencing data of the spleen and the bursa of Fabricius to elucidate the response of the differentially expressed transcriptional profiles and related pathways. AFB1 consumption negatively influenced egg quality, but no obvious organ damage was observed compared to that of the control group. We identified 3918 upregulated and 2415 downregulated genes in the spleen and 231 upregulated and 65 downregulated genes in the bursa of Fabricius. We confirmed that several core genes related to immune and metabolic pathways were activated by AFB1. Furthermore, 42 and 19 differentially expressed miRNAs were found in the spleen and the bursa of Fabricius, respectively. Differentially expressed genes and target genes of differentially expressed miRNAs were mainly associated with cancer progression and immune response. The predicted mRNA-miRNA pathway network illustrated the potential regulatory mechanisms. The present study identified the transcriptional profiles and revealed potential mRNA-miRNA pathway crosstalk. This genetic regulatory network will facilitate the understanding of the immunotoxicity mechanisms of chicken immune organs induced by high concentrations of AFB1.

Host Combats IBDV Infection at Both Protein and RNA Levels

Infectious bursal disease (IBD) is an acute, highly contagious, and immunosuppressive avian disease caused by infectious bursal disease virus (IBDV). In recent years, with the emergence of IBDV variants and recombinant strains, IBDV still threatens the poultry industry worldwide. It seems that the battle between host and IBDV will never end. Thus, it is urgent to develop a more comprehensive and effective strategy for the control of this disease. A better understanding of the mechanisms underlying virus-host interactions would be of help in the development of novel vaccines. Recently, much progress has been made in the understanding of the host response against IBDV infection. If the battle between host and IBDV at the protein level is considered the front line, at the RNA level, it can be taken as a hidden line. The host combats IBDV infection at both the front and hidden lines. Therefore, this review focuses on our current understanding of the host response to IBDV infection at both the protein and RNA levels.

Genome-wide identification of chicken bursae of Fabricius miRNAs in response to very virulent infectious bursal disease virus

Infectious bursal disease virus (IBDV) can cause a highly contagious immunosuppressive disease in young chickens. MicroRNAs (miRNAs) are crucial regulators of gene expression and are involved in the pathogenesis of IBDV infection. To investigate the roles of miRNA in chicken bursae of Fabricius in response to very virulent IBDV (vvIBDV) infection, RNA sequencing was performed to compare the small RNA libraries from uninfected and vvIBDV-infected group which was infected for 3 days. A total of 77 differentially expressed (DE) miRNAs were identified in BF, of which 42 DE miRNAs were upregulated and 35 DE miRNAs were downregulated. A gene ontology analysis showed that genes associated with cellular processes, cells, and binding were enriched. Moreover, pathway analyses suggested that apoptosis, T cell receptor signaling pathways, and chemokine signaling pathways may be activated following vvIBDV infection. In addition, we predicted the target genes of DE miRNAs and constructed an miRNA-mRNA regulatory network. In total, 189 pairs of miRNA-target genes were identified, comprising 67 DE miRNAs and 73 mRNAs. In this network, gga-miR-1684b-3p was identified with the highest fold change, as well as gga-miR-1788-3p and gga-miR-3530-5p showed a high degree of change. The above three miRNAs were considered to play vital roles in vvIBDV-host interactions. This study was the first to perform a comprehensive analysis of DE miRNAs in the bursa of Fabricius in response to vvIBDV infection, and it provided new insights into molecular mechanisms underlying vvIBDV infection and pathogenesis.

Functional prediction of de novo uni-genes from chicken transcriptomic data following infectious bursal disease virus at 3-days post-infection

BACKGROUND

Infectious bursal disease (IBD) is an economically very important issue to the poultry industry and it is one of the major threats to the nation's food security. The pathogen, a highly pathogenic strain of a very virulent IBD virus causes high mortality and immunosuppression in chickens. The importance of understanding the underlying genes that could combat this disease is now of global interest in order to control future outbreaks. We had looked at identified novel genes that could elucidate the pathogenicity of the virus following infection and at possible disease resistance genes present in chickens.

RESULTS

A set of sequences retrieved from IBD virus-infected chickens that did not map to the chicken reference genome were de novo assembled, clustered and analysed. From six inbred chicken lines, we managed to assemble 10,828 uni-transcripts and screened 618 uni-transcripts which were the most significant sequences to known genes, as determined by BLASTX searches. Based on the differentially expressed genes (DEGs) analysis, 12 commonly upregulated and 18 downregulated uni-genes present in all six inbred lines were identified with false discovery rate of q-value < 0.05. Yet, only 9 upregulated and 13 downregulated uni-genes had BLAST hits against the Non-redundant and Swiss-Prot databases. The genome ontology enrichment keywords of these DEGs were associated with immune response, cell signalling and apoptosis. Consequently, the Weighted Gene Correlation Network Analysis R tool was used to predict the functional annotation of the remaining unknown uni-genes with no significant BLAST hits. Interestingly, the functions of the three upregulated uni-genes were predicted to be related to innate immune response, while the five downregulated uni-genes were predicted to be related to cell surface functions. These results further elucidated and supported the current molecular knowledge regarding the pathophysiology of chicken's bursal infected with IBDV.

CONCLUSION

Our data revealed the commonly up- and downregulated novel uni-genes identified to be immune- and extracellular binding-related, respectively. Besides, these novel findings are valuable contributions in improving the current existing integrative chicken transcriptomics annotation and may pave a path towards the control of viral particles especially towards the suppression of IBD and other infectious diseases in chickens.