Abstract Junipers (Juniperus spp.) are woody invasive plants that have caused encroachment problems in the U.S. western rangeland, decreasing forage productivity and biodiversity. A potential solution to this issue is the use of goats in target grazing programs. Selecting goats for juniper consumption is of great interest to rangeland farmers. Therefore, our main objectives were to 1) estimate genetic parameters, including heritability and repeatability, for juniper consumption in Angora (ANG) and Boer x Spanish (BxS) composite population of goats grazing on West Texas rangelands (USA); and, 2) perform a genome-wide association study (GWAS) to identify candidate genes associated with juniper consumption. Fecal samples were collected multiple times between 2008 and 2017 and analyzed using near-infrared reflectance spectroscopy (NIRS) to estimate the percentage of juniper in the diet of goats, using calibration equations previously developed in controlled studies. We analyzed a dataset comprising 2,327 records from 459 ANG and 8,714 records from 1,996 BxS. The pedigree file included 1,601 and 2,426 ANG and BxS goats, respectively. Genomic information for 50,957 single nucleotide polymorphisms (SNP) from 93 ANG and 238 BxS goats was included in the analyses. Sampling date and sex were fitted as categorical fixed effects and age as a linear covariate in the statistical models. Animal direct additive genetic and permanent environment were considered as random effects. Genetic parameters were calculated based on the Restricted Maximum Likelihood method under the single-step genomic BLUP procedure implemented in the BLUPf90+ family of programs. We used de-regressed estimated breeding values and whole-genome sequence data for the GWAS analyses with 8 million SNP for the ANG and 5.3 million SNP for the BxS. The software GCTA was employed for these analyses. Juniper percentage in their diets varied across the years and ranged from 2.14% (± 7.51%) to 14.6% (± 11.0%) for ANG and from 10.5% (± 14.7) to 33.8% (± 18.1) for the BxS population. The heritability estimate for juniper consumption was 0.19 (± 0.038) for ANG and 0.28 (± 0.022) for BxS, indicating that juniper consumption can be increased through selection. The repeatability values were 0.21 for ANG and 0.37 for BxS. The GWAS revealed 1,257 significant SNP located within or close to 220 genes in ANG, and 878 SNP related to 159 genes in BxS, suggesting their association with juniper consumption. These findings suggest that juniper consumption is under genetic control and has a polygenic inheritance pattern influenced by multiple genes of small effect. The heritability and repeatability estimates indicate that genetic improvement of juniper consumption is feasible in both goat populations.
Read full abstract