1. Introduction
Minipuberty of infancy is a transient sex-specific phase of postnatal activation of the reproductive axis. Concentrations of both pituitary gonadotropins, follicle-stimulating hormone (FSH) and luteinizing hormone (LH), start to rise in the first few days after birth in response to disappearance of the suppressive effect of placental hormones, resulting in stimulation of the gonads [
1,
2,
3,
4]. Peak concentrations of FSH in infant girls in this period of life are higher than those of LH, while the subsequent increase in estradiol levels reaches values observed during Tanner stage 4 [
3,
4]. FSH and estradiol levels remain elevated until the ages of 3–4 and 2 years, respectively, while LH decreases to prepubertal levels at around the age of 6 months [
3]. After completing minipuberty, the hypothalamic–pituitary–gonadal axis activity remains in a state of relative quiescence until reactivation at puberty [
1]. The biological significance of minipuberty in girls is still insufficiently understood. It seems to play an important role in the postnatal development of the uterus and breasts [
5]. Hormonal changes associated with minipuberty may also determine sex differences in growth velocity and the pattern of body adipose tissue distribution [
2]. Lastly, minipuberty probably plays a role in forming brain plasticity and in the development of cognitive abilities and emotional competences [
4].
No previous study has investigated a relationship between exposure to subnormal levels of thyroid hormones and the course of minipuberty in girls. However, the results of a nationwide Danish cohort study showed no differences in pubertal development between daughters born to women with thyroid disorders (including hypothyroidism) and female descendants of healthy women [
6]. Unfortunately, the data were not collected by the researchers (the study was register-based), which might have affected their quality. Moreover, the study design made it difficult to differentiate between prevalent and incident cases [
7]. Consequently, some women might have been misdiagnosed with hypothyroidism, and the population of hypothyroid women probably differed in substitution effectiveness. In turn, animal studies assessing genital maturation of female offspring of dams with gestational hypothyroidism provided inconsistent results. Some authors observed differences in the volume, size and weight of ovaries, as well as in the number of secondary follicles between the offspring of females with propylthiouracil-induced gestational hypothyroidism and those of healthy females [
8,
9]. Furthermore, the diameter of Graafian follicles in young female rodents was greater if dams were exposed in utero to propylthiouracil [
10]. However, timing of vaginal opening, a widely used external sign of puberty, was not affected by prenatal expose to propylthiouracil-induced maternal hypothyroidism [
11]. Lastly, female offspring of animals receiving propylthiouracil and melatonin during pregnancy reached puberty, became anestrus, and returned to cyclicity at similar times to contemporary ewe lambs [
12].
Some clinical observations may indirectly suggest an association between gestational hypothyroidism and the course of female minipuberty. Thyroid hypofunction in pregnant women is associated with an increased risk of premature delivery and abnormal (usually high) birth weight [
13,
14]. In turn, premature female infants are characterized by higher gonadotropin and estradiol concentrations than full-term ones [
15,
16], which is suggestive of increased activity of the hypothalamic–pituitary–gonadal axis at the hypothalamic or pituitary level. Moreover, although there are no data for children with high birth weight, infant girls born small for gestational age had higher FSH and sex hormone levels than those born appropriate for gestational age [
17]. Lastly, decreased growth velocity and impaired neurocognitive development are typical features of untreated or undertreated congenital hypothyroidism [
18].
The lack of dedicated human studies and the inconsistent results of animal ones encouraged us to investigate whether maternal hypothyroidism during pregnancy impacts the course of minipuberty in their female offspring. The current study, including only girls, is a parallel study to our recently published work that showed that gestational thyroid hypofunction affected the activity of the reproductive axis and the development of genital organs in infant boys [
19].
2. Results
Four children (two from group 1, one from group 2, and one from group 3) dropped out because they did not participate in the required number of visits. Two children (assigned to group 2 and 3) were withdrawn because they developed asthma and gastroesophageal reflux, requiring chronic pharmacotherapy. One child (from group 1) was withdrawn because 6 months after birth, her mother was diagnosed with Rathke cleft cyst, which was probably present earlier and might have affected the hypothalamic–pituitary–thyroid axis activity during pregnancy. The last child (assigned to group 2) prematurely terminated the study because of going abroad. The results of the remaining 88 infant girls (91.6%) were statistically analyzed.
There were no differences between the study groups in gestational age of delivery, birth order, anthropometric measures (length, weight, weight for length percentile and head circumference), whole-blood thyrotropin concentration between days 3 and 5 after birth, and in the frequency of breastfeeding (
Table 1).
The study groups were comparable with respect to maternal age, education, employment, smoking, body mass index (BMI), and blood pressure (
Table 2). Groups 1 and 2 did not differ statistically in conditions underlying thyroid hypofunction. Thyroid autoimmunity was reported in 15 women (52%) in group 1 and in 14 women (48%) in group 2, partial or total thyroid removal in 9 women (31%) in group 1 and in 8 (28%) women in group 2, and thyroid aplasia/hypoplasia in 3 women (10%) in group 1 and in 4 women (14%) in group 2, whereas genetic defects in the synthesis of thyroid hormones in 2 women (7%) in group 1 and in 3 women (10%) in group 2.
Estradiol was detectable in saliva during the first six months of life in group 1 and during the first 10 months in groups 2 and 3. Stable levels were observed from month 1 to month 4 in group 1 and from month 1 to month 5 in the remaining groups, and decreased thereafter. From month 1 to month 10, estradiol concentrations were lower in group 1 than in groups 2 and 3 (
Table 3).
The remaining steroid hormones were detectable in saliva from month 1 to month 5 (testosterone and androstenedione) and from month 1 to month 12 (dehydroepiandrosterone sulfate [DHEA-S], progesterone, and 17-hydroxyprogesterone). Their concentrations and the detection periods were similar in all study groups and remained stable over the entire detection period (
Table 4,
Table 5,
Table 6,
Table 7 and
Table 8).
FSH was detectable in urine for the first 10 months of life in group 1 and for the first 12 months of life in groups 2 and 3. Stable hormone concentrations were observed from month 1 to month 8 in group 1 and from month 1 to month 10 in the remaining groups. During the first 12 months of life, FSH concentrations were lower in group 1 than groups 2 and 3 (
Table 9).
Urinary LH levels were detectable from month 1 to month 5 in group 1 and from month 1 to month 8 in groups 2 and 3. In all study groups, LH concentrations were stable between month 1 and month 5 and decreased thereafter. During the whole period of detection, urinary LH concentrations were lower in group 1 than in the remaining two groups (
Table 10).
In group 1, ovarian volume remained stable for the first 6 months and decreased thereafter. In groups 2 and 3, ovarian volume was greater between month 2 and month 10 than in months 1, 12, 15, and 18. Between months 2 and 18, ovarian volume was smaller in group 1 than in the remaining groups. There were no differences in ovarian volume between groups 2 and 3 (
Table 11).
There was a decrease in uterine length between month 1 and month 5 in group 1 and between month 1 and month 2 in the remaining groups, with no changes in uterine length thereafter. From month 5 to month 18, uterine length was smaller in group 1 than in the remaining groups. There were no differences in uterine length between groups 2 and 3 (
Table 12).
In all study groups, there was a decrease in breast diameter between month 1 and month 2, with no changes in this diameter thereafter. From month 5 to month 18, breast size was smaller in group 1 than in the remaining groups. There were no differences in breast diameter between groups 2 and 3 (
Table 13).
There were no differences in concentrations of estradiol and gonadotropins between daughters of women with hypothyroidism untreated or inadequately treated during pregnancy resulting from autoimmune thyroid disease (n = 15) and of non-autoimmune origin (
Table 14).
Over the entire detection period, salivary estradiol concentrations positively correlated with urinary LH concentrations (group 1: r values between 0.46 [p = 0.0001] and 0.62 [p < 0.0001] depending on time point [dotp]; group 2: r values between 0.50 [p < 0.0001] and 0.67 [p < 0.0001] dotp; group 3: r values between 0.51 [p < 0.0001] and 0.69 [p < 0.0001] dotp). There were also positive correlations between ovarian volume and concentrations of LH (group 1: r values between 0.29 [p = 0.0424] and 0.43 [p = 0.0008] dotp; group 2: r values between 0.30 [p = 0.0402] and 0.40 [p = 0.0011] dotp; group 3: r values between 0.32 [p = 0.0346] and 0.46 [p = 0.0002] dotp) and FH (group 1: r values between 0.27 [p = 0.0498] and 0.39 [p = 0.0015] dotp; group 2: r values between 0.28 [p = 0.0464] and 0.38 [p = 0.0014] dotp; group 3: r values between 0.28 [p = 0.0471] and 0.40 [p = 0.0012] dotp). Lastly, there were positive correlations between uterine length and estradiol concentration (group 1: r values between 0.28 [p = 0.0424] and 0.41 [p = 0.0008] dotp; group 2: r values between 0.36 [p = 0.0242] and 0.46 [p = 0.0003] dotp; group 3: r values between 0.37 [p = 0.0024] and 0.49 [p = 0.0001] dotp), between breast diameter and estradiol concentration (group 1: r values between 0.26 [p = 0.0495] and 0.38 [p = 0.0228] dotp; group 2: r values between 0.38 [p = 0.0023] and 0.48 [p = 0.0001] dotp; group 3: r values between 0.38 [p = 0.0027] and 0.50 [p < 0.0001] dotp), and in two groups between breast diameter and FSH concentration (group 2: r values between 0.31 [p = 0.0286] and 0.44 [p = 0.0006] dotp; group 3: r values between 0.29 [p = 0.0426] and 0.40 [p = 0.0014] dotp). The remaining correlations were insignificant.
3. Discussion
In the current study, we observed remarkable differences in the detection period for salivary steroid levels in the infant girls. The detection period was longest (12 months) for DHEA-S, progesterone, and 17-hydroxyprogesterone. This finding is in line with the results of our recent studies including boys, in whom these steroids were detectable at all time points, but the study lasted only 12 months [
19,
20]. In turn, detectable concentrations of the remaining steroids in daughters of healthy mothers were found only during the first ten (estradiol) or only during the first five (testosterone and androstenedione) months of life. Interestingly, estradiol levels were higher, while testosterone much lower than salivary concentrations of these hormones in peer boys [
20]. These findings are probably clinically relevant. Unconjugated and protein-unbound steroid hormones easily cross blood capillaries and epithelial cells of salivary acini, and therefore their salivary content reflects the biologically active plasma pool of these hormones [
21,
22]. Among the assessed steroids, only DHEA-S is lipid-insoluble, owing to esterification [
21]. Although its concentration in saliva is over 100 times lower than that of unbound DHEA-S in plasma [
21], relatively high levels in the participants of the current study may be explained by much higher production in comparison with other steroid hormones, and possibly also by paracellular transport through tight junctions. Based on measurements in saliva, some conclusions may be drawn. Firstly, the salivary steroid profile in infancy changes with age, is specific for each steroid hormone, and differs between girls and boys. Thus, reference values for these hormones in the first 12 months of life should consider age and sex. Secondly, undetectable amounts of salivary estradiol, testosterone, and androstenedione in older infants do not have to mean their complete absence. It is probable that these steroid hormones would be detectable in assays with lower cut-off values than in those used in the current study. Lastly, abundant salivation in the first months of life, as well as easy and painless collection, indicates that saliva may be a material of choice for assessment of endogenous steroid production in infancy.
Although some peptides, particularly insulin, are actively transferred from the place of production into the salivary glands and secreted by exocytosis, this phenomenon has not been reported in the case of gonadotropins [
23,
24]. Their sequestration in the salivary glands is delayed in comparison with steroid hormones, and this lag time period may have an impact on the detection of relationships between gonadotropins and other assessed biomarkers [
21]. Thus, gonadotropin levels were assessed in urine, the collection of which is painless and does not generate stress associated with multiple blood collections. Moreover, compared to serum, urine measurements are less subject to variability resulting from pulsatile gonadotropin (mainly FSH) secretion, which is characteristic for normal female minipuberty [
2]. In line with previous observations [
25], we have observed differences in the detection period of both gonadotropins, which was longer for FSH than LH (12 vs. 8 months in case of the offspring of healthy women and 10 vs. 5 months if maternal hypothyroidism was uncontrolled or poorly controlled). The presence of correlations between estradiol in morning saliva and LH in morning urine in all study groups supports the accuracy of the applied methodology.
The major finding of the current study is, however, that untreated or inadequately substituted maternal hypothyroidism affects the course of female minipuberty. Daughters of mothers with thyrotropin levels repeatedly exceeding the upper reference limit for pregnant women recommended by the American Thyroid Association [
26] had lower concentrations of FSH, LH, and estradiol compared with daughters of healthy women. Moreover, gonadotropins and estradiol were detectable for a shorter period of time in the descendants of women with impaired than with normal thyroid function during pregnancy. Interestingly, most mothers with untreated or undertreated thyroid hypofunction were characterized by only mildly increased thyrotropin levels and met the criteria of subclinical hypothyroidism. Thus, it seems that, during pregnancy, all, even mild, cases of hypothyroidism should be detected and adequately managed. Interestingly, impaired activity of the hypothalamic–pituitary–gonadal axis cannot be explained by early postnatal changes in thyroid function because the study groups did not differ in whole-blood thyrotropin levels obtained between days 3 and 5 of life. The obtained results are generally in line with our previous observations concerning sons of hypothyroid mothers [
19]. However, infant girls may be affected to an even greater degree than infant boys because in addition to differences in hormone concentrations, they were characterized by earlier completion of minipuberty (
Table 15). During the first half of pregnancy, the fetus cannot yet produce sufficient amounts of thyroxine and triiodothyronine and is reliant on maternal thyroid hormones [
26]. Thus, decreased thyroid hormone availability may disturb the initial phase of reproductive axis activation in the early and mid-pregnancy, preceding the second phase (minipuberty) and possible modulating its course [
4,
27].
In line with previous observations concerning serum and urine [
16,
28,
29], the current study has shown detectable saliva estradiol concentrations that are higher and present for a longer period of time in daughters of women with normal thyroid function than in the offspring of women with thyroid hypofunction. This finding contrasted with no between-group differences in salivary androgens. However, it is in line with detection of analogous differences in ovarian volume, correlating with gonadotropin concentrations, which suggests impaired activity of the hypothalamic–pituitary–gonadal axis in daughters of untreated or inadequately treated women at the level of ovaries. Interestingly, there were differences between girls and boys in the pattern of sex hormone production, not only in case of descendants of healthy mothers but also if they were born to mothers with thyroid hypofunction uncontrolled or poorly controlled during pregnancy. Irrespectively of the study group, salivary estrogen concentrations were higher and its detection period was longer in infant girls than in infant boys, which contrasted with more pronounced changes in testosterone in boys [
19]. This observation is in agreement with autoptic studies demonstrating larger content of estrogens in the ovaries than in the testes of children suddenly dying during the first two years of life [
30]. Differences in the detection periods and the lack of correlations between salivary estradiol and androgen indicate that decreased estrogen concentrations in daughters of untreated or inadequately treated hypothyroid women result from impaired gonadal estrogen production and/or are secondary to increased activity of aromatase, but cannot be explained by altered androgen production. The association with impaired activation of gonadotropes is supported by positive correlations between salivary estradiol and urinary LH concentrations. However, LH concentrations were detectable in saliva for a shorter period of time than estradiol, and this observation may suggest either high sensitivity of ovarian follicles to this hormone (they may be activated by very low and undetectable LH levels), or the presence of additional activators of estradiol production. Moreover, the presented explanation is in line with observations that in early pregnancy, thyroid hormones (which during this period are almost exclusively of maternal origin) stimulate estradiol secretion by trophoblast cells and increase aromatase activity [
31]. Thus, it cannot be excluded that the impact of low thyroid hormone status on estrogen production persists in the first months after birth.
Differences in activation of the hypothalamic–pituitary–gonadal axis were paralleled by differences in the breast diameter and uterine length, which were greater in children born to either healthy mothers or to women with compensated thyroid function than in daughters of hypothyroid women who were untreated or undertreated during pregnancy. Moreover, both diameters positively correlated with estradiol concentrations, while in the offspring of women with normal gestational thyroid function there were also positive correlations between the breast size and FSH concentrations. Thus, minipuberty is associated with biological effects in the reproductive system, which may be less expressed if maternal hypothyroidism is during pregnancy poorly controlled. Interestingly, we observed an early reduction in uterine length and breast diameter, attributed to a withdrawal of placental estrogens in the first days of life, but less expressed than in the study by Kuiri-Hänninen et al. [
16]. This difference can be explained by a later beginning of our study compared to theirs, which was associated with a relatively long-term interval between cessation of the exposition to maternal estrogens and study commencement. Interestingly, between-group differences in the breast diameter and uterine length were statistically significant only from month 5. This finding indicates that biological effects of estrogens require long-term stimulation of estrogen receptors. Moreover, it may suggest that the hormonal milieu related to minipuberty plays a greater role in the development of target organs than the hormonal milieu associated with pregnancy. It is worth noting that, at the end of the study, both diameters were still smaller in the descendants of untreated or inadequately substituted hypothyroid women than in the remaining girls. The clinical relevance of this finding remains unknown and requires further research. Undoubtedly, it would be reasonable to specify whether this difference persists to the late prepubertal period and is implicated in the final (pubertal) stage of physical maturation.
Another interesting observation resulting from our study is that changes in activity of the hypothalamic–pituitary–ovarian axis and in the size of target organs seem to be a direct consequence of the low thyroid hormone status. Two findings allow us to draw this conclusion. Firstly, female descendants of women who had been effectively supplemented (adequate levothyroxine doses) during gestation were characterized by an unaltered course of minipuberty. Secondly, there were no difference between both groups of hypothyroid women (untreated/inadequately substituted and adequately substituted) in the relative prevalence of disorders leading to thyroid hypofunction. Thus, our findings indicate that potentially unfavorable hormonal and clinical changes associated with minipuberty in female descendants of women with impaired thyroid secretory function may be prevented if hypothyroidism is effectively controlled over the entire period of pregnancy.
Lastly, the obtained results provide evidence that the altered course of female minipuberty is not related to thyroid autoimmunity itself. This observation is in disagreement with the impact of autoimmunity on the pubertal phase of development. Thyroid hypofunction diagnosed in the prepubertal age, which, in this population, is usually secondary to autoimmune thyroid disease [
32], may delay puberty [
33] or result in incomplete isosexual precocity in girls (development of breast and internal genitalia without appearance of pubic hair) [
34]. Moreover, other autoimmune disorders developing before puberty are often complicated by delayed puberty, including juvenile idiopathic arthritis [
35], inflammatory bowel diseases [
36], celiac disease [
37], and type 1 diabetes [
38]. Moreover, the presence of autoimmunity in children was found to be accompanied by slow progression of puberty, isolated delayed menarche in females, and decreased duration and intensity of puberty growth [
35]. Two findings contradict the association between impaired hypothalamic–pituitary–ovarian axis activity in the early childhood and autoimmune thyroid disease. Firstly, there were no differences in hormone concentrations between infant girls born to mothers with thyroid hypofunction of autoimmune and non-autoimmune origin. Secondly, there were also no correlations between concentrations of gonadotropins and steroids and titers of thyroid antibodies. Thus, thyroid peroxidase and thyroglobulin antibodies, despite the ability to penetrate the placental barrier [
39], do not seem to have an impact the course of minipuberty in the offspring.
Some study limitations deserve to be noted. Although the study was sufficiently powered for group comparisons that were statistically significant, the relatively small sample size limits the extent to which our findings can be generalized. Because of its design (a small cohort study), we cannot exclude selection and confounding bias, which might have affected the obtained results. Relatively late recruitment of patients (the end of the first month of life) has not allowed us to investigate the earliest phase of minipuberty (the first few days after birth). Steroid hormone immunoassays used in the current study may theoretically interfere with compounds with structural similarities to the target steroid of the assay [
40], possibly elevating salivary hormone levels. Because the Polish population is iodine-sufficient and selenium-deficient, the effect of gestational hypothyroidism does not have to be the same in populations with inadequate iodine [
41] and sufficient selenium [
42] status. Strict study inclusion and exclusion criteria make it impossible to draw conclusions concerning the course of minipuberty in children born to hypothyroid mothers with comorbidities and/or chronically treated during pregnancy with other drugs. Lastly, taking numerous precautions during study design and data analysis limited, but did not completely eliminate, the regression toward the mean [
43].
4. Materials and Methods
This research was a single-center, prospective, matched, outpatient cohort study, carried out between May 2022 and April 2024. Owing to its nature, clinical trial registration was not required. The study design was reviewed and approved by the local bioethical commission (the Bioethical Committee of the Medical University of Silesia [approval number: PCN/CBN/0052/KB1/17/22, 12 April 2022), and the study was executed in compliance with the Declaration of Helsinki. Written informed consent was obtained from the parents of all infants participating in the study. Methodology and other study details adhere to the Strengthening the Reporting of Observational Studies in Epidemiology guidelines for reporting observational studies.
4.1. Participants
The research included a group of 96 female infants supervised by the principal study investigator (Karolina Kowalcze). The participants were chosen from a larger sample of children meeting the inclusion and exclusion criteria (n = 175) with the aim of matching the groups for women’s age, BMI, and gestational age of delivery. The matching algorithm has been described in our previous study [
19]. Based on maternal thyroid status during pregnancy, the study population was divided into three groups. Groups 1 and 2 encompassed female descendants of women with hypothyroidism. Thirty-two daughters of women who, during pregnancy, had been untreated or insufficiently treated were enrolled into group 1. In this group, circulating thyrotropin concentrations in mothers on two or more gestational measurements during pregnancy (at least six weeks apart) exceeded 4.0 mU/L. Thirty-two girls born to mothers adequately supplemented with levothyroxine were assigned into group 2. Adequate supplementation was defined as all thyrotropin concentrations (three or more) during gestation in the reference range for pregnant women established by the American Thyroid Association [
26] (between the lower limit of normal reduced by 0.4 mU/L and 4.0 mU/L in the first trimester, or between 0.5 and 4.0 mU/L in the remaining trimesters). Group 3 included 32 female descendants of healthy women. To be included into this group, the following criteria were required to be met: (a) thyrotropin concentrations on at least two measurements within the last 12 months (including at least one assessment in the first trimester of gestation) had to be between the lower limit of normal reduced in the first trimester of pregnancy by 0.4 mU/L and 4.0 mU/L; (b) free thyroxine, free triiodothyronine, and thyroid antibodies (thyroid peroxidase and thyroglobulin antibodies), if measured, had to be within normal limits; (c) the mother presented no signs and symptoms of thyroid dysfunction prior to or during gestation; and (d) the mothers had normal thyroid size and structure on ultrasound performed during pregnancy or at the first study visit (if thyroid ultrasound was not performed during pregnancy). The power analysis showed that 26 patients per group were needed to detect 25% between-group differences in gonadotropin and estradiol levels with a power of 80% and a two-tailed alpha level of 5%. Thus, considering the number of participants initiating and completed the study, the study was adequately powered to identify meaningful results.
The infants were excluded if (a) maternal thyrotropin concentration was elevated on only one occasion during gestation, (b) the woman was diagnosed with isolated hypothyroxinemia (free thyroxine levels below 10.2 pmol/L coexisting with thyrotropin concentrations within the reference range), or (c) hormonal results were inconclusive. We also excluded girls if their mothers (a) were diagnosed with other chronic diseases; (b) were admitted to hospital during gestation or postpartum owing to any acute complications; (c) during pregnancy or lactation, received any medication for more than 7 days (except for thyroid hormone preparations and/or vitamin/micronutrient supplements for pregnant and lactating women), or (d) were alcohol- or drug-addicted. The remaining exclusion criteria included the following: birth before gestational age of 36 weeks, a history of birth asphyxia, congenital infections, major congenital anomalies, metabolic diseases, genetic syndromes, other chronic disorders, and any chronic pharmacotherapy.
4.2. Study Design
The infants were followed-up once a month for the first six months of life, once every two months for the following six months, and once every three months until the age of eighteen months (11 visits). At each visit, the parents were asked to assess the child’s health status and medication usage. Moreover, the investigator thoroughly examined the infant and interpreted the results of laboratory and imaging tests (if done). Saliva and urine samples were collected only if the girl was considered healthy and, except for mandatory vaccination, did not receive any medicine for at least 10 days prior to the visit. The results were statistically analyzed only if the infant attended at least 8 of all scheduled visits.
All anthropometric measurements were carried out using standard procedures by the principle investigator. Recumbent length was measured from the top of the head to the sole of the feet with 1 mm precision using a mobile infantometer (Seca, Hamburg, Germany). Head circumference was measured using a flexible measuring tape as the maximum circumference between the supraorbital ridge and the occipital prominence. Weight was determined using a digital infant scale (Seca 834, Hamburg, Germany) with 10 g precision after removing clothing and diapers. Weight-for-length in infants was converted to percentiles using the World Health Organization growth charts. BMI (in mothers) was calculated as weight in kilograms divided by squared length in meters.
Uterine length and ovarian volume were measured by the principle investigator using a high-resolution transducer 5-MHz to 12-MHz (Esaote MyLab Six, Genoa, Italy). Three independent measurements were always taken, and the obtained results were averaged. Uterine length was measured in the long axis from the top of the fundus to the distal aspects of the cervix (external os). The volume of each ovary was calculated using the following formula: 0.52 × length × depth × width. The first two dimensions were determined on parasaggital and the third on transverse images. If both ovaries could be visualized and measured in all three dimensions, the final volume was the average value for both ovaries. Otherwise, the size of the only visible ovary (right or left) was used as the ovarian volume. The number of follicles was also recorded. The mammary gland was assessed according to the method described by Henriksen et al. [
44]. After identifying, by palpation, breast tissue as a firm subcutaneous disc, its diameter was measured to the nearest 1 mm using a caliper (Insize Europe, Zamudio, Spain). Diameters less than 3 mm were regarded as unmeasurable (as equivalent to the diameter of the nipple), and their size was assumed at 1 mm. Both glands were measured, and the mean diameter was the average of these two measurements.
4.3. Laboratory Assays
Urine and saliva samples were obtained between 7.00 and 8.30 a.m. Urine samples were collected using collection bags purchased from Medicavera (Szczecin, Poland). After cleaning and wiping the privates, the ends of the bag were attached on the skin above the privates and in the back of the vagina, while the sides next to each leg. Saliva samples were collected by the principal investigator by aspirating the saliva from the floor of the mouth with a 2 mL syringe. The procedure was performed at least one hour after the last feeding in order to prevent contamination and to limit a possible effect of feeding on hormone secretion [
19,
20]. Saliva collection, lasting for 25–50 s, was well tolerated and did not generate stress for the participants. The collected urine and saliva samples were then kept at −20 °C until analyzed.
All assays were carried out by a person who was unaware of the patient’s assignment. To increase the reliability of the results, the measurements were conducted in duplicate, and the obtained results were averaged. Urine FSH and LH levels were determined using a solid-phase, two-site chemiluminescent immunometric assay (Siemens Healthcare Diagnostics, Erlangen, Germany). Gonadotropin concentrations were then corrected for creatinine to standardize for the concentration of urine collected. Urine creatinine was measured using a modified Jaffe method (Roche Diagnostics, Basel, Switzerland). Concentrations of estradiol, testosterone, androstenedione, DHEA-S, progesterone, and 17-hydroxyprogesterone in saliva were measured by an enzyme-linked immuno-sorbent assay using reagents obtained from BioVendor R&D (Brno, Czech Republic), Diametra (Perugia, Italy), and IBL International (Hamburg, Germany). Limits of detection (LOD) were as follows: 4 pmol/L (estradiol), 10 pmol/L (testosterone), 18 pmol/L (androstenedione), 100 nmol/L (DHEA-S), 16 pmol/L (progesterone), 11 pmol/L (17-hydroxyprogesterone), 0.1 U/L (FSH), and 0.1 U/L (LH). The same methodology was used as in our previous studies [
19,
20].
4.4. Statistical Analysis
All variables were log-transformed to ensure that they closely followed a normal distribution. Within-patient differences at different time points were compared by repeated measures analysis of variance, followed by a post hoc analysis using Tukey’s Honestly Significantly Different test. Comparisons between the groups were performed by repeated measures analysis of variance, with post hoc analysis carried out using Bonferroni’s test. Associations between categorical variables were determined using the chi-square test. Correlations were evaluated using Pearson’s r tests (between two continuous variables), phi coefficient (between one continuous and one categorical variable), and point-biserial (between two categorical variables). The standardized mean difference (SMD) or the standardized difference of proportions (for percentage data) were used to calculate the effect size. In all analyses, differences were regarded as statistically significant if two-tailed p values were below 0.05 and/or 95% confidence intervals did not include the null value.