Mouse Homologue of the Schizophrenia Susceptibility Gene ZNF804A as a Target of Hoxc8
Using a ChIP-cloning technique, we identified a Zinc finger protein 804a (Zfp804a) as one of the putative Hoxc8 downstream target genes. We confirmed binding of Hoxc8 to an intronic region of Zfp804a by ChIP-PCR in F9 cells as well as in mouse embryos. Hoxc8 upregulated Zfp804a mRNA levels and augmented minimal promoter activity in vitro. In E11.5 mouse embryos, Zfp804a and Hoxc8 were coexpressed. Recent genome-wide studies identified Zfp804a (or ZNF804A in humans) as a plausible marker for schizophrenia, leading us to hypothesize that this embryogenic regulatory control might also exert influence in development of complex traits such as psychosis.
Hox genes encode transcription factors containing a 61-amino acid motif that enables binding and transcriptional regulation of target gene expression [1–3]. The mammalian Hox genes have 39 members organized in four clusters (HoxA-D) each located on a different chromosome. Resulting from successive duplications during the evolution of an ancestral Hox complex, each cluster contains a series of paralogues (groups I–XIII). Each Hox gene has a conserved helix-turn-helix DNA-binding domain that recognizes and binds to specific sequence motifs (TAAT/ATTA/TTAT/ATAA) located in promoter sequences, intergenic and intronic regions. Such conserved sequence motif is also found repeatedly several times in the genome underscoring a potentially expansive regulatory landscape of Hox genes [4, 5]. Hox expression domains along the anterior-posterior (A-P) axis is characterized as overlapping, with the order and timing of expression correlating to the position of genes along a cluster. The combinatorial expression of various Hox genes constitutes the proverbial “Hox code”: the genetic program that determines cell fates, that is, morphogenesis, growth, and differentiation during embryonic development . Meanwhile, in adult cells, Hox genes assume some poorly understood oncogenic roles. For example, expression of HOXC5 and HOXC8 is selectively elevated in human cervical cancer cells when compared to normal keratinocytes . In human prostate adenocarcinomas, HOXC8 expression is directly correlated with the loss of tumor differentiation based on Gleason scoring . Moreover, Hox genes have been ascribed with “nonconventional” roles such as postnatal neuronal circuit development  although data is still scant.
To gain better understanding of the mechanisms by which Hox functions during mouse embryogenesis, the identification and functional analysis of target genes are necessary. In this study, we focused our analysis on the downstream targets of Hoxc8, one of three members of the paralogous group VIII, which our group and other laboratories have extensively studied in the past years [6, 10–13]. Several molecular tools have allowed us to identify the developmental regulatory targets of Hoxc8 in a high-throughput manner. Microarray expression analysis by the group of Frank Ruddle has identified 34 genes involved in cell adhesion, migration, metabolism, and apoptosis whose expression levels were differentially changed in mouse embryonic cells overexpressing Hoxc8 . Likewise, our group identified, using proteomics, 15 genes associated with cell motility, protein homeostasis, and metabolism . More recently, modified chromatin immunoprecipitation- (ChIP-) based methods have been developed, such as ChIP-cloning, ChIP Display, Differential Chromatin Scanning, and ChIP-chip and DamID chromatin profiling, permitting the screening for DNA-protein interaction in vivo at a genome-wide scale . We took advantage of a ChIP-cloning gene discovery strategy  to search for novel downstream target genes of Hoxc8 that may be important during mouse development. We identified 12 putative target sites by ChIP-cloning and subjected these to comparative genomics analysis to identify Hoxc8-binding motifs that are conserved across the mammalian lineage leading to the identification of Zfp804a, as a strong candidate target. We further confirm the regulation of Zfp804a by ChIP-PCR, coexpression analysis and a reporter assay.
2. Materials and Methods
2.1. Animal Preparation
In order to get E11.5 embryos, ICR mice were allowed to mate at 6:00 pm, and presence of a vaginal plug in the following morning indicates the 0.5 day post coitum (dpc) time point. After 11 days, the pregnant mice were sacrificed and E11.5 embryos were taken. Maternal and extraembryonic tissues were removed.
2.2. Cell Culture and Transfection
F9 cells were cultured in DMEM media supplemented with 10% of fetal bovine serum (FBS) and 100 μg/ml penicillin-streptomycin (Invitrogen, Carlsbad, CA, USA) in plates that were coated with 0.1% gelatin. The cells were incubated at in a 5% CO2 environment. In 100 mm dishes, upon reaching ~50%, the cells were transfected with 4 μg of pcDNA3 (empty vector control) or pcDNA3:Hoxc8 plasmid DNA using Lipofectamine 2000 system (Invitrogen).
2.3. ChIP Assays
ChIP assays were performed on fresh E11.5 embryos with the head and internal organs removed using a ChIP assay kit (Upstate Biotechnolgy, Lake Placid, NY, USA) with slight modification. To enrich specificity, immunoprecipitation step was performed thrice using a Hoxc8-monoclonal antibody (MMS-266R, Covance Research products, Princeton, NJ, USA). Immunoprecipitation in the absence of the antibody served as negative control. ChIP-DNA was cloned into pBluescriptII SK(+) vector and its inserts were sequenced . We reported our identification of Hoxc8 gene targets using this ChIP-cloning protocol in our previous paper .
In order to verify whether the identified genes are a target of Hoxc8, we PCR-amplified fragments from ChIP-DNA of E11.5 embryos and F9 cells using the primers F, ggagtcacaaccaaatctgc and R, cagggaattttgggcattat. PCR was carried out with a 5-minute hot start at , followed by 33 cycles of denaturation ( for 1 minute), annealing (1 minute at ), and extension (1 minute at ).
2.4. RNA Isolation and RT-PCR
Total cellular RNA was extracted with RNAzol Bee (Tel-test, Inc. Industries Inc., Texas, USA), partitioned with chloroform/isoamyl alcohol and centrifuged at 12 000 rpm at for 15 minutes. RNA was precipitated by adding equal volume of cold iso-propanol and centrifugation at 12 000 rpm at for 15 minutes. Pellet was washed with 70% ethanol, air-dried, and resuspended in DEPC-H2O. RT-PCR was performed on 2 μg of the RNA with the following PCR conditions: 5 minutes hot start at , followed by denaturation at for 1 minute, annealing step at the specific temperatures for 50 seconds (-actin, , Hoxc8, , and Zfp804a, ), and extension step at for 1 minute for 28 cycles. Each reaction was terminated after a final elongation step of for 7 minutes. The sequences of the primers were as follows: β-actin (F, catgtttgagaccttcaacacccc; R, gccatctcctgctcgaagtctag), Hoxc8 (F, cctattacgactgccggttc; R, ttggcggaggatttacagtc), and Zfp804a (F, tcacaagggacccatcaaat; R, ttagcatgttcggctgaatg).
2.5. Effector-Reporter Assay
An effector-reporter assay was done using the Dual-Luciferase Reporter Assay System (Promega, Madison, WI, USA). Zfp804a gene fragment was cloned into the pGL3MP vector  to make the effector-reporter construct. For the luciferase assay, F9 cells/well were seeded into a 6-well plate and cultured until 50% of confluency. Four micrograms of effector (pcDNA3:Hoxc8; pcDNA3:Hoxc8 +siHoxc8) were cotransfected in the cells with 1 μg of the reporter (pGL3-MP; pGL2-MP:Zfp804a) and 25 ng of pRL-TK vectors. After 36 hours, the cells were lysed and luciferase activity was measured using the GLOMAX 20/20 luminometer (Promega).
3. Results and Discussion
Hox genes play a role as a “realisator” of body pattern formation [2, 3]. Since the discovery of Drosophila Hox genes in the early 1980s, there has been an avalanche of studies on their roles in mammalian development. However, while much is already known about Hox gene structure and its function, how it determines the precise formation of the body parts along the A-P axis is still enigmatic. In order to gain further insight into the molecular basis of how various Hox genes control tissue fate, it is necessary that their downstream genes are identified and characterized.
We obtained target genes of Hoxc8 by ChIP-cloning  and further demonstrated regulation of a candidate gene. ChIP-cloning approach is a modification of the ChIP assay which combines immunoprecipitation of sheared chromatin with subsequent DNA cloning and sequencing (Figure 1). In our ChIP-cloning procedures, we used E11.5 mouse embryos because Hoxc8 expression is relatively elevated at this particular stage compared to other time points . After removal of the head and internal organs, the embryo was fixed with formaldehyde to cross-link DNA-protein interactions. After the DNA complexes were pulled down by an anti-Hoxc8 monoclonal antibody, purified ChIPed DNA fragments were cloned and sequenced. Out of 146 clones, we identified 12 candidate genes from 16 sequences, most of which contained Hox-binding motifs (see Table 1). As it is logical to assume that Hox-binding will likely occur in the promoter regions of its target genes, in our experiment, all our cloned DNA fragments were derived from intronic and intergenic regions. A similar observation was reported by Birkaya et al.  showing that 40% of DeltaNp63-ChIP fragments were from the introns. One of the cloned regions of our candidate targets, calpain 10 (NM_011796.2), does not have a Hox-binding motif implying that if this was confirmed as true target, Hoxc8 may likely regulate the gene via indirect binding, probably in collaboration with a cofactor(s), as in the case of TGF- signaling with Hoxc8 [20–22]. Among the candidate genes, we chose to study zinc finger protein 804a (Zfp804a; NM_175513) since this cloned intronic region contains ~8.4 Hox-binding sites per 100 bp. In subsequent in silico analysis, however, out of the 75 Hox-binding motifs present, only 4 were conserved across a wide range of mammalian species: human, chimpanzee, rhesus, orangutan, marmoset, cat, horse, dog, opossum, mouse, rat, and guinea pig.
To further demonstrate Zfp804a-Hoxc8-binding, we validated the interaction by performing ChIP-PCR on cell extracts from Hoxc8-overexpressing F9 murine embryonic teratocarcinoma cell line (in vitro) and E11.5 mouse embryos (in vivo). Located in chromosome 2, the Zfp804a gene contains 4 coding exons with a noncontiguous string of Hoxc8-binding motifs within its first intron. We designed primers to target the area where putative Hox-binding sites are abundant (Figure 2(a)). Zfp804a gene was clearly amplified from the DNA samples that had been precipitated by a commercial monoclonal (Figure 2(b)) and an in-house polyclonal antibody (data not shown) but not by an antimouse IgG. The strongest band of Zfp804a among 12 candidate genes in ChIP-PCR led us to focus our analysis on it (data not shown). The Zfp804a band was also detected from input lane which served as the internal control. Our results, therefore, indicate that the Hoxc8 protein stably and directly binds to the first intron of Zfp804a.
Does this binding lead to a modulation of gene expression? Given the fact that the regulation of Zfp804a is largely unknown, we utilized an unbiased effector-reporter assay by cloning the Zfp804a intronic fragment into the minimal promoter plasmid, pGL3-MP (pGL3-MP:Zfp804a). With transfection into F9 cells and assaying for dual-luciferase activity after 36 hours, enhancement of promoter activity relative to pGL3-MP was seen when Hoxc8 expression plasmid was cotransfected compared with the empty vector (pcDNA3 alone). In addition, we were able to demonstrate that the augmented promoter activity was abolished by adding Hoxc8 siRNA (Figure 2(c)). These data suggest that the intronic fragment of Zfp804a may function as an alternative promoter, and the modulation of activity of the minimal promoter  by such intronic sequence endorses the functionality of the Hoxc8-Zfp804a interaction.
Finally, we studied whether Hoxc8 modulates expression of Zfp804a. For this, we transiently transfected with a Hoxc8 overexpression construct F9 cells having undetectable levels of Hoxc8 and Zfp804a mRNAs. Compared with the empty vector control, the cells overexpressing Hoxc8 showed induction of Zfp804a (Figure 3(a)) indicating that, in this particular in vitro system, expression of Zfp804a was positively regulated by Hoxc8. To prove in vivo, we analyzed levels of Hoxc8 and Zfp804a mRNA transcripts along the A-P axis of E11.5 embryos. From Figure 3(b), Hoxc8 was expressed in the trunk area (regions 3 to 5) except in head whereas Zfp804a mRNA was seen generally through-out the trunk with only a residual signal in the head area (Figure 3(b)) . Thus, the expression patterns on Zfp804a correlate to some extent with that of Hoxc8 in this developmental stage. Since the regulation target genes by Hox is complex and multifactorial [24, 25], we are cautious in interpreting such data since Hox regulates its target in a context-dependent manner and proves to be both a repressor and activator at the same time .
An interesting twist in the relationship of Hoxc8-Zfp804a regulation lies in the recent discovery that its human homologue, ZNF804A, turns out to be a marker for schizophrenia and bipolar disorders . Although ZNF804A is uncharacterized and its functions are still unknown to date, it is interesting that fMRI images of ZNF804A SNPs carriers show altered functional coupling between the dorsolateral prefrontal cortex (DLPFC) and the hippocampus . Since genetic variation in dopaminergic and glutaminergic neurotransmission affects DLPFC or hippocampus connectivity, the examination of ZNF804A and HOXC8 in these neurotransmitter cascades would be interesting in the future. During embryogenesis, Hox gene expression is activated prior to rhombomeric segmentation, and its anterior boundaries, limited to those in the paralogue groups I–IV, are found only in the neuroepithelia that give rise to the hindbrain . Because of this bias, no systematic study on neurobiology of the other Hox paralogues has been done. Therefore, the finding here is rather surprising and very intriguing.
Recently, some HOX genes have been reported to be expressed in adult cells, and misexpression of certain HOX gene(s) can lead to diseases such as cancer . Therefore we examined the expression of Hoxc8 in adult mouse brain. As shown in Figure 3(c), we were also able to verify the expression of Hoxc8 in adult normal brain tissue, cortex, and hippocampus by RT-PCR. Interestingly, Zfp804a was only detected in cortex, but not in hippocampus. Consistent with this expression patterns, Hoxc8 showed a prominent binding to Zfp804a in cortex as well as whole brain (Figure 3(d)).
Given that Hox genes are master transcription factors at the apex of the genetic hierarchy of developmental control, it is tantalizing to regard that their “fetal” molecular networks would extend influence to adulthood and more so, in such complex systems, such as psychosis. A future challenge from our work is to elucidate the molecular mechanisms responsible for this phenomenon and the potential implications of Hox regulatory control in human behavior.
H.J. Chung, H. Min and S.H. Kim are graduate students supported by BK21 (Brain Korea 21) scholarship. The authors thank Kalyani Ruthala for her technical assistance; Drs. Ethan Patterson (Johns Hopkins University), Jinwoong Bok, and Heejei Yoon for the helpful discussions, and JiHyun Ma and Dong-Su Jang for their help in manuscript preparation. This work was supported by Grants 2008-0058561 (2010-0000155) and 2009-0070964 from NRF and partly by 20070401-034-030 from the BioGreen21 Program, RDA, Korea. H. J. Chung and J.-Y. Lee contributed equally to this work.
R. Krumlauf, “Evolution of the vertebrate Hox homeobox genes,” Bioessays, vol. 14, no. 4, pp. 245–252, 1992.View at: Google Scholar
R. Krumlauf, “Hox genes in vertebrate development,” Cell, vol. 78, no. 2, pp. 191–201, 1994.View at: Publisher Site | Google Scholar
W. McGinnis and R. Krumlauf, “Homeobox genes and axial patterning,” Cell, vol. 68, no. 2, pp. 283–302, 1992.View at: Google Scholar
A. Awgulewitsch, C. Bieberich, L. Bogarad, C. Shashikant, and F. H. Ruddle, “Structural analysis of the Hox-3.1 transcription unit and the Hox—3.2-Hox-3.1 intergenic region,” Proceedings of the National Academy of Sciences of the United States of America, vol. 87, no. 16, pp. 6428–6432, 1990.View at: Publisher Site | Google Scholar
H. Hayashi, S. Abdollah, Y. Qiu et al., “The MAD-related protein Smad7 associates with the TGFbeta receptor and functions as an antagonist of TGFbeta signaling,” Cell, vol. 89, no. 7, pp. 1165–1173, 1997.View at: Google Scholar
H. Lei, H. Wang, A. H. Juan, and F. H. Ruddle, “The identification of Hoxc8 target genes,” Proceedings of the National Academy of Sciences of the United States of America, vol. 102, no. 7, pp. 2420–2424, 2005.View at: Publisher Site | Google Scholar
Y. Alami, V. Castronovo, D. Belotti, D. Flagiello, and N. Clausse, “HOXC5 and HOXC8 expression are selectively turned on in human cervical cancer cells compared to normal keratinocytes,” Biochemical and Biophysical Research Communications, vol. 257, no. 3, pp. 738–745, 1999.View at: Publisher Site | Google Scholar
D. Waltregny, Y. Alami, N. Clausse, J. De Leval, and V. Castronovo, “Overexpression of the homeobox gene HOXC8 in human prostate cancer correlates with loss of tumor differentiation,” Prostate, vol. 50, no. 3, pp. 162–169, 2002.View at: Publisher Site | Google Scholar
Y. Narita and F. M. Rijli, “Hox genes in neural patterning and circuit formation in the mouse hindbrain,” Current Topics in Developmental Biology, vol. 88, pp. 139–167, 2009.View at: Publisher Site | Google Scholar
H. J. Chung, M. Kang, and M. H. Kim, “Isolation of putative in vivo hoxc8 downstream target genes using ChIP-cloning method,” Journal of Biomedical Science, vol. 88, pp. 47–53, 2009.View at: Google Scholar
H. Lei, A. H. Juan, M.-S. Kim, and F. H. Ruddle, “Identification of a Hoxc8-regulated transcriptional network in mouse embryo fibroblast cells,” Proceedings of the National Academy of Sciences of the United States of America, vol. 103, no. 27, pp. 10305–10309, 2006.View at: Publisher Site | Google Scholar
H. Lei, A. H. Juan, M.-S. Kim, and F. H. Ruddle, “Mouse naked cuticle 2 (mNkd2) as a direct transcriptional target of Hoxc8 in vivo,” Journal of Experimental Zoology Part A, vol. 307, no. 1, pp. 1–6, 2007.View at: Publisher Site | Google Scholar
H. Min, M. Kang, and M. H. Kim, “A homeotic gene, Hoxc8 regulates the expression of proliferating cell nuclear antgen in NIH3T3 cell,” Journal of Biomedical Science, vol. 13, pp. 239–244, 2007.View at: Google Scholar
Y. Kwon, J. H. Ko, B. Kim, and M. H. Kim, “Analysis of plausible downstream target genes of Hoxc8 in F9 teratocarcinoma cells: putative downstream target genes of Hoxc8,” Molecular Biology Reports, vol. 30, no. 3, pp. 141–148, 2003.View at: Publisher Site | Google Scholar
B. A. T. Rodriguez and T. H. Huang, “Tilling the chromatin landscape: emerging methods for the discovery and profiling of protein-DNA interactions,” Biochemistry and Cell Biology, vol. 83, no. 4, pp. 525–534, 2005.View at: Publisher Site | Google Scholar
M. J. Solomon, P. L. Larsen, and A. Varshavsky, “Mapping protein-DNA interactions in vivo with formaldehyde: evidence that histone H4 is retained on a highly transcribed gene,” Cell, vol. 53, no. 6, pp. 937–947, 1988.View at: Google Scholar
G. Torres, P. D. Frisella, S. J. Yousuf et al., “A ChIP-cloning approach linking SIRT1 to transcriptional modificationof DNA targets,” Biotechniques, vol. 44, no. 7, pp. Pxii–Pxiv, 2008.View at: Publisher Site | Google Scholar
E. S. Patterson, R. C. Addis, M. J. Shamblott, and J. D. Gearhart, “SOX17 directly activates Zfp202 transcription during in vitro endoderm differentiation,” Physiological Genomics, vol. 34, no. 3, pp. 277–284, 2008.View at: Publisher Site | Google Scholar
B. Birkaya, K. Ortt, and S. Sinha, “Novel in vivo targets of Np63 in keratinocytes identified by a modified chromatin immunoprecipitation approach,” BMC Molecular Biology, vol. 8, pp. 1–15, 2007.View at: Publisher Site | Google Scholar
X. Yang, X. Ji, X. Shi, and X. Cao, “Smad1 domains interacting with Hoxc-8 induce osteoblast differentiation,” The Journal of Biological Chemistry, vol. 275, no. 2, pp. 1065–1072, 2000.View at: Publisher Site | Google Scholar
T. G. Hullinger, Q. Pan, H. L. Viswanathan, and M. J. Somerman, “TGFbeta and BMP-2 activation of the OPN promoter: roles of smad- and hox-binding elements,” Experimental Cell Research, vol. 262, no. 1, pp. 69–74, 2001.View at: Publisher Site | Google Scholar
X. Li, S. Nie, C. Chang, T. Qiu, and X. Cao, “Smads oppose Hox transcriptional activities,” Experimental Cell Research, vol. 312, no. 6, pp. 854–864, 2006.View at: Publisher Site | Google Scholar
Y. Kwon, J. Shin, H. W. Park, and M. H. Kim, “Dynamic expression pattern of Hoxc8 during mouse early embryogenesis,” Anatomical Record A, vol. 283, no. 1, pp. 187–192, 2005.View at: Publisher Site | Google Scholar
D. F. B. Miller, B. T. Rogers, A. Kalkbrenner, B. Hamilton, S. L. Holtzman, and T. Kaufman, “Cross-regulation of Hox genes in the Drosophila melanogaster embryo,” Mechanisms of Development, vol. 102, no. 1-2, pp. 3–16, 2001.View at: Publisher Site | Google Scholar
V. S. Chopra and R. K. Mishra, ““Mir”acles in hox gene regulation,” Bioessays, vol. 28, no. 5, pp. 445–448, 2006.View at: Publisher Site | Google Scholar
M. Porsch-Ozcurumez, T. Langmann, S. Heimerl et al., “The zinc finger protein 202 (ZNF202) is a transcriptional repressor of ATP binding cassette transporter A1 (ABCA1) and ABCG1 gene expression and a modulator of cellular lipid efflux,” The Journal of Biological Chemistry, vol. 276, no. 15, pp. 12427–12433, 2001.View at: Publisher Site | Google Scholar
M. C. O'Donovan, N. Craddock, N. Norton et al., “Identification of loci associated with schizophrenia by genome-wide association and follow-up,” Nature Genetics, vol. 40, no. 9, pp. 1053–1055, 2008.View at: Publisher Site | Google Scholar
C. Esslinger, H. Walter, P. Kirsch et al., “Neural mechanisms of a genome-wide supported psychosis variant,” Science, vol. 324, no. 5927, p. 605, 2009.View at: Publisher Site | Google Scholar
M. Davenne, M. K. Maconochie, R. Neun et al., “Hoxa2 and Hoxb2 control dorsoventral patterns of neuronal development in the rostral hindbrain,” Neuron, vol. 22, no. 4, pp. 677–691, 1999.View at: Publisher Site | Google Scholar
K. C. Wang, J. A. Helms, and H. Y. Chang, “Regeneration, repair and remembering identity: the three Rs of Hox gene expression,” Trends in Cell Biology, vol. 19, no. 6, pp. 268–275, 2009.View at: Publisher Site | Google Scholar