This study aimed to identify genomic regions, pathways, and putative candidate genes associated with resistance to gastrointestinal nematode in Santa Ines sheep. The phenotypic information comprised 5529 records from 1703 naturally infected animals. After genomic data quality control, 37,511 SNPs from 589 animals were available. The weighted single-step approach for genome-wide association study was performed to estimate the SNP effects and variances accounted by 10-SNP sliding windows. Confirming the polygenic nature of the studied traits, 20, 22, 21, and 19 genomic windows that explained more than 0.5% of the additive genetic variance were identified for fecal egg counts (FEC), Famacha© (FAM), packed cell volume (PCV), and total plasma protein (TPP), respectively. A total of 81, 122, 106, and 101 protein-coding genes were found in windows associated with FEC, FAM, PCV, and TPP, respectively. Several protein-coding genes related to the immune system and inflammatory response functions were identified within those genomic regions, such as ADCY9, ADRB2, BRAF, CADM1, CCL20, CD70, CREBBP, FNBP1, HTR4, IL16, IL22, IL26, MAPK8, NDFIP1, NLRC3, PAK5, PLCB1, PLCB4, ROCK1, TEK, TNFRSF12A, and VAV1. Functional enrichment analysis by DAVID tool also revealed many significant (P < 0.05) pathways and Gene Ontology terms that could be related to resistance to gastrointestinal nematode in Santa Ines sheep, such as chemokine signaling pathway (oas04062), cAMP signaling pathway (oas04024), cGMP-PKG signaling pathway (Oas04022), platelet activation (Oas04611), Rap1 signaling pathway (oas04015), and oxidoreductase activity, acting on paired donors, with incorporation or reduction of molecular oxygen (GO:0016705). These results contribute to improving the knowledge of the genetic architecture of resistance to gastrointestinal nematode in Santa Ines sheep.

BACKGROUND: The aim of this study was to estimate variance components and to identify genomic regions and pathways associated with resistance to gastrointestinal parasites, particularly Haemonchus contortus, in a breed of sheep adapted to tropical climate. Phenotypes evaluations were performed to verify resistance to gastrointestinal parasites, and were divided into two categories: i) farm phenotypes, assessing body condition score (BCS), degree of anemia assessed by the famacha chart (FAM), fur score (FS) and feces consistency (FC); and ii) lab phenotypes, comprising blood analyses for hematocrit (HCT), white blood cell count (WBC), red blood cell count (RBC), hemoglobin (HGB), platelets (PLT) and transformed (log10) egg per gram of feces (EPGlog). A total of 576 animals were genotyped with the Ovine SNP12k BeadChip (Illumina, Inc.), that contains 12,785 bialleleic SNP markers. The variance components were estimated using a single trait model by single step genomic BLUP procedure.
RESULTS: The overall linkage disequilibrium (LD) mean between pairs of markers measured by r2 was 0.23. The overall LD mean between markers considering windows up to 10 Mb was 0.07. The mean LD between adjacent SNPs across autosomes ranged from 0.02 to 0.10. Heritability estimates were low for EPGlog (0.11), moderate for RBC (0.18), PLT (0.17) HCT (0.20), HGB (0.16) and WBC (0.22), and high for FAM (0.35). A total of 22, 21, 23, 20, 26, 25 and 23 windows for EPGlog for FAM, WBC, RBC, PLT, HCT and HGB traits were identified, respectively. Among the associated windows, 10 were shown to be common to HCT and HGB traits on OAR1, OAR2, OAR3, OAR5, OAR8 and OAR15.
CONCLUSIONS: The traits indicating gastrointestinal parasites resistance presented an adequate genetic variability to respond to selection in Santa Inês breed, and it is expected a higher genetic gain for FAM trait when compared to the others. The level of LD estimated for markers separated by less than 1 Mb indicated that the Ovine SNP12k BeadChip might be a suitable tool for identifying genomic regions associated with traits related to gastrointestinal parasite resistance. Several candidate genes related to immune system development and activation, inflammatory response, regulation of lymphocytes and leukocytes proliferation were found. These genes may help in the selection of animals with higher resistance to parasites.