A phylogenic study of Pinzgau cattle by the methods of immunogenetic analysis]

Phylogenesis of Pinzgau cattle was studied by the method of cluster analysis. Dendrogram of phylogeny of breeds, which participated in creation of Pinzgau breed, was built. The low level of genetic differences between modern populations of Pinzgau breed was determined.     

Genetic diversity among Canadienne, Brown Swiss, Holstein, and Jersey cattle of Canada based on 15 bovine microsatellite markers.

The genetic diversity among Canadienne, Brown Swiss, Holstein, and Jersey cattle was estimated from relationships determined by genotyping 20 distantly related animals in each breed for 15 microsatellites located on separate chromosomes. The Canadienne, Holstein, and Jersey cattle had an average of six alleles per loci compared with five alleles for Brown Swiss. Furthermore, a number of potentially breed-specific alleles were identified. The allele size variance among breeds was similar, but varied considerably among loci. All of the loci studied were equally heterozygous, as were Brown Swiss, Canadienne, and Holstein cattle (0.68-0.69) whereas Jersey cattle showed lower heterozygosity (0.59). The within-breed estimates of genetic distance were greater than zero and significant. The genetic distance between Canadienne and Holstein (0.156), Brown Swiss (0.243), and Jersey (0.235) was negligible, suggesting close relationship. Concurrently, Brown Swiss and Holstein (0.211) cattle also demonstrated close relationship. In contrast, the Jersey breed was genetically distant from the Brown Swiss and Holstein cattle (0.427 and 0.320, respectively). The characterization of Canadienne cattle, as part of the genetic resource conservation effort currently underway in Canada, underscores the difficulty in scientifically establishing unique breeds. Therefore, the need to consider all relevant morphological characteristics and production performance in combination with available cultural, historical, pedigree, and molecular information becomes relevant when identifying breeds for conservation.     

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.     

SNP genotypes reveal breed substructure,selection signatures and highly inbred regions in Piétrain pigs

The Piétrain pig originates from the Belgian village Piétrain some time between 1920 and 1950. Owing to its superior conformation, the Piétrain has spread worldwide since the 1960s. As initial population sizes were limited and close inbreeding was commonplace, the breed’s genetic diversity has been questioned. Therefore, this study examines Piétrain breed substructure, diversity and selection signatures using SNP data in comparison with Duroc, Landrace and Large White populations. Principal component analysis indicated three subpopulations, and F_ST analysis showed that US Piétrains differ most from European Piétrains. Average inbreeding based on runs of homozygosity (ROH) segments larger than 4 Mb ranged between 16.7 and 20.9%. The highest chromosomal inbreeding levels were found on SSC8 (42.7%). ROH islands were found on SSC8, SSC15 and SSC18 in all Piétrain populations, but numerous population-specific ROH islands were also detected. Moreover, a large ROH island on SSC8 (34–126 Mb) appears nearly fixed in all Piétrain populations, with a unique genotype. Chromosomal ROH patterns were similar between Piétrain populations. This study shows that Piétrain populations are genetically diverging, with at least three genetically distinct populations worldwide. Increasing genetic diversity in local Piétrain populations by introgression from other Piétrain populations seems to be only limited. Moreover, a unique 90 Mb region on SSC8 appeared largely fixed in the Piétrain breed, indicating that fixation was already present before the 1960s. We believe that strong selection and inbreeding during breed formation fixed these genomic regions in Piétrains. Finally, we hypothesize that independent coat color selection may have led to large ROH pattern similarities on SSC8 between unrelated pig breeds.     

Genetic diversity analyses reveal first insights into breed‐specific selection signatures within Swiss goat breeds

We used genotype data from the caprine 50k Illumina BeadChip for the assessment of genetic diversity within and between 10 local Swiss goat breeds. Three different cluster methods allowed the goat samples to be assigned to the respective breed groups, whilst the samples of Nera Verzasca and Tessin Grey goats could not be differentiated from each other. The results of the different genetic diversity measures show that Appenzell, Toggenburg, Valais and Booted goats should be prioritized in future conservation activities. Furthermore, we examined runs of homozygosity (ROH) and compared genomic inbreeding coefficients based on ROH (F_ROH) with pedigree‐based inbreeding coefficients (F_PED). The linear relationship between F_ROH and F_PED was confirmed for goats by including samples from the three main breeds (Saanen, Chamois and Toggenburg goats). F_ROH appears to be a suitable measure for describing levels of inbreeding in goat breeds with missing pedigree information. Finally, we derived selection signatures between the breeds. We report a total of 384 putative selection signals. The 25 most significant windows contained genes known for traits such as: coat color variation (MITF, KIT, ASIP), growth (IGF2, IGF2R, HRAS, FGFR3) and milk composition (PITX2). Several other putative genes involved in the formation of populations, which might have been selected for adaptation to the alpine environment, are highlighted. The results provide a contemporary background for the management of genetic diversity in local Swiss goat breeds.     

Analysis of runs of homozygosity and their relationship with inbreeding in five cattle breeds farmed in Italy

Increased inbreeding is an inevitable consequence of selection in livestock populations. The analysis of high‐density single nucleotide polymorphisms (SNPs) facilitates the identification of long and uninterrupted runs of homozygosity (ROH) that can be used to identify chromosomal regions that are identical by descent. In this work, the distribution of ROH of different lengths in five Italian cattle breeds is described. A total of 4095 bulls from five cattle breeds (2093 Italian Holstein, 749 Italian Brown, 364 Piedmontese, 410 Marchigiana and 479 Italian Simmental) were genotyped at 54K SNP loci. ROH were identified and used to estimate molecular inbreeding coefficients (F_ROH), which were compared with inbreeding coefficients estimated from pedigree information (F_PED) and using the genomic relationship matrix (F_GRM). The average number of ROH per animal ranged from 54 ± 7.2 in Piedmontese to 94.6 ± 11.6 in Italian Brown. The highest number of short ROH (related to ancient consanguinity) was found in Piedmontese, followed by Simmental. The Italian Brown and Holstein had a higher proportion of longer ROH distributed across the whole genome, revealing recent inbreeding. The F_PED were moderately correlated with F_ROH > 1 Mb (0.662, 0.700 and 0.669 in Italian Brown, Italian Holstein and Italian Simmental respectively) but poorly correlated with F_GRM (0.134, 0.128 and 0.448 for Italian Brown, Italian Holstein and Italian Simmental respectively). The inclusion of ROH > 8 Mb in the inbreeding calculation improved the correlation of F_ROH with F_PED and F_GRM. ROH are a direct measure of autozygosity at the DNA level and can overcome approximations and errors resulting from incomplete pedigree data. In populations with high linkage disequilibrium (LD) and recent inbreeding (e.g. Italian Holstein and Italian Brown), a medium‐density marker panel, such as the one used here, may provide a good estimate of inbreeding. However, in populations with low LD and ancient inbreeding, marker density would have to be increased to identify short ROH that are identical by descent more precisely.     

Indigenous domestic breeds as reservoirs of genetic diversity: the Argentinean Creole cattle

Contrary to highly selected commercial breeds, indigenous domestic breeds are composed of semi-wild or feral populations subjected to reduced levels of artificial selection. As a consequence, many of these breeds have become locally adapted to a wide range of environments, showing high levels of phenotypic variability and increased fitness under natural conditions. Genetic analyses of three loci associated with milk production (alpha(S1)-casein, kappa-casein and prolactin) and the locus BoLA-DRB3 of the major histocompatibility complex indicated that the Argentinean Creole cattle (ACC), an indigenous breed from South America, maintains high levels of genetic diversity and population structure. In contrast to the commercial Holstein breed, the ACC showed considerable variation in heterozygosity (H(e)) and allelic diversity (A) across populations. As expected, bi-allelic markers showed extensive variation in He whereas the highly polymorphic BoLA-DRB3 showed substantial variation in A, with individual populations having 39-74% of the total number of alleles characterized for the breed. An analysis of molecular variance (AMOVA) of nine populations throughout the distribution range of the ACC revealed that 91.9-94.7% of the total observed variance was explained by differences within populations whereas 5.3-8.1% was the result of differences among populations. In addition, the ACC breed consistently showed higher levels of genetic differentiation among populations than Holstein. Results from this study emphasize the importance of population genetic structure within domestic breeds as an essential component of genetic diversity and suggest that indigenous breeds may be considered important reservoirs of genetic diversity for commercial domestic species.     

Runs of homozygosity and population history in cattle

ABSTRACT: BACKGROUND: Runs of homozygosity (ROH) are contiguous lengths of homozygous genotypes that are present in an individual due to parents transmitting identical haplotypes to their offspring. The extent and frequency of ROHs may inform on the ancestry of an individual and its population. Here we use high density (n = 777,962) bi-allelic SNPs in a range of cattle breed samples to correlate ROH with the pedigree-based inbreeding coefficients and to validate subsequent analyses using 54,001 SNP genotypes. This study provides a first testing of the inference drawn from ROH through comparison with estimates of inbreeding from calculations based on the detailed pedigree data available for several breeds. RESULTS: All animals genotyped on the HD panel displayed at least one ROH that was between 1--5 Mb in length with certain regions of the genome more likely to be involved in a ROH than others. Strong correlations (r = 0.75, p < 0.0001) existed between the pedigree-based inbreeding coefficient and a statistic based on sum of ROH of length > 0.5 KB and suggests that in the absence of an animal's pedigree data, the extent of a genome under ROH may be used to infer aspects of recent population history even from relatively few samples. CONCLUSIONS: Our findings suggest that ROH are frequent across all breeds but differing patterns of ROH length and burden illustrate variations in breed origins and recent management.     

Genomic inbreeding estimation in small populations: evaluation of runs of homozygosity in three local dairy cattle breeds

In the local breeds with small population size, one of the most important problems is the increase of inbreeding coefficient (F). High levels of inbreeding lead to reduced genetic diversity and inbreeding depression. The availability of high-density single nucleotide polymorphism (SNP) arrays has facilitated the quantification of F by genomic markers in farm animals. Runs of homozygosity (ROH) are contiguous lengths of homozygous genotypes and represent an estimate of the degree of autozygosity at genome-wide level. The current study aims to quantify the genomic F derived from ROH (F_ROH) in three local dairy cattle breeds. F_ROH values were compared with F estimated from the genomic relationship matrix (F_GRM), based on the difference between observed v. expected number of homozygous genotypes (F_HOM) and the genomic homozygosity of individual i (F_{MOL i}). The molecular coancestry coefficient (f_{MOL ij}) between individuals i and j was also estimated. Individuals of Cinisara (71), Modicana (72) and Reggiana (168) were genotyped with the 50K v2 Illumina BeadChip. Genotypes from 96 animals of Italian Holstein cattle breed were also included in the analysis. We used a definition of ROH as tracts of homozygous genotypes that were >4 Mb. Among breeds, 3661 ROH were identified. Modicana showed the highest mean number of ROH per individual and the highest value of F_ROH, whereas Reggiana showed the lowest ones. Differences among breeds existed for the ROH lengths. The individuals of Italian Holstein showed high number of short ROH segments, related to ancient consanguinity. Similar results showed the Reggiana with some extreme animals with segments covering 400 Mb and more of genome. Modicana and Cinisara showed similar results between them with the total length of ROH characterized by the presence of large segments. High correlation was found between F_HOM and F_ROH ranged from 0.83 in Reggiana to 0.95 in Cinisara and Modicana. The correlations among F_ROH and other estimated F coefficients were generally lower ranged from 0.45 (F_{MOL i}−F_ROH) in Cinisara to 0.17 (F_GRM−F_ROH) in Modicana. On the basis of our results, recent inbreeding was observed in local breeds, considering that 16 Mb segments are expected to present inbreeding up to three generations ago. Our results showed the necessity of implementing conservation programs to control the rise of inbreeding and coancestry in the three Italian local dairy cattle breeds.     

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.     

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.     

Genetic relationships of Austrian cattle breeds

Methods of phylogenetic analysis were used to study the relationships of ten Central European cattle breeds using gene frequencies at eleven blood group and protein polymorphism loci. The results show remarkable agreement with historical and geographical relationships, but are different from the relationships postulated earlier in the century from studies of skull shape. Lowland cattle of Northern Germany are distinct from the West Alpine breeds, but the Pustertaler and Pinzgauer breeds from the Eastern Alps are more closely related to those lowland breeds than to the other Alpine breeds, suggesting they were brought from the North by Germanic settlers in the early Middle Ages. Thus, a sharp distinction between lowland cattle and all Alpine breeds seems unwarranted. Moreover, we find no evidence to justify a unique taxonomic position for Pustertaler, nor do we find any evidence that Spanish cattle might have influenced the Pinzgau. We find a close genetic relationship of the Murbodner breed to the Fleckvieh (Simmenthal) breeds that is not supported by any known historical relationship; it may indicate a little-known human migration.     

Microsatellite Analysis of Genetic Population Structure of Zebu Cattle (Bos indicus) Breeds from North-Western Region of India

The present study aims to understand the existing genetic diversity and structure of six native cattle breeds (Rathi, Tharparkar, Nagori, Mewati, Gir, and Kankrej) adapted to the north-western arid and semi-arid region of India based on microsatellite loci. Various diversity estimates, mean number of alleles (12.84); effective number of alleles (5.02); gene diversity (0.769), and observed heterozygosity (0.667) reflected the existence of substantial within-breed diversity in all the investigated cattle breeds. Mean estimates of F-statistics: F_IT = 0.144 ± 0.023, F_IS = 0.071 ± 0.021, and F_ST = 0.078 ± 0.014 were significantly different from zero (P < 0.05). The interbreed relationships indicated moderate level of breed differentiation between the six cattle breeds with least differentiation between Kankrej-Mewati pair. The phylogeny structuring further supported close grouping of Kankrej and Mewati breeds. Correspondence analysis plotted Rathi, Tharparkar, and Gir individuals into three separate areas of multivariate space; whereas, Kankrej, Mewati, and Nagori cattle showed low breed specific clustering. This reflected the existence of discrete genetic structure for Tharparkar, Rathi, and Gir, the prominent dairy breeds of the region; whereas, admixture was observed for Kankrej, Mewati, and Nagori individuals.     

Genome-wide analysis in endangered populations: a case study in Barbaresca sheep

Analysis of genomic data is becoming increasingly common in the livestock industry and the findings have been an invaluable resource for effective management of breeding programs in small and endangered populations. In this paper, with the goal of highlighting the potential of genomic analysis for small and endangered populations, genome-wide levels of linkage disequilibrium, measured as the squared correlation coefficient of allele frequencies at a pair of loci, effective population size, runs of homozygosity (ROH) and genetic diversity parameters, were estimated in Barbaresca sheep using Illumina OvineSNP50K array data. Moreover, the breed’s genetic structure and its relationship with other breeds were investigated. Levels of pairwise linkage disequilibrium decreased with increasing distance between single nucleotide polymorphisms. An average correlation coefficient <0.25 was found for markers located up to 50 kb apart. Therefore, these results support the need to use denser single nucleotide polymorphism panels for high power association mapping and genomic selection efficiency in future breeding programs. The estimate of past effective population size ranged from 747 animals 250 generations ago to 28 animals five generations ago, whereas the contemporary effective population size was 25 animals. A total of 637 ROH were identified, most of which were short (67%) and ranged from 1 to 10 Mb. The genetic analyses revealed that the Barbaresca breed tended to display lower variability than other Sicilian breeds. Recent inbreeding was evident, according to the ROH analysis. All the investigated parameters showed a comparatively narrow genetic base and indicated an endangered status for Barbaresca. Multidimensional scaling, model-based clustering, measurement of population differentiation, neighbor networks and haplotype sharing distinguished Barbaresca from other breeds, showed a low level of admixture with the other breeds considered in this study, and indicated clear genetic differences compared with other breeds. Attention should be given to the conservation of Barbaresca due to its critical conservation status. In this context, genomic information may have a crucial role in management of small and endangered populations.     

Unraveling the genetic diversity of Belgian Milk Sheep using medium-density SNP genotypes

The present study focuses on the Belgian Milk Sheep in Flanders (Belgium) and compares its genetic diversity and relationship with the Flemish Sheep, the Friesian Milk Sheep, the French Lacaune dairy sheep and other Northern European breeds. For this study, 94 Belgian Milk Sheep, 23 Flemish Sheep and 22 Friesian Milk Sheep were genotyped with the OvineSNP50 array. In addition, 29 unregistered animals phenotypically similar to Belgian Milk Sheep were genotyped using the 15K ISGC chip. Both Belgian and Friesian Milk Sheep as well as the East Friesian Sheep were found to be less diverse than the other seven breeds included in this study. Genomic inbreeding coefficients based on runs of homozygosity (ROH) were estimated at 14.5, 12.4 and 10.2% for Belgian Milk Sheep, Flemish Sheep and Friesian Milk Sheep respectively. Out of 29 unregistered Belgian Milk Sheep, 28 mapped in the registered Belgian Milk Sheep population. Ancestry analysis, PCA and F_ST calculations showed that Belgian Milk Sheep are more related to Friesian Milk Sheep than to Flemish Sheep, which was contrary to the breeders' expectations. Consequently, breeders may prefer to crossbreed Belgian Milk Sheep with Friesian sheep populations (Friesian Milk Sheep or East Friesian Sheep) in order to increase diversity. This research underlines the usefulness of SNP chip genotyping and ROH analyses for monitoring genetic diversity and studying genetic links in small livestock populations, profiting from internationally available genotypes. As assessment of genetic diversity is vital for long-term breed survival, these results will aid flockbooks to preserve genetic diversity.     

Use of breed-specific single nucleotide polymorphisms to discriminate between Holstein and Jersey dairy cattle breeds

Emphasis on livestock genetic improvement in the past decades has led to commercialization of different breeds of livestock species. Breed validation has become increasingly important to assess the safety and authenticity of livestock products in global and domestic markets. The objective of this study was to evaluate the use of breed-specific single nucleotide polymorphisms (SNPs) in discriminating between Holstein and Jersey dairy cattle breeds. Two separate resource populations were used, including a reference population consisting of 498 Holstein and 83 Jersey bull DNA samples, and a validation population consisting of 260 Holstein and 34 Jersey cow DNA samples. Five Jersey-specific and four Holstein-specific SNPs were identified and genotyped on the reference and validation resource populations. The reference population was used to validate the breed-specific SNPs used in this study and to predict the allocation efficiencies and misclassification probabilities of different combinations of SNPs. Individual animals in the validation population were allocated to either breed based on the presence of breed-specific alleles. It was found that any combination of three breed-specific SNPs had, on average, high breed allocation efficiency of >95% and low misclassification probability of <5%. In conclusion, this study demonstrates a simple, yet effective, method of using breed-specific SNPs to discriminate between Jersey and Holstein cattle breeds.     

Marker-assisted conservation of European cattle breeds: An evaluation

Two methods have been developed for the assessment of conservation priorities on the basis of molecular markers. According to the Weitzman approach, contributions to genetic diversity are derived from genetic distances between populations. Alternatively, diversity within and across populations is optimized by minimizing marker-estimated kinships. We have applied, for the first time, both methods to a comprehensive data set of 69 European cattle breeds, including all cosmopolitan breeds and several local breeds, for which genotypes of 30 microsatellite markers in 25–50 animals per breed have been obtained. Both methods were used to calculate the gain in diversity if a breed was added to a set of nine non-endangered breeds. Weitzman-derived diversities were confounded by genetic drift in isolated populations, which dominates the genetic distances but does not necessarily increase the conservation value of a breed. Marker-estimated kinships across populations were less disturbed by genetic drift than the Weitzman diversities and assigned high conservation values to Mediterranean breeds, which indeed have genetic histories that differ from the non-endangered breeds. Prospects and limitations of marker-assisted decisions on conservation priorities are discussed.     

Genome-wide identification of runs of homozygosity islands and associated genes in local dairy cattle breeds

Runs of homozygosity (ROH) are widely used as predictors of whole-genome inbreeding levels in cattle. They identify regions that have an unfavorable effect on a phenotype when homozygous, but also identify the genes associated with traits of economic interest present in these regions. Here, the distribution of ROH islands and enriched genes within these regions in four dairy cattle breeds were investigated. Cinisara (71), Modicana (72), Reggiana (168) and Italian Holstein (96) individuals were genotyped using the 50K v2 Illumina BeadChip. The genomic regions most commonly associated with ROHs were identified by selecting the top 1% of the single nucleotide polymorphisms (SNPs) most commonly observed in the ROH of each breed. In total, 11 genomic regions were identified in Cinisara and Italian Holstein, and eight in Modicana and Reggiana, indicating an increased ROH frequency level. Generally, ROH islands differed between breeds. The most homozygous region (>45% of individuals with ROH) was found in Modicana on chromosome 6 within a quantitative trail locus affecting milk fat and protein concentrations. We identified between 126 and 347 genes within ROH islands, which are involved in multiple signaling and signal transduction pathways in a wide variety of biological processes. The gene ontology enrichment provided information on possible molecular functions, biological processes and cellular components under selection related to milk production, reproduction, immune response and resistance/susceptibility to infection and diseases. Thus, scanning the genome for ROH could be an alternative strategy to detect genomic regions and genes related to important economic traits.     

Genetic structure of Eurasian cattle (Bos taurus) based on microsatellites: clarification for their breed classification1

We pool three previously published data sets and present population genetic analyses of microsatellite variation in 48 Bos taurus cattle breeds from a wide range of geographical origins in Eurasia, mostly its northern territory. Bayesian model‐based clustering reveals six distinct clusters: besides a single‐population cluster of the Yakutian Cattle from Far Eastern Siberia and a cluster of breeds characteristic of an early origin, the other four major clusters largely correspond to previously defined morphological subgroups of Red Lowland, Lowland Black‐Pied, Longhorned Dairy and North European Polled cattle breeds. The results highlighted past expansion events of the productive breeds such as Danish Red, Angeln, Holstein‐Friesian and Ayrshire in northern and Eastern Europe. Based on genetic assignment of the breeds and the availability of breed information, we provide a preliminary classification of the five breeds that were to date undefined. Furthermore, in the analysis of molecular variance, despite some correspondence between geographical proximity and genetic similarity, the breed classification appears to be a better predictor of genetic structure in the cattle populations (the among‐group variance component: breed classification, 2.47%, P < 0.001; geographical division, 0.77%, P < 0.001).     

The footprint of recent and strong demographic decline in the genomes of Mangalitza pigs

The Mangalitza pig breed has suffered strong population reductions due to competition with more productive cosmopolitan breeds. In the current work, we aimed to investigate the effects of this sustained demographic recession on the genomic diversity of Mangalitza pigs. By using the Porcine Single Nucleotid Polymorphism BeadChip, we have characterized the genome-wide diversity of 350 individuals including 45 Red Mangalitza (number of samples; n=20 from Hungary and n=25 from Romania), 37 Blond Mangalitza, 26 Swallow-belly Mangalitza, 48 Blond Mangalitza × Duroc crossbreds, 5 Bazna swine, 143 pigs from the Hampshire, Duroc, Landrace, Large White and Pietrain breeds and 46 wild boars from Romania (n=18) and Hungary (n=28). Performance of a multidimensional scaling plot showed that Landrace, Large White and Pietrain pigs clustered independently from Mangalitza pigs and Romanian and Hungarian wild boars. The number and total length of ROH (runs of homozygosity), as well as F_ROH coefficients (proportion of the autosomal genome covered ROH) did not show major differences between Mangalitza pigs and other wild and domestic pig populations. However, Romanian and Hungarian Red Mangalitza pigs displayed an increased frequency of very long ROH (>30 Mb) when compared with other porcine breeds. These results indicate that Red Mangalitza pigs underwent recent and strong inbreeding probably as a consequence of severe reductions in census size.