Wei Zhang

Genetic variation facilitates the evolution, environmental adaptability, and biodiversity of organisms. Danish Large White (LW) pigs have more desirable phenotypes compared with local Chinese pigs, which have difficulty adapting to the modern swine industry. However, the genome-wide mutational differences between these pig breeds are yet to be evaluated. Therefore, this study aimed to evaluate genomic variation and identify breed-specific SNPs in Danish LW pigs. Here, 43 LW, 15 Diqing Tibetan (DQZ), and 15 Diannan small-ear (DN) pigs whose genomes were re-sequenced with 5× depth were selected. This was followed by a conjoined analysis of our previous resequencing data of 24 Anqing six-end white (AQ) and six Asian wild (SS) pigs. In total, 39,158,378 SNPs and 13,143,989 insertion–deletions were obtained in all breeds. The variation number of LW pigs was the lowest, with 287,194 breed-specific and 1289 non-synonymous SNPs compared with Chinese breeds. Functional analysis of the breed-specific non-synonymous SNPs indicated that these mutations were mainly associated with the reproductive performance, feed intake, and feed conversion ratio of LW pigs. These findings provide a theoretical basis for genetic improvements in the Chinese swine industry.

Baoshan pigs (BS) are a local breed in Yunnan Province that may face inbreeding owing to its limited population size. To accurately evaluate the inbreeding level of the BS pig population, we used whole-genome resequencing to identify runs of homozygosity (ROH) regions in BS pigs, calculated the inbreeding coefficient based on pedigree and ROH, and screened candidate genes with important economic traits from ROH islands. A total of 22,633,391 SNPS were obtained from the whole genome of BS pigs, and 201 ROHs were detected from 532,450 SNPS after quality control. The number of medium-length ROH (1–5 Mb) was the highest (98.43%), the number of long ROH (>5 Mb) was the lowest (1.57%), and the inbreeding of BS pigs mainly occurred in distant generations. The inbreeding coefficient F_ROH, calculated based on ROH, was 0.018 ± 0.016, and the F_PED, calculated based on the pedigree, was 0.027 ± 0.028, which were positively correlated. Forty ROH islands were identified, containing 507 genes and 891 QTLs. Several genes were associated with growth and development (IGFALS, PTN, DLX5, DKK1, WNT2), meat quality traits (MC3R, ACSM3, ECI1, CD36, ROCK1, CACNA2D1), and reproductive traits (NPW, TSHR, BMP7). This study provides a reference for the protection and utilization of BS pigs.

Copy number variation (CNV) has been widely used to study the evolution of different species. We first discovered different CNVs in 24 Anqingliubai pigs and 6 Asian wild boars using next-generation sequencing at the whole-genome level with 10× depth to understand the relationship between genetic evolution and production traits in wild boars and domestic pigs. A total of 97,489 CNVs were identified and divided into 10,429 copy number variation regions (CNVRs), occupying 32.06% of the porcine genome. Chromosome 1 had the most CNVRs, and chromosome 18 had the least. Ninety-six CNVRs were selected using VST 1% based on the signatures of all CNVRs, and sixty-five genes were identified in the selected regions. These genes were strongly correlated with traits distinguishing groups by enrichment in Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathways, such as growth (CD36), reproduction (CIT, RLN), detoxification (CYP3A29), and fatty acid metabolism (ELOVL6). The QTL overlapping regions were associated with meat traits, growth, and immunity, which was consistent with CNV analysis. Our findings increase the understanding of evolved genome structural variations between wild boars and domestic pigs, and provide new molecular biomarkers to guide breeding and the efficient use of available genetic resources.

Intramuscular fat (IMF) deposition is an important determinant of pork quality and a complex process facilitated by non-coding ceRNAs. In this study, 52 Berkshire × Anqing Sixwhite crossbred pigs were slaughtered to measure eight carcass and pork quality traits. Whole-transcriptome sequencing analysis was performed using longissimus dorsi samples of six low- and high-IMF samples; 34 ceRNA networks, based on 881, 394, 158 differentially expressed (DE) lncRNAs, miRNAs, and mRNAs, were constructed. Following weighted gene co-expression network analysis between the low and high IMF, only one ceRNA, lncRNA4789/miR-381-3p/FABP3, that showed similar DE trend in longissimus dorsi tissue was retained. Dual-luciferase reporter assays further indicated that FABP3 was a direct, functional target of miR-381-3p, where miR-381-3p overexpression inhibited the mRNA and protein expression of FABP3. In addition, overexpressed lncRNA4789 attenuated the effect of miR-381-3p on FABP3 by sponging miR-381-3p. Cell function verification experiment demonstrated that miR-381-3p suppressed IMF deposition by inhibiting preadipocyte cell differentiation and lipid droplet deposition via the suppression of FABP3 expression in the peroxisome proliferator-activated receptor signalling pathway, whereas lncRNA4789 rescued FABP3 expression by sponging miR-381-3p. Our study may aid in identifying novel molecular markers for its optimization in IMF which is of importance in breeding for improving pork quality.