Whole-Genome Resequencing Reveals Selection Signatures of Abigar Cattle for Local Adaptation.

Over time, indigenous cattle breeds have developed disease resistance, heat tolerance, and adaptability to harsh environments. Deciphering the genetic mechanisms underlying adaptive traits is crucial for their improvement and sustainable utilization. For the first time, we performed whole-genome sequencing to unveil the genomic diversity, population structure, and selection signatures of Abigar cattle living in a tropical environment. The population structure analysis revealed that Abigar cattle exhibit high nucleotide diversity and heterozygosity, with low runs of homozygosity and linkage disequilibrium, suggesting a genetic landscape less constrained by inbreeding and enriched by diversity. Using nucleotide diversity (Pi) and population differentiation (FST) selection scan methods, we identified 83 shared genes that are likely associated with tropical adaption. The functional annotation analysis revealed that some of these genes are potentially linked to heat tolerance (HOXC13, DNAJC18, and RXFP2), immune response (IRAK3, MZB1, and STING1), and oxidative stress response (SLC23A1). Given the wider spreading impacts of climate change on cattle production, understanding the genetic mechanisms of adaptation of local breeds becomes crucial to better respond to climate and environmental changes. In this context, our finding establishes a foundation for further research into the mechanisms underpinning cattle adaptation to tropical environments.

Readiness and early response to COVID-19: achievements, challenges and lessons learnt in Ethiopia.

Declaration of the novel coronavirus disease as a Public Health Emergency of International Concern necessitated countries to get ready to respond. Here, we describe key achievements, challenges and lessons learnt during the readiness and early response to COVID-19 in Ethiopia. Readiness activities commenced as early as January 2020 with the activation of a national Public Health Emergency Operations Centre and COVID-19 Incident Management System (IMS) by the Ethiopian Public Health Institute. The COVID-19 IMS conducted rapid risk assessments, developed scenario-based contingency plans, national COVID-19 guidelines and facilitated the enhancement of early warning and monitoring mechanisms. Early activation of a coordination mechanism and strengthening of detection and response capacities contributed to getting the country ready on time and mounting an effective early response. High-level political leadership and commitment led to focused efforts in coordination of response interventions. Health screening, mandatory 14-day quarantine and testing established for all international travellers arriving into the country slowed down the influx of travellers. The International Health Regulations (IHR) capacities in the country served as a good foundation for timely readiness and response. Leveraging on existing IHR capacities in the country built prior to COVID-19 helped slow down the importation and mitigated uncontrolled spread of the disease in the country. Challenges experienced included late operationalisation of a multisectoral coordination platform, shortage of personal protective equipment resulting from global disruption of importation and the huge influx of over 10 000 returnees from different COVID-19-affected countries over a short period of time with resultant constrain on response resources.

Genetic assessment of BoLA-DRB3 polymorphisms by comparing Bangladesh, Ethiopian, and Korean cattle.

Attributable to their major function in pathogen recognition, the use of bovine leukocyte antigens (BoLA) as disease markers in immunological traits in cattle is well established. However, limited report exists on polymorphism of the BoLA gene in zebu cattle breeds by high resolution typing methods. Thus, we used a polymerase chain reaction sequence-based typing (PCR-SBT) method to sequence exon 2 of the BoLA class II DRB3 gene from 100 animals (Boran, n = 13; Sheko, n = 20; Fogera, n = 16; Horro, n = 19), Hanwoo cattle (n = 18) and Bangladesh Red Chittagong zebu (n = 14). Out of the 59 detected alleles, 43 were already deposited under the Immuno Polymorphism Database for major histocompatibility complex (IPD-MHC) while 16 were unique to this study. Assessment of the level of genetic variability at the population and sequence levels with genetic distance in the breeds considered in this study showed that Zebu breeds had a gene diversity score greater than 0.752, nucleotide diversity score greater than 0.152, and mean number of pairwise differences higher than 14, being very comparable to those investigated for other cattle breeds. Regarding neutrality tests analyzed, we investigated that all the breeds except Hanwoo had an excess number of alleles and could be expected from a recent population expansion or genetic hitchhiking. Howbeit, the observed heterozygosity was not significantly (p < 0.05) higher than the expected heterozygosity. The Hardy Weinberg equilibrium (HWE) analysis revealed non-significant excess of heterozygote animals, indicative of plausible over-dominant selection. The pairwise FST values suggested a low genetic variation among all the breeds (FST = 0.056; p < 0.05), besides the rooting from the evolutionary or domestication history of the cattle. No detached clade was observed in the evolutionary divergence study of the BoLA-DRB3 gene, inferred from the phylogenetic tree based on the maximum likelihood model. The investigation herein indicated the clear differences in BoLA-DRB3 gene variability between African and Asian cattle breeds.

Single nucleotide polymorphisms of leptin gene in five Ethiopian indigenous cattle breeds and the Korean Hanwoo breed.

Considering the escalating number of scientific reports on the association between the leptin gene and diverse physiological traits and performance of cattle populations, this study was directed towards identifying SNPs in the leptin gene among five indigenous cattle breeds of Ethiopia. DNA samples were extracted from the nasal swabs of the Ethiopian indigenous cattle breeds: Arsi (n = 18), Horro (n = 20), Begait (n = 21), Boran (n = 19), and Fogera (n = 17) and the Korean Hanwoo (a representative taurine breed) (n = 20), followed by PCR amplification of exon 2 and exon 3 regions of the leptin gene and sequence analysis of the PCR products. Five SNPs, two (generating missense mutations) on exon 2 and three (generating silent mutations) on exon 3 regions, were explicated in this study. Allele frequency and genotype frequency distribution pertaining to the SNPs were recorded for the studied cattle breeds besides the minor allele frequency and deviation from the Hardy-Weinberg equilibrium. Positive FIS index values were recorded for all the markers except SNP2, illustrative of heterozygote deficiency. MEGA X software-based evolutionary divergence analysis of the phylogenetic tree based on the SNP data revealed that the large-sized breeds, Hanwoo, Begait, Boran, and Fogera, were more closely clustered compared to the small-sized Arsi breed. Among the seven haplotypes documented from the various breeds, sequence analysis was suggestive of haplotypes 1 and 2 to be ancestral haplotypes for the leptin gene. This study is envisaged to accelerate molecular breeding programs for the genetic improvement of the Ethiopian cattle breeds.

Genome-wide insights of Ethiopian indigenous sheep populations reveal the population structure related to tail morphology and phylogeography.

BACKGROUND: Ethiopian sheep living in different climatic zones and having contrasting morphologies are a most promising subject of molecular-genetic research. Elucidating their genetic diversity and genetic structure is critical for designing appropriate breeding and conservation strategies. OBJECTIVE: The study was aimed to investigate genome-wide genetic diversity and population structure of eight Ethiopian sheep populations. METHODS: A total of 115 blood samples were collected from four Ethiopian sheep populations that include Washera, Farta and Wollo (short fat-tailed) and Horro (long fat-tailed). DNA was extracted using Quick-DNA™ Miniprep plus kit. All DNA samples were genotyped using Ovine 50 K SNP BeadChip. To infer genetic relationships of Ethiopian sheep at national, continental and global levels, genotype data on four Ethiopian sheep (Adilo, Arsi-Bale, Menz and Black Head Somali) and sheep from east, north, and south Africa, Middle East and Asia were included in the study as reference. RESULTS: Mean genetic diversity of Ethiopian sheep populations ranged from 0.352 ± 0.14 for Horro to 0.379 ± 0.14 for Arsi-Bale sheep. Population structure and principal component analyses of the eight Ethiopian indigenous sheep revealed four distinct genetic cluster groups according to their tail phenotype and geographical distribution. The short fat-tailed sheep did not represent one genetic cluster group. Ethiopian fat-rump sheep share a common genetic background with the Kenyan fat-tailed sheep. CONCLUSION: The results of the present study revealed the principal component and population structure follows a clear pattern of tail morphology and phylogeography. There is clear signature of admixture among the study Ethiopian sheep populations.

Signatures of positive selection underlying beef production traits in Korean cattle breeds.

The difference in the breeding programs and population history may have diversely shaped the genomes of Korean native cattle breeds. In the absence of phenotypic data, comparisons of breeds that have been subjected to different selective pressures can aid to identify genomic regions and genes controlling qualitative and complex traits. In this study to decipher genetic variation and identify evidence of divergent selection, 3 Korean cattle breeds were genotyped using the recently developed high-density GeneSeek Genomic Profiler F250 (GGP-F250) array. The three Korean cattle breeds clustered according to their coat color phenotypes and breeding programs. The Heugu breed reliably showed smaller effective population size at all generations considered. Across the autosomal chromosomes, 113 and 83 annotated genes were identified from Hanwoo-Chikso and Hanwoo-Heugu comparisons, respectively of which 16 genes were shared between the two pairwise comparisons. The most important signals of selection were detected on bovine chromosomes 14 (24.39-25.13 Mb) and 18 (13.34-15.07 Mb), containing genes related to body size, and coat color (XKR4, LYN, PLAG1, SDR16C5, TMEM68, CDH15, MC1R, and GALNS). Some of the candidate genes are also associated with meat quality traits (ACSF3, EIF2B1, BANP, APCDD1, and GALM) and harbor quantitative trait locus (QTL) for beef production traits. Further functional analysis revealed that the candidate genes (DBI, ACSF3, HINT2, GBA2, AGPAT5, SCAP, ELP6, APOB, and RBL1) were involved in gene ontology (GO) terms relevant to meat quality including fatty acid oxidation, biosynthesis, and lipid storage. Candidate genes previously known to affect beef production and quality traits could be used in the beef cattle selection strategies.

Meningococcal Meningitis Outbreaks in the African Meningitis Belt After Meningococcal Serogroup A Conjugate Vaccine Introduction, 2011-2017.

BACKGROUND: In 2010-2017, meningococcal serogroup A conjugate vaccine (MACV) was introduced in 21 African meningitis belt countries. Neisseria meningitidis A epidemics have been eliminated here; however, non-A serogroup epidemics continue. METHODS: We reviewed epidemiological and laboratory World Health Organization data after MACV introduction in 20 countries. Information from the International Coordinating Group documented reactive vaccination. RESULTS: In 2011-2017, 17 outbreaks were reported (31 786 suspected cases from 8 countries, 1-6 outbreaks/year). Outbreaks were of 18-14 542 cases in 113 districts (median 3 districts/outbreak). The most affected countries were Nigeria (17 375 cases) and Niger (9343 cases). Cumulative average attack rates per outbreak were 37-203 cases/100 000 population (median 112). Serogroup C accounted for 11 outbreaks and W for 6. The median proportion of laboratory confirmed cases was 20%. Reactive vaccination was conducted during 14 outbreaks (5.7 million people vaccinated, median response time 36 days). CONCLUSION: Outbreaks due to non-A serogroup meningococci continue to be a significant burden in this region. Until an affordable multivalent conjugate vaccine becomes available, the need for timely reactive vaccination and an emergency vaccine stockpile remains high. Countries must continue to strengthen detection, confirmation, and timeliness of outbreak control measures.

Genomic signatures of high-altitude adaptation in Ethiopian sheep populations.

BACKGROUND: Ethiopian sheep populations such as Arsi-Bale, Horro and Adilo (long fat-tailed, LFT) inhabit mid to high-altitude areas; and Menz sheep (MZ, short fat-tailed) are adapted to cool sub-alpine environments. In contrast, Blackhead Somali sheep (BHS, fat-rumped) thrive well in arid and semi-arid areas characterized by high temperature and low precipitation. The genomic investigation of Ethiopian sheep populations may help to identify genes and biological pathways enable to adapt to the different ecological conditions. OBJECTIVE: To uncover genomic regions and genes showing evidence of positive selection for altitude adaptation in Ethiopian sheep populations. METHODS: A total of 72 animals inhabiting high-versus low-altitude environments were genotyped on an Ovine Infinium HD array (~ 600 K). Pairwise genetic differentiation (Fst) was calculated in sliding windows of 20 SNPs and the upper 1% smoothed Fst values were considered to represent positive selection signatures. Genes within < 25 kb of the most differentiated SNPs were considered as selection candidates. RESULTS: Signatures of selection were detected in genes known to be associated high with altitude adaptation in MZ-BHS pair comparison (PPP1R12A, RELN, PARP2, and DNAH9) and in LFT-BHS pair comparison (VAV3, MSRB3,EIF2AK4, MET, and TACR1). The candidate genes (MITF, FGF5, MTOR, TRHDE, and TUBB3) associated with altitude adaptation and shared between the MZ-BHS and LTF-BHS pair comparisons were also detected as under selection. Further functional analyses reveal that the candidate genes were involved in biological processes and pathways relevant to adaptation under extreme altitudes, including respiratory system development and smoothened signaling pathway. CONCLUSION: The results of the present study could aid in-depth understanding and exploitation of the underlying genetic mechanisms for sheep and other livestock species adaptation to high-altitude environments.

Genome-wide diversity and demographic dynamics of Cameroon goats and their divergence from east African, north African, and Asian conspecifics.

Indigenous goats make significant contributions to Cameroon's national and local economy, but little effort has been devoted to identifying the populations. Here, we assessed the genetic diversity and demographic dynamics of Cameroon goat populations using mitochondrial DNA (two populations) and autosomal markers (four populations) generated with the Caprine 50K SNP chip. To infer genetic relationships at continental and global level, genotype data on six goat populations from Ethiopia and one population each from Egypt, Morocco, Iran, and China were included in the analysis. The mtDNA analysis revealed 83 haplotypes, all belonging to haplogroup A, in Cameroon goats. Four haplotypes were shared between goats found in Cameroon, Mozambique, Namibia, Zimbabwe, Kenya, and Ethiopia. Analysis of autosomal SNPs in Cameroon goats revealed the lowest HO (0.335±0.13) and HE (0.352±0.15) in the North-west Highland and Central Highland populations, respectively. Overall, the highest HO (0.401±0.12) and HE (0.422±0.12) were found for Barki and Iranian goats, respectively. Barki goats had the highest average MAF, while Central Highland Cameroon goats had the lowest. Overall, Cameroon goats demonstrated high FIS. AMOVA revealed that 13.29% of the variation was explained by genetic differences between the six population groups. Low average FST (0.01) suggests intermixing among Cameroon goats. All measures indicated that Cameroon goats are closer to Moroccan goats than to other goat populations. PCA and STRUCTURE analyses poorly differentiated the Cameroon goats, as did genetic distance, Neighbor-Net network, and neighbor-joining tree analyses. The haplotype analysis of mtDNA showed the initial dispersion of goats to Cameroon and central Africa from north-east Africa following the Nile Delta. Whereas, the approximate Bayesian computation indicated Cameroon goats were separated from Moroccan goats after 506 generations in later times (~1518 YA), as supported by the phylogenetic net-work and admixture outputs. Overall, indigenous goats in Cameroon show weak phylogenetic structure, suggesting either extensive intermixing.

Genome-wide association study of carcass weight in commercial Hanwoo cattle.

OBJECTIVE: The objective of the present study was to validate genes and genomic regions associated with carcass weight using a low-density single nucleotide polymorphism (SNP) Chip in Hanwoo cattle breed. METHODS: Commercial Hanwoo steers (n = 220) were genotyped with 20K GeneSeek genomic profiler BeadChip. After applying the quality control of criteria of a call rate ≥90% and minor allele frequency (MAF) ≥0.01, a total of 15,235 autosomal SNPs were left for genome-wide association (GWA) analysis. The GWA tests were performed using single-locus mixed linear model. Age at slaughter was fitted as fixed effect and sire included as a covariate. The level of genome-wide significance was set at 3.28×10-6 (0.05/15,235), corresponding to Bonferroni correction for 15,235 multiple independent tests. RESULTS: By employing EMMAX approach which is based on a mixed linear model and accounts for population stratification and relatedness, we identified 17 and 16 loci significantly (p<0.001) associated with carcass weight for the additive and dominant models, respectively. The second most significant (p = 0.000049) SNP (ARS-BFGL-NGS-28234) on bovine chromosome 4 (BTA4) at 21 Mb had an allele substitution effect of 43.45 kg. Some of the identified regions on BTA2, 6, 14, 22, and 24 were previously reported to be associated with quantitative trait loci for carcass weight in several beef cattle breeds. CONCLUSION: This is the first genome-wide association study using SNP chips on commercial Hanwoo steers, and some of the loci newly identified in this study may help to better DNA markers that determine increased beef production in commercial Hanwoo cattle. Further studies using a larger sample size will allow confirmation of the candidates identified in this study.

Genetic Diversity and Population Structure of Ethiopian Sheep Populations Revealed by High-Density SNP Markers.

Sheep in Ethiopia are adapted to a wide range of environments, including extreme habitats. Elucidating their genetic diversity is critical for improving breeding strategies and mapping quantitative trait loci associated with productivity. To this end, the present study investigated the genetic diversity and population structure of five Ethiopian sheep populations exhibiting distinct phenotypes and sampled from distinct production environments, including arid lowlands and highlands. To investigate the genetic relationships in greater detail and infer population structure of Ethiopian sheep breeds at the continental and global levels, we analyzed genotypic data of selected sheep breeds from the Ovine SNP50K HapMap dataset. All Ethiopian sheep samples were genotyped with Ovine Infinium HD SNP BeadChip (600K). Mean genetic diversity ranged from 0.29 in Arsi-Bale to 0.32 in Menz sheep, while estimates of genetic differentiation among populations ranged from 0.02 to 0.07, indicating low to moderate differentiation. An analysis of molecular variance revealed that 94.62 and 5.38% of the genetic variation was attributable to differences within and among populations, respectively. Our population structure analysis revealed clustering of five Ethiopian sheep populations according to tail phenotype and geographic origin-i.e., short fat-tailed (very cool high-altitude), long fat-tailed (mid to high-altitude), and fat-rumped (arid low-altitude), with clear evidence of admixture between long fat-tailed populations. North African sheep breeds showed higher levels of within-breed diversity, but were less differentiated than breeds from Eastern and Southern Africa. When African breeds were grouped according to geographic origin (North, South, and East), statistically significant differences were detected among groups (regions). A comparison of population structure between Ethiopian and global sheep breeds showed that fat-tailed breeds from Eastern and Southern Africa clustered together, suggesting that these breeds were introduced to the African continent via the Horn and migrated further south.

Genome-wide Single Nucleotide Polymorphism Analyses Reveal Genetic Diversity and Structure of Wild and Domestic Cattle in Bangladesh.

In spite of variation in coat color, size, and production traits among indigenous Bangladeshi cattle populations, genetic differences among most of the populations have not been investigated or exploited. In this study, we used a high-density bovine single nucleotide polymorphism (SNP) 80K Bead Chip derived from Bos indicus breeds to assess genetic diversity and population structure of 2 Bangladeshi zebu cattle populations (red Chittagong, n = 28 and non-descript deshi, n = 28) and a semi-domesticated population (gayal, n = 17). Overall, 95% and 58% of the total SNPs (69,804) showed polymorphisms in the zebu and gayal populations, respectively. Similarly, the average minor allele frequency value was as high 0.29 in zebu and as low as 0.09 in gayal. The mean expected heterozygosity varied from 0.42±0.14 in zebu to 0.148±0.14 in gayal with significant heterozygosity deficiency of 0.06 (F IS) in the latter. Coancestry estimations revealed that the two zebu populations are weakly differentiated, with over 99% of the total genetic variation retained within populations and less than 1% accounted for between populations. Conversely, strong genetic differentiation (FST = 0.33) was observed between zebu and gayal populations. Results of population structure and principal component analyses suggest that gayal is distinct from Bos indicus and that the two zebu populations were weakly structured. This study provides basic information about the genetic diversity and structure of Bangladeshi cattle and the semi-domesticated gayal population that can be used for future appraisal of breed utilization and management strategies.

Genetic Structure of and Evidence for Admixture between Western and Korean Native Pig Breeds Revealed by Single Nucleotide Polymorphisms.

Comprehensive information on genetic diversity and introgression is desirable for the design of rational breed improvement and conservation programs. Despite the concerns regarding the genetic introgression of Western pig breeds into the gene pool of the Korean native pig (KNP), the level of this admixture has not yet been quantified. In the present study, we genotyped 93 animals, representing four Western pig breeds and KNP, using the porcine SNP 60K BeadChip to assess their genetic diversity and to estimate the level of admixture among the breeds. Expected heterozygosity was the lowest in Berkshire (0.31) and highest in Landrace (0.42). Population differentiation (FST) estimates were significantly different (p<0.000), accounting for 27% of the variability among the breeds. The evidence of inbreeding observed in KNP (0.029) and Yorkshire (0.031) may result in deficient heterozygosity. Principal components one (PC1) and two (PC2) explained approximately 35.06% and 25.20% of the variation, respectively, and placed KNP somewhat proximal to the Western pig breeds (Berkshire and Landrace). When K = 2, KNP shared a substantial proportion of ancestry with Western breeds. Similarly, when K = 3, over 86% of the KNP individuals were in the same cluster with Berkshire and Landrace. The linkage disquilbrium (LD) values at r (2) 0.3, the physical distance at which LD decays below a threshold of 0.3, ranged from 72.40 kb in Landrace to 85.86 kb in Yorkshire. Based on our structure analysis, a substantial level of admixture between Western and Korean native pig breeds was observed.

Semi-domesticated and Irreplaceable Genetic Resource Gayal (Bos frontalis) Needs Effective Genetic Conservation in Bangladesh: A Review.

Several studies arduously reported that gayal (Bos frontalis) is an independent bovine species. The population size is shrinking across its distribution. In Bangladesh, it is the only wild relative of domestic cattle and also a less cared animal. Their body size is much bigger than Bangladeshi native cattle and has prominent beef type characters along with the ability to adjust in any adverse environmental conditions. Human interactions and manipulation of biodiversity is affecting the habitats of gayals in recent decades. Besides, the only artificial reproduction center for gayals, Bangladesh Livestock Research Institute (BLRI), has few animals and could not carry out its long term conservation scheme due to a lack of an objective based scientific mission as well as financial support. This indicates that the current population is much more susceptible to stochastic events which might be natural catastrophes, environmental changes or mutations. Further reduction of the population size will sharply reduce genetic diversity. In our recent investigation with 80K indicine single nucleotide polymorphism chip, the F IS (within-population inbreeding) value was reported as 0.061±0.229 and the observed (0.153±0.139) and expected (0.148±0.143) heterozygosities indicated a highly inbred and less diverse gayal population in Bangladesh. Prompt action is needed to tape the genetic information of this semi-domesticated bovine species with considerable sample size and try to investigate its potentials together with native zebu cattle for understanding the large phenotypic variations, improvement and conservation of this valuable creature.

A whole genomic scan to detect selection signatures between Berkshire and Korean native pig breeds.

BACKGROUND: Scanning of the genome for selection signatures between breeds may play important role in understanding the underlie causes for observable phenotypic variations. The discovery of high density single nucleotide polymorphisms (SNPs) provide a useful starting point to perform genome-wide scan in pig populations in order to identify loci/candidate genes underlie phenotypic variation in pig breeds and facilitate genetic improvement programs. However, prior to this study genomic region under selection in commercially selected Berkshire and Korean native pig breeds has never been detected using high density SNP markers. To this end, we have genotyped 45 animals using Porcine SNP60 chip to detect selection signatures in the genome of the two breeds by using the F ST approach. RESULTS: In the comparison of Berkshire and KNP breeds using the FDIST approach, a total of 1108 outlier loci (3.48%) were significantly different from zero at 99% confidence level with 870 of the outlier SNPs displaying high level of genetic differentiation (F ST ≥0.490). The identified candidate genes were involved in a wide array of biological processes and molecular functions. Results revealed that 19 candidate genes were enriched in phosphate metabolism (GO: 0006796; ADCK1, ACYP1, CAMK2D, CDK13, CDK13, ERN1, GALK2, INPP1; MAK, MAP2K5, MAP3K1, MAPK14, P14KB, PIK3C3, PRKC1, PTPRK, RNASEL, THBS1, BRAF, VRK1). We have identified a set of candidate genes under selection and have known to be involved in growth, size and pork quality (CART, AGL, CF7L2, MAP2K5, DLK1, GLI3, CA3 and MC3R), ear morphology and size (HMGA2 and SOX5) stress response (ATF2, MSRB3, TMTC3 and SCAF8) and immune response ( HCST and RYR1). CONCLUSIONS: Some of the genes may be used to facilitate genetic improvement programs. Our results also provide insights for better understanding of the process and influence of breed development on the pattern of genetic variations.

Genetic diversity, population structure and relationships in indigenous cattle populations of Ethiopia and Korean Hanwoo breeds using SNP markers.

In total, 166 individuals from five indigenous Ethiopian cattle populations - Ambo (n = 27), Borana (n = 35), Arsi (n = 30), Horro (n = 36), and Danakil (n = 38) - were genotyped for 8773 single nucleotide polymorphism (SNP) markers to assess genetic diversity, population structure, and relationships. As a representative of taurine breeds, Hanwoo cattle (n = 40) were also included in the study for reference. Among Ethiopian cattle populations, the proportion of SNPs with minor allele frequencies (MAFs) ≥0.05 ranged from 81.63% in Borana to 85.30% in Ambo, with a mean of 83.96% across all populations. The Hanwoo breed showed the highest proportion of polymorphism, with MAFs ≥0.05, accounting for 95.21% of total SNPs. The mean expected heterozygosity varied from 0.370 in Danakil to 0.410 in Hanwoo. The mean genetic differentiation (F ST; 1%) in Ethiopian cattle revealed that within individual variation accounted for approximately 99% of the total genetic variation. As expected, F ST and Reynold genetic distance were greatest between Hanwoo and Ethiopian cattle populations, with average values of 17.62 and 18.50, respectively. The first and second principal components explained approximately 78.33% of the total variation and supported the clustering of the populations according to their historical origins. At K = 2 and 3, a considerable source of variation among cattle is the clustering of the populations into Hanwoo (taurine) and Ethiopian cattle populations. The low estimate of genetic differentiation (F ST) among Ethiopian cattle populations indicated that differentiation among these populations is low, possibly owing to a common historical origin and high gene flow. Genetic distance, phylogenic tree, principal component analysis, and population structure analyses clearly differentiated the cattle population according to their historical origins, and confirmed that Ethiopian cattle populations are genetically distinct from the Hanwoo breed.

Comparison of SNP variation and distribution in indigenous ethiopian and korean cattle (hanwoo) populations.

Although a large number of single nucleotide polymorphisms (SNPs) have been identified from the bovine genome-sequencing project, few of these have been validated at large in Bos indicus breeds. We have genotyped 192 animals, representing 5 cattle populations of Ethiopia, with the Illumina Bovine 8K SNP BeadChip. These include 1 Sanga (Danakil), 3 zebu (Borana, Arsi and Ambo), and 1 zebu × Sanga intermediate (Horro) breeds. The Hanwoo (Bos taurus) was included for comparison purposes. Analysis of 7,045 SNP markers revealed that the mean minor allele frequency (MAF) was 0.23, 0.22, 0.21, 0.21, 0.23, and 0.29 for Ambo, Arsi, Borana, Danakil, Horro, and Hanwoo, respectively. Significant differences of MAF were observed between the indigenous Ethiopian cattle populations and Hanwoo breed (p < 0.001). Across the Ethiopian cattle populations, a common variant MAF (≥0.10 and ≤0.5) accounted for an overall estimated 73.79% of the 7,045 SNPs. The Hanwoo displayed a higher proportion of common variant SNPs (90%). Investigation within Ethiopian cattle populations showed that on average, 16.64% of the markers were monomorphic, but in the Hanwoo breed, only 6% of the markers were monomorphic. Across the sampled Ethiopian cattle populations, the mean observed and expected heterozygosities were 0.314 and 0.313, respectively. The level of SNP variation identified in this particular study highlights that these markers can be potentially used for genetic studies in African cattle breeds.