<title>Background</title> <p>Differentially expressed genes in the early stages of muscle development may be potential candidate genes to improve meat quality and quantity [##UREF##0##1##,##REF##12968692##2##]. Recent studies showed that differential expression of genes during muscle development in prenatal pigs is associated with the differential stages of myogenesis [##REF##15999225##3##,##UREF##1##4##], and expression patterns of these genes differ between pig breeds with different muscle characteristics [##REF##16361485##5##,##REF##17199897##6##].</p> <p>An EST [GenBank: <ext-link ext-link-type="gen" xlink:href="AA063650">AA063650</ext-link>] was identified showing differential expression during the early stages of muscle development [##UREF##2##7##]. In order to identify this EST with expression associated with prenatal muscle development, a full-length cDNA library was constructed using 55-day Chinese Tongcheng pig fetus skeletal muscle [##UREF##2##7##]. From the cDNA library, a 2285-bp cDNA clone was obtained which was identical with this EST. After BLAST search in the non-redundant (nr) sequence database, it was found to be homologous to a clone DKFZp586C1019 [GenBank: <ext-link ext-link-type="gen" xlink:href="AL049397">AL049397</ext-link>] which is part of Homo sapiens <italic>PNAS-4 </italic>gene (also called <italic>C1orf121</italic>) encoding the CGI-146 protein. Recently, a report showed that human <italic>PNAS-4 </italic>is a novel pro-apoptotic protein activated during the early response to DNA damage [##UREF##3##8##]. Furthermore, this gene was also predicted as one of the targets of p53 tumor suppressor via mathematical modeling and quantitative data analysis [##UREF##4##9##]. Presently, Xenopus laevis <italic>PNAS-4 </italic>protein has been purified and confirmed by Western blot [##REF##17630912##10##]. However, the biological functions of <italic>PNAS-4</italic>, especially the role in muscle development, are still unclear.</p> <p>The aim of the present study was to characterize the porcine <italic>PNAS-4 </italic>gene by obtaining its sequence, and gain insight into its potential physiological role in muscle development related to meat production by studying subcellular localization, the temporal prenatal expression pattern, postnatal tissue expression levels and trait association analysis.</p>
<title>Methods</title> <title>cDNA isolation and sequence analysis</title> <p>BLAST in the pig expressed sequence tags (ESTs) database using the human <italic>PNAS-4 </italic>mRNA sequence [GenBank: <ext-link ext-link-type="gen" xlink:href="NM 016076">NM 016076</ext-link>] as probe, revealed three porcine ESTs [GenBank: <ext-link ext-link-type="gen" xlink:href="CN156918">CN156918</ext-link>, <ext-link ext-link-type="gen" xlink:href="CV877591">CV877591</ext-link> and <ext-link ext-link-type="gen" xlink:href="AJ665351">AJ665351</ext-link>]. To obtain the full-length cDNA, primer set GAPF/GAPR (Table ##TAB##2##3##) was used to amplify the unknown nucleotide sequence between ESTs <ext-link ext-link-type="gen" xlink:href="CV877591">CV877591</ext-link> and <ext-link ext-link-type="gen" xlink:href="AJ665351">AJ665351</ext-link>, and the primer 5GSP (Table ##TAB##2##3##) in combination with the Universal Primer A Mix (UPM) supplied in the SMART™ RACE Kit (Clontech Inc, Palo Alto, CA, USA) used to amplify the 5'-end of the gene. Rapid amplification of cDNA ends (RACE) PCR was carried out in a PTC-200 thermal cycler (Bio-Rad, Richmond, CA, USA) with touchdown PCR conditions: 3 min at 95°C; 5 cycles of 94°C for 20 s, 72°C for 3 min; 5 cycles of 94°C for 20 s, 70°C for 30 s; 72°C for 3 min; 32 cycles of 94°C for 20 s, 68°C for 30 s; 72°C for 3 min; 5 min at 72°C. The PCR products were purified by agarose gel electrophoresis and sequenced. Functional elements analysis of the putative protein and UTR were performed using the PROSITE program[##UREF##8##27##], SignalP 3.0 server[##UREF##9##28##] and UTRScan[##UREF##10##29##], respectively.</p> <title>QPCR analysis of gene expression patterns</title> <p>Fifteen tissues (stomach, large intestine, kidney, lymph node of neck, small intestine, liver, heart, longissimus dorsi muscle, lung, gastrocnemius, biceps femoris muscle, brain, backfat, semitendinosus and spleen) were obtained from three Chinese indigenous Wuzhishan (a Chinese mini pig breed) sows. Three pig embryos were collected from every pregnant female during three embryonic periods (33, 65, and 90 dpc) and three postnatal periods (2, 28 days and adult (160 days)). The longissimus dorsi muscle were harvested, mixed, frozen in liquid nitrogen, and stored at -80°C. Furthermore, three Landrace pigs were employed to determine expression differences between western and Chinese pig breeds in prenatal stages. All animal procedures were performed according to guidelines developed by the China Council on Animal Care and protocol approved by Animal Care and Use Committee of Hubei Province, PR China.</p> <p>QPCR was performed with the Chromo4™ Real-Time PCR Detector (Bio-Rad) using gene-specific primers and Taqman probes (Table ##TAB##2##3##). Total RNA (4 <italic>μ</italic>g), extracted from tissues using the Trizol reagent (Invitrogen, Carlsbad, CA, USA) and treated with RNase-free DNase (MBI Fermentas, St. Leon-Rot, Germany). RNA quality was assessed with agarose gel electrophoresis and the ND-1000 Spectrophotometer (NanoDrop Technologies, USA), and was discarded if the 260/280 ratio was not between 1.8 and 2.1. It was reverse-transcribed into cDNA using M-MLV reverse transcriptase (Promega Corp., Madison, WI, USA). cDNA synthesis was performed in duplicate. Both beta-actin gene (<italic>ACTB</italic>) and glyceraldehyde-3-phosphate dehydrogenase (<italic>GAPDH</italic>) were consistently expressed in all tissues of this assay, so we selected the <italic>ACTB</italic>, one of the two housekeeping genes, as an internal control for normalization purposes. The minus RT control with primers for <italic>ACTB </italic>was performed to test if products were derived from DNA. Relative transcript quantification was performed using standard curves generated for <italic>ACTB </italic>and <italic>PNAS-4 </italic>gene from a 10-fold serial dilution of cDNA. Pooled cDNA from a subset of the gastrocnemius samples examined in this study was used to generate the standard curves. In this assay, the efficiency of <italic>ACTB </italic>and <italic>PNAS-4 </italic>gene primers were 96.5% and 96.95%, respectively. The cycling conditions consisted of an initial, single cycle of 95°C for 3 min, followed by 35 cycles of 15 s at 94°C and 1 min at 60°C. Each reaction (in 20 <italic>μ</italic>l) contained 1 × PCR buffer (TaKaRa, Dalian, P. R. C.), 0.4 <italic>μ</italic>M gene-specific primers and 0.2 <italic>μ</italic>M probe, 50 <italic>μ</italic>M of each dNTP, 3 mM MgCl<sub>2</sub>, 2.0 U Taq DNA polymerase (TaKaRa) plus 2 <italic>μ</italic>l normalized template cDNA. All PCR amplifications were performed in triplicate for each RNA sample and gene expression levels were quantified relative to <italic>ACTB </italic>expression using Gene Expression Macro software (Bio-Rad). The results were analyzed using the 2<sup>-ΔΔCt </sup>method described previously[##REF##11846609##30##]. Differences in gene expression between groups were evaluated using Student's t-test and were considered statistically significant at <italic>P </italic>&lt; 0.05.</p> <title>Intracellular distribution of the porcine <italic>PNAS-4 </italic>gene</title> <p>The intracellular distribution of porcine <italic>PNAS-4 </italic>protein was studied by fluorescence and confocal analysis of pig kidney epithelial cells (PK15) cells transiently transfected with P4-pEGFP-C1 construct. The ORF, encoding porcine <italic>PNAS-4 </italic>protein, was amplified by PCR using the primers (P4GFPF/P4GFPR, Table ##TAB##2##3##) designed from its cDNA sequence and subcloned into the <italic>Xho </italic>I-<italic>Hind </italic>III site of the pEGFP-C1 vector (BD Biosciences Clontech, CA, USA).</p> <p>PK15 cells were cultured in Dulbecco's modified Eagle's medium (DMEM) supplemented with 10% fetal calf serum (GIBCO-BRL, Grand Island, NY), 4 mM glutamine, 100 U/ml penicillin, and 0.1 mg/ml streptomycin under humidified air containing 5% CO<sub>2 </sub>at 37°C and seeded onto cover slips in 6-well plates. The PK15 cells were plated at 5 × 10<sup>5</sup>/ml and transfected with P4-pEGFP-C1/N3 constructs (1 <italic>μ</italic>g of plasmid DNA) using Lipofectamine™ 2000 reagent according to the manufacturer's instructions until they were grown to 80–90% confluence.</p> <p>Twenty-four hours after transfection, the cells were incubated at 37°C for 30 min in growth medium containing 200 nM MitoTracker Red CM-H<sub>2</sub>Xros (Molecular Probes, Eugene, Oregon, USA) for mitochondrial labeling. Then the cells were fixed with pre-warmed growth medium containing 3.7% formaldehyde for 15 min at 37°C. After the final washing steps and incubation with 10 <italic>μ</italic>M Hoechst33342 for 10 min, the slides were mounted, sealed and analyzed by Olympus FluoView™ FV1000 Confocal Microscope. The software FV1000 Viewer was used to generate individual fluorescent pictures as well as overlay pictures that demonstrated the relative distribution of the fusion protein.</p> <title>Chromosome mapping by SCH and IMpRH</title> <p>The French SCH panel composed of 27 pig/rodent hybrids [##REF##8697807##31##] was used for the regional assignment of <italic>PNAS-4 </italic>gene on the SSC. To more precisely map genes, the IMpRH panel consisting of 118 hybrid clones obtained by fusing irradiated pig cells with the recipient hamster cell was employed [##REF##9858812##32##]. The PCR was performed in a total volume of 10 <italic>μ</italic>L 1 × PCR buffer (TaKaRa), containing 25 ng panel DNA, 0.3 <italic>μ</italic>M each primer (MapF, MapR, Table ##TAB##2##3##), 75 <italic>μ</italic>M each dNTP, 1.5 mM MgCl<sub>2</sub>, and 1.0 U Taq DNA polymerase (TaKaRa). The PCR cycling conditions were: Initial denaturation at 95°C for 3 min, followed by 32 cycles (20 s at 94°C, 30 s at 65°C and 30 s at 72°C). Statistical analysis of SCH PCR results were performed using the program [##REF##9088711##33##]available for regional assignment of genetic markers. The results of radiation hybrid PCR products were analysed with the IMpRH mapping tool [##REF##10980153##34##].</p> <title>Polymorphism detection and association with meat traits</title> <p>Four Chinese indigenous pig breeds (Xiang (N = 42), Bamaxiang(N = 45), Wuzhishan (N = 43) and Tongcheng (N = 44) pigs) and two western white breeds (Large White (N = 215) and Landrace (N = 15)) were used to analyze genetic variability. Association analyses were performed in our experimental populations that contained three pure breeds, Tongcheng (N = 44), Landrace (N = 15) and Large White<sup>a </sup>(N = 19), and a pair of three-breed cross, Landrace × (Large White × Tongcheng) (N = 44) and Large White × (Landrace × Tongcheng) (N = 35). Growth traits (average daily gain), carcass and meat quality traits (live weight at slaughter, carcass length, dressing percent, loin muscle area, meat color, marbling, drip loss, 0 h loin muscle PH and shear force) were studied. Fat traits studied were percentage of leaf fat, percentage of leaf and caul fat, and backfat thickness between 6th and 7th ribs.</p> <p>Two pairs of primers (SNPF1/SNPR1 located in intron 2 and SNPF2/SNPR2 located in 3'UTR) were designed using the submitted sequences [GenBank: <ext-link ext-link-type="gen" xlink:href="DQ406743">DQ406743</ext-link> and <ext-link ext-link-type="gen" xlink:href="DQ435075">DQ435075</ext-link>]. The PCR products were cloned into the T-easy vector and sequenced to identify the polymorphisms in the porcine <italic>PNAS-4 </italic>gene.</p> <p>Associations between growth traits, carcass traits and meat quality traits on the one hand and genotype on the other hand were studied using a two-step analysis method. Briefly, in the first step, a preliminary general linear model without the genotype information was used to eliminate system effects, including sex, combination and batch, using the SPSS program with the following model:</p> <p></p> <p>where <italic>Y</italic><sub><italic>ijkl </italic></sub>is phenotypic value of the target trait, <italic>μ </italic>is the population mean, <italic>C</italic><sub><italic>i </italic></sub>is the combination (including breed) effect, <italic>B</italic><sub><italic>j </italic></sub>is the batch effect, <italic>S</italic><sub><italic>k </italic></sub>is the sex effect and <italic>e</italic><sub><italic>ijkl </italic></sub>is random error effect for each observation. In addition, the interaction effects between combination and batch (<italic>CB</italic>)<sub><italic>ij</italic></sub>, and between batch and sex (<italic>BS</italic>)<sub><italic>jk</italic></sub>, were also considered. In the second step, the resulting standardized residual values were used to partition the effects of the genotypes using the one-way analysis of variance (ANOVA) model. Simultaneously, significant difference tests between different genotypes were performed with Bonferroni <italic>t </italic>tests [##UREF##11##35##].</p>
<title>Results</title> <title>Molecular cloning and sequence analysis</title> <p>The full-length cDNA of porcine <italic>PNAS-4 </italic>gene consists of 4,059 bp that contains an open reading frame (ORF) of 582-bp encoding a protein of 194 amino acids with a calculated molecular mass of 21.4 kDa and an isolectric point (p<italic>I</italic>) of 4.807 (Figure ##FIG##0##1##). Two polyadenylation signals (AATAAA) were observed at positions 821 to 826 and 2307 to 2312 upstream of the poly (A) stretch. The Brd-Box (AGCTTTA) and GY-Box (GTCTTCC) were found in the 3' untranslated region (UTR) (Figure ##FIG##0##1##). The mRNA sequence of this gene has been submitted to GenBank [GenBank: <ext-link ext-link-type="gen" xlink:href="DQ435075">DQ435075</ext-link>].</p> <p>Comparison of the porcine amino acid sequence with the proteins of seven species reported in GenBank shows that the porcine putative protein possesses a high level of identity (94–98%) with mammals and fowls, and it shares 74% identity with Xenopus laevis throughout the protein sequence. The porcine protein contained several conserved motifs such as C-terminal microbody targeting signal (CMTS), N-myristoylation and a Protein kinase C phosphorylation site (Figure ##FIG##1##2##). In addition, it showed that the porcine <italic>PNAS-4 </italic>encoded a non-secretory protein without signal peptide and transmembrane regions.</p> <title>Spatial and temporal expression pattern</title> <p>The relative expression of the porcine <italic>PNAS-4 </italic>transcript in various pig tissues showed highest expression in skeletal muscle (in decreasing order: gastrocnemius, semitendinosus, longissimus dorsi and biceps femoris) and lymph node. Lower expression levels were determined in the liver, large intestine, backfat, heart, small intestine, stomach, spleen, lung, kidney and brain (Figure ##FIG##2##3##).</p> <p>Figure ##FIG##3##4## showed that the expression level of porcine <italic>PNAS-4 </italic>seemed to decrease gradually during embryonic development of Landrace and Tongcheng pig breeds, which was significant in Landrace pigs between 33 days post conception (dpc) and 90 dpc. This was consistent with the results of semi-quantitative RT-PCR [##UREF##2##7##]. Curiously, after birth, the amount of porcine <italic>PNAS-4 </italic>mRNA in skeletal muscles increased drastically (<italic>P </italic>&lt; 0.01). However, the mRNA expression on the day of weaning (28-day neonate) was lower than the second day after birth and adult (<italic>P </italic>&lt; 0.05). After maturity, the mRNA expression of porcine <italic>PNAS-4 </italic>reached a peak level (<italic>P </italic>&lt; 0.01).</p> <title>Intracellular distribution</title> <p>The CGI-146-GFP fusion proteins localized in the Golgi complex (Figure ##FIG##4##5##), which was identified by the complementation tests using the homo-functional P4-pEGFP-N3 construct. Green fluorescence was detected in the cytoplasm of the control cells transfected with pEGFP vector alone (E-H, Figure ##FIG##4##5##).</p> <title>Chromosome assignment</title> <p>The somatic cell hybrid (SCH) panel revealed that the porcine <italic>PNAS-4 </italic>gene was assigned to SSC10 (q11–q16) (probability of regional localization: 0.9998, error risk &lt; 0.1%). Two-point analysis of INRA-University of Minnesota porcine radiation hybrid (IMpRH) panel revealed that the porcine <italic>PNAS-4 </italic>was closely linked to SW497 (LOD = 13.68, 28 cRs) with 22% retention frequencies. This result was consistent with the SCH panel. The corresponding conserved syntenic region on human chromosome was the end of Homo sapiens chromosome 1 long arm (HSC1q)[##REF##16246521##11##].</p> <title>Polymorphism detection and association with economic traits</title> <p>Comparative sequencing revealed 9 and 5 polymorphisms within two PCR amplified fragments [GenBank: <ext-link ext-link-type="gen" xlink:href="DQ406743">DQ406743</ext-link> and <ext-link ext-link-type="gen" xlink:href="DQ435075">DQ435075</ext-link>] which were parts of <italic>PNAS-4 </italic>genomic DNA. Because of linkage disequilibrium among several polymorphisms, two single nucleotide polymorphisms (SNPs) were selected for further study. The T1400C SNP within intron 2 was detectable after digestion with <italic>Sty </italic>I, resulting in a 762-bp PCR amplicon produced allele T (762 bp) and allele C (549 bp and 213 bp). The other SNP, A1813C within the 3'UTR, located more than 20 kb away from the T1400C SNP in genome, was harbored in a 789-bp PCR amplicon and digested by <italic>Msp </italic>I, resulting in allele A (789 bp) and allele C (606 bp and 183 bp).</p> <p>Analysis of polymorphism revealed that allele C for both polymorphisms (T1400C and A1813C) was prevalent in four Chinese indigenous pig breeds whereas the allele T (T1400C) and the allele A (A1813C) respectively had higher frequencies in Landrace and Large White (Table ##TAB##0##1##). The two SNPs could be linked in Landrace, Large White pigs and a pair of three breed-crosses (data not shown). Therefore, we selected one of the two SNPs for association analysis. Three significant associations were observed between the A1813C SNP genotypes and three fat-related traits (Table ##TAB##1##2##). For the three considered traits, homozygous CC pigs of our experimental population showed higher values than those of the pigs with other genotypes (<italic>P </italic>&lt; 0.05).</p>
<title>Discussion</title> <p>Differential expression of genes during a biological process suggests a direct or indirect relationship between the genes and the biological process[##UREF##5##12##]. An EST [GenBank: <ext-link ext-link-type="gen" xlink:href="AA063650">AA063650</ext-link>] was identified showing differential expression during the early stages of muscle development [##UREF##2##7##]. Therefore, the full-length cDNA of the <italic>PNAS-4 </italic>gene was obtained and analyzed for study its structure and biological functionality. This manuscript reports the cDNA and amino acid sequences of the gene, its mRNA expression profiles, protein localization, and functional analysis in muscle tissue formation and meat production.</p> <title><italic>PNAS-4 </italic>cDNA and amino acid analysis</title> <p>Full-length cDNA sequence analysis revealed high identity of the gene at both the mRNA and the amino acid level with mammalian and avian species. This may indicate both structural and functional conservation of the gene. The Brd-Box and GY-Box within the 3'UTR suggests that the <italic>PNAS-4 </italic>expression was regulated negatively at the post-transcriptional level, both temporal and spatial [##REF##9735368##13##, ####REF##11896390##14##, ##REF##9428421##15####9428421##15##].</p> <title>Functional analysis in muscle tissue related to meat production</title> <p><italic>PNAS-4 </italic>mRNA showed a widespread expression with the highest levels in muscle tissue, and could be detected at all times during the skeletal muscle development. In more detail, the expression of <italic>PNAS-4 </italic>tends to decrease gradually from prenatal day 33 to day 90, an important period for development of muscle fibers and determination of muscle fiber numbers [##UREF##0##1##,##REF##4748487##16##,##UREF##6##17##]. The muscle fiber number appears to increase as a result of the decrease of <italic>PNAS-4 </italic>mRNA expression level during prenatal muscle development, which was consistent with that the <italic>PNAS-4 </italic>is involved in apoptosis[##UREF##3##8##]. However, the expression profiles of Tongcheng pigs and Landrace pigs were not complete similar in early embryonic development, and this kind of difference may be related to breed difference [##REF##16361485##5##]. The <italic>PNAS-4 </italic>mRNA expression decreased faster in early Landrace pig development, and this result may be correlated with the high number of primary fibers found in Landrace compared with Tongcheng pigs and other indigenous Chinese breeds [##REF##17573972##18##, ####REF##8145210##19##, ##REF##7608006##20##, ##REF##8667335##21####8667335##21##]. In other words, there might be a positive correlation between the decrease speed of <italic>PNAS-4 </italic>mRNA expression level in embryonic stage and porcine muscle fiber number. Muscle fiber number is an important determinant of postnatal growth such that pigs with a high fiber number tend to grow faster than pigs with a lower fiber number [##REF##8294285##22##]. This is also supported by the carcass weight of Landrace compared to Tongcheng pigs[##REF##17573972##18##]. After birth, the expression of <italic>PNAS-4 </italic>mRNA in muscle tissue was increased, which could be attributed to an increase in length and girth of the muscle fibres and/or muscle function [##UREF##0##1##]. Consistent with this, the highest expression of the gene was recorded in skeletal muscle of mature pigs. The protein was found to localize in the Golgi complex. It was thus deduced that <italic>PNAS-4 </italic>protein acted more possibly as a regulator than as a structural protein, which was generally in agreement with the opinion that <italic>PNAS-4 </italic>may be a pro-apoptotic gene [##UREF##3##8##].</p> <title><italic>PNAS-4 </italic>as a candidate gene or QTL for fat deposition</title> <p>Porcine <italic>PNAS</italic>-4 gene mapped to q11–16 of SSC10, closely linked to marker SW497. It showed that several quantitative trait loci (QTL) related to growth and carcass traits, such as average daily gain (ADG), Backfat Thickness (BFT), Carcass Weight (CCWT), Ham Weight (HAMWT), are located on this chromosome region in the Pig QTL Database release [##UREF##7##23##, ####REF##12081807##24##, ##REF##16415521##25####16415521##25##]. Consistent with this, significant associations between A1813C SNP genotypes of <italic>PNAS-4 </italic>and three fat-related traits were found in our experimental pig populations. Moreover, the three fat traits values of pigs with homozygous CC were higher than that of the pigs with other genotypes (<italic>P </italic>&lt; 0.05). Our findings have also revealed that allele C of SNP A1813C was prevalent in four Chinese indigenous breeds whereas allele A had higher frequency in Landrace and Large White. Interestingly, western commercial pig breeds (Duroc, Large White and Landrace) are well known for their high growth rate and lean meat percentage [##REF##1644675##26##]. These correlated results suggest that the pigs with allele C in SNP A1813C will have higher fat content. In order to confirm it, further investigations in the use of the gene for commercial breeding are recommended.</p>
<title>Conclusion</title> <p>We cloned and characterized the full-length cDNA of porcine <italic>PNAS-4 </italic>gene. Expression profiles and gene structure suggest that this gene was regulated negatively at the post-transcriptional level, both temporal and spatial. It might play a negative role in modulating muscle fiber numbers, and act as an important factor affecting the muscle mass efficiency. Unexpectedly, association analysis revealed that allele C of SNP A1813C was prevalent in four Chinese indigenous breeds whereas allele A had higher frequency in Landrace and Large White, and the pigs with allele C in SNP A1813C will have higher fat content. It is a potential genetic marker for animal breeding.</p>
<p>This is an Open Access article distributed under the terms of the Creative Commons Attribution License (<ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/2.0"/>), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.</p> <title>Background</title> <p>In a previous screen to identify differentially expressed genes associated with embryonic development, the porcine <italic>PNAS-4 </italic>gene had been found. Considering differentially expressed genes in early stages of muscle development are potential candidate genes to improve meat quality and production efficiency, we determined how porcine <italic>PNAS-4 </italic>gene regulates meat production. Therefore, this gene has been sequenced, expression analyzed and associated with meat production traits.</p> <title>Results</title> <p>We cloned the full-length cDNA of porcine <italic>PNAS-4 </italic>gene encoding a protein of 194 amino acids which was expressed in the Golgi complex. This gene was mapped to chromosome 10, q11–16, in a region of conserved synteny with human chromosome 1 where the human homologous gene was localized. Real-time PCR revealed that <italic>PNAS-4 </italic>mRNA was widely expressed with highest expression levels in skeletal muscle followed by lymph, liver and other tissues, and showed a down-regulated expression pattern during prenatal development while a up-regulated expression pattern after weaning. Association analysis revealed that allele C of SNP A1813C was prevalent in Chinese indigenous breeds whereas A was dominant allele in Landrace and Large White, and the pigs with homozygous CC had a higher fat content than those of the pigs with other genotypes (<italic>P </italic>&lt; 0.05).</p> <title>Conclusion</title> <p>Porcine <italic>PNAS-4 </italic>protein tagged with green fluorescent protein accumulated in the Golgi complex, and its mRNA showed a widespread expression across many tissues and organs in pigs. It may be an important factor affecting the meat production efficiency, because its down-regulated expression pattern during early embryogenesis suggests involvement in increase of muscle fiber number. In addition, the SNP A1813C associated with fat traits might be a genetic marker for molecular-assisted selection in animal breeding.</p>
<title>Abbreviations</title> <p>dpc: days post conception, ESTs: expressed sequence tags, GFP: green fluorescent protein, HSC: homo sapiens chromosome, IMpRH: INRA University of Minnesota porcine radiation hybrid, ORF: open reading frame, PK15: Pig Kidney Epithelial cells, QTL: quantitative trait loci, RACE: rapid amplification of cDNA ends, SCH: somatic cell hybrid, SNPs: single nucleotide polymorphisms, SSC: sus scrofa chromosome, UTR:  untranslated region.</p> <title>Authors' contributions</title> <p>DM performed research and drafted the manuscript. ZM participated in the design and the realization of the experiment. MFWtP provided DNA samples for SNP genotyping, and finalized the writing. XL participated in cell culture and subcellular localization. SY was responsible for the data management. HW contributed to bioinformatics analysis. HLW participated in SNP identification. KL conceived and coordinated the study and corrected the manuscript.</p>
<title>Acknowledgements</title> <p>We are grateful to Dr. Martine Yerle for providing the RH panel. This research was supported by National Natural Science Foundation of China (30400311), Guangdong Natural Science Foundation (07300492), National 973 Project (2006CB102105) and National 863 Project (2006AA10Z135), the State Platform of Technology infrastructure (2005DKA21101), the Project of Science and Technology Innovation Team "Research and Improvement of Domestic Animal Germplasm", Institute of Animal Science, Chinese Academy of Agricultural Science (ywf-td-1).</p>
<fig position="float" id="F1"><label>Figure 1</label><caption><p><bold>Nucleotide and deduced amino acid sequence of porcine <italic>PNAS-4</italic></bold>. The polyadenylation signals (AATAAA) are double-underlined; and the two elements named Brd-Box and GY-Box are boxed, separately. Primer site for 5' RACE assays is single underlined.</p></caption></fig> <fig position="float" id="F2"><label>Figure 2</label><caption><p><bold>Multiple amino acid sequence alignment of porcine CGI-146</bold>. Different amino acid residues among species are presented with boxes. The three highly conserved motifs (N-myristoylation, Protein kinase C phosphorylation and MCTS) are indicated by roman letters (I – III), respectively.</p></caption></fig> <fig position="float" id="F3"><label>Figure 3</label><caption><p><bold>The tissue distribution of porcine <italic>PNAS-4 </italic>mRNA assessed by qPCR</bold>. The values shown in this figure are the averages of three independent experiments. Error bars represent the SD (<italic>n </italic>= 3) of relative mRNA expression levels of <italic>PNAS-4 </italic>normalized to <italic>ACTB</italic>. The values were normalized to endogenous <italic>ACTB </italic>expression and the value of <italic>PNAS-4 </italic>in brain was arbitrarily set to 1.0. In order to tell the detailed difference between the varied tissues, large intestine was used as a control. (* represents <italic>P </italic>&lt; 0.05, ** represents <italic>P </italic>&lt; 0.01).</p></caption></fig> <fig position="float" id="F4"><label>Figure 4</label><caption><p><bold>The developmental profile of porcine <italic>PNAS-4 </italic>mRNA expression</bold>. Porcine <italic>PNAS-4 </italic>mRNA expression levels in skeletal muscle normalized against <italic>ACTB </italic>during embryonic development and the postnatal period, measured by qPCR. Error bars represent SD (<italic>n </italic>= 3). Results were averaged from three independent assays during all stages observed. Symbols used are: *<italic>P </italic>&lt; 0.05; **<italic>P </italic>&lt; 0.01. E33: embryonic day 33; E65: embryonic day 65; E90: embryonic day 90; N2: muscle from 2-day neonate; N28: muscle from 28-day neonate.</p></caption></fig> <fig position="float" id="F5"><label>Figure 5</label><caption><p><bold>Subcellular distribution of CGI-146 protein tagged with GFP</bold>. The CGI-146-<italic>GFP </italic>fusion protein localized in Golgi complex marked with fluorescein isothiocyanate (A). Red fluorescence of mitochondria were stained with the mitochondrial-specific dye, Mito-Tracker Red (B) and blue fluorescence of nuclei were stained with Hoechst33342 (C). The overlay images were produced by merging all three signals together (D). The control micrographs (E-H) showed the GFP in the cytoplasm. Bars represent 30 <italic>μ</italic>M.</p></caption></fig>
<table-wrap position="float" id="T1"><label>Table 1</label><caption><p>Genotypes of porcine SNPs T1400C and A1813C in different pig breeds</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="left">Pig Breeds</td><td align="left">No.</td><td align="left" colspan="3">T1400C (<italic>Sty </italic>I)</td><td align="left" colspan="3">A1813C (<italic>Msp </italic>I)</td></tr><tr><td/><td/><td colspan="6"><hr/></td></tr><tr><td/><td/><td align="left"><italic>TT</italic></td><td align="left"><italic>TC</italic></td><td align="left"><italic>CC</italic></td><td align="left"><italic>AA</italic></td><td align="left"><italic>AC</italic></td><td align="left"><italic>CC</italic></td></tr></thead><tbody><tr><td align="left">Wuzhishan</td><td align="left">43</td><td align="left">0</td><td align="left">0</td><td align="left">43</td><td align="left">7</td><td align="left">18</td><td align="left">18</td></tr><tr><td align="left">Xiang</td><td align="left">42</td><td align="left">0</td><td align="left">21</td><td align="left">21</td><td align="left">8</td><td align="left">13</td><td align="left">21</td></tr><tr><td align="left">Bamaxiang</td><td align="left">45</td><td align="left">0</td><td align="left">14</td><td align="left">31</td><td align="left">0</td><td align="left">0</td><td align="left">45</td></tr><tr><td align="left">Tongcheng</td><td align="left">44</td><td align="left">1</td><td align="left">11</td><td align="left">32</td><td align="left">5</td><td align="left">11</td><td align="left">28</td></tr><tr><td align="left">Landrace</td><td align="left">15</td><td align="left">12</td><td align="left">0</td><td align="left">3</td><td align="left">12</td><td align="left">0</td><td align="left">3</td></tr><tr><td align="left">Large White <sup>a</sup></td><td align="left">19</td><td align="left">15</td><td align="left">0</td><td align="left">4</td><td align="left">15</td><td align="left">0</td><td align="left">4</td></tr><tr><td align="left">Large White <sup>b</sup></td><td align="left">196</td><td align="left">139</td><td align="left">54</td><td align="left">3</td><td align="left">139</td><td align="left">54</td><td align="left">3</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T2"><label>Table 2</label><caption><p>Association analyses between SNP A1813C and three fat-related traits</p></caption><table frame="hsides" rules="groups"><thead><tr><td/><td align="left" colspan="3">SNP A1813C Genotype (mean ± SE)</td><td align="left" colspan="3">P-value</td></tr><tr><td/><td colspan="6"><hr/></td></tr><tr><td align="left">Traits</td><td align="left">AA (67)</td><td align="left">AC (52)</td><td align="left">CC (37)</td><td align="left">AA-AC</td><td align="left">AA-CC</td><td align="left">AC-CC</td></tr></thead><tbody><tr><td align="left">Percentage of leaf fat</td><td align="left">2.794 ± 1.31</td><td align="left">3.259 ± 1.2</td><td align="left">5.261 ± 1.44</td><td align="left">0.847</td><td align="left">0.008**</td><td align="left">0.017*</td></tr><tr><td align="left">Percentage of leaf and caul fat</td><td align="left">5.010 ± 2.29</td><td align="left">5.928 ± 2.6</td><td align="left">9.933 ± 2.96</td><td align="left">0.806</td><td align="left">0.012*</td><td align="left">0.028*</td></tr><tr><td align="left">BFT between 6<sup>th </sup>and 7<sup>th </sup>ribs</td><td align="left">3.257 ± 0.93</td><td align="left">3.199 ± 0.71</td><td align="left">4.262 ± 0.93</td><td align="left">0.183</td><td align="left">0.324</td><td align="left">0.039*</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T3"><label>Table 3</label><caption><p>Primers and probe employed in these experiments</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="left">Function</td><td align="left">Primer name</td><td align="left">Primer and probe sequences (5'-3')</td><td align="left">Tm</td></tr></thead><tbody><tr><td align="left">Cloning</td><td align="left">GAPF</td><td align="left">CTCCAAAGTCACACGCTCAGAAC</td><td align="left">62</td></tr><tr><td/><td align="left">GAPR</td><td align="left">ATAGTCCTACCCCACGAAGAAGC</td><td/></tr><tr><td/><td align="left">5GSP</td><td align="left">GGAAGGGTATGCAGGAGCTGAAGTAGGC</td><td align="left">71</td></tr><tr><td align="left">Subcellular</td><td align="left">P4GFPF</td><td align="left">CGGCTCGAGGGATGGGGGCTAACC</td><td align="left">63</td></tr><tr><td align="left">localization</td><td align="left">P4GFPR</td><td align="left">GAGAAGCTTAGAGTTTGGTGTGGCG</td><td/></tr><tr><td align="left">Real-time</td><td align="left">Real-time F</td><td align="left">CCCAGGAAATGCTTCTGAACTAG</td><td align="left">60</td></tr><tr><td align="left">PCR</td><td align="left">Real-time R</td><td align="left">CCCACAAAGAATCTCTGATAAAGC</td><td/></tr><tr><td/><td align="left">ProbeP4</td><td align="left">AAGCTGTTGTTTTGGGGAGCACTG</td><td/></tr><tr><td/><td align="left"><italic>ACTB </italic>F</td><td align="left">GGATGCAGAAGGAGATCACG</td><td align="left">60</td></tr><tr><td/><td align="left"><italic>ACTB </italic>R</td><td align="left">CTCGTCGTACTCCTGCTTGC</td><td/></tr><tr><td/><td align="left">ProbeBeta</td><td align="left">ATCAAGATCATCGCGCCTCCCGAGC</td><td/></tr><tr><td align="left">Mapping</td><td align="left">MapF</td><td align="left">TGGCAGAGCGGTCGTCACTTAGGC</td><td align="left">65</td></tr><tr><td/><td align="left">MapR</td><td align="left">GAGCAGAATCTCCTCCCACGCACCAT</td><td/></tr><tr><td align="left">SNPs detection</td><td align="left">SNPF1</td><td align="left">CTAGAACCACTCAAACCAAGCAGC</td><td align="left">62</td></tr><tr><td/><td align="left">SNPR1</td><td align="left">ATCAGGCAGGTAAAAGGATAACGG</td><td/></tr><tr><td/><td align="left">SNPF2</td><td align="left">GCCTTCTGAGTAGCAGTATGAGTTG</td><td align="left">62</td></tr><tr><td/><td align="left">SNPR2</td><td align="left">CCTGCGAGAACTGAGAATAATCC</td><td/></tr></tbody></table></table-wrap>
<disp-formula><italic>Y</italic><sub><italic>ijkl </italic></sub>= <italic>μ </italic>+ <italic>C</italic><sub><italic>i </italic></sub>+ <italic>B</italic><sub><italic>j </italic></sub>+ <italic>S</italic><sub><italic>k </italic></sub>+ (<italic>CB</italic>)<sub><italic>ij </italic></sub>+(<italic>BS</italic>)<sub><italic>jk </italic></sub>+ <italic>e</italic><sub><italic>ijkl</italic></sub></disp-formula>
<table-wrap-foot><p><sup>a </sup>DNA samples collected from our experimental pig group.</p><p><sup>b </sup>DNA samples collected from two experimental Large White selection lines originating from one commercial population.</p></table-wrap-foot> <table-wrap-foot><p>* represents <italic>P </italic>&lt; 0.05, ** represents <italic>P </italic>&lt; 0.01</p></table-wrap-foot>
<graphic xlink:href="1471-2156-9-40-1"/> <graphic xlink:href="1471-2156-9-40-2"/> <graphic xlink:href="1471-2156-9-40-3"/> <graphic xlink:href="1471-2156-9-40-4"/> <graphic xlink:href="1471-2156-9-40-5"/>
[{"surname": ["Rehfeldt", "Stickland", "Fiedler", "Wegner"], "given-names": ["C", "NC", "I", "J"], "article-title": ["Environmental and Genetic Factors as Sources of Variation in Skeletal Muscle Fibre Number"], "source": ["Basic Appl Myol"], "year": ["1999"], "volume": ["9"], "fpage": ["235"], "lpage": ["253"]}, {"surname": ["Te Pas", "Pool", "Hulsegge", "Janss"], "given-names": ["MFW", "MH", "I", "LLG"], "article-title": ["Analysis of the differential transcriptome expression profiles during prenatal muscle tissue development."], "source": ["Arch Tierz Dummerstorf"], "year": ["2006"], "volume": ["49"], "fpage": ["110"], "lpage": ["115"]}, {"surname": ["Zhu"], "given-names": ["ZM"], "article-title": ["Identification of pork production trait candidate genes by cDNA macroarray"], "source": ["Department of Animal Science"], "year": ["2004"], "publisher-name": ["Wuhan , Huazhong Agricultural University"]}, {"surname": ["Filippov", "Filippova", "Sinha", "Duerksen-Hughes"], "given-names": ["V", "M", "D", "PJ"], "article-title": ["PNAS-4: A Novel Pro-Apoptotic Gene Activated During the Early Response to DNA Damage [abstract]"], "year": ["2005"]}, {"surname": ["Daniel"], "given-names": ["SB"], "article-title": ["Modelling the p53 gene regulatory network"], "source": ["Centre for Mathematics and Physics in the Life Sciences and Experimental Biology"], "year": ["2006"], "publisher-name": ["London , University of London"]}, {"surname": ["Te Pas", "Cagnazzo", "de Wit", "Priem", "Pool", "Davoli"], "given-names": ["MFW", "M", "AAC", "J", "M", "R"], "article-title": ["Muscle transcriptomes of Duroc and Pietrain pig breeds during prenatal formation of skeletal muscle tissue using microarray technology. "], "source": ["Arch Anim Breed"], "year": ["2005"], "volume": ["48"], "fpage": ["141"], "lpage": ["147"]}, {"surname": ["Swatland"], "given-names": ["HJ"], "article-title": ["Structure and Development of Meat Animals and Poultry"], "source": ["Technomic Publishing Co Inc Lancaster, PA"], "year": ["1994"]}, {"surname": ["Dragos-Wendrich", "Moser", "Bartenschlager", "Reiner", "Geldermann"], "given-names": ["M", "G", "H", "G", "H"], "article-title": ["Linkage and QTL mapping for Sus scrofa chromosome 10"], "source": ["J Anim Breed Genet"], "year": ["2003"], "volume": ["120"], "fpage": ["82"], "lpage": ["88"], "pub-id": ["10.1046/j.0931-2668.2003.00427.x"]}, {"article-title": ["The PROSITE database"]}, {"article-title": ["SignalP 3.0 server"]}, {"article-title": ["UTRscan"]}, {"surname": ["Johnson"], "given-names": ["DE"], "article-title": ["Applied multivariate methods for data analysts"], "source": ["Stamford: Thomson Learning"], "year": ["1998"], "fpage": ["286"], "lpage": ["316"]}]
BMC Genet. 2008 Jun 30; 9:40