Skip to main content

ORIGINAL RESEARCH article

Front. Immunol., 26 November 2021
Sec. Comparative Immunology
This article is part of the Research Topic Porcine Innate Immunity View all 6 articles

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

  • Animal Breeding and Genetics Programme, Institut de Recerca i Tecnologia Agroalimentàries (IRTA), Torre Marimon, Caldes de Montbui, Spain

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.

Introduction

Health-related traits are becoming more and more relevant due to the great challenges that the pig sector has faced in recent years. The current emergence of antibiotic resistance, with pig production as one of the largest consumers of antimicrobials, and the societal demand for healthier livestock products and more sustainable production systems, make necessary to include new traits such as those related to health in the selection programs (1). Thus, health-related traits that have traditionally played a minor role in breeding programs are currently considered relevant traits to take into account.

Several studies have determined an important genetic contribution to the phenotypic variability of health-related traits, describing medium to high heritabilities for some of these traits (27) and identifying genomic regions with candidate genes associated to immunocompetence (7, 8). However, although genome wide association studies (GWAS) have identified single nucleotide polymorphisms (SNPs), or genes, associated with health-related traits, these SNPs generally had a little effect on the phenotypic variation of the traits under study due to their complex nature. Besides, some interesting SNPs could be overlooked due to stringent significance thresholds after multiple test correction (9). Furthermore, these studies do not consider gene-by-gene interactions or the contribution of genetic correlations between related traits.

During the last years, the use of systems biology approaches has allowed the integration of multiple sources of biological information, providing a more complete picture of the functional mechanisms involved in the determination of complex traits, and allowing the identification of candidate genes and genetic variants implicated in their phenotypic diversity (10). One of these approaches, the association weight matrix (AWM), simultaneously considers multiple phenotypes and gene-by-gene interactions by combining GWAS results with network inference algorithms (11, 12). To date, multiple studies in livestock species have implemented this approach to identify gene networks and key regulatory elements implicated in the genetic architecture of meat quality (13, 14), growth (15), gut microbial communities (16), reproductive (11, 17), feed efficiency (18, 19), and milk-related (2023) traits, among others. However, to the best of our knowledge, no studies have applied this methodology for health-related (mainly immunological and haematological) traits in livestock or other animal species.

The present study aimed to identify gene networks, key regulators and candidate genes associated to immuno-competence in pigs by integrating 30 health-related traits, enriched for innate immune phenotypes, using the AWM co-association analysis.

Material and Methods

Animal Material and Phenotypic Parameters

The animal material used in this study was previously reported in Ballester et al. (7). In brief, we used a total of 432 weaned piglets (217 males and 215 females) from a commercial Duroc pig line. The pigs were distributed in six batches and belonged to 134 litters obtained from 132 sows and 22 boars. All animals were raised in the same farm and fed ad libitum with a commercial cereal-based diet. At the time of sample collection (60 ± 8 days of age), all animals were apparently healthy, without any sign of infection.

Details of sampling and laboratory determination of immunological, haematological and stress parameters are described in Ballester et al. (7). In brief, total concentrations of immunoglobulins IgA, IgG and IgM in plasma, and IgA in saliva (IgAsal), were measured by ELISA with commercial kits (Bethyl laboratories Inc., Bionova, Spain), following manufacturer’s instructions. Concentration of acute phase proteins: C-reactive protein (CRP) and haptoglobin (HP) were measured in serum by ELISA kit (Abcam Plc., Spain) and colorimetric assay (Tridelta Development Limited, Ireland), respectively, following manufacturer’s instructions. Total concentration of nitric oxide (NO) in serum was measured by colorimetric assay (Termo Fisher Scientifc, Spain). To determine the percentage of gamma-delta T cells (γδ T cells), peripheral blood mononuclear cells (PBMCs) were separated from heparinized peripheral blood by density-gradient centrifugation with Histopaque-1077 (Sigma, Spain) and 106 PBMCs were stained using the monoclonal antibody APC Rat Anti-Pig γδ T Lymphocytes (MAC320 clone, BD Pharmigen, Spain) and the APC Rat IgG2a κ isotype control (R35-95 clone, BD Pharmigen, Spain). The cells were analysed by flow cytometry using the MACSQuant Analyzer 10 Flow cytometer (Miltenyi Biotec GmbH, Bergisch Gladbach, Germany) and the MACSQuantify sofware v2.6 (Miltenyi Biotec GmbH, Bergisch Gladbach, Germany). Phagocytosis assay was carried out in heparinized whole blood samples incubated with fluorescein (FITC)-labelled opsonized E. coli bacteria by using the Phagotest kit (BD Pharmigen, Spain) according to manufacturer’s instructions. The following percentages of phagocytosis traits in blood: percentage of total phagocytic cells (PHAGO_%); percentage of phagocytic cells among granulocytes (GRANU_PHAGO_%), monocytes (MON_PHAGO_%) and lymphocytes (LYM_PHAGO_%); mean fluorescence in fluorescein isothiocyanate (FITC) among the total phagocytic cells (PHAGO_FITC); and mean fluorescence in FITC among the granulocytes (GRANU_ PHAGO_FITC), monocytes (MON_PHAGO_FITC) and lymphocytes (LYM_PHAGO_FITC) that phagocyte were analysed by flow cytometry (MACSQuant Analyzer 10 Flow cytometer; Miltenyi Biotec GmbH, Bergisch Gladbach, Germany) and the MACSQuantify sofware v2.6 (Miltenyi Biotec GmbH, Bergisch Gladbach, Germany). Haematological traits: haematocrit (HCT), haemoglobin (HB), mean corpuscular volume (MCV), mean corpuscular haemoglobin (MCH), mean corpuscular haemoglobin concentration (MCHC), total number of leukocytes (LEU), eosinophils (EO), lymphocytes (LYM), monocytes (MON), neutrophils (NEU), erythrocytes (ERY) and platelets (PLA) were determined from blood samples by haemogram. Finally, the stress indicators analysed were the neutrophil to lymphocyte ratio (NLR) and the cortisol (CORT) levels in hair, which were measured by ELISA kit (Cusabio Technology LLC., Bionova, Spain). Table S1 summarizes the descriptive statistics of the analysed traits.

SNP Genotyping and Genome Wide Association Analysis

The 432 animals were genotyped with the GGP Porcine HD Array (Illumina, San Diego, CA) using the Infinium HD Assay Ultra protocol (Illumina). The quality control of the 68,516 SNPs was performed using Plink software (24), removing SNPs with a minor allele frequency lower than 5%, SNPs with more than 10% missing genotype data, and SNPs that did not map to the porcine reference genome (Sscrofa11.1 assembly). Thereafter, a subset of 42,641 SNPs was retained for subsequent analysis.

GWAS for the 30 health-related traits was performed using the following mixed linear model with the Genome-wide Complex Trait Analysis (GCTA) software tool (25):

yijk=sexj+bk+ui+slial+eijk

where yijk corresponds to the phenotypic trait (either log-transformed or raw data) of the ith individual of sex j in the kth batch; sexj corresponds to the jth sex effect (2 levels); bk corresponds to the kth batch effect (6 levels) for most traits but for phagocytosis related traits, for which the data of laboratory analysis (12 levels, two by batch) was considered instead; ui is the infinitesimal genetic effect of individual i, with u∼N(0,Gσu2), where G is the genomic relationship matrix calculated using the filtered autosomal SNPs based on the methodology of (25) and σu2 is the additive genetic variance; sli is the genotype (coded as 0,1,2) for the lth SNP, and al is the allele substitution effect of the SNP on the trait under study; and eijk is the residual term.

The GCTA software tool (25) was also used to estimate the heritability value for each one of the 30 phenotypes by performing a restricted maximum likelihood (REML) analysis with sex and batch as fixed effects.

Association Weight Matrix

The results from GWAS analysis were subjected to an AWM approach (11, 12). For the initial step, only those SNPs located in the coding region or within 5 kb of an annotated gene based on the Sscrofa11.1 reference genome assembly were kept. Then, γδ T cells was chosen as key phenotype and those γδ T cells-associated SNPs (P<0.05) were selected. In addition, we also retained those SNPs that were associated (P<0.05) with three or more of the 30 health-related phenotypes and referred to as pleiotropic SNPs. For assessing the gene-gene interactions, the standardized SNP effects across phenotypes were used to infer the gene co-association network using the Partial Correlation and Information Theory (PCIT) approaches (26). In the network, every node represents a gene (or a SNP), whereas every edge connecting two nodes represents a significant gene–gene interaction (based on SNP–SNP co-association). In order to quantify the number of connections as well as to prioritize potential regulators of the gene network we applied an information lossless approach (27) that explored the connectivity of all regulators (transcription factors, miRNAs and lnRNAs) and its target genes. The list of the regulators that was cross-matched with the list of 3,636 genes corresponds to the census of transcription factors reported by Vaquerizas et al. (28). The Cytoscape software (29) was used to visualize the gene network and calculate node centrality values.

Functional and Comparative Analyses

When only the Ensembl ID was available for some of the associated novel genes, they were searched against the database for human orthologues (GRCh38.p13) to extract their putative gene name using the BioMart web-based tool (30).

ClueGO plugin (31) was used to identify the over-represented gene ontology terms, KEGG pathways, and immune functions. The cut-off for considering a significant over-representation (corrected p-value ≤ 0.05) was established by Benjamini and Hochberg multiple-test correction (32). Further categorization of candidate genes was performed using information from the Mouse Genome Database (www.informatics.jax.org) (33) and GeneCards (www.genecards.org) (34).

A comparative analysis between our results and previous published data was performed by retrieving all pig and human quantitative trait loci (QTLs) and association data on Sscrofa11.1 and on GRCh38.p13 (dbSNP Build 153), for health traits from the pigQTL database (35) and from the NHGRI-EBI GWAS Catalog (v1.0.2) (36), respectively.

Results

Gene Co-Association Network Description

From the 42,641 SNPs used in the GWAS analysis, 2,230 SNPs were directly associated (P<0.05) with the key phenotype (γδ T cells), while 8,421 SNPs were associated (P<0.05) with three or more phenotypes. Thus, a total of 9,830 SNPs were obtained after merging both lists of associated SNPs. Two SNPs located in Sus scrofa chromosome (SSC) 13 (rs323669145 and rs81448146) were associated with the maximum number of phenotypes, 13 out of the 30 traits (CRP, HP, NO, CORT, MCHC, LYM, NLR, PHAGO_%, PHAGO_FITC, GRANU_PHAGO_%, GRANU_ PHAGO_FITC, MON_PHAGO_FITC, LYM_PHAGO_%). Those two SNPs were located in intronic regions of the fibroblast growth factor 12 (FGF12) gene. Among other SNPs associated with a high number of phenotypes (n≥10) we also identified rs81284215 on SSC9. Remarkably, this SNP was located in the haemoglobin subunit beta (HBB) gene and was associated with haematological-related phenotypes such as HCT, HB and ERY, as well as with γδ T cells and six out of the eight phagocytosis phenotypes, excluding those related to lymphocytes.

After the annotation step, only 3,544 out of the 9,830 SNPs were retained as they were located in the coding region or 5 kb upstream or downstream of a gene. However, as some SNPs were annotated to more than one gene, a total of 3,650 unique genes were identified. To explore the genetic determinism of the 3,544 SNPs selected from the AWM approach over the 30 phenotypes, their heritability was calculated using these SNPs and compared to the heritabilities obtained using the whole dataset comprised by 42,641 SNPs (Table 1). Heritability calculated with the SNPs subset ranged between 0.238 and 0.791. The highest heritabilities (h2>0.7) were observed for the haematological traits (MCH, MCV, HB, HCT), while CRP, MON, and IgAsal showed low heritabilities (h2<0.3). However, the other immunoglobulin-related traits (IgG, IgA and IgM) had heritabilities above 0.5. In addition, most of the traits (24 out of 30) showed moderate to high heritabilities (h2>0.4). When compared with the heritabilities obtained with the whole dataset, the heritabilities of the SNPs subset showed higher values in all comparisons except one (LYM_PHAGO_FITC), whose values were almost similar (0.374 vs. 0.379). The greatest improvement in the heritability value was observed for HB and HCT, while others, such as the values for CRP, IgG, and the aforementioned LYM_PHAGO_FITC, did not change much. Of note, the heritability of the key phenotype (γδ T cells) was also greatly improved from moderate in the whole dataset (0.439) to high in the subset (0.688). Furthermore, to validate the heritabilities obtained by the reduced dataset of 3,544 SNPs, the heritability values of the 30 phenotypes were calculated by randomly selecting the same number of SNPs from the whole dataset in 100 different iterations. Thus, for each phenotype, the heritability values of the reduced dataset were always higher than the heritability values randomly obtained with each iteration with the exception, again, of the LYM_PHAGO_FITC phenotype.

TABLE 1
www.frontiersin.org

Table 1 Heritability values (h2), and their standard errors (SE), for the immunological, haematological and stress analysed traits obtained using the 3,544 SNPs selected from the AWM approach and the whole dataset (42,641 SNPs).

Genetic relationships among phenotypes based on the normalized additive values of the 3,544 SNPs are shown in Figure 1. The phagocytosis capacity phenotypes (measured through FITC) and the phenotypes for the proportion of phagocytic cells were highly correlated (Figure 1A and Table S2), and clustered together, jointly with NEU and NLR (Figure 1B). Haematological erythrocyte-related measurements had also a high positive correlation, similar to the ones reported for the white blood cell types (Table S2). These phenotypes were separated in a second cluster, with HB, HCT and ERY being sub-clustered together and in close proximity to MCH and MCV, and all leukocyte-related counts but NEU (i.e., LEU, EO, MON and LYM) were grouped together in a differentiated branch (Figure 1B). Associated to them, the last sub-cluster gathered most of the innate immunity traits such as positive correlated acute-phase proteins (CRP and HP) or positive correlated immunoglobulins concentrations (IgA, IgG and IgM) in plasma, but also other phenotypes as CORT, PLA or MCHC. The latter showed a high correlation with IgA measured in saliva (Figure 1A and Table S2). The percentage of γδ T cells, the key phenotype used in the AWM approach, was also enclosed in this last sub-cluster, but far from the rest of phenotypes, since it was not correlated with any phenotype (Figure 1A).

FIGURE 1
www.frontiersin.org

Figure 1 Genetic relationships among phenotypes based on the normalized additive values of the 3,544 SNPs selected for the AWM approach. (A) Heatmap of the correlations estimated by pairwise combinations and (B) Hierarchical cluster analysis among immunity, haematological and stress related traits in pigs.

The gene co-association network was comprised of 3,636 nodes and 451,407 interactions (Figure S1). The topology of the network reflected an average degree of 248.299, a clustering coefficient of 0.347 and an average path length of 2.219, which manifest a great level of connection between the nodes of the network. Of note, the three genes with the highest degree (i.e., the highest number of connections) were REPS2, ACSL3 and WLS, with 893, 858 and 837, respectively. ACSL3 and WLS knockout mice showed abnormal hematopoietic and immune responses, including decreased γδ T cell number in the latter (37).

Key Transcription Factors Regulating Health-Related Traits

After PCIT/network analysis, we identified a total of 246 potential regulators, including three microRNAs, within the 3,636 co-associated genes for the 30 health-related traits. The information lossless approach (12) allowed identifying trios of regulators with maximum connectivity (i.e., interactions) within the network; the top 10 trios are shown in Table 2. Among them, those composed by BRMS1L, MED12L and SUPT3H, and by ARNT2, SUPT3H and TRIM25, were among the top trios that spanned most of the network topology with highest connectivity (1,482 and 1,472 unique interactions, respectively) and minimum redundancy (Table 2). BRMS1L and TRIM25 were the regulators with the highest number of co-associated genes (n=603), followed by SUPT3H (n=580), ARNT2 (n=552), and MED12L (n=466). All these protein-coding genes are transcription factors with repressor (BRMS1L), co-activator (MED12L), or activator functions. For each significant association between a SNP of the five key transcription factors and a trait, the allele substitution effect of the SNP on the respective phenotype is summarised in Table S3.

TABLE 2
www.frontiersin.org

Table 2 Top 10 regulator trios based by the number of unique interactions.

Functional Classification and Identification of Candidate Genes Associated to Immunocompetence

A functional classification considering all the genes (n=1,828) co-associated with the five top regulators was performed to identify biological processes, metabolic pathways and candidate genes directly related to immunocompetence and study their role in immune system processes. The list of immune biological processes and pathways identified is shown in Table S4. Overall, genes for immune cells activation, differentiation and proliferation, regulation of innate and adaptive immune responses, phagocytosis, immunoglobulin and cytokine production and regulation of inflammatory response processes were identified. It is worth to highlight the JAK/STAT pathway, the ERK1/2 and MAPK cascades, and the NOTCH and Phospholipase D signaling pathways among those processes. Furthermore, we also identified other metabolic pathways directly related to immunometabolism such as the AMPK, mTOR, PI3K-Akt and PPAR signaling pathways, which control oxidative phosphorylation, carbohydrate, lipid and peptide metabolic processes, among others. This functional analysis revealed that a total of 589 candidate genes belonged to immune-related functions (Table S4).

In addition, we studied in detail for each top transcription factor those immune processes in which a high number of co-associated genes were involved (Table S5). Figure 2 shows a simplified co-association network representing those immunity-associated genes interacting with the top 5 regulators. Table 3 summarizes the main immune and metabolic processes in which the top transcription factors were associated. It is worthy to note that all the regulators were co-associated with genes related to platelet functions and phospholipase D signaling. Furthermore, among the top five regulators SUPT3H was the transcription factor co-associated with the highest number of pathways and processes related to immunometabolism and haemopoiesis and immune development, and BRMS1L was the only transcription factor co-associated with genes related to the NOTCH signaling pathway. Finally, we identified a specific gene, protein tyrosine phosphatase receptor type C (PTPRC), involved in γδ T cell differentiation and activation that was co-associated with ARNT2, BRMS1L and MED12L.

FIGURE 2
www.frontiersin.org

Figure 2 Simplification of the full gene co-association network formed by 171 nodes and 832 interactions representing the immunity-associated genes (orange) that were correlated with the five top regulator genes (green; ARNT2, BRMS1L, MED12L, SUPT3H, and TRIM25).

TABLE 3
www.frontiersin.org

Table 3 Summary of the main immune and metabolic processes associated with the top transcription factors.

Comparison With Literature Reported Candidate Genes and QTLs for Health-Related Traits

When we studied in detail the top five regulators and their co-associated genes related to immune functions, we identified interesting candidate genes enclosed in genomic regions previously associated with immune phenotypes in studies with humans and pigs (Table S6). For instance, PRKCE, a candidate gene for haematological parameters in humans (3841) and different leukocyte subpopulations in pigs (42), was associated in our study with phagocytosis-related traits. Other genes also associated with phagocytosis-related traits in our population were CD36 and SH2B3. In humans, CD36 has been associated with haematological traits (38, 39), while SH2B3 has been associated with a plethora of immune phenotypes, mainly with leukocyte levels and their types (3840). MYH9 was another candidate gene associated with haematological parameters in humans (3840) together with ANGPT1 (38, 39, 43), where the latter was also associated with the levels of exhaled nitric oxide and IgG glycosylation in serum (44, 45). In our analysis, MYH9 was associated with LYM and phagocytosis-related traits while ANGPT1 was also associated with NO and some haematological traits. The previously mentioned PTPRC (also known as CD45) was associated in humans with leukocyte, lymphocyte and monocyte counts (38, 39, 43) and memory T cell differentiation (46). In our population this gene was associated with LYM, among other traits. Other candidate genes commonly associated with immune cell types in humans were: ATXN1 and CLEC16A (3840, 43), CD4 (46), FTO and KLRK1 (47); while in pigs were AKT2 (48), CLEC16A (42), JMJD1C and RBFOX1 (8), PRKG2 and PAQR3 (49), and SLC30A10 (50). CD4 and KLRK1 are mainly expressed in T cells and Natural Killer cells and they were associated with γδ T cells in our population. Another interesting gene described in association with CRP in humans (51) and with platelet levels in pigs (50) was PDE4B. This gene was associated with haematological and phagocytosis-related traits in our study.

Furthermore, while none of the top five regulators were located in QTL regions associated with health-related traits in our Duroc pig population (7), they were located in QTLs for C3c concentration, HCT, HB, Interleukin (IL)-2 and IL-10 levels, Toll-like receptor (TLR)-2 and TLR-9 levels, and basophil number, among others (Table S6), in other pig populations (35).

Discussion

In the present study, a systems genetic approach using the AWM methodology was used to identify potential regulators, gene-gene interactions and pathways for immunocompetence in pigs considering 30 health-related phenotypes enriched for innate immune phenotypes. We selected γδ T cells as the key phenotype, since these cells play important roles at the interface between innate and adaptive immunity and had the potential to impact other phenotypes analysed in our population. In fact, in a previous study in the same population (7), positive and negative genetic correlations were identified between this phenotype and acute phase proteins and platelets (positive correlations), and IgAsal, PHAGO-FITC and cortisol levels (negative correlations). Similar functions such as response to pathogens and tumours, wound healing, tissue homeostasis, stressed cell clearance, and antigen presentation, among others, have been described for γδ T cells in several species [reviewed in Holderness et al. (52)]. However, compared to humans and mice, this population of cells is strongly enriched in pig blood (53), which also suggest specific immune responses and functions.

After GWAS analysis and the SNP annotation step, only 3,544 SNPs were retained for further analysis. When compared with the heritability values obtained with the whole dataset of SNPs, this subset of 3,544 SNPs helped to explain a greater proportion of the additive variance of the analysed traits. As expected, the heritability of our key phenotype (γδ T cells) was greatly improved from moderate to high after using the reduced dataset, despite other phenotypes such as HB and HCT had a greater improvement and a higher heritability value. This may also be due to the fact that they were the two phenotypes with the highest number of associated SNPs after the key phenotype. Among the associated SNPs with HB and HCT, we identified a pleiotropic SNP (rs81284215) that was associated with more than 10 phenotypes and was located in the HBB gene. It is well established the relationship between HB levels and growth on the early-stages of life (54, 55). Similarly, a positive association between HB and HCT and average daily weight gain has also been described in pigs (56). Therefore, this SNP could be a useful marker to determine the haematological status of our animals, although further analysis will be required to determine the relationship of this polymorphism with weight gain in piglets.

Interestingly, the AWM additive effects correlations of the 3,544 SNPs primarily clustered health-related traits according to their nature and function and were informative enough to capture the genetic associations among the health-related traits in a similar way than the published genetic correlations calculated for the same population (7). Of note, there was a lack of genetic correlation among our key phenotype and the rest of the health-related phenotypes with 225 SNPs exclusively associated with γδ T cells. However, we also identified 534 SNPs associated with γδ T cells and one or more health-related phenotypes, highlighting the genetic interplay between several immune pathways and functions. A possible explanation for these results is the multitasking role of γδ T cells in the immune system with a conventional adaptive role, but with increasing relevance in their influences on innate immunity (57, 58). In fact, it has been proposed a direct but also and indirect role of γδ T cells in promoting monocyte and macrophage differentiation and macrophage recruitment during inflammatory response, as well as in macrophage elimination during the resolution of inflammation (57). Additionally, a specialized role of epithelial γδ T cells in wound repair by promoting platelet and neutrophil accumulation in the limbal vessels has been proposed (59). Overall, these functions are the reflection of the results found during this work that will be discussed below.

In this study, AWM co-association analysis allowed the identification of five key regulators (ARNT2, BRMS1L, MED12L, SUPT3H and TRIM25) for immunocompetence in pigs. It is relevant to note that for most of these regulators (BRMS1L, MED12L, SUPT3H) there is still little information in the literature about their relationship with immunity-related functions, despite all regulators were located within QTL regions associated to immunity phenotypes (pigQTLdb). Among the identified top regulators, BRMS1L and TRIM25 gathered the highest number of associations with other genes. While BRMS1L has been described as a metastatic suppressor gene (60), TRIM25 has been involved in numerous cellular processes including regulation of innate immune response against viral infection (61, 62). Specifically, TRIM25 is involved in the RIG-I/interferon pathway in response to viral RNAs (63). Among the list of genes regulated by TRIM25, we identified two genes (MMP12 and PTPN6) related to regulation of type I interferon-mediated signaling pathway. Another key regulator that also plays a role in the innate immune system was ARNT2. This transcription factor forms heterodimers with the aryl hydrocarbon receptor (AHR), that is expressed in all γδ T cell lymphocyte subsets, being essential for some aspects of their functions [reviewed in Stange and Veldhoen (64)]. This complex regulates the transcription of multiple immune-related genes and the activity and differentiation of phagocytic cells and lymphocytes (65). Accordingly, this transcription factor was associated with leukocyte, lymphocyte and monocyte counts in our population and was co-associated with the PTPRC gene that was related to γδ T cell differentiation. Furthermore, a polymorphism in ARNT2 has been associated with an impaired fungicidal activity by depleting the phagocytic activity of macrophages (66). Notably, in our gene-network we identified ARNT2-regulated genes involved in FcγR-mediated phagocytosis (ACTR3C, AKT2, DOCK1, GAB2, PRKCE, PTPRC) and response to nitric oxide (AIFM1, CCR7, MTR, TXN2) processes. Therefore, these genes are interesting candidate genes to contribute to the phenotypic differences observed for these health-related traits. Consistently with our results, some of these genes (AKT2, PRKCE, PTPRC) were previously associated with health-related traits in humans and/or pigs, and seven of the ten genes have been at least associated with a phagocytosis-related trait in our population. MED12L is a paralog of the MED12 gene, which is required to activate innate immunity genes in Drosophila (67). A recent GWAS study identified a variant within the MED12L gene associated with viral infection response (Epstein-Barr virus (EBV) nuclear antigen (EBNA)) in humans (68). In the network, the MED12L gene was also interacting with genes associated to γδ T cell differentiation (PTPRC) and activation (CXADR and PTPRC). Finally, a genetic study found an association between a polymorphism in the SUPT3H gene with the ADAMTS13 activity, a protein with antithrombotic properties (69).

Except for some genes directly related to the key phenotype (γδ T cell), the five regulators and their 1,828 co-associated genes were involved in a wide variety of biological functions, some of them directly related to immunity but also to immunometabolism. Among them, there were functions related to lymphocytes (differentiation and activation), platelets (activation, degranulation and megakaryocyte differentiation), phagocytosis (FcγR-mediated phagocytosis), NO (response to nitric oxide), and to a lesser extent other functions associated with other immune phenotypes (such as those related to haematological related traits) were identified. Remarkably, for some of these phenotypes (platelets count, NO in serum, and some phagocytosis traits) we did not identify any significantly associated genomic region after performing GWAS (7). Here, it is relevant to highlight that for some of these traits, we identified more than one of the top transcription factors co-associated with genes related to the same functions. Therefore, a plausible explanation for our results would be that due to the complex nature of immunity traits with the involvement of multiple genes and regulators, GWAS analyses are not able to capture all the significant regions associated with their phenotypic variation, since many of the biologically relevant polymorphisms are lost after multiple test corrections. For instance, it is worthy to note the overrepresentation of platelet associated functions, which were commonly found for all the key identified regulators. Apart from being important for coagulation and fibrinolysis processes, platelets are also activators of the immune system and mediators of inflammatory response [reviewed in Trzeciak-Ryczek et al. (70)], therefore participating in multiple functions. In fact, in the hierarchical analysis, platelets phenotype grouped within a heterologous group composed by stress indicators, acute phase proteins and immunoglobulins among others, which may be indicative of its multifunctional nature. Another overrepresented pathway regulated by all the key transcription factors was phospholipase D (PLD) signaling pathway. PLD is an enzyme that catalyses the hydrolysis of phosphatidylcholine, to produce the signal molecule phosphatidic acid (PA). PLD participates in multiple functions of the immune system mediating phagocytosis by macrophages, activation of NADPH oxidation in neutrophils, mast cell degranulation and T cell activation [reviewed in Zhu et al. (71)]. Finally, it is worth mentioning some specific immune functions for the poorly characterized BRMS1L and SUPT3H regulators; for instance, the regulation of the NOTCH signaling pathway identified for BRMS1L. NOTCH signaling plays a pivotal role in T cell fate decision, influencing the differentiation between αβ T cells versus γδ T cells (72, 73). Furthermore, SUPT3H seems to have an active role in immunometabolism by controlling signaling pathways such as AMPK, mTOR and PI3K-Akt. In recent years, immunometabolism has emerged as an important mechanism affecting innate and adaptive immune system. In fact, AMPK and PI3K/Akt/mTOR regulate metabolic pathways important for cell fate and function in the immune system (7476). Overall, these results underline the relevance of the five key transcription factors identified in our network regulating important functions of the immune system in pigs. Furthermore, and taking into account the plethora of different functions that these pathways/molecules perform, we could hypothesize, that depending on the function, either one or the interaction of several key transcription factors will be required.

Among the genes co-associated with the top five regulators, we also identified candidate genes that have been previously described as candidate genes for health-related traits in humans and pigs, some of them associated with similar immunological and haematological traits (ANGPT1, PTPRC, CD4, KLRK1, and FTO). In addition, there were also candidate genes with functions directly related to the associated traits. Indeed, PRKCE, SH2B3 and CD36 were associated in our study with phagocytosis-related traits. PRKCE plays a role in FcγR-mediated phagocytosis (77), while CD36 is a phagocytic receptor that mediates the uptake of different substrates (78, 79). SH2B3 is involved in macrophage function (80), whereas MYH9, which encodes a conventional non-muscle myosin and is required for normal T lymphocyte migration (81), was associated with the quantity of lymphocytes. ANGPT1, a candidate gene in humans, belongs to the angiopoietin family regulating vascular integrity and angiogenesis through endothelium-derived NO release (82), and accordingly, it was also found in association with NO and some haematological traits in our pig population. Also, it is relevant to highlight that the candidate gene PTPRC was associated with the quantity of lymphocytes in blood. This member of the protein tyrosine phosphatase (PTP) family is expressed in all nucleated haematopoietic cells and plays an essential role in T and B cell antigen receptor signaling, and it has been associated to many autoimmune disorders [reviewed in Ulloa-Aguirre et al. (83)]. Finally, PDE4B, which was associated with haematological and phagocytosis-related traits in our study, is essential for the LPS-activated immune response of phagocytic cells (82, 84).

In conclusion, the AWM gene co-association analysis applied in the present study represented an important step to better characterize the genetic architecture of complex traits for immunity in pigs. The network-based approach allowed the identification of five key transcription factors (ARNT2, BRMS1L, MED12L, SUPT3H and TRIM25) regulating pathways directly involved in cell fate and in metabolic and immune functions. Furthermore, we have identified gene-gene interactions and candidate genes, some of which have been previously associated to immunity traits in humans and pigs. This study has been the first one to use AWM for analysing health-related traits in pigs and supports the effectiveness of using network-based approaches to identify potential regulators and candidate genes, representing a major step in understanding the regulatory mechanisms involved in porcine immunity, and reinforcing the use of the pig as biomedical model.

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding authors.

Ethics Statement

All experimental procedures with pigs were performed according to the Spanish Policy for Animal Protection RD 53/2013, which meets the European Union Directive 2010/63/EU about the protection of animals used in experimentation. The experimental protocol was approved by the Ethical Committee of the Institut de Recerca i Tecnologia Agroalimentàries (IRTA).

Author Contributions

RQ and MB designed the study. MB supervised the generation of the animal material used in this work. YR-C, OG-R, MP, RQ and MB performed the sampling. OG-R and MB carried out the laboratory analyses. DC-P, YR-C, RQ and MB analysed the data and interpreted the results. DC-P and MB wrote the manuscript. All authors contributed to the article and approved the submitted version.

Funding

The study was funded by grants AGL2016-75432-R and PID2020-112677RB-C21 awarded by the Spanish Ministry of Economy and Competitiveness (MINECO). YR-C was funded by Marie Skłodowska-Curie grant (P-Sphere) agreement No 6655919 530 (EU). MB was financially supported by a Ramon y Cajal contract (RYC-2013-12573) from the Spanish Ministry of Economy and Competitiveness. The authors belonged to a Consolidated Research Group AGAUR, ref. 2017SGR-1719.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

We gratefully acknowledge the technical staff from IRTA and Selección Batallé S.A., and specially to Josep Reixach, for their collaboration in the experimental protocols at farm and slaughterhouse.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2021.784978/full#supplementary-material

Glossary

www.frontiersin.org

References

1. Merks JWM, Mathur PK, Knol EF. New Phenotypes for New Breeding Goals in Pigs. Animal (2012) 6:535–43. doi: 10.1017/S1751731111002266

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Edfors-Lilja I, Wattrang E, Magnusson U, Fossum C. Genetic Variation in Parameters Reflecting Immune Competence of Swine. Vet Immunol Immunopathol (1994) 40:1–16. doi: 10.1016/0165-2427(94)90011-6

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Henryon M, Heegaard PMH, Nielsen J, Berg P, Juul-Madsen HR. Immunological Traits Have the Potential to Improve Selection of Pigs for Resistance to Clinical and Subclinical Disease. Anim Sci (2006) 82:597–606. doi: 10.1079/ASC200671

CrossRef Full Text | Google Scholar

4. Clapperton M, Glass EJ, Bishop SC. Pig Peripheral Blood Mononuclear Leucocyte Subsets Are Heritable and Genetically Correlated With Performance. Animal (2008) 2:1575–84. doi: 10.1017/S1751731108002929

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Clapperton M, Diack AB, Matika O, Glass EJ, Gladney CD, Mellencamp MA, et al. Traits Associated With Innate and Adaptive Immunity in Pigs: Heritability and Associations With Performance Under Different Health Status Conditions. Genet Sel Evol (2009) 41:54. doi: 10.1186/1297-9686-41-54

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Flori L, Gao Y, Laloë D, Lemonnier G, Leplat J-J, Teillaud A, et al. Immunity Traits in Pigs: Substantial Genetic Variation and Limited Covariation. PloS One (2011) 6:e22717. doi: 10.1371/journal.pone.0022717

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Ballester M, Ramayo-Caldas Y, González-Rodríguez O, Pascual M, Reixach J, Díaz M, et al. Genetic Parameters and Associated Genomic Regions for Global Immunocompetence and Other Health-Related Traits in Pigs. Sci Rep (2020) 10:18462. doi: 10.1038/s41598-020-75417-7

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Bovo S, Mazzoni G, Bertolini F, Schiavo G, Galimberti G, Gallo M, et al. Genome-Wide Association Studies for 30 Haematological and Blood Clinical-Biochemical Traits in Large White Pigs Reveal Genomic Regions Affecting Intermediate Phenotypes. Sci Rep (2019) 9:7003. doi: 10.1038/s41598-019-43297-1

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Rice TK, Schork NJ, Rao DC. Methods for Handling Multiple Testing. Adv Genet (2008) 60:293–308. doi: 10.1016/S0065-2660(07)00412-9

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Civelek M, Lusis AJ. Systems Genetics Approaches to Understand Complex Traits. Nat Rev Genet (2014) 15:34–48. doi: 10.1038/nrg3575

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Fortes MRS, Reverter A, Zhang Y, Collis E, Nagaraj SH, Jonsson NN, et al. Association Weight Matrix for the Genetic Dissection of Puberty in Beef Cattle. Proc Natl Acad Sci USA (2010) 107:13642–7. doi: 10.1073/pnas.1002044107

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Reverter A, Fortes MRS. Association Weight Matrix: A Network-Based Approach Towards Functional Genome-Wide Association Studies. Methods Mol Biol (2013) 1019:437–47. doi: 10.1007/978-1-62703-447-0_20

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Ramayo-Caldas Y, Ballester M, Fortes MRS, Esteve-Codina A, Castelló A, Noguera JL, et al. From SNP Co-Association to RNA Co-Expression: Novel Insights Into Gene Networks for Intramuscular Fatty Acid Composition in Porcine. BMC Genomics (2014) 15:232. doi: 10.1186/1471-2164-15-232

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Ramayo-Caldas Y, Renand G, Ballester M, Saintilan R, Rocha D. Multi-Breed and Multi-Trait Co-Association Analysis of Meat Tenderness and Other Meat Quality Traits in Three French Beef Cattle Breeds. Genet Sel Evol (2016) 48:37. doi: 10.1186/s12711-016-0216-y

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Puig-Oliveras A, Ballester M, Corominas J, Revilla M, Estellé J, Fernández AI, et al. Folch JM. A Co-Association Network Analysis of the Genetic Determination of Pig Conformation, Growth and Fatness. PloS One (2014) 9:e114862. doi: 10.1371/journal.pone.0114862

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Reverter A, Ballester M, Alexandre PA, Mármol-Sánchez E, Dalmau A, Quintanilla R, et al. A Gene Co-Association Network Regulating Gut Microbial Communities in a Duroc Pig Population. Microbiome (2021) 9:52. doi: 10.1186/s40168-020-00994-8

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Kang J-H, Lee E-A, Hong K-C, Kim J-M. Regulatory Gene Network From a Genome-Wide Association Study for Sow Lifetime Productivity Traits. Anim Genet (2018) 49:254–8. doi: 10.1111/age.12640

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Taussat S, Boussaha M, Ramayo-Caldas Y, Martin P, Venot E, Cantalapiedra-Hijar G, et al. Gene Networks for Three Feed Efficiency Criteria Reveal Shared and Specific Biological Processes. Genet Sel Evol (2020) 52:67. doi: 10.1186/s12711-020-00585-z

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Ramayo-Caldas Y, Mármol-Sánchez E, Ballester M, Sánchez JP, González-Prendes R, Amills M, et al. Integrating Genome-Wide Co-Association and Gene Expression to Identify Putative Regulators and Predictors of Feed Efficiency in Pigs. Genet Sel Evol (2019) 51:48. doi: 10.1186/s12711-019-0490-6

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Pegolo S, Dadousis C, Mach N, Ramayo-Caldas Y, Mele M, Conte G, et al. SNP Co-Association and Network Analyses Identify E2F3, KDM5A and BACH2 as Key Regulators of the Bovine Milk Fatty Acid Profile. Sci Rep (2017) 7:17317. doi: 10.1038/s41598-017-17434-7

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Sanchez M-P, Ramayo-Caldas Y, Wolf V, Laithier C, El Jabri M, Michenet A, et al. Sequence-Based GWAS, Network and Pathway Analyses Reveal Genes Co-Associated With Milk Cheese-Making Properties and Milk Composition in Montbéliarde Cows. Genet Sel Evol (2019) 51:34. doi: 10.1186/s12711-019-0473-7

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Pegolo S, Mach N, Ramayo-Caldas Y, Schiavon S, Bittante G, Cecchinato A. Integration of GWAS, Pathway and Network Analyses Reveals Novel Mechanistic Insights Into the Synthesis of Milk Proteins in Dairy Cows. Sci Rep (2018) 8:566. doi: 10.1038/s41598-017-18916-4

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Marete A, Lund MS, Boichard D. Ramayo-Caldas Y. A System-Based Analysis of the Genetic Determinism of Udder Conformation and Health Phenotypes Across Three French Dairy Cattle Breeds. PloS One (2018) 13:e0199931. doi: 10.1371/journal.pone.0199931

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MAR, Bender D, et al. PLINK: A Tool Set for Whole-Genome Association and Population-Based Linkage Analyses. Am J Hum Genet (2007) 81:559–75. doi: 10.1086/519795

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Yang J, Lee SH, Goddard ME, Visscher PM. GCTA: A Tool for Genome-Wide Complex Trait Analysis. Am J Hum Genet (2011) 88:76–82. doi: 10.1016/j.ajhg.2010.11.011

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Reverter A, Chan EKF. Combining Partial Correlation and an Information Theory Approach to the Reversed Engineering of Gene Co-Expression Networks. Bioinformatics (2008) 24:2491–7. doi: 10.1093/bioinformatics/btn482

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Reverter A, Fortes MRS. Breeding and Genetics Symposium: Building Single Nucleotide Polymorphism-Derived Gene Regulatory Networks: Towards Functional Genomewide Association Studies. J Anim Sci (2013) 91:530–6. doi: 10.2527/jas.2012-5780

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Vaquerizas JM, Kummerfeld SK, Teichmann SA. Luscombe NM. A Census of Human Transcription Factors: Function, Expression and Evolution. Nat Rev Genet (2009) 10:252–63. doi: 10.1038/nrg2538

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. Cytoscape: A Software Environment for Integrated Models of Biomolecular Interaction Networks. Genome Res (2003) 13:2498–504. doi: 10.1101/gr.1239303

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Kinsella RJ, Kähäri A, Haider S, Zamora J, Proctor G, Spudich G, et al. Ensembl BioMarts: A Hub for Data Retrieval Across Taxonomic Space. Database (Oxford) (2011) 2011:bar030. doi: 10.1093/database/bar030

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Bindea G, Mlecnik B, Hackl H, Charoentong P, Tosolini M, Kirilovsky A, et al. ClueGO: A Cytoscape Plug-in to Decipher Functionally Grouped Gene Ontology and Pathway Annotation Networks. Bioinformatics (2009) 25:1091–3. doi: 10.1093/bioinformatics/btp101

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Benjamini Y, Hochberg Y. Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. J R Stat Soc Ser B (1995) 57:289–300. doi: 10.1111/j.2517-6161.1995.tb02031.x

CrossRef Full Text | Google Scholar

33. Bult CJ, Blake JA, Smith CL, Kadin JA, Richardson JE. Mouse Genome Database Group. Mouse Genome Database (MGD) 2019. Nucleic Acids Res (2019) 47:D801–6. doi: 10.1093/nar/gky1056

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Stelzer G, Rosen N, Plaschkes I, Zimmerman S, Twik M, Fishilevich S, et al. The GeneCards Suite: From Gene Data Mining to Disease Genome Sequence Analyses. Curr Protoc Bioinform (2016) 54:1.30.1–1.30.33. doi: 10.1002/cpbi.5

CrossRef Full Text | Google Scholar

35. Hu Z-L, Park CA, Reecy JM. Building a Livestock Genetic and Genomic Information Knowledgebase Through Integrative Developments of Animal QTLdb and CorrDB. Nucleic Acids Res (2019) 47:D701–10. doi: 10.1093/nar/gky1084

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Buniello A, MacArthur JAL, Cerezo M, Harris LW, Hayhurst J, Malangone C, et al. The NHGRI-EBI GWAS Catalog of Published Genome-Wide Association Studies, Targeted Arrays and Summary Statistics 2019. Nucleic Acids Res (2019) 47:D1005–12. doi: 10.1093/nar/gky1120

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Augustin I, Gross J, Baumann D, Korn C, Kerr G, Grigoryan T, et al. Loss of Epidermal Evi/Wls Results in a Phenotype Resembling Psoriasiform Dermatitis. J Exp Med (2013) 210:1761–77. doi: 10.1084/jem.20121871

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Chen M-H, Raffield LM, Mousas A, Sakaue S, Huffman JE, Moscati A, et al. Trans-Ethnic and Ancestry-Specific Blood-Cell Genetics in 746,667 Individuals From 5 Global Populations. Cell (2020) 182:1198–1213.e14. doi: 10.1016/j.cell.2020.06.045

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Vuckovic D, Bao EL, Akbari P, Lareau CA, Mousas A, Jiang T, et al. The Polygenic and Monogenic Basis of Blood Traits and Diseases. Cell (2020) 182:1214–1231.e11. doi: 10.1016/j.cell.2020.08.008

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Astle WJ, Elding H, Jiang T, Allen D, Ruklisa D, Mann AL, et al. The Allelic Landscape of Human Blood Cell Trait Variation and Links to Common Complex Disease. Cell (2016) 167:1415–1429.e19. doi: 10.1016/j.cell.2016.10.042

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Hodonsky CJ, Baldassari AR, Bien SA, Raffield LM, Highland HM, Sitlani CM, et al. Ancestry-Specific Associations Identified in Genome-Wide Combined-Phenotype Study of Red Blood Cell Traits Emphasize Benefits of Diversity in Genomics. BMC Genomics (2020) 21:228. doi: 10.1186/s12864-020-6626-9

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Zhang J, Chen JH, Liu XD, Wang HY, Liu XL, Li XY, et al. Genomewide Association Studies for Hematological Traits and T Lymphocyte Subpopulations in a Duroc × Erhualian F Resource Population. J Anim Sci (2016) 94:5028–41. doi: 10.2527/jas.2016-0924

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Kichaev G, Bhatia G, Loh P-R, Gazal S, Burch K, Freund MK, et al. Leveraging Polygenic Functional Enrichment to Improve GWAS Power. Am J Hum Genet (2019) 104:65–75. doi: 10.1016/j.ajhg.2018.11.008

PubMed Abstract | CrossRef Full Text | Google Scholar

44. van der Valk RJ, Duijts L, Timpson NJ, Salam MT, Standl M, Curtin JA, et al. Fraction of Exhaled Nitric Oxide Values in Childhood Are Associated With 17q11.2-Q12 and 17q12-Q21 Variants. J Allergy Clin Immunol (2014) 134:46–55. doi: 10.1016/j.jaci.2013.08.053

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Lauc G, Huffman JE, Pučić M, Zgaga L, Adamczyk B, Mužinić A, et al. Loci Associated With N-Glycosylation of Human Immunoglobulin G Show Pleiotropy With Autoimmune Diseases and Haematological Cancers. PloS Genet (2013) 9:e1003225. doi: 10.1371/journal.pgen.1003225

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Lagou V, Garcia-Perez JE, Smets I, Van Horebeek L, Vandebergh M, Chen L, et al. Genetic Architecture of Adaptive Immune System Identifies Key Immune Regulators. Cell Rep (2018) 25:798–810.e6. doi: 10.1016/j.celrep.2018.09.048

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Roederer M, Quaye L, Mangino M, Beddall MH, Mahnke Y, Chattopadhyay P, et al. The Genetic Architecture of the Human Immune System: A Bioresource for Autoimmunity and Disease Pathogenesis. Cell (2015) 161:387–403. doi: 10.1016/j.cell.2015.02.046

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Lu X, Liu J-F, Gong Y-F, Wang Z-P, Liu Y, Zhang Q. Mapping Quantitative Trait Loci for T Lymphocyte Subpopulations in Peripheral Blood in Swine. BMC Genet (2011) 12:79. doi: 10.1186/1471-2156-12-79

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Reiner G, Fischer R, Hepp S, Berge T, Köhler F, Willems H. Quantitative Trait Loci for White Blood Cell Numbers in Swine. Anim Genet (2008) 39:163–8. doi: 10.1111/j.1365-2052.2008.01700.x

PubMed Abstract | CrossRef Full Text | Google Scholar

50. Ponsuksili S, Reyer H, Trakooljul N, Murani E, Wimmers K. Single- and Bayesian Multi-Marker Genome-Wide Association for Haematological Parameters in Pigs. PloS One (2016) 11:e0159212. doi: 10.1371/journal.pone.0159212

PubMed Abstract | CrossRef Full Text | Google Scholar

51. Han X, Ong J-S, An J, Hewitt AW, Gharahkhani P, MacGregor S. Using Mendelian Randomization to Evaluate the Causal Relationship Between Serum C-Reactive Protein Levels and Age-Related Macular Degeneration. Eur J Epidemiol (2020) 35:139–46. doi: 10.1007/s10654-019-00598-z

PubMed Abstract | CrossRef Full Text | Google Scholar

52. Holderness J, Hedges JF, Ramstead A, Jutila MA. Comparative Biology of γδ T Cell Function in Humans, Mice, and Domestic Animals. Annu Rev Anim Biosci (2013) 1:99–124. doi: 10.1146/annurev-animal-031412-103639

PubMed Abstract | CrossRef Full Text | Google Scholar

53. Mair KH, Sedlak C, Käser T, Pasternak A, Levast B, Gerner W, et al. The Porcine Innate Immune System: An Update. Dev Comp Immunol (2014) 45:321–43. doi: 10.1016/j.dci.2014.03.022

PubMed Abstract | CrossRef Full Text | Google Scholar

54. Gentry JL, Swinkels JW, Lindemann MD, Schrama JW. Effect of Hemoglobin and Immunization Status on Energy Metabolism of Weanling Pigs. J Anim Sci (1997) 75:1032–40. doi: 10.2527/1997.7541032x

PubMed Abstract | CrossRef Full Text | Google Scholar

55. Larson LM, Kubes JN, Ramírez-Luzuriaga MJ, Khishen S, H Shankar A, Prado EL. Effects of Increased Hemoglobin on Child Growth, Development, and Disease: A Systematic Review and Meta-Analysis. Ann N Y Acad Sci (2019) 1450:83–104. doi: 10.1111/nyas.14105

PubMed Abstract | CrossRef Full Text | Google Scholar

56. Bhattarai S, Nielsen JP. Association Between Hematological Status at Weaning and Weight Gain Post-Weaning in Piglets. Livest Sci (2015) 182:64–8. doi: 10.1016/j.livsci.2015.10.017

CrossRef Full Text | Google Scholar

57. Born WK, Reardon CL, O’Brien RL. The Function of γδ T Cells in Innate Immunity. Curr Opin Immunol (2006) 18:31–8. doi: 10.1016/j.coi.2005.11.007

PubMed Abstract | CrossRef Full Text | Google Scholar

58. Vantourout P, Hayday A. Six-Of-the-Best: Unique Contributions of γδ T Cells to Immunology. Nat Rev Immunol (2013) 13:88–100. doi: 10.1038/nri3384

PubMed Abstract | CrossRef Full Text | Google Scholar

59. Li Z, Burns AR, Rumbaut RE, Smith CW. γδ T Cells Are Necessary for Platelet and Neutrophil Accumulation in Limbal Vessels and Efficient Epithelial Repair After Corneal Abrasion. Am J Pathol (2007) 171:838–45. doi: 10.2353/ajpath.2007.070008

PubMed Abstract | CrossRef Full Text | Google Scholar

60. Gong C, Qu S, Lv X-B, Liu B, Tan W, Nie Y, et al. BRMS1L Suppresses Breast Cancer Metastasis by Inducing Epigenetic Silence of FZD10. Nat Commun (2014) 5:5406. doi: 10.1038/ncomms6406

PubMed Abstract | CrossRef Full Text | Google Scholar

61. Martín-Vicente M, Medrano LM, Resino S, García-Sastre A, Martínez I. TRIM25 in the Regulation of the Antiviral Innate Immunity. Front Immunol (2017) 8:1187. doi: 10.3389/fimmu.2017.01187

PubMed Abstract | CrossRef Full Text | Google Scholar

62. Heikel G, Choudhury NR, Michlewski G. The Role of Trim25 in Development, Disease and RNA Metabolism. Biochem Soc Trans (2016) 44:1045–50. doi: 10.1042/BST20160077

PubMed Abstract | CrossRef Full Text | Google Scholar

63. Gack MU, Shin YC, Joo C-H, Urano T, Liang C, Sun L, et al. TRIM25 RING-Finger E3 Ubiquitin Ligase is Essential for RIG-I-Mediated Antiviral Activity. Nature (2007) 446:916–20. doi: 10.1038/nature05732

PubMed Abstract | CrossRef Full Text | Google Scholar

64. Stange J, Veldhoen M. The Aryl Hydrocarbon Receptor in Innate T Cell Immunity. Semin Immunopathol (2013) 35:645–55. doi: 10.1007/s00281-013-0389-1

PubMed Abstract | CrossRef Full Text | Google Scholar

65. Stockinger B, Di Meglio P, Gialitakis M, Duarte JH. The Aryl Hydrocarbon Receptor: Multitasking in the Immune System. Annu Rev Immunol (2014) 32:403–32. doi: 10.1146/annurev-immunol-032713-120245

PubMed Abstract | CrossRef Full Text | Google Scholar

66. Lupiañez CB, Martínez-Bueno M, Sánchez-Maldonado JM, Badiola J, Cunha C, Springer J, et al. Polymorphisms Within the ARNT2 and CX3CR1 Genes Are Associated With the Risk of Developing Invasive Aspergillosis. Infect Immun (2020) 88:e00882–19. doi: 10.1128/IAI.00882-19

PubMed Abstract | CrossRef Full Text | Google Scholar

67. Kuuluvainen E, Hakala H, Havula E, Sahal Estimé M, Rämet M, Hietakangas V, et al. Cyclin-Dependent Kinase 8 Module Expression Profiling Reveals Requirement of Mediator Subunits 12 and 13 for Transcription of Serpent-Dependent Innate Immunity Genes in Drosophila. J Biol Chem (2014) 289:16252–61. doi: 10.1074/jbc.M113.541904

PubMed Abstract | CrossRef Full Text | Google Scholar

68. Kachuri L, Francis SS, Morrison ML, Wendt GA, Bossé Y, Cavazos TB, et al. The Landscape of Host Genetic Factors Involved in Immune Response to Common Viral Infections. Genome Med (2020) 12:93. doi: 10.1186/s13073-020-00790-x

PubMed Abstract | CrossRef Full Text | Google Scholar

69. de Vries PS, Boender J, Sonneveld MAH, Rivadeneira F, Ikram MA, Rottensteiner H, et al. Genetic Variants in the ADAMTS13 and SUPT3H Genes Are Associated With ADAMTS13 Activity. Blood (2015) 125:3949–55. doi: 10.1182/blood-2015-02-629865

PubMed Abstract | CrossRef Full Text | Google Scholar

70. Trzeciak-Ryczek A, Tokarz-Deptuła B, Deptuła W. Platelets – an Important Element of the Immune System. Pol J Vet Sci (2013) 16:407–13. doi: 10.2478/pjvs-2013-0058

PubMed Abstract | CrossRef Full Text | Google Scholar

71. Zhu M, Foreman DP, O’Brien SA, Jin Y, Zhang W. Phospholipase D in TCR-Mediated Signaling and T Cell Activation. J Immunol (2018) 200:2165–73. doi: 10.4049/jimmunol.1701291

PubMed Abstract | CrossRef Full Text | Google Scholar

72. Washburn T, Schweighoffer E, Gridley T, Chang D, Fowlkes BJ, Cado D, et al. Notch Activity Influences the αβ Versus γδ T Cell Lineage Decision. Cell (1997) 88:833–43. doi: 10.1016/s0092-8674(00)81929-7

PubMed Abstract | CrossRef Full Text | Google Scholar

73. Gogoi D, Dar AA, Chiplunkar SV. Involvement of Notch in Activation and Effector Functions of γδ T Cells. J Immunol (2014) 192:2054–62. doi: 10.4049/jimmunol.1300369

PubMed Abstract | CrossRef Full Text | Google Scholar

74. Pearce EL, Pearce EJ. Metabolic Pathways in Immune Cell Activation and Quiescence. Immunity (2013) 38:633–43. doi: 10.1016/j.immuni.2013.04.005

PubMed Abstract | CrossRef Full Text | Google Scholar

75. Makowski L, Chaib M, Rathmell JC. Immunometabolism: From Basic Mechanisms to Translation. Immunol Rev (2020) 295:5–14. doi: 10.1111/imr.12858

PubMed Abstract | CrossRef Full Text | Google Scholar

76. Saxton RA, Sabatini DM. mTOR Signaling in Growth, Metabolism, and Disease. Cell (2017) 168:960–76. doi: 10.1016/j.cell.2017.02.004

PubMed Abstract | CrossRef Full Text | Google Scholar

77. Larsen EC, Ueyama T, Brannock PM, Shirai Y, Saito N, Larsson C, et al. A Role for PKC-ϵ in Fcγr-Mediated Phagocytosis by RAW 264.7 Cells. J Cell Biol (2002) 159:939–44. doi: 10.1083/jcb.200205140

PubMed Abstract | CrossRef Full Text | Google Scholar

78. Fadok VA, Warner ML, Bratton DL, Henson PM. CD36 is Required for Phagocytosis of Apoptotic Cells by Human Macrophages That Use Either a Phosphatidylserine Receptor or the Vitronectin Receptor (Alpha V Beta 3). J Immunol (1998) 161:6250–7.

PubMed Abstract | Google Scholar

79. Grajchen E, Wouters E, van de Haterd B, Haidar M, Hardonnière K, Dierckx T, et al. CD36-Mediated Uptake of Myelin Debris by Macrophages and Microglia Reduces Neuroinflammation. J Neuroinflamm (2020) 17:224. doi: 10.1186/s12974-020-01899-x

CrossRef Full Text | Google Scholar

80. Gueller S, Goodridge HS, Niebuhr B, Xing H, Koren-Michowitz M, Serve H, et al. Adaptor Protein Lnk Inhibits C-Fms-Mediated Macrophage Function. J Leukoc Biol (2010) 88:699–706. doi: 10.1189/jlb.0309185

PubMed Abstract | CrossRef Full Text | Google Scholar

81. Morin NA, Oakes PW, Hyun Y-M, Lee D, Chin YE, Chin EY, et al. Nonmuscle Myosin Heavy Chain IIA Mediates Integrin LFA-1 De-Adhesion During T Lymphocyte Migration. J Exp Med (2008) 205:195–205. doi: 10.1084/jem.20071543

PubMed Abstract | CrossRef Full Text | Google Scholar

82. Babaei S, Teichert-Kuliszewska K, Zhang Q, Jones N, Dumont DJ, Stewart DJ. Angiogenic Actions of Angiopoietin-1 Require Endothelium-Derived Nitric Oxide. Am J Pathol (2003) 162:1927–36. doi: 10.1016/S0002-9440(10)64326-X

PubMed Abstract | CrossRef Full Text | Google Scholar

83. Ulloa-Aguirre A, Zariñán T, Castillo-Badillo JA, Maravillas-Montero JL. Membrane Receptor Signaling in Health and Disease. In: Encyclopedia of Hormones. USA: Elsevier Science (2003). pp. 664–57. doi: 10.1016/B0-12-341103-3/00199-6

CrossRef Full Text | Google Scholar

84. Yang J-X, Hsiung T-C, Weng F-C, Ding S-L, Wu C-P, Conti M, et al. Synergistic Effect of Phosphodiesterase 4 Inhibitor and Serum on Migration of Endotoxin-Stimulated Macrophages. Innate Immun (2018) 24:501–12. doi: 10.1177/1753425918809155

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: immunocompetence, pig, gene networks, transcription factors, systems biology, candidate genes

Citation: Crespo-Piazuelo D, Ramayo-Caldas Y, González-Rodríguez O, Pascual M, Quintanilla R and Ballester M (2021) A Co-Association Network Analysis Reveals Putative Regulators for Health-Related Traits in Pigs. Front. Immunol. 12:784978. doi: 10.3389/fimmu.2021.784978

Received: 28 September 2021; Accepted: 08 November 2021;
Published: 26 November 2021.

Edited by:

Geert Wiegertjes, Wageningen University and Research, Netherlands

Reviewed by:

Yinming Liang, Xinxiang Medical University, China
Apratim Mitra, Eunice Kennedy Shriver National Institute of Child Health and Human Development (NICHD), United States

Copyright © 2021 Crespo-Piazuelo, Ramayo-Caldas, González-Rodríguez, Pascual, Quintanilla and Ballester. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Daniel Crespo-Piazuelo, ZGFuaWVsLmNyZXNwb0BpcnRhLmNhdA==; Maria Ballester, bWFyaWEuYmFsbGVzdGVyQGlydGEuY2F0

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.