Genome-wide analysis of Pax8 binding provides new insights into thyroid functions
© Ruiz-Llorente et al.; licensee BioMed Central Ltd. 2012
Received: 3 January 2012
Accepted: 24 April 2012
Published: 24 April 2012
The transcription factor Pax8 is essential for the differentiation of thyroid cells. However, there are few data on genes transcriptionally regulated by Pax8 other than thyroid-related genes. To better understand the role of Pax8 in the biology of thyroid cells, we obtained transcriptional profiles of Pax8-silenced PCCl3 thyroid cells using whole genome expression arrays and integrated these signals with global cis-regulatory sequencing studies performed by ChIP-Seq analysis
Exhaustive analysis of Pax8 immunoprecipitated peaks demonstrated preferential binding to intragenic regions and CpG-enriched islands, which suggests a role of Pax8 in transcriptional regulation of orphan CpG regions. In addition, ChIP-Seq allowed us to identify Pax8 partners, including proteins involved in tertiary DNA structure (CTCF) and chromatin remodeling (Sp1), and these direct transcriptional interactions were confirmed in vivo. Moreover, both factors modulate Pax8-dependent transcriptional activation of the sodium iodide symporter (Nis) gene promoter. We ultimately combined putative and novel Pax8 binding sites with actual target gene expression regulation to define Pax8-dependent genes. Functional classification suggests that Pax8-regulated genes may be directly involved in important processes of thyroid cell function such as cell proliferation and differentiation, apoptosis, cell polarity, motion and adhesion, and a plethora of DNA/protein-related processes.
Our study provides novel insights into the role of Pax8 in thyroid biology, exerted through transcriptional regulation of important genes involved in critical thyrocyte processes. In addition, we found new transcriptional partners of Pax8, which functionally cooperate with Pax8 in the regulation of thyroid gene transcription. Besides, our data demonstrate preferential location of Pax8 in non-promoter CpG regions. These data point to an orphan CpG island-mediated mechanism that represents a novel role of Pax8 in the transcriptional output of the thyrocyte.
KeywordsPax8 ChIP-Seq Expression arrays CpG island CTCF SP1
Gene regulation has been the subject of intense investigation over the past decades, mainly focusing on detailed characterization of a particular gene or gene family. However, genome-wide mapping of protein-DNA interactions and epigenetic marks is essential for a full understanding of transcriptional regulation. A precise map of binding sites for transcription factors (TFs), core transcriptional machinery, and other DNA-binding proteins is necessary to decode the gene regulatory networks and their contribution to developmental processes and human disease . In fact, regulation of gene expression by TFs is one of the major mechanisms for controlling cell proliferation, differentiation, and function.
To elucidate the mechanism(s) operating in the establishment and maintenance of cell-specific differentiation, we used thyroid epithelial cells as a model system. These cells are the largest cell population of the thyroid gland and express different TFs called Nkx2.1, Foxe1, Hhex and Pax8, which define the thyroid differentiated phenotype [2, 3]. It is well known that these factors bind to the promoter regions of thyroid-specific genes, such as the genes encoding Thyroglobulin (Tg), Thyroperoxidase (Tpo), and the Sodium Iodide Symporter (Nis), thus regulating their expression. Nevertheless, despite the key relevance of these TFs for thyroid biology, few studies have described additional loci that are transcriptionally regulated by the above mentioned TFs, nor have sequences been described to which these factors bind in enhancers, silencers, or boundary elements that could potentially regulate the transcription of genes over large distances.
Among these thyroid TFs, Pax8 is a member of the paired box-containing proteins and is expressed in the thyroid and kidney, and in the central nervous system during development . It plays an essential role in the differentiation of thyroid cells and, according to the phenotype of Pax8 knockout mice, it seems to be responsible for the formation of the follicles of polarized epithelial thyroid cells . Also, the association between mutations of PAX8 and congenital hypothyroidism in humans underlines an important function of this transcription factor in thyroid pathologies . In order to better understand its role in the maintenance of thyroid function, we explored the transcriptional profile of Pax8-silenced thyroid cells, and integrated these signals with global cis-regulatory sequencing studies (chromatin immunoprecipitation followed by sequencing; ChIP-Seq).
The ChIP-Seq strategy allowed us to identify a large number of novel in vivo Pax8 binding sites that were significantly associated with CpG islands or high GC content sequences. Interestingly, immunoprecipitated peaks were mainly located along intronic regions and grouped in distal positions with respect to transcriptional start sites. Consensus sequence screening of these areas suggested Pax8 interaction with several core transcriptional elements (motif ten element, Inr, and BRE), transcription factors belonging to the AP1 family, and trans-elements factors involved in high order chromatin structure (CTCF) and remodeling (Sp1). Co-immunoprecipitation and reporter assays demonstrated both physical binding and transcriptional cooperation between CTCF/Sp1 and Pax8. Combining sequencing and expression array data, we ultimately provided insights into Pax8- transcriptional networks in the differentiated thyroid that predict its involvement in relevant biological processes and pathways.
Genomic features associated with Pax8 binding sites
In order to identify the genome-wide binding patterns of Pax8 in differentiated thyroid cells, we performed ChIP-Seq in PCCl3 rat thyrocytes using IP and non-IP conditions. Prior to massive sequencing, both conditions were interrogated to verify Pax8 binding site enrichment by means of semi-quantitative PCR (Additional file 1). Using this approach, we confirmed DNA immunoprecipitation of Pax8 binding regions in the rat Nis and Tpo promoters (Additional file 1A), as previously described [7–9]. Therefore, we considered both IP and non-IP conditions as useful samples to further identify whole genome Pax8 binding sites by means of high throughput sequencing technology. After sequencing analysis, we obtained 11,613,355 and 12,125,758 raw reads for control and IP conditions, respectively. Of these, 6,714,002 (57.8%) and 6,431,519 (53.0%) fulfilled the ≤2 mismatches quality filter.
To further localize regions of Pax8 enrichment, we identified Pax8 peaks genome-wide. Peak detection analysis using MACS defined 13,151 Pax8-enriched regions with an average length of 681 bp (Additional file 2). Visual inspection of the Pax8 binding sites and the profiling data in a genome browser for well-known Pax8 targets like Nis, Tpo, [7, 9] and WT1 (Wilms' tumour gene 1) , showed Pax8 binding sites close to the 5'-UTRs of these genes as previously described. A detailed analysis of Nis (Slc5a5, Na+/I− symporter), whose transcription status is tightly regulated by Pax8 , showed a significant Pax8 binding site overlapping with the Nis upstream enhancer (Additional file 1B). These findings clearly validated ChIP-Seq as an efficient and powerful technique for mapping Pax8 binding sites in PCCl3 cells.
Pax8 immunoprecipitated regions delineate Pax8 consensus core sequence in vivo
Immunoprecipitation data reveals interaction of Pax8 with various TFs
Main DNA binding motifs overrepresented in Pax8-dependent peaks
Matches in Input (n)
Transcription factor II B (TFIIB) recognition element
Human motif ten element
X gene core promoter element 1
Drosophila motif ten element
Initiator (INR) and downstream promoter element (DPE) with strictly maintained spacing
Nuclear respiratory factor 1
GC-Box factors SP1/GC; Stimulating protein 1, ubiquitous zinc finger transcription factor
GC-Box factors SP1/GC; Stimulating protein 1, ubiquitous zinc finger transcription factor
Zinc finger / POZ domain transcription factor
Expression arrays analysis identifies a wide set of loci regulated by Pax8
We used whole genome expression arrays to identify Pax8-regulated genes by comparing expression profiles of Pax8-silenced PCCl3 cells with both scrambled siRNA-treated and wild type (wt) PCCl3 cells. This last condition was included to consistently integrate both expression array signals and global cis-regulatory sequencing studies into the same experimental conditions. Misinterpretation of expression data due to compensatory effects via Pax8-related paralogues (Pax2 and Pax5) is ruled out, given that both transcription factors are not expressed in thyroid cells.
Regarding the comparison of siPax8-PCCl3 vs. wt PCCl3, 3,035 and 3,354 probes were down and up-regulated in the Pax8-silenced condition, respectively (Additional file 5). A lower number of significant probes was detected for siPax8-PCCl3 vs. siScramble-PCCl3 (797 and 777 probes were down and up-regulated in the Pax8-silenced condition, respectively) (Additional file 5). Statistically significantly differently expressed probes (adjusted p-values <0.005) for both comparisons included 633 down and 565 up-regulated targets (Additional file 5), which represent a set of 849 loci.
Pax8 is involved in controlling key cellular events
FatiScan gene set enrichment analysis
Gene Ontology term
WT vs Pax8
Scramble vs Pax8
Response to external stimulus (GO:0009605)
Response to wounding (GO:0009611)
Cellular component movement (GO:0006928)
Response to hormone stimulus (GO:0009725)
Immune response (GO:0006955)
Cell adhesion (GO:0007155)
Response to steroid hormone stimulus (GO:0048545)
Antigen processing and presentation (GO:0019882)
Cell migration (GO:0016477)
Single functional analysis
We additionally used the FatiGO in silico tool to extract Gene Ontology (GO) terms overrepresented in our down- and up-regulated set of differentially expressed genes. Considering down-regulated probes for each comparison, we observed an enrichment in biological processes related to a wide variety of DNA, RNA, and protein processes (purine and pyrimidine metabolism, response to DNA damage, DNA replication, nucleotide and base exchange repair, mismatch repair and homologous recombination, RNA degradation, and amino acid metabolism), cell response to chemical and stress stimuli, immune response, and p53 and insulin-related pathways (phosphatidyl inositol system and metabolism) (Additional file 9). Concerning GO terms enriched amongst up-regulated probes, it is worth to mention the over-representation of genes involved in biological processes such as immune response, cell response to stimuli, apoptosis and cell death, cell motion/migration/adhesion, and regulation of cell differentiation (Additional file 10).
KEGG pathways associated to Pax8 silencing
Scr. vs Pax8
Wt. vs Pax8
Phosphatidylinositol signaling system
Inositol phosphate metab.
Glycine, serine and threonine metab.
Selenoamino acid metab.
Cysteine and methionine metab.
Arginine and proline metab.
Cell adhesion molecules
Antigen processing and presentation
Autoimmune thyroid disease
Base excision repair
Nucleotide excision repair
p53 signaling pathway
Cell adhesion molecules
Cytosolic DNA- sensing pathway
NOD-like receptor signaling pathway
Toll-like receptor signaling pathway
Chemokine signaling pathway
Renal cell carcinoma
Pathways in cancer
Autoimmune thyroid disease
Antigen processing and presentation
Integrated data reveal a reduced percentage of genes transcriptionally regulated through promoter sequences
Independent validation confirms significant findings defined by ChIP-seq and expression arrays
In addition, Pax8 silencing by means of transient transfection of siRNA was significantly associated with decreased expression levels of these potential targets (Figure 7B, upper pannel), thus demonstrating a direct transcriptional effect of Pax8 on these genes. mRNA expression validation was also done for several genes that were upregulated in the absence of Pax8, including the genes encoding: CCL2, a chemokine involved in thyroid autoimmunity ; S100A4, a calcium-binding protein which plays a role in angiogenesis, extracellular matrix remodelling and tumor microenvironment, and reported to be overexpressed in metastatic papillary thyroid microcarcinomas ; SCNN1G and PADI1, which exert a role in Na+ transport and differentiation in epithelial cells, respectively [29, 30] (Figure 7B, lower panel). In silico analysis of significant IP peaks located along promoter areas of these loci demonstrated Pax8 potential binding sites in 3 out of 4 genes (data not shown). Globally, these findings underscore the efficiency and accuracy of ChIP-Seq and expression array technologies to define a Pax8-dependent gene network, which allowed us to identify biological functions of Pax8 in thyroid cells.
Despite the known relevance of the transcription factor Pax8 for adult thyrocyte physiology, few data have been published concerning Pax8 target genes other than key thyroid-related genes (Tg Tpo, and Nis). The transcriptional output of Pax8 during thyroid development is unknown but essential, given that thyroid follicular precursors are not formed in Pax8 null mouse embryos, which ultimately impairs the formation of follicle structures and thyroid hormone biosynthesis .
With regard to its link to tumour development, Pax8 expression decreases or is lost in follicular thyroid carcinomas as well as in oncogene-transformed thyroid cells . Moreover, several well-known tumour suppressors, including TP53 and WT1, have been defined as Pax8 targets, and cytoplasmic Pax8 staining has been positively associated with tumour size, metastasis, local invasion, recurrence, or persistence in the thyroid . Taking into account all these premises, and in order to better understand the role of Pax8 in the maintenance of thyroid function, we decided to explore the transcriptional profile of Pax8-silenced thyroid PCCl3 cells, and to integrate these signals with genome wide cis-regulatory studies. Thus, our experimental design combined putative and novel Pax8 binding sites with analysis of actual target gene expression regulation, a strategy successfully used for identifying direct targets for other transcription factors [35, 36].
Our unbiased mapping of Pax8 binding sites along the rat genome has identified a large number of DNA sequences that are occupied in living thyrocytes. Moreover, this is the first study addressing in vivo genome-wide mapping of Pax8-DNA binding sites, and the Pax8 consensus binding motif here defined encompasses motifs described by previous reports focused either on single gene regulation [7, 12] or on Paired-box DNA motif characterization [13, 14]. The ChIPSeq approach also led to significant immunoprecipitation of genomic sequences containing CpG islands, as well as CpG dinucleotides. Extensive literature has linked the location of CpG islands and GC-enriched regions to transcriptionally permissive chromatin [37, 38], which could lend support to a relevant role of Pax8 in the transcriptional output of the thyrocyte. About half of all CpG islands self-evidently contain TSSs, while the other half (known as “orphan” CpG islands) are either within or between characterized transcription units and have unknown significance [11, 38]. Despite a lack of association to annotated promoters, “orphan” CpG islands have been associated to transcriptional initiation and dynamic expression during development . In agreement with this, we found significant Pax8 binding to orphan CpG islands in intronic regions and a preferential binding to such islands 10–100 kb upstream or downstream of a transcription start site. In fact, genomic studies indicate that almost half of the human coding genes have alternative promoters  and that transcription factor binding sites (TFBSs) in classically defined promoter regions may represent a minority of genomic binding sites . Moreover, this latter report clearly demonstrated an association between TFBSs and the expression of non-coding RNAs, which could be modulating the expression of the gene encoded by the opposite strand. Less directly, a subset of intergenic H3K4me3 peaks, many of which are likely to correspond to orphan CpG islands, were found to represent TSSs for long non-coding RNAs . Our findings suggest that Pax8 binds orphan CpG islands that could represent alternative promoters of nearby annotated genes  or ncRNAs that regulate gene expression.
Otherwise, Pax8-dependent ChIP-Seq data demonstrated an enrichment of genomic regions with overrepresentation of general transcriptional regulatory elements (Human MTE and Drosophila MTE, Inr-DPE and BRE). MTE constitutes a core promoter element (~20-30 nt downstream of the TSS) associated with RNA polymerase II-mediated transcription [44, 45]. Furthermore, human orphan CpG islands have been associated with RNA polymerase II binding sites . On the other hand, Inr-DPE and BRE elements represent functional binding sites for TFIIB and TFIID (transcription initiation factor IIB and IID, respectively), which are main components of the basal transcription machinery . Interestingly, Jin et al recently described synergistic MTE-Inr-BRE transcriptional modules in more than 9,000 orthologous mouse and human genes . Whereas functional experiments should be performed to demonstrate an interaction of Pax8 with these general core elements, our data underscore the importance of synergistic interactions between core promoter elements and tissue-specific TFs to ultimately modulate gene expression.
Potential Pax8 partners in transcriptional regulation
2Apart from the classical view of TFs interacting with promoter regions, TFs could activate gene expression by interacting with common lineage-specific TFs and/or binding to distal regions (enhancers). Synergistic effects of Pax8 and AP1 proteins have been shown to occur in the regulation of Nis transcription through interaction along the NUE element , and AP1 and PAX proteins also interact to cooperate in the modulation of transcription of other genes . Accordingly, we observed an overrepresentation of binding motifs related to NRF-1 (Nuclear respiratory factor-1), and several AP1 members (c-FOS, BATF3, and c-JUN) were differentially expressed in Pax8-deprived thyroid cells. However, no significant findings were obtained for other transcription factors described to act synergistically with Pax8, such as Nkx2-1 and TAZ/WWTR1 proteins , indicating that this cooperative transcriptional role could be restricted to specific loci rather than representing a global transcription phenomenon in thyroid cells.
Functional studies described in the present paper confirmed physical in vivo interactions between Pax8 and CTCF or Sp1 in thyrocytes. These novel partners were further demonstrated to modulate the effect of Pax8 on the transcription of the NIS gene, thus confirming that these interactions are functionally relevant. Evidence has been accumulating concerning the role of CTCF in the establishment of intra-chromosomal loops which ultimately mediate protein-protein contacts between distal complexes and the general transcription machinery [49, 50]. On the other hand, Sp1 is a ubiquitously expressed transactivator, which physically interacts with several components of TFIIB and TFIID (mentioned above as potential Pax8 interacting proteins) and factors related to epigenetic events, such as histone deacetylases and p300/CBP histone acetyltransferase . Interestingly, several studies have described synergistic interactions between Pax8 and p300 acetyltransferase for enhancing the transcriptional activity of thyroid-related genes [52, 53]. Taking into account this complete transcriptional scenario, our data describe potential interactions of Pax8 with both common TFs and core elements, which could cooperate in chromatin remodeling for transcriptional regulation in thyroid cells.
Identification of biological processes controlled by Pax8 in thyroid cells
Pax8 has been mainly associated to thyroid differentiation and development through its transcriptional role in key thyroid-related genes [54, 55]. At this regard, we observed a downregulation of DIO1 after abolishing Pax8 (Additional file 5), which potentially binds to a critical region for selenocystein insertion in the DIO1 mRNA. Data were recently provided indicating that TSH tightly regulates DIO1 expression in thyroid cells through Pax8-dependent DIO1 mRNA stabilization (S.G. Leoni; unpublished observations). Moreover, gene expression profiling in normal versus malignant thyroid tissues demonstrated a downregulation of DIO1 and DIO2 , which could be linked to Pax8 loss during cancer progression.
Intriguingly, Pax8 modulates the expression of several genes involved in carcinogenesis and thyroid malignancies (phosphatidyl-inositol/insulin and MAPK pathways) and cell cycle processes (CDKN2B CCNB1 and CCNB2, among others) (Additional file 11). These findings are in accordance with previous studies in which Pax8 expression was abolished in the differentiated thyroid cell line FRTL5 [20, 57]. Our data would also explain the biological mechanism underlying the partial decrease in thyrocyte proliferation in response to both IGF-I and TSH (main regulators of thyroid proliferation and differentiation) after both Nkx2.1 and Pax8 mRNA silencing .
DNA-related biological processes involved a plethora of functional categories (replication, repair and metabolism), highlighting the novel finding of Brca1-dependence on Pax8. In this regard, Shih et al described that BRCA1 and BRCA2 germline mutations were twice as common in individuals developing a second non-ovarian carcinoma, with follicular thyroid carcinoma being one of the most frequent secondary tumours . This finding can be of great relevance in the development of sporadic thyroid tumors, given that, as mentioned before, Pax8 expression is decreased or lost in thyroid tumours.
Recent reports have associated the transcription factors Pax2 and Pax5 with increased capabilities for cell motility and adhesion in human cancer [58, 59]. In parallel with these Pax-related functions, we observed significant expression changes of loci involved in cell motion/adhesion, notably the Pax8 effect on NCAM1 (neural cell adhesion molecule 1) transcription. NCAM1 and other components of adherens junctions, such as cadherins, have been described to be essential for maintaining cell polarity and epithelial integrity . Interestingly, Cadherin-16 (Cdh16/Ksp-cadherin) was recently proposed to play a TSH-regulated role in thyroid development , and its expression and promoter activity is controlled by Pax8 [20, 62]. We have not only confirmed transcriptional regulation of Cdh16 by PAX8, but also defined additional PAX8-dependent genes that could be essential for thyroid cell polarity (MYO5b and Rab17, among others). In this regard, germline mutations in MYO5b have been associated with disruption of epithelial cell polarity in MVID (MIM251850) . This role is exerted via its involvement in vesicle trafficking through direct interactions with Rab GTPase proteins, such as RAB11a and RAB8a. Further functional studies should be performed to evaluate potential Myo5b interactions with RAB17, another Rab GTPase protein involved in membrane trafficking and confirmed as a Pax8 target in the present study.
State-of-the-art cis-regulatory sequencing studies have been combined with mRNA silencing and expression arrays to further characterize the functional relevance of TF-interacting DNA regions and thus to define their transcriptional output. In our study, we describe Pax8 as a master regulator of key cellular processes for thyrocyte biology, including cell cycle regulation, DNA repair, replication and metabolism, and cell polarity, and define a large set of genes whose expression is modulated by Pax8. However, only a minor fraction (6.4%) of the Pax8 binding sites identified are close to TSSs and correlate with altered mRNA expression, in agreement with studies carried out on other TFSs (1-10%) [35, 36]. This moderate percentage may be explained by the Pax8 binding site distribution, where most of the binding sites are related with orphan CGI regions. In this regard, our study demonstrates Pax8 binding sites in regions distal to TSSs, preferentially in intronic regions, which highlights a potential role as a distal or alternative transcriptional regulator, although this does not rule out indirect regulation. Distal regulation by Pax8 is supported by the interaction with chromatin remodeling factors such as CTCF and Sp1 described in the present study. Therefore, these findings suggest a new function of Pax8 as a chromatin remodeling factor in thyroid follicullar cells, which should be validated and elucidated in future studies.
Cell culture and plasmids
PCCl3 cells are a continuous line of thyroid follicular cells derived from Fischer rats that express the thyroid-specific genes Tg Tpo, and Nis, as well as the thyroid-specific transcription factors Nkx2.1, Foxe1, and Pax8 . They were grown in Coon’s modified Ham’s F-12 medium supplemented with 5% donor calf serum and a six-hormone mixture . For transfection assays, HeLa cells were used and cultured as described .
The 2,854-bp DNA fragment of the rat Nis promoter (pNIS-2.8) which contains the NUE region with two Pax8 binding sites was cloned in our laboratory . Full length Pax8, Sp1, and CTCF were subcloned respectively in pcDNA3.1+, pBS and pcDNA1 Neo, and have been previously described [66–68].
ChIP samples were prepared from PCCl3 cells as follows: cultures of 10 × 106 cells were cross-linked with 1% formaldehyde for 10 minutes at room temperature. Cross-linking was stopped by the addition of glycine to a final concentration of 125 mM, and cells were washed twice with PBS. The cell pellet was resuspended consecutively in ChIP lysis buffers  and sonicated for 90 minutes (30 seconds high frequency pulsing followed by 30 seconds resting) using the Bioruptor sonicator (Diagenode, Denville, NJ) to produce chromatin fragments of 200–500 bp on average. After isolating the sheared chromatin, we incubated it with Pax8 antibody-coated magnetic beads. To prepare these beads, 100 μl of magnetic sheep anti-rabbit IgG beads (Invitrogen, Carlsbad, CA) were incubated overnight with 10 μg polyclonal anti-mouse Pax8 antibody (Biopat, Milan, Italy), that recognizes also rat Pax8 at 4°C. The following day, the beads were rinsed and added to the sheared chromatin and incubated overnight at 4°C. Samples were then rinsed five times with RIPA buffer, and the antibody was stripped from the beads by incubating in 1% SDS at 65°C for 15 minutes; cross-linking was reversed by incubating overnight at 65°C. The next day, samples were sequentially treated with RNAse A and Proteinase K, phenol-chloroform extracted, ethanol precipitated in the presence of 20 μg glycogen, and resuspended in 50 μl 10 mM Tris pH 8.0. Procedure controls included an input condition, obtained before DNA-protein complex sonication and further used during ChIP-Seq assays as normalization sample, and non-immunoprecipitated DNA (non-IP DNA), which was obtained just prior to Pax8 immunoprecipitation.
Before sequencing, Pax8-IP DNA (IP) was used to confirm enrichment of target DNA fragments (Additional file 13) by means of real time-PCR, using as positive IP controls both the Nis upstream enhancer element (NUE) and Tpo promoter sequences [7, 9]. Negative controls of Pax8 binding to genomic DNA included promoter areas of Gad1 (glutamate decarboxylase 1) and Afm (afamin or alpha-albumin), and a region of the Nis locus that does not bind Pax8. PCR reactions were assembled in triplicate with SYBR Green ER qPCR Supermix (Invitrogen, Carlsbad, CA) and run on an Applied Biosystem 7500 Real Time PCR system. The enrichment of target sequences in ChIP material was calculated relative to the Afm negative control, and normalized to their relative amplification in non-IP DNA.
Illumina high-throughput sequencing
After verifying Pax8 target enrichment, IP and non-IP DNAs were modified for sequencing following the ChIP-Seq manufacturer’s protocol (Illumina, San Diego, CA). Briefly, DNAs were blunted with a combination of T4 DNA polymerase, Klenow polymerase, and T4 PNK. Then, a single 3′-end “A” base was added using Klenow exo (3′-to-5′ exo minus). Adapters provided by Illumina were ligated to the ends of the modified DNA before size selection of 200-bp fragments via polyacrylamide gel electrophoresis followed by extraction. The isolated DNA samples were used as the template for amplification by 18 cycles of PCR, and used for cluster generation on the Illumina Genome Analyzer II. Amplified products were column-purified with the QIAquick PCR Purification Kit (Qiagen, Dusseldorf, Germany) and assayed for quantity and quality with the Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA).
The 50 bp sequence reads were aligned to the rat genome (rn4; NCBI build 4) using the MIRO pipeline (Centre for Genomic Regulation, Barcelona, Spain) allowing 2 mismatches, and aligned tags were converted to BED format and used for identification of binding sites. In order to visualize the data in the University of California Santa Cruz genome browser (http://genome.ucsc.edu), the sequence reads were directionally extended to 300 bp, and for each base pair in the genome the number of overlapping sequence reads was determined and averaged over a 10 bp window. All sequencing data can be downloaded from Gene Expression Omnibus (GEO) under accession number GSE26938.
Peak finding and data analysis
MACS program (Model-based Analysis for ChIP-Seq, v.1.4.1) was used with default parameters to determine enriched Pax8 binding regions using non-IP DNA as control. PeakAnalyzer software  was assessed to identify functional elements proximal to the immunoprecipitated peaks using the annotation of the Ensembl release 63  on the RGSC genome assembly v3.4.
The distance correlation analysis was done with the GenomeInspector tool of the Genomatix suite (Genomatix Software GmbH, Munich, Germany), within +/−10Kb from the middle of the CpG islands and simple nucleotide repeats. CpG islands and simple nucleotide repeat coordinates were obtained from the UCSC Genome Browser .
We used two approaches for motif search in the Pax8 binding sites defined by MACS. The first method, based on MEME-ChIP  and focusing on the central 100 bp portion of each sequence, was used for the 500 peaks with the best FDR. The TOMTOM program was used to compare known transcriptional motifs with the motifs identified by MEME-ChIP . On the other hand, the RegionMiner tool of the Genomatix suite (Genomatix Software GmbH, Munich, Germany) identified the most over-represented motifs, based on the background of occurrences of the transcription factor binding sites (TFBSs) within the whole sequence of the rat genome (rn4; NCBI build 4); we compared these to the Pax8 binding regions identified by MACS. The motifs were ranked by the Z-score/fold change to obtain the most relevant sites .
Polyclonal antibodies (1 μg) were bound to Dynabeads (Invitrogen, Carlsbad, CA) and incubated with 200 μg of nuclear proteins extracted as described  from PCCl3 thyroid cells. The incubation was performed in 300ul of Immunoprecipitation (IP) buffer (20 mM HEPES pH 8, 10 mM KCl, 0,15 mM EGTA, 0,15 mM EDTA, 150 mM NaCl, 0.1% NP-40 with a cocktail of protease inhibitors (Roche, Manhein Germany)). After washing with IP buffer, proteins were eluted in 20 μl of Laemmli sample buffer and boiled for 10 minutes. The immunocomplexes were analyzed by SDS-PAGE and then immunoblotted using anti-Pax8 (BioPat, Milan Italy) and anti-Sp1 (Santa Cruz Biotechnology Inc., Santa Cruz, CA) antibodies. In the case of CTCF and prior to CoIP, we transfected into PCCl3 cells an expression vector containing the full-length human CTCF, as the antibody used (Upstate Biotechnology, Waltham, MASS) recognized the human form more specifically.
Promoter activity assays
HeLa cells were transiently transfected using calcium phosphate with 1 μg of pNIS-2.8 reporter alone or in combination with 0.5 μg of expression vectors for Pax8, Sp1 and CTCF as indicated in the text. The Renilla luciferase-encoding pRL-CMV vector (50 ng) was used to correct for transfection efficiency. Forty-eight hours after transfection cells were harvested, lysed, and analyzed for firefly and renilla luciferase activities by the Dual-Luciferase reporter assay system (Promega, Madison WI). Promoter activity was determined as the ratio between firefly and renilla luciferase and represented as relative luciferase activity. The results were expressed as the mean ± SD of three independent experiments, each performed in triplicate. Data were analyzed with GraphPad Prism (Intuitive Software for Science, San Diego, CA). Statistical significance was determined using an Anova one-way test, and differences were considered significant at a P < 0.05.
The Pax8-dependent gene expression study was performed in differentiated PCCl3 thyroid cells by means of expression arrays (Agilent rat whole genome 44 K arrays). For this purpose, we generated three different conditions to finally establish two main comparisons: wild type vs. Pax8-silenced PCCl3 cells (siPax8 PCCl3), and scrambled siRNA-treated (siScramble PCCl3) vs. Pax8-silenced PCCl3 cells. Given that each comparison was performed using quadruplicates and dye-swaps (Cy3 and Cy5 fluorochromes), our experimental design included sixteen independent competitive hybridizations (Additional file 14).
Transient transfections of PCCl3 cells were performed using Lipofectamine 2000 (Invitrogen, Carlsbad, CA), both for scrambled and for Pax8 siRNA conditions (10 ng siRNA /ml) (Dharmacon, Denver, USA). Pax8 silencing was tested by means of western blotting using a polyclonal Pax8 mouse antibody (Biopat, Milan, Italy) at different time points (24 and 48 hours) after transfection (Additional file 15). Once the 48 hours condition was defined as the best time point for Pax8 silencing, we performed additional transfections to isolate total RNA using TRIzol reagent (Invitrogen, Carlsbad, CA) for each condition considered (siPax8, scrambled siRNA and PCCl3 cells treated with lipofectamine) following the manufacturer’s recommended protocol. RNA quality was evaluated with the Agilent 2100 Bioanalyzer and later amplified and labelled by using the Low RNA Input Linear Amplification Kit PLUS, Two-Color (Agilent Technologies, Palo Alto, CA). Briefly, for each sample 2 μg of total input RNA were amplified in two rounds of amplification following the manufacturer’s instructions. First strand cDNA synthesis and amplification reactions were carried out using random and T7 primers, respectively. During the 2-hour in vitro transcription, Cy3- or Cy5-labeled CTP was incorporated into each amplified RNA (cRNA). Products of the reaction were then purified using RNAeasy mini spin columns (Qiagen, Dusseldorf, Germany). Hybridization and slide and image processing were carried out according to the manufacturer’s instructions (Two-Color Microarray-Based Gene Expression Analysis protocol). In each experiment, 825 ng of contrasting cRNA samples were fragmented at 60°C for 30 min and hybridized at 65°C for 17 hr. The slides were scanned at a 10 μm resolution using the Agilent G2565BA Microarray Scanner (Agilent Technologies, Palo Alto, CA). Signal quantification was carried out with Feature Extraction 9.1 software (Agilent Technologies, Palo Alto, CA), using default analysis parameters for Agilent’s whole rat genome 44 K gene expression arrays. Array data were normalized using loess and quantile methods for normalization within and between arrays, respectively. Differential expression analysis was done using Bioconductor’s limma package. At a later stage, we used the annotate package and the data base rgug4131a.db to obtain the annotations of the rat genome from Agilent. Genes that showed adjusted p-values <0.005 were considered differentially expressed both in wild type vs. siPax8 cells and in siScramble PCCl3 vs. siPax8 PCCl3 cells. Functional analysis of Gene Ontology (GO) terms was carried out using the FatiGO tool and gene set enrichment analysis was performed using FatiScan [76, 77]. All microarray data can be downloaded from the Gene Expression Omnibus (GEO) under accession number GSE26938.
Experimental validation for ChIP-seq and expression array data
Technical validations were performed by means of real-time PCR to verify the Pax8-dependency of 7 loci, which showed significant results for both expression profiling and ChIP-Seq, independently of the peak location along the considered gene. After performing an independent Pax8-chromatin IP, we obtained genomic DNA from IP, non-IP, and input samples, which were further used to amplify specific fragments contained in IP peaks (Additional file 13). The immunoprecipitation ratio for IP peaks was estimated comparing IP versus non-IP amplification values using as normalizing regions those mentioned above (Gad1 and Afm), which were confirmed as negative controls by the ChIP-Seq results.
siRNA transfection was done to obtain cDNA for each of the three conditions initially considered for expression profiling (Pax8 siRNA, scrambled siRNA, and wild type PCCl3 cells). Expression level changes were defined in 7 down-regulated genes and in 4 up-regulated genes by means of real-time PCR for fragments specifically amplifying transcripts of interest (Additional file 16), using GAPDH as a control for target gene expression normalization.
Public database accesion number
Sequencing and gene expression microarray data have been deposited in the GEO database (accession number GSE26938).
Chromatin immunoprecipitation followed by massive sequencing
Model-based Analysis for ChIP-Seq
We thank Dr. Vassart (Université Libre de Bruxelles, Brussels, Belgium), Dr. Angel Pascual (Instituto Investigaciones Biomédicas, Madrid Spain) and Dr. Colin Goding (Ludwig Institute for Cancer Research, UK) for kindly providing the Pax8, Sp1 and CTCF expression vectors, respectively, and Dr. Ronald Hartong for his criticisms and linguistic assistance.
We acknowledge the support of Grants BFU-2010-16025 from the Dirección General de Proyectos de Investigación; RD06/0020/0060 from FIS, Instituto de Salud Carlos III, and S2011/BMD-2328 TIRONET project from the Comunidad de Madrid (Spain).
S. Ruiz-Llorente holds a postdoctoral fellowship of the Instituto de Salud Carlos III (Contrato Postdoctoral de Perfeccionamiento; CD05-0055). A. Sastre-Perona holds a predoctoral fellowship from the Formación Personal Universitario (FPU) program. C. Montero-Conde holds a postdoctoral fellowship of the Spanish Ministry of Science and Innovation (MICINN; BMED2008-0659).
- Farnham PJ: Insights from genomic profiling of transcription factors. Nat Rev Genet. 2009, 10 (9): 605-616. 10.1038/nrg2636.PubMed CentralView ArticlePubMedGoogle Scholar
- Santisteban P, Bernal J: Thyroid development and effect on the nervous system. Rev Endocr Metab Disord. 2005, 6 (3): 217-228. 10.1007/s11154-005-3053-9.View ArticlePubMedGoogle Scholar
- De Felice M, Di Lauro R: Thyroid development and its disorders: genetics and molecular mechanisms. Endocr Rev. 2004, 25 (5): 722-746. 10.1210/er.2003-0028.View ArticlePubMedGoogle Scholar
- Plachov D, Chowdhury K, Walther C, Simon D, Guenet JL, Gruss P: Pax8, a murine paired box gene expressed in the developing excretory system and thyroid gland. Development. 1990, 110 (2): 643-651.PubMedGoogle Scholar
- Mansouri A, Chowdhury K, Gruss P: Follicular cells of the thyroid gland require Pax8 gene function. Nat Genet. 1998, 19 (1): 87-90. 10.1038/ng0598-87.View ArticlePubMedGoogle Scholar
- Macchia PE, Lapi P, Krude H, Pirro MT, Missero C, Chiovato L, Souabni A, Baserga M, Tassi V, Pinchera A: PAX8 mutations associated with congenital hypothyroidism caused by thyroid dysgenesis. Nat Genet. 1998, 19 (1): 83-86. 10.1038/ng0598-83.View ArticlePubMedGoogle Scholar
- Ohno M, Zannini M, Levy O, Carrasco N, di Lauro R: The paired-domain transcription factor Pax8 binds to the upstream enhancer of the rat sodium/iodide symporter gene and participates in both thyroid-specific and cyclic-AMP-dependent transcription. Mol Cell Biol. 1999, 19 (3): 2051-2060.PubMed CentralView ArticlePubMedGoogle Scholar
- Miccadei S, De Leo R, Zammarchi E, Natali PG, Civitareale D: The synergistic activity of thyroid transcription factor 1 and Pax 8 relies on the promoter/enhancer interplay. Mol Endocrinol. 2002, 16 (4): 837-846. 10.1210/me.16.4.837.View ArticlePubMedGoogle Scholar
- Zannini M, Francis-Lang H, Plachov D, Di Lauro R: Pax-8, a paired domain-containing protein, binds to a sequence overlapping the recognition site of a homeodomain and activates transcription from two thyroid-specific promoters. Mol Cell Biol. 1992, 12 (9): 4230-4241.PubMed CentralView ArticlePubMedGoogle Scholar
- Fraizer GC, Shimamura R, Zhang X, Saunders GF: PAX 8 regulates human WT1 transcription through a novel DNA binding site. J Biol Chem. 1997, 272 (49): 30678-30687. 10.1074/jbc.272.49.30678.View ArticlePubMedGoogle Scholar
- Deaton AM, Bird A: CpG islands and the regulation of transcription. Genes Dev. 2011, 25 (10): 1010-1022. 10.1101/gad.2037511.PubMed CentralView ArticlePubMedGoogle Scholar
- Esposito C, Miccadei S, Saiardi A, Civitareale D: PAX 8 activates the enhancer of the human thyroperoxidase gene. Biochem J. 1998, 331 (Pt 1): 37-40.PubMed CentralView ArticlePubMedGoogle Scholar
- Pellizzari L, Fabbro D, Lonigro R, Di Lauro R, Damante G: A network of specific minor-groove contacts is a common characteristic of paired-domain-DNA interactions. Biochem J. 1996, 315 (Pt 2): 363-367.PubMed CentralView ArticlePubMedGoogle Scholar
- Czerny T, Schaffner G, Busslinger M: DNA sequence recognition by Pax proteins: bipartite structure of the paired domain and its binding site. Genes Dev. 1993, 7 (10): 2048-2061. 10.1101/gad.7.10.2048.View ArticlePubMedGoogle Scholar
- Ho Sui SJ, Mortimer JR, Arenillas DJ, Brumm J, Walsh CJ, Kennedy BP, Wasserman WW: oPOSSUM: identification of over-represented transcription factor binding sites in co-expressed genes. Nucleic Acids Res. 2005, 33 (10): 3154-3164. 10.1093/nar/gki624.PubMed CentralView ArticlePubMedGoogle Scholar
- Planque N, Leconte L, Coquelle FM, Benkhelifa S, Martin P, Felder-Schmittbuhl MP, Saule S: Interaction of Maf transcription factors with Pax-6 results in synergistic activation of the glucagon promoter. J Biol Chem. 2001, 276 (38): 35751-35760. 10.1074/jbc.M104523200.View ArticlePubMedGoogle Scholar
- Garcia B, Santisteban P: PI3K is involved in the IGF-I inhibition of TSH-induced sodium/iodide symporter gene expression. Mol Endocrinol. 2002, 16 (2): 342-352. 10.1210/me.16.2.342.View ArticlePubMedGoogle Scholar
- Shih HA, Nathanson KL, Seal S, Collins N, Stratton MR, Rebbeck TR, Weber BL: BRCA1 and BRCA2 mutations in breast cancer families with multiple primary cancers. Clin Cancer Res. 2000, 6 (11): 4259-4264.PubMedGoogle Scholar
- Thedieck C, Kuczyk M, Klingel K, Steiert I, Muller CA, Klein G: Expression of Ksp-cadherin during kidney development and in renal cell carcinoma. Br J Cancer. 2005, 92 (11): 2010-2017. 10.1038/sj.bjc.6602597.PubMed CentralView ArticlePubMedGoogle Scholar
- Di Palma T, Conti A, de Cristofaro T, Scala S, Nitsch L, Zannini M: Identification of novel Pax8 targets in FRTL-5 thyroid cells by gene silencing and expression microarray analysis. PLoS One. 2011, 6 (9): e25162-10.1371/journal.pone.0025162.PubMed CentralView ArticlePubMedGoogle Scholar
- de Cristofaro T, Di Palma T, Fichera I, Lucci V, Parrillo L, De Felice M, Zannini M: An essential role for Pax8 in the transcriptional regulation of Cadherin-16 in thyroid cells. Mol Endocrinol. 2012, 26 (1): 67-78. 10.1210/me.2011-1090.View ArticlePubMedGoogle Scholar
- Zacchi P, Stenmark H, Parton RG, Orioli D, Lim F, Giner A, Mellman I, Zerial M, Murphy C: Rab17 regulates membrane trafficking through apical recycling endosomes in polarized epithelial cells. J Cell Biol. 1998, 140 (5): 1039-1053. 10.1083/jcb.140.5.1039.PubMed CentralView ArticlePubMedGoogle Scholar
- Muller T, Hess MW, Schiefermeier N, Pfaller K, Ebner HL, Heinz-Erian P, Ponstingl H, Partsch J, Rollinghoff B, Kohler H: MYO5B mutations cause microvillus inclusion disease and disrupt epithelial cell polarity. Nat Genet. 2008, 40 (10): 1163-1165. 10.1038/ng.225.View ArticlePubMedGoogle Scholar
- Min J, Zaslavsky A, Fedele G, McLaughlin SK, Reczek EE, De Raedt T, Guney I, Strochlic DE, Macconaill LE, Beroukhim R: An oncogene-tumor suppressor cascade drives metastatic prostate cancer by coordinately activating Ras and nuclear factor-kappaB. Nat Med. 2010, 16 (3): 286-294. 10.1038/nm.2100.PubMed CentralView ArticlePubMedGoogle Scholar
- Bianco AC, Kim BW: Deiodinases: implications of the local control of thyroid hormone action. J Clin Invest. 2006, 116 (10): 2571-2579. 10.1172/JCI29812.PubMed CentralView ArticlePubMedGoogle Scholar
- Nowak RB, Fischer RS, Zoltoski RK, Kuszak JR, Fowler VM: Tropomodulin1 is required for membrane skeleton organization and hexagonal geometry of fiber cells in the mouse lens. J Cell Biol. 2009, 186 (6): 915-928. 10.1083/jcb.200905065.PubMed CentralView ArticlePubMedGoogle Scholar
- Liu C, Papewalis C, Domberg J, Scherbaum WA, Schott M: Chemokines and autoimmune thyroid diseases. Hormone and metabolic research = Hormon- und Stoffwechselforschung = Hormones et metabolisme. 2008, 40 (6): 361-368. 10.1055/s-2008-1073153.View ArticlePubMedGoogle Scholar
- Min HS, Choe G, Kim SW, Park YJ, Park do J, Youn YK, Park SH, Cho BY, Park SY: S100A4 expression is associated with lymph node metastasis in papillary microcarcinoma of the thyroid. Modern pathology : an official journal of the United States and Canadian Academy of Pathology, Inc. 2008, 21 (6): 748-755.View ArticleGoogle Scholar
- Ying S, Dong S, Kawada A, Kojima T, Chavanas S, Mechin MC, Adoue V, Serre G, Simon M, Takahara H: Transcriptional regulation of peptidylarginine deiminase expression in human keratinocytes. J Dermatol Sci. 2009, 53 (1): 2-9. 10.1016/j.jdermsci.2008.09.009.View ArticlePubMedGoogle Scholar
- Riepe FG: Clinical and molecular features of type 1 pseudohypoaldosteronism. Horm Res. 2009, 72 (1): 1-9.View ArticlePubMedGoogle Scholar
- Montero-Conde C, Martin-Campos JM, Lerma E, Gimenez G, Martinez-Guitarte JL, Combalia N, Montaner D, Matias-Guiu X, Dopazo J, de Leiva A: Molecular profiling related to poor prognosis in thyroid carcinoma. Combining gene expression data and biological information. Oncogene. 2008, 27 (11): 1554-1561. 10.1038/sj.onc.1210792.View ArticlePubMedGoogle Scholar
- Stuart ET, Haffner R, Oren M, Gruss P: Loss of p53 function through PAX-mediated transcriptional repression. EMBO J. 1995, 14 (22): 5638-5645.PubMed CentralPubMedGoogle Scholar
- Dehbi M, Pelletier J: PAX8-mediated activation of the wt1 tumor suppressor gene. EMBO J. 1996, 15 (16): 4297-4306.PubMed CentralPubMedGoogle Scholar
- Scouten WT, Patel A, Terrell R, Burch HB, Bernet VJ, Tuttle RM, Francis GL: Cytoplasmic localization of the paired box gene, Pax-8, is found in pediatric thyroid cancer and may be associated with a greater risk of recurrence. Thyroid. 2004, 14 (12): 1037-1046. 10.1089/thy.2004.14.1037.View ArticlePubMedGoogle Scholar
- Scacheri PC, Davis S, Odom DT, Crawford GE, Perkins S, Halawi MJ, Agarwal SK, Marx SJ, Spiegel AM, Meltzer PS: Genome-wide analysis of menin binding provides insights into MEN1 tumorigenesis. PLoS Genet. 2006, 2 (4): e51-10.1371/journal.pgen.0020051.PubMed CentralView ArticlePubMedGoogle Scholar
- Yang A, Zhu Z, Kapranov P, McKeon F, Church GM, Gingeras TR, Struhl K: Relationships between p63 binding, DNA sequence, transcription activity, and biological function in human cells. Mol Cell. 2006, 24 (4): 593-602. 10.1016/j.molcel.2006.10.018.View ArticlePubMedGoogle Scholar
- A M, A AP, AS M, M B, Z M, M H: Evaluation of apoptosis-related genes; Fas (CD94), FasL (CD178) and TRAIL polymorphisms in Iranian multiple sclerosis patients. J Neurol Sci. 2011Google Scholar
- Illingworth RS, Bird AP: CpG islands–'a rough guide'. FEBS Lett. 2009, 583 (11): 1713-1720. 10.1016/j.febslet.2009.04.012.View ArticlePubMedGoogle Scholar
- Illingworth RS, Gruenewald-Schneider U, Webb S, Kerr AR, James KD, Turner DJ, Smith C, Harrison DJ, Andrews R, Bird AP: Orphan CpG islands identify numerous conserved promoters in the mammalian genome. PLoS Genet. 2010, 6 (9)::e1001134PubMed CentralView ArticlePubMedGoogle Scholar
- Kimura K, Wakamatsu A, Suzuki Y, Ota T, Nishikawa T, Yamashita R, Yamamoto J, Sekine M, Tsuritani K, Wakaguri H: Diversification of transcriptional modulation: large-scale identification and characterization of putative alternative promoters of human genes. Genome Res. 2006, 16 (1): 55-65.PubMed CentralView ArticlePubMedGoogle Scholar
- Cawley S, Bekiranov S, Ng HH, Kapranov P, Sekinger EA, Kampa D, Piccolboni A, Sementchenko V, Cheng J, Williams AJ: Unbiased mapping of transcription factor binding sites along human chromosomes 21 and 22 points to widespread regulation of noncoding RNAs. Cell. 2004, 116 (4): 499-509. 10.1016/S0092-8674(04)00127-8.View ArticlePubMedGoogle Scholar
- Guttman M, Amit I, Garber M, French C, Lin MF, Feldser D, Huarte M, Zuk O, Carey BW, Cassady JP: Chromatin signature reveals over a thousand highly conserved large non-coding RNAs in mammals. Nature. 2009, 458 (7235): 223-227. 10.1038/nature07672.PubMed CentralView ArticlePubMedGoogle Scholar
- Maunakea AK, Nagarajan RP, Bilenky M, Ballinger TJ, D'Souza C, Fouse SD, Johnson BE, Hong C, Nielsen C, Zhao Y: Conserved role of intragenic DNA methylation in regulating alternative promoters. Nature. 2010, 466 (7303): 253-257. 10.1038/nature09165.PubMed CentralView ArticlePubMedGoogle Scholar
- Ohler U, Liao GC, Niemann H, Rubin GM: Computational analysis of core promoters in the Drosophila genome. Genome Biol. 2002, 3 (12): RESEARCH0087Google Scholar
- Lim CY, Santoso B, Boulay T, Dong E, Ohler U, Kadonaga JT: The MTE, a new core promoter element for transcription by RNA polymerase II. Genes Dev. 2004, 18 (13): 1606-1617. 10.1101/gad.1193404.PubMed CentralView ArticlePubMedGoogle Scholar
- Juven-Gershon T, Hsu JY, Theisen JW, Kadonaga JT: The RNA polymerase II core promoter - the gateway to transcription. Curr Opin Cell Biol. 2008, 20 (3): 253-259. 10.1016/j.ceb.2008.03.003.PubMed CentralView ArticlePubMedGoogle Scholar
- Jin VX, Singer GA, Agosto-Perez FJ, Liyanarachchi S, Davuluri RV: Genome-wide analysis of core promoter elements from conserved human and mouse orthologous pairs. BMC Bioinforma. 2006, 7: 114-10.1186/1471-2105-7-114.View ArticleGoogle Scholar
- Di Palma T, D'Andrea B, Liguori GL, Liguoro A, de Cristofaro T, Del Prete D, Pappalardo A, Mascia A, Zannini M: TAZ is a coactivator for Pax8 and TTF-1, two transcription factors involved in thyroid differentiation. Exp Cell Res. 2009, 315 (2): 162-175. 10.1016/j.yexcr.2008.10.016.View ArticlePubMedGoogle Scholar
- Lanctot C, Cheutin T, Cremer M, Cavalli G, Cremer T: Dynamic genome architecture in the nuclear space: regulation of gene expression in three dimensions. Nat Rev Genet. 2007, 8 (2): 104-115. 10.1038/nrg2041.View ArticlePubMedGoogle Scholar
- Phillips JE, Corces VG: CTCF: master weaver of the genome. Cell. 2009, 137 (7): 1194-1211. 10.1016/j.cell.2009.06.001.PubMed CentralView ArticlePubMedGoogle Scholar
- Li L, Davie JR: The role of Sp1 and Sp3 in normal and cancer cell biology. Ann Anat. 2010, 192 (5): 275-283. 10.1016/j.aanat.2010.07.010.View ArticlePubMedGoogle Scholar
- De Leo R, Miccadei S, Zammarchi E, Civitareale D: Role for p300 in Pax 8 induction of thyroperoxidase gene expression. J Biol Chem. 2000, 275 (44): 34100-34105.View ArticlePubMedGoogle Scholar
- Grasberger H, Ringkananont U, Lefrancois P, Abramowicz M, Vassart G, Refetoff S: Thyroid transcription factor 1 rescues PAX8/p300 synergism impaired by a natural PAX8 paired domain mutation with dominant negative activity. Mol Endocrinol. 2005, 19 (7): 1779-1791. 10.1210/me.2004-0426.View ArticlePubMedGoogle Scholar
- Pasca di Magliano M, Di Lauro R, Zannini M: Pax8 has a key role in thyroid cell differentiation. Proc Natl Acad Sci U S A. 2000, 97 (24): 13144-13149. 10.1073/pnas.240336397.PubMed CentralView ArticlePubMedGoogle Scholar
- Robson EJ, He SJ, Eccles MR: A PANorama of PAX genes in cancer and development. Nat Rev Cancer. 2006, 6 (1): 52-62. 10.1038/nrc1778.View ArticlePubMedGoogle Scholar
- Arnaldi LA, Borra RC, Maciel RM, Cerutti JM: Gene expression profiles reveal that DCN, DIO1, and DIO2 are underexpressed in benign and malignant thyroid tumors. Thyroid. 2005, 15 (3): 210-221. 10.1089/thy.2005.15.210.View ArticlePubMedGoogle Scholar
- Rossi DL, Acebron A, Santisteban P: Function of the homeo and paired domain proteins TTF-1 and Pax-8 in thyroid cell proliferation. J Biol Chem. 1995, 270 (39): 23139-23142. 10.1074/jbc.270.39.23139.View ArticlePubMedGoogle Scholar
- Buttiglieri S, Deregibus MC, Bravo S, Cassoni P, Chiarle R, Bussolati B, Camussi G: Role of Pax2 in apoptosis resistance and proinvasive phenotype of Kaposi's sarcoma cells. J Biol Chem. 2004, 279 (6): 4136-4143.View ArticlePubMedGoogle Scholar
- Schebesta A, McManus S, Salvagiotto G, Delogu A, Busslinger GA, Busslinger M: Transcription factor Pax5 activates the chromatin of key genes involved in B cell signaling, adhesion, migration, and immune function. Immunity. 2007, 27 (1): 49-63. 10.1016/j.immuni.2007.05.019.View ArticlePubMedGoogle Scholar
- Albelda SM, Oliver PD, Romer LH, Buck CA: EndoCAM: a novel endothelial cell-cell adhesion molecule. J Cell Biol. 1990, 110 (4): 1227-1237. 10.1083/jcb.110.4.1227.View ArticlePubMedGoogle Scholar
- Cali G, Gentile F, Mogavero S, Pallante P, Nitsch R, Ciancia G, Ferraro A, Fusco A, Nitsch L: CDH16/Ksp-Cadherin Is Expressed in the Developing Thyroid Gland and Is Strongly Down-Regulated in Thyroid Carcinomas. Endocrinology. 2011Google Scholar
- de Cristofaro T, Di Palma T, Fichera I, Lucci V, Parrillo L, De Felice M, Zannini M: An Essential Role for Pax8 in the Transcriptional Regulation of Cadherin-16 in Thyroid Cells. Mol Endocrinol. 2011Google Scholar
- Fusco A, Berlingieri MT, Di Fiore PP, Portella G, Grieco M, Vecchio G: One- and two-step transformations of rat thyroid epithelial cells by retroviral oncogenes. Mol Cell Biol. 1987, 7 (9): 3365-3370.PubMed CentralView ArticlePubMedGoogle Scholar
- Costamagna E, Garcia B, Santisteban P: The functional interaction between the paired domain transcription factor Pax8 and Smad3 is involved in transforming growth factor-beta repression of the sodium/iodide symporter gene. J Biol Chem. 2004, 279 (5): 3439-3446.View ArticlePubMedGoogle Scholar
- Landa I, Ruiz-Llorente S, Montero-Conde C, Inglada-Perez L, Schiavi F, Leskela S, Pita G, Milne R, Maravall J, Ramos I: The variant rs1867277 in FOXE1 gene confers thyroid cancer susceptibility through the recruitment of USF1/USF2 transcription factors. PLoS genetics. 2009, 5 (9): e1000637-10.1371/journal.pgen.1000637.PubMed CentralView ArticlePubMedGoogle Scholar
- Vilain C, Rydlewski C, Duprez L, Heinrichs C, Abramowicz M, Malvaux P, Renneboog B, Parma J, Costagliola S, Vassart G: Autosomal dominant transmission of congenital thyroid hypoplasia due to loss-of-function mutation of PAX8. J Clin Endocrinol Metab. 2001, 86 (1): 234-238. 10.1210/jc.86.1.234.PubMedGoogle Scholar
- Courey AJ, Holtzman DA, Jackson SP, Tjian R: Synergistic activation by the glutamine-rich domains of human transcription factor Sp1. Cell. 1989, 59 (5): 827-836. 10.1016/0092-8674(89)90606-5.View ArticlePubMedGoogle Scholar
- Klenova EM, Nicolas RH, Paterson HF, Carne AF, Heath CM, Goodwin GH, Neiman PE, Lobanenkov VV: CTCF, a conserved nuclear factor required for optimal transcriptional activity of the chicken c-myc gene, is an 11-Zn-finger protein differentially expressed in multiple forms. Mol Cell Biol. 1993, 13 (12): 7612-7624.PubMed CentralView ArticlePubMedGoogle Scholar
- Lee TI, Johnstone SE, Young RA: Chromatin immunoprecipitation and microarray-based analysis of protein location. Nat Protoc. 2006, 1 (2): 729-748. 10.1038/nprot.2006.98.PubMed CentralView ArticlePubMedGoogle Scholar
- Salmon-Divon M, Dvinge H, Tammoja K, Bertone P: PeakAnalyzer: genome-wide annotation of chromatin binding and modification loci. BMC Bioinforma. 2010, 11: 415-10.1186/1471-2105-11-415.View ArticleGoogle Scholar
- Flicek P, Amode MR, Barrell D, Beal K, Brent S, Chen Y, Clapham P, Coates G, Fairley S, Fitzgerald S: Ensembl 2011. Nucleic Acids Res. 2011, 39 (Database issue): D800-806.PubMed CentralView ArticlePubMedGoogle Scholar
- Fujita PA, Rhead B, Zweig AS, Hinrichs AS, Karolchik D, Cline MS, Goldman M, Barber GP, Clawson H, Coelho A: The UCSC Genome Browser database: update 2011. Nucleic Acids Res. 2011, 39 (Database issue): D876-882.PubMed CentralView ArticlePubMedGoogle Scholar
- Machanick P, Bailey TL: MEME-ChIP: motif analysis of large DNA datasets. Bioinformatics. 2011, 27 (12): 1696-1697. 10.1093/bioinformatics/btr189.PubMed CentralView ArticlePubMedGoogle Scholar
- Gupta S, Stamatoyannopoulos JA, Bailey TL, Noble WS: Quantifying similarity between motifs. Genome Biol. 2007, 8 (2): R24-10.1186/gb-2007-8-2-r24.PubMed CentralView ArticlePubMedGoogle Scholar
- Perona R, Montaner S, Saniger L, Sanchez-Perez I, Bravo R, Lacal JC: Activation of the nuclear factor-kappaB by Rho, CDC42, and Rac-1 proteins. Genes Dev. 1997, 11 (4): 463-475. 10.1101/gad.11.4.463.View ArticlePubMedGoogle Scholar
- Al-Shahrour F, Minguez P, Tarraga J, Medina I, Alloza E, Montaner D, Dopazo J: FatiGO +: a functional profiling tool for genomic data. Integration of functional annotation, regulatory motifs and interaction data with microarray experiments. Nucleic Acids Res. 2007, 35 (Web Server issue): W91-96.PubMed CentralView ArticlePubMedGoogle Scholar
- Medina I, Carbonell J, Pulido L, Madeira SC, Goetz S, Conesa A, Tarraga J, Pascual-Montano A, Nogales-Cadenas R, Santoyo J: Babelomics: an integrative platform for the analysis of transcriptomics, proteomics and genomic data with advanced functional profiling. Nucleic Acids Res. 2010, 38 (Web Server issue): W210-213.PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.