Genome-wide study refines the quantitative trait locus for number of ribs in a Large White × Minzhu intercross pig population and reveals a new candidate gene

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.

Multi-organ phenotypes in mice lacking latent TGFβ binding protein 2 (LTBP2)

BACKGROUND

Latent TGFβ binding protein-2 (LTBP2) is a fibrillin 1 binding component of the microfibril. LTBP2 is the only LTBP protein that does not bind any isoforms of TGFβ, although it may interfere with the function of other LTBPs or interact with other signaling pathways.

RESULTS

Here, we investigate mice lacking Ltbp2 (Ltbp2-/- ) and identify multiple phenotypes that impact bodyweight and fat mass, and affect bone and skin development. The alterations in skin and bone development are particularly noteworthy since the strength of these tissues is differentially affected by loss of Ltbp2. Interestingly, some tissues that express high levels of Ltbp2, such as the aorta and lung, do not have a developmental or homeostatic phenotype.

CONCLUSIONS

Analysis of these mice show that LTBP2 has complex effects on development through direct effects on the extracellular matrix (ECM) or on signaling pathways that are known to regulate the ECM.

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.

Identification of Candidate Genes for Economically Important Carcass Cutting in Commercial Pigs through GWAS

During the process of pork production, the carcasses of pigs are divided and sold, which provides better economic benefits and market competitiveness for pork production than selling the carcass as a whole. Due to the significant cost of post-slaughter phenotypic measurement, the genetic architecture of tenderloin weight (TLNW) and rib weight (RIBW)-important components of pig carcass economic value-remain unknown. In this study, we conducted genome-wide association studies (GWAS) for TLNW and RIBW traits in a population of 431 Duroc × Landrace × Yorkshire (DLY) pigs. In our study, the most significant single nucleotide polymorphism (SNP) associated with TLNW was identified as ASGA0085853 (3.28 Mb) on Sus scrofa chromosome 12 (SSC12), while for RIBW, it was Affx-1115046258 (172.45 Mb) on SSC13. Through haplotype block analysis, we discovered a novel quantitative trait locus (QTL) associated with TLNW, spanning a 5 kb region on SSC12, and a novel RIBW-associated QTL spanning 1.42 Mb on SSC13. Furthermore, we hypothesized that three candidate genes, TIMP2 and EML1, and SMN1, are associated with TLNW and RIBW, respectively. Our research not only addresses the knowledge gap regarding TLNW, but also serves as a valuable reference for studying RIBW. The identified SNP loci strongly associated with TLNW and RIBW may prove useful for marker-assisted selection in pig breeding programs.

Identification of LTBP2 gene polymorphisms and their association with thoracolumbar vertebrae number, body size, and carcass traits in Dezhou donkeys

The number of thoracolumbar vertebrae in Dezhou donkeys varies from 22 to 24 and is associated with body size and carcass traits. In mammals, the latent transforming growth factor beta binding protein 2 (LTBP2) has been found to have some functions in the development of thoracolumbar vertebrae. The relationship between LTBP2 and TLN (the number of thoracolumbar vertebrae) of Dezhou donkeys is yet to be reported. The purposes of this study are as follows: 1) to quantify the effect of thoracolumbar vertebrae number variation of Dezhou donkeys on body size and carcass trait; 2) to study the distribution of single nucleotide variants (SNVs) in the LTBP2 gene of Dezhou donkeys; and 3) to explore whether these SNVs can be used as candidate sites to study the mechanism of Dezhou donkey muti-thoracolumbar vertebrae development. The TLN, body size, and carcass traits of 392 individuals from a Dezhou donkey breed were recorded. All animals were sequenced for LTBP2 using GBTS liquid chip and 16 SNVs were used for further analysis. We then analyzed the relationship between these SNVs with TLN, body size, and carcass traits. The results showed that: 1) c.5547 + 860 C > T, c.5251 + 281 A > C, c.3769 + 40 C > T, and c.2782 + 3975 A > G were complete genetic linkages and significantly associated with thoracic vertebrae number (TN) (p < 0.05) (wild-type homozygotes had more TN than heterozygotes); 2) c.1381 + 768 T > G and c.1381 + 763 G > T were significantly associated with lumber vertebrae number (LN) (p < 0.05); 3) c.1003 + 704 C > T, c.1003 + 651 C > T, c.1003 + 626 A > G, and c.812 + 22526 T > G were significantly associated with chest circumference (CHC), front carcass weight (CWF), after carcass weight (CWA), and carcass weight (CW) (p < 0.05) (wild-type homozygotes were larger than other genotypes in CHC, CWF, CWA, and CW); and 4) the effect of variation is not consistent in c.565 + 11921 A > G, c.565 + 6840 A > G, c.565 + 3453 C > T, and c.494 + 5808 C > T. These results provide useful information that the polymorphism of LTBP2 is significantly associated with TLN, body size, and carcass traits in Dezhou donkeys, which can serve as a molecule marker to improve donkey production performance.

Genome-wide association study revealed ABCD4 on SSC7 and GREB1L and MIB1 on SSC6 as crucial candidate genes for rib number in Beijing Black pigs

As one of the few animals with variation in the number of rib pairs (RIB), the pig is a good model to study the mechanism of RIB regulation. Quantitative trait loci (QTL) for porcine RIB are present on Sus scrofa chromosome 7 (SSC7). Although several candidate genes exist in this QTL region on SSC7, the causal gene has yet to be verified. Beijing Black pig with 14-17 RIB is a good population for candidate gene mining and 1104 individuals were genotyped using the Illumina Porcine 50K BeadChip. A total of 14 SNPs from 95.49 to 97.78 Mb on SSC7 showed genome-wide significant association with RIB. On SSC7, a locuszoom plot using pairwise linkage disequilibrium displayed the narrowest linkage region encompassing only two genes, ABCD4 and VRTN. In mice, a transcriptome expression profile was obtained using three embryos at E9.5 (the critical period for rib formation). ABCD4 was highly expressed, but no expression of VRTN was detected. On SSC6, there were four genome-wide significant SNPs from 106.42 to 106.92 Mb associated with RIB. GREB1L and MIB1, in this region, were regarded as novel candidate genes. These results revealed a crucial candidate causal gene on SSC7 and novel genes on SSC6 for rib number and provided interesting new insights into its genetic basis.

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.

Genome-Wide Detection of Runs of Homozygosity in Laiwu Pigs Revealed by Sequencing Data

Laiwu pigs, distinguished by their high intramuscular fat of 7-9%, is an indigenous pig breed of China, and recent studies also found that Laiwu pigs showed high resistance to Porcine circovirus type 2. However, with the introduction of commercial varieties, the population of Laiwu pigs has declined, and some lineages have disappeared, which could result in inbreeding. Runs of homozygosity (ROH) can be used as a good measure of individual inbreeding status and is also normally used to detect selection signatures so as to map the candidate genes associated with economically important traits. In this study, we used data from Genotyping by Genome Reducing and Sequencing to investigate the number, length, coverage, and distribution patterns of ROH in 93 Chinese Laiwu pigs and identified genomic regions with a high ROH frequency. The average inbreeding coefficient calculated by pedigree was 0.021, whereas that estimated by all detected ROH segments was 0.133. Covering 13.4% of the whole genome, a total of 7,508 ROH segments longer than 1 Mb were detected, whose average length was 3.76 Mb, and short segments (1-5 Mb) dominated. For individuals, the coverage was in the range between 0.56 and 36.86%. For chromosomes, SSC6 had the largest number (n = 688), and the number of ROH in SSC12 was the lowest (n = 215). Thirteen ROH islands were detected in our study, and 86 genes were found within those regions. Some of these genes were correlated with economically important traits, such as meat quality (ECI1, LRP12, NDUFA4L2, GIL1, and LYZ), immunity capacity (IL23A, STAT2, STAT6, TBK1, IFNG, and ITH2), production (DCSTAMP, RDH16, and GDF11), and reproduction (ODF1 and CDK2). A total of six significant Gene Ontology terms and nine significant Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were identified, most of which were correlated with disease resistance and biosynthesis processes, and one KEGG pathway was related to lipid metabolism. In addition, we aligned all of the ROH islands to the pig quantitative trait loci (QTL) database and finally found eight QTL related to the intramuscular fat trait. These results may help us understand the characteristics of Laiwu pigs and provide insight for future breeding strategies.

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.

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.

A genome-wide association study of reproduction traits in four pig populations with different genetic backgrounds

Objective

Genome-wide association study and two meta-analysis based on GWAS performed to explore the genetic mechanism underlying variation in pig number born alive (NBA) and total number born (TNB).

Methods

Single trait GWAS and two meta-analysis (single-trait meta analysis and multi-trait meta analysis) were used in our study for NBA and TNB on 3,121 Yorkshires from 4 populations, including three different American Yorkshire populations (n = 2,247) and one British Yorkshire populations (n = 874).

Results

The result of single trait GWAS showed that no significant associated single nucleotide polymorphisms (SNPs) were identified. Using single-trait meta analysis and multi-trait meta analysis within populations, 11 significant loci were identified associated with target traits. Spindlin 1, vascular endothelial growth factor A, forkhead box Q1, msh homeobox 1, and LHFPL tetraspan submily member 3 are five functionally plausible candidate genes for NBA and TNB. Compared to the single population GWAS, single-trait Meta analysis can improve the detection power to identify SNPs by integrating information of multiple populations. The multiple-trait analysis reduced the power to detect trait-specific loci but enhanced the power to identify the common loci across traits.

Conclusion

In total, our findings identified novel genes to be validated as candidates for NBA and TNB in pigs. Also, it enabled us to enlarge population size by including multiple populations with different genetic backgrounds and increase the power of GWAS by using meta analysis.

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.

Molecular cloning of WIF1 and HMGA2 reveals ear-preferential expression while uncovering a missense mutation associated with porcine ear size in WIF1

Considerable diversity exists in porcine ear size, which is an important morphological feature of pig breeds. Previously, we localized four crucial candidate genes-high mobility group AT-hook 2 (HMGA2), LEM domain-containing 3 (LEMD3), methionine sulfoxide reductase B3 (MSRB3) and Wnt inhibitory factor 1 (WIF1)-on Sus Scrofa chromosome 5 affecting porcine ear size, then cloned LEMD3 and MSBR3. In this study, we performed rapid amplification of cDNA ends to obtain full-length cDNA sequences of 2338-bp WIF1 and 2998-bp HMGA2. Using quantitative real-time PCR, we revealed that WIF1 expression was highest in ear cartilage of 60-day-old pigs and that this is therefore a better candidate gene for ear size than HMGA2. We further screened coding sequence variants in both genes and identified only one missense mutation (WIF1:c.1167C>G) in a conserved epidermal growth factor-like domain from the mammalian WIF1 protein. The protein-altering mutation was significantly associated with ear size across the Large White × Minzhu hybrid and Beijing Black pig populations. When WIF1:c.1167C>G was included as fixed effect in the model to re-run a genome-wide association study in the Large White × Minzhu intercross population the P-value of the peak SNP on SSC5 from re-running the genome-wide association study dropped from 2.45E-12 to 7.33E-05. Taken together, the WIF1:c.1167C>G could be an important mutation associated with ear size. Our findings provide helpful information for further studies of the molecular mechanisms controlling porcine ear size.

Genome-wide association studies and meta-analysis uncovers new candidate genes for growth and carcass traits in pigs

Genome-wide association studies (GWAS) have been widely used in the genetic dissection of complex traits. As more genomic data is being generated within different commercial or resource pig populations, the challenge which arises is how to collectively investigate the data with the purpose to increase sample size and implicitly the statistical power. This study performs an individual population GWAS, a joint population GWAS and a meta-analysis in three pig F2 populations. D1 is derived from European type breeds (Piétrain, Large White and Landrace), D2 is obtained from an Asian breed (Meishan) and Piétrain, and D3 stems from a European Wild Boar and Piétrain, which is the common founder breed. The traits investigated are average daily gain, backfat thickness, meat to fat ratio and carcass length. The joint and the meta-analysis did not identify additional genomic clusters besides the ones discovered via the individual population GWAS. However, the benefit was an increased mapping resolution which pinpointed to narrower clusters harboring causative variants. The joint analysis identified a higher number of clusters as compared to the meta-analysis; nevertheless, the significance levels and the number of significant variants in the meta-analysis were generally higher. Both types of analysis had similar outputs suggesting that the two strategies can complement each other and that the meta-analysis approach can be a valuable tool whenever access to raw datasets is limited. Overall, a total of 20 genomic clusters that predominantly overlapped at various extents, were identified on chromosomes 2, 7 and 17, many confirming previously identified quantitative trait loci. Several new candidate genes are being proposed and, among them, a strong candidate gene to be taken into account for subsequent analysis is BMP2 (bone morphogenetic protein 2).

A genome scan for selection signatures in Taihu pig breeds using next-generation sequencing

Taihu pig breeds are the most prolific breeds of swine in the world, and they also have superior economic traits, including high resistance to disease, superior meat quality, high resistance to crude feed and a docile temperament. The formation of these phenotypic characteristics is largely a result of long-term artificial or natural selection. Therefore, exploring selection signatures in the genomes of the Taihu pigs will help us to identify porcine genes related to productivity traits, disease and behaviour. In this study, we used both intra-population (Relative Extend Haplotype Homozygosity Test (REHH)) and inter-population (the Cross-Population Extend Haplotype Homozygosity Test (XPEHH); F-STATISTICS, F ST ) methods to detect genomic regions that might be under selection process in Taihu pig breeds. As a result, we found 282 (REHH) and 112 (XPEHH) selection signature candidate regions corresponding to 159.78 Mb (6.15%) and 62.29 Mb (2.40%) genomic regions, respectively. Further investigations of the selection candidate regions revealed that many genes under these genomic regions were related to reproductive traits (such as the TLR9 gene), coat colour (such as the KIT gene) and fat metabolism (such as the CPT1A and MAML3 genes). Furthermore, gene enrichment analyses showed that genes under the selection candidate regions were significantly over-represented in pathways related to diseases, such as autoimmune thyroid and asthma diseases. In conclusion, several candidate genes potentially under positive selection were involved in characteristics of Taihu pig. These results will further allow us to better understand the mechanisms of selection in pig breeding.

Rapid Communication: High-resolution quantitative trait loci analysis identifies LTBP2 encoding latent transforming growth factor beta binding protein 2 associated with thoracic vertebrae number in a large F2 intercross between Landrace and Korean native pigs

Number of vertebrae is associated with body size and meat productivity in pigs. The aim of this study was to identify QTL and associated positional candidate genes affecting the number of thoracic vertebrae (THO). A genomewide association study was conducted in a large resource population derived from an F intercross between Landrace and Korean native pigs using the Porcine SNP 60K BeadChip and the genomewide complex trait analysis (GCTA) program based on a linear mixed-effects model. A total of 38,385 SNP markers from 1,105 F progeny were analyzed for the THO trait after filtering for quality control. A total of 90 genomewide significant SNP markers ( < 1.30 × 10) on SSC 7 covering a 20-Mb region were identified for THO in this study. Several previous studies also mapped QTL for vertebral numbers in this region. The strongest association signals were detected at ASGA0035500 (-value = 4.46 × 10; 103,574,383 bp) and DIAS0000795 (-value = 4.46 × 10; 103,594,753 bp). The QTL region on SSC 7 for THO encompasses and , which are previously described candidate genes for vertebral number variation. To refine the QTL region, a haplotype-based linkage and linkage disequilibrium (LALD) analysis using the DualPHASE program was applied because subsequent conditional association and haplotype block analyses could not resolve the region that contains the 2 loci. The LALD analysis refined the critical region to a 533.9-kb region including ; was located outside the critical region. The gene encoding latent transforming growth factor beta binding protein 2 is involved in bone metabolisms. Based on these data, we propose as a positional candidate gene for THO in pigs. After further functional studies and verification of the association in other independent populations, these results could be useful for optimizing breeding programs that improve THO and other economically important traits in pigs.

Genome-wide genetic structure and differentially selected regions among Landrace, Erhualian, and Meishan pigs using specific-locus amplified fragment sequencing

As typical Chinese indigenous pig breeds, Erhualian and Meishan have been widely used to produce new strain or breed in the world. However, the genetic basis of characteristics of these pig breeds is still limited. Moreover, considering cost and output of sequencing, it is necessary to further develop cost-effective method for pig genome screening. To contribute on this issue, we developed a SLAF-seq (specific-locus amplified fragment sequencing) method for pigs and applied it to analyze the genetic difference among Landrace, Erhualian, and Meishan pigs. A total of 453.75 million reads were produced by SLAF-seq. After quality-control, 165,670 SNPs (single nucleotide polymorphisms) were used in further analysis. The results showed that Landrace had distinct genetic relationship compared to Erhualian (F_ST = 0.5480) and Meishan (F_ST = 0.5800), respectively, while Erhualian and Meishan held the relatively close genetic relationship (F_ST = 0.2335). Furthermore, a genome-wide scanning revealed 268 differentially selected regions (DSRs) with 855 genes and 256 DSRs with 347 genes between Landrace and the two Chinese indigenous pig breeds and between Erhualian and Meishan, respectively. This study provides a new cost-effective method for pig genome study and might contribute to a better understanding on the formation mechanism of genetic difference among pigs with different geographical origins.