BackgroundMarek’s disease (MD) is a commercially important neoplastic disease of chickens caused by the Marek’s disease virus (MDV), a naturally occurring oncogenic alphaherpesvirus. Enhancing MD genetic resistance is desirable to augment current vaccines and other MD control measures. High throughput sequencing was used to profile splenic transcriptomes from individual F1 progeny infected with MDV at 4 days of age from both outbred broilers (meat-type) and inbred layer (egg-type) chicken lines that differed in MD genetic resistance. The resulting information was used to identify SNPs, genes, and biological pathways exhibiting allele-specific expression (ASE) in response to MDV infection in each type of chicken. In addition, we compared and contrasted the results of pathway analyses (ASE and differential expression (DE)) between chicken types to help inform on the biological response to MDV infection.ResultsWith 7 individuals per line and treatment group providing high power, we identified 6,132 single nucleotide polymorphisms (SNPs) in 4,768 genes and 4,528 SNPs in 3,718 genes in broilers and layers, respectively, that exhibited ASE in response to MDV infection. Furthermore, 548 and 434 genes in broilers and layers, respectively, were found to show DE following MDV infection. Comparing the datasets, only 72 SNPs and 850 genes for ASE and 20 genes for DE were common between the two bird types. Although the chicken types used in this study were genetically different, at the pathway level, both TLR receptor and JAK/STAT signaling pathways were enriched as well as exhibiting a high proportion of ASE genes, especially at the beginning of both above mentioned regulatory pathways.ConclusionsRNA sequencing with adequate biological replicates is a powerful approach to identify high confidence SNPs, genes, and pathways that are associated with transcriptional response to MDV infection. In addition, the SNPs exhibiting ASE in response to MDV infection provide a strong foundation for determining the extent to which variation in expression influences MD incidence plus yield genetic markers for genomic selection. However, given the paucity of overlap among ASE SNP sets (broilers vs. layers), it is likely that separate screens need to be incorporated for each population. Finally, comparison of gene lists obtained between these two diverse chicken types indicate the TLR and JAK/STAT signaling are conserved when responding to MDV infection and may be altered by selection of genes exhibiting ASE found at the start of each pathway.

Lymphoid leukosis (LL)-like lymphoma is a low-incidence yet costly and poorly understood disease of domestic chickens. The observed unique characteristics of LL-like lymphomas are that the incidence of the disease is chicken line dependent; pathologically, it appeared to mimic avian leukosis but is free of exogenous ALV infection; inoculation of the nonpathogenic ALV-E or MDV-2 (SB-1) boosts the incidence of the disease; and inoculation of both the nonpathogenic ALV-E and SB-1 escalates it to much higher levels. This study was designed to test the impact of two new ALV-E isolates, recently derived from commercial broiler breeder flocks, in combination with the nonpathogenic SB-1 on LL-like lymphoma incidences in both an experimental egg layer line of chickens and a commercial broiler breeder line of chickens under a controlled condition. Data from this study provided an additional piece of experimental evidence on the potency of nonpathogenic ALV-E, MDV-2, and ALV-E plus MDV-2 in boosting the incidence of LL-like lymphomas in susceptible chickens. This study also generated the first piece of genomic evidence that suggests host transcriptomic variation plays an important role in modulating LL-like lymphoma formation.

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