High-throughput genome sequencing and genotyping have significantly accelerated genetic research. However, the high cost of whole-genome sequencing (WGS) remains a barrier to large-scale studies like genome-wide association studies (GWAS) and genomic prediction. Genotype imputation offers a cost-effective alternative by inferring unobserved variants from lower-density data using haplotype reference panels. In this study, we present the updated Pig Haplotype Reference Panel (PHARP) 4.0, comprising 6449 pig genomes from 154 breeds. PHARP 4.0 encompasses 50.3 million SNPs and 5.8 million indels, making it the largest and most diverse pig reference panel to date. PHARP 4.0 demonstrated superior imputation accuracy compared to existing panels (SWIM, AHC, AGIDB, and PGRP), achieving concordance rates (CR > 0.99) and correlation coefficients (R² > 0.98) in European breeds and improved accuracy in Chinese Jinhua pigs (CR = 0.936, R² = 0.924) when imputing from 80 K SNP chip data to whole-genome sequencing (WGS). We further optimized an RNA-seq-based imputation pipeline by incorporating multiple breeds and applying a 6× sequencing depth filter, achieving CR > 0.95 and R² > 0.90 in European breeds, and a CR of 0.93 with an R² = 0.92 in Chinese Jinhua pigs. Additionally, increasing the specific reference panel size to approximately 400 samples improved the imputation of rare variants. Utilizing PHARP 4.0, we successfully imputed low-density SNP chip data for two GWAS, identifying significant SNPs likely representing causal variants. Overall, PHARP 4.0 serves as a valuable resource for advancing pig genetic research and supporting breeding programs. PHARP 4.0 is an updated pig haplotype reference panel with 6449 genomes from 154 breeds. It demonstrates superior imputation accuracy from chip to WGS data, enables an optimized RNA-seq imputation pipeline, and successfully identifies novel causal variants in GWAS.
Single nucleotide polymorphisms (SNPs) have been widely used in marker-assisted selection to demonstrate the molecular mechanisms underlying hereditary traits by screening candidate genes. Low-coverage whole genome resequencing (lcWGS) followed by genotype imputation represents an efficient strategy for acquiring abundant SNPs with a low cost. Although this strategy has been successfully implemented in animal husbandry, its optimal applicability in aquaculture breeds still need to be explored. Here, we investigated the optimal prephasing and imputation workflow in Paralichthys olivaceus, an economically important bony fish that requires urgent germplasm improvement. A total of 10.3 million high-quality biallelic SNPs were identified and 90 % sites with higher confidence were phased as the reference haplotype panel from whole genome resequencing (18.63x) of 168 P. olivaceus. It is revealed that the combination of prephasing by SHAPEIT and imputation through GLIMPSE1 demonstrated best performance across all three tested chromosomes. Furthermore, 0.5x coverage was identified as the most cost-effective resequencing coverage. Following optimizing an important imputation parameter “effective population number (ne)”, the median Pearson correlation coefficient reached 90 % at 0.5x. Moreover, a group of lcWGS P. olivaceus (1x) was imputed and added to original 168 individuals for population genetic analysis. Distinct genetic differentiation and faster linkage disequilibrium decay were observed in all cultured populations. Strong selection signals were detected on a 3.59 Mb region with 111 candidate genes on chromosome 14 in a cultured breed. Our study provides valuable practical insights for applying genotype imputation in aquaculture species. The identification of SNPs and candidate genes has greatly enriched the germplasm resource of P. olivaceus, offering candidate markers for molecular breeding and providing valuable resources for its preservation and management.
ObjectiveTo address the issues of uneven residual compressive stress distribution and high surface roughness in gears after shot peening, this study investigates the effect of double shot peening on the surface integrity of 18CrNiMo7-6 carburized gear steel.MethodsA shot peening model with random shot distribution was established using Abaqus software. The effects of shot peening parameters such as shot velocity and shot diameter on the residual stress and surface morphology of the 18CrNiMo7-6 carburized gear steel were analyzed for both single and double shot peening.ResultsThe results show that compared with single shot peening, double shot peening leads to more uniform plastic deformation on the material surface. The double shot peening significantly increases the residual compressive stress in the surface layer and shifts the depth position of the maximum residual compressive stress towards the surface, but has no significant effect on the depth of the residual compressive stress layer. Through variance analysis, when the shot coverage is 100%, using a lower shot peening intensity for the second shot peening process improves the uniformity of the surface residual compressive stress more effectively, with the maximum improvement reaching 28.76%.
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.
Pigs not only function as a major meat source worldwide but also are commonly used as an animal model for studying human complex traits. A large haplotype reference panel has been used to facilitate efficient phasing and imputation of relatively sparse genome-wide microarray chips and low-coverage sequencing data. Using the imputed genotypes in the downstream analysis, such as GWASs, TWASs, eQTL mapping and genomic prediction (GS), is beneficial for obtaining novel findings. However, currently, there is still a lack of publicly available and high-quality pig reference panels with large sample sizes and high diversity, which greatly limits the application of genotype imputation in pigs. In response, we built the pig Haplotype Reference Panel (PHARP) database. PHARP provides a reference panel of 2012 pig haplotypes at 34 million SNPs constructed using whole-genome sequence data from more than 49 studies of 71 pig breeds. It also provides Web-based analytical tools that allow researchers to carry out phasing and imputation consistently and efficiently. PHARP is freely accessible at http://alphaindex.zju.edu.cn/PHARP/index.php . We demonstrate its applicability for pig commercial 50 K SNP arrays, by accurately imputing 2.6 billion genotypes at a concordance rate value of 0.971 in 81 Large White pigs (~ 17 × sequencing coverage). We also applied our reference panel to impute the low-density SNP chip into the high-density data for three GWASs and found novel significantly associated SNPs that might be casual variants.
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.
Background Carcass backfat thickness (BFT), carcass lean percentage (CLP) and carcass fat percentage (CFP) are important to the commercial pig industry. Nevertheless, the genetic architecture of BFT, CLP and CFP is still elusive. Here, we performed a genome-wide association study (GWAS) based on specific-locus amplified fragment sequencing (SLAF-seq) to analyze seven fatness-related traits, including five BFTs, CLP, and CFP on 223 four-way crossbred pigs. Results A total of 227, 921 highly consistent single nucleotide polymorphisms (SNPs) evenly distributed throughout the genome were used to perform GWAS. Using the mixed linear model (MLM), a total of 20 SNP loci significantly related to these traits were identified on ten Sus scrofa chromosomes (SSC), of which 10 SNPs were located in previously reported quantitative trait loci (QTL) regions. On SSC7, two SNPs (SSC7:29,503,670 and rs1112937671) for average backfat thickness (ABFT) exceeded 1% and 10% Bonferroni genome-wide significance levels, respectively. These two SNP loci were located within an intron region of the COL21A1 gene, which was a protein-coding gene that played an important role in the porcine backfat deposition by affecting extracellular matrix (ECM) remodeling. In addition, based on the other three significant SNPs on SSC7, five candidate genes, ZNF184 , ZNF391 , HMGA1 , GRM4 and NUDT3 were proposed to influence BFT. On SSC9, two SNPs for backfat thickness at 6–7 ribs (67RBFT) and one SNP for CLP were in the same locus region (19 kb interval). These three SNPs were located in the PGM2L1 gene, which encoded a protein that played an indispensable role in glycogen metabolism, glycolysis and gluconeogenesis as a key enzyme. Finally, one significant SNP on SSC14 for CLP was located within the PLBD2 gene, which participated in the lipid catabolic process. Conclusions A total of two regions on SSC7 and SSC9 and eight potential candidate genes were found for fatness-related traits in pigs. The results of this GWAS based on SLAF-seq will greatly advance our understanding of the genetic architecture of BFT, CLP, and CFP traits. These identified SNP loci and candidate genes might serve as a biological basis for improving the important fatness-related traits of pigs.
Abstract Background Over several decades, a wide range of natural and artificial selection events in response to subtropical environments, intensive pasture and intensive feedlot systems have greatly changed the customary behaviour, appearance, and important economic traits of Shanghai Holstein cattle. In particular, the longevity of the Shanghai Holstein cattle population is generally short, approximately the 2nd to 3rd lactation. In this study, two complementary approaches, integrated haplotype score (iHS) and runs of homozygosity (ROH), were applied for the detection of selection signatures within the genome using genotyping by genome-reduced sequence data from 1092 cows. Results In total, 101 significant iHS genomic regions containing selection signatures encompassing a total of 256 candidate genes were detected. There were 27 significant |iHS| genomic regions with a mean |iHS| score > 2. The average number of ROH per individual was 42.15 ± 25.47, with an average size of 2.95 Mb. The length of 78 % of the detected ROH was within the range of 1–2 MB and 2–4 MB, and 99 % were shorter than 8 Mb. A total of 168 genes were detected in 18 ROH islands (top 1 %) across 16 autosomes, in which each SNP showed a percentage of occurrence > 30 %. There were 160 and 167 genes associated with the 52 candidate regions within health-related QTL intervals and 59 candidate regions within reproduction-related QTL intervals, respectively. Annotation of the regions harbouring clustered |iHS| signals and candidate regions for ROH revealed a panel of interesting candidate genes associated with adaptation and economic traits, such as IL22RA1, CALHM3, ITGA9, NDUFB3, RGS3, SOD2, SNRPA1, ST3GAL4, ALAD, EXOSC10, and MASP2. In a further step, a total of 1472 SNPs in 256 genes were matched with 352 cis-eQTLs in 21 tissues and 27 trans-eQTLs in 6 tissues. For SNPs located in candidate regions for ROH, a total of 108 cis-eQTLs in 13 tissues and 4 trans-eQTLs were found for 1092 SNPs. Eighty-one eGenes were significantly expressed in at least one tissue relevant to a trait (P value < 0.05) and matched the 256 genes detected by iHS. For the 168 significant genes detected by ROH, 47 gene-tissue pairs were significantly associated with at least one of the 37 traits. Conclusions We provide a comprehensive overview of selection signatures in Shanghai Holstein cattle genomes by combining iHS and ROH. Our study provides a list of genes associated with immunity, reproduction and adaptation. For functional annotation, the cGTEx resource was used to interpret SNP-trait associations. The results may facilitate the identification of genes relevant to important economic traits and can help us better understand the biological processes and mechanisms affected by strong ongoing natural or artificial selection in livestock populations.
Laiwu pigs, distinguished by their high intramuscular fat of 7–9%, is an indigenous pig breed of China, and recent studies also found that Laiwu pigs showed high resistance to Porcine circovirus type 2. However, with the introduction of commercial varieties, the population of Laiwu pigs has declined, and some lineages have disappeared, which could result in inbreeding. Runs of homozygosity (ROH) can be used as a good measure of individual inbreeding status and is also normally used to detect selection signatures so as to map the candidate genes associated with economically important traits. In this study, we used data from Genotyping by Genome Reducing and Sequencing to investigate the number, length, coverage, and distribution patterns of ROH in 93 Chinese Laiwu pigs and identified genomic regions with a high ROH frequency. The average inbreeding coefficient calculated by pedigree was 0.021, whereas that estimated by all detected ROH segments was 0.133. Covering 13.4% of the whole genome, a total of 7,508 ROH segments longer than 1 Mb were detected, whose average length was 3.76 Mb, and short segments (1–5 Mb) dominated. For individuals, the coverage was in the range between 0.56 and 36.86%. For chromosomes, SSC6 had the largest number ( n = 688), and the number of ROH in SSC12 was the lowest ( n = 215). Thirteen ROH islands were detected in our study, and 86 genes were found within those regions. Some of these genes were correlated with economically important traits, such as meat quality ( ECI1 , LRP12 , NDUFA4L2 , GIL1 , and LYZ ), immunity capacity ( IL23A, STAT2 , STAT6 , TBK1 , IFNG , and ITH2 ), production ( DCSTAMP , RDH16 , and GDF11 ), and reproduction ( ODF1 and CDK2 ). A total of six significant Gene Ontology terms and nine significant Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were identified, most of which were correlated with disease resistance and biosynthesis processes, and one KEGG pathway was related to lipid metabolism. In addition, we aligned all of the ROH islands to the pig quantitative trait loci (QTL) database and finally found eight QTL related to the intramuscular fat trait. These results may help us understand the characteristics of Laiwu pigs and provide insight for future breeding strategies.
Background China is the country with the most abundant swine genetic resources in the world. Through thousands of years of domestication and natural selection, most of pigs in China have developed unique genetic characteristics. Finding the unique genetic characteristics and modules of each breed is an essential part of their precise conservation. Results In this study, we used the partial least squares method to identify the significant specific SNPs of 19 local Chinese pig breeds and 5 Western pig breeds. A total of 37,514 significant specific SNPs ( p < 0.01) were obtained from these breeds, and the Chinese local pig breed with the most significant SNPs was Hongdenglong (HD), followed by Jiaxing black (JX), Huaibei (HB), Bihu (BH), small Meishan (SMS), Shengxian Hua (SH), Jiangquhai (JQ), Mi (MI), Chunan (CA), Chalu (CL), Jinhualiangtouwu (JHL), Fengjing (FJ), middle Meishan (MMS), Shanzhu (SZ), Pudong white (PD), Dongchuan (DC), Erhualian (EH), Shawutou (SW) and Lanxi Hua (LX) pig. Furthermore, we identified the breeds with the most significant genes, GO terms, pathways, and networks using KOBAS and IPA and then ranked them separately. The results showed that the breeds with the highest number of interaction networks were Hongdenglong (12) and Huaibei (12) pigs. In contrast, the breeds with the lowest interaction networks were Shawutou (4) and Lanxi Hua pigs (3), indicating that Hongdenglong and Huaibei pigs might have the most significant genetic modules in their genome, whereas Shawutou and Lanxi Hua pigs may have the least unique characteristics. To some degree, the identified specific pathways and networks are related to the number of genes and SNPs linked to the specific breeds, but they do not appear to be the same. Most importantly, more significant modules were found to be related to the development and function of the digestive system, regulation of diseases, and metabolism of amino acids in the local Chinese pig breeds, whereas more significant modules were found to be related to the growth rate in the Western pig breeds. Conclusion Our results show that each breed has some relatively unique structural modules and functional characteristics. These modules allow us to better understand the genetic differences among local Chinese and Western pig breeds and therefore implement precise conservation methods. This study could provide a basis for formulating more effective strategies for managing and protecting these genetic resources in the future.
The aim of this study was to detect evidence for signatures of recent selection in the Jinhua pig genome. These results can be useful to better understand the regions under selection in Jinhua pigs and might shed some lights on groups of genes that control production traits. In the present study, we performed extended haplotype homozygosity (EHH) tests to identify significant core regions in 202 Jinhua pigs. A total of 26 161 core regions spanning 636.42 Mb were identified, which occupied approximately 28% of the genome across all autosomes, and 1 158 significant (P<0.01) core haplotypes were selected. Genes in these regions were related to several economically important traits, including meat quality, reproduction, immune responses and exterior traits. A panel of genes including ssc-mir-365-2, KDM8, RABEP2, GSG1L, RHEB, RPH3AL and a signal pathway of PI3K-Akt were detected with the most extreme P-values. The findings in our study could draw a comparatively genome-wide map of selection signature in the pig genome, and also help to detect functional candidate genes under positive selection for further genetic and breeding research in Jinhua and other pigs.
Genetic characterization of Chinese indigenous pig breeds is essential to promote scientific conservation and sustainable development of pigs. Here, we systematically surveyed the genomes of 75 unrelated Diannan small-ear (DSE) pigs from three diverse regions (Yingjiang County, Jinping County, and Sipsongpanna in Yunnan Province) to describe their population structures, genetic diversity, inbreeding coefficients, and selection signatures. First, these individuals were sequenced and genotyped using the genome reducing and sequencing (GGRS) protocol. A total of 438,038 autosomal single-nucleotide polymorphisms (SNPs) were obtained and used for subsequent statistical analysis. The results showed that these DSE pigs were clearly differentiated into three separate clades revealed by the population structure and principal component analysis, which is consistent with their geographical origins. Diannan small-ear pigs owned lower genetic diversity when compared with some other pig breeds, which demonstrated the need to strengthen the conservation strategies for DSE pigs. In addition, the inbreeding coefficients based on runs of homozygosity (ROH) length (F-ROH) were calculated in each ROH length categories, respectively. And the results indicated that the ancient (up to 50 generations ago) inbreeding had greater impacts than recent (within the last five generations) inbreeding within DSE pigs. Some candidate selection signatures within the DSE pig population were detected through the ROH islands and integrated haplotype homozygosity score (iHS) methods. And genes associated with meat quality (COL15A1,RPL3L, andSLC9A3R2), body size (PALM2-AKAP2,NANS,TRAF7, andPACSIN1), adaptability (CLDN9andE4F1), and appetite (GRM4) were identified. These findings can help to understand the genetic characteristics and provide insights into the molecular background of special phenotypes of DSE pigs to promote conservation and sustainability of the breed.
The Jiangquhai (JQ) pig breed is one of the most widely recognized pig populations in China due to its unique and dominant characteristics. In this study, we examined the extent of Linkage disequilibrium (LD) and haplotype block structure of the JQ pig breed, and scanned the blocks for possible genes underlying important QTLs that could either be responsible for some adaptive features in these pigs or might have undergone some selection pressure. We compared some of our results with other Chinese and Western pig breeds. The results show that the JQ breed had the highest total block length (349.73 Mb ≈ 15% of its genome), and the coverage rate of blocks in most of its chromosomes was larger than those of other breeds except for Sus scrofa chromosome 4 (SSC4), SSC6, SSC7, SSC8, SSC10, SSC12, SSC13, SSC14, SSC17, SSC18, and SSCX. Moreover, the JQ breed had more SNPs that were clustered into haplotype blocks than the other breeds examined in this study. Our shared and unique haplotype block analysis revealed that the Hongdenglong (HD) breed had the lowest percentage of shared haplotype blocks while the Shanzhu (SZ) breed had the highest. We found that the JQ breed had an average r2 > 0.2 at SNPs distances 10–20 kb and concluded that about 120,000–240,000 SNPs would be needed for a successful GWAS in the breed. Finally, we detected a total of 88 genes harbored by selected haplotype blocks in the JQ breed, of which only 4 were significantly enriched (p-value ≤ 0.05). These genes were significantly enriched in 2 GO terms (p-value < 0.01), and 2 KEGG pathways (p-value < 0.02). Most of these enriched genes were related to health. Also, most of the overlapping QTLs detected in the haplotype blocks were related to meat and carcass quality, as well as health, with a few of them relating to reproduction and production. These results provide insights into the genetic architecture of some adaptive and meat quality traits observed in the JQ pig breed and also revealed the pattern of LD in the genome of the pig. Our result provides significant guidance for improving the statistical power of GWAS and optimizing the conservation strategy for this JQ pig breed.
Objective: Porcine respiratory disease is one of the most important health problems causing significant economic losses. To understand the genetic basis for susceptibility to swine enzootic pneumonia (EP) in pigs, we detected 102,809 single nucleotide polymorphisms in a total of 249 individuals based on genome-wide sequencing data. Methods: Genome comparison of susceptibility to swine EP in three pig breeds (Jinhua, Erhualian, and Meishan) with two western lines that are considered more resistant (Duroc and Landrace) using cross-population extended haplotype homozygosity and F-statistic (F-ST) statistical approaches identified 691 positively selected genes. Based on quantitative trait loci, gene ontology terms and literature search, we selected 14 candidate genes that have convincible biological functions associated with swine EP or human asthma. Results: Most of these genes were tested by several methods including transcription analysis and candidate genes association study. Among these genes: cytochrome P450 1A1 and catenin beta 1 (CTNNB1) are involved in fertility; transforming growth factor beta receptor 3 plays a role in meat quality traits; Wnt family member 2, CTNNB1 and transcription factor 7 take part in adipogenesis and fat deposition simultaneously; plasminogen activator, urokinase receptor (completely linked to AXL receptor tyrosine kinase, r(2) = 1) plays an essential role in the successful ovulation of matured oocytes in pigs; colipase like 2 (strongly linked to SAM pointed domain containing ETS transcription factor, r(2) = 0.848) is involved in male fertility. Conclusion: These adverse genes susceptible to swine EP may be selected while selecting for economic traits (especially reproduction traits) due to pleiotropic and hitchhiking effect of linked genes. Our study provided a completely new point of view to understand the genetic basis for susceptibility or resistance to swine EP in pigs thereby, provides insight for designing sustainable breed selection programs. Finally, the candidate genes are crucial due to their potential roles in respiratory diseases in a large number of species, including human.
Chinese indigenous pigs in the Taihu Lake region are well known for their high fecundity and other excellent characteristics. To better understand the characteristics of these breeds in this area as well as to provide the government and breeders the molecular basis for formulating a reasonable conservation policy, we explored the structure of haplotype blocks and genetic diversity of the 7 populations which is relevant for the management and conservation of these important genetic resources using next-generation sequencing data. In this study, a total of 131 300 single-nucleotide polymorphisms with minor allele frequencies ⩾0.05 were obtained for further analysis. In general, there are similar within-breed genetic diversities (He, Ho, Pn, Ar) among these 7 pig populations in the Taihu Lake region. Average values for the inbreeding coefficients estimates in the 7 populations are 0.110 (F1), 0.056 (F2), and 0.078 (F3). All the breeds have seen a continuous decline in Ne estimates over time with FJ and SW populations having a very similar curve. Moreover, the Ne of SMS pig breeds were smaller than other Chinese pig breeds, indicating that SMS pig breeds underwent stronger selection pressure than other Chinese pig breeds. The average genetic distances among the 7 populations in the Taihu Lake region were 0.235 (MMS), 0.240 (SMS), 0.269 (EH), 0.248 (MI), 0.221 (FJ), 0.254 (JX), and 0.212 (SW). A summary of the number of haplotype blocks and haplotype diversity was also presented. This study provide a deep understanding of the current situation of conservation in this region, thereby uncovering the pertinent insight to better formulate more reasonable preservation policies for the government departments and breeding planners to follow-up.
本研究旨在探讨杜洛克(D)、大白(Y)和长白(L)3个知名西方品种与金华猪(J)的杂交利用情况.应用"基于基因组简化与测序的基因型判型(GGRS)"技术平台,检测单核苷酸多态性(SNPs)并判定基因型.利用全基因组性状特异的SNP分子标记估算不同群体间的遗传距离,从胴体与肉品质、繁殖、健康和生长4个性状对金华猪不同杂交模式下的杂种优势进行预测.结果表明:共找到繁殖、生长、胴体与肉品质和健康性状潜在相关基因的SNP分别为1835、2174、196057、17200个;在二元杂交组合中,♂杜 × ♀金(D×J)为较优组合;在三元杂交组合中,♂金 × ♀长大(J×LY)为较优组合;在四元杂交组合中,♂杜金 × ♀长大(DJ×LY)为较优组合.本研究应用生物信息学技术对金华猪不同杂交模式的杂交优势进行预测,可大大减少后续的配合力实验,结果可为金华猪的保种和杂交利用提供一定参考.
Jinhua pig, a well-known Chinese indigenous breed, has evolved as a pig breed with excellent meat quality, greater disease resistance, and higher prolificacy. The reduction in the number of Jinhua pigs over the past years has raised concerns about inbreeding. Runs of homozygosity (ROH) along the genome have been applied to quantify individual autozygosity to improve the understanding of inbreeding depression and identify genes associated with traits of interest. Here, we investigated the occurrence and distribution of ROH using next-generation sequencing data to characterize autozygosity in 202 Jinhua pigs, as well as to identify the genomic regions with high ROH frequencies within individuals. The average inbreeding coefficient, based on ROH longer than 1 Mb, was 0.168 +/- 0.052. In total, 18,690 ROH were identified in all individuals, among which shorter segments (1-5 Mb) predominated. Individual ROH autosome coverage ranged from 5.32 to 29.14% in the Jinhua population. On average, approximately 16.8% of the whole genome was covered by ROH segments, with the lowest coverage on SSC11 and the highest coverage on SSC17. A total of 824 SNPs (about 0.5%) and 11 ROH island regions were identified (occurring in over 45% of the samples). Genes associated with reproduction (HOXA3, HOXA7, HOXA10, and HOXA11), meat quality (MYOD1, LPIN3, and CTNNBL1), appetite (NUCB2) and disease resistance traits (MUC4, MUC13, MUC20, LMLN, ITGB5, HEG1, SLC12A8, and MYLK) were identified in ROH islands. Moreover, several quantitative trait loci for ham weight and ham fat thickness were detected. Genes in ROH islands suggested, at least partially, a selection for economic traits and environmental adaptation, and should be subject of future investigation. These findings contribute to the understanding of the effects of environmental and artificial selection in shaping the distribution of functional variants in the pig genome.
To improve development and utilization of Jinhua pig (JH) resources,we designed the crossings of Duroc × J H (D×J),Duroc × Landrace × JH (D×LJ),JH × Landrace × Yorkshire (J×LY),Duroc × JH × Landrace × Yorkshire (DJ×LY) to see how they affected the growth performance,reproductive performance,slaughter performance and meat quality,respectively,of Jinhua pigs.The results indicated that DJ×LY crossbreeding produced higher total litter size,alive litter size and newborn litter weight than the other three crossbreeding designs;D×LJ and DJ×LY had better dressing percentage and lean percentage;J×LY resulted in the highest carcass length,followed by DJ×LY cross;and DJ×LY led to the highest correct color grading and marbling scores.Taken together,all the JH crossbreeding trials in this study could serve as candidate crosses,with the DJ×LY crossbreedding being a more recommended one,for Jinhua pig production.
太湖猪(Sus scrofa)作为我国著名的地方猪品种,具有繁殖力高、耐粗饲、抗病性好、肉质佳等诸多特点,是分子育种和基因组选择的重要遗传资源。本研究基于猪数量性状基因座(quantitative trait locus,QTL)数据库、基因本体(gene ontology, GO)、通路分析(ingenuity pathway analysis, IPA)数据库等生物信息学技术对猪生长性状相关基因进行挖掘,同时筛选出其与西方品种存在品种间差异的特异SNPs,通过基因相关性,SNP特异性,SNP位置等信息综合预测可以通过基因编辑改良太湖地方品种的位点,并对猪基因组单导链RNA(single guide RNA, sgRNA)进行脱靶分析,筛选出CRISPR/Cas9编辑的最优sgRNA,旨在为太湖地方猪生长性状研究提供准确真实的基因编辑位点。在猪基因组中,共检测到127 353 502个前间区序列邻近基序(protospacer adjacent motif, PAM)位点,平均每20 bp存在1个PAM位点,共观察到5种PAM类型,包括TGG、AGG、GGG、CGG和NGG,占比分别为34.52%,32.61%,26.43%,6.43%和0.0082%。最终得到猪生长性状候选基因428个,筛选得到太湖流域猪种二花脸、枫泾、嘉兴黑、米猪、中梅山、沙乌头、小梅山的生长相关特异SNPs数量分别为412、282、372、414、393、279、458个,并为每一个位点综合打分,最后依据脱靶分析,为6个品种的每个生长相关SNP位点分别提供4个最优的sg RNA。本研究为后续高效准确的基因编辑工作提供数据基础。
Inbreeding, which has several causes including genetic drift, population bottlenecks, mating of close relatives and selection, can leave tracts of runs of homozygosity (ROH) along genomes. Recently, decreasing performance on reproductive traits, which might have resulted from inbreeding, has been observed in Chinese pigs. In this study, 830 individuals from Western and Chinese pig breeds were genotyped using the reduced-representation sequencing method. After imputation and quality control, 60 850 high-confidence SNPs were retained for ROH detection. A simulation was performed to explore the reliability of ROH detection with imputed data. Different ROH-related variables were compared between imputed and non-missing genotypes used in ROH detection. Furthermore, ROH islands were evaluated and annotated to find genes influenced by inbreeding in these pigs. The simulation results showed that imputed data with 0.7 as the average missing genotype rate and three heterozygotes allowed in a sliding window have comparable ROH detected compared with data with no missing genotypes. Compared with Western pig breeds, Chinese pigs had more autosomes covered by ROH longer than 5 Mb, indicating higher inbreeding in Chinese pigs in recent times. Genes related to reproduction, immunity, meat quality and adaptability in Chinese pigs and several genes related to growth speed and immunity in Western pigs were observed in short ROH islands. The reproduction-related gene PRM1 was found to be located in the most frequent long ROH island in Chinese pigs, which might explain the decreasing fertility in Chinese pig breeds.