The sheep genome reference sequence: a work in progress.

Until recently, the construction of a reference genome was performed using Sanger sequencing alone. The emergence of next-generation sequencing platforms now means reference genomes may incorporate sequence data generated from a range of sequencing platforms, each of which have different read length, systematic biases and mate-pair characteristics. The objective of this review is to inform the mammalian genomics community about the experimental strategy being pursued by the International Sheep Genomics Consortium (ISGC) to construct the draft reference genome of sheep (Ovis aries). Component activities such as data generation, sequence assembly and annotation are described, along with information concerning the key researchers performing the work. This aims to foster future participation from across the research community through the coordinated activities of the consortium. The review also serves as a 'marker paper' by providing information concerning the pre-publication release of the reference genome. This ensures the ISGC adheres to the framework for data sharing established at the recent Toronto International Data Release Workshop and provides guidelines for data users.

Molecular cytogenetics and gene mapping in sheep (Ovis aries, 2n = 54).

The development of a completely annotated sheep genome sequence is a key need for understanding the phylogenetic relationships and genetic diversity among the many different sheep breeds worldwide and for identifying genes controlling economically and physiologically important traits. The ovine genome sequence assembly will be crucial for developing optimized breeding programs based on highly productive, healthy sheep phenotypes that are adapted to modern breeding and production conditions. Scientists and breeders around the globe have been contributing to this goal by generating genomic and cDNA libraries, performing genome-wide and trait-associated analyses of polymorphism, expression analysis, genome sequencing, and by developing virtual and physical comparative maps. The International Sheep Genomics Consortium (ISGC), an informal network of sheep genomics researchers, is playing a major role in coordinating many of these activities. In addition to serving as an essential tool for monitoring chromosome abnormalities in specific sheep populations, ovine molecular cytogenetics provides physical anchors which link and order genome regions, such as sequence contigs, genes and polymorphic DNA markers to ovine chromosomes. Likewise, molecular cytogenetics can contribute to the process of defining evolutionary breakpoints between related species. The selective expansion of the sheep cytogenetic map, using loci to connect maps and identify chromosome bands, can substantially contribute to improving the quality of the annotated sheep genome sequence and will also accelerate its assembly. Furthermore, identifying major morphological chromosome anomalies and micro-rearrangements, such as gene duplications or deletions, that might occur between different sheep breeds and other Ovis species will also be important to understand the diversity of sheep chromosome structure and its implications for cross-breeding. To date, 566 loci have been assigned to specific chromosome regions in sheep and the new cytogenetic map is presented as part of this review. This review will also summarize the current cytogenomic status of the sheep genome, describe current activities in the sheep cytogenomics research sector, and will discuss the cytogenomics data in context with other major sheep genomics projects.

A radiation hybrid comparative map of ovine chromosome 1 aligned to the virtual sheep genome.

Ovis aries chromosome one (OAR1) is the largest submetacentric chromosome in the sheep genome and is homologous to regions on human chromosomes 1, 2, 3 and 21. Using the USUoRH5000 ovine whole-genome radiation hybrid (RH) panel, we have constructed a RH map of OAR1 comprising 102 framework and 75 placed/binned markers across five linkage groups spanning 3759.43 cR5000, with an average marker density of 21.2 cR5000/marker. The alignment of our OAR1 RH map shows good concordance with the recently developed virtual sheep genome, with fewer than 1.86% discrepancies. A comparative map of OAR1 was constructed by examining the location of RH-mapped orthologues in sheep within the genomes of cow, human, horse and dog. Analysis of the comparative map indicates that conserved syntenies within the five ovine RH linkage groups underwent internal chromosomal rearrangements which, in general, reflect the evolutionary distances between sheep and each of these four species. The ovine RH map presented here integrates all available mapping data and includes new genomic information for ovine chromosome 1.

Genetic mapping of quantitative trait loci for resistance to Haemonchus contortus in sheep.

This paper presents results from a mapping experiment to detect quantitative trait loci (QTL) for resistance to Haemonchus contortus infestation in merino sheep. The primary trait analysed was faecal worm egg count in response to artificial challenge at 6 months of age. In the first stage of the experiment, whole genome linkage analysis was used for broad-scale mapping. The animal resource used was a designed flock comprising 571 individuals from four half-sib families. The average marker spacing was about 20 cM. For the primary trait, 11 QTL (as chromosomal/family combinations) were significant at the 5% chromosome-wide level, with allelic substitution effects of between 0.19 and 0.38 phenotypic standard deviation units. In general, these QTL did not have a significant effect on faecal worm egg count recorded at 13 months of age. In the second stage of the experiment, three promising regions (located on chromosomes 1, 3 and 4) were fine-mapped. This involved typing more closely spaced markers on individuals from the designed flock as well as an additional 495 individuals selected from a related population with a deeper pedigree. Analysis was performed using a linkage disequilibrium-linkage approach, under additive, dominant and multiple QTL models. Of these, the multiple QTL model resulted in the most refined QTL positions, with resolutions of <10 cM achieved for two regions. Because of the moderate size of effect of the QTL, and the apparent age and/or immune status specificity of the QTL, it is suggested that a panel of QTL will be required for significant genetic gains to be achieved within industry via marker-assisted selection.

A high-resolution comparative radiation hybrid map of ovine chromosomal regions that are homologous to human chromosome 6 (HSA6).

In this study, we constructed high-resolution radiation hybrid (RH) and comparative maps of ovine chromosomes or chromosomal segments that are homologous to human chromosome 6 (HSA6). A total of 251 markers were successfully genotyped across the recently developed USUoRH5000 whole-genome panel; 208 of these markers were assigned to five RH linkage groups distributed on three ovine chromosomes (OAR8, 9 and 20). The RH maps have good correspondence with previous chromosome painting data, although a small centromeric region on OAR9 that is homologous to HSA6 had not been previously detected using human chromosome paints on ovine chromosomal spreads. High percentages of the ovine markers were identified as orthologues in the bovine (86.3%), dog (85.8%), horse (69.3%) and human (88.7%) genomes. These maps contribute to investigations in mammalian chromosome evolution and the search for economic trait loci in sheep.

A radiation hybrid map of sheep chromosome 23 based on ovine BAC-end sequences.

More than 375,000 BAC-end sequences (BES) of the CHORI-243 ovine BAC library have been deposited in public databases. blastn searches with these BES against HSA18 revealed 1806 unique and significant hits. We used blastn-anchored BES for an in silico prediction of gene content and chromosome assignment of comparatively mapped ovine BAC clones. Ovine BES were selected at approximately 1.3-Mb intervals of HSA18 and incorporated into a human-sheep comparative map. An ovine 5000-rad whole-genome radiation hybrid panel (USUoRH5000) was typed with 70 markers, all of which mapped to OAR23. The resulting OAR23 RH map included 43 markers derived from BES with high and unique BLAST hits to the sequence of the orthologous HSA18, nine EST-derived markers, 16 microsatellite markers taken from the ovine linkage map and two bovine microsatellite markers. Six new microsatellite markers derived from the 43 mapped BES and the two bovine microsatellite markers were linkage-mapped using the International Mapping Flock (IMF). Thirteen additional microsatellite markers were derived from other ovine BES with high and unique BLAST hits to the sequence of the orthologous HSA18 and also positioned on the ovine linkage map but not incorporated into the OAR23 RH map. This resulted in 24 markers in common and in the same order between the RH and linkage maps. Eight of the BES-derived markers were mapped using fluorescent in situ hybridization (FISH), to thereby align the RH and cytogenetic maps. Comparison of the ovine chromosome 23 RH map with the HSA18 map identified and localized three major breakpoints between HSA18 and OAR23. The positions of these breakpoints were equivalent to those previously shown for syntenic BTA24 and HSA18. This study presents evidence for the usefulness of ovine BES when constructing a high-resolution comprehensive map for a single sheep chromosome. The comparative analysis confirms and refines knowledge about chromosomal conservation and rearrangements between sheep, cattle and human. The constructed RH map demonstrates the resolution and utility of the newly constructed ovine RH panel.

Globally dispersed Y chromosomal haplotypes in wild and domestic sheep.

To date, investigations of genetic diversity and the origins of domestication in sheep have utilised autosomal microsatellites and variation in the mitochondrial genome. We present the first analysis of both domestic and wild sheep using genetic markers residing on the ovine Y chromosome. Analysis of a single nucleotide polymorphism (oY1) in the SRY promoter region revealed that allele A-oY1 was present in all wild bighorn sheep (Ovis canadensis), two subspecies of thinhorn sheep (Ovis dalli), European Mouflon (Ovis musimon) and the Barbary (Ammontragis lervia). A-oY1 also had the highest frequency (71.4%) within 458 domestic sheep drawn from 65 breeds sampled from Africa, Asia, Australia, the Caribbean, Europe, the Middle East and Central Asia. Sequence analysis of a second locus, microsatellite SRYM18, revealed a compound repeat array displaying fixed differences, which identified bighorn and thinhorn sheep as distinct from the European Mouflon and domestic animals. Combined genotypic data identified 11 male-specific haplotypes that represented at least two separate lineages. Investigation of the geographical distribution of each haplotype revealed that one (H6) was both very common and widespread in the global sample of domestic breeds. The remaining haplotypes each displayed more restricted and informative distributions. For example, H5 was likely founded following the domestication of European breeds and was used to trace the recent transportation of animals to both the Caribbean and Australia. A high rate of Y chromosomal dispersal appears to have taken place during the development of domestic sheep as only 12.9% of the total observed variation was partitioned between major geographical regions.

Systematic screening of sheep skin cDNA libraries for microsatellite sequences.

65,000 sheep skin cDNA clones were gridded in high density on to nylon membranes and screened for (CA)n and (GA)n repeat containing clones. 296 dinucleotide repeat-containing clones were identified with approximately 85% non-redundancy. Clones were single-pass 5' sequenced and we compared the Expressed Sequence Tag (EST) sequences to the Swiss-Prot database to ascertain their identity and/or putative function. We then aligned the ESTs against the human genomic sequence to determine the locations of human orthologous sequences. Finally, we developed a subset of polymorphic microsatellite markers and positioned them on the ovine linkage map.

A genome scan for quantitative trait loci affecting resistance to Trichostrongylus colubriformis in sheep.

A genome linkage scan was carried out using a resource flock of 1029 sheep in six half-sib families. The families were offspring of sires derived by crossing divergent lines of sheep selected for response to challenge with the intestinal parasitic nematode Trichostrongylus colubriformis. All animals in the resource flock were phenotypically assessed for worm resistance soon after weaning using a vaccination/challenge regime. After correcting for fixed effects using a least squares linear model the faecal egg count data obtained following the first challenge and the faecal egg count data obtained after the second challenge were designated Trait 1 and Trait 2, respectively. A total of 472 lambs drawn from the phenotypic extremes of the Trait 2 faecal egg count distribution were genotyped with a panel of 133 microsatellite markers covering all 26 sheep autosomes. Detection of quantitative trait loci (QTL) for each of the faecal egg count traits was determined using interval analysis with the Animap program with recombination rates between markers derived from an existing marker map. No chromosomal regions attained genome-wide significance for QTL influencing either of the traits. However, one region attained chromosome-wide significance and five other regions attained point-wise significance for the presence of QTL affecting parasite resistance.

An enhanced linkage map of the sheep genome comprising more than 1000 loci.

A medium-density linkage map of the ovine genome has been developed. Marker data for 550 new loci were generated and merged with the previous sheep linkage map. The new map comprises 1093 markers representing 1062 unique loci (941 anonymous loci, 121 genes) and spans 3500 cM (sex-averaged) for the autosomes and 132 cM (female) on the X chromosome. There is an average spacing of 3.4 cM between autosomal loci and 8.3 cM between highly polymorphic [polymorphic information content (PIC) > or = 0.7] autosomal loci. The largest gap between markers is 32.5 cM, and the number of gaps of > 20 cM between loci, or regions where loci are missing from chromosome ends, has been reduced from 40 in the previous map to 6. Five hundred and seventy-three of the loci can be ordered on a framework map with odds of > 1000 : 1. The sheep linkage map contains strong links to both the cattle and goat maps. Five hundred and seventy-two of the loci positioned on the sheep linkage map have also been mapped by linkage analysis in cattle, and 209 of the loci mapped on the sheep linkage map have also been placed on the goat linkage map. Inspection of ruminant linkage maps indicates that the genomic coverage by the current sheep linkage map is comparable to that of the available cattle maps. The sheep map provides a valuable resource to the international sheep, cattle, and goat gene mapping community.

Increased 15-HPETE production decreases prostacyclin synthase activity during oxidant stress in aortic endothelial cells.

Selenium (Se) is an integral component of glutathione peroxidase and is able to detoxify peroxides that can affect arachidonic acid (AA) metabolism, thereby influencing eicosanoid biosynthesis. This study investigated the effects of oxidant stress, a consequence of Se deficiency, on eicosanoid formation and important key enzyme expression in bovine aortic endothelial cells (BAEC). Bovine aortic endothelial cells cultured in Se-deficient media and stimulated with tumor necrosis factor alpha or H2O2 produced significantly less prostacyclin (PGI(2)) and more 15-hydroxyeicosatetraenoic acid, 15-hydroperoxyeicosatetraenoic acid (15-HPETE), and thromboxane than Se-supplemented BAEC. Additionally, reverse transcription polymerase chain reaction and immunoblotting determined that the mRNA and protein levels of the eicosanoid forming enzymes cyclooxygenase-1 (COX1), cyclooxygenase-2 (COX2), and PGI synthase were not significantly changed. The addition of 15-HPETE to Se-supplemented BAEC inhibited the production of PGI(2) suggesting that the accumulation of lipid hydroperoxides during Se-deficiency may be the underlying factor in the altered eicosanoid production during Se deficiency. Furthermore, inhibition of COX and addition of PGH(2) to Se-deficient or Se-supplemented BAEC still resulted in lower PGI(2) formation by Se-deficient cells. Together, these results suggest that Se deficiency modifies eicosanoid production by affecting the activity of key enzymes, particularly PGI synthase, rather than their transcription or translation.