First survey and functional annotation of prohormone and convertase genes in the pig
BMC Genomics volume 13, Article number: 582 (2012)
The pig is a biomedical model to study human and livestock traits. Many of these traits are controlled by neuropeptides that result from the cleavage of prohormones by prohormone convertases. Only 45 prohormones have been confirmed in the pig. Sequence homology can be ineffective to annotate prohormone genes in sequenced species like the pig due to the multifactorial nature of the prohormone processing. The goal of this study is to undertake the first complete survey of prohormone and prohormone convertases genes in the pig genome. These genes were functionally annotated based on 35 gene expression microarray experiments. The cleavage sites of prohormone sequences into potentially active neuropeptides were predicted.
We identified 95 unique prohormone genes, 2 alternative calcitonin-related sequences, 8 prohormone convertases and 1 cleavage facilitator in the pig genome 10.2 assembly and trace archives. Of these, 11 pig prohormone genes have not been reported in the UniProt, UniGene or Gene databases. These genes are intermedin, cortistatin, insulin-like 5, orexigenic neuropeptide QRFP, prokineticin 2, prolactin-releasing peptide, parathyroid hormone 2, urocortin, urocortin 2, urocortin 3, and urotensin 2-related peptide. In addition, a novel neuropeptide S was identified in the pig genome correcting the previously reported pig sequence that is identical to the rabbit sequence. Most differentially expressed prohormone genes were under-expressed in pigs experiencing immune challenge relative to the un-challenged controls, in non-pregnant relative to pregnant sows, in old relative to young embryos, and in non-neural relative to neural tissues. The cleavage prediction based on human sequences had the best performance with a correct classification rate of cleaved and non-cleaved sites of 92% suggesting that the processing of prohormones in pigs is similar to humans. The cleavage prediction models did not find conclusive evidence supporting the production of the bioactive neuropeptides urocortin 2, urocortin 3, torsin family 2 member A, tachykinin 4, islet amyloid polypeptide, and calcitonin receptor-stimulating peptide 2 in the pig.
The present genomic and functional characterization supports the use of the pig as an effective animal model to gain a deeper understanding of prohormones, prohormone convertases and neuropeptides in biomedical and agricultural research.
In addition to the importance in livestock production, the pig is a well-established biomedical model to study human health due to the genomic, anatomic and physiologic similarities with humans. A wide variety of human health traits including cancer, reproductive health, drug metabolism, wound healing, and cardiovascular disease have been successfully studied using the pig[1–4]. Underlying these and other important traits such as growth and development, feeding, stress, memory and susceptibility to substances abuse are neuropeptides, a class of cell-cell signaling peptides that have paracrine, endocrine, and autocrine effects[5, 6]. Research in selected pig neuropeptides has offered insights into cell transplantation, nervous system diseases, and injury. For example, Yang et al. reported on the inhibitory effect of neuromedin S (NMS) on luteinizing hormone secretion which is mediated via melanocyte-stimulating hormone (MSH) neurons in the hypothalamus-pituitary axis of ovariectomized pigs. Kaminski et al. concluded that hypothalamic peptides, orexin A and orexin B, were involved in the control of food intake, sleep patterns, autonomic and neuroendocrine systems, and may also affect reproductive functions through the influence on the hypothalamic-pituitary-ovarian axis in pigs.
The identification of neuropeptides is more difficult than most proteins because neuropeptides are derived from larger prohormone proteins as a result of complex enzymatic processing. The conversion of the large prohormone to one or multiple smaller neuropeptides involves cleavage by multiple prohormone convertases and additional post-translational modifications such as amidation and glycosylation. This complex processing of prohormones into neuropeptides challenges the identification of neuropeptide genes across genomes solely based on sequence homology to better understood species[5, 6, 10].
Only 40 prohormone and 2 prohormone convertase genes have been empirically confirmed in the pig compared to approximately 100 genes identified in human, rat, mouse, cow and rhesus monkey[11–13]. This partial knowledge of the prohormone, prohormone convertase and associated neuropeptides in the pig is a critical shortcoming, especially considering the important role of pig in biomedical and agricultural research. In addition, few gene expression studies have discussed the expression profile of pig prohormone or prohormone convertase genes. Ross et al. found that estradiol treatment at day 9 of gestation was associated with changes in the expression of the prohormone neuromedin (NMB) in the endometrium of glits. Hausman et al. concluded that the expression of neuropeptide Y (NPY) was down-regulated with age in gilts ranging from 90 to 210 days old.
Understanding the role of neuropeptides in human and livestock traits using the pig as biomedical model requires a comprehensive knowledge of the neuropeptide complement in the recently released pig genome (SScrofa 10.2v18,). This understanding includes the identification of prohormone and prohormone convertase genes, prediction of cleavage sites in prohormones that may result in potentially bioactive neuropeptides, and characterization of gene expression and protein abundance across conditions to gain insights into the role of neuropeptides. A complete survey of the prohormone and prohormone convertase genes in the pig supports the interpretation of gene expression experiments and improves the effectiveness of tandem mass spectrometry studies to identify neuropeptides[17–19]. Functional annotation of these genes can be obtained by the analysis of the large number of gene expression experiments already available[20, 21]. To address the lack of a comprehensive understanding of the prohormone and prohormone convertase genes in the pig, we present the first comprehensive survey and functional annotation of these genes. An all-inclusive catalogue of prohormone and prohormone convertase genes known in other species was used to search complementary pig genome databases. These genes were then characterized by analyzing a large number of gene expression experiments across a wide range of conditions. The potential cleavage sites of prohormones that can result in bioactive neuropeptides were predicted and compared to the cleavages based on known neuropeptide sequences.
Results and discussion
Pig prohormone genes
A comprehensive catalogue of 95 potential pig prohormone genes, 8 prohormone convertase genes and 1 prohormone convertase facilitator gene (7B2) were identified in the pig genome. Table1 lists the genes and the corresponding BLAST matches on the pig Genome, UniProt, Gene and UniGene databases. There are 11 previously unreported (i.e. without empirical confirmation) prohormone genes in the pig and complete sequences where only partially or incomplete sequences have been previously reported. Newly identified genes are intermedin (ADM2), cortistatin (CORT), insulin-like 5 (INSL5), orexigenic neuropeptide QRFP (OX26), prokineticin 2 (PROK2), prolactin-releasing peptide (PRRP), parathyroid hormone 2 (TIP39), urocortin (UCN1), urocortin 2 (UCN2), urocortin 3 (UCN3), and urotensin 2-related peptide (UTS2B). This search also identified two different calcitonin protein entries in public databases that are isoforms of other calcitonin genes. Additional information on the comprehensive catalogue of genes is available [see Additional file1]. The predicted prohormone protein sequences with cleavage sites identification is provided in the NeuroPred format [see Additional file2].
Table2 summarizes the evidence from various repositories supporting the 95 unique pig prohormone genes and 2 alternative sequences detected in this study. The genome-predicted prohormone genes are supported by 66 UniProt entries (Table1) including 47 sequences that have evidence at the protein level of which 39 have complete protein sequence and 8 have partial protein sequence. Additionally, 81, 91 and 19 prohormone genes detected are supported by transcript evidence from UniGene, Gene, and UniProt, respectively. Two genes, apelin (APEL) and PROK2, are supported by evidence in UniProt and by EST records unassigned to a UniGene cluster. UniProt supporting information includes 2 entries from alternatively spliced calcitonin genes, preprocalcitonin gene-related peptide and calcitonin-2, that have evidence at the protein and transcript levels, respectively. However, the calcitonin 2 is not assigned to a UniGene cluster or NCBI Gene identifier because no EST matches the unique region of the reported sequence. The UniProt entry for prothyroliberin (TRH) refers only to the 3-amino acid thyroliberin peptide however, no pig EST has been reported. Only 8 genome predictions do not have supporting evidence in UniProt with 5 of these have supporting transcript evidence. Complete protein sequences were predicted for insulin-like 5 (INSL5) and endothelin-2 (EDN2) and an incomplete protein sequence was predicted for prolactin-releasing peptide (PRRP) genes. The nucleic and protein sequence of EDN2 have been reported but this information is not present in public databases. Of the 23 UniProt predicted sequences, 14 (7) entries have (do not have) transcript support. There are 2 UniProt predicted sequences that correspond to genes that are absent in the pig genome studied. A partial match to UCN1 was found in the trace archives although this gene was not found in the pig genome assembly studied and lacks of external validation. Progonadoliberin-2 (GON2) was not found in the current pig genome assembly although this gene was detected in earlier assembly versions and trace archives. This gene has no current experimental evidence outside homology to other mammalian species. The apparent lack of UCN1 and GON2 in the assembly and fragment evidence of other prohormones is most likely due to poor coverage of the genomic regions where these prohormone genes are located.
At least four calcitonin genes, also known as calcitonin receptor-stimulating peptide genes, were identified with 2 genes exhibiting alternative splicing. The genome assembly permitted the assignment of the UniProt pre-procalcitonin (A6P7L6) and preprocalcitonin gene-related peptide (A6P7L7) entries to the same gene that also produces the UniProt calcitonin (CALC, [Swiss-Prot:P01259]) and calcitonin gene-related peptide (CALCA, [Swiss-Prot:P30880]) peptides, respectively. This alternatively spliced gene corresponds to CALCA gene found in other mammalian species.
The bioinformatics gene prediction pipeline confirmed that the separate Q766Y6 and A0A761 UniProt entries are alternatively spliced variants of the same calcitonin receptor-stimulating peptide 3 (CRSP3) gene as initially reported by Rezaeian et al.. A single genome match was identified for calcitonin receptor-stimulating peptide 2 (CRSP2, Q766Y7). While the calcitonin receptor-stimulating peptide 1 (CRSP1) gene has been reported, the actual gene may be inaccurately assembled in the 10.2 genome release because the region appears to contain a small duplication leading to 2 starting locations. Further supporting this argument, a conserved 19 amino acid region in all calcitonin receptor-stimulating peptide-related protein sequences matched to an additional 5th genome site which was part of a discontinued NCBI Gene entry (Gene ID 100624618). There is insufficient information to conclude whether there is a separate coding gene involved or an assembly-related problem.
Neuropeptide S (NPS) is a potential 12th prohormone gene discovered by the bioinformatics gene prediction pipeline. Two genome matches on different chromosomes for the NPS gene were identified. However, the complete identity of the matched between the predicted sequence and chromosomal locations implied that this was an assembly error rather than a duplication event. Although UniProt has an partial pig NPS entry (B5M997), the genome predicted protein sequence was more similar to the bovine sequence, as expected, than the present partial UniProt pig sequence. The pig NPS protein and nucleic sequences were found to be 100% and 99% identical, respectively, to the rabbit sequence ([GenBank:EU978456]). The similarity between the UniProt pig and rabbit sequences was also evident in the phylogenetic relationships among NPS sequences reported by Yao et al.. These findings call into question the present pig NPS entry in UniProt.
Prohormone genes previously unreported in pig
The genome search identified 11 prohormone genes that do not have empirical confirmation in the UniProt, UniGene or Gene databases (Table1). These genes are ADM2, CORT, INSL5, OX26, PROK2, PRRP, NPS, TIP39, UCN1, UCN2, UCN3, and UTS2B. Only inferred sequences are available for ADM2 and CORT in UniProt and the current pig NPS entry is identical to the rabbit NPS. There is evidence for mammalian homologs of all these genes in UniProt. The protocol followed to identify these genes included a high percentage of identities and similarities with a minimum percentage of mismatches and gaps and conservation of the region encompassing the potential neuropeptide. Intermedin or adrenomedullin 2 (ADM2) is part of the calcitonin family of peptides and has effects similar to those of adrenomedullin (ADML). In humans, intermedin causes hypotension when given peripherally and augments blood pressure and causes sympathetic activation when given to the central nervous system. This neuropeptide induces prolactin release, has anti-diuretic and natriuretic properties and reduces food intake. The amino acid sequences of CORT and somatostatin (SMS) are highly similar and both reduce neuronal activity. In addition, CORT has unique roles such as induction of slow-wave sleep, reduction of locomotor activity, and activation of cation selective currents not responsive to SMS. Although the function of INSL5 is still being determined, high expression in the colon, as well as in the brain and hypothalamus, indicates roles in gut contractility and neuroendocrine signaling. Likewise, the function of OX26 is still being elucidated, although studies in chicken confirm the orexigenic, appetite stimulating activity of this neuropeptide. Takayanagi and Onaka demonstrated that PRRP plays a role in control of energy metabolism and stress response. Prokineticins are involved in tumorigenesis process (prostate, testicles, neuroblastoma, colon, and pancreas) acting as a growth factor for cancer cells, an angiogenic and a chemotactic factor for pro-inflammatory neutrophils. NPS has anxiolytic-like effects (stress reduction) and can induce arousal and wakefulness. TIP39 and the corresponding receptor form a neuromodulator system and the anatomical distribution indicates a role in limbic, endocrine, viscerosensory, and auditory functions. This system has been postulated as potential drug target in anxiety, depression and chronic pain management. Urocortins and their receptors has been found in the central nervous, digestive, reproductive, cardiovascular, immune and endocrine systems, suggesting a variety of roles including cardiovascular activity and cell survival. UTSB2 is a paralog of urotensin 2 (UTS2) that exerts similar biological effects including relaxation of muscles and reduction of blood pressure.
Pig prohormone convertase genes
The sequence of 8 prohormone convertase genes and the 7B2 facilitator gene also known as secretogranin 5 (SCG5) were identified in the pig genome (Table1). The UniProt and Gene databases only had supporting evidence for PCSK1, PCSK2, and 7B2. Six additional prohormone convertase genes (furin, PCSK4, PCSK5, PCSK6, PCSK7, and PCSK9) that were previously unreported or not based on empirical evidence were identified. Only transcript evidence supports the 8 prohormone convertase genes, meanwhile protein evidence is available for 7B2. Dai et al. isolated PCSK1 from the ovary cDNA library of a pregnant sow and Renegar et al. detected PCSK1 in the corpus luteum and brain of pregnant sows. Also, mRNA from PCSK1 and PCSK2 has been identified in the pituitary neurointermediate lobes of pigs. Among the prohormone convertases, furin, PCSK4, PCSK5, PCKS6 and PCSK7 do not have UniGene entries. The present catalogue enhances the currently limited work on pig prohormone convertases.
Functional characterization of the pig prohormone and prohormone convertase genes
Analysis of the large number of microarray gene expression experiments enabled the first comprehensive characterization of the role of prohormone and prohormone convertase genes in biological processes in the pig. The results from these analyses augmented the understanding of the role of these genes on reproduction, health, growth, and other traits of importance to biomedical research and agricultural production.
The query of Affymetrix Porcine Genome Array identified 77 probes representing 56 prohormone and 3 prohormone convertase genes. Table3 lists the total number of differentially expressed probes (P-value < 0.005) within the seven experimental classes considered. A detailed distribution of the differential expression level of each probe and experiment is provided [see Additional file3. A discussion of the findings for the 3 groups with highest number of differentially expressed probes (immune-related, embryo and placenta, and brain and central nervous system) is presented below. Although neuropeptides expressed in the brain and the immune system interact with circulating cytokines to support two-way communications between the brain and immune system, we describe the profiles of prohormones in immune-related tissues separately from the brain and central nervous system tissues to facilitate the interpretation of results.
Several studies have demonstrated that prohormone genes play an important role in pig immune response. This was evidenced by the high number of differentially expressed prohormone and prohormone convertase genes (24 genes) among experiments that evaluated immune-response in blood, spleen, and lymph nodes (Table3). Differentially expressed genes were: ADML, augurin (AUGN), cholecystokinin (CCKN), CRSP3, endothelin-1 (EDN1), galanin (GALA), galanin-like peptide (GALP), progonadoliberin-1 (GON1), insulin-like growth factor I (IGF1), insulin-like growth factor II (IGF2), neuromedin-B (NMB), neuromedin-U (NMU), neuropeptide Y (NPY), platelet-derived growth factor subunit A (PDGFA), proenkephalin-A (PENK), prorelaxin 1 (REL1), secretogranin-1 (SCG1), secretogranin-2 (SCG2), secretogranin-3 (SCG3), SMS, vascular endothelial growth factor C (VEGFC), vascular endothelial growth factor D (VEGFD), PCSK1 and PCSK7.
In general, prohormone genes were under-expressed in pigs under immune challenge relative to the un-challenged controls. AUGN was differentially expressed in two experiments; GSE7313 that profiled lymph nodes and GSE14790 that profiled blood. In GSE14790, 7 day-old pigs were inoculated with porcine circovirus type 2 (PCV2), a virus that is widely spread across pig farms, and gene expression was profiled at 0, 7, 14, 21 and 29 dpi. AUGN was over-expressed in un-inoculated pigs at 29 dpi relative to 7 dpi, regardless of inoculation and relative 21 dpi inoculated pigs (P-value < 2.5 × 10-4). Both contrasts indicate that the expression of AUGN increases with age and this trend is slower in pigs infected with PCV2. In GSE7313, the gene expression of seven week old piglets inoculated with Salmonella Typhimurium was profiled at 8 hours post inoculation (hpi), 24 hpi, 48 hpi, and 21 days post inoculation (dpi). AUGN was over-expressed at 21 dpi relative to 24 hpi and 48 hpi (P-values < 6.8 × 10-5 and 2.7 × 10-6, respectively). Consistent with the differential expression in relation to immune-response observed in this study, AUGN is a putative tumor suppressor gene and is down-regulated in many cancers.
IGF2, a member of the insulin family and is involved in development and growth, was differentially expressed across immune-related experiments. IGF2 was represented by 12 probes in the microarray platform and 6 probes were differentially expressed across experiments. Five probes (Ssc.9365.1.S1_at, Ssc.9365.2.S1_a_at, Ssc.9365.5.A1_at, Ssc.9365.5.S1_at, Ssc.9365.6.S1_x_at) and one probe (Ssc.9365.3.S1_a_at) were differentially expressed in experiments GSE14790 and GSE7314, respectively. In GSE14790, IGF2 was under-expressed in non-inoculated piglets at 7 dpi relative to inoculated pigs at various days (P-value < 1.7 × 10-5, fold change = 0.71). In experiment GSE7314, IGF2 was over-expressed in pigs inoculated with Salmonella choleraesuis at 21 dpi relative to non-inoculated pigs (P-value < 8.1 × 10-4). These results are consistent with reports that IGF2 is down-regulated in pigs immune-challenged with lipopolysaccharide.
SCG1, SCG2, SCG3, members of the secretogranin family, exhibited differential expression among immune-challenge experiments consistent with the known association of these genes with cell activation, cytotoxicity and microbial defense. Probes on all 3 SCGs exhibited differential expression on two immune-related experiments. SCG1 and SCG2 are differentially expressed in GSE14790 while SCG3 was differentially expressed in GSE11787. In GSE11787 SCG3 was under-expressed in inoculated pigs relative to controls (P-value < 1.2 × 10-3, fold change = 0.33). These results are consistent with the lack of synthesis of endogenous granins in rat PC12 cells infected with recombinant vaccinia viruses. In GSE14790, SCG1 and SCG2 were under-expressed in pigs inoculated with PCV2 relative to un-inoculated control pigs (P-value < 1.8 × 10-4, fold change = 0.88 and P-value < 1.5 × 10-5, fold change = 0.93, respectively).
Two members of the vascular endothelial growth factor family, VEGFC and VEGFD, were under-expressed in PCV2 inoculated pigs relative to control pigs (P-value < 1.8 × 10-5, fold change = 0.62) and also under-expressed at early stages (7 dpi) relative to later stages (19 and 29 dpi) in GSE14790. In agreement with these findings, a loss of endothelial growth factor transcription and increase in pro-inflammatory indicators were reported in the endometrial lymphocytes of pigs at sites of fetal arrest.
NPY (probe Ssc.15981.1.A1_at) was under-expressed in PCV2-inoculated pigs relative to control pigs (P-value < 6.6 × 10-4) and, within infection level, NPY was under-expressed at earlier stages relative to 29 dpi in GSE14790. Consistent with these findings, the levels of NPY mRNA decreased in the blood of rats treated with vinblastine, an anti-cancer drug known to decrease the number of white blood cells of the immune system involved in defense. Similarly, NPY was found to decrease in cattle infected with Bovine Spongiform Encephalopathy.
ADML was differentially expressed in GSE14758-D and GSE7314. In GSE14758-D, ADML was under-expressed in the mediastinal lymph nodes of PCV2-infected pigs relative to control pigs at 29 dpi (P-value < 1.7 × 10-3, fold change = 0.6). Whereas, in GSE7314 ADML was over-expressed at 48 dpi in the blood of pigs inoculated with Salmonella choleraesuis relative to controls (P-value 4.8 × 10-3). The latter result is consistent with the up-regulation of ADML gene expression and increases in systemic circulatory concentrations of ADML in response to the onset and progression of trauma, infection, and sepsis. The apparent inconsistency between both experiments may be associated with the differential effects that ADML has on cellular metabolism, immune function, endocrine function, and cardiovascular function.
Of the 3 prohormone convertases available in the microarray platform, PCSK1 and PCSK7 were significantly differentially expressed (P-value < 1.3 × 10-3) and PCSK2 was marginally significantly differentially expressed (P-value < 6.5 × 10-3) in GSE14790. PCSK1 was under-expressed in PCV2-inoculated pigs already at 7 dpi relative to 29 dpi, regardless of inoculation at the later stage (P-value < 5.7 × 10-5). Likewise, PCSK7 is under-expressed in PCV2-inoculated pigs relative to controls already at 7 dpi (P-value < 4.2 × 10-4) and, within controls, PCSK7 was under-expressed at early stages (7 dpi, 21 dpi) relative to 29 dpi (P-value < 4.1 × 10-4). These results are in agreement with similar findings that furin, another prohormone convertase, was dysregulated in the immune cells of advanced human atherosclerotic plaques and imply that prohormone convertase, like prohormone genes, are down regulated under immune challenges.
Embryo and placenta profiling
In GSE18641, IGF2 (probe Ssc.9365.2.S1_a_at) was over-expressed in pregnant sows relative to non-pregnant sows (P-value < 2.7 × 10–3, fold change 1.23). In GSE12705, IGF2 (probes Ssc.9365.4.S1_a_at and Ssc.9365.5.S1_at) was over-expressed in earlier stages (day 11 spherical and day 11 and 12 tubular) relative to later stage (day 12 and 14 filamentous) conceptuses (P-value < 2.1 × 10-4). This IGF2 profile is supported by Pantaleon et al. that showed that IGF2 is needed in order for mouse embryos to progress from early stages to blastocyst stages. Gupta et al.[54, 55] reported that the expression of the embryo survival related gene IGF2 increased with the addition of nonessential amino acids or phytohemaglutinin in pig embryos and blastocysts, respectively.
Both PENK probes were over-expressed in tubular and spherical conceptuses relative to filamentous conceptuses (P-value < 2.6 × 10-6) in experiment GSE12705. This is consistent with results that found PENK mRNA to increase linearly during gestation in the hippocampus of pigs. PTHR was under-expressed in tubular and spherical relative to filamentous conceptuses (P-value < 6.4 × 10-7, fold change = 0.02) in experiment GSE12705. This finding is supported by reports that PTHR is present in higher concentrations in fetal pigs than in sows. VEGFC is a representative of the vascular endothelial growth factor family of prohormones that have an important role in the survival and mitogenesis of endothelial cells and lymphaniogenesis and angiogenesis of embryos. VEGFC was over-expressed in pregnant sows relative to non-pregnant sows (P-value < 7.8 × 10-4) in experiment GSE18641. This finding is supported by a study in the chicken, demonstrating that the chorioallantoic membrane (analogous to the placenta in mammals) contained growth of embryonic microvessels stimulated by VEGFC. The expression profile is also supported by the finding that in mice embryos, VEGFC is required for successful lymphatic vasculature development and lymphatic endothelial cell migration.
Brain and central nervous system
Eleven differentially expressed prohormone genes were identified in experiments concerning the hypothalamus, thyroid, and olfactory bulb (neuroblasts). These genes are Adrenomedullin-5 (ADM5), ADML, C-type natriuretic peptide (ANFC), cocaine and amphetamine regulated transcript protein (CART), IGF1, IGF2, NPY, platelet-derived growth factor subunit A (PDGFA), prodynorphin (PDYN), PTHR, and VEGFC.
ADML was over-expressed in the immortalized porcine olfactory bulb neuroblasts relative to the non-neural epithelial cells (P-value < 2.2 × 10-6, fold change > 10) in experiment GSE16855. This result is supported by a previous study that found that ADML is important for regulation of proliferation and differentiation of neural stem/progenitor cells using the mouse olfactory bulb.
IGF1 was over-expressed in the neuroblasts relative to non-neural epithelial cells (average P-value < 5 × 10-7, fold change > 10) in experiment GSE16855. This result is supported by a study in chickens showing that IGF1 was expressed in the olfactory bulb. Also, IGF2 (probe Ssc.9365.6.S1_x_at) was consistently over-expressed in the hypothalamus of male Iberian pigs relative to all other seven breed-gender combinations (on average, P-value < 2.3 × 10-4, fold change = 2.42) in experiment GSE14739-H[63, 64]. NPY was over-expressed (P-value < 8.1 × 10–4, fold change = 7.94) in neuroblasts relative to non-neuronal cells in GSE16855. This result is consistent with reports that the olfactory bulb exhibit high levels of immunoreactive NPY in the brain of pigs and that NPY may inhibit excitatory neurotransmission in the rat olfactory bulb. VEGFC was over-expressed in neuroblasts relative to non-neuronal cells (P-value < 1.5 × 10-9, fold change > 10) in experiment GSE16855. This result agrees with a 30% increase in dividing neuroblasts in olfactory bulb in culture stimulated with VEGFC compared to controls reported by Le Bras et al.. PTHLH was under-expressed in neuroblasts relative to non-neuronal cells (P-value < 2.6 × 10-4, fold change = 0.20) in GSE16855. This finding is consistent with reports that PTHLH may be a negative regulator in the differentiation of chondrocytes. PDGFA was over-expressed in neuroblasts relative to non-neuronal cells (P-value < 1.2 × 10-4) in experiment GSE16855. Related to this result, Fressinaud et al. reported that platelet-derived growth factors increase the glutamine synthetase activity in astrocytes in the brain.
Prediction of cleavage sites in pig prohormones
All 97 prohormone sequences were used to predict cleavage and confirm the prediction against known or predicted cleavage sites. These sequences were inferred to have 228 cleavage sites that resulting in a 14.6% prevalence rate (proportion of possible sites that are cleaved). Most sites were cleaved at an arginine (R) such that the most frequently cleaved motifs were xxKR (71%), RxxR (34%) and xxRR (41%), where x denotes any amino acid and K denotes lysine. There were 5% (38) C-terminal single R sites that were cleaved without a basic amino acid in the second and fourth positions preceding the cleavage site (P2 or P4 locations, respectively).
The performance of the cleavage prediction models is presented in Table4. The correct classification rate ranged from 82% to 92% indicating that a large proportion of the sites were accurately predicted across all models. The human cleavage prediction models had the best performance for most of the statistics followed by the mammalian model. The Known Motif model provided the highest number of true positive predictions but also the highest number of false positive predictions. The Known Motif model provided the highest sensitivity, 77%, indicating more than three quarters of the cleaved sites were correctly predicted as cleaved. However this model also provided the highest number of false positive predictions. Consequently the Known Motif positive predictive power was 35% indicating that, on average, only 35% of sites predicted to be cleaved are expected to be true cleavage sites.
The human models provided the highest number of true negatives resulting in the best model performance compared to the Known Motif and mammalian models. The human artificial neural network models had approximately 60% positive predictive power indicating that most sites predicted as cleaved are expected to be true positives. Although the human logistic models had lower sensitivity than their artificial neural network counterparts, the differences with the human artificial neural network model were only 4 cleaved and 11 non-cleaved sites. The high performance of the human models suggests that the cleavage of prohormones that result in potential biologically active neuropeptides in the pig is similar to humans. Noteworthy is that the mammalian model was trained on 51 mammalian prohormones that included 8 pig prohormones. This model provided slightly more true positive predictions and a higher sensitivity than the human logistic model. However, the mammalian model had noticeably more false positive predictions than the human logistic model resulting in lower performance in the other accuracy measures.
The comparison of results across models also provides information on the accuracy of the cleavage assignment, prediction accuracy and potential for a gene to produce bioactive peptides. For 10 prohormones, at least 5 of the models did not predict any cleavage site. However, it must be noted that 4 of the prohormones (ANF, GHRL, IGF1 and PDGFD) are likely to have sites cleaved by proteases other than prohormone convertases. For example, ANF is cleaved by corin, serine peptidase (CORIN).
Genes with no predicted cleavage or assigned cleavage that differ from other species can be used to identify proteins are not cleaved to form smaller peptides. There is no evidence for cleavage of UCN2 and UCN3 to produce mature peptides in mammals. All models failed to predict two cleavage sites in TOR2X. The first site, an N-terminal dibasic ‘RK’, is known to be rarely cleaved across species. The second site is a cleavage found in humans that forms alpha- and beta-salusin but this site may not be cleaved in the pig since the pig sequence, like the bovine sequence, only has a single basic site instead of the human dibasic ‘RR’ site. Similarly for TKN4, the genomic prediction and supporting EST data indication a change from an R in other species to a glycine amino acid in the pig sequence that may prevent the formation of the ‘Hemokinin’ peptide.
The pig CRSP2 protein sequence lacks the ‘KR’ and a C-terminal cleavage site that are cleaved in human CALC and CALCB genes to produce Calcitonin gene-related peptide 1 and Calcitonin gene-related peptide 2. Therefore it is unlikely that pig CRSP2 would provide these calcitonin peptides. The assigned cleavages in the RES18 protein are necessary to provide a potential triskadecapeptide peptide reported by Bloomquist et al.. This potential peptide has flanking dibasic cleavage sites in the mouse and rat but this peptide has not been experimentally confirmed. The corresponding region in human, bovine and pig sequences are monobasic and lack common PC cleavage motifs suggesting that these species probably cannot form this peptide.
Examination of the potential cleavage sites in IAPP indicated that a mutation from R to Q in the N-terminal cleavage site is necessary to produce the Islet amyloid polypeptide. Examination of the corresponding ESTs indicated that 2 swine ESTs ([GenBank:AJ649149] and [GenBank:AJ649469]) were 100% identical to the rat genome and consequently invalid sequences. Two other ESTs ([GenBank:EW569366], [GenBank:BF712755]) matched the region that supported the genomic prediction. The predicted protein sequence including the potential cleavages sites of the expected IAPP was less than 80% identical to other mammals sequences compared to typically over 85% identity between the human and most other mammalian sequences. Potter et al. questioned the capability of IAPP to form amyloids after examining the functionality of a synthesized pig sequence based on the [GenBank:BF712755] EST sequence. The predicted prohormone sequence and cleavage prediction results also strongly suggest that the pig is unlikely to be able to form IAPP. This reflects the importance of proteomic studies involving cleavage to first determine that a species can produce a peptide.
The pig is an important biomedical and agricultural research species. Results from the first genome-wide study of pig prohormone and prohormone convertase genes, functional annotation and prediction of prohormone cleavage have been presented. This study was enabled by the availability of the pig genome sequence and of 35 gene expression experiments that evaluated a wide range of conditions in pigs. These results offer more insights into the role of neuropeptides on biological processes such as reproduction, development, growth, and health and support targeted empirical confirmation. The bioinformatics pipeline used in this study can be used to identify prohormones or other sets of genes in species with similar sequence resources. Confirmatory insight into the pig prohormones can be expected from proteomic mass spectrometry studies.
Combining complementary bioinformatic resources, 95 prohormone genes, 8 prohormone convertases and one cleavage facilitator were discovered in the pig genome and raw sequence repositories. We uncovered 11 prohormone genes that have not been previously reported and one potentially incorrectly reported. The high performance of the models used to predict cleavage in the pig prohormones suggests that the prohormone cleavage in pigs is similar to humans. The analysis of 35 gene expression experiments identified various neuropeptide genes differentially expressed in immune-related tissues, embryo and placenta and the central nervous system including AUGN, IGF2, the family of SCG s, NPY, ADM and ADML, PENK, PTHR, and VEGFC. Experiments are required to confirm that the pig does not produce the bioactive neuropeptides UCN2, UCN3, TOR2X, TKN4, IAPP, and CRSP2 as suggested by the cleavage prediction models.
Search for pig prohormone and convertase genes
A registry of approximately 100 candidate mammalian prohormone and convertase genes was built from public sequence databases (including NCBI Gene –release date September 2011, UniGene – release date April 13 2011, and UniProt–release 2011_11 November 16, 2011) and a literature review[11–13, 19, 78–80].
Candidate genes were searched for in the pig genome 10.2 assembly using the bioinformatics pipeline developed by Southey et al.[11, 13]. The protein sequence of each candidate gene in the registry was searched on the pig genome assembly using the TBLASTN, BLASTP and BLASTN programs from NCBI BLAST (version 2.18) with default parameters settings (E-value < 10 and BLOSUM62 scoring matrix) and filtering disabled. In addition, sequences not used in the pig genome assembly (including unassigned genomic regions, whole genome shotgun sequencing and trace archives) were searched when there was no suitable BLAST match to a candidate gene or when the alignment to the genome assembly suggested a missing genomic region. This strategy allowed the annotation of genomic regions that were partly or not included in the assembly.
The BLAST matches were examined based on the alignment score and E-value to identify the most likely matches and genomic location of the corresponding prohormone. The identified pig genomic region that encompassed the BLAST match was further extended approximately 500 base pairs to the 5′ and 3′ ends of the match. Matches were also screened for alignments to multiple homologous prohormone genes that could indicate gene duplication events in the pig genome. The gene parsing tool Wise2 was used to predict the protein sequence within the genome regions detected with BLAST. The genomic region was further extended when only a partial protein sequence was predicted. In this study, Wise2 compared the target protein (preference was given to pig protein sequences, followed by human, cattle and other mammals) to the pig genomic DNA sequence identified by BLAST to infer the gene structure based on a model that includes introns and frameshift errors. Each predicted gene was compared to the UniProt and NCBI Gene databases to assess the accuracy of the prediction based on previously reported pig genes. To further confirm the Wise2 predictions, the protein sequence predicted from the gene model was also compared to the corresponding published mammalian sequences using the multiple sequence alignment tool Clustalw. The multiple sequence alignment maximized the likelihood of identifying homologous genes. The predicted sequences were also searched against the pig entries in the NBCI EST database to confirm the presence of the predicted protein sequence. The pig entries in the NCBI EST database was also used to complete the protein sequence when the genome coverage was incomplete.
Functional annotation of the pig prohormone and convertase genes
A review of the pig microarray gene expression experiments available in the NCBI GEO database indicated that the Affymetrix Porcine Genome Array GPL3533 was the most commonly used platform. The UniGene database was searched for sequences that represent prohormone and prohormone convertase genes. This information was used to identify the probes representing prohormone and prohormone convertase genes in the Affymetrix Porcine Genome platform.
Thirty-five experiments that used the Affymetrix Porcine Genome platform were identified in GEO. Selected experiments had a minimum of 6 microarrays and a maximum of 80 microarrays. The sources and main features of these experiments are provided [see Additional file4]. The wide range of selected microarray experiments available supported a comprehensive characterization of the association of prohormone and associated neuropeptide and convertase genes with various biological processes.
The experiments were grouped into 7 classes: primary immune-response tissues, embryo and placenta, brain and central nervous system, reproduction, muscle, fat, and gut. For experiments encompassing multiple tissues (GSE14739, GSE18359, GSE13528), the samples corresponding to each tissue were grouped and analyzed separately to facilitate the interpretation of results. The number of GEO experiments in each within each class were: immune: 6 (GSE7313, GSE7314, GSE11787, GSE17492, GSE14758-mediastinal lymph nodes, and GSE14790); embryo and placenta: 5 (GSE18467, GSE18641, GSE18343, GSE11853, and GSE12705); brain and nervous system: 5 (GSE16855, GSE12604, GSE14739-hypothalamus, GSE14739-thyroid, and GSE14739-adenohypophsis); reproduction: 2 (GSE11590, and GSE14739-gonads); muscle: 7 (GSE18653, GSE19275, GSE8974, GSE14643, GSE15211, GSE21096, and GSE16348-skeletal muscle); fat: 8 (GSE17309, GSE14373, GSE14739-fat, GSE9333, GSE18359-fat, GSE18359-liver, GSE13528-fat and GSE13528-liver); gut: 2 (GSE14357 and GSE15256).
The gene expression data were pre-processed and normalized using the Affy R package. Steps included the log-2 transformation and GC-robust multichip average normalization of the gene expression measurements. All probes in the platform were analyzed using ANOVA to identify those that exhibited differential expression across the conditions studied. The false discovery rate approach was used to adjust the statistical significance of the differential expression and account for multiple testing across all probes. The normalization, one or two-way ANOVA and multiple test adjustment of the results were done using Beehive.
Prediction of cleavage sites
The location of the cleavage in pig prohormone proteins that would result in potentially active neuropeptides was predicted using NeuroPred. Complete prohormone sequences from UniProt were used to predict cleavage in preference to the predicted sequences. In limited cases, EST sequences were combined with the genomic data and published partial sequences to predict the complete prohormone sequence. For example, for Chromogranin-A (CMGA), three glutamic acids were missing in the genome-based predictions that were present in the corresponding UniProt fragment sequence ([Swiss-Prot:P04404]) and EST sequence [GenBank:EW261315] permitted the prediction of the complete pig CMGA protein sequence. The location of the potential cleavage sites in the pig prohormones were inferred by homology to human data.
Complementary cleavage prediction models trained on confirmed cleavages from mammalian sequences[12, 78, 79] were used to predict cleavages in the pig prohormone sequences. These models included the known motif model that searches for sites with specific combinations of basic amino acid associated to cleavages reported in other species, mammalian logistic regression, and human logistic regression and artificial neural network models based on amino acids only or amino acids combined with the physicochemical properties of amino acids.
Known or predicted cleavage sites on all 97 prohormone sequences were used to assess the performance of the models to predict cleavage. The “observed” cleavage sites known or inferred from homology to other species based on a literature search[11–13, 19, 78–80] were compared to the cleavage sites predicted by the models. The counts of the true positives (number of correctly predicted cleaved sites), true negatives (number of correctly predicted non-cleaved sites), false positives (the number of incorrectly predicted cleaved sites) and false negatives (number of incorrectly predicted non-cleaved sites) or functions of the counts were used to assess the model performance. These measurements were used to compute the correct classification rate (number of correctly predicted sites divided by the total number of all sites), sensitivity (number of true positives divided by the total number of cleaved sites), specificity (number of true negatives divided by the total number of non-cleaved sites), positive predictive power (number of true positives divided by the total number of sites predicted to be cleaved), negative predictive power (number of true negatives divided by the total number of sites predicted to not be cleaved), Mathew’s correlation coefficient between observed and predicted cleavage. The area under the receiver operator characteristic or ROC curve relating sensitivity and 1 - specificity was also calculated where area values lower than 0.7 indicate poor model performance.
Book SA, Bustad LK: The fetal and neonatal pig in biomedical research. J Anim Sci. 1974, 38 (5): 997-1002.
Anzenbacherova E, Anzenbacher P, Svoboda Z, Ulrichova J, Kvetina J, Zoulova J, Perlik F, Martinkova J: Minipig as a model for drug metabolism in man: comparison of in vitro and in vivo metabolism of propafenone. Biomed Pap Med Fac Univ Palacky Olomouc Czech Repub. 2003, 147 (2): 155-159. 10.5507/bp.2003.021.
Sullivan TP, Eaglstein WH, Davis SC, Mertz P: The pig as a model for human wound healing. Wound Repair Regen. 2001, 9 (2): 66-76. 10.1046/j.1524-475x.2001.00066.x.
Kurahashi M, Miyake H, Takagi T, Tashiro S: Changes of lymphatic flow in case of pancreatic duct obstruction in the pig–as a model of pancreatic cancer. J Med Invest. 2004, 51 (1–2): 70-75.
Fricker LD: Neuropeptide-processing enzymes: applications for drug discovery. AAPS J. 2005, 7 (2): E449-55. 10.1208/aapsj070244.
Hook V, Funkelstein L, Lu D, Bark S, Wegrzyn J, Hwang SR: Proteases for processing proneuropeptides into peptide neurotransmitters and hormones. Annu Rev Pharmacol Toxicol. 2008, 48: 393-423. 10.1146/annurev.pharmtox.48.113006.094812.
Schwartz PH, Nethercott H, Kirov II, Ziaeian B, Young MJ, Klassen H: Expression of neurodevelopmental markers by cultured porcine neural precursor cells. Stem Cells. 2005, 23 (9): 1286-1294. 10.1634/stemcells.2004-0306.
Yang G, Su J, Yao Y, Lei Z, Zhang G, Li X: The regulatory mechanism of neuromedin S on luteinizing hormone in pigs. Anim Reprod Sci. 2010, 122 (3–4): 367-374.
Kaminski T, Smolinska N, Nitkiewicz A, Przala J: Expression of orexin receptors 1 (OX1R) and 2 (OX2R) in the porcine hypothalamus during the oestrous cycle. J Physiol Pharmacol. 2010, 61 (3): 363-371.
Southey BR, Sweedler JV, Rodriguez-Zas SL: Prediction of neuropeptide cleavage sites in insects. Bioinformatics. 2008, 24 (6): 815-825. 10.1093/bioinformatics/btn044.
Southey BR, Sweedler JV, Rodriguez-Zas SL: A Python analytical pipeline to identify prohormone precursors and predict prohormone cleavage sites. Front Neuroinform. 2008, 2: 7-
Tegge AN, Southey BR, Sweedler JV, Rodriguez-Zas SL: Comparative analysis of neuropeptide cleavage sites in human, mouse, rat, and cattle. Mamm Genome. 2008, 19 (2): 106-120. 10.1007/s00335-007-9090-9.
Southey BR, Rodriguez-Zas SL, Sweedler JV: Characterization of the prohormone complement in cattle using genomic libraries and cleavage prediction approaches. BMC Genomics. 2009, 10: 228-10.1186/1471-2164-10-228.
Ross JW, Ashworth MD, White FJ, Johnson GA, Ayoubi PJ, DeSilva U, Whitworth KM, Prather RS, Geisert RD: Premature estrogen exposure alters endometrial gene expression to disrupt pregnancy in the pig. Endocrinology. 2007, 148 (10): 4761-4773. 10.1210/en.2007-0599.
Hausman GJ, Barb CR, Dean RG: Patterns of gene expression in pig adipose tissue: insulin-like growth factor system proteins, neuropeptide Y (NPY), NPY receptors, neurotrophic factors and other secreted factors. Domest Anim Endocrinol. 2008, 35 (1): 24-34. 10.1016/j.domaniend.2008.01.004.
Sus scrofa (pig) genome view: http://www.ncbi.nlm.nih.gov/projects/mapview/map_search.cgi?taxid=9823.
Li L, Sweedler JV: Peptides in the brain: mass spectrometry-based measurement approaches and challenges. Annu Rev Anal Chem (Palo Alto Calif). 2008, 1: 451-483. 10.1146/annurev.anchem.1.031207.113053.
Hummon AB, Richmond TA, Verleyen P, Baggerman G, Huybrechts J, Ewing MA, Vierstraete E, Rodriguez-Zas SL, Schoofs L, Robinson GE, Sweedler JV: From the genome to the proteome: uncovering peptides in the Apis brain. Science. 2006, 314 (5799): 647-649. 10.1126/science.1124128.
Xie F, London SE, Southey BR, Annangudi SP, Amare A, Rodriguez-Zas SL, Clayton DF, Sweedler JV: The zebra finch neuropeptidome: prediction, detection and expression. BMC Biol. 2010, 8: 28-10.1186/1741-7007-8-28.
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 (3): 140-148.
Chen H, Li C, Fang M, Zhu M, Li X, Zhou R, Li K, Zhao S: Understanding Haemophilus parasuis infection in porcine spleen through a transcriptomics approach. BMC Genomics. 2009, 10: 64-10.1186/1471-2164-10-64.
Saida K, Uchide T, Zhang S, Adur J, Kawano Y, Ogiso M, Oka S, Takizawa S: Molecular cloning and sequence analysis of the porcine precursor of endothelin-2. J Cardiovasc Pharmacol. 2004, 44 (Suppl 1): S244-7.
Rezaeian AH, Katafuchi T, Yoshizawa M, Hiraiwa N, Saito T, Nishibori M, Hamano K, Minamino N, Yasue H: Genomic organization, expression and evolution of porcine CRSP1, 2, and 3. Cytogenet Genome Res. 2008, 121 (1): 41-49. 10.1159/000124380.
Yao Y, Lin X, Su J, Yang G, Hou Y, Lei Z: Cloning and distribution of neuropeptide S and its receptor in the pig. Neuropeptides. 2009, 43 (6): 465-481. 10.1016/j.npep.2009.09.003.
Hong Y, Hay DL, Quirion R, Poyner DR: The pharmacology of Adrenomedullin 2/Intermedin. Br J Pharmacol. 2012, 166 (1): 110-120. 10.1111/j.1476-5381.2011.01530.x.
Spier AD, de Lecea L: Cortistatin: a member of the somatostatin neuropeptide family with distinct physiological functions. Brain Res Brain Res Rev. 2000, 33 (2–3): 228-241.
Haugaard-Jonsson LM, Hossain MA, Daly NL, Craik DJ, Wade JD, Rosengren KJ: Structure of human insulin-like peptide 5 and characterization of conserved hydrogen bonds and electrostatic interactions within the relaxin framework. Biochem J. 2009, 419 (3): 619-627. 10.1042/BJ20082353.
Ukena K, Tachibana T, Iwakoshi-Ukena E, Saito Y, Minakata H, Kawaguchi R, Osugi T, Tobari Y, Leprince J, Vaudry H, Tsutsui K: Identification, localization, and function of a novel avian hypothalamic neuropeptide, 26RFa, and its cognate receptor, G protein-coupled receptor-103. Endocrinology. 2010, 151 (5): 2255-2264. 10.1210/en.2009-1478.
Takayanagi Y, Onaka T: Roles of prolactin-releasing peptide and RFamide related peptides in the control of stress and food intake. FEBS J. 2010, 277 (24): 4998-5005. 10.1111/j.1742-4658.2010.07932.x.
Monnier J, Samson M: Prokineticins in angiogenesis and cancer. Cancer Lett. 2010, 296 (2): 144-149. 10.1016/j.canlet.2010.06.011.
Reinscheid RK: Neuropeptide S: anatomy, pharmacology, genetics and physiological functions. Results Probl Cell Differ. 2008, 46: 145-158. 10.1007/400_2007_051.
Dobolyi A, Palkovits M, Usdin TB: The TIP39-PTH2 receptor system: unique peptidergic cell groups in the brainstem and their interactions with central regulatory mechanisms. Prog Neurobiol. 2010, 90 (1): 29-59. 10.1016/j.pneurobio.2009.10.017.
Venkatasubramanian S, Newby DE, Lang NN: Urocortins in heart failure. Biochem Pharmacol. 2010, 80 (3): 289-296. 10.1016/j.bcp.2010.03.032.
Vaudry H: Biologically active peptides urotensin II (UII) and urotensin II-related peptide (URP), and to their cognate receptor (UT). Editorial. Peptides. 2008, 29 (5): 647-648. 10.1016/j.peptides.2008.01.010.
Dai G, Smeekens SP, Steiner DF, McMurtry JP, Kwok SC: Characterization of multiple prohormone convertase PC1/3 transcripts in porcine ovary. Biochim Biophys Acta. 1995, 1264 (1): 1-6. 10.1016/0167-4781(95)00111-S.
Renegar RH, Owens CR, Dai G, McMurtry JP, Kwok SC: Expression and localization of prohormone convertase 1/3 (SPC3) in porcine ovary. Mol Reprod Dev. 2000, 57 (4): 361-365. 10.1002/1098-2795(200012)57:4<361::AID-MRD8>3.0.CO;2-T.
Seidah NG, Fournier H, Boileau G, Benjannet S, Rondeau N, Chretien M: The cDNA structure of the porcine pro-hormone convertase PC2 and the comparative processing by PC1 and PC2 of the N-terminal glycopeptide segment of porcine POMC. FEBS Lett. 1992, 310 (3): 235-239. 10.1016/0014-5793(92)81339-N.
Jessop DS: Neuropeptides in the Immune System: Mediators of Stress and Inflammation. Handbook of Neurochemistry and Molecular Neurobiology Part 1. Edited by: Lajtha A, Galoyan A, Besedovsky HO. 2008, New York: Springer, 19-35.
Pampusch MS, Serie JR, Osinski MA, Seybold VS, Murtaugh MP, Brown DR: Expression of nociceptin/OFQ receptor and prepro-nociceptin/OFQ in lymphoid tissues. Peptides. 2000, 21 (12): 1865-1870. 10.1016/S0196-9781(00)00332-6.
Wang Y, Qu L, Uthe JJ, Bearson SM, Kuhar D, Lunney JK, Couture OP, Nettleton D, Dekkers JC, Tuggle CK: Global transcriptional response of porcine mesenteric lymph nodes to Salmonella enterica serovar Typhimurium. Genomics. 2007, 90 (1): 72-84. 10.1016/j.ygeno.2007.03.018.
Tomas A, Fernandes LT, Sanchez A, Segales J: Time course differential gene expression in response to porcine circovirus type 2 subclinical infection. Vet Res. 2010, 41 (1): 12-10.1051/vetres/2009060.
Gonzalez AM, Podvin S, Lin SY, Miller MC, Botfield H, Leadbeater WE, Roberton A, Dang X, Knowling SE, Cardenas-Galindo E, Donahue JE, Stopa EG, Johanson CE, Coimbra R, Eliceiri BP, Baird A: Ecrg4 expression and its product augurin in the choroid plexus: impact on fetal brain development, cerebrospinal fluid homeostasis and neuroprogenitor cell response to CNS injury. Fluids Barriers CNS. 2011, 8 (1): 6-10.1186/2045-8118-8-6.
GSE7314: Expression data from non-infected and Salmonella Choleraesuis infected mesenteric lymph nodes. [http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc = GSE7314],
Solinhac R, Mompart F, Martin P, Robelin D, Pinton P, Iannuccelli E, Lahbib-Mansais Y, Oswald IP, Yerle-Bouissou M: Transcriptomic and nuclear architecture of immune cells after LPS activation. Chromosoma. 2011, 120 (5): 501-520. 10.1007/s00412-011-0328-7.
Radek KA, Lopez-Garcia B, Hupe M, Niesman IR, Elias PM, Taupenot L, Mahata SK, O’Connor DT, Gallo RL: The neuroendocrine peptide catestatin is a cutaneous antimicrobial and induced in the skin after injury. J Invest Dermatol. 2008, 128 (6): 1525-1534. 10.1038/sj.jid.5701225.
Kromer A, Glombik MM, Huttner WB, Gerdes HH: Essential role of the disulfide-bonded loop of chromogranin B for sorting to secretory granules is revealed by expression of a deletion mutant in the absence of endogenous granin synthesis. J Cell Biol. 1998, 140 (6): 1331-1346. 10.1083/jcb.140.6.1331.
Tayade C, Fang Y, Croy BA: A review of gene expression in porcine endometrial lymphocytes, endothelium and trophoblast during pregnancy success and failure. J Reprod Dev. 2007, 53 (3): 455-463. 10.1262/jrd.18170.
Ericsson A, Hemsen A, Lundberg JM, Persson H: Detection of neuropeptide Y-like immunoreactivity and messenger RNA in rat platelets: the effects of vinblastine, reserpine, and dexamethasone on NPY expression in blood cells. Exp Cell Res. 1991, 192 (2): 604-611. 10.1016/0014-4827(91)90082-6.
Almeida LM, Basu U, Williams JL, Moore SS, Guan LL: Microarray analysis in caudal medulla of cattle orally challenged with bovine spongiform encephalopathy. Genet Mol Res. 2011, 10 (4): 3948-3962. 10.4238/2011.October.25.5.
Elsasser TH, Kahl S: Adrenomedullin has multiple roles in disease stress: development and remission of the inflammatory response. Microsc Res Tech. 2002, 57 (2): 120-129. 10.1002/jemt.10058.
Turpeinen H, Raitoharju E, Oksanen A, Oksala N, Levula M, Lyytikainen LP, Jarvinen O, Creemers JW, Kahonen M, Laaksonen R, Pelto-Huikko M, Lehtimaki T, Pesu M: Proprotein convertases in human atherosclerotic plaques: the overexpression of FURIN and its substrate cytokines BAFF and APRIL. Atherosclerosis. 2011, 219 (2): 799-806. 10.1016/j.atherosclerosis.2011.08.011.
Ostrup E, Bauersachs S, Blum H, Wolf E, Hyttel P: Differential endometrial gene expression in pregnant and nonpregnant sows. Biol Reprod. 2010, 83 (2): 277-285. 10.1095/biolreprod.109.082321.
Pantaleon M, Jericho H, Rabnott G, Kaye PL: The role of insulin-like growth factor II and its receptor in mouse preimplantation development. Reprod Fertil Dev. 2003, 15 (1–2): 37-45.
Gupta MK, Uhm SJ, Han DW, Lee HT: Embryo quality and production efficiency of porcine parthenotes is improved by phytohemagglutinin. Mol Reprod Dev. 2007, 74 (4): 435-444. 10.1002/mrd.20547.
Gupta MK, Uhm SJ, Lee SH, Lee HT: Role of nonessential amino acids on porcine embryos produced by parthenogenesis or somatic cell nuclear transfer. Mol Reprod Dev. 2008, 75 (4): 588-597. 10.1002/mrd.20789.
Pittius CW, Ellendorff F, Hollt V, Parvizi N: Ontogenetic development of proenkephalin A and proenkephalin B messenger RNA in fetal pigs. Exp Brain Res. 1987, 69 (1): 208-212.
Anisimov A, Alitalo A, Korpisalo P, Soronen J, Kaijalainen S, Leppanen VM, Jeltsch M, Yla-Herttuala S, Alitalo K: Activated forms of VEGF-C and VEGF-D provide improved vascular function in skeletal muscle. Circ Res. 2009, 104 (11): 1302-1312. 10.1161/CIRCRESAHA.109.197830.
Cao Y, Linden P, Farnebo J, Cao R, Eriksson A, Kumar V, Qi JH, Claesson-Welsh L, Alitalo K: Vascular endothelial growth factor C induces angiogenesis in vivo. Proc Natl Acad Sci U S A. 1998, 95 (24): 14389-14394. 10.1073/pnas.95.24.14389.
Karkkainen MJ, Haiko P, Sainio K, Partanen J, Taipale J, Petrova TV, Jeltsch M, Jackson DG, Talikka M, Rauvala H, Betsholtz C, Alitalo K: Vascular endothelial growth factor C is required for sprouting of the first lymphatic vessels from embryonic veins. Nat Immunol. 2004, 5 (1): 74-80. 10.1038/ni1013.
Uebing-Czipura AU, Dawson HD, Rutherford MS, Scherba G: Transcriptome profile and cytogenetic analysis of immortalized neuronally restricted progenitor cells derived from the porcine olfactory bulb. Anim Biotechnol. 2009, 20 (4): 186-215. 10.1080/10495390903139950.
Vergano-Vera E, Fernandez AP, Hurtado-Chong A, Vicario-Abejon C, Martinez A: Lack of adrenomedullin affects growth and differentiation of adult neural stem/progenitor cells. Cell Tissue Res. 2010, 340 (1): 1-11. 10.1007/s00441-010-0934-3.
Mathonnet M, Comte I, Lalloue F, Ayer-Le Lievre C: Insulin-like growth factor I induced survival of axotomized olfactory neurons in the chick. Neurosci Lett. 2001, 308 (2): 67-70. 10.1016/S0304-3940(01)01715-3.
Perez-Enciso M, Ferraz AL, Ojeda A, Lopez-Bejar M: Impact of breed and sex on porcine endocrine transcriptome: a bayesian biometrical analysis. BMC Genomics. 2009, 10: 89-10.1186/1471-2164-10-89.
Yang B, Bassols A, Saco Y, Perez-Enciso M: Association between plasma metabolites and gene expression profiles in five porcine endocrine tissues. Genet Sel Evol. 2011, 43: 28-10.1186/1297-9686-43-28.
Busch-Sorensen M, Sheikh SP, O’Hare M, Tortora O, Schwartz TW, Gammeltoft S: Regional distribution of neuropeptide Y and its receptor in the porcine central nervous system. J Neurochem. 1989, 52 (5): 1545-1552. 10.1111/j.1471-4159.1989.tb09206.x.
Blakemore LJ, Levenson CW, Trombley PQ: Neuropeptide Y modulates excitatory synaptic transmission in the olfactory bulb. Neuroscience. 2006, 138 (2): 663-674. 10.1016/j.neuroscience.2005.11.027.
Le Bras B, Barallobre MJ, Homman-Ludiye J, Ny A, Wyns S, Tammela T, Haiko P, Karkkainen MJ, Yuan L, Muriel MP, Chatzopoulou E, Breant C, Zalc B, Carmeliet P, Alitalo K, Eichmann A, Thomas JL: VEGF-C is a trophic factor for neural progenitors in the vertebrate embryonic brain. Nat Neurosci. 2006, 9 (3): 340-348. 10.1038/nn1646.
Zenmyo M, Komiya S, Hamada T, Hiraoka K, Suzuki R, Inoue A: P21 and Parathyroid Hormone-Related Peptide in the Growth Plate. Calcif Tissue Int. 2000, 67 (5): 378-381. 10.1007/s002230001153.
Fressinaud C, Weinrauder H, Delaunoy JP, Tholey G, Labourdette G, Sarlieve LL: Glutamine synthetase expression in rat oligodendrocytes in culture: regulation by hormones and growth factors. J Cell Physiol. 1991, 149 (3): 459-468. 10.1002/jcp.1041490315.
Wu Q: The serine protease corin in cardiovascular biology and disease. Front Biosci. 2007, 12: 4179-4190. 10.2741/2379.
Fekete EM, Zorrilla EP: Physiology, pharmacology, and therapeutic relevance of urocortins in mammals: ancient CRF paralogs. Front Neuroendocrinol. 2007, 28 (1): 1-27. 10.1016/j.yfrne.2006.09.002.
Southey BR, Amare A, Zimmerman TA, Rodriguez-Zas SL, Sweedler JV: NeuroPred: a tool to predict cleavage sites in neuropeptide precursors and provide the masses of the resulting peptides. Nucleic Acids Res. 2006, 34 (Web Server issue): W267-W272.
Bloomquist BT, Darlington DN, Mueller GP, Mains RE, Eipper BA: Regulated endocrine-specific protein-18: a short-lived novel glucocorticoid-regulated endocrine protein. Endocrinology. 1994, 135 (6): 2714-2722. 10.1210/en.135.6.2714.
Potter KJ, Abedini A, Marek P, Klimek AM, Butterworth S, Driscoll M, Baker R, Nilsson MR, Warnock GL, Oberholzer J, Bertera S, Trucco M, Korbutt GS, Fraser PE, Raleigh DP, Verchere CB: Islet amyloid deposition limits the viability of human islet grafts but not porcine islet grafts. Proc Natl Acad Sci U S A. 2010, 107 (9): 4305-4310. 10.1073/pnas.0909024107.
Maglott D, Ostell J, Pruitt KD, Tatusova T: Entrez Gene: gene-centered information at NCBI. Nucleic Acids Res. 2011, 39 (Database issue): D52-D57.
Sayers EW, Barrett T, Benson DA, Bolton E, Bryant SH, Canese K, Chetvernin V, Church DM, Dicuccio M, Federhen S, Feolo M, Fingerman IM, Geer LY, Helmberg W, Kapustin Y, Krasnov S, Landsman D, Lipman DJ, Lu Z, Madden TL, Madej T, Maglott DR, Marchler-Bauer A, Miller V, Karsch-Mizrachi I, Ostell J, Panchenko A, Phan L, Pruitt KD, Schuler GD, Sequeira E, Sherry ST, Shumway M, Sirotkin K, Slotta D, Souvorov A, Starchenko G, Tatusova TA, Wagner L, Wang Y, Wilbur WJ, Yaschenko E, Ye J: Database resources of the National Center for Biotechnology Information. Nucleic Acids Res. 2012, 40 (Database issue): D13-D25.
UniProt Consortium: Reorganizing the protein space at the Universal Protein Resource (UniProt). Nucleic Acids Res. 2012, 40 (Database issue): D71-D75.
Southey BR, Rodriguez-Zas SL, Sweedler JV: Prediction of neuropeptide prohormone cleavages with application to RFamides. Peptides. 2006, 27 (5): 1087-1098. 10.1016/j.peptides.2005.07.026.
Amare A, Hummon AB, Southey BR, Zimmerman TA, Rodriguez-Zas SL, Sweedler JV: Bridging neuropeptidomics and genomics with bioinformatics: Prediction of mammalian neuropeptide prohormone processing. J Proteome Res. 2006, 5 (5): 1162-1167. 10.1021/pr0504541.
Delfino KR, Southey BR, Sweedler JV, Rodriguez-Zas SL: Genome-wide census and expression profiling of chicken neuropeptide and prohormone convertase genes. Neuropeptides. 2010, 44 (1): 31-44. 10.1016/j.npep.2009.11.002.
Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25 (17): 3389-3402. 10.1093/nar/25.17.3389.
Birney E, Clamp M, Durbin R: GeneWise and Genomewise. Genome Res. 2004, 14 (5): 988-995. 10.1101/gr.1865504.
Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA, McWilliam H, Valentin F, Wallace IM, Wilm A, Lopez R, Thompson JD, Gibson TJ, Higgins DG: Clustal W and Clustal X version 2.0. Bioinformatics. 2007, 23 (21): 2947-2948. 10.1093/bioinformatics/btm404.
Barrett T, Troup DB, Wilhite SE, Ledoux P, Evangelista C, Kim IF, Tomashevsky M, Marshall KA, Phillippy KH, Sherman PM, Muertter RN, Holko M, Ayanbule O, Yefanov A, Soboleva A: NCBI GEO: archive for functional genomics data sets--10 years on. Nucleic Acids Res. 2011, 39 (Database issue): D1005-D1010.
Gautier L, Cope L, Bolstad BM, Irizarry RA: affy–analysis of Affymetrix GeneChip data at the probe level. Bioinformatics. 2004, 20 (3): 307-315. 10.1093/bioinformatics/btg405.
Benjamini Y, Hochberg Y: Controlling the false discovery rate: A practical and powerful approach to multiple testing. J R Statist Soc Series B. 1995, 57: 289-300.
Beehive: An integrated suite of tools to studymicroarray experiments. [http://stagbeetle.animal.uiuc.edu/Beehive],
The project described was supported by Award Numbers R21DA027548 and P30DA018310 from the National Institute on Drug Abuse (NIDA), Award Number 1R03CA143975 from the National Cancer Institute (NCI), and by Award Number ILLU-538-311 from the United States Department of Agriculture (USDA), Cooperative State Research, Education and Extension Service.
The authors declare that they have no competing interests.
KIP performed the search for prohormone and prohormone convertase genes in the UniProt, Gene, UniGene and ENSEMBL databases, identified the gene probes on the microarray platform, analyzed 35 microarray experiments, contributed to the interpretation of results, and drafted the manuscript. BRS located the prohormone and prohormone convertase genes in the pig genome assembly, trace archives and EST databases, predicted the genes from the genome sequence, compared the prediction to known sequences, contributed to the interpretation of results, and manuscript. JVS obtained funding for the study, contributed to the interpretation of the results, and reviewed the manuscript. SRZ obtained funding for the study, participated in its conception, coordination, interpretation of results, and helped write the manuscript. All authors have read and approved the final version of this manuscript.
Electronic supplementary material
Additional file 1: Table S2. Prohormone and convertase genes identified across multiple pig genome resources. (DOC 79 KB)
Additional file 3: Table S3. Statistical significance P-value corresponding to the comparison between groups within experiment. (DOC 340 KB)
About this article
Cite this article
Porter, K.I., Southey, B.R., Sweedler, J.V. et al. First survey and functional annotation of prohormone and convertase genes in the pig. BMC Genomics 13, 582 (2012). https://doi.org/10.1186/1471-2164-13-582