Anna Ruść

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.

Prevalence of complex vertebral malformation carriers among Polish Holstein-Friesian bulls

An increasing number of Holstein calf births exhibiting vertebral deformations has been detected in Denmark since 1999 by a program monitoring the incidence of genetic diseases. Pedigree analysis demonstrated that the affected calves originated from a family afflicted by an autosomally recessively inherited complex vertebral malformation (CVM) syndrome. To determine the actual carrier frequency of the CVM-determining mutation in a population of Polish Holstein-Friesian (=Polish Black-and-White) cattle, we examined 202 proven bulls (active in 2001–2005) used by 4 domestic artificial insemination companies and 403 unproven bulls (under evaluation for breeding value). Out of the 605 bulls examined, 150 T/G heterozygotes were diagnosed, including 118 that were sons of known CVM carriers. Identification of a gene polymorphism in a bovine solute carrier family 35 member 3, termedSLC35A3, was conducted with the use of a new PCR-SSCP method (polymerase chain reaction - single stranded conformation polymorphism), which - due to its ease of use and high reliability - can be applied in widespread screening programs aimed at reducing the incidence of the CVM defect.

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.

On the origin of mongrels: evolutionary history of free-breeding dogs in Eurasia

Although a large part of the global domestic dog population is free-ranging and free-breeding, knowledge of genetic diversity in these free-breeding dogs (FBDs) and their ancestry relations to pure-breed dogs is limited, and the indigenous status of FBDs in Asia is still uncertain. We analyse genome-wide SNP variability of FBDs across Eurasia, and show that they display weak genetic structure and are genetically distinct from pure-breed dogs rather than constituting an admixture of breeds. Our results suggest that modern European breeds originated locally from European FBDs. East Asian and Arctic breeds show closest affinity to East Asian FBDs, and they both represent the earliest branching lineages in the phylogeny of extant Eurasian dogs. Our biogeographic reconstruction of ancestral distributions indicates a gradual westward expansion of East Asian indigenous dogs to the Middle East and Europe through Central and West Asia, providing evidence for a major expansion that shaped the patterns of genetic differentiation in modern dogs. This expansion was probably secondary and could have led to the replacement of earlier resident populations in Western Eurasia. This could explain why earlier studies based on modern DNA suggest East Asia as the region of dog origin, while ancient DNA and archaeological data point to Western Eurasia.

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.

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.

SNiPORK – A Microarray of SNPs in Candidate Genes Potentially Associated with Pork Yield and Quality – Development and Validation in Commercial Breeds

SNiPORK is an oligonucleotide microarray based on the arrayed primer extension (APEX) technique, allowing genotyping of single nucleotide polymorphisms (SNPs) in genes of interest for pork yield and quality traits. APEX consists of a sequencing reaction primed by an oligonucleotide anchored with its 5′ end to a glass slide and terminating one nucleotide before the polymorphic site. Extension with one fluorescently labeled dideoxynucleotide complementary to the template reveals the polymorphism. Ninety SNPs were selected from those associated directly or potentially with pork traits. Of the 90 SNPs, 5 did not produce a positive signal. For 85 SNPs, 100% repeatiblity was proved by double genotyping of 13 randomly chosen boars. In addition, the accuracy of genotyping was verified in 2 sib-families by a Mendelian inheritance of 49–50 homozygous genotypes from sire to sons. Three genotype discrepancies were found (97% accuracy rate). All inaccurities were confirmed by an alternative method (sequencing and PCR-RFLP assays). Moreover, the exclusion power of the chip was evalueted by an SNP inheritance analysis of unrelated boars within each sib-family. In the validation step, 88 boars (13 Pietrain, 31 Landrace, 16 Large White, 8 Duroc, 7 Hampshire x Pietrain crosses, and 13 other hybrid lines) were screened to validate SNPs. Among the 85 selected SNPs, 12 were found to be monoallelic, the rest showing at least two genotypes for the entire population under study. The primary application of the SNiPORK chip is the simultaneous genotyping of dozens of SNPs to study gene interaction and consequently better understand the genetic background of pork yield and quality. The chip may prospectively be used for evolutionary studies, evaluation of genetic distances between wild and domestic pig breeds, traceability tests, as well as the starting point for developing a platform for identification and paternity analysis.

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.

Assessment of genomic inbreeding in Polish Konik horses

The aim of this study was to assess the inbreeding coefficient of Polish Konik horses based on runs of homozygosity (ROH). Ninety six horses kept in 6 herds located across Poland were genotyped with the use of EquineSNP60 BeadChip (Illumina). SNP markers with a Minor Allele Frequency lower than 0.01 and SNPs assigned to chromosome X or Y were excluded from the study. A total of 50 708 SNPs were included for statistical analysis (SVS software, Golden Helix). The analysis showed that the population is in genetic equilibrium, with He and Ho estimates both equal to 0.3086. Seven categories of Runs of Homozygozity (ROH) length were defined: >0.5, >1, >2, >4, >8, >16, >25 Mb. The genomic inbreeding coefficient derived from ROH (FROH) calculated for each ROH length ranged from 15.96% based on the shortest ROH (>0,5Mb) to 2.71% for the longest ROH (>25Mb). Among individual horses, the inbreeding coefficient ranged from 5.25% to 22.41% (for ROH >1Mb). Analysis of ROH in Polish Koniks allows for more effective management of their inbreeding in the future.

Detection of Brachyspina carriers within Polish Holstein-Friesian bulls.

The aim of this paper was to verify the hypothesis whether carriers of genetic defect Brachyspina occur in the Polish Holstein-Friesian Cattle. PCR method was used to screen 78 Polish Holstein-Friesian bulls. Eight bulls were identified as heterozygotes for 3,3 kb deletion in the FANCI gene - the mutation causing Brachyspina defect. All carriers were sons of 3 sires: Cleitus Jabot, Sandy-Valley Bolton ET and Coyne-Farms Dorcy ET which were descendants of the US sire Sweet Haven Tradition (HOUSAM 1682485). Systematic screening of young bulls having in the pedigree Barchyspina carrier is necessary to prevent spreading of the recessive mutation in the dairy cattle population in Poland.