Ascertaining gene flow patterns in livestock populations of developing countries: a case study in Burkina Faso goat

Investigating the footprint of post-domestication dispersal on the diversity of modern European, African and Asian goats

BACKGROUND

Goats were domesticated in the Fertile Crescent about 10,000 years before present (YBP) and subsequently spread across Eurasia and Africa. This dispersal is expected to generate a gradient of declining genetic diversity with increasing distance from the areas of early livestock management. Previous studies have reported the existence of such genetic cline in European goat populations, but they were based on a limited number of microsatellite markers. Here, we have analyzed data generated by the AdaptMap project and other studies. More specifically, we have used the geographic coordinates and estimates of the observed (Ho) and expected (He) heterozygosities of 1077 European, 1187 African and 617 Asian goats belonging to 38, 43 and 22 different breeds, respectively, to find out whether genetic diversity and distance to Ganj Dareh, a Neolithic settlement in western Iran for which evidence of an early management of domestic goats has been obtained, are significantly correlated.

RESULTS

Principal component and ADMIXTURE analyses revealed an incomplete regional differentiation of European breeds, but two genetic clusters representing Northern Europe and the British-Irish Isles were remarkably differentiated from the remaining European populations. In African breeds, we observed five main clusters: (1) North Africa, (2) West Africa, (3) East Africa, (4) South Africa, and (5) Madagascar. Regarding Asian breeds, three well differentiated West Asian, South Asian and East Asian groups were observed. For European and Asian goats, no strong evidence of significant correlations between Ho and He and distance to Ganj Dareh was found. In contrast, in African breeds we detected a significant gradient of diversity, which decreased with distance to Ganj Dareh.

CONCLUSIONS

The detection of a genetic cline associated with distance to the Ganj Dareh in African but not in European or Asian goat breeds might reflect differences in the post-domestication dispersal process and subsequent migratory movements associated with the management of caprine populations from these three continents.

Ancient autozygous segments subject to positive selection suggest adaptive immune responses in West African cattle

Small-sized and trypanotolerant West African taurine (Bos taurus) cattle are a unique case of human-mediated process of adaptation to a challenging environment. Extensive gene flow with Sahelian zebu (B. indicus), bigger and with some resistance to tick attack, occurred for centuries and allowed the apparition of stable crossbred populations (sanga) having intermediate characteristics. Up to 237 individuals belonging to 10 different taurine, zebu and sanga cattle populations sampled in Benin, Burkina Faso and Niger were typed using the BovineHD BeadChip of Illumina to identify signatures of selection, assessed using three different Extended-Haplotype-Homozygosity-based statistics, overlapping with ancient, originated 1024 or 2048 generations ago, Homozygosity-By-Descent segments in the cattle genome. Candidate genomic regions were defined ensuring their importance within cattle type and using zebu as reference. Functional annotation analysis identified four statistically significant Annotation Clusters in taurine cattle (from ACt1 to ACt4), one (ACs1) in sanga, and another (ACz1) in zebu cattle, fitting well with expectations. ACt1 included genes primarily associated with innate immunity; ACt2 involved bitter taste receptor genes of importance to adaptation to changing environments; ACt3 included 68 genes coding ATP-binding proteins, some of them located on trypanotolerance-related QTL regions, that can partially underlie immune response and the additive mechanism of trypanotolerance; ACt4 was associated with growth and small size (NPPC gene); ACs1 included genes involved in immune response; and ACz1 is related with ectoparasite resistance. Our results provide a new set of genomic areas and candidate genes giving new insights on the genomic impact of adaptation in West African cattle.

Population Structure Assessed Using Microsatellite and SNP Data: An Empirical Comparison in West African Cattle

Simple Summary

Projection of genetic variability on geographic maps is a useful strategy to ascertain population structure and gene flow events when previous genetic information on the scenarios analyzed is not high. Here, we compared the performance of microsatellite sets and Single Nucleotide Polymorphism (SNP) arrays to identify the population structure and between-populations identity in a sample of West African cattle. Large SNP arrays were superior in detecting the population structure due to a more precise assessment of genotypic information of the individuals. However, the projection of genetic parameters on geographical maps was comparable between the SNP and microsatellite data. Geographic-based analyses of genetic variation areuseful inavoiding overinterpretation of the results obtained. Microsatellite markers can still be useful, particularly if the research focuses on non-model organisms or if either the funding or the availability of efficient hardware and software to handle large datasets is limited.

Abstract

A sample of 185 West African cattle belonging to nine different taurine, sanga, and zebu populations was typed using a set of 33 microsatellites and the BovineHD BeadChip of Illumina. The information provided by each type of marker was summarized via clustering methods and principal component analyses (PCA). The aim was to assess differences in performance between both marker types for the identification of population structure and the projection of genetic variability on geographical maps. In general, both microsatellites and Single Nucleotide Polymorphism (SNP) allowed us to differentiate taurine cattle from zebu and sanga cattle, which, in turn, would form a single population. Pearson and Spearman correlation coefficients computed among the admixture coefficients (fitting K = 2) and the eigenvectors corresponding to the first two factors identified using PCA on both microsatellite and SNP data were statistically significant (most of them having p < 0.0001) and high. However, SNP data allowed for a better fine-scale identification of population structure within taurine cattle: Lagunaire cattle from Benin were separated from two different N’Dama cattle samples. Furthermore, when clustering analyses assumed the existence of two parental populations only (K = 2), the SNPs could differentiate a different genetic background in Lagunaire and N’Dama cattle. Although the two N’Dama cattle populations had very different breeding histories, the microsatellite set could not separate the two N’Dama cattle populations. Classic bidimensional dispersion plots constructed using factors identified via PCA gave different shapes for microsatellites and SNPs: plots constructed using microsatellite polymorphism would suggest the existence of weakly differentiated, highly intermingled, subpopulations. However, the projection of the factors identified on synthetic maps gave comparable images. This would suggest that results on population structuring must be interpreted with caution. The geographic projection of genetic variation on synthetic maps avoids interpretations that go beyond the results obtained, particularly when previous information on the analyzed populations is scant. Factors influencing the performance of the projection of genetic parameters on geographic maps, together with restrictions that may affect the election of a given type of markers, are discussed.

Ancient Homozygosity Segments in West African Djallonké Sheep Inform on the Genomic Impact of Livestock Adaptation to the Environment

A sample of Burkina Faso Djallonké (West African Dwarf) sheep was analyzed to identify stretches of homozygous segments (runs of homozygosity; ROH) overlapping with ancient homozygosity-by-descent (HBD) segments. HBD segments were considered ancient if they were likely to be inherited from ancestors living from 1024 to 2048 generations ago, roughly coinciding with the time in which sheep entered into West Africa. It is hypothesized that such homozygous segments can inform on the effect of the sheep genome of human-mediated selection for adaptation to this harsh environment. PLINK analyses allowed to identify a total of 510 ROH segments in 127 different individuals that could be summarized into 124 different ROH. A total of 32,968 HBD segments were identified on 119 individuals using the software ZooRoH. HBD segments inherited from ancestors living 1024 and 2048 generations ago were identified on 61 individuals. The overlap between consensus ROH identified using PLINK and HBD fragments putatively assigned to generations 1024 and 2048 gave 108 genomic areas located on 17 different ovine chromosomes which were considered candidate regions for gene-annotation enrichment analyses. Functional annotation allowed to identify six statistically significant functional clusters involving 50 candidate genes. Cluster 1 was involved in homeostasis and coagulation; functional clusters 2, 3, and 6 were associated to innate immunity, defense against infections, and white blood cells proliferation and migration, respectively; cluster 4 was involved in parasite resistance; and functional cluster 5, formed by 20 genes, was involved in response to stress. The current analysis confirms the importance of genomic areas associated to immunity, disease resistance, and response to stress for adaptation of sheep to the challenging environment of humid Sub-Saharan West Africa.

Genomic scan of selective sweeps in Djallonké (West African Dwarf) sheep shed light on adaptation to harsh environments

The Djallonké (West African Dwarf) sheep is a small-sized haired sheep resulting from a costly evolutionary process of natural adaptation to the harsh environment of West Africa including trypanosome challenge. However, genomic studies carried out in this sheep are scant. In this research, genomic data of 184 Djallonké sheep (and 12 Burkina-Sahel sheep as an outgroup) generated using medium-density SNP Chips were analyzed. Three different statistics (iHS, XP-EHH and nSL) were applied to identify candidate selection sweep regions spanning genes putatively associated with adaptation of sheep to the West African environment. A total of 207 candidate selection sweep regions were defined. Gene-annotation enrichment and functional annotation analyses allowed to identify three statistically significant functional clusters involving 12 candidate genes. Genes included in Functional Clusters associated to selection signatures were mainly related to metabolic response to stress, including regulation of oxidative and metabolic stress and thermotolerance. The bovine chromosomal areas carrying QTLs for cattle trypanotolerance were compared with the regions on which the orthologous functional candidate cattle genes were located. The importance of cattle BTA4 for trypanotolerant response might have been conserved between species. The current research provides new insights on the genomic basis for adaptation and highlights the importance of obtaining information from non-cosmopolite livestock populations managed in harsh environments.

Confronting pastoralists' knowledge of cattle breeds raised in the extensive production systems of Benin with multivariate analyses of morphological traits

Cross-border seasonal livestock movements in West Africa bring into close contact several cattle breeds. In the coastal countries hosting migrating herders from the Sahel, it often affects the genetic variability and geographical distribution of traditional cattle breeds, through their indiscriminate but also intended crossbreeding with larger-framed Sahelian cattle breeds. The need to secure and effectively manage this genetic variability, in order to respond to changing production and market conditions, is widely recognized by the scientific community, livestock herders and policy-makers. This however requires a comprehensive knowledge of the breeds’ characteristics. The indigenous criteria used by pastoralists to characterize and distinguish cattle breeds remain unclear and further validation is required. This study was therefore designed to document and validate herders' knowledge on cattle breeds. From June 2015 to June 2016, 803 cattle herders participated in a phenotypic breed description in seven pastoral communities across the country. Each cattle herder was asked to name and describe morphologically the different cattle breeds in his herd. Subsequently, fifteen body measurements taken on a total of 1401 adult cattle (964 cows and 439 bulls) were submitted to multivariate analyses. Participants distinguished ten different cattle breeds kept in traditional herds according to six primary morphological traits and clearly separated zebuine from taurine breeds. These results were consistent with those of the multivariate analyses of the measured traits. However, herders’ classification approach proved to be more accurate in distinguishing breeds within the zebuine subspecies. Hence, while metric measurements and molecular genetic analyses are promising approaches to fill the knowledge gap on the diversity of local farm animal genetic resources, they should integrate livestock herders’ traditional knowledge for more precision.

Timing and Extent of Inbreeding in African Goats

Genetic characterization of African goats is one of the current priorities in the improvement of goats in the continent. This study contributes to the characterization effort by determining the levels and number of generations to common ancestors ("age") associated with inbreeding in African goat breeds and identifies regions that contain copy number variation mistyped as being homozygous. Illumina 50k single nucleotide polymorphism genotype data for 608 goats from 31 breeds were used to compute the level and age of inbreeding at both local (marker) and global levels (F_G) using a model-based approach based on a hidden Markov model. Runs of homozygosity (ROH) segments detected using the Viterbi algorithm led to ROH-based inbreeding coefficients for all ROH (F_ROH) and for ROH longer than 2 Mb (F_{ROH > 2Mb}). Some of the genomic regions identified as having ROH are likely to be hemizygous regions (copy number deletions) mistyped as homozygous regions. Although the proportion of these miscalled ROH is small and does not substantially affect estimates of levels of inbreeding for individual animals, the inbreeding metrics were adjusted by removing these regions from the ROH. All the inbreeding metrics varied widely across breeds, with overall means of 0.0408, 0.0370, and 0.0691 and medians of 0.0125, 0.0098, and 0.0366 for F_ROH, F_{ROH > 2Mb}, and F_G, respectively. Several breeds (including Menabe and Sofia from Madagascar) had high proportions of recent inbreeding, while Small East African, Ethiopian, and most of the West African breeds (including West African Dwarf) had more ancient inbreeding.

Leveraging Available Resources and Stakeholder Involvement for Improved Productivity of African Livestock in the Era of Genomic Breeding

The African continent is home to diverse populations of livestock breeds adapted to harsh environmental conditions with more than 70% under traditional systems of management. Animal productivity is less than optimal in most cases and is faced with numerous challenges including limited access to adequate nutrition and disease management, poor institutional capacities and lack of adequate government policies and funding to develop the livestock sector. Africa is home to about 1.3 billion people and with increasing demand for animal proteins by an ever growing human population, the current state of livestock productivity creates a significant yield gap for animal products. Although a greater section of the population, especially those living in rural areas depend largely on livestock for their livelihoods; the potential of the sector remains underutilized and therefore unable to contribute significantly to economic development and social wellbeing of the people. With current advances in livestock management practices, breeding technologies and health management, and with inclusion of all stakeholders, African livestock populations can be sustainably developed to close the animal protein gap that exists in the continent. In particular, advances in gene technologies, and application of genomic breeding in many Western countries has resulted in tremendous gains in traits like milk production with the potential that, implementation of genomic selection and other improved practices (nutrition, healthcare, etc.) can lead to rapid improvement in traits of economic importance in African livestock populations. The African livestock populations in the context of this review are limited to cattle, goat, pig, poultry, and sheep, which are mainly exploited for meat, milk, and eggs. This review examines the current state of livestock productivity in Africa, the main challenges faced by the sector, the role of various stakeholders and discusses in-depth strategies that can enable the application of genomic technologies for rapid improvement of livestock traits of economic importance.

Human and ecological determinants of the spatial structure of local breed diversity

Since domestication, a large number of livestock breeds adapted to local conditions have been created by natural and artificial selection, representing one of the most powerful ways in which human groups have constructed niches to meet their need. Although many authors have described local breeds as the result of culturally and environmentally mediated processes, this study, located in mainland Spain, is the first aimed at identifying and quantifying the environmental and human contributions to the spatial structure of local breed diversity, which we refer to as livestock niche. We found that the more similar two provinces were in terms of human population, ecological characteristics, historical ties, and geographic distance, the more similar the composition of local breeds in their territories. Isolation by human population distance showed the strongest effect, followed by isolation by the environment, thus supporting the view of livestock niche as a socio-cultural product adapted to the local environment, in whose construction humans make good use of their ecological and cultural inheritances. These findings provide a useful framework to understand and to envisage the effects of climate change and globalization on local breeds and their livestock niches.

Morphological assessment of the Zebu Bororo (Wodaabé) cattle of Niger in the West African zebu framework

Abstract. A total of 357 adult cows and 29 sires belonging to the long-horned Niger Zebu Bororo cattle population were assessed for 13 body measurements and 11 qualitative traits. Data were jointly analysed with 311 cows and 64 sires belonging to other four West African zebu cattle populations, sampled in Burkina Faso and Benin, representative of both the short-horned and the long-horned West African zebu groups using multivariate statistical methods. Besides the other long-horned zebu breed analysed (Zebu Mbororo of Burkina Faso), Zebu Bororo cattle tended to have the highest mean values for all body measurements. Mahalanobis distance matrices further informed that pairs involving Zebu Bororo cattle had the higher differentiation of the dataset. However, contour plots constructed using eigenvalues computed via principal component analysis (PCA) illustrated a lack of differentiation among West African zebu cattle populations at the body measurements level. Correspondence analysis carried out on the 11 qualitative traits recorded allowed for ascertaining a clear differentiation between the Zebu Bororo and the other zebu cattle populations analysed which, in turn, did not show a clear differentiation at the qualitative type traits level. In our data, Zebu Bororo cattle had in high frequency qualitative features such as dropped ears, lyre-shaped horns and red-pied coat colour that are not frequently present in the other West African zebu populations analysed. A directional selection due to a rough consensus of the stock-keepers may be hypothesised. Performance of further analyses to assess the degree in which such breeding differences may be related to genetic or production differences are advised.

Differences in genetic structure assessed using Y-chromosome and mitochondrial DNA markers do not shape the contributions to diversity in African sires

Up to 173 African sires belonging to 11 different subpopulations representative of four cattle groups were analysed for six Y-specific microsatellite loci and a mitochondrial DNA fragment. Differences in Y-chromosome and mtDNA haplotype structuring were assessed. In addition, the effect of such structuring on contributions to total genetic diversity was assessed. Thirty-five Y-chromosome and 71 mtDNA haplotypes were identified. Most Y-chromosomes analysed (73.4%) were of zebu origin (11 haplotypes). Twenty-two Y-haplotypes (44 samples) belonged to the African taurine subfamily Y2a. All mtDNA haplotypes belonged to the "African" taurine T1 haplogroup with 16 samples and nine haplotypes belonging to a recently identified subhaplogroup (T1e). Median-joining networks showed that Y-chromosome phylogenies were highly reticulated with clear separation between zebu and taurine clusters. Mitochondrial haplotypes showed a clear star-like shape with small number of mutations separating haplotypes. Mitochondrial-based F_ST -statistics computed between cattle groups tended to be statistically non-significant (p > .05). Most F_ST values computed among groups and subpopulations using Y-chromosome markers were statistically significant. AMOVA confirmed that divergence between cattle groups was only significant for Y-chromosome markers (Φ_CT  = 0.209). At the mitochondrial level, African sires resembled an undifferentiated population with individuals explaining 94.3% of the total variance. Whatever the markers considered, the highest contributions to total Nei's gene diversity and allelic richness were found in West African cattle. Genetic structuring had no effect on patterns of contributions to diversity.

Detecting the existence of gene flow between Spanish and North African goats through a coalescent approach

Human-driven migrations are one of the main processes shaping the genetic diversity and population structure of domestic species. However, their magnitude and direction have been rarely analysed in a statistical framework. We aimed to estimate the impact of migration on the population structure of Spanish and African goats. To achieve this goal, we analysed a dataset of 1,472 individuals typed with 23 microsatellites. Population structure of African and Spanish goats was moderate (mean FST = 0.07), with the exception of the Canarian and South African breeds that displayed a significant differentiation when compared to goats from North Africa and Nigeria. Measurement of gene flow with Migrate-n and IMa coalescent genealogy samplers supported the existence of a bidirectional gene flow between African and Spanish goats. Moreover, IMa estimates of the effective number of migrants were remarkably lower than those calculated with Migrate-n and classical approaches. Such discrepancies suggest that recent divergence, rather than extensive gene flow, is the main cause of the weak population structure observed in caprine breeds.

Multivariate characterization of morphological traits in West African cattle sires

Abstract. A total of 183 adult sires belonging to nine West African cattle breeds sampled in 67 villages of Mali, Burkina Faso and Benin were assessed for 16 body measurements and 18 qualitative traits. Within type of cattle (zebu, sanga or taurine), the different breeds analysed showed large differences in body measurements. In general, taurine breeds had lower average values than the zebu breeds while sanga cattle tended to have intermediate values. Principal component analysis identified three factors characterising body measurements. Factor 1 summarised the information provided by those traits characterising the size of the individuals and explained 59.0 % of the variability. Factor 2 tended to gather information characterising the body width and explained 8.0 % of the variation. Less representative, Factor 3 (6.6 % of the variability) had no clear interpretation. Qualitative traits did not allow to distinguish among either cattle groups or breeds. Two Correspondence Analysis Dimensions computed on qualitative traits (explaining 26.2 and 15.5 % of the variability, respectively) did not allow to differentiate between zebu, sanga or taurine cattle breeds. Our results confirm that, in the framework of a general appearance, body measurements are the main criteria for differentiating West African cattle breeds. Furthermore, the current research has not allowed to identify breeding preferences on qualitative type traits in West African cattle sires. Therefore, homogenisation of the appearance of individuals within cattle breed is not expected.

Geographical assessment of body measurements and qualitative traits in West African cattle

A total of 1015 adult cows belonging to nine West African cattle breeds were assessed for 16 body measurements and 18 qualitative traits to ascertain the existence of geographical patterns of variation. Sampling was carried out in 29 different provinces of Mali, Burkina Faso and Benin. For body measurements, taurine breeds took lower average values than the zebu breeds. Sanga cattle took intermediate values. Qualitative traits did not allow to differentiate among cattle groups (taurine, zebu or sanga) or breeds. Principal component analysis identified two factors explaining 56.4 and 9.2 % of the variance for body measurements, respectively. Two correspondence analysis dimensions computed on qualitative traits explained a small proportion of the variability (20.8 and 13.5 %, respectively). Contour plots were constructed using the eigenvalues computed for each individual and either factor or dimension identified; confidence regions calculated confirmed that body measurements clearly differentiated zebu and taurine cattle breeds while qualitative traits did not. Factor 1 was projected on a geographical map, using provinces as nodes, to assess breed-free variation for body measurements. A pattern of continuous variation from the Sahel area southwards was identified. Probably, breeding decisions promoting the crosses between zebu-like and taurine cattle are underlying this geographical pattern of variation. The implementation of selection strategies aiming at the increase of the productivity of native West African taurine cattle breeds while avoiding looses in trypanotolerant ability would be highly advisable.

Assessing introgression of Sahelian zebu genes into native Bos taurus breeds in Burkina Faso

A total of 350 samples were analyzed to estimate zebu gene proportions into two different taurine cattle breeds of Burkina Faso (Lobi and N'Dama) using 38 microsatellites and various statistical methodologies. West African and East African zebu samples were sequentially used as reference parental populations. Furthermore, N'Dama cattle from Congo, the composite South African Bonsmara cattle breed and a pool of European cattle were used successively as second parental populations. Independently of the methodology applied: (a) the use of West African zebu samples gave higher admixture coefficients than the East African zebu; (b) the higher zebu proportions were estimated when the European cattle was used as parental population 2; and (c) the use of the N'Dama population from Congo as parental population 2 gave the more consistent zebu proportion estimates for both the Lobi and the N'Dama breeds. In any case, the zebu admixture proportions estimated were not negligible and were always higher in the N'Dama cattle than in the Lobi cattle of Burkina Faso. This suggested that the introgression of Sahelian zebu genes into the taurine cattle of Southern West Africa can follow a complex pattern that can depend on local agro-ecological features. The current research pointed out that the estimation of admixture coefficients is highly dependent on both the assumptions underlying the methodologies applied and the selection of parental populations. Our analyses suggest that either too high or nil genetic identity between the parental and the expectedly derived populations must be avoided.