Genetic diversity and population structure among six cattle breeds in South Africa using a whole genome SNP panel

Investigation of cosmopolitan and local Italian beef cattle breeds uncover common patterns of heterozygosity

The analysis of livestock heterozygosity is less common compared to the study of homozygous patterns. Heterozygous-Rich Regions (HRRs) may harbor significant loci for functional traits such as immune response, survival rate, and fertility. For this reason, this study was conducted to investigate and characterize the heterozygosity patterns of four beef cattle breeds, which included two cosmopolitan breeds (Limousine and Charolaise) and two local breeds (Sarda and Sardo Bruna). Our analysis identified regions with a high degree of heterozygosity using a consecutive runs approach, the Tajima D test, nucleotide diversity estimation, and Hardy Weinberg equilibrium test. These regions exhibited recurrent heterozygosity peaks and were consistently found on specific chromosomes across all breeds, specifically autosomes 15, 16, 20, and 23. The cosmopolitan and Sardo Bruna breeds also displayed peaks on autosomes 2 and 21, respectively. Thirty-five top runs shared by more than 25% of the populations were identified. These genomic fragments encompassed 18 genes, two of which are directly linked to male fertility, while four are associated with lactation. Two other genes play roles in survival and immune response. Our study also detected a region related to growth and carcass traits in Limousine breed. Our analysis of heterozygosity-rich regions revealed particular segments of the cattle genome linked to various functional traits. It appears that balancing selection is occurring in specific regions within the four examined breeds, and unexpectedly, they are common across cosmopolitan and local breeds. The genes identified hold potential for applications in breeding programs and conservation studies to investigate the phenotypes associated with these heterozygous genotypes. In addition, Tajima D test, Nucleotide diversity, and Hardy Weinberg equilibrium test confirmed the presence of heterozygous fragments found with Heterozygous-Rich Regions analysis.

Genome-wide mining of diversity and evolutionary signatures revealed selective hotspots in Indian Sahiwal cattle

The Sahiwal cattle breed is the best indigenous dairy cattle breed, and it plays a pivotal role in the Indian dairy industry. This is due to its exceptional milk-producing potential, adaptability to local tropical conditions, and its resilience to ticks and diseases. The study aimed to identify selective sweeps and estimate intrapopulation genetic diversity parameters in Sahiwal cattle using ddRAD sequencing-based genotyping data from 82 individuals. After applying filtering criteria, 78,193 high-quality SNPs remained for further analysis. The population exhibited an average minor allele frequency of 0.221 ± 0.119. Genetic diversity metrics, including observed (0.597 ± 0.196) and expected heterozygosity (0.433 ± 0.096), nucleotide diversity (0.327 ± 0.114), the proportion of polymorphic SNPs (0.726), and allelic richness (1.323 ± 0.134), indicated ample genomic diversity within the breed. Furthermore, an effective population size of 74 was observed in the most recent generation. The overall mean linkage disequilibrium (r2) for pairwise SNPs was 0.269 ± 0.057. Moreover, a greater proportion of short Runs of Homozygosity (ROH) segments were observed suggesting that there may be low levels of recent inbreeding in this population. The genomic inbreeding coefficients, computed using different inbreeding estimates (FHOM, FUNI, FROH, and FGROM), ranged from -0.0289 to 0.0725. Subsequently, we found 146 regions undergoing selective sweeps using five distinct statistical tests: Tajima's D, CLR, |iHS|, |iHH12|, and ROH. These regions, located in non-overlapping 500 kb windows, were mapped and revealed various protein-coding genes associated with enhanced immune systems and disease resistance (IFNL3, IRF8, BLK), as well as production traits (NRXN1, PLCE1, GHR). Notably, we identified interleukin 2 (IL2) on Chr17: 35217075-35223276 as a gene linked to tick resistance and uncovered a cluster of genes (HSPA8, UBASH3B, ADAMTS18, CRTAM) associated with heat stress. These findings indicate the evolutionary impact of natural and artificial selection on the environmental adaptation of the Sahiwal cattle population.

Genetic diversity, population structure and origin of the native goats in Central Laos

Maintaining genetic diversity and variation in livestock populations is critical for natural and artificial selection promoting genetic improvement while avoiding problems due to inbreeding. In Laos, there are concerns that there has been a decline in genetic diversity and a rise in inbreeding among native goats in their village-based smallholder system. In this study, we investigated the genetic diversity of Lao native goats in Phin, Songkhone and Sepon districts in Central Laos for the first time using Illumina's Goat SNP50 BeadChip. We also explored the genetic relationships between Lao goats with 163 global goat populations from 36 countries. Our results revealled a close genetic relationship between Lao native goats and Chinese, Mongolian and Pakistani goats, sharing ancestries with Guangfen, Jining Grey and Luoping Yellow breeds (China) and Teddi goats (Pakistan). The observed (Ho) and expected (He) heterozygosity were 0.292 and 0.303 (Laos), 0.288 and 0.288 (Sepon), 0.299 and 0.308 (Phin) and 0.289 and 0.305 (Songkhone), respectively. There was low to moderate genetic differentiation (FST: 0.011-0.043) and negligible inbreeding coefficients (FIS: -0.001 to 0.052) between goat districts. The runs of homozygosity (ROH) had an average length of 5.92-6.85 Mb, with short ROH segments (1-5 Mb length) being the most prevalent (66.34%). Longer ROH segments (20-40 and >40 Mb length categories) were less common, comprising only 4.81% and 1.01%, respectively. Lao goats exhibit moderate genetic diversity, low-inbreeding levels and adequate effective population size. Some genetic distinctions between Lao goats may be explained by geographic and cultural features.

Analysing the genetic diversity and population structure of five native Turkish cattle breeds using SNP data

The conservation and sustainable utilization of cattle genetic resources necessitate a comprehensive understanding of their genetic diversity and population structure. This study provides an analysis of five native Turkish cattle breeds: Anatolian Black (ANB), Turkish Grey (TUR), Anatolian Southern Yellow (ASY), East Anatolian Red (EAR), and South Anatolian Red (SAN) using 50 K SNP data. These breeds were compared with three European breeds, Simmental (SIM), Holstein (HOL), and Jersey (JER), and three Asian Zebu breeds: Arabic Zebu (ZAR), Nelore (NEL), and Red Sindhi (RSI). Genetic diversity indices demonstrated moderate heterogeneity among the breeds, with TUR exhibiting the highest observed heterozygosity (Ho = 0.35). Wright's Fst values indicated significant genetic differentiation, particularly between Turkish breeds and both European (Fst = 0.035-0.071) and Asian breeds (Fst = 0.025-0.150). Principal component analysis distinguished the unique genetic profiles of each breed cluster. Admixture analysis revealed degrees of shared genetic ancestry, suggesting historical gene flow between Turkish, European, and Asian breeds. Analysis of molecular variance (AMOVA) attributed approximately 58% of the variation to population differences. Nei's genetic distances highlighted the closer genetic relatedness within Turkish breeds (distance ranges between 0.032 and 0.046) and suggested a more relative affinity of TUR with European breeds. The study's phylogenetic assessments elucidate the nuanced genetic relationships among these breeds, with runs of homozygosity (ROH) analysis indicating patterns of ancestral relatedness and moderate levels of inbreeding, particularly evident in Turkish breeds. Our findings provide valuable insights into the genetic landscape of Turkish cattle, offering a crucial foundation for informed conservation and breeding strategies aimed at preserving these breeds' genetic integrity and heritage.

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.

Production systems and breeding practices of Begaria cattle breed as input for a community-based breeding program

This study aimed to describe the production systems and breeding practices of Begaria cattle breeds as a resource for community-based conservation and breeding program. A semi-structured questionnaire and focus group discussions were employed to identify and describe the breeding practice of the farmer. Three kebeles: Almahal, Fanguso, and Mankush were considered for the study. The data were analyzed using R software. Cattle were kept for a variety of reasons, including milk, breeding, meat, cash, and saving. It was typical throughout the study area to see herd mixing and keeping together (4.7 ± 0.33 herds on average), herd movement during the dry season, and breeding bull sharing (1.9 ± 1.55 shared bulls). During the summer season, cattle rearing in the areas was hindered by feed and water shortages. The breed's major reproduction constraints were abortion (24 %), repeat breeding (19), and calf mortality (14). The primary reason for keeping Begaria cattle was for milk (29 %), cash (26), meat (22), and breeding (0.22). The average age at first calving (AFC), calving interval (CI) and milk yield per day of the breed were 41.5 ± 1.63, 16.0 ± 0.43 months, and 3.2 ± 0.12 L respectively. The average age of breeding bull selection was 3 ± 0.1 years, and it was based on several criteria, including their phenotype, growth, and white coat colour. Based on the current findings, it is recommended that a multi-trait selection program be designed with full stakeholder participation.

Genetic Distinctness and Diversity of American Aberdeen Cattle Compared to Common Beef Breeds in the United States

American Aberdeen (AD) cattle in the USA descend from an Aberdeen Angus herd originally brought to the Trangie Agricultural Research Centre, New South Wales, AUS. Although put under specific selection pressure for yearling growth rate, AD remain genomically uncharacterized. The objective was to characterize the genetic diversity and structure of purebred and crossbred AD cattle relative to seven common USA beef breeds using available whole-genome SNP data. A total of 1140 animals consisting of 404 purebred (n = 8 types) and 736 admixed individuals (n = 10 types) was used. Genetic diversity metrics, an analysis of molecular variance, and a discriminant analysis of principal components were employed. When linkage disequilibrium was not accounted for, markers influenced basic diversity parameter estimates, especially for AD cattle. Even so, intrapopulation and interpopulation estimates separate AD cattle from other purebred types (e.g., Latter's pairwise FST ranged from 0.1129 to 0.2209), where AD cattle were less heterozygous and had lower allelic richness than other purebred types. The admixed AD-influenced cattle were intermediate to other admixed types for similar parameters. The diversity metrics separation and differences support strong artificial selection pressures during and after AD breed development, shaping the evolution of the breed and making them genomically distinct from similar breeds.

Understanding the underlying genetic mechanisms for age at first calving, inter-calving period and scrotal circumference in Bonsmara cattle

BACKGROUND

Reproduction is a key feature of the sustainability of a species and thus represents an important component in livestock genetic improvement programs. Most reproductive traits are lowly heritable. In order to gain a better understanding of the underlying genetic basis of these traits, a genome-wide association was conducted for age at first calving (AFC), first inter-calving period (ICP) and scrotal circumference (SC) within the South African Bonsmara breed. Phenotypes and genotypes (120,692 single nucleotide polymorphisms (SNPs) post editing) were available on 7,128 South African Bonsmara cattle; the association analyses were undertaken using linear mixed models.

RESULTS

Genomic restricted maximum likelihood analysis of the 7,128 SA Bonsmara cattle yielded genomic heritability's of 0.183 (SE = 0.021) for AFC, 0.207 (SE = 0.022) for ICP and 0.209 (SE = 0.019) for SC. A total of 16, 23 and 51 suggestive (P ≤ 4 × 10-6) SNPs were associated with AFC, ICP and SC, while 11, 11 and 44 significant (P ≤ 4 × 10-7) SNPs were associated with AFC, ICP and SC respectively. A total of 11 quantitative trait loci (QTL) and 11 candidate genes were co-located with these associated SNPs for AFC, with 10 QTL harbouring 11 candidate genes for ICP and 41 QTL containing 40 candidate genes for SC. The QTL identified were close to genes previously associated with carcass, fertility, growth and milk-related traits. The biological pathways influenced by these genes include carbohydrate catabolic processes, cellular development, iron homeostasis, lipid metabolism and storage, immune response, ovarian follicle development and the regulation of DNA transcription and RNA translation.

CONCLUSIONS

This was the first attempt to study the underlying polymorphisms associated with reproduction in South African beef cattle. Genes previously reported in cattle breeds for numerous traits bar AFC, ICP or SC were detected in this study. Over 20 different genes have not been previously reported in beef cattle populations and may have been associated due to the unique genetic composite background of the SA Bonsmara breed.

Population Genomic Evidence for the Diversification of Bellamya aeruginosa in Different River Systems in China

Clarifying the genetic structure can facilitate the understanding of a species evolution history. It is crucial for the management of germplasm resources and providing useful guidance for effective selective breeding. Bellamya is an economically and ecologically important freshwater snail for fish, birds and even humans. Population genetic structures of the Bellamya species, however, were unknown in previous studies. Population genomics approaches with tens to hundreds of thousands of single nucleotide polymorphisms (SNPs) make it possible to detect previously unidentified structures. The population genomic study of seven populations of B. aeruginosa across three river systems (Yellow River, Yangtze River and Pearl River) in China was conducted by SLAF-seq. SLAF-seq obtained a total of 4737 polymorphisms SLAF-tags and 25,999 high-consistency genome-wide SNPs. The population genetic structure showed a clear division among populations from the Yellow River basin (YH and WL) and the Pearl River basin (QSH and LB), as well as population YC from the Yangtze River basin using the SNPs data. However, there existed no distinct population structure using the mitochondrial DNA (mtDNA). The anthropogenic translocation from the Yangtze River basin to the Pearl River basin and the passive dispersion from the Yangtze River basin to the Yellow River basin by flooding have weakened the phylogeographic pattern of B. aeruginosa. The divergence of B. aeruginosa in the three river systems suggests that the anthropogenic dispersal for aquaculture and breeding requires serious consideration of the population structure for the preservation of genetic diversity and effective utilization of germplasm resources.

Revealing Genetic Diversity and Population Structure of Endangered Altay White-Headed Cattle Population Using 100 k SNP Markers

Understanding the genetic basis of native cattle populations that have adapted to the local environment is of great significance for formulating appropriate strategies and programs for genetic improvement and protection. Therefore, it is necessary to understand the genetic diversity and population structure of Altay white-headed cattle so as to meet the current production needs under various environments, carry out continuous genetic improvement, and promote rapid adaptation to changing environments and breeding objectives. A total of 46 individual samples of endangered Xinjiang Altay white-headed cattle were collected in this study, including nine bulls and 37 cows. To collect genotype data, 100 k SNP markers were used, and then studies of genetic diversity, genetic structure, inbreeding degree, and family analysis were carried out. A total of 101,220 SNP loci were detected, and the genotype detection rate for individuals was ≥90%. There were 85,993 SNP loci that passed quality control, of which 93.5% were polymorphic. The average effective allele number was 0.036, the Polymorphism Information Content was 0.304 and the minimum allele frequency was 0.309, the average observed heterozygosity was 0.413, and the average expected heterozygosity was 0.403. The average genetic distance of Idengtical By State (IBS) was 0.3090, there were 461 ROH (genome-length homozygous fragments), 76.1% of which were between 1 and 5 MB in length, and the average inbreeding coefficient was 0.016. The 46 Altay white-headed cattle were divided into their families, and the individual numbers of each family were obviously different. To sum up, the Altay white-headed cattle conservation population had low heterozygosity, a high inbreeding degree, few families, and large differences in the number of individuals in each family, which can easily cause a loss of genetic diversity. In the follow-up seed conservation process, seed selection and matching should be carried out according to the divided families to ensure the long-term protection of Altay white-headed cattle genetic resources.

Selection Signatures in South African Nguni and Bonsmara Cattle Populations Reveal Genes Relating to Environmental Adaptation

Climate change is a major influencing factor in beef production. The greenhouse gases produced from livestock production systems contribute to the overall greenhouse gas emissions. The aim of this study was to identify selection signatures within and between Nguni and Bonsmara cattle in relation to production and adaptation. For this purpose, genomic 150 K single nucleotide polymorphism data from Nguni (n = 231) and Bonsmara (n = 252) cattle in South Africa were used. Extended haplotype homozygosity (EHH) based analysis was executed within each population using integrated haplotype score (iHS). The R package rehh was used for detecting selection signatures across the two populations with cross population EHH (XP-EHH). Total of 121 regions of selection signatures were detected (p < 0.0001) in the Bonsmara and Nguni populations. Several genes relating to DNA methylation, heat stress, feed efficiency and nitrogen metabolism were detected within and between each population. These regions also included QTLs associated with residual feed intake, residual gain, carcass weight, stature and body weight in the Bonsmara, while QTLs associated with conception rate, shear force, tenderness score, juiciness, temperament, heat tolerance, feed efficiency and age at puberty were identified in Nguni. Based on the results of the study it is recommended that the Nguni and Bonsmara be utilized in crossbreeding programs as they have beneficial traits that may allow them to perform better in the presence of climate change. Results of this study coincide with Nguni and Bonsmara breed characteristics and performance, and furthermore support informative crossbreeding programs to enhance livestock productivity in South Africa.

Genetics of Base Coat Colour Variations and Coat Colour-Patterns of the South African Nguni Cattle Investigated Using High-Density SNP Genotypes

Nguni cattle are a Sanga type breed with mixed B. taurus and B. indicus ancestry and proven resistance to ticks, diseases and other harsh conditions of the African geographical landscape. The multi-coloured Nguni coats have found a niche market in the leather industry leading to breeding objectives towards the promotion of such diversity. However, there is limited studies on the genomic architecture underlying the coat colour and patterns hampering any potential breeding and improvement of such trait. This study investigated the genetics of base coat colour, colour-sidedness and the white forehead stripe in Nguni cattle using coat colour phenotyped Nguni cattle and Illumina Bovine HD (770K) genotypes. Base coat colour phenotypes were categorised into eumelanin (n = 45) and pheomelanin (n = 19). Animals were categorised into either colour-sided (n = 46) or non-colour-sided (n = 94) and similarly into presence (n = 15) or absence (n = 67) of white forehead stripe. Genome-wide association tests were conducted using 622,103 quality controlled SNPs and the Efficient Mixed Model Association eXpedited method (EMMAX) implemented in Golden Helix SNP Variation Suite. The genome-wide association studies for base coat colour (eumelanin vs. pheomelanin) resulted into four indicative SNPs on BTA18 and a well-known gene, MC1R, was observed within 1 MB from the indicative SNPs (p < 0.00001) and found to play a role in the melanogenesis (core pathway for melanin production) and the MAPK signalling pathway. GWAS for colour-sidedness resulted in four indicative SNPs, none of which were in close proximity to the KIT candidate gene known for colour-sidedness. GWAS for the white forehead stripe resulted in 17 indicative SNPs on BTA6. Four genes MAPK10, EFNA5, PPP2R3C and PAK1 were found to be associated with the white forehead stripe and were part of the MAPK, adrenergic and Wnt signalling pathways that are synergistically associated with the synthesis of melanin. Overall, our results prove prior knowledge of the role of MC1R in base coat colours in cattle and suggested a different genetic mechanism for forehead stripe phenotypes in Nguni cattle.

Genetic Differentiation among Livestock Breeds—Values for Fst

Simple Summary

The degree of relationship among livestock breeds can be quantified by the F_st statistic, which measures the extent of genetic differentiation between them. An F_st value of 0.1 has often been taken as indicating that two breeds are indeed genetically distinct, but this concept has not been evaluated critically. Here, F_st values have been collated for the six major livestock species: cattle, sheep, goats, pigs, horses, and chickens. These values are remarkably variable both within and between species, demonstrating that F_st > 0.1 is not a reliable criterion for breed distinctiveness. However, the large body of F_st data accumulated in the last 20–30 years represents an untapped database that could contribute to the development of interdisciplinary research involving livestock breeds.

Abstract

(1) Background: The F_st statistic is widely used to characterize between-breed relationships. F_st = 0.1 has frequently been taken as indicating genetic distinctiveness between breeds. This study investigates whether this is justified. (2) Methods: A database was created of 35,080 breed pairs and their corresponding F_st values, deduced from microsatellite and SNP studies covering cattle, sheep, goats, pigs, horses, and chickens. Overall, 6560 (19%) of breed pairs were between breeds located in the same country, 7395 (21%) between breeds of different countries within the same region, 20,563 (59%) between breeds located far apart, and 562 (1%) between a breed and the supposed wild ancestor of the species. (3) Results: General values for between-breed F_st were as follows, cattle: microsatellite 0.06–0.12, SNP 0.08–0.15; sheep: microsatellite 0.06–0.10, SNP 0.06–0.17; horses: microsatellite 0.04–0.11, SNP 0.08–0.12; goats: microsatellite 0.04–0.14, SNP 0.08–0.16; pigs: microsatellite 0.06–0.27, SNP 0.15–0.22; chickens: microsatellite 0.05–0.28, SNP 0.08–0.26. (4) Conclusions: (1) Large amounts of F_st data are available for a substantial proportion of the world’s livestock breeds, (2) the value for between-breed F_st of 0.1 is not appropriate owing to its considerable variability, and (3) accumulated F_st data may have value for interdisciplinary research.

Sustained Effects of Muscle Calpain System Genotypes on Tenderness Phenotypes of South African Beef Bulls during Ageing up to 20 Days

Simple Summary

When searching for genetic markers for the selection of more tender beef, it is important to maintain minimal environmental variation from pre-slaughter, right through to the ageing process, to ensure the accuracy of the obtained phenotypes. This is because beef quality traits have a large environmental component that can greatly alter the characteristics of the meat, which would not reflect a true genetic effect. We propose that variable ageing times are especially important in determining whether markers are associated with tenderization or not. Our analyses included candidate genes for the protein degrading enzyme system for calpains, because they contribute the most to tenderization. We were able to validate these markers in South African beef cattle, where they could be useful for selection. The timing of the collection of tenderness data was critical, as only a few (6/134) genetic markers sustained their association with tenderization over ageing to 20 days. A larger tenderization effect earlier in ageing, as shown here for the capn1_187 and capn1_4751 markers, would decrease the length of ageing. This would not only increase profits, but also decrease the energy needed during the storage and refrigeration of aged beef, decreasing the carbon footprint of beef production.

Abstract

The most important factor that determines beef tenderness is its proteolytic activity, and the balance between calpain-1 protease activity and calpastatin inhibition is especially important, while contributions can also arise from calpain-2 and, possibly, calpain-3. The meat ageing process itself affects these processes. To determine whether genotypes in the calpain–calpastatin system can enhance tenderness through a 20-day ageing period, South African purebred beef bulls (n = 166) were genotyped using the Illumina BovineHD SNP BeadChip through a gene-based association analysis targeting the cast, capn3, capn2 and capn1 genes. The Warner–Bratzler shear force (WBSF) and myofibril fragment length (MFL) of Longissimus thoracis et lumborum (LTL) steaks were evaluated between d 3 and d 20 of ageing, with protease enzyme activity in the first 20 h post-mortem. Although several of the 134 SNPs are associated with tenderness, only seven SNP in the cast, capn2 and capn1 genes sustained genetic associations, additive to the ageing-associated increases in tenderness for at least three of the four ageing periods. While most genomic associations were relatively stable over time, some genotypes within the SNP responded differently to ageing, resulting in altered genomic effects over time. The level of ageing at which genomic associations are performed is an important factor that determines whether SNPs affect tenderness phenotypes.

Genome analyses revealed genetic admixture and selection signatures in Bos indicus

The genomic diversity and relationship among seven diverse cattle breeds viz. Sahiwal, Tharparkar, Gir, Vechur, Ongole, Kangayam and Hariana were investigated in 132 random samples based on high density SNP array comprising > 777 K SNPs. A total of 1993 SNPs (0.25% of the total) having greater power (F_ST ≥ 0.20) to differentiate these cattle populations were identified, and utilized to partition genome of each animal into a predefined number of clusters. The structure of these cattle indicated shared ancestry of dairy breeds viz. Gir, Tharparkar and Sahiwal. Most of the animals (> 76%) of different populations under study except Vechur clustered into their own group of animals called breed. Vechur population retained highest rate of admixture, consistent with its crossing with other breeds. Ongole, Kangayam and Hariana shared comparatively less of their genome (≤ 15%) with other breeds. The study indicated that all seven breeds evolved from their independent ancestry but there was intermixing of these breeds in the recent past. The selection signatures identified between draft (Kangayam) and dairy breeds included several genes like FAM19A2, RAB31P, BEST3, DGKA, AHCY, PIGU and PFKP which are involved in immune response, metabolic pathway, transportation of glucose and sugars, signaling pathways, cellular processes, cell division and glycolysis regulation, respectively. Moreover, these genomic regions also harbour QTLs affecting milk performance traits. The signatures were also identified even between the dairy breeds. In comparison to large-sized cattle, there were significant differences in the number of QTLs affecting production (body weight, growth rate etc.) and morphological traits (height) in short-statured Vechur breed. The presence of HMGA2 gene in the selection signature on chromosome 5 may explain the variations in stature between these cattle.

An Overview of the Use of Genotyping Techniques for Assessing Genetic Diversity in Local Farm Animal Breeds

Simple Summary

The number of local farm animal breeds is declining worldwide. However, these breeds have different degrees of genetic diversity. Measuring genetic diversity is important for the development of conservation strategies and, therefore, various genomic analysis techniques are available. The aim of the present work was to shed light on the use of these techniques in diversity studies of local breeds. In summary, a total of 133 worldwide studies that examined genetic diversity in local cattle, sheep, goat, chicken and pig breeds were reviewed. The results show that over time, almost all available genomic techniques were used and various diversity parameters were calculated. Therefore, the present results provide a comprehensive overview of the application of these techniques in the field of local breeds. This can provide helpful insights into the advancement of the conservation of breeds with high genetic diversity.

Abstract

Globally, many local farm animal breeds are threatened with extinction. However, these breeds contribute to the high amount of genetic diversity required to combat unforeseen future challenges of livestock production systems. To assess genetic diversity, various genotyping techniques have been developed. Based on the respective genomic information, different parameters, e.g., heterozygosity, allele frequencies and inbreeding coefficient, can be measured in order to reveal genetic diversity between and within breeds. The aim of the present work was to shed light on the use of genotyping techniques in the field of local farm animal breeds. Therefore, a total of 133 studies across the world that examined genetic diversity in local cattle, sheep, goat, chicken and pig breeds were reviewed. The results show that diversity of cattle was most often investigated with microsatellite use as the main technique. Furthermore, a large variety of diversity parameters that were calculated with different programs were identified. For 15% of the included studies, the used genotypes are publicly available, and, in 6%, phenotypes were recorded. In conclusion, the present results provide a comprehensive overview of the application of genotyping techniques in the field of local breeds. This can provide helpful insights to advance the conservation of breeds.

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.

Breed Ancestry, Divergence, Admixture, and Selection Patterns of the Simbra Crossbreed

In this study, we evaluated an admixed South African Simbra crossbred population, as well as the Brahman (Indicine) and Simmental (Taurine) ancestor populations to understand their genetic architecture and detect genomic regions showing signatures of selection. Animals were genotyped using the Illumina BovineLD v2 BeadChip (7K). Genomic structure analysis confirmed that the South African Simbra cattle have an admixed genome, composed of 5/8 Taurine and 3/8 Indicine, ensuring that the Simbra genome maintains favorable traits from both breeds. Genomic regions that have been targeted by selection were detected using the linkage disequilibrium-based methods iHS and Rsb. These analyses identified 10 candidate regions that are potentially under strong positive selection, containing genes implicated in cattle health and production (e.g., TRIM63, KCNA10, NCAM1, SMIM5, MIER3, and SLC24A4). These adaptive alleles likely contribute to the biological and cellular functions determining phenotype in the Simbra hybrid cattle breed. Our data suggested that these alleles were introgressed from the breed's original indicine and taurine ancestors. The Simbra breed thus possesses derived parental alleles that combine the superior traits of the founder Brahman and Simmental breeds. These regions and genes might represent good targets for ad-hoc physiological studies, selection of breeding material and eventually even gene editing, for improved traits in modern cattle breeds. This study represents an important step toward developing and improving strategies for selection and population breeding to ultimately contribute meaningfully to the beef production industry.

Genetic diversity and population structure of three native cattle populations in Mozambique

In the present work, the population diversity and structure of three populations of native Mozambican cattle were studied, to develop knowledge that is required for sound conservation and genetic improvement programs of these genetic resources. A total of 228 animals (Landim, Angone, and Tete) were genotyped using the International Dairy and Beef version three (IDBV3) SNP BeadChip array. Population parameters varied within a limited scope, with the average minor allele frequency (MAF) ranging from 0.228 ± 0.154 in the Angone to 0.245 ± 0.145 in the Tete population, while estimates of expected heterozygosities varied from 0.304 ± 0.166 in the Angone to 0.329 ± 0.148 in the Tete population. Low positive (0.065 ± 0.109) inbreeding rates were detected in the three cattle groups. Population structure and admixture analyses indicated low genetic differentiation and various degrees of admixture among the populations. The effective population size has decreased over time and at 12 generations ago ranged between 349 (Tete) and 929 (Landim). The average linkage disequilibrium (LD) of the studied populations ranged from 0.400 ± 0.213 (Tete) to 0.434 ± 0.232 (Landim). The findings of this study will be valuable for formulating management and conservation strategies for indigenous Mozambican cattle populations.

Genetic characterisation of non-descript cattle populations in communal areas of South Africa

Context Indigenous cattle breeds represent an important genetic resource for livelihood of communal-area inhabitants. Indigenous breeds have the ability to withstand harsh climatic conditions, can adapt genetically to poor-quality forages and are resistant to parasites and diseases. These unique traits possessed by indigenous breeds are under threat because of unrestrained crossing with exotic commercial breeds, and this can lead to total loss of a breed. Aims The study was conducted to assess the genetic diversity and population structure of South African non-descript communal beef cattle populations by using 25 microsatellite markers. Methods Unrelated and non-descript animals (n = 150) were sampled from communal areas from five (5) provinces of South Africa, namely, Eastern Cape, KwaZulu–Natal, Limpopo, Mpumalanga and the North West, with 30 samples per breed taken. Six (6) known cattle breeds (n = 180) were used as a reference population. This included Angus, Afrikaner, Bonsmara, Brahman, Drakensberger and the Nguni, with 30 samples per breed. Key results High level of genetic diversity was found across the five non-descript populations, with an average heterozygosity of 75%. The Limpopo population was found to be the most diverse population, with the highest average number of alleles (8.5) and heterozygosity (ranging between observed heterozygosity of 70% and expected heterozygosity of 79%). STRUCTURE software assigned populations (2 ≤ K ≤ 20), with the most probable cluster being at K = 7. The Eastern Cape, KwaZulu–Natal and Limpopo populations had genetic material similar to those possessed by the Nguni and Bonsmara reference populations. Conclusions Results from the study showed that most genetic differentiation occurred within populations rather than among populations, and this might be due to the fact that there is no selection for or against any specific production trait expressed in the populations. Implications The obtained information will serve as a baseline for the development and implementation of sound breeding programs that will assist in controlling the gene flow, so as to lower the possible genetic dilution of the currently available genetic material.

Population demographic history and population structure for Pakistani Nili-Ravi breeding bulls based on SNP genotyping to identify genomic regions associated with male effects for milk yield and body weight

The domestic Nili-Ravi water buffalo (Bubalus bubalis) is the best dairy animal contributing 68% to total milk production in Pakistan. In this study, we identified genome-wide single nucleotide polymorphisms (SNPs) to estimate various population genetic parameters such as diversity, pairwise population differentiation, linkage disequilibrium (LD) distribution and for genome-wide association study for milk yield and body weight traits in the Nili-Ravi dairy bulls that they may pass on to their daughters who are retained for milking purposes. The genotyping by sequencing approach revealed 13,039 reference genome-anchored SNPs with minor allele frequency of 0.05 among 167 buffalos. Population structure analysis revealed that the bulls were grouped into two clusters (K = 2), which indicates the presence of two different lineages in the Pakistani Nili-Ravi water buffalo population, and we showed the extent of admixture of these two lineages in our bull collection. LD analysis revealed 4169 significant SNP associations, with an average LD decay of 90 kb for these buffalo genome. Genome-wide association study involved a multi-locus mixed linear model for milk yield and body weight to identify genome-wide male effects. Our study further illustrates the utility of the genotyping by sequencing approach for identifying genomic regions to uncover additional demographic complexity and to improve the complex dairy traits of the Pakistani Nili-Ravi water buffalo population that would provide the lot of economic benefits to dairy industry.

Evolution and domestication of the Bovini species

Domestication of the Bovini species (taurine cattle, zebu, yak, river buffalo and swamp buffalo) since the early Holocene (ca. 10 000 BCE) has contributed significantly to the development of human civilization. In this study, we review recent literature on the origin and phylogeny, domestication and dispersal of the three major Bos species - taurine cattle, zebu and yak - and their genetic interactions. The global dispersion of taurine and zebu cattle was accompanied by population bottlenecks, which resulted in a marked phylogeographic differentiation of the mitochondrial and Y-chromosomal DNA. The high diversity of European breeds has been shaped through isolation-by-distance, different production objectives, breed formation and the expansion of popular breeds. The overlapping and broad ranges of taurine and zebu cattle led to hybridization with each other and with other bovine species. For instance, Chinese gayal carries zebu mitochondrial DNA; several Indonesian zebu descend from zebu bull × banteng cow crossings; Tibetan cattle and yak have exchanged gene variants; and about 5% of the American bison contain taurine mtDNA. Analysis at the genomic level indicates that introgression may have played a role in environmental adaptation.

Genetic characteristics of Korean Jeju Black cattle with high density single nucleotide polymorphisms

Objective

Conservation and genetic improvement of cattle breeds require information about genetic diversity and population structure of the cattle. In this study, we investigated the genetic diversity and population structure of the three cattle breeds in the Korean peninsula.

Methods

Jeju Black, Hanwoo, Holstein cattle in Korea, together with six foreign breeds were examined. Genetic diversity within the cattle breeds was analyzed with minor allele frequency (MAF), observed and expected heterozygosity (H_O and H_E), inbreeding coefficient (F_IS) and past effective population size. Molecular variance and population structure between the nine breeds were analyzed using a model-based clustering method. Genetic distances between breeds were evaluated with Nei’s genetic distance and Weir and Cockerham’s F_ST.

Results

Our results revealed that Jeju Black cattle had lowest level of heterozygosity (H_E = 0.21) among the studied taurine breeds, and an average MAF of 0.16. The level of inbreeding was −0.076 for Jeju Black, while −0.018 to −0.118 for the other breeds. Principle component analysis and neighbor-joining tree showed a clear separation of Jeju Black cattle from other local (Hanwoo and Japanese cattle) and taurine/indicine cattle breeds in evolutionary process, and a distinct pattern of admixture of Jeju Black cattle having no clustering with other studied populations. The F_ST value between Jeju Black cattle and Hanwoo was 0.106, which was lowest across the pair of breeds ranging from 0.161 to 0.274, indicating some degree of genetic closeness of Jeju Black cattle with Hanwoo. The past effective population size of Jeju Black cattle was very small, i.e. 38 in 13 generation ago, whereas 209 for Hanwoo.

Conclusion

This study indicates genetic uniqueness of Jeju Black cattle. However, a small effective population size of Jeju Black cattle indicates the requirement for an implementation of a sustainable breeding policy to increase the population for genetic improvement and future conservation.

Genome-wide assessment of genetic diversity, linkage disequilibrium and haplotype block structure in Tharparkar cattle breed of India

Knowledge about genetic diversity is very essential for the management and sustainable utilization of livestock genetic resources. In this study, we presented a comprehensive genome-wide analysis of genetic diversity, ROH, inbreeding, linkage disequilibrium, effective population size and haplotype block structure in Tharparkar cattle of India. A total of 24 Tharparkar animals used in this study were genotyped with Illumina BovineSNP50 array. After quality control, 22,825 biallelic SNPs were retained, which were in HWE, MAF > 0.05 and genotyping rate >90%. The overall mean observed (H_O) and expected heterozygosity (H_E) were 0.339 ± 0.156 and 0.325 ± 0.129, respectively. The average minor allele frequency was 0.234 with a standard deviation of ± 0.131. We identified a total of 1832 ROH segments and the highest autosomal coverage of 13.87% was observed on chromosome 23. The genomic inbreeding coefficients estimates by F_ROH, F_HOM, F_GRM and F_UNI were 0.0589, 0.0215, 0.0532 and 0.0160 respectively. The overall mean linkage disequilibrium (LD) for a total of 133,532 pairwise SNPs measured by D' and r² was 0.6452 and 0.1339, respectively. In addition, we observed a gradual decline in effective population size over the past generations.

Genome-wide diversity and global migration patterns in dromedaries follow ancient caravan routes

Dromedaries have been essential for the prosperity of civilizations in arid environments and the dispersal of humans, goods and cultures along ancient, cross-continental trading routes. With increasing desertification their importance as livestock species is rising rapidly, but little is known about their genome-wide diversity and demographic history. As previous studies using few nuclear markers found weak phylogeographic structure, here we detected fine-scale population differentiation in dromedaries across Asia and Africa by adopting a genome-wide approach. Global patterns of effective migration rates revealed pathways of dispersal after domestication, following historic caravan routes like the Silk and Incense Roads. Our results show that a Pleistocene bottleneck and Medieval expansions during the rise of the Ottoman empire have shaped genome-wide diversity in modern dromedaries. By understanding subtle population structure we recognize the value of small, locally adapted populations and appeal for securing genomic diversity for a sustainable utilization of this key desert species.

The genetic heritage of Alpine local cattle breeds using genomic SNP data

BACKGROUND

Assessment of genetic diversity and population structure provides important control metrics to avoid genetic erosion, inbreeding depression and crossbreeding between exotic and locally-adapted cattle breeds since these events can have deleterious consequences and eventually lead to extinction. Historically, the Alpine Arc represents an important pocket of cattle biodiversity with a large number of autochthonous breeds that provide a fundamental source of income for the entire regional economy. By using genotype data from medium-density single nucleotide polymorphism (SNP) arrays, we performed a genome-wide comparative study of 23 cattle populations from the Alpine Arc and three cosmopolitan breeds.

RESULTS

After filtering, we obtained a final genotyping dataset consisting of 30,176 SNPs for 711 individuals. The local breeds showed high or intermediate values of genetic diversity compared to the highly selected cosmopolitan breeds. Patterns of genetic differentiation, multidimensional scaling, admixture analysis and the constructed phylogenetic tree showed convergence, which indicates the presence of gene flow among the breeds according to both geographic origin and historical background. Among the most differentiated breeds, we identified the modern Brown cattle. In spite of admixture events, several local breeds have preserved distinctive characteristics, which is probably due to differences in genetic origin and geographic location.

CONCLUSIONS

This study represents one of the most comprehensive genome-wide analysis of the Alpine cattle breeds to date. Using such a large dataset that includes the majority of the local breeds found in this region, allowed us to expand knowledge on the evaluation and status of Alpine cattle biodiversity. Our results indicate that although many of the analyzed local breeds are listed as endangered, they still harbor a large amount of genetic diversity, even when compared to some cosmopolitan breeds. This finding, together with the reconstruction of the phylogeny and the relationships between these Alpine Arc cattle breeds, provide crucial insights not only into the improvement of genetic stocks but also into the implementation of future conservation strategies.

Genome Wide Assessment of Genetic Variation and Population Distinctiveness of the Pig Family in South Africa

Genetic diversity is of great importance and a prerequisite for genetic improvement and conservation programs in pigs and other livestock populations. The present study provides a genome wide analysis of the genetic variability and population structure of pig populations from different production systems in South Africa relative to global populations. A total of 234 pigs sampled in South Africa and consisting of village (n = 91), commercial (n = 60), indigenous (n = 40), Asian (n = 5) and wild (n = 38) populations were genotyped using Porcine SNP60K BeadChip. In addition, 389 genotypes representing village and commercial pigs from America, Europe, and Asia were accessed from a previous study and used to compare population clustering and relationships of South African pigs with global populations. Moderate heterozygosity levels, ranging from 0.204 for Warthogs to 0.371 for village pigs sampled from Capricorn municipality in Eastern Cape province of South Africa were observed. Principal Component Analysis of the South African pigs resulted in four distinct clusters of (i) Duroc; (ii) Vietnamese; (iii) Bush pig and Warthog and (iv) a cluster with the rest of the commercial (SA Large White and Landrace), village, Wild Boar and indigenous breeds of Koelbroek and Windsnyer. The clustering demonstrated alignment with genetic similarities, geographic location and production systems. The PCA with the global populations also resulted in four clusters that where populated with (i) all the village populations, wild boars, SA indigenous and the large white and landraces; (ii) Durocs (iii) Chinese and Vietnamese pigs and (iv) Warthog and Bush pig. K = 10 (The number of population units) was the most probable ADMIXTURE based clustering, which grouped animals according to their populations with the exception of the village pigs that showed presence of admixture. AMOVA reported 19.92%-98.62% of the genetic variation to be within populations. Sub structuring was observed between South African commercial populations as well as between Indigenous and commercial breeds. Population pairwise F ST analysis showed genetic differentiation (P ≤ 0.05) between the village, commercial and wild populations. A per marker per population pairwise F ST analysis revealed SNPs associated with QTLs for traits such as meat quality, cytoskeletal and muscle development, glucose metabolism processes and growth factors between both domestic populations as well as between wild and domestic breeds. Overall, the study provided a baseline understanding of porcine diversity and an important foundation for porcine genomics of South African populations.

Population structure of indigenous southern African goats based on the Illumina Goat50K SNP panel

In this study, the genetic structure of indigenous Tswana and Swazi goats using the Illumina Goat50K SNP array was investigated. Two South African commercial goat breeds were included to investigate admixture with the indigenous populations in southern Africa. A total of 144 DNA samples including Boer goats (n = 24), Kalahari Red (n = 24), Swazi (n = 48), and Tswana goats (n = 48) were genotyped. Statistical analysis was performed using PLINK version 1.07. Genetic diversity, measured as expected heterozygosity, was estimated at 0.390, 0.398, 0.413, and 0.387 for Boer, Kalahari Red, Tswana, and Swazi goats, respectively. The individual inbreeding coefficient varied from 0.019 ± 0.05 to 0.011 ± 0.06 for the Tswana and Swazi goats, respectively. The Principal component analysis clustered the populations according to geographical origin and breed type. Linkage disequilibrium (LD) for shorter intervals (0-10 kb) ranged from 0.44 to 0.56 and commercial breeds had higher values. Effective population sizes decreased with generations and at the 13th generation ranged between 87 for Boer to 266 for Tswana goats. The Tswana population exhibited the highest level of genetic variation and effective population size, which holds potential for improved production in marginal regions. A national strategy is required to maintain genetic diversity in communal goat production systems through well-structured breeding and conservation programs.

Strategies for Sustainable Use of Indigenous Cattle Genetic Resources in Southern Africa

Indigenous cattle breeds are the most important livestock species in the Southern African Development Community (SADC) region owing to their role in human food, nutrition, income, and social security. Despite the role of these breeds in the household and national economies, they are currently underutilised, their productivity remains low, and populations are faced with extinction. In addition, there are insufficient measures taken to secure their present and future value. The current review highlights strategies for sustainable use of indigenous cattle genetic resources in the region, including the use of novel production and marketing practices, women and youth empowerment, and development of the appropriate capacity building, legislative, and policy structures. At present, the lack of coordination among the different stakeholders still poses a challenge to the implementation of these strategies. To this end, partnerships, collaboration, and stakeholders’ participation are recommended to effectively implement strategies for sustainable use of indigenous cattle breeds.

Genetic variability and relationships in nine South African cattle breeds using microsatellite markers

Genetic variability within and between breeds allows adaptation to a changing environment and consequently prepares producers for the future. Eleven bovine-specific microsatellite markers were used to genotype animals from each of nine South African cattle breeds: Afrikaner (N = 550), Angus (N = 550), Bonsmara (N = 550), Boran (N = 321), Brahman (N = 550), Drakensberger (N = 550), Nguni (N = 550), Simmental (N = 550), and Tuli (N = 311). These breeds were drawn from Bos taurus africanus, Bos taurus, and Bos indicus. Genetic variability estimates included unbiased heterozygosity, effective number of alleles, and inbreeding. Ranges of these parameters were 0.569-0.741, 8.818-11.455, and - 0.001-0.050, respectively. Breed private allele and breed pairwise comparison was also used to characterize the breeds. The analysis of population structure with K = 2 revealed clusters comprised of Sanga-indicine and taurine, while K = 3 included separate clusters of Sanga, indicine, and taurine, and with K = 9 showed the breeds arising from unique progenitor populations. This study broke new ground in molecular cattle genetic diversity by genotyping a large sample size per breed and using a larger number of breeds compared with similar studies that have been conducted in the recent past which have either used a smaller number of breeds or smaller sample sizes but with a larger number of marker loci. Thus, opportunities that arise to explore genetic diversity and relationships in both the livestock and wildlife industries in Southern Africa may capitalize on microsatellite marker databases which remain cost-effective and accessible due to their extensive use for parentage verification.

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

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

Genome-Wide SNP Discovery in Indigenous Cattle Breeds of South Africa

Single nucleotide polymorphism arrays have created new possibilities for performing genome-wide studies to detect genomic regions harboring sequence variants that affect complex traits. However, the majority of validated SNPs for which allele frequencies have been estimated are limited primarily to European breeds. The objective of this study was to perform SNP discovery in three South African indigenous breeds (Afrikaner, Drakensberger, and Nguni) using whole genome sequencing. DNA was extracted from blood and hair samples, quantified and prepared at 50 ng/μl concentration for sequencing at the Agricultural Research Council Biotechnology Platform using an Illumina HiSeq 2500. The fastq files were used to call the variants using the Genome Analysis Tool Kit. A total of 1,678,360 were identified as novel using Run 6 of 1000 Bull Genomes Project. Annotation of the identified variants classified them into functional categories. Within the coding regions, about 30% of the SNPs were non-synonymous substitutions that encode for alternate amino acids. The study of distribution of SNP across the genome identified regions showing notable differences in the densities of SNPs among the breeds and highlighted many regions of functional significance. Gene ontology terms identified genes such as MLANA, SYT10, and CDC42EP5 that have been associated with coat color in mouse, and ADAMS3, DNAJC3, and PAG5 genes have been associated with fertility in cattle. Further analysis of the variants detected 688 candidate selective sweeps (ZHp Z-scores ≤ -4) across all three breeds, of which 223 regions were assigned as being putative selective sweeps (ZHp scores ≤-5). We also identified 96 regions with extremely low ZHp Z-scores (≤-6) in Afrikaner and Nguni. Genes such as KIT and MITF that have been associated with skin pigmentation in cattle and CACNA1C, which has been associated with biopolar disorder in human, were identified in these regions. This study provides the first analysis of sequence data to discover SNPs in indigenous South African cattle breeds. The information will play an important role in our efforts to understand the genetic history of our cattle and in designing appropriate breed improvement programmes.

Evaluation of two bovine SNP genotyping arrays for breed clustering and stratification analysis in well-known taurine and indicine breeds

The present study aimed to evaluate the efficiency of two Bovine SNP genotyping arrays (i.e., 50 K and HD) for breed clustering and stratification related studies in taurine and indicine breeds. The whole-genome SNP data at two densities were assembled into three datasets (A, B and C). Dataset A (N = 213) included 50 K genotypic data for five taurine (Holstein-Friesian, Guernsey, Brown Swiss, Angus and Jersey) and two indicine (Gir and Nellore) breeds. Dataset B (N = 241) included the same breeds with HD density data. Dataset C (N = 299) included 50 K SNP genotypic data for six taurine (Holstein-Friesian, Jersey, Guernsey, Brown Swiss, Angus and Hereford) and six indicine (Hariana, Kankrej, Brahman, Nellore, Sahiwal and Gir) breeds. The analysis was done using ADMIXTURE program (bioinformatics-based) and cross-validation errors and Principal Component Analysis (statistical analysis). The proportion of polymorphic markers and minor allele frequencies were assessed for each breed. The proportion of markers polymorphic was consistently higher in taurine breeds when compared with breeds from indicine group. Minor allele frequency estimates and ADMIXTURE results showed differential patterns for both the lineages. However, no significant increase in the accuracy of genomic clustering was found on moving from 50 K to HD density data.

Genomic Selection and Use of Molecular Tools in Breeding Programs for Indigenous and Crossbred Cattle in Developing Countries: Current Status and Future Prospects

Genomic selection (GS) has resulted in rapid rates of genetic gains especially in dairy cattle in developed countries resulting in a higher proportion of genomically proven young bulls being used in breeding. This success has been undergirded by well-established conventional genetic evaluation systems. Here, the status of GS in terms of the structure of the reference and validation populations, response variables, genomic prediction models, validation methods, and imputation efficiency in breeding programs of developing countries, where smallholder systems predominate and the basic components for conventional breeding are mostly lacking is examined. Also, the application of genomic tools and identification of genome-wide signatures of selection is reviewed. The studies on genomic prediction in developing countries are mostly in dairy and beef cattle usually with small reference populations (500-3,000 animals) and are mostly cows. The input variables tended to be pre-corrected phenotypic records and the small reference populations has made implementation of various Bayesian methods feasible in addition to GBLUP. Multi-trait single-step has been used to incorporate genomic information from foreign bulls, thus GS in developing countries would benefit from collaborations with developed countries, as many dairy sires used are from developed countries where they may have been genotyped and phenotyped. Cross validation approaches have been implemented in most studies resulting in accuracies of 0.20-0.60. Genotyping animals with a mixture of HD and LD chips, followed by imputation to the HD have been implemented with imputation accuracies of 0.74-0.99 reported. This increases the prospects of reducing genotyping costs and hence the cost-effectiveness of GS. Next-generation sequencing and associated technologies have allowed the determination of breed composition, parent verification, genome diversity, and genome-wide selection sweeps. This information can be incorporated into breeding programs aiming to utilize GS. Cost-effective GS in beef cattle in developing countries may involve usage of reproductive technologies (AI and in-vitro fertilization) to efficiently propagate superior genetics from the genomics pipeline. For dairy cattle, sexed semen of genomically proven young bulls could substantially improve profitability thus increase prospects of small holder farmers buying-in into genomic breeding programs.

The role of molecular genetics in livestock production

Genetic variations that lead to easy-to-identify phenotypic changes have always been of interest to livestock breeders since domestication. Molecular genetics has opened up possibilities for identifying these variations and understanding their biological and population effects. Moreover, molecular genetics is part of the most diverse approaches and applications in animal production nowadays, including paternity testing, selection based on genetic variants, diagnostic of genetic diseases, reproductive biotechniques, fraud identification, differentiation of hybrids, parasite identification, genetic evaluation, diversity studies, and genome editing, among others. Therefore, the objective of this review was to describe the different applications of molecular genetics in livestock production, contextualising them with examples and highlighting the importance of the study of these topics and their applications.

SNP markers associated with body size and pelt length in American mink (Neovison vison)

Background

Identification of genes underlying production traits is a key aim of the mink research community. Recent availability of genomic tools have opened the possibility for faster genetic progress in mink breeding. Availability of mink genome assembly allows genome-wide association studies in mink.

Results

In this study, we used genotyping-by-sequencing to obtain single nucleotide polymorphism (SNP) genotypes of 2496 mink. After multiple rounds of filtering, we retained 28,336 high quality SNPs and 2352 individuals for a genome-wide association study (GWAS). We performed the first GWAS for body weight, behavior, along with 10 traits related to fur quality in mink.

Conclusions

Combining association results with existing functional information of genes and mammalian phenotype databases, we proposed WWC3, MAP2K4, SLC7A1 and USP22 as candidate genes for body weight and pelt length in mink.

Electronic supplementary material

The online version of this article (10.1186/s12863-018-0688-6) contains supplementary material, which is available to authorized users.

Genetic Improvement in South African Livestock: Can Genomics Bridge the Gap Between the Developed and Developing Sectors?

South Africa (SA) holds a unique position on the African continent with a rich diversity in terms of available livestock resources, vegetation, climatic regions and cultures. The livestock sector has been characterized by a dual system of a highly developed commercial sector using modern technology vs. a developing sector including emerging and smallholder farmers. Emerging farmers typically aim to join the commercial sector, but lag behind with regard to the use of modern genetic technologies, while smallholder farmers use traditional practices aimed at subsistence. Several factors influence potential application of genomics by the livestock industries, which include available research funding, socio-economic constraints and extension services. State funded Beef and Dairy genomic programs have been established with the aim of building reference populations for genomic selection with most of the potential beneficiaries in the well-developed commercial sector. The structure of the beef, dairy and small stock industries is fragmented and the outcomes of selection strategies are not perceived as an advantage by the processing industry or the consumer. The indigenous and local composites represent approximately 40% of the total beef and sheep populations and present valuable genetic resources. Genomic research has mostly provided insight on genetic biodiversity of these resources, with limited attention to novel phenotypes associated with adaptation or disease tolerance. Genetic improvement of livestock through genomic technology needs to address the role of adapted breeds in challenging environments, increasing reproductive and growth efficiency. National animal recording schemes contributed significantly to progress in the developed sector with regard to genetic evaluations and estimated breeding values (EBV) as a selection tool over the past three decades. The challenge remains on moving the focus to novel traits for increasing efficiency and addressing welfare and environmental issues. Genetic research programs are required that will be directed to bridge the gap between the elite breeders and the developing livestock sector. The aim of this review was to provide a perspective on the dichotomy in the South African livestock sector arguing that a realistic approach to the use of genomics in beef, dairy and small stock is required to ensure sustainable long term genetic progress.

Whole-genome SNP analysis elucidates the genetic structure of Russian cattle and its relationship with Eurasian taurine breeds

BACKGROUND

The origin of native and locally developed Russian cattle breeds is linked to the historical, social, cultural, and climatic features of the diverse geographical regions of Russia. In the present study, we investigated the population structure of nine Russian cattle breeds and their relations to the cattle breeds from around the world to elucidate their origin. Genotyping of single nucleotide polymorphisms (SNPs) in Bestuzhev (n = 26), Russian Black-and-White (n = 21), Kalmyk (n = 14), Kholmogor (n = 25), Kostromsky (n = 20), Red Gorbatov (n = 23), Suksun (n = 20), Yakut (n = 25), and Yaroslavl cattle breeds (n = 21) was done using the Bovine SNP50 BeadChip. SNP profiles from an additional 70 breeds were included in the analysis as references.

RESULTS

The observed heterozygosity levels were quite similar in eight of the nine studied breeds (H_O = 0.337-0.363) except for Yakut (Ho = 0.279). The inbreeding coefficients F_IS ranged from -0.028 for Kalmyk to 0.036 for Russian Black-and-White and were comparable to those of the European breeds. The nine studied Russian breeds exhibited taurine ancestry along the C1 axis of the multidimensional scaling (MDS)-plot, but Yakut was clearly separated from the European taurine breeds on the C2 axis. Neighbor-Net and admixture analyses, discriminated three groups among the studied Russian breeds. Yakut and Kalmyk were assigned to a separate group because of their Turano-Mongolian origin. Russian Black-and-White, Kostromsky and Suksun showed transboundary European ancestry, which originated from the Holstein, Brown Swiss, and Danish Red breeds, respectively. The lowest level of introgression of transboundary breeds was recorded for the Kholmogor, Yaroslavl, Red Gorbatov and Bestuzhev breeds, which can be considered as an authentic genetic resource.

CONCLUSIONS

Whole-genome SNP analysis revealed that Russian native and locally developed breeds have conserved authentic genetic patterns in spite of the considerable influence of Eurasian taurine cattle. In this paper, we provide fundamental genomic information that will contribute to the development of more accurate breed conservation programs and genetic improvement strategies.

Conservation status and historical relatedness of Italian cattle breeds

Background

In the last 50 years, the diversity of cattle breeds has experienced a severe contraction. However, in spite of the growing diffusion of cosmopolite specialized breeds, several local cattle breeds are still farmed in Italy. Genetic characterization of breeds represents an essential step to guide decisions in the management of farm animal genetic resources. The aim of this work was to provide a high-resolution representation of the genome-wide diversity and population structure of Italian local cattle breeds using a medium-density single nucleotide polymorphism (SNP) array.

Results

After quality control filtering, the dataset included 31,013 SNPs for 800 samples from 32 breeds. Our results on the genetic diversity of these breeds agree largely with their recorded history. We observed a low level of genetic diversity, which together with the small size of the effective populations, confirmed that several breeds are threatened with extinction. According to the analysis of runs of homozygosity, evidence of recent inbreeding was strong in some local breeds, such as Garfagnina, Mucca Pisana and Pontremolese. Patterns of genetic differentiation, shared ancestry, admixture events, and the phylogenetic tree, all suggest the presence of gene flow, in particular among breeds that originate from the same geographical area, such as the Sicilian breeds. In spite of the complex admixture events that most Italian cattle breeds have experienced, they have preserved distinctive characteristics and can be clearly discriminated, which is probably due to differences in genetic origin, environment, genetic isolation and inbreeding.

Conclusions

This study is the first exhaustive genome-wide analysis of the diversity of Italian cattle breeds. The results are of significant importance because they will help design and implement conservation strategies. Indeed, efforts to maintain genetic diversity in these breeds are needed. Improvement of systems to record and monitor inbreeding in these breeds may contribute to their in situ conservation and, in view of this, the availability of genomic data is a fundamental resource.

Electronic supplementary material

The online version of this article (10.1186/s12711-018-0406-x) contains supplementary material, which is available to authorized users.

Genetic Diversity of Seven Cattle Breeds Inferred Using Copy Number Variations

Copy number variations (CNVs) comprise deletions, duplications, and insertions found within the genome larger than 50 bp in size. CNVs are thought to be primary role-players in breed formation and adaptation. South Africa boasts a diverse ecology with harsh environmental conditions and a broad spectrum of parasites and diseases that pose challenges to livestock production. This has led to the development of composite cattle breeds which combine the hardiness of Sanga breeds and the production potential of the Taurine breeds. The prevalence of CNVs within these respective breeds of cattle and the prevalence of CNV regions (CNVRs) in their diversity, adaptation and production is however not understood. This study therefore aimed to ascertain the prevalence, diversity, and correlations of CNVRs within cattle breeds used in South Africa. Illumina Bovine SNP50 data and PennCNV were utilized to identify CNVRs within the genome of 287 animals from seven cattle breeds representing Sanga, Taurine, Composite, and cross breeds. Three hundred and fifty six CNVRs of between 36 kb to 4.1 Mb in size were identified. The null hypothesis that one CNVR loci is independent of another was tested using the GENEPOP software. One hunded and two and seven of the CNVRs in the Taurine and Sanga/Composite cattle breeds demonstrated a significant (p ≤ 0.05) association. PANTHER overrepresentation analyses of correlated CNVRs demonstrated significant enrichment of a number of biological processes, molecular functions, cellular components, and protein classes. CNVR genetic variation between and within breed group was measured using phiPT which allows intra-individual variation to be suppressed and hence proved suitable for measuring binary CNVR presence/absence data. Estimate PhiPT within and between breed variance was 2.722 and 0.518 respectively. Pairwise population PhiPT values corresponded with breed type, with Taurine Holstein and Angus breeds demonstrating no between breed CNVR variation. Phylogenetic trees were drawn. CNVRs primarily clustered animals of the same breed type together. This study successfully identified, characterized, and analyzed 356 CNVRs within seven cattle breeds. CNVR correlations were evident, with many more correlations being present among the exotic Taurine breeds. CNVR genetic diversity of Sanga, Taurine and Composite breeds was ascertained with breed types exposed to similar selection pressures demonstrating analogous incidences of CNVRs.

Efficiency of different strategies to mitigate ascertainment bias when using SNP panels in diversity studies

Background

Single nucleotide polymorphism (SNP) panels have been widely used to study genomic variations within and between populations. Methods of SNP discovery have been a matter of debate for their potential of introducing ascertainment bias, and genetic diversity results obtained from the SNP genotype data can be misleading. We used a total of 42 chicken populations where both individual genotyped array data and pool whole genome resequencing (WGS) data were available. We compared allele frequency distributions and genetic diversity measures (expected heterozygosity (H_e), fixation index (F_ST) values, genetic distances and principal components analysis (PCA)) between the two data types. With the array data, we applied different filtering options (SNPs polymorphic in samples of two Gallus gallus wild populations, linkage disequilibrium (LD) based pruning and minor allele frequency (MAF) filtering, and combinations thereof) to assess their potential to mitigate the ascertainment bias.

Results

Rare SNPs were underrepresented in the array data. Array data consistently overestimated H_e compared to WGS data, however, with a similar ranking of the breeds, as demonstrated by Spearman’s rank correlations ranging between 0.956 and 0.985. LD based pruning resulted in a reduced overestimation of H_e compared to the other filters and slightly improved the relationship with the WGS results. The raw array data and those with polymorphic SNPs in the wild samples underestimated pairwise F_ST values between breeds which had low F_ST (<0.15) in the WGS, and overestimated this parameter for high WGS F_ST (>0.15). LD based pruned data underestimated F_ST in a consistent manner. The genetic distance matrix from LD pruned data was more closely related to that of WGS than the other array versions. PCA was rather robust in all array versions, since the population structure on the PCA plot was generally well captured in comparison to the WGS data.

Conclusions

Among the tested filtering strategies, LD based pruning was found to account for the effects of ascertainment bias in the relatively best way, producing results which are most comparable to those obtained from WGS data and therefore is recommended for practical use.

Electronic supplementary material

The online version of this article (doi: 10.1186/s12864-017-4416-9) contains supplementary material, which is available to authorized users.

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.

Evaluation of Bovine High-Density SNP Genotyping Array in Indigenous Dairy Cattle Breeds

In total 52 samples of Sahiwal ( 19 ), Tharparkar ( 17 ), and Gir ( 16 ) were genotyped by using BovineHD SNP chip to analyze minor allele frequency (MAF), genetic diversity, and linkage disequilibrium among these cattle. The common SNPs of BovineHD and 54K SNP Chips were also extracted and evaluated for their performance. Only 40%-50% SNPs of these arrays was found informative for genetic analysis in these cattle breeds. The overall mean of MAF for SNPs of BovineHD SNPChip was 0.248 ± 0.006, 0.241 ± 0.007, and 0.242 ± 0.009 in Sahiwal, Tharparkar and Gir, respectively, while that for 54K SNPs was on lower side. The average Reynold's genetic distance between breeds ranged from 0.042 to 0.055 based on BovineHD Beadchip, and from 0.052 to 0.084 based on 54K SNP Chip. The estimates of genetic diversity based on HD and 54K chips were almost same and, hence, low density chip seems to be good enough to decipher genetic diversity of these cattle breeds. The linkage disequilibrium started decaying (r² < 0.2) at 140 kb inter-marker distance and, hence, a 20K low density customized SNP array from HD chip could be designed for genomic selection in these cattle else the 54K Bead Chip as such will be useful.

Insight into the genetic composition of South African Sanga cattle using SNP data from cattle breeds worldwide

Background

Understanding the history of cattle breeds is important because it provides the basis for developing appropriate selection and breed improvement programs. In this study, patterns of ancestry and admixture in Afrikaner, Nguni, Drakensberger and Bonsmara cattle of South Africa were investigated. We used 50 K single nucleotide polymorphism genotypes that were previously generated for the Afrikaner (n = 36), Nguni (n = 50), Drakensberger (n = 47) and Bonsmara (n = 44) breeds, and for 394 reference animals representing European taurine, African taurine, African zebu and Bos indicus.

Results and discussion

Our findings support previous conclusions that Sanga cattle breeds are composites between African taurine and Bos indicus. Among these breeds, the Afrikaner breed has significantly diverged from its ancestral forebears, probably due to genetic drift and selection to meet breeding objectives of the breed society that enable registration. The Nguni, Drakensberger and Bonsmara breeds are admixed, perhaps unintentionally in the case of Nguni and Drakensberger, but certainly by design in the case of Bonsmara, which was developed through crossbreeding between the Afrikaner, Hereford and Shorthorn breeds.

Conclusions

We established patterns of admixture and ancestry for South African Sanga cattle breeds, which provide a basis for developing appropriate strategies for their genetic improvement.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-016-0266-1) contains supplementary material, which is available to authorized users.