Assessment of Genetic Diversity, Relationships and Structure among Korean Native Cattle Breeds Using Microsatellite Markers

Effects of diets for three growing stages by rumen inocula donors on in vitro rumen fermentation and microbiome

Hanwoo and Jeju Black cattle (Jeju Black) are native breeds of Korean cattle. Jeju Black cattle are recognized as natural monuments and are known to exhibit slower growth rates compared to Hanwoo. While several studies have analyzed the genetic characteristics of these cattle, there has been limited research on the differences in their microbiome. In this study, rumen fluid was obtained from three Hanwoo steers and three Jeju Black steers, and three different diets (total mixed rations [TMRs] for growing, early fattening, and late fattening periods) were used as substrates for in vitro fermentation. The in vitro incubation was conducted for 3 h and 24 h following a 2 × 3 factorial arrangement. After both incubation periods, fermentation characteristics were analyzed, and ruminal microbiome analysis was performed using 16S rRNA gene sequencing, employing both QIIME2 and PICRUSt2. The results revealed significant differences in the ruminal microbiota due to the inoculum effect. At the phylum level, Patescibacteria and Synergistota were found to be enriched in the Jeju Black inoculum-treated group. Additionally, using different inocula also affected the relative abundance of major taxa, including Ruminococcus, Pseudoramibacter, Ruminococcaceae CAG-352, and the [Eubacterium] ruminantium group. These microbial differences induced by the inoculum may have originated from varying levels of domestication between the two subspecies of donor animals, which mainly influenced the fermentation and microbiome features in the early incubation stages, although this was only partially offset afterward. Furthermore, predicted commission numbers of microbial enzymes, some of which are involved in the biosynthesis of secondary metabolites, fatty acids, and alpha amylase, differed based on the inoculum effect. However, these differences may account for only a small proportion of the overall metabolic pathway. Conversely, diets were found to affect protein biosynthesis and its related metabolism, which showed differential abundance in the growing diet and were potentially linked to the growth-promoting effects in beef cattle during the growing period. In conclusion, this study demonstrated that using different inocula significantly affected in vitro fermentation characteristics and microbiome features, mainly in the early stages of incubation, with some effects persisting up to 24 h of incubation.

Genetic diversity of bovine populations raised in Senegal

Background

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

Objective

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

Material and methods

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

Results

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

Conclusion

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

High-Level Gene Flow Restricts Genetic Differentiation in Dairy Cattle Populations in Thailand: Insights from Large-Scale Mt D-Loop Sequencing

Domestication and artificial selection lead to the development of genetically divergent cattle breeds or hybrids that exhibit specific patterns of genetic diversity and population structure. Recently developed mitochondrial markers have allowed investigation of cattle diversity worldwide; however, an extensive study on the population-level genetic diversity and demography of dairy cattle in Thailand is still needed. Mitochondrial D-loop sequences were obtained from 179 individuals (hybrids of Bos taurus and B. indicus) sampled from nine different provinces. Fifty-one haplotypes, of which most were classified in haplogroup "I", were found across all nine populations. All sampled populations showed severely reduced degrees of genetic differentiation, and low nucleotide diversity was observed in populations from central Thailand. Populations that originated from adjacent geographical areas tended to show high gene flow, as revealed by patterns of weak network structuring. Mismatch distribution analysis was suggestive of a stable population, with the recent occurrence of a slight expansion event. The results provide insights into the origins and the genetic relationships among local Thai cattle breeds and will be useful for guiding management of cattle breeding in Thailand.

High-throughput genotype-based population structure analysis of selected buffalo breeds

India is considered as the home tract of some of the best buffalo breeds. However, the genetic structure of the Indian river buffalo is poorly understood. Hence, there is a need to characterize the populations and understand the genetic structure of various buffalo breeds for selection and to design breeding strategies. In this study, we have analyzed genetic variability and population structure of seven buffalo breeds from their respective geographical regions using Axiom Buffalo Genotyping Array. Diversity, as measured by expected heterozygosity, ranged from 0.364 in Surti to 0.384 in Murrah breed, and pair-wise F_ST values revealed the lowest genetic distance between Murrah and Nili-Ravi (0.0022), while the highest between Surti and Pandharpuri (0.030). Principal component analysis and structure analysis unveiled the differentiation of Surti, Pandharpuri, and Jaffarabadi in first two principal components and at K = 4, respectively, while remaining breeds were grouped together as a separate single cluster and admixed. Murrah and Mehsana showed early linkage disequilibrium (LD) decay, while Surti breed showed late decay. In LD blocks to quantitative trait locis (QTLs) concordance analysis, 4.65% of concordance was observed with 873 LD blocks overlapped with 2,330 QTLs. Overall, total 4,090 markers were identified from all LD blocks for six types of traits. Results of this study indicated that these single-nucleotide polymorphism (SNP) markers could differentiate phenotypically distinct breeds like Surti, Pandharpuri, and Jaffarabadi but not others. So, there is a need to develop SNP chip based on SNP markers identified by sequence information of local breeds.

Genetic characterization and population structure of different coat colour variants of Badri cattle

The present study aimed to genetically characterize the Badri cattle and its three colour variants and assess their population structure using 24 microsatellite markers. Out of 96 animals analyzed, 32 each were collected from grey (GVBC), brown (BrVBC) and black (BVBC) colour variants of Badri cattle (BC). The genetic diversity parameters including allele frequencies, observed and effective number of alleles, observed and expected heterozygosity, PIC, Shannon's indices and F-statistics were estimated using POPGENE software. Bottleneck analysis was performed using both qualitative and quantitative approaches. A total of 274 alleles (50 private and 224 shared) were scored for BC, GVBC, BrVBC and BVBC with mean number of 11.417, 9.083, 9.125 and 9.083 alleles, respectively. All populations exhibited average heterozygosity estimate > 0.5 indicating existence of substantial genetic variability, concurrent with revelations from Shannon's indices. Observed mean PIC estimates (> 0.74) were indicative of optimum informativeness of used microsatellite markers. The mean inbreeding estimates (F) in GVBC, BrVBC and BVBC were 0.041, - 0.024 and 0.016, respectively. The pair wise genetic (> 0.91) pointed towards similarity between different colour variant populations. STRUCTURE analysis also revealed clear admixture for the three Badri colour variants indicating absence of genetic differentiation. The present study revealed first-hand information that populations of Badri cattle with different phenotypes with respect to coat colour are genetically related and can be considered as a single breed. The comprehensive knowledge generated for Badri cattle will help in designing breeding plan for its genetic improvement and deciding the conservation priorities.

Genetic Characterization of the Local Pirenaica Cattle for Parentage and Traceability Purposes

Pirenaica is the most important autochthonous cattle breed within the Protected Geographic Indication (PGI) beef quality label in the Basque region, in northern Spain. The short tandem repeats (STRs) are powerful markers to elucidate forensic cases and traceability across the agri-food sector. The main objective of the present work was to study the phylogenetic relationships of Pirenaica cattle and other breeds typically raised in the region and provide the minimum number of STR markers for parentage and traceability purposes. The 30-STR panel recommended by the International Society of Animal Genetics-Food and Agriculture Organization of the United Nations (ISAG-FAO) was compared against other commercial STR panels. The 30-STR panel showed a combined matching probability of 1.89 × 10-25 and a power of exclusion for duos of 0.99998. However, commercial STR panels showed a limited efficiency for a reliable parentage analysis in Pirenaica, and at least a 21-STR panel is needed to reach a power of exclusion of 0.9999. Machine-learning analysis also demonstrated a 95% accuracy in assignments selecting the markers with the highest FST in Pirenaica individuals. Overall, the present study shows the genetic characterization of Pirenaica and its phylogeny compared with other breeds typically raised in the Basque region. Finally, a 21-STR panel with the highest FST markers is proposed for a confident parentage analysis and high traceability.

Microsatellite Diversity and Phylogenetic Relationships among East Eurasian Bos taurus Breeds with an Emphasis on Rare and Ancient Local Cattle

We report the genetic analysis of 18 population samples of animals, which were taken from cattle (Bos taurus) breeds of European and Asian origins. The main strength of our study is the use of rare and ancient native cattle breeds: the Altai, Ukrainian Grey, Tagil, and Buryat ones. The cattle samples studied have different production purposes, belong to various eco-geographic regions, and consequently have distinct farming conditions. In order to clarify the genetic diversity, phylogenetic relationships and historical origin of the studied breeds, we carried out an analysis of the genetic variation of 14 high-variability microsatellite loci at 1168 genotyped animals. High levels of heterozygosity and allelic richness were identified in four of the ancient local breeds, namely the Kalmyk, Tagil, Kyrgyz native, and Buryat breeds. The greatest phylogenetic distances from a common ancestor were observed for the Yakut and Ukrainian Grey breeds, while the Tagil breed showed the smallest difference. By using clustering approaches, we found that the Altai cattle is genetically close to the Kyrgyz one. Moreover, both the Altai and Kyrgyz breeds exposed genetic divergences from other representatives of the Turano-Mongolian type and genetic relationships with the Brown Swiss and Kostroma breeds. This phenomenon can be explained by the extensive use of the Brown Swiss and Kostroma breeds in the breeding and improvement processes for the Kyrgyz breeds, which have been involved in the process of keeping the Altai cattle. Our results can be valuable for conservation and management purposes.

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.

Estimation of environmental effects and genetic parameters of carcass traits on Chikso (Korean brindle cattle)

OBJECTIVE

The purpose of this study was i) to identify the characteristics of carcass traits in Chikso by gender, region, age at slaughter, and coat color using the carcass data collected from the nationwide pedigree information and coat color investigation, and ii) to estimate genetic parameters for breed improvement.

METHODS

A linear model was used to analyze the environmental effects on the carcass traits and to estimate genetic parameters. Analysis of variance was performed using TYPE III sum of squares for the unbalanced data provided by the general linear model procedure. Variance components for genetic parameters was estimated using REMLF90 of the BLUPF90 family programs.

RESULTS

Phenotypic performance of carcass weight (CW), eye muscle area (EMA), and backfat thickness (BF) in Chikso were lower than those of Hanwoo. This is a natural outcome because Hanwoo have undergone significant efforts for improvement at the national level, a phenomenon not observed in Chikso. Another factor influencing the above outcome was the smaller population size of Chikso compared to that of Hanwoo's. The heritabilities of CW, EMA, BF, and marbling score in Chikso were estimated as 0.50, 0.37, 0.35, and 0.53, respectively, which were was higher than those of Hanwoo.

CONCLUSION

Based on the genetic parameters that were estimated in this study, it is expected that the carcass traits will improve when the livestock research institutes at each province conduct small-scale performance tests and the semen is provided to farmers after selecting proven bulls using the state-of-art selection technique such as genomic selection.

Genome-wide linkage disequilibrium and past effective population size in three Korean cattle breeds

The routine collection and use of genomic data are useful for effectively managing breeding programs for endangered populations. Linkage disequilibrium (LD) using high-density DNA markers has been widely used to determine population structures and predict the genomic regions that are associated with economic traits in beef cattle. The extent of LD also provides information about historical events, including past effective population size (N_e ), and it allows inferences on the genetic diversity of breeds. The objective of this study was to estimate the LD and N_e in three Korean cattle breeds that are genetically similar but have different coat colors (Brown, Brindle and Jeju Black Hanwoo). Brindle and Jeju Black are endangered breeds with small populations, whereas Brown Hanwoo is the main breeding population in Korea. DNA samples from these cattle breeds were genotyped using the Illumina BovineSNP50 Bead Chip. We examined 13 cattle breeds, including European taurines, African taurines and indicines, and hybrids to compare their LD values. Brown Hanwoo consistently had the lowest mean LD compared to Jeju Black, Brindle and the other 13 cattle breeds (0.13, 0.19, 0.21 and 0.15-0.22 respectively). The high LD values of Brindle and Jeju Black contributed to small N_e values (53 and 60 respectively), which were distinct from that of Brown Hanwoo (531) for 11 generations ago. The differences in LD and N_e for each breed reflect the breeding strategy applied. The N_e for these endangered cattle breeds remain low; thus, effort is needed to bring them back to a sustainable tract.

Haplogroup Classification of Korean Cattle Breeds Based on Sequence Variations of mtDNA Control Region

Many studies have reported the frequency and distribution of haplogroups among various cattle breeds for verification of their origins and genetic diversity. In this study, 318 complete sequences of the mtDNA control region from four Korean cattle breeds were used for haplogroup classification. 71 polymorphic sites and 66 haplotypes were found in these sequences. Consistent with the genetic patterns in previous reports, four haplogroups (T1, T2, T3, and T4) were identified in Korean cattle breeds. In addition, T1a, T3a, and T3b sub-haplogroups were classified. In the phylogenetic tree, each haplogroup formed an independent cluster. The frequencies of T3, T4, T1 (containing T1a), and T2 were 66%, 16%, 10%, and 8%, respectively. Especially, the T1 haplogroup contained only one haplotype and a sample. All four haplogroups were found in Chikso, Jeju black and Hanwoo. However, only the T3 and T4 haplogroups appeared in Heugu, and most Chikso populations showed a partial of four haplogroups. These results will be useful for stable conservation and efficient management of Korean cattle breeds.

Study of Genetic Diversity among Simmental Cross Cattle in West Sumatra Based on Microsatellite Markers

A study was conducted to assess the genetic diversity among Simmental Cross cattle in West Sumatra using microsatellite DNA markers. A total of 176 individual cattle blood samples was used for obtaining DNA samples. Twelve primers of microsatellite loci as recommended by FAO were used to identify the genetic diversity of the Simmental Cross cattle population. Multiplex DNA fragment analysis method was used for allele identification. All the microsatellite loci in this study were highly polymorphic and all of the identified alleles were able to classify the cattle population into several groups based on their genetic distance. The heterozygosity values of microsatellite loci in this study ranged from 0.556 to 0.782. The polymorphism information content (PIC) value of the 12 observed loci is high (PIC>0.5). The highest PIC value in the Simmental cattle population was 0.893 (locus TGLA53), while the lowest value was 0.529 (locus BM1818). Based on the genetic distance value, the subpopulation of the Simmental Cross-Agam and the Simmental Cross-Limapuluh Kota was exceptionally close to the Simmental Purebred thus indicating that a grading-up process has taken place with the Simmental Purebred. In view of the advantages possessed by the Simmental Cross cattle and the evaluation of the genetic diversity results, a number of subpopulations in this study can be considered as the initial (base) population for the Simmental Cross cattle breeding programs in West Sumatra, Indonesia.

Towards breed formation by island model divergence in Korean cattle

Background

The main cattle breed in Korea is the brown Hanwoo, which has been under artificial selection within a national breeding program for several decades. Varieties of the Hanwoo known as Jeju Black and Chikso were not included in the breeding program and remained isolated from the effects of recent artificial selection advancements. We analysed the Jeju Black and Chikso populations in regards to their genetic variability, state of inbreeding, as well as level of differentiation from the mainland Hanwoo population.

Results

Jeju Black and Chikso were found to have small estimated effective population sizes (N_e) of only 11 and 7, respectively. Despite a small N_e, higher than expected heterozygosity levels were observed (0.303 and 0.306), however, lower allelic richness was found for the two island populations (1.76 and 1.77) compared to the mainland population (1.81). The increase in heterozygosity could be due to environmental disease challenges that promoted maintenance of higher genetic variability; however, no direct proof exists. Increased heterozygosity due to a first generation crossing of genetically different populations is not recorded. The differentiation between the Korean populations had F_ST values between 0.014 and 0.036 which is not as high as the differentiation within European beef or dairy cattle breeds (0.047–0.111). This suggests that the three populations have not separated into independent breeds.

Conclusion

Results agree with an island model of speciation where the brown Hanwoo represents the ancestral breed, whilst the Jeju Black and Chikso diverge from this common ancestor, following different evolutionary trajectories. Nevertheless, differences are minor and whether Jeju Black and Chikso cattle will develop into discrete breeds or reintegrate with the main population has to be seen in the future and will largely depend on human management decisions. This offers a rare opportunity to accompany the development of new breeds but also poses challenges on how to preserve these incipient breeds and ensure their long term viability.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0563-2) contains supplementary material, which is available to authorized users.

Genetic analysis of Thai cattle reveals a Southeast Asian indicine ancestry

Cattle commonly raised in Thailand have characteristics of Bos indicus (zebu). We do not know when or how cattle domestication in Thailand occurred, and so questions remain regarding their origins and relationships to other breeds. We obtained genome-wide SNP genotypic data of 28 bovine individuals sampled from four regions: North (Kho-Khaolampoon), Northeast (Kho-Isaan), Central (Kho-Lan) and South (Kho-Chon) Thailand. These regional varieties have distinctive traits suggestive of breed-like genetic variations. From these data, we confirmed that all four Thai varieties are Bos indicus and that they are distinct from other indicine breeds. Among these Thai cattle, a distinctive ancestry pattern is apparent, which is the purest within Kho-Chon individuals. This ancestral component is only present outside of Thailand among other indicine breeds in Southeast Asia. From this pattern, we conclude that a unique Bos indicus ancestor originated in Southeast Asia, and native Kho-Chon Thai cattle retain the signal of this ancestry with limited admixture of other bovine ancestors.