A within- and across-country assessment of the genomic diversity and autozygosity of South African and eSwatini Nguni cattle.

In southern Africa, the Nguni cattle breed is classified as an indigenous and transboundary animal genetic resource (AnGR) that manifests unique adaptation abilities across distinct agroecological zones. The genetic integrity of various ecotypes is under potential threat due to both indiscriminate crossbreeding and uncontrolled inbreeding. The aim of this study was to assess the genetic diversity and autozygosity that exist both across countries (ES: eSwatini; SA: South Africa) and within countries (SA), between purebred stud animals (SA-S) and research herds (SA-R). Subsets of 96 ES, 96 SA-S, and 96 SA-R genotyped for 40,930 common SNPs were used to study genome-wide profiles of runs of homozygosity (ROH) and heterozygosity (ROHet) as well as inbreeding levels and population structure. The highest percentage (39.8%) of the 2168 ROH segments was 4-8 Mbp in length, whereas 65% of the 935 ROHet segments fell within the 0.5-1 Mbp length category. Inbreeding coefficients indicated positive but low inbreeding (FROH>1Mbp range: 0.025 for SA-S to 0.029 for SA-R). Principal component (PCA) and population structure analyses illustrated genome-level distinctness of (1) the Nguni from global indicine (Boran) and taurine (Hereford) breeds (K = 3), (2) the SA Nguni populations from the ES Nguni population (K = 4), and (3) different Nguni ecotypes within countries (K = 8). Furthermore, greater admixture was observed for the SA-R population compared to purebred SA-S population (shared ancestry = 0.631 ± 0.353 compared to 0.741 ± 0.123), and fewer genomics-defined ES ecotypes were observed than phenotypically (pre)defined. Overall, the results illustrated that genetic uniqueness within the sampled Nguni cattle resulted from both geographic isolation and exposure to different breeding strategies (and, selection pressures). A further loss of genetic variability should be monitored to prevent the endangerment of unique and beneficial ecotypes.

Population structure of indigenous southern African goats based on the Illumina Goat50K SNP panel.

In this study, the genetic structure of indigenous Tswana and Swazi goats using the Illumina Goat50K SNP array was investigated. Two South African commercial goat breeds were included to investigate admixture with the indigenous populations in southern Africa. A total of 144 DNA samples including Boer goats (n = 24), Kalahari Red (n = 24), Swazi (n = 48), and Tswana goats (n = 48) were genotyped. Statistical analysis was performed using PLINK version 1.07. Genetic diversity, measured as expected heterozygosity, was estimated at 0.390, 0.398, 0.413, and 0.387 for Boer, Kalahari Red, Tswana, and Swazi goats, respectively. The individual inbreeding coefficient varied from 0.019 ± 0.05 to 0.011 ± 0.06 for the Tswana and Swazi goats, respectively. The Principal component analysis clustered the populations according to geographical origin and breed type. Linkage disequilibrium (LD) for shorter intervals (0-10 kb) ranged from 0.44 to 0.56 and commercial breeds had higher values. Effective population sizes decreased with generations and at the 13th generation ranged between 87 for Boer to 266 for Tswana goats. The Tswana population exhibited the highest level of genetic variation and effective population size, which holds potential for improved production in marginal regions. A national strategy is required to maintain genetic diversity in communal goat production systems through well-structured breeding and conservation programs.