Next Article in Journal
Development of an Immunochromatographic Test with Recombinant MIC2-MIC3 Fusion Protein for Serological Detection of Toxoplasma gondii
Previous Article in Journal
Polyhexamethylene Biguanide Nanoparticles Inhibit Biofilm Formation by Mastitis-Causing Staphylococcus aureus
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Fine Mapping Identifies Candidate Genes Associated with Swine Inflammation and Necrosis Syndrome

1
Clinic for Swine—Herd Health Management and Molecular Diagnostics, Justus-Liebig-University Giessen, Frankfurter Strasse 112, 35392 Giessen, Germany
2
UEG Hohenlohe, Am Wasen 20, 91567 Herrieden, Germany
*
Author to whom correspondence should be addressed.
Vet. Sci. 2025, 12(5), 508; https://doi.org/10.3390/vetsci12050508
Submission received: 20 March 2025 / Revised: 30 April 2025 / Accepted: 20 May 2025 / Published: 21 May 2025

Simple Summary

Swine inflammation and necrosis syndrome (SINS) is a widespread disease in pigs that causes pain, suffering, and damage, compromising animal health, performance, and consumer demand. In addition to husbandry, feeding, and management, genetics plays a crucial role in the development of and susceptibility to SINS. The aim of this study was to identify functional gene variants in the vicinity of known gene markers to better understand the development of the disease and to explore possibilities for targeted selection against SINS. To this end, DNA samples from 234 already phenotyped piglets were sequenced at a higher resolution. The study revealed 15 significant missense polymorphisms that have functional significance. The genotypes were associated with significant differences in the SINS phenotype (SINS score 3.5 vs. 17.9). These results provide the first evidence of functional genetic markers associated with SINS and open up the possibility of genetic selection against the disease.

Abstract

Swine inflammation and necrosis syndrome (SINS) is a widespread disease in pigs, causing pain, suffering, and damage. Inflammation is documented at different levels based on clinical signs, histopathology, clinical chemistry, metabolomics and transcriptomics. The influence of sow and boar, as well as a heritability of around 0.3, suggest a genetic component to the disease. The aim of the present study was to identify functional single nucleotide polymorphisms (SNPs) in the vicinity of gene markers previously mapped using GWAS. DNA samples were available from 234 already phenotyped piglets. These animals were re-sequenced with additional prior enrichment. The nine selected chromosomal regions cover a total length of 22 Mbp. The genome-wide association study (GWAS) revealed two series with a total of 15 significant missense polymorphisms on chromosomes 11, 14, and 15. The homozygous genotypes of the most discriminating SNPs in series 1 resulted in SINS scores of 3.5 and 17.9, respectively. Despite the partial linkage of the SNPs, interesting candidate genes were defined. The results allow a significant narrowing of the possible candidate genes for understanding the pathogenesis of SINS and for future use in selection breeding to overcome the syndrome. Further studies should be carried out on larger animal populations.

1. Introduction

Improving animal welfare has become an important goal, especially for developed societies [1]. A newly defined but common disease in pigs with significant welfare implications is swine inflammation and necrosis syndrome (SINS) [2,3,4]. Symptoms begin with inflammatory loss of bristles, redness, and swelling. Later, exudative inflammation and even necrosis may occur [2]. Bleeding is also possible. The base of the tail, tip of the tail, ears, face, teats, navel, coronary bands, heels, and claws of newborn and suckling piglets, weaners, and finishers are affected. Macrophages and lymphocytes in the inflamed tails of piglets only a few hours old with intact epidermis show that SINS occurs without biting or technopathy, confirming the endogenous nature of the syndrome [5], although it can be modified by environmental conditions [2]. This endogenous component is also indicated by heritabilities around 0.3 for SINS traits [3]. Furthermore, significant and replicable breed and boar effects have been found [3,6]. SINS is widespread. Studies from Germany, France, and the Netherlands show that SINS affects over 60% of piglets on affected farms [2,3,4,7]. The associated inflammation and necrosis generally cause pain, suffering, and damage to the affected animals [8,9]. Improvements in housing conditions, in particular, the provision of hygienic drinking water from drinking troughs, the creation of opportunities for efficient thermoregulation, the reduction of protein and starch content in favor of crude fiber in feed, and the control of mycotoxins, can help control the syndrome [2,10], but are often difficult to implement in practice.
At the same time, detectable boar effects and moderate inheritance suggest the possibility of selecting against SINS and overcoming the syndrome [2,3]. Using boars that are less susceptible to SINS can result in higher birth and weaning weights and more piglets. At the same time, piglets experience less inflammation, necrosis, biting, and skin damage. This results in significantly improved animal welfare [3]. It is important to look for clues about the genes involved, not only for targeted selection based on functional markers, but also to better understand the pathogenesis of SINS. Studies of the liver transcriptome of piglets with varying degrees of SINS [11,12] allowed the clinical and pathohistological findings to be understood at the level of gene expression. Together with clinical chemistry and metabolomic results [11,13], they confirm a switch in the liver from normal anabolic metabolism to inflammatory and acute phase metabolism in SINS. A first genome-wide association study identified significant but not yet functional SNPs [14]. The significant SNPs associated with SINS traits were distributed across the entire genome with particularly high effects on SSC 7, 9, 11, 12, and 14–17. However, the relatively low coverage of this exploratory study has not yet revealed any functional SNPs.
The aim of the present study was to re-genotype selected regions of the piglets studied by Gerhards et al. [14] at a higher sequencing coverage to search for functional SNPs. Nine regions on eight chromosomes with the highest associations with SINS were selected.

2. Materials and Methods

2.1. Study Design

The study was based on the phenotypes and DNA of 234 piglets from the work of [14]. The phenotypes were used again because they were collected in a complex and precise manner, and no further experimental animals had to be used. The already published animal experiment was carried out in the conventional pig breeding stable of the Oberer Hardthof teaching and research station at Justus Liebig University Giessen under the approval of the authorities in Giessen, Germany with file number 708_M.
The environment was the same for all piglets studied. Sows and suckling piglets were housed together in 4.8 m2 farrowing pens with a plastic slatted floor. The sows were fixed in a farrowing crate. The floor for the sows was a cast iron slatted floor. Nipple drinkers and mother-child basin drinkers ensured the water supply for the animals. The lactating sows were fed according to their specific needs with a spotmix feeder in the trough.
The sow herd was a rotational cross of Topigs with German Landrace. The proportions of the individual genetics in the different sows were not the subject of detailed characterization. On the boar side, three purebred animals were selected that had been characterized for SINS on the basis of their progeny [6]. A Pietrain boar (PI) whose progeny had comparatively high SINS scores was designated high SINS PI (PI−). A second Pietrain boar whose progeny had comparatively low SINS scores was designated low SINS PI (PI+), and a Duroc boar whose progeny had the lowest SINS scores was designated as a low SINS Duroc boar (DU). The quantification of the SINS sensitivity of the boars results from [6]. There was no genetic relationship between these boars. The herd had a performance of 15 live born and 1.4 dead born piglets per sow.
The sows were artificially inseminated with the three boars that differed in terms of SINS. The boars were used in pairs as mixed semen, where the semen from two boars was mixed within one dose. This means that piglets from two different boars were present in each litter at the same time, with the high SINS Pietrain boar always acting as a positive control. This design was used to (i) limit the number of litters and experimental animals, (ii) to minimize environmental effects, (iii) to increase genetic variability within the piglets, (iv) to increase the number of sow–boar combinations, (v) to achieve a manageable number of piglets, and (vi) to minimize the variability of potential pathways involved in SINS sensitivity by using only three boars. At the same time, however, this small number of boars also led to stronger coupling of the SNPS.
The 27 sows produced 27 litters. In total, 234 anatomically normal piglets were used for further analysis. Their SINS phenotype was recorded at their third day of life. The father of the piglets was detected by paternity testing [14] (14 microsatellites in 2 multiplex PCRs, capillary gel electrophoresis) after phenotyping. This kind of order was chosen to ensure blinded data collection. The results of paternity testing revealed 14 mixed litters (from 14 sows) with 77 piglets from the PI+ and 39 piglets from the PI−, as well as 13 mixed litters from 13 sows with 48 piglets from the DU and 70 piglets from the PI−. The higher number of offspring from the PI− boar was due to the fact that it was always used as a control in the mixed semen experiment.

2.2. Clinical Scoring

The SINS phenotypes of the 234 piglets used in this study were taken from the previous study by [14].
In this study, inflammation and necrosis were recorded clinically, as already described in [5]. To ensure comparability with other studies, the piglets were scored on the day 3 of life. In all previous studies, clinical signs were clearly visible at this time, and the piglets were not yet as much exposed to environmental effects like weaners and fatteners. In order to minimize animal load, clinical signs were recorded using a digital camera (Canon EOS DC 8.1V, Canon, Krefeld, Germany) according to a standardized scheme for later detailed evaluation of the images (Windows Media Player, version 12, Microsoft GmbH, Munic, Germany).
Clinical alterations in the tail base and tail tip, the ears, the teats and navel, coronary bands, wall horn, heels, as well as the face were assessed individually. The following clinical characteristics were considered and scored 0 if the sign was not visible or 1 if the sign was visible (Table 1). The tail base and tail tip were independently scored for loss of bristles, swelling, redness, scab formation (tail tip only), rhagades, exudation, necrosis, bleeding (tail tip only), and ring-shaped constrictions (tail tip only). Ears were scored for a shiny skin, the loss of bristles, necrosis, and congested ear veins. The face was scored for the absence or presence of oedema at the eye lids and nose back. Teats were scored for swelling, reddening, scab formation, necrosis, and congested blood vessels. The navel was scored for signs of inflammation in the form of redness or swelling. Claws were scored qualitatively for any signs of inflammation at the coronary bands, walls, and heels. Scores for coronary bands and walls were averaged for front and hind limbs. Because the heels were drastically affected in some piglets, they were scored 0–3 for swelling and 0–3 for bleeding, each in front and hind limbs. The binary scores examined are presented by organ system as a percentage of affected piglets.
All findings were summed up to give an additive body part score (Table 1). This additive body part scores could range from 0 to 12 points. All of these body scores were summed in an unweighted manner to give the SINS score. This resulted in possible SINS scores between 0 and 44 for each piglet. All scores were assigned by two experienced persons using the four-eye principle. An overview of the evaluated phenotypes is given in Table 1. In addition, the SINS score after Z-transformation (ZSINS) was used.

2.3. DNA Extraction and Sequencing

DNA was extracted from docked piglet tails with the smart DNA prep (m) kit (Analytik Jena, Jena, Germany) and quantified using the Qubit dsDNA broad range assay kit (Invitrogen, Thermo Fisher Scientific, Waltham, MA, USA) on the Qubit Flex Fluorometer (Invitrogen, Thermo Fisher Scientific, Waltham, MA, USA). DNA was diluted to a uniform concentration of 50 ng/μL.
During shotgun library preparation, DNA (30–100 ng) was fragmented using the Fragmentase enzyme mix of the Allegro Targeted Genotyping Kit (TECAN/NuGEN, Männedorf, Switzerland). Indexed Illumina libraries were constructed using the Encore Rapid DR Multiplex System 1–96 (TECAN/NuGEN, Männedorf, Switzerland). All libraries were individually amplified for 8 cycles using MyTaq HS Red Mix (Meridian Bioscience, Bioline, Cincinnati, OH, USA) and standard Illumina primers in a 30 µL reaction volume. Additives were used to minimize PCR bias in library amplification (BioStab II PCR Optimizer, Sigma-Aldrich, St. Louis, MO, USA; Herculase II Fusion DNA Polymerase, Agilent, Santa Clara, CA, USA). The final libraries underwent quality control (QC) using agarose gel electrophoresis. If the concentration was slightly low, the libraries were run for up to 5 more PCR cycles. Up to 12 libraries were pooled for enrichment. Compiled pools were cleaned up using AMPure bead purification (using 0.8 Vol) followed by Monarch column purification (New England Biolabs, Ipswich, MA, USA).
Enrichment was performed using a ‘myBaits custom Enrichment kit’ (Daicel Arbor Biosciences, Ann Arbor, MI, USA), following the protocol for ‘most taxa’ and using 65 °C as the hybridization temperature. The only modification was the amplification of the enriched libraries. To avoid PCR bias, we again used the MyTaq HS Red Mix (Meridian Bioscience, Bioline, Cincinnati, OH, USA) with additives as described above. PCR reaction volumes for the enriched pool-DNAs were upscaled to a total volume of 200 µL. PCRs were run for 15 cycles.
Before sequencing, amplified library pools were purified with 0.8 Vol AMPure beads and size-selected on preparative agarose gel (350–650 bp fragments, 250–550 bp inserts).
The final QC of the DNA libraries was carried out using the Agilent Fragment Analyzer (Agilent, Santa Clara, CA, USA) and Qubit Fluorometer (Invitrogen, Thermo Fisher Scientific, Waltham, MA, USA).
Sequencing was performed on Illumina NextSeq 500 /NovaSeq 6000 (Illumina, San Diego, CA, USA). In this process, 150 bp paired-end reads with an average of 4 M read pairs per original sample were generated.
The Sus scrofa chromosomes (SSC) and regions selected for resequencing in this study were based on the SNPs with the highest probabilities of association with SINS traits identified in [14] (Table 2). These SNPs had effects of 8.1 and 8.4 (chromosome-wide) on SSC 9 and 12 and 8.9 to 16.8 (genome-wide) on SSC 7, 11, and 14 to 17. The region was defined around the corresponding SNP, and a range of 0.8 to 2.6 Mb was sequenced for each of these 9 regions.

2.4. OVarFlow Pipeline

For further analysis, the bz2 compressed fastq files received from the sequencing company were decompressed and converted to gz compressed fastq files.
The available raw data were transferred to the open-source workflow OVarFlow. OVarFlow is used for variant discovery of SNPs and insertion/deletion polymorphisms (Indels) in model and non-model organisms [15]. The workflow allows automation, documentation, and reproducibility of the individual analysis steps.
The reference genome and annotation Sscrofa11.1 (GCF_000003025.6; https://www.ncbi.nlm.nih.gov/datasets/taxonomy/9823/ (accessed on 19 May 2025)) were used for the implemented alignment, variant calling, and annotation. The min sequence length was set to 1.
The elements of the remaining analysis steps as well as their programs and versions used in the OVarFlow workflow correspond to those already described in the study by [14].

2.5. Genome-Wide Association Study (GWAS)

The preparation of the data output from OvarFlow for the GWAS and the actual execution of the GWAS was performed in R, version 4.2.1 [16]. RStudio [17] was used as the graphical user interface.
Since the enriched regions only extend to the autosomes, the gonosomes were not analyzed in the following. The annotated VCF file output at the end of the OVarFlow workflow was converted to binary PLINK format using the PLINK package version 1.9 [18,19].
The initial quality control of the genotype data was performed using PLINK. Only variants and individuals that passed these filter criteria were kept for the ongoing GWAS (missingness per marker: <0.01; missingness per individual: <0.1; minor allele frequency (MAF): >0.05; Hardy–Weinberg equilibrium (HWE): p > 0.000001).
A total of 63,362 variants were detected on the autosomes from 234 sequenced pigs.
A principal component analysis (PCA) was performed using TASSEL version 5.2 [20]. TASSEL 5.2 was also used to create a kinship matrix to account for kinship relationships from the offspring of the three different boars in further GWAS analysis.
The actual GWAS was conducted by GAPIT version 3 [21] using the BLINK model [22]. For this purpose, genotype data were converted from PLINK to hapmap format using TASSEL 5.2.
Boar effects, sow (litter), contemporary group, sex, birth weight, and PCA were used as fixed effects.
For PCA, four principal components were considered. Each sow was used only once in the study, so the effects of litter and sow were identical.
For the GWAS, Bonferroni-corrected significance thresholds were calculated genome wide (p ≤ 7.89 × 10−7) and chromosome wide (p ≤ 1.69 × 10−4–2.73 × 10−5) based on the 63,362 SNPs and a significance level of α = 0.05. All of the SNPs presented in the following overcome this significance threshold and can therefore be considered significant.

2.6. Variant Effect Prediction

For the functional annotation of the genetic variants found, information on the genetic coordinates, effect, and impact of each variant was provided. For this purpose, snpEff, version 5.0 [23] and Ensembl Variant Effect Predictor (VEP), release 113 [24] were used. In addition, the Sorting Intolerant From Tolerant (SIFT) score was given for non-synonymous mutations, which predicts whether an amino acid substitution could affect protein function based on sequence homology and the physical properties of the amino acids [25].

2.7. Statistical Analysis of SNP Effects on SINS Traits

Statistical analysis was performed using IBM SPSS version 29 (IBM, Munich, Germany). Body part scores and SINS scores were normally distributed. The analysis was performed using multivariate analysis of variance with the fixed effects of boar*piglet genotype, sex, and batch. The sow (litter) was used as a random effect in the model, and the birth weight of the piglets was used as a covariate. Interactions between effects were included. For the comparative estimation of the effects of boars and genotypes, a modified model was used instead of boar*genotype, with boar and genotype as fixed effects, while keeping all other effects the same. To characterize the boars, this model was modified by excluding the effect of the genotype. All SNPs were 100% linked in one of the two series. This linkage was cross chromosomal, i.e., it was not only haplotypes. This linkage is a consequence of the small number of boars. The genotypes of each SNP were designated as N (norm), and the deviating genotypes were designated as V (variant) according to the reference genome Sscrofa11.1. Due to the 100% agreement of the genotypes of all SNPs in a series, the individual SNPs were not evaluated further, but series 1 and 2, to which all SNPs were assigned, were evaluated. Means, standard errors (SE), and 95% confidence intervals (CI95) were presented. Prevalence was calculated by counting affected and unaffected individuals and presented in a purely descriptive manner. In ANOVA, differences were considered statistically significant at a significance level of α = 0.05. Since two series were tested according to Bonferroni, the lowest significance level was set at p ≤ 0.025. Some values that tended to differ were also shown in order not to lose any indication of possible effects.
The genotypes of the sow and boar were inferred from the genotypes of the piglets.

2.8. Annotation of Potential Candidate Genes

The NCBI Genome Data Viewer was used to search for positional candidate genes using the Sscrofa11.1 reference genome (https://www.ncbi.nlm.nih.gov/gdv/browser/genome/?id=GCF_000003025.6 (accessed on 19 May 2025)). A gene was considered to be a candidate gene if at least one significantly associated functional variant was located within the gene. Functional variants could include missense, splicing, and 5′ UTR variants.
Physiological functions of potential candidate genes were deduced using information from the NCBI database [26], GeneCards [27], PathCards [28], KEGG [29], MGI [30], and DAVID [31].

3. Results

3.1. Phenotypes

The offspring of the three boars differed significantly in their SINS scores (p < 0.001) (Figure 1). The 95% confidence intervals for the offspring of DU, PI+, and PI− were between 9.3 and 11.6, 12.1 and 13.2, and 13.2 and 14.4, respectively.
In line with the significant differences in overall SINS scores, the offspring of the three boars also differed significantly in organ scores (except for claw wall). As expected, the offspring of the DU boar always had the lowest scores and those of the PI− boar always had the highest scores (Table 3). The greatest differences were found in the scores for tail base, tail tip, and teats, where the PI− offspring scores were about twice those of the DU offspring. The PI+ progeny scores were at the level of the DU progeny (tail base, tail tip, coronary bands), at the level of the PI− progeny (ears, heels), or in between (teats). With the exception of tail tip and teat alterations in Du offspring, the prevalence ranged from over 40 to over 90%. Again, Du offspring were always less affected, whereas PI− offspring had the highest prevalence in most cases.

3.2. Single Nucleotide Polymorphisms (SNPs)

The sequenced regions covered about 22 Mbp and revealed a total of 8668 SNPs that were significantly associated with at least one SINS phenotype. Twenty-two of these SNPs were functional and were used for further analysis (Table 4). Of these SNPs, 68.2% were missense variants (Missense) associated with an amino acid exchange in the protein; 22.7% were splice variants (Splice), where altered expression might be expected; and 9.1% were variants in the 5′ untranslated region (5′ UTR), i.e., in the regulatory region in or close to the promoter of a gene. The number of such SNPs on chromosomes 11, 14, and 15 was 10, 6, and 6, respectively. Due to the selected choice of the three boars, there was a strong association between the SNPs, which were grouped into one of two series. Series 1 was associated with SNPs on SSC11, SSC14, and SSC15. The second series was restricted to SSC15. The SNPs were located within or directly in the 5′ or 3′ UTR of defined genes. The SIFT scores of the resulting missense mutations were mostly classified as tolerated (SIFT score > 0.05) and in one case as deleterious (SNP 14_75945482; NUDT13 gene; SIFT score between 0.03–0.04).
The significance of series 1 SNP effects on the phenotypes, expressed as -LOG (p), was 7.5 (GWAS maximum for heels) and 19.2 (GWAS maximum for tail tip). The significance of the series 2 SNPs was 11.3 (GWAS maximum for tail base) (Table 4).
The genotypes of the piglets allowed a clear derivation of the genotypes of the sows and boars (Table 5). For the SNPs of series 1, the DU boar was homozygous for NN, and both Pietrain boars were homozygous for VV. The VV genotype predominated in the sows. The NN genotype only occurred in one sow. The distribution of genotypes in the piglets followed Mendel’s rules. The same picture was observed in series 2 (Table 5), although here the genotype NV was found in the PI− boar. Thus, the genotypes of the SNPs were largely assigned to the boars, and the diversity of the piglet genotypes resulted mainly from the genotypes of the sows.
The favorable genotypes of all SNPs of both series occurred mainly in the progeny of the DU boar (Figure 2, Table 6). With regard to the SNPs of series 1 (Figure 2A, Table 6), the SINS values of the NN genotypes (mean = 3.5) were significantly lower than those of the NV genotypes (mean = 14.6) within the progeny of the DU boar and significantly lower than those of all other boar genotype combinations. At the same time, the scores of the PI− progeny with the VV genotype (mean = 17.9) were significantly higher than those of all other boar genotype combinations. In contrast, the SINS scores of the NV genotypes of the DU and PI boars and the VV genotypes of the PI+ boar did not differ. A similar picture emerged when considering series 2 (Figure 2B, Table 6), the only difference being that the PI− genotypes present were both associated with significantly higher SINS scores than all other boar–genotype combinations. For series 1, no difference could be observed between DU and PI− with respect to the heterozygous genotypes, whereas the VV genotypes differed significantly between PI+ and PI− progeny in favor of the PI+ piglets. In series 2, there was an additional significant disadvantage for the heterozygous genotype of the PI− offspring compared to the DU offspring. In both series, the increase from the NN to the NV genotype dominated, but this could only be reproduced in the DU offspring, as the NN genotype did not occur in the Pietrain offspring. This resulted in an overlap of boar effect and genotype, which could not be clearly separated due to the different distribution of genotypes in the boars. At the same time, however, the effect of increasing SINS scores from NN to NV to VV continued at comparable boar stages (DU, PI−).
The correlation for the SINS score was also largely reproduced at the level of individual body part scores, particularly for tail base, teats, coronary bands, and heels (Table 6).
When genotype and boar effect were simultaneously considered in the saturated model, the largest effect was found for series 1 and 2 (Figure 3). When the SINS scores of the progeny of the most favorable boar (DU) were 20% lower than those of the least favorable boar (PI−), the SINS scores of the NN genotypes were 79% lower than those of the VV genotypes in both series. The heterozygous genotypes were 16% (series 1) and 0% (series 2) below the values of the VV genotypes.

4. Discussion

Inflammation and necrosis of different body parts in pigs (SINS) often occurs at different age stages and affects the health and welfare of the animals [2,3,4,7]. In addition to husbandry and environmental influences [2,5,7], genetic influences on the expression of SINS have been demonstrated at the level of boar effects [6], breed effects [3,6], calculated heritabilities [3], and genomic associations [14] and confirmed using transcriptomic analyses [11,12]. The study by Kuehling et al. [6] showed significant differences between boars of the Duroc and Pietrain breeds in favor of the Duroc boars and between boars within the Pietrain breed. Therefore, a typical Duroc boar, whose progeny were hardly affected by SINS, and two Pietrain boars, whose progeny had the lowest (PI+) or highest (PI+) SINS score in the study by [6], were selected for the present study.
The piglets phenotyped and genotyped accordingly were already part of a genome-wide association study [14] with the aim of generating evidence of the genetic background of the extreme phenotypes. This study was closely accompanied by a transcriptome study [12]. By producing a number of offspring from these individual boars, it should in principle be possible to estimate the effects of specific SNPs. Building on the study by [14], 11 chromosomal regions with particularly high association to the SINS genotype were selected and sequenced with improved sequencing coverage compared to the original study. DNA and phenotypes were taken from the study by [14] in order to save experimental animals.
We deliberately selected only these three extreme boars for the study to reduce the possibility that different lines within a breed might use different pathways that interfere with each other and thus escape detection. However, this also means that there may well be alternative pathways for the regulation of SINS in other individuals and lines. More research is needed in this context. Of course, it was not the aim of the present study to identify differences between the breeds, and any information on differences between the progeny of the three boars can only refer to these three individual boars. Future studies on a larger number of boars are needed to make an overall statement at the breed level.
However, such studies can build directly on the results of the present study and thus save time and resources.
The pronounced chromosome-wide linkage of the SNPs in the present study is also a limitation of the chosen methodological approach. All identified functional SNPs were 100% linked in one of the two genotype series. Therefore, it is not possible to decide which of the SNPs in a series is actually functionally associated with the SINS phenotype. The reason for this is that when a boar is used, all SNPs of homozygous chromosome segments are only passed on to the offspring in a fully linked manner because crossing over has no effect. All linked SNPs will therefore show the same association with the phenotype. Future studies should therefore use several different boars to increase the variability of the SNPs adjacent to the functional SNP, thereby reducing haplotype blocks and breaking the linkage between adjacent SNPs. Nevertheless, in contrast to [14], the present approach identified the first functional SNPs associated with SINS and led for the first time to a manageable number of candidate genes that can be investigated in detail in the future. The association of any SNP in a series with the trait of interest must be seen as the gross effect of the 19 and 3 SNPs in series 1 and 2, respectively. It remains unclear how many of the SNPs are actually directly associated with the target traits and to what extent. The determination of the individual effects of the SNPs is reserved for future work, which would need to be carried out on larger numbers of boars and their progeny in order to separate the individual effects of the SNPs from each other. Future research with higher piglet and boar numbers could attempt to include animal models or haplotype-based GWAS to disentangle these effects.
It must further be taken into account that the phenotypes may be only partially explained by a certain SNP and that a residual variance of unexplained gene effects remains, which ultimately originates from unknown SNPs of the sow and boar.
The alleles of the SNP series were divided into N and V, with the N allele corresponding to the gene bank sequence (Sscrofa11.1) and the V allele to the variant deviating from it. The gene bank sequence originates from a pig of the Duroc breed. The basic assumption was that boars with extreme SINS scores should differ in the allele distribution of the SNPs associated with the phenotype. This hypothesis was confirmed. In relation to the animal species, there is also a clear association between the identified SNPs and the SINS score. However, the effects of the boars and the respective SNPs could not be completely separated. In a complete model, joining boar and SNP effects, however, the effects of the SNP series were estimated to be higher than the boar effects. The comparison of two genotypes (NN vs. NV in DU and NV vs. VV in PI−) also confirmed the directional relationship between SNPs and phenotype. However, it can be assumed that the respective boars should carry further favorable and unfavorable alleles in addition to the identified SNPs with regard to the SINS score, but these could not be detected due to the low effect sizes and the selected chromosome regions. Their identification would require very large numbers of animals.
The 19 SNPs in series 1 cover nine genes. Some of these genes meet the criteria for homologous physiological candidate genes for inflammation. Others show associations with performance parameters in swine or both. Negative side effects of selection for improved productivity that affect animal welfare have been widely reported [32,33]. The genes identified in the present study are interesting candidates for detecting such negative side effects. In addition, they may support the hypothesis that selection for performance may also have led to side effects in the form of inflammation and necrosis observed in SINS. Selection for performance can lead to favored gene variants and their associated selection along with alleles of other genes (selection signature/selective sweeps) [34,35,36]. Such genes can be physiologically unrelated to the performance trait, but involved in other traits [37], for example in the regulation of inflammation. On SSC11, three candidate genes are located within a known sweep region for growth, meatiness, and fertility in the pig.
One of the most interesting candidate genes of the present study is the protein lysine methyltransferase 21C (METTL21C), an important signaling molecule that modulates cellular plasticity in response to stress and ageing by modulating the NF-kB signaling pathway [38,39]. Previously known activities of METTL21C include heat shock protein binding [40,41] and response to oxidative damage during heat stress [38,42]. In addition, numerous studies demonstrate the role of METTL21C in myoblast formation, skeletal muscle development, and differentiation in poultry [39,43,44], which can lead to increased muscle weight and muscle fiber diameter. Steinert et al. [45] showed that the hypertrophic effect of METTL21C is partly due to the inhibition of protein degradation. By regulating the function of the uterus in laying hens, their production performance can be increased [46]. The gene is also associated with adaptability in chickens [47]. As poultry and pigs have undergone a very similar history of selection for conformation and muscle fiber thickness on the one hand, and for reproduction on the other, it is conceivable that a gene variant of METTL21C might have been selected for that has a negative effect on response to oxidative burst. The respiratory burst, which leads to the generation and accumulation of extracellular reactive oxygen species (ROS) to protect against invading pathogens, is one of the characteristic cellular processes observed during inflammation [48,49,50]. However, during oxidative stress, excessive ROS generation leads to imbalances that are associated with numerous pathologies, including cardiovascular disease [51,52,53,54]. Intracellular ROS has been shown to be involved in VEGF-dependent signaling in endothelial cells [55]. Thus, oxidative stress is directly linked to inflammation in a positive feedback loop [56]. Vasculitis plays a key role in the pathogenesis of SINS [5,10], and an intrinsic role of ROS, SOD, and VEGF has been demonstrated by transcriptomic studies [11,12,13]. In sheep and cattle, METTL21C has also been described in relation to stress signaling, oxidative metabolism, and calcium ion transport [57]. In the present study, a missense mutation and a mutation in the 5′ UTR of the SINS-associated METTL21C gene were found in series 1. There was a replacement of the small polar amino acid threonine by the non-polar amino acid methionine. When a missense mutation causes an amino acid change in a protein, it can have different effects depending on where it occurs in the protein and which amino acid is replaced by which other amino acid. It can lead to altered protein function, misfolding, and instability and be directly linked to disease; alternatively, it can be neutral and have no effect [58,59]. A base change in regulatory regions such as promoters or enhancers can affect the binding of transcription factors or other regulatory proteins. This can lead to a change in the expression of the gene, which in turn leads to a change in the amount of the responsible protein. In addition, SNPs in UTR regions can affect the stability and function of the mRNA and thus the amount of protein available [60]. Those substitutions that do not lead to a complete loss of protein function, but only to minor changes in the metabolic machinery, may be of particular interest in this context. For the coiled-coil domain containing 168 (CCDC168) gene, some partially deleterious SNPs were found in a selection scan region on SSC 11 when comparing wild boar and domestic pig sequences. This suggests that this gene is involved in a sweep mechanism, where the frequency of variants in a linked gene is strongly influenced as a result of selection for growth, meat or fertility [35]. It is discussed that CCDC168 may be associated with thermoregulatory ability in cattle and associated reduced fertility [61]. However, the functions of this gene are still largely unknown [62]. Based on the importance of thermoregulation for SINS [2] and the existing sweep effect, a candidate gene role for CCDC168 can also be considered. In the present study, six different missense mutations were identified. Among them were two substitutions of glutamic acid (acidic) by arginine (strongly basic). These substitutions are not considered deleterious, but may contribute to increased susceptibility to unfavorable thermoregulation, which plays a key role in the pathogenesis of SINS [2].
The basic immunoglobulin-like variable motif-containing protein (BIVM) gene has been associated with other selection criteria such as number of teats in Duroc pigs and improved reproduction [34,36]. In the present study, a splice variant was found.
SNPs that influence the splicing of the pre-mRNA by altering the splice site or by influencing the splicing regulators can lead to the formation of different protein variants, so-called isoforms, or to the loss of functional proteins [63]. It remains unclear how BIVM could be involved in the pathogenesis of SINS. However, the association of the gene variants with SINS susceptibility could be due to the fact that it may have arisen as a negative side-effect of selection for reproductive performance and number of teats.
ERCC5 (excision repair 5, endonuclease) is a paralog of BIVM and is involved in Nucleotide excision repair pathway [64]. ERCC5 is important for repairing DNA damage and has been linked to certain cancers and to the skin disease xeroderma pigmentosum [65].
NUDT13 (nudix-hydrolase 13) is involved in nucleotide metabolism. ECD (ecdysoneless cell cycle regulator) is involved in cell cycle regulation and has been detected in the peripheral blood transcriptome of pigs [66]. FAM149B1 (member of the family with sequence similarity 149, member B1) has been shown to be associated with renal function and dysfunction [67]. We did not find any homologous physiological overlap among ERCC5, NUDT13, ECD, FAM149B1, and SINS.
Regarding candidate genes from series 2, myosin IIIB (MYO3B) is associated with performance parameters (loin; fat ratio) in pigs [68,69]. In this context, it has been discussed that this gene may be responsible for the variation of other neighboring genes in the context of selection sweep [70]. In addition, associations with immune response were found [71]. Whether MYO3B has a direct or indirect effect on susceptibility to SINS could not be determined. None of the other candidate genes of series 2 appear to be directly related to inflammation and necrosis.
Current knowledge of the candidate genes and the pathways in which they are involved, particularly in pigs, does not yet allow a complete understanding of the genes in SINS pathogenesis. It is still unclear what effects the identified variants of the amino acid changes of the most interesting candidate genes might have on the phenotypes; however, for such key regulatory genes, we expect small functional differences rather than adverse effects on the target traits. Nevertheless, these results open up many possibilities for future research.
As the most common commercial boar breeds were used in the present study, it is likely that the identified candidate genes also play a role in other individuals of this breed. However, it is also possible that other pathways play an additional or even dominant role in other boars.
To answer this question, a larger study with significantly more boars of the same and other breeds and thus a greater variety of analyzed piglets is needed.
Based on the results of this and future studies, boars could be genotyped for the corresponding SNPs, and sires could be selected for breeding that carry as many favorable gene variants as possible. Therefore, these offspring would be less susceptible to SINS and can also pass this on to their offspring. Improving animal welfare has become an important goal, especially for developed societies [1], and breeding for increased performance has been identified as a major factor contributing to negative side effects on animal welfare [32,33]. This study provides evidence for candidate genes with gene variants that may be responsible for such negative side effects. These side effects occur in the form of SINS in a significant proportion of piglets on many farms, causing considerable pain, suffering, and damage [2]. Studies by several authors show that genetic selection for SINS is possible due to significant boar effects [6] and pronounced heritability [3], and that such selection would also result in significantly improved performance [3,4,10]. It could thus contribute to improved welfare through breeding [1,33,72,73]. It should be emphasized that this would not involve selecting for new gene variants that adapt animals to husbandry conditions, which would be critical in any case [72], but rather selecting against gene variants that have become widespread in populations due to breeding for increased performance and cause negative side effects. The candidate genes presented provide valuable approaches for confirmation in larger studies with a larger number of boars. Once favorable alleles are confirmed, they could be incorporated into practical genomic selection. Corresponding breeding values for boars could be recorded and economically weighted, and breeding targets could be adjusted to support husbandry, feeding, and management measures to combat SINS and improve animal welfare in pig production.

5. Conclusions

The present study demonstrates successful fine mapping by prior enrichment of certain chromosomal regions before sequencing. The analyzed regions were identified in a previous study using GWAS. This provided evidence for functional SNPs in some candidate genes that may contribute to the pathogenesis of SINS due to their homologous physiological function. The location of some of these SNPs in sweep regions may indicate that unfavorable gene variants for predisposition to SINS may have been favored by selection for performance in pigs. The results of this study can now be used as a starting point for further investigations on a larger number of piglets from a significantly larger number of boars. Ultimately, it is possible to assist in the fight against SINS by improving housing, feeding, and management conditions through practical genetic selection.

Author Contributions

Conceptualization, G.R.; methodology, G.R., H.W., K.G. and M.L.; software, H.W. and K.G.; validation, G.R. and K.G.; formal analysis, G.R., H.W. and K.G.; investigation, J.K., K.G. and S.B; resources, G.R. and S.B.; data curation, G.R., J.M. and K.G.; writing—original draft preparation, K.G. and G.R.; writing—review and editing, G.R.; visualization, G.R. and K.G.; supervision, G.R.; project administration, G.R. and K.G.; funding acquisition, G.R. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Tönnies Research Foundation, Rheda, Germany. The publication costs were supported by the “Open-Access-Publikationsfond” of the Justus Liebig University, Giessen, Germany.

Institutional Review Board Statement

No animal experiments were carried out in the present study. We used just phenotypic data and DNA from animals of a preliminary study (Gerhards et al., 2023 [14]) that was carried out under the approval of the authorities of the Justus Liebig University in Giessen, Germany with file number 708_M.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data are available from the corresponding author by reasonable request.

Acknowledgments

The authors wish to thank the team of the Oberer Hardthof swine facility for excellent technical support. We would also like to thank Berthold Fartmann for the informative exchange on genomic data and bioinformatics. Thanks go to the Tönnies Research Foundation, Rheda-Wiedenbrück, Germany for financial support in the realization of this project.

Conflicts of Interest

The authors declare no conflicts of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

Abbreviations

The following abbreviations are used in this manuscript:
ABCB11ATP-binding cassette, subfamily B member 11
BIVMBasic immunoglobulin-like variable motif-containing protein
CCDC168Coiled-coil domain containing 168
CI9595% confidence interval
DULow SINS Duroc boar
ECDEcdysoneless cell cycle regulator
ERCC5Excision repair 5, endonuclease
ERICH2Glutamate rich 2
fForward (5′–3′)
FAM149B1Member of the family with sequence similarity 149, member B1
GWASGenome-wide association study
HWEHardy–Weinberg equilibrium
IndelsInsertion/deletion polymorphisms
LRP2LDL receptor related protein 2
MAFMinor allele frequency
METTL21CProtein lysine methyltransferase 21C
MissenseMissense variant
MODModifier
MYO3BMyosin IIIB
NNorm (e.g., base, triplet, amino acid (AA)) according to reference genome Sscrofa11.1
NNHomozygous genotype of the norm variant
NUDT13Nudix-hydrolase 13
NVHeterozygous genotype
n.s.Not significant
PSignificance
PCAPrincipal component analysis
PIPietrain boar(s)
PI+Low SINS Pietrain boar
PI−High SINS Pietrain boar
QCQuality control
rReverse (3′–5′)
SEStandard error
SIFTSorting Intolerant From Tolerant Score
SINSSwine inflammation and necrosis syndrome
SNPsSingle nucleotide polymorphisms
SpliceSplicing region
SSCSus scrofa chromosome
UTRUntranslated region
VVariant (e.g., base, triplet, amino acid (AA)) deviating from the reference genome Sscrofa11.1
VEPVariant Effect Predictor (Ensembl)
VVHomozygous genotype of the variant
ZSINSSINS score after Z-transformation

References

  1. Li, Y.; Wu, N.; Xu, R.; Li, L.; Zhou, W.; Zhou, X. Empirical analysis of pig welfare levels and their impact on pig breeding efficiency-Based on 773 pig farmers’ survey data. PLoS ONE 2017, 12, e0190108. [Google Scholar] [CrossRef] [PubMed]
  2. Reiner, G.; Kuehling, J.; Loewenstein, F.; Lechner, M.; Becker, S. Swine Inflammation and Necrosis Syndrome (SINS). Animals 2021, 11, 1670. [Google Scholar] [CrossRef] [PubMed]
  3. Leite, N.G.; Knol, E.F.; Nuphaus, S.; Vogelzang, R.; Tsuruta, S.; Wittmann, M.; Lourenco, D. The genetic basis of swine inflammation and necrosis syndrome and its genetic association with post-weaning skin damage and production traits. J. Anim. Sci. 2023, 101, skad067. [Google Scholar] [CrossRef] [PubMed]
  4. Fortune, H.; Micout, S.; Monjouste, A. Évaluation de la prévalence du syndrome inflammatoire et nécrotique porcin dans les troupeaux français. Journées Rech. Porc. 2024, 56, 301–302. [Google Scholar]
  5. Kuehling, J.; Loewenstein, F.; Wenisch, S.; Kressin, M.; Herden, C.; Lechner, M.; Reiner, G. An in-depth diagnostic exploration of an inflammation and necrosis syndrome in a population of newborn piglets. Animal 2021, 15, 100078. [Google Scholar] [CrossRef]
  6. Kuehling, J.; Eisenhofer, K.; Lechner, M.; Becker, S.; Willems, H.; Reiner, G. The effects of boar on susceptibility to swine inflammation and necrosis syndrome in piglets. Porc. Health Manag. 2021, 7, 15. [Google Scholar] [CrossRef]
  7. Koenders-van Gog, K.; Wijnands, T.; Lechner, M.; Reiner, G.; Fink-Gremmels, J. Screening of Piglets for Signs of Inflammation and Necrosis as Early Life Indicators of Animal Health and Welfare Hazards. Animals 2025, 15, 378. [Google Scholar] [CrossRef]
  8. European Food Safety Authority. Statement on the use of animal-based measures to assess the welfare of animals. EFSA J. 2012, 10, 2767. [Google Scholar] [CrossRef]
  9. European Food Safety Authority. Scientific Opinion concerning a Multifactorial approach on the use of animal and non-animal-based measures to assess the welfare of pigs. EFSA J. 2014, 12, 3702. [Google Scholar] [CrossRef]
  10. Reiner, G.; Kuehling, J.; Lechner, M.; Schrade, H.; Saltzmann, J.; Muelling, C.; Dänicke, S.; Loewenstein, F. Swine inflammation and necrosis syndrome is influenced by husbandry and quality of sow in suckling piglets, weaners and fattening pigs. Porc. Health Manag. 2020, 6, 32. [Google Scholar] [CrossRef]
  11. Ringseis, R.; Gessner, D.K.; Loewenstein, F.; Kuehling, J.; Becker, S.; Willems, H.; Lechner, M.; Eder, K.; Reiner, G. Swine Inflammation and Necrosis Syndrome Is Associated with Plasma Metabolites and Liver Transcriptome in Affected Piglets. Animals 2021, 11, 772. [Google Scholar] [CrossRef]
  12. Gerhards, K.; Becker, S.; Kuehling, J.; Lechner, M.; Willems, H.; Ringseis, R.; Reiner, G. Screening for transcriptomic associations with Swine Inflammation and Necrosis Syndrome. BMC Vet. Res. 2025, 21, 26. [Google Scholar] [CrossRef] [PubMed]
  13. Loewenstein, F.; Becker, S.; Kuehling, J.; Schrade, H.; Lechner, M.; Ringseis, R.; Eder, K.; Moritz, A.; Reiner, G. Inflammation and necrosis syndrome is associated with alterations in blood and metabolism in pigs. BMC Vet. Res. 2022, 18, 50. [Google Scholar] [CrossRef]
  14. Gerhards, K.; Becker, S.; Kuehling, J.; Lechner, M.; Bathke, J.; Willems, H.; Reiner, G. GWAS reveals genomic associations with swine inflammation and necrosis syndrome. Mamm. Genome 2023, 34, 586–601. [Google Scholar] [CrossRef]
  15. Bathke, J.; Lühken, G. OVarFlow: A resource optimized GATK 4 based Open source Variant calling workFlow. BMC Bioinform. 2021, 22, 402. [Google Scholar] [CrossRef]
  16. R Core Team. A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2022; Available online: https://www.R-project.org/ (accessed on 19 February 2025).
  17. RStudio Team. RStudio: Integrated Development Environment for R; RStudio PBC: Boston, MA, USA, 2022; Available online: http://www.rstudio.com/ (accessed on 19 January 2025).
  18. Purcell, S.; Neale, B.; Todd-Brown, K.; Thomas, L.; Ferreira, M.A.R.; Bender, D.; Maller, J.; Sklar, P.; Bakker, P.I.W.; de Daly, M.J.; et al. PLINK: A tool set for whole-genome association and population-based linkage analyses. Am. J. Hum. Genet. 2007, 81, 559–575. [Google Scholar] [CrossRef]
  19. Chang, C.C.; Chow, C.C.; Tellier, L.C.; Vattikuti, S.; Purcell, S.M.; Lee, J.J. Second-generation PLINK: Rising to the challenge of larger and richer datasets. Gigascience 2015, 4, 7. [Google Scholar] [CrossRef]
  20. Bradbury, P.J.; Zhang, Z.; Kroon, D.E.; Casstevens, T.M.; Ramdoss, Y.; Buckler, E.S. TASSEL: Software for association mapping of complex traits in diverse samples. Bioinformatics 2007, 23, 2633–2635. [Google Scholar] [CrossRef]
  21. Wang, J.; Zhang, Z. GAPIT Version 3: Boosting Power and Accuracy for Genomic Association and Prediction. Genom. Proteom. Bioinform. 2021, 19, 629–640. [Google Scholar] [CrossRef]
  22. Huang, M.; Liu, X.; Zhou, Y.; Summers, R.M.; Zhang, Z. BLINK: A package for the next level of genome-wide association studies with both individuals and markers in the millions. Gigascience 2019, 8, giy154. [Google Scholar] [CrossRef]
  23. Cingolani, P.; Platts, A.; Le Wang, L.; Coon, M.; Nguyen, T.; Wang, L.; Land, S.J.; Lu, X.; Ruden, D.M. A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff: SNPs in the genome of Drosophila melanogaster strain w1118; iso-2; iso-3. Fly 2012, 6, 80–92. [Google Scholar] [CrossRef] [PubMed]
  24. McLaren, W.; Gil, L.; Hunt, S.E.; Riat, H.S.; Ritchie, G.R.S.; Thormann, A.; Flicek, P.; Cunningham, F. The Ensembl Variant Effect Predictor. Genome Biol. 2016, 17, 122. [Google Scholar] [CrossRef] [PubMed]
  25. Ng, P.C.; Henikoff, S. SIFT: Predicting amino acid changes that affect protein function. Nucleic Acids Res. 2003, 31, 3812–3814. [Google Scholar] [CrossRef]
  26. Brown, G.R.; Hem, V.; Katz, K.S.; Ovetsky, M.; Wallin, C.; Ermolaeva, O.; Tolstoy, I.; Tatusova, T.; Pruitt, K.D.; Maglott, D.R.; et al. Gene: A gene-centered information resource at NCBI. Nucleic Acids Res. 2015, 43, D36–D42. [Google Scholar] [CrossRef]
  27. Safran, M.; Rosen, N.; Twik, M.; BarShir, R.; Iny Stein, T.; Dahary, D.; Fishilevich, S.; Lancet, D. The GeneCards Suite. In Practical Guide to Life Science Databases; Springer: Berlin/Heidelberg, Germany, 2022; pp. 27–56. [Google Scholar]
  28. Belinky, F.; Nativ, N.; Stelzer, G.; Zimmerman, S.; Iny Stein, T.; Safran, M.; Lancet, D. PathCards: Multi-source consolidation of human biological pathways. Database 2015, 2015, bav006. [Google Scholar] [CrossRef]
  29. Kanehisa, M.; Goto, S. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000, 28, 27–30. [Google Scholar] [CrossRef]
  30. Baldarelli, R.M.; Smith, C.L.; Ringwald, M.; Richardson, J.E.; Bult, C.J. Mouse Genome Informatics: An integrated knowledgebase system for the laboratory mouse. Genetics 2024, 227, iyae031. [Google Scholar] [CrossRef]
  31. Sherman, B.T.; Hao, M.; Qiu, J.; Jiao, X.; Baseler, M.W.; Lane, H.C.; Imamichi, T.; Chang, W. DAVID: A web server for functional enrichment analysis and functional annotation of gene lists (2021 update). Nucleic Acids Res. 2022, 50, W216–W221. [Google Scholar] [CrossRef]
  32. Rauw, W.; Kanis, E.; Noordhuizen-Stassen, E.; Grommers, F. Undesirable side effects of selection for high production efficiency in farm animals: A review. Livest. Prod. Sci. 1998, 56, 15–33. [Google Scholar] [CrossRef]
  33. Rydhmer, L.; Lundeheim, N. Chapter 9. Breeding Pigs for Improved Welfare. Welfare of Pigs from Birth to Slaughter; Wageningen Akad. Publ; Éditions Qæ: Wageningen, The Netherlands; Versailles, France, 2008; pp. 243–270. ISBN 9789086866373. [Google Scholar]
  34. Tan, C.; Wu, Z.; Ren, J.; Huang, Z.; Liu, D.; He, X.; Prakapenka, D.; Zhang, R.; Li, N.; Da, Y.; et al. Genome-wide association study and accuracy of genomic prediction for teat number in Duroc pigs using genotyping-by-sequencing. Genet. Sel. Evol. 2017, 49, 35. [Google Scholar] [CrossRef] [PubMed]
  35. Shin, D.; Won, K.-H.; Song, K.-D. In silico approaches to discover the functional impact of non-synonymous single nucleotide polymorphisms in selective sweep regions of the Landrace genome. Asian-Australas. J. Anim. Sci. 2018, 31, 1980–1990. [Google Scholar] [CrossRef] [PubMed]
  36. Bertani, G.R.; Gladney, C.D.; Johnson, R.K.; Pomp, D. Evaluation of gene expression in pigs selected for enhanced reproduction using differential display PCR: II. Anterior pituitary. J. Anim. Sci. 2004, 82, 32–40. [Google Scholar] [CrossRef] [PubMed]
  37. Yoo, C.-K.; Cho, I.-C.; Lee, J.-B.; Jung, E.-J.; Lim, H.-T.; Han, S.-H.; Lee, S.-S.; Ko, M.-S.; Kang, T.; Hwang, J.-H.; et al. QTL analysis of clinical-chemical traits in an F₂ intercross between Landrace and Korean native pigs. Physiol. Genom. 2012, 44, 657–668. [Google Scholar] [CrossRef]
  38. Wiederstein, J.L.; Nolte, H.; Günther, S.; Piller, T.; Baraldo, M.; Kostin, S.; Bloch, W.; Schindler, N.; Sandri, M.; Blaauw, B.; et al. Skeletal Muscle-Specific Methyltransferase METTL21C Trimethylates p97 and Regulates Autophagy-Associated Protein Breakdown. Cell Rep. 2018, 23, 1342–1356. [Google Scholar] [CrossRef]
  39. Huang, J.; Hsu, Y.-H.; Mo, C.; Abreu, E.; Kiel, D.P.; Bonewald, L.F.; Brotto, M.; Karasik, D. METTL21C is a potential pleiotropic gene for osteoporosis and sarcopenia acting through the modulation of the NF-κB signaling pathway. J. Bone Miner. Res. 2014, 29, 1531–1540. [Google Scholar] [CrossRef]
  40. Mayer, M.P.; Bukau, B. Hsp70 chaperones: Cellular functions and molecular mechanism. Cell. Mol. Life Sci. 2005, 62, 670–684. [Google Scholar] [CrossRef]
  41. Cheng, J.; Li, R.; Wang, L.; Zhang, T.; Yuan, G.; Lu, H. Methyltransferase like 21C interaction with Hsc70 affects chicken myoblast differentiation. Ital. J. Anim. Sci. 2023, 22, 605–614. [Google Scholar] [CrossRef]
  42. He, J.; Xia, C.; He, Y.; Pan, D.; Cao, J.; Sun, Y.; Zeng, X. Proteomic responses to oxidative damage in meat from ducks exposed to heat stress. Food Chem. 2019, 295, 129–137. [Google Scholar] [CrossRef]
  43. Yang, G.; Lu, H.; Wang, L.; Zhao, J.; Zeng, W.; Zhang, T. Genome-Wide Identification and Transcriptional Expression of the METTL21C Gene Family in Chicken. Genes 2019, 10, 628. [Google Scholar] [CrossRef]
  44. Wang, S.; Zhao, J.; Wang, L.; Zhang, T.; Zeng, W.; Lu, H. METTL21C mediates lysine trimethylation of IGF2BP1 to regulate chicken myoblast proliferation. Br. Poult. Sci. 2023, 64, 74–80. [Google Scholar] [CrossRef]
  45. Steinert, N.D.; Jorgenson, K.W.; Lin, K.-H.; Hermanson, J.B.; Lemens, J.L.; Hornberger, T.A. A novel method for visualizing in-vivo rates of protein degradation provides insight into how TRIM28 regulates muscle size. iScience 2023, 26, 106526. [Google Scholar] [CrossRef] [PubMed]
  46. Ma, W.; Lu, Y.; Wang, C. Production performance, egg quality, and uterine gene expression for layers as affected by N-Carbamylglutamate supplementation. Front. Vet. Sci. 2023, 10, 1110801. [Google Scholar] [CrossRef]
  47. Fedorova, E.S.; Dementieva, N.V.; Shcherbakov, Y.S.; Stanishevskaya, O.I. Identification of Key Candidate Genes in Runs of Homozygosity of the Genome of Two Chicken Breeds, Associated with Cold Adaptation. Biology 2022, 11, 547. [Google Scholar] [CrossRef]
  48. Filippin, L.I.; Vercelino, R.; Marroni, N.P.; Xavier, R.M. Redox signalling and the inflammatory response in rheumatoid arthritis. Clin. Exp. Immunol. 2008, 152, 415–422. [Google Scholar] [CrossRef]
  49. Lambeth, J.D. NOX enzymes and the biology of reactive oxygen. Nat. Rev. Immunol. 2004, 4, 181–189. [Google Scholar] [CrossRef]
  50. Bulua, A.C.; Simon, A.; Maddipati, R.; Pelletier, M.; Park, H.; Kim, K.-Y.; Sack, M.N.; Kastner, D.L.; Siegel, R.M. Mitochondrial reactive oxygen species promote production of proinflammatory cytokines and are elevated in TNFR1-associated periodic syndrome (TRAPS). J. Exp. Med. 2011, 208, 519–533. [Google Scholar] [CrossRef]
  51. Reuter, S.; Gupta, S.C.; Chaturvedi, M.M.; Aggarwal, B.B. Oxidative stress, inflammation, and cancer: How are they linked? Free. Radic. Biol. Med. 2010, 49, 1603–1616. [Google Scholar] [CrossRef]
  52. Auf dem Keller, U.; Kümin, A.; Braun, S.; Werner, S. Reactive oxygen species and their detoxification in healing skin wounds. J. Investig. Dermatol. Symp. Proc. 2006, 11, 106–111. [Google Scholar] [CrossRef]
  53. Waris, G.; Ahsan, H. Reactive oxygen species: Role in the development of cancer and various chronic conditions. J. Carcinog. 2006, 5, 14. [Google Scholar] [CrossRef]
  54. Zhao, W.; Zhao, T.; Chen, Y.; Ahokas, R.A.; Sun, Y. Reactive oxygen species promote angiogenesis in the infarcted rat heart. Int. J. Exp. Pathol. 2009, 90, 621–629. [Google Scholar] [CrossRef]
  55. Colavitti, R.; Pani, G.; Bedogni, B.; Anzevino, R.; Borrello, S.; Waltenberger, J.; Galeotti, T. Reactive oxygen species as downstream mediators of angiogenic signaling by vascular endothelial growth factor receptor-2/KDR. J. Biol. Chem. 2002, 277, 3101–3108. [Google Scholar] [CrossRef] [PubMed]
  56. Kim, Y.-W.; West, X.Z.; Byzova, T.V. Inflammation and oxidative stress in angiogenesis and vascular disease. J. Mol. Med. 2013, 91, 323–328. [Google Scholar] [CrossRef] [PubMed]
  57. Jerez-Timaure, N.; Gallo, C.; Ramírez-Reveco, A.; Greif, G.; Strobel, P.; Pedro, A.V.F.; Morera, F.J. Early differential gene expression in beef Longissimus thoracis muscles from carcasses with normal (<5.8) and high (5.9) ultimate pH. Meat Sci. 2019, 153, 117–125. [Google Scholar] [CrossRef]
  58. Stefl, S.; Nishi, H.; Petukh, M.; Panchenko, A.R.; Alexov, E. Molecular mechanisms of disease-causing missense mutations. J. Mol. Biol. 2013, 425, 3919–3936. [Google Scholar] [CrossRef]
  59. Zhang, Z.; Miteva, M.A.; Wang, L.; Alexov, E. Analyzing effects of naturally occurring missense mutations. Comput. Math. Methods Med. 2012, 2012, 805827. [Google Scholar] [CrossRef]
  60. Chatterjee, S.; Pal, J.K. Role of 5′- and 3′-untranslated regions of mRNAs in human diseases. Biol. Cell 2009, 101, 251–262. [Google Scholar] [CrossRef]
  61. Huang, Y.; Yuan, C.; Zhao, Y.; Li, C.; Cao, M.; Li, H.; Zhao, Z.; Sun, A.; Basang, W.; Zhu, Y.; et al. Identification and Regulatory Network Analysis of Genes Related to Reproductive Performance in the Hypothalamus and Pituitary of Angus Cattle. Genes 2022, 13, 965. [Google Scholar] [CrossRef]
  62. Hernandez, A.S.; Zayas, G.A.; Rodriguez, E.E.; Sarlo Davila, K.M.; Rafiq, F.; Nunez, A.N.; Titto, C.G.; Mateescu, R.G. Exploring the genetic control of sweat gland characteristics in beef cattle for enhanced heat tolerance. J. Anim. Sci. Biotechnol. 2024, 15, 66. [Google Scholar] [CrossRef]
  63. ElSharawy, A.; Hundrieser, B.; Brosch, M.; Wittig, M.; Huse, K.; Platzer, M.; Becker, A.; Simon, M.; Rosenstiel, P.; Schreiber, S.; et al. Systematic evaluation of the effect of common SNPs on pre-mRNA splicing. Hum. Mutat. 2009, 30, 625–632. [Google Scholar] [CrossRef]
  64. Bai, Z.; Luo, Y.; Tian, L. ERCC5, HES6 and RORA are potential diagnostic markers of coronary artery disease. FEBS Open Bio 2022, 12, 1814–1827. [Google Scholar] [CrossRef]
  65. Zuo, C.; Lv, X.; Liu, T.; Yang, L.; Yang, Z.; Yu, C.; Chen, H. Polymorphisms in ERCC4 and ERCC5 and risk of cancers: Systematic research synopsis, meta-analysis, and epidemiological evidence. Front. Oncol. 2022, 12, 951193. [Google Scholar] [CrossRef] [PubMed]
  66. Maroilley, T.; Lemonnier, G.; Lecardonnel, J.; Esquerré, D.; Ramayo-Caldas, Y.; Mercat, M.J.; Rogel-Gaillard, C.; Estellé, J. Deciphering the genetic regulation of peripheral blood transcriptome in pigs through expression genome-wide association study and allele-specific expression analysis. BMC Genom. 2017, 18, 967. [Google Scholar] [CrossRef] [PubMed]
  67. Shaheen, R.; Jiang, N.; Alzahrani, F.; Ewida, N.; Al-Sheddi, T.; Alobeid, E.; Musaev, D.; Stanley, V.; Hashem, M.; Ibrahim, N.; et al. Bi-allelic Mutations in FAM149B1 Cause Abnormal Primary Cilium and a Range of Ciliopathy Phenotypes in Humans. Am. J. Hum. Genet. 2019, 104, 731–737. [Google Scholar] [CrossRef] [PubMed]
  68. Kim, H.-J.; Son, H.; Sung, J.; Yun, J.M.; Kwon, H.; Cho, B.; Kim, J.-I.; Park, J.-H. A Genome-Wide Association Study on Abdominal Adiposity-Related Traits in Adult Korean Men. Obes. Facts 2022, 15, 590–599. [Google Scholar] [CrossRef]
  69. Qiu, Y.; Zhuang, Z.; Meng, F.; Ruan, D.; Xu, C.; Ma, F.; Peng, L.; Ding, R.; Cai, G.; Yang, M.; et al. Identification of candidate genes associated with carcass component weights in commercial crossbred pigs through a combined GWAS approach. J. Anim. Sci. 2023, 101, skad121. [Google Scholar] [CrossRef]
  70. Saini, T.; Chauhan, A.; Ahmad, S.F.; Kumar, A.; Vaishnav, S.; Singh, S.; Mehrotra, A.; Bhushan, B.; Gaur, G.K.; Dutt, T. Elucidation of population stratifying markers and selective sweeps in crossbred Landlly pig population using genome-wide SNP data. Mamm. Genome 2024, 35, 170–185. [Google Scholar] [CrossRef]
  71. Wang, Q.; Liu, Z.; Zeng, X.; Zheng, Y.; Lan, L.; Wang, X.; Lai, Z.; Hou, X.; Gao, L.; Liang, L.; et al. Integrated analysis of miRNA-mRNA expression of newly emerging swine H3N2 influenza virus cross-species infection with tree shrews. Virol. J. 2024, 21, 4. [Google Scholar] [CrossRef]
  72. Kanis, E.; van den Belt, H.; Groen, A.F.; Schakel, J.; Greef, K.H.d. Breeding for improved welfare in pigs: A conceptual framework and its use in practice. Anim. Sci. 2004, 78, 315–329. [Google Scholar] [CrossRef]
  73. Guy, S.Z.Y.; Thomson, P.C.; Hermesch, S. Selection of pigs for improved coping with health and environmental challenges: Breeding for resistance or tolerance? Front. Genet. 2012, 3, 31610. [Google Scholar] [CrossRef]
Figure 1. SINS scores in offspring from the three boars (DU = low SINS Duroc boar; PI+ = low SINS Pietrain boar, PI− = high SINS Pietrain boar). Number in box: mean; letter in box: groups with different letters were significantly different at p < 0.001; black box: lower 95% confidence interval (CI95) to mean; grey box: mean to upper 95% confidence interval.
Figure 1. SINS scores in offspring from the three boars (DU = low SINS Duroc boar; PI+ = low SINS Pietrain boar, PI− = high SINS Pietrain boar). Number in box: mean; letter in box: groups with different letters were significantly different at p < 0.001; black box: lower 95% confidence interval (CI95) to mean; grey box: mean to upper 95% confidence interval.
Vetsci 12 00508 g001
Figure 2. SINS scores by boar/piglet genotype combinations. (A): Series 1; (B): Series 2. DU = Low SINS Duroc boar; PI+ = Low SINS Pietrain boar, PI− = High SINS Pietrain boar; N: Norm allele as in the reference genome Sscrofa11.1; V: variant allele. Statistically significant differences are observed between columns with different letters (p < 0.05).
Figure 2. SINS scores by boar/piglet genotype combinations. (A): Series 1; (B): Series 2. DU = Low SINS Duroc boar; PI+ = Low SINS Pietrain boar, PI− = High SINS Pietrain boar; N: Norm allele as in the reference genome Sscrofa11.1; V: variant allele. Statistically significant differences are observed between columns with different letters (p < 0.05).
Vetsci 12 00508 g002
Figure 3. Effects of boar and SNP series on the SINS score. Boar: A = Low SINS Duroc; B = PI+ = Low SINS Pietrain; C = PI− = High SINS Pietrain; SNP series: A = NN; B = NV; C = VV; N: Allele as in the reference genome Sscrofa11.1; V: Variant allele. Statistically significant differences are observed between columns with different letters (p < 0.05).
Figure 3. Effects of boar and SNP series on the SINS score. Boar: A = Low SINS Duroc; B = PI+ = Low SINS Pietrain; C = PI− = High SINS Pietrain; SNP series: A = NN; B = NV; C = VV; N: Allele as in the reference genome Sscrofa11.1; V: Variant allele. Statistically significant differences are observed between columns with different letters (p < 0.05).
Vetsci 12 00508 g003
Table 1. Overview of the collected individual characteristics and scores.
Table 1. Overview of the collected individual characteristics and scores.
Skin ShinyNo BristlesSwellingRed-
ness
Scab
Formation
RhagadesExudationNecrosisBleedingRing-
Shaped
Constrictions
Vein
Congestion
Edema Additive
Body Part
Scores
Additive
SINS-
Score
Tail basen.s.0/1 30/10/1n.s.0/10/10/1n.s.n.s.n.s.n.s.0–6Vetsci 12 00508 i001
Tail tipn.s.0/10/10/10/10/10/10/10/10/1n.s.n.s.0–9
Ears0/10/1n.s.n.s.n.s.n.s.n.s.0/1n.s.n.s.0/1n.s.0–4
Face 1n.s.n.s.n.s.n.s.n.s.n.s.n.s.n.s.n.s.n.s.n.s.0/10–2
Teatsn.s.n.s.0/10/10/1n.s.n.s.0/1n.s.n.s.0/1n.s.0–5
Naveln.s.n.s.0/10/1n.s.n.s.n.s.n.s.n.s.n.s.n.s.n.s.0–20–44
Coronary bands 2n.s.n.s.0/10/1n.s.n.s.n.s.n.s.n.s.n.s.n.s.n.s.0–2
Claw wall 2n.s.n.s.0/1n.s.n.s.n.s.n.s.n.s.0/1n.s.n.s.n.s.0–2
Heels 2n.s.n.s.0–3n.s.n.s.n.s.n.s.n.s.0–3n.s.n.s.n.s.0–12
n.s.: Not scored; 1: Eye lids and nose back are separately scored; 2: Hind and front limbs are separately scored; 3: Present (1) or absent (0).
Table 2. Chromosomes and regions selected for the present study based on the significance level of the major SNPs obtained by Gerhards et al. [14].
Table 2. Chromosomes and regions selected for the present study based on the significance level of the major SNPs obtained by Gerhards et al. [14].
SSC 1Start (Mbp)End (Mbp)Length (Mbp)SNP
Name
SNP
Position
SNP Effect
-log(p)
723.526.12.6rs3385099482602536116.3
989.7490.741rs341512035902415778.1
1170.971.91rs3475903338712396798.9
1243.644.40.8rs1112423847447384238.4
1475.576.30.8rs3234888367591819812.7
1491.492.20.8rs3426125619180893413.4
1575772rs3392705827692610616.8
1644451rs3312174554466935811.3
1739.6540.651rs3416286114015712814.8
1 SSC: Sus scrofa chromosome; p: Significance.
Table 3. Mean, standard error, lower and upper 95% confidence intervals, and prevalence in offspring of the three boars (DU; PI+; PI−).
Table 3. Mean, standard error, lower and upper 95% confidence intervals, and prevalence in offspring of the three boars (DU; PI+; PI−).
ScoreBoarMeanSELower CI95Upper CI95Prevalence (%)p
Tail baseDU0.77a 10.160.451.0948.9<0.001
PI+0.92a0.130.671.1752
PI−1.45b0.111.241.6668.5
Tail tipDU0.43a0.150.130.7329.80.012
PI+0.57a0.120.340.8145.3
PI−0.92b0.100.721.1145.4
EarsDU1.57a0.151.281.8778.7<0.001
PI+2.59b0.122.352.8297.3
PI−2.56b0.102.362.7597.2
TeatsDU0.74a0.320.111.3823.40.002
PI+1.59b0.261.082.0965.3
PI−2.10b0.211.682.5260.2
Coronary bandsDU1.04a0.130.781.3161.7<0.001
PI+0.89a0.110.681.1061.3
PI−1.40b0.091.221.5783.3
Claw wallDU1.600.091.421.7883n.s.
PI+1.810.071.671.9696
PI−1.720.061.601.8489.8
HeelsDU5.91a0.285.366.4793.6(0.055)
PI+5.75ab0.225.316.18100
PI−6.42b0.196.056.7899.1
SE: Standard error; p: Significance; DU = Low SINS Duroc boar, PI+ = Low SINS Pietrain boar, PI− = High SINS Pietrain boar; 1: Values within one phenotype with different letters are statistically significant at p < 0.05.
Table 4. Associated functional SNPs by series, position, and phenotype with significance, base exchange, and meaning.
Table 4. Associated functional SNPs by series, position, and phenotype with significance, base exchange, and meaning.
Series 1SSC 2PositionPhenotype 1 3Phenotype 2 3BaseGeneTripletAADirection 6Ex-onEffect 7Meaning (SNPEff)
trait-LOG(p)trait-LOG(p)N 4V 5 N 4V 5N 4V 5
11171015446Heels7.5Tail tip19.2GAMETTL21CCGTCATTMr1MissenseMOD
11171020888Heels7.5Tail tip19.2TCMETTL21C 15′UTRLOW
11171054195Heels7.5Tail tip19.2TCCCDC168CTGCCGERr1MissenseMOD
11171054252Heels7.5Tail tip19.2TCCCDC168ATTACTNSr1MissenseMOD
11171056181Heels7.5Tail tip19.2CTCCDC168ACGATGRHr1MissenseMOD
11171056239Heels7.5Tail tip19.2TCCCDC168CCTCCCRGr1MissenseMOD
11171056968Heels7.5Tail tip19.2CGCCDC168CTCCTGQEr1MissenseMOD
11171057564Heels7.5Tail tip19.2TCCCDC168CTGCCGERr1MissenseMOD
11171114837Heels7.5Tail tip19.2AGBIVM 2SpliceLOW
11171167929Heels7.5Tail tip19.2GAERCC5CGTCATRHf15MissenseMOD
11475945482Heels7.5Tail tip19.2GCNUDT13AGAACARTf2MissenseMOD
11475953186Heels7.5Tail tip19.2AGNUDT13ATGGTGMVf4MissenseMOD
11475956113Heels7.5Tail tip19.2GCNUDT13CGACCARpf7MissenseMOD
11475977767Heels7.5Tail tip19.2CTECD 5SpliceLOW
11475998607Heels7.5Tail tip19.2CTFAM149B1ACAATATIf3MissenseMOD
11476053610Heels7.5Tail tip19.2AGFAM149B1AAGAGGKRf11MissenseMOD
11575617207Tail tip19.2Heels7.5TCLRP2GATGACIVr33MissenseMOD
11575626516Tail tip19.2Heels7.5ACLRP2 41SpliceLOW
11576939828Tail tip19.2Heels7.5GTERICH2promotor 15′UTRLOW
21575405294Tail base8.4 GTABCB11 2SpliceLOW
21576734995Tail base8.4 AGMYO3BAGTGGTSGf30MissenseMOD
21576646899Tail base11.3 AGMYO3B 22SpliceLOW
1: All SNPs were linked in two series; 2: Sus scrofa chromosome; 3: Up to two different phenotypes were significantly associated in the GWAS (see -log(p)); 4: N = Norm relative to the reference genome Sscrofa11.1; 5: V = Variant deviating from the reference genome Sscrofa11.1; 6: direction: f = 5′ to 3′; r = 3′ to 5′; 7: Splice = Splicing region variant, 5′UTR = 5′ Untranslated region variant, Missense = Missense variant.
Table 5. Sow and boar genotypes derived from piglet genotypes.
Table 5. Sow and boar genotypes derived from piglet genotypes.
Boar Sow Piglets
SeriesBreedGenotypeGenotypenGenotypen
1DUNNVV12NV42
NN1NN6
PI+VVVV14VV77
PI−VVVV21VV96
NV5NV13
2DUNNVV12NV42
NN1NN6
PI+VVVV14VV77
PI−NVVV26VV59
NV50
Table 6. SINS scores of the individual body parts by boar/piglet genotype combination and series.
Table 6. SINS scores of the individual body parts by boar/piglet genotype combination and series.
DU/NNDU/NVPI−/NVPI+/VVPI−/VVp
Series 1
Tail base0.15 ± 0.37 a0.73 ± 0.15 a1.11 ± 0.27 b0.88 ± 0.12 a1.5 ± 0.11 b<0.001
Tail tip0.44 ± 0.410.56 ± 0.170.35 ± 0.300.65 ± 0.130.77 ± 0.12n.s.
Ears−0.13 ± 0.30 a1.7 ± 0.13 b1.2 ± 0.22 c2.51 ± 0.1 d2.80 ± 0.09 d<0.001
Teats0.1 ± 0.72 a1.05 ± 0.3 ab1.1 ± 0.53 ab1.62 ± 0.23 b2.32 ± 0.21 c<0.001
Coronary bands−0.11 ± 0.34 a1.15 ± 0.14 bc1.4 ± 0.25 bc0.94 ± 0.11 b1.36 ± 0.10 c<0.001
Claw wall0.15 ± 0.17 a1.71 ± 0.07 b1.86 ± 0.13 b1.82 ± 0.05 b1.71 ± 0.05 b<0.001
Heels1.93 ± 0.7 a6.47 ± 0.29 c4.94 ± 0.51 b5.61 ± 0.22 b6.44 ± 0.20 c<0.001
Series 2
DU/NNDU/NVPI+/VVPI−/NVPI−/VVp
Tail base0.1 ± 0.4 a0.7 ± 0.2 a0.9 ± 0.1 a1.8 ± 0.1 b1.3 ± 0.1 c<0.001
Tail tip0.5 ± 0.40.6 ± 0.20.6 ± 0.10.9 ± 0.10.7 ± 0.1n.s.
Ears−0.1 ± 0.3 a1.7 ± 0.1 b2.5 ± 0.1 c2.7 ± 0.1 c2.7 ± 0.1 c<0.001
Teats0.1 ± 0.8 a1 ± 0.3 ac1.6 ± 0.2 bc2.1 ± 0.3 b2.1 ± 0.3 b0.018
Coronary bands0 ± 0.3 a1.2 ± 0.1 bc0.9 ± 0.1 b1.3 ± 0.1 c1.4 ± 0.1 c<0.001
Claw wall0.1 ± 0.2 a1.7 ± 0.1 b1.8 ± 0.1 b1.7 ± 0.1 b1.8 ± 0.1 b<0.001
Heels2.1 ± 0.7 a6.5 ± 0.3 b5.6 ± 0.2 c6.1 ± 0.3 b6.7 ± 0.3 b<0.001
DU = Low SINS Duroc boar; PI+ = Low SINS Pietrain boar, PI− = High SINS Pietrain boar; N: Allele as in the reference genome Sscrofa11.1, V: Variant allele; p: Significance. Statistically significant differences are observed between values in rows with different letters (p < 0.05).
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Gerhards, K.; Becker, S.; Kühling, J.; Mickan, J.; Lechner, M.; Willems, H.; Reiner, G. Fine Mapping Identifies Candidate Genes Associated with Swine Inflammation and Necrosis Syndrome. Vet. Sci. 2025, 12, 508. https://doi.org/10.3390/vetsci12050508

AMA Style

Gerhards K, Becker S, Kühling J, Mickan J, Lechner M, Willems H, Reiner G. Fine Mapping Identifies Candidate Genes Associated with Swine Inflammation and Necrosis Syndrome. Veterinary Sciences. 2025; 12(5):508. https://doi.org/10.3390/vetsci12050508

Chicago/Turabian Style

Gerhards, Katharina, Sabrina Becker, Josef Kühling, Joel Mickan, Mirjam Lechner, Hermann Willems, and Gerald Reiner. 2025. "Fine Mapping Identifies Candidate Genes Associated with Swine Inflammation and Necrosis Syndrome" Veterinary Sciences 12, no. 5: 508. https://doi.org/10.3390/vetsci12050508

APA Style

Gerhards, K., Becker, S., Kühling, J., Mickan, J., Lechner, M., Willems, H., & Reiner, G. (2025). Fine Mapping Identifies Candidate Genes Associated with Swine Inflammation and Necrosis Syndrome. Veterinary Sciences, 12(5), 508. https://doi.org/10.3390/vetsci12050508

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop