Fine-Mapping of Coccidia-Resistant Quantitative Trait Loci in Chickens

Genomic heterozygosity is associated with parasite abundance, but the effects are not mediated by host condition

Whether, when, and how genetic diversity buffers individuals and populations against infectious disease risk is a critical and open question for understanding wildlife disease and zoonotic disease risk. Several, but not all, studies have found negative relationships between infection and heterozygosity in wildlife. Since they can host multiple zoonotic infections, we sampled a population of wild deer mice (Peromyscus maniculatus), sequenced their genomes, and examined their fecal samples for coccidia and nematode eggs. We analyzed coccidia infection status, abundance, and coinfection status in relation to per-locus and per-individual measures of heterozygosity, as well as identified SNPs associated with infection status. Since heterozygosity might affect host condition, and condition is known to affect immunity, it was included as a co-variate in the per-individual analyses and as response variable in relation to heterozygosity. Not only did coccidia-infected individuals have lower levels of genome-wide per-locus diversity across all metrics, but we found an inverse relationship between genomic diversity and severity of coccidia infection. We also found weaker evidence that coinfected individuals had lower levels of private allelic variation than all other groups. In the per-individual analyses, relationships between heterozygosity and infection were marginal but followed the same negative trends. Condition was negatively correlated with infection, but was not associated with heterozygosity, suggesting that effects of heterozygosity on infection were not mediated by host condition in this system. Association tests identified multiple loci involved in the inflammatory response, with a particular role for NF-κB signaling, supporting previous work on the genetic basis of coccidia resistance. Taken together, we find that increased genome-wide neutral diversity, the presence of specific genetic variants, and improved condition positively impact infection status. Our results underscore the importance of considering host genomic variation as a buffer against infection, especially in systems that can harbor zoonotic diseases.

Supplementary Information

The online version contains supplementary material available at 10.1007/s10682-022-10175-8.

Forward Genetics in Apicomplexa Biology: The Host Side of the Story

Forward genetic approaches have been widely used in parasitology and have proven their power to reveal the complexities of host-parasite interactions in an unbiased fashion. Many aspects of the parasite’s biology, including the identification of virulence factors, replication determinants, antibiotic resistance genes, and other factors required for parasitic life, have been discovered using such strategies. Forward genetic approaches have also been employed to understand host resistance mechanisms to parasitic infection. Here, we will introduce and review all forward genetic approaches that have been used to identify host factors involved with Apicomplexa infections, which include classical genetic screens and QTL mapping, GWAS, ENU mutagenesis, overexpression, RNAi and CRISPR-Cas9 library screens. Collectively, these screens have improved our understanding of host resistance mechanisms, immune regulation, vaccine and drug designs for Apicomplexa parasites. We will also discuss how recent advances in molecular genetics give present opportunities to further explore host-parasite relationships.

Dissecting the Genomic Architecture of Resistance to Eimeria maxima Parasitism in the Chicken

Coccidiosis in poultry, caused by protozoan parasites of the genus Eimeria, is an intestinal disease with substantial economic impact. With the use of anticoccidial drugs under public and political pressure, and the comparatively higher cost of live-attenuated vaccines, an attractive complementary strategy for control is to breed chickens with increased resistance to Eimeria parasitism. Prior infection with Eimeria maxima leads to complete immunity against challenge with homologous strains, but only partial resistance to challenge with antigenically diverse heterologous strains. We investigate the genetic architecture of avian resistance to E. maxima primary infection and heterologous strain secondary challenge using White Leghorn populations of derived inbred lines, C.B12 and 15I, known to differ in susceptibility to the parasite. An intercross population was infected with E. maxima Houghton (H) strain, followed 3 weeks later by E. maxima Weybridge (W) strain challenge, while a backcross population received a single E. maxima W infection. The phenotypes measured were parasite replication (counting fecal oocyst output or qPCR for parasite numbers in intestinal tissue), intestinal lesion score (gross pathology, scale 0-4), and for the backcross only, serum interleukin-10 (IL-10) levels. Birds were genotyped using a high density genome-wide DNA array (600K, Affymetrix). Genome-wide association study located associations on chromosomes 1, 2, 3, and 5 following primary infection in the backcross population, and a suggestive association on chromosome 1 following heterologous E. maxima W challenge in the intercross population. This mapped several megabases away from the quantitative trait locus (QTL) linked to the backcross primary W strain infection, suggesting different underlying mechanisms for the primary- and heterologous secondary- responses. Underlying pathways for those genes located in the respective QTL for resistance to primary infection and protection against heterologous challenge were related mainly to immune response, with IL-10 signaling in the backcross primary infection being the most significant. Additionally, the identified markers associated with IL-10 levels exhibited significant additive genetic variance. We suggest this is a phenotype of interest to the outcome of challenge, being scalable in live birds and negating the requirement for single-bird cages, fecal oocyst counts, or slaughter for sampling (qPCR).

Application of omics technologies for a deeper insight into quali-quantitative production traits in broiler chickens: A review

The poultry industry is continuously facing substantial and different challenges such as the increasing cost of feed ingredients, the European Union's ban of antibiotic as growth promoters, the antimicrobial resistance and the high incidence of muscle myopathies and breast meat abnormalities. In the last decade, there has been an extraordinary development of many genomic techniques able to describe global variation of genes, proteins and metabolites expression level. Proper application of these cutting-edge omics technologies (mainly transcriptomics, proteomics and metabolomics) paves the possibility to understand much useful information about the biological processes and pathways behind different complex traits of chickens. The current review aimed to highlight some important knowledge achieved through the application of omics technologies and proteo-genomics data in the field of feed efficiency, nutrition, meat quality and disease resistance in broiler chickens.

The Evolutionary Biology, Ecology and Epidemiology of Coccidia of Passerine Birds

Coccidia are intracellular parasites of the phylum Apicomplexa that cause a range of pathologies collectively termed coccidiosis. Species of coccidia of commercial importance have been well studied, with the effect of other species on passerine birds receiving increasing attention. In this chapter, we review the literature on coccidia in passerines, with a particular focus on wild populations. The taxonomy and life cycle of passerine coccidia are covered, as is their impact on the health of passerines, their epidemiology and their role in parasite-mediated natural and sexual selection. Coccidia can pose a significant threat to the health of wild passerine populations, and high rates of mortality have been observed in some studies. We examine some of the genetic factors that influence host resistance to coccidia and discuss how these parasites may be important in relation to sexually selected traits. General patterns are beginning to emerge with regard to the epidemiology of the parasites, and the influence of different aspects of the host's ecology on the prevalence and intensity of coccidia is being revealed. We examine these, as well exceptions, in addition to the phenomenon of diurnal oocyst shedding that can bias studies if not accounted for. Finally, we discuss potential future directions for research on coccidia in passerines and the importance of understanding parasite ecology in the management of threatened species.

Genome-wide association study and biological pathway analysis of the Eimeria maxima response in broilers

Background

Coccidiosis is the most common and costly disease in the poultry industry and is caused by protozoans of the Eimeria genus. The current control of coccidiosis, based on the use of anticoccidial drugs and vaccination, faces serious obstacles such as drug resistance and the high costs for the development of efficient vaccines, respectively. Therefore, the current control programs must be expanded with complementary approaches such as the use of genetics to improve the host response to Eimeria infections. Recently, we have performed a large-scale challenge study on Cobb500 broilers using E. maxima for which we investigated variability among animals in response to the challenge. As a follow-up to this challenge study, we performed a genome-wide association study (GWAS) to identify genomic regions underlying variability of the measured traits in the response to Eimeria maxima in broilers. Furthermore, we conducted a post-GWAS functional analysis to increase our biological understanding of the underlying response to Eimeria maxima challenge.

Results

In total, we identified 22 single nucleotide polymorphisms (SNPs) with q value <0.1 distributed across five chromosomes. The highly significant SNPs were associated with body weight gain (three SNPs on GGA5, one SNP on GGA1 and one SNP on GGA3), plasma coloration measured as optical density at wavelengths in the range 465–510 nm (10 SNPs and all on GGA10) and the percentage of β2-globulin in blood plasma (15 SNPs on GGA1 and one SNP on GGA2). Biological pathways related to metabolic processes, cell proliferation, and primary innate immune processes were among the most frequent significantly enriched biological pathways. Furthermore, the network-based analysis produced two networks of high confidence, with one centered on large tumor suppressor kinase 1 (LATS1) and 2 (LATS2) and the second involving the myosin heavy chain 6 (MYH6).

Conclusions

We identified several strong candidate genes and genomic regions associated with traits measured in response to Eimeria maxima in broilers. Furthermore, the post-GWAS functional analysis indicates that biological pathways and networks involved in tissue proliferation and repair along with the primary innate immune response may play the most important role during the early stage of Eimeria maxima infection in broilers.

Electronic supplementary material

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

Detection of two QTL on chicken chromosome 14 for keyhole lymphet haemocyanin

A keyhole lymphet haemocyanin is an antigen which triggers Th1 type of immune response. A QTL for a primary immune response towards keyhole lymphet haemocyanin has been detected on chicken chromosome 14 in three populations. The results from the most recent population were inconsistent and varied depending on the applied QTL detection model. The major goal of the current study was the reanalysis of this data using a 2 QTL model. Additionally, in order to provide more accurate estimates of QTL effects and positions, epistasis between the QTL was considered as a potential important contributor to quantitative traits. Four statistical models were assumed: M1: A model assuming marginal additive effects of two QTL; M2: A model assuming marginal and epistatic additive effects of two QTL; M3: A model assuming marginal additive and dominance effects of two QTL; M4: A model assuming marginal additive and dominance effects of two QTL and all possible pairwise epistases. Two QTL with significant additive and dominance effects were detected on chicken chromosome 14 using model M3. One QTL was detected at 63 cM between MCW0123 and ROS0005, another at 76 cM between ROS0005 and MCW0225/NTN2Lsts1 (FDR = 0.0051). Modelling only additive effects resulted in a significantly worse fit. On the other hand, including epistatic effects did not improve fit significantly. The current study confirms previous reports of the QTL location on GGA14. A notable finding of this study is recognition of two closely related QTL for a keyhole lymphet haemocyanin response at the distal part of chicken chromosome 14.

Identification of parental line specific effects of MLF2 on resistance to coccidiosis in chickens

BACKGROUND

MLF2 was the candidate gene associated with coccidiosis resistance in chickens. Although single marker analysis supported the association between MLF2 and coccidiosis resistance, causative mutation relevant to coccidiosis was not identified yet. Thus, this study suggested segregation analysis of MLF2 haplotype and the association test of the other candidate genes using improved data transformation.

RESULTS

A haplotype probably originated from one parental line was found out of 4 major haplotypes of MLF2. Frequency of this haplotype was 0.2 in parental chickens and its offspring in 12 families. Allele substitution effect of the MLF2 haplotype originated from a specific line was associated with increased body weight and fecal egg count explaining coccidiosis resistance. Nevertheless Box-Cox transformation was able to improve normality; association test did not produce obvious different results compared with analysis with log transformed phenotype.

CONCLUSION

Allele substitution effect analysis and classification of MLF2 haplotype identified the segregation of haplotype associated with coccidiosis resistance. The haplotype originated from a specific parental line was associated with improving disease resistance. Estimating effect of MLF2 haplotype on coccidiosis resistance will provide useful information for selecting animals or lines for future study.

Association of resistance to avian coccidiosis with single nucleotide polymorphisms in the zyxin gene

Our previous genetic studies demonstrated that resistance to avian coccidiosis is linked with microsatellite markers LEI0071 and LEI0101 on chromosome 1. In this study, the associations between parameters of resistance to coccidiosis and single nucleotide polymorphisms (SNP) in 3 candidate genes located between LEI0071 and LEI0101 [zyxin, CD4, and tumor necrosis factor receptor super family 1A (TNFRSF1A)] were determined. The SNP were genotyped in 24 F(1) generation and 290 F(2) generation animals. No SNP were identified in the TNFRSF1A gene, whereas 10 were located in the zyxin gene and 4 in the CD4 gene. At various times following experimental infection of the F(2) generation with Eimeria maxima, BW, fecal oocyst shedding, and plasma levels of carotenoid, nitrite plus nitrate (NO(2)(-) + NO(3)(-)), and interferon-gamma (IFN-gamma) were measured as parameters of resistance. Single marker and haplotype-based tests were applied to determine the associations between the 14 SNP and the parameters of coccidiosis resistance. None of the CD4 SNP were correlated with disease resistance. However, by single marker association, several of the zyxin SNP were significantly associated with carotenoid or NO(2)(-) + NO(3)(-) concentrations. These were the SNP at nucleotide 149 associated with carotenoid at d 3 postinfection (PI), nucleotide 187 with carotenoid at d 6 and 9 PI, and nucleotide 159 with carotenoid between d 3 and 9 PI. In addition, the zyxin SNP at nucleotide 191 was significantly associated with increased levels of NO(2)(-) + NO(3)(-) at d 3 PI. By haplotype association, the zyxin SNP also were found to be highly associated with NO(2)(-) + NO(3)(-) at d 3 PI and increased IFN-gamma at d 6 PI. These results suggest that zyxin is a candidate gene potentially associated with increased resistance to experimental avian coccidiosis.

Microsatellite mapping of QTLs affecting resistance to coccidiosis (Eimeria tenella) in a Fayoumi x White Leghorn cross

Background

Avian coccidiosis is a major parasitic disease of poultry, causing severe economical loss to poultry production by affecting growth and feed efficiency of infected birds. Current control strategies using mainly drugs and more recently vaccination are showing drawbacks and alternative strategies are needed. Using genetic resistance that would limit the negative and very costly effects of the disease would be highly relevant. The purpose of this work was to detect for the first time QTL for disease resistance traits to Eimeria tenella in chicken by performing a genome scan in an F2 cross issued from a resistant Fayoumi line and a susceptible Leghorn line.

Results

The QTL analysis detected 21 chromosome-wide significant QTL for the different traits related to disease resistance (body weight growth, plasma coloration, hematocrit, rectal temperature and lesion) on 6 chromosomes. Out of these, a genome-wide very significant QTL for body weight growth was found on GGA1, five genome-wide significant QTL for body weight growth, plasma coloration and hematocrit and one for plasma coloration were found on GGA1 and GGA6, respectively. Two genome-wide suggestive QTL for plasma coloration and rectal temperature were found on GGA1 and GGA2, respectively. Other chromosme-wide significant QTL were identified on GGA2, GGA3, GGA6, GGA15 and GGA23. Parent-of-origin effects were found for QTL for body weight growth and plasma coloration on GGA1 and GGA3. Several QTL for different resistance phenotypes were identified as co-localized on the same location.

Conclusion

Using an F2 cross from resistant and susceptible chicken lines proved to be a successful strategy to identify QTL for different resistance traits to Eimeria tenella, opening the way for further gene identification and underlying mechanisms and hopefully possibilities for new breeding strategies for resistance to coccidiosis in the chicken. From the QTL regions identified, several candidate genes and relevant pathways linked to innate immune and inflammatory responses were suggested. These results will be combined with functional genomics approaches on the same lines to provide positional candidate genes for resistance loci for coccidiosis. Results suggested also for further analysis, models tackling the complexity of the genetic architecture of these correlated disease resistance traits including potential epistatic effects.

Immunogenomic approaches to study host immunity to enteric pathogens

With increasing consumer demands for safe poultry products, effective control of disease-causing pathogens is becoming a major challenge to the poultry industry. Many chicken pathogens enter the host through the gastrointestinal tract, and over the past few decades, in-feed antibiotics and active vaccination have been the 2 main mechanisms of disease control. However, increasing public concerns are prompting government regulations on the use of growth-promoting drugs in animal production, and the ability of current vaccines to protect against emerging hypervirulent strains of pathogens is becoming an issue. Therefore, there is a need to develop alternative control strategies against poultry pathogens of economic importance as well as to carry out basic research to enhance understanding of host-pathogen interactions at local sites of infection. Effective control strategies against pathogens can only be accomplished by comprehensive analysis of the basic immunobiology of host-pathogen interactions. Recent sequencing of the poultry genome and the availability of several tissue-specific cDNA microarrays are facilitating the rapid application of functional immunogenomic technologies to poultry disease research. Studies using functional genomic, immunology, and bioinformatic approaches have provided novel insights into disease processes and protective immunity to chicken pathogens. In this review, we summarize recent published literature concerning the host response to Eimeria and Salmonella infections with emphasis on our studies using immunogenomic tools to investigate and characterize the mechanisms of avian immunity to these mucosal pathogens. The results clearly indicate that this immunogenomic approach will lead to increased understanding of immune responses to infectious agents that will enable the development of effective prevention strategies against mucosal pathogens.

Mapping of quantitative trait loci affecting resistance/susceptibility to Sarcocystis miescheriana in swine

The outcome of infectious diseases in vertebrates is under genetic control at least to some extent. In swine, e.g., marked differences in resistance/susceptibility to Sarcocystis miescheriana have been shown between Chinese Meishan and European Pietrain pigs, and these differences are associated with high heritabilities. A first step toward the identification of genes and polymorphisms causal for these differences may be the mapping of quantitative trait loci (QTLs). Considering clinical, immunological, and parasitological traits in the above model system, this survey represents the first QTL study on parasite resistance in pigs. QTL mapping was performed in 139 F(2) pigs of a Meishan/Pietrain family infected with S. miescheriana. Fourteen genome-wide significant QTLs were mapped to several chromosomal areas. Among others, major QTLs were identified for bradyzoite numbers in skeletal muscles (F = 17.4; p < 0.001) and for S. miescheriana-specific plasma IgG(2) levels determined 42 days p.i. (F = 20.9; p < 0.001). The QTLs were mapped to different regions of chromosome 7, i.e., to the region of the major histocompatibility complex (bradyzoites) and to an immunoglobulin heavy chain cluster, respectively. These results provide evidence for a direct and causal role for gene variants within these gene clusters (cis-acting) in differences in resistance to S. miescheriana.