Genome-wide association study identifying genetic variants associated with carcass backfat thickness, lean percentage and fat percentage in a four-way crossbred pig population using SLAF-seq technology

Genome-wide association studies reveal the genetic basis of growth and carcass traits in Sichuan Shelduck

China has abundant local duck resource populations, and evaluating the characteristics of these breeds will help improve development and utilization. In this study, we conducted the first investigations of growth and slaughter performance on Sichuan Shelduck (n = 240), an endangered duck local breed. The average body weight is 1497.91 g at 90 d of age. According to the growth curve through data recorded every 2 wk, we observed a low relative growth rate (RGR) for the early growth stage. The RGR shows a decreasing trend with age increasing in the stage from 0 to 56 d of age. The SNP-based heritability estimation showed the growth rate has a relatively high heritability, indicating high genetic stability for this trait. In the correlation analysis, the percentage of leg muscle is positively correlated with the absolute growth rate (AGR) at 28 to 42 d of age, whereas it is negatively correlated with the earlier stages, exhibiting a time-specific correlation result. Additionally, genome-wide association studies (GWAS) identified PCSK6, TOX2, and TOMM7 as potential candidate genes influencing AGR (42-56) and AGR (56-90), while the candidate genes of slaughter traits were PTP4A2, FAM110B, TOX, UBXN2B, and FCHSD2. These results provide an important reference for further understanding the genetic basis of growth and meat production performance of Sichuan Shelduck.

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

BACKGROUND

Carcass traits are essential economic traits in the commercial pig industry. However, the genetic mechanism of carcass traits is still unclear. In this study, we performed a genome-wide association study (GWAS) based on the specific-locus amplified fragment sequencing (SLAF-seq) to study seven carcass traits on 223 four-way intercross pigs, including dressing percentage (DP), number of ribs (RIB), skin thinkness (ST), carcass straight length (CSL), carcass diagonal length (CDL), loin eye width (LEW), and loin eye thickness (LET).

RESULTS

A total of 227,921 high-quality single nucleotide polymorphisms (SNPs) were detected to perform GWAS. A total of 30 SNPs were identified for seven carcass traits using the mixed linear model (MLM) (p < 1.0 × 10- 5), of which 9 SNPs were located in previously reported quantitative trait loci (QTL) regions. The phenotypic variation explained (PVE) by the significant SNPs was from 2.43 to 16.32%. Furthermore, 11 candidate genes (LYPLAL1, EPC1, MATN2, ZFAT, ZBTB10, ZNF704, INHBA, SMYD3, PAK1, SPTBN2, and ACTN3) were found for carcass traits in pigs.

CONCLUSIONS

The GWAS results will improve our understanding of the genetic basis of carcass traits. We hypothesized that the candidate genes associated with these discovered SNPs would offer a biological basis for enhancing the carcass quality of pigs in swine breeding.

Identification of molecular markers and candidate regions associated with grain number per spike in Pubing3228 using SLAF-BSA

Grain number per spike, a pivotal agronomic trait dictating wheat yield, lacks a comprehensive understanding of its underlying mechanism in Pubing3228, despite the identification of certain pertinent genes. Thus, our investigation sought to ascertain molecular markers and candidate regions associated with grain number per spike through a high-density genetic mapping approach that amalgamates site-specific amplified fragment sequencing (SLAF-seq) and bulked segregation analysis (BSA). To facilitate this, we conducted a comparative analysis of two wheat germplasms, Pubing3228 and Jing4839, known to exhibit marked discrepancies in spike shape. By leveraging this methodology, we successfully procured 2,810,474 SLAF tags, subsequently resulting in the identification of 187,489 single nucleotide polymorphisms (SNPs) between the parental strains. We subsequently employed the SNP-index association algorithm alongside the extended distribution (ED) association algorithm to detect regions associated with the trait. The former algorithm identified 24 trait-associated regions, whereas the latter yielded 70. Remarkably, the intersection of these two algorithms led to the identification of 25 trait-associated regions. Amongst these regions, we identified 399 annotated genes, including three genes harboring non-synonymous mutant SNP loci. Notably, the APETALA2 (AP2) transcription factor families, which exhibited a strong correlation with spike type, were also annotated. Given these findings, it is plausible to hypothesize that these genes play a critical role in determining spike shape. In summation, our study contributes significant insights into the genetic foundation of grain number per spike. The molecular markers and candidate regions we have identified can be readily employed for marker-assisted breeding endeavors, ultimately leading to the development of novel wheat cultivars possessing enhanced yield potential. Furthermore, conducting further functional analyses on the identified genes will undoubtedly facilitate a comprehensive elucidation of the underlying mechanisms governing spike development in wheat.

Comparison of differentially expressed genes in longissimus dorsi muscle of Diannan small ears, Wujin and landrace pigs using RNA-seq

Introduction

Pig growth is an important economic trait that involves the co-regulation of multiple genes and related signaling pathways. High-throughput sequencing has become a powerful technology for establishing the transcriptome profiles and can be used to screen genome-wide differentially expressed genes (DEGs). In order to elucidate the molecular mechanism underlying muscle growth, this study adopted RNA sequencing (RNA-seq) to identify and compare DEGs at the genetic level in the longissimus dorsi muscle (LDM) between two indigenous Chinese pig breeds (Diannan small ears [DSE] pig and Wujin pig [WJ]) and one introduced pig breed (Landrace pig [LP]).

Methods

Animals under study were from two Chinese indigenous pig breeds (DSE pig, n = 3; WJ pig, n = 3) and one introduced pig breed (LP, n = 3) were used for RNA sequencing (RNA-seq) to identify and compare the expression levels of DEGs in the LDM. Then, functional annotation, Gene Ontology (GO) enrichment analysis, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis, and Protein-Protein Interaction (PPI) network analysis were performed on these DEGs. Then, functional annotation, Gene Ontology (GO) enrichment analysis, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis, and Protein-Protein Interaction (PPI) network analysis were performed on these DEGs.

Results

The results revealed that for the DSE, WJ, and LP libraries, more than 66, 65, and 71 million clean reads were generated by transcriptome sequencing, respectively. A total of 11,213 genes were identified in the LDM tissue of these pig breeds, of which 7,127 were co-expressed in the muscle tissue of the three samples. In total, 441 and 339 DEGs were identified between DSE vs. WJ and LP vs. DSE in the study, with 254, 193 up-regulated genes and 187, 193 down-regulated genes in DSE compared to WJ and LP. GO analysis and KEGG signaling pathway analysis showed that DEGs are significantly related to contractile fiber, sarcolemma, and dystrophin-associated glycoprotein complex, myofibril, sarcolemma, and myosin II complex, Glycolysis/Gluconeogenesis, Propanoate metabolism, and Pyruvate metabolism, etc. In combination with functional annotation of DEGs, key genes such as ENO3 and JUN were identified by PPI network analysis.

Discussion

In conclusion, the present study revealed key genes including DES, FLNC, PSMD1, PSMD6, PSME4, PSMB4, RPL11, RPL13A, ROS23, RPS29, MYH1, MYL9, MYL12B, TPM1, TPM4, ENO3, PGK1, PKM2, GPI, and the unannotated new gene ENSSSCG00000020769 and related signaling pathways that influence the difference in muscle growth and could provide a theoretical basis for improving pig muscle growth traits in the future.

Identifying Genetic Architecture of Carcass and Meat Quality Traits in a Ningxiang Indigenous Pig Population

Ningxiang pig is a breed renowned for its exceptional meat quality, but it possesses suboptimal carcass traits. To elucidate the genetic architecture of meat quality and carcass traits in Ningxiang pigs, we assessed heritability and executed a genome-wide association study (GWAS) concerning carcass length, backfat thickness, meat color parameters (L.LD, a.LD, b.LD), and pH at two postmortem intervals (45 min and 24 h) within a Ningxiang pig population. Heritability estimates ranged from moderate to high (0.30~0.80) for carcass traits and from low to high (0.11~0.48) for meat quality traits. We identified 21 significant SNPs, the majority of which were situated within previously documented QTL regions. Furthermore, the GRM4 gene emerged as a pleiotropic gene that correlated with carcass length and backfat thickness. The ADGRF1, FKBP5, and PRIM2 genes were associated with carcass length, while the NIPBL gene was linked to backfat thickness. These genes hold the potential for use in selective breeding programs targeting carcass traits in Ningxiang pigs.

Genome-Wide Association Analysis Identifies Genomic Regions and Candidate Genes for Growth and Fatness Traits in Diannan Small-Ear (DSE) Pigs

In the livestock industry, the growth and fatness traits are directly related to production efficiency and economic profits. As for Diannan small-ear (DSE) pigs, a unique indigenous breed, the genetic architecture of growth and fatness traits is still elusive. The aim of this study was to search the genetic loci and candidate genes associated with phenotypic traits in DSE pigs using GWAS based on the Geneseek Porcine 50K SNP Chip data. A total of 22,146 single nucleotide polymorphisms (SNPs) were detected in 265 DSE pigs and used for Genome-wide association studies (GWAS) analysis. Seven SNPs were found to be associated with back height, chest circumference, cannon bone circumference, and backfat thickness at the suggestive significance level. Based on gene annotation results, these seven SNPs were, respectively, mapped to the following candidate genes, VIPR2, SLC10A2, NUCKS1, MCT1, CHCHD3, SMOX, and GPR1, which are mainly involved with adipocyte differentiation, lipid metabolism, skeletal muscle development, and average daily weight gain. Our work offers novel insights into the genetic architecture of economically important traits in swine and may play an important role in breeding using molecular markers in the DSE breed.

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.

Identification of New Genes and Genetic Variant Loci Associated with Breast Muscle Development in the Mini-Cobb F2 Chicken Population Using a Genome-Wide Association Study

Native chicken has become a favorite choice for consumers in many Asian countries recently, not only for its potential nutritional value but also for its deep ties to local food culture. However, low growth performance and limited meat production restrict their economic potential. Conducting a genome-wide association study (GWAS) for chicken-breast muscle development will help identify loci or candidate genes for different traits and potentially provide new insight into this phenotype in chickens and other species. To improve native chicken growth performance, especially breast muscle development, we performed a GWAS to explore the potential genetic mechanisms of breast muscle development in an F2 population constructed by reciprocal crosses between a fast-growing broiler chicken (Cobb500) and a slow-growing native chicken (Daweishan mini chicken). The results showed that 11 SNPs, which exceeded the 10% genome significance level (p = 1.79 × 10-8) were considered associated with breast muscle development traits, where six SNPS, NC_006126.5: g.3138376T>G, NC_006126.5: g.3138452A>G, NC_006088.5: g.73837197A>G, NC_006088.5: g.159574275A>G, NC_006089.5: g.80832197A>G, and NC_006127.5: g.48759869G>T was first identified in this study. In total, 13 genes near the SNPs were chosen as candidate genes, and none of them had previously been studied for their role in breast muscle development. After grouping the F2 population according to partial SNPs, significant differences in breast muscle weight were found among different genotypes (p < 0.05), and the expression levels of ALOX5AP, USPL1, CHRNA9, and EFNA5 among candidate genes were also significantly different (p < 0.05). The results of this study will contribute to the future exploration of the potential genetic mechanisms of breast muscle development in domestic chickens and also support the expansion of the market for native chicken in the world.

Genome-Wide Association Study of Growth Traits in a Four-Way Crossbred Pig Population

Growth traits are crucial economic traits in the commercial pig industry and have a substantial impact on pig production. However, the genetic mechanism of growth traits is not very clear. In this study, we performed a genome-wide association study (GWAS) based on the specific-locus amplified fragment sequencing (SLAF-seq) to analyze ten growth traits on 223 four-way intercross pigs. A total of 227,921 highly consistent single nucleotide polymorphisms (SNPs) uniformly dispersed throughout the entire genome were used to conduct GWAS. A total of 53 SNPs were identified for ten growth traits using the mixed linear model (MLM), of which 18 SNPs were located in previously reported quantitative trait loci (QTL) regions. Two novel QTLs on SSC4 and SSC7 were related to average daily gain from 30 to 60 kg (ADG30–60) and body length (BL), respectively. Furthermore, 13 candidate genes (ATP5O, GHRHR, TRIM55, EIF2AK1, PLEKHA1, BRAP, COL11A2, HMGA1, NHLRC1, SGSM1, NFATC2, MAML1, and PSD3) were found to be associated with growth traits in pigs. The GWAS findings will enhance our comprehension of the genetic architecture of growth traits. We suggested that these detected SNPs and corresponding candidate genes might provide a biological foundation for improving the growth and production performance of pigs in swine breeding.