Genome Sequence Analysis Reveals Selection Signatures in Endangered Trypanotolerant West African Muturu Cattle

Genomic insights into adaptation and inbreeding among Sub-Saharan African cattle from pastoral and agropastoral systems

Background

In Sub-Saharan Africa (SSA), cattle are crucial for socioeconomic stability yet face numerous environmental stressors such as diseases, parasites, and extreme heat within pastoral and agropastoral systems. Despite their significance, gaps remain in understanding how genetic diversity and inbreeding influence traits essential for disease resistance and environmental adaptability. This study examines the genomic adaptations that enable SSA cattle to thrive under these conditions and assesses the impact of inbreeding on such adaptive traits.

Methods

We analyzed genomic data from 113 cattle across four breeds-Kuri, N'dama, Zebu-Fulani, and Zebu-Bororo-employing Runs of Homozygosity (ROH) and Integrated Haplotype Score (iHS) analyses to identify historical and recent genetic selections. Strict quality controls using PLINK software ensured accurate genomic pattern identification related to adaptation and inbreeding.

Results

ROH analysis revealed islands with genes such as RSAD2, CMPK2, and NOTCH1, which are involved in immune response and cellular stress management, highlighting regions of historical selection that have likely provided adaptive advantages in overcoming environmental and pathogenic stresses. In contrast, iHS analysis identified genes under recent selection like HIPK1, involved in stress response regulation, and EPHA5, which plays a crucial role in neural development and synaptic functions, potentially equipping these breeds with novel adaptations to ongoing and emergent environmental challenges.

Conclusion

This research confirms that selective pressures inherent in pastoral and agropastoral systems profoundly influence the genetic structure of SSA cattle. By delineating the genetic bases of key adaptive traits, our study offers crucial insights for targeted breeding programs to enhance cattle resilience and productivity. These findings provide a valuable framework for future genetic improvements and conservation strategies, crucial for sustainable livestock management and economic stability in SSA.

Genomic Reference Resource for African Cattle: Genome Sequences and High-Density Array Variants

The diversity in genome resources is fundamental to designing genomic strategies for local breed improvement and utilisation. These resources also support gene discovery and enhance our understanding of the mechanisms of resilience with applications beyond local breeds. Here, we report the genome sequences of 555 cattle (208 of which comprise new data) and high-density (HD) array genotyping of 1,082 samples (537 new samples) from indigenous African cattle populations. The new sequences have an average genome coverage of ~30X, three times higher than the average (~10X) of the over 300 sequences already in the public domain. Following variant quality checks, we identified approximately 32.3 million sequence variants and 661,943 HD autosomal variants mapped to the Bos taurus reference genome (ARS-UCD1.2). The new datasets were generated as part of the Centre for Tropical Livestock Genetics and Health (CTLGH) Genomic Reference Resource for African Cattle (GRRFAC) initiative, which aspires to facilitate the generation of this livestock resource and hopes for its utilisation for complete indigenous breed characterisation and sustainable global livestock improvement.

The DRB3 gene of the bovine major histocompatibility complex: discovery, diversity and distribution of alleles in commercial breeds of cattle and applications for development of vaccines

The bovine Major Histocompatibility Complex (MHC), also known as the Bovine Leucocyte Antigen (BoLA) complex, is the genomic region that encodes the most important molecules for antigen presentation to initiate immune responses. The first evidence of MHC in bovines pointed to a locus containing 2 antigens, one detected by cytotoxic antiserum (MHC class I) and another studied by mixed lymphocyte culture tests (MHC class II). The most studied gene in the BoLA region is the highly polymorphic BoLA-DRB3, which encodes a β chain with a peptide groove domain involved in antigen presentation for T cells that will develop and co-stimulate cellular and humoral effector responses. BoLA-DRB3 alleles have been associated with outcomes in infectious diseases such as mastitis, trypanosomiasis, and tick loads, and with production traits. To catalog these alleles, 2 nomenclature methods were proposed, and the current use of both systems makes it difficult to list, comprehend and apply these data effectively. In this review we have organized the knowledge available in all of the reports on the frequencies of BoLA-DRB3 alleles. It covers information from studies made in at least 26 countries on more than 30 breeds; studies are lacking in countries that are important producers of cattle livestock. We highlight practical applications of BoLA studies for identification of markers associated with resistance to infectious and parasitic diseases, increased production traits and T cell epitope mapping, in addition to genetic diversity and conservation studies of commercial and creole and locally adapted breeds. Finally, we provide support for the need of studies to discover new BoLA alleles and uncover unknown roles of this locus in production traits.

Whole-genome resequencing provides insights into the diversity and adaptation to desert environment in Xinjiang Mongolian cattle

BACKGROUND

Xinjiang Mongolian cattle is an indigenous breed that inhabits the Taklimakan Desert and is characterized by its small body size. However, the genomic diversity, origin, and genetic basis underlying the adaptation to the desert environment have been poorly studied.

RESULTS

We analyzed patterns of Xinjiang Mongolian cattle genetic variation by sequencing 20 genomes together with seven previously sequenced genomes and comparing them to the 134 genomes of nine representative breeds worldwide. Among the breeds of Bos taurus, we found the highest nucleotide diversity (0.0024) associated with the lower inbreeding coefficient (2.0110-6), the lowest linkage disequilibrium (r2 = 0.3889 at distance of 10 kb), and the highest effective population size (181 at 20 generations ago) in Xinjiang Mongolian cattle. The genomic diversity pattern could be explained by a limited introgression of Bos indicus genes. More importantly, similarly to desert-adapted camel and same-habitat sheep, we also identified signatures of selection including genes, GO terms, and/or KEGG pathways controlling water reabsorption and osmoregulation, metabolic regulation and energy balance, as well as small body size in Xinjiang Mongolian cattle.

CONCLUSIONS

Our results imply that Xinjiang Mongolian cattle might have acquired distinct genomic diversity by virtue of the introgression of Bos indicus, which helps understand the demographic history. The identification of selection signatures can provide novel insights into the genomic basis underlying the adaptation of Xinjiang Mongolian cattle to the desert environment.

Ancestral alleles defined for 70 million cattle variants using a population-based likelihood ratio test

BACKGROUND

The study of ancestral alleles provides insights into the evolutionary history, selection, and genetic structures of a population. In cattle, ancestral alleles are widely used in genetic analyses, including the detection of signatures of selection, determination of breed ancestry, and identification of admixture. Having a comprehensive list of ancestral alleles is expected to improve the accuracy of these genetic analyses. However, the list of ancestral alleles in cattle, especially at the whole genome sequence level, is far from complete. In fact, the current largest list of ancestral alleles (~ 42 million) represents less than 28% of the total number of detected variants in cattle. To address this issue and develop a genomic resource for evolutionary studies, we determined ancestral alleles in cattle by comparing prior derived whole-genome sequence variants to an out-species group using a population-based likelihood ratio test.

RESULTS

Our study determined and makes available the largest list of ancestral alleles in cattle to date (70.1 million) and includes 2.3 million on the X chromosome. There was high concordance (97.6%) of the determined ancestral alleles with those from previous studies when only high-probability ancestral alleles were considered (29.8 million positions) and another 23.5 million high-confidence ancestral alleles were novel, expanding the available reference list to improve the accuracies of genetic analyses involving ancestral alleles. The high concordance of the results with previous studies implies that our approach using genomic sequence variants and a likelihood ratio test to determine ancestral alleles is appropriate.

CONCLUSIONS

Considering the high concordance of ancestral alleles across studies, the ancestral alleles determined in this study including those not previously listed, particularly those with high-probability estimates, may be used for further genetic analyses with reasonable accuracy. Our approach that used predetermined variants in species and the likelihood ratio test to determine ancestral alleles is applicable to other species for which sequence level genotypes are available.

Candidate signatures of positive selection for environmental adaptation in indigenous African cattle: A review

Environmental adaptation traits of indigenous African cattle are increasingly being investigated to respond to the need for sustainable livestock production in the context of unpredictable climatic changes. Several studies have highlighted genomic regions under positive selection probably associated with adaptation to environmental challenges (e.g. heat stress, trypanosomiasis, tick and tick-borne diseases). However, little attention has focused on pinpointing the candidate causative variant(s) controlling the traits. This review compiled information from 22 studies on signatures of positive selection in indigenous African cattle breeds to identify regions under positive selection. We highlight some key candidate genome regions and genes of relevance to the challenges of living in extreme environments (high temperature, high altitude, high infectious disease prevalence). They include candidate genes involved in biological pathways relating to innate and adaptive immunity (e.g. BoLAs, SPAG11, IL1RL2 and GFI1B), heat stress (e.g. HSPs, SOD1 and PRLH) and hypoxia responses (e.g. BDNF and INPP4A). Notably, the highest numbers of candidate regions are found on BTA3, BTA5 and BTA7. They overlap with genes playing roles in several biological functions and pathways. These include but are not limited to growth and feed intake, cell stability, protein stability and sweat gland development. This review may further guide targeted genome studies aiming to assess the importance of candidate causative mutations, within regulatory and protein-coding genome regions, to further understand the biological mechanisms underlying African cattle's unique adaption.

Assessing genetic diversity and defining signatures of positive selection on the genome of dromedary camels from the southeast of the Arabian Peninsula

Dromedary camels (Camelus dromedarius) are members of the Camelini tribe within the Camelidae family. They are distributed throughout North Africa, the Arabian Peninsula and Southeast Asia. This domestic species is characterized by its superior adaptability to the harsh desert environment. In this study, whole autosomal data of 29 dromedary samples from the Southeast Arabian Peninsula in Oman; 10 from Muscat, 14 from Al-Batinah, and 5 from Al-Sharqiya, were investigated to assess their genetic relationship and to define candidate signatures of positive selection. A minimal genetic distinction that separates Muscat dromedaries from the other two populations was observed, with a degree of genetic admixture between them. Using the de-correlated composite of multiple signals (DCMS) approach, a total of 47 candidate regions within the autosomes of these dromedary populations were defined with signatures of positive selection. These candidate regions harbor a total of 154 genes that are mainly associated with functional categories related to immune response, lipid metabolism and energy expenditure, optical and auditory functions, and long-term memory. Different functional genomic variants were called on the candidate regions and respective genes that warrant further investigation to find possible association with the different favorable phenotypes in dromedaries. The output of this study paves the way for further research efforts aimed at defining markers for use in genomic breeding programs, with the goal of conserving the genetic diversity of the species and enhancing its productivity.

Candidate genes associated with heat stress and breeding strategies to relieve its effects in dairy cattle: a deeper insight into the genetic architecture and immune response to heat stress

The need for food products of animal origin is increasing worldwide. Satisfying these needs in a way that has minimal impact on the environment requires cutting-edge technologies and techniques to enhance the genetic quality of cattle. Heat stress (HS), in particular, is affecting dairy cattle with increasing frequency and severity. As future climatic challenges become more evident, identifying dairy cows that are more tolerant to HS will be important for breeding dairy herds that are better adapted to future environmental conditions and for supporting the sustainability of dairy farming. While research into the genetics of HS in the context of the effect of global warming on dairy cattle is gaining momentum, the specific genomic regions involved in heat tolerance are still not well documented. Advances in omics information, QTL mapping, transcriptome profiling and genome-wide association studies (GWAS) have identified genomic regions and variants associated with tolerance to HS. Such studies could provide deeper insights into the genetic basis for response to HS and make an important contribution to future breeding for heat tolerance, which will help to offset the adverse effects of HS in dairy cattle. Overall, there is a great interest in identifying candidate genes and the proportion of genetic variation associated with heat tolerance in dairy cattle, and this area of research is currently very active worldwide. This review provides comprehensive information pertaining to some of the notable recent studies on the genetic architecture of HS in dairy cattle, with particular emphasis on the identified candidate genes associated with heat tolerance in dairy cattle. Since effective breeding programs require optimal knowledge of the impaired immunity and associated health complications caused by HS, the underlying mechanisms by which HS modulates the immune response and renders animals susceptible to various health disorders are explained. In addition, future breeding strategies to relieve HS in dairy cattle and improve their welfare while maintaining milk production are discussed.

Genomic Characterization and Initial Insight into Mastitis-Associated SNP Profiles of Local Latvian Bos taurus Breeds

Latvia has two local Bos taurus breeds-Latvian Brown (LBG) and Latvian Blue (LZG)-characterized by a good adaptation to the local climate, longevity, and high fat and protein contents in milk. Since these are desired traits in the dairy industry, this study investigated the genetic background of the LBG and LZG breeds and identified the genetic factors associated with mastitis. Blood and semen samples were acquired, and whole genome sequencing was then performed to acquire a genomic sequence with at least 35× or 10× coverage. The heterozygosity, nucleotide diversity, and LD analysis indicated that LBG and LZG cows have similar levels of genetic diversity compared to those of other breeds. An analysis of the population structure revealed that each breed clustered together, but the overall differentiation between the breeds was small. The highest genetic variance was observed in the LZG breed compared with the LBG breed. Our results show that SNP rs721295390 is associated with mastitis in the LBG breed, and SNPs rs383806754, chr29:43998719CG>C, and rs462030680 are associated with mastitis in the LZG breed. This study shows that local Latvian LBG and LZG breeds have a pronounced genetic differentiation, with each one suggesting its own mastitis-associated SNP profile.

Genome-wide diversity and admixture of five indigenous cattle populations from the Tigray region of northern Ethiopia

The Tigray region, where we found around eight per cent of the indigenous cattle population of Ethiopia, is considered as the historic centre of the country, with the ancient pre-Aksumite and Aksumite civilisations in contact with the civilisations of the Fertile Crescent and the Indian subcontinent. Here, we used whole genome sequencing data to characterise the genomic diversity, relatedness, and admixture of five cattle populations (Abergelle, Arado, Begait, Erob, and Raya) indigenous to the Tigray region of Ethiopia. We detected 28 to 29 million SNPs and 2.7 to 2.9 million indels in each population, of which 7% of SNPs and 34% of indels were novel. Functional annotation of the variants showed around 0.01% SNPs and 0.22%-0.27% indels in coding regions. Enrichment analysis of genes overlapping missense private SNPs revealed 20 significant GO terms and KEGG pathways that were shared by or specific to breeds. They included important genes associated with morphology (SCN4A, TAS1R2 and KCNG4), milk yield (GABRG1), meat quality (MMRN2, VWC2), feed efficiency (PCDH8 and SLC26A3), immune response (LAMC1, PCDH18, CELSR1, TLR6 and ITGA5), heat resistance (NPFFR1 and HTR7) and genes belonging to the olfactory gene family, which may be related to adaptation to harsh environments. Tigray indigenous cattle are very diverse. Their genome-wide average nucleotide diversity ranged from 0.0035 to 0.0036. The number of heterozygous SNPs was about 0.6-0.7 times higher than homozygous ones. The within-breed average number of ROHs ranged from 777.82 to 1000.45, with the average sum of the length of ROHs ranging from 122.01 Mbp to 163.88 Mbp. The genomic inbreeding coefficients differed among animals and breeds, reaching up to 10% in some Begait and Raya animals. Tigray indigenous cattle shared a common ancestry with Asian indicine (85.6%-88.7%) and African taurine (11.3%-14.1%) cattle, with very small, if any, European taurine introgression. This study identified high within-breed genetic diversity representing an opportunity for breeding improvement programs and, also, significant novel variants that could increase the number of known cattle variants, an important contribution to the knowledge of domestic cattle genetic diversity.

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.

Molecular genotyping reveals mixed bovine and human trypanosomiasis in cattle from West Africa

Background and Aim

Animal trypanosomiasis is a major contributor to agricultural and economic losses, especially in sub-Saharan Africa. We have shown that some animal species expressed genes that are significant players in immune response to bovine trypanosomosis, impeding signs and symptoms of the disease. We hypothesize that such animals are contributors to disease transmission dynamics and severe outcomes. Therefore, this study aims to ascertain trypanosome species diversity in cattle and their potential role as reservoirs for the transmission of human disease.

Materials and Methods

We performed a molecular genotyping of trypanosome internal transcribed spacer 1 (ITS-1) and 18S ribosomal RNA genes on genomic DNA extracts from randomly sampled N'Dama cattle from slaughterhouses in Nigeria. We identified trypanosome species circulating among the animals through polymerase chain reaction and genomic sequencing. We performed multiple sequence alignments as well as conducted a phylogenetic relationship between identified species.

Results

In all, 9 of 127 (7.1%) samples were positively amplified (band sizes ranging from 250 bp to 710 bp), including an isolate with two distinct bands (700 and 710 bp), indicating two trypanosome types. Sequence similarity and homology analysis identified four species, namely: Trypanosoma vivax, Trypanosoma congolense forest type, T. congolense savannah type, and Trypanosoma brucei. Interestingly, one of the bands, additionally verified by nucleotide sequencing, was identified as a human trypanosome (Trypanosoma brucei gambiense), confirming our hypothesis that cattle are potential reservoir hosts for human trypanosomes.

Conclusion

Overall, we observed different trypanosome species in our study area, with animals on the same farm infected with multiple species, which could complicate treatment and disease control strategies. Finding human trypanosome species strengthens the argument that disease transmission dynamics are modulated by other vertebrates, further complicating control programs.

Identification of genomic regions associated with total and progressive sperm motility in Italian Holstein bulls

Sperm motility is directly related to the ability of sperm to move through the female reproductive tract to reach the ovum. Sperm motility is a complex trait that is influenced by environmental and genetic factors and is associated with male fertility, oocyte penetration rate, and reproductive success of cattle. In this study we carried out a GWAS in Italian Holstein bulls to identify candidate regions and genes associated with variations in progressive and total motility (PM and TM, respectively). After quality control, the final data set consisted of 5,960 records from 949 bulls having semen collected in 10 artificial insemination stations and genotyped at 412,737 SNPs (call rate >95%; minor allele frequency >5%). (Co)variance components were estimated using single trait mixed models, and associations between SNPs and phenotypes were assessed using a genomic BLUP approach. Ten windows that explained the greatest percentage of genetic variance were located on Bos taurus autosomes 1, 2, 4, 6, 7, 23, and 26 for TM and Bos taurus autosomes 1, 2, 4, 6, 8, 16, 23, and 26 for PM. A total of 150 genes for TM and 72 genes for PM were identified within these genomic regions. Gene Ontology enrichment analyses identified significant Gene Ontology terms involved with energy homeostasis, membrane functions, sperm-egg interactions, protection against oxidative stress, olfactory receptors, and immune system. There was significant enrichment of quantitative trait loci for fertility, calving ease, immune response, feed intake, and carcass weight within the candidate windows. These results contribute to understanding the architecture of the genetic control of sperm motility and may aid in the development of strategies to identify subfertile bulls and improve reproductive success.

Population Structure and Selection Signatures Underlying Domestication Inferred from Genome-Wide Copy Number Variations in Chinese Indigenous Pigs

Single nucleotide polymorphism was widely used to perform genetic and evolution research in pigs. However, little is known about the effect of copy number variation (CNV) on characteristics in pigs. This study performed a genome-wide comparison of CNVs between Wannan black pigs (WBP) and Asian wild boars (AWB), using whole genome resequencing data. By using Manta, we detected in total 28,720 CNVs that covered approximately 1.98% of the pig genome length. We identified 288 selected CNVs (top 1%) by performing Fst statistics. Functional enrichment analyses for genes located in selected CNVs were found to be muscle related (NDN, TMOD4, SFRP1, and SMYD3), reproduction related (GJA1, CYP26B1, WNT5A, SRD5A2, PTPN11, SPEF2, and CCNB1), residual feed intake (RFI) related (MAP3K5), and ear size related (WIF1). This study provides essential information on selected CNVs in Wannan black pigs for further research on the genetic basis of the complex phenotypic and provides essential information for direction in the protection and utilization of Wannan black pig.

Genomic signatures for drylands adaptation at gene-rich regions in African zebu cattle

Background

Indigenous Sudanese cattle are mainly indicine/zebu (humped) type. They thrive in the harshest dryland environments characterised by high temperatures, long seasonal dry periods, nutritional shortages, and vector disease challenges. Here, we sequenced 60 indigenous Sudanese cattle from six indigenous breeds and analysed the data using three genomic scan approaches to unravel cattle adaptation to the African dryland region.

Results

We identified a set of gene-rich selective sweep regions, detected mostly on chromosomes 5, 7 and 19, shared across African and Gir zebu. These include genes involved in immune response, body size and conformation, and heat stress response. We also identified selective sweep regions unique to Sudanese zebu. Of these, a 250 kb selective sweep on chromosome 16 spans seven genes, including PLCH2, PEX10, PRKCZ, and SKI, which are involved in alternative adaptive metabolic strategies of insulin signalling, glucose homeostasis, and fat metabolism.

Conclusions

Our results suggest that environmental adaptation may involve recent and ancient selection at gene-rich regions, which might be under a common regulatory genetic control, in zebu cattle.

•

Sudanese cattle thrive in the harshest environments of the African drylands.

•

Bos indicus shared selected genes are involved in immune response, conformation, and heat stress response.

•

Sudanese zebu-specific sweep includes genes involved in alternative adaptive metabolic strategies of insulin signalling, glucose homeostasis, and fat metabolism.

•

Environmental adaptation in zebu cattle may involve recent and ancient selection at gene-rich regions, which might be under a common regulatory genetic control.

Assessing genomic diversity and signatures of selection in Pinan cattle using whole-genome sequencing data

BACKGROUND

Crossbreeding is an important way to improve production beef cattle performance. Pinan cattle is a new hybrid cattle obtained from crossing Piedmontese bulls with Nanyang cows. After more than 30 years of cross-breeding, Pinan cattle show a variety of excellent characteristics, including fast growth, early onset of puberty, and good meat quality. In this study, we analyzed the genetic diversity, population structure, and genomic region under the selection of Pinan cattle based on whole-genome sequencing data of 30 Pinan cattle and 169 published cattle genomic data worldwide.  RESULTS: Estimating ancestry composition analysis showed that the composition proportions for our Pinan cattle were mainly Piedmontese and a small amount of Nanyang cattle. The analyses of nucleotide diversity and linkage disequilibrium decay indicated that the genomic diversity of Pinan cattle was higher than that of European cattle and lower than that of Chinese indigenous cattle. De-correlated composite of multiple selection signals, which combines four different statistics including θπ, CLR, F_ST, and XP-EHH, was computed to detect the signatures of selection in the Pinan cattle genome. A total of 83 genes were identified, affecting many economically important traits. Functional annotation revealed that these selected genes were related to immune (BOLA-DQA2, BOLA-DQB, LSM14A, SEC13, and NAALADL2), growth traits (CYP4A11, RPL26, and MYH10), embryo development (REV3L, NT5E, CDX2, KDM6B, and ADAMTS9), hornless traits (C1H21orf62), and climate adaptation (ANTXR2).

CONCLUSION

In this paper, we elucidated the genomic characteristics, ancestry composition, and selective signals related to important economic traits in Pinan cattle. These results will provide the basis for further genetic improvement of Pinan cattle and reference for other hybrid cattle related studies.

Deciphering the Molecular Mechanism Underlying African Animal Trypanosomiasis by Means of the 1000 Bull Genomes Project Genomic Dataset

Simple Summary

Climate change is increasing the risk of spreading vector-borne diseases such as African Animal Trypanosomiasis (AAT), which is causing major economic losses, especially in sub-Saharan African countries. Mainly considering this disease, we have investigated transcriptomic and genomic data from two cattle breeds, namely Boran and N‘Dama, where the former is known for its susceptibility and the latter one for its tolerance to the AAT. Despite the rich literature on this disease, there is still a need to investigate underlying genetic mechanisms to decipher the complex interplay of regulatory SNPs (rSNPs), their corresponding gene expression profiles and the downstream effectors associated with the AAT disease. The findings of this study complement our previous results, which mainly involve the upstream events, including transcription factors (TFs) and their co-operations as well as master regulators. Moreover, our investigation of significant rSNPs and effectors found in the liver, spleen and lymph node tissues of both cattle breeds could enhance the understanding of distinct mechanisms leading to either resistance or susceptibility of cattle breeds.

Abstract

African Animal Trypanosomiasis (AAT) is a neglected tropical disease and spreads by the vector tsetse fly, which carries the infectious Trypanosoma sp. in their saliva. Particularly, this parasitic disease affects the health of livestock, thereby imposing economic constraints on farmers, costing billions of dollars every year, especially in sub-Saharan African countries. Mainly considering the AAT disease as a multistage progression process, we previously performed upstream analysis to identify transcription factors (TFs), their co-operations, over-represented pathways and master regulators. However, downstream analysis, including effectors, corresponding gene expression profiles and their association with the regulatory SNPs (rSNPs), has not yet been established. Therefore, in this study, we aim to investigate the complex interplay of rSNPs, corresponding gene expression and downstream effectors with regard to the AAT disease progression based on two cattle breeds: trypanosusceptible Boran and trypanotolerant N’Dama. Our findings provide mechanistic insights into the effectors involved in the regulation of several signal transduction pathways, thereby differentiating the molecular mechanism with regard to the immune responses of the cattle breeds. The effectors and their associated genes (especially MAPKAPK5, CSK, DOK2, RAC1 and DNMT1) could be promising drug candidates as they orchestrate various downstream regulatory cascades in both cattle breeds.

Selection Signature Analyses Revealed Genes Associated With Adaptation, Production, and Reproduction in Selected Goat Breeds in Kenya

Artificial and natural selection in livestock is expected to leave unique footprints on their genomes. Goat breeds in Kenya have evolved for survival, breeding, and production in various harsh ecological areas, and their genomes are likely to have acquired unique alleles for adaptation to such diverse production environments and other traits of economic importance. To investigate signals of selection for some selected goat breeds in Kenya, Alpine (n = 29), Galla (n = 12), Saanen (n = 24), and Toggenburg (n = 31) were considered. A total of 53,347 single-nucleotide polymorphisms (SNPs) generated using the Illumina GoatSNP50 BeadChip were analyzed. After quality control, 47,663 autosomal single-nucleotide polymorphisms remained for downstream analyses. Several complementary approaches were applied for the following analyses: integrated Haplotype Score (iHS), cross-population-extended haplotype homozygosity (XP-EHH), hapFLK, and FLK. A total of 404 top genomic regions were identified across all the four breeds, based on the four complementary analyses. Out of the 16 identified putative selection signature regions by the intersection of multiple-selective signal analyses, most of the putative regions were found to overlap significantly with the iHS and XP-EHH analyses on chromosomes 3, 4, 10, 15, 22, and 26. These regions were enriched with some genes involved in pathways associated directly or indirectly with environmental adaptation regulating immune responses (e.g., HYAL1 and HYAL3), milk production (e.g., LEPR and PDE4B), and adaptability (e.g., MST1 and PCK). The results revealed few intersect between breeds in genomic selection signature regions. In general, this did not present the typical classic selection signatures as predicted due to the complex nature of the traits. The results support that some various selection pressures (e.g., environmental challenges, artificial selection, and genome admixture challenges) have molded the genome of goat breeds in Kenya. Therefore, the research provides new knowledge on the conservation and utilization of these goat genetic resources in Kenya. In-depth research is needed to detect precise genes connected with adaptation and production in goat breeds in Kenya.

Genome-wide investigations reveal the population structure and selection signatures of Nigerian cattle adaptation in the sub-Saharan tropics

Background

Cattle are considered to be the most desirable livestock by small scale farmers. In Africa, although comprehensive genomic studies have been carried out on cattle, the genetic variations in indigenous cattle from Nigeria have not been fully explored. In this study, genome-wide analysis based on genotyping-by-sequencing (GBS) of 193 Nigerian cattle was used to reveal new insights on the history of West African cattle and their adaptation to the tropical African environment, particularly in sub-Saharan region. 

Results

The GBS data were evaluated against whole-genome sequencing (WGS) data and high rate of variant concordance between the two platforms was evident with high correlated genetic distance matrices genotyped by both methods suggestive of the reliability of GBS applicability in population genetics. The genetic structure of Nigerian cattle was observed to be homogenous and unique from other African cattle populations. Selection analysis for the genomic regions harboring imprints of adaptation revealed genes associated with immune responses, growth and reproduction, efficiency of feeds utilization, and heat tolerance. Our findings depict potential convergent adaptation between African cattle, dogs and humans with adaptive genes SPRY2 and ITGB1BP1 possibly involved in common physiological activities.

Conclusion

The study presents unique genetic patterns of Nigerian cattle which provide new insights on the history of cattle in West Africa based on their population structure and the possibility of parallel adaptation between African cattle, dogs and humans in Africa which require further investigations.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08512-w.

The genetic history of Mayotte and Madagascar cattle breeds mirrors the complex pattern of human exchanges in Western Indian Ocean

Despite their central economic and cultural role, the origin of cattle populations living in Indian Ocean islands still remains poorly documented. Here, we unravel the demographic and adaptive histories of the extant Zebus from the Mayotte and Madagascar islands using high-density SNP genotyping data. We found that these populations are very closely related and both display a predominant indicine ancestry. They diverged in the 16th century at the arrival of European people who transformed the trade network in the area. Their common ancestral cattle population originates from an admixture between an admixed African zebu population and an Indian zebu that occurred around the 12th century at the time of the earliest contacts between human African populations of the Swahili corridor and Austronesian people from Southeast Asia in Comoros and Madagascar. A steep increase in the estimated population sizes from the beginning of the 16th to the 17th century coincides with the expansion of the cattle trade. By carrying out genome scans for recent selection in the two cattle populations from Mayotte and Madagascar, we identified sets of candidate genes involved in biological functions (cancer, skin structure, and UV-protection, nervous system and behavior, organ development, metabolism, and immune response) broadly representative of the physiological adaptation to tropical conditions. Overall, the origin of the cattle populations from Western Indian Ocean islands mirrors the complex history of human migrations and trade in this area.

Small genetic variation affecting mRNA isoforms associated with marbling and meat color in beef cattle

The aim of this study was to identify mRNA isoforms and small genetic variants that may be affecting marbling and beef color in Nellore cattle. Longissimus thoracis muscle samples from 20 bulls with different phenotypes (out of 80 bulls set) for marbling (moderate (n = 10) and low (n = 10) groups) and beef color (desirable (n = 10) and undesirable (n = 9) group) traits were used to perform transcriptomic analysis using RNA sequencing. Fourteen and 15 mRNA isoforms were detected as differentially expressed (DE) (P-value ≤ 0.001) between divergent groups for marbling and meat color traits, respectively. Some of those DE mRNA isoforms have shown sites of splicing modified by small structural variants as single nucleotide variant (SNV), insertion, and/or deletion. Enrichment analysis identified metabolic pathways, such as O2/CO2 exchange in erythrocytes, tyrosine biosynthesis, and phenylalanine degradation. The results obtained suggest potential key regulatory genes associated with these economically important traits for the beef industry and for the consumer.

Genetic diversity and population structure of four Nigerian indigenous cattle breeds

A total of thirty-eight (38) Nigerian indigenous cattle were used to evaluate the genetic relatedness, diversity, and population structure of four indigenous cattle breeds. Blood samples were collected from the experimental animals into vacutainer tubes containing EDTA, and genomic DNA extracted, quantified, and evaluated for integrity. Fourteen (14) microsatellite primers were used for polymerase chain reaction (PCR) and PCR amplification performed under standard conditions followed by electrophoresis in 2.5% Metaphor Agarose gel. Genomic parameter estimates included allele number (Na), observed (H_o) and expected (H_e) heterozygosity, polymorphism information content (PIC), test of Hardy-Weinberg equilibrium, and genetic diversity; pairwise Nei's genetic distance, Wright's F-statistics (FIT, FST, and FIS), and gene flow (Nm); and breed relationship, population structure, and degree of admixture. A total of 112 alleles were detected and mean number of alleles was 4.02 ± 0.190, while mean fixation index was 0.461 ± 0.068. Mean H_o and H_e were 0.352 ± 0.05 and 0.605 ± 0.018, respectively. Pairwise estimates of genetic differentiation, F_ST, were significantly different (p < 0.001) implying distinct breeds. Estimates of Nm were less than 4 but greater than 1, indicating that the cattle breeds do not belong to one panmictic population. Estimates of pairwise genetic distance revealed that White Fulani and Sokoto Gudali were more closely related than Muturu and N'Dama. The results of STRUCTURE, principal coordinate, and phylogenetic analyses revealed four clusters which implies that the breeds were genetically distinct. It is recommended that the four cattle breeds can be used to develop composites with higher genetic potentials for beef production and resistance to endemic diseases and pests. Further efforts should be made to conserve and genetically improve these breeds to meet present and future production and breeding imperatives.

Unraveling Admixture, Inbreeding, and Recent Selection Signatures in West African Indigenous Cattle Populations in Benin

The Dwarf Lagune and the Savannah Somba cattle in Benin are typical representatives of the endangered West African indigenous Shorthorn taurine. The Lagune was previously exported to African and European countries and bred as Dahomey cattle, whereas the Somba contributed to the formation of two indigenous hybrids known as Borgou and Pabli cattle. These breeds are affected by demographic, economic, and environmental pressures in local production systems. Considering current and historical genomic data, we applied a formal test of admixture, estimated admixture proportions, and computed genomic inbreeding coefficients to characterize the five breeds. Subsequently, we unraveled the most recent selection signatures using the cross-population extended haplotype homozygosity approach, based on the current and historical genotypes. Results from principal component analyses and high proportion of Lagune ancestry confirm the Lagune origin of the European Dahomey cattle. Moreover, the Dahomey cattle displayed neither indicine nor European taurine (EUT) background, but they shared on average 40% of autozygosity from common ancestors, dated approximately eight generations ago. The Lagune cattle presented inbreeding coefficients larger than 0.13; however, the Somba and the hybrids (Borgou and Pabli) were less inbred (≤0.08). We detected evidence of admixture in the Somba and Lagune cattle, but they exhibited a similar African taurine (AFT) ancestral proportion (≥96%) to historical populations, respectively. A moderate and stable AFT ancestral proportion (62%) was also inferred for less admixed hybrid cattle including the Pabli. In contrast, the current Borgou samples displayed a lower AFT ancestral proportion (47%) than historical samples (63%). Irrespective of the admixture proportions, the hybrid populations displayed more selection signatures related to economic traits (reproduction, growth, and milk) than the taurine. In contrast, the taurine, especially the Somba, presented several regions known to be associated with adaptive traits (immunity and feed efficiency). The identified subregion of bovine leukocyte antigen (BoLA) class IIb (including DSB and BOLA-DYA) in Somba cattle is interestingly uncommon in other African breeds, suggesting further investigations to understand its association with specific adaptation to endemic diseases in Benin. Overall, our study provides deeper insights into recent evolutionary processes in the Beninese indigenous cattle and their aptitude for conservation and genetic improvement.

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.

Genome-Wide Association Study Using Whole-Genome Sequence Data for Fertility, Health Indicator, and Endoparasite Infection Traits in German Black Pied Cattle

This genome-wide association study (GWAS) aimed to identify sequence variants (SVs) and candidate genes associated with fertility and health in endangered German Black Pied cattle (DSN) based on whole-genome sequence (WGS) data. We used 304 sequenced DSN cattle for the imputation of 1797 genotyped DSN to WGS. The final dataset included 11,413,456 SVs of 1886 cows. Cow traits were calving-to-first service interval (CTFS), non-return after 56 days (NR56), somatic cell score (SCS), fat-to-protein ratio (FPR), and three pre-corrected endoparasite infection traits. We identified 40 SVs above the genome-wide significance and suggestive threshold associated with CTFS and NR56, and three important potential candidate genes (ARHGAP21, MARCH11, and ZNF462). For SCS, most associations were observed on BTA 25. The GWAS revealed 61 SVs, a cluster of 10 candidate genes on BTA 13, and 7 pathways for FPR, including key mediators involved in milk fat synthesis. The strongest associations for gastrointestinal nematode and Dictyocaulus viviparus infections were detected on BTA 8 and 24, respectively. For Fasciola hepatica infections, the strongest associated SVs were located on BTA 4 and 7. We detected 200 genes for endoparasite infection traits, related to 16 pathways involved in host immune response during infection.

Human adaptation, demography and cattle domestication: an overview of the complexity of lactase persistence in Africa

Lactase persistence (LP) is a genetically-determined trait that is prevalent in African, European and Arab populations with a tradition of animal herding and milk consumption. To date, genetic analyses have identified several common variants that are associated with LP. Furthermore, data have indicated that these functional alleles likely have been maintained in pastoralist populations due to the action of recent selection, exemplifying the ongoing evolution of anatomically modern humans. Additionally, demographic history has also played a role in the geographic distribution of LP and associated alleles in Africa. In particular, the migration of ancestral herders and their subsequent admixture with local populations were integral to the spread of LP alleles and the culture of pastoralism across the continent. The timing of these demographic events was often correlated with known major environmental changes and/or the ability of domesticated cattle to resist/avoid infectious diseases. This review summarizes recent advances in our understanding of the genetic basis and evolutionary history of LP, as well as the factors that influenced the origin and spread of pastoralism in Africa.

Population structure, inbreeding and admixture in local cattle populations managed by community-based breeding programs in Burkina Faso

High-throughput genomic markers provide an opportunity to assess important indicators of genetic diversity for populations managed in livestock breeding programs. While well-structured breeding programs are common in developed countries, in developing country situations, especially in West Africa, on-farm performance and pedigree recordings are rare, and thus, genomic markers provide insights to the levels of genetic diversity, inbreeding and introgression by other breeds. In this study, we analysed key population parameters such as population structure, admixture and levels of inbreeding in three neighbouring populations of African taurine and taurine × Zebu crosses managed by community-based breeding programs in the South-West of Burkina Faso. The three populations were pure Baoulé (called Lobi locally) in sedentary production systems, Baoulé x Zebu crossbreds in sedentary systems and Zebu × Baoulé crossbreds in transhumant production systems, respectively. The total sample analysed included 631 animals and 38,207 single nucleotide polymorphisms after quality control. Results of principal component and admixture analyses confirmed the genetic background of two distinct ancestral populations (taurine and zebuine) and levels of admixture in all three breeding populations, including the presumably pure Baoulé group of animals. Inbreeding levels were moderate, compared to European dairy and beef cattle populations and higher than those of Brazilian Nellore cattle. Very few animals with inbreeding levels indicating parent-offspring or full sib mating were observed, and inbreeding levels indicating half sib mating were also rare. For the management of breeding populations, farmers were advised to exchange best young bulls. The crossbreeding levels of presumably pure Baoulé animals are of concern to the breeding program due to the high level of endangerment of pure African taurine cattle populations across West Africa. Future rounds of bull selection in the community-based breeding program will make use of genomic information about admixture levels.

Investigation of ancestral alleles in the Bovinae subfamily

BACKGROUND

In evolutionary theory, divergence and speciation can arise from long periods of reproductive isolation, genetic mutation, selection and environmental adaptation. After divergence, alleles can either persist in their initial state (ancestral allele - AA), co-exist or be replaced by a mutated state (derived alleles -DA). In this study, we aligned whole genome sequences of individuals from the Bovinae subfamily to the cattle reference genome (ARS.UCD-1.2) for defining ancestral alleles necessary for selection signatures study.

RESULTS

Accommodating independent divergent of each lineage from the initial ancestral state, AA were defined based on fixed alleles on at least two groups of yak, bison and gayal-gaur-banteng resulting in ~ 32.4 million variants. Using non-overlapping scanning windows of 10 Kb, we counted the AA observed within taurine and zebu cattle. We focused on the extreme points, regions with top 0. 1% (high count) and regions without any occurrence of AA (null count). High count regions preserved gene functions from ancestral states that are still beneficial in the current condition, while null counts regions were linked to mutated ones. For both cattle, high count regions were associated with basal lipid metabolism, essential for survival of various environmental pressures. Mutated regions were associated to productive traits in taurine, i.e. higher metabolism, cell development and behaviors and in immune response domain for zebu.

CONCLUSIONS

Our findings suggest that retaining and losing AA in some regions are varied and made it species-specific with possibility of overlapping as it depends on the selective pressure they had to experience.

Signatures of Selection in Composite Vrindavani Cattle of India

Vrindavani is an Indian composite cattle breed developed by crossbreeding taurine dairy breeds with native indicine cattle. The constituent breeds were selected for higher milk production and adaptation to the tropical climate. However, the selection response for production and adaptation traits in the Vrindavani genome is not explored. In this study, we provide the first overview of the selection signatures in the Vrindavani genome. A total of 96 Vrindavani cattle were genotyped using the BovineSNP50 BeadChip and the SNP genotype data of its constituent breeds were collected from a public database. Within-breed selection signatures in Vrindavani were investigated using the integrated haplotype score (iHS). The Vrindavani breed was also compared to each of its parental breeds to discover between-population signatures of selection using two approaches, cross-population extended haplotype homozygosity (XP-EHH) and fixation index (F_ST). We identified 11 common regions detected by more than one method harboring genes such as LRP1B, TNNI3K, APOB, CACNA2D1, FAM110B, and SPATA17 associated with production and adaptation. Overall, our results suggested stronger selective pressure on regions responsible for adaptation compared to milk yield.

Genome-wide scan for selection signatures reveals novel insights into the adaptive capacity in local North African cattle

Natural-driven selection is supposed to have left detectable signatures on the genome of North African cattle which are often characterized by the fixation of genetic variants associated with traits under selection pressure and/or an outstanding genetic differentiation with other populations at particular loci. Here, we investigate the population genetic structure and we provide a first outline of potential selection signatures in North African cattle using single nucleotide polymorphism genotyping data. After comparing our data to African, European and indicine cattle populations, we identified 36 genomic regions using three extended haplotype homozygosity statistics and 92 outlier markers based on Bayescan test. The 13 outlier windows detected by at least two approaches, harboured genes (e.g. GH1, ACE, ASIC3, HSPH1, MVD, BCL2, HIGD2A, CBFA2T3) that may be involved in physiological adaptations required to cope with environmental stressors that are typical of the North African area such as infectious diseases, extended drought periods, scarce food supply, oxygen scarcity in the mountainous areas and high-intensity solar radiation. Our data also point to candidate genes involved in transcriptional regulation suggesting that regulatory elements had also a prominent role in North African cattle response to environmental constraints. Our study yields novel insights into the unique adaptive capacity in these endangered populations emphasizing the need for the use of whole genome sequence data to gain a better understanding of the underlying molecular mechanisms.

The mosaic genome of indigenous African cattle as a unique genetic resource for African pastoralism

Cattle pastoralism plays a central role in human livelihood in Africa. However, the genetic history of its success remains unknown. Here, through whole-genome sequence analysis of 172 indigenous African cattle from 16 breeds representative of the main cattle groups, we identify a major taurine × indicine cattle admixture event dated to circa 750-1,050 yr ago, which has shaped the genome of today's cattle in the Horn of Africa. We identify 16 loci linked to African environmental adaptations across crossbred animals showing an excess of taurine or indicine ancestry. These include immune-, heat-tolerance- and reproduction-related genes. Moreover, we identify one highly divergent locus in African taurine cattle, which is putatively linked to trypanotolerance and present in crossbred cattle living in trypanosomosis-infested areas. Our findings indicate that a combination of past taurine and recent indicine admixture-derived genetic resources is at the root of the present success of African pastoralism.

SNP Diversity in CD14 Gene Promoter Suggests Adaptation Footprints in Trypanosome Tolerant N'Dama (Bos taurus) but not in Susceptible White Fulani (Bos indicus) Cattle

Immune response to infections has been shown to be mediated by genetic diversity in pattern recognition receptors, leading to disease tolerance or susceptibility. We elucidated naturally occurring variations within the bovine CD14 gene promoter in trypanosome-tolerant (N'Dama) and susceptible (White Fulani) cattle, with genomic and computational approaches. Blood samples were collected from White Fulani and N'Dama cattle, genomic DNA extracted and the entire promoter region of the CD14 gene amplified by PCR. We sequenced this region and performed in silico computation to identify SNP variants, transcription factor binding sites, as well as micro RNAs in the region. CD14 promoter sequences were compared with the reference bovine genome from the Ensembl database to identify various SNPs. Furthermore, we validated three selected N'Dama specific SNPs using custom Taqman SNP genotyping assay for genetic diversity. In all, we identified a total of 54 and 41 SNPs at the CD14 promoter for N'Dama and White Fulani respectively, including 13 unique SNPs present in N'Dama only. The significantly higher SNP density at the CD14 gene promoter region in N'Dama may be responsible for disease tolerance, possibly an evolutionary adaptation. Our genotype analysis of the three loci selected for validation show that mutant alleles (A/A, C/C, and A/A) were adaptation profiles within disease tolerant N'Dama. A similar observation was made for our haplotype analysis revealing that haplotypes H1 (ACA) and H2 (ACG) were significant combinations within the population. The SNP effect prediction revealed 101 and 89 new transcription factor binding sites in N'Dama and White Fulani, respectively. We conclude that disease tolerant N'Dama possessing higher SNP density at the CD14 gene promoter and the preponderance of mutant alleles potentially confirms the significance of this promoter in immune response, which is lacking in susceptible White Fulani. We, therefore, recommend further in vitro and in vivo study of this observation in infected animals, as the next step for understanding genetic diversity relating to varying disease phenotypes in both breeds.