Robust tests for matched case-control genetic association studies

Background

Infectious disease of livestock continues to be a cause of substantial economic loss and has adverse welfare consequences in both the developing and developed world. New solutions to control disease are needed and research focused on the genetic loci determining variation in immune-related traits has the potential to deliver solutions. However, identifying selectable markers and the causal genes involved in disease resistance and vaccine response is not straightforward. The aims of this study were to locate regions of the bovine genome that control the immune response post immunisation. 195 F2 and backcross Holstein Charolais cattle were immunised with a 40-mer peptide derived from foot-and-mouth disease virus (FMDV). T cell and antibody (IgG1 and IgG2) responses were measured at several time points post immunisation. All experimental animals (F0, F1 and F2, n = 982) were genotyped with 165 microsatellite markers for the genome scan.

Results

Considerable variability in the immune responses across time was observed and sire, dam and age had significant effects on responses at specific time points. There were significant correlations within traits across time, and between IgG1 and IgG2 traits, also some weak correlations were detected between T cell and IgG2 responses. The whole genome scan detected 77 quantitative trait loci (QTL), on 22 chromosomes, including clusters of QTL on BTA 4, 5, 6, 20, 23 and 25. Two QTL reached 5% genome wide significance (on BTA 6 and 24) and one on BTA 20 reached 1% genome wide significance.

Conclusions

A proportion of the variance in the T cell and antibody response post immunisation with an FDMV peptide has a genetic component. Even though the antigen was relatively simple, the humoral and cell mediated responses were clearly under complex genetic control, with the majority of QTL located outside the MHC locus. The results suggest that there may be specific genes or loci that impact on variation in both the primary and secondary immune responses, whereas other loci may be specifically important for early or later phases of the immune response. Future fine mapping of the QTL clusters identified has the potential to reveal the causal variations underlying the variation in immune response observed.     

Genome-wide association study identifies two major loci affecting calving ease and growth-related traits in cattle

Identifying quantitative trait loci (QTL) underlying complex, low-heritability traits is notoriously difficult. Prototypical for such traits, calving ease is an important breeding objective of cattle (Bos taurus)-improving programs. To identify QTL underlying calving ease, we performed a genome-wide association study using estimated breeding values (EBVs) as highly heritable phenotypes for paternal calving ease (pCE) and related traits. The massively structured study population consisted of 1800 bulls of the German Fleckvieh (FV) breed. Two pCE-associated regions on bovine chromosomes (BTA) 14 and 21 (P = 5.72 × 10(-15) and P = 2.27 × 10(-8), respectively) were identified using principal components analysis to correct for population stratification. The two most significantly associated SNPs explain 10% of the EBV variation. Since marker alleles with negative effect on pCE have positive effects on growth-related traits, the QTL may exert their effects on the birthing process through fetal growth traits. The QTL region on BTA14 corresponds to a human chromosome (HSA) region that is associated with growth characteristics. The HSA region corresponding to the BTA21 pCE QTL is maternally imprinted and involved in the Prader-Willi and Angelman syndromes. Resequencing of positional candidate genes on BTA14 revealed a highly significantly (P = 1.96 × 10(-14)) associated polymorphism ablating a polyadenylation signal of the gene encoding ribosomal protein S20 (RPS20). Our study demonstrates the leverage potential of EBVs in unraveling the genetic architecture of lowly heritable traits.     

Quantitative trait loci for fertility traits in Finnish Ayrshire cattle

A whole genome scan was carried out to detect quantitative trait loci (QTL) for fertility traits in Finnish Ayrshire cattle. The mapping population consisted of 12 bulls and 493 sons. Estimated breeding values for days open, fertility treatments, maternal calf mortality and paternal non-return rate were used as phenotypic data. In a granddaughter design, 171 markers were typed on all 29 bovine autosomes. Associations between markers and traits were analysed by multiple marker regression. Multi-trait analyses were carried out with a variance component based approach for the chromosomes and trait combinations, which were observed significant in the regression method. Twenty-two chromosome-wise significant QTL were detected. Several of the detected QTL areas were overlapping with milk production QTL previously identified in the same population. Multi-trait QTL analyses were carried out to test if these effects were due to a pleiotropic QTL affecting fertility and milk yield traits or to linked QTL causing the effects. This distinction could only be made with confidence on BTA1 where a QTL affecting milk yield is linked to a pleiotropic QTL affecting days open and fertility treatments.     

Whole-genome scan for quantitative trait loci associated with birth weight, gestation length and passive immune transfer in a Holstein × Jersey crossbred population

We herein report results from a daughter design genome-scan study aiming to identify quantitative trait loci (QTL) associated with birth weight, direct gestation length and passive immune transfer in a backcross (Holstein × Jersey) × Holstein population. Two-hundred and seventy-six calves, offspring of seven crossbred sires, were genotyped for 161 microsatellite markers distributed along the 29 bovine autosomes. The genome scan was performed through interval mapping using an animal model in order to identify QTL accounting for phenotypic differences between individual animals. Based on significant chi-squared values, we identified putative QTL on BTA7 and BTA14 for gestation length, on BTA2, BTA6 and BTA14 for birth weight and on BTA20 for passive immune transfer. In total, these QTL accounted for 12%, 18% and 1% of the phenotypic variance in gestation length, birth weight and passive immune transfer respectively. We also report results from a supplementary and independent influential grand-daughter Holstein family. In this family, findings on BTA7 and BTA14 for direct gestation length were in agreement with results in the crossbred population. Two other regions on BTA6 and BTA21 putatively underlying QTL for direct gestation length variability were discovered with this analysis.     

Bovine Respiratory Syncytial Virus (BRSV) Pneumonia in Beef Calf Herds Despite Vaccination

The present report describes the clinical, pathological, serological and virological findings in calves from 2 larger Danish beef herds experiencing outbreaks of pneumonia. The calves had been vaccinated with an inactivated bovine respiratory syncytial virus (BRSV) vaccine 2 months prior to the outbreak. The clinical signs comprised nasal discharge, pyrexia, cough and increased respiratory rates. A total of 28 calves died in the 2 herds. The laboratory investigations revealed that BRSV was involved and probably initiated both outbreaks. Furthermore, the serological results suggested that the vaccine induced only sparse levels of antibodies probably due to the presence of maternally derived antibodies at the time of vaccination. Necropsy findings in 5 calves revealed changes typical for infectious pneumonia with involvment of BRSV. In conclusion, vaccination of calves against BRSV in 2 Danish beef herds failed to protect the calves against severe or even fatal BRSV mediated respiratory disease 2 months later.     

Mapping of a quantitative trait locus for beef marbling on bovine chromosome 9 in purebred Japanese black cattle

The goal of this study is to detect quantitative trait loci (QTL) for carcass traits applicable for a DNA-based breeding system in a Japanese Black cattle population. A purebred paternal half-sib family from a commercial line composed of 65 steers was initially analyzed using 188 informative microsatellites giving a 16-cM average interval covering 29 autosomes. A significant QTL for marbling was detected in the centromeric portion of bovine chromosome (BTA) 9. After additional marker genotyping across a larger sample size composed of 169 individuals, this locus was refined to a 20-cM confidence interval between microsatellites BM1227 (24 cM) and DIK2741 (50 cM) at a 1% chromosome-wise threshold. The allele substitution effect between Q and q for a beef marbling standard score (1 to 12 range) on BTA9 was 1.0 (5.7% of total phenotypic variance in QTL contribution in this family). This result provides a primary platform for a marker-assisted selection system of the beef marbling trait within the Japanese Black (Wagyu) cattle population.     

Identification of a 3.7-Mb region for a marbling QTL on bovine chromosome 4 by identical-by-descent and association analysis

QTL mapping for growth and carcass traits was performed using a paternal half-sib family composed of 325 Japanese Black cattle offspring. Nine QTL were detected at the 1% chromosome-wise significance level at a false discovery rate of less than 0.1. These included two QTL for marbling on BTA 4 and 18, two QTL for carcass weight on BTA 14 and 24, two QTL for longissimus muscle area on BTA 1 and 4, two QTL for subcutaneous fat thickness on BTA 1 and 15 and one QTL for rib thickness on BTA 6. Although the marbling QTL on BTA 4 has been replicated with significant linkages in two Japanese Black cattle sires, the three Q (more marbling) haplotypes, each inherited maternally, were apparently different. To compare the three Q haplotypes in more detail, high-density microsatellite markers for the overlapping regions were developed within the 95% CIs (65 markers in 44–78 cM). A detailed haplotype comparison indicated that a small region (<3.7 Mb) around 46 cM was shared between the Qs of the two sires, whose dams were related. An association of this region with marbling was shown by a regression analysis using the local population, in which the two sires were produced and this was confirmed by an association study using a population collected throughout Japan. These results strongly suggest that the marbling QTL on BTA 4 is located in the 3.7-Mb region at around 46 cM.     

Genome wide QTL mapping to identify candidate genes for carcass traits in Hanwoo (Korean Cattle)

Meat quality traits are the most economically important traits affecting the beef industry in Korea. We performed a whole genome quantitative trait locus (QTL) mapping study of carcass data in Hanwoo Korean cattle. Two hundred sixty-six Hanwoo steers from 65 sires were genotyped using a 10K Affymetrix SNP chip. The average SNP interval across the bovine genome was 1.5Mb. Associations between each individual SNP and four carcass traits [carcass weight (CWT), eye muscle area (EMA), back fat thickness (BFT), and marbling (MAR)] were assessed using a linear mixed model of each trait. Combined linkage and linkage disequilibrium analysis (LDLA) detected six potential QTL on BTA04, 06, 13, 16, 17, and 23 at the chromosome-wise level (P<0.05). Two MAR QTL were detected at 52.2 cM of BTA06 and 46.04 cM of BTA17. We identified three genes (ARAP2, LOC539460, and LOC511424) in the QTL region of BTA06 and seven genes (RPS14, SCARB1, LOC782103, BRI3BP, AACS, DHX37, and UBC) in the QTL region of BTA17. One significant QTL for CWT was detected at 100 cM on BTA04 and the corresponding QTL region spanned 1.7 cM from 99.7 to 101.4 cM. For EMA QTL, one significant QTL was detected at 3.9 cM of BTA23 and the most likely QTL interval was 1.4 cM, placing 15 candidate genes in the marker bracket. Finally, two QTL for BFT were identified at 68 cM on BTA13 and 24 cM on BTA16. The LPIN3 gene, which is functionally associated with lipodystrophy in humans, is located in the BFT QTL on BTA13. Thus, two potential candidate genes, acetoacetyl-CoA synthetase (AACS) and lipin (LPIN), were detected in QTL regions on BTA17 for MAR and BTA13 for BFT, respectively. In conclusion, LDLA analysis can be used to detect chromosome regions harboring QTL and candidate genes with a low density SNP panel, yielding relatively narrow confidence intervals regarding location.     

Fine mapping QTL for female fertility on BTA04 and BTA13 in dairy cattle using HD SNP and sequence data

Background

Female fertility is important for the maintenance of the production in a dairy cattle herd. Two QTL regions on BTA04 and on BTA13 previously detected in Nordic Holstein (NH) and validated in the Danish Jersey (DJ) and Nordic Red (NR) were investigated further in the present study to further refine the QTL locations. Refined QTL regions were imputed to the full sequence data. The genes in the regions were then studied to ascertain their possible effect on fertility traits.

Results

BTA04 was screened for number of inseminations (AIS), 56-day non-return rate (NRR), days from first to last insemination (IFL), and the interval from calving to first insemination (ICF) in the range of 38,257,758 to 40,890,784 bp, whereas BTA13 was screened for ICF only in the range from 21,236,959 to 46,150,079 with the HD bovine SNP array for NH, DJ and NR. No markers in the DJ and NR breeds reached significance. By analyzing imputed sequence data the QTL position on BTA04 was narrowed down to two regions in the NH. In these two regions a total of 9 genes were identified. BTA13 was analyzed using sequence data for the NH breed. The highest –log₁₀(P-value) was 19.41 at 33,903,159 bp. Two regions were identified: Region 1: 33,900,143-33,908,994 bp and Region 2: 34,051,815-34,056,728 bp. SNPs within and between these two regions were annotated as intergenic.

Conclusion

Screening BTA04 and BTA13 for female fertility traits in NH, NR and DJ suggested that the QTL for female fertility were specific for NH. A missense mutation in CD36 showed the strongest association with fertility traits on BTA04. The annotated SNPs on BTA13 were all intergenic variants. It is possible that BTA13 at this stage is poorly annotated such that the associated polymorphisms are located in as-yet undiscovered genes. Fertility traits are complex traits as many different biological and physiological factors determine whether a cow is fertile. Therefore it is not expected that there is a simple explanation with an obvious candidate gene but it is more likely a network of genes and intragenic variants that explain the variation of these traits.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-790) contains supplementary material, which is available to authorized users.     

Mapping of a Quantitative Trait Locus for Beef Marbling on Bovine Chromosome 9 in Purebred Japanese Black Cattle

The goal of this study is to detect quantitative trait loci (QTL) for carcass traits applicable for a DNA-based breeding system in a Japanese Black cattle population. A purebred paternal half-sib family from a commercial line composed of 65 steers was initially analyzed using 188 informative microsatellites giving a 16-cM average interval covering 29 autosomes. A significant QTL for marbling was detected in the centromeric portion of bovine chromosome (BTA) 9. After additional marker genotyping across a larger sample size composed of 169 individuals, this locus was refined to a 20-cM confidence interval between microsatellites BM1227 (24 cM) and DIK2741 (50 cM) at a 1% chromosome-wise threshold. The allele substitution effect between Q and q for a beef marbling standard score (1 to 12 range) on BTA9 was 1.0 (5.7% of total phenotypic variance in QTL contribution in this family). This result provides a primary platform for a marker-assisted selection system of the beef marbling trait within the Japanese Black (Wagyu) cattle population.     

Genetic mapping of quantitative trait loci for meat quality and muscle metabolic traits in cattle

A whole‐genome scan was carried out in New Zealand and Australia to detect quantitative trait loci (QTL) for live animal and carcass composition traits and meat quality attributes in cattle. Backcross calves (385 heifers and 398 steers) were generated, with Jersey and Limousin backgrounds. The New Zealand cattle were reared and finished on pasture, whilst Australian cattle were reared on grass and finished on grain for at least 180 days. This paper reports on meat quality traits (tenderness measured as shear force at 4–5 ages on two muscles as well as associated traits of meat colour, pH and cooking loss) and a number of metabolic traits. For meat quality traits, 18 significant QTL (P < 0.05), located in nine linkage groups, were detected on a genome‐wise basis, in combined‐sire (seven QTL) or within‐sire analyses (11 QTL). For metabolic traits, 11 significant QTL (P < 0.05), located in eight linkage groups, were detected on a genome‐wise basis, in combined‐sire (five QTL) or within‐sire analyses (six QTL). BTA2 and BTA3 had QTL for both metabolic traits and meat quality traits. Six significant QTL for meat quality and metabolic traits were found at the proximal end of chromosome 2. BTA2 and BTA29 were the most common chromosomes harbouring QTL for meat quality traits; QTL for improved tenderness were associated with Limousin‐derived and Jersey‐derived alleles on these two chromosomes, respectively.     

Revisiting the quantitative trait loci for milk production traits on BTA6

A parallel association study was performed in two independent cattle populations based on 41 validated, targeted single nucleotide polymorphisms (SNPs) and four microsatellite markers to re-evaluate the multiple quantitative trait loci (QTL) architecture for milk performance on bovine chromosome 6 (BTA6). Two distinct QTL located in the vicinity of the middle region of BTA6, but differing unambiguously regarding their effects on milk composition and yield traits were validated in the German Holstein population. A highly significant association of the protein variant ABCG2 p.Tyr581Ser with milk composition traits reconfirmed the causative molecular relevance of the ABCG2 gene in QTL region 1, whereas in QTL region 2, significant and tentative associations between gene variants RW070 and RW023 (located in the promoter region and exon 9 of the PPARGC1A gene for milk yield traits) were detected. For the German Fleckvieh population, only RW023 showed a tentative association with milk yield traits, whereas those loci with significant effects in German Holsteins (ABCG2 p.Tyr581Ser, RW070) showed fixed alleles. Even though our new data highlight two variants in the PPARGC1A gene (RW023, RW070) in QTL region 2, based on the results of our study, currently no unequivocal conclusion about the causal background of this QTL affecting milk yield traits can be drawn. Notably, the German Holstein and Fleckvieh populations, known for their divergent degree of dairy type, differ substantially in their allele frequencies for the growth-associated NCAPG p.Ile442Met locus.     

Mapping of quantitative trait loci for carcass traits in a Japanese Black (Wagyu) cattle population

To detect quantitative trait loci (QTL) that influence economically important traits in a purebred Japanese Black cattle population, we performed a preliminary genome-wide scan using 187 microsatellite markers across a paternal half-sib family composed of 258 offspring. We located six QTL at the 1% chromosome-wise level on bovine chromosomes (BTA) 4, 6, 13, 14 and 21. A second screen of these six QTL regions using 138 additional paternal offspring half-sib from the same sire, provided further support for five QTL: carcass weight on BTA14 (22-39 cM), one for rib thickness on BTA6 (27-58 cM) and three for beef marbling score (BMS) on BTA4 (59-67 cM), BTA6 (68-89 cM) and BTA21 (75-84 cM). The location of QTL for subcutaneous fat thickness on BTA13 was not supported by the second screen (P > 0.05). We determined that the combined contribution of the three QTLs for BMS was 10.1% of the total variance. The combined phenotypic average of these three Q was significantly different (P < 0.001) from those of other allele combinations. Analysis of additional half-sib families will be necessary to confirm these QTL.     

Whole genome scan to detect quantitative trait loci for bovine milk protein composition

The objective of this study was to perform a whole genome scan to detect quantitative trait loci (QTL) for milk protein composition in 849 Holstein–Friesian cows originating from seven sires. One morning milk sample was analysed for the major milk proteins using capillary zone electrophoresis. A genetic map was constructed with 1341 single nucleotide polymorphisms, covering 2829 centimorgans (cM) and 95% of the cattle genome. The chromosomal regions most significantly related to milk protein composition ( P _genome < 0.05) were found on Bos taurus autosomes (BTA) 6, 11 and 14. The QTL on BTA6 was found at about 80 cM, and affected α_S1-casein, α_S2-casein, β-casein and κ-casein. The QTL on BTA11 was found at 124 cM, and affected β-lactoglobulin, and the QTL on BTA14 was found at 0 cM, and affected protein percentage. The proportion of phenotypic variance explained by the QTL was 3.6% for β-casein and 7.9% for κ-casein on BTA6, 28.3% for β-lactoglobulin on BTA11, and 8.6% for protein percentage on BTA14. The QTL affecting α_S2-casein on BTA6 and 17 showed a significant interaction. We investigated the extent to which the detected QTL affecting milk protein composition could be explained by known polymorphisms in β-casein , κ -casein , β-lactoglobulin and DGAT1 genes. Correction for these polymorphisms decreased the proportion of phenotypic variance explained by the QTL previously found on BTA6, 11 and 14. Thus, several significant QTL affecting milk protein composition were found, of which some QTL could partially be explained by polymorphisms in milk protein genes.     

Confirmation and fine‐mapping of clinical mastitis and somatic cell score QTL in Nordic Holstein cattle

A genome‐wide association study of 2098 progeny‐tested Nordic Holstein bulls genotyped for 36 387 SNPs on 29 autosomes was conducted to confirm and fine‐map quantitative trait loci (QTL) for mastitis traits identified earlier using linkage analysis with sparse microsatellite markers in the same population. We used linear mixed model analysis where a polygenic genetic effect was fitted as a random effect and single SNPs were successively included as fixed effects in the model. We detected 143 SNP‐by‐trait significant associations (P < 0.0001) on 20 chromosomes affecting mastitis‐related traits. Among them, 21 SNP‐by‐trait combinations exceeded the genome‐wide significant threshold. For 12 chromosomes, both the present association study and the previous linkage study detected QTL, and of these, six were in the same chromosomal locations. Strong associations of SNPs with mastitis traits were observed on bovine autosomes 6, 13, 14 and 20. Possible candidate genes for these QTL were identified. Identification of SNPs in linkage disequilibrium with QTL will enable marker‐based selection for mastitis resistance. The candidate genes identified should be further studied to detect candidate polymorphisms underlying these QTL.     

Identification of bovine QTL for growth and carcass traits in Japanese Black cattle by replication and identical-by-descent mapping

To map quantitative trait loci (QTL) for growth and carcass traits in a purebred Japanese Black cattle population, we conducted multiple QTL analyses using 15 paternal half-sib families comprising 7860 offspring. We identified 40 QTL with significant linkages at false discovery rates of less than 0.1, which included 12 for intramuscular fat deposition called marbling and 12 for cold carcass weight or body weight. The QTL each explained 2%–13% of the phenotypic variance. These QTL included many replications and shared hypothetical identical-by-descent (IBD) alleles. The QTL for CW on BTA14 was replicated in five families with significant linkages and in two families with a 1% chromosome-wise significance level. The seven sires shared a 1.1-Mb superior Q haplotype as a hypothetical IBD allele that corresponds to the critical region previously refined by linkage disequilibrium mapping. The QTL for marbling on BTA4 was replicated in two families with significant linkages. The QTL for marbling on BTA6, 7, 9, 10, 20, and 21 and the QTL for body weight on BTA6 were replicated with 1% and/or 5% chromosome-wise significance levels. There were shared IBD Q or q haplotypes in the marbling QTL on BTA4, 6, and 10. The allele substitution effect of these haplotypes ranged from 0.7 to 1.2, and an additive effect between the marbling QTL on BTA6 and 10 was observed in the family examined. The abundant and replicated QTL information will enhance the opportunities for positional cloning of causative genes for the quantitative traits and efficient breeding using marker-assisted selection. Electronic Supplementary Material Electronic Supplementary material is available for this article at and accessible for authorised users.     

A simplified QTL mapping approach for screening and mapping of novel AFLP markers associated with beef marbling

Genome screening of quantitative trait loci (QTL) for a complex trait is usually costly and highly laborious, as it requires a large number of markers spanning the whole genome. Here we present a simplified approach for screening and mapping of QTL-linked markers for beef marbling using a WagyuxLimousin F(2) reference population. This simplified approach involves integration of the amplified fragment length polymorphism (AFLP) with DNA pooling and selective genotyping and comparative bioinformatics tools. AFLP analysis on two high and two low marbling DNA pools yielded ten visually different markers. Among them, four were confirmed based on individual AFLP validation. Sequencing and in silico characterization assigned two of these AFLP markers to bovine chromosomes 1 (BTA1) and 13 (BTA13), which are orthologous to human chromosomes HSA21q22.2 and HSA10p11.23 with both regions harboring QTL for obesity-related phenotypes. Both AFLP markers showed significantly large additive genetic effects (0.28+/-0.11 on BTA1 and 0.54+/-0.21 on BTA13) on beef-marbling score (BMS) (P<0.05). Overall, this approach is less time consuming, inexpensive and in particular, suitable for screening and mapping QTL-linked markers when targeting one or a few complex traits.     

Detection of QTL for milk production traits in cattle by application of a specifically developed marker map of BTA6

Interval mapping was carried out to identify quantitative trait loci (QTL) for milk production traits in five granddaughter design families of the German Holstein population. Fourteen randomly generated markers spanning the whole of BTA6 and six targeted microsatellite markers from BTA6q21-31 were included in the analysis. In one family a QTL with effects on milk fat yield and milk protein yield was mapped to the interval TGLA37-FBN13 (3 CM proximal to FBN13, lodscore 3.22) in the middle part of the chromosome. Although there are several reports about QTL with effects on milk production traits on BTA6 in the literature, a QTL with effects on milk fat and milk protein yield has not been previously described.     

Quantitative trait locus-by-environment interaction for milk yield traits on Bos taurus autosome 6

Genotype-by-environment interactions for production traits in dairy cattle have often been observed, while QTL analyses have focused on detecting genes with general effects on production traits. In this study, a QTL search for genes with environmental interaction for the traits milk yield, protein yield, and fat yield were performed on Bos taurus autosome 6 (BTA6), also including information about the previously investigated candidate genes ABCG2 and OPN. The animals in the study were Norwegian Red. Eighteen grandsires and 716 sires were genotyped for 362 markers on BTA6. Every marker bracket was regarded as a putative QTL position. The effects of the candidate genes and the putative QTL were modeled as a regression on an environmental parameter (herd year), which is based on the predicted herd-year effect for the trait. Two QTL were found to have environmentally dependent effects on milk yield. These QTL were located 3.6 cM upstream and 9.1 cM downstream from ABCG2. No environmentally dependent QTL was found to significantly affect protein or fat yield.     

Mapping of multiple quantitative trait loci affecting bovine spongiform encephalopathy

A whole-genome scan was conducted to map quantitative trait loci (QTL) for BSE resistance or susceptibility. Cows from four half-sib families were included and 173 microsatellite markers were used to construct a 2835-cM (Kosambi) linkage map covering 29 autosomes and the pseudoautosomal region of the sex chromosome. Interval mapping by linear regression was applied and extended to a multiple-QTL analysis approach that used identified QTL on other chromosomes as cofactors to increase mapping power. In the multiple-QTL analysis, two genome-wide significant QTL (BTA17 and X/Y(ps)) and four genome-wide suggestive QTL (BTA1, 6, 13, and 19) were revealed. The QTL identified here using linkage analysis do not overlap with regions previously identified using TDT analysis. One factor that may explain the disparity between the results is that a more extensive data set was used in the present study. Furthermore, methodological differences between TDT and linkage analyses may affect the power of these approaches.