Next Article in Journal
Regulation of Adipose Progenitor Cell Expansion in a Novel Micro-Physiological Model of Human Adipose Tissue Mimicking Fibrotic and Pro-Inflammatory Microenvironments
Next Article in Special Issue
A Mathematical Model of Estradiol Production from Ultrasound Data for Bovine Ovarian Follicles
Previous Article in Journal
Regulation of Conidiogenesis in Aspergillus flavus
Previous Article in Special Issue
The Interplay between Telomeres, Mitochondria, and Chronic Stress Exposure in the Aging Egg
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Cells
Volume 11
Issue 18
10.3390/cells11182797
cells-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

MicroRNA-4516 in Urinary Exosomes as a Biomarker of Premature Ovarian Insufficiency

by
Zobia Umair
1,†,
Mi-Ock Baek
1,2,†,
Jisue Song
3,
Seona An
3,
Seung Joo Chon
3,* and
Mee-Sup Yoon
1,2,4,*
1
Department of Molecular Medicine, School of Medicine, Gachon University, Incheon 21999, Korea
2
Lee Gil Ya Cancer and Diabetes Institute, Gachon University, Incheon 21999, Korea
3
Gachon University Gil Medical Center, Department of Obstetrics and Gynecology, College of Medicine, Gachon University, Incheon 21565, Korea
4
Department of Health Sciences and Technology, GAIHST, Gachon University, Incheon 21999, Korea
*
Authors to whom correspondence should be addressed.
†
These authors contributed equally to this study.
Cells 2022, 11(18), 2797; https://doi.org/10.3390/cells11182797
Submission received: 28 July 2022 / Revised: 26 August 2022 / Accepted: 29 August 2022 / Published: 7 September 2022
(This article belongs to the Special Issue From Development to Death: Molecular Pathways inside the Oocyte)
Browse Figures
Review Reports Versions Notes

Abstract

:
Premature ovarian insufficiency (POI) is a typical disorder of amenorrhea that lasts for a minimum of four months in women < 40 years old and is typically characterized by reduced estrogen levels and elevated serum concentrations of follicle-stimulating hormone. We collected urine samples from two participant cohorts from Gil Hospital of Gachon University (Incheon, Korea): a sequencing cohort of 19 participants (seven patients with POI (POI patients without Turner syndrome), seven patients with Turner syndrome (POI patients with Turner syndrome), and five control individuals (age-matched controls with confirmed ovarian sufficiency)) and a validation cohort of 46 participants (15 patients with POI, 11 patients with Turner syndrome, and 20 control individuals). Among differentially expressed miRNAs, hsa-miR-4516 was significantly upregulated in patients with POI in both cohorts, independent of the presence of Turner syndrome. Moreover, the upregulation of miR-4516 was confirmed in the ovary—but not in the uterus—of a cyclophosphamide and busulfan-induced POI mouse model. This was accompanied by a decrease in STAT3 protein level, a predicted target of miR-4516, via miRTarBase2020. Our study provides compelling evidence that miR-4516 is highly expressed in patients with POI and POI mouse models, suggesting that miR-4516 is a diagnostic marker of POI.

1. Introduction

Premature ovarian insufficiency (POI) is a pathological condition of ovarian reserve exhaustion before the age of 40 years, manifesting as amenorrhea or oligomenorrhea, hypoestrogenism, and elevated serum follicle-stimulating hormone concentrations (>25 mIU/mL) [1]. The prevalence of POI is approximately 1%, with some variation depending on ethnicity [2]. Patients with POI have an increased risk of cardiovascular disease [3], osteoporosis [4], and cognitive impairment [5]. They may also have elevated total and cancer-specific mortality rates [6]. POI is a highly heterogeneous disease and originates from iatrogenic, karyotypic, and genetic factors. However, the etiology of POI remains unknown in approximately 90% of cases. The occurrence of POI in young women of reproductive potential is difficult to predict because clinical symptoms preceding the disease have not been reported. In patients with early ovarian insufficiency, attempts at pregnancy are made through re-transplantation and subsequent in vitro activation treatment after an ex vivo culture by harvesting whole ovaries or ovarian tissues/follicles. However, patients with POI whose ovarian function cannot be restored are unlikely to have successful results with these methods. Therefore, a diagnostic approach is sought to screen for the abnormal deterioration of ovarian function and advise patients with high risks of POI for family planning.
Turner syndrome is the most frequent hereditary cause of POI and is occasionally observed in women with POI, with symptoms typically appearing before menarche [7]. In karyotyping, Turner syndrome is characterized by an entire or partial deletion of one X-chromosome, resulting in oocyte loss throughout childhood [8]. Together with infertility, Turner syndrome is also associated with short height, delayed puberty, ovarian dysgenesis, hypergonadotropic hypogonadism, congenital cardiac abnormalities, endocrine disorders (such as type 1 and type 2 diabetes), osteoporosis, and immunological problems [9]. Most women with Turner syndrome have primary or secondary hypergonadotropic hypogonadism, which necessitates hormone replacement therapy [9]. Management of POI would benefit from a diagnostic marker that distinguishes POI induced by Turner syndrome from POI produced by other causes.
Exosomes are functional vehicles that transport complex proteins, lipids, and nucleic acids between cells [10,11,12,13]. They are endosomal in origin, range in size from 30 to 200 nm, and are secreted from a given parental cell [14]. Exosomes can affect the behavior of nearby cells, their parental cell microenvironment, and the phenotype of distant cells and tissues, thereby causing systemic consequences [15]. Because exosomes may be found in most body fluids, analyzing exosomes in liquid samples collected by using non-invasive or minimally invasive procedures could detect pathological alterations that would otherwise necessitate a tissue biopsy [15]. Exosomes contain numerous microRNAs (miRNAs) [16]; miRNAs are single-stranded, noncoding 17–24-nucleotides-long RNAs that influence post-transcriptional gene silencing by binding to the 3′-untranslated regions (UTRs) or open reading frames of their target mRNAs [17]. These exosomes may include a disease-specific miRNA signature that can be a valuable diagnostic tool [16].
Besides their detection in ovaries, miRNAs derived from liquid biopsy samples, such as blood plasma and serum, have been evaluated for their function in POI [18,19]. Since the various homeostasis mechanisms may reflect changes in the blood and, consequently, changes in miRNA expressions, the urine may represent a better choice as a non-invasive liquid biopsy sample [20]. Moreover, urine samples have additional benefits of ease of handling, low cost, and availability in large volumes. However, miRNA expression profiling in urine as a possible biomarker for POI screening, early identification, or prognosis has not been investigated.
In this study, we identified urinary exosome miRNAs specific for POI, which can be used as potential diagnostic markers in patients with POI. We analyzed the miRNA expression profile of POI patients with or without Turner syndrome, isolated urinary exosomes, selected putative miRNAs, and verified them in a larger patient cohort and a POI mouse model. Considering that early detection of POI is vital for preserving fertility, diagnostic methods for screening patients with POI can improve the management of the condition. Herein, we have also identified miRNA biomarkers that can be potentially used to develop rapid and straightforward diagnostic methods for POI screening.

2. Materials and Methods

2.1. Sample Collection and Characterization

All patients with POI/Turner syndrome and control individuals were women who visited the outpatient clinic in the Department of Gynecology of Gil Hospital, Gachon University (Incheon, Korea), from January 2019 to August 2021. Women younger than 40 years who did not have menstruation for at least six months with persistent hypergonadotropic hypogonadism in four-week intervals were diagnosed with POI. Women diagnosed with POI underwent serum karyotyping to exclude X-chromosome deletions and translocations, which can be inherited and cause POI symptoms in women after age 30. This karyotyping procedure was also performed to screen for Y chromosomes in women under 30. The patients underwent serum fragile X mental retardation 1 (FMR1) gene mutation studies. The duration of amenorrhea was defined as the period of no menstruation, calculated from the day of last menstruation until the first visit to the outpatient clinic. The control group included non-pregnant women with regular menstrual cycles and without abnormal gynecologic ultrasonography findings. None of the patients had a history of drug intake affecting glucose and lipid metabolism and no known medical conditions or diseases, such as polycystic ovarian syndrome, Cushing’s syndrome, congenital adrenal hyperplasia, androgen-secreting tumors, or endometriosis. The patients with POI were classified into two groups: those who were karyotypically confirmed as having Turner syndrome (Turner syndrome group) and those who did not (POI group). Blood samples were collected on the day of POI diagnosis before administration of any hormonal medication to the POI and Turner syndrome groups and the second to the fifth day of the menstrual cycle from the women in the control group. This study initially collected urine samples from seven patients with POI, seven patients with Turner syndrome, and five control individuals in the sequencing cohort. Additional samples were collected from 15 patients with POI, 11 patients with Turner syndrome, and 20 control individuals in the validation cohort. If the participants underwent their menstrual cycles, they were requested to collect morning urine samples on their non-menstrual cycle days.

2.2. Urine Processing and Exosome Isolation

Morning urine of the patients was processed within 4 h after voiding, as no degradation was observed within this period in pre-trials. Approximately 50–100 mL of urine was transferred to tubes and centrifuged for 10 min at 2000× g to remove any cells or cell debris. The supernatant was transferred to fresh tubes and stored at −80 °C or processed to isolate exosomes within a week. Exosomes from the urine were isolated by using an ExoQuick-TC Exosome Precipitation Solution Kit (EXOTC10A-1, System Biosciences, Berkeley, CA, USA). In brief, 250 μL ExoQuick solution was added to 1 mL of concentrated urine (1:4), mixed thoroughly, and incubated overnight at 4 °C. Then the mixture was centrifuged at 1500× g for 30 min at 4 °C, and the supernatant was removed. The exosome pellet was resuspended in 200 μL of PBS.

2.3. RNA Isolation

Total RNA was extracted from exosomes by using TRIzol reagent (15596026, Thermo Fisher Scientific, Waltham, MA, USA), according to the manufacturer’s instructions. RNA quality was assessed by using an Agilent 2100 bioanalyzer containing an RNA 6000 Pico Chip (Agilent Technologies, Amstelveen, The Netherlands), and RNA quantification was performed by using a Nanodrop 2000 Spectrophotometer system (Thermo Fisher Scientific).

2.4. MicroRNA Library Preparation and Sequencing

RNA sequencing was performed by E-Biogen, Inc. (Seoul, Korea). A microRNA library was constructed by using the NEB Next Multiplex Small RNA Library Prep Kit (E7330L, New England BioLabs, Inc., Beverly, MA, USA) according to the manufacturer’s instructions. Briefly, during library construction, 1 μg of total RNA from each sample was used to ligate the adaptors, and cDNA was synthesized by using reverse transcriptase with adaptor-specific primers. PCR was performed to ensure library amplification, and the libraries were cleaned by using a QIAquick PCR Purification Kit (Qiagen, Inc., Hilden, Germany) and AMPure XP beads (Beckman Coulter, Inc., Brea, CA, USA). The yield and size distribution of the small RNA libraries were assessed by using an Agilent 2100 Bioanalyzer for the high-sensitivity DNA assay (Agilent Technologies, Inc., Santa Clara, CA, USA). High-throughput sequences were produced by using a Next Seq500 system with single-end 75 sequencing (Illumina, San Diego, CA, USA).

2.5. Data Analysis

Sequence reads were mapped by using the Bowtie2 software tool to obtain a bam file (alignment file). Mature miRNA sequences were used as a reference for mapping. Read counts mapped on mature miRNA sequences were extracted from the alignment file by using Bedtools (v2.25.0, Salt Lake City, UT, USA) [21] and Bioconductor [22], which use the R statistical programming language (version 3.2.2; R Development Core Team, 2011, Vienna, Austria). Read counts were used to determine the miRNA expression levels. The quantile normalization method was used for comparisons between samples. In the miRNA target study, miRWalk 2.0 [23] was performed. Functional gene classification was performed by using DIANA-miRPath v3.0 (Lamia, Greece) [24].

2.6. Transmission Electron Microscopy

The isolated exosome samples were applied to copper grids coated solely with a thin carbon foil (01340, Ted Pella, Inc. Redding, CA, USA). After allowing the sample to absorb for 2 min and blotting off excess buffer solution onto Whatman filter paper, the samples were stained with 2% (w/v) uranyl acetate for 1 min. Distilled water was added for 1 min to remove the excess uranyl acetate, followed by drying for 15 min. The results were recorded by using a Bio-High Voltage EM system (JEM-1400 Plus at 120 kV and JEM-1000BEF at 1000 kV; JEOL Ltd., Tokyo, Japan) at the Korea Basic Science Institute in the Korea.

2.7. Exosome Particle Size Analysis

The size distribution of the isolated exosomes was evaluated by using a particle analyzer (Litesizer 500; Anton Paar, Graz, Austria). PBS was used to dilute the samples. The size of each sample was measured for 60 s in triplicate, and the captured data were analyzed by using Excel.

2.8. Western Blot Analysis

Exosomes were lysed, using a lysis buffer (9803, Cell Signaling Technology, Danvers, MA, USA), for protein analysis. The supernatant was collected following microcentrifugation at 13,000× g for 10 min and boiled for 5 min in a sodium dodecyl sulfate sample buffer [25]. Proteins in the mouse ovaries were extracted by using a T-PER tissue protein extraction reagent supplemented with protease (05892791001, Roche, Mannheim, Germany) and phosphatase inhibitors (4906837001, Roche, Mannheim, Germany). Ovary tissues (30 mg) were homogenized in an extraction buffer, using steel beads and a TissueLyser II (Qiagen, Hilden, Germany) for 1–2 min at 30 Hz. The lysates were centrifuged at 13,000× g for 10 min at 4 °C to remove debris. The supernatant was mixed with a sodium dodecyl sulfate sample buffer and boiled for 5 min. Proteins were subjected to electrophoresis on sodium dodecyl sulfate–polyacrylamide gels and then transferred to polyvinylidene fluoride membranes (IPVH00010, Millipore, Billerica, MA, USA). The membranes were incubated with the primary antibodies following the manufacturer’s instructions. Horseradish peroxidase-conjugated secondary antibodies were detected by using an Immobilon Western Chemiluminescent HRP Substrate (WBKLS0500, Millipore, Billerica, MA, USA). The band intensity of the Western blot was quantified via densitometry, using ImageJ software. The following antibodies were used: anti-TSG101 (ab125011, 1:1000) and anti-tubulin (ab11304, 1:2000) from Abcam (Cambridge, MA, USA); anti-CD9 (sc-13118, 1:200), anti-CD63 (sc-5275, 1:200), anti-STAT3 (sc-8019, 1:500), and anti-p53 (sc-126, 1:500) from Santa Cruz Biotechnology (Dallas, TX, USA), and anti-Bax (2772, 1:1000) antibody from Cell Signaling Technology. Secondary antibodies (anti-mouse 115-035-003, anti-rabbit 211-002-171) were obtained from Jackson ImmunoResearch Laboratories Inc. (West Grove, PA, USA).

2.9. Quantitative Real-Time Polymerase Chain Reaction

Total RNA was extracted from the exosomes and ovary tissues by using TRIzol reagent (15596026, Thermo Fisher Scientific), according to the manufacturer’s instructions. The miRCURY LNA RT kit (339340, QIAGEN, Hilden, Germany) was used to synthesize cDNA from 5 ng of total RNA by using a 10× miRCURY RT Enzyme Mix with UniSp6 RNA spike-in (QIAGEN) in 20 μL reaction volumes. Quantitative real-time polymerase chain reaction (qRT-PCR) analysis was performed by using a CFX384 C1000 Thermal Cycler (Bio-Rad, Hercules, CA, USA) with a 2× miRCURY SYBER Green Master Mix (QIAGEN). RNA expression levels were normalized against UniSp6 [26]. The sequences of the primers used in this study are listed in Table 1.
For real-time TaqMan analysis of hsa-miR-4516, a Taqman advanced miRNA cDNA synthesis kit (A28007, Applied Biosystems, Foster City, CA, USA) was used to synthesize cDNA for miRNAs by using a multiplex reverse transcription technique with TaqMan® advanced microRNA assays (A25576, Applied Biosystems) for hsa-miR-4516 (478308) and ath-miR-159a (4427975).

2.10. Establishment of A POI Mouse Model

Eight-week-old female pathogen-free-grade C57B6L/N mice were obtained from Orient Bio Inc. (Seongnam, Korea). The mice were kept at a controlled temperature (20–22 °C) and lighting (12/12 h light/dark cycle) and were allowed free access to water and food. Experimental animals were randomly divided into two groups, namely control (n = 7) and cyclophosphamide (CTX)+ busulfan (BUS)-POI (n = 7), and three independent experiments were performed. To induce POI, the mice were intraperitoneally injected with 120 mg/kg CTX and 12 mg BUS/kg weekly for two weeks. The control group was injected with saline. CTX (C0768) and BUS (B2635) were purchased from Sigma-Aldrich (St. Louis, MO, USA).

2.11. Hematoxylin-and-Eosin Staining and SIRIUS Red Staining of the Mouse Ovaries

The ovary tissues were fixed in 4% paraformaldehyde overnight, dehydrated, embedded in paraffin, and cut into 5 μm-thick sections. Every fifth section was stained with hematoxylin and eosin, as described previously [27], and assessed total and atretic follicles from the entire area of a cross-sectioned ovary section from each mouse (Saline, n = 7; CTX-BUS, n = 7). Sirius red staining was performed to determine the presence of collagen fibers in the ovarian tissues. Images of hematoxylin-and-eosin and Sirius-red staining were randomly captured by using a Motic Easyscan Digital Slide Scanner (Motic Hong Kong Limited, Hong Kong, China).

2.12. TUNEL Assay and Ki67 Immunohistochemistry

Apoptosis in ovary tissues was detected by using a DeadEnd™ Fluorometric TUNEL System (G3250, Promega, Fitchburg, WI, USA), following the manufacturer’s instructions. The slides were examined by using an Olympus CKX3-Houn microscope (Olympus, Tokyo, Japan) after preparation with a mounting medium for fluorescence with DAPI (Vector Lab, Newark, CA, USA). The number of TUNEL-positive cells was counted from three random fields (500 μm × 500 μm) from each ovary cross-section and averaged for each mouse. For Ki67 immunohistochemistry, antigen retrieval was performed by incubating the sections in 0.01 M citrate buffer (pH 6.0) and applying intense microwave irradiation for 30 min following deparaffinization and rehydration. Endogenous peroxidase activity was quenched for 10 min with 3% H2O2, and nonspecific binding was blocked for 1 h with 10% normal goat serum. The sections were treated with anti-Ki67 primary antibodies (1:200; ab16667, Abcam) for 24 h at 4 °C. Antigens were visualized by using diaminobenzidine. Finally, the slides were dehydrated and mounted after counterstaining with hematoxylin. Images of Ki67 immunohistochemistry were randomly captured by using a Motic Easyscan Digital Slide Scanner (Motic Hong Kong Limited).

2.13. Statistical Analyses

The baseline clinical characteristics of the included patients were analyzed by using a one-way analysis of variance in IBM SPSS Statistics software, version 28.0. The mean ± standard deviation of at least three independent experiments was calculated (three to seven independent experiments for all figures). Dots represent the individual data points in all the quantitative data graphs. A two-tailed paired Student’s t-test in Excel or GraphPad Prism version 9.0 for Microsoft Windows (GraphPad Software, CA, USA) was used to evaluate the statistical significance of the data when necessary (GraphPad Software, CA, USA). The threshold for statistical significance was set at p < 0.05.

3. Results

3.1. Clinical Characteristics

We selected a sequencing cohort (7 patients with POI (POI patients without Turner syndrome), 7 patients with Turner syndrome (POI patients with Turner syndrome), and 5 control individuals (age-matched controls with confirmed ovarian sufficiency) and a validation cohort (15 patients with POI, 11 patients with Turner syndrome, and 20 control individuals). Table 1 summarizes the clinical features of all the 65 patients in this study. There were no statistically significant differences between groups in age, weight, body mass index, diastolic blood pressure, aspartate aminotransferase, alanine aminotransferase, total cholesterol (TC), triglyceride, low-density lipoprotein cholesterol, or prolactin levels. Patients with Turner syndrome were significantly shorter in height (p < 0.001) and had higher high-density lipoprotein cholesterol (p = 0.013) and follicle-stimulating hormone (p < 0.001) levels (Table 2). Only 2 out of 15 patients with POI patients in the validation cohort and all patients with Turner syndrome had primary amenorrhea (Table 2 and Table 3). Among the patients with secondary amenorrhea, the mean duration of amenorrhea was 37.8 months. A total of 3 of the POI patients of the validation cohort had chromosomal mosaicism, whereas the remaining 12 patients had a normal 46,XX karyotype. None of the patients with POI who did not have Turner syndrome had FMR1 gene mutations (Table 3 and Table 4).

3.2. Identification of Isolated Exosomes in the Urine

To characterize exosomes in human urine, we assessed their morphology and diameter, as well as the presence of exosome-enriched protein markers such as CD9, CD63, and TSG101. Transmission electron microscopy showed that urine exosomes exhibited round or cup-shaped morphologies, with diameters ranging from 95 to 160 nm (Figure 1A). According to nanoparticle tracking analysis, isolated exosomes showed a median diameter of 172.0515 nm, with the highest peak at 5–6% of the expected distribution profiles, as is consistent with pure exosome preparations (Figure 1B). No significant difference was observed in the frequency of exosomes between the groups (Figure 1B). The number of exosomes per unit of urine in the control was comparable to that in patients, as evidenced by the amount of total RNA (Figure 1C). The Western blot analysis revealed the presence of CD9, CD63, and TSG101 (Figure 1D).

3.3. Differential Expression of Urinary Exosome miRNAs in Patients with POI, Patients with Turner Syndrome, and Control Individuals

We performed miRNA deep sequencing by using urinary exosomes from patients with POI (n = 7), patients with Turner syndrome (n = 7), and the control individuals (n = 5) to examine the miRNA expression profiles of urinary exosomes in respective patients. We found 2588 miRNAs in human urinary exosomes from that sample set. The miRNA expression profiles were comparatively analyzed across the three groups (e.g., miRNAs expressed in the POI group versus the control group, the Turner syndrome group versus the control group, the Turner syndrome group versus the POI group, and both POI and Turner syndrome group versus the control group). We identified 33 differentially expressed miRNAs, with differential expression between any given two groups defined as a statistically significant difference in expression (p < 0.05) of at least 1.3 times, with a normalization value of 2 (Table 5).
The heatmap and scatter plots showed differentially expressed miRNAs in POI-, Turner-syndrome-, and control-group-specific clusters (Figure 2A–C). The heatmap showed that hsa-miR-4516 was the most highly upregulated miRNA in POI samples compared to the control individuals. Moreover, hsa-miR-4516 was significantly upregulated in patients with Turner syndrome with lower levels than in patients with POI. To understand the biological roles of the predicted target genes of the differentially expressed miRNAs, a Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis was performed. This process uncovers the primary pathways in which the prospective target genes may be involved. Figure 2D shows the top 20 most significantly enriched pathways. The top five pathways are extracellular matrix–receptor interaction, Hippo signaling pathway, fatty acid biosynthesis, renal cell carcinoma, and proteoglycans in cancer. In addition, a gene category analysis of selected miRNAs demonstrated that most of the selected miRNAs are involved in uterine embryonic development, ovarian follicle development, and antral follicle development; all of these processes are closely related to the development of POI (Figure 2E).

3.4. Validation of the Differentially Expressed miRNAs via qRT-PCR

To further validate the miRNA profiling results of our RNA sequencing analysis, urinary exosomes isolated from the validation cohort (15 patients with POI, 11 patients with Turner syndrome, and 20 control individuals) were used to confirm the expression of selected miRNAs via miRNA sequencing. From the generated heatmap and significance plot, we selected 12 miRNAs that belonged to the following categories: Group 1, miRNAs that were highly differentially expressed in the POI or Turner syndrome group; and Group 2, substantial differences in expression in the POI or Turner syndrome group, but evenly expressed within the same patient group (e.g., the POI group or Turner syndrome group) (Table S1). Quantitative RT-PCR results confirmed the substantial expression of 12 miRNAs in the validation cohort. Comparisons of the relative expression levels of selected miRNAs between the sequencing cohort and the validation cohort are shown in Table 6 and Table 7 and Figure S1A,B. For 10 of the 12 miRNAs selected for the validation analysis, the expression patterns of patients with POI in the validation cohort were the same as the expression patterns of the POI group in the sequencing cohort: hsa-miR-16-5p, hsa-miR-29a-3p, hsa-miR-30b-5p, hsa-miR-99b-3p, hsa-miR-151a-5p, hsa-miR-200b-3p, hsa-miR-423-3p, hsa-miR-941, hsa-miR-4516, and hsa-miR-7847-3p (Table 6 and Figure S1A). Among them, the expression rates of two miRNAs were significantly different between patients with POI and control individuals, as we had originally identified in the sequencing cohort: hsa-miR-4516 was significantly upregulated (p = 0.008), whereas hsa-miR-30b-5p was significantly downregulated (p = 0.026; Table 6 and Figure S1A). Notably, hsa-miR-4516 was significantly upregulated (p = 0.022) in patients with Turner syndrome compared to control individuals to a lesser extent, similar to the RNA sequencing data (Table 7 and Figure S1B). In patients with Turner syndrome, hsa-miR-20a-5p expression was substantially increased compared to the control group (p = 0.447), while hsa-miR-200b-3p expression was substantially decreased (p = 0.296), similar to the sequencing cohort results (Table 7 and Figure S1B). Although the expression pattern was different compared to that of the sequencing cohort, the expressions of hsa-miR-29a-3p (p = 0.021) and hsa-miR-99-3p (p = 0.011) were significantly altered in validation cohort patients with Turner syndrome, compared to the control group (Table 7 and Figure S1B). In conclusion, hsa-miR-4516 expression detected via qRT-PCR analysis of RNA from the validation cohort exactly mapped the RNA sequencing results of the sequencing cohort. Moreover, it was significantly upregulated in patients with POI and Turner syndrome compared to the control group (Figure 3A). The TaqMan PCR assay confirmed the upregulation of hsa-miR-4516 expression in patients with POI and Turner syndrome (Figure 3B).

3.5. Confirmation of miR-4516 Upregulation in the Ovaries of a Chemotherapy-Induced POI Mouse Mode

To evaluate miR-4516 as a biomarker of POI, we established a chemotherapy-induced POI mouse model by injecting the mice with CTX and BUS. There was no apparent death or toxicity in the POI mouse model, but the size of the ovaries was smaller than that in the saline-injected control mice (Figure S2). Compared to the saline-injected control mice, the POI group exhibited significantly reduced follicle numbers, increased ratio of atretic follicles in total follicles, and increased collagen deposition as proven by Sirius red staining (Figure 4A,B). Furthermore, TUNEL-positive cells were evident in the follicles of the POI groups (Figure 4C), whereas Ki67-positive cells were few (Figure 4D), implying that a POI mouse model was successfully established. Next, we assessed the expression of mmu-miR-4516 in the ovaries of the chemotherapy-induced POI mouse model; mmu-miR-4516 expression was significantly increased in the ovaries but not in the uteri in the POI mouse model compared to the control (Figure 4E). These results suggest that upregulation of mmu-miR-4516 expression is associated with pathological changes in the ovaries. To interpret the pathological relevance of mmu-miR-4516 overexpression, we analyzed the targets of miR-4516 by using miRTarBase2020, an experimentally validated miRNA-target interaction database [28]. STAT3 was identified as the target of miR-4516 with solid experimental evidence (Table S2). The 3′-UTR of STAT3 contains two putative hsa-miR-4516 binding sites [29] (Figure S3A). Consistent with this interaction prediction, STAT3 protein expression was significantly decreased in the ovaries of the chemotherapy-induced POI mouse model, together with the increase in the expression of apoptotic markers Bax and p53 (Figure 4F). Taken together, these results suggest that upregulating miR-4516 is closely related to POI pathology in tandem with a decrease in STAT3 expression.

4. Discussion

A longstanding demand exists for non-invasive testing methods for identifying POI biomarkers, which are expected to facilitate convenient screening, early diagnosis, and monitoring therapy responses. In this study, we investigated the clinical utility of urinary exosome-based miRNAs as non-invasive diagnostic and prognostic indicators in patients with POI and Turner syndrome. We first performed miRNA expression profiling in human urinary exosomes of patients with POI and Turner syndrome patients. We identified miR-4516 as a highly upregulated miRNA in the sequencing cohort, validation cohort, and chemotherapy-induced POI mouse model.
The involvement of miRNAs in the development and risk analysis of POI has been previously investigated by using serum and plasma. In contrast to previous studies, we analyzed human urinary exosomes. In this study, 2588 miRNAs were identified from the sequencing cohort, and 33 differentially expressed miRNAs were selected. Among the selected miRNAs, miR-151a-5p, miR-let-7c-5p, and miR-423-5p have been previously associated with POI. Expressions of miR-151a-5p and miR-let-7c-5p were reported in two independent selections of differentially expressed miRNAs in the plasma of patients with POI [18,19], although they changed in opposite trends. A polymorphism in miR-423 (miR-423C>A), a differentially expressed miRNA in our study, has been linked to a higher idiopathic risk of POI in a study involving a Korean cohort [30]. In addition, several previously reported miRNAs, such as miR-23a, miR-27a, miR-22-3p, miR-146a, miR-181a, miR-196a, miR-290-295, and miR-608 [31,32], were not among the selected differentially expressed miRNAs; miR-4516, the highly expressed miRNA in our selection, has not been previously reported in patients with POI or Turner syndrome. In this study, we confirmed the specificity of miR-4516 for POI by using an independent validation cohort and a well-established POI mouse model, thereby supporting miR-4516 as a potential diagnostic biomarker for POI. These discrepancies may have resulted from the different sampling techniques used in this study and the differences in patient grouping. Further investigation with a larger number of patient samples or another independent cohort could help confirm the specificity of miR-4516. In addition, future investigation is required to determine whether miR-4516 levels correlate with various disease-related ovarian dysfunctions such as type 1 diabetes mellitus (T1DM). T1DM shows hypogonadism, reduced ovary size, and decreased oocyte quality [33,34]. Although the mechanisms underlying these effects of diabetes on reproduction are not entirely understood, T1DM seems to impair granulosa cell–oocyte communication, mitochondrial dysfunction during meiosis, apoptosis of cumulus cells, and alteration of DNA methylation status within ovarian follicles [35,36,37,38]. The assessment of whether miR-4516 is a dominantly causative factor of POI would be necessary before using miR-4516 as a diagnostic marker of disease-related ovarian dysfunction.
We also observed inconsistencies between the sequencing and qRT-PCR results. The current diagnostic method of POI could not differentiate the severity or stage of POI, and the cause of POI is highly diverse, which can lead to discrepancies in exosomal miRNA expression between the cohorts. In addition, the sequencing cohort is small and limited compared to the validation cohort, and this may have caused inconsistencies between the cohorts. A future study with a larger cohort size is required for further confirmation.
To our knowledge, this is the first study to implicate hsa-miR-4516 as a highly upregulated miRNA in POI. In addition, miR-4516 has not been previously identified in urinary exosomes, according to a miRTarbase analysis. We predicted STAT3 as a target of miR-4516 in POI, using miRTarBase, an experimentally validated tool (Table S2). There are two complementary regions between hsa-miR-4516 and the 3′-UTR of STAT3; the second miRNA-binding site at 1918–1925 bp of the STAT3 3′-UTR is substantially conserved across different species (Figure S3A) [29]. The cytoplasmic protein STAT3 regulates gene transcription for survival, proliferation, and anti-apoptotic and immune response genes when activated in response to cytokines [39]. We observed that the upregulation of miR-4516 was accompanied by a decrease in STAT3 protein level and an increase in the expression of p53 and Bax in the chemotherapy-induced POI mouse model, thereby implying that miR-4516 might induce apoptosis by controlling STAT3 expression (Supplementary Figure S3B), which is supported by other reports [29,40]. In addition to the role of STAT3 in regulating apoptosis, STAT3 uniquely affects granulosa cells by controlling the coordinated activation of primordial follicles throughout reproductive life. Leukemia inhibitory factor stimulates ovarian follicles in vitro by activating STAT3 and then transcribing and translating SOCS4, which interacts with the cardiac troponin-like cytokine, poly(rC) binding protein 1, and malate dehydrogenase 1 [41]. However, STAT3 possibly contributes to the negative regulation of PTEN-Akt signaling, a critical mediator of the development of primordial and secondary follicles [42], thus leading to apoptosis [43] (Figure S3B). However, we cannot exclude the possibility that decreased STAT3 expression resulted from the side effect of chemotherapy in the POI mice model. Therefore, the decrease in STAT3 protein level in POI transgenic mice models such as Fragile X premutation RNA-knockout mice must be confirmed by further investigation [44].
We identified miRNAs that could differentiate between patients with POI and Turner syndrome. We found that the expressions of hsa-miR-29a-3p and hsa-miR-30b-5p decreased, while that of hsa-miR-4516 increased, in patients with POI, whereas the expression of all three miRNAs increased in patients with Turner syndrome in the validation cohort (Figure 3 and Figure S4). An analysis of the expression patterns of these three miRNAs led us to distinguish patients with POI from those with Turner syndrome. In patients with POI without Turner syndrome, hsa-miR-4516 expression is increased, but the expression of hsa-miR-29a-3p and hsa-miR-30b-5p is decreased; in patients with Turner syndrome, the expression of hsa-miR-4516, hsa-miR-29a-3p, and hsa-miR-30b-5p are all increased. However, we could not examine this possibility with a Turner mouse model including X-monosomic mice (39, XO) that all have a single-sex chromosome because these mice are usually fertile, which is different from their human counterparts [45]. A further expanded human validation cohort is needed to develop a diagnostic method that uses a comparative analysis of the expression of these three miRNAs to screen for POI and Turner syndrome without having to rely on chromosomal examination.
Taken together, we propose that hsa-miR-4516 can be developed as a biomarker of POI within an easy and convenient diagnostic test using non-invasive sampling methods. The evaluation of the miR-4516 levels with the severity of an abnormal decline in ovarian functions is required for its use as a predictive marker for appropriate preparation in advance.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/cells11182797/s1. Table S1. Selected 12 microRNAs from two groups; Table S2. Identification of miR-4516 targets using miRTarBase2020; Figure S1. Comparison of the relative expression patterns of selected microRNAs obtained from the mRNA sequencing and qRT-PCR validation; Figure S2. Representative images of ovaries from saline-injected control (Saline) and cyclophosphamide (CTX) + busulfan (BUS)-injected POI mice (CTX+BUS); Figure S3. Predicted mechanism of the miR-4516-induced regulation of the STAT3 3′-UTR and ovarian function; Figure S4. Relative expression of hsa-miR-29a-3p and hsa-miR-30a-5p in the validation cohorts.

Author Contributions

Conceptualization, S.J.C. and M.-S.Y.; methodology, S.J.C. and M.-S.Y.; data analysis, Z.U., M.-O.B., S.J.C. and M.-S.Y.; investigation, Z.U. and M.-O.B.; establishment of the cohorts and sample collection, J.S., S.A. and S.J.C.; writing—original draft preparation, Z.U., S.J.C., M.-O.B. and M.-S.Y.; writing—review and editing, Z.U., M.-O.B., S.J.C. and M.-S.Y.; supervision, S.J.C. and M.-S.Y.; funding acquisition, M.-O.B., S.J.C. and M.-S.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by a National Research Foundation of Korea (NRF) grant funded by the Korean government (Ministry of Science and ICT; NRF-2020R1C1C1009913, NRF-2021R1A2B5B01002047, NRF-2022R1C1C2004305).

Institutional Review Board Statement

The Ethics Committee of Gil Hospital approved this study under the guidelines and regulations of Gachon University (GDIRB2020-123). All experimental procedures in this study, as well as the protocols for the maintenance and care of the animals, were performed in compliance with the Center of Animal Care and Use, Lee Gil Ya Cancer and Diabetes Institute, Gachon University, and the Institutional Animal Care and Use Committee (permission number: LCDI- 2021-0121).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The datasets GENERATED for this study can be found in the Gene Expression Omnibus, https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE205700 (accessed on 27 July 2022).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Rudnicka, E.; Kruszewska, J.; Klicka, K.; Kowalczyk, J.; Grymowicz, M.; Skorska, J.; Pieta, W.; Smolarczyk, R. Premature ovarian insufficiency-aetiopathology, epidemiology, and diagnostic evaluation. Prz. Menopauzalny 2018, 17, 105–108. [Google Scholar] [CrossRef] [PubMed]
  2. Haller-Kikkatalo, K.; Uibo, R.; Kurg, A.; Salumets, A. The prevalence and phenotypic characteristics of spontaneous premature ovarian failure: A general population registry-based study. Hum. Reprod. 2015, 30, 1229–1238. [Google Scholar] [CrossRef] [PubMed]
  3. Roeters van Lennep, J.E.; Heida, K.Y.; Bots, M.L.; Hoek, A.; collaborators of the Dutch Multidisciplinary Guideline Development Group on Cardiovascular Risk Management after Reproductive, D. Cardiovascular disease risk in women with premature ovarian insufficiency: A systematic review and meta-analysis. Eur. J. Prev. Cardiol. 2016, 23, 178–186. [Google Scholar] [CrossRef] [PubMed]
  4. Asli, I.N.; Fallahian, M.; Seddigh, H.R.; Javadi, H.; Baharfar, N.; Assadi, M. Evaluation of bone mineral density in premature ovarian failure. Hell. J. Nucl. Med. 2010, 13, 261–263. [Google Scholar] [PubMed]
  5. Bove, R.; Secor, E.; Chibnik, L.B.; Barnes, L.L.; Schneider, J.A.; Bennett, D.A.; De Jager, P.L. Age at surgical menopause influences cognitive decline and Alzheimer pathology in older women. Neurology 2014, 82, 222–229. [Google Scholar] [CrossRef]
  6. Wu, X.; Cai, H.; Kallianpur, A.; Li, H.; Yang, G.; Gao, J.; Xiang, Y.B.; Ji, B.T.; Yu, T.; Zheng, W.; et al. Impact of premature ovarian failure on mortality and morbidity among Chinese women. PLoS ONE 2014, 9, e89597. [Google Scholar] [CrossRef]
  7. Chon, S.J.; Umair, Z.; Yoon, M.S. Premature Ovarian Insufficiency: Past, Present, and Future. Front. Cell Dev. Biol 2021, 9, 672890. [Google Scholar] [CrossRef]
  8. Reddy, U.M.; Davis, J.M.; Ren, Z.; Greene, M.F. Opioid Use in Pregnancy, Neonatal Abstinence Syndrome, and Childhood Outcomes: Executive Summary of a Joint Workshop by the Eunice Kennedy Shriver National Institute of Child Health and Human Development, American College of Obstetricians and Gynecologists, American Academy of Pediatrics, Society for Maternal-Fetal Medicine, Centers for Disease Control and Prevention, and the March of Dimes Foundation. Obstet Gynecol 2017, 130, 10–28. [Google Scholar]
  9. Gravholt, C.H.; Viuff, M.H.; Brun, S.; Stochholm, K.; Andersen, N.H. Turner syndrome: Mechanisms and management. Nat. Rev. Endocrinol. 2019, 15, 601–614. [Google Scholar] [CrossRef]
  10. Simpson, R.J.; Lim, J.W.; Moritz, R.L.; Mathivanan, S. Exosomes: Proteomic insights and diagnostic potential. Expert Rev. Proteomics 2009, 6, 267–283. [Google Scholar] [CrossRef]
  11. Valadi, H.; Ekstrom, K.; Bossios, A.; Sjostrand, M.; Lee, J.J.; Lotvall, J.O. Exosome-mediated transfer of mRNAs and microRNAs is a novel mechanism of genetic exchange between cells. Nat. Cell Biol 2007, 9, 654–659. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Vidal, M.; Sainte-Marie, J.; Philippot, J.R.; Bienvenue, A. Asymmetric distribution of phospholipids in the membrane of vesicles released during in vitro maturation of guinea pig reticulocytes: Evidence precluding a role for "aminophospholipid translocase". J. Cell Physiol. 1989, 140, 455–462. [Google Scholar] [CrossRef] [PubMed]
  13. Waldenstrom, A.; Genneback, N.; Hellman, U.; Ronquist, G. Cardiomyocyte microvesicles contain DNA/RNA and convey biological messages to target cells. PLoS ONE 2012, 7, e34653. [Google Scholar] [CrossRef] [PubMed]
  14. Whittle, K.; Kao, S.; Clarke, S.; Grau, G.E.R.; Hosseini-Beheshti, E. Exploring the role of extracellular vesicles and their protein cargo in lung cancer metastasis: A review. Crit. Rev. Oncol. Hematol. 2022, 2022, 103603. [Google Scholar] [CrossRef] [PubMed]
  15. Hu, Q.; Su, H.; Li, J.; Lyon, C.; Tang, W.; Wan, M.; Hu, T.Y. Clinical applications of exosome membrane proteins. Precis. Clin. Med. 2020, 3, 54–66. [Google Scholar] [CrossRef]
  16. Cheng, L.; Sharples, R.A.; Scicluna, B.J.; Hill, A.F. Exosomes provide a protective and enriched source of miRNA for biomarker profiling compared to intracellular and cell-free blood. J. Extracell Vesicles 2014, 3, 743. [Google Scholar] [CrossRef]
  17. Bartel, D.P. MicroRNAs: Genomics, biogenesis, mechanism, and function. Cell 2004, 116, 281–297. [Google Scholar] [CrossRef]
  18. Zhou, Y.; Zhu, Y.; Zhang, S.; Wang, H.; Wang, S.; Yang, X. MicroRNA expression profiles in premature ovarian failure patients and its potential regulate functions. Chin. J. Birth Health Hered. 2011, 19, 20–22. [Google Scholar]
  19. Dang, Y.; Zhao, S.; Qin, Y.; Han, T.; Li, W.; Chen, Z.J. MicroRNA-22-3p is down-regulated in the plasma of Han Chinese patients with premature ovarian failure. Fertil. Steril. 2015, 103, 802–807. [Google Scholar] [CrossRef]
  20. Aftab, M.; Poojary, S.S.; Seshan, V.; Kumar, S.; Agarwal, P.; Tandon, S.; Zutshi, V.; Das, B.C. Urine miRNA signature as a potential non-invasive diagnostic and prognostic biomarker in cervical cancer. Sci. Rep. 2021, 11, 10323. [Google Scholar] [CrossRef]
  21. Quinlan, A.R.; Hall, I.M. BEDTools: A flexible suite of utilities for comparing genomic features. Bioinformatics 2010, 26, 841–842. [Google Scholar] [CrossRef] [PubMed]
  22. Gentleman, R.C.; Carey, V.J.; Bates, D.M.; 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, R80. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Dweep, H.; Sticht, C.; Pandey, P.; Gretz, N. miRWalk--database: Prediction of possible miRNA binding sites by "walking" the genes of three genomes. J. Biomed. Inform. 2011, 44, 839–847. [Google Scholar] [CrossRef]
  24. Vlachos, I.S.; Zagganas, K.; Paraskevopoulou, M.D.; Georgakilas, G.; Karagkouni, D.; Vergoulis, T.; Dalamagas, T.; Hatzigeorgiou, A.G. DIANA-miRPath v3.0: Deciphering microRNA function with experimental support. Nucleic Acids Res. 2015, 43, W460–W466. [Google Scholar] [CrossRef] [PubMed]
  25. Baek, M.O.; Cho, H.J.; Min, D.S.; Choi, C.S.; Yoon, M.S. Self-transducible LRS-UNE-L peptide enhances muscle regeneration. J. Cachexia Sarcopenia Muscle 2022, 13, 1277–1288. [Google Scholar] [CrossRef]
  26. Kinoshita, T.; Yip, K.W.; Spence, T.; Liu, F.-F. MicroRNAs in extracellular vesicles: Potential cancer biomarkers. J. Hum. Genet. 2017, 62, 67–74. [Google Scholar] [CrossRef] [PubMed]
  27. Buigues, A.; Marchante, M.; Herraiz, S.; Pellicer, A. Diminished Ovarian Reserve Chemotherapy-Induced Mouse Model: A Tool for the Preclinical Assessment of New Therapies for Ovarian Damage. Reprod. Sci. 2020, 27, 1609–1619. [Google Scholar] [CrossRef] [PubMed]
  28. Huang, H.Y.; Lin, Y.C.; Li, J.; Huang, K.Y.; Shrestha, S.; Hong, H.C.; Tang, Y.; Chen, Y.G.; Jin, C.N.; Yu, Y.; et al. miRTarBase 2020: Updates to the experimentally validated microRNA-target interaction database. Nucleic Acids Res. 2020, 48, D148–d154. [Google Scholar] [PubMed]
  29. Chowdhari, S.; Saini, N. hsa-miR-4516 mediated downregulation of STAT3/CDK6/UBE2N plays a role in PUVA induced apoptosis in keratinocytes. J. Cell Physiol. 2014, 229, 1630–1638. [Google Scholar] [CrossRef]
  30. Rah, H.; Kim, H.S.; Cha, S.H.; Kim, Y.R.; Lee, W.S.; Ko, J.J.; Kim, N.K. Association of breast cancer–related microRNA polymorphisms with idiopathic primary ovarian insufficiency. Menopause 2015, 22, 437–443. [Google Scholar] [CrossRef]
  31. Guo, Y.; Sun, J.; Lai, D. Role of microRNAs in premature ovarian insufficiency. Reprod. Biol. Endocrinol. 2017, 15, 38. [Google Scholar] [CrossRef] [PubMed]
  32. Zhang, Q.; Sun, H.; Jiang, Y.; Ding, L.; Wu, S.; Fang, T.; Yan, G.; Hu, Y. MicroRNA-181a suppresses mouse granulosa cell proliferation by targeting activin receptor IIA. PLoS ONE 2013, 8, e59667. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Abell, S.K.; Nankervis, A.; Khan, K.S.; Teede, H.J. Type 1 and Type 2 Diabetes Preconception and in Pregnancy: Health Impacts, Influence of Obesity and Lifestyle, and Principles of Management. Semin. Reprod. Med. 2016, 34, 110–120. [Google Scholar] [CrossRef] [PubMed]
  34. Codner, E.; Eyzaguirre, F.C.; Iniguez, G.; Lopez, P.; Perez-Bravo, F.; Torrealba, I.M.; Cassorla, F.; Chilean Group for the Study of Ovarian Function in Type, D. Ovulation rate in adolescents with type 1 diabetes mellitus. Fertil. Steril. 2011, 95, 197–202.e191. [Google Scholar] [CrossRef]
  35. Chang, A.S.; Dale, A.N.; Moley, K.H. Maternal diabetes adversely affects preovulatory oocyte maturation, development, and granulosa cell apoptosis. Endocrinology 2005, 146, 2445–2453. [Google Scholar] [CrossRef]
  36. Diamond, M.P.; Moley, K.H.; Pellicer, A.; Vaughn, W.K.; DeCherney, A.H. Effects of streptozotocin- and alloxan-induced diabetes mellitus on mouse follicular and early embryo development. J. Reprod. Fertil. 1989, 86, 1–10. [Google Scholar] [CrossRef]
  37. Ratchford, A.M.; Chang, A.S.; Chi, M.M.; Sheridan, R.; Moley, K.H. Maternal diabetes adversely affects AMP-activated protein kinase activity and cellular metabolism in murine oocytes. Am. J. Physiol. Endocrinol. Metab. 2007, 293, E1198–E1206. [Google Scholar] [CrossRef]
  38. Wang, Q.; Frolova, A.I.; Purcell, S.; Adastra, K.; Schoeller, E.; Chi, M.M.; Schedl, T.; Moley, K.H. Mitochondrial dysfunction and apoptosis in cumulus cells of type I diabetic mice. PLoS ONE 2010, 5, e15901. [Google Scholar] [CrossRef]
  39. Al Zaid Siddiquee, K.; Turkson, J. STAT3 as a target for inducing apoptosis in solid and hematological tumors. Cell Res. 2008, 18, 254–267. [Google Scholar] [CrossRef]
  40. Li, O.; Li, Z.; Tang, Q.; Li, Y.; Yuan, S.; Shen, Y.; Zhang, Z.; Li, N.; Chu, K.; Lei, G. Long Stress Induced Non-Coding Transcripts 5 (LSINCT5) Promotes Hepatocellular Carcinoma Progression Through Interaction with High-Mobility Group AT-hook 2 and MiR-4516. Med. Sci. Monit. 2018, 24, 8510–8523. [Google Scholar] [CrossRef] [PubMed]
  41. Sutherland, J.M.; Keightley, R.A.; Nixon, B.; Roman, S.D.; Robker, R.L.; Russell, D.L.; McLaughlin, E.A. Suppressor of cytokine signaling 4 (SOCS4): Moderator of ovarian primordial follicle activation. J. Cell Physiol. 2012, 227, 1188–1198. [Google Scholar] [CrossRef] [PubMed]
  42. Kawamura, K.; Cheng, Y.; Suzuki, N.; Deguchi, M.; Sato, Y.; Takae, S.; Ho, C.-h.; Kawamura, N.; Tamura, M.; Hashimoto, S.; et al. Hippo signaling disruption and Akt stimulation of ovarian follicles for infertility treatment. Proc. Natl. Acad. Sci. 2013, 110, 17474–17479. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Weng, Q.; Zhao, M.; Zheng, J.; Yang, L.; Xu, Z.; Zhang, Z.; Wang, J.; Wang, J.; Yang, B.; Richard Lu, Q.; et al. STAT3 dictates β-cell apoptosis by modulating PTEN in streptozocin-induced hyperglycemia. Cell Death Differ. 2020, 27, 130–145. [Google Scholar] [CrossRef] [PubMed]
  44. Lu, C.; Lin, L.; Tan, H.; Wu, H.; Sherman, S.L.; Gao, F.; Jin, P.; Chen, D. Fragile X premutation RNA is sufficient to cause primary ovarian insufficiency in mice. Hum. Mol. Genet. 2012, 21, 5039–5047. [Google Scholar] [CrossRef]
  45. Probst, F.J.; Cooper, M.L.; Cheung, S.W.; Justice, M.J. Genotype, Phenotype, and Karyotype Correlation in the XO Mouse Model of Turner Syndrome. J. Hered. 2008, 99, 512–517. [Google Scholar] [CrossRef]
Cells 11 02797 g001 550
Figure 1. Characterization of exosomes isolated from the urine. (A) Phenotype of urinary exosomes detected via transmission electron microscopy. Scale bar: 100 nm. (B) Frequency of isolated urinary exosomes revealed via nanoparticle tracking analysis. All average particle sizes were calculated by using three independent samples. (C) The amount of total RNA from each sample was measured by using Nanodrop 2000 (ThermoFisher Scientific); patients with premature ovarian insufficiency (POI) = 11, patients with Turner syndrome = 9, control individuals = 12. (D) Expression levels of CD63, CD9, and TSG101 in urinary exosomes were detected via Western blot analysis. Tubulin was used as a negative control.
Figure 1. Characterization of exosomes isolated from the urine. (A) Phenotype of urinary exosomes detected via transmission electron microscopy. Scale bar: 100 nm. (B) Frequency of isolated urinary exosomes revealed via nanoparticle tracking analysis. All average particle sizes were calculated by using three independent samples. (C) The amount of total RNA from each sample was measured by using Nanodrop 2000 (ThermoFisher Scientific); patients with premature ovarian insufficiency (POI) = 11, patients with Turner syndrome = 9, control individuals = 12. (D) Expression levels of CD63, CD9, and TSG101 in urinary exosomes were detected via Western blot analysis. Tubulin was used as a negative control.
Cells 11 02797 g001
Cells 11 02797 g002 550
Figure 2. MicroRNA (miRNA) sequencing profile of urinary exosomes from patients with premature ovarian insufficiency (POI), Turner syndrome, and control individuals. (A) Heatmap with hierarchical clustering of 33 differentially expressed miRNAs (p < 0.05). Red and green colors represent increased or reduced expressions across samples, respectively. (B,C) Scatterplots of miRNA profiles among patients with POI versus the control (B) and patients with Turner syndrome versus the control (C). (D) Top 20 statistically significant enriched pathways in the Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses of the genes targeted by the 33 differentially expressed miRNAs. (E) Gene category analysis for the POI-related target genes of the 33 differentially expressed miRNAs.
Figure 2. MicroRNA (miRNA) sequencing profile of urinary exosomes from patients with premature ovarian insufficiency (POI), Turner syndrome, and control individuals. (A) Heatmap with hierarchical clustering of 33 differentially expressed miRNAs (p < 0.05). Red and green colors represent increased or reduced expressions across samples, respectively. (B,C) Scatterplots of miRNA profiles among patients with POI versus the control (B) and patients with Turner syndrome versus the control (C). (D) Top 20 statistically significant enriched pathways in the Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses of the genes targeted by the 33 differentially expressed miRNAs. (E) Gene category analysis for the POI-related target genes of the 33 differentially expressed miRNAs.
Cells 11 02797 g002
Cells 11 02797 g003 550
Figure 3. Relative expression of hsa-miR-4516 in the validation cohorts. (A) Relative expression of hsa-miR-4516 in the validation cohort, using qRT-PCR (patients with premature ovarian insufficiency (POI) = 15, patients with Turner syndrome = 11, and control individuals = 20). Relative miRNA expression was normalized by using UniSp6 as an internal control. (B) Relative expression of hsa-miR-4516, using real-time Taqman analysis, was normalized by ath-miR-159a as an endogenous control (patients with premature ovarian insufficiency (POI) = 6, patients with Turner syndrome = 5, and control individuals = 3). The p-value was evaluated by comparison to the control group. Asterisks (*) indicate significant differences (p < 0.05) when compared to the control group.
Figure 3. Relative expression of hsa-miR-4516 in the validation cohorts. (A) Relative expression of hsa-miR-4516 in the validation cohort, using qRT-PCR (patients with premature ovarian insufficiency (POI) = 15, patients with Turner syndrome = 11, and control individuals = 20). Relative miRNA expression was normalized by using UniSp6 as an internal control. (B) Relative expression of hsa-miR-4516, using real-time Taqman analysis, was normalized by ath-miR-159a as an endogenous control (patients with premature ovarian insufficiency (POI) = 6, patients with Turner syndrome = 5, and control individuals = 3). The p-value was evaluated by comparison to the control group. Asterisks (*) indicate significant differences (p < 0.05) when compared to the control group.
Cells 11 02797 g003
Cells 11 02797 g004 550
Figure 4. Increase in miR-4516 expression in a chemotherapy-induced premature ovarian insufficiency (POI) mouse model. (A) Representative hematoxylin-and-eosin-stained images of the ovaries from saline- and CTX+BUS-treated mice (n = 7). Scale bar = 200 μm in the upper panel and 80 μm in the lower panel. (B) Representative images of Sirius red staining in the ovaries from saline- and CTX+BUS-treated mice. Scale bar = 200 μm. (C) Representative images of the TUNEL assay in the ovaries from saline- and CTX+BUS-treated mice. Scale bar = 80 μm. TUNEL-positive cells were counted by using ImageJ software (n = 7). (D) Representative immunohistochemistry images of Ki67 expression in the ovaries from saline- and CTX+BUS-treated mice. Scale bar = 80 μm. (E) Relative expression of mmu-miR-4516 in the ovaries and uteri from the saline- and CTX+BUS-treated mice were analyzed by using qRT-PCR. Relative miRNA expression was normalized by using UniSp6 as an internal control (n = 6). (F) Western blot analysis of STAT3, Bax, and p53 expression in the ovaries of saline- and CTX+BUS-treated mice. The relative protein level in (F) was normalized to tubulin levels (mean ± SD; n = 6). Asterisks (*, **) indicate significant differences (p < 0.05, p < 0.01, respectively) when compared to saline-injected mice.
Figure 4. Increase in miR-4516 expression in a chemotherapy-induced premature ovarian insufficiency (POI) mouse model. (A) Representative hematoxylin-and-eosin-stained images of the ovaries from saline- and CTX+BUS-treated mice (n = 7). Scale bar = 200 μm in the upper panel and 80 μm in the lower panel. (B) Representative images of Sirius red staining in the ovaries from saline- and CTX+BUS-treated mice. Scale bar = 200 μm. (C) Representative images of the TUNEL assay in the ovaries from saline- and CTX+BUS-treated mice. Scale bar = 80 μm. TUNEL-positive cells were counted by using ImageJ software (n = 7). (D) Representative immunohistochemistry images of Ki67 expression in the ovaries from saline- and CTX+BUS-treated mice. Scale bar = 80 μm. (E) Relative expression of mmu-miR-4516 in the ovaries and uteri from the saline- and CTX+BUS-treated mice were analyzed by using qRT-PCR. Relative miRNA expression was normalized by using UniSp6 as an internal control (n = 6). (F) Western blot analysis of STAT3, Bax, and p53 expression in the ovaries of saline- and CTX+BUS-treated mice. The relative protein level in (F) was normalized to tubulin levels (mean ± SD; n = 6). Asterisks (*, **) indicate significant differences (p < 0.05, p < 0.01, respectively) when compared to saline-injected mice.
Cells 11 02797 g004
Table
Table 1. Sequences of the microRNA primers used in this study.
Table 1. Sequences of the microRNA primers used in this study.
microRNACat. No.Sequence
hsa-miR-16-5pYP00205702UAGCAGCACGUAAAUAUUGGCG
hsa-miR-20a-5pYP00204292UAAAGUGCUUAUAGUGCAGGUAG
hsa-miR-29a-3pYP00204698UAGCACCAUCUGAAAUCGGUUA
hsa-miR-30b-5pYP00204765UGUAAACAUCCUACACUCAGCU
hsa-miR-99b-3pYP00204064CAAGCUCGUGUCUGUGGGUCCG
hsa-miR-151a-5pYP00204007UCGAGGAGCUCACAGUCUAGU
hsa-miR-200b-3pYP00206071UAAUACUGCCUGGUAAUGAUGA
hsa-miR-423-3pYP00204488AGCUCGGUCUGAGGCCCCUCAGU
hsa-miR-941YP00204574CACCCGGCUGUGUGCACAUGUGC
hsa-miR-4516YP02112882GGGAGAAGGGUCGGGGC
hsa-miR-4492YP02100455GGGGCUGGGCGCGCGCC
hsa-miR-7847-3pYP02105945CGUGGAGGACGAGGAGGAGGC
Abbreviations (Hsa, Homo sapiens; miR, microRNA).
Table
Table 2. Baseline clinical characteristics of the enrolled participants in the three groups.
Table 2. Baseline clinical characteristics of the enrolled participants in the three groups.
Sequencing CohortValidation Cohort
POI (n = 7)Turner
(n = 7)		Control
(n = 5)		p-ValuePOI
(n = 15)		Turner
(n = 11)		Control
(n = 20)		p-Value
Age37.7 ± 8.931.0 ± 6.134.8 ± 5.00.22635.06 ± 9.9234.00 ± 7.3334.05 ± 17.110.052
Height (cm)161.5 ± 6.2148.1 ± 6.2160.9 ± 7.30.002162.39 ± 6.38148.91 ± 5.27163.81 ± 5.30<0.001
Weight (kg)56.9 ± 8.353.3 ± 6.454.8 ± 5.00.62259.00 ± 2.4655.41 ± 9.4551.75 ± 10.410.197
BMI (kg/m2)21.8 ± 2.724.4 ± 3.621.2 ± 1.50.10522.48 ± 3.2719.57 ± 4.2719.24 ± 2.720.340
SBP (mmHg)126.4 ± 3.4135.1 ± 7.9123.3 ± 3.80.003124.31 ± 4.35134.11 ± 11.12122.00 ± 4.18<0.001
DBP (mmHg)75.1 ± 8.480.9 ± 7.380.0 ± 1.70.25876.06 ± 6.1779.89 ± 8.1076.70 ± 4.040.217
AST (IU/L)20.0 ± 7.727.1 ± 16.515.0 ± 2.90.16221.00 ± 5.8927.22 ± 14.3220.85 ± 8.030.224
ALT (IU/L)15.3 ± 6.337.7 ± 40.114.7 ± 4.30.16516.44 ± 7.2335.11 ± 35.5916.50 ± 5.920.022
TC (mg/dL)189.0 ± 39.7202.1 ± 31.3162.2 ± 19.00.101185.81 ± 32.70198.56 ± 20.38185.70 ± 28.250.663
TG (mg/dL)127.9 ± 25.2222.9 ± 245.277.5 ± 22.50.219161.56 ± 124.18141.80 ± 47.27110.30 ± 43.740.296
HDL-C (mg/dL)73.9 ± 18.280.3 ± 17.556.8 ± 14.20.063109.63 ± 41.12105.04 ± 22.1098.95 ± 17.560.013
LDL-C (mg/dL)98.1 ± 31.691.9 ± 9.396.9 ± 11.70.839109.63 ± 41.12105.04 ± 22.1098.95 ± 17.560.534
Prolactin (ng/mL)8.5 ± 2.87.9 ± 2.29.1 ± 1.40.6558.10 ± 2.867.76 ± 2.178.25 ± 2.470.854
TSH (mIU/mL)1.2 ± 0.62.6 ± 1.11.9 ± 0.40.0101.50 ± 1.072.86 ± 1.001.70 ± 1.25<0.001
FSH (mIU/mL)114.8 ± 35.152.8 ± 15.57.1 ± 0.7<0.00192.00 ± 32.9650.46 ± 26.747.2 ± 1.72<0.001
Estradiol (pg/mL)24.6 ± 16.715.7 ± 4.539.6 ± 4.20.00317.95 ± 12.5819.72 ± 8.3781.27 ± 5.27<0.001
Abbreviations: POI, premature ovarian insufficiency; Turner, Turner syndrome; control, control individuals; BMI, body mass index; SBP, systolic blood pressure; DBP, diastolic blood pressure; AST, aspartate aminotransferase; ALT, alanine aminotransferase; TC, total cholesterol; TG, triglyceride; HDL-C, high-density lipoprotein cholesterol; LDL-C, low-density lipoprotein cholesterol; TSH, thyroid-stimulating hormone; FSH, follicle-stimulating hormone).
Table
Table 3. Menstrual history of the patients and serum chromosomal analysis of both premature ovarian insufficiency and Turner syndrome groups in the sequencing cohort.
Table 3. Menstrual history of the patients and serum chromosomal analysis of both premature ovarian insufficiency and Turner syndrome groups in the sequencing cohort.
Menarche (Age)Duration of
Amenorrhea (Months)		KaryotypeFMR
Gene Mutation
POI
1151846,XXNormal
2132447,XXX [39]/46,XX[1]Normal
3152446,XXNormal
4163646,XXNormal
512346,XXNormal
615546,XXNormal
7143046,XXNormal
Turner Syndrome
1no menarcheNA46,X,i(X)(q10),22pstk+[16]/45,X,22pstk+[14]Normal
2no menarcheNA46,X,i(X)(q10)Normal
3no menarcheNA45,XNormal
4no menarcheNA46,X,i(X)(q10)Normal
5no menarcheNA46,X,i(X)(q10)[24]/45,X[6]Normal
6no menarcheNA45,XNormal
7no menarcheNA45,XNormal
Abbreviations: FMR, fragile X mental retardation 1; NA, not available; POI, premature ovarian insufficiency.
Table
Table 4. Menstrual history of the patients and serum chromosomal analysis of both premature ovarian insufficiency and Turner syndrome groups in the validation cohort.
Table 4. Menstrual history of the patients and serum chromosomal analysis of both premature ovarian insufficiency and Turner syndrome groups in the validation cohort.
Menarche (Age)Duration of Amenorrhea (Months)KaryotypeFMR
Gene Mutation
POI
1no menarcheNA46,XXNormal
2151846,XXNormal
315161046,XXNormal
41414646,XXNormal
512167046,XXNormal
61215646,XX,16qh+Normal
71112046,XXNormal
8132447,XXX[39]/46,XX[1]Normal
91412046,XX,1qh+Normal
10152446,XXNormal
11131246,XXNormal
12142446,XXNormal
13123646,XXNormal
14no menarcheNA46,XXNormal
15142246,XXNormal
Turner Syndrome
1no menarcheNA45,XNormal
2no menarcheNA46,X,i(X)(q10)Normal
3no menarcheNA46,X,i(X)(q10)[24]/45,X[6]Normal
4no menarcheNA45,XNormal
5no menarcheNA45,X[21]/46,X,idic(X)(q25)[9]Normal
6no menarcheNA46,X,i(X)(q10)Normal
7no menarcheNA45,XNormal
8no menarcheNA45,X[17]/46,X,i(X)(q10)[13]Normal
9no menarcheNA46,X,i(X)(q10),22pstk+[16]/45,X,22pstk+[14]Normal
10no menarcheNA45,X/46,XXNormal
11no menarcheNA46,X,idic(X)(q22.1)[28]/45,X[22]Normal
Abbreviations: FMR, fragile X mental retardation 1; NA, not available; POI, premature ovarian insufficiency.
Table
Table 5. Differentially expressed microRNAs in the sequencing cohort.
Table 5. Differentially expressed microRNAs in the sequencing cohort.
IDGene SymbolFold Changep-ValueNormalized Data (log2)
P/CT/CT/PP_T/CP/CT/CT/PP_T/CP/CT/CT/PP_T/C
1hsa-miR-451635.5497.2540.20421.4010.3280.0440.3340.4065.41510.5678.2749.835
2hsa-miR-423-3p0.0281.80364.0760.9160.0720.4740.0470.92210.1504.99911.00010.023
3hsa-miR-1273g-5p0.3271.5994.8970.9630.0720.4230.0330.9491.6240.0092.3011.569
4hsa-miR-193b-5p0.0802.20827.5361.1440.2620.3790.0340.8907.9354.2959.0788.129
5hsa-miR-455-5p0.0742.64035.4411.3570.2530.2820.0270.7613.7570.0105.1584.198
6hsa-miR-3184-5p0.0703.18245.4081.6260.1550.2360.0400.6589.2325.39710.9029.933
7hsa-miR-574-5p0.37512.51833.4036.4460.3300.1170.0500.3282.9751.5596.6205.663
8hsa-miR1273d0.6795.5828.2253.1300.3560.0290.0060.2240.7640.2053.2452.410
9hsa-miR-92b-5p0.1856.27433.9013.2290.1580.1500.0490.4242.6670.2335.3174.359
10hsa-miR-45080.3615.37214.8602.8670.4220.1480.0480.4195.0413.5737.4666.560
11hsa-let-7g-5p0.8370.2060.2460.5210.6560.0050.0320.12110.43510.1788.1579.496
12hsa-miR-125b-5p0.8540.2890.3390.5710.8640.0420.4150.48010.50710.2798.7179.700
13hsa-let-7c-5p0.1420.0840.5960.1130.0250.0170.4350.00110.9668.1487.4007.822
14hsa-miR-26a-5p0.4030.2450.6090.3240.0810.0260.2840.00711.85410.5429.82610.228
15hsa-miR-942-5p0.2510.2270.9050.2390.1620.1500.3070.0382.1400.1440.0000.074
16hsa-miR-18a-5p0.0170.0150.9030.0160.0900.0900.3000.0146.0320.1470.0000.075
17hsa-miR-18b-5p0.1510.1500.9940.1510.0790.0790.2670.0112.7340.0090.0000.004
18hsa-miR-1273a0.1060.1050.9940.1050.1770.1760.2690.0483.2530.0090.0000.005
19hsa-miR-125a-5p0.1050.1091.0320.1070.1000.0950.9740.01712.1968.9468.9928.969
20hsa-miR-44520.1990.2031.0200.2010.1020.1040.4630.0182.3370.0080.0370.023
21hsa-miR-183-5p0.1100.1361.2340.1230.1160.1250.8540.0276.7923.6063.9093.766
22hsa-miR-27b-3p0.2850.3041.0650.2940.0760.0770.9210.01710.6598.8488.9398.894
23hsa-miR-92b-3p0.0060.0081.2740.0070.0830.0830.5370.0129.0621.6722.0211.857
24hsa-miR-10a-5p0.0770.1531.9850.1150.1420.1780.3980.04113.89510.19811.18710.776
25hsa-miR-143-3p0.1140.1591.3920.1370.1070.1240.6020.02310.0436.9127.3897.170
26hsa-miR-16-5p0.0070.09913.5750.0530.0380.0580.0950.0047.1080.0103.7732.876
27hsa-miR-200a-3p0.0250.1656.6570.0950.1010.1640.2300.02910.3435.0127.7476.949
28hsa-let-7f-5p0.3480.5541.5930.4510.0260.2560.4470.05610.7299.2069.8789.581
29hsa-miR1260a0.0460.1623.5400.1040.1080.1560.0730.0285.3310.8822.7062.065
30hsa-miR-30d-3p0.2050.8774.2740.5410.2500.8710.0380.3972.2940.0092.1051.408
31hsa-miR-29a-3p0.0010.141248.6950.0710.1360.1970.0930.04310.9940.2058.1637.169
32hsa-miR-30b-5p0.0100.24024.0830.1250.0640.1720.1990.0259.4342.7887.3786.437
33hsa-miR-151a-5p0.1050.5244.9700.3150.0190.2480.0450.0274.1030.8573.1702.434
Fold change is presented as miRNA expression in patients with POI compared to the control individuals. Abbreviations: P, premature ovarian insufficiency; T, Turner syndrome; C, control individuals; P_T, average expression level of premature ovarian insufficiency and Turner syndrome; Hsa, Homo sapiens; miR, microRNA.
Table
Table 6. Comparison of the relative expression patterns of selected microRNAs in patients with premature ovarian insufficiency.
Table 6. Comparison of the relative expression patterns of selected microRNAs in patients with premature ovarian insufficiency.
miRNAmiRNA SequencingValidation Results Obtained Using qRT-PCR
Up/DownFold Changep-ValueUp/DownFold Changep-Value
hsa-miR-16-5pDown0.0070.038Down0.2030.170
hsa-miR-20a-5pDown0.0090.101Up2.4100.361
hsa-miR-29a-3pDown0.0010.136Down0.1860.066
hsa-miR-30b-5pDown0.0100.064Down0.2180.026
hsa-miR-99b-3pDown0.0600.254Down0.5050.185
hsa-miR-151a-5pDown0.1050.019Down0.8130.751
hsa-miR-200b-3pDown0.0020.096Down0.5300.439
hsa-miR-423-3pDown0.0280.072Down0.0890.153
hsa-miR-941Down0.0640.254Down0.2830.065
hsa-miR-4492Up23.0910.344Down0.2760.094
hsa-miR-4516Up35.5490.328Up4.9710.008
hsa-miR-7847-3pUp10.1750.085Up1.4490.568
These patterns were obtained by using miRNA sequencing and validated by using qRT-PCR; Fold change is presented as miRNA expression in patients with POI compared to the control individuals. Abbreviations: Hsa, Homo sapiens; miR, microRNA; qRT-PCR, quantitative real-time polymerase chain reaction.
Table
Table 7. Comparison of the relative expression patterns of selected microRNAs in patients with Turner syndrome.
Table 7. Comparison of the relative expression patterns of selected microRNAs in patients with Turner syndrome.
miRNAmiRNA SequencingValidation Results by qRT-PCR
Up/DownFold Changep-ValueUp/DownFold Changep-Value
hsa-miR-16-5pDown0.0990.058Up1.6730.811
hsa-miR-20a-5pUp1.0140.987Up1.8430.447
hsa-miR-29a-3pDown0.1410.197Up9.7880.021
hsa-miR-30b-5pDown0.2400.284Up5.3800.119
hsa-miR-99b-3pUp14.0810.171Down0.1750.011
hsa-miR-151a-5pDown0.5240.248Up1.0080.991
hsa-miR-200b-3pDown0.6920.671Down0.4390.296
hsa-miR-423-3pUp1.8030.474Down0.3450.363
hsa-miR-941Up5.1380.235Down0.2070.064
hsa-miR-4492Up13.2500.058Down0.2860.107
hsa-miR-4516Up7.2540.044Up2.6630.022
hsa-miR-7847-3pUp12.2520.212Down0.2860.026
These patterns were obtained by using miRNA sequencing and validated by using qRT-PCR; Fold change is presented as miRNA expression in patients with Turner syndrome compared to the control individuals. Abbreviations: Hsa, Homo sapiens; miR, microRNA; qRT-PCR, quantitative real-time polymerase chain reaction.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Umair, Z.; Baek, M.-O.; Song, J.; An, S.; Chon, S.J.; Yoon, M.-S. MicroRNA-4516 in Urinary Exosomes as a Biomarker of Premature Ovarian Insufficiency. Cells 2022, 11, 2797. https://doi.org/10.3390/cells11182797

AMA Style

Umair Z, Baek M-O, Song J, An S, Chon SJ, Yoon M-S. MicroRNA-4516 in Urinary Exosomes as a Biomarker of Premature Ovarian Insufficiency. Cells. 2022; 11(18):2797. https://doi.org/10.3390/cells11182797

Chicago/Turabian Style

Umair, Zobia, Mi-Ock Baek, Jisue Song, Seona An, Seung Joo Chon, and Mee-Sup Yoon. 2022. "MicroRNA-4516 in Urinary Exosomes as a Biomarker of Premature Ovarian Insufficiency" Cells 11, no. 18: 2797. https://doi.org/10.3390/cells11182797

APA Style

Umair, Z., Baek, M. -O., Song, J., An, S., Chon, S. J., & Yoon, M. -S. (2022). MicroRNA-4516 in Urinary Exosomes as a Biomarker of Premature Ovarian Insufficiency. Cells, 11(18), 2797. https://doi.org/10.3390/cells11182797

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