AUTHOR=Chokoe T. C. , Hadebe K. , Muchadeyi F. C. , Nephawe K. A. , Dzomba E. F. , Mphahlele T. D. , Matelele T. C. , Mtileni B. J. 

TITLE=Conservation status and historical relatedness of South African communal indigenous goat populations using a genome-wide single-nucleotide polymorphism marker

JOURNAL=Frontiers in Genetics

VOLUME=13

YEAR=2022

URL=https://www.frontiersin.org/journals/genetics/articles/10.3389/fgene.2022.909472

DOI=10.3389/fgene.2022.909472

ISSN=1664-8021

ABSTRACT=<p>Indigenous goats form the majority of populations in smallholder, low input, low output production systems and are considered an important genetic resource due to their adaptability to different production environments and support of communal farming. Effective population size (<italic>N</italic><sub><italic>e</italic></sub>), inbreeding levels, and the runs of homozygosity (ROHs) are effective tools for exploring the genetic diversity and understanding the demographic history in efforts to support breeding strategies to use and conserve genetic resources. Across populations, the current <italic>N</italic><sub><italic>e</italic></sub> of Gauteng was the lowest at 371 animals, while the historical <italic>N</italic><sub><italic>e</italic></sub> across populations suggests that the ancestor <italic>N</italic><sub><italic>e</italic></sub> has decreased by 53.86%, 44.58%, 42.16%, and 41.16% in Free State (FS), North West (NW), Limpopo (LP), and Gauteng (GP), respectively, over the last 971 generations. Genomic inbreeding levels related to ancient kinship (<italic>F</italic><sub><italic>ROH</italic></sub> &gt; 5 Mb) were highest in FS (0.08 ± 0.09) and lowest in the Eastern Cape (EC) (0.02 ± 0.02). A total of 871 ROH island regions which include important environmental adaptation and hermo-tolerance genes such as <italic>IL10RB, IL23A, FGF9, IGF1, EGR1, MTOR,</italic> and <italic>MAPK3</italic> were identified (occurring in over 20% of the samples) in FS (<italic>n</italic> = 37), GP (<italic>n</italic> = 42), and NW (<italic>n</italic> = 2) populations only. The mean length of ROH across populations was 7.76 Mb and ranged from 1.61 Mb in KwaZulu-Natal (KZN) to 98.05 Mb (GP and NW). The distribution of ROH according to their size showed that the majority (<italic>n</italic> = 1949) of the detected ROH were &gt; 5 Mb in length compared to the other categories. Assuming two hypothetical ancestral populations, the populations from KZN and LP are revealed, supporting PC 1. The genomes of KZN and LP share a common origin but have substantial admixture from the EC and NW populations. The findings revealed that the occurrence of high <italic>N</italic><sub><italic>e</italic></sub> and autozygosity varied largely across breeds in communal indigenous goat populations at recent and ancient events when a genome-wide single-nucleotide polymorphism (SNP) marker was used. The use of Illumina goat SNP50K BeadChip shows that there was a migration route of communal indigenous goat populations from the northern part (LP) of South Africa to the eastern areas of the KZN that confirmed their historical relatedness and coincides with the migration periods of the Bantu nation.</p>