Heifers with positive genetic merit for fertility traits reach puberty earlier and have a greater pregnancy rate than heifers with negative genetic merit for fertility traits.

This study investigated the hypothesis that dairy heifers divergent in genetic merit for fertility traits differ in the age of puberty and reproductive performance. New Zealand's fertility breeding value (FertBV) is the proportion of a sire's daughters expected to calve in the first 42 d of the seasonal calving period. We used the New Zealand national dairy database to identify and select Holstein-Friesian dams with either positive (POS, +5 FertBV, n = 1,334) or negative FertBV (NEG, -5% FertBV, n = 1,662) for insemination with semen from POS or NEG FertBV sires, respectively. The resulting POS and NEG heifers were predicted to have a difference in average FertBV of 10 percentage points. We enrolled 640 heifer calves (POS, n = 324; NEG, n = 316) at 9 d ± 5.4 d (± standard deviation; SD) for the POS calves and 8 d ± 4.4 d old for the NEG calves. Of these, 275 POS and 248 NEG heifers were DNA parent verified and retained for further study. The average FertBV was +5.0% (SD = 0.74) and -5.1% (SD = 1.36) for POS and NEG groups, respectively. Heifers were reared at 2 successive facilities as follows: (1) calf rearing (enrollment to ∼13 wk of age) and (2) grazier, after 13 wk until 22 mo of age. All heifers wore a collar with an activity sensor to monitor estrus events starting at 8 mo of age, and we collected weekly blood samples when individual heifers reached 190 kg of body weight (BW) to measure plasma progesterone concentrations. Puberty was characterized by plasma progesterone concentrations >1 ng/mL in at least 2 of 3 successive weeks. Date of puberty was defined when the first of these samples was >1 ng/mL. Heifers were seasonally bred for 98 d starting at ∼14 mo of age. Transrectal ultrasound was used to confirm pregnancy and combined with activity data to estimate breeding and pregnancy dates. We measured BW every 2 wk, and body condition and stature at 6, 9, 12, and 15 mo of age. The significant FertBV by day interaction for BW was such that the NEG heifers had increasingly greater BW with age. This difference was mirrored with the significant FertBV by month interaction for average daily gain, with the NEG heifers having a greater average daily gain between 9 and 18 mo of age. There was no difference in heifer stature between the POS and NEG heifers. The POS heifers were younger and lighter at puberty, and were at a lesser mature BW, compared with the NEG heifers. As a result, 94 ± 1.6% of the POS and 82 ± 3.2% of the NEG heifers had reached puberty at the start of breeding. The POS heifers were 20% and 11% more likely to be pregnant after 21 d and 42 d of breeding than NEG heifers (relative risk = 1.20, 95% confidence interval of 1.03-1.34; relative risk = 1.11, 95% confidence interval of 1.01-1.16). Results from this experiment support an association between extremes in genetic merit for fertility base on cow traits and heifer reproduction. Our results indicate that heifer puberty and pregnancy rates are affected by genetic merit for fertility traits, and these may be useful phenotypes for genetic selection.

Using genotyping-by-sequencing to predict gender in animals.

Gender assignment errors are common in some animal species and lead to inaccuracies in downstream analyses. Procedures for detecting gender misassignment are available for array-based SNP data but are still being developed for genotyping-by-sequencing (GBS) data. In this study, we describe a method for using GBS data to predict gender using X and Y chromosomal SNPs. From a set of 1286 X chromosomal and 23 Y chromosomal deer (Cervus sp.) SNPs discovered from GBS sequence reads, a prediction model was built using a training dataset of 422 Red deer and validated using a test dataset of 868 Red deer and Wapiti deer. Prediction was based on the proportion of heterozygous genotypes on the X chromosome and the proportion of non-missing genotypes on the Y chromosome observed in each individual. The concordance between recorded gender and predicted gender was 98.6% in the training dataset and 99.3% in the test dataset. The model identified five individuals across both datasets with incorrect recorded gender and was unable to predict gender for another five individuals. Overall, our method predicted gender with a high degree of accuracy and could be used for quality control in gender assignment datasets or for assigning gender when unrecorded, provided a suitable reference genome is available.

Identifying chromosomal selection-sweep regions in facial eczema selection-line animals using an ovine 50K-SNP array.

Facial eczema (FE) is a hepato-mycotoxicosis found mainly in New Zealand sheep and cattle. When genetics was found to be a factor in FE susceptibility, resistant and susceptible selection lines of Romney sheep were established to enable further investigations of this disease trait. Using the Illumina OvineSNP50 BeadChip, we conducted a selection-sweep experiment on these FE genetic lines. Two analytical methods were used to detect selection signals, namely the Peddrift test (Dodds & McEwan, 1997) and fixation index FST (Weir & Hill, 2002). Of 50 975 single nucleotide polymorphism (SNP) markers tested, there were three that showed highly significant allele frequency differences between the resistant and susceptible animals (Peddrift nominal P < 0.000001). These SNP loci are located on chromosomes OAR1, OAR11 and OAR12 that coincide precisely with the three highest genomic FST peaks. In addition, there are nine less significant Peddrift SNPs (nominal P ≤ 0.000009) on OAR6 (n = 2), OAR9 (n = 2), OAR12, OAR19 (n = 2), OAR24 and OAR26. In smoothed FST (five-SNP moving average) plots, the five most prominent peaks are on OAR1, OAR6, OAR7, OAR13 and OAR19. Although these smoothed FST peaks do not coincide with the three most significant Peddrift SNP loci, two (on OAR6 and OAR19) overlap with the set of less significant Peddrift SNPs above. Of these 12 Peddrift SNPs and five smoothed FST regions, none is close to the FE candidate genes catalase and ABCG2; however, two on OAR1 and one on OAR13 fall within suggestive quantitative trait locus regions identified in a previous genome screen experiment. The present studies indicated that there are at least eight genomic regions that underwent a selection sweep in the FE lines.

Towards genomic selection for facial eczema disease tolerance in the New Zealand sheep industry.

Pithomycotoxicosis, more commonly known as facial eczema (FE), is a liver disease that occurs predominantly in New Zealand because of its toxigenic Pithomyces chartarum strains. The first reported case was in sheep in 1887. Since the 1930s, a number of studies have been conducted in an attempt to mitigate the problems FE has on the sheep and dairy industries. The research in these studies included work on fungicide and biological control of the saprophytic fungus, use of different pasture plants to inhibit fungal growth, stock management with respect to pasture fungal spore counts and the use of zinc prophylaxis on animals. The finding that there was a genetic basis in FE sensitivity in sheep prompted research for a genetic approach to mitigation in the form of a diagnostic DNA test for susceptibility to the disease. Recently, we have used the Illumina OvineSNP50 BeadChip to develop a genome-enabled prediction approach to screen for FE-tolerant sheep. Our current best genomic prediction for FE is for the Romney breed and has an accuracy of 0.38. This prediction accuracy is not as high as the individual accuracy gained by an artificial challenge test (0.64). However, it has the advantage of being a non-invasive test and can be provided as part of genomic testing for other traits at minimal cost.

The ovine Booroola fecundity gene (FecB) is linked to markers from a region of human chromosome 4q.

The autosomal Booroola fecundity gene (FecB) mutation in sheep increases ovulation rate and litter size, with associated effects on ovarian physiology and hormone profiles. Analysis of segregation in twelve families (379 female progeny) identified linkage between the mutation, two microsatellite markers (OarAE101 and OarHH55, Zmax > 9.0) and epidermal growth factor (EGF) from human chromosome 4q25 (Zmax > 3.0). The marker OarAE101 was linked to secreted phosphoprotein 1 (SPP1, which maps to chromosome 4q21-23 in man) in the test pedigrees and independent families (Zmax > 9.7). The identification of linkage between the FecB mutation and markers from human chromosome 4q is an important step towards further understanding the control of ovulation rates in mammals.

Mutations in an oocyte-derived growth factor gene (BMP15) cause increased ovulation rate and infertility in a dosage-sensitive manner.

Multiple ovulations are uncommon in humans, cattle and many breeds of sheep. Pituitary gonadotrophins and as yet unidentified ovarian factors precisely regulate follicular development so that, normally, only one follicle is selected to ovulate. The Inverdale (FecXI) sheep, however, carries a naturally occurring X-linked mutation that causes increased ovulation rate and twin and triplet births in heterozygotes (FecXI/FecX+; ref. 1), but primary ovarian failure in homozygotes (FecXI/FecXI; ref. 2). Germ-cell development, formation of the follicle and the earliest stages of follicular growth are normal in FecXI/FecXI sheep, but follicular development beyond the primary stage is impaired. A second family unrelated to the Inverdale sheep also has the same X-linked phenotype (Hanna, FecXH). Crossing FecXI with FecXH animals produces FecXI/FecXH infertile females phenotypically indistinguishable from FecXI/FecXI females. We report here that the FecXI locus maps to an orthologous chromosomal region syntenic to human Xp11.2-11.4, which contains BMP15, encoding bone morphogenetic protein 15 (also known as growth differentiation factor 9B (GDF9B)). Whereas BMP15 is a member of the transforming growth factor beta (TGFbeta) superfamily and is specifically expressed in oocytes, its function is unknown. We show that independent germline point mutations exist in FecXI and FecXH carriers. These findings establish that BMP15 is essential for female fertility and that natural mutations in an ovary-derived factor can cause both increased ovulation rate and infertility phenotypes in a dosage-sensitive manner.

Potential of in-plant intramuscular fat predictions to enable sheep breeders to incorporate consumer preferences in breeding programmes.

The inclusion of eating quality traits in sheep genetic improvement programmes is desirable. Intramuscular fat (IMF) plays a key role in ensuring consumer satisfaction when eating lamb, but genetic progress for IMF is constrained by a lack of routine data collection. This study investigated the potential for IMF predictor traits to substitute for measured IMF in genetic improvement programmes. Carcass and predicted IMF (near-infrared estimated IMF and marbling score) data were available on 10,113 New Zealand lambs, 1678 of which also had measured chemical IMF on a slice of M. longissimus lumborum on which the predictions of IMF had been made. Genetic antagonisms were observed between carcass lean traits and IMF. The genetic correlation between the predictors and measured IMF approached one, indicating that predictors of IMF can be used in genetic improvement programmes. Through using selection indexes, simultaneous increases in IMF and the existing terminal selection index are possible, provided all traits are measured. This study highlights the importance and potential of predicted IMF to achieve genetic improvement in traits of importance to consumers.

Can we have our steak and eat it: The impact of breeding for lowered environmental impact on yield and meat quality in sheep.

Global agreements in place to reduce methane emissions in livestock are a potential threat to food security. Successful but independent breeding strategies for improved production and lower methane are in place. The unanswered questions are whether these strategies can be combined and how they impact one another, physically and economically. The New Zealand economy is largely dependent on pastoral agriculture from grazing ruminants. The sheep industry produces ∼20 million lamb carcasses for export each year primarily from grass. Methane emitted from the fermentation of forage by grazing ruminants accounts for one-third of all New Zealand's greenhouse gas emissions. Here, we use sheep selection lines bred for divergent methane production and large numbers of their relatives to determine the genetic and phenotypic correlations between enteric methane emissions, carcass yield, and meat quality. The primary objectives were to determine whether previously shown physiological differences between methane selection lines (differing by ∼12% in methane) result in a negative impact on meat production and quality by measuring close relatives. The results show no negative effects of breeding for lowered methane on meat and carcass quality. Gross methane emissions were highly correlated with liveweight and measures of carcass weight and negatively correlated with dressing-out percentage and fat yield (GR). Trends were similar but not significant for methane yield (g CH4/kg DMI). Preliminary evidence, to date, shows that breeding for low methane may result in animals with higher lean yields that are economically favorable even before carbon costs and environmental benefits are taken into account. These benefits were seen in animals measured for methane on fixed intakes and require validation on intakes that are allowed to vary.

Heterozygous Inverdale ewes show increased ovulation rate sensitivity to pre-mating nutrition.

This study aimed to determine whether ewes heterozygous (I+) for the Inverdale mutation of the bone morphogenetic protein-15 (BMP15) gene with high natural ovulation rate (OR) show similar sensitivity to nutritional manipulation as non-carriers (++). Increasing pre-mating nutrition results in OR increases in sheep, but whether this effect occurs in ewes with naturally high OR is unknown. Over 2 years, I+ or ++ ewes were given high (ad libitum) or control (maintenance) pasture allowances for 6 weeks prior to mating at a synchronised oestrus, with OR measured 8 days later. The high group increased in weight compared with controls (+5.84kg; P<0.01), accompanied by increased OR (+19%; P<0.01). As well as having higher OR (+45%; P<0.01), I+ ewes responded to increased feed with a larger proportional increase in OR (+27%; P<0.01) compared with the response in ++ ewes (+11%; P<0.05), suggesting an interaction between BMP15 levels and nutritional signals in the follicle to control OR. Although litter size increases only tended to significance (+12%; P=0.06), extra feed resulted in over 50% of I+ ewes giving birth to more than three lambs, compared with 20-31% of I+ ewes on maintenance rations. This information can guide feed management of prolific Inverdale ewes prior to breeding.

A linkage map of sheep chromosome X (OARX) aligned to human chromosome X (HSAX).

We have constructed a genetic linkage map of the sheep X chromosome (OARX) containing 22 new gene loci from across the human X chromosome (HSAX). The female OARX linkage map has a total length of 152.6 cM with average gene spacing of 5.5 cM. Comparison with HSAX confirms one previously reported major breakpoint and inversion, and other minor rearrangements between OARX and HSAX. Comparison of the linkage map with sheep sequence data OAR 1.0 reveals a different arrangement of markers on the q arm, which may more accurately reflect the genuine arrangement of this region.

A genome-screen experiment to detect quantitative trait loci affecting resistance to facial eczema disease in sheep.

Facial eczema (FE) is a secondary photosensitization disease arising from liver cirrhosis caused by the mycotoxin sporidesmin. The disease affects sheep, cattle, deer and goats, and costs the New Zealand sheep industry alone an estimated NZ$63M annually. A long-term sustainable solution to this century-old FE problem is to breed for disease-resistant animals by marker-assisted selection. As a step towards finding a diagnostic DNA test for FE sensitivity, we have conducted a genome-scan experiment to screen for quantitative trait loci (QTL) affecting this trait in Romney sheep. Four F(1) sires, obtained from reciprocal matings of FE resistant and susceptible selection-line animals, were used to generate four outcross families. The resulting half-sib progeny were artificially challenged with sporidesmin to phenotype their FE traits measured in terms of their serum levels of liver-specific enzymes, namely gamma-glutamyl transferase and glutamate dehydrogenase. In a primary screen using selective genotyping on extreme progeny of each family, a total of 244 DNA markers uniformly distributed over all 26 ovine autosomes (with an autosomal genome coverage of 79-91%) were tested for linkage to the FE traits. Data were analysed using Haley-Knott regression. The primary screen detected one significant and one suggestive QTL on chromosomes 3 and 8 respectively. Both the significant and suggestive QTL were followed up in a secondary screen where all progeny were genotyped and analysed; the QTL on chromosome 3 was significant in this analysis.

Quantitative trait loci for live animal and carcass composition traits in Jersey and Limousin back-cross cattle finished on pasture or feedlot.

A quantitative trait locus (QTL) study was carried out in two countries, recording live animal and carcass composition traits. Back-cross calves (385 heifers and 398 steers) were generated, with Jersey and Limousin breed backgrounds. The New Zealand cattle were reared on pasture to carcass weights averaging 229 kg, whilst the Australian cattle were reared on grass and finished on grain (for at least 180 days) to carcass weights averaging 335 kg. From 11 live animal traits and 31 carcass composition traits respectively, 5 and 22 QTL were detected in combined-sire analyses, which were significant (P < 0.05) on a genome-wise basis. Fourteen significant traits for carcass composition QTL were on chromosome 2 and these were traits associated with muscling and fatness. This chromosome carried a variant myostatin allele (F94L), segregating from the Limousin ancestry. Despite very different cattle management systems between the two countries, the two populations had a large number of QTL in common. Of the 18 traits which were common to both countries, and which had significant QTL at the genome-wise level, eight were significant in both countries.

Generation of a preliminary bovine gene atlas, using expression clustering to annotate gene function.

Genes whose products function in a common biological process are often co-regulated. When regulation occurs at the transcriptional level, co-expressed genes can be detected globally by expression arrays or by sequencing non-normalized cDNA libraries. We examined bovine gene expression in 27 tissues using non-normalized cDNA library sequencing. Contigs were generated from expressed sequence tags whose sequences overlapped. Contigs containing a minimum of five expressed sequence tags were ordered via a hierarchical clustering process, where the distance between the contigs represents their expression pattern similarity across tissues. Gene ontology terms associated with the genes in each cluster showed that co-clustered genes encoded proteins involved in a common biological process. This process can be used to annotate genes of unknown function in the cluster. Gene expression was compared between bovine and human tissues; there were significant correlations between species for each tissue, with the exception of thyroid and placenta. Tissues were also clustered based on the genes they express; tissues with similar physiological functions clustered closely. Based on this information, we generated the first preliminary gene atlas of the bovine genome. Genes with similar expression patterns were clustered, and genes with a common function co-clustered. This method can be used to annotate genes of unknown function in the bovine genome.

Heritability of ram mating success in multi-sire breeding situations.

Multi-sire mating of a mob of ewes is commonly used in commercial sheep production systems. However, ram mating success (defined as the number of lambs sired by an individual) can vary between rams in the mating group. If this trait was repeatable and heritable, selection of rams capable of siring larger numbers of lambs could reduce the number of rams required for mating and ultimately lead to increased genetic gain. However, genetic correlations with other productive traits, such as growth and female fertility, could influence the potential for ram mating success to be used as a selection trait. In order to investigate this trait, parentage records (including accuracy of sire assignment) from 15 commercial ram breeding flocks of various breeds were utilised to examine the repeatability and heritability of ram mating success in multi-sire mating groups. In addition, genetic and phenotypic correlations with growth and female fertility traits were estimated using ASReml. The final model used for the ram mating success traits included age of the ram and mating group as fixed effects. Older rams (3+years old) had 15% to 20% greater mating success than younger rams (1 or 2 years of age). Increasing the stringency of the criteria for inclusion of both an individual lamb, based on accuracy of sire assignment, or a whole mating group, based on how many lambs had an assigned sire, increased repeatability and heritability estimates of the ram mating success traits examined. With the most stringent criteria employed, where assignment of sire accuracy was >0.95 and the total number of lambs in the progeny group that failed to have a sire assigned was<0.05, repeatability and heritability for loge(number of lambs) was 0.40±0.09 and 0.26±0.12, respectively. For proportion of lambs sired, repeatability and heritability were both 0.30±0.09. The two ram mating traits (loge(nlamb) and proportion) were highly correlated, both phenotypically and genetically (0.88±0.01 and 0.94±0.06, respectively). Both phenotypic and genetic correlations between ram mating success and growth and other female fertility traits were low and non-significant. In conclusion, there is scope to select rams capable of producing high numbers of progeny and thus increase selection pressure on rams to increase genetic gain.

Identifying genes for intestinal nematode resistance using transcriptional profiling.

Gene expression was compared between resistant and susceptible Perendale lambs that had either never been exposed to gastrointestinal nematode challenge (had a naïve immune system with respect to parasites) or had been naturally challenged on pasture with nematodes. Only a small number of genes were differentially expressed between the naive resistant and susceptible animals, but many genes were differentially expressed between the resistant and susceptible challenged animals. The differentially expressed genes were involved in a variety of biological processes, most notably the immune response, the stress response and gene regulation via chromatin remodelling. The transcriptional profiling experiments also detected gene expression differences in the Ovar-DQA1 gene between resistant and susceptible challenged animals. A null allele of this gene was demonstrated to be associated with susceptibility to gastrointestinal parasites in some, but not all populations. This allele is not thought to be causal for susceptibility.

Comparative mapping of sheep chromosome 2q.

Sheep chromosome 2q (OAR2q), which is homologous with human chromosome 2q (HSA2q), and cattle chromosome 2 (BTA2), is known to contain several loci contributing to carcass traits. However, the chromosomal rearrangements differentiating these chromosomes among the three species have not yet been determined and thus precise correspondences between the locations of sheep and human genes are not known. Twenty-six genes from HSA2q (2q21.1-->2q36) have been assigned to OAR2q by genetic linkage mapping to refine this area of the sheep genome. Seventy-six genes were initially selected from HSA2q. Sixty-eight percent of the PCR primer sets designed for these genes amplified successfully in sheep, and 34% amplified polymorphic products. Part of the proximal arm of OAR2q was found to be inverted compared with HSA2q. The breakpoint has been localised near the growth differentiation factor 8 gene (GDF8), spanning 380 kb between the positions of the hypothetical protein (FLJ20160) (HSA2:191008944-191075046) and glutaminase (GLS) (HSA2:191453847-191538510) (Build36.1).

Characterization of an X-chromosomal non-mosaic monosomy (59, X0) dairy heifer detected using routinely available single nucleotide polymorphism genotype data.

Evidence exists from a range of species on the impact of karyotype abnormalities on reproductive performance. Despite this, cytogenetic analyses of cattle, especially females, are not routinely undertaken. Genome-wide single nucleotide polymorphism (SNP) genotype data are now, however, routinely being generated in many species globally at a relatively low cost. The objective of the present study was to evaluate the potential of routinely available SNP genotype data to identify sex-chromosome aberrations using X chromosome monosomy 59,X0 as a case study for illustration. A single 2.5-yr old Holstein-Friesian heifer was detected with a mean allelic intensity of SNP on the X chromosome almost 17 standard deviations less than the mean of other genotyped females ( = 103,326). Following cytogenetic analysis (10 replicates by karyotyping and a further 140 by FISH), the female was deduced to be a non-mosaic 59,X0. The female had never produced a calf and, although gross examination revealed no physical abnormalities, she was smaller in size than expected based on her breed and age. Given the age of the animal at slaughter, the uterus and uterine tubes appeared immature and inactive. The oviduct appeared normal while the single ovary present contained a markedly reduced number of follicles. There was, however, some evidence of prior ovulation and formation of corpora lutea. The approach proposed in the present study to identify allosome aneuploidy from routinely available genotype data can be used to screen for such abnormalities at no additional cost to the breeder or producer.

A putative autosomal gene increasing ovulation rate in Romney sheep.

Ovulation rates were measured in 547 progeny of 24 rams in a Romney flock with a long history of high prolificacy. These sheep were from the same family line and the distribution of ovulation rates suggests the presence of a segregating major gene (FecW) that increases prolificacy. The phenotype differs from those previously described for major genes affecting prolificacy in sheep. The putative gene shows autosomal inheritance and one copy increases ovulation rate by 0.8-1.0 eggs per ewe ovulating. To date, we have found no evidence of infertility among putative homozygous ewes, as described in some autosomal major genes for prolificacy.

Linkage disequilibrium in domestic sheep.

The last decade has seen a dramatic increase in the number of livestock QTL mapping studies. The next challenge awaiting livestock geneticists is to determine the actual genes responsible for variation of economically important traits. With the advent of high density single nucleotide polymorphism (SNP) maps, it may be possible to fine map genes by exploiting linkage disequilibrium between genes of interest and adjacent markers. However, the extent of linkage disequilibrium (LD) is generally unknown for livestock populations. In this article microsatellite genotype data are used to assess the extent of LD in two populations of domestic sheep. High levels of LD were found to extend for tens of centimorgans and declined as a function of marker distance. However, LD was also frequently observed between unlinked markers. The prospects for LD mapping in livestock appear encouraging provided that type I error can be minimized. Properties of the multiallelic LD coefficient D' were also explored. D' was found to be significantly related to marker heterozygosity, although the relationship did not appear to unduly influence the overall conclusions. Of potentially greater concern was the observation that D' may be skewed when rare alleles are present. It is recommended that the statistical significance of LD is used in conjunction with coefficients such as D' to determine the true extent of LD.

A new gene mapping resource: interspecies hybrids between Père David's deer (Elaphurus davidianus) and red deer (Cervus elaphus).

Three male F1 hybrids between Père David's deer and red deer were mated to red deer to produce 143 backcross calves. The pedigrees are a rare example of a fertile hybrid between evolutionarily divergent species. We examined the use of these families for genetic mapping of evolutionarily conserved (Type I) loci by testing for genetic linkage between five species-specific protein variants and 12 conserved DNA probes. Two probes were homologous, and the remainder syntenic, to the protein coding loci in cattle or humans. Using six restriction enzymes, each DNA probe detected one or more restriction fragments specific to Père David's deer. Linkage analyses among the species-specific variants placed the loci into four linkage groups within which linkage between adjacent loci and gene order was supported by a LOD > 3. The linkage groups were (HPX, HBB)-FSHB-ACP2, LDHA-CD5-IGF2, BMP3-(GC, ALB)-(KIT, PDGFRA) and LDLR-C3-FGF1. Southern and protein analysis of LDHA and ALB provided identical segregation data. These linkage groups were consistent with the cattle gene map and provide new information for comparing the gene maps of ruminants, humans and mice. The deer hybrids are an important new resource that can contribute to the comparative analysis of the mammalian genome.