Document Detail

Genome wide scan for somatic cell counts in holstein bulls.
Jump to Full Text
MedLine Citation:
PMID:  21645296     Owner:  NLM     Status:  PubMed-not-MEDLINE    
BACKGROUND: Mastitis is the most costly disease for dairy production, and control of the disease is often difficult, due to its multi-factorial nature. Susceptibility to mastitis is under partial genetic control and the industry uses indirect selection for decreased concentrations of somatic cells in milk to reduce mastitis.
METHODS: A genome-wide scan was performed to identify genomic regions associated with deregressed estimated breeding values (EBVs) for somatic cell counts (SCC) in Holstein bulls. In total 1183 proven bulls of the Italian of Holstein population, were genotyped with the BovineSNP50 BeadChip (Illumina, San Diego, CA) and a whole genome association analysis was performed using the R package GenABEL.
RESULTS: Two chromosomal regions showed association with SCC, a region on chromosome 14 with high significance (P < 5x10-6) and a region on chromosome 6 with moderate significance (P < 5x10-5).
CONCLUSIONS: Two regions with effects on SCC have been identified with good statistical support. A further study of these candidate regions will be performed to verify the results and identify the causal mutations.
Giulietta Minozzi; Ezequiel L Nicolazzi; Francesco Strozzi; Alessandra Stella; Riccardo Negrini; Paolo Ajmone-Marsan; John L Williams
Related Documents :
3162536 - Mild and severe muscular dystrophy associated with deletions in xp21 of the human x chr...
9731526 - Dysferlin, a novel skeletal muscle gene, is mutated in miyoshi myopathy and limb girdle...
9781046 - Localisation of the gene for a dominant congenital spinal muscular atrophy predominantl...
3010296 - Dna linkage analysis of x chromosome-linked chronic granulomatous disease.
7759096 - Characterization of the human gene for microfibril-associated glycoprotein (mfap2), ass...
8536986 - Comparative analysis of qtls affecting plant height and maturity across the poaceae, in...
Publication Detail:
Type:  Journal Article     Date:  2011-06-03
Journal Detail:
Title:  BMC proceedings     Volume:  5 Suppl 4     ISSN:  1753-6561     ISO Abbreviation:  BMC Proc     Publication Date:  2011  
Date Detail:
Created Date:  2011-06-07     Completed Date:  2011-07-14     Revised Date:  2011-07-28    
Medline Journal Info:
Nlm Unique ID:  101316936     Medline TA:  BMC Proc     Country:  England    
Other Details:
Languages:  eng     Pagination:  S17     Citation Subset:  -    
Parco Tecnologico Padano, Via Einstein, Polo Universitario, Lodi 26900, Italy.
Export Citation:
APA/MLA Format     Download EndNote     Download BibTex
MeSH Terms

From MEDLINE®/PubMed®, a database of the U.S. National Library of Medicine

Full Text
Journal Information
Journal ID (nlm-ta): BMC Proc
ISSN: 1753-6561
Publisher: BioMed Central
Article Information
Download PDF
Copyright ©2011 Minozzi et al; licensee BioMed Central Ltd.
collection publication date: Year: 2011
Electronic publication date: Day: 3 Month: 6 Year: 2011
Volume: 5 Issue: Suppl 4
First Page: S17 Last Page: S17
ID: 3108211
PubMed Id: 21645296
Publisher Id: 1753-6561-5-S4-S17
DOI: 10.1186/1753-6561-5-S4-S17

Genome wide scan for somatic cell counts in holstein bulls
Giulietta Minozzi1 Email:
Ezequiel L Nicolazzi2 Email:
Francesco Strozzi1 Email:
Alessandra Stella13 Email:
Riccardo Negrini2 Email:
Paolo Ajmone-Marsan2 Email:
John L Williams1 Email:
1Parco Tecnologico Padano, Via Einstein, Polo Universitario, Lodi 26900, Italy
2Istituto di Zootecnica, Facoltà di Agraria, Università Cattolica del Sacro Cuore, Piacenza, 29122 Italy
3IBBA, CNR, Lodi, 26900, Italy


Mastitis, an inflammation of the mammary gland caused by an infection with a range of bacteria, is the most costly disease for dairy production. Control of mastitis is difficult due to its multi-factorial nature. Susceptibility to mastitis is under partial genetic control and the industry uses selection on a correlated trait (somatic cells score in milk), to reduce mastitis incidence in the population. Over the last few years, several studies have identified genetic loci putatively associated with somatic cell counts or clinical mastitis [1,2]. The availability of the bovine genome sequence and high density genotyping panels of single nucleotide polymorphisms has allowed a considerable number of bulls worldwide to be genotyped for genomic evaluation and selection. Furthermore, this information can be used to perform association studies with high precision at genome-wide level. The work reported here used genotypic data from the genomic selection project to perform a genome-wide scan with the objective of identifying genomic regions associated to deregressed estimated breeding values (DR-EBVs) for somatic cell counts (SCC) in Holstein bulls.


The bulls chosen for the genome wide association study were selected from among the 3155 animals progeny tested in Italy with DNA samples available. All these bulls will be used by the Italian National breeders association of Holstein Frisian Cattle (ANAFI) to perform national genomic evaluations.

Selection criteria used for association studies were intended to obtain: i) bulls with high selection index reliability (PFT > 0.75%); and ii) as low relationships between animals in the dataset as possible by trying to keep as many families (father – son couples) as possible. Among the 3155 bulls with biological material available 2109 bulls had appropriate criteria to be included in the study, 1183 of which had been already genotyped with the Bovine 50K SNP chip (Illumina Inc, San Diego).

Phenotype: deregressed EBV for SCC

The EBV for SCC had a mean of 98.73 ± 5.3 for the 2109 bulls, and a mean 98.77 ± 6.3 in the cohort of 1183 animals included in the study. Furthermore, deregressed EBVs (DR-EBVs) had mean of 0 and a standard deviation of 5. The DR-EBVs and reliabilities for somatic cell counts were derived from a reduced animal model for single records on a single trait. The strategy used to estimate the deregressed estimates was a simplified version of the algorithm of Jamrozik et al. [3], appropriate to a single trait reduced animal model.

Statistical analysis

Genome-wide association analysis was performed with the GenABEL package in R using a three step GRAMMAR-CG approach, (Genome wide Association using Mixed Model and Regression - Genomic Control) [4,5]. Uncorrected p-values of P < 5 x 10-7 were accepted to represent very strong proof of genome-wide association, while p-values between 5 x 10-7 and 5 x 10-5 were considered as moderately significant associations.

Genotyping and quality control filters

A total of 1183 progeny tested bulls were genotyped with the BovineSNP50 BeadChip (Illumina, San Diego, CA). Genotype quality assurance was performed within the R statistical environment using the GenABEL package (“check.marker” function) [6]. SNPs were checked for marker call rate (>5%) and minor allele frequency (<5%): markers missing 5% of data and with MAF of less than 5% were removed. Genotyping efficiency of samples was also verified, thus, samples with more than 5% missing data were removed. Classical Multi Dimension Scaling (MDS) was used to explore population substructure and to verify the genetic homogeneity of the dataset prior to analysis.

Results and discussion
Quality control

Following quality control checks, 641 markers were excluded because of low call rate and 11404 markers were excluded because of low minor allele frequency. Furthermore, markers on the sex chromosomes were removed from the analysis. A total of 8 samples were removed because of low call rate and other 2 were eliminated because of high autosomal heterozygosity (FDR < 1%). Mean heterozygosity of the dataset after quality check was 0.33 ± 0.01, while the samples removed had heterozygosity higher than 0.63, indicating possible sample contamination. No samples were removed due to high IBS (Identity By State). Mean IBS was 0.70 ± 0.01, based on 2000 autosomal markers, while the threshold for IBS was set to > 0.95. No outliers were identified by Classical Multi Dimension Scaling (MDS).

After quality controls, the final dataset used in the following association analysis contained 1173 samples and 41209 Genome wide SNPs.

Association analysis

Two chromosomal regions showed associations with SCC, a region on chromosome 14 with high significance (P < 5x10-6) and a region on chromosome 6 with moderate significance (P < 5x10-5). These two chromosome regions should be further tested to confirm these associations and to potentially identify the causative variations that affect this trait.

A recent review of QTL reported on chromosome 14 [7] identified 10 QTL for disease traits as mastitis, seven of which were related to somatic cell score [8-11]. Interestingly the SNP identified in this study and located on chromosome 14, is actually within 1Mb from the QTL identified by Kaupe et al. [7] to be associated to somatic cell score, but significantly distant from all other chromosomal regions that harbor QTL for clinical mastitis [7].

Furthermore Nilsen et al, [2] characterized a region of chromosome 6 in which QTL for clinical mastitis had been identified, and found the Mucin 7 gene to be significantly associated with the trait. Mucin7 is located close to the casein cluster on chromosome 6. However, the SNP located on chromosome 6 obtained in this study is more than 2 Mb distant from the casein cluster, indicating that different genes could be involved. To confirm the results found in the current study, both SNP identified will be tested in a second independent set of animals.

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

GM and JLW wrote the manuscript. JLW and PAM conceived the study. ELN extracted the phenotypic and pedigree information and defined the selection criteria. RN coordinated the sample collection and DNA preparation. GM performed statistical analysis. AS and FS contributed to the data analysis. All authors revised the paper, and approved the final manuscript.


Gerald Jansen is acknowledged for estimating deregressed proofs .Stefano Biffani is acknowledged for valuable suggestions and collaboration. This work was carried out through the SelMol and ProZoo projects.

This article has been published as part of BMC Proceedings Volume 5 Supplement 4, 2011: Proceedings of the International Symposium on Animal Genomics for Animal Health (AGAH 2010). The full contents of the supplement are available online at

Klungland H,Sabry A,Heringstad B,Olsen HG,Gomez-Raya L,Våge DI,Olsaker I,Ødegård J,Klemetsdal G,Schulman N,Vilkki J,Ruane J,Aasland M,Rønningen K,Lien S,Quantitative trait loci affecting clinical mastitis and somatic cell count in dairy cattleMamm GenomeYear: 2001128374210.1007/s00335001-2081-311845286
Nilsen H,Olsen HG,Hayes B,Nome T,Sehested E,Svendsen M,Meuwissen TH,Lien S,Characterization of a QTL region affecting clinical mastitis and protein yield on BTA6Anim GenetYear: 2009407011210.1111/j.1365-2052.2009.01908.x19466933
Jamrozik J,Schaeffer LR,Jansen GB,Approximate accuracies of prediction from random regression modelsLivestYear: 200066859210.1016/S0301-6226(00)00158-5
Amin N,van Duijn CM,Aulchenko YS,A genomic background based method for association analysis in related individualsPLoS ONEYear: 2007212e127410.1371/journal.pone.000127418060068
Aulchenko YS,de Koning DJ,Haley C,Genomewide rapid association using mixed model and regression: a fast and simple method for genome-wide pedigree-based quantitative trait loci association analysisGeneticsYear: 20071775778510.1534/genetics.107.07561417660554
Aulchenko YS,Ripke S,Isaacs A,van Duijn CM,GenABEL: an R package for genome-wide association analysisBioinformaticsYear: 2007231294610.1093/bioinformatics/btm10817384015
Wibowo TA,Gaskins CT,Newberry RC,Thorgaard GH,Michal JJ,Jiang Z,Genome assembly anchored QTL map of bovine chromosome 14Int J Biol SciYear: 200844061410.3923/ijb.2008.406.42019043607
Ashwell MS,Da Y,VanRaden PM,Rexroad CE Jr,Miller RH,Detection of putative loci affecting conformational type traits in an elite population of United States Holsteins using microsatellite markersJ Dairy SciYear: 1998811120112510.3168/jds.S0022-0302(98)75674-79594401
Kaupe B,Brandt H,Prinzenberg EM,Erhardt G,Joint analysis of the influence of CYP11B1 and DGAT1 genetic variation on milk production, somatic cell score, conformation, reproduction, and productive lifespan in German Holstein cattleJ Anim SciYear: 200785112110.2527/jas.2005-75317179535
Rodriguez-Zas SL,Southey BR,Heyen DW,Lewin HA,Interval and composite interval mapping of somatic cell score, yield, and components of milk in dairy cattleJ Dairy SciYear: 20028530819110.3168/jds.S0022-0302(02)74395-612487475
Rupp R,Lagriffoul G,Astruc JM,Barillet F,Genetic parameters for milk somatic cell scores and relationships with production traits in French Lacaune dairy sheepJ Dairy SciYear: 2003861476148110.3168/jds.S0022-0302(03)73732-112741573

Article Categories:
  • Proceedings

Conference: International Symposium on Animal Genomics for Animal Health (AGAH 2010). Paris, France. 31 May - 2 June 2010.

Previous Document:  Genome-wide association studies and genetic architecture of common human diseases.
Next Document:  Across-line SNP association study for (innate) immune and behavioral traits in laying hens.