Genomic diversity and population structure of Mexican and Spanish bovine Lidia breed

The Lidia bovine breed is distinguished for its low genetic exchangeability given its selection on aggressive behavior, its management uniqueness and its subdivided structure. In this study, we present a comprehensive genome‐wide analysis of genetic diversity, population structure and admixture of 468 animals from Mexican and Spanish Lidia breed populations and 64 samples belonging to 10 Spanish native and American‐creole breeds using 37 148 single nucleotide polymorphisms. We found similar average inbreeding values in the Lidia breed, with different distributions within groups; variability of inbreeding values among Spanish lineages was significant and no differences were found among the Mexican sub‐populations. Together, the high F_IS of the lineages and the behavior of the runs of homozygosity are consequences of the lineage's small effective population sizes, contributing to their inbreeding increase. Population admixture analysis discarded any influence on the genetic structure of the Lidia populations from the Spanish native and American‐creole breeds. In addition, both Lidia populations depicted different genetic origins, with the exception of some Mexican individuals whose origins traced back to recent Spanish importations.     

Immunogenetic and population genetic analyses of Iberian cattle

Blood samples were collected from more than 100 animals in each of 2 Spanish cattle breeds (Retinto and De Lidia), 2 Portuguese breeds (Alentejana and Mertolenga), and American Longhorn cattle. All samples for the 4 Iberian breeds were tested for 20 polymorphic systems; American Longhorn were tested for 19 of the 20. For each breed an average inbreeding coefficient was estimated by a comparison of the observed and expected heterozygosity at 7 or 8 codominant systems tested. All breeds had positive values but only 3 breeds had estimates of inbreeding that were statistically significantly different from 0: De Lidia with = 0.17, Retinto with = 0.08 and Mertolenga with f = 0.05. The De Lidia breed especially may be suffering from inbreeding depression since this high value is greater than expected if all of the animals were progeny of half-sib matings. Genetic distances were calculated from the gene frequency data on these 5 breeds plus 9 other European breeds. Analyses of these distances show a closely related group of the 4 Iberian breeds and American Longhorn, confirming the close relationships among the Iberian breeds and the Iberian, probably Portuguese, origin of American Longhorn cattle.     

Ancestral matrilineages and mitochondrial DNA diversity of the Lidia cattle breed

To clarify the genetic ancestry and the mitochondrial DNA (mtDNA) diversity of the Lidia cattle breed, a 521-bp D-loop fragment was sequenced in 527 animals belonging to 70 herds distributed across 29 lineages. The mtDNA diversity recorded was similar to that seen for Middle Eastern breeds and greater than that recorded for the majority of European breeds. Haplotype T3 was the most common (81%), followed by the African T1 haplotype (17%); very low frequencies were recorded for haplotypes T and T2. The results agree with there being two major ancestral lines for the Lidia breed, European and African, similar to that seen for other Mediterranean breeds. A wide range of variation in haplotype frequencies was seen between the examined lineages. Haplotype T3 was present in all those analysed; in five it was the only one present, and in only one lineage (Miura) was its frequency lower than that of T1. T1*, a haplotype reported in Criollo breeds and to date in only a single European breed (the Retinta breed from Spain), was found in a single animal belonging to the Concha y Sierra lineage. Network analysis of the Lidia breed revealed the presence of two major haplotypes: T3 and T1. The Lidia breed appears to be more closely related to prehistoric Iberian and Italian than to British aurochs.     

A genome‐wide scan for selection signatures in Yorkshire and Landrace pigs based on sequencing data

Pigs have experienced dramatic selection due to domestication, which has led to many different phenotypes when compared to their wild counterparts, especially in the last several decades. Currently, genome‐wide scans in both cattle and humans showing positive selection footprints have been investigated. However, few studies have focused on porcine selection footprints, particularly on a genome‐wide scale. Surveying for selection footprints across porcine genomes can be quite valuable for revealing the genetic mechanisms of phenotypic diversity. Here, we employed a medium sequencing depth (5–20x/site per individual, on average) approach called genotyping by genome reducing and sequencing (GGRS) to detect genome‐wide selection signatures of two domestic pig breeds (Yorkshire and Landrace) that have been under intensive selection for traits of muscle development, growth and behavior. The relative extended haplotype homozygosity test, which identifies selection signatures by measuring the characteristics of haplotypes’ frequency distribution within a single population, was also applied to identify potential positively selected regions. As a result, signatures of positive selection were found in each breed. However, most selection signatures were population specific and related to genomic regions containing genes for biological categories including brain development, metabolism, growth and olfaction. Furthermore, the result of the gene set enrichment analysis indicated that selected regions of the two breeds presented a different over‐representation of genes in the Gene Ontology annotations and Kyoto Encyclopedia of Genes and Genomes pathways. Our results revealed a genome‐wide map of selection footprints in pigs and may help us better understand the mechanisms of selection in pig breeding.     

Change in certain forms of aggressive behavior and the concentration of monoamines in the brain during selection of wild rats for taming

Norway rats have been selected during 20 generations by the absence of aggressive reaction to man (tamed rats). From 7 up to 20th generations of selection, different forms of aggressive behaviour (reaction to glove, intermale, shock-induced aggression and predatory aggression) were studied, and the level of noradrenaline, serotonin and its metabolite 5-hydroxyindoleacetic acid was determined in the brain. In the absence of aggressive reaction to glove in tamed rats, the shock-induced aggression considerably decreased while the predatory aggressiveness (mouse-killing behaviour) and intermale aggressiveness did not change. Beginning from 15-16th generation of selection, a higher level of the 5-hydroxyindoleacetic acid in the hypothalamus was established, in the 20th generation an increased content of serotonin was revealed in the hypothalamus and the midbrain. In some generations of selection an increased level of noradrenaline in the hypothalamus in comparison to wild rats was observed. A conclusion is made that the selection of animals by taming unequally influences different kinds of aggressiveness and is accompanied by inherited consolidated reorganization of the monoamine brain systems.     

Genetic diversity and genomic signatures of selection among cattle breeds from Siberia,eastern and northern Europe

Domestication in the near eastern region had a major impact on the gene pool of humpless taurine cattle (Bos taurus). As a result of subsequent natural and artificial selection, hundreds of different breeds have evolved, displaying a broad range of phenotypic traits. Here, 10 Eurasian B. taurus breeds from different biogeographic and production conditions, which exhibit different demographic histories and have been under artificial selection at various intensities, were investigated using the Illumina BovineSNP50 panel to understand their genetic diversity and population structure. In addition, we scanned genomes from eight breeds for signatures of diversifying selection. Our population structure analysis indicated six distinct breed groups, the most divergent being the Yakutian cattle from Siberia. Selection signals were shared (experimental P‐value < 0.01) with more than four breeds on chromosomes 6, 7, 13, 16 and 22. The strongest selection signals in the Yakutian cattle were found on chromosomes 7 and 21, where a miRNA gene and genes related to immune system processes are respectively located. In general, genomic regions indicating selection overlapped with known QTL associated with milk production (e.g. on chromosome 19), reproduction (e.g. on chromosome 24) and meat quality (e.g. on chromosome 7). The selection map created in this study shows that native cattle breeds and their genetic resources represent unique material for future breeding.     

Association of the expression levels in the skeletal muscle and a SNP in the CDC10 gene with growth‐related traits in Japanese Black beef cattle

Growth performance, as well as marbling, is the main breeding objective in Japanese Black (JB) cattle, the major beef breed in Japan. The septin 7 (CDC10) gene, involved in cellular proliferation, is located within a genomic region of a quantitative trait locus for growth‐related traits. In this study, we first showed that the expression levels of the CDC10 gene in the skeletal muscle were higher in JB steers with extremely high growth performance than in JB steers with extremely low growth, using real‐time PCR. Further, a single nucleotide polymorphism (SNP), NC_007302.5:g.63264949G>C, was detected in the promoter region of the CDC10 gene and genotyped in three Japanese cattle breeds (known as ‘Wagyu’ in Japan) and the Brown Swiss dairy cattle breed. All four cattle populations showed a moderate genetic diversity at the SNP of the CDC10 gene. An association analysis indicated that the SNP was associated with growth‐related traits in JB cattle. These findings suggest possible effects of the expression levels in the skeletal muscle and the SNP of the CDC10 gene on growth‐related traits in JB cattle. The CDC10 SNP may be useful for effective marker‐assisted selection to increase beef productivity in JB beef cattle.     

Genetic diversity and population structure of 20 North European cattle breeds

Blood samples were collected from 743 animals from 15 indigenous, 2 old imported, and 3 commercial North European cattle breeds. The samples were analyzed for 11 erythrocyte antigen systems, 8 proteins, and 10 microsatellites, and used to assess inter- and intrabreed genetic variation and genetic population structures. The microsatellites BoLA-DRBP1 and CSSM66 were nonneutral markers according to the Ewens-Watterson test, suggesting some kind of selection imposed on these loci. North European cattle breeds displayed generally similar levels of multilocus heterozygosity and allelic diversity. However, allelic diversity has been reduced in several breeds, which was explained by limited effective population sizes over the course of man-directed breed development and demographic bottlenecks of indigenous breeds. A tree showing genetic relationships between breeds was constructed from a matrix of random drift-based genetic distance estimates. The breeds were classified on the basis of the tree topology into four major breed groups, defined as Northern indigenous breeds, Southern breeds, Ayrshire and Friesian breeds, and Jersey. Grouping of Nordic breeds was supported by documented breed history and geographical divisions of native breeding regions of indigenous cattle. Divergence estimates between Icelandic cattle and indigenous breeds suggested a separation time of more than 1,000 years between Icelandic cattle and Norwegian native breeds, a finding consistent with historical evidence.     

Skeletal muscle transcriptional profiles in two Italian beef breeds, <Emphasis Type="Italic">Chianina</Emphasis> and <Emphasis Type="Italic">Maremmana</Emphasis>, reveal breed specific variation

Chianina and Maremmana breeds play an important role in the Italian cattle meat market. The Chianina breed is an ancient breed principally raised for draught. Now this breed is the worldwide recognized producer of top quality beef, tasteful and tender, specifically the famous “Florentine steak”. The Maremmana characterized by a massive skeletal structure, is a rustic cattle breed selected for adaptability to the marshy land of the Maremma region. We used a high throughput mRNA sequencing to analyze gene expression in muscle tissues of two Italian cattle breeds, Maremmana (MM) and Chianina (CN) with different selection history. We aim to examine the specific genetic contribution of each breed to meat production and quality, comparing the skeletal muscle tissue from Maremmana and Chianina. Most of the differentially expressed genes were grouped in the Glycolysis/Gluconeogenesis pathways. The rate and the extent of post-mortem energy metabolism have a critical effect on the conversion of muscle to meat. Furthermore, we aim at discovering the differences in nucleotide variation between the two breeds which might be attributable to the different history of selection/divergence. In this work we could emphasize the involvement of pathways of post-mortem energy metabolism. Moreover, we detected a collection of coding SNPs which could offer new genomic resources to improve phenotypic selection in livestock breeding program.     

Contribution of Lidia cattle breed historical castes to the paternal genetic stock of Spain

The main objective of this work was to determine whether the five founding castes defined in the Lidia cattle breed actually have an important contribution to the Spanish paternal genetic stock as well as to the paternal genetic origin support. A total of 1300 Bos taurus male individuals were genotyped for five microsatellites (INRA189, UMN0103, UMN0307, BM861 and BYM1) and one indel (ZFY10). Microsatellite and indel alleles were combined into haplotypes, identifying a total of 38 haplotypes, 11 of them belonging to haplogroup Y1 and 27 to haplogroup Y2. Ten different haplotypes were found in the Lidia cattle breed, with five being exclusive to this breed. Our results agree with different male genetic stocks in the Lidia breed: one hypothetically representing the ancient Iberian bovine genetic stock (Gallardo, Navarra and Cabrera castes and some encastes from Vistahermosa) and a second one that is the result of the more recent breeding strategy of choosing the most aggressive individuals from traditional herds (including some Vistahermosa encastes and the Vazqueña caste). In terms of conservation, it would be better to not consider this breed as a unit but to consider the caste, or even better the encaste, as the target of putative conservation efforts.     

Application of Selection Mapping to Identify Genomic Regions Associated with Dairy Production in Sheep

In Europe, especially in Mediterranean areas, the sheep has been traditionally exploited as a dual purpose species, with income from both meat and milk. Modernization of husbandry methods and the establishment of breeding schemes focused on milk production have led to the development of “dairy breeds.” This study investigated selective sweeps specifically related to dairy production in sheep by searching for regions commonly identified in different European dairy breeds. With this aim, genotypes from 44,545 SNP markers covering the sheep autosomes were analysed in both European dairy and non-dairy sheep breeds using two approaches: (i) identification of genomic regions showing extreme genetic differentiation between each dairy breed and a closely related non-dairy breed, and (ii) identification of regions with reduced variation (heterozygosity) in the dairy breeds using two methods. Regions detected in at least two breeds (breed pairs) by the two approaches (genetic differentiation and at least one of the heterozygosity-based analyses) were labeled as core candidate convergence regions and further investigated for candidate genes. Following this approach six regions were detected. For some of them, strong candidate genes have been proposed (e.g. ABCG2, SPP1), whereas some other genes designated as candidates based on their association with sheep and cattle dairy traits (e.g. LALBA, DGAT1A) were not associated with a detectable sweep signal. Few of the identified regions were coincident with QTL previously reported in sheep, although many of them corresponded to orthologous regions in cattle where QTL for dairy traits have been identified. Due to the limited number of QTL studies reported in sheep compared with cattle, the results illustrate the potential value of selection mapping to identify genomic regions associated with dairy traits in sheep.     

Comparative Study on the Genetic Diversity of GHR Gene in Tibetan Cattle and Holstein Cows

Due to the phenotype-based artificial selection in domestic cattle, the underlying functional genes may be indirectly selected and show decreasing diversity in theory. The growth hormone receptor (GHR) gene has been widely proposed to significantly associate with critical economic traits in cattle. In the present study, we comparatively studied the genetic diversity of GHR in Tibetan cattle (a traditional unselected breed, n = 93) and Chinese Holstein cow (the intensively selected breed, n = 94). The Tibetan yak (n = 38) was also included as an outgroup breed. A total of 21 variants were detected by sequencing 1279 bp genomic fragments encompassing the largest exon 9. Twelve haplotypes (H1～H12) constructed by 15 coding SNPs were presented as a star-like network profile, in which haplotype H2 was located at the central position and almost occupied by Tibetan yaks. Furthermore, H2 was also identical to the formerly reported sequence specific to African cattle. Only haplotype H5 was simultaneously shared by all three breeds. Tibetan cattle showed higher nucleotide diversity (0.00215 ± 0.00015) and haplotype diversity (0.678 ± 0.026) than Holstein cow. Conclusively, we found Tibetan cattle have retained relatively high genetic variation of GHR. The predominant presence of African cattle specific H2 in the outgroup yak breed would highlight its ancestral relationship, which may be used as one informative molecular marker in the phylogenetic studies.     

Signatures of Diversifying Selection in European Pig Breeds

Following domestication, livestock breeds have experienced intense selection pressures for the development of desirable traits. This has resulted in a large diversity of breeds that display variation in many phenotypic traits, such as coat colour, muscle composition, early maturity, growth rate, body size, reproduction, and behaviour. To better understand the relationship between genomic composition and phenotypic diversity arising from breed development, the genomes of 13 traditional and commercial European pig breeds were scanned for signatures of diversifying selection using the Porcine60K SNP chip, applying a between-population (differentiation) approach. Signatures of diversifying selection between breeds were found in genomic regions associated with traits related to breed standard criteria, such as coat colour and ear morphology. Amino acid differences in the EDNRB gene appear to be associated with one of these signatures, and variation in the KITLG gene may be associated with another. Other selection signals were found in genomic regions including QTLs and genes associated with production traits such as reproduction, growth, and fat deposition. Some selection signatures were associated with regions showing evidence of introgression from Asian breeds. When the European breeds were compared with wild boar, genomic regions with high levels of differentiation harboured genes related to bone formation, growth, and fat deposition.     

Footprints of recent selection and variability in breed composition in the Göttingen Minipig genome

The Göttingen Minipig (GMP) developed at the University of Göttingen is a synthetic breed that is widely used in medical research and toxicology. It combines the high fertility of the Vietnamese potbellied pig, the low body weight of the Minnesota Minipig and the white coat colour of the German Landrace pig. The aim of this study was to find genomic regions that may have undergone selection since the creation of the breed in the 1960s. Therefore, the whole genome was screened for footprints of recent selection based on single nucleotide polymorphism (SNP) genotypes from the Illumina Porcine SNP60 BeadChip using two methods: the extended haplotype homozygosity (EHH) test and the estimation of the genomic proportion of the three original breeds at each SNP using a Bayesian approach. Local deviations from the average genome‐wide breed composition were tested with a permutation‐based empirical test. Results for a comprehensive whole‐genome scan for both methods are presented. Several regions showing the highest P‐values in the EHH test are related to breeding goals relevant in the GMP, such as growth (SOCS2, TXN, DDR2 and GRB10 genes) and white colour (PRLR gene). Additionally, the calculated proportion of the founder breeds diverged significantly in many regions from the pedigree‐based expectations and the genome average. The results provide a genome‐wide map of selection signatures in the GMP, which leads to a better understanding of selection that took place over the last decades in GMP breed development.     

Comparative selection signature analyses identify genomic footprints in Reggiana cattle,the traditional breed of the Parmigiano-Reggiano cheese production system

Reggiana is an autochthonous cattle breed reared mainly in the province of Reggio Emilia, located in the North of Italy. Reggiana cattle (originally a triple-purpose population largely diffused in the North of Italy) are characterised by a typical solid red coat colour. About 2500 cows of this breed are currently registered to its herd book. Reggiana is now considered a dual-purpose breed even if it is almost completely dedicated to the production of a mono-breed branded Protected Designation of Origin Parmigiano-Reggiano cheese, which is the main driver of the sustainable conservation of this local genetic resource. In this study, we provided the first overview of genomic footprints that characterise Reggiana and define the diversity of this local cattle breed. A total of 168 Reggiana sires (all bulls born over 35 years for which semen was available) and other 3321 sires from 3 cosmopolitan breeds (Brown, Holstein and Simmental) were genotyped with the Illumina BovineSNP50 panel. ADMIXTURE analysis suggested that Reggiana breed might have been influenced, at least in part, by the other three breeds included in this study. Selection signatures in the Reggiana genome were identified using three statistical approaches based on allele frequency differences among populations or on properties of haplotypes segregating in the populations (fixation index (F_ST); integrated haplotype score; cross-population extended haplotype homozygosity). We identified several regions under peculiar selection in the Reggiana breed, particularly on bovine chromosome (BTA) 6 in the KIT gene region, that is known to be involved in coat colour pattern distribution, and within the region of the LAP3, NCAPG and LCORL genes, that are associated with stature, conformation and carcass traits. Another already known region that includes the PLAG1 gene (BTA14), associated with conformation traits, showed a selection signature in the Reggiana cattle. On BTA18, a signal of selection included the MC1R gene that causes the red coat colour in cattle. Other selection sweeps were in regions, with high density of quantitative trait loci for milk production traits (on BTA20) and in several other large regions that might have contributed to shape and define the Reggiana genome (on BTA17 and BTA29). All these results, overall, indicate that the Reggiana genome might still contain several signs of its multipurpose and non-specialised utilisation, as already described for other local cattle populations, in addition to footprints derived by its ancestral origin and by its adaptation to the specialised Parmigiano-Reggiano cheese production system.     

Polymorphism of Bovine Prolactin Gene: Microsatellites,PCR–RFLP

In the samples of Russian Ayrshire and Gorbatov Red cattle breeds, distribution of frequencies of prolactin (PRL) gene alleles generated due to the presence of polymorphic RsaI site in exon 3 were studied. In the breeds, the frequencies of the Ballele of the PRLgene (with RsaI(+) site) detected by the PCR–RFLP method were 14.1 and 8.6%, respectively. In Black Pied, Ayrshire and Gorbatov Red cattle breeds, variation of the microsatellite dinucleotide repeat in the regulatory region of the gene PRLwas also studied. Gorbatov Red breed was monomorphic at the microsatellite locus with the only allele 164 bp in length. Two alleles (164 bp and 162 bp) were detected in the other breeds studied. The frequencies of 164-bp allele of the microsatellite locus were 93.7 and 90.0% in Black Pied and Ayrshire breeds, respectively. In Gorbatov Red breed of dairy type with good beef qualities and low milk-fat yield, lower level of heterozygosity for PRLgene was demonstrated compared to Ayrshire and Black Pied breeds that have high milk-fat yield. In three cattle breeds, higher mean estimate of polymorphism information content of PCR–RFLP in exon 3 (PIC = 0.21) was revealed compared with the same estimate (PIC = 0.09) for the microsatellite locus variability in the regulatory region of the PRLgene. Characteristics of allele Bdistribution of thePRLgene in the representatives of the Bovidae family are considered.     

Genome‐wide variability and selection signatures in Italian island cattle breeds

In the present study, a sample of 88 animals belonging to four local (Modicana, Sarda, Sardo‐Bruna and Sardo‐Modicana) and one cosmopolitan (Italian Brown Swiss) cattle breeds were genotyped with a medium density SNP beadchip and compared to investigate their genetic diversity and the existence of selection signatures. A total of 43 012 SNPs distributed across all 29 autosomal chromosomes were retained after data quality control. Basic population statistics, Wright fixation index and runs of homozygosity (ROH) analyses confirmed that the Italian Brown Swiss genome was shaped mainly by selection, as underlined by the low values of heterozygosity and minor allele frequency. As expected, local cattle exhibited a large within‐breed genetic heterogeneity. The F_ST comparison revealing the largest number of significant SNPs was Sardo‐Bruna vs. Sardo‐Modicana, whereas the smallest was observed for Italian Brown Swiss vs. Sardo‐Modicana. Modicana exhibited the largest number of detected ROHs, whereas the smallest was observed for Sardo‐Modicana. Signatures of selection were detected in genomic regions that harbor genes involved in milk production traits for Italian Brown Swiss and fitness traits for local breeds. According to the results of multi‐dimensional scaling and the admixture analysis the Sardo‐Bruna is more similar to the Sarda than to the Italian Brown Swiss breed. Moreover, the Sardo‐Modicana is genetically closer to the Modicana than to the Sarda breed. Results of the present work confirm the usefulness of single nucleotide polymorphisms in deciphering the genetic architecture of livestock breeds.     

Characterization of the Russian beef cattle breed gene pools using inter simple sequence repeat DNA analysis (ISSR analysis)

The gene pools of beef cattle breeds bred in Russia were characterized on the basis of inter simple sequence repeat DNA analysis (ISSR analysis). Samples of Aberdeen Angus, Kalmyk, and Kazakh Whitehead breeds from Russia, as well as of Hereford breed, hybrids of Kazakh Whitehead and Hereford breeds, and Kazakh Whitehead breed from the Republic of Kazakhstan, were examined. In the examined breeds, 27 AG-ISSR fragments were identified, 25 of which were polymorphic. The examined breeds were different both in the fragment profiles (the presence/absence of individual ISSR fragments) and in their frequencies. It was demonstrated that the hybrid animals lacked some ISSR fragments that were present with high frequencies in parental forms, suggesting considerable genome rearrangement in the hybrid animals (at the regions of microsatellite localization) in crossings of the individuals from different breeds. The level of genetic diversity in Russian beef breeds was consistent with the values typical of farmed populations (breeds). The genetic diversity parameters assessed by applying Nei’s gene diversity index and the Shannon index varied from 0.0218 to 0.0605 and from 0.0225 to 0.0819, respectively. The highest Shannon index value was detected in the Kalmyk breed (0.0837) and Kazakh Whitehead breed from Russia (0.0819), and the highest level of Nei’s gene diversity index was found in the Kalmyk breed (0.0562) and in both populations of the Kazakh Whitehead breed (0.0509 and 0.0605). The high level of genetic similarity (according to Nei) was revealed between Russian beef cattle breeds and Hereford cattle: 0.839 (for the Kazakh Whitehead breed from Russia) and 0.769 (for the Kalmyk breed).     

Diversity of the TLR4 Immunity Receptor in Czech Native Cattle Breeds Revealed Using the Pacific Biosciences Sequencing Platform

The allelic variants of immunity genes in historical breeds likely reflect local infection pressure and therefore represent a reservoir for breeding. Screening to determine the diversity of the Toll-like receptor gene TLR4 was conducted in two conserved cattle breeds: Czech Red and Czech Red Pied. High-throughput sequencing of pooled PCR amplicons using the PacBio platform revealed polymorphisms, which were subsequently confirmed via genotyping techniques. Eight SNPs found in coding and adjacent regions were grouped into 18 haplotypes, representing a significant portion of the known diversity in the global breed panel and presumably exceeding diversity in production populations. Notably, the ancient Czech Red breed appeared to possess greater haplotype diversity than the Czech Red Pied breed, a Simmental variant, although the haplotype frequencies might have been distorted by significant crossbreeding and bottlenecks in the history of Czech Red cattle. The differences in haplotype frequencies validated the phenotypic distinctness of the local breeds. Due to the availability of Czech Red Pied production herds, the effect of intensive breeding on TLR diversity can be evaluated in this model. The advantages of the Pacific Biosciences technology for the resequencing of long PCR fragments with subsequent direct phasing were independently validated.