A. M. Crawford

A second-generation linkage map of the sheep genome.

A genetic map of Ovis aries (haploid n = 27) was developed with 519 markers (504 microsatellites) spanning ∼3063 cM in 26 autosomal linkage groups and 127 cM (female specific) of the X Chromosome (Chr). Genotypic data were merged from the IMF flock (Crawford et al., Genetics 140, 703, 1995) and the USDA mapping flock. Seventy-three percent (370/504) of the microsatellite markers on the map are common to the USDA-ARS MARC cattle linkage map, with 27 of the common markers derived from sheep. The number of common markers per homologous linkage group ranges from 5 to 22 and spans a total of 2866 cM (sex average) in sheep and 2817 cM in cattle. Marker order within a linkage group was consistent between the two species with limited exceptions. The reported translocation between the telomeric end of bovine Chr 9 (BTA 9) and BTA 14 to form ovine Chr 9 is represented by a 15-cM region containing 5 common markers. The significant genomic conservation of marker order will allow use of linkage maps in both species to facilitate the search for quantitative trait loci (QTLs) in cattle and sheep.

Mapping of the melatonin receptor 1a (MTNR1A) gene in pigs, sheep, and cattle.

and Implications Human and sheep Melatonin receptor 1a (MTNR1A) gene information was used to clone a portion of the coding region of this gene in pigs, and to identify polymorphisms of the gene for its assignment to both the genetic linkage and physical maps. MTNR1A maps to pig chromosome 17, establishing a new region of conserved synteny between this chromosome and human chromosome 4. Furthermore, we have assigned MTNR1A to bovine chromosome 27 and sheep chromosome 26. The addition of genes like MTNR1A to livestock genome maps allows questions about evolutionary events and the genetic basis for quantitative traits in livestock to be addressed.

Using comparative genomics to reorder the human genome sequence into a virtual sheep genome

BackgroundIs it possible to construct an accurate and detailed subgene-level map of a genome using bacterial artificial chromosome (BAC) end sequences, a sparse marker map, and the sequences of other genomes?ResultsA sheep BAC library, CHORI-243, was constructed and the BAC end sequences were determined and mapped with high sensitivity and low specificity onto the frameworks of the human, dog, and cow genomes. To maximize genome coverage, the coordinates of all BAC end sequence hits to the cow and dog genomes were also converted to the equivalent human genome coordinates. The 84,624 sheep BACs (about 5.4-fold genome coverage) with paired ends in the correct orientation (tail-to-tail) and spacing, combined with information from sheep BAC comparative genome contigs (CGCs) built separately on the dog and cow genomes, were used to construct 1,172 sheep BAC-CGCs, covering 91.2% of the human genome. Clustered non-tail-to-tail and outsize BACs located close to the ends of many BAC-CGCs linked BAC-CGCs covering about 70% of the genome to at least one other BAC-CGC on the same chromosome. Using the BAC-CGCs, the intrachromosomal and interchromosomal BAC-CGC linkage information, human/cow and vertebrate synteny, and the sheep marker map, a virtual sheep genome was constructed. To identify BACs potentially located in gaps between BAC-CGCs, an additional set of 55,668 sheep BACs were positioned on the sheep genome with lower confidence. A coordinate conversion process allowed us to transfer human genes and other genome features to the virtual sheep genome to display on a sheep genome browser.ConclusionWe demonstrate that limited sequencing of BACs combined with positioning on a well assembled genome and integrating locations from other less well assembled genomes can yield extensive, detailed subgene-level maps of mammalian genomes, for which genomic resources are currently limited.

Gene expression profiling of Naïve sheep genetically resistant and susceptible to gastrointestinal nematodes

BackgroundGastrointestinal nematodes constitute a major cause of morbidity and mortality in grazing ruminants. Individual animals or breeds, however, are known to differ in their resistance to infection. Gene expression profiling allows us to examine large numbers of transcripts simultaneously in order to identify those transcripts that contribute to an animals susceptibility or resistance.ResultsWith the goal of identifying genes with a differential pattern of expression between sheep genetically resistant and susceptible to gastrointestinal nematodes, a 20,000 spot ovine cDNA microarray was constructed. This array was used to interrogate the expression of 9,238 known genes in duodenum tissue of four resistant and four susceptible female lambs. Naïve animals were used in order to look at genes that were differentially expressed in the absence of infection with gastrointestinal nematodes. Forty one unique known genes were identified that were differentially expressed between the resistant and susceptible animals. Northern blotting of a selection of the genes confirmed differential expression. The differentially expressed genes had a variety of functions, although many genes relating to the stress response and response to stimulus were more highly expressed in the susceptible animals.ConclusionWe have constructed the first reported ovine microarray and used this array to examine gene expression in lambs genetically resistant and susceptible to gastrointestinal nematode infection. This study indicates that susceptible animals appear to be generating a hyper-sensitive immune response to non-nematode challenges. The gastrointestinal tract of susceptible animals is therefore under stress and compromised even in the absence of gastrointestinal nematodes. These factors may contribute to the genetic susceptibility of these animals.

Discovery of quantitative trait loci for resistance to parasitic nematode infection in sheep: I. Analysis of outcross pedigrees

BackgroundCurrently most pastoral farmers rely on anthelmintic drenches to control gastrointestinal parasitic nematodes in sheep. Resistance to anthelmintics is rapidly increasing in nematode populations such that on some farms none of the drench families are now completely effective. It is well established that host resistance to nematode infection is a moderately heritable trait. This study was undertaken to identify regions of the genome, quantitative trait loci (QTL) that contain genes affecting resistance to parasitic nematodes.ResultsRams obtained from crossing nematode parasite resistant and susceptible selection lines were used to derive five large half-sib families comprising between 348 and 101 offspring per sire. Total offspring comprised 940 lambs. Extensive measurements for a range of parasite burden and immune function traits in all offspring allowed each lamb in each pedigree to be ranked for relative resistance to nematode parasites.Initially the 22 most resistant and 22 most susceptible progeny from each pedigree were used in a genome scan that used 203 microsatellite markers spread across all sheep autosomes. This study identified 9 chromosomes with regions showing sufficient linkage to warrant the genotyping of all offspring. After genotyping all offspring with markers covering Chromosomes 1, 3, 4, 5, 8, 12, 13, 22 and 23, the telomeric end of chromosome 8 was identified as having a significant QTL for parasite resistance as measured by the number of Trichostrongylus spp. adults in the abomasum and small intestine at the end of the second parasite challenge. Two further QTL for associated immune function traits of total serum IgE and T. colubiformis specific serum IgG, at the end of the second parasite challenge, were identified on chromosome 23.ConclusionDespite parasite resistance being a moderately heritable trait, this large study was able to identify only a single significant QTL associated with it. The QTL concerned adult parasite burdens at the end of the second parasite challenge when the lambs were approximately 6 months old. Our failure to discover more QTL suggests that most of the genes controlling this trait are of relatively small effect. The large number of suggestive QTL discovered (more than one per family per trait than would be expected by chance) also supports this conclusion.

A physical map of the bovine genome

BackgroundCattle 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.ResultsA 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.ConclusionFurther 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 paradox of MHC-DRB exon/intron evolution: α-helix and β-sheet encoding regions diverge while hypervariable intronic simple repeats coevolve with β-sheet codons

Twenty-one different caprine and 13 ovine MHC-DRB exon 2 sequences were determined including part of the adjacent introns containing simple repetitive (gt)n(ga)m elements. The positions for highly polymorphic DRB amino acids vary slightly among ungulates and other mammals. From man and mouse to ungulates the basic (gt)n(ga)m structure is fixed in evolution for 7 × 107 years whereas ample variations exist in the tandem (gt)n and (ga)m dinucleotides and especially their “degenerated” derivatives. Phylogenetic trees for the α-helices and β-pleated sheets of the ungulate DRB sequences suggest different evolutionary histories. In hoofed animals as well as in humans DRB β-sheet encoding sequences and adjacent intronic repeats can be assembled into virtually identical groups suggesting coevolution of noncoding as well as coding DNA. In contrast a-helices and C-terminal parts of the first DRB domain evolve distinctly. In the absence of a defined mechanism causing specific, site-directed mutations, double-recombination or gene-conversion-like events would readily explain this fact. The role of the intronic simple (gt)n(ga)m repeat is discussed with respect to these genetic exchange mechanisms during evolution.

Microsatellite Evolution: Testing the Ascertainment Bias Hypothesis

Abstract. Previous studies suggest the median allele length of microsatellites is longest in the species from which the markers were derived, suggesting that an ascertainment bias was operating. We have examined whether the size distribution of microsatellite alleles between sheep and cattle is source dependent using a set of 472 microsatellites that can be amplified in both species. For those markers that were polymorphic in both species we report a significantly greater number of markers (P < 0.001) with longer median allele sizes in sheep, regardless of microsatellite origin. This finding suggests that any ascertainment bias operating during microsatellite selection is only a minor contributor to the variation observed.

Microsatellites and associated repetitive elements in the sheep genome

To determine the frequency and clustering of a variety of simple di-and trinucleotide repeats, an Artiodactyl short interspersed element (SINE), an ovine satellite repeat, and a human Alu 1 repeat were used to screen a random selection of cosmids containing inserts of ovine genomic DNA. In total, 197 individual cosmids were digested with EcoRI and the fragments separated on 0.7% agarose gels. Southern blots of these gels were then sequentially probed with (AC)7, (CT)9, and (CAC)6 oligonucleotides, and the repeats described above. The frequency at which (AC)1, (CT)n, and (CAC)n repeats were found in the cosmids indicated that they occurred at average intervals of 65 kb, 367 kb, and 213 kb respectively within the ovine genome. The Artiodactyl SINE was the most common, occurring at an average interval of 20 kb. No human Alu 1 sequences were detected. There was a significant positive association between the (AC)n and the Artiodactyl SINE. This association is quite strong as there was significant clustering of the two repeats both within cosmids and also within the EcoRI fragments of the digested genomic fragments. With the exception of the sheep satellite sequence, which occurs in tandem arrays, none of the other repeats showed significant clustering within the 41-kb (average size) cosmid inserts. The first 25 ovine microsatellites we characterized had an average polymorphic information content (PIC) of 0.65. The different microsatellite types, containing either perfect, imperfect, or compound repeats, had similar average PICs of 0.64, 0.65, and 0.66 respectively. There was a weak regression relationship (R2(adj)%=21.9) between the length of the longest uninterrupted dinucleotide repeat in the largest allele and the PIC of the microsatellite.

Identification and genetic mapping of random amplified polymorphic DNA (RAPD) markers to the sheep genome

The random amplified polymorphic DNA (RAPD) assay utilizes the polymerase chain reaction (PCR) and short primers of arbitrary nucleotide sequence to amplify DNA. In this study, the RAPD assay was used to identify and map polymorphic markers in the AgResearch International Mapping Flock (IMF) sheep pedigrees. Sires and dams of eight of the full-sib IMF pedigrees were screened with 131 different 10-mer oligonucleotide primers. An average of 85 RAPD polymorphisms was identified between each parental pair, and 53 markers were contributed to the AgResearch IMF collaboration. Forty-five of the RAPD markers were mapped in the AgResearch IMF genetic linkage map, and at least one marker was located on 17 of the 26 autosomes and both sex chromosomes. Three lines of evidence were used to check for the homology of scored polymorphisms in different pedigrees, pedigree evaluation, segregation analysis, and Southern blot analysis. These results demonstrate that the RAPD assay is a powerful approach for identifying polymorphisms that can be used as markers for constructing a sheep genetic linkage map.