Identification of a two-marker-haplotype on Bos taurus autosome 18 associated with somatic cell score in German Holstein cattle
© Brand et al; licensee BioMed Central Ltd. 2009
Received: 15 June 2009
Accepted: 2 September 2009
Published: 2 September 2009
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 on Bos 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 marker BMS833 (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 haplotype 249-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-haplotype 249-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 and Homo sapiens autosome 19 (HSA19) was constructed to take advantage of the more advanced gene annotation of the human genome (Additional file 1: 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), and poliovirus 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 . For CDH1 and PVRL2 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 except CDH1 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. For CDH1, only the genomic sequence spanning exon 13 to 15 was analysed. Primer information and the genomic position for each primer are given in additional file 2: 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).
Descriptive statistics for daughter yield deviations
Number of Phenotypes
The marker set included a total of 28 markers covering the telomeric region of BTA18 from CDH1 to DIK4013. 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.
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
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 markers BMON117, DIK4672, BMS833, DIK4232 and BB710, 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 file 1: 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 .
where y is a vector of phenotypic observations (DYD) for sires, μ is the overall mean, MGi is the fixed effect of the marker genotype i, aj is the random polygenic effect of animal j and eij 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 . For CALM3 [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 file 4: polymorphisms): ten SNP and one 12 bp deletion within the assumed promoter region.
For HAMP [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 of HAMP 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 file 4: polymorphisms), one SNP within intron 1 and two SNP within the first 110 nucleotides downstream of HAMP transcript.
For CDH1, 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 file 4: polymorphisms) two of them were located in exon 13 and cause an amino acid substitution.
For PVRL2 [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 file 4: polymorphisms). The examination of the promoter region (1 kb upstream of transcription start) of PVRL2 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 file 1: 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 file 1: marker table).
Association analyses of marker genotypes with SCS in six German Holstein halfsib families
Position in cM
Analysis of maternally inherited two-marker haplotypes and two-marker-interval gametes
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
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
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-haplotype 249-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 (Figure 1). 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 marker TGLA227. 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 marker TGLA227 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 marker BMS833 (p = 0.004) is associated with SCS and PVRL2_c.392G>A (p = 0.017) and CALM3_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, as BMS833 is one of the flanking markers of this interval. PVRL2_c.392G>A and CALM3_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 of PVRL2_c.392G>A. For CALM3_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 of BMS833 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 (Figure 2). 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 (Figure 2, Figure 3). 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 of BMS833 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 between 249-G (0.0378 (± 0.0018)) and 253-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 that 249-G (p = 0.027) is associated with SCS at the nominal 5% significance level in the German Holstein population. The weak association of PVRL2_c.392G>A genotype with SCS within the families we investigated and the association of maternally inherited two-marker-haplotype 249-G with SCS indicates that PVRL2_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 markers BMS833 and PVRL2_c.392G>A it still remains unclear whether there are two mutations (two QTL), one located near marker BMS833 (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 that HAMP and PVRL2 were selected within the region harbouring at least one QTL for SCS. Particularly PVRL2 still remains interesting. On the one hand, the association of the PVRL2_c.392G>A genotype and that of the 249-G haplotype suggests that another polymorphism within PVRL2 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 that PVRL2 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-haplotype 249-G of interval 12-13 (BB710 - PVRL2_c.392G>A) with SCS indicates that micosatellite marker BB710 could be a suitable candidate-marker for MAS, but association of microsatellite marker BB710 with SCS has to be verified. To confirm the association of the two-marker-haplotype 249-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 different BB710 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.
- Rupp R, Boichard D: Genetics of resistance to mastitis in dairy cattle. Vet Res. 2003, 34: 671-688. 10.1051/vetres:2003020.View ArticlePubMedGoogle Scholar
- 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 ArticlePubMedGoogle Scholar
- Kalm E: Development of cattle breeding strategies in Europe. Arch Tierz. 2002, 45: 5-12.Google Scholar
- 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. 10.1186/1297-9686-36-2-163.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.Google Scholar
- Holmberg M, Andersson-Eklund L: Quantitative Trait Loci Affecting Health Traits in Swedish Dairy Cattle. J Dairy Sci. 2004, 87: 2653-2659.View ArticlePubMedGoogle Scholar
- 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 ArticlePubMedGoogle Scholar
- 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. 10.3168/jds.2007-0177.View ArticlePubMedGoogle Scholar
- 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 ArticlePubMedGoogle Scholar
- 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 ArticlePubMedGoogle Scholar
- 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.Google Scholar
- 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. 10.3168/jds.2008-1742.View ArticlePubMedGoogle Scholar
- 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.Google Scholar
- Pfitzer G: Signal Transduction in Smooth Muscle: Invited Review: Regulation of myosin phosphorylation in smooth muscle. J Appl Physiol. 2001, 91: 497-503.PubMedGoogle Scholar
- Walsh MP: Calmodulin and the regulation of smooth muscle contraction. Mol Cell Biochem. 1994, 135: 21-41. 10.1007/BF00925958.View ArticlePubMedGoogle Scholar
- Park CH, Valore EV, Waring AJ, Ganz T: Hepcidin, a Urinary Antimicrobial Peptide Synthesized in the Liver. J Biol Chem. 2001, 276: 7806-7810. 10.1074/jbc.M008922200.View ArticlePubMedGoogle Scholar
- 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. 10.1182/blood-2005-06-2259.PubMed CentralView ArticlePubMedGoogle Scholar
- Ganz T: Molecular Control of Iron Transport. J Am Soc Nephrol. 2007, 18: 394-400. 10.1681/ASN.2006070802.View ArticlePubMedGoogle Scholar
- Ivanov DB, Philippova MP, Tkachuk VA: Structure and Functions of Classical Cadherins. Biochemistry (Mosc). 2001, 66: 1174-1186. 10.1023/A:1012445316415.View ArticleGoogle Scholar
- Ebnet K: Organization of multiprotein complexes at cell-cell junctions. Histochem Cell Biol. 2008, 130: 1-20. 10.1007/s00418-008-0418-7.PubMed CentralView ArticlePubMedGoogle Scholar
- 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. 10.1146/annurev.cellbio.24.110707.175339.View ArticlePubMedGoogle Scholar
- Niessen CM: Tight Junctions/Adherens Junctions: Basic Structure and Function. J Invest Dermatol. 2007, 127: 2525-2532. 10.1038/sj.jid.5700865.View ArticlePubMedGoogle Scholar
- Ireton K: Entry of the bacterial pathogen Listeria monocytogenes into mammalian cells. Cell Microbiol. 2007, 9: 1365-1375. 10.1111/j.1462-5822.2007.00933.x.View ArticlePubMedGoogle Scholar
- Pizarro-Cerdá J, Cossart P: Bacterial Adhesion and Entry into Host Cells. Cell. 2006, 124: 715-727. 10.1016/j.cell.2006.02.012.View ArticlePubMedGoogle Scholar
- 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.Google Scholar
- 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 ArticlePubMedGoogle Scholar
- Cattle Genome Mapping Project. [http://www.marc.usda.gov/genome/cattle/cattle.html]
- Vincze T, Posfai J, Roberts RJ: NEBcutter: a program to cleave DNA with restriction enzymes. Nucl Acids Res. 2003, 31: 3688-3691. 10.1093/nar/gkg526.PubMed CentralView ArticlePubMedGoogle Scholar
- NEBcutter V2.0. [http://tools.neb.com/NEBcutter2/index.php]
- Green P, Falls K, Crooks S: Documentation of CRI-MAP, Version 2.4. 1990, Washington University School of Medicine, St. LouisGoogle Scholar
- 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-10.1186/1471-2164-6-74.PubMed CentralView ArticlePubMedGoogle Scholar
- 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. 10.1101/gr.2741704.PubMed CentralView ArticlePubMedGoogle Scholar
- 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. 10.1016/j.ygeno.2004.12.007.View ArticlePubMedGoogle Scholar
- NCBI Map Viewer. [http://www.ncbi.nlm.nih.gov/mapview/]
- 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. 10.1046/j.1365-2052.2000.00641.x.View ArticlePubMedGoogle Scholar
- Boehnke M, Lunetta K, Hauser E, Lange K, Uro J, VanderStoep J: RHMAP: Statistical Package for Multipoint Radiation Hybrid Mapping Version 3.0. 1996, University of Michigan, Ann Arbor, MichiganGoogle Scholar
- Baes C, Reinsch N: TIGER: A software system for fine-mapping quantitative trait loci. Arch Tierz. 2008, 51: 402-412.Google Scholar
- Reinsch N: A multiple-species, multiple-project database for genotypes at codominant loci. J Anim Breed Genet. 1999, 116: 425-435. 10.1046/j.1439-0388.1999.00217.x.View ArticleGoogle Scholar
- Meuwissen THE, Goddard ME: Fine mapping of quantitative trait loci using linkage disequilibria with closely linked marker loci. Genetics. 2000, 155: 421-430.PubMed CentralPubMedGoogle Scholar
- Meuwissen THE, Goddard ME: Prediction of identity by descent probabilities from marker-haplotypes. Genet Sel Evol. 2001, 33: 605-634. 10.1186/1297-9686-33-6-605.PubMed CentralView ArticlePubMedGoogle Scholar
- 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 CentralPubMedGoogle Scholar
- Baes C, Reinsch N: Computing the condensed conditional gametic QTL relationship matrix and its inverse. Arch Tierz. 2007, 50: 294-308.Google Scholar
- Gilmour AR, Gogel BJ, Cullis BR, Thompson R: ASReml User Guide Release 2.0. 2006, VSN International Ltd, Hemel Hempstead, HP1 1ES, UKGoogle Scholar
- Lander ES, Botstein D: Mapping Mendelian Factors Underlying Quantitative Traits Using Rflp Linkage Maps. Genetics. 1989, 121: 185-199.PubMed CentralPubMedGoogle Scholar
- Ott J: Analysis of Human Genetic Linkage. Revised edition. 1991, Baltimore Maryland: John Hopkins University Press, 65-68.Google Scholar
- NCBI HomePage. [http://www.ncbi.nlm.nih.gov/]
- NCBI BLAST. [http://blast.ncbi.nlm.nih.gov/Blast.cgi]
- NCBI Trace Archives. [http://www.ncbi.nlm.nih.gov/Traces/home/]
- 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. 10.1126/science.1169588.PubMed CentralView ArticleGoogle Scholar
- 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. 10.1111/j.1365-2052.1997.00115.x.View ArticleGoogle Scholar
- 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. 10.3168/jds.2007-0585.View ArticlePubMedGoogle Scholar
- 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. 10.1111/j.1365-2052.2005.01337.x.View ArticlePubMedGoogle Scholar
- 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. 10.1046/j.1439-0388.2000.00263.x.View ArticleGoogle Scholar
- Liu KJ, Muse SV: PowerMarker: an integrated analysis environment for genetic marker analysis. Bioinformatics. 2005, 21: 2128-2129. 10.1093/bioinformatics/bti282.View ArticlePubMedGoogle Scholar
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.