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

Discover QTLs for the level of blood components in Shaoxing duck using GWAS and haplotype sharing analysis

Blood composition is indicative of health-related traits such as immunity and metabolism. The use of molecular genetics to investigate alterations in these attributes in laying ducks is a novel approach. Our objective was to employ genome - wide association studies (GWAS) and haplotype - sharing analysis to identify genomic regions and potential genes associated with 11 blood components in Shaoxing ducks. Our findings revealed 35 SNPs and 1 SNP associated with low-density lipoprotein cholesterol (LDL) and globulin (GLB), respectively. We identified 36 putative candidate genes for the LDL trait in close proximity to major QTLs and key loci. Based on their biochemical and physiological properties, TRA2A, NPY, ARHGEF26, DHX36, and AADAC are the strongest putative candidate genes. Through linkage disequilibrium analysis and haplotype sharing analysis, we identified three haplotypes and one haplotype, respectively, that were significantly linked with LDL and GLB. These haplotypes could be selected as potential candidate haplotypes for molecular breeding of Shaoxing ducks. Additionally, we utilized a bootstrap test to verify the reliability of GWAS with small experimental samples. The test can be accessed at https://github.com/xuwenwu24/Bootstrap-test.

Unravelling the Genetic Architecture of Serum Biochemical Indicators in Sheep

Serum biochemical indicators serve as vital proxies that reflect the physiological state and functions of different organs. The genetic parameters and molecular mechanisms underlying serum biochemical indicators of sheep (Ovis aries) have not been well understood. Therefore, the aim of the present study was to identify the genetic architecture and genomic loci underlying ten serum biochemical indicators in sheep, including alanine transaminase, aspartate transferase, lactate dehydrogenase, cholesterol, glucose, phosphorus, calcium, creatinine, urea and total protein levels. We implemented genetic parameter estimations and GWASs for each trait in 422 Akkaraman lambs. Overall, low to moderate heritability estimates were found in the range of 0.14-0.55. Additionally, low to high genetic correlations were observed among traits. In total, 23 SNP loci were associated with serum biochemical indicators leading to 19 genes. These were SPTA1, MGST2, CACUL1, IGFBP7, PARD3, PHB1, SLC15A5, TRIM35, RGS6, NUP93, CNTNAP2, SLC7A11, B3GALT5, DPP10, HST2ST1, NRP1, LRP1B, MAP3K9 and ENSOARG00020040484.1, as well as LOC101103187, LOC101117162, LOC105611309 and LOC101118029. To our knowledge, these data provide the first associations between SPTA1 and serum cholesterol and between ENSOARG00020040484.1 and serum glucose. The current findings provide a comprehensive inventory of the relationships between serum biochemical parameters, genetic variants and disease-relevant characteristics. This information may facilitate the identification of therapeutic targets and fluid biomarkers and establish a strong framework for comprehending the pathobiology of complex diseases as well as providing targets for sheep genetic improvement programs.

Copy number variants selected during pig domestication inferred from whole genome resequencing

Over extended periods of natural and artificial selection, China has developed numerous exceptional pig breeds. Deciphering the germplasm characteristics of these breeds is crucial for their preservation and utilization. While many studies have employed single nucleotide polymorphism (SNP) analysis to investigate the local pig germplasm characteristics, copy number variation (CNV), another significant type of genetic variation, has been less explored in understanding pig resources. In this study, we examined the CNVs of 18 Wanbei pigs (WBP) using whole genome resequencing data with an average depth of 12.61. We identified a total of 8,783 CNVs (~30.07 Mb, 1.20% of the pig genome) in WBP, including 8,427 deletions and 356 duplications. Utilizing fixation index (Fst), we determined that 164 CNVs were within the top 1% of the Fst value and defined as under selection. Functional enrichment analyses of the genes associated with these selected CNVs revealed genes linked to reproduction (SPATA6, CFAP43, CFTR, BPTF), growth and development (NR6A1, SMYD3, VIPR2), and immunity (PARD3, FYB2). This study enhances our understanding of the genomic characteristics of the Wanbei pig and offers a theoretical foundation for the future breeding of this breed.

Genome-wide association study identifies multiple loci influencing duck serum biochemical indicators in the laying period

1. Laying performance is an important economic trait in poultry. The blood is essential in transporting nutrients to the yolk and albumen and is necessary for egg formation.2. This study calculated the phenotypic relationships of duck egg quality, egg production efficiency and 22 serum parameters in the egg-laying stage. Using a variety of methodologies, a genome-wide association study (GWAS) was carried out to uncover the genetic foundations of the 22 serum biochemical markers of laying ducks.3. Spearman correlation coefficients between the egg production (226-329 per day) and the serum parameters were all weak, being less than 0.3. This analysis was done on 22 serum parameters, with total protein (TP), total triglycerides (TG), calcium (Ca) and phosphorous (P) having the highest correlation coefficients (r = 0.56-0.88). The coefficients for blood markers, such as total cholesterol (CHOL), total bilirubin (TBIL), low-density lipoprotein cholesterol (LDL-C) and high-density lipoprotein cholesterol (HDL-C) varied from 0.70-0.94.4. Based on single-marker single-trait genome-wide analyses by a mixed linear model program of EMMAX, nine candidate genes were associated with enzyme traits (AST/ALT aspartate transaminase/glutamic-pyruvic transaminase, creatine kinase) and 19 candidate genes were associated with metabolism and protein-related serum parameters (glucose, total bile acid, uric acid (UA), albumin (ALB).5. The mvLMM (multivariate linear mixed model) of GEMMA software was used to carry out multiple trait integrated GWAS. Two candidate genes were found in the TP-TG-CA-P analysis and seven candidate genes in the CHOL_LDL-C_HDL-C_TBIL study. There was a high genetic correlation between the two groups.

Investigation of ABO Gene Variants across More Than 60 Pig Breeds and Populations and Other Suidae Species Using Whole-Genome Sequencing Datasets

Polymorphisms in the human ABO gene determine the major blood classification system based on the three well-known forms: A; B; and O. In pigs that carry only two main alleles in this gene (A and O), we still need to obtain a more comprehensive distribution of variants, which could also impact its function. In this study, we mined more than 500 whole-genome sequencing datasets to obtain information on the ABO gene in different Suidae species, pig breeds, and populations and provide (i) a comprehensive distribution of the A and O alleles, (ii) evolutionary relationships of ABO gene sequences across Suidae species, and (iii) an exploratory evaluation of the effect of the different ABO gene variants on production traits and blood-related parameters in Italian Large White pigs. We confirmed that allele O is likely under balancing selection, present in all Sus species investigated, without being fixed in any of them. We reported a novel structural variant in perfect linkage disequilibrium with allele O that made it possible to estimate the evolutionary time window of occurrence of this functional allele. We also identified two single nucleotide polymorphisms that were suggestively associated with plasma magnesium levels in pigs. Other studies can also be constructed over our results to further evaluate the effect of this gene on economically relevant traits and basic biological functions.

Genome-wide association study using a single-step approach for teat number in Duroc, Landrace and Yorkshire pigs in Korea

We investigated the genetic basis of teat number in sows, which is an important factor in their reproductive performance. We collected genotyping data from 20 353 pigs of three breeds (Duroc, Landrace and Yorkshire) using the Porcine SNP60K Bead Chip, and analyzed phenotypic data from 240 603 pigs. The heritability values of total teat number were 0.33 ± 0.02, 0.51 ± 0.01 and 0.50 ± 0.01 in Duroc, Landrace and Yorkshire pigs, respectively. A genome-wide association study was used to identify significant chromosomal regions associated with teat number in SSC7 and SSC9 in Duroc pig, SSC3, SSC7 and SSC18 in Landrace pig, and SSC7, SSC8 and SSC10 in Yorkshire pig. Among the markers, MARC0038565, located between the vertnin (VRTN) and synapse differentiation-inducing 1-like (SYNDIG1L) genes, showed the strongest association in the Duroc pig and was significant in all breeds. In Landrace and Yorkshire pigs, the most significant markers were located within the apoptosis resistant E3 ubiquitin protein ligase 1 (AREL1) and latent transforming growth factor beta-binding protein 2 (LTBP2) genes in SSC7, respectively. VRTN is a candidate gene regulating the teat number. Most markers were located in SSC7, indicating their significance in determining teat number and their potential as valuable genomic selection targets for improving this trait. Extensive linkage disequilibrium blocks were identified in SSC7, supporting their use in genomic selection strategies. Our study provides valuable insights into the genetic architecture of teat numbers in pigs, and helps identify candidate genes and genomic regions that may contribute to this economically important trait.

Genome-wide population structure, homozygosity, and heterozygosity patterns of Nero Siciliano pig in the framework of Italian and cosmopolitan breeds

Analysis of genomic data is becoming more and more common for the effective management of livestock breeding programmes, even in the case of local populations. In this work, the genome-wide data of Nero Siciliano pig breed were compared to that of wild boar, Italian local and cosmopolitan breeds to investigate its genetic structure, and runs of homozygosity (ROH) and heterozygosity patterns. The Nero Siciliano has been reported to have the highest rate of genetic diversity among the Italian breeds, and a genetic variability comparable to that of the cosmopolitan breeds. Analyses of genomic structure and relationships underlined its proximity to wild boar, and an internal substructure probably linked to different family lines. The breed showed a low value of inbreeding estimated from ROH, and the highest diversity index among the Italian breeds, even if lower than that of the cosmopolitans. Four ROH islands in three chromosomes (SSC8, SSC11, and SSC14) and one heterozygosity-rich region (SSC1) were identified in Nero Siciliano, highlighting genomic regions related to productive QTL. Across breeds, SSC8 and SSC14 were the chromosomes with most ROH islands, with Mora Romagnola and wild boar showing the highest level of autozygosity. Chromosomes SSC2, SSC6, SSC8 and SSC13 showed the majority of runs of heterozygosity regions, mainly found in the cosmopolitan pig breeds, which reported several genes associated with health-related QTL. The outlined results can help to better identify the genomic profile of this local breed in order to plan matings, maintain adequate internal diversity and exploit the production system.

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.

Expression profiles of the CD274 and PLEKHH2 gene and association of its polymorphism with hematologic parameters in sheep

CD274 and PLEKHH2 genes have been identified as immune- and multiple diseases-related genes, and have recently garnered significant interest. However, their role in regulating immune functions in sheep remains largely unexplored. In this study, we aimed to investigate the effects of polymorphisms in CD274 and PLEKHH2 on hematologic parameters in 915 sheep. Our results showed that the CD274 and PLEKHH2 genes were most highly expressed in the spleen and tail fat, respectively, as determined by qRT-PCR. We also identified a G to A mutation (g 0.11858 G > A) in the exon 4 region of CD274, and a C to G mutation (g 0.38384 C > G) in the intron 8 region of PLEKH2. Association analysis revealed that CD274 g 0.11858 G > A was significantly associated with RBC, HCT, MCHC, and MCV (P < 0.05), while PLEKHH2 g 0.38384 C > G was significantly associated with HCT, MPV, MCHC, and MCV (P < 0.05). These results suggest that CD274 and PLEKHH2 genes may play a role in regulating blood physiological indicators and could be potential functional candidates for influencing immune traits in sheep breeding programs.

Comparison of two multi-trait association testing methods and sequence-based fine mapping of six additive QTL in Swiss Large White pigs

BACKGROUND

Genetic correlations between complex traits suggest that pleiotropic variants contribute to trait variation. Genome-wide association studies (GWAS) aim to uncover the genetic underpinnings of traits. Multivariate association testing and the meta-analysis of summary statistics from single-trait GWAS enable detecting variants associated with multiple phenotypes. In this study, we used array-derived genotypes and phenotypes for 24 reproduction, production, and conformation traits to explore differences between the two methods and used imputed sequence variant genotypes to fine-map six quantitative trait loci (QTL).

RESULTS

We considered genotypes at 44,733 SNPs for 5,753 pigs from the Swiss Large White breed that had deregressed breeding values for 24 traits. Single-trait association analyses revealed eleven QTL that affected 15 traits. Multi-trait association testing and the meta-analysis of the single-trait GWAS revealed between 3 and 6 QTL, respectively, in three groups of traits. The multi-trait methods revealed three loci that were not detected in the single-trait GWAS. Four QTL that were identified in the single-trait GWAS, remained undetected in the multi-trait analyses. To pinpoint candidate causal variants for the QTL, we imputed the array-derived genotypes to the sequence level using a sequenced reference panel consisting of 421 pigs. This approach provided genotypes at 16 million imputed sequence variants with a mean accuracy of imputation of 0.94. The fine-mapping of six QTL with imputed sequence variant genotypes revealed four previously proposed causal mutations among the top variants.

CONCLUSIONS

Our findings in a medium-size cohort of pigs suggest that multivariate association testing and the meta-analysis of summary statistics from single-trait GWAS provide very similar results. Although multi-trait association methods provide a useful overview of pleiotropic loci segregating in mapping populations, the investigation of single-trait association studies is still advised, as multi-trait methods may miss QTL that are uncovered in single-trait GWAS.

Comparison of blood gas parameters, ions, and glucose concentration in polish Holstein-Friesian Dairy cows at different milk production levels

Genetic selection for increased milk yield has been a key driver of dairy intensification. The modern dairy cow produces much higher amounts of milk than the cattle of several years ago, and this may have an influence on hematologic values at different stages of lactation and on cows with different levels of milk production. The purpose of the study was to investigate the variations in blood parameters such as Ht, tHb, sO₂, FO₂Hb, FCOHb, FMetHb, FHHb, pH, pCO₂, pO₂, standard HCO₃-, actual HCO₃-, BE, BE ecf, ctCO₂, BO₂, p50, and ctO₂ in cows at different milk production levels. In addition, ions such as Na+, K+ , Ca++, Ca++ (7.4), and Cl-, and AnGap and glucose were examined. Our findings indicated that differences in the examined blood parameters between low and high-production dairy cattle do exist. The most apparent differences were connected with blood pH (p < 0.01), oxygen metabolism (Ht, tHb, sO2, FO2Hb; p < 0.01), and glucose utilization (p < 0.01) The results confirm that the parameters connected with blood oxygen metabolism and glucose metabolism increase significantly in high-production animals. In conclusion, reference values should be considered in light of the lactation stage and level of milk production, because these might influence how changes should be interpreted. The main limitation of the study is the delay to analysis. However, the blood was properly stored (4C), thus changes were delayed. Anyway, it is very hard in the field practice to perform it within 5 min after the blood collection and according to studies it has low impact on clinical outcomes.

Identification of transcriptional regulatory variants in pig duodenum, liver, and muscle tissues

BACKGROUND

In humans and livestock species, genome-wide association studies (GWAS) have been applied to study the association between variants distributed across the genome and a phenotype of interest. To discover genetic polymorphisms affecting the duodenum, liver, and muscle transcriptomes of 300 pigs from 3 different breeds (Duroc, Landrace, and Large White), we performed expression GWAS between 25,315,878 polymorphisms and the expression of 13,891 genes in duodenum, 12,748 genes in liver, and 11,617 genes in muscle.

RESULTS

More than 9.68 × 1011 association tests were performed, yielding 14,096,080 significantly associated variants, which were grouped in 26,414 expression quantitative trait locus (eQTL) regions. Over 56% of the variants were within 1 Mb of their associated gene. In addition to the 100-kb region upstream of the transcription start site, we identified the importance of the 100-kb region downstream of the 3'UTR for gene regulation, as most of the cis-regulatory variants were located within these 2 regions. We also observed 39,874 hotspot regulatory polymorphisms associated with the expression of 10 or more genes that could modify the protein structure or the expression of a regulator gene. In addition, 2 motifs (5'-GATCCNGYGTTGCYG-3' and a poly(A) sequence) were enriched across the 3 tissues within the neighboring sequences of the most significant single-nucleotide polymorphisms in each cis-eQTL region.

CONCLUSIONS

The 14 million significant associations obtained in this study are publicly available and have enabled the identification of expression-associated cis-, trans-, and hotspot regulatory variants within and across tissues, thus shedding light on the molecular mechanisms of regulatory variations that shape end-trait phenotypes.

Genome-wide association study reveals the genetic determinism of serum biochemical indicators in ducks

BACKGROUND

The serum is rich in nutrients and plays an essential role in electrolyte and acid-base balance, maintaining cellular homeostasis. In addition, serum parameters have been commonly used as essential biomarkers for clinical diagnosis. However, little is known about the genetic mechanism of the serum parameters in ducks.

RESULTS

This study measured 18 serum parameters in 320 samples of the F₂ segregating population generated by Mallard × Pekin duck. The phenotypic correlations showed a high correlation between LDH, HBDH, AST, and ALT (0.59-0.99), and higher coefficients were also observed among TP, ALB, HDL-C, and CHO (0.46-0.87). And then, we performed the GWAS to reveal the genetic basis of the 18 serum biochemical parameters in ducks. Fourteen candidate protein-coding genes were identified with enzyme traits (AST, ALP, LDH, HBDH), and 3 protein-coding genes were associated with metabolism and protein-related serum parameters (UA, TG). Moreover, the expression levels of the above candidate protein-coding genes in different stages of breast muscle and different tissues were analyzed. Furthermore, the genes located within the high-LD region (r² > 0.4 and - log₁₀(P) < 4) neighboring the significant locus also remained. Finally, 86 putative protein-coding genes were used for GO and KEGG enrichment analysis, the enzyme-linked receptor protein signaling pathway and ErbB signaling pathway deserve further focus.

CONCLUSIONS

The obtained results can contribute to new insights into blood metabolism and provide new genetic biomarkers for application in duck breeding programs.

Assessing Population Structure and Signatures of Selection in Wanbei Pigs Using Whole Genome Resequencing Data

Wanbei pig (WBP) is one of the indigenous pig resources in China and has many germplasm characteristics. However, research on its genome is lacking. To assess the genomic variation, population structure, and selection signatures, we resequenced 18 WBP for the first time and performed a comprehensive analysis with resequenced data of 10 Asian wild boars. In total, 590.03 Gb of data and approximately 41 million variants were obtained. Polymorphism level (θπ) ratio and genetic differentiation (fixation index)-based cross approaches were applied, and 539 regions, which harbored 176 genes, were selected. Functional analysis of the selected genes revealed that they were associated with lipid metabolism (SCP2, APOA1, APOA4, APOC3, CD36, BCL6, ADCY8), backfat thickness (PLAG1, CACNA2D1), muscle (MYOG), and reproduction (CABS1). Overall, our results provide a valuable resource for characterizing the uniqueness of WBP and a basis for future breeding.

Estimates of genomic heritability and genome-wide association studies for blood parameters in Akkaraman sheep

The aim of this study was to estimate genomic heritability and the impact that genetic backgrounds have on blood parameters in Akkaraman sheep by conducting genome-wide association studies and regional heritability mapping analysis. Genomic heritability estimates for blood parameters ranged from 0.00 to 0.55, indicating that measured phenotypes have a low to moderate heritability. A total of 7 genome- and 13 chromosome-wide significant SNPs were associated with phenotypic changes in 15 blood parameters tested. Accordingly, SCN7A, SCN9A, MYADM-like, CCDC67, ITGA9, MGAT5, SLC19A1, AMPH, NTRK2, MSRA, SLC35F3, SIRT6, CREB3L3, and NAV3 genes as well as three undefined regions (LOC101117887, LOC106991526 and LOC105608461) were suggested as candidates. Most of the identified genes were involved in basic biological processes that are essential to immune system function and cellular growth; specific functions include cellular transport, histone deacetylation, cell differentiation, erythropoiesis, and endocytosis. The top significant SNP for HCT, MCH, and MCHC was found within a genomic region mainly populated by the MYADM-like gene family. This region was previously suggested to be under historical selection pressure in many sheep breeds from various parts of the world. These results have implications on animal breeding program studies due to the effect that the genetic background has on blood parameters, which underlying many productive and wellness related traits.

Lost pigs of Angola: Whole genome sequencing reveals unique regions of selection with emphasis on metabolism and feed efficiency

Angola, in the western coast of Africa, has been through dramatic social events that have led to the near-disappearance of native swine populations, and the recent introduction of European exotic breeds has also contributed to the erosion of this native swine repertoire. In an effort to investigate the genetic basis of native pigs in Angola (ANG) we have generated whole genomes from animals of a remote local pig population in Huambo province, which we have compared with 78 genomes of European and Asian pig breeds as well as European and Asian wild boars that are currently in public domain. Analyses of population structure showed that ANG pigs grouped within the European cluster and were clearly separated from Asian pig breeds. Pairwise F_ST ranged from 0.14 to 0.26, ANG pigs display lower levels of genetic differentiation towards European breeds. Finally, we have identified candidate regions for selection using a complementary approach based on various methods. All results suggest that selection towards feed efficiency and metabolism has occurred. Moreover, all analysis identified CDKAL1 gene, which is related with insulin and cholesterol metabolism, as a candidate gene overlapping signatures of selection unique to ANG pigs. This study presents the first assessment of the genetic relationship between ANG pigs and other world breeds and uncovers selection signatures that may indicate adaptation features unique to this important genetic resource.

Genetic Influence on Frequencies of Myeloid-Derived Cell Subpopulations in Mouse

Differences in frequencies of blood cell subpopulations were reported to influence the course of infections, atopic and autoimmune diseases, and cancer. We have discovered a unique mouse strain B10.O20 containing extremely high frequency of myeloid-derived cells (MDC) in spleen. B10.O20 carries 3.6% of genes of the strain O20 on the C57BL/10 genetic background. It contains much higher frequency of CD11b⁺Gr1⁺ cells in spleen than both its parents. B10.O20 carries O20-derived segments on chromosomes 1, 15, 17, and 18. Their linkage with frequencies of blood cell subpopulations in spleen was tested in F₂ hybrids between B10.O20 and C57BL/10. We found 3 novel loci controlling MDC frequencies: Mydc1, 2, and 3 on chromosomes 1, 15, and 17, respectively, and a locus controlling relative spleen weight (Rsw1) that co-localizes with Mydc3 and also influences proportion of white and red pulp in spleen. Mydc1 controls numbers of CD11b⁺Gr1⁺ cells. Interaction of Mydc2 and Mydc3 regulates frequency of CD11b⁺Gr1⁺ cells and neutrophils (Gr1⁺Siglec-F^- cells from CD11b⁺ cells). Interestingly, Mydc3/Rsw1 is orthologous with human segment 6q21 that was shown previously to determine counts of white blood cells. Bioinformatics analysis of genomic sequence of the chromosomal segments bearing these loci revealed polymorphisms between O20 and C57BL/10 that change RNA stability and genes’ functions, and we examined expression of relevant genes. This identified potential candidate genes Smap1, Vps52, Tnxb, and Rab44. Definition of genetic control of MDC can help to personalize therapy of diseases influenced by these cells.

Comparison of the choice of animals for re-sequencing in two maternal pig lines

Next-generation sequencing is a promising approach for the detection of causal variants within previously identified quantitative trait loci. Because of the costs of re-sequencing experiments, this application is currently mainly restricted to subsets of animals from already genotyped populations. Imputation from a lower to a higher marker density could represent a useful complementary approach. An analysis of the literature shows that several strategies are available to select animals for re-sequencing. This study demonstrates an animal selection workflow under practical conditions. Our approach considers different data sources and limited resources such as budget and availability of sampling material. The workflow combines previously described approaches and makes use of genotype and pedigree information from a Landrace and Large White population. Genotypes were phased and haplotypes were accurately estimated with AlphaPhase. Then, AlphaSeqOpt was used to optimize selection of animals for re-sequencing, reflecting the existing diversity of haplotypes. AlphaSeqOpt and ENDOG were used to select individuals based on pedigree information and by taking into account key animals that represent the genetic diversity of the populations. After the best selection criteria were determined, a subset of 57 animals was selected for subsequent re-sequencing. In order to evaluate and assess the advantage of this procedure, imputation accuracy was assessed by setting a set of single nucleotide polymorphism (SNP) chip genotypes to missing. Accuracy values were compared to those of alternative selection scenarios and the results showed the clear benefits of a targeted selection within this practical-driven approach. Especially imputation of low-frequency markers benefits from the combined approach described here. Accuracy was increased by up to 12% compared to a randomized or exclusively haplotype-based selection of sequencing candidates.

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

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.

Environmentally enriched housing conditions affect pig welfare, immune system and gut microbiota in early life

BACKGROUND

Conventional pig housing and management conditions are associated with gastrointestinal pathophysiology and disease susceptibility in early life. Developing new strategies to reduce both therapeutic and prophylactic antibiotic use is urgent for the sustainable swine production globally. To this end, housing methodology providing effective environmental enrichment could be a promising alternative approach to reduce antibiotic usage, as it has been proven to positively influence pig welfare and immune status and reduce susceptibility to infections. It is, however, poorly understood how this enriched housing affects systemic and local pulmonary immune status and gut microbiota colonization during early life. In the present study, we compared the effects of two housing conditions, i.e., conventional housing: (CH) versus enriched housing (EH), on immune status and gut microbiota from birth until 61 days of age.

RESULTS

The expected benefits of enrichment on pig welfare were confirmed as EH pigs showed more positive behaviour, less aggression behaviour during the weaning transition and better human animal relation during the post weaning phase. Regarding the pigs' immune status, EH pigs had higher values of haemoglobin and mean corpuscular volume in haematological profiles and higher percentages of T cells and cytotoxic T cells in peripheral blood. Furthermore, EH pigs showed higher ex vivo secretion of IL1ß and TNF-α after lipopolysaccharide stimulation of whole blood than CH pigs. The structure of the developing faecal microbiota of CH and EH pigs significantly differed as early as day 12 with an increase in the relative abundance of several bacterial groups known to be involved in the production of short chain fatty acids, such as Prevotella_2, Christensenellaceae_R_7_group and Ruminococcus gauvreauii group. Furthermore, the main difference between both housing conditions post weaning was that on day 61, CH pigs had significantly larger inter-individual variation of ileal and colonic microbiota than EH pigs. In addition to housing, other intrinsic factors (e.g., sex) were associated with gut microbiota development and immune competence.

CONCLUSIONS

In addition to the known welfare benefits for pigs, environmentally enriched housing also positively drives important aspects of the development of the immune system and the establishment of gut microbiota in early life. Consequently, EH may contribute to increasing productivity of pigs and reducing antibiotic use.

Genotyping and Whole-Genome Resequencing of Welsh Sheep Breeds Reveal Candidate Genes and Variants for Adaptation to Local Environment and Socioeconomic Traits

Background

Advances in genetic tools applied to livestock breeding has prompted research into the previously neglected breeds adapted to harsh local environments. One such group is the Welsh mountain sheep breeds, which can be farmed at altitudes of 300 m above sea level but are considered to have a low productive value because of their poor wool quality and small carcass size. This is contrary to the lowland breeds which are more suited to wool and meat production qualities, but do not fare well on upland pasture. Herein, medium-density genotyping data from 317 individuals representing 15 Welsh sheep breeds were used alongside the whole-genome resequencing data of 14 breeds from the same set to scan for the signatures of selection and candidate genetic variants using haplotype- and SNP-based approaches.

Results

Haplotype-based selection scan performed on the genotyping data pointed to a strong selection in the regions of GBA3, PPARGC1A, APOB, and PPP1R16B genes in the upland breeds, and RNF24, PANK2, and MUC15 in the lowland breeds. SNP-based selection scan performed on the resequencing data pointed to the missense mutations under putative selection relating to a local adaptation in the upland breeds with functions such as angiogenesis (VASH1), anti-oxidation (RWDD1), cell stress (HSPA5), membrane transport (ABCA13 and SLC22A7), and insulin signaling (PTPN1 and GIGFY1). By contrast, genes containing candidate missense mutations in the lowland breeds are related to cell cycle (CDK5RAP2), cell adhesion (CDHR3), and coat color (MC1R).

Conclusion

We found new variants in genes with potentially functional consequences to the adaptation of local sheep to their environments in Wales. Knowledge of these variations is important for improving the adaptative qualities of UK and world sheep breeds through a marker-assisted selection.

Single-marker and haplotype-based genome-wide association studies for the number of teats in two heavy pig breeds

The number of teats is a reproductive‐related trait of great economic relevance as it affects the mothering ability of the sows and thus the number of properly weaned piglets. Moreover, genetic improvement of this trait is fundamental to parallelly help the selection for increased litter size. We present the results of single‐marker and haplotypes‐based genome‐wide association studies for the number of teats in two large cohorts of heavy pig breeds (Italian Large White and Italian Landrace) including 3990 animals genotyped with the 70K GGP Porcine BeadChip and other 1927 animals genotyped with the Illumina PorcineSNP60 BeadChip. In the Italian Large White population, genome scans identified three genome regions (SSC7, SSC10, and SSC12) that confirmed the involvement of the VRTN gene (as we previously reported) and highlighted additional loci known to affect teat counts, including the FRMD4A and HOXB1 gene regions. A different picture emerged in the Italian Landrace population, with a total of 12 genome regions in eight chromosomes (SSC3, SSC6, SSC8, SSC11, SSC13, SSC14, SSC15, and SSC16) mainly detected via the haplotype‐based genome scan. The most relevant QTL was close to the ARL4C gene on SSC15. Markers in the VRTN gene region were not significant in the Italian Landrace breed. The use of both single‐marker and haplotype‐based genome‐wide association analyses can be helpful to exploit and dissect the genome of the pigs of different populations. Overall, the obtained results supported the polygenic nature of the investigated trait and better elucidated its genetic architecture in Italian heavy pigs.

Genomic Regions and Candidate Genes Linked to Capped Hock in Pig

Capped hock affects the exterior of pedigree pigs, making them unsalable and resulting in a negative impact on the efficiency of pig-breeding centers. The purpose of this paper was to carry out pilot studies aimed at finding genomic regions and genes linked to the capped hock in pigs. The studies were carried out on Landrace pigs (n = 75) and Duroc pigs (n = 70). To identify genomic regions linked to capped hock in pigs, we used smoothing FST statistics. Genotyping was performed with GeneSeek^® GGP Porcine HD Genomic Profiler v1 (Illumina Inc, San Diego, CA, USA). The research results showed 70 SNPs linked to capped hock in Landrace (38 SNPs) and Duroc (32 SNPs). The identified regions overlapped with QTLs related with health traits (blood parameters) and meat and carcass traits (fatness). In total, 31 genes were identified (i.e., 17 genes in Landrace, 14 genes in Durocs). Three genes appeared in both the Landrace and Duroc groups, including A2ML1 (SSC5), ROBO2 (SSC13), and MSI1 (SSC14). We identified genomic regions directly or indirectly linked to capped hock, which thus might contribute to identifying genetic variants and using them as genetic markers in pig breeding.

Eosinophils may serve as CEA-secreting cells for allergic bronchopulmonary aspergillosis (ABPA) patients

Allergic bronchopulmonary aspergillosis (ABPA) is a condition characterized by an exaggerated response of the immune system to the fungus Aspergillus. This study aimed to assess the relationship between carcinoembryonic antigen (CEA) and eosinophils in ABPA patients. We describes a case of a 50-year-old patient who was diagnosed with ABPA presenting with high level of CEA and eosinophils. Besides,we used immunohistochemistry and immunofluorescence to identify eosinophils and CEA in sections which were obtained by Endobronchial ultrasound-guided transbronchial lung biopsy aspiration (EBUS-TBLB). The sections were then visualized using confocal microscopy. We also retrospectively analyzed a cohort of 37 ABPA patients between January 2013 and December 2019 in our hospital. We found the patient whose serum CEA levels were consistent with eosinophils during the follow-up (r = 0.929, P = 0.022). The positive expression of CEA and abnormal expression of eosinophils was higher in the ABPA tissue compared to the normal lung tissue. The co-localization was represented as pixels containing both red and green color in the image (with various shades of orange and yellow) which signified that eosinophils were immunohistochemically positive for CEA. Patients with higher levels of eosinophils had higher levels of CEA in the serum (P < 0.001). The results of Pearson correlation analysis showed that the levels of eosinophils were positively correlated with serum CEA levels (r = 0.459 and r = 0.506, P = 0.004 and P = 0.001). Serum CEA level is elevated in ABPA patients. The elevated serum CEA level was shown to be normalized after treatment. Increased CEA levels in ABPA patients may be positively correlated with eosinophil levels, and eosinophils may be served as CEA-secreting cells in patients with ABPA.

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.

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.

Haplotype-based genome-wide association studies reveal new loci for haematological and clinical-biochemical parameters in Large White pigs

We report haplotype-based GWASs for 33 blood parameters measured in 843 Italian Large White pigs. In the single-trait analysis, a total of 30 QTL for number of basophils, six erythrocyte traits (haemoglobin, haematocrit, mean corpuscular haemoglobin, mean corpuscular haemoglobin concentration, mean corpuscular volume and red blood cell count) and two clinical-biochemical traits (alkaline phosphatase and Ca²⁺ contents) were identified. In the multiple-trait analysis, a total of five QTL affected three different clusters of traits. Only four of these QTL were already reported in the single-marker and multi-marker GWASs we previously carried out on the same pig population. QTL on SSC11 and SSC17 showed effects on multiple traits. These results further dissected the genetic architecture of parameters that could be used as proxies in breeding programmes for more complex traits. In addition, these results might help to better define the pig as an animal model for several blood-related biological functions.

Genome-Wide Association Study Using Fix-Length Haplotypes and Network Analysis Revealed New Candidate Genes for Nematode Resistance and Body Weight in Blackface Lambs

The objectives of this study were to identify genomic regions by Bayesian methods (BayesA, BayesB, or BayesN) that fit fixed-length haplotypes or SNPs using GenSel. Covariates for haplo-type alleles of five lengths (125, 250, 500 kb, 1 or 2 Mb) were generated, and rare haplotypes were removed at three thresholds (1, 5, or 10%). Subsequently, we performed gene network analyses to investigate the biological processes shared by genes that were identified for the same across traits. Genotypes at 41,034 SNPs that were common on OvineSNP50 panel were phased for 751 Scottish Blackface (SBF) lambs. This is the first study to quantify the proportion of genetic variance using haplotypes across the whole genome in an SBF population. The genetic variance explained of haplotype-based GWAS was higher than that of SNP-based GWAS in across traits studied. In this population, fitting 500kb haplotypes with a 1% frequency threshold resulted in the highest proportion of genetic variance explained for nematode resistance and fitting 2Mb haplotypes with a 10% frequency threshold improved genetic variance explained for body weight comparable to fitting SNPs by BayesB. Candidate genes, including CXCR4, STAT4, CCL1, CCL2, CCL3, CCL5, CCL8, CCL16, CCL18, CARMIL2, and HSPA14 were identified for nematode resistance and ADH5, PPP3CA, and FABP4 for body weight traits. Network analysis provided annotation results linking to all identified candidate genes. This study supported previous results from GWAS of nematode resistance and body weight and revealed additional regions in the ovine genome associated with these economically important traits. These results suggest that network analysis can provide new information regarding biological mechanisms and genes leading to complex phenotypes, like nematode resistance and body weight of lamb.

Genome-wide association studies for the number of teats and teat asymmetry patterns in Large White pigs

The number of teats is a morphological trait that influences the mothering ability of the sows and thus their reproduction performances. In this study, we carried out GWASs for the total number of teats and other 12 related parameters in 821 Italian Large White heavy pigs. All pigs were genotyped with the Illumina PorcineSNP60 BeadChip array. For four investigated parameters (total number of teats, the number of teats of the left line, the number of teats of the right line and the maximum number of teats comparing the two sides), significant markers were identified on SSC7, in the region of the vertnin (VRTN) gene. Significant markers for the numbers of posterior teats and the absolute difference between anterior and posterior teat numbers were consistently identified on SSC6. The most significant SNP for these parameters was an intron variant in the TOX high mobility group box family member 3 (TOX3) gene. For the other four parameters (absolute difference between the two sides; anterior teats; the ratio between the posterior and the anterior number of teats; and the absence or the presence of extra teats) only suggestively significant markers were identified on several other chromosomes. This study further supported the role of the VRTN gene region in affecting the recorded variability of the number of teats in the Italian Large White pig population and identified a genomic region potentially affecting the biological mechanisms controlling the developmental programme of morphological features in pigs.

Conservation Genomic Analysis of the Croatian Indigenous Black Slavonian and Turopolje Pig Breeds

The majority of the nearly 400 existing local pig breeds are adapted to specific environments and human needs. The demand for large production quantities and the industrialized pig production have caused a rapid decline of many local pig breeds in recent decades. Black Slavonian pig and Turopolje pig, the latter highly threatened, are the two Croatian local indigenous breeds typically grown in extensive or semi-intensive systems. In order to guide a long-term breeding program to prevent the disappearance of these breeds, we analyzed their genetic diversity, inbreeding level and relationship with other local breeds across the world, as well as modern breeds and several wild populations, using high throughput genomic data obtained using the Illumina Infinium PorcineSNP60 v2 BeadChip. Multidimensional scaling analysis positioned Black Slavonian pigs close to the UK/North American breeds, while the Turopolje pig clustered within the Mediterranean breeds. Turopolje pig showed a very high inbreeding level (F_ROH>4Mb = 0.400 and F_ROH>8Mb = 0.332) that considerably exceeded the level of full-sib mating, while Black Slavonian pig showed much lower inbreeding (F_ROH>4Mb = 0.098 and F_ROH>8Mb = 0.074), indicating a planned mating strategy. In Croatian local breeds we identified several genome regions showing adaptive selection signals that were not present in commercial breeds. The results obtained in this study reflect the current genetic status and breeding management of the two Croatian indigenous local breeds. Given the small populations of both breeds, a controlled management activity has been implemented in Black Slavonian pigs since their commercial value has been recognized. In contrast, the extremely high inbreeding level observed in Turopolje pig argues for an urgent conservation plan with a long-term, diversity-oriented breeding program.

Genome-wide association studies for iris pigmentation and heterochromia patterns in Large White pigs

Eye colour genetics have been extensively studied in humans since the rediscovery of Mendel's laws. This trait was first interpreted using simplistic genetic models but soon it was realised that it is more complex. In this study, we analysed eye colour variability in a Large White pig population (n = 897) and report the results of GWASs based on several comparisons including pigs having four main eye colour categories (three with both pigmented eyes of different brown grades: pale, 17.9%; medium, 14.8%; and dark, 54.3%; another one with both eyes completely depigmented, 3.8%) and heterochromia patterns (heterochromia iridis - depigmented iris sectors in pigmented irises, 3.2%; heterochromia iridum - one whole eye iris of depigmented phenotype and the other eye with the iris completely pigmented, 5.9%). Pigs were genotyped with the Illumina PorcineSNP60 BeadChip and GEMMA was used for the association analyses. The results indicated that SLC45A2 (on chromosome 16, SSC16), EDNRB (SSC11) and KITLG (SSC5) affect the different grades of brown pigmentation of the eyes, the bilateral eye depigmentation defect and the heterochromia iridis defect recorded in this white pig population respectively. These genes are involved in several mechanisms affecting pigmentation. Significant associations for the eye depigmented patterns were also identified for SNPs on two SSC4 regions (including two candidate genes: NOTCH2 and PREX2) and on SSC6, SSC8 and SSC14 (including COL17A1 as candidate gene). This study provided useful information to understand eye pigmentation mechanisms, further valuing the pig as animal model to study complex phenotypes in humans.

Genome-Wide Epistatic Interaction Networks Affecting Feed Efficiency in Duroc and Landrace Pigs

Interactions among genomic loci have often been overlooked in genome-wide association studies, revealing the combinatorial effects of variants on phenotype or disease manifestation. Unexplained genetic variance, interactions among causal genes of small effects, and biological pathways could be identified using a network biology approach. The main objective of this study was to determine the genome-wide epistatic variants affecting feed efficiency traits [feed conversion ratio (FCR) and residual feed intake (RFI)] based on weighted interaction SNP hub (WISH-R) method. Herein, we detected highly interconnected epistatic SNP modules, pathways, and potential biomarkers for the FCR and RFI in Duroc and Landrace purebreds considering the whole population, and separately for low and high feed efficient groups. Highly interacting SNP modules in Duroc (1,247 SNPs) and Landrace (1,215 SNPs) across the population and for low feed efficient (Duroc-80 SNPs, Landrace-146 SNPs) and high feed efficient group (Duroc-198 SNPs, Landrace-232 SNPs) for FCR and RFI were identified. Gene and pathway analyses identified ABL1, MAP3K4, MAP3K5, SEMA6A, KITLG, and KAT2B from chromosomes 1, 2, 5, and 13 underlying ErbB, Ras, Rap1, thyroid hormone, axon guidance pathways in Duroc. GABBR2, GNA12, and PRKCG genes from chromosomes 1, 3, and 6 pointed towards thyroid hormone, cGMP-PKG and cAMP pathways in Landrace. From Duroc low feed efficient group, the TPK1 gene was found involved with thiamine metabolism, whereas PARD6G, DLG2, CRB1 were involved with the hippo signaling pathway in high feed efficient group. PLOD1 and SETD7 genes were involved with lysine degradation in low feed efficient group in Landrace, while high feed efficient group pointed to genes underpinning valine, leucine, isoleucine degradation, and fatty acid elongation. Some SNPs and genes identified are known for their association with feed efficiency, others are novel and potentially provide new avenues for further research. Further validation of epistatic SNPs and genes identified here in a larger cohort would help to establish a framework for modelling epistatic variance in future methods of genomic prediction, increasing the accuracy of estimated genetic merit for FE and helping the pig breeding industry.

Integrative analysis of blood and gut microbiota data suggests a non-alcoholic fatty liver disease (NAFLD)-related disorder in French SLAdd minipigs

Minipigs are a group of small-sized swine lines, which show a broad range of phenotype variation and which often tend to be obese. The SLA^dd (DD) minipig line was created by the NIH and selected as homozygous at the SLA locus. It was brought to France more than 30 years ago and maintained inbred ever since. In this report, we characterized the physiological status of a herd of French DD pigs by measuring intermediate phenotypes from blood and faeces and by using Large White (LW) pigs as controls. Three datasets were produced, i.e. complete blood counts (CBCs), microarray-based blood transcriptome, and faecal microbiota obtained by 16S rRNA sequencing. CBCs and expression profiles suggested a non-alcoholic fatty liver disease (NAFLD)-related pathology associated to comorbid cardiac diseases. The characterization of 16S sequencing data was less straightforward, suggesting only a potential weak link to obesity. The integration of the datasets identified several fine-scale associations between CBCs, gene expression, and faecal microbiota composition. NAFLD is a common cause of chronic liver disease in Western countries and is linked to obesity, type 2 diabetes mellitus and cardiac pathologies. Here we show that the French DD herd is potentially affected by this syndrome.

A machine learning approach for the identification of population-informative markers from high-throughput genotyping data: application to several pig breeds

Single nucleotide polymorphisms (SNPs) able to describe population differences can be used for important applications in livestock, including breed assignment of individual animals, authentication of mono-breed products and parentage verification among several other applications. To identify the most discriminating SNPs among thousands of markers in the available commercial SNP chip tools, several methods have been used. Random forest (RF) is a machine learning technique that has been proposed for this purpose. In this study, we used RF to analyse PorcineSNP60 BeadChip array genotyping data obtained from a total of 2737 pigs of 7 Italian pig breeds (3 cosmopolitan-derived breeds: Italian Large White, Italian Duroc and Italian Landrace, and 4 autochthonous breeds: Apulo-Calabrese, Casertana, Cinta Senese and Nero Siciliano) to identify breed informative and reduced SNP panels using the mean decrease in the Gini Index and the Mean Decrease in Accuracy parameters with stability evaluation. Other reduced informative SNP panels were obtained using Delta, Fixation index and principal component analysis statistics, and their performances were compared with those obtained using the RF-defined panels using the RF classification method and its derived Out Of Bag rates and correct prediction proportions. Therefore, the performances of a total of six reduced panels were evaluated. The correct assignment of the animals to its breed was close to 100% for all tested approaches. Porcine chromosome 8 harboured the largest number of selected SNPs across all panels. Many SNPs were included in genomic regions in which previous studies identified signatures of selection or genes (e.g. ESR1, KITL and LCORL) that could contribute to explain, at least in part, phenotypically or economically relevant traits that might differentiate cosmopolitan and autochthonous pig breeds. Random forest used as preselection statistics highlighted informative SNPs that were not the same as those identified by other methods. This might be due to specific features of this machine learning methodology. It will be interesting to explore if the adaptation of RF methods for the identification of selection signature regions could be able to describe population-specific features that are not captured by other approaches.