- Open Access
Regulating role of fetal thyroid hormones on placental mitochondrial DNA methylation: epidemiological evidence from the ENVIRONAGE birth cohort study
Clinical Epigeneticsvolume 9, Article number: 66 (2017)
Fetal development largely depends on thyroid hormone availability and proper placental function with an important role played by placental mitochondria. The biological mechanisms by which thyroid hormones exert their effects on mitochondrial function are not well understood. We investigated the role of fetal thyroid hormones on placental mitochondrial DNA (mtDNA) content and mtDNA methylation. We collected placental tissue and cord blood from 305 mother–child pairs that were enrolled between February 2010 and June 2014 in the ENVIRONAGE (ENVIRonmental influence ON early AGEing) birth cohort (province of Limburg, Belgium). Placental mtDNA content was determined by qPCR and placental mtDNA methylation by bisulfite-pyrosequencing in two regions, i.e., the D-loop control region and 12S ribosomal RNA (MT-RNR1). The levels of free thyroid hormones (FT3, FT4) and thyroid-stimulating hormone (TSH) were measured in cord blood.
Cord blood FT3 and FT4 were inversely associated with placental mtDNA methylation at the MT-RNR1 (p ≤ 0.01) and D-loop (p ≤ 0.05) regions, whereas a positive association was observed for both hormones with placental mtDNA content (p ≤ 0.04). Assuming causality, we estimated that MT-RNR1 and D-loop methylation mediated, respectively, 77% [indirect effect +14.61% (95% CI 2.64 to 27.98%, p = 0.01)] and 47% [indirect effect +8.60% (95% CI 1.23 to 16.50%, p = 0.02] of the positive association between FT3 and placental mtDNA content. Mediation models with FT4 gave similar results but the estimated effect proportions were smaller compared with those of FT3 (54% and 24%, respectively).
We showed that epigenetic modification at specific loci of the mitochondrial genome could intervene with the thyroid-dependent regulation of mitochondrial DNA copy numbers.
Thyroid hormones are known to have profound effects on mitochondrial energetics and biogenesis [1, 2]. Under the control of thyroid-stimulating hormone (TSH), the thyroid gland produces thyroxine (T4), the major form of thyroid hormone, and triiodothyronine (T3), the active form. On account of its biological activity, unbound or free T3 (FT3) regulates mitochondriogenesis most likely in two different cell compartments: (i) binding of FT3 to thyroid receptors (TRs) that consecutively bind to response elements in the nucleus, activating the expression of mitochondria-related genes such as peroxisome proliferator-activated receptor γ-coactivator1α (PPARGC1A) or (ii) direct binding to specific TRs in the mitochondrial matrix stimulating transcription of the mitochondrial genome .
The placenta is a relevant tissue to investigate the interplay between thyroid hormones and mitochondria. This metabolically active organ contains a high number of mitochondria and regulates nutrient, oxygen, and hormonal transfer, allowing for fetal growth. It is also well established that thyroid hormones are critical for placental  and fetal development . FT3 stimulates the production of factors that control trophoblast growth and development . The placenta expresses two types of iodothyronine deiodinases (D2 and D3) that are capable of metabolizing FT4 to FT3 , and thus, plays an important role in thyroid hormone homeostasis.
Although some studies are questioning the presence of DNA methylation in the mitochondrial genome [8, 9], in the last few years, mtDNA methylation has emerged as a next-generation biomarker that might be a diagnostic tool for several diseases such as cardiovascular disease . Here, we examined two specific mtDNA loci, i.e., the displacement loop (D-loop) and MT-RNR1. DNA methylation levels in the D-loop control region are suggested to play an important role in modulating either the replication or transcription of mtDNA since nearly the entire mitochondrial genome transcribes from this region . Our other studied locus MT-RNR1 encodes for the 12S rRNA protein which is critical for normal function and integrity of the mitochondrial ribosome . In this mother–newborn study, we aimed to elucidate possible mechanisms involved in the regulation of mitochondrial biogenesis by investigating the association between fetal FT3, FT4, and TSH with mtDNA content and mtDNA methylation in placental tissue. We hypothesized that methylation at specific loci of the mitochondrial genome could interfere with thyroid hormone-dependent regulation of mitochondrial biogenesis.
Mother–newborn characteristics and demographics
The study included 305 mother–child pairs (mean maternal age, 29.1 yr; range, 18–42 years). Demographic and prenatal lifestyle factors are reported in Table 1. Briefly, mean (10th–90th percentile) pre-pregnancy body mass index (BMI) of the participating mothers was 24.1 (19.7–29.7) kg/m2. Fifty-two mothers (17.1%) reported to have smoked during pregnancy and smoked on average 7.8 cigarettes per day. Most women (64.9%, n = 198) never smoked cigarettes. The majority (>50%) of the mothers were highly educated. Nearly half of the newborn population were boys (n = 151; 49.5%). The overall mean gestational age was 39.2 weeks (10th–90th percentile 18–41) and included a vast majority of primiparous (51.8%, n = 158) or secundiparous (36.7%, n = 112) newborns.
Levels of fetal thyroid hormones, placental mtDNA methylation, and mtDNA content
The geometric means (10th–90th percentile) of thyroid hormone levels in cord blood were 2.63 (2.15–3.22) pmol/L for FT3, 15.66 (13.26–18.53) pmol/L for FT4, and 11.65 (5.16–20.59) mU/L for TSH (Table 2). A positive correlation was observed between the FT3 and FT4 values (r = 0.25; p < 0.0001) and between the FT3 and TSH values (r = 0.19; p = 0.0008) (Fig. 1).
Table 2 also shows mtDNA methylation and mtDNA content values measured in placental tissue. Placental mtDNA methylation averaged (10th–90th percentile) 9.51% (4.46–14.67) for MT-RNR1 and 3.61% (1.94–5.24) for D-loop. The methylation levels of MT-RNR1 and D-loop were strongly correlated (r = 0.85, p < 0.0001) (Fig. 1 and Additional file 1: Figure S1). The geometric mean of placental mtDNA content was 1.11 (0.49–2.49) (unitless). In this subset, and as previously published , we report a strong inverse correlation between placental mtDNA methylation and mtDNA content (r = −0.63, p < 0.0001 for MT-RNR1 and r = −0.40, p < 0.0001 for D-loop) (Fig. 1).
Association of fetal thyroid hormones with placental mtDNA methylation and mtDNA content
Both fetal FT3 and FT4 were inversely correlated with mtDNA methylation at the MT-RNR1 and D-loop region, whereas a positive correlation was observed for both hormones with mtDNA content in placental tissue (see Fig. 1 for correlations). Fetal TSH levels did not correlate with either placental mtDNA methylation or mtDNA content.
The associations remained significant after adjustment for maternal age, pre-pregnancy BMI, gestational age, newborn’s sex, smoking status, parity, maternal education, ethnicity, and cord plasma insulin level. A 10th–90th percentile increment (53%) in cord blood FT3 (log10 values) was associated with a lowering in absolute methylation of −1.50% (95% CI −2.70 to −0.30%, p = 0.01) for MT-RNR1 and −0.46% (95% CI −0.83 to −0.09%, p = 0.02) for D-loop. A similar association, but with smaller estimates, was observed for FT4 (Table 3). On the other hand, a 10th–90th percentile increment in cord blood FT3 (53%) or FT4 (13%) was associated with a higher placental mtDNA content (relative change) of +20.35% (95% CI 0.47 to 44.17%, p = 0.04) for FT3 and +11.04% (95% CI 4.00 to 18.54%, p = 0.002) for FT4.
We performed mediation analysis to estimate the proportion of the associations between cord blood thyroid hormones and placental mtDNA content that might be mediated by changes in mtDNA methylation if the underlying causal assumptions of the mediation analysis are valid. We selected both MT-RNR1 and D-loop as potential mediators because methylation levels at both regions were significantly associated with FT3 and FT4 as well as with placental mtDNA content (Fig. 1). While adjusting for the aforementioned variables, we estimated that MT-RNR1 methylation mediated 77% [indirect effect +14.61% (95% CI 2.64 to 27.98%, p = 0.01)] and D-loop methylation mediated 47% [indirect effect +8.60% (95% CI 1.23 to 16.50%, p = 0.02)] of the positive association between FT3 and placental mtDNA content (Fig. 2). Mediation models with FT4 gave similar results but the estimated effect proportions were smaller compared with those of FT3 (54 and 24%, respectively) (Additional file 1: Figure S2).
Certain CpGs were more associated with fetal thyroid hormones compared to others, underscoring the importance of methylation levels at specific CpG (Additional file 2: Table S1). For example, FT3 was associated with methylation levels at all three CpG sites of D-loop (the transcription start site of the mitochondrial genome).
Excluding women who had undergone a cesarean section (n = 13) slightly increased our effect estimates for FT3 (Additional file 2: Table S2), whereas excluding small for gestational age infants (n = 22) slightly decreased our effect estimates for FT3 (Additional file 2: Table S3). However, we verified whether the association between fetal thyroid hormones and placental mtDNA content was still present in the initial group of 547 mother–child pairs (flowchart Additional file 1: Figure S3). Adjusting for the aforementioned variables, the association between cord blood FT3 and placental mtDNA content was stronger compared with the reported association in the main analysis, i.e., relative change of +29.22% (95% CI 16.22 to 43.67%, p < 0.0001) for a 10th–90th percentile increment in cord blood FT3, whereas the association for FT4 was less (+6.65%, 95% CI 2.23 to 11.26%, p = 0.003).
Furthermore, additional adjustments for cord blood plasma estradiol (n = 304), passive indoor tobacco smoke exposure (n = 295), alcohol consumption (n = 295), or pH of the arterial cord blood (indicator of hypoxemia) (n = 264), did not alter our associations between fetal thyroid hormones and mtDNA content or mtDNA methylation (data not shown).
Mitochondrial biogenesis relies on a tightly coordinated process between the nuclear and mitochondrial genome. For example, the nuclear gene PPARGC1A is a central regulator of mitochondrial gene expression and biogenesis that is also controlled by FT3 . We observed an inverse association between fetal FT3 levels and placental PPARGC1A promoter methylation, i.e., −2.31% (95% CI −4.04 to −0.59, p = 0.009) for a 10th–90th percentile increment in FT3 (Additional file 2: Table S1). No association was observed between placental PPARGC1A promoter methylation and mtDNA content (p = 0.23) or mtDNA methylation (D-loop, p = 0.54; MT-RNR1, p = 0.74).
In this study, we report for the first time associations between fetal thyroid hormones and epigenetic modification at specific loci in the mitochondrial genome that could, at least in part, mediate the fetal thyroid-dependent regulation of mitochondrial biogenesis in placental tissue.
From mid-gestation onwards, the fetus starts secreting small amounts of thyroid hormones in addition to the transplacental supply of maternal thyroid hormones . The optimal concentration of fetal thyroid hormones at each stage of development is maintained by the placenta, in conjunction with the fetal thyroid gland, liver, and kidneys . The placenta expresses two types of iodothyronine deiodinases (D2 and D3) that are capable of metabolizing FT3 and FT4 , and thus, plays an important role in thyroid hormone homeostasis. Besides regulating fetal development , thyroid hormones have profound effects on mitochondrial energetics and biogenesis [1, 2].
Despite the spectacular progress in the knowledge of the T4-T3 nuclear pathway, a clear answer to how it regulates mitochondrial biogenesis is lacking. While it is known that an upregulation of nuclear-encoded respiratory genes including PPARGC1A occurs within hours after injection of T3 in hypothyroid rats, other non-genomic direct effects are detectable within minutes suggesting another mode of action [16, 17]. T3 binds directly to specific receptors inside the mitochondria. A 43 kDa c-Erb A alpha1 protein (p43), located exclusively in the mitochondrial matrix, acts as a T3-dependent transcription factor and specifically binds to four mitochondrial DNA sequences with a high similarity to nuclear T3 response elements [18, 19]. Overexpression of p43 increases mitochondrial genome transcription and protein synthesis, stimulating mitochondrial biogenesis in a T3-dependent manner . Interestingly, of the four mitochondrial DNA sequences, two response elements are located in the D-loop region and one in the MT-RNR1 gene. We hypothesized that DNA methylation at these mitochondrial hotspots could interfere with thyroid hormone-dependent regulation of mitochondrial biogenesis. First, we have shown that fetal thyroid hormones, especially FT3, are positively associated with placental mtDNA content and inversely associated with placental methylation levels at the D-loop and MT-RNR1 region. Next, we underscored a substantial mediating role of placental mtDNA methylation between the association of fetal thyroid hormones and placental mtDNA content. Hence, we postulate that high DNA methylation levels in the mitochondrial genome are related to conformational or structural changes making the mtDNA less accessible to proteins and transcription factors such as the T3-dependent transcription factor p43. Previously, we have shown that exposure to particulate matter (PM2.5) during pregnancy is associated not only with changes in fetal thyroid hormones  but also with higher levels of methylation of the D-loop and the MT-RNR1 sequence, affecting mtDNA content in placental tissue . An increase in mtDNA methylation was observed in blood leukocytes of steel workers exposed to air pollution  and in an experimental study of human cultured cells treated with ethidium bromide . The ethidium-bromide-exposed cells recovered from mtDNA depletion and showed increased overall methylation levels indicating that mtDNA was less packed with proteins during active mtDNA replication .
Our mediation analysis revealed that the effect of FT3 on mtDNA content is mediated for 77% by MT-RNR1 and for 54% by D-loop methylation levels, meaning that there are alternative routes, most likely nuclear pathways, by which thyroid hormones exert their action on mitochondrial biogenesis. T3 controls the expression of PPARGC1A , a transcriptional co-activator of several nuclear-encoded transcription factors including mitochondrial transcription factor A (TFAM) that regulates mitochondrial biogenesis  (Fig. 3). In addition to specific mitochondrial actions of FT3, we also observed an inverse association between placental PPARGC1A promoter methylation and fetal FT3 levels. Most likely this is indicative of upregulated PPARGC1A mRNA levels in the presence of high FT3 concentrations since promoter hypomethylation is usually associated with increased gene expression due to unwound DNA that is available for transcription factors. Interestingly, our interrogated region is a putative binding site for CREB transcription factors. The transcriptional activity of CREB is critical for the establishment and maintenance of energy homeostasis in mice neonates  and appears to be involved in the regulation of PPARGC1A. It is noteworthy to mention that the promoter region of PPARGC1A was largely unmethylated in cord blood but not in placental tissue , suggesting differential regulation of PPARGC1A in cord blood and placental tissue. Even though nuclear-encoded genes that are involved in mitochondrial function are regulated by epigenetic mechanisms , it would be interesting to determine the direct functional significance at mRNA expression level in response to FT3 levels. Unfortunately, we did not measure mRNA expression due to lack of suitable RNA samples and consider this a limitation of the study. We also have no clear information on “threshold” changes in mtDNA methylation that lead to functional effects of mitochondriogenesis. The mitochondrial genome, as far as we observed, shows a far less DNA methylation level (range 0 to 40%) compared with nuclear DNA (range 0 to 100%). Hence, small changes in mtDNA methylation might still have functional effects. It has been shown that small changes in single CpG sites and regional methylation changes interfere with gene/protein expression  by altering the affinity of transcription factors to their binding sites or chromosome looping events . Even when small changes in methylation would not impact functional effects, mtDNA methylation can provide invaluable help in developing new and useful markers of exposure and disease . However, our findings should first be confirmed in an experimental study and we need to investigate the health consequences in later life before we can translate it into a clinical meaning.
We are aware of the issue of placental tissue heterogeneity but it is difficult to assess the impact of this variation on our methylation results. To our knowledge, there is only one study that explored genome-wide DNA methylation patterns in the two main cell types of the placenta (cytotrophoblasts and fibroblasts) . These cell types did not demonstrate genome-wide differentiation in DNA methylation but some specific genes showed differential promoter methylation. Although these authors demonstrated that the methylation profile of the placenta is mainly driven by cytotrophoblasts, there is a possibility that our epigenetic data is confounded by variation in cell type distributions. To overcome this obstacle, one should correct for differences in tissue composition using algorithms that estimate cell type proportions based on DNA methylation signatures derived from Illumina microarrays [31, 32]. However, statistical adjustments of DNA methylation data for cell distributions require the availability of reference epigenomes for the component cell types created by cell sorting . Currently, there is no algorithm available to estimate the cell type proportions in placental tissue. Another limitation of our study is that we were not able to measure thyroid hormones in placental tissue but we used cord blood levels as a proxy since it circulates through the fetal side of the placenta.
Given our epidemiological findings and other experimental research data, it seems that there exist coordinated events between mtDNA methylation, mtDNA content, and thyroid hormones. We summarized our findings in a simplified scheme (Fig. 3) indicating that epigenetic modification at specific loci of the mitochondrial genome could intervene with thyroid-dependent regulation of mitochondrial biogenesis. Our findings could contribute to further epidemiological understanding of mitochondrial disorders. Whether alterations in mitochondrial function or newborn’s thyroid hormone levels have health consequences later in life should be elucidated.
Within the on-going ENVIRONAGE birth cohort (ENVIRonmental influence ON early AGEing) , we conducted our investigation in a group of 589 singleton pregnancies for which the mothers agreed withdrawal of cord blood, the collection of the placenta after delivery, and the use of information from their medical files. The placenta could not be collected for ten newborns, four placentas had insufficient DNA yield, 12 placentas had missing measurements of mtDNA content, and 16 mothers with thyroid gland complications were excluded. Because mtDNA methylation was measured only in a subset of this large study sample (n = 547), we ended with a final sample size of 305 newborns for the main analysis (see flowchart; Additional file 1: Figure S3).
Mother-child pairs were recruited from February 2010 to June 2014 at the East-Limburg Hospital in Genk (Belgium) following procedures approved by the Ethical Committee of Hasselt University and the East-Limburg Hospital and according to the principles outlined in the Helsinki Declaration for investigation of human subjects as previously described in details . Briefly, written informed consent was obtained from eligible participants before delivery. Questionnaires and medical records consulted after birth provided information on maternal age, maternal education, smoking status, ethnicity, pre-pregnancy BMI, gestational age, newborn’s sex, Apgar scores, birth weight and length, parity, and ultrasonographic data. Information about maternal tobacco smoke exposure was obtained by asking whether mothers smoked during pregnancy, whether they smoked at any time during their life, or if they never smoked in their life. All neonates were assessed for congenital anomalies immediately after birth and were considered healthy with an Apgar score after 5 min ranging between 7 and 10. No neonate was delivered in the Neonatal Intensive Care Unit. The ENVIRONAGE birth cohort generally consists of mothers with normal pregnancies without complications and with healthy neonates.
Cord blood collection and thyroid hormones measurements
Immediately after delivery, the umbilical cord was clamped and cord blood was drawn in plastic BD Vacutainer® Lithium Heparin Tubes (BD, Franklin Lakes, NJ, USA). Blood tubes were centrifuged at 3200 rpm for 15 min to retrieve plasma which was instantly kept at −80 °C. FT4 (pmol/L), FT3 (pmol/L), TSH (mU/L), and insulin were measured in plasma using an electro-chemiluminescence immunoassay using the Modular E170 automatic analyzer (Roche, Basel, Switzerland) at the clinical lab of East-Limburg Hospital.
Placentas were deep-frozen within 10 min of delivery, and afterwards, placental specimens were taken for DNA extraction after minimally thawing of the placentas. We took villous tissue (1 to 2 cm3) at a fixed location from the fetal side of the placenta, approximately 1–1.5 cm below the chorio-amniotic membrane, and preserved the biopsies at −80 °C . Genomic DNA was isolated from placental tissue using the QIAamp DNA mini kit (Qiagen, Inc., Venlo, the Netherlands) and stored at −80 °C until further use.
DNA methylation analyses
In 305 placentas, we performed DNA methylation analysis by highly quantitative bisulfite-PCR pyrosequencing as previously described . Briefly, bisulfite conversions were performed using 1 μg of extracted genomic DNA with the EZ-96 DNA methylation Gold kit (Zymo Research, Orange, CA, USA) according to the manufacturer’s instructions. We interrogated CpG sites within specific regions of the mitochondrial genome (MT-RNR1 and D-loop) and the promoter of PPARGC1A as described by Byun et al.  and Janssen et al. . Detailed information regarding primer sequences is given in Additional file 2: Table S4. PCR amplification of regions of interest prior to pyrosequencing was performed in a total reaction volume of 30 μl, using 15 μl GoTaq Hot Start Green Master Mix (Promega, Madison, WI, USA), 10 pmol forward primer, 10 pmol reverse primer, 1 μl bisulfite-treated genomic DNA, and water. PCR products were purified and sequenced by pyrosequencing using the PyroMark Q96 MD Pyrosequencing System (Qiagen, Inc., Germantown, MD, USA). The degree of methylation was expressed as the percentage of methylated cytosines over the sum of methylated and unmethylated cytosines. Samples were run in duplicate on two different plates from which the average methylation levels were used. The coefficient of variation was 3.7% for MT-RNR1, 5.8% for D-loop, and 8.1% for PPARGC1A. The efficiency of the bisulfite-conversion process was assessed using non-CpG cytosine residues within the sequence. We used 0% (PSQ-T oligo: 5′-TTGCGATACAACGGGAACAAACGTTGAATTC-3′) and 100% (PSQ-C oligo: 5′-TTGCGATACGACGGGAACAAACGTTGAATTC-3′) DNA methylation control oligos. The sequencing primer for the control oligo was 5′-AACGTTTGTTCCCGT-3′. We mixed the PSQ-C oligo (or PSQ-T oligo) with the sequencing oligo in PyroMark Annealing Buffer (Qiagen, Inc., Valencia, CA, USA) and performed pyrosequencing with the sequencing entry C/TGTAT. The between- and within-placenta variability, exemplified by the intra-class correlation coefficient, was evaluated in a subset of 19 placentas and was 58 vs. 42% (p = 0.009) for MT-RNR1, 61 vs. 39% (p = 0.01) for the D-loop region and 64 vs. 36% (p = 0.005) for PPARGC1A .
mtDNA content analysis
mtDNA content was measured in placental tissue by determining the ratio of two mitochondrial gene copy numbers (MTF3212/R3319 and MT-ND1) to two single-copy nuclear control genes (RPLP0 and ACTB) using a quantitative real-time polymerase chain reaction (qPCR) assay as previously described  but with small modification. Briefly, 2.5 μl diluted genomic DNA (5 ng/μl) was added to 7.5 μl mastermix consisting of Fast SYBR® Green I dye 2× (5 μl/reaction), forward and reverse primer (each 0.3 μl/reaction), and RNase free water (1.9 μl/reaction). Primer sequences (Additional file 2: Table S5) were diluted to a final concentration of 300 nM in the master mix. Samples were run in triplicate in 384-well format. Real-time PCR was performed using the 7900HT Fast Real-Time PCR System (Applied Biosystems, Foster City, CA, USA) with the following thermal cycling profile: 20 s at 95 °C (activation), followed by 40 cycles of 1 s at 95 °C (denaturation) and 20 s at 60 °C (annealing/extension), ending with melting curve analysis (15 s at 95 °C, 15 s at 60 °C, 15 s at 95 °C). qBase software (Biogazelle, Zwijnaarde, BE) was used to normalize data and correct for run-to-run differences .
For database management and statistical analysis, we used the SAS software program (version 9.4; SAS Institute Inc., Cary, NC, USA). mtDNA content and thyroid hormone levels were log10-transformed to improve normality. For each subject, we measured in placental tissue methylation levels of CpG sites at two regions of the mitochondrial genome and at the promoter region of PPARGC1A. The pyrosequencing-based DNA methylation analysis produced a methylation value (%) for each CpG of MT-RNR1 (two CpGs), the D-loop region (three CpGs), and PPARGC1A (three CpGs). In the main analysis, we used the average methylation levels of the different CpGs. Pearson correlation coefficients were calculated between the different thyroid hormone levels in cord blood (FT3, FT4, and TSH), mtDNA methylation (MT-RNR1 and D-loop), and mtDNA content in placental tissue using R software packages. We performed multiple linear regression to determine the association between thyroid hormones and placental mtDNA methylation (both MT-RNR1 and D-loop) and between thyroid hormones and placental mtDNA content. Thyroid hormones were fitted as linear variables in the models, and effect estimates on mtDNA methylation and mtDNA content were calculated for a 10th–90th percentile increment in thyroid hormones, which corresponds to a 52% change in FT3, a 12% change in FT4, and a 84% change in TSH. All models were adjusted for maternal age, pre-pregnancy BMI, gestational age, newborn’s sex, smoking status, parity, maternal education, ethnicity, and cord plasma insulin level. The Shapiro-Wilk statistic and Q-Q plots of the residuals were used to test the assumptions of all linear models.
We used mediation analysis to investigate potential associations that may underlie the relation between the exposure variable (FT3 or FT4) and the continuous outcome variable (mtDNA content) by examining how they relate to a third variable, the mediator (mtDNA methylation) . We accomplished this by decomposing the total effect into a direct effect (DE; exposure effect on outcome at a fixed level of the mediator) and an indirect effect (IE; exposure effect on outcome that operates through the mediator).
To underscore the importance of single CpG sites, we performed separate linear regression analyses between levels of FT3 or FT4 and the methylation levels at each CpG site of MT-RNR1 and D-loop. Furthermore, it is known that cord blood thyroid hormone levels are influenced not only by different external factors such as fetal stress exemplified by cesarean section  but also by other factors such as cord plasma estradiol , passive smoking , alcohol consumption , and pH of arterial cord blood . Besides adjusting the main model with the aforementioned variables, we performed additional analyses in which we excluded women who had undergone a cesarean section or excluding infants born with a birth weight less than the 10th percentile. Lastly, we used linear regression analyses to investigate the association between fetal thyroid hormones and placental promoter methylation of the nuclear gene PPARGC1A.
Harper M-E, Seifert EL. Thyroid hormone effects on mitochondrial energetics. Thyroid. 2008;18(2):145–56.
Goglia F, Moreno M, Lanni A. Action of thyroid hormones at the cellular level: the mitochondrial target. FEBS Lett. 1999;452(3):115–20.
Cheng S-Y, Leonard JL, Davis PJ. Molecular aspects of thyroid hormone actions. Endocr Rev. 2010;31(2):139–70.
Chan SY, Vasilopoulou E, Kilby MD. The role of the placenta in thyroid hormone delivery to the fetus. Nat Clin Pract End Met. 2009;5(1):45–54.
Burrow GN, Fisher DA, Larsen PR. Maternal and fetal thyroid function. N Engl J Med. 1994;331(16):1072–8.
Maruo T, Matsuo H, Mochizuki M. Thyroid hormone as a biological amplifier of differentiated trophoblast function in early pregnancy. Acta Endocrinol (Copenh). 1991;125(1):58–66.
Koopdonk-Kool JM, de Vijlder JJ, Veenboer GJ, Ris-Stalpers C, Kok JH, Vulsma T, et al. Type II and type III deiodinase activity in human placenta as a function of gestational age. J Clin Endocrinol Metab. 1996;81(6):2154–8.
Liu B, Du Q, Chen L, Fu G, Li S, Fu L, et al. CpG methylation patterns of human mitochondrial DNA. Sci Rep. 2016;6:23421.
Hong EE, Okitsu CY, Smith AD, Hsieh C-L. Regionally specific and genome-wide analyses conclusively demonstrate the absence of CpG methylation in human mitochondrial DNA. Mol Cell Biol. 2013;33(14):2683–90.
Baccarelli AA, Byun H-M. Platelet mitochondrial DNA methylation: a potential new marker of cardiovascular disease. Clin Epigenetics. 2015;7(1):1–9.
Aloni Y, Attardi G. Expression of the mitochondrial genome in HeLa cells. II. Evidence for complete transcription of mitochondrial DNA. J Mol Biol. 1971;55(2):251–67.
Metodiev MD, Lesko N, Park CB, Cámara Y, Shi Y, Wibom R, et al. Methylation of 12S rRNA is necessary for in vivo stability of the small subunit of the mammalian mitochondrial ribosome. Cell Metab. 2009;9(4):386–97.
Janssen BG, Byun H-M, Gyselaers W, Lefebvre W, Baccarelli AA, Nawrot TS. Placental mitochondrial methylation and exposure to airborne particulate matter in the early life environment: an ENVIRONAGE birth cohort study. Epigenetics. 2015;10(6):536–44.
Ventura-Clapier R, Garnier A, Veksler V. Transcriptional control of mitochondrial biogenesis: the central role of PGC-1α. Cardiovasc Res. 2008;79(2):208–17.
Thorpe-Beeston JG, Nicolaides KH, Felton CV, Butler J, McGregor AM. Maturation of the secretion of thyroid hormone and thyroid-stimulating hormone in the fetus. N Engl J Med. 1991;324(8):532–6.
Weitzel JM, Radtke C, Seitz HJ. Two thyroid hormone-mediated gene expression patterns in vivo identified by cDNA expression arrays in rat. Nucleic Acids Res. 2001;29(24):5148–55.
Hiroi Y, Kim HH, Ying H, Furuya F, Huang Z, Simoncini T, et al. Rapid nongenomic actions of thyroid hormone. Proc Natl Acad Sci U S A. 2006;103(38):14104–9.
Wrutniak-Cabello C, Casas F, Cabello G. Thyroid hormone action in mitochondria. J Mol Endocrinol. 2001;26(1):67–77.
Casas F, Rochard P, Rodier A, Cassar-Malek I, Marchal-Victorion S, Wiesner RJ, et al. A variant form of the nuclear triiodothyronine receptor c-ErbAalpha1 plays a direct role in regulation of mitochondrial RNA synthesis. Mol Cell Biol. 1999;19(12):7913–24.
Janssen BG, Saenen ND, Roels HA, Madhloum N, Gyselaers W, Lefebvre W, et al. Fetal thyroid function and in utero exposure to fine particle air pollution: a birth cohort study. Environ Health Perspect. 2017;125(4):699–705.
Byun H-M, Panni T, Motta V, Hou L, Nordio F, Apostoli P, et al. Effects of airborne pollutants on mitochondrial DNA methylation. Part Fibre Toxicol. 2013;10(1):18.
Rebelo AP, Williams SL, Moraes CT. In vivo methylation of mtDNA reveals the dynamics of protein-mtDNA interactions. Nucleic Acids Res. 2009;37(20):6701–15.
Lee HC, Wei YH. Mitochondrial biogenesis and mitochondrial DNA maintenance of mammalian cells under oxidative stress. Int J Biochem Cell Biol. 2005;37(4):822–34.
Wang ND, Finegold MJ, Bradley A, Ou CN, Abdelsayed SV, Wilde MD, et al. Impaired energy homeostasis in C/EBP alpha knockout mice. Science. 1995;269(5227):1108–12.
Janssen BG, Byun H-M, Cox B, Gyselaers W, Izzi B, Baccarelli AA, et al. Variation of DNA methylation in candidate age-related targets on the mitochondrial-telomere axis in cord blood and placenta. Placenta. 2014;35(9):665–72.
Barrès R, Osler ME, Yan J, Rune A, Fritz T, Caidahl K, et al. Non-CpG methylation of the PGC-1α promoter through DNMT3B controls mitochondrial density. Cell Metab. 2009;10(3):189–98.
Izzi B, Francois I, Labarque V, Thys C, Wittevrongel C, Devriendt K, et al. Methylation defect in imprinted genes detected in patients with an Albright’s hereditary osteodystrophy like phenotype and platelet Gs hypofunction. PLoS ONE. 2012;7(6):e38579.
Xu J, Pope SD, Jazirehi AR, Attema JL, Papathanasiou P, Watts JA, et al. Pioneer factor interactions and unmethylated CpG dinucleotides mark silent tissue-specific enhancers in embryonic stem cells. Proc Natl Acad Sci U S A. 2007;104(30):12377–82.
Lambertini L, Byun H-M. Mitochondrial epigenetics and environmental exposure. Curr Envir Health Rpt. 2016;3(3):214–24.
Grigoriu A, Ferreira JC, Choufani S, Baczyk D, Kingdom J, Weksberg R. Cell specific patterns of methylation in the human placenta. Epigenetics. 2011;6(3):368–79.
Houseman EA, Accomando WP, Koestler DC, Christensen BC, Marsit CJ, Nelson HH, et al. DNA methylation arrays as surrogate measures of cell mixture distribution. BMC Bioinformatics. 2012;13:86.
Jaffe AE, Irizarry RA. Accounting for cellular heterogeneity is critical in epigenome-wide association studies. Genome Biol. 2014;15(2):1–9.
Michels KB, Binder AM, Dedeurwaerder S, Epstein CB, Greally JM, Gut I, et al. Recommendations for the design and analysis of epigenome-wide association studies. Nat Meth. 2013;10(10):949–55.
Janssen BG, Madhloum N, Gyselaers W, Bijnens E, Clemente DB, Cox B, et al. Cohort Profile: The ENVIRonmental influence ON early AGEing (ENVIRONAGE): a Birth Cohort Study. Int J Epidemiol. 2017; doi: 10.1093/ije/dyw269.
Janssen BG, Munters E, Pieters N, Smeets K, Cox B, Cuypers A, et al. Placental mitochondrial DNA content and particulate air pollution during in utero life. Environ Health Perspect. 2012;120(9):1346–52.
Hellemans J, Mortier G, De Paepe A, Speleman F, Vandesompele J. qBase relative quantification framework and software for management and automated analysis of real-time quantitative PCR data. Genome Biol. 2007;8(2):R19.
Valeri L, Vanderweele TJ. Mediation analysis allowing for exposure-mediator interactions and causal interpretation: theoretical assumptions and implementation with SAS and SPSS macros. Psychol Methods. 2013;18(2):137–50.
Gupta A, Srivastava S, Bhatnagar A. Cord blood thyroid stimulating hormone level—interpretation in light of perinatal factors. Indian Pediatr. 2014;51(1):32–6.
Lv P-P, Meng Y, Lv M, Feng C, Liu Y, Li J-Y, et al. Altered thyroid hormone profile in offspring after exposure to high estradiol environment during the first trimester of pregnancy: a cross-sectional study. BMC Med. 2014;12:240.
Soldin OP, Goughenour BE, Gilbert SZ, Landy HJ, Soldin SJ. Thyroid hormone levels associated with active and passive cigarette smoking. Thyroid. 2009;19(8):817–23.
Herbstman J, Apelberg BJ, Witter FR, Panny S, Goldman LR. Maternal, infant, and delivery factors associated with neonatal thyroid hormone status. Thyroid. 2008;18(1):67–76.
Chan LY-S, Leung TN, Lau TK. Influences of perinatal factors on cord blood thyroid-stimulating hormone level. Acta Obstet Gynecol Scand. 2001;80(11):1014–8.
The authors thank the participating mothers and neonates, as well as the staff of the maternity ward, midwives, and the staff of the clinical laboratory of East-Limburg Hospital in Genk.
The ENVIRONAGE birth cohort is supported by the European Research Council (ERC-2012-StG.310898), by the Flemish Scientific Fund (FWO, G.0.733.15.N), and the Special Research Fund (BOF) of Hasselt University. This work was also supported by funding from the National Institute of Environmental Health Sciences (R21ES022694 and R01ES021733).
Availability of data and materials
The datasets during and/or analyzed during the current study available from the corresponding author on reasonable request.
TSN coordinates the ENVIRONAGE birth cohort and designed the current study together with BGJ and AAB. WG and BGJ gave guidance to the midwives and did the quality control of the database. BGJ performed the experiments with the help of HMB and BGJ carried out the statistical analysis. JP performed the thyroid hormone analysis. BGJ, HMB, AAB, and TSN did the interpretation of the data. BGJ wrote the first draft of the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
All procedures were according to the principles outlined in the Helsinki Declaration for investigation of human subjects and approved by the Ethical Committee of Hasselt University and the East-Limburg Hospital. All participants provided written informed consent.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional figures. Figure S1. Correlation plot between MT-RNR1 and D-loop mtDNA methylation levels. The displayed methylation levels are absolute percentages. Figure S2. Estimated proportion of effects of FT4 exposure on mtDNA content mediated by mtDNA methylation. Figure S3. Flowchart depicting the selection for arriving at the final study sample either for placental mtDNA content measurements (n = 547) or for mtDNA methylation analysis (n = 305). (DOCX 520 kb)
Additional tables. Table S1. Associations of cord blood thyroid hormones FT3 and FT4 with CpG-specific mtDNA methylation and PPARGC1A promoter methylation. Table S2. Associations of fetal thyroid hormones with placental mtDNA methylation and mtDNA content while excluding women who had undergone a cesarean section (n = 13). Table S3. Associations of fetal thyroid hormones with placental mtDNA methylation and mtDNA content while excluding small for gestational age infants (n = 22). Table S4. Bisulfite-pyrosequencing primer sequence information based upon Assembly GRCh37/hg19 of the UCSC genome browser. Table S5. Mitochondrial and nuclear primer sequence information based upon Assembly GRCh37/hg19 of the UCSC genome browser. (DOCX 54 kb)