Next Article in Journal
Systems Drug Discovery for Diffuse Large B Cell Lymphoma Based on Pathogenic Molecular Mechanism via Big Data Mining and Deep Learning Method
Next Article in Special Issue
Changes in the Metabolic Profile of Melatonin Synthesis-Related Indoles during Post-Embryonic Development of the Turkey Pineal Organ
Previous Article in Journal
To Die or Not to Die: Cell Death in Biology and Disease
Previous Article in Special Issue
Use of Melatonin in Cancer Treatment: Where Are We?
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 23
Issue 12
10.3390/ijms23126731
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Effects of Melatonin on the Transcriptome of Human Granulosa Cells, Fertilization and Blastocyst Formation

by
Isao Tamura
,
Hiroshi Tamura
*,
Mai Kawamoto-Jozaki
,
Yuichiro Shirafuta
,
Taishi Fujimura
,
Yumiko Doi-Tanaka
,
Yumiko Mihara
,
Toshiaki Taketani
and
Norihiro Sugino
Department of Obstetrics and Gynecology, Graduate School of Medicine, Yamaguchi University, Minamikogushi 1-1-1, Ube 755-8505, Japan
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2022, 23(12), 6731; https://doi.org/10.3390/ijms23126731
Submission received: 9 May 2022 / Revised: 14 June 2022 / Accepted: 14 June 2022 / Published: 16 June 2022
(This article belongs to the Special Issue Pleiotropic Benefits of Melatonin: From Basic Mechanisms to Disease)
Browse Figures
Versions Notes

Abstract

:
Melatonin is a promising reagent that can improve assisted reproductive technology (ART) outcomes in infertility patients. However, melatonin is not effective for all infertile patients, and it remains unclear for which patients melatonin would be effective. This study examined the effects of melatonin on ART outcomes and examined its mechanisms. Melatonin increased the fertilization rate in patients whose fertilization rates in the previous cycle were less than 50%, but not in patients whose fertilization rates were more than 50% in the previous cycle. Melatonin increased the blastocyst formation rate in patients whose embryo development rates in the previous cycle were less than 50%, but not in patients whose embryo development rates were more than 50% in the previous cycle. To clarify its mechanisms, transcriptome changes by melatonin treatment in granulosa cells (GCs) of the patients were examined by RNA-sequence. Melatonin treatment altered the transcriptomes of GCs of patients with poor ART outcomes so that they were similar to the transcriptomes of patients with good ART outcomes. The altered genes were associated with the inhibition of cell death and T-cell activity, and the activation of steroidogenesis and angiogenesis. Melatonin treatment was effective for patients with poor fertilization rates and poor embryo development rates in the previous ART cycle. Melatonin alters the GCs transcriptome and, thus, their functions, and this could improve the oocyte quality, leading to good ART outcomes.

1. Introduction

Melatonin (N-acetyl-5-methoxytryptamine) is a hormone secreted by the pineal gland in a circadian manner. Its secretion is regulated by light and dark environments, and it functions in the regulation of circadian rhythms and sleep [1]. Melatonin is a direct free radical scavenger with more powerful antioxidant potential than conventional antioxidants such as vitamins C and E, mannitol, and glutathione [2]. In recent years, remarkable advances have been made in assisted reproductive technology (ART) for infertility treatment [3,4]. However, satisfactory conception rates have not been achieved. The major cause of this is attributed to problems in the quality of oocytes [5]. It is thought that oxidative stress caused by reactive oxygen species (ROS) in the follicle is one of the critical factors determining oocyte quality, because oocytes under high oxidative stress show poor oocyte quality with poor fertilization ability and poor embryo development [6,7]. We previously demonstrated that melatonin protects oocytes and granulosa cells (GCs) from apoptosis by scavenging ROS produced within follicles during ovulation [6,8,9,10,11]. Our previous clinical study showed that oral melatonin administration for infertile women increases intrafollicular melatonin concentrations, reduces intrafollicular oxidative stress, and elevates fertilization and pregnancy rates by ART [11]. After our initial study, other reports also showed that melatonin administration improves the clinical outcome of ART as shown by increases in the number of mature oocytes, the fertilization rate, and the number of high-quality embryos [12,13,14,15,16]. Therefore, melatonin is a promising reagent that can improve ART outcomes by reducing intrafollicular oxidative stress. As a result, we have used melatonin to treat a number of patients who failed to become pregnant by ART. However, while melatonin treatment improved the fertilization or embryo development in some cases, it had no effect in other cases. In other words, melatonin is not effective for all infertile patients, and it remains unclear for which patients melatonin would be effective. Therefore, it is necessary to identify the characteristics of patients whose ART results may be improved by melatonin treatment.
The mechanisms by which melatonin improves oocyte quality have not been well clarified. The quality and development potential of oocytes are linked to the status of GCs [17,18]. Because melatonin protects GCs from apoptosis by scavenging ROS, melatonin may improve oocyte quality by altering the cellular functions of GCs. To investigate the altered cellular functions, it is important to examine the effects of melatonin administration on gene expressions in GCs. Although previous reports showed that melatonin alters the gene expression of human GCs, these studies were carried out only in vitro culture of GCs [19,20,21]. Therefore, it remains unclear whether oral melatonin administration alters gene expression of human GCs in vivo. This led us to examine the effects of melatonin on the transcriptome of GCs, which we expected would clarify how melatonin treatment improves ART outcomes.
In this study, we retrospectively examined patients whose ART outcomes were improved by melatonin treatment. Furthermore, to clarify its mechanisms, we obtained GCs in the ART cycle before and during melatonin treatment from the same patients and analyzed the transcriptome changes by RNA-sequence.

2. Results

2.1. Association between Previous Fertilization Rate and the Effect of Melatonin

To retrospectively evaluate the effects of melatonin on fertilization rates, the fertilization rates of 75 women were compared between the cycle before melatonin treatment and the cycle under melatonin treatment. Overall, fertilization rates were not significantly different before and under melatonin treatment (Figure 1, left graph). We divided the patients into two groups according to their fertilization rate before melatonin treatment; 19 patients whose fertilization rate before melatonin treatment was less than 50% (Poor fertilization rate cases), and 56 patients whose fertilization rate before melatonin treatment was more than 50% (Good fertilization rate cases). The backgrounds of the two groups (Table 1) were not significantly different.
Then, the effect of melatonin on the fertilization rate was compared, respectively. Melatonin significantly improved the fertilization rate in poor fertilization rate cases (Figure 1, middle graph), but no further improvement was observed in good fertilization rate cases (Figure 1, right graph).

2.2. Association between Previous Embryo Development Rate (EDR) and the Effect of Melatonin

We next evaluated the effect of melatonin on embryo development in the 56 cases with good fertilization rate, because they had a sufficient number of fertilized oocytes to analyze the subsequent embryo development. Overall, the formation rate of good quality embryos was not significantly different before and under melatonin treatment (Figure 2A, left graphs). We divided the patients into two groups according to their EDR before melatonin treatment; 41 cases with poor EDR (either the formation rate of good quality cleavage-stage embryo or the formation rate of blastocysts was less than 50% before melatonin treatment), and 15 cases with good EDR. The backgrounds of the two groups (Table 1) were not significantly different. The formation rate of good quality embryos was not significantly different before and under melatonin treatment in both good and poor EDR cases (Figure 2A, middle and right graphs).
Because the developmental capacity of embryos is more associated with the formation of blastocysts than with the quality of cleavage-stage embryos [22,23,24], we evaluated the effect of melatonin on the blastocyst formation rate. Because 18 of the 56 cases underwent embryo transfer at the cleavage-stage and had no embryos to be cultured until the blastocyst stage, they were excluded from the analysis. In total, the blastocyst formation rate was evaluated in 38 cases. Overall, the blastocyst formation rate was not significantly different before and under melatonin treatment (Figure 2B, left graph). Melatonin significantly improved the blastocyst formation rate in poor EDR cases (Figure 2B, middle graph), but no further improvement was observed in good EDR cases (Figure 2B, right graph).
We evaluated the implantation rate and pregnancy rate after melatonin treatment. In poor fertilization rate cases, 12 cases (22 embryos) underwent embryo transfer in the cycle during the melatonin treatment and 1 case became pregnant (implantation rate; 4.5%, pregnancy rate; 8.3%). In poor EDR cases, 25 cases (38 embryos) underwent embryo transfer in the cycle during the melatonin treatment and 4 cases became pregnant (implantation rate; 10.5%, pregnancy rate; 16.0%). In good EDR cases, 11 cases (16 embryos) underwent embryo transfer in the cycle during the melatonin treatment and 1 case became pregnant (implantation rate; 6.3%, pregnancy rate; 9.1%). There were no significant differences in the implantation rate and pregnancy rate between three groups. Therefore, it was indicated that implantation rate and pregnancy rate were not associated with the improvement in fertilization or blastocyst formation rate in this study.

2.3. Effect of Melatonin Treatment on Transcriptome Profile in GCs

For this experiment, we obtained GCs in the cycles before and during melatonin treatment from three patients whose blastocyst formation was significantly improved by melatonin treatment (Table 2, Cases No. 1–3). For controls, we obtained GCs from four patients whose blastocyst formation rate was more than 50% and who became pregnant without melatonin treatment (Table 2, Cases No. 4–7). These 10 samples were subjected to RNA-sequence analysis. As shown in Figure 3A, a hierarchical clustering analysis showed that GCs before melatonin treatment (Cases No. 1–3; indicated by blue boxes) were classified into different clusters of those with good outcomes without melatonin treatment (Cases No. 4–7; indicated by yellow boxes). GCs under the melatonin treatment (Cases No. 1–3; indicated by green boxes) were classified into the same clusters as GCs with good outcomes without melatonin treatment. These results suggest that melatonin treatment altered the transcriptomes of GCs of patients with poor ART outcomes so they were similar to the transcriptomes of patients with good ART outcomes.
Melatonin treatment down-regulated 97 genes and up-regulated 152 genes (Table S1). The top 20 down- or up-regulated genes according to the p-values are defined as the top differentially expressed genes (DEGs) (Table 3). The expression levels of several top DEGs (COX7B, CXCL, TRAF6 and CENPV) were examined by real time RT-PCR, which validated that melatonin treatment significantly altered their expression levels (Figure S1). According to a gene ontology analysis, the down-regulated genes were associated with cell death (circled in red) and T-cell activation (circled in blue) (Figure 3B, Table S2), while the up-regulated genes were associated with steroidogenesis (circled in green) and angiogenesis (circled in yellow) (Figure 3C, Table S2). The changes in mRNA expression level of the representative genes in each function are shown in Figure 3D,E. Although we performed a pathway analysis, there were no enriched terms in up- or down-regulated genes.

3. Discussion

This study revealed that melatonin treatment is effective for patients who had a poor fertilization rate or poor embryo development rate (EDR) in the previous ART cycle. Furthermore, this is the first study to show that melatonin treatment alters the transcriptome in GCs, which may contribute to improving ART outcomes.
We and others previously reported the favorable effects of melatonin administration on ART outcomes [11,12,13,14,15,16]. However, melatonin is not always effective for all patients, and it remains unclear for which patients melatonin would be effective. The present results reveal that melatonin increases the fertilization rate in patients whose fertilization rates in the previous cycle were less than 50%, but not in patients whose fertilization rates were more than 50% in the previous cycle. In addition, melatonin increases the blastocyst formation rate in patients whose EDR in the previous cycle was less than 50%, but not in patients whose EDR was more than 50% in the previous cycle. ART outcomes are closely associated with oocyte quality. One of critical factors determining oocyte quality is oxidative stress caused by ROS in the follicle, because excessive ROS is cytotoxic to oocytes and GCs [6,25]. In fact, the oocytes retrieved from the follicles under high oxidative stress show poor quality, and poor ability of fertilization and embryo development [7,11]. We previously reported that melatonin administration is effective in decreasing oxidative stress and improving oocyte quality [11]. Therefore, it is likely that the cases with poor fertilization rate or poor EDR were under oxidative stress and that melatonin treatment could improve ART outcomes by scavenging ROS. These facts appear to explain why melatonin treatment was effective only for patients with poor fertilization rate and poor EDR. In fact, we found that melatonin has protective effects on mouse GCs in the presence of oxidative stress, but not in its absence [26]. Taken together, our results provide important insights on the selection of patients for melatonin treatment. Although melatonin treatment improved the blastocyst formation rate, the formation rate of good quality cleavage stage embryos was not affected. This is not surprising because the development capacity of embryos is more closely associated with blastocyst formation than with the quality of cleavage-stage embryos [22,23]. Only 30% of cleavage-stage embryos can reach the blastocyst stage [24]. These facts indicate that the developmental capacity of embryos should be determined by monitoring the developmental process up to the blastocyst stage. Therefore, it is reasonable to conclude that melatonin improves embryo development because melatonin increased the blastocyst formation rate.
It is well known that GCs are involved in the regulation of oocyte quality [17,18]. Because melatonin protects GCs from apoptosis by scavenging ROS [6,8,9,10,11], it may improve oocyte quality by altering the cellular functions of GCs. Therefore, it is important to examine whether melatonin treatment alters the gene expressions of GCs. Several studies have shown that supplementing culture medium with melatonin alters the gene expressions of VEGF, StAR and CYP19A1 in human GCs, which were collected at oocyte retrieval [20,21,27]. However, these GCs no longer exist in the follicles and are not associated with oocytes. Therefore, it is unclear whether these in vitro effects of melatonin are also true in GCs within the follicles. While melatonin administration in vivo has previously been shown to affect gene expressions in mouse GCs, such as Sirt1 and Myo10 [28,29], the present study is the first to show it also affects the transcriptome in human GCs. Our analysis revealed that melatonin administration altered the transcriptome in human GCs of patients with poor ART outcomes into that of patients with good ART outcomes, indicating that melatonin improves ART outcomes by altering the expressions of a number of genes in human GCs. In other words, melatonin improves the cellular functions of GCs by reducing oxidative stress.
Gene ontology analysis revealed that one of the cellular functions associated with down-regulated genes by melatonin treatment was cell death, including the genes such as RASSF2 and TNFRSF6B (Figure 3D). These genes are associated with cell death because their expressions increase under oxidative stress, which induces apoptosis [30,31,32]. These results are consistent with the previously reported function of melatonin in protecting GCs from apoptosis by scavenging ROS [6,8,9,10,11]. Another cellular function associated with down-regulated genes was T-cell activation, including the gene CCL5 (Figure 3D). The effect of melatonin on immunomodulation has not been well clarified. CCL5 is induced by oxidative stress and is involved in the chemotaxis and activation of T lymphocytes [33]. This suggests that by decreasing CCL5 expression in GCs, melatonin suppresses the accumulation of activated T lymphocytes around GCs. In fact, melatonin suppresses both oxidative stress-induced CCL5 expression in bovine mammary epithelial cells [34] and the chemotaxis of leukocytes [35]. Given that apoptosis is induced by the accumulation of cytotoxic T-cells [36], melatonin may suppress the cellular function of GCs that activates T-cells, which contributes to the suppression of GC apoptosis.
Among the genes up-regulated by melatonin, some seem to be associated with angiogenesis. Angiogenic factors secreted by GCs play essential roles in angiogenesis during ovulation [37,38]. We reported that the expression of vascular endothelial growth factor (Vegf) increases after the LH surge in rat GCs [39]. The vascularization status of ovarian follicles affects the quality of oocytes by regulating the oxygen supply. An adequate blood supply increases the intrafollicular oxygen levels [40]. A decrease in the intrafollicular oxygen level due to poor vascularization at the perifollicular region causes metabolic disorders and chromosomal abnormalities in oocytes, leading to poor oocyte quality [40,41]. Therefore, angiogenesis is an essential event determining oocyte quality. In fact, the concentrations of angiogenic factors, including VEGF, are higher in follicular fluid with mature oocytes than in that with immature oocytes [42,43]. In addition, oocytes from the follicle with higher perifollicular vascularity show higher rates of fertilization, embryo development and implantation than those with poor vascularity [44,45]. Therefore, it is likely that melatonin improves the angiogenetic function of GCs, which contributes to improving oocyte quality. In fact, melatonin increases VEGF expression in cultured human GCs, supporting the involvement of melatonin in angiogenesis [20].
Another cellular function associated with up-regulated genes is steroidogenesis. The relationship between melatonin and steroidogenesis has been reported in cultured human GCs [19,21,27]. We previously reported that the cytoprotective effects of melatonin as an antioxidant are essential for progesterone production in human GCs [46]. The present study found that melatonin up-regulated the genes associated with steroidogenesis, such as ACOX2, CYP2R1 and HNF1A (Figure 3D). These genes are involved in the metabolism of cholesterol [47,48,49], which is a precursor of progesterone. Progesterone production from GCs is suggested to be important for oocyte maturation during ovulation [50]. Therefore, it is likely that melatonin treatment improves progesterone production in GCs during ovulation, which may improve oocyte quality.
In addition to the role of melatonin as a direct radical scavenger, melatonin can affect gene expressions via the membrane melatonin receptors (MT1 and MT2), both of which are expressed in GCs and oocytes [51,52]. Therefore, there is a possibility that melatonin alters the cellular functions of GCs by activating the signaling pathway downstream of melatonin receptors. However, considering that melatonin treatment is effective for patients with poor fertilization or poor EDR, and that their oocytes are under higher oxidative stress [7,11], it is more likely that melatonin improves the cellular functions of GCs by scavenging ROS.
In conclusion, the present study showed that melatonin treatment is effective for patients with poor fertilization rate and poor ER in the previous ART cycle. This will help to identify patients that would most benefit from melatonin treatment. Melatonin reduces oxidative stress and protects GCs in patients with poor ART outcomes. This improves the cellular functions of GCs by altering the transcriptomes of GCs, which may contribute to improving oocyte quality, leading to good ART outcomes.

4. Materials and Methods

4.1. Study Population

This retrospective study included 75 female patients who had previously failed to become pregnant by conventional in vitro fertilization (cIVF) or intracytoplasmic sperm injection (ICSI), followed by embryo transfer at Yamaguchi University Hospital. They were treated with melatonin in the next ART cycle. Informed consent was obtained from all the patients in this study. The study design was reviewed and approved by the institutional review board of Yamaguchi University Hospital.

4.2. ART Procedure

Controlled ovarian hyperstimulation (COH) was performed using standard gonadotropin releasing hormone agonist (GnRHa)/follicle-stimulating hormone (FSH) protocols, as we reported previously [53]. Nasal spray GnRHa (900 μg/day) was given from the mid-luteal phase in the previous cycle to continuously suppress pituitary gonadotropin secretion until the injection of HCG (10,000 IU). COH was initiated from the 2nd day of the IVF-ET cycle by the injection of 225 IU FSH for 3 days, followed by a daily injection of 150 IU HMG. When more than three leading follicles reached 18 mm or more, HCG was injected for ovulation induction. Oocyte retrieval was carried out 35 h after HCG injection. For patients who showed poor response to COH in the previous IVF cycle, they underwent short GnRH agonist protocol. In this protocol, COH and GnRHa were simultaneously initiated from the 2nd day of the IVF-ET cycle. 300 IU FSH or 300 IU HMG was daily injected until the day of HCG injection. The semen samples were collected by masturbation and the motile sperm was collected by the swim-up technique. Oocytes were fertilized by either cIVF or ICSI with Piezo-assisted ICSI system (Prime Tech Ltd., Tokyo, Japan). If the sperm concentration was low (motile sperm concentration is less than 5 × 106/mL after swim-up) or there was a history of fertilization failure (previous fertilization rate under cIVF was less than 25%), all retrieved oocytes were fertilized by ICSI. Fertilization was confirmed on day 1 (17–19 h after insemination) by the presence of two pronuclei. The quality of the cleavage-stage embryo was assessed on day 2 according to Veeck’s criteria. Grade 1 and grade 2 embryos were defined as good quality embryos, as reported previously [54,55]. Blastocyst formation was assessed on day 5. A fresh embryo was transferred on either day 2 (cleavage stage) or day 5 (blastocyst stage).

4.3. Melatonin Treatment

Patients who failed to become pregnant in the previous ART cycle due to poor fertilization rate, poor embryo development rate (EDR) (low formation rate of the good quality embryo or low formation rate of blastocyst), or repeated ART failure were administered melatonin orally in the next cycle, as reported previously [11]. They were given a 3 mg tablet of melatonin orally at 22:00 h from the fifth day of the previous menstrual cycle until the day of oocyte retrieval. For patients who were treated with a short GnRH agonist protocol, melatonin was simultaneously started on day 1 of COH.

4.4. Outcomes

The numbers of retrieved oocytes, fertilized oocytes, cleavage-stage embryos, and their quality, and blastocysts were evaluated in each cycle. Fertilization rate was expressed as the ratio of the number of fertilized oocytes to the number of retrieved oocytes. The formation rates of good quality embryos and blastocysts were expressed as a ratio of the number of good quality embryos (grade1 and grade 2) or blastocysts to the number of fertilized oocytes. If the embryo was transferred at the cleavage stage, the blastocyst formation rate was expressed as the ratio of the number of blastocysts to the number of fertilized oocytes minus the number of transferred embryos. Cases where embryos were transferred at the cleavage stage were excluded from the analysis of blastocyst formation rate. Outcomes of the cycles before and during melatonin treatment were compared. Cases in which different fertilization methods (cIVF vs. ICSI) were used between the cycles before and under melatonin treatment were excluded from the analysis. Implantation rate was expressed as the ratio of the number of gestational sacs to the number of transferred embryos. Pregnancy rate was expressed as the ratio of the number of pregnant cases to the number of cases who underwent fresh embryo transfer.

4.5. Transcriptome Analysis of Human GCs

Human GCs obtained from seven patients at oocyte retrieval were used for the transcriptome analysis. For three of the seven patients, in which blastocyst formation was not observed before melatonin treatment but was significantly improved under melatonin treatment, GCs were obtained in the cycles before and under melatonin treatment. The other four patients had good ART outcomes (their blastocyst formation rate was more than 50% and they became pregnant by the transfer of these embryos). Their GCs were used as controls of good ART outcomes. The ages and blastocyst formation rates of the seven patients are shown in Table S1. GCs were collected from follicular fluids of IVF-ET patients as reported previously [46]. After the removal of the oocyte–cumulus complex, cells were carefully aspirated from the follicular fluid and subsequently separated from red blood cells with Ficoll-Paque (Amersham Pharmacia Biotech AB, Uppsala, Sweden) by centrifuging at 500× g for 20 min. Total RNA was isolated from the cultured cells with a RNeasy® Mini Kit (Qiagen, Valencia, CA, USA). In total, 10 samples (3 before-melatonin treated GCs, 3 under-melatonin treated GCs, and 4 control GCs) were subjected to RNA-sequence analysis, as reported previously [56,57,58]. mRNA was purified with oligo dT beads (NEBNext Poly (A) mRNA magnet Isolation Module, New England Biolabs, NEB, Beverly, MA, USA). Complementary DNA (cDNA) libraries for Illumina sequencing were generated with NEBNext Ultra II RNA library Prep kit (NEB) and NEBNextplex Oligos for Illumina. The confirmed libraries were sequenced with an Illumina Next-seq DNA sequencer with a 75-bp pair-end cycle sequencing kit (Illumina, San Diego, CA, USA). To produce the raq bcl, or base call files, quality assessment and image analyses were performed using Next-seq packaging software Real Time Analysis, and bcl2fastaq Conversion Software v2.19 (Illumina) was used for de-multiplexing the samples. Reads with more than two ambiguous nucleotides and reads with quality scores less than 20 as calculated by the Phred program were removed using CLC Genomics Workbench software (ver. 8.01, Qiagen, Valencia, CA, USA). Long reads with more than 1000 nucleotides and short reads with fewer than 20 nucleotides were also discarded. Trimmed reads were mapped to the human reference genome GRCh37 in default settings. The numbers of the identified transcript and sequenced paired end reads in each sample are shown in Table S3. Clustering analyses were performed using the R package “Cluster” [59]. Gene expression values were calculated as “transcripts per million” (TPM). We added 0.5 to the TPM value before the following calculation. Genes whose TPM increased more than 1.4-fold by melatonin treatment and whose p-values were < 0.05 by paired t-test were defined as genes that were up-regulated by melatonin treatment. Similarly, genes whose TPM decreased more than 1.4-fold by melatonin treatment and whose p-values were < 0.05 were defined as genes that were down-regulated by melatonin treatment. The top 20 down- or up-regulated genes according to the p-values are defined as the top differentially expressed genes (DEGs) (Table 3). Gene ontology was analyzed with DAVID Bioinformatics Resources v.6.8 to determine the functions of these genes [60]. Gene ontology terms with p-values < 0.05 were considered as enriched terms. These terms were summarized by removing redundancy and plotted using reduce and visualize gene ontology (REVIGO) with Allowed Similarity as “medium (0.5)” [61,62].

4.6. Real-Time RT-PCR

The relative mRNA expressions of several top DEGs (COX7B, CXCL, TRAF6 and CENPV) were examined by real-time RT-PCR as reported previously [56] with sequence-specific primer sets (Table S4). GAPDH was used as an internal control.

4.7. Statistical Analysis

Mann–Whitney U test or Chi-squared test was applied to analyze the difference in backgrounds of the two groups. Wilcoxon’s rank-sum test was used to evaluate the differences in ART outcomes before and under melatonin treatment. All statistical analyses were performed using SPSS for Windows version 11 (SPSS Inc., Chicago, IL, USA). Differences were considered significant at p-values < 0.05.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/ijms23126731/s1

Author Contributions

Conceptualization, I.T. and H.T.; methodology, I.T. and Y.S.; formal analysis, I.T., M.K.-J., T.F., Y.D.-T., Y.M. and T.T.; writing—original draft preparation, I.T.; writing—review and editing, N.S. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported in part by JSPS KAKENHI (Grant Number 20K09645, 20K18191, 20H03825, 21K16816, 21K09542, 21K09517).

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki, and approved by the Institutional Review Board of Yamaguchi University Hospital. (protocol code; H27-190, date of approval; 24 February 2016).

Informed Consent Statement

Not applicable.

Data Availability Statement

RNA-sequence data were deposited in the Gene Expression Omnibus (number GSE198742).

Acknowledgments

The authors thank the Open Research Facilities of the Institute of Gene Research of Yamaguchi University Science Research Center for RNA-sequence analysis.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhao, D.; Yu, Y.; Shen, Y.; Liu, Q.; Zhao, Z.; Sharma, R.; Reiter, R.J. Melatonin Synthesis and Function: Evolutionary History in Animals and Plants. Front. Endocrinol. 2019, 10, 249. [Google Scholar] [CrossRef] [PubMed]
  2. Tamura, H.; Jozaki, M.; Tanabe, M.; Shirafuta, Y.; Mihara, Y.; Shinagawa, M.; Tamura, I.; Maekawa, R.; Sato, S.; Taketani, T.; et al. Importance of Melatonin in Assisted Reproductive Technology and Ovarian Aging. Int. J. Mol. Sci. 2020, 21, 1135. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Jiang, L.; Qian, Y.; Chen, X.; Ji, X.; Ou, S.; Li, R.; Yang, D.; Li, Y. Effect of early rescue ICSI and split IVF-ICSI in preventing low fertilization rate during the first ART cycle: A real-world retrospective cohort study. Reprod. Med. Biol. 2022, 21, e12420. [Google Scholar] [CrossRef] [PubMed]
  4. Ishihara, O.; Jwa, S.C.; Kuwahara, A.; Katagiri, Y.; Kuwabara, Y.; Hamatani, T.; Harada, M.; Ichikawa, T. Assisted reproductive technology in Japan: A summary report for 2017 by the Ethics Committee of the Japan Society of Obstetrics and Gynecology. Reprod. Med. Biol. 2020, 19, 3–12. [Google Scholar] [CrossRef] [Green Version]
  5. Lee, J. Is age-related increase of chromosome segregation errors in mammalian oocytes caused by cohesin deterioration? Reprod. Med. Biol. 2020, 19, 32–41. [Google Scholar] [CrossRef] [Green Version]
  6. Tamura, H.; Takasaki, A.; Taketani, T.; Tanabe, M.; Kizuka, F.; Lee, L.; Tamura, I.; Maekawa, R.; Asada, H.; Yamagata, Y.; et al. Melatonin as a free radical scavenger in the ovarian follicle. Endocr. J. 2013, 60, 1–13. [Google Scholar] [CrossRef] [Green Version]
  7. Seino, T.; Saito, H.; Kaneko, T.; Takahashi, T.; Kawachiya, S.; Kurachi, H. Eight-hydroxy-2′-deoxyguanosine in granulosa cells is correlated with the quality of oocytes and embryos in an in vitro fertilization-embryo transfer program. Fertil. Steril. 2002, 77, 1184–1190. [Google Scholar] [CrossRef]
  8. Tamura, H.; Takasaki, A.; Taketani, T.; Tanabe, M.; Kizuka, F.; Lee, L.; Tamura, I.; Maekawa, R.; Aasada, H.; Yamagata, Y.; et al. The role of melatonin as an antioxidant in the follicle. J. Ovarian Res. 2012, 5, 5. [Google Scholar] [CrossRef] [Green Version]
  9. Tamura, H.; Takasaki, A.; Taketani, T.; Tanabe, M.; Lee, L.; Tamura, I.; Maekawa, R.; Aasada, H.; Yamagata, Y.; Sugino, N. Melatonin and female reproduction. J. Obstet. Gynaecol. Res. 2014, 40, 1–11. [Google Scholar] [CrossRef]
  10. Tamura, H.; Nakamura, Y.; Korkmaz, A.; Manchester, L.C.; Tan, D.X.; Sugino, N.; Reiter, R.J. Melatonin and the ovary: Physiological and pathophysiological implications. Fertil. Steril. 2009, 92, 328–343. [Google Scholar] [CrossRef]
  11. Tamura, H.; Takasaki, A.; Miwa, I.; Taniguchi, K.; Maekawa, R.; Asada, H.; Taketani, T.; Matsuoka, A.; Yamagata, Y.; Shimamura, K.; et al. Oxidative stress impairs oocyte quality and melatonin protects oocytes from free radical damage and improves fertilization rate. J. Pineal Res. 2008, 44, 280–287. [Google Scholar] [CrossRef] [PubMed]
  12. Batioglu, A.S.; Sahin, U.; Gurlek, B.; Ozturk, N.; Unsal, E. The efficacy of melatonin administration on oocyte quality. Gynecol. Endocrinol. 2012, 28, 91–93. [Google Scholar] [CrossRef] [PubMed]
  13. Eryilmaz, O.G.; Devran, A.; Sarikaya, E.; Aksakal, F.N.; Mollamahmutoglu, L.; Cicek, N. Melatonin improves the oocyte and the embryo in IVF patients with sleep disturbances, but does not improve the sleeping problems. J. Assist. Reprod. Genet. 2011, 28, 815–820. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Jahromi, B.N.; Sadeghi, S.; Alipour, S.; Parsanezhad, M.E.; Alamdarloo, S.M. Effect of Melatonin on the Outcome of Assisted Reproductive Technique Cycles in Women with Diminished Ovarian Reserve: A Double-Blinded Randomized Clinical Trial. Iran J. Med. Sci. 2017, 42, 73–78. [Google Scholar]
  15. Nishihara, T.; Hashimoto, S.; Ito, K.; Nakaoka, Y.; Matsumoto, K.; Hosoi, Y.; Morimoto, Y. Oral melatonin supplementation improves oocyte and embryo quality in women undergoing in vitro fertilization-embryo transfer. Gynecol. Endocrinol. 2014, 30, 359–362. [Google Scholar] [CrossRef]
  16. Espino, J.; Macedo, M.; Lozano, G.; Ortiz, A.; Rodriguez, C.; Rodriguez, A.B.; Bejarano, I. Impact of Melatonin Supplementation in Women with Unexplained Infertility Undergoing Fertility Treatment. Antioxidants 2019, 8, 338. [Google Scholar] [CrossRef] [Green Version]
  17. Alam, M.H.; Miyano, T. Interaction between growing oocytes and granulosa cells in vitro. Reprod. Med. Biol. 2020, 19, 13–23. [Google Scholar] [CrossRef] [Green Version]
  18. Regan, S.L.P.; Knight, P.G.; Yovich, J.L.; Leung, Y.; Arfuso, F.; Dharmarajan, A. Granulosa Cell Apoptosis in the Ovarian Follicle—A Changing View. Front. Endocrinol. 2018, 9, 61. [Google Scholar] [CrossRef]
  19. Fang, L.; Li, Y.; Wang, S.; Yu, Y.; Li, Y.; Guo, Y.; Yan, Y.; Sun, Y.P. Melatonin induces progesterone production in human granulosa-lutein cells through upregulation of StAR expression. Aging 2019, 11, 9013–9024. [Google Scholar] [CrossRef]
  20. Li, Y.; Fang, L.; Zhang, R.; Wang, S.; Li, Y.; Yan, Y.; Yu, Y.; Cheng, J.C.; Sun, Y.P. Melatonin stimulates VEGF expression in human granulosa-lutein cells: A potential mechanism for the pathogenesis of ovarian hyperstimulation syndrome. Mol. Cell Endocrinol. 2020, 518, 110981. [Google Scholar] [CrossRef]
  21. Cheng, J.C.; Fang, L.; Li, Y.; Wang, S.; Li, Y.; Yan, Y.; Jia, Q.; Wu, Z.; Wang, Z.; Han, X.; et al. Melatonin stimulates aromatase expression and estradiol production in human granulosa-lutein cells: Relevance for high serum estradiol levels in patients with ovarian hyperstimulation syndrome. Exp. Mol. Med. 2020, 52, 1341–1350. [Google Scholar] [CrossRef] [PubMed]
  22. Gardner, D.K.; Schoolcraft, W.B.; Wagley, L.; Schlenker, T.; Stevens, J.; Hesla, J. A prospective randomized trial of blastocyst culture and transfer in in-vitro fertilization. Hum. Reprod. 1998, 13, 3434–3440. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Balaban, B.; Urman, B.; Sertac, A.; Alatas, C.; Aksoy, S.; Mercan, R. Blastocyst quality affects the success of blastocyst-stage embryo transfer. Fertil. Steril. 2000, 74, 282–287. [Google Scholar] [CrossRef]
  24. Dumoulin, J.C.; Coonen, E.; Bras, M.; van Wissen, L.C.; Ignoul-Vanvuchelen, R.; Bergers-Jansen, J.M.; Derhaag, J.G.; Geraedts, J.P.; Evers, J.L. Comparison of in-vitro development of embryos originating from either conventional in-vitro fertilization or intracytoplasmic sperm injection. Hum. Reprod. 2000, 15, 402–409. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. He, C.; Wang, J.; Zhang, Z.; Yang, M.; Li, Y.; Tian, X.; Ma, T.; Tao, J.; Zhu, K.; Song, Y.; et al. Mitochondria Synthesize Melatonin to Ameliorate Its Function and Improve Mice Oocyte’s Quality under in Vitro Conditions. Int. J. Mol. Sci. 2016, 17, 939. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Tanabe, M.; Tamura, H.; Taketani, T.; Okada, M.; Lee, L.; Tamura, I.; Maekawa, R.; Asada, H.; Yamagata, Y.; Sugino, N. Melatonin protects the integrity of granulosa cells by reducing oxidative stress in nuclei, mitochondria, and plasma membranes in mice. J. Reprod. Dev. 2015, 61, 35–41. [Google Scholar] [CrossRef] [Green Version]
  27. Yu, K.; Wang, R.X.; Li, M.H.; Sun, T.C.; Zhou, Y.W.; Li, Y.Y.; Sun, L.H.; Zhang, B.L.; Lian, Z.X.; Xue, S.G.; et al. Melatonin Reduces Androgen Production and Upregulates Heme Oxygenase-1 Expression in Granulosa Cells from PCOS Patients with Hypoestrogenia and Hyperandrogenia. Oxid. Med. Cell Longev. 2019, 2019, 8218650. [Google Scholar] [CrossRef]
  28. Shen, M.; Cao, Y.; Jiang, Y.; Wei, Y.; Liu, H. Melatonin protects mouse granulosa cells against oxidative damage by inhibiting FOXO1-mediated autophagy: Implication of an antioxidation-independent mechanism. Redox. Biol. 2018, 18, 138–157. [Google Scholar] [CrossRef]
  29. Zhang, H.; Li, C.; Wen, D.; Li, R.; Lu, S.; Xu, R.; Tang, Y.; Sun, Y.; Zhao, X.; Pan, M. Melatonin improves the quality of maternally aged oocytes by maintaining intercellular communication and antioxidant metabolite supply. Redox. Biol. 2021, 49, 102215. [Google Scholar] [CrossRef]
  30. Vos, M.D.; Ellis, C.A.; Elam, C.; Ulku, A.S.; Taylor, B.J.; Clark, G.J. RASSF2 is a novel K-Ras-specific effector and potential tumor suppressor. J. Biol. Chem. 2003, 278, 28045–28051. [Google Scholar] [CrossRef] [Green Version]
  31. Hudlikar, R.R.; Sargsyan, D.; Wu, R.; Su, S.; Zheng, M.; Kong, A.N. Triterpenoid corosolic acid modulates global CpG methylation and transcriptome of tumor promotor TPA induced mouse epidermal JB6 P+ cells. Chem. Biol. Interact. 2020, 321, 109025. [Google Scholar] [CrossRef] [PubMed]
  32. Hsieh, S.L.; Lin, W.W. Decoy receptor 3: An endogenous immunomodulator in cancer growth and inflammatory reactions. J. Biomed. Sci. 2017, 24, 39. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Chan, O.; Burke, J.D.; Gao, D.F.; Fish, E.N. The chemokine CCL5 regulates glucose uptake and AMP kinase signaling in activated T cells to facilitate chemotaxis. J. Biol. Chem. 2012, 287, 29406–29416. [Google Scholar] [CrossRef] [Green Version]
  34. Yu, G.M.; Kubota, H.; Okita, M.; Maeda, T. The anti-inflammatory and antioxidant effects of melatonin on LPS-stimulated bovine mammary epithelial cells. PLoS ONE 2017, 12, e0178525. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Kepka, M.; Szwejser, E.; Pijanowski, L.; Verburg-van Kemenade, B.M.; Chadzinska, M. A role for melatonin in maintaining the pro- and anti-inflammatory balance by influencing leukocyte migration and apoptosis in carp. Dev. Comp. Immunol. 2015, 53, 179–190. [Google Scholar] [CrossRef]
  36. Zhang, N.; Hartig, H.; Dzhagalov, I.; Draper, D.; He, Y.W. The role of apoptosis in the development and function of T lymphocytes. Cell Res. 2005, 15, 749–769. [Google Scholar] [CrossRef]
  37. Bender, H.R.; Trau, H.A.; Duffy, D.M. Placental Growth Factor Is Required for Ovulation, Luteinization, and Angiogenesis in Primate Ovulatory Follicles. Endocrinology 2018, 159, 710–722. [Google Scholar] [CrossRef]
  38. Kizuka, F.; Tokuda, N.; Takagi, K.; Adachi, Y.; Lee, L.; Tamura, I.; Maekawa, R.; Taketani, T.; Tamura, H.; Suzuki, T.; et al. Involvement of bone marrow-derived vascular progenitor cells in neovascularization during formation of the corpus luteum in mice. Biol. Reprod. 2012, 87, 55. [Google Scholar] [CrossRef]
  39. Shinagawa, M.; Tamura, I.; Maekawa, R.; Sato, S.; Shirafuta, Y.; Mihara, Y.; Okada-Matsumoto, M.; Taketani, T.; Asada, H.; Tamura, H.; et al. C/EBPbeta regulates Vegf gene expression in granulosa cells undergoing luteinization during ovulation in female rats. Sci. Rep. 2019, 9, 714. [Google Scholar] [CrossRef] [Green Version]
  40. Van Blerkom, J.; Antczak, M.; Schrader, R. The developmental potential of the human oocyte is related to the dissolved oxygen content of follicular fluid: Association with vascular endothelial growth factor levels and perifollicular blood flow characteristics. Hum. Reprod. 1997, 12, 1047–1055. [Google Scholar] [CrossRef]
  41. Van Blerkom, J. Intrafollicular influences on human oocyte developmental competence: Perifollicular vascularity, oocyte metabolism and mitochondrial function. Hum. Reprod. 2000, 15 (Suppl. 2), 173–188. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Bokal, E.V.; Vrtovec, H.M.; Virant Klun, I.; Verdenik, I. Prolonged HCG action affects angiogenic substances and improves follicular maturation, oocyte quality and fertilization competence in patients with polycystic ovarian syndrome. Hum. Reprod. 2005, 20, 1562–1568. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Malamitsi-Puchner, A.; Sarandakou, A.; Baka, S.G.; Tziotis, J.; Rizos, D.; Hassiakos, D.; Creatsas, G. Concentrations of angiogenic factors in follicular fluid and oocyte-cumulus complex culture medium from women undergoing in vitro fertilization: Association with oocyte maturity and fertilization. Fertil. Steril. 2001, 76, 98–101. [Google Scholar] [CrossRef]
  44. Naredi, N.; Singh, S.K.; Sharma, R. Does Perifollicular Vascularity on the Day of Oocyte Retrieval Affect Pregnancy Outcome in an In Vitro Fertilization Cycle? J. Hum. Reprod. Sci. 2017, 10, 281–287. [Google Scholar] [CrossRef] [PubMed]
  45. Bhal, P.S.; Pugh, N.D.; Chui, D.K.; Gregory, L.; Walker, S.M.; Shaw, R.W. The use of transvaginal power Doppler ultrasonography to evaluate the relationship between perifollicular vascularity and outcome in in-vitro fertilization treatment cycles. Hum. Reprod. 1999, 14, 939–945. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Taketani, T.; Tamura, H.; Takasaki, A.; Lee, L.; Kizuka, F.; Tamura, I.; Taniguchi, K.; Maekawa, R.; Asada, H.; Shimamura, K.; et al. Protective role of melatonin in progesterone production by human luteal cells. J. Pineal Res. 2011, 51, 207–213. [Google Scholar] [CrossRef]
  47. Miller, W.L. Genetic disorders of Vitamin D biosynthesis and degradation. J. Steroid. Biochem. Mol. Biol. 2017, 165, 101–108. [Google Scholar] [CrossRef]
  48. Nenicu, A.; Luers, G.H.; Kovacs, W.; David, M.; Zimmer, A.; Bergmann, M.; Baumgart-Vogt, E. Peroxisomes in human and mouse testis: Differential expression of peroxisomal proteins in germ cells and distinct somatic cell types of the testis. Biol. Reprod. 2007, 77, 1060–1072. [Google Scholar] [CrossRef] [Green Version]
  49. Shih, D.Q.; Bussen, M.; Sehayek, E.; Ananthanarayanan, M.; Shneider, B.L.; Suchy, F.J.; Shefer, S.; Bollileni, J.S.; Gonzalez, F.J.; Breslow, J.L.; et al. Hepatocyte nuclear factor-1alpha is an essential regulator of bile acid and plasma cholesterol metabolism. Nat. Genet. 2001, 27, 375–382. [Google Scholar] [CrossRef]
  50. Salehnia, M.; Zavareh, S. The effects of progesterone on oocyte maturation and embryo development. Int. J. Fertil. Steril. 2013, 7, 74–81. [Google Scholar]
  51. Pang, Y.; Zhao, S.; Sun, Y.; Jiang, X.; Hao, H.; Du, W.; Zhu, H. Protective effects of melatonin on the in vitro developmental competence of bovine oocytes. Anim. Sci. J. 2018, 89, 648–660. [Google Scholar] [CrossRef] [PubMed]
  52. Fang, Y.; Zhang, J.; Li, Y.; Guo, X.; Li, J.; Zhong, R.; Zhang, X. Melatonin-induced demethylation of antioxidant genes increases antioxidant capacity through RORalpha in cumulus cells of prepubertal lambs. Free Radic. Biol. Med. 2019, 131, 173–183. [Google Scholar] [CrossRef] [PubMed]
  53. Tamura, I.; Kawamoto-Jozaki, M.; Fujimura, T.; Doi-Tanaka, Y.; Takagi, H.; Shirafuta, Y.; Mihara, Y.; Taketani, T.; Tamura, H.; Sugino, N. Relationship between follicular size and developmental capacity of oocytes under controlled ovarian hyperstimulation in assisted reproductive technologies. Reprod. Med. Biol. 2021, 20, 299–304. [Google Scholar] [CrossRef] [PubMed]
  54. Nomura, M.; Iwase, A.; Furui, K.; Kitagawa, T.; Matsui, Y.; Yoshikawa, M.; Kikkawa, F. Preferable correlation to blastocyst development and pregnancy rates with a new embryo grading system specific for day 3 embryos. J. Assist. Reprod. Genet. 2007, 24, 23–28. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Tomsu, M.; Sharma, V.; Miller, D. Embryo quality and IVF treatment outcomes may correlate with different sperm comet assay parameters. Hum. Reprod. 2002, 17, 1856–1862. [Google Scholar] [CrossRef] [Green Version]
  56. Tamura, I.; Maekawa, R.; Jozaki, K.; Ohkawa, Y.; Takagi, H.; Doi-Tanaka, Y.; Shirafuta, Y.; Mihara, Y.; Taketani, T.; Sato, S.; et al. Transcription factor C/EBPbeta induces genome-wide H3K27ac and upregulates gene expression during decidualization of human endometrial stromal cells. Mol. Cell Endocrinol. 2021, 520, 111085. [Google Scholar] [CrossRef]
  57. Tamura, I.; Ohkawa, Y.; Sato, T.; Suyama, M.; Jozaki, K.; Okada, M.; Lee, L.; Maekawa, R.; Asada, H.; Sato, S.; et al. Genome-wide analysis of histone modifications in human endometrial stromal cells. Mol. Endocrinol. 2014, 28, 1656–1669. [Google Scholar] [CrossRef] [Green Version]
  58. Watanabe, K.; Yamamoto, S.; Sakaguti, S.; Isayama, K.; Oka, M.; Nagano, H.; Mizukami, Y. A novel somatic mutation of SIN3A detected in breast cancer by whole-exome sequencing enhances cell proliferation through ERalpha expression. Sci. Rep. 2018, 8, 16000. [Google Scholar] [CrossRef] [Green Version]
  59. Yu, G.; Wang, L.G.; Han, Y.; He, Q.Y. clusterProfiler: An R package for comparing biological themes among gene clusters. OMICS 2012, 16, 284–287. [Google Scholar] [CrossRef]
  60. Huang da, W.; Sherman, B.T.; Lempicki, R.A. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat. Protoc. 2009, 4, 44–57. [Google Scholar] [CrossRef]
  61. Supek, F.; Bosnjak, M.; Skunca, N.; Smuc, T. REVIGO summarizes and visualizes long lists of gene ontology terms. PLoS ONE 2011, 6, e21800. [Google Scholar]
  62. Shirafuta, Y.; Tamura, I.; Ohkawa, Y.; Maekawa, R.; Doi-Tanaka, Y.; Takagi, H.; Mihara, Y.; Shinagawa, M.; Taketani, T.; Sato, S.; et al. Integrated Analysis of Transcriptome and Histone Modifications in Granulosa Cells During Ovulation in Female Mice. Endocrinology 2021, 162, bqab128. [Google Scholar] [CrossRef] [PubMed]
Ijms 23 06731 g001 550
Figure 1. Association between previous fertilization rate and the effect of melatonin. Effect of melatonin on the fertilization rate. The fertilization rates in the cycles before and under melatonin treatment are shown in the violin plot. Left graph: Effect of melatonin in total 75 cases. Middle graph: Effect of melatonin in 19 cases whose fertilization rate before melatonin treatment was less than 50% (Poor fertilization rate cases). Right graph: Effect of melatonin in 56 cases whose fertilization rate before melatonin treatment was more than 50% (Good fertilization rate cases). * p < 0.01 vs. before melatonin.
Figure 1. Association between previous fertilization rate and the effect of melatonin. Effect of melatonin on the fertilization rate. The fertilization rates in the cycles before and under melatonin treatment are shown in the violin plot. Left graph: Effect of melatonin in total 75 cases. Middle graph: Effect of melatonin in 19 cases whose fertilization rate before melatonin treatment was less than 50% (Poor fertilization rate cases). Right graph: Effect of melatonin in 56 cases whose fertilization rate before melatonin treatment was more than 50% (Good fertilization rate cases). * p < 0.01 vs. before melatonin.
Ijms 23 06731 g001
Ijms 23 06731 g002 550
Figure 2. Association between previous embryo development rate (EDR) and the effect of melatonin. (A) Effect of melatonin on the formation rate of good quality embryos. The formation rates of good quality embryos in the cycles before and under melatonin treatment are shown in the violin plot. Left graph: Effect of melatonin in total 56 cases. Middle graph: Effect of melatonin in 41 cases with poor embryo development rate (EDR) (either the formation rate of good quality cleavage-stage embryos or the formation rate of blastocysts before melatonin treatment was less than 50%). Right graph: Effect of melatonin in 15 cases with good EDR. (B) Effect of melatonin on the formation rate of blastocysts. The formation rates of blastocysts in the cycles before and under melatonin treatment are shown in the violin plot. Left graph: Effect of melatonin in total 38 cases. Middle graph: Effect of melatonin in 30 cases with poor EDR. Right graph: Effect of melatonin in 8 cases with good EDR. * p < 0.05 vs. before melatonin.
Figure 2. Association between previous embryo development rate (EDR) and the effect of melatonin. (A) Effect of melatonin on the formation rate of good quality embryos. The formation rates of good quality embryos in the cycles before and under melatonin treatment are shown in the violin plot. Left graph: Effect of melatonin in total 56 cases. Middle graph: Effect of melatonin in 41 cases with poor embryo development rate (EDR) (either the formation rate of good quality cleavage-stage embryos or the formation rate of blastocysts before melatonin treatment was less than 50%). Right graph: Effect of melatonin in 15 cases with good EDR. (B) Effect of melatonin on the formation rate of blastocysts. The formation rates of blastocysts in the cycles before and under melatonin treatment are shown in the violin plot. Left graph: Effect of melatonin in total 38 cases. Middle graph: Effect of melatonin in 30 cases with poor EDR. Right graph: Effect of melatonin in 8 cases with good EDR. * p < 0.05 vs. before melatonin.
Ijms 23 06731 g002
Ijms 23 06731 g003a 550Ijms 23 06731 g003b 550
Figure 3. Effect of melatonin treatment on transcriptome profile in GCs. (A) Hierarchical cluster analysis comparing the transcriptome in GCs. Each color indicates the cases before melatonin treatment (Cases No. 1-3, blue), under melatonin treatment (Cases No. 1-3, green) and the cases with good ART outcomes without melatonin treatment (Cases No. 4-7, yellow). Distances of gene expression pattern are indicated as height. (B) Enriched gene ontology terms and p-values in down-regulated genes by melatonin treatment. The terms circled in red are associated with cell death, and the terms circled in blue are associated with T-cell activation. (C) Enriched gene ontology terms and p-values in up-regulated genes by melatonin treatment. The terms circled in green are associated with steroidogenesis, and the terms circled in yellow are associated with angiogenesis. (D) Changes in mRNA expression level of the representative down-regulated genes by melatonin treatment. Expression values before and under melatonin treatment are shown as dot plots of TPM determined by RNA-sequence. Blue line and dots; Case 1. Orange line and dots; Case 2. Gray line and dots; Case 3. Gene names and associated cellular functions are shown above the graphs. (E) Changes in mRNA expression level of the representative up-regulated genes by melatonin treatment. Expression values before and under melatonin treatment are shown as dot plots of TPM determined by RNA-sequence. Blue line and dots; Case 1. Orange line and dots; Case 2. Gray line and dots; Case 3. Gene symbols and associated cellular functions are shown above the graphs.
Figure 3. Effect of melatonin treatment on transcriptome profile in GCs. (A) Hierarchical cluster analysis comparing the transcriptome in GCs. Each color indicates the cases before melatonin treatment (Cases No. 1-3, blue), under melatonin treatment (Cases No. 1-3, green) and the cases with good ART outcomes without melatonin treatment (Cases No. 4-7, yellow). Distances of gene expression pattern are indicated as height. (B) Enriched gene ontology terms and p-values in down-regulated genes by melatonin treatment. The terms circled in red are associated with cell death, and the terms circled in blue are associated with T-cell activation. (C) Enriched gene ontology terms and p-values in up-regulated genes by melatonin treatment. The terms circled in green are associated with steroidogenesis, and the terms circled in yellow are associated with angiogenesis. (D) Changes in mRNA expression level of the representative down-regulated genes by melatonin treatment. Expression values before and under melatonin treatment are shown as dot plots of TPM determined by RNA-sequence. Blue line and dots; Case 1. Orange line and dots; Case 2. Gray line and dots; Case 3. Gene names and associated cellular functions are shown above the graphs. (E) Changes in mRNA expression level of the representative up-regulated genes by melatonin treatment. Expression values before and under melatonin treatment are shown as dot plots of TPM determined by RNA-sequence. Blue line and dots; Case 1. Orange line and dots; Case 2. Gray line and dots; Case 3. Gene symbols and associated cellular functions are shown above the graphs.
Ijms 23 06731 g003aIjms 23 06731 g003b
Table
Table 1. Backgrounds of patients.
Table 1. Backgrounds of patients.
Poor Fertilization Rate CasesGood Fertilization Rate CasesPoor EDR CasesGood EDR Cases
No. of patients19564115
Patient age (mean ± SD)35.4± 3.737.3 ± 4.336.9 ± 4.238.3 ± 4.1
Patient’s complications
Myoma/adenomyosis (%)15.816.126.712.2
Endometriosis (%)26.316.113.317.1
PCOS (%)5.35.40.07.3
Tubal factor (%)21.121.413.324.4
Unilateral obstruction (%)10.512.513.312.2
Bilateral obstruction (%)10.58.90.012.2
Male factor (%)15.816.126.712.2
Asthenozoospermia (%)5.38.913.37.3
Oligospermia (%)10.57.113.34.9
EDR, embryo development rate
Table
Table 2. Information on cases provided for transcriptome analysis.
Table 2. Information on cases provided for transcriptome analysis.
Blastocyst Formation Rate (%)GnRH Agonist ProtocolFertilization
Case No. AgeBefore Melatonin TreatmentUnder Melatonin TreatmentBefore Melatonin TreatmentUnder Melatonin TreatmentBefore Melatonin TreatmentUnder Melatonin Treatment
132 040.0long protocollong protocolICSIICSI
233 028.9long protocollong protocolcIVFcIVF
336 0100.0long protocollong protocolICSIICSI
435 61.5 long protocol cIVF
536 50.0 long protocol cIVF
634 50.0 short protocol cIVF
737 100.0 long protocol cIVF
long protocol, GnRHa (900 μg/day) was given from the mid-luteal phase in the previous cycle until two days before oocyte retrieval.
cIVF, conventional in vitro fertilizationICSI, intracytoplasmic sperm injection
Table
Table 3. List of top DEGs.
Table 3. List of top DEGs.
Down-Regulated GenesUp-Regulated Genes
Genep ValueGenep Value
NUP62CL1.08 × 10−03RP11-211N8.33.88 × 10−04
TRAF61.40 × 10−03SNRPGP29.07 × 10−04
ZNF5011.98 × 10−03RPS17P59.07 × 10−04
RP11-355I22.72.26 × 10−03RPS11P79.07 × 10−04
RP11-3N2.133.29 × 10−03CTD-2013M15.19.94 × 10−04
ZDHHC155.43 × 10−03COX7B1.06 × 10−03
CENPV5.82 × 10−03RP4-738P15.11.36 × 10−03
RP11-564C4.66.78 × 10−03RP11-787I22.31.38 × 10−03
CD2446.83 × 10−03RP11-874J12.41.39 × 10−03
AC078883.37.24 × 10−03CTC-507E2.11.68 × 10−03
ZNF8357.36 × 10−03SMIM112.09 × 10−03
RP11-63E5.67.54 × 10−03GJB73.60 × 10−03
MMP97.63 × 10−03PTTG14.90 × 10−03
RP11-768F21.18.44 × 10−03CXCL65.12 × 10−03
AC093668.28.57 × 10−03RP11-315I20.15.28 × 10−03
RP11-73M18.71.09 × 10−02ACOX25.66 × 10−03
AC061975.71.14 × 10−02MAGOH26.78 × 10−03
AC090673.21.14 × 10−02RP11-440L14.16.78 × 10−03
CTD-2514K5.41.14 × 10−02ENPP7P107.43 × 10−03
MRPS21P31.14 × 10−02RPL7AP307.66 × 10−03
RN7SL474P1.14 × 10−02RP11-514P8.77.91 × 10−03
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Tamura, I.; Tamura, H.; Kawamoto-Jozaki, M.; Shirafuta, Y.; Fujimura, T.; Doi-Tanaka, Y.; Mihara, Y.; Taketani, T.; Sugino, N. Effects of Melatonin on the Transcriptome of Human Granulosa Cells, Fertilization and Blastocyst Formation. Int. J. Mol. Sci. 2022, 23, 6731. https://doi.org/10.3390/ijms23126731

AMA Style

Tamura I, Tamura H, Kawamoto-Jozaki M, Shirafuta Y, Fujimura T, Doi-Tanaka Y, Mihara Y, Taketani T, Sugino N. Effects of Melatonin on the Transcriptome of Human Granulosa Cells, Fertilization and Blastocyst Formation. International Journal of Molecular Sciences. 2022; 23(12):6731. https://doi.org/10.3390/ijms23126731

Chicago/Turabian Style

Tamura, Isao, Hiroshi Tamura, Mai Kawamoto-Jozaki, Yuichiro Shirafuta, Taishi Fujimura, Yumiko Doi-Tanaka, Yumiko Mihara, Toshiaki Taketani, and Norihiro Sugino. 2022. "Effects of Melatonin on the Transcriptome of Human Granulosa Cells, Fertilization and Blastocyst Formation" International Journal of Molecular Sciences 23, no. 12: 6731. https://doi.org/10.3390/ijms23126731

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop