Search

Domesticated sheep have adapted to contrasting and extreme environments and continue to play important roles in local community-based economies throughout Africa. Here we review the Neolithic migrations of thin-tailed sheep and the later introductions of fat-tailed sheep into eastern Africa. According to contemporary pictorial evidence, the latter occurred in Egypt not before the Ptolemaic period (305-25 BCE). We further describe the more recent history of sheep in Egypt, the Maghreb, west and central Africa, central-east Africa, and southern Africa. We also present a comprehensive molecular survey based on the analysis of 50 K SNP genotypes for 59 African breeds contributed by several laboratories. We propose that gene flow and import of fat-tailed sheep have partially overwritten the diversity profile created by the initial migration. We found a genetic contrast between sheep north and south of the Sahara and a west-east contrast of thin- and fat-tailed sheep. There is no close relationship between African and central and east Asian fat-tailed breeds, whereas we observe within Africa only a modest effect of tail types on breed relationships., (© 2024 The Author(s). Animal Genetics published by John Wiley & Sons Ltd on behalf of Stichting International Foundation for Animal Genetics.)

South Africa boasts a diverse range of pig populations, encompassing intensively raised commercial breeds, as well as indigenous and village pigs reared under low-input production systems. The aim of this study was to investigate how natural and artificial selection have shaped the genomic landscape of South African pig populations sampled from different genetic backgrounds and production systems. For this purpose, the integrated haplotype score (iHS), as well as cross population extended haplotype homozygosity (XP-EHH) and Lewontin and Krakauer's extension of the Fst statistic based on haplotype information (HapFLK) were utilised. Our results revealed several population-specific signatures of selection associated with the different production systems. The importance of natural selection in village populations was highlighted, as the majority of genomic regions under selection were identified in these populations. Regions under natural and artificial selection causing the distinct genetic footprints of these populations also allow for the identification of genes and pathways that may influence production and adaptation. In the context of intensively raised commercial pig breeds (Large White, Kolbroek, and Windsnyer), the identified regions included quantitative loci (QTLs) associated with economically important traits. For example, meat and carcass QTLs were prevalent in all the populations, showing the potential of village and indigenous populations' ability to be managed and improved for such traits. Results of this study therefore increase our understanding of the intricate interplay between selection pressures, genomic adaptations, and desirable traits within South African pig populations.

In this study, we investigated how IBD patterns shared between individuals of the same breed could be informative of its admixture level, with the underlying assumption that the most admixed breeds, i.e. the least genetically isolated, should have a much more fragmented genome. We considered 111 goat breeds (i.e. 2501 individuals) and 156 sheep breeds (i.e. 3304 individuals) from Europe, Africa and Asia, for which beadchip SNP genotypes had been performed. We inferred the breed's level of admixture from: (i) the proportion of the genome shared by breed's members (i.e. "genetic integrity level" assessed from ADMIXTURE software analyses), and (ii) the "AV index" (calculated from Reynolds' genetic distances), used as a proxy for the "genetic distinctiveness". In both goat and sheep datasets, the statistical analyses (comparison of means, Spearman correlations, LM and GAM models) revealed that the most genetically isolated breeds, also showed IBD profiles made up of more shared IBD segments, which were also longer. These results pave the way for further research that could lead to the development of admixture indicators, based on the characterization of intra-breed shared IBD segments, particularly effective as they would be independent of the knowledge of the whole genetic landscape in which the breeds evolve. Finally, by highlighting the fragmentation experienced by the genomes subjected to crossbreeding carried out over the last few generations, the study reminds us of the need to preserve local breeds and the integrity of their adaptive architectures that have been shaped over the centuries., (© 2023. The Author(s).)

The African Goat Improvement Network (AGIN) is a collaborative group of scientists focused on genetic improvement of goats in small holder communities across the African continent. The group emerged from a series of workshops focused on enhancing goat productivity and sustainability. Discussions began in 2011 at the inaugural workshop held in Nairobi, Kenya. The goals of this diverse group were to: improve indigenous goat production in Africa; characterize existing goat populations and to facilitate germplasm preservation where appropriate; and to genomic approaches to better understand adaptation. The long-term goal was to develop cost-effective strategies to apply genomics to improve productivity of small holder farmers without sacrificing adaptation. Genome-wide information on genetic variation enabled genetic diversity studies, facilitated improved germplasm preservation decisions, and provided information necessary to initiate large scale genetic improvement programs. These improvements were partially implemented through a series of community-based breeding programs that engaged and empowered local small farmers, especially women, to promote sustainability of the production system. As with many international collaborative efforts, the AGIN work serves as a platform for human capacity development. This paper chronicles the evolution of the collaborative approach leading to the current AGIN organization and describes how it builds capacity for sustained research and development long after the initial program funds are gone. It is unique in its effectiveness for simultaneous, multi-level capacity building for researchers, students, farmers and communities, and local and regional government officials. The positive impact of AGIN capacity building has been felt by participants from developing, as well as developed country partners., Competing Interests: Author TSS was employed by Acceligen Inc. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2023 Van Tassell, Rosen, Woodward-Greene, Silverstein, Huson, Sölkner, Boettcher, Rothschild, Mészáros, Nakimbugwe, Gondwe, Muchadeyi, Nandolo, Mulindwa, Banda, Kaumbata, Getachew, Haile, Soudre, Ouédraogo, Rischkowsky, Mwai, Dzomba, Nash, Abegaz, Masiga, Wurzinger, Sayre, Stella, Tosser-Klopp and Sonstegard.)

Background: To enhance and extend the knowledge about the global historical and phylogenetic relationships between Merino and Merino-derived breeds, 19 populations were genotyped with the OvineSNP50 BeadChip specifically for this study, while an additional 23 populations from the publicly available genotypes were retrieved. Three complementary statistical tests, Rsb (extended haplotype homozygosity between-populations), XP-EHH (cross-population extended haplotype homozygosity), and runs of homozygosity (ROH) islands were applied to identify genomic variants with potential impact on the adaptability of Merino genetic type in two contrasting climate zones., Results: The results indicate that a large part of the Merino's genetic relatedness and admixture patterns are explained by their genetic background and/or geographic origin, followed by local admixture. Multi-dimensional scaling, Neighbor-Net, Admixture, and TREEMIX analyses consistently provided evidence of the role of Australian, Rambouillet and German strains in the extensive gene introgression into the other Merino and Merino-derived breeds. The close relationship between Iberian Merinos and other South-western European breeds is consistent with the Iberian origin of the Merino genetic type, with traces from previous contributions of other Mediterranean stocks. Using Rsb and XP-EHH approaches, signatures of selection were detected spanning four genomic regions located on Ovis aries chromosomes (OAR) 1, 6 and 16, whereas two genomic regions on OAR6, that partially overlapped with the previous ones, were highlighted by ROH islands. Overall, the three approaches identified 106 candidate genes putatively under selection. Among them, genes related to immune response were identified via the gene interaction network. In addition, several candidate genes were found, such as LEKR1, LCORL, GHR, RBPJ, BMPR1B, PPARGC1A, and PRKAA1, related to morphological, growth and reproductive traits, adaptive thermogenesis, and hypoxia responses., Conclusions: To the best of our knowledge, this is the first comprehensive dataset that includes most of the Merino and Merino-derived sheep breeds raised in different regions of the world. The results provide an in-depth picture of the genetic makeup of the current Merino and Merino-derived breeds, highlighting the possible selection pressures associated with the combined effect of anthropic and environmental factors. The study underlines the importance of Merino genetic types as invaluable resources of possible adaptive diversity in the context of the occurring climate changes., (© 2023. The Author(s).)

Merino sheep are a breed of choice across the world, popularly kept for their wool and mutton value. They are often reared as a pure breed or used in crossbreeding and are a common component in synthetic breed development. This study evaluated genetic diversity, population structure, and breed divergence in 279 animals of Merino and Merino-based sheep breeds in South Africa using the Illumina Ovine SNP 50K BeadChip. The sheep breeds analysed included the three Merino-derived breeds of Dohne Merino ( n = 50); Meatmaster ( n = 47); and Afrino ( n = 52) and five presumed ancestral populations of Merinos (Merino ( n = 46); South African Merino ( n = 10); and South African Mutton Merino ( n = 8)); and the non-Merino founding breeds of Damara ( n = 20); Ronderib Afrikaner ( n = 17); and Nguni ( n = 29). Highest genetic diversity values were observed in the Dohne Merino (DM), with H o = 0.39 ± 0.01, followed by the Meatmaster and South African Merino (SAM), with H o = 0.37 ± 0.03. The level of inbreeding ranged from 0.0 ± 0.02 (DM) to 0.27 ± 0.05 (Nguni). Analysis of molecular variance (AMOVA) showed high within-population variance (>80%) across all population categories. The first principal component (PC1) separated the Merino, South African Mutton Merino (SAMM), DM, and Afrino (AFR) from the Meatmaster, Damara, Nguni, and Ronderib Afrikaner (RDA). PC2 aligned each Merino-derived breed with its presumed ancestors and separated the SAMM from the Merino and SAM. The iHS analysis yielded selection sweeps across the AFR (12 sweeps), Meatmaster (four sweeps), and DM (29 sweeps). Hair/wool trait genes such as FGF12 ; metabolic genes of ICA1 , NXPH1 , and GPR171 ; and immune response genes of IL22 , IL26 , IFNAR1 , and IL10RB were reported. Other genes include HMGA , which was observed as selection signatures in other populations; WNT5A , important in the development of the skeleton and mammary glands; ANTXR2 , associated with adaptation to variation in climatic conditions; and BMP2 , which has been reported as strongly selected in both fat-tailed and thin-tailed sheep. The DM vs. SAMM shared all six sweep regions on chromosomes 1, 10, and 11 with AFR vs. SAMM. Genes such as FGF12 on OAR 1:191.3-194.7 Mb and MAP2K4 on OAR 11:28.6-31.3 Mb were observed. The selection sweep on chromosome 10 region 28.6-30.3 Mb harbouring the RXFP2 for polledness was shared between the DM vs. Merino, the Meatmaster vs. Merino, and the Meatmaster vs. Nguni. The DM vs. Merino and the Meatmaster vs. Merino also shared an Rsb-based selection sweep on chromosome 1 region 268.5-269.9 Mb associated with the Calpain gene, CAPN7. The study demonstrated some genetic similarities between the Merino and Merino-derived breeds emanating from common founding populations and some divergence driven by breed-specific selection goals. Overall, information regarding the evolution of these composite breeds from their founding population will guide future breed improvement programs and management and conservation efforts., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2023 Dzomba, Van Der Nest, Mthembu, Soma, Snyman, Chimonyo and Muchadeyi.)

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 (F ROH>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., (© 2022. The Author(s), under exclusive licence to Springer Nature B.V.)

Growth and carcass quality are economically important traits in goat production. This study investigated differentially expressed genes from the caprine pituitary gland transcriptome of South African indigenous goat breeds of varying growth performances and carcass quality parameters. Tissues were harvested from the pituitary gland of three South African Boer goats and three village ecotype goats all raised under similar conditions simulating intensive commercial production systems. Three additional tissues were harvested from village ecotype goats that were raised extensively on village farms. Between breed differences were investigated by comparing differential gene expression among three South African Boer and three village goats that were both raised under intensive commercial production system at a research farm. Within-breed differences were investigated by comparing differential gene expression among three village goats raised under extensive conditions (on-farm in Pella, S.A. village farming community) and three village goats raised under intensive commercial production system (at ARC research farm in Pretoria, South Africa. Total RNA was isolated from the pituitary gland of 36-week-old animals ( n = 9 ) and sequenced individually in triplicates. An average of 28,298,512 trimmed, and quality-controlled reads/animal were mapped to the goat genome ( Capra_hircus.ARS1.94 ) using HiSat2 software. Transcript assembly and quantification yielded 104 differentially expressed genes for village goats raised under extensive system and 62 for village goats raised under the intensive production system at the false discovery rate (FRD) of ≤0.05 and a fold change of ≥2. Growth-related genes such as POU3F4 and TSHZ1 were highly expressed within breeds raised under both production systems. Conversely, growth-related genes such as FGFR2 and SMPX genes were highly expressed between breeds raised under similar production systems. Ballgown analysis revealed a high expression of GH1 and IGF1 in the intensively raised compared to extensively raised goats. Both genes were also highly expressed in the village goats when compared to the Boer. The differential gene expression data provided insights into genes and molecular mechanisms associated with growth and growth development in goats., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Ncube, Dzomba, Rosen, Schroeder, Van Tassell and Muchadeyi.)

The Nguni cattle of South Africa are a Sanga breed, characterized by many eco-types and research populations that have been established in an effort to conserve the diversity within the breed. The aim of this study was to investigate the overall genetic diversity as well as similarities and differences within and between two conservation herds of the South African Nguni Cattle. Mean LD ( r 2 ) estimates were 0.413 ± 0.219 for Bartlow Combine and 0.402 ± 0.209 for Kokstad. Genome-wide average LD ( r 2 ) decreased with increasing genetic marker distance for both populations from an average of 0.76 ± 0.28 and 0.77 ± 0.27 at 0-1 kb bin to 0.31 ± 0.13 and 0.32 ± 0.13 at 900-1000 kb bin in Bartlow Combine and Kokstad populations, respectively. Variation in LD levels across autosomes was observed in both populations. The results showed higher levels of LD than previously reported in Nguni field populations and other South African breeds, especially at shorter marker distances of less than 20 kb. A total number of 77,305 and 66,237 haplotype blocks covering a total of 1570.09 Mb (61.99% genome coverage) and 1367.42 Mb (53.96% genome coverage) were detected in Bartlow Combine and Kokstad populations, respectively. A total of 18,449 haploblocks were shared between the two populations while 58,856 and 47,788 haploblocks were unique to Bartlow Combine and Kokstad populations, respectively. Effective population size (N e ) results demonstrated a rapid decrease in N e across generations for both Bartlow Combine and Kokstad conservation herds. Two complementary methods, integrated haplotype score (iHS) and Extend Haplotype Homozygosity Test (XP-EHH), were implemented in this study to detect the selection signatures in the two herds. A total of 553 and 166 selected regions were identified in Bartlow Combine and Kokstad populations, respectively. DAVID and GO terms analysis of the regions under selection reported genes/QTLs associated with fertility, carcass weight, coat colour, immune response, and eye area pigmentation. Some genes, such as HCAR1, GNAI1, PIK3R3, WNT3, RAB5A, BOLA-N (Class IB MHC Antigen QA-2-Related), BOLA (Class IB MHC Antigen QA-2-Related), and Rab-8B, etc., were found in regions under selection in this study. Overall, the study implied reduced genetic diversity in the two herds calling for corrective measures to maintain the diversity of the South African Nguni cattle. This study presented a comprehensive analysis of the genomic architecture of South African Nguni cattle populations, providing essential genetic information of utility in the management of conservation flocks.

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 ( N e ), 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 N e of Gauteng was the lowest at 371 animals, while the historical N e across populations suggests that the ancestor N e 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 ( F ROH > 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 IL10RB, IL23A, FGF9, IGF1, EGR1, MTOR, and MAPK3 were identified (occurring in over 20% of the samples) in FS ( n = 37), GP ( n = 42), and NW ( n = 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 ( n = 1949) of the detected ROH were > 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 N e 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., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Chokoe, Hadebe, Muchadeyi, Nephawe, Dzomba, Mphahlele, Matelele and Mtileni.)

Nguni cattle are a Sanga type breed with mixed B. taurus and B. indicus ancestry and proven resistance to ticks, diseases and other harsh conditions of the African geographical landscape. The multi-coloured Nguni coats have found a niche market in the leather industry leading to breeding objectives towards the promotion of such diversity. However, there is limited studies on the genomic architecture underlying the coat colour and patterns hampering any potential breeding and improvement of such trait. This study investigated the genetics of base coat colour, colour-sidedness and the white forehead stripe in Nguni cattle using coat colour phenotyped Nguni cattle and Illumina Bovine HD (770K) genotypes. Base coat colour phenotypes were categorised into eumelanin ( n = 45) and pheomelanin ( n = 19). Animals were categorised into either colour-sided ( n = 46) or non-colour-sided ( n = 94) and similarly into presence ( n = 15) or absence ( n = 67) of white forehead stripe. Genome-wide association tests were conducted using 622,103 quality controlled SNPs and the Efficient Mixed Model Association eXpedited method (EMMAX) implemented in Golden Helix SNP Variation Suite. The genome-wide association studies for base coat colour (eumelanin vs. pheomelanin) resulted into four indicative SNPs on BTA18 and a well-known gene, MC1R, was observed within 1 MB from the indicative SNPs (p < 0.00001) and found to play a role in the melanogenesis (core pathway for melanin production) and the MAPK signalling pathway. GWAS for colour-sidedness resulted in four indicative SNPs, none of which were in close proximity to the KIT candidate gene known for colour-sidedness. GWAS for the white forehead stripe resulted in 17 indicative SNPs on BTA6. Four genes MAPK10, EFNA5, PPP2R3C and PAK1 were found to be associated with the white forehead stripe and were part of the MAPK, adrenergic and Wnt signalling pathways that are synergistically associated with the synthesis of melanin. Overall, our results prove prior knowledge of the role of MC1R in base coat colours in cattle and suggested a different genetic mechanism for forehead stripe phenotypes in Nguni cattle., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Kunene, Muchadeyi, Hadebe, Mészáros, Sölkner, Dugmore and Dzomba.)

Processed meat is a target in meat adulteration for economic gain. This study demonstrates a molecular and bioinformatics diagnostic pipeline, utilizing the mitochondrial 16S ribosomal RNA (rRNA) gene, to determine processed meat product mislabeling through Next-Generation Sequencing. Nine pure meat samples were collected and artificially mixed at different ratios to verify the specificity and sensitivity of the pipeline. Processed meat products (n = 155), namely, minced meat, biltong, burger patties, and sausages, were collected across South Africa. Sequencing was performed using the Illumina MiSeq sequencing platform. Each sample had paired-end reads with a length of ±300 bp. Quality control and filtering was performed using BBDuk (version 37.90a). Each sample had an average of 134,000 reads aligned to the mitochondrial genomes using BBMap v37.90. All species in the artificial DNA mixtures were detected. Processed meat samples had reads that mapped to the Bos (90% and above) genus, with traces of reads mapping to Sus and Ovis (2-5%) genus. Sausage samples showed the highest level of contamination with 46% of the samples having mixtures of beef, pork, or mutton in one sample. This method can be used to authenticate meat products, investigate, and manage any form of mislabeling.

Carcass quality includes a battery of essential economic meat traits that play a significant role in influencing farmer breed preferences. A preliminary study was undertaken to investigate the carcass quality and the associated genomic regions in a small nucleus of animals that are representative of South African goat genetic resources. Samples of the South African Boer ( n = 14), Northern Cape Speckled ( n = 14), Eastern Cape Xhosa Lob ear ( n = 12), Nguni/Mbuzi ( n = 13), and Village ( n = 20) were genotyped using the Illumina goat SNP50K and were phenotyped for carcass quality traits. SA Boer goats had heavier warm and cold carcass weights (17.2 ± 2.3 kg and 16.3 ± 2.3 kg). Pella village goats raised under an intensive system had significantly ( p < 0.05) heavier warm and cold carcass weights (9.9 ± 1.1 kg and 9.2 ± 1.2 kg) compared to the village goats that are raised extensively (9.1 ± 2.0 kg and 8.4 ± 1.9). A total of 40 SNPs located on chromosomes 6, 10, 12, 13, 19, and 21 were significantly associated with carcass traits at (-log10 [ p < 0.05]). Candidate genes that were associated with carcass characteristics ( GADD45G , IGF2R , GAS1 , VAV3 , CAPN8 , CAPN7 , CAPN2 , GHSR , COLQ , MRAS , and POU1F1 ) were also observed. Results from this study will inform larger future studies that will ultimately find use in breed improvement programs.

A major obstacle in applying genomic selection (GS) to uniquely adapted local breeds in less-developed countries has been the cost of genotyping at high densities of single-nucleotide polymorphisms (SNP). Cost reduction can be achieved by imputing genotypes from lower to higher densities. Locally adapted breeds tend to be admixed and exhibit a high degree of genomic heterogeneity thus necessitating the optimization of SNP selection for downstream imputation. The aim of this study was to quantify the achievable imputation accuracy for a sample of 1,135 South African (SA) Drakensberger cattle using several custom-derived lower-density panels varying in both SNP density and how the SNP were selected. From a pool of 120,608 genotyped SNP, subsets of SNP were chosen (1) at random, (2) with even genomic dispersion, (3) by maximizing the mean minor allele frequency (MAF), (4) using a combined score of MAF and linkage disequilibrium (LD), (5) using a partitioning-around-medoids (PAM) algorithm, and finally (6) using a hierarchical LD-based clustering algorithm. Imputation accuracy to higher density improved as SNP density increased; animal-wise imputation accuracy defined as the within-animal correlation between the imputed and actual alleles ranged from 0.625 to 0.990 when 2,500 randomly selected SNP were chosen vs. a range of 0.918 to 0.999 when 50,000 randomly selected SNP were used. At a panel density of 10,000 SNP, the mean (standard deviation) animal-wise allele concordance rate was 0.976 (0.018) vs. 0.982 (0.014) when the worst (i.e., random) as opposed to the best (i.e., combination of MAF and LD) SNP selection strategy was employed. A difference of 0.071 units was observed between the mean correlation-based accuracy of imputed SNP categorized as low (0.01 < MAF ≤ 0.1) vs. high MAF (0.4 < MAF ≤ 0.5). Greater mean imputation accuracy was achieved for SNP located on autosomal extremes when these regions were populated with more SNP. The presented results suggested that genotype imputation can be a practical cost-saving strategy for indigenous breeds such as the SA Drakensberger. Based on the results, a genotyping panel consisting of ~10,000 SNP selected based on a combination of MAF and LD would suffice in achieving a <3% imputation error rate for a breed characterized by genomic admixture on the condition that these SNP are selected based on breed-specific selection criteria., (© The Author(s) 2021. Published by Oxford University Press on behalf of the American Society of Animal Science. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

Ungulates are a group of hoofed animals that have long interacted with humans as essential sources of food, labor, clothing, and transportation. These consist of domesticated, feral, and wild species raised in a wide range of habitats and biomes. Given the diverse and extreme environments inhabited by ungulates, unique adaptive traits are fundamental for fitness. The documentation of genes that underlie their genomic signatures of selection is crucial in this regard. The increasing availability of advanced sequencing technologies has seen the rapid growth of ungulate genomic resources, which offers an exceptional opportunity to understand their adaptive evolution. Here, we summarize the current knowledge on evolutionary genetic signatures underlying the adaptations of ungulates to different habitats.

A metagenomic approach was used to study the gut microbiome of Haemonchus contortus field strains and that of its predilection site, the abomasum of Dohne Merino sheep. The abomasum contents and H. contortus were collected from 10 naturally infected Dohne Merino sheep. The H. contortus specimens were classified and sexually differentiated using morphometric characters and was further confirmed through molecular identification. We investigated differences and similarities between the bacterial composition of the adult male and female H. contortus gut microbiomes, which were both dominated by bacteria from the Escherichia, Shigella, Vibrio and Halomonas genera. Major abundance variations were identified between the shared adult male and female H. contortus microbiomes. The results also revealed that Succiniclasticum, Rikenellaceae RC9 gut group and Candidatus Saccharimonas were the predominant genera in the Dohne Merino abomasum. This study provides insight into the highly diverse bacterial composition of the H. contortus gut microbiome and the Dohne Merino abomasum which needs to be studied further to explore the complex interactions of different gastrointestinal nematode microbiomes with the host.

In this study, we evaluated an admixed South African Simbra crossbred population, as well as the Brahman (Indicine) and Simmental (Taurine) ancestor populations to understand their genetic architecture and detect genomic regions showing signatures of selection. Animals were genotyped using the Illumina BovineLD v2 BeadChip (7K). Genomic structure analysis confirmed that the South African Simbra cattle have an admixed genome, composed of 5/8 Taurine and 3/8 Indicine, ensuring that the Simbra genome maintains favorable traits from both breeds. Genomic regions that have been targeted by selection were detected using the linkage disequilibrium-based methods iHS and Rsb. These analyses identified 10 candidate regions that are potentially under strong positive selection, containing genes implicated in cattle health and production (e.g., TRIM63, KCNA10, NCAM1, SMIM5, MIER3, and SLC24A4). These adaptive alleles likely contribute to the biological and cellular functions determining phenotype in the Simbra hybrid cattle breed. Our data suggested that these alleles were introgressed from the breed's original indicine and taurine ancestors. The Simbra breed thus possesses derived parental alleles that combine the superior traits of the founder Brahman and Simmental breeds. These regions and genes might represent good targets for ad-hoc physiological studies, selection of breeding material and eventually even gene editing, for improved traits in modern cattle breeds. This study represents an important step toward developing and improving strategies for selection and population breeding to ultimately contribute meaningfully to the beef production industry., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 van der Nest, Hlongwane, Hadebe, Chan, van der Merwe, De Vos, Greyling, Kooverjee, Soma, Dzomba, Bradfield and Muchadeyi.)

Background: Population history, production system and within-breed selection pressure impacts the genome architecture resulting in reduced genetic diversity and increased frequency of runs of homozygosity islands. This study tested the hypothesis that production systems geared towards specific traits of importance or natural or artificial selection pressures influenced the occurrence and distribution of runs of homozygosity (ROH) in the South African sheep population. The Illumina OvineSNP50 BeadChip was used to genotype 400 sheep belonging to 13 breeds from South Africa representing mutton, pelt and mutton and wool dual-purpose breeds, including indigenous non-descript breeds that are reared by smallholder farmers. To get more insight into the autozygosity and distribution of ROH islands of South African breeds relative to global populations, 623 genotypes of sheep from worldwide populations were included in the analysis. Runs of homozygosity were computed at cut-offs of 1-6 Mb, 6-12 Mb, 12-24 Mb, 24-48 Mb and > 48 Mb, using the R package detectRUNS. The Golden Helix SVS program was used to investigate the ROH islands., Results: A total of 121,399 ROH with mean number of ROH per animal per breed ranging from 800 (African White Dorper) to 15,097 (Australian Poll Dorset) were obtained. Analysis of the distribution of ROH according to their size showed that, for all breeds, the majority of the detected ROH were in the short (1-6 Mb) category (88.2%). Most animals had no ROH > 48 Mb. Of the South African breeds, the Nguni and the Blackhead Persian displayed high ROH based inbreeding (F ROH ) of 0.31 ± 0.05 and 0.31 ± 0.04, respectively. Highest incidence of common runs per SNP across breeds was observed on chromosome 10 with over 250 incidences of common ROHs. Mean proportion of SNPs per breed per ROH island ranged from 0.02 ± 0.15 (island ROH224 on chromosome 23) to 0.13 ± 0.29 (island ROH175 on chromosome 15). Seventeen (17) of the islands had SNPs observed in single populations (unique ROH islands). The MacArthur Merino (MCM) population had five unique ROH islands followed by Blackhead Persian and Nguni with three each whilst the South African Mutton Merino, SA Merino, White Vital Swakara, Karakul, Dorset Horn and Chinese Merino each had one unique ROH island. Genes within ROH islands were associated with predominantly metabolic and immune response traits and predomestic selection for traits such as presence or absence of horns., Conclusions: Overall, the frequency and patterns of distribution of ROH observed in this study corresponds to the breed history and implied selection pressures exposed to the sheep populations under study.

South Africa has a diverse array of phenotypically distinct and locally adapted sheep breeds that have been developed for different production systems ranging from mutton to wool and pelt, and some dual-purpose and nondescript breeds kept by smallholder farmers. This study investigated genetic diversity, population genetic structure and divergence between South African sheep breeds in order to gain an insight into breed history and genomic architecture aligned to breeding goals and production systems. The Illumina OvineSNP50 BeadChip was used to genotype 400 sheep belonging to 14 breeds representing mutton, pelt and mutton and wool dual-purpose breeds. Nguni sheep were included as a representative of indigenous nondescript breeds that are reared by smallholder farmers. Seeking a clearer understanding of the genetic diversity of South African breeds relative to global populations, 623 genotypes of sheep from worldwide populations were included in the analysis. These sheep breeds included six African, two Asian and eight European breeds. Across breeds, genetic diversity ranged from observed heterozygosity (H 0 ) = 0.26 ± 0.02 in Namaqua Afrikaner to H 0  = 0.38 ± 0.01 in Dohne Merino. The overall mean H 0 was 0.35 ± 0.04. The African and Asian populations were the most inbred populations with F IS ranging from 0.17 ± 0.05 in Grey Swakara and Ronderib Afrikaner sheep to 0.34 ± 0.07 in the Namaqua Afrikaner. The South African Dohne Merino (F IS  = 0.03 ± 0.01), SA Merino (F IS  = 0.05 ± 0.04) and Afrino (F IS  = 0.09 ± 0.02) and other global Merino-derived breeds were the least inbred. The first principal component explained 27.7% of the variation and separated the fat- and rump-tailed sheep (i.e. Swakara, Nguni, Blackhead Persian, Ethiopian Menzi, Meatmaster) from the Merino and Merino-derived breeds and the Dorset Horn. The second principal component separated the Merino and Merino-derived breeds from the English breed of Dorset Horn. Overall, South African indigenous breeds clustered together with indigenous breeds from other African and Asian countries. The optimal admixture cluster (K = 20) revealed various sources of within- and amongst-breed genomic variation associated with production purpose, adaptation and history of the breeds. The Blackhead Persian, Nguni and Namaqua Afrikaner breeds differed significantly from other breeds, particularly with the South African Mutton Merino and Dorset Horn. Breed-differentiating SNPs were observed within genomic regions associated with growth, adaptation and reproduction. Genes such as RAB44, associated with growth and meat/carcass traits, differentiated the Blackhead Persian from the Dorset Horn and South African Mutton Merino. The MAP2 and HRAS genes, which are associated with immune traits involving Toll-like receptors and Chemokine signalling pathways, differentiated the Nguni from the Dorset Horn. The current results give insight into the current status of the sheep genetic resources of South Africa relative to the global sheep population, highlighting both genetic similarities as well as divergence associated with production systems, geographical distribution and local adaptation., (© 2020 Stichting International Foundation for Animal Genetics.)

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.

Genetic diversity is of great importance and a prerequisite for genetic improvement and conservation programs in pigs and other livestock populations. The present study provides a genome wide analysis of the genetic variability and population structure of pig populations from different production systems in South Africa relative to global populations. A total of 234 pigs sampled in South Africa and consisting of village ( n = 91), commercial ( n = 60), indigenous ( n = 40), Asian ( n = 5) and wild ( n = 38) populations were genotyped using Porcine SNP60K BeadChip. In addition, 389 genotypes representing village and commercial pigs from America, Europe, and Asia were accessed from a previous study and used to compare population clustering and relationships of South African pigs with global populations. Moderate heterozygosity levels, ranging from 0.204 for Warthogs to 0.371 for village pigs sampled from Capricorn municipality in Eastern Cape province of South Africa were observed. Principal Component Analysis of the South African pigs resulted in four distinct clusters of (i) Duroc; (ii) Vietnamese; (iii) Bush pig and Warthog and (iv) a cluster with the rest of the commercial (SA Large White and Landrace), village, Wild Boar and indigenous breeds of Koelbroek and Windsnyer. The clustering demonstrated alignment with genetic similarities, geographic location and production systems. The PCA with the global populations also resulted in four clusters that where populated with (i) all the village populations, wild boars, SA indigenous and the large white and landraces; (ii) Durocs (iii) Chinese and Vietnamese pigs and (iv) Warthog and Bush pig. K = 10 (The number of population units) was the most probable ADMIXTURE based clustering, which grouped animals according to their populations with the exception of the village pigs that showed presence of admixture. AMOVA reported 19.92%-98.62% of the genetic variation to be within populations. Sub structuring was observed between South African commercial populations as well as between Indigenous and commercial breeds. Population pairwise F ST analysis showed genetic differentiation ( P ≤ 0.05) between the village, commercial and wild populations. A per marker per population pairwise F ST analysis revealed SNPs associated with QTLs for traits such as meat quality, cytoskeletal and muscle development, glucose metabolism processes and growth factors between both domestic populations as well as between wild and domestic breeds. Overall, the study provided a baseline understanding of porcine diversity and an important foundation for porcine genomics of South African populations., (Copyright © 2020 Hlongwane, Hadebe, Soma, Dzomba and Muchadeyi.)

Goats play a major role in poor marginalized communities of South Africa for food security and socio-economic purposes. Majority of the goats are raised in villages with poor infrastructure and resources, therefore facing challenges that affect growth performance which leads to low mature weights. Investigating growth profiles will shed light on growth performances and will aid in goat improvement and selection. This study investigated the growth profiles and genomic structure of SA indigenous breeds raised in different production systems to unravel the genetic potential of indigenous goat populations. Live weights and morphological body measurements were collected from a total of 83 kids representing the commercial meat-producing SA Boer (n = 14); the indigenous veld goats (IVG) of NC Skilder (n = 14), Mbuzi (n = 13), and Xhosa lob (n = 14) raised under intensive systems; and nondescript village goat populations (n = 14) raised in intensive and others (n = 14) raised in extensive production systems. The remaining 72 of 83 phenotyped goats were genotyped using the Illumina Caprine SNP50K BeadChip. The SA Boer had a higher weight (28.96 ± 0.30 kg) gain as compared to other populations. The Mbuzi population was the smallest (14.83 ± 0.33 kg), while the village goats raised in Pella Village were relatively smaller (17.55 ± 0.37 kg) than those raised on the research farm (19.55 ± 0.36 kg). The study concluded that both genetics and management systems can lead to improved growth performance in goat production. The outputs of this study can be used to identify suitable breeds and potential genotypes for optimal growth and establish optimal goat management systems suitable for communal farmers for improved productivity.

Nematodes of the genus Ascaridia are known to infect many species of birds and result in fatal diseases. A. galli damages the intestinal mucosa of chickens leading to blood loss, secondary infection and occasionally the obstruction of small intestines due to high worm burden. This study investigated the gene expression profiles in chickens from two different provinces of South Africa naturally exposed to A. galli infestations and tested either positive or negative for the parasite. The study further investigated gene expression profiles of the A. galli infected duodenum, jejunum and ileum tissues of the small intestines. The A. galli positive intestines displayed hypertrophy of the intestinal villi with accumulation of inflammatory cells and necrosis of the crypts of Lieberkühn glands, lesions that were absent in the uninfected intestines. Total RNA isolated from small intestines of infected and non-infected intestines was sequenced using Illumina HiSeq technology to generate up to 23,856,130 reads. Between any two-way comparisons of the intestines, 277 and 190 transcripts were significantly expressed in Limpopo and KwaZulu-Natal (KZN) chickens, respectively. Gene ontology analysis of the differentially expressed genes (DEGs) revealed an enrichment of genes reported to function in the immune response, defense response, inflammatory response and cell signalling genes. T cell receptor signalling pathways and arachidonic acid metabolism pathways were among the most significantly impacted pathways. Overall, the study provided insights into adaptative mechanisms for chickens extensively raised in parasite infected environments., (© 2019 The Authors.)

[This corrects the article DOI: 10.3389/fgene.2018.00163.].

Background: Bacillus endophyticus is a soil plant-endophytic bacterium, while B. anthracis is the causative agent of anthrax. The virulence factors of B. anthracis are the plasmid encoded tripartite toxins (pXO1) and poly-γ-glutamic acid (PGA) capsule (pXO2). B. endophyticus isolated alongside B. anthracis from animals that died of anthrax in Northern Cape Province (NCP), South Africa, harbored polyglutamate genes. The study compared the characteristics of B. anthracis and B. endophyticus with other Bacillus species with a focus on the presence of the PGA capsule or/and unbound PGA. The morphology and whole genome sequence analysis of B. endophyticus strains and B. anthracis were compared., Results: In conventional microbiology, B. endophyticus showed gram-positive round-shaped rods in single/short chains, which were endospore-forming, non-motile, non-haemolytic with white and dry colonies, and γ-phage resistant. B. anthracis was differentiated from B. endophyticus based on the latter's box-shaped rods in pairs/long chains, white-grey and slimy colonies, encapsulated and γ-phage susceptible. The study identified a PGA polyglutamate synthase operon that consisted of pgsBCA, γ-glutamyltranspeptidase (ggt) and pgsE in B. endophyticus genomes., Conclusions: PGA regions of B. anthracis contain capBCADE genes located in the pXO2 required for capsulation formation, while B. endophyticus contain the pgsBCAE genes in the chromosome. Whole genome and microbiology analysis identified B. endophyticus, as a non-capsuled endospore-forming bacterium that consists of PGA required for biosynthesis. B. endophyticus strains do not synthesize surface associated PGA, therefore capsule visualization of B. anthracis is a key diagnostic characteristic. The study highlights the significance of using whole genome shotgun sequencing to identify virulence and other important genes that might be present amongst unknown samples from natural outbreaks. None of the B. anthracis related plasmids or virulence genes were found in the B. endophyticus genomes.

Copy number variations (CNVs) comprise deletions, duplications, and insertions found within the genome larger than 50 bp in size. CNVs are thought to be primary role-players in breed formation and adaptation. South Africa boasts a diverse ecology with harsh environmental conditions and a broad spectrum of parasites and diseases that pose challenges to livestock production. This has led to the development of composite cattle breeds which combine the hardiness of Sanga breeds and the production potential of the Taurine breeds. The prevalence of CNVs within these respective breeds of cattle and the prevalence of CNV regions (CNVRs) in their diversity, adaptation and production is however not understood. This study therefore aimed to ascertain the prevalence, diversity, and correlations of CNVRs within cattle breeds used in South Africa. Illumina Bovine SNP50 data and PennCNV were utilized to identify CNVRs within the genome of 287 animals from seven cattle breeds representing Sanga, Taurine, Composite, and cross breeds. Three hundred and fifty six CNVRs of between 36 kb to 4.1 Mb in size were identified. The null hypothesis that one CNVR loci is independent of another was tested using the GENEPOP software. One hunded and two and seven of the CNVRs in the Taurine and Sanga/Composite cattle breeds demonstrated a significant ( p ≤ 0.05) association. PANTHER overrepresentation analyses of correlated CNVRs demonstrated significant enrichment of a number of biological processes, molecular functions, cellular components, and protein classes. CNVR genetic variation between and within breed group was measured using phiPT which allows intra-individual variation to be suppressed and hence proved suitable for measuring binary CNVR presence/absence data. Estimate PhiPT within and between breed variance was 2.722 and 0.518 respectively. Pairwise population PhiPT values corresponded with breed type, with Taurine Holstein and Angus breeds demonstrating no between breed CNVR variation. Phylogenetic trees were drawn. CNVRs primarily clustered animals of the same breed type together. This study successfully identified, characterized, and analyzed 356 CNVRs within seven cattle breeds. CNVR correlations were evident, with many more correlations being present among the exotic Taurine breeds. CNVR genetic diversity of Sanga, Taurine and Composite breeds was ascertained with breed types exposed to similar selection pressures demonstrating analogous incidences of CNVRs.

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.

Expansion of goat improvement programs requires exploration of the factors that influence the production system and breeding initiatives. Characterization of goat breeds or populations is crucial in providing information on prevalent goat types and their attributes and may suffice as a guideline on conservation, development, and selection for improved productivity. This study investigated the existing village goat production system and phenotypic diversity of the different village populations from four South African provinces. The study further investigated the use of phenotypic attributes to classify goats to breeds or populations. Data was collected from 142 households in 26 villages of the Eastern Cape (n = 2 villages), KwaZulu-Natal (n = 6 villages), Limpopo (n = 13 villages), and North West (n = 5 villages) provinces through a survey. Individual interviews and focus group discussions revealed that the mean goat herd size per household was least in Limpopo at 13.2 ± 12.40 and highest in Eastern Cape (34.18 ± 28.36). Flocks had more (p < 0.05) adults than kids, and the distribution of breeding animals was biased towards does and less bucks. Goats were kept mainly for meat, for selling, and for ritual ceremonies. The goat production system was mainly scavenging. Goat health was a major challenge across households and villages. Qualitative traits such coat, horn, ears, and wattle characteristics were recorded for populations of village goats (n = 319) and a feral Tankwa breed (n = 50). The dominant coat pattern was plain (74.53%) with black as the most common coat color (31.98%). Across breeds, a majority (88.08%) of the goats had horns, and 7.59% had wattles while 56.64% had beard. Adult goats (<3 years; n = 398) were further analyzed for five quantitative traits of chest girth, height, length, and pin bone and there were significant (p < 0.05) breed differences in all. A stepwise discriminatory procedure was used to rank the quantitative traits according to their discriminatory power to separate breeds or populations. Significant traits were then subjected to canonical discriminant analysis for principle component analysis. Based on the quantitative traits, the Tankwa, Xhosa, and Tswana goats formed their own cluster separated from commercial meat type breeds and the Venda and Zulu ecotypes. The discriminant function analysis correctly classified 90.41% of the Zulu goats and 82.93% of the Xhosa village populations. None of the Savanna goats were correctly classified. The study demonstrated diversity in village goat populations and production systems, which would allow for within population selection in genetic improvement programs. The heterogeneity in the phenotypic traits in the village goats is reflective of the role of village production systems in the maintenance of animal diversity in local populations.

Goats are an important source of livelihood especially in smallholder communities. Infections with gastrointestinal nematodes (GIN) remain the most prevalent parasitic diseases affecting small ruminants. The study was conducted to assess management, the level of knowledge and control of gastrointestinal nematodes. Surveys were conducted in Chipinge, Shurugwi, Binga, Tsholotsho and Matobo districts, representing the five natural/agro-ecological regions (NR) in Zimbabwe. Data was collected in 135 households using a pre-tested semi-structured questionnaire. Results indicated that goats were ranked the most important livestock species, with high flock sizes in NR IV and V. Partitioning of roles was such that the adult males were involved in decision-making while females and children were involved in day-to-day management of animals. Farmers showed low levels of input use, with natural pasture (98.4%) being the main feed source and indigenous breeds (73.2%) being kept. Farmers ranked food and financial benefits as the main reasons for keeping goats. Gastrointestinal nematodes ranked the highest as the most common disease, with majority of farmers (57%) not controlling or treating animals and 63% of farmers not having knowledge on the spread of GIN. Access to veterinary services, anthelmintic class used and breeds used by the farmers had the highest effects on parasitic infections in households. Farmer education is required for capacitation of farmer in terms of disease prevention and control so as to improve goat production.

Background: Understanding the history of cattle breeds is important because it provides the basis for developing appropriate selection and breed improvement programs. In this study, patterns of ancestry and admixture in Afrikaner, Nguni, Drakensberger and Bonsmara cattle of South Africa were investigated. We used 50 K single nucleotide polymorphism genotypes that were previously generated for the Afrikaner (n = 36), Nguni (n = 50), Drakensberger (n = 47) and Bonsmara (n = 44) breeds, and for 394 reference animals representing European taurine, African taurine, African zebu and Bos indicus., Results and Discussion: Our findings support previous conclusions that Sanga cattle breeds are composites between African taurine and Bos indicus. Among these breeds, the Afrikaner breed has significantly diverged from its ancestral forebears, probably due to genetic drift and selection to meet breeding objectives of the breed society that enable registration. The Nguni, Drakensberger and Bonsmara breeds are admixed, perhaps unintentionally in the case of Nguni and Drakensberger, but certainly by design in the case of Bonsmara, which was developed through crossbreeding between the Afrikaner, Hereford and Shorthorn breeds., Conclusions: We established patterns of admixture and ancestry for South African Sanga cattle breeds, which provide a basis for developing appropriate strategies for their genetic improvement.

A longitudinal study was conducted in low-input low-output farming systems to determine the prevalence of gastrointestinal parasitic infections in different age groups, sex and associated risk factors in goats. A total of 580 indigenous goats were randomly selected in areas representing the five agro-ecological regions of Zimbabwe in the dry and wet seasons. Blood and faecal samples were collected from each animal and egg/oocyst per gram of faeces (epg/opg), larval culture, and packed cell volumes (PCV) were determined. Factors affecting parasitic infections were evaluated. Highest prevalence was determined for Eimeria oocysts (43%), strongyles (31%) and lower levels in trematodes and cestodes. Parasites identified were Haemonchus, Strongyloides and Oesophagostomum. Area, season, sex and age significantly influenced patterns of gastrointestinal infections ( P   <  0.05). Cannonical correlations indicated that parasite species composition varied by area and impacts of risk factors also differed. Risk of infection was very high for goats sampled in Natural regions (NR) I, II, III (OR = 6.6-8.2; P  < 0.05) as compared to those in NR IV and V. Highest helminths and Eimeria infections were observed in the wet vs. dry season ( P   <  0.05). Young animals were more susceptible to parasitic infections ( P   <  0.05). Prevalence was higher in males than females, with odds of infection for males being almost three times to that for females ( P   <  0.0001). Knowledge concerning gastrointestinal helminth biology and epidemiological infection patterns caused by these parasites is essential in the development of appropriate control strategies and this has a potential to reduce production losses.

Introduction: Bacillus anthracis is the causative agent of anthrax, a disease endemic in regions of Northern Cape Province and Kruger National Park of South Africa. Accurate identification of virulent B. anthracis is essential but challenging due to its close relationship with other members of B. cereus group. This study characterized B. anthracis and Bacillus species that were recovered from animals and the environment where animals died of anthrax symptoms in southern Africa using a polyphasic approach., Methodology: For this purpose, 3 B. anthracis and 10 Bacillus isolates were subjected to microbiology tests, BiologOmniLog identification system (Biolog), 16S ribosomal RNA (rRNA) sequence analysis, polymerase chain reaction (PCR) detection of protective antigen (pag) and capsule (cap) regions, and real-time PCR using hybridization probes targeting chromosomal, pag, and capC genes., Results: The Bacillus isolates were non-hemolytic, non-motile, and susceptible to penicillin, which is typical of B. anthracis, but resistant to gamma phage, unlike typical B. anthracis. The Biolog system and 16S rRNA gene sequence analysis identified most of the Bacillus isolates as B. endophyticus (7 of 10). Conventional PCR revealed that most of the Bacillus isolates contained capBCA gene regions. This highlights the limitation of the specificity of conventional PCR and the fact that the real-time PCR is more specific and reliable for anthrax diagnosis., Conclusions: Real-time PCR, 16S rRNA sequencing, and confirmatory microbiology tests including phage resistance distinguished Bacillus isolates from B. anthracis in this study. Identification of B. anthracis should be done using a polyphasic approach.

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., (© 2016 Stichting International Foundation for Animal Genetics.)

The control of gastrointestinal nematodes (GIN) is mainly based on the use of drugs, grazing management, use of copper oxide wire particles and bioactive forages. Resistance to anthelmintic drugs in small ruminants is documented worldwide. Host genetic resistance to parasites, has been increasingly used as a complementary control strategy, along with the conventional intervention methods mentioned above. Genetic diversity in resistance to GIN has been well studied in experimental and commercial flocks in temperate climates and more developed economies. However, there are very few report outputs from the more extensive low-input/output smallholder systems in developing and emerging countries. Furthermore, results on quantitative trait loci (QTL) associated with nematode resistance from various studies have not always been consistent, mainly due to the different nematodes studied, different host breeds, ages, climates, natural infections versus artificial challenges, infection level at sampling periods, among others. The increasing use of genetic markers (Single Nucleotide Polymorphisms, SNPs) in GWAS or the use of whole genome sequence data and a plethora of analytic methods offer the potential to identify loci or regions associated nematode resistance. Genomic selection as a genome-wide level method overcomes the need to identify candidate genes. Benefits in genomic selection are now being realised in dairy cattle and sheep under commercial settings in the more advanced countries. However, despite the commercial benefits of using these tools, there are practical problems associated with incorporating the use of marker-assisted selection or genomic selection in low-input/output smallholder farming systems breeding schemes. Unlike anthelmintic resistance, there is no empirical evidence suggesting that nematodes will evolve rapidly in response to resistant hosts. The strategy of nematode control has evolved to a more practical manipulation of host-parasite equilibrium in grazing systems by implementation of various strategies, in which improvement of genetic resistance of small ruminant should be included. Therefore, selection for resistant hosts can be considered as one of the sustainable control strategy, although it will be most effective when used to complement other control strategies such as grazing management and improving efficiency of anthelmintics currently., (Copyright © 2016 The Authors. Published by Elsevier B.V. All rights reserved.)

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., (© 2016 Stichting International Foundation for Animal Genetics.)

The majority of rural households in developing countries own village chickens that are reared under traditional scavenging systems with few inputs and exposure to various parasitic infestations. Understanding of the village chicken farming system and its influence on helminth infestation is a prerequisite for optimal prevention and control strategies. This study investigated the village chicken production system and associated gastrointestinal parasites in 87 households from Limpopo (n = 39) and KwaZulu-Natal (n = 48) provinces of South Africa. A total of 191 village chicken faecal samples and 145 intestines were collected to determine the prevalence of gastrointestinal parasites in villages of Limpopo and KwaZulu-Natal provinces, respectively. The faecal floatation analysis of samples from Limpopo and KwaZulu-Natal provinces indicated infestations by Ascaridia galli (18.77%), Heterakis gallinarum (15.56%) and Capillaria spp. (4.00%); tapeworms Choanotaenia infundibulum (2.10%) and Raillietina cesticillus (6.00%) and Eimeria spp. (29.46%). Mixed infestations were observed in five (4.90%) samples from Limpopo province and in only four (4.49%) from KwaZulu-Natal province, of which 1.12% were a mixture of C. infundibulum and Eimeria spp. and 3.37% a combination of H. gallinarum and Eimeria spp. In Limpopo, 2.94% of the chickens were positive for H. gallinarum and Eimeria spp., whilst 0.98% had A. galli and Capillaria spp. infestations. Further investigation is needed to understand the impact of gastrointestinal parasites on village chicken health and production and develop appropriate intervention and control strategies feasible for smallholder farmers.

Africa is host to diverse and locally adapted cattle breeds that are expected to survive the harsh and extreme tropical environments associated with diseases and parasite infections, heat stress and episodes of feed and water scarcity. Genomic copy number variations (CNVs) are considered to be primary role players in cattle breed formation and adaptation where isolation and genetic drift together with subsequent mutations have created an enormous diversity of local populations. CNVs are modifications in DNA structure comprising deletions, duplications and insertions that are >1 kb in size. Despite attracting much attention, the frequency and pattern of bovine CNV events, especially in African cattle breeds, are for the most part largely unknown. Characterization of genetic variation in the indigenous cattle of Africa will be a vital step toward dissecting the molecular mechanisms underlying phenotypic variation and local adaptation. This review therefore aims to describe the current knowledge regarding bovine CNVs and the implications and potentials they encompass for dissecting genetic adaptation and the genotypic skeleton of tropical African cattle populations., (© 2015 Stichting International Foundation for Animal Genetics.)

The present study investigated the seroprevalence of antibodies to Ehrlichia ruminantium and the associated risk factors in goats from five different farming provinces of South Africa. Sera collected from 686 goats of the commercial meat type (n=179), mohair type (n=9), non-descript indigenous goats from Eastern Cape (n=56), KwaZulu-Natal (n=209), Limpopo (n=111), North West (n=61) and Northern Cape (n=11) provinces and a feral Tankwa goat (n=50) were tested for the presence of immunoglobulin G (IgG) antibodies to antigens of E. ruminantium using the indirect fluorescent-antibody test (IFAT). Fifty two percent of these goats had ticks. The overall seroprevalence of antibodies to E. ruminantium was 64.87% (445/686) with the highest seroprevalence reported for Limpopo (95.50%) and lowest for Northern Cape (20.29%). Highest seroprevalence for antibodies to E. ruminantium was observed in goats from endemic regions (76.09%), and from smallholder production systems (89.54%). High seroprevalence was also observed in non-descript indigenous goats (85.04%), adult goat (69.62%), in does (67.46%) and goats infested with ticks (85.79%). The logistic model showed a gradient of increasing risk for commercial meat type Savanna (OR=3.681; CI=1.335-10.149) and non-descript indigenous (OR=3.466; CI=1.57-7.645) compared to Boer goats and for goats from the smallholder production system (OR=2.582; CI=1.182-5.639) and those with ticks (OR=3.587; CI=2.105-6.112). Results from this study showed that E. ruminantium infections were prevalent but were widely and unevenly distributed throughout South Africa. Findings from the study facilitate identification and mapping of risk areas for heartwater and its endeminicity in South Africa and should be taken into consideration for future disease control strategies and local goat improvement programs., (Copyright © 2016 Elsevier B.V. All rights reserved.)

Ascaridia galli is one of the most common nematode affecting chickens. This study characterized A. galli parasites collected from South African village chickens of Limpopo (n=18) and KwaZulu-Natal (n=22) provinces using the 510bp sequences of cytochrome C oxidase subunit 1 gene of the mitochondrial DNA. Fourteen and 12 polymorphic sites were observed for Limpopo and KwaZulu-Natal sequences, respectively. Six haplotypes were observed in total. Haplotype diversity was high and ranged from 0.749 for Limpopo province to 0.758 for KwaZulu-Natal province isolates. There was no genetic differentiation between A. galli from Limpopo and KwaZulu-Natal provinces. The six South African haplotypes were unique compared to those published in the GeneBank sampled from Hy-line chickens raised under organic farming in Denmark. The utility of cytochrome C oxidase subunit 1 gene as a potential genetic marker for studying A. galli in village chicken populations is presented., (Copyright © 2015 Elsevier B.V. All rights reserved.)

Knowledge on the extent of linkage disequilibrium (LD) in livestock populations is essential to determine the minimum distance between markers required for effective coverage when conducting genome-wide association studies (GWAS). This study evaluated the extent of LD, persistence of allelic phase and effective population size (Ne) for four Sanga cattle breeds in South Africa including the Afrikaner (n = 44), Nguni (n = 54), Drakensberger (n = 47), and Bonsmara breeds (n = 46), using Angus (n = 31) and Holstein (n = 29) as reference populations. We found that moderate LD extends up to inter-marker distances of 40-60 kb in Angus (0.21) and Holstein (0.21) and up to 100 kb in Afrikaner (0.20). This suggests that genomic selection and association studies performed within these breeds using an average inter-marker r (2)≥ 0.20 would require about 30,000-50,000 SNPs. However, r (2)≥ 0.20 extended only up to 10-20 kb in the Nguni and Drakensberger and 20-40 kb in the Bonsmara indicating that 75,000 to 150,000 SNPs would be necessary for GWAS in these breeds. Correlation between alleles at contiguous loci indicated that phase was not strongly preserved between breeds. This suggests the need for breed-specific reference populations in which a much greater density of markers should be scored to identify breed specific haplotypes which may then be imputed into multi-breed commercial populations. Analysis of effective population size based on the extent of LD, revealed Ne = 95 (Nguni), Ne = 87 (Drakensberger), Ne = 77 (Bonsmara), and Ne = 41 (Afrikaner). Results of this study form the basis for implementation of genomic selection programs in the Sanga breeds of South Africa.

Background: The detection of selection signatures in breeds of livestock species can contribute to the identification of regions of the genome that are, or have been, functionally important and, as a consequence, have been targeted by selection., Methods: This study used two approaches to detect signatures of selection within and between six cattle breeds in South Africa, including Afrikaner (n = 44), Nguni (n = 54), Drakensberger (n = 47), Bonsmara (n = 44), Angus (n = 31) and Holstein (n = 29). The first approach was based on the detection of genomic regions in which haplotypes have been driven towards complete fixation within breeds. The second approach identified regions of the genome that had very different allele frequencies between populations (F ST)., Results and Discussion: Forty-seven candidate genomic regions were identified as harbouring putative signatures of selection using both methods. Twelve of these candidate selected regions were shared among the breeds and ten were validated by previous studies. Thirty-three of these regions were successfully annotated and candidate genes were identified. Among these genes the keratin genes (KRT222, KRT24, KRT25, KRT26, and KRT27) and one heat shock protein gene (HSPB9) on chromosome 19 between 42,896,570 and 42,897,840 bp were detected for the Nguni breed. These genes were previously associated with adaptation to tropical environments in Zebu cattle. In addition, a number of candidate genes associated with the nervous system (WNT5B, FMOD, PRELP, and ATP2B), immune response (CYM, CDC6, and CDK10), production (MTPN, IGFBP4, TGFB1, and AJAP1) and reproductive performance (ADIPOR2, OVOS2, and RBBP8) were also detected as being under selection., Conclusions: The results presented here provide a foundation for detecting mutations that underlie genetic variation of traits that have economic importance for cattle breeds in South Africa.

Bacillus anthracis is a Gram-positive bacterium that causes anthrax, mainly in herbivores through exotoxins and capsule produced on plasmids, pXO1 and pXO2. This paper compares the whole-genome sequences of two B. anthracis strains from an endemic region and a sporadic outbreak in South Africa. Sequencing was done using next-generation sequencing technologies., (Copyright © 2015 Lekota et al.)

Background: Copy number variations (CNVs) are modifications in DNA structure comprising of deletions, duplications, insertions and complex multi-site variants. Although CNVs are proven to be involved in a variety of phenotypic discrepancies, the full extent and consequence of CNVs is yet to be understood. To date, no such genomic characterization has been performed in indigenous South African Nguni cattle. Nguni cattle are recognized for their ability to sustain harsh environmental conditions while exhibiting enhanced resistance to disease and parasites and are thought to comprise of up to nine different ecotypes., Methods: Illumina BovineSNP50 Beadchip data was utilized to investigate genomic population structure and the prevalence of CNVs in 492 South African Nguni cattle. PLINK, ADMIXTURE, R, gPLINK and Haploview software was utilized for quality control, population structure and haplotype block determination. PennCNV hidden Markov model identified CNVs and genes contained within and 10 Mb downstream from reported CNVs. PANTHER and Ensembl databases were subsequently utilized for gene annotation analyses., Results: Population structure analyses on Nguni cattle revealed 5 sub-populations with a possible sub-structure evident at K equal to 8. Four hundred and thirty three CNVs that formed 334 CNVRs ranging from 30 kb to 1 Mb in size are reported. Only 231 of the 492 animals demonstrated CNVRs. Two hundred and eighty nine genes were observed within CNVRs identified. Of these 149, 28, 44, 2 and 14 genes were unique to sub-populations A, B, C, D and E respectively. Gene ontology analyses demonstrated a number of pathways to be represented by respective genes, including immune response, response to abiotic stress and biological regulation processess., Conclusions: CNVs may explain part of the phenotypic diversity and the enhanced adaptation evident in Nguni cattle. Genes involved in a number of cellular components, biological processes and molecular functions are reported within CNVRs identified. The significance of such CNVRs and the possible effect thereof needs to be ascertained and may hold interesting insight into the functional and adaptive consequence of CNVs in cattle.

This study investigated the haplo-block structure, haplotype sharing, and diversity in extensively raised chicken populations of Southern Africa. Two hundred ninety village chickens from Malawi (N = 30), South Africa (N = 132), and Zimbabwe (N = 128) were included in the study, from which 649, 2104, and 2442 haplo-blocks were observed, respectively. The majority of haplo-blocks were smaller than 25 kb in size and only five blocks were more than 2000 kb in size. The low chromosomal coverage of haplo-blocks observed across the genome suggests that multiple recombination events fragmented the ancestral haplo-blocks into smaller sizes. Haplo-block sharing was observed between populations with 2325 haplo-blocks common between Zimbabwe and Malawi and 2689 between South Africa and Zimbabwe. Haplotype sharing allows transferability of genomic tools between these extensively raised chicken populations of Southern Africa. The unique haplo-blocks could have originated from isolated evolution taking place in specific agro-ecological zones. Quantitative trait loci analysis revealed that genes related to body composition were spanned by these haplo-blocks. Body composition traits are important for village chicken populations, which have to harness poor quality feed obtained from the environment to meet their maintenance and production needs.

Extensively raised village chickens are considered a valuable source of biodiversity, with genetic variability developed over thousands of years that ought to be characterized and utilized. Surveys that can reveal a population's genetic structure and provide an insight into its demographic history will give valuable information that can be used to manage and conserve important indigenous animal genetic resources. This study reports population diversity and structure, linkage disequilibrium and effective population sizes of Southern African village chickens and conservation flocks from South Africa. DNA samples from 312 chickens from South African village and conservation flocks (n = 146), Malawi (n = 30) and Zimbabwe (n = 136) were genotyped using the Illumina iSelect chicken SNP60K BeadChip. Population genetic structure analysis distinguished the four conservation flocks from the village chicken populations. Of the four flocks, the Ovambo clustered closer to the village chickens particularly those sampled from South Africa. Clustering of the village chickens followed a geographic gradient whereby South African chickens were closer to those from Zimbabwe than to chickens from Malawi. Different conservation flocks seemed to have maintained different components of the ancestral genomes with a higher proportion of village chicken diversity found in the Ovambo population. Overall population LD averaged over chromosomes ranged from 0.03 ± 0.07 to 0.58 ± 0.41 and averaged 0.15 ± 0.16. Higher LD, ranging from 0.29 to 0.36, was observed between SNP markers that were less than 10 kb apart in the conservation flocks. LD in the conservation flocks steadily decreased to 0.15 (PK) and 0.24 (VD) at SNP marker interval of 500 kb. Genomewide LD decay in the village chickens from Malawi, Zimbabwe and South Africa followed a similar trend as the conservation flocks although the mean LD values for the investigated SNP intervals were lower. The results suggest low effective population sizes particularly in the conservation flocks. The utility and limitations of the iselect chicken SNP60K in village chicken populations is discussed.

Information about genetic diversity and population structure among cattle breeds is essential for genetic improvement, understanding of environmental adaptation as well as utilization and conservation of cattle breeds. This study investigated genetic diversity and the population structure among six cattle breeds in South African (SA) including Afrikaner (n = 44), Nguni (n = 54), Drakensberger (n = 47), Bonsmara (n = 44), Angus (n = 31), and Holstein (n = 29). Genetic diversity within cattle breeds was analyzed using three measures of genetic diversity namely allelic richness (AR), expected heterozygosity (He) and inbreeding coefficient (f). Genetic distances between breed pairs were evaluated using Nei's genetic distance. Population structure was assessed using model-based clustering (ADMIXTURE). Results of this study revealed that the allelic richness ranged from 1.88 (Afrikaner) to 1.73 (Nguni). Afrikaner cattle had the lowest level of genetic diversity (He = 0.24) and the Drakensberger cattle (He = 0.30) had the highest level of genetic variation among indigenous and locally-developed cattle breeds. The level of inbreeding was lower across the studied cattle breeds. As expected the average genetic distance was the greatest between indigenous cattle breeds and Bos taurus cattle breeds but the lowest among indigenous and locally-developed breeds. Model-based clustering revealed some level of admixture among indigenous and locally-developed breeds and supported the clustering of the breeds according to their history of origin. The results of this study provided useful insight regarding genetic structure of SA cattle breeds.

Enterotoxigenic Escherichia coli (ETEC) and shiga toxin E. coli (STEC) are important causes of colibacillosis in piglets. Recently, enteroaggregative E. coli heat-stable enterotoxin 1 (EAST-1) has been implicated in pig diarrhoea. This study investigated the prevalence of enterotoxin [heat-labile toxins (LT), heat-stable toxin a (STa), heat-stable toxin b (STb)], shiga toxins (Stx1, Stx2, Stx2e), enteroaggregative heat-stable E. coli (EAST-1), associated fimbriae (F4, F5, F6, F41, F18ab, F18ac) and non-fimbrial adhesins [adhesin involved in diffuse adherence 1 (AIDA-1), attaching and effacing factor, porcine attaching- and effacing-associated factor] in South African pigs. A total of 263 E. coli strains were isolated from Landrace (n = 24), Large White (n = 126), Duroc (n = 28) and indigenous (n = 85) breeds of piglets aged between 9 and 136 days. PCR was used in the analysis. Virulent genes were detected in 40.3% of the isolates, of which 18.6, 0.4 and 17.5% were classified as ETEC, STEC and enteroaggregative E. coli (EAEC), respectively. Individual genes were found in the following proportions: STb (19.01%), LT (0.4%), STa (3.4%), St2xe (1.1%) and EAST-1 (20.2%) toxins. None of the tested fimbriae were detected in ETEC and STEC isolates. About one third of the ETEC and STEC isolates was tested negative for both fimbrial and non-fimbrial adhesins. Twenty-five pathotypes from ETEC-, EAEC- and STEC-positive strains were identified. Pathotypes EAST-1 (30.2%), STb (13.2%) and STb/AIDA-1 (10.4%) were most prevalent. The study provided insight on possible causes of colibacillosis in South African pigs.