Resistant and susceptible chicken lines show distinctive responses to Newcastle disease virus infection in the lung transcriptome

Breed-specific responses to coccidiosis in chickens: identification of intestinal bacteria linked to disease resistance

BACKGROUND

Coccidiosis, caused by Eimeria parasites, is a major enteric disease in poultry, significantly impacting animal health, production performance, and welfare. This disease imposes a substantial economic burden, costing the global poultry industry up to $13 billion annually. However, effective mitigation strategies for coccidiosis remain elusive. While different chicken breeds exhibit varying resistance to coccidiosis, no commensal bacteria have been directly linked to this resistance.

METHODS

To assess relative resistance of different breeds to coccidiosis, 10-day-old Fayoumi M5.1, Leghorn Ghs6, and Cobb chickens were challenged with 50,000 sporulated Eimeria maxima oocysts or mock-infected. Body weight changes, small intestinal lesions, and fecal oocyst shedding were evaluated on d 17. Ileal and cecal digesta were collected from individual animals on d 17 and subjected to microbiome analysis using 16S rRNA gene sequencing.  RESULTS: Fayoumi M5.1 chickens showed the lowest growth retardation, intestinal lesion score, fecal oocyst shedding, and pathobiont proliferation compared to Ghs6 and Cobb chickens. The intestinal microbiota of M5.1 chickens also differed markedly from the other two breeds under both healthy and coccidiosis conditions. Notably, group A Lactobacillus and Ligilactobacillus salivarius were the least prevalent in both the ileum and cecum of healthy M5.1 chickens, but became highly enriched and comparable to Ghs6 and Cobb chickens in response to coccidiosis. Conversely, Weissella, Staphylococcus gallinarum, and Enterococcus durans/hirae were more abundant in the ileum of healthy M5.1 chickens than in the other two breeds. Despite being reduced by Eimeria, these bacteria retained higher abundance in M5.1 chickens compared to the other breeds.

CONCLUSIONS

Fayoumi M5.1 chickens exhibit greater resistance to coccidiosis than Leghorn Ghs6 layers and Cobb broilers. Several commensal bacteria, including group A Lactobacillus, L. salivarius, Weissella, S. gallinarum, and E. durans/hirae, are differentially enriched in Fayoumi M5.1 chickens with strong correlation with coccidiosis resistance. These bacteria hold potential as probiotics for coccidiosis mitigation.

Different MHC class I cell surface expression levels in diverse chicken lines, associations with B blood group, and proposed relationship to antigen-binding repertoire

The Major Histocompatibility Complex (MHC) is a cluster of genes with primarily immune-related functions. The MHC class I genes are responsible for self- versus non-self-recognition and viral antigen presentation to T lymphocytes. The chicken MHC class I protein binds its cognate antigen(s) over a repertoire spectrum ranging from promiscuous (generalist) to fastidious (specialist). The MHC class I protein expression level at the cell surface is inversely related to the promiscuity of its peptide-binding repertoire. In our study, erythrocytes from 6 diverse and highly inbred lines of chickens, a closed broiler line, and a highly advanced intercross line were evaluated for MHC class I antigen expression level by flow cytometry using monoclonal antibodies to chicken MHC class I molecules. In chickens, the B blood group antigens include the MHC class I antigen expressed from the MHC. Thus, the B blood group has historically been used as a genetic marker for Marek's Disease virus response. Erythrocytes of the inbred lines were blood typed by serology. The B21 blood type is widely recognized as relatively resistant to Marek's disease and regarded as an MHC class I generalist with low MHC class I expression. The Spanish line, which types serologically as B21.1 (similar to B21), was the lowest MHC class I expressing line. The two sublines (B5.1 and B15.2) of the Fayoumi breed, which significantly differed in their MHC class I expression, also differ in response to multiple pathogens. These defined genetic lines of chickens, with distinct MHC class I expression levels, provide an excellent platform to further interrogate the hypothesis of high or low MHC class I expression (antigenic specialists or generalists, respectively) determining diverse responses to pathogens.

Integrative study of chicken lung transcriptome to understand the host immune response during Newcastle disease virus challenge

Introduction

Newcastle disease is one of the significant issues in the poultry industry, having catastrophic effects worldwide. The lung is one of the essential organs which harbours Bronchus-associated lymphoid tissue and plays a vital role in the immune response. Leghorn and Fayoumi breeds are known to have differences in resistance to Newcastle disease. Along with genes and long non-coding RNAs (lncRNAs) are also known to regulate various biological pathways through gene regulation.

Methods

This study analysed the lung transcriptome data and identified the role of genes and long non-coding RNAs in differential immune resistance. The computational pipeline, FHSpipe, as used in our previous studies on analysis of harderian gland and trachea transcriptome was used to identify genes and lncRNAs. This was followed by differential expression analysis, functional annotation of genes and lncRNAs, identification of transcription factors, microRNAs and finally validation using qRT-PCR.

Results and discussion

A total of 8219 novel lncRNAs were identified. Of them, 1263 lncRNAs and 281 genes were differentially expressed. About 66 genes were annotated with either an immune-related GO term or pathway, and 12 were annotated with both. In challenge and breed-based analysis, most of these genes were upregulated in Fayoumi compared to Leghorn, and in timepoint-based analysis, Leghorn challenge chicken showed downregulation between time points. A similar trend was observed in the expression of lncRNAs. Co-expression analysis has revealed several lncRNAs co-expressing with immune genes with a positive correlation. Several genes annotated with non-immune pathways, including metabolism, signal transduction, transport of small molecules, extracellular matrix organization, developmental biology and cellular processes, were also impacted. With this, we can understand that Fayoumi chicken showed upregulated immune genes and positive cis-lncRNAs during both the non-challenged and NDV-challenge conditions, even without viral transcripts in the tissue. This finding shows that these immune-annotated genes and coexpressing cis-lncRNAs play a significant role in Fayoumi being comparatively resistant to NDV compared to Leghorn. Our study affirms and expands upon the outcomes of previous studies and highlights the crucial role of lncRNAs during the immune response to NDV.

Conclusion

This analysis clearly shows the differences in the gene expression patterns and lncRNA co-expression with the genes between Leghorn and Fayoumi, indicating that the lncRNAs and co-expressing genes might potentially have a role in differentiating these breeds. We hypothesise that these genes and lncRNAs play a vital role in the higher resistance of Fayoumi to NDV than Leghorn. This study can pave the way for future studies to unravel the biological mechanism behind the regulation of immune-related genes.

Changes in the Transcriptome Profile in Young Chickens after Infection with LaSota Newcastle Disease Virus

Understanding gene expression changes in chicks after vaccination against Newcastle Disease (ND) can reveal vaccine biomarkers. There are limited data on chicks' early immune response after ND vaccination. Two trials focused on this knowledge gap. In experiment one, 42 13-day-old specific-pathogen-free (SPF) chicks were used. Harderian glands (Hgs) and tracheas (Tcs) from five birds per group were sampled at 12, 24, and 48 h post-vaccination (hpv) to evaluate the gene transcription levels by RNA sequencing (RNA-seq) and RT-qPCR. The results of RNA-seq were compared by glmFTest, while results of RT-qPCR were compared by t-test. With RNA-seq, a significant up-regulation of interferon-related genes along with JAK-STAT signaling pathway regulation was observed in the Hgs at 24 hpv. None of the differentially expressed genes (DEGs) identified by RNA-seq were positive for RT-qPCR. Experiment 2 used 112 SPF and commercial chickens that were 1 day old and 14 days old. Only the commercial birds had maternal antibodies for Newcastle Disease virus (NDV). By RNA-seq, 20 core DEGs associated with innate immunity and viral genome replication inhibition were identified. Genes previously unlinked to NDV response, such as USP41, were identified. This research present genes with potential as immunity biomarkers for vaccines, yet further investigation is needed to correlate the core gene expression with viral shedding post-vaccination.

Molecular Detection and Serological Investigation of Newcastle Disease in Intensive, Semi-Intensive, and Backyard Production Systems in Central and Southwestern Areas of Ethiopia

Purpose

The purpose of this research is to detect Newcastle disease virus and to assess the seropositivity among backyard, semi-intensive, and intensive farms located in central and southwestern areas of Ethiopia.

Material and Methods

A total of 239 oropharyngeal and cloacal swab samples were collected from symptomatic birds found in Holeta, Burayu, Jimma towns as well as Seka Chekorsa and Nadhigibe woredas of Jimma Zone. In addition, ninety blood samples were collected from wing veins of unvaccinated birds found in the study areas of Jimma zone. Side-by-side information related to risk factors estimated to contribute to the susceptibility of the disease was collected by interviewing owners of sampled birds. Reverse transcription polymerase-chain reaction (RT-PCR) was conducted to detect NDV. Likewise, Enzyme-linked immunosorbent assay (ELISA) was performed to determine the seropositivity of ND.

Results

The proportion of samples where NDV was detected was 24.6%. Similarly, 68.9% of the sampled birds were seropositive. It was observed that adult birds were more likely to encounter the disease than youngs (OR = 11.6; 95% CI: 4.0-33.3; P = 0.000). Birds owned by respondents who leave diseased birds in the flock were more likely infected (OR = 6.2; 95% CI: 1.8-21.2; P=0.004) as compared to those isolated and mode of disposal of dead chicken significantly affect exposure (OR = 0.13; 95% CI: 0.10-4.88; P = 0.044). Likewise, access to veterinary services highly likely reduces susceptibility to the disease (OR = 12.4; 95% CI: 3.2-46.9; P = 0.000). It was also found that birds farmed intensively were the most at risk (OR = 2.8; 95% CI: 0.58-13.71; P = 0.199).

Conclusion

Detection of ND from a significant proportion of sampled birds and their high seropositivity percentage revealed the circulation of the virus in the study areas.

Integrated analysis of genes and long non-coding RNAs in trachea transcriptome to decipher the host response during Newcastle disease challenge in different breeds of chicken

Newcastle disease is a highly infectious economically devastating disease caused by Newcastle disease Virus in Chicken (Gallus gallus). Leghorn and Fayoumi are two breeds which show differential resistance patterns towards NDV. This study aims to identify the differentially expressed genes and lncRNAs during NDV challenge which could play a potential role in this differential resistance pattern. A total of 552 genes and 1580 lncRNAs were found to be differentially expressing. Of them, 52 genes were annotated with both Immune related pathways and Gene ontologies. We found that most of these genes were upregulated in Leghorn between normal and challenged chicken but several were down regulated between different timepoints after NDV challenge, while Fayoumi showed no such downregulation. We also observed that higher number of positively correlating lncRNAs was found to be downregulated along with these genes. This shows that although Leghorn is showing higher number of differentially expressed genes in challenged than in non-challenged, most of them were downregulated during the disease between different timepoints. With this we hypothesize that the downregulation of immune related genes and co-expressing lncRNAs could play a significant role behind the Leghorn being comparatively susceptible breed than Fayoumi. The computational pipeline is available at https://github.com/Venky2804/FHSpipe.

Development and function of chicken XCR1+ conventional dendritic cells

Conventional dendritic cells (cDCs) are antigen-presenting cells (APCs) that play a central role in linking innate and adaptive immunity. cDCs have been well described in a number of different mammalian species, but remain poorly characterised in the chicken. In this study, we use previously described chicken cDC specific reagents, a novel gene-edited chicken line and single-cell RNA sequencing (scRNAseq) to characterise chicken splenic cDCs. In contrast to mammals, scRNAseq analysis indicates that the chicken spleen contains a single, chemokine receptor XCR1 expressing, cDC subset. By sexual maturity the XCR1+ cDC population is the most abundant mononuclear phagocyte cell subset in the chicken spleen. scRNAseq analysis revealed substantial heterogeneity within the chicken splenic XCR1+ cDC population. Immature MHC class II (MHCII)LOW XCR1+ cDCs expressed a range of viral resistance genes. Maturation to MHCIIHIGH XCR1+ cDCs was associated with reduced expression of anti-viral gene expression and increased expression of genes related to antigen presentation via the MHCII and cross-presentation pathways. To visualise and transiently ablate chicken XCR1+ cDCs in situ, we generated XCR1-iCaspase9-RFP chickens using a CRISPR-Cas9 knockin transgenesis approach to precisely edit the XCR1 locus, replacing the XCR1 coding region with genes for a fluorescent protein (TagRFP), and inducible Caspase 9. After inducible ablation, the chicken spleen is initially repopulated by immature CD1.1+ XCR1+ cDCs. XCR1+ cDCs are abundant in the splenic red pulp, in close association with CD8+ T-cells. Knockout of XCR1 prevented this clustering of cDCs with CD8+ T-cells. Taken together these data indicate a conserved role for chicken and mammalian XCR1+ cDCs in driving CD8+ T-cells responses.

Whole transcriptome analysis in offspring whose fathers were exposed to a developmental insult: a novel avian model

Although the effects of paternal exposure to insults on the offspring received limited attention in the past, it is currently gaining interest especially after understanding the mechanisms which may mediate such exposure effects. In the current study, the well-controlled avian model (Fayoumi) was utilized to investigate the effects of paternal exposure to the developmental insult, chlorpyrifos on the offspring's gene expression via mRNA and small RNA sequencing. Numerous mRNA gene expression changes were detected in the offspring after paternal exposure to the developmental insult, especially in genes related to neurogenesis, learning and memory. qPCR analysis of several genes, that were significantly changed in mRNA sequencing, confirmed the results obtained in mRNA sequencing. On the other hand, small RNA sequencing did not identify significant microRNA genes expression changes in the offspring after paternal exposure to the developmental insult. The effects of the paternal exposure were more pronounced in the female offspring compared to the male offspring. The results identified expression alterations in major genes (some of which were pertinent to the functional changes observed in other forms of early developmental exposure) after paternal insult exposure and provided a direction for future studies involving the most affected genes.

Comparison of Immune-Related Gene Expression in Two Chicken Breeds Following Infectious Bronchitis Virus Vaccination

This study aims to identify the immune-related genes and the corresponding biological pathways following infectious bronchitis virus vaccination in Taiwan Country and White Leghorn chicken breeds. Transcriptomic analyses of the spleen of these two breeds were conducted by next-generation sequencing. Compared to White Leghorn chicken, Taiwan Country chicken showed a significantly higher level of anti-infectious bronchitis virus (IBV) antibodies at 14 and 21 days pos vaccination. At 7 days post vaccination, in the Taiwan Country chicken, higher expression of mitogen-activated protein kinase 10, Major histocompatibility complex class 1, and V-set pre-B cell surrogate light chain 3 were found. In contrast, the White Leghorn chicken had a high expression of interleukin 4 induced 1, interleukin 6, and interleukin 22 receptor subunit alpha 2. These findings have highlighted the variations in immune induction between chickens with distinct genetic background and provided biological pathways and specific genes involved in immune responses against live attenuated IBV vaccine.

Parental Preconception and Pre-Hatch Exposure to a Developmental Insult Alters Offspring's Gene Expression and Epigenetic Regulations: An Avian Model

Parental exposure to insults was initially considered safe if stopped before conception. In the present investigation, paternal or maternal preconception exposure to the neuroteratogen chlorpyrifos was investigated in a well-controlled avian model (Fayoumi) and compared to pre-hatch exposure focusing on molecular alterations. The investigation included the analysis of several neurogenesis, neurotransmission, epigenetic and microRNA genes. A significant decrease in the vesicular acetylcholine transporter (SLC18A3) expression was detected in the female offspring in the three investigated models: paternal (57.7%, p < 0.05), maternal (36%, p < 0.05) and pre-hatch (35.6%, p < 0.05). Paternal exposure to chlorpyrifos also led to a significant increase in brain-derived neurotrophic factor (BDNF) gene expression mainly in the female offspring (27.6%, p < 0.005), while its targeting microRNA, miR-10a, was similarly decreased in both female (50.5%, p < 0.05) and male (56%, p < 0.05) offspring. Doublecortin's (DCX) targeting microRNA, miR-29a, was decreased in the offspring after maternal preconception exposure to chlorpyrifos (39.8%, p < 0.05). Finally, pre-hatch exposure to chlorpyrifos led to a significant increase in protein kinase C beta (PKCß; 44.1%, p < 0.05), methyl-CpG-binding domain protein 2 (MBD2; 44%, p < 0.01) and 3 (MBD3; 33%, p < 0.05) genes expression in the offspring. Although extensive studies are required to establish a mechanism-phenotype relationship, it should be noted that the current investigation does not include phenotype assessment in the offspring.

A comparative evaluation of transcriptome changes in lung and spleen tissues of chickens infected with velogenic and mesogenic Avian Orthoavulavirus 1

Newcastle disease is an acute, highly contagious disease responsible for severe economic losses to the poultry industry worldwide. Clinical assessment of different pathotypes of AOaV-1 strains is well-elucidated in chickens. However, a paucity of data exists for a comparative assessment of avian innate immune responses in birds after infection with two different pathotypes of AOaV-1. We compared early immune responses in chickens infected with a duck-originated velogenic strain (high virulent: genotype VII) and a pigeon-originated mesogenic stain (moderate virulent; genotype VI). Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) identified 4737 differentially expressed genes (DEGs) in the transcriptional profiles of lung and spleen tissues of chickens infected with both pathotypes. More DEGs were expressed in spleen tissue infected with velogenic strain compared to spleen or lung exposed to mesogenic strain. An enriched expression was observed for genes involved in metabolic processes and cellular components, including innate immune-associated signaling pathways. Most DEGs were involved in RIG-I, Toll-like, NF-Kappa B, and MAPK signaling pathways to activate interferon-stimulated genes (ISGs). This study provided a comparative insight into complicated molecular mechanisms and associated DEGs involved in early immune responses of birds to two different AOaV-1 strains.

Improved biomarker discovery through a plot twist in transcriptomic data analysis

Background

Transcriptomic analysis is crucial for understanding the functional elements of the genome, with the classic method consisting of screening transcriptomics datasets for differentially expressed genes (DEGs). Additionally, since 2005, weighted gene co-expression network analysis (WGCNA) has emerged as a powerful method to explore relationships between genes. However, an approach combining both methods, i.e., filtering the transcriptome dataset by DEGs or other criteria, followed by WGCNA (DEGs + WGCNA), has become common. This is of concern because such approach can affect the resulting underlying architecture of the network under analysis and lead to wrong conclusions. Here, we explore a plot twist to transcriptome data analysis: applying WGCNA to exploit entire datasets without affecting the topology of the network, followed with the strength and relative simplicity of DEG analysis (WGCNA + DEGs). We tested WGCNA + DEGs against DEGs + WGCNA to publicly available transcriptomics data in one of the most transcriptomically complex tissues and delicate processes: vertebrate gonads undergoing sex differentiation. We further validate the general applicability of our approach through analysis of datasets from three distinct model systems: European sea bass, mouse, and human.

Results

In all cases, WGCNA + DEGs clearly outperformed DEGs + WGCNA. First, the network model fit and node connectivity measures and other network statistics improved. The gene lists filtered by each method were different, the number of modules associated with the trait of interest and key genes retained increased, and GO terms of biological processes provided a more nuanced representation of the biological question under consideration. Lastly, WGCNA + DEGs facilitated biomarker discovery.

Conclusions

We propose that building a co-expression network from an entire dataset, and only thereafter filtering by DEGs, should be the method to use in transcriptomic studies, regardless of biological system, species, or question being considered.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01398-w.

Host response to successive challenges with lentogenic and velogenic Newcastle disease virus in local chickens of Ghana

Newcastle disease (ND) is a highly contagious viral disease that constantly threatens poultry production. The velogenic (highly virulent) form of ND inflicts the most damage and can lead to 100% mortality in unvaccinated village chicken flocks. This study sought to characterize responses of local chickens in Ghana after challenging them with lentogenic and velogenic Newcastle disease virus (NDV) strains. At 4 wk of age, chicks were challenged with lentogenic NDV. Traits measured were pre- and post-lentogenic infection growth rates (GR), viral load at 2 and 6 d post-lentogenic infection (DPI), viral clearance rate and antibody levels at 10 DPI. Subsequently, the chickens were naturally exposed to velogenic NDV (vNDV) after anti-NDV antibody titers had waned to levels ≤1:1,700. Body weights and blood samples were again collected for analysis. Finally, chickens were euthanized and lesion scores (LS) across tissues were recorded. Post-velogenic exposure GR; antibody levels at 21 and 34 days post-velogenic exposure (DPE); LS for trachea, proventriculus, intestines, and cecal tonsils; and average LS across tissues were measured. Variance components and heritabilities were estimated for all traits using univariate animal models. Mean pre- and post-lentogenic NDV infection GRs were 6.26 g/day and 7.93 g/day, respectively, but mean post-velogenic NDV exposure GR was −1.96 g/day. Mean lesion scores ranged from 0.52 (trachea) to 1.33 (intestine), with males having significantly higher (P < 0.05) lesion scores compared to females. Heritability estimates for the lentogenic NDV trial traits ranged from moderate (0.23) to high (0.55) whereas those for the vNDV natural exposure trial were very low (≤ 0.08). Therefore, in contrast to the vNDV exposure trial, differences in the traits measured in the lentogenic challenge were more affected by genetics and thus selection for these traits may be more feasible compared to those following vNDV exposure. Our results can form the basis for identifying local chickens with improved resilience in the face of NDV infection for selective breeding to improve productivity.

Slc6a13 Deficiency Attenuates Pasteurella multocida Infection-Induced Inflammation via Glycine-Inflammasome Signaling

We have previously demonstrated that Slc6a13-deficient (Slc6a13-/-; KO) mice are resistant to P. multocida infection, which might be in connection with macrophage-mediated inflammation; however, the specific metabolic mechanism is still enigmatic. Here we reproduce the less sensitive to P. multocida infection in overall survival assays as well as reduced bacterial loads, tissue lesions, and inflammation of lungs in KO mice. The transcriptome sequencing analysis of wild-type (WT) and KO mice shows a large number of differentially expressed genes that are enriched in amino acid metabolism by functional analysis. Of note, glycine levels are substantially increased in the lungs of KO mice with or without P. multocida infection in comparison to the WT controls. Interestingly, exogenous glycine supplementation alleviates P. multocida infection-induced inflammation. Mechanistically, glycine reduces the production of inflammatory cytokines in macrophages by blocking the activation of inflammasome (NALP1, NLRP3, NLRC4, AIM2, and Caspase-1). Together, Slc6a13 deficiency attenuates P. multocida infection through lessening the excessive inflammatory responses of macrophages involving glycine-inflammasome signaling.

Transcriptional atlas analysis from multiple tissues reveals the expression specificity patterns in beef cattle

BACKGROUND

A comprehensive analysis of gene expression profiling across tissues can provide necessary information for an in-depth understanding of their biological functions. We performed a large-scale gene expression analysis and generated a high-resolution atlas of the transcriptome in beef cattle.

RESULTS

Our transcriptome atlas was generated from 135 bovine tissues in adult beef cattle, covering 51 tissue types of major organ systems (e.g., muscular system, digestive system, immune system, reproductive system). Approximately 94.76% of sequencing reads were successfully mapped to the reference genome assembly ARS-UCD1.2. We detected a total of 60,488 transcripts, and 32% of them were not reported before. We identified 2654 housekeeping genes (HKGs) and 477 tissue-specific genes (TSGs) across tissues. Using weighted gene co-expression network analysis, we obtained 24 modules with 237 hub genes (HUBGs). Functional enrichment analysis showed that HKGs mainly maintain the basic biological activities of cells, while TSGs were involved in tissue differentiation and specific physiological processes. HKGs in bovine tissues were more conserved in terms of expression pattern as compared to TSGs and HUBGs among multiple species. Finally, we obtained a subset of tissue-specific differentially expressed genes (DEGs) between beef and dairy cattle and several functional pathways, which may be involved in production and health traits.

CONCLUSIONS

We generated a large-scale gene expression atlas across the major tissues in beef cattle, providing valuable information for enhancing genome assembly and annotation. HKGs, TSGs, and HUBGs further contribute to better understanding the biology and evolution of multiple tissues in cattle. DEGs between beef and dairy cattle also fill in the knowledge gaps about differential transcriptome regulation of bovine tissues underlying economically important traits.

Differential gene expression reveals host factors for viral shedding variation in mallards (Anas platyrhynchos) infected with low-pathogenic avian influenza virus

Intraspecific variation in pathogen shedding impacts disease transmission dynamics; therefore, understanding the host factors associated with individual variation in pathogen shedding is key to controlling and preventing outbreaks. In this study, ileum and bursa of Fabricius tissues of wild-bred mallards (Anas platyrhynchos) infected with low-pathogenic avian influenza (LPAIV) were evaluated at various post-infection time points to determine genetic host factors associated with intraspecific variation in viral shedding. By analysing transcriptome sequencing data (RNA-seq), we found that LPAIV-infected wild-bred mallards do not exhibit differential gene expression compared to uninfected birds, but that gene expression was associated with cloacal viral shedding quantity early in the infection. In both tissues, immune gene expression was higher in high/moderate shedding birds compared to low shedding birds, and significant positive relationships with viral shedding were observed. In the ileum, expression for host genes involved in viral cell entry was lower in low shedders compared to moderate shedders at 1 day post-infection (DPI), and expression for host genes promoting viral replication was higher in high shedders compared to low shedders at 2 DPI. Our findings indicate that viral shedding is a key factor for gene expression differences in LPAIV-infected wild-bred mallards, and the genes identified in this study could be important for understanding the molecular mechanisms driving intraspecific variation in pathogen shedding.

Genome-Wide Analysis of Alternative Splicing during Host-Virus Interactions in Chicken

The chicken is a model animal for the study of evolution, immunity and development. In addition to their use as a model organism, chickens also represent an important agricultural product. Pathogen invasion has already been shown to modulate the expression of hundreds of genes, but the role of alternative splicing in avian virus infection remains unclear. We used RNA-seq data to analyze virus-induced changes in the alternative splicing of Gallus gallus, and found that a large number of alternative splicing events were induced by virus infection both in vivo and in vitro. Virus-responsive alternative splicing events preferentially occurred in genes involved in metabolism and transport. Many of the alternatively spliced transcripts were also expressed from genes with a function relating to splicing or immune response, suggesting a potential impact of virus infection on pre-mRNA splicing and immune gene regulation. Moreover, exon skipping was the most frequent AS event in chickens during virus infection. This is the first report describing a genome-wide analysis of alternative splicing in chicken and contributes to the genomic resources available for studying host-virus interaction in this species. Our analysis fills an important knowledge gap in understanding the extent of genome-wide alternative splicing dynamics occurring during avian virus infection and provides the impetus for the further exploration of AS in chicken defense signaling and homeostasis.

Heritability Estimate for Antibody Response to Vaccination and Survival to a Newcastle Disease Infection of Native chicken in a Low-Input Production System

The Newcastle disease virus (NDV) is the deadliest chicken pathogen in low-input village poultry, and selecting for NDV resistance has been recommended as a sustainable strategy in backyard poultry production systems. However, selecting for disease resistance needs precision data from either a big population sample size or on many generations with good pedigree records for effective prediction of heritability (h²) and breeding values of the foundation stock. Such conditions are almost impossible to meet in low-input backyard production systems. This study aimed at proposing a realistic method for estimating the heritability of the immune response to vaccination and survival of NDV infection in village poultry production to inform a breeding strategy for ND resistance in Cameroon. A 1 and 3% selection intensity of cocks and hens for higher antibody (ab) response (ABR) to vaccination followed by progeny selection of chickens who survived an experimental NDV infection was conducted from an initial population of 1,702 chickens. The selection induced an increase of 1012.47units/ml (p<0.01) of the NDV antibody of the progeny as well as an effective survival rate (ESR) increase of 11.75%. Three methods were used to estimate the heritability (h2) of NDV antibody response to vaccination. h2 was low irrespective of the method with estimates of 0.2227, 0.2442, and 0.2839 for the breeder’s equation method, the graphical method, and the full-sib/half-sib nested design, respectively. The mortality rate of infected chickens was high (86%). The antibody response to selection was not influenced by sex and genetic type even though the opposite was observed (p<0.05) for the ESR to NDV infection with naked neck chickens recording an ESR of 14% against 2.25% for the normal feather type. A very low heritability (0.0891) was observed for the survival against NDV infection. We confirm the evidence of disease resistance and the effect of selection for antibody response to vaccination on the improvement of the survival against NDV disease. Although the full sib/half sib nested design is more appropriate in case of availability of pedigree information, the direct methods are still useful in case of unavailability of full pedigree information. It is recommended that gene expression analysis should be prioritized for disease-resistance assessment and selection of native breeds of poultry.

Indicators of the molecular pathogenesis of virulent Newcastle disease virus in chickens revealed by transcriptomic profiling of spleen

Newcastle disease virus (NDV) has caused significant outbreaks in South-East Asia, particularly in Indonesia in recent years. Recently emerged genotype VII NDVs (NDV-GVII) have shifted their tropism from gastrointestinal/respiratory tropism to a lymphotropic virus, invading lymphoid organs including spleen and bursa of Fabricius to cause profound lymphoid depletion. In this study, we aimed to identify candidate genes and biological pathways that contribute to the disease caused by this velogenic NDV-GVII. A transcriptomic analysis based on RNA-Seq of spleen was performed in chickens challenged with NDV-GVII and a control group. In total, 6361 genes were differentially expressed that included 3506 up-regulated genes and 2855 down-regulated genes. Real-Time PCR of ten selected genes validated the RNA-Seq results as the correlation between them is 0.98. Functional and network analysis of Differentially Expressed Genes (DEGs) showed altered regulation of ElF2 signalling, mTOR signalling, proliferation of cells of the lymphoid system, signalling by Rho family GTPases and synaptogenesis signalling in spleen. We have also identified modified expression of IFIT5, PI3K, AGT and PLP1 genes in NDV-GVII infected chickens. Our findings in activation of autophagy-mediated cell death, lymphotropic and synaptogenesis signalling pathways provide new insights into the molecular pathogenesis of this newly emerged NDV-GVII.

Transcriptomic analysis to infer key molecular players involved during host response to NDV challenge in Gallus gallus (Leghorn & Fayoumi)

Long non-coding RNAs (lncRNAs) are the transcripts of length longer than 200 nucleotides. They are involved in the regulation of various biological activities. Leghorn and Fayoumi breeds of Gallus gallus were known to be having differential resistance against Newcastle Disease Virus (NDV) infection. Differentially expressed genes which were thought to be involved in this pattern of resistance were already studied. Here we report the analysis of the transcriptomic data of Harderian gland of Gallus gallus for studying the lncRNAs involved in regulation of these genes. Using bioinformatics approaches, a total of 37,411 lncRNAs were extracted and 359 lncRNAs were differentially expressing. Functional annotation using co-expression analysis revealed the involvement of lncRNAs in the regulation of various pathways. We also identified 1232 quantitative trait loci (QTLs) associated with the genes interacting with lncRNA. Additionally, we identified the role of lncRNAs as putative micro RNA precursors, and the interaction of differentially expressed Genes with transcription factors and micro RNAs. Our study revealed the role of lncRNAs during host response against NDV infection which would facilitate future experiments in unravelling regulatory mechanisms of development in the genetic improvement of the susceptible breeds of Gallus gallus.

Distinct transcriptomic response to Newcastle disease virus infection during heat stress in chicken tracheal epithelial tissue

Newcastle disease (ND) has a great impact on poultry health and welfare with its most virulent (velogenic) strain. In addition, issues exacerbated by the increase in global temperatures necessitates a greater understanding of the host immune response when facing a combination of biotic and abiotic stress factors in poultry production. Previous investigations have revealed that the host immune response is tissue-specific. The goal of this study was to identify genes and/or signaling pathways associated with immune response to NDV (Newcastle disease virus) in the trachea, an essential organ where NDV replicate after the infection, by profiling the tissue specific transcriptome response in two genetically distinct inbred chicken lines when exposed to both abiotic and biotic stressors. Fayoumis appear to be able to respond more effectively (lower viral titer, higher antibody levels, immune gene up-regulation) and earlier than Leghorns. Our results suggest NDV infection in Fayoumis appears to elicit proinflammatory processes, and pathways such as the inhibition of cell viability, cell proliferation of lymphocytes, and transactivation of RNA, more rapidly than in Leghorns. These differences in immune response converge at later timepoints which may indicate that Leghorns eventually regulate its immune response to infection. The profiling of the gene expression response in the trachea adds to our understanding of the chicken host response to NDV infection and heat stress on a whole genome level and provides potential candidate genes and signaling pathways for further investigation into the characterization of the time-specific and pathway specific responses in Fayoumis and Leghorns.

Genomic Diversity and Evolution of Quasispecies in Newcastle Disease Virus Infections

Newcastle disease virus (NDV) infections are well known to harbour quasispecies, due to the error-prone nature of the RNA polymerase. Quasispecies variants in the fusion cleavage site of the virus are known to significantly change its virulence. However, little is known about the genomic patterns of diversity and selection in NDV viral swarms. We analyse deep sequencing data from in vitro and in vivo NDV infections to uncover the genomic patterns of diversity and the signatures of selection within NDV swarms. Variants in viruses from in vitro samples are mostly localised in non-coding regions and 3′ and 5′ untranslated regions (3′UTRs or 5′UTRs), while in vivo samples contain an order of magnitude more variants. We find different patterns of genomic divergence and diversity among NDV genotypes, as well as differences in the genomic distribution of intra-host variants among in vitro and in vivo infections of the same strain. The frequency spectrum shows clear signatures of intra-host purifying selection in vivo on the matrix protein (M) coding gene and positive or diversifying selection on nucleocapsid (NP) and haemagglutinin-neuraminidase (HN). The comparison between within-host polymorphisms and phylogenetic divergence reveals complex patterns of selective pressure on the NDV genome at between- and within-host level. The M sequence is strongly constrained both between and within hosts, fusion protein (F) coding gene is under intra-host positive selection, and NP and HN show contrasting patterns: HN RNA sequence is positively selected between hosts while its protein sequence is positively selected within hosts, and NP is under intra-host positive selection at the RNA level and negative selection at the protein level.

Patterns of RNA Editing in Newcastle Disease Virus Infections

The expression of accessory non-structural proteins V and W in Newcastle disease virus (NDV) infections depends on RNA editing. These proteins are derived from frameshifts of the sequence coding for the P protein via co-transcriptional insertion of one or two guanines in the mRNA. However, a larger number of guanines can be inserted with lower frequencies. We analysed data from deep RNA sequencing of samples from in vitro and in vivo NDV infections to uncover the patterns of mRNA editing in NDV. The distribution of insertions is well described by a simple Markov model of polymerase stuttering, providing strong quantitative confirmation of the molecular process hypothesised by Kolakofsky and collaborators three decades ago. Our results suggest that the probability that the NDV polymerase would stutter is about 0.45 initially, and 0.3 for further subsequent insertions. The latter probability is approximately independent of the number of previous insertions, the host cell, and viral strain. However, in LaSota infections, we also observe deviations from the predicted V/W ratio of about 3:1 according to this model, which could be attributed to deviations from this stuttering model or to further mechanisms downregulating the abundance of W protein.

Transcriptome Analysis Reveals Inhibitory Effects of Lentogenic Newcastle Disease Virus on Cell Survival and Immune Function in Spleen of Commercial Layer Chicks

As a major infectious disease in chickens, Newcastle disease virus (NDV) causes considerable economic losses in the poultry industry, especially in developing countries where there is limited access to effective vaccination. Therefore, enhancing resistance to the virus in commercial chickens through breeding is a promising way to promote poultry production. In this study, we investigated gene expression changes at 2 and 6 days post inoculation (dpi) at day 21 with a lentogenic NDV in a commercial egg-laying chicken hybrid using RNA sequencing analysis. By comparing NDV-challenged and non-challenged groups, 526 differentially expressed genes (DEGs) (false discovery rate (FDR) < 0.05) were identified at 2 dpi, and only 36 at 6 dpi. For the DEGs at 2 dpi, Ingenuity Pathway Analysis predicted inhibition of multiple signaling pathways in response to NDV that regulate immune cell development and activity, neurogenesis, and angiogenesis. Up-regulation of interferon induced protein with tetratricopeptide repeats 5 (IFIT5) in response to NDV was consistent between the current and most previous studies. Sprouty RTK signaling antagonist 1 (SPRY1), a DEG in the current study, is in a significant quantitative trait locus associated with virus load at 6 dpi in the same population. These identified pathways and DEGs provide potential targets to further study breeding strategy to enhance NDV resistance in chickens.

Transcriptomic Analysis of Chicken Lungs Infected With Avian and Bovine Pasteurella multocida Serotype A

Pasteurella multocida (P. multocida) is a common animal pathogen responsible for many animal diseases. Strains from different hosts exhibit disparate degrees of effect in other species. Here, we characterize an avian P. multocida serogroup A strain (PmQ) showing high lethality to chickens and a bovine P. multocida serogroup A strain (PmCQ2) with no lethality to chickens. We used RNA-seq to profile the transcriptomes of chicken lungs infected with PmQ and PmCQ2. A total of 1,649 differentially expressed genes (DEGs) due to PmQ infection (831 upregulated genes and 818 downregulated genes) and 1427 DEGs (633 upregulated genes and 794 downregulated genes) due to PmCQ2 infection were identified. Functional analysis of these DEGs demonstrated that the TNF signaling pathway, the toll-like receptor signaling pathway, complement and coagulation cascades, and cytokine–cytokine receptor interaction were both enriched in PmQ and PmCQ2 infection. STAT and apoptosis signaling pathways were uniquely enriched by PmQ infection, and the NOD-like receptor signaling pathway was enriched only by PmCQ2 infection. Cell-type enrichment analysis of the transcriptomes showed that immune cells, including macrophages and granulocytes, were enriched in both infection groups. Collectively, our study profiled the transcriptomic response of chicken lungs infected with P. multocida and provided valuable information to understand the chicken responses to P. multocida infection.

Usutu Virus Infection of Embryonated Chicken Eggs and a Chicken Embryo-Derived Primary Cell Line

Usutu virus (USUV) is a mosquito-borne flavivirus, closely related to the West Nile virus (WNV). Similar to WNV, USUV may cause infections in humans, with occasional, but sometimes severe, neurological complications. Further, USUV can be highly pathogenic in wild and captive birds and its circulation in Europe has given rise to substantial avian death. Adequate study models of this virus are still lacking but are critically needed to understand its pathogenesis and virulence spectrum. The chicken embryo is a low-cost, easy-to-manipulate and ethically acceptable model that closely reflects mammalian fetal development and allows immune response investigations, drug screening, and high-throughput virus production for vaccine development. While former studies suggested that this model was refractory to USUV infection, we unexpectedly found that high doses of four phylogenetically distinct USUV strains caused embryonic lethality. By employing immunohistochemistry and quantitative reverse transcriptase-polymerase chain reaction, we demonstrated that USUV was widely distributed in embryonic tissues, including the brain, retina, and feather follicles. We then successfully developed a primary cell line from the chorioallantoic membrane that was permissive to the virus without the need for viral adaptation. We believe the future use of these models would foster a significant understanding of USUV-induced neuropathogenesis and immune response and allow the future development of drugs and vaccines against USUV.

Novel Combined Tissue Transcriptome Analysis After Lentogenic Newcastle Disease Virus Challenge in Inbred Chicken Lines of Differential Resistance

Disease has large negative impacts on poultry production. A more comprehensive understanding of host-pathogen interaction can lead to new and improved strategies to maintain health. In particular, host genetic factors can lead to a more effective response to pathogens, hereafter termed resistance. Fayoumi and Leghorn chicken lines have demonstrated relative resistance and susceptibility, respectively, to the Newcastle disease virus (NDV) vaccine strain and many other pathogens. This biological model was used to better understand the host response to a vaccine strain of NDV across three tissues and time points, using RNA-seq. Analyzing the Harderian gland, trachea, and lung tissues together using weighted gene co-expression network analysis (WGCNA) identified important genes that were co-expressed and associated with parameters including: genetic line, days post-infection (dpi), challenge status, sex, and tissue. Pathways and driver genes, such as EIF2AK2, MPEG1, and TNFSF13B, associated with challenge status, dpi, and genetic line were of particular interest as candidates for disease resistance. Overall, by jointly analyzing the three tissues, this study identified genes and gene networks that led to a more comprehensive understanding of the whole animal response to lentogenic NDV than that obtained by analyzing the tissues individually.

Development of innate immunity in chicken embryos and newly hatched chicks: a disease control perspective

Newly hatched chickens are confronted by a wide array of pathogenic microbes because their adaptive immune defences have limited capabilities to control these pathogens. In such circumstances, and within this age group, innate responses provide a degree of protection. Moreover, as the adaptive immune system is relatively naïve to foreign antigens, synergy with innate defences is critical. This review presents knowledge on the ontogeny of innate immunity in chickens pre-hatch and early post-hatch and provides insights into possible interventions to modulate innate responses early in the life of the bird. As in other vertebrate species, the chicken innate immune system which include cellular mediators, cytokine and chemokine repertoires and molecules involved in antigen detection, develop early in life. Comparison of innate immune systems in newly hatched chickens and mature birds has revealed differences in magnitude and quality, but responses in younger chickens can be boosted using innate immune system modulators. Functional expression of pattern recognition receptors and several defence molecules by innate immune system cells of embryos and newly hatched chicks suggests that innate responses can be modulated at this stage of development to combat pathogens. Improved understanding of innate immune system ontogeny and functionality in chickens is critical for the implementation of sound and safe interventions to provide long-term protection against pathogens. Next-generation tools for studying genetic and epigenetic regulation of genes, functional metagenomics and gene knockouts can be used in the future to explore and dissect the contributions of signalling pathways of innate immunity and to devise more efficacious disease control strategies.

Insights into the chicken bursa of fabricius response to Newcastle disease virus at 48 and 72 hours post-infection through RNA-seq

Newcastle disease virus (NDV) causes significant economic losses to the poultry industry worldwide. As a lymphoid organ, the bursa of Fabricius (BF) plays a pivotal role in destroying invading pathogens. Virulent NDV strains can cause rapid atrophy of the BF; however, there is limited knowledge regarding the BF innate immune response to NDV infection. In this study, we used the virulent NDV strain F48E9 to infect four-week-old chickens and found atrophy of the BF, with severe damage and high NDV viral loads after NDV infection in dying chickens. To better understand the interactions between the host and NDV, we compared the transcriptional profiles at 48 and 72 h following infection with the virulent NDV strain F48E9 using RNA-seq. We identified a total of 1498 differentially expressed genes (DEGs), which were enriched in a variety of biological processes and pathways according to Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses. The enriched pathways were associated with innate immune and inflammatory responses as well as metabolism-related signalling pathways. Excessive inflammatory and innate immune responses induced by the NDV strain may be related to severe BF damage. The global survey of changes in gene expression performed herein provides new insights into complicated molecular mechanisms underlying the interaction between NDV and chickens and will enable the use of new strategies to protect chickens against NDV.

Specific genes related to nucleopolyhedrovirus in Bombyx mori susceptible and near-isogenic resistant strains through comparative transcriptome analysis

Bombyx mori nucleopolyhedrovirus (BmNPV) is one of the primary pathogens that causes severe economic losses to sericulture. Comparative transcriptomics analysis has been widely applied to explore the antiviral mechanism in resistant strains. Here, to identify genes involved in BmNPV infection, we identified differentially expressed genes (DEGs) and performed weighted gene co-expression network analysis (WGCNA) between two Bombyx mori strains: strain 871 (susceptible to BmNPV infection) and the near-isogenic strain 871C (resistant to BmNPV). Our results showed that 400 genes were associated with resistance in strain 871C, and 76 genes were related to susceptibility in strain 871. In addition, the correlation analysis of DEGs and WGCNA showed that 40 genes related to resistance were highly expressed in the resistant strain. Among them, gene BGIBMGA004291 was the most noticeable. We further identified the effect of gene BGIBMGA004291, which encoded a multiprotein bridge factor 2 (MBF2) family member (MBF2-10), on viral infection in cells. Our data suggested that MBF2-10 inhibited viral infection. Taken together, this study showed specific module trait correlations related to viral infection in strains 871 and 871C, and we identified a resistance-related gene. These findings suggested promising candidate genes with antiviral activity, aiding in the analysis of the antiviral molecular mechanisms in resistant strains.

Innate Immune Genes Associated With Newcastle Disease Virus Load in Chick Embryos From Inbred and Outbred Lines

Newcastle disease virus (NDV) causes substantial economic losses to smallholder farmers in low- and middle-income countries with high levels of morbidity and mortality in poultry flocks. Previous investigations have suggested differing levels of susceptibility to NDV between specific inbred lines and amongst breeds of chickens, however, the mechanisms contributing to this remain poorly understood. Studies have shown that some of these differences in levels of susceptibility to NDV infection may be accounted for by variability in the innate immune response amongst various breeds of poultry to NDV infection. Recent studies, in inbred Fayoumi and Leghorn lines, uncovered conserved, breed-dependent, and subline-dependent responses. To better understand the role of innate immune genes in engendering a protective immune response, we assessed the transcriptional responses to NDV of three highly outbred Tanzanian local chicken ecotypes, the Kuchi, the Morogoro Medium, and the Ching’wekwe. Hierarchical clustering and principal coordinate analysis of the gene expression profiles of 21-day old chick embryos infected with NDV clustered in an ecotype-dependent manner and was consistent with the relative viral loads for each of the three ecotypes. The Kuchi and Morogoro Medium exhibit significantly higher viral loads than the Ching’wekwe. The results show that the outbred ecotypes with increased levels of expression of CCL4, NOS2, and SOCS1 also had higher viral loads. The higher expression of SOCS1 is inconsistent with the expression in inbred lines. These differences may uncover new mechanisms or pathways in these populations that may have otherwise been overlooked when examining the response in highly inbred lines. Taken together, our findings provide insights on the specific conserved and differentially expressed innate immune-related genes involved the response of highly outbred chicken lines to NDV. This also suggests that several of the specific innate immunity related genes identified in the current investigation may serve as markers for the selection of chickens with reduced susceptibility to NDV.

Conserved, breed-dependent, and subline-dependent innate immune responses of Fayoumi and Leghorn chicken embryos to Newcastle disease virus infection

Newcastle disease virus (NDV) is a threat to the global poultry industry, but particularly for smallholder farmers in low- and middle-income countries. Previous reports suggest that some breeds of chickens are less susceptible to NDV infection, however, the mechanisms contributing to this are unknown. We here examined the comparative transcriptional responses of innate immune genes to NDV infection in inbred sublines of the Fayoumi and Leghorn breeds known to differ in their relative susceptibility to infection as well as at the microchromosome bearing the major histocompatability complex (MHC) locus. The analysis identified a set of five core genes, Mx1, IRF1, IRF7, STAT1, and SOCS1, that are up-regulated regardless of subline. Several genes were differentially expressed in a breed- or subline-dependent manner. The breed-dependent response involved TLR3, NOS2, LITAF, and IFIH1 in the Fayoumi versus IL8, CAMP, and CCL4 in the Leghorn. Further analysis identified subline-dependent differences in the pro-inflammatory response within the Fayoumi breed that are likely influenced by the MHC. These results have identified conserved, breed-dependent, and subline-dependent innate immune responses to NDV infection in chickens, and provide a strong framework for the future characterization of the specific roles of genes and pathways that influence the susceptibility of chickens to NDV infection.

Genetics and Genomic Regions Affecting Response to Newcastle Disease Virus Infection under Heat Stress in Layer Chickens

Newcastle disease virus (NDV) is a highly contagious avian pathogen that poses a tremendous threat to poultry producers in endemic zones due to its epidemic potential. To investigate host genetic resistance to NDV while under the effects of heat stress, a genome-wide association study (GWAS) was performed on Hy-Line Brown layer chickens that were challenged with NDV while under high ambient temperature to identify regions associated with host viral titer, circulating anti-NDV antibody titer, and body weight change. A single nucleotide polymorphism (SNP) on chromosome 1 was associated with viral titer at two days post-infection (dpi), while 30 SNPs spanning a quantitative trait loci (QTL) on chromosome 24 were associated with viral titer at 6 dpi. Immune related genes, such as CAMK1d and CCDC3 on chromosome 1, associated with viral titer at 2 dpi, and TIRAP, ETS1, and KIRREL3, associated with viral titer at 6 dpi, were located in two QTL regions for viral titer that were identified in this study. This study identified genomic regions and candidate genes that are associated with response to NDV during heat stress in Hy-Line Brown layer chickens. Regions identified for viral titer on chromosome 1 and 24, at 2 and 6 dpi, respectively, included several genes that have key roles in regulating the immune response.

Immune genes are hotspots of shared positive selection across birds and mammals

Consistent patterns of positive selection in functionally similar genes can suggest a common selective pressure across a group of species. We use alignments of orthologous protein-coding genes from 39 species of birds to estimate parameters related to positive selection for 11,000 genes conserved across birds. We show that functional pathways related to the immune system, recombination, lipid metabolism, and phototransduction are enriched for positively selected genes. By comparing our results with mammalian data, we find a significant enrichment for positively selected genes shared between taxa, and that these shared selected genes are enriched for viral immune pathways. Using pathogen-challenge transcriptome data, we show that genes up-regulated in response to pathogens are also enriched for positively selected genes. Together, our results suggest that pathogens, particularly viruses, consistently target the same genes across divergent clades, and that these genes are hotspots of host-pathogen conflict over deep evolutionary time.

Integrated Proteomic and Transcriptomic Analysis of Differential Expression of Chicken Lung Tissue in Response to NDV Infection during Heat Stress

Newcastle disease virus (NDV) is a devastating worldwide poultry pathogen with major implications for global food security. In this study, two highly inbred and genetically distinct chicken lines, Fayoumis and Leghorns, were exposed to a lentogenic strain of NDV, while under the effects of heat stress, in order to understand the genetic mechanisms of resistance during high ambient temperatures. Fayoumis, which are relatively more resistant to pathogens than Leghorns, had larger numbers of differentially expressed genes (DEGs) during the early stages of infection when compared to Leghorns and subsequently down-regulated their immune response at the latter stages to return to homeostasis. Leghorns had very few DEGs across all observed time points, with the majority of DEGs involved with metabolic and glucose-related functions. Proteomic analysis corroborates findings made within Leghorns, while also identifying interesting candidate genes missed by expression profiling. Poor correlation between changes observed in the proteomic and transcriptomic datasets highlights the potential importance of integrative approaches to understand the mechanisms of disease response. Overall, this study provides novel insights into global protein and expression profiles of these two genetic lines, and provides potential genetic targets involved with NDV resistance during heat stress in poultry.

Association of Candidate Genes with Response to Heat and Newcastle Disease Virus

Newcastle disease is considered the number one disease constraint to poultry production in low and middle-income countries, however poultry that is raised in resource-poor areas often experience multiple environmental challenges. Heat stress has a negative impact on production, and immune response to pathogens can be negatively modulated by heat stress. Candidate genes and regions chosen for this study were based on previously reported associations with response to immune stimulants, pathogens, or heat, including: TLR3, TLR7, MX, MHC-B (major histocompatibility complex, gene complex), IFI27L2, SLC5A1, HSPB1, HSPA2, HSPA8, IFRD1, IL18R1, IL1R1, AP2A2, and TOLLIP. Chickens of a commercial egg-laying line were infected with a lentogenic strain of NDV (Newcastle disease virus); half the birds were maintained at thermoneutral temperature and the other half were exposed to high ambient temperature before the NDV challenge and throughout the remainder of the study. Phenotypic responses to heat, to NDV, or to heat + NDV were measured. Selected SNPs (single nucleotide polymorphisms) within 14 target genes or regions were genotyped; and genotype effects on phenotypic responses to NDV or heat + NDV were tested in each individual treatment group and the combined groups. Seventeen significant haplotype effects, among seven genes and seven phenotypes, were detected for response to NDV or heat or NDV + heat. These findings identify specific genetic variants that are associated with response to heat and/or NDV which may be useful in the genetic improvement of chickens to perform favorably when faced with pathogens and heat stress.

Novel insights into the host immune response of chicken Harderian gland tissue during Newcastle disease virus infection and heat treatment

Background

Newcastle disease virus, in its most pathogenic form, threatens the livelihood of rural poultry farmers where there is a limited infrastructure and service for vaccinations to prevent outbreaks of the virus. Previously reported studies on the host response to Newcastle disease in chickens have not examined the disease under abiotic stressors, such as heat, which commonly experienced by chickens in regions such as Africa. The objective of this study was to elucidate the underlying biological mechanisms that contribute to disease resistance in chickens to the Newcastle disease virus while under the effects of heat stress.

Results

Differential gene expression analysis identified genes differentially expressed between treated and non-treated birds across three time points (2, 6, and 10 days post-infection) in Fayoumi and Leghorn birds. Across the three time points, Fayoumi had very few genes differentially expressed between treated and non-treated groups at 2 and 6 days post-infection. However, 202 genes were differentially expressed at 10 days post-infection. Alternatively, Leghorn had very few genes differentially expressed at 2 and 10 days post-infection but had 167 differentially expressed genes at 6 days post-infection. Very few differentially expressed genes were shared between the two genetic lines, and pathway analysis found unique signaling pathways specific to each genetic line. Fayoumi had significantly lower viral load, higher viral clearance, higher anti-NDV antibody levels, and fewer viral transcripts detected compared to Leghorns. Fayoumis activated immune related pathways including SAPK/JNK and p38 MAPK signaling pathways at earlier time points, while Leghorn would activate these same pathways at a later time. Further analysis revealed activation of the GP6 signaling pathway that may be responsible for the susceptible Leghorn response.

Conclusions

The findings in this study confirmed our hypothesis that the Fayoumi line was more resistant to Newcastle disease virus infection compared to the Leghorn line. Within line and interaction analysis demonstrated substantial differences in response patterns between the two genetic lines that was not observed from the within line contrasts. This study has provided novel insights into the transcriptome response of the Harderian gland tissue during Newcastle disease virus infection while under heat stress utilizing a unique resistant and susceptible model.

Electronic supplementary material

The online version of this article (10.1186/s12917-018-1583-0) contains supplementary material, which is available to authorized users.

Genetic Analysis of a Commercial Egg Laying Line Challenged With Newcastle Disease Virus

In low income countries, chickens play a vital role in daily life. They provide a critical source of protein through egg production and meat. Newcastle disease, caused by avian paramyxovirus type 1, has been ranked as the most devastating disease for scavenging chickens in Africa and Asia. High mortality among flocks infected with velogenic strains leads to a devastating loss of dietary protein and buying power for rural households. Improving the genetic resistance of chickens to Newcastle Disease virus (NDV), in addition to vaccination, is a practical target for improvement of poultry production in low income countries. Because response to NDV has a component of genetic control, it can be influenced through selective breeding. Adding genomic information to a breeding program can increase the amount of genetic progress per generation. In this study, we challenged a commercial egg-laying line with a lentogenic strain of NDV, measured phenotypic responses, collected genotypes, and associated genotypes with phenotypes. Collected phenotypes included viral load at 2 and 6 days post-infection (dpi), antibody levels pre-challenge and 10 dpi, and growth rates pre- and post-challenge. Six suggestive QTL associated with response to NDV and/or growth were identified, including novel and known QTL confirming previously reported associations with related traits. Additionally, previous RNA-seq analysis provided support for several of the genes located in or near the identified QTL. Considering the trend of negative genetic correlation between antibody and Newcastle Disease tolerance (growth under disease) and estimates of moderate to high heritability, we provide evidence that these NDV response traits can be influenced through selective breeding. Producing chickens that perform favorably in challenging environments will ultimately increase the supply of quality protein for human consumption.

What Makes the Harderian Gland Transcriptome Different From Other Chicken Immune Tissues? A Gene Expression Comparative Analysis

The Harderian gland is a sparsely characterized immune tissue known to play an important role in local immunity. The function of the Harderian gland, however, is not clearly defined. Measuring the expression of all genes using RNA-seq enables the identification of genes, pathways, or networks of interest. Our relative RNA-seq expression analysis compared the chicken Harderian gland transcriptome to other important primary and secondary immune tissues including the bursa of Fabricius, thymus, and spleen of non-challenged birds. A total of 2,386 transcripts were identified as highly expressed in the Harderian gland. Gene set enrichment showed the importance of G-protein coupled receptor signaling and several immune pathways. Among the genes highly expressed in the Harderian gland were 48 miRNAs, a category of genetic elements involved in regulation of gene expression. Several identified miRNAs have immune related functions. This analysis gives insight to the unique immune processes inherent in the Harderian gland.

Novel analysis of the Harderian gland transcriptome response to Newcastle disease virus in two inbred chicken lines

Behind each eye of the chicken resides a unique lymph tissue, the Harderian gland, for which RNA sequencing (RNA-seq) analysis is novel. We characterized the response of this tissue to Newcastle disease virus (NDV) in two inbred lines with different susceptibility to NDV across three time points. Three-week-old relatively resistant (Fayoumi) and relatively susceptible (Leghorn) birds were inoculated with a high-titered (10⁷EID₅₀) La Sota strain of NDV via an oculonasal route. At 2, 6, and 10 days post infection (dpi) Harderian glands were collected and analyzed via RNA-seq. The Fayoumi had significantly more detectable viral transcripts in the Harderian gland at 2 dpi than the Leghorn, but cleared the virus by 6 dpi. At all three time points, few genes were declared differentially expressed (DE) between the challenged and nonchallenged birds, except for the Leghorns at 6 dpi, and these DE genes were predicted to activate an adaptive immune response. Relative to the Leghorn, the Fayoumi was predicted to activate more immune pathways in both challenged and nonchallenged birds suggesting a more elevated immune system in the Fayoumis under homeostatic conditions. Overall, this study helped characterize the function of this important tissue and its response to NDV.

Deep Sequencing-Based Transcriptome Profiling Reveals Avian Interferon-Stimulated Genes and Provides Comprehensive Insight into Newcastle Disease Virus-Induced Host Responses

Newcastle disease virus (NDV) is an avian paramyxovirus that causes significant economic losses to the poultry industry worldwide, with variations in NDV pathogenicity due to the differences in virulence between strains. However, there is limited knowledge regarding the avian innate immune response to NDV infection. In this study, transcriptional profiles were obtained from chick embryo fibroblasts (CEFs) that were infected with the highly virulent NDV Herts/33 strain or the nonvirulent LaSota strain using RNA-seq. This yielded 8433 transcripts that were associated with NDV infection. This list of candidate genes was then further examined using Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses. It showed a high enrichment in the areas of cellular components and metabolic processes, with the cellular components possibly being associated with NDV pathogenicity. Among these 8433 transcripts, 3616 transcripts associated with interferon-stimulated genes (ISGs) were obtained; these transcripts are involved in metabolic processes, including protein phosphorylation and protein modification. These results provide further insight into the identification of genes that are involved in NDV infection. The global survey of changes in gene expression performed herein provides new insights into the complicated molecular mechanisms underlying virus and host interactions and will enable the use of new strategies to protect chickens against this virus.

Transcriptional Innate Immune Response of the Developing Chicken Embryo to Newcastle Disease Virus Infection

Traditional approaches to assess the immune response of chickens to infection are through animal trials, which are expensive, require enhanced biosecurity, compromise welfare, and are frequently influenced by confounding variables. Since the chicken embryo becomes immunocompetent prior to hatch, we here characterized the transcriptional response of selected innate immune genes to Newcastle disease virus (NDV) infection in chicken embryos at days 10, 14, and 18 of embryonic development. The results suggest that the innate immune response 72 h after challenge of 18-day chicken embryo is both consistent and robust. The expression of CCL5, Mx1, and TLR3 in lung tissues of NDV challenged chicken embryos from the outbred Kuroiler and Tanzanian local ecotype lines showed that their expression was several orders of magnitude higher in the Kuroiler than in the local ecotypes. Next, the expression patterns of three additional innate-immunity related genes, IL-8, IRF-1, and STAT1, were examined in the highly congenic Fayoumi (M5.1 and M15.2) and Leghorn (Ghs6 and Ghs13) sublines that differ only at the microchromosome bearing the major histocompatibility locus. The results show that the Ghs13 Leghorn subline had a consistently higher expression of all genes except IL-8 and expression seemed to be subline-dependent rather than breed-dependent, suggesting that the innate immune response of chicken embryos to NDV infection may be genetically controlled by the MHC-locus. Taken together, the results suggest that the chicken embryo may represent a promising model to studying the patterns and sources of variation of the avian innate immune response to infection with NDV and related pathogens.