Cross-breed comparisons identified a critical 591-kb region for bovine carcass weight QTL (CW-2) on chromosome 6 and the Ile-442-Met substitution in NCAPG as a positional candidate
© Setoguchi et al; licensee BioMed Central Ltd. 2009
Received: 17 February 2009
Accepted: 4 August 2009
Published: 4 August 2009
Growth-related traits have been mapped on bovine chromosome 6 (BTA 6) in various bovine breed populations. We previously mapped a significant quantitative trait locus (QTL) for carcass and body weight (CW-2) between 38 and 55 cM on BTA 6 using a Japanese Black half-sib family. Additional QTL mapping studies detected four QTL for body or carcass weight that overlapped with CW-2 in Japanese Black and Japanese Brown half-sib families. To map the region in greater detail, we applied cross-breed comparisons of haplotypes that have been shown to be powerful in canine.
We used 38 microsatellite markers to search for a shared Q (increasing carcass and/or body weight) haplotype within the 17-cM CW-2 region among five sires. Linkage disequilibrium mapping using maternal alleles of the offspring showed that an 815-kb shared Q haplotype was associated with body or carcass weight in both breeds. The addition of 43 single nucleotide polymorphism (SNP) markers narrowed the region to 591 kb containing 4 genes. The SNP changing Ile-442 to Met in NCAPG (chromosome condensation protein G) was significantly associated with carcass weight (p < 1.2 × 10-11) in a large Japanese Black population as well as in the five families. The Q allele of the SNP was also associated with a larger longissimus muscle area and thinner subcutaneous fat thickness in steers of all five families, indicating that the CW-2 locus is pleiotropic and favorable for marker-assisted selection of beef cattle.
A 591-kb critical region for CW-2 was identified. The SNP changing Ile-442 to Met in NCAPG (chromosome condensation protein G) can be used as a positional candidate of CW-2 for marker-assisted selection.
Body size is one of the most visible animal characteristics and many genes can affect body size. In cattle, body size is correlated with meat quantity, an economically important trait that varies within and across breeds. We previously performed bovine quantitative trait locus (QTL) mapping for growth and carcass traits using Japanese Black paternal half-sib families constructed from a commercial population [1–3]. In these studies, a carcass weight QTL, CW-1 on bovine chromosome 14 (BTA 14), was detected in five families with significant linkages and successfully narrowed down to a 1.1-Mb region by identical-by-descent mapping  and linkage disequilibrium (LD) mapping . In contrast with CW-1, another significant QTL for carcass and body weight, designated CW-2, on BTA 6 was replicated at a 1% chromosome-wise significance level, but the two significantly segregating sires had no apparent shared Q haplotypes . LD mapping only narrowed the region down to 13.1 Mb (see Results). We thus changed strategies to narrow down the CW-2 region. We recently detected carcass weight QTL in regions that overlapped with CW-2 in Japanese Brown populations. The Japanese Black breed was established in 1948, basically from indigenous populations of the Japanese Islands, while the Japanese Brown breed originated from indigenous populations of the Korean Peninsula. Approximately 100 years ago, the Korean cattle were imported to Japan and crossed several times with Simmental bulls, followed by the establishment of the Japanese Brown breed in 1948. Therefore, it can be estimated that the two indigenous populations have survived separately for thousands of years. It will be interesting if a common carcass weight QTL is present across the entire Bos taurus population. Indeed, the CW-2 region was repeatedly highlighted for its association with growth-related traits in a mixed breed population , postweaning growth in a beef cattle population , and weight and body length at birth in a Charolais × Holstein cross-bred population . These data indicate that CW-2 may be shared across breeds. If this is the case, cross-breed comparisons might be useful for fine-mapping to identify a shared hypothetical identical-by-descent haplotype.
In dogs, cross-breed comparisons have been proposed for fine-mapping to pinpoint disease-related genes . Canine LD patterns reflect two bottlenecks in dog history (early domestication and recent breed creation) and long-range breed-specific haplotype blocks retain the underlying short-range ancestral haplotype blocks, suggesting that genetic risk factors are shared across breeds. This situation provides for an efficient mapping strategy: initial mapping within breeds and subsequent fine-mapping by cross-breed comparisons. The strategy was validated in coat color and hair ridge studies , in which cross-breed comparisons successfully refined the location of the coat color locus from 800 kb to 100 kb within the pigmentation-related gene. A recent study in cattle also revealed footprints of ancestral LD at short distances (< 10 kb) and these ancestral blocks are organized into larger blocks of a few hundred kilobases within a breed .
Here, we demonstrate that cross-breed comparisons are also efficient for fine-mapping of QTL in cattle. The CW-2 locus was narrowed down from a 17-cM region to a 591-kb region shared by the Japanese Black and Japanese Brown breeds. Furthermore, a single nucleotide polymorphism (SNP) changing Ile-442 to Met in NCAPG (chromosome condensation protein G) was identified as a positional candidate. The NCAPG-LCORL locus was recently mapped as a locus associated with human adult height [11, 12]. The findings of the present study indicate that genetic variations in NCAPG may affect body size in cattle as well as in humans.
Carcass and/or body weight QTL in Japanese Black and Japanese Brown cattle
Identical-by-state haplotype associated with carcass and body weight
In Japanese Brown families, offspring from Sires D and E were combined to maximize the number of offspring. Carcass weight was higher in DIK4852-DIK9016 (p = 0.011) and in DIK9019-MNB202 (p = 0.0052 at DIK9019–DIK9020) in the maternal Q3-harboring steers than in the non-Q3- harboring steers among the offspring with the paternally inherited q haplotype (Figure 3B). The two regions were in linkage disequilibrium. Among 161 steers that inherited the paternal q haplotype, 50 harbored maternal Q3 in DIK4852-DIK9016, 31 of which harbored maternal Q3 in DIK9019–DIK9020. On the other hand, all of the steers that harbored maternal Q3 in DIK9019–DIK9020 also harbored maternal Q3 in DIK4852-DIK9016. The DIK4852-DIK9016 region was included in the narrowed-down CW-2 region and corresponded to one of the possible IBS regions, whereas the DIK9019-MNB202 region was included in a homozygous region in Sire C. An association of the DIK4852-DIK9016 region with carcass weight was also observed among offspring with the paternally inherited Q haplotype.
These results strongly suggest that the causal alleles in the three QTL haplotypes Q1, Q2, and Q3 are identical, and that the IBS region between DIK4852 and DIK9017 might be responsible for the QTL.
Ile-442-Met Mutation as a positional candidate
The possible IBS region between DIK4852 and DIK9017 corresponds to the 17.2 to 17.9-Mb region on human chromosome 4 that contains six genes (LAP3, MED28, FAM184B, C6H4orf30, NCAPG, and LCORL) [13–15]. To examine the region in detail, SNPs were searched for by resequencing Sire A in the coding regions of the six orthologous bovine genes (LAP3, MED28, LOC523874, C6H4orf30, NCAPG, and LOC540095) and genotyped in the five sires (Figure 2B; Additional file 3). This ensured and refined the IBS Q region to an 879-kb interval flanked by SNP-0 and BTC-071108**. In addition, the telomeric part of the region from BTA-09966** to BTC-071108** was homozygous in Sire C, resulting in the narrowing of the CW-2 region to a 591-kb interval flanked by SNP-0 and BTA-09966** (Figure 2B; Additional file 3).
To assess the association of these SNPs with carcass weight, 187 steers selected from the 4.7% extremes of 7990 Japanese Black steers were genotyped. SNP-9 of NCAPG had the strongest association (p < 1.2 × 10-11, Fisher's exact test for independence; frequency of G-allele, 0.14; Figure 4A; Additional file 3), and it was in strong linkage disequilibrium with the four intronic SNPs (Figure 4BC). These data indicated that SNP-9 is the most appropriate CW-2 marker among the examined SNPs as a candidate for the causative sequence variation, although three other genes are in the same LD block and cannot be excluded as candidate genes for CW-2.
Pleiotropy of CW-2
QTL mapping studies detected QTL for the longissimus muscle area (LMA) and subcutaneous fat thickness (SFT) in the regions overlapping with CW-2 at less than 5% chromosome-wise significance levels in families of Sires A and B (Additional file 4). A QTL for LMA was also detected in Sire D's family at the 1% chromosome-wise significance level (Additional file 4). These QTL peaks were also extinguished in a fixed effect QTL model accounting for the effect of SNP-9 (Additional file 4). Therefore the SNP-9 genotypes were tested for LMA and SFT in each family (Figure 5). The G/G- or G/T-genotypes had a significantly greater average LMA than the T/T-genotype in the steers of the five families and the cows of Sire D's offspring (p < 0.01, t-test), although the correlation between carcass weight and LMA traits was not high (r = 0.26–0.54). On the other hand, the G/G-genotype had a significantly lower average SFT than the T/T or G/T-genotypes in the steers of the five families (p < 0.039, t-test). This is a particularly interesting finding because the correlation between carcass weight and SFT traits was r = 0.18–0.39 in each family (Additional file 5), meaning that greater carcass weight was correlated with thicker SFT. For example, a Q allele of another carcass weight QTL CW-1 on BTA 14 is associated with greater carcass weight and thicker SFT . In CW-2, however, the SNP-9 genotype was associated with thinner SFT in the steers of the five families, indicating that the CW-2 locus is pleiotropic. The results further support the notion that the CW-2 QTL in the Japanese Black population is the same as the QTL in a respective position in the Japanese Brown population.
Here we show that cross-breed comparisons are useful for fine-mapping of the QTL in cattle, as recently demonstrated in the fine-mapping of a monogenic trait locus in dogs . We detected a carcass weight QTL on BTA 6 in two breed populations, Japanese Black and Japanese Brown. The genetic distance between the two cattle breeds can be thousands of years, because their ancestors have survived separately on the Japanese Islands and the Korean Peninsula. Linkage disequilibrium mapping using the two breeds narrowed the CW-2 QTL region to a 591-kb section. Comparison of the Q haplotypes inherited in different lineages is, therefore, a useful strategy for narrowing down a QTL.
One of the critical issues should be whether the CW-2 Q in the Japanese Black population is identical to the QTL in the Japanese Brown population. In this study, the effects of the CW-2 on carcass weight, LMA, and SFT traits were similar in the two breed populations. Namely, the CW-2 Q allele was associated with greater carcass weight, larger LMA, and thinner SFT. Importantly, although greater carcass weight is generally associated with a thicker SFT, in both populations the CW-2 Q allele was associated with a thinner SFT. In addition, QTL for growth-related traits were mapped around the CW-2 region in various cattle populations [5–7], suggesting that the CW-2 Q allele is widely distributed in Bos taurus populations.
The 591-kb critical region contains four genes, LOC523874, C6H4orf30, NCAPG, and LOC540095 (the bovine ortholog of LCORL). The NCAPG-LCORL region was recently identified as a QTL for human adult height [11, 12]. An SNP causing Ile-442 to Met in NCAPG was detected as a positional candidate for the quantitative trait nucleotide for CW-2. The NCAPG gene encodes chromosome condensation subunit G, which is a catalytic subunit of the mammalian condensin I complex. Disruption of the NCAPG ortholog (cnd3) causes defective chromosome condensation in fission yeast , while in Drosophila, the NCAPG ortholog (Cap-G) is essential for chromosome condensation in the metaphase of single, unreplicated sister chromatids, and also has a role during the interphase in regulating heterochromatic gene expression . In Hela cells, NCAPG interacts with HSF2 to mediate HSPA1A (hsp70i) bookmarking . The NCAPG expression is upregulated in more aggressive metastatic melanomas than in less aggressive primary melanomas . These findings indicate that NCAPG may affect cell proliferation and growth through regulating cell cycle and chromosome condensation. Further studies are required to reveal the molecular mechanism that links NCAPG to growth and body size.
We identified a 591-kb critical region for CW-2. Among the four genes located within the region, NCAPG is a candidate for the causative gene because it has an amino acid substitution that is significantly associated with carcass weight. The SNP for the amino acid substitution provides a useful CW-2 marker for marker-assisted selection. The pleiotropic characteristic of CW-2, especially the effect on SFT, is clearly different from a previously fine-mapped carcass weight QTL, CW-1 . Therefore, the two QTL probably affect carcass weight in different ways and both will be useful for breeding beef cattle.
DNA samples and phenotype data
Paternal half-sib families of Sires A through E were constructed from carcass data and pedigree records collected by the Japan Wagyu Register Association (Kyoto, Japan) and Kumamoto Union of Livestock and Agriculture Cooperation (Kumamoto, Japan). Offspring of each sire were reared in different herds in a prefecture. The offspring of each sire were collected over a period ranging from 2 to 5 years. The average slaughter age was 878 days for the offspring of Japanese Black Sires A through C and 750 days for the offspring of Japanese Brown Sires D and E. Sire DNA was obtained from semen. Offspring DNA samples were collected from adipose tissues around the kidney at the slaughterhouses or blood at the individual farmer's houses. Six traits were analysed in this study: body weight at slaughter (Sire C only), cold carcass weight, LMA, rib thickness (thickness of a muscle layer in a rib of beef), SFT, and marbling. These traits were systematically measured by certified graders and recorded at the slaughterhouses in Japan.
where y is the vector of phenotypic value, X is the design matrix of fixed effects composed of sex, slaughter year, age (day), and probability of having the Q phase at a given location(Prob(Q)), b is the vector of fixed effects, and e is the residual error. b was estimated by the least squares method. An F-statistic value at each position was calculated from the residual sum square regressed with Prob(Q), and the total residual sum square without Prob(Q). The analysis was performed at 2-cM intervals along each chromosome. To evaluate whether the QTL effect was well estimated, the information content was calculated as a variance of Prob(Q) divided by 0.25, which was the maximum possible value of Prob(Q) . The allele substitution effect from q to Q was calculated as an estimator of the cofactor for Prob(Q) in the b. The contribution ratio was calculated as a proportion of the trait variance explained by the paternal allele substitution from q to Q. A threshold for significance of the F-statistic value was obtained by 10,000 random permutations of the phenotypic data . The 95% CI of the QTL locations was calculated using the bootstrapping method . Briefly, a set of offspring was chosen so as to be the same number as the original half-sibs by resampling from the original half-sibs. Resampling was repeated 10,000 times. The position of the F-statistic peak in each bootstrapping was collected. The CI was determined by the distribution of the peaks. Therefore, the CI may be fragmented into separated regions, and not a single contiguous region.
Microsatellite development and genotyping
Microsatellites were searched for in the genomic sequences and the primers were designed using Primer 3 . Twenty microsatellite markers from DIK9001 to DIK9020 were developed in this study. The markers were anchored to the BAC fingerprint map  by a BLAST search (BLASTN) against bovine BAC sequences (HTGS)  or by screening bovine BAC libraries, RPCI-42  and CHORI-240 . Marker information, such as primer sequences, genomic positions, and locations in the BAC fingerprint map, is shown in Additional file 2. Genotyping was performed using polymerase chain reaction (PCR) with a fluorescent-labeled reverse primer, followed by electrophoresis using ABI 3730 DNA analyzer (Applied Biosystems, Foster City, CA) and analysis using GeneMapper software (Applied Biosystems). The sires and their offspring were genotyped to determine the phase of the sires' chromosomes.
Linkage disequilibrium mapping
Linkage disequilibrium mapping was performed using the maternal alleles of the offspring steers. First, offspring were divided into two groups according to paternally inherited haplotypes (Q or q) of the CW-2 region (38–55 cM). When a crossover was observed or suspected within the region in the paternal haplotype, the offspring were eliminated. The numbers of paternal Q- and q-inherited steers were 126 and 117 in Sire C's family, 106 and 94 in Sire D's family, and 104 and 67 in Sire E's family, respectively. The maternal allele was determined by subtracting the paternal allele from the genotype. The markers for which the number of alleles was less than four (Sire C family) or three (families of Sires D and E), or heterozygosity was less than 0.4 were eliminated. The effect of pairwise Q alleles on body or carcass weight was evaluated in each group using a t-test. Phenotypic values for carcass and body weight were adjusted based on slaughter year and age in the family. The average of the adjusted carcass or body weight was compared between the steers that harbored two consecutive Q alleles in the maternal alleles and the other steers, in each group. In Sire C family, the test was performed for Q1 and Q2 alleles in the regions where corresponding Q frequencies were higher than 5%. In Japanese Brown families, offspring from Sires D and E were combined to maximize the number of offspring. Carcass weight was not adjusted between families, because the effect of sire was not significant (p = 0.21). Frequencies of two consecutive Q3 alleles in the maternal alleles were more than 4.3% in the entire region.
SNP discovery and genotyping
Six genes (LAP3, MED28, LOC523874, C6H4orf30, NCAPG, and LOC540095) were located between DIK4852 and DIK9017. SNPs in the coding regions of the genes were searched for by resequencing Sire A. SNPs chosen from a 120 K SNP collection  and a commercially available bovine 50 K SNP chip (Illumina, San Diego, CA) were also used. The sires and their offspring were genotyped to determine the phases of the sires' chromosomes. SNPs were genotyped by direct sequencing of the genomic PCR products. As for SNP (Indel) 19, PCR was performed using a fluorescent-labeled reverse primer and genotyped in the same way as a microsatellite marker. SNP information, such as alleles, genomic positions, and primer sequences used for amplification and sequencing, is shown in Additional file 3.
A population consisting of 7990 Japanese Black steers was constructed from DNA samples collected at one slaughterhouse over 7 years. These cattle were reared in different herds throughout Japan. To maximize the power and minimize the pedigree effect of a selected population, steers were selected from the 4.7% extremes, from many different sires, and did not include more than 5 half-sibs in each extreme. The selected population consisted of 90 offspring (> 570 kg) from 41 sires and 97 offspring (< 410 kg) from 52 sires, among which 18 sires were common. The 187 steers were genotyped with SNPs. The association of each SNP with carcass weight was assessed by Fisher's exact test for independence using a 2 × 2 contingency table consisting of an observed number of each allele in each extreme population.
Calculation of LD
SNPs were confirmed not to have deviated from Hardy-Weinberg equilibrium (p > 0.05). Haplotype frequencies consisting of two SNPs were estimated using the expectation-maximization algorithm of Dempster et al. . An LD coefficient (r2) between two SNPs was calculated from the estimated haplotype frequencies.
The authors thank Y. Kaneuchi and A. Narushima for technical assistance. The work was partly supported by the Ministry of Agriculture, Forestry, and Fishery, Japan, and by the Japan Racing and Livestock Promotion Foundation.
- Mizoshita K, Watanabe T, Hayashi H, Kubota C, Yamakuchi H, Todoroki J, Sugimoto Y: Quantitative trait loci analysis for growth and carcass traits in a half-sib family of purebred Japanese Black (Wagyu) cattle. J Anim Sci. 2004, 82: 3415-3420.PubMedGoogle Scholar
- Mizoguchi Y, Watanabe T, Fujinaka K, Iwamoto E, Sugimoto Y: Mapping of quantitative trait loci for carcass traits in a Japanese Black (Wagyu) cattle population. Anim Genet. 2006, 37: 51-54. 10.1111/j.1365-2052.2005.01367.x.View ArticlePubMedGoogle Scholar
- Takasuga A, Watanabe T, Mizoguchi Y, Hirano T, Ihara N, Takano A, Yokouchi K, Fujikawa A, Chiba K, Kobayashi N, et al: Identification of bovine QTL for growth and carcass traits in Japanese Black cattle by replication and identical-by-descent mapping. Mamm Genome. 2007, 18: 125-136. 10.1007/s00335-006-0096-5.View ArticlePubMedGoogle Scholar
- Mizoshita K, Takano A, Watanabe T, Takasuga A, Sugimoto Y: Identification of a 11-Mb region for a carcass weight QTL on bovine Chromosome 14. Mamm Genome. 2005, 16: 532-537. 10.1007/s00335-005-0024-0.View ArticlePubMedGoogle Scholar
- Allan MF, Thallman RM, Cushman RA, Echternkamp SE, White SN, Kuehn LA, Casas E, Smith TP: Association of a single nucleotide polymorphism in SPP1 with growth traits and twinning in a cattle population selected for twinning rate. J Anim Sci. 2007, 85: 341-347. 10.2527/jas.2006-460.View ArticlePubMedGoogle Scholar
- White SN, Casas E, Allan MF, Keele JW, Snelling WM, Wheeler TL, Shackelford SD, Koohmaraie M, Smith TP: Evaluation in beef cattle of six deoxyribonucleic acid markers developed for dairy traits reveals an osteopontin polymorphism associated with postweaning growth. J Anim Sci. 2007, 85: 1-10. 10.2527/jas.2006-314.View ArticlePubMedGoogle Scholar
- Gutierrez-Gil B, Williams JL, Homer D, Burton D, Haley CS, Wiener P: Search for QTL affecting growth and carcass traits in a cross population of beef and dairy cattle. J Anim Sci. 2009, 87: 24-36. 10.2527/jas.2008-0922.View ArticlePubMedGoogle Scholar
- Lindblad-Toh K, Wade CM, Mikkelsen TS, Karlsson EK, Jaffe DB, Kamal M, Clamp M, Chang JL, Kulbokas EJ, Zody MC, et al: Genome sequence, comparative analysis and haplotype structure of the domestic dog. Nature. 2005, 438: 803-819. 10.1038/nature04338.View ArticlePubMedGoogle Scholar
- Karlsson EK, Baranowska I, Wade CM, Salmon Hillbertz NH, Zody MC, Anderson N, Biagi TM, Patterson N, Pielberg GR, Kulbokas EJ, et al: Efficient mapping of mendelian traits in dogs through genome-wide association. Nat Genet. 2007, 39: 1321-1328. 10.1038/ng.2007.10.View ArticlePubMedGoogle Scholar
- Gautier M, Faraut T, Moazami-Goudarzi K, Navratil V, Foglio M, Grohs C, Boland A, Garnier JG, Boichard D, Lathrop GM, et al: Genetic and haplotypic structure in 14 European and African cattle breeds. Genetics. 2007, 177: 1059-1070. 10.1534/genetics.107.075804.PubMed CentralView ArticlePubMedGoogle Scholar
- Weedon MN, Lango H, Lindgren CM, Wallace C, Evans DM, Mangino M, Freathy RM, Perry JR, Stevens S, Hall AS, et al: Genome-wide association analysis identifies 20 loci that influence adult height. Nat Genet. 2008, 40: 575-583. 10.1038/ng.121.PubMed CentralView ArticlePubMedGoogle Scholar
- Gudbjartsson DF, Walters GB, Thorleifsson G, Stefansson H, Halldorsson BV, Zusmanovich P, Sulem P, Thorlacius S, Gylfason A, Steinberg S, et al: Many sequence variants affecting diversity of adult human height. Nat Genet. 2008, 40: 609-615. 10.1038/ng.122.View ArticlePubMedGoogle Scholar
- Weikard R, Goldammer T, Laurent P, Womack JE, Kuehn C: A gene-based high-resolution comparative radiation hybrid map as a framework for genome sequence assembly of a bovine chromosome 6 region associated with QTL for growth, body composition, and milk performance traits. BMC Genomics. 2006, 7: 53-10.1186/1471-2164-7-53.PubMed CentralView ArticlePubMedGoogle Scholar
- Nilsen H, Hayes B, Berg PR, Roseth A, Sundsaasen KK, Nilsen K, Lien S: Construction of a dense SNP map for bovine chromosome 6 to assist the assembly of the bovine genome sequence. Anim Genet. 2008, 39: 97-104. 10.1111/j.1365-2052.2007.01686.x.View ArticlePubMedGoogle Scholar
- Human (Homo sapiens) Genome Browser Gateway. [http://genome.ucsc.edu/cgi-bin/hgGateway]
- Sutani T, Yuasa T, Tomonaga T, Dohmae N, Takio K, Yanagida M: Fission yeast condensin complex: essential roles of non-SMC subunits for condensation and Cdc2 phosphorylation of Cut3/SMC4. Genes Dev. 1999, 13: 2271-2283. 10.1101/gad.13.17.2271.PubMed CentralView ArticlePubMedGoogle Scholar
- Dej KJ, Ahn C, Orr-Weaver TL: Mutations in the Drosophila condensin subunit dCAP-G: defining the role of condensin for chromosome condensation in mitosis and gene expression in interphase. Genetics. 2004, 168: 895-906. 10.1534/genetics.104.030908.PubMed CentralView ArticlePubMedGoogle Scholar
- Xing H, Wilkerson DC, Mayhew CN, Lubert EJ, Skaggs HS, Goodson ML, Hong Y, Park-Sarge OK, Sarge KD: Mechanism of hsp70i gene bookmarking. Science. 2005, 307: 421-423. 10.1126/science.1106478.View ArticlePubMedGoogle Scholar
- Ryu B, Kim DS, Deluca AM, Alani RM: Comprehensive expression profiling of tumor cell lines identifies molecular signatures of melanoma progression. PLoS ONE. 2007, 2: e594-10.1371/journal.pone.0000594.PubMed CentralView ArticlePubMedGoogle Scholar
- Ihara N, Takasuga A, Mizoshita K, Takeda H, Sugimoto M, Mizoguchi Y, Hirano T, Itoh T, Watanabe T, Reed KM, et al: 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
- Knott SA, Elsen JM, Haley CS: Methods for multiple-marker mapping of quantitative trait loci in half-sib populations. Theor Appl Genet. 1996, 93: 71-80. 10.1007/BF00225729.View ArticlePubMedGoogle Scholar
- Seaton G, Haley CS, Knott SA, Kearsey M, Visscher PM: QTL Express: mapping quantitative trait loci in simple and complex pedigrees. Bioinformatics. 2002, 18: 339-340. 10.1093/bioinformatics/18.2.339.View ArticlePubMedGoogle Scholar
- Knott SA, Marklund L, Haley CS, Andersson K, Davies W, Ellegren H, Fredholm M, Hansson I, Hoyheim B, Lundström K, et al: Multiple marker mapping of quantitative trait loci in a cross between outbred wild boar and large white pigs. Genetics. 1998, 149: 1069-1080.PubMed CentralPubMedGoogle Scholar
- Churchill GA, Doerge RW: Empirical threshold values for quantitative trait mapping. Genetics. 1994, 138: 963-971.PubMed CentralPubMedGoogle Scholar
- Visscher PM, Thompson R, Haley CS: Confidence intervals in QTL mapping by bootstrapping. Genetics. 1996, 143: 1013-1020.PubMed CentralPubMedGoogle Scholar
- Primer 3. [http://frodo.wi.mit.edu/cgi-bin/primer3/primer3_www.cgi]
- A BAC fingerprint map of the bovine genome. [http://www.bcgsc.ca/platform/mapping/bovine]
- NCBI/BLAST Home. [http://blast.ncbi.nlm.nih.gov/Blast.cgi]
- RPCI – 42 Male Bovine BAC Library. [http://bacpac.chori.org/mbovine42.htm]
- CHORI-240 Bovine BAC Library (Male). [http://bacpac.chori.org/bovine240.htm]
- Bovine Genome Project. [http://www.hgsc.bcm.tmc.edu/projects/bovine/]
- Dempster A, Laird N, Rubin D: Maximum likelihood estimation from incomplete data via the EM algorithm. J R Stat Soc. 1977, 39: 1-38.Google 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.