Progeny-testing of full-sibs IBD in a SSC2 QTL region highlights epistatic interactions for fatness traits in pigs

Identification and validation of a regulatory mutation upstream of the BMP2 gene associated with carcass length in pigs

Background

Carcass length is very important for body size and meat production for swine, thus understanding the genetic mechanisms that underly this trait is of great significance in genetic improvement programs for pigs. Although many quantitative trait loci (QTL) have been detected in pigs, very few have been fine-mapped to the level of the causal mutations. The aim of this study was to identify potential causal single nucleotide polymorphisms (SNPs) for carcass length by integrating a genome-wide association study (GWAS) and functional assays.

Results

Here, we present a GWAS in a commercial Duroc × (Landrace × Yorkshire) (DLY) population that reveals a prominent association signal (P = 4.49E−07) on pig chromosome 17 for carcass length, which was further validated in two other DLY populations. Within the detected 1 Mb region, the BMP2 gene stood out as the most likely causal candidate because of its functions in bone growth and development. Whole-genome gene expression studies showed that the BMP2 gene was differentially expressed in the cartilage tissues of pigs with extreme carcass length. Then, we genotyped an additional 267 SNPs in 500 selected DLY pigs, followed by further whole-genome SNP imputation, combined with deep genome resequencing data on multiple pig breeds. Reassociation analyses using genotyped and imputed SNP data revealed that the rs320706814 SNP, located approximately 123 kb upstream of the BMP2 gene, was the strongest candidate causal mutation, with a large association with carcass length, with a ~ 4.2 cm difference in length across all three DLY populations (N = 1501; P = 3.66E−29). This SNP segregated in all parental lines of the DLY (Duroc, Large White and Landrace) and was also associated with a significant effect on body length in 299 pure Yorkshire pigs (P = 9.2E−4), which indicates that it has a major value for commercial breeding. Functional assays showed that this SNP is likely located within an enhancer and may affect the binding affinity of transcription factors, thereby regulating BMP2 gene expression.

Conclusions

Taken together, these results suggest that the rs320706814 SNP on pig chromosome 17 is a putative causal mutation for carcass length in the widely used DLY pigs and has great value in breeding for body size in pigs.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-021-00689-0.

Porcine ubiquitin-like 5 (UBL5) gene: genomic organization, polymorphisms, mRNA cloning, splicing variants and association study

Ubiquitin-like 5 (UBL5), which is supposed to be involved in regulation of feed intake, energy metabolism, obesity and type 2 diabetes, is located at position 62.1 cM on the pig chromosome 2 region harbouring quantitative trait loci for carcass and meat quality. The 4,354 bp genomic sequence (FR798948) of the porcine gene encompassing the promoter and entire gene was cloned by polymerase chain reaction. Comparative sequencing revealed 13 polymorphisms in noncoding regions. Synthesis of full-length cDNA sequences using rapid amplification of 5' and 3' ends showed three splice variants. Variants 1 and 2 differ in transcription length for the untranslated part of exon 1 with deduced protein of 73 amino acid (aa) residues and 100 % identities between human, mouse and other species. Variant 3, with 4 bp deletion at the 3' end of exon 2, encodes a truncated protein with 28 aa residues. In a Wild boar×Meishan F2 population (n = 334) with 47 recorded traits, loci FR798948:g.2788G>A and FR798948:g.2141T>C were associated at nominal P < 0.05 with fat deposition, growth and fattening and muscling but after adjustment for multiple testing (Benjamini and Hochberg, J R Stat Soc B 57:289-300, 1995) only eight fat deposition traits showed suggestive association with FR798948:g.2788G>A at adjusted P < 0.10. In a Meishan×Large White (MLW) cross (n = 562) with six trait records available, FR798948:g.2141T>C showed suggestive association with growth (adjusted P = 0.0690). As association mapping conducted in the outbred MLW population is more precise than in the three generation F2 population the UBL5 gene tends to be associated with growth rather than with fat accretion.

Fine mapping of fatness QTL on porcine chromosome X and analyses of three positional candidate genes

BACKGROUND

Porcine chromosome X harbors four QTL strongly affecting backfat thickness (BFT), ham weight (HW), intramuscular fat content (IMF) and loin eye area (LEA). The confidence intervals (CI) of these QTL overlap and span more than 30 cM, or approximately 80 Mb. This study therefore attempts to fine map these QTL by joint analysis of two large-scale F₂ populations (Large White × Meishan and White Duroc × Erhualian constructed by INRA and JXAU respectively) and furthermore, to determine whether these QTL are caused by mutations in three positional candidate genes (ACSL4, SERPINA7 and IRS4) involved in lipid biosynthesis.

RESULTS

A female-specific linkage map with an average distance of 2 cM between markers in the initial QTL interval (SW2456-SW1943) was created and used here. The CI of QTL for BFT, HW and LEA were narrowed down to 6-7 cM, resulting from the joint analysis. For IMF, two linked QTL were revealed in the INRA population but not in the JXAU population, causing a wider CI (13 cM) for IMF QTL. Linkage analyses using two subsets of INRA F₁ dam families demonstrate that the BFT and HW QTL were segregating in the Meishan pigs. Moreover, haplotype comparisons between these dams suggest that within the refined QTL region, the recombination coldspot (~34 Mb) flanked by markers MCSE3F14 and UMNP1218 is unlikely to contain QTL genes. Two SNPs in the ACSL4 gene were identified and showed significant association with BFT and HW, but they and the known polymorphisms in the other two genes are unlikely to be causal mutations.

CONCLUSION

The candidate QTL regions have been greatly reduced and the QTL are most likely located downstream of the recombination coldspot. The segregation of SSCX QTL for BFT and HW within Meishan breed provides an opportunity for us to make effective use of Meishan chromosome X in crossbreeding. Further studies should attempt to identify the impact of additional DNA sequence (e.g. CNV) and expression variation in the three genes or their surrounding genes on these traits.

Association mapping of quantitative trait loci for carcass and meat quality traits at the central part of chromosome 2 in Italian Large White pigs

Association mapping of the central part of porcine chromosome 2 harboring QTLs for carcass and meat quality traits was performed with 17 gene-tagged SNPs located between 44.0 and 77.5 Mb on a physical map (Sscrofa10.2) in Italian Large White pigs. For the analyzed animals records of estimated breeding values for average daily gain, back fat thickness, lean cuts, ham weight, feed conversion ratio, pH1, pHu, CIE L*, CIE a*, CIE b* and drip loss were available. A significant QTL for fat deposition (adjusted P=0.0081) and pH1 (adjusted P=0.0972) to MYOD1 at position 44.4 Mb and a QTL for growth and meatiness (adjusted P=0.0238-0.0601) to UBL5 at position 68.9 Mb were mapped. These results from association mapping are much more accurate than those from linkage mapping and facilitate further search for position candidate genes and causative mutations needed for application of markers through marker assisted selection.