Large-scale association study for structural soundness and leg locomotion traits in the pig

Deep learning based landmark detection for measuring hock and knee angles in sows

This paper presents a visual deep learning approach to automatically determine hock and knee angles from sow images. Lameness is the second largest reason for culling of breeding herd females and relies on human observers to provide visual scoring for detection which can be slow, subjective, and inconsistent. A deep learning model classified and detected ten and two key body landmarks from the side and rear profile images, respectively (mean average precision = 0.94). Trigonometric-based formulae were derived to calculate hock and knee angles using the features extracted from the imagery. Automated angle measurements were compared with manual results from each image (average root mean square error [RMSE] = 4.13°), where all correlation slopes (average R 2 = 0.84) were statistically different from zero (P < 0.05); all automated measurements were in statistical agreement with manually collected measurements using the Bland-Altman procedure. This approach will be of interest to animal geneticists, scientists, and practitioners for obtaining objective angle measurements that can be factored into gilt replacement criteria to optimize sow breeding units.

Assessing Genetic Diversity and Searching for Selection Signatures by Comparison between the Indigenous Livni and Duroc Breeds in Local Livestock of the Central Region of Russia

Indigenous pig breeds are mainly associated with the adaptive capacity that is necessary to respond adequately to climate change, food security, and livelihood needs, and natural resources conservation. Livni pigs are an indigenous fat-type breed farmed in a single farm in the Orel region and located in the Central European part of the Russian Federation. To determine the genomic regions and genes that are affected by artificial selection, we conducted the comparative study of two pig breeds with different breeding histories and breeding objectives, i.e., the native fat-type Livni and meat-type Duroc breeds using the Porcine GGP HD BeadChip, which contains ~80,000 SNPs. To check the Livni pigs for possible admixture, the Landrace and the Large White breeds were included into the study of genetic diversity as these breeds participated in the formation of the Livni pigs. We observed the highest level of genetic diversity in Livni pigs compared to commercial breeds (UHE = 0.409 vs. 0.319–0.359, p < 0.001; AR = 1.995 vs. 1.894–1.964, p < 0.001). A slight excess of heterozygotes was found in all of the breeds. We identified 291 candidate genes, which were localized within the regions under putative selection, including 22 and 228 genes, which were specific for Livni and Duroc breeds, respectively, and 41 genes common for both breeds. A detailed analysis of the molecular functions identified the genes, which were related to the formation of meat and fat traits, and adaptation to environmental stress, including extreme temperatures, which were different between breeds. Our research results are useful for conservation and sustainable breeding of Livni breed, which shows a high level of genetic diversity. This makes Livni one of the valuable national pig genetic resources.

Genetic diversity and population structure of six autochthonous pig breeds from Croatia, Serbia, and Slovenia

BACKGROUND

The importance of local breeds as genetic reservoirs of valuable genetic variation is well established. Pig breeding in Central and South-Eastern Europe has a long tradition that led to the formation of several local pig breeds. In the present study, genetic diversity parameters were analysed in six autochthonous pig breeds from Slovenia, Croatia and Serbia (Banija spotted, Black Slavonian, Turopolje pig, Swallow-bellied Mangalitsa, Moravka and Krskopolje pig). Animals from each of these breeds were genotyped using microsatellites and single nucleotide polymorphisms (SNPs). The results obtained with these two marker systems and those based on pedigree data were compared. In addition, we estimated inbreeding levels based on the distribution of runs of homozygosity (ROH) and identified genomic regions under selection pressure using ROH islands and the integrated haplotype score (iHS).

RESULTS

The lowest heterozygosity values calculated from microsatellite and SNP data were observed in the Turopolje pig. The observed heterozygosity was higher than the expected heterozygosity in the Black Slavonian, Moravka and Turopolje pig. Both types of markers allowed us to distinguish clusters of individuals belonging to each breed. The analysis of admixture between breeds revealed potential gene flow between the Mangalitsa and Moravka, and between the Mangalitsa and Black Slavonian, but no introgression events were detected in the Banija spotted and Turopolje pig. The distribution of ROH across the genome was not uniform. Analysis of the ROH islands identified genomic regions with an extremely high frequency of shared ROH within the Swallow-bellied Mangalitsa, which harboured genes associated with cholesterol biosynthesis, fatty acid metabolism and daily weight gain. The iHS approach to detect signatures of selection revealed candidate regions containing genes with potential roles in reproduction traits and disease resistance.

CONCLUSIONS

Based on the estimation of population parameters obtained from three data sets, we showed the existence of relationships among the six pig breeds analysed here. Analysis of the distribution of ROH allowed us to estimate the level of inbreeding and the extent of homozygous regions in these breeds. The iHS analysis revealed genomic regions potentially associated with phenotypic traits and allowed the detection of genomic regions under selection pressure.

An association study on imputed whole-genome resequencing from high-throughput sequencing data for body traits in crossbred pigs

Body traits are important economic factors in the pig industry. Genome-wide association studies (GWASs) have been widely applied using high-density genotype data to detect QTL in pigs. The aim of the present study was to detect the genetic variants significantly associated with body traits in crossbred pigs using the Illumina Porcine SNP50 BeadChip and imputed whole-genome sequence data. A set of seven body traits - body length, body height, chest circumference, cannon bone circumference, leg buttock circumference, back fat thickness and loin muscle depth - were measured. Moderate to high heritabilities were obtained for most traits (from 0.14 to 0.46), and significant genetic and phenotypic correlations among them were observed. GWAS identified 714 significantly associated SNPs located at 39 regions on all autosomes for body traits, and a total of seven functionally related candidate genes: PIK3CD, HOXA, PCGF2, CHST11, COL2A1, BMI1 and OSR2. Functional enrichment analysis revealed that candidate genes were enriched in the estrogen signaling pathway, embryonic skeletal system morphogenesis and embryonic skeletal system development. These results aim to uncover the genetic mechanisms underlying body development and marker-assisted selection programs focusing on body traits in pigs.

Genome-Wide Detection of Genetic Loci and Candidate Genes for Body Conformation Traits in Duroc × Landrace × Yorkshire Crossbred Pigs

The Duroc × (Landrace × Yorkshire) hybrid pigs (DLY) are the most popular commercial pigs, providing consumers with the largest source of pork. In order to gain more insights into the genetic architecture of economically important traits in pigs, we performed a genome-wide association study (GWAS) using the GeneSeek Porcine 50 K SNP Chip to map the genetic markers and genes associated with body conformation traits (BCT) in 311 DLY pigs. The quantitative traits analyzed included body weight (BW), carcass length (CL), body length (BL), body height (BH), and body mass index (BMI). BMI was defined as BMI_CL, BMI_BL, and BMI_BH, respectively, based on CL, BL, and BH phenotypic data. We identified 82 SNPs for the seven traits by GEMMA-based and FarmCPU-based GWASs. Both methods detected two quantitative trait loci (QTL) on SSC8 and SSC17 for body conformation traits. Several candidate genes (such as TNFAIP3, KDM4C, HSPG2, BMP2, PLCB4, and GRM5) were found to be associated with body weight and body conformation traits in pigs. Notably, the BMP2 gene had pleiotropic effects on CL, BL, BH, BMI_CL, and BMI_BL and is proposed as a strong candidate gene for body size due to its involvement in growth and bone development. Furthermore, gene set enrichment analysis indicated that most of the pathway terms are associated with regulation of cell growth, negative regulation of cell population proliferation, and chondrocyte differentiation. We anticipate that these results further advance our understanding of the genetic architecture of body conformation traits in the popular commercial DLY pigs and provide new insights into the genetic architecture of BMI in pigs.

A Genome-Wide Association Study on Feed Efficiency Related Traits in Landrace Pigs

Feed efficiency (FE) traits in pigs are of utmost economic importance. Genetic improvement of FE related traits in pigs might significantly reduce production cost and energy consumption. Hence, our study aimed at identifying SNPs and candidate genes associated with FE related traits, including feed conversion ratio (FCR), average daily gain (ADG), average daily feed intake (ADFI), and residual feed intake (RFI). A genome-wide association study (GWAS) was performed for the four FE related traits in 296 Landrace pigs genotyped with PorcineSNP50 BeadChip. Two different single-trait methods, single SNP linear model GWAS (LM-GWAS) and single-step GWAS (ssGWAS), were implemented. Our results showed that the two methods showed high consistency with respect to SNP identification. A total of 32 common significant SNPs associated with the four FE related traits were identified. Bioinformatics analysis revealed eight common QTL regions, of which three QTL regions related to ADFI and RFI traits were overlapped. Gene ontology analysis revealed six common candidate genes (PRELID2, GPER1, PDX1, TEX2, PLCL2, ICAM2) relevant for the four FE related traits. These genes are involved in the processes of fat synthesis and decomposition, lipid transport process, insulin metabolism, among others. Our results provide, new insights into the genetic mechanisms and candidate function genes of FE related traits in pigs. However, further investigations to validate these results are warranted.

Genome-wide association study of bone mineral density trait among three pig breeds

Leg weakness (LW) issues are a great concern for pig breeding industry. And it also has a serious impact on animal welfare. To dissect the genetic architecture of limb-and-hoof firmness in commercial pigs, a genome-wide association study was conducted on bone mineral density (BMD) in three sow populations, including Duroc, Landrace and Yorkshire. The BMD data were obtained by ultrasound technology from 812 pigs (including Duroc 115, Landrace 243 and Yorkshire 454). In addition, all pigs were genotyped using genome-by-sequencing and a total of 224 162 single-nucleotide polymorphisms (SNPs) were obtained. After quality control, 218 141 SNPs were used for subsequent genome-wide association analysis. Nine significant associations were identified on chromosomes 3, 5, 6, 7, 9, 10, 12 and 18 that passed Bonferroni correction threshold of 0.05/(total SNP numbers). The most significant locus that associated with BMD (P value = 1.92e-14) was detected at approximately 41.7 Mb on SSC6 (SSC stands for Sus scrofa chromosome). CUL7, PTK7, SRF, VEGFA, RHEB, PRKAR1A and TPO that are located near the lead SNP of significant loci were highlighted as functionally plausible candidate genes for sow limb-and-hoof firmness. Moreover, we also applied a new method to measure the BMD data of pigs by ultrasound technology. The results provide an insight into the genetic architecture of LW and can also help to improve animal welfare in pigs.

Genome-wide comparative analysis ofbone morphogenetic proteins: genomic structure, phylogeny, and expression patterns in the golden pompano,Trachinotus ovatus(Linnaeus, 1758)

Bone morphogenetic proteins (BMPs) play important roles in various physiological processes, especially during the formation and maintenance of various organs. In this study, we first obtained and characterized twenty BMP genes from the Trachinotus ovatus genome (designated as ToBMPs). Sequence alignment and phylogenetic analysis both indicated that the predicted amino acid sequences of ToBMP were highly conserved with corresponding homologs of other species. Moreover, a comparative analysis was performed with seven representative vertebrate genomes and found difference in number of BMP3 genes in different species, which three members, BMP3a, BMP3b-1, and BMP3b-2, existed in diploid T. ovatus, but there were four and two members in tetraploidized Cyprinus carpio (BMP3a-1, BMP3a-2, BMP3b-1, and BMP3b-2) and diploid Danio rerio (BMP3a and BMP3b), respectively. The amino acid alignment and genomic structure analysis of ToBMP3 also suggested that the BMP3 gene had expanded in T. ovatus. Furthermore, tissue expression patterns were assessed for the small intestine, liver, white muscle, brain, spleen, fin, gill, head kidney, stomach, blood, and gonads. It was discovered that BMP1, BMP2, BMP3a, BMP4, BMP6, BMP7b, BMP11, and BMP16 were ubiquitously expressed in all the tissues tested. To study the regulatory function of BMP in response to the intake of different types of food, the expression changes in BMP mRNAs were detected by qRT-PCR, and the results showed that the majority of the BMP genes had the highest mRNA levels in the small intestine and liver after ingesting pelleted feed. Our data provide a useful resource for further studies on how paralogous genes may have different expression profiles in T. ovatus.

Evaluation of genotype quality parameters for SowPro90, a new genotyping array for swine1

Understanding early predictors of sow fertility has the potential to improve genomic predictions. A custom SNP array (SowPro90 produced by Affymetrix) was developed to include genetic variants overlapping quantitative trait loci for age at puberty, one of the earliest indicators of sow fertility, as well as variants related to innate and adaptive immunity. The polymorphisms included in the custom genotyping array were identified using multiple genomic approaches including deep genomic and transcriptomic sequencing and genome-wide associations. Animals from research and commercial populations (n = 2,586) were genotyped for 103,476 SNPs included in SowPro90. To assess the quality of data generated, genotype concordance was evaluated between the SowPro90 and Porcine SNP60 BeadArray using a subset of common SNP (n = 44,708) and animals (n = 277). The mean genotype concordance rate per SNP was 98.4%. Differences in distribution of data quality were observed between the platforms indicating the need for platform specific thresholds for quality parameters. The optimal thresholds for SowPro90 (≥97% SNP and ≥93% sample call rate) were obtained by analyzing the data quality distribution and genotype concordance per SNP across platforms. At ≥97% SNP call rate, there were 42,151 SNPs (94.3%) retained with a mean genotype concordance of 98.6% across platforms. Similarly, ≥94% SNPs and ≥85% sample call rates were established as thresholds for Porcine SNP60 BeadArray. At ≥94% SNPs call rate, there were 41,043 SNPs (91.8%) retained with a mean genotype concordance of 98.6% across platforms. Final evaluation of SowPro90 array content (n = 103,476) at ≥97% SNPs and ≥93% sample call rates allowed retention of 89,040 SNPs (86%) for downstream analysis. The findings and strategy for quality control could be helpful in identifying consistent, high-quality genotypes for genomic evaluations, especially when integrating genotype data from different platforms.

New world goat populations are a genetically diverse reservoir for future use

Western hemisphere goats have European, African and Central Asian origins, and some local or rare breeds are reported to be adapted to their environments and economically important. By-in-large these genetic resources have not been quantified. Using 50 K SNP genotypes of 244 animals from 12 goat populations in United States, Costa Rica, Brazil and Argentina, we evaluated the genetic diversity, population structure and selective sweeps documenting goat migration to the "New World". Our findings suggest the concept of breed, particularly among "locally adapted" breeds, is not a meaningful way to characterize goat populations. The USA Spanish goats were found to be an important genetic reservoir, sharing genomic composition with the wild ancestor and with specialized breeds (e.g. Angora, Lamancha and Saanen). Results suggest goats in the Americas have substantial genetic diversity to use in selection and promote environmental adaptation or product driven specialization. These findings highlight the importance of maintaining goat conservation programs and suggest an awaiting reservoir of genetic diversity for breeding and research while simultaneously discarding concerns about breed designations.

Objective evaluation of female feet and leg joint conformation at time of selection and post first parity in swine1

Feet and legs of replacement females were objectively evaluated at selection, i.e., approximately 150 d of age (n = 319) and post first parity, i.e., any time after weaning of first litter and before second parturition (n = 277) to 1) compare feet and leg joint angle ranges between selection and post first parity; 2) identify feet and leg joint angle differences between selection and first 3 wk of second gestation; 3) identify feet and leg joint angle differences between farms and gestation days during second gestation; and 4) obtain genetic variance components for conformation angles for the two time points measured. Angles for carpal joint (knee), metacarpophalangeal joint (front pastern), metatarsophalangeal joint (rear pastern), tarsal joint (hock), and rear stance were measured using image analysis software. Between selection and post first parity, significant differences were observed for all joints measured (P < 0.05). Knee, front and rear pastern angles were less (more flexion), and hock angles were greater (less flexion) as age progressed (P < 0.05), while the rear stance pattern was less (feet further under center) at selection than post first parity (only including measures during first 3 wk of second gestation). Only using post first parity leg conformation information, farm was a significant source of variation for front and rear pasterns and rear stance angle measurements (P < 0.05). Knee angle was less (more flexion; P < 0.05) as gestation age progressed. Heritability estimates were low to moderate (0.04-0.35) for all traits measured across time points. Genetic correlations between the same joints at different time points were high (>0.8) between the front leg joints and low (<0.2) between the rear leg joints. High genetic correlations between time points indicate that the trait can be considered the same at either time point, and low genetic correlations indicate that the trait at different time points should be considered as two separate traits. Minimal change in the front leg suggests conformation traits that remain between selection and post first parity, while larger changes in rear leg indicate that rear leg conformation traits should be evaluated at multiple time periods.

Genetic variants of major genes contributing to phosphate and calcium homeostasis and their association with serum parameters in pigs

Calcium and phosphorus are irreplaceable components of life. Tracking the fate of calcium and phosphorus in organisms deserves high attention due to their relevance in bone metabolism and subsequently animal health. Indeed, bone serves as reservoir for calcium and phosphorus, whose formation and resorption follow specific molecular routes including hormones, receptors, and transcription factors. The objective of the study was to analyze the genetic variation of major components driving mineral utilization such as calcitonin receptor, calcium sensing receptor, fibroblast growth factor 23 (FGF23), parathyroid hormone receptor, osteopontin, stanniocalcin 1, RAF-type zinc finger domain containing 1 (TRAFD1), and vitamin D receptor. A German Landrace pig population (n = 360) was used to perform an association analysis between selected single nucleotide polymorphisms (SNP) and relevant serum parameters (calcium, phosphorus, calcium/phosphorus ratio, alkaline phosphatase). Analyzed SNPs in FGF23 (rs710498025) and TRAFD1 (rs345195312) were significantly (p ≤ 0.05) associated with the serum calcium/phosphorus ratio and serum phosphorus levels, respectively. This might represent a modulation of the homeostatic balance between calcium and phosphorus. Furthermore, TRAFD1 is known to be involved in skeletal disorders which emphasize its link to phosphorus utilization and immune system. However, none of the analyzed genetic variants of these major regulators of phosphate and calcium homeostasis showed significant associations after correction for multiple testing (q value > 0.05). Thus, minor contributors as well as unknown and yet to be elucidated regulators of mineral homeostasis need to be characterized towards the implementation of improved phosphorus efficiency in pig breeding programs.

Characterization and spatiotemporal expression analysis of nine bone morphogenetic protein family genes during intermuscular bone development in blunt snout bream

Intermuscular bones (IBs) only exist in the myosepta of lower teleosts and its molecular mechanism remains to be clarified. Bone morphogenetic proteins (BMPs) have been demonstrated to be involved in various physiological processes, including bone and cartilage formation. In this study, we firstly obtained and characterized nine bmp genes for Megalobrama amblycephala, which belongs to Cyprinidae and have a certain amount of IBs. Sequence alignment and phylogenetic analysis both documented that the mature proteins of M. amblycephala bmp genes were highly conserved with other corresponding homologs, respectively, indicating that the function of each bmp gene has been conserved throughout evolution. As a step to characterize potential involvement of bmp genes in IB formation and development, spatiotemporal expressions of nine bmp genes (bmp2a, bmp2b, bmp3, bmp4, bmp5, bmp7b, bmp8a, bmp14 and bmp16) were investigated during the key development stages of IBs. During the ossification process from stage I (the IBs haven't emerged) to stage IV (all of the IBs ossified in the tail with the mature morphology), the expression profiles revealed that bmp16 was the most abundant transcript while bmp4 had the lowest abundance. The mRNA levels of bmp3, bmp4, bmp5 and bmp8a increased significantly at stage II, suggesting their roles in stimulating IB formation. The expression of bmp7b reached the highest level at stage III (the rapid period of IB development), suggesting potential involvement of bmp7b in promoting osteoblast differentiation. With the exception of bmp7b and bmp16, most bmp genes appeared a significant increase at IB maturation phase (stage IV), which means that they may play important roles in maintenance of IB morphogenesis. Spatial tissue distribution of bmp genes showed that most bmp genes were observed at the highest level in developing IBs at one year old fish. Spatiotemporal expression patterns suggest the potential key roles of these bmp genes in IBs formation and maintenance in fish, being as possible promoters or inhibitors.

Effect of rubber flooring on group-housed sows' gait and claw and skin lesions

This study evaluated the influence of floor type on sow welfare in terms of lameness, claw lesions, and skin lesions. In a 2 × 3 factorial design, we have investigated the effect of rubber coverings on concrete floors and the effect of 3 levels of dietary zinc supplementation on locomotion and claw and skin lesions in group-housed sows. Six groups of 21 ± 4 hybrid sows were monitored during 3 successive reproductive cycles. The sows were group housed from d 28 after insemination (d 0) until 1 wk before expected farrowing date (d 108) in pens with either exposed concrete floors or concrete floors covered with rubber in part of the lying area and the fully slatted area. During each reproductive cycle, locomotion and skin lesions were assessed 4 times (d 28, 50, 108, and 140) and claw lesions were assessed twice (d 50 and 140). Results are given as least squares means ± SE. Locomotion and claw scores were given in millimeters, on analog scales of 150 and 160 mm, respectively. Here, we report on the effect of floor type, which did not interact with dietary zinc concentration ( > 0.10 for all variables). At move to group (d 28) and mid gestation (d 50), no differences between floor treatments were seen in locomotion ( > 0.10). At the end of gestation (d 108), sows housed on rubber flooring scored 9.9 ± 4.1 mm better on gait ( < 0.001). Regarding claw disorders, both parameters "heel overgrowth and erosion" (difference of 4.6 ± 1.8 mm; = 0.01) and "heel-sole crack" (difference of 3.1 ± 1.5 mm; = 0.04) scores were better for sows on rubber flooring at mid gestation (d 50). However, sows on rubber flooring scored worse for "vertical cracks in the wall horn" (difference of 3.4 ± 1.7 mm; = 0.04). At the end of lactation (d 140), both "white line" (difference of 2.9 ± 1 mm; = 0.02) and "claw length" (difference of 4.7 ± 1.4 mm; < 0.001) had better scores on rubber flooring. No differences for skin lesions were observed between floor treatments. The improved scores for gait toward the end of gestation and some types of claw disorders at mid gestation suggest that rubber flooring in group housing has a beneficial effect on the overall leg health of sows. The documented increase in vertical cracks in the wall horn at d 50 requires further investigation.

Locomotion Disorders and Skin and Claw Lesions in Gestating Sows Housed in Dynamic versus Static Groups

Lameness and lesions to the skin and claws of sows in group housing are commonly occurring indicators of reduced welfare. Typically, these problems are more common in group housing than in individual housing systems. Group management type (dynamic versus static) and stage of gestation influence the behavior of the animals, which in turn influences the occurrence of these problems. The present study compared prevalence, incidence and mean scores of lameness and skin and claw lesions in static versus dynamic group housed sows at different stages of gestation during three consecutive reproductive cycles. A total of 10 Belgian sow herds were monitored; 5 in which dynamic groups and 5 in which static groups were utilized. All sows were visually assessed for lameness and skin lesions three times per cycle and the claws of the hind limbs were assessed once per cycle. Lameness and claw lesions were assessed using visual analogue scales. Static groups, in comparison with dynamic groups, demonstrated lower lameness scores (P<0.05) and decreased skin lesion prevalence (24.9 vs. 47.3%, P<0.05) at the end of gestation. There was no difference between treatment group regarding claw lesion prevalence with 75.5% of sows demonstrating claw lesions regardless of group management. Prevalences of lameness (22.4 vs. 8.9%, P<0.05) and skin lesions (46.6 vs. 4.4%, P<0.05) were highest during the group-housed phase compared to the individually housed phases. Although the prevalence of lameness and skin lesions did not differ three days after grouping versus at the end of the group-housing phase, their incidence peaked during the first three days after moving from the insemination stalls to the group. In conclusion, the first three days after grouping was the most risky period for lameness incidence, but there was no significant difference between static or dynamic group management.

Genetic trends of conformation traits and genetic correlations to osteochondrosis in boars

The objective of our study was to investigate the heritabilities and genetic correlations between traits from a linear exterior assessment system and osteochondrosis (OC) measured by computed tomography (CT), and in addition, to study the genetic trend in a population where the conformation traits have been included in the breeding goal. The data material consisted of phenotypes from a total of 4571 Norsvin Landrace test boars. At the end of the test period, all boars were subjected to a detailed exterior assessment system. Within 10 days of the assessment, the boars were CT scanned for measuring OC. The total score of osteochondrosis (OCT), used in this study, is the sum of phenotypes from the assessment on the medial and lateral condyles at the distal end of both the humerus and the femur of the right and the left leg of the boar based on images from CT. The exterior assessment traits included in the study were; 'front leg knee' (FKNE), 'front leg pasterns' (FPAS), 'front leg stance' (FSTA), 'front leg twisted pasterns' (FFLK), 'hind leg stance', 'hind leg pasterns' (HPAS), 'hind leg standing under' (HSTU), 'hind leg small inner toe', 'dipped back', 'arched back' (ARCH) and 'waddling hindquarters' (WADL). The estimation of (co)variance components and breeding values were performed using bivariate animal genetic models. Breeding values for HSTU, HPAS, FPAS, WADL and OCT traits were additional outputs from the same bivariate analyses. The lowest heritability was found for FFLK (h 2 FFLK=0.05), whereas FPAS was estimated to have the highest heritability (h 2 FPAS=0.36), and OCT demonstrating a heritability of 0.29. Significant genetic correlations were found between several traits; the strongest correlation was between FSTA and FFLK (0.94), which was followed by the correlation between FPAS and FKNE (0.69). The traits ARCH and FSTA had significant genetic correlations to OCT, whereas all other genetic correlations between OCT and the conformation traits were low and not significantly different from 0. Our study shows positive genetic trends for the conformation traits included in the breeding goal. In general, low genetic correlations between conformation traits and OC were observed in our study.

Longitudinal gait development and variability of growing pigs reared on three different floor types

Biomechanical investigation into locomotor pathology in commercial pigs is lacking despite this being a major concern for the industry. Different floor types are used in modern, intensive pig production systems at different stages of the pigs' production cycle. The general perception holds that slatted and/or hard solid concrete surfaces are inferior to soft straw-covered floors regarding healthy musculoskeletal development. Previous studies have compared pigs housed on different floor types using clinical, subjective assessment of leg weakness and lameness. However, reliability studies generally report a low repeatability of clinical lameness scoring. The objective of this study was to quantitatively assess the long-term effect of pen floors, reflected in the biomechanical gait characteristics and associated welfare of the pigs. A cohort of 24 pigs housed on one of three different floor types was followed from 37 to 90 kg average liveweight, with gait analysis (motion capture) starting at 63 kg. The three floor types were fully slatted concrete, partly slatted concrete and deep straw-bedded surfaces, all located within the same building. Pigs underwent five repeated camera-based motion captures, 7 to 10 days apart, during which 3D coordinate data of reflective skin markers attached to leg anatomical landmarks were collected. Pigs walked on the same solid concrete walkway during captures. One-way ANOVA and repeated measures ANOVA were used to analyse the gait data. Results revealed changes over time in the spatiotemporal gait pattern which were similar in magnitude and direction for the pigs from different floor types. Significant increases in elbow joint flexion with age were observed in all pigs (P⩽0.050; +6°). There were few differences between floor groups, except for the step-to-stride ratio in the hind legs being more irregular in pigs housed on partly slatted floors (P=0.012; 3.6 times higher s.d.) compared with those on 5 to 10 cm straw-bedding in all pen areas. As the level of clinical problems was generally low in this cohort, it may be that floors elicit problems only when there is a primary predisposing factor increasing weakness in susceptible tissues.

Genome wide analysis reveals single nucleotide polymorphisms associated with fatness and putative novel copy number variants in three pig breeds

BACKGROUND

Obesity, excess fat tissue in the body, can underlie a variety of medical complaints including heart disease, stroke and cancer. The pig is an excellent model organism for the study of various human disorders, including obesity, as well as being the foremost agricultural species. In order to identify genetic variants associated with fatness, we used a selective genomic approach sampling DNA from animals at the extreme ends of the fat and lean spectrum using estimated breeding values derived from a total population size of over 70,000 animals. DNA from 3 breeds (Sire Line Large White, Duroc and a white Pietrain composite line (Titan)) was used to interrogate the Illumina Porcine SNP60 Genotyping Beadchip in order to identify significant associations in terms of single nucleotide polymorphisms (SNPs) and copy number variants (CNVs).

RESULTS

By sampling animals at each end of the fat/lean EBV (estimate breeding value) spectrum the whole population could be assessed using less than 300 animals, without losing statistical power. Indeed, several significant SNPs (at the 5% genome wide significance level) were discovered, 4 of these linked to genes with ontologies that had previously been correlated with fatness (NTS, FABP6, SST and NR3C2). Quantitative analysis of the data identified putative CNV regions containing genes whose ontology suggested fatness related functions (MCHR1, PPARα, SLC5A1 and SLC5A4).

CONCLUSIONS

Selective genotyping of EBVs at either end of the phenotypic spectrum proved to be a cost effective means of identifying SNPs and CNVs associated with fatness and with estimated major effects in a large population of animals.

Articular osteochondrosis: a comparison of naturally-occurring human and animal disease

BACKGROUND

Osteochondrosis (OC) is a common developmental orthopedic disease affecting both humans and animals. Despite increasing recognition of this disease among children and adolescents, its pathogenesis is incompletely understood because clinical signs are often not apparent until lesions have progressed to end-stage, and examination of cadaveric early lesions is not feasible. In contrast, both naturally-occurring and surgically-induced animal models of disease have been extensively studied, most notably in horses and swine, species in which OC is recognized to have profound health and economic implications. The potential for a translational model of human OC has not been recognized in the existing human literature.

OBJECTIVE

The purpose of this review is to highlight the similarities in signalment, predilection sites and clinical presentation of naturally-occurring OC in humans and animals and to propose a common pathogenesis for this condition across species.

STUDY DESIGN

Review.

METHODS

The published human and veterinary literature for the various manifestations of OC was reviewed. Peer-reviewed original scientific articles and species-specific review articles accessible in PubMed (US National Library of Medicine) were eligible for inclusion.

RESULTS

A broad range of similarities exists between OC affecting humans and animals, including predilection sites, clinical presentation, radiographic/MRI changes, and histological appearance of the end-stage lesion, suggesting a shared pathogenesis across species.

CONCLUSION

This proposed shared pathogenesis for OC between species implies that naturally-occurring and surgically-induced models of OC in animals may be useful in determining risk factors and for testing new diagnostic and therapeutic interventions that can be used in humans.

Gene expression profiling of articular cartilage reveals functional pathways and networks of candidate genes for osteochondrosis in pigs

Osteochondrosis (OC) is a joint disorder that frequently causes leg weakness in growing pigs, resulting in welfare problems and economic losses. We aimed to detect molecular pathways relevant to the emergence of the disease and to identify candidate genes for the liability to the disorder. Therefore, we compared microarray-based expression patterns of articular cartilage with (n=11) and without (n=11) histologically diagnosed OC lesions obtained from discordant sib-pairs. A total of 1,564 genes were found with different transcript abundance [differentially expressed (DE) genes] at q≤0.05. To further identify candidate genes, we integrated data from quantitative trait loci (QTL) and genome-wide association (GWA) studies with the expression analysis. We detected 317 DE genes within the QTL confidence intervals, of which 26 DE genes also overlapped GWA regions. Ingenuity Pathway Analysis suggests a pathogenic role of immune response, angiogenesis, and synthesis of extracellular matrix pathways for OC. These processes could facilitate the emergence of defects. But they may also promote the degradation of articular cartilage and the worsening of the disease. A functional network was derived that comprised genes with functional and positional clues of their role in bone and cartilage metabolisms and development, including extracellular matrix genes (e.g., LOX, OGN, and ASPN), angiogenesis genes (e.g., ANGPTL4 and PDGFA), and immune response genes (e.g., ICAM1, AZGP1, C1QB, C1QC, PDE4B, and CDA). The study identified molecular processes linked to OC and several genes with positional, genetic-statistical, and functional evidence for their role in the emergence of articular cartilage lesions and the liability to OC.

KRT8, FAF1 and PTH1R gene polymorphisms are associated with leg weakness traits in pigs

The liability to lesions of dysfunctions of bone and joints in pigs, summarized as leg weakness and mostly expressed as osteochondrosis, is an animal welfare and economic issue in pig production. The objective of this study was to identify polymorphisms in the functional and positional candidate genes keratin 8 (KRT8), Fas-associated factor 1 (FAF1) and parathyroid hormone type I receptor (PTH1R) and to evaluate their association with leg weakness traits. Therefore, osteochondrosis lesions were scored in animals of a Duroc × Pietrain F2 population (DuPi; n = 310) and commercial herds of the breed Large White (n = 298). In addition, bone mineralization traits were observed in DuPi population. SNPs were identified in genes KRT8 (g.8,039G > A), FAF1 (g.380,914T > C) and PTH1R (c.1,672C > T). KRT8 showed significant association with bone mineral density and content (P ≤ 0.05). FAF1 was association with OC lesions score of all joints inspected (P ≤ 0.05). PTH1R showed significant dominance effects on OC lesion scores of the distal femur articular cartilage (P = 0.01) and epiphysis of the distal ulna (P = 0.05) as well as sums of scores of all joints (OCsum, P = 0.04) and assignment to groups of either severely or gently affected animals (OCcat, P = 0.01). This study reveals clear genetic-statistical evidence for a link of KRT8, FAF1 and PTH1R with some of leg weakness related traits in pigs.

A genome wide association study for backfat thickness in Italian Large White pigs highlights new regions affecting fat deposition including neuronal genes

Background

Carcass fatness is an important trait in most pig breeding programs. Following market requests, breeding plans for fresh pork consumption are usually designed to reduce carcass fat content and increase lean meat deposition. However, the Italian pig industry is mainly devoted to the production of Protected Designation of Origin dry cured hams: pigs are slaughtered at around 160 kg of live weight and the breeding goal aims at maintaining fat coverage, measured as backfat thickness to avoid excessive desiccation of the hams. This objective has shaped the genetic pool of Italian heavy pig breeds for a few decades. In this study we applied a selective genotyping approach within a population of ~ 12,000 performance tested Italian Large White pigs. Within this population, we selectively genotyped 304 pigs with extreme and divergent backfat thickness estimated breeding value by the Illumina PorcineSNP60 BeadChip and performed a genome wide association study to identify loci associated to this trait.

Results

We identified 4 single nucleotide polymorphisms with P≤5.0E-07 and additional 119 ones with 5.0E-07<P≤5.0E-05. These markers were located throughout all chromosomes. The largest numbers were found on porcine chromosomes 6 and 9 (n=15), 4 (n=13), and 7 (n=12) while the most significant marker was located on chromosome 18. Twenty-two single nucleotide polymorphisms were in intronic regions of genes already recognized by the Pre-Ensembl Sscrofa10.2 assembly. Gene Ontology analysis indicated an enrichment of Gene Ontology terms associated with nervous system development and regulation in concordance with results of large genome wide association studies for human obesity.

Conclusions

Further investigations are needed to evaluate the effects of the identified single nucleotide polymorphisms associated with backfat thickness on other traits as a pre-requisite for practical applications in breeding programs. Reported results could improve our understanding of the biology of fat metabolism and deposition that could also be relevant for other mammalian species including humans, confirming the role of neuronal genes on obesity.

Genetics of behavioural adaptation of livestock to farming conditions

Behavioural adaptation of farm animals to environmental changes contributes to high levels of production under a wide range of farming conditions, from highly controlled indoor systems to harsh outdoor systems. The genetic variation in livestock behaviour is considerable. Animals and genotypes with a larger behavioural capacity for adaptation may cope more readily with varying farming conditions than those with a lower capacity for adaptation. This capacity should be exploited when the aim is to use a limited number of species extensively across the world. The genetics of behavioural traits is understood to some extent, but it is seldom accounted for in breeding programmes. This review summarizes the estimates of genetic parameters for behavioural traits in cattle, pigs, poultry and fish. On the basis of the major studies performed in the last two decades, we focus the review on traits of common interest in the four species. These concern the behavioural responses to both acute and chronic stressors in the physical environment (feed, temperature, etc.) and those in the social environment (other group members, progeny, humans). The genetic strategies used to improve the behavioural capacity for adaptation of animals differ between species. There is a greater emphasis on responses to acute environmental stress in fish and birds, and on responses to chronic social stress in mammals.

Association and expression study of MMP3, TGFβ1 and COL10A1 as candidate genes for leg weakness-related traits in pigs

The present study was aimed to determine the association between metalloproteinase 3 (MMP3), transforming growth factor beta 1 (TGFβ1) and collagen type X alpha I (COL10A1) gene polymorphisms with traits related to leg weakness in pigs. Three hundred Duroc × Pietrain cross breds (DuPi) and 299 pigs of a commercial population (CP) were used for the experiment. DuPi animals were examined for 10 different traits describing leg and feet structure, osteochondrosis (OC) scores and bone density status. Data of OC score at condylus medialis humeri, condylus medialis femoris and distal epiphysis ulna regions of CP were used for association analysis. Significant association (P < 0.05) was found for MMP3 SNP (g.158 C>T) with OC at head of femur and bone mineral density in the DuPi population. Association (P < 0.05) was found between SNP of TGFβ1 (g.180 G>A) with rear leg score and the principle component denoting both OC and feet and leg scores in the DuPi population. No association was found between COL10A1 (g.72 C>T) and leg weakness related traits. The associations of SNPs with OC traits could not be confirmed in the commercial population. Expression analysis of the three candidate genes was performed to compare between healthy and OC. TGFβ1 was found to be highly expressed (P < 0.05) in the OC compared to healthy cartilages, but no significant different expressions were observed for MMP3 and COL10A1 genes. The present finding suggested that TGFβ1 and MMP3 genes variants have an effect on some of the leg weakness related traits.

A whole-genome association study for pig reproductive traits

A whole-genome association study was performed for reproductive traits in commercial sows using the PorcineSNP60 BeadChip and Bayesian statistical methods. The traits included total number born (TNB), number born alive (NBA), number of stillborn (SB), number of mummified foetuses at birth (MUM) and gestation length (GL) in each of the first three parities. We report the associations of informative QTL and the genes within the QTL for each reproductive trait in different parities. These results provide evidence of gene effects having temporal impacts on reproductive traits in different parities. Many QTL identified in this study are new for pig reproductive traits. Around 48% of total genes located in the identified QTL regions were predicted to be involved in placental functions. The genomic regions containing genes important for foetal developmental (e.g. MEF2C) and uterine functions (e.g. PLSCR4) were associated with TNB and NBA in the first two parities. Similarly, QTL in other foetal developmental (e.g. HNRNPD and AHR) and placental (e.g. RELL1 and CD96) genes were associated with SB and MUM in different parities. The QTL with genes related to utero-placental blood flow (e.g. VEGFA) and hematopoiesis (e.g. MAFB) were associated with GL differences among sows in this population. Pathway analyses using genes within QTL identified some modest underlying biological pathways, which are interesting candidates (e.g. the nucleotide metabolism pathway for SB) for pig reproductive traits in different parities. Further validation studies on large populations are warranted to improve our understanding of the complex genetic architecture for pig reproductive traits.

A functional mutation at position -155 in porcine APOE promoter affects gene expression

BACKGROUND

Apolipoprotein E, a component of the plasma lipoproteins, plays an important role in the transport and metabolism of cholesterol and other lipids. Three single nucleotide polymorphisms (SNPs) -491A>T, -219T>G and +113G>C in the regulatory region of human apolipoprotein E gene (APOE) change the promoter activity and are associated with a wide variety of disorders including Alzheimer disease (AD). Functional SNPs in porcine APOE gene 5' regulatory region have not been explored.

RESULTS

We examined SNPs within this region (from -831 to +855), and the analysis revealed that the T>A SNP at position -155 among these three polymorphism sites (-440, -155, +501) was found to exert a marked influence on the transcription of the porcine APOE gene. Electrophoretic mobility shift assays showed that the binding affinity of oligonucletides containing the -155A to transcription factor(s) was stronger than that of the -155T. Transient transfection assays and quantitative real-time PCR results revealed that the -155T>A variant enhanced the activity of the APOE promoter and was associated with increased APOE mRNA levels in vivo.

CONCLUSIONS

These data suggest that the -155T>A mutation in the promoter region of the porcine APOE gene is an important functional variant. The results provided new insights into aspects of pig genetics and might also facilitate the application of pigs in biomedical studies addressing important human diseases.

Genome-wide association study identifies Loci for body composition and structural soundness traits in pigs

BACKGROUND

The recent completion of the swine genome sequencing project and development of a high density porcine SNP array has made genome-wide association (GWA) studies feasible in pigs.

METHODOLOGY/PRINCIPAL FINDINGS

Using Illumina's PorcineSNP60 BeadChip, we performed a pilot GWA study in 820 commercial female pigs phenotyped for backfat, loin muscle area, body conformation in addition to feet and leg (FL) structural soundness traits. A total of 51,385 SNPs were jointly fitted using Bayesian techniques as random effects in a mixture model that assumed a known large proportion (99.5%) of SNPs had zero effect. SNP annotations were implemented through the Sus scrofa Build 9 available from pig Ensembl. We discovered a number of candidate chromosomal regions, and some of them corresponded to QTL regions previously reported. We not only have identified some well-known candidate genes for the traits of interest, such as MC4R (for backfat) and IGF2 (for loin muscle area), but also obtained novel promising genes, including CHCHD3 (for backfat), BMP2 (for loin muscle area, body size and several FL structure traits), and some HOXA family genes (for overall leg action). The candidate regions responsible for body conformation and FL structure soundness did not overlap greatly which implied that these traits were controlled by different genes. Functional clustering analyses classified the genes into categories related to bone and cartilage development, muscle growth and development or the insulin pathway suggesting the traits are regulated by common pathways or gene networks that exert roles at different spatial and temporal stages.

CONCLUSIONS/SIGNIFICANCE

This study is one of the earliest GWA reports on important quantitative traits in pigs, and the findings will contribute to the further biological function analysis of the identified candidate genes and potential utilization of them in marker assisted selection.

The pig genome project has plenty to squeal about

Significant progress on pig genetics and genomics research has been witnessed in recent years due to the integration of advanced molecular biology techniques, bioinformatics and computational biology, and the collaborative efforts of researchers in the swine genomics community. Progress on expanding the linkage map has slowed down, but the efforts have created a higher-resolution physical map integrating the clone map and BAC end sequence. The number of QTL mapped is still growing and most of the updated QTL mapping results are available through PigQTLdb. Additionally, expression studies using high-throughput microarrays and other gene expression techniques have made significant advancements. The number of identified non-coding RNAs is rapidly increasing and their exact regulatory functions are being explored. A publishable draft (build 10) of the swine genome sequence was available for the pig genomics community by the end of December 2010. Build 9 of the porcine genome is currently available with Ensembl annotation; manual annotation is ongoing. These drafts provide useful tools for such endeavors as comparative genomics and SNP scans for fine QTL mapping. A recent community-wide effort to create a 60K porcine SNP chip has greatly facilitated whole-genome association analyses, haplotype block construction and linkage disequilibrium mapping, which can contribute to whole-genome selection. The future 'systems biology' that integrates and optimizes the information from all research levels can enhance the pig community's understanding of the full complexity of the porcine genome. These recent technological advances and where they may lead are reviewed.

Whole-genome association analyses for lifetime reproductive traits in the pig

Profits for commercial pork producers vary in part because of sow productivity or sow productive life (SPL) and replacement costs. During the last decade, culling rates of sows have increased to more than 50% in the United States. Both SPL and culling rates are influenced by genetic and nongenetic factors. A whole-genome association study was conducted for pig lifetime reproductive traits, including lifetime total number born (LTNB), lifetime number born alive (LNBA), removal parity, and the ratio between lifetime nonproductive days and herd life. The proportion of phenotypic variance explained by markers was 0.15 for LTNB and LNBA, 0.12 for removal parity, and 0.06 for the ratio between lifetime nonproductive days and herd life. Several informative QTL regions (e.g., 14 QTL regions for LTNB) and genes within the regions (e.g., SLC22A18 on SSC2 for LTNB) were associated with lifetime reproductive traits in this study. Genes associated with LTNB and LNBA were similar, reflecting the high genetic correlation (0.99 ± 0.003) between these traits. Functional annotation revealed that many genes at the associated regions are expressed in reproductive tissues. For instance, the SLC22A18 gene on SSC2 associated with LTNB has been shown to be expressed in the placenta of mice. Many of the QTL regions showing associations coincided with previously identified QTL for fat deposition. This reinforces the role of fat regulation for lifetime reproductive traits. Overall, this whole-genome association study provides a list of genomic locations and markers associated with pig lifetime reproductive traits that could be considered for SPL in future studies.

DNA sequence polymorphisms in a panel of eight candidate bovine imprinted genes and their association with performance traits in Irish Holstein-Friesian cattle

Background

Studies in mice and humans have shown that imprinted genes, whereby expression from one of the two parentally inherited alleles is attenuated or completely silenced, have a major effect on mammalian growth, metabolism and physiology. More recently, investigations in livestock species indicate that genes subject to this type of epigenetic regulation contribute to, or are associated with, several performance traits, most notably muscle mass and fat deposition. In the present study, a candidate gene approach was adopted to assess 17 validated single nucleotide polymorphisms (SNPs) and their association with a range of performance traits in 848 progeny-tested Irish Holstein-Friesian artificial insemination sires. These SNPs are located proximal to, or within, the bovine orthologs of eight genes (CALCR, GRB10, PEG3, PHLDA2, RASGRF1, TSPAN32, ZIM2 and ZNF215) that have been shown to be imprinted in cattle or in at least one other mammalian species (i.e. human/mouse/pig/sheep).

Results

Heterozygosities for all SNPs analysed ranged from 0.09 to 0.46 and significant deviations from Hardy-Weinberg proportions (P ≤ 0.01) were observed at four loci. Phenotypic associations (P ≤ 0.05) were observed between nine SNPs proximal to, or within, six of the eight analysed genes and a number of performance traits evaluated, including milk protein percentage, somatic cell count, culled cow and progeny carcass weight, angularity, body conditioning score, progeny carcass conformation, body depth, rump angle, rump width, animal stature, calving difficulty, gestation length and calf perinatal mortality. Notably, SNPs within the imprinted paternally expressed gene 3 (PEG3) gene cluster were associated (P ≤ 0.05) with calving, calf performance and fertility traits, while a single SNP in the zinc finger protein 215 gene (ZNF215) was associated with milk protein percentage (P ≤ 0.05), progeny carcass weight (P ≤ 0.05), culled cow carcass weight (P ≤ 0.01), angularity (P ≤ 0.01), body depth (P ≤ 0.01), rump width (P ≤ 0.01) and animal stature (P ≤ 0.01).

Conclusions

Of the eight candidate bovine imprinted genes assessed, DNA sequence polymorphisms in six of these genes (CALCR, GRB10, PEG3, RASGRF1, ZIM2 and ZNF215) displayed associations with several of the phenotypes included for analyses. The genotype-phenotype associations detected here are further supported by the biological function of these six genes, each of which plays important roles in mammalian growth, development and physiology. The associations between SNPs within the imprinted PEG3 gene cluster and traits related to calving, calf performance and gestation length suggest that this domain on chromosome 18 may play a role regulating pre-natal growth and development and fertility. SNPs within the bovine ZNF215 gene were associated with bovine growth and body conformation traits and studies in humans have revealed that the human ZNF215 ortholog belongs to the imprinted gene cluster associated with Beckwith-Wiedemann syndrome--a genetic disorder characterised by growth abnormalities. Similarly, the data presented here suggest that the ZNF215 gene may have an important role in regulating bovine growth. Collectively, our results support previous work showing that (candidate) imprinted genes/loci contribute to heritable variation in bovine performance traits and suggest that DNA sequence polymorphisms within these genes/loci represents an important reservoir of genomic markers for future genetic improvement of dairy and beef cattle populations.