Effect of VRTN gene polymorphisms on Duroc pig production and carcass traits, and their genetic relationships

Genomic background of biotypes related to growth, carcass and meat quality traits in Duroc pigs based on principal component analysis

As the swine industry continues to explore pork quality traits alongside growth, feed efficiency and carcass leanness traits, it becomes imperative to understand their underlying genetic relationships. Due to this increase in the number of desirable traits, animal breeders must also consider methods to efficiently perform direct genetic changes for each trait and evaluate alternative selection indexes with different sets of phenotypic measurements. Principal component analysis (PCA) and genome-wide association studies (GWAS) can be combined to understand the genetic architecture and biological mechanisms by defining biological types (biotypes) that relate these valuable traits. Therefore, the main objectives of this study were to: (1) estimate genomic-based genetic parameters; (2) define animal biotypes utilizing PCA; and (3) utilize GWAS to link the biotypes to candidate genes and quantitative trait loci (QTL). The phenotypic dataset included 2583 phenotypic records from female Duroc pigs from a terminal sire line. The pedigree file contained 193,764 animals and the genotype file included 21,309 animals with 35,651 single nucleotide polymorphisms (SNPs). Eight principal components (PCs), accounting for a total of 99.7% of the population variation, were defined for three growth, eight conventional carcass, 10 pork quality and 18 novel carcass traits. The eight biotypes defined from the PCs were found to be related to growth rate, maturity, meat quality and body structure, which were then related to candidate genes. Of the 175 candidate genes found, six of them [LDHA (SSC1), PIK3C3 (SSC6), PRKAG3 (SSC15), VRTN (SSC7), DLST (SSC7) and PAPPA (SSC1)] related to four PCs were found to be associated with previously defined QTL, linking the biotypes with biological processes involved with muscle growth, fat deposition, glycogen levels and skeletal development. Further functional analyses helped to make connections between biotypes, relating them through common KEGG pathways and gene ontology (GO) terms. These findings contribute to a better understanding of the genetic relationships between growth, carcass and meat quality traits in Duroc pigs, enabling breeders to better understand the biological mechanisms underlying the phenotypic expression of these traits.

Advancements in Genetic Marker Exploration for Livestock Vertebral Traits with a Focus on China

In livestock breeding, the number of vertebrae has gained significant attention due to its impact on carcass quality and quantity. Variations in vertebral traits have been observed across different animal species and breeds, with a strong correlation to growth and meat production. Furthermore, vertebral traits are classified as quantitative characteristics. Molecular marker techniques, such as marker-assisted selection (MAS), have emerged as efficient tools to identify genetic markers associated with vertebral traits. In the current review, we highlight some key potential genes and their polymorphisms that play pivotal roles in controlling vertebral traits (development, length, and number) in various livestock species, including pigs, donkeys, and sheep. Specific genetic variants within these genes have been linked to vertebral development, number, and length, offering valuable insights into the genetic mechanisms governing vertebral traits. This knowledge has significant implications for selective breeding strategies to enhance structural characteristics and meat quantity and quality in livestock, ultimately improving the efficiency and quality of the animal husbandry industry.

Genome-wide association study reveals new QTL and functional candidate genes for the number of ribs and carcass length in pigs

The number of ribs (NR) and carcass length (CL) are important economic traits in pig breeding programs. Pigs with a higher NR and longer CL produce greater pork yields. In the present study, Suhuai pigs with NR and CL phenotypes were genotyped using the Neogen^® GGP Porcine 80 K SNP array to identify the QTL affecting NR and CL and dissect the candidate genes for the two traits. The SNP-chip data was imputed to the whole-genome sequence (iWGS) to increase the probability of identifying causal variants. Through genome-wide association studies (GWAS) based on both chip and iWGS data, significant SNPs were detected on Sus scrofa chromosome (SSC) 1, SSC4 and SSC7 for NR and on SSC5, SSC16 and SSC17 for CL. Moreover, two SNPs (H3GA0022644 and WU_10.2_7_103460706) on SSC7 detected in chip-based GWAS were significantly associated with both NR and CL. Through Bayes fine mapping, one reported QTL for NR on SSC7 and two reported QTL for CL on SSC17 were verified, and two new QTL (SSC1: 14.05-15.84 Mb and SSC4: 64.83-66.59 Mb) affecting NR and two new QTL (SSC5: 58.31-59.84 Mb and SSC16: 22.98-23.43 Mb) affecting CL were detected. According to the biological functions of genes, MTHFD1L on SSC1 and SULF1 on SSC4 are novel functional candidate genes for NR, and EMP1 on SSC5 and EGFLAM on SSC16 are novel functional candidate genes for CL. Overall, our findings provide a basis for identifying new causal genes and mutations affecting NR and CL.

Genome-Wide Association Study Reveals Additive and Non-Additive Effects on Growth Traits in Duroc Pigs

Growth rate plays a critical role in the pig industry and is related to quantitative traits controlled by many genes. Here, we aimed to identify causative mutations and candidate genes responsible for pig growth traits. In this study, 2360 Duroc pigs were used to detect significant additive, dominance, and epistatic effects associated with growth traits. As a result, a total number of 32 significant SNPs for additive or dominance effects were found to be associated with various factors, including adjusted age at a specified weight (AGE), average daily gain (ADG), backfat thickness (BF), and loin muscle depth (LMD). In addition, the detected additive significant SNPs explained 2.49%, 3.02%, 3.18%, and 1.96% of the deregressed estimated breeding value (DEBV) variance for AGE, ADG, BF, and LMD, respectively, while significant dominance SNPs could explain 2.24%, 13.26%, and 4.08% of AGE, BF, and LMD, respectively. Meanwhile, a total of 805 significant epistatic effects SNPs were associated with one of ADG, AGE, and LMD, from which 11 sub-networks were constructed. In total, 46 potential genes involved in muscle development, fat deposition, and regulation of cell growth were considered as candidates for growth traits, including CD55 and NRIP1 for AGE and ADG, TRIP11 and MIS2 for BF, and VRTN and ZEB2 for LMD, respectively. Generally, in this study, we detected both new and reported variants and potential candidate genes for growth traits of Duroc pigs, which might to be taken into account in future molecular breeding programs to improve the growth performance of pigs.

Genome-wide association study revealed ABCD4 on SSC7 and GREB1L and MIB1 on SSC6 as crucial candidate genes for rib number in Beijing Black pigs

As one of the few animals with variation in the number of rib pairs (RIB), the pig is a good model to study the mechanism of RIB regulation. Quantitative trait loci (QTL) for porcine RIB are present on Sus scrofa chromosome 7 (SSC7). Although several candidate genes exist in this QTL region on SSC7, the causal gene has yet to be verified. Beijing Black pig with 14-17 RIB is a good population for candidate gene mining and 1104 individuals were genotyped using the Illumina Porcine 50K BeadChip. A total of 14 SNPs from 95.49 to 97.78 Mb on SSC7 showed genome-wide significant association with RIB. On SSC7, a locuszoom plot using pairwise linkage disequilibrium displayed the narrowest linkage region encompassing only two genes, ABCD4 and VRTN. In mice, a transcriptome expression profile was obtained using three embryos at E9.5 (the critical period for rib formation). ABCD4 was highly expressed, but no expression of VRTN was detected. On SSC6, there were four genome-wide significant SNPs from 106.42 to 106.92 Mb associated with RIB. GREB1L and MIB1, in this region, were regarded as novel candidate genes. These results revealed a crucial candidate causal gene on SSC7 and novel genes on SSC6 for rib number and provided interesting new insights into its genetic basis.

Expression analysis and single-nucleotide polymorphisms of SYNDIG1L and UNC13C genes associated with thoracic vertebral numbers in sheep (Ovis aries)

The objective of the current study was to analyze expression levels of synapse differentiation inducing 1-like (SYNDIG1L) and unc-13 homolog C (UNC13C) genes in different tissues, while single-nucleotide polymorphisms (SNPs) of two genes were associated with multiple thoracic vertebrae traits in both Small-tailed Han sheep (STH) and Sunite sheep (SNT). The expression levels of SYNDIG1L and UNC13C were analyzed in the brain, cerebellum, heart, liver, spleen, lung, kidney, adrenal gland, uterine horn, longissimus muscle, and abdominal adipose tissues of two sheep breeds with different thoracic vertebral number (TVN) sheep (T13 groups and T14 groups) by real-time quantitative polymerase chain reaction (RT-qPCR). Meanwhile, the polymorphisms of UNC13C gene g.52919279C>T and SYNDIG1L gene g.82573325C>A in T14 and T13 were genotyped by the Sequenom MassARRAY^® SNP assay, and association analysis was performed with the TVN. The results demonstrated that UNC13C gene was extensively expressed in 11 tissues. The expression of UNC13C gene in longissimus muscle of T14 groups of STH was significantly higher than that of T13 groups (P<0.05). SYNDIG1L gene was overexpressed in brain and cerebellum tissues, and the expression level of UNC13C gene in the brain and cerebellum of T13 groups in SNT was significantly higher than that of T14 groups (P<0.01). Association analysis showed that SNPs found in the UNC13C gene had no significant effects on TVN for both two genes. The polymorphism of SYNDIG1L g.82573325C>A was significantly correlated with the TVN in both STH (P<0.05) and SNT (P<0.01). Taken together, the SYNDIG1L gene was related to thoracic vertebral development, and this variation may be potentially used as a molecular marker to select the multiple thoracic vertebrae in sheep.

Estimation of Genetic Parameters for Pork Quality, Novel Carcass, Primal-Cut and Growth Traits in Duroc Pigs

Simple Summary

There is a growing interest in worldwide swine breeding programs to genetically select for pork quality and primal cuts in addition to the traditional growth and carcass leanness traits. Accurate population genetic parameters are needed to estimate correlated responses to selection and incorporate novel traits in the selection objective. Therefore, we estimated heritabilities and genetic correlations for 39 pork quality, growth and carcass traits in Duroc pigs. In general, moderate and favorable genetic correlations were observed between pork quality (e.g., loin color and marbling scores) and carcass traits. Additionally, moderate to low correlations were found among pork quality, growth and carcass traits. Our findings suggest that pig breeders can successfully incorporate pork quality and novel carcass traits in the selection objectives without undesirable impacts upon growth rate and carcass leanness.

Abstract

More recently, swine breeding programs have aimed to include pork quality and novel carcass (e.g., specific primal cuts such as the Boston butt or belly that are not commonly used in selection indexes) and belly traits together with growth, feed efficiency and carcass leanness in the selection indexes of terminal-sire lines, in order to efficiently produce pork with improved quality at a low cost to consumers. In this context, the success of genetic selection for such traits relies on accurate estimates of heritabilities and genetic correlations between traits. The objective of this study was to estimate genetic parameters for 39 traits in Duroc pigs (three growth, eight conventional carcass (commonly measured production traits; e.g., backfat depth), 10 pork quality and 18 novel carcass traits). Phenotypic measurements were collected on 2583 purebred Duroc gilts, and the variance components were estimated using both univariate and bivariate models and REML procedures. Moderate to high heritability estimates were found for most traits, while genetic correlations tended to be low to moderate overall. Moderate to high genetic correlations were found between growth, primal-cuts and novel carcass traits, while low to moderate correlations were found between pork quality and growth and carcass traits. Some genetic antagonisms were observed, but they are of low to moderate magnitude. This indicates that genetic progress can be achieved for all traits when using an adequate selection index.

Polymorphism of VRTN Gene g.20311_20312ins291 Was Associated with the Number of Ribs, Carcass Diagonal Length and Cannon Bone Circumference in Suhuai Pigs

Simple Summary

An increase in the number of ribs (RIB) could improve carcass length (CL) and body size. Cannon bone circumference (CBC) is a pivotal body size trait, and a large CBC could enhance the capacity to bear excessive body weight, vigorous exercise, and resistance to injuries. Several researchers showed that the vertnin (VRTN) gene g.20311_20312ins291 (NC_010449.5 7: g.20311_20312ins291) is an important variant that is related to RIB and CL of Western pigs. However, it is unknown whether this variant could affect the CBC of pigs. Our study showed that this variant was significantly associated with RIB, carcass diagonal length (CDL), and CBC in Suhuai pigs; therefore, it could be used as a potential molecular marker for improving RIB, CDL, and CBC in this breed.

Abstract

The vertnin (VRTN) gene g.20311_20312ins291 was reported as an important variant related to the number of ribs (RIB), and the ins/ins genotype was advantageous for improving RIB of Western pigs. The purpose of this study was to determine whether the VRTN gene g.20311_20312ins291 influences RIB, carcass traits, and body size traits, including cannon bone circumference (CBC) in Chinese Suhuai pigs. We found that the VRTN gene g.20311_20312ins291 was polymorphic in Suhuai fattening pigs and gilts. The polymorphism of g.20311_20312ins291 was significantly associated with RIB and CDL in Suhuai fattening pigs (p < 0.01), whereas this variant had no influence on carcass weight (CWT). There was a tendency of association between this variant and carcass straight length (CSL) in Suhuai fattening pigs (p = 0.06). The polymorphism of g.20311_20312ins291 was also significantly associated with CBC in Suhuai gilts (p = 0.04). Furthermore, CBC was positively genetically correlated with body length (0.22, p < 0.01) and body weight (0.15, p < 0.01). Our results indicated that the VRTN gene g.20311_20312ins291 could be used as a potential marker for improving RIB, CDL, and CBC in Suhuai pigs.

Genome-wide association study reveals novel loci associated with body size and carcass yields in Pekin ducks

Background

Pekin duck products have become popular in Asia over recent decades and account for an increasing market share. However, the genetic mechanisms affecting carcass growth in Pekin ducks remain unknown. This study aimed to identify quantitative trait loci affecting body size and carcass yields in Pekin ducks.

Results

We measured 18 carcass traits in 639 Pekin ducks and performed genotyping using genotyping-by-sequencing (GBS). Loci-based association analysis detected 37 significant loci for the 17 traits. Thirty-seven identified candidate genes were involved in many biological processes. One single nucleotide polymorphism (SNP) (Chr1_140105435 A > T) located in the intron of the ATPase phospholipid transporting 11A gene (ATP11A) attained genome-wide significance associated with five weight traits. Eight SNPs were significantly associated with three body size traits, including the candidate gene plexin domain containing 2 (PLXDC2) associated with breast width and tensin 3 (TNS3) associated with fossil bone length. Only two SNPs were significantly associated with foot weight and four SNPs were significantly associated with heart weight. In the gene-based analysis, three genes (LOC101791418, TUBGCP3 (encoding tubulin gamma complex-associated protein 3), and ATP11A) were associated with four traits (42-day body weight, eviscerated weight, half-eviscerated weight, and leg muscle weight percentage). However, no loci were significantly associated with leg muscle weight in this study.

Conclusions

The novel results of this study improve our understanding of the genetic mechanisms regulating body growth in ducks and thus provide a genetic basis for breeding programs aimed at maximizing the economic potential of Pekin ducks.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5379-1) contains supplementary material, which is available to authorized users.

Effects of increased vertebral number on carcass weight in PIC pigs

Variation of the vertebral number is associated with carcass traits in pigs. However, results from different populations do not match well with others, especially for carcass weight. Therefore, effects of increased vertebral number on carcass weight were investigated by analyzing the relationship between two loci multi-vertebra causal loci (NR6A1 g.748 C > T and VRTN g.20311_20312ins291) and carcass weight in PIC pigs. Results from the association study between vertebral number and carcass weight showed that increased thoracic number had negative effects on carcass weight, but the results were not statistically significant. Further, VRTN Ins/Ins genotype increased more than one thoracic than that of Wt/Wt genotype on average in this PIC population. Meanwhile, there was a significant negative effect of VRTN Ins on carcass weight (P < 0.05). Thus, our results suggested negative effect of increased thoracic number on carcass weight in PIC pigs.

Genome-wide association study in an F2 Duroc x Pietrain resource population for economically important meat quality and carcass traits

Meat quality is essential for consumer acceptance, it ultimately impacts pork production profitability and it is subject to genetic control. The objective of this study was to map genomic regions associated with economically important meat quality and carcass traits. We performed a genome-wide association (GWA) analysis to map regions associated with 38 meat quality and carcass traits recorded for 948 F2 pigs from the Michigan State University Duroc × Pietrain resource population. The F0, F1, and 336 F2 pigs were genotyped with the Illumina Porcine SNP60 BeadChip, while the remaining F2 pigs were genotyped with the GeneSeek Genomic Profiler for Porcine Low Desnisty (LD) chip, and imputed with high accuracy ( = 0.97). Altogether the genomic dataset comprised 1,019 animals and 44,911 SNP. A Gaussian linear mixed model was fitted to estimate the breeding values and the variance components. A linear transformation was performed to estimate the marker effects and variances. Type I error rate was controlled at a False Discovery Rate of 5%. Seven putative QTL found in this study were previously reported in other studies. Two novel QTL associated with tenderness (TEN) were located on SSC3 [135.6:137.5Mb; False Discovery rate (FDR) < 0.03] and SSC5 (67.3:69.1Mb; FDR < 0.02). The QTL region identified on SSC15 includes Protein Kinase AMP-activated ɣ 3-subunit gene (), which has been associated with 24-h pH (pH24), drip loss (DL) and cook yield (CY). Also, novel candidate genes were identified for TEN in the region on SSC5 [A Kinase (PRKA) Anchor Protein 3 (], and for tenth rib backfat thickness (BF10) [Carnitine O-Acetyltransferase ()] on SSC1. The association of gene polymorphisms with pork quality traits has been reported for several pig populations. However, there are no SNP for this gene on the chip used, thus we genotyped the animals for 2 non-synonymous variants ( and ). We then performed a GWA conditioning on the genotype of both SNP and was associated with pH24, DL, protein content (PRO) and CY ( < 0.004) and T30N with Juiciness, TEN, shear force, pH24, PRO, and CY < 0.04). Finally, we performed a GWA conditioning on the genotype of the SNP peak detected in this study, and T30N remained associated only with PRO ( < 0.02). Therefore, in this study we identified 2 novel QTL regions, suggest 2 novel candidate genes, and conclude that other SNP in PRKAG3 or nearby gene(s) explain the observed associations on SSC15 in this population.

Effects of vertebral number variations on carcass traits and genotyping of Vertnin candidate gene in Kazakh sheep

Objective

The vertebral number is associated with body length and carcass traits, which represents an economically important trait in farm animals. The variation of vertebral number has been observed in a few mammalian species. However, the variation of vertebral number and quantitative trait loci in sheep breeds have not been well addressed.

Methods

In our investigation, the information including gender, age, carcass weight, carcass length and the number of thoracic and lumbar vertebrae from 624 China Kazakh sheep was collected. The effect of vertebral number variation on carcass weight and carcass length was estimated by general linear model. Further, the polymorphic sites of Vertnin (VRTN) gene were identified by sequencing, and the association of the genotype and vertebral number variation was analyzed by the one-way analysis of variance model.

Results

The variation of thoracolumbar vertebrae number in Kazakh sheep (18 to 20) was smaller than that in Texel sheep (17 to 21). The individuals with 19 thoracolumbar vertebrae (T13L6) were dominant in Kazakh sheep (79.2%). The association study showed that the numbers of thoracolumbar vertebrae were positively correlated with the carcass length and carcass weight, statistically significant with carcass length. To investigate the association of thoracolumbar vertebrae number with VRTN gene, we genotyped the VRTN gene. A total of 9 polymorphic sites were detected and only a single nucleotide polymorphism (SNP) (rs426367238) was suggested to associate with thoracic vertebral number statistically.

Conclusion

The variation of thoracolumbar vertebrae number positively associated with the carcass length and carcass weight, especially with the carcass length. VRTN gene polymorphism of the SNP (rs426367238) with significant effect on thoracic vertebral number could be as a candidate marker to further evaluate its role in influence of thoracolumbar vertebral number.

SNP- and haplotype-based genome-wide association studies for growth, carcass, and meat quality traits in a Duroc multigenerational population

Background

The aim of the present study was to compare the power of single nucleotide polymorphism (SNP)-based genome-wide association study (GWAS) and haplotype-based GWAS for quantitative trait loci (QTL) detection, and to detect novel candidate genes affecting economically important traits in a purebred Duroc population comprising seven-generation pedigree. First, we performed a simulation analysis using real genotype data of this population to compare the power (based on the null hypothesis) of the two methods. We then performed GWAS using both methods and real phenotype data comprising 52 traits, which included growth, carcass, and meat quality traits.

Results

In total, 836 animals were genotyped using the Illumina PorcineSNP60 BeadChip and 14 customized SNPs from regions of known candidate genes related to the traits of interest. The power of SNP-based GWAS was greater than that of haplotype-based GWAS in a simulation analysis. In real data analysis, a larger number of significant regions was obtained by SNP-based GWAS than by haplotype-based GWAS. For SNP-based GWAS, 23 genome-wide significant SNP regions were detected for 17 traits, and 120 genome-wide suggestive SNP regions were detected for 27 traits. For haplotype-based GWAS, 6 genome-wide significant SNP regions were detected for four traits, and 11 genome-wide suggestive SNP regions were detected for eight traits. All genome-wide significant SNP regions detected by haplotype-based GWAS were located in regions also detected by SNP-based GWAS. Four regions detected by SNP-based GWAS were significantly associated with multiple traits: on Sus scrofa chromosome (SSC) 1 at 304 Mb; and on SSC7 at 35–39 Mb, 41–42 Mb, and 103 Mb. The vertnin gene (VRTN) in particular, was located on SSC7 at 103 Mb and was significantly associated with vertebrae number and carcass lengths. Mapped QTL regions contain some candidate genes involved in skeletal formation (FUBP3; far upstream element binding protein 3) and fat deposition (METTL3; methyltransferase like 3).

Conclusion

Our results show that a multigenerational pig population is useful for detecting QTL, which are typically segregated in a purebred population. In addition, a novel significant region could be detected by SNP-based GWAS as opposed to haplotype-based GWAS.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-016-0368-3) contains supplementary material, which is available to authorized users.

Possible introgression of the VRTN mutation increasing vertebral number, carcass length and teat number from Chinese pigs into European pigs

Vertnin (VRTN) variants have been associated with the number of thoracic vertebrae in European pigs, but the association has not been evidenced in Chinese indigenous pigs. In this study, we first performed a genome-wide association study in Chinese Erhualian pigs using one VRTN candidate causative mutation and the Illumina Porcine 60K SNP Beadchips. The VRTN mutation is significantly associated with thoracic vertebral number in this population. We further show that the VRTN mutation has pleiotropic and desirable effects on teat number and carcass (body) length across four diverse populations, including Erhualian, White Duroc × Erhualian F₂ population, Duroc and Landrace pigs. No association was observed between VRTN genotype and growth and fatness traits in these populations. Therefore, testing for the VRTN mutation in pig breeding schemes would not only increase the number of vertebrae and nipples, but also enlarge body size without undesirable effects on growth and fatness traits, consequently improving pork production. Further, by using whole-genome sequence data, we show that the VRTN mutation was possibly introgressed from Chinese pigs into European pigs. Our results provide another example showing that introgressed Chinese genes greatly contributed to the development and production of modern European pig breeds.