Research article | Open | Published:
Regulation of gene expression in the bovine blastocyst by colony stimulating factor 2
BMC Research Notesvolume 9, Article number: 250 (2016)
Colony stimulating factor 2 can have multiple effects on the function of the preimplantation embryo that include increased potential to develop to the blastocyst stage, reduced apoptosis, and enhanced ability of inner cell mass (ICM) to remain pluripotent after culture. The objective of the current experiment was to identify genes regulated by CSF2 in the ICM and trophectoderm (TE) of the bovine blastocyst with the goal of identifying possible molecular pathways by which CSF2 increases developmental competence for survival. Embryos were produced in vitro and cultured from Day 6 to 8 in serum-free medium containing 10 ng/ml recombinant bovine CSF2 or vehicle. Blastocysts were harvested at Day 8 and ICM separated from TE by magnetic-activated cell sorting. RNA was purified and used to prepare amplified cDNA, which was then subjected to high-throughput sequencing using the SOLiD 4.0 system. Three pools of amplified cDNA were analyzed per treatment.
The number of genes whose expression was regulated by CSF2, using P < 0.05 and >1.5-fold difference as cut-offs, was 945 in the ICM (242 upregulated by CSF2 and 703 downregulated) and 886 in the TE (401 upregulated by CSF2 and 485 downregulated). Only 49 genes were regulated in a similar manner by CSF2 in both cell types. The three significant annotation clusters in which genes regulated by ICM were overrepresented were related to membrane signaling. Genes downregulated by CSF2 in ICM were overrepresented in several pathways including those for ERK and AKT signaling. The only significant annotation cluster containing an overrepresentation of genes regulated by CSF2 in TE was for secreted or extracellular proteins. In addition, genes downregulated in TE were overrepresented in TGFβ and Nanog pathways.
Differentiation of the blastocyst is such that, by Day 8 after fertilization, the ICM and TE respond differently to CSF2. Analysis of the genes regulated by CSF2 in ICM and TE are suggestive that CSF2 reinforces developmental fate and function of both cell lineages.
Maternally-derived signals called embryokines can alter developmental phenotype of the preimplantation embryo immediately or later during pre-natal or post-natal life [1, 2]. CSF2 is one of the best-studied embryokines and has been shown to affect developmental processes in preimplantation embryo of the cow, mouse, human, and pig [2, 3]. In the cow, the species studied here, treatment of embryos with CSF2 from Day 5–7 of development can increase the proportion of embryos becoming blastocysts [4–6], block apoptosis of morula-stage embryos in response to heat shock , improve ability of isolated inner cell mass (ICM) to survive passage while remaining in a pluripotent state , increase the proportion of embryos developing to term after transfer into recipient females [8, 9] and result in calves with increased growth rate after birth .
The objective of the current experiment was to identify genes regulated by CSF2 in the ICM and trophectoderm (TE) of the bovine blastocyst with the goal of identifying possible molecular pathways by which CSF2 increases developmental competence for survival.
Data were obtained as part of an experiment on gene expression in isolated ICM and TE of the bovine blastocyst that was reported earlier . The earlier paper contained data from control embryos only and the current report details additional information about differences in gene expression between control embryos and embryos treated with CSF2 from Day 6 to 8 of development.
Embryo production, isolation of ICM and TE, and RNA-Seq
These procedures were described in detail earlier . Briefly, bovine embryos were produced by fertilizing slaughterhouse-derived oocytes with Percoll-purified spermatozoa from a pool of frozen-thawed semen from three bulls. Fertilization proceeded for 18–20 h and then embryos were cultured in 50 µl SOF-BE1 medium  at 38.5 °C in a humidified atmosphere of 5 % CO2 and 5 % O2 with the balance N2. At Day 6, an additional 5 µl culture medium was added containing either medium alone (control) or, for the CSF2 group, 110 ng/ml recombinant bovine CSF2 (CIBA-GEIGY, Basel, Switzerland) to achieve a final CSF2 concentration of 10 ng/ml CSF2. Blastocysts were harvested at Day 8 after fertilization and used to prepare preparations of ICM and TE using magnetic-activated cell sorting as described . Three separate pools of TE and ICM for each treatment were obtained. Each pool was prepared using 88–102 blastocysts.
Total RNA isolated from each pool of embryonic cells was used to prepare amplified cDNA using the Ovation RNA-Seq kit (NuGen Technology, San Carlos, CA, USA). Barcoded fragment libraries were constructed using the SOLiD™ v4 fragment library kit. Amplified bead-bound libraries were clonally amplified by emulsion PCR according to the Applied Biosystems SOLiD™ 4 Systems Templated Bead Preparation Guide (Applied Biosystems, Foster City, CA, USA). Clonally amplified template DNA were P2-enriched and extended with a bead linker by terminal transferase. Approximately 600–700 M beads were deposited on each slide and sequenced using ‘sequencing by ligation’ chemistry and the 50 × 5 bp protocol on the SOLiD™ v4 sequencer (Applied Biosystems) at the Interdisciplinary Center for Biotechnology Research, University of Florida. Results were obtained as color space fasta files. Processing of reads was described earlier . A total of 64.4 % or 595 million reads were mapped successfully. Data were deposited as Submission DRA000504 in the DDBJ Sequence Read Archive at http://www.ddbj.nig.ac.jp/index-e.html.
Analysis of RNA-Seq data
The number of mapped reads for each individual gene was counted using the HTSeq tool (http://www-huber.embl.de/users/anders/HTSeq/doc/overview.html) with intersection-nonempty mode. The DESeq package  in R was used to determine treatment effects on gene expression. Treatment effects were determined by making the following comparisons: ICM-control vs ICM-CSF2 and TE-control vs TE-CSF2. Two separate sets of differentially expressed genes were identified. The first included those with a P value of <0.05 and a minimum fold-change of >1.5 or <0.66 (CSF2/control) and the second included the same fold-change cut-off and a P < 0.05 after adjustment for a 2 % false discovery rate using the Benjamini-Hochberg procedure .
The differentially expressed genes in TE and ICM were analyzed by the Database for Annotation, Visualization and Integrated Discovery (DAVID; (DAVID Bioinformatics Resources 6.7, http://david.abcc.ncifcrf.gov/) . Ambiguous gene identifiers were submitted to the gene conversion tool to create DAVID gene names. The DAVID database was queried to identify functional annotation terms and clusters in which differentially expressed genes were overrepresented. Clusters of annotation terms with Benjamini-Hochberg corrected P values of <0.05 were identified. Lists of upregulated and downregulated genes in ICM and TE were also separately analyzed for overrepresentation of genes in SuperPaths using GeneAnalytics (https://ga.genecards.org). SuperPaths are a collection of 1073 unified pathways generated from multiple pathway sources using nearest neighbor graphing and hierarchical clustering . All SuperPaths with an adjusted P value of <0.0001 are reported as well as selected paths with an adjusted P value of <0.05. Lists of differentially expressed genes in ICM and TE were separately uploaded to Ingenuity Pathway Analysis (IPA, Qiagen, Redwood City, CA, USA) for 1) clustering of differentially expressed genes into diseases and functions predicted to be changed as a result of differential gene expression and 2) identification of putative transcriptional regulators of differentially expressed genes. A function was predicted to be changed when the activation z-score ≥2.0. For upstream regulators, transcriptional regulators that were predicted to be activated or inhibited by CSF2 (i.e., activation z-score ≥2.0) are reported.
Genes regulated by CSF2
The effects of CSF2 on gene expression are described in Additional file 1. The number of genes whose expression was regulated by CSF2, using P < 0.05 and >1.5-fold or <0.66-fold difference as cut-offs, was 945 in the ICM (242 upregulated by CSF2 and 703 downregulated) and 886 in the TE (401 upregulated by CSF2 and 485 downregulated). A total of 74 genes were regulated by CSF2 in both ICM and TE (Additional file 1). Of these, 9 were upregulated by CSF2 in both cell types, 40 were downregulated by CSF2 in both cell types, 13 were upregulated by CSF2 in ICM and downregulated in TE and 12 were downregulated in ICM and upregulated in TE.
When a more stringent criterion was applied, an adjusted P value of <0.05 after performing a Benjamini-Hochberg correction for a false discovery rate of 2 %, CSF2 changed expression of 25 genes in the ICM (8 upregulated and 17 downregulated by CSF2) and 23 genes in the TE (14 upregulated and 9 downregulated). The list of these genes for ICM is shown in Table 1. The largest increase in expression caused by CSF2 was for FRMPD3 (14.7-fold), NDNF (13.3-fold), LAPTM5 (12.3-fold) and an ortholog of OLFR1239 (8.5-fold). The largest decrease in expression was for a bovine ortholog of MYADM (undetectable in cells from CSF2-treated embryos), olfactory receptor, family 5, subfamily T like (0.11control/CSF2 or a 19.9-fold decrease), CMYA5 (19.7-fold decrease) and another family 5 olfactory receptor (OR5T2; 8.8-fold decrease).
Results for genes regulated by CSF2 in TE after correcting for multiple testing are summarized in Table 2. The largest increase in expression caused by CSF2 was for CLIC6 (undetectable in the absence of CSF2 and present in CSF2-treated TE), the olfactory receptor OR5AN1 (16.2-fold), CSN3 (13.1-fold; also upregulated 4.6-fold in ICM) and a uncharacterized gene (11.5-fold). The largest decrease was for a gene orthologous to interferon-regulatory factor (19.2-fold decrease), NPY1R (17.8-fold decrease), PAQR7 (17.0-fold decrease) and VAV3 (4.8-fold decrease).
Functional annotation of differentially expressed genes
A total of 671 differentially expressed genes in the ICM was annotated by DAVID. There were three annotation clusters having Benjamini-Hochberg adjusted P values <0.05 (Additional file 2). These included a cluster of 194 genes related to membrane signaling (enrichment score = 4.26), another of 254 genes related to membrane and extracellular space (enrichment score = 2.94) and a cluster of 30 genes related to sensory perception (enrichment score = 2.36). Of note was the large number of differentially expressed genes associated with G-protein coupled receptor protein signaling pathway (67 genes) and olfactory receptor activity (48 genes, adjusted P value = 0.013). CSF2 increased expression of some olfactory genes and decreased expression of others. For example, LOC100294956 (olfactory receptor 4A47) was increased 8.48-fold and OR5T2 was decreased 9.09-fold (Table 1). The absolute number of reads for most putative olfactory receptor genes was low (<5).
For differentially expressed genes in the TE, 648 differentially expressed genes were annotated by DAVID. The only cluster of pathways with a Benjamini-Hochberg adjusted P value <0.05 was for 108 genes related to secreted or extracellular proteins (enrichment score = 4.22) (Additional file 2).
GeneAnalytics was used to identify SuperPaths in which genes regulated by CSF2 were overrepresented. Using a P value of 0.001 as a cutoff, there were no SuperPaths identified for genes upregulated by CSF2 in ICM. Among genes downregulated by CSF2 were those in the ERK signaling, integrin pathway, visual cycle and AKT signaling (Table 3). There were 11 Superpaths in which genes upregulated by CSF2 in TE were overrepresented (Table 4). Among these were those for apoptotic pathways in synovial fibroblasts, G protein signaling Ras family GTPases in kinase cascades, and ERK signaling. There were four SuperPaths containing genes downregulated by CSF2 in TE with P < 0.001 including for CREB pathway, actin nucleation by ARP-WASP complex, Akt signaling, and ERK signaling (Table 5). In addition, there were two SuperPaths related to pluripotency and differentiation of the blastocyst (TGF-β, Nanog in mammalian ESC pluripotency) that contained an overrepresentation of downregulated genes in the TE at P < 0.05.
Sets of genes that were differentially expressed in ICM and TE were separately subjected to analysis by IPA to determine 1) diseases and functions predicted to be activated or inhibited as a result of changes in gene expression caused by CSF2 and 2) transcriptional regulators that were postulated to be activated or inhibited in response to CSF2. For ICM, a total of 22 functions (many overlapping) were predicted to change in response to alterations in gene expression caused by CSF2 (Additional file 3). In all cases, the function was predicted to decrease. Among these functions were 9 related to phagocytosis or cell movement, 2 related to changes in cell morphology, and 3 related to immune responses. For TE, there were 6 functions predicted to change in response to CSF2 and all were predicted to decrease. These functions were related to antibody response, quantity of cells and cell migration, and dermatitis. Three transcription factors were predicted to be activated by CSF2 in ICM (HLX, IKZF1, and ZNF217) while 13 were predicted to be inhibited (HSF1, STAT3, SRF, STAT4, CREM, SP1, CREB1, FOXM1, ERG, ETS1, SMARCA4, CREBBP and CTNNB1) (Additional file 3). Only two transcription factors were predicted to be activated by CSF2 in TE (SATB1 and NFATC2) and none were predicted to be inhibited.
Comparison to genes regulated by CSF2 in bovine morulae and mouse blastocysts
The set of genes regulated by CSF2 in the ICM and TE were compared to sets of genes previously reported to be regulated by CSF2 in embryos. In an earlier experiment with bovine embryos, Loureiro el al.  identified 67 genes in Day 6 morula that were upregulated by CSF2 treatment at Day 5 and 92 genes that were downregulated. Of these 159 genes, only 7 were regulated by CSF2 in a similar manner for Day 8 blastocysts. Of genes upregulated by CSF2 in morula, SLC26A5 was upregulated in ICM, ANKRD37 was upregulated in TE, and PYGO1 was upregulated in ICM and TE. Of genes downregulated in morula, DTX3 and VGLL1 were downregulated in ICM and NEDD4 and PHGDH were downregulated in TE.
In another experiment, Chin et al.  identified 16 genes related to stress responses that were downregulated in mouse blastocysts treated with CSF2. Two of these genes, DFF5 and HSPH1, were also downregulated in ICM in the current experiment.
CSF2 can have multiple effects on the function of the developing bovine embryo. These include increased potential to develop to the blastocyst stage [4–6], reduced apoptosis , and enhanced ability of ICM to remain pluripotent after culture . As a result of these or other actions, embryos treated with CSF2 have increased ability to develop to term after transfer to recipients [8, 9]. Accordingly, CSF2 could be a useful molecule for commercial production of bovine embryos based on in vitro fertilization and culture. In the current study, it was shown that effects of CSF2 on embryonic function are accompanied by changes in gene expression in both the ICM and TE of the blastocyst. Moreover, the actions of CSF2 on gene expression were largely distinct for the two cell types. Of 1782 genes regulated by CSF2 > 1.5-fold, only 49 genes were regulated by CSF2 in the same direction in both cell types. Thus, differentiation of the blastocyst is such that, by Day 8 after fertilization, the ICM and TE respond differently to maternal signals. Analysis of the genes regulated by CSF2 in ICM and TE are suggestive that CSF2 reinforces developmental fate and function of both cell lineages. Actions of CSF2 on the ICM reduced expression of genes in signaling pathways involved in differentiation while actions of CSF2 on the TE reduced expression of genes in pluripotency pathways and may have increased capacity of the TE to interact with the maternal system through upregulation of genes involved in extracellular processes.
Many of the genes regulated by CSF2 in the ICM are involved in signal transduction. The three significant annotation clusters in which genes regulated by CSF2 were overrepresented are related to membrane signaling. The overall picture is for downregulation of genes involved in signal transduction, particularly for the ERK/MAPK and AKT signaling pathways. Repression of MAPK signaling by CSF2 may facilitate maintenance of the ICM in a pluripotent condition. In cattle, inhibition of MAPK signaling increases numbers of cells positive for NANOG while reducing expression of the TE gene IFNT and cells positive for the hypoblast marker GATA6 [18, 19]. Inhibition of MAPK also promotes formation of naïve-type induced pluripotent stem cells in the cow . Downregulation of genes involved in AKT signaling would counteract the downregulation of MAPK signaling because this pathway can itself inhibit MAPK signaling . However, inhibition of PI3K (upstream regulator of AKT) did not affect expression of pluripotency-related genes in the bovine blastocyst although ICM numbers were reduced . CSF2 can signal through the phosphatidylinositol 3-kinase/AKT pathway  and downregulation of expression of genes in the AKT pathway by CSF2 could reflect feedback inhibition. Three transcription factors that are involved in actions of CSF2 in other cells, CTTNB1, SP1, and STAT3 [24, 25], were also predicted to be inhibited in the ICM by CSF2.
Another indication that CSF2 promotes maintenance of the ICM in a pluripotent state was the observation that the transcription factor ETS1 was predicted to be inhibited in the ICM. ETS1 promotes differentiation of human trophoblast  and a related molecule, ETS2, is involved in transcriptional regulation of trophoblast specific genes like IFNT and TKDP1 . Another transcriptional regulator predicted to be inhibited by CSF2 in the ICM, STAT3, is important for formation of the ICM and expression of pluripotency markers in cattle  while another, CREBBP participates in downregulation of the trophoblast gene IFNT .
The TE is the cell type in the blastocyst that interacts most closely with the maternal reproductive tract, both because of its location at the outer boundaries of the embryo and because tight junctions between TE cells  limit movement of molecules from the ICM to the extra-embryonic space. CSF2 may facilitate interactions between the TE and maternal cells because genes regulated by CSF2 in the TE included a large cluster of genes related to secreted or extracellular proteins. Among the genes whose expression in TE was downregulated by CSF2 was MUC1, which encodes for a mucin that inhibits adhesion to cell surfaces. MUC1 is present on human trophoblast and overexpression in the Jar trophoblast cell line reduced invasiveness of the cells . Perhaps inhibition of MUC1 expression by CSF2 increases the competence of the TE to interact with the endometrium. In addition, MUC1 can participate in cell signaling because its cytoplasmic tail can regulate β-catenin stability .
One gene whose expression in TE was increased by CSF2 was FOS. The promoter for the maternal recognition of pregnancy protein, IFNT, contains an AP1 site  and deletions in that site can decrease IFNT1 gene expression . Treatment of embryos with CSF2 from Day 5–7 caused an increase in IFNT accumulation in the uterus at Day 15 of pregnancy for male embryos  and it is possible that regulation of FOS is involved in this action of CSF2. If so, effects of CSF2 could be sex-specific because CSF2 decreased IFNT accumulation in the uterus of cows that received female embryos .
CSF2 may also promote differentiation of the TE. Unlike for the mouse, the TE of the bovine embryo only gradually becomes committed to the TE fate and, when aggregated with 8-cell embryos, some TE cells from the day 7 blastocyst can revert to ICM . Several genes downregulated by CSF2 in TE were identified as members of the Nanog and TGF-β signaling pathways. NANOG is a well-described promoter of pluripotency in the ICM  and TGF-β signaling has been reported to promote pluripotency in human embryonic stem cells  but to reduce epiblast formation in the mouse embryo .
One of the surprising set of signaling genes regulated by CSF2 in the ICM was that for olfactory receptors. Some of these genes were increased in expression by CSF2 and others decreased. Several olfactory receptors were also regulated by CSF2 in the TE including OR5AN1, which was increased 16-fold by CSF2 (Table 2). Olfactory receptors are a large gene family. A total of 1071 olfactory receptor related sequences were identified in the bovine genome including 881 functional genes with the remainder being pseudogenes or partial sequences . The fact that specific olfactory receptors are regulated by CSF2 implies that olfactory receptor signaling may be important for embryonic development and that actions of specific odorants on the embryo could be modified by CSF2.
CSF2 is best known for its role in the immune system, where in addition to promoting formation of granulocytes and macrophages from hematopoietic stem cells, it also controls neutrophil function and macrophage differentiation [39, 40]. Regulation of genes involved in immune function by CSF2 in the embryo accounts for the observation that many of the functions predicted to change in response to CSF2 were related to immune function.
The finding that the genes regulated by CSF2 in the ICM were distinct from those regulated by CSF2 in the TE strengthens the idea that transcriptional regulation by CSF2 depends on cell type. In addition, the genes regulated by CSF2 in the ICM and TE were largely different than those found to be regulated by CSF2 in the bovine morula . The 7 genes that were regulated similarly by CSF2 in morula and blastocysts include 3 upregulated genes (ANKRD37, a putative nuclear localization signal, PYGO1, a coactivator of transcription of WNT-regulated genes, and SLC26A5, an incomplete anion transporter involved in changing cell shape) and 4 downregulated genes (DTX3 and NEDD4, both involved in ubiquitination, PHGDH, an enzyme involved in serine synthesis, and VGLL1, a coactivator of transcription factors). In the bovine, activation of WNT signaling either reduced [41, 42] or increased  competence of embryos to become blastocysts and preferentially reduced number of TE cells as compared to ICM cells . Deubiquitinating enzymes are required for blastocyst formation in mice and pigs . Serine is probably not a major energy source for the bovine embryo because it is not taken up by blastocysts to a large degree  and addition of the amino acid to culture medium of bovine embryos had no effect on development to the blastocyst stage .
CSF2 can cause differential actions on gene expression in the ICM and TE. One interpretation of changes in gene expression is that CSF2 reinforces the lineage commitment decisions in ICM and TE and enhances ability of the TE to communicate with the maternal reproductive tract. Additional functional studies will be required to test these hypotheses.
inner cell mass
Hansen PJ, Denicol AC, Dobbs KB. Maternal embryokines that regulate development of the bovine preimplantation embryo. Turkish J Vet Anim Sci. 2014;38:589–98.
Hansen PJ, Dobbs KB, Denicol AC, Siqueira LG. Sex and the preimplantation embryo: implications of sexual dimorphism in the preimplantation period for maternal programming of embryonic development. Cell Tiss Res. 2016;363:237–47.
Hansen PJ, Dobbs KB, Denicol AC. Programming of the preimplantation embryo by the embryokine colony stimulating factor 2. Anim Reprod Sci. 2014;149:59–66.
de Moraes AA, Hansen PJ. Granulocyte-macrophage colony-stimulating factor promotes development of in vitro produced bovine embryos. Biol Reprod. 1997;57:1060–5.
Loureiro B, Bonilla L, Block J, Fear JM, Bonilla AQ, Hansen PJ. Colony-stimulating factor 2 (CSF-2) improves development and post transfer survival of bovine embryos produced in vitro. Endocrinology. 2009;150:5046–54.
Dobbs KB, Khan FA, Sakatani M, Moss JI, Ozawa M, Ealy AD, Hansen PJ. Regulation of pluripotency of inner cell mass and growth and differentiation of trophectoderm of the bovine embryo by colony stimulating factor 2. Biol Reprod. 2013;89:141.
Loureiro B, Oliveira LJ, Favoreto MG, Hansen PJ. Colony-stimulating factor 2 inhibits induction of apoptosis in the bovine preimplantation embryo. Am J Reprod Immunol. 2011;65:578–88.
Loureiro B, Bonilla L, Block J, Fear JM, Bonilla AQ, Hansen PJ. Colony-stimulating factor 2 (CSF-2) improves development and post transfer survival of bovine embryos produced in vitro. Endocrinology. 2009;150:5046–54.
Denicol AC, Block J, Kelley DE, Pohler KG, Dobbs KB, Mortensen CJ, Ortega MS, Hansen PJ. The WNT signaling antagonist Dickkopf-1 directs lineage commitment and promotes survival of the preimplantation embryo. FASEB J. 2014;28:3975–86.
Kannampuzha-Francis J, Denicol AC, Loureiro B, Kaniyamattam K, Ortega MS, Hansen PJ. Exposure to colony stimulating factor 2 during preimplantation development increases postnatal growth in cattle. Mol Reprod Dev. 2015;82:892–7.
Ozawa M, Sakatani M, Yao J, Shanker S, Yu F, Yamashita R, Wakabayashi S, Nakai K, Dobbs KB, Sudano MJ, Farmerie WG, Hansen PJ. Global gene expression of the inner cell mass and trophectoderm of the bovine blastocyst. BMC Dev Biol. 2012;12:33.
Fields SD, Hansen PJ, Ealy AD. Fibroblast growth factor requirement for in vitro development of bovine embryos. Theriogenology. 2011;75:1466–75.
Anders S, Huber W. Differential expression analysis for sequence count data. Genome Biol. 2010;11:R106.
Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J Royal Stat Soc. 1995;57:289–300.
Huang DW, Sherman BT, Lempicki RA. Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 2009;37:1–13.
Belinky F, Nativ N, Stelzer G, Zimmerman S, Iny Stein T, Safran M, Lancet D. PathCards: multi-source consolidation of human biological pathways. Database 2015; doi:10.1093/database/bav006
Chin PY, Macpherson AM, Thompson JG, Lane M, Robertson SA. Stress response genes are suppressed in mouse preimplantation embryos by granulocyte-macrophage colony-stimulating factor (GM-CSF). Hum Reprod. 2009;24:2997–3009.
Kuijk EW, van Tol LT, Van de Velde H, Wubbolts R, Welling M, Geijsen N, Roelen BA. The roles of FGF and MAP kinase signaling in the segregation of the epiblast and hypoblast cell lineages in bovine and human embryos. Development. 2012;139:871–82.
Brinkhof B, van Tol HT, Groot Koerkamp MJ, Riemers FM, IJzer SG, Mashayekhi K, Haagsman HP, Roelen BA. A mRNA landscape of bovine embryos after standard and MAPK-inhibited culture conditions: a comparative analysis. BMC Genom. 2015;16:277.
Kawaguchi T, Tsukiyama T, Kimura K, Matsuyama S, Minami N, Yamada M, Imai H. Generation of naïve bovine induced pluripotent stem cells using piggybac transposition of doxycycline-inducible transcription factors. PLoS One. 2015;10:e0135403.
Dalton S. Signaling networks in human pluripotent stem cells. Curr Opin Cell Biol. 2013;25:241–6.
Meng F, Forrester-Gauntlett B, Turner P, Henderson H, Oback B. Signal inhibition reveals JAK/STAT3 pathway as critical for bovine inner cell mass development. Biol Reprod. 2015;93:132.
Jeong W, Kim J, Bazer FW, Song G. Proliferation-stimulating effect of colony stimulating factor 2 on porcine trophectoderm cells is mediated by activation of phosphatidylinositol 3-kinase and extracellular signal-regulated kinase 1/2 mitogen-activated protein kinase. PLoS One. 2014;9:e88731.
Mangan JK, Tantravahi RV, Rane SG, Reddy EP. Granulocyte colony-stimulating factor-induced upregulation of Jak3 transcription during granulocytic differentiation is mediated by the cooperative action of Sp1 and Stat3. Oncogene. 2006;25:2489–99.
Brown AL, Salerno DG, Sadras T, Engler GA, Kok CH, Wilkinson CR, Samaraweera SE, Sadlon TJ, Perugini M, Lewis ID, Gonda TJ, D’Andrea RJ. The GM-CSF receptor utilizes β-catenin and Tcf4 to specify macrophage lineage differentiation. Differentiation. 2012;83:47–59.
Kessler CA, Stanek JW, Stringer KF, Handwerger S. ETS1 induces human trophoblast differentiation. Endocrinology. 2015;156:1851–9.
Chakrabarty A, Roberts MR. Ets-2 and C/EBP-beta are important mediators of ovine trophoblast Kunitz domain protein-1 gene expression in trophoblast. BMC Mol Biol. 2007;8:14.
Sakurai T, Bai H, Bai R, Sato D, Arai M, Okuda K, Ideta A, Aoyagi Y, Godkin JD, Imakawa K. Down-regulation of interferon tau gene transcription with a transcription factor EOMES. Mol Reprod Dev. 2013;80:371–83.
Eckert JJ, Fleming TP. Tight junction biogenesis during early development. Biochim Biophys Acta. 2008;1778:717–28.
Shyu MK, Lin MC, Liu CH, Fu YR, Shih JC, Lee CN, Chen HY, Huang J, Huang MC, Hsieh FJ. MUC1 expression is increased during human placental development and suppresses trophoblast-like cell invasion in vitro. Biol Reprod. 2008;79:233–9.
Poh TW, Bradley JM, Mukherjee P, Gendler SJ. Lack of Muc1-regulated & β-catenin stability results in aberrant expansion of CD11b + Gr1 + myeloid-derived suppressor cells from the bone marrow. Cancer Res. 2009;69:3554–62.
Ezashi T, Roberts RM. Regulation of interferon-τ (IFN-τ) gene promoters by growth factors that target the Ets-2 composite enhancer: a possible model for maternal control of IFN-τ production by the conceptus during early pregnancy. Endocrinology. 2004;145:4452–60.
Dobbs KB, Gagné D, Fournier E, Dufort I, Robert C, Block J, Sirard MA, Bonilla L, Ealy AD, Loureiro B, Hansen PJ. Sexual dimorphism in developmental programming of the bovine preimplantation embryo caused by colony-stimulating factor 2. Biol Reprod. 2014;91:80.
Berg DK, Smith CS, Pearton DJ, Wells DN, Broadhurst R, Donnison M, Pfeffer PL. Trophectoderm lineage determination in cattle. Dev Cell. 2011;20:244–55.
Artus J, Chazaud C. A close look at the mammalian blastocyst: epiblast and primitive endoderm formation. Cell Mol Life Sci. 2014;71:3327–38.
James D, Levine AJ, Besser D, Hemmati-Brivanlou A. TGFbeta/activin/nodal signaling is necessary for the maintenance of pluripotency in human embryonic stem cells. Development. 2005;132:1273–82.
Ghimire S, Heindryckx B, Van der Jeught M, Neupane J, O’Leary T, Lierman S, De Vos WH, Chuva de Sousa Lopes S, Deroo T, De Sutter P. Inhibition of transforming growth factor β signaling promotes epiblast formation in mouse embryos. Stem Cells Dev. 2015;24:497–506.
Lee K, Nguyen DT, Choi M, Cha SY, Kim JH, Dadi H, Seo HG, Seo K, Chun T, Park C. Analysis of cattle olfactory subgenome: the first detail study on the characteristics of the complete olfactory receptor repertoire of a ruminant. BMC Genom. 2013;14:596.
Pitrak DL. Effects of granulocyte colony-stimulating factor and granulocyte-macrophage colony-stimulating factor on the bactericidal functions of neutrophils. Curr Opin Hematol. 1997;4:183–90.
Page AV, Liles WC. Colony-stimulating factors in the prevention and management of infectious diseases. Infect Dis Clin North Am. 2011;25:803–17.
Aparicio IM, Garcia-Herreros M, Fair T, Lonergan P. Identification and regulation of glycogen synthase kinase-3 during bovine embryo development. Reproduction. 2010;140:83–92.
Denicol AC, Dobbs KB, McLean KM, Carambula SF, Loureiro B, Hansen PJ. Canonical WNT signaling regulates development of bovine embryos to the blastocyst stage. Sci Rep. 2013;3:1266.
Mtango NR, Latham KE, Sutovsky P. Deubiquitinating enzymes in oocyte maturation, fertilization and preimplantation embryo development. Adv Exp Med Biol. 2014;759:89–110.
Jung YG, Sakata T, Lee ES, Fukui Y. Amino acid metabolism of bovine blastocysts derived from parthenogenetically activated or in vitro fertilized oocytes. Reprod Fertil Dev. 1998;10:279–87.
Lee ES, Fukui Y, Lee BC, Lim JM, Hwang WS. Promoting effect of amino acids added to a chemically defined medium on blastocyst formation and blastomeres proliferation of bovine embryos cultured in vitro. Anim Reprod Sci. 2004;84:257–67.
All authors except PJH performed experiments, data were analyzed by MO and PJH. All authors read and approved the final manuscript.
Authors thank Savita Shanker, JiQiang Yao, and Fahong Yu and the University of Florida Interdisciplinary Center for Biotechnology Research, for performing RNA Seq and bioinformatics, William Rembert, for collecting oocytes, the Chernin family and Central Packing (Center Hill, Florida), for donating ovarian tissue, Scott Randel of Southeastern Semen (Wellborn Florida), for donating semen, and Novartis (Basle, Switzerland), for donation of CSF2.
The authors declare that they have no competing interests.
This work was supported by Agriculture and Food Research Initiative Grant Nos. 2009-65203-05732 and 2011-67015-30688 from United States Dept. of Agriculture National Institute of Food and Agriculture, the Research Opportunity Fund, University of Florida Research Foundation, and NIH Grant R03 HD080855.