Whole-Genome Analyses of Korean Native and Holstein Cattle Breeds by Massively Parallel Sequencing

A main goal of cattle genomics is to identify DNA differences that account for variations in economically important traits. In this study, we performed whole-genome analyses of three important cattle breeds in Korea—Hanwoo, Jeju Heugu, and Korean Holstein—using the Illumina HiSeq 2000 sequencing platform. We achieved 25.5-, 29.6-, and 29.5-fold coverage of the Hanwoo, Jeju Heugu, and Korean Holstein genomes, respectively, and identified a total of 10.4 million single nucleotide polymorphisms (SNPs), of which 54.12% were found to be novel. We also detected 1,063,267 insertions–deletions (InDels) across the genomes (78.92% novel). Annotations of the datasets identified a total of 31,503 nonsynonymous SNPs and 859 frameshift InDels that could affect phenotypic variations in traits of interest. Furthermore, genome-wide copy number variation regions (CNVRs) were detected by comparing the Hanwoo, Jeju Heugu, and previously published Chikso genomes against that of Korean Holstein. A total of 992, 284, and 1881 CNVRs, respectively, were detected throughout the genome. Moreover, 53, 65, 45, and 82 putative regions of homozygosity (ROH) were identified in Hanwoo, Jeju Heugu, Chikso, and Korean Holstein respectively. The results of this study provide a valuable foundation for further investigations to dissect the molecular mechanisms underlying variation in economically important traits in cattle and to develop genetic markers for use in cattle breeding.     

Genome-wide analyses of the Jeju,Thoroughbred, and Jeju crossbred horse populations using the high density SNP array

The Jeju horse is an indigenous Korean horse breed that is currently registered with the Food and Agriculture Organization of the United Nations. However, there is severe lack of genomic studies on Jeju horse. This study was conducted to investigate genetic characteristics of horses including Jeju horse, Thoroughbred and Jeju crossbred (Jeju?×?Thoroughbred) populations. We compared the genomes of three horse populations using the Equine SNP70 Beadchip array. Short-range Linkage disequilibrium was the highest in Thoroughbred, whereas r² values were lowest in Jeju horse. Expected heterozygosity was the highest in Jeju crossbred (0.351), followed by the Thoroughbred (0.337) and Jeju horse (0.311). The level of inbreeding was slightly higher in Thoroughbred (??0.009) than in Jeju crossbred (??0.035) and Jeju horse (??0.038). F_ST value was the highest between Jeju horse and Thoroughbred (0.113), whereas Jeju crossbred and Thoroughbred showed the lowest value (0.031). The genetic relationship was further assessed by principal component analysis, suggesting that Jeju crossbred is more genetically similar to Thoroughbred than Jeju horse population. Additionally, we detected potential selection signatures, for example, in loci located on LCORL/NCAPG and PROP1 genes that are known to influence body. Genome-wide analyses of the three horse populations showed that all the breeds had somewhat a low level of inbreeding within each population. In the population structure analysis, we found that Jeju crossbred was genetically closer to Thoroughbred than Jeju horse. Furthermore, we identified several signatures of selection which might be associated with traits of interest. To our current knowledge, this study is the first genomic research, analyzing genetic relationships of Jeju horse, Thoroughbred and Jeju crossbred.     

Performance of different SNP panels for parentage testing in two East Asian cattle breeds

The International Society for Animal Genetics (ISAG) proposed a panel of single nucleotide polymorphisms (SNPs) for parentage testing in cattle (a core panel of 100 SNPs and an additional list of 100 SNPs). However, markers specific to East Asian taurine cattle breeds were not included, and no information is available as to whether the ISAG panel performs adequately for these breeds. We tested ISAG's core (100 SNP) and full (200 SNP) panels on two East Asian taurine breeds: the Korean Hanwoo and the Japanese Wagyu, the latter from the Australian herd. Even though the power of exclusion was high at 0.99 for both ISAG panels, the core panel performed poorly with 3.01% false‐positive assignments in the Hanwoo population and 3.57% in the Wagyu. The full ISAG panel identified all sire–offspring relations correctly in both populations with 0.02% of relations wrongly excluded in the Hanwoo population. Based on these results, we created and tested two population‐specific marker panels: one for the Wagyu population, which showed no false‐positive assignments with either 100 or 200 SNPs, and a second panel for the Hanwoo, which still had some false‐positive assignments with 100 SNPs but no false positives using 200 SNPs. In conclusion, for parentage assignment in East Asian cattle breeds, only the full ISAG panel is adequate for parentage testing. If fewer markers should be used, it is advisable to use population‐specific markers rather than the ISAG panel.     

Characterization of experimental cattle populations from Argentina with a low-density SNP genotyping panel

The objectives of the present experiment were to evaluate a low-density SNP array designed for the molecular characterisation of gene banks and to assess the genetic diversity and population structure of beef cattle herds from an Argentinean research station. Forty-nine animals from three breeds (Angus, Hereford, and Argentinean Creole) were genotyped using the multi-species IMAGE001 60-K SNP array (10 K for cattle). Genotypes of other 19 cattle populations from Argentina, other American countries, and Europe were included in the study. Of special interest was the characterization of the Argentinean Creole, the only autochthonous cattle breed in the country. Due to the merging of different datasets, approximately 5 K SNPs were effectively used. Genetic differentiation (F_ST), principal component analysis, neighbour-joining tree of Reynolds distances and ancestry analysis showed that autochthonous American breeds are clearly differentiated, but all have genetic influences of Iberian cattle. The analysed herds of Argentinean Creole showed no evidence of recent admixture and represent a unique genetic pool within local American breeds. An experimental herd and the local commercial Hereford population have also diverged, probably due to the influence of current selection objectives in the breed. Our results illustrate the utility of using low-cost, low density SNP arrays in the evaluation of animal genetic resources. This type of panels could become a very useful resource in developing countries, where most endangered cattle breeds are located. The results also reinforce the importance of experimental herds as reservoir of genetic diversity, particularly in the case of local breeds under-represented in traditional production systems.     

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.     

Signatures of selection in five Italian cattle breeds detected by a 54K SNP panel

In this study we used a medium density panel of SNP markers to perform population genetic analysis in five Italian cattle breeds. The BovineSNP50 BeadChip was used to genotype a total of 2,935 bulls of Piedmontese, Marchigiana, Italian Holstein, Italian Brown and Italian Pezzata Rossa breeds. To determine a genome-wide pattern of positive selection we mapped the F _st values against genome location. The highest F _st peaks were obtained on BTA6 and BTA13 where some candidate genes are located. We identified selection signatures peculiar of each breed which suggest selection for genes involved in milk or meat traits. The genetic structure was investigated by using a multidimensional scaling of the genetic distance matrix and a Bayesian approach implemented in the STRUCTURE software. The genotyping data showed a clear partitioning of the cattle genetic diversity into distinct breeds if a number of clusters equal to the number of populations were given. Assuming a lower number of clusters beef breeds group together. Both methods showed all five breeds separated in well defined clusters and the Bayesian approach assigned individuals to the breed of origin. The work is of interest not only because it enriches the knowledge on the process of evolution but also because the results generated could have implications for selective breeding programs.     

Genome‐wide detection of signatures of selection in Korean Hanwoo cattle

The Korean Hanwoo cattle have been intensively selected for production traits, especially high intramuscular fat content. It is believed that ancient crossings between different breeds contributed to forming the Hanwoo, but little is known about the genomic differences and similarities between other cattle breeds and the Hanwoo. In this work, cattle breeds were grouped by origin into four types and used for comparisons: the Europeans (represented by six breeds), zebu (Nelore), African taurine (N'Dama) and Hanwoo. All animals had genotypes for around 680 000 SNPs after quality control of genotypes. Average heterozygosity was lower in Nelore and N'Dama (0.22 and 0.21 respectively) than in Europeans (0.26–0.31, with Shorthorn as outlier at 0.24) and Hanwoo (0.29). Pairwise F_ST analyses demonstrated that Hanwoo are more related to European cattle than to Nelore, with N'Dama in an intermediate position. This finding was corroborated by principal components and unsupervised hierarchical clustering. Using genome‐wide smoothed F_ST, 55 genomic regions potentially under positive selection in Hanwoo were identified. Among these, 29 were regions also detected in previous studies. Twenty‐four regions were exclusive to Hanwoo, and a number of other regions were shared with one or two of the other groups. These regions overlap a number of genes that are related to immune, reproduction and fatty acid metabolism pathways. Further analyses are needed to better characterize the ancestry of the Hanwoo cattle and to define the genes responsible to the identified selection peaks.     

Genomic Diversity in Pig (Sus scrofa) and its Comparison with Human and other Livestock

We have reviewed the current pig (Sus scrofa) genomic diversity within and between sites and compared them with human and other livestock. The current Porcine 60K single nucleotide polymorphism (SNP) panel has an average SNP distance in a range of 30 - 40 kb. Most of genetic variation was distributed within populations, and only a small proportion of them existed between populations. The average heterozygosity was lower in pig than in human and other livestock. Genetic inbreeding coefficient (F(IS)), population differentiation (F(ST)), and Nei's genetic distance between populations were much larger in pig than in human and other livestock. Higher average genetic distance existed between European and Asian populations than between European or between Asian populations. Asian breeds harboured much larger variability and higher average heterozygosity than European breeds. The samples of wild boar that have been analyzed displayed more extensive genetic variation than domestic breeds. The average linkage disequilibrium (LD) in improved pig breeds extended to 1 - 3 cM, much larger than that in human (~ 30 kb) and cattle (~ 100 kb), but smaller than that in sheep (~ 10 cM). European breeds showed greater LD that decayed more slowly than Asian breeds. We briefly discuss some processes for maintaining genomic diversity in pig, including migration, introgression, selection, and drift. We conclude that, due to the long time of domestication, the pig possesses lower heterozygosity, higher F(IS), and larger LD compared with human and cattle. This implies that a smaller effective population size and less informative markers are needed in pig for genome wide association studies.     

Genome-wide assessment of genetic diversity and population structure insights into admixture and introgression in Chinese indigenous cattle

Background

China exhibits a great diversity of ecosystems and abundant cattle resources, with nearly 30 million cattle from 53 indigenous breeds reared in specific geographic regions. To explore the genetic diversity and population structure of Chinese indigenous cattle, a population genetic analysis at both the individual and population levels was conducted and the admixture analysis was performed. We genotyped 572 samples from 20 Chinese indigenous cattle breeds using GeneSeek Genomic Profiler Bovine LD (GGP-LD, 30?K) and downloaded the published data of 77 samples from 4 worldwide commercial breeds genotyped with Illumina BovineSNP50 Beadchip (SNP50, 50?K).

Results

In principal component analysis (PCA) and neighbour-joining (NJ) tree analysis, samples of the same breeds were grouped together, leading to clear separation from other breeds. And Chinese indigenous cattle were clustered into two groups of southern and northern breeds, originated from Asian Bos indicus lineage and Eurasian Bos taurus lineage, respectively. In STRUCTURE K?=?2, a clear transition occurred from the northern breeds to the southern breeds. Additionally, the northern breeds contained a smaller Eurasian taurine (62.5%) descent proportion than that reported previously (more than 90%). In STRUCTURE K?=?3, a distinct descent was detected in the southern Chinese breeds, which could reflect a long-term selection history of Chinese indigenous cattle. The results from TreeMix and f3 statistic provided the evidence of an admixture history between southern breeds and northern breeds.

Conclusions

Consistent with the observed geographical distributions, Chinese indigenous cattle were divided into two genetic clusters, northern indigenous cattle and southern indigenous cattle. Three improved breeds in the northern area also exhibited northern indigenous ancestry. We found that the breeds distributed in the northern China showed more southern lineage introgression than previously reported. Central-located populations appeared to the admixture between southern and northern lineages, and introgression events from European cattle were observed in Luxi Cattle, Qinchuan Cattle and Jinnan Cattle. The study revealed the population structures and levels of admixture pattern among Chinese indigenous cattle, shedding light on the origin and evolutionary history of these breeds.

     

Endangered Pinzgauer cattle subtype Jochberger Hummel are genetically distinct

Autochthonous cattle breeds are an important part of cultural heritage and reservoir of genetic diversity. Usually, such breeds have been selected over centuries and reflect adaptation to a specific local environment and human demands. However, owing to low effective population size in conjunction with increased inbreeding and genetic drift, many of these lineages are threatened with extinction. The Jochberger Hummel (‘Jochberger Bumblebee’) is such an endangered subtype of the Pinzgauer cattle originating from a single polled female calf. To evaluate the suspected uniqueness of this subtype and to assess whether it should be kept separate or managed together with the Pinzgauer cattle as one population, I examined the genetic diversity in a set of 844 cattle (Angus, Charolais, Holstein, Jochberger Hummel, Pinzgauer, Uckermaerker and Tux-Zillertaler) using principal component and biogeographical ancestry analysis. The analysis showed that the Jochberger Hummel is still a distinct subtype of Pinzgauer cattle with less than 5% admixture and a low inbreeding coefficient.     

Inference of population structure of purebred dairy and beef cattle using high-density genotype data

Information on the genetic diversity and population structure of cattle breeds is useful when deciding the most optimal, for example, crossbreeding strategies to improve phenotypic performance by exploiting heterosis. The present study investigated the genetic diversity and population structure of the most prominent dairy and beef breeds used in Ireland. Illumina high-density genotypes (777 962 single nucleotide polymorphisms; SNPs) were available on 4623 purebred bulls from nine breeds; Angus (n=430), Belgian Blue (n=298), Charolais (n=893), Hereford (n=327), Holstein-Friesian (n=1261), Jersey (n=75), Limousin (n=943), Montbéliarde (n=33) and Simmental (n=363). Principal component analysis revealed that Angus, Hereford, and Jersey formed non-overlapping clusters, representing distinct populations. In contrast, overlapping clusters suggested geographical proximity of origin and genetic similarity between Limousin, Simmental and Montbéliarde and to a lesser extent between Holstein, Friesian and Belgian Blue. The observed SNP heterozygosity averaged across all loci was 0.379. The Belgian Blue had the greatest mean observed heterozygosity (H_O=0.389) among individuals within breed while the Holstein-Friesian and Jersey populations had the lowest mean heterozygosity (H_O=0.370 and 0.376, respectively). The correlation between the genomic-based and pedigree-based inbreeding coefficients was weak (r=0.171; P<0.001). Mean genomic inbreeding estimates were greatest for Jersey (0.173) and least for Hereford (0.051). The pair-wise breed fixation index (F_st) ranged from 0.049 (Limousin and Charolais) to 0.165 (Hereford and Jersey). In conclusion, substantial genetic variation exists among breeds commercially used in Ireland. Thus custom-mating strategies would be successful in maximising the exploitation of heterosis in crossbreeding strategies.     

Insights into the Genetic History of French Cattle from Dense SNP Data on 47 Worldwide Breeds

Background

Modern cattle originate from populations of the wild extinct aurochs through a few domestication events which occurred about 8,000 years ago. Newly domesticated populations subsequently spread worldwide following breeder migration routes. The resulting complex historical origins associated with both natural and artificial selection have led to the differentiation of numerous different cattle breeds displaying a broad phenotypic variety over a short period of time.

Methodology/Principal Findings

This study gives a detailed assessment of cattle genetic diversity based on 1,121 individuals sampled in 47 populations from different parts of the world (with a special focus on French cattle) genotyped for 44,706 autosomal SNPs. The analyzed data set consisted of new genotypes for 296 individuals representing 14 French cattle breeds which were combined to those available from three previously published studies. After characterizing SNP polymorphism in the different populations, we performed a detailed analysis of genetic structure at both the individual and population levels. We further searched for spatial patterns of genetic diversity among 23 European populations, most of them being of French origin, under the recently developed spatial Principal Component analysis framework.

Conclusions/Significance

Overall, such high throughput genotyping data confirmed a clear partitioning of the cattle genetic diversity into distinct breeds. In addition, patterns of differentiation among the three main groups of populations—the African taurine, the European taurine and zebus—may provide some additional support for three distinct domestication centres. Finally, among the European cattle breeds investigated, spatial patterns of genetic diversity were found in good agreement with the two main migration routes towards France, initially postulated based on archeological evidence.     

Genetic relationships between Japanese native and commercial breeds using 70 chicken autosomal SNP genotypes by the DigiTag2 assay

Recently, single nucleotide polymorphisms (SNPs) have been used to identify genes or genomic regions responsible for economic traits, including genetic diseases in domestic animals, and to examine genetic diversity of populations. In this study, we genotyped 70 chicken autosomal SNPs using DigiTag2 assay to understand the genetic structure of the Japanese native chicken breeds Satsumadori and Ingie, and the relationship of these breeds with other established breeds, Rhode Island Red (RIR), commercial broiler and layer. Five breeds, each consisting of approximately 20 chickens, were subjected to the assay, revealing the following: Average expected heterozygosities of broiler, Satsumadori, RIR, layer and Ingie were 0.265, 0.254, 0.244, 0.179 and 0.176, respectively. Phylogenetic analysis using the concatenated 70 autosomal SNP genotypes distinguished all chickens and formed clusters of chickens belonging to the respective breeds. In addition, the 2-D scatter plot of the first two principal components was consistent with the phylogenic tree. Taken together with the pairwise F(st) distances, broiler and RIR were closely positioned near each other, while Ingie was positioned far from the other breeds. Structure analysis revealed that the probable number of genetic clusters (K) was six and four with maximum likelihood and ΔK values, respectively. The clustering with maximum likelihood revealed that, in addition to the clustering of the other five breeds, the Satsumadori was subdivided into two genetic clusters. The clustering with ΔK value indicated that the broiler and Rhode Island Red were assigned to the same genetic cluster.     

Characterization of the Genetic Diversity and Population Structure for the Yellow Cattle in Taiwan Based on Microsatellite Markers

In recent years, the population size of Taiwan yellow cattle has drastically declined, even become endangered. A preservation project, Taiwan Yellow Cattle Genetic Preservation Project (TYCGPP), was carried out at the Livestock Research Institute (LRI) Hengchun branch (1988–present). An analysis of intra- and inter- population variability was performed to be the first step to preserve this precious genetic resource. In this work, a total number of 140 individuals selected from the five Taiwan yellow cattle populations were analyzed using 12 microsatellite markers (loci). These markers determined the level of genetic variation within and among populations as well as the phylogenetic structure. The total number of alleles detected (122, 10.28 per locus) and the expected heterozygosity (0.712) indicated that these five populations had a high level of genetic variability. Bayesian cluster analysis showed that the most likely number of groups was 2 (K = 2). Genetic differentiation among clusters was moderate (F _ST = 0.095). The result of AMOVA showed that yellow cattle in Taiwan had maintained a high level of within-population genetic differentiation (91%), the remainder being accounted for by differentiation among subpopulations (4%), and by differentiation among regions (5%). The results of STRUCTURE and principal component analysis (PCA) revealed two divergent clusters. The individual unrooted phylogenetic tree showed that some Kinmen yellow cattle in the Hengchun facility (KMHC individuals) were overlapped with Taiwan yellow cattle (TW) and Taiwan yellow cattle Hengchun (HC) populations. Also, they were overlapped with Kinmen × Taiwan (KT) and Kinmen yellow cattle (KM) populations. It is possible that KMHC kept similar phenotypic characteristics and analogous genotypes between TW and KM. A significant inbreeding coefficient (F _IS = 0.185; P < 0.01) was detected, suggesting a medium level of inbreeding for yellow cattle in Taiwan. The hypothesis that yellow cattle in Taiwan were derived from two different clusters was also supported by the phylogenetic tree constructed by the UPGMA, indicating that the yellow cattle in Taiwan and in Kinmen should be treated as two different management units. This result will be applied to maintain a good level of genetic variability and rusticity (stress-resistance) and to avoid further inbreeding for yellow cattle population in Taiwan.     

Genetic characterization of Argentine and Bolivian Creole cattle breeds assessed through microsatellites

In the present report, the polymorphisms from 9 microsatellites were used to assess genetic diversity and relationships in 4 Creole cattle breeds from Argentina and Bolivia, 4 European taurine breeds, and 2 American zebu populations. The Creole populations display a relatively high level of genetic variation as estimated by allelic diversity and heterozygosity, whereas the British breeds displayed reduced levels of genetic diversity. The analysis of molecular variance indicated that 7.8% of variance can be explained by differences among taurine and zebu breeds. Consistent with these results, the first principal component (PC), which comprised the 40% of the total variance, clearly distinguishes these 2 groups. In addition, all constructed phylogenetic trees cluster together Nelore and Brahman breeds with robust bootstrap values. Only 1% of variance was due to difference between American Creole and European taurine cattle. Although this secondary split was supported by the classical genetic distance and the second PC (15%), the topology of trees is not particularly robust. The presence of zebu-specific alleles in Creole cattle allowed estimating a moderate degree of zebu admixture. When these data were compared with mitochondrial and Y chromosomal studies, a clear pattern of male-mediated introgression was revealed. The results presented here contribute to the understanding of origin and history of the American Creole cattle.     

Massively parallel sequencing of Chikso (Korean brindle cattle) to discover genome-wide SNPs and InDels

Since the completion of the bovine sequencing projects, a substantial number of genetic variations such as single nucleotide polymorphisms have become available across the cattle genome. Recently, cataloguing such genetic variations has been accelerated using massively parallel sequencing technology. However, most of the recent studies have been concentrated on European Bos taurus cattle breeds, resulting in a severe lack of knowledge for valuable native cattle genetic resources worldwide. Here, we present the first whole-genome sequencing results for an endangered Korean native cattle breed, Chikso, using the Illumina HiSeq 2,000 sequencing platform. The genome of a Chikso bull was sequenced to approximately 25.3-fold coverage with 98.8% of the bovine reference genome sequence (UMD 3.1) covered. In total, 5,874,026 single nucleotide polymorphisms and 551,363 insertion/deletions were identified across all 29 autosomes and the X-chromosome, of which 45% and 75% were previously unknown, respectively. Most of the variations (92.7% of single nucleotide polymorphisms and 92.9% of insertion/deletions) were located in intergenic and intron regions. A total of 16,273 single nucleotide polymorphisms causing missense mutations were detected in 7,111 genes throughout the genome, which could potentially contribute to variation in economically important traits in Chikso. This study provides a valuable resource for further investigations of the genetic mechanisms underlying traits of interest in cattle, and for the development of improved genomics-based breeding tools.     

Characterization of bovine MHC DRB3 diversity in Latin American Creole cattle breeds

In cattle, bovine leukocyte antigens (BoLAs) have been extensively used as markers for diseases and immunological traits. However, none of the highly adapted Latin American Creole breeds have been characterized for BoLA gene polymorphism by high resolution typing methods. In this work, we sequenced exon 2 of the BoLA class II DRB3 gene from 179 cattle (113 Bolivian Yacumeño cattle and 66 Colombian Hartón del Valle cattle breeds) using a polymerase chain reaction sequence-based typing (PCR-SBT) method. We identified 36 previously reported alleles and three novel alleles. Thirty-five (32 reported and three new) and 24 alleles (22 reported and two new) were detected in Yacumeño and Hartón del Valle breeds, respectively. Interestingly, Latin American Creole cattle showed a high degree of gene diversity despite their small population sizes, and 10 alleles including three new alleles were found only in these two Creole breeds. We next compared the degree of genetic variability at the population and sequence levels and the genetic distance in the two breeds with those previously reported in five other breeds: Holstein, Japanese Shorthorn, Japanese Black, Jersey, and Hanwoo. Both Creole breeds presented gene diversity higher than 0.90, a nucleotide diversity higher than 0.07, and mean number of pairwise differences higher than 19, indicating that Creole cattle had similar genetic diversity at BoLA-DRB3 to the other breeds. A neutrality test showed that the high degree of genetic variability may be maintained by balancing selection. The F_ST index and the exact G test showed significant differences across all cattle populations (F_ST = 0.0478; p < 0.001). Results from the principal components analysis and the phylogenetic tree showed that Yacumeño and Hartón del Valle breeds were closely related to each other. Collectively, our results suggest that the high level of genetic diversity could be explained by the multiple origins of the Creole germplasm (European, African and Indicus), and this diversity might be maintained by balancing selection.     

Genome-wide genetic diversity,population structure and admixture analysis in African and Asian cattle breeds

Knowledge about genetic diversity and population structure is useful for designing effective strategies to improve the production, management and conservation of farm animal genetic resources. Here, we present a comprehensive genome-wide analysis of genetic diversity, population structure and admixture based on 244 animals sampled from 10 cattle populations in Asia and Africa and genotyped for 69 903 autosomal single-nucleotide polymorphisms (SNPs) mainly derived from the indicine breed. Principal component analysis, STRUCTURE and distance analysis from high-density SNP data clearly revealed that the largest genetic difference occurred between the two domestic lineages (taurine and indicine), whereas Ethiopian cattle populations represent a mosaic of the humped zebu and taurine. Estimation of the genetic influence of zebu and taurine revealed that Ethiopian cattle were characterized by considerable levels of introgression from South Asian zebu, whereas Bangladeshi populations shared very low taurine ancestry. The relationships among Ethiopian cattle populations reflect their history of origin and admixture rather than phenotype-based distinctions. The high within-individual genetic variability observed in Ethiopian cattle represents an untapped opportunity for adaptation to changing environments and for implementation of within-breed genetic improvement schemes. Our results provide a basis for future applications of genome-wide SNP data to exploit the unique genetic makeup of indigenous cattle breeds and to facilitate their improvement and conservation.     

Whole-Genome Resequencing Analysis of Hanwoo and Yanbian Cattle to Identify Genome-Wide SNPs and Signatures of Selection

Over the last 30 years, Hanwoo has been selectively bred to improve economically important traits. Hanwoo is currently the representative Korean native beef cattle breed, and it is believed that it shared an ancestor with a Chinese breed, Yanbian cattle, until the last century. However, these two breeds have experienced different selection pressures during recent decades. Here, we whole-genome sequenced 10 animals each of Hanwoo and Yanbian cattle (20 total) using the Illumina HiSeq 2000 sequencer. A total of approximately 3.12 and 3.07 billion sequence reads were mapped to the bovine reference sequence assembly (UMD 3.1) at an average of approximately 10.71- and 10.53-fold coverage for Hanwoo and Yanbian cattle, respectively. A total of 17,936,399 single nucleotide polymorphisms (SNPs) were yielded, of which 22.3% were found to be novel. By annotating the SNPs, we further retrieved numerous nonsynonymous SNPs that may be associated with traits of interest in cattle. Furthermore, we performed whole-genome screening to detect signatures of selection throughout the genome. We located several promising selective sweeps that are potentially responsible for economically important traits in cattle; the PPP1R12A gene is an example of a gene that potentially affects intramuscular fat content. These discoveries provide valuable genomic information regarding potential genomic markers that could predict traits of interest for breeding programs of these cattle breeds.     

Genetic characterization of southwestern European bovine breeds: a historical and biogeographical reassessment with a set of 16 microsatellites

The origin of Iberian cattle has been suggested by some authors to be the product of European and north African cattle entrances during the last few thousands of years. However, these hypotheses were mainly based on morphological similarities. This study analyzed 889 unrelated individuals from 15 representative Iberian breeds and 3 French breeds for 16 microsatellite loci. Statistical tests were used to calculate interpopulation genetic distances (D(A)) and principal components analysis (PCA). To visualize the geographical distribution of the genetic differentiation between Iberian cattle breeds, data from the PCA analysis were used to construct synthetic maps. Genetic similarity among neighboring Iberian breeds is mainly caused by gene flow. However, recent demographic fluctuations and reproductive isolation in Alistana, Mirandesa, and Tudanca has increased genetic drift, which may be the main cause for the relatively high differentiation of these populations. The synthetic maps constructed with the first and second PCs revealed (1) a large differentiation between Northern Iberian breeds rather than between more geographically distant breeds, and (2) a clear east-west gradient that may be related with the model of demic diffusion of agriculture. Finally, we detected no strong evidence for an African genetic influence in the Iberian cattle breeds analyzed in this study.