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

ABCD4 is associated with mammary gland development in mammals

BACKGROUND

Mammary gland development is a critical process in mammals, crucial for their reproductive success and offspring nourishment. However, the functional roles of key candidate genes associated with teat number, including ABCD4, VRTN, PROX2, and DLST, in this developmental process remain elusive. To address this gap in knowledge, we conducted an in-depth investigation into the dynamic expression patterns, functional implications, and regulatory networks of these candidate genes during mouse mammary gland development.

RESULTS

In this study, the spatial and temporal patterns of key genes were characterized in mammary gland development. Using time-series single-cell data, we uncovered differences in the expression of A bcd4, Vrtn, Prox2, and Dlst in cell population of the mammary gland during embryonic and adult stages, while Vrtn was not detected in any cells. We found that only overexpression and knockdown of Abcd4 could inhibit proliferation and promote apoptosis of HC11 mammary epithelial cells, whereas Prox2 and Dlst had no significant effect on these cells. Using RNA-seq and qPCR, further analysis revealed that Abcd4 can induce widespread changes in the expression levels of genes involved in mammary gland development, such as Igfbp3, Ccl5, Tlr2, and Prlr, which were primarily associated with the MAPK, JAK-STAT, and PI3K-AKT pathways by functional enrichment.

CONCLUSIONS

These findings revealed ABCD4 as a candidate gene pivotal for regulating mammary gland development and lactation during pregnancy by influencing PRLR expression.

Detection and Analysis of Genes Affecting the Number of Thoracic Vertebrae in Licha Black Pigs

The number of thoracic vertebrae (NTV) in pigs is an important economic trait that significantly influences pork production. While the Licha black pig is a well-known Chinese pig breed with multiple thoracic vertebrae, the genetic mechanism is still unknown. Here, we performed a selective signal analysis on the genome of Licha black pigs, comparing individuals with 15 NTV versus those with 16 NTV to better understand functional genes associated with NTV. A total of 2265 selection signal regions were detected across the genome, including 527 genes and 1073 QTL that overlapped with the selection signal regions. Functional enrichment analysis revealed that LRP5 and SP5 genes were involved in biological processes such as bone morphogenesis and Wnt protein binding. Furthermore, three genes, LRP8, DEF6, and SCUBE3, associated with osteoblast differentiation and bone formation, were located within or close to the QTL related to bone development and vertebrae number. These five genes were hypothesized to be potential candidates for regulating the NTV trait in Licha black pigs. Our findings revealed several candidate genes that play crucial roles in NTV regulation and provide a theoretical foundation to understand the genetic mechanism of the NTV trait in pig breeding.

Effects of reference population size and structure on genomic prediction of maternal traits in two pig lines using whole-genome sequence-, high-density- and combined annotation-dependent depletion genotypes

The aim of this study was to investigate the reference population size required to obtain substantial prediction accuracy within- and across-lines and the effect of using a multi-line reference population for genomic predictions of maternal traits in pigs. The data consisted of two nucleus pig populations, one pure-bred Landrace (L) and one Synthetic (S) Yorkshire/Large White line. All animals were genotyped with up to 30 K animals in each line, and all had records on maternal traits. Prediction accuracy was tested with three different marker data sets: High-density SNP (HD), whole genome sequence (WGS), and markers derived from WGS based on pig combined annotation dependent depletion-score (pCADD). Also, two different genomic prediction methods (GBLUP and Bayes GC) were compared for four maternal traits; total number piglets born (TNB), total number of stillborn piglets (STB), Shoulder Lesion Score and Body Condition Score. The main results from this study showed that a reference population of 3 K-6 K animals for within-line prediction generally was sufficient to achieve high prediction accuracy. However, when the number of animals in the reference population was increased to 30 K, the prediction accuracy significantly increased for the traits TNB and STB. For multi-line prediction accuracy, the accuracy was most dependent on the number of within-line animals in the reference data. The S-line provided a generally higher prediction accuracy compared to the L-line. Using pCADD scores to reduce the number of markers from WGS data in combination with the GBLUP method generally reduced prediction accuracies relative to GBLUP using HD genotypes. The BayesGC method benefited from a large reference population and was less dependent on the different genotype marker datasets to achieve a high prediction accuracy.

Genomic background of biotypes related to growth, carcass and meat quality traits in Duroc pigs based on principal component analysis

As the swine industry continues to explore pork quality traits alongside growth, feed efficiency and carcass leanness traits, it becomes imperative to understand their underlying genetic relationships. Due to this increase in the number of desirable traits, animal breeders must also consider methods to efficiently perform direct genetic changes for each trait and evaluate alternative selection indexes with different sets of phenotypic measurements. Principal component analysis (PCA) and genome-wide association studies (GWAS) can be combined to understand the genetic architecture and biological mechanisms by defining biological types (biotypes) that relate these valuable traits. Therefore, the main objectives of this study were to: (1) estimate genomic-based genetic parameters; (2) define animal biotypes utilizing PCA; and (3) utilize GWAS to link the biotypes to candidate genes and quantitative trait loci (QTL). The phenotypic dataset included 2583 phenotypic records from female Duroc pigs from a terminal sire line. The pedigree file contained 193,764 animals and the genotype file included 21,309 animals with 35,651 single nucleotide polymorphisms (SNPs). Eight principal components (PCs), accounting for a total of 99.7% of the population variation, were defined for three growth, eight conventional carcass, 10 pork quality and 18 novel carcass traits. The eight biotypes defined from the PCs were found to be related to growth rate, maturity, meat quality and body structure, which were then related to candidate genes. Of the 175 candidate genes found, six of them [LDHA (SSC1), PIK3C3 (SSC6), PRKAG3 (SSC15), VRTN (SSC7), DLST (SSC7) and PAPPA (SSC1)] related to four PCs were found to be associated with previously defined QTL, linking the biotypes with biological processes involved with muscle growth, fat deposition, glycogen levels and skeletal development. Further functional analyses helped to make connections between biotypes, relating them through common KEGG pathways and gene ontology (GO) terms. These findings contribute to a better understanding of the genetic relationships between growth, carcass and meat quality traits in Duroc pigs, enabling breeders to better understand the biological mechanisms underlying the phenotypic expression of these traits.

Whole-Genome Resequencing of Ujimqin Sheep Identifies Genes Associated with Vertebral Number

The number of vertebrae is a crucial economic trait that can significantly impact the carcass length and meat production in animals. However, our understanding of the quantitative trait loci (QTLs) and candidate genes associated with the vertebral number in sheep (Ovis aries) remains limited. To identify these candidate genes and QTLs, we collected 73 Ujimqin sheep with increased numbers of vertebrae (T13L7, T14L6, and T14L7) and 23 sheep with normal numbers of vertebrae (T13L6). Through high-throughput genome resequencing, we obtained a total of 24,130,801 effective single-nucleotide polymorphisms (SNPs). By conducting a selective-sweep analysis, we discovered that the most significantly selective region was located on chromosome 7. Within this region, we identified several genes, including VRTN, SYNDIG1L, LTBP2, and ABCD4, known to regulate the spinal development and morphology. Further, a genome-wide association study (GWAS) performed on sheep with increased and normal vertebral numbers confirmed that ABCD4 is a candidate gene for determining the number of vertebrae in sheep. Additionally, the most significant SNP on chromosome 7 was identified as a candidate QTL. Moreover, we detected two missense mutations in the ABCD4 gene; one of these mutations (Chr7: 89393414, C > T) at position 22 leads to the conversion of arginine (Arg) to glutamine (Gln), which is expected to negatively affect the protein's function. Notably, a transcriptome expression profile in mouse embryonic development revealed that ABCD4 is highly expressed during the critical period of vertebral formation (4.5-7.5 days). Our study highlights ABCD4 as a potential major gene influencing the number of vertebrae in Ujimqin sheep, with promising prospects for future genome-assisted breeding improvements in sheep.

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.

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.

Young SINEs in pig genomes impact gene regulation, genetic diversity, and complex traits

Transposable elements (TEs) are a major source of genetic polymorphisms and play a role in chromatin architecture, gene regulatory networks, and genomic evolution. However, their functional role in pigs and contributions to complex traits are largely unknown. We created a catalog of TEs (n = 3,087,929) in pigs and found that young SINEs were predominantly silenced by histone modifications, DNA methylation, and decreased accessibility. However, some transcripts from active young SINEs showed high tissue-specificity, as confirmed by analyzing 3570 RNA-seq samples. We also detected 211,067 dimorphic SINEs in 374 individuals, including 340 population-specific ones associated with local adaptation. Mapping these dimorphic SINEs to genome-wide associations of 97 complex traits in pigs, we found 54 candidate genes (e.g., ANK2 and VRTN) that might be mediated by TEs. Our findings highlight the important roles of young SINEs and provide a supplement for genotype-to-phenotype associations and modern breeding in pigs.

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.

Genome-wide association study reveals new QTL and functional candidate genes for the number of ribs and carcass length in pigs

The number of ribs (NR) and carcass length (CL) are important economic traits in pig breeding programs. Pigs with a higher NR and longer CL produce greater pork yields. In the present study, Suhuai pigs with NR and CL phenotypes were genotyped using the Neogen^® GGP Porcine 80 K SNP array to identify the QTL affecting NR and CL and dissect the candidate genes for the two traits. The SNP-chip data was imputed to the whole-genome sequence (iWGS) to increase the probability of identifying causal variants. Through genome-wide association studies (GWAS) based on both chip and iWGS data, significant SNPs were detected on Sus scrofa chromosome (SSC) 1, SSC4 and SSC7 for NR and on SSC5, SSC16 and SSC17 for CL. Moreover, two SNPs (H3GA0022644 and WU_10.2_7_103460706) on SSC7 detected in chip-based GWAS were significantly associated with both NR and CL. Through Bayes fine mapping, one reported QTL for NR on SSC7 and two reported QTL for CL on SSC17 were verified, and two new QTL (SSC1: 14.05-15.84 Mb and SSC4: 64.83-66.59 Mb) affecting NR and two new QTL (SSC5: 58.31-59.84 Mb and SSC16: 22.98-23.43 Mb) affecting CL were detected. According to the biological functions of genes, MTHFD1L on SSC1 and SULF1 on SSC4 are novel functional candidate genes for NR, and EMP1 on SSC5 and EGFLAM on SSC16 are novel functional candidate genes for CL. Overall, our findings provide a basis for identifying new causal genes and mutations affecting NR and CL.

A high-throughput study of visceral organs in CT-scanned pigs

It has been debated whether intensive selection for growth and carcass yield in pig breeding programmes can affect the size of internal organs, and thereby reduce the animal’s ability to handle stress and increase the risk of sudden deaths. To explore the respiratory and circulatory system in pigs, a deep learning based computational pipeline was built to extract the size of lungs and hearts from CT-scan images. This pipeline was applied on CT images from 11,000 boar selection candidates acquired during the last decade. Further, heart and lung volumes were analysed genetically and correlated with production traits. Both heart and lung volumes were heritable, with h² estimated to 0.35 and 0.34, respectively, in Landrace, and 0.28 and 0.4 in Duroc. Both volumes were positively correlated with lean meat percentage, and lung volume was negatively genetically correlated with growth (r_g = − 0.48 ± 0.07 for Landrace and r_g = − 0.44 ± 0.07 for Duroc). The main findings suggest that the current pig breeding programs could, as an indirect response to selection, affect the size of hearts- and lungs. The presented methods can be used to monitor the development of internal organs in the future.

Multi-Thoracolumbar Variations and NR6A1 Gene Polymorphisms Potentially Associated with Body Size and Carcass Traits of Dezhou Donkey

The number of thoracolumbar vertebrae is a quantitative trait positively correlated with the economic traits of livestock. More thoracolumbar vertebrae individuals could genetically be used to improve the livestock population, as more thoracolumbar vertebrae means a longer carcass, which could bring more meat production. Nuclear receptor subfamily 6 group A member 1 (NR6A1) is considered a strong candidate gene for effecting the number of vertebrae in livestock. The purposes of this study are as follows: (a) Analyzing the effect of TLN variation on body size and carcass traits of Dezhou donkey; (b) Studying the distribution of seven single nucleotide variants (SNVs) in NR6A1 gene of Dezhou donkey; (c) Exploring the relationship between latent SNVs and TLN, the body size and carcass traits. We examined the thoracic and lumbar vertebrae number and seven SNVs in NR6A1 gene of 455 Dezhou donkeys, and analyzed the relationships between them. Five types of thoracolumbar combinations (T17L5 (individual with 17 thoracic and five lumbar vertebrae) 2.4%, T18L5 75.8%, T19L5 1.1%, T17L6 11.9%, and T18L6 8.8%) of Dezhou donkeys were detected in this study. For one thoracolumbar vertebra added, the body length of Dezhou donkey increases by 3 cm and the carcass weight increases by 6 kg. Seven SNVs (g.18093100G > T, g.18094587G > T, g.18106043G > T, g.18108764G > T, g.18110615T > G, g.18112000C > T and g.18114954T > G) of the NR6A1 gene were found to have a significant association with the TLN, body size and carcass traits of Dezhou donkey (p < 0.05), respectively. For instance, g.18114954C > T is significantly associated with lumber vertebrae number, the total number of thoracolumbar, and carcass weight, and individuals with TT genotype had significantly larger value than CC genotype (p < 0.05). Using these 7SNVs, 16 different haplotypes were estimated. Compared to Hap3Hap3, individuals homozygous for Hap2Hap2 showed significantly longer length in one thoracic spine (STL), the total thoracic vertebrae and one thoracolumbar spine. Our study will not only extend the understanding of genetic variation in the NR6A1 gene of Dezhou donkey, but also provide useful information for marker assisted selection in donkey breeding program.

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.

Evidence for and localization of proposed causative variants in cattle and pig genomes

BACKGROUND

This paper reviews the localization of published potential causative variants in contemporary pig and cattle reference genomes, and the evidence for their causality. In spite of the difficulties inherent to the identification of causative variants from genetic mapping and genome-wide association studies, researchers in animal genetics have proposed putative causative variants for several traits relevant to livestock breeding.

RESULTS

For this review, we read the literature that supports potential causative variants in 13 genes (ABCG2, DGAT1, GHR, IGF2, MC4R, MSTN, NR6A1, PHGK1, PRKAG3, PLRL, RYR1, SYNGR2 and VRTN) in cattle and pigs, and localized them in contemporary reference genomes. We review the evidence for their causality, by aiming to separate the evidence for the locus, the proposed causative gene and the proposed causative variant, and report the bioinformatic searches and tactics needed to localize the sequence variants in the cattle or pig genome.

CONCLUSIONS

Taken together, there is usually good evidence for the association at the locus level, some evidence for a specific causative gene at eight of the loci, and some experimental evidence for a specific causative variant at six of the loci. We recommend that researchers who report new potential causative variants use referenced coordinate systems, show local sequence context, and submit variants to repositories.

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.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Genome scanning reveals novel candidate genes for vertebral and teat number in the Beijing Black Pig

Porcine vertebral and teat numbers are variable and important economic traits in pig production. However, the quantitative trait loci (QTL) and candidate genes for both of these traits in the Beijing Black Pig are not yet known. In the present study, number of vertebrae and number of teats were obtained for 891 individuals of the Beijing Black Pig and genotyped using the Illumina Porcine 50 K BeadChip. Genome scanning was performed to detect associated variants and candidate genes for both traits using a genome-wide association study by tassel software. For vertebral number, 15 significant SNPs were located on SSC7. According to linkage disequilibrium analysis on SSC7, a haplotype block of 221 kb from 97.4 to 97.6 Mb was shown to contain a good candidate gene ABCD4. Interestingly, on SSC12, we recorded a novel QTL containing three significant SNPs and 34 annotated genes from 24.0 to 25.7 Mb for vertebral number. Of the 34 genes, nine Hox family genes (HOXB 1-7, 9, and 13) were found to be good candidate genes. Using the 34 genes, a gene ontology analysis was performed to detect enrichment of anterior/posterior pattern specification. For teat number, a novel chromosome-wide significant QTL was identified on SSC10. In this QTL region, one significant SNP was identified. The nearest gene, NTRK2, was regarded as a candidate gene. The present results expanded the QTL for vertebral and teat numbers and provided useful molecular markers for breeding in the Beijing Black Pig population.

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.

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 Analysis Identified BMPR1A as a Novel Candidate Gene Affecting the Number of Thoracic Vertebrae in a Large White × Minzhu Intercross Pig Population

Simple Summary

The number of thoracic vertebrae (NTV) and number of vertebrae (NV) varies among pig breeds with a high correlation of about 0.8. It is important to discover variants associated with the NTV by considering the effect of the NV in pig. The results suggest that regulation variants on SSC7 might play crucial roles in the NTV and the FOS on SSC7 should be further studied as a critical candidate gene. In addition, BMPR1A was identified as a novel candidate gene affecting the NTV in pigs.

Abstract

The number of vertebrae (NV), especially the number of thoracic vertebrae (NTV), varies among pig breeds. The NTV is controlled by vertebral segmentation and the number of somites during embryonic development. Although there is a high correlation between the NTV and NV, studies on a fixed NV have mainly considered the absolute numbers of thoracic vertebrae instead of vertebral segmentation. Therefore, this study aimed to discover variants associated with the NTV by considering the effect of the NV in pigs. The NTV and NV of 542 F2 individuals from a Large White × Minzhu pig crossbreed were recorded. All animals were genotyped for VRTN g.19034 A > C, LTBP2 c.4481A > C, and 37 missense or splice variants previously reported in a 951-kb interval on SSC7 and 147 single nucleotide polymorphisms (SNPs) on SSC14. To identify NTV-associated SNPs, we firstly performed a genome-wide association study (GWAS) using the Q + K (population structure + kinship matrix) model in TASSEL. With the NV as a covariate, the obtained data were used to identify the SNPs with the most significant genome-wide association with the NTV by performing a GWAS on a PorcineSNP60K Genotyping BeadChip. Finally, a conditional GWAS was performed by fixing this SNP. The GWAS showed that 31 SNPs on SSC7 have significant genome-wide associations with the NTV. No missense or splice variants were found to be associated with the NTV significantly. A linkage disequilibrium analysis suggested the existence of quantitative trait loci (QTL) in a 479-Kb region on SSC7, which contained a critical candidate gene FOS for the NTV in pigs. Subsequently, a conditional GWAS was performed by fixing M1GA0010658, the most significant of these SNPs. Two SNPs in BMPR1A were found to have significant genome-wide associations and a significant dominant effect. The leading SNP, S14_87859370, accounted for 3.86% of the phenotypic variance. Our study uncovered that regulation variants in FOS on SSC7 and in BMPR1A on SSC14 might play important roles in controlling the NTV, and thus these genetic factors may be harnessed for increasing the NTV in pigs.

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.