NRF2, a member of the NFE2 family of transcription factors, is not essential for murine erythropoiesis, growth, and development (original) (raw)

Abstract

The locus control region of the β-globin gene is composed of four erythroid-specific hypersensitive sites. Hypersensitive site 2 has been shown to be a powerful enhancer and contains a tandem repeat sequence for the transcription factors AP1 and NFE2 (activating protein 1 and nuclear factor erythroid 2, respectively). The human NRF2 (NFE2 related factor 2) has been isolated by bacterial expression screening using this core sequence as a probe. p45-NFE2, NRF1, and NRF2 belong to the CNC (“cap ‘n’ collar”) subfamily of the basic region–leucine zipper transcription factors, which exhibits strong homology at specific regions such as the “CNC” and the DNA binding and leucine zipper domains. Although the erythroid-specific p45-NFE2 has been implicated in globin gene regulation, p45-NFE2 null mice succumb to bleedings due to lack of platelets and those that survive exhibit only a mild anemia. To determine the function of NRF2, which we found to be widely expressed in vivo, we have characterized the genomic structure of the mouse NRF2 gene, disrupted the Nrf2 gene by homologous recombination in mouse embryonic stem cells (ES cells), and generated NRF2−/− mice. Homozygous mutant mice developed normally, were not anemic, reached adulthood, and reproduced. Our studies indicate that NRF2 is dispensable for mouse development.

Keywords: locus control region, hypersensitive site, gene knockout, CNC-bZip protein

A wealth of knowledge has been accumulated by the study of the human β-globin gene expression generated since the characterization of the locus control region (LCR) (1, 2, 3), which consists of four erythroid specific DNase I hypersensitive sites (HS1–HS4) and is responsible for the erythroid and developmental stage specific expression of the globin genes (4, 5, 6, 7, 8, 9). Work centered on hypersensitive site 2 (HS2) carried out through in vitro and_in vivo_ footprinting experiments (10, 11) identified a tandem AP1-NFE2 (activating protein 1 and nuclear factor erythroid 2, respectively) repeat sequence (GCTGAGTCATGATGAGTCA). Experiments using transient transfections and transgenic mice (10, 12, 13, 14) suggested that this dimeric motif was essential for strong enhancer activity and that the proteins that bind to this sequence play a critical role in the expression of β-globin and β-globin-like genes in humans.

Work undertaken by several groups to clone the transcription factors that bind to this NFE2 site so far identified three genes that are members of the novel CNC (“cap ‘n’ collar”) subfamily of basic region–leucine zipper (bZIP) transcription factors. CNC (15) is a homeotic gene involved in the development of the head and neck structure in Drosophila. NFE2 is a heterodimer consisting of the hematopoietic-specific p45-NFE2 subunit (16) and a ubiquitously expressed p18-MAF subunit (17). Two p45-NFE2-related genes_NRF1_ (or LCR-F1) and NRF2 (18, 19, 20) have also been cloned in our laboratory and others. These two genes are not restricted to the erythroid cells; they are widely expressed in different tissues, though with different levels of expression. p45-NFE2, NRF1, and NRF2 constitute a new subfamily transcription factors, for they share homology at several protein domains including the CNC and bZIP. Most notably their DNA binding domains are virtually identical in amino acid sequence with only conservative amino acid substitutions. Although these three genes are mapped to three different chromosomes (21), each is situated within loci that contain a member of gene families such as integrin, collagen, and Hox genes. NRF1 and NRF2 have been shown to be capable of transactivating reporter genes in cell culture (10, 20, 22), but their function in vivo are still unknown.

Targeted disruption of p45-NFE2 (23, 24) results in blockage of maturation of megakaryocytes, but only mild impairment of erythropoieis. Mice lacking p18-MAF show no detectable phenotype and are indistinguishable from their normal littermates (25). Although NRF1 and NRF2 are not up-regulated in the _p45-NFE2_null mice (24), they may complement the lack of p45-NFE2. To investigate the function of NRF2 in vivo, we used homologous recombination in embryonic stem cells (26) to generate NRF2 mutant mice.

MATERIALS AND METHODS

cDNA Library Screening.

The human NRF2 cDNA was used to screen a mouse lung cDNA library in λZAPII (Stratagene). Several clones were purified, excised with helper phage, and recircularized to generate subclones in Bluescript SK− phagemid vectors (Stratagene). To determine the sequence of the 5′ end of cDNA, rapid amplification of cDNA 5′ ends (RACE) was carried out according to the manufacturer’s instructions (GIBCO/BRL). poly(A) RNA was isolated from mouse kidney using the FastTrack 2.0 kit (Invitrogen) and used as a template for the synthesis of the first strand cDNA.

Genomic Clones Isolation.

Several phage clones were isolated from a 129SVJ library in λFIX II (Stratagene) through several rounds of screening, first with a mouse mNRF2 cDNA probe and then with oligonucleotides. The whole inserts from the phage clones were subcloned into pBluescript KSII at the _Not_I site to facilitate their mapping.

Sequence Analysis.

Sequence analysis was performed using the chain termination methods with cycle sequencing kits (Promega and Amersham) and the sequence analysis package from Genetics Computer Group (Madison, WI).

In Situ Hybridization.

Tissues from 15.5-day-old embryos of normal or knockout mice were harvested and fixed immediately in PBS containing 4% paraformaldehyde. After fixation for 24 hr at 4°C, the specimens were washed in PBS, dehydrated through a series of ethanol washes, cleared in xylene, and embedded in paraffin. Sections were cut at 7 microns and mounted on Superfrost/plus (Fisher Scientific) treated slides. Riboprobes were prepared from an _Eco_RV–_Hin_dIII fragment of the cDNA that is in the region deleted by the targeting vector (see below). This fragment is subcloned in pBluescript KSII. [35S]UTP-labeled riboprobes were synthesized with T3 or T7 RNA polymerase from linearized plasmid DNA. Hybridization and washes were carried out according to established protocols (27). The slides were then dipped in NTB2 emulsion (Kodak), and the emulsion was exposed for 4–7 days at 4°C. After development, the slides were stained with eosin and hematoxylin for histology and then mounted with Permount. Microphotography was performed using a Zeiss Axiophot microscope.

Targeting Vector Construction.

The targeting vector was based on the pPNT vector (28), which contains both the positive selection neomycin resistance and the negative selection thymidine kinase markers under the control of the phosphoglycerate kinase promoter and polyadenylylation signal. The 5′ homology region consists of a 5.5-kb _Nhe_I–_Eco_RV fragment that encompasses the 3′ end of intron 1 to part of exon 4. This 5′ fragment was first fused in-frame to a LacZ reporter cassette containing the last intron and polyadenylylation signal of protamine 1 (29) and is subcloned into the _Not_I–_Xho_I sites of pPNT. The 3′ homology region consists of a 6.2-kb _Xba_I–_Nhe_I fragment downstream of the NRF2 gene. The construct is linearized with_Not_I before electroporation.

ES Cell Line and Electroporation.

The JM1 ES cells derived from 129SVJ mice were cultured on embryonic fibroblast feeder layers in DMEM supplemented with 15% FBS/0.1 mM of 2-mercaptoethanol/2 mM glutamine/0.1 mM nonessential amino acid/purified leukemia inhibiting factor (GIBCO/BRL) at 103 units/ml. Cells (n = 107) were electroporated with 20 μg of linearized targeting construct DNA using an electroporator (Bio-Rad) in 1 ml of complete ES cell culture medium at 240 V and 25 mF. Transfected cells were first cultured in nonselective medium overnight, and then subjected to selection with G418 at 200 μg/ml (GIBCO) and ganciclovir at 2 μM (Syntex, Palo Alto, CA). Single colonies were picked after 10–12 days, and cells were expanded in 24-well dishes. Aliquots of each clone were frozen for stock and the rest were expanded for DNA analysis. The ES cells with homologous recombination were injected into blastocysts, and chimeric animals were then bred to produce F1 mice. Southern blots were used to identify germ-line transmission, and the heterozygous animals were interbred to produce homozygous knockout animals.

Northern and Southern Blot Analyses.

Northern and Southern blots analyses were performed according to standard procedures (30). Total RNA was extracted from intestine, kidney, stomach, and testis using RNAzol B (Biotex Laboratories, Edmonton, Canada), and poly(A) RNA was subsequently isolated. Multiple tissue RNA blots and mouse embryo blots were purchased from CLONTECH.

Hematologic Analysis.

Blood was collected from adult mice by retroorbital puncture. The hematocrit, hemoglobin concentration, mean corpuscular volume, mean corpuscular hemoglobin, reticulocyte count were done (CVD, Sacramento, CA).

Reverse Transcription–PCR Analysis.

Total RNA was used as template for the first strand cDNA synthesis using Superscript II reverse transcriptase (GIBCO/BRL) and random hexamers. Primer pairs spanning at least one intron were used for PCR analysis to distinguish amplification products of cDNA from those of genomic DNA.

RESULTS

mNRF2 cDNA.

Sequence analysis of 13 clones isolated from a mouse lung cDNA library revealed that they were all cDNA fragments derived from the same messenger RNA. Two of the overlapping clones covered 2005 bases, excluding the poly(A) tail. RACE using gene specific primers on exon 2 extended the final clone to 2380 bases. After we obtained the mouse sequence, the mNRF2 sequence was reported by Chui et al. (31). Our sequence data is similar to that published one with several differences. (i) The translation start site in our sequence is extended by 16 amino acids at the amino end. We believe that these 16 amino acids are translated as they are present in a human cDNA sequence of Nrf2 (GenBank accession no. T39964T39964) (Fig. 1A). The difference between our translation sequence and that of Chui et al. (31) could be due to a G nucleotide missing in their sequence. This G nucleotide is present in both our cDNA and genomic sequences. (ii) The 5′ untranslated region of our cDNA, determined from the analysis of six RACE clones, is shorter than that of Chui et al. (31), whose sequence extends for approximately 105 nucleotides further. This difference is unexplained. This extension is present in our genomic sequence with a few mismatches (Fig. 1B).

Figure 1.

Figure 1

Sequence of NRF2. (A) Comparison of the putative translation start site of NRF2. The translation of the murine cDNA at the exon 1 and exon 2 junction is compared with the translation of the human NRF2 sequence T39964T39964 found in the EST data base. The extra amino acids are italicized, the methionine described in previous publications is in boldface type, and the G nucleotide underlined is missing in entry U20532 (GenBank) (31). (B) Alignment of the 5′ ends of the cDNAs (Nrf2cDNA and U20532) and genomic (Nrf2GEN) sequences of the first exon of mNrf2. The numbers above the cDNA represent the number of RACE clones that start at that base. ∗, Putative translation start sites derived from our sequence; •, the published translation start site (31). The sequence of intron 1 is between parentheses with gt representing the splice donor junction; sequence differences are underlined; the missing nucleotide G in U20532 is in boldface italicized type in our cDNA and genomic sequences.

Genomic Structure of m_NRF2._

The transcriptional unit of mNRF2 was mapped within 30 kb segment of genomic DNA. The NRF2 gene consists of five exons interrupted by four introns, the first of which spans ≈25 kb. The sequences at all the exon–intron boundaries were determined and are shown in Table1. The DNA sequence was generated by a combination of subcloning and primer walking. The promoter of_NRF2_ is very GC rich and contains neither a TATA box nor a CCAAT box. Several SP1, EKLF, and AP-2 binding sites are found within the first 500 bp upstream of exon 1. Similar to the _p45-NFE2_locus (32), a gene coding for an heterogenous nuclear ribonucleoprotein (hnRNP) is found 10 kb downstream of NRF2_with its transcriptional direction opposite that of_NRF2. This is further evidence suggesting that the CNC-bZIP family of gene evolved from one ancestral gene.

Table 1.

Sequences at the exon–intron boundaries of murine Nrf2

Exon no. Exon size Sequences at exon–intron junctions (splice acceptor//splice donor)
1 184 Promoter //TCCCAGCAGgtgctgccc
2 267 taaacatagGACATGGAT//TACTCCCAGgtacactcg
3 90 ttcctccagGTTGCCCAC//GACCATGAGgtataaaaa
4 168 gcatttcagTCGCTTGCC//GAATTACAGgtaagagag
5 1671 ctttcctagTGTCTTAAT//ACTTAAACT_ACTTTGG.._

m_NRF2_ Expression in Tissues.

Northern blot analysis showed that like the human NRF2, m_NRF2_ is also expressed in a wide range of tissues at different levels in different tissues (Fig. 2 Left). Northern blot analysis of mouse fetal RNA shows significant amounts of NRF2 mRNA from developmental stage of day 7 to 17. (Fig. 2 Right). In addition, reverse transcription–PCR analysis of ES cells RNA revealed expression of NRF2, indicating that_NRF2_ is expressed from cells at the blastocyst stage (data not shown). In situ hybridization (Fig.3) using 35S-labeled riboprobe on section of 15.5-day-old embryos shows that although NRF2 is expressed at low levels in many tissues, it is expressed at high levels in several organs such as the liver, lungs, and kidneys. The highest level of expression is found in the cells lining the digestive tract, extending from the esophagus to the small and large intestines. In the central nervous system, NRF2 is expressed in the lateral ventricle and its medial wall and the choroid plexus in the fourth ventricle, which is the site of the future pituitary. NRF2 expression can also be found in the olfactory epithelium in the nasal cavity, the thyroid and submandibular glands, and the brown fat layer in the back. On higher magnification, NRF2 expression is most pronounced in the luminal cells of the stomach and intestine, the lining of the bronchi and alveoli and the renal tubules (Fig. 4).

Figure 2.

Figure 2

Expression of the mNRF2. Adult mouse multiple tissue RNA blots (Left) and mouse embryos RNA blots (Right) hybridized with mNRF2 and β-actin.

Figure 3.

Figure 3

Localization of NRF2 mRNA in the mouse embryo.In situ hybridization was performed on neighboring sections using antisense (Left) and sense (Right) probe. Two saggital sections are shown under darkfield illumination. Li, liver; Lu, lung; K, kidney; I, intestines; S, stomach; Oe, esophagus entrance (nasopharynx); C, choroid plexus; Sm, submandibular gland; O, Olfactory epithelium; T, testis; Bf, brown fat; > points to the thyroid.

Figure 4.

Figure 4

Localization of Nrf2 mRNA in the day 16 embryo developing stomach (A), intestine (B), kidney (C), and lung (D) under darkfield illumination. m, Muscle; l, lumen; k, kidney; a, adrenal gland; b, bronchus; v, blood vessel. (Bar = 100 μm.)

Production of _NRF2_−/− Mice.

A replacement type targeting vector was constructed as shown (Fig.5A). The 5′ and 3′ homologous genomic sequences were 5.5 kb and 6.2 kb long, respectively. In this targeting vector, a 4.2-kb segment of DNA containing part of exon 4 and all of exon 5 of NRF2 is replaced by a 5.5-kb fragment of DNA coding for the bacterial gene LacZ followed by the positive selection neomycin resistance cassette. This replacement vector was engineered to delete the carboxyl 457 amino acid of the NRF2 molecule and to splice in the bacterial reporter gene LacZ. This design would effectively nullify the NRF2 gene function, since the CNC bZIP regions are deleted, and the construct provides a means of monitoring the_NRF2_ promoter activity. DNA from ES cell clones was analyzed by Southern blotting to screen for homologous recombinants (Fig. 5 B and C). Of the 60 clones analyzed, 10_NRF2_+/− heterozygous ES clones were identified, thus a targeting efficiency of 1 in 6 was achieved in JM1 ES cells.

Figure 5.

Figure 5

(A) Genomic sequence and targeting constructs with the resultant NRF2 knockout gene. (B) Southern blot analysis of ES cells. An external probe_Xba_I–_Hin_dIII fragment gave a 15.5-kb band for the wild-type and 16.5-kb band for the homologous recombinant. (C) Southern blot analysis of F2 offspring. A_Stu_I–_Eco_RV probe gave a 2.4-kb band for the wild type and a 1.5-kb band for the recombinant.

Chimeric mice were generated, and germ-line transmission was detected. Heterozygous F1 mice appeared normal and were interbred to generate homozygous _NRF2_−/− mice. Of 63 offspring, 17 (27%) were homozygous for the disrupted allele and 17 (27%) were wild type, indicating a normal Mendelian inheritance pattern.

Analysis of _NRF2_−/− Mice.

The homozygous mice were normal in appearance and showed no discernible difference in their development or in their behavior compared with their wild-type littermates. At 6 weeks of age, wild-type,NRF2+/−, and _NRF2_−/− were killed for gross and histological examination. No obvious anatomical differences in the major organs were seen. Histological examination of the stomach, intestine, and kidney did not reveal any abnormality in the heterozygous and homozygous animals (data not shown). Hematological indices of _NRF2_−/− mice were all within the range of heterozygous and wild-type littermates (data not shown). Attempts at staining adult tissues with 5-bromo-4-chloro-3-indolyl β-d-galactoside (X-Gal) failed to detect any β-galactosidase activity. Homozygous and heterozygous embryos isolated at 13.5 and 15.5 days of gestation appeared to have developed normally and synchronously with wild-type littermates up to 9 months of age.

RNA Analysis of _NRF2_−/− Mice.

In situ hybridizations were carried out on paraffin section of adult stomach and intestine from wild-type and_NRF2_−/− mice using a riboprobe synthesized from a fragment of cDNA that was deleted in the targeted allele. No hybridizations were detected on sections of stomach and intestine from_NRF2_−/− mice (data not shown). poly(A) RNA from intestine, kidney, stomach, and testis were isolated from wild-type, NRF2+/−, and_NRF2_−/− mice and analyzed by Northern blotting. A 500-bp cDNA probe (_Bsr_DI–_Eco_RV) hybridized to 2.5-kb message on the NRF2+/+ RNA, to 2.5- and 4.0-kb messages in the NRF2+/− RNA, and only to the 4.0-kb message in the _NRF2_−/− RNA (Fig. 6). A lacZ probe only hybridized to this 4.0-kb message (data not shown). This larger 4.0-kb transcript produced from the mutated allele appears to be unstable, as it hybridized much less intensely than the 2.5 kb RNA, and was barely discernible in the RNA from some tissues.

Figure 6.

Figure 6

Northern blot analysis of RNA from wild-type, NRF2+/−, and NRF2−/− mice. The NRF2 message is 2.5 kb, and the NRF2-LacZ chimeric message is 4 kb.

DISCUSSION

Our analysis on DNA from the targeted ES cells as well as the DNA from the ES cells derived mice shows that we have disrupted the_NRF2_ gene. Several facts indicate that the modification in the NRF2 gene results in a null mutation. The LacZ cassette linked to the neomycin cassette in the targeting construct replaced part of exon 4 and the whole of exon 5, which contain the CNC, DNA binding, and leucine zipper domain, leaving about 125 amino acids of the amino end of the NRF2 protein linked in-frame to the LacZ reporter gene. This extensive deletion is borne out by Northern blot analysis, reverse transcription–PCR, and in situ hybridization. From these observations, we conclude that no functional_NRF2_ protein could possibly be made in the Nrf2 null mice.

It was unfortunate that the lacZ gene fused to the_NRF2_ gene failed to generate any detectable β-galactosidase activity for it could have provided a means of monitoring the expression of NRF2 in the mice throughout their development. A test construct using a cDNA that contained exon 2 to part of exon 4 fused to the same LacZ cassette driven by the cytomegalovirus promoter successfully showed β-galactosidase activity upon X-Gal staining when transiently transfected in 293 kidney cells. It is possible that the 16 extra amino acids coded by exon 1 in addition to the one coded by exon 2 to exon 4 created enough steric bulk to abolish the activity of the chimeric lacZ molecule. Alternately, as shown by the Northern blot data, chimeric NRF2-LacZ transcript is much less abundant than the wild type and is unstable.

As NRF2 is expressed widely, and even in ES cells, one might expect that the homozygous null mice would exhibit defect in embryogenesis. Yet, interbreeding of heterozygous animals resulted in a normal Mendelian ratio of homozygous null offspring, which are indistinguishable from their heterozygous and normal littermates. They are fertile and produce normal litter size. These results show that NRF2 gene function is not necessary for mouse development, growth, and fertility.

It is interesting that using a p45-NFE2 probe to screen a chicken cDNA library, Itoh and coworkers (33) isolated a cDNA, termed_Ech_, that encodes a bZIP protein more homologous to NRF2 than p45-NFE2. Of the 40 of 42 ECH clones Itoh and coworkers (33) identified, all were homologous to NRF2 and none to_p45-NFE2_. They also showed that this chicken_NRF2_-like gene is highly expressed in hematopoietic cells as well as some other tissues. It heterodimerizes with MAF protein and transactivates through the NFE2 motifs. Perhaps NRF2 serves important hematopoietic functions in the avian system, but is replaced by p45-NFE2 or NRF1 in mammals through evolution.

Although the NRF2 null mice are viable and exhibit an outwardly normal phenotype, we cannot rule out that lack of NRF2 results in subtle changes that do not produce a recognizable phenotype under laboratory conditions, but will only reveal changes that can be detected when subjected to specific challenges or when the mice age. NRF2 absence may also have consequences in compound knock-out animals involving other CNC-bZIP genes. Experiments are in progress to address these issues.

Acknowledgments

We thank Dr. Jefferson Chan for reviewing the manuscript, Dr. T. S. B. Yen for reviewing the sections from the mice, Mei-Chi Cheung for her help in the art work, and Lisa Woldin for editorial assistance. The work is supported in part by National Institutes of Health Grants DK16666 and DK50267. Y.W.K. is an Investigator of the Howard Hughes Medical Institute.

Footnotes

The publication costs of this article were defrayed in part by page charge payment. This article must therefore be hereby marked “_advertisement_” in accordance with 18 U.S.C. §1734 solely to indicate this fact.

Abbreviations: CNC, cap ‘n’ collar; bZIP, basic region–leucine zipper; RACE, rapid amplification of cDNA 5′ ends; ES cells, embryonic stem cells.

Data deposition: The sequences reported in this paper have been deposited in the GenBank data base (accession nos. U70474U70474 andU70475U70475).

References