Mapping quantitative trait loci for T lymphocyte subpopulations in peripheral blood in swine

Genome wide association study of thyroid hormone levels following challenge with porcine reproductive and respiratory syndrome virus

Introduction: Porcine reproductive and respiratory syndrome virus (PRRSV) causes respiratory disease in piglets and reproductive disease in sows. Piglet and fetal serum thyroid hormone (i.e., T3 and T4) levels decrease rapidly in response to Porcine reproductive and respiratory syndrome virus infection. However, the genetic control of T3 and T4 levels during infection is not completely understood. Our objective was to estimate genetic parameters and identify quantitative trait loci (QTL) for absolute T3 and/or T4 levels of piglets and fetuses challenged with Porcine reproductive and respiratory syndrome virus. Methods: Sera from 5-week-old pigs (N = 1792) at 11 days post inoculation (DPI) with Porcine reproductive and respiratory syndrome virus were assayed for T3 levels (piglet_T3). Sera from fetuses (N = 1,267) at 12 or 21 days post maternal inoculation (DPMI) with Porcine reproductive and respiratory syndrome virus of sows (N = 145) in late gestation were assayed for T3 (fetal_T3) and T4 (fetal_T4) levels. Animals were genotyped using 60 K Illumina or 650 K Affymetrix single nucleotide polymorphism (SNP) panels. Heritabilities, phenotypic correlations, and genetic correlations were estimated using ASREML; genome wide association studies were performed for each trait separately using Julia for Whole-genome Analysis Software (JWAS). Results: All three traits were low to moderately heritable (10%-16%). Phenotypic and genetic correlations of piglet_T3 levels with weight gain (0-42 DPI) were 0.26 ± 0.03 and 0.67 ± 0.14, respectively. Nine significant quantitative trait loci were identified for piglet_T3, on Sus scrofa chromosomes (SSC) 3, 4, 5, 6, 7, 14, 15, and 17, and collectively explaining 30% of the genetic variation (GV), with the largest quantitative trait loci identified on SSC5, explaining 15% of the genetic variation. Three significant quantitative trait loci were identified for fetal_T3 on SSC1 and SSC4, which collectively explained 10% of the genetic variation. Five significant quantitative trait loci were identified for fetal_T4 on SSC1, 6, 10, 13, and 15, which collectively explained 14% of the genetic variation. Several putative immune-related candidate genes were identified, including CD247, IRF8, and MAPK8. Discussion: Thyroid hormone levels following Porcine reproductive and respiratory syndrome virus infection were heritable and had positive genetic correlations with growth rate. Multiple quantitative trait loci with moderate effects were identified for T3 and T4 levels during challenge with Porcine reproductive and respiratory syndrome virus and candidate genes were identified, including several immune-related genes. These results advance our understanding of growth effects of both piglet and fetal response to Porcine reproductive and respiratory syndrome virus infection, revealing factors associated with genomic control of host resilience.

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.

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.

Genome-wide association studies for 30 haematological and blood clinical-biochemical traits in Large White pigs reveal genomic regions affecting intermediate phenotypes

Haematological and clinical-biochemical parameters are considered indicators of the physiological/health status of animals and might serve as intermediate phenotypes to link physiological aspects to production and disease resistance traits. The dissection of the genetic variability affecting these phenotypes might be useful to describe the resilience of the animals and to support the usefulness of the pig as animal model. Here, we analysed 15 haematological and 15 clinical-biochemical traits in 843 Italian Large White pigs, via three genome-wide association scan approaches (single-trait, multi-trait and Bayesian). We identified 52 quantitative trait loci (QTLs) associated with 29 out of 30 analysed blood parameters, with the most significant QTL identified on porcine chromosome 14 for basophil count. Some QTL regions harbour genes that may be the obvious candidates: QTLs for cholesterol parameters identified genes (ADCY8, APOB, ATG5, CDKAL1, PCSK5, PRL and SOX6) that are directly involved in cholesterol metabolism; other QTLs highlighted genes encoding the enzymes being measured [ALT (known also as GPT) and AST (known also as GOT)]. Moreover, the multivariate approach strengthened the association results for several candidate genes. The obtained results can contribute to define new measurable phenotypes that could be applied in breeding programs as proxies for more complex traits.

Genomic variant in porcine TNFRSF1A gene and its effects on TNF signaling pathway in vitro

Our initial genome-wide association study (GWAS) revealed the presence of single nucleotide polymorphisms (SNPs) related to immune traits in the tumor necrosis factor receptor superfamily member 1A (TNFRSF1A) gene, suggesting the association of this gene with immune function in pigs. To better understand the immune functions of the TNFRSF1A gene, SNPs within the TNFRSF1A gene were detected by sequencing. One SNP (c.1394C > T) in exon 6 of TNFRSF1A was identified, and association analysis in two pig populations was subsequently performed. The results showed that this SNP was significantly associated with CD4^-CD8^-CD3^-, CD4⁺CD8^-CD3⁺, and CD4⁺/CD8⁺ (P = 0.0038, P = 0.0007, and P = 0.0076, respectively). Based on quantitative real-time polymerase chain reaction (RT-qPCR), the TNFRSF1A mRNA was shown to be widely expressed in six different tissues. Finally, functional verification of the TNFRSF1A gene was performed in vitro to better understand its role. RNAi was used to generate a porcine PK15 cell line with a silenced TNFRSF1A gene, and a vector was also constructed to assess overexpression of TNFRSF1A. RT-qPCR was then used to detect changes in the expression levels of five critical genes. Our results indicated that TNFRSF1A activated the TNF signaling pathway and inhibited the NFκB signaling pathway in vitro. These findings provide evidence for an immune-related regulatory function for porcine TNFRSF1A.

Identification of porcine CTLA4 gene polymorphism and their association with piglet diarrhea and performance traits

The objective of this study was to evaluate the association between the cytotoxic T-lymphocyte-associated antigen 4 (CTLA4) gene and piglet diarrhea. In this study, the mRNA expression of the CTLA4 gene increased significantly in IPEC-J2 cells after Escherichia coli K88 infection. Single nucleotide polymorphisms (SNPs) located in the 5' flanking region (SNPs g.107281989C>T) and 3'-untranslated region (3'-UTR; SNPs g.107288753C>A) were identified, and they were in linkage disequilibrium in both Min pigs and the Landrace population. Association analysis showed that Landrace piglets with a TT or AA genotype had a lower diarrhea index, and AA animals had higher average daily gain when compared to CC pigs, respectively (p < 0.05). However, the relationship between SNPs and diarrhea and performance traits in the Min population was not significant. Haplotype analysis indicated that the TC haplotype had the lowest diarrhea index. The 5' flanking deletion assay suggested that SNP g.107281989C>T was a molecular marker instead of the functional marker. This research demonstrated that genetic variances in the CTLA4 gene had significant effects on Landrace piglet diarrhea resistance.

A genome-wide association study of fetal response to type 2 porcine reproductive and respiratory syndrome virus challenge

Control of porcine reproductive and respiratory syndrome (PRRS) is economically important for the swine industry worldwide. As current PRRS vaccines do not completely protect against heterologous challenge, alternative means of control, including enhanced genetic resilience, are needed. For reproductive PRRS, the genetic basis of fetal response to PRRS virus (PRRSV) infection is poorly understood. Genome-wide association studies (GWAS) were done here using data from 928 fetuses from pregnant gilts experimentally challenged with type 2 PRRSV. Fetuses were assessed for viral load in thymus (VLT), viral load in endometrium (VLE), fetal death (FD) and fetal viability (FV), and genotyped at a medium density. Collectively, 21 candidate genomic regions were found associated with these traits, seven of which overlap with previously reported QTLs for pig health and reproduction. A comparison with ongoing and related transcriptomic analyses of fetal response to PRRSV infection found differentially expressed genes within 18 candidate regions. Some of these genes have immune system functions, and have been reported to contribute to host response to PRRSV infection. The results provide new evidence about the genetic basis of fetal response to PRRSV challenge, and may ultimately lead to alternative control strategies to reduce the impact of reproductive PRRS.

Single-nucleotide polymorphisms in CD8A and their associations with T lymphocyte subpopulations in pig

Findings from previous studies suggested that the cluster of the differentiation 8 alpha (CD8A) gene plays a prominent role in human T lymphocyte subpopulations. However, the evidence from pig population is still rare. To determine whether the important role of the CD8A gene is conserved in pig, a candidate gene analysis was performed herein through genotype-phenotype associations. Five single-nucleotide polymorphisms (SNPs) locating in the regulatory region of porcine CD8A gene were detected and tested for association analysis with seven T lymphocyte subpopulations (proportion of CD4(-)CD8(-), CD4(+)CD8(+), CD4(+)CD8(-), CD4(-)CD8(+), CD4(+), CD8(+), and the ratio of CD4(+) to CD8(+) T cells in peripheral blood) in 382 Large White piglets. After Bonferroni correction for multiple testing, four SNPs were significantly associated with some or all of the seven T lymphocyte subpopulations. Analyses of pairwise D' measures of linkage disequilibrium between all SNPs were also explored. Two haplotype blocks was inferred and the association study on haplotype level revealed similar effects on T lymphocyte subpopulations. In addition, the tissue-specific RNA expression pattern and electrophoretic mobility shift assay offered further explanation of the link between the CD8A gene with porcine T lymphocyte subpopulations. The findings presented here provide strong evidence for associations of CD8A variants with T lymphocyte subpopulations and may be applied in porcine breeding programs.

A Whole Genome Association Study to Detect Single Nucleotide Polymorphisms for Blood Components (Immunity) in a Cross between Korean Native Pig and Yorkshire

The purpose of this study was to detect significant SNPs for blood components that were related to immunity using high single nucleotide polymorphism (SNP) density panels in a Korean native pig (KNP)×Yorkshire (YK) cross population. A reciprocal design of KNP×YK produced 249 F₂ individuals that were genotyped for a total of 46,865 available SNPs in the Illumina porcine 60K beadchip. To perform whole genome association analysis (WGA), phenotypes were regressed on each SNP under a simple linear regression model after adjustment for sex and slaughter age. To set up a significance threshold, 0.1% point-wise p value from F distribution was used for each SNP test. Among the significant SNPs for a trait, the best set of SNP markers were determined using a stepwise regression procedure with the rates of inclusion and exclusion of each SNP out of the model at 0.001 level. A total of 54 SNPs were detected; 10, 6, 4, 4, 5, 4, 5, 10, and 6 SNPs for neutrophil, lymphocyte, monocyte, eosinophil, basophil, atypical lymph, immunoglobulin, insulin, and insulin-like growth factor-I, respectively. Each set of significant SNPs per trait explained 24 to 42% of phenotypic variance. Several pleiotropic SNPs were detected on SSCs 4, 13, 14 and 15.

Association of the porcine cluster of differentiation 4 gene with T lymphocyte subpopulations and its expression in immune tissues

Cluster of differentiation 4 (CD4) is mainly expressed on CD4⁺ T cells, which plays an important role in immune response. The aim of this study was to detect the association between polymorphisms of the CD4 gene and T lymphocyte subpopulations in pigs, and to investigate the effects of genetic variation on the CD4 gene expression level in immune tissues. Five missense mutations in the CD4 gene were identified using DNA pooling sequencing assays, and two main haplotypes (CCTCC and AGCTG) in strong linkage disequilibrium (with frequencies of 50.26% and 46.34%, respectively) were detected in the population of Large White pigs. Our results indicated that the five SNPs and the two haplotypes were significantly associated with the proportions of CD4⁻CD8⁻, CD4⁺CD8⁺, CD4⁺CD8⁻, CD4⁺ and CD4⁺/CD8⁺ in peripheral blood (p<0.05). Gene expression analysis showed the mRNA level of the CD4 gene in thymus was significantly higher than that in lymph node and spleen (p<0.05). However, no significant difference was observed between animals with CCTCC/CCTCC genotype and animals with AGCTG/AGCTG genotype in the three immune tissues (p>0.05). These results indicate that the CD4 gene may influence T lymphocyte subpopulations and can be considered as a candidate gene affecting immunity in pigs.

Genome-wide association study for cytokines and immunoglobulin G in swine

Increased disease resistance through improved immune capacity would be beneficial for the welfare and productivity of farm animals. To identify genomic regions responsible for immune capacity traits in swine, a genome-wide association study was conducted. In total, 675 pigs were included. At 21 days of age, all piglets were vaccinated with modified live classical swine fever vaccine. Blood samples were sampled when the piglets were 20 and 35 days of age, respectively. Four traits, including Interferon-gamma (IFN-γ) and Interleukin 10 (IL-10) levels, the ratio of IFN-γ to IL-10 and Immunoglobulin G (IgG) blocking percentage to CSFV in serum were measured. All the samples were genotyped for 62,163 single nucleotide polymorphisms (SNP) using the Illumina porcineSNP60k BeadChip. After quality control, 46,079 SNPs were selected for association tests based on a single-locus regression model. To tackle the issue of multiple testing, 10,000 permutations were performed to determine the chromosome-wise and genome-wise significance level. In total, 32 SNPs with chromosome-wise significance level (including 4 SNPs with genome-wise significance level) were identified. These SNPs account for 3.23% to 13.81% of the total phenotypic variance individually. For the four traits, the numbers of significant SNPs range from 5 to 15, which jointly account for 37.52%, 82.94%, 26.74% and 24.16% of the total phenotypic variance of IFN-γ, IL-10, IFN-γ/IL-10, and IgG, respectively. Several significant SNPs are located within the QTL regions reported in previous studies. Furthermore, several significant SNPs fall into the regions which harbour a number of known immunity-related genes. Results herein lay a preliminary foundation for further identifying the causal mutations affecting swine immune capacity in follow-up studies.

Genome-wide association study for T lymphocyte subpopulations in swine

Background

Lymphocytes act as a major component of the adaptive immune system, taking very crucial responsibility for immunity. Differences in proportions of T-cell subpopulations in peripheral blood among individuals under same conditions provide evidence of genetic control on these traits, but little is known about the genetic mechanism of them, especially in swine. Identification of the genetic control on these variants may help the genetic improvement of immune capacity through selection.

Results

To identify genomic regions responsible for these immune traits in swine, a genome-wide association study was conducted. A total of 675 pigs of three breeds were involved in the study. At 21 days of age, all individuals were vaccinated with modified live classical swine fever vaccine. Blood samples were collected when the piglets were 20 and 35 days of age, respectively. Seven traits, including the proportions of CD4+, CD8+, CD4+CD8+, CD4+CD8−, CD4−CD8+, CD4−CD8− and the ratio of CD4+ to CD8+ T cells were measured at the two ages. All the samples were genotyped for 62,163 single nucleotide polymorphisms (SNP) using the Illumina porcineSNP60k BeadChip. 40833 SNPs were selected after quality control for association tests between SNPs and each immune trait considered based on a single-locus regression model. To tackle the issue of multiple testing in GWAS, 10,000 permutations were performed to determine the chromosome-wise and genome-wise significance levels of association tests. In total, 61 SNPs with chromosome-wise significance level and 3 SNPs with genome-wise significance level were identified. 27 significant SNPs were located within the immune-related QTL regions reported in previous studies. Furthermore, several significant SNPs fell into the regions harboring known immunity-related genes, 14 of them fell into the regions which harbor some known T cell-related genes.

Conclusions

Our study demonstrated that genome-wide association studies would be a feasible way for revealing the potential genetics variants affecting T-cell subpopulations. Results herein lay a preliminary foundation for further identifying the causal mutations underlying swine immune capacity in follow-up studies.