| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
,
Department of Molecular Physiology and Biophysics, Center for Human Genetics Research, Vanderbilt University School of Medicine, 519 Light Hall, 2215 Garland Avenue, Nashville, Tennessee 37232,1 Department of Biological Sciences, Vanderbilt University, 2401 Stevenson Center, 465 21st Avenue South, Nashville, Tennessee 372322
Received 28 August 2006/ Returned for modification 16 October 2006/ Accepted 22 January 2007
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
The expression of Bmp2 during development has been extensively characterized, and it plays many roles in pattern formation and morphogenesis, cell fate determination, and/or differentiation. Bmp2 is expressed dynamically in numerous embryonic locations, such as kidneys, hair follicles, tooth buds, and gut epithelium (4, 16, 43, 51, 52). The involvement of Bmp2 in limb patterning is evidenced by its expression in the apical ectodermal ridge (AER) and limb bud mesenchyme (28, 50). From 8.5 to 9.5 days postcoitus (dpc), Bmp2 is transcribed in extraembryonic tissues and heart mesoderm (51, 75), and the early embryonic lethality observed in Bmp2 knockout mice is probably due to defects in these tissues (75). There is extensive evidence for the involvement of Bmp2 in bone and cartilage formation. In developing endochondral bones, Bmp2 is expressed in hypertrophic chondrocytes of the growth plates and also in the osteogenic perichondrium (52, 60, 65). Bmp2 can also stimulate signaling cascades that are required for osteoblast differentiation (9). The recently reported genetic linkage and association of osteoporosis (OP) to the human BMP2 gene also supports its role in human bone homeostasis and disease (67). In addition, Bmp2 is of tremendous interest as a therapeutic agent for OP, fracture healing, and other ailments for which increased bone formation would be beneficial (9, 55, 69).
Despite the potential importance of Bmp2 for the proper development of bone and other tissues, little is known about the molecular mechanisms that control when and where Bmp2 is expressed in vivo. Previous reports have described the regulatory functions of Bmp2 promoter fragments in cell lines (1, 18, 23, 27, 68, 76), but these provide only limited information regarding how Bmp2 transcription is controlled in vivo. A transgene containing a 2.7-kilobase Bmp2 promoter fragment has been reported to drive reporter gene expression in skeletal perichondrium and chondrocytes in vivo (18), suggesting that some important elements for Bmp2 skeletal expression lie close to the promoter. However, the same transgene failed to confer many other sites of Bmp2 expression (S. Harris, personal communication). On the other hand, several lines of evidence suggest Bmp2 might be controlled by multiple, long-range regulatory elements outside the promoter region. First, the complex patterns of Bmp2 expression are reminiscent of those for Bmp5, Gdf6, and dpp, which are other BMP genes that each have multiple, distant cis-regulatory elements (5, 12, 13, 53). Second, Bmp2 is flanked by "stable gene deserts" (i.e., large genomic regions that are devoid of other genes yet conserved across species) (58), suggesting that evolution has preserved an arrangement of distant cis-regulatory elements flanking Bmp2.
To test the hypothesis that Bmp2 is regulated by distant cis-acting sequences in vivo, we used bacterial artificial chromosome (BAC) transgenes to survey a region of approximately 400 kilobases surrounding the mouse Bmp2 gene. To do this we engineered Bmp2 BAC clones to contain a lacZ reporter cassette and then tested the BACs in transgenic embryos for their ability to direct normal patterns of Bmp2 expression in mid-gestation mouse embryos. Here we show that individual BAC clones that extend far 5' or 3' to the Bmp2 gene can direct distinct subsets of Bmp2 expression patterns in vivo, strongly suggesting that different cis-acting regulatory sequences are located in separate genomic regions flanking either side of Bmp2. Surprisingly, cis-regulatory sequences that control Bmp2 expression in hypertrophic chondrocytes are very far from those controlling expression in osteoblast progenitors, which are in a distant 3' flanking region. We also found that in the osteogenic perichondrium of the endochondral bones, Bmp2 expression is closely preceded temporally by the key osteoblast factor Runx2. By testing additional Bmp2 BACs with engineered deletions across the 3' region, we more closely refined the location of the osteoblast progenitor regulatory element, and we could also parse distinct elements controlling Bmp2 expression in tooth buds, hair placodes, kidney, and other tissues into four separate genomic regions 3' to Bmp2. Finally, using heterologous minigenes, we identified a conserved osteoblast progenitor Bmp2 enhancer located approximately 156 kb from the promoter. Taken together, these findings demonstrate that numerous and distant regulatory sequences direct Bmp2 expression during development and that Bmp2 transcription in differentiating osteoblasts is controlled by a distant 3' enhancer.
| MATERIALS AND METHODS |
|---|
|
|
|---|
BAC IRES-ß-geo modifications. An internal ribosome entry site (IRES)-ß-geo cassette was inserted into each Bmp2 BAC, in place of Bmp2 exon 3 mature-region-coding sequences, using the homologous recombination technique of Lee et al. (47) as follows. The plasmid pIBG-FTet was generated by ligating the IRES-ß-geo-SV40pA cassette from pGT1.8 (54) and an FRT-flanked tetracycline resistance cassette into a modified pBluescript II SK(+) backbone. The recombination cassette was constructed by subcloning a 50-bp 5' recombination arm overlapping part of Bmp2 intron 2 and part of exon 3 and a 50-bp 3' recombination arm containing part of exon 3 into pIBG-FTet such that the recombination arms flanked the IRES-ß-geo-FRT-Tet-FRT cassette. The forward strand (relative to Bmp2) sequences of the 50-bp homology arms were as follows: for the 5' arm, TTGCACCAAGATGAACACAGCTGGTCACAGATAAGGCCATTGCTAGTGACT; for the 3' arm, GTGTCGTTAGCACAGCAAGAATAAATAAATAAATATATATATTTTAGAAAC. Both recombination arms were created by annealing polyacrylamide gel electrophoresis (PAGE)-purified oligonucleotides designed to allow overhanging ends for direct ligation to pIBG-FTet. The final cassette with recombination arms was digested from the vector, gel purified, and then recombined with Bmp2 BACs as described previously (47). Successful recombinants were selected by tetracycline resistance. The tetracycline resistance gene was then removed by L-arabinose promoter-driven FLpe recombinase excision (47, 53). Correctly modified BACs were verified by conventional and pulsed-field gel analysis of restriction digests to confirm expected banding patterns, as well as direct BAC sequencing.
Bmp2 3' deletion lacZ-BAC modifications.
Four 3' deletion lacZ-BACs were engineered using the mouse 3' BAC (RP23-409L24) that contains the lacZ insertion (as described above) and the galK homologous recombination technique of Warming et al. (72) as follows. First, the galK gene was used as a selection cassette and inserted into the 3' lacZ-BAC. Selection cassettes used for targeting the 3' lacZ-BAC deletions were generated by PCR amplification of the galK gene with PAGE-purified primers, which also contain homology flanking the genomic regions deleted. The primer sequences (5'
3') used are as follows (homology arms are in uppercase, and sequences that anneal within galK gene are in lowercase): for deletion 1 (D1), TGGGCACAGAGAGCAGAGCTAACAGTTGAGTGGGGAGTGTCTAAGAATTCcctgttgacaattaatcatcggca (forward) and CAATTGTATATACCCAGTGACTTTCCCTGATTCACTGGGGAGACATGGGAtcagcactgtcctgctcctt (reverse); for D2, TGTGTGTGTGTGTGTGCTTTTTTAGACTAAGGTTGGCAATGTTCCTAAGGcctgttgacaattaatcatcggca (forward) and AATGGGGAGGTGGACTTGATAGATGTGTACACTCTGCCTATTGTGTACATtcagcactgtcctgctcctt (reverse); for D3, CACAATATCAAAATGTGTCAATGGTATGAAGTAGTAGAGCTCATAATGACcctgttgacaattaatcatcggca (forward) and AGTAAGATAGGTGCACATAGAGAAAGGAGAATTGGGACAGACGGACtcagcactgtcctgctcctt (reverse); for D4, CCAAGCCCAGCTGCTTCTGACTGGGTATTCTCTGGTGACCTGTTCTTTGAcctgttgacaattaatcatcggca (forward) and ATTTGCGGCCGCTAATACGACTCACTATAGGGAGAGGATCCGCGGAATTCtcagcactgtcctgctcctt (reverse). Note that the 3' homology arms used to generate D4 have homology to sequences present in the BAC vector. The final galK cassette with recombination arms was gel purified and then recombined with 3' lacZ-BAC as described previously (72). Successful recombinants were selected for the ability to use galactose as a carbon source in minimal medium. Second, the galK cassettes were deleted via recombination, by using synthetic DNA fragments made by annealing PAGE-purified oligonucleotides that were homologous to sequences flanking the galK insertion site. The forward-strand (relative to Bmp2) sequences of the annealed oligonucleotides used for substitutions are as follows: for D1, TGGGCACAGAGAGCAGAGCTAACAGTTGAGTGGGGAGTGTCTAAGAATTCTCCCATGTCTCCCCAGTGAATCAGGGAAAGTCACTGGGTATATACAATTG; for D2, TGTGTGTGTGTGTGTGCTTTTTTAGACTAAGGTTGGCAATGTTCCTAAGGATGTACACAATAGGCAGAGTGTACACATCTATCAAGTCCACCTCCCCATT; for D3, CACAATATCAAAATGTGTCAATGGTATGAAGTAGTAGAGCTCATAATGACATGAGTCCGTCTGTCCCAATTCTCCTTTCTCTATGTGCACCTATCTTACT; for D4, CCAAGCCCAGCTGCTTCTGACTGGGTATTCTCTGGTGACCTGTTCTTTGAGAATTCCGCGGATCCTCTCCCTATAGTGAGTCGTATTAGCGGCCGCAAAT. Annealed oligonucleotides with recombination arms were recombined with 3' galK+ deletion lacZ-BACs as described previously (72). Successful galK-negative recombinants were selected for resistance to 2-deoxygalactose in minimal medium with glycerol as the carbon source. Correctly modified deletion BACs were verified by conventional and pulsed-field gel analysis of restriction digests to confirm expected banding patterns, as well as direct BAC sequencing.
ECR-Hsp promoter-lacZ transgene construction. Evolutionarily conserved region (ECR)-Hsp68 promoter-lacZ transgenes were made as follows. PCR products corresponding to genomic sequences for ECR-1 + 2 (4.5 kb), ECR-1 (656 bp), and ECR-2 (358 bp) were amplified using 3' BAC DNA (RP23-409L24) and the following primer sequences: for ECR-1 + 2, AGAGAGAGCTAGATGGTGGAC (forward) and TGCCTTAGTCCTCTCCAGAG (reverse); for ECR-1, CGCTGCTTTTGAATGCTT (forward) and TAGACCCAAGATGGAGTATTAG (reverse); for ECR-2, TGTACTTACTCGATTGTGTGC (forward) and TGCCTTAGTCCTCTCCAGAG (reverse). Each PCR product was ligated into the pGEM-T Easy vector (Promega). A NotI restriction fragment containing each ECR fragment was ligated into pSfi-Hsp68lacZ in the forward orientation (13, 53, 63) to create ECR-Hsp68 promoter-lacZ transgenes. All ECR plasmids were verified by analysis of restriction digests and direct sequencing.
Transgenic mice. Bmp2 BAC DNAs were purified according to established techniques and used for pronuclear injection of C57BL/6J x DBA/2J F1 hybrid embryos (12). All BACs were injected as uncut circular DNAs. Circular BACs have been shown to be as efficient as linearized BACs for generating transgenic mice (24, 53). Plasmids used for microinjection were prepared as follows. Plasmid DNA was isolated using the QIAGEN Endo-Free Maxi plasmid isolation kit. One hundred micrograms of plasmid DNA was linearized with SalI, and the enzyme was removed using Micropure EZ restriction enzyme removal columns (Millipore). The filtrate containing the DNA was then transferred to Montage PCR filter units (Millipore) and purified according to the manufacturer's protocol. Purified and linearized DNAs were then dialyzed against 3 liters of microinjection buffer (10 mM Tris-HCl [pH 7.4], 0.15 mM EDTA [pH 8.0]) using 10,000-molecular-weight-cutoff Slide-A-Lyzer dialysis cassettes (Pierce). DNAs were further dialyzed and concentrated with 30,000-molecular-weight-cutoff Centriprep centrifugal filter devices (Millipore) prior to microinjection. Injections and oviduct transfers were performed by the Vanderbilt Transgenic Core Facility using standard techniques in accordance with protocols approved by the Vanderbilt University Institutional Animal Care and Use Committee. Transgenic embryos or weanlings were verified by PCR from yolk sac or tail DNAs. All of the BAC-transgenic embryos and lines analyzed for lacZ expression in this study were positive for the ß-geo insertion and BAC vector, as determined by both conventional and real-time PCRs.
Bmp2 BAC transgene copy number determination.
Genomic DNAs were harvested from mouse embryo yolk sacs or tail samples from live-born transgenic mice by using standard genomic DNA isolation methods. BAC transgene copy number was estimated using TaqMan-based (Applied Biosystems Inc.) real-time PCR on a 7900HT platform (Applied Biosystems Inc.). All reactions were performed according to the manufacturer's instructions, using 20 ng of mouse genomic DNA. Custom TaqMan primer/probe assays for the neomycin (Neo) resistance gene (ß-geo insertion) and the chloramphenicol (Cam) cassette (BAC vector) were provided by Applied Biosystems (Neo and Cam primer/probe sequences are available upon request). For all mouse genomic DNA samples, the Jun gene (GenBank accession no. NM_010591.1; Applied Biosystems Inc. assay identification no. Mm00495062_s1) was used as an internal control. Copy number estimates for either Neo or Cam were derived using the triplicate
CT value relative to the internal Jun genomic control and compared to the
CT values of standard-curve samples derived from C57BL/6J x DBA/2J F1 hybrid genomic DNAs spiked with known quantities of BAC DNAs diluted to copy number equivalents (1, 2, 4, 8, and 16 BAC copies per diploid genome).
Preparation of agarose-embedded high-molecular-weight DNA from BAC-transgenic embryos. Mouse embryonic fibroblasts (MEFs) used for the generation of agarose-embedded high-molecular-weight DNA were isolated and cultured as follows. Embryonic day 13.5 embryos were generated by crossing BAC-transgenic males with wild-type CD-1 females. Embryos were dissected at room temperature in sterile phosphate-buffered saline (PBS), and the head and viscera were removed. The head was X-Gal (5-bromo-4-chloro-3-indolyl-ß-D-galactopyranoside) stained in parallel to identify transgenic embryos. Individual embryos were then forced through an 18-gauge needle into a 10-cm tissue culture dish containing 10 ml of MEF medium (Dulbecco's modified Eagle medium supplemented with 4.5 g/liter glucose, 2 mM sodium pyruvate, 0.5 g/liter L-glutamine, 10% fetal bovine serum, 200 units/liter of penicillin G, 200 µg/liter streptomycin sulfate, and 0.5 µg/liter of amphotericin B). Cells were incubated at 37°C in a 5% CO2 atmosphere and fed every 2 days with MEF medium until they reached confluence. Confluent cells were passaged at 1:3 dilutions with trypsin and propagated in MEF medium until sufficient cell numbers were attained. Confluent MEF cells were washed twice in sterile 1x PBS, disassociated with a rubber policeman, and resuspended in sterile cell suspension buffer (0.1 M EDTA [pH 8.0], 10 mM Tris-HCl [pH 7.6], 20 mM NaCl) at a concentration of approximately 7 x 106 cells/ml. Cells were embedded in 0.75% low-melting-point agarose (Invitrogen) using agarose plug molds (Bio-Rad). Approximately 10 plugs were then placed in 10 ml of lysis buffer (1% sodium dodecyl sulfate, 2 mg/ml proteinase K [Roche], 0.5 M EDTA [pH 8.0]) and incubated overnight in a 50°C water bath with no agitation. Plugs were sequentially washed in 50 ml of the following buffers at room temperature for 30 min each with agitation: twice in 0.1 EX buffer (0.1 M EDTA [pH 8.0], 0.01% Triton X-100), once in 10:10 buffer (10 mM Tris-HCl [pH 8.0], 10 mM EDTA [pH 8.0]), once in proteinase K inactivation buffer (10 mM Tris-HCl [pH 8.0], 10 mM EDTA [pH 8.0], 1 mM phenylmethylsulfonyl fluoride), twice in 10:10 buffer, twice in 0.5 EX buffer (0.5 M EDTA [pH 8.0], 0.01% Triton X-100), and once in 0.5 M EDTA (pH 8.0). The plugs were then stored until use at 4°C in 0.5 M EDTA (pH 8.0).
Southern analysis of high-molecular-weight transgenic DNA.
Agarose plugs (isolated as described above) were washed twice in 50 ml of TEX buffer (10 mM Tris-HCl [pH 8.0], 1 mM EDTA [pH 8.0], 0.01% Triton X-100) and once in 50 ml of Tris-EDTA (pH 8.0) at room temperature with agitation. Plugs were then transferred to 2-ml screw cap tubes (one plug per tube) and equilibrated for 30 min at room temperature in 2 ml of 1x restriction enzyme buffer with agitation. For NotI enzyme, 10 mM spermidine trihydrochloride was also included in the buffer. Once equilibrated, the solution was replaced with 800 µl of 1x restriction enzyme buffer containing 160 units of NotI or FseI (NEB) and incubated for 6 to 8 h at 4°C on a Nutator rotator to allow the enzyme to infiltrate the agarose plug. After 6 to 8 h, the tubes were transferred to a 37°C dry incubator and incubated for 12 to 16 h with agitation. Plugs were then washed twice in 2 ml TEX buffer at 4°C for 30 min each with agitation and equilibrated at room temperature in 2 ml of 0.5x Tris-borate-EDTA gel electrophoresis buffer. Restriction-digested high-molecular-weight DNA was then resolved by pulsed-field gel electrophoresis in 1% agarose-0.5x Tris-borate-EDTA for 20 h at 15°C (6 V/cm; 6- to 80-s switch time). DNA fragments were depurinated with 0.25 M HCl for 20 min and transferred by downward alkaline capillary transfer onto a Zeta-Probe GT membrane (Bio-Rad) with 0.4 N NaOH for 24 h. The membrane was neutralized with 0.2 M Tris-HCl (pH 7.5)-2x SSC (1x SSC is 0.15 M NaCl plus 0.015 M sodium citrate) for 10 min and baked for 30 min at 80°C. Probe used in hybridizations was generated from a 4.7-kb XbaI fragment containing the IRES-ß-geo cassette from pIBG-FTet (described above) and was labeled using Ready-to-Go labeling beads (GE/Amersham) and [
-32P]dCTP (GE/Amersham). Membranes were washed twice with sterile H2O and hybridized for 3 h with Rapid-hyb buffer (GE/Amersham) according to the manufacturer's instructions. Membranes were exposed to a Kodak phosphorimaging screen for 16 h and imaged using a Pharos FX imaging system and software (Bio-Rad).
X-Gal staining and histology. Whole-mount X-Gal staining was performed as described previously (13, 53). Specimens that were counterstained with alizarin red were first X-Gal stained and then processed for alizarin red staining as follows. Postfixed X-Gal-stained specimens were rinsed in H2O and then cleared in 1.0% KOH for 4 to 6 h at room temperature. Ossifying regions of bone were stained with 0.05 mg/ml alizarin red S (Sigma) in 1.0% KOH for 3 h to overnight. Specimens were then rinsed several times in H2O and cleared in glycerol. Whole-mount X-Gal-stained specimens were prepared for histology by dehydrating through a graded ethanol series, clearing in Citrisolv (Fisher), and embedding in paraffin. Samples were sectioned at 10 µm and collected on Superfrost-Plus slides (Fisher). Sections for histology were counterstained with eosin or nuclear fast red.
In situ hybridization. In situ hybridizations were performed with probes specific for mouse Bmp2 (52), Runx2, and Spp1. The Bmp2 probe template was kindly provided by B. L. M. Hogan. Mouse Runx2 and Spp1 probe templates were made by PCR amplifying fragments of the 3' untranslated region of either gene and ligating them into the pGEM T-easy vector (Promega). The 3' untranslated regions were inspected for nonrepetitive regions suitable for probe design by using the UCSC Genome Browser (35). Primers were as follows: for Runx2, AAAGCTTGCAGAACTCTTAGAATGA (forward) and ATCATATTAAAAAGCCAAGCACAAG (reverse); for Spp1, ACATGAAGAGCGGTGAGTCTAAG (forward) and AAATGCAGTGGCCGTTTG (reverse). Single-stranded antisense and sense RNA probes were labeled with digoxigenin by using either Sp6, T7, or T3 polymerase as described previously (29). Tissues were fixed in 4% paraformaldehyde-1x PBS at 4°C overnight, followed by dehydration through a graded methanol series. Specimens were then cleared in Citrisolv (Fisher) and embedded in paraffin. Samples were sectioned at 10 µm and collected on Superfrost-Plus slides (Fisher). Sections were deparaffinized in Citrisolv and rehydrated in 1x PBS. Tissues were permeabilized by incubation in 3 mg/ml proteinase K-1x PBS for 10 min. Probe hybridizations were performed overnight at 63°C in 40% formamide-5x SSC (pH 7.0) with a probe concentration of 100 ng/ml (Spp1), 150 ng/ml (Bmp2), or 200 ng/ml (Runx2). Posthybridization stringency washes were performed in 2x SSC and 0.2x SSC at 63°C for 1 h. Sections were then blocked with blocking reagent (Roche) and incubated in the presence of a 1:5,000 dilution of antidigoxigenin-alkaline phosphatase Fab fragments (Roche) for 1 h at room temperature. Sections were allowed to develop in an nitroblue tetrazolium-5-bromo-4-chloro-3-indolylphosphate solution containing 6% polyvinyl alcohol. Slides were washed and mounted using aqueous medium.
Immunohistochemistry.
Immunohistochemistry for Runx2 (also known as polyomavirus enhancer binding protein 2
[PEBP2
], Cbfa1, or AML1) was performed using a rabbit polyclonal antibody (PEBP2
A M-70; Santa Cruz Biotechnology, Inc.) (30) specific for Runx2/PEBP2
. Whole-mount X-Gal staining was performed on 10% neutral buffered formalin-fixed embryos. Whole-mount X-Gal-stained specimens were washed with PBS and postfixed in neutral buffered formalin for 3 h at 4°C. Specimens were dehydrated through a graded ethanol series. The material was then cleared in Citrisolv and embedded in paraffin. Samples were sectioned at 10 µm and collected on Superfrost-Plus slides (Fisher). Sections were deparaffinized in Citrisolv and rehydrated into H2O. Sections were placed in 10 mM sodium citrate, slowly heated until boiling, and then cooled to room temperature to aid in antigen retrieval. Endogenous peroxidase activity was quenched by placing the sections in 3% hydrogen peroxide-100% methanol for 30 min. Sections were blocked in a solution containing 3 mg/ml normal goat serum and 1% bovine sheep albumin (Jackson ImmunoResearch, Inc.) and then incubated in the presence of a 1:50 dilution of anti-PEBP2
A/Runx2 antibody for 30 min at room temperature. Secondary detection of rabbit immunoglobulin G was performed using a 1:200 dilution of biotinylated goat anti-rabbit immunoglobulin G (Vector Laboratories). Indirect enzyme-linked immunostaining was performed using a Vectastain ABC-diaminobenzidine substrate kit (Vector Laboratories). Sections were dehydrated and mounted using nonaqueous medium.
Microscopy and imaging. Sections from X-Gal staining, in situ hybridization, and immunohistochemistry experiments were analyzed and photographed on an Olympus BX51 microscope using bright-field optics and a digital camera.
| RESULTS |
|---|
|
|
|---|
|
In general, for each BAC the staining patterns at 15.5 dpc were mostly identical for multiple embryos/lines generated from the same construct. However, the overall lacZ expression intensity varied such that in some embryos staining was faint or undetectable in tissues where Bmp2 expression is relatively low. This variation could in theory be due to either copy number effects, fragmentation of the BAC causing loss of cis-regulatory elements, or repression due to site-specific position effects. Since BAC transgenes are thought to be fairly resistant to position effects and it has been reported that BAC-transgenic lines that have multiple copies are very likely to contain full-length copies of the transgene (24), we therefore estimated BAC transgene copy number by using quantitative real-time PCR. We confirmed our results using assays for both the ß-geo gene and the BAC vector (see Materials and Methods). This revealed that embryos or lines that had several (
3 or more) copies of the transgene almost always had robust expression patterns that were identical to those of the other high-copy founders, while weakly stained embryos almost invariably had only one or two transgene copies (Fig. 1 and data not shown). However, real-time PCR has limited ability to confirm whether a transgenic embryo or line contains a full-length BAC transgene insertion. Therefore, we performed Southern blot analysis on high-molecular-weight DNAs isolated from both 5' and 3' BAC-carrying stable lines used to generate the embryos in Fig. 1b to f and h to l, using rare restriction enzyme sites present in the BACs and a transgene-specific probe (see Materials and Methods; see also Fig. S1 in the supplemental material). This analysis revealed that both 5' and 3' BAC-carrying stable lines most likely contained one or more intact BAC copies, as evidenced by strong bands at the expected size of the BAC transgenes (see Fig. S1 in the supplemental material). These data are consistent with previous evidence that multicopy transgenes tend to integrate as tandem concatemers (8, 19). Our results are very consistent with the idea that BACs are highly resistant to position effects and that multiple-copy BAC transgenes usually contain one or more intact BAC copies that drive highly reproducible expression (24).
Since we confirmed that our 5' and 3' BAC lines each had multiple-copy integrations and that embryos from these lines had identical and complete staining patterns at 15.5 dpc as listed in Fig. 1, we analyzed staining in these lines at additional time points. Progeny embryos from representative stable lines for each Bmp2 lacZ-BAC were analyzed for lacZ expression from 9.5 dpc to 17.5 dpc. Whole-mount X-Gal-stained 9.5-, 11.5-, 13.5-, and 15.5-dpc embryo specimens for each Bmp2 lacZ-BAC are shown in Fig. 1.
The overall staining patterns for each BAC were quite distinct, suggesting differences in the distribution of cis-acting regulatory sequences surrounding the Bmp2 gene in the genomic region(s) covered by the two BACs. Similar patterns of lacZ expression were consistently observed in transgenic embryos across multiple independent transgene integration events at 15.5 dpc (Fig. 1a, e, g, k, and m). These consistent lacZ patterns closely reflect normal sites of endogenous Bmp2 expression as previously reported and/or reported here (see below), suggesting that ectopic or position site effects on BAC transgene expression were very minimal or nonexistent. While many lacZ expression patterns were specific to either BAC construct, some sites of lacZ expression were obviously similar for each BAC and reflect sites of endogenous Bmp2 expression. For example, at 11.5 dpc both BACs drive lacZ expression in the AER and in limb bud mesenchyme (28) (Fig. 1c and i, insets). In addition, both 5' BAC and 3' BAC transgenic embryos display lacZ expression in vibrissae (whisker hair follicles) at 15.5 dpc (Fig. 1e and k) in hair shaft precursor or precortical cells adjacent to the dermal papilla (see Fig. S2 in the supplemental material) (4, 43). Likewise, both Bmp2 lacZ-BACs consistently directed lacZ expression in the retinal pigmented epithelium (16), ventral footpads (28), gut mesentery, pituitary gland (71), and dorsal tongue papillae (33) at 15.5 dpc (Fig. 1a and data not shown; see Fig. S1 in the supplemental material). Additionally, expression was also seen in the hyaline cartilage of the nonossifying costal or ventral ribs in 5' (Fig. 1d to g) and 3' (Fig. 1j to m) BAC embryos. These data suggest that cis-regulatory elements sufficient to drive several aspects of embryonic Bmp2 expression are located in the region shared by both BACs (i.e., between kb 2.7 and +53.7), which also contains the Bmp2 transcription unit.
However, many lacZ expression domains were associated only with the 5' BAC or the 3' BAC but not both, probably because distant 5' or 3' regions are required to upregulate Bmp2 in specific embryonic locations. For example, embryos from stable 5' BAC lines displayed lacZ expression in heart mesoderm at 9.5 dpc (Fig. 1b), while transgenic embryos from a 3' BAC line had no observable X-Gal staining in the heart primordia at this stage (Fig. 1h). These data suggest that the genomic region present upstream of kb 2.7 contains sequences required for Bmp2 transcription in heart mesoderm.
Other patterns consistently observed in 5' BAC embryos but not in 3' BAC-transgenic mice include expression in the ventral trunk and extraembryonic tissues at 9.5 dpc (Fig. 1b) and in vasculature, dorsal tail bud, and eye lens and in a segmented pattern along the trunk at 11.5 dpc (Fig. 1c). Expression in eye lens and vasculature persisted at 13.5 dpc (Fig. 1d). By 15.5 dpc, staining in thymus, lung, and liver was clearly seen, as well as expression in all endochondral bones (Fig. 1e to g; see Fig. S2 in the supplemental material). In all 12 transgenic embryos/lines carrying the 5' BAC, expression in developing endochondral bones appeared to be concentrated to the hypertrophic chondrocyte zones (Fig. 1a and e to g). This was first noticeable in the proximal limb bones, vertebrae, and dorsal ribs at approximately 13.5 dpc (Fig. 1d). Closer inspection of 15.5-dpc whole-mount-stained 5' BAC embryos also consistently revealed staining in digit tips, choroid plexus, skeletal muscle, adrenal glands, synovial joints, and faintly in skin epithelium (Fig. 1 and data not shown; also see Fig. 2 and 3).
|
|
Additionally, 3' BAC embryos expressed lacZ in the placodes of the forming whisker follicles at 12.5 dpc and in pelage (coat) hair follicles clustered towards the dorsal midline at 13.5 dpc (not shown). Hair placode expression continued throughout the surface ectoderm at 15.5 dpc (Fig. 1k and m). Hair placode expression is also consistent with previously reported Bmp2 expression (51). Interestingly, since both 5' BAC and 3' BAC clones drive expression in whisker follicle hair shaft precortical cells at 15.5 dpc, this indicates that separate regulatory sequences control Bmp2 expression early in the epithelial hair placodes versus later in hair shaft precursor cells. As indicated for regulatory elements controlling Bmp2 expression in osteoblast progenitors, midbrain, and kidney, sequences controlling expression in hair follicle placodes are likely located between kb +53.7 and +207.1 3' to the promoter (Fig. 1a).
In summary, the analysis of 5' and 3' BACs indicated that many sites of embryonic Bmp2 expression require cis-regulatory elements that are more than 2.7 kb 5' or more than 53 kb 3' to the promoter. Regulatory elements for expression in hypertrophic chondrocytes, synovial joints, liver, lung, gut, and other tissues are in the 5' region, whereas elements controlling expression in osteoblast progenitors, tooth buds, hair follicles, kidney, midbrain, and other locations are in a distant 3' region; notably, only a minority of the distinct expression sites at 15.5 dpc were common to both BACs and are thus likely to be controlled by more proximal elements. Taken together, our data strongly suggest the presence of numerous enhancer elements in the 5' and 3' genomic regions unique to either Bmp2 lacZ-BAC, and some can affect promoter activity from large distances.
Bmp2 lacZ-BAC transgenes recapitulate endogenous expression in numerous tissues in vivo. In order to determine if the Bmp2 lacZ-BAC transgenes were able to direct expression patterns similarly to endogenous Bmp2 mRNA, we performed in situ hybridizations on wild-type embryo sections and compared them to corresponding X-Gal-stained transgenic specimens. Figure 2 shows a comparison of lacZ directed by either the 5' or 3' BAC clone to endogenous Bmp2 mRNA in several developing organs. Close similarity of 5' BAC expression to endogenous Bmp2 was seen in the gut, adrenal gland, and gonads at 15.5 dpc (Fig. 2a to i). Expression in gut epithelium was similar to previous reports of Bmp2 expression (4, 26) (Fig. 2a and b). Expression of Bmp2 in adrenal gland has not been previously reported. We found that the 5' BAC drives robust lacZ expression similarly to Bmp2 mRNA in the adrenal gland outer cortex layer (Fig. 2d and e). Similar patterns of lacZ and endogenous Bmp2 were also observed in the liver and in lung, where lacZ-positive cells were localized to proximal airway epithelium and mesothelium at 15.5 dpc (see Fig. S2 in the supplemental material). Bmp2 is also expressed in mouse gonads at midgestation (7) (Fig. 2g and h). Previous studies indicate similarity between our 5' BAC lacZ patterns and those for Bmp2 mRNA in digit tips, choroid plexus, and thymus (not shown) (22, 25, 28). We also detected weak 5' BAC expression and Bmp2 mRNA in skin epithelium adjacent to hair follicle placodes (not shown).
In 3' BAC embryos, we observed close similarities of lacZ expression to Bmp2 mRNA patterns in kidneys, mammary gland, and pelage hair follicle placodes (Fig. 2j to r). Expression in the kidney (Fig. 2j and k) was localized to comma- and s-shaped bodies (16). In the mammary gland, expression was prevalent in the developing secretory epithelium (Fig. 2m and n). Bmp2 and lacZ expression were prominent in the epithelial bud of hair follicle placodes at 15.5 dpc (Fig. 2p and q).
Sites of Bmp2 expression in numerous skeletal tissues were also recapitulated by BAC-transgenic mice (Fig. 3). The 5' lacZ-BAC drives lacZ expression in synovial joints, in a manner very similar to that for endogenous Bmp2, within the middle or deep articular layer of 17.5-dpc digit joints (Fig. 3a and b). The 5' BAC also directed expression in the hypertrophic cartilage of the vertebral bodies (Fig. 3d). In contrast, the 3' lacZ-BAC drives expression in the intervertebral discs of the spine and vertebral perichondrium and adjacent connective tissue (Fig. 3e and h). Both lacZ and Bmp2 mRNA are localized to fibrous cartilaginous cells of the annulus fibrosus surrounding the nucleus pulposus (Fig. 3e, f, h, and i). Bmp2 protein expression has also been reported to occur in the annuli of adult, senescent mice (70). In 15.5-dpc developing incisor buds, 3' BAC-directed expression in the enamel knot is very similar to that for endogenous Bmp2 mRNA (Fig. 3k and l). These patterns were also similar in molar tooth bud enamel knot (data not shown). Weak X-Gal staining and Bmp2 mRNA were also observed in the dental papilla (Fig. 3k and l). Our findings corroborate previous reports of Bmp2 mRNA expression in tooth buds (4, 51). Taken together, these results suggest that both BAC transgenes accurately reflect overlapping subsets of normal Bmp2 expression and that regulatory elements in 5' and 3' flanking regions within the Bmp2 locus are endogenous Bmp2 regulatory elements.
Distinct cis-regulatory regions drive Bmp2 expression in developing bone in vivo. The endochondral skeleton, which is composed of limb bones, the chondrocranium (cartilaginous bones of the skull base), and axial components (the vertebrae and ribs), undergoes ossification via cartilage templates that prefigure the future bone elements. In contrast, the intramembranous skeleton is composed of the cranial vault (calvarial bones) and viscerocranium (facial bones) and forms via the direct conversion of mesenchyme to bone, thus bypassing a cartilage template. Although these two ossification processes are distinct, differentiating osteoblasts are present in all endochondral and intramembranous bones. Embryonic Bmp2 expression has been previously documented in both the perichondrium and hypertrophic chondrocytes of endochondral long bones and also in intramembranous bones (10, 52, 60, 65). Given the significance of Bmp2 for osteoblast differentiation, we examined the skeletal expression of both BACs in more detail.
Figure 4 shows images of X-Gal-stained skeletons from 5' lacZ-BAC and 3' lacZ-BAC 17.5-dpc transgenic embryos generated from the stable lines described above, counterstained with alizarin red to highlight ossifying regions of bone. In 17.5-dpc 5' BAC embryos, strong staining is present in the endochondrally derived occipital bone (Fig. 4a). No X-Gal staining is present in the intramembranous parietal bone of the calvarium or zygomatic bone of the facial skeleton in 5' BAC embryos (Fig. 4a). In contrast, lacZ-positive cells are present in all intramembranous bones in 3' BAC embryos at 17.5 dpc (Fig. 4b), such as the parietal and zygomatic bones of the skull, as well as the occipital bone. In the radii and ulnae of 5' BAC embryos, stripes of X-Gal staining are localized to the ends of the ossifying bone collar (Fig. 4c). However, in 3' BAC embryos, X-Gal staining extends throughout the entire length of the bone collar (Fig. 4d). These data suggest that the 3' BAC can direct expression in cells in ossifying regions of both intramembranous and endochondral bones, while the 5' BAC directs expression only in endochondral bones.
|
|
To confirm whether these lacZ expression patterns matched endogenous Bmp2 expression at 17.5 dpc, we performed in situ hybridizations for Bmp2. In age-matched metacarpal sections of nontransgenic embryos, Bmp2 transcripts are present in both the hypertrophic chondrocytes and adjacent perichondrium/bone collar (Fig. 5g). These are very similar to the combined patterns of both lacZ expression domains observed in Fig. 5a and d. The apparent low levels of endogenous expression in comparison to X-Gal-stained sections may reflect difficulties in detecting low Bmp2 mRNA levels by in situ hybridization, whereas the enzymatic staining procedure used to detect lacZ is very robust. In the 17.5-dpc radius Bmp2 mRNA was detected in hypertrophic chondrocytes located in the growth plate (Fig. 5h and h'), similar to 5' BAC lacZ expression. Numerous Bmp2-positive cells are also distributed throughout the bone collar surrounding the marrow cavity (Fig. 5h and h') and in perichondrium flanking the hypertrophic cartilage zones (Fig. 5h'). Therefore, while both Bmp2 lacZ-BACs can confer endogenous patterns of Bmp2 expression in developing endochodral bones, they each direct complementary patterns that together reflect the complete Bmp2 regulatory domains at 17.5 dpc.
Since the 3' BAC could also direct expression in intramembranous bones, X-Gal-stained sections through the orbital plate of frontal bones from 15.5-dpc 3' BAC embryos (Fig. 5k) were compared to corresponding Bmp2 in situ hybridizations of wild-type embryos (Fig. 5l). This confirmed that lacZ-positive cells are distributed throughout the network of intramembranous bone in 3' BAC embryos (Fig. 5k), similar to endogenous Bmp2 mRNA (Fig. 5l) (10). We also observed lacZ expression throughout the intramembranous calvarium and mandible at 15.5 dpc (data not shown).
Clearly, each BAC directed a different subset of normal patterns in bone; the 5' BAC drives expression in hypertrophic chondrocytes of endochondral bones, and the 3' BAC drives expression in differentiating osteoblasts present in the perichondrium/bone collar and intramembranous bones. Surprisingly, this implies that regulatory sequences controlling embryonic Bmp2 expression in differentiating osteoblasts are located at greater than 53 kb from the Bmp2 promoter.
Relationship between Bmp2 transcription and osteogenic markers in vivo.
The conversion of the perichondrium to an osteoblast-containing bone collar is intimately linked to osteoblast differentiation, as well as molecular and morphological changes in adjacent chondrocytes, such as proliferation and subsequent hypertrophy (34). To correlate Bmp2 and 3' BAC transgene expression to the progression of osteoblast differentiation in vivo, we compared both to Runx2 and Spp1 expression. Runx2 (also named Cbfa1/PEBP2
A) is one of the earliest markers of, and is essential for, osteoblast differentiation (15, 42, 57). Runx2 is expressed in the mesenchyme surrounding future skeletal elements, in osteogenic cells within the perichondrium and bone collar, and in both prehypertrophic and hypertrophic chondrocytes at midgestation developmental stages (15, 38, 57). Spp1 (also named Bsp/bone sialoprotein/osteopontin) is expressed in osteoblasts and hypertrophic chondrocytes (20, 32). Its expression coincides with collagenous and noncollagenous matrix production, and in osteoblasts it is indicative of a bone matrix-producing phenotype (2). In addition, both Runx2 and Bmp2 are implicated in osteoblast-specific control of Spp1 transcription (15, 42, 46).
In 15.5-dpc radius condensations, Runx2 transcripts are present in cells within the perichondrium/bone collar flanking the hypertrophic and prehypertrophic zone (Fig. 6a and a '). In the osteogenic inner perichondrial layer, domains of Runx2 expression closely resemble those for endogenous Bmp2 mRNA (Fig. 6b and b') as well as 3' BAC-directed (Fig. 6c and c') lacZ expression at 15.5 dpc, as the latter both label perichondrial cells flanking the prehypertrophic and hypertrophic zones. However, the inner perichodrium region surrounding the nonhypertrophic proliferative and resting zones appears faintly positive for Runx2 (Fig. 6a, arrow) but not Bmp2-lacZ or endogenous Bmp2 transcripts. Runx2 is also transcribed in prehypertrophic chondrocytes (Fig. 6a and a'). In addition, endogenous Bmp2 mRNA is detected in the hypertrophic zone (Fig. 6b and b'). In contrast, Spp1 mRNA is undetected in perichondrial cells adjacent to the prehypertrophic chondrocytes (Fig. 6d and d') but is upregulated in bone matrix-producing osteoblasts flanking the hypertrophic zone and lining the nascent bone collar (Fig. 6d and d') and in late-stage hypertrophic chondrocytes. These data suggest that both Bmp2 and Runx2 are transcribed in osteoblast progenitors within the inner perichondrial layer flanking the prehypertrophic zone, prior to onset of Spp1 production by these cells. In addition, Bmp2 and Runx2 expression persists as these cells upregulate Spp1 and transition to a more mature osteoblastic phenotype capable of bone matrix production and bone collar formation (Fig. 6g).
|
Numerous cis-regulatory sequences are dispersed throughout the distant 3' region. In order to further refine the location of cis-acting sequences throughout the distant 3' region (kb +53.7 to +207.1), including that of the osteoblast progenitor element, four BACs were engineered to contain deletions in the 3' region and these were tested for ability to drive lacZ expression at midgestation time points. Using the galK homologous recombination system in bacteria (72), serial and nonoverlapping deletions of approximately 40 kilobases in size were engineered into the full-length Bmp2-lacZ 3' BAC (described above). We refer to these deletion BACs as D1-, D2-, D3-, and D4-BACs (see Fig. 8 for the locations of the deletions relative to the Bmp2 transcription start site [23]). This system allowed us to generate "seamless" contiguous deletions for systematically analyzing the cis-regulatory functions of the entire 153.5-kb 3' region (from kb +53.7 to +207.1). As described above, the 3' deletion lacZ-BAC DNAs were purified and injected into one-cell-stage mouse embryos. For each deletion BAC, between 3 and 10 transgenic founder embryos were analyzed directly for lacZ expression at 15.5 dpc. Each founder embryo represents an independent transgene integration event.
|
|
The 3' deletion BACs revealed other regions that are also essential for Bmp2-lacZ gene expression in various tissues (Fig. 8) and that were lacZ positive in full-length 3' lacZ-BAC embryos (Fig. 1a and k to m) at 15.5 dpc. For example, lacZ expression in the kidney s-shaped bodies is absent in D1-BAC embryos (Fig. 8; see Fig. S3 in the supplemental material). In addition, X-Gal staining was lacking in the tooth bud enamel knot (Fig. 8; see Fig. S3 in the supplemental material) and mammary glands (Fig. 8; and data not shown) in all D2-BAC embryos. Expression was also observed in the intervertebral discs (Fig. 3e and h) and midbrain regions of full-length 3' lacZ-BAC embryos (see Fig. S2 in the supplemental material) and in at least one or more each of the D1-, D2-, and D3-BAC embryos (Fig. 8 and data not shown); however, no expression was seen in the intervertebral discs in any of the 10 D4-BAC embryos analyzed (Fig. 8 and data not shown).
We analyzed lacZ transgene copy number for all of the 3' deletion BAC founder embryos (Fig. 7, insets). Figure 7 demonstrates the general correlation between strength of X-Gal staining and copy number, (e.g., hair follicle placodes). Notably, of the embryos shown in Fig. 7, X-Gal staining is relatively weaker in the D1-BAC embryo, which also has the fewest BAC copies (
4) yet still has detectable lacZ in bone and hair follicles. Comparison of the D3-BAC and D4-BAC embryos (both having an estimated copy number of
6) demonstrates that while staining in hair follicles and digits is almost identical for both, no X-Gal staining can be seen in the long bones of the D3-BAC embryo (Fig. 7d and i). Taken together, these data strongly suggest that numerous cis-regulatory sequences are dispersed throughout the distant 3' region and that a critical osteoblast progenitor-specific Bmp2 control element is located between kb +132.8 and +168.8.
Bmp2 expression in differentiating osteoblasts is controlled by a distant 3' enhancer. Because the deletion BAC analysis suggested that the region covered by D3 contains an osteoblast regulatory element, we next asked whether or not highly conserved sequences present within this region were capable of osteoblast progenitor-specific enhancer function. The entire D3 region covers approximately 36 kilobases of sequence (mouse February 2006 assembly; chr2, 133377031 to 133413037) (35). Based on inspection of cross-species conservation data on the UCSC genome browser (35), we focused on two regions that had the highest PhastCons conservation scores (64) in the D3 interval. These are located at approximately +156 and +160 kb from the previously reported Bmp2 transcription start site (23). Both had PhastCons scores of greater than 200, while no other regions in the D3 interval had scores over 89 (data not shown). These sequences were named ECR-1 and ECR-2 (Fig. 9a). ECR-1 is conserved across mammals and chicken, while ECR-2 conservation is detected only in mammals, based on MultiZ cross-species sequence alignments (6) (Fig. 9a); neither had alignment to frog (Xenopus tropicalis) or teleost fish (Tetraodon) (Fig. 9a). For both ECRs, the values for length of alignment versus percent identity for mouse/human and mouse/chicken comparisons were generated using the ECR browser (59) (Fig. 9a).
|
Multiple transgenic embryos were analyzed for each construct at 15.5 dpc. Some embryos had localized staining in the brain, vasculature, or other soft tissues which was not consistent between any two embryos; therefore, these were likely ectopic patterns due to position effects. Most embryos had staining in the eye, which is a previously known cryptic effect of the Hsp68-lacZ vector backbone (63). However, for both the ECR-1 + 2 and ECR1 constructs, a majority of transgenic embryos (10/12 and 7/8, respectively) had obvious lacZ expression in osteoblast progenitors present in developing bone at 15.5 dpc, while no ECR-2 embryos (n = 4) had any consistent lacZ expression. One ECR-1 + 2 and two representative ECR-1 embryos are shown in Fig. 9b to d. Specific X-Gal staining can be seen in the ossifying regions of the limb long bones and ribs (Fig. 9b to d). Oblique histological sections through the radius indicated numerous lacZ-positive osteoblast progenitors in the perichondrium flanking the prehypertrophic zone and in the bone collar (Fig. 9e to g). In addition, X-Gal-stained osteoblasts were present in the intramembranous mandible (Fig. 9h to j).
We noted that in some ECR transgenic embryos, X-Gal staining in osteoblast progenitors appeared mosaic (Fig. 9d, g, and j). These data might be partly explained by somatic mosaicism of the transgene in some founder embryos or position effect variegation; alternatively, we cannot rule out that other sequences are involved in activating uniform Bmp2 expression in differentiating osteoblasts. Nevertheless, these data strongly indicate that ECR-1 contains a conserved enhancer that activates Bmp2 expression in differentiating osteoblasts of both endochondral and intramembranous bones and that this element is located 156.3 kilobases from the Bmp2 promoter.
| DISCUSSION |
|---|
|
|
|---|
Each deletion specifically ablated expression in certain tissues, corroborating the full-length 3' BAC data and strongly suggesting that these endogenous Bmp2 expression patterns require modular cis-regulatory elements. By testing highly conserved sequences within the D3 region, we identified an element that can activate a heterologous promoter in osteoblast progenitors in vivo, suggesting that it normally functions as a tissue-specific enhancer. Thus, Bmp2 expression in osteoblast progenitors is controlled by a distant cis-regulatory module.
These data support the notion that Bmp2 cis-regulatory elements have evolved independently to permit new roles for Bmp2 in different structures and cell types. The large intergenic sequences flanking Bmp2 may have facilitated this process. In support of this, the large intergenic regions flanking Bmp2 show conserved synteny across vertebrates, including mammals and chicken (58). In addition, the osteoblast progenitor control sequence is highly conserved between these two clades; however, it appears to be absent in organisms that also have a bony skeleton, such as frog and teleost fish. This suggests that this sequence may have evolved in amniotes. Interestingly, endogenous Bmp2 expression in the perichondrium is similar in both mouse (data presented here) and chicken (