Skip to main content
Download PDF

Identification and characterization of structural variants related to meat quality in pigs using chromosome-level genome assemblies

BMC Genomics volume 25, Article number: 299 (2024) Cite this article

Abstract

Background

Many studies have been performed to identify various genomic loci and genes associated with the meat quality in pigs. However, the full genetic architecture of the trait still remains unclear in part because of the lack of accurate identification of related structural variations (SVs) which resulted from the shortage of target breeds, the limitations of sequencing data, and the incompleteness of genome assemblies. The recent generation of a new pig breed with superior meat quality, called Nanchukmacdon, and its chromosome-level genome assembly (the NCMD assembly) has provided new opportunities.

Results

By applying assembly-based SV calling approaches to various genome assemblies of pigs including Nanchukmacdon, the impact of SVs on meat quality was investigated. Especially, by checking the commonality of SVs with other pig breeds, a total of 13,819 Nanchukmacdon-specific SVs (NSVs) were identified, which have a potential effect on the unique meat quality of Nanchukmacdon. The regulatory potentials of NSVs for the expression of nearby genes were further examined using transcriptome- and epigenome-based analyses in different tissues.

Conclusions

Whole-genome comparisons based on chromosome-level genome assemblies have led to the discovery of SVs affecting meat quality in pigs, and their regulatory potentials were analyzed. The identified NSVs will provide new insights regarding genetic architectures underlying the meat quality in pigs. Finally, this study confirms the utility of chromosome-level genome assemblies and multi-omics analysis to enhance the understanding of unique phenotypes.

Peer Review reports

Background

Meat quality has been an important feature in the livestock industry throughout human history. The importance of the trait has increased in the last 50 years with the increased demand for high-quality meats, and the trait is currently a considerably important characteristic in a related industry [1, 2]. Therefore, the livestock industry has continuously developed new breeds that can provide higher-quality meat to meet the demands of customers [3]. As the genomic sequences and annotations of many pig breeds have been revealed [4], it is important to understand the genetic architectures regulating the trait to make a new breed with better meat quality. Based on the maturity of sequencing technologies, diverse genomic studies have been performed to identify genomic loci affecting meat quality in pigs. Especially, several population-level studies using short read data, such as population genomic analyses [5,6,7] and genome-wide association studies [8,9,10], have identified plentiful regions and genes associated with the traits at the whole-genome level.

However, it is still not enough to fully uncover the genetic architectures for the trait because the previous studies mainly focused on small variants such as single nucleotide polymorphisms, small insertions and deletions. Structural variations (SVs), which are large variants involving more than 50 base pairs, also contribute to phenotypic alterations [11,12,13]. However, understanding the effect of SVs on meat quality in pigs has lagged because it has been hard to accurately discover SVs due to the limited availability of sequencing data and the incompleteness of genome assemblies especially for a pig breed specifically inbred for higher meat quality [14, 15]. In this situation, the high-quality and chromosome-level genome assemblies of such pig breeds can be used as valuable resources to accurately identify SVs [16, 17] even in repetitive and complex regions [18].

We previously reconstructed a chromosome-level genome assembly [19], called the NCMD assembly, for Nanchukmacdon which is a new pig breed derived by mating three different commercial pig breeds (Korean native pig, Duroc and Landrace) with outstanding levels of intramuscular fat deposition and redness in meat compared to other pig breeds [19, 20]. Therefore, the high-quality and chromosome-level NCMD assembly can provide new opportunities for fully understanding the impact of SVs on meat quality in pigs.

In this study, we discovered SVs in diverse pig breeds including Nanchukmacdon by comparing the prebuilt NCMD assembly against different pig genome assemblies. Furthermore, we identified several candidate Nanchukmacdon-specific SVs (NSVs) responsible for the superior meat quality of Nanchukmacdon. In addition, the analysis using multi-omics data enabled us to interrogate the regulatory roles of the NSVs for the expression of adjacent genes related to meat quality. This study will establish a foundation for understanding the role of SVs on breed-specific traits, and act as an important resource for future development of breeds with superior meat quality.

Methods

Collinearity comparison of the chromosome-level assemblies of pig breeds

For comparing genome assemblies of pig breeds, the genome assemblies of the Nanchukmacdon (NCMD), Duroc (Sscrofa11.1) and Meishan (MSCAAS v1) were obtained from the NCBI database (accession GCA_031306245.1, GCF_000003025.6, and GCA_017957985.1). Whole-genome sequence alignments of the NCMD assembly against the MSCAAS v1 and Sscrofa11.1 (the Duroc pig) assembly were next constructed using LASTZ (v.1.04.00) [21]. Synteny blocks among them were then constructed by the synteny block detection program in InferCars [22] with a resolution of 300 Kbp. Breakpoint regions, which are genomic regions between two adjacent genomic blocks belonging to different synteny blocks, were identified in each assembly using the in-house script. The syntenic relationships and breakpoint regions between the two assemblies were visualized using mySyntenyPortal [23] and jcvi (https://github.com/tanghaibao/jcvi).

For analyzing the mapping patterns of the sequencing reads at the breakpoint regions, we obtained long and short reads of Nanchukmacdon and Meishan from the previous studies [19, 24]. The short and long reads were first mapped against the corresponding assembly (the NCMD assembly for the reads of Nanchukmacdon, and the MSCAAS v1 assembly for the reads of Meishan). Then, the reads mapped to the breakpoint regions of each assembly were next mapped to the other two assemblies of different pig breeds. The short and long reads were mapped using BWA MEM [25] and minimap2 (v.2.17) [26], respectively. Also. the Hi-C reads of Nanchukmacdon were also obtained from a recent study [19] and aligned to the other two assemblies, and the Hi-C contact maps at the breakpoint regions in each assembly were constructed using the HiC-Pro pipeline (v.2.11.4) [27] with default parameters and 100 Kbp bins. The Hi-C contact maps were visualized using HiCPlotter (v.0.8.1) [28].

To assess the functional effect of genome rearrangements, the breakpoint regions were annotated using SnpEff (v.4.3t) [29] with “-ud 5000” parameter. For the annotation of breakpoint regions, we used the pig reference (Sscrofa11.1.101) and Nanchukmacdon gene annotation information which were downloaded from the NCBI database and a previous paper [19], respectively. We defined the breakpoint regions affecting exons or 3/5 prime UTRs as exonic breakpoint regions.

Variant analysis

To call variants based on genome assemblies, whole-genome pairwise alignment was performed for each of the NCMD and Landrace_pig_v1 assembly (Landrace pig) against the Sscrofa11.1 genome assembly (Duroc pig) using minimap2 (v.2.17) with "-a -x asm5" parameters. For that, the public genome assembly of Landrace was additionally obtained from the NCBI database (GCA_001700215.1). Using the alignment results, structural variants (SVs) were called in each assembly with SVIM-asm (v.1.0.3) [30]. We excluded the breakend variants in this study. To identify candidate Nanchukmacdon-specific variants (NSVs), we filtered out the variants found in the NCMD assembly whose positions on the Sscrofa11.1 assembly overlapped those found in the Landrace_pig_v1 assembly. The called NSVs were annotated using VEP [31] with the pig reference annotation information (Sscrofa11.1.101) downloaded from NCBI database. We defined the stop gained/lost, frameshift variant, inframe insertion/deletion, coding sequence variant, and 3/5 prime UTR variant as exonic variants, splice donor/acceptor/region variant and intron variant as intronic variants, and upstream/downstream gene variant and intergenic variant as intergenic variants. GO enrichment analysis for the NSVs was conducted using g:Profiler [32] with g:SCS threshold 0.05.

NSV genotyping analysis

To genotype NSVs, we obtained the whole genome sequencing data of 10, 6, and 13 samples of Nanchukmacdon, Duroc, and Landrace, respectively, from the NCBI database (Additional file 2: Table S1). The quality of the collected whole genome sequencing data was examined using FastQC (v.0.11.9) [33]. Low-quality reads and adaptor sequences in the reads were removed using NGStoolkit IlluQC.pl (v.2.3.3) [34]. When there were no adapter sequences in the reads, “N A -s 20 -l 70” parameters were used for NGStoolkit to filter only low-quality reads with low-quality base percentages (Phred quality score < 20) larger than 30%. If adapter sequences were present, “2 A -s 20 -l 70” parameters were used. Low-quality bases (Phred quality score < 20) were next trimmed at the 3' end of the reads, and the trimmed reads shorter than 45 bases were removed using NGStoolkit TrimmingReads.pl with “-q 20 -n 45” parameters. If adapter sequences still remained when checked with FastQC, additional trimming was performed using TrimGalore (v.0.6.0) [35] with the following parameters: “-q 20 –length 45 –paired –illumina”. The cleaned reads for each sample were mapped against the pig reference assembly (Sscrofa11.1) using BWA MEM (v.0.7.17) [25]. The NSVs are genotyped using the mapping results with graphtyper (v.2.7.5) [36]. We only used the aggregated genotype calls and filtered out the calls not annotated as PASS for further analysis.

PCA analysis

Based on the genotyped NSVs in autosomes, a principal component analysis (PCA) was performed using the Genome-wide Complex Trait Analysis (GCTA) tool (v.1.91.4) [37]. The genetic relationship matrix for pairs of individuals was estimated using “–make-grm –autosome –autosome-num 18” parameters. Based on the matrix, eigenvalues and eigenvectors were calculated using the “–pca 3” parameter in GCTA. The results were visualized using the ggplot2 R package [38].

Integrative multi-omics analysis for NSV flanking genes

For the differential gene expression (DGE) analysis, we first collected the raw RNA sequencing (RNA-seq) data from 15, 12, and 17 samples of three different pig breeds for adipose, liver and muscle tissues, respectively. The raw RNA-seq data of the three tissues of Nanchukmacdon were obtained from previous studies [20, 39]. The raw RNA-seq data of the other two pig breeds (Duroc and Landrace) were obtained from the NCBI database (Additional file 2: Table S2). For each tissue, a separate DGE analysis was performed for Nanchukmacdon against each of Duroc and Landrace using the RNA-seq data. The DGE analyses were performed by following the procedure in a previous study [20]. In brief, the quality of RNA-seq reads was checked using FastQC [40], and they were trimmed by Trimmomatic [41]. The trimmed RNA-seq reads of each pig breed were aligned to the pig reference assembly (Sscrofa11.1) using HISAT2 [42] with default parameters. For each gene in the pig reference assembly, the number of mapped RNA-seq reads was counted using FeatureCount [43], and the DGE analysis between two pig breeds was performed using the DESeq2 R package [44]. Genes with |log2FC|> = 1 and an adjusted P-value < 0.05 were identified as the final differentially expressed genes (DEGs). Finally, we identified the common DEGs of Nanchukmacdon against the two breeds for each tissue.

For the epigenome analysis, we obtained the H3K4me1, H3K27ac and H3K4me3 histone modification profile and chromatin state information of pigs from a previous study [45]. For visualization, we used histone modification profile data of the P348 sample.

Results

Collinearity comparison with chromosome-level assemblies of pig breeds

We first compared the NCMD assembly against chromosome-level assemblies of two different pig breeds, such as MSCAAS v1 for Meishan and Sscrofa11.1 for Duroc by constructing synteny. The compared assemblies were overall collinear with only four and eight breakpoint regions in the NCMD assembly against the Sscrofa11.1 and MSCAAS v1 assembly, respectively (Fig. 1a).

Fig. 1
figure 1

Comparison of the NCMD assembly with assemblies of other pig breeds. a Syntenic relationships of Nanchukmacdon chromosomes containing breakpoint regions against the Duroc (Sscrofa11.1) and Meishan (MSCAAS v1) chromosomes. The syntenic relationships were visualized using mySyntenyPortal (http://bioinfo.konkuk.ac.kr/mySyntenyPortal_RDA/publish/Nanchukmacdon_2Breed). b, c Examples of breakpoint regions in the Nanchukmacdon genome against Meishan (b) and Duroc (c) genome with Hi-C heatmaps in the regions. Colored regions in each assembly represent the regions corresponding to the breakpoint regions in the assembly. Green triangles pinpoint irregular Hi-C contact frequencies near breakpoint regions. d Paralogous genes flanking in the breakpoint regions between the Nanchukmacdon and Duroc genome. Green and red boxes respectively represent the genes in forward and reverse strand in each genome, respectively

Full size image

To check whether the identified breakpoint regions are the results of misassembly in the NCMD assembly or not, the Hi-C contact maps were created by mapping the Hi-C reads of Nanchukmacdon to each of the assemblies, and the generated Hi-C contact maps were compared. As shown in Fig. 1b-c, normal contact frequencies were observed in all breakpoint regions in the NCMD assembly, while abnormal patterns were found in the corresponding breakpoint regions in the other two assemblies. Mapping patterns of the short and long Nanchukmacdon reads against the three assemblies were also investigated (Additional file 1: Figure S1-S5). Both short and long reads were properly mapped in the breakpoint regions of the NCMD assembly only, but not in those of the other two assemblies. Consistent patterns were also observed when the reads of Meishan were mapped (Additional file 1: Figure S6). These results indicate that the observed genome rearrangements are not artifacts generated by the misassembly of the NCMD assembly.

Six and 35 genes were located within 5 Kbp upstream and downstream of the breakpoint regions against Duroc and Meishan, respectively (Additional file 2: Table S3-S4). Most of the breakpoints were in intronic and intergenic regions, and only one (1.45%) and eight (5.00%) exonic regions were affected by the breakpoint regions against Duroc and Meishan, respectively (Additional file 2: Table S5). Interestingly, some of these genes were paralogous to each other. For example, the igLON Family Member 5 (igLON5; LOC100512905) and igLON5-like (LOC110261043) gene were located close to the boundary of breakpoint regions against Duroc (Fig. 1d). The breakpoint region against Meishan also contained several olfactory receptor genes (Additional file 1: Figure S7).

Identification of structural variants related to Nanchukmacdon-specific phenotypes

We also used the NCMD assembly to identify genomic features related to the superior phenotype of Nanchukmacdon compared to its related breeds (Duroc, Landrace) including higher levels of intramuscular fat deposition and reddish meat color [20]. Here, putative structural variants (SVs) of Nanchukmacdon and Landrace were identified against Duroc based on whole-genome pairwise comparisons (Methods). A total of 30,784 and 62,373 of SVs in length \(\ge\) 50 bp were detected in the genome of Nanchukmacdon and Landrace, respectively (Table 1). Most of the SVs were insertions or deletions in both breeds. As described in the previous genomic variant studies [46, 47], the SVs were not randomly distributed in both breeds, and more SVs were observed in the telomere regions in both breeds (Fig. 2a).

Table 1 Statistics of structural variants in Nanchukmacdon and Landrace genome. NSV: Nanchukmacdon-specific variant, DEL: Deletion, INS: Insertion, DUP: Duplication, INV: Inversion
Full size table
Fig. 2
figure 2

Structural variants (SVs) in Nanchukmacdon and Landrace. a Distribution of SVs in the Nanchukmacdon and Landrace genomes. NSVs represent Nanchukmacdon-specific SVs. b Annotation result for NSVs. c PCA result of 29 samples of three different pig breeds based on NSVs. d Enriched biological process gene ontology terms for differentially expressed genes related to NSVs in the backfat, liver, and muscle tissue

Full size image

To determine the SVs affecting the superior phenotypes of Nanchukmacdon, we found the candidate Nanchukmacdon-specific SVs (NSVs), which were observed only in the Nanchukmacdon genome (Method). We identified 16,898 NSVs including 7,012 deletions, 70 duplications, 9,785 insertions, and 31 inversions (Fig. 2a, Table 1). The NSVs were genotyped using whole genome sequencing data of total 29 samples of Nanchukmacdon, Duroc, Landrace (Additional file 2: Table S6). An average of 7,051 (41.73%) NSVs were genotyped in each sample, and the number of genotyped NSVs of Nanchukmacdon was larger than those of the other breeds. In PCA performed using the genotyped NSVs in autosomes, clear separation between Nanchukmacdon and the other breeds was identified by the first PC (Fig. 2c). To further identify the separation came from the differences of the rate of genotyped NSVs among breeds, we performed the PCA analysis using the NSVs with minor allele frequency not lower than 0.05 and missing rate not larger than 0.1. The clear division between Nanchukmacdon and the other breeds was consistently identified in the PCA using the filtered NSVs (Additional file 1: Figure S8).

When annotating those NSVs, a total of 4,793 genes contained the NSVs within 5 Kbp upstream and downstream of their gene body. We also found that only 2.09% (734) NSVs were mapped within exonic regions. The remaining 64.27% (22,613) and 32.91% (11,581) were in intronic and intergenic regions, respectively (Fig. 2b, Additional file 2: Table S7). For the genes adjacent to or containing the NSVs, differential gene expression analyses were followed. In this analysis, 414, 473 and 558 genes were differentially expressed only in Nanchukmacdon (not both in Duroc and Landrace) in the backfat, liver and muscle tissue of Nanchukmacdon, respectively (Additional file 2: Table S8).

For those differentially expressed genes (DEGs) in each tissue, the gene ontology (GO) enrichment analysis was performed (Additional file 2: Table S9). For biological process terms, the muscle and backfat tissue had the highest and lowest number of enriched terms, respectively (Fig. 2d). In the backfat tissue, the DEGs were mainly involved in several development-related terms, such as developmental process (GO:0032502), anatomical structure development (GO:0048856) and system development (GO:0048731). Interestingly, most of the development-related terms were commonly identified in all three tissues. In addition, we also identified many tissue-specific enriched GO terms in liver and muscle, such as cell adhesion (Cell adhesion (GO:0007155), cell–cell adhesion (GO:0098609)), actin filament (Actin filament-based process (GO:0030029), actin cytoskeleton organization (GO:0030036), actin filament network formation (GO:0051639)), and phosphorylation process (Phosphorylation (GO:0016310), phosphorus metabolic process (GO:0006793)) for liver, and macromolecule modification (Macromolecule modification (GO:0043412), protein modification process (GO:0036211)), cell development (Cell differentiation (GO:0030154), regulation of cell size (GO:0008361), regulation of anatomical structure size (GO:0090066), regulation of cellular component size (GO:0032535)) and blood circulatory system development (Circulatory system development (GO:0072359), blood circulation (GO:0008015), and blood vessel development (GO:0001568)) for muscle. These results suggest that the NSVs can also shape the superior phenotypes of Nanchukmacdon in a tissue-specific manner.

Regulatory potentials of NSVs affecting the expression of nearby genes

Using chromatin state data obtained from a published paper [45], we examined the regulatory potentials of NSVs affecting the expression of nearby genes related to the specific phenotypes of Nanchukmacdon. We found that 643 NSVs near DEGs were located in actual regulatory regions (Additional file 2: Table S10). For example, an NSV with a length of 293 bp (Deletion at Sscrofa11.1: X:53,609,703–53,609,995) was discovered in the first intron of the androgen receptor (AR) gene which is involved in the diverse developmental terms (Fig. 3a). This NSV was mapped flanking to the active promoter region (TssAFlnk) in all three tissues (Additional file 2: Table S10). Also, as shown in Fig. 3b, strong signals were observed from all three histone modification marks containing the H3K4me3 signal, an indicator of the active promoter, in all three tissues. The RNA-seq data revealed that the AR gene was commonly down-regulated in all three tissues only in Nanchukmacdon, not in the other two breeds. These results indicate that the deletion could inhibit the AR gene expression by disrupting promoter activity.

Fig. 3
figure 3

Example NSVs potentially affecting expression of the AR (a) and SPAG9 (b) gene. Red colored lines highlight the location of NSVs. Area charts with blue, yellow and green color indicate H3K4me1, H3K4me3, and H3K27ac signals of corresponding regions respectively

Full size image

We also identified an NSV with a length of 304 bp (Deletion at Sscrofa11.1:12:27,201,206–27,201,509) in the last exon of the sperm associated antigen 9 (SPAG9) gene which was down-regulated only in the muscle tissue of Nanchukmacdon (Fig. 3b). The NSV was located in the enhancer region (EnhG) with relatively high H3K4me1 and H3K27ac signal only in the muscle tissue (Fig. 3b, Additional file 2: Table S10). Because those marks are indicators of an active enhancer, the deletion in that region by the NSV could prevent a normal expression of the SPAG9 gene only in the muscle tissue of Nanchukmacdon. The genotypes of both NSVs were not fully identified in the three breeds.

Discussion

In this study, we investigated the role of SVs on the meat quality trait by comparing the genomes of Nanchukmacdon and other pig breeds. Nanchukmacdon breed shows various unique meat quality phenotypes compared to other commercial pig breeds. Especially, the breed shows outstanding levels of intramuscular fat feature which is one of the most important factors determining the quality of meat [48]. Moreover, the recent construction of a chromosome-level genome assembly of Nanchukmacdon, called the NCMD assembly [19], has provided the full genomic structure of the breed. The high-quality genome assembly enables to accurately identify large and complex SVs [16, 17]. These make Nanchukmacdon a good animal model for studying the impacts of SVs on meat quality.

Here, using the NCMD assembly, we confirmed two different scales of SVs within pig species and the related biological features. For example, we found several large-scale genomic rearrangements among pig breeds through the collinearity analysis of chromosome-level assemblies based on the synteny. Especially, even though Duroc genome assembly was used to construct the NCMD assembly, we successfully revealed the rearrangement events between Nanchukmacdon and Duroc. In addition, we identified that there are several paralogous genes, such as igLON, close to the breakpoint regions. It suggests that the non-allelic homologous recombination (NAHR) between paralogous genes is one of the mechanisms of genomic rearrangements in pig species. Also, we identified 30,784 SVs including deletions, duplications, insertions and inversions in the Nanchukmacdon genome against Duroc pig. In particular, we discovered 16,898 NSVs that may contribute to the unique meat quality phenotypes of the Nanchukmacdon breed through comparison with the SVs of the Landrace breed.

Furthermore, the impacts of NSVs on the phenotypes were further interrogated with diverse multi-omics data. We identified the 414, 473 and 558 DEGs containing or close to the NSVs in the backfat, liver and muscle tissue of Nanchukmacdon against the European domestic pigs. Especially, these DEGs were related to several development-related GO terms in all three tissues. The development-related terms have been reported to affect intramuscular fat (IMF) traits in many previous studies [49,50,51]. We also identified the terms that appeared in a specific tissue, for example, cell adhesion, actin filament, and phosphorylation process for liver, and macromolecular modification, cell and blood circulation system development for muscle tissue. The actin filaments perform diverse cellular functions including cell adhesion and directly or indirectly regulate lipid synthesis and metabolism [52]. Also, the phosphorylation has been reported to be involved in the hepatic lipid metabolism [53]. The blood circulation also has been reported to affect IMF traits in chickens [54]. Interestingly, while lipid biosynthesis mainly occurs in the adipose tissue of pigs, lipids are mainly synthesized in liver tissue and transported via blood circulation in chickens [54, 55]. The significant GO terms related to lipid metabolism in liver and blood circulation in muscle suggest that the lipid biosynthesis mechanism like chickens can exist in Nanchukmacdon and affect its high IMF content. The cell development-related terms also seem to be involved in the superior meat quality, such as reddish color [56].

We also found the impact of NSVs on the expression of flanking genes investigated with the pig histone modification profile data and confirmed that several NSVs have functional regulatory effects on flanking genes, for example AR and SPAG9. The AR gene involved in the diverse development terms is commonly down-regulated in the all three tissues of Nanchukmacdon compared to the different pig breeds. The AR gene generates the Androgen receptor protein (AR), and the AR signaling pathway regulates the body fat mass [57, 58] and muscle development [59, 60]. Especially, the AR-knockout mouse models show lipogenesis and increased fat mass [60]. In addition, the sex hormones including androgen, which are ligands of the AR protein, are also reported to be associated with the amount and distribution of body fat [57, 61, 62]. SPAG9 promotes the MAPK signaling pathway which regulates lipid deposition and metabolism [63,64,65]. Also, it interacts with the NRP1 gene involved in adipogenesis [66]. We suggest that the NSVs in these genes could shape the outstanding levels of intramuscular fat deposition and redness in meat of Nanchukmacdon by affecting expression of their nearby genes.

Unfortunately, the genotyping of the NSVs using short reads was not possible in this study, which may be caused by the limitation of short reads for genotyping [67]. In addition, the specificity of NSVs on Nanchukmacdon was not experimentally validated. Thus, the genotyping using other types of sequencing data, such as long reads, and the experimental validation of NSVs by PCR or gene expression assays will be performed as a future research direction.

Conclusion

In this study, we performed the assembly-based variant calling using the chromosome-level genome assemblies of multiple pig breeds, and identified various SVs and their diverse biological features related to genome structure and rearrangements in pigs. Furthermore, the potential SVs shaping the unique meat phenotypes of Nanchukmacdon, called NSVs, were detected through comparative analysis of SVs among different pig breeds. The impact of the NSVs on the meat quality trait was further investigated by analyses using diverse multi-omics data. These findings contribute to a better understanding of the biological mechanisms of SVs underlying the meat quality trait in pigs. Also, integrative analyses of the assembly-based variant calling with diverse multi-omics data show the great potential to reveal SVs related to unique and breed-specific phenotypes.

Availability of data and materials

For comparing genome assemblies of different pig breeds, the genome assembly of the Meishan pig (MSCAAS v1) was obtained from the NCBI database (accession GCA_017957985.1). The public genome assembly of Duroc and Landrace was obtained from the NCBI database (accession GCF_000003025.6 and GCA_001700215.1, respectively), and the chromosome-level NCMD assembly was obtained from a recent study [16]. The called NSVs were annotated using VEP [27] with the pig reference annotation information (Sscrofa11.1.101) downloaded from NCBI database.

Abbreviations

SV:

Structural variant

NSV:

Nanchukmacdon-specific structural variant

log2FC:

Log2 fold change

IMF:

Intramuscular fat

DEG:

Differentially expressed gene

DGE analysis:

Differential gene expression analysis

NAHR:

Non-allelic homologous recombination

References

  1. Bonnet C, Bouamra-Mechemache Z, Requillart V, Treich N. Viewpoint: Regulating meat consumption to improve health, the environment and animal welfare. Food Policy. 2020;97:Article 101847.

  2. Barbosa MLV, Hauser-Davis RA, Siciliano S, Dias TLP, Alves RRN, Costa-Neto EM. Historical Shark Meat Consumption and Trade Trends in a Global Richness Hotspot. Ethnobiol Lett. 2019;10(1):97–103.

    Article  Google Scholar 

  3. Mortimer SI, Przybylski W. Breeding strategies for improving meat quality. Meat Quality Genetic and Environmental Factors. Boca Raton: CRC Press; 2016. p. 257–85.

  4. Warr A, Affara N, Aken B, Beiki H, Bickhart DM, Billis K, Chow W, Eory L, Finlayson HA, Flicek P, et al. An improved pig reference genome sequence to enable pig genetics and genomics research. Gigascience. 2020;9(6):giaa051.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Lee D, Cho M, Hong WY, Lim D, Kim HC, Cho YM, Jeong JY, Choi BH, Ko Y, Kim J. Evolutionary Analyses of Hanwoo (Korean Cattle)-Specific Single-Nucleotide Polymorphisms and Genes Using Whole-Genome Resequencing Data of a Hanwoo Population. Mol Cells. 2016;39(9):692–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Mei C, Wang H, Liao Q, Wang L, Cheng G, Wang H, Zhao C, Zhao S, Song J, Guang X, et al. Genetic Architecture and Selection of Chinese Cattle Revealed by Whole Genome Resequencing. Mol Biol Evol. 2018;35(3):688–99.

    Article  CAS  PubMed  Google Scholar 

  7. Lee D, Lee J, Heo KN, Kwon K, Moon Y, Lim D, Lee KT, Kim J. Population analysis of the Korean native duck using whole-genome sequencing data. BMC Genomics. 2020;21(1):554.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Bedhane M, van der Werf J, Gondro C, Duijvesteijn N, Lim D, Park B, Park MN, Hee RS, Clark S. Genome-Wide Association Study of Meat Quality Traits in Hanwoo Beef Cattle Using Imputed Whole-Genome Sequence Data. Front Genet. 2019;10:1235.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Gao G, Gao N, Li S, Kuang W, Zhu L, Jiang W, Yu W, Guo J, Li Z, Yang C, Zhao Y. Genome-Wide Association Study of Meat Quality Traits in a Three-Way Crossbred Commercial Pig Population. Front Genet. 2021;12:614087.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Guo Q, Huang L, Bai H, Wang Z, Bi Y, Chen G, Jiang Y, Chang G. Genome-Wide Association Study of Potential Meat Quality Trait Loci in Ducks. Genes (Basel). 2022;13(6):986.

    Article  CAS  PubMed  Google Scholar 

  11. Wold J, Koepfli KP, Galla SJ, Eccles D, Hogg CJ, Le Lec MF, Guhlin J, Santure AW, Steeves TE. Expanding the conservation genomics toolbox: Incorporating structural variants to enhance genomic studies for species of conservation concern. Mol Ecol. 2021;30(23):5949–65.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Collins RL, Brand H, Karczewski KJ, Zhao X, Alfoldi J, Francioli LC, Khera AV, Lowther C, Gauthier LD, Wang H, et al. A structural variation reference for medical and population genetics. Nature. 2020;581(7809):444–51.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Weischenfeldt J, Symmons O, Spitz F, Korbel JO. Phenotypic impact of genomic structural variation: insights from and for human disease. Nat Rev Genet. 2013;14(2):125–38.

    Article  CAS  PubMed  Google Scholar 

  14. Tattini L, D’Aurizio R, Magi A. Detection of Genomic Structural Variants from Next-Generation Sequencing Data. Front Bioeng Biotechnol. 2015;3:92.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Mahmoud M, Gobet N, Cruz-Davalos DI, Mounier N, Dessimoz C, Sedlazeck FJ. Structural variant calling: the long and the short of it. Genome Biol. 2019;20(1):246.

    Article  PubMed  PubMed Central  Google Scholar 

  16. De Coster W, Van Broeckhoven C. Newest Methods for Detecting Structural Variations. Trends Biotechnol. 2019;37(9):973–82.

    Article  PubMed  Google Scholar 

  17. Lee H, Kim J, Lee J. Benchmarking datasets for assembly-based variant calling using high-fidelity long reads. BMC Genomics. 2023;24(1):148.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Wala JA, Bandopadhayay P, Greenwald NF, O’Rourke R, Sharpe T, Stewart C, Schumacher S, Li Y, Weischenfeldt J, Yao X, et al. SvABA: genome-wide detection of structural variants and indels by local assembly. Genome Res. 2018;28(4):581–91.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Kwon D, Park N, Wy S, Lee D, Chai HH, Cho IC, Lee J, Kwon K, Kim H, Moon Y, et al. A chromosome-level genome assembly of the Korean crossbred pig Nanchukmacdon (Sus scrofa). Sci Data. 2023;10(1):761.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Arora D, Srikanth K, Lee J, Lee D, Park N, Wy S, Kim H, Park JE, Chai HH, Lim D, et al. Integration of multi-omics approaches for functional characterization of muscle related selective sweep genes in Nanchukmacdon. Sci Rep. 2021;11(1):7219.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Harris RS: Improved pairwise alignment of genomic DNA: The Pennsylvania State University; 2007.

  22. Ma J, Zhang LX, Suh BB, Raney BJ, Burhans RC, Kent WJ, Blanchette M, Haussler D, Miller W. Reconstructing contiguous regions of an ancestral genome. Genome Res. 2006;16(12):1557–65.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Lee J, Lee D, Sim M, Kwon D, Kim J, Ko Y, Kim J. mySyntenyPortal: an application package to construct websites for synteny block analysis. BMC Bioinformatics. 2018;19(1):216.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Zhou R, Li ST, Yao WY, Xie CD, Chen Z, Zeng ZJ, Wang D, Xu K, Shen ZJ, Mu Y, et al. The Meishan pig genome reveals structural variation-mediated gene expression and phenotypic divergence underlying Asian pig domestication. Mol Ecol Resour. 2021;21(6):2077–92.

    Article  CAS  PubMed  Google Scholar 

  25. Li H, Durbin R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009;25(14):1754–60.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Li H. Minimap2: pairwise alignment for nucleotide sequences. Bioinformatics. 2018;34(18):3094–100.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Servant N, Varoquaux N, Lajoie BR, Viara E, Chen CJ, Vert JP, Heard E, Dekker J, Barillot E. HiC-Pro: an optimized and flexible pipeline for Hi-C data processing. Genome Biol. 2015;16:259.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Akdemir KC, Chin L. HiCPlotter integrates genomic data with interaction matrices. Genome Biol. 2015;16:198.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Cingolani P. Variant Annotation and Functional Prediction: SnpEff. Methods Mol Biol. 2022;2493:289–314.

    Article  CAS  PubMed  Google Scholar 

  30. Heller D, Vingron M. SVIM-asm: Structural variant detection from haploid and diploid genome assemblies. Bioinformatics. 2020;36(22–23):5519–21.

    CAS  PubMed Central  Google Scholar 

  31. McLaren W, Gil L, Hunt SE, Riat HS, Ritchie GR, Thormann A, Flicek P, Cunningham F. The Ensembl Variant Effect Predictor. Genome Biol. 2016;17(1):122.

    Article  PubMed  PubMed Central  Google Scholar 

  32. Reimand J, Kull M, Peterson H, Hansen J, Vilo J. g Profiler–a web-based toolset for functional profiling of gene lists from large-scale experiments. Nucleic Acids Res. 2007;35(Web Server issue):193–200.

    Article  Google Scholar 

  33. Andrews S: FastQC: a quality control tool for high throughput sequence data. In.: Babraham Bioinformatics, Babraham Institute, Cambridge, United Kingdom; 2010.

  34. Patel RK, Jain M. NGS QC Toolkit: a toolkit for quality control of next generation sequencing data. PLoS ONE. 2012;7(2):e30619.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Krueger F: Trim Galore: a wrapper tool around Cutadapt and FastQC to consistently apply quality and adapter trimming to FastQ files, with some extra functionality for MspI-digested RRBS-type (Reduced Representation Bisufite-Seq) libraries. URL http://www.bioinformaticsbabrahamacuk/projects/trim_galore/(Accessed date: 28/04/2016) 2012.

  36. Eggertsson HP, Jonsson H, Kristmundsdottir S, Hjartarson E, Kehr B, Masson G, Zink F, Hjorleifsson KE, Jonasdottir A, Jonasdottir A, et al. Graphtyper enables population-scale genotyping using pangenome graphs. Nat Genet. 2017;49(11):1654–60.

    Article  CAS  PubMed  Google Scholar 

  37. Yang J, Lee SH, Goddard ME, Visscher PM. GCTA: a tool for genome-wide complex trait analysis. Am J Hum Genet. 2011;88(1):76–82.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Ginestet C. ggplot2: Elegant Graphics for Data Analysis. J R Stat Soc a Stat. 2011;174:245–245.

    Article  Google Scholar 

  39. Arora D, Park J-E, Lim D, Choi B-H, Cho I-C, Srikanth K, Kim J, Park W. Comparative Methylation and RNA-Seq Expression Analysis in CpG Context to Identify Genes Involved in Backfat vs Liver Diversification in Nanchukmacdon Pig. BMC Genomics. 2021;22(1):801.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Andrew S. A quality control tool for high throughput sequence data. Fast QC. 2010;532:1.

    Google Scholar 

  41. Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30(15):2114–20.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Kim D, Langmead B, Salzberg SL. HISAT: a fast spliced aligner with low memory requirements. Nat Methods. 2015;12(4):357–60.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Liao Y, Smyth GK, Shi W. featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics. 2014;30(7):923–30.

    Article  CAS  PubMed  Google Scholar 

  44. Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15(12):550.

    Article  PubMed  PubMed Central  Google Scholar 

  45. Pan ZY, Yao YL, Yin HW, Cai ZX, Wang Y, Bai LJ, Kern C, Halstead M, Chanthavixay G, Trakooljul N, et al. Pig genome functional annotation enhances the biological interpretation of complex traits and human disease. Nat Commun. 2021;12(1):5848.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Audano PA, Sulovari A, Graves-Lindsay TA, Cantsilieris S, Sorensen M, Welch AE, Dougherty ML, Nelson BJ, Shah A, Dutcher SK, et al. Characterizing the Major Structural Variant Alleles of the Human Genome. Cell. 2019;176(3):663–75.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. He Y, Luo X, Zhou B, Hu T, Meng X, Audano PA, Kronenberg ZN, Eichler EE, Jin J, Guo Y, et al. Long-read assembly of the Chinese rhesus macaque genome and identification of ape-specific structural variants. Nat Commun. 2019;10(1):4233.

    Article  PubMed  PubMed Central  Google Scholar 

  48. Khan R, Li A, Raza SHA. Editorial: Genetic Regulation of Meat Quality Traits in Livestock Species. Front Genet. 2022;13:1092562.

    Article  PubMed  Google Scholar 

  49. Lee HJ, Jang M, Kim H, Kwak W, Park W, Hwang JY, Lee CK, Jang GW, Park MN. Kim HC et al Comparative Transcriptome Analysis of Adipose Tissues Reveals that ECM-Receptor Interaction Is Involved in the Depot-Specific Adipogenesis in Cattle. Plos One. 2013;8(6):e66267.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. Ding RR, Yang M, Quan JP, Li SY, Zhuang ZW, Zhou SP, Zheng EQ, Hong LJ, Li ZC, Cai GY, et al. Single-Locus and Multi-Locus Genome-Wide Association Studies for Intramuscular Fat in Duroc Pigs. Frontiers in Genetics. 2019;10:619.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Ovilo C, Benitez R, Fernandez A, Nunez Y, Ayuso M, Fernandez AI, Rodriguez C, Isabel B, Rey AI, Lopez-Bote C, Silio L. Longissimus dorsi transcriptome analysis of purebred and crossbred Iberian pigs differing in muscle characteristics. Bmc Genomics. 2014;15:413.

    Article  PubMed  PubMed Central  Google Scholar 

  52. Liu GY, Li JY, Wu CY. Reciprocal regulation of actin filaments and cellular metabolism. Eur J Cell Biol. 2022;101(4):151281.

    Article  CAS  PubMed  Google Scholar 

  53. Ramatchandirin B, Pearah A, He L. Regulation of Liver Glucose and Lipid Metabolism by Transcriptional Factors and Coactivators. Life (Basel). 2023;13(2):515.

    CAS  PubMed  Google Scholar 

  54. Cui HX, Liu RR, Zhao GP, Zheng MQ, Chen JL, Wen J. Identification of differentially expressed genes and pathways for intramuscular fat deposition in pectoralis major tissues of fast-and slow-growing chickens. Bmc Genomics. 2012;13:213.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  55. Leveille GA. In vitro hepatic lipogenesis in the hen and chick. Comp Biochem Physiol. 1969;28(1):431–5.

    Article  CAS  PubMed  Google Scholar 

  56. Listrat A, Lebret B, Louveau I, Astruc T, Bonnet M, Lefaucheur L, Picard B, Bugeon J. How muscle structure and composition influence meat and flesh quality. The Scientific World Journal. 2016;2016:3182746.

    Article  PubMed  PubMed Central  Google Scholar 

  57. Venkatesh VS, Grossmann M, Zajac JD, Davey RA. The role of the androgen receptor in the pathogenesis of obesity and its utility as a target for obesity treatments. Obes Rev. 2022;23(6):e13429.

    Article  PubMed  PubMed Central  Google Scholar 

  58. McInnes KJ, Smith LB, Hunger NI, Saunders PT, Andrew R, Walker BR. Deletion of the androgen receptor in adipose tissue in male mice elevates retinol binding protein 4 and reveals independent effects on visceral fat mass and on glucose homeostasis. Diabetes. 2012;61(5):1072–81.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  59. Zhang Q, Lee HG, Han JA, Kang SK, Lee NK, Baik M, Choi YJ. Differentially expressed proteins associated with myogenesis and adipogenesis in skeletal muscle and adipose tissue between bulls and steers. Mol Biol Rep. 2012;39(2):953–60.

    Article  CAS  PubMed  Google Scholar 

  60. Rizk J, Sahu R, Duteil D. An overview on androgen-mediated actions in skeletal muscle and adipose tissue. Steroids. 2023;199: 109306.

    Article  CAS  PubMed  Google Scholar 

  61. Hamilton EJ, Gianatti E, Strauss BJ, Wentworth J, Lim-Joon D, Bolton D, Zajac JD, Grossmann M. Increase in visceral and subcutaneous abdominal fat in men with prostate cancer treated with androgen deprivation therapy. Clin Endocrinol (Oxf). 2011;74(3):377–83.

    Article  CAS  PubMed  Google Scholar 

  62. Jones TH. Testosterone deficiency: a risk factor for cardiovascular disease? Trends Endocrinol Metab. 2010;21(8):496–503.

    Article  CAS  PubMed  Google Scholar 

  63. Xiao C, Li M, Huang Q, Si-Tu J. SPAG9 promotes prostate cancer proliferation and metastasis via MAPK signaling pathway. Am J Transl Res. 2019;11(8):5249–60.

    CAS  PubMed  PubMed Central  Google Scholar 

  64. Shen Y, Zhao Y, Zheng D, Chang X, Ju S, Guo L. Effects of orexin A on GLUT4 expression and lipid content via MAPK signaling in 3T3-L1 adipocytes. J Steroid Biochem Mol Biol. 2013;138:376–83.

    Article  CAS  PubMed  Google Scholar 

  65. Gubern A, Barcelo-Torns M, Barneda D, Lopez JM, Masgrau R, Picatoste F, Chalfant CE, Balsinde J, Balboa MA, Claro E. JNK and ceramide kinase govern the biogenesis of lipid droplets through activation of group IVA phospholipase A2. J Biol Chem. 2009;284(47):32359–69.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  66. Yu Y, Uchida-Fukuhara Y, Weng Y, He Y, Ikegame M, Wang Z, Yoshida K, Okamura H, Qiu L. Neuropilin 1 (NRP1) Positively Regulates Adipogenic Differentiation in C3H10T1/2 Cells. Int J Mol Sci. 2023;24(8):7394.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  67. De Coster W, Weissensteiner MH, Sedlazeck FJ. Towards population-scale long-read sequencing. Nat Rev Genet. 2021;22(9):572–87.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This work was supported by Agenda (project No. PJ01334302) of the National Institute of Animal Science funded by the Rural Development Administration (RDA) and the Ministry of Science and ICT (NRF-2021M3H9A2097134, NRF-2022R1F1A1065159).

Author information

Authors and Affiliations

  1. Department of Biomedical Science and Engineering, Konkuk University, Seoul, 05029, Republic of Korea

    Daehong Kwon, Nayoung Park, Suyeon Wy, Daehwan Lee, Jongin Lee, Kisang Kwon, Heesun Kim, Youngbeen Moon, Juyeon Kim & Jaebum Kim

  2. Animal Genomics and Bioinformatics Division, National Institute of Animal Science, RDA, Wanju, 55365, Republic of Korea

    Woncheoul Park & Han-Ha Chai

  3. Subtropical Livestock Research Institute, National Institute of Animal Science, RDA, Jeju, 63242, Republic of Korea

    In-Cheol Cho

Authors
  1. Daehong Kwon
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Nayoung Park
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Suyeon Wy
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Daehwan Lee
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Woncheoul Park
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Han-Ha Chai
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. In-Cheol Cho
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Jongin Lee
    View author publications

    You can also search for this author in PubMed Google Scholar

  9. Kisang Kwon
    View author publications

    You can also search for this author in PubMed Google Scholar

  10. Heesun Kim
    View author publications

    You can also search for this author in PubMed Google Scholar

  11. Youngbeen Moon
    View author publications

    You can also search for this author in PubMed Google Scholar

  12. Juyeon Kim
    View author publications

    You can also search for this author in PubMed Google Scholar

  13. Jaebum Kim
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

J.B.K. conceived and supervised the study. W.C.P, I.C.C and H.H.C prepared the sequencing data. J.B.K, D.H.K., N.Y.P., S.Y.W., D.H.L. and J.I.L. designed analyses. D.H.K., N.Y.P., S.Y.W., H.S.K., D.H.L., J.I.L. and Y.B.M., J.Y.K., K.S.K. performed analyses. J.B.K, D.H.K., N.Y.P. and S.Y.W. interpreted results. J.B.K, D.H.K, N.Y.P., S.Y.W. and H.S.K. wrote the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Jaebum Kim.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1. 

Supplementary Figures S1-S8

Additional file 2. 

Supplemental Tables S1-S10

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Kwon, D., Park, N., Wy, S. et al. Identification and characterization of structural variants related to meat quality in pigs using chromosome-level genome assemblies. BMC Genomics 25, 299 (2024). https://doi.org/10.1186/s12864-024-10225-1

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12864-024-10225-1

Keywords

BMC Genomics

ISSN: 1471-2164

Contact us