Search

Background: Backfat serves as a vital fat reservoir in pigs, and its excessive accumulation will adversely impact pig growth performance, farming efficiency, and pork quality. The aim of this research is to integrate assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) and RNA sequencing (RNA-seq) to explore the molecular mechanisms underlying porcine backfat deposition., Results: ATAC-seq analysis identified 568 genes originating from 698 regions exhibiting differential accessibility, which were significantly enriched in pathways pertinent to adipocyte differentiation and lipid metabolism. Besides, a total of 283 transcription factors (TFs) were identified by motif analysis. RNA-seq analysis revealed 978 differentially expressed genes (DEGs), which were enriched in pathways related to energy metabolism, cell cycle and signal transduction. The integration of ATAC-seq and RNA-seq data indicates that DEG expression levels are associated with chromatin accessibility. This comprehensive study highlights the involvement of critical pathways, including the Wnt signaling pathway, Jak-STAT signaling pathway, and fatty acid degradation, in the regulation of backfat deposition. Through rigorous analysis, we identified several candidate genes (LEP, CTBP2, EHHADH, OSMR, TCF7L2, BCL2, FGF1, UCP2, CCND1, TIMP1, and VDR) as potentially significant contributors to backfat deposition. Additionally, we constructed TF-TF and TF-target gene regulatory networks and identified a series of potential TFs related to backfat deposition (FOS, STAT3, SMAD3, and ESR1)., Conclusions: This study represents the first application of ATAC-seq and RNA-seq, affording a novel perspective into the mechanisms underlying backfat deposition and providing invaluable resources for the enhancement of pig breeding programs., (© 2024. The Author(s).)

Background: Teat number is a vital reproductive trait in sows, crucial for providing immunity and nutrition to piglets during lactation. However, "missing heritability" in Single Nucleotide Polymorphism (SNP)-based Genome-Wide Association Studies (GWAS) has led to an increasing focus on structural variations in the genetic analysis of complex biological traits., Results: In this study, we generated a comprehensive CNV map in a population of French Yorkshire pigs (n = 644) and identified 429 CNVRs. Notably, 44% (189 CNVRs) of these were detected for the first time. Subsequently, we conducted GWAS for teat number in the French Yorkshire pig population using both 80K chip and its imputed data, as well as a GWAS analysis based on CNV regions (CNVR). Interestingly, 80K chip GWAS identified two SNPs located on Sus scrofa chromosome 5 (SSC5) that were simultaneously associated with Total Teat Number (TTN), Left Teat Number (LTN), and Right Teat Number (RTN). The leading SNP (WU_10.2_5_76130558) explained 3.33%, 2.69%, and 2.67% of the phenotypic variance for TTN, LTN, and RTN, respectively. Moreover, through imputed GWAS, we successfully identified 30 genetic variants associated with TTN located within the 73.22 -73.30 Mb region on SSC5. The two SNPs identified in the 80K chip GWAS were also located in this region. In addition, CNVR-based GWAS revealed three significant CNVRs associated with TTN. Finally, through gene annotation, we pinpointed two candidate genes, TRIM66 and PRICKLE1, which are related to diverse processes such as breast cancer and abnormal vertebral development., Conclusions: Our research provides an in-depth analysis of the complex genetic structure underlying teat number, contributing to the genetic enhancement of sows with improved reproductive performance and, ultimately, bolstering the economic benefits of swine production enterprises., (© 2024. The Author(s).)

Background: Integration of high throughput DNA genotyping and RNA-sequencing data enables the discovery of genomic regions that regulate gene expression, known as expression quantitative trait loci (eQTL). In pigs, efforts to date have been mainly focused on purebred lines for traits with commercial relevance as such growth and meat quality. However, little is known on genetic variants and mechanisms associated with the robustness of an animal, thus its overall health status. Here, the liver, lung, spleen, and muscle transcriptomes of 100 three-way crossbred female finishers were studied, with the aim of identifying novel eQTL regulatory regions and transcription factors (TFs) associated with regulation of porcine metabolism and health-related traits., Results: An expression genome-wide association study with 535,896 genotypes and the expression of 12,680 genes in liver, 13,310 genes in lung, 12,650 genes in spleen, and 12,595 genes in muscle resulted in 4,293, 10,630, 4,533, and 6,871 eQTL regions for each of these tissues, respectively. Although only a small fraction of the eQTLs were annotated as cis-eQTLs, these presented a higher number of polymorphisms per region and significantly stronger associations with their target gene compared to trans-eQTLs. Between 20 and 115 eQTL hotspots were identified across the four tissues. Interestingly, these were all enriched for immune-related biological processes. In spleen, two TFs were identified: ERF and ZNF45, with key roles in regulation of gene expression., Conclusions: This study provides a comprehensive analysis with more than 26,000 eQTL regions identified that are now publicly available. The genomic regions and their variants were mostly associated with tissue-specific regulatory roles. However, some shared regions provide new insights into the complex regulation of genes and their interactions that are involved with important traits related to metabolism and immunity., (© 2024. The Author(s).)

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., (© 2024. The Author(s).)

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., (© 2024. The Author(s).)

Background: Skeletal muscle development plays a crucial role in yield and quality of pork; however, this process is influenced by various factors. In this study, we employed whole-genome bisulfite sequencing (WGBS) and transcriptome sequencing to comprehensively investigate the longissimus dorsi muscle (LDM), aiming to identify key genes that impact the growth and development of Duroc pigs with different average daily gains (ADGs)., Results: Eight pigs were selected and divided into two groups based on ADGs: H (774.89 g) group and L (658.77 g) group. Each pair of the H and L groups were half-siblings. The results of methylation sequencing revealed 2631 differentially methylated genes (DMGs) involved in metabolic processes, signalling, insulin secretion, and other biological activities. Furthermore, a joint analysis was conducted on these DMGs and the differentially expressed genes (DEGs) obtained from transcriptome sequencing of the same individual. This analysis identified 316 differentially methylated and differentially expressed genes (DMEGs), including 18 DMEGs in promoter regions and 294 DMEGs in gene body regions. Finally, LPAR1 and MEF2C were selected as candidate genes associated with muscle development. Bisulfite sequencing PCR (BSP) and quantitative real-time PCR (qRT-PCR) revealed that the promoter region of LPAR1 exhibited significantly lower methylation levels (P < 0.05) and greater expression levels (P < 0.05) in the H group than in the L group. Additionally, hypermethylation was observed in the gene body region of MEF2C, as was a low expression level, in the H group (P < 0.05)., Conclusions: These results suggest that the differences in the ADGs of Duroc pigs fed the same diet may be influenced by the methylation levels and expression levels of genes related to skeletal muscle development., (© 2024. The Author(s).)

Background: The identification of low-frequency haplotypes, never observed in homozygous state in a population, is considered informative on the presence of potentially harmful alleles (candidate alleles), putatively involved in inbreeding depression. Although identification of candidate alleles is challenging, studies analyzing the dynamics of potentially harmful alleles are lacking. A pedigree of the highly endangered Gochu Asturcelta pig breed, including 471 individuals belonging to 51 different families with at least 5 offspring each, was genotyped using the Axiom PigHDv1 Array (658,692 SNPs). Analyses were carried out on four different cohorts defined according to pedigree depth and at the whole population (WP) level., Results: The 4,470 Linkage Blocks (LB) identified in the Base Population (10 individuals), gathered a total of 16,981 alleles in the WP. Up to 5,466 (32%) haplotypes were statistically considered candidate alleles, 3,995 of them (73%) having one copy only. The number of alleles and candidate alleles varied across cohorts according to sample size. Up to 4,610 of the alleles identified in the WP (27% of the total) were present in one cohort only. Parentage analysis identified a total of 67,742 parent-offspring incompatibilities. The number of mismatches varied according to family size. Parent-offspring inconsistencies were identified in 98.2% of the candidate alleles and 100% of the LB in which they were located. Segregation analyses informed that most potential candidate alleles appeared de novo in the pedigree. Only 17 candidate alleles were identified in the boar, sow, and paternal and maternal grandparents and were considered segregants., Conclusions: Our results suggest that neither mutation nor recombination are the major forces causing the apparition of candidate alleles. Their occurrence is more likely caused by Allele-Drop-In events due to SNP calling errors. New alleles appear when wrongly called SNPs are used to construct haplotypes. The presence of candidate alleles in either parents or grandparents of the carrier individuals does not ensure that they are true alleles. Minimum Allele Frequency thresholds may remove informative alleles. Only fully segregant candidate alleles should be considered potentially harmful alleles. A set of 16 candidate genes, potentially involved in inbreeding depression, is described., (© 2024. The Author(s).)

Background: Heat stress (HS) poses significant threats to the sustainability of livestock production. Genetically improving heat tolerance could enhance animal welfare and minimize production losses during HS events. Measuring phenotypic indicators of HS response and understanding their genetic background are crucial steps to optimize breeding schemes for improved climatic resilience. The identification of genomic regions and candidate genes influencing the traits of interest, including variants with pleiotropic effects, enables the refinement of genotyping panels used to perform genomic prediction of breeding values and contributes to unraveling the biological mechanisms influencing heat stress response. Therefore, the main objectives of this study were to identify genomic regions, candidate genes, and potential pleiotropic variants significantly associated with indicators of HS response in lactating sows using imputed whole-genome sequence (WGS) data. Phenotypic records for 18 traits and genomic information from 1,645 lactating sows were available for the study. The genotypes from the PorcineSNP50K panel containing 50,703 single nucleotide polymorphisms (SNPs) were imputed to WGS and after quality control, 1,622 animals and 7,065,922 SNPs were included in the analyses., Results: A total of 1,388 unique SNPs located on sixteen chromosomes were found to be associated with 11 traits. Twenty gene ontology terms and 11 biological pathways were shown to be associated with variability in ear skin temperature, shoulder skin temperature, rump skin temperature, tail skin temperature, respiration rate, panting score, vaginal temperature automatically measured every 10 min, vaginal temperature measured at 0800 h, hair density score, body condition score, and ear area. Seven, five, six, two, seven, 15, and 14 genes with potential pleiotropic effects were identified for indicators of skin temperature, vaginal temperature, animal temperature, respiration rate, thermoregulatory traits, anatomical traits, and all traits, respectively., Conclusions: Physiological and anatomical indicators of HS response in lactating sows are heritable but highly polygenic. The candidate genes found are associated with important gene ontology terms and biological pathways related to heat shock protein activities, immune response, and cellular oxidative stress. Many of the candidate genes with pleiotropic effects are involved in catalytic activities to reduce cell damage from oxidative stress and cellular mechanisms related to immune response., (© 2024. The Author(s).)

Background: The health and size of the testes are crucial for boar fertility. Testicular development is tightly regulated by epigenetics. N6-methyladenosine (m6A) modification is a prevalent internal modification on mRNA and plays an important role in development. The mRNA m6A methylation in boar testicular development still needs to be investigated., Results: Using the MeRIP-seq technique, we identify and profile m6A modification in boar testes between piglets and adults. The results showed 7783 distinct m6A peaks in piglets and 6590 distinct m6A peaks in adults, with 2,471 peaks shared between the two groups. Enrichment of GO and KEGG analysis reveal dynamic m6A methylation in various biological processes and signalling pathways. Meanwhile, we conjointly analyzed differentially methylated and expressed genes in boar testes before and after sexual maturity, and reproductive related genes (TLE4, TSSK3, TSSK6, C11ORF94, PATZ1, PHLPP1 and PAQR7) were identified. Functional enrichment analysis showed that differential genes are associated with important biological functions, including regulation of growth and development, regulation of metabolic processes and protein catabolic processes., Conclusion: The results demonstrate that m6A methylation, differential expression and the related signalling pathways are crucial for boar testicular development. These results suggest a role for m6A modification in boar testicular development and provided a resource for future studies on m6A function in boar testicular development., (© 2024. The Author(s).)

Background: The pig (Sus Scrofa) is one of the oldest domesticated livestock species that has undergone extensive improvement through modern breeding. European breeds have advantages in lean meat development and highly-productive body type, whereas Asian breeds possess extraordinary fat deposition and reproductive performance. Consequently, Eurasian breeds have been extensively used to develop modern commercial breeds for fast-growing and high prolificacy. However, limited by the sequencing technology, the genome architecture of some nascent developed breeds and the human-mediated impact on their genomes are still unknown., Results: Through whole-genome analysis of 178 individuals from an Asian locally developed pig breed, Beijing Black pig, and its two ancestors from two different continents, we found the pervasive inconsistent gene trees and species trees across the genome of Beijing Black pig, which suggests its introgressive hybrid origin. Interestingly, we discovered that this developed breed has more genetic relationships with European pigs and an unexpected introgression from Asian pigs to this breed, which indicated that human-mediated introgression could form the porcine genome architecture in a completely different type compared to native introgression. We identified 554 genomic regions occupied 63.30 Mb with signals of introgression from the Asian ancestry to Beijing Black pig, and the genes in these regions enriched in pathways associated with meat quality, fertility, and disease-resistant. Additionally, a proportion of 7.77% of genomic regions were recognized as regions that have been under selection. Moreover, combined with the results of a genome-wide association study for meat quality traits in the 1537 Beijing Black pig population, two important candidate genes related to meat quality traits were identified. DNAJC6 is related to intramuscular fat content and fat deposition, and RUFY4 is related to meat pH and tenderness., Conclusions: Our research provides insight for analyzing the origins of nascent developed breeds and genome-wide selection remaining in the developed breeds mediated by humans during modern breeding., (© 2024. The Author(s).)

Background: Many studies have been performed to identify various genomic loci and genes associated with the meat quality in pigs. However, the full genetic architecture of the trait still remains unclear in part because of the lack of accurate identification of related structural variations (SVs) which resulted from the shortage of target breeds, the limitations of sequencing data, and the incompleteness of genome assemblies. The recent generation of a new pig breed with superior meat quality, called Nanchukmacdon, and its chromosome-level genome assembly (the NCMD assembly) has provided new opportunities., Results: By applying assembly-based SV calling approaches to various genome assemblies of pigs including Nanchukmacdon, the impact of SVs on meat quality was investigated. Especially, by checking the commonality of SVs with other pig breeds, a total of 13,819 Nanchukmacdon-specific SVs (NSVs) were identified, which have a potential effect on the unique meat quality of Nanchukmacdon. The regulatory potentials of NSVs for the expression of nearby genes were further examined using transcriptome- and epigenome-based analyses in different tissues., Conclusions: Whole-genome comparisons based on chromosome-level genome assemblies have led to the discovery of SVs affecting meat quality in pigs, and their regulatory potentials were analyzed. The identified NSVs will provide new insights regarding genetic architectures underlying the meat quality in pigs. Finally, this study confirms the utility of chromosome-level genome assemblies and multi-omics analysis to enhance the understanding of unique phenotypes., (© 2024. The Author(s).)

Background: Muscle growth post-birth relies on muscle fiber number and size. Myofibre number, metabolic and contractile capacities are established pre-birth during prenatal myogenesis. The aim of this study was to identify genes involved in skeletal muscle development in cattle, sheep, and pigs - livestock., Results: The cattle analysis showed significant differences in 5043 genes during the 135-280 dpc period. In sheep, 444 genes differed significantly during the 70-120 dpc period. Pigs had 905 significantly different genes for the 63-91 dpc period.The biological processes and KEGG pathway enrichment results in each species individually indicated that DEGs in cattle were significantly enriched in regulation of cell proliferation, cell division, focal adhesion, ECM-receptor interaction, and signaling pathways (PI3K-Akt, PPAR, MAPK, AMPK, Ras, Rap1); in sheep - positive regulation of fibroblast proliferation, negative regulation of endothelial cell proliferation, focal adhesion, ECM-receptor interaction, insulin resistance, and signaling pathways (PI3K-Akt, HIF-1, prolactin, Rap1, PPAR); in pigs - regulation of striated muscle tissue development, collagen fibril organization, positive regulation of insulin secretion, focal adhesion, ECM-receptor interaction, and signaling pathways (PPAR, FoxO, HIF-1, AMPK). Among the DEGs common for studied animal species, 45 common genes were identified. Based on these, a protein-protein interaction network was created and three significant modules critical for skeletal muscle myogenesis were found, with the most significant module A containing four recognized hub genes - EGFR, VEGFA, CDH1, and CAV1. Using the miRWALK and TF2DNA databases, miRNAs (bta-miR-2374 and bta-miR-744) and transcription factors (CEBPB, KLF15, RELA, ZNF143, ZBTB48, and REST) associated with hub genes were detected. Analysis of GO term and KEGG pathways showed that such processes are related to myogenesis and associated with module A: positive regulation of MAP kinase activity, vascular endothelial growth factor receptor, insulin-like growth factor binding, focal adhesion, and signaling pathways (PI3K-Akt, HIF-1, Rap1, Ras, MAPK)., Conclusions: The identified genes, common to the prenatal developmental period of skeletal muscle in livestock, are critical for later muscle development, including its growth by hypertrophy. They regulate valuable economic characteristics. Enhancing and breeding animals according to the recognized genes seems essential for breeders to achieve superior gains in high-quality muscle mass., (© 2024. The Author(s).)

Background: Cultured porcine cell lines are powerful tools for functional genomics and in vitro phenotypic testing of candidate causal variants. However, to be utilised for genomic or variant interrogation assays, the genome sequence and structure of cultured cell lines must be realised. In this work, we called variants and used read coverage in combination with within-sample allele frequency to detect potential aneuploidy in two immortalised porcine kidney epithelial (PK15) cell lines and in a pig embryonic fibroblast line., Results: We compared two PK15 cultured cells samples: a new American Type Culture Collection (ATCC) sample and one that has been utilised and passaged within the laboratory for an extended period (> 10 years). Read coverage and within-sample allele frequencies showed that several chromosomes are fully or partially aneuploid in both PK15 lines, including potential trisomy of chromosome 4 and tetrasomy of chromosome 17. The older PK15 line showed evidence of additional structural variation and potentially clonal variation. By comparison, the pig embryonic fibroblast line was free from the gross aneuploidies seen in the PK15s., Conclusions: Our results show that the PK15 cell lines examined have aneuploidies and complex structural variants in their genomes. We propose that screening for aneuploidy should be considered for cell lines, and discuss implications for livestock genomics., (© 2024. The Author(s).)

Background: Shaziling pig is a well-known indigenous breed in China who has superior meat quality traits. However, the genetic mechanism and genomic evidence underlying meat quality characteristics of Shaziling pigs are still unclear. To explore and investigate the germplasm characteristics of Shaziling pigs, we totally analyzed 67 individual's whole genome sequencing data for the first time (20 Shaziling pigs [S], 20 Dabasha pigs [DBS], 11 Yorkshire pigs [Y], 10 Berkshire pigs [BKX], 5 Basha pigs [BS] and 1 Warthog)., Results: A total of 2,538,577 SNPs with high quality were detected and 9 candidate genes which was specifically selected in S and shared in S to DBS were precisely mined and screened using an integrated analysis strategy of identity-by-descent (IBD) and selective sweep. Of them, dickkopf WNT signaling pathway inhibitor 2 (DKK2), the antagonist of Wnt signaling pathway, was the most promising candidate gene which was not only identified an association of palmitic acid and palmitoleic acid quantitative trait locus in PigQTLdb, but also specifically selected in S compared to other 48 Chinese local pigs of 12 populations and 39 foreign pigs of 4 populations. Subsequently, a mutation at 12,726-bp of DKK2 intron 1 (g.114874954 A > C) was identified associated with intramuscular fat content using method of PCR-RFLP in 21 different pig populations. We observed DKK2 specifically expressed in adipose tissues. Overexpression of DKK2 decreased the content of triglyceride, fatty acid synthase and expression of relevant genes of adipogenic and Wnt signaling pathway, while interference of DKK2 got contrary effect during adipogenesis differentiation of porcine preadipocytes and 3T3-L1 cells., Conclusions: Our findings provide an analysis strategy for mining functional genes of important economic traits and provide fundamental data and molecular evidence for improving pig meat quality traits and molecular breeding., (© 2023. The Author(s).)

Background: Mapping expression quantitative trait loci (eQTLs) in skeletal muscle tissue in pigs is crucial for understanding the relationship between genetic variation and phenotypic expression of carcass traits in meat animals. Therefore, the primary objective of this study was to evaluate the impact of different sets of single nucleotide polymorphisms (SNP), including scenarios removing SNPs pruned for linkage disequilibrium (LD) and SNPs derived from SNP chip arrays and RNA-seq data from liver, brain, and skeletal muscle tissues, on the identification of eQTLs in the Longissimus lumborum tissue, associated with carcass and body composition traits in Large White pigs. The SNPs identified from muscle mRNA were combined with SNPs identified in the brain and liver tissue transcriptomes, as well as SNPs from the GGP Porcine 50 K SNP chip array. Cis- and trans-eQTLs were identified based on the skeletal muscle gene expression level, followed by functional genomic analyses and statistical associations with carcass and body composition traits in Large White pigs., Results: The number of cis- and trans-eQTLs identified across different sets of SNPs (scenarios) ranged from 261 to 2,539 and from 29 to 13,721, respectively. Furthermore, 6,180 genes were modulated by eQTLs in at least one of the scenarios evaluated. The eQTLs identified were not significantly associated with carcass and body composition traits but were significantly enriched for many traits in the "Meat and Carcass" type QTL. The scenarios with the highest number of cis- (n = 304) and trans- (n = 5,993) modulated genes were the unpruned and LD-pruned SNP set scenarios identified from the muscle transcriptome. These genes include 84 transcription factor coding genes., Conclusions: After LD pruning, the set of SNPs identified based on the transcriptome of the skeletal muscle tissue of pigs resulted in the highest number of genes modulated by eQTLs. Most eQTLs are of the trans type and are associated with genes influencing complex traits in pigs, such as transcription factors and enhancers. Furthermore, the incorporation of SNPs from other genomic regions to the set of SNPs identified in the porcine skeletal muscle transcriptome contributed to the identification of eQTLs that had not been identified based on the porcine skeletal muscle transcriptome alone., (© 2023. The Author(s).)

Background: Jianli pig, a renowned indigenous breed in China, has the characteristics of a two-end black (TEB) coat color, excellent meat quality, strong adaptability and increased prolificacy. However, there is limited information available regarding the genetic diversity, population structure and genomic regions under selection of Jianli pig. On the other hand, the genetic mechanism of TEB coat color has remained largely unknown., Results: In this study, the whole genome resequencing of 30 Jianli pigs within a context of 153 individuals representing 13 diverse breeds was performed. The population structure analysis revealed that Jianli pigs have close genetic relationships with the Tongcheng pig breed, their geographical neighbors. Three methods (observed heterozygosity, expected heterozygosity, and runs of homozygosity) implied a relatively high level of genetic diversity and, a low inbreeding coefficient in Jianli compared with other pigs. We used Fst and XP-EHH to detect the selection signatures in Jianli pigs compared with Asian wild boar. A total of 451 candidate genes influencing meat quality (CREBBP, ADCY9, EEPD1 and HDAC9), reproduction (ESR1 and FANCA), and coat color (EDNRB, MITF and MC1R), were detected by gene annotation analysis. Finally, to fine-map the genomic region for the two-end black (TEB) coat color phenotype in Jianli pigs, we performed three signature selection methods between the TEB coat color and no-TEB coat color pig breeds. The current study, further confirmed that the EDNRB gene is a candidate gene for TEB color phenotype found in Chinese pigs, including Jinhua pigs, and the haplotype harboring 25 SNPs in the EDNRB gene may promote the formation of TEB coat color. Further ATAC-seq and luciferase reporter assays of these regions suggest that the 25-SNPs region was a strong candidate causative mutation that regulates the TEB coat color phenotype by altering enhancer function., Conclusion: Our results advanced the understanding of the genetic mechanism behind artificial selection, and provided further resources for the protection and breeding improvement of Jianli pigs., (© 2024. The Author(s).)

Background: Chinese indigenous pigs in Yunnan exhibit considerable phenotypic diversity, but their population structure and the biological interpretation of signatures of artificial selection require further investigation. To uncover population genetic diversity, migration events, and artificial selection signatures in Chinese domestic pigs, we sampled 111 Yunnan pigs from four breeds in Yunnan which is considered to be one of the centres of livestock domestication in China, and genotyped them using Illumina Porcine SNP60K BeadChip. We then leveraged multiple bioinformatics database tools to further investigate the signatures and associated complex traits., Results: Population structure and migration analyses showed that Diannanxiaoer pigs had different genetic backgrounds from other Yunnan pigs, and Gaoligongshan may undergone the migration events from Baoshan and Saba pigs. Intriguingly, we identified a possible common target of sharing artificial selection on a 265.09 kb region on chromosome 5 in Yunnan indigenous pigs, and the genes on this region were associated with cardiovascular and immune systems. We also detected several candidate genes correlated with dietary adaptation, body size (e.g., PASCIN1, GRM4, ITPR2), and reproductive performance. In addition, the breed-sharing gene MMP16 was identified to be a human-mediated gene. Multiple lines of evidence at the mammalian genome, transcriptome, and phenome levels further supported the evidence for the causality between MMP16 variants and the metabolic diseases, brain development, and cartilage tissues in Chinese pigs. Our results suggested that the suppression of MMP16 would directly lead to inactivity and insensitivity of neuronal activity and skeletal development in Chinese indigenous pigs., Conclusion: In this study, the population genetic analyses and identification of artificial selection signatures of Yunnan indigenous pigs help to build an understanding of the effect of human-mediated selection mechanisms on phenotypic traits in Chinese indigenous pigs. Further studies are needed to fully characterize the process of human-mediated genes and biological mechanisms., (© 2023. The Author(s).)

Background: Eurasian pigs have undergone lineage admixture throughout history. It has been confirmed that the genes of indigenous pig breeds in China have been introduced into Western commercial pigs, providing genetic materials for breeding Western pigs. Pigs in Taihu Lake region (TL), such as the Meishan pig and Erhualian pig, serve as typical representatives of indigenous pig breeds in China due to their high reproductive performances. These pigs have also been imported into European countries in 1970 and 1980 s. They have played a positive role in improving the reproductive performances in European commercial pigs such as French Large White pigs (FLW). However, it is currently unclear if the lineage of TL pigs have been introgressed into the Danish Large White pigs (DLW), which are also known for their high reproductive performances in European pigs. To systematically identify genomic regions in which TL pigs have introgressed into DLW pigs and their physiological functions, we collected the re-sequencing data from 304 Eurasian pigs, to identify shared haplotypes between DLW and TL pigs., Results: The findings revealed the presence of introgressed genomic regions from TL pigs in the genome of DLW pigs indeed. The genes annotated within these regions were found to be mainly enriched in neurodevelopmental pathways. Furthermore, we found that the 115 kb region located in SSC16 exhibited highly shared haplotypes between TL and DLW pigs. The major haplotype of TL pigs in this region could significantly improve reproductive performances in various pig populations. Around this genomic region, NDUFS4 gene was highly expressed and showed differential expression in multiple reproductive tissues between extremely high and low farrowing Erhualian pigs. This suggested that NDUFS4 gene could be an important candidate causal gene responsible for affecting the reproductive performances of DLW pigs., Conclusions: Our study has furthered our knowledge of the pattern of introgression from TL into DLW pigs and the potential effects on the fertility of DLW pigs., (© 2023. The Author(s).)

Background: Feed efficiency is a research priority to support a sustainable meat production. It is recognized as a complex trait that integrates multiple biological pathways orchestrated in and by various tissues. This study aims to determine networks between biological entities to explain inter-individual variation of feed efficiency in growing pigs., Results: The feed conversion ratio (FCR), a measure of feed efficiency, and its two component traits, average daily gain and average daily feed intake, were obtained from 47 growing pigs from a divergent selection for residual feed intake and fed high-starch or high-fat high-fiber diets during 58 days. Datasets of transcriptomics (60 k porcine microarray) in the whole blood and metabolomics (1H-NMR analysis and target gas chromatography) in plasma were available for all pigs at the end of the trial. A weighted gene co-expression network was built from the transcriptomics dataset, resulting in 33 modules of co-expressed molecular probes. The eigengenes of eight of these modules were significantly ([Formula: see text]) or tended to be ([Formula: see text]) correlated to FCR. Great homogeneity in the enriched biological pathways was observed in these modules, suggesting co-expressed and co-regulated constitutive genes. They were mainly enriched in genes participating to immune and defense-related processes, and to a lesser extent, to translation, cell development or learning. They were also generally associated with growth rate and percentage of lean mass. In the whole network, only one module composed of genes participating to the response to substances, was significantly associated with daily feed intake and body adiposity. The plasma profiles in circulating metabolites and in fatty acids were summarized by weighted linear combinations using a dimensionality reduction method. Close association was thus found between a module composed of co-expressed genes participating to T cell receptor signaling and cell development process in the whole blood and related to FCR, and the circulating concentrations of polyunsaturated fatty acids in plasma., Conclusion: These systemic approaches have highlighted networks of entities driving key biological processes involved in the phenotypic difference in feed efficiency between animals. Connecting transcriptomics and metabolic levels together had some additional benefits., (© 2023. The Author(s).)

Background: Immune traits are considered to serve as potential biomarkers for pig's health. Medium to high heritabilities have been observed for some of the immune traits suggesting genetic variability of these phenotypes. Consideration of previously established genetic correlations between immune traits can be used to identify pleiotropic genetic markers. Therefore, genome-wide association study (GWAS) approaches are required to explore the joint genetic foundation for health biomarkers. Usually, GWAS explores phenotypes in a univariate (uv), trait-by-trait manner. Besides two uv GWAS methods, four multivariate (mv) GWAS approaches were applied on combinations out of 22 immune traits for Landrace (LR) and Large White (LW) pig lines., Results: In total 433 (LR: 351, LW: 82) associations were identified with the uv approach implemented in PLINK and a Bayesian linear regression uv approach (BIMBAM) software. Single Nucleotide Polymorphisms (SNPs) that were identified with both uv approaches (n = 32) were mostly associated with immune traits such as haptoglobin, red blood cell characteristics and cytokines, and were located in protein-coding genes. Mv GWAS approaches detected 647 associations for different mv immune trait combinations which were summarized to 133 Quantitative Trait Loci (QTL). SNPs for different trait combinations (n = 66) were detected with more than one mv method. Most of these SNPs are associated with red blood cell related immune trait combinations. Functional annotation of these QTL revealed 453 immune-relevant protein-coding genes. With uv methods shared markers were not observed between the breeds, whereas mv approaches were able to detect two conjoint SNPs for LR and LW. Due to unmapped positions for these markers, their functional annotation was not clarified., Conclusions: This study evaluated the joint genetic background of immune traits in LR and LW piglets through the application of various uv and mv GWAS approaches. In comparison to uv methods, mv methodologies identified more significant associations, which might reflect the pleiotropic background of the immune system more accurately. In genetic research of complex traits, the SNP effects are generally small. Furthermore, one genetic variant can affect several correlated immune traits at the same time, termed pleiotropy. As mv GWAS methods consider strong dependencies among traits, the power to detect SNPs can be boosted. Both methods revealed immune-relevant potential candidate genes. Our results indicate that one single test is not able to detect all the different types of genetic effects in the most powerful manner and therefore, the methods should be applied complementary., (© 2023. BioMed Central Ltd., part of Springer Nature.)

Background: Numerous circular RNAs (circRNAs) have been recently identified in porcine tissues and cell types. Nevertheless, their significance in porcine spleen development is yet unelucidated. Herein, we reported an extensive overlook of circRNA expression profile during spleen development in Meishan pigs., Results: Overall, 39,641 circRNAs were identified from 6,914 host genes. Among them, many circRNAs are up- or down-regulated at different time points of pig spleen development. Using WGCNA analysis, we revealed two essential modules for protein-coding genes and circRNAs. Subsequent correlation analysis revealed 67 circRNAs/co-expressed genes that participated in immnue-associated networks. Furthermore, a competing endogenous RNA (ceRNA) network analysis of circRNAs revealed that 12 circRNAs modulated CD226, MBD2, SAMD3, SIT1, SRP14, SYTL3 gene expressions via acting as miRNA sponges. Moreover, the circRNA&#95;21767/miR-202-3p axis regulated SIT1 expression in a ceRNA manner, which is critical for the immune-based regulation of spleen development in Meishan pigs., Conclusions: Overall, our results demonstrated that the circRNAs were differentially expressed during different stages of porcine spleen development, meanwhile the circRNAs interacted with immune-related genes in a ceRNA-based fashion. Moreover, we presented biomedical researchers with RNAseqTools, a user-friendly and powerful software for the visualization of transcriptome profile data., (© 2023. BioMed Central Ltd., part of Springer Nature.)

Background: This study investigated if the allele effect of a given single nucleotide polymorphism (SNP) for crossbred performance in pigs estimated in a genomic prediction model differs depending on its breed-of-origin, and how these are related to estimated effects for purebred performance. Results: SNP-allele substitution effects were estimated for a commonly used SNP panel using a genomic best linear unbiased prediction model with breed-specific partial relationship matrices. Estimated breeding values for purebred and crossbred performance were converted to SNP-allele effects by breed-of-origin. Differences between purebred and crossbred, and between breeds-of-origin were evaluated by comparing percentage of variance explained by genomic regions for back fat thickness (BF), average daily gain (ADG), and residual feed intake (RFI). From ten regions explaining most additive genetic variance for crossbred performance, 1 to 5 regions also appeared in the top ten for purebred performance. The proportion of genetic variance explained by a genomic region and the estimated effect of a haplotype in such a region were different depending upon the breed-of-origin. To illustrate underlying mechanisms, we evaluated the estimated effects across breeds-of-origin for haplotypes associated to the melanocortin 4 receptor (MC4R) gene, and for the MC4Rsnp itself which is a missense mutation with a known effect on BF and ADG. Although estimated allele substitution effects of the MC4Rsnp mutation were very similar across breeds, explained genetic variance of haplotypes associated to the MC4R gene using a SNP panel that does not include the mutation, was considerably lower in one of the breeds where the allele frequency of the mutation was the lowest. Conclusions: Similar regions explaining similar additive genetic variance were observed across purebred and crossbred performance. Moreover, there was some overlap across breeds-of-origin between regions that explained relatively large proportions of genetic variance for crossbred performance; albeit that the actual proportion of variance deviated across breeds-of-origin. Results based on a missense mutation in MC4R confirmed that even if a causal locus has similar effects across breeds-of-origin, estimated effects and explained variance in its region using a commonly used SNP panel can strongly depend on the allele frequency of the underlying causal mutation. [ABSTRACT FROM AUTHOR]

Background: Efforts to improve sustainability in livestock production systems have focused on two objectives: investigating the genetic control of immune function as it pertains to robustness and disease resistance, and finding predictive markers for use in breeding programs. In this context, the peripheral blood transcriptome represents an important source of biological information about an individual's health and immunological status, and has been proposed for use as an intermediate phenotype to measure immune capacity. The objective of this work was to study the genetic architecture of variation in gene expression in the blood of healthy young pigs using two approaches: an expression genome-wide association study (eGWAS) and allele-specific expression (ASE) analysis. Results: The blood transcriptomes of 60-day-old Large White pigs were analyzed by expression microarrays for eGWAS (242 animals) and by RNA-Seq for ASE analysis (38 animals). Using eGWAS, the expression levels of 1901 genes were found to be associated with expression quantitative trait loci (eQTLs). We recovered 2839 local and 1752 distant associations (Single Nucleotide Polymorphism or SNP located less or more than 1 Mb from expression probe, respectively). ASE analyses confirmed the extensive cis-regulation of gene transcription in blood, and revealed allelic imbalance in 2286 SNPs, which affected 763 genes. eQTLs and ASE-genes were widely distributed on all chromosomes. By analyzing mutually overlapping eGWAS results, we were able to describe putative regulatory networks, which were further refined using ASE data. At the functional level, genes with genetically controlled expression that were detected by eGWAS and/or ASE analyses were significantly enriched in biological processes related to RNA processing and immune function. Indeed, numerous distant and local regulatory relationships were detected within the major histocompatibility complex region on chromosome 7, revealing ASE for most class I and II genes. Conclusions: This study represents, to the best of our knowledge, the first genome-wide map of the genetic control of gene expression in porcine peripheral blood. These results represent an interesting resource for the identification of genetic markers and blood biomarkers associated with variations in immunity traits in pigs, as well as any other complex traits for which blood is an appropriate surrogate tissue. [ABSTRACT FROM AUTHOR]

Background: The intramuscular fat content (IMF) refers to the amount of fat within muscles, including the sum of phospholipids mainly found in cell membranes, triglycerides and cholesterol, and is determined both by hyperplasia and hypertrophy of adipocyte during the development of pigs. The IMF content is an important economic trait that is genetically controlled by multiple genes. The Laiwu pig is an indigenous fatty pig breed distributed in North China, characterized by excessively higher level of IMF content (9%~12%), therefore, is suitable for the identification of genes controlling IMF variations. To identify genes underlying IMF deposition, we performed genome-wide transcriptome and methylome analyses on longissimus dorsi (LD) muscle in Laiwu pigs across four developmental stages. Results: A total of 22,524 expressed genes were detected and 1158 differentially expressed genes (DEGs) were hierarchically clustered in the LD muscle over four developmental stages from 60 d to 400 d. These genes were significantly clustered into four temporal expression profiles and genes participating in fat cell differentiation and lipid biosynthesis processes were identified. From 120 d to 240 d, the period with the maximum IMF deposition rate, the lipid biosynthesis related genes (FOSL1, FAM213B and G0S2), transcription factors (TFs) (EGR1, KLF5, SREBF2, TP53 and TWIST1) and enriched pathways (steroid biosynthesis and fatty acid biosynthesis) were revealed; and fat biosynthesis relevant genes showing differences in DNA methylation in gene body or intergenic region were detected, such as FASN, PVALB, ID2, SH3PXD2B and EGR1. Conclusions: This study provides a comprehensive landscape of transcriptome of the LD muscle in Laiwu pigs ranging from 60 to 400 days old and methylome of the LD muscle in 120 d and 240 d Laiwu pigs. A set of candidate genes and TFs involved in fat biosynthesis process were identified, which were probably responsible for IMF deposition. The results from this study would provide a reference for the identification of genes controlling IMF variation and for exploring molecular mechanisms underlying IMF deposition in pigs. [ABSTRACT FROM AUTHOR]

Background: The Ashanti Dwarf Pig (ADP) of Ghana is an endangered pig breed with hardy and disease resistant traits. Characterisation of animal genetic resources provides relevant data for their conservation and sustainable use for food security and economic development. We investigated the origin and phylogenetic status of the local ADP of Ghana and their crosses with modern commercial breeds based on mtDNA, MC1R, Y-chromosome sequence polymorphisms, and genome-wide SNP genotyping. Results: The study involved 164 local pigs sampled from the three agro-ecological zones of Ghana. Analyses of the mitochondrial D-loop region and Y-chromosome sequences revealed both European and Asian genetic signatures, with differences between the geographical zones. Black coat colour is the most predominant within the breed, with black MC1R alleles of both Asian and European origin. European alleles for spotting are present at a low frequency in the sample set, and may account for the occurrence of spotted piglets in some APD litters. PCA analysis of SNP data revealed a strong location and breed effect on clustering of local Ghanaian pigs. On a global level, Ghanaian local pigs cluster closely with European pigs of commercial origin, but we identified intervals via FST analyses that may elucidate loci for ADP specific traits. Conclusions: The presence of both European and Asian contributions, with differences between geographical zones probably reflects trading and colonial influences. Understanding the effects of admixture on important adaptive and economic traits of the ADP and other local breeds in Africa is critical for developing sustainable conservation programmes to prevent the decline of these genetic resources. [ABSTRACT FROM AUTHOR]

Background: One of the most critical periods for the loss of pig embryos is the 12th day of gestation when implantation begins. Recent studies have shown that non-coding RNAs (ncRNAs) play important regulatory roles during pregnancy. Circular RNAs (circRNAs) are a kind of ubiquitously expressed ncRNAs that can directly regulate the binding proteins or regulate the expression of target genes by adsorbing micro RNAs (miRNA)., Results: We used the Illumina Novaseq6,000 technology to analyze the circRNA expression profile in the endometrium of three Erhualian (EH12) and three Yorkshire (YK12) pigs on day 12 of gestation. Overall, a total of 22,108 circRNAs were identified. Of these, 4051 circRNAs were specific to EH12 and 5889 circRNAs were specific to YK12, indicating a high level of breed specificity. Further analysis showed that there were 641 significant differentially expressed circRNAs (SDEcircRNAs) in EH12 compared with YK12 (FDR < 0.05). Functional enrichment of differential circRNA host genes revealed many pathways and genes associated with reproduction and regulation of embryo development. Network analysis of circRNA-miRNA interactions further supported the idea that circRNAs act as sponges for miRNAs to regulate gene expression. The prediction of differential circRNA binding proteins further explored the potential regulatory pathways of circRNAs. Analysis of SDEcircRNAs suggested a possible reason for the difference in embryo survival between the two breeds at the peri-implantation stage., Conclusions: Together, these data suggest that circRNAs are abundantly expressed in the endometrium during the peri-implantation period in pigs and are important regulators of related genes. The results of this study will help to further understand the differences in molecular pathways between the two breeds during the critical implantation period of pregnancy, and will help to provide insight into the molecular mechanisms that contribute to the establishment of pregnancy and embryo loss in pigs., (© 2023. The Author(s).)

Background: Feed restriction occurs frequently during pig growth, either due to economic reasons or stressful environmental conditions. Local breeds are suggested to have better tolerance to periods of feed restriction. However, the mechanisms underlying the response to feed restriction in different breeds is largely unknown. The aims of the present study were (1) to compare the blood transcriptome profile in response to feed restriction and refeeding of two contrasted breeds, Large White (LW), which has been selected for high performance, and Creole (CR), which is adapted to tropical conditions, and (2) to investigate the effect of a moderate feed restriction and refeeding on whole blood transcriptome. Analysis of blood transcriptome allows to study the response to feed restriction and refeeding in a dynamic way. RNAseq was performed on blood samples of growing LW and CR pigs at two time points: after 3 weeks of feed restriction and after 3 weeks of refeeding. The data was compared with samples from control animals offered the same diet on an ad libitum basis throughout the whole experiment., Results: In terms of performance (body weight and feed efficiency), CR pigs were less impacted by feed restriction than LW. The transcriptional response to feed restriction and refeeding between CR and LW was contrasted both in terms of number of DEGs and enriched pathways. CR demonstrated a stronger transcriptional response to feed restriction whereas LW had a stronger response to refeeding. Differences in the transcriptional response to feed restriction between CR and LW were related to cell stress response (Aldosterone Signalling, Protein ubiquitination, Unfolded Protein Signalling) whereas after refeeding, differences were linked to thermogenesis, metabolic pathways and cell proliferation (p38 MAPK, ERK/MAPK pathway). In both breeds, transcriptional changes related to the immune response were found after restriction and refeeding., Conclusions: Altogether, the present study indicates that blood transcriptomics can be a useful tool to study differential genetic response to feed restriction in a dynamic way. The results indicate a differential response of blood gene expression to feed restriction and refeeding between breeds, affecting biological pathways that are in accordance with performance and thermoregulatory results., (© 2023. The Author(s).)

Background: Alternative splicing is an important step in gene expression, generating multiple isoforms for the same genes and greatly expanding the diversity of proteomes. Genetic variation in alternative splicing contributes to phenotypic diversity in natural populations. However, the genetic basis of variation in alternative splicing in livestock including pigs remains poorly understood., Results: In this study, using a Duroc x Pietrain F2 pig population, we performed genome-wide analysis of alternative splicing estimated from stranded RNA-Seq data in skeletal muscle. We characterized the genetic architecture of alternative splicing and compared its basic features with those of overall gene expression. We detected a large number of novel alternative splicing events that were not previously annotated. We found heritability of quantitative alternative splicing scores (percent spliced in or PSI) to be lower than that of overall gene expression. In addition, heritabilities showed little correlation between alternative splicing and overall gene expression. We mapped expression QTLs (eQTLs) and splice QTLs (sQTLs) and found them to be largely non-overlapping. Finally, we integrated sQTL mapping with phenotype QTL (pQTL mapping to identify potential mediator of pQTL effect by alternative splicing., Conclusions: Our results suggest that regulatory variation exists at multiple levels and that their genetic controls are distinct, offering opportunities for genetic improvement., (© 2023. The Author(s).)

Background: N 6 -methyladenosine (m 6 A) refers to the methylation modification of N 6 position of RNA adenine, a dynamic reversible RNA epigenetic modification that plays an important regulatory role in a variety of life processes. In this study, we used MeRIP-Seq and RNA-Seq of the longissimus dorsi (LD) muscle of adult (QA) and newborn (QN) Queshan Black pigs to screen key genes with m 6 A modification involved in muscle growth by bioinformatics analysis., Results: A total of 23,445 and 25,465 m 6 A peaks were found in the whole genomes of QA and QN, respectively. Among them, 613 methylation peaks were significantly different (DMPs) and 579 genes were defined as differentially methylated genes (DMGs). Compared with the QN group, there were 1,874 significantly differentially expressed genes (DEGs) in QA group, including 620 up-regulated and 1,254 down-regulated genes. In order to investigate the relationship between m 6 A and mRNA expression in the muscle of Queshan Black pigs at different periods, a combined analysis of MeRIP-Seq and RNA-Seq showed that 88 genes were significantly different at both levels. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes results showed that DEGs and DMGs were mainly involved in skeletal muscle tissue development, FoxO signaling pathway, MAPK signaling pathway, insulin signaling pathway, PI3K-Akt signaling pathway, and Wnt signaling pathway. Four DEGs (IGF1R, CCND2, MYOD1 and FOS) and four DMGs (CCND2, PHKB, BIN1 and FUT2), which are closely related to skeletal muscle development, were selected as candidate genes for verification, and the results were consistent with the sequencing results, which indicated the reliability of the sequencing results., Conclusions: These results lay the foundation for understanding the specific regulatory mechanisms of growth in Queshan Black pigs, and provide theoretical references for further research on the role of m 6 A in muscle development and breed optimization selection., (© 2023. The Author(s).)

Background: De novo mutations arising in the germline are a source of genetic variation and their discovery broadens our understanding of genetic disorders and evolutionary patterns. Although the number of de novo single nucleotide variants (dnSNVs) has been studied in a number of species, relatively little is known about the occurrence of de novo structural variants (dnSVs). In this study, we investigated 37 deeply sequenced pig trios from two commercial lines to identify dnSVs present in the offspring. The identified dnSVs were characterised by identifying their parent of origin, their functional annotations and characterizing sequence homology at the breakpoints., Results: We identified four swine germline dnSVs, all located in intronic regions of protein-coding genes. Our conservative, first estimate of the swine germline dnSV rate is 0.108 (95% CI 0.038-0.255) per generation (one dnSV per nine offspring), detected using short-read sequencing. Two detected dnSVs are clusters of mutations. Mutation cluster 1 contains a de novo duplication, a dnSNV and a de novo deletion. Mutation cluster 2 contains a de novo deletion and three de novo duplications, of which one is inverted. Mutation cluster 2 is 25 kb in size, whereas mutation cluster 1 (197 bp) and the other two individual dnSVs (64 and 573 bp) are smaller. Only mutation cluster 2 could be phased and is located on the paternal haplotype. Mutation cluster 2 originates from both micro-homology as well as non-homology mutation mechanisms, where mutation cluster 1 and the other two dnSVs are caused by mutation mechanisms lacking sequence homology. The 64 bp deletion and mutation cluster 1 were validated through PCR. Lastly, the 64 bp deletion and the 573 bp duplication were validated in sequenced offspring of probands with three generations of sequence data., Conclusions: Our estimate of 0.108 dnSVs per generation in the swine germline is conservative, due to our small sample size and restricted possibilities of dnSV detection from short-read sequencing. The current study highlights the complexity of dnSVs and shows the potential of breeding programs for pigs and livestock species in general, to provide a suitable population structure for identification and characterisation of dnSVs., (© 2023. The Author(s).)

Background: As a diverse and abundant class of endogenous RNAs, circular RNAs (circRNAs) participate in various biological processes including cell proliferation and apoptosis. Nevertheless, few researchers have investigated the role of circRNAs in muscle development in cultivated pigs., Results: In this study, we used RNA-seq to construct circRNA expression profiles in skeletal muscle of Jinfen White pigs at the age of 1, 90, and 180 days. Among the 16,990 identified circRNAs, 584 circRNAs were differentially expressed. Moreover, the enrichment analysis of DE circRNA host genes showed that they were mainly involved in muscle contraction, muscle organ development and muscle system processes, as well as AMPK and cAMP-related signal pathways. We also constructed a circRNA-miRNA-mRNA co-expression network to find key circRNAs which many involved in the regulation of porcine skeletal muscle development through the competitive endogenous RNA (ceRNA) mechanism. It is noteworthy that circ&#95;0018595/miR-1343/PGM1 axis may play a regulatory role in the development of porcine skeletal muscle., Conclusions: This study identified the circRNAs and present the circRNA expression profile in the development of pigs, revealed that DE circRNA host genes participate in different cell fates and enriched the porcine ceRNA network. Thus, this work will become a valuable resource for further in-depth study of the regulatory mechanism of circRNA in the development of porcine skeletal muscle., (© 2023. The Author(s).)

Background: Genetic correlations between complex traits suggest that pleiotropic variants contribute to trait variation. Genome-wide association studies (GWAS) aim to uncover the genetic underpinnings of traits. Multivariate association testing and the meta-analysis of summary statistics from single-trait GWAS enable detecting variants associated with multiple phenotypes. In this study, we used array-derived genotypes and phenotypes for 24 reproduction, production, and conformation traits to explore differences between the two methods and used imputed sequence variant genotypes to fine-map six quantitative trait loci (QTL)., Results: We considered genotypes at 44,733 SNPs for 5,753 pigs from the Swiss Large White breed that had deregressed breeding values for 24 traits. Single-trait association analyses revealed eleven QTL that affected 15 traits. Multi-trait association testing and the meta-analysis of the single-trait GWAS revealed between 3 and 6 QTL, respectively, in three groups of traits. The multi-trait methods revealed three loci that were not detected in the single-trait GWAS. Four QTL that were identified in the single-trait GWAS, remained undetected in the multi-trait analyses. To pinpoint candidate causal variants for the QTL, we imputed the array-derived genotypes to the sequence level using a sequenced reference panel consisting of 421 pigs. This approach provided genotypes at 16 million imputed sequence variants with a mean accuracy of imputation of 0.94. The fine-mapping of six QTL with imputed sequence variant genotypes revealed four previously proposed causal mutations among the top variants., Conclusions: Our findings in a medium-size cohort of pigs suggest that multivariate association testing and the meta-analysis of summary statistics from single-trait GWAS provide very similar results. Although multi-trait association methods provide a useful overview of pleiotropic loci segregating in mapping populations, the investigation of single-trait association studies is still advised, as multi-trait methods may miss QTL that are uncovered in single-trait GWAS., (© 2023. The Author(s).)

Background: In many traits, not only individual trait levels are under genetic control, but also the variation around that level. In other words, genotypes do not only differ in mean, but also in (residual) variation around the genotypic mean. New statistical methods facilitate gaining knowledge on the genetic architecture of complex traits such as phenotypic variability. Here we study litter size (total number born) and its variation in a Large White pig population using a Double Hierarchical Generalized Linear model, and perform a genome-wide association study using a Bayesian method. Results: In total, 10 significant single nucleotide polymorphisms (SNPs) were detected for total number born (TNB) and 9 SNPs for variability of TNB (varTNB). Those SNPs explained 0.83 % of genetic variance in TNB and 1.44 % in varTNB. The most significant SNP for TNB was detected on Sus scrofa chromosome (SSC) 11. A possible candidate gene for TNB is ENOX1, which is involved in cell growth and survival. On SSC7, two possible candidate genes for varTNB are located. The first gene is coding a swine heat shock protein 90 (HSPCB = Hsp90), which is a well-studied gene stabilizing morphological traits in Drosophila and Arabidopsis. The second gene is VEGFA, which is activated in angiogenesis and vasculogenesis in the fetus. Furthermore, the genetic correlation between additive genetic effects on TNB and on its variation was 0.49. This indicates that the current selection to increase TNB will also increase the varTNB. Conclusions: To the best of our knowledge, this is the first study reporting SNPs associated with variation of a trait in pigs. Detected genomic regions associated with varTNB can be used in genomic selection to decrease varTNB, which is highly desirable to avoid very small or very large litters in pigs. However, the percentage of variance explained by those regions was small. The SNPs detected in this study can be used as indication for regions in the Sus scrofa genome involved in maintaining low variability of litter size, but further studies are needed to identify the causative loci. [ABSTRACT FROM AUTHOR]

Background: Pigs (Sus scrofa) provide relevant biomedical models to dissect complex diseases due to their anatomical, genetic, and physiological similarities with humans. Aberrant DNA methylation has been linked to many of these diseases and is associated with gene expression; however, the functional similarities and differences between porcine and human DNA methylation patterns are largely unknown. Methods: DNA and RNA was isolated from eight tissue samples (fat, heart, kidney, liver, lung, lymph node, muscle, and spleen) from the adult female Duroc utilized for the pig genome sequencing project. Reduced representation bisulfite sequencing (RRBS) and RNA-seq were performed on an Illumina HiSeq2000. RRBS reads were aligned using BSseeker2, and only sites with a minimum depth of 10 reads were used for methylation analysis. RNA-seq reads were aligned using Tophat, and expression analysis was performed using Cufflinks. In addition, SNP calling was performed using GATK for targeted control and whole genome sequencing reads for CpG site validation and allelic expression analysis, respectively. Results: Analysis on the influence of DNA variation in methylation calling revealed a reduced effectiveness of WGS datasets in covering CpG rich regions, as well as the usefulness of a targeted control library for SNP detection. Analysis of over 500,000 CpG sites demonstrated genome wide methylation patterns similar to those observed in humans, including reduced methylation within CpG islands and at transcription start sites (TSS), X chromosome inactivation, and anticorrelation of TSS CpG methylation with gene expression. In addition, a positive correlation between TSS CpG density and expression, and a negative correlation between TSS TpG density and expression were demonstrated. Low but non-random non-CpG methylation (<1%) was also detected in all non-neuronal somatic tissues, with differences in tissue clustering observed based on CpG and non-CpG methylation patterns. Finally, allele specific expression analysis revealed enrichment of genes involved in metabolic and regulatory processes. Discussion: These results provide transcriptional and DNA methylation datasets for the biomedical community that are directly relatable to current genomic resources. In addition, the correlation between TSS CpG density and expression suggests increased mutation rates at CpG sites play a significant role in adaptive evolution by reducing CpG density at TSS over time, resulting in higher methylation levels in these regions and more permanent changes to lower gene expression. This is proposed to occur predominantly through deamination of 5-methylcytosine to thymidine, resulting in the replacement of CpG with TpG sites in these regions, as indicated by the increased TSS TpG density observed in non-expressed genes, resulting in a negative correlation between expression and TSS TpG density. Conclusions: This study provides baseline methylation and gene transcription profiles for a healthy adult pig, reports similar patterns to those observed in humans, and supports future porcine studies related to human disease and development. Additionally, the observed reduced CpG and increased TpG density at TSS of lowly expressed genes suggests DNA methylation plays a significant role in adaptive evolution through more permanent changes to lower gene expression. [ABSTRACT FROM AUTHOR]

Background: The growth and development of skeletal muscle directly impacts the quantity and quality of pork production. Chinese indigenous pig breeds and exotic species vary greatly in terms of muscle production and performance traits. We present transcriptome profiles of 110 skeletal muscle samples from Tongcheng (TC) and Yorkshire (YK) pigs at 11 developmental periods (30, 40, 55, 63, 70, 90, and 105 days of gestation, and 0, 1, 3, and 5 weeks of age) using digital gene expression on Solexa/Illumina's Genome Analyzer platform to investigate the differences in prenatal and postnatal skeletal muscle between the two breeds. Results: Muscle morphological changes indicate the importance of primary fiber formation from 30 to 40 dpc (days post coitus), and secondary fiber formation from 55 to 70 dpc. We screened 4,331 differentially expressed genes in TC and 2,259 in YK (log2 ratio >1 and probability >0.7). Cluster analysis showed different gene expression patterns between TC and YK pigs. The transcripts were annotated in terms of Gene Ontology related to muscle development. We found that the genes CXCL10, EIF2B5, PSMA6, FBXO32, and LOC100622249 played vital roles in the muscle regulatory networks in the TC breed, whereas the genes SGCD, ENG, THBD, AQP4, and BTG2 played dominant roles in the YK breed. These genes showed breed-specific and development-dependent differential expression patterns. Furthermore, 984 genes were identified in myogenesis. A heat map showed that significantly enriched pathways (FDR <0.05) had stage-specific functional regulatory mechanisms. Finally, the differentially expressed genes from our sequencing results were confirmed by real-time quantitative polymerase chain reaction. Conclusions: This study detected many functional genes and showed differences in the molecular mechanisms of skeletal muscle development between TC and YK pigs. TC pigs showed slower muscle growth and more complicated genetic regulation than YK pigs. Many differentially expressed genes showed breed-specific expression patterns. Our data provide a better understanding of skeletal muscle developmental differences and valuable information for improving pork quality. [ABSTRACT FROM AUTHOR]

Background: Unraveling the genetic mechanisms associated with reduced gene flow between genetically differentiated populations is key to understand speciation. Different types of structural variations (SVs) have been found as a source of genetic diversity in a wide range of species. Previous studies provided detailed knowledge on the potential evolutionary role of SVs, especially copy number variations (CNVs), between well diverged species of e.g. primates. However, our understanding of their significance during ongoing speciation processes is limited due to the lack of CNV data from closely related species. The genus Sus (pig and its close relatives) which started to diverge ~4 Mya presents an excellent model for studying the role of CNVs during ongoing speciation. Results: In this study, we identified 1408 CNV regions (CNVRs) across the genus Sus. These CNVRs encompass 624 genes and were found to evolve ~2.5 times faster than single nucleotide polymorphisms (SNPs). The majority of these copy number variable genes are olfactory receptors (ORs) known to play a prominent role in food foraging and mate recognition in Sus. Phylogenetic analyses, including novel Bayesian analysis, based on CNVRs that overlap ORs retain the well-accepted topology of the genus Sus whereas CNVRs overlapping genes other than ORs show evidence for random drift and/or admixture. Conclusion: We hypothesize that inter-specific variation in copy number of ORs provided the means for rapid adaptation to different environments during the diversification of the genus Sus in the Pliocene. Furthermore, these regions might have acted as barriers preventing massive gene flow between these species during the multiple hybridization events that took place later in the Pleistocene suggesting a possible prominent role of ORs in the ongoing Sus speciation. [ABSTRACT FROM AUTHOR]

Background: BMPR1B (Bone morphogenetic protein receptor type-1B) is a receptor in the bone morphogenetic protein (BMP) family and has been identified as a candidate gene for reproductive traits in pigs. Our previous study in Taihu pigs found a specific estrogen response element (ERE) in the first intron of the BMPR1B gene that is associated with the number born alive trait. However, little is known about the mechanism by which the ERE regulates the expression of BMPR1B in the endometrium., Results: Here, a 15-bp InDel (insertion/deletion) (AGCCAGAAAGGAGGA) was identified as a unique variation in Taihu pigs, and was shown to be responsible for the binding of the type I receptor of estrogen (ESR1) to the ERE using dual-luciferase assays. Four BMPR1B transcripts (T1, T2, T3, and T4) were identified by 5' RACE in endometrial tissue. Expression of T3 and T4 in the endometrium of Meishan pigs was significantly higher than in Duroc pigs during pregnancy. Luciferase assays showed that three distinct BMPR1B promoters may drive expression of T1, T3, and T4. Interestingly, ERE-mediated enhancement of T4 promoter activity significantly increased expression of Transcript T4 in the endometrium of Taihu pigs (P < 0.05). In contrast, the ERE inhibited activity of the T3 promoter and decreased expression of the T3 transcript in the Duroc background (P < 0.05). In summary, we identified a 15-bp InDel in the Taihu ERE that can be used as a molecular marker for the number born alive trait, characterized the 5' untranslated regions (UTRs) of BMPR1B transcripts in the endometrium, and determined how the transcripts are processed by alternative splicing events., Conclusions: Our results provide a foundation for understanding the transcriptional regulation of BMPR1B and its contributions to the unique breeding prolificacy characteristics of Taihu pigs., (© 2022. The Author(s).)

Background: Average backfat thickness (BFT) is a critical complex trait in pig and an important indicator for fat deposition and lean rate. Usually, genome-wide association study (GWAS) was used to discover quantitative trait loci (QTLs) of BFT in a single population. However, the power of GWAS is limited by sample size in a single population. Alternatively, meta-analysis of GWAS (metaGWAS) is an attractive method to increase the statistical power by integrating data from multiple breeds and populations. The aim of this study is to identify shared genetic characterization of BFT across breeds in pigs via metaGWAS.  RESULTS: In this study, we performed metaGWAS on BFT using 15,353 pigs (5,143 Duroc, 7,275 Yorkshire, and 2,935 Landrace) from 19 populations. We detected 40 genome-wide significant SNPs (Bonferroni corrected P < 0.05) and defined five breed-shared QTLs in across-breed metaGWAS. Markers within the five QTL regions explained 7 ~ 9% additive genetic variance and showed strong heritability enrichment. Furthermore, by integrating information from multiple bioinformatics databases, we annotated 46 candidate genes located in the five QTLs. Among them, three important (MC4R, PPARD, and SLC27A1) and seven suggestive candidate genes (PHLPP1, NUDT3, ILRUN, RELCH, KCNQ5, ITPR3, and U3) were identified., Conclusion: QTLs and candidate genes underlying BFT across breeds were identified via metaGWAS from multiple populations. Our findings contribute to the understanding of the genetic architecture of BFT and the regulating mechanism underlying fat deposition in pigs., (© 2022. The Author(s).)

Background: The transcriptional changes around zygotic genome activation (ZGA) in preimplantation embryos are critical for studying mechanisms of embryonic developmental arrest and searching for key transcription factors. However, studies on the transcription profile of porcine ZGA are limited., Results: In this study, we performed RNA sequencing in porcine in vivo developed (IVV) and somatic cell nuclear transfer (SCNT) embryo at different stages and compared the transcriptional activity of porcine embryos with mouse, bovine and human embryos. The results showed that the transcriptome map of the early porcine embryos was significantly changed at the 4-cell stage, and 5821 differentially expressed genes (DEGs) in SCNT embryos failed to be reprogrammed or activated during ZGA, which mainly enrichment to metabolic pathways. c-MYC was identified as the highest expressed transcription factor during ZGA. By treating with 10,058-F4, an inhibitor of c-MYC, the cleavage rate (38.33 ± 3.4%) and blastocyst rate (23.33 ± 4.3%) of porcine embryos were significantly lower than those of the control group (50.82 ± 2.7% and 34.43 ± 1.9%). Cross-species analysis of transcriptome during ZGA showed that pigs and bovines had the highest similarity coefficient in biological processes. KEGG pathway analysis indicated that there were 10 co-shared pathways in the four species., Conclusions: Our results reveal that embryos with impaired developmental competence may be arrested at an early stage of development. c-MYC helps promote ZGA and preimplantation embryonic development in pigs. Pigs and bovines have the highest coefficient of similarity in biological processes during ZGA. This study provides an important reference for further studying the reprogramming regulatory mechanism of porcine embryos during ZGA., (© 2022. The Author(s).)

Background The two main genetic types in Iberian pig production show important phenotypic differences in growth, fattening and tissue composition since early developmental stages. The objective of this work was the evaluation of muscle transcriptome profile in piglets of both genetic types, in order to identify genes, pathways and regulatory factors responsible for their phenotypic differences. Contemporary families coming from pure Iberian pigs (IB) or from crossing with Duroc boars (DUxIB) were generated. Piglets (14 from each genetic type) were slaughtered at weaning (28 days) and longissimus dorsi was sampled for composition and gene expression studies. RNA was obtained and hybridized to Affymetrix Porcine Genechip expression arrays. Results Loin muscle chemical composition showed significant differences between genetic types in intramuscular fat content (6.1% vs. 4.3% in IB and DUxIB animals, respectively, P = 0.009) and in saturated (P = 0.019) and monounsaturated fatty acid proportions (P = 0.044). The statistical analysis of gene expression data allowed the identification of 256 differentially expressed (DE) genes between genetic types (FDR < 0.10), 102 upregulated in IB and 154 upregulated in DUxIB. Transcript differences were validated for a subset of DE genes by qPCR. We observed alteration in biological functions related to extracellular matrix function and organization, cellular adhesion, muscle growth, lipid metabolism and proteolysis. Candidate genes with known effects on muscle growth were found among the DE genes upregulated in DUxIB. Genes related to lipid metabolism and proteolysis were found among those upregulated in IB. Regulatory factors (RF) potentially involved in the expression differences were identified by calculating the regulatory impact factors. Twenty-nine RF were found, some of them with known relationship with tissue development (MSTN, SIX4, IRX3), adipogenesis (CEBPD, PPARGC1B), or extracellular matrix processes (MAX, MXI1). Correlation among the expression of these RF and DE genes show relevant differences between genetic types. Conclusion These results provide valuable information about genetic mechanisms determining the phenotypic differences on growth and meat quality between the genetic types studied, mainly related to the development and function of the extracellular matrix and also to some metabolic processes as proteolysis and lipid metabolism. Transcription factors and regulatory mechanisms are proposed for these altered biological functions. [ABSTRACT FROM AUTHOR]

Background Obesity, excess fat tissue in the body, can underlie a variety of medical complaints including heart disease, stroke and cancer. The pig is an excellent model organism for the study of various human disorders, including obesity, as well as being the foremost agricultural species. In order to identify genetic variants associated with fatness, we used a selective genomic approach sampling DNA from animals at the extreme ends of the fat and lean spectrum using estimated breeding values derived from a total population size of over 70,000 animals. DNA from 3 breeds (Sire Line Large White, Duroc and a white Pietrain composite line (Titan)) was used to interrogate the Illumina Porcine SNP60 Genotyping Beadchip in order to identify significant associations in terms of single nucleotide polymorphisms (SNPs) and copy number variants (CNVs). Results By sampling animals at each end of the fat/lean EBV (estimate breeding value) spectrum the whole population could be assessed using less than 300 animals, without losing statistical power. Indeed, several significant SNPs (at the 5% genome wide significance level) were discovered, 4 of these linked to genes with ontologies that had previously been correlated with fatness (NTS, FABP6, SST and NR3C2). Quantitative analysis of the data identified putative CNV regions containing genes whose ontology suggested fatness related functions (MCHR1, PPARα, SLC5A1 and SLC5A4). Conclusions Selective genotyping of EBVs at either end of the phenotypic spectrum proved to be a cost effective means of identifying SNPs and CNVs associated with fatness and with estimated major effects in a large population of animals. [ABSTRACT FROM AUTHOR]

Background: Structural variants and tandem repeats are relevant sources of genomic variation that are not routinely analyzed in genome wide association studies mainly due to challenging identification and genotyping. Here, we profiled these variants via state-of-the-art strategies in the founder animals of four F 2 pig crosses using whole-genome sequence data (20x coverage). The variants were compared at a founder level with the commonly screened SNPs and small indels. At the F 2 level, we carried out an association study using imputed structural variants and tandem repeats with four growth and carcass traits followed by a comparison with a previously conducted SNPs and small indels based association study., Results: A total of 13,201 high confidence structural variants and 103,730 polymorphic tandem repeats (with a repeat length of 2-20 bp) were profiled in the founders. We observed a moderate to high (r from 0.48 to 0.57) level of co-localization between SNPs or small indels and structural variants or tandem repeats. In the association step 56.56% of the significant variants were not in high LD with significantly associated SNPs and small indels identified for the same traits in the earlier study and thus presumably not tagged in case of a standard association study. For the four growth and carcass traits investigated, many of the already proposed candidate genes in our previous studies were confirmed and additional ones were identified. Interestingly, a common pattern on how structural variants or tandem repeats regulate the phenotypic traits emerged. Many of the significant variants were embedded or nearby long non-coding RNAs drawing attention to their functional importance. Through which specific mechanisms the identified long non-coding RNAs and their associated structural variants or tandem repeats contribute to quantitative trait variation will need further investigation., Conclusions: The current study provides insights into the characteristics of structural variants and tandem repeats and their role in association studies. A systematic incorporation of these variants into genome wide association studies is advised. While not of immediate interest for genomic prediction purposes, this will be particularly beneficial for elucidating biological mechanisms driving the complex trait variation., (© 2022. The Author(s).)

Background: Carcass backfat thickness (BFT), carcass lean percentage (CLP) and carcass fat percentage (CFP) are important to the commercial pig industry. Nevertheless, the genetic architecture of BFT, CLP and CFP is still elusive. Here, we performed a genome-wide association study (GWAS) based on specific-locus amplified fragment sequencing (SLAF-seq) to analyze seven fatness-related traits, including five BFTs, CLP, and CFP on 223 four-way crossbred pigs., Results: A total of 227, 921 highly consistent single nucleotide polymorphisms (SNPs) evenly distributed throughout the genome were used to perform GWAS. Using the mixed linear model (MLM), a total of 20 SNP loci significantly related to these traits were identified on ten Sus scrofa chromosomes (SSC), of which 10 SNPs were located in previously reported quantitative trait loci (QTL) regions. On SSC7, two SNPs (SSC7:29,503,670 and rs1112937671) for average backfat thickness (ABFT) exceeded 1% and 10% Bonferroni genome-wide significance levels, respectively. These two SNP loci were located within an intron region of the COL21A1 gene, which was a protein-coding gene that played an important role in the porcine backfat deposition by affecting extracellular matrix (ECM) remodeling. In addition, based on the other three significant SNPs on SSC7, five candidate genes, ZNF184, ZNF391, HMGA1, GRM4 and NUDT3 were proposed to influence BFT. On SSC9, two SNPs for backfat thickness at 6-7 ribs (67RBFT) and one SNP for CLP were in the same locus region (19 kb interval). These three SNPs were located in the PGM2L1 gene, which encoded a protein that played an indispensable role in glycogen metabolism, glycolysis and gluconeogenesis as a key enzyme. Finally, one significant SNP on SSC14 for CLP was located within the PLBD2 gene, which participated in the lipid catabolic process., Conclusions: A total of two regions on SSC7 and SSC9 and eight potential candidate genes were found for fatness-related traits in pigs. The results of this GWAS based on SLAF-seq will greatly advance our understanding of the genetic architecture of BFT, CLP, and CFP traits. These identified SNP loci and candidate genes might serve as a biological basis for improving the important fatness-related traits of pigs., (© 2022. The Author(s).)

Background: Carcass traits are important in pig breeding programs for improving pork production. Understanding the genetic variants underlies complex phenotypes can help explain trait variation in pigs. In this study, we integrated a weighted single-step genome-wide association study (wssGWAS) and copy number variation (CNV) analyses to map genetic variations and genes associated with loin muscle area (LMA), loin muscle depth (LMD) and lean meat percentage (LMP) in Duroc pigs., Results: Firstly, we performed a genome-wide analysis for CNV detection using GeneSeek Porcine SNP50 Bead chip data of 3770 pigs. A total of 11,100 CNVs were detected, which were aggregated by overlapping 695 CNV regions (CNVRs). Next, we investigated CNVs of pigs from the same population by whole-genome resequencing. A genome-wide analysis of 21 pigs revealed 23,856 CNVRs that were further divided into three categories (851 gain, 22,279 loss, and 726 mixed), which covered 190.8 Mb (~ 8.42%) of the pig autosomal genome. Further, the identified CNVRs were used to determine an overall validation rate of 68.5% for the CNV detection accuracy of chip data. CNVR association analyses identified one CNVR associated with LMA, one with LMD and eight with LMP after applying stringent Bonferroni correction. The wssGWAS identified eight, six and five regions explaining more than 1% of the additive genetic variance for LMA, LMD and LMP, respectively. The CNVR analyses and wssGWAS identified five common regions, of which three regions were associated with LMA and two with LMP. Four genes (DOK7, ARAP1, ELMO2 and SLC13A3) were highlighted as promising candidates according to their function., Conclusions: We determined an overall validation rate for the CNV detection accuracy of low-density chip data and constructed a genomic CNV map for Duroc pigs using resequencing, thereby proving a value genetic variation resource for pig genome research. Furthermore, our study utilized a composite genetic strategy for complex traits in pigs, which will contribute to the study for elucidating the genetic architecture that may be influenced and regulated by multiple forms of variations., (© 2022. The Author(s).)

Background: Quantitative real time PCR (qPCR) is a powerful tool to evaluate mRNA expression level. However, reliable qPCR results require normalization with validated reference gene(s). In this study, we investigated stable reference genes in seven tissues according to four developmental stages in minipigs. Six candidate reference genes and one target gene (ACE2) were selected and qPCR was performed. BestKeeper, geNorm, NormFinder, and delta Ct method through the RefFinder web-based tool were used to evaluate the stability of candidate reference genes. To verify the selected stable genes, relative expression of ACE2 was calculated and compared with each other., Results: As a result, HPRT1 and 18S genes had lower SD value, while HMBS and GAPDH genes had higher SD value in all samples. Using statistical algorithms, HPRT1 was the most stable gene, followed by 18S, β-actin, B2M, GAPDH, and HMBS. In intestine, all candidate reference genes exhibited similar patterns of ACE2 gene expression over time, whereas in liver, lung, and kidney, gene expression pattern normalized with stable reference genes differed from those normalized with less stable genes. When normalized with the most stable genes, the expression levels of ACE2 in minipigs highly increased in intestine and kidney at PND28, which is consistent with the ACE2 expression pattern in humans., Conclusions: We suggest that HPRT1 and 18S are good choices for analyzing all these samples across the seven tissues and four developmental stages. However, this study can be a reference literature for gene expression experiments using minipig because reference gene should be validated and chosen according to experimental conditions., (© 2022. The Author(s).)

Background: Past selection events left footprints in the genome of domestic animals, which can be traced back by stretches of homozygous genotypes, designated as runs of homozygosity (ROHs). The analysis of common ROH regions within groups or populations displaying potential signatures of selection requires high-quality SNP data as well as carefully adjusted ROH-defining parameters. In this study, we used a simultaneous testing of rule- and model-based approaches to perform strategic ROH calling in genomic data from different pig populations to detect genomic regions under selection for specific phenotypes., Results: Our ROH analysis using a rule-based approach offered by PLINK, as well as a model-based approach run by RZooRoH demonstrated a high efficiency of both methods. It underlined the importance of providing a high-quality SNP set as input as well as adjusting parameters based on dataset and population for ROH calling. Particularly, ROHs ≤ 20 kb were called in a high frequency by both tools, but to some extent covered different gene sets in subsequent analysis of ROH regions common for investigated pig groups. Phenotype associated ROH analysis resulted in regions under potential selection characterizing heritage pig breeds, known to harbour a long-established breeding history. In particular, the selection focus on fitness-related traits was underlined by various ROHs harbouring disease resistance or tolerance-associated genes. Moreover, we identified potential selection signatures associated with ear morphology, which confirmed known candidate genes as well as uncovered a missense mutation in the ABCA6 gene potentially supporting ear cartilage formation., Conclusions: The results of this study highlight the strengths and unique features of rule- and model-based approaches as well as demonstrate their potential for ROH analysis in animal populations. We provide a workflow for ROH detection, evaluating the major steps from filtering for high-quality SNP sets to intersecting ROH regions. Formula-based estimations defining ROHs for rule-based method show its limits, particularly for efficient detection of smaller ROHs. Moreover, we emphasize the role of ROH detection for the identification of potential footprints of selection in pigs, displaying their breed-specific characteristics or favourable phenotypes., (© 2022. The Author(s).)

Background Immune traits (ITs) are potentially relevant criteria to characterize an individual's immune response. Our aim was to investigate whether the peripheral blood transcriptome can provide a significant and comprehensive view of IT variations in pig. Results Sixty-day-old Large White pigs classified as extreme for in vitro production of IL2, IL10, IFNγ and TNFα, phagocytosis activity, in vivo CD4-/CD8+ or TCRγδ + cell counts, and anti- Mycoplasma antibody levels were chosen to perform a blood transcriptome analysis with a porcine generic array enriched with immunity-related genes. Differentially expressed (DE) genes for in vitro production of IL2 and IL10, phagocytosis activity and CD4-/CD8+ cell counts were identified. Gene set enrichment analysis revealed a significant overrepresentation of immune response functions. To validate the microarray-based results, a subset of DE genes was confirmed by RT-qPCR. An independent set of 74 animals was used to validate the covariation between gene expression levels and ITs. Five potential gene biomarkers were found for prediction of IL2 (RALGDS), phagocytosis (ALOX12) or CD4- /CD8+ cell count (GNLY, KLRG1 and CX3CR1). On average, these biomarkers performed with a sensitivity of 79% and a specificity of 86%. Conclusions Our results confirmed that gene expression profiling in blood represents a relevant molecular phenotype to refine ITs in pig and to identify potential biomarkers that can provide new insights into immune response analysis. [ABSTRACT FROM AUTHOR]

Background: Runs of homozygosity (ROH) are continuous homozygous regions typically located in the DNA sequence of diploid organisms. Identifications of ROH that lead to reduced performance can provide valuable insight into the genetic architecture of complex traits. Here, we systematically investigated the population genetic structure of five Anhui indigenous pig breeds (AHIPs), and compared them to those of five Western commercial pig breeds (WECPs). Furthermore, we examined the occurrence and distribution of ROHs in the five AHIPs and estimated the inbreeding coefficients based on the ROHs (F ROH ) and homozygosity (F HOM ). Finally, we identified genomic regions with high frequencies of ROHs and annotated candidate genes contained therein., Results: The WECPs and AHIPs were clearly differentiated into two separate clades consistent with their geographical origins, as revealed by the population structure and principal component analysis. We identified 13,530 ROHs across all individuals, of which 4,555 and 8,975 ROHs were unique to AHIPs and WECPs, respectively. Most ROHs identified in our study were short (< 10 Mb) or medium (10-20 Mb) in length. WECPs had significantly higher numbers of short ROHs, and AHIPs generally had longer ROHs. F ROH values were significantly lower in AHIPs than in WECPs, indicating that breed improvement and conservation programmes were successful in AHIPs. On average, F ROH and F HOM values were highly correlated (0.952-0.991) in AHIPs and WECPs. A total of 27 regions had a high frequency of ROHs and contained 17 key candidate genes associated with economically important traits in pigs. Among these, nine candidate genes (CCNT2, EGR2, MYL3, CDH13, PROX1, FLVCR1, SETD2, FGF18, and FGF20) found in WECPs were related to muscular and skeletal development, whereas eight candidate genes (CSN1S1, SULT1E1, TJP1, ZNF366, LIPC, MCEE, STAP1, and DUSP) found in AHIPs were associated with health, reproduction, and fatness traits., Conclusion: Our findings provide a useful reference for the selection and assortative mating of pig breeds, laying the groundwork for future research on the population genetic structures of AHIPs, ultimately helping protect these local varieties., (© 2022. The Author(s).)

Background: The success of marker assisted selection depends on the amount of linkage disequilibrium (LD) across the genome. To implement marker assisted selection in the swine breeding industry, information about extent and degree of LD is essential. The objective of this study is to estimate LD in four US breeds of pigs (Duroc, Hampshire, Landrace, and Yorkshire) and subsequently calculate persistence of phase among them using a 60 k SNP panel. In addition, we report LD when using only a fraction of the available markers, to estimate persistence of LD over distance. Results: Average r² between adjacent SNP across all chromosomes was 0.36 for Landrace, 0.39 for Yorkshire, 0.44 for Hampshire and 0.46 for Duroc. For markers 1 Mb apart, r² ranged from 0.15 for Landrace to 0.20 for Hampshire. Reducing the marker panel to 10% of its original density, average r² ranged between 0.20 for Landrace to 0.25 for Duroc. We also estimated persistence of phase as a measure of prediction reliability of markers in one breed by those in another and found that markers less than 10 kb apart could be predicted with a maximal accuracy of 0.92 for Landrace with Yorkshire. Conclusions: Our estimates of LD, although in good agreement with previous reports, are more comprehensive and based on a larger panel of markers. Our estimates also confirmed earlier findings reporting higher LD in pigs than in American Holstein cattle, especially at increasing marker distances (> 1 Mb). High average LD (r² > 0.4) between adjacent SNP found in this study is an important precursor for the implementation of marker assisted selection within a livestock species. Results of this study are relevant to the US purebred pig industry and critical for the design of programs of whole genome marker assisted evaluation and selection. In addition, results indicate that a more cost efficient implementation of marker assisted selection using low density panels with genotype imputation, would be... [ABSTRACT FROM AUTHOR]

Background: Male infertility is an increasing problem in all domestic species including man. Localization and identification of genes involved in defects causing male infertility provide valuable information of specific events in sperm development. Sperm development is a complex process, where diploid spermatogonia develop into haploid, highly specialized spermatozoa. Correct expression and function of various genes and their protein products are required for production of fertile sperm. We have identified an infertility defect in Finnish Yorkshire boars caused by spermatogenic arrest. The aim of this study was to locate the disease associated region using genome wide screen with the PorcineSNP60 Beadchip and identify the causal mutation by candidate gene approach. Results: In the Finnish Yorkshire pig population the spermatogenic arrest (SA) defect appears to be of genetic origin and causes severe degeneration of germ cells and total absence of spermatozoa. Genome wide scan with the PorcineSNP60 Beadchip localized the SA defect to porcine chromosome 12 in a 2 Mbp region. Sequencing of a candidate gene Tex14 revealed a 51 bp insertion within exon 27, which caused differential splicing of the exon and created a premature translation stop codon. The expression of Tex14 was markedly down regulated in the testis of a SA affected boar compared to control boars and no protein product was identified by Western blotting. The SA insertion sequence was also found within intron 27 in all analyzed animals, thus the insertion appears to be a possible duplication event. Conclusion: In this study we report the identification of a causal mutation for infertility caused by spermatogenic arrest at an early meiotic phase. Our results highlight the role of TEX14 specifically in spermatogenesis and the importance of specific genomic remodeling events as causes for inherited defects. [ABSTRACT FROM AUTHOR]