- Research article
- Open Access
Analysis of porcine MUC4 gene as a candidate gene for prolificacy QTL on SSC13 in an Iberian × Meishan F2 population
© Balcells et al; licensee BioMed Central Ltd. 2011
- Received: 4 July 2011
- Accepted: 31 October 2011
- Published: 31 October 2011
Reproductive traits, such as prolificacy, are of great interest to the pig industry. Better understanding of their genetic architecture should help to increase the efficiency of pig productivity through the implementation of marker assisted selection (MAS) programmes.
The Mucin 4 (MUC4) gene has been evaluated as a candidate gene for a prolificacy QTL described in an Iberian × Meishan (Ib × Me) F2 intercross. For association analyses, two previously described SNPs (DQ124298:g.243A>G and DQ124298:g.344A>G) were genotyped in 347 pigs from the Ib × Me population. QTL for the number of piglets born alive (NBA) and for the total number of piglets born (TNB) were confirmed on SSC13 at positions 44 cM and 51 cM, respectively. The MUC4 gene was successfully located within the confidence intervals of both QTL. Only DQ124298:g.344A>G MUC4 polymorphism was significantly associated with both NBA and TNB (P-value < 0.05) with favourable effects coming from the Meishan origin. MUC4 expression level was determined in F2 sows displaying extreme phenotypes for the number of embryos (NE) at 30-32 days of gestation. Differences in the uterine expression of MUC4 were found between high (NE ≥ 13) and low (NE ≤ 11) prolificacy sows. Overall, MUC4 expression in high prolificacy sows was almost two-fold increased compared with low prolificacy sows.
Our data suggest that MUC4 could play an important role in the establishment of an optimal uterine environment that would increase embryonic survival during pig gestation.
- Quantitative Trait Locus
- Quantitative Trait Locus Effect
- Significant Quantitative Trait Locus
- MUC4 Expression
- MUC4 Gene
Prolificacy traits have been widely explored during the last decades as a potential tool for increasing efficiency of sow productivity in the pig industry. Genetic improvement programmes have achieved moderate gains in prolificacy related traits owing to their low heritability, late expression in life and sex limitation . Increasing knowledge on the genetic architecture of prolificacy traits would provide new tools to improve the efficiency of genetic selection by implementing marker assisted selection (MAS).
So far a relatively low number of quantitative trait loci (QTL) for prolificacy traits reaching the genome-wide significance level have been identified [2, 3]. The most significant QTL affecting the number of piglets born alive (NBA) and the total number of piglets born (TNB) were described by Noguera et al.  in the same resource population as used in the present study, an Iberian (Ib) by Meishan (Me) F2 intercross. A number of epistatic QTL were also detected, thus indicating that the genetic architecture of reproductive traits is built as a complex network of interactions throughout the genome. Some of these epistatic QTL were further confirmed and putative candidate interacting genes were identified .
Porcine chromosome 13 (SSC13) harbours the most significant QTL for TNB and NBA . The Mucin 4 (MUC4) gene is located within the confidence interval of prolificacy QTL. Mucins are large heterodimeric glycoproteins commonly located on apical surfaces of many wet-surfaced epithelia that play a key role in the lubrication and protection of the uterine mucosa [5–7]. They have been shown to present anti-adhesive and anti-recognition properties which are necessary to protect the endometrium from the binding and invasion of the trophoectoderm [8, 9]. A role of MUC4 has been pointed out in rodents and pigs during pregnancy although its expression during the peri-implantational period varies depending on the type of implantation in each species. In mice and rats, which have an invasive type of implantation, MUC4 expression is downregulated to generate the receptive state for uterine implantation [8–12]. Conversely, in pigs, where a non-invasive epitheliochorial placental attachment takes place, MUC4 is upregulated in the uterus . A protective role has been suggested for MUC4 owing to the fact that it is localized on the endometrium epithelium blocking the access of different substrates to the cell surface . The endometrium is then protected from proteolytic activity of porcine conceptus  and from microbial invasion  resulting in better uterine conditions for embryo development. In pigs, the disruption of the uterine microenvironment could affect embryo viability which could lead to prenatal mortality rates ranging from 20 to 46% . The improvement of the uterine microenvironment would increase embryonic survival and, in consequence, the number of piglets born alive.
In humans, polymorphisms in the MUC4 nucleotide sequence have been significantly associated with the development of endometriosis and endometriosis related infertility . However, no association with implantation failure has been detected . In livestock species, the genetic association of MUC4 gene variants with reproductive traits has not yet been explored. In pigs, polymorphisms in the MUC4 gene were shown to be in linkage disequilibrium with susceptibility/resistance to Enterotoxigenic Escherichia coli (ETEC) F4ab/ac infection .
In the current study, we have examined the porcine MUC4 gene as a functional and positional candidate gene to explain the prolificacy QTL previously identified on SSC13 in the Ib × Me population .
Refinement of SSC13 QTL for NBA and TNB
The resulting linkage map for SSC13 in the Ib × Me population was as follows (distance is indicated in centimorgan (cM)): S0076 (0.0) - ITIH3 (18.7) - SWR1008 (27.5) - MUC4 (39.1) - SW398 (55.8) - SW2440 (78.5) - SW769 (91.5). The position of the MUC4 gene is in agreement with the pig genome sequence (Sus scrofa, Ensembl release 64 - September 2011).
Results of QTL analyses, association tests and marker assisted association tests for prolificacy traits
QTL model (model 1)
MUC4 association (model 2)
QTL + MUC4 association model (model 3)1
a MUC4 (SE)4
P -QTL + MUC4
Candidate gene association analyses
Allelic frequencies of the DQ124298:g.243A>G and DQ124298:g.344A>G MUC4 polymorphisms in the Me × Ib F2 population
The results obtained with association studies must be interpreted with caution due to the fact that the extensive linkage disequilibrium described in F2 crosses can bias the estimations. For this reason, in order to detect possible false positives obtained in association studies with DQ124298:g.344A>G SNP, a marker assisted association test (MAAT) including the QTL effect in the association test  was performed with model 3. Results from the MAAT are summarized in Table 1. The DQ124298:g.344A>G SNP effect on NBA and TNB remained significant when the QTL effect was considered which is in agreement with the association studies analyses. Nevertheless, the significance of the QTL disappeared when the DQ124298:g.344A>G SNP was included in the model. MAAT shows that the DQ124298:g.344A>G SNP genotype explains better the effects on NBA and TNB and confirms the association detected with model 2.
Expression analysis of the porcine MUC4 gene
Statistical association between the DQ124298:g.344A>G mutation in the porcine MUC4 gene and prolificacy related traits has been reported. No association was found, however, between the DQ124298:g.243A>G SNP and the prolificacy related traits although this SNP is located 100bp upstream of the DQ124298:g.344A>G mutation. The difference in association analysis detected between both MUC4 SNPs indicates that they are not in linkage disequilibrium. The MUC4 gene successfully mapped within the confidence interval of the SSC13 QTL affecting NBA and TNB traits previously described by Noguera et al.  using the same resource population. The favourable effect of DQ124298:g.344A>G SNP was inherited from the Me maternal breed, as expected, since Me is more prolific than Ib. The effect of the DQ124298:g.344A>G MUC4 SNP was stronger for NBA than for TNB suggesting that the MUC4 gene could be related to the embryonic survival in the uterus. This hypothesis is supported by the differential expression of the MUC4 gene found in the uterus at 30-32 days of gestation, where the expression was two-fold higher in high than in low prolificacy sows. Nevertheless, further studies must be performed in other unrelated porcine populations to validate the association which was obtained in this study.
MUC4 has been identified as a potential regulator of placentation in pigs . In porcine uterine surface epithelium, MUC4 gene expression increases during the first 18 days of gestation  whereas in rats a reduction in MUC4 gene expression in the uterus is associated with the period of implantation . Differences in MUC4 expression profiles could be explained by the different placentation types of rodents and pigs. Pigs have a non-invasive type of placentation and MUC4 is thought to play a role in protecting the uterus from erosion by the embryo [13, 15]. MUC4 and MUC1 are the major mucin proteins expressed in the endometrial epithelium [22, 23]. Human MUC1 and MUC4 present highly polymorphic sites with a variable number of tandem repeats (VNTRs). MUC1 VNTR variants have been related to alterations in both T-antigen presentation and in the local immune response in cancer [17, 24, 25]. These results suggest that mucins may play a role in the immunological processes that take place during the implantation period essential to ensure the correct establishment of maternal-foetal tolerance .
MUC4 polymorphism DQ124298:g.344A>G is associated with litter size in the Ib × Me population and, moreover, MUC4 is differentially expressed regarding the number of embryos in uterus at 30 days of gestation. These results suggest that MUC4 may participate in the establishment of an optimal uterine environment essential for adequate embryo development during the early stages of gestation and increase litter size in pigs.
Animal material and phenotypic measurements
An F2 population was generated by crossing 3 Ib males from the Guadyerbas line (Dehesón del Encinar, Toledo, Spain) with 18 Me females (Domaine du Magneraud, INRA, France). A total of 8 boars and 97 sows from the F1 generation were mated to obtain the F2 progeny in the Nova Genètica S.A experimental farm (Lleida, Spain). All animals were obtained according to the European animal experimentation ethics law approved by the Ethical and Care Committee at IRTA.
Measurements of sow reproduction traits including TNB and NBA were recorded in 255 F2 sows during 4 successive parities. In the fifth parity, sows were slaughtered at 30-32 days of gestation, when the uterus samples were recollected and the number of embryos (NE) was recorded.
Genotyping of MUC4 polymorphisms
List of primer sequences used for typing MUC4 polymorphisms and quantitative PCR
Primer Sequence (5' →3')
QTL and association analyses
The linkage map of porcine chromosome 13 (SSC13) was constructed with the Build option of CRIMAP 2.4 software . Overall, seven markers were used: five microsatellites that had been previously described by , one SNP at the ITIH3 gene (dbSNP accession number ss315834911) and the DQ124298:g.344A>G MUC4 SNP described by .
Three models were used to analyse the Ib × Me F2 population data: (1) a QTL model, (2) an association model and (3) a QTL + association model to perform the marker assisted association test (MAAT) proposed by .
where y ijk was the phenotypic data vector for NBA or TNB; H i and O j were the fixed effects for year-season and parturition order, respectively; u k was the random polygenic effect of each individual; p k was the environmental permanent effect of the sow; a was the QTL additive effect; C a was the probability of the individual being homozygous for Ib alleles minus the probability of being homozygous for Me alleles at the QTL position of interest; and e ijkl was the random residual term. The dominance effect was not included in the model because the likelihood ratio test performed indicated that a model with only additive QTL effect fitted better.
where y ijk was the vector containing the phenotypic data for NBA and TNB and λ ik was a variable that indicated the number of copies (0, 1, or 2) of the kth allele presented by each individual.
All analyses were performed with Qxpack software . QTL scans were analysed every cM and nominal P-values were calculated with the maximum likelihood ratio test, assuming a χ2 distribution of the likelihood ratio test with degrees of freedom calculated as the difference between the number of parameters in the alternative and in the null models.
Expression analysis of the MUC4 gene
Uterine expression of the MUC4 gene was measured by reverse transcription quantitative real time PCR (RT-qPCR). Total RNA from uterus samples was extracted by means of the RiboPure™ kit (Ambion, Applied Biosystems). One microgram of total RNA in 40 μl reaction was reverse transcribed with the High Capacity cDNA Transcription Kit (Applied Biosystems). Primers for MUC4 (Table 3) were designed with Primer Express® 2.0 software (Applied Biosystems, Warrington, UK). Sus scrofa hypoxanthine phosphoribosyltransferase 1 (HPRT1) was used as a reference gene for normalization [for primer sequences see 30]. qPCR reactions were performed in triplicate in a 20 μl final volume including 2X FastStart SYBR Green Master (Roche), 0.3 μM of each primer and 5 μl of the cDNA diluted twenty times on an ABI PRISM® 7900HT sequence detection system (Applied Biosystems, Warrington, UK). Thermal conditions were 95°C for 10 min followed by 40 cycles of 95°C for 15 s and 60°C for 1 min. Dissociation curve analyses were performed in order to detect unspecific amplifications. PCR efficiencies of MUC4 and HPRT1 genes were calculated to validate the use of the 2-ΔΔCt method . The MUC4 gene expression level was measured individually in 36 Ib × Me F2 sows from 32 different litters that were classified into two groups according to the NE at slaughter (30-32 days of pregnancy): low (NE ≤ 11, n = 20) and high (NE ≥ 13, n = 16). Gene expression data were log10 transformed and analysed by a t-test using the Statistical Analysis System (Statistics, V 9.1.3; SAS Institute, Inc., Cary, NC). The significance threshold was set at α < 0.05.
The authors are indebted to L. Varona, M. Arqué, J. Tarrés, M. Fina, and the staff of Nova Genètica, particularly E. Ramells, F. Márquez, R. Malé, F. Rovira, and I. Riart, for cooperating in the experimental protocol. The authors gratefully acknowledge the contributions of the INRA (France) and the CIA El Dehesón del Encinar (Spain) in providing the purebred Meishan sows and Iberian boars, respectively. The authors also thank J. Casellas, V. Martínez and S. Cirera for critical comments on the manuscript. This research was funded in part by Project AGL2004-08368-C03 and by the Consolider-Ingenio 2010 Programme (CSD2007-00036), both from the Spanish Ministry of Science and Innovation. IB is the recipient of a PIF PhD fellowship from Universitat Autònoma de Barcelona.
- Spotter A, Distl O: Genetic approaches to the improvement of fertility traits in the pig. Vet J. 2006, 172: 234-247. 10.1016/j.tvjl.2005.11.013.View ArticlePubMedGoogle Scholar
- Tribout T, Iannuccelli N, Druet T, Gilbert H, Riquet J, Gueblez R, Mercat MJ, Bidanel JP, Milan D, Le Roy P: Detection of quantitative trait loci for reproduction and production traits in Large White and French Landrace pig populations. Genet Sel Evol. 2008, 40: 61-78.PubMed CentralPubMedGoogle Scholar
- Noguera JL, Rodriguez C, Varona L, Tomas A, Munoz G, Ramirez O, Barragan C, Arque M, Bidanel JP, Amills M, Ovilo C, Sanchez A: A bi-dimensional genome scan for prolificacy traits in pigs shows the existence of multiple epistatic QTL. BMC Genomics. 2009, 10: 636-10.1186/1471-2164-10-636.PubMed CentralView ArticlePubMedGoogle Scholar
- Fernandez-Rodriguez A, Rodriguez C, Varona L, Balcells I, Noguera JL, Ovilo C, Fernandez AI: Analysis of candidate genes underlying two epistatic quantitative trait loci on SSC12 affecting litter size in pig. Anim Genet. 2010, 41: 73-80. 10.1111/j.1365-2052.2009.01966.x.View ArticlePubMedGoogle Scholar
- Moniaux N, Nollet S, Porchet N, Degand P, Laine A, Aubert JP: Complete sequence of the human mucin MUC4: a putative cell membrane-associated mucin. Biochem J. 1999, 338: 325-333. 10.1042/0264-6021:3380325.PubMed CentralView ArticlePubMedGoogle Scholar
- Lagow E, DeSouza MM, Carson DD: Mammalian reproductive tract mucins. Hum Reprod Update. 1999, 5: 280-292. 10.1093/humupd/5.4.280.View ArticlePubMedGoogle Scholar
- Govindarajan B, Gipson IK: Membrane-tethered mucins have multiple functions on the ocular surface. Exp Eye Res. 2010, 90: 655-663. 10.1016/j.exer.2010.02.014.PubMed CentralView ArticlePubMedGoogle Scholar
- McNeer RR, Carraway CA, Fregien NL, Carraway KL: Characterization of the expression and steroid hormone control of sialomucin complex in the rat uterus: implications for uterine receptivity. J Cell Physiol. 1998, 176: 110-119. 10.1002/(SICI)1097-4652(199807)176:1<110::AID-JCP13>3.0.CO;2-B.View ArticlePubMedGoogle Scholar
- McNeer RR, Huang D, Fregien NL, Carraway KL: Sialomucin complex in the rat respiratory tract: a model for its role in epithelial protection. Biochem J. 1998, 330: 737-744.PubMed CentralView ArticlePubMedGoogle Scholar
- Idris N, Carraway KL: Sialomucin complex (Muc4) expression in the rat female reproductive tract. Biol Reprod. 1999, 61: 1431-1438. 10.1095/biolreprod61.6.1431.View ArticlePubMedGoogle Scholar
- Idris N, Carraway KL: Regulation of sialomucin complex/Muc4 expression in rat uterine luminal epithelial cells by transforming growth factor-beta: implications for blastocyst implantation. J Cell Physiol. 2000, 185: 310-316. 10.1002/1097-4652(200011)185:2<310::AID-JCP16>3.0.CO;2-W.View ArticlePubMedGoogle Scholar
- Carraway KL, Idris N: Regulation of sialomucin complex/Muc4 in the female rat reproductive tract. Biochem Soc Trans. 2001, 29: 162-166. 10.1042/BST0290162.View ArticlePubMedGoogle Scholar
- Ferrell AD, Malayer JR, Carraway KL, Geisert RD: Sialomucin complex (Muc4) expression in porcine endometrium during the oestrous cycle and early pregnancy. Reprod Domest Anim. 2003, 38: 63-65. 10.1046/j.1439-0531.2003.00405.x.View ArticlePubMedGoogle Scholar
- Komatsu M, Carraway CA, Fregien NL, Carraway KL: Reversible disruption of cell-matrix and cell-cell interactions by overexpression of sialomucin complex. J Biol Chem. 1997, 272: 33245-33254. 10.1074/jbc.272.52.33245.View ArticlePubMedGoogle Scholar
- Carraway KL, Perez A, Idris N, Jepson S, Arango M, Komatsu M, Haq B, Price-Schiavi SA, Zhang J, Carraway CA: Muc4/sialomucin complex, the intramembrane ErbB2 ligand, in cancer and epithelia: to protect and to survive. Prog Nucleic Acid Res Mol Biol. 2002, 71: 149-185.View ArticlePubMedGoogle Scholar
- Ross JW, Ashworth MD, Stein DR, Couture OP, Tuggle CK, Geisert RD: Identification of differential gene expression during porcine conceptus rapid trophoblastic elongation and attachment to uterine luminal epithelium. Physiol Genomics. 2009, 36: 140-148.View ArticlePubMedGoogle Scholar
- Chang CY, Chang HW, Chen CM, Lin CY, Chen CP, Lai CH, Lin WY, Liu HP, Sheu JJ, Tsai FJ: MUC4 gene polymorphisms associate with endometriosis development and endometriosis-related infertility. BMC Med. 2011, 9: 19-10.1186/1741-7015-9-19.PubMed CentralView ArticlePubMedGoogle Scholar
- Koscinski I, Viville S, Porchet N, Bernigaud A, Escande F, Defossez A, Buisine MP: MUC4 gene polymorphism and expression in women with implantation failure. Hum Reprod. 2006, 21: 2238-2245. 10.1093/humrep/del189.View ArticlePubMedGoogle Scholar
- Peng QL, Ren J, Yan XM, Huang X, Tang H, Wang YZ, Zhang B, Huang LS: The g.243A>G mutation in intron 17 of MUC4 is significantly associated with susceptibility/resistance to ETEC F4ab/ac infection in pigs. Anim Genet. 2007, 38: 397-400. 10.1111/j.1365-2052.2007.01608.x.View ArticlePubMedGoogle Scholar
- Zhao H, Rothschild MF, Fernando RL, Dekkers JC: Tests of candidate genes in breed cross populations for QTL mapping in livestock. Mamm Genome. 2003, 14: 472-482. 10.1007/s00335-002-2215-y.View ArticlePubMedGoogle Scholar
- Ostrup E, Bauersachs S, Blum H, Wolf E, Hyttel P: Differential endometrial gene expression in pregnant and nonpregnant sows. Biol Reprod. 2010, 83: 277-285. 10.1095/biolreprod.109.082321.View ArticlePubMedGoogle Scholar
- Gollub EG, Goswami S, Kouba D, Marom Z: Regulation of mucin gene expression in secretory epithelial cells. Biochem Biophys Res Commun. 1993, 197: 667-673. 10.1006/bbrc.1993.2531.View ArticlePubMedGoogle Scholar
- Audie JP, Tetaert D, Pigny P, Buisine MP, Janin A, Aubert JP, Porchet N, Boersma A: Mucin gene expression in the human endocervix. Hum Reprod. 1995, 10: 98-102. 10.1093/humrep/10.1.98.View ArticlePubMedGoogle Scholar
- Fowler JC, Teixeira AS, Vinall LE, Swallow DM: Hypervariability of the membrane-associated mucin and cancer marker MUC1. Hum Genet. 2003, 113: 473-479. 10.1007/s00439-003-1011-8.View ArticlePubMedGoogle Scholar
- Santos-Silva F, Fonseca A, Caffrey T, Carvalho F, Mesquita P, Reis C, Almeida R, David L, Hollingsworth MA: Thomsen-Friedenreich antigen expression in gastric carcinomas is associated with MUC1 mucin VNTR polymorphism. Glycobiology. 2005, 15: 511-517.View ArticlePubMedGoogle Scholar
- Guerin LR, Prins JR, Robertson SA: Regulatory T-cells and immune tolerance in pregnancy: a new target for infertility treatment?. Hum Reprod Update. 2009, 15: 517-535. 10.1093/humupd/dmp004.PubMed CentralView ArticlePubMedGoogle Scholar
- Ronaghi M, Pettersson B, Uhlen M, Nyren P: PCR-introduced loop structure as primer in DNA sequencing. BioTechniques. 1998, 25: 876-8. 880-2, 884PubMedGoogle Scholar
- Green P, Falls K, Crooks S: Documentation of CRI-MAP. [http://compgen.rutgers.edu/old/multimap/crimap/index.html]
- Perez-Enciso M, Misztal I: Qxpak: a versatile mixed model application for genetical genomics and QTL analyses. Bioinformatics. 2004, 20: 2792-2798. 10.1093/bioinformatics/bth331.View ArticlePubMedGoogle Scholar
- Fernandes LT, Tomas A, Bensaid A, Perez-Enciso M, Sibila M, Sanchez A, Segales J: Exploratory study on the transcriptional profile of pigs subclinically infected with porcine circovirus type 2. Anim Biotechnol. 2009, 20: 96-109. 10.1080/10495390902885785.View ArticlePubMedGoogle Scholar
- Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2001, 25: 402-408. 10.1006/meth.2001.1262.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.