Selection signature analysis and genome-wide divergence of South African Merino breeds from their founders

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.

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

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.

Runs of homozygosity analysis of South African sheep breeds from various production systems investigated using OvineSNP50k data

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.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07314-2.

The genomic architecture of South African mutton, pelt, dual-purpose and nondescript sheep breeds relative to global sheep populations

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.26 ± 0.02 in Namaqua Afrikaner to H₀  = 0.38 ± 0.01 in Dohne Merino. The overall mean H₀ 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.

Relationship between population genomic structure and growth profiles of South African goats under different production systems

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.

Gene expression profiles of the small intestine of village chickens from an Ascaridia galli infested environment

•

The transcriptome of chickens from parasite infested environment was sequenced.

•

Different genes were reported between A. galli infected and non-infected chickens.

•

Upregulated immune and inflammatory response genes are associated with fighting parasites.

•

T cell receptor signalling and arachidonic acid metabolism pathways were impacted.

•

Different segments of the intestines differed in gene expression and associated pathways.

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.

Management and control of gastrointestinal nematodes in communal goat farms in Zimbabwe

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.

Breeding for resistance to gastrointestinal nematodes - the potential in low-input/output small ruminant production systems

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.

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

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

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

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

Population genetic structure of Ascaridia galli of extensively raised chickens of South Africa

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.

Haplo-block structure of Southern African village chicken populations inferred using genome-wide SNP data

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.

Mitochondrial DNA D-loop sequences suggest a Southeast Asian and Indian origin of Zimbabwean village chickens

This study sought to assess mitochondrial DNA (mtDNA) diversity and phylogeographic structure of chickens from five agro-ecological zones of Zimbabwe. Furthermore, chickens from Zimbabwe were compared with populations from other geographical regions (Malawi, Sudan and Germany) and other management systems (broiler and layer purebred lines). Finally, haplotypes of these animals were aligned to chicken sequences, taken from GenBank, that reflected populations of presumed centres of domestication. A 455-bp fragment of the mtDNA D-loop region was sequenced in 283 chickens of 14 populations. Thirty-two variable sites that defined 34 haplotypes were observed. In Zimbabwean chickens, diversity within ecotypes accounted for 96.8% of the variation, indicating little differentiation between ecotypes. The 34 haplotypes clustered into three clades that corresponded to (i) Zimbabwean and Malawian chickens, (ii) broiler and layer purebred lines and Northwest European chickens, and (iii) a mixture of chickens from Zimbabwe, Sudan, Northwest Europe and the purebred lines. Diversity among clades explained more than 80% of the total variation. Results indicated the existence of two distinct maternal lineages evenly distributed among the five Zimbabwean chicken ecotypes. For one of these lineages, chickens from Zimbabwe and Malawi shared major haplotypes with chicken populations that have a Southeast Asian background. The second maternal lineage, probably from the Indian subcontinent, was common to the five Zimbabwean chicken ecotypes, Sudanese and Northwest European chickens as well as purebred broiler and layer chicken lines. A third maternal lineage excluded Zimbabwean and other African chickens and clustered with haplotypes presumably originating from South China.

Variation in village chicken production systems among agro-ecological zones of Zimbabwe

The degree to which village chickens are integrated in the smallholder farming systems differs depending on the socio-economic, cultural and biological factors within each system. The objective of this study was to characterise the village chicken farming systems and identify possible threats to, and opportunities for, local chickens in the agro-ecological zones of Zimbabwe. A pre-tested questionnaire was administered to households randomly selected from five districts, Risitu (n=97), Hurungwe (n=56), Gutu (n=77), Gokwe-South (n=104) and Beitbridge (n=37) in eco-zones I-V, respectively. Age of head of household averaged 47 years (SD = 14.3). Land holdings per household averaged 4.82 ha (SD = 3.6). Overall, 17.7 percent of the households ranked livestock as the major source of income compared to 70.8 percent who ranked crops as the main contributor. Chicken flock size averaged 16.7 (SD = 12.4), and the highest flock sizes were observed in eco-zones I and IV. Households owning cattle, goats and other livestock assigned less important ranks to chickens. Chickens were usedmainly for the provision of meat and eggs whilst the use of chicken feathers and investment were uncommon practises. Results indicate that more support is necessary for village chickens in the non-cropping regions of the country.

Absence of population substructuring in Zimbabwe chicken ecotypes inferred using microsatellite analysis

The objective of this study was to investigate the population structure of village chickens found in the five agro-ecological zones of Zimbabwe. Twenty-nine microsatellites were genotyped for chickens randomly selected from 13 populations, including the five eco-zones of Zimbabwe (n = 238), Malawi (n = 60), Sudan (n = 48) and six purebred lines (n = 180). A total of 280 alleles were observed in the 13 populations. Forty-eight of these alleles were unique to the Zimbabwe chicken ecotypes. The average number (+/-SD) of alleles/locus was 9.7 +/- 5.10. The overall heterozygote deficiency in the Zimbabwe chickens (F(IT) +/- SE) was 0.08 +/- 0.01, over 90% of which was due to within-ecotype deficit (F(IS)). Small Nei's standard genetic distances ranging from 0.02 to 0.05 were observed between Zimbabwe ecotypes compared with an average of 0.6 between purebred lines. The structure software program was used to cluster individuals to 2 </= K </= 7 assumed clusters. The most probable clustering was found at K = 6. Ninety-seven of 100 structure runs were identical, in which Malawi, Sudan and purebred lines split out as independent clusters and the five Zimbabwe ecotypes clustered into one population. The within-ecotype marker-estimated kinships (mean = 0.13) differed only slightly from the between-ecotype estimates. Results from this study lead to a rejection of the hypothesis that village chickens are substructured across agro-ecological zones but indicated high genetic diversity within the Zimbabwe chicken population.

Village chicken flock dynamics and the contribution of chickens to household livelihoods in a smallholder farming area in Zimbabwe

The development of successful production strategies for poultry rearing depends on an accurate description of village chicken production systems. In Rushinga District of Zimbabwe, participatory rural appraisals (PRAs) followed by checklists and intensive case studies were carried out in three villages. The role of chickens in the livelihoods of households was evaluated. Flock dynamics were monitored monthly for 24 months. Women who were resident on the farm headed 19% of the households. A household comprised 4.8 +/- 2.5 members with arable land of approximately 2.6 ha. In addition to chickens, households grew maize, cotton and sunflower and kept large animals. Flocks ranging from 1 to 50 village chickens per household were reared under a scavenging system of management with suboptimal housing, inadequate feeding and poor health care. The use of ethno-veterinary medicine was common in treating sick chickens. The largest flock sizes were observed in the hot-wet season. Over 90% of an average of 5.4 entries/household per month were from hatched chicks. Mortality claimed an average of 80% of the total exits. Chicken production potential (CPP), which defined the proportion of chickens that could be utilized by a household, averaged 50%. Chicken production efficiency (CPE) was approximately 15% of the CPP. Egg consumption patterns were low and similar across seasons.