- Research article
- Open Access
Thyroid hormone receptor orthologues from invertebrate species with emphasis on Schistosoma mansoni
BMC Evolutionary Biology volume 7, Article number: 150 (2007)
Thyroid hormone receptors (TRs) function as molecular switches in response to thyroid hormone to regulate gene transcription. TRs were previously believed to be present only in chordates.
We isolated two TR genes from the Schistosoma mansoni and identified TR orthologues from other invertebrates: the platyhelminths, S. japonium and Schmidtea mediterranea, the mollusc, Lottia gigantean and the arthropod Daphnia pulex. Phylogenetic analysis of the DNA binding domain and/or ligand binding domain shows that invertebrate and vertebrate TRs cluster together, TRs from the vertebrates and from the jawless vertebrate (lamprey) clustered within separate subgroups, Platyhelminth TRs cluster outside of the vertebrate TR subgroups and that the schistosome TRs and S. mediterranea TRs clustered within separate subgroups.
Alignment of the C-terminus of the A/B domain revealed a conserved TR-specific motif, termed TR 'N-terminus signature sequence', with a consensus sequence of (G/P)YIPSY(M/L)XXXGPE(D/E)X.
Heterodimer formation between S. mansoni TRs and SmRXR1 suggests that the invertebrate TR protein gained the ability to form a heterodimer with RXR. ESMA analysis showed that SmTRα could bind to a conserved DNA core motif as a monomer or homodimer.
Vertebrate TR genes originated from a common ancestor of the Bilateria. TR genes underwent duplication independently in the Protostomia and Deuterostomia. The duplication of TRs in deuterostomes occurred after the split of jawless and jawed vertebrates. In protostomes, TR genes underwent duplication in Platyhelminths, occurring independently in trematode and turbellarian lineages. Using S. mansoni TRs as an example, invertebrate TRs exhibited the ability to form a dimer with RXR prior to the emergence of the vertebrate TRs and were able to bind to vertebrate TR core DNA elements as a monomer or homodimer.
Thyroid hormones (TH) play important roles in growth, development and metabolism in vertebrates. TH is synthesized in the thyroid gland under the control of thyroid-stimulating hormone (TSH) secreted by the pituitary. TSH secretion is controlled by thyrotropin-releasing hormone (TRH) which is secreted from the hypothalamus. THs are lipophilic molecules able to passively cross the membrane and bind to its receptor, the thyroid hormone receptor (TR). TRs belong to a superfamily of transcription factors called nuclear receptor (NR) superfamily, based on protein sequence similarities, structural motifs and functionality . TRs function as a molecular switch in response to the thyroid hormones T3 or T4 to activate or repress gene transcription depending on the promoter context and thyroid hormone binding status . The typical nuclear receptor contains an N-terminal A/B domain, a conserved C domain (DNA binding domain, DBD), a D domain (hinge region) and a moderately conserved E domain (ligand binding domain, LBD). The most conserved DBD contains two zinc finger motifs (CI and CII). Like all NRs, TRs regulate transcription through its binding to the promoter region of a target gene by the DBD and they activate or repress mRNA synthesis through co-regulators bound to the LBD . The specific target DNA sequence to which NRs bind is called a hormone response element (HRE). The typical HRE is a direct, inverted or everted repeat or palindrome of the DNA sequence AGGTCA. TRs can bind to the HRE as a monomer, a homodimer or as a heterodimer with RXR, another member of nuclear receptor superfamily which contributes to the specificity of the TR. TH binds to the LBD of TR which results in a conformational change in the C-terminus of the receptor. Corepressors then dissociated from the TR allowing coactivators to bind to the C-terminus of the TR in a hormone-dependent manner. TR and the coactivator complex activate the expression of the target gene [3, 4].
TR was previously believed to be an innovation of chordates as the genomes of insects (Drosophila and mosquito) and nematodes (Caenorhabditis elegans and C. briggus) do not contain TR genes [5–8]. Recently, we identified two thyroid receptor homologues in the flatworm Schistosoma mansoni , one of which was found in the S. mansoni EST database [6, 10]. The presence of TR homologues in S. mansoni demonstrated that the TR orthologue genes are present outside of chordates. However, it is still unclear whether these prostostome TRs possess the same functional domains as in vertebrate TRs. Another question is whether the TR orthologue is present in other invertebrates or just in the platyhelminth lineage? Answers to these questions will help to understand the origin of TR genes and evolution of the function of vertebrate thyroid hormone network. To begin to address these questions, we isolated cDNAs of S. mansoni TRs and demonstrated that TRs in platyhelminths are highly conserved not only in sequence similarity, but also in gene organization, protein-protein interaction and in DNA-binding ability. Furthermore, we mined the available genome data and demonstrated that TR orthologues are present in different invertebrate animals but not in Porifera or Cnidaria. Phylogenetic analysis showed that the TR orthologue likely originated from a common ancestor of the Bilateria.
TR orthologue genes in invertebrate animals
By an extensive search of available databases, predicted genes encoding TR orthologues were found in different invertebrate animals using the conserved DBD as a query. They include two genes in each of the platyhelminth species evaluated, the turbellarian Schmidtea mediterranea and the trematodes, Schistosoma mansoni and S. japonium, one each from the mollusc Lottia gigantean (owl limpet) and the crustacean Daphnia pulex (water flea) (Fig. 1 and see additional file 1). No TR orthologues were identified in any species of the Radiata.
Sequence alignment analysis
Alignment of DBD amino acid sequences show that the invertebrate TR orthologues are highly conserved (Fig. 1). All TRs (from chordate and non-chordate species) possess an identical P-box sequence, which determines the specificity of target DNA binding. The methionine at position 16 of the N-terminus of the DBD (highlighted in red in Fig. 1) was non-chordate TR-specific. Instead at this position an isoleucine is conserved in all chordate TRs examined to date. Two amino acids (His and Pro) at positions 33–34 of N-terminal of DBD were chordate TR-specific (highlighted blue in Fig. 1). The second zinc finger is less conserved in non-chordate TR orthologues compared to the first zinc finger, but still shows a high degree of conservation (Fig. 1). The T-box and A-box which are located in the C-terminal extension (CTE) of the DBD are also highly conserved in non-chordate TR orthologues. They show 75–83% and 57–100% identity to the chordate T-box and A-box, respectively. These regions in the water flea TR orthologue are highly divergent, probably due to a recombination event that happened in this region (see below). An aspartic acid in the vertebrate TRs at the 8th position of N-terminus of the T-box is highly conserved for jawed vertebrate TRs, a glutamic acid at this position is conserved in other TRs. Although the N-terminal A/B domain of nuclear receptors is divergent, there is a highly conserved motif found at the 3' end of the TR A/B domain, we termed the 'N-terminal signature sequence' (NTSS). The TR NTSS is a consensus sequence of (G/P)YIPSY(M/L)XXXGPE(D/E) (Fig. 1). In non-chordate TR orthologues, the sequence is PYIPSYMXXXGPEEP; while in chordates, three amino acids are deleted generating a sequence of GYIPYL(D/E)KE(P/Q/L). The deleted three amino acids in chordate NTSS were deduced from the non-chordate sea urchin TR (SpTR). Since the deuterostome SpTR shares the same NTSS with that of protostome TRs, this suggests that the non-chordate NTSS is the primitive state.
Functional domains of non-chordate TRs
To determine the functional domains of a non-chordate TR orthologue, cDNAs encoding the entire open reading frame (ORF) of the two S. mansoni TRs (SmTRα and SmTRβ) were isolated, sequenced and the sequences deposited in GenBank (SmTRα [GenBank:AY395038, AY395059–AY395063, AH013464) and SmTRβ [GenBank:AY395039]).
SmTRα and SmTRβ encode proteins of 1115 amino acids and 847 amino acids, respectively. Each protein exhibits a modular structure characteristic of the nuclear receptor superfamily with an A/B domain, a conserved DBD, a hinge region and a LBD (Fig. 2 and 3). Alignment of the sequences demonstrated that the LBD is less conserved. Helices 1–2 in the LBD are highly divergent, similar to that of other reported schistosome NRs [11–19]. Although the LBD is less conserved (Fig. 2), the consensus motif I (from helix 3 to helix 6) and the consensus motif II (from the middle of helix 8 to the middle of helix 10) are conserved  (Fig. 2). A putative AF2 activation domain core (AF2-AD) (helix 12) is present in both SmTRs (Fig. 2). It exhibits a high degree of conservation (represented by YLHELF and YFHELF in SmTRα and SmTRβ, respectively) in comparison to the common consensus AF2-AD core structure of ΦΦXEΦΦ, where Φ denotes a hydrophobic residue [12, 21, 22].
The two consensus LBD motifs (motif I and motif II) were examined in other non-chordate TR orthologues when their DBD-containing contig contained these regions (Fig. 4). Two amino acids in motif I (a glycine at position 21 and an aspartic acid at position 29) and the three amino acids in motif II (alanine at position 3, a leucine at position 9 and an isoleucine or leucine at postion 45) were found to be invertebrate TR-specific (Fig. 4).
Phylogenetic analysis of DBD sequences shows that all invertebrate TRs cluster with the vertebrate TR subgroup (Fig. 5), suggesting that they originated from a common ancestor gene. The two TR orthologues in the Platyhelminths, the turbellarian S. mediterranea and the trematodes (S. mansoni and S. japonicum) clustered outside of the vertebrate TR subgroups suggesting that they underwent duplication after the split of Prostostomia and Deuterostomia. Furthermore, each of the two duplicated trematode TR homologues cluster together within a subgroup, while the two turbellarian (S. mediterranea) TR homologues clustered within another subgroup. This result suggests that the TR duplication in the Platyhelminths occurred after the split of the trematodes and the turbellarians.
A high resolution phylogenetic tree was constructed using the DBD and LBD sequences (Fig. 6). Since the only full length LBD sequence was identified from S. mansoni TRs, the S. mansoni TRs were used to represent the protostome TRs. The results show that all jawed vertebrate TRs clustered within a group including TRα and TRβ subgroups. The two jawless vertebrate lamprey TRs clustered outside of the jawed vertebrate TRs. This suggests that the vertebrate TR gene underwent duplication after the split of jawless and jawed vertebrates. The urochordate Ciona TR and the two duplicated trematode SmTRs cluster outside of the vertebrate TR groups. These results confirmed those obtained from the DBD analysis which suggested that the TR gene duplication occurred independently in vertebrates and invertebrates.
The gene organization of the TR homologue genes in S. mansoni (SmTRα and SmTRβ) was identified by alignment of the cDNA sequence with the genomic sequence (Fig. 3). The genomic sequences were obtained by screening S. mansoni bacterial artificial chromosome (BAC) libraries (SmBAC1 or CHOR103) and sequencing the BAC DNA , and by Blast searching WTSI S. mansoni WGS database . Three clones (SmBAC1 3A6, 4G3 and 9D10) for SmTRα and 3 clones (CHOR103 3I15, 7K22 and 21F18) for SmTRβ were identified by screening the BAC libraries. A gDNA contig (Sm00.scaff00056.0030.1) for SmTRα and a gDNA contig (Sm00.scaff00631.0050) for SmTRβ were obtained from WTSI S. mansoni WGS database. The Contig Sm00.scaff00056.0030.1 contains a partial coding region for SmTRα, BAC DNAs (SmBAC1 3A6 and 4G3) were sequenced to generate a contig which contained full coverage of SmTRα cDNA sequence.
The ORF of SmTRα is encoded by seven exons spanning over 32 kb and the ORF of SmTRβ is encoded by six exons spanning over 21 kb (Fig. 3). All splice donor and acceptor sites for the two genes fit the GT-AG role. The A/B domain, hinge and E domains (LBD) are each encoded by 2–3 exons, respectively. The DBD of both proteins is encoded by one exon (Fig. 3). The 5' UTR and 3'UTR in both SmTRα and SmTRβ were not determined due to a short 3'UTR in SmTRα which seems not to be full length and due to the genomic contig not covering the 5'- and 3'UTR for SmTRβ.
It is known that there are four conserved splice junctions in nuclear receptor genes [3, 4, 9, 20]. The first conserved junction position (JP1) is within the DBD encoding region, the second (JP2) is in the region which encodes the 3' end of the DBD, the third (JP3) is in the LBD motif I (also known as the LBD signature sequence, Tτ) encoding region  and the fourth one (JP4) is in the LBD motif II encoding region . JP1 exhibits diversity among nuclear receptor genes but is conserved for certain gene groups including vertebrate TR genes . JP2 is highly conserved usually in the region encoding the fourth amino acid of the hinge region, but most of the Drosophila NR genes have lost this junction . JP3 is NR gene-specific  and JP4 is at the same position in all vertebrate NR genes . All identified non-chordate TR orthologues possess the four JPs, except JP1 was lost in the schistsome TR genes (SmTRα, SmTRβ, SjTRα and SjTRβ) and JP2 was lost in DpTR. JP3 of all invertebrate TR genes is at the same position as that found in vertebrate TR genes (NR1A) (Fig. 7). JP4 in all invertebrate TRs is located at the same conserved position as found in all analyzed NRs . The conserved splice junctions in these invertebrate TRs suggest that the organization of the TR genes was ancient and has been maintained throughout evolution. It adds further support to the notion that vertebrate and invertebrate TR orthologues originated from a common ancestor gene.
One property of vertebrate TRs is the ability to form a heterodimer with RXR . To determine if a non-chordate TR would form a heterodimer with RXR, the interaction between SmRXR1 with SmTRα and SmTRβ was tested by GST pull-down experiments. Either GST-SmTRα-LBD or GST-SmTRβ-LBD fusion protein was incubated with in vitro translated 35S-SmRXR1. GST-bound resin was used as a negative control to assess non-specific background binding. The results showed that both SmTRα-LBD and SmTRβ-LBD could form a heterodimer with SmRXR1 in vitro. The results demonstrate that non-chordate TRs are capable of forming a heterodimer with RXR (Fig. 8).
Electrphoretic mobility shift assays (EMSAs) were performed to determine DNA binding specificity of SmTRα and SmTRβ. A DNA element containing the half-site AGGTCA, a direct repeat of the half-site spaced with 0–5 nucleic acids (DR0-DR5) and palindrome repeat of the half-site not separated by nucleic acids (Pal0) were employed. A gel shift was observed when γ-32P labelled half-site, DR0-DR5 or Pal0 were added to SmTRα. SmTRα could bind to all the tested oligonucleotides either as a monomer or as a homodimer (Fig. 9A). SmTRβ demonstrated a weak gel shift with the same labelled oligonucleotides as template (Fig. 9B). Heterodimers of either SmTRα or SmTRβ with SmRXR1 did not cause a gel shift (Fig. 9A and 9B). SmRXR1 alone could bind to the tested oligonucleotides (Fig. 9A and 9B) as we previously demonstrated .
The above result showed that SmTRα could bind to a half site, DR0-DR5 and Pal0 element as a monomer or a homodimer. To test whether the DBD domain (which contains a weaker-dimer interface in the D-box) could bind to the oligonucleotides as a homodimer, in vitro synthesized SmTRα (Gln182-Ala288) (containing 20 aa at the 5' end of the DBD, the DBD and 40 aa at the 3' end of the DBD) was tested. SmTRα (Gln182-Ala288) bound to all tested oligonucleotides as either a monomer or a homodimer (Fig. 10). The results suggest that there is a dimer interface located in DBD region of SmTRα [25–27]. Specificity was demonstrated by the ability of cold specific competitor but not the non-specific competitor to prevent SmTRα binding to the oligonucleotide template.
An extensive search of whole genomic sequence (WGS) databases extracted from NCBI and online resources revealed the presence of TR orthologues in the Platyhelminths, Mollusca, Crustacea, and Echinodermata. However, no TR orthologues nor sequences encoding a TR P-box sequence of CEGCKG followed by the amino acid sequence FFRR (CEGCKGFFRR), which is unique to NR subfamily 1 including TR, was found in sponges (Reniera sp) nor cnidarians (Hydra magnipapillata and Nematostella vectensis). The results suggested that the TR orthologue gene originated in a common ancestor of the Bilateria. The gene duplication of TRs was deduced from a phylogenetic analysis employing DBD and LBD sequences. The analyses demonstrated that the TR orthologue genes underwent duplication independently in invertebrates and vertebrates. Furthermore, the TR orthologues in the Platyhelminths underwent duplication after the split of the turbellaria and trematodes. In vertebrates, the TR genes underwent a duplication event in lampreys, a jawless vertebrate, after the lamprey split from a common ancestor of jawed vertebrates. This result is consistent with the duplication of Hox gene clusters in lampreys . Although the A/B domain of NR is highly divergent, we identified a conserved NTSS (N-terminal signature sequence, (G/P)YIPSY(M/L)XXXGPE(D/E)) located in the C-terminus of the A/B domain of TRs. Analysis of other NRs (data not shown) demonstrates this motif is found only in TRs suggesting it is TR-specific. The NTSS of TRs was found to be sequence-specific for the chordate and for the non-chordate TRs, three amino acids present in non-chordate TRs are missing in chordate NTSS. In the DBD, a methionine at position 16 of the N-terminus was non-chordate-specific, while two amino acids inserted at positions 33–34 of the N-terminus were chordate-specific, These differences which are conserved, support the notion that TR genes duplicated independently in non-chordate and chordate species.
That the JP3 conserved in all invertebrate TRs are also conserved in vertebrate TR genes support the notion that the vertebrate and invertebrate TR orthologue genes originated from a common ancestor gene. Further, that the four JPs are conserved in both invertebrate and vertebrate TR orthologue genes suggest that the gene organization of vertebrate TR genes is ancient and has been maintained through out evolution of TRs.
GST pull-down experiments showed that SmTRs exhibits similarity to vertebrate TR proteins as they both can form a heterodimer with RXR (SmRXR1). RXRs have been characterized in other invertebrates such as Mollusca and Arthropoda [29, 30]. The formation of SmTRs with SmRXR1 indicates that other invertebrate TRs can possibly form a heterodimer with RXR. The ability of vertebrate NR to bind to a DNA element on a target gene was described as the 1–2–3–4–5 rule [31, 32]. Vertebrate TRs are known to bind to half-sites as a monomer, to DR0 as a heterodimer with RXR  and to DR4 as a homodimer or heterodimer [32, 34–43]. The binding of the TR monomer or homodimer is weak, which is due to a rapid dissociation of the TR-TRE complex . Our EMSA show that SmTRα can bind to a half-site as a monomer and bind to DR0-DR5 as a homodimer. The results suggested that TR gained the ability to bind to conserved half-site repeats in a common ancestor, but did not subsequently evolve a distinct specificity as regards spacing between half sites. Although the DBD sequence between SmTRα and SmTRβ are highly similar, SmTRβ showed weak binding to the DNA elements compared to that of SmTRα. This suggests that other domains may have a role in determining the ability of SmTRβ to bind to the DNA element, such as has been shown for another schistosome nuclear receptor, SmNR1 . Although both SmTRs could form a heterodimer with SmRXR1 as determined in a GST pull-down assay, SmTRα/SmRXR1 or SmTRβ/SmRXR1 did not bind to the DNA elements tested. This may due to the flanking or spacing of nucleotides between the HRE, as it is known that these nucleotides also determine the ability of the protein to bind to it. Identification of the HRE in a target gene will help to understand the evolution of TR DNA binding properties. Mammalian TRs are ligand-dependent transcription factors, they repress basal transcription through association with a variety of corepressors in the absence of T3. A structural change occurs when T3 binds to the TR that results in the release of natural corepressors and results in the recruitment of coactivators to activate gene transcription . Determination of whether invertebrate TR can bind to T3 will help to understand the evolution of the thyroid hormone pathway.
We have identified ortholouges of TR in various invertebrate phyla. Sequence analysis demonstrates that there is a 'N-terminal signature sequence' (NTSS) located in the A/B domain -(N-terminal region) specific to TRs. Phylogenetic analysis of conserved DBD sequences showed that TR orthologues originated from a common ancestor of the Bilateria which underwent gene duplication independently in the Protostomia and Deuterostomia. Using the schistosome TR orthologues as representative of the invertebrate TR genes, we have demonstrated that they can form heterodimers with RXR and bind as monomers and homodimers to consensus vertebrate DNA half sites and the direct repeat of the half sites DR0 to DR5.
cDNAs containing the entire open reading frame (ORF) of two S. mansoni TRs (SmTRα and SmTRβ) were isolated by a PCR strategy using a S. mansoni paried adult worm cDNA library (pBluescript SK (+/-) phagemid) pool as template DNA. The PCR primers for one end (either 5' or 3' end) were designed according to a fragment encoding the previously identified DBD region of each gene . The primer for the other end (either 5' or 3' end) was a vector universal primer (M13-Rev and T3, or M13-For and T7 primers). PCR products were separated on 1.2% agarose gels, ligated into pCR2.1 TOPO vector (Invitrogen) and sequenced. After the correct fragments were identified, the cDNA sequence containing the 5' UTR, ORF and 3' UTR were obtained by PCR. Each cDNA was shown to be related to a single mRNA species by sequencing the PCR products obtained from single-stranded cDNA using primers located within the 5'UTR and 3'UTR of each gene.
Whole genomic sequences (WGS) were extracted from the GenBank public ftp site  and imported into the StarBlast program (DNASTAR) to build a local database  which was searched by tblastn using the sequence of the DBD of SmTRα as the query. Any sequence that contained a zinc finger structure of the DBD (Cys-X2-Cys-X13-Cys-X2-Cys or Cys-X5-Cys-X9-Cys-X2-Cys) was retained. Sequence walking was carried out to assemble the contigs. Website databases of GenBank (nr, EST_human, EST_mouse and EST_other databases) Genbank blast , European Bioinformatics Institute  and Swiss-Prot  were also mined by tblasn or blastp using the same query sequence as above .
BAC library screening and BAC DNA sequencing
S. mansoni BAC clones containing SmTRα and SmTRβ were identified by screening S. mansoni SmBAC1 or CHOR-1 BAC library with methods previously described . For BAC DNA sequencing, the BAC clone was grown in 100 ml LB medium (12.5 μg/ml choramphenicol) and the BAC DNA was purified using Plasmid Midi kit (Qiagen) and sequenced on an ABI-377 automatic sequencing machine.
Phylogenetic tree construction
Phylogenetic trees were constructed from deduced sequences of the DBD and the sequence of DBD plus LBD, respectively. Sequences were aligned with ClustalW . Phylogenetic analysis of the data set was carried out using Maximum Likelihood method under Jones-Taylor-Thornton (JTT) substitution model  with a gamma distribution of rates between sites (eight categories, parameter alpha, estimated by the program) using PHYML (v2.4.4))  with support values obtained by bootstrapping a 100 replicates. The same data set was also tested by Bayesian inference using MrBAYES v3.1.1 with a JTT amino acid replacement model + gamma rates . The trees were started randomly; four simultaneous Markov chains were run for 5 million generations for the DBD data set and 3 million generations for the DBD+LBD data set, respectively. The trees were sampled every 100 generations, the Bayesian posterior probabilities (PPs) were calculated using a Markov chain Monte Carlo (MCMC) sampling approach implemented in MrBAYES v3.1.1, with a burn-in value setting at 12,500 for DBD data set and 7,500 generations for the DBD+LBD data set, respectively.
cDNAs encoding part of the hinge region plus the LBD domain of SmTRα (Val739-Glu1115) and SmTRβ (Pro487-Cys847) were inserted into pGEX-4T-1 vectors to form pGEX-SmTRα(LBD) and pGEX-SmTRβ(LBD). cDNA encoding SmRXR1 was previously inserted into pCITE-4a vector to form pCITE-SmRXR1 . E. coli AD 494 (DE3) pLys S competent cells (Novagen) were transformed with pGEX-SmTRα(LBD) or pGEX-SmTRβ(LBD) and the GST fusion proteins were purified by passage over a glutathione-Sepharose column according to standard protocols. pCITE-SmRXR1 was transcribed and translated using the Single Tube Protein System (Novagen) following the manufacture's protocol to produce 35S labeled protein. For pull-down experiments, a 50 μl reaction that contained 2 μl of the in vitro translation reaction, SmTRα(LBD) or SmTRβ(LBD) GST fusion protein, or GST protein (negative control) affixed to glutathione-Sepharose beads (about 2 μg) and binding buffer (50 mM Tris-HCl, pH 7.5, 100 mM NaCl, 10% glycerol, 0.15% Nonidet P40) was used . The reaction was incubated overnight at 4°C and then washed three times with binding buffer. The bound proteins were analyzed by 10% SDS-PAGE and autoradiography.
Electrphoretic mobility shift assay (EMSAs)
cDNA encoding SmTRα, SmTRβ and SmTRα(Gln182-Ala288) (containing 20 aa at the 5' end of the DBD, DBD and 40 aa at the 3' end of the DBD) was inserted into pCITE-4a vector to form pCITE-SmTRα, pCITE-SmTRβ and pCITE-SmTRα(Gln182-Ala288). The protein was produced in vitro using TNT quick coupled transcription/translation system (Promega). The following complementary single-stranded oligonucleotides containing consensus half-sites AGGTCA were synthesized:
half-site: 5'-GTACCGTAAGGTCA CTCGCGT-3',
DR0: 5'-CCGTAAGGTCAAGGTCA CTCG-3',
DR1: 5'-CCGTAAGGTCA CAGGTCA CTCG-3',
DR2: 5'-CCGTAAGGTCA CAAGGTCA CTCG-3',
DR3: 5'-CCGTAAGGTCA CAGAGGTCA CTCG-3',
DR4: 5'-CCGTAAGGTCA CAGGAGGTCA CTCG-3',
DR5: 5'-CCGTAAGGTCA CCAGGAGGTCA CTCG-3'.
Pal0: 5'-CGCAAGGTCATGACCT CG-3' . One strand of each oligonucleotide was annealed after incubation at 100°C for 3 minutes to its complementary oligonucleotide and then labeled with T4 polynucleotide kinase and γ-32P adenosine triphosphate. The binding reactions were incubated on ice for 40 minutes in 15 μl reaction mixture containing 40,000 cpm of probe, 3 μl in vitro translation reaction, 3 μl 5 × buffer (20% glycerol, 5 mM MgCl2, 2.5 mM EDTA, 2.5 mM DTT, 250 mM NaCl, 50 mM Tris-HCl (pH 7.5) and 0.25 mg/ml poly(dI-dC)poly (dI-dC)). The reaction was separated on 6% (v/v) native polyacrylamide gel containing 2.5% glycerol in 1 × TBE buffer at 4°C . Gels were dried, exposed to x-ray film and autoradiographed.
Bacterial artificial chromosome
Ligand binding domain
N-terminal signature sequence
Thyroid hormone receptor
Whole genomic sequence.
Moras D, Gronemeyer H: The nuclear receptor ligand-binding domain: structure and function. Curr Opin Cell Biol. 1998, 10 (3): 384-391. 10.1016/S0955-0674(98)80015-X.
Baniahmad A, Kohne AC, Renkawitz R: A transferable silencing domain is present in the thyroid hormone receptor, in the v-erbA oncogene product and in the retinoic acid receptor. The EMBO journal. 1992, 11 (3): 1015-1023.
Gronemeyer H, Laudet V: Transcription factors 3: nuclear receptors. Protein Profile. 1995, 2 (11): 1173-1308.
Laudet V, Gronemeyer H: The Nuclear Receptor Factsbook. 2001, Academic Press
Adams MD, Celniker SE, Holt RA, Evans CA, Gocayne JD, Amanatides PG, Scherer SE, Li PW, Hoskins RA, Galle RF: The genome sequence of Drosophila melanogaster. Science. 2000, 287 (5461): 2185-2195. 10.1126/science.287.5461.2185.
Bertrand S, Brunet FG, Escriva H, Parmentier G, Laudet V, Robinson-Rechavi M: Evolutionary genomics of nuclear receptors: from twenty-five ancestral genes to derived endocrine systems. Mol Biol Evol. 2004, 21 (10): 1923-1937. 10.1093/molbev/msh200.
Holt RA, Subramanian GM, Halpern A, Sutton GG, Charlab R, Nusskern DR, Wincker P, Clark AG, Ribeiro JM, Wides R: The genome sequence of the malaria mosquito Anopheles gambiae. Science. 2002, 298 (5591): 129-149. 10.1126/science.1076181.
Sluder AE, Maina CV: Nuclear receptors in nematodes: themes and variations. Trends Genet. 2001, 17 (4): 206-213. 10.1016/S0168-9525(01)02242-9.
Wu W, Niles EG, El-Sayed N, Berriman M, LoVerde PT: Schistosoma mansoni (Platyhelminthes, Trematoda) nuclear receptors: Sixteen new members and a novel subfamily. Gene. 2006, 366 (2): 303-315. 10.1016/j.gene.2005.09.013.
Verjovski-Almeida S, DeMarco R, Martins EA, Guimaraes PE, Ojopi EP, Paquola AC, Piazza JP, Nishiyama MY, Kitajima JP, Adamson RE: Transcriptome analysis of the acoelomate human parasite Schistosoma mansoni. Nat Genet. 2003, 35 (2): 148-157. 10.1038/ng1237.
de Mendonca RL, Bouton D, Bertin B, Escriva H, Noel C, Vanacker JM, Cornette J, Laudet V, Pierce RJ: A functionally conserved member of the FTZ-F1 nuclear receptor family from Schistosoma mansoni. Eur J Biochem. 2002, 269 (22): 5700-5711. 10.1046/j.1432-1033.2002.03287.x.
de Mendonca RL, Escriva H, Bouton D, Zelus D, Vanacker JM, Bonnelye E, Cornette J, Pierce RJ, Laudet V: Structural and functional divergence of a nuclear receptor of the RXR family from the trematode parasite Schistosoma mansoni. Eur J Biochem. 2000, 267 (11): 3208-3219. 10.1046/j.1432-1327.2000.01344.x.
Freebern WJ, Niles EG, LoVerde PT: RXR-2, a member of the retinoid x receptor family in Schistosoma mansoni. Gene. 1999, 233 (1–2): 33-38. 10.1016/S0378-1119(99)00161-4.
Freebern WJ, Osman A, Niles EG, Christen L, LoVerde PT: Identification of a cDNA encoding a retinoid X receptor homologue from Schistosoma mansoni. Evidence for a role in female-specific gene expression. J Biol Chem. 1999, 274 (8): 4577-4585. 10.1074/jbc.274.8.4577.
Hu R, Wu W, Niles EG, LoVerde PT: Isolation and characterization of Schistosoma mansoni constitutive androstane receptor. Mol Biochem Parasitol. 2006, 148 (1): 31-43. 10.1016/j.molbiopara.2006.02.017.
Hu R, Wu W, Niles EG, Loverde PT: SmTR2/4, a Schistosoma mansoni homologue of TR2/TR4 orphan nuclear receptor. Int J Parasitol. 2006, 36: 1113-1122. 10.1016/j.ijpara.2006.06.003.
Lu C, Wu W, Niles EG, Loverde PT: Identification and characterization of a novel fushi tarazu factor 1 (FTZ-F1) nuclear receptor in Schistosoma mansoni. Mol Biochem Parasitol. 2006, 150 (1): 25-36. 10.1016/j.molbiopara.2006.06.005.
Wu W, Niles EG, Hirai H, LoVerde PT: Evolution of a novel subfamily of nuclear receptors with members that each contain two DNA binding domains. BMC evolutionary biology. 2007, 7: 27-10.1186/1471-2148-7-27.
Wu W, Niles EG, Hirai H, LoVerde PT: Identification and characterization of a nuclear receptor subfamily I member in the Platyhelminth Schistosoma mansoni (SmNR1). Febs J. 2007, 274 (2): 390-405. 10.1111/j.1742-4658.2006.05587.x.
Zhang Z, Burch PE, Cooney AJ, Lanz RB, Pereira FA, Wu J, Gibbs RA, Weinstock G, Wheeler DA: Genomic analysis of the nuclear receptor family: new insights into structure, regulation, and evolution from the rat genome. Genome Res. 2004, 14 (4): 580-590. 10.1101/gr.2160004.
Nagy L, Kao HY, Love JD, Li C, Banayo E, Gooch JT, Krishna V, Chatterjee K, Evans RM, Schwabe JW: Mechanism of corepressor binding and release from nuclear hormone receptors. Genes Dev. 1999, 13 (24): 3209-3216. 10.1101/gad.13.24.3209.
Perissi V, Staszewski LM, McInerney EM, Kurokawa R, Krones A, Rose DW, Lambert MH, Milburn MV, Glass CK, Rosenfeld MG: Molecular determinants of nuclear receptor-corepressor interaction. Genes Dev. 1999, 13 (24): 3198-3208. 10.1101/gad.13.24.3198.
Le Paslier MC, Pierce RJ, Merlin F, Hirai H, Wu WJ, Williams DL, Johnston D, LoVerde PT, Le Paslier D: Construction and characterization of a Schistosoma mansoni bacterial artificial chromosome library. Genomics. 2000, 65 (2): 87-94. 10.1006/geno.2000.6147.
WTSI S. mansoni WGS database. [http://www.sanger.ac.uk/Projects/S_mansoni/]
IJpenberg A, Jeannin E, Wahli W, Desvergne B: Polarity and specific sequence requirements of peroxisome proliferator-activated receptor (PPAR)/retinoid X receptor heterodimer binding to DNA. A functional analysis of the malic enzyme gene PPAR response element. J Biol Chem. 1997, 272 (32): 20108-20117. 10.1074/jbc.272.32.20108.
Predki PF, Zamble D, Sarkar B, Giguere V: Ordered binding of retinoic acid and retinoid-X receptors to asymmetric response elements involves determinants adjacent to the DNA-binding domain. Mol Endocrinol. 1994, 8 (1): 31-39. 10.1210/me.8.1.31.
Wilson TE, Paulsen RE, Padgett KA, Milbrandt J: Participation of non-zinc finger residues in DNA binding by two nuclear orphan receptors. Science. 1992, 256 (5053): 107-110. 10.1126/science.1314418.
Fried C, Prohaska SJ, Stadler PF: Independent Hox-cluster duplications in lampreys. Journal of experimental zoology Part B. 2003, 299 (1): 18-25. 10.1002/jez.b.37.
Guo X, Xu Q, Harmon MA, Jin X, Laudet V, Mangelsdorf DJ, Palmer MJ: Isolation of two functional retinoid X receptor subtypes from the Ixodid tick, Amblyomma americanum (L.). Mol Cell Endocrinol. 1998, 139 (1–2): 45-60. 10.1016/S0303-7207(98)00073-2.
Oro AE, McKeown M, Evans RM: Relationship between the product of the Drosophila ultraspiracle locus and the vertebrate retinoid X receptor. Nature. 1990, 347 (6290): 298-301. 10.1038/347298a0.
Shulemovich K, Dimaculangan DD, Katz D, Lazar MA: DNA bending by thyroid hormone receptor: influence of half-site spacing and RXR. Nucleic Acids Res. 1995, 23 (5): 811-818. 10.1093/nar/23.5.811.
Umesono K, Murakami KK, Thompson CC, Evans RM: Direct repeats as selective response elements for the thyroid hormone, retinoic acid, and vitamin D3 receptors. Cell. 1991, 65 (7): 1255-1266. 10.1016/0092-8674(91)90020-Y.
Shin DJ, Plateroti M, Samarut J, Osborne TF: Two uniquely arranged thyroid hormone response elements in the far upstream 5' flanking region confer direct thyroid hormone regulation to the murine cholesterol 7alpha hydroxylase gene. Nucleic Acids Res. 2006, 34 (14): 3853-3861. 10.1093/nar/gkl506.
Andersson ML, Nordstrom K, Demczuk S, Harbers M, Vennstrom B: Thyroid hormone alters the DNA binding properties of chicken thyroid hormone receptors alpha and beta. Nucleic Acids Res. 1992, 20 (18): 4803-4810. 10.1093/nar/20.18.4803.
Baniahmad A, Steiner C, Kohne AC, Renkawitz R: Modular structure of a chicken lysozyme silencer: involvement of an unusual thyroid hormone receptor binding site. Cell. 1990, 61 (3): 505-514. 10.1016/0092-8674(90)90532-J.
Glass CK, Holloway JM, Devary OV, Rosenfeld MG: The thyroid hormone receptor binds with opposite transcriptional effects to a common sequence motif in thyroid hormone and estrogen response elements. Cell. 1988, 54 (3): 313-323. 10.1016/0092-8674(88)90194-8.
Kurokawa R, Yu VC, Naar A, Kyakumoto S, Han Z, Silverman S, Rosenfeld MG, Glass CK: Differential orientations of the DNA-binding domain and carboxy-terminal dimerization interface regulate binding site selection by nuclear receptor heterodimers. Genes Dev. 1993, 7 (7B): 1423-1435. 10.1101/gad.7.7b.1423.
Miyamoto T, Suzuki S, DeGroot LJ: High affinity and specificity of dimeric binding of thyroid hormone receptors to DNA and their ligand-dependent dissociation. Mol Endocrinol. 1993, 7 (2): 224-231. 10.1210/me.7.2.224.
Naar AM, Boutin JM, Lipkin SM, Yu VC, Holloway JM, Glass CK, Rosenfeld MG: The orientation and spacing of core DNA-binding motifs dictate selective transcriptional responses to three nuclear receptors. Cell. 1991, 65 (7): 1267-1279. 10.1016/0092-8674(91)90021-P.
Ribeiro RC, Kushner PJ, Apriletti JW, West BL, Baxter JD: Thyroid hormone alters in vitro DNA binding of monomers and dimers of thyroid hormone receptors. Mol Endocrinol. 1992, 6 (7): 1142-1152. 10.1210/me.6.7.1142.
Sap J, Munoz A, Schmitt J, Stunnenberg H, Vennstrom B: Repression of transcription mediated at a thyroid hormone response element by the v-erb-A oncogene product. Nature. 1989, 340 (6230): 242-244. 10.1038/340242a0.
Umesono K, Evans RM: Determinants of target gene specificity for steroid/thyroid hormone receptors. Cell. 1989, 57 (7): 1139-1146. 10.1016/0092-8674(89)90051-2.
Yen PM, Sugawara A, Chin WW: Triiodothyronine (T3) differentially affects T3-receptor/retinoic acid receptor and T3-receptor/retinoid X receptor heterodimer binding to DNA. J Biol Chem. 1992, 267 (32): 23248-23252.
Desvergne B: How do thyroid hormone receptors bind to structurally diverse response elements?. Mol Cell Endocrinol. 1994, 100 (1–2): 125-131. 10.1016/0303-7207(94)90291-7.
GenBank public ftp site. [ftp://ftp.ncbi.nlm.nih.gov/pub/TraceDB/]
European Bioinformatics Institute. [http://www.ebi.ac.uk/blast2]
Jones DT, Taylor WR, Thornton JM: The rapid generation of mutation data matrices from protein sequences. Comput Appl Biosci. 1992, 8 (3): 275-282.
Guindon S, Gascuel O: A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol. 2003, 52 (5): 696-704. 10.1080/10635150390235520.
Huelsenbeck JP, Ronquist F: MRBAYES: Bayesian inference of phylogenetic trees. Bioinformatics (Oxford, England). 2001, 17 (8): 754-755. 10.1093/bioinformatics/17.8.754.
Wagner RL, Apriletti JW, McGrath ME, West BL, Baxter JD, Fletterick RJ: A structural role for hormone in the thyroid hormone receptor. Nature. 1995, 378 (6558): 690-697. 10.1038/378690a0.
This research was supported by NIH grants AI046762 and D43TW006580. We thank the genome sequencing centers for making their data publicly available.
WW performed the experiments and analysis. He contributed to the writing and is the primary author of the manuscript. EGN was responsible for the design and analysis of the experiments. He contributed to the writing of the manuscript. PTL was responsible for the overall design, analysis and interpretation of the results. He contributed to the writing and preparation of the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Wu, W., Niles, E.G. & LoVerde, P.T. Thyroid hormone receptor orthologues from invertebrate species with emphasis on Schistosoma mansoni. BMC Evol Biol 7, 150 (2007) doi:10.1186/1471-2148-7-150
- Bacterial Artificial Chromosome
- Bacterial Artificial Chromosome Clone
- Bacterial Artificial Chromosome Library
- Thyroid Hormone Receptor
- Schistosoma Mansoni