Statement of ethics
The study was conducted in accordance with the Declaration of Helsinki and was approved by the Medical Ethics Committee of Ganzhou Hospital of Traditional Chinese Medicine (approval no. GZSZYYKYLL20240092, Date: 2024-12-20).
All materials, reagents, assay kits, and equipment used in this protocol are listed in the Table of Materials.
Study design
This retrospective case-control study included 62 patients with PCOS diagnosed at the gynecology outpatient clinic of Ganzhou Hospital of Traditional Chinese Medicine between June 2024 and August 2025 and 71 healthy women of reproductive age examined at the hospital’s health management center during the same period. Controls were selected to achieve group-level comparability in age; BMI was not used as a matching variable and was evaluated as a baseline characteristic.
Study subjects
Inclusion criteria: The inclusion criteria for the PCOS group were as follows: (1) age >18 years; (2) diagnosis of PCOS according to the modified Rotterdam criteria endorsed by the 2023 International Evidence-based Guideline, requiring at least two of the following three features after exclusion of related disorders: oligo-ovulation or anovulation, defined as menstrual cycles longer than 35 days or fewer than eight cycles per year; clinical hyperandrogenism and/or biochemical hyperandrogenism, with biochemical hyperandrogenism defined as a serum total T concentration above the assay-specific upper reference limit; and polycystic ovarian morphology, defined as at least 20 follicles in one ovary or an ovarian volume of at least 10 mL when complete follicle counting was not feasible14; (3) complete electronic medical records; (4) no hormone replacement therapy, hypoglycemic agents, or immunosuppressants during the preceding 3 months; and (5) availability of a fasting serum sample collected at enrollment and stored in the institutional biobank.
The inclusion criteria for the control group were as follows: (1) attendance for a health examination during the same period; (2) regular menstrual cycles; (3) normal ovarian morphology on transvaginal ultrasonography; (4) no history of endocrine disorders, gynecological disease, infertility, or clinically significant hepatic or renal dysfunction; and (5) availability of a fasting serum sample collected on the day of examination and stored in the institutional biobank.
The exclusion criteria were as follows: (1) endometriosis, uterine fibroids, ovarian tumors, pelvic inflammatory disease, or other gynecological disorders; (2) thyroid dysfunction, Cushing syndrome, hyperprolactinemia, congenital adrenal hyperplasia, or other endocrine disorder that could account for ovulatory dysfunction or hyperandrogenism; (3) severe cardiac, hepatic, renal, or pulmonary dysfunction, malignant tumors, autoimmune diseases, or hematologic disorders; (4) pregnancy, delivery, or abortion within the preceding 6 months; (5) use of glucocorticoids, oral contraceptives, or other medications affecting glucose or lipid metabolism during the preceding 3 months; (6) psychiatric disorders or cognitive impairment; (7) missing key clinical information; and (8) severe serum hemolysis, defined as a hemoglobin concentration >2.0 g/L, lipemia, jaundice, more than one freeze–thaw cycle, or interrupted sample storage. Thyroid dysfunction, hyperprolactinemia, and non-classic congenital adrenal hyperplasia were excluded using thyroid-stimulating hormone, serum prolactin, and 17-hydroxyprogesterone measurements, respectively; follicle-stimulating hormone was reviewed where appropriate, while Cushing syndrome or an adrenal tumor was investigated when clinically indicated. The same pre-analytical procedure was applied to all clinical samples, and no sample included in the primary analysis underwent more than one freeze–thaw cycle.
Serum sample collection and storage
Fasting venous blood (10 mL) was collected in the morning into anticoagulant-free tubes, allowed to clot at room temperature for 30 min, and centrifuged at 1,505 x g for 15 min at 4 °C. Serum was divided into 200 µL aliquots and stored at −80 °C until analysis. The same pre-analytical procedure was applied to all clinical samples, and no sample included in the primary analysis underwent more than one freeze–thaw cycle.
Preparation of samples for methodological validation
Ten pooled serum samples from healthy volunteers were used to evaluate the effects of centrifugation speed (1,505, 3,000, and 4,000 x g for 15 min), freeze–thaw frequency (0, 1, 3, and 5 cycles), storage temperature (4 °C, −20 °C, and −80 °C for 7 days), storage duration (1, 3, 6, and 12 months at −80 °C), and hemolysis (<0.1, 0.1–0.5, and 0.5–2.0 g/L hemoglobin).
Laboratory testing
All laboratory assays were performed by the same trained technician under standardized environmental conditions and according to prespecified operating procedures. Serum miRNAs were quantified by qRT-PCR; ADPN, LEP, and RBP4 were measured by sandwich ELISA; T and IL-6 were measured by chemiluminescence-based immunoassays; and glucose-metabolic variables were obtained using routine automated laboratory methods.
Serum miRNA quantification by qRT-PCR
Total RNA was extracted from 200 µL of serum using TRIzol reagent. Spectrophotometric RNA concentration and A260/A280 purity data were not available for the archived serum RNA preparations. qRT-PCR assay acceptability was therefore assessed by the presence of a single specific melting-curve peak and concordance among three technical replicates. After chloroform phase separation, isopropanol precipitation, and two washes with 75% ethanol, the RNA pellet was dissolved in enzyme-free water. The 20 µL reaction system contained 4 µL of 5x buffer, 1 µL of reverse transcriptase, 1 µL of 10 µM stem-loop primer, and 5 µL of RNA template. The reaction conditions were 16 °C for 30 min, 42 °C for 30 min, and 85 °C for 5 s. The reverse-transcribed products were stored at −20 °C. Quantitative PCR was carried out using the TaKaRa SYBR Premix Ex Taq II kit. Each 20 µL reaction contained 10 µL of 2× SYBR mixture, 0.8 µL each of forward and reverse primers (10 µM), and 2 µL of cDNA template. Amplification was performed on an Applied Biosystems 7500 real-time PCR system under the following conditions: pre-denaturation at 95 °C for 30 s; 40 cycles of 95 °C for 5 s and 60 °C for 30 s; melting curve analysis was then performed to confirm amplification specificity. U6 snRNA served as the internal reference gene, and relative expression levels of miR-146a, miR-642a, and miR-2861 were calculated using the 2-ΔΔCt method. Each sample was tested in triplicate. The sequences of all primers used in this study are listed in Table 1.
Serum protein marker detection
Serum ADPN, LEP, and RBP4 concentrations were measured by double-antibody sandwich ELISA using the Human Total ADPN/Acrp30 Quantikine ELISA Kit (analytical sensitivity, 0.891 ng/mL; assay range, 3.9–250 ng/mL), Human Leptin Quantikine ELISA Kit (analytical sensitivity, 7.8 pg/mL; assay range, 15.6–1,000 pg/mL), and Human RBP4 Quantikine ELISA Kit (analytical sensitivity, 0.628 ng/mL; assay range, 1.56–100 ng/mL), respectively. Each assay was specific for the corresponding human analyte, and serum samples were diluted according to the manufacturer’s instructions. Absorbance was measured at 450 nm, and concentrations were calculated using four-parameter logistic standard curves. The intra-assay or inter-assay CVs were 2.5% or 3.2% for ADPN, 2.7% or 3.4% for LEP, and 2.6% or 3.3% for RBP4.
T measurement
Serum T was measured by electrochemiluminescence immunoassay using a Roche Cobas e601 automated immunoassay analyzer.
Glucose-metabolic and inflammatory measurements
Fasting plasma glucose, fasting insulin, and HbA1c values were obtained from same-day routine clinical laboratory records rather than from the archived serum aliquots used for the miRNA and metabolic protein assays. Fasting plasma glucose was measured by the glucose oxidase method, fasting insulin by immunoturbidimetry, and HbA1c by high-performance liquid chromatography. HOMA-IR was calculated as [fasting plasma glucose (mmol/L) × fasting insulin (mIU/L)]/22.5. Serum IL-6 was measured by chemiluminescence immunoassay.
Methodological validation
The alternative-method comparisons were performed to determine whether the primary measurements were robust to assay platform and were not used to generate the principal between-group results. Expression levels of the three miRNAs and five biochemical markers were measured under the different processing conditions described above, and the coefficient of variation for each indicator was calculated as CV = (standard deviation/mean) x 100%. A CV <15% was considered acceptable for evaluating the stability of assay results under different sample processing and storage conditions. Within-day precision was assessed by repeated measurement of the same pooled serum sample 10 times on a single day, whereas between-day precision was determined from measurements performed once daily over five consecutive days. For each indicator, the mean, standard deviation, and corresponding CV were calculated. For analytical agreement testing, five biochemical analytes—ADPN, LEP, RBP4, IL-6, and T—were measured in the validation samples using the primary assay and a Roche chemiluminescence-based immunoassay platform. Bland–Altman analysis was used to calculate the 95% limits of agreement. The primary qRT-PCR assay was used for all clinical miRNA measurements, whereas the miScript SYBR Green PCR Kit was used only as an orthogonal method in the analytical agreement experiment; Pearson correlation coefficients were calculated between the two qRT-PCR systems.
Statistical analysis
Data analysis was performed using statistical software. All continuous variables were first tested for normality with the Shapiro-Wilk test. Data with a normal distribution are expressed as mean ± standard deviation. Comparisons between two groups were made with the independent-samples t test, while comparisons among multiple groups were performed using one-way analysis of variance, followed by the LSD-t test for pairwise comparisons. Data not conforming to a normal distribution are presented as median (interquartile range) [M (P25, P75)]. Comparisons between two groups were then performed with the Mann-Whitney U test, and comparisons among multiple groups with the Kruskal-Wallis H test. Categorical data were presented as frequencies and percentages. We used the χ2 test or Fisher’s exact test to compare the groups. Correlations between variables were assessed using Pearson correlation analysis or Spearman rank correlation analysis according to data distribution. As an exploratory analysis, Pearson correlations between each miRNA and the six metabolic markers were repeated within non-obese (BMI <28 kg/m2) and obese (BMI ≥28 kg/m2) PCOS subgroups. Correlation coefficients were compared using Fisher’s r-to-z transformation, and the resulting between-subgroup p values were adjusted using the Benjamini–Hochberg false-discovery-rate procedure. We evaluated diagnostic performance using receiver operating characteristic (ROC) curves and the area under the curve (AUC). In addition, we analyzed combined detection using a logistic regression model. All tests were two-sided, and p < 0.05 was considered statistically significant.