Unveiling the genetic structure of pig population in a Himalayan state Uttarakhand through microsatellite and mitochondrial DNA analyses.

This study traced the maternal lineage of the domestic swine populations using mitochondrial DNA control region markers and genetic diversity using microsatellite markers in Uttarakhand, an Indian state situated at the foothills of the world's youngest (geo-dynamically sensitive) mountain system, "the Himalayas". Analysis of 68 maternally unrelated individuals revealed 20 haplotypes. The maternal signature of the Pacific, Southeast Asian, European, and ubiquitously distributed Chinese haplotypes was present in Uttarakhand's domestic pig population. The D3 haplotype reported in wild pigs from North India was also identified in 47 domestic samples. A unique gene pool, UKD (Uttarakhand Domestic), as another lineage specific to this region has been proposed. Genotypes were analyzed, using 13 sets of microsatellite markers. The observed (Ho) and expected (He) heterozygosities were 0.83 ± 0.02 and 0.84 ± 0.01, respectively. The average polymorphic information content value of 0.83 ± 0.01 indicated the high informativeness of the marker. The overall mean FIS value for all the microsatellite markers was low (F = 0.04, P < 0.01). Seven loci deviated from Hardy-Weinberg equilibrium (HWE) at a significant level (p < 0.05). Two clusters were identified, indicating overlapping populations. These results suggested that though belonging to different maternal lineages, the traditional management practices in Uttarakhand have allowed for genetic mixing and the sharing of genetic material among pig populations. It could contribute to increased genetic diversity but might also result in the loss of distinct genetic characteristics or breed purity of the local breeds if not carefully managed.

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.

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.

Integrated Proteomics Analysis of Baseline Protein Expression in Pig Tissues.

The availability of an increasingly large amount of public proteomics data sets presents an opportunity for performing combined analyses to generate comprehensive organism-wide protein expression maps across different organisms and biological conditions. Sus scrofa, a domestic pig, is a model organism relevant for food production and for human biomedical research. Here, we reanalyzed 14 public proteomics data sets from the PRIDE database coming from pig tissues to assess baseline (without any biological perturbation) protein abundance in 14 organs, encompassing a total of 20 healthy tissues from 128 samples. The analysis involved the quantification of protein abundance in 599 mass spectrometry runs. We compared protein expression patterns among different pig organs and examined the distribution of proteins across these organs. Then, we studied how protein abundances were compared across different data sets and studied the tissue specificity of the detected proteins. Of particular interest, we conducted a comparative analysis of protein expression between pig and human tissues, revealing a high degree of correlation in protein expression among orthologs, particularly in brain, kidney, heart, and liver samples. We have integrated the protein expression results into the Expression Atlas resource for easy access and visualization of the protein expression data individually or alongside gene expression data.

Natural selection on feralization genes contributed to the invasive spread of wild pigs throughout the United States.

Despite a long presence in the contiguous United States (US), the distribution of invasive wild pigs (Sus scrofa × domesticus) has expanded rapidly since the 1980s, suggesting a more recent evolutionary shift towards greater invasiveness. Contemporary populations of wild pigs represent exoferal hybrid descendants of domestic pigs and European wild boar, with such hybridization expected to enrich genetic diversity and increase the adaptive potential of populations. Our objective was to characterize how genetic enrichment through hybridization increases the invasiveness of populations by identifying signals of selection and the ancestral origins of selected loci. Our study focused on invasive wild pigs within Great Smoky Mountains National Park, which represents a hybrid population descendent from the admixture of established populations of feral pigs and an introduction of European wild boar to North America. Accordingly, we genotyped 881 wild pigs with multiple high-density single-nucleotide polymorphism (SNP) arrays. We found 233 markers under putative selection spread over 79 regions across 16 out of 18 autosomes, which contained genes involved in traits affecting feralization. Among these, genes were found to be related to skull formation and neurogenesis, with two genes, TYRP1 and TYR, also encoding for crucial melanogenesis enzymes. The most common haplotypes associated with regions under selection for the Great Smoky Mountains population were also common among other populations throughout the region, indicating a key role of putatively selective variants in the fitness of invasive populations. Interestingly, many of these haplotypes were absent among European wild boar reference genotypes, indicating feralization through genetic adaptation.

Gut-associated microbes are present and active in the pig nasal cavity.

The nasal microbiota is a key contributor to animal health, and characterizing the nasal microbiota composition is an important step towards elucidating the role of its different members. Efforts to characterize the nasal microbiota composition of domestic pigs and other farm animals frequently report the presence of bacteria that are typically found in the gut, including many anaerobes from the Bacteroidales and Clostridiales orders. However, the in vivo role of these gut-microbiota associated taxa is currently unclear. Here, we tackled this issue by examining the prevalence, origin, and activity of these taxa in the nasal microbiota of piglets. First, analysis of the nasal microbiota of farm piglets sampled in this study, as well as various publicly available data sets, revealed that gut-microbiota associated taxa indeed constitute a substantial fraction of the pig nasal microbiota that is highly variable across individual animals. Second, comparison of herd-matched nasal and rectal samples at amplicon sequencing variant (ASV) level showed that these taxa are largely shared in the nasal and rectal microbiota, suggesting a common origin driven presumably by the transfer of fecal matter. Third, surgical sampling of the inner nasal tract showed that gut-microbiota associated taxa are found throughout the nasal cavity, indicating that these taxa do not stem from contaminations introduced during sampling with conventional nasal swabs. Finally, analysis of cDNA from the 16S rRNA gene in these nasal samples indicated that gut-microbiota associated taxa are indeed active in the pig nasal cavity. This study shows that gut-microbiota associated taxa are not only present, but also active, in the nasal cavity of domestic pigs, and paves the way for future efforts to elucidate the function of these taxa within the nasal microbiota.

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.

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.

A novel QTL region for pH and meat color in Duroc pigs.

One of the most important processes that occur during the transformation of muscle to meat is the pH decline as a consequence of the post-mortem metabolism of muscle tissue. Abnormal pH declines lead to pork defects such as pale, soft, and exudative meat. There is genetic variance for ultimate pH and the role of some genes on this phenotype is well established. After conducting a genome-wide association study on ultimate pH using 526 purebred Duroc pigs, we identified associated regions on Sus scrofa chromosomes (SSC) 3, 8, and 15. Functional candidate genes in these regions included PRKAG3 and PHKG1. The SSC8 region, at 71.6 Mb, was novel and, although no candidate causative gene could be identified, it may have regulatory effects. Subsequent analysis on 828 pigs from the same population confirmed the impact of the three associated regions on pH and meat color. We detected no interaction between the three regions. Further investigations are necessary to unravel the functional significance of the novel genomic region at SSC8. These variants could be used as markers in marker-assisted selection for improving meat quality.

Genome-wide association study identified five quantitative trait loci and two candidate genes for digestive traits in Suhuai pigs.

This work aimed to identify markers and candidate genes underlying porcine digestive traits. In total, 331 pigs were genotyped by 80 K Chip data or 50 K Chip data. For apparent neutral detergent fiber digestibility, a total of 19 and 21 candidate single nucleotide polymorphisms (SNP) were respectively identified using a genome-wide efficient mixed-model association algorithm and linkage-disequilibrium adjusted kinship. Among them, three quantitative trait locus (QTL) regions were identified. For apparent acid detergent fiber digestibility, a total of 16 and 17 SNPs were identified by these two methods, respectively. Of these, three QTL regions were also identified. Moreover, two candidate genes (MST1 and LATS1), which are functionally related to intestinal homeostasis and health, were detected near these significant SNPs. Taken together, our results could provide a basis for deeper research on digestive traits in pigs.

Glyphosate-induced changes in the expression of galanin and GALR1, GALR2 and GALR3 receptors in the porcine small intestine wall.

Glyphosate is the active ingredient of glyphosate-based herbicides and the most commonly used pesticide in the world. The goal of the present study was to verify whether low doses of glyphosate (equivalent to the environmental exposure) evoke changes in galanin expression in intramural neurons in the small intestine in pigs and to quantitatively determine changes in the level of galanin receptor encoding mRNA (GALR1, GALR2, GALR3) in the small intestine wall. The experiment was conducted on 15 sexually immature gilts divided into three study groups: control (C)-animals receiving empty gelatin capsules; experimental 1 (G1)-animals receiving a low dose of glyphosate (0.05 mg/kg b.w./day); experimental 2 (G2)-animals receiving a higher dose of glyphosate (0.5 mg/kg b.w./day) orally in gelatine capsules for 28 days. Glyphosate ingestion led to an increase in the number of GAL-like immunoreactive intramural neurons in the porcine small intestine. The results of RT-PCR showed a significant increase in the expression of mRNA, which encodes the GAL-receptors in the ileum, a decreased expression in the duodenum and no significant changes in the jejunum. Additionally, intoxication with glyphosate increased the expression of SOD2-encoding mRNA in the duodenum and decreased it in the jejunum and ileum, but it did not affect SOD1 expression. The results suggest that it may be a consequence of the cytotoxic and/or neurotoxic properties of glyphosate and/or its ability to induce oxidative stress.

Resistance to African swine fever virus among African domestic pigs appears to be associated with a distinct polymorphic signature in the RelA gene and upregulation of RelA transcription.

African swine fever virus (ASFV) is a highly contagious and fatal hemorrhagic disease of domestic pigs, which poses a major threat to the swine industry worldwide. Studies have shown that indigenous African pigs tolerate ASFV infection better than European pigs. The porcine v-rel avian reticuloendotheliosis viral oncogene homolog A (RelA) encoding a p65 kD protein, a major subunit of the NF-kB transcription factor, plays important roles in controlling both innate and adaptive immunity during infection with ASFV. In the present study, RelA genes from ASFV-surviving and symptomatic pigs were sequenced and found to contain polymorphisms revealing two discrete RelA amino acid sequences. One was found in the surviving pigs, and the other in symptomatic pigs. In total, 16 nonsynonymous SNPs (nsSNPs) resulting in codon changes were identified using bioinformatics software (SIFT and Polyphen v2) and web-based tools (MutPre and PredictSNP). Seven nsSNPs (P374-S, T448-S, P462-R, V464-P, Q478-H, L495-E, and P499-Q) were predicted to alter RelA protein function and stability, while 5 of these (P374-S, T448-S, P462-R, L495-E, and Q499-P) were predicted as disease-related SNPs.Additionally, the inflammatory cytokine levels of IFN-α, IL-10, and TNF-α at both the protein and the mRNA transcript levels were measured using ELISA and Real-Time PCR, respectively. The resulting data was used in correlation analysis to assess the association between cytokine levels and the RelA gene expression. Higher levels of IFN-α and detectable levels of IL-10 protein and RelA mRNA were observed in surviving pigs compared to healthy (non-infected). A positive correlation of IFN-α cytokine levels with RelA mRNA expression was also obtained. In conclusion, 7 polymorphic events in the coding region of the RelA gene may contribute to the tolerance of ASFV in pigs.

A chromosome-level genome of Chenghua pig provides new insights into the domestication and local adaptation of pigs.

As a country with abundant genetic resources of pigs, the domestication history of pigs in China and the adaptive evolution of Chinese pig breeds at different latitudes have rarely been elucidated at the genome-wide level. To fill this gap, we first assembled a high-quality chromosome-level genome of the Chenghua pig and used it as a benchmark to analyse the genomes of 272 samples from three genera of three continents. The divergence of the three species belonging to three genera, Phacochoerus africanus, Potamochoerus porcus, and Sus scrofa, was assessed. The introgression of pig breeds redefined that the migration routes were basically from southern China to central and southwestern China, then spread to eastern China, arrived in northern China, and finally reached Europe. The domestication of pigs in China occurred ∼12,000 years ago, earlier than the available Chinese archaeological domestication evidence. In addition, FBN1 and NR6A1 were identified in our study as candidate genes related to extreme skin thickness differences in Eurasian pig breeds and adaptive evolution at different latitudes in Chinese pig breeds, respectively. Our study provides a new resource for the pig genomic pool and refines our understanding of pig genetic diversity, domestication, migration, and adaptive evolution at different latitudes.

Full genome sequence analysis of African swine fever virus isolates from Cameroon.

African swine fever (ASF) is a devastating disease of domestic pigs that has spread across the globe since its introduction into Georgia in 2007. The etiological agent is a large double-stranded DNA virus with a genome of 170 to 180 kb in length depending on the isolate. Much of the differences in genome length between isolates are due to variations in the copy number of five different multigene families that are encoded in repetitive regions that are towards the termini of the covalently closed ends of the genome. Molecular epidemiology of African swine fever virus (ASFV) is primarily based on Sanger sequencing of a few conserved and variable regions, but due to the stability of the dsDNA genome changes in the variable regions occur relatively slowly. Observations in Europe and Asia have shown that changes in other genetic loci can occur and that this could be useful in molecular tracking. ASFV has been circulating in Western Africa for at least forty years. It is therefore reasonable to assume that changes may have accumulated in regions of the genome other than the standard targets over the years. At present only one full genome sequence is available for an isolate from Western Africa, that of a highly virulent isolate collected from Benin during an outbreak in 1997. In Cameroon, ASFV was first reported in 1981 and outbreaks have been reported to the present day and is considered endemic. Here we report three full genome sequences from Cameroon isolates of 1982, 1994 and 2018 outbreaks and identify novel single nucleotide polymorphisms and insertion-deletions that may prove useful for molecular epidemiology studies in Western Africa and beyond.

Duplex Droplet Digital PCR Assay for Quantification of Hepatitis E Virus in Food.

Hepatitis E virus (HEV) represents an emerging risk in industrialized countries where the consumption of contaminated food plays a pivotal role. Quantitative real-time RT-PCR (RT-qPCR) is one of the most suitable methods for the detection and quantification of viruses in food. Nevertheless, quantification using RT-qPCR has limitations. Droplet digital PCR (ddPCR) provides the precise quantification of nucleic acids without the need for a standard curve and a reduction in the effect on virus quantification due to the presence of inhibitors. The objectives of the present work were (i) to develop a method for the absolute quantification of HEV in swine tissues based on ddPCR technology and provide internal process control for recovery assessment and (ii) to evaluate the performance of the method by analyzing a selection of naturally contaminated wild boar muscle samples previously tested using RT-qPCR. The method was optimized using a set of in vitro synthesized HEV RNA and quantified dsDNA. The limit of detection of the developed ddPCR assay was 0.34 genome copies/µL. The analysis of the wild boar samples confirmed the validity of the ddPCR assay. The duplex ddPCR method showed no reduction in efficiency compared to individual assays. The method developed in the present study could represent a sensitive assay for the detection and absolute quantification of HEV RNA in food samples with the advantage of presenting the co-amplification of internal process control.

Elucidation of population stratifying markers and selective sweeps in crossbred Landlly pig population using genome-wide SNP data.

The present study was aimed at the identification of population stratifying markers from the commercial porcine SNP 60K array and elucidate the genome-wide selective sweeps in the crossbred Landlly pig population. Original genotyping data, generated on Landlly pigs, was merged in various combinations with global suid breeds that were grouped as exotic (global pig breeds excluding Indian and Chinese), Chinese (Chinese pig breeds only), and outgroup pig populations. Post quality control, the genome-wide SNPs were ranked for their stratifying power within each dataset in TRES (using three different criteria) and FIFS programs and top-ranked SNPs (0.5K, 1K, 2K, 3K, and 4K densities) were selected. PCA plots were used to assess the stratification power of low-density panels. Selective sweeps were elucidated in the Landlly population using intra- and inter-population haplotype statistics. Additionally, Tajima's D-statistics were calculated to determine the status of balancing selection in the Landlly population. PCA plots showed 0.5K marker density to effectively stratify Landlly from other pig populations. The A-score in DAPC program revealed the Delta statistic of marker selection to outperform other methods (informativeness and FST methods) and that 3000-marker density was suitable for stratification of Landlly animals from exotic pig populations. The results from selective sweep analysis revealed the Landlly population to be under selection for mammary (NAV2), reproductive efficiency (JMY, SERGEF, and MAP3K20), body conformation (FHIT, WNT2, ASRB, DMGDH, and BHMT), feed efficiency (CSRNP1 and ADRA1A), and immunity (U6, MYO3B, RBMS3, and FAM78B) traits. More than two methods suggested sweeps for immunity and feed efficiency traits, thus giving a strong indication for selection in this direction. The study is the first of its kind in Indian pig breeds with a comparison against global breeds. In conclusion, 500 markers were able to effectively stratify the breeds. Different traits under selective sweeps (natural or artificial selection) can be exploited for further improvement.

Towards identification of new genetic determinants for post-weaning diarrhea in piglets.

Post-weaning diarrhea in pigs is a considerable challenge in the pig farming industry due to its effect on animal welfare and production costs, as well as the large volume of antibiotics, which are used to treat diarrhea in pigs after weaning. Previous studies have revealed loci on SSC6 and SSC13 associated with susceptibility to specific diarrhea causing pathogens. This study aimed to identify new genetic loci for resistance to diarrhea based on phenotypic data. In depth clinical characterization of diarrhea was performed in 257 pigs belonging to two herds during the first 14 days post weaning. The daily diarrhea assessments were used for the classification of pigs into case and control groups. Pigs were assigned to case and control groups based only on the incidence of diarrhea in the second week of the study in order to differentiate between differences in etiology. Genome-wide association studies and metabolomics association analysis were performed in order to identify new biological determinants for diarrhea susceptibility. With the present work, we revealed a new locus for diarrhea resistance on SSC16. Furthermore, studies of metabolomics in the same pigs revealed one metabolite associated with diarrhea.

Phenotypic correlations of carpal gland diverticular number with production traits and its genome-wide association analysis in multiple pig populations.

Pig carpal glands play crucial roles in territorial recognition, reproductive behavior, and information exchange; however, their effects on production traits and underlying genetic mechanisms remain unclear. In this study, 1028 pigs from six populations were counted for the carpal gland diverticular numbers (CGDNs) on the left (CGDNL) and right (CGDNR) legs, and their carcass and meat quality traits were assessed. The CGDNs were significantly different among the populations, and Licha Black pigs had a lower CGDN than the Bama Xiang breed. It was also significantly different between sexes, with males having more diverticula than females (p ≤ 0.0391). Moreover, the number was asymmetric, with CGDNR being significantly higher than CGDNL. Notably, CGDNs was significantly correlated with each other in phenotype and genetics and with 24-h pH, 24-h meat color score, 24-h marbling score, fat content, moisture content, sodium salt content, and saturated fatty acid content in phenotype. Furthermore, genome-wide association analyses identified seven SNPs in association with CGDNs at a 5% genome-wide significance level, all of which were located in a 1.78-Mb (35.347-37.129 Mb) region on chromosome 1. CNC10010837 and CNC10010840 were the top SNPs: both had an additive effect of 0.789 ± 0.120 on CGDNR with p = 8.31E-10. These findings provide important insights into the functions and underlying genetic mechanisms of swine carpal glands.

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

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