Mapping QTL affecting resistance to Marek's disease in an F6 advanced intercross population of commercial layer chickens

Research Note: Phenotypic trends for the multigenerational advanced intercross of the Virginia body weight lines of chickens

Random samples from generation S41 of the Virginia high and low 8-week body weight lines formed the base population for producing a multigenerational reciprocal intercross population. Although genetic mapping from this intercross has been reported, lacking are phenotypic trends across multiple generations. Here, we provide phenotypic information for the parental base population, the F1 reciprocal cross, and subsequent segregating recombinant generations F2 to F17. Heterosis for the selected trait in the F1 was negative for both reciprocal crosses. Phenotypic correlations for the selected trait in the recombinant generations were essentially nil for both males and females as was percent sexual dimorphism and coefficients of variation.

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.

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.

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.

The Chicken A and E Blood Systems Arise from Genetic Variation in and around the Regulators of Complement Activation Region

The tightly linked A and E blood alloantigen systems are 2 of 13 blood systems identified in chickens. Reported herein are studies showing that the genes encoding A and E alloantigens map within or near to the chicken regulator of complement activation (RCA) gene cluster, a region syntenic with the human RCA. Genome-wide association studies, sequence analysis, and sequence-derived single-nucleotide polymorphism information for known A and/or E system alleles show that the most likely candidate gene for the A blood system is C4BPM gene (complement component 4 binding protein, membrane). Cosegregation of single-nucleotide polymorphism–defined C4BPM haplotypes and blood system A alleles defined by alloantisera provide a link between chicken blood system A and C4BPM. The best match for the E blood system is the avian equivalent of FCAMR (Fc fragment of IgA and IgM receptor). C4BPM is located within the chicken RCA on chicken microchromosome 26 and is separated from FCAMR by 89 kbp. The genetic variation observed at C4BPM and FCAMR could affect the chicken complement system and differentially guide immune responses to infectious diseases.

The Diverse Major Histocompatibility Complex Haplotypes of a Common Commercial Chicken Line and Their Effect on Marek's Disease Virus Pathogenesis and Tumorigenesis

The major histocompatibility complex (MHC) is crucial for appropriate immune responses against invading pathogens. Chickens possess a single predominantly-expressed class I molecule with strong associations between disease resistance and MHC haplotype. For Marek's disease virus (MDV) infections of chickens, the MHC haplotype is one of the major determinants of genetic resistance and susceptibility. VALO specific pathogen free (SPF) chickens are widely used in biomedical research and vaccine production. While valuable findings originate from MDV infections of VALO SPF chickens, their MHC haplotypes and associated disease resistance remained elusive. In this study, we used several typing systems to show that VALO SPF chickens possess MHC haplotypes that include B9, B9:02, B15, B19 and B21 at various frequencies. Moreover, we associate the MHC haplotypes to MDV-induced disease and lymphoma formation and found that B15 homozygotes had the lowest tumor incidence while B21 homozygotes had the lowest number of organs with tumors. Finally, we found transmission at variable levels to all contact birds except B15/B21 heterozygotes. These data have immediate implications for the use of VALO SPF chickens and eggs in the life sciences and add another piece to the puzzle of the chicken MHC complex and its role in infections with this oncogenic herpesvirus.

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.

From Chickens to Humans: The Importance of Peptide Repertoires for MHC Class I Alleles

In humans, killer immunoglobulin-like receptors (KIRs), expressed on natural killer (NK) and thymus-derived (T) cells, and their ligands, primarily the classical class I molecules of the major histocompatibility complex (MHC) expressed on nearly all cells, are both polymorphic. The variation of this receptor-ligand interaction, based on which alleles have been inherited, is known to play crucial roles in resistance to infectious disease, autoimmunity, and reproduction in humans. However, not all the variation in response is inherited, since KIR binding can be affected by a portion of the peptide bound to the class I molecules, with the particular peptide presented affecting the NK response. The extent to which the large multigene family of chicken immunoglobulin-like receptors (ChIRs) is involved in functions similar to KIRs is suspected but not proven. However, much is understood about the two MHC-I molecules encoded in the chicken MHC. The BF2 molecule is expressed at a high level and is thought to be the predominant ligand of cytotoxic T lymphocytes (CTLs), while the BF1 molecule is expressed at a much lower level if at all and is thought to be primarily a ligand for NK cells. Recently, a hierarchy of BF2 alleles with a suite of correlated properties has been defined, from those expressed at a high level on the cell surface but with a narrow range of bound peptides to those expressed at a lower level on the cell surface but with a very wide repertoire of bound peptides. Interestingly, there is a similar hierarchy for human class I alleles, although the hierarchy is not as wide. It is a question whether KIRs and ChIRs recognize class I molecules with bound peptide in a similar way, and whether fastidious to promiscuous hierarchy of class I molecules affect both T and NK cell function. Such effects might be different from those predicted by the similarities of peptide-binding based on peptide motifs, as enshrined in the idea of supertypes. Since the size of peptide repertoire can be very different for alleles with similar peptide motifs from the same supertype, the relative importance of these two properties may be testable.

Selection for Favorable Health Traits: A Potential Approach to Cope with Diseases in Farm Animals

Disease is a global problem for animal farming industries causing tremendous economic losses (>USD 220 billion over the last decade) and serious animal welfare issues. The limitations and deficiencies of current non-selection disease control methods (e.g., vaccination, treatment, eradication strategy, genome editing, and probiotics) make it difficult to effectively, economically, and permanently eliminate the adverse influences of disease in the farm animals. These limitations and deficiencies drive animal breeders to be more concerned and committed to dealing with health problems in farm animals by selecting animals with favorable health traits. Both genetic selection and genomic selection contribute to improving the health of farm animals by selecting certain health traits (e.g., disease tolerance, disease resistance, and immune response), although both of them face some challenges. The objective of this review was to comprehensively review the potential of selecting health traits in coping with issues caused by diseases in farm animals. Within this review, we highlighted that selecting health traits can be applied as a method of disease control to help animal agriculture industries to cope with the adverse influences caused by diseases in farm animals. Certainly, the genetic/genomic selection solution cannot solve all the disease problems in farm animals. Therefore, management, vaccination, culling, medical treatment, and other measures must accompany selection solution to reduce the adverse impact of farm animal diseases on profitability and animal welfare.

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.

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.

Comparative transcriptomics of genetically divergent lines of chickens in response to Marek's disease virus challenge at cytolytic phase

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

Characterization of Copy Number Variation's Potential Role in Marek's Disease

Marek’s Disease (MD) is a highly contagious pathogenic and oncogenic disease primarily affecting chickens. Chicken Lines 63 and 72, as well as their recombinant congenic strains (RCS) with varied susceptibility to MD, are ideal models to study the complex mechanisms of genetic resistance to MD. In this study, we investigated copy number variation (CNV) in these inbred chicken lines using the Affymetrix Axiom HD 600 K SNP genotyping array. We detected 393 CNV segments across all ten chicken lines, of which 12 CNVs were specifically identified in Line 72. We then assessed genetic structure based on CNV and observed markedly different patterns. Finally, we validated two deletion events in Line 72 and correlated them with genes expression using qPCR and RNA-seq, respectively. Our combined results indicated that these two CNV deletions were likely to contribute to MD susceptibility.

Genome-wide association studies of immune, disease and production traits in indigenous chicken ecotypes

Background

The majority of chickens in sub-Saharan Africa are indigenous ecotypes, well adapted to the local environment and raised in scavenging production systems. Although they are generally resilient to disease challenge, routine vaccination and biosecurity measures are rarely applied and infectious diseases remain a major cause of mortality and reduced productivity. Management and genetic improvement programmes are hampered by lack of routine data recording. Selective breeding based on genomic technologies may provide the means to enhance sustainability. In this study, we investigated the genetic architecture of antibody response to four major infectious diseases [infectious bursal disease (IBDV), Marek’s disease (MDV), fowl typhoid (SG), fowl cholera (PM)] and resistance to Eimeria and cestode parasitism, along with two production traits [body weight and body condition score (BCS)] in two distinct indigenous Ethiopian chicken ecotypes. We conducted variance component analyses, genome-wide association studies, and pathway and selective sweep analyses.

Results

The large majority of birds was found to have antibody titres for all pathogens and were infected with both parasites, suggesting almost universal exposure. We derived significant moderate to high heritabilities for IBDV, MDV and PM antibody titres, cestodes infestation, body weight and BCS. We identified single nucleotide polymorphisms (SNPs) with genome-wide significance for each trait. Based on these associations, we identified for each trait, pathways, networks and functional gene clusters that include plausible candidate genes. Selective sweep analyses revealed a locus on chromosome 18 associated with viral antibody titres and resistance to Eimeria parasitism that is within a positive selection signal. We found no significant genetic correlations between production, immune and disease traits, implying that selection for altered antibody response and/or disease resistance will not affect production.

Conclusions

We confirmed the presence of genetic variability and identified SNPs significantly associated with immune, disease and production traits in indigenous village chickens. Results underpin the feasibility of concomitant genetic improvement for enhanced antibody response, resistance to parasitism and productivity within and across indigenous chicken ecotypes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-016-0252-7) contains supplementary material, which is available to authorized users.

Dual transcriptomics of virus-host interactions: comparing two Pacific oyster families presenting contrasted susceptibility to ostreid herpesvirus 1

BACKGROUND

Massive mortality outbreaks affecting Pacific oyster (Crassostrea gigas) spat in various countries have been associated with the detection of a herpesvirus called ostreid herpesvirus type 1 (OsHV-1). However, few studies have been performed to understand and follow viral gene expression, as it has been done in vertebrate herpesviruses. In this work, experimental infection trials of C. gigas spat with OsHV-1 were conducted in order to test the susceptibility of several bi-parental oyster families to this virus and to analyze host-pathogen interactions using in vivo transcriptomic approaches.

RESULTS

The divergent response of these oyster families in terms of mortality confirmed that susceptibility to OsHV-1 infection has a significant genetic component. Two families with contrasted survival rates were selected. A total of 39 viral genes and five host genes were monitored by real-time PCR. Initial results provided information on (i) the virus cycle of OsHV-1 based on the kinetics of viral DNA replication and transcription and (ii) host defense mechanisms against the virus.

CONCLUSIONS

In the two selected families, the detected amounts of viral DNA and RNA were significantly different. This result suggests that Pacific oysters are genetically diverse in terms of their susceptibility to OsHV-1 infection. This contrasted susceptibility was associated with dissimilar host gene expression profiles. Moreover, the present study showed a positive correlation between viral DNA amounts and the level of expression of selected oyster genes.

Fine-mapping QTLs in advanced intercross lines and other outbred populations

Quantitative genetic studies in model organisms, particularly in mice, have been extremely successful in identifying chromosomal regions that are associated with a wide variety of behavioral and other traits. However, it is now widely understood that identification of the underlying genes will be far more challenging. In the last few years, a variety of populations have been utilized in an effort to more finely map these chromosomal regions with the goal of identifying specific genes. The common property of these newer populations is that linkage disequilibrium spans relatively short distances, which permits fine-scale mapping resolution. This review focuses on advanced intercross lines (AILs) which are the simplest such population. As originally proposed in 1995 by Darvasi and Soller, an AIL is the product of intercrossing two inbred strains beyond the F2 generation. Unlike recombinant inbred strains, AILs are maintained as outbred populations; brother-sister matings are specifically avoided. Each generation of intercrossing beyond the F2 further degrades linkage disequilibrium between adjacent makers, which allows for fine-scale mapping of quantitative trait loci (QTLs). Advances in genotyping technology and techniques for the statistical analysis of AILs have permitted rapid advances in the application of AILs. We review some of the analytical issues and available software, including QTLRel, EMMA, EMMAX, GEMMA, TASSEL, GRAMMAR, WOMBAT, Mendel, and others.

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.

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

A challenge test following inoculation with a standard amount of a vv+ strain of the Marek's disease (MD) virus in multiple lines and multiple generations of egg type chicken and the corresponding phenotypic trend are described. This program significantly reduced mortality of progeny from selected sires for three to 11 generations in eight of the nine elite lines studied herein. In brown egg lines, a retrospective analysis of DNA indicated an association between the blood type B (major histocompatibility complex) of the sire and the MD mortality in the challenge of its progeny. As a result of the multigeneration stock amplification and crossbreeding processes used in the commercial breeding industry, improvement in survival after challenge at the elite level will translate to improved welfare for millions of birds at the commercial production level.

Genome-wide copy number variant analysis in inbred chickens lines with different susceptibility to Marek's disease

Breeding of genetically resistant chickens to Marek’s disease (MD) is a vital strategy to poultry health. To find the markers underlying the genetic resistance to MD, copy number variation (CNV) was examined in inbred MD-resistant and -susceptible chicken lines. A total of 45 CNVs were found in four lines of chickens, and 28 were potentially involved in immune response and cell proliferation, etc. Importantly, two CNVs related with MD resistance were transmitted to descendent recombinant congenic lines that differ in susceptibility to MD. Our findings may lead to better strategies for genetic improvement of disease resistance in poultry.

Toward integrative genomics study of genetic resistance to Salmonella and Campylobacter intestinal colonization in fowl

Salmonella enterica serotypes Enteritidis and Typhimurium and Campylobacter jejuni are responsible for most cases of food poisoning in Europe. These bacteria do not cause severe disease symptoms in chicken, but they are easily propagated by symptomless chicken carriers which cannot be easily isolated. This animal tolerance is detrimental to food safety. In this particular case, increasing animal's resistance is not sufficient, since some animals considered as resistant are able to carry bacteria during several weeks without displaying disease symptoms. We review studies aimed at evaluating the resistance of chicken to Salmonella and Campylobacter intestinal colonization, either a few days or several weeks after infection. While studies of the genetic control of Campylobacter colonization are only beginning, mostly due to technical difficulties in infection protocols, genetic studies of Salmonella colonization have been conducted for now more than 20 years. They have initially reported an estimation of the genetic parameters associated with resistance to Salmonella colonization and are now aimed at identifying the genomic regions controlling variation of this trait in experimental lines and commercial populations. With the advent of high-throughput genomics, we are closer than ever to identify the true genes controlling resistance to Enterobacteria colonization in chicken. The comparison of genes involved in early resistance to intestinal colonization with genes controlling resistance to bacteria persistence several weeks after infection (i.e., carrier-state) should soon highlight the differences between the molecular mechanisms underlying those two distinct phenotypes. It will also be highly interesting to compare the genes or genomic regions controlling Campylobacter and Salmonella, in order to evaluate the feasibility of a selection conducted on both bacteria simultaneously.

Whole-genome scan for signatures of recent selection reveals loci associated with important traits in White Leghorn chickens

Chicken is considered to be an excellent model for genetic studies of phenotypic and genomic evolution, with large effective population size, specialized commercial lines, and strong human-driven selection. High-density chicken SNP chips can help to achieve a better understanding of the selection mechanisms in artificially selected populations. We performed the genome-wide tests for the selection signature in 385 White Leghorn hens and mapped positively selected regions to the genome annotations. Ten QTL related to egg production, egg quality, growth, and disease resistance traits were selected for extended haplotype homozygosity tests to give a brief overview of recent selection signatures in chicken QTL. We also reported 185 candidate genes/CDSs showing top P-values and slower decay of haplotype homozygosities. Some of these genes seemed to have significant effects on important economical traits, and most of them have not been reported in chickens. The current study provides a genome-wide map of linkage disequilibrium extents and distributions and selection footprints in the chicken genome. A panel of genes, including PRL, NCKX1, NRF1, LHX2, and SFRP1 associated with egg production, metabolism traits, and response to illumination were identified. In addition, there were more genes identified that have not yet been reported in chickens, and our results provide new clues for further study.

MHC class I of saltwater crocodiles (Crocodylus porosus): polymorphism and balancing selection

Saltwater crocodiles are in high demand for the production of luxury fashion items. However, their susceptibility to disease incurs substantial losses and it is hoped to be able to genetically select these animals for disease resistance. So far, this has only been enabled by phenotypic selection. Investigating the major histocompatibility complex (MHC) could provide insight into the ability of an individual to respond to pathogens acting as a selective pressure on the host. Here, we assessed genetic diversity and a role of selection in shaping the diversity of MHC class I exon 3 among 42 saltwater crocodiles from nine river basins in the Northern Territory, Australia. We generated 640 sequences using cloning and sequencing methods and identified 43 MHC variants among them. Phylogenetic analyses clustered these variants into two major clades, which may suggest two gene lineages. We found the number of variants within an individual varying between one and seven, indicating that there are at least four gene loci in this species. Selection detection analyses revealed an elevated ratio of nonsynonymous to synonymous substitutions (mean = 1.152 per codon), suggesting balancing selection. Population differentiation analyses revealed that the MHC did not show structuring among the river basins, and there were some shared variants among them. This may be a result of possible gene flow and/or similar selection pressures among populations. These findings provide background knowledge to identify potential MHC markers, which could be used for selecting genetically variable individuals for future disease associations. All MHC class I exon 3 sequences reported in this paper were submitted to the GenBank database with following accession numbers: HQ008785-HQ008789, HQ008791-HQ008798, HQ008808-HQ008815, HQ008824, HQ008826-HQ008830, HQ008835, HQ008839, HQ008842-HQ008850, and JX023536-JX023540.

A genome-wide SNP scan reveals two loci associated with the chicken resistance to Marek's disease

Marek's disease (MD) is a neoplastic disease in chickens, caused by the Marek's disease virus (MDV). To investigate host genetic resistance to MD, we conducted a genome-wide association study (GWAS) on 67 MDV-infected chickens based on a case and control design, including 57 susceptible chickens in the case group and 10 resistant chickens as controls. After searching 38 655 valid genomic markers, two SNPs were found to be associated with host resistance to MD. One SNP, rs14527240, reaching chromosome-wide significance level (P < 0.01) was located in the SPARC-related modular calcium-binding 1 (SMOC1) gene on GGA5. The other one, GGaluGA156129, reaching genome-wide significance (P < 0.05), was located in the protein tyrosine phosphatase, non-receptor type 3 (PTPN3) gene on GGA2. In addition, expression patterns of these two genes in spleens were detected by qPCR. The expression of SMOC1 was significantly up-regulated (P < 0.05), whereas the expression of PTNP3 did not show significance when the case group was compared with the control group. Up-regulation of SMOC1 in susceptible spleens suggests its important roles in MD tumorigenesis. This is the first study to investigate MD-resistant loci, and it demonstrates the power of GWASs for mapping genes associated with MD resistance.

Genome-Wide Identification and Quantification of cis- and trans-Regulated Genes Responding to Marek's Disease Virus Infection via Analysis of Allele-Specific Expression

Marek’s disease (MD) is a commercially important neoplastic disease of chickens caused by Marek’s disease virus (MDV), a naturally occurring oncogenic alphaherpesvirus. Selecting for increased genetic resistance to MD is a control strategy that can augment vaccinal control measures. To identify high-confidence candidate MD resistance genes, we conducted a genome-wide screen for allele-specific expression (ASE) amongst F₁ progeny of two inbred chicken lines that differ substantially in MD resistance. High throughput sequencing was initially used to profile transcriptomes from pools of uninfected and infected individuals at 4 days post-infection to identify any genes showing ASE in response to MDV infection. RNA sequencing identified 22,655 single nucleotide polymorphisms (SNPs) of which 5,360 in 3,773 genes exhibited significant allelic imbalance. Illumina GoldenGate assays were subsequently used to quantify regulatory variation controlled at the gene (cis) and elsewhere in the genome (trans) by examining differences in expression between F₁ individuals and artificial F₁ RNA pools over six time periods in 1,536 of the most significant SNPs identified by RNA sequencing. Allelic imbalance as a result of cis-regulatory changes was confirmed in 861 of the 1,233 GoldenGate assays successfully examined. Furthermore we have identified seven genes that display trans-regulation only in infected animals and ∼500 SNP that show a complex interaction between cis- and trans-regulatory changes. Our results indicate ASE analyses are a powerful approach to identify regulatory variation responsible for differences in transcript abundance in genes underlying complex traits. And the genes with SNPs exhibiting ASE provide a strong foundation to further investigate the causative polymorphisms and genetic mechanisms for MD resistance. Finally, the methods used here for identifying specific genes and SNPs have practical implications for applying marker-assisted selection to complex traits that are difficult to measure in agricultural species, when expression differences are expected to control a portion of the phenotypic variance.

Systems analysis of immune responses in Marek's disease virus-infected chickens identifies a gene involved in susceptibility and highlights a possible novel pathogenicity mechanism

Marek's disease virus (MDV) is a highly contagious oncogenic alphaherpesvirus that causes disease that is both a cancer model and a continuing threat to the world's poultry industry. This comprehensive gene expression study analyzes the host response to infection in both resistant and susceptible lines of chickens and inherent expression differences between the two lines following the infection of the host. A novel pathogenicity mechanism, involving the downregulation of genes containing HIC1 transcription factor binding sites as early as 4 days postinfection, was suggested from this analysis. HIC1 drives antitumor mechanisms, suggesting that MDV infection switches off genes involved in antitumor regulation several days before the expression of the MDV oncogene meq. The comparison of the gene expression data to previous QTL data identified several genes as candidates for involvement in resistance to MD. One of these genes, IRG1, was confirmed by single nucleotide polymorphism analysis to be involved in susceptibility. Its precise mechanism remains to be elucidated, although the analysis of gene expression data suggests it has a role in apoptosis. Understanding which genes are involved in susceptibility/resistance to MD and defining the pathological mechanisms of the disease gives us a much greater ability to try to reduce the incidence of this virus, which is costly to the poultry industry in terms of both animal welfare and economics.

Temporal transcriptome changes induced by MDV in Marek's disease-resistant and -susceptible inbred chickens

Background

Marek's disease (MD) is a lymphoproliferative disease in chickens caused by Marek's disease virus (MDV) and characterized by T cell lymphoma and infiltration of lymphoid cells into various organs such as liver, spleen, peripheral nerves and muscle. Resistance to MD and disease risk have long been thought to be influenced both by genetic and environmental factors, the combination of which contributes to the observed outcome in an individual. We hypothesize that after MDV infection, genes related to MD-resistance or -susceptibility may exhibit different trends in transcriptional activity in chicken lines having a varying degree of resistance to MD.

Results

In order to study the mechanisms of resistance and susceptibility to MD, we performed genome-wide temporal expression analysis in spleen tissues from MD-resistant line 6₃, susceptible line 7₂ and recombinant congenic strain M (RCS-M) that has a phenotype intermediate between lines 6₃ and 7₂ after MDV infection. Three time points of the MDV life cycle in chicken were selected for study: 5 days post infection (dpi), 10dpi and 21dpi, representing the early cytolytic, latent and late cytolytic stages, respectively. We observed similar gene expression profiles at the three time points in line 6₃ and RCS-M chickens that are both different from line 7₂. Pathway analysis using Ingenuity Pathway Analysis (IPA) showed that MDV can broadly influence the chickens irrespective of whether they are resistant or susceptible to MD. However, some pathways like cardiac arrhythmia and cardiovascular disease were found to be affected only in line 7₂; while some networks related to cell-mediated immune response and antigen presentation were enriched only in line 6₃ and RCS-M. We identified 78 and 30 candidate genes associated with MD resistance, at 10 and 21dpi respectively, by considering genes having the same trend of expression change after MDV infection in lines 6₃ and RCS-M. On the other hand, by considering genes with the same trend of expression change after MDV infection in lines 7₂ and RCS-M, we identified 78 and 43 genes at 10 and 21dpi, respectively, which may be associated with MD-susceptibility.

Conclusions

By testing temporal transcriptome changes using three representative chicken lines with different resistance to MD, we identified 108 candidate genes for MD-resistance and 121 candidate genes for MD-susceptibility over the three time points. Genes included in our resistance or susceptibility genes lists that are also involved in more than 5 biofunctions, such as CD8α, IL8, USP18, and CTLA4, are considered to be important genes involved in MD-resistance or -susceptibility. We were also able to identify several biofunctions related with immune response that we believe play an important role in MD-resistance.

Immune responses against Marek's disease virus

It is more than a century since Marek's disease (MD) was first reported in chickens and since then there have been concerted efforts to better understand this disease, its causative agent and various approaches for control of this disease. Recently, there have been several outbreaks of the disease in various regions, due to the evolving nature of MD virus (MDV), which necessitates the implementation of improved prophylactic approaches. It is therefore essential to better understand the interactions between chickens and the virus. The chicken immune system is directly involved in controlling the entry and the spread of the virus. It employs two distinct but interrelated mechanisms to tackle viral invasion. Innate defense mechanisms comprise secretion of soluble factors as well as cells such as macrophages and natural killer cells as the first line of defense. These innate responses provide the adaptive arm of the immune system including antibody- and cell-mediated immune responses to be tailored more specifically against MDV. In addition to the immune system, genetic and epigenetic mechanisms contribute to the outcome of MDV infection in chickens. This review discusses our current understanding of immune responses elicited against MDV and genetic factors that contribute to the nature of the response.

Comparative study of tumorigenesis and tumor immunity in invertebrates and nonmammalian vertebrates

Despite intense study in mammals, the different roles played by the immune system in detecting (immunosurveillance), controlling and remodeling (immunoediting) neoplasia, and perhaps in metastasis are not fully understood. In this review, I will present evidence of neoplasia and invasive malignancy, as well as tumor immunity in invertebrates and nonmammalian vertebrates. I will also present a comparative and evolutionary view of the complex interactions between neoplasia and the host immune system. Overall, I wish to go beyond the too simplistic dichotomy between invertebrates with innate immunity that are only affected with benign neoplasia and vertebrates with adaptive immunity that are affected by metastatic malignancies or cancer.