Single nucleotide polymorphism discovery in bovine liver using RNA-seq technology

Background RNA-seq is a useful next-generation sequencing (NGS) technology that has been widely used to understand mammalian transcriptome architecture and function. In this study, a breed-specific RNA-seq experiment was utilized to detect putative single nucleotide polymorphisms (SNPs) in liver tissue of young bulls of the Polish Red, Polish Holstein-Friesian (HF) and Hereford breeds, and to understand the genomic variation in the three cattle breeds that may reflect differences in production traits. Results The RNA-seq experiment on bovine liver produced 107,114,4072 raw paired-end reads, with an average of approximately 60 million paired-end reads per library. Breed-wise, a total of 345.06, 290.04 and 436.03 million paired-end reads were obtained from the Polish Red, Polish HF, and Hereford breeds, respectively. Burrows-Wheeler Aligner (BWA) read alignments showed that 81.35%, 82.81% and 84.21% of the mapped sequencing reads were properly paired to the Polish Red, Polish HF, and Hereford breeds, respectively. This study identified 5,641,401 SNPs and insertion and deletion (indel) positions expressed in the bovine liver with an average of 313,411 SNPs and indel per young bull. Following the removal of the indel mutations, a total of 195,3804, 152,7120 and 205,3184 raw SNPs expressed in bovine liver were identified for the Polish Red, Polish HF, and Hereford breeds, respectively. Breed-wise, three highly reliable breed-specific SNP-databases (SNP-dbs) with 31,562, 24,945 and 28,194 SNP records were constructed for the Polish Red, Polish HF, and Hereford breeds, respectively. Using a combination of stringent parameters of a minimum depth of ≥10 mapping reads that support the polymorphic nucleotide base and 100% SNP ratio, 4,368, 3,780 and 3,800 SNP records were detected in the Polish Red, Polish HF, and Hereford breeds, respectively. The SNP detections using RNA-seq data were successfully validated by kompetitive allele-specific PCR (KASPTM) SNP genotyping assay. The comprehensive QTL/CG analysis of 110 QTL/CG with RNA-seq data identified 20 monomorphic SNP hit loci (CARTPT, GAD1, GDF5, GHRH, GHRL, GRB10, IGFBPL1, IGFL1, LEP, LHX4, MC4R, MSTN, NKAIN1, PLAG1, POU1F1, SDR16C5, SH2B2, TOX, UCP3 and WNT10B) in all three cattle breeds. However, six SNP loci (CCSER1, GHR, KCNIP4, MTSS1, EGFR and NSMCE2) were identified as highly polymorphic among the cattle breeds. Conclusions This study identified breed-specific SNPs with greater SNP ratio and excellent mapping coverage, as well as monomorphic and highly polymorphic putative SNP loci within QTL/CGs of bovine liver tissue. A breed-specific SNP-db constructed for bovine liver yielded nearly six million SNPs. In addition, a KASPTM SNP genotyping assay, as a reliable cost-effective method, successfully validated the breed-specific putative SNPs originating from the RNA-seq experiments.


Results
mRNA sequencing and read alignment mRNA sequencing of bovine liver at single-nucleotide resolution was carried out using two biological replicates of poly(A)-enriched mRNA of young bulls aged 6, 9, and 12 months from three cattle breeds. These mRNA samples were first converted into barcoded strand-specific dUTP RNA-seq libraries, followed by HT sequencing on the Illumina NextSeq 500 sequencer. The HT sequencing produced a total of 107,114,4072 raw paired-end reads with a length of 156 bases. The reads were de-multiplexed to assign reads to each sequenced sample according to its index. The FASTQ sequence dataset of each library (Table 1) was submitted to the NCBI-SRA database with NCBI-SRA experiment number SRS1296732 (https://www.ncbi. nlm.nih.gov/sra?linkname=bioproject_sra_all&from_uid=312148).
Using the Burrows-Wheeler Aligner (BWA) program under default conditions, the breedspecific liver transcripts were mapped to the bovine reference genome (UMD3.1 assembly plus Y chromosome). Our read alignment results showed that 98.55% sequencing reads (105,567,9630) were successfully aligned to the UMD3.1 bovine reference genome. Furthermore, the acquired BWA mapping profile results showed that 81.35%, 82.81% and 84.21% of the mapped sequencing reads were properly paired for the Polish Red, Polish HF, and Hereford breeds, respectively (Tables 2-4).

SNP discoveries in cattle breeds
Breed-specific raw SNP-db records. Using the SAMtool package, a total of 5,641,401 (~5.6 million) breed-specific SNPs and indel positions expressed in the bovine liver were detected with the RNA-seq reads, with an average of 313,411 (~0.31 million) SNPs and indels per young bull (Table 5). Breed-wise, this raw SNP-db comprised 1,995,571 (35.4%), 1,556,048 (27.6%), and 2,089,782 (37%) SNPs and indels for the Polish Red, Polish HF, and Hereford breeds, respectively. Following the removal of the indel mutations, a total of 1,953,804 (35.3%), 1,527,120 (27.6%), and 2,053,184 (37.1%) raw SNPs expressed in bovine liver were recovered from the Polish Red, Polish HF, and Hereford breeds, respectively. Using the SAMtool package, single-base substitutions (SNPs) and small indels were also identified. In this study, a total of 41,767, 36,604 and 28,934 indel mutations were identified in the Polish Red, Polish HF, and Hereford breeds, respectively. SNP distribution in Venn plot. In our initial SNPs analysis, a stringent filtering parameter of read count with a minimum depth of !5 SNP reads that support the polymorphic nucleotide base and existed in both replicates was utilized to allow the identification of approximately 0.8 million SNPs among the three cattle breeds as shown in Venn diagram (Fig 1).
Breed-specific SNP-db records. For the detection of breed-specific putative SNPs expressed in bovine liver, only the records in the raw SNP-db (S1-S18 Tables) were combined into one file to construct a highly reliable SNP-db with 84,701 SNP hit records. Three highly reliable breed-specific SNP-dbs comprising 31,562 (37.27%), 24,945 (29.45%) and 28,194  Error removals: During SNP detection using SAMtool, some records were observed as more than one SNP mutation. Such records were considered as error records and were excluded from the SNP-db. In our study, a total of 4381, 1164 and 1202 such error records in the Polish Red, Polish HF and Hereford breeds were observed and excluded. After removal of the error records, a total of 27,182, 23,781 and 26,992 SNP-db records were recovered from the Polish Red, Polish HF and Hereford breeds, respectively (S22-S24 Tables).
SNP filtering. We utilized the stringent parameter of a minimum depth of 10 SNP reads that support the polymorphic nucleotide base with a SNP ratio of 100%, because the SNP filtering criteria of !10 SNP reads with a SNP ratio of 100% could cover and explain all the HT SNP variations compared with the SNP filtering criteria of !10 SNP reads with a SNP ratio of 90% [9].
Initially, for the SNP filtering analysis, we utilized stringent parameters with a minimum depth of !10 SNP reads that support the polymorphic nucleotide base and identified 15,197,11,346 and 12,455 SNP records for the Polish Red, Polish HF, and Hereford breeds, respectively (S25-S27 Tables). Similarly, by utilizing stringent parameters of 100% SNP ratios, we identified 10,206, 9,684 and 9,778 SNP records for the Polish Red, Polish HF, and Hereford breeds, respectively (S28-S30 Tables). Finally, the combination of both stringent filtering parameters of !10 reads, and a 100% SNP ratio, yielded a total of 4,368, 3,780 and 3,800 SNPs records for the Polish Red, Polish HF, and Hereford breeds, respectively (S31-S33 Tables).
Breed comparison. Comparison of the breed-specific SNP records resulted in the identification of common and unique SNPs. For each breed, three comparisons were made to find common and unique SNPs within the investigated breeds. For example, i) a comparison between breed-1 (Polish Red) with breed-2 (Polish HF); ii) a comparison between breed-1 and breed-3 (Hereford); and iii) a comparison of breed-1 to both breed-2 and breed-3.
In this study, a total of 50 and 81 unique SNP loci were identified in the Polish Red breed that were not present in either the Polish HF or Hereford breeds, respectively (S34 and S35 Tables). In addition, only six unique SNP records were identified in the Polish Red breed that was not detected in either Polish HF or Hereford breeds (S36 Table).
Similarly, the SNP data comparison of the Polish HF breed to other breeds resulted in the identification of 7 and 22 unique SNPs that were not present in the Polish Red or Hereford breeds, respectively (S37 and S38 Tables, respectively). However, no single unique SNP record was identified in the Polish HF breed that was not present in either the Polish Red or Hereford breeds, respectively (S39 Table). Finally, the SNP data comparison of the Hereford breed to other breeds identified a total of 80 and 41 unique SNPs that were not detected in Polish Red or Polish HF breeds, respectively (S40 and S41 Tables). Furthermore, a total of nine unique SNP records of the Hereford breed were identified, which were not present in either the Polish Red or Polish HF breeds (S42 Table).
De novo SNPs. Deep and extensive SNP analysis provided evidence that certain regions of the Bos taurus genome were still unknown (base = N, according to recent mapping to UMD3.1). Using the SNP filtering criteria (!10 SNP reads with SNP ratio of 100%), a total of 217, 193 and 265 best candidates of de novo SNP reads were identified for the Polish Red, Polish HF, and Hereford breeds, respectively (S43-S45 Tables).

Breed-specific SNP discovery and QTL/CG analysis
For the QTL/CG analysis on identified breed-specific SNP data, a set of 110 QTL/CG loci (http://www.animalgenome.org/cgi-bin/QTLdb/index) (S46 Table)  investigated for each breed and ages of the young bulls (S47-S52 Tables). We summarized, for each breed, the chromosomal locations and SNP locations of identified putative SNP loci hits of RNA-seq data on 110 potential CGs from the bovine QTL-db (S47-S49 Tables). In S50-S52 Identification of highly polymorphic breed-specific SNP loci in the bovine QTL/CG db. Highly polymorphic breed-specific SNP loci were summarized as top-10 and top-20 SNP hits (S50-S52 Tables). Among the top 10 SNP hits, five CGs (CCSER1, GHR, KCNIP4, MTSS1, EGFR and NSMCE2) were highly polymorphic to all three cattle breeds. Furthermore, SNP hits in the IGF-I CG were highly polymorphic to the Polish Red and Polish HF breeds, SNP hits in the CAST CG were highly polymorphic to the Polish HF and Hereford breeds, and SNP hits in the LEPR CG were highly polymorphic to the Polish Red and Hereford breeds. In addition, SNP hits in the IGF2R and MAP3K5 CG loci were highly polymorphic to the Polish Red breed, SNP hits in the INSIG1 and ATP6V1B2 CG loci were highly polymorphic to the Polish HF breed, and SNP hits in the PRLR, IGFBP2 CG loci were highly polymorphic to the Hereford breed.  Phylogenetic analysis of the breed-specific SNP-db Additionally, a phylogenetic tree was constructed using a JTT matrix-based model with superior log maximum likelihood values [28] and molecular evolutionary genetics analysis version 7 (MEGA7) software [29] to examine the close relationship of the investigated bulls within and between breeds. Results showed that all three cattle breeds were separated from each other, but clustered together for each breed. The percentage of trees in which the associated taxa clustered together is shown next to the branches (Fig 2).

Breed-specific SNP validation
Selection of nine breed-specific putative SNPs. A breed-specific SNP validation experiment was carried out using a subset of nine putative SNPs derived from the Polish Red (n = 2), Polish HF (n = 4) and Hereford (n = 3) breeds, on single-plex KASP TM genotyping assays (LGC Genomics) based on fluorescently labelled allele-specific PCR primers ( Table 6). The selection of nine putative breed-specific SNPs was performed according to the RNA-seq experimental results (S53 Table). We selected three SNPs that were uniquely expressed in only one cattle breed (CTNS and P4HA2 SNP gene loci specific to the Polish Red breed and IQGAP2 SNP gene locus specific to the Hereford breed), and the remaining six SNP gene loci that were specific to all the breeds (GHR, IGF2R, IGF2BP3, IGFBP4 SNP gene loci specific to the Polish HF breed, and GHR and IGF2 SNP gene loci specific to the Hereford breed, respectively). Based on the RNA-seq experimental results, detailed information of the primer design of the nine selected breed-specific SNP loci, such as the candidate genes and their symbols, genome   Table. The selected subset of the nine breed-specific putative SNP markers from the Polish Red, Polish HF and Hereford breeds worked well in KASP TM SNP genotyping assay and did not reveal either non-amplification or ambiguous clustering, except for a few samples due to poor DNA quality (Figs 3-11).
KASP TM SNP assay analysis. Initially, an additional analysis for the estimation of a direct relationship between breed and the detected SNP polymorphisms was performed using the PROC MIXED SAS 9.2 package with age as a random effect. Two SNP markers, IGFBP4 (Polish HF) and IGF2 (Hereford) genes were identified as homozygous at the SNP loci, and therefore were excluded for further statistical SNP validation analysis (S55 Table). The results of the remaining seven SNP loci did not show any significant association between age and SNPs selected for the validation. However, there were significant and highly significant associations   Table). Furthermore, a chi-square analysis showed significant differences in genotypes and allele frequencies for the CTNS gene SNP locus specific to the Polish Red breed, IGF2R gene SNP locus specific to the Polish HF breed, and the GHR and IQGAP2 genes SNP loci specific to the Hereford breed (S57 Table). Statistical analysis using Genepop software. Using the Fisher's Exact Probability test, the genetic differentiation of SNP alleles and SNP genotype results revealed significant differences in SNP allele frequencies for the CTNS, IGF2R, GHR (Polish HF breed), IGF2BP3, GHR (Hereford breed) and IQGAP2 SNP loci (S58 Table). Regarding the genotype frequencies, significant differences were observed in the CTNS, P4HA2, IGF2R, GHR (Polish HF breed), IGF2BP3, GHR (Hereford breed), and IQGAP2 SNP loci (S59 Table) in all investigated cattle breeds. Moreover, the genetic differentiation comparison of the SNP alleles (S60 Table) and genotypes (S61 Table) in all SNP loci and all investigated cattle breeds were also performed, using Fisher's Exact Probability test and Fisher's Exact G test. Results showed highly significant differences in the SNP alleles (S60 Table) and SNP genotypes (S61 Table) among breeds to validate the breed-specific SNP markers. For SNP validation statistics using a Markov chain method, the selected SNP markers were further examined by testing the deviation from Hardy Weinberg equilibrium (HWE) for each SNP locus (S62 Table), and for each breed population (S63 Table). The results based on genetic differentiation of the investigated SNPs across all loci showed significant differences in the CTNS, IGF2R, IGF2BP3, GHR (Polish HF), and IQGAP2 SNP loci (S62 Table). Furthermore, the results based on genetic differentiation of the investigated SNPs across all breeds revealed significant differences among them (S63 Table). Discussion mRNA sequencing and read alignment RNA-seq technology has great potential in identifying genetic variation at many loci, with respect to SNP polymorphisms and gene expression patterns across different organ tissues. In the present study, we identified approximately 107 million raw paired-end reads with an average length of 156 bases in bovine liver tissue compared to 113 million raw paired-end reads detected in the bovine pituitary gland [9]. This variation in the RNA-seq yield from different organ tissues might be due to technical variation, such as differences in the quality and quantity of the RNA recovered during tissue-specific sample preparations, batch effects in library preparation [30,31], flow cell and lane effects caused by the Illumina sequencing platform, or adapter bias [32,33]. However, both the liver and pituitary gland tissues revealed higher percentages (99.4% and 98.5%, respectively) of read alignment to the UMD3.1 bovine reference genome. Furthermore, using the BWA program under default conditions, showed that 94.39%, 93.04% and 83.46% of the mapped sequencing reads of the pituitary gland tissues were properly paired for the Polish Red, Polish HF, and Hereford breeds, respectively [9]. The FASTQ sequences of both tissues were submitted with the NCBI-SRA experiment number SRS1296732 (http://www.ncbi.nlm.nih.gov/sra?linkname=bioproject_sra_all&from_ uid=312148). Currently, there are 24 SRA records (FASTQ sequence dataset) for bovine liver tissue at the SRA NCBI database [34] (https://trace.ncbi.nlm.nih.gov/Traces/sra/sra.cgi?view= studies&term=(bos%20taurus%20liver)%20NOT%20cluster_dbgap%5BPROP%5D)).

SNP discovery in cattle breeds
One of the potential applications of the RNA-seq technology is to identify several thousands of SNPs and as well as to construct the tissue-specific SNP datasets. In our study, we constructed breed-specific SNP-db records for bovine liver tissue, with more than 5.6 million SNP records compared to 13.7 million SNP-db records for the bovine pituitary gland tissue [9]. The differences in the SNP data yield in tissue-specific RNA-seq experiments in the bovine liver and pituitary gland might be due to pre-and post-quality assessment of read alignments, i.e., quality control (QC) metrics before alignment [35]. In most SNP discoveries based on bovine RNA-seq studies [7][8][9][16][17][18][19][20], selection of stringent parameters in the SNP filtering process is undoubtedly a critical task for identifying the most reliable novel putative SNPs. Using the stringent parameter of a minimum depth of 10 SNP reads that support the polymorphic nucleotide base with a SNP ratio of 100%, we identified 20,573, 31,978 and 30,052 breed-specific SNP records for bovine pituitary gland in the Polish Red, Polish HF and Hereford breeds, respectively [9]. However, in the current study, we identified 31,562, 24,945 and 28,194 breedspecific SNP records for bovine liver tissue in the Polish Red, Polish HF and Hereford cattle breeds, respectively.
It should be emphasized that, in cattle breeding practices, identification of breed-specific gene-associated SNPs can serve as suitable markers for trait-associated studies and can be effectively utilized in genomic selection (GS) programmes [36]. In our study, we investigated 110 QTL/CG loci to identify novel putative gene-associated SNPs within a breed, with respect to the animal's age, and between breeds. Results based on within-breed QTL/CG analysis revealed that approximately one third of the selected QTL/CG loci were monomorphic in all investigated breeds. These findings indicate that the dairy and beef cattle breeds might have certain genetic selection signatures due to fixation of certain genotypes (monomorphic SNPs). It is noteworthy that, between breeds, the QTL/CG analysis identified CCSER1, GHR, KCNIP4, MTSS1, EGFR and NSMCE2 gene loci as highly polymorphic. Moreover, analysis of the 76 QTL/CGs in bovine pituitary gland RNA-seq data demonstrated that KCNIP4, CCSER1, DPP6, MAP3K5 and GHR genes loci were highly polymorphic in all the investigated cattle breeds [9]. These results indicate that there is a still an abundant genetic variation in the Polish dairy and beef cattle breeds due to the high degrees of polymorphism in the CGs of economically important traits, which could be exploited in GS and other breeding programmes.
In cattle, many studies have utilized QTL/CG analysis with respect to whole genome sequence (WGS) data, genome-wide association study (GWAS) data and expression quantitative trait loci (eQTL) data, to investigate the causal relationships of variant-trait and variantgene expression, respectively [37][38][39][40]. Using human and bovine data, Seo et al. [5] have successfully demonstrated the application of the association test approach based on RNA-seq analysis to identify trait associated genes (TAGs).
In this study, we have validated and revealed that breeds of the investigated bulls were separated from each other using an ITT matrix based maximum likelihood model [28] and MEGA 7 [29]. With regards to bovine pituitary gland tissue, we successfully utilized another phylogenetic analysis method using the SNPphylo model [41], which revealed that the three cattle breeds were clustered together for each breed and were separated from each other [9].
Finally, the utilization of KASP TM genotyping assay technology in our SNP validation experiment has confirmed the presence of seven breed-specific SNP markers, which are expressed in bovine liver tissue. In addition, a similar SNP experiment has been successfully performed to validate the expression of the seven SNP markers in bovine pituitary tissue [9]. Thus, the SNP experiment using KASP TM genotyping assay technology is an effective and reliable method to validate RNA-seq in various bovine organ tissues.

Conclusions
Breed-specific SNP discovery using NGS-based RNA-seq in bovine liver tissue has been reported to provide a global view of the complexity of the bovine liver transcriptome. Our study has developed three breed-specific SNP-dbs based on expressed genes in the bovine liver, which might provide valuable resources for trait-associated genomic and genome-wide association studies. Our study has demonstrated the utility of QTL/CG analysis on RNA-seq SNP data to identify putative trait-associated SNPs from the bovine QTL-db. Our phylogenetic analyses have shown that all three cattle breeds were separated from each other long ago and that each breed is represented as unique genetic clusters. The transcriptome sequencing (RNA-seq) technique employed in the current study is similar to the Genotyping-by-Sequencing (GBS) and it has good potential for improving the accuracy of genomic selection because it involves several millions of SNPs covering the entire transcriptome, which increases the chance of identifying QTLs or quantitative trait nucleotides (QTNs) in linkage disequilibrium with SNPs. Furthermore, these results suggest that there are wide differences in the liver transcriptomes between the breeds, which could be useful to study the mechanisms underlying genetic variability in meat quality and other production traits in cattle.

Experimental design
The bovine liver tissue samples were collected from 18 young bulls of three cattle breeds ( Table 6) stationed at the Institute of Genetics and Animal Breeding (IGAB), Polish Academy of Science (PAS), Jastrzębiec, Poland. The investigated animals were selected randomly, and after the purchase at birth, they were reared in the experimental farm of IGAB, PAS, Jastrzębiec, Poland, with uniform environmental and feeding conditions. After slaughtering, the collected liver tissues were immediately kept in liquid nitrogen, and stored at -80˚C. All procedures involving animals were conducted in accordance with the guiding principles for the care and use of research animals. The investigating research materials were approved by the local ethics commission of IGAB, PAS, Jastrzębiec, Poland (permission No. 3/2005). The experimental designs, with detailed laboratory procedures including isolation of total RNA, library preparation and NGS sequencing using the Illumina NextSeq 500 High Output/300 cycle platform and kits (Illumina), were performed as previously described in a recent study [9].

Bioinformatics analysis
Complete bioinformatics analysis, including read alignment of RNA-seq data to the reference genome, breed-specific SNP detections, QTL/CG analysis, phylogenetic analysis, SNP validation by KASP TM genotyping assay, was performed as previously described in a recent study [9]. However, i) the SNP filtering criterion, ii) the number of QTL/CG loci to investigate the RNA-seq SNP-db, and iii) SNP markers selected for SNP validation by KASP TM genotyping assay were slightly modified from our recent bovine pituitary gland RNA-seq study [9]. SNP filtering criterion. To construct the bovine liver SNP-db, stringent parameters of a minimum depth of 10 SNP reads that support the polymorphic nucleotide base with a SNP ratio of 100%, were carried out using Microsoft office excel program in the following manner: • Stringent parameter of SNP filtering with minimum depth of 10 SNP reads that support the polymorphic nucleotide base, • Stringent parameter of SNP filtering with a SNP ratio of 100%, • Combining the stringent parameters of a minimum depth of 10 SNP reads that support the polymorphic nucleotide base with a SNP ratio of 100%.
Using Microsoft office excel program, the bovine liver SNP-dbs of the 18 young bulls were combined according to breed and were further trimmed to one SNP-db that was highly reliable and specific to the three breeds.

Breed-specific SNP discovery and QTL/CG analysis
Based on a publicly available animal QTL database (http://www.animalgenome.org/cgi-bin/ QTLdb/index), a total of 110 potential QTL/CGs (S46 Table) for bovine body growth and developmental trait were included to investigate the RNA-seq SNP-db of bovine liver using Microsoft Office Excel.
Phylogenetic analysis. The phylogenetic analysis of the breed-specific SNP-dbs of bovine liver transcriptome was inferred by using the Maximum Likelihood method based on the JTT matrix-based model [28]. The phylogenetic tree(s) of the 18 bull samples representing the Polish Red, Polish HF and Hereford breeds were constructed, firstly with heuristic search by applying Neighbor-Join and BioNJ algorithms to a matrix of pairwise distances estimated using a JTT model, and then by selecting the topology with superior log likelihood value. Finally, the evolutionary analyses were conducted in MEGA7 [29].

SNP validation by KASP TM genotyping assay
The primer sequences of nine breed-specific SNPs from the Polish Red (n = 2), Polish HF (n = 4), Hereford (n = 3) breeds selected for KASP TM SNP genotyping assay to test the SNP validation experiment are presented in Table 7. The complete methodological procedure of the KASP TM genotyping assay, including statistical analysis using Genepop software (http:// genepop.curtin.edu.au) was performed as previously described in a recent study [9].
Supporting information S1   Table. RNA-seq SNP-db of bovine liver tissue of young bull-9 of Polish HF cattle aged 9 months.
(XLSX) S10 Table. RNA-seq SNP-db of bovine liver tissue of young bull-10 of Polish HF cattle aged 9 months. (XLSX) S11 Table. RNA-seq SNP-db of bovine liver tissue of young bull-11 of Polish HF cattle aged 12 months.
(XLSX) S12 Table. RNA-seq SNP-db of bovine liver tissue of young bull-12 of Polish HF cattle aged 12 months.