A review of the development of chicken lines to resolve genes determining resistance to diseases

Epigenetic Factor MicroRNAs Likely Mediate Vaccine Protection Efficacy against Lymphomas in Response to Tumor Virus Infection in Chickens through Target Gene Involved Signaling Pathways

Epigenetic factors, including microRNAs (miRNAs), play an important role in affecting gene expression and, therefore, are involved in various biological processes including immunity protection against tumors. Marek's disease (MD) is a highly contagious disease of chickens caused by the MD virus (MDV). MD has been primarily controlled by vaccinations. MD vaccine efficacy might, in part, be dependent on modulations of a complex set of factors including host epigenetic factors. This study was designed to identify differentially expressed miRNAs in the primary lymphoid organ, bursae of Fabricius, in response to MD vaccination followed by MDV challenge in two genetically divergent inbred lines of White Leghorns. Small RNA sequencing and bioinformatic analyses of the small RNA sequence reads identified hundreds of miRNAs among all the treatment groups. A small portion of the identified miRNAs was differentially expressed within each of the four treatment groups, which were HVT or CVI988/Rispens vaccinated line 63-resistant birds and line 72-susceptible birds. A direct comparison between the resistant line 63 and susceptible line 72 groups vaccinated with HVT followed by MDV challenge identified five differentially expressed miRNAs. Gene Ontology analysis of the target genes of those five miRNAs revealed that those target genes, in addition to various GO terms, are involved in multiple signaling pathways including MAPK, TGF-β, ErbB, and EGFR1 signaling pathways. The general functions of those pathways reportedly play important roles in oncogenesis, anti-cancer immunity, cancer cell migration, and metastatic progression. Therefore, it is highly likely that those miRNAs may, in part, influence vaccine protection through the pathways.

Effects of DHAV-3 infection on innate immunity, antioxidant capacity, and lipid metabolism in ducks with different DHAV-3 susceptibilities

The aim of the experiment was to evaluate the status of innate immunity, oxidative status and lipid accumulation in ducklings exhibiting varying susceptibilities to DHAV-3 infection. In the experiment, ducklings with different DHAV-3 susceptibilities were used. Samples were collected at 6, 12, 15, and 24 h post infection (hpi), with 5 samples per time point. Plasma biochemistry, antioxidant enzyme activities, lipid content of liver and kidney were detected in the experiment. Elevated plasma level of total bilirubin, direct bilirubin, and creatinine indicated the injury of liver and kidney in susceptible ducklings (P < 0.05). The histopathological sections showed the injury in kidney. During the infection time, there was an increase in the concentrations of reactive oxygen species and oxidative damage markers (malondialdehyde and nitric oxide) in plasma of susceptible ducklings, particularly at 24 hpi (P < 0.05). Compared with the resistant ducklings, DHAV-3 infection resulted in a significant increase in the plasma total triglyceride (TG) level and a decrease in glucose level in susceptible ducklings. Gene expression of the innate immune response was both investigated in liver and kidney. In resistant ducklings, the expressions levels of pattern recognition receptors RIG-I, MDA5 remained constant. In contrast, the gene expressions peaked at 24 hpi in the susceptible ducklings. DHAV-3 infection promoted the expression of IFN, IL6, IL12β, caspase-8 or caspase-9 in both liver and kidney of susceptible ducklings. In conclusion, DHAV-3 infection led to the mobilization of antioxidant defenses, alterations in lipid metabolism, and oxidative stress in susceptible ducklings during DHAV-3 infection.

Immunoglobulin-like receptors in chickens: identification, functional characterization, and renaming to cluster homolog of immunoglobulin-like receptors

The cluster homolog of immunoglobulin-like receptors (CHIRs), previously known as the "chicken homolog of immunogloublin-like receptors," represents is a large group of transmembrane glycoproteins that direct the immune response. However, the full repertoire of putatively activating, inhibitory, or dual function CHIRA, CHIRB, and CHIRAB on chickens' immune responses is poorly understood. Herein, the study objective was to determine the genes encoding CHIR proteins and predict their function by searching canonical protein structure. A bioinformatics pipeline based on previous work was employed to search for the CHIRs from the newly updated broiler and layer genomes. The categorization into CHIRA, CHIRB, and CHIRAB types was assigned through motif searches, multiple sequence alignment, and phylogeny. In total, 150 protein-encoding genes on Chromosome 31 were identified as CHIRs. Gene members of each functional group (CHIRA, CHIRB, CHIRAB) were classified in accordance with previously recognized proteins. The genes were renamed to "cluster homolog of immunoglobulin-like receptors" (CHIRs) to allow for the naming of orthologous genes in other avian species. Additionally, expression analysis of the classified CHIRs across various reinforces their importance as immune regulators and activation in inflammatory tissues. Furthermore, over 1,000 diverse and rare CHIRs variants associated with differential Marek's disease response (P < 0.05) emphasize the impact of CHIRs on shaping avian immune responses in diverse contexts. The practical applications of these findings encompass advancing immunology, improving poultry health management, optimizing breeding programs for disease resistance, and enhancing overall animal health through a deeper understanding of the roles and functions of CHIRA, CHIRB, and CHIRAB types in avian immune responses.

Transgenerational epigenetic inheritance and immunity in chickens that vary in Marek's disease resistance

Marek's disease virus (MDV), a naturally oncogenic, highly contagious alpha herpesvirus, induces a T cell lymphoma in chickens that causes severe economic loss. Marek's disease (MD) outcome in an individual is attributed to genetic and environmental factors. Further investigation of the host-virus interaction mechanisms that impact MD resistance is needed to achieve greater MD control. This study analyzed genome-wide DNA methylation patterns in 2 highly inbred parental lines 63 and 72 and 5 recombinant congenic strains (RCS) C, L, M, N, and X strains from those parents. Lines 63 and 72, are MD resistant and susceptible, respectively, whereas the RCS have different combinations of 87.5% Line 63 and 12.5% Line 72. Our DNA methylation cluster showed a strong association with MD incidence. Differentially methylated regions (DMRs) between the parental lines and the 5 RCS were captured. MD-resistant and MD-susceptible markers of DNA methylation were identified as transgenerational epigenetic inheritable. In addition, the growth of v-src DNA tumors and antibody response against sheep red blood cells differed among the 2 parental lines and the RCS. Overall, our results provide very solid evidence that DNA methylation patterns are transgenerational epigenetic inheritance (TEI) in chickens and also play a vital role in MD tumorigenesis and other immune responses; the specific methylated regions may be important modulators of general immunity.

Epidemiology, pathogenesis, and diagnosis of Aleutian disease caused by Aleutian mink disease virus: A literature review with a perspective of genomic breeding for disease control in American mink (Neogale vison)

Aleutian disease (AD) is a multi-systemic infectious disease in American mink (Neogale vison) caused by the Aleutian mink disease virus (AMDV). Commonly referred to as mink plasmacytosis, AD is an economically significant disease in mink-breeding countries. Aleutian disease mainly induces weight loss, lower fertility, and dropped pelt quality in adults and can result in acute interstitial pneumonia with high mortality rates in kits. In this review, we employed the scientific literature on AD over the last 70 years to discuss the historical and contemporary status of AD outbreaks and seroprevalence in mink farming countries. We also explained different forms of AD and the differences between the pathogenicity of the virus in kits and adults. The application of the available AD serological tests in AD control strategies was argued. We explained how selection programs could help AD control and proposed different approaches to selecting animals for building AD-tolerant herds. The advantages of genomic selection for AD tolerance over traditional breeding strategies were discussed in detail. We also explained how genomic selection could help AD control by selecting tolerant animals for the next generation based on genome-wide single nucleotide polymorphisms (SNP) data and the challenges of implementing genomic selection for AD tolerance in the mink industry. This review collected the information required for designing successful breeding programs for AD tolerance. Examples of the application of information are presented, and data gaps are highlighted. We showed that AD tolerance is necessary to be among the traits that animals are selected for in the mink industry.

Avian Leukosis: Will We Be Able to Get Rid of It?

Avian leukosis viruses (ALVs) have been virtually eradicated from commercial poultry. However, some niches remain as pockets from which this group of viruses may reemerge and induce economic losses. Such is the case of fancy, hobby, backyard chickens and indigenous or native breeds, which are not as strictly inspected as commercial poultry and which have been found to harbor ALVs. In addition, the genome of both poultry and of several gamebird species contain endogenous retroviral sequences. Circumstances that support keeping up surveillance include the detection of several ALV natural recombinants between exogenous and endogenous ALV-related sequences which, combined with the well-known ability of retroviruses to mutate, facilitate the emergence of escape mutants. The subgroup most prevalent nowadays, ALV-J, has emerged as a multi-recombinant which uses a different receptor from the previously known subgroups, greatly increasing its cell tropism and pathogenicity and making it more transmissible. In this review we describe the ALVs, their different subgroups and which receptor they use to infect the cell, their routes of transmission and their presence in different bird collectivities, and the immune response against them. We analyze the different systems to control them, from vaccination to the progress made editing the bird genome to generate mutated ALV receptors or selecting certain haplotypes.

The immune cell landscape and response of Marek's disease resistant and susceptible chickens infected with Marek's disease virus

Genetically resistant or susceptible chickens to Marek's disease (MD) have been widely used models to identify the molecular determinants of these phenotypes. However, these prior studies lacked the basic identification and understanding of immune cell types that could be translated toward improved MD control. To gain insights into specific immune cell types and their responses to Marek's disease virus (MDV) infection, we used single-cell RNA sequencing (scRNAseq) on splenic cells from MD resistant and susceptible birds. In total, 14,378 cells formed clusters that identified various immune cell types. Lymphocytes, specifically T cell subtypes, were the most abundant with significant proportional changes in some subtypes upon infection. The largest number of differentially expressed genes (DEG) response was seen in granulocytes, while macrophage DEGs differed in directionality by subtype and line. Among the most DEG in almost all immune cell types were granzyme and granulysin, both associated with cell-perforating processes. Protein interactive network analyses revealed multiple overlapping canonical pathways within both lymphoid and myeloid cell lineages. This initial estimation of the chicken immune cell type landscape and its accompanying response will greatly aid efforts in identifying specific cell types and improving our knowledge of host response to viral infection.

Targeted Ablation of Exon 2 of the Avian Leukosis Virus-A (ALV-A) Receptor Gene in a Chicken Fibroblast Cell Line by CRISPR Abrogates ALV-A Infection

The U.S. Department of Agriculture Avian Disease and Oncology Laboratory currently relies on live birds of specific genetic backgrounds for producing chicken-embryo fibroblasts that are used for the diagnosis and subtyping of field isolates associated with avian leukosis virus (ALV) outbreaks. As an alternative to maintaining live animals for this purpose, we are currently developing cell lines capable of achieving the same result by ablation of the entry receptors utilized by ALV strains. We used CRISPR-Cas9 on the cell fibroblast-derived cell line DF-1 to disrupt the tva gene, which encodes the receptor required for binding and entry of ALV-A into cells. We ultimately identified seven DF-1 clones that had biallelic and homozygous indels at the Cas9 target site, exon 2 of tva. When tested in vitro for their ability to host ALV-A, the five clones that had frameshift mutations that disrupted the Tva protein were unable to support ALV-A replication. This result clearly demonstrates that modified cell lines can be used as part of a battery of tests to determine ALV subtype for isolate characterization, thus eliminating the need for live birds.

Small RNA deep sequencing revealed microRNAs' involvement in modulating cellular senescence and immortalization state

Unlike rodent cells, spontaneous immortalization of avian cells and human cells is a very rare event. According to patent publications and current literature, there are no more than 4 spontaneously immortalized chicken embryo fibroblast (CEF) cell lines established up to date. One of those cell lines is ADOL (Avian Disease and Oncology Laboratory) ZS-1 cell line, which was established by continuous passaging of the CEFs derived from the specific pathogen free (SPF) 0.TVB*S1 (commonly known as rapid feathering susceptible or RFS) genetic line of chickens. The RFS genetic line of chickens was developed and has been maintained on the SPF chicken farm of USDA-ARS facility, ADOL, in East Lansing, Michigan, which is known as one of a few lines of chickens that are free of any known avian endogenous virus genes. To explore potential roles that epigenetic factors may play in modulating cellular senescence processes and spontaneous immortalization state, total RNAs extracted from samples of the RFS primary CEFs, RFS CEFs reached the 21st passage, and the ZS-1 cells were subjected to small RNA sequencing. Collectively, a total of 531 miRNAs was identified in the 3 types of samples. In contrast to the primary CEF samples, 50 miRNAs were identified with significantly differential expression only in the 21st passage samples; a different subset of 63 differentially expressed miRNAs was identified only in the ZS-1 samples; the majority of differentially expressed miRNAs identified in both the 21st passage CEF and the ZS-1 samples were more or less directionally consistent. Gene Ontology analysis results suggested that the epigenetic factor, miRNAs, plays a role in modulating the cellular senescence and spontaneous immortalization processes through various bioprocesses and key pathways including ErbB and MAPK signaling pathways. These findings provided the experimental and bioinformatic evidence for a better understanding on the epigenetic factor of miRNAs in association with cellular senescence and spontaneous immortalization process in avian cells.

Short communication: diversity of endogenous avian leukosis virus subgroup E elements in 11 chicken breeds

Avian leukosis virus subgroup E (ALVE) as a kind of endogenous retroviruses extensively exists in chicken genome. The insertion of ALVE has some effects on chicken production traits and appearance. Most of the work on ALVEs has been done with commercial breeds. We present here an investigation of ALVE elements in seven Chinese domestic breeds and four standard breeds. Firstly, we established an ALVE insertion site dataset by using the obsERVer pipeline to identify ALVEs from whole-genome sequence data of eleven chicken breeds, seven Chinese domestic breeds, including Beijing You (BY), Dongxiang (DX), Luxi Game (LX), Shouguang (SG), Silkie (SK), Tibetan (TB) and Wenchang (WC), four standard breeds, including White Leghorn (WL), White Plymouth Rock (WR), Cornish (CS), and Rhode Island Red (RIR). A total of 37 ALVE insertion sites were identified and 23 of them were novel. Most of these insertion sites were distributed in intergenic regions and introns. We then used locus-specific PCR to validate the insertion sites in an expanded population with 18~60 individuals in each breed. The results showed that all predicted integration sites in 11 breeds were verified by PCR. Some ALVE insertion sites were breeds specific, and 16 out of 23 novel ALVEs were found in only one Chinese domestic chicken breed. We randomly selected three ALVE insertions including ALVE_CAU005, ALVE_ros127, and ALVE_ros276, and obtained their insertion sequences by long-range PCR and Sanger sequencing. The insertion sequences were all 7525 bp, which were full-length ALVE insertion and all of them were highly homologous to ALVE1 with similarity of 99%. Our study identified the distribution of ALVE in 11 chicken breeds, which expands the current research on ALVE in Chinese domestic breeds.

The Role of Dendritic Cells in the Host Response to Marek’s Disease Virus (MDV) as Shown by Transcriptomic Analysis of Susceptible and Resistant Birds

Despite the successful control of highly contagious tumorigenic Marek’s disease (MD) by vaccination, a continuous increase in MD virus (MDV) virulence over recent decades has put emphasis on the development of more MD-resistant chickens. The cell types and genes involved in resistance therefore need to be recognized. The virus is primarily lymphotropic, but research should also focus on innate immunity, as innate immune cells are among the first to encounter MDV. Our previous study on MDV–macrophage interaction revealed significant differences between MHC-congenic lines 61 (MD-resistant) and 72 (MD-susceptible). To investigate the role of dendritic cells (DCs) in MD resistance, bone-marrow-derived DCs from these lines were infected with MDV in vitro. They were then characterized by cell sorting, and the respective transcriptomes analysed by RNA-seq. The differential expression (DE) of genes revealed a strong immune activation in DCs of the susceptible line, although an inherent immune supremacy was shown by the resistant line, including a significant expression of tumour-suppressor miRNA, gga-mir-124a, in line 61 control birds. Enrichment analysis of DE genes revealed high expression of an oncogenic transcription factor, AP-1, in the susceptible line following MDV challenge. This research highlights genes and pathways that may play a role in DCs in determining resistance or susceptibility to MDV infection.

Genetic variation in chicken interferon signalling pathway genes in research lines showing differential viral resistance

Avian viruses of economic interest are a significant burden on the poultry industry, affecting production traits and resulting in mortality. Furthermore, the zoonosis of avian viruses risks pandemics developing in humans. Vaccination is the most common method of controlling viruses; however current vaccines often lack cross-protection against multiple strains of each virus. The mutagenicity of these viruses has also led to virulent strains emerging that can overcome the protection offered by vaccines. Breeding chickens with a more robust innate immune response may help in tackling current and emerging viruses. Understanding the genetic evolution of different lines will thus provide a useful tool in helping the host in the fight against pathogens. This study focuses on the interferon genes and their receptors in different chicken lines that are known to be more resistant or susceptible to particular avian viruses. Comparing genotypic differences in these core immune genes between the chicken lines may explain the phenotypic differences observed and aid the identification of causative variations. The relative resistance/susceptibility of each line to viruses of interest (Marek's disease virus, infectious bursal disease, infectious bronchitis virus and avian influenza virus) has previously been determined. Here we identify single nucleotide polymorphisms in interferons and downstream genes. Functional prediction tools were used to identify variants that may be affecting protein structure, mRNA secondary structure or transcription factor and micro-RNA binding sites. These variants were then considered in the context of the research lines and their distribution between phenotypes. We highlight 60 variants of interest in the interferon pathway genes that may account for susceptibility/resistance to viral pathogens.

MDV-induced differential microRNA expression in the primary lymphoid organ of thymus

Marek's disease virus (MDV), a highly contagious cell associated virus, is the etiological agent of Marek's disease (MD), a lymphoproliferative and neuropathic disease of domestic chickens. Clinical signs of MD include transient paralysis, bursal/thymic atrophy, and T cell lymphomas. MicroRNAs (miRNAs) are short single-stranded non-coding RNAs that regulate gene expression by transcriptional suppression or mRNA degradation. Herpesviruses, including MDV, encode for miRNAs that are known to play essential roles in viral pathogenicity, oncogenesis, and evasion of immune responses. In this study, we performed miRNA sequencing in thymuses of control and MDV-infected chickens of MD-resistant (6₃) and susceptible (7₂) lines at 21 days post infection (dpi). The thymus is a lymphoid organ that undergoes severe atrophy due to MDV-induced apoptotic mediated destruction of T cells. Sequence analysis identified 658 total chicken miRNAs in the thymuses of control and MDV-infected birds of both lines. Of these, 453 were novel and 205 were known microRNAs. All novel miRNAs mapped to chicken genome with no sequence homology to existing miRNAs in the chicken miRbase. Comparative analysis between the thymuses of control and infected birds of resistant and susceptible lines identified 78 differentially expressed microRNAs that might provide insights into mechanisms of thymus atrophy.

The Microbiota-Gut-Brain Axis: Gut Microbiota Modulates Conspecific Aggression in Diversely Selected Laying Hens

The gut microbiota plays an important role in regulating brain function, influencing psychological and emotional stability. The correlations between conspecific aggression, gut microbiota, and physiological homeostasis were further studied in inbred laying chicken lines, 63 and 72, which were diversely selected for Marek’s disease, and they also behave differently in aggression. Ten sixty-week-old hens from each line were sampled for blood, brain, and cecal content. Neurotransmitters, cytokines, corticosterone, and heterophil/lymphocyte ratios were determined. Cecal microbiota compositions were determined by bacterial 16s rRNA sequencing, and functional predictions were performed. Our data showed that the central serotonin and tryptophan levels were higher in line 63 compared to line 72 (p < 0.05). Plasma corticosterone, heterophil/lymphocyte ratios, and central norepinephrine were lower in line 63 (p < 0.05). The level of tumor necrosis factor α tended to be higher in line 63. Faecalibacterium, Oscillibacter, Butyricicoccus, and Bacteriodes were enriched in line 63 birds, while Clostridiales vadin BB60, Alistipes, Mollicutes RF39 were dominated in line 72. From the predicted bacterial functional genes, the kynurenine pathway was upregulated in line 72. These results suggested a functional linkage of the line differences in serotonergic activity, stress response, innate immunity, and gut microbiota populations.

Marek's Disease Virus Telomeric Integration Profiles of Neoplastic Host Tissues Reveal Unbiased Chromosomal Selection and Loss of Cellular Diversity during Tumorigenesis

The avian α-herpesvirus known as Marek's disease virus (MDV) linearly integrates its genomic DNA into host telomeres during infection. The resulting disease, Marek's disease (MD), is characterized by virally-induced lymphomas with high mortality. The temporal dynamics of MDV-positive (MDV+) transformed cells and expansion of MD lymphomas remain targets for further understanding. It also remains to be determined whether specific host chromosomal sites of MDV telomere integration confer an advantage to MDV-transformed cells during tumorigenesis. We applied MDV-specific fluorescence in situ hybridization (MDV FISH) to investigate virus-host cytogenomic interactions within and among a total of 37 gonad lymphomas and neoplastic splenic samples in birds infected with virulent MDV. We also determined single-cell, chromosome-specific MDV integration profiles within and among transformed tissue samples, including multiple samples from the same bird. Most mitotically-dividing cells within neoplastic samples had the cytogenomic phenotype of 'MDV telomere-integrated only', and tissue-specific, temporal changes in phenotype frequencies were detected. Transformed cell populations composing gonad lymphomas exhibited significantly lower diversity, in terms of heterogeneity of MDV integration profiles, at the latest stages of tumorigenesis (>50 days post-infection (dpi)). We further report high interindividual and lower intraindividual variation in MDV integration profiles of lymphoma cells. There was no evidence of integration hotspots into a specific host chromosome(s). Collectively, our data suggests that very few transformed MDV+ T cell populations present earlier in MDV-induced lymphomas (32-50 dpi), survive, and expand to become the dominant clonal population in more advanced MD lymphomas (51-62 dpi) and establish metastatic lymphomas.

Genetic Diversity of MHC B-F/B-L Region in 21 Chicken Populations

The chicken major histocompatibility complex (MHC) on chromosome 16 is the most polymorphic region across the whole genome, and also an ideal model for genetic diversity investigation. The MHC B-F/B-L region is 92 kb in length with high GC content consisting of 18 genes and one pseudogene (Blec4), which plays important roles in immune response. To evaluate polymorphism of the Chinese indigenous chickens as well as to analyze the effect of selection to genetic diversity, we used WaferGen platform to identify sequence variants of the B-F/B-L region in 21 chicken populations, including the Red Jungle Fowl (RJF), Cornish (CS), White Leghorns (WLs), 16 Chinese domestic breeds, and two well-known inbred lines 63 and 72. A total of 3,319 single nucleotide polymorphism (SNPs) and 181 INDELs in the B-F/B-L region were identified among 21 populations, of which 2,057 SNPs (62%) and 159 INDELs (88%) were novel. Most of the variants were within the intron and the flanking regions. The average variation density was 36 SNPs and 2 INDELs per kb, indicating dramatical high diversity of this region. Furthermore, BF2 was identified as the hypervariable genes with 67 SNPs per kb. Chinese domestic populations showed higher diversity than the WLs and CS. The indigenous breeds, Nandan Yao (NY), Xishuangbanna Game (XG), Gushi (GS), and Xiayan (XY) chickens, were the top four with the highest density of SNPs and INDELs. The highly inbred lines 63 and 72 have the lowest diversity, which might be resulted from a long-term intense selection for decades. Collectively, we refined the genetic map of chicken MHC B-F/B-L region, and illustrated genetic diversity of 21 chicken populations. Abundant genetic variants were identified, which not only strikingly expanded the current Ensembl SNP database, but also provided comprehensive data for researchers to further investigate association between variants in MHC and immune traits.

Chickens as a simple system for scientific discovery: The example of the MHC

Chickens have played many roles in human societies over thousands of years, most recently as an important model species for scientific discovery, particularly for embryology, virology and immunology. In the last few decades, biomedical models like mice have become the most important model organism for understanding the mechanisms of disease, but for the study of outbred populations, they have many limitations. Research on humans directly addresses many questions about disease, but frank experiments into mechanisms are limited by practicality and ethics. For research into all levels of disease simultaneously, chickens combine many of the advantages of humans and of mice, and could provide an independent, integrated and overarching system to validate and/or challenge the dogmas that have arisen from current biomedical research. Moreover, some important systems are simpler in chickens than in typical mammals. An example is the major histocompatibility complex (MHC) that encodes the classical MHC molecules, which play crucial roles in the innate and adaptive immune systems. Compared to the large and complex MHCs of typical mammals, the chicken MHC is compact and simple, with single dominantly-expressed MHC molecules that can determine the response to infectious pathogens. As a result, some fundamental principles have been easier to discover in chickens, with the importance of generalist and specialist MHC alleles being the latest example.

Transcriptomes of an Array of Chicken Ovary, Intestinal, and Immune Cells and Tissues

While the chicken (Gallus gallus) is the most consumed agricultural animal worldwide, the chicken transcriptome remains understudied. We have characterized the transcriptome of 10 cell and tissue types from the chicken using RNA-seq, spanning intestinal tissues (ileum, jejunum, proximal cecum), immune cells (B cells, bursa, macrophages, monocytes, spleen T cells, thymus), and reproductive tissue (ovary). We detected 17,872 genes and 24,812 transcripts across all cell and tissue types, representing 73% and 63% of the current gene annotation, respectively. Further quantification of RNA transcript biotypes revealed protein-coding and lncRNAs specific to an individual cell/tissue type. Each cell/tissue type also has an average of around 1.2 isoforms per gene, however, they all have at least one gene with at least 11 isoforms. Differential expression analysis revealed a large number of differentially expressed genes between tissues of the same category (immune and intestinal). Many of these differentially expressed genes in immune cells were involved in cellular processes relating to differentiation and cell metabolism as well as basic functions of immune cells such as cell adhesion and signal transduction. The differential expressed genes of the different segments of the chicken intestine (jejunum, ileum, proximal cecum) correlated to the metabolic processes in nutrient digestion and absorption. These data should provide a valuable resource in understanding the chicken genome.

Transcriptomic Analysis of Inbred Chicken Lines Reveals Infectious Bursal Disease Severity Is Associated with Greater Bursal Inflammation In Vivo and More Rapid Induction of Pro-Inflammatory Responses in Primary Bursal Cells Stimulated Ex Vivo

In order to better understand differences in the outcome of infectious bursal disease virus (IBDV) infection, we inoculated a very virulent (vv) strain into White Leghorn chickens of inbred line W that was previously reported to experience over 24% flock mortality, and three inbred lines (15I, C.B4 and 0) that were previously reported to display no mortality. Within each experimental group, some individuals experienced more severe disease than others but line 15I birds experienced milder disease based on average clinical scores, percentage of birds with gross pathology, average bursal lesion scores and average peak bursal virus titre. RNA-Seq analysis revealed that more severe disease in line W was associated with significant up-regulation of pathways involved in inflammation, cytoskeletal regulation by Rho GTPases, nicotinic acetylcholine receptor signaling, and Wnt signaling in the bursa compared to line 15I. Primary bursal cell populations isolated from uninfected line W birds contained a significantly greater percentage of KUL01+ macrophages than cells isolated from line 15I birds (p < 0.01) and, when stimulated ex vivo with LPS, showed more rapid up-regulation of pro-inflammatory gene expression than those from line 15I birds. We hypothesize that a more rapid induction of pro-inflammatory cytokine responses in bursal cells following IBDV infection leads to more severe disease in line W birds than in line 15I.

Comparison of Marek's disease virus challenge strains and bird types for vaccine licensing

Marek's disease virus (MDV) is an important poultry pathogen which is controlled through widespread vaccination with avirulent and attenuated strains, but continued evolution of field viruses to higher virulence has required ongoing improvement of available vaccine strains, and these vaccine strains also offer an attractive platform for designing recombinant vector vaccines with cross-protection against MDV and additional pathogens.  Recent reports of failures in vaccine licensing trials of positive controls to reach appropriately high levels of MD incidence prompted us to evaluate possible combinations of outbred specific pathogen-free (SPF) layer lines and alternative virulent challenge strains which could provide more consistent models for serotype-3 vectored vaccine development.  Choice of layer line and virulent MDV challenge strain each contributed to the ability of a challenge model to reach 80 percent virulence in unvaccinated positive control groups in the majority of trials without overwhelming serotype-3 vectored vaccine protection in vaccinated groups.  Conversely, reducing challenge virus dose by a factor of four, or vaccine dose by half, had no consistent effect across these models.  Although MDV strain 617A had the most potential as an alternative to strains that are currently approved for licensing trials, no combination of layer line and challenge virus consistently met the goals for a successful challenge model in all study replicates, indicating that high variability is an inherent difficulty in MDV challenge studies, at least when outbred birds are used.

The Chicken MHC: Insights into Genetic Resistance, Immunity, and Inflammation Following Infectious Bronchitis Virus Infections

The chicken immune system has provided an immense contribution to basic immunology knowledge by establishing major landmarks and discoveries that defined concepts widely used today. One of many special features on chickens is the presence of a compact and simple major histocompatibility complex (MHC). Despite its simplicity, the chicken MHC maintains the essential counterpart genes of the mammalian MHC, allowing for a strong association to be detected between the MHC and resistance or susceptibility to infectious diseases. This association has been widely studied for several poultry infectious diseases, including infectious bronchitis. In addition to the MHC and its linked genes, other non-MHC loci may play a role in the mechanisms underlying such resistance. It has been reported that innate immune responses, such as macrophage function and inflammation, might be some of the factors driving resistance or susceptibility, consequently influencing the disease outcome in an individual or a population. Information about innate immunity and genetic resistance can be helpful in developing effective preventative measures for diseases such as infectious bronchitis, to which a systemic antibody response is often not associated with disease protection. In this review, we summarize the importance of the chicken MHC in poultry disease resistance, particularly to infectious bronchitis virus (IBV) infections and the role played by innate immunity and inflammation on disease outcome. We highlight how future studies focusing on the MHC and non-MHC genes can potentially bring clarity to observed resistance in some chicken B haplotype lines.

Genome-wide characterization of copy number variations in the host genome in genetic resistance to Marek's disease using next generation sequencing

Background

Marek’s disease (MD) is a highly neoplastic disease primarily affecting chickens, and remains as a chronic infectious disease that threatens the poultry industry. Copy number variation (CNV) has been examined in many species and is recognized as a major source of genetic variation that directly contributes to phenotypic variation such as resistance to infectious diseases. Two highly inbred chicken lines, 6₃ (MD-resistant) and 7₂ (MD-susceptible), as well as their F₁ generation and six recombinant congenic strains (RCSs) with varied susceptibility to MD, are considered as ideal models to identify the complex mechanisms of genetic and molecular resistance to MD.

Results

In the present study, to unravel the potential genetic mechanisms underlying resistance to MD, we performed a genome-wide CNV detection using next generation sequencing on the inbred chicken lines with the assistance of CNVnator. As a result, a total of 1649 CNV regions (CNVRs) were successfully identified after merging all the nine datasets, of which 90 CNVRs were overlapped across all the chicken lines. Within these shared regions, 1360 harbored genes were identified. In addition, 55 and 44 CNVRs with 62 and 57 harbored genes were specifically identified in line 6₃ and 7₂, respectively. Bioinformatics analysis showed that the nearby genes were significantly enriched in 36 GO terms and 6 KEGG pathways including JAK/STAT signaling pathway. Ten CNVRs (nine deletions and one duplication) involved in 10 disease-related genes were selected for validation by using quantitative real-time PCR (qPCR), all of which were successfully confirmed. Finally, qPCR was also used to validate two deletion events in line 7₂ that were definitely normal in line 6₃. One high-confidence gene, IRF2 was identified as the most promising candidate gene underlying resistance and susceptibility to MD in view of its function and overlaps with data from previous study.

Conclusions

Our findings provide valuable insights for understanding the genetic mechanism of resistance to MD and the identified gene and pathway could be considered as the subject of further functional characterization.

Association of LEI0258 Marker Alleles and Susceptibility to Virulent Newcastle Disease Virus Infection in Kuroiler, Sasso, and Local Tanzanian Chicken Embryos

Newcastle disease (ND) control by vaccination and an institution of biosecurity measures is less feasible in backyard chicken in developing countries. Therefore, an alternative disease control strategy like the genetic selection of less susceptible chicken genotypes is a promising option. In the present study, genetic polymorphism of LEIO258 marker and association with susceptibility to virulent Newcastle disease virus (NDV) infection in Kuroilers, Sasso, and local Tanzanian chicken embryos were investigated. Samples from high (15%) and less (15%) susceptible cohorts were genotyped by sequencing of LEI0258 marker. A total of 75 DNA sequences comprised of 29 Kuroiler, 29 local Tanzanian chickens, and 17 Sasso were analyzed. Neighbor-joining phylogenetic trees were constructed to depict the clustering of LEI0258 marker alleles and relationship with susceptibility. Alleles with frequency ≥3 were considered for association with susceptibility by the use of the inference technique. The present findings suggest that some LEI0258 marker genetic polymorphisms apart from LEI0258 marker allelic based on sizes may be linked with chicken MHC-B haplotypes that confer chickens variability in resistance or susceptibility to infections. Furthermore, these results demonstrate the presence of relationship between LEI0258 marker polymorphisms and variations in chicken susceptibility to NDV infection, which could be utilized in breeding programs designed to improve chicken disease resistance.

Evaluation of Protective Efficacy When Combining Turkey Herpesvirus-Vector Vaccines

Turkey herpesvirus (HVT) is widely used as a vaccine against Marek's disease in chickens and recently as a vector for foreign genes from infectious bursal disease virus, Newcastle disease (ND) virus, infectious laryngotracheitis (ILT) virus, and avian influenza virus. Advantages of HVT-vector vaccines are that the vaccines do not contain live respiratory viruses or live infectious bursal disease virus able to replicate and cause disease or embryo mortality, they can be administered at hatch or in ovo, and they are relatively insensitive to interference from maternally derived antibodies. As producers have tried to combine HVT-vector vaccines to protect against additional diseases, reports have indicated that applying two vectored vaccines using the same HVT vector is reported to reduce the efficacy of one or both vaccines. To confirm this interference, we evaluated commercial vaccines from multiple companies, including products with inserts designed to protect against ND, infectious ILT, and infectious bursal disease (IBD). Using a standard dosage, we found that the ILT product was most severely affected by the addition of other vaccines, as demonstrated by a significant increase in clinical signs, significant decrease in weight gain, and increase in quantity of challenge virus observed from tracheal swabs collected from Days 3-5 postchallenge. The ND and IBD products were also affected by the addition of other vaccines, although in most cases differences compared to vaccination with the vector alone were not statistically significant. This study demonstrates the importance of following manufacturer guidelines and the need for validating alternative strategies to benefit from the high level of protection offered by vector vaccines.

Pathogen transmission from vaccinated hosts can cause dose-dependent reduction in virulence

Many livestock and human vaccines are leaky because they block symptoms but do not prevent infection or onward transmission. This leakiness is concerning because it increases vaccination coverage required to prevent disease spread and can promote evolution of increased pathogen virulence. Despite leakiness, vaccination may reduce pathogen load, affecting disease transmission dynamics. However, the impacts on post-transmission disease development and infectiousness in contact individuals are unknown. Here, we use transmission experiments involving Marek disease virus (MDV) in chickens to show that vaccination with a leaky vaccine substantially reduces viral load in both vaccinated individuals and unvaccinated contact individuals they infect. Consequently, contact birds are less likely to develop disease symptoms or die, show less severe symptoms, and shed less infectious virus themselves, when infected by vaccinated birds. These results highlight that even partial vaccination with a leaky vaccine can have unforeseen positive consequences in controlling the spread and symptoms of disease.

Marek's disease vaccines-induced differential expression of known and novel microRNAs in primary lymphoid organ bursae of White Leghorn

Marek’s disease (MD) is a contagious disease of domestic chickens caused by MD viruses. MD has been controlled primarily by vaccinations, yet sporadic outbreaks of MD take place worldwide. Commonly used MD vaccines include HVT, SB-1 and CVI988/Rispens and their efficacies are reportedly dependent of multiple factors including host genetics. Our previous studies showed protective efficacy of a MD vaccine can differ drastically from one chicken line to the next. Advanced understanding on the underlying genetic and epigenetic factors that modulate vaccine efficacy would greatly improve the strategy in design and development of more potent vaccines. Two highly inbred lines of White Leghorn were inoculated with HVT and CVI988/Rispens. Bursa samples were taken 26 days post-vaccination and subjected to small RNA sequencing analysis to profile microRNAs (miRNA). A total of 589 and 519 miRNAs was identified in one line, known as line 6₃, 490 and 630 miRNAs were identified in the other, known as line 7₂, in response to HVT or CVI988/Rispens inoculation, respectively. HVT and CVI988/Rispens induced mutually exclusive 4 and 13 differentially expressed (DE) miRNAs in line 6₃ birds in contrast to a non-vaccinated group of the same line. HVT failed to induce any DE miRNA and CVI988/Rispens induced a single DE miRNA in line 7₂ birds. Thousands of target genes for the DE miRNAs were predicted, which were enriched in a variety of gene ontology terms and pathways. This finding suggests the epigenetic factor, microRNA, is highly likely involved in modulating vaccine protective efficacy in chicken.

MicroRNA profiling in the bursae of Marek's disease virus-infected resistant and susceptible chicken lines

Marek's disease (MD) is a lymphoproliferative disease of domestic chickens caused by a cell-associated oncogenic alpha-herpesvirus, Marek's disease virus (MDV). Clinical signs of MD include bursal/thymic atrophy, neurologic disorders, and T cell lymphomas. MiRNAs play key roles in regulation of gene expression by targeting translational suppression or mRNA degradation. MDV encodes miRNAs that are associated with viral pathogenicity and oncogenesis. In this study, we performed miRNA sequencing in the bursal tissues, non-tumorous but viral-induced atrophied lymphoid organ, from control and infected MD-resistant and susceptible chickens at 21 days post infection. In addition to some known miRNAs, a minimum of 300 novel miRNAs were identified in each group that mapped to the chicken genome with no sequence homology to existing miRNAs in chicken miRbase. Comparative analysis identified 54 deferentially expressed miRNAs between the chicken lines that might shed light on underlying mechanism of bursal atrophy and resistance or susceptibility to MD.

MicroRNA gga-miR-200a-3p modulates immune response via MAPK signaling pathway in chicken afflicted with necrotic enteritis

MicroRNAs (miRNAs) are small non-coding RNAs that contribute to host immune response as post-transcriptional regulation. The current study investigated the biological role of the chicken (Gallus gallus) microRNA-200a-3p (gga-miR-200a-3p), using 2 necrotic enteritis (NE) afflicted genetically disparate chicken lines, 6.3 and 7.2, as well as the mechanisms underlying the fundamental signaling pathways in chicken. The expression of gga-miR-200a-3p in the intestinal mucosal layer of NE-induced chickens, was found to be upregulated during NE infection in the disease-susceptible chicken line 7.2. To validate the target genes, we performed an overexpression analysis of gga-miR-200a-3p using chemically synthesized oligonucleotides identical to gga-miR-200a-3p, reporter gene analysis including luciferase reporter assay, and a dual fluorescence reporter assay in cultured HD11 chicken macrophage cell lines. Gga-miR-200a-3p was observed to be a direct transcriptional repressor of ZAK, MAP2K4, and TGFβ2 that are involved in mitogen-activated protein kinase (MAPK) pathway by targeting the 3′-UTR of their transcripts. Besides, gga-miR-200a-3p may indirectly affect the expression of protein kinases including p38 and ERK1/2 at both transcriptional and translational levels, suggesting that this miRNA may function as an important regulator of the MAPK signaling pathway. Proinflammatory cytokines consisting of IL-1β, IFN-γ, IL-12p40, IL-17A, and LITAF belonging to Th1 and Th17-type cytokines, were upregulated upon gga-miR-200a-3p overexpression. These findings have enhanced our knowledge of the immune function of gga-miR-200a-3p mediating the chicken immune response via regulation of the MAPK signaling pathway and indicate that this miRNA may serve as an important biomarker of diseases in domestic animals.

Differential expression of immune-related genes in the bursa of Fabricius of two inbred chicken lines following infection with very virulent infectious bursal disease virus

Among different inbred chickens' lines, we previously showed that lines P and N of Institute for Animal Health, Compton, UK are the most susceptible and the least affected lines, respectively, following infection with very virulent infectious bursal disease virus (vvIBDV). In this study, the differential expressions of 29 different immune-related genes were characterized. Although, birds from both lines succumbed to infection, line P showed greater bursal lesion scores and higher viral copy numbers compared to line N. Interestingly, line N showed greater down-regulation of B cell related genes (BLNK, TNFSF13B and CD72) compared to line P. While up-regulation of T-cell related genes (CD86 and CTLA4) and Th1 associated cytokines (IFNG, IL2, IL12A and IL15) were documented in both lines, the expression levels of these genes were different in the two lines. Meanwhile, the expression of IFN-related genes IFNB, STAT1, and IRF10, but not IRF5, were up-regulated in both lines. The expression of pro-inflammatory cytokines (IL1B, IL6, IL18, and IL17) and chemokines (CXCLi2, CCL4, CCL5 and CCR5) were up-regulated in both lines with greater increase documented in line P compared to line N. Strikingly, the expression of IL12B was detected only in line P whilst the expression of IL15RA was detected only in line N. In conclusion, the bursal immunopathology of IBDV correlates more with expression of proinflammatory response related genes and does not related to expression of B-cell related genes.

Parent-of-origin has no detectable effect on survival days of Marek's disease virus infected White Leghorns

Marek's Disease (MD) is a neoplastic disease of chickens and remains as a chronic infectious disease that threatens the poultry industry. Improving genetic resistance to MD in poultry is an important long-term goal, which would significantly augment the current control measures against MD and eventually reduce the annual economic loss. In this study, survival patterns of F2 birds from 2 reciprocal crosses were compared to examine possible difference in survival between the reciprocal crosses in response to MD virus (MDV) challenge. A total of 246 and 224 F2 chicks derived from reciprocal crosses of lines 63 × 72 and lines 72 × 63, respectively, were sampled from an MDV challenge trial and survival days were recorded from the MDV-inoculation date to the end of experiment. Statistical analyses, including Principal Component Analysis (PCA) followed by a cox-regression model, showed there was no significant difference in survival days between reciprocal crosses (P > 0.05). To the best of our knowledge, this is the first MD survival study on reciprocal crosses of 2 genetically diversified lines of chickens differing in MD resistance. This report documented the experimental evidence that the genetic lineage of grandparental (maternal or paternal) effect on survival days was minimal, if present at all.

Linc-GALMD1 Regulates Viral Gene Expression in the Chicken

A rapidly increasing number of reports on dysregulated long intergenic non-coding RNA (lincRNA) expression across numerous types of cancers indicates that aberrant lincRNA expression may be a major contributor to tumorigenesis. Marek’s disease (MD) is a T cell lymphoma of chickens induced by Marek’s disease virus (MDV). Although we have investigated the roles of lincRNAs in bursa tissue of MDV-infected chickens in previous studies, the molecular mechanisms of lincRNA functions in T cells remain poorly understood. In the present study, Linc-GALMD1 was identified from CD4+ T cells and MSB1 cells, and its expression was significantly downregulated in MD-resistant line of birds in response to MDV challenge. Furthermore, loss-of-function experiments indicated that linc-GALMD1 significantly affected the expression of 290 genes in trans. Through integrated analysis of differentially expressed genes (DEGs) induced by MDV and linc-GALMD1, we found that IGLL1 gene expression levels had a positive correlation with the degree of MD infection and could potentially serve as an indicator for clinical diagnosis of MD. Moreover, an interaction between MDV and linc-GALMD1 was also observed. Accordingly, chicken embryonic fibroblast cells were inoculated with MDV with and without the linc-GALMD1 knockdown, and the data showed that linc-GALMD1 could repress MDV gene expression during the course of MDV infection. These findings uncovered a role of linc-GALMD1 as a viral gene regulator and suggested a function of linc-GALMD1 contributing to tumor suppression by coordinating expression of MDV genes and tumor-related genes and regulating immune responses to MDV infection.

Endogenous Avian Leukosis Virus in Combination with Serotype 2 Marek's Disease Virus Significantly Boosted the Incidence of Lymphoid Leukosis-Like Bursal Lymphomas in Susceptible Chickens

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.

In 2010, sporadic cases of avian leukosis virus (ALV)-like bursal lymphoma, also known as spontaneous lymphoid leukosis (LL)-like tumors, were identified in two commercial broiler breeder flocks in the absence of exogenous ALV infection. Two individual ALV subgroup E (ALV-E) field strains, designated AF227 and AF229, were isolated from two different breeder farms. The role of these ALV-E field isolates in development of and the potential joint impact in conjunction with a Marek’s disease virus (MDV) vaccine (SB-1) were further characterized in chickens of an experimental line and commercial broiler breeders. The experimental line 0.TVB*S1, commonly known as the rapid feathering-susceptible (RFS) line, of chickens lacks all endogenous ALV and is fully susceptible to all subgroups of ALV, including ALV-E. Spontaneous LL-like tumors occurred following infection with AF227, AF229, and a reference ALV-E strain, RAV60, in RFS chickens. Vaccination with serotype 2 MDV, SB-1, in addition to AF227 or AF229 inoculation, significantly enhanced the spontaneous LL-like tumor incidence in the RFS chickens. The spontaneous LL-like tumor incidence jumped from 14% by AF227 alone to 42 to 43% by AF227 in combination with SB-1 in the RFS chickens under controlled conditions. RNA-sequencing analysis of the LL-like lymphomas and nonmalignant bursa tissues of the RFS line of birds identified hundreds of differentially expressed genes that are reportedly involved in key biological processes and pathways, including signaling and signal transduction pathways. The data from this study suggested that both ALV-E and MDV-2 play an important role in enhancement of the spontaneous LL-like tumors in susceptible chickens. The underlying mechanism may be complex and involved in many chicken genes and pathways, including signal transduction pathways and immune system processes, in addition to reported viral genes.

IMPORTANCE 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.

RNA Sequencing revealed differentially expressed genes functionally associated with immunity and tumor suppression during latent phase infection of a vv + MDV in chickens

Very virulent plus Marek's disease (MD) virus (vv + MDV) induces tumors in relatively resistant lines of chickens and early mortality in highly susceptible lines of chickens. The vv + MDV also triggers a series of cellular responses in both types of chickens. We challenged birds sampled from a highly inbred chicken line (line 6₃) that is relatively resistant to MD and from another inbred line (line 7₂) that is highly susceptible to MD with a vv + MDV. RNA-sequencing analysis was performed with samples extracted from spleen tissues taken at 10-day and 21-day post infection (dpi). A total of 64 and 106 differentially expressed genes was identified in response to the vv + MDV challenge at latent phase in the resistant and susceptible lines of chickens, respectively. Direct comparisons between samples of the two lines identified 90 and 126 differentially expressed genes for control and MDV challenged groups, respectively. The differentially expressed gene profiles illustrated that intensive defense responses were significantly induced by vv + MDV at 10 dpi and 21 dpi but with slight changes in the resistant line. In contrast, vv + MDV induced a measurable suppression of gene expression associated with host defense at 10 dpi but followed by an apparent activation of the defense response at 21 dpi in the susceptible line of chickens. The observed difference in gene expression between the two genetic lines of chickens in response to MDV challenge during the latent phase provided a piece of indirect evidence that time points for MDV reactivation differ between the genetic lines of chickens with different levels of genetic resistance to MD. Early MDV reactivation might be necessary and potent to host defense system readiness for damage control of tumorigenesis and disease progression, which consequently results in measurable differences in phenotypic characteristics including early mortality (8 to 20 dpi) and tumor incidence between the resistant and susceptible lines of chickens. Combining differential gene expression patterns with reported GO function terms and quantitative trait loci, a total of 27 top genes was selected as highly promising candidate genes for genetic resistance to MD. These genes are functionally involved with virus process (F13A1 and HSP90AB1), immunity (ABCB1LB, RGS5, C10ORF58, OSF-2, MMP7, CXCL12, GAL1, GAL2, GAL7, HVCN1, PDE4D, IL4I1, PARP9, EOMES, MPEG1, PDK4, CCLI10, K60 and FST), and tumor suppression (ADAMTS2, LXN, ARRDC3, WNT7A, CLDN1 and HPGD). It is anticipated that these findings will facilitate advancement in the fundamental understanding on mechanisms of genetic resistance to MD. In addition, such advancement may also provide insights on tumor virus-induced tumorigenesis in general and help the research community recognize MD study may serve as a good model for oncology study involving tumor viruses.

Allele-Specific Expression of CD4+ T Cells in Response to Marek's Disease Virus Infection

Marek's disease (MD) is a T cell lymphoma disease induced by Marek's disease virus (MDV), a highly oncogenic α herpesvirus primarily affecting chickens. MD is a chronic infectious disease that threatens the poultry industry. However, the mechanisms of genetic resistance for MD are complex and not completely understood. In this study, to identify high-confidence candidate genes of MD genetic resistance, high throughput sequencing (RNA-seq) was used to obtain transcriptomic data of CD4+ T cells isolated from MDV-infected and non-infected groups of two reciprocal crosses of individuals mating by two highly inbred chicken lines (63 MD-resistant and 72 MD-susceptible). After RNA-seq analysis with two biological replicates in each group, we identified 61 and 123 single nucleotide polymorphisms (SNPs) (false discovery rate (FDR) < 0.05) annotated in 39 and 132 genes in intercrosses 63 × 72 and 72 × 63, respectively, which exhibited allele-specific expression (ASE) in response to MDV infection. Similarly, we identified 62 and 79 SNPs annotated in 66 and 96 genes in infected and non-infected groups, respectively. We identified 534 and 1543 differentially expressed genes (DEGs) (FDR < 0.05) related to MDV infection in intercrosses 63 × 72 and 72 × 63, respectively. We also identified 328 and 20 DEGs in infected and non-infected groups, respectively. The qRT-PCR using seven DEGs further verified our results of RNA-seq analysis. The qRT-PCR of 11 important ASE genes was performed for gene functional validation in CD4+ T cells and tumors. Combining the analyses, six genes (MCL1, SLC43A2, PDE3B, ADAM33, BLB1, and DMB2), especially MCL1, were highlighted as the candidate genes with the potential to be involved in MDV infection. Gene-set enrichment analysis revealed that many ASE genes are linked to T cell activation, T cell receptor (TCR), B cell receptor (BCR), ERK/MAPK, and PI3K/AKT-mTOR signaling pathways, which play potentially important roles in MDV infection. Our approach underlines the importance of comprehensive functional studies for gaining valuable biological insight into the genetic factors behind MD and other complex traits, and our findings provide additional insights into the mechanisms of MD and disease resistance breeding in poultry.

Avian Expression Patterns and Genomic Mapping Implicate Leptin in Digestion and TNF in Immunity, Suggesting That Their Interacting Adipokine Role Has Been Acquired Only in Mammals

In mammals, leptin and tumor-necrosis factor (TNF) are prominent interacting adipokines mediating appetite control and insulin sensitivity. While TNF pleiotropically functions in immune defense and cell survival, leptin is largely confined to signaling energy stores in adipocytes. Knowledge about the function of avian leptin and TNF is limited and they are absent or lowly expressed in adipose, respectively. Employing radiation-hybrid mapping and FISH-TSA, we mapped TNF and its syntenic genes to chicken chromosome 16 within the major histocompatibility complex (MHC) region. This mapping position suggests that avian TNF has a role in regulating immune response. To test its possible interaction with leptin within the immune system and beyond, we compared the transcription patterns of TNF, leptin and their cognate receptors obtained by meta-analysis of GenBank RNA-seq data. While expression of leptin and its receptor (LEPR) were detected in the brain and digestive tract, TNF and its receptor mRNAs were primarily found in viral-infected and LPS-treated leukocytes. We confirmed leptin expression in the duodenum by immunohistochemistry staining. Altogether, we suggest that whereas leptin and TNF interact as adipokines in mammals, in birds, they have distinct roles. Thus, the interaction between leptin and TNF may be unique to mammals.

Towards a mechanistic understanding of the synergistic response induced by bivalent Marek's disease vaccines to prevent lymphomas

BACKGROUND

Marek's disease (MD) is a lymphoproliferative disease of chickens caused by Marek's disease virus (MDV), an oncogenic α-herpesvirus. Since 1970, MD has been controlled by widespread vaccination; however, more effective MD vaccines are needed to counter more virulent MDV strains. The bivalent vaccine combination of SB-1 and herpesvirus of turkey (HVT) strain FC126 has been widely used. Nonetheless, the mechanism(s) underlying this synergistic effect has not been investigated.

METHODS

Three experiments were conducted where SB-1 or HVT were administered as monovalent or bivalent vaccines to newly hatched chickens, then challenged five days later with MDV. In Experiment 1, levels of MDV replication in PBMCs were measured over time, and tumor incidence and vaccinal protection determined. In Experiment 2, MDV and vaccine strains replication levels in lymphoid organs were measured at 1, 5, 10, and 14 days post-challenge (DPC). In Experiment 3, to verify that the bursa was necessary for HVT protection, a subset of chicks were bursectomized and these birds plus controls were similarly vaccinated and challenged, and the levels of vaccinal protection determined.

RESULTS

The efficacy of bivalent SB-1 + HVT surpasses that of either SB-1 or HVT monovalent vaccines in controlling the level of pathogenic MDV in PBMCs until the end of the study, and this correlated with the ability to inhibit tumor formation. SB-1 replication in the spleen increased from 1 to 14 DPC, while HVT replicated only in the bursa at 1 DPC. The bursa was necessary for immune protection induced by HVT vaccine.

CONCLUSION

Synergy of SB-1 and HVT vaccines is due to additive influences of the individual vaccines acting at different times and target organs. And the bursa is vital for HVT to replicate and induce immune protection.

RNA-seq analysis of viral gene expression in the skin of Marek's disease virus infected chickens

Marek's disease virus (MDV), a highly cell-associated oncogenic avian α-herpesvirus, is the causative agent of malignant transformation of T cells in domestic chickens. The latently infected CD4⁺CD8^- T cells carry the virus through the blood stream and establish lymphomas in the skin, visceral organs and peripheral nerves. The feather follicle epithelium (FFE) is the only anatomical site where fully infectious enveloped virions are produced and eventually disseminated into the environment to infect contact birds. Therefore, skin and FFE play a critical role as being the common source of re-infection of birds sharing the same habitat. The molecular mechanism involved in the replication and assembly of MDV in the FFE leading to the production and release of cell-free infectious virus particles is unknown and to date no viral or host gene has been implicated in the process. To examine alterations in the expression pattern of viral genes, we performed RNA-seq on the skin samples of Marek's disease virus-infected susceptible chickens at 10, 20, and 30 days post infection. For comparative analysis of the expression patterns of viral genes between the skin and spleen of the MD-susceptible and resistant lines, Real-Time RT-PCR was employed. In total, RNA-seq based analysis identified 42 viral genes that were differentially expressed in the skin of infected birds. Majority of the identified genes are involved in DNA replication, capsid, tegument, and envelop formation. Comparative analysis between the skin and spleen of MD-susceptible and resistant chicken lines, revealed significantly higher expression of the genes in the skin of either lines than the spleen. Furthermore, much higher expression of the genes was observed in the skin of the susceptible line than the resistant line.

Marek's Disease Virus Infection Induced Mitochondria Changes in Chickens

Mitochondria are crucial cellular organelles in eukaryotes and participate in many cell processes including immune response, growth development, and tumorigenesis. Marek’s disease (MD), caused by an avian alpha-herpesvirus Marek’s disease virus (MDV), is characterized with lymphomas and immunosuppression. In this research, we hypothesize that mitochondria may play roles in response to MDV infection. To test it, mitochondrial DNA (mtDNA) abundance and gene expression in immune organs were examined in two well-defined and highly inbred lines of chickens, the MD-susceptible line 7₂ and the MD-resistant line 6₃. We found that mitochondrial DNA contents decreased significantly at the transformation phase in spleen of the MD-susceptible line 7₂ birds in contrast to the MD-resistant line 6₃. The mtDNA-genes and the nucleus-genes relevant to mtDNA maintenance and transcription, however, were significantly up-regulated. Interestingly, we found that POLG2 might play a potential role that led to the imbalance of mtDNA copy number and gene expression alteration. MDV infection induced imbalance of mitochondrial contents and gene expression, demonstrating the indispensability of mitochondria in virus-induced cell transformation and subsequent lymphoma formation, such as MD development in chicken. This is the first report on relationship between virus infection and mitochondria in chicken, which provides important insights into the understanding on pathogenesis and tumorigenesis due to viral infection.

Identification of Marek's disease virus genes associated with virulence of US strains

Marek's disease virus (MDV) is the most well-cited example of vaccine-driven virulence evolution. MDV induces a lymphoproliferative disease in chickens, which is currently controlled by widespread vaccination of flocks. Unfortunately, Marek's disease (MD) vaccines, while effective in preventing tumours, do not prevent viral replication and mutation, which has been hypothesized as the major driving force for increased MDV virulence of field strains during the past 40 years in US commercial flocks. To limit future virulence increases, there is interest in characterizing MDV strain genomes collected over the years and associating genetic variations with variation in virulence. In this study, we characterized 70 MDV genomes with known virulence by complete or targeted DNA sequencing, and identified genetic variants that showed association with virulence. Our results revealed a number of MDV genes as would be expected for a complex trait. In addition, phylogenetic analysis revealed a clear separation of strains that varied by virulence. Interestingly, high virulence isolates from the same farms persisted over years despite eradication attempts, which has implications on control efforts. Given the growing ability to bioengineer the MDV genome, it should be feasible to experimentally test whether these individual variants influence virulence markers alone or combinations. Once validated, these markers may provide an alternative to live bird testing for evaluating virulence of new MDV field strains.

Avian Sarcoma and Leukosis Virus Envelope Glycoproteins Evolve to Broaden Receptor Usage Under Pressure from Entry Competitors †

The subgroup A through E avian sarcoma and leukosis viruses (ASLV(A) through ASLV(E)) are a group of highly related alpharetroviruses that have evolved their envelope glycoproteins to use different receptors to enable efficient virus entry due to host resistance and/or to expand host range. Previously, we demonstrated that ASLV(A) in the presence of a competitor to the subgroup A Tva receptor, SUA-rIgG immunoadhesin, evolved to use other receptor options. The selected mutant virus, RCASBP(A)Δ155–160, modestly expanded its use of the Tvb and Tvc receptors and possibly other cell surface proteins while maintaining the binding affinity to Tva. In this study, we further evolved the Δ155–160 virus with the genetic selection pressure of a soluble form of the Tva receptor that should force the loss of Tva binding affinity in the presence of the Δ155–160 mutation. Viable ASLVs were selected that acquired additional mutations in the Δ155–160 Env hypervariable regions that significantly broadened receptor usage to include Tvb and Tvc as well as retaining the use of Tva as a receptor determined by receptor interference assays. A similar deletion in the hr1 hypervariable region of the subgroup C ASLV glycoproteins evolved to broaden receptor usage when selected on Tvc-negative cells.

Reverse Engineering Provides Insights on the Evolution of Subgroups A to E Avian Sarcoma and Leukosis Virus Receptor Specificity

The initial step of retrovirus entry—the interaction between the virus envelope glycoprotein trimer and a cellular receptor—is complex, involving multiple, noncontiguous determinants in both proteins that specify receptor choice, binding affinity and the ability to trigger conformational changes in the viral glycoproteins. Despite the complexity of this interaction, retroviruses have the ability to evolve the structure of their envelope glycoproteins to use a different cellular protein as receptors. The highly homologous subgroup A to E Avian Sarcoma and Leukosis Virus (ASLV) glycoproteins belong to the group of class 1 viral fusion proteins with a two-step triggering mechanism that allows experimental access to intermediate structures during the fusion process. We and others have taken advantage of replication-competent ASLVs and exploited genetic selection strategies to force the ASLVs to naturally evolve and acquire envelope glycoprotein mutations to escape the pressure on virus entry and still yield a functional replicating virus. This approach allows for the simultaneous selection of multiple mutations in multiple functional domains of the envelope glycoprotein that may be required to yield a functional virus. Here, we review the ASLV family and experimental system and the reverse engineering approaches used to understand the evolution of ASLV receptor usage.

Deciphering desirable immune responses from disease models with resistant and susceptible chickens

Coccidiosis and necrotic enteritis (NE) are among the most significant diseases affecting the poultry industry. These diseases have become more prominent in the wake of policies to reduce the use of antibiotics in animal production. This has led to more research focused on better understanding the immune system and its responses to pathogen challenge, and thus developing informed strategies to exploit immune responses that can support enhanced disease resistance and growth performance. Some chicken breeds and lines show greater resistance or susceptibility to various diseases, and thus these birds maybe able to shed light on immune processes or pathways that contribute to the more resistant/susceptible state. This review attempts to identify potentially important genes that show some consistency in (relative) up or downregulation in key tissues between the resistant and susceptible chickens. For coccidiosis and NE, relative downregulation of IL-10 and (slightly less consistently) upregulation of IFN-γ appear to be features of more resistant birds. Data for IFN-α, IL-12, and IL-17D are currently less consistent. Gene expression data from NE studies have identified some potentially interesting, perhaps less well understood, immune-related genes (e.g., TCF12, BCL2, IRF2, TRAF3, TAB3, etc.,) that maybe associated with the resistant and/or susceptible phenotype. Salmonella and Campylobacter are important foodborne pathogens harbored by the chicken intestinal tract, while infectious bursal disease and infectious bronchitis are also important viral diseases of poultry. We, therefore, consider whether there are consistent features from resistant/susceptible disease models with these pathogens that relate to findings from the coccidiosis and NE studies. It is not anticipated that ideal immune responses to these pathogens will be identical but rather that consistent elements maybe identified that could help inform breeding or alternative strategies to support general disease resistance and enhanced (and efficient) flock productivity.

Vaccination and Host Marek's Disease-Resistance Genotype Significantly Reduce Oncogenic Gallid alphaherpesvirus 2 Telomere Integration in Host Birds

Marek's disease (MD) is an infectious disease characterized by lymphomas and high mortality in susceptible chickens. The causative and ubiquitous alpha-herpesvirus known as MD virus (MDV) integrates into host telomeres during early infection through latency, known to be an important phase for oncogenic transformation. Herein, we sought to determine the influence of vaccination and host genetics on the temporal dynamics of MDV-host genome interactions. We studied integration profiles using 2 MD vaccines that vary in protective efficacy in 2 genetic lines that differ in MD resistance/susceptibility. Virus integration of both oncogenic MDV and vaccine strains was observed in both MD susceptible and resistant birds, however, the lines differed in their dynamic telomere-integration profiles. Notably, the resistant host genotype exhibited a smaller percentage of replicating cells with the virus telomere-integrated only phenotype as compared to the susceptible genotype. Vaccination with Rispens, the most protective MD vaccine, also reduced the establishment of the virus telomere-integrated only phenotype, suggesting a significant role of the phenotype in MD lymphoma development. The effect of Rispens vaccination was most dramatic in the susceptible genotype. These results suggest important connections between vaccinal immunity, MDV telomere integration, virus-induced oncogenesis, and virus-host genome interactions in the context of host genetics and disease susceptibility.

Marek's disease vaccine activates chicken macrophages

To provide insights into the role of innate immune responses in vaccine-mediated protection, we investigated the effect of Marek's disease (MD) vaccine, CVI988/Rispens, on the expression patterns of selected genes associated with activation of macrophages in MD-resistant and MD-susceptible chicken lines. Upregulation of interferon γ, interleukin (IL)-1β, IL-8, and IL-12 at different days post-inoculation (dpi) revealed activation of macrophages in both chicken lines. A strong immune response was induced in cecal tonsils of the susceptible line at 5 dpi. The highest transcriptional activities were observed in spleen tissues of the resistant line at 3 dpi. No increase in the population of CD3⁺ T cells was observed in duodenum of vaccinated birds at 5 dpi indicating a lack of involvement of the adaptive immune system in the transcriptional profiling of the tested genes. There was, however, an increase in the number of macrophages in the duodenum of vaccinated birds. The CVI988/Rispens antigen was detected in the duodenum and cecal tonsils of the susceptible line at 5 dpi but not in the resistant line. This study sheds light on the role of macrophages in vaccine-mediated protection against MD and on the possible development of new recombinant vaccines with enhanced innate immune system activation properties.

Genetic assessment of inbred chicken lines indicates genomic signatures of resistance to Marek's disease

Background

Marek’s disease (MD) is a highly contagious pathogenic and oncogenic disease primarily affecting chickens. However, the mechanisms of genetic resistance for MD are complex and not fully understood. MD-resistant line 6₃ and MD-susceptible line 7₂ are two highly inbred progenitor lines of White Leghorn. Recombinant Congenic Strains (RCS) were developed from these two lines, which show varied susceptibility to MD.

Results

We investigated genetic structure and genomic signatures across the genome, including the line 6₃ and line 7₂, six RCSs, and two reciprocally crossed flocks between the lines 6₃ and 7₂ (F1 6₃ × 7₂ and F1 7₂ × 6₃) using Affymetrix® Axiom® HD 600 K genotyping array. We observed 18 chickens from RCS lines were specifically clustered into resistance sub-groups distributed around line 6₃. Additionally, homozygosity analysis was employed to explore potential genetic components related to MD resistance, while runs of homozygosity (ROH) are regions of the genome where the identical haplotypes are inherited from each parent. We found several genes including SIK, SOX1, LIG4, SIK1 and TNFSF13B were contained in ROH region identified in resistant group (line 6₃ and RCS), and these genes have been reported that are contribute to immunology and survival. Based on F_ST based population differential analysis, we also identified important genes related to cell death and anti-apoptosis, including AKT1, API5, CDH13, CFDP and USP15, which could be involved in divergent selection during inbreeding process.

Conclusions

Our findings offer valuable insights for understanding the genetic mechanism of resistance to MD and the identified genes could be considered as candidate biomarkers in further evaluation.

Electronic supplementary material

The online version of this article (10.1186/s40104-018-0281-x) contains supplementary material, which is available to authorized users.

Evaluation of genetic resistance to Salmonella Pullorum in three chicken lines

Resistance to diseases varies considerably among populations of the same species and can be ascribed to both genetic and environmental factors. Salmonella Pullorum (SP) is responsible for significant losses in the poultry industry, especially in developing countries. To better understand SP resistance in chicken populations with different genetic backgrounds, we orally challenged 3 chicken lines with SP-a highly selected commercial breed (Rhode Island Red, RIR), a local Chinese chicken (Beijing You, BY), and a synthetic layer line (dwarf, DW)-at 4 d of age. Two traits related to SP resistance, survival, and bacterial carriage in the spleen were evaluated after infection. Survival rates were recorded up to 40 d of age when all chickens still alive were killed to verify the presence of SP in the spleen to determine carrier state. Mortalities for RIR, BY, and DW chicks were 25.1%, 8.3%, and 22.7%, respectively, and the corresponding carrier-states in the spleens were 17.9%, 0.6%, and 15.8%. Survival and carrier-state heritabilities were estimated using an animal threshold model. Survival heritability was 0.197, 0.091, and 0.167 in RIR, BY, and DW populations, respectively, and the heritabilities of carrier state for DW and RIR were 0.32 and 0.16, respectively. This is the first time that the heritability of the SP carrier state has been evaluated in chickens. Our study provides experimental evidence that chickens with various genetic background exhibited significantly different SP-resistant activities and heritabilities. These results may be useful for selecting lines with better disease resistance.

Development of reliable techniques for the differential diagnosis of avian tumour viruses by immunohistochemistry and polymerase chain reaction from formalin-fixed paraffin-embedded tissue sections

A variety of techniques have been developed as diagnostic tools for the differential diagnosis of tumours produced by Marek's disease virus from those induced by avian leukosis virus and reticuloendotheliosis virus. However, most current techniques are unreliable when used in formalin-fixed paraffin-embedded (FFPE) tissues, which often is the only sample type available for definitive diagnosis. A collection of tumours was generated by the inoculation of different strains of Marek's disease virus, reticuloendotheliosis virus or avian leukosis virus singularly or in combination. FFPE tissue sections from tumour and non-tumour tissues were analysed by optimized immunohistochemistry (IHC) techniques and traditional as well as quantitative polymerase chain reaction (PCR) with newly designed primers ideal for DNA fragmented by fixation. IHC and PCR results were highly sensitive and specific in tissues from single-infected birds. Virus quantity was higher in tumours compared to non-tumour spleens from Marek's disease (MD) virus-infected birds. Thus, using FFPE sections alone may be sufficient for the diagnosis of MD by demonstration of high quantities of viral antigens or genome in tumour cells, along with the absence of other tumour viruses by traditional PCR, and if standard criteria are met based on clinical history and histology. IHC furthermore allowed detection of the specific cells that were infected with different viruses in tumours from birds that had been inoculated simultaneously with multiple viruses. Following validation with field samples, these new protocols can be applied for both diagnostic and research purposes to help accurately identify avian tumour viruses in routine FFPE tissue sections.