Fine-scale population structure analysis of seven local Swiss sheep breeds using genome-wide SNP data

As part of the global sheep Hapmap project, 24 individuals from each of seven indigenous Swiss sheep breeds (Bundner Oberländer sheep (BOS), Engadine Red sheep (ERS), Swiss Black-Brown Mountain sheep (SBS), Swiss Mirror sheep (SMS), Swiss White Alpine (SWA) sheep, Valais Blacknose sheep (VBS) and Valais Red sheep (VRS)), were genotyped using Illumina's Ovine SNP50 BeadChip. In total, 167 animals were subjected to a detailed analysis for genetic diversity using 45193 informative single nucleotide polymorphisms. The results of the phylogenetic analyses supported the known proximity between populations such as VBS and VRS or SMS and SWA. Average genomic relatedness within a breed was found to be 12 percent (BOS), 5 percent (ERS), 9 percent (SBS), 10 percent (SMS), 9 percent (SWA), 12 percent (VBS) and 20 percent (VRS). Furthermore, genomic relationships between breeds were found for single individuals from SWA and SMS, VRS and VBS as well as VRS and BOS. In addition, seven out of 40 indicated parent-offspring pairs could not be confirmed. These results were further supported by results from the genome-wide population cluster analysis. This study provides a better understanding of fine-scale population structures within and between Swiss sheep breeds. This relevant information will help to increase the conservation activities of the local Swiss sheep breed

A Genome-Wide Association Study Reveals Loci Influencing Height and Other Conformation Traits in Horses

The molecular analysis of genes influencing human height has been notoriously difficult. Genome-wide association studies (GWAS) for height in humans based on tens of thousands to hundreds of thousands of samples so far revealed ∼200 loci for human height explaining only 20% of the heritability. In domestic animals isolated populations with a greatly reduced genetic heterogeneity facilitate a more efficient analysis of complex traits. We performed a genome-wide association study on 1,077 Franches-Montagnes (FM) horses using ∼40,000 SNPs. Our study revealed two QTL for height at withers on chromosomes 3 and 9. The association signal on chromosome 3 is close to the LCORL/NCAPG genes. The association signal on chromosome 9 is close to the ZFAT gene. Both loci have already been shown to influence height in humans. Interestingly, there are very large intergenic regions at the association signals. The two detected QTL together explain ∼18.2% of the heritable variation of height in horses. However, another large fraction of the variance for height in horses results from ECA 1 (11.0%), although the association analysis did not reveal significantly associated SNPs on this chromosome. The QTL region on ECA 3 associated with height at withers was also significantly associated with wither height, conformation of legs, ventral border of mandible, correctness of gaits, and expression of the head. The region on ECA 9 associated with height at withers was also associated with wither height, length of croup and length of back. In addition to these two QTL regions on ECA 3 and ECA 9 we detected another QTL on ECA 6 for correctness of gaits. Our study highlights the value of domestic animal populations for the genetic analysis of complex traits

Strong signatures of selection in the domestic pig genome

Domestication of wild boar (Sus scrofa) and subsequent selection have resulted in dramatic phenotypic changes in domestic pigs for a number of traits, including behavior, body composition, reproduction, and coat color. Here we have used whole-genome resequencing to reveal some of the loci that underlie phenotypic evolution in European domestic pigs. Selective sweep analyses revealed strong signatures of selection at three loci harboring quantitative trait loci that explain a considerable part of one of the most characteristic morphological changes in the domestic pig—the elongation of the back and an increased number of vertebrae. The three loci were associated with the NR6A1, PLAG1, and LCORL genes. The latter two have repeatedly been associated with loci controlling stature in other domestic animals and in humans. Most European domestic pigs are homozygous for the same haplotype at these three loci. We found an excess of derived nonsynonymous substitutions in domestic pigs, most likely reflecting both positive selection and relaxed purifying selection after domestication. Our analysis of structural variation revealed four duplications at the KIT locus that were exclusively present in white or white-spotted pigs, carrying the Dominant white, Patch, or Belt alleles. This discovery illustrates how structural changes have contributed to rapid phenotypic evolution in domestic animals and how alleles in domestic animals may evolve by the accumulation of multiple causative mutations as a response to strong directional selection

Runs of homozygosity and signatures of selection: a comparison among eight local Swiss sheep breeds.

A dataset consisting of 787 animals with high-density SNP chip genotypes (346 774 SNPs) and 939 animals with medium-density SNP chip genotypes (33 828 SNPs) from eight indigenous Swiss sheep breeds was analyzed to characterize population structure, quantify genomic inbreeding based on runs of homozygosity and identify selection signatures. In concordance with the recent known history of these breeds, the highest genetic diversity was observed in Engadine Red sheep and the lowest in Valais Blacknose sheep. Correlation between F and F was around 0.50 and thereby lower than that found in similar studies in cattle. Mean F estimates from medium-density data and HD data were highly correlated (0.95). Signatures of selection and candidate gene analysis revealed that the most prominent signatures of selection were found in the proximity of genes associated with body size (NCAPG, LCORL, LAP3, SPP1, PLAG1, ALOX12, TP53), litter size (SPP1), milk production (ABCG2, SPP1), coat color (KIT, ASIP, TBX3) and horn status (RXFP2). For the Valais Blacknose sheep, the private signatures in proximity of genes/QTL influencing body size, coat color and fatty acid composition were confirmed based on runs of homozygosity analysis. These private signatures underline the genetic uniqueness of the Valais Blacknose sheep breed. In conclusion, we identified differences in the genetic make-up of Swiss sheep breeds, and we present relevant candidate genes responsible for breed differentiation in locally adapted breeds

Genome-wide association studies based on sequence-derived genotypes reveal new QTL associated with conformation and performance traits in the Franches-Montagnes horse breed

To identify novel quantitative trait loci (QTL) within horses, we performed genome-wide association studies (GWAS) based on sequence-level genotypes for conformation and performance traits in the Franches-Montagnes (FM) horse breed. Sequence-level genotypes of FM horses were derived by re-sequencing 30 key founders and imputing 50K data of genotyped horses. In total, we included 1077 FM horses genotyped for ~4 million SNPs and their respective de-regressed breeding values of the traits in the analysis. Based on this dataset, we identified a total of 14 QTL associated with 18 conformation traits and one performance trait. Therefore, our results suggest that the application of sequence-derived genotypes increases the power to identify novel QTL which were not identified previously based on 50K SNP chip data

A complex structural variant at the KIT locus in cattle with the Pinzgauer spotting pattern.

A specific white spotting phenotype, termed finching or line-backed spotting, is known for all Pinzgauer cattle and occurs occasionally in Tux-Zillertaler cattle, two Austrian breeds. The so-called Pinzgauer spotting is inherited as an autosomal incompletely dominant trait. A genome-wide association study using 27 white spotted and 16 solid-coloured Tux-Zillertaler cattle, based on 777k SNP data, revealed a strong signal on chromosome 6 at the KIT locus. Haplotype analyses defined a critical interval of 122 kb downstream of the KIT coding region. Whole-genome sequencing of a Pinzgauer cattle and comparison to 338 control genomes revealed a complex structural variant consisting of a 9.4-kb deletion and an inversely inserted duplication of 1.5 kb fused to a 310-kb duplicated segment from chromosome 4. A diagnostic PCR was developed for straightforward genotyping of carriers for this structural variant (KIT ) and confirmed that the variant allele was present in all Pinzgauer and most of the white spotted Tux-Zillertaler cattle. In addition, we detected the variant in all Slovenian Cika, British Gloucester and Spanish Berrenda en negro cattle with similar spotting patterns. Interestingly, the KIT variant occurs in some white spotted animals of the Swiss breeds Evolèner and Eringer. The introgression of the KIT variant confirms admixture and the reported historical relationship of these short-headed breeds with Austrian Tux-Zillertaler and suggests a mutation event, occurring before breed formation

Genome-wide association study in Brown Swiss for udder traits based on sequence data

Identification of QTL, especially of causative variants within QTL is still challenging. Higher SNP densities aid the identification and fine-mapping of QTL. Based on imputed sequence data we performed GWAS for udder traits in Brown Swiss cattle. The GWAS was performed using a mixed-model approach with deregressed breeding values as phenotypes. The traits investigated included: udder depth (UD), fore udder attachment (FUA), rear udder width (RUW), rear udder height (RUH), fore udder length (FUL), and central ligament (CL). We found significant associations on BTA 3 (UD, FUA), BTA 5 (UD), BTA 17 (FUL, RUW, CL) and BTA 20 (FUA). A single gene was located in the significantly associated regions on BTA 5 (ABCC9) and BTA 20 (HCN1). The region on BTA 17 spans almost 3 Mb and includes 74 genes (maximal region for all the traits combined) and the region on BTA 3 includes 91 genes across 3 Mb. We also looked for associated missense variants in these intervals. Neither for ABCC9 nor for HCN1 we could identify such a variant. On BTA 17 we identified 2 missense variants that were significantly associated with CL. On BTA 3, 11 missense variants were significantly associated with UD and/or FUA. The advantages of using imputed sequence data compared to SNP chip genotypes are mainly through the inclusion of potential causative variants

Genome-wide association study in Brown Swiss for fertility traits based on sequence data

Fertility in dairy cattle has declined continuously during the last years. In Switzerland low fertility is one of the most important culling reason. Therefore this trait complex is of great economic importance. But fertility traits are functional traits with low heritabilities. We performed GWAS in Brown Swiss and Holstein cattle for the fertility traits: non-return rate after 56 days in heifers (NRH) and cows (NRC); days between first and last insemination in heifers (FLIH) and cows (FLIC); and days to first service (DFS). We used deregressed breeding values as input phenotypes and imputed whole-genome sequence dosages as genotypes. In Holstein (2,336 to 3,257 individuals) we could not identify any significant region for any of the investigated traits. For Brown Swiss cattle (1,392 to 3,278 individuals) however we could identify a QTL on BTA 17 for NRH and FLIH. This locus also showed a suggestive association with NRC and FLIC. The region significantly associated to NRH and FLIH comprises 27 genes. Among the associated variants two non-synonymous variants in the genes GAS2L and ASCC2 could be identified. These two variants could be potentially causative ariants