Genomewide association studies for hematological traits and T lymphocyte subpopulations in a Duroc × Erhualian F resource population

Identification of SNPs and Candidate Genes Associated with Monocyte/Lymphocyte Ratio and Neutrophil/Lymphocyte Ratio in Duroc × Erhualian F2 Population

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 F2 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.

An atlas of expression quantitative trait loci of microRNAs in longissimus muscle of eight-way crossbred pigs

MicroRNAs (miRNAs) are key regulators of myocyte development and traits, yet insight into the genetic basis of variation in miRNA expression is still limited. Here, we present a systematic analysis of expression quantitative trait loci (eQTL) for miRNA profiling in longissimus muscle of pigs from an eight-breed crossed heterogeneous population. By integrating the whole-genome sequencing and miRNAomics data, we map 54 cis- and 292 trans-eQTLs at high resolution that are associated with the expression of 54 and 92 miRNAs, respectively. Twenty-three trans-acting loci are identified to affect the expression of nine myomiRs (known muscle-specific miRNAs). MiRNAs in mammalian conserved miRNA clusters are found to be subjected to regulation by shared cis-eQTLs, while the expression of mature miRNA-5p/-3p counterparts is more likely to be regulated by different cis-eQTLs. Fine mapping and bioinformatics analyses pinpoint the peak cis-eSNP of miR-4331-5p, rs344650810, which is located in its seed region, as a causal variant for the changes in expression and function of this miRNA. Additionally, rs344650810 is significantly (P < 0.01) correlated with the density and percentage of type I muscle fibers. Altogether, this study provides a comprehensive atlas of miRNA-eQTLs in porcine skeletal muscle and new insights into regulatory mechanisms of miRNA expression.

Genetic architecture of innate and adaptive immune cells in pigs

Pig industry is facing new challenges that make necessary to reorient breeding programs to produce more robust and resilient pig populations. The aim of the present work was to study the genetic determinism of lymphocyte subpopulations in the peripheral blood of pigs and identify genomic regions and biomarkers associated to them. For this purpose, we stained peripheral blood mononuclear cells to measure ten immune-cell-related traits including the relative abundance of different populations of lymphocytes, the proportions of CD4⁺ T cells and CD8⁺ T cells, and the ratio of CD4⁺/CD8⁺ T cells from 391 healthy Duroc piglets aged 8 weeks. Medium to high heritabilities were observed for the ten immune-cell-related traits and significant genetic correlations were obtained between the proportion of some lymphocytes populations. A genome-wide association study pointed out 32 SNPs located at four chromosomal regions on pig chromosomes SSC3, SSC5, SSC8, and SSCX as significantly associated to T-helper cells, memory T-helper cells and γδ T cells. Several genes previously identified in human association studies for the same or related traits were located in the associated regions, and were proposed as candidate genes to explain the variation of T cell populations such as CD4, CD8A, CD8B, KLRC2, RMND5A and VPS24. The transcriptome analysis of whole blood samples from animals with extreme proportions of γδ T, T-helper and memory T-helper cells identified differentially expressed genes (CAPG, TCF7L1, KLRD1 and CD4) located into the associated regions. In addition, differentially expressed genes specific of different T cells subpopulations were identified such as SOX13 and WC1 genes for γδ T cells. Our results enhance the knowledge about the genetic control of lymphocyte traits that could be considered to optimize the induction of immune responses to vaccines against pathogens. Furthermore, they open the possibility of applying effective selection programs for improving immunocompetence in pigs and support the use of the pig as a very reliable human biomedical model.

Assessment of Hematologic and Biochemical Parameters for Healthy Commercial Pigs in China

Hematologic and biochemical data are useful for indicating disease diagnosis and growth performance in swine. However, the assessment of these parameters in healthy commercial pigs is rare in China. Thus, blood samples were collected from 107 nursery pigs and 87 sows and were analyzed for 25 hematologic and 14 biochemical variables. After the rejection of the outliers and the detection of the data distribution, the correlations between the blood parameters were analyzed and the hematologic/biochemical RIs were preliminarily established using the 95% percentile RI. Correlation analysis showed that albumin was the hub parameter among the blood parameters investigated, and genes overlapping with key correlated variables were discovered. Most of the hematologic and biochemical parameters were significantly different between nursery pigs and sows. The 95% RIs of white blood cells and red blood cells were 7.18–24.52 × 109/L and 5.62–7.84 × 1012/L, respectively, for nursery pigs, but 9.34–23.84 × 109/L and 4.98–8.29 × 1012/L for sows. The 95% RIs of total protein and albumin were 43.16–61.23 g/dL and 19.35–37.86 g/dL, respectively, for nursery pigs, but 64.96–88.68 g/dL and 31.91–43.28 g/dL for sows. In conclusion, our study highlights the variability in blood parameters between nursery pigs and sows and provides fundamental data for the health monitoring of commercial pigs in China.

Assessment of DNA methylation in porcine immune cells reveals novel regulatory elements associated with cell-specific gene expression and immune capacity traits

Background

Genetics studies in the porcine immune system have enhanced selection practices for disease resistance phenotypes and increased the efficacy of porcine models in biomedical research; however limited functional annotation of the porcine immunome has hindered progress on both fronts. Among epigenetic mechanisms that regulate gene expression, DNA methylation is the most ubiquitous modification made to the DNA molecule and influences transcription factor binding as well as gene and phenotype expression. Human and mouse DNA methylation studies have improved mapping of regulatory elements in these species, but comparable studies in the pig have been limited in scope.

Results

We performed whole-genome bisulfite sequencing to assess DNA methylation patterns in nine pig immune cell populations: CD21⁺ and CD21⁻ B cells, four T cell fractions (CD4⁺, CD8⁺, CD8⁺CD4⁺, and SWC6γδ⁺), natural killer and myeloid cells, and neutrophils. We identified 54,391 cell differentially methylated regions (cDMRs), and clustering by cDMR methylation rate grouped samples by cell lineage. 32,737 cDMRs were classified as cell lowly methylated regions (cLMRs) in at least one cell type, and cLMRs were broadly enriched in genes and regions of intermediate CpG density. We observed strong correlations between differential methylation and expression across immune cell populations, with cell-specific low methylation disproportionately impacting genes exhibiting enriched gene expression in the same cell type. Motif analysis of cLMRs revealed cell type-specific enrichment of transcription factor binding motifs, indicating that cell-specific methylation patterns may influence accessibility by trans-acting factors. Lastly, cDMRs were enriched for immune capacity GWAS SNPs, and many such overlaps occurred within genes known to influence immune cell development and function (CD8B, NDRG1).

Conclusion

Our DNA methylation data improve functional annotation of the porcine genome through characterization of epigenomic regulatory patterns that contribute to immune cell identity and function, and increase the potential for identifying mechanistic links between genotype and phenotype.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08773-5.

A Co-Association Network Analysis Reveals Putative Regulators for Health-Related Traits in Pigs

In recent years, the increase in awareness of antimicrobial resistance together with the societal demand of healthier meat products have driven attention to health-related traits in livestock production. Previous studies have reported medium to high heritabilities for these traits and described genomic regions associated with them. Despite its genetic component, health- and immunity-related traits are complex and its study by association analysis with genomic markers may be missing some information. To analyse multiple phenotypes and gene-by-gene interactions, systems biology approaches, such as the association weight matrix (AWM), allows combining genome wide association study results with network inference algorithms. The present study aimed to identify gene networks, key regulators and candidate genes associated to immunocompetence in pigs by integrating multiple health-related traits, enriched for innate immune phenotypes, using the AWM approach. The co-association network analysis unveiled a network comprised of 3,636 nodes (genes) and 451,407 edges (interactions), including a total of 246 regulators. From these, five genes (ARNT2, BRMS1L, MED12L, SUPT3H and TRIM25) were selected as key regulators as they were associated with the maximum number of genes with the minimum overlapping (1,827 genes in total). The five regulators were involved in pathways related to immunity such as lymphocyte differentiation and activation, platelet activation and degranulation, megakaryocyte differentiation, FcγR-mediated phagocytosis and response to nitric oxide, among others, but also in immunometabolism. Furthermore, we identified genes co-associated with the key regulators previously reported as candidate genes (e.g., ANGPT1, CD4, CD36, DOCK1, PDE4B, PRKCE, PTPRC and SH2B3) for immunity traits in humans and pigs, but also new candidate ones (e.g., ACSL3, CXADR, HBB, MMP12, PTPN6, WLS) that were not previously described. The co-association analysis revealed new regulators associated with health-related traits in pigs. This approach also identified gene-by-gene interactions and candidate genes involved in pathways related to cell fate and metabolic and immune functions. Our results shed new light in the regulatory mechanisms involved in pig immunity and reinforce the use of the pig as biomedical model.

Genome-wide associations for immune traits in two maternal pig lines

Background

In recent years, animal welfare and health has become more and more important in pig breeding. So far, numerous parameters have been considered as important biomarkers, especially in the immune reaction and inflammation. Previous studies have shown moderate to high heritabilities in most of these traits. However, the genetic background of health and robustness of pigs needs to be extensively clarified. The objective of this study was to identify genomic regions with a biological relevance for the immunocompetence of piglets. Genome-wide Association Studies (GWAS) in 535 Landrace (LR) and 461 Large White (LW) piglets were performed, investigating 20 immune relevant traits. Besides the health indicators of the complete and differential blood count, eight different cytokines and haptoglobin were recorded in all piglets and their biological dams to capture mediating processes and acute phase reactions. Additionally, all animals were genotyped using the Illumina PorcineSNP60v2 BeadChip.

Results

In summary, GWAS detected 25 genome-wide and 452 chromosome-wide significant SNPs associated with 17 immune relevant traits in the two maternal pig lines LR and LW. Only small differences were observed considering the maternal immune records as covariate within the statistical model. Furthermore, the study identified across- and within-breed differences as well as relevant candidate genes. In LR more significant associations and related candidate genes were detected, compared with LW. The results detected in LR and LW are partly in accordance with previously identified quantitative trait loci (QTL) regions. In addition, promising novel genomic regions were identified which might be of interest for further detailed analysis. Especially putative pleiotropic regions on SSC5, SSC12, SSC15, SSC16 and SSC17 are of major interest with regard to the interacting structure of the immune system. The comparison with already identified QTL gives indications on interactions with traits affecting piglet survival and also production traits.

Conclusion

In conclusion, results suggest a polygenic and breed-specific background of immune relevant traits. The current study provides knowledge about regions with biological relevance for health and immune traits. Identified markers and putative pleiotropic regions provide first indications in the context of balancing a breeding-based modification of the porcine immune system.

Supplementary Information

The online version contains supplementary material available at (10.1186/s12864-021-07997-1).

Genome-wide association studies for production, respiratory disease, and immune-related traits in Landrace pigs

Identification of a quantitative trait locus (QTL) related to a chronic respiratory disease such as Mycoplasmal pneumonia of swine (MPS) and immune-related traits is important for the genetic improvement of disease resistance in pigs. The objective of this study was to detect a novel QTL for a total of 22 production, respiratory disease, and immune-related traits in Landrace pigs. A total of 874 Landrace purebred pigs, which were selected based on MPS resistance, were genotyped using the Illumina PorcineSNP60 BeadChip. We performed single nucleotide polymorphism (SNP)-based and haplotype-based genome-wide association studies (GWAS) to detect a novel QTL and to evaluate the possibility of a pleiotropic QTL for these traits. SNP-based GWAS detected a total of six significant regions in backfat thickness, ratio of granular leucocytes to lymphatic cells, plasma concentration of cortisol at different ages, and complement alternative pathway activity in serum. The significant region detected by haplotype-based GWAS was overlapped across the region detected by SNP-based GWAS. Most of these detected QTL regions were novel regions with some candidate genes located in them. With regard to a pleiotropic QTL among traits, only three of these detected QTL regions overlapped among traits, and many detected regions independently affected the traits.

Genetic parameters and associated genomic regions for global immunocompetence and other health-related traits in pigs

The inclusion of health-related traits, or functionally associated genetic markers, in pig breeding programs could contribute to produce more robust and disease resistant animals. The aim of the present work was to study the genetic determinism and genomic regions associated to global immunocompetence and health in a Duroc pig population. For this purpose, a set of 30 health-related traits covering immune (mainly innate), haematological, and stress parameters were measured in 432 healthy Duroc piglets aged 8 weeks. Moderate to high heritabilities were obtained for most traits and significant genetic correlations among them were observed. A genome wide association study pointed out 31 significantly associated SNPs at whole-genome level, located in six chromosomal regions on pig chromosomes SSC4, SSC6, SSC17 and SSCX, for IgG, γδ T-cells, C-reactive protein, lymphocytes phagocytic capacity, total number of lymphocytes, mean corpuscular volume and mean corpuscular haemoglobin. A total of 16 promising functionally-related candidate genes, including CRP, NFATC2, PRDX1, SLA, ST3GAL1, and VPS4A, have been proposed to explain the variation of immune and haematological traits. Our results enhance the knowledge of the genetic control of traits related with immunity and support the possibility of applying effective selection programs to improve immunocompetence in pigs.

Assessment of Autozygosity Derived From Runs of Homozygosity in Jinhua Pigs Disclosed by Sequencing Data

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.