Genetic analysis of teat number in pigs reveals some developmental pathways independent of vertebra number and several loci which only affect a specific side

Genome-wide detection of multiple variants associated with teat number in French Yorkshire pigs

BACKGROUND

Teat number is a vital reproductive trait in sows, crucial for providing immunity and nutrition to piglets during lactation. However, "missing heritability" in Single Nucleotide Polymorphism (SNP)-based Genome-Wide Association Studies (GWAS) has led to an increasing focus on structural variations in the genetic analysis of complex biological traits.

RESULTS

In this study, we generated a comprehensive CNV map in a population of French Yorkshire pigs (n = 644) and identified 429 CNVRs. Notably, 44% (189 CNVRs) of these were detected for the first time. Subsequently, we conducted GWAS for teat number in the French Yorkshire pig population using both 80K chip and its imputed data, as well as a GWAS analysis based on CNV regions (CNVR). Interestingly, 80K chip GWAS identified two SNPs located on Sus scrofa chromosome 5 (SSC5) that were simultaneously associated with Total Teat Number (TTN), Left Teat Number (LTN), and Right Teat Number (RTN). The leading SNP (WU_10.2_5_76130558) explained 3.33%, 2.69%, and 2.67% of the phenotypic variance for TTN, LTN, and RTN, respectively. Moreover, through imputed GWAS, we successfully identified 30 genetic variants associated with TTN located within the 73.22 -73.30 Mb region on SSC5. The two SNPs identified in the 80K chip GWAS were also located in this region. In addition, CNVR-based GWAS revealed three significant CNVRs associated with TTN. Finally, through gene annotation, we pinpointed two candidate genes, TRIM66 and PRICKLE1, which are related to diverse processes such as breast cancer and abnormal vertebral development.

CONCLUSIONS

Our research provides an in-depth analysis of the complex genetic structure underlying teat number, contributing to the genetic enhancement of sows with improved reproductive performance and, ultimately, bolstering the economic benefits of swine production enterprises.

Genome-Wide Association Study Meta-Analysis Elucidates Genetic Structure and Identifies Candidate Genes of Teat Number Traits in Pigs

The teat number is a pivotal reproductive trait that significantly influences the survival rate of piglets. A meta-analysis is a robust instrument, enhancing the universality of research findings and improving statistical power by increasing the sample size. This study aimed to identify universal candidate genes associated with teat number traits using a genome-wide association study (GWAS) meta-analysis with three breeds. We identified 21 chromosome threshold significant single-nucleotide polymorphisms (SNPs) associated with five teat number traits in single-breed and cross-breed meta-GWAS analyses. Using a co-localization analysis of expression quantitative trait loci and GWAS loci, we detected four unique genes that were co-localized with cross-breed GWAS loci associated with teat number traits. Through a meta-analysis and integrative analysis, we identified more reliable candidate genes associated with multiple-breed teat number traits. Our research provides new information for exploring the genetic mechanism affecting pig teat number for breeding selection and improvement.

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

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

CYP24A1 is associated with fetal mummification in pigs

Mummified piglets are among the leading causes of fertility loss and severely hamper reproductive performance in pigs. However, the contributions of genomic variation to the emergence of mummified piglets (MUM) have rarely been studied. This study aims to (1) elucidate the genetic architecture of MUM in sows of parity 1 - 3 using a single-step genome-wide association study (ssGWAS). The ssGWAS involved genotyping-by-sequencing of Large White and Landrace pig breeds. (2) Explore the biological role of the candidate genes at the cellular level. A total of 185 and 48 genome-wide significant SNPs are associated with MUM in Large White and Landrace pigs, explaining 0.01-36.52% genetic variance for different significant loci, respectively. All the significant SNPs are parity-specific, and the numerous, consecutive significant loci likely generated the nine significant peaks in different parities. Multiple candidate genes (including CYP24A1, FBXO30, and ARHGEF28) are associated with fetal congenital and maternal diseases. Collectively, CYP24A1 regulation contributes to steady-state levels of embryo development genes. CYP24A1 is involved in reproduction and, immune and gestational disorders. Thus, it is associated with known newborn death traits and MUM in Large White sows. Altogether, these results improve the current understanding of the genetic architecture of MUM and expand the knowledge on genetic variations for selecting against mummified piglets in pig breeding.

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

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

The impact of functional teat number on reproductive throughput in swine

The objective was to evaluate the impact of functional teat number on reproductive throughput in swine. Data included 735 multiparous Landrace × Large White F1 females. Sow underlined traits consisted of total teat number (TT), functional teat number (FT), nonfunctional teat number (NFT), and number of functional mammary glands (FMG). Weaning traits were calculated for both the biological and the nurse dam. For the biological dam, litter size at weaning (LSW) included a sow's biological piglets regardless of cross-fostering. For nurse dam, number weaned (NW) included the piglets a sow weaned. For the biological dam, piglet survival (PS) was calculated as litter size at weaning / (total number born × 100). Linear regression estimates were calculated in RStudio v. 1.1.456 and variance components were estimated using GIBBS3F90. Average total number born, number born alive, TT, FT, NFT, and FMG were 14.22, 13.12, 14.43, 13.96, 0.42, and 10.7, respectively. An increase in one FT enhanced (P < 0.05) LSW by 0.32 piglets and NW by 0.33 piglets. Similarly, an increase in one FT improved (P < 0.05) PS by 1.63% and reduced (P < 0.05) preweaning mortality by 2.73%. However, an increase in one FT reduced (P < 0.05) average piglet weaning weight (WW) for biological and nurse dams by 35 and 94 g, respectively. Yet an increase in one FT enhanced (P < 0.05) litter weaning weight (LWW) for biological and nurse dams by 1.3 and 1.5 kg, respectively. Heritability estimates for TT, FT, NFT, FMG, WW, LWW, LSW, and PS were 0.25, 0.22, 0.53, 0.18, 0.21, 0.22, 0.16, and 0.18, respectively. Genetic correlation estimates between FT with TT, NFT, and FMG were 0.79, 0.09, and 0.28, respectively. Estimated genetic correlations between TT with WW, LWW, LSW, and PS were 0.37, 0.38, 0.11, and -0.19, respectively. Genetic correlation estimates between FT with WW, LWW, LSW, and PS were 0.44, 0.49, 0.39, and 0.35, respectively. Results suggest increasing functional teat number would enhance both piglet survival and reproductive throughput.

Genetic Parameter Estimates for Teat and Mammary Traits in Commercial Sows

The objective was to evaluate the genetics of sow teat and mammary traits at farrowing and at weaning. Data were recorded on 3099 Landrace × Large White F1 sows. Underline traits included the total teat number (TT), the functional teat number (FT), the non-functional teat number (NFT), the damaged teat number (DT), and the number of functional mammary glands (FMG). Variance components were estimated using AIREMLF90. Means for TT, FT, and NFT at farrowing were 14.93, 13.90, and 1.03, respectively. Heritability estimates for TT, FT, and NFT ranged from 0.18 to 0.37, 0.16 to 0.28, and 0.14 to 0.18, respectively. Estimates of heritability for DT and FMG at weaning were 0.03 and 0.06, respectively. Estimated genetic correlations between FT with TT and NFT were 0.68 to 0.78 and -0.19 to -0.57, respectively. Genetic correlation estimates between TT, FT, and NFT with the number weaned were 0.25, 0.50, and -0.38, respectively. An increase of one TT and FT enhanced (p < 0.05) the number weaned by 0.14 to 0.16 and 0.18 to 0.27 piglets, respectively. The results suggest that genetically increasing the number of functional teats on a sow at farrowing would improve the number of piglets at weaning.

Identification of ROH Islands Conserved through Generations in Pigs Belonging to the Nero Lucano Breed

The recovery of Nero Lucano (NL) pigs in the Basilicata region (Southern Italy) started in 2001 with the collaboration of several public authorities in order to preserve native breeds that can play a significant economic role both due to their remarkable ability to adapt to difficult environments and the value of typical products from their area of origin. In this study, by using the Illumina Porcine SNP60 BeadChip, we compared the genetic structures of NL pigs reared in a single farm in two different periods separated by a time interval corresponding to at least three generations. The results showed an increase in the percentage of polymorphic loci, a decrease in the inbreeding coefficient calculated according to ROH genome coverage (F_ROH), a reduction in the number of ROH longer than 16 Mb and an increase in ROH with a length between 2 and 4 Mb, highlighting a picture of improved genetic variability. In addition, ROH island analysis in the two groups allowed us to identify five conserved regions, located on chromosomes 1, 4, 8, 14 and 15, containing genes involved in biological processes affecting immune response, reproduction and production traits. Only the conserved ROH island on chromosome 14 contains markers which, according to the literature, are associated with QTLs affecting thoracic vertebra number, teat number, gestation length, age at puberty and mean platelet volume.

Genome-Wide Association Study Identifies the Crucial Candidate Genes for Teat Number in Crossbred Commercial Pigs

The number of teats is a crucial reproductive trait with significant economic implications on maternal capacity and litter size. Consequently, improving this trait is essential to facilitate genetic selection for increased litter size. In this study, we performed a genome-wide association study (GWAS) of the number of teats in a three-way crossbred commercial Duroc × (Landrace × Yorkshire) (DLY) pig population comprising 1518 animals genotyped with the 50K BeadChip. Our analysis identified crucial quantitative trait loci (QTL) for the number of teats, containing the ABCD4 and VRTN genes on porcine chromosome 7. Our results establish SNP variants of ABCD4 and VRTN as new molecular markers for improving the number of teats in DLY pigs. Furthermore, the most significant noteworthy single nucleotide polymorphism (SNP) (7_97568284) was identified within the ABCD4 gene, exhibiting a significant association with the total teat number traits. This SNP accounted for a substantial proportion of the genetic variance, explaining 6.64% of the observed variation. These findings reveal a novel gene on SSC7 for the number of teats and provide a deeper understanding of the genetic mechanisms underlying reproductive traits.

Estimating genetics of body dimensions and activity levels in pigs using automated pose estimation

Pig breeding is changing rapidly due to technological progress and socio-ecological factors. New precision livestock farming technologies such as computer vision systems are crucial for automated phenotyping on a large scale for novel traits, as pigs’ robustness and behavior are gaining importance in breeding goals. However, individual identification, data processing and the availability of adequate (open source) software currently pose the main hurdles. The overall goal of this study was to expand pig weighing with automated measurements of body dimensions and activity levels using an automated video-analytic system: DeepLabCut. Furthermore, these data were coupled with pedigree information to estimate genetic parameters for breeding programs. We analyzed 7428 recordings over the fattening period of 1556 finishing pigs (Piétrain sire x crossbred dam) with two-week intervals between recordings on the same pig. We were able to accurately estimate relevant body parts with an average tracking error of 3.3 cm. Body metrics extracted from video images were highly heritable (61–74%) and significantly genetically correlated with average daily gain (r_g = 0.81–0.92). Activity traits were low to moderately heritable (22–35%) and showed low genetic correlations with production traits and physical abnormalities. We demonstrated a simple and cost-efficient method to extract body dimension parameters and activity traits. These traits were estimated to be heritable, and hence, can be selected on. These findings are valuable for (pig) breeding organizations, as they offer a method to automatically phenotype new production and behavioral traits on an individual level.

Welfare of pigs on farm

This scientific opinion focuses on the welfare of pigs on farm, and is based on literature and expert opinion. All pig categories were assessed: gilts and dry sows, farrowing and lactating sows, suckling piglets, weaners, rearing pigs and boars. The most relevant husbandry systems used in Europe are described. For each system, highly relevant welfare consequences were identified, as well as related animal-based measures (ABMs), and hazards leading to the welfare consequences. Moreover, measures to prevent or correct the hazards and/or mitigate the welfare consequences are recommended. Recommendations are also provided on quantitative or qualitative criteria to answer specific questions on the welfare of pigs related to tail biting and related to the European Citizen's Initiative 'End the Cage Age'. For example, the AHAW Panel recommends how to mitigate group stress when dry sows and gilts are grouped immediately after weaning or in early pregnancy. Results of a comparative qualitative assessment suggested that long-stemmed or long-cut straw, hay or haylage is the most suitable material for nest-building. A period of time will be needed for staff and animals to adapt to housing lactating sows and their piglets in farrowing pens (as opposed to crates) before achieving stable welfare outcomes. The panel recommends a minimum available space to the lactating sow to ensure piglet welfare (measured by live-born piglet mortality). Among the main risk factors for tail biting are space allowance, types of flooring, air quality, health status and diet composition, while weaning age was not associated directly with tail biting in later life. The relationship between the availability of space and growth rate, lying behaviour and tail biting in rearing pigs is quantified and presented. Finally, the panel suggests a set of ABMs to use at slaughter for monitoring on-farm welfare of cull sows and rearing pigs.

Genetic aspects of piglet survival and related traits: a review

In livestock, mortality in general, and mortality of the young, is societal worries and is economically relevant for farm efficiency. Genetic change is cumulative; if it exists for survival of the young and genetic merit can be estimated with sufficient accuracy, it can help alleviate the pressure of mortality. Lack of survival is a moving target; livestock production is in continuous change and labor shortage is a given. There is now ample evidence of clear genetic variance and of models able to provide genomic predictions with enough accuracy for selection response. Underlying traits such as birth weight, uniformity in birth weight, gestation length, number of teats, and farrowing duration all show genetic variation and support selection for survival or, alternatively, be selected for on their own merit.

A Whole Genome Sequencing-Based Genome-Wide Association Study Reveals the Potential Associations of Teat Number in Qingping Pigs

Teat number plays an important role in the reproductive performance of sows and the growth of piglets. However, the quantitative trait loci (QTLs) and candidate genes for the teat number-related traits in Qingping pigs remain unknown. In this study, we performed GWAS based on whole-genome single-nucleotide polymorphisms (SNPs) and insertions/deletions (Indels) for the total number of teats and five other related traits in 100 Qingping pigs. SNPs and Indels of all 100 pigs were genotyped using 10× whole genome resequencing. GWAS using General Linear Models (GLM) detected a total of 28 SNPs and 45 Indels as peak markers for these six traits. We also performed GWAS for the absolute difference between left and right teat number (ADIFF) using Fixed and random model Circulating Probability Unification (FarmCPU). The most strongly associated SNP and Indel with a distance of 562,788 bp were significantly associated with ADIFF in both GLM and FarmCPU models. In the 1-Mb regions of the most strongly associated SNP and Indel, there were five annotated genes, including TRIML1, TRIML2, ZFP42, FAT1 and MTNR1A. We also highlighted TBX3 as an interesting candidate gene for SSC14. Enrichment analysis of candidate genes suggested the Wnt signaling pathway may contribute to teat number-related traits. This study expanded significant marker-trait associations for teat number and provided useful molecular markers and candidate genes for teat number improvement in the breeding of sows.

Genome-Wide Association Study for Body Length, Body Height, and Total Teat Number in Large White Pigs

Body length, body height, and total teat number are economically important traits in pig breeding, as these traits are usually associated with the growth, reproductivity, and longevity potential of piglets. Here, we report a genetic analysis of these traits using a population comprising 2,068 Large White pigs. A genotyping-by-sequencing (GBS) approach was used to provide high-density genome-wide SNP discovery and genotyping. Univariate and bivariate animal models were used to estimate heritability and genetic correlations. The results showed that heritability estimates for body length, body height, and total teat number were 0.25 ± 0.04, 0.11 ± 0.03, and 0.22 ± 0.04, respectively. The genetic correlation between body length and body height exhibited a strongly positive correlation (0.63 ± 0.15), while a positive but low genetic correlation was observed between total teat number and body length. Furthermore, we used two different genome-wide association study (GWAS) approaches: single-locus GWAS and weighted single-step GWAS (WssGWAS), to identify candidate genes for these traits. Single-locus GWAS detected 76, 13, and 29 significant single-nucleotide polymorphisms (SNPs) associated with body length, body height, and total teat number. Notably, the most significant SNP (S17_15781294), which is located 20 kb downstream of the BMP2 gene, explained 9.09% of the genetic variance for body length traits, and it also explained 9.57% of the genetic variance for body height traits. In addition, another significant SNP (S7_97595973), which is located in the ABCD4 gene, explained 8.92% of the genetic variance for total teat number traits. GWAS results for these traits identified some candidate genomic regions, such as SSC6: 14.96–15.02 Mb, SSC7: 97.18–98.18 Mb, SSC14: 128.29–131.15 Mb, SSC17: 15.39–17.27 Mb, and SSC17: 22.04–24.15 Mb, providing a starting point for further inheritance research. Most quantitative trait loci were detected by single-locus GWAS and WssGWAS. These findings reveal the complexity of the genetic mechanism of the three traits and provide guidance for subsequent genetic improvement through genome selection.

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

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

Expression analysis and single-nucleotide polymorphisms of SYNDIG1L and UNC13C genes associated with thoracic vertebral numbers in sheep (Ovis aries)

The objective of the current study was to analyze expression levels of synapse differentiation inducing 1-like (SYNDIG1L) and unc-13 homolog C (UNC13C) genes in different tissues, while single-nucleotide polymorphisms (SNPs) of two genes were associated with multiple thoracic vertebrae traits in both Small-tailed Han sheep (STH) and Sunite sheep (SNT). The expression levels of SYNDIG1L and UNC13C were analyzed in the brain, cerebellum, heart, liver, spleen, lung, kidney, adrenal gland, uterine horn, longissimus muscle, and abdominal adipose tissues of two sheep breeds with different thoracic vertebral number (TVN) sheep (T13 groups and T14 groups) by real-time quantitative polymerase chain reaction (RT-qPCR). Meanwhile, the polymorphisms of UNC13C gene g.52919279C>T and SYNDIG1L gene g.82573325C>A in T14 and T13 were genotyped by the Sequenom MassARRAY^® SNP assay, and association analysis was performed with the TVN. The results demonstrated that UNC13C gene was extensively expressed in 11 tissues. The expression of UNC13C gene in longissimus muscle of T14 groups of STH was significantly higher than that of T13 groups (P<0.05). SYNDIG1L gene was overexpressed in brain and cerebellum tissues, and the expression level of UNC13C gene in the brain and cerebellum of T13 groups in SNT was significantly higher than that of T14 groups (P<0.01). Association analysis showed that SNPs found in the UNC13C gene had no significant effects on TVN for both two genes. The polymorphism of SYNDIG1L g.82573325C>A was significantly correlated with the TVN in both STH (P<0.05) and SNT (P<0.01). Taken together, the SYNDIG1L gene was related to thoracic vertebral development, and this variation may be potentially used as a molecular marker to select the multiple thoracic vertebrae in sheep.

Revealing New Candidate Genes for Teat Number Relevant Traits in Duroc Pigs Using Genome-Wide Association Studies

Simple Summary

Number of teats is very important for lactating sows. We conducted genome-wide association studies (GWAS) and estimated the genetic parameters for traits related to teat number. Results showed that there were nine and 22 SNPs exceeding genome-wide significance and suggestive significance levels, respectively. Eighteen genes annotated near them were concentrated on chromosomes 7 and 10. Among them, three new candidate genes were located on the genomic regions around the significant SNPs. Our findings provide new insight into investigating the complex genetic mechanism of traits related to teat number in pigs.

Abstract

The number of teats is related to the nursing ability of sows. In the present study, we conducted genome-wide association studies (GWAS) for traits related to teat number in Duroc pig population. Two mixed models, one for counted and another for binary phenotypic traits, were employed to analyze seven traits: the right (RTN), left (LTN), and total (TTN) teat numbers; maximum teat number on a side (MAX); left minus right side teat number (LR); the absolute value of LR (ALR); and the presence of symmetry between left and right teat numbers (SLR). We identified 11, 1, 4, 13, and 9 significant SNPs associated with traits RTN, LTN, MAX, TTN, and SLR, respectively. One significant SNP (MARC0038565) was found to be simultaneous associated with RTN, LTN, MAX and TTN. Two annotated genes (VRTN and SYNDIG1L) were located in genomic region around this SNP. Three significant SNPs were shown to be associated with TTN, RTN and MAX traits. Seven significant SNPs were simultaneously detected in two traits of TTN and RTN. Other two SNPs were only identified in TTN. These 13 SNPs were clustered in the genomic region between 96.10—98.09 Mb on chromosome 7. Moreover, nine significant SNPs were shown to be significantly associated with SLR. In total, four and 22 SNPs surpassed genome-wide significance and suggestive significance levels, respectively. Among candidate genes annotated, eight genes have documented association with the teat number relevant traits. Out of them, DPF3 genes on Sus scrofa chromosome (SSC) 7 and the NRP1 gene on SSC 10 were new candidate genes identified in this study. Our findings demonstrate the genetic mechanism of teat number relevant traits and provide a reference to further improve reproductive performances in practical pig breeding programs.

A genome-wide association study for the number of teats in European rabbits (Oryctolagus cuniculus) identifies several candidate genes affecting this trait

In the European rabbit (Oryctolagus cuniculus), a polytocous livestock species, the number of teats indirectly impacts the doe reproduction efficiency and, in turn, the sustainable production of rabbit meat. In this study, we carried out a genome-wide association study (GWAS) for the total number of teats in 247 Italian White does included in the Italian White rabbit breed selection program, by applying a selective genotyping approach. Does had either 8 (n = 121) or 10 teats (n = 126). All rabbits were genotyped with the Affymetrix Axiom OrcunSNP Array. Genomic data from the two extreme groups of rabbits were also analysed with the single-marker fixation index statistic and combined with the GWAS results. The GWAS identified 50 significant SNPs and the fixation index analysis identified a total of 20 SNPs that trespassed the 99.98th percentile threshold, 19 of which confirmed the GWAS results. The most significant SNP (P = 4.31 × 10^-11 ) was located on OCU1, close to the NUDT2 gene, a breast carcinoma cells proliferation promoter. Another significant SNP identified as candidate gene NR6A1, which is well known to play an important role in affecting the correlated number of vertebrae in pigs. Other significant markers were close to candidate genes involved in determining body length in mice. Markers associated with increased number of teats could be included in selection programmes to speed up the improvement for this trait in rabbit lines that need to increase maternal performances.

Genome-wide association analyses identify known and novel loci for teat number in Duroc pigs using single-locus and multi-locus models

BACKGROUND

More teats are necessary for sows to nurse larger litters to provide immunity and nutrient for piglets prior to weaning. Previous studies have reported the strong effect of an insertion mutation in the Vertebrae Development Associated (VRTN) gene on Sus scrofa chromosome 7 (SSC7) that increased the number of thoracic vertebrae and teat number in pigs. We used genome-wide association studies (GWAS) to map genetic markers and genes associated with teat number in two Duroc pig populations with different genetic backgrounds. A single marker method and several multi-locus methods were utilized. A meta-analysis that combined the effects and P-values of 34,681 single nucleotide polymorphisms (SNPs) that were common in the results of single marker GWAS of American and Canadian Duroc pigs was conducted. We also performed association tests between the VRTN insertion and teat number in the same populations.

RESULTS

A total of 97 SNPs were found to be associated with teat number. Among these, six, eight and seven SNPs were consistently detected with two, three and four multi-locus methods, respectively. Seven SNPs were concordantly identified between single marker and multi-locus methods. Moreover, 26 SNPs were newly found by multi-locus methods to be associated with teat number. Notably, we detected one consistent quantitative trait locus (QTL) on SSC7 for teat number using single-locus and meta-analysis of GWAS and the top SNP (rs692640845) explained 8.68% phenotypic variance of teat number in the Canadian Duroc pigs. The associations between the VRTN insertion and teat number in two Duroc pig populations were substantially weaker. Further analysis revealed that the effect of VRTN on teat number may be mediated by its LD with the true causal mutation.

CONCLUSIONS

Our study suggested that VRTN insertion may not be a strong or the only candidate causal mutation for the QTL on SSC7 for teat number in the analyzed Duroc pig populations. The combination of single-locus and multi-locus GWAS detected additional SNPs that were absent using only one model. The identified SNPs will be useful for the genetic improvement of teat number in pigs by assigning higher weights to associated SNPs in genomic selection.

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

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

Estimation of economic weights for number of teats and sperm quality traits in pigs

Accurate predictive modelling facilitates efficient and effective trait selection in animal breeding and can decrease costs while maximizing profits when raising economically important animals. The objective of this study was to extend a previously developed bioeconomic model and computer program to calculate the marginal economic values (MEVs) and economic weights (EWs) for direct and maternal pig traits affected by new reproductive traits, namely the number of sow functional teats (NFTs) and boar sperm quality traits (SQTs) that included sperm volume, sperm concentration, motility percentage and percentage of abnormal spermatozoa. The MEV of NFTs represented the cost differences between naturally and artificially reared piglets until weaning and the cost differences between naturally and artificially reared finished animals. The MEVs of SQTs expressed the saved costs for artificial insemination, assuming a decreased price per insemination dose when improving the SQTs. The absolute and relative EWs for the newly defined complex of traits in the breeding objectives for pig breeds involved in the Czech national three-way crossing system (Czech Large White [CLW], Czech Landrace [CL] and Pietrain [PN]) were calculated using gene flow methods. The NFT trait was included only for dam breeds, and the relative EW averaged 3.6% of the total economic importance based on the genetic standard deviations of all 19 simultaneously evaluated traits in CLW and CL breeds. The relative EWs of the four SQTs comprised 2.0% of the total economic importance of the 19 traits in the CLW and CL dam breeds and 8% of the total economic importance of the 18 traits in the PN sire breed. Therefore, inclusion of the NFTs for dam breeds and SQTs for sire breeds in the breeding goal is recommended to aid in obtaining ideal outcomes with optimal economic values.

A Survey of Copy Number Variation in the Porcine Genome Detected From Whole-Genome Sequence

Copy number variations (CNVs) are gains and losses of large regions of genomic sequence between individuals of a species. Although CNVs have been associated with various phenotypic traits in humans and other species, the extent to which CNVs impact phenotypic variation remains unclear. In swine, as well as many other species, relatively little is understood about the frequency of CNV in the genome, sizes, locations, and other chromosomal properties. In this work, we identified and characterized CNV by utilizing whole-genome sequence from 240 members of an intensely phenotyped experimental swine herd at the U.S. Meat Animal Research Center (USMARC). These animals included all 24 of the purebred founding boars (12 Duroc and 12 Landrace), 48 of the founding Yorkshire-Landrace composite sows, 109 composite animals from generations 4 through 9, 29 composite animals from generation 15, and 30 purebred industry boars (15 Landrace and 15 Yorkshire) used as sires in generations 10 through 15. Using a combination of split reads, paired-end mapping, and read depth approaches, we identified a total of 3,538 copy number variable regions (CNVRs), including 1,820 novel CNVRs not reported in previous studies. The CNVRs covered 0.94% of the porcine genome and overlapped 1,401 genes. Gene ontology analysis identified that CNV-overlapped genes were enriched for functions related to organism development. Additionally, CNVRs overlapped with many known quantitative trait loci (QTL). In particular, analysis of QTL previously identified in the USMARC herd showed that CNVRs were most overlapped with reproductive traits, such as age of puberty and ovulation rate, and CNVRs were significantly enriched for reproductive QTL.

A QTL for Number of Teats Shows Breed Specific Effects on Number of Vertebrae in Pigs: Bridging the Gap Between Molecular and Quantitative Genetics

Modern breeding schemes for livestock species accumulate a large amount of genotype and phenotype data which can be used for genome-wide association studies (GWAS). Many chromosomal regions harboring effects on quantitative traits have been reported from these studies, but the underlying causative mutations remain mostly undetected. In this study, we combine large genotype and phenotype data available from a commercial pig breeding scheme for three different breeds (Duroc, Landrace, and Large White) to pinpoint functional variation for a region on porcine chromosome 7 affecting number of teats (NTE). Our results show that refining trait definition by counting number of vertebrae (NVE) and ribs (RIB) helps to reduce noise from other genetic variation and increases heritability from 0.28 up to 0.62 NVE and 0.78 RIB in Duroc. However, in Landrace, the effect of the same QTL on NTE mainly affects NVE and not RIB, which is reflected in reduced heritability for RIB (0.24) compared to NVE (0.59). Further, differences in allele frequencies and accuracy of rib counting influence genetic parameters. Correction for the top SNP does not detect any other QTL effect on NTE, NVE, or RIB in Landrace or Duroc. At the molecular level, haplotypes derived from 660K SNP data detects a core haplotype of seven SNPs in Duroc. Sequence analysis of 16 Duroc animals shows that two functional mutations of the Vertnin (VRTN) gene known to increase number of thoracic vertebrae (ribs) reside on this haplotype. In Landrace, the linkage disequilibrium (LD) extends over a region of more than 3 Mb also containing both VRTN mutations. Here, other modifying loci are expected to cause the breed-specific effect. Additional variants found on the wildtype haplotype surrounding the VRTN region in all sequenced Landrace animals point toward breed specific differences which are expected to be present also across the whole genome. This Landrace specific haplotype contains two missense mutations in the ABCD4 gene, one of which is expected to have a negative effect on the protein function. Together, the integration of largescale genotype, phenotype and sequence data shows exemplarily how population parameters are influenced by underlying variation at the molecular level.

Genetic selection against intrauterine growth retardation in piglets: a problem at the piglet level with a solution at the sow level

Background

In polytocous livestock species, litter size and offspring weight act antagonistically; in modern pig breeds, selection for increased litter size has resulted in lower mean birth weights, an increased number of small piglets and an increased number of those affected by varying degrees of intrauterine growth retardation (IUGR). IUGR poses life-long challenges, both mental, with morphological brain changes and altered cognition, and physical, such as immaturity of organs, reduced colostrum intake and weight gain. In pigs, head morphology of newborn piglets is a good phenotypic marker for identifying such compromised piglets. Growth retardation could be considered as a property of the dam, in part due to either uterine capacity or insufficiency. A novel approach to this issue is to consider the proportion of IUGR-affected piglets in a litter as an indirect measure of uterine capacity. However, uterine capacity or sufficiency cannot be equated solely to litter size and thus is a trait difficult to measure on farm.

Results

A total of 21,159 Landrace × Large White or Landrace × White Duroc piglets (born over 52 weeks) with recorded head morphology and birth weights were followed from birth until death or weaning. At the piglet level, the estimated heritability for IUGR (as defined by head morphology) was low at 0.01 ± 0.01. Piglet direct genetic effects of birth weight (h² = 0.07 ± 0.02) were strongly negatively correlated with head morphology (− 0.93), in that IUGR-affected piglets tended to have lower birth weights. At the sow level, analysis of the proportion of IUGR-affected piglets in a litter gave a heritability of 0.20 ± 0.06, with high and negative genetic correlations of the proportion of IUGR-affected piglets with average offspring birth weight (− 0.90) and with the proportion of piglets surviving until 24 h (− 0.80).

Conclusions

This suggests that the proportion of IUGR-affected piglets in a litter is a suitable indirect measure of uterine capacity for inclusion in breeding programmes that aim at reducing IUGR in piglets and improving piglet survival.