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.

The 'heritability' of domestication and its functional partitioning in the pig.

We propose to estimate the proportion of variance explained by regression on genome-wide markers (or genomic heritability) when wild/domestic status is considered the phenotype of interest. This approach differs from the standard Fst in that it can accommodate genetic similarity between individuals in a general form. We apply this strategy to complete genome data from 47 wild and domestic pigs from Asia and Europe. When we partitioned the total genomic variance into components associated to subsets of single nucleotide polymorphisms (SNPs) defined in terms of their annotation, we found that potentially deleterious non-synonymous mutations (9566 SNPs) explained as much genetic variance as the whole set of 25 million SNPs. This suggests that domestication may have affected protein sequence to a larger extent than regulatory or other kinds of mutations. A pathway-guided analysis revealed ovarian steroidogenesis and leptin signaling as highly relevant in domestication. The genomic regression approach proposed in this study revealed molecular processes not apparent through typical differentiation statistics. We propose that at least some of these processes are likely new discoveries because domestication is a dynamic process of genetic selection, which may not be completely characterized by a static metric like Fst. Nevertheless, and despite some particularly influential mutation types or pathways, our analyses tend to rule out a simplistic genetic basis for the domestication process: neither a single pathway nor a unique set of SNPs can explain the process as a whole.

Diversity and linkage disequilibrium in farmed Tasmanian Atlantic salmon.

Farmed Atlantic salmon (Salmo salar) is a globally important production species, including in Australia where breeding and selection has been in progress since the 1960s. The recent development of SNP genotyping platforms means genome-wide association and genomic prediction can now be implemented to speed genetic gain. As a precursor, this study collected genotypes at 218 132 SNPs in 777 fish from a Tasmanian breeding population to assess levels of genetic diversity, the strength of linkage disequilibrium (LD) and imputation accuracy. Genetic diversity in Tasmanian Atlantic salmon was lower than observed within European populations when compared using four diversity metrics. The distribution of allele frequencies also showed a clear difference, with the Tasmanian animals carrying an excess of low minor allele frequency variants. The strength of observed LD was high at short distances (<25 kb) and remained above background for marker pairs separated by large chromosomal distances (hundreds of kb), in sharp contrast to the European Atlantic salmon tested. Genotypes were used to evaluate the accuracy of imputation from low density (0.5 to 5 K) up to increased density SNP sets (78 K). This revealed high imputation accuracies (0.89-0.97), suggesting that the use of low density SNP sets will be a successful approach for genomic prediction in this population. The long-range LD, comparatively low genetic diversity and high imputation accuracy in Tasmanian salmon is consistent with known aspects of their population history, which involved a small founding population and an absence of subsequent introgression. The findings of this study represent an important first step towards the design of methods to apply genomics in this economically important population.

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.

Brachygnathia, cardiomegaly and renal hypoplasia syndrome (BCRHS) in Merino sheep maps to a 1.1-megabase region on ovine chromosome OAR2.

A genome scan was conducted to map the autosomal recessive lethal disorder brachygnathia, cardiomegaly and renal hypoplasia syndrome (BCRHS) in Poll Merino sheep. The scan involved 10 affected and 27 unaffected animals from a single Poll Merino/Merino sheep flock, which were genotyped with the Illumina Ovine SNP50 BeadChip. Association and homozygosity mapping analyses located the disorder in a region comprising 20 consecutive SNPs spanning 1.1 Mb towards the distal end of chromosome OAR2. All affected animals and none of the unaffected animals were homozygous for the associated haplotype in this region. These results provide the basis for identifying the causative mutation(s) and should enable the development of a DNA test to identify carriers in the Poll Merino sheep population. Understanding the molecular control of BCRHS may provide insight into the fundamental genetic control and regulation of the affected organ systems.

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.

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.

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.

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.

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.

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.

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.

A frameshift mutation in MC1R and a high frequency of somatic reversions cause black spotting in pigs.

Black spotting on a red or white background in pigs is determined by the E(P) allele at the MC1R/Extension locus. A previous comparison of partial MC1R sequences revealed that E(P) shares a missense mutation (D121N) with the E(D2) allele for dominant black color. Sequence analysis of the entire coding region now reveals a second mutation in the form of a 2-bp insertion at codon 23 (nt67insCC). This mutation expands a tract of six C nucleotides to eight and introduces a premature stop codon at position 56. This frameshift mutation is expected to cause a recessive red color, which was in fact observed in some breeds with the E(P) allele present (Tamworth and Hereford). RT-PCR analyses were conducted using skin samples taken from both spotted and background areas of spotted pigs. The background red area had transcript only from the mutant nt67insCC MC1R allele, whereas the black spot also contained a transcript without the 2-bp insertion. This indicates that black spots are due to somatic reversion events that restore the frame and MC1R function. The phenotypic expression of the E(P) allele is highly variable and the associated coat color ranges from red, red with black spots, white with black spots, to almost completely solid black. In several breeds of pigs the phenotypic manifestation of this allele has been modified by selection for or against black spots.

The origin of the domestic pig: independent domestication and subsequent introgression.

The domestic pig originates from the Eurasian wild boar (Sus scrofa). We have sequenced mitochondrial DNA and nuclear genes from wild and domestic pigs from Asia and Europe. Clear evidence was obtained for domestication to have occurred independently from wild boar subspecies in Europe and Asia. The time since divergence of the ancestral forms was estimated at approximately 500,000 years, well before domestication approximately 9,000 years ago. Historical records indicate that Asian pigs were introduced into Europe during the 18th and early 19th centuries. We found molecular evidence for this introgression and the data indicated a hybrid origin of some major "European" pig breeds. The study is an advance in pig genetics and has important implications for the maintenance and utilization of genetic diversity in this livestock species.

Melanocortin receptor 1 (MC1R) mutations and coat color in pigs.

The melanocortin receptor 1 (MC1R) plays a central role in regulation of eumelanin (black/brown) and phaeomelanin (red/yellow) synthesis within the mammalian melanocyte and is encoded by the classical Extension (E) coat color locus. Sequence analysis of MC1R from seven porcine breeds revealed a total of four allelic variants corresponding to five different E alleles. The European wild boar possessed a unique MC1R allele that we believe is required for the expression of a wild-type coat color. Two different MC1R alleles were associated with the dominant black color in pigs. MC1R*2 was found in European Large Black and Chinese Meishan pigs and exhibited two missense mutations compared with the wild-type sequence. Comparative data strongly suggest that one of these, L99P, may form a constitutively active receptor. MC1R*3 was associated with the black color in the Hampshire breed and involved a single missense mutation D121N. This same MC1R variant was also associated with EP, which results in black spots on a white or red background. Two different missense mutations were identified in recessive red (e/e) animals. One of these, A240T, occurs at a highly conserved position, making it a strong candidate for disruption of receptor function.

Enrichment of microsatellites from the citrus genome using biotinylated oligonucleotide sequences bound to streptavidin-coated magnetic particles.

A method is described for the rapid isolation of microsatellite sequences using a biotin-labeled oligonucleotide attached to streptavidin-coated magnetic particles. The oligonucleotide "hook" in solution hybridizes to complementary single-stranded lengths of genomic DNA onto which have been engineered specific PCR priming sites. The final product is an enriched library of microsatellites of defined sequences. The method is applicable to any genome and in principle is adaptable to the rapid isolation of both repetitive as well as genic sequences. It is illustrated by the isolation of trinucleotide repeat (TAA)n sequences from the citrus genome.

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.