A study of vertebra number in pigs confirms the association of vertnin and reveals additional QTL

Genome-Wide Association Studies of Growth Trait Heterosis in Crossbred Meat Rabbits

The application of heterosis can not only effectively improve the disease resistance and meat quality of livestock, but also significantly enhance the reproduction and growth of livestock and poultry. We conducted genome-wide association studies using data from F2 crossbred meat rabbits to screen out candidate genes with significant dominant effects associated with economic trait variation. High-throughput sequencing technology was used to obtain SNPs covering the whole genome to evaluate the homozygosity of the population genome, and analyze the number, length, frequency, and distribution of ROHs in the population. Candidate genes related to economic traits of meat rabbits were searched based on high-frequency ROH regions. After quality control filtering of genotype data, 380 F2 crossbred rabbits were identified with 78,579 SNPs and 42,018 ROHs on the autosomes. The fitting of the Logistic growth curve model showed that 49-day-old rabbits were a growth inflection point. Then, through genome-wide association studies, 10 SNP loci and seven growth trait candidate genes were found to be significantly related to body weight in meat rabbits at 84 days of age. In addition, we revealed the functional roles and locations of 20 candidate genes in the high-frequency ROH region associated with economic traits in meat rabbits. This study identified potential genes associated with growth and development in the high-frequency ROH region of meat rabbits. In this study, the identified candidate genes can be used as molecular markers for assisted selection in meat rabbits. At the same time, the inbreeding situation based on ROH assessment can provide reference for breeding and breeding preservation of meat rabbits.

Advancements in Genetic Marker Exploration for Livestock Vertebral Traits with a Focus on China

In livestock breeding, the number of vertebrae has gained significant attention due to its impact on carcass quality and quantity. Variations in vertebral traits have been observed across different animal species and breeds, with a strong correlation to growth and meat production. Furthermore, vertebral traits are classified as quantitative characteristics. Molecular marker techniques, such as marker-assisted selection (MAS), have emerged as efficient tools to identify genetic markers associated with vertebral traits. In the current review, we highlight some key potential genes and their polymorphisms that play pivotal roles in controlling vertebral traits (development, length, and number) in various livestock species, including pigs, donkeys, and sheep. Specific genetic variants within these genes have been linked to vertebral development, number, and length, offering valuable insights into the genetic mechanisms governing vertebral traits. This knowledge has significant implications for selective breeding strategies to enhance structural characteristics and meat quantity and quality in livestock, ultimately improving the efficiency and quality of the animal husbandry industry.

A meta-analysis of genetic and phenotypic diversity of European local pig breeds reveals genomic regions associated with breed differentiation for production traits

BACKGROUND

Intense selection of modern pig breeds has resulted in genetic improvement of production traits while the performance of local pig breeds has remained lower. As local pig breeds have been bred in extensive systems, they have adapted to specific environmental conditions, resulting in a rich genotypic and phenotypic diversity. This study is based on European local pig breeds that have been genetically characterized using DNA-pool sequencing data and phenotypically characterized using breed level phenotypes related to stature, fatness, growth, and reproductive performance traits. These data were analyzed using a dedicated approach to detect signatures of selection linked to phenotypic traits in order to uncover potential candidate genes that may underlie adaptation to specific environments.

RESULTS

Analysis of the genetic data of European pig breeds revealed four main axes of genetic variation represented by the Iberian and three modern breeds (i.e. Large White, Landrace, and Duroc). In addition, breeds clustered according to their geographical origin, for example French Gascon and Basque breeds, Italian Apulo Calabrese and Casertana breeds, Spanish Iberian, and Portuguese Alentejano breeds. Principal component analysis of the phenotypic data distinguished the larger and leaner breeds with better growth potential and reproductive performance from the smaller and fatter breeds with low growth and reproductive efficiency. Linking the signatures of selection with phenotype identified 16 significant genomic regions associated with stature, 24 with fatness, 2 with growth, and 192 with reproduction. Among them, several regions contained candidate genes with possible biological effects on stature, fatness, growth, and reproductive performance traits. For example, strong associations were found for stature in two regions containing, respectively, the ANXA4 and ANTXR1 genes, for fatness in a region containing the DNMT3A and POMC genes and for reproductive performance in a region containing the HSD17B7 gene.

CONCLUSIONS

In this study on European local pig breeds, we used a dedicated approach for detecting signatures of selection that were supported by phenotypic data at the breed level to identify potential candidate genes that may have adapted to different living environments and production systems.

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.

Integrated Analysis of Transcriptome Expression Profiles Reveals miRNA-326-NKX3.2-Regulated Porcine Chondrocyte Differentiation

The porcine body length trait is an essential factor affecting meat production and reproductive performance. It is evident that the development/lengthening of individual vertebrae is one of the main reasons for increases in body length; however, the underlying molecular mechanism remains unclear. In this study, RNA-seq analysis was used to profile the transcriptome (lncRNA, mRNA, and miRNA) of the thoracic intervertebral cartilage (TIC) at two time points (1 and 4 months) during vertebral column development in Yorkshire (Y) and Wuzhishan pigs (W). There were four groups: 1- (Y1) and 4-month-old (Y4) Yorkshire pigs and 1- (W1) and 4-month-old (W4) Wuzhishan pigs. In total, 161, 275, 86, and 126 differentially expressed (DE) lncRNAs, 1478, 2643, 404, and 750 DE genes (DEGs), and 74,51, 34, and 23 DE miRNAs (DE miRNAs) were identified in the Y4 vs. Y1, W4 vs. W1, Y4 vs. W4, and Y1 vs. W1 comparisons, respectively. Functional analysis of these DE transcripts (DETs) demonstrated that they had participated in various biological processes, such as cellular component organization or biogenesis, the developmental process, the metabolic process, bone development, and cartilage development. The crucial bone development-related candidate genes NK3 Homeobox 2 (NKX3.2), Wnt ligand secretion mediator (WLS), gremlin 1 (GREM1), fibroblast growth factor receptor 3 (FGFR3), hematopoietically expressed homeobox (HHEX), (collagen type XI alpha 1 chain (COL11A1), and Wnt Family Member 16 (WNT16)) were further identified by functional analysis. Moreover, lncRNA, miRNA, and gene interaction networks were constructed; a total of 55 lncRNAs, 6 miRNAs, and 7 genes formed lncRNA-gene, miRNA-gene, and lncRNA-miRNA-gene pairs, respectively. The aim was to demonstrate that coding and non-coding genes may co-regulate porcine spine development through interaction networks. NKX3.2 was identified as being specifically expressed in cartilage tissues, and it delayed chondrocyte differentiation. miRNA-326 regulated chondrocyte differentiation by targeting NKX3.2. The present study provides the first non-coding RNA and gene expression profiles in the porcine TIC, constructs the lncRNA-miRNA-gene interaction networks, and confirms the function of NKX3.2 in vertebral column development. These findings contribute to the understanding of the potential molecular mechanisms regulating pig vertebral column development. They expand our knowledge about the differences in body length between different pig species and provide a foundation for future studies.

Association of HOXC8 Genetic Polymorphisms with Multi-Vertebral Number and Carcass Weight in Dezhou Donkey

An increase in the number of vertebrae can significantly affect the meat production performance of livestock, thus increasing carcass weight, which is of great importance for livestock production. The homeobox gene C8 (HOXC8) has been identified as an essential candidate gene for regulating vertebral development. However, it has not been researched on the Dezhou donkey. This study aimed to verify the Dezhou donkey HOXC8 gene's polymorphisms and assess their effects on multiple vertebral numbers and carcass weight. In this study, the entire HOXC8 gene of the Dezhou donkey was sequenced, SNPs at the whole gene level were identified, and typing was accomplished utilizing a targeted sequencing genotype detection technique (GBTS). Then, a general linear model was used to perform an association study of HOXC8 gene polymorphism loci, multiple vertebral numbers, and carcass weight for screening candidate markers that can be used for molecular breeding of Dezhou donkeys. These findings revealed that HOXC8 included 12 SNPs, all unique mutant loci. The HOXC8 g.15179224C>T was significantly negatively associated with carcass weight (CW) and lumbar vertebrae length (LL) (p < 0.05). The g.15179674G>A locus was shown to be significantly positively associated with the number of lumbar vertebrae (LN) (p < 0.05). The phylogenetic tree constructed for the Dezhou donkey HOXC8 gene and seven other species revealed that the HOXC8 gene was highly conserved during animal evolution but differed markedly among distantly related animals. The results suggest that HOXC8 is a vital gene affecting multiple vertebral numbers and carcass weight in Dezhou donkeys, and the two loci g.15179224C>T and g.15179674G>A may be potential genetic markers for screening and breeding of new strains of high-quality and high-yielding Dezhou donkeys.

Spinal fixation after laminectomy in pigs prevents postoperative spinal cord injury

Background

A safe, effective, and ethically sound animal model is essential for preclinical research to investigate spinal medical devices. We report the initial failure of a porcine spinal survival model and a potential solution by fixating the spine.

Methods

Eleven female Dutch Landrace pigs underwent a spinal lumbar interlaminar decompression with durotomy and were randomized for implantation of a medical device or control group. Magnetic resonance imaging (MRI) was performed before termination.

Results

Neurological deficits were observed in 6 out of the first 8 animals. Three of these animals were terminated prematurely because they reached the predefined humane endpoint. Spinal cord compression and myelopathy was observed on postoperative MRI imaging. We hypothesized postoperative spinal instability with epidural hematoma, inherent to the biology of the model, and subsequent spinal cord injury as a potential cause. In the subsequent 3 animals, we fixated the spine with Lubra plates. All these animals recovered without neurological deficits. The extent of spinal cord compression on MRI was variable across animals and did not seem to correspond well with neurological outcome.

Conclusion

This study shows that in a porcine in vivo model of interlaminar decompression and durotomy, fixation of the spine after lumbar interlaminar decompression is feasible and may improve neurological outcomes. Additional research is necessary to evaluate this hypothesis.

Genome‑wide association study and genomic prediction for growth traits in yellow-plumage chicken using genotyping-by-sequencing

Background

Growth traits are of great importance for poultry breeding and production and have been the topic of extensive investigation, with many quantitative trait loci (QTL) detected. However, due to their complex genetic background, few causative genes have been confirmed and the underlying molecular mechanisms remain unclear, thus limiting our understanding of QTL and their potential use for the genetic improvement of poultry. Therefore, deciphering the genetic architecture is a promising avenue for optimising genomic prediction strategies and exploiting genomic information for commercial breeding. The objectives of this study were to: (1) conduct a genome-wide association study to identify key genetic factors and explore the polygenicity of chicken growth traits; (2) investigate the efficiency of genomic prediction in broilers; and (3) evaluate genomic predictions that harness genomic features.

Results

We identified five significant QTL, including one on chromosome 4 with major effects and four on chromosomes 1, 2, 17, and 27 with minor effects, accounting for 14.5 to 34.1% and 0.2 to 2.6% of the genomic additive genetic variance, respectively, and 23.3 to 46.7% and 0.6 to 4.5% of the observed predictive accuracy of breeding values, respectively. Further analysis showed that the QTL with minor effects collectively had a considerable influence, reflecting the polygenicity of the genetic background. The accuracy of genomic best linear unbiased predictions (BLUP) was improved by 22.0 to 70.3% compared to that of the conventional pedigree-based BLUP model. The genomic feature BLUP model further improved the observed prediction accuracy by 13.8 to 15.2% compared to the genomic BLUP model.

Conclusions

A major QTL and four minor QTL were identified for growth traits; the remaining variance was due to QTL effects that were too small to be detected. The genomic BLUP and genomic feature BLUP models yielded considerably higher prediction accuracy compared to the pedigree-based BLUP model. This study revealed the polygenicity of growth traits in yellow-plumage chickens and demonstrated that the predictive ability can be greatly improved by using genomic information and related features.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-021-00672-9.

A further survey of the quantitative trait loci affecting swine body size and carcass traits in five related pig populations

Breeding for good meat quality performance while maintaining large body size and desirable carcass traits has been the major challenge for modern swine selective breeding. To address this goal, in the present work we studied five related populations produced by two commercial breeds (Berkshire and Duroc) and two Chinese breeds (Licha black pig and Lulai black pig). A single-trait GWAS performed on 20 body size and carcass traits using a self-developed China Chip-1 porcine SNP50K BeadChip identified 11 genome-wide significant QTL on nine chromosomes and 22 suggestive QTL on 15 chromosomes. For the 11 genome-wide significant QTL, eight were detected in at least two populations, and the rest were population-specific and only mapped in Shanxia black pig. Most of the genome-wide significant QTL were pleiotropic; for example, the QTL around 75.65 Mb on SSC4 was associated with four traits at genome-wide significance level. After screening the genes within 50 kb of the top SNP for each genome-wide significant QTL, NR6A1 and VRTN were chosen as candidate genes for vertebrae number; PLAG1 and BMP2 were identified as candidate genes for body size; and MC4R was the strong candidate gene for body weight. The four genes have been reported as candidates for thoracic vertebrae number, lumbar vertebrae number, carcass length and body weight respectively in previous studies. The effects of VRTN on thoracic vertebrae number, carcass length and body length have been verified in Shanxia black pig. Therefore, the VRTN genotype could be used in gene-assisted selection, and this could accelerate genetic improvement of body size and carcass traits in Shanxia black pig.

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.

Perosomus Elumbis in Piglets: Pathological, Radiological and Cytogenetic Findings

Perosomus elumbis (PE) is a rare congenital condition characterized by agenesis of the lumbar, sacral and coccygeal vertebrae. Perosomus elumbis has rarely been reported in literature as morphological description of singles or few cases. Here we report the first extensive description of eight cases of PE detected in two consecutive litters from the same parents of Casertana pig breed. In August 2018, eight piglets were investigated for multiple malformations. All malformed animals, but one, died in the first day of life. The survivor piglet died at 23 days of age. Pathological, radiological and cytogenetic examination was performed. Furthermore, a farm epidemiological investigation was carried out to investigate the percentage of piglets born dead or with malformations in 2018. The radiological and pathological exams showed skeletal abnormalities at the spinal cord level and visceral malformations. Cytogenetic investigations showed a normal chromosome arrangement. Finally, epidemiological investigation revealed a low prevalence of malformations in newborn pigs, equal to 0.5% of the total birth rate of the farm. Our findings report the first extensive description of PE cases in pigs and suggest an underestimation of this malformation in veterinary medicine. Our findings also suggest a specific genetic etiological basis as cause of PE in pigs and exclude chromosomal abnormalities. Further studies will be performed to confirm this hypothesis.

Combined caudal retrocostal and lateral ultrasound-guided approach for transversus abdominis plane injection: A descriptive pilot study in pig cadavers

Transversus abdominis plane (TAP) block is a regional anesthetic technique used to desensitize the abdominal wall in several species. This study aimed to describe the anatomical characteristics of the abdominal wall and to identify a feasible approach for an US-guided TAP injection that would result in adequate staining of the relevant nerves in the abdominal wall in pig cadavers. Fresh cadavers from five Landrace pigs (age, 12 weeks; body weight, 35.5 ± 1.6 kg) were used. One pig (n = 1) was anatomically dissected, and four pigs (n = 4; i.e., 8 hemiabdomens) were used for TAP injections and evaluation of dye spread. The volume of 0.3 mL/kg/injection point of methylene blue was injected bilaterally. In the caudal retrocostal approach, the injection was performed ventral to the most caudal part of the costal arch. In the lateral approach, the injection was performed between the last rib and iliac crest. A needle was inserted in plane for the caudal retrocostal and the lateral approach caudocranially and craniocaudally, respectively. Successful staining was defined as presence of dye on the nerve for a length of >1 cm in its entire circumference. The TAP was found between different muscle layers in the described anatomical regions. In the caudal retrocostal approach the TAP was found between the external abdominal oblique and transversus abdominis muscle bellies. In the lateral approach the TAP was found between the internal abdominal oblique and transversus abdominis muscles. The approach combining lateral and caudal retrocostal injections at the studied volume stained a median of 5 (3–6) target nerves from the fourth-last thoracic nerve to L2 (six nerves). Combined caudal retrocostal and lateral TAP injections of 0.3 mL/kg/injection point, resulted in staining of target nerve branches which supply the periumbilical and caudal abdominal wall in pig cadavers.

Description of different ultrasound guided transversus abdominis plane block in pig cadavers - A pilot study

OBJECTIVE

To develop an ultrasound (US)-guided transversus abdominis plane (TAP) block in pigs.

STUDY DESIGN

Prospective, descriptive, experimental study.

ANIMALS

A total of 26 hemi-abdomens belonging to 13 Seghers Hybrid cadavers.

METHODS

The study consisted of two parts. In part 1, a preliminary anatomical understanding of the abdominal wall innervation in pigs was established (two hemi-abdomens of one pig). Part 2 was divided into three phases, and a US-guided TAP technique using methylene blue dye (0.3 mL kg^-1 for each hemi-abdomen) was developed. In chronological order, a pilot study (phase I) was conducted to establish good injection points (four hemi-abdomens). In phase II a two-point injection technique (10 hemi-abdomens) was performed. A cranial injection was made at two-thirds of the distance between the xyphoid process and the iliac crest, immediately ventral to the rib arch. A caudal injection was performed ventral to the last rib. In phase III a three-point injection technique was performed (10 hemi-abdomens) with an extra injection point halfway between the cranial and caudal injection point. Staining of the different nerves was recorded during dissection.

RESULTS

In part 1 the nerves innervating the pig's abdominal wall, thoracic vertebrae 12-16 (T12-16) and lumbar vertebrae 1-3 (L1-3), were identified. In part 2, data from the pilot study, adequate nerve staining with the two- or three-point injection technique, were obtained in T12 (29 versus 45% respectively), T13 (29 versus 100%), T14 (73 versus 45%), T15 (75 versus 78%), T16 (58 versus 78%), L1 (100 versus 75%), L2 (88 versus 84%) and L3 (23 versus 0%).

CONCLUSIONS AND CLINICAL RELEVANCE

US-guided TAP block can be used in swine but only a moderate success rate for adequate nerve staining was achieved in this study. Further studies are necessary to determine a correct injection volume and assess clinical intra- and postoperative efficacy.

Image-guided study of swine anatomy as a tool for urologic surgery research and training

PURPOSE

To describe the anatomy of the swine urinary system using computed tomography and to discuss the role of this animal as an experimental model for urological procedures.

METHODS

Three male Landrace pigs underwent computed tomography and the anatomy of the urinary system and renal circulation was analyzed and described.

RESULTS

In all animals, 2 kidneys, 2 ureters and one bladder were identified. Each kidney presented a single renal artery vascularization, with a mean diameter on the right of 4.45 and 5.31 mm on the left (p < 0.0001) and single renal vein drainage, with a mean diameter on the right of 5.78 and 5.82 mm on the left (p = 0.0336). The average renal length was 9.85 cm on the right and 10.30 cm on the left (p < 0.0001). The average renal volume was 113.70 cm3 on the right and 109.70 cm3 on the left (p < 0.0001). The average length of the ureter was 19.78 cm on the right and 22.08 cm on the left (p < 0.0001). The average bladder volume was 423.70 cm3.

CONCLUSIONS

The data obtained show similarities with human anatomy, suggesting the viability of the swine model for planning preclinical trials, basic research, refinement in experimental surgery and surgical training for urological procedures.

A multi-breed GWAS for morphometric traits in four Beninese indigenous cattle breeds reveals loci associated with conformation, carcass and adaptive traits

BACKGROUND

Specific adaptive features including disease resistance and growth abilities in harsh environments are attributed to indigenous cattle breeds of Benin, but these breeds are endangered due to crossbreeding. So far, there is a lack of systematic trait recording, being the basis for breed characterizations, and for structured breeding program designs aiming on conservation. Bridging this gap, own phenotyping for morphological traits considered measurements for height at withers (HAW), sacrum height (SH), heart girth (HG), hip width (HW), body length (BL) and ear length (EL), including 449 cattle from the four indigenous Benin breeds Lagune, Somba, Borgou and Pabli. In order to utilize recent genomic tools for breed characterizations and genetic evaluations, phenotypes for novel traits were merged with high-density SNP marker data. Multi-breed genetic parameter estimations and genome-wide association studies (GWAS) for the six morphometric traits were carried out. Continuatively, we aimed on inferring genomic regions and functional loci potentially associated with conformation, carcass and adaptive traits.

RESULTS

SNP-based heritability estimates for the morphometric traits ranged between 0.46 ± 0.14 (HG) and 0.74 ± 0.13 (HW). Phenotypic and genetic correlations ranged from 0.25 ± 0.05 (HW-BL) to 0.89 ± 0.01 (HAW-SH), and from 0.14 ± 0.10 (HW-BL) to 0.85 ± 0.02 (HAW-SH), respectively. Three genome-wide and 25 chromosome-wide significant SNP positioned on different chromosomes were detected, located in very close chromosomal distance (±25 kb) to 15 genes (or located within the genes). The genes PIK3R6 and PIK3R1 showed direct functional associations with height and body size. We inferred the potential candidate genes VEPH1, CNTNAP5, GYPC for conformation, growth and carcass traits including body weight and body fat deposition. According to their functional annotations, detected potential candidate genes were associated with stress or immune response (genes PTAFR, PBRM1, ADAMTS12) and with feed efficiency (genes MEGF11 SLC16A4, CCDC117).

CONCLUSIONS

Accurate measurements contributed to large SNP heritabilities for some morphological traits, even for a small mixed-breed sample size. Multi-breed GWAS detected different loci associated with conformation or carcass traits. The identified potential candidate genes for immune response or feed efficiency indicators reflect the evolutionary development and adaptability features of the breeds.

Polymorphism of VRTN Gene g.20311_20312ins291 Was Associated with the Number of Ribs, Carcass Diagonal Length and Cannon Bone Circumference in Suhuai Pigs

Simple Summary

An increase in the number of ribs (RIB) could improve carcass length (CL) and body size. Cannon bone circumference (CBC) is a pivotal body size trait, and a large CBC could enhance the capacity to bear excessive body weight, vigorous exercise, and resistance to injuries. Several researchers showed that the vertnin (VRTN) gene g.20311_20312ins291 (NC_010449.5 7: g.20311_20312ins291) is an important variant that is related to RIB and CL of Western pigs. However, it is unknown whether this variant could affect the CBC of pigs. Our study showed that this variant was significantly associated with RIB, carcass diagonal length (CDL), and CBC in Suhuai pigs; therefore, it could be used as a potential molecular marker for improving RIB, CDL, and CBC in this breed.

Abstract

The vertnin (VRTN) gene g.20311_20312ins291 was reported as an important variant related to the number of ribs (RIB), and the ins/ins genotype was advantageous for improving RIB of Western pigs. The purpose of this study was to determine whether the VRTN gene g.20311_20312ins291 influences RIB, carcass traits, and body size traits, including cannon bone circumference (CBC) in Chinese Suhuai pigs. We found that the VRTN gene g.20311_20312ins291 was polymorphic in Suhuai fattening pigs and gilts. The polymorphism of g.20311_20312ins291 was significantly associated with RIB and CDL in Suhuai fattening pigs (p < 0.01), whereas this variant had no influence on carcass weight (CWT). There was a tendency of association between this variant and carcass straight length (CSL) in Suhuai fattening pigs (p = 0.06). The polymorphism of g.20311_20312ins291 was also significantly associated with CBC in Suhuai gilts (p = 0.04). Furthermore, CBC was positively genetically correlated with body length (0.22, p < 0.01) and body weight (0.15, p < 0.01). Our results indicated that the VRTN gene g.20311_20312ins291 could be used as a potential marker for improving RIB, CDL, and CBC in Suhuai pigs.

Insight into unique somitogenesis of yak (Bos grunniens) with one additional thoracic vertebra

Background

The yak is a species of livestock which is crucial for local communities of the Qinghai-Tibet Plateau and adjacent regions and naturally owns one more thoracic vertebra than cattle. Recently, a sub-population of yak termed as the Jinchuan yak has been identified with over half its members own a thoracolumbar vertebral formula of T15L5 instead of the natural T14L5 arrangement. The novel T15L5 positioning is a preferred genetic trait leading to enhanced meat and milk production. Selective breeding of this trait would have great agricultural value and exploration of the molecular mechanisms underlying this trait would both accelerate this process and provide us insight into the development and regulation of somitogenesis.

Results

Here we investigated the genetic background of the Jinchuan yak through resequencing fifteen individuals, comprising five T15L5 individuals and ten T14L5 individuals with an average sequencing depth of > 10X, whose thoracolumbar vertebral formulae were confirmed by anatomical observation. Principal component analysis, linkage disequilibrium analysis, phylogenetic analysis, and selective sweep analysis were carried out to explore Jinchuan yak’s genetic background. Three hundred and thirty candidate markers were identified as associated with the additional thoracic vertebrae and target sequencing was used to validate seven carefully selected markers in an additional 51 Jinchuan yaks. The accuracies of predicting 15 thoracic vertebrae and 20 thoracolumbar vertebrae with these 7 markers were 100.00 and 33.33% despite they both could only represent 20% of all possible genetic diversity. Two genes, PPP2R2B and TBLR1, were found to harbour the most candidate markers associated with the trait and likely contribute to the unique somitic number and identity according to their reported roles in the mechanism of somitogenesis.

Conclusions

Our findings provide a clear depiction of the Jinchuan yak’s genetic background and a solid foundation for marker-assistant selection. Further exploitation of this unique population and trait could be promoted with the aid of our genomic resource.

GWAS for Meat and Carcass Traits Using Imputed Sequence Level Genotypes in Pooled F2-Designs in Pigs

In order to gain insight into the genetic architecture of economically important traits in pigs and to derive suitable genetic markers to improve these traits in breeding programs, many studies have been conducted to map quantitative trait loci. Shortcomings of these studies were low mapping resolution, large confidence intervals for quantitative trait loci-positions and large linkage disequilibrium blocks. Here, we overcome these shortcomings by pooling four large F2 designs to produce smaller linkage disequilibrium blocks and by resequencing the founder generation at high coverage and the F1 generation at low coverage for subsequent imputation of the F2 generation to whole genome sequencing marker density. This lead to the discovery of more than 32 million variants, 8 million of which have not been previously reported. The pooling of the four F2 designs enabled us to perform a joint genome-wide association study, which lead to the identification of numerous significantly associated variant clusters on chromosomes 1, 2, 4, 7, 17 and 18 for the growth and carcass traits average daily gain, back fat thickness, meat fat ratio, and carcass length. We could not only confirm previously reported, but also discovered new quantitative trait loci. As a result, several new candidate genes are discussed, among them BMP2 (bone morphogenetic protein 2), which we recently discovered in a related study. Variant effect prediction revealed that 15 high impact variants for the traits back fat thickness, meat fat ratio and carcass length were among the statistically significantly associated variants.

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.

Whole-Genome Resequencing Reveals Loci Associated With Thoracic Vertebrae Number in Sheep

The number of vertebrae, especially thoracic vertebrae, is an important economic trait that may influence carcass length and meat production in animals. However, the genetic basis of vertebrae number in sheep is still poorly understood. To detect the candidate genes, 400 increased number of thoracic vertebrae (T14L6) and 200 normal (T13L6) Kazakh sheep were collected. We generated and sequenced 60 pools of genomic DNA (each pool prepared by mixing genomic DNA from 10 sheep with the same thoracic traits), with an average depth of coverage of 25.65×. We identified a total of 42,075,402 SNPs and 11 putatively selected genomic regions, including the VRTN gene and the HoxA gene family that regulate vertebral development. The most prominent areas of selective elimination were located in a region of chromosome 7, including VRTN, which regulates spinal development and morphology. Further investigation indicated that the expression level of the VRTN gene during fetal development was significantly higher in sheep with more thoracic vertebrae than in those with a normal number of thoracic vertebrae. A genome-wide comparison between sheep with increased and normal numbers of thoracic vertebrae showed that the VRTN gene is the major selection locus for the number of thoracic vertebrae in sheep and has the potential to be utilized in sheep breeding in the future.

Validation of a methodology for neuro-urological and lumbosacral stimulation studies in domestic pigs: a humanlike animal model

OBJECTIVE

Spinal cord injuries (SCIs) result in loss of movement and sensory feedback, but also organ dysfunction. Nearly all patients with complete SCI lose bladder control and are prone to kidney failure if intermittent catheterization is not performed. Electrical stimulation of sacral spinal roots was initially considered to be a promising approach for restoring continence and micturition control, but many patients are discouraged by the need for surgical deafferentation as it could lead to a loss of sensory functions and reflexes. Nevertheless, recent research findings highlight the renewed interest in spinal cord stimulation (SCS). It is thought that synergic recruitment of spinal fibers could be achieved by stimulating the spinal neural networks involved in regulating physiological processes. Paradoxically, most of these recent studies focused on locomotor issues, while few addressed visceral dysfunction. This could at least partially be attributed to the lack of methodological tools. In this study, the authors aim to fill this gap by presenting a comprehensive method for investigating the potential of SCS to restore visceral functions in domestic pigs, a large-animal model considered to be a close approximation to humans.

METHODS

This methodology was tested in 7 female pigs (Landrace pig breed, 45-60 kg, 4 months old) during acute experiments. A combination of morphine and propofol was used for anesthesia when transurethral catheterization and lumbosacral laminectomy (L4-S4) were performed. At the end of the operation, spinal root stimulation (L6-S5) and urodynamic recordings were performed to compare the evoked responses with those observed intraoperatively in humans.

RESULTS

Nervous excitability was preserved despite long-term anesthesia (mean 8.43 ± 1.5 hours). Transurethral catheterization and conventional laminectomy were possible while motor responses (gluteus muscle monitoring) were unaffected throughout the procedure. Consistent detrusor (approximately 25 cm H2O) and sphincter responses were obtained, whereas spinal root stimulation elicited detrusor and external urethral sphincter co-contractions similar to those observed intraoperatively in humans.

CONCLUSIONS

Pigs represent an ideal model for SCS studies aimed at visceral function investigation and restoration because of the close similarities between female domestic pigs and humans, both in terms of anatomical structure and experimental techniques implemented. This article provides methodological keys for conducting experiments with equipment routinely used in clinical practice.

Distribution, Morphological Characterization, and Resiniferatoxin-Susceptibility of Sensory Neurons That Innervate Rat Perirenal Adipose Tissue

Perirenal adipose tissue (PrAT) is a visceral adipose tissue involved in the pathogenesis of obesity and cardiovascular diseases via neural pathways. However, the origins, morphological characterization, and resiniferatoxin (RTX)-susceptibility of sensory neurons that innervate rat PrAT are yet unclear. Using neural tracing, an injection of DiI (1,1′-dioctadecyl-3,3,3′,3′-tetramethylindocarbocyanine perchlorate) into PrAT revealed that sensory neurons that innervate PrAT reside in T9-L3 dorsal root ganglia (DRG). Peak labeling occurred in T13 and L1 DRGs. Two distinct peaks were observed in cross-sectional areas of the labeled soma, and the mean cross-sectional area was 717.1 ± 27.7 μm2. Immunofluorescence staining for transient receptor potential cation channel subfamily V member 1 (TRPV1) separated DiI-positive neurons into three subpopulations: small TRPV1-negative, small TRPV1-positive, and large TRPV1-negative. Furthermore, the injection of RTX into PrAT reduced labeled cells by 36.7% where TRPV1-positive cells were the main target of RTX denervation. These novel findings provide a structural basis for future TRPV1-dependent and TRPV1-independent studies on the sensory innervation of PrAT, which may be of interest for future therapeutic obesity treatment and intervention.

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.

Genomic analysis reveals genes affecting distinct phenotypes among different Chinese and western pig breeds

The differences in artificial and natural selection have been some of the factors contributing to phenotypic diversity between Chinese and western pigs. Here, 830 individuals from western and Chinese pig breeds were genotyped using the reduced-representation genotyping method. First, we identified the selection signatures for different pig breeds. By comparing Chinese pigs and western pigs along the first principal component, the growth gene IGF1R; the immune genes IL1R1, IL1RL1, DUSP10, RAC3 and SWAP70; the meat quality-related gene SNORA50 and the olfactory gene OR1F1 were identified as candidate differentiated targets. Further, along a principal component separating Pudong White pigs from others, a potential causal gene for coat colour (EDNRB) was discovered. In addition, the divergent signatures evaluated by F_st within Chinese pig breeds found genes associated with the phenotypic features of coat colour, meat quality and feed efficiency among these indigenous pigs. Second, admixture and genomic introgression analysis were performed. Shan pigs have introgressed genes from Berkshire, Yorkshire and Hongdenglong pigs. The results of introgression mapping showed that this introgression conferred adaption to the local environment and coat colour of Chinese pigs and the superior productivity of western pigs.

VRTN is Required for the Development of Thoracic Vertebrae in Mammals

Vertnin (VRTN) variants are associated with thoracic vertebral number (TVN) in pigs. However, the biological function of VRTN remains poorly understood. Here we first conducted a range of experiments to demonstrate that VRTN is a responsible gene for TVN and two causative variants in the regulatory region of VRTN additively regulate TVN. Then, we show that VRTN is a novel DNA-binding transcription factor as it localizes exclusively in the nucleus, binds to DNA on a genome-wide scale and regulates the transcription of a set of genes that harbor VRTN binding motifs. Next, we illustrate that VRTN is essential for the development of thoracic vertebrae. Vrtn-null embryos display somitogenesis defect with the failure of axial rotation and fewer somites at the thoracic somite stage. Half of Vrtn heterozygous mice show abnormal spinal development with fewer thoracic vertebrae and ribs than their wild-type littermates. Lastly, we reveal that VRTN could modulate somite segmentation via the Notch signaling pathway. The findings advance our understanding of the mechanisms underlying the development of thoracic vertebrate in mammals, and VRTN causative variants provide a robust tool to improve pork production by selecting the alleles increasing the number of thoracic vertebrae and ribs.

Genome-wide association study in an F2 Duroc x Pietrain resource population for economically important meat quality and carcass traits

Meat quality is essential for consumer acceptance, it ultimately impacts pork production profitability and it is subject to genetic control. The objective of this study was to map genomic regions associated with economically important meat quality and carcass traits. We performed a genome-wide association (GWA) analysis to map regions associated with 38 meat quality and carcass traits recorded for 948 F2 pigs from the Michigan State University Duroc × Pietrain resource population. The F0, F1, and 336 F2 pigs were genotyped with the Illumina Porcine SNP60 BeadChip, while the remaining F2 pigs were genotyped with the GeneSeek Genomic Profiler for Porcine Low Desnisty (LD) chip, and imputed with high accuracy ( = 0.97). Altogether the genomic dataset comprised 1,019 animals and 44,911 SNP. A Gaussian linear mixed model was fitted to estimate the breeding values and the variance components. A linear transformation was performed to estimate the marker effects and variances. Type I error rate was controlled at a False Discovery Rate of 5%. Seven putative QTL found in this study were previously reported in other studies. Two novel QTL associated with tenderness (TEN) were located on SSC3 [135.6:137.5Mb; False Discovery rate (FDR) < 0.03] and SSC5 (67.3:69.1Mb; FDR < 0.02). The QTL region identified on SSC15 includes Protein Kinase AMP-activated ɣ 3-subunit gene (), which has been associated with 24-h pH (pH24), drip loss (DL) and cook yield (CY). Also, novel candidate genes were identified for TEN in the region on SSC5 [A Kinase (PRKA) Anchor Protein 3 (], and for tenth rib backfat thickness (BF10) [Carnitine O-Acetyltransferase ()] on SSC1. The association of gene polymorphisms with pork quality traits has been reported for several pig populations. However, there are no SNP for this gene on the chip used, thus we genotyped the animals for 2 non-synonymous variants ( and ). We then performed a GWA conditioning on the genotype of both SNP and was associated with pH24, DL, protein content (PRO) and CY ( < 0.004) and T30N with Juiciness, TEN, shear force, pH24, PRO, and CY < 0.04). Finally, we performed a GWA conditioning on the genotype of the SNP peak detected in this study, and T30N remained associated only with PRO ( < 0.02). Therefore, in this study we identified 2 novel QTL regions, suggest 2 novel candidate genes, and conclude that other SNP in PRKAG3 or nearby gene(s) explain the observed associations on SSC15 in this population.

Impact of direct epispinal stimulation on bladder and bowel functions in pigs: A feasibility study

AIMS

This study assesses the potential of epispinal (subdural) stimulation application in the treatment of urinary and bowel neurological disorders. Acute experiments were performed on a large animal model - the domestic pig - to develop a new methodology facilitating future results and technology transfers to human.

METHODS

After rectal and bladder catheterization, four Landrace pigs (45-50 kg) underwent spinal cord surgery - that is, lumbosacral incision, laminectomy [L4-S4], dural opening and microsurgical arachnoid dissection. Three successive electrical stimulation sessions were carried out: 1) nerve roots stimulation, 2) epispinal stimulation with a matrix electrode, 3) epispinal stimulation with a small diameter needle electrode. Changes in rectal and bladder pressures were monitored throughout the various procedures to identify spinal areas inducing responses while evaluating the influence of electrode contacts size in the measured responses amplitudes.

RESULTS

An interesting area was identified in the upper portion of the spinal myelomeres (ie, spinal cord segment delimited by two successive pairs of spinal roots) directly adjoining root with best pressures (either rectal or vesical). Significant responses (up to 40 cmH₂ O) were also obtained with a needle electrode. Furthermore, bowel evacuation was triggered in one of the animals. Despite the use of smaller electrode contacts, no detrusor or rectum selective responses were observed in none of the sessions.

CONCLUSION

This study showed, for the first time, that epispinal stimulation causes significant detrusor and rectal responses in pigs and allows considering further studies with the objective of treating urinary and rectal disorders in spinal cord injury patients.

Effects of vertebral number variations on carcass traits and genotyping of Vertnin candidate gene in Kazakh sheep

Objective

The vertebral number is associated with body length and carcass traits, which represents an economically important trait in farm animals. The variation of vertebral number has been observed in a few mammalian species. However, the variation of vertebral number and quantitative trait loci in sheep breeds have not been well addressed.

Methods

In our investigation, the information including gender, age, carcass weight, carcass length and the number of thoracic and lumbar vertebrae from 624 China Kazakh sheep was collected. The effect of vertebral number variation on carcass weight and carcass length was estimated by general linear model. Further, the polymorphic sites of Vertnin (VRTN) gene were identified by sequencing, and the association of the genotype and vertebral number variation was analyzed by the one-way analysis of variance model.

Results

The variation of thoracolumbar vertebrae number in Kazakh sheep (18 to 20) was smaller than that in Texel sheep (17 to 21). The individuals with 19 thoracolumbar vertebrae (T13L6) were dominant in Kazakh sheep (79.2%). The association study showed that the numbers of thoracolumbar vertebrae were positively correlated with the carcass length and carcass weight, statistically significant with carcass length. To investigate the association of thoracolumbar vertebrae number with VRTN gene, we genotyped the VRTN gene. A total of 9 polymorphic sites were detected and only a single nucleotide polymorphism (SNP) (rs426367238) was suggested to associate with thoracic vertebral number statistically.

Conclusion

The variation of thoracolumbar vertebrae number positively associated with the carcass length and carcass weight, especially with the carcass length. VRTN gene polymorphism of the SNP (rs426367238) with significant effect on thoracic vertebral number could be as a candidate marker to further evaluate its role in influence of thoracolumbar vertebral number.

Genome-wide association study for conformation traits in three Danish pig breeds

Background

Selection for sound conformation has been widely used as a primary approach to reduce lameness and leg weakness in pigs. Identification of genomic regions that affect conformation traits would help to improve selection accuracy for these lowly to moderately heritable traits. Our objective was to identify genetic factors that underlie leg and back conformation traits in three Danish pig breeds by performing a genome-wide association study followed by meta-analyses.

Methods

Data on four conformation traits (front leg, back, hind leg and overall conformation) for three Danish pig breeds (23,898 Landrace, 24,130 Yorkshire and 16,524 Duroc pigs) were used for association analyses. Estimated effects of single nucleotide polymorphisms (SNPs) from single-trait association analyses were combined in two meta-analyses: (1) a within-breed meta-analysis for multiple traits to examine if there are pleiotropic genetic variants within a breed; and (2) an across-breed meta-analysis for a single trait to examine if the same quantitative trait loci (QTL) segregate across breeds. SNP annotation was implemented through Sus scrofa Build 10.2 on Ensembl to search for candidate genes.

Results

Among the 14, 12 and 13 QTL that were detected in the single-trait association analyses for the three breeds, the most significant SNPs explained 2, 2.3 and 11.4% of genetic variance for back quality in Landrace, overall conformation in Yorkshire and back quality in Duroc, respectively. Several candidate genes for these QTL were also identified, i.e. LRPPRC, WRAP73, VRTN and PPARD likely control conformation traits through the regulation of bone and muscle development, and IGF2BP2, GH1, CCND2 and MSH2 can have an influence through growth-related processes. Meta-analyses not only confirmed many significant SNPs from single-trait analyses with higher significance levels, but also detected several additional associated SNPs and suggested QTL with possible pleiotropic effects.

Conclusions

Our results imply that conformation traits are complex and may be partly controlled by genes that are involved in bone and skeleton development, muscle and fat metabolism, and growth processes. A reliable list of QTL and candidate genes was provided that can be used in fine-mapping and marker assisted selection to improve conformation traits in pigs.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-017-0289-2) contains supplementary material, which is available to authorized users.

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

Background

Number of functional teats is an important trait in commercial swine production. As litter size increases, the number of teats must also increase to supply nutrition to all piglets. Therefore, a genome-wide association analysis was conducted to identify genomic regions that affect this trait in a commercial swine population. Genotypic data from the Illumina Porcine SNP60v1 BeadChip were available for 2951 animals with total teat number (TTN) records. A subset of these animals (n = 1828) had number of teats on each side recorded. From this information, the following traits were derived: number of teats on the left (LTN) and right side (RTN), maximum number of teats on a side (MAX), difference between LTN and RTN (L − R) and absolute value of L − R (DIF). Bayes C option of GENSEL (version 4.61) and 1-Mb windows were implemented. Identified regions that explained more than 1.5% of the genomic variation were tested in a larger group of animals (n = 5453) to estimate additive genetic effects.

Results

Marker heritabilities were highest for TTN (0.233), intermediate for individual side counts (0.088 to 0.115) and virtually nil for difference traits (0.002 for L − R and 0.006 for DIF). Each copy of the VRTN mutant allele increased teat count by 0.35 (TTN), 0.16 (LTN and RTN) and 0.19 (MAX). 15, 18, 13 and 18 one-Mb windows were detected that explained more than 1.0% of the genomic variation for TTN, LTN, RTN, and MAX, respectively. These regions cumulatively accounted for over 50% of the genomic variation of LTN, RTN and MAX, but only 30% of that of TTN. Sus scrofa chromosome SSC10:52 Mb was associated with all four count traits, while SSC10:60 and SSC14:54 Mb were associated with three count traits. Thirty-three SNPs accounted for nearly 39% of the additive genetic variation in the validation dataset. No effect of piglet sex or percentage of males in litter was detected, but birth weight was positively correlated with TTN.

Conclusions

Teat number is a heritable trait and use of genetic markers would expedite selection progress. Exploiting genetic variation associated with teat counts on each side would enhance selection focused on total teat counts. These results confirm QTL on SSC4, seven and ten and identify a novel QTL on SSC14.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-016-0282-1) contains supplementary material, which is available to authorized users.