- Open Access
Transcriptome analysis of the effect of pyrroloquinoline quinone disodium (PQQ·Na2) on reproductive performance in sows during gestation and lactation
Journal of Animal Science and Biotechnologyvolume 10, Article number: 62 (2019)
Pyrroloquinoline quinone (PQQ), which is a water soluble, thermo-stable triglyceride-quinone, is widely distributed in nature and characterized as a mammalian vitamin-like redox cofactor. The objective of this study was to investigate the effects of pyrroloquinoline quinone disodium (PQQ·Na2) on reproductive performance in sows.
Dietary supplementation with PQQ·Na2 significantly increased the total number of piglets born, the number of piglets born alive and the born alive litter weight. It also increased the antioxidant status in the placenta, plasma and milk. The concentration of NO was significantly increased in the plasma and placenta. RNA-seq analysis showed that 462 unigenes were differentially expressed between the control (Con) treatment and PQQ treatment groups. Among these unigenes, 199 were upregulated, while 263 unigenes were downregulated. The assigned functions of the unigenes covered a broad range of GO categories. Reproduction (27, 7.03%) and the reproduction process (27, 7.03%) were assigned to the biological process category. By matching DEGs to the KEGG database, we identified 29 pathways.
In conclusion, dietary supplementation with PQQ·Na2 in gestating and lactating sows had positive effects on their reproductive performance.
Sow reproductive performance affects the production level and provides an economic benefit to the pig industry. Litter size is an important reproductive trait and a critical indicator of sow reproduction performance . Maintaining optimal reproductive and litter performance is essential for meeting economic targets in commercial pig production . Increasing litter size has long been a goal of pig breeders and producers; prolificacy is of great interest to the pig industry . Maternal nutrition has substantial implications for fetal health. Enhancing reproductive performance through nutrition and management strategies in gestation and lactation sows have been of research interest for several decades . The placental tissue is the only site for contact between the fetus and the mother during pregnancy; thus, the tissue is closely related to the health and development of the fetus . Maternal conditions have been demonstrated to affect the placental morphology, blood flow, fetomaternal exchanges, and endocrine function  . In addition, the placenta plays a pivotal role in maternal nutrient supply and metabolic waste removal, protection against bacterial and viral infections, and production of hormones supporting pregnancy .
Pyrroloquinoline quinone (PQQ), which is a water soluble, thermo-stable triglyceride-quinone . Initially identified as a novel cofactor of various bacterial dehydrogenases , PQQ is an essential animal nutrient. PQQ-deficient animals display a variety of illnesses . PQQ has attracted considerable attention, as it is important for mammalian growth, development, reproduction and immune function . PQQ is an effective antioxidant, protecting mitochondria against oxidative stress-induced lipid peroxidation, protein carbonyl formation and inactivation of the mitochondrial respiratory chain . Although PQQ is not biosynthesized in mammals, trace amounts of PQQ have been found in human and rat tissues at picomolar to nanomolar levels, and an especially large amount has been found in human milk .
Because of its versatile functions, PQQ disodium (PQQ·Na2) salt has been authorized in Canada as a Natural Health Product, providing 20 mg PQQ·Na2 salt per day as an antioxidant for the maintenance of good health . On August 13, 2018, the European Commission issued regulations (EU) 2018/1122 approving pyrroquinone sodium salt as a new type of food. The European commission defines pyrroquinone sodium as a dietary supplement. Although PQQ has a positive effect on reproductive performance, its mechanism is not clear. Therefore, in this study, we used the Illumina HiSeqTM2500 platform to perform a large-scale transcriptome analysis of the placenta of sows. The next-generation RNA sequencing (RNA-seq) platform has emerged as the method of choice for studying transcriptomes . Second generation sequencing technology, also called RNA sequencing (RNA-seq), is powerful for gene identification, comparative gene expression analysis and investigation of the functional complexity of the transcriptome . In recent years, an RNA-seq approach has been widely used in animals for novel gene identification and differentially expressed gene (DEGs) analysis, because it is high throughput, low cost, covers a multitude of low abundance genes, and has high sensitivity. Although PQQ has gained interest in medicine and health-related research in recent decades, the usefulness of PQQ has not yet been fully demonstrated in animal agriculture, especially in the pig industry. In addition, there are no published data on the effects of PQQ on reproductive performance in sows. Therefore, the objective of this study was to test the efficacy of dietary PQQ·Na2 supplementation on reproductive performance. We also analyzed the RNA-seq of placentas, which revealed the genes that may be involved in placental development and function, thus playing a role in determining litter size. Furthermore, the data collected were used to establish the relationship between PQQ and the placenta of sows, as well as reproductive performance. We hope the results of this study could lay a foundation for the further study of PQQ in the pig industry.
Animals and management
A total of 40 crossbred (Landrace×Large White crossed with Duroc boar) multiparity gestation sows with an average parity of 4.3 were used in the study. Forty sows were allotted to 2 dietary treatments after breeding. One group was the control sows, which were fed a corn-soybean meal control diet (Con treatment, n = 20), and the other group was the treatment sows, fed a control diet with 20 mg/kg PQQ·Na2 after breeding and through gestation (PQQ treatment, n = 20). The PQQ·Na2 (purity, ≥ 98%) was synthesized by chemical reactions. It was diluted with corn starch to a concentration of 1 g/kg mixture before being mixed into the diet. Based on the known range of PQQ in foods , we inferred that the concentration of PQQ in the basal diet was less than 0.01 mg/kg. The sows were kept in single crates (0.6 m × 2.0 m) from insemination to day 110 of gestation. On d 110 of gestation, sows were transported to the farrowing facility, where they were placed in individual farrowing crates (2.4 m × 2.4 m). Each crate had steel mesh floors with a heat lamp for newborn pigs. The crates were mounted over a solid concrete floor, and manure was removed manually each day. The farrowing room temperature was maintained at approximately 18 to 20 °C. Births were watched, but the observers interfered as little as possible in the farrowing process. The protocols used in this experiment were approved by the Northeast Agricultural University Institutional Animal Care and Use Committee. All animal experimental diets (Table 1) were formulated to meet or exceed the recommended nutrient requirements of the NRC (2012). From d 1 of gestation until d 90 of gestation, all sows were fed 2.5 kg of the gestation diet daily. From d 91 of gestation, all sows were fed 4.5 kg of the gestation diet daily. The amounts of parturition feed provided to each sow at d 112, 113 and 114 of gestation were 2.0, 1.5 and 1.0 kg, respectively.
Sample and data collection
Blood was collected from the ear vein of a random subset of sows (n = 8 per treatment) at d 90 of gestation and d 21 of lactation. The blood was centrifuged at 3000×g for 15 min to obtain the plasma, and the plasma was stored at − 20 °C until analysis. Eight sows per treatment were randomly selected and marked for milk sample collection during lactation. Colostrum was collected within 6 h of farrowing (d 0 of lactation). Approximately 30 to 50 mL of milk was collected from all functional mammary glands using a mechanical milk pump after the injection of 1 mL oxytocin. The samples were immediately stored at − 20 °C until analysis. Placenta allantochorion tissue samples were collected immediately during parturition to preserve RNA stability for mRNA analysis. A section of samples was stored at − 20 °C, and another section was snap-frozen in liquid nitrogen for further analysis. Sow back-fat thickness was measured at d 0 and 90 of gestation, within 24 h of farrowing (d 0) and d 21 of lactation (n = 8 per treatment). Back-fat thickness was measured at the P2 position (left side of the 10th rib and 6 cm lateral to the spine) by digital B-ultrasonography (Kaixin, Xuzhou, China). At farrowing, the number of piglets born, litter birth weight and individual birth weights were recorded.
Evaluation of antioxidant enzyme activity and the concentration of nitric oxide
Superoxide dismutase (SOD), glutathione peroxidase (GSH-Px) and catalase (CAT) enzyme activities in the milk, plasma and placenta were determined using commercially available kits (Nanjing Jiancheng Bioengineering Institute, Nanjing, China) in accordance with our previous study. The results of the measurements were expressed as U/mL in plasma and milk and as U/mg protein in placenta. Lipid peroxidation in the plasma, milk and placenta was determined by measuring the amounts of malondialdehyde (MDA) through the thiobarbituric acid method using commercially available kits (Nanjing Jiancheng Bioengineering Institute, Nanjing, China). The results of the measurements were expressed as nmol/mL in plasma and as nmol/mg protein in placenta. Nitric oxide (NO) and inducible NOS (iNOS) in the plasma and placenta of sows were determined using assay kits obtained from Jiancheng Biochemistry (Nanjing, China). The colostrum was analyzed for lactose, protein, fat, and total solids with a fully automatic milk analyzer (Milko ScanTM FT+ Analyzer, Foss). The milk samples were analyzed for immunoglobulin G (IgG), immunoglobulin A (IgA), and immunoglobulin M (IgM) using immunoglobulin-specific kits (Jinma Biological Engineering Co., Ltd., Shanghai, China).
RNA extraction, cDNA library construction and sequencing
Total RNA was isolated from the placenta using TRIzol reagent (Invitrogen, USA) according to the manufacturer’s instructions. The extracted RNA was treated with DNase I (Takara Biotechnology, China) for 45 min at 37 °C to remove residual DNA. The RNA concentration and integrity were measured using an UltrasecTM 2100 pro UV/Visible Spectrophotometer (Amersham Biosciences, Uppsala, Sweden) and gel electrophoresis. Equal amounts of high-quality RNA from each specimen were pooled for RNA-Seq library construction. A cDNA library was prepared with a TruSeq RNA sample preparation kit following the manufacturer’s instructions (Illumina) and sequenced on an Illumina HiSeqTM 2500 platform in 100 pair-ended mode (Biomarker Technologies).
De novo transcriptome assembly and functional annotation
To obtain clean reads, the raw reads were filtered by removing the adapter, poly-N and low-quality sequences. De novo assembly was performed using the Trinity method . The K-mer and group pairs distance were set at 25 and 300, respectively, while the other parameters were set at default levels. Based on their overlap regions, the short reads were assembled into longer contigs, which were then clustered and further assembled into unigenes with the paired-end information. For gene functional annotation, all of the assembled transcripts were aligned to the publicly available protein databases, including the National Center for Biotechnology Information (NCBI), nonredundant protein (Nr), the Swiss-Prot protein, Gene Ontology (GO) (http://wego.genomics.org.cn/cgi-bin/wego/index.pl), Clusters of Orthologous Groups (COG), and the Kyoto Encyclopedia of Genes and Genomes (KEGG) (http://www.genome.jp/kegg/kegg2.html) using BLASTx analysis with a cut-off E-value of 10− 5.
Differentially expressed genes (DEGs) analysis
Fragments per kilobase of transcript per million fragments mapped (FPKM) was calculated to represent the expression abundance of the unigenes. FPKM may reflect the molar concentration of a transcript by normalizing for RNA length and for the total read number. DEGs between PQQ treated and control samples were identified by EBSeq. An FDR (false discovery rate) < 0.05 and |fold change (FC)| ≥ 2 was set as the threshold for significantly different expression
Quantitative real-time PCR (qRT-PCR)
Total RNA from each sample was converted into cDNA using the Prime Script RT reagent Kit (TaKaRa Bio Catalog), and the cDNA was used for qRT-PCR. GAPDH was used as an internal control gene, and it did not respond to dietary treatments. The primer sequences are shown in Table 2. qRT-PCR was performed using the SYBR Green I Kit (TaKaRa Bio Catalog). For analyses, using an ABI PRISM 7500 SDS thermal cycler, PCRs were performed with 2.0 mL of first-strand cDNA and 0.4 mL of forward and reverse primers in a final volume of 20 mL. Samples were centrifuged briefly and run on the PCR machine using the default fast program (1 cycle at 95 °C for 30 s, 40 cycles of 95 °C for 5 s and 60 °C for 34 s). All of the PCRs were performed in triplicate. The relative gene expression levels were determined using the 2–ΔΔCt method.
All data analyses were performed with SPSS 19.0 software (IBM-SPSS Inc., Chicago, Illinois, USA). The data were analyzed by One-way analysis of variance (ANOVA), and multiple comparisons were analyzed with the Tukey’s test in SPSS. The individual sow and her litter were defined as the experimental unit. The results were presented as the mean values and the standard error of the mean (SEM). In all statistical tests used, P < 0.05 was considered significantly different.
The results pertaining to the reproductive performance of sows are shown in Table 3. The back-fat thickness of sows at d 0 and 90 of gestation and at weaning did not differ between the two dietary treatments (P > 0.05). Dietary PQQ·Na2 supplementation in the gestation and lactation diets had no effect on back-fat thickness gain during gestation or loss during lactation (P > 0.05). The total piglets born, number of piglets born alive and born alive litter weight were significantly increased (P < 0.05) by dietary PQQ·Na2 supplementation during gestation and lactation. The number of piglets weaned, litter weaning weight and piglet weaning weight were not affected by dietary PQQ·Na2 during gestation and lactation. (P > 0.05).
Antioxidant status and the concentration of nitric oxide
The effects of dietary PQQ·Na2 supplementation during gestation and lactation on the antioxidant capacity in the placenta of sows are shown in Fig. 1. PQQ·Na2 dietary supplementation in gestation and lactation significantly increased the activities of SOD (P = 0.011) and GSH-Px (P = 0.032). The MDA (P = 0.018) activity in the placenta after PQQ·Na2 supplementation was significantly lower than the control treatment. However, the CAT activity of the placenta was not significantly different (P > 0.05). The results of the antioxidant status in the plasma of sows are presented in Fig. 2. On d 90 of gestation, the activities of CAT (P = 0.029), SOD (P = 0.041) and GSH-Px (P = 0.030) were significantly increased by PQQ·Na2 supplementation. The MDA (P > 0.05) activity of the placenta was not significantly changed. On d 21 of lactation, the SOD (P = 0.027) and GSH-Px (P = 0.023) activities were significantly increased and the MDA (P = 0.023) activity was significantly decreased by PQQ·Na2 supplementation. The concentration of NO (P = 0.020) and iNOS (P = 0.037) were significantly increased in plasma on d 90 of gestation by PQQ·Na2 supplementation are shown in Fig. 3. PQQ·Na2 dietary supplementation in gestation and lactation significantly increased the concentration of NO (P = 0.034) and iNOS (P = 0.016) in the placenta of sows.
Figure 4 shows the protein, lactose, fat and total milk solids content of the sow milk. The concentrations of protein (P = 0.011) and total solids (P = 0.040) were significantly increased by PQQ·Na2 supplementation. The concentrations of fat and lactose were not significant (P > 0.05). The SOD (P = 0.011) and GSH-Px (P = 0.016) activity in milk from the PQQ sows was higher than that in milk from the Con sows, as shown in Fig. 5. The effects of dietary PQQ·Na2 supplementation during gestation and lactation on immunoglobulin concentrations in the colostrum of the sows are presented in Fig. 6. The concentrations of IgA, IgG and IgM were significantly (P < 0.05) increased in the colostrum with PQQ·Na2 supplementation.
High-throughput transcriptome sequencing and de novo assembly
To understand the molecular basis of the difference in the reproductive performance between Con and PQQ, the placenta was used to build 6 libraries for high-throughput sequencing. We obtained a total of 54.18 Gb of raw data for the 6 samples (Table 4). We discarded low-quality reads, which contained adapters and unknown or low-quality bases, and after stringent quality checks and data cleaning, the clean reads were obtained (Table 4). The GC (guanine + cytosine) contents of these samples were 51.86–53.06%, with an average of 52.51%. The average Q20 and Q3 percentages reached 97.85% and 94.51%, respectively (Table 4).
Gene annotation and functional classification
All unigenes were aligned to 7 protein databases, including COG, GO, KEGG, KOG, Swiss-Prot, and Nr, using BLASTx with an E-value threshold of 10− 5 and Pfam using HMMER with an E-value threshold of 10− 10. As shown in Table 5, of 462 unigenes annotated, 448 (96.97%) unigenes had significant BLASTx matches in the Nr database. Based on comparison against the Swiss-Prot database, 370 (80.09%) unigenes had significant matches. In the Pfam and GO databases, 415 (89.93%) and 384 (83.12%) unigenes were also found to have significant matches, respectively, and 316 (68.40%) unigenes were similar to proteins in the KEGG database. To further evaluate the completeness of our transcriptome library and the effectiveness of our annotation process, we searched the annotated sequences for genes with COG (cluster of orthologous groups) classifications, and 145 unigenes were assigned to the COG classification (Fig. 7). Among the 25 COG categories, the cluster for “General function prediction only” (24, 15.19%) represented the largest group, followed by “Carbohydrate transport and metabolism” (20, 12.66%) and “Posttranslational modification, protein turnover, chaperones” (19, 12.03%). Gene ontology (GO) was also used to classify the functions of the predicted unigenes. Based on the sequence homology, 384 sequences were categorized into 61 functional groups (Fig. 6). The assigned functions of the unigenes covered a broad range of GO categories. The unigenes were assigned to three main categories, including the cellular component, molecular function, and biological process categories (Fig. 8). Reproduction (27, 7.03%) and reproduction process (27, 7.03%) were assigned to the biological process category.
Analysis of the differentially expressed unigenes (DEGs) by RNA-seq
The gene expression levels of DEGs in the placenta were measured using the calculated values of the RPKM parameter (reads per kilobase of exon per million reads mapped). The TMM (trimmed mean of M-values) method was used to standardize the read counts. Thereafter, the differentially expressed unigene analysis was performed using EB Seq. An absolute value of log2 fold change (FC) ≥2 and a false discovery rate (FDR) value of < 0.05 were employed to identify the DEGs and to explore the gene expression levels of DEGs in the placenta between the Con treatment and PQQ treatment. As shown in Additional file 1: Table S1, 462 unigenes were differentially expressed between the Con treatment and PQQ treatment. Among these unigenes, 199 were upregulated, while 263 unigenes were downregulated. By matching DEGs to the KEGG database, we identified 29 pathways (P < 0.05, Table 6).
Validation of differentially expressed genes by qRT-PCR
To validate the results of RNA-seq analysis, qRT-PCR was employed to determine the relative expression of 8 genes in the placenta (Fig. 9), including SOD1, IL6, IL8, NOS2, CDX2, CCN1, GCLC and CALM. IL6, IL8 and CCN1 were significantly downregulated. SOD1, NOS2, CDX2, GCLC and CALM were upregulated by PQQ·Na2 dietary supplementation (P < 0.05), which was consistent with the data from the RNA-seq analysis. Differences in the magnitude of fold-change values were likely due to differences in detection sensitivity of the two methods. In addition, correlation analysis demonstrated that the values of log2 (fold-change) obtained from RNA-seq and qRT-PCR were significantly correlated (R2 = 0.94). Thus, our RNA-seq analysis results are valid.
The use of PQQ in nutrition is increasingly being discussed in the literature [19, 20]. Jonscher KR et al.  reported that PQQ treatment led to a significant increase in placental weight and placental surface area. Yin et al.  suggested that the dietary supplementation of 1.50 mg/kg PQQ·Na2 is the lowest functional dose to improve the growth performance for weaned pigs and the expression of the jejunal tight junction protein ZO-1 was significantly higher in pigs with PQQ·Na2 supplementation. Zhang et al.  suggested that dietary 1.6 mg/kg PQQ·Na2 supplementation during gestation and lactation of female rats can significantly increase the number of viable foetuses per litter, born alive litter weight and the mRNA expression levels of CAT, GPX2 and SOD1 in the placenta. Steinberg et al.  observed that 0.4 mg PQQ/kg in an amino acid-based diet optimized reproduction. To our knowledge, this is the first study to examine the effects of dietary PQQ·Na2 supplementation during gestation and lactation in sows. Several reports have investigated DEGs from reproductive tissues, such as the endometrium, in pigs . DEGs in the placenta have rarely been demonstrated using RNA-seq.
During pregnancy and lactation in sows, high energy and oxygen levels are required to satisfy increasing metabolic burdens for fetal growth, placenta development and milk production, which could cause elevated ROS production  Oxidative stress results from increased production of reactive oxygen species (ROS) or a decrease in antioxidant defense. Oxidative damage is a strong indicator of the health status and wellbeing of animals . A recent study showed that pregnant sows had elevated oxidative stress during late gestation and lactation, which was responsible for impaired milk production, reproductive performance, and longevity of sows . The protective role of an antioxidant against oxidative stress in sows has been clearly demonstrated , and dietary addition of antioxidants can reduce oxidative stress and improve the reproductive performance of sows . Thus, dietary antioxidant concentrations need to be added in sow diets, especially to prevent excessive oxidative stress during gestation and lactation. In this study, we found that the antioxidant status of sows, including in the plasma, placenta and milk, was partially improved and oxidative stress markers were partially reduced by dietary PQQ·Na2 supplementation. SOD is known to serve a protective function for the elimination of reactive free radicals and, therefore, it represents an important antioxidant defense in nearly all cells exposed to oxygen. GSH forms an important part of the nonenzymatic antioxidants . Similar to other sulfhydryl-containing products, GSH also has regulatory and protective roles in the body, as it establishes the defenses of the body against tissue injury due to chemicals through its ROS scavenging, cell viability and membrane-stabilizing effects . In this study, dietary PQQ·Na2 supplementation can increase the antioxidant enzyme activities of SOD and GSH-Px in the placenta, plasma and milk. A previous study has demonstrated that SOD activity can be upregulated by PQQ·Na2 treatment , and this is consistent with our study results. The antioxidant enzymes SOD and CAT are considered the first line of defense against ROS. In this study, CAT was increased in the plasma (90 d of gestation). MDA is an end product of free-radical chain reactions and lipid peroxidation , so it is frequently used in the measurement of lipid peroxide levels, and it correlates well with the degree of lipid peroxidation . In the present study, the MDA levels in the placenta and plasma (21 d of lactation) was significantly decreased by PQQ·Na2 supplementation. PQQ acts as an antioxidant by scavenging O2− and protects mitochondria from oxidative stress-induced damage. In our study, the total piglets born, number of piglets born alive and born alive litter weight were significantly increased. Evidence has shown that supplementation with antioxidants, such as selenium, vitamin E and vitamin C, improves antioxidant status and reproductive performance in sows [35, 36]. Similarly, a growing number of studies have demonstrated that certain functional substances, such as chitosan , resveratrol , and isoflavone , alleviate oxidative stress and improve the reproductive performance of sows. Our results are in alignment with these studies and demonstrate that dietary PQQ·Na2 exerts a beneficial role in antioxidant defense and the reproductive performance of sows.
Colostrum is of great importance for the growth and development of piglets during and after lactation . Because neither creep feed nor milk replacer was used for the suckling piglets in this study, the sow milk served asthe sole source of nutrients and antioxidants for the piglets. Therefore, litter performance reflected the nutrient composition and antioxidant status of the colostrum and sow milk only. There is evidence that sow’s colostrum and milk include various antioxidants, including SOD, GSH-Px and GSH [25, 40]. In the present study, SOD and GSH-Px activity in the colostrum was significantly increased by dietary PQQ·Na2 during gestation and lactation. Milk antioxidants, which provide antioxidant protection to suckling piglets, may be as important as nutrients or immunological factors in protecting the health of the neonatal piglet . Colostrum is the main external resource providing piglets with nutrients and maternal immune molecules. Newborn piglets can hardly obtain passive immunity from the maternal blood during the fetal period because of the special epitheliochorial structure of the pig placenta. Before their own immune system is fully developed, colostrum is the sole external resource which provides piglets with nutrients and maternal immune molecules . In the present study, the concentrations of protein and total solids were increased. The high total solids and protein content in colostrum is mostly due to immunoglobulin . Immunoglobulin in colostrum, mainly IgG, provides humoral immune protection for the newborn piglet until its own immune system has sufficiently matured to respond to antigens . Additionally, immunoglobulin in the colostrum and milk could increase susceptibility to infection in newborn animals, not only in the immediate postnatal period but also after weaning . IgA is the main immunoglobulin of milk and could protect piglets against local pathogens, commensal bacteria and food antigens in the digestive tract . In our study, the concentrations of IgG, IgA and IgM in the colostrum were significantly increased by PQQ·Na2 supplementation, which might aid the humoral and mucosal immunity of piglets.
Nitric oxide (NO) is a pleiotropic regulator and is pivotal to numerous biological processes, including vasodilation, neurotransmission, and macrophage-mediated immunity . NO plays an important role in regulating placental-fetal blood flow, contributing to maternal systemic vasodilatation during pregnancy and regulating uterine and fetal placental blood flow . It has been reported that antioxidant activity inhibits lipid peroxidation, increases nitric oxide (NO) production, reduces oxidation of low-density lipoproteins, and preserves superoxide dismutase (SOD) activity . In the present study, PQQ increased antioxidants and iNOS in the placenta and plasma, which can increase NO. The increased NO plays an important part in improving placental vascular function and promoting the nutrient supply to the fetus. NO plays an important role in vasodilatation and regulates uterine blood flow, promoting the transfer of nutrients to the fetus . It has reported that NO is a key regulator of angiogenesis and embryogenesis as well as placental and fetal growth . In addition, NO might positively influence the final litter size, as it also induces follicle rupture in vitro in rabbit ovaries .
PQQ has attracted considerable attention, as it is important for mammalian growth, development, reproduction and immune function. However, the mechanism by which PQQ is beneficial to reproduction is not clear. Thus, we performed an RNA-seq analysis on the placenta of sows. We found 462 DEGs, including 199 upregulated genes and 263 downregulated genes. To validate the DEGs identified by RNA-seq analysis, we confirmed the expression levels of 8 genes by qRT-PCR. Comparison of the results obtained using the two methods revealed similar trends, confirming their validity.
The placenta is an important endocrine organ throughout pregnancy and markedly affects fetal health by supplying nutrients. Among the 25 COG categories, we found that the cluster for “General function prediction only” (24, 15.19%) represented the largest group, followed by “Carbohydrate transport and metabolism” (20, 12.66%). Effects of inadequate dietary protein and carbohydrate intake on maternal venous plasma amino acid profiles occurred from early pregnancy on, probably with impacts on placental growth and function and eventually on fetal development . By matching DEGs to the KEGG database, we identified 29 pathways. From these pathways, the “Steroid biosynthesis”, “TGF-beta signaling pathway”, “Jak-STAT signaling pathway”, “Insulin signaling pathway”, “Hippo signaling pathway” and “Folate biosynthesis” are closely related to reproductive performance. Ovarian and placental steroids are essential for the maintenance of pregnancy. In some mammals it is evident that the placenta is responsible for the production of steroids . The TGF-beta superfamily has been observed to oppose trophoblastic migration, suggesting an intricate balance of growth factor expression that is fundamental for placental health and embryo development . Tojo reported that the entire TGF-β signal transduction pathway is essential to proper development and embryonic survival . The JAK-STAT pathway is frequently activated and is indispensable and pivotal in many biological processes, including immunity and the inflammatory response . In previous studies we showed that PQQ regulates intracellular JAK/STAT signaling pathway activation. Eskouhie Tchaparian et al. reported changes in gene expression patterns and transcriptional networks that respond to dietary PQQ restriction or pharmacological administration, they found JAK/STAT pathways seem particularly influenced by PQQ . Zhang reported that the insulin signaling pathway was down-regulated in the placenta of women with gestational diabetes mellitus . Dysregulation of Hippo signaling component genes can result in embryonic lethality . Folates are needed for fetal growth and placental development, since they activate cell growth and biosynthetic processes that are essential during pregnancy .
GO analyses showed that 27 DEGs were involved in reproduction. Among these DEGs, the ovo-like zinc finger 2 (OVOL2), distal-less homeobox 1 (DLX1), distal-less homeobox 2 (DLX2), distal-less homeobox 5 (DLX5), msh homeobox 2 (MSX2), caudal type homeobox 2 (CDX2), nitric oxide synthase 2 (NOS2) and calmodulin (CALM) were upregulated. Unezaki reported that ovol2 function is required for endothelial cell growth during heart development and angiogenesis of extraembryonic and embryonic vessels . The Dlx genes encode a family of transcription factors with important roles in patterning and differentiation during vertebrate embryogenesis . Previous studies indicate that MSX2 plays an important role in mammalian embryonic diapause [63, 64]. Sakurai suggested that CDX2 is essential for early development and gene expression and is involved in differentiation of the inner cell mass and trophectoderm lineages in embryos . Kwon showed, by comparing sperm proteins from different litter sizes, that CALM was highly expressed in high swine litter sizes and was positively related to litter size . Furthermore, antioxidant genes, such as SOD1, GCLC and DHCR24, were also upregulated by PQQ·Na2 supplementation. In sows, supplementation improves antioxidant levels and alleviates oxidative stress effectively, which are beneficial to litter size and piglet growth. DHCR24 exerts cytoprotective effects against endoplasmic reticulum stress by eliminating ROS. DHCR24 can scavenge hydrogen peroxide (H2O2), protecting cells from oxidative stress-induced apoptosis . GCLC is an important part of GSH, which is an extremely important antioxidant. It not only scavenges free radicals but also maintains the redox-sensitive active sites of many enzymes from an oxidized form to a reduced form . The supply of glucose to the embryo from the maternal circulation is important for normal metabolism and growth, as glucose constitutes the main energy substrate during embryogenesis . The solute carrier (SLC) family (including SLC1A1, SLC7A4, SLC7A10 and SLC19A1) was upregulated, and many nutrient carriers and growth factors decrease as the dam is exposed to stress. Oxidative stress, defined as an imbalance between the production of free radicals and reactive metabolites, is closely related to inflammation. GO analysis identified genes involved in the inflammatory response (including IL6, IL8, IL11 and CCN1) that were downregulated. KEEG pathway analysis showed that the cytokine-cytokine receptor interaction pathway was significantly reduced, which indicated that the inflammatory state of the placenta was alleviated by maternal PQQ·Na2 supplementation. A previous study showed that PQQ can reduce the expression of inflammatory cytokine genes  which was consistent with our study.
In conclusion, our results have shown that dietary 20 mg/kg PQQ·Na2 supplementation during gestation and lactation in sows can significantly increase the total piglets born, number of piglets born alive and born alive litter weight. It also increased antioxidant levels in the placenta, plasma and milk. The concentration of NO was significantly increased in the plasma and placenta. RNA-seq analysis showed that 462 unigenes were differentially expressed between Con treatment and PQQ treatment. Among these unigenes, 199 were upregulated, while 263 unigenes were downregulated. By matching DEGs to the KEGG database, we identified 29 pathways. These provide a theoretical basis to further explore the effect of PQQ on the reproductive performance mechanism of sows. The present study can provide a scientific basis for dietary PQQ·Na2 supplementation in sows.
Availability of data and materials
The datasets produced and/or analyzed during the current study are available from the corresponding author on reasonable request.
Analysis of variance
Cellular communication network factor 1
Caudal type homeobox 2
Caudal type homeobox 2
Clusters of Orthologous Groups
Differentially expressed genes
Distal-less homeobox 1
Distal-less homeobox 2
Distal-less homeobox 5
False discovery rate
Fragments per kilobase of transcript per million fragments mapped
Glutamate-cysteine ligase catalytic subunit
Kyoto Encyclopedia of Genes and Genomes
msh homeobox 2
National Center for Biotechnology Information
Nitric oxide synthase 2
Ovo like zinc finger 2
- PQQ·Na2 :
Pyrroloquinoline quinone disodium
Quantitative real-time PCR
Reactive oxygen species
Standard error of the mean
Solute carrier family 19 member 1
Solute carrier family 1 member 1
Solute carrier family 7 member 10
Solute carrier family 7 member 4
Zhang S, Yang J, Wang L, Li Z, Pang P, Li F. SLA-11 mutations are associated with litter size traits in large white and Chinese DIV pigs. Anim Biotechnol. 2018. https://doi.org/10.1080/10495398.2018.1471401.
Balcells I, Castelló A, Mercadé A, et al. Analysis of porcine MUC4 gene as a candidate gene for prolificacy QTL on SSC13 in an Iberian × Meishan F2 ovulation. BMC Genetics. 2011;12(1):93.
Huang L, Yin ZJ, Feng YF, et al. Identification and differential expression of microRNAs in the ovaries of pigs (\r, Sus scrofa\r, ) with high and low litter sizes. Anim Genet. 2016. https://doi.org/10.1111/age.12452.
Ren P, Yang XJ, Railton R, et al. Effects of different levels of feed intake during four short periods of gestation and housing systems on sows and litter performance. Anim Reprod Sci. https://doi.org/10.1016/j.anireprosci.2017.11.001.
Leddy MA, Power ML, Schulkin J. The impact of maternal obesity on maternal and fetal health. Rev Obstet Gynecol. 2008;1(4):170–8.
Tarrade A, Panchenko P, Junien C, et al. Placental contribution to nutritional programming of health and diseases: epigenetics and sexual dimorphism. J Exp Biol. 2015;218(1):50–8.
Lee D , Nam J , Kim S W , et al. Proteomic analysis of reproduction proteins involved in litter size from porcine placenta. Journal of the Agricultural Chemical Society of Japan, 2015. doi: https://doi.org/10.1080/09168451.2015.1039478.
Zhang Y, Feustel PJ, Kimelberg HK. Neuroprotection by pyrroloquinoline quinone (PQQ) in reversible middle cerebral artery occlusion in the adult rat. Brain Res. 2006;1094(1):200–6.
Killgore J, Smidt C, Duich L, et al. Nutritional importance of pyrroloquinoline quinone. Science. 1989;245(4920):850–2.
Akagawa M , Nakano M , Ikemoto K . Recent progress in studies on the health benefits of pyrroloquinoline quinone. Journal of the Agricultural Chemical Society of Japan, 2016, 80(1):10.
Ikemoto K, Mori S, Mukai K. Synthesis and crystal structure of pyrroloquinoline quinol (PQQH2) and pyrroloquinoline quinone (PQQ). Acta Crystallogr Sect B: Struct Sci. 2017;73(3):489–97.
Hwang P, Willoughby DS. Mechanisms behind Pyrroloquinoline Quinone supplementation on skeletal muscle mitochondrial biogenesis: possible synergistic effects with exercise. J Am Coll Nutr. 2018:1–11.
Mitchell AE, Jones AD, Mercer RS, et al. Characterization of Pyrroloquinoline Quinone amino acid derivatives by electrospray ionization mass spectrometry and detection in human Milk. Anal Biochem. 1999;269(2):317–25.
Samuel KG, Zhang HJ, Wang J, et al. Effects of dietary pyrroloquinoline quinone disodium on growth performance, carcass yield and antioxidant status of broiler chicks. Animal. 2014;9(3):8.
O’Loughlin A, Lynn DJ, Mark MG, et al. Transcriptomic analysis of the stress response to weaning at housing in bovine leukocytes using RNA-seq technology. BMC Genomics. 2012;13(1):250.
Gao W, Sun HX, Xiao H, et al. Combining metabolomics and transcriptomics to characterize tanshinone biosynthesis inSalvia miltiorrhiza. BMC Genomics. 2014;15(1):73.
Noji N, Nakamura T, Kitahata N, et al. Simple and sensitive method for Pyrroloquinoline Quinone (PQQ) analysis in various foods using liquid chromatography/electrospray-ionization tandem mass spectrometry. J Agric Food Chem. 2007;55(18):7258–63.
Grabherr MG, Haas BJ, Yassour M, et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011;29(7):644–52.
Wang J, Zhang HJ, Xu L, et al. Dietary supplementation of pyrroloquinoline quinone disodium protects against oxidative stress and liver damage in laying hens fed an oxidized sunflower oil-added diet. Animal. 2016;10(07):1129–36.
Harris CB, Chowanadisai W, Mishchuk DO, et al. Dietary pyrroloquinoline quinone (PQQ) alters indicators of inflammation and mitochondrial-related metabolism in human subjects[J]. J Nutr Biochem. 2013;24(12):2076–84.
Jonscher KR, Stewart MS, Alfonso-Garcia A, et al. Early PQQ supplementation has persistent long-term protective effects on developmental programming of hepatic lipotoxicity and inflammation in obese mice. FASEB J. 2017;31(4):1434–48.
Yin X, Ming D, Bai L, et al. Effects of pyrroloquinoline quinone supplementation on growth performance and small intestine characteristics in weaned pigs. J Anim Sci. 2018. https://doi.org/10.1093/jas/sky387.
Zhang B, Yang W, Zhang H, He S, Meng Q, Chen Z, Shan A. Effect of pyrroloquinoline quinone disodium in female rats during gestating and lactating on reproductive performance and the intestinal barrier functions in the progeny. Br. J. Nutr. 2019;121(7). https://doi.org/10.1017/S0007114519000047.
Steinberg FM. Dietary pyrroloquinoline quinone : growhth and immune response in BALB/c mice. J Nutr. 1994;124:744.
Chen X, Li A, Chen W, et al. Differential gene expression in uterine endometrium during implantation in pigs. Biol Reprod. 2015;92(2):52–2.
Hu YJ, Gao KG, Zheng CT, et al. Effect of dietary supplementation with glycitein during late pregnancy and lactation on antioxidative indices and performance of primiparous sows. J Anim Sci. 2015;93(5):2246.
Agarwal A, Makker K, Sharma R. Clinical relevance of oxidative stress in male factor infertility: an update. Am J Reprod Immunol. 2010;59(1):2–11.
Lapointe J. Mitochondria as promising targets for nutritional interventions aiming to improve performance and longevity of sows. J Anim Physiol Anim Nutr. 2014;98(5):809–21.
Agarwal A, Aponte-Mellado A, Premkumar BJ, et al. The effects of oxidative stress on female reproduction: a review. Reprod Biol Endocrinol. 2012;10(1):49.
Cheung CCC, Zheng GJ, Li AMY, et al. Relationships between tissue concentrations of polycyclic aromatic hydrocarbons and antioxidative responses of marine mussels, Perna viridis. Aquatic Toxicology (Amsterdam). 2001;52(3–4):189–203.
Szabó S. Role of sulfhydryls and early vascular lesion in gastric mucosal injury. Acta Physiol Hung. 1984;64(3–4):203–14.
Zhang Q, Ding M, Gao XR, et al. Pyrroloquinoline quinone rescues hippocampal neurons from glutamate-induced cell death through activation of Nrf2 and up-regulation of antioxidant genes. Genet Mol Res. 2012;11(3):2652–64.
Tomita M, Okuyama T, Kawai S. Determination of malonaldehyde in oxidized biological materials by high-performance liquid chromatography. J Chromatogr A. 1990;515(5):391–7.
Wei B, Nie S, Meng Q, et al. Effects of l-carnitine and/or maize distillers dried grains with solubles in diets of gestating and lactating sows on the intestinal barrier functions of their offspring. Br J Nutr. 2016;116(03):459–69.
Chen J, Han JH, Guan WT, et al. Selenium and vitamin E in sow diets: I. effect on antioxidant status and reproductive performance in multiparous sows. Anim Feed Sci Tech. 2016;221:111–23.
Pinellisaavedra A. Vitamin E in immunity and reproductive performance in pigs. Reprod Nutr Dev. 2003;43(5):397.
Xie C, Wu X, Long C, et al. Chitosan oligosaccharide affects antioxidant defense capacity and placental amino acids transport of sows. BMC Vet Res. 2016;12(1):243.
Qingwei, Meng, Gaoqiang, et al. Dietary resveratrol improves antioxidant status of sows and piglets and regulates antioxidant gene expression in placenta by Keap1-Nrf2 pathway and Sirt1. Journal of livestock and biotechnology, 2018(3):639–651.
Hurley WL, Theil PK. Perspectives on immunoglobulins in colostrum and milk[J]. Nutrients. 2011;3(4):442.
Lipkoprzybylska J, Kankofer M. Antioxidant defence of colostrum and milk in consecutive lactations in sows[J]. Ir Vet J. 2012;65(1):4-4.
Tan L, Wei T, Yuan A, et al. Dietary supplementation of Astragalus polysaccharides enhanced immune components and growth factors EGF and IGF-1 in sow colostrum. J Immunol Res. 2017;2017(10):1–6.
Wang LS, Su BC, Shi Z, et al. Dietary supplementation with maize distillers dried grains with solubles during late gestation and lactation: effects on sow and litter performance, and on colostrum and milk composition. Anim Feed Sci Technol. 2013;179(1–4):149–53.
Liu ST, Hou WX, Cheng SY, et al. Effects of dietary citric acid on performance, digestibility of calcium and phosphorus, milk composition and immunoglobulin in sows during late gestation and lactation. Anim Feed Sci Technol. 2014;191:67–75.
Varley MA, Towle AMA. Artificial rearing of piglets: the administration of two sources of immunoglobulins after birth. Anim Prod. 2010;43(1):121–6.
Brandtzaeg P. The Mucosal Immune System and Its Integration with the Mammary Glands. J Pediatr. 2010;156(2 suppl):S8–15.
Rao CV. Nitric oxide signaling in colon cancer chemoprevention. Mutat Res. 2004;555(1–2):107–19.
Sladek SM, Magness RR, Conrad KP. Nitric oxide and pregnancy. Am J Phys. 1997;272(2 Pt 2):R441–63.
Mason RP. Mechanisms of plaque stabilization for the dihydropyridine calcium channel blocker amlodipine: review of the evidence. Atherosclerosis. 2002;165(2):191–9.
Moncada S, Higgs EA. Molecular mechanisms and therapeutic strategies related to nitric-oxide. FASEB J. 1995;9(13):1319.
Wu G, Bazer FW, Cudd TA, et al. Maternal nutrition and fetal development. Early Hum Dev. 1980;4(2):99–120.
Yamauchi J, Miyazaki T, Iwasaki S, et al. Effects of nitric oxide on ovulation and ovarian steroidogenesis and prostaglandin production in the rabbit. Endocrinology. 1997;138(9):3630–7.
Strakovsky RS, Zhou D, Pan YX. A low-protein diet during gestation in rats activates the placental mammalian amino acid response pathway and programs the growth capacity of offspring. J Nutr. 2010;140(12):2116–20.
Braun BC, Zschockelt L, Dehnhard M, et al. Progesterone and estradiol in cat placenta--biosynthesis and tissue concentration. J Steroid Biochem Mol Biol. 2012;132(3–5):295–302.
Nikuei P , Malekzadeh K , Rajaei M , et al. The imbalance in expression of angiogenic and anti-angiogenic factors as candidate predictive biomarker in preeclampsia. 2015.
Tojo M, Takebe A, Takahashi S, et al. Smad7-deficient mice show growth retardation with reduced viability. J Biochem. 2012;151(6):621–31.
O’Shea JJ, Schwartz DM, Villarino AV, et al. The JAK-STAT pathway: impact on human disease and therapeutic intervention*. Annu Rev Med. 2015;66(1):311–28.
Tchaparian E, Marshal L, Cutler G, et al. Identification of transcriptional networks responding to pyrroloquinoline quinone dietary supplementation and their influence on thioredoxin expression, and the JAK/STAT and MAPK pathways. Biochem J. 2010;429(3):515–26.
Zhang W, Su R, Lin L, et al. ARHGEF11 affecting the placental insulin signaling pathway in fetal macrosomia of normal glucose tolerance pregnant women. Placenta. . https://doi.org/10.1016/j.placenta.2017.12.010.
Halder G, Johnson RL. Hippo signaling: growth control and beyond. Development. 2011;138(1):9–22.
Castaño E, Caviedes L, Hirsch S, et al. Folate Transporters in Placentas from Preterm Newborns and Their Relation to Cord Blood Folate and Vitamin B12 Levels. Plos One. 2017;12(1):e0170389.
Unezaki S, Horai R, Sudo K, Iwakura Y, Ito S. Ovol2/Movo, a homologue of Drosophila ovo, is required for angiogenesis, heart formation and placental development in mice. Genes Cells. 2007;12(6):13.
Coubrough ML, Bendall AJ. Impaired nuclear import of mammalian Dlx4 proteins as a consequence of rapid sequence divergence. Exp Cell Res. 2006;312(19):3880–91.
Cha J, Sun X, Bartos A. A new role for muscle segment homeobox genes in mammalian embryonic diapause. Open Biol. 2013;3(4):130035.
Renfree MB, Shaw G. Embryo-endometrial interactions during early development after embryonic diapause in the marsupial tammar wallaby. Int J Dev Biol. 2014;58(2–4):175–81.
Sakurai N , Takahashi K , Emura N , et al. The Necessity of OCT-4 and CDX2 for Early Development and Gene Expression Involved in Differentiation of Inner Cell Mass and Trophectoderm Lineages in Bovine Embryos. Cellular Reprogramming, 2016:cell.2015.0081.
Kwon WS, Rahman MS, Lee JS, et al. Discovery of predictive biomarkers for litter size in boar spermatozoa. Mol Cell Proteomics. 2015;14(5):1230.
Lu X, Kambe F, Cao X, et al. 3beta-Hydroxysteroid-delta24 reductase is a hydrogen peroxide scavenger, protecting cells from oxidative stress-induced apoptosis. Endocrinology. 2008;149(7):3267.
Jian Z, Mei P, Ki K, et al. Fucoxanthin enhances the level of reduced glutathione via the Nrf2-mediated pathway in human keratinocytes. Marine Drugs. 2014;12(7):4214–30.
Jensen VFH, AM ML, Lykkesfeldt J, et al. Effect of maternal hypoglycaemia during gestation on materno-foetal nutrient transfer and embryo-foetal development: Evidence from experimental studies focused primarily on the rat. Reprod Toxicol. 2018. https://doi.org/10.1016/j.reprotox.2018.01.007.
The authors would like to thank the Mudanjiang Dawan Animal Husbandry Company with limited liability for animal management.
This work was supported by the National Key Research and Development Plan of China (2016YFD0501207), the China Agriculture Research System (CARS-36).
Ethics approval and consent to participate
The protocols used in this experiment were approved by the Northeast Agricultural University Institutional Animal Care and Use Committee. All animal experimental procedures were approved by the Ethical and Animal.
Welfare Committee of Heilongjiang Province, China.
Consent for publication
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Table S1. Analysis of the differentially expressed unigenes (DEGs) by RNA-seq. (DOC 644 kb)