Expressions of mitochondria-related genes in pregnant women with subclinical hypothyroidism, and expressions of miRNAs in maternal and cord blood
Thyroid Research volume 16, Article number: 38 (2023)
Subclinical hypothyroidism in pregnancy and definition by upper thyrotropin (TSH) cutoff are controversial. As mitochondria are influenced by thyroid hormones, the purpose in this study was to measure expression of mitochondria-related genes in euthyroid and subclinical hypothyroid pregnant women to obtain more knowledge of potential metabolic consequences of maternal subclinical hypothyroidism. In addition, we wished to test if applied TSH-cutoff significantly changed our results of expressed gene-levels. Moreover, we aimed to identify potential microRNA-biomarkers for subclinical hypothyroidism – markers that could be traced to offspring as well.
From a cohort of at-term pregnant women undergoing planned cesarean section, 77 women had expression levels of the mitochondria-related genes Peroxisome Proliferator-activated Receptor-γ coactivator-1β (PGC-1β), mitochondrial Transcription Factor A (TFAM), Superoxide Dismutase 2 (SOD2) and Nuclear Respiratory Factor 2 (NRF-2) determined by qPCR from blood sampled in prior to delivery. Two TSH-cutoff levels defining subclinical hypothyroidism (> 3.0 and > 3.7 mIU/L) were applied for the procession of results, generating two data analyses of the same cohort. In 22 pairwise maternal-cord samples (subclinical hypothyroid/euthyroid-rate 0.5, TSH-cutoff > 3.0 mIU/L), microRNA-expressions (miRNA) were analyzed.
All gene expressions were lower in the subclinical hypothyroid group regardless of applied TSH-cutoff, but insignificant except for PGC-1β at TSH cutoff > 3.0 mIU/L. Two miRNAs (hsa-let-7d-3p and hsa-miR-345-5p) were upregulated in blood from women and offspring (cord blood) with subclinical hypothyroidism.
A trend towards decreased mitochondrial gene expressions in subclinical hypothyroidism were demonstrated. The miRNAs hsa-let-7d-3p and hsa-miR-345-5p might be potential markers of maternal subclinical hypothyroidism. However, larger studies are needed to verify the findings.
Subclinical hypothyroidism (SCH) is characterized by increased thyrotropin (TSH) and normal levels of the thyroid hormones thyroxine (T4) and triiodothyronine (T3). Thyroid hormones increase fat mobilization, stimulate insulin-dependent glucose uptake, gluconeogenesis and glycogenolysis . By genomic and non-genomic actions, thyroid hormones regulate important cellular and metabolic processes, including mitochondrial function [2, 3]. On the genomic (nuclear) level, thyroid hormones regulate thyroid hormone nuclear receptors and thyroid receptor transcription activity . These effects can be measured within hours to days . On the non-genomic level, thyroid hormones more directly regulate mitochondrial energetics, mitochondrial biogenesis, induction of mitochondrial DNA and non-genomic actions in mitochondria . These effects can be measured within seconds to minutes [2, 3]. Mitochondrial function can be evaluated by all these effects and not only by measuring energetics, which traditionally has been the most common approach. However, thyroid hormone regulation of mitochondria holds many unanswered and complex questions, and thus a detailed description of regulation mechanisms is yet to be detected .
Examples of thyroid regulated mitochondria-related transcripts (genomic level) are the co-activators from the PGC-1 family [5,6,7,8], the transcription factors Nuclear Respiratory Factor 1 (NRF-1) or 2 (NRF-2) [5, 9], mitochondrial Transcription Factor A (TFAM) [5, 6] and the antioxidant enzymes Superoxide Dismutase 1–3 (SOD1-3) [10, 11]. Expression of these mitochondria-related genes are all linked to mitochondrial biogenesis and function, and therefore, thyroid dysfunction may be associated with changes in these. Even subtle thyroid changes within normal reference values may lead to mitochondrial changes . For instance, Kristensen et al. demonstrated that hemi-thyroidectomy significantly increased TSH (within normal range), and concomitantly influenced mitochondrial energetics in terms of increasing the mitochondrial membrane potential in mononuclear blood cells . Also, the mitochondria-related genes peroxisome proliferator-activated receptor-γ coactivator-1α (PGC-1α) and superoxide dismutase 2 (SOD2) were downregulated . Findings like this indicate a necessity to gain more knowledge of subclinical hypothyroid state and the association with cellular metabolism.
The number of mitochondria vary, dependent on tissue sort and energy demands [10, 15]. Likewise, thyroid hormone levels are regulated in different tissues by the expressions of different deiodinases. In placental and fetal tissue, the type 3 deiodinase (D3) is predominantly expressed and inactivates thyroid hormones, while type 2 deiodinase (D2) is also expressed in the placenta, where it deiodinates T4 to the active hormone T3 . These mechanisms make pregnancy and thyroid function very interesting, as thyroid hormone exposure to the fetus seems very adequately tuned. As thyroid hormones regulate mitochondria, it is therefore worth taking a glance towards cellular processes concerning the “pregnancy load” on the thyroid function, as thyroid”normality” in pregnancy has been difficult to define and to understand the impact of.
In short, recommendations for the definition of subclinical hypothyroidism in pregnancy have in recent years been eagerly debated and changed from in 2011 (by the American Thyroid Association) and in 2012 (by the Endocrine Society) a definition in the third trimester of TSH > 3.0 mIU/L [16, 17], over TSH > 3.5 mIU/L in 2014 (by the European Thyroid Association) , to in 2017 local established references, with a maximum of 4.0 mIU/L (by the American Thyroid Association) . Moreover, the clinical importance and the need for treatment of subclinical hypothyroidism have been discussed . Studies are difficult to compare, as TSH levels are dependent on factors such as ethnicity, iodine supply, age, laboratory methods  and thyroid autoimmunity  as well as goitrogens, for instance from smoking .
In Denmark, the Danish Society of Endocrinologists in 2018 recommended to use TSH > 3.7mIU/L as the upper cutoff for the definition of subclinical hypothyroidism in the third trimester, in the absence of local laboratory pregnancy references . However, this threshold was recently changed to TSH > 3.5 mIU/L  which underlines that definitions and clinicians´ perception are persistently changing, also within populations. In relation to the controversial debate concerning the clinical significance of sub-hypothyroid state in pregnancy, it seems increasingly important to gain more knowledge about thyroid impact, also on the cellular level.
To our knowledge, few studies have investigated mitochondrial changes in subclinical thyroid disorder in pregnancy . As subclinical hypothyroidism may be rather common in pregnancy [12, 26], it seems relevant to assess thyroid regulated mitochondria-related gene expressions in subclinical hypothyroid pregnant women. Optimizing our understanding of metabolic effects, and thereby thyroid normality in pregnancy, could contribute to a more reliable definition.
Other regulators of mitochondria-related genes and of thyroid function are microRNAs (miRNA). A miRNA is a 19–25 nucleotides long non-coding RNA-strand that regulate gene expression at the posttranscriptional level, by binding to messenger RNA (mRNA) . In this way, they either suppress mRNA translation or induce mRNA degradation . Each miRNA can target several mRNAs, as one mRNA can be the target of several miRNAs [27, 28], and in this way miRNAs regulate both normal and pathological processes  in cell proliferation, differentiation, apoptosis and in metabolism . Concerning the thyroid, miRNA research has primarily been centered around cancer. However, evidence exists that miRNAs are essential for thyroid gland development, thyroid hormone production, and for the regulation of locally expressed deiodinases in different tissues which adjust the bioavailability of thyroid hormones, and thereby their effects . Knowledge of miRNA-expressions related to subclinical hypothyroidism is sparse.
Our hypothesis for this study was that a mitochondrial dysfunction could be present in subclinical hypothyroid pregnant women, reflected by a downregulated level of mitochondria related gene transcripts. Moreover, that a potential mitochondrial dysfunction could be traced in cord blood at delivery by biomarkers (miRNA) that would be present in the mothers, too.
Aims of the present study were to
Compare mitochondrial biogenesis in euthyroid and subclinical hypothyroid pregnant women by assessing levels of TFAM, PGC-1β, NRF-2 and SOD2.
Investigate if changing TSH-cutoff markedly affects results of gene-expression levels of TFAM, PGC-1β, NRF-1 and SOD2 in euthyroid vs. subclinical hypothyroid pregnant women.
Compare mean levels of expressed miRNAs in euthyroid and subclinical hypothyroid pregnant women, and their offspring, and assess a potential relation to mitochondrial biogenesis.
In the present study (Fig. 1), medical journals for all planned cesarean sections were in advance reviewed to define exclusion and inclusion criteria. Afterwards, third trimester pregnant women were recruited by an interview and blood samples prior to a planned cesarean section. After delivery, their offspring were enrolled by cord blood samples. The women were included at Obstetrics Dept., Naestved Hospital, Denmark from January – February 2014 (pilot study) and again from June 2014—July 2015 (main study period).
Inclusion criteria: maternal age of at least eighteen years, an expected healthy, singleton pregnancy verified by a routine ultrasound scan in the second trimester and maternal health without any diseases or medication causing adverse effects to the fetus.
Exclusion criteria: multiple pregnancies, medical treated thyroid disease, diabetes (including gestational), preeclampsia or any other disease suspected to cause adverse effect to the fetus. Lack of informed verbal and written consent from both parents led to exclusion, too.
In the study period, a total of 180 planned cesarean sections were performed. Sampling logistics and laboratory capacity (n = 31), challenged family situation and missing full parental consent (n = 14) as well as positive exclusion criteria (n = 39) were the cause that only 96 women were invited to participate. Of these, 68 accepted (70.83%). In the pilot study, 9 were included, obtaining a total of 77 participating women. According to guidelines at inclusion time (TSH > 3.0 mIU/L) , 24 women fulfilled the criteria for SCH, and 13 according to recent Danish guidelines (TSH > 3.7 mIU/L) [19, 24]. The rest were euthyroid. According to the different guidelines, study results were analyzed as cohorts SCH1 (subclinical hypothyroid group 1) and EU1 (euthyroid group 1) for TSH cutoff 3.0 mIU/L and SCH2 (subclinical hypothyroid group 2) and EU2 (euthyroid group 2) for TSH cutoff 3.7 mIU/L.
Maternal and cord metabolic analyses
A 3.5 mL heparinized blood sample was collected for same-day analysis of thyroid status, glucose, cholesterol and triglycerides. A serum gel tube (3.5 mL) was collected for thyroid peroxidase antibodies (anti-TPO) measurements.
Thyrotropin (TSH), free triiodthyronin (fT3) and free thyroxine (fT4) were analyzed by Siemens Dimension Vista System by an electrochemical luminescent immunoassay based on LOCI-technology (2008 Siemens Healthcare Diagnostics). Third trimester normal values for maternal TSH was defined as TSH: 0.3—3.0 mIU/L (EU1) or 0.3—3.7 mIU/L (EU2). The lowest value followed the standard used at the Dept. of Biochemistry at Naestved Hospital, and the upper value was in accordance with the different guidelines [16, 24]. Maternal normal values for fT3 and fT4 followed the references used at Naestved Hospital: fT4 = 8.5 – 26.0 pmol/L and fT3 = 2.7—6.1 pmol/L. Anti-TPO was measured on Kryptor by TRACE-technology (2005 Brahms Kryptor) with a detection limit of 11 kU/l. An anti-TPO > 60 kU/l was considered a positive test. To compare metabolism, glucose, total-cholesterol, low-density lipoprotein cholesterol (LDL), high-density lipoprotein cholesterol (HDL) and triglycerides were analyzed by Siemens Dimension Vista System (2008 Siemens Healthcare Diagnostics) by photometry.
Mitochondria-related RNA sequencing
For reverse transcription by quantitative Polymerase Chain Reaction (qPCR), a 3 mL Tempus™ Blood RNA Tube (Life Technologies, Denmark) containing a reagent that instantly stabilizes intracellular RNA, was used. Samples were subsequently frozen and stored at -80 °C until analysis. The laboratory technician was blinded to which thyroid groups the samples belonged. Sample collection and handling were performed according to the manufacturer´s recommendation.
RNA was purified using PerfectPure™ RNA Blood kit (5 Prime, AH Diagnostics) in accordance with the manufacturer’s recommendations. The integrity of RNA was characterized by the RNA integrity number (RIN) measured on an Agilent 2100 Bioanalyzer (Agilent Technologies, Copenhagen, Denmark). For analysis, a RIN value above 7.0 was considered acceptable. Whole blood gene expressions of PGC-1β, TFAM, NRF-2 and SOD2 were examined. Table 1 presents the primer sequences used.
Reverse transcription was performed using Quanti-Tect Reverse Transcription Kit (Qiagen, Copenhagen, Denmark). A SYBR Green based qRT-PCR of PGC-1β, TFAM, NRF-2 and SOD2 mRNA was measured on a LightCycler 480 (Roche). To correct for potential variation in cDNA loading and quantity, the measured gene transcripts levels were normalized to the expression of ribosomal 18 S RNA. The qRT-PCR reactions were performed according to the LightCycler standard protocol. The qPCR reaction mixture included 1 × LightCycler FastStart DNA Master PLUS SYBR Green I (Roche, Copenhagen, Denmark) and 1 μl mRNA specific- and 0.5 μl of the 18 S primers, 1 μmol/l respectively. For each reaction 5 μl template cDNA and sterile water were added in a total reaction volume of 20 μl. Cycling conditions were 95 °C for 10 min, followed by 45 cycles at 95 °C for 10 s, and 60 °C for 1 min. All quantitative RT-PCR measurements were performed in triplicates. Melting curves were completed for the control of unspecific DNA amplification after each run. Unit of measurement for gene expressions was arbitrary units (a.u.).
Some results did not reach acceptable quality, and were excluded. The excluded samples were from 2 women with TSH > 4 mIU/L, and from 2 with TSH < 2 mIU/L, BMI-range 25.1–27.6. Therefore, only 73 samples presented full gene profile results (EU1 n = 51, SCH1 n = 22, EU2 n = 62, SCH2 n = 11), and for TFAM one result could be added to euthyroid-groups (EU1 n = 52, EU2 n = 63).
The amount obtainable of cord blood were too low for qPCR analysis. This was due to a very sparse amount of blood left for the Tempus tube after collection of blood for standard pH-measurements (2–4 mL), for thyroid function-tests (3.5 mL), anti-TPO measurements (3.5 mL), for flow cytometry (4 mL) (results published in [26, 32]) and for miRNA analyses (4 mL). Seldom, blood was left in the cord after collection for these other tubes. If more blood was left, it was very difficult to do more sampling due to coagulation.
Maternal and cord micro-RNA panel
For miRNA-screening, only a minor number of participants were chosen due to cost of analyses. Twenty-two maternal samples (11 euthyroid, 11 subclinical hypothyroid) and the associated cord samples were chosen randomly, as they fulfilled quality criteria concerning sufficient blood volume and no signs of hemolysis by the manufacturer´s primary quality check.
One EDTA-plasma fraction was sampled for miRNA-analysis and subsequently frozen and stored at -80 °C until shipping for analysis.
The real-time PCR panel analysis was performed by Exiqon the following way: Each RNA sample was successfully reverse transcribed (RT) into cDNA and run on the miRCURY LNA™ Universal RT microRNA PCR Human Panel I + II (752 assays). For normalization of data, the average of the assays detected in all samples was applied as it was the most stable normalizer. Numbers of miRNA present were detected and the quantification cycle (Cq) value of the global mean for each of the samples was identified.
Each individual amplification product on PCR panels was scrutinized by melting curve analysis, calculation of amplification efficiency and comparison of Cq value to background level in a negative control sample. For quality control, two types of RNA-spike-ins kits (Exiqon) were used: For RNA isolation control, UniSp2, UniSp4, UniSp5 (Exiqon) were added to the purification to detect any differences in extraction efficiency. For cDNA synthesis control, UniSp6 (Exiqon) was added in the RT reaction, giving the opportunity to evaluate it. In addition, a DNA spike-in (UniSp3, Exiqon) was present on all panels, to indicate inhibitions at the qPCR level by deviations in this reaction.
Exiqon was blinded to which thyroid groups the samples belonged.
Two-tailed t-test was used to compare data with a normal distribution. Otherwise, results were handled non-parametrically (non-normality confirmed by Shapiro-Wilks test). Wilcoxon Rank-sum test was used to compare level of differences. Spearman´s rho ρ was used to access correlations. A p-value < 0.050 was considered statistically significant. Data processes were carried out using STATA version 15.0 for Windows (StataCorp LP, College station, TX, USA).
Regarding miRNA, Exiqon analysed the data as follows: A Principal Component Analysis (PCA) was used to reduce the dimension of the large data set and identify the miRNAs with the largest variation across samples. The most differentially expressed miRNAs were analysed by t-test and a Benjamini–Hochberg correction. The analyses were made in accordance with the definition of SCH at that time (TSH > 3.0 mIU/L).
A calculation with 90% strength at the 5% significance level (alpha) specified recruitment of 242 euthyroid and 48 subclinical hypothyroid women to show a significant decrease in subclinical hypothyroid group TFAM. Calculations were based on a mean value of 1.70 in euthyroid and 1.19 in subclinical hypothyroid group, and a common SD of 1.00. For PGC-1β, NFR-2 and SOD2, recruitment of 252 euthyroid and 50 subclinical hypothyroid women was necessary to demonstrate a significant difference. Calculations were based on a mean value of 4.70 in the euthyroid and 3.80 in the subclinical hypothyroid group, and a common SD of 2.00. Calculations were based on a rate of subclinical hypothyroidism /euthyroidism = 0.20 (TSH cutoff 3.7 mU/L).
Calculations were performed retrospectively, as this study was performed in supplement to a flow cytometric evaluation of mitochondrial function (primary outcome) of the same cohort (unpublished results ), preceded by a method study . Unfortunately, it was not possible to obtain the calculated power of the study, due to health circumstances in the research team, which interrupted the main study prematurely.
Maternal metabolic profiles
Table 2 depicts relevant maternal data. Age, BMI and gestational age were comparable between cohorts. In general, other metabolic parameters such as blood glucose level, total cholesterol, LDL, HDL and triglycerides did not differ (Table 2). Thus, other metabolic differences that could affect our results on the mitochondrial level were not present.
Mitochondria-related gene expression in relation to thyroid status
Mitochondria-related gene expression levels in euthyroid and subclinical hypothyroid pregnant women are depicted by the different reference standards of TSH in Fig. 2A (TSH cutoff 3.0 mIU/L) and Fig. 2B (TSH cutoff 3.7 mIU/L).
Irrespective of TSH cutoff, every gene showed a lower expression level in pregnant women with subclinical hypothyroidism compared with euthyroids. However, apart from PGC-1β expression in group EU1 vs. SCH1 (p = 0.033), results were non-significant (Fig. 2A-B). Notably, the mitochondria-related genes were positively correlated with each other (Table 3).
In supplementary section, maternal TSH is plotted against the different gene expressions (Fig. S1).
In the miRNA-study, only the TSH-cutoff of 3.0 mIU/L was used to define euthyroidism vs. subclinical hypothyroidism, as the analyses were made by Exiqon according to the definition of subclinical hypothyroidism at that time . Apart from maternal TSH (p < 0.001), thyroid profiles did not differ between euthyroid and subclinical hypothyroid women or offspring (Fig. 3A-B).
Mean maternal age did not differ significantly, and in both cohorts, the women delivered at mean gestational age of 39 + 1 (Table 4). The women with subclinical hypothyroidism were leaner, but not significant. Metabolic profile differed a bit from the profile in the genomic study, but there were no significant differences in blood glucose or lipids of euthyroid and subclinical hypothyroid women (Table 4).
Total cholesterol in cords varied from < 1.3 to 1.8 mmol/l, LDL from unmeasurable to 0.9 mmol/l, HDL from 0.4 to 0.8 mmol/l and triglycerids from < 0.1 to 0.3 mmol/l with no difference between groups (data not presented). Birthweights were very alike in the two groups, and rate of gender in the two groups were not significantly different (Table 4).
In all of the samples, 219 identical miRNAs could be presented. In each sample, an average of 407 miRNAs was detectable. When comparing maternal subclinical hypothyroid and euthyroid group, 30 miRNAs were differentially expressed by t-test, but none significantly after performing the Benjamini-Hochberg (BH) correction. For cord samples, 29 miRNAs were differentially expressed by t-test, and one persisted to be significant after the BH-correction. The 20 most differentially expressed miRNAs in maternal respective cord samples are depicted in Tables 5 and 6.
In other words, we were able to demonstrate biomarker (miRNA) differences related to TSH level by t-test. However, after BH-adjustment, only one (hsa-miR-340-3p) remained significant. Interestingly, two miRNAs (hsa-let-7d-3p and hsa-miR-345-5p) were equally upregulated in subclinical hypothyroid maternal and cord blood, which suggest they may be markers of SCH.
Mitochondria-related gene expression in relation to thyroid status
The results of mitochondria-related gene expressions in relation to thyroid status presented here were insignificant, except for PGC-1β expression at TSH-cutoff 3.0 mIU/L. The insignificant findings in this study are most likely caused by the too small number of participants, which increases the risk of type-II errors. However, as the decrease in absolute terms in gene expressions in subclinical hypothyroid cohort was persistent—irrespective of TSH-cutoff applied—we suggest that a larger study could be appropriate, to investigate these very interesting findings more thoroughly.
To support the theory that subclinical hypothyroidism is associated with decreased expression of mitochondria-related genes, other studies have demonstrated noteworthy changes by overt hypothyroidism. Sinha et al. demonstrated decreased expression of the mitochondria-related genes NRF-1 α, PGC-1α and TFAM in brain tissue of hypothyroid rodents , and Sagliocchi et al. demonstrated decreased SOD2 levels in hypothyroid muscles compared to control muscles, and oppositely increased SOD2 levels in hyperthyroid muscles . Overall, these are in accordance with our findings.
Another noteworthy finding in the present study is the demonstration of positive correlations between the investigated mitochondria-related genes, confirming they are part of a common regulatory pathway. To our knowledge, only few studies have investigated the correlation between regulators of mitochondrial function , but Fabricius et al. had similar results with positive correlations between expressed mitochondria-related genes within specific tissue types . They investigated expressions in smooth muscle tissue (uterus), white adipose tissue and blood of 17 persons and found different levels of expressions in different tissues. PGC-1β and TFAM were most pronouncedly expressed in smooth muscle tissue, while NRF-1 and NRF-2 were expressed lowest in adipose tissue. In blood, all genes had a low, but similar expression  which makes it ideal for comparative studies where the objective is to detect small changes in expression between groups. This, and the factor that blood is easy obtainable, have led us to use blood in the present study.
PGC-1β was chosen as it is a co-activator that respond rapidly and robust to thyroid hormones by binding to different transcription factors that induce mitochondrial biogenesis  and energy metabolism in terms of regulation of beta-oxidation of fatty acids and oxidative phosphorylation in mitochondria . Moreover, it regulates targets of mitochondrial fusion and fission, which is fundamental for mitochondrial repair and cell fate . Functionally, it appears related to PGC-1α but are probably more energy inducing as its activation leads to a higher level of coupled respiration .
NRF-2 is a DNA-binding transcription factor implicated in activating cytochrome oxidase expression, and controlling nucleus-encoded subunits of cytochrome oxidase . Together, NRF-2 and NRF-1 influence the expression of nuclear genes related to mitochondrial function . It is well known that NRF-1 is a target of regulation for every member of the PGC-1 family, but studies in mouse tissue and cultured cells suggest NRF-2 as a target, too . For that reason, it seemed important to choose NRF-2 in the present study.
SOD2 is involved in intracellular response to oxidative stress. In most cells, it is present in abundant amounts where it neutralizes reactive oxygen species (ROS), by accepting unpaired electrons [10, 11]. SOD2 differs from the other Superoxide Dismutases by its localization inside the mitochondrial matrix .
By the choice of transcription factors that regulate different compartments related to mitochondria, and a co-activator that regulate at least one of these, we hoped to demonstrate a common trend of changes in response to increasing TSH, and in absolute terms we succeeded in this.
In the present study, we demonstrated a general variety of expressed miRNAs in maternal vs. cord samples. These different expressions were most likely caused by the different composition of maternal vs. cord blood. Cord blood contain stem cells, and another metabolism is present in cord blood compared to adult blood . This different metabolism was reflected by the very different levels of thyroid hormones and lipid status in cord blood.
We demonstrated a tendency of different miRNA-expressions related to thyroid status, and perhaps identified a potential biomarker for maternal subclinical hypothyroidism (hsa-miR-340-3p). Interestingly, two potential biomarkers for subclinical hypothyroid disease in the woman and in her offspring were identified (hsa-let-7d-3p and hsa-miR-345-5p), implicating a potential effect on the offspring to maternal subclinical hypothyroid disease. In a recent overview of miRNAs linked to autoimmune thyroid disease, none of these miRNAs were described . On the contrary, one of those isoforms (hsa-let7d) has been described to be highly expressed by normal functioning thyroid glands , which is why this upregulation in subclinical hypothyroid state is very interesting. However, Massolt et al.  investigated the expression of a panel of 384 selected miRNAs in different thyroid states of 13 thyroidectomized patients in treatment due to differentiated thyroid cancer. Some had suppressed TSH, and some had a very high TSH. In contrast to our study, they did not find any differences in expression of miRNAs in relation to thyroid status . This could be caused by their smaller number of study participants, or perhaps by choice of other markers. Larger studies of miRNAs as biomarkers for changed metabolism in relation to thyroid disease would be appropriate to gain more knowledge about this field.
We searched the microPIR , a database that contains information about promoter interactions. The miRNAs hsa-let-7d-3p and hsa-miR-345-5p were not identified in it. However, every miRNA with accessible information had the specific mitochondria-related genes, the current paper assesses, among its targets (or the “family members” NRF-1, SOD1, SOD3), except from hsa-miR-551a. We also searched the TargetScanHuman database  and the miRBase  for a relation between the most noteworthy miRNAs in this study (hsa-let-7d-3p, hsa-miR-345-5p, hsa-miR-340-3p) and interactions with the examinated gene expressions in this study, but with no positive results.
In future studies, it could be beneficial to use miRNA mimics or inhibitors in cellular models to test if the identified miRNAs are directly associated with subclinical hypothyroidism, or if their effects are merely through regulation of gene targets already known.
In the microPIR database , PGC-1β was not listed as a target of regulation by miRNAs. Instead, PGC-1β is known to be a coactivator and a regulator of miRNAs . This is in accordance with a study of thyroid tumor cells, with a high content of functional mitochondria, where evidence was found that coactivators of the PGC-1 family control mitochondrial function by regulation of miRNA expression .
Strengths to the study were comparable study groups in terms of gestational age and delivery method as well as metabolic profiles. Moreover, in the miRNA-study, the ratio of offspring gender was acceptable, factors that otherwise may influence miRNA expression .
Weaknesses were primarily lack of power, but also supplementary analyzes of cord blood quality could have been performed locally . However, Exiqon is experienced in miRNA services, and therefore the plasma quality tests were run by their standards. Regarding collection and freezing, the samples were handled in a standardized way. However, variation in time delay from delivery to possibility of sampling after the cord clamping, could not be avoided, even though the standardization should minimize this aspect.
Minor weaknesses were the increased rate of smoking in the subclinically hypothyroid group compared to the euthyroid group in the miRNA-study, as smoking might modulate expressions. However, smoking in pregnancy is associated with a lower TSH than in non-smokers , which we have also demonstrated in our material .
To evaluate the overall interaction between maternal and fetal environment, it would have been beneficial (had the study been larger) to add measurements of gene- and miRNA expressions in placental tissue, and gene expressions in cord blood, too.
Although non-significant, this study demonstrates a decrease in absolute terms in expression-level of the mitochondria-related genes PGC-1β, TFAM, NRF-2 and SOD2 in third trimester pregnant women with subclinical hypothyroidism. Two potential miRNA markers of subclinical hypothyroidism in maternal and cord blood (offspring) and one potential maternal miRNA marker were identified, and a variety of miRNAs with mitochondria-regulating genes as targets had different expressions in subclinical hypothyroid vs. euthyroid maternal samples. This could indicate changes on the cellular and metabolic level by increased TSH, even though the measured levels of thyroid hormones are considered normal in subclinical hypothyroidism. Therefore, this study may facilitate new research on thyroid regulation of mitochondria on the cellular level and in relation to heritage, and we encourage further research due to the small sample size in our study. Expression of metabolic biomarkers could be helpful in establishing the correct cutoff for subclinical hypothyroid disease in pregnancy, as they could help to understand normal and pathological processes and thereby determine the biological normal cutoff. In this way, clinical decision making of “to treat or not to treat” a patient could be supported.
Availability of data and materials
The data that support the findings of this study are available from the main author, but restrictions apply to the availability of these data, which were used under license for the current study, and so are not publicly available. Data are however available from the authors upon reasonable request and with permission of the Danish Data Protection Agency. The Exicon Project Report, describing miRNA data analysis, is also available on request.
Quantitative Polymerase Chain Reaction
Peroxisome Proliferator-activated Receptor-γ coactivator-1β
Mitochondrial Transcription Factor A
Superoxide Dismutase 2
Nuclear Respiratory Factor 2
Cicatiello AG, Di Girolamo D, Dentice M. Metabolic effects of the intracellular regulation of thyroid hormone: old players, new concepts. Front Endocrinol (Lausanne). 2018;9:474. https://doi.org/10.3389/fendo.2018.00474.
Cioffi F, Senese R, Lanni A, Goglia F. Thyroid hormones and mitochondria: with a brief look at derivatives and analogues. Mol Cell Endocrinol. 2013;379(1–2):51–61. https://doi.org/10.1016/j.mce.2013.06.006.
Cioffi F, Giacco A, Goglia F, Silvestri E. Bioenergetic aspects of mitochondrial actions of thyroid hormones. Cells. 2022;11(6). https://doi.org/10.3390/cells11060997.
Cheng SY, Leonard JL, Davis PJ. Molecular aspects of thyroid hormone actions. Endocr Rev. 2010;31(2):139–70. https://doi.org/10.1210/er.2009-0007.
Scarpulla RC. Metabolic control of mitochondrial biogenesis through the PGC-1 family regulatory network. Biochim Biophys Acta. 2011;1813(7):1269–78. https://doi.org/10.1016/j.bbamcr.2010.09.019.
Fabricius MH, Wilms LK, Larsen J, Pedersen PL, Anthonsen S, Kvetny J. Measure of expression of mitochondrial related genes in human mononuclear blood cells, adipose white tissue and smooth muscle cells. Clin Chim Acta. 2010;411(9–10):749–53. https://doi.org/10.1016/j.cca.2010.02.005.
Shao D, Liu Y, Liu X, Zhu L, Cui Y, Cui A, et al. PGC-1 beta-regulated mitochondrial biogenesis and function in myotubes is mediated by NRF-1 and ERR alpha. Mitochondrion. 2010;10(5):516–27. https://doi.org/10.1016/j.mito.2010.05.012.
Weitzel JM, Iwen KA. Coordination of mitochondrial biogenesis by thyroid hormone. Mol Cell Endocrinol. 2011;342(1–2):1–7. https://doi.org/10.1016/j.mce.2011.05.009.
Klinge CM. Estrogens regulate life and death in mitochondria. J Bioenerg Biomembr. 2017;49(4):307–24. https://doi.org/10.1007/s10863-017-9704-1.
Malik AN, Czajka A. Is mitochondrial DNA content a potential biomarker of mitochondrial dysfunction? Mitochondrion. 2013;13(5):481–92. https://doi.org/10.1016/j.mito.2012.10.011.
Flynn JM, Melov S. SOD2 in mitochondrial dysfunction and neurodegeneration. Free Radic Biol Med. 2013;62:4–12. https://doi.org/10.1016/j.freeradbiomed.2013.05.027.
Feldthusen AD, Larsen J, Pedersen PL, Toft Kristensen T, Kvetny J. Pregnancy-induced alterations in mitochondrial function in euthyroid pregnant women and pregnant women with subclinical hypothyroidism; relation to adverse outcome. J Clin Transl Endocrinol. 2014;1(1):e13–7. https://doi.org/10.1016/j.jcte.2013.12.003.
Toft Kristensen T, Larsen J, Pedersen PL, Feldthusen AD, Ellervik C, Jelstrup S, et al. Persistent cellular metabolic changes after hemithyroidectomy for benign euthyroid goiter. Eur Thyroid J. 2014;3(1):10–6. https://doi.org/10.1159/000357943.
Kristensen TT, Pedersen PL, Larsen J, Feldthusen AD, Jelstrup S, Ellervik C. Hormone and metabolic research = Hormon- und Stoffwechselforschung = Hormones et metabolisme. 2019;51(2):127–33. https://doi.org/10.1055/a-0822-3066.
Harper ME, Seifert EL. Thyroid hormone effects on mitochondrial energetics. Thyroid. 2008;18(2):145–56. https://doi.org/10.1089/thy.2007.0250.
De Groot L, Abalovich M, Alexander EK, Amino N, Barbour L, Cobin RH, et al. Management of thyroid dysfunction during pregnancy and postpartum: an Endocrine Society clinical practice guideline. J Clin Endocrinol Metab. 2012;97(8):2543–65. https://doi.org/10.1210/jc.2011-2803.
Stagnaro-Green A, Abalovich M, Alexander E, Azizi F, Mestman J, Negro R, et al. Guidelines of the American Thyroid Association for the diagnosis and management of thyroid disease during pregnancy and postpartum. Thyroid. 2011;21(10):1081–125. https://doi.org/10.1089/thy.2011.0087.
Lazarus J, Brown RS, Daumerie C, Hubalewska-Dydejczyk A, Negro R, Vaidya B. 2014 European thyroid association guidelines for the management of subclinical hypothyroidism in pregnancy and in children. Eur Thyroid J. 2014;3(2):76–94. https://doi.org/10.1159/000362597.
Alexander EK, Pearce EN, Brent GA, Brown RS, Chen H, Dosiou C, et al. 2017 Guidelines of the American Thyroid Association for the diagnosis and management of thyroid disease during pregnancy and the postpartum. Thyroid. 2017;27(3):315–89. https://doi.org/10.1089/thy.2016.0457.
Tsakiridis I, Giouleka S, Kourtis A, Mamopoulos A, Athanasiadis A, Dagklis T. Thyroid disease in pregnancy: a descriptive review of guidelines. Obstet Gynecol Surv. 2022;77(1):45–62. https://doi.org/10.1097/ogx.0000000000000960.
Ittermann T, Khattak RM, Nauck M, Cordova CM, Volzke H. Shift of the TSH reference range with improved iodine supply in Northeast Germany. Eur J Endocrinol. 2015;172(3):261–7. https://doi.org/10.1530/EJE-14-0898.
Bliddal S, Derakhshan A, Xiao Y, Chen LM, Männistö T, Ashoor G, et al. Association of thyroid peroxidase antibodies and thyroglobulin antibodies with thyroid function in pregnancy: an individual participant data meta-analysis. Thyroid. 2022;32(7):828–40. https://doi.org/10.1089/thy.2022.0083.
Laurberg P, Cerqueira C, Ovesen L, Rasmussen LB, Perrild H, Andersen S, et al. Iodine intake as a determinant of thyroid disorders in populations. Best Pract Res Clin Endocrinol Metab. 2010;24(1):13–27. https://doi.org/10.1016/j.beem.2009.08.013.
Juul AB, et al. Society of Danish endocrinologists, National guideline:Thyroid disease in pregnancy and infertility. Was provided at: https://www.endocrinologydk/nbv/thyroideasygdomme/thyroideasygdom-ved-graviditet-og-infertiliet/. 2018.
Grovesen DL PJ, Pedersen IB, Nielsen HS, Feldt-Rasmussen U, Nielsen T, Nygård B, Juul AB, Boas M, Biddal S, Brix T. Society of Danish endocrinologists, National guideline:Thyroid disease in pregnancy and infertility. Now provided at: https://www.endocrinologydk/nbv/thyroideasygdomme/thyroideasygdom-ved-graviditet-og-infertiliet/. 2022
Stryhn JKG, Larsen J, Pedersen PL, Gæde PH. Subclinical hypothyroidism in pregnancy – assessment of offspring thyroid status and mitochondrial robustness to stress. Re-submitted to Scandinavian Journal of Clinical and Laboratory Investigation. 2023.
Boufraqech M, Klubo-Gwiezdzinska J, Kebebew E. MicroRNAs in the thyroid. Best Pract Res Clin Endocrinol Metab. 2016;30(5):603–19. https://doi.org/10.1016/j.beem.2016.10.001.
Ramirez-Moya J, Santisteban P. miRNA-Directed Regulation of the Main Signaling Pathways in Thyroid Cancer. Front Endocrinol (Lausanne). 2019;10:430. https://doi.org/10.3389/fendo.2019.00430.
Lizarraga D, Huen K, Combs M, Escudero-Fung M, Eskenazi B, Holland N. miRNAs differentially expressed by next-generation sequencing in cord blood buffy coat samples of boys and girls. Epigenomics. 2016;8(12):1619–35. https://doi.org/10.2217/epi-2016-0031.
Morales-Prieto DM, Ospina-Prieto S, Chaiwangyen W, Schoenleben M, Markert UR. Pregnancy-associated miRNA-clusters. J Reprod Immunol. 2013;97(1):51–61. https://doi.org/10.1016/j.jri.2012.11.001.
Aranda A. MicroRNAs and thyroid hormone action. Mol Cell Endocrinol. 2021;525:111175. https://doi.org/10.1016/j.mce.2021.111175.
Stryhn JKG, Larsen J, Pedersen PL, Feldthusen AD, Kvetny J, Gaede PH. Mitochondrial energetics and contents evaluated by flow cytometry in human maternal and umbilical cord blood. Scand J Clin Lab Invest. 2020;80(5):351–9. https://doi.org/10.1080/00365513.2020.1768584.
Campbell CT, Kolesar JE, Kaufman BA. Mitochondrial transcription factor A regulates mitochondrial transcription initiation, DNA packaging, and genome copy number. Biochim Biophys Acta. 2012;1819(9–10):921–9. https://doi.org/10.1016/j.bbagrm.2012.03.002.
Kristensen TT, Pedersen PL, Larsen J, Feldthusen AD, Jelstrup S, Ellervik C. Reduced Gene Expression of Peroxisome Proliferator-Activated Receptor Gamma Coactivator-1α in Whole Blood in Euthyroid Patients One Year After Hemithyroidectomy for Benign Euthyroid Goiter. Horm Metab Res. 2019;51(2):127–33. https://doi.org/10.1055/a-0822-3066.
Piriyapongsa J, Bootchai C, Ngamphiw C, Tongsima S. microPIR: an integrated database of microRNA target sites within human promoter sequences. PLoS One. 2012;7(3):e33888. https://doi.org/10.1371/journal.pone.0033888.
Sinha RA, Pathak A, Mohan V, Babu S, Pal A, Khare D, et al. Evidence of a bigenomic regulation of mitochondrial gene expression by thyroid hormone during rat brain development. Biochem Biophys Res Commun. 2010;397(3):548–52. https://doi.org/10.1016/j.bbrc.2010.05.154.
Sagliocchi S, Cicatiello AG, Di Cicco E, Ambrosio R, Miro C, Di Girolamo D, et al. The thyroid hormone activating enzyme, type 2 deiodinase, induces myogenic differentiation by regulating mitochondrial metabolism and reducing oxidative stress. Redox Biol. 2019;24:101228. https://doi.org/10.1016/j.redox.2019.101228.
Dumont S, Le Pennec S, Donnart A, Teusan R, Steenman M, Chevalier C, et al. Transcriptional orchestration of mitochondrial homeostasis in a cellular model of PGC-1-related coactivator-dependent thyroid tumor. Oncotarget. 2018;9(22):15883–94. https://doi.org/10.18632/oncotarget.24633.
Borengasser SJ, Faske J, Kang P, Blackburn ML, Badger TM, Shankar K. In utero exposure to prepregnancy maternal obesity and postweaning high-fat diet impair regulators of mitochondrial dynamics in rat placenta and offspring. Physiol Genomics. 2014;46(23):841–50. https://doi.org/10.1152/physiolgenomics.00059.2014.
Martínez-Hernández R, Marazuela M. MicroRNAs in autoimmune thyroid diseases and their role as biomarkers. Best Pract Res Clin Endocrinol Metab. 2023;37(2):101741. https://doi.org/10.1016/j.beem.2023.101741.
Fuziwara CS, Kimura ET. MicroRNAs in thyroid development, function and tumorigenesis. Mol Cell Endocrinol. 2017;456:44–50. https://doi.org/10.1016/j.mce.2016.12.017.
Massolt ET, Chaker L, Visser TJ, Gillis AJM, Dorssers LCJ, Beukhof CM, et al. Serum microRNA profiles in athyroid patients on and off levothyroxine therapy. PLoS One. 2018;13(4):e0194259. https://doi.org/10.1371/journal.pone.0194259.
Agarwal V, Bell GW, Nam JW, Bartel DP. Predicting effective microRNA target sites in mammalian mRNAs. Elife. 2015;4. https://doi.org/10.7554/eLife.05005.
Kozomara A, Birgaoanu M, Griffiths-Jones S. miRBase: from microRNA sequences to function. Nucleic Acids Res. 2019;47(D1):D155–62. https://doi.org/10.1093/nar/gky1141.
Ibarra A, Vega-Guedes B, Brito-Casillas Y, Wagner AM. Diabetes in pregnancy and MicroRNAs: promises and limitations in their clinical application. Noncoding RNA. 2018;4(4). https://doi.org/10.3390/ncrna4040032.
Dypas LB, Gutzkow KB, Olsen AK, Duale N. MiRNA profiles in blood plasma from mother-child duos in human biobanks and the implication of sample quality: Circulating miRNAs as potential early markers of child health. PLoS One. 2020;15(4):e0231040. https://doi.org/10.1371/journal.pone.0231040.
McDonald SD, Walker MC, Ohlsson A, Murphy KE, Beyene J, Perkins SL. The effect of tobacco exposure on maternal and fetal thyroid function. Eur J Obstet Gynecol Reprod Biol. 2008;140(1):38–42. https://doi.org/10.1016/j.ejogrb.2008.02.010.
We would like to thank professor and senior consultant Jan Kvetny for inspiration and help to initiate this study. Jan Kvetny was a dedicated endocrinologist that participated in the team as the project was launched, but are now retired.
Moreover, we would like to thank laboratory technician Jette Ellehauge for her skilled efforts in doing the qPCR-analyses. Also, we would like to thank Jane Clements for English proofreading.
The study was funded by Region Zealand Research Fund and Naestved, Slagelse, Ringsted (NSR) Hospital Research Fund. Both funds sponsored materials, laboratory costs and salary to the main author during the main study period. Moreover, the Dept. of Gynecology and Obstetrics at Slagelse Hospital provided financial support for leave of absence for 2 months of writing.
Ethics approval and consent to participate
The study was approved by the Danish Data Protection Agency, and the Regional Ethic Committee of Zealand (SJ-361). Every adult participant received oral and written information. Both parents had to accept inclusion of the child to the study by written consent. The study conformed to the Declaration of Helsinki. The study was registered in Clinical Trials with identification number NCT02061111.
Consent for publication
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Stryhn, J.K.G., Larsen, J., Pedersen, P.L. et al. Expressions of mitochondria-related genes in pregnant women with subclinical hypothyroidism, and expressions of miRNAs in maternal and cord blood. Thyroid Res 16, 38 (2023). https://doi.org/10.1186/s13044-023-00180-6