Skip to main content

Differentially expressed genes in cotyledon of ewes fed mycotoxins

Abstract

Background

Ergot alkaloids (E+) are mycotoxins produced by the endophytic fungus, Epichloë coenophiala, in tall fescue that are associated with ergotism in animals. Exposure to ergot alkaloids during gestation reduces fetal weight and placental mass in sheep. These reductions are related to vasoconstrictive effects of ergot alkaloids and potential alterations in nutrient transport to the fetus. Cotyledon samples were obtained from eight ewes that were fed E+ (n = 4; E+/E+) or E- (endophyte-free without ergot alkaloids; n = 4; E−/E-) seed during both mid (d 35 to 85) and late (d 85–133) gestation to assess differentially expressed genes associated with ergot alkaloid induced reductions in placental mass and fetal weight, and discover potential adaptive mechanisms to alter nutrient supply to fetus.

Results

Ewes fed E+/E+ fescue seed during both mid and late gestation had 20% reduction in fetal body weight and 33% reduction in cotyledon mass compared to controls (E−/E-). Over 13,000 genes were identified with 110 upregulated and 33 downregulated. Four genes had a |log2FC| > 5 for ewes consuming E+/E+ treatment compared to controls: LECT2, SLC22A9, APOC3, and MBL2. REViGO revealed clusters of upregulated genes associated glucose, carbohydrates, lipid, protein, macromolecular and cellular metabolism, regulation of wound healing and response to starvation. For downregulated genes, no clusters were present, but all enriched GO terms were associated with anion and monocarboxylic acid transport. The complement and coagulation cascade and the peroxisome proliferator-activated receptor signaling pathway were found to be enriched for ewes consuming E+/E+ treatment.

Conclusions

Consumption of ergot alkaloids during gestation altered the cotyledonary transcriptome specifically related to macronutrient metabolism, wound healing and starvation. These results show that ergot alkaloid exposure upregulates genes involved in nutrient metabolism to supply the fetus with additional substrates in attempts to rescue fetal growth.

Background

Tall fescue [Lolium arundinaceum (Schreb.) Darbysh; Schedonorus phoenix (Scop.) Holub] is the primary cool season perennial grass utilized in the eastern U.S. occupying more than 14 million ha [1]. The majority of tall fescue contains an endophyte (Epichloë coenophiala), which produces ergot alkaloids (i.e. ergovaline, ergovalinine, etc. [2];, mycotoxins associated with ergotism in humans and animals. The endophyte is beneficial to the plant and aids in establishment, persistence, and drought tolerance [1]; however, ingestion of ergot alkaloids by grazing livestock results in fescue toxicosis, a syndrome which reduces animal growth [3] and reproductive performance [4]. Previous research shows that exposure to ergovaline and ergovalinine during late (d 86 to parturition) gestation alters placental development [5, 6] and vasoactivity of umbilical arteries [7], and results in asymmetrical growth and intrauterine growth restriction of the fetus [8]; however, little is known about the mechanisms associated with these changes in the placenta.

Intrauterine growth restriction (IUGR) is the inability of a fetus to reach its expected growth potential in-utero, which leads to 10-fold higher rates of morbidity and mortality, the primary cause of stillbirth in humans [9,10,11]. IUGR is characterized by asymmetrical growth in which brain development is spared at the expense of overall growth [12]. The most common cause of IUGR is placental insufficiency where reductions in placental mass or impaired placental transport capacity lead to an inefficient transfer of substrates to the growing fetus [13]. In most mammalian species, uterine and umbilical blood flows are directly correlated with placental and fetal weights [14, 15] and reductions in uterine blood flow can lead to stunted fetal growth and an increase in offspring morbidity and mortality [16, 17]. Conditions that negatively influence fetal growth, such as maternal undernutrition or increased fetal number, are often associated with reductions in nutrient uptake, blood flow, and placental angiogenesis, the formation of new vascular beds [18, 19]. Angiogenesis and vasodilation of the uterine and placental vessels are primary mechanisms that work to increase blood flow during gestation. Furthermore, establishment of efficient fetal and maternal vascular beds within the placenta is crucial to support the rapid increase in placental blood flow, especially during the later stages of gestation when the majority of fetal growth occurs [14, 15].

Placental insufficiency and subsequent IUGR can be induced in gestating sheep through the consumption of ergovaline, an ergot alkaloid found in endophyte-infected tall fescue known to cause systematic vasoconstriction [6, 8]. Ergovaline binds to 5-hydroxytryptamine receptor 2A (5HT2a) and α-2 adrenergic receptors located throughout the cardiovascular system. These receptors are also located in the umbilical and uterine arteries of pregnant sheep and cows [20,21,22]. Vasoconstrictive activity associated with ergot alkaloid exposure has also been reported in the uterine vessels of heifers, uterine arteries of pregnant sheep, and umbilical arteries in sheep [7, 20, 23]. The resulting vasoconstriction reduces uteroplacental blood flow and limits fetal growth [6]. Sheep are utilized as a model for placental insufficiency and IUGR in humans due to similarities in placental structure and parallel IUGR pathologies after birth and into adulthood [24,25,26,27]. Proper placental vascularization and function depends on the finely tuned regulation of gene expression throughout gestation [28]. Aberrant gene expression in the placenta indicates altered molecular pathways which can provide insight into the causes of diseases such as preeclampsia and IUGR. This makes RNA-sequencing useful in detecting altered pathways associated with diseases states during pregnancy [29, 30]. It was hypothesized that reduction in placental and fetal mass was likely due to the vasoconstrictive effects exerted on uterine and placental blood flow that are commonly associated with ergot alkaloid exposure [20]. Several known angiogenic and vascular growth factors were investigated in the cotyledon tissues of ergovaline-fed ewes, but no differences were found that could account for the overall reduction in placental mass [6]. Therefore, the objective of this study was to assess differentially expressed genes using RNA-Seq associated with ergot alkaloid induced reductions in placentome and fetal weight and discover potential adaptive mechanisms to alter nutrient supply to fetus.

Results

Placentomes

Total placentome weight was 30% lower (P < 0.05) for E+/E+ ewes compared to E−/E-; however, total placentome number and caruncle weight did not differ due to ergot alkaloid treatment (Table 1). Exposure to ergot alkaloids during gestation reduced (P < 0.05) cotyledon weight by 34% compared to E−/E-. Fetal weight was 20% lower (P < 0.05) for E+/E+ ewes compared to E−/E-. Ratios of placental efficiency did not differ between ergot alkaloid treatments.

Table 1 Placental and fetal weights at d 133 of gestation for ewes fed endophyte-infected, ergot alkaloid containing (E+/E+) or endophyte-free, no ergot alkaloids (E−/E-) tall fescue seed during mid and late gestation

PCoA and RNAseq data

A Principle Coordinates Analysis (PCoA) was conducted to visualize variation in the samples (Fig. 1). Principle component 1 (PC1) accounted for 29.39% of the variation in the samples and separated the samples based on fescue treatment. Principle component 2 (PC2) accounted for 20.22% of the variation. PC2 accounts for variation in the control samples whereas the E+/E+ samples remain closely associated with one another. In total, 13,572 genes were identified through mapping to the ovine genome across all samples (Table 2). Out of these, 110 genes were upregulated (FDR < 0.1) with 15 having a Log2FC between 0 and 1. The remaining 95 upregulated genes had a Log2FC ≥ 1. A total of 33 genes that were downregulated with 17 having a Log2FC between 0 and − 1. The remaining 16 downregulated genes had a Log2FC ≤ − 1. A volcano plot for all DEGs detected in cotyledon tissues for ewes on E+/E+ fescue treatment compared to E−/E- is presented in Fig. 1. Black dots represent an FDR ≥ 0.10 or |log2FC| < 1 while colored dots represent genes determined to be differentially expressed with a |log2FC| ≥ 1 and an FDR < 0.10.

Fig. 1
figure1

Principle Coordinates Analysis (PCoA) plot for RNAseq results in cotyledon tissues for ewes on E+/E+ fescue treatment compared to E−/E- at d 133 of gestation (a) and volcano plot for all DEGs detected in cotyledon tissues for ewes on E+/E+ fescue treatment compared to E−/E- at d 133 of gestation (b). Black dots represent an FDR ≥ 0.10 or |log2FC| < 1. Colored dots represent genes determined to be differentially expressed with a |log2FC| > 1 and an FDR < 0.10. Differences in color represent varying degrees of FC

Table 2 Summary of reads mapped to the ovine genome

Genes of interest

Genes of interest with an FDR < 0.01 and |log2FC| > 3 are presented in Table 3. Four genes had a |log2FC| > 5 in the cotyledon for ewes consuming E+/E+ treatment compared to E−/E- treatment. These included leukocyte cell derived chemotaxin 2 (LECT2), solute carrier family 22 member 9 (SLC22A9), apolipoprotein C3 (APOC3), and mannose binding lectin 2 (MBL2). In some cases, multiple genes within a given gene family were identified as being upregulated in the cotyledon tissue of E+/E+ treated ewes compared to E−/E- treatment. Five apolipoproteins (APOs) were upregulated including APOC3, APOC4, APOA1, APOB, and APOH. Several serpin (SERPIN) genes were also identified including SERPINA3–8, SERPINF2, SERPINA1, and SERPIND1. Inter-alpha-trypsin inhibitor heavy chains 1–3 (ITIH1, ITIH2, ITIH3) and the cytochrome P450 gene, CYP2E1, were also found to be upregulated in E+/E+ treatment compared to the control.

Table 3 Candidate genes from cotyledon tissues for ewes on E+/E+ fescue treatment compared to E−/E- at d 133 of gestation. Only genes with a |log2FC| > 3 and a FDR < 0.01 are displayed

The five genes assessed using qPCR were found to be similar in their calculated Log2FC compared to RNA sequencing results (Fig. 2). The expression of LECT2 and CYP2E1 was found to be up-regulated (P < 0.05) in E+/E+ compared to E−/E-. Additionally, STRA6 expression was down-regulated (P < 0.05) in E+/E+ compared to E−/E-. There was no statistical difference based on treatment for the two other genes tested, AGT and FABP1.

Fig. 2
figure2

Differential expression of genes measured by qPCR in cotyledon of E+/E+ ewes compared to E−/E- (n = 4/treatment)

Figure 3 shows the average LogCPM normalized to GAPDH for genes with high levels of expression that have not previously been identified in ovine cotyledon samples. These genes included TIMP metallopeptidase inhibitor 2 and 3 (TIMP2, TIMP3), carboxypeptidase X, M14 family member 1 (CPXM1), calpain 6 (CAPN6), and placentally expressed transcript 1 (PLET1).

Fig. 3
figure3

The average LogCPM normalized to GAPDH across all cotyledon samples regardless of treatment including vascular and growth factors and genes with high levels of expression that have not previously been identified in ovine cotyledon samples

Gene ontology and KEGG pathway analysis

REViGO generates a cluster representation of enriched biological process GO terms among upregulated genes. The REViGO analyses can be seen in Fig. 4 for upregulated and downregulated genes. A closer proximity denotes a closer relationship between the terms and an increase in size indicates fold-change. For upregulated genes, the most tightly associated cluster included responses to glucose, carbohydrates, and xenobiotic stimuli. There were two GO term clusters associated with lipid and protein metabolism. GO terms representing the response to starvation, regulation of wound healing, and liver development were also present but fairly independent of other terms. For downregulated genes, no clusters were present, but all enriched GO terms were associated with transport. This included anion, monocarboxylic and organic acid, and drug transmembrane transport. There was an enrichment for the regulation of peptide breakdown in E+/E+ cotyledon tissue compared to the control tissues. There was a fold enrichment of 15 or more for several types of peptidase and endopeptidase activity which included both regulation and inhibition of activity. The molecular functions of protease binding, carboxylic ester hydrolase activity, and carbohydrate derivative binding were also highly enriched in the gene set. As shown Table 4, KEGG pathway analysis identified two distinct pathways that were upregulated in E+/E+ treated samples. The complement and coagulation cascade was highly enriched and included plasminogen activator tissue type (PLAT), plasminogen (PLG), paraoxonase 3 (PON3), mannose binding lectin 2 (MBL2), complements C9, C6, and C1S, fibrinogen alpha chain (FGA), fibrinogen gamma chain (FGG), and serpins SERPINF2, SERPINC1, SERPINA1 and SERPIND1. The peroxisome proliferator-activated receptor (PPAR) signaling pathway was also enriched and included apolipoprotein A2 (APOA2), apolipoprotein A5 (APOA5), apolipoprotein C3 (APOC3), enoyl-CoA hydratase and 3-hydroxyacyl CoA dehydrogenase (EHHADH), fatty acid binding protein 1 (FABP1), adiponectin (ADIPOQ), and phosphoenolpyruvate carboxykinase 1 (PCK1).

Fig. 4
figure4

A scatterplot generated with REViGO shows a cluster representation of enriched biological process Gene Ontology (GO) term among upregulated genes (a) and downregulated genes (b)

Table 4 KEGG pathway enrichment analysis using David (v 6.8) showed enrichment in the complement and coagulation cascade and PPAR signaling pathway for ewes on E+/E+ fescue compared to E−/E- treatment

Discussion

Exposure to ergot alkaloids during mid and late gestation reduced placentome and fetal weights at near-term, d 133 of an average 145 d gestation length. The reductions in the placentome weight were a direct result of reductions in the cotyledon, the fetal side of the placentome, mass and not of the caruncle, the maternal side. Estimates of placental sufficiency showed no difference in the fetal to cotyledon weight ratio, which suggests that the reduction in cotyledon weight limited growth of the fetus. During the first half of gestation, the ovine placenta experiences rapid proliferative growth and peaks in weight around d 80 after which a period of remodeling begins where connective tissue in the core of the fetal villi is replaced by vascular beds [31, 32]. This angiogenic vascular development in the cotyledon results in a 12-fold increase in capillary density, a decrease in capillary size, and continues throughout late gestation to support rapid fetal growth [33]. Previous research documents the extensive list of vascular and growth factors necessary during this process throughout mid and late gestation [15, 32,33,34,35] and differences in gene expression have been reported in cases of undernutrition [36], hyperthermia [37], and hypoxia [38]. Ergot alkaloids are known vasoconstrictors [39, 40] that bind to 5HT2A serotonergic receptors. Previous examinations of known vascular factors involved in placental growth did not reveal any differences due to ergot alkaloid exposure [6] and therefore this RNA-sequencing project was developed to identify differentially expressed genes associated with ergot alkaloid exposure and discover potential adaptive mechanisms that alter nutrient supply to fetus.

The RNA-Seq results identified 95 genes that were up-regulated and 16 down-regulated in response to feeding mycotoxins. Interestingly, LECT2 was the most upregulated gene with a 7-fold increase in the cotyledon samples of E+/E+ treated ewes. It is described as an energy-sensing hepatokine that is preferentially expressed in human adult and fetal liver cells but is secreted into the bloodstream [41]. Several studies have associated LECT2 with glucose metabolism, insulin resistance and inflammation [42,43,44]. Knockout mouse studies show that loss of LECT2 increases insulin signaling in skeletal muscle [41]. To our knowledge, LECT2 expression has not been previously identified in placental tissues nor associated with IUGR. Therefore, further research is warranted to investigate the role of LECT2 in cases of mycotoxin exposure and irregular placental development.

Ewes on E+/E+ fescue seed treatment during gestation experienced an upregulation in genes associated with responses to carbohydrates, glucose and carbohydrate biosynthesis. The main carbohydrate utilized by the placenta is glucose which is transported by GLUT transporters through the process of facilitated diffusion [45]. Under normal conditions the majority of glucose is derived from maternal circulation, but glucose can be derived from fetal blood if maternal concentrations are lacking [46]. Because of dependence on the concentration gradient, changes in the placental uptake and consumption of glucose is directly related to changes in maternal concentrations [47]. The upregulation in genes associated with responses to carbohydrates and glucose may indicate an increase in maternal concentrations of glucose reaching the placenta. However, plasma glucose concentrations were not different based on fescue seed treatment in the associated study [6]. During late gestation, placental consumption of glucose is up to 10 times higher than the fetus for the production of ATP and other sugars and carbohydrates including polyols [48]. The polyol pathways are highly active in the early placenta and are closely associated with the pentose phosphate pathway, which supports rapid cell proliferation [49, 50]. In sheep, additional carbohydrates such as lactate and fructose are produced and shuttled to umbilical and uterine circulations. Sheep with growth restricted placentae exhibit a reduction in overall glucose consumption on a placental weight basis, but an increase in the conversion of glucose (or other substrates) to lactate, much of which is shuttled to the fetus [51]. A similar scenario is seen in humans in which up to 22% of glucose is converted to lactate under normoxic conditions [52]. Because lactate is unable to be used by the placenta, it is speculated that this may be an attempt to set aside resources for the fetus [49]. The increase in carbohydrate biosynthesis suggests a similar scenario in which glucose or other substrates are converted to lactate or fructose at an increased rate in ewes on E+/E+ fescue seed treatment.

Triglyceride and cellular lipid metabolism were also upregulated in E+/E+ cotyledon samples compared to controls. Lipids and free fatty acids are required for growth and development of both the fetus and placenta but transport across the placenta appears minimal for most livestock species [53]. In sheep, the placenta may hydrolyze esterified lipids and then desaturate and/or elongate them to supply essential fatty acids to the fetus [53]. These processes, along with the synthesis of lipids from glucose and ketoacids, may provide necessary fats for both placental and fetal tissues during gestation [45]. Research also shows that volatile fatty acids (VFAs) produced through the process of rumen fermentation are utilized by the placenta and transported to the fetus in both cows and sheep [53]. It is hypothesized that the placental utilization of VFAs may work to generate ATP or synthesize fatty acids and, in adverse conditions, lipids could be utilized as an alternate energy source if insufficient glucose is available [48, 54]. In cases of maternal undernutrition in sheep, there is an increase in placental fatty acid transporters [55] and placental lipid metabolism is known to be altered in cases of IUGR [56]. Additionally, KEGG pathway analysis determined that PPAR signaling was upregulated in the cotyledon of E+/E+ treated ewes. Recent research has elucidated critical roles for PPARs in placental development and the pathophysiology of IUGR. PPARs are ligand-activated transcription factors that regulate gene expression in a variety of tissues and play a critical role in placental lipid metabolism [57]. Human placentae were found to have increased PPAR expression levels in cases of preeclampsia and IUGR when compared to controls and it has been suggested this may be an adaptive response to compensate for insufficient placental development [58]. APOC3 is found within the PPAR pathway and had a 5-fold greater expression in cotyledon samples from E+/E+ treated ewes. APOC3 is associated with hypertriglyceridemia and elevates plasma triglyceride levels by preventing clearance of very-low-density lipoproteins (VLDLs) and high-density lipoproteins (HDLs) but limited research is available denoting its presence or role in placental function.

A primary function of the placenta is to facilitate nutrient exchange to the fetus through the use of transport systems. Transporters are in place to mediate the transfer of endogenous compounds across the placental barrier but, in some cases, transporters accept specific exogenous substrates such as drugs or xenobiotics if they are present [59]. Such transporters are identified as drug transporters and include monocarboxylate transporters (MCTs) and anion transporters [59]. MCTs generally work to transport lactate, much of which is produced by placental metabolism when glucose is in short supply [60]. The placenta also uses anion transporters, such as organic anion transporters (OATs), which have been discovered in almost all barrier epithelia within the body [61]. Both MCTs and OATs also function to transport a wide range of drugs and toxins in addition to their endogenous substrates [59]. Ewes on E+/E+ fescue seed treatment experienced a consistent downregulation in a variety of transport systems, including monocarboxylic acid and anion transport, within the cotyledon compared to ewes on the control treatment. In contrast, solute carrier family 22 member 9 (SLC22A9; also known as OAT7) experienced a 5-fold increase in the cotyledon from E+/E+ ewes and was previously thought to be liver-specific [62]. It has been hypothesized that ergot alkaloids or metabolites may cross the placental barrier [7] but no research is available documenting the presence of these compounds in the fetus or which transport systems may be involved if placental transport exists. It was previously proposed that the downregulation of certain transport systems may augment fetal exposure to environmental toxicants during pregnancy [63]. A similar mechanism may be at work to limit fetal exposure to ergot alkaloids if placental transport occurs.

During late gestation, the primary function of the placenta is to facilitate nutrient exchange to the fetus in order to keep pace with exponential fetal growth [33]. REViGO revealed several upregulated GO clusters associated with macromolecule metabolism, including protein and lipid metabolism, in the cotyledon of ewes exposed to E+ fescue seed during gestation compared to those ewes on E- treatment. In addition to facilitating nutrient exchange to the fetus, the placenta is a metabolically active organ with specific energy and nutrient requirements. The placenta also utilizes, produces, and converts amino acids throughout gestation [53]. The placental protein turnover rate was previously reported at 60% per day in sheep [64]. Such high protein turnover rates are heavily dependent on the amino-acid availability in maternal circulation. During the initial period of placental growth in humans, the placenta is the primary recipient of maternally circulating amino acids. As fetal growth accelerates during late gestation, the fetus becomes the primary beneficiary of amino acids [65]. In sheep, the placenta is a net consumer of glutamate, serine, valine, leucine and isoleucine and shuttles greater concentrations of glutamine, methionine and glycine towards the fetus compared to levels found in maternal uterine circulation suggesting significant utilization and/or conversion of amino acids [66]. The placenta experiences a high level of protein synthesis during late gestation due to the changes in placental morphology. It was hypothesized that this period would result in a high rate of protein turnover [45, 67]. Nutritional stress in sheep has been shown to alter fetoplacental amino acid metabolism and cycling and in some cases changes in metabolism appear permanent [68, 69]. The upregulation in protein metabolism and proteolysis in the cotyledon of E+/E+ fescue treated ewes may indicate increased placental remodeling or a higher rate of protein turnover.

The KEGG pathway analysis determined the complement and coagulation cascade pathway was highly upregulated in the cotyledon of E+/E+ treated ewes. While the complement and coagulation networks are distinct, several key crossover points keep them linked and both must be tightly regulated during pregnancy. Increased thrombin production, which has been associated with preeclampsia and IUGR, may be a result of increased activation of the coagulation cascade beyond that of normal pregnancy [70]. Additionally, several coagulation components are necessary for vascular differentiation which is often impaired in the placentae of preeclamptic women [70]. Regulation of the innate and adaptive immune response through the complement system is necessary for successful placental and fetal development [71]. While some degree of activation is required, early embryonic loss and IUGR are often associated with increased complement activation [72]. The complement system has three routes of activation: the classical pathway, alternate pathway, and mannose-binding lectin pathway. The cotyledon samples of ewes exposed to E+/E+ fescue seed treatment had a 5-fold upregulation of MBL. Mannose binding lectin functions as a critical part of the innate immune system and works as a first line of defense against microorganisms [72]. Under normal physiological conditions, MBL does not recognize an organism’s own tissues. However, in cases of cellular hypoxia, altered cell surface glycosylation can stimulate MBL expression and activate the complement system [73]. In this study, ergot alkaloid induced vasoconstriction may result in placental hypoxia that stimulates MBL and the complement cascade.

Several genes presented with a high rate of expression that have not previously been identified in ovine cotyledonary tissues. TIMP2 and TIMP3 are members of the TIMP metallopeptidase inhibitor family that are involved in degradation of extracellular matrix and suppress endothelial cell proliferation. CPXM1 is a member of the carboxypeptidase X, M14 family and placental expression has not been reported previously but it has been shown to be a positive regulator of adipogenesis [74]. CAPN6 is a member of the calpain family which are calcium-dependent cysteine proteases shown to be highly expressed in the placenta where they are involved in tumorigenesis by promoting angiogenesis [75]. Placentally expressed transcript 1 (PLET1) has been characterized in the pig where it was shown to be involved in trophoblastic elongation of the conceptus [76]. While these genes did not show differences based on treatment, their high level of expression not previously reported in sheep was worth noting.

The research presented here was our initial exploration to identify the genes associated with reduced cotyledon mass observed when feeding mycotoxins containing ergovaline and ergovalinine. Using RNA-Seq methods, we identified several genes (LECT2, SLC22A9, APOC3, MBL2) that were highly (log2FC > 5; FDR < 0.01) up-regulated in response to mycotoxin exposure and one (STRA6) that was down-regulated (log2FC > 3; FDR < 0.01). Comparisons between RNA-Seq and qPCR showed 86% agreement in fold-change but genes with lower fold-change were not significant for qPCR. These comparisons are similar to those reported by Everaert and co-workers [77] who reported 85% agreement between RNA-Seq and qPCR comparisons of gene expression fold-change but genes with lower expression had inconsistent expression measurements. We acknowledge that this study had a very small sample size (n = 3 biological replicates per treatment) and that this is a limitation of the study. These results identified novel genes that may be involved in mycotoxin exposure and further research with larger numbers is warranted to examine how the cotyledon tissue responds to ergot alkaloid induced vasoconstriction.

Conclusion

The present study demonstrates that the cotyledon transcriptome was altered in response to ergot alkaloid exposure during mid and late gestation compared to control ewes. Upregulated genes were in clusters for glucose, carbohydrates, lipid, protein, macromolecular and cellular metabolism, regulation of wound healing and response to starvation. For downregulated genes, no clusters were present, but all enriched GO terms were associated with anion and monocarboxylic acid transport. The complement and coagulation cascade and the peroxisome proliferator-activated receptor signaling pathway were found to be enriched for ewes consuming E+/E+ treatment. These results show that ergot alkaloid exposure upregulates genes involved in nutrient metabolism to supply the fetus with additional substrates in attempts to rescue fetal growth.

Methods

Experimental design

Placental samples utilized in this study were collected from a larger study and detailed experimental design information is available [6]. Mature Suffolk ewes (n = 57; 4.8 yr of age; 82 kg avg. BW), naïve to endophyte-infected tall fescue, were purchased from a private farm and transported to Clemson University Small Ruminant Facility 90 d prior to the start of the experiment. Ewes were synchronized using an intravaginal controlled internal drug release (CIDR) insert (Eazi-Breed CIDR, Zoetis Animal Health) for 7 d. Upon CIDR removal, ewes were given prostaglandin F2α (12.5 mg i.m.; Lutalyse; Pfizer, New York, NY) and turned in with a purebred Suffolk ram. The ram was fitted with a marking harness and crayon that was changed weekly. Ewes were checked twice daily for marks to estimate breeding date and confirmed pregnant by transrectal ultrasonography on d 30 of gestation. Pregnant Suffolk ewes (n = 32; n = 8/treatment) were randomly assigned to E+ (endophyte-infected; 1.77 mg/head/day of ergovaline and ergovalinine) or E- (endophyte-free; 0 mg/head/day of ergovaline and ergovalinine) tall fescue seed during MID (d 35–85) and/or LATE (d 86–133) gestation which provided four dietary treatments: E−/E-, E+/E-, E−/E+, and E+/E+. Ewes were individually penned into stalls (1.8 × 0.5 × 0.91 m) after 0700 and individually fed their respective treatment diet for 90 min. After individual feeding, ewes were removed from stalls and kept in a large pen (10–12 hd/pen) with ad libitum access to water and minerals (Purina Sheep Mineral, Land O’Lakes Inc., Arden Hills, MN), and with access to inside and outside areas devoid of forage or hay. For this study, a subsample (n = 4/treatment) was selected from the E−/E- and E+/E+ treatments that had twin fetuses and were representative of treatment means for further examination of cotyledon by RNA-seq. Ewe was the experimental unit for this study.

Sample collection

On d 133 of gestation, ewes were transported to Godley-Snell Research Facility (14.3 km) at 0730 for terminal surgery. Each ewe was given an intravenous injection of Ketamine (10 mg/kg) and Diazepam (0.25 mg/kg) upon arrival for sedation and intubation. Ewes were intubated with a 10 mm endotracheal tube and placed on 4–5% isoflurane with 1–2 L/min of O2 for induction. Upon successful anesthetization, ewes were maintained at 3–5% isoflurane with 1–2 L/min of O2 and placed on a ventilator at 15–20 breaths per minute. Corneal reflex and heart rate were monitored to confirm adequate levels of anesthesia. The abdominal area was shaved, scrubbed with chlorhexidine, and ewes were subjected to a mid-ventral laparotomy. Once the uterus was exposed, an incision was made in the uterine wall to collect fetuses and a placentome of the type B morphology [78] was selected adjacent to the initial incision. Placentomes were immediately separated into the cotyledon and caruncle portions, and the cotyledon was flash frozen in liquid nitrogen and stored at − 80 °C for subsequent RNA extraction. Once removed, fetal lambs were euthanized with a 3 mL intracardiac injection of Beuthanasia-D Special (Merck Animal Health, Madison, New Jersey). Each fetus was towel dried and fetal weight was collected. Once all fetuses were removed, ewes were euthanized with an intravenous injection of 20 mL Beuthanasia-D Special. Placental and fetal data were analyzed using ANOVA with fescue treatment in the model.

Sample preparation and RNA sequencing analysis

Total cellular RNA was extracted using Trizol reagent and purified using the PureLink Mini RNA kit (Invitrogen, Thermo Fisher Scientific, Waltham, MA) according to manufacturer instructions. RNA quantity was estimated using a NanoDrop 1000 Spectrophotometer (ThermoFisher). RNA integrity was determined using an Agilent 2100 Bioanalyzer (Agilent, Santa Clara, CA). All samples submitted for RNAseq had a RIN value of 9.3 or greater. RNA samples were shipped on dry ice to LCSciences (Houston, TX) for mRNA sequencing. Library preparation was completed using the Illumina Small RNA Sample Preparation Kit (Illumina, San Diego, CA USA). Samples were sequenced to an average sequencing depth of 10 M reads using an Illuminia HiSeq 2500 sequencing platform at 50 bp single end reads. Downsampling to ensure uniform coverage between samples was done using seqtk v1.3-r106 (https://github.com/lh3/seqtk). FastQC v0.11.7 was used to check quality of the samples [79]. Any low-quality bases and adapter sequences were trimmed using Trimmomatic v0.38 [80]. After trimming, 84.8% of the total trimmed reads across all samples were aligned to the ovine reference genome (ftp://ftp.ensembl.org/pub/release-96/fasta/ovis_aries/) using GSNAP v2018-07-04 [81]. Subread’s feature Counts v1.6.2 software was utilized to count uniquely mapped reads that aligned to known genes in the reference genome [82]. Raw gene counts were input to edgeR v3.22.5 and genes with at least 4 samples having at least 1 count per million (cpm) were kept and these counts were normalized [83]. A principal coordinates analysis plot was constructed showing one sample from each treatment to be an outlier. Sequencing and mapping data can be seen in Table 1. These samples were removed leaving n = 3/treatment as seen in Fig. 1. Raw counts from the subset samples were rerun through edgeR, filtering out genes with less than 3 samples having a cpm of less than 1, and used to calculate differentially expressed genes (DEG) in E−/E- vs E+/E+. A volcano plot was generated from the data. Genes were considered differentially expressed if their false discovery rate (FDR) was < 0.05 and their |log2FC| was greater than or equal to 1.

Gene ontology (GO) and KEGG pathway analysis

Gene ontology (GO) enrichment analysis was conducted using Panther 14.1 (http://pantherdb.org/). DEGs with a |log2FC| greater than or equal to 1 were used to analyze upregulated and downregulated gene sets separately. Bos taurus was utilized as a reference in the absence of the sheep genome. Panther GO-Slim was used to assess biological process, molecular function, and cellular component. Enriched GO terms were semantically clustered and visualized using REViGO. Upregulated and downregulated genes were analyzed separately. David (v 6.8) was used to determine KEGG pathway enrichment.

RT-qPCR validation

Gene expression analysis was conducted using quantitative real-time RT-PCR methods according to Duckett and co-workers [84] to serve as a comparison for RNAseq results. Aliquots from the same RNA samples submitted for RNAseq were converted to cDNA using qScript cDNA SuperMix (Quanta Bio, Beverly, MA) according to manufacturer instructions. Real-time PCR was performed using a QuantStudio3 (Applied Biosystems, Thermo Fisher) and Perfecta (Quanta Bio, Beverly, MA) SYBR green according to the manufacturer’s directions. An initial hold of 2 min at 95 °C was followed by 40 cycles of 95 °C for 15 s and 60 °C for 30 s. Primers for leukocyte cell derived chemotaxin 2 (LECT2), fatty acid binding protein 1 (FABP1), angiotensinogen (AGT), cytochrome p450 2E1 (CYP2E1) and stimulated by retinoic acid 6 (STRA6) were designed using PrimerQuest (Integrated DNA Technologies, IDT, Coralville, IA) and used to compare to RNA-Seq. Glyceraldehyde 3-phosphate dehydrogenase (GAPDH), βactin (bACT), and cyclophilin (CYC) were tested for stability using RefFinder for the selection of the most stable housekeeping gene [85]. The most stable housekeeping gene, GAPDH, was used for normalization. Normalized CT values (ΔCT = CT, gene − CT, GAPDH) were calculated for each sample. The relative quantification was the ratio of the target gene to internal control genes using the ΔΔCT method. Normalized CT values were analyzed using an ANOVA with fescue treatment in the model. Genes in cotyledon tissue from ewes fed E+/E+ fescue were considered differentially expressed relative to control at a cut of P < 0.05 for qPCR.

Availability of data and materials

RNA-Seq data are published and freely available in Gene Expression Omnibus GEO, accession number GSE155975, https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE155975.

Abbreviations

5HT2a:

5-hydroxytryptamine receptor 2A

ADIPOQ:

adiponectin

APO:

apolipoprotein

APOA2:

apolipoprotein A2

APOA5:

apolipoprotein A5

APOC3:

apolipoprotein C3

APOC4:

apolipoprotein C4

APOA1:

apolipoprotein A1

APOB:

apolipoprotein B

APOH:

apolipoprotein H

ATP:

Adenosine triphosphate

BW:

body weight

C1S:

Complement C1S

C6:

Complement C6

C9:

Complement C9

CPM:

counts per million

CYP3A4:

cytochrome P450 3A4

CYP2E1:

cytochrome P450 2E1

DEG:

differentially expressed genes

EHHADH:

Enoyl-CoA Hydratase And 3-Hydroxyacyl CoA Dehydrogenase

FABP1:

Fatty acid binding protein 1

FC:

Fold change

FDR:

false discovery rate

FGA:

Fibrinogen Alpha Chain

FGG:

Fibrinogen Gamma Chain

GAPDH:

Glyceraldehyde 3-phosphate dehydrogenase

GO:

gene ontology

IUGR:

Intrauterine growth restriction

ITIH:

Inter-alpha-trypsin inhibitor heavy chains

LECT2:

leukocyte cell derived chemotaxin 2

mRNA:

messanger RNA

MBL2:

mannose binding lectin 2

MCT:

monocarboxylate transporters

PON3:

paraoxonase 3

OAT:

organic anion transporters

PCK:

Phosphoenolpyruvate Carboxykinase 1

PCoA:

Principle Coordinates Analysis

PLAT:

Plasminogen Activator Tissue Type

PLG:

Plasminogen

PPAR:

peroxisome proliferator-activated receptor

RIN:

RNA integrity number

RNA:

ribonucleic acid

RNAseq:

RNA sequencing

SERPINA1:

Serpin Family A Member 1

SERPINA3–8:

Serpin Family A Member 3–8

SERPINC1:

Serpin Family C Member 1

SERPIND1:

Serpin Family D Member 1

SERPINF2:

Serpin Family F Member 2

1SLC22A9:

solute carrier family 22 member 9

TMR:

total mixed ration

References

  1. 1.

    Stuedemann JA, Hoveland CS. Fescue endophyte: history and impact on animal agriculture. J Prod Agric. 1988;1:39–44.

    Article  Google Scholar 

  2. 2.

    Young CA, Charlton ND, Takach JE, Swoboda G, Trammell MA, Huhman DV, Hopkins AA. Characterization of Epichloe coenophiala within the US: are all tall fescue endophytes created equal? Front Chem. 2014;2:95. https://doi.org/10.3389/fchem.2014.00095.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  3. 3.

    Hoveland CS, Schmidt SP, King CC Jr, Odom JW, Clark EM, McGuire JA, et al. Steer performance and association of Acremonium coenophialum fungal endophyte on tall fescue pastures. Agron J. 1983;75:821–4.

    Article  Google Scholar 

  4. 4.

    Peters CW, Grigsby KN, Aldrich CG, Paterson JA, Lipsey RJ, Kerley MS, Garner GB. Performance, forage utilization, and ergovaline consumption by beef cows grazing endophyte fungus-infected tall fescue, endophyte fungus-free tall fescue, and orchardgrass pastures. J Anim Sci. 1992;70:1550–61.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  5. 5.

    Duckett SK, Andrae JG, Pratt SL. Exposure to ergot alkaloids during gestation reduces fetal growth in sheep. Front Chem. 2014;2:68. https://doi.org/10.3389/fchem.2014.00068.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  6. 6.

    Britt JL, Greene MA, Bridges WC, Klotz JL, Aiken GE, Andrae JG, et al. Ergot alkaloid exposure during gestation alters. I. Maternal characteristics and placental development of pregnant ewes1. J Anim Sci. 2019;97(4):1874–90.

    PubMed  PubMed Central  Article  Google Scholar 

  7. 7.

    Klotz JL, Britt JL, Miller MF, Snider MA, Aiken GE, Long NM, et al. Ergot alkaloid exposure during gestation alters: II. Uterine and umbilical artery vasoactivity1. J Anim Sci. 2019;97(4):1891–902.

    PubMed  PubMed Central  Article  Google Scholar 

  8. 8.

    Greene MA, Britt JL, Powell RR, Feltus FA, Bridges WC, Bruce T, et al. Ergot alkaloid exposure during gestation alters: 3. Fetal Growth, Muscle Fiber Development and miRNA Transcriptome. J Anim Sci. 2019;97(4):1874–90. https://doi.org/10.1093/jas/skz068.

  9. 9.

    Sharma D, Sharma P, Shastri S. Genetic, metabolic and endocrine aspect of intrauterine growth restriction: an update. J Matern Fetal Neonatal Med. 2017;30(19):2263–75.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  10. 10.

    Froen JF, Gardosi JO, Thurmann A, Francis A, Stray-Pedersen B. Restricted fetal growth in sudden intrauterine unexplained death. Acta Obstet Gynecol Scand. 2004;83(9):801–7.

    PubMed  Article  PubMed Central  Google Scholar 

  11. 11.

    McIntire DD, Bloom SL, Casey BM, Leveno KJ. Birth weight in relation to morbidity and mortality among newborn infants. N Engl J Med. 1999;340(16):1234–8.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  12. 12.

    Cohen E, Baerts W, van Bel F. Brain-sparing in intrauterine growth restriction: considerations for the neonatologist. Neonatology. 2015;108(4):269–76.

    PubMed  Article  PubMed Central  Google Scholar 

  13. 13.

    Morrison JL. Sheep models of intrauterine growth restriction: fetal adaptations and consequences. Clin Exp Pharmacol Physiol. 2008;35(7):730–43.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  14. 14.

    Reynolds LP, Borowicz PP, Vonnahme KA, Johnson ML, Grazul-Bilska AT, Wallace JM, et al. Animal models of placental angiogenesis. Placenta. 2005;26(10):689–708.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  15. 15.

    Reynolds LP, Borowicz PP, Vonnahme KA, Johnson ML, Grazul-Bilska AT, Redmer DA, et al. Placental angiogenesis in sheep models of compromised pregnancy. J Physiol. 2005;565(Pt 1):43–58.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  16. 16.

    Lang U, Baker RS, Braems G, Zygmunt M, Kunzel W, Clark KE. Uterine blood flow--a determinant of fetal growth. Eur J Obstet Gynecol Reprod Biol. 2003;110(Suppl 1):S55–61.

    PubMed  Article  PubMed Central  Google Scholar 

  17. 17.

    Marconi AM, Ronzoni S, Vailati S, Bozzetti P, Morabito A, Battaglia FC. Neonatal morbidity and mortality in intrauterine growth restricted (IUGR) pregnancies is predicated upon prenatal diagnosis of clinical severity. Reprod Sci. 2009;16(4):373–9.

    PubMed  Article  PubMed Central  Google Scholar 

  18. 18.

    Reynolds LP, Caton JS, Redmer DA, Grazul-Bilska AT, Vonnahme KA, Borowicz PP, et al. Evidence for altered placental blood flow and vascularity in compromised pregnancies. J Physiol. 2006;572(Pt 1):51–8.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  19. 19.

    Wallace JM, Regnault TR, Limesand SW, Hay WW Jr, Anthony RV. Investigating the causes of low birth weight in contrasting ovine paradigms. J Physiol. 2005;565(Pt 1):19–26.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  20. 20.

    Dyer DC. Evidence that ergovaline acts on serotonin receptors. Life Sci. 1993;53(14):PL223–8.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  21. 21.

    Klotz JL, Brown KR, Xue Y, Matthews JC, Boling JA, Burris WR, et al. Alterations in serotonin receptor-induced contractility of bovine lateral saphenous vein in cattle grazing endophyte-infected tall fescue. J Anim Sci. 2012;90(2):682–93.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  22. 22.

    Klotz JL, Aiken GE, Johnson JM, Brown KR, Bush LP, Strickland JR. Antagonism of lateral saphenous vein serotonin receptors from steers grazing endophyte-free, wild-type, or novel endophyte-infected tall fescue. J Anim Sci. 2013;91(9):4492–500.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  23. 23.

    Poole DH, Lyons SE, Poole RK, Poore MH. Ergot alkaloids induce vasoconstriction of bovine uterine and ovarian blood vessels. J Anim Sci. 2018;96(11):4812–22. https://doi.org/10.1093/jas/sky328.

  24. 24.

    Barry JS, Anthony RV. The pregnant sheep as a model for human pregnancy. Theriogenology. 2008;69(1):55–67.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  25. 25.

    Leiser R, Krebs C, Ebert B, Dantzer V. Placental vascular corrosion cast studies: a comparison between ruminants and humans. Microsc Res Tech. 1997;38(1–2):76–87.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  26. 26.

    Beede KA, Limesand SW, Petersen JL, Yates DT. Real supermodels wear wool: summarizing the impact of the pregnant sheep as an animal model for adaptive fetal programming. Anim Front. 2019;9(3):34–43.

    PubMed  PubMed Central  Article  Google Scholar 

  27. 27.

    Hay WW Jr, Brown LD, Rozance PJ, Wesolowski SR, Limesand SW. Challenges in nourishing the intrauterine growth-restricted foetus - lessons learned from studies in the intrauterine growth-restricted foetal sheep. Acta Paediatr. 2016;105(8):881–9.

    PubMed  PubMed Central  Article  Google Scholar 

  28. 28.

    Mikheev AM, Nabekura T, Kaddoumi A, Bammler TK, Govindarajan R, Hebert MF, et al. Profiling gene expression in human placentae of different gestational ages: an OPRU network and UW SCOR study. Reprod Sci. 2008;15(9):866–77.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  29. 29.

    Majewska M, Lipka A, Paukszto L, Jastrzebski JP, Myszczynski K, Gowkielewicz M, et al. Transcriptome profile of the human placenta. Funct Integr Genomics. 2017;17(5):551–63.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  30. 30.

    Cox B, Leavey K, Nosi U, Wong F, Kingdom J. Placental transcriptome in development and pathology: expression, function, and methods of analysis. Am J Obstet Gynecol. 2015;213(4 Suppl):S138–51.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  31. 31.

    Stegeman HJ. A study of the maturation of the placenta in sheep. Acta Morphol Neerl Scand. 1972;10(4):400.

    CAS  PubMed  PubMed Central  Google Scholar 

  32. 32.

    Borowicz PP, Arnold DR, Johnson ML, Grazul-Bilska AT, Redmer DA, Reynolds LP. Placental growth throughout the last two thirds of pregnancy in sheep: vascular development and angiogenic factor expression. Biol Reprod. 2007;76(2):259–67.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  33. 33.

    Reynolds LP, Borowicz PP, Caton JS, Vonnahme KA, Luther JS, Buchanan DS, et al. Uteroplacental vascular development and placental function: an update. Int J Dev Biol. 2010;54(2–3):355–66.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  34. 34.

    Carr DJ, David AL, Aitken RP, Milne JS, Borowicz PP, Wallace JM, et al. Placental vascularity and markers of angiogenesis in relation to prenatal growth status in overnourished adolescent ewes. Placenta. 2016;46:79–86.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  35. 35.

    Reynolds LP, Redmer DA. Angiogenesis in the placenta. Biol Reprod. 2001;64(4):1033–40.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  36. 36.

    Redmer DA, Milne JS, Aitken RP, Johnson ML, Borowicz PP, Reynolds LP, et al. Decreasing maternal nutrient intake during the final third of pregnancy in previously overnourished adolescent sheep: effects on maternal nutrient partitioning and feto-placental development. Placenta. 2012;33(2):114–21.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  37. 37.

    Regnault TR, Galan HL, Parker TA, Anthony RV. Placental development in normal and compromised pregnancies-- a review. Placenta. 2002;23 Suppl A:S119–29.

  38. 38.

    Zhang S, Barker P, Botting KJ, Roberts CT, McMillan CM, McMillen IC, et al. Early restriction of placental growth results in placental structural and gene expression changes in late gestation independent of fetal hypoxemia. Physiol Rep. 2016;4(23):e13049. https://doi.org/10.14814/phy2.13049.

  39. 39.

    Strickland JR, Looper ML, Matthews JC, Rosenkrans CF Jr, Flythe MD, Brown KR. Board-invited review: St. Anthony’s fire in livestock: causes, mechanisms, and potential solutions. J Anim Sci. 2011;89(5):1603–26. https://doi.org/10.2527/jas.2010-3478.

    CAS  Article  PubMed  Google Scholar 

  40. 40.

    Klotz JL. Activities and effects of ergot alkaloids on livestock physiology and production. Toxins (Basel). 2015;7(8):2801–21. https://doi.org/10.3390/toxins7082801.

    CAS  Article  Google Scholar 

  41. 41.

    Yamagoe S, Mizuno S, Suzuki K. Molecular cloning of human and bovine LECT2 having neutrophil chemotactic activity and its specific expression in the liver. Biochim Biophys Acta. 1998;1396(1):105–13. https://doi.org/10.1016/s0167-4781(97)00181-4.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  42. 42.

    Lan F, Misu H, Chikamoto K, Takayama H, Jujychi A, Nihri K, et al. LECT2 functions as a hepatokine that links obesity to skeletal muscle insulin resistance. Diabetes. 2014;63:1649–64. https://doi.org/10.2337/db13-07828.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  43. 43.

    Jung TW, Chung YH, Kim HC, El-Aty AMA, Jeong JH. LECT2 promotes inflammation and insulin resistance in adipocytes via P38 pathways. J Mol Endocrinol. 2018;61(1):37–45. https://doi.org/10.1530/JME-17-0267.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  44. 44.

    Yoo HJ, Hwang SY, Choi JH, Lee HJ, Chung HS, Seo JA, et al. Association of leukocyte cell-derived chemotaxin 2 (LECT2) with NAFLD, metabolic syndrome and atherosclerosis. PLoS One. 2017;12(4):e0174717. https://doi.org/10.1371/journal.pone.0174717.ecollection.

    Article  PubMed  PubMed Central  Google Scholar 

  45. 45.

    Bell AW, Ehrhardt RA. Regulation of placental nutrient transport and implications for fetal growth. Nutr Res Rev. 2002;15(2):211–30.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  46. 46.

    Simmons MA, Battaglia FC, Meschia G. Placental transfer of glucose. J Dev Physiol. 1979;1(3):227–43.

    CAS  PubMed  PubMed Central  Google Scholar 

  47. 47.

    Hay WW Jr, Molina RA, DiGiacomo JE, Meschia G. Model of placental glucose consumption and glucose transfer. Am J Phys. 1990;258(3 Pt 2):R569–77.

    CAS  Google Scholar 

  48. 48.

    Vaughan OR, Fowden AL. Placental metabolism: substrate requirements and the response to stress. Reprod Domest Anim. 2016;51(Suppl 2):25–35. https://doi.org/10.1111/rda.12797.

  49. 49.

    Jauniaux E, Hempstock J, Teng C, Battaglia FC, Burton GJ. Polyol concentrations in the fluid compartments of the human conceptus during the first trimester of pregnancy: maintenance of redox potential in a low oxygen environment. J Clin Endocrinol Metab. 2005;90(2):1171–5.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  50. 50.

    Burton GJ, Fowden AL. The placenta: a multifaceted, transient organ. Philos Trans R Soc Lond Ser B Biol Sci. 2015;370(1663):20140066.

    Article  Google Scholar 

  51. 51.

    Owens JA, Falconer J, Robinson JS. Effect of restriction of placental growth on fetal and utero-placental metabolism. J Dev Physiol. 1987;9(3):225–38.

    CAS  PubMed  PubMed Central  Google Scholar 

  52. 52.

    Schneider H. Placental oxygen consumption. Part II: in vitro studies--a review. Placenta. 2000;21(Suppl A):S38–44.

    PubMed  Article  PubMed Central  Google Scholar 

  53. 53.

    Vaughan OR, Fowden AL. Placental metabolism: substrate requirements and the response to stress. Reprod D0mest Anim. 2016;51(Suppl 2):25–35.

    CAS  Article  Google Scholar 

  54. 54.

    Christie WW, Noble RC. Fatty acid biosynthesis in sheep placenta and maternal and fetal adipose tissue. Biol Neonate. 1982;42(1–2):79–86.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  55. 55.

    Ma Y, Zhu MJ, Uthlaut AB, Nijland MJ, Nathanielsz PW, Hess BW, et al. Upregulation of growth signaling and nutrient transporters in cotyledons of early to mid-gestational nutrient restricted ewes. Placenta. 2011;32(3):255–63.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  56. 56.

    Cetin I, Alvino G. Intrauterine growth restriction: implications for placental metabolism and transport. A review. Placenta. 2009;30(Suppl A):S77–82.

    PubMed  Article  CAS  PubMed Central  Google Scholar 

  57. 57.

    Xu Y, Wang Q, Cook TJ, Knipp GT. Effect of placental fatty acid metabolism and regulation by peroxisome proliferator activated receptor on pregnancy and fetal outcomes. J Pharm Sci. 2007;96(10):2582–606.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  58. 58.

    Holdsworth-Carson SJ, Lim R, Mitton A, Whitehead C, Rice GE, Permezel M, et al. Peroxisome proliferator-activated receptors are altered in pathologies of the human placenta: gestational diabetes mellitus, intrauterine growth restriction and preeclampsia. Placenta. 2010;31(3):222–9.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  59. 59.

    Ganapathy V, Prasad PD. Role of transporters in placental transfer of drugs. Toxicol Appl Pharmacol. 2005;207(2 Suppl):381–7.

    PubMed  Article  CAS  PubMed Central  Google Scholar 

  60. 60.

    Nagai A, Takebe K, Nio-Kobayashi J, Takahashi-Iwanaga H, Iwanaga T. Cellular expression of the monocarboxylate transporter (MCT) family in the placenta of mice. Placenta. 2010;31(2):126–33.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  61. 61.

    Nigam SK, Bush KT, Martovetsky G, Ahn SY, Liu HC, Richard E, et al. The organic anion transporter (OAT) family: a systems biology perspective. Physiol Rev. 2015;95(1):83–123.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  62. 62.

    Jungst C, Klein K, Eloranta JJ, Kullak-Ublick GA. Bile acids Downregulate the human hepatic organic anion transporter 7. J Hepatol. 2012;56:S141–S2.

    Article  Google Scholar 

  63. 63.

    Fontes KN, Reginatto MW, Silva NL, et al. Dysregulation of placental ABC transporters in a murine model of malaria-induced preterm labor. Sci Rep. 2019;9:11488 https://doi.org/10.1038/s41598-019-47865-3.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  64. 64.

    Young M. Protein turnover rate in early life. Acta Paediatr Acad Sci Hung. 1982;23(1):99–117.

    CAS  PubMed  PubMed Central  Google Scholar 

  65. 65.

    Chien PF. Investigations of protein metabolism in human pregnancy: the term foetus and placenta studied using stable isotope labelled amino-acids. Clin Nutr. 1991;10(Suppl):70–6.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  66. 66.

    Chung M, Teng C, Timmerman M, Meschia G, Battaglia FC. Production and utilization of amino acids by ovine placenta in vivo. Am J Phys. 1998;274(1):E13–22.

    CAS  Article  Google Scholar 

  67. 67.

    Vaughan OR, Rosario FJ, Powell TL, Jansson T. Regulation of placental amino acid transport and fetal growth. In: Huckle WR, editor. Molecular biology of placental development and disease 2017. p. 217–21.

  68. 68.

    Kwon H, Ford SP, Bazer FW, Spencer TE, Nathanielsz PW, Nijland MJ, et al. Maternal nutrient restriction reduces concentrations of amino acids and polyamines in ovine maternal and fetal plasma and fetal fluids. Biol Reprod. 2004;71(3):901–8.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  69. 69.

    Liechty EA, Kelley J, Lemons JA. Effect of fasting on uteroplacental amino acid metabolism in the pregnant sheep. Biol Neonate. 1991;60(3–4):207–14.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  70. 70.

    Klaitman V, Beer-Wiesel R, Rafaeli T, Mazor M, Erez O. The role of the coagulation system in preterm parturition. 2013. In: Perterm Birth [Internet]. IntechOpen. Available from: https://www.intechopen.com/books/preterm-birth/the-role-of-the-coagulation-system-in-preterm-parturition.

  71. 71.

    Regal JF, Gilbert JS, Burwick RM. The complement system and adverse pregnancy outcomes. Mol Immunol. 2015;67(1):56–70.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  72. 72.

    Takahashi K, Ezekowitz RA. The role of the mannose-binding lectin in innate immunity. Clin Infect Dis. 2005;41(Suppl 7):S440–4.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  73. 73.

    Collard CD, Vakeva A, Morrissey MA, Agah A, Rollins SA, Reenstra WR, et al. Complement activation after oxidative stress: role of the lectin complement pathway. Am J Pathol. 2000;156(5):1549–56.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  74. 74.

    Kim YH, Barclay JL, He J, Luo X, O’Neill HM, Keshvari S, et al. Identification of carboxypeptidase X (CPX)-1as a positivie regulator of adipogenesis. FASEB J. 2016;30(7):2528–40. https://doi.org/10.1096/fj.201500107R.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  75. 75.

    Rho SB, Byun HJ, Park SY, Chun T. Calpain 6 supports tumorigenesis by inhibiting apoptosis and facilitating angiogenesis. Cancer Lett. 2008;271(12):306–13 https://doi.org/10.1016/j.canlet.2008.06.020.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  76. 76.

    Zhao SH, Tuggle CK. Linkage mapping and expression analyses of a novel gene, placentally expressed transcript 1 (PLET1) in the pig. Anim Gen. 2003;35(1):72–7 https://doi.org/10.1046/j.0268-9146.2003.01075.x.

    Article  Google Scholar 

  77. 77.

    Everaert C, Luypaert M, Maag JLV, Cheng QX, Dinger ME, Hellemans J, Mestdagh P. Benchmarking of RNA-sequencing analysis workflows using whole-transcriptome RT-qPCR expression data. Sci Rep 7:1559. Doi:https://doi.org/10.1038/s41598-017-01617-3.

  78. 78.

    Vatnick I, Schoknecht PA, Darrigrand R, Bell AW. Growth and metabolism of the placenta after unilateral fetectomy in twin pregnant ewes. J Dev Physiol. 1991;15(6):351–6.

    CAS  PubMed  PubMed Central  Google Scholar 

  79. 79.

    Andrews S. http://www.bioinformatics.babraham.ac.uk/projects/fastqc2010. Available from: http://www.bioinformatics.babraham.ac.uk/projects/fastqc.

  80. 80.

    Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30(15):2114–20.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  81. 81.

    Wu TD, Nacu S. Fast and SNP-tolerant detection of complex variants and splicing in short reads. Bioinformatics. 2010;26(7):873–81.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  82. 82.

    Liao Y, Smyth GK, Shi W. featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics. 2014;30(7):923–30.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  83. 83.

    Robinson MD, McCarthy DJ, Smyth GK. edgeR: a bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26(1):139–40.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  84. 84.

    Duckett SK, Furusho-Garcia I, Rico JE, McFadden JW. Flaxseed oil or n-7 fatty acid-enhanced fish oil supplementation alters fatty acid composition, plasma insulin and serum ceramide concentrations, and gene expression in lambs. Lipids. 2019;54:389–99. https://doi.org/10.1002/lipd.12156.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  85. 85.

    Xie F, Xiao P, Chen D, Xu L, Zhang B. miRDeepFinder: a miRNA analysis tool for deep sequencing of plant small RNAs. Plant Mol Biol. 2012;80:75–84.

Download references

Acknowledgments

Appreciation is expressed to: M. C. Miller for animal handling, R. Smith for laboratory assistance, and the AVS 4220 Fetal Fescue Research Team for assistance.

Consent to participate

Not applicable.

Funding

Technical Contribution No. 6774 of the Clemson University Experiment Station. This research was supported by USDA Agriculture and Food Research Initiative Competitive Grant no. 2015–67015-23218. Clemson University Genomics and Bioinformatics Facility is supported by Grant P20GM109094 an Institutional Development Award (IDeA) from the National Institute of General Medical Sciences of the National Institutes of Health. The funders had no role in the design of the study or collection, analysis, and interpretation of data and in writing the manuscript.

Author information

Affiliations

Authors

Contributions

JLB performed RNA extractions, analyzed the DEG data with guidance from RN, interpreted data, and wrote the first draft of the manuscript. RN offered additional guidance concerning bioinformatics data processing. SKD conceived the original research study and was the principle investigator. All authors have read and approved the manuscript.

Corresponding author

Correspondence to S. K. Duckett.

Ethics declarations

Ethics approval

All animal experimental procedures were reviewed and approved by the Clemson University Institutional Animal Care and Use Committee (AUP 2014–081).

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Britt, J.L., Noorai, R.E. & Duckett, S.K. Differentially expressed genes in cotyledon of ewes fed mycotoxins. BMC Genomics 21, 680 (2020). https://doi.org/10.1186/s12864-020-07074-z

Download citation

Keywords

  • Sheep
  • Ergot alkaloid
  • Mycotoxin
  • Cotyledon
  • Vasoconstriction
  • Nutrient transport