Abstract. A production herd of Czech Simmental cattle (Czech Red Pied, CRP), the conserved subpopulation of this breed, and the ancient local breed Czech Red cattle (CR) were screened for diversity in the antibacterial toll-like receptors (TLRs), which are members of the innate immune system. Polymerase chain reaction (PCR) amplicons of TLR1, TLR2, TLR4, TLR5, and TLR6 from pooled DNA samples were sequenced with PacBio technology, with 3–5× coverage per gene per animal. To increase the reliability of variant detection, the gDNA pools were sequenced in parallel with the Illumina X-ten platform at low coverage (60× per gene). The diversity in conserved CRP and CR was similar to the diversity in conserved and modern CRP, representing 76.4 % and 70.9 % of its variants, respectively. Sixty-eight (54.4 %) polymorphisms in the five TLR genes were shared by the two breeds, whereas 38 (30.4 %) were specific to the production herd of CRP; 4 (3.2 %) were specific to the broad CRP population; 7 (5.6 %) were present in both conserved populations; 5 (4.0 %) were present solely for the conserved CRP; and 3 (2.4 %) were restricted to CR. Consequently, gene pool erosion related to intensive breeding did not occur in Czech Simmental cattle. Similarly, no considerable consequences were found from known bottlenecks in the history of Czech Red cattle. On the other hand, the distinctness of the conserved populations and their potential for resistance breeding were only moderate. This relationship might be transferable to other non-abundant historical cattle breeds that are conserved as genetic resources. The estimates of polymorphism impact using Variant Effect Predictor and SIFT software tools allowed for the identification of candidate single-nucleotide polymorphisms (SNPs) for association studies related to infection resistance and targeted breeding. Knowledge of TLR-gene diversity present in Czech Simmental populations may aid in the potential transfer of variant characteristics from other breeds.

Hybrid resequencing of the antibacterial innate immune genes coding for toll-like receptors, namely TLR1, TLR2, TLR4, TLR5, and TLR6, using HiSeq and PacBio technologies of pooled population samples of Czech Simmental (Czech Red) cattle allowed us to determine haplotypes formed by the polymorphisms present. Directly determined haplotypes within the range of the large proximal amplicon in TLR2 formed two clusters in the network tree graph. The distribution of the statistically reconstructed haplotypes based on individual genotyping of the present SNPs was consistent. Similarly, the statistically reconstructed haplotypes in TLR5 and TLR6 formed two clusters. The trend of bimodal distribution was also observed in TLR4, while the limited diversity of TLR1 did not allow for any conclusion. The observed bimodal distribution is consistent with earlier reports for cattle populations worldwide. The stability of this phenomenon cannot be ascribed to historical origin but rather to a long-term effect of balancing selection. The equilibrium might be based on two different essential functions performed by the TLR genes or their products. The formation of two kinds of heterodimers by the TLR2 product, namely, TLR2/TLR1 and TLR2/TLR6 with different ligand specificities, is considered to be a particular case. On the other hand, the better expression of the bimodal groups in the 5′-proximal SNPs supports the localization of the selection targets in the upstream regulatory regions or the functional interactions in the proximal part of the transcripts.

Driven by co-evolution with pathogens, host immunity continuously adapts to optimize defence against pathogens within a given environment. Recent advances in genetics, genomics and transcriptomics have enabled a more detailed investigation into how immunogenetic variation shapes the diversity of immune responses seen across domestic and wild animal species. However, a deeper understanding of the diverse