Comprehensive inbred variation discovery in Bama pigs using de novo assemblies

Assembly of novel sequences for Chinese domestic pigs reveals new genes and regulatory variants providing new insights into their diversity

There is an increasing understanding that a reference genome representing an individual cannot capture all the gene repertoire of a species. Here, we conduct a population-scale missing sequences detection of Chinese domestic pigs using whole-genome sequencing data from 534 individuals. We identify 132.41 Mb of sequences absent in the reference assembly, including eight novel genes. In particular, the breeds spread in Chinese high-altitude regions perform significantly different frequencies of new sequences in promoters than other breeds. Furthermore, we dissect the role of non-coding variants and identify a novel sequence inserted in the 3'UTR of the FMO3 gene, which may be associated with the intramuscular fat phenotype. This novel sequence could be a candidate marker for meat quality. Our study provides a comprehensive overview of the missing sequences in Chinese domestic pigs and indicates that this dataset is a valuable resource for understanding the diversity and biology of pigs.

Integrated Analysis of the Transcriptome and Microbial Diversity in the Intestine of Miniature Pig Obesity Model

Obesity, a key contributor to metabolic disorders, necessitates an in-depth understanding of its pathogenesis and prerequisites for prevention. Guangxi Bama miniature pig (GBM) offers an apt model for obesity-related studies. In this research, we used transcriptomics and 16S rRNA gene sequencing to discern the differentially expressed genes (DEGs) within intestinal (jejunum, ileum, and colon) tissues and variations in microbial communities in intestinal contents of GBM subjected to normal diets (ND) and high-fat, high-carbohydrate diets (HFHCD). After a feeding duration of 26 weeks, the HFHCD-fed experimental group demonstrated notable increases in backfat thickness, BMI, abnormal blood glucose metabolism, and blood lipid levels alongside the escalated serum expression of pro-inflammatory factors and a marked decline in intestinal health status when compared to the ND group. Transcriptomic analysis revealed a total of 1669 DEGs, of which 27 had similar differences in three intestinal segments across different groups, including five immune related genes: COL6A6, CYP1A1, EIF2AK2, NMI, and LGALS3B. Further, we found significant changes in the microbiota composition, with a significant decrease in beneficial bacterial populations within the HFHCD group. Finally, the results of integrated analysis of microbial diversity with transcriptomics show a positive link between certain microbial abundance (Solibacillus, norank_f__Saccharimonadaceae, Candidatus_Saccharimonas, and unclassified_f__Butyricicoccaceae) and changes in gene expression (COL6A6 and NMI). Overall, HFHCD appears to co-contribute to the initiation and progression of obesity in GBM by aggravating inflammatory responses, disrupting immune homeostasis, and creating imbalances in intestinal flora.

Uncovering the candidate genes related to sheep body weight using multi-trait genome-wide association analysis

In sheep, body weight is an economically important trait. This study sought to map genetic loci related to weaning weight and yearling weight. To this end, a single-trait and multi-trait genome-wide association study (GWAS) was performed using a high-density 600 K single nucleotide polymorphism (SNP) chip. The results showed that 43 and 56 SNPs were significantly associated with weaning weight and yearling weight, respectively. A region associated with both weaning and yearling traits (OARX: 6.74-7.04 Mb) was identified, suggesting that the same genes could play a role in regulating both these traits. This region was found to contain three genes (TBL1X, SHROOM2 and GPR143). The most significant SNP was Affx-281066395, located at 6.94 Mb (p = 1.70 × 10-17), corresponding to the SHROOM2 gene. We also identified 93 novel SNPs elated to sheep weight using multi-trait GWAS analysis. A new genomic region (OAR10: 76.04-77.23 Mb) with 22 significant SNPs were discovered. Combining transcriptomic data from multiple tissues and genomic data in sheep, we found the HINT1, ASB11 and GPR143 genes may involve in sheep body weight. So, multi-omic anlaysis is a valuable strategy identifying candidate genes related to body weight.

Genome-Wide Detection of Copy Number Variations and Evaluation of Candidate Copy Number Polymorphism Genes Associated With Complex Traits of Pigs

Copy number variation (CNV) has been considered to be an important source of genetic variation for important phenotypic traits of livestock. In this study, we performed whole-genome CNV detection on Suhuai (SH) (n = 23), Chinese Min Zhu (MZ) (n = 11), and Large White (LW) (n = 12) pigs based on next-generation sequencing data. The copy number variation regions (CNVRs) were annotated and analyzed, and 10,885, 10,836, and 10,917 CNVRs were detected in LW, MZ, and SH pigs, respectively. Some CNVRs have been randomly selected for verification of the variation type by real-time PCR. We found that SH and LW pigs are closely related, while MZ pigs are distantly related to the SH and LW pigs by CNVR-based genetic structure, PCA, VST, and QTL analyses. A total of 14 known genes annotated in CNVRs were unique for LW pigs. Among them, the cyclin T2 (CCNT2) is involved in cell proliferation and the cell cycle. The FA Complementation Group M (FANCM) is involved in defective DNA repair and reproductive cell development. Ten known genes annotated in 47 CNVRs were unique for MZ pigs. The genes included glycerol-3-phosphate acyltransferase 3 (GPAT3) is involved in fat synthesis and is essential to forming the glycerol triphosphate. Glutathione S-transferase mu 4 (GSTM4) gene plays an important role in detoxification. Eleven known genes annotated in 23 CNVRs were unique for SH pigs. Neuroligin 4 X-linked (NLGN4X) and Neuroligin 4 Y-linked (NLGN4Y) are involved with nerve disorders and nerve signal transmission. IgLON family member 5 (IGLON5) is related to autoimmunity and neural activities. The unique characteristics of LW, MZ, and SH pigs are related to these genes with CNV polymorphisms. These findings provide important information for the identification of candidate genes in the molecular breeding of pigs.

Gene Cloning, Tissue Expression Profiles and Antiviral Activities of Interferon-β from Two Chinese Miniature Pig Breeds

The porcine interferon (PoIFN) complex represents an ideal model for studying IFN evolution which has resulted from viral pressure during domestication. Bama and Banna miniature pigs are the two Chinese miniature pig breeds that have been developed as laboratory animal models for studying virus infection, pathogenesis, and vaccine evaluation. However, the PoIFN complex of such miniature pig breeds remains to be studied. In the present study, we cloned PoIFN-β genes from Bama and Banna miniature pigs, detected their PoIFN-β tissue expression profiles, prepared recombinant PoIFN-β (rPoIFN-β) using the E. coli expression system, and measured their antiviral activities against three different pig viruses. At the amino acid sequence level, PoIFN-βs of the two miniature pig breeds were identical, which shared 100% identity with that of Congjiang Xiang pigs, 99.4–100% identity with that of domestic pigs, and 99.5% identity with that of three species of African wild boars. The tissue expression profiles of PoIFN-β mRNA differed not only between the two miniature pig breeds but between miniature pigs and domestic pigs as well. The four promoter domains of PoIFN-β of the two miniature pig breeds were identical with that of humans, domestic pigs, and three species of African wild boars. The recombinant PoIFN-β prepared from the two miniature pig breeds showed dose-dependent pre-infection and post-infection antiviral activities against vesicular stomatitis virus, porcine respiratory and reproductive syndrome virus, and pig pseudorabies virus. This study provided evidence for the high sequence conservation of PoIFN-β genes within the Suidae family with different tissue expression profiles and antiviral activities.

Gene cloning, prokaryotic expression and antiviral activities of interferon-αω from Chinese Bama miniature pigs

Porcine interferon (PoIFN) complex represents an ideal model for studying IFN evolution that resulted from viral pressure during domestication. IFN-αω is an emergent subtype of type I IFNs which has been primarily characterized in domestic pigs. In this study, the PoIFN-αω cDNA was cloned from Chinese Bama miniature pigs by RT-PCR, and its tissue expression profile was analyzed by real-time RT-PCR. The cDNA was expressed in Escherichia coli as a His-tagged protein and purified by nickel affinity chromatography. The antiviral activities of recombinant PoIFN-αω (rPoIFN-αω) against four different pig viruses were measured using cytopathic effect (CPE) inhibition assay. Although the PoIFN-αω sequence of Bama miniature pigs was identical to that of domestic pigs, the tissue expression profiles differed significantly between the two pig species. The rPoIFN-αω showed dose-dependent pre-infection antiviral activities against porcine pseudorabies virus, vesicular stomatitis virus and porcine reproductive and respiratory syndrome virus, but not against porcine circovirus type 2. When used as treatment post infection with the three viruses, rPoIFN-αω showed the efficacy in decreasing CPE in the infected cells in a time-dependent manner. Therefore, the expressed rPoIFN-αω could be used as an antiviral agent against pig virus infections.

Analysis of long non-coding RNAs in skeletal muscle of Bama Xiang pigs in response to heat stress

The aim of this study was to identify differentially expressed long non-coding RNAs (lncRNAs) molecules and predict their target genes related to muscle development and lipid metabolism in longissimus dorsi (LD) muscles of Bama Xiang pigs under constant heat stress. Ten male Bama Xiang pigs with an average initial body weight of 14 kg were randomly divided into control group (22°C) and heat stress (35 °C) group. The experiment lasted for 28 days. All the pigs were slaughtered at the end of the experiment, and LD muscles were collected for muscle quality analysis and transcriptome sequencing. Heat stress reduced meat quality of Bama Xiang pigs. lncRNAs in LD were identified systematically by deep RNA sequencing between the two groups. The results showed that 365 lncRNAs from the LD were identified, including 128 intergenic lncRNAs, 82 intronic lncRNAs, and 155 anti-sense lncRNAs. The differences lie in transcript of length, number of exons and wider size distribution, and expression level per KB fragment in three subtypes of lncRNAs. The three types of transposable elements coverage, including Line/L1, SINE/tRNA, and LTR/ERVL-MaLR, are the highest in mRNA and the three subtypes of lncRNAs in pigs. lncRNAs and mRNAs were different in comparison of features. The results predicted the target genes of the significant differentially expressed lncRNAs related to muscle development and lipid metabolism. This is the first study to expand the knowledge about muscle-related lncRNAs biology in Bama Xiang pigs under heat stress and will contribute to the development of alleviating the adverse effects of heat stress on pork quality targeting lncRNAs.