Imputation-Based Whole-Genome Sequence Association Study Reveals Constant and Novel Loci for Hematological Traits in a Large-Scale Swine F2 Resource Population

Genome-wide association studies for loin muscle area, loin muscle depth and backfat thickness in DLY pigs

This study aimed at identifying genes associated with loin muscle area (LMA), loin muscle depth (LMD) and backfat thickness (BFT). We performed single-trait and multi-trait genome-wide association studies (GWASs) after genotyping 685 Duroc × (Landrace × Yorkshire) (DLY) pigs using the Geneseek Porcine 50K SNP chip. In the single-trait GWASs, we identified two, eight and two significant SNPs associated with LMA, LMD and BFT, respectively, and searched genes within the 1 Mb region near the significant SNPs with relevant functions as candidate genes. Consequently, we identified one (DOCK5), three (PID1, PITX2, ELOVL6) and three (CCR1, PARP14, CASR) promising candidate genes for LMA, LMD and BFT, respectively. Moreover, the multi-trait GWAS identified four significant SNPs associated with the three traits. In conclusion, the GWAS analysis of LMA, LMD and BFT in a DLY pig population identified several associated SNPs and candidate genes, further deepening our understanding of the genetic basis of these traits, and they may be useful for marker-assisted selection to improve the three traits in DLY pigs.

Uncovering the candidate genes related to sheep body weight using multi-trait genome-wide association analysis

In sheep, body weight is an economically important trait. This study sought to map genetic loci related to weaning weight and yearling weight. To this end, a single-trait and multi-trait genome-wide association study (GWAS) was performed using a high-density 600 K single nucleotide polymorphism (SNP) chip. The results showed that 43 and 56 SNPs were significantly associated with weaning weight and yearling weight, respectively. A region associated with both weaning and yearling traits (OARX: 6.74-7.04 Mb) was identified, suggesting that the same genes could play a role in regulating both these traits. This region was found to contain three genes (TBL1X, SHROOM2 and GPR143). The most significant SNP was Affx-281066395, located at 6.94 Mb (p = 1.70 × 10-17), corresponding to the SHROOM2 gene. We also identified 93 novel SNPs elated to sheep weight using multi-trait GWAS analysis. A new genomic region (OAR10: 76.04-77.23 Mb) with 22 significant SNPs were discovered. Combining transcriptomic data from multiple tissues and genomic data in sheep, we found the HINT1, ASB11 and GPR143 genes may involve in sheep body weight. So, multi-omic anlaysis is a valuable strategy identifying candidate genes related to body weight.

Functional genomics of reproduction in pigs: Are we there yet?

Reproductive failure is the main reason for culling females in swine herds and is both a financial and sustainability issue. Because reproductive traits are complex and lowly to moderately heritable, genomic selection within populations can achieve substantial genetic gain in reproductive efficiency. A better understanding of the physiological components affecting the expression of these traits will facilitate greater understanding of the genes affecting reproductive traits and is necessary to improve and optimize management strategies to maximize reproductive success of gilts and sows. Large-scale genotyping with single-nucleotide polymorphism (SNP) arrays are used for genome-wide association studies (GWAS) and have facilitated identification of positional candidate genes. Transcriptomic data can be used to weight SNP for GWAS and could lead to previously unidentified candidate genes. Resequencing and fine mapping of candidate genes are necessary to identify putative functional variants and some of these have been incorporated into new genotyping arrays. Sequence imputation and genotype by sequence are newer strategies that could reveal novel functional mutations. In this study, these approaches are discussed. Advantages and limitations are highlighted where additional research is needed.

Integrated Single-Trait and Multi-Trait GWASs Reveal the Genetic Architecture of Internal Organ Weight in Pigs

Internal organ weight is an essential indicator of growth status as it reflects the level of growth and development in pigs. However, the associated genetic architecture has not been well explored because phenotypes are difficult to obtain. Herein, we performed single-trait and multi-trait genome-wide association studies (GWASs) to map the genetic markers and genes associated with six internal organ weight traits (including heart weight, liver weight, spleen weight, lung weight, kidney weight, and stomach weight) in 1518 three-way crossbred commercial pigs. In summation, single-trait GWASs identified a total of 24 significant single- nucleotide polymorphisms (SNPs) and 5 promising candidate genes, namely, TPK1, POU6F2, PBX3, UNC5C, and BMPR1B, as being associated with the six internal organ weight traits analyzed. Multi-trait GWAS identified four SNPs with polymorphisms localized on the APK1, ANO6, and UNC5C genes and improved the statistical efficacy of single-trait GWASs. Furthermore, our study was the first to use GWASs to identify SNPs associated with stomach weight in pigs. In conclusion, our exploration of the genetic architecture of internal organ weights helps us better understand growth traits, and the key SNPs identified could play a potential role in animal breeding programs.

Distinct traces of mixed ancestry in western commercial pig genomes following gene flow from Chinese indigenous breeds

Studying gene flow between different livestock breeds will benefit the discovery of genes related to production traits and provide insight into human historical breeding. Chinese pigs have played an indispensable role in the breeding of Western commercial pigs. However, the differences in the timing and volume of the contribution of pigs from different Chinese regions to Western pigs are not yet apparent. In this paper, we combine the whole-genome sequencing data of 592 pigs from different studies and illustrate patterns of gene flow from Chinese pigs into Western commercial pigs. We describe introgression patterns from four distinct Chinese indigenous groups into five Western commercial groups. There were considerable differences in the number and length of the putative introgressed segments from Chinese pig groups that contributed to Western commercial pig breeds. The contribution of pigs from different Chinese geographical locations to a given western commercial breed varied more than that from a specific Chinese pig group to different Western commercial breeds, implying admixture within Europe after introgression. Within different Western commercial lines from the same breed, the introgression patterns from a given Chinese pig group seemed highly conserved, suggesting that introgression of Chinese pigs into Western commercial pig breeds mainly occurred at an early stage of breed formation. Finally, based on analyses of introgression signals, allele frequencies, and selection footprints, we identified a ∼2.65 Mb Chinese-derived haplotype under selection in Duroc pigs (CHR14: 95.68-98.33 Mb). Functional and phenotypic studies demonstrate that this PRKG1 haplotype is related to backfat and loin depth in Duroc pigs. Overall, we demonstrate that the introgression history of domestic pigs is complex and that Western commercial pigs contain distinct traces of mixed ancestry, likely derived from various Chinese pig breeds.

Imputation to whole-genome sequence and its use in genome-wide association studies for pork colour traits in crossbred and purebred pigs

Imputed whole-genome sequence (WGS) has been proposed to improve genome-wide association studies (GWAS), since all causative mutations responsible for phenotypic variation are expected to be present in the data. This approach was applied on a large number of purebred (PB) and crossbred (CB) pigs for 18 pork color traits to evaluate the impact of using imputed WGS relative to medium-density marker panels. The traits included Minolta A*, B*, and L* for fat (FCOL), quadriceps femoris muscle (QFCOL), thawed loin muscle (TMCOL), fresh ham gluteus medius (GMCOL), ham iliopsoas muscle (ICOL), and longissimus dorsi muscle on the fresh loin (FMCOL). Sequence variants were imputed from a medium-density marker panel (61K for CBs and 50K for PBs) in all genotyped pigs using BeagleV5.0. We obtained high imputation accuracy (average of 0.97 for PBs and 0.91 for CBs). GWAS were conducted for three datasets: 954 CBs and 891 PBs, and the combined CBs and PBs. For most traits, no significant associations were detected, regardless of panel density or population type. However, quantitative trait loci (QTL) regions were only found for a few traits including TMCOL Minolta A* and GMCOL Minolta B* (CBs), FMCOL Minolta B*, FMCOL Minolta L*, and ICOL Minolta B* (PBs) and FMCOL Minolta A*, FMCOL Minolta B*, GMCOL Minolta B*, and ICOL Minolta B* (Combined dataset). More QTL regions were identified with WGS (n = 58) relative to medium-density marker panels (n = 22). Most of the QTL were linked to previously reported QTLs or candidate genes that have been previously reported to be associated with meat quality, pH and pork color; e.g., VIL1, PRKAG3, TTLL4, and SLC11A1, USP37. CTDSP1 gene on SSC15 has not been previously associated with meat color traits in pigs. The findings suggest any added value of WGS was only for detecting novel QTL regions when the sample size is sufficiently large as with the Combined dataset in this study. The percentage of phenotypic variance explained by the most significant SNPs also increased with WGS compared with medium-density panels. The results provide additional insights into identification of a number of candidate regions and genes for pork color traits in different pig populations.

Genome-Wide Association Study Identified a Quantitative Trait Locus and Two Candidate Genes on Sus scrofa Chromosome 2 Affecting Vulvar Traits of Suhuai Pigs

Vulvar size and angle are meaningful traits in pig production. Sows with abnormal vulva generally show reproductive disorders. In order to excavate candidate loci and genes associated with pig’s vulvar traits, 270 Suhuai pigs with vulvar phenotype were genotyped by a porcine single nucleotide polymorphisms (SNP) Chip. Then, Chip data were imputed using resequenced data of 30 Suhuai pigs as a reference panel. Next, we estimated the heritability and performed a genome-wide association study (GWAS) for vulvar traits. The heritabilities for the traits vulvar length (VL), vulvar width (VW) and vulvar angle (VA) in this pig population were 0.23, 0.32 and 0.22, respectively. GWAS based on Chip data identified nine significant SNPs on the Sus scrofa chromosomes (SSC) 2, 7, 9 and 13 for VL or VW. GWAS based on imputed data identified 11 new quantitative trait loci (QTL) on SSC1, 2, 7, 8, 9, 11, 13, 16 and 17 for VL or VW. The most significant QTL for VL on SSC2 were refined to a 3.48–3.97 Mb region using linkage disequilibrium and linkage analysis (LDLA). In this refined region, FGF19 and CCND1, involved in the development of the reproductive tract, cell growth and vulvar cancer, could be new candidate genes affecting VL. Our results provided potential genetic markers for the breeding of vulvar traits in pigs and deepened the understanding of the genetic mechanism of vulvar traits.

Scan of the endogenous retrovirus sequences across the swine genome and survey of their copy number variation and sequence diversity among various Chinese and Western pig breeds

In pig-to-human xenotransplantation, the transmission risk of porcine endogenous retroviruses (PERVs) is of great concern. However, the distribution of PERVs in pig genomes, their genetic variation among Eurasian pigs, and their evolutionary history remain unclear. We scanned PERVs in the current pig reference genome (assembly Build 11.1), and identified 36 long complete or near-complete PERVs (lcPERVs) and 23 short incomplete PERVs (siPERVs). Besides three known PERVs (PERV-A, -B, and -C), four novel types (PERV-JX1, -JX2, -JX3, and -JX4) were detected in this study. According to evolutionary analyses, the newly discovered PERVs were more ancient, and PERV-Bs probably experienced a bottleneck ~0.5 million years ago (Ma). By analyzing 63 high-quality porcine whole-genome resequencing data, we found that the PERV copy numbers in Chinese pigs were lower (32.0±4.0) than in Western pigs (49.1±6.5). Additionally, the PERV sequence diversity was lower in Chinese pigs than in Western pigs. Regarding the lcPERV copy numbers, PERV-A and -JX2 in Western pigs were higher than in Chinese pigs. Notably, Bama Xiang (BMX) pigs had the lowest PERV copy number (27.8±5.1), and a BMX individual had no PERV-C and the lowest PERV copy number (23), suggesting that BMX pigs were more suitable for screening and/or modification as xenograft donors. Furthermore, we identified 451 PERV transposon insertion polymorphisms (TIPs), of which 86 were shared by all 10 Chinese and Western pig breeds. Our findings provide systematic insights into the genomic distribution, variation, evolution, and possible biological function of PERVs.

Genome-Wide Association Analysis Reveals Genetic Loci and Candidate Genes for Chest, Abdominal, and Waist Circumferences in Two Duroc Pig Populations

Chest circumference (CC), abdominal circumference (AC), and waist circumference (WC) are regarded as important indicators for improving economic traits because they can reflect the growth and physiological status in pigs. However, the genetic architecture of CC, AC, and WC is still elusive. Here, we performed single-trait and multi-trait genome-wide association studies (GWASs) for CC, AC, and WC in 2,206 American origin Duroc (AOD) and 2,082 Canadian origin Duroc (COD) pigs. As a result, one novel quantitative trait locus (QTL) on Sus scrofa chromosome (SSC) one was associated with CC and AC in COD pigs, which spans 6.92 Mb (from 170.06 to 176.98 Mb). Moreover, multi-trait GWAS identified 21 significant SNPs associated with the three conformation traits, indicating the multi-trait GWAS is a powerful statistical approach that uncovers pleiotropic locus. Finally, the three candidate genes (ITGA11, TLE3, and GALC) were selected that may play a role in the conformation traits. Further bioinformatics analysis indicated that the candidate genes for the three conformation traits mainly participated in sphingolipid metabolism and lysosome pathways. For all we know, this study was the first GWAS for WC in pigs. In general, our findings further reveal the genetic architecture of CC, AC, and WC, which may offer a useful reference for improving the conformation traits in pigs.

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

An imputed whole-genome sequence-based GWAS approach pinpoints causal mutations for complex traits in a specific swine population

Sequencing-based genome-wide association studies (GWAS) have facilitated the identification of causal associations between genetic variants and traits in diverse species. However, it is cost-prohibitive for the majority of research groups to sequence a large number of samples. Here, we carried out genotype imputation to increase the density of single nucleotide polymorphisms in a large-scale Swine F₂ population using a reference panel including 117 individuals, followed by a series of GWAS analyses. The imputation accuracies reached 0.89 and 0.86 for allelic concordance and correlation, respectively. A quantitative trait nucleotide (QTN) affecting the chest vertebrate was detected directly, while the investigation of another QTN affecting the residual glucose failed due to the presence of similar haplotypes carrying wild-type and mutant allelesin the reference panel used in this study. A high imputation accuracy was confirmed by Sanger sequencing technology for the most significant loci. Two candidate genes, CPNE5 and MYH3, affecting meat-related traits were proposed. Collectively, we illustrated four scenarios in imputation-based GWAS that may be encountered by researchers, and our results will provide an extensive reference for future genotype imputation-based GWAS analyses in the future.

Accelerated deciphering of the genetic architecture of agricultural economic traits in pigs using a low-coverage whole-genome sequencing strategy

BACKGROUND

Uncovering the genetic architecture of economic traits in pigs is important for agricultural breeding. However, high-density haplotype reference panels are unavailable in most agricultural species, limiting accurate genotype imputation in large populations. Moreover, the infinitesimal model of quantitative traits implies that weak association signals tend to be spread across most of the genome, further complicating the genetic analysis. Hence, there is a need to develop new methods for sequencing large cohorts without large reference panels.

RESULTS

We describe a Tn5-based highly accurate, cost- and time-efficient, low-coverage sequencing method to obtain 11.3 million whole-genome single-nucleotide polymorphisms in 2,869 Duroc boars at a mean depth of 0.73×. On the basis of these single-nucleotide polymorphisms, a genome-wide association study was performed, resulting in 14 quantitative trait loci (QTLs) for 7 of 21 important agricultural traits in pigs. These QTLs harbour genes, such as ABCD4 for total teat number and HMGA1 for back fat thickness, and provided a starting point for further investigation. The inheritance models of the different traits varied greatly. Most follow the minor-polygene model, but this can be attributed to different reasons, such as the shaping of genetic architecture by artificial selection for this population and sufficiently interconnected minor gene regulatory networks.

CONCLUSIONS

Genome-wide association study results for 21 important agricultural traits identified 14 QTLs/genes and showed their genetic architectures, providing guidance for genetic improvement harnessing genomic features. The Tn5-based low-coverage sequencing method can be applied to large-scale genome studies for any species without a good reference panel and can be used for agricultural breeding.

A further look at quantitative trait loci for growth and fatness traits in a White Duroc × Erhualian F3 intercross population

Genetic analysis of porcine growth and fatness traits is beneficial to the swine industry and provides a reference to understand human obesity. Here, we obtained 29 growth and fatness traits for 473 individuals from a White Duroc × Erhualian F₃ intercross population. Basic statistical analyses showed that: (1) Positive correlations between different-stage body weights were detected, the shorter the time interval the stronger the correlation. (2) Strong correlations existed in the paired fatness traits. (3) With the growth of age, the correlation between fatness and body weight was increasing. All pigs were genotyped by Illumina 50 K SNP chips and their whole-genome genotypes were imputed referred to 109 re-sequencing data. We performed common and imputation-based GWASs for these traits. Two genome-wide significant loci on swine chromosome (SSC) 4 and 7 were repeatedly detected. The strongest association (P = 3.24 × 10^-19) was detected at 31.96 Mb on SSC7 for leaf fat weight. On this locus, seven major haplotypes were identified, of which two were novel and had an increasing-fatness effect. In the imputation-based GWAS, three new loci were identified. Our findings provide further insights into and enhance our understanding of genetic mechanism of porcine growth and fat deposition.

Identification of Major Loci and Candidate Genes for Meat Production-Related Traits in Broilers

Background

Carcass traits are crucial characteristics of broilers. However, the underlying genetic mechanisms are not well understood. In the current study, significant loci and major-effect candidate genes affecting nine carcass traits related to meat production were analyzed in 873 purebred broilers using an imputation-based genome-wide association study.

Results

The heritability estimates of nine carcass traits, including carcass weight, thigh muscle weight, and thigh muscle percentage, were moderate to high and ranged from 0.21 to 0.39. Twelve genome-wide significant SNPs and 118 suggestively significant SNPs of 546,656 autosomal variants were associated with carcass traits. All SNPs for six weight traits (body weight at 42 days of age, carcass weight, eviscerated weight, whole thigh weight, thigh weight, and thigh muscle weight) were clustered around the 24.08 Kb region (GGA24: 5.73–5.75 Mb) and contained only one candidate gene (DRD2). The most significant SNP, rs15226023, accounted for 4.85–7.71% of the estimated genetic variance of the six weight traits. The remaining SNPs for carcass composition traits (whole thigh percentage and thigh percentage) were clustered around the 42.52 Kb region (GGA3: 53.03–53.08 Mb) and contained only one candidate gene (ADGRG6). The most significant SNP in this region, rs13571431, accounted for 11.89–13.56% of the estimated genetic variance of two carcass composition traits. Some degree of genetic differentiation in ADGRG6 between large and small breeds was observed.

Conclusion

We identified one 24.08 Kb region for weight traits and one 42.52 Kb region for thigh-related carcass traits. DRD2 was the major-effect candidate gene for weight traits, and ADGRG6 was the major-effect candidate gene for carcass composition traits. Our results supply essential information for causative mutation identification of carcass traits in broilers.

Introgression of Eastern Chinese and Southern Chinese haplotypes contributes to the improvement of fertility and immunity in European modern pigs

Background

Pigs were domesticated independently from European and Asian wild boars nearly 10,000 years ago. Chinese indigenous pigs have been historically introduced to improve Europe local pigs. However, the geographic origin and biological functions of introgressed Chinese genes in modern European pig breeds remain largely unknown.

Results

Here we explored whole-genome sequencing data from 266 Eurasian wild boars and domestic pigs to produce a fine-scale map of introgression between French Large White (FLW) and Chinese pigs. We show that FLW pigs had historical admixture with both Southern Chinese (SCN) and Eastern Chinese (ECN) pigs ∼200–300 years ago. Moreover, a set of SCN haplotypes was shown to be beneficial for improving disease resistance and ECN haplotypes are favorable for improved reproductive performance in FLW pigs. In addition, we confirm human-mediated introgression events at the AHR locus, at which the haplotype of most likely ECN origin contributes to increased fertility of FLW pigs.

Conclusions

This study advances our understanding of the breeding history of global domestic pigs and highlights the importance of artificial introgression in the formation of phenotypic characteristics in domestic animals.

Genomic diversity revealed by whole-genome sequencing in three Danish commercial pig breeds

Whole-genome sequencing of 217 animals from three Danish commercial pig breeds (Duroc, Landrace [LL], and Yorkshire [YY]) was performed. Twenty-six million single-nucleotide polymorphisms (SNPs) and 8 million insertions or deletions (indels) were uncovered. Among the SNPs, 493,099 variants were located in coding sequences, and 29,430 were predicted to have a high functional impact such as gain or loss of stop codon. Using the whole-genome sequence dataset as the reference, the imputation accuracy for pigs genotyped with high-density SNP chips was examined. The overall average imputation accuracy for all biallelic variants (SNP and indel) was 0.69, while it was 0.83 for variants with minor allele frequency > 0.1. This study provides whole-genome reference data to impute SNP chip-genotyped animals for further studies to fine map quantitative trait loci as well as improving the prediction accuracy in genomic selection. Signatures of selection were identified both through analyses of fixation and differentiation to reveal selective sweeps that may have had prominent roles during breed development or subsequent divergent selection. However, the fixation indices did not indicate a strong divergence among these three breeds. In LL and YY, the integrated haplotype score identified genomic regions under recent selection. These regions contained genes for olfactory receptors and oxidoreductases. Olfactory receptor genes that might have played a major role in the domestication were previously reported to have been under selection in several species including cattle and swine.

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.

New insights into the associations among feed efficiency, metabolizable efficiency traits and related QTL regions in broiler chickens

Background

Improving the feed efficiency would increase profitability for producers while also reducing the environmental footprint of livestock production. This study was conducted to investigate the relationships among feed efficiency traits and metabolizable efficiency traits in 180 male broilers. Significant loci and genes affecting the metabolizable efficiency traits were explored with an imputation-based genome-wide association study. The traits measured or calculated comprised three growth traits, five feed efficiency related traits, and nine metabolizable efficiency traits.

Results

The residual feed intake (RFI) showed moderate to high and positive phenotypic correlations with eight other traits measured, including average daily feed intake (ADFI), dry excreta weight (DEW), gross energy excretion (GEE), crude protein excretion (CPE), metabolizable dry matter (MDM), nitrogen corrected apparent metabolizable energy (AMEn), abdominal fat weight (AbF), and percentage of abdominal fat (AbP). Greater correlations were observed between growth traits and the feed conversion ratio (FCR) than RFI. In addition, the RFI, FCR, ADFI, DEW, GEE, CPE, MDM, AMEn, AbF, and AbP were lower in low-RFI birds than high-RFI birds (P < 0.01 or P < 0.05), whereas the coefficients of MDM and MCP of low-RFI birds were greater than those of high-RFI birds (P < 0.01). Five narrow QTLs for metabolizable efficiency traits were detected, including one 82.46-kb region for DEW and GEE on Gallus gallus chromosome (GGA) 26, one 120.13-kb region for MDM and AMEn on GGA1, one 691.25-kb region for the coefficients of MDM and AMEn on GGA5, one region for the coefficients of MDM and MCP on GGA2 (103.45–103.53 Mb), and one 690.50-kb region for the coefficient of MCP on GGA14. Linkage disequilibrium (LD) analysis indicated that the five regions contained high LD blocks, as well as the genes chromosome 26 C6orf106 homolog (C26H6orf106), LOC396098, SH3 and multiple ankyrin repeat domains 2 (SHANK2), ETS homologous factor (EHF), and histamine receptor H3-like (HRH3L), which are known to be involved in the regulation of neurodevelopment, cell proliferation and differentiation, and food intake.

Conclusions

Selection for low RFI significantly decreased chicken feed intake, excreta output, and abdominal fat deposition, and increased nutrient digestibility without changing the weight gain. Five novel QTL regions involved in the control of metabolizable efficiency in chickens were identified. These results, combined through nutritional and genetic approaches, should facilitate novel insights into improving feed efficiency in poultry and other species.

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.

Landscape of Loci and Candidate Genes for Muscle Fatty Acid Composition in Pigs Revealed by Multiple Population Association Analysis

Genome wide association analyses in diverse populations can identify complex trait loci that are specifically present in one population or shared across multiple populations, which help to better understand the genetic architecture of complex traits in a broader genetic context. In this study, we conducted genome-wide association studies and meta-analysis for 38 fatty acid composition traits with 12–19 million imputed genome sequence SNPs in 2446 pigs from six populations, encompassing White Duroc × Erhualian F₂, Sutai, Duroc-Landrace-Yorkshire (DLY) three-way cross, Laiwu, Erhualian, and Bamaxiang pigs that were originally genotyped with 60 K or 1.4 million single nucleotide polymorphism (SNP) chips. The analyses uncovered 285 lead SNPs (P < 5 × 10^-8), among which 78 locate more than 1 Mb to the lead chip SNPs were considered as novel, largely augmented the landscape of loci for porcine muscle fatty acid composition. Meta-analysis enhanced the association significance at loci near FADS2, ABCD2, ELOVL5, ELOVL6, ELOVL7, SCD, and THRSP genes, suggesting possible existence of population shared mutations underlying these loci. Further haplotype analysis at SCD loci identified a shared 3.7 kb haplotype in F₂, Sutai and DLY pigs showing consistent effects of decreasing C18:0 contents in the three populations. In contrast, at FASN loci, we found an Erhualian specific haplotype explaining the population specific association signals in Erhualian pigs. This study refines our understanding on landscape of loci and candidate genes for fatty acid composition traits of pigs.

Rediscover and Refine QTLs for Pig Scrotal Hernia by Increasing a Specially Designed F3 Population and Using Whole-Genome Sequence Imputation Technology

Pig scrotal hernia is one of the most common congenital defects triggered by both genetic and environmental factors, leading to severe economic loss as well as poor animal welfare in the pig industry. Identification and implementation of genomic regions controlling scrotal hernia in breeding is of great appeal to reduce incidences of hernia in pig production. The aim of this study was to identify such regions or molecular markers affecting scrotal hernia in pigs. First of all, we summarized and analyzed the results of some international teams on scrotal hernia and designed a specially population which contains 246 male individuals. We then performed genome-wide association study (GWAS) in this specially designed population using two scenarios, i.e., the target panel data before and after imputation, which contain 42,365 SNPs and 18,756,672 SNPs, respectively. In addition, a series of methods including genetic differentiation analysis, linkage disequilibrium and linkage analysis (LDLA), and haplotype sharing analysis were appropriate to provide for further analysis to identify the potential gene underlying the QTL. The GWAS in this report detected a highly significant region affecting scrotal hernia within a 24.8Mb region (114.1-138.9Mb) on SSC8. And the result of genetic differentiation analysis also showed a strong genetic differentiation signal between 116.1 and 132.7Mb on SSC8. In addition, the QTL interval was refined to 2.99Mb by combining LDLA and genetic differentiation analysis. Finally, two susceptibility haplotypes were identified through haplotype sharing analysis, with one potential causal gene in it. Our study provided deeper insights into the genetic architecture of pig scrotal hernia and contributed to further fine-mapping and characterize haplotype and gene that influence scrotal hernia in pigs.

Genome-wide association and evolutionary analyses reveal the formation of swine facial wrinkles in Chinese Erhualian pigs

Wrinkles are uneven concave-convex folds, ridges or creases in skin. Facial wrinkles appear in head, typically increasing along with aging. However in several Chinese indigenous pigs, such as Erhualian pigs, rich facial wrinkles have been generated during the growth stages as one of their breed characteristics. To investigate the genetic basis underlying the development of swine facial wrinkles, we estimated the folding extent of facial wrinkles in a herd of Erhualian pigs (n=332), and then conducted genome-wide association studies and multi-trait meta-analysis for facial wrinkles using 60K porcine chips. We found that facial wrinkles had high heritability estimates of ~0.7 in Erhualian pigs. Notably, only one genome-wide significant QTL was detected at 34.8 Mb on porcine chromosome 7. The most significant SNP rs80983858 located at the 3255-bp downstream of candidate gene GRM4, and the G allele was of benefit to increase facial wrinkles. Evolutionary and selection analyses suggested that the haplotypes containing G allele were under artificial selection, which was consistent with local animal sacrificial custom praying for longevity. Our findings made important clues for further deciphering the molecular mechanism of swine facial wrinkles formation, and shed light on the research of skin wrinkle development in human or other mammals.