{'fraction_non_alphanumeric': 0.03094691168963757, 'fraction_numerical': 0.021918070444104133, 'mean_word_length': 5.003639149588201, 'pattern_counts': {'":': 0, '<': 4, '<?xml version=': 0, '>': 0, 'https://': 1, 'lorem ipsum': 0, 'www.': 1, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 1, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Genetic architecture and major genes for backfat thickness in pig lines of diverse genetic backgrounds

Background Backfat thickness is an important carcass composition trait for pork production and is commonly included in swine breeding programmes. In this paper, we report the results of a large genome-wide association study for backfat thickness using data from eight lines of diverse genetic backgrounds. Methods Data comprised 275,590 pigs from eight lines with diverse genetic backgrounds (breeds included Large White, Landrace, Pietrain, Hampshire, Duroc, and synthetic lines) genotyped and imputed for 71,324 single-nucleotide polymorphisms (SNPs). For each line, we estimated SNP associations using a univariate linear mixed model that accounted for genomic relationships. SNPs with significant associations were identified using a threshold of p < 10–6 and used to define genomic regions of interest. The proportion of genetic variance explained by a genomic region was estimated using a ridge regression model. Results We found significant associations with backfat thickness for 264 SNPs across 27 genomic regions. Six genomic regions were detected in three or more lines. The average estimate of the SNP-based heritability was 0.48, with estimates by line ranging from 0.30 to 0.58. The genomic regions jointly explained from 3.2 to 19.5% of the additive genetic variance of backfat thickness within a line. Individual genomic regions explained up to 8.0% of the additive genetic variance of backfat thickness within a line. Some of these 27 genomic regions also explained up to 1.6% of the additive genetic variance in lines for which the genomic region was not statistically significant. We identified 64 candidate genes with annotated functions that can be related to fat metabolism, including well-studied genes such as MC4R, IGF2, and LEPR, and more novel candidate genes such as DHCR7, FGF23, MEDAG, DGKI, and PTN. Conclusions Our results confirm the polygenic architecture of backfat thickness and the role of genes involved in energy homeostasis, adipogenesis, fatty acid metabolism, and insulin signalling pathways for fat deposition in pigs. The results also suggest that several less well-understood metabolic pathways contribute to backfat development, such as those of phosphate, calcium, and vitamin D homeostasis. Supplementary Information The online version contains supplementary material available at 10.1186/s12711-021-00671-w.

growth efficiency and lean meat content [2]. This is typically achieved by including backfat thickness in the economic index for selection within pig lines. Backfat thickness is a good indirect predictor of overall body fat content, can be measured on the live animal by ultrasound, and has a high heritability [3][4][5][6].
Over the last two decades, there has been great interest in identifying candidate genes that regulate backfat thickness. Prominent genes that were discovered by linkage analysis include IGF2 [7][8][9], MC4R [10,11], and LEPR [12]. Since then, more than 1400 quantitative trait loci (QTL) related to backfat thickness have been reported (https:// www. anima lgeno me. org/ QTLdb). Results from these studies showed that backfat thickness is a polygenic trait that is regulated by a large number of small-effect variants. With the advent of single-nucleotide polymorphism (SNP) genotyping arrays, gene expression analyses, and other high-throughput genotyping technologies, many more candidate genes for backfat thickness have been reported that are involved in very diverse biological functions and metabolic pathways, such as: adipogenesis [13,14]; lipid metabolism (biosynthesis, absorption, transport, catabolism and homeostasis) pathways, including those related to fatty acids and triglycerides [13,15,16]; regulation of feed intake and energy homeostasis, through hormone-mediated responses [17][18][19][20] or even taste perception [21]; the adipocytokine signalling pathway [17,19]; the vitamin D metabolic pathway [13]; and nervous system development and regulation [22].
The accumulation of evidence for the association of genomic regions with backfat thickness across diverse genetic backgrounds could disentangle which of the reported QTL represent the most prevalent genes and pathways that underlie backfat deposition. In turn, it could also be hypothesized that less prevalent and population-specific associations may derive from variants with larger effects that have been (nearly) fixed in intensely selected populations. In this study, we performed a large genome-wide association study (GWAS) for backfat thickness in eight pig breeding lines of diverse genetic backgrounds, with ~ 15,000 to ~ 55,000 pigs each, for a total of 275,590 pigs. Our main objectives were to determine the genetic architecture of backfat thickness and to identify the main genes and pathways that underlie its genetic variance.

Data
Data comprised 278,112 purebred pigs from eight lines (A to H) of diverse genetic backgrounds (Table 1)  Backfat thickness was measured by ultrasound in the live pigs at about 145 days of age at the tenth rib. Phenotype values were preadjusted for non-genetic effects (contemporary group, litter, and weight) by line. In total, 2522 outlier preadjusted phenotype values, defined as those outside ± 3 standard deviations of the mean within line, were excluded, and 275,590 records remained for further analyses. Pigs were genotyped with either the GGP-Porcine LD BeadChip with 15 k SNPs or the GGP-Porcine HD BeadChip with 50-80 k SNPs (GeneSeek, Lincoln, NE). We used SNPs that mapped to autosomes based on the reference genome version Sscrofa11.1 and excluded SNPs with a call rate lower than 0.95 and a minor allele frequency lower than 0.01. We also excluded individuals with more than 10% missing genotypes. The remaining SNP genotypes were imputed using multi-locus iterative peeling with the AlphaPeel software [23]. Table 1 summarises the number of individuals and SNPs per line that remained after filtering.

Genome-wide association study
For each line, we estimated SNP associations by fitting a univariate linear mixed model that accounted for the genomic relationship matrix as: where y is the vector of preadjusted phenotypes, x i is the vector of genotypes for the i th SNP coded as 0 and 2 if homozygous for either allele or 1 if heterozygous, β i is the additive effect of the i th SNP on the trait, u ∼ N (0, σ 2 u K) is the vector of polygenic effects with the covariance matrix equal to the product of the polygenic additive variance σ 2 u and the genomic relationship matrix K , and e is a vector of uncorrelated residuals. The genomic relationship matrix K was calculated using centred nonstandardized SNP genotypes. We used the GEMMA 0.96 software [24] to fit the model. To assess that the GWAS did not have an increased rate of false positives, we inspected the distribution of the p-values in quantilequantile (Q-Q) plots using the qqman R package [25]. We applied a Bonferroni correction for multiple tests and considered SNP associations with a p-value less than 10 -6 , as significant.
For each line, we defined genomic regions of interest that harboured significant SNPs by applying 0.5-Mb flanking regions downstream and upstream of the significant SNP. The genomic regions that overlapped because they arose from nearby significant SNPs within the same line were merged into a single genomic region and those that overlapped across lines were also merged into a single larger genomic region to facilitate comparison across lines.

SNP-based heritability and genetic variance partitioning by genomic region
To estimate the SNP-based heritability and the genetic variance explained by each genomic region, we fitted a ridge regression model, as implemented in AlphaBayes [26], which uses a Bayesian approach with a Gaussian prior for the SNP effects, a flat prior for the intercept, and a scaled-inverse chi-squared prior for the residual variance. Posterior samples of the SNP effects within each genomic region were obtained from 60,000 Markovchain Monte Carlo iterations after a burn-in period of 20,000 iterations. In each iteration, the total additive genetic variance was calculated as the variance of breeding values across all individuals. The breeding value of the j th individual was calculated as BV j = x j β , where x j is the vector of genotypes of all SNPs of the j th individual, and β is the vector of corresponding SNP effects. For each genomic region, regional breeding values were calculated for all individuals using only the subset of SNPs in each genomic region. The variance of the breeding values obtained for each genomic region was calculated and divided by the total additive genetic variance to estimate the proportion of the additive genetic variance explained by the genomic region. The SNP-based heritability was calculated as the total additive genetic variance divided by the phenotypic variance. All breeding values, variances, and variance ratios were calculated in each iteration to obtain posterior distributions for the proportion of the SNP-based heritability and the genetic variance explained by each genomic region. We summarised these posterior distributions by reporting the median value.

Functional candidate genes and previously reported QTL
The genes located within each genomic region were extracted using the BioMart tool of the Ensembl Genome Browser (Ensembl Genes 100). In order to detect potential functional candidate genes, gene annotation was retrieved from databases of the Gene Ontology project and the KEGG Pathway Database integrated in the Enrichr gene analysis [27].
Data belonging to previously reported QTL that mapped to the reference genome version Sscrofa11.1 were downloaded from the Animal QTLdb [28] (February 2021). QTL entries for traits related to backfat thickness (e.g., average backfat thickness, backfat thickness at last rib, or backfat thickness at last lumbar vertebrae), fat metabolism and deposition (e.g., obesity index, intramuscular fat or triglycerides level), and feed efficiency (e.g., daily feed intake or feed conversion ratio) were selected. Only entries for QTL that were shorter than 5 Mb and that overlapped with the genomic regions found in our GWAS were retained. Enrichment of the genomic regions for QTL terms was tested using the hypergeometric test approach implemented in the GALLO package [29], where the number of QTL entries in the genomic regions identified by the GWAS was compared with the total number of QTL entries for the same term along the whole genome.

Results
We found significant genome-wide associations with backfat thickness for 264 SNPs in 27 genomic regions, of which six were detected in three or more lines. Genomewide associations by line are shown in Fig. 1. The significant SNPs (p < 10 -6 ) and their location in the genome are in Table 2. In general, estimates of SNP effects were low to moderate, but a small fraction of SNPs had larger effects of up to 0.55 mm (0.30 additive genetic standard deviation units) [see Additional file 1: Figure S1]. The estimates of SNP effects were largely consistent across lines. The correlations of estimates of SNP effects between lines were positive (0.05 to 0.18 when all SNPs were considered, Fig. 2), and especially high for SNPs that were significant in at least one line (0.22 to 0.70, Fig. 3). The Q-Q plots for each line are in Figure S2 [see Additional file 1: Figure S2  Within the 27 genomic regions, we identified 1219 positional candidate genes, of which 64 are annotated to functions that can be related to fat metabolism ( Table 4). The functional annotation of these 64 candidate genes supports a role for energy homeostasis genes in regulating backfat thickness development in pigs, such as MC4R on SSC1 at ~ 160.8 Mb (significant association in five lines). Other candidate genes with fat-related functions were identified, such as MEDAG on SSC11 at ~ 7.5 Mb, which is involved in adipocyte differentiation and showed a significant association in four lines. While this region did not show significant associations in the other lines, it explained 1.6% of the additive genetic variance in one other of the lines for which the genomic region was not statistically significant. Another example is the genomic region that contains the IGF2 gene on SSC2 at ~ 1.5 Mb, which is involved in regulation of fat deposition. This region was significantly associated with backfat thickness in only one line but explained 1.1% of the additive genetic variance in at least one other line. The results also suggest that genes involved in phosphate, calcium, and vitamin D homeostasis pathways, such as FGF23 on SSC5 at ~ 66.0 Mb, contribute to backfat thickness development. Table 4 also shows 207 previously published QTL entries for 21 fatness and feed efficiency traits that overlapped 20 of the regions with significant SNPs. The detected genomic regions were enriched for previous QTL entries for average backfat thickness, backfat thickness at last rib, daily feed intake, and lean meat percentage (p < 10 -4 ). For seven of the regions we found no previously reported QTL, and for 12 of the regions we found previously reported QTL for fat metabolism and deposition or for feed efficiency traits but not for backfat thickness traits. These 12 regions showed a significant association in only one of the lines and, in general, explained a low proportion of the genetic variance.

Discussion
To our knowledge, this is the largest-to-date GWAS for backfat thickness in pigs. We report results from eight large populations, which ranged from ~ 15,000 to ~ 55,000 genotyped pigs and differed in breed of origin and selection history. Large sample sizes are required for high power of GWAS and, thus, this dataset provides valuable insight into the genetic architecture of backfat thickness and the main genes and pathways that underlie its genetic variance. We found significant associations for 27 genomic regions, of which one region was detected in five out of the eight lines, three regions in four lines and two regions in three lines. Moreover, some of these genomic regions explained a relatively large proportion of the additive genetic variance of backfat thickness in lines for which the GWAS detected no significant association. The genomic regions that were detected in five or four of the lines contained candidate genes MC4R, MEDAG, FGF23, DGKI, and PTN. Together with the candidate genes that were found in the other genomic regions, the results support the involvement of energy homeostasis, adipogenesis, fatty acid metabolism, and insulin signalling pathways, and suggest the contribution of other metabolic pathways, which are less well understood, to genetic variation for backfat thickness in pigs, such as the phosphate, calcium, and vitamin D homeostasis pathways.
In the light of these findings, we will focus our discussion on: (1) the genetic architecture of backfat thickness, (2) the role of energy homeostasis genes on backfat thickness, (3) the role of adipogenesis, fatty acid metabolism, and insulin signalling genes on backfat thickness, and (4) the role of phosphate, calcium, and vitamin D homeostasis genes on backfat thickness.

Genetic architecture of backfat thickness
The results of the detected SNP associations and genetic variance partitioning confirm the polygenic architecture of backfat thickness, with many loci with small individual effects and only a small fraction of SNPs with larger effects. The effect of the significant SNPs was largely maintained across the eight studied lines. Our moderate-to-high estimates for SNP-based heritability are in the range of previous SNP-based estimates in purebred Duroc lines (0.31 [30] or 0.37 [19]), Landrace (0.47) [31], Large White (0.35) [31], and Pietrain (0.39) [31], and even in Pietrain crosses with Large White x Landrace (0.45), Meishan (0.73) and wild boar (0.42) [32]. Table 3 Estimates of SNP-based genetic parameters and percentage of additive genetic variance of backfat thickness explained by the significant genomic regions in each line a Italic type indicates a significant association of the genomic region found in the GWAS with data from that line We detected several genomic regions that significantly affected backfat thickness and individually explained up to 8.0% of the genetic variance for the trait. Genomic regions that were detected in more than one line or that explained a large proportion of genetic variance generally overlapped with entries for backfat-related QTL with evidence of enrichment and with candidate genes with plausible annotated functions. The fact that a genomic region was detected in a single line or in multiple lines was not always related to the proportion of genetic variance explained by the region in those lines. Moreover, some genomic regions explained a relevant proportion of genetic variance in lines for which that genomic region was not significant. Previous reports in a Duroc population estimated that a single genomic region, in particular the region on SSC6 where the LEPR gene is located, could explain up to 19.8% of the additive genetic variance of backfat thickness [17]. We observed no instances of Table 4 Functional candidate genes and previously reported QTL for backfat-related traits in the significant genomic regions BFT average backfat thickness, ADAR adipocyte area, AFW abdominal fat weight, BF10R backfat thickness at tenth rib, BF34R backfat thickness between the third and fourth rib, BF67R backfat thickness between the sixth and seventh rib, BFFR backfat thickness at first rib, BFLL backfat thickness at last lumbar, BFLR backfat thickness at last rib, BFMB backfat thickness at mid-back, BFMD backfat thickness above muscle dorsi, BFR backfat thickness at rump, BFW backfat weight, DFI daily feed intake, FCR feed conversion ratio; IMF intramuscular fat content, LMP lean meat percentage, OI obesity index, RFI residual feed intake, SSFT shoulder subcutaneous fat thickness, TGL triglyceride level. Within parentheses: number of entries if there were more than one. *Enriched QTL traits (p < 10 -4 ) any single genomic region that explained such a high percentage of the genetic variance. In that same population, the genomic region on SSC1, which includes the MC4R gene, explained 1.1% of the additive genetic variance [17]. These two regions explained up to 0.5% (LEPR region) and 5.3% (MC4R region) of the additive genetic variance in the lines studied here. After accounting for significant genomic regions from the GWAS, the majority of the additive genetic variance remained distributed across non-significantly associated genomic regions; in the genetic variance partitioning analysis, the residual polygenic term due to these non-significant SNPs explained from 80.5 to 96.8% of the additive genetic variance in each line. Thus, most of the genetic variance is explained by minor loci that were not detected in the GWAS. Taken together, these results confirm that backfat thickness has a polygenic architecture, although some major genes that agree with previous studies [13,17,22] contribute large proportions of the genetic variance of the trait in some lines. In that regard, shifting towards an omnigenic model [33] may provide a more suitable conceptualisation of the genetic architecture of backfat thickness.

SNP-based genetic parameters
Previous studies have estimated that dominance effects account for 4 to 15% of the phenotypic variance, and imprinting effects for 1 to 3% [30,31]. Indeed, major genes such as LEPR and IGF2 have been reported to have dominance [19,34] and imprinting [8] effects, respectively, on backfat thickness. While non-additive effects may be of interest for understanding genetic or physiological mechanisms, they have limited applications in breeding practices [35]. Statistical additive variance captures a fraction of the non-additive effects. Thus, our study focused on the additive variance, because it is the fraction of genetic variance that is most commonly targeted in GWAS and most useful for directional selection.

Energy homeostasis genes
The GWAS results support the role of energy homeostasis genes for genetic variation in backfat thickness in pigs. Significant SNPs in the genomic region on SSC1 at ~ 160 Mb were found in five lines. The functional candidate gene MC4R is located in this region. Energy homeostasis in mammals is a feedback system that balances energy intake and expenditure. The melanocortin-4 receptor that is encoded by MC4R has been described as a critical coordinator of mammalian energy homeostasis and body weight [36]. Mutations in MC4R are well known in pigs and are involved in regulating appetite [11,37]. The role of mutations in MC4R in human monogenic obesity was described in the late 1990s [38,39] and has also been described to affect variation in fatness, growth, and feed intake in different pig breeds [10,20,37,[40][41][42].
Another key regulator of feed intake and energy homeostasis is the LEPR gene. A mutation in this gene was described in an Iberian × Landrace cross [43] that reduced leptin signalling and resulted in greater feed intake and therefore greater carcass fat content [44]. Similar findings were found in Duroc pigs [19] and in a Duroc × Landrace/Large White cross [37]. In fact, studies in pigs confirmed that serum leptin concentration is an effective predictor of fat accumulation [4]. We only found one significant SNP near this gene in one of the eight lines, on SSC6 at 147.5 Mb, which is only 0.6 Mb from the LEPR gene at 146. 8-146.9 Mb. Non-significance of this region in the other lines could be the result of intense selection for feed efficiency in the studied lines [34]. A study on signatures of selection in a Duroc line that was selected for increased intramuscular fat content, with a correlated response for backfat thickness, revealed greater extended haplotype homozygosity in this region compared to a control line [45]. The candidate genes JAK1 and LEPROT also map to this region. In fact, the significant SNP at 147.5 Mb is located in an intron of the JAK1 gene. The annotated functions of JAK1 and its effect on backfat thickness are difficult to disentangle from the effects of other genes, such as LEPR or LEPROT.

Adipogenesis, fatty acid metabolism, and insulin signalling genes
The results of the GWAS also support that other pathways are involved in the development of backfat in pigs, such as adipogenesis, fatty acid metabolism, and insulin signalling. Adipogenesis involves differentiation from pre-adipocytes to adipocytes. The region on SSC11 at ~ 8 Mb, which was identified in four lines, contains the adipogenic gene MEDAG, which promotes adipocyte differentiation and lipid accumulation in mature adipocytes [46] and was shown to be upregulated in fat compared to lean pigs [14]. The same region also contains the gene ALOX5AP, which is involved in subcutaneous fat deposition in pigs [14].
The gene PTN, which is in the region on SSC18 at ~ 10 Mb that was identified in four lines, has a role in a signalling pathway that negatively regulates adipogenesis [47]. Recent in vitro studies in mice demonstrated that PTN plays an essential role in the dynamics of adipose lipid turnover and plasticity, as it preserves insulin sensitivity and regulates energy metabolism and thermogenesis [48]. The gene DGKI, which maps to this same region, has been reported to be under positive selection in polar bears and could be related to the development of corporal fat to provide thermal isolation [49].
Other genes with annotated functions related to adipogenesis that were implicated in this GWAS include genes that have been previously associated with body fat content in other species, such as LRP5 (body fat distribution in humans) [50], BRSK2 (abdominal fat in chicken) [51], and DOLK (subcutaneous fat in lambs) [52]. Of these, LRP5 and BRSK2 are in the 0.00-4.82 Mb genomic region on SSC2. Although we detected a significant association for this region in only one line, several GWAS on backfat thickness in pigs have revealed a significant association of this genomic region with backfat thickness, average daily gain, and meat-to-fat ratio in diverse genetic backgrounds, from F 2 populations derived from breeds such as Pietrain, Large White and Landrace [32] to crosses of Iberian pigs with Landrace, Pietrain, and Duroc [13], and many QTL reports support these findings. This region is gene-rich and includes many candidate genes, such as the INS gene, which encodes insulin that regulates blood glucose levels, promotes cell fat storage, and regulates the activity of enzymes that intervene in lipid metabolism [53], and the IGF2 gene, which encodes the insulin-like growth factor 2, and is widely considered as a major candidate gene for muscle mass and fat deposition in pigs [7-9, 22, 32, 54, 55]. However, it has been suggested that other genes in the same region could have an effect on backfat thickness independent of IGF2, such as CTSD, which encodes a protease [55], and genes related to fatty acid metabolism, such as CPT1A, which is involved in the oxidation of long-chain fatty acids, which are the main storage lipids that form backfat. The FADS1, FADS2, and FADS3 genes, which encode fatty acid desaturases, are located near this region (at ~ 9.6-9.7 Mb). Genes related to fatty acid oxidation, such as PLAAT3 (at ~ 8.4 Mb), and fibroblast growth factor genes, such as FGF19 [13,32], have also been considered as potential candidate genes for growth and fat deposition traits in this region.
The PLIN5 gene, which is located in the significant SSC2 region at ~ 74.3 Mb, is involved in the control of intracellular lipid deposition and some results indicate that it may be involved in regulation of the expression of hormone-sensitive lipase [56]. Other genes of the perilipin family have also been associated to differences in backfat thickness [57].
Across the genomic regions we identified candidate genes with functions in bile acid metabolism. Bile, which is predominantly formed by steroid bile acids synthesized from cholesterol, breaks down fat into monoglyceride and fatty acids that can be absorbed by the digestive tract. The SOAT2 gene [58] is involved in biliary cholesterol metabolism and the SLC51B [59], ATP8B1 [60], and AKR1D1 [61] genes are involved in bile formation. This is not the first study that pointed to bile acid metabolism genes as candidates for backfat thickness in pigs through the mechanism of lipid absorption in the intestine [15], although the BAAT gene that was proposed in that previous study was not in any of the significant genomic regions identified in our study.

Phosphate, calcium and vitamin D homeostasis genes
Due to the biological complexity of polygenic traits, some metabolic pathways that have been less explored could nonetheless also have an impact on backfat thickness. For instance, calcium [62,63] and phosphate [64,65] have been linked to adipocyte differentiation and lipid metabolism in human and rats. We found significant SNPs in the genomic region on SSC5 at ~ 66 Mb in four lines. The candidate gene FGF23 is located in this region. This gene is responsible for phosphate homeostasis through a pathway that involves feedback regulation by phosphate, calcium, and vitamin D [66][67][68]. To our knowledge, this is the first GWAS that detects a significant association of the FGF23 genomic region with backfat thickness in pigs, although this gene was previously linked to mineral utilization and homeostasis in Landrace pigs (not significant after correction for multiple testing) in relation with bone development [69]. However, there is some evidence for a potential mechanistic link between FGF23 and adiposity. Lean adipose tissue secretes adiponectin, which causes a significant reduction in the expression of FGF23 in osteocytes [70], while expanded adipose tissue secretes leptin, which increases FGF23 expression in osteocytes [71]. In humans, clinical studies support a potential role of FGF23 signalling in the metabolic status of individuals, including insulin resistance, dyslipidemia, and obesity [72,73]. For instance, it has been reported that the level of FGF23 in blood was higher in obese compared to normal-weight adolescents [74] and that this level was positively correlated with fat mass and triglyceride levels [75].
The genomic region on SSC18 at ~ 10 Mb that includes the PTN gene was detected in the same four lines as the SSC5 genomic region that contains FGF23. As well as having a role in regulation of adipogenesis (as discussed above) [47], the PTN gene is also involved in vitamin D-dependent regulation of calcium and phosphate homeostasis [76]. To our knowledge, there is, however, no evidence for any interaction between PTN and FGF23 [77].
The DHCR7 gene, which is located in the genomic region on SSC2 at ~ 2.4 Mb close to the IGF2 and other genes, has also been associated with backfat thickness in different pig populations [13,32]. The enzyme encoded by DHCR7 catalyses the conversion of 7-dehydrocholesterol to cholesterol, the final step in the production of cholesterol. 7-dehydrocholesterol is also a precursor for vitamin D and, therefore, DHCR7 plays an important role in vitamin D metabolism in humans [78,79]. However, the role of vitamin D in adiposity is unclear. Meta-analyses in humans have shown that, although low vitamin D levels are commonly observed in obese people, which is likely due to sequestration of the fat-soluble vitamin in adipose tissue, vitamin D supplementation did not consistently reduce body weight in clinical trials [80,81]. Although this gene is in a genomic region that contains other candidate genes with more plausible annotated functions, the characterization and validation of DHCR7 and other genes such as FGF23 and PTN with functions related to phosphate, calcium and vitamin D homeostasis could shed new light on genetic variation for backfat thickness.

Conclusions
Our GWAS results obtained on 275,590 pigs from lines with diverse genetic backgrounds confirmed the polygenic architecture of backfat thickness and the importance of genes associated with energy homeostasis, adipogenesis, fatty acid metabolism, and insulin signalling pathways for fat deposition in pigs. The results also suggested that genes involved in phosphate, calcium, and vitamin D homeostasis contribute to backfat development. While the association with backfat thickness of genes such as MC4R, IGF2, and LEPR has been studied during the last decades, the genomic regions detected here also contained more novel candidate genes, such as DHCR7, FGF23, MEDAG, DGKI, and PTN. We quantified that these and other genomic regions could individually contribute up to 8.0% of the genetic variance in the studied lines. The characterization of genes with annotated functions that are not well understood is challenging but can shed new light on the genetic and physiological mechanisms that control adiposity. Further research on these candidate genes is encouraged in order to identify putative causal genomic variants that contribute to the genetic variance in backfat thickness in pigs and to assess their potential application in swine breeding programmes.