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.

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.

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.

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.

[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<or=0.01 within one of the traits and the frequency of homozygotes with a minor allele of at least 0.1. Among 10 possible permutations, statistically significant associations were found only for a combination of SNP markers in TNNT3 x SFRSI genes and, interestingly, for all combinations with SNP located within the calpastatin (CAST) gene commonly known as a gene influencing pork quality traits. This study also found that CAST allele A (which is beneficial for pork tenderness) is also favorable for growth rate. This effect is clearly increased with additive alleles C of myogenic factor MYF6 and troponin T3 (TNNT3) and is decreased with each allele Tfrom serine rich splicing factor 1 (SFRS1) gene. For meat content, the most favorable genotype of calpastatin gene was AC, the effect of which was generally increased with each C allele of MYF6 and TNNT3 and decreased with each allele T from SFRSl and PKLR (puruvate kinase) genes, respectively. The calpastatin AC genotype seemed to be beneficial for selection index, although in combination with the CC genotype of troponin T3, calpastatin genotypes AA and CC reached higher values. In the case of the combination of genotypes TNNT3 x SFRS1, the most preferable for all analyzed traits is the CC genotype of TNNT3, especially in combination with CTor CC genotypes in the SFRS1 gene. We conclude that searching for interaction effects between candidate SNPs (even of moderate influence) may lead to interesting and valuable findings enabled better understanding and applications of SNPs in pork yield and quality improvement programmes.

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.

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.

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.