SNP panel for evaluation of genetic variability and relationship in roe deer (Capreolus capreolus).

Blood samples from forty-six roe deer ( Capreolus capreolus) acquired during officially approved hunting in six hunting divisions throughout Poland were used to isolate the genomic DNA. All individuals were genotyped by MD_Bovine BeadChip (Illumina) for 46.750 Single Nucleotide Polymorphism (SNP) markers. SNPs of inappropriate clusters, with a marker call rate lower than 90% and with a minor allele frequency (MAF) lower than 0.01, located on sex chromosomes and mitochondrial DNA, were removed. Altogether, 21.033 SNP markers were included for further analysis. Observed and expected heterozygosity amounted to 0.098 and 0.119, respectively. Among 21.033 markers, a panel of 148 SNPs were selected for relationship analysis. They were unlinked and had a MAF higher than 0.2. This set of SNPs showed a probability of parentage exclusion of 1.29x10 -6 and 2.37x10 -19 for one, and two known parents, respectively. The probability of identity was estimated at 1.8x10 -40. The probabilities obtained in this study are sufficient for the monitoring and effective management of the genetic diversity of roe deer in Poland and is a cost-effective complementary tool for forensic applications.

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.

Genome-wide association study for semen volume and total number of sperm in Holstein-Friesian bulls.

In artificial insemination industry bulls producing high volume of semen with relatively high concentration of sperm are very desirable since they ensure stable production of commercial straws especially in case of top bulls. The aim of the study was to screen the entire bull genome to identify markers and candidate genes underlying semen volume (SV) and total number of sperm (TNS) in ejaculate produced by Holstein-Friesian bulls. Data on semen production were retrieved from records of AI center and included a population of 877 Holstein-Friesian bulls. Each bull was genotyped using the Illumina BovineSNP50 BeadChip. Genome-wide association analysis was performed with the use of GoldenHelix SVS7 software. An additive model for Linear Regression Analysis was used to estimate the effect of SNP marker for SV and TNS. After Bonferroni correction, 3 markers located on chromosome 22 reached the highest significance (rs41625599, rs41584616, rs42012507) for both traits. In the vicinity of these significant markers 3 genes are located (DCP1A, SFMBT1, TMEM110). Moreover, marker rs110109069 located on chromosome 25 was significantly associated with TNS and marker rs42438348 located on chromosome 10 has been found to be associated with SV. Some additional candidate genes were suggested to be potentially involved in analyzed traits (GALC, PRKCD, PHF7, TLR9, SPATA7). Identifying SNPs associated with the lower total number of sperm may be very useful for early recognition of a young sire as less suitable for effective semen production in artificial insemination centers.

Genome-wide association study for sperm concentration in Holstein-Friesian bulls.

The aim of the study was to screen the entire bull genome to identify markers and candidate genes underlying sperm concentration. The analysed data set originates from a population of 877 Polish Holstein-Friesian bulls. Based on sperm concentration value, two extreme groups of bulls were created: Low (L, n = 126) and High (H, n = 140). Each bull was genotyped using the Illumina BovineSNP50 BeadChip. Genome-wide association analysis was performed with the use of GoldenHelix SVS7 software. An additive model with a Cohran-Armitage test, Correlation/Trend adjusted by a Bonferroni test, was used to estimate the effect of SNP marker for sperm concentration. Thirteen markers reached genome-wide significance. The most significant SNPs were located on chromosome 3 (rs109154964 and rs108965556), 14 (rs41621145) and 18 (rs41615539), in the close vicinity of protein arginine methyltransferase 6 (PRMT6), Sel1 repeat containing 1 (SELRC1), triple QxxK/R motif containing (TRIQK) and zinc finger homeobox 3 (ZFHX3) genes, respectively. For three other candidate genes located close to significant markers (within a distance of ca 1 Mb), namely histone deacetylase 9 (HDAC9), an inhibitor of DNA binding 2 (ID2) and glutathione S-transferase theta 1 (GSTT1), their potential role in the production of male germ cells was confirmed in earlier studies. Six additional candidate genes (Vav3, GSTM1, CDK5, NOS3, PDP1 and GAL3ST1) were suspected of being significantly associated with sperm concentration or semen biochemistry. Our results indicate the genetic complexity of sperm concentration but also open the possibility for finding causal polymorphism useful in marker-assisted selection.

Genome-wide association study for poor sperm motility in Holstein-Friesian bulls.

The aim of the study was to screen the whole bull genome to identify markers and candidate genes underlying poor sperm motility. The analyzed data set originates from the Polish Holstein-Friesian bull population and consists of 41 Case and 279 Control bulls (selected from 1581 bulls). The most distinguishing trait of case group was very poor sperm motility (average 25.61%) when compared to control samples (average 72.95%). Each bull was genotyped using the Illumina BovineSNP50 BeadChip. Genome-wide association analysis was performed with the use of GoldenHelix SVS7 software. An additive model with a Cohran-Armitage test, Correlation/Trend adjusted by Bonferroni test were used to estimate the effect of Single Nucleotide Polymorphism (SNP) marker for poor sperm motility. Markers (n=34) reached genome-wide significance. The most significant SNP were located on chromosome 24 (rs110876480), 5 (rs110827324 and rs29011704), and 1 (rs110596818), in the close vicinity of melanocortin 4 receptor (MC4R), PDZ domain containing ring finger 4 (PDZRN4) and ethanolamine kinase 1 (ETNK1), olfactory receptor 5K3-like (LOC785875) genes, respectively. For five other candidate genes located close to significant markers (in distance of ca. 1 Mb), namely alkaline phosphatase, liver/bone/kidney (ALPL), tripartite motif containing 36 (TRIM36), 3-hydroxyisobutyrate dehygrogenase (HIBADH), kelch-like 1 (KLHL1), protein kinase C, beta (PRKCB), their potential role in sperm motility was confirmed in the earlier studies. Five additional candidate genes, cystic fibrosis transmembrane conductance regulator (CFTR), insulin-like growth factor 1 receptor (IGF1R), steroid-5-alpha-reductase, alpha polypeptide 2 (SRD5A2), cation channel, sperm associated 1 (CATSPER1) calpain 1 (mu/I) large subunit (CAPN1) were suggested to be significantly associated with sperm motility or semen biochemistry. Results of the present study indicate there is a genetic complexity of poor sperm motility but also indicate there might be a causal polymorphism useful in marker-assisted selection. Identifying genomic regions associated with poor sperm motility may be very important for early recognition of a young sire as unsuitable for effective semen production in artificial insemination centers.

Towards an integrated approach to study SNPs and expression of candidate genes associated with milk protein biosynthesis.

MilkProtChip is oligonucleotide microarray allowing bovine genotyping based on single nucleotide polymorphisms (SNPs) in genes influencing milk protein biosynthesis. A total of 71 SNPs in 42 genes were selected as associated with milk protein biosynthesis. Genotyping of about 300 animals of Polish Black-and-White cattle showed that SNPs in acyl-CoA: 1,2-diacylglycerol O-transferase (DGAT1), lactoferrin (LTF), casein kappa (CSN3) and growth hormone receptor (GHR) genes were associated with several milk performance traits. Analysis of correlations between SNPs and milk production traits showed that SNPs in single genes rarely affect the investigated traits. Only 4 of 42 investigated single SNPs had impact on milk production traits while 22 combinations of paired SNPs in these genes had impact. Positive effect SNP combinations in two genes can be a result of additive effect on these SNPs on the same traits or effect of genes interaction. The MilkBovExp chip representing 90 genes encoding transcription factors expressed in the bovine mammary gland and/or involved in mammary gland signaling pathways was designed for further investigation of impact of gene expression and/or its encoded products on milk traits performance.

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.