You are currently viewing a beta version of our website. If you spot anything unusual, kindly let us know.

Preprint

Article

Genomic regions associated with resistance to gastrointestinal parasites in Australian Merino sheep

Altmetrics

Downloads

109

Views

Comments

A peer-reviewed article of this preprint also exists.

This version is not peer-reviewed

Abstract

The objective of this study was to identify genomic regions and genes associated with resistance to gastrointestinal nematodes in Australian Merino sheep in Uruguay, using the single-step GWAS methodology (ssGWAS), which is based on genomic estimated breeding values (GEBVs) obtained from a combination of pedigree, genomic and phenotypic data. This methodology converts GEBVs into SNP effects. Analysis included 26,638 animals with fecal egg count (FEC) records obtained in two independent parasitic cycles (FEC1 and FEC2), and 1,700 50K SNP genotypes. Comparison of genomic regions was based on genetic variances (gVar(%)) explained by non-overlapping regions of 20 SNPs. For FEC1 and FEC2, 18 and 22 genomic windows exceeded the significance threshold (gVar(%) ≥ 0.22%), respectively. The genomic regions with strong associations with FEC1 were located on chromosomes OAR 2, 6, 11, 21 and 25, and for FEC2 on OAR 5, 6 and 11. The proportion of genetic variance attributed to the top windows was 0.83% and 1.9% for FEC1 and FEC2, respectively. The 33 candidate genes shared between the two traits were subjected to enrichment analysis, revealing a marked enrichment in biological processes related to immune system functions. These results contribute to the understanding of the genetics underlying gastrointestinal parasite resistance and its implications for other productive and welfare traits in animal breeding programs.

Keywords:

Subject: Biology and Life Sciences - Animal Science, Veterinary Science and Zoology

1. Introduction

The Australian Merino breed is a wool breed traditionally bred in Uruguay in extensive production systems, mainly located in the basalt region. Currently, this breed represents about 40% of the country’s sheep stock and is the main breed in fine and superfine wool production. One of the main problems affecting sheep production is infection by gastrointestinal parasites (GIP) [1]. The decrease in production is a consequence of growth retardation, decreased weight gain, fleece weight and wool quality, and increased mortality [2,3], with Haemonchus contortus and Trichostrongylus colubriformis being the most prevalent parasites in Uruguay [4].

Given the problem of parasitism and the reported anthelmintic resistance [5,6,7], one alternative control method involves the breeding of animals genetically resistant to GIP. The selection criterion used to assess resistance is the parasite fecal egg count (FEC), which is a moderately heritable trait [8] included in the National Genetic Evaluations (NGE) in Uruguay. Currently, the NGE protocol includes the recording of two counts: the initial FEC (FEC1) shortly after weaning (7-9 months of age), and the subsequent FEC (FEC2) at 10-14 months of age [3].

Genome-wide association studies (GWAS) have been implemented in sheep populations to identify genes potentially associated with GIP resistance [9,10,11,12,13]. Unlike other methodologies that are limited to using genomic information exclusively from phenotyped animals, ssGWAS integrates genomic, pedigree and phenotypic data from both genotyped and non-genotyped animals [14], thus considering population structure [15]. The ssGWAS procedure combines traditional pedigree relationships with those derived from genetic markers, and by conversion of GEBVs to marker effects and weights[14] .

The objectives of our study were: (1) to estimate variance and heritability components for FEC1 and FEC2 traits in the Australian Merino sheep population in Uruguay; (2) to identify genomic regions and candidate genes associated with each trait by ssGWAS; and (3) to explore genes associated with FEC1 and FEC2, providing insight into the biological mechanisms underlying resistance to GIP in sheep.

2. Materials and Methods

2.1. Phenotypic Data

FEC determination was conducted according to the current protocol used for recording phenotypes in the NGE [3]. Fecal samples were collected from two independent natural parasitic cycles separated by an anthelmintic treatment (Figure 1). Samples were collected post-weaning, the mean age at recording (days) for FEC1 and FEC2 was 273 (± 69) and 341 (± 62), respectively. This scheme allowed collecting samples in two seasons with climatic conditions that favor the prevalence of different GPI. A fecal sample was obtained from each individual directly from the rectum and the modified Mc Master technique with a sensitivity of 100 eggs per gram of feces was used to estimate the FEC [16].

To establish the optimal time for individual FEC1 and FEC2 sampling, FEC was monitored in a sample of animals until an average load of 500 was obtained. Counts were transformed to natural logarithm as described by Ciappesoni et al. [8] due to their non-normal distribution (LogFEC = Log_e (FEC+100)). In this study, we refer to logFEC1 and logFEC2 as FEC1 and FEC2, respectively.

A total of 26,638 animals born between 2001 and 2020 and belonging to 13 farms had FEC1 records. Among these, 18,971 animals also had FEC2 records (Table 1).

2.2. Genomic and Pedigree Data

Genomic DNA was extracted from blood samples as described by Carracelas et al. [10] . A total of 1702 individuals were genotyped with the GeneSeek® Genomic Profiler™ Ovine 50K panel (GGP, 43,705 SNPs). Genomic data quality control was performed using preGSf90 (Aguilar et al., 2014). SNPs with call rate below 90%, with minor allele frequency (MAF) less than 1%, monomorphic SNPs and animals with call rate less than 90% were removed. Finally, 38,268 SNPs for 1697 sheep were used in the analysis.

The pedigree file was corrected using SeekParentF90 [17], which detects incompatibilities based on Mendelian conflict counts, as described in Wiggans et al. [18].

2.3. Statistical Analysis

Genetic parameters were estimated using methods based on pedigree relationships (BLUP) and pedigree-genomic models (ssGBLUP) [19]. A univariate model was conducted to estimate the variance components and heritabilities for FEC1 and FEC2 employing the AIREMLF90 software from the BLUPF90 family of programs [17]. Additionally, genetic correlations between FEC1 and FEC2 were estimated.

The following univariate model was used:

y = X b + Z u + e

Where y is the vector of phenotypes for FEC1 or FEC2; X and Z are incidence matrices for fixed and random effects, respectively; b is the vector of fixed effects including 494 contemporary groups (year of birth, sex, flock), dam age (three levels: 2, 3, and > 4 years), type of birth (two levels: single or multiple), and lamb age at FEC1 or FEC2 recording as a covariate; u is the vector of random additive genetic effects, e is the vector of residual effects.

In BLUP estimates, the random effects were modeled as u ~ N(0, A

σ_{a}^{2}

) and e ~ N(0, I

σ_{e}^{2}

), where

A

represents the numerator relationship matrix, I is the identity matrix,

σ_{a}^{2}

stands for the additive genetic variance, and

σ_{e}^{2}

is the residual variance. In the ssGBLUP model, the numerator relationship matrix (

A^{- 1}

) utilized in BLUP is substituted with the

H^{- 1}

matrix.

H^{- 1} = A^{- 1} + [\begin{matrix} 0 & 0 \\ 0 & G^{- 1} - A_{22}^{- 1} \end{matrix}]

Where

A^{- 1}

and

A_{22}^{- 1}

are the inverse of the pedigree relationship matrix for all animals and for genotyped animals, respectively, and

G^{- 1}

is the inverse of the genomic relationship matrix. The matrix

G

was constructed as described by Van Raden [20]:

G = Z D Z^{'} q

Where

Z

is the SNP incidence matrix adjusted for allele frequencies,

D

is a weight matrix for SNP (initially D = I), and q is a weighting factor. The SNP effects and weighting factor were derived using an iterative process described by Wang et al. [14]. In this study, a single iteration was used as there was no significant change in SNP effects with additional iterations.

The percentage of genetic variance explained by region was calculated as:

\frac{V a r (a_{i})}{σ_{a}^{2}} x 100 % = \frac{V a r (Σ_{j = 1}^{20} Z_{j} {\hat{u}}_{j})}{σ_{a}^{2}} x 100 %

where

a_{i}

is the genetic value of the i-th region (20 contiguous SNPs),

σ_{a}^{2}

is the total genetic variance,

Z_{j}

is the gene content vector of the j-th SNP for all individuals, and

{\hat{u}}_{j}

is the effect of the j-th SNP marker within the i-th region [21].

Variance components and heritabilities were estimated by AIREMLF90. Heritability (

h^{2}

) was calculated as

h^{2} = \frac{σ_{a}^{2}}{σ_{p}^{2}}

and the total phenotypic variance (

σ_{p}^{2}

) as the sum of the additive genetic variance (

σ_{a}^{2}

) and the residual variance (

σ_{e}^{2}

2.4. Single-Step GWAS Analysis

ssGWAS is a two-step iterative procedure: 1) prediction of GEBVs using ssGBLUP, and 2) prediction of SNP effects based on GEBV. The detailed algorithm was described by Wang et al. [14].

The ssGWAS analysis was done independently for each trait and was done sequentially using RENUMF90 for general dataset preparation, PREGSF90 for quality control and generation of clean genotypes, BLUPF90, and POSTGSF90 for the prediction of breeding values and SNP effects, respectively. These programs are part of the BLUPF90 software family, and for this study, the step-by-step tutorial reported by Masuda [22] was followed.

2.5. Identification of Candidate Genes and Functional Enrichment Analysis

Assuming that all windows explain the same proportion of genetic variance, the proportion of genetic variance explained by each of the 2519 1-Mb windows, including the 38,268 SNPs in the sheep genome, was 0.04%. Therefore, windows that explained at least 0.22% of the genetic variance, which is 5 times higher than expected (0.045%*5 = 0.22%), were considered to contain putative QTL [23,24]. Regions representing 0.22% or more of the genetic variance

σ_{a}^{2}

were defined as significant regions. SNPs within these regions were identified and mapped onto the Oar_v3.1 sheep genome assembly using the Ensembl database [25]. A range of 5kb upstream and downstream of the variant position were considered to identify candidate genes.

Functional enrichment analysis was conducted on the set of common candidate genes for both traits using DAVID (https://david.ncifcrf.gov/tools.jsp accessed on 10 February 2024) [26]. Gene Ontology (GO) terms with a p-value ≤ 0.05 were reported as significant terms.

3. Results

3.1. Variance Components and Heritabilities

The variance components and heritabilities for GIP resistance in Australian Merino sheep are presented in Table 2. BLUP and ssGBLUP estimates of the heritabilities for FEC1 and FEC2 were close to 0.19. On the other hand, the genetic correlation between both traits was 0.88 (±0.03).

3.2. Genome-Wide Association Analysis

Figure 2 shows the Manhattan plots for GIP resistance traits. The genomic regions that explained the largest genetic variance for FEC1 were located on chromosomes 2, 6, 11, 21 and 25, and for FEC2 they were on chromosomes 5, 6 and 11. The top windows (most significant) explained 0.9% and 2% of the genetic variance for FEC1 and FEC2, respectively (Table 3).

Eighteen and 22 windows with genetic variance greater than the significance threshold of 0.22 gVar(%) were identified, which included 316 and 376 SNPs related to FEC1 and FEC2 traits, respectively. Positional candidate genes close to SNPs (≤ 0.5 Mbps) were identified using the Ovis aries 3.1. reference genome map and a total of 67 and 63 genes were mapped for FEC 1 and FEC 2, respectively. In total, 33 genes were shared between both traits (Figure 3). Details of the common positional candidate genes for the FEC 1 and FEC 2 traits are summarized in Table 4.

3.3. Enrichment Analysis

The ssGWAS results were complemented with a gene ontology (GO) enrichment analysis, which revealed 17 significantly enriched GO terms (p ≤ 0.05). Among these, ten were related to biological processes, two to molecular functions and five to cellular components. The KEGG pathway enrichment analysis revealed three enriched metabolic pathways for the set of analyzed genes (Figure 4). Details of the enriched GO categories and the metabolic pathways involved are shown in Table 5.

4. Discussion

Variance component estimates obtained using the traditional pedigree-based approach (BLUP) were like those obtained using the ssGBLUP procedure, as well as the estimated heritability values (0.19 vs. 0.20). Medium to low heritabilities for resistance to GIP in sheep [27,28] and for Australian Merino [8] have already been reported. Genetic parameter estimates obtained using ssGBLUP are known to be less biased and more accurate [29,30,31] since relationships between animals are better estimated [32]. In our case, the results were similar and because there were no changes in the genetic base, the same additive variance is expected when including the genomic coefficients, as reported by Forni et al. [29].

The strong genetic correlation between FEC1 and FEC2 (0.88) suggests that they can be considered as the same trait genetically, even when these traits were measured at different ages and correspond to two different parasitic cycles and seasons of the year (Figure 1), in which animals could have been exposed to different parasites. This high estimate is in agreement with other studies that explored the genetic association between FEC recorded at different ages in different breeds (i.e. 0.85 Romney; 0.82 Katahdin) [33,34]. A high genetic correlation of 0.74 between FEC by Strongyles and Nematodirus was also reported by Pacheco et al. [35].

On the other hand, ssGWAS revealed that the genomic regions reported as significant explain only 7 and 6% of the genetic variance for FEC1 and FEC2, respectively. These small variances suggest that resistance to GIP is a polygenic trait with a large number of variants involved in the resistance mechanism [10,36,37]. Significant regions associated with GIP resistance have been previously reported on chromosomes 2, 6, 18 and 24 in several Australian sheep populations including the Merino breed [38] and GIP-related QTL regions are also known [39].

In this study, seventeen positional candidate genes were identified on the OAR 2 for the FEC1 trait: AMER3, BMP1, CCAR2, CYFIP1, FAM160B2, HERC2, HR, NIPA1, NIPA2, PDLIM2, PIWIL2, POLR3D, PPP3CC, PTPN18, SORBS3, TUBGCP5, XPO7. Among these, several are involved in the mechanisms of innate and adaptive immunity in mammals, such as HERC2 and CYFIP1 that are also involved in cytokine signaling [38,40]. The PDLIM2 gene has been associated in transcriptomic studies with the immune system and reproduction in sheep [41]. PTPN18 is involved in the B-cell receptor signaling pathway, being involved in differentiation, proliferation and immunoglobulin (Ig) production [38], while other reports relate it to pigmentation in Merino sheep [42]. The SORBS3, PPP3CC and PIWIL2 genes have been linked to growth and wool quality traits [43]. PPP3CC has been associated with heat tolerance in cattle [44], while PIWIL2 has been linked to reproductive traits in pigs [45]. In addition, these genes have been reported in selection signature studies, revealing their involvement in parasite resistance of Slovakian sheep populations [46]. Al Kalaldeh et al. [38] also reported the association of the HERC2, NIPA1, NIPA2, CYFIP1, TUBGCP5, PTPN18 and AMER3 genes with GIP resistance in sheep and their involvement in immune system mechanisms. In addition, the CAST gene (OAR 5) was identified in the ssGWAS for FEC2 and it is known to have relevance to traits such as muscle production, carcass and meat quality in sheep [47,48].

On the other hand, 12 candidate genes significantly associated with the FEC1 or FEC2 traits were identified in OAR 6: DCAF16, FAM184B, GPRIN3, HERC3, HERC5, IBSP, LAP3, LCORL, MED28, NCAPG, PIGY and PYURF, all of them except LAP3 were detected by Al Kalaldeh et al. [38] in significant genomic regions for GIP resistance in sheep. Particularly the FAM184B gene has also been reported to be strongly associated with body size and weight in livestock [49,50,51], milk production [52] and reproductive traits in sheep [53].

The MED28 gene is involved in milk production in ewes [43,52] and is also related to liveweight [54]. In addition, this gene has been linked to both pre and postnatal body weights in ewes [49,54]. On the other hand, LAP3 also contributes to growth, milk production and feed efficiency traits in sheep [55]. The GPRIN3 gene has been linked to prolificacy, litter size [56] and temperament in sheep [57].

Five candidate genes were identified on OAR 8: FAM120B, SMOC2, TBP, THBS2, WDR27, and WDR27. The THBS2 gene is related to growth and development traits [58] and high prolificacy in sheep [59], while the WDR27 gene has been associated with the cow milk fat trait [60].

On chromosome 11, a total of 43 candidate genes for FEC1 and FEC2 were identified, including ACADVL, ACAP1, ALOX12, ALOX15, ARRB2, ASGR1, BCL6B, C11H17orf49, C1QBP, AMTA2, CCDC42, CD68, CXCL16, DLG4, DNAH2, DVL2, INCA1, KIF1C, LOC101115381, LOC101121185, MED11, MINK1, MPDU1, MYH10, NDEL1, NUP88, ODF4, PELP1, PIK3R5, PLD2, RABEP1, RPAIN, SENP3, SLC16A13, SLC25A35, SPAG7, STX8, TNFSF12, USP43, WRAP53, WSCD1, ZFP3 and ZMYND15. The ALOX15 gene is known to play a role in Haemonchus contortus infection in sheep [61]. Other studies have reported QTLs for Trichostrongylus colubriformis egg count in the Merino breed in this genomic region [62], as well as QTLs for height, weight and parasite resistance in the region between 32.13 and 32.19 Mbp [63]. In addition, selection signatures for resistance to Haemonchus contortus in sheep and goats have been reported by Estrada-Reyes et al. [64]. While the CTNNA gene, located on chromosome 25, was significantly associated with FEC1 and FEC2 traits and is related to growth traits [65] and brucellosis resistance in sheep [66].

The GO analysis (Figure 4) revealed multiple GO terms related to host defense mechanisms against pathogen including lipoxin A4 biosynthetic process (GO:2001303). Lipoxins have been reported to be endogenous anti-inflammatory molecules involved in reducing excessive tissue damage and chronic inflammation [67]. These lipoxins can be synthesized from arachidonic acid, which justifies the overrepresentation of biological processes such as lipoxygenase pathway (GO:0019372) and arachidonic acid metabolic process (GO:0019369), as well as arachidonate 15-lipoxygenase activity and linoleate 13S-lipoxygenase activity molecular functions (GO:0050473, GO:001616).

Cytokines, such as IL-12, are critical for host resistance to many pathogens but can also be detrimental when expressed in an uncontrolled manner [68,69], so it makes biological sense that GO term negative regulation of interleukin-12 production (GO:0032695) is overrepresented. The role of lipoxins in mediating the immune response has been studied against parasitic pathogens [70] and in other diseases [67]. Another GO term related to host defense against pathogens is the GO term linoleic acid metabolic process (GO:0043651). The activation of linoleic acid metabolism in macrophages in bacterias promotes pathogen killing [71]. In that sense, one of the most overrepresented genes is the ALOX15 gene, which is related to the oxidation of arachidonic acid and the production of anti-inflammatory lipoxins [72] and it was identified as a differentially expressed gene associated with sheep resistance to the nematode Teladorsagia circumcincta [73].

As with all infections, parasites cause inflammation, involving mobilization, proliferation and recruitment of leukocytes to the affected area. This trafficking of immune cells and the effector functions of these cells can serve to control the pathogen or exacerbate the pathology [74]. In this regard, one of the most overrepresented metabolic pathways is the oas04062: Chemokine signaling pathway. Two of the functions of Chemokines are to attract immune cells to sites of inflammation [75] and to guide the migration of neurons and other migratory cells [76].

In summary, the enrichment analysis shows that the candidate genes enriched biological processes and molecular functions related to the metabolism of linoleic acid, which is a metabolic precursor of arachidonic acid. Free arachidonic acid and its metabolites promote and modulate the type 2 immune response, playing a crucial role in GIP resistance through the action of eosinophils, basophils and mast cells [77].

5. Conclusions

In this study, genetic parameters were estimated to evaluate resistance to GIP in two independent parasite cycles, each separated by anthelmintic treatment. A total of 18 and 22 genomic regions were identified that showed significant association with FEC1 and FEC2, respectively. We report positional candidate genes for both cycles using ssGWAS in Australian Merino sheep, some of which are novel for these traits. Our study reveals a set of candidate genes that share mechanisms related to immune response, body size and weight, as well as genes associated with reproductive traits. In summary, our findings provide a basis for future genomic research and could contribute significantly to breeding programs.

Author Contributions

B.V.: Investigation, Data Curation, Formal analysis, Writing—Original draft, Writing—Review & Editing, B.C.: Resources, Methodology, Writing—Review & Editing, P.P.: Resources, Methodology, Writing—Review & Editing, E.N.: Conceptualization, Methodology, Investigation, Writing—Review & Editing, Project administration, Funding acquisition, E. VL.: Writing—Review & Editing, Funding acquisition, G.C.: Conceptualization, Writing—Review & Editing, Project administration, Funding acquisition. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Instituto Nacional de Investigación Agropecuaria (INIA_CL_40 and INIA_CL_38), CSIC_I+D_2018_287 (Comisión Sectorial de Investigación Científica, CSIC, Universidad de la República) and the European Union’s Horizon 2020 research and innovation program (SMARTER, agreement no. 772787).

Institutional Review Board Statement

The animal study protocol was approved by the Animal Ethics Committee of INIA (approval number INIA_2018.2, July 2018) and approved by the Honorary Commission for Animal Experimentation (CHEA, protocols N° 021130-006307-11 and 020300-000653-18-ID, CHEA 701, Universidad de la República).

Informed Consent Statement

Not applicable.

Data Availability Statement

Restrictions apply to the availability of these data. Data were obtained by INIA and are available from the authors with INIA’s permission.

Acknowledgments

We thank Ignacio Aguilar for his help using BLUPF90 programs, the Australian Merino Breeders Association of Uruguay and Estación Experimental Facultad de Agronomía Salto and to the INIA staff who collaborated in collecting data for this research, especially to Carlos Monzalvo for managing the phenotypic database.

Conflicts of Interest

The authors declare no conflicts of interest.

References

Zajac, A.M. Gastrointestinal Nematodes of Small Ruminants: Life Cycle, Anthelmintics, and Diagnosis. Vet. Clin. North Am. Food Anim. Pract. 2006, 22, 529–541. [Google Scholar] [CrossRef] [PubMed]
Zajac, A.M.; Garza, J. Biology, Epidemiology, and Control of Gastrointestinal Nematodes of Small Ruminants. Vet. Clin. North Am. Food Anim. Pract. 2020, 36, 73–87. [Google Scholar] [CrossRef] [PubMed]
Castells, D. Evaluación de Resistencia de Ovinos Corriedale a Los Nematodos En Uruguay: Heredabilidad y Correlaciones Genéticas El Recuento De Huevos De Nematodos Características Productivas. Tesis de maestría, UDELAR: Montevideo, Uruguay, 2008.
Nari, A.; Cardozo, H.; Berdié, J.; Canábez, F.; Bawden, R. Dinámica de población para nematodes gastrointestinales de ovinos en Uruguay: Trabajo asistido por UNDP y FAO. Vet. Montev. 1977, 14, 11–23. [Google Scholar]
Mederos, A.E.; Ramos, Z.; Banchero, G.E. First Report of Monepantel Haemonchus Contortus Resistance on Sheep Farms in Uruguay. Parasit. Vectors 2014, 7, 598. [Google Scholar] [CrossRef] [PubMed]
Roeber, F.; Jex, A.R.; Gasser, R.B. Impact of Gastrointestinal Parasitic Nematodes of Sheep, and the Role of Advanced Molecular Tools for Exploring Epidemiology and Drug Resistance - an Australian Perspective. Parasit. Vectors 2013, 6, 153. [Google Scholar] [CrossRef]
Waller, P.J.; Echevarria, F.; Eddi, C.; Maciel, S.; Nari, A.; Hansen, J.W. The Prevalence of Anthelmintic Resistance in Nematode Parasites of Sheep in Southern Latin America: General Overview. Vet. Parasitol. 1996, 62, 181–187. [Google Scholar] [CrossRef]
Ciappesoni, G.; Goldberg, V.; Gimeno, D. Estimates of Genetic Parameters for Worm Resistance, Wool and Growth Traits in Merino Sheep of Uruguay. Livest. Sci. 2013, 157, 65–74. [Google Scholar] [CrossRef]
Arzik, Y.; Kizilaslan, M.; White, S.N.; Piel, L.M.W.; Çınar, M.U. Genomic Analysis of Gastrointestinal Parasite Resistance in Akkaraman Sheep. Genes 2022, 13, 2177. [Google Scholar] [CrossRef]
Carracelas, B.; Navajas, E.A.; Vera, B.; Ciappesoni, G. Genome-Wide Association Study of Parasite Resistance to Gastrointestinal Nematodes in Corriedale Sheep. Genes 2022, 13, 1548. [Google Scholar] [CrossRef]
Estrada-Reyes, Z.M.; Ogunade, I.M.; Pech-Cervantes, A.A.; Terrill, T.H. Copy Number Variant-Based Genome Wide Association Study Reveals Immune-Related Genes Associated with Parasite Resistance in a Heritage Sheep Breed from the United States. Parasite Immunol. 2022, 44, e12943. [Google Scholar] [CrossRef]
Pacheco, A.; Conington, J.; Corripio-Miyar, Y.; Frew, D.; Banos, G.; McNeilly, T.N. Genetic Profile of Adaptive Immune Traits and Relationships with Parasite Resistance and Productivity in Scottish Blackface Sheep. animal 2024, 18, 101061. [Google Scholar] [CrossRef]
Thorne, J.W.; Redden, R.; Bowdridge, S.A.; Becker, G.M.; Stegemiller, M.R.; Murdoch, B.M. Genome-Wide Analysis of Sheep Artificially or Naturally Infected with Gastrointestinal Nematodes. Genes 2023, 14, 1342. [Google Scholar] [CrossRef]
Wang, H.; Misztal, I.; Aguilar, I.; Legarra, A.; Muir, W.M. Genome-Wide Association Mapping Including Phenotypes from Relatives without Genotypes. Genet. Res. 2012, 94, 73–83. [Google Scholar] [CrossRef]
Mancin, E.; Lourenco, D.; Bermann, M.; Mantovani, R.; Misztal, I. Accounting for Population Structure and Phenotypes From Relatives in Association Mapping for Farm Animals: A Simulation Study. Front. Genet. 2021, 12. [Google Scholar] [CrossRef] [PubMed]
Whitlock, H.V. Some Modifications of the McMaster Helminth Egg-Counting Technique and Apparatus. 1948, 21, 177–180. [Google Scholar]
Misztal, I.; Tsuruta, S.; Lourenco, D.; Masuda, Y.; Aguilar, I. BLUPF90 Family of Programs Available online: University of Georgia.
Wiggans, G.R.; VanRaden, P.M.; Bacheller, L.R.; Tooker, M.E.; Hutchison, J.L.; Cooper, T.A.; Sonstegard, T.S. Selection and Management of DNA Markers for Use in Genomic Evaluation. J. Dairy Sci. 2010, 93, 2287–2292. [Google Scholar] [CrossRef] [PubMed]
Aguilar, I.; Misztal, I.; Johnson, D.L.; Legarra, A.; Tsuruta, S.; Lawlor, T.J. Hot Topic: A Unified Approach to Utilize Phenotypic, Full Pedigree, and Genomic Information for Genetic Evaluation of Holstein Final Score1. J. Dairy Sci. 2010, 93, 743–752. [Google Scholar] [CrossRef] [PubMed]
VanRaden, P.M. Efficient Methods to Compute Genomic Predictions. J. Dairy Sci. 2008, 91, 4414–4423. [Google Scholar] [CrossRef] [PubMed]
Wang, H.; Misztal, I.; Aguilar, I.; Legarra, A.; Fernando, R.L.; Vitezica, Z.; Okimoto, R.; Wing, T.; Hawken, R.; Muir, W.M. Genome-Wide Association Mapping Including Phenotypes from Relatives without Genotypes in a Single-Step (ssGWAS) for 6-Week Body Weight in Broiler Chickens. Front. Genet. 2014, 5. [Google Scholar] [CrossRef]
Yutaka Masuda Available online:. Available online: https://masuday.github.io/ (accessed on 12 February 2024).
Sollero, B.P.; Junqueira, V.S.; Gomes, C.C.G.; Caetano, A.R.; Cardoso, F.F. Tag SNP Selection for Prediction of Tick Resistance in Brazilian Braford and Hereford Cattle Breeds Using Bayesian Methods. Genet. Sel. Evol. 2017, 49, 49. [Google Scholar] [CrossRef]
Onteru, S.K.; Gorbach, D.M.; Young, J.M.; Garrick, D.J.; Dekkers, J.C.M.; Rothschild, M.F. Whole Genome Association Studies of Residual Feed Intake and Related Traits in the Pig. PLOS ONE 2013, 8, e61756. [Google Scholar] [CrossRef] [PubMed]
International Sheep Genomics Consortium; Archibald, A.L.; Cockett, N.E.; Dalrymple, B.P.; Faraut, T.; Kijas, J.W.; Maddox, J.F.; McEwan, J.C.; Hutton Oddy, V.; Raadsma, H.W.; et al. The Sheep Genome Reference Sequence: A Work in Progress. Anim. Genet. 2010, 41, 449–453. [Google Scholar] [CrossRef] [PubMed]
Sherman, B.T.; Hao, M.; Qiu, J.; Jiao, X.; Baseler, M.W.; Lane, H.C.; Imamichi, T.; Chang, W. DAVID: A Web Server for Functional Enrichment Analysis and Functional Annotation of Gene Lists (2021 Update). Nucleic Acids Res. 2022, 50, W216–W221. [Google Scholar] [CrossRef] [PubMed]
Gutiérrez-Gil, B.; Pérez, J.; de la Fuente, L.F.; Meana, A.; Martínez-Valladares, M.; San Primitivo, F.; Rojo-Vázquez, F.A.; Arranz, J.J. Genetic Parameters for Resistance to Trichostrongylid Infection in Dairy Sheep. Animal 2010, 4, 505–512. [Google Scholar] [CrossRef] [PubMed]
Safari, E.; Fogarty, N.M.; Gilmour, A.R. A Review of Genetic Parameter Estimates for Wool, Growth, Meat and Reproduction Traits in Sheep. Livest. Prod. Sci. 2005, 92, 271–289. [Google Scholar] [CrossRef]
Forni, S.; Aguilar, I.; Misztal, I. Different Genomic Relationship Matrices for Single-Step Analysis Using Phenotypic, Pedigree and Genomic Information. Genet. Sel. Evol. 2011, 43, 1. [Google Scholar] [CrossRef] [PubMed]
Gordo, D.G.M.; Espigolan, R.; Tonussi, R.L.; Júnior, G.A.F.; Bresolin, T.; Magalhães, A.F.B.; Feitosa, F.L.; Baldi, F.; Carvalheiro, R.; Tonhati, H.; et al. Genetic Parameter Estimates for Carcass Traits and Visual Scores Including or Not Genomic Information1. J. Anim. Sci. 2016, 94, 1821–1826. [Google Scholar] [CrossRef] [PubMed]
Tsuruta, S.; Misztal, I.; Lourenco, D. a. L.; Lawlor, T.J. Assigning Unknown Parent Groups to Reduce Bias in Genomic Evaluations of Final Score in US Holsteins. J. Dairy Sci. 2014, 97, 5814–5821. [Google Scholar] [CrossRef]
Lourenco, D.A.L.; Tsuruta, S.; Fragomeni, B.O.; Masuda, Y.; Aguilar, I.; Legarra, A.; Bertrand, J.K.; Amen, T.S.; Wang, L.; Moser, D.W.; et al. Genetic Evaluation Using Single-Step Genomic Best Linear Unbiased Predictor in American Angus1. J. Anim. Sci. 2015, 93, 2653–2662. [Google Scholar] [CrossRef]
Morris, C.; Wheeler, M. Correlated Responses Following Genetic Selection to Change Faecal Egg Count in Romney. Proc. N. Z. Soc. Anim. Prod. 2010, 70, 229–234. [Google Scholar]
Ngere, L.; Burke, J.M.; Morgan, J.L.M.; Miller, J.E.; Notter, D.R. Genetic Parameters for Fecal Egg Counts and Their Relationship with Body Weights in Katahdin Lambs. J. Anim. Sci. 2018, 96, 1590–1599. [Google Scholar] [CrossRef] [PubMed]
Pacheco, A.; McNeilly, T.N.; Banos, G.; Conington, J. Genetic Parameters of Animal Traits Associated with Coccidian and Nematode Parasite Load and Growth in Scottish Blackface Sheep. Animal 2021, 15, 100185. [Google Scholar] [CrossRef] [PubMed]
Kemper, K.E.; Emery, D.L.; Bishop, S.C.; Oddy, H.; Hayes, B.J.; Dominik, S.; Henshall, J.M.; Goddard, M.E. The Distribution of SNP Marker Effects for Faecal Worm Egg Count in Sheep, and the Feasibility of Using These Markers to Predict Genetic Merit for Resistance to Worm Infections. Genet. Res. 2011, 93, 203–219. [Google Scholar] [CrossRef]
Riggio, V.; Matika, O.; Pong-Wong, R.; Stear, M.J.; Bishop, S.C. Genome-Wide Association and Regional Heritability Mapping to Identify Loci Underlying Variation in Nematode Resistance and Body Weight in Scottish Blackface Lambs. Heredity 2013, 110, 420–429. [Google Scholar] [CrossRef]
Al Kalaldeh, M.; Gibson, J.; Lee, S.H.; Gondro, C.; van der Werf, J.H.J. Detection of Genomic Regions Underlying Resistance to Gastrointestinal Parasites in Australian Sheep. Genet. Sel. Evol. GSE 2019, 51, 37. [Google Scholar] [CrossRef]
Cunha, S.M.F.; Lam, S.; Mallard, B.; Karrow, N.A.; Cánovas, Á. Genomic Regions Associated with Resistance to Gastrointestinal Nematode Parasites in Sheep—A Review. Genes 2024, 15, 187. [Google Scholar] [CrossRef] [PubMed]
Abied, A.; Bagadi, A.; Bordbar, F.; Pu, Y.; Augustino, S.M.A.; Xue, X.; Xing, F.; Gebreselassie, G.; Han, J.-L.; Mwacharo, J.M.; et al. Genomic Diversity, Population Structure, and Signature of Selection in Five Chinese Native Sheep Breeds Adapted to Extreme Environments. Genes 2020, 11, 494. [Google Scholar] [CrossRef]
Pokharel, K.; Peippo, J.; Weldenegodguad, M.; Honkatukia, M.; Li, M.-H.; Kantanen, J. Transcriptome Analysis Reveals the Importance of the Immune System during Early Pregnancy in Sheep (Ovis Aries) 2019, 673558.
García-Gámez, E.; Reverter, A.; Whan, V.; McWilliam, S.M.; Arranz, J.J.; Kijas, J. Using Regulatory and Epistatic Networks to Extend the Findings of a Genome Scan: Identifying the Gene Drivers of Pigmentation in Merino Sheep. PLoS ONE 2011, 6, e21158. [Google Scholar] [CrossRef]
Mohammadi, H.; Rafat, S.A.; Moradi Shahrbabak, H.; Shodja, J.; Moradi, M.H. Genome-Wide Association Study and Gene Ontology for Growth and Wool Characteristics in Zandi Sheep. J. Livest. Sci. Technol. 2020, 8, 45–55. [Google Scholar] [CrossRef]
Li, R.; Li, C.; Chen, H.; Li, R.; Chong, Q.; Xiao, H.; Chen, S. Genome-Wide Scan of Selection Signatures in Dehong Humped Cattle for Heat Tolerance and Disease Resistance. Anim. Genet. 2020, 51, 292–299. [Google Scholar] [CrossRef]
Fischer, D.; Laiho, A.; Gyenesei, A.; Sironen, A. Identification of Reproduction-Related Gene Polymorphisms Using Whole Transcriptome Sequencing in the Large White Pig Population. G3 GenesGenomesGenetics 2015, 5, 1351–1360. [Google Scholar] [CrossRef] [PubMed]
Mészárosová, M.; Mészáros, G.; Moravčíková, N.; Pavlík, I.; Margetín, M.; Kasarda, R. Within- and between-Breed Selection Signatures in the Original and Improved Valachian Sheep. Animals 2022, 12, 1346. [Google Scholar] [CrossRef] [PubMed]
Greguła-Kania, M.; Gruszecki, T.M.; Junkuszew, A.; Juszczuk-Kubiak, E.; Florek, M. Association of CAST Gene Polymorphism with Carcass Value and Meat Quality in Two Synthetic Lines of Sheep. Meat Sci. 2019, 154, 69–74. [Google Scholar] [CrossRef] [PubMed]
Sutikno, *!!! REPLACE !!!*; Yamin, M.; Sumantri, C. Association of Polymorphisms Calpastatin Gene with Body Weight of Local Sheep in Jonggol, Indonesia. Media Peternakan Vol 34, No 1 2011. [Google Scholar] [CrossRef]
Al-Mamun, H.A.; Kwan, P.; Clark, S.A.; Ferdosi, M.H.; Tellam, R.; Gondro, C. Genome-Wide Association Study of Body Weight in Australian Merino Sheep Reveals an Orthologous Region on OAR6 to Human and Bovine Genomic Regions Affecting Height and Weight. Genet. Sel. Evol. 2015, 47, 66. [Google Scholar] [CrossRef] [PubMed]
Bolormaa, S.; Hayes, B.J.; van der Werf, J.H.J.; Pethick, D.; Goddard, M.E.; Daetwyler, H.D. Detailed Phenotyping Identifies Genes with Pleiotropic Effects on Body Composition. BMC Genomics 2016, 17, 224. [Google Scholar] [CrossRef] [PubMed]
Yuan, Z.; Sunduimijid, B.; Xiang, R.; Behrendt, R.; Knight, M.I.; Mason, B.A.; Reich, C.M.; Prowse-Wilkins, C.; Vander Jagt, C.J.; Chamberlain, A.J.; et al. Expression Quantitative Trait Loci in Sheep Liver and Muscle Contribute to Variations in Meat Traits. Genet. Sel. Evol. 2021, 53, 8. [Google Scholar] [CrossRef] [PubMed]
Lukic, B.; Curik, I.; Drzaic, I.; Galić, V.; Shihabi, M.; Vostry, L.; Cubric-Curik, V. Genomic Signatures of Selection, Local Adaptation and Production Type Characterisation of East Adriatic Sheep Breeds. J. Anim. Sci. Biotechnol. 2023, 14, 142. [Google Scholar] [CrossRef] [PubMed]
Gong, Y.; Chen, Q.; He, X.; Wang, X.; He, X.; Wang, Y.; Pan, Z.; Chu, M.; Di, R. Association Analyses between Single Nucleotide Polymorphisms in ZFAT, FBN1, FAM184B Genes and Litter Size of Xinggao Mutton Sheep. Anim. Open Access J. MDPI 2023, 13, 3639. [Google Scholar] [CrossRef]
Ramos, Z.; Garrick, D.J.; Blair, H.T.; Vera, B.; Ciappesoni, G.; Kenyon, P.R. Genomic Regions Associated with Wool, Growth and Reproduction Traits in Uruguayan Merino Sheep. Genes 2023, 14, 167. [Google Scholar] [CrossRef]
Yurchenko, A.A.; Deniskova, T.E.; Yudin, N.S.; Dotsev, A.V.; Khamiruev, T.N.; Selionova, M.I.; Egorov, S.V.; Reyer, H.; Wimmers, K.; Brem, G.; et al. High-Density Genotyping Reveals Signatures of Selection Related to Acclimation and Economically Important Traits in 15 Local Sheep Breeds from Russia. BMC Genomics 2019, 20, 294. [Google Scholar] [CrossRef] [PubMed]
Tao, L.; He, X.; Jiang, Y.; Liu, Y.; Ouyang, Y.; Shen, Y.; Hong, Q.; Chu, M. Genome-Wide Analyses Reveal Genetic Convergence of Prolificacy between Goats and Sheep. Genes 2021, 12, 480. [Google Scholar] [CrossRef] [PubMed]
Romaniuk, E.; Vera, B.; Peraza, P.; Ciappesoni, G.; Damián, J.P.; Van Lier, E. Identification of Candidate Genes and Pathways Linked to the Temperament Trait in Sheep. Genes 2024, 15, 229. [Google Scholar] [CrossRef]
Zhu, M.; Yang, Y.; Yang, H.; Zhao, Z.; Zhang, H.; Blair, H.T.; Zheng, W.; Wang, M.; Fang, C.; Yu, Q.; et al. Whole-Genome Resequencing of the Native Sheep Provides Insights into the Microevolution and Identifies Genes Associated with Reproduction Traits. BMC Genomics 2023, 24, 392. [Google Scholar] [CrossRef] [PubMed]
Chen, H.Y.; Shen, H.; Jia, B.; Zhang, Y.S.; Wang, X.H.; Zeng, X.C. Differential Gene Expression in Ovaries of Qira Black Sheep and Hetian Sheep Using RNA-Seq Technique. PLoS ONE 2015, 10, e0120170. [Google Scholar] [CrossRef] [PubMed]
Peters, S.O.; Kızılkaya, K.; Ibeagha-Awemu, E.M.; Sinecen, M.; Zhao, X. Comparative Accuracies of Genetic Values Predicted for Economically Important Milk Traits, Genome-Wide Association, and Linkage Disequilibrium Patterns of Canadian Holstein Cows. J. Dairy Sci. 2021, 104, 1900–1916. [Google Scholar] [CrossRef] [PubMed]
Guo, Z.; González, J.F.; Hernandez, J.N.; McNeilly, T.N.; Corripio-Miyar, Y.; Frew, D.; Morrison, T.; Yu, P.; Li, R.W. Possible Mechanisms of Host Resistance to Haemonchus Contortus Infection in Sheep Breeds Native to the Canary Islands. Sci. Rep. 2016, 6, 26200. [Google Scholar] [CrossRef] [PubMed]
Beh, K.J.; Hulme, D.J.; Callaghan, M.J.; Leish, Z.; Lenane, I.; Windon, R.G.; Maddox, J.F. A Genome Scan for Quantitative Trait Loci Affecting Resistance to Trichostrongylus Colubriformis in Sheep. Anim. Genet. 2002, 33, 97–106. [Google Scholar] [CrossRef]
Niciura, S.C.M.; Benavides, M.V.; Okino, C.H.; Ibelli, A.M.G.; Minho, A.P.; Esteves, S.N.; Chagas, A.C. de S. Genome-Wide Association Study for Haemonchus Contortus Resistance in Morada Nova Sheep. Pathog. Basel Switz. 2022, 11, 939. [Google Scholar] [CrossRef]
Estrada-Reyes, Z.M.; Tsukahara, Y.; Amadeu, R.R.; Goetsch, A.L.; Gipson, T.A.; Sahlu, T.; Puchala, R.; Wang, Z.; Hart, S.P.; Mateescu, R.G. Signatures of Selection for Resistance to Haemonchus Contortus in Sheep and Goats. BMC Genomics 2019, 20, 735. [Google Scholar] [CrossRef]
Zhao, L.; Li, F.; Yuan, L.; Zhang, X.; Zhang, D.; Li, X.; Zhang, Y.; Zhao, Y.; Song, Q.; Wang, J.; et al. Expression of Ovine CTNNA3 and CAP2 Genes and Their Association with Growth Traits. Gene 2022, 807, 145949. [Google Scholar] [CrossRef] [PubMed]
Li, X.; Wu, Q.; Zhang, X.; Li, C.; Zhang, D.; Li, G.; Zhang, Y.; Zhao, Y.; Shi, Z.; Wang, W.; et al. Whole-Genome Resequencing to Study Brucellosis Susceptibility in Sheep. Front. Genet. 2021, 12. [Google Scholar] [CrossRef] [PubMed]
Chandrasekharan, J.A.; Sharma-Walia, N. Lipoxins: Nature’s Way to Resolve Inflammation. J. Inflamm. Res. 2015, 8, 181–192. [Google Scholar] [CrossRef] [PubMed]
Aliberti, J.; Serhan, C.; Sher, A. Parasite-Induced Lipoxin A4 Is an Endogenous Regulator of IL-12 Production and Immunopathology in Toxoplasma Gondii Infection. J. Exp. Med. 2002, 196, 1253–1262. [Google Scholar] [CrossRef]
Eyles, J.L.; Metcalf, D.; Grusby, M.J.; Hilton, D.J.; Starr, R. Regulación Negativa de La Señalización de Interleucina-12 Por El Supresor de La Señalización de Citocina-1 *. J. Biol. Chem. 2002, 277, 43735–43740. [Google Scholar] [CrossRef] [PubMed]
Wei, F.; Gong, W.; Wang, J.; Yang, Y.; Liu, J.; Wang, Y.; Cao, J. Role of the Lipoxin A4 Receptor in the Development of Neutrophil Extracellular Traps in Leishmania Infantum Infection. Parasit. Vectors 2019, 12, 275. [Google Scholar] [CrossRef]
Yan, B.; Fung, K.; Ye, S.; Lai, P.-M.; Xin Wei, Y.; Sze, K.-H.; Yang, D.; Gao, P.; Yi-Tsun Kao, R. Linoleic Acid Metabolism Activation in Macrophages Promotes the Clearing of Intracellular Staphylococcus Aureus. Chem. Sci. 2022, 13, 12445–12460. [Google Scholar] [CrossRef] [PubMed]
Samuelsson, B.; Dahlén, S.E.; Lindgren, J.A.; Rouzer, C.A.; Serhan, C.N. Leukotrienes and Lipoxins: Structures, Biosynthesis, and Biological Effects. Science 1987, 237, 1171–1176. [Google Scholar] [CrossRef] [PubMed]
Wilkie, H.; Xu, S.; Gossner, A.; Hopkins, J. Variable Exon Usage of Differentially-Expressed Genes Associated with Resistance of Sheep to Teladorsagia Circumcincta. Vet. Parasitol. 2015, 212, 206–213. [Google Scholar] [CrossRef]
McGovern, K.E.; Wilson, E.H. Role of Chemokines and Trafficking of Immune Cells in Parasitic Infections. Curr. Immunol. Rev. 2013, 9, 157–168. [Google Scholar] [CrossRef]
Oppenheim, J.J.; Zachariae, C.O.C.; Mukaida, N.; Matsushima, K. Properties of the Novel Proinflammatory Supergene “Intercrine” Cytokine Family. Annu. Rev. Immunol. 1991, 9, 617–648. [Google Scholar] [CrossRef] [PubMed]
Wang, J.; Knaut, H. Chemokine Signaling in Development and Disease. Dev. Camb. Engl. 2014, 141, 4199–4205. [Google Scholar] [CrossRef] [PubMed]
Tallima, H.; El Ridi, R. Arachidonic Acid: Physiological Roles and Potential Health Benefits – A Review. J. Adv. Res. 2017, 11, 33–41. [Google Scholar] [CrossRef]

Figure 1. Sampling scheme for fecal egg counts (FEC) in two independent parasitic cycles and in different seasonal periods.

Figure 2. Manhattan plots depict the genetic variance explained (%) by 20 adjacent SNP windows for FEC1 (a) and FEC2 (b) in Australian Merino sheep. Each dot represents a window, with the percentage of additive genetic variance explained by each window. The horizontal line indicates the suggestive threshold of 0.22 of the gVar(%).

Figure 3. List of positional candidate genes identified in genomic regions that explain >0.22 gVar(%) for FEC1 and FEC2 (a). The black dot represents that this gene was identified in the analysis for each trait. The bar plot shows the number of candidate genes in common for both traits at the intersection (b).

Figure 4. Gene ontology (GO) term enrichment analysis of genes in common associated with FEC1 and FEC2 traits in Australian Merino sheep. Categories included biological process (BP), cellular component (CC), molecular function (MF). On the x-axis the -log₁₀ of the adjusted p-value (<0.05) and on the y-axis the GO term.

Table 1. Descriptive statistics for fecal egg counts (LogFEC1 and LogFEC2).

Trait	N	Mean	SD^a	Min^b	Max^c
LogFEC1	26,638	6.63	1.14	4.60	10.51
LogFEC2	18,971	6.64	1.13	4.60	10.87

^a SD: standard deviation; ^b Min: minimum; ^c Max: maximum.

Table 2. Additive genetic variances, residuals, and heritability for FEC1 and FEC2 in Australian Merino sheep.

Trait	Method	$σ_{a}^{2}$	$σ_{e}^{2}$	h²
FEC1	BLUP	0.15 ± 0.01	0.66 ± 0.01	0.19 ± 0.01
	SSGBLUP	0.16 ± 0.01	0.66 ± 0.01	0.20 ± 0.01
FEC2	BLUP	0.15 ± 0.01	0.66 ± 0.01	0.19 ± 0.02
	SSGBLUP	0.16 ± 0.01	0.66 ± 0.01	0.19 ± 0.02

σ_{a}^{2}

= additive genetic variance;

σ_{e}^{2}

=residual variance;

h^{2}

= heritability.

Table 3. Chromosome, location, proportion of genetic variance, and candidate genes within the top 10 windows associated with the FEC 1 and FEC 2 traits in Australian Merino sheep.

Trait	Chr ^a	gVar(%) ^b	Windows Bounds (pb)	Candidate Genes
FEC1	6	0.83	35,429,074 - 35,554,875	GPRIN3
	6	0.73	37,072,532 - 37,134,706	LAP3, MED28, FAM184B
	11	0.61	26,119,722 - 26,407,016	MINK1, PLD2, ZMYND15, CXCL16, MED11, ARRB2, PELP1, ALOX15
	11	0.57	27,002,322 - 27,807,356	DNAH2, LOC101121185, SLC25A35, ODF4, NDEL1, MYH10
	25	0.49	21,324,145 - 22,654,093	CTNNA3
	11	0.36	25,298,833 - 25,633,178	WSCD1
	6	0.30	36,066,911 - 36,197,551	HERC3, PYURF, PIGY, HERC5
	21	0.30	42,654,067 - 42,734,621	SYVN1, SPDYC, CAPN1
	2	0.29	112,528,931 - 113,355,547	HERC2, NIPA1, NIPA2, CYFIP1, TUBGCP5, PTPN18, AMER3
	11	0.28	27,811,348 - 28,475,634	CCDC42, PIK3R5, STX8, USP43
FEC2	11	1.94	26,119,722 - 26,407,016	MINK1, PLD2, ZMYND15, CXCL16, MED11, ARRB2, PELP1, ALOX15
	5	1.13	93,463,383 - 93,482,724	CAST
	11	0.90	27,002,322 - 27,807,356	DNAH2, LOC101121185, SLC25A35, ODF4, NDEL1, MYH10
	11	0.77	27,811,348 - 28,475,634	CCDC42, PIK3R5, STX8, USP43
	11	0.76	25,298,833 - 25,633,178	WSCD1
	5	0.68	93,437,720 - 93,463,383	CAST
	6	0.64	37,072,532 - 37,134,706	LAP3, MED28, FAM184B
	6	0.53	35,429,074 - 35,554,875	GPRIN3
	11	0.52	25,744,499 - 26,098,771	C1QBP, RPAIN, NUP88, RABEP1, ZFP3, KIF1C, INCA1, CAMTA2, SPAG7
	6	0.51	37,483,582 - 37,665,116	-

^aChr = chromosome number. ^bgVar(%) = proportion of genetic variance represented by each region comprising 20 SNPs.

Table 4. Positional candidate genes in common for FEC1 and FEC2.

Chr	Pos (bp)	Candidate genes
6	35,511,497 - 37,257,065	FAM184B, MED28, LAP3, GPRIN3
11	25,298,833 - 28,475,634	USP43, STX8, PIK3R5, CCDC42, MYH10, NDEL1, ODF4, SLC25A35, LOC101121185, DNAH2, ALOX15, PELP1, ARRB2, CXCL16, MED11, ZMYND15, PLD2, MINK1, CAMTA2, SPAG7, INCA1, KIF1C, ZFP3, RABEP1, NUP88, C1QBP, RPAIN, WSCD1
25	22,143,113 - 22,654,093	CTNNA3

Table 5. Gene ontology (GO) terms such as biological processes, molecular functions, and KEGG pathways of candidate genes associated to FEC 1 and FEC 2 traits in Australian Merino sheep.

Category	Term	Genes	p-value	-log₁₀(p-value)
GOTERM_BP	GO:2001303~lipoxin A4 biosynthetic process	ALOX15, LOC101121185	0.005	2.289
GOTERM_BP	GO:0001764~neuron migration	NDEL1, ENSOARG00000002682, MYH10	0.007	2.140
GOTERM_BP	GO:0019372~lipoxygenase pathway	ALOX15, LOC101121185	0.010	1.989
GOTERM_BP	GO:0043651~linoleic acid metabolic process	ALOX15, LOC101121185	0.010	1.989
GOTERM_BP	GO:0007097~nuclear migration	ENSOARG00000002682, MYH10	0.014	1.865
GOTERM_BP	GO:0032695~negative regulation of interleukin-12 production	ARRB2, C1QBP	0.020	1.690
GOTERM_BP	GO:0019369~arachidonic acid metabolic process	ALOX15, LOC101121185	0.020	1.690
GOTERM_BP	GO:0030334~regulation of cell migration	MINK1, ENSOARG00000002682	0.075	1.128
GOTERM_BP	GO:0007286~spermatid development	CCDC42, ZMYND15	0.082	1.084
GOTERM_BP	GO:0007018~microtubule-based movement	DNAH2, KIF1C	0.093	1.029
GOTERM_CC	GO:0070847~core mediator complex	MED11, MED28	0.020	1.690
GOTERM_CC	GO:0005737~cytoplasm	NDEL1, LOC101121185, CTNNA3, ENSOARG00000002682, PELP1, INCA1, ZMYND15, ENSOARG00000007066	0.037	1.427
GOTERM_CC	GO:0030139~endocytic vesicle	RABEP1, ARRB2	0.049	1.313
GOTERM_CC	GO:0005819~spindle	NDEL1, MYH10	0.081	1.092
GOTERM_CC	GO:0005829~cytosol	ALOX15, LOC101121185, PIK3R5, STX8, LAP3, ARRB2, C1QBP, MYH10	0.082	1.086
GOTERM_MF	GO:0050473~arachidonate 15-lipoxygenase activity	ALOX15, LOC101121185	0.005	2.301
GOTERM_MF	GO:0016165~linoleate 13S-lipoxygenase activity	ALOX15, LOC101121185	0.007	2.176
KEGG_PATHWAY	oas04062:Chemokine signaling pathway	PIK3R5, ARRB2, CXCL16	0.033	1.486
KEGG_PATHWAY	oas04814:Motor proteins	DNAH2, MYH10, KIF1C	0.033	1.478
KEGG_PATHWAY	oas04144:Endocytosis	RABEP1, PLD2, ARRB2	0.055	1.260

Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

© 2024 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.

Submitted:

03 June 2024

Posted:

04 June 2024

You are already at the latest version

Alerts

Abstract

Keywords:

Subject: Biology and Life Sciences - Animal Science, Veterinary Science and Zoology

1. Introduction