Does a Little Difference Make a Big Difference? Bovine ß-Casein A1 and A2 Variants and Human Health-An Update.

For over 20 years, bovine beta-casein has been a subject of increasing scientific interest because its genetic A1 variant during gastrointestinal digestion releases opioid-like peptide ß-casomorphin-7 (ß-CM-7). Since ß-CM-7 is involved in the dysregulation of many physiological processes, there is a growing discussion of whether the consumption of the ß-casein A1 variant has an influence on human health. In the last decade, the number of papers dealing with this problem has substantially increased. The newest clinical studies on humans showed a negative effect of variant A1 on serum glutathione level, digestive well-being, cognitive performance score in children, and mood score in women. Scientific reports in this field can affect the policies of dairy cattle breeders and the milk industry, leading to the elimination of allele A1 in dairy cattle populations and promoting milk products based on milk from cows with the A2A2 genotype. More scientific proof, especially in well-designed clinical studies, is necessary to determine whether a little difference in the ß-casein amino acid sequence negatively affects the health of milk consumers.

Genomic Footprints of Recovery in the European Bison.

After extinction in the wild in the beginning of the 20th century, the European bison has been successfully recovered in 2 distinct genetic lines from only 12 and 7 captive founders. We here aimed at characterizing the levels of realized inbreeding in these 2 restored lines to provide empirical insights into the genomic footprints left by population recovery from a small number of founders. To that end, we genotyped 183 European bison born over the last 40 years with the Illumina BovineHD beadchip that contained 22 602 informative autosomal single-nucleotide polymorphisms after data filtering. We then identified homozygous-by-descent (HBD) segments and classified them into different age-related classes relying on a model-based approach. As expected, we observed that the strong and recent founder effect experienced by the 2 lines resulted in very high levels of recent inbreeding and in the presence of long HBD tracks (up to 120 Mb). These long HBD tracks were associated with ancestors living approximately from 4 to 32 generations in the past, suggesting that inbreeding accumulated over multiple generations after the bottleneck. The contribution to inbreeding of the most recent groups of ancestors was however found to be decreasing in both lines. In addition, comparison of Lowland individuals born at different time periods showed that the levels of inbreeding tended to stabilize, HBD segments being shorter in animals born more recently which indicates efficient control of inbreeding. Monitoring HBD segment lengths over generations may thus be viewed as a valuable genomic diagnostic tool for populations in conservation or recovery programs.

Effects of a novel missense polymorphism within the SIGLEC5 gene on fertility traits in Holstein-Friesian cattle.

The aim of the study was to find functional polymorphism within two exons of the SIGLEC5 (sialic acid-binding Ig-like lectin-5) gene and to examine its effects on the production and fertility traits of cows and bulls. Two hundred seventytwo Holstein-Friesian cows and 574 bulls were included in the study. Novel missense polymorphism (A > G) within exon 3 causing substitution of amino acid arginine by glutamate in position 260 of SIGLEC5 protein (R260Q) was identified by sequencing and digestion by restriction enzyme Msp I. Basic production and fertility traits of cows and estimated breeding values (EBV) of bulls were analysed. The study demonstrated a significant association of SIGLEC5 R260Q polymorphism with days open and calving interval in cows as well as with breeding value for calving interval in bulls. An opposite effect of SIGLEC5 alleles for production and fertility traits was observed: the allele G increased the breeding value for the protein yield, while the allele A increased the breeding value for the calving interval. The current study suggests the involvement of SIGLEC5 R260Q polymorphism in biological processes related to fertility traits. This finding can be applied as a biomarker for a genetic improvement programme in Holstein-Friesian cattle.

The assessment of inter-individual variation of whole-genome DNA sequence in 32 cows.

Despite the growing number of sequenced bovine genomes, the knowledge of the population-wide variation of sequences remains limited. In many studies, statistical methodology was not applied in order to relate findings in the sequenced samples to a population-wide level. Our goal was to assess the population-wide variation in DNA sequence based on whole-genome sequences of 32 Holstein-Friesian cows. The number of SNPs significantly varied across individuals. The number of identified SNPs increased with coverage, following a logarithmic curve. A total of 15,272,427 SNPs were identified, 99.16 % of them being bi-allelic. Missense SNPs were classified into three categories based on their genomic location: housekeeping genes, genes undergoing strong selection, and genes neutral to selection. The number of missense SNPs was significantly higher within genes neutral to selection than in the other two categories. The number of variants located within 3'UTR and 5'UTR regions was also significantly different across gene families. Moreover, the number of insertions and deletions differed significantly among cows varying between 261,712 and 330,103 insertions and from 271,398 to 343,649 deletions. Results not only demonstrate inter-individual variation in the number of SNPs and indels but also show that the number of missense SNPs differs across genes representing different functional backgrounds.

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.

Molecular Diversity of the Casein Gene Cluster in Bovidae: Insights from SNP Microarray Analysis.

The casein gene cluster spans 250 to 350 kb across mammalian species and is flanked by non-coding DNA with largely unknown functions. These regions likely harbor elements regulating the expression of the 4 casein genes. In Bovidae, this cluster is well studied in domestic cattle and to a lesser extent in zebu and water buffalo. This study used a cattle-specific SNP microarray to analyze 12 Bovidae taxa and estimate casein gene cluster variability across 5 bovid subfamilies. Genotyping identified 126 SNPs covering the entire casein gene cluster and 2 Mb of upstream and downstream regions. Dairy cattle, watusi, and zebu showed the highest polymorphism: 63.7-68.2% in the 5'-upstream region, 35.6-40.0% in the casein cluster, and 40.4-89.4% in the 3'-downstream region. Among wild bovids, only a 'semi-aquatic' lechwe revealed high polymorphism similar to cattle. Other species exhibited lower variability, ranging from 9.1-27.3% in the 5'-upstream, 8.9-20.0% in the casein, and 4.2-10.6% in the 3'-downstream regions. For the first time, genome variability data were obtained for impala, waterbuck, and lechwe. It appears that higher variability in cattle's casein gene cluster may relate to its intense expression. This study confirms the effectiveness of cattle-derived microarrays for genotyping Bovidae.

Milk from cows of different ß-casein genotypes as a source of ß-casomorphin-7.

The aim of this study was to quantify ß-casomorphin-7 in raw, hydrolyzed and processed milk in different stages of the cow lactation. The obtained results lead to the following conclusion: the highest amount of ß-casomorphin-7 released from the hydrolyzed and processed milk is related to the ß-casein A1 allele, irrespective of a lactation period. Some traces of ß-casomorphin-7 in milk from cows with the ß-casein A2 variant are probably a result of the acid hydrolysis of ß-casein during its digestion with pepsin. It has been shown that processing of raw milk at high temperatures affects, in a slight degree, the differences between ß-casomorphins-7 originating from different ß-casein genotypes. The obtained results suggest a possibility to provide a new nutritional factor for milk quality based on the content of ß-casomorphin-7 liberated in vivo from milk digested by a mixture of the gastrointestinal enzymes.

High-frequency marker haplotypes in the genomic selection of dairy cattle.

The aim of this study was to predict the genomic breeding value (DGV) of production, selected conformation and reproductive traits, and somatic cell score of dairy cattle in Poland using high-frequency marker haplotypes. The dataset consisted of phenotypic, genotypic, and pedigree data of 1216 Polish Holstein-Friesian bulls. The genotypic data consisted of 54,000 single-nucleotide polymorphisms (SNPs). The data were divided into two subsets: a test dataset (n = 1064) and a validation dataset (n = 152). Genotypic data were selected using three criteria: the percentage of missing genotypes, minor allele frequency, and linkage disequilibrium. The purpose of the data selection was to identify blocks of SNPs that were then used for the construction of haplotypes. Only haplotypes with a frequency higher than 25% were selected. DGV was predicted using four variants of a linear model with random haplotype effects and deregressed breeding values as the response variables. The accuracy of genomic prediction was checked by comparing DGVs with estimated breeding values (EBVs) using two methods: Pearson's correlations and the regression of EBV on DGV. The use of high-frequency haplotypes showed a tendency to underestimate DGVs. None of the models tested was clearly superior with regard to the traits studied. DGVs of production and conformation traits as well as somatic cell score (medium or high heritability traits) were more accurate than those estimated for fertility traits (low heritability traits).

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.

Identification of genome regions determining semen quality in Holstein-Friesian bulls using information theory.

Use of information theory can be an alternative statistical approach to detect genome regions and candidate genes that are associated with livestock traits. The aim of this study was to verify the validity of the SNPs effects on some semen quality variables of bulls using entropy analysis. Records from 288 Holstein-Friesian bulls from one AI station were included. The following semen quality variables were analyzed: CASA kinematic variables of sperm (total motility, average path velocity, straight line velocity, curvilinear velocity, amplitude of lateral head displacement, beat cross frequency, straightness, linearity), sperm membrane integrity (plazmolema, mitochondrial function), sperm ATP content. Molecular data included 48,192 SNPs. After filtering (call rate = 0.95 and MAF = 0.05), 34,794 SNPs were included in the entropy analysis. The entropy and conditional entropy were estimated for each SNP. Conditional entropy quantifies the remaining uncertainty about values of the variable with the knowledge of SNP. The most informative SNPs for each variable were determined. The computations were performed using the R statistical package. A majority of the loci had relatively small contributions. The most informative SNPs for all variables were mainly located on chromosomes: 3, 4, 5 and 16. The results from the study indicate that important genome regions and candidate genes that determine semen quality variables in bulls are located on a number of chromosomes. Some detected clusters of SNPs were located in RNA (U6 and 5S_rRNA) for all the variables for which analysis occurred. Associations between PARK2 as well GALNT13 genes and some semen characteristics were also detected.

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, termed SLC35A3, 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.

Polymorphism of bovine beta-casein and its potential effect on human health.

Proteins in bovine milk are a common source of bioactive peptides. The peptides are released by the digestion of caseins and whey proteins. In vitro the bioactive peptide beta-casomorphin 7 (BCM-7) is yielded by the successive gastrointestinal proteolytic digestion of bovine beta-casein variants A1 and B, but this was not seen in variant A2. In hydrolysed milk with variant A1 of beta-casein, BCM-7 level is 4-fold higher than in A2 milk. Variants A1 and A2 of beta-casein are common among many dairy cattle breeds. A1 is the most frequent in Holstein-Friesian (0.310-0.660), Ayrshire (0.432-0.720) and Red (0.710) cattle. In contrast, a high frequency of A2 is observed in Guernsey (0.880-0.970) and Jersey (0.490-0.721) cattle. BCM-7 may play a role in the aetiology of human diseases. Epidemiological evidence from New Zealand claims that consumption of beta-casein A1 is associated with higher national mortality rates from ischaemic heart disease. It seems that the populations that consume milk containing high levels of beta-casein A2 have a lower incidence of cardiovascular disease and type 1 diabetes. BCM-7 has also been suggested as a possible cause of sudden infant death syndrome. In addition, neurological disorders, such as autism and schizophrenia, seem to be associated with milk consumption and a higher level of BCM-7. Therefore, careful attention should be paid to that protein polymorphism, and deeper research is needed to verify the range and nature of its interactions with the human gastrointestinal tract and whole organism.

Effectiveness of a program aimed at the elimination of BLAD-carrier bulls from Polish Holstein-Friesian cattle.

The molecular basis of BLAD is the D128G mutation of the gene coding for the CD18 subunit of beta-2 integrin. This mutation is lethal, since homozygous (BL/BL) animals die before they reach sexual maturity. In the 1990s, BLAD was the most widespread genetic disease in HF cattle worldwide. The aim of the present study was to determine the frequency of BLAD carriers among 4645 young breeding bulls in Poland in 1995-2006. The frequency of carriers of the mutated allele showed a clear decreasing trend. The highest frequency (7.9%) was recorded while implementing the BLAD control program (1995-1997). Regular monitoring has enabled a great reduction of this threat to the tested population. Today only sporadic cases of BL/TL heterozygotes are reported (ca. 0.8% in 2004-2006).

Characterization of the CHORI-240 BAC clones containing the bovine CSN1S1, CSN2, STATH, CSN1S2 and CSN3 genes.

Nineteen BAC clones were identified by hybridization of the bovine genomic BAC library CHORI-240 with mixed CSN1S1- and CSN3-specific probes. Two of the clones were shown to contain the genes CSN1S1, CSN1S2, CSN2, STATH and CSN3, and five were proved to include the genes CSN2, STATH, CSN1S2 and CSN3. These data showed that the BAC contig was established for the whole casein cluster, including all known five genes.

Genome-wide association study for sperm membrane integrity in frozen-thawed semen of Holstein-Friesian bulls.

The aim of the study was to screen the entire bull genome to identify SNP markers and propose candidate genes potentially involved in the variation of sperm membrane integrity in Holstein-Friesian bulls. Two hundred eighty eight bulls kept in one AI center were included in the study. Each bull was genotyped for 54.001 Single Nucleotide Polymorpisms (SNP) by the Illumina BovineSNP50 BeadChip. Commercial straws of frozen-thawed semen were used for the evaluation of sperm plasma membrane integrity (SYBR-14/PI staining) and sperm mitochondrial function (JC1/PI staining). An additive model for Linear Regression Analysis was applied to estimate the effect of SNP marker for sperm membrane integrity (by the use of GoldenHelix SVS7 software). Five significant markers (encompassing 2,2 MB region located on chromosome 6) for SYBR-14/PI were found. Among them one marker-rs41570391 passed Bonferroni correction test. Within approximately 3 Mb genomic region including significant markers three candidate genes: SGMS2 (Sphingomyelin Synthase 2), TET2 (Methylcytosine dioxygenase 2) and GSTCD genes (Gluthatione S-transferase C terminal domain) were proposed as potentially involved in sperm membrane integrity in frozen-thawed semen of Holstein-Friesian bulls.

Novel Graphical Analyses of Runs of Homozygosity among Species and Livestock Breeds.

Runs of homozygosity (ROH), uninterrupted stretches of homozygous genotypes resulting from parents transmitting identical haplotypes to their offspring, have emerged as informative genome-wide estimates of autozygosity (inbreeding). We used genomic profiles based on 698 K single nucleotide polymorphisms (SNPs) from nine breeds of domestic cattle (Bos taurus) and the European bison (Bison bonasus) to investigate how ROH distributions can be compared within and among species. We focused on two length classes: 0.5-15 Mb to investigate ancient events and >15 Mb to address recent events (approximately three generations). For each length class, we chose a few chromosomes with a high number of ROH, calculated the percentage of times a SNP appeared in a ROH, and plotted the results. We selected areas with distinct patterns including regions where (1) all groups revealed an increase or decrease of ROH, (2) bison differed from cattle, (3) one cattle breed or groups of breeds differed (e.g., dairy versus meat cattle). Examination of these regions in the cattle genome showed genes potentially important for natural and human-induced selection, concerning, for example, meat and milk quality, metabolism, growth, and immune function. The comparative methodology presented here permits visual identification of regions of interest for selection, breeding programs, and conservation.

Diversifying Selection Between Pure-Breed and Free-Breeding Dogs Inferred from Genome-Wide SNP Analysis.

Domesticated species are often composed of distinct populations differing in the character and strength of artificial and natural selection pressures, providing a valuable model to study adaptation. In contrast to pure-breed dogs that constitute artificially maintained inbred lines, free-ranging dogs are typically free-breeding, i.e., unrestrained in mate choice. Many traits in free-breeding dogs (FBDs) may be under similar natural and sexual selection conditions to wild canids, while relaxation of sexual selection is expected in pure-breed dogs. We used a Bayesian approach with strict false-positive control criteria to identify FST-outlier SNPs between FBDs and either European or East Asian breeds, based on 167,989 autosomal SNPs. By identifying outlier SNPs located within coding genes, we found four candidate genes under diversifying selection shared by these two comparisons. Three of them are associated with the Hedgehog (HH) signaling pathway regulating vertebrate morphogenesis. A comparison between FBDs and East Asian breeds also revealed diversifying selection on the BBS6 gene, which was earlier shown to cause snout shortening and dental crowding via disrupted HH signaling. Our results suggest that relaxation of natural and sexual selection in pure-breed dogs as opposed to FBDs could have led to mild changes in regulation of the HH signaling pathway. HH inhibits adhesion and the migration of neural crest cells from the neural tube, and minor deficits of these cells during embryonic development have been proposed as the underlying cause of "domestication syndrome." This suggests that the process of breed formation involved the same genetic and developmental pathways as the process of domestication.

Nucleotide sequence polymorphism within exon 4 of the bovine prolactin gene and its associations with milk performance traits.

Prolactin plays an important regulatory function in mammary gland development, milk secretion, and expression of milk protein genes. Hence the PRL gene is a potential quantitative trait locus and genetic marker of production traits in dairy cattle. We analysed the sequence of the PRL gene to investigate whether mutations in this sequence might be responsible for quantitative variations in milk yield and composition. Using SSCP and direct sequencing, we detected six single-nucleotide polymorphisms within a 294-bp prolactin gene fragment involving exon 4. All detected mutations were silent with respect to the amino acid sequence of the protein. PCR-RFLP genotyping of SNP 8398 R (RsaI) was used to assess allele frequencies in 186 Black-and-White cows (0.113 and 0.887 for A and G, respectively) and in 138 Jersey cows (0.706 and 0.294 for A and G, respectively). Black-and-White cows with genotype AG showed the highest milk yield, while cows with genotype GG showed the highest fat content.

No incidence of DUMPS carriers in Polish dairy cattle.

DUMPS (Deficiency of Uridine Monophosphate Synthase) is a hereditary recessive disorder in Holstein cattle causing early embryo mortality during its implantation in the uterus. The only way to avoid the economic losses is early detection of DUMPS carriers. Because American Holstein semen has been intensively imported to Poland since 1970, there was a risk that DUMPS could have spread in Polish dairy cattle. In our study, 2209 dairy cattle of the Polish Holstein breed have been screened by the DNA test. The dominant group was young bulls entering the testing program (1171) and proven bulls (781). They represented all sires entering Polish breeding programs between 1999 and 2003. Also, 257 sire dams were included in the screening program. No DUMPS carrier has been found. Our results then indicate that the population of dairy cattle reared in Poland is free from DUMPS. Because of the economical significance of the DUMPS mutation and its recessive mode of inheritance, attention has to be paid to any case of a bull having in his origin any known DUMPS carrier. Such a bull should be tested and if positive eliminated from the active population. Also, young bulls (testing bulls) should be screened for DUMPS if in their progeny a high incidence of embryo mortality is observed and their genealogy cannot exclude their relatedness to any DUMPS carriers.

MilkProtChip--a microarray of SNPs in candidate genes associated with milk protein biosynthesis--development and validation.

MilkProtChip is an oligonucleotide microarray based on the arrayed primer extension (APEX) technique, allowing genotyping of single nucleotide polymorphisms (SNPs) in genes of interest for bovine milk protein biosynthesis. 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. The extension with one fluorescently labeled dideoxy nucleotide complementary to the template reveals the polymorphism. A total of 75 SNPs were selected among those associated directly or potentially with milk protein content. Among the 75 SNPs, 4 did not produce a positive signal. Most of the remaining SNPs produced a signal for both strands, except for 4 (one strand). In the validation step, 12 Polish Holstein bulls, 1 Polish Red bull, 1 bison (Bison bonasus), 11 Jersey cows and 25 Polish Holstein cows were screened to validate SNPs. Among the 71 selected SNPs--26 were found monoallelic, the rest showing at least two genotypes for the entire population under study. All the animals were earlier genotyped for 2-5 SNPs by PCR-RFLP and PCR sequencing and all showed complete concordance with APEX genotyping. APEX reactions showed relatively high signal frequencies: more than 0.9, 0.9-0.8 and below 0.8, for 65, 4 and 2 DNA samples, respectively. The primary application of the MilkProtChip is the simultaneous genotyping of dozens of SNPs to reveal and clarify the genetic background of milk protein biosynthesis. The chip may possibly be used for dairy cattle identification and paternity analysis, evolutionary studies, the evaluation of genetic distances between wild and domestic cattle breeds and the domestication history of bovine species.