Improving the outcome of a Marek's disease challenge in multiple lines of egg type chickens

Purinergic signaling during Marek's disease in chickens

Purinergic receptors (PRs) have been reported as potential therapeutic targets for many viral infections including herpesviruses, which urges the investigation into their role in Marek's disease (MD), a herpesvirus induced cancer in chickens that is an important pathogen for the poultry industry. MD is caused by MD virus (MDV) that has a similar viral life cycle as human varicella zoster virus in that it is shed from infected epithelial skin cells and enters the host through the respiratory route. In this report, PR responses during natural MDV infection and disease progression was examined in MD-resistant white Leghorns (WL) and MD-susceptible Pure Columbian (PC) chickens during natural infection. Whole lung lavage cells (WLLC) and liver tissue samples were collected from chickens infected but showing no clinical signs of MD (Infected) or presenting with clinical disease (Diseased). RNA was extracted followed by RT-qPCR analysis with gene specific primers against members of the P1, P2X, and P2Y PR families. Differential expression (p < 0.05) was observed in breed and disease conditions. Some PRs showed tissue specific expression (P1A1, P2X1, and P2X6 in WLLC) whereas others responded to MDV infection only in MD-susceptible (PC) chickens (P1A2A, P2X1, P2X5, P2X7). P2Y PRs had differential expression in both chicken lines in response to MDV infection and MD progression. This study is the first to our knowledge to examine PR responses during MDV infection and disease progression. These results suggest PR signaling may an important area of research for MDV replication and MD.

Sex Differences in Response to Marek's Disease: Mapping Quantitative Trait Loci Regions (QTLRs) to the Z Chromosome

Marek's Disease (MD) has a significant impact on both the global poultry economy and animal welfare. The disease pathology can include neurological damage and tumour formation. Sexual dimorphism in immunity and known higher susceptibility of females to MD makes the chicken Z chromosome (GGZ) a particularly attractive target to study the chicken MD response. Previously, we used a Hy-Line F₆ population from a full-sib advanced intercross line to map MD QTL regions (QTLRs) on all chicken autosomes. Here, we mapped MD QTLRs on GGZ in the previously utilized F₆ population with individual genotypes and phenotypes, and in eight elite commercial egg production lines with daughter-tested sires and selective DNA pooling (SDP). Four MD QTLRs were found from each analysis. Some of these QTLRs overlap regions from previous reports. All QTLRs were tested by individuals from the same eight lines used in the SDP and genotyped with markers located within and around the QTLRs. All QTLRs were confirmed. The results exemplify the complexity of MD resistance in chickens and the complex distribution of p-values and Linkage Disequilibrium (LD) pattern and their effect on localization of the causative elements. Considering the fragments and interdigitated LD blocks while using LD to aid localization of causative elements, one must look beyond the non-significant markers, for possible distant markers and blocks in high LD with the significant block. The QTLRs found here may explain at least part of the gender differences in MD tolerance, and provide targets for mitigating the effects of MD.

Genetic Analyses of Response of Local Ghanaian Tanzanian Chicken Ecotypes to a Natural Challenge with Velogenic Newcastle Disease Virus

Simple Summary

Newcastle disease (ND) is a global threat to poultry production and often has a major impact on chicken welfare and the livelihoods of rural poultry farmers. We exposed unvaccinated Ghanaian and Tanzanian local chicken ecotypes to velogenic Newcastle disease virus strains, and measured response traits to understand the genetic basis of ND. We identified heritable ND response traits and revealed differences in survival between Ghanaian and Tanzanian local chicken ecotype birds. Our findings indicate that velogenic ND resistance could be improved through selective breeding of local chicken ecotypes in regions where the disease is endemic.

Abstract

Newcastle disease is a devastating poultry disease that often causes significant economic losses in poultry in the developing countries of Africa, Asia, as well as South and Central America. Velogenic Newcastle disease virus (NDV) outbreaks are associated with high mortalities, which can threaten household livelihoods, especially in the rural areas, and lead to loss of high-quality proteins in the form of meat and eggs, as well as household purchasing power. In this study, we exposed unvaccinated Ghanaian and Tanzanian chickens of six local ecotypes to velogenic NDV strains, measured NDV response traits, sequenced their DNA on a genotyping-by-sequencing platform, and performed variance component analyses. The collected phenotypes included: growth rates (pre- and post-exposure); lesion scores (gross lesion severity) in the trachea, proventriculus, intestine, and cecal tonsils; natural antibody levels; anti-NDV antibody levels at 7 days post exposure (dpe); tear and cloacal viral load at 2, 4, and 6 dpe; and survival time. Heritability estimates were low to moderate, ranging from 0.11 for average lesion scores to 0.36 for pre-exposure growth rate. Heritability estimates for survival time were 0.23 and 0.27 for the Tanzanian and Ghanaian ecotypes, respectively. Similar heritability estimates were observed when data were analyzed either separately or combined for the two countries. Survival time was genetically negatively correlated with lesion scores and with viral load. Results suggested that response to mesogenic or velogenic NDV of these local chicken ecotypes could be improved by selective breeding. Chickens that are more resilient to velogenic NDV can improve household livelihoods in developing countries.

Why breed disease-resilient livestock, and how?

BACKGROUND

Fighting and controlling epidemic and endemic diseases represents a considerable cost to livestock production. Much research is dedicated to breeding disease resilient livestock, but this is not yet a common objective in practical breeding programs. In this paper, we investigate how future breeding programs may benefit from recent research on disease resilience.

MAIN BODY

We define disease resilience in terms of its component traits resistance (R: the ability of a host animal to limit within-host pathogen load (PL)) and tolerance (T: the ability of an infected host to limit the damage caused by a given PL), and model the host's production performance as a reaction norm on PL, depending on R and T. Based on this, we derive equations for the economic values of resilience and its component traits. A case study on porcine respiratory and reproductive syndrome (PRRS) in pigs illustrates that the economic value of increasing production in infectious conditions through selection for R and T can be more than three times higher than by selection for production in disease-free conditions. Although this reaction norm model of resilience is helpful for quantifying its relationship to its component traits, its parameters are difficult and expensive to quantify. We consider the consequences of ignoring R and T in breeding programs that measure resilience as production in infectious conditions with unknown PL-particularly, the risk that the genetic correlation between R and T is unfavourable (antagonistic) and that a trade-off between them neutralizes the resilience improvement. We describe four approaches to avoid such antagonisms: (1) by producing sufficient PL records to estimate this correlation and check for antagonisms-if found, continue routine PL recording, and if not found, shift to cheaper proxies for PL; (2) by selection on quantitative trait loci (QTL) known to influence both R and T in favourable ways; (3) by rapidly modifying towards near-complete resistance or tolerance, (4) by re-defining resilience as the animal's capacity to resist (or recover from) the perturbation caused by an infection, measured as temporal deviations of production traits in within-host longitudinal data series.

CONCLUSIONS

All four alternatives offer promising options for genetic improvement of disease resilience, and most rely on technological and methodological developments and innovation in automated data generation.

Mapping QTL Associated with Resistance to Avian Oncogenic Marek's Disease Virus (MDV) Reveals Major Candidate Genes and Variants

Marek’s disease (MD) represents a significant global economic and animal welfare issue. Marek’s disease virus (MDV) is a highly contagious oncogenic and highly immune-suppressive α-herpes virus, which infects chickens, causing neurological effects and tumour formation. Though partially controlled by vaccination, MD continues to have a profound impact on animal health and on the poultry industry. Genetic selection provides an alternative and complementary method to vaccination. However, even after years of study, the genetic mechanisms underlying resistance to MDV remain poorly understood. The Major Histocompatability Complex (MHC) is known to play a role in disease resistance, along with a handful of other non-MHC genes. In this study, one of the largest to date, we used a multi-facetted approach to identify quantitative trait locus regions (QTLR) influencing resistance to MDV, including an F₆ population from a full-sib advanced intercross line (FSIL) between two elite commercial layer lines differing in resistance to MDV, RNA-seq information from virus challenged chicks, and genome wide association study (GWAS) from multiple commercial lines. Candidate genomic elements residing in the QTLR were further tested for association with offspring mortality in the face of MDV challenge in eight pure lines of elite egg-layer birds. Thirty-eight QTLR were found on 19 chicken chromosomes. Candidate genes, microRNAs, long non-coding RNAs and potentially functional mutations were identified in these regions. Association tests were carried out in 26 of the QTLR, using eight pure lines of elite egg-layer birds. Numerous candidate genomic elements were strongly associated with MD resistance. Genomic regions significantly associated with resistance to MDV were mapped and candidate genes identified. Various QTLR elements were shown to have a strong genetic association with resistance. These results provide a large number of significant targets for mitigating the effects of MDV infection on both poultry health and the economy, whether by means of selective breeding, improved vaccine design, or gene-editing technologies.

Advances in methodologies for detecting MHC-B variability in chickens

The chicken major histocompatibility B complex (MHC-B) region is of great interest owing to its very strong association with resistance to many diseases. Variation in the MHC-B was initially identified by hemagglutination of red blood cells with specific alloantisera. New technologies, developed to identify variation in biological materials, have been applied to the chicken MHC. Protein variation encoded by the MHC genes was examined by immunoprecipitation and 2-dimensional gel electrophoresis. Increased availability of DNA probes, PCR, and sequencing resulted in the application of DNA-based methods for MHC detection. The chicken reference genome, completed in 2004, allowed further refinements in DNA methods that enabled more rapid examination of MHC variation and extended such analyses to include very diverse chicken populations. This review progresses from the inception of MHC-B identification to the present, describing multiple methods, plus their advantages and disadvantages.

Genomic selection for the improvement of antibody response to Newcastle disease and avian influenza virus in chickens

Newcastle disease (ND) and avian influenza (AI) are the most feared diseases in the poultry industry worldwide. They can cause flock mortality up to 100%, resulting in a catastrophic economic loss. This is the first study to investigate the feasibility of genomic selection for antibody response to Newcastle disease virus (Ab-NDV) and antibody response to Avian Influenza virus (Ab-AIV) in chickens. The data were collected from a crossbred population. Breeding values for Ab-NDV and Ab-AIV were estimated using a pedigree-based best linear unbiased prediction model (BLUP) and a genomic best linear unbiased prediction model (GBLUP). Single-trait and multiple-trait analyses were implemented. According to the analysis using the pedigree-based model, the heritability for Ab-NDV estimated from the single-trait and multiple-trait models was 0.478 and 0.487, respectively. The heritability for Ab-AIV estimated from the two models was 0.301 and 0.291, respectively. The estimated genetic correlation between the two traits was 0.438. A four-fold cross-validation was used to assess the accuracy of the estimated breeding values (EBV) in the two validation scenarios. In the family sample scenario each half-sib family is randomly allocated to one of four subsets and in the random sample scenario the individuals are randomly divided into four subsets. In the family sample scenario, compared with the pedigree-based model, the accuracy of the genomic prediction increased from 0.086 to 0.237 for Ab-NDV and from 0.080 to 0.347 for Ab-AIV. In the random sample scenario, the accuracy was improved from 0.389 to 0.427 for Ab-NDV and from 0.281 to 0.367 for Ab-AIV. The multiple-trait GBLUP model led to a slightly higher accuracy of genomic prediction for both traits. These results indicate that genomic selection for antibody response to ND and AI in chickens is promising.

Genetic variation within the Mx gene of commercially selected chicken lines reveals multiple haplotypes, recombination and a protein under selection pressure

The Mx protein is one of the best-characterized interferon-stimulated antiviral mediators. Mx homologs have been identified in most vertebrates examined; however, their location within the cell, their level of activity, and the viruses they inhibit vary widely. Recent studies have demonstrated multiple Mx alleles in chickens and some reports have suggested a specific variant (S631N) within exon 14 confers antiviral activity. In the current study, the complete genome of nine elite egg-layer type lines were sequenced and multiple variants of the Mx gene identified. Within the coding region and upstream putative promoter region 36 SNP variants were identified, producing a total of 12 unique haplotypes. Each elite line contained from one to four haplotypes, with many of these haplotypes being found in only one line. Observation of changes in haplotype frequency over generations, as well as recombination, suggested some unknown selection pressure on the Mx gene. Trait association analysis with either individual SNP or haplotypes showed a significant effect of Mx haplotype on several egg production related traits, and on mortality following Marek's disease virus challenge in some lines. Examination of the location of the various SNP within the protein suggests synonymous SNP tend to be found within structural or enzymatic regions of the protein, while non-synonymous SNP are located in less well defined regions. The putative resistance variant N631 was found in five of the 12 haplotypes with an overall frequency of 47% across the nine lines. Two Mx recombinants were identified within the elite populations, indicating that novel variation can arise and be maintained within intensively selected lines. Collectively, these results suggest the conflicting reports in the literature describing the impact of the different SNP on chicken Mx function may be due to the varying context of haplotypes present in the populations studied.

Genome-wide association study for Marek's disease mortality in layer chickens

A genome-wide association study (GWAS) using Bayesian variable selection was performed to determine genomic regions associated with mortality due to Marek's disease virus (MDV) infection in layers. Mortality (%) under experimental disease challenge (500 plaque-forming units of a very virulent plus MDV strain) was recorded for progeny groups (average 15.5 birds; range 3 to 30) of 253 genotyped sires from four generations of a brown-egg layer line. An additional generation of 43 sires with progeny data was used to validate results. Sires were genotyped with a 42K Illumina single-nucleotide polymorphism (SNP) chip. Methods BayesB (pi = 0.995) and BayesCpi, with or without weighting residuals by the size of progeny groups were applied. The proportion of genetic variance contributed by SNPs within each 1-megabase (Mb) genomic region was quantified. Average mortality was 33% but differed significantly between generations. Genetic markers explained about 11% of phenotypic variation in mortality. Correlations between genomic estimated breeding values and percentage of progeny mortality for the validation generation (sons of individuals in training) were 0.12, 0.17, 0.02, and 0.16 for BayesB, weighted BayesB, BayesCpi, and weighted BayesCpi, respectively, when using the whole genome, and 0.03, 0.20, -0.06, and 0.14, when using only SNP from the 10, 1-Mb regions, explaining the largest proportion of genetic variance according to each method. Results suggest that regions on chromosomes 2, 3, 4, 9, 15, 18, and 21 are associated with Marek's disease resistance and can be used for selection and that accounting for the size of progeny groups has a large impact on correct localization of such genomic regions.