Assessment of genetic diversity and population structure in cultured Australian Pacific oysters.

The recent development of Pacific oyster (Crassostrea gigas) SNP genotyping arrays has allowed detailed characterisation of genetic diversity and population structure within and between oyster populations. It also raises the potential of harnessing genomic selection for genetic improvement in oyster breeding programmes. The aim of this study was to characterise a breeding population of Australian oysters through genotyping and analysis of 18 027 SNPs, followed by comparison with genotypes of oyster sampled from Europe and Asia. This revealed that the Australian populations had similar population diversity (HE ) to oysters from New Zealand, the British Isles, France and Japan. Population divergence was assessed using PCA of genetic distance and revealed that Australian oysters were distinct from all other populations tested. Australian Pacific oysters originate from planned introductions sourced from three Japanese populations. Approximately 95% of these introductions were from geographically, and potentially genetically, distinct populations from the Nagasaki oysters assessed in this study. Finally, in preparation for the application of genomic selection in oyster breeding programmes, the strength of LD was evaluated and subsets of loci were tested for their ability to accurately infer relationships. Weak LD was observed on average; however, SNP subsets were shown to accurately reconstitute a genomic relationship matrix constructed using all loci. This suggests that low-density SNP panels may have utility in the Australian population tested, and the findings represent an important first step towards the design and implementation of genomic approaches for applied breeding in Pacific oysters.

Selection signature analysis reveals genes associated with tail type in Chinese indigenous sheep.

Fat-tailed sheep have commercial value because consumers prefer high-protein and low-fat food and producers care about feed conversion rate. However, fat-tailed sheep still have some scientific significance, as the fat tail is commonly regarded as a characteristic of environmental adaptability. Finding the candidate genes associated with fat tail formation is essential for breeding and conservation. To identify these candidate genes, we applied FST and hapFLK approaches in fat- and thin-tailed sheep with available 50K SNP genotype data. These two methods found 6.24 Mb of overlapped regions and 43 genes that may associated with fat tail development. Gene annotation showed that HOXA11, BMP2, PPP1CC, SP3, SP9, WDR92, PROKR1 and ETAA1 may play important roles in fat tail formation. These findings provide insight into tail fat development and a guide for molecular breeding and conservation.

Beta-globin gene evolution in the ruminants: evidence for an ancient origin of sheep haplotype B.

Domestic sheep (Ovis aries) can be divided into two groups with significantly different responses to hypoxic environments, determined by two allelic beta-globin haplotypes. Haplotype A is very similar to the goat beta-globin locus, whereas haplotype B has a deletion spanning four globin genes, including beta-C globin, which encodes a globin with high oxygen affinity. We surveyed the beta-globin locus using resequencing data from 70 domestic sheep from 42 worldwide breeds and three Ovis canadensis and two Ovis dalli individuals. Haplotype B has an allele frequency of 71.4% in O. aries and was homozygous (BB) in all five wild sheep. This shared ancestry indicates haplotype B is at least 2-3 million years old. Approximately 40 kb of the sequence flanking the ~37-kb haplotype B deletion had unexpectedly low identity between haplotypes A and B. Phylogenetic analysis showed that the divergent region of sheep haplotype B is remarkably distinct from the beta-globin loci in goat and cattle but still groups with the Ruminantia. We hypothesize that this divergent ~40-kb region in haplotype B may be from an unknown ancestral ruminant and was maintained in the lineage to O. aries, but not other Bovidae, evolving independently of haplotype A. Alternatively, the ~40-kb sequence in haplotype B was more recently acquired by an ancestor of sheep from an unknown non-Bovidae ruminant, replacing part of haplotype A. Haplotype B has a lower nucleotide diversity than does haplotype A, suggesting a recent bottleneck, whereas the higher frequency of haplotype B suggests a subsequent spread through the global population of O. aries.

Genomewide association for a dominant pigmentation gene in sheep.

Most published genomewide association studies (GWAS) in sheep have investigated recessively inherited monogenic traits. The objective here was to assess the feasibility of performing GWAS for a dominant trait for which the genetic basis was already known. A total of 42 Manchega and Rasa Aragonesa sheep that segregate solid black or white coat pigmentation were genotyped using the SNP50 BeadChip. Previous analysis in Manchegas demonstrated a complete association between the pigmentation trait and alleles of the MC1R gene, setting an a priori expectation for GWAS. Multiple methods were used to identify and quantify the strength of population substructure between black and white animals, before allelic association testing was performed for 49,034 SNPs. Following correction for substructure, GWAS identified the most strongly associated SNP (s26449) was also the closest to the MC1R gene. The finding was strongly supported by the permutation tree-based random forest (RF) analysis. Importantly, GWAS identified unlinked SNP with only slightly lower p-values than for s26449. Random forest analysis indicated these were false positives, suggesting interpretation based on both approaches was beneficial. The results indicate that a combined analytical approach can be successful in studies where a modest number of animals are available and substantial population stratification exists.

Consistent divergence times and allele sharing measured from cross-species application of SNP chips developed for three domestic species.

Recent advances in technology facilitated development of large sets of genetic markers for many taxa, though most often model or domestic organisms. Cross-species application of genomic technologies may allow for rapid marker discovery in wild relatives of taxa with well-developed resources. We investigated returns from cross-species application of three commercially available SNP chips (the OvineSNP50, BovineSNP50 and EquineSNP50 BeadChips) as a function of divergence time between the domestic source species and wild target species. Across all three chips, we observed a consistent linear decrease in call rate (~1.5% per million years), while retention of polymorphisms showed an exponential decay. These results will allow researchers to predict the expected amplification rate and polymorphism of cross-species application for their taxa of interest, as well as provide a resource for estimating divergence times.

Accuracy of genotype imputation in sheep breeds.

Although genomic selection offers the prospect of improving the rate of genetic gain in meat, wool and dairy sheep breeding programs, the key constraint is likely to be the cost of genotyping. Potentially, this constraint can be overcome by genotyping selection candidates for a low density (low cost) panel of SNPs with sparse genotype coverage, imputing a much higher density of SNP genotypes using a densely genotyped reference population. These imputed genotypes would then be used with a prediction equation to produce genomic estimated breeding values. In the future, it may also be desirable to impute very dense marker genotypes or even whole genome re-sequence data from moderate density SNP panels. Such a strategy could lead to an accurate prediction of genomic estimated breeding values across breeds, for example. We used genotypes from 48 640 (50K) SNPs genotyped in four sheep breeds to investigate both the accuracy of imputation of the 50K SNPs from low density SNP panels, as well as prospects for imputing very dense or whole genome re-sequence data from the 50K SNPs (by leaving out a small number of the 50K SNPs at random). Accuracy of imputation was low if the sparse panel had less than 5000 (5K) markers. Across breeds, it was clear that the accuracy of imputing from sparse marker panels to 50K was higher if the genetic diversity within a breed was lower, such that relationships among animals in that breed were higher. The accuracy of imputation from sparse genotypes to 50K genotypes was higher when the imputation was performed within breed rather than when pooling all the data, despite the fact that the pooled reference set was much larger. For Border Leicesters, Poll Dorsets and White Suffolks, 5K sparse genotypes were sufficient to impute 50K with 80% accuracy. For Merinos, the accuracy of imputing 50K from 5K was lower at 71%, despite a large number of animals with full genotypes (2215) being used as a reference. For all breeds, the relationship of individuals to the reference explained up to 64% of the variation in accuracy of imputation, demonstrating that accuracy of imputation can be increased if sires and other ancestors of the individuals to be imputed are included in the reference population. The accuracy of imputation could also be increased if pedigree information was available and was used in tracking inheritance of large chromosome segments within families. In our study, we only considered methods of imputation based on population-wide linkage disequilibrium (largely because the pedigree for some of the populations was incomplete). Finally, in the scenarios designed to mimic imputation of high density or whole genome re-sequence data from the 50K panel, the accuracy of imputation was much higher (86-96%). This is promising, suggesting that in silico genome re-sequencing is possible in sheep if a suitable pool of key ancestors is sequenced for each breed.

A genome-wide set of SNPs detects population substructure and long range linkage disequilibrium in wild sheep.

The development of genomic resources for wild species is still in its infancy. However, cross-species utilization of technologies developed for their domestic counterparts has the potential to unlock the genomes of organisms that currently lack genomic resources. Here, we apply the OvineSNP50 BeadChip, developed for domestic sheep, to two related wild ungulate species: the bighorn sheep (Ovis canadensis) and the thinhorn sheep (Ovis dalli). Over 95% of the domestic sheep markers were successfully genotyped in a sample of fifty-two bighorn sheep while over 90% were genotyped in two thinhorn sheep. Pooling the results from both species identified 868 single-nucleotide polymorphisms (SNPs), 570 were detected in bighorn sheep, while 330 SNPs were identified in thinhorn sheep. The total panel of SNPs was able to discriminate between the two species, assign population of origin for bighorn sheep and detect known relationship classes within one population of bighorn sheep. Using an informative subset of these SNPs (n=308), we examined the extent of genome-wide linkage disequilibrium (LD) within one population of bighorn sheep and found that high levels of LD persist over 4 Mb.

Haplogroup relationships between domestic and wild sheep resolved using a mitogenome panel.

Five haplogroups have been identified in domestic sheep through global surveys of mitochondrial (mt) sequence variation, however these group classifications are often based on small fragments of the complete mtDNA sequence; partial control region or the cytochrome B gene. This study presents the complete mitogenome from representatives of each haplogroup identified in domestic sheep, plus a sample of their wild relatives. Comparison of the sequence successfully resolved the relationships between each haplogroup and provided insight into the relationship with wild sheep. The five haplogroups were characterised as branching independently, a radiation that shared a common ancestor 920,000 ± 190,000 years ago based on protein coding sequence. The utility of various mtDNA components to inform the true relationship between sheep was also examined with Bayesian, maximum likelihood and partitioned Bremmer support analyses. The control region was found to be the mtDNA component, which contributed the highest amount of support to the tree generated using the complete data set. This study provides the nucleus of a mtDNA mitogenome panel, which can be used to assess additional mitogenomes and serve as a reference set to evaluate small fragments of the mtDNA.

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.

Technical note: High fidelity of whole-genome amplified sheep (Ovis aries) deoxyribonucleic acid using a high-density single nucleotide polymorphism array-based genotyping platform.

Advances in high-throughput genotyping technologies have afforded researchers the opportunity to study ever-increasing numbers of SNP in animal genomes. However, many studies encounter difficulties in obtaining sufficient quantities of high-quality DNA for such analyses, particularly when the source biological material is limited or degraded. The recent development of in vitro whole-genome amplification approaches has permitted researchers to circumvent these challenges by increasing the amount of usable DNA in normally small-quantity samples. Here, we assess the performance of whole-genome amplification products generated from ovine genomic DNA using a high-throughput SNP genotyping platform, the newly developed Illumina ovineSNP50 BeadChip. Our results demonstrate a high genotype call rate for conventional genomic DNA and whole-genome amplified genomic DNA. The data also reveal an exceptionally high concordance rate ( > or = 99%) between the genotypes generated from whole-genome amplified products and their conventional genomic DNA counterparts. This study supports the use of whole-genome amplification as a viable solution for the analysis of high-density SNP genotypic data using compromised or limited starting material.

Evaluation of 16 loci to examine the cross-species utility of single nucleotide polymorphism arrays.

Large collections of single nucleotide polymorphisms (SNPs) have recently been identified from a number of livestock genomes. This raises the possibility that SNP arrays might be useful for analysis in related species for which few genetic markers are currently available. To address the likely success of such an approach, the aim of this study was to examine the threshold number and position of flanking mutations which act to prevent genotype calls being produced. Sequence diversity was measured across 16 loci containing SNPs known either to work successfully between species or fail between species. In pairwise comparisons between domestic and wild sheep, sequence divergence surrounding working SNP assays was significantly lower than that surrounding non-functional assays. In addition, the location of flanking mismatches tended to be closer to the target SNP in loci that failed to generate genotype calls across species. The magnitude of sequence divergence observed for both working and non-functional assays was compared with the divergence separating domestic sheep from European Mouflon, African Barbary, goat and cattle. The results suggest that the utility of SNP arrays for analysis of shared polymorphism will be restricted to closely related pairs of species. Analysis across more divergent species will, however, be successful for other objectives, such as the identification of the ancestral state of SNPs.

Re-sequencing regions of the ovine Y chromosome in domestic and wild sheep reveals novel paternal haplotypes.

The male-specific region of the ovine Y chromosome (MSY) remains poorly characterized, yet sequence variants from this region have the potential to reveal the wild progenitor of domestic sheep or examples of domestic and wild paternal introgression. The 5' promoter region of the sex-determining gene SRY was re-sequenced using a subset of wild sheep including bighorn (Ovis canadensis), thinhorn (Ovis dalli spp.), urial (Ovis vignei), argali (Ovis ammon), mouflon (Ovis musimon) and domestic sheep (Ovis aries). Seven novel SNPs (oY2-oY8) were revealed; these were polymorphic between but not within species. Re-sequencing and fragment analysis was applied to the MSY microsatellite SRYM18. It contains a complex compound repeat structure and sequencing of three novel size fragments revealed that a pentanucleotide element remained fixed, whilst a dinucleotide element displayed variability within species. Comparison of the sequence between species revealed that urial and argali sheep grouped more closely to the mouflon and domestic breeds than the pachyceriforms (bighorn and thinhorn). SNP and microsatellite data were combined to define six previously undetected haplotypes. Analysis revealed the mouflon as the only species to share a haplotype with domestic sheep, consistent with its status as a feral domesticate that has undergone male-mediated exchange with domestic animals. A comparison of the remaining wild species and domestic sheep revealed that O. aries is free from signatures of wild sheep introgression.

Mitochondrial haplotypes reveal a strong genetic structure for three Indian sheep breeds.

This survey represents the first characterization of mitochondrial DNA diversity within three breeds of Indian sheep (two strains of the Deccani breed, as well as the Bannur and Garole breeds) from different geographic regions and with divergent phenotypic characteristics. A 1061-bp fragment of the mitochondrial genome spanning the control region, a portion of the 12S rRNA gene and the complete phenyl tRNA gene, was sequenced from 73 animals and compared with the corresponding published sequence from European and Asian breeds and the European Mouflon (Ovis musimon). Analysis of all 156 sequences revealed 73 haplotypes, 52 of which belonged to the Indian breeds. The three Indian breeds had no haplotypes in common, but one Indian haplotype was shared with European and other Asian breeds. The highest nucleotide and haplotype diversity was observed in the Bannur breed (0.00355 and 0.981 respectively), while the minimum was in the Sangamneri strain of the Deccani breed (0.00167 and 0.882 respectively). All 52 Indian haplotypes belonged to mitochondrial lineage A. Therefore, these Indian sheep are distinct from other Asian and European breeds studied so far. The relationships among the haplotypes showed strong breed structure and almost no introgression among these Indian breeds, consistent with Indian sheep husbandry, which discourages genetic exchange between breeds. These results have implications for the conservation of India's ovine biodiversity and suggest a common origin for the breeds investigated.

The effect of genetic variation of the retinoic acid receptor-related orphan receptor C gene on fatness in cattle.

Genotypes at the retinoic acid receptor-related orphan receptor C (RORC) gene were associated with fatness in 1750 cattle. Ten SNPs were genotyped in RORC and the adjacent gene leucine-rich repeat neuronal 6D (LRRN6D) to map the QTL, 7 of which are in a 4.2-kb sequence around the ligand-binding domain of the RORC gene. Of the 29 inferred haplotypes for these SNPs, 2 have a combined frequency of 54.6% while the top 5 haplotypes have a combined frequency of 85.3%. The average D' value of linkage disequilibrium was 0.92 although the average r2 was a low 0.18. The RORC:g.3290T>G SNP had the strongest association with marbling. The inferred haplotypes were significantly associated with marbling and the difference between the most divergent haplotypes was 0.35 sigma(p) of marbling and 0.28 sigma(p) of rump fat, explaining the previously reported QTL effect. cDNA for RORC were sequenced and 2 new alternative transcripts were found. Fetal tissue shows 40 times greater transcription of RORC than adult tissue. The highest expression in fetal tissue was found in liver and kidney, but in adults the longissimus muscle had the greatest expression of the tissues tested.

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.

Sequence diversity and rates of molecular evolution between sheep and cattle genes.

Experiments that aim to identify genes of importance in sheep are currently inhibited by a paucity of genomic resources. One approach, therefore, is to exploit the wealth of data and associated capabilities becoming available for the bovine genome. Cross-species application of microarrays and comparative sequencing to identify single nucleotide polymorphisms are two possibilities; however, both are dependant on the level of nucleotide sequence similarity between the two species. This study used 120 gene orthologues consisting of over 60 kb of aligned sequence to estimate the gene diversity between cattle and sheep. Less than 3% of protein-coding nucleotide positions were found to be different, indicating that the prospect for successfully using cross-species strategies is high. Substitution at synonymous sites ranged between 6.9 and 7.7% (+/- 0.3%), and was higher than at non-synonymous sites (1.4-1.7 +/- 0.1%). The relative rate test was used to determine whether the observed mutation rates were constant between the two lineages. While the rate at synonymous sites appeared constant, the rate at non-synonymous sites was significantly higher within the caprinae lineage (sheep) when compared with bovinae (cattle; chi2 = 10.03; d.f. = 1, P < 0.01). This is the first demonstration that variable rates of molecular evolution may be present within the family Bovidae.

Mitochondrial sequence reveals high levels of gene flow between breeds of domestic sheep from Asia and Europe.

Sequence variation present within the mitochondrial genome was used to investigate genetic diversity within sheep breeds from Asia and Europe. Comparison of 2027 bp of sequence from 121 animals revealed 44 phylogenetically informative nucleotide positions and a single insertion/deletion. A total of 57 haplotypes were observed which formed two distinct clades. Type A haplotypes were found in breeds from Asia (India, Indonesia, Mongolia, and Tibet), while type B haplotypes were observed at the highest frequency in breeds sourced from Europe (nine breeds from Austria, Aland, Finland, Spain, and northwestern Russia). The distribution of haplotypes indicates sheep appear to have the weakest population structure and the highest rate of intercontinental dispersal of any domestic animal reported to date. Only 2.7% of the sequence variation observed was partitioned between continents, which is lower than both goat (approximately 10%) and cattle (approximately 50%). Diagnostic restriction fragment length polymorphism polymerase chain reaction (RFLP-PCR) tests which distinguish type A and B haplotypes were used to test an additional 223 animals from 17 breeds of European and Asian origin. A mixture of the two lineages was found in every breed except Suffolk and the Indian Garole, indicating introgression has played a major part during breed development and subsequent selection.

Nucleotide diversity on the ovine Y chromosome.

To investigate the impact of male-mediated introgression during the evolution of sheep breeds, a sequencing approach was used to identify single nucleotide polymorphisms (SNPs) from the male-specific region of the ovine Y chromosome (MSY). A total of 4380 bp, which comprised nine fragments from five MSY genes was sequenced within a panel of 14 males from seven breeds. Sequence alignment identified a single segregating site, an A/G SNP located approximately 1685 bp upstream of the ovine SRY gene. The resulting estimation of nucleotide diversity (piY = 0.90 +/- 0.50 x 10(-4)) falls towards the lower end of estimates from other species. This was compared with the nucleotide diversity estimated from the autosomal component of the genome. Sequence analysis of 2933 bp amplified from eight autosomal genes revealed a nucleotide diversity (piA = 2.15 +/- 0.27 x 10(-3)) higher than previously reported for sheep. Following adjustment for the contrasting influence of effective population size and a male biased mutation rate, comparison revealed that approximately 10% of the expected nucleotide diversity is present on the ovine Y chromosome.

Canine models of ocular disease: outcross breedings define a dominant disorder present in the English mastiff and bull mastiff dog breeds.

Progressive retinal atrophies (PRA) are a heterogeneous group of inherited eye diseases common to both dogs and man. Over 100 individual canine breeds display some sort of retinal degeneration, making the dog an extremely valuable resource both for finding the genetic determinants of inherited blindness and for developing naturally occurring animal models that mimic human disease. Progressive retinal atrophies within the English mastiff displayed an ambiguous mode of inheritance. By conducting outcross matings between affected English mastiffs and normal animals from other breeds, the mode of inheritance was confirmed as dominant. This directed candidate gene analysis and led to identification of two synonymous mutations and one nonsynonymous mutation within the canine rhodopsin gene. The nonsynonymous mutation (T4R) is the cause of PRA in the English mastiff, and a test was developed to investigate its presence in 17 additional breeds. Testing of PRA-affected animals from 16 breeds revealed that none carry the T4R mutation, indicating a different cause of PRA. Analysis of two affected bull mastiffs revealed one heterozygote (+/T4R) and one homozygous normal individual (+/+). These findings suggest that the genetic origin of PRA is often breed specific and underline the value of outcross mating to circumvent problems that act to mask the mode of inheritance.