Detection of selection signatures within candidate regions underlying trypanotolerance in outbred cattle populations

Assessment of Genomic Diversity and Selective Pressures in Crossbred Dairy Cattle of Pakistan

Improving the low productivity levels of native cattle breeds in smallholder farming systems is a pressing concern in Pakistan. Crossbreeding high milk-yielding holstein friesian (HF) breed with the adaptability and heat tolerance of Sahiwal cattle has resulted in offspring that are well-suited to local conditions and exhibit improved milk yield. The exploration of how desirable traits in crossbred dairy cattle are selected has not yet been investigated. This study aims to provide the first overview of the selective pressures on the genome of crossbred dairy cattle in Pakistan. A total of eighty-one crossbred, thirty-two HF and twenty-four Sahiwal cattle were genotyped, and additional SNP genotype data for HF and Sahiwal were collected from a public database to equate the sample size in each group. Within-breed selection signatures in crossbreds were investigated using the integrated haplotype score. Crossbreds were also compared to each of their parental breeds to discover between-population signatures of selection using two approaches: cross-population extended haplotype homozygosity and fixation index. We identified several overlapping genes associated with production, immunity, and adaptation traits, including U6, TMEM41B, B4GALT7, 5S_rRNA, RBM27, POU4F3, NSD1, PRELID1, RGS14, SLC34A1, TMED9, B4GALT7, OR2AK3, OR2T16, OR2T60, OR2L3, and CTNNA1. Our results suggest that regions responsible for milk traits have generally experienced stronger selective pressure than others.

Genetic differentiation and structuration of the Gobra zebu cattle breeds reared in Senegal

The Gobra zebu genetic breeding program has resulted in the genetic improvement of a new population. This population gained genetic characteristics that set them apart from the other cattle populations reared in Senegal. The cause of these differences might be the reproductive isolation and selection to which this population of the "Centre de Recherches Zootechniques" of Dahra has been subjected since the 1950s. This study aimed to assess the genetic differentiation and structuration of this population compared to the main cattle breeds used in Senegal. A total of 180 individuals, selected from the Gobra selection nucleus and bovine populations from four main breeds in Senegal, were included in this study. We used a panel of 21 microsatellite markers among those recommended by the Food Agriculture Organization, to conduct the molecular genotyping of our sampled populations. The basic genetic parameters of differentiation and structuration were calculated using various bioinformatics software. The results of this study, particularly the degrees of genetic differentiation (Fst), the coefficient of genetic homogeneity (Gst), and the gene flow (Nm), show a significant genetic differentiation of the Gobra from the station compared to the other populations studied. Structuring and phylogeny analyses reveal a micro-structuring within the Gobra population as a novelty. This micro-structuring clearly identifies the Gobra individuals from Dahra's station among the other Gobra populations studied. The main causes of these observations would be reproductive isolation and the selection pressure exerted on this population for several decades.

Genetic diversity of bovine populations raised in Senegal

Background

The Gobra zebu and N'dama taurine cattle breeds are important genetic animal resources for Senegal. For several decades, genetic breeding programmes have been devoted to them at the Centre de Recherches Zootechniques de Dahra and Kolda. Since then, these animals have been subjected to mass selection, mainly in closed selection nuclei.

Objective

This study aims to assess the genetic diversity within these selection nuclei in order to orient future selection strategies.

Material and methods

The study was carried out on the Gobra zebu and N'dama taurine populations from selection nuclei of Dahra and Kolda respectively, which were compared to 5 other populations of the main cattle breeds in Senegal. One hundred eighty (180) animals were genotyped with 21 microsatellite markers recommended by the Food and Agriculture Organisation.

Results

All populations were found to be polymorphic with a PIC of over 55%. However, animals from the CRZ‐Dahra (indigenous) and CRZ‐Kolda stations had the lowest mean heterozygosity (0.643 and 0.591 respectively). The other populations had an average heterozygosity between 0.650 and 0.737.

Conclusion

The cattle populations maintained at the different CRZs show a lower genetic diversity than the other populations described in our study. The main reasons for this are reproductive isolation and selection pressure on these populations.

Identification and Characterization of Copy Number Variations Regions in West African Taurine Cattle

A total of 106 West African taurine cattle belonging to the Lagunaire breed of Benin (33), the N’Dama population of Burkina Faso (48), and N’Dama cattle sampled in Congo (25) were analyzed for Copy Number Variations (CNVs) using the BovineHDBeadChip of Illumina and two different CNV calling programs: PennCNV and QuantiSNP. Furthermore, 89 West African zebu samples (Bororo cattle of Mali and Zebu Peul sampled in Benin and Burkina Faso) were used as an outgroup to ensure that analyses reflect the taurine cattle genomic background. Analyses identified 307 taurine-specific CNV regions (CNVRs), covering about 56 Mb on all bovine autosomes. Gene annotation enrichment analysis identified a total of 840 candidate genes on 168 taurine-specific CNVRs. Three different statistically significant functional term annotation clusters (from ACt1 to ACt3) involved in the immune function were identified: ACt1 includes genes encoding lipocalins, proteins involved in the modulation of immune response and allergy; ACt2 includes genes encoding coding B-box-type zinc finger proteins and butyrophilins, involved in innate immune processes; and Act3 includes genes encoding lectin receptors, involved in the inflammatory responses to pathogens and B- and T-cell differentiation. The overlap between taurine-specific CNVRs and QTL regions associated with trypanotolerant response and tick-resistance was relatively low, suggesting that the mechanisms underlying such traits may not be determined by CNV alterations. However, four taurine-specific CNVRs overlapped with QTL regions associated with both traits on BTA23, therefore suggesting that CNV alterations in major histocompatibility complex (MHC) genes can partially explain the existence of genetic mechanisms shared between trypanotolerance and tick resistance in cattle. This research contributes to the understanding of the genomic features of West African taurine cattle.

Local Ancestry to Identify Selection in Response to Trypanosome Infection in Baoulé x Zebu Crossbred Cattle in Burkina Faso

The genomes of crossbred (admixed) individuals are a mosaic of ancestral haplotypes formed by recombination in each generation. The proportion of these ancestral haplotypes in certain genomic regions can be responsible for either susceptibility or tolerance against pathogens, and for performances in production traits. Using a medium-density genomic marker panel from the Illumina Bovine SNP50 BeadChip, we estimated individual admixture proportions for Baoulé x Zebu crossbred cattle in Burkina Faso, which were tested for trypanosome infection by direct ELISA from blood samples. Furthermore, we calculated local ancestry deviation from average for each SNP across 29 autosomes to identify potential regions under selection in the trypanotolerant Baoulé cattle and their crossbreds. We identified significant deviation from the local average ancestry (above 5 and 10% genome-wide thresholds) on chromosomes 8 and 19 in the positive animals, while the negative ones showed higher deviation on chromosomes 6, 19, 21, and 22. Some candidate genes on chromosome 6 (PDGFRA) and chromosome 19 (CDC6) have been found associated to trypanotolerance in West African taurines. Screening for F_ST outliers in trypanosome positive/negative animals we detected seven variants putatively under selection. Finally, we identified a minimum set of highly ancestry informative markers for routine admixture testing. The results of this study contribute to a better understanding of the genetic basis of trypanotolerance in Baoulé cattle and their crossbreeds. Furthermore, we provide a small informative marker set to monitor admixture in this valuable indigenous breed. As such, our results are important for conserving the genetic uniqueness and trypanotolerance of Baoulé cattle, as well as for the improvement of Baoulé and Zebu crossbreds in specific community-based breeding programs.

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.

A multi-breed GWAS for morphometric traits in four Beninese indigenous cattle breeds reveals loci associated with conformation, carcass and adaptive traits

BACKGROUND

Specific adaptive features including disease resistance and growth abilities in harsh environments are attributed to indigenous cattle breeds of Benin, but these breeds are endangered due to crossbreeding. So far, there is a lack of systematic trait recording, being the basis for breed characterizations, and for structured breeding program designs aiming on conservation. Bridging this gap, own phenotyping for morphological traits considered measurements for height at withers (HAW), sacrum height (SH), heart girth (HG), hip width (HW), body length (BL) and ear length (EL), including 449 cattle from the four indigenous Benin breeds Lagune, Somba, Borgou and Pabli. In order to utilize recent genomic tools for breed characterizations and genetic evaluations, phenotypes for novel traits were merged with high-density SNP marker data. Multi-breed genetic parameter estimations and genome-wide association studies (GWAS) for the six morphometric traits were carried out. Continuatively, we aimed on inferring genomic regions and functional loci potentially associated with conformation, carcass and adaptive traits.

RESULTS

SNP-based heritability estimates for the morphometric traits ranged between 0.46 ± 0.14 (HG) and 0.74 ± 0.13 (HW). Phenotypic and genetic correlations ranged from 0.25 ± 0.05 (HW-BL) to 0.89 ± 0.01 (HAW-SH), and from 0.14 ± 0.10 (HW-BL) to 0.85 ± 0.02 (HAW-SH), respectively. Three genome-wide and 25 chromosome-wide significant SNP positioned on different chromosomes were detected, located in very close chromosomal distance (±25 kb) to 15 genes (or located within the genes). The genes PIK3R6 and PIK3R1 showed direct functional associations with height and body size. We inferred the potential candidate genes VEPH1, CNTNAP5, GYPC for conformation, growth and carcass traits including body weight and body fat deposition. According to their functional annotations, detected potential candidate genes were associated with stress or immune response (genes PTAFR, PBRM1, ADAMTS12) and with feed efficiency (genes MEGF11 SLC16A4, CCDC117).

CONCLUSIONS

Accurate measurements contributed to large SNP heritabilities for some morphological traits, even for a small mixed-breed sample size. Multi-breed GWAS detected different loci associated with conformation or carcass traits. The identified potential candidate genes for immune response or feed efficiency indicators reflect the evolutionary development and adaptability features of the breeds.

The role of agro-ecological factors and transboundary transhumance in shaping the genetic diversity in four indigenous cattle populations of Benin

The indigenous cattle population of Benin is a diverse mix of taurine and hybrid breeds shaped by diverse ecological and climatic conditions with eight agro-ecological zones (AEZ). Presumably, the taurine breeds face current endangerment due to ongoing indicine introgression following climate change and transboundary transhumance. The aim of the study was to investigate the genetic diversity and population structure of the indigenous breeds Lagune, Somba, Pabli and Borgou considering spatial agro-ecological and socio-economic factors (transhumance) based on 50k SNP and microsatellite data. Among the four sampled breeds, six genetic clusters were identified using model-free (discriminant analysis of principal components) and model-based (TESS and ADMIXTURE) methods separating taurine from hybrid breeds. Results based on an extension with publicly available historic SNP data sets from taurine and indicine West African cattle and additional outgroups provided additional insight into changes of genetic structure in the sampled breeds over time. Both taurine breeds, Somba and Lagune, showed a stable foundation but also spatially limited partial indicine introgression associated with transhumance leading to high genetic diversity. In addition, we found evidence for spatial diversity and changes in genetic structure over time in the Borgou breed in comparison of our samples with the historic samples which could be explained by potential continuous indicine introgression into the Borgou breed in two sample regions. Results for the Pabli breed do not conclusively point to full absorbance by the Borgou in comparison with all available Borgou samples. Further research is needed in this regard.

Genome-wide detection of signatures of selection in three Valdostana cattle populations

The Valdostana is a local dual purpose cattle breed developed in Italy. Three populations are recognized within this breed, based on coat colour, production level, morphology and temperament: Valdostana Red Pied (VPR), Valdostana Black Pied (VPN) and Valdostana Chestnut (VCA). Here, we investigated putative genomic regions under selection among these three populations using the Bovine 50K SNP array by combining three different statistical methods based either on allele frequencies (F_ST ) or extended haplotype homozygosity (iHS and Rsb). In total, 8, 5 and 8 chromosomes harbouring 13, 13 and 16 genomic regions potentially under selection were identified by at least two approaches in VPR, VPN and VCA, respectively. Most of these candidate regions were population-specific but we found one common genomic region spanning 2.38 Mb on BTA06 which either overlaps or is located close to runs of homozygosity islands detected in the three populations. This region included inter alia two well-known genes: KDR, a well-established coat colour gene, and CLOCK, which plays a central role in positive regulation of inflammatory response and in the regulation of the mammalian circadian rhythm. The other candidate regions identified harboured genes associated mainly with milk and meat traits as well as genes involved in immune response/inflammation or associated with behavioural traits. This last category of genes was mainly identified in VCA, which is selected for fighting ability. Overall, our results provide, for the first time, a glimpse into regions of the genome targeted by selection in Valdostana cattle. Finally, this study illustrates the relevance of using multiple complementary approaches to identify genomic regions putatively under selection in livestock.

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.

Analysis of pooled genome sequences from Djallonke and Sahelian sheep of Ghana reveals co-localisation of regions of reduced heterozygosity with candidate genes for disease resistance and adaptation to a tropical environment

BACKGROUND

The Djallonke sheep is well adapted to harsh environmental conditions, and is relatively resistant to Haemonchosis and resilient to animal trypanosomiasis. The larger Sahelian sheep, which cohabit the same region, is less well adapted to these disease challenges. Haemonchosis and Trypanosomiasis collectively cost the worldwide animal industry billions of dollars in production losses annually.

RESULTS

Here, we separately sequenced and then pooled according to breed the genomes from five unrelated individuals from each of the Djallonke and Sahelian sheep breeds (sourced from Ghana), at greater than 22-fold combined coverage for each breed. A total of approximately 404 million (97%) and 343 million (97%) sequence reads from the Djallonke and Sahelian breeds respectively, were successfully mapped to the sheep reference genome Oar v3.1. We identified approximately 11.1 million and 10.9 million single nucleotide polymorphisms (SNPs) in the Djallonke and Sahelian breeds, with approximately 15 and 16% respectively of these not previously reported in sheep. Multiple regions of reduced heterozygosity were also found; 70 co-localised within genomic regions harbouring genes that mediate disease resistance, immune response and adaptation in sheep or cattle. Thirty- three of the regions of reduced heterozygosity co-localised with previously reported genes for resistance to haemonchosis and trypanosomiasis.

CONCLUSIONS

Our analyses suggest that these regions of reduced heterozygosity may be signatures of selection for these economically important diseases.

Haematology of N'Dama and West African Shorthorn cattle herds under natural Trypanosoma vivax challenge in Ghana

Background: Animal trypanosomosis is a major cause of economic loss in livestock production in Africa. A suggested control measure is to use breeds with traits of trypanotolerance. The study examines the effect of natural Trypanosoma vivax challenge on haematological parameters in two trypanotolerant cattle [N’Dama and West African Shorthorn (WASH)] herds.

Methods: Trypanosoma vivax-specific primers were used to diagnose T. vivax infection in an N’Dama herd at Cape Coast in southern Ghana and a WASH herd at Chegbani in northern Ghana from May to July 2011 in a cross-sectional study. Levels of haematological parameters comprising packed cell volume (PCV), haemoglobin (Hb) concentration and red blood cell (RBC) and total white blood cell (TWBC) counts; differential WBC counts (neutrophils, lymphocytes, eosinophils, monocytes and basophils); and RBC indices of mean corpuscular volume (MCV), mean corpuscular haemoglobin (MCH) and mean corpuscular haemoglobin concentration (MCHC) were determined in blood samples and then compared between infected and uninfected cattle.

Results: We found that haematological indices for infected and uninfected animals in both breeds were within the normal range. However, the mean PCV values for T. vivax-infected WASH and N’Dama were lower in infected compared to uninfected animals. The difference was significant ( p< 0.05) in N’Dama but not in WASH.

Conclusion: Despite the presence of infection by T. vivax, N’Dama and WASH cattle maintained their haematological parameters within acceptable normal ranges,  which confirms their trypanotolerant trait. This highlights the need for low-input traditional African farmers in medium, high and severe tsetse challenge areas to be educated on the advantages of N’Dama and WASH breeds to increase their utilization in integrated tsetse and trypanosomosis control programmes.

Combatting African Animal Trypanosomiasis (AAT) in livestock: The potential role of trypanotolerance

African Animal Trypanosomiasis (AAT) is endemic in at least 37 of the 54 countries in Africa. It is estimated to cause direct and indirect losses to the livestock production industry in excess of US$ 4.5 billion per annum. A century of intervention has yielded limited success, owing largely to the extraordinary complexity of the host-parasite interaction. Trypanotolerance, which refers to the inherent ability of some African livestock breeds, notably Djallonke sheep, N'Dama cattle and West African Dwarf goats, to withstand a trypanosomiasis challenge and still remain productive without any form of therapy, is an economically sustainable option for combatting this disease. Yet trypanotolerance has not been adequately exploited in the fight against AAT. In this review, we describe new insights into the genetic basis of trypanotolerance and discuss the potential of exploring this phenomenon as an integral part of the solution for AAT, particularly, in the context of African animal production systems.

A Meta-Assembly of Selection Signatures in Cattle

Since domestication, significant genetic improvement has been achieved for many traits of commercial importance in cattle, including adaptation, appearance and production. In response to such intense selection pressures, the bovine genome has undergone changes at the underlying regions of functional genetic variants, which are termed “selection signatures”. This article reviews 64 recent (2009–2015) investigations testing genomic diversity for departure from neutrality in worldwide cattle populations. In particular, we constructed a meta-assembly of 16,158 selection signatures for individual breeds and their archetype groups (European, African, Zebu and composite) from 56 genome-wide scans representing 70,743 animals of 90 pure and crossbred cattle breeds. Meta-selection-scores (MSS) were computed by combining published results at every given locus, within a sliding window span. MSS were adjusted for common samples across studies and were weighted for significance thresholds across and within studies. Published selection signatures show extensive coverage across the bovine genome, however, the meta-assembly provides a consensus profile of 263 genomic regions of which 141 were unique (113 were breed-specific) and 122 were shared across cattle archetypes. The most prominent peaks of MSS represent regions under selection across multiple populations and harboured genes of known major effects (coat color, polledness and muscle hypertrophy) and genes known to influence polygenic traits (stature, adaptation, feed efficiency, immunity, behaviour, reproduction, beef and dairy production). As the first meta-assembly of selection signatures, it offers novel insights about the hotspots of selective sweeps in the bovine genome, and this method could equally be applied to other species.

Tolerance to Trypanosomatids: A Threat, or a Key for Disease Elimination?

So far, research on trypanosomatid infections has been driven by 'disease by disease' approaches, leading to different concepts and control strategies. It is, however, increasingly clear that they share common features such as the ability to generate long-lasting asymptomatic infections in their mammalian hosts. Trypanotolerance, long integrated in animal African trypanosomiasis control, historically refers to the ability of cattle breeds to limit Trypanosoma infection and pathology, but has only recently been recognized in humans. Whilst trypanotolerance is absent from the vocabulary on leishmaniasis and Chagas disease, asymptomatic infections also occur. We review the concept of trypanotolerance across the trypanosomatids and discuss the importance of asymptomatic carriage in the current context of elimination.

Genetic diversity and phylogenetic relationships in local cattle breeds of Senegal based on autosomal microsatellite markers

Aim:

In Senegal, uncontrolled cross-breeding of cattle breeds and changes in production systems are assumed to lead to an increase of gene flow between populations. This might constitute a relevant threat to livestock improvement. Therewith, this study was carried out to assess the current genetic diversity and the phylogenetic relationships of the four native Senegalese cattle breeds (Gobra zebu, Maure zebu, Djakoré, and N’Dama).

Methods:

Genomic DNA was isolated from blood samples of 120 unrelated animals collected from three agro-ecological areas of Senegal according to their phenotypic traits. Genotyping was done using 11 specific highly polymorphic microsatellite makers recommended by Food and Agriculture Organization. The basic measures of genetic variation and phylogenetic trees were computed using bioinformatics’ software.

Results:

A total of 115 alleles were identified with a number of alleles (Na) at one locus ranging from 6 to 16. All loci were polymorphic with a mean polymorphic information content of 0.76. The mean allelic richness (Rs) lay within the narrow range of 5.14 in N’Dama taurine to 6.10 in Gobra zebu. While, the expected heterozygosity (H_E) per breed was high in general with an overall mean of 0.76±0.04. Generally, the heterozygote deficiency (F_IS) of 0.073±0.026 was relatively due to inbreeding among these cattle breeds or the occurrence of population substructure. The high values of allelic and gene diversity showed that Senegalese native cattle breeds represented an important reservoir of genetic variation. The genetic distances and clustering trees concluded that the N’Dama cattle were most distinct among the investigated cattle populations. So, the principal component analyses showed qualitatively that there was an intensive genetic admixture between the Gobra zebu and Maure zebu breeds.

Conclusions:

The broad genetic diversity in Senegalese cattle breeds will allow for greater opportunities for improvement of productivity and adaptation relative to global changes. For the development of sustainable breeding and crossbreeding programs of Senegalese local breeds, effective management is needed towards genetic selection and transhumance to ensure their long-term survival.

Trypanosomosis: potential driver of selection in African cattle

Trypanosomosis is a serious cause of reduction in productivity of cattle in tsetse-fly infested areas. Baoule and other local Taurine cattle breeds in Burkina Faso are trypanotolerant. Zebuine cattle, which are also kept there are susceptible to trypanosomosis but bigger in body size. Farmers have continuously been intercrossing Baoule and Zebu animals to increase production and disease tolerance. The aim of this study was to compare levels of zebuine and taurine admixture in genomic regions potentially involved in trypanotolerance with background admixture of composites to identify differences in allelic frequencies of tolerant and non-tolerant animals. The study was conducted on 214 animals (90 Baoule, 90 Zebu, and 34 composites), genotyped with 25 microsatellites across the genome and with 155 SNPs in 23 candidate regions. Degrees of admixture of composites were analyzed for microsatellite and SNP data separately. Average Baoule admixture based on microsatellites across the genomes of the Baoule- Zebu composites was 0.31, which was smaller than the average Baoule admixture in the trypanosomosis candidate regions of 0.37 (P = 0.15). Fixation index F_ST measured in the overall genome based on microsatellites or with SNPs from candidate regions indicates strong differentiation between breeds. Nine out of 23 regions had F_ST ≥ 0.20 calculated from haplotypes or individual SNPs. The levels of admixture were significantly different from background admixture, as revealed by microsatellite data, for six out of the nine regions. Five out of the six regions showed an excess of Baoule ancestry. Information about best levels of breed composition would be useful for future breeding ctivities, aiming at trypanotolerant animals with higher productive capacity.

Adaptive admixture in the West African bovine hybrid zone: insight from the Borgou population

Understanding the adaptive response to environmental fluctuations represents a central issue in evolutionary biology. Population admixture between divergent ancestries has often been considered as an efficient short-term adaptation strategy. Cattle populations from the West African Bos taurus × Bos indicus hybrid zone represent a valuable resource to characterize the effect of such adaptive admixture at the genome level. We here provide a detailed assessment of the global and local genome ancestries of the Borgou breed, one of the most representative cattle of this hybrid zone. We analysed a large data set consisting of 38,100 SNPs genotyped on 203 Borgou and 591 individuals representative of all the different cattle ancestries. At the global genomic level, we show that Borgou is a stabilized admixed breed whose origin (c. 130 years ago) traces back to the great African rinderpest pandemic, several centuries after the last admixture events, the West African zebus originate from (c. 500 years ago). To identify footprints of adaptive admixture, we combined the identification of signatures of selection and the functional annotation of the underlying genes using systems biology tools. The detection of the SILV coat coloration gene likely under artificial selection may be viewed as a validation of our approach. Overall, our results suggest that the long-term presence of pathogens and the intermediate environmental conditions are the main acting selective pressures. Our analytical framework can be extended to other model or nonmodel species to understand the process that shapes the patterns of genetic variability in hybrid zones.

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.

Molecular ecology and selection in the drought-related Asr gene polymorphisms in wild and cultivated common bean (Phaseolus vulgaris L.)

BACKGROUND

The abscisic acid (ABA) pathway plays an important role in the plants' reaction to drought stress and ABA-stress response (Asr) genes are important in controlling this process. In this sense, we accessed nucleotide diversity at two candidate genes for drought tolerance (Asr1 and Asr2), involved in an ABA signaling pathway, in the reference collection of cultivated common bean (Phaseolus vulgaris L.) and a core collection of wild common bean accessions.

RESULTS

Our wild population samples covered a range of mesic (semi-arid) to very dry (desert) habitats, while our cultivated samples presented a wide spectrum of drought tolerance. Both genes showed very different patterns of nucleotide variation. Asr1 exhibited very low nucleotide diversity relative to the neutral reference loci that were previously surveyed in these populations. This suggests that strong purifying selection has been acting on this gene. In contrast, Asr2 exhibited higher levels of nucleotide diversity, which is indicative of adaptive selection. These patterns were more notable in wild beans than in cultivated common beans indicting that natural selection has played a role over long time periods compared to farmer selection since domestication.

CONCLUSIONS

Together these results suggested the importance of Asr1 in the context of drought tolerance, and constitute the first steps towards an association study between genetic polymorphism of this gene family and variation in drought tolerance traits. Furthermore, one of our major successes was to find that wild common bean is a reservoir of genetic variation and selection signatures at Asr genes, which may be useful for breeding drought tolerance in cultivated common bean.

A quasi-exclusive European ancestry in the Senepol tropical cattle breed highlights the importance of the slick locus in tropical adaptation

Background

The Senepol cattle breed (SEN) was created in the early XX^th century from a presumed cross between a European (EUT) breed (Red Poll) and a West African taurine (AFT) breed (N’Dama). Well adapted to tropical conditions, it is also believed trypanotolerant according to its putative AFT ancestry. However, such origins needed to be verified to define relevant husbandry practices and the genetic background underlying such adaptation needed to be characterized.

Methodology/Principal Findings

We genotyped 153 SEN individuals on 47,365 SNPs and combined the resulting data with those available on 18 other populations representative of EUT, AFT and Zebu (ZEB) cattle. We found on average 89% EUT, 10.4% ZEB and 0.6% AFT ancestries in the SEN genome. We further looked for footprints of recent selection using standard tests based on the extent of haplotype homozygosity. We underlined i) three footprints on chromosome (BTA) 01, two of which are within or close to the polled locus underlying the absence of horns and ii) one footprint on BTA20 within the slick hair coat locus, involved in thermotolerance. Annotation of these regions allowed us to propose three candidate genes to explain the observed signals (TIAM1, GRIK1 and RAI14).

Conclusions/Significance

Our results do not support the accepted concept about the AFT origin of SEN breed. Initial AFT ancestry (if any) might have been counter-selected in early generations due to breeding objectives oriented in particular toward meat production and hornless phenotype. Therefore, SEN animals are likely susceptible to African trypanosomes which questions the importation of SEN within the West African tsetse belt, as promoted by some breeding societies. Besides, our results revealed that SEN breed is predominantly a EUT breed well adapted to tropical conditions and confirmed the importance in thermotolerance of the slick locus.

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

Background

Introgression of Sahel livestock genes southwards in West Africa may be favoured by human activity and the increase of the duration of the dry seasons since the 1970’s. The aim of this study is to assess the gene flow patterns in Burkina Faso goat and to ascertain the most likely factors influencing geographic patterns of genetic variation in the Burkina Faso goat population.

Results

A total of 520 goat were sampled in 23 different locations of Burkina Faso and genotyped for a set of 19 microsatellites. Data deposited in the Dryad repository: http://dx.doi.org/10.5061/dryad.41h46j37. Although overall differentiation is poor (F_ST = 0.067 ± 0.003), the goat population of Burkina Faso is far from being homogeneous. Barrier analysis pointed out the existence of: a) genetic discontinuities in the Central and Southeast Burkina Faso; and b) genetic differences within the goat sampled in the Sahel or the Sudan areas of Burkina Faso. Principal component analysis and admixture proportion scores were computed for each population sampled and used to construct interpolation maps. Furthermore, Population Graph analysis revealed that the Sahel and the Sudan environmental areas of Burkina Faso were connected through a significant number of extended edges, which would be consistent with the hypothesis of long-distance dispersal. Genetic variation of Burkina Faso goat followed a geographic-related pattern. This pattern of variation is likely to be related to the presence of vectors of African animal trypanosomosis. Partial Mantel test identified the present Northern limit of trypanosome vectors as the most significant landscape boundary influencing the genetic variability of Burkina Faso goat (p = 0.008). The contribution of Sahel goat genes to the goat populations in the Northern and Eastern parts of the Sudan-Sahel area of Burkina Faso was substantial. The presence of perennial streams explains the existence of trypanosome vectors. The South half of the Nakambé river (Southern Ouagadougou) and the Mouhoun river loop determined, respectively, the Eastern and Northern limits for the expansion of Sahelian goat genes. Furthermore, results from partial Mantel test suggest that the introgression of Sahelian goat genes into Djallonké goat using human-influenced genetic corridors has a limited influence when compared to the biological boundary defined by the northern limits for the distribution of the tsetse fly. However, the genetic differences found between the goat sampled in Bobo Dioulasso and the other populations located in the Sudan area of Burkina Faso may be explained by the broad goat trade favoured by the main road of the country.

Conclusions

The current analysis clearly suggests that genetic variation in Burkina Faso goat: a) follows a North to South clinal; and b) is affected by the distribution of the tsetse fly that imposes a limit to the Sahelian goat expansion due to their trypanosusceptibility. Here we show how extensive surveys on livestock populations can be useful to indirectly assess the consequences of climate change and human action in developing countries.

Genomic scan of selective sweeps in thin and fat tail sheep breeds for identifying of candidate regions associated with fat deposition

BACKGROUND

Identification of genomic regions that have been targets of selection for phenotypic traits is one of the most important and challenging areas of research in animal genetics. However, currently there are relatively few genomic regions identified that have been subject to positive selection. In this study, a genome-wide scan using ~50,000 Single Nucleotide Polymorphisms (SNPs) was performed in an attempt to identify genomic regions associated with fat deposition in fat-tail breeds. This trait and its modification are very important in those countries grazing these breeds.

RESULTS

Two independent experiments using either Iranian or Ovine HapMap genotyping data contrasted thin and fat tail breeds. Population differentiation using FST in Iranian thin and fat tail breeds revealed seven genomic regions. Almost all of these regions overlapped with QTLs that had previously been identified as affecting fat and carcass yield traits in beef and dairy cattle. Study of selection sweep signatures using FST in thin and fat tail breeds sampled from the Ovine HapMap project confirmed three of these regions located on Chromosomes 5, 7 and X. We found increased homozygosity in these regions in favour of fat tail breeds on chromosome 5 and X and in favour of thin tail breeds on chromosome 7.

CONCLUSIONS

In this study, we were able to identify three novel regions associated with fat deposition in thin and fat tail sheep breeds. Two of these were associated with an increase of homozygosity in the fat tail breeds which would be consistent with selection for mutations affecting fat tail size several thousand years after domestication.

Farmers' perceptions on trypanosomosis and trypanotolerance character of the taurine Sheko

In the humid and subhumid tropics, trypanosomosis is an economically important zoonotic protozoan disease of the commonly kept farm animal species and their wild relatives. For example, more than 20% of the humid western and southwestern Ethiopia, which is home to more than 14 million heads of cattle, is under varying levels of trypanosomosis risk. Our study was, therefore, initiated to document farmers' perception on trypanosomosis and Sheko's trypanotolerance character. Our findings showed that trypanosomosis was the most frequently reported cattle disease in the Bench Maji Zone. Accordingly, 76.7% of the farmers reported the epidemiological importance of trypanosomosis, and they also noted that trypanosomosis on average accounted for 63.0% of annualized cattle death. The reported signs of trypanosomosis and trypanotolerance indicators were consistent with literature reports. Moreover, 66.7% of the farmers reported Sheko's trypanotolerance character. In the course of time, smallholder farmers have developed ethnoveterinary practices that are mainly used to prevent the landing of vector flies on the animal. Wet and warm seasons of the year, i.e., spring and, to some extent, the beginning of summer and autumn, were reported as peak periods of trypanosomosis risk. Therefore, this showed the need for incorporating farmers' knowledge in trypanosomosis control programs.

Footprints of selection in the ancestral admixture of a New World Creole cattle breed

Admixed populations represent attractive biological models to study adaptive selection. Originating from several waves of recent introduction from European (EUT), African (AFT) and zebus (ZEB) cattle, New World Creole cattle allow investigating the response to tropical environmental challenges of these three ancestries. We here provide a detailed assessment of their genetic contributions to the Creole breed from Guadeloupe (CGU). We subsequently look for footprints of selection by combining results from tests based on the extent of haplotype homozygosity and the identification of excess/deficiency of local ancestry. To tackle these issues, 140 CGU individuals and 25 Brahman zebus from Martinique were genotyped at 44 057 SNPs. These data were combined to those available on 23 populations representative of EUT, AFT or ZEB. We found average proportions of 26.1%, 36.0% and 37.9% of EUT, AFT and ZEB ancestries in the CGU genome indicating a higher level of African and zebu ancestries than suggested by historical records. We further identified 23 genomic regions displaying strong signal of selection, most of them being characterized by an excess of ZEB local ancestry. Among the candidate gene underlying these regions, several are associated with reproductive functions (RXFP2, PMEPA1, IGFBP3, KDR, PPP1R8, TBXA2R and SLC7A5) and metabolism (PDE1B and CYP46A1). Finally, two genes (CENTD3 and SAMD12) are involved in cellular signalization of immune response. This study illustrates the relevance of admixed populations to identify footprints of selection by combining several tests straightforward to implement on large data sets.

A whole genome Bayesian scan for adaptive genetic divergence in West African cattle

BACKGROUND

The recent settlement of cattle in West Africa after several waves of migration from remote centres of domestication has imposed dramatic changes in their environmental conditions, in particular through exposure to new pathogens. West African cattle populations thus represent an appealing model to unravel the genome response to adaptation to tropical conditions. The purpose of this study was to identify footprints of adaptive selection at the whole genome level in a newly collected data set comprising 36,320 SNPs genotyped in 9 West African cattle populations.

RESULTS

After a detailed analysis of population structure, we performed a scan for SNP differentiation via a previously proposed Bayesian procedure including extensions to improve the detection of loci under selection. Based on these results we identified 53 genomic regions and 42 strong candidate genes. Their physiological functions were mainly related to immune response (MHC region which was found under strong balancing selection, CD79A, CXCR4, DLK1, RFX3, SEMA4A, TICAM1 and TRIM21), nervous system (NEUROD6, OLFM2, MAGI1, SEMA4A and HTR4) and skin and hair properties (EDNRB, TRSP1 and KRTAP8-1).

CONCLUSION

The main possible underlying selective pressures may be related to climatic conditions but also to the host response to pathogens such as Trypanosoma(sp). Overall, these results might open the way towards the identification of important variants involved in adaptation to tropical conditions and in particular to resistance to tropical infectious diseases.

Divergent antimicrobial peptide (AMP) and acute phase protein (APP) responses to Trypanosoma congolense infection in trypanotolerant and trypanosusceptible cattle

African animal trypanosomiasis (AAT) is endemic across Sub-Saharan African and is a major constraint to livestock production. The ability of certain cattle breeds to remain productive despite infection is known as trypanotolerance; however, the underlying immune mechanisms contributing to this trait remain poorly understood. Antimicrobial peptides (AMPs) and acute phase proteins (APPs) are evolutionarily conserved effector molecules of the innate immune system that have important roles in the resolution of infection and activation of the adaptive immune response. Expression levels of AMP genes (TAP, LAP, BNBD4, DEFB1, DEFB5 and LEAP2) and APP genes (HP, CP, AGP, LBP, SAA3 and CRP) were investigated using real time quantitative reverse transcription PCR (qRT-PCR) in peripheral blood mononuclear cells (PBMC) isolated from two breeds of African cattle (trypanotolerant N'Dama and trypanosusceptible Boran), experimentally infected with Trypanosoma congolense. Haptoglobin and serum amyloid A (SAA) were also measured in plasma using quantitative protein assays. Results demonstrated that tracheal antimicrobial peptide (TAP) gene expression increased by 32-fold in Boran, compared to only 3-fold in N'Dama, by 14 days post-infection (dpi) and rising to 136-fold at 29 dpi in Boran, compared to 47-fold in N'Dama (P<0.05). Protein expression levels of SAA are elevated in N'Dama, rising to 163 microg/ml at 14 dpi compared with 72 microg/ml in Boran. The SAA expression profile mirrors the wave of parasitaemia detected in N'Dama. Seven single nucleotide polymorphisms (SNPs) were identified in the promoter regions of the SAA3 and SAA4 genes, which are predicted to affect transcription factor binding and thereby contributing to the differential patterns of expression detected between the breeds. Whereas elevated TAP expression is a conserved component of the innate immune response to infection in both breeds, higher SAA expression levels may contribute toward trypanotolerance in N'Dama.

Population genetic structure of wild and farmed rusa deer (Cervus timorensis russa) in New-Caledonia inferred from polymorphic microsatellite loci

Historical records indicate that 12 rusa deer (Cervus timorensis russa) were introduced in New-Caledonia during the 1870s. We used eight polymorphic microsatellite DNA loci to assess the genetic differentiation and diversity of farmed and wild deer populations. Past genetic bottlenecks were detected in both sub-populations, although higher genetic diversity was maintained in farmed populations, probably due to the regular introduction of reproducers from wild populations and from other farms. The genetic structure of farmed and wild populations differed significantly. There was a significant isolation by distance for wild populations, whereas farmed populations were significantly differentiated between farms independently from their geographical proximity. Wild rusa deer consisted of small populations (with effective population sizes ranging between 7 and 19 individuals depending on the methods used), with a low parent-offspring dispersion range (0.20-2.02 km). Genetic tools and direct observations provided congruent estimates of dispersion and population sizes. We discuss the relevance of our results for management purposes.