Andrea M. Santangelo,1,
Viviana Bumaschny,1
María Elena Avale,1
James L. Smart,2
Malcolm J. Low,2,3,4 and
Marcelo Rubinstein1,4,5*
Instituto de Investigaciones en Ingeniería Genética y Biología Molecular, CONICET and Departamento Fisiología, Biología Molecular y Celular, FCEyN, Universidad de Buenos Aires, Buenos Aires, Argentina,1 Vollum Institute,2 Department of Behavioral Neuroscience,3 Center for the Study of Weight Regulation and Associated Disorders, Oregon Health & Science University, Portland, Oregon,4 Centrode Estudios Científicos, Valdivia, Chile5
Received 25 October 2004/ Returned for modification 7 January 2005/ Accepted 18 January 2005
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
-melanocyte stimulating hormone (
-MSH) and ß-endorphin. In the brain, a specialized population of neurons concentrated in the hypothalamic arcuate nucleus also express Pomc, giving rise to
-, ß-, and
-melanocortins and the potent analgesic opioid peptide ß-endorphin (46, 50). During the last few years, POMC hypothalamic neurons have received a great deal of attention because they express receptors for the adipostatic hormone leptin and play a central role in the control of energy homeostasis and body weight regulation (10, 52, 56). Once leptin crosses the blood-brain barrier, it stimulates both the transcription of the POMC gene and the activity of POMC arcuate neurons, leading to the release of melanocortins which, in turn, stimulate central melanocortin receptors to decrease food intake and increase metabolic rate (4, 52). The importance of the central melanocortin pathway in feeding behavior is clearly observed in mice and humans with null Pomc homozygous mutations, which display hyperphagia and early-onset obesity (8, 29, 30, 65). Even though total POMC deficiency is very rare in humans, POMC is a strong candidate gene to predispose to familial obesity. Several independent genome-wide scans for quantitative trait loci (QTL) have found a highly significant genetic linkage between a relatively narrow region in chromosome 2 containing the POMC locus and obesity-related traits (11, 15, 47). However, polymorphisms in the coding sequences of POMC that alter the structure or function of POMC peptides apparently do not account for this correlation (25), suggesting the alternative possibility that mutations in noncoding regulatory sequences may alter the level of POMC RNA transcripts and consequently the concentration of POMC peptides in the brain.
Transcriptional regulation of Pomc has been studied primarily in pituitary corticotrophs and melanotrophs because it has been possible to combine in vitro transfection in immortalized pituitary cell lines (40) with in vivo expression in transgenic mice (33, 34). Cell-specific expression of the rat or mouse Pomc genes in the pituitary is controlled by the combinatorial presence of cis-acting elements localized within approximately 400 nucleotides adjacent to the transcriptional start site (33, 34). Several transcription factors including Nurr77, NeuroD1/BETA2 heterodimers, Pitx1, SP1, and Tpit have been implicated in the pituitary regulation of Pomc expression or the determination of corticotroph and melanotroph lineages (31, 34, 43, 44, 45, 60). In contrast, little is known about enhancer regions and transcription factors that regulate Pomc expression in arcuate neurons. In previous studies, we demonstrated that the 400-bp proximal Pomc promoter region is incapable of directing reporter gene expression to the POMC arcuate neurons of transgenic mice (49). However, the inclusion of additional genomic sequences suggested that cis-acting regions controlling Pomc expression in the brain are localized within 11 kb of distal 5' flanking sequences upstream of the Pomc promoter (13, 41, 66). Combining a functional expression analysis in transgenic mice together with in silico phylogenetic footprinting, we report here the identification of two conserved enhancers that are necessary and sufficient for transgene expression in POMC arcuate neurons. Our results also indicate that the Pomc promoter has a modular architecture that can allow independent transcriptional control of Pomc expression in the brain and the pituitary by a distal and proximal set of enhancers, respectively. These results advance the understanding of the molecular identity of POMC neurons and might contribute to a further determination of the role that POMC plays in feeding centers of the brain and the genetics of obesity.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Production of transgenic mice. Transgenic mice were generated by pronuclear microinjection of B6CBF2 or B6D2F2 zygotes as described previously (66). Microinjected zygotes were transferred to the oviduct of NIH or B6CB pseudopregnant females. Transgenic pups were identified by tail genomic DNA PCR as described previously (66). For transgenes 1 to 4 and 6 to 9, primers M329 (5'-GAAGTACGTCATGGGTCACT-3') and M330 (5'-AGCTCCCTCTTGAACTCTAG-3'), which identify a heterologous insertion present in exon 3 of the transgenes (66), were used. For transgenes 5 and 10 to 12, primers TK (5'-CAGCGTCTTGTCATTGGCG-3') and GH (5'-AGTGGTTCGGGGAGTTGGG-3'), which amplify a 170-bp product from the TKGH sequence, were used, and for transgene 13, primers HP3 and PPH9 (5'-TCAAGGGCAAAGAGGAATCA-3'), which amplify an 800-bp, human-specific product, were used. The numbers of independent transgenic lines (in parentheses) analyzed for each transgene were as follows: transgene 1 (1), transgene 2 (2), transgene 3 (2), transgene 4 (1), transgene 5 (3), transgene 6 (1), transgene 7 (2), transgene 8 (2), transgene 9 (5), transgene 10 (2), transgene 11 (1), transgene 12 (5), transgene 13 (3). All animal procedures were conducted in accordance with the Guide for the Care and Use of Laboratory Animals by the U.S. Public Health Service.
Immunohistochemistry. Mice were perfused with 4% paraformaldehyde (PFA) in KPBS (0.9% NaCl, 16 mM K2HPO4, 3.6 mM KH2PO4), and brains were excised, postfixed in 4% PFA-KPBS overnight at 4°C, and sectioned (50 µm) with Vibratome 1000 (Ted Pella, Redding, Calif.). Brain slices were treated with 1% H2O2 in KPBS for 20 min, washed twice with KPBS, and incubated overnight at 4°C with the primary antibodies diluted in KPBS-0.3% Triton X-100 and 2% normal goat serum. The primary antibodies were rabbit polyclonal anti-ACTH-IC-1 (1:1,000; National Hormone and Pituitary Program, National Institutes of Health [NIH]), anti-EGFP ab290 (1:10,000; AbCam), and anti-hGH-IC-3 (1:1,000; National Hormone and Pituitary Program, NIH) antibodies. The next day, slices were washed in KPBS and incubated with biotinylated anti-rabbit immunoglobulin G antibody (Vector) diluted 1:200 in KPBS-0.3% Triton X-100 for 2 h at room temperature (RT). After washing in KPBS, slices were incubated with avidin/biotin-horseradish peroxidase complex (Vectastain Elite ABC kit; Vector) for 1 h at RT, washed in KPBS, and developed with 25 mg of diaminobenzidine (DAB; Sigma)/ml and 0.05% H2O2 in TBS (150 mM NaCl, 50 mM Tris-HCl, pH 7.5). Alternatively, reactions were developed with a blue chromogen (SG kit; Vector). Stained slices were then mounted onto gelatin-coated slides.
Immunofluorescence. To check for colocalization of EGFP and POMC, brain slices containing EGFP-expressing cells were subjected to an anti-ACTH immunohistochemistry as described above up to the second antibody incubation and then slices were incubated in streptavidin conjugated with Cy3 (1:800; Jackson Immunoresearch) in KPBS for 1 h at RT. Subsequently slices were washed in KPBS and mounted onto slides with FluorSave reagent (Calbiochem) and photographed with a microscope (Leica, Wetzlar, Germany) under UV light by using different filters to check for colocalization of green (EGFP) and red (ACTH-Cy3) fluorescence. The protocol for pituitaries was similar except that, prior to immunohistochemistry, perfused pituitaries were dehydrated in an ascending series of sucrose concentrations (10, 20, and 30% in PBS), frozen in embedding medium (Tissue-Tek), and sectioned with a cryostat (Microm HM 505N; Micron, Heidelberg, Germany) and then postfixed in 4% PFA.
Double in situ hybridization and immunohistochemistry. To colocalize endogenous Pomc mRNA and the human GH (hGH) transgene product, nonradioactive POMC in situ hybridization was performed and followed by anti-hGH immunohistochemistry in 12-well plastic plates. A 600-bp mouse POMC exon 3 riboprobe was synthesized using a DIG-labeling kit (Roche) according to the instructions of the manufacturer. PFA-fixed, Vibratome-sectioned brain slices were treated with 0.1 M triethanolamine (TEA), pH 8.0, for 3 min, followed by an incubation in 0.0025% acetic anhydride in 0.1 M TEA for 10 min to block positive charges, and washed twice in 2x SSC (0.3 M NaCl, 0.03 M sodium citrate, pH 7.2). Brain slices were incubated in hybridization solution (50% formamide, 200 mM NaCl, 1x Denhardt's solution, 1 mM EDTA, 10 mM Tris-HCl [pH 8.0], 10% dextran sulfate) for 2 h at 57°C. The riboprobe was heated at 65°C for 5 min (in 100 µl of water with 5 µg of tRNA) and added to the hybridization solution of each well. Hybridization followed overnight at 57°C. The next day, slices were washed for 10 s in 2x SSC and subjected to RNase A digestion (20 µg of RNase A/ml, 50 mM NaCl, 10 mM Tris-HCl [pH 8.0], 1 mM EDTA [pH 8.0]) at 37°C for 30 min, de-salted in 1x and 0.25x SSC for 30 min at 60°C. Slices were washed in buffer 1 (0.1 M Tris-HCl [pH 7.5], 0.15 M NaCl) and incubated with anti-DIG mouse monoclonal antibody (1:5,000; Jackson ImmunoResearch) in buffer 1 containing 0.2% normal goat serum for 2 h at 37°C. Slides were then washed once in buffer 1 and once in KPBS and incubated with goat biotinylated anti-mouse immunoglobulin G antibody (1:200, Vector) for 2 h at RT and then developed with DAB as described above. After that, slides were washed in KPBS and then incubated with anti-hGH (1:1,000) and finally developed with benzidine hydrochloride (Sigma) as described in reference 18.
Radioactive in situ hybridization of pituitaries.
Pituitaries were fresh frozen in embedding medium (Tissue-Tek), and 12-µm-thick sections were collected on gelatinized glass slides with a cryostat microtome (Microm HM 505N; Micron). A 176-nucleotide riboprobe labeled with [
-35S]UTP was synthesized from a plasmid containing part of exon 5 of the hGH gene by using the MAXIscript kit (Ambion). The in situ hybridization was performed as described previously (66).
X-Gal staining. The 50-µm-thick Vibratome brain slices of transgenic mice carrying construct 11 were stained with 1 mg of 5-bromo-4-chloro-3-indolyl-ß-D-glucuronic acid (X-Gal)/ml in staining solution (PBS [pH 7.3] containing 2.12 mg of potassium ferrocyanide/ml, 1.64 mg of potassium ferricyanide/ml, 2 mM MgCl2, 0.01% sodium deoxycolate, and 0.02% NP-40) for 4 h at 37°C. Stained slices were then subjected to anti-ACTH immunohistochemistry as described above.
Isolation and partial mapping of bovine nPE sequences. A genomic clone containing bovine nPE sequences was isolated from a bovine genomic phage library (number BL1015, average insert size of 40 kb; Clontech) by a PCR-based strategy (27) using degenerate primers for a conserved segment of nPE1: PPH7 (5'-GACTGAGCTGAGTGCCTGT-3') and PPH3 (5'-ANGAATGCNGAGTTCTCCA-3'). From this clone, partial sequences of nPE1 and nPE2 were amplified by PCR using primers PPH7 and PPH3 for nPE1 and PPH8 (5'-GGATAAAAGCNGTCTCAAGG-3') and PPH9 for nPE2, subcloned into pGEM-T (Promega), and sequenced. The same clone also contained POMC exon 1 as confirmed by performing PCR using primers specific for the exon 1 of the bovine POMC gene (5-' GGAAGTCCACTCAACGTC-3' and 5'-CTGCGCCCTTACCTGTCA-3'), which were designed based on a published sequence (GenBank accession number J00016).
Sequences and programs. POMC loci sequences from human, mouse, rat, chimpanzee, chick, and zebrafish were retrieved from the Ensembl website (http://www.ensembl.org). Partial dog sequences were retrieved from the National Center for Biotechnology Information (NCBI) Trace database (http://www.ncbi.nih.gov/Traces/trace.cgi). Takifugu rubripes and Xenopus tropicalis sequences were retrieved from the Joint Genome Institute website (http://www.jgi.doe.gov). Human and mouse sequences were compared using the DOTTER program (55) downloaded from http://www.cgb.ki.se/cgb/groups/sonnhammer/Dotter.html, and multiple sequences were compared using the PipMaker program (available at http://bio.cse.psu.edu/pipmaker). Multiple alignments were done using the program Clustal W (58) at http://www.ebi.ac.uk/clustalw. Identification of conserved and aligned transcription factor binding sites was performed using the program rVista2.0 (37) at http://www.dcode.org.
| RESULTS |
|---|
|
|
|---|
|
|
The Pomc neuron-specific region contains two highly conserved elements in mouse and human. Among vertebrate genomes, regulatory elements are often more conserved than average intergenic regions (16, 57). Therefore, in order to identify candidate enhancers within the 13- to 9-kb mouse Pomc fragment, we sequenced this 4-kb region and compared it to the publicly available human POMC locus sequence. Figure 3A shows a global visualization of local sequence alignments (DOTTER program) (55) between the regions from kb 13 to 9 and kb 11 to 7 from the mouse and human POMC genes, respectively. Two regions of approximately 600 and 150 bp with high identities between the two species were readily identifiable. We named these two putative neuronal POMC enhancers nPE1 and nPE2. In addition to the conservation of nucleotide sequence, the distances from each element to the transcriptional start were also similar, indicating that the overall genomic organization at the POMC locus is conserved between these two species (Fig. 3B).
|
|
|
|
nPE1 and nPE2 are conserved in mammals but not in other vertebrates. Since the POMC gene is expressed in the ventromedial hypothalamic neurons of all vertebrate species studied to date, we investigated whether the neuronal regulatory elements nPE1 and nPE2 are conserved in different vertebrate orders. We searched the sequenced POMC loci from the completed genome projects of rat, chimpanzee, and chicken and the teleost fishes pufferfish (T. rubripes) and zebrafish (Danio rerio), as well as the ongoing genome projects of the dog and the frog X. tropicalis (see Materials and Methods). Figure 7A shows a multiple percentage-identity plot (PIP) (53) comparing the human POMC locus to the complete mouse, rat, chimpanzee, chicken, frog, pufferfish and zebra fish POMC loci, as well as to a partial dog sequence. All mammalian sequences showed extensive identity to the human POMC locus, including several conserved noncoding blocks. Among these, nPE1 and nPE2 displayed a particularly greater degree of conservation, evidenced by long stretches and a high percentage of sequence identity relative to those of other less conserved intergenic or intronic regions. The chimpanzee POMC locus was almost 100% identical to the human locus in coding and noncoding regions, as expected given the close evolutionary distance between these two primate species. In contrast, nonmammalian POMC genes did not show significant blocks of sequence identity to the human POMC locus, except for exons 2 and 3 (chicken and Xenopus) and exon 3 (fishes). These data show that both nPE1 and nPE2 are strongly conserved among mammals but cannot be recognized in birds, amphibians, or teleost fishes. Figure 7B shows a multiple alignment of nPE1 and nPE2 of four different mammals, including partial nPE sequences that we amplified by PCR from a bovine genomic phage clone containing the POMC gene. The alignment shows several blocks of high sequence identity, which might indicate possible binding sites for transcription factors. nPE1 had an overall mouse and human similarity of 76%, with a 40-bp core of 100% identity, whereas mouse and human nPE2 had 90% similarity (138 bp out of 153 bp are identical). Interestingly, the similarities between mouse and human exons 1, 2, and 3 of 64, 87, and 82%, respectively, were not higher than the interspecies identity for nPE1 and nPE2.
|
Within nPE2, we identified a sequence that showed high similarity to the consensus binding site for the homeobox gene Nkx6.1 and the POU domain gene Brn2.0. Contiguous to this site, there is a sequence that carries overlapping matches for chicken ovoalbumin upstream promoter (COUP) and estrogen-related receptor alpha (ERR
), two transcription factors that belong to the orphan nuclear receptor family. Another canonical ERR
site was located further downstream.
| DISCUSSION |
|---|
|
|
|---|
In the arcuate nucleus of the mouse hypothalamus, approximately 3,000 POMC-expressing neurons are intermingled with neurons that express other neuropeptides or monoamines, including NPY/AGRP, ghrelin, somatostatin, or dopamine (13, 14, 51). The specific set of neuropeptides expressed in each of these neuronal subpopulations is presumably determined by particular combinations of transcription factors acting during ontogenesis and subsequently in the maintenance of established phenotypes. nPE enhancers probably play an important role in these processes, acting as recognition sites for the assembly of transcriptional complexes that participate in the remodeling of chromatin structure at the POMC locus and the expression of POMC in arcuate neurons. Although the mouse and human regions containing nPE1 and nPE2 proved to be necessary and sufficient to drive reporter gene expression to POMC arcuate neurons of transgenic mice, they were not able to isolate transgenes from the transcriptional influences present at the insertion sites. In consequence, most transgenic pedigrees that showed reporter expression in arcuate POMC neurons also displayed variable patterns of ectopic transgene expression.
Interestingly, our transgenic mouse analysis demonstrated that nPE1 and nPE2 do not need to be simultaneously present to direct transgene expression to POMC arcuate neurons, suggesting some level of functional redundancy between these two enhancers. Nevertheless, it is possible that each nPE enhancer plays a unique role in POMC transcriptional regulation. For example, leptin receptors localized in POMC arcuate neurons signal via phosphorylated STAT3 and deficits in this signal transduction pathway within these neurons correlate with very low levels of POMC mRNA (2, 17). A recent study performed in vitro using a heterologous cell line showed that STAT3 might increase POMC transcription by interacting with a noncanonical STAT3 binding site located immediately upstream of the TATAA box (39). Although it is still not known whether phosphorylation of STAT3 activates Pomc expression directly, nPE1 contains a STAT3 canonical site that is conserved and aligned across mammals and could participate in the leptin regulation of Pomc expression. Further studies are needed to address this possibility.
Modular architecture of the POMC gene.
Our results showed that genomic mouse and human fragments containing nPE1 and nPE2 fused to heterologous minimal promoters direct the expression of transgenes to POMC arcuate neurons in the absence of detectable reporter expression in pituitary melanotrophs or corticotrophs (transgenes 5 and 13). On the other hand, simultaneous deletion of nPE1 and nPE2 in the context of the entire POMC transcriptional unit and flanking regions abolishes reporter gene expression in the arcuate nucleus while maintaining authentic transgenic expression in pituitary POMC cells (transgene 9). Together, these results indicate that two independent sets of enhancers control the cell-specific expression of POMC: a distal module acting in the brain and a proximal module acting in the pituitary. Whether these two modular regions cross talk under particular physiological circumstances to regulate POMC expression in the pituitary or the brain remains to be investigated. Our previous studies demonstrated that all necessary elements to provide cell-specific and hormonally regulated expression of transgenes in pituitary POMC cells are located in approximately 300 bp upstream of the rat or mouse TATAA box (33, 34). Moreover, several transcription factors involved in pituitary POMC expression have been identified, including the homeobox protein Pitx1, the T-box factor Tpit/Tbx19, and basic helix-loop-helix proteins of the NeuroD1 subfamily (31, 32, 44). In contrast, the transcription factors controlling POMC expression in the brain are completely unknown. In this regard, nPE1 and nPE2 could provide a valuable platform for the discovery of such factors. Predictive analysis of binding sites based solely on computer algorithms that identify canonical sequences generally yields a high number of false positives (63) because transcription factors bind to short sequence motifs (6 to 12 bp) that occur very frequently in a genome and different functional DNA elements for the same transcription factor admit a considerable level of degeneracy. However, novel bioinformatics tools have been recently developed that select for canonical binding sites that are aligned within highly conserved sequences of two or more species. Using one of these tools, rVista2.0 (37), we were able to identify five potential transcription factor binding sites along the 600 bp of nPE1 and three sites along the 150 bp of nPE2 that complied with such stringent criteria. The high selective pressure that these sequences withstood during more than 65 million years of evolutionary distance existing between mice and men suggests that these sites play a key role in the commitment of the POMC hypothalamic lineage and the cell-specific and hormone-regulated expression of Pomc in the brain. In fact, nPE1 contains a conserved and aligned site carrying two consecutive canonical sequences for Brn 4.0 and OCT-1, two POU domain transcription factors that participate in the final differentiation of glucagon-expressing pancreatic
cells (26), and GnRH-expressing hypothalamic neurons, respectively (9). This site is adjacent to a consensus sequence for NERF/Elf-2, a member of the Ets transcription factor family. The close proximity of these two sites suggests the existence of a unique combination of a POU domain protein present in POMC neurons that would heterodimerize with an Ets transcription factor to control POMC expression. Both POU domain and Ets factors have been shown to form heterodimers to activate gene expression. In fact, physical interaction of the pituitary-specific POU homeodomain protein Pit-1 with Ets-1 has been shown to be critical to establish the lactotroph pituitary lineage and to induce prolactin gene expression (5). nPE1 also contains a perfect canonical site for STAT3 whose potential significance in the leptin regulation of POMC expression and the central control of satiety mechanisms have been addressed above. In addition, two separate conserved CREB-like sites are present in nPE1, suggesting that transcription factors from the leucine zipper family could also participate in the neuronal regulation of POMC expression.
Within nPE2, there is a conserved and aligned homeodomain binding site carrying a canonical sequence for Nkx6.1 and Brn2.0, two homeodomain transcription factors that participate in the final differentiation of pancreatic ß cells (28) or hypothalamic neuropeptide neurons expressing oxytocin, vasopressin, and CRH (21), respectively. Binding of a homeobox or POU domain factor to this site may be critical to induce neuron-specific expression of POMC. Adjacent to this site, there are overlapping canonical sequences for COUP and ERR
, two transcription factors that belong to the orphan nuclear receptor family. The close proximity of these two binding sites suggests that an orphan nuclear receptor may physically interact with a homeodomain factor to repress POMC expression. Such a type of transcriptional repression has been observed between the homeodomain protein Pit-1 and different members of the orphan nuclear receptor factor family (20). An additional ERR
site is located further downstream at the 3' end of nPE2. ERR
participates in the control of energy balance by regulating the mitochondrial ß-oxidation of fatty acids (62) as well as the development of brown and white adipocytes (38, 64). ERR
-induced transcriptional regulation of hypothalamic POMC could provide another level for the control of energy homeostasis by regulating food intake in concert with leptin and NPY activation of POMC neurons. Taken together, it is tempting to speculate that some, if not all, of these highly conserved and aligned binding sites are the cis-acting platforms where unique clusters of transcription factors present in hypothalamic POMC neurons are recruited to specifically transactivate POMC. Additional functional and biochemical assays will have to be performed in order to finally characterize the cis- and trans-acting partners that control neuronal POMC expression.
The independent modular architecture of the POMC gene may also have important implications in the study of genetic predisposition to familial obesity. Total POMC deficiency causes a complicated syndrome that includes adrenal insufficiency, skin and hair pigmentation deficits, and early-onset obesity probably due to the lack of ACTH and peripheral and central melanocortins, respectively (8, 29, 65). Although homozygous null mutations of POMC are extremely rare in humans, genome-wide QTL scans have demonstrated a genetic linkage between obesity-related traits and a region in human chromosome 2 containing the POMC locus with a reported logarithm of odds score of 7.5, the highest ever found in a human genetic study of obesity-related traits (reviewed in reference 12). This linkage first determined in a population of Mexican-Americans has been reproduced in other studies that analyzed individuals of Caucasian French and African-American origins (23, 25, 47). Although the strongest candidate gene for this chromosome 2 QTL is POMC, no loss-of-function mutations were found within the POMC transcriptional unit (25). Therefore, it is tempting to speculate about the possibility that polymorphisms or mutations in nPE1 and nPE2 could decrease the transcription of POMC and thereby the neural content of POMC peptides, impairing central energy homeostatic mechanisms and predisposing to excessive weight gain without altering the levels or physiological actions of pituitary-released POMC products. Thus, nPE regions are natural candidates to screen for sequence polymorphisms that influence human obesity.
Phylogenetic footprinting to identify neuronal POMC enhancers. Until recently, it was difficult to predict the localization of transcriptional regulatory elements based solely on nucleotide sequence analysis. However, the recent availability of genomic sequences from an increasing number of vertebrate species has allowed the identification of regulatory regions by comparing homologous intergenic regions from different species, a technique termed phylogenetic footprinting (16, 36, 57). This strategy is based on the idea that, as a result of selective pressure, mutations accumulate faster in nonfunctional nucleotide residues whereas functional residues remain more conserved. Therefore, multiple alignment comparisons between orthologous gene loci allow the discovery of potential regulatory sequences embedded in large intergenic regions.
In the study reported here, we combined phylogenetic footprinting with expression analysis in transgenic mice to identify two POMC neuronal enhancers. We found nPE1 and nPE2 in representative species of the mammalian orders Rodentia, Primates, Carnivora, and Artiodactyla, and it is reasonable to predict that nPEs will be found in all mammals. Conversely, we failed to identify conserved nPE enhancers outside Mammalia in spite of the fact that POMC is expressed in the ventromedial hypothalamus of the chicken and teleost fishes like trout, goldfish, and zebrafish (7, 19, 35, 59). Even though there are several examples of enhancers conserved among mammals, birds, and fishes (1, 6, 22, 48, 61), the only conserved regions we observed after comparing the human POMC locus with orthologous genomic sequences from other nonmammalian vertebrates were coding exons 2 and 3. Nevertheless, given that hypothalamic expression of POMC is a conserved vertebrate feature and that core recognition sites for transcription factors are usually very short (6 to 10 bp), it is likely that a similar set of transcription factors control POMC expression in all vertebrates but that the accumulation of neutral mutations throughout evolution prevented us from identifying neuronal POMC enhancers by global alignment programs. The example of POMC illustrates the fact that for each gene, a different set of species might be necessary to obtain useful results by phylogenetic footprinting. In this case, the human X mouse comparison provided an optimal degree of sequence conservation for the identification of functional enhancers while the human X chimpanzee pair was too similar and the sequences of birds, amphibians, or fishes was too different.
| ACKNOWLEDGMENTS |
|---|
This work was supported in part by a Fogarty International Research Collaborative Award TW01233 (M.J.L. and M.R.), NIH grant DK68400 (M.J.L. and M.R.), International Research Scholar Grant of the Howard Hughes Medical Institute (M.R.), Agencia Nacional de Promoción Científica y Tecnológica (M.R.), JS Guggenheim Foundation (M.R.) and Universidad de Buenos Aires (M.R.). F. S. J. de Souza, A. M. Santangelo, V. F. Bumaschny, and M. E. Avale are recipients of fellowships from CONICET, Argentina.
M.J.L. declares a financial interest in transgenic models utilizing the POMC neural regulatory regions, and this potential conflict of interest has been disclosed to and managed by the OHSU Office for Integrity in Research.
| FOOTNOTES |
|---|
F.S.J.D.S. and A.M.S. contributed equally to this work. ![]()
| REFERENCES |
|---|
|
|
|---|
2. Bates, S. H., W. H. Stearns, T. A. Dundon, M. Schubert, A. W. Tso, Y. Wang, A. S. Banks, H. J. Lavery, A. K. Haq, E. Maratos-Flier, B. G. Neel, M. W. Schwartz, and M. G. Myers, Jr. 2003. STAT3 signalling is required for leptin regulation of energy balance but not reproduction. Nature 421:856-859.[CrossRef][Medline]
3. Batterham, R. L., M. A. Cowley, C. J. Small, H. Herzog, M. A. Cohen, C. L. Dakin, A. M. Wren, A. E. Brynes, M. J. Low, M. A. Ghatei, R. D. Cone, and S. R. Bloom. 2002. Gut hormone PYY(3-36) physiologically inhibits food intake. Nature 418:650-654.[CrossRef][Medline]
4. Bjorbaek, C., and B. B. Kahn. 2004. Leptin signaling in the central nervous system and the periphery. Recent Prog. Horm. Res. 59:305-331.
5. Bradford, A. P., K. S. Brodsky, S. E. Diamond, L. C. Kuhn, Y. Liu, and A. Gutierrez-Hartmann. 2000. The Pit-1 homeodomain and beta-domain interact with Ets-1 and modulate synergistic activation of the rat prolactin promoter. J. Biol. Chem. 275:3100-3106.
6. Brenner, S., B. Venkatesh, W. H. Yap, C. F. Chou, A. Tay, S. Ponniah, Y. Wang, and Y. H. Tan. 2002. Conserved regulation of the lymphocyte-specific expression of lck in the Fugu and mammals. Proc. Natl. Acad. Sci. USA 99:2936-2941.
7. Cerda-Reverter, J. M., H. B. Schioth, and R. E. Peter. 2003. The central melanocortin system regulates food intake in goldfish. Regul. Pept. 115:101-113.[CrossRef][Medline]
8. Challis, B. G., A. P. Coll, G. S. Yeo, S. B. Pinnock, S. L. Dickson, R. R. Thresher, J. Dixon, D. Zahn, J. J. Rochford, A. White, R. L. Oliver, G. Millington, S. A. Aparicio, W. H. Colledge, A. P. Russ, M. B. Carlton, and S. O'Rahilly. 2004. Mice lacking pro-opiomelanocortin are sensitive to high-fat feeding but respond normally to the acute anorectic effects of peptide-YY(3-36). Proc. Natl. Acad. Sci. USA 101:4695-4700.
9. Clark, M. E., and P. L. Mellon. 1995. The POU homeodomain transcription factor Oct-1 is essential for activity of the gonadotropin-releasing hormone neuron-specific enhancer. Mol. Cell. Biol. 15:6169-6177.[Abstract]
10. Coll, A. P., I. S. Farooqi, B. G. Challis, G. S. Yeo, and S. O'Rahilly. 2004. Proopiomelanocortin and energy balance: insights from human and murine genetics. J. Clin. Endocrinol. Metab. 89:2557-2562.
11. Comuzzie, A. G., J. E. Hixson, L. Almasy, B. D. Mitchell, M. C. Mahaney, T. D. Dyer, M. P. Stern, J. W. MacCluer, and J. Blangero. 1997. A major quantitative trait locus determining serum leptin levels and fat mass is located on human chromosome 2. Nat. Genet. 15:273-276.[CrossRef][Medline]
12. Comuzzie, A. G. 2002. The emerging pattern of the genetic contribution to human obesity. Best. Pract. Res. Clin. Endocrinol. Metab. 16:611-621.[CrossRef][Medline]
13. Cowley, M. A., J. L. Smart, M. Rubinstein, M. G. Cerdan, S. Diano, T. L. Horvath, R. D. Cone, and M. J. Low. 2001. Leptin activates anorexigenic POMC neurons through a neural network in the arcuate nucleus. Nature 411:480-484.[CrossRef][Medline]
14. Cowley, M. A., R. G. Smith, S. Diano, M. Tschop, N. Pronchuk, K. L. Grove, C. J. Strasburger, M. Bidlingmaier, M. Esterman, M. L. Heiman, L. M. Garcia-Segura, E. A. Nillni, P. Mendez, M. J. Low, P. Sotonyi, J. M. Friedman, H. Liu, S. Pinto, W. F. Colmers, R. D. Cone, and T. L. Horvath. 2003. The distribution and mechanism of action of ghrelin in the CNS demonstrates a novel hypothalamic circuit regulating energy homeostasis. Neuron 37:649-661.[CrossRef][Medline]
15. Delplanque, J., M. Barat-Houari, C. Dina, P. Gallina, K. Clement, B. Guy-Grand, F. Vasseur, P. Boutin, and P. Froguel. 2000. Linkage and association studies between the proopiomelanocortin (POMC) gene and obesity in caucasian families. Diabetologia 43:1554-1557.[CrossRef][Medline]
16. Fickett, J. W., and W. W. Wasserman. 2000. Discovery and modeling of transcriptional regulatory regions. Curr. Opin. Biotechnol. 11:19-24.[CrossRef][Medline]
17. Gao, Q., M. J. Wolfgang, S. Neschen, K. Morino, T. L. Horvath, G. I. Shulman, and X. Y. Fu. 2004. Disruption of neural signal transducer and activator of transcription 3 causes obesity, diabetes, infertility, and thermal dysregulation. Proc. Natl. Acad. Sci. USA 101:4661-4666.
18. Gelman, D. M., D. Noain, M. E. Avale, V. Otero, M. J. Low, and M. Rubinstein. 2003. Transgenic mice engineered to target Cre/loxP-mediated DNA recombination into catecholaminergic neurons. Genesis 36:196-202.[CrossRef][Medline]
19. Gerets, H. H., K. Peeters, L. Arckens, F. Vandesande, and L. R. Berghman. 2000. Sequence and distribution of pro-opiomelanocortin in the pituitary and the brain of the chicken (Gallus gallus). J. Comp. Neurol. 417:250-262.[CrossRef][Medline]
20. Gonzalez, M. M., and C. Carlberg. 2002. Cross-repression, a functional consequence of the physical interaction of non-liganded nuclear receptors and POU domain transcription factors. J. Biol. Chem. 277:18501-18509.
21. Goshu, E., H. Jin, J. Lovejoy, J. F. Marion, J. L. Michaud, and C. M. Fan. 2004. Sim2 contributes to neuroendocrine hormone gene expression in the anterior hypothalamus. Mol. Endocrinol. 18:1251-1262.
22. Göttgens, B., L. M. Barton, J. G. Gilbert, A. J. Bench, M. J. Sanchez, S. Bahn, S. Mistry, D. Grafham, A. McMurray, M. Vaudin, E. Amaya, D. R. Bentley, A. R. Green, and A. M. Sinclair. 2000. Analysis of vertebrate SCL loci identifies conserved enhancers. Nat. Biotechnol. 18:181-186.[CrossRef][Medline]
23. Hager, J., C. Dina, S. Francke, S. Dubois, M. Houari, V. Vatin, E. Vaillant, N. Lorentz, A. Basdevant, K. Clement, B. Guy-Grand, and P. Froguel. 1998. A genome-wide scan for human obesity genes reveals a major susceptibility locus on chromosome 10. Nat. Genet. 20:304-308.[CrossRef][Medline]
24. Heisler, L. K., M. A. Cowley, L. H. Tecott, W. Fan, M. J. Low, J. L. Smart, M. Rubinstein, J. B. Tatro, J. N. Marcus, H. Holstege, C. E. Lee, R. D. Cone, and J. K. Elmquist. 2002. Activation of central melanocortin pathways by fenfluramine. Science 297:609-611.
25. Hixson, J. E., L. Almasy, S. Cole, S. Birnbaum, B. D. Mitchell, M. C. Mahaney, M. P. Stern, J. W. MacCluer, J. Blangero, and A. G. Comuzzie. 1999. Normal variation in leptin levels is associated with polymorphisms in the proopiomelanocortin gene, POMC. J. Clin. Endocrinol. Metab. 84:3187-3191.
26. Hussain, M. A., C. P. Miller, and J. F. Habener. 2002. Brn-4 transcription factor expression targeted to the early developing mouse pancreas induces ectopic glucagon gene expression in insulin-producing beta cells. J. Biol. Chem. 277:16028-16032.
27. Israel, D. I. 1993. A PCR-based method for high stringency screening of DNA libraries. Nucleic Acids Res. 21:2627-2631.
28. Iype, T., D. G. Taylor, S. M. Ziesmann, J. C. Garmey, H. Watada, and R. G. Mirmira. 2004. The transcriptional repressor Nkx6.1 also functions as a deoxyribonucleic acid context-dependent transcriptional activator during pancreatic beta-cell differentiation: evidence for feedback activation of the nkx6.1 gene by Nkx6.1. Mol. Endocrinol. 18:1363-1375.
29. Krude, H., H. Biebermann, W. Luck, R. Horn, G. Brabant, and A. Gruters. 1998. Severe early-onset obesity, adrenal insufficiency and red hair pigmentation caused by POMC mutations in humans. Nat. Genet. 19:155-157.[CrossRef][Medline]
30. Krude, H., H. Biebermann, D. Schnabel, M. Z. Tansek, P. Theunissen, P. E. Mullis, and A. Gruters. 2003. Obesity due to proopiomelanocortin deficiency: three new cases and treatment trials with thyroid hormone and ACTH4-10. J. Clin. Endocrinol. Metab. 88:4633-4640.
31. Lamolet, B., A. M. Pulichino, T. Lamonerie, Y. Gauthier, T. Brue, A. Enjalbert, and J. Drouin. 2001. A pituitary cell-restricted T box factor, Tpit, activates POMC transcription in cooperation with Pitx homeoproteins. Cell 104:849-859.[CrossRef][Medline]
32. Lamonerie, T., J. J. Tremblay, C. Lanctot, M. Therrien, Y. Gauthier, and J. Drouin. 1996. Ptx1, a bicoid-related homeo box transcription factor involved in transcription of the pro-opiomelanocortin gene. Genes Dev. 10:1284-1295.
33. Liu, B., G. D. Hammer, M. Rubinstein, M. Mortrud, and M. J. Low. 1992. Identification of DNA elements cooperatively activating proopiomelanocortin gene expression in the pituitary glands of transgenic mice. Mol. Cell. Biol. 12:3978-3990.
34. Liu, B., M. Mortrud, and M. J. Low. 1995. DNA elements with AT-rich core sequences direct pituitary cell-specific expression of the pro-opiomelanocortin gene in transgenic mice. Biochem. J. 312:827-832.
35. Liu, N. A., H. Huang, Z. Yang, W. Herzog, M. Hammerschmidt, S. Lin, and S. Melmed. 2003. Pituitary corticotroph ontogeny and regulation in transgenic zebrafish. Mol. Endocrinol. 17:959-966.
36. Loots, G. G., R. M. Locksley, C. M. Blankespoor, Z. E. Wang, W. Miller, E. M. Rubin, and K. A. Frazer. 2000. Identification of a coordinate regulator of interleukins 4, 13, and 5 by cross-species sequence comparisons. Science 288:136-140.
37. Loots, G. G., and I. Ovcharenko. 2004. rVISTA 2.0: evolutionary analysis of transcription factor binding sites. Nucleic Acids Res. 32:W217-W221.
38. Luo, J., R. Sladek, J. Carrier, J.-A. Bader, D. Richard, and V. Giguère. 2003. Reduced fat mass in mice lacking orphan nuclear receptor estrogen-related receptor
. Mol. Cell. Biol. 23:7947-7956.
39. Münzberg, H., L. Huo, E. A. Nillni, A. N. Hollenberg, and C. Bjorbaek. 2003. Role of signal transducer and activator of transcription 3 in regulation of hypothalamic proopiomelanocortin gene expression by leptin. Endocrinology 144:2121-2131.
40. Newell-Price, J. 2003. Proopiomelanocortin gene expression and DNA methylation: implications for Cushing's syndrome and beyond. J. Endocrinol. 177:365-372.[Abstract]
41. Overstreet, L. S., S. T. Hentges, V. F. Bumaschny, F. S. de Souza, J. L. Smart, A. M. Santangelo, M. J. Low, G. L. Westbrook, and M. Rubinstein. 2004. A transgenic marker for newly born granule cells in dentate gyrus. J. Neurosci. 24:3251-3259.
42. Pfeffer, P. L., M. Bouchard, and M. Busslinger. 2000. Pax2 and homeodomain proteins cooperatively regulate a 435 bp enhancer of the mouse Pax5 gene at the midbrain-hindbrain boundary. Development 127:1017-1028.[Abstract]
43. Philips, A., M. Maira, A. Mullick, M. Chamberland, S. Lesage, P. Hugo, and J. Drouin. 1997. Antagonism between Nur77 and glucocorticoid receptor for control of transcription. Mol. Cell. Biol. 17:5952-5959.[Abstract]
44. Poulin, G., B. Turgeon, and J. Drouin. 1997. NeuroD1/beta2 contributes to cell-specific transcription of the proopiomelanocortin gene. Mol. Cell. Biol. 17:6673-6682.[Abstract]
45. Pulichino, A. M., S. Vallette-Kasic, J. P. Tsai, C. Couture, Y. Gauthier, and J. Drouin. 2003. Tpit determines alternate fates during pituitary cell differentiation. Genes Dev. 17:738-747.
46. Raffin-Sanson, M. L., Y. de Keyzer, and X. Bertagna. 2003. Proopiomelanocortin, a polypeptide precursor with multiple functions: from physiology to pathological conditions. Eur. J. Endocrinol. 149:79-90.[Abstract]
47. Rotimi, C. N., A. G. Comuzzie, W. L. Lowe, A. Luke, J. Blangero, and R. S. Cooper. 1999. The quantitative trait locus on chromosome 2 for serum leptin levels is confirmed in African-Americans. Diabetes 48:643-644.[Medline]
48. Rowitch, D. H., Y. Echelard, P. S. Danielian, K. Gellner, S. Brenner, and A. P. McMahon. 1998. Identification of an evolutionarily conserved 110 base-pair cis-acting regulatory sequence that governs Wnt-1 expression in the murine neural plate. Development 125:2735-2746.[Abstract]
49. Rubinstein, M., M. Mortrud, B. Liu, and M. J. Low. 1993. Rat and mouse proopiomelanocortin gene sequences target tissue-specific expression to the pituitary gland but not to the hypothalamus of transgenic mice. Neuroendocrinology 58:373-380.[Medline]
50. Rubinstein, M., J. S. Mogil, M. Japon, E. C. Chan, R. G. Allen, and M. J. Low. 1996. Absence of opioid stress-induced analgesia in mice lacking beta-endorphin by site-directed mutagenesis. Proc. Natl. Acad. Sci. USA 93:3995-4000.
51. Saper, C. B., T. C. Chou, and J. K. Elmquist. 2002. The need to feed: homeostatic and hedonic control of eating. Neuron 36:199-211.[CrossRef][Medline]
52. Schwartz, M. W., S. C. Woods, D. Porte, Jr., R. J. Seeley, and D. G. Baskin. 2000. Central nervous system control of food intake. Nature 404:661-671.[Medline]
53. Schwartz, S., Z. Zhang, K. A. Frazer, A. Smit, C. Riemer, J. Bouck, R. Gibbs, R. Hardison, and W. Miller. 2000. PipMakera web server for aligning two genomic DNA sequences. Genome Res. 10:577-586.
54. Scully, K. M., and M. G. Rosenfeld. 2002. Pituitary development: regulatory codes in mammalian organogenesis. Science 295:2231-2235.
55. Sonnhammer, E. L., and R. Durbin. 1995. A dot-matrix program with dynamic threshold control suited for genomic DNA and protein sequence analysis. Gene 167:GC1-GC10.[CrossRef][Medline]
56. Spiegelman, B. M., and J. S. Flier. 2001. Obesity and the regulation of energy balance. Cell 104:531-543.[CrossRef][Medline]
57. Thomas, J. W., and J. W. Touchman. 2002. Vertebrate genome sequencing: building a backbone for comparative genomics. Trends Genet. 18:104-108.[CrossRef][Medline]
58. Thompson, J. D., D. G. Higgins, and T. J. Gibson. 1994. CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 22:4673-4680.
59. Tollemer, H., C. A. Teitsma, J. Leprince, T. Bailhache, F. Vandesande, O. Kah, M. C. Tonon, and H. Vaudry. 1999. Immunohistochemical localization and biochemical characterization of two novel decapeptides derived from POMC-A in the trout hypothalamus. Cell Tissue Res. 295:409-417.[CrossRef][Medline]
60. Tremblay, J. J., C. Lanctot, and J. Drouin. 1998. The pan-pituitary activator of transcription, Ptx1 (pituitary homeobox 1), acts in synergy with SF-1 and Pit1 and is an upstream regulator of the Lim-homeodomain gene Lim3/Lhx3. Mol. Endocrinol. 12:428-441.
61. Uchikawa, M., Y. Ishida, T. Takemoto, Y. Kamachi, and H. Kondoh. 2003. Functional analysis of chicken Sox2 enhancers highlights an array of diverse regulatory elements that are conserved in mammals. Dev. Cell 4:509-519.[CrossRef][Medline]
62. Vega, R. B., and D. P. Kelly. 1997. A role for estrogen-related receptor alpha in the control of mitochondrial fatty acid beta-oxidation during brown adipocyte differentiation. J. Biol. Chem. 272:31693-31699.
63. Wasserman, W. W., and A. Sandelin. 2004. Applied bioinformatics for the identification of regulatory elements. Nat. Rev. Genet. 5:276-287.[CrossRef][Medline]
64. Willy, P. J., I. R. Murray, J. Qian, B. B. Busch, W. C. Stevens, Jr., R. Martin, R. Mohan, S. Zhou, P. Ordentlich, P. Wei, D. W. Sapp, R. A. Horlick, R. A. Heyman, and I. G. Schulman. 2004. Regulation of PPARgamma coactivator 1alpha (PGC-1alpha) signaling by an estrogen-related receptor alpha (ERRalpha) ligand. Proc. Natl. Acad. Sci. USA 101:8912-8917.