Transcriptional profiling of host gene expression in chicken liver tissues infected with oncogenic Marek's disease virus

The Efficacy of a Live Attenuated TW I-Type Infectious Bronchitis Virus Vaccine Candidate

Infectious bronchitis (IB) is a highly contagious avian disease caused by infection with infectious bronchitis virus (IBV), which seriously affects the development of the global poultry industry. The distribution of TW I-type IBV in China has increased in recent years, becoming a widespread genotype. We previously isolated a TW I-type IBV strain termed CK/CH/GD/GZ14 in 2014, but its pathogenicity and possibility for vaccine development were not explored. Therefore, this research aimed to develop a live-attenuated virus vaccine based on the CK/CH/GD/GZ14 strain. The wild type IBV CK/CH/GD/GZ14 strain was serially passaged in SPF embryos for 145 generations. The morbidity and mortality rate of wild-type strain in 14 day-old chickens is 100% and 80% respectively, while the morbidity rate in the attenuated strain was 20% in the 95th and 105th generations and there was no death. Histopathological observations showed that the pathogenicity of the 95th and 105th generations in chickens was significantly weakened. Further challenge experiments confirmed that the attenuated CK/CH/GD/GZ14 strain in the 95th and 105th generations could resist CK/CH/GD/GZ14 (5th generation) infection and the protection rate was 80%. Tracheal cilia stagnation, virus shedding, and viral load experiments confirmed that the 95th and 105th generations provide good immune protection in chickens, and the immunogenicity of the 105th generation is better than that of the 95th generation. These data suggest that the attenuated CK/CH/GD/GZ14 strain in the 105th generation may be applied as a vaccine candidate against TW I-type IBV.

Atmospheric H2S triggers immune damage by activating the TLR-7/MyD88/NF-κB pathway and NLRP3 inflammasome in broiler thymus

Atmospheric hydrogen sulfide (H₂S) is a highly toxic air pollutant that has a negative effect on human health and animal welfare. The immunotoxicity of H₂S has been explored previously, but its mechanism still needs to be clarified, especially in chickens. To further evaluate the immunotoxicity of H₂S, 1-day-old broilers were recruited and exposed to atmospheric H₂S for 42 days of age. Our results showed that H₂S significantly reduced the thymus index and the CD4⁺ and CD8⁺ T-lymphocyte numbers and that it also changed the CD4⁺/CD8⁺ ratio. The morphological analysis showed that H₂S incrassated the medulla and generated inflammatory infiltration. In addition, it caused the mitochondria to swell and the chromatin to condense, and destroyed nuclear structures were observed. We also conducted bioinformation and transcriptomic analyses to delve the mechanism of H₂S toxicity in chicken thymus. We measured 172 differently expression genes (DEGs) after H₂S exposure and further filtrated the DEGs that are related to inflammation and cell death that play a critical role in immune function. We concluded that H₂S significantly increased IL-1β, IL-4 and IL-10 levels, whereas it downregulated IL-12 and IFN-γ. This study confirmed that H₂S triggered the thymus inflammatory response and caused a Th1/Th2 imbalance. Moreover, our results demonstrated that H₂S triggered the TLR-7/MyD88/NF-κB pathway to promote NLRP3 inflammasome activation. In conclusion, atmospheric H₂S actives the TLR-7/MyD88/NF-κB pathway and the NLRP3 inflammasome to promote an inflammatory response, which then causes tissues damage in broiler thymus. These results provide new insights for unveiling the immunotoxic effects of H₂S.

Transcriptome analysis of responses to bluetongue virus infection in Aedes albopictus cells

Background

Bluetongue virus (BTV) causes a disease among wild and domesticated ruminants which is not contagious, but which is transmitted by biting midges of the Culicoides species. BTV can induce an intense cytopathic effect (CPE) in mammalian cells after infection, although Culicoides- or mosquito-derived cell cultures cause non-lytic infection with BTV without CPE. However, little is known about the transcriptome changes in Aedes albopictus cells infected with BTV.

Methods

Transcriptome sequencing was used to identify the expression pattern of mRNA transcripts in A. albopictus cells infected with BTV, given the absence of the Culicoides genome sequence. Bioinformatics analyses were performed to examine the biological functions of the differentially expressed genes. Subsequently, quantitative reverse transcription–polymerase chain reaction was utilized to validate the sequencing data.

Results

In total, 51,850,205 raw reads were generated from the BTV infection group and 51,852,293 from the control group. A total of 5769 unigenes were common to both groups; only 779 unigenes existed exclusively in the infection group and 607 in the control group. In total, 380 differentially expressed genes were identified, 362 of which were up-regulated and 18 of which were down-regulated. Bioinformatics analyses revealed that the differentially expressed genes mainly participated in endocytosis, FoxO, MAPK, dorso-ventral axis formation, insulin resistance, Hippo, and JAK-STAT signaling pathways.

Conclusion

This study represents the first attempt to investigate transcriptome-wide dysregulation in A. albopictus cells infected with BTV. The understanding of BTV pathogenesis and virus–vector interaction will be improved by global transcriptome profiling.

Electronic supplementary material

The online version of this article (10.1186/s12866-019-1498-3) contains supplementary material, which is available to authorized users.

The molecular basis for host responses to Marek's disease viruses integrated with different retro-viral long terminal repeat

Integration of retro-viral long terminal repeat (LTR) into the Marek's disease virus (MDV) genome can occur both in co-cultivation cell cultures and naturally in dual infected chickens. It is clear that the LTR insert is associated with the pathogenicity of MDV. The objective of this study was to compare the host responses to MDV with a different retro-viral LTR insert. Gene-chip containing chicken genome was employed to investigate the gene transcription profile of chicken embryo fibroblasts cells, and 795 genes were differentially expressed in chicken embryo fibroblasts infected with GX0101 with a reticuloendotheliosis virus LTR insert as compared to GX0101-ALV-LTR significantly. The differentially expressed genes were mostly associated with the regulation of transcription and the development of multiple organs. Based on the bio functions of the differential genes, infection of GX0101 was predicated with a greater development disorder of multiple systems, resulting in higher growth retardation, mortality, tumorigenicity, and immunosuppression in chickens than GX0101-ALV-LTR. Collectively, our results provided valuable insights into elucidation of the possible relationship between retro-viral LTR insert and the observed phenotypes caused by MDV recombinant viruses.

Avian Influenza Virus Subtype H9N2 Affects Intestinal Microbiota, Barrier Structure Injury, and Inflammatory Intestinal Disease in the Chicken Ileum

Avian influenza virus subtype H9N2 (H9N2 AIV) has caused significant losses to the poultry industry due to the high mortality associated with secondary infections attributable to E. coli. This study tries to address the underlying secondary mechanisms after H9N2 AIV infection. Initially, nine day-old specific pathogen-free chickens were assigned to control (uninfected) and H9N2-infected groups, respectively. Using Illumina sequencing, histological examination, and quantitative real-time PCR, it was found that H9N2 AIV caused intestinal microbiota disorder, injury, and inflammatory damage to the intestinal mucosa. Notably, the genera Escherichia, especially E. coli, significantly increased (p < 0.01) at five days post-infection (dpi), while Lactobacillus, Enterococcus, and other probiotic organisms were significantly reduced (p < 0.01). Simultaneously, the mRNA expression of tight junction proteins (ZO-1, claudin 3, and occludin), TFF2, and Muc2 were significantly reduced (p < 0.01), indicating the destruction of the intestinal epithelial cell tight junctions and the damage of mucin layer construction. Moreover, the mRNA expression of proinflammatory cytokines IFN-γ, IL-22, IFN-α, and IL-17A in intestinal epithelial cells were significantly upregulated, resulting in the inflammatory response and intestinal injury. Our findings may provide a theoretical basis for observed gastroenteritis-like symptoms such as diarrhea and secondary E. coli infection following H9N2 AIV infection.

Common microRNA⁻mRNA Interactions in Different Newcastle Disease Virus-Infected Chicken Embryonic Visceral Tissues

To investigate the roles and explore the altered expression of microRNAs (miRNAs) and mRNAs in chicken embryos in response to Newcastle disease virus (NDV) infection, deep sequencing was performed. Then, a conjoint analysis of small RNA-seq and mRNA-seq was performed to screen interactional miRNA⁻mRNA pairs during NDV infection. In total, 15 and 17 up- and downregulated miRNAs were identified that potentially targeted 4279 and 6080 mRNAs in NDV-infected chicken embryonic tissues, respectively; in addition, 595 upregulated and 480 downregulated mRNAs were identified. The conjoint analysis of the obtained data identified 1069 miRNA⁻mRNA pairs. Among these pairs, 130 pairs were related to immune or inflammatory responses. The relationship between gga-miR-203a and its target transglutaminase 2 (TGM2) was confirmed using a dual-luciferase reporter system and a real time quantitative polymerase chain reaction (RT-qPCR) assay. Overall, the discovery of miRNAs, mRNAs, and their potential pairing relationships, which may be involved in the regulation of NDV infection, will facilitate our understanding of the complex regulatory relationship between the host and the virus.

Dynamic Changes in the Splenic Transcriptome of Chickens during the Early Infection and Progress of Marek's Disease

Gallid alphaherpesvirus 2 (GaHV2) is an oncogenic avian herpesvirus inducing Marek’s disease (MD) and rapid-onset T-cell lymphomas. To reveal molecular events in MD pathogenesis and tumorigenesis, the dynamic splenic transcriptome of GaHV2-infected chickens during early infection and pathogenic phases has been determined utilizing RNA-seq. Based on the significant differentially expressed genes (DEGs), analysis of gene ontology, KEGG pathway and protein-protein interaction network has demonstrated that the molecular events happening during GaHV2 infection are highly relevant to the disease course. In the ‘Cornell Model’ description of MD, innate immune responses and inflammatory responses were established at early cytolytic phase but persisted until lymphoma formation. Humoral immunity in contrast began to play a role firstly in the intestinal system and started at late cytolytic phase. Neurological damage caused by GaHV2 is first seen in early cytolytic phase and is then sustained throughout the following phases over a long time period. During the proliferative phase many pathways associated with transcription and/or translation were significantly enriched, reflecting the cell transformation and lymphoma formation. Our work provides an overall view of host responses to GaHV2 infection and offers a meaningful basis for further studies of MD biology.

Comparative transcriptomics of genetically divergent lines of chickens in response to Marek's disease virus challenge at cytolytic phase

Marek's disease (MD), caused by Marek's disease virus (MDV), remains an economically significant threat to the poultry industry worldwide. Genetic resistance to MD is a promising alternative strategy to augment current control measures (vaccination and management). However, only a few functional genes reportedly conferring MD resistance have been identified. Here, we performed a comparative transcriptomics analysis of two highly inbred yet genetically divergent lines of chickens (line 63 and 72) that are resistant and susceptible to MD, respectively, in response to a very virulent plus strain of MDV (vv+MDV) challenge at cytolytic phase. A total of 203 DEGs in response to MDV challenge were identified in the two lines. Of these, 96 DEGs were in common for both lines, in addition to 36 and 71 DEGs that were specific for line 63 and 72, respectively. Functional enrichment analysis results showed the DEGs were significantly enriched in GO terms and pathways associated with immune response. Especially, the four DEGs, FGA, ALB, FN1, and F13A1 that reportedly facilitate virus invasion or immunosuppression, were found to be significantly up-regulated in the susceptible line 72 but down-regulated in the resistant line 63 birds. These results provide new resources for future studies to further elucidate the genetic mechanism conferring MD resistance.

Transcriptional and Bioinformatic Analysis Provide a Relationship between Host Response Changes to Marek's Disease Viruses Infection and an Integrated Long Terminal Repeat

GX0101, Marek's disease virus (MDV) strain with a long terminal repeat (LTR) insert of reticuloendotheliosis virus (REV), was isolated from CVI988/Rispens vaccinated birds showing tumors. We have constructed a LTR deleted strain GX0101ΔLTR in our previous study. To compare the host responses to GX0101 and GX0101ΔLTR, chicken embryo fibroblasts (CEF) cells were infected with two MDV strains and a gene-chip containing chicken genome was employed to examine gene transcription changes in host cells in the present study. Of the 42,368 chicken transcripts on the chip, there were 2199 genes that differentially expressed in CEF infected with GX0101 compared to GX0101ΔLTR significantly. Differentially expressed genes were distributed to 25 possible gene networks according to their intermolecular connections and were annotated to 56 pathways. The insertion of REV LTR showed the greatest influence on cancer formation and metastasis, followed with immune changes, atherosclerosis, and nervous system disorders in MDV-infected CEF cells. Based on these bio functions, GX0101 infection was predicated with a greater growth and survival inhibition but lower oncogenicity in chickens than GX0101ΔLTR, at least in the acute phase of infection. In summary, the insertion of REV LTR altered the expression of host genes in response to MDV infection, possibly resulting in novel phenotypic properties in chickens. Our study has provided the evidence of retroviral insertional changes of host responses to herpesvirus infection for the first time, which will promote to elucidation of the possible relationship between the LTR insertion and the observed phenotypes.

Transcriptional analysis of host responses to Marek's disease virus infection in chicken thymus

Marek's disease virus (MDV) is a cell-associated alpha-herpesvirus that causes T-cell lymphomas and nervous disorders in chickens. Different from other lymphoid organs, the thymus is the site of T-cell maturation and differentiation. However, the transcriptional response to MDV infection in the chicken thymus is still not known. In this study, we performed genome-wide expression analysis in thymus tissues of RB1B-infected chickens at different time points to investigate the molecular mechanisms of MDV pathogenesis. The number of differentially expressed genes with 2-fold or higher changes (>2) are as follows: 1,250 genes (7 dpi), 834 genes (14 dpi), 1,958 genes (21 dpi), and 2,306 genes (28 dpi). Gene ontology enrichment analysis revealed that the upregulated genes were involved in immune and inflammatory response at 7 dpi; angiogenesis, cytoskeleton organization, cell adhesion, and signal transduction showed different expressions at 21 and 28 dpi. The expression pattern of 18 randomly selected genes was confirmed by real-time RT-PCR. Several differently expressed host genes associated with tumor development are discussed. We identified the global host-gene expression pattern in the thymus of chickens that responded to MDV infection. The present data may provide groundwork for future investigation in the biology and pathogenesis of MDV.

Characterization of cytotoxicity-related gene expression in response to virulent Marek's disease virus infection in the bursa of Fabricius

Cell-mediated cytotoxic responses are critical for control of Marek's disease virus (MDV) infection and tumour development. However, the mechanisms of virus clearance mediated by cytotoxic responses in the bursa of Fabricius of chickens during MDV infection are not fully understood. In this study, the host cytotoxic responses during MDV infection in the bursa were investigated by examining the expression of genes in the cell lysis pathways. Partial up-regulation existed in the expression of the important cytolytic molecule granzyme A (GzmA), Fas, NK lysin and DNA repair enzyme Ape1, whereas little or no expression appeared in other cytolytic molecules, including perforin (PFN) and Fas ligand (FasL), and molecules involved in DNA repair and apoptosis in the bursa during MDV infection. These results suggest that less sustained cytotoxic activities are generated in the bursa of MDV-infected chickens. The findings of this study provide a more detailed insight into the host cytotoxic responses to MDV infection.

Marek's disease virus infection induces widespread differential chromatin marks in inbred chicken lines

Background

Marek’s disease (MD) is a neoplastic disease in chickens caused by the MD virus (MDV). Successful vaccine development against MD has resulted in increased virulence of MDV and the understanding of genetic resistance to the disease is, therefore, crucial to long-term control strategies. Also, epigenetic factors are believed to be one of the major determinants of disease response.

Results

Here, we carried out comprehensive analyses of the epigenetic landscape induced by MDV, utilizing genome-wide histone H3 lysine 4 and lysine 27 trimethylation maps from chicken lines with varying resistance to MD. Differential chromatin marks were observed on genes previously implicated in the disease such as MX1 and CTLA-4 and also on genes reported in other cancers including IGF2BP1 and GAL. We detected bivalent domains on immune-related transcriptional regulators BCL6, CITED2 and EGR1, which underwent dynamic changes in both lines as a result of MDV infection. In addition, putative roles for GAL in the mechanism of MD progression were revealed.

Conclusion

Our results confirm the presence of widespread epigenetic differences induced by MD in chicken lines with different levels of genetic resistance. A majority of observed epigenetic changes were indicative of increased levels of viral infection in the susceptible line symptomatic of lowered immunocompetence in these birds caused by early cytolytic infection. The GAL system that has known anti-proliferative effects in other cancers is also revealed to be potentially involved in MD progression. Our study provides further insight into the mechanisms of MD progression while revealing a complex landscape of epigenetic regulatory mechanisms that varies depending on host factors.