Multi-omic characterization of allele-specific regulatory variation in hybrid pigs.

Hybrid mapping is a powerful approach to efficiently identify and characterize genes regulated through mechanisms in cis. In this study, using reciprocal crosses of the phenotypically divergent Duroc and Lulai pig breeds, we perform a comprehensive multi-omic characterization of regulatory variation across the brain, liver, muscle, and placenta through four developmental stages. We produce one of the largest multi-omic datasets in pigs to date, including 16 whole genome sequenced individuals, as well as 48 whole genome bisulfite sequencing, 168 ATAC-Seq and 168 RNA-Seq samples. We develop a read count-based method to reliably assess allele-specific methylation, chromatin accessibility, and RNA expression. We show that tissue specificity was much stronger than developmental stage specificity in all of DNA methylation, chromatin accessibility, and gene expression. We identify 573 genes showing allele specific expression, including those influenced by parent-of-origin as well as allele genotype effects. We integrate methylation, chromatin accessibility, and gene expression data to show that allele specific expression can be explained in great part by allele specific methylation and/or chromatin accessibility. This study provides a comprehensive characterization of regulatory variation across multiple tissues and developmental stages in pigs.

Single-Nucleus RNA-Seq Reveals Spermatogonial Stem Cell Developmental Pattern in Shaziling Pigs.

Normal testicular development ensures the process of spermatogenesis, which is a complex biological process. The sustained high productivity of spermatogenesis throughout life is predominantly attributable to the constant proliferation and differentiation of spermatogonial stem cells (SSCs). The self-renewal and differentiation processes of SSCs are strictly regulated by the SSC niche. Therefore, understanding the developmental pattern of SSCs is crucial for spermatogenesis. The Shaziling pig is a medium-sized indigenous pig breed originating from central China. It is renowned for its superior meat quality and early male sexual maturity. The spermatogenic ability of the boars is of great economic importance to the pig industry. To investigate testicular development, particularly the pattern of SSC development in Shaziling pigs, we used single-cell transcriptomics to identify gene expression patterns in 82,027 individual cells from nine Shaziling pig testes at three key postnatal developmental stages. We generated an unbiased cell developmental atlas of Shaziling pig testicular tissues. We elucidated the complex processes involved in the development of SSCs within their niche in the Shaziling pig. Specifically, we identified potential marker genes and cellular signaling pathways that regulate SSC self-renewal and maintenance. Additionally, we proposed potential novel marker genes for SSCs that could be used for SSC isolation and sorting in Shaziling pigs. Furthermore, by immunofluorescence staining of testicular tissues of different developmental ages using marker proteins (UCHL1 and KIT), the developmental pattern of the spermatogonia of Shaziling pigs was intensively studied. Our research enhances the comprehension of the development of SSCs and provides a valuable reference for breeding Shaziling pigs.

Genome-Wide Identification of Specific Genetic Loci Common to Sheep and Goat.

Sheep and goat may become carriers of some zoonotic diseases. They are important livestock and experimental model animals for human beings. The fast and accurate identification of genetic materials originating from sheep and goat can prevent and inhibit the spread of some zoonotic diseases, monitor market product quality, and maintain the stability of animal husbandry and food industries. This study proposed a methodology for identifying sheep and goat common specific sites from a genome-wide perspective. A total of 150 specific sites were selected from three data sources, including the coding sequences of single copy genes from nine species (sheep, goat, cow, pig, dog, horse, human, mouse, and chicken), the dbSNPs for these species, and human 100-way alignment data. These 150 sites exhibited low intraspecific heterogeneity in the resequencing data of 1450 samples from five species (sheep, goat, cow, pig, and chicken) and high interspecific divergence in the human 100-way alignment data after quality control. The results were proven to be reliable at the data level. Using the process proposed in this study, specific sites of other species can be screened, and genome-level species identification can be performed using the screened sites.

Transcriptomic Analyses Reveal the Effects of Walnut Kernel Cake on Adipose Deposition in Pigs.

With the rising cost of animal feed protein, finding affordable and effective substitutes is crucial. Walnut kernel cake, a polyphenol-, fiber-, protein- and fat-rich byproduct of walnut oil extraction, has been underexplored as a potential protein replacement in pig feed. In this study, we found that feeding large Diqing Tibetan pigs walnut kernel cake promoted adipose deposition and improved pork quality during pig growth. Transcriptome analysis revealed the upregulation of genes ANGPTL8, CCNP, ETV4, and TRIB3, associated with adipose deposition. Pathway analysis highlighted enrichment in adipose deposition-related pathways, including PPAR, insulin, PI3K-Akt, Wnt, and MAPK signaling. Further analysis identified DEGs (differentially expressed genes) positively correlated with adipose-related traits, such as PER2 and PTGES. Single-cell transcriptome data pointed to the specific expression of CD248 and PTGES in adipocyte progenitor/stem cells (APSCs), pivotal for adipocyte differentiation and adipose deposition regulation. This study demonstrates walnut kernel cake's potential to substitute soybean cake in pig feed, providing high-quality protein and promoting adipose deposition. It offers insights into feed protein replacement, human functional food, fat metabolism, and related diseases, with marker genes and pathways supporting pig breeding and pork quality improvement.

Placental Transcriptome Analysis in Connection with Low Litter Birth Weight Phenotype (LBWP) Sows.

It is possible to identify sub-populations of sows in every pig herd that consistently give birth to low birth weight (BW) piglets, irrespective of the litter size. A previous study from our group demonstrated that placental development is a main factor affecting the litter birth weight phenotype (LBWP) in sows, thereby impacting the BW of entire litters, but the biological and molecular pathways behind this phenomenon are largely unknown. The aim of this study was to investigate the differential gene expression in placental tissues at day 30 of gestation between low LBWP (LLBWP) vs. high LBWP (HLBWP) sows from a purebred Large White maternal line. Using mRNA sequencing, we found 45 differentially expressed genes (DEGs) in placental tissues of LLBWP and HLBWP sows. Furthermore, (GO) enrichment of upregulated DEGs predicted that there were two biological processes significantly related to cornification and regulation of cell population proliferation. To better understand the molecular interaction between cell proliferation and cornification, we conducted transcriptional factor binding site (TFBS) prediction analysis. The results indicated that a highly significant TFBS was located at the 5' upstream of all four upregulated genes (CDSN, DSG3, KLK14, KRT17), recognized by transcription factors EGR4 and FOSL1. Our findings provide novel insight into how transcriptional regulation of two different biological processes interact in placental tissues of LLBWP sows.

Integration of transcriptome and machine learning to identify the potential key genes and regulatory networks affecting drip loss in pork.

Low level of drip loss (DL) is an important quality characteristic of meat with high economic value. However, the key genes and regulatory networks contributing to DL in pork remain largely unknown. To accurately identify the key genes affecting DL in muscles postmortem, 12 Duroc × (Landrace × Yorkshire) pigs with extremely high (n = 6, H group) and low (n = 6, L group) DL at both 24 and 48 h postmortem were selected for transcriptome sequencing. The analysis of differentially expressed genes and weighted gene co-expression network analysis (WGCNA) were performed to find the overlapping genes using the transcriptome data, and functional enrichment and protein-protein interaction (PPI) network analysis were conducted using the overlapping genes. Moreover, we used machine learning to identify the key genes and regulatory networks related to DL based on the interactive genes of the PPI network. Finally, nine potential key genes (IRS1, ESR1, HSPA6, INSR, SPOP, MSTN, LGALS4, MYLK2, and FRMD4B) mainly associated with the MAPK signaling pathway, the insulin signaling pathway, and the calcium signaling pathway were identified, and a single-gene set enrichment analysis (GSEA) was performed to further annotate the functions of these potential key genes. The GSEA results showed that these genes are mainly related to ubiquitin-mediated proteolysis and oxidative reactions. Taken together, our results indicate that the potential key genes influencing DL are mainly related to insulin signaling mediated differences in glycolysis and ubiquitin-mediated changes in muscle structure and improve the understanding of gene expression and regulation related to DL and contribute to future molecular breeding for improving pork quality.

A low level of drip loss (DL) is critical for the economic value of pork. However, the genetic basis underlying DL remains unclear. In this study, pigs with extremely high and low DL at both 24 and 48 h postmortem were selected, and total RNA from longissimus dorsi (LD) muscles was extracted for transcriptome sequencing. Subsequently, a variety of analytical methods, were integrated to identify the potential key genes and pathways affecting DL. As a result, nine potential key genes (IRS1, ESR1, HSPA6, INSR, SPOP, MSTN, LGALS4, MYLK2, and FRMD4B) mainly associated with the MAPK signaling pathway, insulin signaling pathway, and calcium signaling pathway, were identified, and these genes are primarily related to ubiquitin-mediated proteolysis and oxidation reactions. This study contributes new evidence for elucidating the molecular mechanism of DL and provides potential target genes for precise genetic improvement of DL.

Transcriptional Regulation Analysis Provides Insight into the Function of GSK3ß Gene in Diannan Small-Ear Pig Spermatogenesis.

Glycogen synthase kinase-3ß (GSK3ß) not only plays a crucial role in regulating sperm maturation but also is pivotal in orchestrating the acrosome reaction. Here, we integrated single-molecule long-read and short-read sequencing to comprehensively examine GSK3ß expression patterns in adult Diannan small-ear pig (DSE) testes. We identified the most important transcript ENSSSCT00000039364 of GSK3ß, obtaining its full-length coding sequence (CDS) spanning 1263 bp. Gene structure analysis located GSK3ß on pig chromosome 13 with 12 exons. Protein structure analysis reflected that GSK3ß consisted of 420 amino acids containing PKc-like conserved domains. Phylogenetic analysis underscored the evolutionary conservation and homology of GSK3ß across different mammalian species. The evaluation of the protein interaction network, KEGG, and GO pathways implied that GSK3ß interacted with 50 proteins, predominantly involved in the Wnt signaling pathway, papillomavirus infection, hippo signaling pathway, hepatocellular carcinoma, gastric cancer, colorectal cancer, breast cancer, endometrial cancer, basal cell carcinoma, and Alzheimer's disease. Functional annotation identified that GSK3ß was involved in thirteen GOs, including six molecular functions and seven biological processes. ceRNA network analysis suggested that DSE GSK3ß was regulated by 11 miRNA targets. Furthermore, qPCR expression analysis across 15 tissues highlighted that GSK3ß was highly expressed in the testis. Subcellular localization analysis indicated that the majority of the GSK3ß protein was located in the cytoplasm of ST (swine testis) cells, with a small amount detected in the nucleus. Overall, our findings shed new light on GSK3ß's role in DSE reproduction, providing a foundation for further functional studies of GSK3ß function.

Integrating Multiple Database Resources to Elucidate the Gene Flow in Southeast Asian Pig Populations.

The domestic pig (Sus scrofa) and its subfamilies have experienced long-term and extensive gene flow, particularly in Southeast Asia. Here, we analyzed 236 pigs, focusing on Yunnan indigenous, European commercial, East Asian, and Southeast Asian breeds, using the Pig Genomics Reference Panel (PGRP v1) of Pig Genotype-Tissue Expression (PigGTEx) to investigate gene flow and associated complex traits by integrating multiple database resources. In this study, we discovered evidence of admixtures from European pigs into the genome of Yunnan indigenous pigs. Additionally, we hypothesized that a potential conceptual gene flow route that may have contributed to the genetic composition of the Diannan small-ear pig is a gene exchange from the Vietnamese pig. Based on the most stringent gene introgression scan using the fd statistic, we identified three specific loci on chromosome 8, ranging from 51.65 to 52.45 Mb, which exhibited strong signatures of selection and harbored the NAF1, NPY1R, and NPY5R genes. These genes are associated with complex traits, such as fat mass, immunity, and litter weight, in pigs, as supported by multiple bio-functionalization databases. We utilized multiple databases to explore the potential dynamics of genetic exchange in Southeast Asian pig populations and elucidated specific gene functionalities.

Unveiling the Influence of Copy Number Variations on Genetic Diversity and Adaptive Evolution in China's Native Pig Breeds via Whole-Genome Resequencing.

Copy number variations (CNVs) critically influence individual genetic diversity and phenotypic traits. In this study, we employed whole-genome resequencing technology to conduct an in-depth analysis of 50 pigs from five local swine populations [Rongchang pig (RC), Wuzhishan pig (WZS), Tibetan pig (T), Yorkshire (YL) and Landrace (LR)], aiming to assess their genetic potential and explore their prospects in the field of animal model applications. We identified a total of 96,466 CNVs, which were subsequently integrated into 7112 non-redundant CNVRs, encompassing 1.3% of the swine genome. Functional enrichment analysis of the genes within these CNVRs revealed significant associations with sensory perception, energy metabolism, and neural-related pathways. Further selective scan analyses of the local pig breeds RC, T, WZS, along with YL and LR, uncovered that for the RC variety, the genes PLA2G10 and ABCA8 were found to be closely related to fat metabolism and cardiovascular health. In the T breed, the genes NCF2 and CSGALNACT1 were associated with immune response and connective tissue characteristics. As for the WZS breed, the genes PLIN4 and CPB2 were primarily linked to fat storage and anti-inflammatory responses. In summary, this research underscores the pivotal role of CNVs in fostering the diversity and adaptive evolution of pig breeds while also offering valuable insights for further exploration of the advantageous genetic traits inherent to China's local pig breeds. This facilitates the creation of experimental animal models tailored to the specific characteristics of these breeds, contributing to the advancement of livestock and biomedical research.

Full-Length Transcriptome Analysis of Skeletal Muscle of Jiangquan Black Pig at Different Developmental Stages.

Skeletal muscle grows in response to a combination of genetic and environmental factors, and its growth and development influence the quality of pork. Elucidating the molecular mechanisms regulating the growth and development of skeletal muscle is of great significance to both animal husbandry and farm management. The Jiangquan black pig is an excellent pig breed based on the original Yimeng black pig, importing the genes of the Duroc pig for meat traits, and cultivated through years of scientific selection and breeding. In this study, full-length transcriptome sequencing was performed on three growth stages of Jiangquan black pigs, aiming to study the developmental changes in Jiangquan black pigs at different developmental stages at the molecular level and to screen the key genes affecting the growth of skeletal muscle in Jiangquan black pigs. We performed an enrichment analysis of genes showing differential expression and constructed a protein-protein interaction network with the aim of identifying core genes involved in the development of Jiangquan black pigs. Notably, genes such as TNNI2, TMOD4, PLDIM3, MYOZ1, and MYH1 may be potential regulators of muscle development in Jiangquan black pigs. Our results contribute to the understanding of the molecular mechanisms of skeletal muscle development in this pig breed, which will facilitate molecular breeding efforts and the development of pig breeds to meet the needs of the livestock industry.

Genetic parameters for novel climatic resilience indicators derived from automatically-recorded vaginal temperature in lactating sows under heat stress conditions.

BACKGROUND: Longitudinal records of automatically-recorded vaginal temperature (TV) could be a key source of data for deriving novel indicators of climatic resilience (CR) for breeding more resilient pigs, especially during lactation when sows are at an increased risk of suffering from heat stress (HS). Therefore, we derived 15 CR indicators based on the variability in TV in lactating sows and estimated their genetic parameters. We also investigated their genetic relationship with sows' key reproductive traits. RESULTS: The heritability estimates of the CR traits ranged from 0.000 ± 0.000 for slope for decreased rate of TV (SlopeDe) to 0.291 ± 0.047 for sum of TV values below the HS threshold (HSUB). Moderate to high genetic correlations (from 0.508 ± 0.056 to 0.998 ± 0.137) and Spearman rank correlations (from 0.431 to 1.000) between genomic estimated breeding values (GEBV) were observed for five CR indicators, i.e. HS duration (HSD), the normalized median multiplied by normalized variance (Nor_medvar), the highest TV value of each measurement day for each individual (MaxTv), and the sum of the TV values above (HSUA) and below (HSUB) the HS threshold. These five CR indicators were lowly to moderately genetically correlated with shoulder skin surface temperature (from 0.139 ± 0.008 to 0.478 ± 0.048) and respiration rate (from 0.079 ± 0.011 to 0.502 ± 0.098). The genetic correlations between these five selected CR indicators and sow reproductive performance traits ranged from - 0.733 to - 0.175 for total number of piglets born alive, from - 0.733 to - 0.175 for total number of piglets born, and from - 0.434 to - 0.169 for number of pigs weaned. The individuals with the highest GEBV (most climate-sensitive) had higher mean skin surface temperature, respiration rate (RR), panting score (PS), and hair density, but had lower mean body condition scores compared to those with the lowest GEBV (most climate-resilient). CONCLUSIONS: Most of the CR indicators evaluated are heritable with substantial additive genetic variance. Five of them, i.e. HSD, MaxTv, HSUA, HSUB, and Nor_medvar share similar underlying genetic mechanisms. In addition, individuals with higher CR indicators are more likely to exhibit better HS-related physiological responses, higher body condition scores, and improved reproductive performance under hot conditions. These findings highlight the potential benefits of genetically selecting more heat-tolerant individuals based on CR indicators.

Genome-wide identification of quantitative trait loci and candidate genes for seven carcass traits in a four-way intercross porcine population.

BACKGROUND: Carcass traits are essential economic traits in the commercial pig industry. However, the genetic mechanism of carcass traits is still unclear. In this study, we performed a genome-wide association study (GWAS) based on the specific-locus amplified fragment sequencing (SLAF-seq) to study seven carcass traits on 223 four-way intercross pigs, including dressing percentage (DP), number of ribs (RIB), skin thinkness (ST), carcass straight length (CSL), carcass diagonal length (CDL), loin eye width (LEW), and loin eye thickness (LET). RESULTS: A total of 227,921 high-quality single nucleotide polymorphisms (SNPs) were detected to perform GWAS. A total of 30 SNPs were identified for seven carcass traits using the mixed linear model (MLM) (p < 1.0 × 10- 5), of which 9 SNPs were located in previously reported quantitative trait loci (QTL) regions. The phenotypic variation explained (PVE) by the significant SNPs was from 2.43 to 16.32%. Furthermore, 11 candidate genes (LYPLAL1, EPC1, MATN2, ZFAT, ZBTB10, ZNF704, INHBA, SMYD3, PAK1, SPTBN2, and ACTN3) were found for carcass traits in pigs. CONCLUSIONS: The GWAS results will improve our understanding of the genetic basis of carcass traits. We hypothesized that the candidate genes associated with these discovered SNPs would offer a biological basis for enhancing the carcass quality of pigs in swine breeding.

Whole-genome SNP allele frequency differences between Tibetan and Large white pigs reveal genes associated with skeletal muscle growth.

BACKGROUND: The skeletal muscle growth rate and body size of Tibetan pigs (TIB) are lower than Large white pigs (LW). However, the underlying genetic basis attributing to these differences remains uncertain. To address this knowledge gap, the present study employed whole-genome sequencing of TIB (slow growth) and LW (fast growth) individuals, and integrated with existing NCBI sequencing datasets of TIB and LW individuals, enabling the identification of a comprehensive set of genetic variations for each breed. The specific and predominant SNPs in the TIB and LW populations were detected by using a cutoff value of 0.50 for SNP allele frequency and absolute allele frequency differences (△AF) between the TIB and LW populations. RESULTS: A total of 21,767,938 SNPs were retrieved from 44 TIB and 29 LW genomes. The analysis detected 2,893,106 (13.29%) and 813,310 (3.74%) specific and predominant SNPs in the TIB and LW populations, and annotated to 24,560 genes. Further GO analysis revealed 291 genes involved in biological processes related to striated and/or skeletal muscle differentiation, proliferation, hypertrophy, regulation of striated muscle cell differentiation and proliferation, and myoblast differentiation and fusion. These 291 genes included crucial regulators of muscle cell determination, proliferation, differentiation, and hypertrophy, such as members of the Myogenic regulatory factors (MRF) (MYOD, MYF5, MYOG, MYF6) and Myocyte enhancer factor 2 (MEF2) (MEF2A, MEF2C, MEF2D) families, as well as muscle growth inhibitors (MSTN, ACVR1, and SMAD1); KEGG pathway analysis revealed 106 and 20 genes were found in muscle growth related positive and negative regulatory signaling pathways. Notably, genes critical for protein synthesis, such as MTOR, IGF1, IGF1R, IRS1, INSR, and RPS6KA6, were implicated in these pathways. CONCLUSION: This study employed an effective methodology to rigorously identify the potential genes associated with skeletal muscle development. A substantial number of SNPs and genes that potentially play roles in the divergence observed in skeletal muscle growth between the TIB and LW breeds were identified. These findings offer valuable insights into the genetic underpinnings of skeletal muscle development and present opportunities for enhancing meat production through pig breeding.

Haplotype-resolved assembly of a pig genome using single-sperm sequencing.

Single gamete cell sequencing together with long-read sequencing can reliably produce chromosome-level phased genomes. In this study, we employed PacBio HiFi and Hi-C sequencing on a male Landrace pig, coupled with single-sperm sequencing of its 102 sperm cells. A haplotype assembly method was developed based on long-read sequencing and sperm-phased markers. The chromosome-level phased assembly showed higher phasing accuracy than methods that rely only on HiFi reads. The use of single-sperm sequencing data enabled the construction of a genetic map, successfully mapping the sperm motility trait to a specific region on chromosome 1 (105.40-110.70 Mb). Furthermore, with the assistance of Y chromosome-bearing sperm data, 26.16 Mb Y chromosome sequences were assembled. We report a reliable approach for assembling chromosome-level phased genomes and reveal the potential of sperm population in basic biology research and sperm phenotype research.

Exploring the Interplay between the Hologenome and Complex Traits in Bovine and Porcine Animals Using Genome-Wide Association Analysis.

Genome-wide association studies (GWAS) significantly enhance our ability to identify trait-associated genomic variants by considering the host genome. Moreover, the hologenome refers to the host organism's collective genetic material and its associated microbiome. In this study, we utilized the hologenome framework, called Hologenome-wide association studies (HWAS), to dissect the architecture of complex traits, including milk yield, methane emissions, rumen physiology in cattle, and gut microbial composition in pigs. We employed four statistical models: (1) GWAS, (2) Microbial GWAS (M-GWAS), (3) HWAS-CG (hologenome interaction estimated using COvariance between Random Effects Genome-based restricted maximum likelihood (CORE-GREML)), and (4) HWAS-H (hologenome interaction estimated using the Hadamard product method). We applied Bonferroni correction to interpret the significant associations in the complex traits. The GWAS and M-GWAS detected one and sixteen significant SNPs for milk yield traits, respectively, whereas the HWAS-CG and HWAS-H each identified eight SNPs. Moreover, HWAS-CG revealed four, and the remaining models identified three SNPs each for methane emissions traits. The GWAS and HWAS-CG detected one and three SNPs for rumen physiology traits, respectively. For the pigs' gut microbial composition traits, the GWAS, M-GWAS, HWAS-CG, and HWAS-H identified 14, 16, 13, and 12 SNPs, respectively. We further explored these associations through SNP annotation and by analyzing biological processes and functional pathways. Additionally, we integrated our GWA results with expression quantitative trait locus (eQTL) data using transcriptome-wide association studies (TWAS) and summary-based Mendelian randomization (SMR) methods for a more comprehensive understanding of SNP-trait associations. Our study revealed hologenomic variability in agriculturally important traits, enhancing our understanding of host-microbiome interactions.

Predictive ability of multi-population genomic prediction methods of phenotypes for reproduction traits in Chinese and Austrian pigs.

BACKGROUND: Multi-population genomic prediction can rapidly expand the size of the reference population and improve genomic prediction ability. Machine learning (ML) algorithms have shown advantages in single-population genomic prediction of phenotypes. However, few studies have explored the effectiveness of ML methods for multi-population genomic prediction. RESULTS: In this study, 3720 Yorkshire pigs from Austria and four breeding farms in China were used, and single-trait genomic best linear unbiased prediction (ST-GBLUP), multitrait GBLUP (MT-GBLUP), Bayesian Horseshoe (BayesHE), and three ML methods (support vector regression (SVR), kernel ridge regression (KRR) and AdaBoost.R2) were compared to explore the optimal method for joint genomic prediction of phenotypes of Chinese and Austrian pigs through 10 replicates of fivefold cross-validation. In this study, we tested the performance of different methods in two scenarios: (i) including only one Austrian population and one Chinese pig population that were genetically linked based on principal component analysis (PCA) (designated as the "two-population scenario") and (ii) adding reference populations that are unrelated based on PCA to the above two populations (designated as the "multi-population scenario"). Our results show that, the use of MT-GBLUP in the two-population scenario resulted in an improvement of 7.1% in predictive ability compared to ST-GBLUP, while the use of SVR and KKR yielded improvements in predictive ability of 4.5 and 5.3%, respectively, compared to MT-GBLUP. SVR and KRR also yielded lower mean square errors (MSE) in most population and trait combinations. In the multi-population scenario, improvements in predictive ability of 29.7, 24.4 and 11.1% were obtained compared to ST-GBLUP when using, respectively, SVR, KRR, and AdaBoost.R2. However, compared to MT-GBLUP, the potential of ML methods to improve predictive ability was not demonstrated. CONCLUSIONS: Our study demonstrates that ML algorithms can achieve better prediction performance than multitrait GBLUP models in multi-population genomic prediction of phenotypes when the populations have similar genetic backgrounds; however, when reference populations that are unrelated based on PCA are added, the ML methods did not show a benefit. When the number of populations increased, only MT-GBLUP improved predictive ability in both validation populations, while the other methods showed improvement in only one population.

Mapping and functional characterization of structural variation in 1060 pig genomes.

BACKGROUND: Structural variations (SVs) have significant impacts on complex phenotypes by rearranging large amounts of DNA sequence. RESULTS: We present a comprehensive SV catalog based on the whole-genome sequence of 1060 pigs (Sus scrofa) representing 101 breeds, covering 9.6% of the pig genome. This catalog includes 42,487 deletions, 37,913 mobile element insertions, 3308 duplications, 1664 inversions, and 45,184 break ends. Estimates of breed ancestry and hybridization using genotyped SVs align well with those from single nucleotide polymorphisms. Geographically stratified deletions are observed, along with known duplications of the KIT gene, responsible for white coat color in European pigs. Additionally, we identify a recent SINE element insertion in MYO5A transcripts of European pigs, potentially influencing alternative splicing patterns and coat color alterations. Furthermore, a Yorkshire-specific copy number gain within ABCG2 is found, impacting chromatin interactions and gene expression across multiple tissues over a stretch of genomic region of ~200 kb. Preliminary investigations into SV's impact on gene expression and traits using the Pig Genotype-Tissue Expression (PigGTEx) data reveal SV associations with regulatory variants and gene-trait pairs. For instance, a 51-bp deletion is linked to the lead eQTL of the lipid metabolism regulating gene FADS3, whose expression in embryo may affect loin muscle area, as revealed by our transcriptome-wide association studies. CONCLUSIONS: This SV catalog serves as a valuable resource for studying diversity, evolutionary history, and functional shaping of the pig genome by processes like domestication, trait-based breeding, and adaptive evolution.

Novel off-Targeting Events Identified after Genome Wide Analysis of CRISPR-Cas Edited Pigs.

CRISPR-Cas technology has transformed our ability to introduce targeted modifications, allowing unconventional animal models such as pigs to model human diseases and improve its value for food production. The main concern with using the technology is the possibility of introducing unwanted modifications in the genome. In this study, we illustrate a pipeline to comprehensively identify off-targeting events on a global scale in the genome of three different gene-edited pig models. Whole genome sequencing paired with an off-targeting prediction software tool filtered off-targeting events amongst natural variations present in gene-edited pigs. This pipeline confirmed two known off-targeting events in IGH knockout pigs, AR and RBFOX1, and identified other presumably off-targeted loci. Independent validation of the off-targeting events using other gene-edited DNA confirmed two novel off-targeting events in RAG2/IL2RG knockout pig models. This unique strategy offers a novel tool to detect off-targeting events in genetically heterogeneous species after genome editing.

Residual network improves the prediction accuracy of genomic selection.

Genetic improvement of complex traits in animal and plant breeding depends on the efficient and accurate estimation of breeding values. Deep learning methods have been shown to be not superior over traditional genomic selection (GS) methods, partially due to the degradation problem (i.e. with the increase of the model depth, the performance of the deeper model deteriorates). Since the deep learning method residual network (ResNet) is designed to solve gradient degradation, we examined its performance and factors related to its prediction accuracy in GS. Here we compared the prediction accuracy of conventional genomic best linear unbiased prediction, Bayesian methods (BayesA, BayesB, BayesC, and Bayesian Lasso), and two deep learning methods, convolutional neural network and ResNet, on three datasets (wheat, simulated and real pig data). ResNet outperformed other methods in both Pearson's correlation coefficient (PCC) and mean squared error (MSE) on the wheat and simulated data. For the pig backfat depth trait, ResNet still had the lowest MSE, whereas Bayesian Lasso had the highest PCC. We further clustered the pig data into four groups and, on one separated group, ResNet had the highest prediction accuracy (both PCC and MSE). Transfer learning was adopted and capable of enhancing the performance of both convolutional neural network and ResNet. Taken together, our findings indicate that ResNet could improve GS prediction accuracy, affected potentially by factors such as the genetic architecture of complex traits, data volume, and heterogeneity.

Study on the Characteristics of Coarse Feeding Tolerance of Ding'an Pigs: Phenotypic and Candidate Genes Identification.

Ding'an (DA) pig, a prominent local breed in Hainan Province, exhibits notable advantages in coarse feeding tolerance and high-quality meat. To explore the potential genetic mechanism of coarse feeding tolerance in DA pigs, 60-day-old full sibling pairs of DA and DLY (Duroc-Landrace-Yorkshire) pigs were subjected to fed normal (5%) and high (10%) crude fiber diets for 56 days, respectively. The findings showed that increasing the crude fiber level had no impact on the apparent digestibility of crude fiber, intramuscular fat, and marbling scores in DA pigs, whereas these factors were significantly reduced in DLY pigs (p < 0.05). Through differential expression analysis and Weighted Gene Co-expression Network Analysis (WGCNA) of the colonic mucosal transcriptome data, 65 and 482 candidate genes with coarse feeding tolerance in DA pigs were identified, respectively. Joint analysis screened four key candidate genes, including LDHB, MLC1, LSG1, and ESM1, potentially serving as key regulated genes for coarse feeding tolerance. Functional analysis revealed that the most significant pathway enriched in differential genes associated with coarse feeding tolerance in Ding'an pigs was the signaling receptor binding. The results hold substantial significance for advancing our understanding of the genetic mechanisms governing coarse feeding tolerance in Ding'an pigs.