Microarray analysis of differential gene expression profiles in blood cells of naturally BLV-infected and uninfected Holstein-Friesian cows.

The aim of the present study was to examine gene expression changes in response to bovine leukemia virus (BLV) infection, in an effort to determine genes that take a part in molecular events leading to persistent lymphocytosis (PL), and to better define genes involved in host response to BLV infection. Using bovine 70-mer oligonucleotide spotted microarrays (BLOPlus) and qRT-PCR validation, we studied global gene expression profiles in blood cells in vivo of 12 naturally BLV-infected Polish Holstein cows, and 12 BLV non-infected controls of the same breed and reared in herds with high BLV seroprevalence. With an arbitrary cut-off value of 1.5-fold change in gene expression, we identified the down-regulation of 212 genes (M value ≤-0.585) and the up-regulation of 158 genes (M value of ≥0.585) at 1% false discovery rate in BLV-positive animals in comparison to the BLV-negative group. The gene set enrichment analysis demonstrated that the differentially expressed (DE) genes could be classified to diverse biological processes, including immune response of host blood cells. Interestingly, our data indicated the potential involvement of the innate immunity, including complement system activation, NK-cell cytotoxicity and TREM-1 signaling, during the BLV-induced pathogenesis. We showed the occurrence of numerous regulatory processes that are targeted by BLV-infection. We also suggest that a complex network of interrelated pathways is disturbed, causing the interruption of the control of B-cell proliferation and programmed cell death.

Polymorphism and expression of the tumor necrosis factor receptor II gene in cows infected with the bovine leukemia virus.

A single T>C nucleotide polymorphism (rs42686850) of bovine tumor necrosis factor receptor type II gene (TNF-RII) is located within a sequence with allele-specific affinity to bind E2F transcription factors, considered pivotal in the regulation of cell cycle and cell proliferation. The objective of the study was to determine the effect of this SNP and BLV infection on the TNF-RII gene expression at the mRNA and protein levels in peripheral blood mononuclear cells (PBMC). We noted that analyzed TNF-RII gene polymorphism influenced the expression of the TNF-RII gene at the mRNA level but only in BLV-positive cows. Concurrently, no statistically significant association was found between gene polymorphism and TNF-RII expression at the protein level. However, we found a significant effect of BLV infection status on the amount of TNF-RII mRNA and the percentage of PBMC expressing TNF-RII. These results show an unclear effect of considered T>C polymorphism on TNF-RII gene expression in bovine leukocytes and they suggest the involvement of BLV in modifying the TNF-RII expression in BLV-infected cows potentially implying the EBL (Enzootic Bovine Leukosis) associated pathogenesis.

Association of polymorphism within LTF gene promoter with lactoferrin concentration in milk of Holstein cows.

This study analyzed the association between single nucleotide polymorphism (A/C) in position -28 located in the TATA box of LTF gene and the lactoferrin concentration in bovine milk secreted by healthy and infected udders. Out of 241, 69 cows were selected into the experimental group and were divided into 3 groups according to mean value of somatic cell count (SCC): I < 180,000 cells/mL, II: 180,000-350,000 cells/mL and III > 350,000 cells/mL. In each SCC group, three LTF genotypes: AA, AC and CC were identified by PCR-SSCP method. A total of 604 milk samples were collected monthly and lactoferrin concentration was measured by ELISA. The 1-way ANOVA within SCC groups was performed to estimate association of -28 A/C genotypes with mean lactoferrin concentration per lactation. In the group of healthy cows (< 180,000 cells/mL) LTF concentration in milk cows with the AA genotype (107.58 ± 17.92 µg/mL) was significantly higher than in homozygotes CC (52.09 ± 19.01 µg/mL). Unexpectedly, in cows with elevated SCC (> 350,000 cells/mL) we observed a significant opposite relationship (207.21 ± 28.50 in CC vs 115.0 ± 28.6 µg/mL in AA). We hypothesized that a promoter with allele C, which cannot be recognized as a TATA sequence is becoming more accessible for other transcription factors, which may induce alternative LTF gene expression. We assume that our results demonstrate a very interesting effect of differential gene expression depending on polymorphism in a key regulatory motif (TATA box) and also on the health status of mammary tissues.

[Associations between 60 SNPs identified by APEX microarray and growth rate, meatiness and selection index in boars].

A total of 312 boars (201 Landrace and 111 Large White) were genotyped with a custom-made low throughput genotyping microarray (called SNiPORK) based on array primer extension (APEX) technology. The results were used to association studies between genotyped SNP markers and daily gains, meat content and selection index. Among the 60 SNP markers analyzed, 14 of them showed statistically significant associations between the genotype and the level of at least one trait. In order to find extremely beneficial or unfavorable intergenic diplotype combinations, 5 SNP markers were selected: CASTA499C, MYF6 T255C, PKLR C384T, SFRSI C1146T and TNNT3T 153C, which showed statistically significant associations at P

Screening 52 single nucleotide polymorphisms for extreme value of glycolytic potential and drip loss in pigs.

Summary Low heritability of meat quality traits and the lack of their systematic registration in breeding programs have encouraged the search for single nucleotide polymorphisms (SNPs) located within genes coding the proteins involved in muscle and fat metabolism. In this report, a panel of 52 SNPs was used to find which alleles and genotypes are more/less frequent in groups of pigs differentiated by extreme value of glycolytic potential (GP) and drip loss (DL). The analysis was carried out in 52 fatteners (chosen from 246 pigs), of which 28 were Landrace and 27 Landrace x Yorkshire. Two designs were performed: I, fatteners were divided into two groups showing extreme value of GP (<125 versus >145), II, fatteners were divided into two groups showing extreme value of DL (<6.0 versus >6.0). Allele frequency differences between the phenotypic groups of extreme GL or DL were not influenced by the breed. The frequency of 52 SNPs alleles for each of group was calculated and a chi-squared test was used to estimate the significance of differences in allele frequencies between alternative groups in each experimental design. Three SNPs (DECR1, PPARGC1, MC4R) and another two (CYP21, SFRS1) showed significant differences between groups of extreme GP and DL, respectively. To exemplify and validate potential associations of candidate SNPs for GP and DL, 293 fatteners representing three commercial breeds/crosses (95 Landrace, 66 Landrace x Yorkshire and 132 Landrace x Yorkshire x Duroc were genotyped for DECR1 and CYP21 by PCR-RFLP assays. DECR1 showed significant associations with GP in Landrace and Landrace x Yorkshire x Duroc fatteners. CYP21 showed significant associations with DL in all breeds/crosses. Interestingly, the CYP21 polymorphism revealed adverse associations trend in Landrace x Yorkshire x Duroc pigs in comparison to Landrace and Landrace x Yorkshire fatteners.

Transcriptome analysis of boar spermatozoa with different freezability using RNA-Seq.

Semen freezability is associated with genetic markers, and there is a diverse set of sperm transcripts that have been attributed to various cellular functions. RNA-Seq was performed to compare the transcript profiles of spermatozoa from boars with different semen freezability. We examined ejaculates from the Polish large white (PLW) boars that were classified as having good and poor semen freezability (GSF and PSF, respectively; n = 3 boars per group) by assessing post-thaw motility characteristics, mitochondrial membrane potential, plasma membrane and acrosome integrity. Total RNA was isolated from fresh spermatozoa from boars of the GSF and PSF groups and subjected to RNA-Seq (Illumina NextSeq 500 platform). Transcript abundance was assessed with the DESeq2, DESeq, and EdgeR Bioconductor R packages, and varying numbers of differentially expressed gene (DEG) transcripts were detected in the spermatozoa of each boar. Using RNA-Seq, we identified several genes associated with inflammation and apoptosis (FOS, NFATC3, ITGAL, EAF2 and ZDHHC14), spermatogenesis (FGF-14 and BAMBI), autophagy (RAB33B), protein phosphorylation (PTPRU and PTPN2) and energy metabolism (ND6 and ACADM) that were predominantly up-regulated in poor freezability ejaculates. Quantitative reverse transcription polymerase chain reaction (RT-qPCR) validated the transcript expression levels detected by RNA-Seq and thus confirmed the reliability of this technique. Subsequent validation with western blotting showed that the expression of three proteins was in accordance with the transcript abundance. Overall, we demonstrated that the up-regulation of the DEG transcripts in spermatozoa was associated with poor semen freezability. We suggest that spermatozoa transcriptome profiling provides a foundation to further elucidate the relevance of sperm-related transcripts on cryo-survival. The sperm-related transcripts, namely FOS, NFATC3, EAF2, BAMBI, PTPRU, PTPN2, ND6 and ACADM, are potential markers for predicting the freezability of boar semen.

Microarray of SNPs for diverse applications in commercial pig breeding.

Modern pig production needs new tools for fast, reliable, more effective breeding. In the present paper we present a chip containing 45 SNP (Single Nucleotide Polymorphisms) which enables the determining of 1 genetic disease (PSS-Porcine Stress Syndrome), 4 QTLs genes: PRKAG3, CAST, MC4R and ESR, which together with the remaining SNPs create a panel useful in marker-assisted selection and veterinary control. The SNPs were genotyped using the PCR-APEX (Arrayed Primer Extension) technique. Special attention is paid to evaluation of the 45 SNP chip as an alternative approach to parentage and identity control. Based on allele frequency estimations, for a sample of 88 individuals of commercial pig lines, the probabilities that a randomly chosen candidate parent would be excluded from paternity or maternity were estimated to be 99.9% when genotypes of both parents and a progeny were known, and 98% when the genotypes of only one parent and a piglet were available. The marker set presented here also reached a probability of identity in the order of 10(-16), which allows for unequivocal discrimination of animals or their products among billions of individuals. Further improvements for upcoming chip versions were also considered.

Total RNA quality in boar spermatozoa with different freezability.

In this study the quality of total RNA, isolated from fresh spermatozoa, was compared between boars with good and poor semen freezability (GSF and PSF, respectively). Semen from 3 boars with GSF exhibited significantly higher total motility, mitochondrial function, plasma membrane integrity and reduced lipid peroxidation compared with 3 boars with PSF after cryo- preservation. There were variations in the quality of RNA isolated from spermatozoa of boars with GSF and PSF. Boars with GSF exhibited mainly full-length, intact RNA, whereas substantial amounts of degraded RNA were detected in spermatozoa from boars with PSF. Further under- standing of the biological relevance of RNAs in sperm function is critical to improve the freezabil- ity of boar semen.

Single nucleotide polymorphism in the promoter region of the lactoferrin gene and its associations with milk performance traits in Polish Holstein-Friesian cows.

Bovine lactoferrin (LTF) is a multifunctional small glycoprotein found in milk acting mainly as a defense factor in the mammary gland. Many polymorphisms have been found in the bovine LTF gene but almost none were considered as genetic markers of production traits in dairy cattle. In this study, the promoter fragment of LTF gene containing mutation (G/C) in position +32 has been amplified by PCR followed by genotyping by the SSCP and RFLP method. 358 Polish Holstein-Friesian cows were screened, giving the following frequency of genotypes: 0.628, 0.313 and 0.059 for GG, GC and CC, respectively. GLM (General Linear Model) analysis was applied to evaluate the associations of lactoferrin with milk performance traits, including SCC - somatic cell count. It was found that CC cows show significantly higher (P < or = 0.01) protein content in milk in comparison with GG cows. The values of other milk performance traits were also higher but at non-significant levels. SCC in milk was the lowest in CC cows, but also at a non-significant level.

Genome-wide association study for host response to bovine leukemia virus in Holstein cows.

The mechanisms of leukemogenesis induced by bovine leukemia virus (BLV) and the processes underlying the phenomenon of differential host response to BLV infection still remain poorly understood. The aim of the study was to screen the entire cattle genome to identify markers and candidate genes that might be involved in host response to bovine leukemia virus infection. A genome-wide association study was performed using Holstein cows naturally infected by BLV. A data set included 43 cows (BLV positive) and 30 cows (BLV negative) genotyped for 54,609 SNP markers (Illumina Bovine SNP50 BeadChip). The BLV status of cows was determined by serum ELISA, nested-PCR and hematological counts. Linear Regression Analysis with a False Discovery Rate and kinship matrix (computed on the autosomal SNPs) was calculated to find out which SNP markers significantly differentiate BLV-positive and BLV-negative cows. Nine markers reached genome-wide significance. The most significant SNPs were located on chromosomes 23 (rs41583098), 3 (rs109405425, rs110785500) and 8 (rs43564499) in close vicinity of a patatin-like phospholipase domain containing 1 (PNPLA1); adaptor-related protein complex 4, beta 1 subunit (AP4B1); tripartite motif-containing 45 (TRIM45) and cell division cycle associated 2 (CDCA2) genes, respectively. Furthermore, a list of 41 candidate genes was composed based on their proximity to significant markers (within a distance of ca. 1 Mb) and functional involvement in processes potentially underlying BLV-induced pathogenesis. In conclusion, it was demonstrated that host response to BLV infection involves nine sub-regions of the cattle genome (represented by 9 SNP markers), containing many genes which, based on the literature, could be involved to enzootic bovine leukemia progression. New group of promising candidate genes associated with the host response to BLV infection were identified and could therefore be a target for future studies. The functions of candidate genes surrounding significant SNP markers imply that there is no single regulatory process that is solely targeted by BLV infection, but rather the network of interrelated pathways is deregulated, leading to the disruption of the control of B-cell proliferation and programmed cell death.

Evaluation of reference genes for qRT-PCR gene expression studies in whole blood samples from healthy and leukemia-virus infected cattle.

Real-time quantitative reverse transcription PCR (qRT-PCR) is the method of choice to investigate the alterations in gene expression involved with BLV pathogenesis. However, the reliability of qRT-PCR data critically depends on proper normalization to a set of stably expressed reference genes. The aim of the study was to validate the expression stability of ten candidate reference genes in RNA isolated from whole blood cells of BLV-negative and BLV-positive cows with hematological abnormalities. The rankings of the candidate genes according to their expression stability were calculated using BestKeeper, NormFinder and qBase(PLUS) software with implemented geNorm(PLUS) algorithm. The results showed that two genes are sufficient for normalization of qRT-PCR studies in whole blood RNA isolated from cows infected with BLV. According to geNorm, UCHL5 and RPLP0 were the best choice, but taking into account possible intergroup variation, NormFinder recommended RPLP0 and B2M as a most suitable pair. The overall ranking based on the geometric mean of the ranking numbers from each method separately showed UCHL5, RPLP0 and TBP as the most stable candidate reference genes. In addition, all three methods unanimously pointed at the commonly used ACTB and GAPDH as the least stable genes. These results further emphasize the need to accurately validate candidate reference genes before use in gene expression qRT-PCR studies.

Effect of new SNP within bovine prolactin gene enhancer region on expression in the pituitary gland.

A new single nucleotide polymorphism was revealed using PCR-SSCP and sequencing methods within the bovine prolactin distal promoter region described as a functional enhancer. The A-->G transition at position -1043 abolishes the recognition site for Hsp92II restriction endonuclease, allowing for PCR-RFLP genotyping. The application of real-time PCR revealed that the prolactin gene expression level in the pituitary was higher in cattle with the AA genotype than in those with the GG genotype. EMSA analysis, however, showed increased nuclear protein binding to the sequence variant with G, suggesting a possible inhibition event, in which the transcription factors Pit1, Oct1, and YY1 could be involved.

Associations between milk performance traits in Holstein cows and 16 candidate SNPs identified by arrayed primer extension (APEX) microarray.

An oligonucleotide microarray-which allows for parallel genotyping of many SNPs in genes involved in cow milk protein biosynthesis-was used to identify which of the 16 candidate SNPs are associated with milk performance traits in Holstein cows. Four hundred cows were genotyped by the developed and validated microarray. Significant associations were found between four single SNPs, namely DGAT1 (acyloCoA:diacylglycerol acyltransferase), LTF (lactoferrin), CSN3 (kappa-casein), and GHR (growth hormone receptor) and with fat and protein yield and percentage. Many significant associations between combined genotypes (two SNPs) and milk performance traits were found. The associations between the combined genotypes DGAT1/LTF and DGAT1/LEPTIN analyzed traits are presented as examples. The microarray based on APEX (Arrayed Primer Extension) is a fast and reliable method for multiple SNP analysis of potential application in marker-assisted selection. After further development, the chip may prospectively be used for dairy cattle paternity analysis and evolutionary studies.