Construction and characterization of a new bovine bacterial artificial chromosome library with 10 genome-equivalent coverage

Internal Duplications of DH, JH, and C Region Genes Create an Unusual IgH Gene Locus in Cattle

It has been suspected for many years that cattle possess two functional IgH gene loci, located on Bos taurus autosome (BTA) 21 and BTA11, respectively. In this study, based on fluorescence in situ hybridization and additional experiments, we showed that all functional bovine IgH genes were located on BTA21, and only a truncated μCH2 exon was present on BTA11. By sequencing of seven bacterial artificial chromosome clones screened from a Hostein cow bacterial artificial chromosome library, we generated a 678-kb continuous genomic sequence covering the bovine IGHV, IGHD, IGHJ, and IGHC genes, which are organized as IGHVn-IGHDn-IGHJn-IGHM1-(IGHDP-IGHV3-IGHDn)3-IGHJn-IGHM2-IGHD-IGHG3-IGHG1-IGHG2-IGHE-IGHA. Although both of two functional IGHM genes, IGHM1 and IGHM2, can be expressed via independent VDJ recombinations, the IGHM2 can also be expressed through class switch recombination. Likely because more IGHD segments can be involved in the expression of IGHM2, the IGHM2 gene was shown to be dominantly expressed in most tissues throughout different developmental stages. Based on the length and identity of the coding sequence, the 23 IGHD segments identified in the locus could be divided into nine subgroups (termed IGHD1 to IGHD9). Except two members of IGHD9 (14 nt in size), all other functional IGHD segments are longer than 30 nt, with the IGHD8 gene (149 bp) to be the longest. These remarkably long germline IGHD segments play a pivotal role in generating the exceptionally great H chain CDR 3 length variability in cattle.

First steps: bovine genomics in historical perspective

The vision of Morris Soller was instrumental in launching the field of bovine genomics. This study is a review of the early years of bovine gene mapping leading up to the sequencing and assembly of the bovine genome in 2009. A historical perspective of parasexual, linkage and physical mapping is provided with a focus on the contribution of these maps to the eventual assignment and orientation of genes and sequence to cattle chromosomes.

Construction and preliminary characterization of a river buffalo bacterial artificial chromosome library

River buffalo genome analyses have advanced significantly in the last decade, and the genome sequence of Bubalus bubalis will be available shortly. Nonetheless, large-insert DNA library resources such as bacterial artificial chromosomes (BAC) are still required for validation and accurate assembly of the genome sequence. We constructed a river buffalo BAC library containing 52,224 clones with an average insert size of 97 kb, representing 1.7 × coverage of the genome. This genomic resource for river buffalo will facilitate further studies in this economically important species allowing for instance, whole genome physical mapping and isolation of genes and gene clusters, contributing to the elucidation of gene organization and identification of regulatory elements.

Construction of a llama bacterial artificial chromosome library with approximately 9-fold genome equivalent coverage

The Ilama is an important agricultural livestock in much of South America. The llama is increasing in popularity in the United States as a companion animal. Little work has been done to improve llama production using modern technology. A paucity of information is available regarding the llama genome. We report the construction of a llama bacterial artificial chromosome (BAC) library of about 196,224 clones in the vector pECBAC1. Using flow cytometry and bovine, human, mouse, and chicken as controls, we determined the llama genome size to be 2.4 × 10⁹ bp. The average insert size of the library is 137.8 kb corresponding to approximately 9-fold genome coverage. Further studies are needed to further characterize the library and llama genome. We anticipate that this new library will help facilitate future genomic studies in the llama.

A 4 Mb high resolution BAC contig on bovine chromosome 1q12 and comparative analysis with human chromosome 21q22

The bovine RPCI-42 BAC library was screened to construct a sequence-ready ~4 Mb single contig of 92 BAC clones on BTA 1q12. The contig covers the region between the genes KRTAP8P1 and CLIC6. This genomic segment in cattle is of special interest as it contains the dominant gene responsible for the hornless or polled phenotype in cattle. The construction of the BAC contig was initiated by screening the bovine BAC library with heterologous cDNA probes derived from 12 human genes of the syntenic region on HSA 21q22. Contig building was facilitated by BAC end sequencing and chromosome walking. During the construction of the contig, 165 BAC end sequences and 109 single-copy STS markers were generated. For comparative mapping of 25 HSA 21q22 genes, genomic PCR primers were designed from bovine EST sequences and the gene-associated STSs mapped on the contig. Furthermore, bovine BAC end sequence comparisons against the human genome sequence revealed significant matches to HSA 21q22 and allowed the in silico mapping of two new genes in cattle. In total, 31 orthologues of human genes located on HSA 21q22 were directly mapped within the bovine BAC contig, of which 16 genes have been cloned and mapped for the first time in cattle. In contrast to the existing comparative bovine–human RH maps of this region, these results provide a better alignment and reveal a completely conserved gene order in this 4 Mb segment between cattle, human and mouse. The mapping of known polled linked BTA 1q12 microsatellite markers allowed the integration of the physical contig map with existing linkage maps of this region and also determined the exact order of these markers for the first time. Our physical map and transcript map may be useful for positional cloning of the putative polled gene in cattle. The nucleotide sequence data reported in this paper have been submitted to EMBL and have been assigned Accession Numbers AJ698510–AJ698674.

Development and characterization of a high density SNP genotyping assay for cattle

The success of genome-wide association (GWA) studies for the detection of sequence variation affecting complex traits in human has spurred interest in the use of large-scale high-density single nucleotide polymorphism (SNP) genotyping for the identification of quantitative trait loci (QTL) and for marker-assisted selection in model and agricultural species. A cost-effective and efficient approach for the development of a custom genotyping assay interrogating 54,001 SNP loci to support GWA applications in cattle is described. A novel algorithm for achieving a compressed inter-marker interval distribution proved remarkably successful, with median interval of 37 kb and maximum predicted gap of <350 kb. The assay was tested on a panel of 576 animals from 21 cattle breeds and six outgroup species and revealed that from 39,765 to 46,492 SNP are polymorphic within individual breeds (average minor allele frequency (MAF) ranging from 0.24 to 0.27). The assay also identified 79 putative copy number variants in cattle. Utility for GWA was demonstrated by localizing known variation for coat color and the presence/absence of horns to their correct genomic locations. The combination of SNP selection and the novel spacing algorithm allows an efficient approach for the development of high-density genotyping platforms in species having full or even moderate quality draft sequence. Aspects of the approach can be exploited in species which lack an available genome sequence. The BovineSNP50 assay described here is commercially available from Illumina and provides a robust platform for mapping disease genes and QTL in cattle.

A multiway analysis for identifying high integrity bovine BACs

BACKGROUND

In large genomics projects involving many different types of analyses of bacterial artificial chromosomes (BACs), such as fingerprinting, end sequencing (BES) and full BAC sequencing there are many opportunities for the identities of BACs to become confused. However, by comparing the results from the different analyses, inconsistencies can be identified and a set of high integrity BACs preferred for future research can be defined.

RESULTS

The location of each bovine BAC in the BAC fingerprint-based genome map and in the genome assembly were compared based on the reported BESs, and for a smaller number of BACs the full sequence. BACs with consistent positions in all three datasets, or if the full sequence was not available, for both the fingerprint map and BES-based alignments, were deemed to be correctly positioned. BACs with consistent BES-based and fingerprint-based locations, but with conflicting locations based on the fully sequenced BAC, appeared to have been misidentified during sequencing, and included a number of apparently swapped BACs. Inconsistencies between BES-based and fingerprint map positions identified thirty one plates from the CHORI-240 library that appear to have suffered substantial systematic problems during the end-sequencing of the BACs. No systematic problems were identified in the fingerprinting of the BACs. Analysis of BACs overlapping in the assembly identified a small overrepresentation of clones with substantial overlap in the library and a substantial enrichment of highly overlapping BACs on the same plate in the CHORI-240 library. More than half of these BACs appear to have been present as duplicates on the original BAC-library plates and thus should be avoided in subsequent projects.

CONCLUSION

Our analysis shows that approximately 95% of the bovine CHORI-240 library clones with both a BAC fingerprint and two BESs mapping to the genome in the expected orientations (approximately 27% of all BACs) have consistent locations in the BAC fingerprint map and the genome assembly. We have developed a broadly applicable methodology for checking the integrity of BAC-based datasets even where only incomplete and partially assembled genomic sequence is available.

Characterization of the bovine pseudoautosomal boundary: Documenting the evolutionary history of mammalian sex chromosomes

Here, we report the sequence characterization of the bovine pseudoautosomal boundary (PAB) and its neighborhood. We demonstrate that it maps to the 5' end of the GPR143 gene, which has concomitantly lost upstream noncoding exons on the Y chromosome. We show that the bovine PAB was created approximately 20.7 million years ago by illegitimate intrachromatid recombination between inverted, ruminant-specific Bov-tA repeats. Accordingly, we demonstrate that cattle share their PAB with all other examined ruminants including sheep, but not with cetaceans or more distantly related mammals. We provide evidence that, since its creation, the ancestral ruminant PAB has been displaced by attrition, which occurs at variable rates in different species, and that it is capable of retreat by attrition erasure. We have estimated the ratio of male to female mutation rates in the Bovidae family as approximately 1.7, and we provide evidence that the mutation rate is higher in the recombining pseudoautosomal region than in the adjacent, nonrecombining gonosome-specific sequences.

Characterization of a 320-kb region containing the HEXA gene on bovine chromosome 10 and analysis of its association with BSE susceptibility

Bovine spongiform encephalopathy (BSE) belongs to a group of neurodegenerative diseases known as transmissible prion diseases. Recently, variants in the promoter region of the prion protein (PRNP) gene have been shown to have a considerable effect on the susceptibility to BSE. However, a previous genome scan revealed other putative BSE-susceptibility loci. Here, we analysed such a region on BTA10, which contains the functional candidate gene HEXA. Three hundred and twenty kilobases that, besides HEXA, also contain ARIH1, BRUNOL6 and PARP6 were characterized and screened for polymorphisms. Genotyping of 38 SNPs in Holstein-Friesian animals from the UK (350 diseased and 270 controls) revealed two intronic SNPs that were associated with BSE incidence, with experiment-wise P-values of 3.5 x 10(-3) and 7.7 x 10(-3) respectively. Both SNPs were in strong linkage disequilibrium and the rare alleles had a protective effect. These alleles were contained in a haplotype dubbed 'UK-protective' that was significantly overrepresented in the controls with a permuted P-value of 2 x 10(-3). An association study in German Holstein animals (73 diseased and 627 controls) revealed an opposite effect of the 'UK-protective' haplotype in this population, i.e. it was overrepresented in the diseased animals, although not significant after correction for multiple testing. These findings indicate a causal variant for BSE susceptibility on BTA10 in linkage disequilibrium with the markers studied. Candidate gene analyses of the surrounding region and additional association studies will help to clarify the origin of the protective effects and to identify causal variants for BSE susceptibility on BTA10.

A physical map of the bovine genome

A new physical map of the bovine genome has been constructed by integrating data from genetic and radiation hybrid maps, and a new bovine BAC map, with the bovine genome draft assembly.

Background

Cattle are important agriculturally and relevant as a model organism. Previously described genetic and radiation hybrid (RH) maps of the bovine genome have been used to identify genomic regions and genes affecting specific traits. Application of these maps to identify influential genetic polymorphisms will be enhanced by integration with each other and with bacterial artificial chromosome (BAC) libraries. The BAC libraries and clone maps are essential for the hybrid clone-by-clone/whole-genome shotgun sequencing approach taken by the bovine genome sequencing project.

Results

A bovine BAC map was constructed with HindIII restriction digest fragments of 290,797 BAC clones from animals of three different breeds. Comparative mapping of 422,522 BAC end sequences assisted with BAC map ordering and assembly. Genotypes and pedigree from two genetic maps and marker scores from three whole-genome RH panels were consolidated on a 17,254-marker composite map. Sequence similarity allowed integrating the BAC and composite maps with the bovine draft assembly (Btau3.1), establishing a comprehensive resource describing the bovine genome. Agreement between the marker and BAC maps and the draft assembly is high, although discrepancies exist. The composite and BAC maps are more similar than either is to the draft assembly.

Conclusion

Further refinement of the maps and greater integration into the genome assembly process may contribute to a high quality assembly. The maps provide resources to associate phenotypic variation with underlying genomic variation, and are crucial resources for understanding the biology underpinning this important ruminant species so closely associated with humans.

The role of the bovine growth hormone receptor and prolactin receptor genes in milk, fat and protein production in Finnish Ayrshire dairy cattle

We herein report new evidence that the QTL effect on chromosome 20 in Finnish Ayrshire can be explained by variation in two distinct genes, growth hormone receptor (GHR) and prolactin receptor (PRLR). In a previous study in Holstein-Friesian dairy cattle an F279Y polymorphism in the transmembrane domain of GHR was found to be associated with an effect on milk yield and composition. The result of our multimarker regression analysis suggests that in Finnish Ayrshire two QTL segregate on the chromosomal region including GHR and PRLR. By sequencing the coding sequences of GHR and PRLR and the sequence of three GHR promoters from the pooled samples of individuals of known QTL genotype, we identified two substitutions that were associated with milk production traits: the previously reported F-to-Y substitution in the transmembrane domain of GHR and an S-to-N substitution in the signal peptide of PRLR. The results provide strong evidence that the effect of PRLR S18N polymorphism is distinct from the GHR F279Y effect. In particular, the GHR F279Y has the highest influence on protein percentage and fat percentage while PRLR S18N markedly influences protein and fat yield. Furthermore, an interaction between the two loci is suggested.

Evidence that bovine forebrain embryonic zinc finger-like gene influences immune response associated with mastitis resistance

Mastitis, a mammary gland inflammation in response to bacterial infection, is a major problem in the dairy industry. We found that cows susceptible to mastitis have a three-base insertion in a glycine-coding stretch of the gene for forebrain embryonic zinc finger-like (FEZL), a transcription factor with a role in neuronal development. Mastitis induces FEZL expression in mammary glands, and induced FEZL promotes expression of the axon-attracting molecule semaphorin 5A (SEMA5A) through a GCAG sequence. FEZL also induces SEMA5A expression in susceptible cattle but at a lower level than in resistant cattle. Enhanced SEMA5A induces expression of at least nine genes related to the host's immune response, including TNF-alpha and IL-8. We propose that susceptibility to mastitis results from an impaired immune response due to the lower transcription activity of susceptible FEZL. Our results provide an avenue to select for genetic improvement of resistance to mastitis and suggest that the FEZL-SEMA5A pathway might control both neuronal development and innate immunity.

Identification of a missense mutation in the bovine ABCG2 gene with a major effect on the QTL on chromosome 6 affecting milk yield and composition in Holstein cattle

We previously localized a quantitative trait locus (QTL) on chromosome 6 affecting milk fat and protein concentration to a 4-cM confidence interval, centered on the microsatellite BM143. We characterized the genes and sequence variation in this region and identified common haplotypes spanning five polymorphic sites in the genes IBSP, SPP1, PKD2, and ABCG2 for two sires heterozygous for this QTL. Expression of SPP1 and ABCG2 in the bovine mammary gland increased from parturition through lactation. SPP1 and all the coding exons of ABCG2 and PKD2 were sequenced for these two sires. The single nucleotide change capable of encoding a substitution of tyrosine-581 to serine (Y581S) in the ABCG2 transporter was the only polymorphism corresponding to the segregation status of all 3 heterozygous and 15 homozygous sires for the QTL in the Israeli and U.S. Holstein populations. The allele substitution fixed effects on the genetic evaluations of 335 Israeli sires were -341 kg milk, +0.16% fat, and +0.13% protein (F-value = 200). No other polymorphism gave significant effect for fat and protein concentration in models that also included Y581S. The allele substitution effects on the genetic evaluations of 670 cows, daughters of two heterozygous sires, were -226 kg milk, 0.09% fat, and 0.08% protein (F-value = 394), with partial dominance towards the 581S homozygotes. We therefore propose that Y581S in ABCG2 is the causative site for this QTL.

Molecular characterization and chromosomal assignment of the bovine glycinamide ribonucleotide formyltransferase (GART) gene on cattle chromosome 1q12.1-q12.2

The mammalian glycinamide ribonucleotide formyltransferase (GART) genes encode a trifunctional polypeptide involved in the de novo purine biosynthesis. We isolated a bacterial artificial chromosome (BAC) clone containing the bovine GART gene and determined the complete DNA sequence of the BAC clone. Cloning and characterization of the bovine GART gene revealed that the bovine gene consists of 23 exons spanning approximately 27 kb. RT-PCR amplification of bovine GART in different organs showed the expression of two GART transcripts in cattle similar to human and mouse. The GART transcripts encode two proteins of 1010 and 433 amino acids, respectively. Eleven single nucleotide polymorphisms (SNPs) were detected in a mutation scan of 24 unrelated animals of three different cattle breeds, including one SNP that affects the amino acid sequence of GART. The chromosomal localization of the gene was determined by fluorescence in situ hybridization. Comparative genome analysis between cattle, human and mouse indicates that the chromosomal location of the bovine GART gene is in agreement with a previously published mapping report.

A refined radiation hybrid map of the telomeric region of bovine chromosome 18q25-q26 compared with human chromosome 19q13

Genome-wide scans have mapped economically important quantitative trait loci (QTL) for mastitis susceptibility in dairy cattle at the telomeric end of bovine chromosome 18 (BTA18). In order to increase the density of markers in this chromosomal region and to improve breakpoint resolution in the human-bovine comparative map, this study describes the chromosomal assignment of seven newly developed gene-associated markers and five microsatellites and eight previously mapped sequence tagged site markers near these QTL. The orientation of KCNJ14, BAX, CD37, NKG7, LIM2, PRKCG, TNNT1, MGC2705, RPL28, EPN1, ZNF582, ZIM2, STK13, ZNF132 and SLC27A5 on the 3000-rad radiation hybrid (RH) map of BTA18 is homologous to the organization found on the corresponding 10 Mbp of human chromosome 19q (HSA19q). The resulting bovine RH map with a length of 20.9 cR spans over about 11 cM on the bovine linkage map. The location of KCNJ14 and SLC27A5 flanking the RH map on BTA18q25-26 has been confirmed by fluorescence in situ hybridization. The data of this refined human-bovine comparative map should improve selection of candidate genes for mastitis susceptibility in dairy cattle.

An expanded comparative map of bovine chromosome 27 targeting dairy form QTL regions

At present, the density of genes on the bovine maps is extremely limited and current resolution of the human-bovine comparative map is insufficient for selection of candidate genes controlling many economic traits of interest in dairy cattle. This study describes the chromosomal mapping of 10 selected gene-associated markers to bovine linkage and radiation hybrid maps to improve the breakpoint resolution in the human-bovine comparative map near two previously identified quantitative trait loci for the linear type trait, dairy form. Two regions of conserved synteny not previously described are reported between the telomeric region of bovine chromosome 27 (BTA27) and human chromosome 3 (HSA3) p24 region and between the HSA4q34.1 region and BTA8. These data increase the number of genes positioned on the bovine gene maps, refine the human-bovine comparative map, and should improve the efficiency of candidate gene selection for the dairy form trait in cattle.

Assessment of the gene content of the chromosomal regions flanking bovine DGAT1

As a first step towards verifying the candidate status of DGAT1 as the causal gene for milk fat percentage in cattle, we constructed a bovine BAC contig spanning 576 kb of the chromosomal region containing DGAT1. High content of NotI sites (21 within the contig) indicated that the region is gene-rich. Twenty-three genes neighboring DGAT1 were mapped, including two bovine cDNA sequences that have no orthologous sequences within the NCBI sequence databases. On average, 2015 bp for each of the 23 neighboring genes were sequenced and entered into EMBL. Likewise, 10 new STS markers were established by BAC-end sequencing. Within the genes and STS markers, 55 polymorphisms were discovered. These will form the basis of future linkage disequilibrium studies to test whether any genes neighboring DGAT1 are associated with variation in milk fat percentage, thereby testing the candidate status of DGAT1.

Evaluation of single-nucleotide polymorphisms in CAPN1 for association with meat tenderness in cattle

Micromolar calcium activated neutral protease (CAPN1) was evaluated as a candidate gene for a quantitative trait locus (QTL) on BTA29 affecting meat tenderness by characterization of nucleotide sequence variation in the gene. Single-nucleotide polymorphisms (SNP) were identified by sequencing all 22 exons and 19 of the 21 introns in two sires (Piedmontese x Angus located at the U.S. Meat Animal Research Center in Clay Center, NE; Jersey x Limousin located at AgResearch in New Zealand) of independent resource populations previously shown to be segregating meat tenderness QTL on BTA29. The majority of the 38 SNP were found in introns or were synonymous substitutions in the coding regions, with two exceptions. Exons 14 and 9 contained SNP that were predicted to alter the protein sequence by the substitution of isoleucine for valine in Domain III of the protein, and alanine for glycine in Domain II of the protein. The resource populations were genotyped for these two SNP in addition to six intronic polymorphisms and two silent substitutions. Analysis of genotypes and shear force values in both populations revealed a difference between paternal CAPN1 alleles in which the allele encoding isoleucine at position 530 and glycine at position 316 associated with decreased meat tenderness (increased shear force values) relative to the allele encoding valine at position 530 and alanine at position 316 (P < 0.05). The association of maternal alleles with meat tenderness phenotypes is consistent with the hypothesis of CAPN1 as the gene underlying the QTL effect in two independent resource populations and presents the possibility of using these markers for selective breeding to reduce the numbers of animals with unfavorable meat tenderness traits.

Identification of the single base change causing the callipyge muscle hypertrophy phenotype, the only known example of polar overdominance in mammals

A small genetic region near the telomere of ovine chromosome 18 was previously shown to carry the mutation causing the callipyge muscle hypertrophy phenotype in sheep. Expression of this phenotype is the only known case in mammals of paternal polar overdominance gene action. A region surrounding two positional candidate genes was sequenced in animals of known genotype. Mutation detection focused on an inbred ram of callipyge phenotype postulated to have inherited chromosome segments identical-by-descent with exception of the mutated position. In support of this hypothesis, this inbred ram was homozygous over 210 Kb of sequence, except for a single heterozygous base position. This single polymorphism was genotyped in multiple families segregating the callipyge locus (CLPG), providing 100% concordance with animals of known CLPG genotype, and was unique to descendants of the founder animal. The mutation lies in a region of high homology among mouse, sheep, cattle, and humans, but not in any previously identified expressed transcript. A substantial open reading frame exists in the sheep sequence surrounding the mutation, although this frame is not conserved among species. Initial functional analysis indicates sequence encompassing the mutation is part of a novel transcript expressed in sheep fetal muscle we have named CLPG1.

Parallel construction of orthologous sequence-ready clone contig maps in multiple species

Comparison is a fundamental tool for analyzing DNA sequence. Interspecies sequence comparison is particularly powerful for inferring genome function and is based on the simple premise that conserved sequences are likely to be important. Thus, the comparison of a genomic sequence with its orthologous counterpart from another species is increasingly becoming an integral component of genome analysis. In ideal situations, such comparisons are performed with orthologous sequences from multiple species. To facilitate multispecies comparative sequence analysis, a robust and scalable strategy for simultaneously constructing sequence-ready bacterial artificial chromosome (BAC) contig maps from targeted genomic regions has been developed. Central to this approach is the generation and utilization of "universal" oligonucleotide-based hybridization probes ("overgo" probes), which are designed from sequences that are highly conserved between distantly related species. Large collections of these probes are used en masse to screen BAC libraries from multiple species in parallel, with the isolated clones assembled into physical contig maps. To validate the effectiveness of this strategy, efforts were focused on the construction of BAC-based physical maps from multiple mammalian species (chimpanzee, baboon, cat, dog, cow, and pig). Using available human and mouse genomic sequence and a newly developed computer program to design the requisite probes, sequence-ready maps were constructed in all species for a series of targeted regions totaling approximately 16 Mb in the human genome. The described approach can be used to facilitate the multispecies comparative sequencing of targeted genomic regions and can be adapted for constructing BAC contig maps in other vertebrates.

A high-resolution comparative RH map of the proximal part of bovine chromosome 1

Current comparative maps between human chromosome 21 and the proximal part of cattle chromosome 1 are insufficient to define chromosomal rearrangements because of the low density of mapped genes in the bovine genome. The recently completed sequence of human chromosome 21 facilitates the detailed comparative analysis of corresponding segments on BTA1. In this study eight bovine bacterial artificial chromosome (BAC) clones containing bovine orthologues of human chromosome 21 genes, i.e. GRIK1, CLDN8, TIAM1, HUNK, SYNJ1, OLIG2, IL10RB, and KCNE2 were physically assigned by fluorescence in situ hybridization (FISH) to BTA1q12.1-q12.2. Sequence tagged site (STS) markers derived from these clones were mapped on the 3000 rad Roslin/Cambridge bovine radiation hybrid (RH) panel. In addition to these eight novel markers, 17 known markers from previously published BTA1 linkage or RH maps were also mapped on the Roslin/Cambridge bovine RH panel resulting in an integrated map with 25 markers of 355.4 cR(3000) length. The human-cattle genome comparison revealed the existence of three chromosomal breakpoints and two probable inversions in this region.

Association of a lysine-232/alanine polymorphism in a bovine gene encoding acyl-CoA:diacylglycerol acyltransferase (DGAT1) with variation at a quantitative trait locus for milk fat content

DGAT1 encodes diacylglycerol O-acyltransferase (EC ), a microsomal enzyme that catalyzes the final step of triglyceride synthesis. It became a functional candidate gene for lactation traits after studies indicated that mice lacking both copies of DGAT1 are completely devoid of milk secretion, most likely because of deficient triglyceride synthesis in the mammary gland. Our mapping studies placed DGAT1 close to the region of a quantitative trait locus (QTL) on bovine chromosome 14 for variation in fat content of milk. Sequencing of DGAT1 from pooled DNA revealed significant frequency shifts at several variable positions between groups of animals with high and low breeding values for milk fat content in different breeds (Holstein-Friesian, Fleckvieh, and Braunvieh). Among the variants was a nonconservative substitution of lysine by alanine (K232A), with the lysine-encoding allele being associated with higher milk fat content. Haplotype analysis indicated the lysine variant to be ancestral. Two animals that were typed heterozygous (Qq) at the QTL based on marker-assisted QTL-genotyping were heterozygous for the K232A substitution, whereas 14 animals that are most likely qq at the QTL were homozygous for the alanine-encoding allele. An independent association study in Fleckvieh animals confirmed the positive effect of the lysine variant on milk fat content. We consider the nonconservative K232A substitution to be directly responsible for the QTL variation, although our genetic studies cannot provide formal proof.

Positional candidate cloning of a QTL in dairy cattle: identification of a missense mutation in the bovine DGAT1 gene with major effect on milk yield and composition

We recently mapped a quantitative trait locus (QTL) with a major effect on milk composition--particularly fat content--to the centromeric end of bovine chromosome 14. We subsequently exploited linkage disequilibrium to refine the map position of this QTL to a 3-cM chromosome interval bounded by microsatellite markers BULGE13 and BULGE09. We herein report the positional candidate cloning of this QTL, involving (1) the construction of a BAC contig spanning the corresponding marker interval, (2) the demonstration that a very strong candidate gene, acylCoA:diacylglycerol acyltransferase (DGAT1), maps to that contig, and (3) the identification of a nonconservative K232A substitution in the DGAT1 gene with a major effect on milk fat content and other milk characteristics.

Development of microsatellite markers and comparative mapping for bovine chromosome 19

Previous research has mapped an ovulation rate quantitative trait locus (QTL) to bovine chromosome 19. In an effort to enhance comparative mapping information and develop additional markers for refined QTL mapping, microsatellite markers were developed in a targeted approach. A bovine bacterial artificial chromosome (BAC) library was screened for loci with either known or predicted locations on bovine chromosome 19. An average of 6.4 positive BAC were identified per screened locus. A total of 10 microsatellite markers were developed for five targeted loci with heterozygosity of 7-83% in a sample of reference family parents. The newly developed markers were typed on reference families along with four previously mapped marker loci and used to create a linkage map. Comparison of locus order between human and cattle provides support for previously observed rearrangement. One of the mapped loci myotubularin related protein 4 (MTMR4) potentially extends the proximal boundary of a conserved linkage group.

Nuclear transfer technology in mammalian cloning

The past several years have witnessed remarkable progress in mammalian cloning using nuclear transfer (NT). Until 1997 and the announcement of the successful cloning of sheep from adult mammary gland or fetal fibroblast cells, our working assumption was that cloning by NT could only be accomplished with relatively undifferentiated embryonic cells. Indeed, live offspring were first produced by NT over 15 years ago from totipotent, embryonic blastomeres derived from early cleavage-stage embryos. However, once begun, the progression to somatic cell cloning or NT employing differentiated cells as the source of donor nuclei was meteoric, initially involving differentiated embryonic cell cultures in sheep in 1996 and quickly thereafter, fetal or adult somatic cells in sheep, cow, mouse, goat, and pig. Several recent reviews provide a background for and discussion of these successes. Here we will focus on the potential uses of reproductive cloning along with recent activities in the field and a discussion concerning current interests in human reproductive and therapeutic cloning.

Partial deletion of the bovine ED1 gene causes anhidrotic ectodermal dysplasia in cattle

Anhidrotic ectodermal dysplasia (ED1) is characterized by hypotrichosis, reduced number of sweat glands, and incisior anodontia in human, mouse, and cattle. In affected humans and mice, mutations in the ED1 gene coding for ectodysplasin 1 are found. Ectodysplasin 1 is a novel trimeric transmembrane protein with an extracellular TNF-like signaling domain that is believed to be involved in the formation of hair follicles and tooth buds during fetal development. We report the construction of a 480-kb BAC contig harboring the complete bovine ED1 gene on BTA Xq22-Xq24. Physical mapping and sequence analysis of the coding parts of the ED1 gene revealed that a large genomic region including exon 3 of the ED1 gene is deleted in cattle with anhidrotic ectodermal dysplasia in a family of German Holstein cattle with three affected maternal half sibs.

Mapping of expressed sequence tags from a porcine early embryonic cDNA library

The goal of this study was to identify and map genes expressed during the elongation phase of embryogenesis in swine. Expressed sequence tags were analysed from a previously described porcine cDNA library prepared from elongating swine embryos. Average insert length of randomly selected clones was approximately 600 bp, with a range from < 100 to > 2500 bp. Single-pass, coding strand sequences from 1132 independent clones were compared with the GenBank non-redundant (nr) database via BLASTN analysis to identify potential porcine homologous of known genes. Among these sequences, 781 (69%) showed significant (score > 300) homology to non- mitochondrial sequences previously deposited in GenBank. Sequences matching interleucin 1 beta and thymosin beta 10 were most frequently observed (24 and 18 clones, respectively), in addition to matches with 310 other distinct genes. No significant match in the GenBank nr database was obtained for 303 sequences. Analysis demonstrated that 151 (50%) had open reading frames (ORF) extending at least 50 codons from the first base of the clone insert. Genetic markers were developed and used to map a subset of 17 genes, selected on the basis of function or of the ability to design primers that successfully amplified porcine genomic DNA, to 10 different porcine chromosomes, providing a set of mapped markers corresponding to genes expressed during conceptus elongation.