Gene expression of bovine embryos developing at the air-liquid interface on oviductal epithelial cells (ALI-BOEC)
- Vera A. van der Weijden†1,
- Shuai Chen†2,
- Stefan Bauersachs1,
- Susanne E. Ulbrich1 and
- Jennifer Schoen2Email authorView ORCID ID profile
© The Author(s). 2017
Received: 22 June 2017
Accepted: 9 November 2017
Published: 25 November 2017
We recently developed an air-liquid interface long-term culture of differentiated bovine oviductal epithelial cells (ALI-BOEC). This ex vivo oviduct epithelium is capable of supporting embryo development in co-culture up to the blastocyst stage without addition of embryo culture medium. However, blastocyst rates in co-culture were markedly lower than in conventional in vitro embryo production procedures. In the present study, we assessed target gene expression of ALI-BOEC derived embryos to test their similarity to embryos from conventional in vitro embryo culture. We screened previously published data from developing bovine embryos and selected 41 genes which are either differentially expressed during embryo development, or reflect differences between various in vitro culture conditions or in vitro and in vivo embryos. Target gene expression was measured in 8-cell embryos and blastocysts using a 48.48 Dynamic Array™ on a Biomark HD instrument. For comparison with the ALI-BOEC system, we generated embryos by two different standard IVP protocols. The culture conditions lead to differential gene expression in both 8-cell embryos and blastocysts. Across the expression of all target genes the embryos developing on ALI-BOEC did not depart from conventional IVP embryos. These first results prove that gene expression in ALI-BOEC embryos is not largely aberrant. However, there was no clear indication for a more in vivo-like target gene expression of these embryos. This calls for further optimization of the ALI-BOEC system to increase its efficiency both quantitatively and qualitatively.
Establishment and maintenance of a viable embryo and subsequent implantation depends on a complex, dynamic network of endocrine signalling pathways and local factors produced by the embryo and the female reproductive tract . These early embryo-maternal interactions are subtle and difficult to investigate in vivo, especially in monotocous species like cattle. Therefore, reliable in vitro tools faithfully mimicking the early embryonic environment are needed. The simple, columnar-shaped oviduct epithelium builds the first maternal “contact zone” for the early embryo. However, when grown under standard culture conditions, oviduct epithelial cells (OEC) lose their polarized phenotype and dedifferentiate both morphologically and functionally [2, 3]. We recently developed an air-liquid interface (ALI) culture system for the long-term culture of differentiated bovine oviductal epithelial cells (ALI-BOEC) . The ALI-BOEC model forms tissue-like epithelial layers, which actively generate an oviductal fluid surrogate on its apical cell surface. Proteomic analyses revealed the abundance of >3000 proteins with distinct similarities to in vivo oviductal fluid. Initial co-culture experiments with living zygotes have proven that co-cultured embryos develop up to the blastocyst stage. The embryonic development is solely supported by the ALI-BOEC milieu, as no embryo culture medium was supplemented. In the present study, we assessed the stage-dependent gene expression of 41 target genes in embryos produced in this ALI-BOEC system. Selected target genes represented different functional categories (imprinted genes and DNA methylation; stem cell markers and differentiation; apoptosis; embryo development, quality and competence; metabolism). For comparison, we used embryos from two different conventional in vitro embryo production (IVP) protocols. We hypothesized that embryos successfully developing in the environment of the ALI-BOEC system exhibit similar gene expression compared to embryos produced in conventional bovine IVP procedures.
Foetal bovine serum, DMEM/Ham’s F12, HEPES, L-glutamine, sodium pyruvate, amphotericin B and penicillin/streptomycin was purchased from Merck Millipore (Billerica, MA, USA). If not otherwise indicated all other reagents were purchased from Sigma Aldrich (St. Louis, MO, USA).
Bovine ovaries and oviducts were by-products from local slaughterhouses (Teterower Fleisch GmbH, Teterow, Germany; Biopark, Güstrow, Germany). All samples were collected within 15 min after slaughter. Ovaries were maintained at 38.5 °C in pre-warmed DPBS with 1% penicillin, while oviducts were transported to the laboratory on ice.
In vitro production of embryos
The IVP procedure was performed as described previously (; experiment 1, protocol 2). In trial 1 (collection of 8-cell embryos; n = 115 oocytes) cleavage rate was 77.39% (n = 89). In trial 2 (collection of blastocysts; n = 95 oocytes) cleavage rate was 72.63% (n = 69) and blastocyst rate 47.83% (n = 33).
Commercial serum-free IVP suite
The IVP procedure was performed using the IVF Bioscience media suite for bovine embryo production according to the manufacturer’s instructions (IVF Bioscience, Falmouth, UK, cat. # 61008, 61002, 61004, 61003, and 61001).
In trial 1 (n = 34 oocytes), the cleavage rate was 97.06% (n = 33). In trial 2 (n = 73 oocytes), the cleavage rate was 91.78% (n = 67) and the blastocyst rate 49.25% (n = 33).
The protocol for culturing ALI-BOEC on hanging inserts has recently been reported  and histological samples to verify differentiation of the cell culture were prepared following published protocols . A scheme of the culture system and a representative histological picture of ALI-BOEC are shown in Fig. 1a and b. The maturation and fertilization of oocytes followed the procedure described in paragraph “1) Standard IVP”. 100 oocytes were applied for trial 1 [cleavage rate 60% (n = 60)]. Another 100 oocytes was applied in trial 2 [cleavage rate 72% (n = 72) and blastocyst rate 9.72% (n = 7)].
Selection of target genes
Forty-one target genes were selected based on previously published transcriptomics and qRT-PCR data of developing bovine embryos. The PubMed database was used to search papers with the following keywords: ‘bovine embryo development’, ‘transcriptomics’, ‘gene expression’, ‘in vivo’ and ‘in vitro fertilization’. The selection criteria included the availability of gene expression data of different days of embryo development, or differences between various in vitro culture conditions or in vivo. Twenty-six publications (Additional file 1: Table S1) were used for a final selection of 41 target genes. The selected genes were subdivided into five functional categories, i.e. ‘imprinted genes and DNA methylation’, ‘stem cell markers and differentiation’, ‘apoptosis’, ‘embryo development, quality and competence’, and ‘metabolism’.
cDNA synthesis and specific target amplification
The CellsDirect™ One-Step qRT-PCR Kit (ThermoFisher Scientific, Waltham, MA, USA, cat. # 11753100) was used for cDNA synthesis and specific-target amplification (STA), as described previously  with minor modifications. In brief, an STA master mix was prepared, consisting of 5 μl 2× Reaction Mix, 0.2 μl CellsDirect Enzyme Mix, 2.5 μl primer mix, 0.2 μl SUPERase• In™ RNase Inhibitor (20 U/μl) (ThermoFisher Scientific, Waltham, MA, USA, cat. # AM2694), and 1 × TE buffer (ThermoFisher Scientific, Waltham, MA, USA, cat. # 12090015). All primers were ordered from Microsynth (Balgach, Switzerland), and product specificity was assessed based on product size by gel electrophoresis. Gene descriptions and primer sequences are available in the Additional file 2: Table S2. Nine microliter of the STA master mix was added to the frozen embryos, which were then homogenized by means of pipetting. Both reverse transcription and STA were performed in a thermal cycler by incubation for 15 min at 50 °C, followed by 2 min incubation at 95 °C. The STA was performed by 18 cycles of 15 s at 95 °C and 4 min at 60 °C. Reactions were cleaned from residual primers by a treatment with an Exonuclease I master mix, consisting of 0.8 μl Exonuclease I (20 U/μl) (ThermoFisher Scientific, Waltham, MA, USA, cat. # EN0581), 0.4 μl 10× Exonuclease I Reaction Buffer, and 2.8 μl nuclease free water. Four microliter of Exonuclease master mix was added to each STA product and incubated at 37 °C for 15 min followed by heat-inactivation of the enzyme at 80 °C for 15 min.
Biomark delta gene assay
The STA samples were used for gene expression analysis using a Biomark HD instrument. The Sample Pre-Mix consisted of 3 μl 2× TaqMan Gene Expression Master Mix (Applied Biosystems, Foster City, CA, USA, cat. # PN 4369016), 0.3 μl 20× DNA Binding Dye Sample Loading Reagent (Fluidigm, San Francisco, CA, USA, cat. # PN 100–0388), 0.3 μl 20 × EvaGreen DNA binding dye (Biotium, Fremont, CA, USA, cat. # PN 31000), and 0.9 μl TE buffer. The Sample Pre-Mix was combined with 1.5 μl of 10× diluted cleaned-up STA cDNA. The Assay Mix was prepared for 48 primer pairs, consisting of 3 μl 2× Assay Loading Reagent (Fluidigm, San Francisco, CA, USA, cat. # PN 85000736), 0.3 μl TE buffer, and 2.7 μl of 20 μM of Forward and Reverse Primer Mix. One 48.48 Dynamic Array™ (Fluidigm, San Francisco, CA, USA, cat. # BMK-M-48.48) chip was loaded and run according to manufacturer instructions as described in the Fluidigm Advanced Development Protocol 14. After a hot-start of 120 s at 50 °C and 600 s at 95 °C, the amplification was performed by 40 cycles of 15 s at 95 °C and 60 s at 60 °C. A melting curve was generated by a temperature increase from 60 °C to 95 °C with increments of 1 °C/s.
The Fluidigm Real-Time PCR Analysis Software was used for quality control of the experiment and for validation of product specificity. A cut-off of Ct = 25 was set as limit of detection. To allow analysis of expression patterns of very low abundant genes a Ct value of 25 was assigned if gene expression levels were lower than the limit of detection. Genes with low expression levels in some samples were PTGS2, FADS1, and HSPA1A (Additional file 3: Table S4; specific samples with expression levels lower than the limit of detection highlighted in blue). Using the geNorm algorithm within the GenEx6 software, the geometric mean of all reference genes, i.e. UBB, H3F3A, YWHAZ, GAPDH, and SDHA, was found to be the most stable reference . The normalized expression (∆Ct) values were log2 transformed and centred for both PCA plot and heatmap, which were processed in R Studio. In addition, the differential gene expression from 8-cell embryos to blastocysts was displayed in a heatmap using the log2 transformed fold changes (∆∆Ct).The statistical analysis of differentially expressed genes (DEGs) was performed on ∆Ct values in IBM SPSS Statistics 23 using one-way ANOVA with Tukey post-hoc test. Graphs of log2 fold-changes were made in GraphPad Prism 7.02.
Results and discussion
The differentiation status of ALI-BOEC was verified by their polarized morphology and presence of ciliated and secretory sub-populations (Fig. 1b). The overall cleavage rate in ALI-BOEC co-culture (ALI, 66%) was comparable to the standard IVP procedure (S, 75.01%), and lower than the cleavage rate reached with the commercial serum-free media suit (SF, 94.42%). However, the blastocyst rate in ALI-BOEC co-culture (9.72%) was considerably lower than in either IVP-S (47.83%) or the commercial IVP-SF system (49.25%). This may be caused by the fact that in vivo bovine zygotes stay in the oviduct for 3–5 days and only develop up to the 8–16 cell stage before entering the uterus . As the oviductal secretions are greatly influenced by hormones during the oestrous cycle , the lack of hormonal stimulation during co-culture may also affect the competence of ALI-BOEC to adequately support embryo development.
Surprisingly, the mRNA expression of OVGP1, a marker gene for OEC , was detected in all 8-cell embryo groups. Initially chosen as a target gene to notice any contamination of the ALI-BOEC embryos with epithelial cells, the expression of OVGP1 mRNA was significantly higher in embryos cultured in commercial serum-free media compared to ALI and was also detectable in standard IVP embryos. The results indicate that a) OVGP1 mRNA is actively expressed in IVP embryos and b) OEC do not necessarily contribute to this expression. The presence and significance of OVGP1 expression in IVP embryos requires further analysis.
The differentially expressed genes in blastocysts derived from the different IVP systems either have a function in embryo development, quality, and competence (CCL26, CDH1, NID2, and IFNAR1) [22–26] or in embryonic metabolism (SLC2A5, SREBP1, SERPINE1) [27–29]. Eighteen of our target genes were previously reported to be differentially expressed between in vivo and in vitro blastocysts (Additional file 5: Table S3) and could therefore be indicative for more in vivo-like embryonic development. However, only three of these genes are among the differentially expressed genes in our analysis (CDH1, NID2 and SLC2A5). ALI-OEC blastocysts show significant lower NID2 expression compared to both conventional IVP groups. In vivo embryos are reported to have lower expression of NID2 compared to embryos produced in conventional IVP systems as well (Additional file 5: Table S3). Likewise, IVP-SF embryos, displayed differential expression of SLC2A5, a fructose transporter, compared to both other groups. SLC2A5 has previously been reported to be higher expressed in vivo than in in vitro embryos. IVP-SF embryos, however, showed a lower expression of this gene. CDH1 expression was significantly different only between ALI and IVC-S embryos.
The ALI-BOEC co-culture system was much less efficient in supporting blastocyst formation than optimized conventional IVP procedures. The different culture conditions lead to differential gene expression in both 8-cell embryos and blastocysts. However, across the expression of all target genes, the embryos developing on ALI-BOEC did not clearly depart from conventional IVP embryos. Our results neither hint for largely aberrant, nor for more in vivo-like gene expression of embryos produced in co-culture with ALI-BOEC. To further optimize the ALI-BOEC system we propose to develop a dynamic hormonal (progesterone) stimulation protocol mimicking the hormonal environment in vivo. Furthermore, the establishment of a sequential culture system of oviductal and uterine epithelial cells might increase the efficiency of the production system not only quantitatively, but also qualitatively in view of its competence to support more in vivo-like embryonic development.
We are grateful to Bianka Drawert and Caterina Poeppel, FBN Dummerstorf Germany, for their technical assistance. We thank Peter Nestorov for sharing his useful insights using the Biomark Delta Gene Assay for low-input purposes. Data produced and analysed in this paper were generated in collaboration with the Genetic Diversity Centre (GDC), ETH Zurich. The authors are active participants of the COST Action CA16119 (In vitro 3-D total cell guidance and fitness).
This study was partly funded by the Schweizerischer Nationalfonds (SNF, 31003A_159734). The publication of this article was funded by the Open Access Fund of the Leibniz Association.
Availability of data and materials
The datasets supporting the conclusions of this article are included within the article and its additional files.
VvdW performed the qRT-PCR experiments, analysed the data and drafted the manuscript. SC carried out the cell and embryo culture experiments, helped with data analysis and drafting the manuscript. SB supervised the data analysis and bioinformatics. SU and JS conceived of the study, supervised the experiments and finalized the manuscript. All authors read, edited and approved the final manuscript.
Ethics approval and consent to participate
Not applicable. Animal materials used in this study were by-products collected in a local commercial slaughterhouse.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.
- Fazeli A. Maternal communication with gametes and embryos. Theriogenology. 2008;70:1182–7.View ArticlePubMedGoogle Scholar
- Danesh MS, Sharbati J, Einspanier R, Gabler C. mRNA expression pattern of selected candidate genes differs in bovine oviductal epithelial cells in vitro compared with the in vivo state and during cell culture passages. Reprod Biol Endocrinol. 2016;14:44.View ArticleGoogle Scholar
- Ulbrich SE, Zitta K, Hiendleder S, Wolf E. In vitro systems for intercepting early embryo-maternal cross-talk in the bovine oviduct. Theriogenology. 2010;73:802–16.View ArticlePubMedGoogle Scholar
- Chen S, Palma-Vera SE, Langhammer M, Galuska SP, Braun BC, Krause E, Lucas-Hahn A, Schoen J. An air-liquid interphase approach for modeling the early embryo-maternal contact zone. Sci Rep. 2017;7:42298.View ArticlePubMedPubMed CentralGoogle Scholar
- Alm H, Torner H, Kanitz W, Roschlau K. Influence of oocyte recovery method, in vitro fertilization method and serum source on embryonic development of in vitro matured bovine oocytes. Archiv Fur Tierzucht-Archives of Animal Breeding. 2008;51:224–34.View ArticleGoogle Scholar
- Chen S, Einspanier R, Schoen J. Long-term culture of primary porcine oviduct epithelial cells: validation of a comprehensive in vitro model for reproductive science. Theriogenology. 2013;80:862–9.View ArticlePubMedGoogle Scholar
- Nestorov P, Hotz HR, Liu Z, Peters AH. Dynamic expression of chromatin modifiers during developmental transitions in mouse preimplantation embryos. Sci Rep. 2015;5:14347.View ArticlePubMedPubMed CentralGoogle Scholar
- Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, De Paepe A, Speleman F. Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol. 2002;3:RESEARCH0034.View ArticlePubMedPubMed CentralGoogle Scholar
- Hackett AJ, Durnford R, Mapletoft RJ, Marcus GJ. Location and status of embryos in the genital-tract of Superovulated cows 4 to 6 days after insemination. Theriogenology. 1993;40:1147–53.View ArticleGoogle Scholar
- Lamy J, Labas V, Harichaux G, Tsikis G, Mermillod P, Saint-Dizier M. Regulation of the bovine oviductal fluid proteome. Reproduction. 2016;152:629–44.View ArticlePubMedGoogle Scholar
- Graf A, Krebs S, Zakhartchenko V, Schwalb B, Blum H, Wolf E. Fine mapping of genome activation in bovine embryos by RNA sequencing. Proc Natl Acad Sci U S A. 2014;111:4139–44.View ArticlePubMedPubMed CentralGoogle Scholar
- Azzi S, Sas TC, Koudou Y, Le Bouc Y, Souberbielle JC, Dargent-Molina P, Netchine I, Charles MA. Degree of methylation of ZAC1 (PLAGL1) is associated with prenatal and post-natal growth in healthy infants of the EDEN mother child cohort. Epigenetics. 2014;9:338–45.View ArticlePubMedGoogle Scholar
- Hajian M, Hosseini SM, Ostadhosseini S, Nasr-Esfahani MH. Targeting the transforming growth factor-beta signaling during pre-implantation development in embryos of cattle, sheep and goats. Growth Factors. 2016;34:141–8.View ArticlePubMedGoogle Scholar
- Heo JY, Jing K, Song KS, Seo KS, Park JH, Kim JS, Jung YJ, Hur GM, Jo DY, Kweon GR, et al. Downregulation of APE1/Ref-1 is involved in the senescence of mesenchymal stem cells. Stem Cells. 2009;27:1455–62.View ArticlePubMedGoogle Scholar
- Kobayashi Y, Yamamoto Y, Kageyama S, Hirayama H, Kimura K, Okuda K. Regulation of bovine oviductal NO synthesis by follicular steroids and prostaglandins. Reproduction. 2016;151:577–87.View ArticlePubMedGoogle Scholar
- Rizos D, Lonergan P, Boland MP, Arroyo-Garcia R, Pintado B, de la Fuente J, Gutierrez-Adan A. Analysis of differential messenger RNA expression between bovine blastocysts produced in different culture systems: implications for blastocyst quality. Biol Reprod. 2002;66:589–95.View ArticlePubMedGoogle Scholar
- Tesfaye D, Lonergan P, Hoelker M, Rings F, Nganvongpanit K, Havlicek V, Besenfelder U, Jennen D, Tholen E, Schellander K. Suppression of connexin 43 and E-cadherin transcripts in in vitro derived bovine embryos following culture in vitro or in vivo in the homologous bovine oviduct. Mol Reprod Dev. 2007;74:978–88.View ArticlePubMedGoogle Scholar
- Mitra S, Izumi T, Boldogh I, Bhakat KK, Chattopadhyay R, Szczesny B. Intracellular trafficking and regulation of mammalian AP-endonuclease 1 (APE1), an essential DNA repair protein. DNA Repair (Amst). 2007;6:461–9.View ArticleGoogle Scholar
- Son MY, Choi H, Han YM, Cho YS. Unveiling the critical role of REX1 in the regulation of human stem cell pluripotency. Stem Cells. 2013;31:2374–87.View ArticlePubMedGoogle Scholar
- Shamseddine AA, Airola MV, Hannun YA. Roles and regulation of neutral sphingomyelinase-2 in cellular and pathological processes. Adv Biol Regul. 2015;57:24–41.View ArticlePubMedGoogle Scholar
- Lagow E, DeSouza MM, Carson DD. Mammalian reproductive tract mucins. Hum Reprod Update. 1999;5:280–92.View ArticlePubMedGoogle Scholar
- Clemente M, Lopez-Vidriero I, O'Gaora P, Mehta JP, Forde N, Gutierrez-Adan A, Lonergan P, Rizos D. Transcriptome changes at the initiation of elongation in the bovine conceptus. Biol Reprod. 2011;85:285–95.View ArticlePubMedGoogle Scholar
- Imakawa K, Tamura K, Lee RS, Ji Y, Kogo H, Sakai S, Christenson RK. Temporal expression of type I interferon receptor in the peri-implantation ovine extra-embryonic membranes: demonstration that human IFNalpha can bind to this receptor. Endocr J. 2002;49:195–205.View ArticlePubMedGoogle Scholar
- Larue L, Ohsugi M, Hirchenhain J, Kemler R. E-cadherin null mutant embryos fail to form a trophectoderm epithelium. Proc Natl Acad Sci U S A. 1994;91:8263–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Marchand M, Horcajadas JA, Esteban FJ, McElroy SL, Fisher SJ, Giudice LC. Transcriptomic signature of trophoblast differentiation in a human embryonic stem cell model. Biol Reprod. 2011;84:1258–71.View ArticlePubMedGoogle Scholar
- Smith SL, Everts RE, Tian XC, Du F, Sung LY, Rodriguez-Zas SL, Jeong BS, Renard JP, Lewin HA, Yang X. Global gene expression profiles reveal significant nuclear reprogramming by the blastocyst stage after cloning. Proc Natl Acad Sci U S A. 2005;102:17582–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Cagnone G, Sirard MA. The impact of exposure to serum lipids during in vitro culture on the transcriptome of bovine blastocysts. Theriogenology. 2014;81:712–22. e711-713View ArticlePubMedGoogle Scholar
- Lonergan P, Rizos D, Gutierrez-Adan A, Moreira PM, Pintado B, de la Fuente J, Boland MP. Temporal divergence in the pattern of messenger RNA expression in bovine embryos cultured from the zygote to blastocyst stage in vitro or in vivo. Biol Reprod. 2003;69:1424–31.View ArticlePubMedGoogle Scholar
- Vergnes L, Chin RG, de Aguiar Vallim T, Fong LG, Osborne TF, Young SG, Reue K. SREBP-2-deficient and hypomorphic mice reveal roles for SREBP-2 in embryonic development and SREBP-1c expression. J Lipid Res. 2016;57:410–21.View ArticlePubMedPubMed CentralGoogle Scholar