QTL modulating ear size and erectness in pigs

Ear size and erectness are important conformation measurements in pigs. An F(2) population established by crossing European Large White (small, erect ears) with Chinese Meishan (large, flop ears) was used to study the genetic influence of the two ear traits for the first time. A linkage map incorporating 152 markers on 18 autosomal chromosomes was utilised in a genome scan for QTL. Significant QTL were found on SSC1, 5, 7, 9 and 12 for the two traits. The QTL on SSC5 and SSC7 had major effects and were significant at the genome-wide level (P < 0.01). The QTL on SSC1 for ear erectness also had a major effect and was genome-wide significant (P < 0.01). The 95% confidence interval (CI) of the ear size QTL on SSC5 spanned only 4 cM. The QTL on SSC7 for the two ear traits each had a CI of <20 cM, and their positions overlapped with those of the major QTL affecting subcutaneous fat depths on the same chromosome. This study provides insights on the complex genetic influences underlying pig ear traits and will facilitate positional candidate gene analysis to identify causative DNA variants.     

A genome scan for quantitative trait loci affecting three ear traits in a White Duroc × Chinese Erhualian resource population

Chinese Erhualian pigs have larger and floppier ears compared with White Duroc pigs (small, half- or fully-pricked ears). To identify quantitative trait loci (QTL) for ear weight and area as well as erectness, a genome-wide scan with 194 microsatellites was performed in a White Duroc × Chinese Erhualian resource population (>1000 F₂ animals). Twenty-three genome-wide significant QTL and 12 suggestive QTL were identified. All QTL for ear erectness and size detected in two previous studies, bar two on SSC6 and 9, were confirmed here. The 1% genome-wide significant QTL at 70 cM on SSC5 and at 58 cM on SSC7 have profound and pleiotropic effects on the three ear traits, with Erhualian alleles increasing weight and area but decreasing erectness. Notably, the 95% confidence interval of the QTL for weight and area on SSC7 spanned only 3 cM. New QTL reaching 1% genome-wide significance were found on SSC8 (at 37 cM) for all three ear traits, on SSC4 and 16 for weight and area, and on SSCX for area. Unexpectedly, Erhualian alleles at these loci were associated with lighter and smaller or erect ear. Some new suggestive QTL were also found on other chromosome regions. Almost all the QTL for weight and area had essentially additive effects, while the QTL for erectness on SSC2, 5 and 7 showed not only additive effects but also partial dominance effects of Erhualian alleles. The two most significant QTL on SSC7 and SSC5 could be promising targets for fine mapping and identification of the causative mutations.     

All sheeps and sizes: a genetic investigation of mature body size across sheep breeds reveals a polygenic nature

Mature body size is genetically correlated with growth rate, an important economic trait in the sheep industry. Mature body size has been studied extensively in humans as well as cattle and other domestic animal populations but not in sheep. Six‐hundred and sixteen ewes, across 22 breeds, were measured for 28 linear measurements representing various skeletal parts. PCA from these measures generated principal components 1 and 2 which represented 66 and 7% of the phenotypic variation respectively. Two‐hundred and twenty sheep were genotyped on the Illumina Ovine HD beadchip for a GWAS investigating mature body size and linear body measurements. Forty‐six (Bonferroni P < 0.05) SNP associations across 14 chromosomes were identified utilizing principal component 1, representing overall body size, revealing mature body size to have fewer loci of large effect than other domestic species such as dogs and horses. Genome‐wide associations for individual linear measures identified major quantitative trait loci for withers height and ear length. Withers height was associated (Bonferroni P < 0.05) with 12 SNPs across six chromosomes whereas ear length was associated with a single locus on chromosome 3, containing MSRB3. This analysis identified several loci known to be associated with mature body size in other species such as NCAPG, LCORL, and HMGA2. Mature body size is more polygenic in sheep than other domesticated species, making the development of genomic selection for the trait the most efficient option for maintaining or reducing mature body size in sheep.     

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.     

Linked genetic variants on chromosome 10 control ear morphology and body mass among dog breeds

Background

The domestic dog is a rich resource for mapping the genetic components of phenotypic variation due to its unique population history involving strong artificial selection. Genome-wide association studies have revealed a number of chromosomal regions where genetic variation associates with morphological characters that typify dog breeds. A region on chromosome 10 is among those with the highest levels of genetic differentiation between dog breeds and is associated with body mass and ear morphology, a common motif of animal domestication. We characterised variation in this region to uncover haplotype structure and identify candidate functional variants.

Results

We first identified SNPs that strongly associate with body mass and ear type by comparing sequence variation in a 3 Mb region between 19 breeds with a variety of phenotypes. We next genotyped a subset of 123 candidate SNPs in 288 samples from 46 breeds to identify the variants most highly associated with phenotype and infer haplotype structure. A cluster of SNPs that associate strongly with the drop ear phenotype is located within a narrow interval downstream of the gene MSRB3, which is involved in human hearing. These SNPs are in strong genetic linkage with another set of variants that correlate with body mass within the gene HMGA2, which affects human height. In addition we find evidence that this region has been under selection during dog domestication, and identify a cluster of SNPs within MSRB3 that are highly differentiated between dogs and wolves.

Conclusions

We characterise genetically linked variants that potentially influence ear type and body mass in dog breeds, both key traits that have been modified by selective breeding that may also be important for domestication. The finding that variants on long haplotypes have effects on more than one trait suggests that genetic linkage can be an important determinant of the phenotypic response to selection in domestic animals.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1702-2) contains supplementary material, which is available to authorized users.     

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.     

Genome-wide detection of copy number variants in European autochthonous and commercial pig breeds by whole-genome sequencing of DNA pools identified breed-characterising copy number states

In this study, we identified copy number variants (CNVs) in 19 European autochthonous pig breeds and in two commercial breeds (Italian Large White and Italian Duroc) that represent important genetic resources for this species. The genome of 725 pigs was sequenced using a breed-specific DNA pooling approach (30–35 animals per pool) obtaining an average depth per pool of 42×. This approach maximised CNV discovery as well as the related copy number states characterising, on average, the analysed breeds. By mining more than 17.5 billion reads, we identified a total of 9592 CNVs (~683 CNVs per breed) and 3710 CNV regions (CNVRs; 1.15% of the reference pig genome), with an average of 77 CNVRs per breed that were considered as private. A few CNVRs were analysed in more detail, together with other information derived from sequencing data. For example, the CNVR encompassing the KIT gene was associated with coat colour phenotypes in the analysed breeds, confirming the role of the multiple copies in determining breed-specific coat colours. The CNVR covering the MSRB3 gene was associated with ear size in most breeds. The CNVRs affecting the ELOVL6 and ZNF622 genes were private features observed in the Lithuanian Indigenous Wattle and in the Turopolje pig breeds respectively. Overall, the genome variability unravelled here can explain part of the genetic diversity among breeds and might contribute to explain their origin, history and adaptation to a variety of production systems.     

The evolutionary history of the DMRT3 ‘Gait keeper’ haplotype

A previous study revealed a strong association between the DMRT3:Ser301STOP mutation in horses and alternate gaits as well as performance in harness racing. Several follow‐up studies have confirmed a high frequency of the mutation in gaited horse breeds and an effect on gait quality. The aim of this study was to determine when and where the mutation arose, to identify additional potential causal mutations and to determine the coalescence time for contemporary haplotypes carrying the stop mutation. We utilized sequences from 89 horses representing 26 breeds to identify 102 SNPs encompassing the DMRT3 gene that are in strong linkage disequilibrium with the stop mutation. These 102 SNPs were genotyped in an additional 382 horses representing 72 breeds, and we identified 14 unique haplotypes. The results provided conclusive evidence that DMRT3:Ser301STOP is causal, as no other sequence polymorphisms showed an equally strong association to locomotion traits. The low sequence diversity among mutant chromosomes demonstrated that they must have diverged from a common ancestral sequence within the last 10 000 years. Thus, the mutation occurred either just before domestication or more likely some time after domestication and then spread across the world as a result of selection on locomotion traits.     

A genome-wide association study identifies a new variant associated with word reading fluency in Chinese children

Reading disability exhibited defects in different cognitive domains, including word reading fluency, word reading accuracy, phonological awareness, rapid automatized naming and morphological awareness. To identify the genetic basis of Chinese reading disability, we conducted a genome-wide association study (GWAS) of the cognitive traits related to Chinese reading disability in 2284 unrelated Chinese children. Among the traits analyzed in the present GWAS, we detected one genome-wide significant association (p < 5 × 10⁻⁸) on word reading fluency for one SNP on 4p16.2, within EVC genes (rs6446395, p = 7.33 × 10⁻¹⁰). Rs6446395 also showed significant association with Chinese character reading accuracy (p = 2.95 × 10⁻⁴), phonological awareness (p = 7.11 × 10⁻³) and rapid automatized naming (p = 4.71 × 10⁻³), implying multiple effects of this variant. The eQTL data showed that rs6446395 affected EVC expression in the cerebellum. Gene-based analyses identified a gene (PRDM10) to be associated with word reading fluency at the genome-wide level. Our study discovered a new candidate susceptibility variant for reading ability and provided new insights into the genetics of developmental dyslexia in Chinese children.     

A genome-wide significant association on chromosome 15 for congenital entropion in Swiss White Alpine sheep

Entropion is a known congenital disorder in sheep presumed to be heritable but no causative genetic variant has been reported. Affected lambs show a variable inward rolling of the lower eyelids leading to blindness in severe cases. In Switzerland, the Swiss White Alpine (SWA) breed showed a significantly higher prevalence for entropion than other breeds. A GWAS using 150 SWA sheep (90 affected lambs and 60 controls), based on 600k SNP data, revealed a genome-wide significant signal on chromosome 15. The 0.2 Mb associated region contains functional candidate genes, SMTNL1 and CTNND1. Pathogenic variants in human CTNND1 cause blepharocheilodontic syndrome 2, a rare disorder including eyelid anomalies, and SMTNL1 regulates contraction and relaxation of skeletal and smooth muscle. WGS of a single entropion-affected lamb revealed two private missense variants in SMTNL1 and CTNND1. Subsequent genotyping of both variants in 231 phenotyped SWA sheep was performed. The SMTNL1 variant p.(Asp452Asn) affects an evolutionary conserved residue within an important domain and represents a rare allele, which occurred also in controls. The p.(Glu943Lys) variant in CTNND1 represents a common variant unlikely to cause entropion as the mutant allele occurred more frequently in non-affected sheep. Therefore, we propose that these protein-changing variants are unlikely to explain the phenotype. Additionally, WGS of three further disconcordant pairs of full siblings was carried out but revealed no obvious causative variant. Finally, we conclude that entropion represents a more complex disease caused by different non-coding regulatory variants.     

Detection of nitrogen deficiency QTL in juvenile wild barley introgression linesgrowing in a hydroponic system

Background

Identification of genomic regions that have been targets of selection for phenotypic traits is one of the most important and challenging areas of research in animal genetics. However, currently there are relatively few genomic regions identified that have been subject to positive selection. In this study, a genome-wide scan using ~50,000 Single Nucleotide Polymorphisms (SNPs) was performed in an attempt to identify genomic regions associated with fat deposition in fat-tail breeds. This trait and its modification are very important in those countries grazing these breeds.

Results

Two independent experiments using either Iranian or Ovine HapMap genotyping data contrasted thin and fat tail breeds. Population differentiation using F_ST in Iranian thin and fat tail breeds revealed seven genomic regions. Almost all of these regions overlapped with QTLs that had previously been identified as affecting fat and carcass yield traits in beef and dairy cattle. Study of selection sweep signatures using F_ST in thin and fat tail breeds sampled from the Ovine HapMap project confirmed three of these regions located on Chromosomes 5, 7 and X. We found increased homozygosity in these regions in favour of fat tail breeds on chromosome 5 and X and in favour of thin tail breeds on chromosome 7.

Conclusions

In this study, we were able to identify three novel regions associated with fat deposition in thin and fat tail sheep breeds. Two of these were associated with an increase of homozygosity in the fat tail breeds which would be consistent with selection for mutations affecting fat tail size several thousand years after domestication.     

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_PED and F_ROH was around 0.50 and thereby lower than that found in similar studies in cattle. Mean F_ROH 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.     

Investigation of the genetic architecture of a bone carcass weight QTL on BTA6

A previous analysis of an F₂/Backcross Charolais × Holstein cross population identified the presence of a highly significant QTL on chromosome 6 (BTA6) affecting the proportion of bone in the carcass. Two closely linked QTL affected birth weight (BW) and body length at birth (BBL). In this report, the marker density around the QTL on BTA6 was increased, adding four additional microsatellite markers across the chromosome and 46 SNPs within the target QTL confidence interval. Of the SNPs, 26 were in positional candidate genes and the remaining 20 provided an even distribution of markers in the target QTL region. As a bone‐related trait, the sum of the bone weight for all the left fore‐ and hindquarter joints of the carcass was analysed. We also studied the BW and BBL. Analyses of the data substantially reduced the QTL confidence interval. No strong evidence was found that the QTL for the three traits studied are different, and we conclude that the results are consistent with a single pleiotropic QTL influencing the three traits, with the largest effects on the proportion of bone in the carcass. The analyses also suggest that none of the SNPs tested is the sole causative variant of the QTL effects. Specifically, the SNP in the NCAPG gene previously reported as a causal mutation for foetal growth and carcass traits in other cattle populations was excluded as the causal mutation for the QTL reported here. Polymorphisms located in other previously identified candidate genes including SPP1, ABCG2, IBSP, MEPE and PPARGC1A were also excluded. The results suggest that SNP51_BTA‐119876 is the polymorphism in strongest linkage disequilibrium with the causal mutation(s). Further research is required to identify the causal variant(s) associated with this bone‐related QTL.     

Genetic effects of PRNP gene insertion/deletion (indel) on phenotypic traits in sheep

Prion protein (PRNP) gene is well known for affecting mammal transmissible spongiform encephalopathies (TSE), and is also reported to regulate phenotypic traits (e.g. growth traits) in healthy ruminants. To identify the insertion/deletion (indel) variations of the PRNP gene and evaluate their effects on growth traits, 768 healthy individuals from five sheep breeds located in China and Mongolia were identified and analyzed. Herein, four novel indel polymorphisms, namely, Intron-1-insertion-7bp (I1-7bp), Intron-2-insertion-15bp (I2-15bp), Intron-2-insertion-19bp (I2-19bp), and 3′ UTR-insertion-7bp (3′ UTR-7bp), were found in the sheep PRNP gene. In five analyzed breeds, the minor allelic frequencies (MAF) of the above indels were in the range of 0.008 to 0.986 (I1-7bp), 0.113 to 0.336 (I2-15bp), 0.281 to 0.510 (I2-19bp), and 0.040 to 0.238 (3′ UTR-7bp). Additionally, there were 15 haplotypes and the haplotype ‘I_I2-15bp-D_3’UTR-7bp-D_I2-19bp-D_I1-7bp’ had the highest frequency, which varied from 0.464 to 0.629 in five breeds. Moreover, association analysis revealed that all novel indel polymorphisms were significantly associated with 13 different growth traits (P < 0.05). Particularly, the influences of I2-15bp on chest width (P = 0.001) in Small Tail Han sheep (ewe), 3′ UTR-7bp on chest circumference (P = 0.003) in Hu sheep, and I2-19bp on tail length (P = 0.001) in Tong sheep, were highly significant (P < 0.01). These findings may be a further step toward the detection of indel-based typing within and across sheep breeds, and of promising target loci for accelerating the progress of marker-assisted selection in sheep breeding.     

SNP detection and haplotype analysis in partial sequence of MSTN gene in sheep

Polymerase chain reaction (PCR) products of the MSTN gene amplified from sixty sheep of nine Chinese indigenous sheep breeds and one imported sheep breed were sequenced to identify the single-nucleotide polymorphisms (SNPs) in a 378-bp fragment including intron 2 and exon 3 of the MSTN gene. A total of fifteen SNPs (A1937C, T1942G, C1956T, A1972C, A1990G, A2008C, A2011G, C2019T, A2025C, A2027C, T2085G, T2173C, C2198T, C2210T and C2213T) were detected among the sixty sequenced individuals and they were all located in intron 2. Twelve haplotypes were identified from these fifteen SNPs, of which haplotype I (CGTCGCGTCCGCTTT) and VIII (ATCAAAACAATTCCC) were the two major and basic ones with frequencies of 12.25% and 77.80%, respectively. Haplotype VIII was distributed in all sheep breeds and all individuals of the meat or meat-wool type sheep breeds were homozygous with respect to this haplotype. This suggests that haplotype VIII might be related to meat production traits in sheep. Haplotype I was only distributed in the fur, lambskin type and fur-meat type sheep breeds. This suggests that haplotype I may have some relationship with fur traits in sheep.     

Resequencing and signatures of selection scan in two Siberian native sheep breeds point to candidate genetic variants for adaptation and economically important traits

Russian sheep breeds represent an important economic asset by providing meat and wool, whilst being adapted to extreme climates. By resequencing two Russian breeds from Siberia: Tuva (n = 20) and Baikal (n = 20); and comparing them with a European (UK) sheep outgroup (n = 14), 41 million variants were called, and signatures of selection were identified. High‐frequency missense mutations on top of selection peaks were found in genes related to immunity (LOC101109746) in the Baikal breed and wool traits (IDUA), cell differentiation (GLIS1) and fat deposition (AADACL3) in the Tuva breed. In addition, genes found under selection owing to haplotype frequency changes were related to wool traits (DSC2), parasite resistance (CLCA1), insulin receptor pathway (SOCS6) and DNA repair (DDB2) in the Baikal breed, and vision (GPR179) in the Tuva breed. Our results present candidate genes and SNPs for future selection programmes, which are necessary to maintain and increase socioeconomic gain from Siberian breeds.     

Genome-wide association study for calving performance using high-density genotypes in dairy and beef cattle

Background

Calving difficulty and perinatal mortality are prevalent in modern-day cattle production systems. It is well-established that there is a genetic component to both traits, yet little is known about their underlying genomic architecture, particularly in beef breeds. Therefore, we performed a genome-wide association study using high-density genotypes to elucidate the genomic architecture of these traits and to identify regions of the bovine genome associated with them.

Results

Genomic regions associated with calving difficulty (direct and maternal) and perinatal mortality were detected using two statistical approaches: (1) single-SNP (single nucleotide polymorphism) regression and (2) a Bayesian approach. Data included high-density genotypes on 770 Holstein-Friesian, 927 Charolais and 963 Limousin bulls. Several novel or previously identified genomic regions were detected but associations differed by breed. For example, two genomic associations, one each on chromosomes 18 and 2 explained 2.49 % and 3.13 % of the genetic variance in direct calving difficulty in the Holstein-Friesian and Charolais populations, respectively. Imputed Holstein-Friesian sequence data was used to refine the genomic regions responsible for significant associations. Several candidate genes on chromosome 18 were identified and four highly significant missense variants were detected within three of these genes (SIGLEC12, CTU1, and ZNF615). Nevertheless, only CTU1 contained a missense variant with a putative impact on direct calving difficulty based on SIFT (0.06) and Polyphen (0.95) scores. Using imputed sequence data, we refined a genomic region on chromosome 4 associated with maternal calving difficulty in the Holstein-Friesian population and found the strongest association with an intronic variant in the PCLO gene. A meta-analysis was performed across the three breeds for each calving performance trait to identify common variants associated with these traits in the three breeds. Our results suggest that a portion of the genetic variation in calving performance is common to all three breeds.

Conclusion

The genomic architecture of calving performance is complex and mainly influenced by many polymorphisms of small effect. We identified several associations of moderate effect size but the majority were breed-specific, indicating that breed-specific alleles exist for calving performance or that the linkage phase between genotyped allele and causal mutation varies between breeds.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-015-0126-4) contains supplementary material, which is available to authorized users.