KRT8, FAF1 and PTH1R gene polymorphisms are associated with leg weakness traits in pigs

Genetic basis of phenotypic convergence in pig terminal sires using pathway-based selection signature detection methods

The primary purpose of genetic improvement in lean pig breeds is to enhance production performance. Owing to their similar breeding directions, Duroc and Pietrain pigs are ideal models for investigating the phenotypic convergence underlying artificial selection. However, most important economic traits are controlled by a polygenic basis, so traditional strategies for detecting selection signatures may not fully reveal the genetic basis of complex traits. The pathway-based gene network analysis method utilizes each pathway as a unit, overcoming the limitations of traditional strategies for detecting selection signatures by revealing the selection of complex biological processes. Here, we utilized 13 122 398 high-quality SNPs from whole-genome sequencing data of 48 Pietrain pigs, 156 Duroc pigs and 36 European wild boars to detect selective signatures. After calculating FST and iHS scores, we integrated the pathway information and utilized the r/bioconductor graphite and signet packages to construct gene networks, identify subnets and uncover candidate genes underlying selection. Using the traditional strategy, a total of 47 genomic regions exhibiting parallel selection were identified. The enriched genes, including INO80, FZR1, LEPR and FAF1, may be associated with reproduction, fat deposition and skeletal development. Using the pathway-based selection signatures detection method, we identified two significant biological pathways and eight potential candidate genes underlying parallel selection, such as VTN, FN1 and ITGAV. This study presents a novel strategy for investigating the genetic basis of complex traits and elucidating the phenotypic convergence underlying artificial selection, by integrating traditional selection signature methods with pathway-based gene network analysis.

Structural variant landscapes reveal convergent signatures of evolution in sheep and goats

BACKGROUND

Sheep and goats have undergone domestication and improvement to produce similar phenotypes, which have been greatly impacted by structural variants (SVs). Here, we report a high-quality chromosome-level reference genome of Asiatic mouflon, and implement a comprehensive analysis of SVs in 897 genomes of worldwide wild and domestic populations of sheep and goats to reveal genetic signatures underlying convergent evolution.

RESULTS

We characterize the SV landscapes in terms of genetic diversity, chromosomal distribution and their links with genes, QTLs and transposable elements, and examine their impacts on regulatory elements. We identify several novel SVs and annotate corresponding genes (e.g., BMPR1B, BMPR2, RALYL, COL21A1, and LRP1B) associated with important production traits such as fertility, meat and milk production, and wool/hair fineness. We detect signatures of selection involving the parallel evolution of orthologous SV-associated genes during domestication, local environmental adaptation, and improvement. In particular, we find that fecundity traits experienced convergent selection targeting the gene BMPR1B, with the DEL00067921 deletion explaining ~10.4% of the phenotypic variation observed in goats.

CONCLUSIONS

Our results provide new insights into the convergent evolution of SVs and serve as a rich resource for the future improvement of sheep, goats, and related livestock.

Analysis of 206 whole-genome resequencing reveals selection signatures associated with breed-specific traits in Hu sheep

As an invaluable Chinese sheep germplasm resource, Hu sheep are renowned for their high fertility and beautiful wavy lambskins. Their distinctive characteristics have evolved over time through a combination of artificial and natural selection. Identifying selection signatures in Hu sheep can provide a straightforward insight into the mechanism of selection and further uncover the candidate genes associated with breed-specific traits subject to selection. Here, we conducted whole-genome resequencing on 206 Hu sheep individuals, each with an approximate 6-fold depth of coverage. And then we employed three complementary approaches, including composite likelihood ratio, integrated haplotype homozygosity score and the detection of runs of homozygosity, to detect selection signatures. In total, 10 candidate genomic regions displaying selection signatures were simultaneously identified by multiple methods, spanning 88.54 Mb. After annotating, these genomic regions harbored collectively 92 unique genes. Interestingly, 32 candidate genes associated with reproduction were distributed in nine genomic regions detected. Out of them, two stood out as star candidates: BMPR1B and GNRH2, both of which have documented associations with fertility, and a HOXA gene cluster (HOXA1-5, HOXA9, HOXA10, HOXA11 and HOXA13) had also been linked to fertility. Additionally, we identified other genes that are related to hair follicle development (LAMTOR3, EEF1A2), ear size (HOXA1, KCNQ2), fat tail formation (HOXA10, HOXA11), growth and development (FAF1, CCNDBP1, GJB2, GJA3), fat deposition (ACOXL, JAZF1, HOXA3, HOXA4, HOXA5, EBF4), immune (UBR1, FASTKD5) and feed intake (DAPP1, RNF17, NPBWR2). Our results offer novel insights into the genetic mechanisms underlying the selection of breed-specific traits in Hu sheep and provide a reference for sheep genetic improvement programs.

Genome-Wide Association Studies for Body Conformation Traits in Korean Holstein Population

The objective of this study was to identify quantitative trait loci (QTL) and nearby candidate genes that influence body conformation traits. Phenotypic data for 24 body conformation traits were collected from a population of 2329 Korean Holstein cattle, and all animals were genotyped using the 50 K Illumina bovine SNP chip. A total of 24 genome-wide significant SNPs associated with 24 body conformation traits were identified by genome-wide association analysis. The selection of the most promising candidate genes was based on gene ontology (GO) terms and the previously identified functions that influence various body conformation traits as determined in our study. These genes include KCNA1, RYBP, PTH1R, TMIE, and GNAI3 for body traits; ANGPT1 for rump traits; MALRD1, INHBA, and HOXA13 for feet and leg traits; and CDK1, RHOBTB1, and SLC17A1 for udder traits, respectively. These findings contribute to our understanding of the genetic basis of body conformation traits in this population and pave the way for future breeding strategies aimed at enhancing desirable traits in dairy cattle.

Across-breed genetic investigation of canine hip dysplasia, elbow dysplasia, and anterior cruciate ligament rupture using whole-genome sequencing

Here, we report the use of genome-wide association study (GWAS) for the analysis of canine whole-genome sequencing (WGS) repository data using breed phenotypes. Single-nucleotide polymorphisms (SNPs) were called from WGS data from 648 dogs that included 119 breeds from the Dog10K Genomes Project. Next, we assigned breed phenotypes for hip dysplasia (Orthopedic Foundation for Animals (OFA) HD, n = 230 dogs from 27 breeds; hospital HD, n = 279 dogs from 38 breeds), elbow dysplasia (ED, n = 230 dogs from 27 breeds), and anterior cruciate ligament rupture (ACL rupture, n = 279 dogs from 38 breeds), the three most important canine spontaneous complex orthopedic diseases. Substantial morbidity is common with these diseases. Previous within- and between-breed GWAS for HD, ED, and ACL rupture using array SNPs have identified disease-associated loci. Individual disease phenotypes are lacking in repository data. There is a critical knowledge gap regarding the optimal approach to undertake categorical GWAS without individual phenotypes. We considered four GWAS approaches: a classical linear mixed model, a haplotype-based model, a binary case-control model, and a weighted least squares model using SNP average allelic frequency. We found that categorical GWAS was able to validate HD candidate loci. Additionally, we discovered novel candidate loci and genes for all three diseases, including FBX025, IL1A, IL1B, COL27A1, SPRED2 (HD), UGDH, FAF1 (ED), TGIF2 (ED & ACL rupture), and IL22, IL26, CSMD1, LDHA, and TNS1 (ACL rupture). Therefore, categorical GWAS of ancestral dog populations may contribute to the understanding of any disease for which breed epidemiological risk data are available, including diseases for which GWAS has not been performed and candidate loci remain elusive.

Genome-Wide Estimates of Runs of Homozygosity, Heterozygosity, and Genetic Load in Two Chinese Indigenous Goat Breeds

Runs of homozygosity (ROH) and heterozygosity (ROHet) are windows into population demographic history and adaptive evolution. Numerous studies have shown that deleterious mutations are enriched in the ROH of humans, pigs, cattle, and chickens. However, the relationship of deleterious variants to ROH and the pattern of ROHet in goats have been largely understudied. Here, 240 Guangfeng and Ganxi goats from Jiangxi Province, China, were genotyped using the Illumina GoatSNP50 BeadChip and genome-wide ROH, ROHet, and genetic load analyses were performed in the context of 32 global goat breeds. The classes with the highest percentage of ROH and ROHet were 0.5–2 Mb and 0.5–1 Mb, respectively. The results of inbreeding coefficients (based on SNP and ROH) and ROHet measurements showed that Guangfeng goats had higher genetic variability than most Chinese goats, while Ganxi goats had a high degree of inbreeding, even exceeding that of commercial goat breeds. Next, the predicted damaging homozygotes were more enriched in long ROHs, especially in Guangfeng goats. Therefore, we suggest that information on damaging alleles should also be incorporated into the design of breeding and conservation programs. A list of genes related to fecundity, growth, and environmental adaptation were identified in the ROH hotspots of two Jiangxi goats. A sense-related ROH hotspot (chromosome 12: 50.55–50.81 Mb) was shared across global goat breeds and may have undergone selection prior to goat domestication. Furthermore, an identical ROHet hotspot (chromosome 1: 132.21–132.54 Mb) containing two genes associated with embryonic development (STAG1 and PCCB) was detected in domestic goat breeds worldwide. Tajima’s D and BetaScan2 statistics indicated that this region may be caused by long-term balancing selection. These findings not only provide guidance for the design of conservation strategies for Jiangxi goat breeds but also enrich our understanding of the adaptive evolution of goats.

Discovery of selection-driven genetic differences of Duroc, Landrace, and Yorkshire pig breeds by EigenGWAS and Fst analyses

Pigs are one of the earliest domesticated animals and multiple breeds have been developed to meet the various demands of consumers. EigenGWAS is a novel strategy to identify candidate genes that underlying population genetic differences and to infer candidate regions under selection as well. In this study, EigenGWAS and F_st analyses were performed using the public re-sequencing data of three typical commercial pig breeds, Duroc, Landrace and Yorkshire. The intersection of genome-wide significant SNPs detected by EigenGWAS and top-ranked 1% SNPs of F_st results were treated as signals under selection. Using the data of all three breeds, 3062 signals under selection were detected and the nearby genomic regions within 300 kb upstream and downstream covered 6.54% of whole genome. Pairs of breeds were analysed along with the pathway analysis. The gene function enrichment results indicated that many candidate genes located in the genomic regions of the signals under selection were associated with biological processes related to growth, metabolism, reproduction, sensory perception, etc. Among the candidate genes, the FSHB, AHR, PTHLH, KDR and FST genes were reported to be associated with reproductive performance; the KIT, KITLG, MITF, MC1R and EDNRB genes were previously identified to affect coat colour; the RETREG1, TXNIP, BMP5, PPARD and RBP4 genes were reported to be associated with lipid metabolism and growth traits. The identified genetic differences across the three commercial breeds will advance understanding of the artificial selection history of pigs and the signals under selection will suggest potential uses in pig genomic breeding programmes.

Identification of mRNAs Related to Tibial Cartilage Development of Yorkshire Piglets

Cartilage dysplasia is one of the important reasons for the weakness of pig limbs and hooves. Porcine rickets with weak limbs and hooves bring huge economic losses to the pig industry. However, research on the development of pig cartilage is lacking. This study investigated the key genes and molecular mechanisms involved in cartilage development via an RNA-seq technique. Samples of proximal tibia cartilage were collected from three normal piglets with 1 day, 14 days, and 28 days of age, respectively, and then these samples were divided into two comparison groups (1-day vs. 14-day group, 14-day vs. 28-day group). Through the transcriptome analysis, 108 differentially expressed genes (DEGs), such as FORL2, were obtained from 1-day vs. 14-day comparison group, and 3602 DEGs were obtained from 14-day vs. 28-day comparison group, including SOX9, BMP6, and MMP13. The gene ontology (GO) functional and KEGG pathway enrichment revealed that many functions of DEGs were related to bone development. The pathways of DEGs from Day 1 vs. Day 14 were mainly enriched in mineral absorption, but the DEGs of Day 14 vs. Day 28 were enriched in osteoclast differentiation. Then, the expression patterns of six candidate genes were verified via qPCR. In conclusion, candidate genes affecting cartilage development in Yorkshire pigs were obtained by transcriptome analysis, and the clues showed that Day 14 to Day 28 is a more active and extensive period in cartilage developments, which played a key role in revealing the molecular mechanism of pig cartilage development basis, also compensating for vacancies in cartilage research.

Genetic variants of major genes contributing to phosphate and calcium homeostasis and their association with serum parameters in pigs

Calcium and phosphorus are irreplaceable components of life. Tracking the fate of calcium and phosphorus in organisms deserves high attention due to their relevance in bone metabolism and subsequently animal health. Indeed, bone serves as reservoir for calcium and phosphorus, whose formation and resorption follow specific molecular routes including hormones, receptors, and transcription factors. The objective of the study was to analyze the genetic variation of major components driving mineral utilization such as calcitonin receptor, calcium sensing receptor, fibroblast growth factor 23 (FGF23), parathyroid hormone receptor, osteopontin, stanniocalcin 1, RAF-type zinc finger domain containing 1 (TRAFD1), and vitamin D receptor. A German Landrace pig population (n = 360) was used to perform an association analysis between selected single nucleotide polymorphisms (SNP) and relevant serum parameters (calcium, phosphorus, calcium/phosphorus ratio, alkaline phosphatase). Analyzed SNPs in FGF23 (rs710498025) and TRAFD1 (rs345195312) were significantly (p ≤ 0.05) associated with the serum calcium/phosphorus ratio and serum phosphorus levels, respectively. This might represent a modulation of the homeostatic balance between calcium and phosphorus. Furthermore, TRAFD1 is known to be involved in skeletal disorders which emphasize its link to phosphorus utilization and immune system. However, none of the analyzed genetic variants of these major regulators of phosphate and calcium homeostasis showed significant associations after correction for multiple testing (q value > 0.05). Thus, minor contributors as well as unknown and yet to be elucidated regulators of mineral homeostasis need to be characterized towards the implementation of improved phosphorus efficiency in pig breeding programs.

Genome-wide QTL mapping results for regional DXA body composition and bone mineral density traits in pigs

Abstract. In a previous study, genome-wide mapping of quantitative trait loci (QTL) for five body composition traits, three bone mineral traits and live weight was performed using whole-body dual-energy X-ray absorptiometry (DXA) data. Since QTL for bone mineral traits were rare, the current study aimed to clarify whether the mapping results were influenced by the analysed body regions. Thus, the same material (551 pigs) and methods as in the whole-body QTL mapping study were used. However, for evaluation of the DXA scans, we manually defined two body regions: (i) from the last ribs to the pelvis (A) and (ii) including the pelvis and the hind limbs (P). Since live weight was not affected by the regional analysis, it was omitted from the QTL mapping design. Our results show an overall high consistency of mapping results especially for body composition traits. Two thirds of the initial whole-body QTL are significant for both A and P. Possible causes for the still low number of bone mineral QTL and the lower consistency found for these traits are discussed. For body composition traits, the data presented here show high genome-wide Pearson correlations between mapping results that are based on DXA scans with the time-saving whole-body standard setting and mapping results for DXA data that were obtained by time-consuming manual definition of the regions of interest. However, our results also suggest that whole-body or regional DXA scans might generally be less suitable for mapping of bone mineral traits in pigs. An analysis of single reference bones could be more useful.

Polymorphisms in ten candidate genes are associated with conformational and locomotive traits in Spanish Purebred horses

The Spanish Purebred horses, also known as Andalusian horses, compete to the highest standards in international dressage events. Gait and conformation could be used as early selection criteria to detect young horses with promising dressage ability. Although the genetic background of equine size variation has been recently uncovered, the genetic basis of horse conformational and locomotive traits is not known, hampered by the complex genetic architecture underlying quantitative traits and the lack of phenotypic data. The aim of this study was to validate the loci associated with size in 144 Spanish Purebred horses, and to seek novel associations between loci previously associated with the development of osteochondrosis (OC) lesions and 20 conformational and locomotive traits. Ten loci were associated with different conformational and locomotive traits (LCORL/NCAPG, HMGA2, USP31, MECR, COL24A1, MGP, FAM184B, PTH1R, KLF3 and SGK1), and the LCORL/NCAPG association with size in the Spanish Purebred horse was validated. Except for HMGA2, all polymorphisms seem to influence both the prevalence of OC lesions and morphological characters, supporting the link between conformation and OC. Also, the implication of most genes in either immune and inflammatory responses and cellular growth, or ossification processes, reinforces the role that these mechanisms have in the aetiology of OC, as well as their reflection on the general conformation of the individual. These polymorphisms could be used in marker-assisted selection (MAS) programmes to improve desirable conformational traits, but taking into account their possible detrimental effect on OC prevalence.

Non-invasive methods for the determination of body and carcass composition in livestock: dual-energy X-ray absorptiometry, computed tomography, magnetic resonance imaging and ultrasound: invited review

The ability to accurately measure body or carcass composition is important for performance testing, grading and finally selection or payment of meat-producing animals. Advances especially in non-invasive techniques are mainly based on the development of electronic and computer-driven methods in order to provide objective phenotypic data. The preference for a specific technique depends on the target animal species or carcass, combined with technical and practical aspects such as accuracy, reliability, cost, portability, speed, ease of use, safety and for in vivo measurements the need for fixation or sedation. The techniques rely on specific device-driven signals, which interact with tissues in the body or carcass at the atomic or molecular level, resulting in secondary or attenuated signals detected by the instruments and analyzed quantitatively. The electromagnetic signal produced by the instrument may originate from mechanical energy such as sound waves (ultrasound - US), 'photon' radiation (X-ray-computed tomography - CT, dual-energy X-ray absorptiometry - DXA) or radio frequency waves (magnetic resonance imaging - MRI). The signals detected by the corresponding instruments are processed to measure, for example, tissue depths, areas, volumes or distributions of fat, muscle (water, protein) and partly bone or bone mineral. Among the above techniques, CT is the most accurate one followed by MRI and DXA, whereas US can be used for all sizes of farm animal species even under field conditions. CT, MRI and US can provide volume data, whereas only DXA delivers immediate whole-body composition results without (2D) image manipulation. A combination of simple US and more expensive CT, MRI or DXA might be applied for farm animal selection programs in a stepwise approach.

Genome-wide QTL mapping of nine body composition and bone mineral density traits in pigs

Background

Since the pig is one of the most important livestock animals worldwide, mapping loci that are associated with economically important traits and/or traits that influence animal welfare is extremely relevant for efficient future pig breeding. Therefore, the purpose of this study was a genome-wide mapping of quantitative trait loci (QTL) associated with nine body composition and bone mineral traits: absolute (Fat, Lean) and percentage (FatPC, LeanPC) fat and lean mass, live weight (Weight), soft tissue X-ray attenuation coefficient (R), absolute (BMC) and percentage (BMCPC) bone mineral content and bone mineral density (BMD).

Methods

Data on the nine traits investigated were obtained by Dual-energy X-ray absorptiometry for 551 pigs that were between 160 and 200 days old. In addition, all pigs were genotyped using Illumina’s PorcineSNP60 Genotyping BeadChip. Based on these data, a genome-wide combined linkage and linkage disequilibrium analysis was conducted. Thus, we used 44 611 sliding windows that each consisted of 20 adjacent single nucleotide polymorphisms (SNPs). For the middle of each sliding window a variance component analysis was carried out using ASReml. The underlying mixed linear model included random QTL and polygenic effects, with fixed effects of sex, housing, season and age.

Results

Using a Bonferroni-corrected genome-wide significance threshold of P < 0.001, significant peaks were identified for all traits except BMCPC. Overall, we identified 72 QTL on 16 chromosomes, of which 24 were significantly associated with one trait only and the remaining with more than one trait. For example, a QTL on chromosome 2 included the highest peak across the genome for four traits (Fat, FatPC, LeanPC and R). The nearby gene, ZNF608, is known to be associated with body mass index in humans and involved in starvation in Drosophila, which makes it an extremely good candidate gene for this QTL.

Conclusions

Our QTL mapping approach identified 72 QTL, some of which confirmed results of previous studies in pigs. However, we also detected significant associations that have not been published before and were able to identify a number of new and promising candidate genes, such as ZNF608.

Electronic supplementary material

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

Gene expression profiling of articular cartilage reveals functional pathways and networks of candidate genes for osteochondrosis in pigs

Osteochondrosis (OC) is a joint disorder that frequently causes leg weakness in growing pigs, resulting in welfare problems and economic losses. We aimed to detect molecular pathways relevant to the emergence of the disease and to identify candidate genes for the liability to the disorder. Therefore, we compared microarray-based expression patterns of articular cartilage with (n=11) and without (n=11) histologically diagnosed OC lesions obtained from discordant sib-pairs. A total of 1,564 genes were found with different transcript abundance [differentially expressed (DE) genes] at q≤0.05. To further identify candidate genes, we integrated data from quantitative trait loci (QTL) and genome-wide association (GWA) studies with the expression analysis. We detected 317 DE genes within the QTL confidence intervals, of which 26 DE genes also overlapped GWA regions. Ingenuity Pathway Analysis suggests a pathogenic role of immune response, angiogenesis, and synthesis of extracellular matrix pathways for OC. These processes could facilitate the emergence of defects. But they may also promote the degradation of articular cartilage and the worsening of the disease. A functional network was derived that comprised genes with functional and positional clues of their role in bone and cartilage metabolisms and development, including extracellular matrix genes (e.g., LOX, OGN, and ASPN), angiogenesis genes (e.g., ANGPTL4 and PDGFA), and immune response genes (e.g., ICAM1, AZGP1, C1QB, C1QC, PDE4B, and CDA). The study identified molecular processes linked to OC and several genes with positional, genetic-statistical, and functional evidence for their role in the emergence of articular cartilage lesions and the liability to OC.