Landscape genomics and pathway analysis to understand genetic adaptation of South African indigenous goat populations.

In Africa, extensively raised livestock populations in most smallholder farming communities are exposed to harsh and heterogeneous climatic conditions and disease pathogens that they adapt to in order to survive. Majority of these livestock species, including goats, are of non-descript and uncharacterized breeds and their response to natural selection presented by heterogeneous environments is still unresolved. This study investigated genetic diversity and its association with environmental and geographic conditions in 194 South African indigenous goats from different geographic locations genotyped on the Illumina goat SNP50K panel. Population structure analysis revealed a homogeneous genetic cluster of the Tankwa goats, restricted to the Northern Cape province. Overall, the Boer, Kalahari Red, and Savanna showed a wide geographic spread of shared genetic components, whereas the village ecotypes revealed a longitudinal distribution. The relative importance of environmental factors on genetic variation of goat populations was assessed using redundancy analysis (RDA). Climatic and geographic variables explained 22% of the total variation while climatic variables alone accounted for 17% of the diversity. Geographic variables solitarily explained 1% of the total variation. The first axis (Model I) of the RDA analysis revealed 329 outlier SNPs. Landscape genomic approaches of spatial analysis method (SAM) identified a total of 843 (1.75%) SNPs, while latent factor mixed models (LFMM) identified 714 (1.48%) SNPs significantly associated with environmental variables. Significant markers were within genes involved in biological functions potentially important for environmental adaptation. Overall, the study suggested environmental factors to have some effect in shaping the genetic variation of South African indigenous goat populations. Loci observed to be significant and under selection may be responsible for the adaption of the goat populations to local production systems.

Population genomic structure and linkage disequilibrium analysis of South African goat breeds using genome-wide SNP data.

The sustainability of goat farming in marginal areas of southern Africa depends on local breeds that are adapted to specific agro-ecological conditions. Unimproved non-descript goats are the main genetic resources used for the development of commercial meat-type breeds of South Africa. Little is known about genetic diversity and the genetics of adaptation of these indigenous goat populations. This study investigated the genetic diversity, population structure and breed relations, linkage disequilibrium, effective population size and persistence of gametic phase in goat populations of South Africa. Three locally developed meat-type breeds of the Boer (n = 33), Savanna (n = 31), Kalahari Red (n = 40), a feral breed of Tankwa (n = 25) and unimproved non-descript village ecotypes (n = 110) from four goat-producing provinces of the Eastern Cape, KwaZulu-Natal, Limpopo and North West were assessed using the Illumina Goat 50K SNP Bead Chip assay. The proportion of SNPs with minor allele frequencies >0.05 ranged from 84.22% in the Tankwa to 97.58% in the Xhosa ecotype, with a mean of 0.32 ± 0.13 across populations. Principal components analysis, admixture and pairwise FST identified Tankwa as a genetically distinct population and supported clustering of the populations according to their historical origins. Genome-wide FST identified 101 markers potentially under positive selection in the Tankwa. Average linkage disequilibrium was highest in the Tankwa (r(2)  = 0.25 ± 0.26) and lowest in the village ecotypes (r(2) range = 0.09 ± 0.12 to 0.11 ± 0.14). We observed an effective population size of <150 for all populations 13 generations ago. The estimated correlations for all breed pairs were lower than 0.80 at marker distances >100 kb with the exception of those in Savanna and Tswana populations. This study highlights the high level of genetic diversity in South African indigenous goats as well as the utility of the genome-wide SNP marker panels in genetic studies of these populations.

Targeted high-throughput growth hormone 1 gene sequencing reveals high within-breed genetic diversity in South African goats.

This study assessed the genetic diversity in the growth hormone 1 gene (GH1) within and between South African goat breeds. Polymerase chain reaction-targeted gene amplification together with Illumina MiSeq next-generation sequencing (NGS) was used to generate the full length (2.54 kb) of the growth hormone 1 gene and screen for SNPs in the South African Boer (SAB) (n = 17), Tankwa (n = 15) and South African village (n = 35) goat populations. A range of 27-58 SNPs per population were observed. Mutations resulting in amino acid changes were observed at exons 2 and 5. Higher within-breed diversity of 97.37% was observed within the population category consisting of SA village ecotypes and the Tankwa goats. Highest pairwise FST values ranging from 0.148 to 0.356 were observed between the SAB and both the South African village and Tankwa feral goat populations. Phylogenetic analysis indicated nine genetic clusters, which reflected close relationships between the South African populations and the other international breeds with the exception of the Italian Sarda breeds. Results imply greater potential for within-population selection programs, particularly with SA village goats.