Genetic mapping of quantitative trait loci affecting susceptibility to Marek's disease virus induced tumors in F2 intercross chickens

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.

A Wider and Deeper Peptide-Binding Groove for the Class I Molecules from B15 Compared with B19 Chickens Correlates with Relative Resistance to Marek's Disease

The chicken MHC is known to confer decisive resistance or susceptibility to various economically important pathogens, including the iconic oncogenic herpesvirus that causes Marek's disease (MD). Only one classical class I gene, BF2, is expressed at a high level in chickens, so it was relatively easy to discern a hierarchy from well-expressed thermostable fastidious specialist alleles to promiscuous generalist alleles that are less stable and expressed less on the cell surface. The class I molecule BF2*1901 is better expressed and more thermostable than the closely related BF2*1501, but the peptide motif was not simpler as expected. In this study, we confirm for newly developed chicken lines that the chicken MHC haplotype B15 confers resistance to MD compared with B19. Using gas phase sequencing and immunopeptidomics, we find that BF2*1901 binds a greater variety of amino acids in some anchor positions than does BF2*1501. However, by x-ray crystallography, we find that the peptide-binding groove of BF2*1901 is narrower and shallower. Although the self-peptides that bound to BF2*1901 may appear more various than those of BF2*1501, the structures show that the wider and deeper peptide-binding groove of BF2*1501 allows stronger binding and thus more peptides overall, correlating with the expected hierarchies for expression level, thermostability, and MD resistance. Our study provides a reasonable explanation for greater promiscuity for BF2*1501 compared with BF2*1901, corresponding to the difference in resistance to MD.

Functional immune diversity in reindeer reveals a high Arctic population at risk

Climate changes the geographic range of both species as well as pathogens, causing a potential increase in the vulnerability of populations or species with limited genetic diversity. With advances in high throughput sequencing (HTS) technologies, we can now define functional expressed genetic diversity of wild species at a larger scale and identify populations at risk. Previous studies have used genomic DNA to define major histocompatibility complex (MHC) class II diversity in reindeer. Varying numbers of expressed genes found in many ungulates strongly argues for using cDNA in MHC typing strategies to ensure that diversity estimates relate to functional genes. We have used available reindeer genomes to identify candidate genes and established an HTS approach to define expressed MHC class I and class II diversity. To capture a broad diversity we included samples from wild reindeer from Southern Norway, semi-domesticated reindeer from Northern Norway and reindeer from the high Artic archipelago Svalbard. Our data show a medium MHC diversity in semi-domesticated and wild Norwegian mainland reindeer, and low MHC diversity reindeer in Svalbard reindeer. The low immune diversity in Svalbard reindeer provides a potential risk if the pathogenic pressure changes in response to altered environmental conditions due to climate change, or increased human-related activity.

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.

Seed protein content and its relationships with agronomic traits in pigeonpea is controlled by both main and epistatic effects QTLs

The genetic architecture of seed protein content (SPC) and its relationships to agronomic traits in pigeonpea is poorly understood. Accordingly, five F₂ populations segregating for SPC and four agronomic traits (seed weight (SW), seed yield (SY), growth habit (GH) and days to first flowering (DFF)) were phenotyped and genotyped using genotyping-by-sequencing approach. Five high-density population-specific genetic maps were constructed with an average inter-marker distance of 1.6 to 3.5 cM, and subsequently, integrated into a consensus map with average marker spacing of 1.6 cM. Based on analysis of phenotyping data and genotyping data, 192 main effect QTLs (M-QTLs) with phenotypic variation explained (PVE) of 0.7 to 91.3% were detected for the five traits across the five populations. Major effect (PVE ≥ 10%) M-QTLs included 14 M-QTLs for SPC, 16 M-QTLs for SW, 17 M-QTLs for SY, 19 M-QTLs for GH and 24 M-QTLs for DFF. Also, 573 epistatic QTLs (E-QTLs) were detected with PVE ranging from 6.3 to 99.4% across traits and populations. Colocalization of M-QTLs and E-QTLs explained the genetic basis of the significant (P < 0.05) correlations of SPC with SW, SY, DFF and GH. The nature of genetic architecture of SPC and its relationship with agronomic traits suggest that genomics-assisted breeding targeting genome-wide variations would be effective for the simultaneous improvement of SPC and other important traits.

Genetic differences in host infectivity affect disease spread and survival in epidemics

Survival during an epidemic is partly determined by host genetics. While quantitative genetic studies typically consider survival as an indicator for disease resistance (an individual's propensity to avoid becoming infected or diseased), mortality rates of populations undergoing an epidemic are also affected by endurance (the propensity of diseased individual to survive the infection) and infectivity (i.e. the propensity of an infected individual to transmit disease). Few studies have demonstrated genetic variation in disease endurance, and no study has demonstrated genetic variation in host infectivity, despite strong evidence for considerable phenotypic variation in this trait. Here we propose an experimental design and statistical models for estimating genetic diversity in all three host traits. Using an infection model in fish we provide, for the first time, direct evidence for genetic variation in host infectivity, in addition to variation in resistance and endurance. We also demonstrate how genetic differences in these three traits contribute to survival. Our results imply that animals can evolve different disease response types affecting epidemic survival rates, with important implications for understanding and controlling epidemics.

Macrophages from Susceptible and Resistant Chicken Lines have Different Transcriptomes following Marek's Disease Virus Infection

Despite successful control by vaccination, Marek’s disease (MD) has continued evolving to greater virulence over recent years. To control MD, selection and breeding of MD-resistant chickens might be a suitable option. MHC-congenic inbred chicken lines, 6₁ and 7₂, are highly resistant and susceptible to MD, respectively, but the cellular and genetic basis for these phenotypes is unknown. Marek’s disease virus (MDV) infects macrophages, B-cells, and activated T-cells in vivo. This study investigates the cellular basis of resistance to MD in vitro with the hypothesis that resistance is determined by cells active during the innate immune response. Chicken bone marrow-derived macrophages from lines 6₁ and 7₂ were infected with MDV in vitro. Flow cytometry showed that a higher percentage of macrophages were infected in line 7₂ than in line 6₁. A transcriptomic study followed by in silico functional analysis of differentially expressed genes was then carried out between the two lines pre- and post-infection. Analysis supports the hypothesis that macrophages from susceptible and resistant chicken lines display a marked difference in their transcriptome following MDV infection. Resistance to infection, differential activation of biological pathways, and suppression of oncogenic potential are among host defense strategies identified in macrophages from resistant chickens.

Genome-wide association analysis of intra-specific QTL associated with the resistance for enteric septicemia of catfish

Disease resistance is one of the most important traits for aquaculture industry. For catfish industry, enteric septicemia of catfish (ESC), caused by the bacterial pathogen Edwardsiella ictaluri, is the most severe disease, causing enormous economic losses every year. In this study, we used three channel catfish families with 900 individuals (300 fish per family) and the 690K catfish SNP array, and conducted a genome-wide association study to detect the quantitative trait loci (QTL) associated with ESC resistance. Three significant QTL, with two of located on LG1 and one on LG26, and three suggestive QTL located on LG1, LG3, and LG21, respectively, were identified to be associated with ESC resistance. With a well-assembled- and -annotated reference genome sequence, genes around the involved QTL regions were identified. Among these genes, 37 genes had known functions in immunity, which may be involved in ESC resistance. Notably, nlrc3 and nlrp12 identified here were also found in QTL regions of ESC resistance in the channel catfish × blue catfish interspecific hybrid system, suggesting this QTL was operating within both intra-specific channel catfish populations and interspecific hybrid backcross populations. Many of the genes of the Class I MHC pathway, for mediated antigen processing and presentation, were found in the QTL regions. The positional correlation found in this study and the expressional correlation found in previous studies indicated that Class I MHC pathway was significantly associated with ESC resistance. This study validated one QTL previously identified using the second and fourth generation of the interspecific hybrid backcross progenies, and identified five additional QTL among channel catfish families. Taken together, it appears that there are only a few major QTL for ESC disease resistance, making marker-assisted selection an effective approach for genetic improvements of ESC resistance.

Relationships between genetic change and infectious disease in domestic livestock

The relations between genetic change in domestic livestock and infectious disease (including both its epidemiology and the animal's reaction to it) are examined. The overall picture is confusing because there are different, and seemingly unrelated, ways of considering the issue. An attempt is made to put these together into a more general framework. Four processes of particular interest are distinguished and discussed in more detail: (i) the way a population's genetic potential for immunocompetence can be changed by breeding, (ii) the way an animal's immunocompetence is influenced by that animal's production potential, in combination with the environmental resources that are available to it at a given time, (iii) the way the disease status of an animal (and a population of animals) is influenced by its immunocompetence, and (iv) the way the production level of an animal is influenced by activation of its immune system. Ultimately, all four processes influence the realised level of production.

This comes down to four questions that need to be addressed: (i) can we use genetic variation in immunocompetence in animal breeding? (ii) does a higher production potential (today's direction of breeding) have a negative impact on immunocompetence? (iii) does improved immunocompetence result in improved health? and (iv) how large is the negative impact of disease on production?

Bursal transcriptome profiling of different inbred chicken lines reveals key differentially expressed genes at 3 days post-infection with very virulent infectious bursal disease virus

Infectious bursal disease is a highly contagious disease in the poultry industry and causes immunosuppression in chickens. Genome-wide regulations of immune response genes of inbred chickens with different genetic backgrounds, following very virulent infectious bursal disease virus (vvIBDV) infection are poorly characterized. Therefore, this study aims to analyse the bursal tissue transcriptome of six inbred chicken lines 6, 7, 15, N, O and P following infection with vvIBDV strain UK661 using strand-specific next-generation sequencing, by highlighting important genes and pathways involved in the infected chicken during peak infection at 3 days post-infection. All infected chickens succumbed to the infection without major variations among the different lines. However, based on the viral loads and bursal lesion scoring, lines P and 6 can be considered as the most susceptible lines, while lines 15 and N were regarded as the least affected lines. Transcriptome profiling of the bursa identified 4588 genes to be differentially expressed, with 2985 upregulated and 1642 downregulated genes, in which these genes were commonly or uniquely detected in all or several infected lines. Genes that were upregulated are primarily pro-inflammatory cytokines, chemokines and IFN-related. Various genes that are associated with B-cell functions and genes related to apoptosis were downregulated, together with the genes involved in p53 signalling. In conclusion, bursal transcriptome profiles of different inbred lines showed differential expressions of pro-inflammatory cytokines and chemokines, Th1 cytokines, JAK-STAT signalling genes, MAPK signalling genes, and their related pathways following vvIBDV infection.

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.

Validation of Alternative Transcript Splicing in Chicken Lines that Differ in Genetic Resistance to Marek's Disease

Utilizing RNA-seq data, 1,574 candidate genes with alternative splicing were previously identified between two chicken lines that differ in Marek's disease (MD) genetic resistance under control and Marek's disease virus infection conditions. After filtering out 1,530 genes with splice variants in the first or last exon, 44 genes were screened for possible exon loss or gain using PCR and gel electrophoresis. Consequently, 7 genes exhibited visually detectable differential expression of splice variants between lines 6 (MD resistant) and 7 (MD susceptible), and the resultant PCR products verified by DNA sequencing. Birds from inbred line 6 have transcripts that preferentially retain an exon compared to line 7 chickens for ITGB2, SGPL1, and COMMD5. Birds from inbred line 7 have alleles that preferentially retain an exon compared to line 6 for MOCS2. CCBL2 exon 1a is absent and ATAD1 exon 2 is truncated by 87 nucleotides in transcripts expressed by line 7 compared to those from line 6. For CHTF18, line 6 transcripts have an indel mutation with 7 additional nucleotides in exon 21 compared to line 7. The current study validates 7 genes with alternatively spliced isomers between the two chicken lines, which helps provide potential underlying mechanisms for the phenotypic differences.

Virus and host genomic, molecular, and cellular interactions during Marek's disease pathogenesis and oncogenesis

Marek's Disease Virus (MDV) is a chicken alphaherpesvirus that causes paralysis, chronic wasting, blindness, and fatal lymphoma development in infected, susceptible host birds. This disease and its protective vaccines are highly relevant research targets, given their enormous impact within the poultry industry. Further, Marek's disease (MD) serves as a valuable model for the investigation of oncogenic viruses and herpesvirus patterns of viral latency and persistence--as pertinent to human health as to poultry health. The objectives of this article are to review MDV interactions with its host from a variety of genomic, molecular, and cellular perspectives. In particular, we focus on cytogenetic studies, which precisely assess the physical status of the MDV genome in the context of the chicken host genome. Combined, the cytogenetic and genomic research indicates that MDV-host genome interactions, specifically integration of the virus into the host telomeres, is a key feature of the virus life cycle, contributing to the viral achievement of latency, transformation, and reactivation of lytic replication. We present a model that outlines the variety of virus-host interactions, at the multiple levels, and with regard to the disease states.

A Novel Statistical Model to Estimate Host Genetic Effects Affecting Disease Transmission

There is increasing recognition that genetic diversity can affect the spread of diseases, potentially affecting plant and livestock disease control as well as the emergence of human disease outbreaks. Nevertheless, even though computational tools can guide the control of infectious diseases, few epidemiological models can simultaneously accommodate the inherent individual heterogeneity in multiple infectious disease traits influencing disease transmission, such as the frequently modeled propensity to become infected and infectivity, which describes the host ability to transmit the infection to susceptible individuals. Furthermore, current quantitative genetic models fail to fully capture the heritable variation in host infectivity, mainly because they cannot accommodate the nonlinear infection dynamics underlying epidemiological data. We present in this article a novel statistical model and an inference method to estimate genetic parameters associated with both host susceptibility and infectivity. Our methodology combines quantitative genetic models of social interactions with stochastic processes to model the random, nonlinear, and dynamic nature of infections and uses adaptive Bayesian computational techniques to estimate the model parameters. Results using simulated epidemic data show that our model can accurately estimate heritabilities and genetic risks not only of susceptibility but also of infectivity, therefore exploring a trait whose heritable variation is currently ignored in disease genetics and can greatly influence the spread of infectious diseases. Our proposed methodology offers potential impacts in areas such as livestock disease control through selective breeding and also in predicting and controlling the emergence of disease outbreaks in human populations.

Genome-wide identification of copy number variations between two chicken lines that differ in genetic resistance to Marek's disease

Background

Copy number variation (CNV) is a major source of genome polymorphism that directly contributes to phenotypic variation such as resistance to infectious diseases. Lines 6₃ and 7₂ are two highly inbred experimental chicken lines that differ greatly in susceptibility to Marek’s disease (MD), and have been used extensively in efforts to identify the genetic and molecular basis for genetic resistance to MD. Using next generation sequencing, we present a genome-wide assessment of CNVs that are potentially associated with genetic resistance to MD.

Methods

Three chickens randomly selected from each line were sequenced to an average depth of 20×. Two popular software, CNVnator and Pindel, were used to call genomic CNVs separately. The results were combined to obtain a union set of genomic CNVs in the two chicken lines.

Results

A total of 5,680 CNV regions (CNVRs) were identified after merging the two datasets, of which 1,546 and 1,866 were specific to the MD resistant or susceptible line, respectively. Over half of the line-specific CNVRs were shared by 2 or more chickens, reflecting the reduced diversity in both inbred lines. The CNVRs fixed in the susceptible lines were significantly enriched in genes involved in MAPK signaling pathway. We also found 67 CNVRs overlapping with 62 genes previously shown to be strong candidates of the underlying genes responsible for the susceptibility to MD.

Conclusions

Our findings provide new insights into the genetic architecture of the two chicken lines and additional evidence that MAPK signaling pathway may play an important role in host response to MD virus infection. The rich source of line-specific CNVs is valuable for future disease-related association studies in the two chicken lines.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2080-5) contains supplementary material, which is available to authorized users.

Fine mapping of QTL and genomic prediction using allele-specific expression SNPs demonstrates that the complex trait of genetic resistance to Marek's disease is predominantly determined by transcriptional regulation

Background

Marek’s disease (MD) is a lymphoproliferative disease of poultry induced by Marek’s disease virus (MDV), a highly oncogenic alphaherpesvirus. Identifying the underlying genes conferring MD genetic resistance is desired for more efficacious control measures including genomic selection, which requires accurately identified genetic markers throughout the chicken genome.

Methods

Hypothesizing that variants located in transcriptional regulatory regions are the main mechanism underlying this complex trait, a genome-wide association study was conducted by genotyping a ~1,000 bird MD resource population derived from experimental inbred layers with SNPs containing 1,824 previously identified allele-specific expression (ASE) SNPs in response to MDV infection as well as 3,097 random SNPs equally spaced throughout the chicken genome. Based on the calculated associations, genomic predictions were determined for 200 roosters and selected sires had their progeny tested for Marek’s disease incidence.

Results

Our analyses indicate that these ASE SNPs account for more than 83 % of the genetic variance and exhibit nearly all the highest associations. To validate these findings, 200 roosters had their genetic merit predicted from the ASE SNPs only, and the top 30 and bottom 30 ranked roosters were reciprocally mated to random hens. The resulting progeny showed that after only one generation of bidirectional selection, there was a 22 % difference in MD incidence and this approach gave a 125 % increase in accuracy compared to current pedigree-based estimates.

Conclusions

We conclude that variation in transcriptional regulation is the major driving cause for genetic resistance to MD, and ASE SNPs identify the underlying genes and are sufficiently linked to the causative polymorphisms that they can be used for accurate genomic prediction as well as help define the underlying molecular basis. Furthermore, this approach should be applicable to other complex traits.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2016-0) contains supplementary material, which is available to authorized users.

Expression levels of MHC class I molecules are inversely correlated with promiscuity of peptide binding

Highly polymorphic major histocompatibility complex (MHC) molecules are at the heart of adaptive immune responses, playing crucial roles in many kinds of disease and in vaccination. We report that breadth of peptide presentation and level of cell surface expression of class I molecules are inversely correlated in both chickens and humans. This relationship correlates with protective responses against infectious pathogens including Marek's disease virus leading to lethal tumours in chickens and human immunodeficiency virus infection progressing to AIDS in humans. We propose that differences in peptide binding repertoire define two groups of MHC class I molecules strategically evolved as generalists and specialists for different modes of pathogen resistance. We suggest that differences in cell surface expression level ensure the development of optimal peripheral T cell responses. The inverse relationship of peptide repertoire and expression is evidently a fundamental property of MHC molecules, with ramifications extending beyond immunology and medicine to evolutionary biology and conservation.

Establishment of six homozygous MHC-B haplotype populations associated with susceptibility to Marek's disease in Chinese specific pathogen-free BWEL chickens

The highly polymorphic chicken major histocompatibility complex (MHC) is associated with different levels of immunologic responses to certain avian pathogens. MHC-B haplotype chickens are an important genetic resource for studying the genetic determination of pathogen resistance and susceptibility. The BWEL chicken population is the only specific pathogen-free (SPF) chickens bred and developed by the State Center of Poultry Genetic Resources of Laboratory Animals in China. In this study, we successfully established six homozygous MHC-B haplotype populations from the BWEL chickens using microsatellite marker technology, named as BW/G(1, 2, 3, 5, 6, 7) lines, and their molecular genotypes were matched to six serologically defined MHC-B haplotypes, B13, B15, B2, B5, B21 and B19, respectively. The sequences of BF genes exons 2 and 3 from four successive generations (F1-F4) of the BW/G(n) lines were completely consistent with those of serologically defined MHC-B haplotypes. Subsequently, six BW/G(n) line specific allo-antisera were prepared by immunization with red blood cells (RBCs) and hemagglutination tests results showed the BW/G(n) SPF chickens could be serologically differentiated. Additionally, susceptibility to Marek's disease (MD) in the BW/G3 (B2 haplotype) and BW/G7 (B19 haplotype) lines were determined by comparing mortality, macroscopic and histopathological lesions, and viral loads in feather pulp. The BW/G7 line showed greater genetic susceptibility to the very virulent MD virus (MDV) strain than the BW/G3 line. The establishment of MHC-B haplotype chicken populations associated with susceptibility to MD will be helpful for studying host immune responses and further developing the more effective vaccines in the context of MHC specificities, and they are also very useful for an understanding of MHC genes architecture and function.

Genome-wide patterns of genetic variation in two domestic chickens

Domestic chickens are excellent models for investigating the genetic basis of phenotypic diversity, as numerous phenotypic changes in physiology, morphology, and behavior in chickens have been artificially selected. Genomic study is required to study genome-wide patterns of DNA variation for dissecting the genetic basis of phenotypic traits. We sequenced the genomes of the Silkie and the Taiwanese native chicken L2 at ∼23- and 25-fold average coverage depth, respectively, using Illumina sequencing. The reads were mapped onto the chicken reference genome (including 5.1% Ns) to 92.32% genome coverage for the two breeds. Using a stringent filter, we identified ∼7.6 million single-nucleotide polymorphisms (SNPs) and 8,839 copy number variations (CNVs) in the mapped regions; 42% of the SNPs have not found in other chickens before. Among the 68,906 SNPs annotated in the chicken sequence assembly, 27,852 were nonsynonymous SNPs located in 13,537 genes. We also identified hundreds of shared and divergent structural and copy number variants in intronic and intergenic regions and in coding regions in the two breeds. Functional enrichments of identified genetic variants were discussed. Radical nsSNP-containing immunity genes were enriched in the QTL regions associated with some economic traits for both breeds. Moreover, genetic changes involved in selective sweeps were detected. From the selective sweeps identified in our two breeds, several genes associated with growth, appetite, and metabolic regulation were identified. Our study provides a framework for genetic and genomic research of domestic chickens and facilitates the domestic chicken as an avian model for genomic, biomedical, and evolutionary studies.

Applied animal genomics: results from the field

Genomic selection (GS) is the use of statistical methods to estimate the genetic merit of a genotyped animal based on prediction equations derived from large ancestral populations with both phenotypes and genotypes. It has revolutionized the dairy cattle breeding industry and has been implemented with varying degrees of success in other animal breeding programs, including swine, poultry, and beef cattle. The findings of empirical field studies applying GS to the breeding sectors of these main animal protein industries are reviewed. Several translational considerations must be addressed before implementing GS in genetic improvement programs. These include determining and obtaining economically relevant phenotypes and determining the optimal size of the training population, cost-effective genotyping strategies, the practicality of field implementation, and the relative costs versus the benefits of the realized rate of genetic gain. GS may additionally change the optimal breeding scheme design, and studies that address this consideration are also reviewed briefly.

Transcriptional profiling of mEq-dependent genes in Marek's disease resistant and susceptible inbred chicken lines

Marek’s disease (MD) is an economically significant disease in chickens caused by the highly oncogenic Marek’s disease virus (MDV). Understanding the genes and biological pathways that confer MD genetic resistance should lead towards the development of more disease resistant commercial poultry flocks or improved MD vaccines. MDV mEq, a bZIP transcription factor, is largely attributed to viral oncogenicity though only a few host target genes have been described, which has impeded our understanding of MDV-induced tumorigenesis. Given the importance of mEq in MDV-induced pathogenesis, we explored the role of mEq in genetic resistance to MDV. Using global transcriptome analysis and cells from MD resistant or susceptible birds, we compared the response to infection with either wild type MDV or a nononcogenic recombinant lacking mEq. As a result, we identified a number of specific genes and pathways associated with either MD resistance or susceptibility. Additionally, integrating prior information from ChIP-seq, microarray analysis, and SNPs exhibiting allele-specific expression (ASE) in response to MDV infection, we were able to provide evidence for 24 genes that are polymorphic within mEq binding sites are likely to account for gene expression in an allele-specific manner and potentially for the underlying genetic differences in MD incidence.

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.

Expression pattern of genes of RLR-mediated antiviral pathway in different-breed chicken response to Marek's disease virus infection

It has been known that the chicken's resistance to disease was affected by chicken's genetic background. And RLR-mediated antiviral pathway plays an important role in detection of viral RNA. However, little is known about the interaction of genetic background with RLR-mediated antiviral pathway in chicken against MDV infection. In this study, we adopted economic line-AA broilers and native Erlang mountainous chickens for being infected with MDV. Upon infection with MDV, the expression of MDA-5 was upregulated in two-breed chickens at 4, 7, and 21 d.p.i. It is indicated that MDA-5 might be involved in detecting MDV in chicken. Interestingly, the expression of IRF-3 and IFN-β genes was decreased in spleen and thymus of broilers at 21 d.p.i, but it was upregulated in immune tissues of Erlang mountainous chickens. And the genome load of MDV in spleen of broiler is significantly higher than that in Erlang mountainous chickens. Meanwhile, we observed that the death of broiler mainly also occurred in this phase. Collectively, these present results demonstrated that the expression patters of IRF-3 and IFN-β genes in chicken against MDV infection might be affected by the genetic background which sequently influence the resistance of chicken response to MDV.

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.

Comparison and contrast of genes and biological pathways responding to Marek's disease virus infection using allele-specific expression and differential expression in broiler and layer chickens

BACKGROUND

Marek's disease (MD) is a commercially important neoplastic disease of chickens caused by the Marek's disease virus (MDV), a naturally occurring oncogenic alphaherpesvirus. Enhancing MD genetic resistance is desirable to augment current vaccines and other MD control measures. High throughput sequencing was used to profile splenic transcriptomes from individual F1 progeny infected with MDV at 4 days of age from both outbred broilers (meat-type) and inbred layer (egg-type) chicken lines that differed in MD genetic resistance. The resulting information was used to identify SNPs, genes, and biological pathways exhibiting allele-specific expression (ASE) in response to MDV infection in each type of chicken. In addition, we compared and contrasted the results of pathway analyses (ASE and differential expression (DE)) between chicken types to help inform on the biological response to MDV infection.

RESULTS

With 7 individuals per line and treatment group providing high power, we identified 6,132 single nucleotide polymorphisms (SNPs) in 4,768 genes and 4,528 SNPs in 3,718 genes in broilers and layers, respectively, that exhibited ASE in response to MDV infection. Furthermore, 548 and 434 genes in broilers and layers, respectively, were found to show DE following MDV infection. Comparing the datasets, only 72 SNPs and 850 genes for ASE and 20 genes for DE were common between the two bird types. Although the chicken types used in this study were genetically different, at the pathway level, both TLR receptor and JAK/STAT signaling pathways were enriched as well as exhibiting a high proportion of ASE genes, especially at the beginning of both above mentioned regulatory pathways.

CONCLUSIONS

RNA sequencing with adequate biological replicates is a powerful approach to identify high confidence SNPs, genes, and pathways that are associated with transcriptional response to MDV infection. In addition, the SNPs exhibiting ASE in response to MDV infection provide a strong foundation for determining the extent to which variation in expression influences MD incidence plus yield genetic markers for genomic selection. However, given the paucity of overlap among ASE SNP sets (broilers vs. layers), it is likely that separate screens need to be incorporated for each population. Finally, comparison of gene lists obtained between these two diverse chicken types indicate the TLR and JAK/STAT signaling are conserved when responding to MDV infection and may be altered by selection of genes exhibiting ASE found at the start of each pathway.

An assessment of opportunities to dissect host genetic variation in resistance to infectious diseases in livestock

This paper reviews the evidence for host genetic variation in resistance to infectious diseases for a wide variety of diseases of economic importance in poultry, cattle, pig, sheep and Atlantic salmon. Further, it develops a method of ranking each disease in terms of its overall impact, and combines this ranking with published evidence for host genetic variation and information on the current state of genomic tools in each host species. The outcome is an overall ranking of the amenability of each disease to genomic studies that dissect host genetic variation in resistance. Six disease-based assessment criteria were defined: industry concern, economic impact, public concern, threat to food safety or zoonotic potential, impact on animal welfare and threat to international trade barriers. For each category, a subjective score was assigned to each disease according to the relative strength of evidence, impact, concern or threat posed by that particular disease, and the scores were summed across categories. Evidence for host genetic variation in resistance was determined from available published data, including breed comparison, heritability studies, quantitative trait loci (QTL) studies, evidence of candidate genes with significant effects, data on pathogen sequence and on host gene expression analyses. In total, 16 poultry diseases, 13 cattle diseases, nine pig diseases, 11 sheep diseases and three Atlantic salmon diseases were assessed. The top-ranking diseases or pathogens, i.e. those most amenable to studies dissecting host genetic variation, were Salmonella in poultry, bovine mastitis, Marek's disease and coccidiosis, both in poultry. The top-ranking diseases or pathogens in pigs, sheep and Atlantic salmon were Escherichia coli, mastitis and infectious pancreatic necrosis, respectively. These rankings summarise the current state of knowledge for each disease and broadly, although not entirely, reflect current international research efforts. They will alter as more information becomes available and as genome tools become more sophisticated for each species. It is suggested that this approach could be used to rank diseases from other perspectives as well, e.g. in terms of disease control strategies.

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.

Genetics and vaccine efficacy: host genetic variation affecting Marek's disease vaccine efficacy in White Leghorn chickens

Marek's disease (MD) is a T-cell lymphoma disease of domestic chickens induced by MD virus (MDV), a naturally oncogenic and highly contagious cell-associated α-herpesvirus. Earlier reports have shown that the MHC haplotype as well as non-MHC genes are responsible for genetic resistance to MD. The MHC was also shown to affect efficiency of vaccine response. Using specific-pathogen-free chickens from a series of 19 recombinant congenic strains and their 2 progenitor lines (lines 6(3) and 7(2)), vaccine challenge experiments were conducted to examine the effect of host genetic variation on vaccine efficacy. The 21 inbred lines of White Leghorns share the same B*2 MHC haplotype and the genome of each recombinant congenic strain differs by a random 1/8 sample of the susceptible donor line (7(2)) genome. Chickens from each of the lines were divided into 2 groups. One was vaccinated with turkey herpesvirus strain FC126 at the day of hatch and the other was treated as a nonvaccinated control. Chickens of both groups were inoculated with a very virulent plus strain of MDV on the fifth day posthatch. Analyses of the MD data showed that the genetic line significantly influenced MD incidence and days of survival post-MDV infection after vaccination of chickens (P<0.01). The protective indices against MD varied greatly among the lines with a range of 0 up to 84%. This is the first evidence that non-MHC host genetic variation significantly affects MD vaccine efficacy in chickens in a designed prospective study.

Evidence of major genes affecting resistance to bacterial cold water disease in rainbow trout using Bayesian methods of segregation analysis

Bacterial cold water disease (BCWD) causes significant economic loss in salmonid aquaculture. We previously detected genetic variation for BCWD resistance in our rainbow trout population, and a family-based selection program to improve resistance was initiated at the National Center for Cool and Cold Water Aquaculture (NCCCWA). This study investigated evidence of major trait loci affecting BCWD resistance using only phenotypic data (without using genetic markers) and Bayesian methods of segregation analysis (BMSA). A total of 10,603 juvenile fish from 101 full-sib families corresponding to 3 generations (2005, 2007, and 2009 hatch years) of the NCCCWA population were challenged by intraperitoneal injection with Flavobacterium psychrophilum, the bacterium that causes BCWD. The results from single- and multiple-QTL models of BMSA suggest that 6 to 10 QTL explaining 83 to 89% of phenotypic variance with either codominant or dominant disease-resistant alleles plus polygenic effects may underlie the genetic architecture of BCWD resistance. This study also highlights the importance of polygenic background effects in the genetic variation of BCWD resistance. The polygenic heritability on the observed scale of survival status is slightly larger than that previously reported for rainbow trout BCWD resistance. These findings provide the basis for designing informative crosses for QTL mapping and carrying out genome scans for QTL affecting BCWD resistance in rainbow trout.

Mapping quantitative trait loci affecting susceptibility to Marek's disease virus in a backcross population of layer chickens

Marek's disease (MD), caused by the oncogenic MD avian herpes virus (MDV), is a major source of economic losses to the poultry industry. A reciprocal backcross (BC) population (total 2052 individuals) was generated by crossing two partially inbred commercial Leghorn layer lines known to differ in MDV resistance, measured as survival time after challenge with a (vv+) MDV. QTL affecting resistance were identified by selective DNA pooling using a panel of 198 microsatellite markers covering two-thirds of the chicken genome. Data for each BC were analyzed separately, and as a combined data set. Markers showing significant association with resistance generally appeared in blocks of two or three, separated by blocks of nonsignificant markers. Defined this way, 15 chromosomal regions (QTLR) affecting MDV resistance, distributed among 10 chromosomes (GGA 1, 2, 3, 4, 5, 7, 8, 9, 15, and Z), were identified. The identified QTLR include one gene and three QTL associated with resistance in previous studies of other lines, and three additional QTL associated with resistance in previous studies of the present lines. These QTL could be used in marker-assisted selection (MAS) programs for MDV resistance and as a platform for high-resolution mapping and positional cloning of the resistance genes.

Functional genomics of the chicken--a model organism

Since the sequencing of the genome and the development of high-throughput tools for the exploration of functional elements of the genome, the chicken has reached model organism status. Functional genomics focuses on understanding the function and regulation of genes and gene products on a global or genome-wide scale. Systems biology attempts to integrate functional information derived from multiple high-content data sets into a holistic view of all biological processes within a cell or organism. Generation of a large collection ( approximately 600K) of chicken expressed sequence tags, representing most tissues and developmental stages, has enabled the construction of high-density microarrays for transcriptional profiling. Comprehensive analysis of this large expressed sequence tag collection and a set of approximately 20K full-length cDNA sequences indicate that the transcriptome of the chicken represents approximately 20,000 genes. Furthermore, comparative analyses of these sequences have facilitated functional annotation of the genome and the creation of several bioinformatic resources for the chicken. Recently, about 20 papers have been published on transcriptional profiling with DNA microarrays in chicken tissues under various conditions. Proteomics is another powerful high-throughput tool currently used for examining the dynamics of protein expression in chicken tissues and fluids. Computational analyses of the chicken genome are providing new insight into the evolution of gene families in birds and other organisms. Abundant functional genomic resources now support large-scale analyses in the chicken and will facilitate identification of transcriptional mechanisms, gene networks, and metabolic or regulatory pathways that will ultimately determine the phenotype of the bird. New technologies such as marker-assisted selection, transgenics, and RNA interference offer the opportunity to modify the phenotype of the chicken to fit defined production goals. This review focuses on functional genomics in the chicken and provides a road map for large-scale exploration of the chicken genome.

Livestock genomics: bridging the gap between mice and men

Dissecting the genetic control of variation in complex traits, such as disease resistance and agricultural-product quality, remains very challenging. Farm animals are now well placed to bridge the gap between human biology and traditional model species. Livestock species share with model species the benefits of controlled breeding, and their biology is often much closer to that of humans. Genetic research in model species focuses on differences between homogenous lines, whereas genetic research in humans focuses on genetic variation within populations. Livestock genetics has the strengths of both human and model-species genetics because researchers can exploit both the abundant genetic variation between divergent breeds and the variation that is segregating within breeds. Therefore, livestock genomics fills the void where the genetics of model species proves intractable or where model species are not a good proxy for human biology.