Quantitative trait loci analysis for leg weakness-related traits in a Duroc × Pietrain crossbred population

Unravelling novel and pleiotropic genes for cannon bone circumference and bone mineral density in Yorkshire pigs

Leg weakness is a prevalent health condition in pig farms. The augmentation of cannon bone circumference and bone mineral density can effectively improve limb strength in pigs and alleviate leg weakness. This study measured forelimb cannon bone circumference (fCBC) and rear limb cannon bone circumference (rCBC) using an inelastic tapeline and rear limb metatarsal area bone mineral density (raBMD) using a dual-energy X-ray absorptiometry bone density scanner. The samples of Yorkshire castrated boars were genotyped using a 50K single-nucleotide polymorphism (SNP) array. The SNP-chip data were imputed to the level of whole-genome sequencing data (iWGS). This study used iWGS data to perform genome-wide association studies and identified novel significant SNPs associated with fCBC on SSC6, SSC12, and SSC13, rCBC on SSC12 and SSC14, and raBMD on SSC7. Based on the high phenotypic and genetic correlations between CBC and raBMD, multi-trait meta-analysis was performed to identify pleiotropic SNPs. A significant potential pleiotropic quantitative trait locus (QTL) regulating both CBC and raBMD was identified on SSC15. Bayes fine mapping was used to establish the confidence intervals for these novel QTLs with the most refined confidence interval narrowed down to 56 kb (15.11 to 15.17 Mb on SSC12 for fCBC). Furthermore, the confidence interval for the potential pleiotropic QTL on SSC15 in the meta-analysis was narrowed down to 7.45 kb (137.55 to137.56 Mb on SSC15). Based on the biological functions of genes, the following genes were identified as novel regulatory candidates for different phenotypes: DDX42, MYSM1, FTSJ3, and MECOM for fCBC; SMURF2, and STC1 for rCBC; RGMA for raBMD. Additionally, RAMP1, which was determined to be located 23.68 kb upstream of the confidence interval of the QTL on SSC15 in the meta-analysis, was identified as a potential pleiotropic candidate gene regulating both CBC and raBMD. These findings offered valuable insights for identifying pathogenic genes and elucidating the genetic mechanisms underlying CBC and BMD.

A combined GWAS approach reveals key loci for socially-affected traits in Yorkshire pigs

Socially affected traits in pigs are controlled by direct genetic effects and social genetic effects, which can make elucidation of their genetic architecture challenging. We evaluated the genetic basis of direct genetic effects and social genetic effects by combining single-locus and haplotype-based GWAS on imputed whole-genome sequences. Nineteen SNPs and 25 haplotype loci are identified for direct genetic effects on four traits: average daily feed intake, average daily gain, days to 100 kg and time in feeder per day. Nineteen SNPs and 11 haplotype loci are identified for social genetic effects on average daily feed intake, average daily gain, days to 100 kg and feeding speed. Two significant SNPs from single-locus GWAS (SSC6:18,635,874 and SSC6:18,635,895) are shared by a significant haplotype locus with haplotype alleles 'GGG' for both direct genetic effects and social genetic effects in average daily feed intake. A candidate gene, MT3, which is involved in growth, nervous, and immune processes, is identified. We demonstrate the genetic differences between direct genetic effects and social genetic effects and provide an anchor for investigating the genetic architecture underlying direct genetic effects and social genetic effects on socially affected traits in pigs.

Genome-wide association study of bone mineral density trait among three pig breeds

Leg weakness (LW) issues are a great concern for pig breeding industry. And it also has a serious impact on animal welfare. To dissect the genetic architecture of limb-and-hoof firmness in commercial pigs, a genome-wide association study was conducted on bone mineral density (BMD) in three sow populations, including Duroc, Landrace and Yorkshire. The BMD data were obtained by ultrasound technology from 812 pigs (including Duroc 115, Landrace 243 and Yorkshire 454). In addition, all pigs were genotyped using genome-by-sequencing and a total of 224 162 single-nucleotide polymorphisms (SNPs) were obtained. After quality control, 218 141 SNPs were used for subsequent genome-wide association analysis. Nine significant associations were identified on chromosomes 3, 5, 6, 7, 9, 10, 12 and 18 that passed Bonferroni correction threshold of 0.05/(total SNP numbers). The most significant locus that associated with BMD (P value = 1.92e-14) was detected at approximately 41.7 Mb on SSC6 (SSC stands for Sus scrofa chromosome). CUL7, PTK7, SRF, VEGFA, RHEB, PRKAR1A and TPO that are located near the lead SNP of significant loci were highlighted as functionally plausible candidate genes for sow limb-and-hoof firmness. Moreover, we also applied a new method to measure the BMD data of pigs by ultrasound technology. The results provide an insight into the genetic architecture of LW and can also help to improve animal welfare in pigs.

Effect of Cranial Flexion of Pelvic Limbs on Interlaminar Length of the Lumbosacral Space in Sternally and Laterally Recumbent Juvenile Duroc and Adult Yucatan Pigs

Epidural puncture in swine is technically challenging. Several combinations of limb and body positions have been suggested to increase lumbosacral interlaminar space (LSS) and lumbosacral angle (LSA). This study investigated whether cranial hyperflexion of pelvic limbs increased LSS and LSA in laterally and sternally recumbent juvenile Duroc and adult Yucatan pigs and assessed which position produced the largest LSS. Juvenile Duroc (n = 7) and adult Yucatan (n = 7) pigs were euthanized and randomly placed in 4 positions: sternal with neutral limbs, sternal with cranially hyperflexed limbs, lateral with neutral limbs, and lateral with hyperflexed limbs. LSS and LSA were measured on transverse axial CT images of the spine and compared by using multivariate ANOVA and the Student t test. In both age groups, LSS was greater in lateral flexed (juvenile, 7.0 ± 0.7 mm; adult, 15.9 ± 1.1 mm) and sternal flexed (juvenile, 7.5 ± 1 mm; adult, 17.1 ± 1.1 mm) positions than in lateral neutral (juvenile, 5.4 ± 0.9 mm; adult, 9.6 ± 1.6 mm) position. In addition, in both age groups, LSS and LSA in lateral neutral position were smaller than lateral flexed, sternal neutral, and sternal flexed positions. In adults, LSS was greater in lateral flexed and sternal flexed than in sternal neutral position. Hyperflexion of pelvic limbs increases LSS and LSA in sternally recumbent adult Yucatan pigs and laterally recumbent adult Yucatan and juvenile Duroc swine. Increased LSS from positioning pigs with pelvic limbs flexed in sternal or lateral recumbence may facilitate epidural puncture compared with neutral limb positioning.

Whole genome SNPs discovery in Nero Siciliano pig

Autochthonous pig breeds represent an important genetic reserve to be utilized mainly for the production of typical products. To explore its genetic variability, here we present for the first time whole genome sequencing data and SNPs discovered in a male domestic Nero Siciliano pig compared to the last pig reference genome Sus scrofa11.1.A total of 346.8 million paired reads were generated by sequencing. After quality control, 99.03% of the reads were mapped to the reference genome, and over 11 million variants were detected.Additionally, we evaluated sequence diversity in 21 fitness-related loci selected based on their biological function and/or their proximity to relevant QTLs. We focused on genes that have been related to environmental adaptation and reproductive traits in previous studies regarding local breeds. A total of 6,747 variants were identified resulting in a rate of 1 variant every ~276 bases. Among these variants 1,132 were novel to the dbSNP151 database. This study represents a first step in the genetic characterization of Nero Siciliano pig and also provides a platform for future comparative studies between this and other swine breeds.

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.

Balancing selection at a premature stop mutation in the myostatin gene underlies a recessive leg weakness syndrome in pigs

Balancing selection provides a plausible explanation for the maintenance of deleterious alleles at moderate frequency in livestock, including lethal recessives exhibiting heterozygous advantage in carriers. In the current study, a leg weakness syndrome causing mortality of piglets in a commercial line showed monogenic recessive inheritance, and a region on chromosome 15 associated with the syndrome was identified by homozygosity mapping. Whole genome resequencing of cases and controls identified a mutation causing a premature stop codon within exon 3 of the porcine Myostatin (MSTN) gene, similar to those causing a double-muscling phenotype observed in several mammalian species. The MSTN mutation was in Hardy-Weinberg equilibrium in the population at birth, but significantly distorted amongst animals still in the herd at 110 kg, due to an absence of homozygous mutant genotypes. In heterozygous form, the MSTN mutation was associated with a major increase in muscle depth and decrease in fat depth, suggesting that the deleterious allele was maintained at moderate frequency due to heterozygous advantage (allele frequency, q = 0.22). Knockout of the porcine MSTN by gene editing has previously been linked to problems of low piglet survival and lameness. This MSTN mutation is an example of putative balancing selection in livestock, providing a plausible explanation for the lack of disrupting MSTN mutations in pigs despite many generations of selection for lean growth.

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.

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.

Whole-genome sequence analysis reveals differences in population management and selection of European low-input pig breeds

BACKGROUND

A major concern in conservation genetics is to maintain the genetic diversity of populations. Genetic variation in livestock species is threatened by the progressive marginalisation of local breeds in benefit of high-output pigs worldwide. We used high-density SNP and re-sequencing data to assess genetic diversity of local pig breeds from Europe. In addition, we re-sequenced pigs from commercial breeds to identify potential candidate mutations responsible for phenotypic divergence among these groups of breeds.

RESULTS

Our results point out some local breeds with low genetic diversity, whose genome shows a high proportion of regions of homozygosis (>50%) and that harbour a large number of potentially damaging mutations. We also observed a high correlation between genetic diversity estimates using high-density SNP data and Next Generation Sequencing data (r = 0.96 at individual level). The study of non-synonymous SNPs that were fixed in commercial breeds and also in any local breed, but with different allele, revealed 99 non-synonymous SNPs affecting 65 genes. Candidate mutations that may underlie differences in the adaptation to the environment were exemplified by the genes AZGP1 and TAS2R40. We also observed that highly productive breeds may have lost advantageous genotypes within genes involve in immune response--e.g. IL12RB2 and STAB1-, probably as a result of strong artificial in the intensive production systems in pig.

CONCLUSIONS

The high correlation between genetic diversity computed with the 60K SNP and whole genome re-sequence data indicates that the Porcine 60K SNP Beadchip provides reliable estimates of genomic diversity in European pig populations despite the expected bias. Moreover, this analysis gave insights for strategies to the genetic characterization of local breeds. The comparison between re-sequenced local pigs and re-sequenced commercial pigs made it possible to report candidate mutations to be responsible for phenotypic divergence among those groups of breeds. This study highlights the importance of low input breeds as a valuable genetic reservoir for the pig production industry. However, the high levels of ROHs, inbreeding and potentially damaging mutations emphasize the importance of the genetic characterization of local breeds to preserve their genomic variability.

Emerging genetic basis of osteochondritis dissecans

Genome-wide association studies (GWAS) provide an unbiased approach in the identification of genes that increase the risk for osteochondritis dissecans (OCD). Recent GWAS in humans, horses, and pigs are reviewed and genes identified. The identified genes tended to cluster with respect to function and biologic processes. GWAS in humans are a critical next step in the effort to provide a better understanding of the causes of OCD, which will, in turn, allow preventive strategies for treatment of adolescents and young adults who are at risk for the development of degenerative joint disease due to the effects of OCD.

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.

An f2 pig resource population as a model for genetic studies of obesity and obesity-related diseases in humans: design and genetic parameters

Obesity is a rising worldwide public health problem. Difficulties to precisely measure various obesity traits and the genetic heterogeneity in human have been major impediments to completely disentangle genetic factors causing obesity. The pig is a relevant model for studying human obesity and obesity-related (OOR) traits. Using founder breeds divergent with respect to obesity traits we have created an F2 pig resource population (454 pigs), which has been intensively phenotyped for 36 OOR traits. The main rationale for our study is to characterize the genetic architecture of OOR traits in the F2 pig design, by estimating heritabilities, genetic, and phenotypic correlations using mixed- and multi-trait BLUP animal models. Our analyses revealed high coefficients of variation (15–42%) and moderate to high heritabilities (0.22–0.81) in fatness traits, showing large phenotypic and genetic variation in the F2 population, respectively. This fulfills the purpose of creating a resource population divergent for OOR traits. Strong genetic correlations were found between weight and lean mass at dual-energy x-ray absorptiometry scanning (0.56–0.97). Weight and conformation also showed strong genetic correlations with slaughter traits (e.g., r_g between abdominal circumference and leaf fat at slaughtering: 0.66). Genetic correlations between fat-related traits and the glucose level vary between 0.35 and 0.74 and show a strong correlation between adipose tissue and impaired glucose metabolism. Our power calculations showed a minimum of 80% power for QTL detection for all phenotypes. We revealed genetic correlations at population level, for the first time, for several difficult to measure and novel OOR traits and diseases. The results underpin the potential of the established F2 pig resource population for further genomic, systems genetics, and functional investigations to unravel the genetic background of OOR traits.

In vivo body composition in autochthonous and conventional pig breeding groups by dual-energy X-ray absorptiometry and magnetic resonance imaging under special consideration of Cerdo Ibérico

The improvement of carcass quality is one of the main breeding goals in pig production. To select appropriate breeding animals, it is of major concern to exactly and reliably analyze the body composition in vivo. Therefore, the objective of the study was to examine whether the combination of dual-energy X-ray absorptiometry (DXA) and magnetic resonance imaging (MRI) offers the opportunity to reliably analyze quantitative and qualitative body composition characteristics of different pig breeding groups in vivo. In this study, a total of 77 pigs were studied by DXA and MRI at an average age of 154 days. The pigs originated from different autochthonous or conventional breeds or crossbreeds and were grouped into six breed types: Cerdo Ibérico (Ib); Duroc × Ib (Du_Ib); White Sow Lines (WSL, including German Landrace and German Large White); Hampshire/Pietrain (Pi_Ha, including Hampshire, Pietrain × Hampshire (PiHa) and Pietrain × PiHa); Pietrain/Duroc (Pi_Du, including Pietrain × Duroc (PiDu) and Pietrain × PiDu); crossbred WSL (PiDu_WSL, including Pietrain × WSL and PiDu × WSL). A whole-body scan was performed by DXA with a GE Lunar DPX-IQ in order to measure the amount and percentage of fat tissue (FM; %FM), lean tissue (LM; %LM) and bone mineral, whereas a Siemens Magnetom Open with a large body coil was used for MRI in the thorax region between 13th and 14th vertebrae in order to measure the area of the loin (LA) and the above back fat area (FA) of both body sides. A GLM procedure using SAS 9.2 was used to analyze the data. As expected, the native breed Ib followed by Du_Ib crossbreeds showed the highest %FM (27.2%, 25.0%) combined with the smallest LA (46.2 cm2, 73.6 cm2), whereas Ib had the lowest BW at an average age of 154 days. Pigs with Pi_Ha origin presented the least %FM (12.4%) and largest LA (99.5 cm2). The WSL and PiDu_WSL showed an intermediate body composition. Therefore, it could be concluded that DXA and MRI and especially their combination are very suitable methods to reliably identify differences in body composition and carcass traits among different pig lines in vivo.

Association and expression study of MMP3, TGFβ1 and COL10A1 as candidate genes for leg weakness-related traits in pigs

The present study was aimed to determine the association between metalloproteinase 3 (MMP3), transforming growth factor beta 1 (TGFβ1) and collagen type X alpha I (COL10A1) gene polymorphisms with traits related to leg weakness in pigs. Three hundred Duroc × Pietrain cross breds (DuPi) and 299 pigs of a commercial population (CP) were used for the experiment. DuPi animals were examined for 10 different traits describing leg and feet structure, osteochondrosis (OC) scores and bone density status. Data of OC score at condylus medialis humeri, condylus medialis femoris and distal epiphysis ulna regions of CP were used for association analysis. Significant association (P < 0.05) was found for MMP3 SNP (g.158 C>T) with OC at head of femur and bone mineral density in the DuPi population. Association (P < 0.05) was found between SNP of TGFβ1 (g.180 G>A) with rear leg score and the principle component denoting both OC and feet and leg scores in the DuPi population. No association was found between COL10A1 (g.72 C>T) and leg weakness related traits. The associations of SNPs with OC traits could not be confirmed in the commercial population. Expression analysis of the three candidate genes was performed to compare between healthy and OC. TGFβ1 was found to be highly expressed (P < 0.05) in the OC compared to healthy cartilages, but no significant different expressions were observed for MMP3 and COL10A1 genes. The present finding suggested that TGFβ1 and MMP3 genes variants have an effect on some of the leg weakness related traits.