Transcriptomic analyses of host-virus interactions during in vitro infection with wild-type and glycoprotein g-deficient (ΔgG) strains of ILTV in primary and continuous cell cultures

Infectious laryngotracheitis (ILT) remains a significant concern for the poultry industry worldwide due to its impact on animal welfare and its substantial economic consequences. The disease is caused by the alphaherpesvirus, infectious laryngotracheitis virus (ILTV). This study investigated in vitro host-virus interactions of a glycoprotein G (gG) deletion mutant vaccine strain of ILTV (ΔgG ILTV), and its parent wild-type strain (CSW-1 ILTV). Inoculations were performed separately for the two strains of ILTV using both a primary (chicken embryonic kidney, CEK) and a continuous culture (leghorn male hepatoma, LMH) of chicken cells. Transcriptome analysis was performed at 12 hours post infection. Each cell-type displayed distinct effects on host and viral gene transcription, with a greater number of viral and host genes differentially transcribed in CEK cells and LMH cells, respectively. Both cell-types infected with either strain demonstrated enrichment of pathways related to signalling, and gene ontologies (GO) associated with chemotaxis. Infection with either strain upregulated both SOCS proteins and certain proto-oncogenes, which may contribute to prolonged viral persistence by promoting immunosuppression and preventing apoptosis, respectively. Patterns of gene transcription related to cytokines, chemokines, endosomal TLRs, and interferon responses, as well as pathways associated with histone acetylation, transport, and extracellular matrix organization were similar within each cell type, regardless of the viral strain. In CEK cells, GO terms and pathways were downregulated uniquely after CSW-1 ILTV infection, indicating a viral-strain specific effect in this cell-type. Overall, this study highlights that the observed differences in host and ILTV gene transcription in vitro were more strongly influenced by the cell-types used rather than the presence or absence of gG. This underscores the importance of cell-line selection in studying host-virus interactions and interpreting experimental results.

Transcriptome sequencing reveals non-coding RNAs respond to porcine reproductive and respiratory syndrome virus and Haemophilus parasuis co-infection in Kele piglets

Co-infection with porcine reproductive and respiratory syndrome virus (PRRSV) and Haemophilus parasuis (HPS) has severely restricted the healthy development of pig breeding. Exploring disease resistance of non-coding RNAs in pigs co-infected with PRRSV and HPS is therefore critical to complement and elucidate the molecular mechanisms of disease resistance in Kele piglets and to innovate the use of local pig germplasm resources in China. RNA-seq of lungs from Kele piglets with single-infection of PRRSV or HPS and co-infection of both pathogens was performed. Two hundred and twenty-five differentially expressed long non-coding RNAs (DElncRNAs) and 30 DEmicroRNAs (DEmiRNAs) were identified and characterized in the PRRSV and HPS co-infection (PRRSV-HPS) group. Compared with the single-infection groups, 146 unique DElncRNAs, 17 unique DEmiRNAs, and 206 target differentially expressed genes (DEGs) were identified in the PRRSV-HPS group. The expression patterns of 20 DEmiRNAs and DElncRNAs confirmed by real-time quantitative polymerase chain reaction (RT-qPCR) were consistent with those determined by high-throughput sequencing. In the PRRSV-HPS group, the target DEGs were enriched in eight immune Gene Ontology terms relating to two unique DEmiRNAs and 16 DElncRNAs, and the unique target DEGs participated the host immune response to pathogens infection by affecting 15 immune-related Kyoto Encyclopedia of Genes and Genomes enrichment pathways. Notably, competitive endogenous RNA (ceRNA) networks of different groups were constructed, and the ssc-miR-671-5p miRNA was validated as a potential regulatory factor to regulate DTX4 and AEBP1 genes to achieve innate antiviral effects and inhibit pulmonary fibrosis by dual-luciferase reporter assays. These results provided insight into further study on the molecular mechanisms of resistance to PRRSV and HPS co-infection in Kele piglets.

Host Responses Following Infection with Canadian-Origin Wildtype and Vaccine Revertant Infectious Laryngotracheitis Virus

Infectious laryngotracheitis (ILT) is caused by Gallid herpesvirus-1 (GaHV-1) or infectious laryngotracheitis virus (ILTV) and was first described in Canadian poultry flocks. In Canada, ILTV infection is endemic in backyard flocks, and commercial poultry encounters ILT outbreaks sporadically. A common practice to control ILT is the use of live attenuated vaccines. However, outbreaks still occur in poultry flocks globally due to ILTV vaccine strains reverting to virulence and emergence of new ILTV strains due to recombination in addition to circulating wildtype strains. Recent studies reported that most of the ILT outbreaks in Canada were induced by the chicken-embryo-origin (CEO) live attenuated vaccine revertant strains with the involvement of a small percentage of wildtype ILTV. It is not known if the host responses induced by these two ILTV strains are different. The objective of the study was to compare the host responses elicited by CEO revertant and wildtype ILTV strains in chickens. We infected 3-week-old specific pathogen-free chickens with the two types of ILTV isolates and subsequently evaluated the severity of clinical and pathological manifestations, in addition to host responses. We observed that both of the isolates show high pathogenicity by inducing several clinical and pathological manifestations. A significant recruitment of immune cells at both 3 and 7 days post-infection (dpi) was observed in the tracheal mucosa and the lung tissues of the infected chickens with wildtype and CEO vaccine revertant ILTV isolates when compared to uninfected controls. Overall, this study provides a better understanding of the mechanism of host responses against ILTV infection.

Network Meta-Analysis of Chicken Microarray Data Following Avian Influenza Challenge—A Comparison of Highly and Lowly Pathogenic Strains

The current bioinformatics study was undertaken to analyze the transcriptome of chicken (Gallus gallus) after influenza A virus challenge. A meta-analysis was carried out to explore the host expression response after challenge with lowly pathogenic avian influenza (LPAI) (H1N1, H2N3, H5N2, H5N3 and H9N2) and with highly pathogenic avian influenza (HPAI) H5N1 strains. To do so, ten microarray datasets obtained from the Gene Expression Omnibus (GEO) database were normalized and meta-analyzed for the LPAI and HPAI host response individually. Different undirected networks were constructed and their metrics determined e.g., degree centrality, closeness centrality, harmonic centrality, subgraph centrality and eigenvector centrality. The results showed that, based on criteria of centrality, the CMTR1, EPSTI1, RNF213, HERC4L, IFIT5 and LY96 genes were the most significant during HPAI challenge, with PARD6G, HMG20A, PEX14, RNF151 and TLK1L having the lowest values. However, for LPAI challenge, ZDHHC9, IMMP2L, COX7C, RBM18, DCTN3, and NDUFB1 genes had the largest values for aforementioned criteria, with GTF3C5, DROSHA, ATRX, RFWD2, MED23 and SEC23B genes having the lowest values. The results of this study can be used as a basis for future development of treatments/preventions of the effects of avian influenza in chicken.

Cytokine gene transcription in the trachea, Harderian gland, and trigeminal ganglia of chickens inoculated with virulent infectious laryngotracheitis virus (ILTV) strain

The objective of this study was to determine how cytokine transcription profiles correlate with patterns of infectious laryngotracheitis virus (ILTV) replication in the trachea, Harderian gland, and trigeminal ganglia during the early and late stages of infection after intratracheal inoculation. Viral genomes and transcripts were detected in the trachea and Harderian gland but not in trigeminal ganglia. The onset of viral replication in the trachea was detected at day one post-infection and peaked by day three post-infection. The peak of pro-inflammatory (CXCLi2, IL-1β, IFN-γ) and anti-inflammatory (IL-13, IL-10) cytokine gene transcription, 5 days post-infection, coincided with the increased recruitment of inflammatory cells, extensive tissue damage, and limiting of virus replication in the trachea. In contrast, transcription of the IFN-β gene in the trachea remained unaffected suggesting that ILTV infection blocks type I interferon responses. In the Harderian gland, the most evident transcription change was the early and transient upregulation of the IFN-γ gene at 1 day post-infection, which suggests that the Harderian gland is prepared to rapidly respond to ILTV infection. Overall, results from this study suggest that regulation of Th1 effector cells and macrophage activity by Th1/2 cytokines was pertinent to maintain a balanced immune response capable of providing an adequate Th1-mediated protective immunity, while sustaining some immune homeostasis in preparation for the regeneration of the tracheal mucosa.

Genome-wide detection of CNVs associated with beak deformity in chickens using high-density 600K SNP arrays

Beak deformity (crossed beaks) is found in several indigenous chicken breeds including Beijing-You studied here. Birds with deformed beaks have reduced feed intake and poor production performance. Recently, copy number variation (CNV) has been examined in many species and is recognized as a source of genetic variation, especially for disease phenotypes. In this study, to unravel the genetic mechanisms underlying beak deformity, we performed genome-wide CNV detection using Affymetrix chicken high-density 600K data on 48 deformed-beak and 48 normal birds using penncnv. As a result, two and eight CNV regions (CNVRs) covering 0.32 and 2.45 Mb respectively on autosomes were identified in deformed-beak and normal birds respectively. Further RT-qPCR studies validated nine of the 10 CNVRs. The ratios of six CNVRs were significantly different between deformed-beak and normal birds (P < 0.01). Within these six regions, three and 21 known genes were identified in deformed-beak and normal birds respectively. Bioinformatics analysis showed that these genes were enriched in six GO terms and one KEGG pathway. Five candidate genes in the CNVRs were further validated using RT-qPCR. The expression of LRIG2 (leucine rich repeats and immunoglobulin like domains 2) was lower in birds with deformed beaks (P < 0.01). Therefore, the LRIG2 gene could be considered a key factor in view of its known functions and its potential roles in beak deformity. Overall, our results will be helpful for future investigations of the genomic structural variations underlying beak deformity in chickens.

Accurate characterization of the IFITM locus using MiSeq and PacBio sequencing shows genetic variation in Galliformes

Background

Interferon inducible transmembrane (IFITM) proteins are effectors of the immune system widely characterized for their role in restricting infection by diverse enveloped and non-enveloped viruses. The chicken IFITM (chIFITM) genes are clustered on chromosome 5 and to date four genes have been annotated, namely chIFITM1, chIFITM3, chIFITM5 and chIFITM10. However, due to poor assembly of this locus in the Gallus Gallus v4 genome, accurate characterization has so far proven problematic. Recently, a new chicken reference genome assembly Gallus Gallus v5 was generated using Sanger, 454, Illumina and PacBio sequencing technologies identifying considerable differences in the chIFITM locus over the previous genome releases.

Methods

We re-sequenced the locus using both Illumina MiSeq and PacBio RS II sequencing technologies and we mapped RNA-seq data from the European Nucleotide Archive (ENA) to this finalized chIFITM locus. Using SureSelect probes capture probes designed to the finalized chIFITM locus, we sequenced the locus of a different chicken breed, namely a White Leghorn, and a turkey.

Results

We confirmed the Gallus Gallus v5 consensus except for two insertions of 5 and 1 base pair within the chIFITM3 and B4GALNT4 genes, respectively, and a single base pair deletion within the B4GALNT4 gene. The pull down revealed a single amino acid substitution of A63V in the CIL domain of IFITM2 compared to Red Jungle fowl and 13, 13 and 11 differences between IFITM1, 2 and 3 of chickens and turkeys, respectively. RNA-seq shows chIFITM2 and chIFITM3 expression in numerous tissue types of different chicken breeds and avian cell lines, while the expression of the putative chIFITM1 is limited to the testis, caecum and ileum tissues.

Conclusions

Locus resequencing using these capture probes and RNA-seq based expression analysis will allow the further characterization of genetic diversity within Galliformes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3801-8) contains supplementary material, which is available to authorized users.

Expression of perforin, granzyme A and Fas ligand mRNA in caecal tissues upon Eimeria tenella infection of naïve and immune chickens

Cytotoxic cells of the immune system may kill infected or transformed host cells via the perforin/granzyme or the Fas ligand (FasL) pathways. The purpose of this study was to determine mRNA expression of perforin, granzyme A and FasL in Eimeria tenella-infected tissues at primary infection and infection of immune chickens as an indirect measure of cytotoxic cell activity. Chickens were rendered immune by repeated E. tenella infections, which were manifested as an absence of clinical signs or pathological lesions and significantly reduced oocyst production upon challenge infection. During primary E. tenella infection, perforin, granzyme A and FasL mRNA expression in caecal tissue was significantly increased at 10 days after infection, compared to uninfected birds. In contrast, at infection of immune birds, perforin and granzyme A mRNA expression in caecal tissue was significantly increased during the early stages of E. tenella challenge infection, days 1-4, which coincided with a substantial reduction of parasite replication in these birds. These results indicate the activation of cytotoxic pathways in immune birds and support a role for cytotoxic T cells in the protection against Eimeria infections.

Infection of Goose with Genotype VIId Newcastle Disease Virus of Goose Origin Elicits Strong Immune Responses at Early Stage

Newcastle disease (ND), caused by virulent strains of Newcastle disease virus (NDV), is a highly contagious disease of birds that is responsible for heavy economic losses for the poultry industry worldwide. However, little is known about host-virus interactions in waterfowl, goose. In this study, we aim to characterize the host immune response in goose, based on the previous reports on the host response to NDV in chickens. Here, we evaluated viral replication and mRNA expression of 27 immune-related genes in 10 tissues of geese challenged with a genotype VIId NDV strain of goose origin (go/CH/LHLJ/1/06). The virus showed early replication, especially in digestive and immune tissues. The expression profiles showed up-regulation of Toll-like receptor (TLR)1–3, 5, 7, and 15, avian β-defensin (AvBD) 5–7, 10, 12, and 16, cytokines [interleukin (IL)-8, IL-18, IL-1β, and interferon-γ], inducible NO synthase (iNOS), and MHC class I in some tissues of geese in response to NDV. In contrast, NDV infection suppressed expression of AvBD1 in cecal tonsil of geese. Moreover, we observed a highly positive correlation between viral replication and host mRNA expressions of TLR1-5 and 7, AvBD4-6, 10, and 12, all the cytokines measured, MHC class I, FAS ligand, and iNOS, mainly at 72 h post-infection. Taken together, these results demonstrated that NDV infection induces strong innate immune responses and intense inflammatory responses at early stage in goose which may associate with the viral pathogenesis.

Spleen Transcriptome Profile of Muscovy Ducklings in Response to Infection With Muscovy Duck Reovirus

Muscovy duck reovirus (MDRV) causes high morbidity and mortality in ducklings. However, the molecular basis for pathogenesis of this virus remains poorly understood, and the complete genome sequence of Muscovy duck is lacking. Here we report the transcriptome profile of Muscovy ducks in response to MDRV infection. RNA sequencing technology was employed to obtain a representative complement of transcripts from the spleen of ducklings, and then differential gene expression was analyzed between MDRV-YB strain infected ducks and noninfected ducks. This analysis generated 65,536 unigenes. Of these, 6458 genes were found to be significantly differentially expressed between the infected and noninfected groups. The symptom and pathology of ducks infected with MDRV-YB was correlated with the greater proportion of decreased expression genes (4906) than increased expression (1552) level. Gene ontology analysis assigned differentially expressed genes to the categories: "biological processes," "cellular functions," and "molecular functions." Differentially expressed genes involved in the innate immune system were analyzed further, and 128 of these genes showed decreased expression and 86 showed increased expression between the infected and noninfected groups. These genes represented the Janus kinase-signal transducer and activator of transcription signaling pathway, and the retinoic acid-inducible gene I (RIG-I)-like and Toll-like receptor (TLR) signaling pathways and included interferon (IFN) α, IFNγ, interleukin 6, RIG-I, and TLR4. The data were verified by SYBR fluorescence quantitative polymerase chain reaction (SYBR-qPCR). Our findings offer new insight into the host immune response to MDRV infection.

Transcriptomic Characterization of Innate and Acquired Immune Responses in Red-Legged Partridges (Alectoris rufa): A Resource for Immunoecology and Robustness Selection

Present and future challenges for wild partridge populations include the resistance against possible disease transmission after restocking with captive-reared individuals, and the need to cope with the stress prompted by new dynamic and challenging scenarios. Selection of individuals with the best immune ability may be a good strategy to improve general immunity, and hence adaptation to stress. In this study, non-infectious challenges with phytohemagglutinin (PHA) and sheep red blood cells allowed the classification of red-legged partridges (Alectoris rufa) according to their overall immune responses (IR). Skin from the area of injection of PHA and spleen, both from animals showing extreme high and low IR, were selected to investigate the transcriptional profiles underlying the different ability to cope with pathogens and external aggressions. RNA-seq yielded 97 million raw reads from eight sequencing libraries and approximately 84% of the processed reads were mapped to the reference chicken genome. Differential expression analysis identified 1488 up- and 107 down-regulated loci in individuals with high IR versus low IR. Partridges displaying higher innate IR show an enhanced activation of host defence gene pathways complemented with a tightly controlled desensitization that facilitates the return to cellular homeostasis. These findings indicate that the immune system ability to respond to aggressions (either diseases or stress produced by environmental changes) involves extensive transcriptional and post-transcriptional regulations, and expand our understanding on the molecular mechanisms of the avian immune system, opening the possibility of improving disease resistance or robustness using genome assisted selection (GAS) approaches for increased IR in partridges by using genes such as AVN or BF2 as markers. This study provides the first transcriptome sequencing data of the Alectoris genus, a resource for molecular ecology that enables integration of genomic tools in further studies.

Methylome Analysis in Chickens Immunized with Infectious Laryngotracheitis Vaccine

In this study we investigated the methylome of chickens immunized with Infectious laryngotracheitis (ILT) vaccine derived from chicken embryos. Methyl-CpG binding domain protein-enriched genome sequencing (MBD-Seq) method was employed in the detection of the 1,155 differentially methylated regions (DMRs) across the entire genome. After validation, we ascertained the genomic DMRs distribution and annotated them regarding genes, transcription start sites (TSS) and CpG islands. We found that global DNA methylation decreased in vaccinated birds, presenting 704 hypomethylated and 451 hypermethylated DMRs, respectively. Additionally, we performed an enrichment analysis detecting gene networks, in which cancer and RNA post-transcriptional modification appeared in the first place, followed by humoral immune response, immunological disease and inflammatory disease. The top four identified canonical pathways were EIF2 signaling, regulation of EIF4 and p70S6K signaling, axonal guidance signaling and mTOR signaling, providing new insight regarding the mechanisms of ILT etiology. Lastly, the association between DNA methylation and differentially expressed genes was examined, and detected negative correlation in seventeen of the eighteen genes.