Association Between Neonatal Thyroid Function and Anogenital Distance from Birth to 48 Months of Age

Background: Evidence from animal studies has indicated that neonatal thyroid function is vital for the reproductive development. Anogenital distance (AGD), a sensitive biomarker of the fetal hormonal milieu, can be used to predict adult reproductive disorders. However, few human studies have examined the association between neonatal thyroid function and AGD. We aimed to explore their associations in a birth cohort study.

Methods: Concentrations of thyroid stimulating hormone (TSH) and thyroid hormones (THs), including total triiodothyronine (TT₃), total thyroxine (TT₄), free triiodothyronine (FT₃), and free thyroxine (FT₄) were measured in cord plasma in the Shanghai-Minhang Birth Cohort. The offspring AGD (AGD_AP [anus–penis] and AGD_AS [anus–scrotum] for boys and AGD_AC [anus–clitoris] and AGD_AF [anus–fourchette] for girls), body weight and anogenital index (AGI = AGD/weight [mm/kg]) were obtained at each follow-up visit. In total, 344 children (194 boys and 150 girls) with cord plasma concentrations of THs and TSH and at least one AGD measurement at birth and at 6, 12, and 48 months of age were included. Multiple linear regression and generalized estimating equation (GEE) models were used to examine the associations of cord plasma concentrations of THs and TSH with AGI.

Results: Multiple linear regression models showed inverse associations of TT₄, FT₃, and FT₄ with female AGI, although statistical significance was only reached at birth, 6 and 48 months of age. These associations were also found in GEE models: higher TT₄ and FT₄ concentrations were associated with lower AGI_AC (TT₄: β = -0.27, 95% CI: -0.50, -0.03 for middle vs. lowest tertile; FT₄: β = -0.38, 95% CI: -0.61, -0.16 for middle and β = -0.30, 95% CI: -0.55, -0.04 for highest vs. lowest tertile). Besides, girls with the highest tertile of FT₃ concentrations had lower AGI_AF than those with the lowest tertile (the highest vs. lowest tertile: β = -0.22, 95% CI: -0.36, -0.08). Positive associations between TSH and AGI at birth and at 12 months of age were observed in boys.

Conclusions: This study provides further evidence on the effects of neonatal thyroid function on reproductive development at an early life stage.

1 Introduction

Accumulating evidence from animal studies has indicated that neonatal thyroid function influences gonadal differentiation and reproductive function (1–4). Administration of large doses of thyroxine to neonatal female rats has been shown to delay vaginal opening and first estrus (5). Neonatal 3,5,3’-triiodothyronine excess in male rats was associated with a decrease in adult testis size (6), and might promote the apoptosis of germ cells in the neonatal testis (7). In humans, thyroid dysfunction is associated with decreased semen quality, decreased sexual activity, menstrual disturbances, and infertility (8–11). However, limited studies have evaluated the influence of thyroid function starting from birth on reproductive development since the information on reproductive development in early life is difficult to obtain, and decades may lapse between neonatal thyroid function and some manifest reproductive end points such as puberty development.

Anogenital distance (AGD), the distance from the center of the anus to the genital tubercle, has been widely recognized as a sensitive biomarker of the fetal hormonal milieu, and a measure of reproductive toxicity in animal models (12). Measurements of AGD at birth and during early childhood are also known to track through adulthood and can be used to predict various reproductive health disorders in adult humans (13, 14). Shorter AGD has been found to be associated with poor sperm quality (15), infertility (16), and lower testosterone concentrations (16) in men and associated with increased risk of endometriomas and deep infiltrating endometriosis in women (17). To date, only one human study has assessed the associations of thyroid stimulating hormone (TSH), free triiodothyronine (FT₃), and free thyroxine (FT₄) concentrations in umbilical cord serum with AGD at birth (18), indicating that thyroid function may be involved in male gonadal development in early life stages. However, its longitudinal effects at a later age are unknown.

In the present study, we examined the longitudinal associations between neonatal thyroid function, as reflected by the concentrations of TSH and thyroid hormones (THs, including total triiodothyronine (TT₃), total thyroxine (TT₄), FT₃, and FT₄)) in cord plasma, and repeated AGD measurements from birth to 48 months of age.

2 Materials and Methods

2.1 Study Design and Population

The Shanghai-Minhang Birth Cohort Study (S-MBCS) is an ongoing prospective cohort study designed to determine the distributions of a wide spectrum of maternal environmental exposures and examine their effects on pregnant women and their children (19, 20). From April to December 2012, pregnant women were recruited during their first prenatal care visit (at 12–16 weeks of gestation) at Maternal and Child Health Hospital of Minhang district in Shanghai, China. Briefly, women were eligible for inclusion if they were registered residents of Shanghai, had no history of hospital-diagnosed chronic disease, planned to be deliver their baby in the study hospital, and were willing to attend specified interviews during pregnancy and after delivery. In total, 1,292 eligible pregnant women agreed to participate in this study and completed a structured questionnaire at enrollment.

At delivery, 1,225 participants gave singleton live births, including 667 (54.4%) boys and 558 (45.6%) girls. A total of 611 women voluntarily provided cord blood samples during delivery. Due to limited funding, concentrations of THs and TSH were measured in a subset of 348 cord blood samples of infants with complete information at delivery, sufficient cord plasma volume (≥1 ml), and at least one follow-up at 12 or 48 months of age. Among them, two AGD metrics were measured for 186, 153, 154, and 163 boys, and for 144, 115, 105, and 124 girls at birth, 6, 12 and 48 months of age, respectively. Ultimately, 344 children (194 boys and 150 girls) who both have cord plasma concentrations of THs and TSH and at least one AGD measurement at the four time points were included in the present study, as shown in Figure 1.

FIGURE 1

Figure 1 Study population of the present study from the Shanghai-Minhang Birth Cohort Study.

The study protocol was approved by the ethical committee of Shanghai Institute of Planned Parenthood Research. All women provided informed consent for themselves and their children before participating in the S-MBCS and postnatal follow-up visits.

2.2 Measurement of Neonatal Thyroid Function

Cord blood samples were collected from the umbilical vein at the time of birth. The samples were centrifuged to retrieve plasma, which was immediately frozen at −80°C until shipment to the clinical laboratory of the affiliated hospital of Shanghai Institute of Planned Parenthood Research for the measurement of THs and TSH concentrations. An Electrochemiluminescence immunoassay with kits obtained from Roche Cobas e601 analyzer was performed to measure the concentrations of TSH, TT₄, FT₄, TT₃, and FT₃ (21). The analytical sensitivity was 0.3 nmol/L, 5.4 nmol/L, 0.4 pmol/L, 0.3 pmol/L, and 0.27 μIU/ml for TT₃, TT₄, FT₃, FT₄, and TSH, respectively. We also measured thyroid peroxidase antibody (TPOAb) concentration. TPOAb positivity was defined as a TPOAb concentration of ≥ 34 IU/L.

2.3 Measurement of AGD

The methods and procedures for the measurement of AGD have been described in detail elsewhere (22–24). In girls, AGD_AC was measured from the center of the anus to the anterior surface of the clitoral hood, and AGD_AF was measured from the center of the anus to the posterior end of the fourchette. In boys, AGD_AP was measured from the center of the anus to the anterior base of the penis where the penile tissue met the pubic bone, and AGD_AS was measured from the center of the anus to the posterior base of the scrotum where the skin changed from rugate to smooth. All children were placed in a dorsal decubitus position with both legs held back in a frog leg position to obtain accurate AGD measurements. Before each follow-up visit, we conducted standardized training to ensure the accuracy of AGD measurements. AGD measurements were conducted using the same Vernier caliper with increments of 0.1 mm in children during their first 3 days of life and at 6 months ± 2 weeks, 12 months ± 2 weeks, and 48 months ± 2 weeks of age by the same four trained examiners. To evaluate the inter-examiner variability, two examiners took independent measurements on each newborn on the same day using the same method in 15 girls and 15 boys. The intra-class correlation coefficients of AGD_AP, AGD_AS, AGD_AC and AGD_AF were 0.836, 0.731, 0.624, and 0.722, respectively, indicating moderate-to-good inter-rater reliability for the AGD measurements (25). The examiners who participated in the AGD measurements had no knowledge regarding the neonatal THs or TSH concentrations.

2.4 Covariates

The demographic characteristics, health conditions, medical and reproductive histories, and lifestyle factors (e.g., smoking and drinking habits) of pregnant women and their partners were obtained at enrollment through a structured questionnaire. After delivery, information on infant sex, birth weight, date of birth, mode of delivery, and gestational age was extracted from the electronic medical records of the study hospitals. Feeding practices and additional information were obtained from questionnaires filled by the mother or another caregiver at each postnatal follow-up visit.

2.5 Statistical Analyses

The distributions of participant characteristics were summarized as mean (standard deviation, SD) for continuous variables or counts and percentages for categorical variables. The distributions of TT₃, TT₄, FT₃, and FT₄ concentrations were normal, whereas that of TSH concentrations were skewed. Thus, TSH concentrations were log₁₀-transformed for further analyses. Since AGD was reported to be dependent on body weight, we used the anogenital index [AGI = AGD/weight (mm/kg)] as a weight-normalized index of AGD, as reported previously (26). We also used mean (SD), and percentiles to describe the distributions of AGI at birth and at 6, 12, and 48 months of age and to describe the distributions of neonatal TT₃, TT₄, FT₃, FT₄, and TSH concentrations.

2.5.1 Main Analyses of the Association Between Neonatal THs and TSH and AGI

A generalized additive model (GAM) was used to examine the nonlinearity assumptions between neonatal THs and log₁₀-transformed TSH concentrations and AGI. As non-linear associations were observed from GAM, the concentrations of TT₃, TT₄, FT₃, FT₄, and log₁₀-transformed TSH were treated as categorical variables (by tertiles) in the analyses. We first performed multiple linear regression models at each visit to assess the associations of THs and TSH concentrations with AGI. To take advantage of the longitudinal design and repeated measurements of AGD, generalized estimating equation (GEE) models were then applied to estimate β coefficients and 95% confidence intervals (CIs) for neonatal THs and TSH concentrations in relation to AGI after accounting for correlations of repeated AGD measures that had been obtained from four visits. We also added child age at follow-up into the GEE model and entered the interaction terms between THs and TSH concentrations (categorical) and child age (categorically) to assess differences in the associations over time. Since most interaction terms were not statistically significant (P > 0.1, Table S1), overall estimates were provided for AGI at 0–4 years of age.

Potential confounders were selected based on prior knowledge (27) or the change-in-estimate principle. For the latter criterion, the covariates were also retained in the final analyses if they caused a ≥10% change in the effect estimates for the associations of THs and TSH concentrations with AGI. Maternal education was also considered as a covariate since it is a good measure of socio-economic status (28) which can also indicate some often unobservable characteristics (29), and women’s socio-economic status were associated with the health endpoints of children (29, 30). Potential confounders finally considered in the models were maternal age at conception (years), maternal education (high school or below/college or above), maternal pre-pregnancy body mass index (BMI: <18.5, 18.5–24, and ≥24 kg/m²), maternal pre-pregnancy passive smoking (yes/no), paternal alcohol consumption before conception (yes/no), and gestational weeks (weeks).

2.5.2 Sensitivity Analyses

First, in addition to AGI, we used AGD as the dependent variable, and adjusted for weight-for-length z-scores in the GEE models to facilitate comparison with other studies (18). Second, TPOAb positivity is associated with a higher risk of developing autoimmune thyroiditis in children and adolescents (31), which may have adverse effects on reproductive outcomes (32). We thus excluded the subjects with TPOAb positivity (four girls and three boys) to test the robustness of the results. Third, to reduce the potential influence of incomplete neonatal thyroid development, we excluded subjects with gestational age <37 weeks. In addition, considering the significant variation in neonatal TSH and THs concentrations across different modes of delivery and maternal BMI (33–35), we re-ran the analyses in children who were delivered vaginally and whose mothers were of normal weight (18.5 < BMI < 24 kg/m²). We also re-ran the final analyses after excluding six children whose mothers reported prior or current diagnoses of any thyroid-related disease (e.g., autoimmune thyroid diseases, thyroid nodules, thyroid cancer, and thyroiditis)/any related medication use at enrollment to evaluate whether maternal thyroid-related disease/related medication use could affect the reliability of our results.

SAS 9.4 (SAS Institute Inc., Cary, NC, USA) was used for statistical analyses. A two-tailed value of p < 0.05 was considered statistically significant.

3 Results

3.1 Participant Characteristics

The characteristics of the 344 mother–child pairs are presented in Table 1. The mean (SD) age of mothers at conception was 28.09 (3.38) years and the mean (SD) gestational age was 39.67 (1.18) weeks. Most of mothers were Han Chinese (97.66%), nulliparous (86.88%), received a college degree or higher (78.13%), reported a monthly per capita household income >4000 RMB (79.53%), and had normal weight before pregnancy (72.78%). In addition, 58.31% mothers reported that they were not exposed to passive smoking before pregnancy, and 67.44% fathers reported no alcohol consumption within 3 months before conception. All newborns had a normal 5-min Apgar score (≥8). There was no statistically significant difference in demographic characteristics between the included mother–child pairs and those excluded (Table S2).

TABLE 1

Table 1 Descriptive characteristics of the study participants.

3.2 Distributions of AGI

All AGI measures were normally distributed. In girls, the mean (SD) AGI_AC was 8.75 (1.26), 4.52 (0.92), 4.06 (0.73), and 3.25 (0.67) mm/kg and the mean (SD) AGI_AF were 2.53 (0.70), 1.80 (0.58), 1.71 (0.52), and 1.78 (0.42) mm/kg at birth, 6 months of age, 12 months of age, and 48 months of age, respectively. In boys, the mean (SD) AGI_AP was 11.73 (1.54), 7.50 (1.25), 6.49 (1.09), and 5.29 (0.68) mm/kg and the mean (SD) AGI_AS was 4.32 (1.13), 3.10 (1.19), 2.79 (0.87), and 2.49 (0.53) mm/kg at birth, 6 months of age, 12 months of age, and 48 months of age, respectively (Table 2). The distributions of AGD from birth to 48 months of age are also shown in Figure S1.

TABLE 2

Table 2 Distributions of AGI and neonatal THs and TSH concentrations.

3.3 Distributions and Determinants of THs and TSH

In all subjects, the mean concentrations of TT₃, TT₄, FT₃, FT₄ and TSH in cord plasma were 0.86 nmol/L, and 94.12 nmol/L, 1.82 pmol/L, 14.22 pmol/L, and 6.47 uIU/mL, respectively, as shown in Table 2.

In univariate analyses (Table S3), children whose parents were younger than 25 years, whose fathers drank alcohol within 3 months before conception, and who were delivered vaginally (vs. cesarean section) had significantly higher neonatal TSH concentrations. Additionally, children born to mothers who were underweight (vs. normal weight) and gave birth via vaginal delivery had lower FT₃ and FT₄ concentrations. Further, children born to mothers who were exposed to passive smoking before pregnancy had higher FT₃ and FT₄ concentrations and children born to fathers who drank alcohol within 3 months before conception had higher FT₃ concentrations.

3.4 Associations of Neonatal THs and TSH Concentrations With AGI in Girls

In general, multiple linear regression models showed that AGI_AC and AGI_AF in girls tended to be lower in the higher tertile groups of TT₄, FT₄, and FT₃ in cord plasma than those of girls in the lowest tertile, although statistical significance was not reached for all time points (Figure 2). Girls who had higher concentrations of TT₄ had lower AGI_AC at birth (β = -0.56, 95% CI: -1.08, -0.04 for highest vs. lowest tertiles) and 6 months of age (β = -0.57, 95% CI: -1.03, -0.11 for middle vs. lowest tertiles), and lower AGI_AF at 6 months of age (β = -0.35, 95% CI: -0.65, -0.06 for middle vs. lowest tertiles). Higher FT₄ concentrations were associated with lower AGI_AC at birth (β = -0.58, 95% CI: -1.11, -0.05 for middle and β = -0.74, 95% CI: -1.26, -0.22 for highest vs. lowest tertiles) and 48 months of age (β = -0.35, 95% CI: -0.67, -0.02 for middle vs. lowest tertiles), and with lower AGI_AF at 6 months of age (β = -0.35, 95% CI: -0.65, -0.05 for middle vs. lowest tertiles) and 48 months of age (β = -0.22, 95% CI: -0.43, -0.02 for middle vs. lowest tertiles). In addition, girls in the highest tertile of FT₃ concentrations had lower AGI_AF at birth (β = -0.39, 95% CI: -0.68, -0.10) than those in the lowest tertile.

FIGURE 2

Figure 2 Associations between neonatal TH and TSH concentrations and AGI_AC and AGI_AF in girls from birth to 48 months of age. In multiple linear regression models. (A) AGI_AC [anus-clitoris] (B) AGI_AF [anus-fourchette]. All models adjusted for maternal age, maternal education, maternal pre-pregnancy BMI, gestational weeks, maternal passive smaoking, and paternal alcohol comsumption.

Similar inverse associations of TT₄, FT₄, and FT₃ with AGI in girls were also found in the GEE models (Table 3). Higher neonatal TT₄ and FT₄ concentrations were associated with lower AGI_AC (TT₄: β = -0.27, 95% CI: -0.50, -0.03 for middle vs. lowest tertiles; FT₄: β = -0.38, 95% CI: -0.61, -0.16 for middle and β = -0.30, 95% CI: -0.55, -0.04 for highest vs. lowest tertiles). Additionally, girls in the highest tertile of FT₃ concentrations had lower AGI_AF than those in the lowest tertile (β = -0.22, 95% CI: -0.36, -0.08).

TABLE 3

Table 3 Associations between TH and TSH concentrations in cord plasma and AGI from birth to 48 months of age in GEE models.

3.5 Associations of Neonatal THs and TSH Concentrations With AGI in Boys

For THs, the pattern of the results in multiple linear regression models (Figure 3) was inconsistent from birth to 48 months of age, although there were a positive association between TT₃ concentrations and AGI_AP at 6 months of age (β = 0.61, 95% CI: 0.10, 1.11 for highest vs. lowest tertiles) and an inverse association between TT₄ concentrations and AGI_AP at 6 months of age (β = -0.69, 95% CI: -1.19, -0.19 for middle vs. lowest tertiles).

FIGURE 3

Figure 3 Associations between neonatal TH and TSH concentrations and AGI_AP and AGI_AS in boys from birth to 48 months of age. in multiple linear regression models. (A) AGI_AP [anus-penis] (B) AGI_AS [anus-scrotum]. All models adjusted for maternal age, maternal education, maternal pre-pregnancy BMI, gestational weeks, maternal passive smaoking, and paternal alcohol comsumption.

Boys in the higher tertile of log₁₀-transformed TSH had higher AGI_AS at birth (β = 0.53, 95% CI: 0.13, 0.92 for middle and β = 0.50, 95% CI: 0.09, 0.91 for highest vs. lowest tertiles) and 12 months of age (β = 0.41, 95% CI: 0.04, 0.79 for highest vs. lowest tertiles; Figure 3). The positive association between TSH concentrations and AGI_AS remained in GEE models (β = 0.23, 95% CI: 0.0005, 0.46 for highest vs. lowest tertiles; Table 3), but the interaction term between TSH and child age at follow-up was statistically significant, indicating that the association between TSH and AGI_AS may change over time (P < 0.1; Table S2).

3.6 Sensitivity Analyses

We conducted several sensitivity analyses using GEE models. When we used AGD instead of AGI as the dependent variable, inverse associations between concentrations of TT₄, FT₃, and FT₄ and AGD were also observed in girls, but the association between TSH concentrations and AGD in boys became non-significant in GEE models (Table S4). The results were virtually identical in sensitivity analyses after excluding infants with TPOAb positivity, with gestational age <37 weeks, or those whose mothers reported prior or current diagnoses of thyroid-related disease/related medication use, and after restricting girls whose mothers had a normal weight (Tables S5–S8). Similar patterns were also observed in girls who were delivered vaginally, although some associations were not statistically significant, possibly due to the small sample size (Table S9).

4 Discussion

The present study is the first longitudinal assessment of the association between neonatal thyroid function and AGD from birth to 48 months of age. Higher neonatal concentrations of TT₄, FT₃, and FT₄ were associated with lower AGI in girls at 0–48 months of age; however, no consistent pattern for THs was found in boys. Positive associations between TSH concentrations and AGI at birth and at 12 months of age were observed in boys.

Numerous epidemiological studies have reported that alterations in neonatal THs and TSH concentrations are associated with infant neurodevelopment (21, 36) and physical growth (18, 37). However, evidence on the effects of THs on reproductive development is limited. Only one human study has assessed the effects of FT₃, FT₄, and TSH concentrations in cord serum on AGD at birth (18). In female newborns, the overall patterns of association between higher FT₃ and FT₄ and shorter AGD were observed in both our study and the study by Liu et al. (18), although the association was not statistically significant in the study by Liu et al. These differences are likely due to different neonatal TH concentrations (mean comparison: 1.82 vs. 1.63 for FT₃; 14.22 vs. 15.61 for FT₄). The concentrations of THs in the previous study were measured in cord serum and not in cord plasma, and the laboratory method used was also different (18). In line with the study by Liu et al. (18), we found a positive association between TSH concentration and AGI_AS in male newborns. We also found a positive association between TSH and AGI_AS at 12 months of age. Given that TSH concentrations at birth are more likely to be affected by factors, such as the timing of blood collection and temperature (38), further studies are necessary to confirm these results.

The inverse associations between THs and female AGI might be explained by the effects of THs on sex steroid hormone synthesis (1). Female AGD development during the fetal masculinization window was also suggested to be affected by androgen concentrations (39). In in vitro studies, THs along with other gonadotropic hormones could inhibit the excessive production of androgens in theca cells isolated from medium-sized follicles (40). Animal studies have also reported that THs could decrease the expression of steroid 5α-reductase type 2, one androgen-related gene, playing anti-androgen effects in the ovary tissue of female S. tropicalis frogs (41). Another mechanism that affecting female AGD through estrogenic pathways has also been proposed (39). THs could also up-regulate the expression of the estrogen receptor alpha and increase the SHBG concentrations in females, resulting in increased serum estrogen concentrations (42, 43). Male reproductive development is intricately dependent on fetal androgen action. Increased fetal androgen concentrations are the main cause for longer AGD in males (44). It is reported that increased TSH concentrations are associated with decreased sex-hormone binding globulin (SHBG) concentrations in males (42), which in turn, will increase the bioavailability of testosterone (43).

Our study has several strengths. First, this longitudinal study used repeated measurements of AGD at multiple ages of early life to examine the association between neonatal thyroid function and AGD, which enabled us to assess the average effect of neonatal thyroid function on AGD development from birth to 48 months of age. Second, the present study was based on a well-designed birth cohort study, which enabled us to obtain information on a wide range of covariates. In addition, the findings were strengthened when similar associations were observed through various analytic strategies.

Several limitations should also be acknowledged in the present study. First, the AGD metrics at each follow-up visit were measured only once, which may have resulted in misclassification. However, it has been reported that the reliability coefficient increased only slightly when the average of two or three repeated measurements for AGD was used (14). In addition, the examiners were blinded to data on neonatal thyroid function, and the misclassification would be non-differential and would bias the associations towards null. Second, the sample size was relatively small, and a total of 80 associations between neonatal thyroid function and AGI were tested using multiple linear regression models in the present study, which may have increased the type I error. However, consistent patterns of inverse associations between neonatal THs and female AGI in GEE models were also observed, which may be less prone to the issue of multiple comparisons when the estimates were calculated for the overall effects from birth to 48 months of age. Third, circulating THs and TSH of both maternal and fetal origin are present in the fetus from the second trimester onward (45). We did not measure maternal TH and TSH concentrations; thus, the contribution of maternal TH and TSH status to our findings cannot be evaluated. Finally, because we measured THs at birth rather than during the fetal masculinization window, reverse causality could probably occur considering that sex steroids can also affect the TH axis. However, the effects of THs on sex steroids synthesis and action are supported by considerable data while the evidence for the opposite direction of crosstalk is much more restricted (1).

In conclusion, inverse associations between neonatal THs and female AGI and a positive association between neonatal TSH and male AGI were indicated. These findings provide additional evidence for the effects of neonatal thyroid function on reproductive development during early childhood.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics Statement

The studies involving human participants were reviewed and approved by the ethical committee of Shanghai Institute of Planned Parenthood Research. Written informed consent to participate in this study was provided by the participants’ legal guardian/next of kin.

Author Contributions

ML: data collection, formal analysis, investigation, writing-original draft, and writing-review and editing. HL: design of the data collection instruments, data collection, and writing-review and editing. GHF laboratory analyses, and writing-review and editing. ZW: design of the data collection instruments, data collection, funding acquisition, and writing-review and editing. XS: writing-review and editing. AC: writing-review and editing. MM: study conception and design, project supervision, funding acquisition, and writing-review and editing. WY: study conception and design, project supervision, and writing-review and editing. All authors accept accountability for the overall work. All authors contributed to the article and approved the submitted version.

Funding

This work was supported by grants from the National key research and development program (2016YFC1000505); National Natural Science Foundation of China (22076123, 81903346); the Science and Technology Commission of Shanghai Municipality (20ZR1448000); Shanghai Sailing Program (18YF1420500).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fendo.2021.736505/full#supplementary-material

References

1. Duarte-Guterman P, Navarro-Martin L, Trudeau VL. Mechanisms of Crosstalk Between Endocrine Systems: Regulation of Sex Steroid Hormone Synthesis and Action by Thyroid Hormones. Gen Comp Endocrinol (2014) 203:69–85. doi: 10.1016/j.ygcen.2014.03.015

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Wagner MS, Wajner SM, Maia AL. The Role of Thyroid Hormone in Testicular Development and Function. J Endocrinol (2008) 199:351–65. doi: 10.1677/JOE-08-0218

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Dijkstra G, de Rooij DG, de Jong FH, van den Hurk R. Effect of Hypothyroidism on Ovarian Follicular Development, Granulosa Cell Proliferation and Peripheral Hormone Levels in the Prepubertal Rat. Eur J Endocrinol (1996) 134:649–54. doi: 10.1530/eje.0.1340649

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Romano RM, Gomes SN, Cardoso NC, Schiessl L, Romano MA, Oliveira CA. New Insights for Male Infertility Revealed by Alterations in Spermatic Function and Differential Testicular Expression of Thyroid-Related Genes. Endocrine (2017) 55:607–17. doi: 10.1007/s12020-016-0952-3

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Soliman FA, Reineke EP. Influence of Variations in Environmental Temperature and Thyroid Status on Sexual Function in Young Female Mice. Am J Physiol (1952) 168:400–5. doi: 10.1152/ajplegacy.1952.168.2.400

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Holsberger DR, Kiesewetter SE, Cooke PS. Regulation of Neonatal Sertoli Cell Development by Thyroid Hormone Receptor Alpha1. Biol Reprod (2005) 73:396–403. doi: 10.1095/biolreprod.105.041426

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Rijntjes E, Wientjes AT, Swarts HJ, de Rooij DG, Teerds KJ. Dietary-Induced Hyperthyroidism Marginally Affects Neonatal Testicular Development. J Androl (2008) 29:643–53. doi: 10.2164/jandrol.108.005108

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Lotti F, Maseroli E, Fralassi N, Degl'Innocenti S, Boni L, Baldi E, et al. Is Thyroid Hormones Evaluation of Clinical Value in the Work-Up of Males of Infertile Couples? Hum Reprod (2016) 31:518–29. doi: 10.1093/humrep/dev338

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Krassas GE, Pontikides N. Male Reproductive Function in Relation With Thyroid Alterations. Best Pract Res Clin Endocrinol Metab (2004) 18:183–95. doi: 10.1016/j.beem.2004.03.003

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Doufas AG, Mastorakos G. The Hypothalamic-Pituitary-Thyroid Axis and the Female Reproductive System. Ann NY Acad Sci (2000) 900:65–76. doi: 10.1111/j.1749-6632.2000.tb06217.x

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Krassas GE, Poppe K, Glinoer D. Thyroid Function and Human Reproductive Health. Endocr Rev (2010) 31:702–55. doi: 10.1210/er.2009-0041

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Welsh M, Saunders PT, Fisken M, Scott HM, Hutchison GR, Smith LB, et al. Identification in Rats of a Programming Window for Reproductive Tract Masculinization, Disruption of Which Leads to Hypospadias and Cryptorchidism. J Clin Invest (2008) 118:1479–90. doi: 10.1172/JCI34241

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Kahn LG, Philippat C, Nakayama SF, Slama R, Trasande L. Endocrine-Disrupting Chemicals: Implications for Human Health. Lancet Diabetes Endocrinol (2020) 8:703–18. doi: 10.1016/S2213-8587(20)30129-7

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Priskorn L, Petersen JH, Jorgensen N, Kyhl HB, Andersen MS, Main KM, et al. Anogenital Distance as a Phenotypic Signature Through Infancy. Pediatr Res (2018) 83:573–9. doi: 10.1038/pr.2017.287

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Mendiola J, Stahlhut RW, Jorgensen N, Liu F, Swan SH. Shorter Anogenital Distance Predicts Poorer Semen Quality in Young Men in Rochester, New York. Environ Health Persp (2011) 119:958–63. doi: 10.1289/ehp.1103421

CrossRef Full Text | Google Scholar

16. Eisenberg ML, Jensen TK, Walters RC, Skakkebaek NE, Lipshultz LI. The Relationship Between Anogenital Distance and Reproductive Hormone Levels in Adult Men. J Urol (2012) 187:594–8. doi: 10.1016/j.juro.2011.10.041

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Mendiola J, Sanchez-Ferrer ML, Jimenez-Velazquez R, Canovas-Lopez L, Hernandez-Penalver AI, Corbalan-Biyang S, et al. Endometriomas and Deep Infiltrating Endometriosis in Adulthood Are Strongly Associated With Anogenital Distance, A Biomarker for Prenatal Hormonal Environment. Hum Reprod (2016) 31:2377–83. doi: 10.1093/humrep/dew163

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Liu RJ, Xu XJ, Zhang YL, Zheng XB, Kim SS, Dietrich KN, et al. Thyroid Hormone Status in Umbilical Cord Serum Is Positively Associated With Male Anogenital Distance. J Clin Endocr Metab (2016) 101:3378–85. doi: 10.1210/jc.2015-3872

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Luan M, Liang H, Yang F, Yuan W, Chen A, Liu X, et al. Prenatal Polybrominated Diphenyl Ethers Exposure and Anogenital Distance in Boys From a Shanghai Birth Cohort. Int J Hyg Environ Health (2019) 222:513–23. doi: 10.1016/j.ijheh.2019.01.008

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Tian Y, Zhou Y, Miao M, Wang Z, Yuan W, Liu X, et al. Determinants of Plasma Concentrations of Perfluoroalkyl and Polyfluoroalkyl Substances in Pregnant Women From a Birth Cohort in Shanghai, China. Environ Int (2018) 119:165–73. doi: 10.1016/j.envint.2018.06.015

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Ji HL, Liang H, Wang ZL, Miao MH, Wang X, Zhang XT, et al. Associations of Prenatal Exposures to Low Levels of Polybrominated Diphenyl Ether (PBDE) With Thyroid Hormones in Cord Plasma and Neurobehavioral Development in Children at 2 and 4 Years. Environ Int (2019) 131:105010. doi: 10.1016/j.envint.2019.105010

CrossRef Full Text | Google Scholar

22. Xia RL, Jin LM, Li DK, Liang H, Yang F, Chen JP, et al. Association Between Paternal Alcohol Consumption Before Conception and Anogenital Distance of Offspring. Alcohol Clin Exp Res (2018) 42:735–42. doi: 10.1111/acer.13595

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Sathyanarayana S, Beard L, Zhou C, Grady R. Measurement and Correlates of Ano-Genital Distance in Healthy, Newborn Infants. Int J Androl (2010) 33:317–23. doi: 10.1111/j.1365-2605.2009.01044.x

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Salazar-Martinez E, Romano-Riquer P, Yanez-Marquez E, Longnecker MP, Hernandez-Avila M. Anogenital Distance in Human Male and Female Newborns: A Descriptive, Cross-Sectional Study. Environ Health (2004) 3:8. doi: 10.1186/1476-069X-3-8

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Sun XW, Liu C, Wang ZL, Yang F, Liang H, Miao MH, et al. Prenatal Exposure to Residential PM2.5 and Anogenital Distance in Infants at Birth: A Birth Cohort Study From Shanghai, China. Environ Pollut (2020) 264:114684. doi: 10.1016/j.envpol.2020.114684

CrossRef Full Text | Google Scholar

26. Swan SH, Main KM, Liu F, Stewart SL, Kruse RL, Calafat AM, et al. Decrease in Anogenital Distance Among Male Infants With Prenatal Phthalate Exposure. Environ Health Perspect (2005) 113:1056–61. doi: 10.1289/ehp.8100

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Swan SH, Sathyanarayana S, Barrett ES, Janssen S, Liu F, Nguyen RH, et al. First Trimester Phthalate Exposure and Anogenital Distance in Newborns. Hum Reprod (2015) 30:963–72. doi: 10.1093/humrep/deu363

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Winkleby MA, Jatulis DE, Frank E, Fortmann SP. Socioeconomic-Status and Health - How Education, Income, and Occupation Contribute to Risk-Factors for Cardiovascular-Disease. Am J Public Health (1992) 82:816–20. doi: 10.2105/Ajph.82.6.816

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Maternal Socio-Economic Status and the Well-Being of the Next Generation. (2017). Available at: https://www.iza.org/publications/dp/10714/maternal-socio-economic-status-and-the-well-being-of-the-next-generations (Accessed August 21, 2021).

Google Scholar

30. Luo ZC, Wilkins R, Kramer MS, Canadi FIHSG. Effect of Neighbourhood Income and Maternal Education on Birth Outcomes: A Population-Based Study. Can Med Assoc J (2006) 174:1415–21. doi: 10.1503/cmaj.051096

CrossRef Full Text | Google Scholar

31. Svensson J, Lindberg B, Ericsson UB, Ofsson P, Jonsson B, Ivarsson SA. Thyroid Autoantibodies in Cord Blood Sera From Children and Adolescents With Autoimmune Thyroiditis. Thyroid (2006) 16:79–83. doi: 10.1089/thy.2006.16.79

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Plowden TC, Connell MT, Hill MJ, Mendola P, Kim K, Nobles CJ, et al. Family History of Autoimmune Disease in Relation to Time-to-Pregnancy, Pregnancy Loss, and Live Birth Rate. J Transl Autoimmun (2020) 3:100059. doi: 10.1016/j.jtauto.2020.100059

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Derakhshan A, Shu H, Broeren MAC, de Poortere RA, Wikstrom S, Peeters RP, et al. Reference Ranges and Determinants of Thyroid Function During Early Pregnancy: The SELMA Study. J Clin Endocrinol Metab (2018) 103:3548–56. doi: 10.1210/jc.2018-00890

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Herbstman J, Apelberg BJ, Witter FR, Panny S, Goldman LR. Maternal, Infant, and Delivery Factors Associated With Neonatal Thyroid Hormone Status. Thyroid (2008) 18:67–76. doi: 10.1089/thy.2007.0180

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Fan P, Luo ZC, Tang N, Wang W, Liu Z, Zhang J, et al. Advanced Maternal Age, Mode of Delivery, and Thyroid Hormone Levels in Chinese Newborns. Front Endocrinol (Lausanne) (2019) 10:913. doi: 10.3389/fendo.2019.00913

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Williams FL, Watson J, Ogston SA, Visser TJ, Hume R, Willatts P. Maternal and Umbilical Cord Levels of T4, FT4, TSH, TPOAb, and TgAb in Term Infants and Neurodevelopmental Outcome at 5.5 Years. J Clin Endocrinol Metab (2013) 98:829–38. doi: 10.1210/jc.2012-3572

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Shields BM, Knight BA, Hill A, Hattersley AT, Vaidya B. Fetal Thyroid Hormone Level at Birth Is Associated With Fetal Growth. J Clin Endocrinol Metab (2011) 96:E934–8. doi: 10.1210/jc.2010-2814

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Ryckman KK, Berberich SL, Shchelochkov OA, Cook DE, Murray JC. Clinical and Environmental Influences on Metabolic Biomarkers Collected for Newborn Screening. Clin Biochem (2013) 46:133–8. doi: 10.1016/j.clinbiochem.2012.09.013

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Schwartz CL, Christiansen S, Vinggaard AM, Axelstad M, Hass U, Svingen T. Anogenital Distance as a Toxicological or Clinical Marker for Fetal Androgen Action and Risk for Reproductive Disorders. Arch Toxicol (2019) 93:253–72. doi: 10.1007/s00204-018-2350-5

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Gregoraszczuk EL, Skalka M. Thyroid Hormone as a Regulator of Basal and Human Chorionic Gonadotrophin-Stimulated Steroidogenesis by Cultured Porcine Theca and Granulosa Cells Isolated at Different Stages of the Follicular Phase. Reprod Fertil Dev (1996) 8:961–7. doi: 10.1071/rd9960961

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Campbell DEK, Langlois VS. Thyroid Hormones and Androgens Differentially Regulate Gene Expression in Testes and Ovaries of Sexually Mature Silurana Tropicalis. Gen Comp Endocrinol (2018) 267:172–82. doi: 10.1016/j.ygcen.2018.07.001

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Kjaergaard AD, Marouli E, Papadopoulou A, Deloukas P, Kus A, Sterenborg R, et al. Thyroid Function, Sex Hormones and Sexual Function: A Mendelian Randomization Study. Eur J Epidemiol (2021) 36:335–44. doi: 10.1007/s10654-021-00721-z

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Abalovich M, Levalle O, Hermes R, Scaglia H, Aranda C, Zylbersztein C, et al. Hypothalamic-Pituitary-Testicular Axis and Seminal Parameters in Hyperthyroid Males. Thyroid (1999) 9:857–63. doi: 10.1089/thy.1999.9.857

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Shah R, Alshaikh B, Schall JI, Kelly A, Ford E, Zemel BS, et al. Endocrine-Sensitive Physical Endpoints in Newborns: Ranges and Predictors. Pediatr Res (2021) 89:660–6. doi: 10.1038/s41390-020-0950-2

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Vulsma T, Gons MH, de Vijlder JJ. Maternal-Fetal Transfer of Thyroxine in Congenital Hypothyroidism Due to a Total Organification Defect or Thyroid Agenesis. N Engl J Med (1989) 321:13–6. doi: 10.1056/nejm198907063210103

PubMed Abstract | CrossRef Full Text | Google Scholar