Identifying Selection Signatures for Backfat Thickness in Yorkshire Pigs Highlights New Regions Affecting Fat Metabolism

Gender difference in the association of OSBPL8 polymorphisms with nephrolithiasis within a Chinese cohort

BACKGROUND

Kidney stone disease (KSD) is a common disorder of the urinary system and is closely related to genetic polymorphisms. However, the relationship between OSBPL8 polymorphisms and kidney stones has not been thoroughly investigated.

METHODS

Six OSBPL8 polymorphisms (rs17042391,rs17042409,rs4761431,rs7303892,rs4761434, and rs17042390) were analyzed in a Chinese case-control cohort containing 923 nephrolithiasis patients and 945 healthy controls.The association of these OSBPL8 gene polymorphisms with KSD susceptibility was analyzed using logistic regression, and examined by calculating the odds ratios (ORs) and 95% confidence intervals (CIs).

RESULTS

The OSBPL8 polymorphisms (GG for rs17042391, rs17042409, rs4761431, rs7303892; AA for rs4761434; and G for rs17042390) were significantly associated with a decreased risk of KSD in females.The protected alleles (G allele of rs17042391, G allele of rs17042409, G allele of rs4761431, A allele of rs4761434, and G allele of rs17042390) were related to decreased BMI levels in KSD patients; female patients with these alleles also exhibited lower BMI, HDL, and LDL levels,the G allele of rs7303892 was linked to reduced serum cholesterol levels in these females.Additionally, the haplotype ACAAGA was associated with decreased KSD risk in females, but haplotype GGGGAG presented an opposing effect.

CONCLUSION

Our research shows that the OSBPL8 gene polymorphisms reduced the risk of KSD in females, and were also associated with lipid-related metabolic traits.

A genome-wide association study identified candidate regions and genes for commercial traits in a Landrace population

Backfat thickness (BFT) and feed conversion ratio (FCR) are important commercial traits in the pig industry. With the increasing demand for human health and meat production, identifying functional genomic regions and genes associated with these commercial traits is critical for enhancing production efficiency. In this research, we conducted a genome-wide association study (GWAS) on a Landrace population comprising 4,295 individuals with chip data for BFT and FCR. Our analysis revealed a total of 118 genome-wide significant signals located on chromosomes SSC1, SSC2, SSC7, SSC12, and SSC13, respectively. Furthermore, we identified 10 potential regions associated with the two traits and annotated the genes within these regions. In addition, enrichment analysis was also performed. Notably, candidate genes such as SHANK2, KCNQ1, and ABL1 were found to be associated with BFT, whereas NAP1L4, LSP1, and PPFIA1 genes were related to the FCR. Our findings provide valuable insights into the genetic architecture of these two traits and offer guidance for future pig breeding efforts.

Characterization and Function Analysis of miRNA Editing during Fat Deposition in Chinese Indigenous Ningxiang Pigs

This study aimed to identify active miRNA editing sites during adipose development in Ningxiang pigs and analyze their characteristics and functions. Based on small RNA-seq data from the subcutaneous adipose tissues of Ningxiang pigs at four stages-30 days (piglet), 90 days (nursery), 150 days (early fattening), and 210 days (late fattening)-we constructed a developmental map of miRNA editing in the adipose tissues of Ningxiang pigs. A total of 505 miRNA editing sites were identified using the revised pipeline, with C-to-U editing types being the most prevalent, followed by U-to-C, A-to-G, and G-to-U. Importantly, these four types of miRNA editing exhibited base preferences. The number of editing sites showed obvious differences among age groups, with the highest occurrence of miRNA editing events observed at 90 days of age and the lowest at 150 days of age. A total of nine miRNA editing sites were identified in the miRNA seed region, with significant differences in editing levels (p < 0.05) located in ssc-miR-23a, ssc-miR-27a, ssc-miR-30b-5p, ssc-miR-15a, ssc-miR-497, ssc-miR-15b, and ssc-miR-425-5p, respectively. Target gene prediction and KEGG enrichment analyses indicated that the editing of miR-497 might potentially regulate fat deposition by inhibiting adipose synthesis via influencing target binding. These results provide new insights into the regulatory mechanism of pig fat deposition.

The Role of SOCS3 in Regulating Meat Quality in Jinhua Pigs

Meat quality is an important economic trait that influences the development of the pig industry. Skeletal muscle development and glycolytic potential (GP) are two crucial aspects that significantly impact meat quality. It has been reported that abnormal skeletal muscle development and high glycogen content results in low meat quality. However, the genetic mechanisms underlying these factors are still unclear. Compared with intensive pig breeds, Chinese indigenous pig breeds, such as the Jinhua pig, express superior meat quality characteristics. The differences in the meat quality traits between Jinhua and intensive pig breeds make them suitable for uncovering the genetic mechanisms that regulate meat quality traits. In this study, the Jinhua pig breed and five intensive pig breeds, including Duroc, Landrace, Yorkshire, Berkshire, and Pietrain pig breeds, were selected as experimental materials. First, the F_ST and XP-EHH methods were used to screen the selective signatures on the genome in the Jinhua population. Then, combined with RNA-Seq data, the study further confirmed that SOCS3 could be a key candidate gene that influences meat quality by mediating myoblast proliferation and glycometabolism because of the down-regulated expression of SOCS3 in Jinhua pigs compared with Landrace pigs. Finally, through SOCS3 knockout (KO) and overexpression (OE) experiments in mouse C2C12 cells, the results showed that SOCS3 regulated the cell proliferation of myoblasts. Moreover, SOCS3 is involved in regulating glucose uptake by the IRS1/PI3K/AKT signaling pathway. Overall, these findings provide a basis for the genetic improvement of meat quality traits in the pig industry.

The role of genetic selection and climatic factors in the dispersal of anatomically modern humans out of Africa

The evolutionarily recent dispersal of anatomically modern humans (AMH) out of Africa (OoA) and across Eurasia provides a unique opportunity to examine the impacts of genetic selection as humans adapted to multiple new environments. Analysis of ancient Eurasian genomic datasets (~1,000 to 45,000 y old) reveals signatures of strong selection, including at least 57 hard sweeps after the initial AMH movement OoA, which have been obscured in modern populations by extensive admixture during the Holocene. The spatiotemporal patterns of these hard sweeps provide a means to reconstruct early AMH population dispersals OoA. We identify a previously unsuspected extended period of genetic adaptation lasting ~30,000 y, potentially in the Arabian Peninsula area, prior to a major Neandertal genetic introgression and subsequent rapid dispersal across Eurasia as far as Australia. Consistent functional targets of selection initiated during this period, which we term the Arabian Standstill, include loci involved in the regulation of fat storage, neural development, skin physiology, and cilia function. Similar adaptive signatures are also evident in introgressed archaic hominin loci and modern Arctic human groups, and we suggest that this signal represents selection for cold adaptation. Surprisingly, many of the candidate selected loci across these groups appear to directly interact and coordinately regulate biological processes, with a number associated with major modern diseases including the ciliopathies, metabolic syndrome, and neurodegenerative disorders. This expands the potential for ancestral human adaptation to directly impact modern diseases, providing a platform for evolutionary medicine.

Genome-wide association study reveals genetic loci and candidate genes for meat quality traits in a four-way crossbred pig population

Meat quality traits (MQTs) have gained more attention from breeders due to their increasing economic value in the commercial pig industry. In this genome-wide association study (GWAS), 223 four-way intercross pigs were genotyped using the specific-locus amplified fragment sequencing (SLAF-seq) and phenotyped for PH at 45 min post mortem (PH45), meat color score (MC), marbling score (MA), water loss rate (WL), drip loss (DL) in the longissimus muscle, and cooking loss (CL) in the psoas major muscle. A total of 227, 921 filtered single nucleotide polymorphisms (SNPs) evenly distributed across the entire genome were detected to perform GWAS. A total of 64 SNPs were identified for six meat quality traits using the mixed linear model (MLM), of which 24 SNPs were located in previously reported QTL regions. The phenotypic variation explained (PVE) by the significant SNPs was from 2.43% to 16.32%. The genomic heritability estimates based on SNP for six meat-quality traits were low to moderate (0.07-0.47) being the lowest for CL and the highest for DL. A total of 30 genes located within 10 kb upstream or downstream of these significant SNPs were found. Furthermore, several candidate genes for MQTs were detected, including pH45 (GRM8), MC (ANKRD6), MA (MACROD2 and ABCG1), WL (TMEM50A), CL (PIP4K2A) and DL (CDYL2, CHL1, ABCA4, ZAG and SLC1A2). This study provided substantial new evidence for several candidate genes to participate in different pork quality traits. The identification of these SNPs and candidate genes provided a basis for molecular marker-assisted breeding and improvement of pork quality traits.

Integrating genome-wide association studies and population genomics analysis reveals the genetic architecture of growth and backfat traits in pigs

Growth and fat deposition are complex traits, which can affect economical income in the pig industry. Due to the intensive artificial selection, a significant genetic improvement has been observed for growth and fat deposition in pigs. Here, we first investigated genomic-wide association studies (GWAS) and population genomics (e.g., selection signature) to explore the genetic basis of such complex traits in two Large White pig lines (n = 3,727) with the GeneSeek GGP Porcine HD array (n = 50,915 SNPs). Ten genetic variants were identified to be associated with growth and fatness traits in two Large White pig lines from different genetic backgrounds by performing both within-population GWAS and cross-population GWAS analyses. These ten significant loci represented eight candidate genes, i.e., NRG4, BATF3, IRS2, ANO1, ANO9, RNF152, KCNQ5, and EYA2. One of them, ANO1 gene was simultaneously identified for both two lines in BF100 trait. Compared to single-population GWAS, cross-population GWAS was less effective for identifying SNPs with population-specific effect, but more powerful for detecting SNPs with population-shared effects. We further detected genomic regions specifically selected in each of two populations, but did not observe a significant enrichment for the heritability of growth and backfat traits in such regions. In summary, the candidate genes will provide an insight into the understanding of the genetic architecture of growth-related traits and backfat thickness, and may have a potential use in the genomic breeding programs in pigs.

Identifying selection signatures for litter size in Guangxi Bama Xiang pigs

Litter size is an important economic trait in pig production. However, the genetic mechanisms underlying varying litter size in Guangxi Bama Xiang pigs remain unknown. To identify selection signatures for litter size in Guangxi Bama Xiang pigs, we obtained 297 Illumina PorcineSNP50 BeadChip array data and the average born number (ABN) from parity one to nine in Guangxi Bama Xiang pigs. Fixation index (Fst) methods were used to identify the selection signature of the litter size, and three phenotypic gradient differential population pairs (according to the ABN) in individuals were used to reduce the false positives of signature selections. Single nucleotide polymorphisms (SNPs) were identified in the VEGFA promoter and exons. The general linear model was used to analyse the differences in distinct genotypes after they were typed using three-round multiplex PCR technology. Finally, the transcriptome factor and CpG island in the VEGFA promoter were predicted. A total of 328, 328 and 317 significant loci were identified in the 1st, 2nd and 3rd population pairs, respectively. After removing the false positives, 25 SNPs were defined as the selection signatures in relation to litter size. Ten (VEGFA, USP49, USP25, SRPK1, SLC26A8, RPL10A, PPARD, MAPK14, HMGA1 and CHRDL2) out of 52 genes in the selection regions were annotated as the candidate genes of litter size, respectively, VEGFA. There were no SNPs in the VEGFA exon region, but we obtained three SNPs (rs786889605, rs343769603 and rs323942424) in the VEGFA promoter regions. The ABN in CC was significantly higher than that in TT in rs786889605, and the ABN in TT was significantly lower than that in GG in rs323942424. Meanwhile, the mutation of the VEGFA promoter result in the loss of Sp1 and NF-1 and the formation of Oct-1. In summary, we obtained ten candidate genes, and two mutations in the VEGFA promoter that could be important potential molecular biomarkers for litter size in Bama Xiang pigs.

Identifying Candidate Genes for Short Gestation Length Trait in Chinese Qingping Pigs by Whole-Genome Resequencing and RNA Sequencing

Gestation length is a complex polygenic trait that affects pig fetal development. The Qingping (QP) pig, a Chinese native black pig breed, is characterized by short gestation length. However, the genetic architecture of short gestation length is still not clear. The present study aimed to explore the genetic architecture of short gestation length in QP pigs. In this study, selective sweep analyses were performed to detect selective sweep signatures for short gestation length traits between 100 QP pigs and 219 pigs from 15 other breeds. In addition, differentially expressed genes for the short gestation length between QP pigs and Large White pigs were detected by RNA sequencing. Comparing candidate genes from these methods with known genes for preterm birth in the database, we obtained 111 candidate genes that were known preterm birth genes. Prioritizing other candidate genes, 839 novel prioritized candidate genes were found to have significant functional similarity to preterm birth genes. In particular, we highlighted EGFR, which was the most prioritized novel candidate relative to preterm birth genes. Experimental validations in placental and porcine trophectoderm cells suggest that EGFR is highly expressed in the QP pigs with short gestation length and could regulate the NF-κΒ pathway and downstream expression of PTGS2. These findings comprehensively identified candidate genes for short gestation length trait at the genomic and transcriptomic levels. These candidate genes provide an important new resource for further investigation and genetic improvement of gestation length.

Trait-specific Selection Signature Detection Reveals Novel Loci of Meat Quality in Large White Pigs

In past decades, meat quality traits have been shaped by human-driven selection in the process of genetic improvement programs. Exploring the potential genetic basis of artificial selection and mapping functional candidate genes for economic traits are of great significance in genetic improvement of pigs. In this study, we focus on investigating the genetic basis of five meat quality traits, including intramuscular fat content (IMF), drip loss, water binding capacity, pH at 45 min (pH45min), and ultimate pH (pH24h). Through making phenotypic gradient differential population pairs, Wright’s fixation index (F_ST) and the cross-population extended haplotype homozogysity (XPEHH) were applied to detect selection signatures for these five traits. Finally, a total of 427 and 307 trait-specific selection signatures were revealed by F_ST and XPEHH, respectively. Further bioinformatics analysis indicates that some genes, such as USF1, NDUFS2, PIGM, IGSF8, CASQ1, and ACBD6, overlapping with the trait-specific selection signatures are responsible for the phenotypes including fat metabolism and muscle development. Among them, a series of promising trait-specific selection signatures that were detected in the high IMF subpopulation are located in the region of 93544042-95179724bp on SSC4, and the genes harboring in this region are all related to lipids and muscle development. Overall, these candidate genes of meat quality traits identified in this analysis may provide some fundamental information for further exploring the genetic basis of this complex trait.

Genomic Analyses Revealed the Genetic Difference and Potential Selection Genes of Growth Traits in Two Duroc Lines

Duroc pigs are famous for their high growth rate, feed conversion efficiency, and lean meat percentage. Given that they have been subjected to artificial selection and breeding in multiple countries, various lines with obvious differences in production performance have formed. In this study, we genotyped 3,770 American Duroc (AD) pigs and 2,098 Canadian Duroc (CD) pigs using the GeneSeek Porcine SNP50 Beadchip to dissect the genetic differences and potential selection genes of growth traits in these two Duroc pig lines. Population structure detection showed that there were significant genetic differences between the two Duroc pig lines. Hence, we performed F _ST and cross-population extended haplotype homozygosity (XP-EHH) analyses between the two lines. As a result, we identified 38 annotated genes that were significantly enriched in the gland development pathway in the AD line, and 61 annotated genes that were significantly enriched in the immune-related pathway in the CD line. For three growth traits including backfat thickness (BFT), loin muscle depth (LMD), and loin muscle area (LMA), we then performed selection signature detection at 5 and 10% levels within the line and identified different selected regions and a series of candidate genes that are involved in lipid metabolism and skeletal muscle development or repair, such as IRX3, EBF2, WNT10B, TLR2, PITX3, and SGCD. The differences in selected regions and genes between the two lines may be the cause of the differences in growth traits. Our study suggests significant genetic differences between the AD and CD lines, which provide a theoretical basis for selecting different Duroc lines as sires for different needs.

Whole Transcriptome Analysis Reveals a Potential Regulatory Mechanism of LncRNA-FNIP2/miR-24-3p/FNIP2 Axis in Chicken Adipogenesis

Lipid biosynthesis is a complex process, which is regulated by multiple factors including lncRNA. However, the role of lncRNA in chicken abdominal fat accumulation is still unclear. In this research, we collected liver tissues from six high abdominal fat rate Sanhuang broilers and six low abdominal fat rate Sanhuang broilers to perform lncRNA sequencing and small RNA sequencing. A total of 2,265 lncRNAs, 245 miRNAs, and 5,315 mRNAs were differently expressed. Among of them, 1,136 differently expressed genes were enriched in the metabolic process. A total of 36 differently expressed genes, which were considered as differently expressed lncRNAs' targets, were enriched in the metabolic process. In addition, we also found out that eight differently expressed miRNAs could target 19 differently expressed genes. FNIP2 and PEX5L were shared in a cis-regulatory network and a differently expressed miRNA target relationship network. LncRNA-FNIP2/miR-24-3p/FNIP2 axis was considered as a potential candidate that may participate in lipid synthesis. Experimentally, the objective reality of lncRNA-FNIP2/miR-24-3p/FNIP2 axis was clarified and the regulation effect of lncRNA-FNIP2/miR-24-3p/FNIP2 axis on synthesis was validated. In brief, our study reveals a potential novel regulatory mechanism that lncRNA-FNIP2/miR-24-3p/FNIP2 axis was considered as being involved in lipid synthesis during chicken adipogenesis in liver.

Detection of selection signatures for response to Aleutian mink disease virus infection in American mink

Aleutian disease (AD) is the most significant health issue for farmed American mink. The objective of this study was to identify the genomic regions subjected to selection for response to infection with Aleutian mink disease virus (AMDV) in American mink using genotyping by sequencing (GBS) data. A total of 225 black mink were inoculated with AMDV and genotyped using a GBS assay based on the sequencing of ApeKI-digested libraries. Five AD-characterized phenotypes were used to assign animals to pairwise groups. Signatures of selection were detected using integrated measurement of fixation index (F_ST) and nucleotide diversity (θπ), that were validated by haplotype-based (hap-FLK) test. The total of 99 putatively selected regions harbouring 63 genes were detected in different groups. The gene ontology revealed numerous genes related to immune response (e.g. TRAF3IP2, WDR7, SWAP70, CBFB, and GPR65), liver development (e.g. SULF2, SRSF5) and reproduction process (e.g. FBXO5, CatSperβ, CATSPER4, and IGF2R). The hapFLK test supported two strongly selected regions that contained five candidate genes related to immune response, virus–host interaction, reproduction and liver regeneration. This study provided the first map of putative selection signals of response to AMDV infection in American mink, bringing new insights into genomic regions controlling the AD phenotypes.

The effect of divergent selection for intramuscular fat on the domestic rabbit genome

An experiment of divergent selection for intramuscular fat was carried out at Universitat Politècnica de València. The high response of selection in intramuscular fat content, after nine generations of selection, and a multidimensional scaling analysis showed a high degree of genomic differentiation between the two divergent populations. Therefore, local genomic differences could link genomic regions, encompassing selective sweeps, to the trait used as selection criterion. In this sense, the aim of this study was to identify genomic regions related to intramuscular fat through three methods for detection of selection signatures and to generate a list of candidate genes. The methods implemented in this study were Wright's fixation index, cross population composite likelihood ratio and cross population - extended haplotype homozygosity. Genomic data came from the 9th generation of the two populations divergently selected, 237 from Low line and 240 from High line. A high single nucleotide polymorphism (SNP) density array, Affymetrix Axiom OrcunSNP Array (around 200k SNPs), was used for genotyping samples. Several genomic regions distributed along rabbit chromosomes (OCU) were identified as signatures of selection (SNPs having a value above cut-off of 1%) within each method. In contrast, 8 genomic regions, harbouring 80 SNPs (OCU1, OCU3, OCU6, OCU7, OCU16 and OCU17), were identified by at least 2 methods and none by the 3 methods. In general, our results suggest that intramuscular fat selection influenced multiple genomic regions which can be a consequence of either only selection effect or the combined effect of selection and genetic drift. In addition, 73 genes were retrieved from the 8 selection signatures. After functional and enrichment analyses, the main genes into the selection signatures linked to energy, fatty acids, carbohydrates and lipid metabolic processes were ACER2, PLIN2, DENND4C, RPS6, RRAGA (OCU1), ST8SIA6, VIM (OCU16), RORA, GANC and PLA2G4B (OCU17). This genomic scan is the first study using rabbits from a divergent selection experiment. Our results pointed out a large polygenic component of the intramuscular fat content. Besides, promising positional candidate genes would be analysed in further studies in order to bear out their contributions to this trait and their feasible implications for rabbit breeding programmes.