Genome‐wide association study to identify genomic regions affecting prolificacy in Lori‐Bakhtiari sheep

Several causative mutations in candidate genes affecting prolificacy have been detected in various sheep breeds. A genome‐wide association study was performed on estimated breeding values for litter size in Lori‐Bakhtiari sheep. Prolific ewes with twinning records and others with only singleton records were genotyped using the medium‐density Illumina Ovine SNP50 array. Four single nucleotide polymorphisms (SNPs) associated with litter size were identified on chromosomes 3, 6 and 22. The region on sheep chromosome 3 between 75 739 167 and 75 745 152 bp included two significant SNPs (s52383.1 and OAR3_80038014_X.1) in high linkage disequilibrium with each other. The region that surrounds these SNPs contains a novel putative candidate gene: luteinizing hormone/choriogonadotropin receptor (LHCGR), known to be involved in ovarian steroidogenesis and organism‐specific biosystem pathways in sheep. Known prolificacy genes BMPR1B, BMP15 and GDF9 were not associated with litter size in Lori‐Bakhtiari sheep, suggesting that other biological mechanisms could be responsible for the trait's variation in this breed.     

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^PINZ) 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^PINZ variant occurs in some white spotted animals of the Swiss breeds Evolèner and Eringer. The introgression of the KIT^PINZ 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.     

Searching new signals for production traits through gene‐based association analysis in three Italian cattle breeds

Genome‐wide association studies (GWAS) have been widely applied to disentangle the genetic basis of complex traits. In cattle breeds, classical GWAS approaches with medium‐density marker panels are far from conclusive, especially for complex traits. This is due to the intrinsic limitations of GWAS and the assumptions that are made to step from the association signals to the functional variations. Here, we applied a gene‐based strategy to prioritize genotype–phenotype associations found for milk production and quality traits with classical approaches in three Italian dairy cattle breeds with different sample sizes (Italian Brown n = 745; Italian Holstein n = 2058; Italian Simmental n = 477). Although classical regression on single markers revealed only a single genome‐wide significant genotype–phenotype association, for Italian Holstein, the gene‐based approach identified specific genes in each breed that are associated with milk physiology and mammary gland development. As no standard method has yet been established to step from variation to functional units (i.e., genes), the strategy proposed here may contribute to revealing new genes that play significant roles in complex traits, such as those investigated here, amplifying low association signals using a gene‐centric approach.     

Genome‐wide association analysis identifies the genetic basis of fat deposition in the tails of sheep (Ovis aries)

Fat‐tailed sheep (Ovis aries) can survive in harsh environments and satisfy human's intake of dietary fat. However, the animals require more feed, which increases the cost of farming. Thus, most farmers currently prefer thin‐tailed, short‐tailed or docked sheep. To date, the molecular mechanism of the formation of fat tails in sheep has not been completely elucidated. Here, we conducted a genome‐wide association study using phenotypes and genotypes (the Ovine Infinium HD SNP BeadChip genotype data) of two breeds of contrasting tail types (78 Small‐tailed and 78 Large‐tailed Han sheep breeds) to identify functional genes and variants associated with fat deposition. We identified four significantly (rs416433540, rs409848439, rs408118325 and rs402128848) and three approximately associated autosomal SNPs (rs401248376, rs402445895 and rs416201901). Gene annotation indicated that the surrounding genes (CREB1, STEAP4, CTBP1 and RIP140, also known as NRIP1) function in lipid storage or fat cell regulation. Furthermore, through an X‐chromosome‐wide association analysis, we detected significantly associated SNPs in the OARX: 88–89 Mb region, which could be a strong candidate genomic region for fat deposition in tails of sheep. Our results represent a new genomic resource for sheep genetics and breeding. In addition, the findings provide novel insights into genetic mechanisms of fat deposition in the tail of sheep and other mammals.     

Screening for Booroola (FecB) and Galway (FecX) mutations in Indian sheep

This study reports the status of the Booroola (FecB) and Galway (FecX^G) mutations in Indian sheep breeds. The Kendrapada sheep (n = 46) was genotyped for the presence of FecB and FecX^G mutations, while the Garole (n = 34), Malpura (n = 30), and Decanni sheep (n = 15) for the FecX^G mutation. The FecB and FecX^G genotyping was carried out by forced restriction fragment length polymorphism PCR technique. In the present study, FecB mutation was discovered in the Kendrapada sheep of Orissa, which is now the second prolific sheep of India after the Garole. Out of 46 individuals of Kendrapada sheep, 26 were homozygous (BB), 15 heterozygous (B+) and 5 non-carriers (++) for the FecB mutation. The frequency of the FecB allele in this sample was about 0.73. Results indicated that the frequency of the FecB mutation is high, but the gene is not fixed in the population as reported in Garole sheep. None of sheep breeds carried the FecX^G mutation. The discovery of the FecB mutation in Kendrapada sheep will facilitate the use of FecB allele in improving the prolificacy of non-prolific sheep breeds of India.     

Deletion variant near ZNF389 is associated with control of ovine lentivirus in multiple sheep flocks

Ovine lentivirus (OvLV) is a macrophage‐tropic lentivirus found in many countries that causes interstitial pneumonia, mastitis, arthritis and cachexia in sheep. There is no preventive vaccine and no cure, but breed differences suggest marker‐assisted selective breeding might improve odds of infection and control of OvLV post‐infection. Although variants in TMEM154 have consistent association with odds of infection, no variant in any gene has been associated with host control of OvLV post‐infection in multiple animal sets. Proviral concentration is a live‐animal diagnostic measure of OvLV control post‐infection related to severity of OvLV‐induced lesions. A recent genome‐wide association study identified a region including four zinc finger genes associated with proviral concentration in one Rambouillet flock. To refine this region, we tested additional variants and identified a small insertion/deletion variant near ZNF389 that showed consistent association with proviral concentration in three animal sets (P < 0.05). These animal sets contained Rambouillet, Polypay and crossbred sheep from multiple locations and management conditions. Strikingly, one flock had exceptionally high prevalence (>87%, including yearlings) and mean proviral concentration (>950 copies/μg), possibly due to needle sharing. The best estimate of proviral concentration by genotype, obtained from all 1310 OvLV‐positive animals tested, showed insertion homozygotes had less than half the proviral concentration of other genotypes (P < 0.0001). Future work will test additional breeds, management conditions and viral subtypes, and identify functional properties of the haplotype this deletion variant tracks. To our knowledge, this is the first genetic variant consistently associated with host control of OvLV post‐infection in multiple sheep flocks.     

Mapping the four‐horned locus and testing the polled locus in three Chinese sheep breeds

Four‐horned sheep are an ideal animal model for illuminating the genetic basis of horn development. The objective of this study was to locate the genetic region responsible for the four‐horned phenotype and to verify a previously reported polled locus in three Chinese breeds. A genome‐wide association study (GWAS) was performed using 34 two‐horned and 32 four‐horned sheep from three Chinese indigenous breeds: Altay, Mongolian and Sishui Fur sheep. The top two significant single nucleotide polymorphisms (SNPs) associated with the four‐horned phenotype were both located in a region spanning positions 132.6 to 132.7 Mb on sheep chromosome 2. Similar locations for the four‐horned trait were previously identified in Jacob, Navajo‐Churro, Damara and Sishui Fur sheep, suggesting a common genetic component underlying the four‐horned phenotype. The two identified SNPs were both downstream of the metaxin 2 (MTX2) gene and the HOXD gene cluster. For the top SNP—OAR2:g.132619300G>A—the strong associations of the AA and AG genotypes with the four‐horned phenotype and the GG genotype with the two‐horned phenotype indicated the dominant inheritance of the four‐horned trait. No significant SNPs for the polled phenotype were identified in the GWAS analysis, and a PCR analysis for the detection of the 1.8‐kb insertion associated with polled sheep in other breeds failed to verify the association with polledness in the three Chinese breeds. This study supports the hypothesis that two different loci are responsible for horn existence and number. This study contributes to the understanding of the molecular regulation of horn development and enriches the knowledge of qualitative traits in domestic animals.     

Identification of copy number variations in three Chinese horse breeds using 70K single nucleotide polymorphism BeadChip array

Copy number variation (CNV), an essential form of genetic variation, has been increasingly recognized as one promising genetic marker in the analysis of animal genomes. Here, we used the Equine 70K single nucleotide polymorphism genotyping array for the genome‐wide detection of CNVs in 96 horses from three diverse Chinese breeds: Debao pony (DB), Mongolian horse (MG) and Yili horse (YL). A total of 287 CNVs were determined and merged into 122 CNV regions (CNVRs) ranging from 199 bp to 2344 kb in size and distributed in a heterogeneous manner on chromosomes. These CNVRs were integrated with seven existing reports to generate a composite genome‐wide dataset of 1558 equine CNVRs, revealing 69 (56.6%) novel CNVRs. The majority (69.7%) of the 122 CNVRs overlapped with 438 genes, whereas 30.3% were located in intergenic regions. Most of these genes were associated with common CNVRs, which were shared by divergent horse breeds. As many as 60, 42 and 91 genes overlapping with the breed‐specific ss were identified in DB, MG and YL respectively. Among these genes, FGF11, SPEM1, PPARG, CIDEB, HIVEP1 and GALR may have potential relevance to breed‐specific traits. These findings provide valuable information for understanding the equine genome and facilitating association studies of economically important traits with equine CNVRs in the future.     

Polyceraty (multi‐horns) in Damara sheep maps to ovine chromosome 2

Polyceraty (presence of multiple horns) is rare in modern day ungulates. Although not found in wild sheep, polyceraty does occur in a small number of domestic sheep breeds covering a wide geographical region. Damara are fat‐tailed hair sheep, from the south‐western region of Africa, which display polyceraty, with horn number ranging from zero to four. We conducted a genome‐wide association study for horn number with 43 Damara genotyped with 606 006 SNP markers. The analysis revealed a region with multiple significant SNPs on ovine chromosome 2, in a location different from the mutation for polled in sheep on chromosome 10. The causal mutation for polyceraty was not identified; however, the region associated with polyceraty spans nine HOXD genes, which are critical in embryonic development of appendages. Mutations in HOXD genes are implicated in polydactly phenotypes in mice and humans. There was no evidence for epistatic interactions contributing to polyceraty. This is the first report on the genetic mechanisms underlying polyceraty in the under‐studied Damara.     

Genome‐wide association reveals the locus responsible for four‐horned ruminant

Phenotypic variability in horn characteristics, such as their size, number and shape, offers the opportunity to elucidate the molecular basis of horn development. The objective of this study was to map the genetic determinant controlling the production of four horns in two breeds, Jacob sheep and Navajo‐Churro, and examine whether an eyelid abnormality occurring in the same populations is related. Genome‐wide association mapping was performed using 125 animals from the two breeds that contain two‐ and four‐horned individuals. A case–control design analysis of 570 712 SNPs genotyped with the ovine HD SNP Beadchip revealed a strong association signal on sheep chromosome 2. The 10 most strongly associated SNPs were all located in a region spanning Mb positions 131.9–132.6, indicating the genetic architecture underpinning the production of four horns is likely to involve a single gene. The closest genes to the most strongly associated marker (OAR2_132568092) were MTX2 and the HOXD cluster, located approximately 93 Kb and 251 Kb upstream respectively. The occurrence of an eyelid malformation across both breeds was restricted to polled animals and those carrying more than two horns. This suggests the eyelid abnormality may be associated with departures from the normal developmental production of two‐horned animals and that the two conditions are developmentally linked. This study demonstrated the presence of separate loci responsible for the polled and four‐horned phenotypes in sheep and advanced our understanding of the complexity that underpins horn morphology in ruminants.     

Genetic factors controlling wool shedding in a composite Easycare sheep flock

Historically, sheep have been selectively bred for desirable traits including wool characteristics. However, recent moves towards extensive farming and reduced farm labour have seen a renewed interest in Easycare breeds. The aim of this study was to quantify the underlying genetic architecture of wool shedding in an Easycare flock. Wool shedding scores were collected from 565 pedigreed commercial Easycare sheep from 2002 to 2010. The wool scoring system was based on a 10‐point (0–9) scale, with score 0 for animals retaining full fleece and 9 for those completely shedding. DNA was sampled from 200 animals of which 48 with extreme phenotypes were genotyped using a 50‐k SNP chip. Three genetic analyses were performed: heritability analysis, complex segregation analysis to test for a major gene hypothesis and a genome‐wide association study to map regions in the genome affecting the trait. Phenotypes were treated as a continuous or binary variable and categories. High estimates of heritability (0.80 when treated as a continuous, 0.65–0.75 as binary and 0.75 as categories) for shedding were obtained from linear mixed model analyses. Complex segregation analysis gave similar estimates (0.80 ± 0.06) to those above with additional evidence for a major gene with dominance effects. Mixed model association analyses identified four significant (P < 0.05) SNPs. Further analyses of these four SNPs in all 200 animals revealed that one of the SNPs displayed dominance effects similar to those obtained from the complex segregation analyses. In summary, we found strong genetic control for wool shedding, demonstrated the possibility of a single putative dominant gene controlling this trait and identified four SNPs that may be in partial linkage disequilibrium with gene(s) controlling shedding.     

Genome‐wide association analysis for quantitative trait loci influencing Warner–Bratzler shear force in five taurine cattle breeds

We performed a genome‐wide association study for Warner–Bratzler shear force (WBSF), a measure of meat tenderness, by genotyping 3360 animals from five breeds with 54 790 BovineSNP50 and 96 putative single‐nucleotide polymorphisms (SNPs) within μ‐calpain [HUGO nomenclature calpain 1, (mu/I) large subunit; CAPN1] and calpastatin (CAST). Within‐ and across‐breed analyses estimated SNP allele substitution effects (ASEs) by genomic best linear unbiased prediction (GBLUP) and variance components by restricted maximum likelihood under an animal model incorporating a genomic relationship matrix. GBLUP estimates of ASEs from the across‐breed analysis were moderately correlated (0.31–0.66) with those from the individual within‐breed analyses, indicating that prediction equations for molecular estimates of breeding value developed from across‐breed analyses should be effective for genomic selection within breeds. We identified 79 genomic regions associated with WBSF in at least three breeds, but only eight were detected in all five breeds, suggesting that the within‐breed analyses were underpowered, that different quantitative trait loci (QTL) underlie variation between breeds or that the BovineSNP50 SNP density is insufficient to detect common QTL among breeds. In the across‐breed analysis, CAPN1 was followed by CAST as the most strongly associated WBSF QTL genome‐wide, and associations with both were detected in all five breeds. We show that none of the four commercialized CAST and CAPN1 SNP diagnostics are causal for associations with WBSF, and we putatively fine‐map the CAPN1 causal mutation to a 4581‐bp region. We estimate that variation in CAST and CAPN1 explains 1.02 and 1.85% of the phenotypic variation in WBSF respectively.     

Progressive retinal atrophy in Shetland sheepdog is associated with a mutation in the CNGA1 gene

Progressive retinal atrophy (PRA) is the collective name of a class of hereditary retinal dystrophies in the dog and is often described as the equivalent of retinitis pigmentosa in humans. PRA is characterized by visual impairment due to degeneration of the photoreceptors in the retina, usually leading to blindness. PRA has been reported in dogs from more than 100 breeds and can be genetically heterogeneous both between and within breeds. The disease can be subdivided by age at onset and rate of progression. Using genome‐wide association with 15 Shetland Sheepdog (Sheltie) cases and 14 controls, we identified a novel PRA locus on CFA13 (P_raw = 8.55 × 10^?7, P_genome = 1.7 × 10^?4). CNGA1, which is known to be involved in human cases of retinitis pigmentosa, was located within the associated region and was considered a likely candidate gene. Sequencing of this gene identified a 4‐bp deletion in exon 9 (c.1752_1755delAACT), leading to a frameshift and a premature stop codon. The study indicated genetic heterogeneity as the mutation was present in all PRA‐affected individuals in one large family of Shelties, whereas some other cases in the studied Sheltie population were not associated with this CNGA1 mutation. To our knowledge, this is the first report of a mutation in CNGA1 causing PRA in dogs.     

Genome‐wide association study reveals novel genes for the ear size in sheep (Ovis aries)

Variations in ear size can be observed in livestock such as sheep; however, the genetic basis of variable ear size in sheep is still poorly understood. To investigate causative genes associated with ear size in sheep, a genome‐wide association study was performed in 115 adult Duolang sheep with different‐sized floppy ears using the Ovine Infinium HD BeadChip. We found 38 significant SNPs at the genome‐wide or chromosome‐wise 5% significance level after Bonferroni correction. The most significant association (P = 1.61 × 10^?6) was found at SNP rs402740419, located in the DCC gene, which plays a critical role in ear development. Also, we observed two additional significant SNPs, rs407891215 in PTPRD and rs407769095 in SOX5, both of which are functionally associated with ear developmental processes. Our results are useful for future sheep breeding and provide insights into the genetic basis of ear size development in sheep and other livestock.     

A genome‐wide significant association on chromosome 2 for footrot resistance/susceptibility in Swiss White Alpine sheep

Footrot is one of the most important causes of lameness in global sheep populations and is characterized by a bacterial infection of the interdigital skin. As a multifactorial disease, its clinical representation depends not only on pathogen factors and environmental components but also on the individual resistance/susceptibility of the host. A genetic component has been shown in previous studies; however, so far no causative genetic variant influencing the risk of developing footrot has been identified. In this study, we genotyped 373 Swiss White Alpine sheep, using the ovine high‐density 600k SNP chip, in order to run a DNA‐based comparison of individuals with known clinical footrot status. We performed a case–control genome‐wide association study, which revealed a genome‐wide significant association for SNP rs418747104 on ovine chromosome 2 at 81.2 Mb. The three best associated SNP markers were located at the MPDZ gene, which codes for the multiple PDZ domain crumbs cell polarity complex component protein, also known as multi‐PDZ domain protein 1 (MUPP1). This protein is possibly involved in maintaining the barrier function and integrity of tight junctions. Therefore, we speculate that individuals carrying MPDZ variants may differ in their footrot resistance/susceptibility due to modified horn and interdigital skin integrity. In conclusion, our study reveals that MPDZ might represent a functional candidate gene, and further research is needed to explore its role in footrot affected sheep.     

Molecular cloning of WIF1 and HMGA2 reveals ear‐preferential expression while uncovering a missense mutation associated with porcine ear size in WIF1

Considerable diversity exists in porcine ear size, which is an important morphological feature of pig breeds. Previously, we localized four crucial candidate genes—high mobility group AT‐hook 2 (HMGA2), LEM domain‐containing 3 (LEMD3), methionine sulfoxide reductase B3 (MSRB3) and Wnt inhibitory factor 1 (WIF1)—on Sus Scrofa chromosome 5 affecting porcine ear size, then cloned LEMD3 and MSBR3. In this study, we performed rapid amplification of cDNA ends to obtain full‐length cDNA sequences of 2338‐bp WIF1 and 2998‐bp HMGA2. Using quantitative real‐time PCR, we revealed that WIF1 expression was highest in ear cartilage of 60‐day‐old pigs and that this is therefore a better candidate gene for ear size than HMGA2. We further screened coding sequence variants in both genes and identified only one missense mutation (WIF1:c.1167C>G) in a conserved epidermal growth factor‐like domain from the mammalian WIF1 protein. The protein‐altering mutation was significantly associated with ear size across the Large White × Minzhu hybrid and Beijing Black pig populations. When WIF1:c.1167C>G was included as fixed effect in the model to re‐run a genome‐wide association study in the Large White × Minzhu intercross population the P‐value of the peak SNP on SSC5 from re‐running the genome‐wide association study dropped from 2.45E‐12 to 7.33E‐05. Taken together, the WIF1:c.1167C>G could be an important mutation associated with ear size. Our findings provide helpful information for further studies of the molecular mechanisms controlling porcine ear size.     

Late‐onset progressive retinal atrophy in the Gordon and Irish Setter breeds is associated with a frameshift mutation in C2orf71

Progressive retinal atrophy (PRA) in dogs is characterised by the degeneration of the photoreceptor cells of the retina, resulting in vision loss and eventually complete blindness. The condition affects more than 100 dog breeds and is known to be genetically heterogeneous between breeds. Around 14 mutations have now been identified that are associated with PRA in around 49 breeds, but for the majority of breeds the mutation(s) responsible have yet to be identified. Using genome‐wide association with 16 Gordon Setter PRA cases and 22 controls, we identified a novel PRA locus, termed rod–cone degeneration 4 (rcd4), on CFA17 (P_raw = 2.22 × 10^?8, P_genome = 2.00 × 10^?5), where a 3.2‐Mb region was homozygous within cases. A frameshift mutation was identified in C2orf71, a gene located within this region. This variant was homozygous in 19 of 21 PRA cases and was at a frequency of approximately 0.37 in the Gordon Setter population. Approximately 10% of cases in our study (2 of 21) are not associated with this C2orf71 mutation, indicating that PRA in this breed is genetically heterogeneous and caused by at least two mutations. This variant is also present in a number of Irish Setter dogs with PRA and has an estimated allele frequency of 0.26 in the breed. The function of C2orf71 remains unknown, but it is important for retinal development and function and has previously been associated with autosomal recessive retinitis pigmentosa in humans.