- Open Access
Multi-tissue gene-expression analysis in a mouse model of thyroid hormone resistance
Genome Biology volume 5, Article number: R31 (2004)
Resistance to thyroid hormone (RTH) is caused by mutations of the thyroid hormone receptor β (TRβ) gene. To understand the transcriptional program underlying TRβ mutant-induced phenotypic expression of RTH, cDNA microarrays were used to profile the expression of 11,500 genes in a mouse model of human RTH.
We analyzed transcript levels in cerebellum, heart and white adipose tissue from a knock-in mouse (TRβPV/PV mouse) that harbors a human mutation (referred to as PV) and faithfully reproduces human RTH. Because TRβPV/PV mice have elevated thyroid hormone (T3), to define T3-responsive genes in the context of normal TRβ, we also analyzed T3 effects in hyperthyroid wild-type gender-matched littermates. Microarray analysis revealed 163 genes responsive to T3 treatment and 187 genes differentially expressed between TRβPV/PV mice and wild-type littermates. Both the magnitude and gene make-up of the transcriptional response varied widely across tissues and conditions. We identified genes modulated in T3-dependent PV-independent, T3- and PV-dependent, and T3-independent PV-dependent pathways that illuminated the biological consequences of PV action in vivo. Most T3-responsive genes that were dysregulated in the heart and white adipose tissue of TRβPV/PV mice were repressed in T3-treated wild-type mice and upregulated in TRβPV/PV mice, suggesting the inappropriate activation of T3-suppressed genes in RTH.
Comprehensive multi-tissue gene-expression analysis uncovered complex multiple signaling pathways that mediate the molecular actions of TRβ mutants in vivo. In particular, the T3-independent mutant-dependent genomic response unveiled the contribution of a novel 'change-of-function' of TRβ mutants to the pathogenesis of RTH. Thus, the molecular actions of TRβ mutants are more complex than previously envisioned.
Thyroid hormone (T3) regulates growth, development and differentiation. These actions are mediated by high-affinity thyroid hormone receptors (TRs) that bind T3 and localize to the nucleus, where they regulate transcription of target genes. Four T3-binding TR isoforms, β1, β2, β3 and α1, are derived from two TR genes (β and α genes) by alternative splicing of the primary transcripts [1–3]. The expression of these TR isoforms is tissue dependent and developmentally regulated [1, 2, 4]. As transcription factors, the TRs regulate gene expression by binding thyroid-hormone response elements (TREs) in the regulatory domains of target genes which can confer specificity for TR isoforms . The transcriptional activity of the TRs is also modulated by a host of co-repressors and co-activators .
Mutations in TRβ that affect its ability to bind T3 or interact with co-repressors results in the syndrome known as resistance to thyroid hormone (RTH) [7–9]. RTH is a dominantly inherited abnormality that manifests as an impaired sensitivity of thyroid hormone-responsive tissues to circulating T3 [7, 9]. Clinical diagnosis of RTH typically recognizes elevated levels of thyroid hormone associated with nonsuppressible thyroid-stimulating hormone (TSH) [7, 9]. The clinical presentation of the disease is hypothyroidism, with symptoms such as delayed growth, cognitive dysfunction, and hypercholesterolemia, and, concurrently, signs consistent with hyperthyroidism, including tachycardia, weight loss, attention deficit-hyperactivity disorder and advanced bone age. The hypothyroid-like effects are presumably the consequence of mutant TRβ interference with, or inhibition of, normal T3 signaling pathways, whereas the signs reflective of hyperthyroidism result from the elevated T3 driving the activity of the TRα1 isoform .
To investigate the physiological consequences of a mutant β receptor in the germline, a mouse model expressing a TRβ mutant was created using homologous recombination and the Cre/loxP system . The targeted TRβ mutation (referred to as PV) was the same as that from a patient with severe RTH whose symptoms included elevated T3 and T4, nonsuppressible TSH, goiter, tachycardia and short stature . The PV mutant was shown in vitro to lack both T3-binding activity and transcriptional capacity and to strongly interfere with the transcriptional activity of normal TRs . The PV mutant was found expressed in all mouse T3-target tissues examined, including cerebrum, cerebellum, pituitary, liver, brown and white adipocytes, heart, muscle, lung, spleen and kidney . TRβPV mice exhibited goiter, delayed bone development, impaired weight gain, hearing defects and hypercholesterolemia, all of which are reminiscent of the clinical presentation of RTH in humans [11, 13, 14].
The availability of this mouse model provides an extraordinary opportunity to understand the molecular actions of TRβ mutants and to elucidate the affected cellular pathways in vivo. To this end, we applied high-density cDNA microarrays to identify PV-affected genes in the heart, cerebellum and white adipose tissue (WAT) of TRβPV/PV mice. Concurrently, we analyzed the expression profiles of wild-type littermates treated with T3 to uncover T3-responsive genes in the context of normal TRβ in a model of pharmacologic hyperthyroidism. Across the three tissues analyzed, we identified (by twofold change or more) 163 distinct genes responsive to T3 treatment and 187 genes differentially expressed between TRβPV/PV mice and wild-type littermates.
The expression patterns of these genes showed a diverse transcriptional response comprising genes with tissue-specific expression and genes with similar or contrasting expression in multiple tissues. Category analysis of expression patterns identified genes that were modulated by T3-dependent PV-independent, T3- and PV-dependent, and T3-independent PV-dependent mechanisms. Intra-tissue expression analysis provided a biological glimpse of the adverse effects of PV in a tissue-dependent manner. Finally, hierarchical clustering of multi-tissue gene-expression patterns revealed evidence of discrete biological pathways including immune response, lipogenesis and cell-cycle inhibition that are modulated in multiple tissues of hyperthyroid and T3-resistant TRβPV/PV mice. These findings provide a molecular framework for understanding the variability in tissue sensitivity of RTH and provide insight into signaling pathways of mutant TRβ in vivo.
RTH, caused by mutation of the TRβ gene, is the consequence of abnormal transcription of thyroid-hormone-responsive genes in T3 target tissues. Using a mouse cDNA microarray containing 11,500 gene probes (representing over 10,000 unique named and unnamed expressed sequence tags (ESTs)), we profiled gene-expression patterns in mouse models of hyperthyroidism and RTH to reveal T3-responsive (hyperthyroid-associated) and mutant PV-dysregulated genes. The three T3-target tissues profiled in this study were the cerebellum, heart and WAT, all of which are important T3-target tissues in which T3-responsive genes have yet to be comprehensively analyzed. The experimental comparisons in this study included: wild-type mice versus TRβPV/PV littermates (the RTH group); wild-type mice injected with T3 versus those receiving saline injection only (the iT3 group); and wild-type mice versus a second set of wild-type animals selected in the same way (the control group).
In the RTH experiments, we identified changes in gene expression resulting from the in vivo action of PV and/or elevated T3 characteristic of clinical RTH. The iT3 experiments revealed genes that are responsive to increased levels of T3 in pharmacologic hyperthyroidism. These experiments allowed us to analyze the expression profiles in RTH that are the effects due to the increased levels of circulating thyroid hormone and thus to understand which TR-subtype-dependent transcriptional mechanisms govern gene expression. The purpose of the control group was to identify genes that naturally fluctuate as the result of normal inter-animal differences . Genes identified as differentially expressed in the control group were subsequently censored in a tissue-dependent manner. For each experimental comparison, replicate microarray hybridizations were performed with reciprocal labeling (dye-swapping) to minimize technical noise (see Materials and methods).
Identification and characterization of T3-responsive genes
Differentially expressed genes (outliers) were selected based on the 2.0r criterion whereby the fold change detected on each of two reciprocally labeled arrays (r designation) was greater than twofold (up or down, but in a reciprocal manner) in at least one experimental comparison. The average expression ratios were calculated as described (see Materials and methods) and are used herein as the expression ratio measurements. To evaluate the probe-level reproducibility of the microarray data, we examined the variability of the expression ratios for all genes represented by multiple probes on the array and identified as outliers in multiple experimental comparisons (see Additional data file 1). A high degree of expression-ratio reproducibility was observed, even for changes less than twofold, lending confidence to our microarray results and outlier selection criteria. Notably, the 2.0r criterion was applied to achieve a high level of stringency in detecting gene outliers. However, this approach precludes the identification of genes that are differentially expressed at levels below 2.0r detection (for example, 1.9-fold change) thus obscuring the distinction of tissue-specific outliers and excluding potentially important distinguishing gene cassettes. In the light of this, we compiled a dataset that lists all genes passing the 2.0r criterion in at least one experimental comparison and identifies all instances (in other experimental comparisons) where these genes pass a less stringent threshold of 1.2r (see Additional data files 2 and 3). Thus, a transcript that shows a minimum change of twofold in one comparison can be subsequently screened for smaller changes, or categorized as having virtually no change in expression (that is, less than 1.2-fold), in other comparisons.
To gauge the accuracy of our methodology in correctly identifying T3-regulated genes, we searched for the 2.0r outliers across the three tissues for genes or gene products previously known to be regulated by T3. Representative genes are listed in Additional data file 4 and include 10 genes whose mRNA expression have been reported previously to be regulated by T3 and five genes that are known to be affected by T3 at the protein level. Two genes categorized as effectors are also listed. These include transthyretin (Ttr), a high-affinity serum thyroid-hormone-binding protein that transports thyroid hormone to target tissues, thus facilitating thyroid hormone distribution , and the retinoid X receptor α (Rxra), a member of the nuclear receptor superfamily that heterodimerizes with TRs and modulates the transcriptional activity of TRs .
To confirm the validity of the microarray results, four outliers (two up- and two downregulated genes) each were selected from the cerebellum and WAT for the determination of mRNA expression by reverse transcription PCR (RT-PCR) analysis. The patterns and fold changes in the expression of adenomatous polyposis coli (APC), somatostatin, carbonic anhydrase 4 and tenascin C in the cerebellum (Figure 1a) and those of lysyl oxidase, carbonic anhydrase 4, enolase 3β and creatine kinase in the WAT (Figure 1b) determined by RT-PCR were all consistent with those obtained by the arrays. Even though there was a slight variability in the magnitude of the response between the two methods, the patterns of response were in concordance.
Variability in tissue-dependent transcriptional response in mice treated with T3 and in TRβPV/PVmice
The number of genes that are differentially expressed in a tissue following T3 treatment can be viewed as an indirect measure of pharmacodynamic effect. For example, tissues that show a large number of gene outliers in response to T3 treatment would suggest higher sensitivity to T3 than tissues responding with a few outliers. This has been found to be true for chemotherapeutic effects, and for estrogen responsiveness . This principle can therefore be used as a metric for assessing tissue sensitivity to T3 treatment or the PV mutant. A breakdown of the number of 2.0r outliers identified in each tissue of each comparison is shown in Table 1. The number of outliers varied markedly between tissue types. For example, in the WAT, 58 and 172 outliers were detected in each of the comparisons (that is, with versus without T3 treatment (iT3 mice) and wild-type versus TRβPV/PV mice (RTH group), respectively), whereas in the cerebellum only 21 were found in the iT3 group and none was found in the RTH group. In the heart, 98 outliers were found in the iT3 group and 17 in the RTH group. Therefore, adult cerebellum appears to be the tissue least sensitive to T3 treatment and the action of PV, whereas the heart is a particularly responsive tissue in our model of hyperthyroidism, and white adipose is the most sensitive tissue in our model of RTH.
Intriguingly, the distribution patterns of up- and downregulated genes showed a consistent trend in two tissues that distinguished the hyperthyroid and RTH conditions (Table 1). In both the heart and WAT, the predominant response to T3 administration (that is, iT3 group) was mostly downregulation of expression (71% and 76% of the outliers, respectively). By contrast, in the TRβPV/PV mice we observed higher proportions of transcriptionally induced genes. In the heart and WAT of RTH animals, 65% and 53% of outliers, respectively, were upregulated genes. Together, these data support the intriguing view that the hyperthyroid phenotype is largely mediated through T3-induced suppression of gene expression, whereas RTH is underscored by an inappropriate increase in transcriptional activity in T3 target tissues. These data strongly suggest that target-tissue response to T3 administration and the RTH genotype are distinctly different.
To further characterize the target-tissue transcriptional response, we analyzed each tissue for its degree of transcriptional uniqueness. We found that in the T3-treated mice, 24% (5/21), 24% (24/98), and 22% (13/58) of outliers in the cerebellum, heart and WAT, respectively, were tissue specific (that is, differentially expressed by twofold or more in only one tissue, and not more than 1.2-fold in any other tissue). In the RTH mice, 0% (0/0), 53% (9/17), and 47% (81/172) of outliers in the cerebellum, heart and WAT, respectively, were detected in only one tissue type (data not shown). These observations reveal a considerable degree of transcriptional specificity that is dependent on tissue type, consistent with the notion that the transcriptional behavior of TRs depends on the molecular composition of the cell type with respect to the availability and levels of cofactors that modulate TR transcriptional activity [14, 18].
Gene classification by transcriptional response
The pathological hallmark of RTH is elevated thyroid hormone associated with nonsuppressible TSH. Even though the phenotypic consequences of elevated thyroid hormone in RTH patients are known, it is not clear what genes mediate the phenotypes. Furthermore, RTH may manifest as hyperthyroidism in one tissue, and, simultaneously, hypothyroidism in another . Dissecting this variable phenotype into the hyperthyroid and non-hyperthyroid components at the molecular level requires a comprehensive knowledge of the hyperthyroid-related genes in each tissue, and the expression patterns of these genes in the context of RTH. Therefore, we performed a tissue-by-tissue analysis cross-comparing the genomic effects of iT3 and RTH to segregate genes according to cognate response patterns.
Four categories of response patterns were discerned and labeled A, B, C and D, as shown in Figure 2. The outliers in category A (brown bars in Figure 2) are the genes that showed the same directionality of change in both T3-treated and TRβPV/PV mice (that is, up in both or down in both). That the mutant TRβ does not alter the transcriptional response of genes responsive to elevated T3 suggests their potential involvement in the hyperthyroid phenotype. As shown in Figure 2, two genes were identified in the cerebellum, five in the heart, and 39 in the WAT.
Category B (orange bars) includes the outlier genes that showed reversed patterns of change between iT3 and RTH mice (that is, up in iT3 but down in RTH, or vice versa). Presumably, these genes, which are responsive to increased T3 levels (in iT3 mice), are expressed in the opposite direction in TRβPV/PV mice as a consequence of mutant TRβ activity. One such outlier was identified in the cerebellum, 25 in the heart, and 33 in the WAT. Interestingly, 84% (21/25) and 79% (26/33) of these genes in the heart and WAT, respectively, were downregulated in response to T3 treatment and concomitantly upregulated in RTH, further supporting the view that an important molecular aspect of RTH is the abnormal expression of otherwise T3-repressed genes. A mechanistic explanation is that these category B genes are regulated by T3 mostly (or exclusively) through the TRβ receptor.
The majority of the outliers identified in this study, however, were found exclusively in either iT3 (category C, blue bars) or RTH (category D, pink bars) mice. Genes in category C showed a greater than twofold change in iT3 mice and less than 1.2-fold change in RTH mice. Therefore, these genes were responsive to increased T3 in the iT3 group but insensitive to elevated T3 in the presence of the mutant TRβ receptor. In this category, we identified 18, 72 and 27 such genes in the cerebellum, heart and WAT, respectively. Conversely, the outliers in category D were detected in RTH but not in the iT3 group. Compared with category C, a relatively smaller number of such genes were identified in the cerebellum (0 genes) and the heart (3 genes), but a significantly larger number of outliers (116 genes) were detected in the WAT. Taken together, these data show that the expression patterns of iT3 and RTH responsive genes are not only tissue-dependent, but also can differ markedly between the hyperthyroid and RTH states within the same tissue.
Hierarchical clustering identifies biological pathways associated with hyperthyroid (iT3) and RTH phenotypes
Hierarchical clustering of gene-expression patterns can provide a robust, composite view of cellular pathways operative in a biological system as evidenced by the coordinate expression of functionally related genes. To gain insight into the pathways most perturbed in our models, we analyzed the expression patterns of genes differentially expressed in one or more tissues of hyperthyroid (iT3) or RTH mice. In Figure 3, the expression patterns of these genes were hierarchically clustered and the functionality of genes within clusters was analyzed using Gene Ontology (GO) terms. We identified three gene clusters with ostensible biological implications: A, an immunity cluster: B, a lipogenesis cluster; and C, a cell cycle/growth inhibitory cluster. The immunity cluster is largely characterized by a trans-tissue downregulation of genes involved in immune response pathways or immune-cell biology (A, black bar, Figure 3). These genes are repressed in one or more tissues of T3-treated and/or TRβPV/PV mice and are all consistently downregulated in the WAT of TRβPV/PV mice. These include several HLA class II antigen genes (H2-Ab1, H2-Aa, and H2-Eb1), chemotactic factor genes (Ccl22 and Ccl19), and genes involved in lymphocyte activation (Lcp2, Ptprcap and Ms4a1) and adherence (Sell), strongly suggesting an immunomodulatory response in tissues of both hyperthyroid and TRβPV/PV mice.
In the lipogenesis cluster (B, blue bar, Figure 2), six of eight genes have clearly defined roles in fatty-acid and lipid metabolism and include Fasn, Acly, Gpd1, Elovl6, Slc25a1, and a gene named 'similar to acetyl-coenzyme A carboxylase, clone IMAGE: 5151139' which has 99% identity at the protein level with rat acetyl-coenzyme A carboxylase. These genes are consistently repressed by T3 in the heart and WAT, and expressed at reduced levels in the WAT of TRβPV/PV mice, with half of the genes being induced in the heart of TRβPV/PV mice (blue bar, Figure 3). These patterns are consistent with the negative regulation of lipogenesis in the hyperthyroid heart and WAT. In RTH mice this suggests suppression of lipogenesis in the WAT with simultaneous increase in lipogenesis in the heart.
The smaller cell-cycle/growth-inhibitory cluster is composed of several genes known to inhibit growth or cell-cycle progression or promote apoptosis (C, pink bars, Figure 3). These genes include Bnip3, Gadd45a and Cdkn1a (p21). Shown adjacent to this gene cluster is the similarly expressed gene, APC (adenomatous polyposis coli), whose expression is also associated with growth inhibition (that is, via suppression of the Wnt signaling pathway which we recently showed to be negatively regulated by T3 ). A hallmark of these genes (with the exception of Gadd45a) is increased expression in all three tissues of the hyperthyroid mouse. Notably, p21 is also upregulated in the heart and WAT of TRβPV/PV mice, suggesting its induction in these tissues in response to high circulating T3. These findings raise the possibility that induction of growth-inhibitory mechanisms in multiple T3 target tissues may contribute to the pathogenesis of RTH.
Gene-expression patterns in hyperthyroid and RTH hearts are distinct
The heart is highly sensitive to thyroid hormone and the primary mode of T3 action is a direct influence on cardiac gene expression . Increased heart rate (tachycardia) and myocardial contractility are clinical features common to both hyperthyroidism and RTH [21, 22]. As the predominant TR isoform in the heart is TRα1, and RTH patients have normal expression levels of TRα1 together with elevated thyroid hormone, it has been postulated that the cardiac effects common to both hyperthyroidism and RTH are mediated by TRα1 signaling [20, 23]. To examine this hypothesis at the molecular level, we compared the genomic effects of induced hyperthyroidism and RTH in the mouse heart. In total, we identified 105 gene outliers that showed twofold change or more in the iT3 or RTH groups (see Additional data file 5). Notably, this list does not include several genes previously implicated in T3-mediated cardiac effects such as those for the α- and β-myosin heavy chains, SERCA2 and HCN2, as they were not represented on our microarray. To focus on the set of most relevant genes in our list, we used the GO database of biological and functional gene classifications  to define genes known or expected to have a role in cardiac muscle contraction.
Hence, we classified 23 genes into one or more of the following categories (Table 2): muscle-related; calcium-ion binding (calcium-ion release and retrieval is central to muscle contraction and relaxation); ATP binding and electron transport (ATP is the immediate source of energy that powers muscle contraction); mitochondrion related (cardiac muscle has the richest supply of mitochondria and T3 directly boosts energy metabolism in mitochondria via gene transcription).
We hypothesized that the expression patterns of these genes would be similar in the iT3 and RTH hearts as they are presumably regulated by TRα1 and responsive in the context of elevated T3. Surprisingly, however, this analysis indicated that the iT3 and RTH hearts are molecularly distinct. Of the 23 genes identified, none showed the same or similar response patterns. Eighteen genes, which ranged from fourfold induction to 20-fold suppression in the iT3 heart, showed less than a 1.2-fold change in RTH (that is, category C genes, Figure 2). Moreover, the remaining five genes were expressed in opposite directions (that is, category B genes, Figure 2). These include the gene for myosin light chain, regulatory A (> 20-fold repressed in iT3 and > 6.6-fold increased in RTH) and the sarcolipin gene (> 4.8-fold repressed in iT3 and > 3.6-fold increased in RTH).
This expression pattern trend was also apparent in the remaining 82 non-muscle-contraction-related outliers (see Additional data file 5). Of these genes, 54 showed from 2- to 11-fold change in response to T3 injection and concurrently no change in RTH mice, while only three showed a greater then twofold change in RTH with no concurrent change in iT3 mice. Of the 25 genes in this group that showed change in both iT3 and RTH mice, only five showed change in the same direction (that is, candidates for TRα1-dependent transcription, category A genes, Figure 2), whereas the remaining 20 genes showed opposing expression patterns. Such marked differences between iT3 and RTH hearts suggest a greater role for TRβ in the heart than previously envisioned, or interference of TRα1-mediated transcription by TRβPV or its downstream effects.
Functional interpretation of expression profiles in the WAT of RTH mice
Thyroid hormone is an important regulator of adipose tissue development and metabolism. T3 is thought to exert its effects on fat cells through transcriptional regulation by both α1 and β receptors . Although a number of T3 target genes in adipose tissues have been described [26, 27], little is known of the transcriptional response of adipocyte T3 target genes in RTH. Accordingly, we compared WAT expression patterns in the iT3 and RTH groups. In total, we identified 215 genes that showed twofold change or more in either iT3 or RTH (see complete list in Additional data file 6). Assignment of the named genes to GO biological processes (tier 3) revealed four predominant biological classes: metabolism (n = 45), cell growth and maintenance (n = 24), cell communication (n = 24) and response to external stimuli (n = 16). The expression patterns of genes in these categories provide a molecular framework for interpreting the biological properties of WAT in RTH. The length limitation of this paper permits only description of a selected subset of genes. For example, within the cell communication class we identified a number of genes upregulated in RTH that are involved in cell-cell adhesion via desmosomal junctions (Table 3). Desmosomes are major intercellular adhesive junctions that provide strong mechanical attachments between adjacent cells and are thought to have a role in tissue morphogenesis and differentiation . The adhesive core of the desmosome is comprised of proteins of the cadherin family, namely desmogleins and desmocollins. In the WAT of RTH animals, we observed a marked upregulation of transcript levels of cadherin 1, desmoglein 2, desmocollin 2 and the CEA-related cell adhesion molecule 1. Desmosomes are anchored to intermediate filaments of keratin in the cytoplasm. Here, we also observed increased transcript levels of the keratins Krt1-13, Krt1-19 and Krt2-6a. In addition, the most highly upregulated transcript in the WAT of RTH mice was a gene called 5730453H04Rik, which shares 95% identity (over 480 amino acids) with human desmoplakin, a major protein of desmosomes involved in the anchoring of keratin intermediate filaments to desmosomes.
Within the metabolism class, we identified a number of genes negatively regulated in the WAT of iT3 and RTH mice that are directly involved in fatty-acid and lipid metabolism (Table 3). These genes include those previously discovered in the lipogenesis cluster of Figure 3 (that is, Fasn, Acly, Gpd1, Elovl6, Slc25a1 and IMAGE: 5151139). Fatty-acid biosynthesis is known to use large amounts of NADPH. As a consequence, adipose tissue is known to express high levels of the enzymes of the pentose phosphate pathway that generate NADPH. Hexose-6-phosphate dehydrogenase and phosphogluconate dehydrogenase are two of the three major enzymes in the pentose phosphate pathway, and both were negatively regulated in the WAT of RTH animals.
The genes assigned to the class response to external stimuli were, with only one exception, all downregulated in RTH animals, and in some cases, in iT3 mice (Table 3). All the genes identified here are involved in immune response or immune-cell biology and overlap with the previously identified immunity-gene cluster of Figure 3. They include a number of chemokines, HLA antigens and other markers of lymphoid lineages.
Taken together, these data suggest that the WAT of RTH mice is characterized by increased cell-cell adhesion via desmosomal structures, and they recapitulate our initial findings of a possible hyperthyroid-associated inhibition of lipogenesis and a repression of immunomodulatory signaling.
Gene-expression patterns in the cerebellum of iT3 and RTH mice
Thyroid hormones are essential for normal brain development and regulate neuronal proliferation, migration and synaptogenesis in the cerebellum [29, 30]. In our study, the cerebellum was the least responsive tissue in both iT3 and RTH mice, with only 13 genes and no genes, respectively showing twofold or greater changes. The complete list of genes with all GO classifications is shown in Additional data file 7. The most highly induced gene (> threefold) in the cerebellum of T3-treated mice was the T3/T4-binding protein, transthyretin, which is synthesized in large amounts in the choroid plexus and is essential for the transport of thyroid hormones from the blood to the brain . Increased expression of transthyretin may represent a positive feedback mechanism for enhancing the effect of thyroid hormone in the brain. mRNA for the retinoid X receptor (RXR) was also induced more than threefold in response to T3. RXR is a member of the nuclear receptor superfamily that heterodimerizes with TRs and modulates the transcriptional activity of TRs . T3-induced upregulation of RXR may reflect the relative importance of TR-RXR heterodimerization for T3 action in the cerebellum. Also induced by T3 treatment were the genes APC and dickkopf homolog 3, which are known or suspected antagonists of the Wnt signaling pathway. Both T3 and Wnt pathways are known to have important roles in neuronal growth and synaptogenesis [32–35]. We recently presented genetic and biochemical evidence that thyroid hormone can inhibit the Wnt signaling pathway . Thus, T3-mediated inhibition of Wnt signaling may be an important aspect of cerebellar development and function. The physiologic impact of the activation of these genes in hyperthyroidism warrants further investigation.
In this study, we sought to define the molecular basis of the target-tissue phenotype for hereditary TRβ mutations using a global gene expression approach. Specifically, we used expression profiling as a complex readout for the in vivo actions of T3 and PV in target tissues. The animal model was a murine line with the TRβ disrupted in the identical manner as found in the human condition (TRβPV mice). This TRβPV mouse faithfully reproduces human RTH [11, 13, 14]. To study the in vivo effects of T3 in the presence of normal TRβ, we used a mouse model of pharmacologic hyperthyroidism by treating the wild-type siblings with T3 (iT3 mice). This model aided in the delineation of genes associated with the hyperthyroid and RTH phenotype. Gene-expression patterns were analyzed in three T3 target tissues - cerebellum, heart and WAT - previously shown to have marked variation in metabolic patterns after T3 treatment.
Several interesting observations emerged from our study. First, we found that the normal physiologic effect of T3 on specific tissues is predominantly downregulation of gene expression. In the T3-treated mice, the majority of expression outliers in the heart (71%) and WAT (76%) were repressed. The same phenomenon has been observed in other mouse in vivo expression studies, albeit using less comprehensive array technologies. Feng et al.  identified 55 outliers in the livers of hypothyroid mice injected with T3; of these, 14 (25%) were upregulated and 41 (75%) were downregulated. Similarly, in a mouse model of T3-induced involution of a thyrotrope tumor, Wood et al.  identified 47 T3-responsive genes: seven (15%) were upregulated whereas 40 (85%) were downregulated. Thus, taken together, these studies indicate that T3 predominantly suppresses gene expression in responsive tissues. As would be predicted, in the RTH model, the majority of genes in the heart (63%) and, to a lesser extent in the WAT (53%), were upregulated. Of the 25 genes in the heart and 33 genes in WAT that showed opposite expression patterns in RTH (PV) and iT3 mice, 84% and 79% in heart and WAT, respectively, were repressed by T3 treatment and induced in RTH. Together, these data suggest that certain RTH (and perhaps hypothyroid) tissue phenotypes may be due predominantly to the inappropriate induction of genes normally repressed by T3. Using the number of differentially expressed genes as a semi-quantitative measure of pharmacodynamic effect, we found that the cerebellum is by far the least responsive in both iT3 and RTH, whereas the heart is the most responsive in iT3 (approximately fivefold more than cerebellum and approximately twofold more than WAT; Table 1) and the WAT is the most responsive in RTH (around 10-fold more than heart, while the cerebellum showed no changes in gene expression greater than twofold; Table 1).
Category analysis of expression patterns in iT3 and RTH mice shed light on the transcriptional consequences of the PV mutant (Figure 2). Category A consists of genes for which the T3-induced response was not affected by the expression of PV, indicating that the responses of these genes to increased T3 are likely to be exclusively mediated by TRα1. However, we have recently demonstrated that PV can form inactive heterodimers with TRα1, thereby suppressing the transcriptional activity of TRα1 . This suggests the possibility that genes in category A may, in some instances, be modulated by TR-independent pathways: a notion supported by recent studies demonstrating that T3 can act in alternate pathways independent of TRs . Categories B and C consist of T3-response genes mediated by TRs, but their transcriptional responses are altered by PV. In category B, the directionality of response in iT3 and RTH is reversed. These genes may be of particular importance in the phenotypic manifestations of RTH. Notable examples are the TSHβ and the α-glycoprotein common subunit genes that are repressed in the mouse pituitary following T3 administration . However, in the pituitary of TRβPV/PV mice, the expression of the PV mutant leads to the induction of these two genes, resulting in the highly elevated serum TSH .
Genes in category C were responsive to T3 treatment yet insensitive to hyperthyroid levels of T3 in the presence of PV. These genes are not differentially expressed between untreated euthyroid wild-type mice and TRβPV/PV littermates and thus may reflect changes specific to pharmacologic hyperthyroidism. Alternatively, these genes may be primary targets of TRα1 and their expression levels in TRβPV/PV mice mimic euthyroid levels as a result of partial transcriptional interference by PV. Some of the interesting genes in this category include cytochrome c oxidase, subunit VIII and Rous sarcoma oncogene (Table 2).
One particularly interesting finding concerns the genes in category D. These genes did not respond to T3 treatment, but showed altered expression in the presence of PV. WAT had a particularly high number of genes in this category (53% (116/215) compared to heart (3%) (3/105)). The modulation of these genes in RTH could arise from the systemic effects of T3 resistance. Alternatively, the transcriptional response of these genes could be mediated by PV, either directly or indirectly, in a T3-independent manner. For example, we have previously shown that PV can heterodimerize with the 9-cis retinoacid receptor, RXR . Heterodimerization of a mutant TRα with other nuclear receptors could simultaneously alter other hormone signaling pathways in RTH in a T3-independent manner.
The experimental design of this study afforded a unique opportunity to investigate the involvement of cellular pathways that may contribute to the pathophysiology of mutant TRβ action in the intact animal. Through pattern analysis of genes identified as differentially expressed in iT3 or TRβPV/PV mice in one or more tissues, we identified several gene clusters with clear functional associations that may have pathological implications in thyroid hormone disease (Figure 1). Then, on a tissue-by-tissue basis, we were able to focus on expression patterns and pathways most relevant to the tissue type.
The lipogenesis cluster in Figure 3 consists of six genes (Fasn, Acly, Gpd1, Elovl6, Slc25a1 and the homolog of rat Acac) that were repressed in both iT3 and RTH mice (category A genes; Table 1). These genes are known to have important roles in fatty-acid and lipid synthesis. This finding presumably reflects the known function of T3 as a regulator of lipogenesis and lipolysis in adipocytes  and, specifically, as a transcriptional modulator of Fasn, Acly and Gpd1 [40–42]. Whereas T3 is known to induce lipogenesis in the liver , downregulation of fatty-acid synthase and acetyl-CoA carboxylase genes was reported in brown adipose tissue of T3-treated hypothyroid rats . This is consistent with the physiologic findings that hyperthyroid rats display diminished fatty-acid synthesis in brown and white adipose tissues . That all six of these genes were downregulated in the WAT of both iT3 and RTH mice suggests that the T3-mediated inhibition of lipogenesis in this tissue is independent of TRβ and may be the main targets of TRα1 in this tissue. This possibility is supported by the observations of multiple TR transcriptional response elements in the promoters of both Fas and Gpd1 [43, 44] and the observation that mice heterozygous for a TRα1 point mutation (resulting in a dominant-negative TRα1 mutant with elevated T3 and T4) showed a more than fourfold increase in body fat compared to wild-type mice . The suppression of these genes in both T3-treated and TRβPV/PV mice could represent an important mechanism in the phenotype of body-fat reduction common to both hyperthyroidism and RTH.
We identified in Figure 3 a small cluster of genes with similar expression patterns that are composed of cell-cycle inhibitory and apoptosis-promoting genes. Cdkn1a (encoding p21, a major cyclin-dependent kinase inhibitor that induces G1 arrest), APC (a negative regulator of Wnt-induced cell proliferation), and Bnip3 (a potent promoter of apoptosis) are all induced by T3 treatment in the cerebellum, heart and WAT. The gene for the growth inhibitor Gadd45a is induced in the heart of iT3 mice and the WAT of RTH mice, in a fashion similar to that of Cdkn1a. Though the mitogenic activity of thyroid hormone has been documented in vitro [46, 47], little is known of the cellular mechanisms underlying the antiproliferative [48, 49] and apoptotic effects [50, 51] of T3. Recently, we and others have shown that T3 treatment of cultured cells can induce expression of the p53 tumor suppressor gene [19, 52]. Several of the genes in this T3-inducible cluster are known transcriptional targets of p53, including Gadd45a , Cdkn1a  and APC , suggesting that a coordinated transcriptional cassette supportive of growth arrest induced by T3 is likely to be mediated directly and/or indirectly by TRα1. The induction of an antiproliferative/tumor suppressor pathway is also consistent with our previous findings that the Wnt pathway is coordinately repressed by T3 in rat pituitary cells . v-erbA is a viral oncogene whose oncogenic function is the inhibition of the thyroid hormone receptor through a dominant-negative effect. These observations support our hypothesis that a major mechanism for the transforming function of v-erbA is its cellular inhibition of T3/TRα1 effects, resulting in coordinate upregulation of oncogenic pathways (that is, the Wnt signaling pathway), and repression of antiproliferation and tumor suppressor signals (for example, p53, Gadd45a, Cdkn1a).
Finally, we discovered an overall repression in the expression of immune-related genes in iT3 and RTH mice (Figure 3). To further understand the functions of these genes, we analyzed the expression of all outlier genes from the Figure 3 dataset with known roles in immune response or immune-cell biology. The expression profiles of 28 genes meeting these criteria are shown in Figure 4. Interestingly, all but six of these 28 genes showed reduced expression in the WAT of TRβPV/PV mice, and half were repressed in the hearts of iT3 mice. That these genes represent lineage-specific markers of a variety of immune cells, and that they are consistently reduced in two tissues (the heart (iT3) and WAT (RTH)), suggests that at least a fraction of these genes may reflect a general lack of lymphocytes in these tissues of hyperthyroid and RTH mice. Indeed, it has been observed that hypothyroid mouse strains have significantly reduced numbers of pro-B, pre-B, and B cells, and treatment of these mice with T4 increases the number of B-lineage cells . In addition, increases in numbers of splenic and thymic T cells and splenic NK cells have also been observed in mice treated long term with T3 . Intriguingly, mouse primary adipocytes themselves are known to be capable of expressing multiple chemokines and inflammatory mediators . As shown in Figure 4, we also observed reduced expression of chemokines CCL22, CCL6 and CXCL13 (which together are chemotactic for B cells, T cells, monocytes and macrophages) and induction of Socs5 (a suppressor of cytokine signaling) in the WAT of RTH mice. The repression of these genes could contribute to an attenuation of lymphocytic residence in this tissue.
The heart is a major target of thyroid hormone action, and T3 has a direct effect on heart rate and contractile function. Both TRα1 and TRβ are expressed in cardiac myocytes; however, TRα1 is the predominant subtype, and recent reports have suggested a primary role for TRα1 in mediating T3 effects in the heart [10, 59, 60]. It has been hypothesized that, in cardiac tissue from TRβPV/PV mice, the primary effects of the elevated T3 would be mediated through heightened T3/TRα1 signaling and that the expression response would be similar, if not identical, to that found in the myocardium of wild-type mice with induced hyperthyroidism by exogenous administration of T3. Contrary to this hypothesis, we observed a large degree of discordant cardiac gene expression when we compared the responsive genes of T3-treated animals having wild-type TRα1 and TRβ with those of TRβPV/PV mice with wild-type TRα1 (Table 2). The hearts of T3-treated animals showed expression of genes indicative of an altered workload. When heart workload increases or decreases, there is a switch in energy source from predominantly the oxidation of fatty acids to the oxidation of glucose . Three key regulators facilitate this switch: malonyl-CoA decarboxylase (MCD), uncoupling protein 3 (UCP3) and pyruvate dehydrogenase kinase 4 (PDK4). In our study, both MCD (which regulates β-oxidation of long-chain fatty acids) and UCP3 (which uncouples the oxidative phosphorylation of ADP) were upregulated at the mRNA level in the hearts of T3-treated animals, with concommitant downregulation of all the genes of the lipogenesis cluster (Figure 3). In contrast, the RTH heart did not show induction of MCD or UCP3, but rather shows increased mRNA levels for the lipogenic genes Fasn, Elovl6 and Slc25a1, suggesting a different metabolic signaling in the RTH heart. Recently, Swanson et al. reported that the TRβPV/PV mice, under euthyroid conditions, have decreased heart rate and contractile function, demonstrating that despite a lower expression of TRβ as compared with TRα1 in the mouse heart, the homozygous PV mutant can negatively interfere with TRα1 signaling in this organ . Our genomic results are in agreement with these physiologic findings and provide a glimpse at the gene-by-gene tissue response to T3 stress.
The clinical manifestations of RTH and hyperthyroidism are well documented, and although the genetic cause of RTH is known, the target-tissue transcriptional responses and cellular pathways affected are far from being understood. In our study we have used expression profiling to dissect the molecular consequences of a dominant-negative TRβ homozygous mutation in vivo, and in the process, develop a better understanding of TRα1 effects in the whole animal. Our results show that T3 primarily acts to repress gene expression, and that TRβ has a greater modulating effect in the heart than originally thought. Moreover, we identified novel physiologic candidates for more subtle T3 action such as changes in immune-gene expression, and in the induction of antiproliferative genes. Lastly, our analysis further confirms that the relative levels of TR isoforms lead to dramatic differential effects on gene expression. The results permit the identification of genetic elements that contribute to downstream isoform effects in individual tissues [59, 60].
Materials and methods
Mouse strains and treatment
The animal protocols used in the present study have been approved by the National Cancer Institute Animal Care and Use Committee. The TRβPV mice contain a cytosine insertion in exon 10 of the TRβ1 gene at nucleotide position 1,642 of the TRβ1 cDNA. This mutation leads to a frameshift of the carboxy-terminal 14 amino acids of TRβ1 [11, 12]. TRβPV mice were prepared via homologous recombination and maintained on a mixed background of 129/Sv × C57BL/6J . TRβ+/+ and TRβPV/PV male siblings at 8-10 weeks of age were used in this study. Another group of wild-type (TRβ+/+) male littermates of TRβPV/PV received daily intraperitoneal injections of T3 (5 μg/mouse) or the same volume of vehicle (phosphate-buffered saline) for 7 days. On the eighth day the animals were sacrificed by CO2 exposure and tissues (cerebellum, heart and white adipose) were harvested.
Total RNAs were extracted from liquid nitrogen-frozen tissue specimens using Trizol reagent according to the manufacturer's instructions (Gibco-BRL). To minimize differences in gene expression owing to inter-animal variability, for each tissue and animal category, 5 μg of total RNA from four animals were pooled before RNA amplification. RNA amplification and quality assessment was carried out as previously described .
Microarray design and manufacture
The mouse cDNA microarrays contained 11,520 sequences derived from various sequence-verified clone collections including the Incyte GEM1 mouse clone set (8,700; Incyte), the Research Genetics mouse sequence-verified set (2,600; Research Genetics,), and a private collection of 200 mouse ESTs (a gift from Lothar Heninghausen, NCI, Bethesda, MD). All arrays were manufactured at the NCI Microarray Facility at the Advanced Technology Center, Gaithersburg, MD. Arrays were printed using an OmniGrid microarrayer (GeneMachines) on poly-L-lysine-coated-glass according to Eisen and Brown .
Microarray hybridization and scanning
For each experimental comparison, duplicate microarray hybridizations were performed according to a reciprocal strategy termed dye-swapping. In this approach, the Cy3 and Cy5 labeling scheme for the two arrays are reversed such that a true differentially expressed gene having a ratio of 2.0 (that is, a 2-fold red (Cy5-biased) spot) on one array, for example, would have a reciprocal ratio of approximately 0.5 (that is, a 2-fold green (Cy3-biased) spot) on the other array. Notably, this form of experimental replication is more rigorous than conventional duplication (that is, true replicates using the same labeling scheme) as it takes into account a recognized form of experimental noise in microarray technology characterized by the reproducible bias of one fluorophore over the other at low signal intensity features (termed dye bias). Three micrograms amplified RNA were used to generate Cy3- or Cy5-labeled cDNA for each array hybridization. cDNA labeling reaction and array hybridization were carried out essentially as described . Microarrays were scanned on a GenePix 4000A microarray scanner (Axon Instruments) to generate 16-bit TIFF images of Cy3 and Cy5 signal intensities. The images were analyzed using GenePix Pro 3.0 microarray analysis software to measure fluorescence signals and format data for database deposition.
Outlier selection and data analysis
All array data were deposited in the NCI-CIT microarray database  where Cy3 and Cy5 signals were median normalized and expression ratios (Cy5/Cy3) were calculated. The following selection criteria were applied. All spots having a mean signal (after background subtraction) less than twice that of background in both Cy3 and Cy5 channels were systematically excluded from the dataset to eliminate measurements based on poorly detected features. The data were also filtered to exclude spots flagged as missing or corrupt in one or both arrays of a dye-swap pair. Remaining outliers were then selected by the 2.0r criterion, whereby an expression ratio must demonstrate a twofold or greater change on each array of a pair of reciprocal dye-swap hybridizations. We also demanded that the ratios must be in opposite expression directions (that is, > 2.0 on one array, and < 0.5 on the other) to exclude 'false-positive outliers' that showed a reproducible bias towards either red or green fluorescence and thus failed to reciprocate in the dye-swap array pairs. We next calculated the average expression ratio from the dye-swap pairs (after taking the reciprocal ratio from one of the two arrays) and report the average expression ratios in all analyses herein. For reporting genes by name, IMAGE clone IDs corresponding to the microarray probe sequences were used to extract UniGene Cluster IDs and names (UniGene Build 129). For genes represented by multiple probes on the array, a single probe was used (that is, the probe that showed the greatest transcriptional variation) to report the expression ratio so as to prevent gene number biasing owing to redundant probes.
The expression of selected genes was also analyzed by quantitative RT-PCR using the following pairs of primers:
APC forward primer 5'-GCTGACATCTGTGCTGTGGATGG-3';
APC reverse primer 5'-TCCTTAAAGCTGCTGCACTTCCC-3';
somatostatin forward primer 5'-TGCATCGTCCTGGCTTTGG-3';
somatostatin reverse primer 5'-GAGTTAAGGAAGAGATATGGG-3';
tenascin C forward primer 5'-TTCGTGTGTTCGCCATCTTG-3';
tenascin C reverse primer 5'-GTGTGAGGTCGATGGTGGT-3';
carbonic anhydrase 4 forward primer 5'-TCACTGCTAGGACAAAGGTG-3';
carbonic anhydrase 4 reverse primer 5'-AGAGGTTGAATGGGGTTTGG-3';
lysyl oxidase forward primer 5'-CTACATCCAGGCTTCCACG-3';
lysyl oxidase reverse primer 5'-TCTCCTCTGTGTGTTGGCAT-3';
enolase 3β, muscle forward primer 5'-TCAAGGGTTACACTCTGCCT-3';
enolase 3β, muscle reverse primer 5'-TCGTTCACCCACACAAGGAA-3';
creatine kinase, muscle forward primer 5'-TTCCTTGTGTGGGTGAACGA-3';
creatine kinase, muscle reverse primer 5'-TTTTCCAGCTTCTTCTCCATC-3';
GAPDH forward primer 5'-ACATCATCCCTGCATCCACT-3';
GAPDH reverse primer 5'-GTCCTCAGTGAGCCCAAG-3'.
Total RNA (5 μg) prepared as described above was used in the quantitative RT-PCR reaction. cDNA was prepared by SuperScript II reverse transcriptase I (Invitrogen) in the presence of poly(dT) primer. Preliminary experiments were carried out to define the conditions under which the amplification by PCR was in the linear range. The final conditions used for each were: 94°C for 30 sec, 60°C for 30 sec, 72°C for 30 sec followed by 72°C for 2 min. A total of 26 cycles were used for the amplification of APC, carbonic anhydrase 4, lysyl oxidase and somatostatin cDNA. A total of 32 cycles were used for the amplification of tenascin C, enolase 3β (muscle) and creatine kinase (muscle) cDNA. The PCR products were analyzed by 2% agarose gel electrophoresis and detected by ethidium bromide staining. The intensities of the bands were quantified by Eagle Eye (Stratagene). The expression of GAPDH was determined similarly (22 cycles) as an internal control. The band intensities of the quest genes were normalized against the control. The data are expressed as mean value ± SEM (n = 3 independent determinations). Differences between groups were examined for statistical significance using Student's t-test. p < 0.05 is considered statistically significant.
Additional data files
The following additional data files are available with the online version of this article and also at the authors' website : an Excel table (Additional data file 1) of all redundant microarray probes and their average expression ratio measurements; an Excel file (Additional data file 2) containing a comprehensive filtered gene list of average expression ratios; a figure (Additional data file 3) showing a hierarchical clustergram of 2.0r outlier genes with full annotation; a Word table (Additional data file 4) listing the responsive genes known to be modulated by T3; an Excel table (Additional data file 5) listing genes differentially expressed in the heart of iT3 and/or RTH mice, while another table (Additional data file 6) lists genes differentially expressed in the white adipose tissue (WAT) of iT3 and/or RTH mice. Also provided is an Excel file (Additional data file 7) listing genes that are differentially expressed in the cerebellum of iT3 and/or RTH mice; and, finally, a zipped text file (Additional data file 8) containing the raw data.
Cheng SY: Multiple mechanisms for regulation of the transcriptional activity of thyroid hormone receptors. Rev Endocr Metab Disord. 2000, 1: 9-18. 10.1023/A:1010052101214.
Yen PM: Physiological and molecular basis of thyroid hormone action. Physiol Rev. 2001, 81: 1097-1142.
Harvey CB, Williams GR: Mechanism of thyroid hormone action. Thyroid. 2002, 12: 441-446. 10.1089/105072502760143791.
O'Shea PJ, Williams GR: Insight into the physiological actions of thyroid hormone receptors from genetically modified mice. J Endocrinol. 2002, 175: 553-570.
Forrest D, Vennstrom B: Functions of thyroid hormone receptors in mice. Thyroid. 2000, 10: 41-52.
McKenna NJ, Lanz RB, O'Malley BW: Nuclear receptor coregulators: cellular and molecular biology. Endocrinol Rev. 1999, 20: 321-344. 10.1210/er.20.3.321.
Refetoff S, Weiss RE, Usala SJ: The syndromes of resistance to thyroid hormone. Endocr Rev. 1993, 14: 348-399. 10.1210/er.14.3.348.
Yoh SM, Chatterjee VK, Privalsky ML: Thyroid hormone resistance syndrome manifests as an aberrant interaction between mutant T3 receptors and transcriptional corepressors. Mol Endocrinol. 1997, 11: 470-480. 10.1210/me.11.4.470.
Weiss RE, Refetoff S: Resistance to thyroid hormone. Rev Endocr Metab Disord. 2000, 1: 97-108. 10.1023/A:1010072605757.
Zhang XY, Kaneshige M, Kamiya Y, Kaneshige K, McPhie P, Cheng SY: Differential expression of thyroid hormone receptor isoforms dictates the dominant negative activity of mutant beta receptor. Mol Endocr. 2002, 16: 2077-2092. 10.1210/me.2002-0080.
Kaneshige M, Kaneshige K, Zhu X, Dace A, Garrett L, Carter TA, Kazlauskaite R, Pankratz DG, Wynshaw-Boris A, Refetoff S, et al: Mice with a targeted mutation in the thyroid hormone beta receptor gene exhibit impaired growth and resistance to thyroid hormone. Proc Natl Acad Sci USA. 2000, 97: 13209-13214. 10.1073/pnas.230285997.
Parrilla R, Mixson AJ, McPherson JA, McClaskey JH, Weintraub BD: Characterization of seven novel mutations of the c-erbA beta gene in unrelated kindreds with generalized thyroid hormone resistance. Evidence for two "hot spot" regions of the ligand binding domain. J Clin Invest. 1991, 88: 2123-2130.
Griffith A, Szymko Y, Kaneshige M, Quinonez R, Kaneshige K, Heintz K, Masteroianni M, Kelly M, Cheng SY: Knock-in mouse model for resistance to thyroid hormone (RTH): an RTH mutation in the thyroid hormone receptor β gene disrupts cochlear morphogenesis. J Assoc Res Otolaryngol. 2002, 3: 279-288. 10.1007/s101620010092.
Kamiya Y, Zhang XY, Ying H, Kato Y, Willingham MC, Xu J, O'Malley BW, Cheng SY: Modulation by steroid receptor coactivator-1 of target-tissue responsiveness in resistance to thyroid hormone. Endocrinology. 2003, 144: 4144-4153. 10.1210/en.2003-0239.
Pritchard CC, Hsu L, Delrow J, Nelson PS: Project normal: defining normal variance in mouse gene expression. Proc Natl Acad Sci USA. 2001, 98: 13266-13271. 10.1073/pnas.221465998.
Robbins J: Thyroid hormone transport proteins and the physiology of hormone binding. In The Thyroid. Edited by: Braverman LE, Utiger RD. 2000, Philadelphia, PA: Lippincott William and Wilkins, 105-120. 8
Sotiriou C, Powles TJ, Dowsett M, Jazaeri AA, Feldman AL, Assersohn L, Gadisetti C, Libutti SK, Liu ET: Gene expression profiles derived from fine needle aspiration correlate with response to systemic chemotherapy in breast cancer. Breast Cancer Res. 2002, 4: R3-10.1186/bcr433.
Sadow PM, Chassande O, Gauthier K, Samarut J, Xu J, O'Malley BW, Weiss RE: Specificity of thyroid hormone receptor subtype and steroid receptor coactivator-1 on thyroid hormone action. Am J Physiol Endocrinol Metab. 2003, 284: E36-E46.
Miller L, Park K, Guo Q, Alkharouf N, Malek R, Lee N, Liu E, Cheng SY: Silencing of Wnt signaling and activation of multiple metabolic pathways in response to thyroid hormone-stimulated cell proliferation. Mol Cell Biol. 2001, 21: 6626-6639. 10.1128/MCB.21.19.6626-6639.2001.
Dillmann WH, Gloss BR: The role of thyroid hormone receptors in the heart. Methods Mol Biol. 2002, 202: 55-70. 10.1385/1-59259-174-4:55.
Fadel BM, Ellahham S, Ringel MD, Lindsay J, Wartofsky L, Burman KD: Hyperthyroid heart disease. Clin Cardiol. 2000, 23: 402-408.
Refetoff S, Weiss RE, Usala SJ: The syndromes of resistance to thyroid hormone. Endocr Rev. 1993, 14: 348-399. 10.1210/er.14.3.348.
Kahaly GJ, Matthews CH, Mohr-Kahaly S, Richards CA, Chatterjee VK: Cardiac involvement in thyroid hormone resistance. J Clin Endocrinol Metab. 2002, 87: 204-212. 10.1210/jc.87.1.204.
Gene Ontology Consortium. [http://www.geneontology.org]
Reyne Y, Nougues J, Cambon B, Viguerie-Bascands N, Casteilla L: Expression of c-erbA alpha, c-erbA beta and Rev-erbA alpha mRNA during the conversion of brown adipose tissue into white adipose tissue. Mol Cell Endocrinol. 1996, 116: 59-65. 10.1016/0303-7207(95)03696-2.
Viguerie N, Millet L, Avizou S, Vidal H, Larrouy D, Langin D: Regulation of human adipocyte gene expression by thyroid hormone. J Clin Endocrinol Metab. 2002, 87: 630-634. 10.1210/jc.87.2.630.
Pucci E, Chiovato L, Pinchera A: Thyroid and lipid metabolism. Int J Obes Relat Metab Disord. 2000, S109-S112. 24 Suppl 2
Huber O: Structure and function of desmosomal proteins and their role in development and disease. Cell Mol Life Sci. 2003, 60: 1872-1890. 10.1007/s00018-003-3050-7.
Koibuchi N, Chin WW: Thyroid hormone action and brain development. Trends Endocrinol Metab. 2000, 11: 123-128. 10.1016/S1043-2760(00)00238-1.
Poguet AL, Legrand C, Feng X, Yen PM, Meltzer P, Samarut J, Flamant F: Microarray analysis of knockout mice identifies cyclin D2 as a possible mediator for the action of thyroid hormone during the postnatal development of the cerebellum. Dev Biol. 2003, 254: 188-199. 10.1016/S0012-1606(02)00039-8.
Chanoine JP, Braverman LE: The role of transthyretin in the transport of thyroid hormone to cerebrospinal fluid and brain. Acta Med Austriaca. 1992, 19 Suppl 1: 25-28.
Vincent J, Legrand C, Rabie A, Legrand J: Effects of thyroid hormone on synaptogenesis in the molecular layer of the developing rat cerebellum. J Physiol (Paris). 1982, 78: 729-738.
Madeira MD, Paula-Barbosa MM: Reorganization of mossy fiber synapses in male and female hypothyroid rats: a stereological study. J Comp Neurol. 1993, 337: 334-352.
Hall AC, Lucas FR, Salinas PC: Axonal remodeling and synaptic differentiation in the cerebellum is regulated by WNT-7a signaling. Cell. 2000, 100: 525-535. 10.1016/S0092-8674(00)80689-3.
Salinas PC: Wnt factors in axonal remodelling and synaptogenesis. Biochem Soc Symp. 1999, 65: 101-109.
Feng X, Jiang Y, Meltzer P, Yen PM: Thyroid hormone regulation of hepatic genes in vivo detected by complementary DNA microarray. Mol Endocrinol. 2000, 14: 947-955. 10.1210/me.14.7.947.
Wood WM, Sarapura VD, Dowding JM, Woodmansee WW, Haakinson DJ, Gordon DF, Ridgway EC: Early gene expression changes preceding thyroid hormone-induced involution of a thyrotrope tumor. Endocrinology. 2002, 143: 347-359. 10.1210/en.143.2.347.
Davis PJ, Tillmann HC, Davis FB, Wehling M: Comparison of the mechanisms of nongenomic actions of thyroid hormone and steroid hormones. J Endocrinol Invest. 2002, 25: 377-388.
Shupnik MA, Chin WW, Habener JF, Ridgway EC: Transcriptional regulation of the thyrotropin subunit genes by thyroid hormone. J Biol Chem. 1985, 260: 2900-2903.
Goodridge AG: Regulation of the gene for fatty acid synthase. Fed Proc. 1986, 45: 2399-2405.
Gharbi-Chihi J, Facchinetti T, Berge-Lefranc JL, Bonne J, Torresani J: Triiodothyronine control of ATP-citrate lyase and malic enzyme during differentiation of a murine preadipocyte cell line. Horm Metab Res. 1991, 23: 423-427.
Sugisaki T, Noguchi T, Beamer WG, Kozak LP: Genetic hypothyroid mice: normal cerebellar morphology but altered glycerol-3-phosphate dehydrogenase in Bergmann glia. J Neurosci. 1991, 11: 2614-2621.
Xiong S, Chirala SS, Hsu MH, Wakil SJ: Identification of thyroid hormone response elements in the human fatty acid synthase promoter. Proc Natl Acad Sci USA. 1998, 95: 12260-12265. 10.1073/pnas.95.21.12260.
Weitzel JM, Kutz S, Radtke C, Grott S, Seitz HJ: Hormonal regulation of multiple promoters of the rat mitochondrial glycerol-3-phosphate dehydrogenase gene: identification of a complex hormone-response element in the ubiquitous promoter B. Eur J Biochem. 2001, 268: 4095-4103. 10.1046/j.1432-1327.2001.02332.x.
Liu YY, Schultz JJ, Brent GA: A thyroid hormone receptor alpha gene mutation (P398H) is associated with visceral adiposity and impaired catecholamine-stimulated lipolysis in mice. J Biol Chem. 2003, 278: 38913-38920. 10.1074/jbc.M306120200.
Barrera-Hernandez G, Park KS, Dace A, Zhan Q, Cheng SY: Thyroid hormone-induced cell proliferation in GC cells is mediated by changes in G1 cyclin/cyclin-dependent kinase levels and activity. Endocrinology. 1999, 140: 5267-5274. 10.1210/en.140.11.5267.
Pibiri M, Ledda-Columbano GM, Cossu C, Simbula G, Menegazzi M, Shinozuka H, Columbano A: Cyclin D1 is an early target in hepatocyte proliferation induced by thyroid hormone (T3). FASEB J. 2001, 15: 1006-1013. 10.1096/fj.00-0416com.
Gonzalez-Sancho JM, Figueroa A, Lopez-Barahona M, Lopez E, Beug H, Munoz A: Inhibition of proliferation and expression of T1 and cyclin D1 genes by thyroid hormone in mammary epithelial cells. Mol Carcinog. 2002, 34: 25-34. 10.1002/mc.10046. (Erratum in: Mol Carcinog 2002, 34:164)
Perez-Juste G, Aranda A: The cyclin-dependent kinase inhibitor p27(Kip1) is involved in thyroid hormone-mediated neuronal differentiation. J Biol Chem. 1999, 274: 5026-5031. 10.1074/jbc.274.8.5026.
Nishikawa A, Kaiho M, Yoshizato K: Cell death in the anuran tadpole tail: thyroid hormone induces keratinization and tail-specific growth inhibition of epidermal cells. Dev Biol. 1989, 131: 337-344.
James RA, Sarapura VD, Bruns C, Raulf F, Dowding JM, Gordon DF, Wood WM, Ridgway EC: Thyroid hormone-induced expression of specific somatostatin receptor subtypes correlates with involution of the TtT-97 murine thyrotrope tumor. Endocrinology. 1997, 138: 719-724. 10.1210/en.138.2.719.
Dinda S, Sanchez A, Moudgil V: Estrogen-like effects of thyroid hormone on the regulation of tumor suppressor proteins, p53 and retinoblastoma, in breast cancer cells. Oncogene. 2002, 21: 761-768. 10.1038/sj.onc.1205136.
Kastan MB, Zhan Q, el-Deiry WS, Carrier F, Jacks T, Walsh WV, Plunkett BS, Vogelstein B, Fornace AJ: A mammalian cell cycle checkpoint pathway utilizing p53 and GADD45 is defective in ataxia-telangiectasia. Cell. 1992, 71: 587-597.
el-Deiry WS, Tokino T, Velculescu VE, Levy DB, Parsons R, Trent JM, Lin D, Mercer WE, Kinzler KW, Vogelstein B: WAF1, a potential mediator of p53 tumor suppression. Cell. 1993, 75: 817-825.
Jaiswal AS, Narayan S: p53-dependent transcriptional regulation of the APC promoter in colon cancer cells treated with DNA alkylating agents. J Biol Chem. 2001, 276: 18193-18199. 10.1074/jbc.M101298200.
Foster MP, Montecino-Rodriguez E, Dorshkind K: Proliferation of bone marrow pro-B cells is dependent on stimulation by the pituitary/thyroid axis. J Immunol. 1999, 163: 5883-5890.
Watanabe K, Iwatani Y, Hidaka Y, Watanabe M, Amino N: Long-term effects of thyroid hormone on lymphocyte subsets in spleens and thymuses of mice. Endocr J. 1995, 42: 661-668.
Ruan H, Zarnowski MJ, Cushman SW, Lodish HF: Standard isolation of primary adipose cells from mouse epididymal fat pads induces inflammatory mediators and down-regulates adipocyte genes. J Biol Chem. 2003, 278: 47585-47593. 10.1074/jbc.M305257200.
Gloss B, Trost S, Bluhm W, Swanson E, Clark R, Winkfein R, Janzen K, Giles W, Chassande O, Samarut J, Dillmann W: Cardiac ion channel expression and contractile function in mice with deletion of thyroid hormone receptor alpha or beta. Endocrinology. 2001, 142: 544-550. 10.1210/en.142.2.544.
Gloss B, Sayen MR, Trost SU, Bluhm WF, Meyer M, Swanson EA, Usala SJ, Dillmann WH: Altered cardiac phenotype in transgenic mice carrying the delta337 threonine thyroid hormone receptor beta mutant derived from the S family. Endocrinology. 1999, 140: 897-902. 10.1210/en.140.2.897.
Taegtmeyer H, Cohen DM: Overestimating glycolysis in rat heart. Am J Physiol Endocrinol Metab. 2002, 283: E1102-1103. (Author reply E1103-1104)
Swanson EA, Gloss B, Belke DD, Kaneshige M, Cheng SY, Dillmann WH: Cardiac expression and function of thyroid hormone receptor beta and its PV mutant. Endocrinology. 2003, 144: 4820-4825. 10.1210/en.2003-0522.
Eisen MB, Brown PO: DNA arrays for analysis of gene expression. Methods Enzymol. 1999, 303: 179-205.
Wang E, Miller LD, Ohnmacht GA, Liu ET, Marincola FM: High-fidelity mRNA amplification for gene profiling. Nat Biotechnol. 2000, 18: 457-459. 10.1038/74546.
National Cancer Institute: Center for Cancer Research. [http://nciarray.nci.nih.gov]
Genome Institute of Singapore: supplementary information. [http://www.gis.a-star.edu.sg/homepage/toolssup.jsp]
Electronic supplementary material
Additional data file 1: An Excel table of all redundant microarray probes and their average expression ratio measurements (XLS 31 KB)
Additional data file 2: An Excel file containing a comprehensive filtered gene list of average expression ratios (XLS 181 KB)
Additional data file 3: A figure showing a hierarchical clustergram of 2.0r outlier genes with full annotation (PDF 53 KB)
Additional data file 4: A Word table listing the responsive genes known to be modulated by T3 (DOC 34 KB)
Additional data file 5: An Excel table listing genes differentially expressed in the heart of iT3 and/or RTH mice (XLS 72 KB)
Additional data file 6: An Excel table lists genes differentially expressed in the white adipose tissue (WAT) of iT3 and/or RTH mice (XLS 130 KB)
Additional data file 7: An Excel file listing genes that are differentially expressed in the cerebellum of iT3 and/or RTH mice (XLS 28 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
About this article
Cite this article
Miller, L.D., McPhie, P., Suzuki, H. et al. Multi-tissue gene-expression analysis in a mouse model of thyroid hormone resistance. Genome Biol 5, R31 (2004). https://doi.org/10.1186/gb-2004-5-5-r31
- Thyroid Hormone
- White Adipose Tissue
- Adenomatous Polyposis Coli
- Additional Data File
- Average Expression Ratio