Prerpints.org logo
Preprint
Article

Comprehensive Atlas of Alternative Splicing Reveals NSRP1 Promoting Adipogenesis through CCDC18

Altmetrics

Downloads

108

Views

28

Comments

0

A peer-reviewed article of this preprint also exists.

supplementary.zip (60.91KB )

Submitted:

19 January 2024

Posted:

23 January 2024

You are already at the latest version

Alerts
Abstract
Alternative splicing (AS) plays a crucial role in regulating gene expression, function and diversity. However, limited reports exist on identification and comparation of AS in Eastern and Western pigs. Here, we analyzed 243 transcriptome data from eight tissues, integrating information on transcription factors (TFs), selection signals, splicing factors (SFs), and QTLs to comprehensively study alternative splicing events (ASEs) in pigs. Five ASE types were identified, with Mutually Exclusive Exon (MXE) and Skipped Exon (SE) ASEs being the most prevalent. A significant portion of genes with ASEs (ASGs) showed conservation across all eight tissues (63.21% - 76.13% per tissue). Differentially alternative splicing genes (DASGs) and differentially expressed genes (DEGs) exhibited tissue specificity, with blood and adipose tissues having more DASGs. Functional enrichment analysis revealed coDASG_DEGs in adipose were enriched in pathways associated with adipose deposition and immune inflammation, while coDASG_DEGs in blood were enriched in pathways related to immune inflammation and metabolism. Adipose deposition in Eastern pigs might be linked to down-regulation of immune-inflammation-related pathways and reduced insulin resistance. The TFs, selection signals, and SFs appeared to regulate ASEs. Notably, ARID4A (TF), NSRP1 (SF), ANKRD12, IFT74, KIAA2026, CCDC18, NEXN, PPIG, and ROCK1 genes in adipose tissue showed potential regulatory effects on adipose-deposition traits. And NSRP1 could promote adipogenesis by regulating alternative splicing and expression of CCDC18. Conducting an in-depth investigation into AS, this study has successfully identified key marker genes essential for pig genetic breeding and the enhancement of meat quality, which will play important roles in promoting the diversity of pork quality and meeting market demand.
Keywords: 
Subject: Biology and Life Sciences  -   Endocrinology and Metabolism

1. Introduction

As an essential livestock animal, pigs provide abundant food resources for human beings and hold significant symbolic importance in various cultures. The Eastern and Western pig populations, as the main subtypes, have distinct ecological environments and developmental histories [1]. Eastern pigs are mainly distributed in Asia, including China, Japan, South Korea and other countries, and their history of communication with humans can be traced back tens of thousands of years ago [2,3]. On the other hand, Western pigs originated in Europe and have spread to other continents through human exploration and expansion [4]. Throughout their long domestication process, these two pig types have developed distinct genetic characteristics and ecological habits [4,5]. For example, Western pigs have a faster growth rate, while Eastern pigs have stronger disease resistance and fat content [6,7,8].
The AS of pre-RNA is a crucial transcriptional regulatory mechanism involving splicing mRNA precursors of genes in different ways. This process allows a single gene to generate multiple transcripts, translating diverse proteins and significantly enriching the transcriptome and proteome diversity [9,10]. AS plays a pivotal role in organism growth, development, and cell differentiation [11] and is widespread in mammals [12]. In humans, 40-60% of genes undergo AS, enhancing the coding capacity of the genome [13]. Furthermore, AS is linked to numerous disease factors [14], underscoring its critical importance. There are many factors that influence AS, such as transcription factors, selection signals, splicing factors, etc [15,16,17,18].
Studying AS in pigs will deepen our understanding of pig gene functionality and regulatory mechanisms and sheds light on their relationships and evolutionary processes with other species. Moreover, as society progresses and concerns about food safety and environmental protection grow, pigs have garnered increased attention as a major livestock animal [19]. Investigating AS characteristics in Eastern and Western pigs can offer valuable insights into pig domestication and breeding, leading to improvements in productivity and adaptability and contributing to the sustainable development of the pork industry. Researchers have identified widespread pre-RNA AS occurrences in multiple pig tissues [20,21,22]. However, the differences in AS among various tissues of Eastern and Western pigs and their molecular regulatory mechanisms remain unclear.
In this study, we collected a total of 243 transcriptome datasets from eight tissues (adipose, blood, heart, lung, kidney, muscle, ovary, and spleen) in Eastern and Western pigs from public databases. By integrating information on transcription factors (TFs), selection signals, splicing factors (SFs), and QTLs, we comprehensively compared AS characteristics between the two pig populations. The results revealed tissue-specific differences in both differentially alternative splicing genes (DASGs) and differentially expressed genes (DEGs), potentially influenced by TFs, selection signals, and SFs. Furthermore, by incorporating QTLs information, we identified marker genes associated with meat quality traits such as backfat thickness and intermuscular fat deposition. Our research aims to provide scientific evidence for the conservation and rational utilization of Eastern and Western pig genetic resources. Additionally, we seek to establish a new theoretical foundation for pig genetic breeding and the improvement of meat quality traits, thereby contributing to the sustainable development of the livestock industry and deepening our understanding of pig gene function and regulatory mechanisms.

2. Results

2.1. Genome-wide identification of ASEs in different tissues

To investigate the differences in alternative splicing events (ASEs) across various tissues of Eastern and Western pigs, we conducted identification and differential analysis of five types of ASEs in the eight tissues, including A3SS, A5SS, MXE, RI, and SE (Figure 1A). The results revealed that MXE and SE types of ASEs and differentially alternative splicing events (DASEs) had the highest occurrence in each tissue. Among them, SE ASEs were the most abundant, while MXE had the highest number of DASEs (Figure 1B). In terms of different tissues, blood and adipose tissues exhibited a higher number of ASEs and DASEs. Additionally, these two tissues had the highest proportion of DASEs, with adipose showing 1.5 times (3.51% / 2.34%) more than that of blood (Figure 1B). We further counted the number of genes with alternative splicing (ASGs) and found that the percentage of genes undergoing AS in each tissue ranged from 32.59% to 39.25% of the genome-wide genes (Figure 1B), indicating that the number of genes undergoing AS in each tissue was relatively stable. Moreover, we observed that the proportion of DASGs was highest in adipose and blood tissues (33.43% and 25.74%, respectively), with adipose showing 1.3 times more than blood (Figure 1B). This finding suggested that AS plays a significant role in adipose deposition of pigs. Furthermore, we investigated conserved genes with AS (ASGs) in different tissues and identified 7917 common ASGs present in all eight tissues, accounting for 63.21% to 76.13% of ASGs in each tissue (Figure 1C). This finding suggested that genes undergoing AS in various tissues display a certain level of conservation.

2.2. Tissue-specific DASGs and DEGs

We further investigated the tissue-specificity of differentially alternative splicing genes (DASGs) and differentially expressed genes (DEGs) in each tissue. The results revealed that, in the comparison between Western and Eastern pigs, both DASGs and DEGs exhibited tissue-specificity, with adipose and blood tissues showing higher tissue-specificity for both categories. Among all tissues, the ovary had the highest tissue-specificity for DEGs (65.28%), but a relatively lower tissue-specificity for DASGs (17.29%), suggesting that gene expression in the ovary might be less influenced by AS (Table 1). Furthermore, we calculated the ratio of DASGs to DEGs in each tissue and found that, on average, each gene underwent 4.98, 3.02, 0.97, 1.11, 0.88, 0.79, 0.43, and 2.78 ASEs in adipose, blood, heart, kidney, lung, muscle, ovary, and spleen tissues, respectively. These results further indicated that genes in adipose and blood tissues were more likely to undergo AS compared to other tissues (Table 1), and we focused on these two tissues in the following research. We identified 840 and 960 DEGs for adipose and blood, respectively, and filtered down to 83 and 164 coDASG_DEGs for adipose and blood, respectively (Table 1). These coDASG_DEGs may be significantly influenced by ASEs and represent promising candidates for further investigation.

2.3. Functional analysis of DEGs association with ASEs

To investigate the impact of alternative splicing events (ASEs) on the biological functions in adipose and blood tissues, we performed KEGG pathway enrichment analysis on coDASG_DEGs specific to these two tissues. In adipose tissue, we identified a total of 29 significantly enriched pathways (P < 0.05). Among them, the Ras signaling pathway, FoxO signaling pathway, PI3K-Akt signaling pathway, insulin resistance, MAPK signaling pathway, and AMPK signaling pathway were related to adipose deposition, involving genes such as IGF1, MLXIPL, NTRK2, IL6, KIT, TNFSF10, PFKFB1, and PAK3 (Figure 2A). Notably, genes IGF1 and PFKFB1 in the AMPK signaling pathway were significantly upregulated in Eastern pigs, while genes MLXIPL and IL6 in the insulin resistance pathway were significantly downregulated (Figure 2A). Furthermore, we observed that genes related to immune inflammation pathways, such as Non-alcoholic fatty liver disease, Legionellosis, IL-17 signaling pathway, Rheumatoid arthritis, Amoebiasis, and TNF signaling pathway were significantly downregulated, involving MLXIPL, CXCL2, KIT, NDUFS4, and IL6 (Figure 2A). Moreover, we conducted Protein-Protein Interaction (PPI) analysis on coDASG_DEGs using the STRING database. We applied K-means clustering and divided the adipose tissue into two clusters. Cluster 1 in adipose tissue was enriched in insulin stimulus-related pathways (Figure 2B), further suggesting a possible association between adipose deposition in Eastern pigs and insulin resistance.
In blood, we identified a total of seven significantly enriched KEGG pathways (P < 0.05), mainly related to immune inflammation and metabolism. All coDASG_DEGs in pathways such as Porphyrin metabolism, Rheumatoid arthritis, Cytosolic DNA-sensing pathway, and Autophagy – other showed significant downregulation in Eastern pigs, involving ATG4A, ATP6V1E1, BLVRB, CCL5, FECH, GABARAP, HMBS, IL15, IL18, IRF7, and UROS (Figure 2C). Similarly, we conducted PPI analysis on coDASG_DEGs using the STRING database for blood tissue. We applied K-means clustering and divided the blood tissue into three clusters. Cluster 1 was enriched in immune inflammation-related pathways, while Cluster 2 and Cluster 3 were primarily enriched in metabolism-related pathways (Figure 2D). These pathways and genes suggested differences in disease resistance and metabolic levels between Eastern and Western pigs.

2.4. Effects of transcription factors on alternative splicing

Researches show that transcription factors (TFs) can regulate the splicing process of genes [15,16]. To investigate the regulatory role of TFs on gene AS, we conducted a joint analysis of DEGs, DASGs, and TFs. We identified 63 and 75 TFs in the DEGs (coDEG_TFs) of adipose and blood, respectively (Figure 3A,B), accounting for 7.50% and 7.81% of adipose and blood DEGs. Among them, in adipose tissue, 12 coDEG_TFs also exhibited differentially alternative splicing between Eastern and Western pigs, including ARID4A, BBX, BCL6B, CREBRF, FOSL2, HMBOX1, MLXIPL, NFIB, NR1D2, SON, SOX5, and ZNF148 (Figure 3A). In terms of blood tissue, 16 coDEG_TFs also showed differentially alternative splicing between Eastern and Western pigs, including ETV7, IRF1, IRF7, KMT2C, LITAF, MAX, MGA, NCOR2, NFAT5, PRDM11, SOX13, TET2, TFDP2, USF2, YBX3, and ZMIZ1 (Figure 3A).
We conducted Pearson correlation analysis to study the regulatory effect of these coDEG_TFs on coDASG_DEGs gene expression. Using the criteria r > 0.9 and P < 3.4E-14, we identified 110 TF-gene pairs, involving 24 coDEG_TFs and 38 (45.78%) coDASG_DEGs in adipose tissue (Figure 3C). In blood tissue, we found 350 TF-gene pairs, involving 46 coDEG_TFs and 86 (52.44%) coDASG_DEGs (Figure 3D). These results indicate that approximately half of the coDASG_DEGs are significantly influenced by TFs, revealing the important regulatory role of TFs in gene AS.

2.5. The impact of genes under potential selection on alternative splicing

The differences in ASEs between Eastern and Western pigs may be influenced by long-term artificial selection during domestication and breeding. Therefore, we also investigated the impact of selection signals on AS. Based on the previously identified genes under selection (SGs) in Eastern and Western pigs [23], we further explored their influence on AS. In both adipose and blood tissues, we identified 20 genes under selection intersecting with coDASG_DEGs (Figure 4). Among these, five TFs (SOX5, ZNF148, CREBRF, ARID4A, SON) were involved in adipose tissue (Figure 4A), and two TFs (MGA, KMT2C) were involved in blood tissue (Figure 4B).
In adipose tissue, we found that SOX5, ENSSSCG00000014414, PPP1R12A, ZNF148, TNFSF10, ROCK1, NKD2, RUSC2, and CXCL2 were under selection in Eastern pigs (Figure 4A). Among them, SOX5, ENSSSCG00000014414, PPP1R12A, ZNF148, TNFSF10, and ROCK1 were significantly upregulated in Eastern pig's adipose tissue (Figure 4A), potentially influencing adipose deposition in pigs. In blood tissue, MACF1, IL1RL1, MGA, TTC14, ANKRD17, PDLIM7, FGD3, NPRL3, ENSSSCG00000014060, ENSSSCG00000039214, CYB5A, and HAGH were under selection in Eastern pigs (Figure 4B). Among them, MACF1, IL1RL1, MGA, TTC14, and ANKRD17 were significantly upregulated in Eastern pig's blood tissue (Figure 4B), potentially affecting disease resistance in pigs.

2.6. Effects of splicing factors on alternative splicing

Splicing factors (SFs) are protein factors involved in the RNA precursor splicing process [24]. To study the regulatory role of SFs on alternative splicing, we identified five SFs (SERP1, MBNL3, NOVA1, NSRP1, and SRSF12) from DEGs in adipose tissue. Then, we conducted Pearson correlation analysis between the gene expression levels of these five SFs and coDASG_DEGs. The results showed that out of 415 SF-coDASG_DEG pairs, 263 (63.37%) pairs exhibited significant correlations (|r| > 0.324, P < 0.05, Figure 5). It indicates that these five SFs play a crucial role in the alternative splicing of DEGs in the adipose tissue of Eastern and Western pigs.
Furthermore, among the 83 coDASG_DEGs, we found that 74 (89.16%) genes exhibited consistent correlation patterns with these five SFs, with 44 genes showing positive correlations and 30 genes showing negative correlations. It suggested that these SFs may play a co-regulatory role in the process of AS (Figure 5). After sorting the 415 SF-coDASG_DEG pairs in descending order of correlation coefficient, we found that the top 18 pairs of genes were significantly and highly positively correlated with the splicing factor NSRP1 (r > 0.86, P < 6.104E-12, Figure 5). These genes included PPIG, IFT74, AKAP9, CCDC18, NEXN, ROCK1, AHI1, ENSSSCG00000037142, EIF5B, ARID4A, CCDC191, ANKRD12, GOLGA4, TWISTNB, KIAA2026, GKAP1, SMC4, and ENSSSCG00000009128 (Figure 5), indicating that the SF NSRP1 may have significant regulatory effects on adipose AS and pig adipose deposition.

2.7. Splicing factor NSRP1 and correlated genes impact adipose-deposition traits

To investigate the potential role of the splicing factor NSRP1 in adipose deposition, we downloaded QTLs data for pigs from PigQTLdb (https://www.animalgenome.org/cgi-bin/QTLdb/SS/index). We then conducted a combined analysis using NSRP1 and the 18 genes highly correlated with the QTLs information. As a result, 1,057 QTLs were identified, of which 235 were associated with adipose-deposition traits. These QTLs were linked to the 18 genes (excluding ENSSSCG00000037142), and were associated with 30 traits, including abdominal fat weight, adipocyte diameter, backfat thickness, and intramuscular fat content (Figure 6A, Table S1). Further analysis revealed that all these genes in adipose-related QTLs were highly expressed in the Eastern pigs (Figure 6B). Moreover, as many as 14 genes were found in the Average backfat thickness QTL, further indicating a potential link between NSRP1 and its related genes with adipose deposition. Notably, five genes, namely ANKRD12, ARID4A, IFT74, KIAA2026, and ROCK1, appeared in 12 or more QTLs, suggesting their potentially significant roles in adipose-deposition traits (Figure 6A). Additionally, ANKRD12, CCDC18, NEXN, PPIG, and ROCK1 were associated with the intramuscular fat content QTL (Figure 6A), implying their potential importance in intermuscular fat deposition and meat quality improvement. These findings suggested a possible regulatory network involving NSRP1 and its correlated genes that influence pigs' adipose deposition and meat quality traits.

2.8. Splicing factor NSRP1 regulates adipogenesis

To evaluate the influence of NSRP1 on pre-adipocyte proliferation and differentiation, we conducted interference experiments by transfecting three siRNAs targeting NSRP1. The results indicated that si2 and si3 exhibited the highest interference efficiency during both proliferation and differentiation stages (Figure 7A). Consequently, NSRP1 knockdown was achieved through transfection with si2 and si3 (Figure 7B). Subsequently, CCK-8 assays were performed to assess the rate of cell proliferation. Over time, the proliferation rate accelerated following si-NSRP1 treatment, as demonstrated by the results of the CCK-8 assay (Figure 7C). RT-qPCR and Western blot analyses confirmed an increase in the expression of proliferation marker genes, Ki67 and PCNA, during pre-adipocyte proliferation post si-NSRP1 treatment (Figure 7D,E). EdU assay results exhibited a parallel trend, indicating an increase in the percentage of EdU-positive cells post si-NSRP1 treatment (Figure 7F,G). Moreover, RT-qPCR analysis validated a reduction in the expression of differentiation marker genes, PPARG and CEBPA, during pre-adipocyte differentiation following si-NSRP1 treatment (Figure 7H). Western blot analysis further substantiated after the successful interference of NSRP1, demonstrating a decrease in the expression of differentiation marker genes PPARG and CEBPA (Figure 7I). Oil Red O staining revealed a significant reduction in lipid droplet formation following si-NSRP1 treatment (Figure 7J). In summary, following the downregulation of NSRP1 expression, there was an observed acceleration in the proliferation rate of pre-adipocytes, coupled with a diminished capacity for lipid droplet formation. This implies that the NSRP1 has the capability to inhibit the proliferation of pre-adipocytes, and additionally, it plays a role in promoting lipid droplet formation during the adipogenic stage.

2.9. NSRP1 promotes adipogenesis by regulating AS and expression of CCDC18

CCDC18 was a coDASG_DEG that was significantly highly expressed in Eastern pigs, which was found to be highly positively correlated with NSRP1 (r = 0.96, P < 1.00E-16). To validate differential alternative splicing events in CCDC18 identified through transcriptome data analysis between Eastern and Western pig breeds, we extracted adipose tissues from the representative Western breed, Duroc; and the Eastern breed, Luchuan pig. Semiquantitative RT-qPCR experiments were conducted, revealing significant differences in the alternative splicing events of CCDC18 between Eastern and Western pig breeds. Notably, the Percent Spliced In (PSI) values in Luchuan pig were markedly higher than those in Duroc (Figure 8A). Furthermore, we validated whether NSRP1, as the alternative splicing factor, could regulate the occurrence of alternative splicing in CCDC18. Semiquantitative results demonstrated a significant decrease in the PSI values of CCDC18 when NSRP1 was downregulated in pre-adipocytes (Figure 8B). This suggested that NSRP1 indeed has the capability to regulate the strength of alternative splicing events in CCDC18. Afterward, we further verified whether NSRP1 could regulate the expression of CCDC18. The RT-qPCR results revealed a significant inhibition in the mRNA of CCDC18 during both cell proliferation and differentiation stages following si-NSRP1 treatment (Figure 8C). The above results indicate that NSRP1 serves as the alternative splicing regulatory factor, not only modulating the occurrence of alternative splicing in CCDC18 but also concurrently regulating the expression of CCDC18.
Therefore, we further explored the regulatory role of CCDC18 in pre-adipocyte proliferation and differentiation. We synthesized three siRNAs targeting CCDC18 and individually transfected them into pre-adipocytes. The interference efficiency of siRNAs was detected by RT-qPCR, and the results showed that the expression of CCDC18 was significantly inhibited by all three siRNAs (Figure 8D). Among these siRNAs, si1 and si2 were selected to verify the function of CCDC18 in cell proliferation and differentiation. RT-qPCR showed that the expression of proliferative marker genes (Ki67 and PCNA) increased, and the expression of differentiation marker genes (PPARG and CEBPA) decreased (Figure 8E). The EdU assay results showed an increase in the percentage of EdU-positive cells (Figure 8F,G). Oil red O staining revealed a significantly reduced in lipid droplet formation (Figure 8H). In conclusion, CCDC18 inhibits pre-adipocyte proliferation and promotes differentiation, which is consistent with the function of NSRP1.

3. Discussion

As a classical type of gene expression regulatory mechanism, AS plays a crucial role in trait formation [9,25,26,27]. However, there have been few studies on the regulation of phenotypic differences between Eastern and Western pigs by gene alternative splicing. In this study, we utilized transcriptomic data from multiple tissues of Eastern and Western pigs, combined with information on TFs, selection signals, SFs, and QTLs, to comprehensively analyze the differences in AS and its regulatory mechanisms in various tissues of Eastern and Western pigs. Furthermore, we identified marker genes that influence traits such as pig adipose deposition and meat quality.
The study found that the most abundant types of ASEs in all tissues were MEX and SE, consistent with other studies findings [28,29]. The common ASGs in the eight tissues accounted for 63.21%-76.13% of each tissue, indicating a certain degree of conservation in the genes undergoing ASEs across tissues. Additionally, we observed tissue-specificity in both DASGs and DEGs, in line with previous studies [30,31,32]. Among them, adipose and blood tissues showed the highest tissue-specificity in DASGs and DEGs, and these two tissues exhibited the highest average number of ASEs per gene. There are many immune cells in the blood, which are related to the disease resistance of pigs. Disease resistance and fat content are two traits that differ greatly between Eastern and Western pigs [6,7,8]. The above results indicated that AS might ultimately lead to the differences between the two traits of Eastern and Western pigs through its effects on these two tissues.
The coDASG_DEGs in adipose tissue were mainly enriched in pathways related to adipose deposition, immune inflammation, and insulin resistance. Previous researches have shown that pathways such as the Ras signaling pathway, FoxO signaling pathway, PI3K-Akt signaling pathway, insulin resistance, MAPK signaling pathway, and AMPK signaling pathway play essential regulatory roles in adipose deposition [33,34,35,36,37]. Inflammation response and the immune system are closely related to adipose deposition, as they can influence the function, metabolism, and differentiation of adipocytes, thereby affecting adipose deposition [38,39,40]. Studies have demonstrated that inflammation response and insulin resistance downregulation can promote adipose deposition [39,41]. In this study, the downregulation of immune inflammation-related pathways and insulin resistance pathways in Eastern pigs may also promote adipose deposition. These pathways and related genes may be vital in regulating adipose deposition in Eastern pigs. However, further research is needed to explore how these pathways and genes interact and jointly regulate pig adipose deposition.
This study found that TFs played a significant regulatory role in gene AS. TFs are a class of proteins that regulate gene transcription by binding to specific sequences on DNA to initiate or suppress the transcription process [42,43]. They can also regulate splicing processes through interactions with splicing regulatory factors [44,45] or RNA-binding proteins [45]. This study revealed that approximately half of the coDASG_DEGs are significantly influenced by TFs, indicating the important regulatory role of TFs in gene AS. Notably, we identified a TF, ARID4A, in adipose tissue. This TF is under selection in Western pigs but shows high expression in Eastern pigs. Studies have shown that the ARID4A protein can form complexes with proteins such as SIN3A and HDAC1, participating in histone deacetylation modifications, regulating chromatin structure and accessibility, and thereby controlling gene expression and cellular functions [46,47]. ARID4A can influence cell proliferation and division processes by regulating the expression of cell cycle-related genes [48,49]. Moreover, ARID4A is considered a potential tumor suppressor gene [50,51]. However, research on the role of ARID4A in pig adipose deposition is limited.
In this study, we found a strong positive correlation (r > 0.9, P < 3.40E-14) between ARID4A and genes such as AHI1, AKAP9, ANKRD12, CCDC18, ENSSSCG00000037142, IFT74, KIAA2026, NEXN, PPIG, ROCK1, SMC4, TWISTNB, BAZ1A, EIF5B, NSRP1, and BBX. Among these genes, AHI1, AKAP9, ANKRD12, CCDC18, IFT74, NEXN, etc. have been shown to play an important role in adipose deposition or obesity [52,53,54,55,56,57]. Moreover, AHI1, AKAP9, ANKRD12, CCDC18, IFT74, KIAA2026, NEXN, PPIG, ROCK1, SMC4, TWISTNB, EIF5B, and NSRP1 were also identified as potentially related to adipose deposition in this study. ARID4A is also associated with many QTLs related to adipose deposition traits, suggesting that this transcription factor may have an important regulatory role in pig adipose deposition, which requires further investigation.
The phenotypic differences between Western and Eastern pigs are influenced by artificial selection during domestication [1,4,5]. In this study, we further identified 20 candidate SGs in both adipose and blood tissues, which intersect with coDASG_DEGs, indicating that these genes may have been driven by artificial selection during domestication and breeding. Specifically, in adipose tissue, we found that the SOX5, ENSSSCG00000014414, PPP1R12A, ZNF148, TNFSF10, and ROCK1 genes were under selection in Eastern pigs and significantly upregulated in adipose tissue, suggesting their potential influence on adipose deposition. Moreover, previous research has shown that SOX5, PPP1R12A, ZNF148, TNFSF10, and ROCK1 genes regulate adipose deposition processes [58,59,60,61,62,63]. In blood, the MACF1, IL1RL1, MGA, TTC14, ANKRD17, PDLIM7, FGD3, NPRL3, ENSSSCG00000014060, ENSSSCG00000039214, CYB5A, and HAGH genes were under selection in Eastern pigs, and the MACF1, IL1RL1, MGA, TTC14, and ANKRD17 genes were significantly upregulated in Eastern pig blood. Previous reports have implicated these genes in affecting animal immunity and inflammation [64,65,66,67,68,69,70]. They might also play a crucial role in influencing the disease resistance of pigs. These findings suggested that the AS of these candidate SGs might have been subjected to artificial selection during pig domestication, and they could contribute to the phenotypic differences observed between Eastern and Western pigs in adipose deposition and disease resistance.
SFs are a class of proteins that play a crucial role in the process of AS, acting as regulators and mediators [24]. This study found that the SF NSRP1 played an important role in promoting lipid droplet formation during the adipogenic stage. NSRP1 is a protein that regulates the formation and function of spliceosomes [71]. The spliceosome is a complex involved in RNA splicing regulation, and its function is closely related to RNA splicing and post-transcriptional regulation [72]. Therefore, NSRP1 may influence the splicing regulation of genes related to lipid metabolism, thereby affecting adipose deposition. In this study, we discovered that the splicing factor NSRP1, along with 17 highly correlated genes (r > 0.86, P < 6.104E-12), is located in QTLs associated with pig adipose deposition. Among them, ANKRD12, ARID4A, IFT74, KIAA2026, and ROCK1 are located in multiple fat-related QTLs, and ANKRD12, CCDC18, NEXN, PPIG, and ROCK1, are present in the intramuscular fat content QTL. Reports have already showed that ANKRD12, NEXN, IFT74 and ROCK1 play an important role in adipose deposition or obesity [56,73,74,75]. While CCDC18 is a new marker gene for fatty liver in chicken [76], its role in fat deposition has been rarely reported. The present study showed that CCDC18 was highly positively correlated with ARID4A (SG, TF, r = 0.94, P < 1.00E-16) and NSRP1 (SF, r = 0.96, P < 1.00E-16), implying that CCDC18 might be regulated by transcription factor ARID4A and splicing factor NSRP1. Furthermore, this study verified for the first time through a series of experiments that NSRP1 promotes adipogenesis by regulating alternative splicing and expression of CCDC18. However, the role of the transcription factor ARID4A in this process needs further study.

4. Materials and Methods

4.1. Data collection

Transcriptome data of 243 samples from eight tissues of Eastern (Est) and Western (Wst) pigs were obtained from the NCBI-SRA database. The samples comprised adipose (37: Est = 11, Wst = 26), blood (59: Est = 10, Wst = 49), heart (18: Est = 9, Wst = 9), kidney (11: Est = 6, Wst = 5), lung (9: Est = 4, Wst = 5), muscle (55: Est = 34, Wst = 21), ovary (24: Est = 8, Wst = 16), and spleen (30: Est = 25, Wst = 5). For detailed information about the data, please refer to Table S2.

4.2. Transcriptome data quality control and alignment

Firstly, the collected sequencing data underwent quality control using the fastp software to eliminate adapter sequences and low-quality reads [77]. Subsequently, the fastqc software was employed to assess the quality of the clean data obtained after fastp quality control, ensuring that the data met the required quality standards before proceeding with the subsequent alignment [78]. Next, the Sscrofa 11.1 version of the reference genome and annotation files for pigs were downloaded from the Ensembl website (https://ftp.ensembl.org/pub/release-97/), and an index was constructed using STAR [79]. Finally, the STAR software was used to align the quality-assured clean data generated from the fastp quality control process.

4.3. Identification and differential analysis of alternative splicing events

The output bam files from STAR were processed using rMATS to analyze alternative splicing events (ASEs) [80]. We investigated five types of AS, namely alternative 5' splice site (A5SS), alternative 3' splice site (A3SS), mutually exclusive exon (MXE), skipped exon (SE), and retained intron (RI). By applying a cutoff of FDR < 0.05 and |IncLevelDifference| > 5%, we detected differentially alternative splicing events (DASEs). Subsequently, we referred to the genes identified from differentially alternative splicing as Differentially Alternative Splicing Genes (DASGs).

4.4. Identification of differentially expressed genes

First, we utilized featureCounts (v2.0.0) to compute gene expression levels and removed genes with an average count of less than one across all samples [81]. Subsequently, we employed the DESeq2 (v1.38.3) software to normalize gene expression levels and conducted further analysis to identify differentially expressed genes (DEGs) [82]. For this analysis, we applied the criteria |log2(fold change)| > 1 (|log2FC| > 1) and P < 0.05 to select DEGs.

4.5. Gene functional enrichment analysis

We then identified the shared genes between DASGs and DEGs, referred to as coDASG_DEGs. These genes may be subject to regulation by AS and display significant changes in expression. Subsequently, we conducted KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway enrichment analysis on coDASG_DEGs using the R package clusterProfiler (v4.6.2) to explore their biological functions [83]. Pathways with P < 0.05 were considered significantly enriched for these genes.

4.6. Protein-protein interaction network analysis

We analyzed the Protein-Protein Interaction (PPI) of the coDASG_DEGs using the STRING database (https://string-db.org). Employing K-means clustering, we divided the coDASG_DEGs from each tissue into different clusters. Subsequently, we conducted functional enrichment analysis on genes within each cluster using clusterProfiler (v4.6.2) to study their biological significance.

4.7. Integrated analysis of transcription factors, selection signals, splicing factors, and QTLs

We obtained a list of pig transcription factors (TFs) from AnimalTFDB v4.0 (http://bioinfo.life.hust.edu.cn/AnimalTFDB4/#/) and filtered for TFs among the DEGs in each tissue. Subsequently, we conducted a correlation analysis between the identified TFs and coDASG_DEGs to investigate their regulatory role in significant splicing differences between Eastern and Western pigs. We also retrieved the list of genes under selective sweep (SGs) from our previous publication for Eastern and Western pigs [23]. By intersecting this list with coDASG_DEGs, we studied potential coDASG_DEGs influenced by artificial selection during domestication. Splicing factors (SFs) play a regulatory role in the AS process, we further screened for SFs among the DEGs. Then, we performed a correlation analysis between these identified splicing factors and coDASG_DEGs to study their regulatory role in significant splicing differences between Eastern and Western pigs. Lastly, we downloaded pig QTLs information from PigQTLdb (https://www.animalgenome.org/cgi-bin/QTLdb/SS/index) and conducted functional analysis on the coDASG_DEGs highly correlated with SFs. All correlation analyses were performed using the R programming language through Pearson correlation analysis, and a P-value less than 0.05 was considered statistically significant.

4.8. Adipocyte culture and induced differentiation

The stromal vascular fraction cells (SVF cells) from Bama pig were generously provided by Dr. Yangli Pei at Foshan University. The cells were cultured in DMEM supplemented with 10% fetal bovine serum (FBS) and 1% penicillin/streptomycin (Gibco, Grand Island, NY, USA). Upon reaching full confluence and allowing fusion for a day or two, the growth medium was replaced with induction medium to initiate cell differentiation. The cells were cultured until the sixth day, at which point the maintenance medium was changed, and the culture was continued until the lipid droplets reached maturity, with the fluid being refreshed every two days. The induction medium consisted of 1 mg/ml insulin, 0.5 mmol/L isobutyl-1-methylxanthine (IBMX), 1 mmol/L dexamethasone, 1 μmol/L rosiglitazone, and 10% FBS. The maintenance medium included 1 mg/ml insulin and 10% FBS. Transient transfections were carried out using Lipofectamine 3000 (Invitrogen, USA), following the manufacturer′s instructions.

4.9. RNA extraction and RT-qPCR

Total RNA was extracted using Trizol (Invitrogen, Shanghai, China) following the manufacturer′s instructions. The quality and quantity of RNA were assessed using the NanoDrop 2000 (Thermo Fisher Scientific, Massachusetts, USA). Quantitative real-time PCR (RT-qPCR) was conducted using Fast ChamQ Universal SYBR qPCR Master Mix (Vazyme, Nanjing, China) for mRNA on an ABI Step One Plus Real-Time PCR system (Applied Biosystems, USA). The 2-ΔΔCt method was employed to analyze the relative expression levels of mRNA. β-actin was used as endogenous controls to normalize the expression of mRNA. The sequence information of primers used for RT-qPCR is shown in Table S3.

4.10. Western blot analysis

Proteins were extracted by RIPA buffer (Thermo Scientific, Massachusetts, USA) supplemented with phosphorylase inhibitor (Roche 5892791001, Basel, Switzerland) and protease inhibitor (Roche 04693132001, Basel, Switzerland). The concentration of obtained protein was measured by the BCA kit (Beyotime, China). The proteins were separated in 10% sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) gels (EpiZyme, Shanghai, China) and transferred onto 0.45 µm Hybridization Nitrocellulose Filter (NC) membrane (Merck, New Jersey, USA), and then probed with antibodies following standard procedures. These membranes were blocked with 5% skim milk at room temperature for 2 h, and then incubated with primary antibodies at 4°C overnight. Subsequently, the membranes were incubated with secondary antibodies at room temperature for 1 h. The following antibodies were used in the present work: NSRP1 (1:500, 21360-1-AP, Proteintech, China), GAPDH (1:50000, 60004-1-Ig, Proteintech, China), β-actin (1:20000, 66009-1-Ig, Proteintech, China), KI67 (ab16667; 1:1000; Abcam, UK), PCNA (1:5000, 10205-2-AP, Proteintech Group, China), PPARG (1:1000, 16643-1-AP, Proteintech Group, China), CEBPA (1:500, 18311-1-AP, Proteintech Group, China). Secondary antibodies: Goat Anti-Rabbit (ZB-2301, 1:1000, ZSGB-BIO, China) and Goat Anti-Mouse (ZB-2305, 1:1000, ZSGB-BIO, China).

4.11. Cell Counting Kit-8 proliferation assay

SVF cells were seeded into 96-well plates, and the proliferation of cells was assessed at 0 h, 24 h, 48 h and 72 h after transfection using the Cell Counting Kit-8 (CCK-8) (Beyotime C0038, Beijing, China). After incubation for 1 h, the absorbance at 450 nm was measured using a microplate reader.

4.12. 5-Ethynyl-2′-deoxyuridine (EdU) staining

Cells were seeded in 12-well plates and cultured until they reached 50% confluency. After that, they were transfected and allowed to incubate for 48 hours. Following the 48-hour transfection period, EdU staining was performed using the BeyoClickTM EdU Cell Proliferation Kit (Beyotime, China). Briefly, the cells were fixed in 4% paraformaldehyde and permeabilized with 0.5% Triton X-100. Next, the cells were stained with Click Additive Solution in the dark for 30 minutes, and the nuclei were counterstained with 4’,6-diamidino-2-phenylindole (DAPI) solution. Images were acquired using a Nikon ECLIPSE Ti microscope, and the ImageJ software was employed to calculate the proportion of EdU-positive cells.

4.13. Oil Red O Staining

After induction of adipogenesis, cells were washed three times with PBS (Gibco, Carlsbad, CA, USA) and then fixed for 30 min in 4% paraformaldehyde. The samples were rinsed twice with 60% isopropanol and dried for 30 min before being treated with 1 mL of the oil red O dye working solution. A microscope was used to observe the oil red O staining after adding 1 mL PBS (Gibco, Carlsbad, CA, USA) to the culture plate.

4.14. Semiquantitative RT-PCR analysis of alternative splicing events

PCR products were separated by 1.5% agarose gel in 1×TAE buffer for 40 min at 120 V. Quantification of gels was performed by densitometry using ImageJ software (National Institutes of Health) for analysis. The sequences of the primers (Sangon Biotech) were as Table S4.

4.15. Statistical analysis

Alternative splicing events were automatically detected and quantified using the percent-spliced-in (PSI, C) metric based on long (L) and short (S) forms of splicing events presents (equation shown below). Briefly, a PSI value was given according to the ratio of the long form on total form present (short form and long form) to characterize inclusion of exon.
C=L/(L+S)
The results are represented as the mean ± SD. Statistical analyses of the differences between groups were performed using Student’s t-test. Statistical significance was set at *P < 0.05, **P < 0.01 and ***P < 0.001.

5. Conclusions

In the present study, we found that AS might play an important role in fat deposition and disease resistance in pigs, which might be affected by transcription factors, selection signals and splicing factors. We identified the TF ARID4A and the SF NSRP1 in adipose tissue, along with genes such as ANKRD12, IFT74, KIAA2026, CCDC18, NEXN, PPIG, and ROCK1 might play critical roles in adipose-deposition traits. For the first time, we further confirmed that NSRP1 regulates adipogenesis by regulating alternative splicing and expression of CCDC18 through experimental verification. These findings contribute to a better understanding of the formation and regulatory mechanisms of meat quality traits. The identified marker genes can serve as molecular markers for pig genetic breeding and meat quality improvement practices, providing a scientific basis for enhancing and optimizing the pork industry.

Supplementary Materials

The following supporting information can be downloaded at the website of this paper posted on Preprints.org, Table S1: QTLs information that were associated with adipose-deposition traits; Table S2: The detailed information of the transcriptome data; Table S3: Primers used for RT-qPCR; Table S4: Primers used for semiquantitative RT-PCR.

Author Contributions

Z.T. conceived the study. L.L. and X.L. performed bioinformatics analyses. W.W. and W.L. were responsible for experimental verification work. L.L. and W.W. wrote the original draft. Z.T., G.Y., A.A.A. and H.H. revised the manuscript. All authors reviewed and approved the final manuscript.

Funding

This work was supported by the National Natural Science Foundation of China (32002151, U23A20229), Science and Technology Major Project of Shenzhen (KJZD20230923115003006), sustainable development special project from Shenzhen (KCXFZ20201221173213037) and Science and Technology Program project Yazhou Bay Seed Laboratory (B21HJ0211).

Institutional Review Board Statement

The animal procedures were conducted in adherence to the guidelines established by the Institutional Animal Care and Use Committee of the Agricultural Genomics Institute at Shenzhen, Chinese Academy of Agricultural Sciences.

Informed Consent Statement

Not applicable.

Data Availability Statement

The RNA-seq data was obtained from the NCBI-SRA database (https://www.ncbi.nlm.nih.gov/sra/?term=). The detailed information of the transcriptome data has been showed in Table S2.

Acknowledgments

We are very grateful to Dr. Yangli Pei from Foshan University for generously providing stromal vascular fraction cells (SVF cells) of Bama pig.

Conflicts of Interest

The authors declare that there are no conflicts of interest to disclose.

References

  1. White, S. , From globalized pig breeds to capitalist pigs: a study in animal cultures and evolutionary history. Environmental History 2011. [Google Scholar] [CrossRef]
  2. Zhang, L.; Tian, H.; Shi, H.; Pan, S.; Qin, X.; Pan, N.; Dangal, S. R. , Methane emissions from livestock in East Asia during 1961− 2019. Ecosystem Health and Sustainability 2021, 7, 1918024. [Google Scholar] [CrossRef]
  3. Choi, S. K.; Lee, J.-E.; Kim, Y.-J.; Min, M.-S.; Voloshina, I.; Myslenkov, A.; Oh, J. G.; Kim, T.-H.; Markov, N.; Seryodkin, I. , Genetic structure of wild boar (Sus scrofa) populations from East Asia based on microsatellite loci analyses. BMC genetics 2014, 15, 1–10. [Google Scholar] [CrossRef] [PubMed]
  4. Larson, G.; Dobney, K.; Albarella, U.; Fang, M.; Matisoo-Smith, E.; Robins, J.; Lowden, S.; Finlayson, H.; Brand, T.; Willerslev, E. , Worldwide phylogeography of wild boar reveals multiple centers of pig domestication. Science 2005, 307, 1618–1621. [Google Scholar] [CrossRef] [PubMed]
  5. Evin, A.; Cucchi, T.; Cardini, A.; Vidarsdottir, U. S.; Larson, G.; Dobney, K. , The long and winding road: identifying pig domestication through molar size and shape. Journal of Archaeological Science 2013, 40, 735–743. [Google Scholar] [CrossRef]
  6. Liu, Y.; Fu, Y.; Yang, Y.; Yi, G.; Lian, J.; Xie, B.; Yao, Y.; Chen, M.; Niu, Y.; Liu, L.; Wang, L.; Zhang, Y.; Fan, X.; Tang, Y.; Yuan, P.; Zhu, M.; Li, Q.; Zhang, S.; Chen, Y.; Wang, B.; He, J.; Lu, D.; Liachko, I.; Sullivan, S. T.; Pang, B.; Chen, Y.; He, X.; Li, K.; Tang, Z. , Integration of multi-omics data reveals cis-regulatory variants that are associated with phenotypic differentiation of eastern from western pigs. Genet Sel Evol 2022, 54, 62. [Google Scholar] [CrossRef]
  7. Deng, L.; Li, W.; Liu, W.; Liu, Y.; Xie, B.; Groenen, M. A. M.; Madsen, O.; Yang, X.; Tang, Z. , Integrative metabolomic and transcriptomic analysis reveals difference in glucose and lipid metabolism in the longissimus muscle of Luchuan and Duroc pigs. Front Genet 2023, 14, 1128033. [Google Scholar] [CrossRef]
  8. Yang, Y.; Liang, G.; Niu, G.; Zhang, Y.; Zhou, R.; Wang, Y.; Mu, Y.; Tang, Z.; Li, K. , Comparative analysis of DNA methylome and transcriptome of skeletal muscle in lean-, obese-, and mini-type pigs. Sci Rep 2017, 7, 39883. [Google Scholar] [CrossRef]
  9. Wright, C. J.; Smith, C. W.; Jiggins, C. D. , Alternative splicing as a source of phenotypic diversity. Nature Reviews Genetics 2022, 23, 697–710. [Google Scholar] [CrossRef]
  10. Chen, H.; Gao, F.; He, M.; Ding, X. F.; Wong, A. M.; Sze, S. C.; Yu, A. C.; Sun, T.; Chan, A. W. H.; Wang, X. , Long-read RNA sequencing identifies alternative splice variants in hepatocellular carcinoma and tumor-specific isoforms. Hepatology 2019, 70, 1011–1025. [Google Scholar] [CrossRef]
  11. Baralle, F. E.; Giudice, J. , Alternative splicing as a regulator of development and tissue identity. Nature reviews Molecular cell biology 2017, 18, 437–451. [Google Scholar] [CrossRef]
  12. Mazin, P. V.; Khaitovich, P.; Cardoso-Moreira, M.; Kaessmann, H. , Alternative splicing during mammalian organ development. Nature genetics 2021, 53, 925–934. [Google Scholar] [CrossRef] [PubMed]
  13. Modrek, B.; Lee, C. , A genomic view of alternative splicing. Nature genetics 2002, 30, 13–19. [Google Scholar] [CrossRef] [PubMed]
  14. Scotti, M. M.; Swanson, M. S. , RNA mis-splicing in disease. Nature Reviews Genetics 2016, 17, 19–32. [Google Scholar] [CrossRef] [PubMed]
  15. Han, H.; Braunschweig, U.; Gonatopoulos-Pournatzis, T.; Weatheritt, R. J.; Hirsch, C. L.; Ha, K. C. H.; Radovani, E.; Nabeel-Shah, S.; Sterne-Weiler, T.; Wang, J.; O'Hanlon, D.; Pan, Q.; Ray, D.; Zheng, H.; Vizeacoumar, F.; Datti, A.; Magomedova, L.; Cummins, C. L.; Hughes, T. R.; Greenblatt, J. F.; Wrana, J. L.; Moffat, J.; Blencowe, B. J. , Multilayered Control of Alternative Splicing Regulatory Networks by Transcription Factors. Mol Cell 2017, 65, 539–553 e7. [Google Scholar] [CrossRef] [PubMed]
  16. Bielli, P.; Busa, R.; Di Stasi, S. M.; Munoz, M. J.; Botti, F.; Kornblihtt, A. R.; Sette, C. , The transcription factor FBI-1 inhibits SAM68-mediated BCL-X alternative splicing and apoptosis. EMBO Rep 2014, 15, 419–27. [Google Scholar] [CrossRef] [PubMed]
  17. Long, J. C.; Caceres, J. F. , The SR protein family of splicing factors: master regulators of gene expression. Biochem J 2009, 417, 15–27. [Google Scholar] [CrossRef] [PubMed]
  18. Chen, K.; Baxter, T.; Muir, W. M.; Groenen, M. A.; Schook, L. B. , Genetic resources, genome mapping and evolutionary genomics of the pig (Sus scrofa). Int J Biol Sci 2007, 3, 153–65. [Google Scholar] [CrossRef] [PubMed]
  19. Li, P. J. , Exponential growth, animal welfare, environmental and food safety impact: The case of China’s livestock production. Journal of Agricultural and Environmental Ethics 2009, 22, 217–240. [Google Scholar] [CrossRef]
  20. Nygard, A.-B.; Cirera, S.; Gilchrist, M. J.; Gorodkin, J.; Jørgensen, C. B.; Fredholm, M. , A study of alternative splicing in the pig. BMC Research Notes 2010, 3, 1–12. [Google Scholar] [CrossRef]
  21. Hao, W.; Yang, Z.; Sun, Y.; Li, J.; Zhang, D.; Liu, D.; Yang, X. , Characterization of alternative splicing events in porcine skeletal muscles with different intramuscular Fat contents. Biomolecules 2022, 12, 154. [Google Scholar] [CrossRef] [PubMed]
  22. Shu, Z.; Wang, L.; Wang, J.; Zhang, L.; Hou, X.; Yan, H.; Wang, L. , Integrative Analysis of Nanopore and Illumina Sequencing Reveals Alternative Splicing Complexity in Pig Longissimus Dorsi Muscle. Frontiers in Genetics 2022, 13, 877646. [Google Scholar] [CrossRef] [PubMed]
  23. Liu, Y.; Fu, Y.; Yang, Y.; Yi, G.; Lian, J.; Xie, B.; Yao, Y.; Chen, M.; Niu, Y.; Liu, L. , Integration of multi-omics data reveals cis-regulatory variants that are associated with phenotypic differentiation of eastern from western pigs. Genetics Selection Evolution 2022, 54, 62. [Google Scholar] [CrossRef]
  24. Long, J. C.; Caceres, J. F. , The SR protein family of splicing factors: master regulators of gene expression. Biochemical Journal 2009, 417, 15–27. [Google Scholar] [CrossRef] [PubMed]
  25. Bhadra, M.; Howell, P.; Dutta, S.; Heintz, C.; Mair, W. B. , Alternative splicing in aging and longevity. Human genetics 2020, 139, 357–369. [Google Scholar] [CrossRef]
  26. Jacobs, A.; Elmer, K. R. , Alternative splicing and gene expression play contrasting roles in the parallel phenotypic evolution of a salmonid fish. Molecular Ecology 2021, 30, 4955–4969. [Google Scholar] [CrossRef]
  27. Bush, S. J.; Chen, L.; Tovar-Corona, J. M.; Urrutia, A. O. , Alternative splicing and the evolution of phenotypic novelty. Philosophical Transactions of the Royal Society B: Biological Sciences 2017, 372, 20150474. [Google Scholar] [CrossRef]
  28. Yang, Y.; Yuan, H.; Yang, Q.; Cai, Y.; Ren, Y.; Li, Y.; Gao, C.; Zhao, S. , Post-transcriptional regulation through alternative splicing in the lungs of Tibetan pigs under hypoxia. Gene 2022, 819, 146268. [Google Scholar] [CrossRef]
  29. Tan, B.; Zeng, J.; Meng, F.; Wang, S.; Xiao, L.; Zhao, X.; Hong, L.; Zheng, E.; Wu, Z.; Li, Z. , Comprehensive analysis of pre-mRNA alternative splicing regulated by m6A methylation in pig oxidative and glycolytic skeletal muscles. BMC genomics 2022, 23, 1–14. [Google Scholar] [CrossRef]
  30. Jiang, W.; Chen, L. , Tissue specificity of gene expression evolves across mammal species. Journal of Computational Biology 2022, 29, 880–891. [Google Scholar] [CrossRef]
  31. Mosthaf, L.; Grako, K.; Dull, T.; Coussens, L.; Ullrich, A.; McClain, D. , Functionally distinct insulin receptors generated by tissue-specific alternative splicing. The EMBO journal 1990, 9, 2409–2413. [Google Scholar] [CrossRef] [PubMed]
  32. Qu, A.; Bai, Y.; Zhang, X.; Zeng, J.; Pu, F.; Wu, L.; Xu, P.; Zhou, T. , Tissue-Specific Analysis of Alternative Splicing Events and Differential Isoform Expression in Large Yellow Croaker (Larimichthys crocea) After Cryptocaryon irritans Infection. Marine Biotechnology 2022, 24, 640–654. [Google Scholar] [CrossRef] [PubMed]
  33. Yu, T.; Tian, X.; Li, D.; He, Y.; Yang, P.; Cheng, Y.; Zhao, X.; Sun, J.; Yang, G. , Transcriptome, proteome and metabolome analysis provide insights on fat deposition and meat quality in pig. Food Research International 2023, 166, 112550. [Google Scholar] [CrossRef] [PubMed]
  34. Shen, J.; Hao, Z.; Luo, Y.; Zhen, H.; Liu, Y.; Wang, J.; Hu, J.; Liu, X.; Li, S.; Zhao, Z. , Deep small RNA sequencing reveals important miRNAs related to muscle development and intramuscular fat deposition in longissimus dorsi muscle from different goat breeds. Frontiers in Veterinary Science 2022, 9, 911166. [Google Scholar] [CrossRef] [PubMed]
  35. Wang, H.; Wang, X.; Li, M.; Wang, S.; Chen, Q.; Lu, S. , Identification of key sex-specific pathways and genes in the subcutaneous adipose tissue from pigs using WGCNA method. BMC Genomic Data 2022, 23, 35. [Google Scholar] [CrossRef] [PubMed]
  36. Jiao, Y.; Liang, X.; Hou, J.; Aisa, Y.; Wu, H.; Zhang, Z.; Nuermaimaiti, N.; Zhao, Y.; Jiang, S.; Guan, Y. , Adenovirus type 36 regulates adipose stem cell differentiation and glucolipid metabolism through the PI3K/Akt/FoxO1/PPARγ signaling pathway. Lipids in Health and Disease 2019, 18, 1–12. [Google Scholar] [CrossRef] [PubMed]
  37. Chen, H.; Peng, T.; Shang, H.; Shang, X.; Zhao, X.; Qu, M.; Song, X. , RNA-Seq analysis reveals the potential molecular mechanisms of puerarin on intramuscular fat deposition in heat-stressed beef cattle. Frontiers in Nutrition 2022, 9, 817557. [Google Scholar] [CrossRef]
  38. Poklukar, K.; Čandek-Potokar, M.; Batorek Lukač, N.; Tomažin, U.; Škrlep, M. , Lipid deposition and metabolism in local and modern pig breeds: A review. Animals 2020, 10, 424. [Google Scholar] [CrossRef]
  39. Ahmed, B.; Sultana, R.; Greene, M. W. , Adipose tissue and insulin resistance in obese. Biomedicine & Pharmacotherapy 2021, 137, 111315. [Google Scholar]
  40. Nedunchezhiyan, U.; Varughese, I.; Sun, A. R.; Wu, X.; Crawford, R.; Prasadam, I. , Obesity, inflammation, and immune system in osteoarthritis. Front Immunol 2022, 13, 907750. [Google Scholar] [CrossRef]
  41. Zhao, G.-N.; Tian, Z.-W.; Tian, T.; Zhu, Z.-P.; Zhao, W.-J.; Tian, H.; Cheng, X.; Hu, F.-J.; Hu, M.-L.; Tian, S. , TMBIM1 is an inhibitor of adipogenesis and its depletion promotes adipocyte hyperplasia and improves obesity-related metabolic disease. Cell metabolism 2021, 33, 1640–1654. e8. [Google Scholar] [CrossRef] [PubMed]
  42. Stadhouders, R.; Filion, G. J.; Graf, T. , Transcription factors and 3D genome conformation in cell-fate decisions. Nature 2019, 569, 345–354. [Google Scholar] [CrossRef] [PubMed]
  43. Andrews, G.; Fan, K.; Pratt, H. E.; Phalke, N.; §, Z. C.; Karlsson, E. K.; Lindblad-Toh, K.; Gazal, S.; Moore, J. E.; Weng, Z. , Mammalian evolution of human cis-regulatory elements and transcription factor binding sites. Science 2023, 380, eabn7930. [Google Scholar] [CrossRef] [PubMed]
  44. Bielli, P.; Busà, R.; Di Stasi, S. M.; Munoz, M. J.; Botti, F.; Kornblihtt, A. R.; Sette, C. , The transcription factor FBI-1 inhibits SAM 68-mediated BCL-X alternative splicing and apoptosis. EMBO reports 2014, 15, 419–427. [Google Scholar] [CrossRef] [PubMed]
  45. Han, H.; Braunschweig, U.; Gonatopoulos-Pournatzis, T.; Weatheritt, R. J.; Hirsch, C. L.; Ha, K. C.; Radovani, E.; Nabeel-Shah, S.; Sterne-Weiler, T.; Wang, J. , Multilayered control of alternative splicing regulatory networks by transcription factors. Molecular cell 2017, 65, 539–553. e7. [Google Scholar] [CrossRef] [PubMed]
  46. Farooqi, A. A.; Fayyaz, S.; Poltronieri, P.; Calin, G.; Mallardo, M. In Epigenetic deregulation in cancer: Enzyme players and non-coding RNAs, Seminars in cancer biology, 2022; Elsevier: 2022; pp 197-207.
  47. Milazzo, G.; Mercatelli, D.; Di Muzio, G.; Triboli, L.; De Rosa, P.; Perini, G.; Giorgi, F. M. , Histone deacetylases (HDACs): evolution, specificity, role in transcriptional complexes, and pharmacological actionability. Genes 2020, 11, 556. [Google Scholar] [CrossRef] [PubMed]
  48. Delcuve, G. P.; Khan, D. H.; Davie, J. R. , Targeting class I histone deacetylases in cancer therapy. Expert opinion on therapeutic targets 2013, 17, 29–41. [Google Scholar] [CrossRef] [PubMed]
  49. rivas, P. D.; Papavassiliou, A. G., Transcriptional corepressors in cancer: emerging targets for therapeutic intervention. In Wiley Online Library: 2013; Vol. 119, pp 1120-1128.
  50. Wu, M.-Y.; Eldin, K. W.; Beaudet, A. L. , Identification of chromatin remodeling genes Arid4a and Arid4b as leukemia suppressor genes. Journal of the National Cancer Institute 2008, 100, 1247–1259. [Google Scholar] [CrossRef]
  51. Liang, Y. K.; Han, Z. D.; Lu, J. M.; Liu, Z. Z.; Zhuo, Y. J.; Zhu, X. J.; Chen, J. X.; Ye, J. H.; Liang, Y. X.; He, H. C. , Downregulation of ARID4A and ARID4B promote tumor progression and directly regulated by microRNA-30d in patient with prostate cancer. Journal of Cellular Biochemistry 2018, 119, 7245–7255. [Google Scholar] [CrossRef]
  52. Huang, L.; Teng, D.; Wang, H.; Sheng, G.; Liu, T. , Association of copy number variation in the AHI1 gene with risk of obesity in the Chinese population. European journal of endocrinology 2012, 166, 727–734. [Google Scholar] [CrossRef]
  53. Cho, Y.; Gutierrez, L.; Bordonaro, M.; Russo, D.; Anzelmi, F.; Hooven, J. T.; Cerra, C.; Lazarova, D. L. , Effects of propolis and gamma-cyclodextrin on intestinal neoplasia in normal weight and obese mice. Cancer Medicine 2016, 5, 2448–2458. [Google Scholar] [CrossRef] [PubMed]
  54. Deaton, A. M.; Dubey, A.; Ward, L. D.; Dornbos, P.; Flannick, J.; Consortium, A.-T. D.-G.; Yee, E.; Ticau, S.; Noetzli, L.; Parker, M. M. , Rare loss of function variants in the hepatokine gene INHBE protect from abdominal obesity. Nature Communications 2022, 13, 4319. [Google Scholar] [CrossRef] [PubMed]
  55. Tan, X.; Liu, R.; Xing, S.; Zhang, Y.; Li, Q.; Zheng, M.; Zhao, G.; Wen, J. , Genome-wide detection of key genes and epigenetic markers for chicken fatty liver. International Journal of Molecular Sciences 2020, 21, 1800. [Google Scholar] [CrossRef] [PubMed]
  56. Zappaterra, M.; Gioiosa, S.; Chillemi, G.; Zambonelli, P.; Davoli, R. , Muscle transcriptome analysis identifies genes involved in ciliogenesis and the molecular cascade associated with intramuscular fat content in Large White heavy pigs. PLoS One 2020, 15, e0233372. [Google Scholar] [CrossRef] [PubMed]
  57. Pan, C.; Yang, C.; Wang, S.; Ma, Y. , Identifying key genes and functionally enriched pathways of diverse adipose tissue types in cattle. Frontiers in Genetics 2022, 13, 790690. [Google Scholar] [CrossRef] [PubMed]
  58. Peng, Y.; Shi, X.-e.; Huang, K.-l.; Yao, X.-p.; Chen, F.-f.; Li, X.; Yang, G.-s. , Knock-down Sox5 suppresses porcine adipogenesis through BMP R-Smads signal pathway. Biochemical and Biophysical Research Communications 2020, 527, 574–580. [Google Scholar] [CrossRef] [PubMed]
  59. Li, Z.; Liu, X.; Li, Y.; Wang, W.; Wang, N.; Xiao, F.; Gao, H.; Guo, H.; Li, H.; Wang, S. , Chicken C/EBPζ gene: Expression profiles, association analysis, and identification of functional variants for abdominal fat. Domestic Animal Endocrinology 2021, 76, 106631. [Google Scholar] [CrossRef]
  60. Jiang, H.; Gao, Q.; Che, X.; Zhu, L.; Zhang, Z.; Chen, Y.; Dai, Y. , Inhibition of penile tunica albuginea myofibroblasts activity by adipose-derived stem cells. Experimental and Therapeutic Medicine 2017, 14, 5149–5156. [Google Scholar] [CrossRef]
  61. Dou, Y.; Qi, K.; Liu, Y.; Li, C.; Song, C.; Wei, Y.; Zhang, Z.; Li, X.; Wang, K.; Li, X. , Identification and Functional Prediction of Long Non-Coding RNA in Longissimus Dorsi Muscle of Queshan Black and Large White Pigs. Genes 2023, 14, 197. [Google Scholar] [CrossRef]
  62. Mota, L. F.; Lopes, F. B.; Fernandes Júnior, G. A.; Rosa, G. J.; Magalhães, A. F.; Carvalheiro, R.; Albuquerque, L. G. , Genome-wide scan highlights the role of candidate genes on phenotypic plasticity for age at first calving in Nellore heifers. Scientific reports 2020, 10, 6481. [Google Scholar] [CrossRef]
  63. Liu, X.; Wei, S.; Deng, S.; Li, D.; Liu, K.; Shan, B.; Shao, Y.; Wei, W.; Chen, J.; Zhang, L. , Genome-wide identification and comparison of mRNA s, lnc RNA s and circ RNA s in porcine intramuscular, subcutaneous, retroperitoneal and mesenteric adipose tissues. Animal genetics 2019, 50, 228–241. [Google Scholar] [CrossRef]
  64. Hu, L.; Xiao, Y.; Xiong, Z.; Zhao, F.; Yin, C.; Zhang, Y.; Su, P.; Li, D.; Chen, Z.; Ma, X. In MACF1, versatility in tissue-specific function and in human disease, Seminars in cell & developmental biology, 2017; Elsevier: 2017; pp 3-8.
  65. Akhabir, L.; Sandford, A. , Genetics of interleukin 1 receptor-like 1 in immune and inflammatory diseases. Current genomics 2010, 11, 591–606. [Google Scholar] [CrossRef] [PubMed]
  66. Shaw, C. H.; Palti, Y.; Gao, G.; Wiens, G. D. , Evolution of interleukin-1 receptor-like 1 and its role in rainbow trout (Oncorhynchus mykiss) resistance to Flavobacterium psychrophilum. The Journal of Immunology 2017, 198 (1_Supplement), 201.22–201.22. [Google Scholar] [CrossRef]
  67. Chen, Y.; Zeng, B.; Shi, P.; Xiao, H.; Chen, S. , Comparative analysis of the Liver and Spleen transcriptomes between Holstein and Yunnan humped cattle. Animals 2019, 9, 527. [Google Scholar] [CrossRef] [PubMed]
  68. Sun, L.; Li, M.; Deng, L.; Niu, Y.; Tang, Y.; Wang, Y.; Guo, L. , MGA Mutation as a novel biomarker for immune checkpoint therapies in non-squamous non-small cell lung cancer. Frontiers in pharmacology 2021, 12, 625593. [Google Scholar] [CrossRef]
  69. Wang, J.; Shivakumar, S.; Barker, K.; Tang, Y.; Wallstrom, G.; Park, J. G.; Tsay, J.-C. J.; Pass, H. I.; Rom, W. N.; LaBaer, J. , Comparative study of autoantibody responses between lung adenocarcinoma and benign pulmonary nodules. Journal of Thoracic Oncology 2016, 11, 334–345. [Google Scholar] [CrossRef] [PubMed]
  70. Zhou, G.; Jiang, N.; Zhang, W.; Guo, S.; Xin, G. , Biomarker identification in membranous nephropathy using a long non-coding RNA-mediated competitive endogenous RNA network. Interdisciplinary Sciences: Computational Life Sciences 2021, 13, 615–623. [Google Scholar] [CrossRef]
  71. Kim, C.-H.; Park, S.-M.; Lee, S.-j.; Kim, Y.-D.; Jang, S.-H.; Woo, S.-M.; Kwon, T.-K.; Park, Z.-Y.; Chung, I.-J.; Kim, H.-R. , NSrp70 is a lymphocyte-essential splicing factor that controls thymocyte development. Nucleic Acids Research 2021, 49, 5760–5778. [Google Scholar] [CrossRef]
  72. Shi, Y. , Mechanistic insights into precursor messenger RNA splicing by the spliceosome. Nature reviews Molecular cell biology 2017, 18, 655–670. [Google Scholar] [CrossRef]
  73. Jiang, H.; Gao, Q.; Che, X.; Zhu, L.; Zhang, Z.; Chen, Y.; Dai, Y. , Inhibition of penile tunica albuginea myofibroblasts activity by adipose-derived stem cells. Exp Ther Med 2017, 14, 5149–5156. [Google Scholar] [CrossRef]
  74. Deaton, A. M.; Dubey, A.; Ward, L. D.; Dornbos, P.; Flannick, J.; Consortium, A.-T. D. G.; Yee, E.; Ticau, S.; Noetzli, L.; Parker, M. M.; Hoffing, R. A.; Willis, C.; Plekan, M. E.; Holleman, A. M.; Hinkle, G.; Fitzgerald, K.; Vaishnaw, A. K.; Nioi, P. , Rare loss of function variants in the hepatokine gene INHBE protect from abdominal obesity. Nat Commun 2022, 13, 4319. [Google Scholar] [CrossRef] [PubMed]
  75. Pan, C.; Yang, C.; Wang, S.; Ma, Y. , Identifying Key Genes and Functionally Enriched Pathways of Diverse Adipose Tissue Types in Cattle. Front Genet 2022, 13, 790690. [Google Scholar] [CrossRef] [PubMed]
  76. Tan, X.; Liu, R.; Xing, S.; Zhang, Y.; Li, Q.; Zheng, M.; Zhao, G.; Wen, J. , Genome-Wide Detection of Key Genes and Epigenetic Markers for Chicken Fatty Liver. Int J Mol Sci 2020, 21. [Google Scholar] [CrossRef]
  77. Chen, S.; Zhou, Y.; Chen, Y.; Gu, J. , fastp: an ultra-fast all-in-one FASTQ preprocessor. Bioinformatics 2018, 34, i884–i890. [Google Scholar] [CrossRef] [PubMed]
  78. Andrews, S., FastQC: a quality control tool for high throughput sequence data. In Babraham Bioinformatics, Babraham Institute, Cambridge, United Kingdom: 2010.
  79. Dobin, A.; Davis, C. A.; Schlesinger, F.; Drenkow, J.; Zaleski, C.; Jha, S.; Batut, P.; Chaisson, M.; Gingeras, T. R. , STAR: ultrafast universal RNA-seq aligner. Bioinformatics 2013, 29, 15–21. [Google Scholar] [CrossRef] [PubMed]
  80. Shen, S.; Park, J. W.; Lu, Z.-x.; Lin, L.; Henry, M. D.; Wu, Y. N.; Zhou, Q.; Xing, Y. , rMATS: robust and flexible detection of differential alternative splicing from replicate RNA-Seq data. Proceedings of the National Academy of Sciences 2014, 111, E5593–E5601. [Google Scholar] [CrossRef] [PubMed]
  81. Liao, Y.; Smyth, G. K.; Shi, W. , featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics 2014, 30, 923–930. [Google Scholar] [CrossRef]
  82. Love, M. I.; Huber, W.; Anders, S. , Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome biology 2014, 15, 1–21. [Google Scholar] [CrossRef]
  83. Yu, G.; Wang, L.-G.; Han, Y.; He, Q.-Y. , clusterProfiler: an R package for comparing biological themes among gene clusters. Omics: a journal of integrative biology 2012, 16, 284–287. [Google Scholar] [CrossRef]
Preprints 96858 g001
Figure 1. Identification of alternative splicing across tissues. (A) Experimental design and flowchart. DEGs = Differentially Expressed Genes; DASGs = Differentially Alternative Splicing Genes; coDASG_DEGs = Common genes between DASGs and DEGs; PPI = Protein-Protein Interaction; A5SS = Alternative 5’ Splice Site; A3SS = Alternative 3’ Splice Site; MXE = Mutually Exclusive Exon; SE = Skipped Exon; RI = Retained Intron; TFs =Transcription Factors; SGs = Genes under Selection; SFs = Splicing Factors; QTLs = Quantitative Trait Loci. (B) Alternative splicing events and genes in multiple tissues. ASEs = Alternative Splicing Events; DASEs = Differentially Alternative Splicing Events; ASGs = Alternative Splicing Genes; DASGs = Differentially Alternative Splicing Genes; ASGs/% = The proportion of ASGs in all genes of the pig genome; The values before and after the slash of A3SS, A5SS, MXE, RI and SE represent the number of DASEs and ASEs, respectively; The numbers in parentheses of DASEs/ASEs(%) indicate the percentage of DASEs in ASEs; The numbers in parentheses of DASGs/ASGs(%) indicate the percentage of DASGs in ASGs. (C) Conservation of alternative splicing across tissues. The panel above indicates the overlapped ASGs among different tissues. The bottom panel represents statistics on the number of ASGs in different tissues.
Figure 1. Identification of alternative splicing across tissues. (A) Experimental design and flowchart. DEGs = Differentially Expressed Genes; DASGs = Differentially Alternative Splicing Genes; coDASG_DEGs = Common genes between DASGs and DEGs; PPI = Protein-Protein Interaction; A5SS = Alternative 5’ Splice Site; A3SS = Alternative 3’ Splice Site; MXE = Mutually Exclusive Exon; SE = Skipped Exon; RI = Retained Intron; TFs =Transcription Factors; SGs = Genes under Selection; SFs = Splicing Factors; QTLs = Quantitative Trait Loci. (B) Alternative splicing events and genes in multiple tissues. ASEs = Alternative Splicing Events; DASEs = Differentially Alternative Splicing Events; ASGs = Alternative Splicing Genes; DASGs = Differentially Alternative Splicing Genes; ASGs/% = The proportion of ASGs in all genes of the pig genome; The values before and after the slash of A3SS, A5SS, MXE, RI and SE represent the number of DASEs and ASEs, respectively; The numbers in parentheses of DASEs/ASEs(%) indicate the percentage of DASEs in ASEs; The numbers in parentheses of DASGs/ASGs(%) indicate the percentage of DASGs in ASGs. (C) Conservation of alternative splicing across tissues. The panel above indicates the overlapped ASGs among different tissues. The bottom panel represents statistics on the number of ASGs in different tissues.
Preprints 96858 g001
Preprints 96858 g002
Figure 2. Functional enrichment analysis of the coDASG_DEGs. (A) Significant KEGG pathway of coDASG_DEGs in adipose (P < 0.05). (B) PPI analysis of coDASG_DEGs in adipose tissue. (C) Significant KEGG pathway of coDASG_DEGs in blood (P < 0.05). (D) PPI analysis of coDASG_DEGs in blood. The left panel represents the PPI predicted with the coDASG_DEGs. These genes are divided and colored into different clusters according to K-means. The right panel shows the biological process terms predicted with the genes in different clusters. The color of the cluster on the right panel is the same as the color of the gene on the left panel.
Figure 2. Functional enrichment analysis of the coDASG_DEGs. (A) Significant KEGG pathway of coDASG_DEGs in adipose (P < 0.05). (B) PPI analysis of coDASG_DEGs in adipose tissue. (C) Significant KEGG pathway of coDASG_DEGs in blood (P < 0.05). (D) PPI analysis of coDASG_DEGs in blood. The left panel represents the PPI predicted with the coDASG_DEGs. These genes are divided and colored into different clusters according to K-means. The right panel shows the biological process terms predicted with the genes in different clusters. The color of the cluster on the right panel is the same as the color of the gene on the left panel.
Preprints 96858 g002
Preprints 96858 g003
Figure 3. Combined analysis of the coDASG_DEGs and TFs. (A, B) TFs in DEGs of adipose (A) and blood (B). The red and green bars denote the up-regulation and down-regulation of TF in the Eastern pig, respectively. The red font indicates that the differentially expressed TF also undergoes differentially alternative splicing. (C, D) The interaction network of TF and the coDASG_DEGs in adipose (C) and blood (D). TFs and coDASG_DEGs with a correlation coefficient (r) greater than 0.9 and a P value less than 3.4E-14 were screened and made into an interaction network. The red and blue dots separately represent TF and genes.
Figure 3. Combined analysis of the coDASG_DEGs and TFs. (A, B) TFs in DEGs of adipose (A) and blood (B). The red and green bars denote the up-regulation and down-regulation of TF in the Eastern pig, respectively. The red font indicates that the differentially expressed TF also undergoes differentially alternative splicing. (C, D) The interaction network of TF and the coDASG_DEGs in adipose (C) and blood (D). TFs and coDASG_DEGs with a correlation coefficient (r) greater than 0.9 and a P value less than 3.4E-14 were screened and made into an interaction network. The red and blue dots separately represent TF and genes.
Preprints 96858 g003
Preprints 96858 g004
Figure 4. Genes under selection (SGs) in coDASG_DEGs of adipose and blood tissues. (A, B) The overlapped genes of SGs and coDASG_DEGs in adipose (A) and blood (B). Est and Wst represent SGs in Eastern and Western pigs, respectively. Red dots with positive values and green dots with negative values indicate the genes are up-regulated and down-regulated in Eastern pigs, respectively. The genes in red font are TFs.
Figure 4. Genes under selection (SGs) in coDASG_DEGs of adipose and blood tissues. (A, B) The overlapped genes of SGs and coDASG_DEGs in adipose (A) and blood (B). Est and Wst represent SGs in Eastern and Western pigs, respectively. Red dots with positive values and green dots with negative values indicate the genes are up-regulated and down-regulated in Eastern pigs, respectively. The genes in red font are TFs.
Preprints 96858 g004
Preprints 96858 g005
Figure 5. Pearson correlation analysis of SFs and coDASG_DEGs. Genes with red and green backgrounds were positively and negatively correlated with SFs, respectively. The number in the figure represents the correlation coefficient. The positive and negative values represent positive and negative correlations, accompanied by the red and green bars. And the length of the bars represents the absolute value of the correlation coefficient. The number in bold font indicates a significant correlation (P < 0.05).
Figure 5. Pearson correlation analysis of SFs and coDASG_DEGs. Genes with red and green backgrounds were positively and negatively correlated with SFs, respectively. The number in the figure represents the correlation coefficient. The positive and negative values represent positive and negative correlations, accompanied by the red and green bars. And the length of the bars represents the absolute value of the correlation coefficient. The number in bold font indicates a significant correlation (P < 0.05).
Preprints 96858 g005
Preprints 96858 g006
Figure 6. Comparison of NSRP1 and its related genes with QTLs related to adipose traits. (A) The information of QTLs and the related genes. The genes in the figure are splicing factor NSRP1 and the genes highly correlated with it (r > 0.86, P < 6.104E-12). Green and rose-red squares represent genes present in the corresponding QTLs. The numbers on the top indicate the number of QTLs in which a specific gene is present, and the numbers on the right indicate the number of genes present in a specific QTL. (B) Comparison of the expression level of these genes in Eastern and Western pigs. Est: Eastern, Wst: Western. ** represents the P-value is less than 0.01, *** means the P-value is less than 0.001.
Figure 6. Comparison of NSRP1 and its related genes with QTLs related to adipose traits. (A) The information of QTLs and the related genes. The genes in the figure are splicing factor NSRP1 and the genes highly correlated with it (r > 0.86, P < 6.104E-12). Green and rose-red squares represent genes present in the corresponding QTLs. The numbers on the top indicate the number of QTLs in which a specific gene is present, and the numbers on the right indicate the number of genes present in a specific QTL. (B) Comparison of the expression level of these genes in Eastern and Western pigs. Est: Eastern, Wst: Western. ** represents the P-value is less than 0.01, *** means the P-value is less than 0.001.
Preprints 96858 g006
Preprints 96858 g007
Figure 7. Effect of NSRP1 on proliferation and differentiation of pre-adipocyte and lipid droplet formation. (A) The RT-qPCR results demonstrate the interference efficiency of three distinct small interfering RNAs targeting NSRP1 during both cell proliferation (top panel) and differentiation (bottom panel) stages. (B) Western blot experiment verified the interference efficiency of si-NSRP1. (C) Cell proliferation rates in response to downregulated NSRP1 expression were investigated through CCK8 assay. (D) The mRNA expression levels of pre-adipocyte proliferation marker genes (Ki67 and PCNA) were evaluated following the downregulation of NSRP1 expression. (E) The protein expression levels of pre-adipocyte proliferation marker genes (Ki67 and PCNA) were evaluated following the downregulation of NSRP1 expression. (F) EdU assay was carried out after si-NSRP1 transfection for 24 h. Cells undergoing DNA replication were stained by EdU (red) and cell nuclei were stained with DAPI (blue). Scale bar, 200 µm. (G) The percentage of DNA+ nuclei in Figure F was quantified using ImageJ. (H) The mRNA expression levels of pre-adipocyte differentiation (adipogenesis) marker genes (PPARG and CEBPA) were evaluated following the downregulation of NSRP1 expression. (I) The protein expression levels of pre-adipocyte differentiation (adipogenesis) marker genes (PPARG and CEBPA) were evaluated following the downregulation of NSRP1 expression. (J) The cells were subjected to adipogenic differentiation for 10 days downregulating NSRP1 expression, and the assessment of lipid droplet formation was conducted using Oil Red O staining. Scale bar, 200 µm.
Figure 7. Effect of NSRP1 on proliferation and differentiation of pre-adipocyte and lipid droplet formation. (A) The RT-qPCR results demonstrate the interference efficiency of three distinct small interfering RNAs targeting NSRP1 during both cell proliferation (top panel) and differentiation (bottom panel) stages. (B) Western blot experiment verified the interference efficiency of si-NSRP1. (C) Cell proliferation rates in response to downregulated NSRP1 expression were investigated through CCK8 assay. (D) The mRNA expression levels of pre-adipocyte proliferation marker genes (Ki67 and PCNA) were evaluated following the downregulation of NSRP1 expression. (E) The protein expression levels of pre-adipocyte proliferation marker genes (Ki67 and PCNA) were evaluated following the downregulation of NSRP1 expression. (F) EdU assay was carried out after si-NSRP1 transfection for 24 h. Cells undergoing DNA replication were stained by EdU (red) and cell nuclei were stained with DAPI (blue). Scale bar, 200 µm. (G) The percentage of DNA+ nuclei in Figure F was quantified using ImageJ. (H) The mRNA expression levels of pre-adipocyte differentiation (adipogenesis) marker genes (PPARG and CEBPA) were evaluated following the downregulation of NSRP1 expression. (I) The protein expression levels of pre-adipocyte differentiation (adipogenesis) marker genes (PPARG and CEBPA) were evaluated following the downregulation of NSRP1 expression. (J) The cells were subjected to adipogenic differentiation for 10 days downregulating NSRP1 expression, and the assessment of lipid droplet formation was conducted using Oil Red O staining. Scale bar, 200 µm.
Preprints 96858 g007
Preprints 96858 g008
Figure 8. NSRP1 modulates adipogenesis through the regulation of AS and expression of CCDC18. (A) Semiquantitative RT-PCR analyses were conducted to assess splicing changes of CCDC18 in representative Western pig breed (Duroc) and Eastern pig breed (Luchuan) adipose tissues. (B) Semiquantitative RT-PCR analyses were performed to evaluate splicing changes of CCDC18 following downregulation of NSRP1 expression. Left panel, the agarose gel electrophoresis results. Right panel, the quantification of alternative splicing levels through ImageJ processing of grayscale values. (C) The RT-qPCR experiment verified the effect of downregulating NSRP1 on CCDC18. Left panel, proliferation stage. Right panel, differentiation stage. (D) The RT-qPCR results demonstrate the interference efficiency of three distinct small interfering RNAs targeting CCDC18. (E) The mRNA expression levels of cell proliferation and differentiation marker genes after si-CDDC18. (F) EdU assay was carried out after transfection si-CCDC18 for 24 h. Cells undergoing DNA replication were stained by EdU (red) and cell nuclei were stained with DAPI (blue). Scale bar 200 µm. (G) The percentage of DNA+ nuclei in Figure F was quantified using ImageJ. (H) The cells were subjected to adipogenic differentiation for 10 days downregulating NSRP1 expression, and the assessment of lipid droplet formation was conducted using Oil Red O staining. Scale bar, 200 µm.
Figure 8. NSRP1 modulates adipogenesis through the regulation of AS and expression of CCDC18. (A) Semiquantitative RT-PCR analyses were conducted to assess splicing changes of CCDC18 in representative Western pig breed (Duroc) and Eastern pig breed (Luchuan) adipose tissues. (B) Semiquantitative RT-PCR analyses were performed to evaluate splicing changes of CCDC18 following downregulation of NSRP1 expression. Left panel, the agarose gel electrophoresis results. Right panel, the quantification of alternative splicing levels through ImageJ processing of grayscale values. (C) The RT-qPCR experiment verified the effect of downregulating NSRP1 on CCDC18. Left panel, proliferation stage. Right panel, differentiation stage. (D) The RT-qPCR results demonstrate the interference efficiency of three distinct small interfering RNAs targeting CCDC18. (E) The mRNA expression levels of cell proliferation and differentiation marker genes after si-CDDC18. (F) EdU assay was carried out after transfection si-CCDC18 for 24 h. Cells undergoing DNA replication were stained by EdU (red) and cell nuclei were stained with DAPI (blue). Scale bar 200 µm. (G) The percentage of DNA+ nuclei in Figure F was quantified using ImageJ. (H) The cells were subjected to adipogenic differentiation for 10 days downregulating NSRP1 expression, and the assessment of lipid droplet formation was conducted using Oil Red O staining. Scale bar, 200 µm.
Preprints 96858 g008
Table 1. Tissue-specific differentially alternative splicing genes and differentially expressed genes.
Table 1. Tissue-specific differentially alternative splicing genes and differentially expressed genes.
Adipose Blood Heart Kidney Lung Muscle Ovary Spleen
DASGs 4187 2898 173 242 315 244 584 1721
SDASGs 1740 909 41 63 80 56 101 280
SDASGs/DASGs 41.56% 31.37% 23.70% 26.03% 25.40% 22.95% 17.29% 16.27%
DEGs 840 960 178 218 358 310 1348 620
SDEGs 504 604 69 109 210 161 880 311
SDEGs/DEGs 60% 62.92% 38.76% 50% 58.66% 51.94% 65.28% 50.16%
coDASG_DEGs 83 164 0 4 4 2 15 65
DASGs=Differentially Alternative Splicing Genes; SDASGs=Tissue-specific DASGs; DEGs=Differentially Expressed Genes; SDEGs=Tissue-specific DEGs; DASGs: FDR < 0.05, |IncLevelDifference| > 5%; DEGs: P < 0.05, |Fold change| > 2; coDASG_DEGs: Common genes between DASGs and DEGs.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

facebook logotwitter logolinkedin logo
MDPI Initiatives
Important Links
Subscribe

Disclaimer

Terms of Use

Privacy Policy

© 2024 MDPI (Basel, Switzerland) unless otherwise stated