Quantitative trait loci (QTL) analysis in a Meishan x Göttingen cross population

Genome-wide association analysis reveals 6 copy number variations associated with the number of cervical vertebrae in Pekin ducks

As a critical developmental stage in vertebrates, the vertebral column formation process is under strict control; however, we observed variations in the number of cervical vertebrae in duck populations in our previous study. Here, we further explored the variations in the number of vertebrae in two duck populations: 421 Pekin duck × mallard F2 ducks and 850 Pekin ducks. Using resequencing data of 125 Pekin ducks with different numbers of cervical vertebrae and 352 Pekin duck × mallard F2 ducks with different numbers of thoracic vertebrae, we detected whole-genome copy number variations (CNVs) and implemented a genome-wide association study (GWAS) to identify the genetic variants related to the traits. The findings verified the existence of variations in the number of cervical vertebrae in duck populations. The number of cervical vertebrae in most ducks was 15, while that in a small number of the ducks was 14 or 16. The number of cervical vertebrae had a positive influence on the neck production, and one cervical vertebra addition could increase 11 g or 2 cm of duck neck. Genome-wide CNV association analysis identified six CNVs associated with the number of cervical vertebrae, and the associated CNV regions covered 15 genes which included WNT10A and WNT6. These findings improve our understanding of the variations in the number of vertebrae in ducks and lay a foundation for future duck breeding.

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.

Genomic evaluation and genome-wide association studies for total number of teats in a combined American and Danish Yorkshire pig populations selected in China

Joint genomic evaluation by combining data recordings and genomic information from different pig herds and populations is of interest for pig breeding companies because the efficiency of genomic selection (GS) could be further improved. In this work, an efficient strategy of joint genomic evaluation combining data from multiple pig populations is investigated. Total Teat Number (TTN), a trait that is equally recorded on 13 060 American Yorkshire (AY) populations (~14.68 teats) and 10 060 Danish Yorkshire (DY) pigs (~14.29 teats), was used to explore the feasibility and accuracy of GS combining datasets from different populations. We first estimated the genetic correlation (rg) of TTN between AY and DY pig populations (rg=0.79, se=0.23). Then we employed the genome-wide association study (GWAS) to identify QTL regions that are significantly associated with TTN and investigate the genetic architecture of TTN in different populations. Our results suggested that the genomic regions controlling TTN are slight different in the two Yorkshire populations, where the candidate QTL regions were on SSC 7 and SSC 8 for AY population and on SSC 7 for DY population. Finally, we explored an optimal way of genomic prediction for TTN via three different Genomic Best Linear Unbiased Prediction (GBLUP) models and we concluded that when TTN across populations are regarded as different, but correlated, traits in a multi-trait model, predictive abilities for both Yorkshire populations improve. As a conclusion, joint genomic evaluation for target traits in multiple pig populations is feasible in practice and more accurate, provided a proper model is used.

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

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

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

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

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.

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.

Differences in gene expression profiles for subcutaneous adipose, liver, and skeletal muscle tissues between Meishan and Landrace pigs with different backfat thicknesses

Backfat thickness is one of the most important traits of commercially raised pigs. Meishan pigs are renowned for having thicker backfat than Landrace pigs. To examine the genetic factors responsible for the differences, we first produced female crossbred pig lines by mating Landrace (L) × Large White (W) × Duroc (D) females (LWD) with Landrace (L) or Meishan (M) boars (i.e., LWD × L = LWDL for Landrace offspring and LWD × M = LWDM for the Meishan offspring). We confirmed that LWDM pigs indeed had a thicker backfat than LWDL pigs. Next, we performed gene expression microarray analysis in both genetic lines to examine differentially expressed genes (DEGs) in energy metabolism-related tissues, subcutaneous adipose (fat), liver, and longissimus dorsi muscle tissues. We analyzed the annotation of DEGs (2-fold cutoff) to functionally categorize them by Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathways. The number of DEGs in muscle tissues of both lines was much less than that in fat and liver tissues, indicating that DEGs in muscle tissues may not contribute much to differences in backfat thickness. In contrast, several genes related to muscle (in fat tissue) and lipid metabolism (in liver tissue) were more upregulated in LWDM pigs than LWDL pigs, indicating that those DEGs might be responsible for differences in backfat thickness. The different genome-wide gene expression profiles in the fat, liver, and muscle tissues between genetic lines can provide useful information for pig breeders.

Whole-genome resequencing reveals candidate mutations for pig prolificacy

Changes in pig fertility have occurred as a result of domestication, but are not understood at the level of genetic variation. To identify variations potentially responsible for prolificacy, we sequenced the genomes of the highly prolific Taihu pig breed and four control breeds. Genes involved in embryogenesis and morphogenesis were targeted in the Taihu pig, consistent with the morphological differences observed between the Taihu pig and others during pregnancy. Additionally, excessive functional non-coding mutations have been specifically fixed or nearly fixed in the Taihu pig. We focused attention on an oestrogen response element (ERE) within the first intron of the bone morphogenetic protein receptor type-1B gene (BMPR1B) that overlaps with a known quantitative trait locus (QTL) for pig fecundity. Using 242 pigs from 30 different breeds, we confirmed that the genotype of the ERE was nearly fixed in the Taihu pig. ERE function was assessed by luciferase assays, examination of histological sections, chromatin immunoprecipitation, quantitative polymerase chain reactions, and western blots. The results suggest that the ERE may control pig prolificacy via the cis-regulation of BMPR1B expression. This study provides new insight into changes in reproductive performance and highlights the role of non-coding mutations in generating phenotypic diversity between breeds.

Effects of increased vertebral number on carcass weight in PIC pigs

Variation of the vertebral number is associated with carcass traits in pigs. However, results from different populations do not match well with others, especially for carcass weight. Therefore, effects of increased vertebral number on carcass weight were investigated by analyzing the relationship between two loci multi-vertebra causal loci (NR6A1 g.748 C > T and VRTN g.20311_20312ins291) and carcass weight in PIC pigs. Results from the association study between vertebral number and carcass weight showed that increased thoracic number had negative effects on carcass weight, but the results were not statistically significant. Further, VRTN Ins/Ins genotype increased more than one thoracic than that of Wt/Wt genotype on average in this PIC population. Meanwhile, there was a significant negative effect of VRTN Ins on carcass weight (P < 0.05). Thus, our results suggested negative effect of increased thoracic number on carcass weight in PIC pigs.

Imputation-Based Whole-Genome Sequence Association Study Rediscovered the Missing QTL for Lumbar Number in Sutai Pigs

Resequencing a number of individuals of various breeds as reference population and imputing the whole-genome sequences of individuals that were genotyped with medium-density chips to perform an association study is a very efficient strategy. Previously, we performed a genome-wide association study (GWAS) of lumbar number using 60K SNPs from the porcine Illumina chips in 418 Sutai pigs and did not detect any significant signals. Therefore, we imputed the whole-genome sequences of 418 Sutai individuals from 403 deeply resequenced reference individuals and performed association tests. We identified a quantitative trait locus (QTL) for lumbar number in SSC1 with a P value of 9.01E-18 that was close to the potential causative gene of NR6A1. The result of conditioning on the top SNP association test indicated that only one QTL was responsible for this trait in SSC1. The linkage disequilibrium (LD) drop test result for the condition of the reported potential causative mutation (c.575T > C missense mutation of NR6A1) indicated that this mutation was probably not the underlying mutation that affected lumbar number in our study. As the first trial of imputed whole-genome sequence GWAS in swine, this approach can be also powerful to investigate complex traits in pig like in human and cattle.

Identification of loci affecting teat number by genome-wide association studies on three pig populations

OBJECTIVE

Three genome-wide association studies (GWAS) and a meta-analysis of GWAS were conducted to explore the genetic mechanisms underlying variation in pig teat number.

METHODS

We performed three GWAS and a meta-analysis for teat number on three pig populations, including a White Duroc×Erhualian F₂ resource population (n = 1,743), a Chinese Erhualian pig population (n = 320) and a Chinese Sutai pig population (n = 383).

RESULTS

We detected 24 single nucleotide polymorphisms (SNPs) that surpassed the genome-wide significant level on Sus Scrofa chromosomes (SSC) 1, 7, and 12 in the F₂ resource population, corresponding to four loci for pig teat number. We highlighted vertnin (VRTN) and lysine demethylase 6B (KDM6B) as two interesting candidate genes at the loci on SSC7 and SSC12. No significant associated SNPs were identified in the meta-analysis of GWAS.

CONCLUSION

The results verified the complex genetic architecture of pig teat number. The causative variants for teat number may be different in the three populations.

Possible introgression of the VRTN mutation increasing vertebral number, carcass length and teat number from Chinese pigs into European pigs

Vertnin (VRTN) variants have been associated with the number of thoracic vertebrae in European pigs, but the association has not been evidenced in Chinese indigenous pigs. In this study, we first performed a genome-wide association study in Chinese Erhualian pigs using one VRTN candidate causative mutation and the Illumina Porcine 60K SNP Beadchips. The VRTN mutation is significantly associated with thoracic vertebral number in this population. We further show that the VRTN mutation has pleiotropic and desirable effects on teat number and carcass (body) length across four diverse populations, including Erhualian, White Duroc × Erhualian F₂ population, Duroc and Landrace pigs. No association was observed between VRTN genotype and growth and fatness traits in these populations. Therefore, testing for the VRTN mutation in pig breeding schemes would not only increase the number of vertebrae and nipples, but also enlarge body size without undesirable effects on growth and fatness traits, consequently improving pork production. Further, by using whole-genome sequence data, we show that the VRTN mutation was possibly introgressed from Chinese pigs into European pigs. Our results provide another example showing that introgressed Chinese genes greatly contributed to the development and production of modern European pig breeds.

Genome-wide association analyses reveal significant loci and strong candidate genes for growth and fatness traits in two pig populations

BACKGROUND

Recently, genome-wide association studies (GWAS) have been reported on various pig traits. We performed a GWAS to analyze 22 traits related to growth and fatness on two pig populations: a White Duroc × Erhualian F2 intercross population and a Chinese Sutai half-sib population.

RESULTS

We identified 14 and 39 loci that displayed significant associations with growth and fatness traits at the genome-wide level and chromosome-wide level, respectively. The strongest association was between a 750 kb region on SSC7 (SSC for Sus scrofa) and backfat thickness at the first rib. This region had pleiotropic effects on both fatness and growth traits in F2 animals and contained a promising candidate gene HMGA1 (high mobility group AT-hook 1). Unexpectedly, population genetic analysis revealed that the allele at this locus that reduces fatness and increases growth is derived from Chinese indigenous pigs and segregates in multiple Chinese breeds. The second strongest association was between the region around 82.85 Mb on SSC4 and average backfat thickness. PLAG1 (pleiomorphic adenoma gene 1), a gene under strong selection in European domestic pigs, is proximal to the top SNP and stands out as a strong candidate gene. On SSC2, a locus that significantly affects fatness traits mapped to the region around the IGF2 (insulin-like growth factor 2) gene but its non-imprinting inheritance excluded IGF2 as a candidate gene. A significant locus was also detected within a recombination cold spot that spans more than 30 Mb on SSCX, which hampered the identification of plausible candidate genes. Notably, no genome-wide significant locus was shared by the two experimental populations; different loci were observed that had both constant and time-specific effects on growth traits at different stages, which illustrates the complex genetic architecture of these traits.

CONCLUSIONS

We confirm several previously reported QTL and provide a list of novel loci for porcine growth and fatness traits in two experimental populations with Chinese Taihu and Western pigs as common founders. We showed that distinct loci exist for these traits in the two populations and identified HMGA1 and PLAG1 as strong candidate genes on SSC7 and SSC4, respectively.

Feed intake, average daily gain, feed efficiency, and real-time ultrasound traits in Duroc pigs: II. Genomewide association

Efficient use of feed resources has become a clear challenge for the U.S. pork industry as feed costs continue to be the largest variable expense. The availability of the Illumina Porcine60K BeadChip has greatly facilitated whole-genome association studies to identify chromosomal regions harboring genes influencing those traits. The current study aimed at identifying genomic regions associated with variation in feed efficiency and several production traits in a Duroc terminal sire population, including ADFI, ADG, feed conversion ratio, residual feed intake (RFI), real-time ultrasound back fat thickness (BF), ultrasound muscle depth, intramuscular fat content (IMF), birth weight (BW at birth), and weaning weight (BW at weaning). Single trait association analyses were performed using Bayes B models with 35,140 SNP on 18 autosomes after quality control. Significance of nonoverlapping 1-Mb length windows (n = 2,380) were tested across 3 QTL inference methods: posterior distribution of windows variances from Monte Carlo Markov Chain, naive Bayes factor, and nonparametric bootstrapping. Genes within the informative QTL regions for the traits were annotated. A region ranging from166 to 140 Mb (4-Mb length) on SSC 1, approximately 8 Mb upstream of the MC4R gene, was significantly associated with ADFI, ADG, and BF, where SOCS6 and DOK6 are proposed as the most likely candidate genes. Another region affecting BW at weaning was identified on SSC 4 (84-85 Mb), harboring genes previously found to influence both human and cattle height: PLAG1, CHCHD7, RDHE2 (or SDR16C5), MOS, RPS20, LYN, and PENK. No QTL were identified for RFI, IMF, and BW at birth. In conclusion, we have identified several genomic regions associated with traits affecting nutrient utilization that could be considered for future genomic prediction to improve feed utilization.

Effect of VRTN gene polymorphisms on Duroc pig production and carcass traits, and their genetic relationships

The thoracic vertebral number is associated with body length and carcass traits, and represents one of the most important traits in the pig industry. Recent studies have shown that vertnin (VRTN) gene is associated with variations in the vertebral number in commercial European pigs. However, the genetic relationships and effect of this VRTN gene in pig production and carcass traits remain uncertain. Therefore, we investigated the genetic relationships among traits such as vertebral numbers, carcass weight and length-related traits, and meat production traits, and the effect of VRTN gene polymorphisms on these traits in a Duroc purebred population selected for its meat production traits. Highly positive genetic correlations were obtained between the thoracic vertebral numbers and length-related traits (0.56 to 0.84), whereas low correlations were obtained with production traits and carcass weight (-0.16 to 0.05). VRTN gene polymorphisms indicated that the number of thoracic vertebrae and length-related traits were significantly associated with the VRTN genotype, but had no significant effect on production traits and carcass weight. The results indicate that VRTN gene may be used as an effective selection marker to obtain pigs with high thoracic vertebral numbers and length-related traits, without adversely affecting meat production traits.

A genome-wide association study reveals dominance effects on number of teats in pigs

Dominance has been suggested as one of the genetic mechanisms explaining heterosis. However, using traditional quantitative genetic methods it is difficult to obtain accurate estimates of dominance effects. With the availability of dense SNP (Single Nucleotide Polymorphism) panels, we now have new opportunities for the detection and use of dominance at individual loci. Thus, the aim of this study was to detect additive and dominance effects on number of teats (NT), specifically to investigate the importance of dominance in a Landrace-based population of pigs. In total, 1,550 animals, genotyped for 32,911 SNPs, were used in single SNP analysis. SNPs with a significant genetic effect were tested for their mode of gene action being additive, dominant or a combination. In total, 21 SNPs were associated with NT, located in three regions with additive (SSC6, 7 and 12) and one region with dominant effects (SSC4). Estimates of additive effects ranged from 0.24 to 0.29 teats. The dominance effect of the QTL located on SSC4 was negative (−0.26 teats). The additive variance of the four QTLs together explained 7.37% of the total phenotypic variance. The dominance variance of the four QTLs together explained 1.82% of the total phenotypic variance, which corresponds to one-fourth of the variance explained by additive effects. The results suggest that dominance effects play a relevant role in the genetic architecture of NT. The QTL region on SSC7 contains the most promising candidate gene: VRTN. This gene has been suggested to be related to the number of vertebrae, a trait correlated with NT.

Genome-wide association analysis to identify SNP markers affecting teat numbers in an F2 intercross population between Landrace and Korean native pigs

Most reproductive traits have low heritability and are greatly affected by environmental factors. Teat number and litter size are traits related to the reproduction ability of pigs. To identify quantitative trait loci (QTLs) for teat number traits, a genome-wide association study (GWAS) was conducted using an F2 intercross between Landrace and Korean native pigs. Genotype analysis was performed using the porcine SNP 60 K beadchip. The GWAS was performed using a mixed-effects model and linear regression approach. When a genome-wide threshold was determined using the Bonferroni method (P = 1.61 × 10(-6)), 38 single nucleotide polymorphism (SNP) markers in pig chromosome 7 (SSC7) were significantly associated with three teat number traits (total teat number, left teat number, and right teat number). Among these, SNPs in 5 genes (HDDC3, LOC100156276, LOC100155863, ANPEP, SCAMP2) were selected for further study based primarily on their statistical significance. A significant association was detected in SCAMP2 g.25280 G>A for total teat number (P = 2.0 × 10(-12)), HDDC3 g.1319 G>A SNP for left teat number (P = 2.3 × 10(-7)), and SCAMP2 g.14198 G>A for right teat number (P = 4.7 × 10(-12)). These results provide valuable information about the selective breeding for desirable teat numbers in pigs.

The Possibility of TBC1D21 as a Candidate Gene for Teat Numbers in Pigs

Based on a quantitative traits locus (QTL) study using a F₂ intercross between Landrace and Korean native pigs, a significant QTL affecting teat numbers in SSC7 was identified. The strong positional candidate gene, TBC1D21, was selected due to its biological function for epithelial mesenchymal cell development. Sequence analysis revealed six single nucleotide polymorphisms (SNPs) in the TBC1D21 gene. Among these, two SNP markers, one silent mutation (SNP01) for g.13,050A>G and one missense mutation (SNP04) for c.829A>T (S277C), were genotyped and they showed significant associations with teat number traits (p value = 6.38E-05 for SNP01 and p value = 1.06E-07 for SNP04 with total teat numbers). Further functional validation of these SNPs could give valuable information for understanding the teat number variation in pigs.

A Whole Genome Association Study on Meat Quality Traits Using High Density SNP Chips in a Cross between Korean Native Pig and Landrace

A whole genome association (WGA) study was performed to detect significant polymorphisms for meat quality traits in an F₂ cross population (N = 478) that were generated with Korean native pig sires and Landrace dams in National Livestock Research Institute, Songwhan, Korea. The animals were genotyped using Illumina porcine 60k SNP beadchips, in which a set of 46,865 SNPs were available for the WGA analyses on ten carcass quality traits; live weight, crude protein, crude lipids, crude ash, water holding capacity, drip loss, shear force, CIE L, CIE a and CIE b. Phenotypes were regressed on additive and dominance effects for each SNP using a simple linear regression model, after adjusting for sex, sire and slaughter stage as fixed effects. With the significant SNPs for each trait (p<0.001), a stepwise regression procedure was applied to determine the best set of SNPs with the additive and/or dominance effects. A total of 106 SNPs, or quantitative trait loci (QTL) were detected, and about 32 to 66% of the total phenotypic variation was explained by the significant SNPs for each trait. The QTL were identified in most porcine chromosomes (SSCs), in which majority of the QTL were detected in SSCs 1, 2, 12, 13, 14 and 16. Several QTL clusters were identified on SSCs 12, 16 and 17, and a cluster of QTL influencing crude protein, crude lipid, drip loss, shear force, CIE a and CIE b were located between 20 and 29 Mb of SSC12. A pleiotropic QTL for drip loss, CIE L and CIE b was also detected on SSC16. These QTL need to be validated in commercial pig populations for genetic improvement in meat quality via marker-assisted selection.

High-resolution association mapping of number of teats in pigs reveals regions controlling vertebral development

BACKGROUND

Selection pressure on the number of teats has been applied to be able to provide enough teats for the increase in litter size in pigs. Although many QTL were reported, they cover large chromosomal regions and the functional mutations and their underlying biological mechanisms have not yet been identified. To gain a better insight in the genetic architecture of the trait number of teats, we performed a genome-wide association study by genotyping 936 Large White pigs using the Illumina PorcineSNP60 Beadchip. The analysis is based on deregressed breeding values to account for the dense family structure and a Bayesian approach for estimation of the SNP effects.

RESULTS

The genome-wide association study resulted in 212 significant SNPs. In total, 39 QTL regions were defined including 170 SNPs on 13 Sus scrofa chromosomes (SSC) of which 5 regions on SSC7, 9, 10, 12 and 14 were highly significant. All significantly associated regions together explain 9.5% of the genetic variance where a QTL on SSC7 explains the most genetic variance (2.5%). For the five highly significant QTL regions, a search for candidate genes was performed. The most convincing candidate genes were VRTN and Prox2 on SSC7, MPP7, ARMC4, and MKX on SSC10, and vertebrae δ-EF1 on SSC12. All three QTL contain candidate genes which are known to be associated with vertebral development. In the new QTL regions on SSC9 and SSC14, no obvious candidate genes were identified.

CONCLUSIONS

Five major QTL were found at high resolution on SSC7, 9, 10, 12, and 14 of which the QTL on SSC9 and SSC14 are the first ones to be reported on these chromosomes. The significant SNPs found in this study could be used in selection to increase number of teats in pigs, so that the increasing number of live-born piglets can be nursed by the sow. This study points to common genetic mechanisms regulating number of vertebrae and number of teats.

A further look at porcine chromosome 7 reveals VRTN variants associated with vertebral number in Chinese and Western pigs

The number of vertebrae is an economically important trait that affects carcass length and meat production in pigs. A major quantitative trait locus (QTL) for thoracic vertebral number has been repeatedly identified on pig chromosome (SSC) 7. To dissect the genetic basis of the major locus, we herein genotyped a large sample of animals from 3 experimental populations of Chinese and Western origins using 60K DNA chips. Genome-wide association studies consistently identified the locus across the 3 populations and mapped the locus to a 947-Kb region on SSC7. An identical-by-descent sharing assay refined the locus to a 100-Kb segment that harbors only two genes including VRTN and SYNDIG1L. Of them, VRNT has been proposed as a strong candidate of the major locus in Western modern breeds. Further, we resequenced the VRTN gene using DNA samples of 35 parental animals with known QTL genotypes by progeny testing. Concordance tests revealed 4 candidate causal variants as their genotypes showed the perfect segregation with QTL genotypes of the tested animals. An integrative analysis of evolutional constraints and functional elements supported two VRTN variants in a complete linkage disequilibrium phase as the most likely causal mutations. The promising variants significantly affect the number of thoracic vertebrae (one vertebra) in large scale outbred animals, and are segregating at rather high frequencies in Western pigs and at relatively low frequencies in a number of Chinese breeds. Altogether, we show that VRTN variants are significantly associated with the number of thoracic vertebrae in both Chinese and Western pigs. The finding advances our understanding of the genetic architecture of the vertebral number in pigs. Furthermore, our finding is of economical importance as it provides a robust breeding tool for the improvement of vertebral number and meat production in both Chinese indigenous pigs and Western present-day commercial pigs.

Detection of quantitative trait loci for teat number and female reproductive traits in Meishan × Large White F2 pigs

A quantitative trait locus (QTL) analysis of female reproductive data from a three-generation experimental cross between Meishan (MS) and Large White (LW) pig breeds is presented. Six F1 boars and 23 F1 sows, progeny of six LW boars and six MS sows, produced 573 F2 females and 530 F2 males. Six traits, i.e. teat number (TN), age at puberty (AP), ovulation rate (OR), weight at mating (WTM), number of viable embryos (NVE) and embryo survival (ES) at 30 days of gestation were analysed. Animals were genotyped for a total of 137 markers covering the entire porcine genome. Analyses were carried out based on interval mapping methods, using a line-cross (LC) regression and a half-full sib (HFS) maximum likelihood test. Genome-wide (GW) highly significant (P < 0.001) QTL were detected for WTM on SSC 7 and for AP on SSC 13. They explained, respectively, 14.5% and 8.9% of the trait phenotypic variance. Other GW significant (P < 0.05) QTL were detected for TN on SSC 3, 7, 8, 16 and 17, for OR on SSC 4 and 5, and for ES on SSC 9. Two additional chromosome-wide significant (P < 0.05) QTL were detected for TN, three for WTM, four for AP, three for OR, three for NVE and two for ES. With the exception of the two above-mentioned loci, the QTL explained from 1.2% to 4.6% of trait phenotypic variance. QTL alleles were in most cases not fixed in the grand-parental populations and Meishan alleles were not systematically associated with higher reproductive performance.

The Human Obesity Gene Map: The 2003 Update

This is the tenth update of the human obesity gene map, incorporating published results up to the end of October 2003 and continuing the previous format. Evidence from single-gene mutation obesity cases, Mendelian disorders exhibiting obesity as a clinical feature, quantitative trait loci (QTLs) from human genome-wide scans and animal crossbreeding experiments, and association and linkage studies with candidate genes and other markers is reviewed. Transgenic and knockout murine models relevant to obesity are also incorporated (N = 55). As of October 2003, 41 Mendelian syndromes relevant to human obesity have been mapped to a genomic region, and causal genes or strong candidates have been identified for most of these syndromes. QTLs reported from animal models currently number 183. There are 208 human QTLs for obesity phenotypes from genome-wide scans and candidate regions in targeted studies. A total of 35 genomic regions harbor QTLs replicated among two to five studies. Attempts to relate DNA sequence variation in specific genes to obesity phenotypes continue to grow, with 272 studies reporting positive associations with 90 candidate genes. Fifteen such candidate genes are supported by at least five positive studies. The obesity gene map shows putative loci on all chromosomes except Y. Overall, more than 430 genes, markers, and chromosomal regions have been associated or linked with human obesity phenotypes. The electronic version of the map with links to useful sites can be found at http://obesitygene.pbrc.edu.

The Human Obesity Gene Map: The 2004 Update

This paper presents the eleventh update of the human obesity gene map, which incorporates published results up to the end of October 2004. Evidence from single-gene mutation obesity cases, Mendelian disorders exhibiting obesity as a clinical feature, transgenic and knockout murine models relevant to obesity, quantitative trait loci (QTLs) from animal cross-breeding experiments, association studies with candidate genes, and linkages from genome scans is reviewed. As of October 2004, 173 human obesity cases due to single-gene mutations in 10 different genes have been reported, and 49 loci related to Mendelian syndromes relevant to human obesity have been mapped to a genomic region, and causal genes or strong candidates have been identified for most of these syndromes. There are 166 genes which, when mutated or expressed as transgenes in the mouse, result in phenotypes that affect body weight and adiposity. The number of QTLs reported from animal models currently reaches 221. The number of human obesity QTLs derived from genome scans continues to grow, and we have now 204 QTLs for obesity-related phenotypes from 50 genome-wide scans. A total of 38 genomic regions harbor QTLs replicated among two to four studies. The number of studies reporting associations between DNA sequence variation in specific genes and obesity phenotypes has also increased considerably with 358 findings of positive associations with 113 candidate genes. Among them, 18 genes are supported by at least five positive studies. The obesity gene map shows putative loci on all chromosomes except Y. Overall, >600 genes, markers, and chromosomal regions have been associated or linked with human obesity phenotypes. The electronic version of the map with links to useful publications and genomic and other relevant sites can be found at http://obesitygene.pbrc.edu.

Association of the porcine transforming growth factor beta type I receptor (TGFBR1) gene with growth and carcass traits

BACKGROUND

Growth and carcass traits are of great economic importance in livestock production. A large number of quantitative trait loci (QTL) have been identified for growth and carcass traits on porcine chromosome one (SSC1). A key positional candidate for this chromosomal region is TGFBR1 (transforming growth factor beta type I receptor). This gene plays a key role in inherited disorders at cardiovascular, craniofacial, neurocognitive, and skeletal development in mammals.

RESULTS

In this study, 27 polymorphic SNPs in the porcine TGFBR1 gene were identified on the University of Illinois Yorkshire × Meishan resource population. Three SNPs (SNP3, SNP43, SNP64) representing major polymorphic patterns of the 27 SNPs in F1 and F0 individuals of the Illinois population were selected for analyses of QTL association and genetic diversity. An association analysis for growth and carcass traits was completed using these three representative SNPs in the Illinois population with 298 F2 individuals and a large commercial population of 1008 animals. The results indicate that the TGFBR1 gene polymorphism (SNP64) is significantly associated (p < 0.05) with growth rates including average daily gains between birth and 56 kg (p = 0.049), between 5.5 and 56 kg (p = 0.024), between 35 and 56 kg (p = 0.021). Significant associations (p < 0.05) were also identified between TGFBR1 gene polymorphisms (SNP3/SNP43) and carcass traits including loin-eye-area (p = 0.022) in the Illinois population, and back-fat thickness (p = 0.0009), lean percentage (p = 0.0023) and muscle color (p = 0.021) in the commercial population. These three SNPs were also used to genotype a diverse panel of 130 animals representing 11 pig breeds. Alleles SNP3_T and SNP43_G were fixed in Pietrain and Sinclair pig breeds. SNP64_G allele was uniquely identified in Chinese Meishan pigs. Strong evidence of association (p < 0.01) between both SNP3 and SNP64 alleles and reproductive traits including gestation length and number of corpora lutea were also observed in the Illinois population.

CONCLUSION

This study gives the first evidence of association between the porcine TGFBR1 gene and traits of economic importance and provides support for using TGFBR1 markers for pig breeding and selection programs. The genetic diversities in different pig breeds would be helpful to understand the genetic background and migration of the porcine TGFBR1 gene.

Mapping and fine mapping of quantitative trait loci for the number of vertebrae in a White Duroc × Chinese Erhualian intercross resource population

The number of vertebrae is associated with body size and meat production in pigs. To identify quantitative trait loci (QTL) for the number of vertebrae, phenotypic values were measured in 1029 individuals from a White Duroc × Chinese Erhualian intercross F(2) population. A whole genome scan was performed with 194 microsatellite markers in the F(2) population. Four genome-wide significant QTL and eight chromosome-wide significant QTL for the number of vertebrae were identified on pig chromosomes (SSC) 1, 2, 6, 7, 10 and 12. The most significant QTL was detected on SSC7 with a confidence interval of 1 cM, explaining 42.32% of the phenotypic variance in the thoracic vertebral number. The significant QTL on SSC1, 2 and 7 confirmed previous reports. A panel of 276 animals representing seven Western and Chinese breeds was genotyped with 34 microsatellite markers in the SSC7 QTL region. No obvious selective sweep effect was observed in the tested breeds, indicating that intensive selection for enlarged body size in Western commercial breeds did not wipe out the genetic variability in the QTL region. The Q alleles for increased vertebral number originated from both Chinese Erhualian and White Duroc founder animals. A haplotype block of approximately 900 kb was found to be shared by all Q-bearing chromosomes of F(1) sires except for one distinct Q chromosome. The critical region harbours the newly reported VRTN gene associated with vertebral number. Further investigations are required to confirm whether VRTN or two other positional candidate genes, PROX2 and FOS, cause the QTL effect.

Changes in gene expression in a porcine preadipocyte cell line during differentiation

Adipocyte differentiation plays an important role in the formation of fat tissues in pigs and affects meat quality and productivity. Clarification of the nature of the pig genes that participate in adipocyte differentiation will provide a clue to the regulation of fat content and thickness in pig carcases by dietary control; it will also help to find target genes for exploring potentially useful polymorphisms for molecular breeding aimed at fat traits. We constructed a DNA oligomer microarray based on pig transcripts, and we used the array to investigate time-dependent changes in gene expression in the PSPA porcine preadipocyte cell line during differentiation into adipocytes. We selected genes with markedly altered expression (at least fivefold difference in comparison with expression in undifferentiated cells) and classified them into five groups according to gene expression pattern. In the early stage after stimulation of adipocyte differentiation, we observed up-regulation of many genes encoding proteins involved in regulating cell proliferation and transcription. Among the probes corresponding to transcripts that showed marked changes in expression, 27 were located within previously reported QTL regions for traits related to adipose tissues. These results will be valuable resources for finding the genes responsible for fat-related traits that have been identified in previous studies using various pig resource families.

Detection of quantitative trait loci for growth- and fatness-related traits in a large-scale White Duroc × Erhualian intercross pig population

Growth and fatness are economically important traits in pigs. In this study, a genome scan was performed to detect quantitative trait loci (QTL) for 14 growth and fatness traits related to body weight, backfat thickness and fat weight in a large-scale White Duroc × Erhualian F(2) intercross. A total of 76 genome-wide significant QTL were mapped to 16 chromosomes. The most significant QTL was found on pig chromosome (SSC) 7 for fatness with unexpectedly small confidence intervals of ∼2 cM, providing an excellent starting point to identify causal variants. Common QTL for both fatness and growth traits were found on SSC4, 5, 7 and 8, and shared QTL for fat deposition were detected on SSC1, 2 and X. Time-series analysis of QTL for body weight at six growth stages revealed the continuously significant effects of the QTL on SSC4 at the fattening period and the temporal-specific expression of the QTL on SSC7 at the foetus and fattening stages. For fatness traits, Chinese Erhualian alleles were associated with increased fat deposition except that at the major QTL on SSC7. For growth traits, most of White Duroc alleles enhanced growth rates except for those at three significant QTL on SSC6, 7 and 9. The results confirmed many previously reported QTL and also detected novel QTL, revealing the complexity of the genetic basis of growth and fatness in pigs.

Mapping QTL in the porcine MHC region affecting fatness and growth traits in a Meishan/Large White composite population

A number of studies have mapped QTL regulating porcine fatness and growth traits to the region of the major histocompatibility complex (MHC) on porcine chromosome 7 using various experimental crosses. The QTL results from crosses using the Chinese Meishan (MS) (slow growing and fat) are particularly interesting because the MS alleles have been found to be associated with increased growth rate and reduced backfat depth. We investigated these QTL further in a composite population derived previously over eight generations by intercrossing Meishan and the European Large White breeds. Genotype information from 32 markers in a 15cM target region was used in linkage and association analyses. A two-step variance component analysis identified QTL for three growth-related traits, explaining 19 ∼ 24% of the phenotypic variance with a confidence interval of 4 cM in the target region. SNP association analyses found that ss181128966 and ss181128924 within the QTL interval were strongly associated with the growth traits. Only weak signals for an effect on backfat depth were found in the association and linkage analyses, possibly because of past directional selection in the composite population.

Identification of a second gene associated with variation in vertebral number in domestic pigs

Background

The number of vertebrae in pigs varies and is associated with body size. Wild boars have 19 vertebrae, but European commercial breeds for pork production have 20 to 23 vertebrae. We previously identified two quantitative trait loci (QTLs) for number of vertebrae on Sus scrofa chromosomes (SSC) 1 and 7, and reported that an orphan nuclear receptor, NR6A1, was located at the QTL on SSC1. At the NR6A1 locus, wild boars and Asian local breed pigs had the wild-type allele and European commercial-breed pigs had an allele associated with increased numbers of vertebrae (number-increase allele).

Results

Here, we performed a map-based study to define the other QTL, on SSC7, for which we detected genetic diversity in European commercial breeds. Haplotype analysis with microsatellite markers revealed a 41-kb conserved region within all the number-increase alleles in the present study. We also developed single nucleotide polymorphisms (SNPs) in the 450-kb region around the QTL and used them for a linkage disequilibrium analysis and an association study in 199 independent animals. Three haplotype blocks were detected, and SNPs in the 41-kb region presented the highest associations with the number of vertebrae. This region encodes an uncharacterized hypothetical protein that is not a member of any other known gene family. Orthologs appear to exist not only in mammals but also birds and fish. This gene, which we have named vertnin (VRTN) is a candidate for the gene associated with variation in vertebral number. In pigs, the number-increase allele was expressed more abundantly than the wild-type allele in embryos. Among candidate polymorphisms, there is an insertion of a SINE element (PRE1) into the intron of the Q allele as well as the SNPs in the promoter region.

Conclusions

Genetic diversity of VRTN is the suspected cause of the heterogeneity of the number of vertebrae in commercial-breed pigs, so the polymorphism information should be directly useful for assessing the genetic ability of individual animals. The number-increase allele of swine VRTN was suggested to add an additional thoracic segment to the animal. Functional analysis of VRTN may provide novel findings in the areas of developmental biology.

Molecular characterization and association with carcass traits of the porcine SLC39A7 gene

In this study, the molecular characterization and potential association of SLC39A7 gene with carcass traits were investigated in pigs. The sequence of SLC39A7 cDNA was obtained by in silico cloning and RT-polymerase chain reaction (PCR). Two transcripts, variant 1 (2398 bp) and variant 2 (2088 bp), of the SLC39A7 gene were identified. Expression analysis of SLC39A7 in 10 different tissues by RT-PCR showed that variant 1 was ubiquitously expressed in all tissues analysed, but variant 2 was not found in fat tissue. The cDNA regions of variant 1 and 2 were organized in seven and eight exons respectively. A c.205G>A substitution in exon 3, which changes a codon for glycine into a codon for arginine, (p.Gly69Arg) and a c.1138-216T>C substitution in intron 6 were detected by PCR-HpaII-restriction fragment length polymorphisms (RFLP) and PCR-cofI-RFLP respectively. Significant differences were found in the allele frequencies of c.205G>A among six Chinese indigenous pig breeds and two commercial pig breeds. Linkage analysis showed that the c.205G>A polymorphism within the SLC39A7 gene was closely linked to the marker Sw1856 on pig chromosome 7 in a Large White x Meishan F(2) resource population. The QTL and association studies between polymorphisms of the SLC39A7 gene and carcass traits were carried out. Significant associations of the SLC39A7 polymorphisms with backfat thickness at thorax-waist (p < 0.05), average backfat thickness (p < 0.05) and leaf fat weight (p < 0.01) were found. Additional F-drop test or marker assisted association analyses also supported the association of the mutation in SLC39A7 with the above traits. Together, the present study provided the useful information for the characterization of SLC39A7 gene and potential association with carcass traits in pigs.

Using microarrays to identify positional candidate genes for QTL: the case study of ACTH response in pigs

BACKGROUND

Microarray studies can supplement QTL studies by suggesting potential candidate genes in the QTL regions, which by themselves are too large to provide a limited selection of candidate genes. Here we provide a case study where we explore ways to integrate QTL data and microarray data for the pig, which has only a partial genome sequence. We outline various procedures to localize differentially expressed genes on the pig genome and link this with information on published QTL. The starting point is a set of 237 differentially expressed cDNA clones in adrenal tissue from two pig breeds, before and after treatment with adrenocorticotropic hormone (ACTH).

RESULTS

Different approaches to localize the differentially expressed (DE) genes to the pig genome showed different levels of success and a clear lack of concordance for some genes between the various approaches. For a focused analysis on 12 genes, overlapping QTL from the public domain were presented. Also, differentially expressed genes underlying QTL for ACTH response were described. Using the latest version of the draft sequence, the differentially expressed genes were mapped to the pig genome. This enabled co-location of DE genes and previously studied QTL regions, but the draft genome sequence is still incomplete and will contain many errors. A further step to explore links between DE genes and QTL at the pathway level was largely unsuccessful due to the lack of annotation of the pig genome. This could be improved by further comparative mapping analyses but this would be time consuming.

CONCLUSION

This paper provides a case study for the integration of QTL data and microarray data for a species with limited genome sequence information and annotation. The results illustrate the challenges that must be addressed but also provide a roadmap for future work that is applicable to other non-model species.

Molecular genetics of the swine major histocompatibility complex, the SLA complex

The swine major histocompatibility complex (MHC) or swine leukocyte antigen (SLA) complex is one of the most gene-dense regions in the swine genome. It consists of three major gene clusters, the SLA class I, class III and class II regions, that span approximately 1.1, 0.7 and 0.5Mb, respectively, making the swine MHC the smallest among mammalian MHC so far examined and the only one known to span the centromere. This review summarizes recent updates to the Immuno Polymorphism Database-MHC (IPD-MHC) website (http://www.ebi.ac.uk/ipd/mhc/sla/) which serves as the repository for maintaining a list of all SLA recognized genes and their allelic sequences. It reviews the expression of SLA proteins on cell subsets and their role in antigen presentation and regulating immune responses. It concludes by discussing the role of SLA genes in swine models of transplantation, xenotransplantation, cancer and allergy and in swine production traits and responses to infectious disease and vaccines.

Genome-wide QTL mapping for three traits related to teat number in a White Duroc x Erhualian pig resource population

BACKGROUND

Teat number is an important fertility trait for pig production, reflecting the mothering ability of sows. It is also a discrete and often canalized trait presenting bilateral symmetry with minor differences between the two sides, providing a potential power to evaluate fluctuating asymmetry and developmental instability. The knowledge of its genetic control is still limited. In this study, a genome-wide scan was performed with 183 microsatellites covering the pig genome to identify quantitative trait loci (QTL) for three traits related to teat number including the total teat number (TTN), the teat number at the left (LTN) and right (RTN) sides in a large scale White Duroc x Erhualian resource population.

RESULTS

A sex-average linkage map with a total length of 2350.3 cM and an average marker interval of 12.84 cM was constructed. Eleven genome-wide significant QTL for TTN were detected on 8 autosomes including pig chromosomes (SSC) 1, 3, 4, 5, 6, 7, 8 and 12. Six suggestive QTL for this trait were detected on SSC6, 9, 13, 14 and 16. Eight chromosomal regions each on SSC1, 3, 4, 5, 6, 7, 8 and 12 showed significant associations with LTN. These regions were also evidenced as significant QTL for RTN except for those on SSC6 and SSC8. The most significant QTL for the 3 traits were all located on SSC7. Erhualian alleles at most of the identified QTL had positive additive effects except for three QTL on SSC1 and SSC7, at which White Duroc alleles increased teat numbers. On SSC1, 6, 9, 13 and 16, significant dominance effects were observed on TTN, and predominant imprinting effect on TTN was only detected on SSC12.

CONCLUSION

The results not only confirmed the QTL regions from previous experiments, but also identified five new QTL for the total teat number in swine. Minor differences between the QTL regions responsible for LTN and RTN were validated. Further fine mapping should be focused on consistently identified regions with small confidence intervals, such as those on SSC1, SSC7 and SSC12.

Genetic evidence for the introgression of Western NR6A1 haplotype into Chinese Licha breed associated with increased vertebral number

There is evidence that NR6A1 is a strong candidate for being a causal gene underlying vertebral number in pigs. The Licha Black is one of the leanest Chinese indigenous pig breeds, having an average vertebral number of 21.5. The introgression of Western germplasm into Licha Black, resulting in increased vertebral number, has been assumed but is not confirmed. This study detected allele frequencies of the NR6A1 causative mutation (c.575T>C) in 519 pigs from three Western and seven Chinese breeds including Licha Black, and evaluated the genetic variation in a 650-kb region containing NR6A1 in the 10 breeds. Allele T for increased vertebral number was fixed in Western breeds. In contrast, this allele was very rare in most of the Chinese native breeds. Notably, the T allele was present in the Licha Black at a rather higher frequency (0.585) and in the Laiwu at lower frequency (0.250). As expected, selection pressure has wiped out the genetic variability in the 650 kb region in Western breeds. Conversely, Chinese indigenous breeds showed a high degree of genetic variability in this region. However, the Licha Black displayed dramatically reduced heterozygosity at the loci proximal to the causative mutation. Moreover, a high proportion (45.9%) of Licha Black pigs and a small number (21%) of Laiwu pigs had the Western NR6A1 haplotype, and the two breeds showed closer relationships with Western commercial breeds than other Chinese breeds in the phylogenic tree. When the results are taken together, this study supports the assumption that the Western NR6A1 haplotypes were introduced into Licha Black and possibly Laiwu and are associated with increased vertebral number.

Structure and polymorphism analysis of transforming growth factor beta receptor 1 (TGFBR1) in pigs

Many quantitative trait loci (QTL) for growth and reproductive traits have been detected on the porcine chromosome region 1qter (SSC1qter), making it one of the most important genomic regions for pig breeding. SSC1q corresponds to human chromosome 9, on which lies transforming growth factor beta receptor 1 (TGFBR1). We cloned the porcine TGFBR1 cDNA and gene (as a candidate for QTL) and analyzed the gene structure and polymorphism. Porcine TGFBR1 consists of 9 exons and 8 introns. Intron 2 is alternatively spliced at the acceptor site, resulting in two kinds of mRNA, with putative open reading frames of 1500 and 1512 bp in length. The shorter one encodes 499 amino acid residues. The amino acid sequence has 96.2 and 97.2% sequence similarity to those of human and bovine TGFBR1, respectively. The sequence similarity between porcine and murine TGFBR1 is 95.6%. We detected three single-nucleotide substitutions in exons 1, 2, and 7. Those in exons 1 and 7 are nonsynonymous substitutions resulting in Pro8Ser and Ile413Val substitutions, respectively.