Background: Genome-wide association studies (GWAS) are an effective way to explore genotype–phenotype associations in humans, animals, and plants. Various GWAS methods have been developed based on different genetic or statistical assumptions. However, no single method is optimal for all traits and, for many traits, the putative single nucleotide polymorphisms (SNPs) that are detected by the different methods do not entirely overlap due to the diversity of the genetic architecture of complex traits. Therefore, multi-tool-based GWAS strategies that combine different methods have been increasingly employed. To take this one step further, we propose an ensemble-like GWAS strategy (E-GWAS) that statistically integrates GWAS results from different single GWAS methods. Results: E-GWAS was compared with various single GWAS methods using simulated phenotype traits with different genetic architectures. E-GWAS performed stably across traits with different genetic architectures and effectively controlled the number of false positive genetic variants detected without decreasing the number of true positive variants. In addition, its performance could be further improved by using a bin-merged strategy and the addition of more distinct single GWAS methods. Our results show that the numbers of true and false positive SNPs detected by the E-GWAS strategy slightly increased and decreased, respectively, with increasing bin size and when the number and the diversity of individual GWAS methods that were integrated in E-GWAS increased, the latter being more effective than the bin-merged strategy. The E-GWAS strategy was also applied to a real dataset to study backfat thickness in a pig population, and 10 candidate genes related to this trait and expressed in adipose-associated tissues were identified. Conclusions: Using both simulated and real datasets, we show that E-GWAS is a reliable and robust strategy that effectively integrates the GWAS results of different methods and reduces the number of false positive SNPs without decreasing that of true positive SNPs.

Appropriate growth curves can reflect more sophisticated growth patterns of animals than body weight, and thus, the identification of genes and variants related to the growth curve parameter traits contributes to revealing the fine growth and development characteristics of livestock. However, the ability of single genome-wide association analysis (GWAS) and transcriptome analyses to identify valuable genes and variants is limited. In this study, based on genome and transcriptome data, the growth curve parameter traits of hybrid pigs were analyzed, and a set of genes and variants were identified. The Gompertz–Laird growth curve model was optimized to reveal the growth pattern of F₂ individuals of Duroc × Erhualian pigs over four time points. Five growth parameters were estimated, including initial body weight ${(W}_0)$, instantaneous growth rate per day (L), coefficient of relative growth or maturing index (k), body weight at inflection point ${(W}_i)$, and average growth rate (GR). These five parameters were subjected to a genome-wide association study, differential gene expression analysis, and weighted gene co-expression network analysis (WGCNA). In the study, 336 pigs were genotyped, and 39,494 SNP markers were used for each pig in the analysis. Thirty of these pigs were also included in the transcriptomics analysis. Based on genome and transcriptome data, the integrated analyses identified five putative SNPs (including INRA0056566 on chromosome X, DRGA0004151 on chromosome 3, INRA0056460 on chromosome X, H3GA0049324 on chromosome 17, and H3GA0037747 on chromosome 13) and 15 candidate genes (PDGFA, VEGFD, CSPP1, EFHC1, PIK3C3, ZZZ3, GCC2, MAPK14, ZPR1, ISG15, ANG, CEBPD, ZHX3, CTBP2, and MYNN). The functional analysis indicated that these candidate genes played important roles in cell division and differentiation, development and aging, and skeletal muscle and fat formation. Our results provide insight into the genetic mechanisms underlying the growth and development of hybrid pigs and offer a theoretical basis for genomic breeding.

The development of multi-omics has increased the likelihood of further improving genomic prediction (GP) of complex traits. Gene expression data can directly reflect the genotype effect, and thus, they are widely used for GP. Generally, the gene expression data are integrated into multiple random effect models as independent data layers or used to replace genotype data for genomic prediction. In this study, we integrated pedigree, genotype, and gene expression data into the single-step method and investigated the effects of this integration on prediction accuracy. The integrated single-step method improved the genomic prediction accuracy of more than 90% of the 54 traits in the Duroc × Erhualian F₂ pig population dataset. On average, the prediction accuracy of the single-step method integrating gene expression data was 20.6% and 11.8% higher than that of the pedigree-based best linear unbiased prediction (ABLUP) and genome-based best linear unbiased prediction (GBLUP) when the weighting factor (w) was set as 0, and it was 5.3% higher than that of the single-step best linear unbiased prediction (ssBLUP) under different w values. Overall, the analyses confirmed that the integration of gene expression data into a single-step method could effectively improve genomic prediction accuracy. Our findings enrich the application of multi-omics data to genomic prediction and provide a valuable reference for integrating multi-omics data into the genomic prediction model.

Heterosis has been extensively used for pig genetic breeding and production, but the genetic basis of heterosis remains largely elusive. Crossbreeding between commercial and native breeds provides a good model to parse the genetic basis of heterosis. This study uses Duhua hybrid pigs, a crossbreed of Duroc and Liangguang small spotted pigs, as materials to explore the genetic basis underlying heterosis related to growth traits at the genomic level. The mid-parent heterosis (MPH) analysis showed heterosis of this Duhua offspring on growth traits. In this study, we examined the impact of additive and dominance effects on 100 AGE (age adjusted to 100 kg) and 100 BF (backfat thickness adjusted to 100 kg) of Duhua hybrid pigs. Meanwhile, we successfully identified SNPs associated with growth traits through both additive and dominance GWASs (genome-wide association studies). These findings will facilitate the subsequent in-depth studies of heterosis in the growth traits of Duhua pigs.

Understanding the pig immune function is crucial for disease-resistant breeding and potentially for human health research due to shared immune system features. Immune cell ratios, like monocyte/lymphocyte ratio (MLR) and neutrophil/lymphocyte ratio (NLR), offer a more comprehensive view of immune status compared to individual cell counts. However, research on pig immune cell ratios remains limited. This study investigated MLR and NLR in a Duroc × Erhualian F₂ resource population. Heritability analysis revealed high values (0.649 and 0.688 for MLR and NLR, respectively), suggesting a strong genetic component. Furthermore, we employed an ensemble-like GWAS (E-GWAS) strategy and functional annotation analysis to identify 11 MLR-associated and 6 NLR-associated candidate genes. These genes were significantly enriched in immune-related biological processes. These findings provide novel genetic markers and candidate genes associated with porcine immunity, thereby providing valuable insights for addressing biosecurity and animal welfare concerns in the pig industry.