Previous Article | Next Article ![]()
Molecular and Cellular Biology, December 2008, p. 7394-7401, Vol. 28, No. 24
0270-7306/08/$08.00+0 doi:10.1128/MCB.01087-08
Copyright © 2008, American Society for Microbiology. All Rights Reserved.
,
Hematopoiesis Section, Genetics and Molecular Biology Branch, NHGRI, NIH, 49 Convent Drive, Bethesda, Maryland 20892-4442,1 Department of Biochemistry and Molecular Biology, The George Washington University Medical Center, 2300 Eye Street NW, Washington, DC 20037,2 Duke University Medical Center, DUMC Box 3912, Durham, North Carolina 27710,3 Duke Institute for Genome Sciences and Policy 2177B CIEMAS, 101 Science Drive, DUMC Box 3382, Durham, North Carolina 27708,4 Department of Pediatrics, Yale University School of Medicine, 333 Cedar Street, Box 208064, New Haven, Connecticut 06520-80645
Received 10 July 2008/ Returned for modification 3 September 2008/ Accepted 1 October 2008
|
|
|---|
|
|
|---|
Studies of global gene expression using RNA extracted from whole E13.5 wild-type and Eklf–/– fetal liver cells (16, 20, 31, 35) demonstrated that expression of the membrane structural protein mRNAs ankyrin, β-spectrin, dematin, and band 3 are significantly downregulated in EKLF-deficient fetal liver, as is the globin chaperone, alpha-hemoglobin stabilizing protein. The pleiotropic defects observed in Eklf–/– embryos led us to hypothesize that the defect in Eklf–/– embryos was the result of a primary failure of definitive erythropoiesis, as opposed to individual defects in the mature erythrocytes themselves. In this report, we demonstrate that E13.5 Eklf–/– erythroid progenitor cells fail to undergo terminal differentiation due to a failure of cell cycle progression. To identify potential EKLF target genes responsible for the block in erythropoiesis, we performed transcriptional profiling with RNA extracted from homogeneous populations of wild-type and Eklf–/– fetal liver early erythroid progenitor cells. Our analysis revealed significant dysregulation of >3,000 genes. Ingenuity Pathways Analysis (IPA; Ingenuity Systems, Mountain View, CA) of the transcriptional profiles identified significant perturbation of a cell cycle and DNA replication pathway and identified the transcription factor, E2F2, which regulates the G1-to-S-phase transition, as the most significantly affected gene in the pathway. We confirmed that Eklf–/– early erythroid progenitor cells were perturbed at the G1-to-S-phase transition. Chromatin immunoprecipitation (ChIP) analysis demonstrated EKLF occupancy at the proximal E2f2 promoter in vivo. EKLF was required to confer DNase I sensitivity across the E2f2 locus. These results demonstrate that E2f2 is a direct target of EKLF and that EKLF deficiency leads to cell cycle perturbation and defective terminal erythroid differentiation.
|
|
|---|
Flow cytometric analyses. Fetal livers from E13.5 wild-type and Eklf–/– embryos were dissociated to a single cell suspension and stained with 7-amino-actinomycin D and anti-CD71 and Ter119 antibodies (BD Pharmingen, San Jose, CA). Cell populations were isolated by using a FACSAria flow cytometer running FACSDiva software (BD Biosciences, San Jose, CA). Cell cycle analysis was performed by using propidium iodide staining solution (BD Pharmingen) on a FACSCalibur flow cytometer (BD Biosciences) running FlowJo software (TreeStar, Inc., Ashland, OR). Each experiment was performed a minimum of four times, and the data were analyzed by using a Student t test.
Colony-forming assays. Sorted R1-R5 cells from E13.5 wild-type fetal livers and R1 and R2 cells from E13.5 Eklf–/– fetal livers were suspended at 106 or 104 cells/ml and plated in duplicate in either 1.5 ml of Methocult M3334 (CFU-E) or 1.5 ml of Methocult GF M3434 (BFU-E and CFU-GM; Stem Cell Technologies, Vancouver, British Columbia, Canada). CFU-E colonies containing 12 to 20 cells were counted 2 to 3 days postplating (wild type) or 3 to 5 days postplating (Eklf–/–). BFU-E colonies containing 10,000 to 20,000 cells were counted 10 to 12 days postplating (wild type) and 13 to 15 days postplating (Eklf–/–). CFU-GM colonies were counted from the same plates 10 to 12 days postplating for both the wild-type and Eklf–/– cell populations. Each experiment was performed a minimum of four times, and the data were analyzed by using a Student t test.
Microarray analyses. RNA was prepared from R1 and R2 cells sorted from E13.5 wild-type and Eklf–/– fetal liver cells using TRIzol reagent (Invitrogen, Carlsbad, CA) according to the manufacturer's instructions. RNA quality was verified with an Agilent 2100 bioanalyzer (Agilent Technologies, Santa Clara, CA). GeneChip (Affymetrix, Santa Clara, CA) target probes from three independent Eklf–/– and three wild-type RNA samples were prepared according to Affymetrix instructions starting with 10 µg of RNA. Double-stranded cDNA was synthesized by using a T7-linked oligo(dT) primer, followed by second-strand synthesis. Biotin-labeled cRNA, produced by in vitro transcription, was synthesized and subsequently fragmented. The fragmented cRNA was hybridized to Affymetrix GeneChip Mouse Genome 430 2.0 arrays at 45°C for 16 h, followed by washing and staining with streptavidin-phycoerythrin (Molecular Probes, Carlsbad, CA). Signal amplification was performed with a biotinylated anti-streptavidin antibody (Vector Laboratories, Burlingame, CA). Arrays were scanned with an Affymetrix GeneChip Scanner, and hybridization patterns were detected as light emitted from the fluorescent reporter groups incorporated into the target and hybridized to oligonucleotide probes. Signal intensity measurements computed in the Affymetrix Microarray Analysis Suite 5.0 served as a relative indicator of the level of expression. Scaling factors were also computed for each array based on an arbitrary target intensity of 500. All raw data files (CEL) are available (http://data.genome.duke.edu/EKLFDef) and have been submitted to the GEO database. These experiments comply with MIAME (Minimum Information About a Microarray Experiment) standards.
Genespring 6.1 (Agilent) was used to perform initial data analysis. Affymetrix MAS 5.0 signal intensities were normalized per chip and per gene by taking each measurement divided by the 50th percentile of all measurements in that sample, and each gene was divided by the median of its measurements in all samples. Based on the data from three independent analyses of each condition and a threshold of twofold change in expression relative to the control, a two-way analysis of variance with a P value cutoff of 0.05 (applying a Bonferroni correction) was performed. The expression of each gene was reported as the ratio of the value obtained for each condition relative to control conditions after data normalization. Affymetrix probe identifications were imported into the IPA software and categorized based on location; cellular components; and the reported or suggested biochemical, biologic, and molecular functions. Genes were also mapped to networks in the Ingenuity database and ranked by score. The score is the probability that a collection of genes equal to or greater than the number in a network could be achieved by chance alone.
Expression analyses.
Quantitative real-time PCR was performed to confirm the levels of mRNA observed in the transcriptional profiling using sequence-specific oligonucleotide primers (see Table S1 in the supplemental material) designed to amplify
150-bp fragments, each spanning an intron. PCRs were performed with Sybr green in an iCycler (Bio-Rad, Hercules, CA) to measure the fluorescence intensity. The values obtained for target gene expression were normalized to beta-actin and were expressed relative to the expression in control samples. For the threshold cycle (CT) calculations, the 2–
CT formula was used: 
CT = [CT(target) – CT(actin)]experimental sample - [CT(target) – CT(actin)]control sample.
Western blot analysis was performed with wild-type and Eklf–/– E13.5 total fetal liver proteins as described previously (35). Blots were probed with a rat anti-mouse E2F2 polyclonal antibody (sc-633) and a goat anti-human actin polyclonal antibody (sc-1616), both from Santa Cruz Biotechnology (Santa Cruz, CA). Western blots were exposed to an imaging plate, scanned using a Storm 860 scanner, and quantified by using scanning densitometry with ImageQuant software (both from Molecular Dynamics, Inc., Sunnyvale, CA).
ChIP analysis. We performed ChIP on chromatin obtained from sorted R1+R2 fetal liver cells from E13.5 HA-EKLF-TAP-tagged mice (43). Cross-linked cells were processed by using the Magna ChIP A assay kit (catalog no. 17-610; Millipore, Billerica, MA) according to the manufacturer's instructions. Chromatin was immunoprecipitated with a mouse monoclonal anti-HA probe (F-7, sc-7329X; Santa Cruz). Nonimmune control samples were incubated with protein A beads but without antibody. The PCR primers are listed in Table S2 in the supplemental material. We performed the analyses on five independent chromatin preparations from sorted cells. Each sample was analyzed in triplicate.
DNase I sensitivity assay.
A DNase I sensitivity assay of the E2f2 locus was performed by using the PCR-based method described by Bernat et al. (5) and Dorschner et al. (14). Chromatin was prepared from
107 nuclei from sorted R1+R2 fetal liver cells obtained from E13.5 wild-type or Eklf–/– embryos as previously described (39). Nuclei were exposed to increasing amounts of DNase I (0.0 to 4.0 g/ml; Worthington Biochemical Corp., Lakewood, NJ) and incubated at 37°C for 10 min before DNA extraction. Southern blot analysis with a probe for a relatively DNase I resistant region of the keratin locus was performed to identify samples with equivalent levels of DNase I digestion. DNA from untreated, "low DNase" (the first sample with evidence of DNase I activity) and "high DNase" (the first sample with an absence of high-molecular-weight DNA) was analyzed by quantitative real-time PCR using a set of 287 primer pairs (average amplicon size of 280 bp, average gap of 150 bp) covering the E2f2 locus and 50 kb of sequence up- and downstream of the gene. The DNA was amplified on a 384-well format ABI Prism 7900 Fast Real-Time PCR instrument (Applied Biosystems, Foster City, CA). The difference in CT between DNase-treated and untreated template for each primer pair was calculated as
CT = CT(treated) – CT(zero), which serves as a measure of the relative DNase sensitivity of the amplified region. Three or four independent chromatin preparations were analyzed for each primer pair in triplicate. The strict statistical criteria described by Bernat et al. (5) and Dorschner et al. (14) were used to compare the DNase I sensitivity of regions. The mean CT and standard deviation for each point in the wild-type and Eklf–/– samples were calculated and analyzed by a Student t test. A DNase I-sensitive region is defined as a region with significantly more digestion (P < 0.01).
|
|
|---|
![]() View larger version (46K): [in a new window] |
FIG. 1. (A) Fluorescence-activated cell sorting analysis of E13.5 fetal liver cells from wild-type and Eklf–/– embryos. The cells were stained with anti-CD71 (transferrin receptor; y axis) and TER119 (x axis) antibodies. The sorting windows for the R1-R5 populations are indicated by the labeled ovals. (B) Morphology of representative Wright-Giemsa-stained cells in the sorted R1-R5 populations. Wild-type R1 was composed of 98% large blast cells and 2% other cells, wild-type R2 was composed of 99% smaller blast cells and 1% other cells, wild-type R3 was composed of 99% smaller hemoglobin-staining cells with compacted nuclei and 1% other cells, wild-type R4 was composed of 99% cells with dark hemoglobin staining and decentralized nuclei and 1% other cells, and wild-type R5 was composed of 99% reticulocytes and 1% other cells. Eklf–/– R1 was composed of 97.5% large blast cells and 2.5% other cells, while Eklf–/– R2 was composed of 97% smaller blast cells and 3% other cells. WT, wild type.
|
![]() View larger version (32K): [in a new window] |
FIG. 2. Analysis of erythroid colony forming cells in E13.5 fetal liver cells from wild-type and Eklf–/– embryos. Fetal liver cells were sorted using the windows depicted in Fig. 1. The R1 and R2 populations were isolated and suspended in semisolid medium. The black bars represent the average number of colonies from wild-type cells (with the standard deviation), and the open bars represent the average number of colonies from Eklf–/– cells (with the standard deviation). (A) Frequency and total number of BFU-E colonies. BFU-E colonies were counted on days 10 to 12 (wild type) and days 13 to 15 (Eklf–/–) of culture. (B) Frequency and total number of CFU-E colonies. CFU-E colonies were counted on day 2 to 3 (WT) and day 3 to 5 (Eklf–/–) of culture. The P values for the significant differences between wild type and Eklf–/– are shown. NS, not significant. There was no significant difference in the frequency or absolute number of CFU-GM in wild-type and Eklf–/– fetal liver cells. The total numbers of colonies were calculated as follows: (number of colonies/10,000 R1 or R2 cells) x (number of R1 or R2 cells/total fetal liver cells) x the total number fetal liver cells. WT, wild type.
|
For an unbiased analysis of our transcript profile, we performed IPA on the more inclusive set of 3,199 genes showing dysregulated expression at the P < 0.05 level. The Affymetrix probe lists were overlaid on a cellular pathway map based on location; cellular components; and reported or suggested biochemical, biologic, and molecular functions in the Ingenuity Systems (Redwood City, CA) Pathways Knowledge Base (winter 2004 release containing 20,000 genes). The genes with significantly altered expression were mapped to genetic networks and then ranked by score. The score is the probability that the collection of probe genes equal to or greater than the number represented in a network could be achieved by chance alone. A score of
3 indicates a 1/1,000 chance that the focus genes are in a network due to random chance. The pathway with the highest IPA score was cell cycle and DNA replication, which generated a network involving 35 focus genes with significantly dysregulated expression (score = 27; Fig. 3). At a node of this network was E2F2, a transcription factor and a critical regulator of cell proliferation by its control of genes regulating S-phase entry and DNA synthesis (1, 30), whose expression was significantly decreased in Eklf–/– R1+R2 fetal liver cells (P < 0.037). In four paired, independent RNA samples, RT-PCR analysis verified that the level of E2f2 mRNA in Eklf–/– R1+R2 cells was decreased to 4.8% ± 3.5% of the level seen in wild-type R1+R2 cells (P < 0.001) (Fig. 4A). Consistent with the mRNA levels, Western blot analysis demonstrated that E2F2 protein was undetectable in Eklf–/– R1+R2 fetal liver cells (Fig. 4B; the band immediately beneath the E2F2 band is a nonspecific band recognized by the SC-633 antibody as described in the original reference) (29). PhosphorImager analysis of these data confirmed a fivefold reduction in E2F2 in Eklf–/– R1+R2 fetal liver cells (Fig. 4C). An independent analysis of the microarray data using GeneSifter software (VizX Labs, Seattle, WA) identified the same pathway and target gene.
![]() View larger version (26K): [in a new window] |
FIG. 3. Network of differentially expressed genes in the cell cycle control pathway in sorted R1/R2 Eklf–/– E13.5 fetal liver cells. IPA was performed on microarray data comparing R1 and R2 mRNA isolated from wild-type and Eklf–/– E13.5 fetal liver cells. Using 3,199 differentially expressed genes as the focus gene set, the highest-scoring pathway was cell cycle and DNA replication (score = 27; scores >3 indicate a significant chance that a biologically relevant network has been identified). Red nodes indicate upregulation in Eklf–/– cells; green nodes indicate downregulation in Eklf–/– cells. Arrows and lines denote interactions between specific genes within the network. A, activation; E, expression regulation; P, phosphorylation; PP, protein-protein interaction; PD, protein-DNA interaction; PR, protein-RNA interaction; RB, regulation of binding; L, proteolysis.
|
![]() View larger version (42K): [in a new window] |
FIG. 4. Verification of decreased levels of E2f2 in sorted R1/R2 cells from E13.5 wild-type and Eklf–/– fetal livers. (A) Quantitative RT-PCR analysis of E2f2 mRNA in sorted R1/R2 cells from E13.5 wild-type and Eklf–/– (–/–) embryos. The mean mRNA levels (with the standard deviation) are related to the level in wild-type cells (designated 100%). (B) Representative Western blot of E2F2 protein levels in sorted R1/R2 cells from E13.5 wild-type, heterozygous Eklf+/– (+/–) and Eklf–/– (–/–) embryos. Actin levels are shown as controls. (C) Relative E2F2 protein levels (with the standard deviation) in sorted R1/R2 cells from E13.5 wild-type, heterozygous Eklf+/– (+/–), and Eklf–/– (–/–) embryos. The mean protein levels of each genotype were quantified by using a Molecular Dynamics Storm 860 scanner and are normalized to the level in wild-type cells (designated 100%). WT, wild type.
|
![]() View larger version (21K): [in a new window] |
FIG. 5. Analysis of cell cycle in sorted R1 and R2 E13.5 fetal liver cells from wild-type and Eklf–/– embryos. Fetal liver cells were sorted into R1 and R2 populations using the windows depicted in Fig. 1 and analyzed for DNA content with propidium iodide. The bars represent the average percentage of cells (with the standard deviation; y axis) at each stage of the cell cycle (x axis). Black bars represent cells sorted from wild-type fetal liver, and the open bars represent cells sorted from Eklf–/– fetal liver. (Left panel) Cell cycle analysis of R1 cells; (right panel) cell cycle analysis of R2 cells. WT, wild type.
|
-globin and keratin promoters (data not shown). Regions 1 to 3 contain three (region 1), two overlapping (region 2), or two distinct (region 3) NCNCNCCCN motifs. These regions were significantly more enriched than the -C region (P < 0.002), indicating EKLF occupancy. The two most proximal regions (region 4, containing one NCNCNCCCN motif, and region 5, containing two overlapping NCNCNCCCN motifs) were significantly enriched compared to both the -C region and regions 1 to 3 (all P < 0.008), demonstrating higher levels of EKLF occupancy (Fig. 6).
![]() View larger version (10K): [in a new window] |
FIG. 6. ChIP analysis of the E2f2 promoter region in sorted R1/R2 E13.5 fetal liver cells from wild-type embryos expressing HA-EKLF. Chromatin fragments were precipitated with anti-HA antibody and amplified with six sets of primers from the E2f2 promoter region. Regions 1 to 5 contain one or more consensus EKLF binding motifs (NCNCNCCCN). The negative control region (-C) does not contain an EKLF binding motif. The location of the primers flanking each region is indicated. The dark gray bars represent the enrichment of each sequence relative to the -C region (gray bar; designated as 1.0).
|
![]() View larger version (36K): [in a new window] |
FIG. 7. DNase I sensitivity analysis of the mouse E2f2 locus in sorted R1/R2 E13.5 fetal liver cells from wild-type and Eklf–/– embryos. DNA was extracted from the nuclei of sorted cells treated with or without DNase I. The DNA was amplified with a series of 287 primers spanning the E2f2 gene plus 50 kb of 5' and 3' flanking DNA by using a real-time PCR assay. The top of the figure shows the location of the E2f2 and the flanking loci on mouse chromosome 4. The paired vertical bars represent the DNase I sensitivity at each point in this region. DNase I sensitivity was calculated by subtracting the average number of cycles needed to amplify a specific amount of product in untreated DNA from the average number of cycles needed to amplify the same amount of product in the DNase-treated DNA. The results are the means of triplicate analyses of each point in a composite of three or four (depending on the region) independent experiments. The blue bars represent the signals in DNA extracted from wild-type chromatin. The black bars represent the signals in DNA extracted from Eklf–/– chromatin. The mean level of DNase digestion in the region indicated by bracket A is similar in wild-type and Eklf–/– chromatin. The mean level of DNase digestion in the region indicated by bracket B is significantly greater in wild-type chromatin compared to Eklf–/– chromatin (P < 0.001). WT, wild type.
|
|
|
|---|
2-fold differences in expression between unfractionated wild-type and Eklf–/– fetal liver mRNA were expressed at
4-fold-higher levels in wild-type R3+R4+R5 cells compared to wild-type R1+R2 cells. Our observations support a model of defective definitive erythropoiesis in Eklf–/– embryos, resulting in the limited production of abnormal erythrocytes at insufficient numbers to support fetal development. With quantitatively more BFU-E and CFU-E than wild type, the severe anemia of E13.5 Eklf–/– embryos is not due to a deficit of erythroid progenitor cells. Our unbiased analysis of the entire gene expression profile suggested that cell cycle defects were the major contributor to the failure of definitive erythropoiesis and that the block in cell cycle progression was in part due to a lack of E2F2. During early G1, an unphosphorylated retinoblastoma tumor suppressor protein (pRb) sequesters E2F proteins in the cytosol, preventing E2F from transactivating target genes, such as the G1 cyclins A and E. Cyclins A and E associate with cyclin-dependent kinases to phosphorylate pRb, causing a release of E2F. Free E2F then transactivates target gene transcription, stimulating a feed-forward mechanism of cyclin-cyclin-dependent kinase activation and pRb phosphorylation, pushing the cell cycle into S phase (17). Dirlam et al. (12) have demonstrated that E2F2 is the major pRb-associated member of the E2F family active during terminal erythroid maturation.
Our model predicts that EKLF acts upstream of E2f2. Consistent with this model, we also observed significant decreases in the mRNA level of Rbl1, Rbl2, E2f4, E2f1, and the E2F target genes cyclin A and cyclin E in Eklf–/– R1 and R2 fetal liver cells. Additional evidence is provided by the observation that E2f1–/– E2f2–/– mice have impaired erythropoiesis due to delayed S-phase progression coupled with a severe decrease in the level of β-globin mRNA (30), similar to Eklf–/– embryos. Rb–/– mice also exhibit defective definitive erythropoiesis, marked by impaired cell cycle exit and failure to enucleate (11, 23, 36), whereas E2f4–/– mice have an abnormal expansion of erythroid progenitors and abnormal cell cycle progression (21, 26). Although E2F2 is expressed in many cell types, its role(s) in erythropoiesis appears to be especially critical. We propose that EKLF binding is responsible for maintaining the E2f2 locus in an active chromatin state for transcription in early erythroid cells. In turn, E2F2 protein transactivates other pathway members, allowing progression of early erythroid progenitor cells, including BFU-E and/or CFU-E, through the G1/S transition in the divisions prior to terminal erythroid differentiation.
Published ahead of print on 13 October 2008. ![]()
Supplemental material for this article may be found at http://mcb.asm.org/. ![]()
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»