The somatic cell score (SCS) is implemented in routine sire evaluations in many countries as an indicator trait for udder health. Somatic cell score is highly correlated with clinical mastitis, and in the German Holstein population quantitative trait loci (QTL) for SCS have been repeatedly mapped onBos taurus autosome 18 (BTA18). In the present study, we report a refined analysis of previously detected QTL regions on BTA18 with the aim of identifying marker and marker haplotypes in linkage disequilibrium with SCS. A combined linkage and linkage disequilibrium approach was implemented, and association analyses of marker genotypes and maternally inherited two-marker-haplotypes were conducted to identify marker and haplotypes in linkage disequilibrium with a locus affecting SCS in the German Holstein population.
We detected a genome-wide significant QTL within marker interval 9 (HAMP_c.366+109G>A -BMS833) in the middle to telomeric region on BTA18 and a second putative QTL in marker interval 12-13 (BB710 -PVRL2_c.392G>A). Association analyses with genotypes of markers flanking the most likely QTL positions revealed the microsatellite markerBMS833 (interval 9) to be associated with a locus affecting SCS within the families investigated. A further analysis of maternally inherited two-marker haplotypes and effects of maternally inherited two-marker-interval gametes indicated haplotype249-G in marker interval 12-13 (BB710 -PVRL2_c.392G>A) to be associated with SCS in the German Holstein population.
Our results confirmed previous QTL mapping results for SCS and support the hypothesis that more than one locus presumably affects udder health in the middle to telomeric region of BTA18. However, a subsequent investigation of the reported QTL regions is necessary to verify the two-QTL hypothesis and confirm the association of two-marker-haplotype249-G in marker interval 12-13 (BB710 -PVRL2_c.392G>A) with SCS. For this purpose, higher marker density and multiple-trait and multiple-QTL models are required to narrow down the position of the causal mutation or mutations affecting SCS in German Holstein cattle.
Udder health, somatic cell score and subclinical and clinical mastitis remain major challenges for the economy of milk production in respect to milk production efficiency and animal health and welfare. Several studies have attempted to identify chromosomal regions, genes and polymorphisms that influence udder health in order to improve breeding strategies. SCS has been used as an indicator of udder health and is implemented in routine sire evaluations in many countries . SCS has a low to medium heritability (h2 = 0.15; ) and a strong correlation to mastitis in the German Holstein population (rg = 0.84; ). However, selection on low SCS as well as on decreased mastitis incidence is hampered by three aspects: first the low heritability of SCS and liability to mastitis, second the difficulties in recording mastitis related data and third by potential population-wide antagonisms between milk production traits (milk, fat and protein yield) and udder health [1,2]. Recently, marker assisted selection (MAS) has been determined as a promising tool to improve current selection strategies based on phenotypic data . MAS implements genetic marker information of confirmed QTL regions to identify individuals with favourable genetic background concerning the trait of interest. Thus, the confirmation and fine mapping of known QTL regions and the estimation of QTL effects will advance the use of MAS.
For clinical mastitis (CM) or SCS, QTL have been detected on nearly all autosomes  and several studies repeatedly detected QTL for SCS or CM in the telomeric region of BTA18 [5–12]. In addition, Kühn et al.  have shown in a proof-of-principle approach that information of 5 markers located in the telomeric region of BTA18 indeed enabled successful MAS, which identified halfsib heifers prior to first calving that exhibited significant differences in SCS after parturition.
The intention of this study was to further analyse the middle to telomeric region of BTA18 with the aim of identifying markers and marker haplotypes in linkage disequilibrium (LD) with SCS in German Holstein cattle to improve MAS for udder health. Therefore, we increased the marker density in the telomeric region on BTA18 and selected four functional candidate genes within the QTL regions reported by [5,7,11]. Polymorphisms detected within these candidate genes were used as additional markers for fine mapping previously identified QTL regions and to analyse effects of candidate gene polymorphisms on SCS in the German Holstein. In the present study, we detected a genome-wide significant QTL in the middle to telomeric region on BTA18. Furthermore, we analysed effects of maternally inherited marker haplotypes and identified a two-marker-haplotype associated with SCS in German Holstein cattle.
Selection of candidate genes
Based on the positional information derived from the previously mentioned QTL regions identified by [5,7,11] and preliminary results from microarray experiments, an intensive literature research of genes bearing potential function in innate immune defence, immune defence, mammary gland development or udder morphology was performed. To enable the investigation of positional candidate genes within the telomeric region on BTA18, a comparative map for BTA18 andHomo sapiens autosome 19 (HSA19) was constructed to take advantage of the more advanced gene annotation of the human genome (Additional file1: marker table). Four candidate genes were selected:calmodulin 3 (phosphorylase kinase, delta) (CALM3),Hepcidin antimicrobial peptide (HAMP),cadherin 1, type 1, E-Cadherin (ephitelial) (CDH1), andpoliovirus receptor-related 2 (herpesvirus entry mediator B) (PVRL2). Calmodulin 3 is a ubiquitously expressed Ca2+-binding protein that is involved in many Ca2+ modulated signal pathways. It was selected due to its influence on smooth muscle contraction [14,15] and the potential effects on milk leaking and milk flow, which is a trait with a substantial genetic correlation to SCS (rg = 0.4; ).HAMP was selected as an immunological candidate gene because of its antifungal and antibacterial activity , the toll-like receptor-4 dependent induction by bacterial pathogens in myeloid cells  and to a lesser extent due to its function as a key regulator of iron metabolism . ForCDH1 andPVRL2 initial results from own microarray experiments indicated their differential expression in clinically unaffected heifers with different predisposition to udder infection. Additionally, both genes are involved in cell-cell junctions that have a strong impact on tissue development. E-Cadherin, a Ca2+-dependent cell-cell adhesion molecule, is involved in tissue and organ development as part of the cadherin-catenin-complex [19,20], and poliovirus receptor-related 2, a Ca2+-independent immunoglobulin-like cell adhesion molecule, is involved in the organisation of intercellular junctions as part of the nectin-afadin-complex [20–22]. Furthermore E-Cadherin can serve as a receptor for pathogens and is involved in their internalization [23,24].
Screening for polymorphisms
Polymorphisms within the candidate genes were detected by comparative sequencing of genomic DNA from heifers selected for QTL alleles associated with high or low SCS. The heifers originate from most likely QTL heterozygous sires selected from the German Holstein population based on marker information regarding a confirmed QTL for SCS on BTA18. The selection strategy for the heifers and their phenotypes are described by Kühn et al.  in detail. For sequencing, we selected three pairs of halfsibs, where one daughter inherited the SCS increasing paternal chromosomal region (q) and the other inherited the SCS decreasing paternal chromosomal region (Q). In addition, two animals were screened for variants which originated from a genetically divergent Charolais × Holstein F2 cross background.
For each gene exceptCDH1 the entire gene (large introns excluded), 800 to 1500 bp of the promoter region and up to 500 bp downstream of the transcripts were investigated. ForCDH1, only the genomic sequence spanning exon 13 to 15 was analysed. Primer information and the genomic position for each primer are given in additional file2: primer table sequencing. Sequencing was performed by amplification of genomic DNA and subsequent sequence analysis with the DYEnamic ET Terminator Cycle Sequencing reaction and the MegaBACE™1000 DNA Analysis System (GE Healthcare, Munich, Germany). For evaluation of polymorphisms BioEdit 126.96.36.199  was used. Essentially, polymerase chain reaction (PCR) primers were used for sequencing. Additional sequencing primers were used only for longer PCR fragments or PCR fragments that were difficult to sequence.
The pedigree material used for genotyping included a total of 1,054 animals originating from six paternal halfsib families. Some of the animals are a subset of the granddaughter designs previously described by [7,11]. Numbers of sons per grandsire ranged from 60 to 353, with an average family size of 175 sons. The German genetic evaluation center (VIT) in Verden, Germany provided additional pedigree information including non-genotyped ancestors of genotyped animals (7,627 animals).
The phenotype information for SCS was provided by the VIT as daughter yield deviations (DYD) for the first lactation (Table1). SCS is the log2 transformed somatic cell count (log2 (somatic cell count/100000) + 3). DYD for SCS were calculated based on a random regression test day model . The reliabilities associated with the DYD were expressed as the number of effective daughter contributions (EDC) as described by . DYD and EDC for genotyped animals were obtained from the official release of the April 2008 routine genetic evaluation.
Descriptive statistics for daughter yield deviations
Number of Phenotypes
Arithmetic mean and standard deviation as well as lowest and highest daughter yield deviation for SCS considering all animals are given. SCS is the log2 transformed somatic cell count (log2 (somatic cell count/100000) + 3).
The marker set included a total of 28 markers covering the telomeric region of BTA18 fromCDH1 toDIK4013. Six of the 28 markers were already genotyped within previous QTL mapping studies [5,7,11] and include an erythrocyte antigen marker. The other 22 markers were newly selected and genotyped. Fifteen new microsatellite and seven new single nucleotide polymorphism (SNP) markers were chosen based on the information of the putative QTL positions reported by [5,7,11]. The fifteen new microsatellite markers were selected using the bovine linkage map of the United States Department of Agriculture's Meat Animal Research Center (MARC USDA) . SNP markers were selected from all detected polymorphisms within candidate genes based on allele frequencies, position (coding- or non-coding region), effect of the SNP (synonymous or nonsynonymous) and whether they are in linkage disequilibrium to each other.
Microsatellite markers were genotyped by PCR or Multiplex-PCR with fluorescence labelled primers followed by a fragment length analysis using the MegaBACE™1000 DNA Analysis System and MegaBACE Fragment Profiler Version 1.2 software (GE Healthcare, Munich, Germany) (Additional file3: primer table genotyping). The genotyping methods used for detection of SNP were PCR-restriction fragment length polymorphism (RFLP) assays forCDH1_c.2102C>T [NCBI dbSNP: rs41862198],HAMP_c.366+109G>A [NCBI dbSNP: ss141026745],PVRL2_c.-1268G>C [NCBI dbSNP: ss141026747],PVRL2_c.392G>A [NCBI dbSNP: rs41884977],CALM3_c.3+1795C>T [NCBI dbSNP: ss141026788] and a multiplex pyrosequencing assay forHAMP_c.86+430G>A [NCBI dbSNP: ss141026740] andCALM3_c.3+1678C>T [NCBI dbSNP: ss141026787] (Table2). The enzymes used for detection of RFLPs were identified with the NEBcutter V2.0 webtool [28,29] and PCR primers were designed with primer analysis software Oligo 4.1 (National Bioscience Inc., Plymouth, MN, USA). For the RFLP assays, a PCR specific for each SNP was used to amplify genomic DNA, and PCR products were subsequently incubated for 8 h with SNP specific restriction enzymes. For visualization of the RFLP a 2.5% agarose gel was used.
Primer and enzymes used for genotyping of SNP
5'-CAT AGA CAA CCA GAA CAA AGA C
5'-TGG ACC TCT GGG GAG ACT G
(Fermentas Life Sciences)
5'-AGA CAC CCA CTT TCC CAT C
5'-AGC TCC ACA GTC TCT TCT C
(Fermentas Life Sciences)
5'-AAT GCC AGT CAA TCA CAG TCT C
5'-GGA TTC TAC ACC CGC TGC TC
(Fermentas Life Sciences)
5'-TTC CTC AAA CTG TCT TAT CTG G
5'-GTG TAG TTG CCC TCG TCC TC
(Fermentas Life Sciences)
5'-TTG AGA GAA AAC CAG CAG AC
5'-CCA GGC AGC AGT GTT AGA
(New England Biolabs Inc.)
5'-GAG CCC TCC CTG AGT GCT TC
5'-AGC GGC TGC CTG TTC TCC
5'-AGG ATG GCT GCA CAC
5'-AAA AGA TGG TGG GAG AGT AAT GG
5'-CCT CTG CAC TTG CCT GTA AGA CTT
5'-CCA AAT AGG TCA AAT AAC A
Summary of methods and primers used for genotyping SNP. Forward primers for pyrosequencing were biotyinlated at the 5'-end for purification of single strand DNA.
Primer design for pyrosequencing was performed with Pyrosequencing™ Assay Design Software (Biotage AB, Uppsala, Sweden). A SNP specific PCR was used for amplification of genomic DNA and the products of both PCR were merged and analysed with the PSQ™HS 96A pyrosequencing system (Biotage AB, Uppsala, Sweden) in a multiplex run.
The genetic linkage map was calculated based on a refined marker order using CRIMAP software . The marker order on BTA18 used for the calculation was evaluated in two steps. First, information from published linkage maps [31,32], RH-maps , and human and bovine sequence-assemblies  were merged and compared to own linkage-mapping results. Second, in a region including the markersBMON117,DIK4672,BMS833,DIK4232 andBB710, no unequivocal marker order was obtained. Thus, the marker order in this region was verified by RH-mapping using the 12000 rad whole-genome radiation hybrid panel  and RH-MAP3.0 software . For some marker groups, a recombination rate of zero was calculated. To avoid technical difficulties arising in the calculation of transmitting probabilities the marker spacing was set to small values greater than zero (Additional file1: marker table).
A combined linkage and linkage disequilibrium analysis (LALD) was performed using the software system TIGER . TIGER is a UNIX script linking several individual Fortran programmes to perform combined linkage and linkage disequilibrium analysis. Six steps are implemented in the script. First, allele frequencies and transmitting probabilities for each putative QTL position are calculated using BIGMAP . The putative QTL positions were considered as the midpoint of each marker interval, resulting in a total of 27 putative QTL positions. Second, the identical by descent (IBD) sub-matrices for each putative QTL position are computed based on the gene dropping procedure described by [39–41] and third, the IBD sub-matrices are tested for positive definiteness and inverted. The software program COBRA  computes a condensed gametic relationship matrix and its inverse at each putative QTL position for the calculation of gametic effects. Transmitting probabilities and IBD sub-matrices are used for the set up of the condensed gametic relationship matrix. Finally, an LALD analysis is performed analysing every putative QTL position with restricted maximum likelihood (REML) methods applied in ASReml . A detailed description of the QTL mapping procedure and the model applied in ASReml is given by .
Analyses were conducted using a likelihood ratio test, where the REML of the full model was compared with the REML of the model missing the QTL effect. Chromosome-wide and genome-wide significance thresholds were determined as restricted log likelihood ratio (RLRT) equivalents of logarithm of odds (LOD) units  where a LOD > 2 indicates chromosome-wide (RLRT = 9.2) and a LOD > 3 indicates genome-wide significance (RLRT = 13.8) . Confidence intervals were estimated using the LOD drop-off method described by .
To investigate the association of candidate gene polymorphisms and markers flanking interval 9 and interval 12-13 with SCS, a mixed model including a random polygenic effect and the fixed effect of marker genotypes was applied in ASReml:
wherey is a vector of phenotypic observations (DYD) for sires,μ is the overall mean,MGi is the fixed effect of the marker genotypei,aj is the random polygenic effect of animalj andeij is the random residual. The polygenic effect that accounts for the family structure of the population was estimated using an extended pedigree of non-genotyped ancestors of genotyped animals including a total of 7,627 animals. To account for multiple testing, a 5% experiment-wise significance threshold was obtained by Bonferroni correction of the nominal p-value assuming a 5% Type 1 error (pexp = 0.0057).
Analysis of maternally inherited two-marker haplotypes and two-marker-interval gametes
Due to the limited number of sires, paternally inherited chromosomes could have a strong impact on genotype effects estimated in their offspring. To exclude these effects, maternally inherited two-marker-intervals including flanking markers of most likely QTL position were investigated. The most probable linkage phases of genotyped sires were calculated to determine the maternally inherited haplotypes using BIGMAP. The TIGER software system was then applied to estimate maternally inherited gametic effects for the putative QTL positions in interval 9 and interval 12 based on IBD sub-matrices. Finally, the estimated gametic effects of maternally inherited two-marker-interval gametes were analysed and plotted using SAS software (SAS Institute, Cary, NC, USA).
To verify differences in estimated two-marker-interval gamete effects, we performed an association analysis including a fixed maternally inherited two-marker-haplotype effect independent of the IBD sub-matrices. For this purpose, the same model was used as for marker genotypes (1), except that the fixed genotype effect was replaced by the fixed maternally inherited two-marker-haplotype effect. Two-marker-haplotypes with the highest and lowest mean for IBD gametic effects were tested against the total of all other haplotypes within respective intervals. For interval 12-13 additionally the most frequent two-marker-haplotype was investigated.
Screening for polymorphisms
Sequence analyses were based on the provisional reference sequences obtained from NCBI . ForCALM3 [GenBank: NM_001046249; GeneID: 520277], the entire gene (10728 bp) was resequenced including ~1.3 kb upstream from the transcription start (assumedly promoter region). A total of eleven polymorphisms were detected (Additional file4: polymorphisms): ten SNP and one 12 bp deletion within the assumed promoter region.
ForHAMP [GenBank: NM_001114508; GeneID: 512301], in silico analyses revealed that this gene is not annotated in Btau4.0 but has been annotated in NCBI Build3.1. The sequence ofHAMP is still located on BTA18 in Btau4.0 [GenBank: NC_007316.3: 476114 bp-477570 bp]. The whole gene (2777 bp) was resequenced, and three SNP were detected (Additional file4: polymorphisms), one SNP within intron 1 and two SNP within the first 110 nucleotides downstream ofHAMP transcript.
ForCDH1, only the coding sequence was provided as provisional mRNA reference sequence [GenBank: NM_001002763; GeneID: 282637]. Due to discrepancies between the mRNA reference sequence and the genome assembly sequence, only the genomic reference sequence was used for sequence comparisons. The genomic sequence including exon 13 to 15 (1673 bp) was investigated and six SNP were detected (Additional file4: polymorphisms) two of them were located in exon 13 and cause an amino acid substitution.
ForPVRL2 [GenBank: NM_001075210.1; GeneID: 505580], the whole gene, excluding 13181 bp of intron 2 and 670 bp of intron 4, was resequenced (11766 bp). A total of 17 SNP and two length polymorphisms were identified (Additional file4: polymorphisms). The examination of the promoter region (1 kb upstream of transcription start) ofPVRL2 revealed a gap of 367 bp with missing sequence information in the genome assembly Btau4.0. By a BLAST search  a whole genome shotgun sequence (WGS)-Trace [Trace Archive: 1632871820] that overlaps the gap was identified within the NCBI Trace Archives . This WGS-Trace was incorporated in the reference sequence and sequence information was confirmed by resequencing. Coordinates of polymorphisms upstream of the translation initiation codon ATG are indicated according to the updated sequence [GenBank: FJ829796].
Markers and Map
The genetic marker map covered 50 cM of the telomeric region of BTA18 (Additional file1: marker table). Marker intervals ranged from 0.05 cM to 7.9 cM with an average marker interval of 1.85 cM. The number of alleles for each marker ranged from two for SNP to 34 for the erythrocyte antigen marker. The marker order is in good agreement with previously published linkage- and RH-maps [31–33] and no discrepancies to the bovine sequence assembly Btau4.0  were observed (Additional file1: marker table).
The restricted Log Likelihood Ratio profile of the combined LALD analysis is shown in Figure1. Two peaks exceeding the genome-wide significance level can be observed. The maximum of the first peak is located in marker interval 9 at position 71.775 cM withHAMP_c.366+109G>A andBMS833 as flanking markers. The maximum of the second peak is located in marker interval 12 at position 77.6 cM with the flanking markersBB710 andPVRL2_c.-1268G>C. The LOD drop-off method was used to estimate confidence intervals for each of the peaks. An approximate 96% confidence interval included marker intervals 3 to 9 for the first peak and the marker intervals 2 to 15 for the second peak.
An association analysis between the candidate gene polymorphisms and SCS was performed to evaluate the influence of the selected candidate genes on the variation in SCS in German Holstein cattle. The markersBMS833 andBB710 were also included in the association analysis, as they are flanking markers of the maxima observed in the LALD test statistic. The association analyses (Table3) indicated a significant effect of theBMS833 genotype (p = 0.004) on SCS within the 5% experiment-wise significance threshold (pexp = 0.0057).PVRL2_c.392G>A (p = 0.017) andCALM3_c.3+1678C>T (p = 0.055) approached nominal significance, but were not significant at the 5% experiment-wise significance level. For the genotypes of SNP markersCDH1_c.2102C>T,HAMP_c.366+109G>A,HAMP_c.86+430G>A,PVRL2_c.-1268G>C andCALM3_c.3+1795C>T no significant effects on SCS were observed.
Association analyses of marker genotypes with SCS in six German Holstein halfsib families
Position in cM
Summary of results from association analyses for marker genotypes based on a mixed model (1) including a fixed genotype effect and a random polygenic effect. The experiment-wise 5% significance level is pexp = 0.0057.
Analysis of maternally inherited two-marker haplotypes and two-marker-interval gametes
To exclude any specific sire gamete effects that occur in a paternal halfsib design, maternally inherited gametic effects for two-marker-interval gametes were analysed and plotted for the putative QTL positions in interval 9 and interval 12-13. In interval 9 (HAMP_c.366+109G>A -BMS833), six maternally inherited two-marker-allele combinations occurred. Two-marker-interval gametes includingBMS833 alleles115 and119 were excluded, because they both occurred only once (Table4; Figure2A). Gametes carrying theBMS833 allele117 in interval 9 have a mean estimated maternally inherited gametic effect of 0.022 (± 0.003) and 0.028 (± 0.003), respectively, whereas gametes carrying allele113 of markerBMS833 have mean effects of - 0.0006 (± 0.0014) and - 0.0059 (± 0.0022), respectively. Investigating the gametic effect of each flanking marker of interval 9 alone (Figure2B), revealed that the main differences in gametic effects in the two-marker-interval results from the discrimination by microsatellite marker alleles. The difference in the mean effects for two-marker-interval gametes carrying the two alleles ofHAMP_c.366+109G>A is 0.009, whereas the difference between gametes carrying alleles113 and117 of markerBMS833 is 0.0268.
Estimated effects on SCS for maternally inherited two-marker-gametes in interval 9 and interval 12-13
Number of Animals
Mean of Estimates
Standard Error Mean
A - 113
A - 117
G - 113
G - 117
12 - 13
249 - G
12 - 13
249 - A
12 - 13
251 - G
12 - 13
251 - A
12 - 13
253 - G
12 - 13
253 - A
12 - 13
255 - G
12 - 13
255 - A
12 - 13
257 - G
12 - 13
257 - A
Summary of means for gametic effects estimated for maternally inherited two-marker-interval gametes considering IBD coefficients. Haplotype frequencies were determined by direct counting.
For interval 12-13 (BB710 -PVRL2_c.392G>A) we selectedPVRL2_c.392G>A (interval 13) as a flanking marker. Both polymorphisms,PVRL2_c.-1268G>C (interval 12) andPVRL2_c.392G>A (interval 13) are located in close vicinity within thePVRL2 gene and showed a high linkage disequilibrium (r2 = 0.68). In interval 12-13 (BB710 -PVRL2_c.392G>A), two-marker-interval gametes with twelve different two-marker-allele combinations were observed. Two-marker-interval gametes carrying the twoBB710 alleles243 and245 were excluded, because they both only occurred once. In addition, the frequencies for allele combinations253-G and257-G were smaller than 1% and therefore gametes carrying these two allele combinations are not included in Figure3 (Table4). Analogous to the situation in interval 9, the SNP alleles ofPVRL2_c.392G>A seem to have only a small influence on discriminating estimated two-marker-interval gamete effects (Figure3B), because the difference in the mean effects discriminated by the alleles is 0.0078. For microsatellite marker alleles (Figure3A), two-marker-interval gametes carrying the allele249 have mean effects of 0.0378 (± 0.0018) and 0.0176 (± 0.0009) and gametes carrying allele253 have mean effects of - 0.0150 (± 0.0021) and - 0.0253 (± 0.0012). The biggest differences in the mean estimated gametic effects was observed for the two two-marker-interval gametes249-G (0.0378 (± 0.0018)) and253-A (-0.0253 (± 0.0012)). Thus, the difference in the mean maternally inherited gametic effects equals 0.0631, which is equivalent to 0.16 phenotypic standard deviations.
To validate the differences observed between maternally inherited two-marker-interval gametes, we performed a direct association analysis for maternally inherited two-marker-haplotypes without considering IBD coefficients in interval 9 (HAMP_c.366+109G>A -BMS833) and interval 12-13 (BB710 -PVRL2_c.392G>A) (Table5). Only the two-marker-haplotypes with the highest and lowest mean for IBD gametic effects were tested against the total of all other haplotypes. For interval 12-13, we also tested255-A, because it was the most frequent haplotype (frequency at 40%) in the data set. In interval 9, no significant association of the maternally inherited haplotypesG-113 andG-117 were observed, whereas haplotype249-G (p = 0.027) in interval 12-13 seems to be associated with SCS in the German Holstein population.
Association analyses of maternally inherited two-marker-haplotypes with SCS
G - 113
G - 117
12 - 13
249 - G
12 - 13
253 - A
12 - 13
255 - A
Summary of results from an association analyses for maternally inherited two-marker-haplotypes based on a model (1) including a fixed two-marker-haplotype effect and a random polygenic effect. Individual two-marker-haplotypes were compared to all other maternally inherited haplotypes in the data set. Haplotype frequencies were determined by direct counting.
In an approach to identify marker and marker haplotypes affecting SCS in the German Holstein population, we analysed association of marker genotypes and maternally inherited two-marker-intervals flanking putative QTL positions detected in an LALD analysis for SCS on BTA18, and showed that the two-marker-haplotype249-G of interval 12-13 (BB710 -PVRL2_c.392G>A) is in LD with SCS in the German Holstein population.
Initially, we detected a genome-wide significant QTL for SCS on BTA18 (Figure1). The maxima of the QTL test statistic were in interval 9 (HAMP_c.366+109G>A -BMS833) and in interval 12 (BB710 -PVRL2_c.-1268G>C). In the same region telomeric on BTA18, several studies reported QTL for SCS. Kühn et al. , Brink  and Xu et al.  all identified a QTL for SCS in the German Holstein population at the telomeric end of BTA18 near markerTGLA227. Schulman et al.  reported a QTL for SCS as well as for mastitis at the telomeric end of BTA18 in Finnish Ayrshire cattle, and Ashwell et al.  localized a QTL for SCS in the middle to telomeric region on BTA18 at marker BM2078 in US Holstein, whereas Lund et al.  reported a QTL for SCS in Finnish Ayrshire, Swedish Red and White, and Danish Red with the maximum of the test statistic in the middle of BTA18 between the markers ILSTS002 and BMS2639. Interestingly Lund et al.  as well as Xu et al.  indicated that there might be more than one QTL for SCS on BTA18, but in both studies no significant evidence could be provided for a second QTL. Due to the lower marker density in previous studies and different approaches to detect QTL, a direct comparison between our results and previously mentioned studies is impeded. Nevertheless, the confidence intervals of the maxima observed in our LALD analysis did not include markerTGLA227 at the telomeric end of BTA18, and the QTL position reported by Lund et al.  as well as the assumption of a second QTL by Xu et al.  further in the middle of BTA18, indicated that we identified a second QTL for SCS in German Holstein cattle and possibly discovered a third QTL on BTA18 in our studies. Further indications for a third QTL on BTA18 in German Holstein cattle arise from association and haplotype analyses but no formal proof for a third QTL is given in our studies.
Subsequent to QTL mapping, we performed an association analysis of marker genotypes to verify the association of candidate gene polymorphisms and flanking markers of the most likely QTL positions observed in LALD analysis with SCS. Our results indicated that the microsatellite markerBMS833 (p = 0.004) is associated with SCS andPVRL2_c.392G>A (p = 0.017) andCALM3_c.3+1678C>T (p = 0.055) showed a respective tendency of association. Corresponding to the results of our LALD analysis, these results confirmed the position of one QTL for SCS in interval 9, asBMS833 is one of the flanking markers of this interval.PVRL2_c.392G>A andCALM3_c.3+1678C>T are flanking markers of the intervals 13 (PVRL2_c.-1268G>C -PVRL2_c.392G>A), 14 (PVRL2_c.392G>A - DIK3014) and intervals 16 (DIK4234 -CALM3_c.3+1678C>T) and 17 (CALM3_c.3+1678C>T -CALM3_c.3+1795C>T). The weak association observed for these markers is also in accordance to our LALD analysis, because interval 13 and 14 are within the genome-wide significance threshold in LALD analysis. Additionally, the second maximum observed in our LALD analysis in interval 12 is in an immediately adjacent interval and the putative QTL position in interval 12 is approximately 1 MB upstream ofPVRL2_c.392G>A. ForCALM3_c.3+1678C>T, only interval 16 (DIK4234 -CALM3_c.3+1678C>T) is within the chromosome-wide significance level in our LALD test statistic confirming the weaker association with SCS observed. To further test the results of our LALD analysis and the association ofBMS833 genotypes with SCS on a population wide level, we investigated maternally inherited two-marker-intervals for interval 9 (HAMP_c.366+109G>A -BMS833) as well as for interval 12-13 (BB710 -PVRL2_c.392G>A). First we analysed effects of maternally inherited two-marker-interval gametes estimated based on IBD coefficients and second we performed a direct association analysis for maternally inherited two-marker-haplotypes of interval 9 and interval 12-13 without considering IBD coefficients.
For interval 9 (HAMP_c.366+109G>A -BMS833), maternally inherited two-marker-interval gametes showed differences in estimated effects for SCS (Figure2). The variance of estimated gametic effects was higher for two-marker-interval gametes in interval 9 compared to two-marker-interval gametes in interval 12-13 and single marker analyses revealed that the microsatellite marker alleles are the main force in discriminating the effects of maternally inherited two-marker-interval gametes (Figure2, Figure3). Association analysis for maternally inherited two-marker-haplotypes of interval 9 showed that none of the maternally inherited two-marker-haplotypes in interval 9 (HAMP_c.366+109G>A -BMS833) are in LD with SCS in the German Holstein population. However, association analysis ofBMS833 genotypes showed an association with SCS in our half-sib design indicating that the association is due to linkage but not linkage disequilibrium of the BMS833 locus with the causal mutation affecting SCS.
For interval 12-13 (BB710 -PVRL2_c.392G>A), the biggest difference in gametic effects estimated for maternally inherited two-marker-interval gametes was observed between249-G (0.0378 (± 0.0018)) and253-A (- 0.0253 (± 0.0012)), where a positive mean of estimates indicates an unfavourable effect on SCS (high number of cells) and a negative mean of estimates indicates a favourable effect on SCS (low number of cells). Association analyses for the maternally inherited two-marker-haplotypes of interval 12-13 showed that249-G (p = 0.027) is associated with SCS at the nominal 5% significance level in the German Holstein population. The weak association ofPVRL2_c.392G>A genotype with SCS within the families we investigated and the association of maternally inherited two-marker-haplotype249-G with SCS indicates thatPVRL2_c.392G>A is not the causal mutation affecting SCS in German Holstein cattle, but the causal mutation has to be located near or within marker interval 12-13 (BB710 -PVRL2_c.392G>A). Combining results obtained for markersBMS833 andPVRL2_c.392G>A it still remains unclear whether there are two mutations (two QTL), one located near markerBMS833 (interval 9) and one near or within interval 12-13, or only one mutation, presumably in interval 12-13, affecting SCS in the middle part of BTA18.
The candidate gene polymorphisms we investigated were not the causal mutations affecting SCS in German Holstein cattle. However, the results of association analyses for single marker genotypes and maternally inherited two-marker-haplotypes indicated thatHAMP andPVRL2 were selected within the region harbouring at least one QTL for SCS. ParticularlyPVRL2 still remains interesting. On the one hand, the association of thePVRL2_c.392G>A genotype and that of the249-G haplotype suggests that another polymorphism withinPVRL2 is the causal mutation affecting SCS in the German Holstein population. On the other hand,PVRL2 was selected as a gene with possible impact on mammary gland development or udder morphology and several studies have detected QTL for udder conformation on BTA18 [51,52]. Therefore, it is also possible thatPVRL2 does not directly affect SCS but affects udder conformation traits like udder depth or fore udder attachment that are correlated with SCS .
In summary, our results suggest that the chromosomal region including interval 9 (HAMP_c.366+109G>A -BMS833) and interval 12-13 (BB710 -PVRL2_c.392G>A), in the middle to telomeric region on BTA18 has a strong impact on SCS in the German Holstein population. The analyses of maternally inherited two-marker-interval gamete effects and the association of the maternally inherited two-marker-haplotype249-G of interval 12-13 (BB710 -PVRL2_c.392G>A) with SCS indicates that micosatellite markerBB710 could be a suitable candidate-marker for MAS, but association of microsatellite markerBB710 with SCS has to be verified. To confirm the association of the two-marker-haplotype249-G with SCS and approve the hypothesis of two QTL in this region a further investigation is necessary. Thus, a mapping of udder conformation traits including a multiple-trait and multiple-QTL model might be useful to verify the existence of two QTL, and whether they are both directly affecting SCS or one is affecting a correlated trait. Likewise, a higher marker density within this region has to be achieved and families segregating for differentBB710 alleles have to be identified. Hence, it might be useful to cover the region including interval 9 and interval 12-13 with equally distributed SNP to narrow down the position of the casual mutation or mutations affecting SCS in the German Holstein population by a further fine mapping approach.
The authors would like to thank M. Schwerin and A. Hartmann for information on unpublished microarry data and the colleagues within the FUGATO M.A.S.net project for fruitful discussions. The financial support of the German Federal Ministry of Education and Research (BMBF) (Projekt FUGATO M.A.S.net, FKZ 0313390A) and the Development Association for Biotechnology Research (FBF) e.V., Bonn, is gratefully acknowledged. In addition, we thank F. Reinhardt, Z. Liu and the Vereinigte Informationssysteme Verden w.V. (VIT) for assistance with the selection of animals and for providing data.
Research Unit Molecular Biology, Research Institute for the Biology of Farm Animals
Research Unit Genetics and Biometry, Research Institute for the Biology of Farm Animals
Rupp R, Boichard D:Genetics of resistance to mastitis in dairy cattle.Vet Res 2003,34:671–688.View ArticlePubMed
Hinrichs D, Stamer E, Junge W, Kalm E:Genetic analyses of mastitis data using animal threshold models and genetic correlation with production traits.J Dairy Sci 2005,88:2260–2268.View ArticlePubMed
Kalm E:Development of cattle breeding strategies in Europe.Arch Tierz 2002,45:5–12.
Khatkar MS, Thomson PC, Tammen I, Raadsma HW:Quantitative trait loci mapping in dairy cattle: review and meta-analysis.Genet Sel Evol 2004,36:163–190. View ArticlePubMed
Brink M:Ein Beitrag zur Feinkartierung von QTL-Regionen für Eutergesundheit beim Rind.Schriftenreihe des Instituts für Tierzucht und Tierhaltung der Christian-Albrechts-Universität zu Kiel 2003,142:1–131.
Holmberg M, Andersson-Eklund L:Quantitative Trait Loci Affecting Health Traits in Swedish Dairy Cattle.J Dairy Sci 2004,87:2653–2659.View ArticlePubMed
Kühn C, Bennewitz J, Reinsch N, Xu N, Thomsen H, Looft C, Brockmann GA, Schwerin M, Weimann C, Hiendleder S, Erhardt G, Medjugorac I, Forster M, Brenig B, Reinhardt F, Reents R, Russ I, Averdunk G, Blumel J, Kalm E:Quantitative Trait Loci Mapping of Functional Traits in the German Holstein Cattle Population.J Dairy Sci 2003,86:360–368.View ArticlePubMed
Lund MS, Sahana G, Andersson-Eklund L, Hastings N, Fernandez A, Schulman N, Thomsen B, Viitala S, Williams JL, Sabry A, Viinalass H, Vilkki J:Joint Analysis of Quantitative Trait Loci for Clinical Mastitis and Somatic Cell Score on Five Chromosomes in Three Nordic Dairy Cattle Breeds.J Dairy Sci 2007,90:5282–5290.View ArticlePubMed
Schrooten C, Bovenhuis H, Coppieters W, Van Arendonk JAM:Whole Genome Scan to Detect Quantitative Trait Loci for Conformation and Functional Traits in Dairy Cattle.J Dairy Sci 2000,83:795–806.View ArticlePubMed
Schulman NF, Viitala SM, de Koning DJ, Virta J, Maki-Tanila A, Vilkki JH:Quantitative Trait Loci for Health Traits in Finnish Ayrshire Cattle.J Dairy Sci 2004,87:443–449.View ArticlePubMed
Xu N, Paul S, Bennewitz J, Reinsch N, Thaller G, Reinhardt F, Kühn C, Schwerin M, Erhardt G, Weimann C, Thomsen H, Mishra S, Kalm E:Confirmation of quantitative trait loci for somatic cell score on bovine chromosome 18 in the German Holstein.Arch Tierz 2006,49:111–119.
Baes C, Brand B, Mayer M, Kuhn C, Liu Z, Reinhardt F, Reinsch N:Refined positioning of a quantitative trait locus affecting somatic cell score on chromosome 18 in the German Holstein using linkage disequilibrium.J Dairy Sci 2009,92:4046–4054.View ArticlePubMed
Kühn C, Reinhardt F, Schwerin M:Marker assisted selection of heifers improved milk somatic cell count compared to selection on conventional pedigree breeding values.Arch Tierz 2008,51:23–32.
Pfitzer G:Signal Transduction in Smooth Muscle: Invited Review: Regulation of myosin phosphorylation in smooth muscle.J Appl Physiol 2001,91:497–503.PubMed
Walsh MP:Calmodulin and the regulation of smooth muscle contraction.Mol Cell Biochem 1994,135:21–41.View ArticlePubMed
Park CH, Valore EV, Waring AJ, Ganz T:Hepcidin, a Urinary Antimicrobial Peptide Synthesized in the Liver.J Biol Chem 2001,276:7806–7810.View ArticlePubMed
Peyssonnaux C, Zinkernagel AS, Datta V, Lauth X, Johnson RS, Nizet V:TLR4-dependent hepcidin expression by myeloid cells in response to bacterial pathogens.Blood 2006,107:3727–3732.View ArticlePubMed
Ganz T:Molecular Control of Iron Transport.J Am Soc Nephrol 2007,18:394–400.View ArticlePubMed
Ivanov DB, Philippova MP, Tkachuk VA:Structure and Functions of Classical Cadherins.Biochemistry (Mosc) 2001,66:1174–1186.View Article
Ebnet K:Organization of multiprotein complexes at cell-cell junctions.Histochem Cell Biol 2008,130:1–20.View ArticlePubMed
Takai Y, Ikeda W, Ogita H, Rikitake Y:The Immunoglobulin-Like Cell Adhesion Molecule Nectin and Its Associated Protein Afadin.Annu Rev Cell Dev Biol 2008,24:309–342.View ArticlePubMed
Ireton K:Entry of the bacterial pathogen Listeria monocytogenes into mammalian cells.Cell Microbiol 2007,9:1365–1375.View ArticlePubMed
Pizarro-Cerdá J, Cossart P:Bacterial Adhesion and Entry into Host Cells.Cell 2006,124:715–727.View ArticlePubMed
Hall TA:BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT.Nucleic Acids Symp Ser 1999,41:95–98.
Liu Z, Reinhardt F, Bunger A, Reents R:Derivation and Calculation of Approximate Reliabilities and Daughter Yield-Deviations of a Random Regression Test-Day Model for Genetic Evaluation of Dairy Cattle.J Dairy Sci 2004,87:1896–1907.View ArticlePubMed
Green P, Falls K, Crooks S:Documentation of CRI-MAP, Version 2.4 Washington University School of Medicine, St. Louis 1990.
Snelling W, Casas E, Stone R, Keele J, Harhay G, Bennett G, Smith T:Linkage mapping bovine EST-based SNP.BMC Genomics 2005,6:74. View ArticlePubMed
Ihara N, Takasuga A, Mizoshita K, Takeda H, Sugimoto M, Mizoguchi Y, Hirano T, Itoh T, Watanabe T, Reed KM, Snelling WM, Kappes SM, Beattie CW, Bennett GL, Sugimoto Y:A Comprehensive Genetic Map of the Cattle Genome Based on 3802 Microsatellites.Genome Res 2004,14:1987–1998.View ArticlePubMed
Itoh T, Watanabe T, Ihara N, Mariani P, Beattie CW, Sugimoto Y, Takasuga A:A comprehensive radiation hybrid map of the bovine genome comprising 5593 loci.Genomics 2005,85:413–424.View ArticlePubMed
Rexroad CE, Owens EK, Johnson JS, Womack JE:A 12 000 rad whole genome radiation hybrid panel for high resolution mapping in cattle: characterization of the centromeric end of chromosome 1.Anim Genet 2000,31:262–265.View ArticlePubMed
Boehnke M, Lunetta K, Hauser E, Lange K, Uro J, VanderStoep J:RHMAP: Statistical Package for Multipoint Radiation Hybrid Mapping Version 3.0 University of Michigan, Ann Arbor, Michigan 1996.
Baes C, Reinsch N:TIGER: A software system for fine-mapping quantitative trait loci.Arch Tierz 2008,51:402–412.
Reinsch N:A multiple-species, multiple-project database for genotypes at codominant loci.J Anim Breed Genet 1999,116:425–435.View Article
Meuwissen THE, Goddard ME:Fine mapping of quantitative trait loci using linkage disequilibria with closely linked marker loci.Genetics 2000,155:421–430.PubMed
Meuwissen THE, Goddard ME:Prediction of identity by descent probabilities from marker-haplotypes.Genet Sel Evol 2001,33:605–634. View ArticlePubMed
Meuwissen THE, Karlsen A, Lien S, Olsaker I, Goddard ME:Fine mapping of a quantitative trait locus for twinning rate using combined linkage and linkage disequilibrium mapping.Genetics 2002,161:373–379.PubMed
Baes C, Reinsch N:Computing the condensed conditional gametic QTL relationship matrix and its inverse.Arch Tierz 2007,50:294–308.
Gilmour AR, Gogel BJ, Cullis BR, Thompson R:ASReml User Guide Release 2.0 VSN International Ltd, Hemel Hempstead, HP1 1ES, UK 2006.
The Bovine Genome Sequencing and Analysis Consortium, Elsik CG, Tellam RL, Worley KC:The Genome Sequence of Taurine Cattle: A Window to Ruminant Biology and Evolution.Science 2009,324:522–528.View ArticlePubMed
Ashwell M, Rexroad CEJ, Miller R, Vanraden P, Da Y:Detection of loci affecting milk production and health traits in an elite US Holstein population using microsatellite markers.Anim Genet 1997,28:216–222.View Article
Kolbehdari D, Wang Z, Grant JR, Murdoch B, Prasad A, Xiu Z, Marques E, Stothard P, Moore SS:A Whole-Genome Scan to Map Quantitative Trait Loci for Conformation and Functional Traits in Canadian Holstein Bulls.J Dairy Sci 2008,91:2844–2856.View ArticlePubMed
Schnabel RD, Sonstegard TS, Taylor JF, Ashwell MS:Whole-genome scan to detect QTL for milk production, conformation, fertility and functional traits in two US Holstein families.Anim Genet 2005,36:408–416.View ArticlePubMed
Thomsen H, Reinsch N, Xu N, Looft C, Grupe S, Kuhn C, Brockmann GA, Schwerin M, Leyhe-Horn B, Hiendleder S, Erhardt G, Medjugorac I, Russ I, Forster M, Brenig B, Reinhardt F, Reents R, Blumel J, Averdunk G, Kalm E:A male bovine linkage map for the ADR granddaughter design.J Anim Breed Genet 2000,117:289–306.View Article
Liu KJ, Muse SV:PowerMarker: an integrated analysis environment for genetic marker analysis.Bioinformatics 2005,21:2128–2129.View ArticlePubMed
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.