- Open Access
Sex impacts cardiac function and the proteome response to thyroid hormone in aged mice
Proteome Science volume 18, Article number: 11 (2020)
Sex and age have substantial influence on thyroid function. Sex influences the risk and clinical expression of thyroid disorders (TDs), with age a proposed trigger for the development of TDs. Cardiac function is affected by thyroid hormone levels with gender differences. Accordingly, we investigated the proteomic changes involved in sex based cardiac responses to thyroid dysfunction in elderly mice.
Aged (18–20 months) male and female C57BL/6 mice were fed diets to create euthyroid, hypothyroid, or hyperthyroid states. Serial echocardiographs were performed to assess heart function. Proteomic changes in cardiac protein profiles were assessed by 2-D DIGE and LC-MS/MS, and a subset confirmed by immunoblotting.
Serial echocardiographs showed ventricular function remained unchanged regardless of treatment. Heart rate and size increased (hyperthyroid) or decreased (hypothyroid) independent of sex. Pairwise comparison between the six groups identified 55 proteins (≥ 1.5-fold difference and p < 0.1). Compared to same-sex controls 26/55 protein changes were in the female hypothyroid heart, whereas 15/55 protein changes were identified in the male hypothyroid, and male and female hyperthyroid heart. The proteins mapped to oxidative phosphorylation, tissue remodeling and inflammatory response pathways.
We identified both predicted and novel proteins with gender specific differential expression in response to thyroid hormone status, providing a catalogue of proteins associated with thyroid dysfunction. Pursuit of these proteins and their involvement in cardiac function will expand our understanding of mechanisms involved in sex-based cardiac response to thyroid dysfunction.
Thyroid hormone (TH) modulates myocardial function and cellular architecture through multiple mechanisms. Acting as a ligand, thyroid hormone can bind cardiomyocyte nuclear receptors and regulate transcription of genes responsible for myocardial contractility and systolic function [1, 2]. Thyroid hormone can also act extranuclear, activating cardiac ion channels and mitochondrial biogenesis [3, 4]. Disruptions in thyroid hormone homeostasis, even those defined as subclinical , can have substantial impact on cardiac morphology and function. Clinically, two important demographic factors independently modulate impact of thyroid dysfunction: age and sex . However, as age and sex modify numerous factors, which impose cardiovascular risk, their interaction with different forms of thyroid hormone disruption are difficult to determine in a human population. The interaction between sex and age requires clarification, as thyroid hormone abnormalities occur at higher rates in elderly and females compared respectively to young and males [7, 8].
The aging population exhibits incidences for diabetes and atherosclerosis which exceed rates in younger cohorts . These factors complicate interpretation of data from aging human populations with regards to thyroid hormone disturbances and its direct impact on the heart. In controlled mouse models without diabetes or atherosclerosis, aging alone promotes myocardial hypertrophy and cardiac dysfunction , emphasizing the importance of further examination of the impact of thyroid hormone disturbance on the aging heart.
In relation to sex and cardiovascular diseases, several studies highlight the importance of sex differences in the incidence and clinical manifestation observed between males and females [11,12,13]. Sex hormones and TH modulate each other through multiple mechanisms affecting gene transcription and signal transducing pathways [14,15,16] necessitating the need to better understand the role sex plays in modulating the impact of thyroid hormone on the aging heart. We address this issue using a proteomic approach to investigate the role sex and age play in the TH altered heart in the aged mouse.
This investigation conforms to the Guide for the Care and Use of Laboratory Animals  and all animal protocols were reviewed and approved by Institutional Animal Care and Use Committee at Seattle Children’s Research Institute.
The male and female C57/BL6 mice, aged 18–20 months old, used in this study were obtained from Charles Rivers Laboratories (Wilmington, MA). The mice were divided into six groups: male euthyroid (MCON), male hypothyroid (MPTU), male hyperthyroid (MTH), female euthyroid (FCON), female hypothyroid (FPTU), female hyperthyroid (FTH). The mice were fed one of the following diets: 1) control iodine diet (TD.97350, CON), 2) hypothyroid − 0.15% propylthiouracil (PTU) diet (TD.95125), or 3) hyperthyroid - 0.1% thyroid powder (TH) diet (TD.150769) purchased from Envigo (Madison, WI). The thyroid powder (T6384) was purchased from Sigma-Aldrich (St. Louis, MO), and dosage was based on work by Thomas et al. . All mice were fed for 4 weeks prior to sacrifice. Water was available ab libitum.
T4 serum ELISA assays
Serum samples were collected prior to commencement of diet regime by submandibular bleeds and clearing by centrifugation. Final serum samples were collected at time of sacrifices and heart collection. The T4 levels in serum were determined using the mouse T4/Thyroxine ELISA kit (LS-F25787) from LifeSpan BioSciences, Inc., (Seattle, WA). (n = 5 animals/sex/treatment).
As previously performed in our laboratory  serial echocardiograms were conducted prior to thyroid hormone intervention and after four weeks. Briefly, mice were initially sedated with 3% isoflurane in O2 at a flow of 1 LPM and placed in a supine position at which time the isoflurane is reduced to 1.5% administered via a small nose cone. ECG leads were placed for simultaneous ECG monitoring during image acquisition. Echocardiographic images were performed with a Vevo 2100 machine using a MS250 or MS550 transducer (VisualSonics, Inc., Toronto, Canada). M-Mode measurements at the midpapillary level of the left ventricle (LV) were performed at end-diastole (LVEDD) and end-systole (LVESD) to determine LV function via the fractional shortening [(LVEDD-LVESD)/LVEDD * 100] in a parasternal short axis mode for at least three heart beats. The investigators performing and interpreting the echocardiograms were blinded to the treatments. (n = 6–9 animals/sex/treatment).
Following the final echocardiogram, mice were euthanized, the hearts extracted, rinsed in PBS buffer and frozen in liquid nitrogen. The left tibial bone was dissected out and measured using a Vernier caliper.
2D-gel and LC-MS analysis
Crushed heart tissue was sent to Applied Biomics (Hayward, CA) for 2D-DIGE and LC-MS/MS. In short, samples were labeled with CyDye DIGE fluors (size and charge matched) and subjected to IEF and SDS-PAGE electrophoresis. After electrophoresis, the gel was scanned using Typhoon TRIO (GE Healthcare, Chicago, IL). Each scan reveals one of the CyDye signals. The ImageQuant (version 6) (GE Healthcare, Chicago, IL) software is used to generate the image presentation data including the single and overlay images. A comparative analysis of all spots was performed using DeCyder software (version 6.5) (GE Healthcare, Chicago, IL). The protein expression ratio cut-off of 1.5-fold difference and p ≤ 0.1 was used in comparing same-sex control and treatment groups or males to females of same treatment. Protein spots of interest were picked from the 2D gel with the Ettan Spot Picker (Amersham Biosciences, Little Chalfont, UK). The 2D gel spots underwent in-gel trypsin digestion to generate peptides. The peptides were subjected to MALDI-TOF MS and TOF/TOF tandem MS/MS performed on an AB SCIEX TOF/TOF™ 5800 System (AB SCIEX, Framingham, MA). The resulting peptide mass and associated fragmentation spectra were submitted to GPS Explorer workstation equipped with MASCOT search engine (Matrix science) to search the database of National Center for Biotechnology Information non-redundant (NCBInr) or Swiss-Prot database. Candidates with either protein score C.I.% or Ion C.I.% greater than 95 were considered significant. An n = 3 animals/sex/treatment was used for these experiments. The mass spectrometry data have been deposited to the ProteomeXchange Consortium (http://proteomecentral.proteomeexchange.org) via the PRIDE  partner repository with the dataset identifier PXD022557.
Protein extraction and Western blotting
Heart tissue was homogenized in RIPA buffer (20 mM Tris-HCl (pH 7.5), 150 mM NaCl, 1MM Na2EDTA, 0.1% SDS, 1% sodium deoxycholate, 1X protease/phosphatase inhibitor (Thermo Fisher Scientific)). Thirty micrograms of total protein extract were resolved using SDS-PAGE and electroblotted onto PVDF membranes. Standard immunoblotting protocols were followed. Horseradish peroxidase secondary antibodies were used and visualized with enhanced chemiluminescence using the ChemiDoc-It imager™ (AnalytiKjena, Germany). Membranes were stripped by washing them 2 × 15 min with 100 mM DTT, 2% (wt/vol) SDS, and 62.5 mM Tris·HCl, pH 6.7, at 70 °C, followed by three 10-min washes with TBS. Efficiency of protein transfer to PVDF membrane was determined using Pierce Reversible Protein Stain Kit for PVDF Membranes (no. 24585) (Thermo Fisher Scientific, Rockford, IL). Densitometry analyses were performed using Image J and Quantity One (Bio-Rad, Hercules, CA). All immunoblots were performed in duplicate. The primary antibodies used in this study were as follows: α 1B-glycoprotein (AF8537) purchased from R&D Systems (Minneapolis, MN) and ACAA2 (GTX115417) purchased from GeneTex (Irvine, CA). An n = 4 animals/sex/treatment was used for these experiments.
Bioinformatics and motif analysis
The gene-annotation enrichment analysis software, DAVID Bioinformatics Resources 6.8 [21, 22], was utilized to acquire molecular function (MF) and KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway  information to protein sequences obtained. A search tool for the retrieval of interacting genes/proteins (STRING, http://string-db.org/) was used to construct protein-protein interactions between the proteins. The searchable database of regulatory elements and their inferred target genes, GeneHancer embedded in GeneCards (www.genecards.org) [24, 25], was used to ascertain potential thyroid hormone receptor binding sites.
All reported values are expressed as mean ± standard error of the mean. Significance between baseline and post-treatment values was assessed using two-way ANOVA with Tukey post hoc analysis where appropriate. One-way ANOVA was performed for values comparing treatment to same-sex control group. For immunoblots comparing control and experimental groups of a single sex we performed one-way ANOVA with Tukey post hoc analysis; and two-way ANOVA with Tukey post hoc analysis on blots containing both sexes . Hypothyroid and hyperthyroid treatment groups were not compared to each other.
Criterion for significance was p < 0.05 for all comparisons. We preformed statistical analysis using GraphPad Prism version 7.03 (GraphPad Software).
Serum thyroid hormone levels
Efficacy of the mouse diets to induce the hypo- or hyperthyroid state was confirmed by mouse T4/thyroxine ELISA. Total T4 (TT4) serum levels were assayed prior to start and at completion of the 4-week diet regime (Fig. 1). No significant sex differences were found in the T4 serum levels of the control mice, FCON or MCON, prior to or post treatment. TH levels were altered by the diet regime. The TH deprivation diet resulted in an average decrease in TT4 levels of 17.8 ± 3.4 ng/ml and 24.5 ± 2 ng/ml in females (FPTU) and males (MPTU), respectively. Whereas, the TH diet resulted in an average increase in TT4 levels of 59 ± 2.7 ng/ml and 53.4 ± 6.1 ng/ml TT4 in females (FTH) and males (MTH), respectively.
Morphometric and echocardiographic data
We performed serial echocardiography to assess changes in left ventricular size and function compared to baseline values (Table 1) in response to treatment. LVWPD, left ventricular posterior wall diameter, an important determinant of diastolic stiffness and pressure, significantly increased for the male by 0.42 ± 0.09 mm (p = 0.02) and for the female by 0.29 ± 0.08 mm (p = 0.042) in mice fed the hyperthyroid diet. There was no significant change observed for the other functional echocardiographic parameters between pre- and post-treatment.
The male hypothyroid mice experienced a significant decrease of 83 ± 26 beats per minute (bpm) (p < 0.006) final heart rate compared to baseline; the effect did not reach significance for the females. However, both the male and female hyperthyroid mice had significantly increased heart rates. The male mice experienced a 102.8 ± 11 bpm (p < 0.0002), and the female mice a 92 ± 25 bpm (p < 0.002) increase in final heart rate compared to baseline.
Serial body weight measurements provided an assessment of systemic effect of thyroid hormone state. The male hypothyroid mice experienced a significant weight loss of 3.3 ± 0.5 g (p < 0.0001); whereas no change in the females was observed. In contrast, the hyperthyroid females experienced significant weight gain of 2.6 ± 0.5 g (p < 0.02), whereas the males remained constant.
An examination of the excised heart weight found that the male and female hypo- and hyperthyroid mice hearts had significant weight differences compared to their same sex control. The male and female hypothyroid heart weighed 33 ± 11 mg (p = 0.0025) and 31 ± 4 mg (p = 0.0021), respectively, less compared to same-sex controls. In contrast, the male and female hyperthyroid heart weighed 77 ± 8 mg (p = 0.0001) and 77 ± 10 mg (p = 0.0001), respectively, more compared to same-sex controls. The normalization of heart weight to tibia length reflected the same changes observed by heart weight measurements alone.
2D-DIGE and LC MS/MS was performed to analyze TH affected proteomic alterations in aged mouse hearts. Three biological replicates for each treatment in both sexes were used in the 2D-DIGE, allowing for the following pairwise comparisons: FCON to FPTU, FCON to FTH, MCON to MPTU, MCON to MTH, FTH to MTH and FPTU to MPTU. Image analysis of the 2D-DIGE revealed 66 protein spots showing a change in average abundance of ≥1.5 fold (p ≤ 0.1) between at least one of the six pairwise comparisons. Representative images of the 2D-DIGE gels is shown in Fig. 2. LC MS/MS analysis of the 66 spots identified 55 proteins (Table 2). Protein fold change ratios are depicted by a color gradient from black - representing the largest downregulation observed - to light gray – representing the largest upregulation. Hatch boxes represent no significant change in expression ratio. The largest number of proteins affected by TH status was observed between the FCON to FPTU comparison at 26/55 proteins. The FCON to FTH, MCON to MPTU, MCON to MTH comparisons each identified 15/ 55 proteins affected by TH state. The FTH to MTH comparison identified 8/55 proteins different between the sexes; and the FPTU to MPTU comparison identified 3/55. Interestingly, few of the same proteins were altered between the sexes for the same treatment.
Ontology and protein-protein interaction (PPI) network
To gain insights into molecular functions of the processes altered we performed a GO enrichment analysis using the DAVID software (Table 3). A cutoff p value of 0.05 was used in creating the table. The molecular functions identified for the 55 proteins by the DAVID software involved oxidoreductase and structural roles. A KEGG pathway analysis also identified predominantly cellular energetics pathways (Table 3). These results support our data as TH is known have a predominant involvement in energy metabolism and cell architecture.
To explore the functional connectivity between the proteins the STRING software was used to construct a protein-protein interaction network based on an edge confidence ≥0.4 (i.e. medium to high confidence) (Fig. 3). All 55 proteins (nodes) were analyzed and 176 pairwise associations identified. Not all proteins possessed identifiable associations within this dataset with 5 proteins unlinked. The line thickness indicates the strength of data supporting the interaction between nodes. Proteins involved in energy metabolic pathways made up the largest cluster (unenclosed). Smaller clusters involved proteins associated with structural remodeling (black square) or inflammatory response roles (black circle). Of interest, the protein changes in the female hypothyroid heart represented primarily energy-metabolism proteins; whereas, the male hypothyroid proteome expressed more structural-related and inflammatory-response protein changes. In the hyperthyroid hearts, male and female, the predominant changes involved the energy-metabolism proteins.
The GeneCard database was utilized to assess if the identified proteins altered by TH status possessed potential TH receptor binding sites. The database search identified 29/55 (53%) possessing potential TH receptor binding sites (Fig. 3).
We validated the 2D-DIGE MS/MS findings with immunoblotting for two proteins. The candidate proteins were chosen for separate reasons. A1BG (A1BG_MOUSE), alpha 1-B glycoprotein exhibited the greatest fold difference between sexes in the hyperthyroid state; and ACAA2 (THIM_MOUSE), mitochondrial 3-ketoacyl-CoA thiolase was chosen because it plays an important role in fatty acid oxidation in the heart.
First, we examined the changes between treatments in each sex separately (Fig. 4). The top panels represent 2D-DIGE images highlighting the A1BG or ACAA2 spots. The encircled spots in Fig. 4a for FCON and FTH reveal multiple isoforms of A1BG. These spots are absent in FPTU. The A1BG spots are also absent in the male A1BG 2D-DIGE panels (Fig. 4b). The immunoblots confirmed the protein expression changes due to TH status observed by 2D-DIGE MS/MS. The FPTU heart exhibited a 2-fold A1BG expression decrease (p = 0.003) compared to FCON. The undetectable levels observed in the MCON mice hearts did not change with diet. Figure 4c and d show the changes in the ACAA2 protein. Like the 2D-DIGE, the ACAA2 immunoblot shows expression decreased 1.6-fold in the MTH hearts compared to MCON (p = 0.0065), with no observable shift in the female hearts.
For comparisons between males and females, western blots with all 6 groups (male and female: eu, hypo- and hyperthyroid) were generated (Fig. 5). There was a significant difference in expression of the A1BG protein between MCON and FCON hearts (p < 0.0001), due the lack of detectable A1BG in the MCON hearts. The hyperthyroid diet did not alter this disparity between male and females; however, the FPTU hearts A1BG protein levels decreased to levels observed in the MCON and MPTU hearts. ACAA2 expression was similar between MCON and FCON. The hypothyroid diet did not alter this observation, and although the MTH ACAA2 expression trended lower than the in FTH hearts it did not reach significance (p < 0.09).
We provide the first detailed proteomic analysis of thyroid impact on the heart in male and female aged mice. Previous studies evaluating thyroid effects in the hearts of elderly mice specifically excluded females or did not provide analyses according to sex [26,27,28,29]. This is an important gap in knowledge as thyroid hormone interacts with sex hormones at multiple levels ranging from transcriptional to post-translational. We show TH state elicits differential sex-based proteomic changes, especially in the female hypothyroid heart.
Thyroid hormone is known to modulate cardiac phenotype. While hyperthyroidism generally induces myocardial hypertrophy, hypothyroidism can cause a dilated cardiomyopathy . These conditions are highly prevalent in aging populations [10, 30]. The impact of therapy can be highly variable, particularly on the heart . Thus, the principal objective of this study was to assess the proteomic changes due to TH status in the heart for both sexes. Echocardiography allows us to make longitudinal analyses of thyroid effect on cardiac function and morphology. Both sexes appeared exquisitely sensitive to the hyperthyroid state, which induced significant elevations in heart rate. The chronotropic response occurred concomitant with increasing left posterior wall thickness, but with no change in functional parameters such as fractional shortening and ejection fraction. Although, we could not directly measure stroke volume, the increase in heart rate along with stable cardiac chamber size and ejection fraction suggests that TH induced a high cardiac output state. In contrast, the hypothyroid state caused a significant decrease in heart rate in males; however, there was limited chronotropic effect in females. With similar logic and no change in echocardiographic derived dimensions, males presumably had a low cardiac output, while females exhibited no overall change due hypothyroidism. Thus, the data suggest that the hemodynamics of the female aged hearts are more resistant to the effects of low TH levels.
Altered TT4 serum levels, including echocardiographic and morphometric data confirmed that we successfully altered thyroid state in the heart with our perturbations, establishing the basis for our proteomic analyses. We then further defined the proteome for each sex in our aging hearts according to thyroid state. Surprisingly, prior to our analysis, detailed proteomic analyses have not been performed in hearts subjected to thyroid treatment.
A prior investigation  reported that sex modulates TH action in both 12 and 20 months old C5/BL/7 mice. Those investigators studied effects of thyroid state on multiple systemic parameters such as body temperature and weight, as well as muscle strength and locomotor coordination. However, they did not perform end-organ specific evaluation other than heart rate. Additionally, molecular analyses were not included. Those authors also alluded to their finding of sex-based differences in TT4 serum levels in response to treatment as a potential reason for phenotypic differences. However, we found no such differences in our aged mice. The TT4 responses to treatments were virtually identical in the two sexes. This discrepancy could be related to differences in protocol; mainly as we administered daily thyronine/thyroxine in chow, and they provided T4 by intraperitoneal injection every 48 h. Regardless, sex related effects, noted in our study, would not be due to differences in serum TT4 levels.
2D-DIGE in conjunction with mass spectrometry identified 55 cardiac proteins which responded to alterations in thyroid status. Searches of pathway analyses databases, STRING and DAVID, indicated that they are involved in processes relating to mitochondrial function, structural and inflammatory responses, each of which is integral to the function and remodeling of the heart. Additionally, a search of the GeneHancer database identified 29/55 (53%) of the proteins as potential TH receptor targets. This agrees with Ayers et al.  report that ~ 50% of TH induced genes display THRβ binding sites.
Prior studies in our laboratory have shown that both aging and chronic hypothyroidism separately impairs myocardial fatty acid oxidation in male rats [33, 34]. Accordingly, we chose the ACAA2 protein to validate our 2D-DIGE MS/MS data given this enzyme regulates the critical final step of mitochondrial fatty acid oxidation. We found that hyperthyroidism decreased ACAA2 in aged hearts only in males. Similarly, a prior study showed decreased ACAA2 protein expression in male rat hyperthyroid liver, supporting our findings .
The female hypothyroid hearts appear to undergo a more dynamic response compared to their male counterparts, with 26/55 affected proteins compared 15/55 observed for the male hearts. Many of the proteins altered are members of the mitochondrial respiratory chain, such as QCR1, a subunit of the ubiquinol-cytochrome c reductase complex, and the NDUB or NDUS proteins which participate in the mitochondrial membrane respiratory chain NADH dehydrogenase (Complex 1). This together with our physiological data suggest the female aged heart is more accommodating to the hypothyroid state. Therefore, the finding of more metabolic proteins undergoing altered expression in female hypothyroid heart may yield insight into potential mechanistic difference between sexes.
Inflammatory response proteins were also altered by thyroid status, many of which are glycoproteins. For example, A1BG, a novel thyroid hormone target, is a plasma glycoprotein of unknown function. Glycoproteins play key roles in inflammatory diseases, are known to be altered in cardiomyopathies and thyroid disorders, and have the potential to be exploited for diagnosis and/or treatment targets [36, 37]. Glycoproteins identified in this study, such as Serpina3K (SPA3K) and A1BG, have been shown to be differentially expressed in myocardial pathology [38, 39]. In this study A1BG protein was absent in the MCON hearts and highly expressed in the FCON. The FPTU A1BG expression plunged to the non-detectable levels of MCON and MPTU. Thyroid dysfunction is primarily diagnosed by the levels of thyroid stimulating hormone (TSH), a regulator of TH production. However, bioactivity of TSH is related to glycosylation modifications and age ; current clinical assays do not distinguish between TSH glycosylated forms. The expression of varied modified TSH protein in the aged population may account for discrepancies in diagnosis and prevalence observed in studies of aging patients, emphasizing the need for identifying other diagnostic targets. Xiao et al. reported on zinc-α2- glycoprotein increase and its correlation with TH levels in hyperthyroid patients; however, that study did not analyse male and females separately . We have identified several glycoproteins altered in the aged mouse by TH state and differentially expressed between the sexes. Hence, our findings show the importance of examining both males and females for the identification and use of potential biomarkers towards better discernment of subclinical thyroid states or indications for treatment.
TH regulates cardiovascular remodelling, affecting size, shape and function in the heart [42, 43]. Our study identified several cardiac remodelling proteins further validating our data. Thyroid disorders affect the myocardium similarly to that observed by the aging process, such as changes to cardiac contractility, cardiac output and systemic vascular resistance . Additionally, it is well recognized that structural and functional differences exist between male and female, such as anatomical size  and electrical activity . We identified myosin heavy chain α (MHC6), for instance, a major protein in cardiac muscle and contractile function that is known to be regulated by TH via the TH hormone nuclear receptor sites in its promoter region . Other structural proteins identified included regulatory myosin light chain 2 and 3 which interact directly with MHC6. Structure-related proteins identified, such as alpha B-crystallin interact with desmin/ actin cytoskeletal complexes, and are predicted to play a protective role during stress conditions . The structural proteins identified in addition to being altered by TH state also showed preferential sex expression. Further investigations of these proteins and their cellular interactions will aid in better understanding any sexual dimorphisms.
An interesting disparity was observed in our study in that the echocardiographic measurements showed modest phenotypic differences between males and females compared to the key differences observed for the proteomic changes. One explanation is that the proteins we observed altered by TH status are predominantly involved in myocardial substrate utilization as opposed to structural or hypertrophic changes. Supporting this rationale, previous work in our lab showed TH treated aged male mouse heart experienced significant alterations in substrate utilization with fatty acids, but exhibited limited or unchanged cardiac functions compared to baseline or aged untreated [31, 33, 34]. Thus, we would need to do specific studies on cardiac substrate utilization for fatty acids to define these sex-differences.
Limitations of the study
First, the proteins involved in thyroid hormone regulation identified by the 2D-DIGE and LC-MS/MS method used in this study are by no means meant to be an exhaustive list. Second, as a screening study, it is not within the scope of the present research to determine if the changes in protein expression levels are due to direct or indirect modulation of the proteins by thyroid hormone. Third, small sample size limits the statistical power of our data . Fourth, it is important to note a major difference between aging female mice and postmenopausal aged women; aging female mice do not have the equivalent of menopause. Aging female mice become acyclic by 11–16 months and develop estrogen deficiency but not as profoundly as aged women .
In this study we demonstrate that thyroid hormone differentially regulates the cardiac proteomic profile in an age and sex influenced manner. Although thyroid dysfunction is more prevalent in women previous animal studies have excluded females. In our study, the aged female heart had a greater number of proteomic changes in response to thyroid hormone state. The identification of numerous sex-based differentially expressed proteins altered by thyroid hormone status shows the importance of examining both sexes. The integration of sex-specific data will yield a better understanding of thyroid hormone regulation in the aged heart.
Availability of data and materials
All data used to support the findings of this study are available upon request to the corresponding author. The dataset generated in this study has been deposited in the ProteomeXchange Cosortium (http://proteomecentral.proteomeexchange.org) via the PRIDE partner repository with the dataset identifier PXD022557.
Two dimensional difference gel electrophoresis
liquid chromatography with tandem mass spectrometry
Alpha 1-B glycoprotein
Mitochondrial 3-ketoacyl-CoA thiolase
Fazio S, Palmieri EA, Lombardi G, Biondi B. Effects of thyroid hormone on the cardiovascular system. Recent Prog Horm Res. 2004;59:31–50.
Cheng SY, Leonard JL, Davis PJ. Molecular aspects of thyroid hormone actions. Endocr Rev. 2010;31(2):139–70.
Marin-Garcia J. Thyroid hormone and myocardial mitochondrial biogenesis. Vasc Pharmacol. 2010;52(3–4):120–30.
Purtell K, Roepke TK, Abbott GW. Cardiac arrhythmia and thyroid dysfunction: a novel genetic link. Int J Biochem Cell Biol. 2010;42(11):1767–70.
Klein I, Ojamaa K. Thyroid hormone and the cardiovascular system. N Engl J Med. 2001;344(7):501–9.
Morganti S, Ceda GP, Saccani M, Milli B, Ugolotti D, Prampolini R, Maggio M, Valenti G, Ceresini G. Thyroid disease in the elderly: sex-related differences in clinical expression. J Endocrinol Invest. 2005;28(11 Suppl Proceedings):101–4.
Surks MI, Hollowell JG. Age-specific distribution of serum thyrotropin and antithyroid antibodies in the US population: implications for the prevalence of subclinical hypothyroidism. J Clin Endocrinol Metab. 2007;92(12):4575–82.
Suzuki S, Nishio S, Takeda T, Komatsu M. Gender-specific regulation of response to thyroid hormone in aging. Thyroid Res. 2012;5(1):1.
Stout RW. Diabetes, Atherosclerosis, and Aging. Diabetes Care. 1990;13(Supplement 2):20–3.
Dai DF, Chen T, Johnson SC, Szeto H, Rabinovitch PS. Cardiac aging: from molecular mechanisms to significance in human health and disease. Antioxid Redox Signal. 2012;16(12):1492–526.
Maas AH, Appelman YE. Gender differences in coronary heart disease. Neth Heart J. 2010;18(12):598–602.
Subramanya V, Zhao D, Ouyang P, Ying W, Vaidya D, Ndumele CE, Heckbert SR, Budoff MJ, Post WS, Michos ED. Association of endogenous sex hormone levels with coronary artery calcium progression among post-menopausal women in the multi-ethnic study of atherosclerosis (MESA). J Cardiovasc Comput Tomogr. 2019;13(1):41–7.
Mehta LS, Beckie TM, DeVon HA, Grines CL, Krumholz HM, Johnson MN, Lindley KJ, Vaccarino V, Wang TY, Watson KE, et al. Acute myocardial infarction in women: a scientific statement from the American Heart Association. Circulation. 2016;133(9):916–47.
Hernandez A. Thyroid hormone role and economy in the developing testis. Vitam Horm. 2018;106:473–500.
Mazer NA. Interaction of estrogen therapy and thyroid hormone replacement in postmenopausal women. Thyroid. 2004;14(Suppl 1):S27–34.
Santin AP, Furlanetto TW. Role of estrogen in thyroid function and growth regulation. J Thyroid Res. 2011;2011:875125.
National Research Council. Guide for the Care and Use of Laboratory Animals. Washington, DC: National Academies Press; 2011.
Thomas TA, Kuzman JA, Anderson BE, Andersen SM, Schlenker EH, Holder MS, Gerdes AM. Thyroid hormones induce unique and potentially beneficial changes in cardiac myocyte shape in hypertensive rats near heart failure. Am J Physiol Heart Circ Physiol. 2005;288(5):H2118–22.
Zhu WZ, El-Nachef D, Yang X, Ledee D, Olson AK. O-GlcNAc Transferase promotes compensated cardiac function and protein kinase a O-GlcNAcylation during early and established pathological hypertrophy from pressure overload. J Am Heart Assoc. 2019;8(11):e011260.
Perez-Riverol Y, Csordas A, Bai J, Bernal-Llinares M, Hewapathirana S, Kundu DJ, Inuganti A, Griss J, Mayer G, Eisenacher M, et al. The PRIDE database and related tools and resources in 2019: improving support for quantification data. Nucleic Acids Res. 2019;47(D1):D442–50.
da Huang W, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4(1):44–57.
da Huang W, Sherman BT, Lempicki RA. Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 2009;37(1):1–13.
Kanehisa M, Sato Y, Furumichi M, Morishima K, Tanabe M. New approach for understanding genome variations in KEGG. Nucleic Acids Res. 2019;47(D1):D590–5.
Stelzer G, Rosen N, Plaschkes I, Zimmerman S, Twik M, Fishilevich S, Stein TI, Nudel R, Lieder I, Mazor Y, et al. The GeneCards Suite: From Gene Data Mining to Disease Genome Sequence Analyses. Curr Protoc Bioinformatics. 2016;54:1.30.1-1. 30.33.
Fishilevich S, Nudel R, Rappaport N, Hadar R, Plaschkes I, Iny Stein T, Rosen N, Kohn A, Twik M, Safran M, et al. GeneHancer: genome-wide integration of enhancers and target genes in GeneCards. Database (Oxford). 2017;2017.
do Imperio GE, Ramos IP, Santiago LA, Pereira GF, dos Santos Almeida NA, Fuziwara CS, Pazos-Moura CC, Kimura ET, Olivares EL, Ortiga-Carvalho TM. The Impact of a Non-Functional Thyroid Receptor Beta upon Triiodotironine-Induced Cardiac Hypertrophy in Mice. Cell Physiol Biochem. 2015;37(2):477–90.
Weiss RE, Korcarz C, Chassande O, Cua K, Sadow PM, Koo E, Samarut J, Lang R. Thyroid hormone and cardiac function in mice deficient in thyroid hormone receptor-alpha or -beta: an echocardiograph study. Am J Physiol Endocrinol Metab. 2002;283(3):E428–35.
Wadosky KM, Berthiaume JM, Tang W, Zungu M, Portman MA, Gerdes AM, Willis MS. MuRF1 mono-ubiquitinates TRalpha to inhibit T3-induced cardiac hypertrophy in vivo. J Mol Endocrinol. 2016;56(3):273–90.
Chen YF, Redetzke RA, Said S, Beyer AJ, Gerdes AM. Changes in left ventricular function and remodeling after myocardial infarction in hypothyroid rats. Am J Physiol Heart Circ Physiol. 2010;298(1):H259–62.
Barbesino G. Thyroid function changes in the elderly and their relationship to cardiovascular health: a mini-review. Gerontology. 2019;65(1):1–8.
Rakov H, Engels K, Hones GS, Brix K, Kohrle J, Moeller LC, Zwanziger D, Fuhrer D. Sex-specific phenotypes of hyperthyroidism and hypothyroidism in aged mice. Biol Sex Differ. 2017;8(1):38.
Ayers S, Switnicki MP, Angajala A, Lammel J, Arumanayagam AS, Webb P. Genome-wide binding patterns of thyroid hormone receptor beta. PLoS One. 2014;9(2):e81186.
Hyyti OM, Ledee D, Ning XH, Ge M, Portman MA. Aging impairs myocardial fatty acid and ketone oxidation and modifies cardiac functional and metabolic responses to insulin in mice. Am J Physiol Heart Circ Physiol. 2010;299(3):H868–75.
Ledee D, Portman MA, Kajimoto M, Isern N, Olson AK. Thyroid hormone reverses aging-induced myocardial fatty acid oxidation defects and improves the response to acutely increased afterload. PLoS One. 2013;8(6):e65532.
Severino V, Locker J, Ledda-Columbano GM, Columbano A, Parente A, Chambery A. Proteomic characterization of early changes induced by triiodothyronine in rat liver. J Proteome Res. 2011;10(7):3212–24.
Debus JD, Milting H, Brodehl A, Kassner A, Anselmetti D, Gummert J, Gaertner-Rommel A. In vitro analysis of arrhythmogenic cardiomyopathy associated desmoglein-2 (DSG2) mutations reveals diverse glycosylation patterns. J Mol Cell Cardiol. 2019;129:303–13.
Ufret-Vincenty CA, Baro DJ, Lederer WJ, Rockman HA, Quinones LE, Santana LF. Role of sodium channel deglycosylation in the genesis of cardiac arrhythmias in heart failure. J Biol Chem. 2001;276(30):28197–203.
Yang S, Chen L, Sun S, Shah P, Yang W, Zhang B, Zhang Z, Chan DW, Kass DA, van Eyk JE, et al. Glycoproteins identified from heart failure and treatment models. Proteomics. 2015;15(2–3):567–79.
McDonough CW, Gong Y, Padmanabhan S, Burkley B, Langaee TY, Melander O, Pepine CJ, Dominiczak AF, Cooper-Dehoff RM, Johnson JA. Pharmacogenomic association of nonsynonymous SNPs in SIGLEC12, A1BG, and the selectin region and cardiovascular outcomes. Hypertension. 2013;62(1):48–54.
Zabczynska M, Kozlowska K, Pochec E. Glycosylation in the Thyroid Gland: Vital Aspects of Glycoprotein Function in Thyrocyte Physiology and Thyroid Disorders. Int J Mol Sci. 2018;19(9):2792.
Xiao XH, Qi XY, Li JY, Wang YB, Wang YD, Liao ZZ, Yang J, Ran L, Wen GB, Liu JH. Serum zinc-alpha2-glycoprotein levels are elevated and correlated with thyroid hormone in newly diagnosed hyperthyroidism. BMC Endocr Disord. 2019;19(1):12.
Rajagopalan V, Gerdes AM. Role of thyroid hormones in ventricular remodeling. Curr Heart Fail Rep. 2015;12(2):141–9.
Cokkinos DV, Chryssanthopoulos S. Thyroid hormones and cardiac remodeling. Heart Fail Rev. 2016;21(4):365–72.
Wallace K, Hofmann MT. Thyroid dysfunction: how to manage overt and subclinical disease in older patients. Geriatrics. 1998;53(4):32–8 41.
Rabbia F, Grosso T, Cat Genova G, Conterno A, De Vito B, Mulatero P, Chiandussi L, Veglio F. Assessing resting heart rate in adolescents: determinants and correlates. J Hum Hypertens. 2002;16(5):327–32.
Surawicz B, Parikh SR. Differences between ventricular repolarization in men and women: description, mechanism and implications. Ann Noninvasive Electrocardiol. 2003;8(4):333–40.
Flink IL, Morkin E. Interaction of thyroid hormone receptors with strong and weak cis-acting elements in the human alpha-myosin heavy chain gene promoter. J Biol Chem. 1990;265(19):11233–7.
Bennardini F, Wrzosek A, Chiesi M. Alpha B-crystallin in cardiac tissue. Association with actin and desmin filaments. Circ Res. 1992;71(2):288–94.
Altman DG, Bland JM. Absence of evidence is not evidence of absence. Aust Vet J. 1996;74(4):311.
Nelson JF, Felicio LS, Randall PK, Sims C, Finch CE. A longitudinal study of estrous cyclicity in aging C57BL/6J mice: I. cycle frequency, length and vaginal cytology. Biol Reprod. 1982;27(2):327–39.
The authors thank Margaret Sedensky for her thoughtful comments toward improving our manuscript.
Research reported in this publication was supported by the National Heart, Lung, and Blood Institute of the National Institutes of Health under award number NIH K01HL120886 to Ledee and NIH R01HL122546 to Olson.
Ethics approval and consent to participate
Consent for publication
All authors declare that there are no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Zhu, W.Z., Olson, A., Portman, M. et al. Sex impacts cardiac function and the proteome response to thyroid hormone in aged mice. Proteome Sci 18, 11 (2020). https://doi.org/10.1186/s12953-020-00167-3
- Thyroid hormone