- Open Access
Genome-wide microRNA expression profiling in placentae from frozen-thawed blastocyst transfer
Clinical Epigenetics volume 9, Article number: 79 (2017)
Frozen-thawed embryo transfer (FET) is increasingly available for the improvement of the success rate of assisted reproductive technologies other than fresh embryo transfer (ET). There have been numerous findings that FET provides better obstetric and perinatal outcomes. However, the birth weight of infants conceived using FET is heavier than that of those conceived via ET. In addition, some reports have suggested that FET is associated with perinatal diseases such as placenta accreta and pregnancy-induced hypertension (PIH).
In this study, we compared the microRNA (miRNA) expression profiles in term placentae derived from FET, ET, and spontaneous pregnancy (SP). We identified four miRNAs, miR-130a-3p, miR-149-5p, miR-423-5p, and miR-487b-3p, that were significantly downregulated in FET placentae compared with those from SP and ET. We found that DNA methylation of MEG3-DMR, not but IG-DMR, was associated with miRNA expression of the DLK1-DIO3 imprinted domain in the human placenta. In functional analyses, GO terms and signaling pathways related to positive regulation of gene expression, growth, development, cell migration, and type II diabetes mellitus (T2DM) were enriched.
This study supports the hypothesis that the process of FET may increase exposure of epigenome to external influences.
Recently, due to dramatic increases in success rates with frozen-thawed embryo transfer (FET), obstetric and perinatal outcomes have improved [1, 2]. A large Japanese registry study indicated that FET reduced the incidences of small for gestational age (SGA), low birthweight (LBW), and prematurity when compared with fresh embryo transfer (ET) . In other large studies and comprehensive meta-analyses, FET was associated with increased success rates and decreased obstetric and perinatal complications compared to ET [1, 4, 5]. However, singleton FET is adversely associated with increased risks of macrosomia and large for gestational age (LGA) and higher incidences of placenta accreta and pregnancy-induced hypertension (PIH) [3, 6].
For mice, it has been reported that FET affects the expression of imprinted genes and DNA methylation in both the fetus and placenta . With assisted reproductive technologies (ART), including FET, major epigenetic reprogramming events occur and may expose the epigenome to external influences, preventing proper embryogenesis and placentogenesis [8,9,10]. In humans, epigenetic alterations of the placenta in utero may increase the risk of developing various diseases later in life such as diabetes mellitus (DM) .
MicroRNAs (miRNAs), approximately 22 nucleotides in length, are evolutionarily well-conserved non-coding RNAs that regulate gene expression by binding specific sites of 3′ or 5′ untranslated regions of target genes. They promote target mRNA degradation and translational inhibition, and are considered a mode of epigenetic regulation . Abnormal expression of miRNAs has been observed in pregnancy-associated diseases and recurrent pregnancy loss [13, 14]. Several pregnancy-related miRNAs cluster in the large human genomic regions located in imprinted chromosomes 14q32 (C14MC) and 19q13.42 (C19MC) [15, 16]. Parent-of-origin-specific methylation imprints are established during oogenesis or spermatogenesis and stably maintained through fertilization and subsequent embryonic and placental development in mammals .
In this study we, for the first time, performed genome-wide miRNA analysis of term placentae derived from FET, ET, and spontaneous pregnancy (SP). We found that expression of miR-130a-3p, miR-149-5p, miR-423-5p, and miR-487b-3p was decreased in the term placentae derived from FET compared with those derived from ET or SP. Functional analyses suggest that these miRNAs are involved in an increase of neonatal birthweight and potentiality lead to diseases later in life such as DM. This study supports our hypothesis that the process of FET may increase exposure of the epigenome to external influences.
Pregnant women who conceived by SP, FET or ET were recruited from one private hospital. Patient information is shown in Additional file 1: Table S1. There were no severe perinatal complications. Controlled ovarian stimulation was performed with combination of a GnRH agonist, letrozole/GnRH agonist, GnRH antagonist or clomiphene citrate, and recombinant FSH, hMG, and hCG. All inseminations for FET and ET were performed by intracytoplasmic sperm injection (ICSI). After ICSI, embryos were cultured in Quinn’s Advantage Medium at 37 °C for 5 or 6 days in an atmosphere of 6% CO2, 5% O2, and 89% N2 under humidified conditions. Blastocysts were vitrified using a Cryotop (Kitazato BioPharma, Fuji, Japan) and then were immersed and stored in liquid nitrogen (LN2) until embryo transfer. To warm vitrified blastocysts, the Cryotop was removed from the LN2 and instantly immersed into 3 ml of the medium containing 1.0 M sucrose at 37 °C for 1 min. All patients who underwent FET used the hormone replacement cycle (HRC). Under the HRC, patients were prepared for ET using conjugated estrogens or transdermal estradiol. Placentae were collected from singleton pregnancies after term delivery. Two 1-cm3 cubes of the placental attachment of the umbilical cord were excised from term placenta and were washed in PBS at least five times to remove blood cells. For RNA extraction, one was homogenized in 5 ml of Isogen (Nippon Gene, Toyama, Japan) and was preserved immediately. For DNA extraction, the other was immediately transferred into a 1.5 ml tube. All samples were stored at −80 °C until use. We used adjusted birthweight by subtracting the value of the birthweight from the value of the fetal growth curve considering the Japanese gestational age-, gender-, and parity-specific growth standard .
RNA extraction and DNA extraction
Total RNA was isolated from placental tissue using an miRNeasy Mini Kit (Qiagen, Valencia, CA, USA). RNA was quantified using a NanoDrop ND-1000 spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA). The integrity of the RNA was assessed on an Agilent TapeStation 2200 system (Agilent Technologies, Santa Clara, CA, USA). Purified RNA aliquots were stored at −80 °C until further processing. Genomic DNA was obtained from placental tissue using the standard extraction protocol and stored at −30 °C until use.
For microarray analysis, we randomly picked four representative placentae (two male and two female infants) derived from each group. The Agilent SurePrint G3 Human miRNA 8 × 60k (release 21.0) microarray platform for 2588 mature miRNAs was used. Total RNA (100 ng) was labeled and hybridized. The microarray slides were scanned on an Agilent G4900DA SureScan Microarray Scanner System and data were extracted from the scanned images using Agilent Feature Extraction software version 188.8.131.52. Raw microarray data were deposited in the NCBI Gene Expression Omnibus repository (GSE85270). Subsequent analysis was carried out in R (http://www.r-project.org) using functions and packages collected in the Bioconductor project . Raw microarray data were filtered by detection in all samples. Selected miRNA values (n = 471) were transformed logarithmically and were normalized by the quantile method.
Quantitative reverse transcription polymerase chain reaction
Quantitative RT-PCR (qRT-PCR) analysis was performed using TaqMan MicroRNA Assays (Applied Biosystems, Foster City, CA, USA). Briefly, 5 ng of total RNA was used for the synthesis of the first strand cDNA using a TaqMan MicroRNA Reverse Transcription Kit (Applied Biosystems; assay IDs are listed in Additional file 2: Table S2). qRT-PCR was performed on a StepOne Real-time PCR System (Applied Biosystems) using TaqMan Fast Advanced Master Mix (Applied Biosystems). The fold change for miRNA, relative to RNU44 snoRNA, was calculated using the 2−ΔΔCt method . All RT-PCR reactions were performed in duplicate.
miRNA-predicted target genes and gene ontology (GO) and pathway analyses
Putative target genes of differentially expressed miRNAs were predicted using six available bioinformatics algorithms (miRWalk, DIANA-microT4, miRanda, miRDB, PICTAR2, and TargetScan) with miRWalk2.0 . To reduce the number of false positives, only target genes that were predicted by at least four of the six programs were selected and used for further investigation. The biological annotation and the potential pathway were analyzed using DAVID, version 6.7  and KEGG pathway enrichment analysis , respectively.
Combined bisulfite restriction analysis coupled with the Agilent 2200 TapeStation platform (Bio-COBRA)
Bisulfite treatment of genomic DNA was carried out using an EZ DNA methylation-Gold Kit (Zymo Research, Orange, CA, USA). PCR was performed using TaKaRa EpiTaq HS (Takara Biomedical, Otsu, Japan). The PCR reaction was performed in a volume of 20 μl. Primer and PCR conditions are shown in Additional file 3: Table S3. Of the 20 μl amplicons, 10 μl was digested with 5 units of BstUI, HinfI (New England Biolabs, Hitchin, UK), and TaqI (Nippon Gene) for IG-DMR, MEG3-DMR, and C19MC-DMR for 16 h, respectively. The digested amplicons were separated and visualized using an Agilent 2200 TapeStation system.
Statistical analyses were performed using R or JMP Pro 11.2.0 (SAS Institute, Cary, NC, USA). Statistical significance among the clinical characteristics, miRNA expression levels quantified by qRT-PCR and DNA methylation levels of the three groups was determined using the Steel-Dwass test at p < 0.05. Significance in microarray analysis among the three groups was determined by one-way ANOVA with the Tukey-Kramer post hoc test. Correlations were calculated using the Pearson correlation coefficient. Fisher’s exact test with the Benjamini-Hochberg procedure was used to calculate the p values for enriched gene ontology (GO) terms in the DAVID analysis and pathways in the KEGG analysis.
We collected 108 placentae, including 28 SP, 64 FET, and 16 ET placentae. The clinical characteristics of these pregnancies are summarized in Table 1. Adjusted birthweight (p = 0.0092) and placental weight (p = 0.0237) in the FET group were significantly heavier than in the SP group but not the ET group. Mothers were older in the FET (p < 0.0001) and ET (p = 0.0111) groups than in the SP group. The LGA rate for FET was significantly different (p = 0.0150) from SP, but not ET.
miRNA microarray analysis
To investigate whether placental miRNAs were associated with the characteristics of FET, comprehensive microarray analysis was performed using total RNAs from four placentae each derived from SP, FET, and ET (Fig. 1). Among the 471 human miRNAs detected by the microarray, which targets 2588 mature miRNAs, miRNAs with 2-fold changes between each 2 groups and p < 0.05 were identified as differentially expressed. We identified 94 miRNAs significant differentially expressed (46 upregulated and 48 downregulated), between FET and SP placentae, and 40 (18 upregulated and 22 downregulated), between FET and ET placentae. However, we found none between SP and ET placentae. We constructed a heat-map of all differentially expressed miRNAs among the 3 groups and 12 samples separated into 2 groups (4 FETs/1 ET, and 4 SPs/3 ETs) based on the miRNA expression patterns (Fig. 1d). A total of 39 miRNAs, 18 upregulated and 21 downregulated, were commonly differentially expressed in the FET placentae compared with those from SP and ET (Table 2). Notably, among these 95 miRNAs, 18 were located in 3 imprinted regions. Thirteen of them belonged to C19MC, 4 to C14MC, and 1 to the IGF2 gene (Additional file 4: Table S4).
miRNA expression validation by qRT-PCR
To confirm the microarray results, we conducted qRT-PCR analysis of the expression levels of 21 common differentially expressed miRNAs that had raw signal data with values higher than 50 in all samples of the microarray using 96 placentae. These miRNAs consisted of 6 of the 18 that were upregulated and 15 of the 21 downregulated ones in microarray analysis of the FET placentae. Of these miRNAs, four (miR-130a-3p, miR-149-5p, miR-423-5p, and miR-487b-3p) were significantly downregulated in FET placentae compared with those from SP and ET. miR-193b-3p was significantly downregulated in FET and ET placentae compared with SP placentae. Nine (miR-125a-5p, miR-224-3p, miR-331-3p, miR-365a-3p, miR-495-3p, miR-518b, miR-518f-3p, miR-543, and miR-7977) were significantly downregulated in FET placentae compared with SP, but not ET (Fig. 2 and Additional file 5: Figure S1).
To investigate the relations between miRNA expression and adjusted birthweight and placental weight, we calculated the correlation coefficients of the 21 miRNAs validated by qRT-PCR. There were weak inverse correlations between miR-193b-3p expression and not only placental weight (R = −0.29, p = 0.0047) but also adjusted birthweight (R = −0.29, p = 0.0050). For the other 20 miRNAs, there was no significant correlation between miRNA expression and adjusted birthweight or placental weight (Additional file 6: Figure S2 and Additional file 7: Table S5). In four FET placenta-specific miRNAs, there was no significant correlation between miRNA expression and maternal age, whereas there were weak correlations between expression levels of miR-125a-5p, miR-224-3p, miR-331-3p, miR-365a-3p, miR-518b, miR-518f-3p, and miR-543 and maternal age (Additional file 7: Table S5).
miRNA-predicted target genes and functional analysis
To investigate the pathways in which miRNAs were involved, we searched for potential target genes of the four differentially expressed miRNAs in the FET placentae using six prediction programs, GO terms, and signaling pathways related to the target genes. The target predictions for the four miRNAs showed 4243 putative target sites in 2686 genes. GO enrichment analysis based on the biological process showed that these predicted target genes were clustered into 553 GO terms, among which 161 were significant (corrected p value < 0.05) (Additional file 8: Table S6). The top 20 enrichment GO terms are shown in Table 3. Among these, six were associated with positive regulation of transcription, four with the metabolic system, four with the biosynthetic process, and cell migration and cell motion GO terms were enriched. KEGG analysis of these predicted target genes revealed 52 pathways (Additional file 9: Table S7), including 32 significant ones. The top 20 enrichment pathways are shown in Table 4. The signaling, endocrine, gap junction, and focal adhesion pathways, as well as pathways related to some cancers, and the type II DM (T2DM) pathway, were enriched (corrected p value = 0.047, Additional file 9: Table S7).
DNA methylation analysis
We identified 17 differentially expressed miRNAs in imprinted C14MC and C19MC among the three groups (Fig. 1d and Additional file 4: Table S4). To determine whether the downregulation of these miRNAs could be explained by the DNA methylation of imprinted differentially methylated regions (DMRs), we examined the DNA methylation of IG-DMR and MEG3-DMR, located in C14MC and C19MC-DMR, using Bio-COBRA (Fig. 3). Although the median methylation levels of IG-DMR and C19MC-DMR were not significantly different among the three groups, that of MEG3-DMR in the FET placentae was significantly increased 4.86% (p = 0.0053) compared to the SP placentae, but not the ET ones (increased 4.85%, p = 0.0979).
With significant improvements of vitrification, the number of FET procedures and success rate have increased in Japan. However, the safety of FET with respect to perinatal complications and the offspring remains unknown. miRNAs are epigenetic mechanisms for the prediction and detection of placenta-mediated pregnancy complications . In this study, we performed genome-wide miRNA analysis comparing term placentae derived from FET, ET, and SP and confirmed four FET placenta-specific miRNAs.
In microarray screening, since the difference between FET and SP was larger than that between FET and SP and there was none between ET and SP, the FET procedure might affect miRNA expression more than ET. The 39 miRNAs commonly differentially expressed in the FET placentae might be involved in conditions such as overgrowth and perinatal complications.
Of the 21 miRNAs that were commonly differentially expressed in the FET placentae compared with those from SP and ET, 14 showed significant differential expression in FET by qRT-PCR. Four of these 14 miRNAs were specifically downregulated in FET. They might contribute to increased birthweight and placental weight; however, their underlying mechanism is unclear. The downregulation of these four miRNAs might have been affected by the FET procedure, not maternal age, although maternal age in the FET and ET groups was higher than for SP and the expression levels of miRNAs decrease with aging . We identified GO terms and pathways related to the four FET-specific miRNAs, as potential causal pathways in the pathogenesis of overgrowth and placenta accreta. Although pathways related to some cancers were enriched, the genes predicted by these four miRNAs might be involved in growth, but not related to cancers themselves. Although there were no differences in the outcomes except for birthweight and placental weight among the three groups, these GO terms and pathways might cause perinatal complications easily.
That the other nine miRNAs showed significantly different expression between FET and SP, but not ET, might be due to the fact that the difference of the miRNA expression levels between FET and SP was greater than between ET and SP. Seven of the nine miRNAs that showed significantly different expression only between FET and SP might be affected by maternal age rather than the FET procedure. The other two might be involved in overgrowth, although they tended to be decreased in FET compared with ET and showed no significant difference.
miR-193b-3p was decreased in FET and ET placentae compared to those from SP and was related to overgrowth. The decreased expression of miR-193b-3p might be due to procedures common to FET and ET, including ovulation-induction methods, insemination, culture system, and embryo transfer. The lack of a significant difference in birthweight and placental weight between ET and SP might be due to factors that compensate for the influence of miR-193b-3p and complement the increases of fetal and placental growth in ET. For 20 miRNAs, each miRNA alone might not affect overgrowth in FET, nevertheless overgrowth might be a synergistic effect of a miRNAs combination. In fact, miRNA binding is partially complementary to multiple mRNA sequences and some miRNAs downregulate large numbers of target mRNAs . In other words, expression of one mRNA is regulated by multiple miRNAs. Therefore, it seems that miRNAs and mRNAs form a complex gene expression network. Moreover, though all human protein-coding genes might not be regulated by miRNAs, more than 30% of human protein-coding genes contain at least one miRNA-binding site . For these reasons, there might be no good correlation between the expression levels of the four miRNAs and the birthweight and placental weight in this study.
In microarray analysis, 13 detectable miRNAs of C19MC and 4 of C14MC were decreased in the FET placentae. Since some C19MC and C14MC miRNAs are associated with pregnancy-related complications , the miRNAs we detected might also be. The methylation in MEG3-DMR was increased in FET placentae and might play a role in the regulation of C14MC miRNA expression. On the other hand, the methylation of C19MC miRNAs did not differ among the groups and other mechanisms such as histone modifications or unknown promoters might regulate C19MC miRNAs. The methylation of MEG3-DMR was increased in FET and ET placentae compared to those from SP, but the lack of a significant difference between ET and SP suggested that the effect for not only miRNA expression but also DNA methylation of FET might be greater than for ET. Many studies have suggested that there is an increased incidence of rare imprinting disorders, including Beckwith–Wiedemann syndrome (BWS)  and Silver-Russell syndrome (SRS) , associated with ART . MEG3-DMR, which is the location of the MEG3 promoter, encodes large non-coding RNA, as well as H19-DMR and KCNQ1OT1-DMR involved in SRS and BWS, respectively, might be easily affected by ART. Meg3 is involved in growth and development although Meg3-mutant mice exhibit very unique inheritance . Those results and ours may suggest that some imprinted genes or domains are labile and readily changeable.
In animal studies, the in vitro culture methods and compositions of culture media might induce male large offspring syndrome in ruminants [32, 33]. The effects of ART, including in vitro culture, can be maintained, or more pronounced even after the freezing and thawing of embryos, as embryo vitrification aggravated the loss of methylation of H19-DMR in an IVF group with vitrification in mice . Both ovulation-induction treatment and embryo culture increase the perturbation of genomic imprinting [34, 35]. Recently, a randomized study showed that different media for culture of human in vitro fertilization (IVF) embryos affect the birthweights of newborns [36, 37]. The developmental origins of health and disease paradigm postulate that suboptimal growth early in life can program changes that have a lifelong effect on health, increasing the risks for various diseases [38, 39]. Interestingly, not only LBW but also high birthweight infants are at increased risk for future development of type 2 diabetes mellitus (T2DM) . We found a decrease of expression of placental miR-487b-3p, which is located in C14MC, and the pathway associated with T2DM. Recently, islets . Abnormalities occurring in the placenta such as placental insufficiency, placental infarction, and placental angiogenesis might indirectly have adverse effects on the fetus [42,43,44]. Importantly, the role of environmental epigenetics in FET suggests the possibility of leading to diseases later in life such as DM .
There are some limitations to this study. First, this was a study of only Japanese term placentae with no involvement of severe perinatal complications. Second, we could collect only 16 ET placentae. Third, we performed genome-wide miRNA assay comparing three groups. Although the microarray technologies for high-throughput analysis are well established, they cannot detect unknown miRNAs. However, a strength of this study is that both patients with FET and fresh ET underwent the same ICSI protocol at one ART institute in spite of the different backgrounds of the patients.
In conclusion, our miRNA data support the hypothesis that major epigenetic events taking place during early embryogenesis and the process of FET may increase exposure of the epigenome to external influences. This could lead to being disproportionately large at birth and placenta-mediated pregnancy complications linked to perinatal and adult diseases. Our knowledge that the epigenomes of gametes and newly fertilized embryos are susceptible to environmentally induced epigenetic changes has particularly important implications as changes in lifestyle and modes of reproduction may have long-term implications for human health that are not yet fully appreciated.
Assisted reproductive technology
Chromosome 14 microRNA cluster
Chromosome 19 microRNA cluster
Combined bisulfite restriction analysis
Database for Annotation, Visualization and Integrated Discovery
Frozen-thawed embryo transfer
Follicular stimulating hormone
Human chorionic gonadotropin
Human menopausal gonadotropin
Hormone replacement cycle
Intracytoplasmic sperm injection
In vitro fertilization
Kyoto Encyclopedia of Genes and Genomes
Large for gestational age
- LN2 :
Mitogen-activated protein kinase
Recombinant follicular stimulating hormone
Reverse transcription polymerase chain reaction
Small for gestational age
Type 2 diabetes mellitus
Vascular endothelial growth factor
Maheshwari A, Pandey S, Shetty A, Hamilton M, Bhattacharya S. Obstetric and perinatal outcomes in singleton pregnancies resulting from the transfer of frozen thawed versus fresh embryos generated through in vitro fertilization treatment: a systematic review and meta-analysis. Fertil Steril. 2012;98(2):368–77. e361-369
Shapiro BS, Daneshmand ST, Garner FC, Aguirre M, Hudson C. Clinical rationale for cryopreservation of entire embryo cohorts in lieu of fresh transfer. Fertil Steril. 2014;102(1):3–9.
Ishihara O, Araki R, Kuwahara A, Itakura A, Saito H, Adamson GD. Impact of frozen-thawed single-blastocyst transfer on maternal and neonatal outcome: an analysis of 277,042 single-embryo transfer cycles from 2008 to 2010 in Japan. Fertil Steril. 2014;101(1):128–33.
Shih W, Rushford DD, Bourne H, Garrett C, McBain JC, Healy DL, Baker HW. Factors affecting low birthweight after assisted reproduction technology: difference between transfer of fresh and cryopreserved embryos suggests an adverse effect of oocyte collection. Hum Reprod. 2008;23(7):1644–53.
Pinborg A, Loft A, Aaris Henningsen AK, Rasmussen S, Andersen AN. Infant outcome of 957 singletons born after frozen embryo replacement: the Danish National Cohort Study 1995–2006. Fertil Steril. 2010;94(4):1320–7.
Barsky M, St Marie P, Rahil T, Markenson GR, Sites CK. Are perinatal outcomes affected by blastocyst vitrification and warming? Am J Obstet Gynecol. 2016;215:603.e1–5.
Wang Z, Xu L, He F. Embryo vitrification affects the methylation of the H19/Igf2 differentially methylated domain and the expression of H19 and Igf2. Fertil Steril. 2010;93(8):2729–33.
Guo H, Zhu P, Yan L, Li R, Hu B, Lian Y, Yan J, Ren X, Lin S, Li J, et al. The DNA methylation landscape of human early embryos. Nature. 2014;511(7511):606–10.
Smith ZD, Chan MM, Humm KC, Karnik R, Mekhoubad S, Regev A, Eggan K, Meissner A. DNA methylation dynamics of the human preimplantation embryo. Nature. 2014;511(7511):611–5.
Okae H, Chiba H, Hiura H, Hamada H, Sato A, Utsunomiya T, Kikuchi H, Yoshida H, Tanaka A, Suyama M, et al. Genome-wide analysis of DNA methylation dynamics during early human development. PLoS Genet. 2014;10(12):e1004868.
Rinaudo P, Wang E. Fetal programming and metabolic syndrome. Annu Rev Physiol. 2012;74:107–30.
Ha M, Kim VN. Regulation of microRNA biogenesis. Nat Rev Mol Cell Biol. 2014;15(8):509–24.
Santamaria X, Taylor H. MicroRNA and gynecological reproductive diseases. Fertil Steril. 2014;101(6):1545–51.
Mouillet JF, Ouyang Y, Coyne CB, Sadovsky Y. MicroRNAs in placental health and disease. Am J Obstet Gynecol. 2015;213(4 Suppl):S163–72.
Gu Y, Sun J, Groome LJ, Wang Y. Differential miRNA expression profiles between the first and third trimester human placentas. Am J Physiol Endocrinol Metab. 2013;304(8):E836–43.
Morales-Prieto DM, Ospina-Prieto S, Chaiwangyen W, Schoenleben M, Markert UR. Pregnancy-associated miRNA-clusters. J Reprod Immunol. 2013;97(1):51–61.
Barlow DP, Bartolomei MS. Genomic imprinting in mammals. Cold Spring Harbor Perspec Biol. 2014;6(2):a018382.
Itabashi K, Miura F, Uehara R, Nakamura Y. New Japanese neonatal anthropometric charts for gestational age at birth. Pediatr Int. 2014;56(5):702–8.
Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, Ellis B, Gautier L, Ge Y, Gentry J, et al. Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004;5(10):R80.
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(−Delta Delta C(T)) Method. Methods. 2001;25(4):402–8.
Dweep H, Gretz N. miRWalk2.0. a comprehensive atlas of microRNA-target interactions. Nat Methods. 2015;12(8):697.
Huang da W, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4(1):44–57.
Altermann E, Klaenhammer TR. PathwayVoyager: pathway mapping using the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. BMC Genomics. 2005;6:60.
Fu G, Brkic J, Hayder H, Peng C. MicroRNAs in Human Placental Development and Pregnancy Complications. Int J Mol Sci. 2013;14(3):5519–44.
Noren Hooten N, Fitzpatrick M, Wood WH 3rd, De S, Ejiogu N, Zhang Y, Mattison JA, Becker KG, Zonderman AB, Evans MK. Age-related changes in microRNA levels in serum. Aging. 2013;5(10):725–40.
Lim LP, Lau NC, Garrett-Engele P, Grimson A, Schelter JM, Castle J, Bartel DP, Linsley PS, Johnson JM. Microarray analysis shows that some microRNAs downregulate large numbers of target mRNAs. Nature. 2005;433(7027):769–73.
Lewis BP, Burge CB, Bartel DP. Conserved seed pairing, often flanked by adenosines, indicates that thousands of human genes are microRNA targets. Cell. 2005;120(1):15–20.
Choufani S, Shuman C, Weksberg R. Beckwith-Wiedemann syndrome. Am J Med Genet C Sem Med Genet. 2010;154C(3):343–54.
Eggermann T. Russell-Silver syndrome. Am J Med Genet C Sem Med Genet. 2010;154C(3):355–64.
Kitamura A, Miyauchi N, Hamada H, Hiura H, Chiba H, Okae H, Sato A, John RM, Arima T. Epigenetic alterations in sperm associated with male infertility. Congenit Anom. 2015;55(3):133–44.
Takahashi N, Okamoto A, Kobayashi R, Shirai M, Obata Y, Ogawa H, Sotomaru Y, Kono T. Deletion of Gtl2, imprinted non-coding RNA, with its differentially methylated region induces lethal parent-origin-dependent defects in mice. Hum Mol Genet. 2009;18(10):1879–88.
Young LE, Sinclair KD, Wilmut I. Large offspring syndrome in cattle and sheep. Rev Reprod. 1998;3(3):155–63.
Chen Z, Robbins KM, Wells KD, Rivera RM. Large offspring syndrome: a bovine model for the human loss-of-imprinting overgrowth syndrome Beckwith-Wiedemann. Epigenetics. 2013;8(6):591–601.
Market-Velker BA, Fernandes AD, Mann MR. Side-by-side comparison of five commercial media systems in a mouse model: suboptimal in vitro culture interferes with imprint maintenance. Biol Reprod. 2010;83(6):938–50.
Market-Velker BA, Zhang L, Magri LS, Bonvissuto AC, Mann MR. Dual effects of superovulation: loss of maternal and paternal imprinted methylation in a dose-dependent manner. Hum Mol Genet. 2010;19(1):36–51.
Dumoulin JC, Land JA, Van Montfoort AP, Nelissen EC, Coonen E, Derhaag JG, Schreurs IL, Dunselman GA, Kester AD, Geraedts JP, et al. Effect of in vitro culture of human embryos on birthweight of newborns. Hum Reprod. 2010;25(3):605–12.
Nelissen EC, Van Montfoort AP, Coonen E, Derhaag JG, Geraedts JP, Smits LJ, Land JA, Evers JL, Dumoulin JC. Further evidence that culture media affect perinatal outcome: findings after transfer of fresh and cryopreserved embryos. Hum Reprod. 2012;27(7):1966–76.
Bateson P, Barker D, Clutton-Brock T, Deb D, D'Udine B, Foley RA, Gluckman P, Godfrey K, Kirkwood T, Lahr MM, et al. Developmental plasticity and human health. Nature. 2004;430(6998):419–21.
Gluckman PD, Hanson MA. Living with the past: evolution, development, and patterns of disease. Science. 2004;305(5691):1733–6.
Wei JN, Sung FC, Li CY, Chang CH, Lin RS, Lin CC, Chiang CC, Chuang LM. Low birth weight and high birth weight infants are both at an increased risk to have type 2 diabetes among schoolchildren in taiwan. Diabetes Care. 2003;26(2):343–8.
Kameswaran V, Bramswig NC, McKenna LB, Penn M, Schug J, Hand NJ, Chen Y, Choi I, Vourekas A, Won KJ, et al. Epigenetic regulation of the DLK1-MEG3 microRNA cluster in human type 2 diabetic islets. Cell Metab. 2014;19(1):135–45.
Pardi G, Cetin I, Marconi AM, Lanfranchi A, Bozzetti P, Ferrazzi E, Buscaglia M, Battaglia FC. Diagnostic value of blood sampling in fetuses with growth retardation. N Engl J Med. 1993;328(10):692–6.
Naeye RL. Placental infarction leading to fetal or neonatal death. A prospective study. Obstet Gynecol. 1977;50(5):583–8.
Ahmed A, Perkins J. Angiogenesis and intrauterine growth restriction. Bailliere Best Pract Res Clin Obstet Gynaecol. 2000;14(6):981–98.
Kyle UG, Pichard C. The Dutch famine of 1944–1945: a pathophysiological model of long-term consequences of wasting disease. Curr Opin Clin Nutr Metab Care. 2006;9(4):388–94.
Gardner DK, Schoolcraft WB. In vitro culture of human blastocyst. In: Jansen R, Carnforth MD, editors. Towards Reproductive Certainty: Infertility and Genetics Beyond. UK: Parthenon Press; 1999. p. 377–88.
We would like to thank all the members of our laboratory for their technical assistance, support and valuable suggestions. In particular, we thank Dr. R. M. John (Cardiff School of Biosciences) for comments on the manuscript.
This work was supported by Grants-in-Aid for Scientific Research (KAKENHI) (21028003, 23013003, 23390385), and the Takeda Foundation (TA).
Availability of data and materials
The datasets generated and analyzed during the current study are available in the NCBI Gene Expression Omnibus repository, GEO Series accession number: GSE85270.
Ethics approval and consent to participate
This study was approved by the Ethics Committee of the Sendai ART Clinic and Tohoku University School of Medicine. All samples were obtained after receiving written informed consent. All experiments handling cells and tissues were conducted in accordance with the principles outlined in the Declaration of Helsinki. Written informed consent was provided by all patients.
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.
Patients’ information. Blastocysts were graded using a blastocyst scoring system . FSH follicular stimulating hormone, HMG human menopausal gonadotropin. (XLS 51 kb)
List of TaqMan probe ID numbers. (XLS 29 kb)
Primer sequences and PCR conditions in this study. (XLS 29 kb)
Differentially expressed miRNAs with ≥2-fold changes and p < 0.05 among the three groups. (XLS 45 kb)
Validation of placental miRNAs expression using qRT-PCR. The boxplots show the expression levels of miR-197-5p (a), miR-4697-5p (b), miR-4721 (c), miR-5006-5p (d), miR-575 (e), miR-6893-5p (f), miR-125a-5p (g), miR-1260b (h), miR-224-3p (i), miR-331-3p (j), miR-365a-3p (k), miR-495-3p (l), miR-518b (m), miR-518f-3p (n), miR-543 (o) and miR-7977 (p). Data were normalized to RNU44 expression and are presented as boxplots with whiskers. The upper and lower limits of the boxes present the 75th and 25th percentiles, respectively. The upper and lower whiskers represent the maximum and minimum values that are no more than 1.5 times the span of the interquartile range (range of the values between the 25th and the 75th percentiles). The circles indicate the outliers. The median is indicated by the line in each box. They were analyzed using the Steel-Dwass test and considered statistically significant when p < 0.05. (PDF 1684 kb)
Correlation between miR-193b-3p expression and adjusted birthweight or placental weight. (a) The miR-193b-3p expression level in the placenta was correlated with adjusted birthweight. (b) The miR-193b-3p expression level in the placenta was correlated with placental weight. Black circles, white circles, and triangles indicate SP, FET and ET samples, respectively. (PDF 893 kb)
Pearson correlation coefficients between miRNA expression and adjusted birthweight, placental weight or maternal age. (XLS 33 kb)
GO term enrichment analysis of target genes. (XLS 272 kb)
KEGG enrichment analysis of target genes. (XLS 45 kb)
About this article
Cite this article
Hiura, H., Hattori, H., Kobayashi, N. et al. Genome-wide microRNA expression profiling in placentae from frozen-thawed blastocyst transfer. Clin Epigenet 9, 79 (2017). https://doi.org/10.1186/s13148-017-0379-6
- Assisted reproductive technologies (ART)
- Frozen-thawed embryo transfer (FET)
- MicroRNAs (miRNAs)
- Real-time PCR