Characterization of sexual maturity-associated N6-methyladenosine in boar testes.

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.

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

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

Reaction-norm animal model analysis of average daily gain heat tolerance in purebred Duroc pigs.

The present study aimed to genetically improve growth performance under high-heat environments by specifically designing a reaction-norm animal model (RNAM) for purebred Duroc pigs in Japan. A total of 54,750 records of average daily gain (ADG) measured for pigs reared at four farms in different prefectures were analyzed. Estimated maximum daily temperatures at the respective farm locations were used to calculate the average cumulative thermal load (TL). The TL values served as an indicator of high-heat environments for pigs. The plausible cumulative period length and threshold temperature for calculating TL were determined to be 8 weeks until just before shipping and 25°C, respectively. Variance components were estimated via RNAM analysis using TL as a linear covariate. The estimated additive genetic variances under both responsive and non-responsive to TL were found to be significant. Moreover, the estimated heritability of ADG ranged from 0.38 to 0.73 for TL values of 0-8. These results suggest that the RNAM developed holds the potential for improving the genetic ability of growth under high-heat environments in pigs.

Detection and Analysis of Antidiarrheal Genes and Immune Factors in Various Shanghai Pig Breeds.

The aim of this study was to identify effective genetic markers for the Antigen Processing Associated Transporter 1 (TAP1), α (1,2) Fucosyltransferase 1 (FUT1), Natural Resistance Associated Macrophage Protein 1 (NRAMP1), Mucin 4 (MUC4) and Mucin 13 (MUC13) diarrhea-resistance genes in the local pig breeds, namely Shanghai white pigs, Fengjing pigs, Shawutou pigs, Meishan pigs and Pudong white pigs, to provide a reference for the characterization of local pig breed resources in Shanghai. Polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLR) and sequence sequencing were applied to analyze the polymorphisms of the above genes and to explore the effects on the immunity of Shanghai local pig breeds in conjunction with some immunity factors. The results showed that both TAP1 and MUC4 genes had antidiarrheal genotype GG in the five pig breeds, AG and GG genotypes of the FUT1 gene were detected in Pudong white pigs, AA antidiarrheal genes of the NRAMP1 gene were detected in Meishan pigs, the AB type of the NRAMP1 gene was detected in Pudong white pigs, and antidiarrheal genotype GG of the MUC13 gene was only detected in Shanghai white pigs. The MUC13 antidiarrhea genotype GG was only detected in Shanghai white pigs. The TAP1 gene was moderately polymorphic in Shanghai white pigs, Fengjing pigs, Shawutou pigs, Meishan pigs and Pudong white pigs, among which TAP1 in Shanghai white pigs and Shawutou pigs did not satisfy the Hardy-Weinberg equilibrium. The FUT1 gene of Pudong white pigs was in a state of low polymorphism. NRAMP1 of Meishan pigs and Pudong white pigs was in a state of moderate polymorphism, which did not satisfy the Hardy-Weinberg equilibrium. The MUC4 genes of Shanghai white pigs and Pudong white pigs were in a state of low polymorphism, and the MUC4 genes of Fengjing pigs and Shawutou pigs were in a state of moderate polymorphism, and the MUC4 genes of Fengjing pigs and Pudong white pigs did not satisfy the Hardy-Weinberg equilibrium. The MUC13 gene of Shanghai white pigs and Pudong white pigs was in a state of moderate polymorphism. Meishan pigs had higher levels of IL-2, IL-10, IgG and TNF-α, and Pudong white pigs had higher levels of IL-12 than the other pigs. The level of interleukin 12 (IL-12) was significantly higher in the AA genotype of the MUC13 gene of Shanghai white pigs than in the AG genotype. The indicator of tumor necrosis factor alpha (TNF-α) in the AA genotype of the TAP1 gene of Fengjing pigs was significantly higher than that of the GG and AG genotypes. The indicator of IL-12 in the AG genotype of the Shawutou pig TAP1 gene was significantly higher than that of the GG genotype. The level of TNF-α in the AA genotype of the NRAMP1 gene of Meishan pigs was markedly higher than that of the AB genotype. The IL-2 level of the AG type of the FUT1 gene was obviously higher than that of the GG type of Pudong white pigs, the IL-2 level of the AA type of the MUC4 gene was dramatically higher than that of the AG type, and the IgG level of the GG type of the MUC13 gene was apparently higher than that of the AG type. The results of this study are of great significance in guiding the antidiarrhea breeding and molecular selection of Shanghai white pigs, Fengjing pigs, Shawutou pigs, Meishan pigs and Pudong white pigs and laying the foundation for future antidiarrhea breeding of various local pig breeds in Shanghai.

Integrated analysis of the DNA methylome and RNA transcriptome during the development of skeletal muscle in Duroc pigs.

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.

Genomic regions, candidate genes, and pleiotropic variants associated with physiological and anatomical indicators of heat stress response in lactating sows.

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.

Population dynamics of potentially harmful haplotypes: a pedigree analysis.

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.

Impact of the leptin receptor gene on pig performance and quality traits.

The recessive T allele of the missense polymorphism rs709596309 C > T of the leptin receptor gene is associated with intramuscular fat. However, its overall impact on pork production is still partial. In this work, we investigated the all-round effects of the TT genotype on lean growth efficiency and carcass, meat and fat quality using data from an experiment that compared the performance of 48 TT and 48 C- (24 CT and 24 CC) Duroc barrows. The TT pigs were less efficient for lean growth than the C- pigs. Although heavier, their carcasses had less lean content, were shorter and had lighter loins. Apart from increasing marbling and saturated fatty acid content, changes caused by the TT genotype in meat and fat quality are likely not enough to be perceived by consumers. The effect on visual marbling score exceeded that on intramuscular fat content, which suggests a direct influence of the T allele on the pattern of fat distribution in muscle. With current low-protein diets, the T allele is expected to be cost-effective only in niche markets where a very high level of marbling is critical.

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

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

Identification of Breed-Specific SNPs of Danish Large White Pig in Comparison with Four Chinese Local Pig Breed Genomes.

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

Whole-Transcriptome Analysis Sheds Light on the Biological Contexts of Intramuscular Fat Deposition in Ningxiang Pigs.

The quality of pork is significantly impacted by intramuscular fat (IMF). However, the regulatory mechanism of IMF depositions remains unclear. We performed whole-transcriptome sequencing of the longissimus dorsi muscle (IMF) from the high (5.1 ± 0.08) and low (2.9 ± 0.51) IMF groups (%) to elucidate potential mechanisms. In summary, 285 differentially expressed genes (DEGs), 14 differentially expressed miRNAs (DEMIs), 83 differentially expressed lncRNAs (DELs), and 79 differentially expressed circRNAs (DECs) were identified. DEGs were widely associated with IMF deposition and liposome differentiation. Furthermore, competing endogenous RNA (ceRNA) regulatory networks were constructed through co-differential expression analyses, which included circRNA-miRNA-mRNA (containing 6 DEMIs, 6 DEGs, 47 DECs) and lncRNA-miRNA-mRNA (containing 6 DEMIs, 6 DEGs, 36 DELs) regulatory networks. The circRNAs sus-TRPM7_0005, sus-MTUS1_0004, the lncRNAs SMSTRG.4269.1, and MSTRG.7983.2 regulate the expression of six lipid metabolism-related target genes, including PLCB1, BAD, and GADD45G, through the binding sites of 2-4068, miR-7134-3p, and miR-190a. For instance, MSTRG.4269.1 regulates its targets PLCB1 and BAD via miRNA 2_4068. Meanwhile, sus-TRPM7_0005 controls its target LRP5 through ssc-miR-7134-3P. These findings indicate molecular regulatory networks that could potentially be applied for the marker-assisted selection of IMF to enhance pork quality.

Gene expression supports a single origin of horns and antlers in hoofed mammals.

Horns, antlers, and other bony cranial appendages of even-toed hoofed mammals (ruminant artiodactyls) challenge traditional morphological homology assessments. Cranial appendages all share a permanent bone portion with family-specific integument coverings, but homology determination depends on whether the integument covering is an essential component or a secondary elaboration of each structure. To enhance morphological homology assessments, we tested whether juvenile cattle horn bud transcriptomes share homologous gene expression patterns with deer antlers relative to pig outgroup tissues, treating the integument covering as a secondary elaboration. We uncovered differentially expressed genes that support horn and antler homology, potentially distinguish them from non-cranial-appendage bone and other tissues, and highlight the importance of phylogenetic outgroups in homology assessments. Furthermore, we found differentially expressed genes that could support a shared cranial neural crest origin for horns and antlers and expression patterns that refine our understanding of the timing of horn and antler differentiation.

Detection and Analysis of Genes Affecting the Number of Thoracic Vertebrae in Licha Black Pigs.

The number of thoracic vertebrae (NTV) in pigs is an important economic trait that significantly influences pork production. While the Licha black pig is a well-known Chinese pig breed with multiple thoracic vertebrae, the genetic mechanism is still unknown. Here, we performed a selective signal analysis on the genome of Licha black pigs, comparing individuals with 15 NTV versus those with 16 NTV to better understand functional genes associated with NTV. A total of 2265 selection signal regions were detected across the genome, including 527 genes and 1073 QTL that overlapped with the selection signal regions. Functional enrichment analysis revealed that LRP5 and SP5 genes were involved in biological processes such as bone morphogenesis and Wnt protein binding. Furthermore, three genes, LRP8, DEF6, and SCUBE3, associated with osteoblast differentiation and bone formation, were located within or close to the QTL related to bone development and vertebrae number. These five genes were hypothesized to be potential candidates for regulating the NTV trait in Licha black pigs. Our findings revealed several candidate genes that play crucial roles in NTV regulation and provide a theoretical foundation to understand the genetic mechanism of the NTV trait in pig breeding.

Genetic diversity, growth and heart function of Auckland Island pigs, a potential source for organ xenotransplantation.

One of the prerequisites for successful organ xenotransplantation is a reasonable size match between the porcine organ and the recipient's organ to be replaced. Therefore, the selection of a suitable genetic background of source pigs is important. In this study, we investigated body and organ growth, cardiac function, and genetic diversity of a colony of Auckland Island pigs established at the Center for Innovative Medical Models (CiMM), LMU Munich. Male and female Auckland Island pig kidney cells (selected to be free of porcine endogenous retrovirus C) were imported from New Zealand, and founder animals were established by somatic cell nuclear transfer (SCNT). Morphologically, Auckland Island pigs have smaller body stature compared to many domestic pig breeds, rendering their organ dimensions well-suited for human transplantation. Furthermore, echocardiography assessments of Auckland Island pig hearts indicated normal structure and functioning across various age groups throughout the study. Single nucleotide polymorphism (SNP) analysis revealed higher runs of homozygosity (ROH) in Auckland Island pigs compared to other domestic pig breeds and demonstrated that the entire locus coding the swine leukocyte antigens (SLAs) was homozygous. Based on these findings, Auckland Island pigs represent a promising genetic background for organ xenotransplantation.

Effects of foetal size, sex and developmental stage on adaptive transcriptional responses of skeletal muscle to intrauterine growth restriction in pigs.

Intrauterine growth restriction (IUGR) occurs both in humans and domestic species. It has a particularly high incidence in pigs, and is a leading cause of neonatal morbidity and mortality as well as impaired postnatal growth. A key feature of IUGR is impaired muscle development, resulting in decreased meat quality. Understanding the developmental origins of IUGR, particularly at the molecular level, is important for developing effective strategies to mitigate its economic impact on the pig industry and animal welfare. The aim of this study was to characterise transcriptional profiles in the muscle of growth restricted pig foetuses at different gestational days (GD; gestational length ~ 115 days), focusing on selected genes (related to development, tissue injury and metabolism) that were previously identified as dysregulated in muscle of GD90 fetuses. Muscle samples were collected from the lightest foetus (L) and the sex-matched foetus with weight closest to the litter average (AW) from each of 22 Landrace x Large White litters corresponding to GD45 (n = 6), GD60 (n = 8) or GD90 (n = 8), followed by analyses, using RT-PCR and protein immunohistochemistry, of selected gene targets. Expression of the developmental genes, MYOD, RET and ACTN3 were markedly lower, whereas MSTN expression was higher, in the muscle of L relative to AW littermates beginning on GD45. Levels of all tissue injury-associated transcripts analysed (F5, PLG, KNG1, SELL, CCL16) were increased in L muscle on GD60 and, most prominently, on GD90. Among genes involved in metabolic regulation, KLB was expressed at higher levels in L than AW littermates beginning on GD60, whereas both IGFBP1 and AHSG were higher in L littermates on GD90 but only in males. Furthermore, the expression of genes specifically involved in lipid, hexose sugar or iron metabolism increased or, in the case of UCP3, decreased in L littermates on GD60 (UCP3, APOB, ALDOB) or GD90 (PNPLA3, TF), albeit in the case of ALDOB this only involved females. In conclusion, marked dysregulation of genes with critical roles in development in L foetuses can be observed from GD45, whereas for a majority of transcripts associated with tissue injury and metabolism differences between L and AW foetuses were apparent by GD60 or only at GD90, thus identifying different developmental windows for different types of adaptive responses to IUGR in the muscle of porcine foetuses.

Integrating Single-Cell and Spatial Transcriptomics Reveals Heterogeneity of Early Pig Skin Development and a Subpopulation with Hair Placode Formation.

The dermis and epidermis, crucial structural layers of the skin, encompass appendages, hair follicles (HFs), and intricate cellular heterogeneity. However, an integrated spatiotemporal transcriptomic atlas of embryonic skin has not yet been described and would be invaluable for studying skin-related diseases in humans. Here, single-cell and spatial transcriptomic analyses are performed on skin samples of normal and hairless fetal pigs across four developmental periods. The cross-species comparison of skin cells illustrated that the pig epidermis is more representative of the human epidermis than mice epidermis. Moreover, Phenome-wide association study analysis revealed that the conserved genes between pigs and humans are strongly associated with human skin-related diseases. In the epidermis, two lineage differentiation trajectories describe hair follicle (HF) morphogenesis and epidermal development. By comparing normal and hairless fetal pigs, it is found that the hair placode (Pc), the most characteristic initial structure in HFs, arises from progenitor-like OGN+/UCHL1+ cells. These progenitors appear earlier in development than the previously described early Pc cells and exhibit abnormal proliferation and migration during differentiation in hairless pigs. The study provides a valuable resource for in-depth insights into HF development, which may serve as a key reference atlas for studying human skin disease etiology using porcine models.

Constructing eRNA-mediated gene regulatory networks to explore the genetic basis of muscle and fat-relevant traits in pigs.

BACKGROUND: Enhancer RNAs (eRNAs) play a crucial role in transcriptional regulation. While significant progress has been made in understanding epigenetic regulation mediated by eRNAs, research on the construction of eRNA-mediated gene regulatory networks (eGRN) and the identification of critical network components that influence complex traits is lacking. RESULTS: Here, employing the pig as a model, we conducted a comprehensive study using H3K27ac histone ChIP-seq and RNA-seq data to construct eRNA expression profiles from multiple tissues of two distinct pig breeds, namely Enshi Black (ES) and Duroc. In addition to revealing the regulatory landscape of eRNAs at the tissue level, we developed an innovative network construction and refinement method by integrating RNA-seq, ChIP-seq, genome-wide association study (GWAS) signals and enhancer-modulating effects of single nucleotide polymorphisms (SNPs) measured by self-transcribing active regulatory region sequencing (STARR-seq) experiments. Using this approach, we unraveled eGRN that significantly influence the growth and development of muscle and fat tissues, and identified several novel genes that affect adipocyte differentiation in a cell line model. CONCLUSIONS: Our work not only provides novel insights into the genetic basis of economic pig traits, but also offers a generalizable approach to elucidate the eRNA-mediated transcriptional regulation underlying a wide spectrum of complex traits for diverse organisms.

Genomic evidence for human-mediated introgressive hybridization and selection in the developed breed.

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.

Genetic introgression from commercial European pigs to the indigenous Chinese Lijiang breed and associated changes in phenotypes.

BACKGROUND: Gene flow is crucial for enhancing economic traits of livestock. In China, breeders have used hybridization strategies for decades to improve livestock performance. Here, we performed whole-genome sequencing of a native Chinese Lijiang pig (LJP) breed. By integrating previously published data, we explored the genetic structure and introgression of genetic components from commercial European pigs (EP) into the LJP, and examined the impact of this introgression on phenotypic traits. RESULTS: Our analysis revealed significant introgression of EP breeds into the LJP and other domestic pig breeds in China. Using a haplotype-based approach, we quantified introgression levels and compared EP to LJP and other Chinese domestic pigs. The results show that EP introgression is widely prevalent in Chinese domestic pigs, although there are significant differences between breeds. We propose that LJP could potentially act as a mediator for the transmission of EP haplotypes. We also examined the correlation between EP introgression and the number of thoracic vertebrae in LJP and identified VRTN and STUM as candidate genes for this trait. CONCLUSIONS: Our study provides evidence of introgressed European haplotypes in the LJP breed and describes the potential role of EP introgression on phenotypic changes of this indigenous breed.