Identification of causative genes or genetic variants associated with phenotype traits benefits the genetic improvement of animals. CISH plays a role in immunity and growth, however, the upstream transcriptional factors of porcine CISH and the genetic variations in these factors remain unclear. In this study, we firstly identified the minimal core promoter of porcine CISH and confirmed the existence of STATx binding sites. Overexpression and RT-qPCR demonstrated STAT5A increased CISH transcriptional activity (P < 0.01) and mRNA expression (P < 0.01), while GATA1 inhibited CISH transcriptional activity (P < 0.01) and the following mRNA expression (P < 0.05 or P < 0.01). Then, the putative functional genetic variations of porcine STAT5A were screened and a PCR-SSCP was established for genotype g.508A>C and g.566C>T. Population genetic analysis showed the A allele frequency of g.508A>C and C allele frequency of g.566C>T was 0.61 and 0.94 in Min pigs, respectively, while these two alleles were fixed in the Landrace population. Statistical analysis showed that Min piglets with CC genotype at g.566C>T or Hap1: AC had higher 28-day body weight, 35-day body weight, and ADG than TC or Hap3: CT animals (P < 0.05, P < 0.05). Further luciferase activity assay demonstrated that the activity of g.508A>C in the C allele was lower than the A allele (P < 0.05). Collectively, the present study demonstrated that STAT5A positively regulated porcine CISH transcription, and SNP g.566C>T in the STAT5A was associated with the Min piglet growth trait.

cis-Splicing of adjacent genes (cis-SAGe) has been involved in multiple physiological and pathological processes in humans. However, to the best of our knowledge, there is no report of cis-SAGe in adipogenic regulation. In this study, a cis-SAGe product, BCL2L2-PABPN1 (BP), was characterized in fat tissue of pigs with RT-PCR and RACE method. BP is an in-frame fusion product composed of 333 aa and all the functional domains of both parents. BP is highly conserved among species and rich in splicing variants. BP was found to promote proliferation and inhibit differentiation of primary porcine preadipocytes. A total of 3074/44 differentially expressed mRNAs (DEmRs)/known miRNAs (DEmiRs) were identified in porcine preadipocytes overexpressing BP through RNA-Seq analysis. Both DEmRs and target genes of DEmiRs were involved in various fat-related pathways with MAPK and PI3K-Akt being the top enriched. PPP2CB, EGFR, Wnt5A and EHHADH were hub genes among the fat-related pathways identified. Moreover, ssc-miR-339-3p was found to be critical for BP regulating adipogenesis through integrated analysis of mRNA and miRNA data. The results highlight the role of cis-SAGe in adipogenesis and contribute to further revealing the mechanisms underlying fat deposition, which will be conductive to human obesity control.

Meat quality is one of the most important economic traits in pig breeding and production. Intramuscular fat (IMF) is a major factor that improves meat quality. To better understand the alternative splicing (AS) events underlying meat quality, long-read isoform sequencing (Iso-seq) was used to identify differential (D)AS events between the longissimus thoracis (LT) and semitendinosus (ST), which differ in IMF content, together with short-read RNA-seq. Through Iso-seq analysis, we identified a total of 56,789 novel transcripts covering protein-coding genes, lncRNA, and fusion transcripts that were not previously annotated in pigs. We also identified 456,965 AS events, among which 3930 were DAS events, corresponding to 2364 unique genes. Through integrative analysis of Iso-seq and RNA-seq, we identified 1174 differentially expressed genes (DEGs), among which 122 were DAS genes, i.e., DE-DAS genes. There are 12 overlapped pathways between the top 20 DEGs and DE-DAS genes, as revealed by KEGG (Kyoto Encyclopedia of Genes and Genomes) analysis, indicating that DE-DAS genes play important roles in the differential phenotype of LT and ST. Further analysis showed that upregulated DE-DAS genes are more important than downregulated ones in IMF deposition. Fatty acid degradation and the PPAR (peroxisome proliferator-activated receptor) signaling pathway were found to be the most important pathways regulating the differential fat deposition of the two muscles. The results update the existing porcine genome annotations and provide data for the in-depth exploration of the mechanisms underlying meat quality and IMF deposition.

The cAMP response element-binding protein (CREB), a basic leucine zipper transcription factor, is involved in the activation of numerous genes in a variety of cell types. The CREB gene is rich in alternative splicing (AS) events. However, studies on the AS of CREB genes in pigs are limited, and few reports have compared the roles of isoforms in activating gene expression. Here, five AS transcripts, V1-5, were characterized by RT-PCR and two, V3 and V5, were new identifications. Both V1 and V2 have all the functional domains of the CREB protein, with similar tissue expression profiles and mRNA stability, suggesting that they have similar roles. The transcriptional transactivation activities of four isoforms encoding complete polypeptides were analyzed on the expression of the B-cell CLL/lymphoma 2-like protein 2 and the poly (A)-binding protein, nuclear 1 genes with a dual-luciferase reporter system, and differential activities were observed. Both V1 and V2 have promoting effects, but their roles are gene-specific. V3 has no effect on the promoter of the two genes, while V4 functions as a repressor. The mechanisms underlying the differential roles of V1 and V2 were analyzed with RNA-seq, and the genes specifically regulated by V1 and V2 were identified. These results will contribute to further revealing the role of CREB and to analyzing the significance of AS in genes.

The Min pig (Sus scrofa) is a well-known indigenous breed in China. One of its main advantages over European breeds is its high meat quality. Additionally, different cuts of pig also show some different traits of meat quality. To explore the underlying mechanism responsible for the differences of meat quality between different breeds or cuts, the longissimus dorsi muscle (LM) and the biceps femoris muscle (BF) from Min and Large White pigs were investigated using transcriptome analysis. The gene expression profiling identified 1371 differentially expressed genes (DEGs) between LM muscles from Min and Large White pigs, and 114 DEGs between LM and BF muscles from the same Min pigs. Gene Ontology (GO) enrichment of biological functions and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed that the gene products were mainly involved in the IRS1/Akt/FoxO1 signaling pathway, adenosine 5'-monophosphate-activated protein kinase (AMPK) cascade effects, lipid metabolism and amino acid metabolism pathway. Such pathways contributed to fatty acid metabolism, intramuscular fat deposition, and skeletal muscle growth in Min pig. These results give an insight into the mechanisms underlying the formation of skeletal muscle and provide candidate genes for improving meat quality. It will contribute to improving meat quality of pigs through molecular breeding.

Matrix metalloproteinase 9 (MMP9) plays critical roles in multiple biological processes, such as reproduction, cell proliferation and differentiation, and host defenses. The aim of this study was to evaluate whether MMP9 is a candidate gene for resistance to diarrhea in piglets. In this study, quantitative real-time PCR was used to analyze the expression of MMP9 mRNA in different tissues of specific pathogen-free piglets. MMP9 was expressed in all the tissues (heart, liver, spleen, lung, kidney, stomach, duodenum, jejunum, ileum, and colon) analyzed. An association analysis between MMP9 polymorphisms and piglet diarrhea score and performance traits were performed in Min (Chinese indigenous breed) and Landrace populations. In the statistical analysis, at the g.48178429 G>A locus, AA piglets had a lower diarrhea score than that of GA in the Min population (P < .05), whereas GG had higher day-35 body weight and average daily gain (ADG) than AA in the Landrace breed (P < .05). At the rs336583561 locus, Min piglets with the GG genotype have a lower diarrhea score than AG piglets (P < .05). At g.48184777C>T, CC animals have higher body weight than TC Landrace piglets (P < .05 or P < .01). A 5' flanking deletion assay indicated that g.48178429 G>A was not located in the MMP9 promoter region. Our results suggest that the A allele at the g.48178429 G>A locus and the G allele at rs336583561 are resistance alleles in Min pigs. Before these markers are used in pig breeding programs, more studies in larger populations are needed.

Variations in Cytokine inducible SH2-containing protein (CISH) gene influence human susceptibility to common infectious diseases, but little is known about CISH in swine. The objectives of this study were to 1) determine porcine CISH (pCISH) mRNA expression level in different tissues of piglets, 2) predict putative functional genetic variations within pCISH, 3) investigate the association between a identified variation in the 3'UTR and piglets phenotype traits in Min (n = 226) and Landrace (n = 186) population, and explore the function of this variation. Results of quantitative PCR showed pCISH mRNA expressed in all the collected tissues with higher level in lung and ileum than colon (p < 0.05). In-silico analysis indicated none of the functional ns-SNPs existed in pCISH coding region. Results from the characterizing of 3'UTR presented a novel 12-bp insertion/deletion (indel) mutation. Statistical analysis demonstrated that this 12-bp indel associated with piglets diarrhea score in the Landrace population, and animals with AA genotype (12-bp insertion) presented lower diarrhea score when compared with BB (p < 0.05) or AB (p < 0.01) carriers. The in vitro study indicated that the luciferase activity of reconstruct plasmid psiCHECK-2-CISH-AA or psiCHECK-2-CISH-BB was significantly lower than the negative control (p < 0.05), and luciferase activity of psiCHECK-2-CISH-AA was higher than that of the psiCHECK-2-CISH-BB (p < 0.05). Although results herein suggested the 12-bp indel might affect Landrace piglet susceptibility to diarrhea, further association studies in more populations are needed before this preliminary finding could be used for pig breeding.

Matrix metalloproteinase 7 (MMP7) is involved in the degradation of extracellular matrix in disease processes and therefore plays an important role in host disease resistance/susceptibility. To better understanding the effects of porcine MMP7 (pMMP7) on piglets diarrhea trait, we characterized pMMP7 gene, identified genetic variations in pMMP7 and explored the relationship between pMMP7 polymorphisms and piglets diarrhea in Min pig and Landrace populations. The complete coding sequence of pMMP7 is 804 bp encoding a protein of 267 amino acids. Sequence alignment showed that the identity between pMMP7 and human MMP7 was approximately 80%. The expression of pMMP7 in the gut of healthy piglets were weak and the distribution of the pMMP7-EGFP fusion protein was observed mainly in the cytoplasm. After the identification of 21 genetic variations in 5' flanking region and exons, Hae III and Eco72 Ⅰ PCR-RFLP were established to genotype SNP rs327380117 and rs329429922, respectively. Statistical analysis indicated that Landrace piglets with a TT genotype at rs327380117 had a lower diarrhea score and day-14 wt than TC piglets (p < 0.05); the diarrhea score of AA Landrace animals with rs329429922 was lower than that of GG individuals (p < 0.05). The findings presented here contribute to the understanding of the biological function of pMMP7 and may provide new molecular markers for pig breeding.

Long-chain acyl-CoA synthetase 1 (ACSL1) plays an important role in fatty acid metabolism and fat deposition. The transcription of the ACSL1 gene is regulated specifically among cells and physiological processes, and transcriptional regulation of ACSL1 in adipogenesis remains elusive. Here, we characterize transcription factors (TFs) associated with adipogenesis in the porcine ACSL1 gene. CCAAT-enhancer binding protein (C/EBP)α, a well-known adipogenic marker, was found to enhance the expression of the ACSL1 gene via binding two tandem motifs in the promoter. Further, we demonstrate that ACSL1 mediates C/EBPα effects on adipogenesis in preadipocytes cultured from subcutaneous fat tissue of pigs via gain- and loss-of-function analyses. The cAMP-response element binding protein, another TF involved in adipogenesis, was also identified in the regulation of ACSL1 gene expression. Additionally, single nucleotide polymorphisms (SNPs) were screened in the promoter of ACSL1 among four breeds including the Chinese indigenous Min, and Duroc, Berkshire, and Yorkshire pigs through sequencing of PCR products. Two tightly linked SNPs, -517G>T and -311T>G, were found exclusively in Min pigs. The haplotype mutation decreases promoter activity in PK-15 and ST cells, and in vivo the expression of ACSL1, illustrating a possible role in adipogenesis regulated by C/EBPα/ACSL1 axis. Additionally, a total of 24 alternative splicing transcripts were identified, indicating the complexity of alternative splicing in the ACSL1 gene. The results will contribute to further revealing the regulatory mechanisms of ACSL1 during adipogenesis and to the characterization of molecular markers for selection of fat deposition in pigs.

IFN-induced protein with tetratricopeptide repeats 2 (IFIT2) plays important roles in host defense against viral infection as revealed by studies in humans and mice. However, little is known on porcine IFIT2 (pIFIT2). Here, we performed molecular cloning, expression profile, and transcriptional regulation analysis of pIFIT2. pIFIT2 gene, located on chromosome 14, is composed of two exons and have a complete coding sequence of 1407 bp. The encoded polypeptide, 468 aa in length, has three tetratricopeptide repeat motifs. pIFIT2 gene was unevenly distributed in all eleven tissues studied with the most abundance in spleen. Poly(I:C) treatment notably strongly upregulated the mRNA level and promoter activity of pIFIT2 gene. Upstream sequence of 1759 bp from the start codon which was assigned +1 here has promoter activity, and deltaEF1 acts as transcription repressor through binding to sequences at position - 1774 to - 1764. Minimal promoter region exists within nucleotide position - 162 and - 126. Two adjacent interferon-stimulated response elements (ISREs) and two nuclear factor (NF)-κB binding sites were identified within position - 310 and - 126. The ISRE elements act alone and in synergy with the one closer to start codon having more strength, so do the NF-κB binding sites. Synergistic effect was also found between the ISRE and NF-κB binding sites. Additionally, a third ISRE element was identified within position - 1661 to - 1579. These findings will contribute to clarifying the antiviral effect and underlying mechanisms of pIFIT2.

The large scale genome wide association studies (GWAS) have identified approximately 80 single nucleotide polymorphisms (SNPs) conferring susceptibility to type 2 diabetes (T2D). However, most of these loci have not been replicated in diverse populations and much genetic heterogeneity has been observed across ethnic groups. We tested 28 SNPs previously found to be associated with T2D by GWAS in a Mongolian sample of Northern China (497 diagnosed with T2D and 469 controls) for association with T2D and diabetes related quantitative traits. We replicated T2D association of 11 SNPs, namely, rs7578326 (IRS1), rs1531343 (HMGA2), rs8042680 (PRC1), rs7578597 (THADA), rs1333051 (CDKN2), rs6723108 (TMEM163), rs163182 and rs2237897 (KCNQ1), rs1387153 (MTNR1B), rs243021 (BCL11A), and rs10229583 (PAX4) in our sample. Further, we showed that risk allele of the strongest T2D associated SNP in our sample, rs757832 (IRS1), is associated with increased level of TG. We observed substantial difference of T2D risk allele frequency between the Mongolian sample and the 1000G Caucasian sample for a few SNPs, including rs6723108 (TMEM163) whose risk allele reaches near fixation in the Mongolian sample. Further study of genetic architecture of these variants in susceptibility of T2D is needed to understand the role of these variants in heterogeneous populations.

The Toll-like receptor 5 (TLR5) recognizes flagellin of Gram-positive and -negative bacteria and plays an important role in the host defense system. Here, we surveyed single nucleotide polymorphisms (SNPs) in the coding sequence of the porcine TLR5 gene in 83 individuals from five pig breeds, these including Chinese local populations and Western commercial pig breeds. A total of 19 medium polymorphic SNPs (0.25 < PIC < 0.5) were identified, three of which were missense mutations that clustered within the extracellular domain of TLR5. One of the non-synonymous SNPs fell within a 228-amino acid region which has been shown to be important for flagellin recognition. Four SNPs were only found with high frequencies in Oriental pig breeds. The 19 SNPs were found in 30 haplotypes, one of which segregated at high frequency in all samples. Compared with Western pig breeds, Chinese local populations had higher genetic diversity and more haplotypes. Tajima's test showed no evidence for deviation from neutrality. The data provide useful information for future genetic marker characterization by means of disease association analysis and/or stimulating the mutation carrier with relevant ligands.

It is of positive significance to explore the mechanism of antioxidant and metabolic response of Canna indica under Cr stress mediated by rhizosphere niche. However, the mechanisms of recruitment and interaction of rhizosphere microorganisms in plants still need to be fully understood. This study combined physiology, microbiology, and metabolomics, revealing the interaction between C. indica and rhizosphere microorganisms under Cr stress. The results showed that Cr stress increased the content of malondialdehyde (MDA) and oxygen-free radicals (ROS) in plants. At the same time, the activities of antioxidant enzymes (SOD, POD, and APX) and the contents of glutathione (GSH) and soluble sugar were increased. In addition, Cr stress decreased the α diversity index of C. indica rhizosphere bacterial community and changed its community structure. The dominant bacteria, namely, Actinobacteriota, Proteobacteria, and Chloroflexi accounted for 75.16% of the total sequence. At the same time, with the extension of stress time, the colonization amount of rhizosphere-dominant bacteria increased significantly, and the metabolites secreted by roots were associated with the formation characteristics of Proteobacteria, Actinobacteria, Bacteroidetes, and other specific bacteria. Five critical metabolic pathways were identified by metabolome analysis, involving 79 differentially expressed metabolites, which were divided into 15 categories, mainly including lipids, terpenoids, and flavonoids. In conclusion, this study revealed the recruitment and interaction response mechanism between C. indica and rhizosphere bacteria under Cr stress through multi-omics methods, providing the theoretical basis for the remediation of Cr-contaminated soil.

The present study aimed to search for functional mutations within the promoter of porcine STAT3 and to provide causative genetic variants associated with piglet diarrhea. We firstly confirmed that STAT3 expressed higher in the small intestine than in the spleen, stomach and large intestine of SPF piglets, respectively (P < 0.05). Then, 10 genetic variations in the porcine STAT3 promoter region was identified by direct sequencing. Among them, three mutations SNP1: g.-870 G>A, SNP2: g.-584 A>C and a 6-bp Indel in the promoter region that displayed significant differential transcriptional activities were identified. Association analyses showed that SNP1: g.-870 G>A was significantly associated with piglet diarrhea (P < 0.05) and the GG animals had lower diarrhea score than AA piglets (P < 0.01) in both Min and Landrace population. Further functional analysis revealed that E2F6 repressed the transcriptional efficiency of STAT3 in vitro, by binding the G allele of SNP1. The present study suggested that SNP1: g.-870 G>A was a piglet diarrhea-associated variant that directly affected binding with E2F6, leading to changes in STAT3 transcription which might partially contribute to piglet diarrhea susceptibility or resistance.

Chimeric RNA was considered a special marker of cancer. However, recent studies have demonstrated that chimeric RNAs also exist in non-cancerous cells and tissues. Here, we analyzed and predicted jointly 49 chimeric RNAs by Star-Fusion and FusionMap. One chimeric RNA, we named TNNI2-ACTA1, and its eight transcript variants were identified by reverse transcriptase-polymerase chain reaction. The overexpression of TNNI2-ACTA1 V1 inhibited the proliferation of porcine skeletal muscle satellite cells through down-regulating the mRNA expression levels of cell cycle-related genes cyclinD1. However, as parental genes, there is no such effect in the TNNI2 and ACTA1. To explore the underlying mechanism for this phenomenon, we used RNA-seq to profile the transcriptomes of PSCs with overexpression. Compared with the negative control group, 1,592 differentially expressed genes (DEGs) were upregulated and 1,077 DEGs downregulated in TNNI2 group; 1,226 DEGs were upregulated and 902 DEGs downregulated in ACTA1 group; and 13 DEGs were upregulated and 16 DEGs downregulated in TNNI2-ACTA1 V1 group, respectively. Compared with the parental gene groups, three specific genes were enriched in the TNNI2-ACTA1 V1 group (NCOA3, Radixin, and DDR2). These three genes may be the key to TNNI2-ACTA1 V1 regulating cell proliferation. Taken together, our study explores the role of chimeric RNAs in normal tissues. In addition, our study as the first research provides the foundation for the mechanism of chimeric RNAs regulating porcine skeletal muscle growth.

Chimeric RNA is a crucial target for tumor diagnosis and drug therapy, also having its unique biological role in normal tissues. TNNI2-ACTA1-V1 (TA-V1), a chimeric RNA discovered by our laboratory in porcine muscle tissue, can inhibit the proliferation of Porcine Skeletal Muscle Satellite Cells (PSCs). The regulatory mechanism of TA-V1 in PSCs remains unclear, but we speculate that NCOA3, DDR2 and RDX may be the target genes of TA-V1. In this study, we explored the effects of NCOA3, DDR2 and RDX on cell viability and cell proliferation by CCK-8 assay, EdU staining and flow cytometry. Furthermore, the regulatory pathway of proliferation in PSCs mediated by TA-V1 through NCOA3 or CyclinD1 was elucidated by co-transfection and co-immunoprecipitation (Co-IP). The results revealed that overexpression of NCOA3 significantly increased cell viability and the expression level of CyclinD1, and also promotes cell proliferation by changing cells from the G1 phase to the S phase. In addition, inhibiting the expression of NCOA3 substantially reduced cell viability and inhibited cell proliferation. Overexpression of DDR2 and RDX had no significant effect on cell viability and proliferation. Co-transfection experiments showed that NCOA3 could rescue the proliferation inhibition of PSCs caused by TA-V1. Co-IP assay indicated that TA-V1 directly interacts with NCOA3. Our study explores the hypothesis that TA-V1 directly regulates NCOA3, indirectly regulating CyclinD1, thereby regulating PSCs proliferation. We provide new putative mechanisms of porcine skeletal muscle growth and lay the foundation for the study of chimeric RNA in normal tissues.