Research | Open | Published:
Single nucleotide polymorphisms in candidate genes are significantly associated with resistance to Haemonchus contortus infection in goats
Journal of Animal Science and Biotechnologyvolume 10, Article number: 30 (2019)
Haemonchosis is a major economic problem in goat production in humid, tropical and subtropical regions. The disease is caused by an abomasal nematode, Haemonchus contortus, which is highly pathogenic in small ruminants. The aim of this study was to identifying single-nucleotide polymorphisms (SNP) that were associated with fecal egg counts (FEC) and could be used as markers to identify resistance to H. contortus in goats.
Ten novel variants in the CIITA, ATP2A3, HSPA8, STAT5B, ESYT1, and SERPING1 genes were associated with FEC in goats with a nominal significance level of P < 0.05. Two missense mutation in the exon region of the caprine CIITA gene resulted in replacement of arginine with cysteine at position 9473550 (R9473550C) and aspartic acid with glutamic acid at position 9473870 (D9473870E). Chinese goat breeds had significantly higher FEC than Bangladeshi goat breeds within their respective genotypes. Polymorphism information content (PIC), effective allele number (Ne), and heterozygosity (He) were greatest for the STAT5B_197_A > G SNP locus in all goat breeds. Pairwise coefficients of linkage disequilibrium (D´, r2) revealed complete LD (r2 = 1) between significant SNP polymorphisms in CIITA and SERPING1 and strong LD (r2 = 0.93 and 0.98) between polymorphisms in HSPA8 and ATP2A3, respectively. Correlation coefficient (r) between FEC and body weight (BW) was significantly positive (r = 0.56***, P < 0.001) but that between FEC and packed cell volume (PCV) was negatively significant (r = − 0.47**, P < 0.01) in the total population of goats. On the other hand, correlation coefficient (r) between BW and PCV was not significant in total population of goats. Association analysis revealed that haplotypes within ATP2A3, HSPA8, and SERPING1 were significantly associated with FEC. Quantitative real-time PCR revealed that the relative expression of mRNA was higher (P < 0.001) for resistant, compared to susceptible, groups of goats for all candidate genes except CIITA.
This study identified SNP markers that can potentially be used in marker-assisted selection programs to develop goat breeds that are resistant to H. contortus.
Livestock, and especially small ruminants, provide one of the principal means to improve living standards in many developing countries and play a critical role in both the national economy and livelihoods in rural communities [1, 2]. Among small ruminant species, goats are particularly susceptible to gastrointestinal nematode (GIN) parasites in humid, tropical and subtropical regions where goats are reared under natural grazing conditions and are especially affected by Haemonchus contortus (H. contortus), haematophagous abomasal nematode parasite [3, 4]. Haemonchosis causes significant economic losses due to negative effects on production and increases in the cost of anthelmintic treatment [5,6,7].
Livestock producers and breeders are therefore trying to utilize genetic variation in susceptibility to parasitic infection within and among breeds to develop lines that are resistant to H. contortus [8,9,10]. The advance of molecular genetics has provided an opportunity to identify genes and associated pathways which provide a better understanding of the physiological processes associated with targeted genes and can potentially be used as genetic markers in marker-assisted selection programs . Several candidate genes associated with parasite susceptibility have been identified and partially characterized in different farm animal species. However, only a few studies have addressed candidate genes and polymorphisms associated with resistance to GIN in goat [6, 12, 13].
Gene expression studies revealed that resistant animals expressed hundreds of genes involved in immune function . Genetic variants that affect phenotypic variation are generally identified by association analysis . Such association studies have already been done in different animals including human [16, 17]. Major histocompatibility complex class II molecules (CIITA), the endoplasmic Ca2+ pump (ATP2A3), heat stress protein 70 (HSPA8), serpin peptidase inhibitor clade G member 1 (SERPING1), signal transducer activator of transcription 5B (STAT5B), and extended synaptotagmin-1 (ESYT1) have been shown to impact immune function in humans and several livestock species. Polymorphisms in these genes have been shown to be associated with autoimmune diseases in human , mice , and other farm animal species such as sheep, goat, yak, buffalo, and chicken [16, 19, 20], but comparable association studies involving resistance to GIN and, particularly, to Haemonchus contortus in goats have not been reported. The present research was therefore undertaken to search for associations between single-nucleotide polymorphism (SNP) in CIITA, ATP2A3, HSPA8, STAT5B, ESYT1 and SERPING1 and resistance to H. contortus using fecal egg counts (FEC) as an indicator of the level of parasite infection.
All experiment were carried out in accordance with the Guide for the Care and Use of Laboratory Animals (Ministry of Science and Technology of China, 2006) and approved by the Standing Committee of the Hubei People’s Congress, and the Ethics Committee of Huazhong Agricultural University, Wuhan, China (Permission Number 4200896859).
A total of 507 goats were used for the study and included 32 Yichang White (YCW), 56 Nanjiang Yellow (NJY), 37 Enshi Black (ESB), and 155 crossbred goats (CCB; produced by the mating YCW, NJY, and ESB goats) from China and 227 Black Bengal goats (BBG) from Bangladesh. Experimental animals were all grazed on native pastures and had not been dewormed for at least 1 year to maximize the likelihood of natural parasite infection. Goats were selected from three locations: Enshi city (30o N, 109o E) and Yichang city (31o N, 111o E) in Hubei province in southern China and the Natore district (24o N, 89o E) in western Bangladesh. The animal pedigree was ensured from the record of animals maintained by the farmers as well as experimental record on the animals, especially in Bangladesh.
Phenotypic data and sample collection
The data were collected at the end of rainy season, from June to August 2015, when parasitic infection was maximum due to high moisture levels and ambient temperatures . Body weights were determined and fecal and peripheral blood samples were collected from each goat. The FEC was determined using a modified McMaster’s technique [21, 22]. Blood was collected using venoject tubes coated with ethylene diaminetetraacetic acid (EDTA) as an anti-coagulant and transferred to a laboratory for DNA extraction using the standard phenol–chloroform method  and determination of packed cell volume (PCV) using the capillary micro-hematocrit method .
Targeted re-sequencing and identification of polymorphisms
A total of 84 candidate genes were identified for targeted re-sequencing and overall 129 polymorphisms (the number of SNP ranged from 1 to 6 for each of the genes) were identified in goats that were involved in immune system and body defense mechanisms. The candidate genes (NCBI Gen Bank association numbers presented in Additional file 1: Table S1) were selected based on analysis of the global list of sheep Entrez Gene IDs in the bovine KEGG database using KEGGARRAY. For the amplification of polymerase chain reaction (PCR) product, oligonucleotide primers were designed for partial regions of genes in a panel of eight unrelated goats selected from different breeds located in major environmental areas under present study. The final PCR reaction volume was 50 μL consisting of a 2-μL template containing 50 ng of genomic DNA, 1.5 μL of each primer, 20 μL of double-distilled H2O, and 25 μL of Taq premix (TaKaRa, Dalian. China). For all primers, samples were incubated at 94 °C for 10 min and then amplified for 35 cycles; each cycle consisted of denaturation at 95 °C for 30 s, annealing at 58 °C for 30 s, an extension at 72 °C for the 50 s, and the final extension at 72 °C for 5 min. The PCR products were checked by electrophoresis on 1% agarose gel to confirm amplification and sequenced using an ABI 3730XL genetic analyzer (Tsingke Biological Technology Co. Ltd., Wuhan, China). Sequences were generated from both ends, and single nucleotide polymorphisms were identified using Codon Code Aligner version 8.0.1 and BioEdit software  and confirmed by manual inspection.
After identification of the candidate SNP, the animals were genotyped using a Kompetitive allele-specific PCR (KASPar) assay based on FRET chemistry (KBiosciences, LGC Genomics. UK). Two forward primers specific to each SNP allele were designed and combined with the respective proprietary tail sequence complementing the FAM or HEX fluorescence reporting system . A common reveres primer was designed for each genotyping assay. Thermal cycling conditions are presented in Additional file 1: Table S2. Genotype calling was based on an endpoint allele discrimination module to measure fluorescence intensities recorded for each of two alleles. The emission data for all the samples on the plate were plotted on corresponding X and Y axis for each allele, and the genotypes were called based on distinct clustering . Touchdown PCR (TD-PCR) assays were used to amplify the variable regions of the different DNA sequences and allow genotyping of each of the target SNP locus. The TD-PCR protocol used a total of 36 cycles of amplification, each including two steps of extension and annealing processes. Touchdown PCR protocol (61–55 TD) was used for genotyping of CIITA and ATP2A3 and protocol (68–62 TD) was used for genotyping of HSPA8, STAT5B, ESYT1 and SERPING1 (Additional file 1: Table S2).
Artificial challenge trial
An artificial challenge with H. contortus was performed at the Lao Gao Huang goat farm, Yichang, China to asses mRNA expression for each candidate gene in goats characterized as resistant or susceptible to H. contortus. An initial sample of 81 Yichang White goats of both sexes with ages of 9 to 12 months was selected from a herd of 450 goats and FEC, circulating haemoglobin levels (Hg), and PCV were determined for those animals. Goats were then divided into two groups based on their FEC. Goats with low FEC (FEC < 500) were assigned to the resistant group and goats with high FEC (FEC > 500) were assigned to the susceptible group. From these groups, 4 resistant goats and 4 susceptible goats with an average age of 305 ± 10 days were selected for artificial challenge. Experimental goats were moved to drylot, treated with Ivermectin (0.25 mg/kg body weight), and, 15 days later, with levamisole (8 mg/kg body weight) [27, 28]. In dry lot, goats were fed GIN-free fresh grass, and a concentrate mixture was provided twice per day at a rate of 1.5% of body weight; pens was cleaned regularly. Complete deworming was confirmed by determination of FEC. When all experimental goats achieved a FEC of zero, resistant and susceptible goats were inoculated with a single dose of 5000 third-stage (L3) larvae of H. contortus. The day of inoculation was considered to be day 0, and FEC, PCV, Hg, and body weight (BW) were determined weekly for 6 weeks after inoculation. At day 42, all experimental goats were sacrificed, contents of the abomasum and small intestine were collected, and mucosa was thoroughly washed into the same containers. Worms from the abomasum and small intestine were then counted.
RNA extraction, cDNA synthesis and qRT-PCR
Approximately 100 g of abomasal tissue including lymph node from the hemorrhagic part of the abomasal wall were collected and total RNA was extracted using TRlzol (Invitrogen, Carlsbad, CA, USA) following the manufacturer’s recommended procedure. A Nano Drop ND2000 spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA) was used to assess the quality of extracted RNA that was utilized to mRNA expression for ATP2A3, SERPING1, CIITA, HSPA8, ESYT1, and STAT5B by quantitative real-time PCR (qRT-PCR) using cDNA generated from approximately 1 μg of RNA using a standard PrimeScipt™ RT reagent kit with gDNA eraser (Perfect Real Time, TAKARA Bio, Inc.). The qRT-PCR was performed using a CFX-96 Bio-Rad thermal cycler with SYBR green real-time PCR master mix (Toyobo Co., Ltd., Osaka, Japan). The qRT-PCR protocol was a single cycle of denaturation at 95°C for 5 min, 45 cycles of denaturation at 95 °C for 20 s, annealing at 58 °C for 20 s, and extension at 72 °C for 15 s. β-actin (ACTB) was used as a housekeeping gene to normalize the samples. Primers used for qRT-PCR are shown in Additional file 1: Table S3.
The FEC were not normally distributed and exhibited positive skewness. A logarithmic transformation [log10 (FEC + 25)] was therefore applied before analysis [9, 29]. Descriptive statistics (arithmetic means, variances, standard errors, etc.; ) were derived using SAS (Version 9; SAS Institute, Inc., Cary, NC USA). No transformation was applied to PCV or body weight.
Genotypic and allelic frequencies were determined for each candidate genes in each population . Correspondence of observed genotypic frequencies to expectations based on Hardy-Weinberg equilibrium (HWE) was assessed by chi-square (χ2) tests in POPGENE software (Version 3.2;) . Population genetic diversity indices [gene heterozygosity (He), polymorphism information content (PIC), and effective allele numbers (Ne)] were calculated using PICcalculator an online program available at https://www.liverpool.ac.uk/~kempsj/pic.html [33,34,35]. Heploview blocks and linkage disequilibrium (LD) measures (D′, r2) were derived for the significant SNP loci using genotypes for all goats in this study and the default settings in the Heploview software. Descriptive statistics were estimated using the partition–ligation, combination–subdivision expectation maximization algorithm  in Heploview software [11, 37].
Relationships between SNP and FEC were tested using R software  and a linear model:
where y is the FEC, μ is the overall population mean, B and V are fixed effects of breed and location associated with the FEC trait, G is one of the 129 SNP which were tested for an association with FEC, and e is residual error . Multiple comparisons between genotypes of significant SNP and fecal egg count (FEC) in all goat population were carried out using R software (Version 3.0.2).
Evaluation of gene expression in the artificial challenge trial
Differences in gene expression between resistant and susceptible goats in the challenge trial were tested using the 2-ΔΔCT method in SAS [40, 41]. Student’s t-tests and a significance level of P < 0.05 were used to compare gene expression in resistant and susceptible goats using Graph Pad Software Prism7 (San Diego, CA USA).
Based on FEC, 237 out of 507 goats were infected with H. contortus, resulting in an overall infection rate as 46.75%. Average infection rates were 68.9% for Chinese goat breeds and 19.4% for BBG goats. Least-squares means for FEC and PCV differed (P < 0.001) between goats from China and Bangladesh. Least-squares means for FEC and PCV also differed (P < 0.001 and P < 0.05, respectively) among breeds, but not between locations, in China (Table 1). Across the five goat breeds evaluated in this study, transformed FEC were highest (2.53 ± 0.10 epg) for ESB goats and lowest (1.58 ± 0.04 epg) for BBG goats. However, the maximum FEC (15,600 epg) was observed in a Chinese crossbred goat (Table 1). Male goats had higher FEC than female goats (P < 0.01), but did not differ from female goats in PCV; this result was consistent across countries In contrast to results for FEC, body weights were similar among breeds, locations, and sexes (Table 1).
Identification of polymorphisms
Out of 129 polymorphisms from 84 different genes, 10 polymorphisms in 6 genes (CIITA, ATP2A3, HSPA8, STAT5B, ESYT1 and SERPING1) had significant associations with FEC at a nominal significance level of P < 0.05. The rest of polymorphisms (119 polymorphisms in 78 genes) had non-significant associations with the studied phenotype (FEC) are presented in Additional file 1: Table S5. Two non-synonymous mutations in CIITA resulted in replacement of arginine with cysteine at position 9473550 and aspartic acid with glutamic acid at position 9473870 (Table 2). The remaining eight SNP were located in non-coding regions and included a SNP in the 3´UTR of STAT5B (A/G at position 41961106), two SNP in the 3´UTR of ATP2A3 (A/G at position 24358932 and A/C at position 24358402), and five SNP in various intron regions (Table 2).
Genotypic and allelic frequencies
Across the entire sample of 507 goats (Table 3), frequencies of CC homozygotes at CIITA_161_C > T (0.98) and TT homozygotes at CIITA_481_A > T (0.98) were very high. Frequencies of GG, GG, and AA homozygotes were also relatively high at ESYT1_559_G > C, HSPA8_1024_A > G, and HSPA8_1064_A > G (0.75, 0.69, and 0.65, respectively). Genotypic frequencies for SNP polymorphisms in STAT5B, SERPING1, and ATP2A3 were more balanced than those of the other genes considered in the study. Frequencies of the C allele at CIITA_161_C > T and T allele at CIITA_481_A > T were very high (0.99). Frequencies of the G allele at ESYT1_559_G > C and HSPA8_1024_A > G and the A allele at HSPA8_1064_A > G likewise exceeded 0.80, but frequencies of the most common allele at the remaining loci were moderate, ranging from 0.52 to 0.76 (Table 3). When Chinese and Bangladeshi goat breeds were considered as separate populations, similar trends were observed for allelic and genotypic frequencies in all candidate genes (Table 4).
Chi-square tests revealed that departures of observed genotypic frequencies from HWE expectations were significant for CIITA_161_C > T, CIITA_481_A > T, HSPA8_1024_A > G, HSPA8_1064_A > G, ESYT1_559_G > C, and STAT5B_197_A > G in Chinese goats and CIITA_161_C > T, CIITA_481_A > T, and ATP2A3_150_A > G in Bangladeshi goats (Table 4). Significant departures from HWE conditions for CIITA_161_C > T and CIITA_481_A > T in both populations are likely spurious results, reflecting very high allelic frequencies for the most common allele and limited sample sizes. Genotypic frequencies at other loci did not differ from HWE expectations
Population genetic diversity
Indicators of population genetic diversity such as the He, Ne, and PIC were highest for STAT5B_197_A > G and highest for YCW goats (0.50, 2.00, and 0.38, respectively) followed by CCB, BBG, ESB, and NJY goats. The lowest values for He, Ne, and PIC were observed for CIITA in all goat breeds. Values for He, Ne, and PIC were relatively high for SERPING1 and ATP2A3 in all goats breeds but were consistently moderate for the remaining loci (Table 5). Based on PIC, STAT5B, SERPING1, and ATP2A3 were categorized as having medium genetic diversity (0.25 < PIC < 0.50) for all goat breeds. CIITA, HSPA8, and ESYT1 were categorized as having low genetic diversity (PIC < 0.25) for all breeds except CCB and NJY for HSPA8_1024_A > G and HSPA8_1064_A > G and YCW for ESYT1_559_G > C, which were classified as having medium genetic diversity (Table 5).
Population association analysis
Association analysis confirmed that 10 SNP from six genes (CIITA, HSPA8, STAT5B, ESYT1, ATP2A3, and SERPING1) were significantly associated with FEC, with nominal significance levels of P < 0.05 to P < 0.01 (Table 6). Highly significant (P < 0.01) associations were observed for HSPA8_1024_A > G, HSPA8_1064_A > G, SERPING1_312_C > T, SERPING1 615_G > T, and ATP2A3_680_A > C (Table 6). Large differences were observed between TT and CC genotypes (P = 0.003) and TT and TC genotypes (P = 0.004) at SERPING1_312_C > T, between TT and GG genotypes (P = 0.002) and TT and GT genotypes (P = 0.005) at SERPING1_615_G > T, and between CC and AA genotypes (P = 0.012) at ATP2A3_680_A > C. Differences between some others genotypes also reached lower levels of significance (P < 0.05) (Table 6).
Means and standard errors for FEC, PCV, and body weight (BW) for each genotype for candidate genes with significant effects on FEC are presented in Table 7. Across all goat breeds, generally undesirable effects (i.e., higher FEC) were associated with the less-frequent allele at all loci except STAT5B_197_A > G. Dominance relationships for these unfavorable alleles varied among loci. Additive gene action was indicated only for ATP2A3_680_A > C and, perhaps, ATP2A3_150_A > G. Allelic frequencies were similar for the two alleles at STAT5B_197_A > G, with an apparent recessive favorable effect of the G allele. Particularly large and apparently dominant unfavorable effects of rare minor alleles were observed for both CIITA loci. Significant differences among SNP genotypes in PCV were not observed for any of the loci. Effects of genotype on body weight were detected only for the two ATP2A3 loci, with approximately additive effects for ATP2A3_150_A > G, but apparent dominance of the A allele for ATP2A3_680_A > C.
The FEC for Chinese goat breeds were higher than those of Bangladeshi goats (Table 1, Additional file 1: Table S4). There were significant differences for FEC among all the genotypes within significant SNP loci for both Chinese and Bangladeshi goat breeds except for genotypes within HSPA8 and CIITA in Bangladeshi goats.
Correlation among measured traits
Correlation coefficients (Table 8) between PCV and FEC were negative and significant in Chinese goats (r = − 0.45*; P < 0.05), Bangladeshi goats (r = − 0.44***; P < 0.001) and also in total population of goats (r = − 0.47; P < 0.001). On the other hand, Correlation coefficient between PCV and BW was non-significant in Chinese goats, Bangladeshi goats and total population of goats. The correlations between BW and FEC was significantly positive (r = 0.66*; P < 0.05) in Chinese goats, but non-significant and negative (− 0.27) in Bangladeshi goats. In total population of goats, the correlation between BW and FEC was highly significant and positive (r = 0.56; P < 0.001) (Table 8).
Linkage disequilibrium (LD) and haplotype association analysis
Coefficients of linkage disequilibrium (D´, r2) were determined for each pair of SNP. Linkage disequilibrium was present between HSPA8_1024_A > G and HSPA8_1064_A > G, CIITA_161_C > T and CIITA_481_A > T, SERPING1_312_C > T and SERPING1_615_G > T, and ATP2A3_150_A > G and ATP2A3_680_A > C (Fig. 1). In particular, SNP in the intronic region of SERPING1 and the exon region of CIITA were in complete LD (D´ = 100 and r2 = 1). The SNP variants in the intronic region of HSPA8 and the 3´ UTR region of ATP2A3 also expressed strong LD (D´ ≥ 93 and r2 ≥ 0.93; Fig. 1). High levels of LD between SNP within these genes suggest that significant effects of these SNP in Table 7 were likely associated with common causal regions.
Associations between 13 haplotypes at CIITA, HSPA8, SERPING1, and ATP2A3 and FEC were examined (Table 9). Haplotypes in ATP2A3 were associated (P < 0.05) with FEC in YCW, NJY, and CCB goats. Associations with FEC (P < 0.05) were also observed for haplotypes in HSPA8 in NJY goats and SERPING1 in CCB goats. The CCB goats that were heterozygous for the CG and TC chromosome segments had lower FEC compared to those that were homozygous for the TT chromosome segment, but this relationship was not evident for the other breeds of Chinese goats. No significant difference (P > 0.05) were found between any of the haplotypes and FEC in BBG goats.
Artificial challenge trial
The FEC increased gradually from day 21 through days 42 and was higher in susceptible goats compared to resistant goats during that period (Fig. 2a). Additionally, the body weight of resistant goats increased throughout the experimental periods whereas that of susceptible goats was stable during the experimental period (Fig. 2b). PCV and Hg level began to decline at days 28 and 14, respectively, in susceptible goats and continued to decline through day 42 (Fig. 2c and d). On the other hand, PCV and Hg values in resistant goats did not decline as rapidly as those in susceptible goats (Fig. 2c and d).
Differential gene expression analysis
Relative levels of expression of mRNA in resistant goats were higher than those in susceptible goats for HSPA8, ESYT1 and SERPING1 (all P < 0.001) and for ATP2A3 and STAT5B (both P < 0.0001) (Fig. 3). Resistant and susceptible goats did not differ in expression of CIITA.
The FEC is an important phenotypic marker for resistance to H. contortus in sheep and goats [42,43,44]. Host genetics significantly affect FEC and can be used to assess the level of H. contortus infection . The current study revealed that Black Bengal goats of Bangladesh had lower FEC than Chinese goat breeds, supporting previous work . Variation among the breeds for FEC may have reflected both genetic and non-genetic differences, because the Bangladeshi and Chinese breeds were evaluated separately in their home countries. However, goats in the two locations were maintained under similar environmental condition, feces were collected at the same period, and FEC were determined using the same method, supporting the hypothesis that the breeds differed genetically in resistance. This is the first time that goats from the two countries have been evaluated for resistance to H. contortus using comparable experimental protocols.
Ten novel SNP in CIITA, ATP2A3, HSPA8, STAT5B, ESYT1, and SERPING1 genes were present in both Chinese and Bangladeshi goats and may provide a basis for identification of genetic markers associated with H. contortus resistance . CIITA, ATP2A3, HSPA8, STAT5B, ESYT1, and SERPING1 have all been shown to affect regulation of the immune system in humans [46,47,48], laboratory mice , and a few livestock species such as yak , chicken , and sheep . CIITA produces MHC class II proteins that are found on surfaces of several vertebrate immune cells and regulate immune responses to various pathogens. A role for this gene in immune function has been clearly demonstrated in mice and humans but only rarely been reported in lower vertebrates such as fishes . ATP2A3 encodes SERCA Ca2+-ATPase enzymes which play a fundamental role in maintaining intracellular homeostasis by supporting pumping of Ca2+ into endoplasmic reticulum of muscle cells. Expression of ATP2A3 has also been reported to be dramatically reduced in human colon, breast, and lung cancers [48, 51]. Mutations in SERPING1 cause hereditary angioedema in humans . HSPA8 has been reported have to anti-tumor effects by inducing chemokine production from tumor cells and activation of chemo-attracted dendritic cells via the TLR4 pathway in mice . Extended synaptotagmin-1 (ESYT1) acts as a Ca2+-regulated lipid-transfer protein and membrane-fusion regulator, which facilitates intracellular signaling and expression of immune responses to viruses . Polymorphisms in STAT5B have been associated with production and growth traits in chicken and cattle [19, 53]. The only reported association between polymorphisms in these genes and GIN resistance in goats was a positive association involving a polymorphism in ATP2A3 . However, there are several reports [6, 12, 13] of associations between polymorphism in DRB-1, IGF-1, IL2, IL13, and IFNG and GIN resistance in goats.
Population genetic parameters can be used to characterize genetic diversity within and among populations. The PIC is an indicator of the extent of polymorphism in a population and the value of PIC has been classified as high, intermediate, and low levels of polymorphism for PIC > 0.50, 0.25 < PIC < 0.50, and PIC < 0.25, respectively [34, 53, 55]. Results of the current study revealed low to intermediate levels of polymorphism for all genes except CIITA in Black Bengal goats from Bangladesh. Similar or somewhat higher levels of genetic diversity were observed for Chinese goats compared to Bangladeshi, goats. These results indicated that there was adequate genetic diversity for selection to be effective in controlling GIN infections for all genes except, perhaps, CIITA. The SNP polymorphisms in CIITA resulted in amino acid substitutions in the gene product and therefore had potential to have a functional association with GIN resistance. Given the low frequency of the deleterious allele at CIITA, these polymorphisms may reflect the presence of rare, deleterious mutations. Three significant mutations were detected in the 3′ UTR regions of STAT5B and ATP2A3 and may alter gene regulation and impact gene expression through mechanisms that disrupt miRNA binding . Previous studies have shown that mutations in introns of bovine, caprine, swine and human genes were significantly associated with performance [16, 57,58,59]. Synonymous mutations can therefore still be used as genetic markers if they are in LD with polymorphisms in functional genes.
Linkage disequilibrium (LD) plays a vital role in mapping genes that affect complex diseases and identifying association among genetic markers and functional genes . Understanding LD among SNP also avoids redundant inferences involving non-independent genetic markers. Result of this study indicated that four pairs of variants were in significant LD with each other. The two polymorphisms in CIITA and the two polymorphisms in SERPING1 were each in complete LD (r2 = 1), and, within each gene, represented by only two unique haplotypes. High levels of LD were also observed for polymorphisms in ATP2A3 (r2 = 0.93) and HSPA8 (r2 = 0.98). Polymorphisms with r2 > 0.33 are generally considered to be in relatively strong LD  and commonly inherited together. When r2 values are large, haplotype analysis is preferred to the analysis of individual SNP variants . The LD analysis revealed that haplotype differences in ATP2A3, HSPA8, and SERPING1 were significantly associated with FEC.
The artificial challenge trial revealed significantly higher expression levels for all candidate genes except CIITA in GIN-resistant, compared to GIN-susceptible, goats. Resistant and susceptible goats were identified by screening a large sample of goats exposed to natural GIN infection. Animals were then dewormed and re-infected with H. contortus. Results therefore addressed the consistency of indicators of parasite resistance and subsequent levels of gene expression in different infection cycles. Goats that had lower FEC, indicating greater resistance to infection, under natural grazing conditions appeared much more resistant, based on indicators of parasite resistance, and also had greater expression of five of the six candidate genes in response to a subsequent controlled GIN infection. The high repeatability for indicators of parasite resistance across infection cycles was consistent with results reported for Pelibuey hair sheep , and differences in expression of genes involved in immune function have been reported between resistant and susceptible sheep breeds . The results of present study show that the majority of candidate genes selected for this study were differentially expressed in resistant and susceptible goats, further supporting the potential value of SNP in these genes as markers for GIN resistance in goats.
Population genetic parameter, LD among SNP markers, identification of non-synonymous mutations in candidate genes, differences in relative gene expression between resistant and susceptible goats, and associations with FEC involving both individual SNP genotypes and SNP haplotypes can be used to advance our understanding of options to utilize selective breeding and molecular markers to improve resistance to H. contortus and other GIN in goats. Ten SNP within six candidate genes were associated with FEC and provided a suite of potential molecular markers for further study and possible use in screening individuals for resistance to H. contortus.
Black Bengal goats
Fecal egg count
Kyoto encyclopedia of genes and genomes
Larval stage 3
Logarithmic transformation fecal egg count
Packed cell volume
Polymorphism information content
Single nucleotide polymorphisms
Bacha A. Study on prevalence of gastrointestinal nematodes and coccidian parasites affecting cattle in west arsi zone, Ormia regional state, Ethiopia. Glob J Anim Sci Res. 2014;3:77–86.
Ahmed J, Duguma A, Regassa D, Belina D, Jilo R. Gastrointestinal nematode parasites of small ruminants and anthelmintics efficacy test in sheep of Haramaya District, Eastern Ethiopia. Anim Vet Sci. 2017;5:39–44.
Traoré A, Notter DR, Soudre A, Kaboré A, Álvarez I, Fernández I. Resistance to gastrointestinal parasite infection in Djallonké sheep. Animal. 2017;11:1354–62.
Omar A, Alam M, Faruque M, Mondal M, Raihan M, Adetula A. Variability in fecal egg count of Haemonchus contortus infection to native goat breeds of China and Bangladesh under natural grazing condition. Progress Agric. 2017;27:473–81.
Kim E-S, Sonstegard TS, da Silva MV, Gasbarre LC, Van Tassell CP. Genome-wide scan of gastrointestinal nematode resistance in closed Angus population selected for minimized influence of MHC. PLoS One. 2015;10:0119380.
Bressani F, Tizioto P, Giglioti R, Meirelles S, Coutinho R, Benvenuti C. Single nucleotide polymorphisms in candidate genes associated with gastrointestinal nematode infection in goats. Gene Molec Res. 2014;13:8530–6.
Howell SB, Burke JM, Miller JE, Terrill TH, Valencia E, Williams MJ. Prevalence of anthelmintic resistance on sheep and goat farms in the southeastern United States. J Amer Vet Med Assoc. 2008;233:1913–9.
Mandonnet N, Aumont G, Fleury J, Arquet R, Varo H, Gruner L. Assessment of genetic variability of resistance to gastrointestinal nematode parasites in creole goats in the humid tropics. J Anim Sci. 2001;79:1706–12.
Rout P, Chauhan K, Matika O, Bishop S. Exploring the genetic resistance to natural gastrointestinal nematode infection in Indian goats. Vet Parasitol. 2011;180:315–22.
Chauhan K, Rout P, Singh P, Mandal A, Singh S, Roy R. Genetic resistance of Barbari and Jamunapari kids to natural infection with gastrointestinal nematodes. Trop Anim Health Prod. 2003;35:397–408.
Su L, Mei S, Tao H, Peng X, Sun X, Wu H. Identification of the promoter region and genetic mutations of the porcine GALP gene. Mol Biol Rep. 2013;40:2821–7.
Alim M, Fu Y, Wu Z. Zhao S, Cao J. Single nucleotide polymorphisms of toll-like receptors and association with Haemonchus contortus infection in goats. Pak Vet J. 2016;36:286–91.
Asif AR, Qadri S, Yuhua F, Alim M, Wu Z, Ijaz N. Single nucleotide polymorphisms in DRB1, IGF1 and ILs associated with fecal egg count confers resistance against Haemonchus contortus infection in goats. Pak J Agric Sci. 2016;53.
Diez-Tascón C, Keane OM, Wilson T, Zadissa A, Hyndman DL, Baird DB. Microarray analysis of selection lines from outbred populations to identify genes involved with nematode parasite resistance in sheep. Physiol Genomics. 2005;21:59–69.
Lee J-B, Yoo C-K, Park H-B, Cho I-C, Lim H-T. Association of the single nucleotide polymorphisms in RUNX1, DYRK1A, and KCNJ15 with blood related traits in pigs. Asia Austral J Anim Sci. 2016;29:1675.
Lan D, Lin B, Xiong X, Yang X, Li J. Identification and characteristics analysis of toll-like receptors family genes in yak. Gene Genom. 2016;38:429–38.
Desnues B, Macedo AB, Roussel-Queval A, Bonnardel J, Henri S, Demaria O. TLR8 on dendritic cells and TLR9 on B cells restrain TLR7-mediated spontaneous autoimmunity in C57BL/6 mice. Proc Nati Acad Sci USA. 2014;111:1497–502.
Dai Y, Condorelli G, Mehta JL. Scavenger receptors and non-coding RNAs: relevance in atherogenesis. Cardiovasc Res. 2015;109:24–33.
Niknafs S, Javaremi AN, Sadeghi M. Single nucleotide polymorphisms in BMPR-IB and STAT5B genes and their association with growth and reproductive traits in chicken. Song J Sci Technol. 2014;36:137–142.
Rout P, Kaushik R, Ramachandran N. Differential expression pattern of heat shock protein 70 gene in tissues and heat stress phenotypes in goats during peak heat stress period. Cell Stress Chaperones. 2016;21:645–51.
Zajac AM, Conboy GA. Veterinary Clinical Parasitology. 2012, Chichester, West Sussex, UK: Wiley-Blackwell, 8.
Gordon HM, Whitlock H. A new technique for counting nematode eggs in sheep faeces. J Coun Sci Indus Res. 1939;12:50–2.
Sambrook J, Russell DW. Molecular cloning: a laboratory manual. 3rd ed. New York: Cold Spring Harbor Laboratory Press; 2001.
Jain N. Normal values in blood of laboratory, fur-bearing, and miscellaneous zoo, domestic, and wild animals. Schalm’s Veterinary Hematology 4th ed. NC Jain Philadelphia, Lea & Febiger; 1986. p. 337.
Hall TA, editor. BioEdit: a user-friendly biological sequence alignment editor and analysis program for windows 95/98/NT. Nucleic Acids Symp Ser. 1999;41:1979–2000.
Periasamy K, Pichler R, Poli M, Cristel S, Cetrá B, Medus D, et al. Candidate gene approach for parasite resistance in sheep–variation in immune pathway genes and association with fecal egg count. PLoS One. 2014;9:88337.
Bambou J-C, Chevrotière C, Varo H, Arquet R, Kooyman FN, Mandonnet N. Serum antibody responses in creole kids experimentally infected with Haemonchus contortus. Vet Parasitol. 2008;158:311–8.
Kurukulasuriya MS, Silva P, Dematawewa C, Ariyaratne H, Rajapakshe R. Responses of Sri Lankan indigenous goats and their Jamnapari crosses to artificial challenge with Haemonchus contortus. Vet Parasitol Reg Stud Rep. 2018;11:41–8.
González-Garduño R, Mendoza-de Gives P, Torres-Hernández G. Variability in the fecal egg count and the parasitic burden of hair sheep after grazing in nematode infected paddocks. Pes Vet Brasil. 2013;33:469–75.
Chiejina S, Musongong G, Fakae B, Behnke J, Ngongeh L, Wakelin D. The modulatory influence of Trypanosoma brucei on challenge infection with Haemonchus contortus in Nigerian west African dwarf goats segregated into weak and strong responders to the nematode. Vet Parasitol. 2005;128:29–40.
Falconer D, Mackay T. Introduction to quantitative genetics. 4th ed. London: Longman Scientific and Technical; 1996.
Yeh F, Yang R, Boyle T. Microsoft windows-based free ware for population genetic analysis (Release 1.31). Molecular Biology and Biotechnology Center, University of Alberta, Edmonton, Canada; 1999.
Nei M, Roychoudhury AK. Sampling variances of heterozygosity and genetic distance. Gene. 1974;76:379–90.
Botstein D, White RL, Skolnick M, Davis RW. Construction of a genetic linkage map in man using restriction fragment length polymorphisms. Amer J Human Genet. 1980;32:314.
Nagy S, Poczai P, Cernák I, Gorji AM, Hegedűs G, Taller J. PICcalc: an online program to calculate polymorphic information content for molecular genetic studies. Biochem Genet. 2012;50:670–2.
Wang A, Zhang Y, Li M, Lan X, Wang J, Chen H. SNP identification in FBXO32 gene and their associations with growth traits in cattle. Gene. 2013;515:181–6.
Yong Y, Lin H. SHEsis, a powerful software platform for analyses of linkage disequilibrium, haplotype construction, and genetic association at polymorphism loci. Cell Res. 2005;15:97–8.
Archana P, Aleena J, Pragna P, Vidya M, Niyas A, Bagath M. Role of heat shock proteins in livestock adaptation to heat stress. J Dairy Vet Anim Res. 2017;5:00127.
Davis BW, Seabury CM, Brashear WA, Li G, Roelke-Parker M, Murphy WJ. Mechanisms underlying mammalian hybrid sterility in two feline interspecies models. Mol Biol Evol. 2015;32:2534–46.
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2−ΔΔCT method. Methods. 2001;25:402–8.
Ling D. SASqPCR: robust and rapid analysis of RT-qPCR data in SAS. PLoS One. 2012;7:29788.
Amarante A, Bricarello P, Rocha R, Gennari S. Resistance of Santa Ines, Suffolk and Ile de France sheep to naturally acquired gastrointestinal nematode infections. Vet Parasitol. 2004;120:91–106.
Bricarello P, Gennari S, Oliveira-Sequeira T, Vaz C, de Gonçalves IG, Echevarria F. Response of Corriedale and Crioula lanada sheep to artificial primary infection with Haemonchus contortus. Vet Res Commun. 2002;26:447–57.
Vanimisetti H, Greiner S, Zajac A, Notter D. Performance of hair sheep composite breeds: resistance of lambs to Haemonchus contortus. J Anim Sci. 2004;82:595–604.
Zhang Z-R, Liu Y-P, Yao Y-G, Jiang X-S, Du H-R, Zhu Q. Identification and association of the single nucleotide polymorphisms in calpain3 (CAPN3) gene with carcass traits in chickens. BMC Genet. 2009;10:10.
Su M, Chen H, Zeng Y, Chen T, Chen J, Jiang L. Association of genetic variation of CIITA and NTCP with chronic hepatitis B virus infection in Han Chinese populations. Hepa Month. 2017;17.
López-Lera A, Pernia O, López-Trascasa M. Expression of the SERPING1 gene is not regulated by promoter hypermethylation in peripheral blood mononuclear cells from patients with hereditary angioedema due to C1-inhibitor deficiency. Orph J Rare Dis. 2014;9:103.
Flores-Peredo L, Rodriguez G, Zarain-Herzberg A. Induction of cell differentiation activates transcription of the sarco/endoplasmic reticulum calcium-ATPase 3 gene (ATP2A3) in gastric and colon cancer cells. Mole Carcino. 2017;56:735–50.
Singh K, Singh S, Ganguly I, Nachiappan RK, Ganguly A, Venkataramanan R. Association of heat stress protein 90 and 70 gene polymorphism with adaptability traits in Indian sheep (Ovis aries). Cell Stress Chaperones. 2017;5:675–684.
Liu Y, Meng Y, Wang Q, Sha Z. Class II, major histocompatibility complex, transactivator (CIITA) in channel catfish: identification and expression patterns responding to different pathogens. Mol Biol Rep. 2012;39:11041–50.
Izquierdo-Torres E, Rodríguez G, Meneses-Morales I, Zarain-Herzberg A. ATP2A3 gene as an important player for resveratrol anticancer activity in breast cancer cells. Molec Carcino. 2017;56:1703.
El Kasmi I, Khadivjam B, Lackman M, Duron J, Bonneil E, Thibault P. Extended synaptotagmin 1 interacts with herpes simplex virus 1 glycoprotein M and negatively modulates virus-induced membrane fusion. J Virol. 2018;92:01281–17.
Capomaccio S, Milanesi M, Bomba L, Cappelli K, Nicolazzi EL, Williams JL. Searching new signals for production traits through gene-based association analysis in three Italian cattle breeds. Anim Genet. 2015;46:361–70.
Bhuiyan AA, Ni P, Elokil AA, Li J, Bhuyan AA, Zhao S, Du X. The polymorphisms in ATPase sarco / endoplasmic reticulum Ca2+ transporting-3geneassociated with gastrointestinal nematode infection in goat. Proc World Congr Gen Appl Livest Prod. 2018;2:145.
Wang J, Li G, Elzo MA, Yan L, Chen S, Jia X. A novel single nucleotide polymorphism of the POU1F1 gene associated with meat quality traits in rabbits. Ann Anim Sci. 2015;15:611–20.
Skeeles LE, Fleming JL, Mahler KL, Toland AE. The impact of 3′ UTR variants on differential expression of candidate cancer susceptibility genes. PLoS One. 2013;8:58609.
Mačeková S, Bernasovský I, Gabriková D, Bôžiková A, Bernasovská J, Boroňová I. Association of the FTO rs9939609 polymorphism with obesity in Roma/gypsy population. Amer J Phys Anthropol. 2012;147:30–4.
Msalya G, Shimogiri T, Nishitani K, Okamoto S, Kawabe K, Minesawa M. Indels within promoter and intron 1 of bovine prion protein gene modulate the gene expression levels in the medulla oblongata of two Japanese cattle breeds. Anim Genet. 2010;41:218–21.
Gao Y, Zhang YH, Zhang S, Li F, Wang S, Dai L. Association of A-FABP gene polymorphism in intron 1 with meat quality traits in Junmu no. 1 white swine. Gene. 2011;487:170–3.
Hazelett DJ, Conti DV, Han Y, Al Olama AA, Easton D, Eeles RA. Reducing GWAS complexity. Cell Cycle. 2016;15:22–4.
Ren G, Huang Y-Z, Wei T-B, Liu J-X, Lan X-Y, Lei C-Z. Linkage disequilibrium and haplotype distribution of the bovine LHX4 gene in relation to growth. Gene. 2014;538:354–60.
Daly MJ, Rioux JD, Schaffner SF, Hudson TJ, Lander ES. High-resolution haplotype structure in the human genome. Nat Genet. 2001;29:229.
Ojeda-Robertos N, Torres-Acosta J, González-Garduño R, Notter D. Phenotypic expression of parasite susceptibility to Haemonchus contortus in pelibuey sheep. Vet Parasitol. 2017;239:57–61.
MacKinnon K, Burton J, Zajac A, Notter D. Microarray analysis reveals difference in gene expression profiles of hair and wool sheep infected with Haemonchus contortus. Vet Immunol Immunopathol. 2009;130:210–20.
We gratefully thank all the people whose hard work in laboratories and in the field in China and Bangladesh contributed to this work.
This work was supported by the National Nature Science Foundation of China (grant no. 31402040) and by a CRP project (grant no. 16087 for China and grant no. 16073 for Bangladesh) from the International Atomic Energy Agency.
Availability of data and materials
All data generated or analyzed during this study are included in this published article.
Ethics approval and consent to participate
All experiment were carried out in accordance with the Guide for the Care and Use of Laboratory Animals and approved by the Standing Committee of the Hubei People’s Congress, and the Ethics Committee of Huazhong Agricultural University, Wuhan, China.
Consent for publication
The authors declare that they have no competing interests.
Table S1. Primer sequences and PCR conditions for genotyping SNP in caprine candidate genes. Table S2. Cycling condition for Touchdown PCR (TD-PCR) genotyping of SNP. Table S3. Primer information for measurement of expression of mRNA by quantitative reverse-transcriptase polymerase chain reaction. Table S4. Least-squares means and standard errors for FEC, PCV and BW for each SNP genotype in Chinese and Bangladeshi goats. Table S5. Results of association and multiple comparison test P-values of non-significant polymorphisms. (DOCX 57 kb)