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

BACKGROUND: Increasing resilience is a priority in modern pig breeding. Recent research shows that general resilience can be quantified via variability in longitudinal data. The collection of such longitudinal data on weight, feed intake and feeding behaviour in pigs has been facilitated by the development of technologies such as automated feeding stations. The goal of this study was to investigate resilience traits, which were estimated as deviations from longitudinal weight, feed intake and feeding behaviour data during the finishing phase. A dataset with 324,207 records between the age of 95 and 155 days on 5,939 Piétrain pigs with known pedigree and genomic information was used. We provided guidelines for a rigid quality control of longitudinal body weight data, as we found that outliers can significantly affect results. Gompertz growth curve analysis, linear modelling and trajectory analyses were used for quantifying resilience traits. RESULTS: To our knowledge, this is the first study comparing resilience traits from longitudinal body weight, feed intake and feeding behaviour data in pigs. We demonstrated that the resilience traits are lowly to moderately heritable for deviations in body weight (h2 = 2.9%-20.2%), in feed intake (9.4%-23.3%) and in feeding behaviour (16.2%-28.3%). Additionally, these traits have good predictive abilities in cross-validation analyses. Deviations in individual body weight and feed intake trajectories are highly correlated (rg = 0.78) with low to moderate favourable genetic correlations with feed conversion ratio (rg = 0.39-0.49). Lastly, we showed that some resilience traits, such as the natural logarithm of variances of observed versus predicted body weights (lnvarweight), are more robust to lower observation frequencies and are repeatable over three different time periods of the finishing phase. CONCLUSIONS: Our results will help future studies investigating resilience traits and resilience-related traits. Moreover, our study provides first results on standardization of quality control and efficient data sampling from automated feeding station data. Our findings will be valuable for breeding organizations as they offer evidence that pigs' general resilience can be selected on with good accuracy. Moreover, this methodology might be extended to other species to quantify resilience based on longitudinal data.

Genome-wide association studies for additive and dominance effects for body composition traits in commercial crossbred Piétrain pigs.

Fat depth (FD) and muscle depth (MD) are economically important traits and used to estimate carcass lean content (LMP), which is one of the main breeding objectives in pig breeding programmes. We assessed the genetic architectures of body composition traits for additive and dominance effects in commercial crossbred Piétrain pigs using both 50 K array and sequence genotypes. We first performed a genome-wide association study (GWAS) using single-marker association analysis with a false discovery rate of 0.1. Then, we estimated the additive and dominance effects of the most significant variant in the quantitative trait loci (QTL) regions. It was investigated whether the use of whole-genome sequence (WGS) will improve the QTL detection (both additive and dominance) with a higher power compared with lower density SNP arrays. Our results showed that more QTL regions were detected by WGS compared with 50 K array (n = 54 vs. n = 17). Of the novel associated regions associated with FD and LMP and detected by WGS, the most pronounced peak was on SSC13, situated at ~116-118, 121-127 and 129-134 Mbp. Additionally, we found that only additive effects contributed to the genetic architecture of the analysed traits and no significant dominance effects were found for the tested SNPs at QTL regions, regardless of panel density. The associated SNPs are located in or near several relevant candidate genes. Of these genes, GABRR2, GALR1, RNGTT, CDH20 and MC4R have been previously reported as being associated with fat deposition traits. However, the genes on SSC1 (ZNF292, ORC3, CNR1, SRSF12, MDN1, TSHZ1, RELCH and RNF152) and SSC18 (TTC26 and KIAA1549) have not been reported previously to our best knowledge. Our current findings provide insights into the genomic regions influencing composition traits in Piétrain pigs.

Longitudinal individual predictions from irregular repeated measurements data.

Intensive longitudinal data can be used to explore important associations and patterns between various types of inputs and outcomes. Nonlinear relations and irregular measurement occasions can pose problems to develop an accurate model for these kinds of data. This paper focuses on the development, fitting and evaluation of a prediction model with irregular intensive longitudinal data. A three-step process for developing a prediction tool for (daily) monitoring and prediction is outlined and illustrated for intensive weight measurements in piglets. Step 1 addresses a nonlinear relation in the data by developing and applying a normalizing transformation. Step 2 addresses the intermittent nature of the time points by aligning the measurement times to a common time grid with a broken-stick model. Step 3 addresses the prediction problem by selecting and evaluating inputs and covariates in the model to obtain accurate predictions. The final model predicts future outcomes accurately, while allowing for nonlinearities between input and output and for different measurement histories of individuals. The methodology described can be used to develop a tool to deal with intensive irregular longitudinal data that uses the available information in an optimal way. The resulting tool demonstrated to perform well for piglet weight prediction and can be adapted to many different applications.

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.

CRISPR/Cas9-editing of KISS1 to generate pigs with hypogonadotropic hypogonadism as a castration free trait.

Introduction: Most male pigs are surgically castrated to avoid puberty-derived boar taint and aggressiveness. However, this surgical intervention represents a welfare concern in swine production. Disrupting porcine KISS1 is hypothesized to delay or abolish puberty by inducing variable hypogonadotropism and thus preventing the need for castration. Methods: To test this hypothesis, we generated the first KISS1-edited large animal using CRISPR/Cas9-ribonucleoproteins and single-stranded donor oligonucleotides. The targeted region preceded the sequence encoding a conserved core motif of kisspeptin. Genome editors were intracytoplasmically injected into 684 swine zygotes and transferred to 19 hormonally synchronized surrogate sows. In nine litters, 49 American Yorkshire and 20 Duroc liveborn piglets were naturally farrowed. Results: Thirty-five of these pigs bore KISS1-disruptive alleles ranging in frequency from 5% to 97% and did not phenotypically differ from their wild-type counterparts. In contrast, four KISS1-edited pigs (two boars and two gilts) with disruptive allele frequencies of 96% and 100% demonstrated full hypogonadotropism, infantile reproductive tracts, and failed to reach sexual maturity. Change in body weight during development was unaffected by editing KISS1. Founder pigs partially carrying KISS1-disruptive alleles were bred resulting in a total of 53 KISS1 +/+, 60 KISS1 +/-, and 34 KISS1 -/- F1 liveborn piglets, confirming germline transmission. Discussion: Results demonstrate that a high proportion of KISS1 alleles in pigs must be disrupted before variation in gonadotropin secretion is observed, suggesting that even a small amount of kisspeptin ligand is sufficient to confer proper sexual development and puberty in pigs. Follow-on studies will evaluate fertility restoration in KISS1 KO breeding stock to fully realize the potential of KISS1 gene edits to eliminate the need for surgical castration.

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.

Identification of protein-damaging mutations in 10 swine taste receptors and 191 appetite-reward genes.

BACKGROUND: Taste receptors (TASRs) are essential for the body's recognition of chemical compounds. In the tongue, TASRs sense the sweet and umami and the toxin-related bitter taste thus promoting a particular eating behaviour. Moreover, their relevance in other organs is now becoming evident. In the intestine, they regulate nutrient absorption and gut motility. Upon ligand binding, TASRs activate the appetite-reward circuitry to signal the nervous system and keep body homeostasis. With the aim to identify genetic variation in the swine TASRs and in the genes from the appetite and the reward pathways, we have sequenced the exons of 201 TASRs and appetite-reward genes from 304 pigs belonging to ten breeds, wild boars and to two phenotypically extreme groups from a F2 resource with data on growth and fat deposition. RESULTS: We identified 2,766 coding variants 395 of which were predicted to have a strong impact on protein sequence and function. 334 variants were present in only one breed and at predicted alternative allele frequency (pAAF) ≥ 0.1. The Asian pigs and the wild boars showed the largest proportion of breed specific variants. We also compared the pAAF of the two F2 groups and found that variants in TAS2R39 and CD36 display significant differences suggesting that these genes could influence growth and fat deposition. We developed a 128-variant genotyping assay and confirmed 57 of these variants. CONCLUSIONS: We have identified thousands of variants affecting TASRs as well as genes involved in the appetite and the reward mechanisms. Some of these genes have been already associated to taste preferences, appetite or behaviour in humans and mouse. We have also detected indications of a potential relationship of some of these genes with growth and fat deposition, which could have been caused by changes in taste preferences, appetite or reward and ultimately impact on food intake. A genotyping array with 57 variants in 31 of these genes is now available for genotyping and start elucidating the impact of genetic variation in these genes on pig biology and breeding.

Identification of genetic variation in the swine toll-like receptors and development of a porcine TLR genotyping array.

BACKGROUND: Toll-like receptors (TLR) are crucial in innate immunity for the recognition of a broad range of microbial pathogens and are expressed in multiple cell types. There are 10 TLR genes described in the pig genome. RESULTS: With a twofold objective i.e. to catalogue genetic variants in porcine TLR genes and develop a genotyping array for genetic association studies on immune-related traits, we combined targeted sub-genome enrichment and high-throughput sequencing to sequence the 10 porcine TLR genes in 266 pigs from 10 breeds and wild boars using a DNA-pooling strategy. We identified 306 single nucleotide variants across the 10 TLR and 11 populations, 87 of which were novel. One hundred and forty-seven positions i.e. six stop-gains and 141 non-synonymous substitutions were predicted to alter the protein sequence. Three positions were unique to a single breed with alternative allele frequencies equal to or higher than 0.5. We designed a genotyping array for future applications in genetic association studies, with a selection of 126 variants based on their predicted impact on protein sequence. Since TLR4, TLR7 and TLR9 were underrepresented in this selection, we also included three variants that were located in the 3'UTR of these genes. We tested the array by genotyping 214 of the 266 sequenced pigs. We found that 93 variants that involved the 10 TLR genes were polymorphic in these animals. Twelve of these variants were novel. Furthermore, seven known variants that are associated with immune-related phenotypes are present on the array and can thus be used to test such associations in additional populations. CONCLUSIONS: We identified genetic variations that potentially have an impact on the protein sequence of porcine TLR. A genotyping array with 80 non-synonymous, 10 synonymous and three 3'UTR polymorphisms in the 10 TLR genes is now available for association studies in swine populations with measures on immune-related traits.

Genetic and phenotypic correlations between performance traits with meat quality and carcass characteristics in commercial crossbred pigs.

Genetic correlations between performance traits with meat quality and carcass traits were estimated on 6,408 commercial crossbred pigs with performance traits recorded in production systems with 2,100 of them having meat quality and carcass measurements. Significant fixed effects (company, sex and batch), covariates (birth weight, cold carcass weight, and age), random effects (additive, litter and maternal) were fitted in the statistical models. A series of pairwise bivariate analyses were implemented in ASREML to estimate heritability, phenotypic, and genetic correlations between performance traits (n = 9) with meat quality (n = 25) and carcass (n = 19) traits. The animals had a pedigree compromised of 9,439 animals over 15 generations. Performance traits had low-to-moderate heritabilities (±SE), ranged from 0.07±0.13 to 0.45±0.07 for weaning weight, and ultrasound backfat depth, respectively. Genetic correlations between performance and carcass traits were moderate to high. The results indicate that: (a) selection for birth weight may increase drip loss, lightness of longissimus dorsi, and gluteus medius muscles but may reduce fat depth; (b) selection for nursery weight can be valuable for increasing both quantity and quality traits; (c) selection for increased daily gain may increase the carcass weight and most of the primal cuts. These findings suggest that deterioration of pork quality may have occurred over many generations through the selection for less backfat thickness, and feed efficiency, but selection for growth had no adverse effects on pork quality. Low-to-moderate heritabilities for performance traits indicate that they could be improved using traditional selection or genomic selection. The estimated genetic parameters for performance, carcass and meat quality traits may be incorporated into the breeding programs that emphasize product quality in these Canadian swine populations.