Immunity traits in pigs: substantial genetic variation and limited covariation

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.

Multivariate genome-wide associations for immune traits in two maternal pig lines

BACKGROUND

Immune traits are considered to serve as potential biomarkers for pig's health. Medium to high heritabilities have been observed for some of the immune traits suggesting genetic variability of these phenotypes. Consideration of previously established genetic correlations between immune traits can be used to identify pleiotropic genetic markers. Therefore, genome-wide association study (GWAS) approaches are required to explore the joint genetic foundation for health biomarkers. Usually, GWAS explores phenotypes in a univariate (uv), trait-by-trait manner. Besides two uv GWAS methods, four multivariate (mv) GWAS approaches were applied on combinations out of 22 immune traits for Landrace (LR) and Large White (LW) pig lines.

RESULTS

In total 433 (LR: 351, LW: 82) associations were identified with the uv approach implemented in PLINK and a Bayesian linear regression uv approach (BIMBAM) software. Single Nucleotide Polymorphisms (SNPs) that were identified with both uv approaches (n = 32) were mostly associated with immune traits such as haptoglobin, red blood cell characteristics and cytokines, and were located in protein-coding genes. Mv GWAS approaches detected 647 associations for different mv immune trait combinations which were summarized to 133 Quantitative Trait Loci (QTL). SNPs for different trait combinations (n = 66) were detected with more than one mv method. Most of these SNPs are associated with red blood cell related immune trait combinations. Functional annotation of these QTL revealed 453 immune-relevant protein-coding genes. With uv methods shared markers were not observed between the breeds, whereas mv approaches were able to detect two conjoint SNPs for LR and LW. Due to unmapped positions for these markers, their functional annotation was not clarified.

CONCLUSIONS

This study evaluated the joint genetic background of immune traits in LR and LW piglets through the application of various uv and mv GWAS approaches. In comparison to uv methods, mv methodologies identified more significant associations, which might reflect the pleiotropic background of the immune system more accurately. In genetic research of complex traits, the SNP effects are generally small. Furthermore, one genetic variant can affect several correlated immune traits at the same time, termed pleiotropy. As mv GWAS methods consider strong dependencies among traits, the power to detect SNPs can be boosted. Both methods revealed immune-relevant potential candidate genes. Our results indicate that one single test is not able to detect all the different types of genetic effects in the most powerful manner and therefore, the methods should be applied complementary.

A promising resilience parameter for breeding: the use of weight and feed trajectories in growing pigs

BACKGROUND

Increasing resilience is a priority in modern pig breeding. Recent research shows that general resilience can be quantified via variability in longitudinal data. The collection of such longitudinal data on weight, feed intake and feeding behaviour in pigs has been facilitated by the development of technologies such as automated feeding stations. The goal of this study was to investigate resilience traits, which were estimated as deviations from longitudinal weight, feed intake and feeding behaviour data during the finishing phase. A dataset with 324,207 records between the age of 95 and 155 days on 5,939 Piétrain pigs with known pedigree and genomic information was used. We provided guidelines for a rigid quality control of longitudinal body weight data, as we found that outliers can significantly affect results. Gompertz growth curve analysis, linear modelling and trajectory analyses were used for quantifying resilience traits.

RESULTS

To our knowledge, this is the first study comparing resilience traits from longitudinal body weight, feed intake and feeding behaviour data in pigs. We demonstrated that the resilience traits are lowly to moderately heritable for deviations in body weight (h² = 2.9%-20.2%), in feed intake (9.4%-23.3%) and in feeding behaviour (16.2%-28.3%). Additionally, these traits have good predictive abilities in cross-validation analyses. Deviations in individual body weight and feed intake trajectories are highly correlated (r_g = 0.78) with low to moderate favourable genetic correlations with feed conversion ratio (r_g = 0.39-0.49). Lastly, we showed that some resilience traits, such as the natural logarithm of variances of observed versus predicted body weights (lnvar_weight), are more robust to lower observation frequencies and are repeatable over three different time periods of the finishing phase.

CONCLUSIONS

Our results will help future studies investigating resilience traits and resilience-related traits. Moreover, our study provides first results on standardization of quality control and efficient data sampling from automated feeding station data. Our findings will be valuable for breeding organizations as they offer evidence that pigs' general resilience can be selected on with good accuracy. Moreover, this methodology might be extended to other species to quantify resilience based on longitudinal data.

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.

Genetic parameters of immune traits for Landrace and Large White pig breeds

Improving the immunocompetence towards pathogens represents a desirable objective of breeding strategies to increase resilience. However, the immune system is complex and the genetic foundation of the underlying components is not yet clarified. In the present study, we focused on 22 blood parameters of 1,144 Landrace (LR) and Large White (LW) piglets at the age of 6-7 weeks. The immune profiles covered immune cells, red blood cell characteristics and cytokines. Genetic parameters based on pedigree information along with possible environmental effects were estimated. Litter effects play an important role in the expression of immune parameters of their young progenies. Hence, litter impacts on the piglet's immune profile including the immune parameters of the dam itself were investigated by different models. To incorporate the complexity of the immune network, the data were further investigated with a principal component analysis. Immune traits showed low to high breed-specific heritabilities (h² ). Strong positive r_g were estimated among red blood cell characteristics (0.77-0.99) and among cytokines (0.48-0.99). Neutrophils and lymphocytes illustrated a high negative r_g (-0.96 to -0.98). The litter impact on piglet's immunity was examined and strengthened already observed breed differences. In LR, h² (0.22-0.15) and litter effect (c² ) (0.52-0.44) for IFN-γ decreased after statistical consideration of maternal impact. In LW, a decrease in h² (0.32-0.18) for IFN-γ and an increase in c² (0.54-0.56) were observed. Here, sufficient correlations were detected within various immune traits and functional biological networks of principal components. Most immune traits are heritable and are promising to cover global breed-specific immunocompetence in pigs. The analysis of immune traits has to be extended in order to find an optimal range and to characterize relationships between immunity and performance to gain an improved immune system without accidental losses in productivity.

Application of Gene Editing Technology in Resistance Breeding of Livestock

As a new genetic engineering technology, gene editing can precisely modify the specific gene sequence of the organism’s genome. In the last 10 years, with the rapid development of gene editing technology, zinc-finger nucleases (ZFNs), transcription activator-like endonucleases (TALENs), and CRISPR/Cas9 systems have been applied to modify endogenous genes in organisms accurately. Now, gene editing technology has been used in mice, zebrafish, pigs, cattle, goats, sheep, rabbits, monkeys, and other species. Breeding for disease-resistance in agricultural animals tends to be a difficult task for traditional breeding, but gene editing technology has made this easier. In this work, we overview the development and application of gene editing technology in the resistance breeding of livestock. Also, we further discuss the prospects and outlooks of gene editing technology in disease-resistance breeding.

Plasma Cytokines and Birth Weight as Biomarkers of Vaccine-Induced Humoral Responses in Piglets

Failure to mount an effective immune response to vaccination leaves individuals at risk for infection and can compromise herd immunity. Vaccine unresponsiveness can range from poor responses “low responders” to a failure to seroconvert “non-responders.” Biomarkers of vaccine unresponsiveness, particularly those measured at the time of vaccination, could facilitate more strategic vaccination programs. We previously reported that pro-inflammatory cytokine signaling within peripheral blood mononuclear cells, elevated plasma interferon-gamma (IFNγ), and low birth weight correlated with vaccine-induced serum IgG titers in piglets that were below the threshold of detectable seroconversion (vaccine non-responders). These observations suggested that plasma IFNγ concentration and birth weight might serve as pre-vaccination biomarkers of vaccine unresponsiveness. To test this hypothesis, piglets (n = 67) from a different production facility were vaccinated with the same commercial Mycoplasma hyopneumoniae bacterin (RespiSure-One) to determine if there was a consistent and significant association between vaccine-induced serum IgG titers and either plasma cytokine concentrations or birth weight. All piglets seroconverted following vaccination with significantly less variability in vaccine-induced serum IgG titers than observed in the previous vaccine trial. Piglets exhibited highly variable birth weights and plasma cytokine concentrations prior to vaccination, but there were no significant associations (p > 0.05) between these variables and vaccine-induced serum IgG titers. There were significant (p < 0.001) differences in plasma IFNγ concentrations among individual litters (n = 6), and plasma IFNγ concentrations decreased in all pigs from birth to 63-days of age. One of the six litters (n = 11 piglets) exhibited significantly elevated plasma IFNγ concentrations during the first 3 weeks of life (p < 0.001) and at the time of vaccination (p < 0.01). This litter, however, had similar vaccine-induced serum IgG titers when compared to the other piglets in this study. Collectively the two studies indicate that while plasma cytokines and birth weight can be associated with vaccine non-responsiveness, their temporal and individual variation, as well as the complexity of the vaccine responsiveness phenotype, make them inconsistent biomarkers for predicting the less extreme phenotype of vaccine low responders.

Breeding for disease resilience: opportunities to manage polymicrobial challenge and improve commercial performance in the pig industry

Disease resilience, defined as an animal's ability to maintain productive performance in the face of infection, provides opportunities to manage the polymicrobial challenge common in pig production. Disease resilience can deliver a number of benefits, including more sustainable production as well as improved animal health and the potential for reduced antimicrobial use. However, little progress has been made to date in the application of disease resilience in breeding programs due to a number of factors, including (1) confusion around definitions of disease resilience and its component traits disease resistance and tolerance, and (2) the difficulty in characterizing such a complex trait consisting of multiple biological functions and dynamic elements of rates of response and recovery from infection. Accordingly, this review refines the definitions of disease resistance, tolerance, and resilience based on previous studies to help improve the understanding and application of these breeding goals and traits under different scenarios. We also describe and summarize results from a "natural disease challenge model" designed to provide inputs for selection of disease resilience. The next steps for managing polymicrobial challenges faced by the pig industry will include the development of large-scale multi-omics data, new phenotyping technologies, and mathematical and statistical methods adapted to these data. Genome editing to produce pigs resistant to major diseases may complement selection for disease resilience along with continued efforts in the more traditional areas of biosecurity, vaccination and treatment. Altogether genomic approaches provide exciting opportunities for the pig industry to overcome the challenges provided by hard-to-manage diseases as well as new environmental challenges associated with climate change.

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.

Leveraging host-genetics and gut microbiota to determine immunocompetence in pigs

BACKGROUND

The gut microbiota influences host performance playing a relevant role in homeostasis and function of the immune system. The aim of the present work was to identify microbial signatures linked to immunity traits and to characterize the contribution of host-genome and gut microbiota to the immunocompetence in healthy pigs.

RESULTS

To achieve this goal, we undertook a combination of network, mixed model and microbial-wide association studies (MWAS) for 21 immunity traits and the relative abundance of gut bacterial communities in 389 pigs genotyped for 70K SNPs. The heritability (h2; proportion of phenotypic variance explained by the host genetics) and microbiability (m2; proportion of variance explained by the microbial composition) showed similar values for most of the analyzed immunity traits, except for both IgM and IgG in plasma that was dominated by the host genetics, and the haptoglobin in serum which was the trait with larger m2 (0.275) compared to h2 (0.138). Results from the MWAS suggested a polymicrobial nature of the immunocompetence in pigs and revealed associations between pigs gut microbiota composition and 15 of the analyzed traits. The lymphocytes phagocytic capacity (quantified as mean fluorescence) and the total number of monocytes in blood were the traits associated with the largest number of taxa (6 taxa). Among the associations identified by MWAS, 30% were confirmed by an information theory network approach. The strongest confirmed associations were between Fibrobacter and phagocytic capacity of lymphocytes (r = 0.37), followed by correlations between Streptococcus and the percentage of phagocytic lymphocytes (r = -0.34) and between Megasphaera and serum concentration of haptoglobin (r = 0.26). In the interaction network, Streptococcus and percentage of phagocytic lymphocytes were the keystone bacterial and immune-trait, respectively.

CONCLUSIONS

Overall, our findings reveal an important connection between gut microbiota composition and immunity traits in pigs, and highlight the need to consider both sources of information, host genome and microbial levels, to accurately characterize immunocompetence in pigs.

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

Single Nucleotide Polymorphisms of Immunity-Related Genes and Their Effects on Immunophenotypes in Different Pig Breeds

Enhancing resistance and tolerance to pathogens remains an important selection objective in the production of livestock animals. Single nucleotide polymorphisms (SNPs) vary gene expression at the transcriptional level, influencing an individual's immune regulation and susceptibility to diseases. In this study, we investigated the distribution of SNP sites in immune-related genes and their correlations with cell surface markers of immune cells within purebred (Taiwan black, Duroc, Landrace and Yorkshire) and crossbred (Landrace-Yorkshire) pigs. Thirty-nine SNPs of immune-related genes, including 11 cytokines, 5 chemokines and 23 Toll-like receptors (TLRs) (interferon-α and γ (IFN-α, γ), tumor necrosis factor-α (TNF-α), granulocyte-macrophage colony-stimulating factor (GM-CSF), Monocyte chemoattractant protein-1 (MCP-1) and TLR3, TLR4, TLR7, TLR8, and TLR9) were selected, and the percentages of positive cells with five cell surface markers of CD4, CD8, CD80/86, MHCI, and MHCII were analyzed. There were 28 SNPs that were significantly different among breeds, particularly between Landrace and Taiwan black. For instance, the frequency of SNP1 IFN-α -235A/G in Taiwan black and Landrace was 11.11% and 96.15%, respectively. In addition, 18 SNPs significantly correlated with the expression of cell surface markers, including CD4, CD8, CD80/86, and MHCII. The percentage of CD4+ (39.27%) in SNP33 TLR-8 543C/C was significantly higher than those in A/C (24.34%), at p < 0.05. Together, our findings show that Taiwan black pigs had a unique genotype distribution, whereas Landrace and Yorkshire had a more similar genotype distribution. Thus, an understanding of the genetic uniqueness of each breed could help to identify functionally important SNPs in immunoregulation.

Proliferation of peripheral blood mononuclear cells from healthy piglets after mitogen stimulation as indicators of disease resilience

Disease resilience refers to the productivity of an animal under disease. Given the high biosecurity of pig nucleus herds, traits that can be measured on healthy pigs and that are genetically correlated with disease resilience, that is, genetic indicator traits, offer a strategy to select for disease resilience. Our objective was to evaluate mitogen stimulation assays (MSAs) on peripheral blood mononuclear cells (PBMCs) from young healthy pigs as genetic indicators for disease resilience. Data were from a natural disease challenge in which batches of 60 or 75 naïve Yorkshire × Landrace piglets were introduced every 3 wk into a continuous flow barn that was seeded with multiple diseases. In this environment, disease resilience traits, including growth, treatment, and mortality rates, were recorded on 3,136 pigs that were genotyped with a high-density marker panel. PBMCs from 882 of these pigs from 19 batches were isolated from whole blood collected prior to the disease challenge and stimulated with five mitogens: concanavalin A (ConA), phytohemagglutinin (PHA), pokeweed mitogen (PWM), lipopolysaccharide (LPS), and phorbol myristate acetate (PMA). The proliferation of cells was evaluated at 48, 72, and 96 h and compared with unstimulated samples (rest count). Heritabilities of cell proliferation were estimated using a model with batch as a fixed effect and covariates of entry age; rest count; complete blood count proportions of lymphocytes, monocytes, eosinophils, and basophils; and pen, litter, and animal genetics as random effects. Heritability estimates were highest for response to ConA (0.30 ± 0.09, 0.28 ± 0.10, 0.17 ± 0.10, and 0.25 ±0.10 at 48, 72, and 96 h after stimulation and for area under the curve across the three time points, respectively). Estimates were in a similar range for response to PHA and PMA but low for PWM and LPS. Responses to ConA, PHA, and PMA were moderately genetically correlated with several disease resilience traits and in the expected direction, but individual estimates were not significantly different from zero due to large SEs. In conclusion, although validation is needed, MSAss, in particular based on ConA, show promise as genetic indicator traits for disease resilience.

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.

Links between fecal microbiota and the response to vaccination against influenza A virus in pigs

This study describes the associations between fecal microbiota and vaccine response variability in pigs, using 98 piglets vaccinated against the influenza A virus at 28 days of age (D28) with a booster at D49. Immune response to the vaccine is measured at D49, D56, D63, and D146 by serum levels of IAV-specific IgG and assays of hemagglutination inhibition (HAI). Analysis of the pre-vaccination microbiota characterized by 16S rRNA gene sequencing of fecal DNA reveals a higher vaccine response in piglets with a richer microbiota, and shows that 23 operational taxonomic units (OTUs) are differentially abundant between high and low IAV-specific IgG producers at D63. A stronger immune response is linked with OTUs assigned to the genus Prevotella and family Muribaculaceae, and a weaker response is linked with OTUs assigned to the genera Helicobacter and Escherichia-Shigella. A set of 81 OTUs accurately predicts IAV-specific IgG and HAI titer levels at all time points, highlighting early and late associations between pre-vaccination fecal microbiota composition and immune response to the vaccine.

High Heritabilities for Antibiotic Usage Show Potential to Breed for Disease Resistance in Finishing Pigs

The use of antimicrobials in animal production is under public debate, mainly due to the risk of transfer of resistance to pathogenic bacteria in humans. Therefore, measures have been taken during the last few decades to reduce antibiotic usage in animals, for instance, by national monitoring programmes and by improving animal health management. Although some initiatives exist in molecular genetic selection, quantitative genetic selection of animals towards decreased antibiotic usage is an underexplored area to reduce antibiotic usage. However, this strategy could yield cumulative effects. In this study, we derived new phenotypes from on-farm parenteral antibiotic records from commercially grown crossbred finishing pigs used in the progeny test of Piétrain terminal sires to investigate the heritability of antibiotics usage. Parenteral antibiotic records, production parameters and pedigree records of 2238 full-sib pens from two experimental farms in Belgium between 2014 and 2020 were analysed. Heritability estimates were moderate (18–44%) for phenotypes derived from all antibiotic treatments, and low (1–15%) for phenotypes derived from treatments against respiratory diseases only. Moreover, genetic correlations between these new phenotypes and mortality were low to moderate (0.08–0.60) and no strong adverse genetic correlations with production traits were found. The high heritabilities and favourable genetic correlations suggest these new phenotypes, derived from on-farm antibiotics records, to be promising for inclusion in future pig breeding programs to breed for a decrease in antibiotics usage.

Influence of genetics and the pre-vaccination blood transcriptome on the variability of antibody levels after vaccination against Mycoplasma hyopneumoniae in pigs

Background

The impact of individual genetic and genomic variations on immune responses is an emerging lever investigated in vaccination strategies. In our study, we used genetic and pre-vaccination blood transcriptomic data to study vaccine effectiveness in pigs.

Results

A cohort of 182 Large White pigs was vaccinated against Mycoplasma hyopneumoniae (M. hyo) at weaning (28 days of age), with a booster 21 days later. Vaccine response was assessed by measuring seric M. hyo antibodies (Ab) at 0 (vaccination day), 21 (booster day), 28, 35, and 118 days post-vaccination (dpv). Inter-individual variability of M. hyo Ab levels was observed at all time points and the corresponding heritabilities ranged from 0.46 to 0.57. Ab persistence was higher in females than in males. Genome-wide association studies with a 658 K SNP panel revealed two genomic regions associated with variations of M. hyo Ab levels at 21 dpv at positions where immunity-related genes have been mapped, DAB2IP on chromosome 1, and ASAP1, CYRIB and GSDMC on chromosome 4. We studied covariations of Ab responses with the pre-vaccination blood transcriptome obtained by RNA-Seq for a subset of 82 pigs. Weighted gene correlation network and differential expression analyses between pigs that differed in Ab responses highlighted biological functions that were enriched in heme biosynthesis and platelet activation for low response at 21 dpv, innate antiviral immunity and dendritic cells for high response at 28 and 35 dpv, and cell adhesion and extracellular matrix for high response at 118 dpv. Sparse partial least squares discriminant analysis identified 101 genes that efficiently predicted divergent responders at all time points. We found weak negative correlations of M. hyo Ab levels with body weight traits, which revealed a trade-off that needs to be further explored.

Conclusions

We confirmed the influence of the host genetics on vaccine effectiveness to M. hyo and provided evidence that the pre-vaccination blood transcriptome co-varies with the Ab response. Our results highlight that both genetic markers and blood biomarkers could be used as potential predictors of vaccine response levels and more studies are required to assess whether they can be exploited in breeding programs.

Genetic parameters of animal traits associated with coccidian and nematode parasite load and growth in Scottish Blackface sheep

Gastrointestinal parasitism is a global problem for grazing ruminants which can be addressed in a sustainable way through breeding animals to be more resistant to disease. This study estimates the genetic parameters of common and new disease phenotypes associated with natural nematode and coccidian infection in Scottish Blackface sheep to underpin future genetic improvement strategies for parasite control. Data on faecal egg counts (FEC) from different species of strongyle parasites and faecal oocyst counts (FOC) from coccidian parasites were collected on 3-month-old lambs together with a faecal soiling score in the breech area dagginess (DAG) and live weight (LWT). Faecal count data were obtained for Strongyles (FEC_S), Nematodirus (FEC_N) and Coccidia (FOC). Data from 3 731 lambs sampled between 2011 and 2017 were included. Faecal egg counts and DAG records were log-transformed prior to analysis. Data were analysed using linear mixed models. Average age at sampling was 92 days with a mean LWT of 24.5 kg. Faecal soiling was not evident in 69% of lambs. Coccidia were the most prevalent parasite (99.5%), while Strongyles and Nematodirus had a prevalence of 95.4% and 72.7%, respectively. Heritability estimates (±SE) were 0.16 ± 0.03, 0.17 ± 0.03, 0.09 ± 0.03, 0.09 ± 0.03 and 0.33 ± 0.04 for FEC_S, FEC_N, FOC, DAG and LWT, respectively. Strongyles faecal egg count had a strong and positive genetic correlation with FEC_N (0.74 ± 0.09) and a moderate positive correlation with FOC (0.39 ± 0.15) while DAG was negatively genetically correlated with LWT (-0.33 ± 0.15). The significant positive genetic correlations between FEC_S, FEC_N and FOC at 3 months of age show that co-selection of sheep for resistance to these different parasites is feasible. Selection for increased resistance to parasite infection is not expected to adversely affect live BW, as no significant antagonistic genetic correlations were found between LWT and FEC. There were significant antagonistic phenotypic and genetic relationships between DAG and LWT being -0.19 ± 0.02 and -0.33 ± 0.15, respectively, indicating that the expression of the manifestation of disease in lambs may be a more meaningful indicator of the impact of parasite burden on productivity.

A 1.6 Mb region on SSC2 is associated with antibody response to classical swine fever vaccination in a mixed pig population

Classical Swine Fever (CSF) is a contagious viral disease of pigs which is endemic in several parts of the world, including India. Prophylactic vaccination using live attenuated vaccine is the preferred method of control. However, there is significant inter-individual variation in the antibody response to vaccination. In this study, we measured the E2 antibody blocking percentage after 21 days of CSF vaccination in a mixed pig population consisting of Landrace, indigenous Ghurrah pigs, and their crossbreds. A Genome Wide Association Study (GWAS) carried out using single-SNP and haplotype based methods detected a 1.6 Mb region on SSC2 (28.92-30.52 Mb) as significantly associated with antibody response to CSF vaccination. The significant region and 1 Mb flanking sequences encompass 3 genes - EIF3M, DNAJC24 and ARL14EP, which code for proteins involved in Pestivirus replication and host immune response system. Our results combined with previous studies on immune response of pigs present this region as a suitable candidate for future functional investigations.

Improvement of Disease Resistance in Livestock: Application of Immunogenomics and CRISPR/Cas9 Technology

Disease occurrence adversely affects livestock production and animal welfare, and have an impact on both human health and public perception of food-animals production. Combined efforts from farmers, animal scientists, and veterinarians have been continuing to explore the effective disease control approaches for the production of safe animal-originated food. Implementing the immunogenomics, along with genome editing technology, has been considering as the key approach for safe food-animal production through the improvement of the host genetic resistance. Next-generation sequencing, as a cutting-edge technique, enables the production of high throughput transcriptomic and genomic profiles resulted from host-pathogen interactions. Immunogenomics combine the transcriptomic and genomic data that links to host resistance to disease, and predict the potential candidate genes and their genomic locations. Genome editing, which involves insertion, deletion, or modification of one or more genes in the DNA sequence, is advancing rapidly and may be poised to become a commercial reality faster than it has thought. The clustered regulatory interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein 9 (Cas9) [CRISPR/Cas9] system has recently emerged as a powerful tool for genome editing in agricultural food production including livestock disease management. CRISPR/Cas9 mediated insertion of NRAMP1 gene for producing tuberculosis resistant cattle, and deletion of CD163 gene for producing porcine reproductive and respiratory syndrome (PRRS) resistant pigs are two groundbreaking applications of genome editing in livestock. In this review, we have highlighted the technological advances of livestock immunogenomics and the principles and scopes of application of CRISPR/Cas9-mediated targeted genome editing in animal breeding for disease resistance.

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.

The genetic basis of natural antibody titers of young healthy pigs and relationships with disease resilience

Background

Disease resilience is the ability to maintain performance under pathogen exposure but is difficult to select for because breeding populations are raised under high health. Selection for resilience requires a trait that is heritable, easy to measure on healthy animals, and genetically correlated with resilience. Natural antibodies (NAb) are important parts of the innate immune system and are found to be heritable and associated with disease susceptibility in dairy cattle and poultry. Our objective was to investigate NAb and total IgG in blood of healthy, young pigs as potential indicator traits for disease resilience.

Results

Data were from Yorkshire x Landrace pigs, with IgG and IgM NAb (four antigens) and total IgG measured by ELISA in blood plasma collected ~ 1 week after weaning, prior to their exposure to a natural polymicrobial challenge. Heritability estimates were lower for IgG NAb (0.12 to 0.24, + 0.05) and for total IgG (0.19 + 0.05) than for IgM NAb (0.33 to 0.53, + 0.07) but maternal effects were larger for IgG NAb (0.41 to 0.52, + 0.03) and for total IgG (0.19 + 0.05) than for IgM NAb (0.00 to 0.10, + 0.04). Phenotypically, IgM NAb titers were moderately correlated with each other (average 0.60), as were IgG NAb titers (average 0.42), but correlations between IgM and IgG NAb titers were weak (average 0.09). Phenotypic correlations of total IgG were moderate with NAb IgG (average 0.46) but weak with NAb IgM (average 0.01). Estimates of genetic correlations among NAb showed similar patterns but with small SE, with estimates averaging 0.76 among IgG NAb, 0.63 among IgM NAb, 0.17 between IgG and IgM NAb, 0.64 between total IgG and IgG NAb, and 0.13 between total IgG and IgM NAb. Phenotypically, pigs that survived had slightly higher levels of NAb and total IgG than pigs that died. Genetically, higher levels of NAb tended to be associated with greater disease resilience based on lower mortality and fewer parenteral antibiotic treatments. Genome-wide association analyses for NAb titers identified several genomic regions, with several candidate genes for immune response.

Conclusions

Levels of NAb in blood of healthy young piglets are heritable and potential genetic indicators of resilience to polymicrobial disease.

Maternal effect on body measurement and meat production traits in purebred Duroc pigs

We investigated maternal effect on nine body measurement traits (body height, body length, front width (FW), chest width (CW), hind width (HW), chest depth, chest girth (CHG), front cannon circumference (FCC) and rear cannon circumference (RCC)) measured at the end of performance testing and five meat production traits (ages at the start and end of performance testing (D30 and D105), average daily gain (ADG), backfat thickness and loin muscle area) in purebred Duroc pigs. Genetic parameters for each trait were estimated by using six single-trait models with and without common litter environmental effect, maternal genetic effect and direct-maternal genetic correlation. The value of Akaike's information criterion was lowest with the model including direct additive genetic and common litter environmental effects for 10 traits. The estimated proportion of common litter environmental variance to phenotypic variance was approximately ≥0.1 for D30, D105, ADG, FW, CW, HW, CHG, FCC and RCC. Using a model without common litter environmental effect would overestimate the direct heritability of most traits. Standard errors of estimated genetic parameters tended to be larger in models including maternal genetic effect. The results indicate that a compromise could be made for accurate genetic parameter estimation for body measurement traits, as well as meat production traits, in pigs by considering common litter environmental effect.

The Role of Innate Immune Response and Microbiome in Resilience of Dairy Cattle to Disease: The Mastitis Model

Animal health is affected by many factors such as metabolic stress, the immune system, and epidemiological features that interconnect. The immune system has evolved along with the phylogenetic evolution as a highly refined sensing and response system, poised to react against diverse infectious and non-infectious stressors for better survival and adaptation. It is now known that high genetic merit for milk yield is correlated with a defective control of the inflammatory response, underlying the occurrence of several production diseases. This is evident in the mastitis model where high-yielding dairy cows show high disease prevalence of the mammary gland with reduced effectiveness of the innate immune system and poor control over the inflammatory response to microbial agents. There is growing evidence of epigenetic effects on innate immunity genes underlying the response to common microbial agents. The aforementioned agents, along with other non-infectious stressors, can give rise to abnormal activation of the innate immune system, underlying serious disease conditions, and affecting milk yield. Furthermore, the microbiome also plays a role in shaping immune functions and disease resistance as a whole. Accordingly, proper modulation of the microbiome can be pivotal to successful disease control strategies. These strategies can benefit from a fundamental re-appraisal of native cattle breeds as models of disease resistance based on successful coping of both infectious and non-infectious stressors.

Genetic Basis of Blood-Based Traits and Their Relationship With Performance and Environment in Beef Cattle at Weaning

The objectives of this study were to explore the usefulness of blood-based traits as indicators of health and performance in beef cattle at weaning and identify the genetic basis underlying the different blood parameters obtained from complete blood counts (CBCs). Disease costs represent one of the main factors determining profitability in animal production. Previous research has observed associations between blood cell counts and an animal’s health status in some species. CBC were recorded from approximately 570 Angus based, crossbred beef calves at weaning born between 2015 and 2016 and raised on toxic or novel tall fescue. The calves (N = ∼600) were genotyped at a density of 50k SNPs and the genotypes (N = 1160) were imputed to a density of 270k SNPs. Genetic parameters were estimated for 15 blood and 4 production. Finally, with the objective of identifying the genetic basis underlying the different blood-based traits, genome-wide association studies (GWAS) were performed for all traits. Heritability estimates ranged from 0.11 to 0.60, and generally weak phenotypic correlations and strong genetic correlations were observed among blood-based traits only. Genome-wide association study identified ninety-one 1-Mb windows that accounted for 0.5% or more of the estimated genetic variance for at least 1 trait with 21 windows overlapping across two or more traits (explaining more than 0.5% of estimated genetic variance for two or more traits). Five candidate genes have been identified in the most interesting overlapping regions related to blood-based traits. Overall, this study represents one of the first efforts represented in scientific literature to identify the genetic basis of blood cell traits in beef cattle. The results presented in this study allow us to conclude that: (1) blood-based traits have weak phenotypic correlations but strong genetic correlations among themselves. (2) Blood-based traits have moderate to high heritability. (3) There is evidence of an important overlap of genetic control among similar blood-based traits which will allow for their use in improvement programs in beef cattle.

Exploring Phenotypes for Disease Resilience in Pigs Using Complete Blood Count Data From a Natural Disease Challenge Model

Disease resilience is a valuable trait to help manage infectious diseases in livestock. It is anticipated that improved disease resilience will sustainably increase production efficiency, as resilient animals maintain their performance in the face of infection. The objective of this study was to identify phenotypes related to disease resilience using complete blood count (CBC) data from a wean-to-finish natural disease challenge model, established to mimic the disease pressure caused by many common pathogens at the commercial level of pig production. In total, 2433 F1 crossbred (Landrace × Yorkshire) barrows that went through the natural disease challenge model were classified into four groups (resilient, average, susceptible, and dead) based on their divergent responses in terms of growth and individual treatment. Three sets of blood samples for CBC analysis were drawn at 2-weeks before, and at 2- and 6-weeks after the challenge: Blood 1, Blood 3, and Blood 4 respectively. CBC of Blood 1 taken from healthy pigs before challenge did not show differences between groups. However, resilient animals were found to be primed to initiate a faster adaptive immune response and recover earlier following infection, with greater increases of lymphocyte concentration from Blood 1 to Blood 3 and for hemoglobin concentration and hematocrit from Blood 3 to Blood 4, but a lower neutrophil concentration from Blood 3 to Blood 4 than in susceptible and dead animals (FDR < 0.05). The CBC traits in response to the challenge were found to be heritable and genetically correlated with growth and treatment, which may indicate the potential for developing CBC under disease or commercial conditions as a phenotype in commercial systems as part of developing predictions for disease resilience.

Towards the quantitative characterisation of piglets' robustness to weaning: a modelling approach

Weaning is a critical transition phase in swine production in which piglets must cope with different stressors that may affect their health. During this period, the prophylactic use of antibiotics is still frequent to limit piglet morbidity, which raises both economic and public health concerns such as the appearance of antimicrobial-resistant microbes. With the interest of developing tools for assisting health and management decisions around weaning, it is key to provide robustness indexes that inform on the animals' capacity to endure the challenges associated with weaning. This work aimed at developing a modelling approach for facilitating the quantification of piglet resilience to weaning. A total of 325 Large White pigs weaned at 28 days of age were monitored and further housed and fed conventionally during the post-weaning period without antibiotic administration. Body weight and diarrhoea scores were recorded before and after weaning, and blood was sampled at weaning and 1 week later for collecting haematological data. A dynamic model was constructed based on the Gompertz-Makeham law to describe live weight trajectories during the first 75 days after weaning, following the rationale that the animal response is partitioned in two time windows (a perturbation and a recovery window). Model calibration was performed for each animal. Our results show that the transition time between the two time windows, as well as the weight trajectories are characteristic for each individual. The model captured the weight dynamics of animals at different degrees of perturbation, with an average coefficient of determination of 0.99, and a concordance correlation coefficient of 0.99. The utility of the model is that it provides biologically meaningful parameters that inform on the amplitude and length of perturbation, and the rate of animal recovery. Our rationale is that the dynamics of weight inform on the capability of the animal to cope with the weaning disturbance. Indeed, there were significant correlations between model parameters and individual diarrhoea scores and haematological traits. Overall, the parameters of our model can be useful for constructing weaning robustness indexes by using exclusively the growth curves. We foresee that this modelling approach will provide a step forward in the quantitative characterisation of robustness.

PBMCs transcriptome profiles identified breed-specific transcriptome signatures for PRRSV vaccination in German Landrace and Pietrain pigs

Porcine reproductive and respiratory syndrome (PRRS) is a devastating viral disease affecting the swine industry worldwide. Genetic variation in host immunity has been considered as one of the potential determinants to improve the immunocompetence, thereby resistance to PRRS. Therefore, the present study aimed to investigate the breed difference in innate immune response to PRRSV vaccination between German Landrace (DL) and Pietrain (Pi) pigs. We analyzed microarray-based transcriptome profiles of peripheral blood mononuclear cells (PBMCs) collected before (0 h) and 24 h after PRRSV vaccination from purebred DL and Pi pigs with three biological replicates. In total 4,269 transcripts were identified to be differentially expressed in PBMCs in at least any of four tested contrast pairs (i.e. DL-24h vs. DL-0h, Pi-24h vs. Pi-0h, DL-0h vs. Pi-0h and DL-24h vs. Pi-24h). The number of vaccine-induced differentially expressed genes (DEGs) was much higher (2,459) in DL pigs than that of Pi pigs (291). After 24 h of PRRSV vaccination, 1,046 genes were differentially expressed in PMBCs of DL pigs compared to that of Pi (DL-24h vs. Pi-24h), indicating the breed differences in vaccine responsiveness. The top biological pathways significantly affected by DEGs of both breeds were linked to immune response functions. The network enrichment analysis identified ADAM17, STAT1, MMS19, RPA2, BAD, UCHL5 and APC as potential regulatory genes for the functional network of PRRSV vaccine response specific for DL; while FOXO3, IRF2, ADRBK1, FHL3, PPP2CB and NCOA6 were found to be the most potential hubs of Pi specific transcriptome network. In conclusion, our data provided insights of breed-specific host transcriptome responses to PRRSV vaccination which might contribute in better understanding of PPRS resistance in pigs.

Identification of resilient sows in Porcine Reproductive and Respiratory Syndrome virus infected farms

The identification of resilient sows can improve reproductive performance in farms exposed to multiple challenges. A common challenge is the porcine reproductive and respiratory syndrome virus (PRRSV). A key issue to deal with disease resilience is to set up a feasible phenotyping strategy. Our aim was to develop a phenotyping criterion to discriminate susceptible from resilient sows in PRRSV-infected farms. A total of 517 Landrace x Large White gilts were classified as resilient (R) or susceptible (S) to PRRSV virus, following vaccination with MLV-PRRSV at 6-7 wk of age, in a PRRSV negative multiplication farm. Female piglets were phenotyped as R if their serum was negative to PRRSV at 7 and 21 d post-vaccination (DPV) or as S if their serum was positive at 7 and/or 21 DPV. Amongst them, 382 gilts were transferred to a PRRSV-positive production farm, where the number of piglets born alive (NBA), stillborn (NSB), mummified (NMU), lost (NLP=NSB+NMU) and total born (NTB = NBA+NLP) were recorded for almost three years. Data were collected during two periods according to the PRRSV farm health status, which were confirmed as either PRRSV-positive stable (endemic) or inestable (epidemic). Analyses were carried out under a Bayesian approach. The heritability for the resilience criterion was estimated using a threshold model. A linear (for NTB and NBA) and a binomial model (for NSB, NMU and NLP) on the resilience criterion by the farm health status were used to assess the difference between R and S sows. The heritability of the resilience criterion was 0.46 (SD 0.06). The probability of a piglet being lost was greater (≥0.97) in S than in R litters, regardless of whether the delivery occurred during a PRRSV outbreak (20.5% vs 17.0%) or not (15.8% vs 13.7%). The lower piglet mortality rate in R sows was due to NSB, in the endemic phase (13.0% vs 15.0% of NTB, with a posterior probability of 98% of S sows showing higher NSB than R sows), and to NMU, in the epidemic phase (4.0% vs 8.4% of NTB, with a posterior probability of >99% of S sows showing higher NMU than R sows). During a PRRSV outbreak, the S sows were twice as likely to give birth to a mummified piglet as compared to R sows. These findings provide evidence that the described phenotyping scheme has a potential use as a PRRSV resilience criterion.

Invited review: Piglet survival: benefits of the immunocompetence

Piglet mortality has a negative impact on animal welfare and public acceptance. Moreover, the number of weaned piglets per sow mainly determines the profitability of piglet production. Increased litter sizes are associated with lower birth weights and piglet survival. Decreased survival rates and performance of piglets make the control of diseases and infections within pig production even more crucial. Consequently, selection for immunocompetence becomes an important key aspect within modern breeding programmes. However, the phenotypic recording of immune traits is difficult and expensive to realize within farm routines. Even though immune traits show genetic variability, only few examples exist on their respective suitability within a breeding programme and their relationships to economically important production traits. The analysis of immune traits for an evaluation of immunocompetence to gain a generally improved immune response is promising. Generally, in-depth knowledge of the genetic background of the immune system is needed to gain helpful insights about its possible incorporation into breeding programmes. Possible physiological drawbacks for enhanced immunocompetence must be considered with regards to the allocation theory and possible trade-offs between the immune system and performance. This review aims to discuss the relationships between the immunocompetence of the pig, piglet survival as well as the potential of these traits to be included into a breeding strategy for improved robustness.

Methodologies for Assessing Disease Tolerance in Pigs

Features of intensive farming can seriously threaten pig homeostasis, well-being and productivity. Disease tolerance of an organism is the adaptive ability in preserving homeostasis and at the same time limiting the detrimental impact that infection can inflict on its health and performance without affecting pathogen burden per se. While disease resistance (DRs) can be assessed measuring appropriately the pathogen burden within the host, the tolerance cannot be quantified easily. Indeed, it requires the assessment of the changes in performance as well as the changes in pathogen burden. In this paper, special attention is given to criteria required to standardize methodologies for assessing disease tolerance (DT) in respect of infectious diseases in pigs. The concept is applied to different areas of expertise and specific examples are given. The basic physiological mechanisms of DT are reviewed. Disease tolerance pathways, genetics of the tolerance-related traits, stress and disease tolerance, and role of metabolic stress in DT are described. In addition, methodologies based on monitoring of growth and reproductive performance, welfare, emotional affective states, sickness behavior for assessment of disease tolerance, and methodologies based on the relationship between environmental challenges and disease tolerance are considered. Automated Precision Livestock Farming technologies available for monitoring performance, health and welfare-related measures in pig farms, and their limitations regarding DT in pigs are also presented. Since defining standardized methodologies for assessing DT is a serious challenge for biologists, animal scientists and veterinarians, this work should contribute to improvement of health, welfare and production in pigs.

Toward the development of a one-dose classical swine fever subunit vaccine: antigen titration, immunity onset, and duration of immunity

Highly contagious classical swine fever (CSF) remains a major trade and health problem in the pig industry, resulting in large economic losses worldwide. In CSF-endemic countries, attenuated CSF virus (CSFV) vaccines have been routinely used to control the disease. However, eradication of CSFV in a geographical area would require permanent reduction to zero presence of the virus. It is therefore of paramount importance to develop a safe, potent, and non-infectious CSF vaccine. We have previously reported on a cost-effective CSF E2 subunit vaccine, KNB-E2, which can protect against CSF symptoms in a single dose containing 75 µg of recombinant CSFV glycoprotein E2. In this study, we report on a series of animal studies undertaken to elucidate further the efficacy of KNB-E2. We found that pigs vaccinated with a single KNB-E2 dose containing 25 µg of recombinant CSFV glycoprotein E2 were protected from clinical symptoms of CSF. In addition, KNB-E2-mediated reduction of CSF symptoms was observed at two weeks post-vaccination and the vaccinated pigs continued to exhibit reduced CSF clinical signs when virus challenged at two months and four months post-vaccination. These results suggest that KNB-E2 effectively reduces CSF clinical signs, indicating the potential of this vaccine for safely minimizing CSF-related losses.

Immune-associated traits measured in milk of Holstein-Friesian cows as proxies for blood serum measurements

Previous work has highlighted that immune-associated (IA) traits measurable in blood are associated with health, productivity, and reproduction in dairy cows. The aim of the present study was to determine relationships between IA traits measured in blood serum and those simultaneously measured in milk as well as their association with disease phenotypes. All animals were Holstein-Friesian cows from the Langhill research herd (n = 546) housed at the SRUC Dairy Research Centre in Scotland. Milk and serum samples were collected on 20 separate occasions between July 2010 and March 2015 and analyzed by ELISA for haptoglobin (Hp), tumor necrosis factor-α (TNF-α), and natural antibodies binding keyhole limpet hemocyanin (NAb_KLH) and lipopolysaccharide (NAb_LPS). Data were analyzed using mixed linear models that included pedigree information. Analyses revealed positive phenotypic correlations between milk and serum NAb (0.59 ≤ r ≤ 0.77), Hp (r = 0.37), and TNF-α (r = 0.12). Milk and serum NAb were also found to have a strong genetic correlation (0.81 ≤ r ≤ 0.94) and were genetically correlated with cow lameness (0.66 and 0.79 for milk NAb_KLH and serum NAb_LPS, respectively). Clinical mastitis was found to be phenotypically correlated with both milk and serum Hp (0.09 ≤ r ≤ 0.23). Serum Hp was also strongly genetically correlated with other cellular IA traits such as percent NKp46⁺ (a natural killer cell marker; 0.35) and percent peripheral blood mononuclear cells (PBMC; -0.90). Similarly, genetic correlations were found to exist between serum TNF-α and percent NKp46⁺ (0.22), percent PBMC (0.41), and percent lymphocytes (0.47). Excluding serum Hp, all milk and serum IA traits were repeatable, ranging from 0.11 (milk Hp) to 0.43 (serum NAb_LPS). Between-animal variation was highest in milk and serum NAb (0.34-0.43) and significant estimates of heritability were also observed in milk and serum NAb (0.17-0.37). Our findings show that certain IA traits, such as NAb_KLH and NAb_LPS, found in milk and serum are strongly correlated and highlight the potential of using routinely collected milk samples as a less invasive and cost-effective source of informative data for predictive modeling of animal IA traits.

Deciphering the genetic regulation of peripheral blood transcriptome in pigs through expression genome-wide association study and allele-specific expression analysis

BACKGROUND

Efforts to improve sustainability in livestock production systems have focused on two objectives: investigating the genetic control of immune function as it pertains to robustness and disease resistance, and finding predictive markers for use in breeding programs. In this context, the peripheral blood transcriptome represents an important source of biological information about an individual's health and immunological status, and has been proposed for use as an intermediate phenotype to measure immune capacity. The objective of this work was to study the genetic architecture of variation in gene expression in the blood of healthy young pigs using two approaches: an expression genome-wide association study (eGWAS) and allele-specific expression (ASE) analysis.

RESULTS

The blood transcriptomes of 60-day-old Large White pigs were analyzed by expression microarrays for eGWAS (242 animals) and by RNA-Seq for ASE analysis (38 animals). Using eGWAS, the expression levels of 1901 genes were found to be associated with expression quantitative trait loci (eQTLs). We recovered 2839 local and 1752 distant associations (Single Nucleotide Polymorphism or SNP located less or more than 1 Mb from expression probe, respectively). ASE analyses confirmed the extensive cis-regulation of gene transcription in blood, and revealed allelic imbalance in 2286 SNPs, which affected 763 genes. eQTLs and ASE-genes were widely distributed on all chromosomes. By analyzing mutually overlapping eGWAS results, we were able to describe putative regulatory networks, which were further refined using ASE data. At the functional level, genes with genetically controlled expression that were detected by eGWAS and/or ASE analyses were significantly enriched in biological processes related to RNA processing and immune function. Indeed, numerous distant and local regulatory relationships were detected within the major histocompatibility complex region on chromosome 7, revealing ASE for most class I and II genes.

CONCLUSIONS

This study represents, to the best of our knowledge, the first genome-wide map of the genetic control of gene expression in porcine peripheral blood. These results represent an interesting resource for the identification of genetic markers and blood biomarkers associated with variations in immunity traits in pigs, as well as any other complex traits for which blood is an appropriate surrogate tissue.

Transcriptome profile of lung dendritic cells after in vitro porcine reproductive and respiratory syndrome virus (PRRSV) infection

The porcine reproductive and respiratory syndrome (PRRS) is an infectious disease that leads to high financial and production losses in the global swine industry. The pathogenesis of this disease is dependent on a multitude of factors, and its control remains problematic. The immune system generally defends against infectious diseases, especially dendritic cells (DCs), which play a crucial role in the activation of the immune response after viral infections. However, the understanding of the immune response and the genetic impact on the immune response to PRRS virus (PRRSV) remains incomplete. In light of this, we investigated the regulation of the host immune response to PRRSV in porcine lung DCs using RNA-sequencing (RNA-Seq). Lung DCs from two different pig breeds (Pietrain and Duroc) were collected before (0 hours) and during various periods of infection (3, 6, 9, 12, and 24 hours post infection (hpi)). RNA-Seq analysis revealed a total of 20,396 predicted porcine genes, which included breed-specific differentially expressed immune genes. Pietrain and Duroc infected lung DCs showed opposite gene expression courses during the first time points post infection. Duroc lung DCs reacted more strongly and distinctly than Pietrain lung DCs during these periods (3, 6, 9, 12 hpi). Additionally, cluster analysis revealed time-dependent co-expressed groups of genes that were involved in immune-relevant pathways. Key clusters and pathways were identified, which help to explain the biological and functional background of lung DCs post PRRSV infection and suggest IL-1β1 as an important candidate gene. RNA-Seq was also used to characterize the viral replication of PRRSV for each breed. PRRSV was able to infect and to replicate differently in lung DCs between the two mentioned breeds. These results could be useful in investigations on immunity traits in pig breeding and enhancing the health of pigs.

Intestinal toxicity of the type B trichothecene mycotoxin fusarenon-X: whole transcriptome profiling reveals new signaling pathways

The few data available on fusarenon-X (FX) do not support the derivation of health-based guidance values, although preliminary results suggest higher toxicity than other regulated trichothecenes. Using histo-morphological analysis and whole transcriptome profiling, this study was designed to obtain a global view of the intestinal alterations induced by FX. Deoxynivalenol (DON) served as a benchmark. FX induced more severe histological alterations than DON. Inflammation was the hallmark of the molecular toxicity of both mycotoxins. The benchmark doses for the up-regulation of key inflammatory genes by FX were 4- to 45-fold higher than the previously reported values for DON. The transcriptome analysis revealed that both mycotoxins down-regulated the peroxisome proliferator-activated receptor (PPAR) and liver X receptor - retinoid X receptor (LXR-RXR) signaling pathways that control lipid metabolism. Interestingly, several pathways, including VDR/RXR activation, ephrin receptor signaling, and GNRH signaling, were specific to FX and thus discriminated the transcriptomic fingerprints of the two mycotoxins. These results demonstrate that FX induces more potent intestinal inflammation than DON. Moreover, although the mechanisms of toxicity of both mycotoxins are similar in many ways, this study emphasize specific pathways targeted by each mycotoxin, highlighting the need for specific mechanism-based risk assessments of Fusarium mycotoxins.

Acute social stress-induced immunomodulation in pigs high and low responders to ACTH

Pig husbandry is known as an intensive breeding system, piglets being submitted to multiple stressful events such as early weaning, successive mixing, crowding and shipping. These stressors are thought to impair immune defences and might contribute, at least partly, to the prophylactic use of antibiotics. Robustness was recently defined as the ability of an individual to express a high-production potential in a wide variety of environmental conditions. Increasing robustness thus appears as a valuable option to improve resilience to stressors and could be obtained by selecting piglets upon their adrenocortical activity. In this study, we aimed at depicting the consequences of an acute social stress on the immune capacity of piglets genetically selected upon divergent hypothalamic-pituitary-adrenocortical (HPA) axis activity. For this purpose, we monitored neuroendocrine and immune parameters, in high- (HPA^hi) and low- (HPA^lo) responders to ACTH, just before and immediately after a one-hour mixing with unfamiliar conspecifics. As expected, stressed piglets displayed higher levels of circulating cortisol and norepinephrine. Blood cell count analysis combined to flow cytometry revealed a stress-induced leukocyte mobilization in the bloodstream with a specific recruitment of CD8α⁺ lymphocytes. Besides, one-hour mixing decreased LPS-induced IL-8 and TNFα secretions in whole-blood assays (WBA) and reduced mononuclear cell phagocytosis. Altogether, our data demonstrate that acute social stress alters immune competence of piglets from both groups, and bring new insights in favour of good farming practices. While for most parameters high- and low-responders to ACTH behaved similarly, HPA^hi piglets displayed higher number of CD4⁺ CD8α^- T cells, as well as increased cytokine production in WBA (LPS-induced TNFα and PIL-induced IL-8), which could confer them increased resistance to pathogens. Finally, a principal component analysis including all parameters highlighted that overall stress effects were less pronounced on piglets with a strong HPA axis. Thus, selection upon adrenocortical axis activity seems to reduce the magnitude of response to stress and appears as a good tool to increase piglet robustness.

Cellular and humoral immunity in a wild mammal: Variation with age & sex and association with overwinter survival

Immune defenses are expected to be crucial for survival under the considerable parasite pressures experienced by wild animals. However, our understanding of the association between immunity and fitness in nature remains limited due to both the complexity of the vertebrate immune system and the often‐limited availability of immune reagents in nonmodel organisms. Here, we use methods and reagents developed by veterinary researchers for domestic ungulates on blood samples collected from a wild Soay sheep population, to evaluate an unusually broad panel of immune parameters. Our evaluation included different innate and acquired immune cell types as well as nematode parasite‐specific antibodies of different isotypes. We test how these markers correlate with one another, how they vary with age‐group and sex, and, crucially, whether they predict overwinter survival either within or among demographic groups. We found anticipated patterns of variation in markers with age, associated with immune development, and once these age trends were accounted for, correlations among our 11 immune markers were generally weak. We found that females had higher proportions of naïve T cells and gamma–delta T cells than males, independent of age, while our other markers did not differ between sexes. Only one of our 11 markers predicted overwinter survival: sheep with higher plasma levels of anti‐nematode IgG antibodies were significantly more likely to survive the subsequent high mortality winter, independent of age, sex, or weight. This supports a previous finding from this study system using a different set of samples and shows that circulating antibody levels against ecologically relevant parasites in natural systems represent an important parameter of immune function and may be under strong natural selection. Our data provide rare insights into patterns of variation among age‐ and sex groups in different T‐cell subsets and antibody levels in the wild, and suggest that certain types of immune response—notably those likely to be repeatable within individuals and linked to resistance to ecologically relevant parasites—may be most informative for research into the links between immunity and fitness under natural conditions.

MicroRNA Transcriptome of Poly I:C-Stimulated Peripheral Blood Mononuclear Cells Reveals Evidence for MicroRNAs in Regulating Host Response to RNA Viruses in Pigs

MicroRNAs (miRNAs) are one family of small noncoding RNAs that function to modulate the activity of specific mRNA targets in animals. To understand the role of miRNAs in regulating genes involved in the host immune response to RNA viruses, we profiled and characterized the miRNAs of swine peripheral blood mononuclear cells (PBMC) stimulated with poly I:C, a synthetic dsRNA analog, by miRNA-sequencing (miRNA-seq). We identified a total of 905 miRNAs, of which 503 miRNAs were firstly exploited herein with no annotation in the latest miRBase 21.0. Expression analysis demonstrated that poly I:C stimulation can elicit significantly differentially expressed (DE) miRNAs in Dapulian (n = 20), one Chinese indigenous breed, as well as Landrace (n = 23). By integrating the mRNA expression profiles of the same sample with miRNA profiles, we carried out function analyses of the target genes of these DE miRNAs, with the results indicating that target genes were most enriched in some immune-related pathways and gene ontology (GO) terms, suggesting that DE miRNAs play an important role in the regulation of host to poly I:C stimulation. Furthermore, we also detected 43 and 61 significantly DE miRNAs between the two breeds in the control sample groups and poly I:C stimulation groups, respectively, which may be involved in regulation of the different characteristics of the two breeds. This study describes for the first time the PBMC miRNA transcriptomic response to poly I:C stimulation in pigs, which not only contributes to a broad view of the pig miRNAome but improves our understanding of miRNA function in regulating host immune response to RNA viruses.

Resistance to infectious diseases is a heritable trait in rabbits

Selection for disease resistance is a powerful way to improve the health status of herds and to reduce the use of antibiotics. The objectives of this study were to estimate 1) the genetic parameters for simple visually assessed disease syndromes and for a composite trait of resistance to infectious disease including all syndromes and 2) their genetic correlations with production traits in a rabbit population. Disease symptoms were recorded in the selection herds of 2 commercial paternal rabbit lines during weighing at the end of the test (63 and 70 d of age, respectively). Causes of mortality occurring before these dates were also recorded. Seven disease traits were analyzed: 3 elementary traits visually assessed by technicians on farm (diarrhea, various digestive syndromes, and respiratory syndromes), 2 composite traits (all digestive syndromes and all infectious syndromes), and 2 mortality traits (digestive mortality and infectious mortality). Each animal was assigned only 1 disease trait, corresponding to the main syndrome ( = 153,400). Four production traits were also recorded: live weight the day before the end of test on most animals ( = 137,860) and cold carcass weight, carcass yield, and perirenal fat percentage of the carcass on a subset of slaughtered animals ( = 13,765). Records on both lines were analyzed simultaneously using bivariate linear animal models after validation of consistency with threshold models applied to logit-transformed traits. The heritabilities were low for disease traits, from 0.01 ± 0.002 for various digestive syndromes to 0.04 ± 0.004 for infectious mortality, and moderate to high for production traits. The genetic correlations between digestive syndromes were high and positive, whereas digestive and respiratory syndromes were slightly negatively correlated. The genetic correlations between the composite infectious disease trait and digestive or respiratory syndromes were moderate. Genetic correlations between disease and production traits were favorable. Our results indicate that it is possible to select rabbits using visually assessed disease syndromes without the need for a trade-off between health and production traits. Using a composite criterion that includes all infectious syndromes is easy to implement and heritable and is, therefore, a promising way to improve the general disease resistance in livestock species.

Evolutionary Diversity of Suid Herpesvirus 1 Based on Ul44 Partial Sequences

OBJECTIVE

The aim of this study was to use partial Ul44 sequences (glycoprotein C) of Suid herpesvirus 1 to examine the evolution and dynamics of the virus in different periods and hosts.

METHODS

Phylogenetic trees were constructed using the software MrBayes after analysis in the software jModelTest to evaluate the best phylogenetic models. The software SplitsTree 4.0 was used to create phylogenetic networks, and the BEAST program was used to generate data on phylogeography. Replication kinetics and serum neutralization tests were applied to tree strains from different phylogenetic groups.

RESULTS

Ul44 sequences derived from domestic swine and wild swine clustered in different clades and had different selective pressures depending on the host. We found no differences in replication kinetics and serum neutralization tests in the strains tested. Data show that the evolution of herpesviruses is complex, and different genetic groups may be evolving at different rates. Ul44 is an important marker for molecular evolution and epidemiology studies, but it is not useful for biological information.

Genetic relationships between measures of sexual development, boar taint, health, and aggressiveness in pigs

Breeding intact boars is a promising alternative to surgical castration of piglets. Genetic selection should enable farmers to solve problems due to boar taint and aggressiveness while taking into account potential consequences on other traits of interest. The aim of the study was to estimate genetic relations between sexual development, boar taint, health, and aggressiveness. About 1,600 Pietrain (purebred) or Pietrain × Large White (crossbred) boars were raised in a testing station. Blood samples were collected at about 105 kg BW for measuring sex hormones (testosterone and estradiol) and indicators of the inflammatory status (C-reactive protein [CRP], pig major acute-phase protein [pigMAP], and blood formula). Animals were slaughtered 9 d later and measured for boar taint compounds present in fat (androstenone and skatole) and skin lesions on carcass, an indicator of aggressiveness. For both genetic types, heritability was moderate for sex hormones (from 0.17 to 0.29) and skatole (0.24 for purebred and 0.37 for crossbred) and high for androstenone (0.63 and 0.70 for purebred and crossbred, respectively). Genetic correlations between sex hormones and boar taint compounds were moderate to high (from 0.31 to 0.95). Heritability was moderate for CRP (0.24 and 0.46 for purebred and crossbred, respectively) and very low for pigMAP (0.06 and 0.05 for purebred and crossbred, respectively. Numbers of leukocytes had moderate to high heritabilities according to the genetic type (from 0.21 to 0.52). Heritability of skin lesions was moderate for both genetic types (0.31). Genetic correlations were negative between sex hormones and inflammatory measures (from -0.46 to -0.05), positive between testosterone and number of lesions (0.43 and 0.53 for purebred and crossbred, respectively), and low between androstenone and lesions (-0.06 and -0.17 for purebred and crossbred, respectively). Overall, both breeds of pigs had very similar estimations of heritabilities, but estimates of genetic correlations were different for some pairs of traits. It would be possible to select boars based on their plasma concentration of sex hormones to decrease boar taint and aggressiveness without important consequences on the immune response. However, because of the strong links between boar taint and reproductive function, the possible consequences on the reproductive performance should be evaluated.

Transcriptomic Analysis Identifies Candidate Genes and Gene Sets Controlling the Response of Porcine Peripheral Blood Mononuclear Cells to Poly I:C Stimulation

Polyinosinic-polycytidylic acid (poly I:C), a synthetic dsRNA analog, has been demonstrated to have stimulatory effects similar to viral dsRNA. To gain deep knowledge of the host transcriptional response of pigs to poly I:C stimulation, in the present study, we cultured and stimulated peripheral blood mononuclear cells (PBMC) of piglets of one Chinese indigenous breed (Dapulian) and one modern commercial breed (Landrace) with poly I:C, and compared their transcriptional profiling using RNA-sequencing (RNA-seq). Our results indicated that poly I:C stimulation can elicit significantly differentially expressed (DE) genes in Dapulian (g = 290) as well as Landrace (g = 85). We also performed gene set analysis using the Gene Set Enrichment Analysis (GSEA) package, and identified some significantly enriched gene sets in Dapulian (g = 18) and Landrace (g = 21). Most of the shared DE genes and gene sets were immune-related, and may play crucial rules in the immune response of poly I:C stimulation. In addition, we detected large sets of significantly DE genes and enriched gene sets when comparing the gene expression profile between the two breeds, including control and poly I:C stimulation groups. Besides immune-related functions, some of the DE genes and gene sets between the two breeds were involved in development and growth of various tissues, which may be correlated with the different characteristics of the two breeds. The DE genes and gene sets detected herein provide crucial information towards understanding the immune regulation of antiviral responses, and the molecular mechanisms of different genetic resistance to viral infection, in modern and indigenous pigs.

Maternal body size influences offspring immune configuration in an oviparous snake

Like most ectothermic vertebrates, keelback snakes (Tropidonophis mairii) do not exhibit parental care. Thus, offspring must possess an immune system capable of dealing with challenges such as pathogens, without assistance from an attendant parent. We know very little about immune system characteristics of neonatal reptiles, including the magnitude of heritability and other maternal influences. To identify sources of variation in circulating white blood cell (WBC) concentrations and differentials, we examined blood smears from 246 hatchling snakes and their field-caught mothers. WBC concentrations were lower in hatchlings than in adults, and hatchlings had more basophils and fewer azurophils than adults. A hatchling keelback's WBC differential was also influenced by its sex and body size. Although hatchling WBC measures exhibited negligible heritability, they were strongly influenced by maternal body size and parasite infection (but not by maternal body condition, relative clutch mass or time in captivity). Larger mothers produced offspring with more azurophils and fewer lymphocytes. The mechanisms and consequences of WBC variation are currently unknown, but if these maternal effects enhance offspring fitness, the impact of maternal body size on reproductive success may be greater than expected simply from allometric increases in the numbers and sizes of progeny.

Time course of the response to ACTH in pig: biological and transcriptomic study

Background

HPA axis plays a major role in physiological homeostasis. It is also involved in stress and adaptive response to the environment. In farm animals in general and specifically in pigs, breeding strategies have highly favored production traits such as lean growth rate, feed efficiency and prolificacy at the cost of robustness. On the hypothesis that the HPA axis could contribute to the trade-off between robustness and production traits, we have designed this experiment to explore individual variation in the biological response to the main stress hormone, cortisol, in pigs. We used ACTH injections to trigger production of cortisol in 120 juvenile Large White (LW) pigs from 28 litters and the kinetics of the response was measured with biological variables and whole blood gene expression at 4 time points. A multilevel statistical analysis was used to take into account the longitudinal aspect of the data.

Results

Cortisol level reached its peak 1 h after ACTH injection. White blood cell composition was modified with a decrease of lymphocytes and monocytes and an increase of granulocytes (FDR<0.05). Basal level of cortisol was correlated with birth and weaning weights. Microarray analysis identified 65 unique genes of which expression responded to the injection of ACTH (adjusted P<0.05). These genes were classified into 4 clusters with distinctive kinetics in response to ACTH injection. The first cluster identified genes strongly correlated to cortisol and previously reported as being regulated by glucocorticoids. In particular, DDIT4, DUSP1, FKBP5, IL7R, NFKBIA, PER1, RGS2 and RHOB were shown to be connected to each other by the glucocorticoid receptor NR3C1. Most of the differentially expressed genes that encode transcription factors have not been described yet as being important in transcription networks involved in stress response. Their co-expression may mean co-regulation and they could thus provide new patterns of biomarkers of the individual sensitivity to cortisol.

Conclusions

We identified 65 genes as biological markers of HPA axis activation at the gene expression level. These genes might be candidates for a better understanding of the molecular mechanisms of the stress response.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2118-8) contains supplementary material, which is available to authorized users.