- About this Journal ·
- Abstracting and Indexing ·
- Aims and Scope ·
- Annual Issues ·
- Article Processing Charges ·
- Articles in Press ·
- Author Guidelines ·
- Bibliographic Information ·
- Citations to this Journal ·
- Contact Information ·
- Editorial Board ·
- Editorial Workflow ·
- Free eTOC Alerts ·
- Publication Ethics ·
- Reviewers Acknowledgment ·
- Submit a Manuscript ·
- Subscription Information ·
- Table of Contents
Evidence-Based Complementary and Alternative Medicine
Volume 2011 (2011), Article ID 508689, 8 pages
Gene Expression Profile of the Hypothalamus in DNP-KLH Immunized Mice Following Electroacupuncture Stimulation
1Department of Physiology, College of Oriental Medicine, Kyung Hee University, #1 Hoegi-dong, Dongdaemoon-gu, Seoul 130-701, Republic of Korea
2BK21 Oriental Medical Science Center, Republic of Korea
3Acupuncture & Meridian Science Research Center, Republic of Korea
4Department of Acupuncture & Moxibustion, College of Oriental Medicine, Kyung Hee University, #1 Hoeki-dong, Dongdaemoon-gu, Seoul, 130-701, Republic of Korea
5Department of Oriental Gynecology, College of Oriental Medicine, Republic of Korea
6Department of Physiology, College of Medicine, Kyung Hee University, Seoul, Republic of Korea
Received 5 August 2009; Accepted 25 November 2009
Copyright © 2011 Sun Kwang Kim et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Clinical evidence indicates that electroacupuncture (EA) is effective for allergic disorder. Recent animal studies have shown that EA treatment reduces levels of IgE and Th2 cytokines in BALB/c mice immunized with 2,4-dinitrophenylated keyhole limpet protein (DNP-KLH). The hypothalamus, a brain center of the neural-immune system, is known to be activated by EA stimulation. This study was performed to identify and characterize the differentially expressed genes in the hypothalamus of DNP-KLH immunized mice that were stimulated with EA or only restrained. To this aim, we conducted a microarray analysis to evaluate the global gene expression profiles, using the hypothalamic RNA samples taken from three groups of mice: (i) normal control group (no treatments); (ii) IMH group (DNP-KLH immunization + restraint); and (iii) IMEA group (immunization + EA stimulation). The microarray analysis revealed that total 39 genes were altered in their expression levels by EA treatment. Ten genes, including T-cell receptor alpha variable region family 13 subfamily 1 (Tcra-V13.1), heat shock protein 1B (Hspa1b) and 2′–5′ oligoadenylate synthetase 1F (Oas1f), were up-regulated in the IMEA group when compared with the IMH group. In contrast, 29 genes, including decay accelerating factor 2 (Daf2), NAD(P)H dehydrogenase, quinone 1 (Nqo1) and programmed cell death 1 ligand 2 (Pdcd1lg2) were down-regulated in the IMEA group as compared with the IMH group. These results suggest that EA treatment can modulate immune response in DNP-KLH immunized mice by regulating expression levels of genes that are associated with innate immune, cellular defense and/or other kinds of immune system in the hypothalamus.
Several clinical studies have previously suggested that electroacupuncture (EA) is effective for the treatment of various types of allergic disorders, including asthma and chronic urticaria [1–4]. It is generally accepted that hyperproduction of IgE in response to stimulation by Th2 cytokines (i.e., IL-4, IL-5 and IL-13) promotes the development of these allergic disorders [5, 6]. Our previous study demonstrated that successive EA at the ST36 acupoint reduces serum levels of IgE in BALB/c mice that have been immunized with 2,4-dinitrophenylated keyhole limpet protein (DNP-KLH) via suppression of the production of Th2 cytokines . Furthermore, it has been shown that these effects of EA are mediated, at least in part, by -adrenergic receptors  and are acupoint (ST36) specific .
The hypothalamus is the primary center for neuroendocrine-immune modulation . The paraventricular nucleus of the hypothalamus (PVN) represents an integral part of the neuroendocrine circuit that modulates immune function. Moreover, stimulation of the ventral nucleus of the hypothalamus (VMH) stimulation enhances the immune function associated with the lateral hypothalamus (LH) and ventral tegmental area (VTA) of the brain. Interestingly, previous brain imaging studies of animals and humans have shown that EA treatment affects the neuronal activity of the hypothalamus [11–16]. Although it has also been reported that EA influences the immunomodulatory effects of the hypothalamus [17, 18], the studies that have been conducted to date have only investigated the effects of EA on normal rats.
Oligonucleotide microarray analysis is the most powerful tool for functional genomics that provides direct information about mRNA expression levels from a large number of genes. This method has been used to find the genes correlated with the pathogenesis of diseases and to assess the drug actions on the diseases [19, 20]. To examine if and how EA treatment modulates immune function of the hypothalamus at transcriptional level in a Th2-skewed condition, we conducted a microarray analysis on the hypothalamus of DNP-KLH immunized mice that were stimulated with EA.
Eight-week-old female BALB/c mice were purchased from Orient Bio (Sungnam, Korea) and housed in an air-controlled, pathogen-free animal facility under a 12-h light/dark cycle at 23 ± 2°C. Mice were provided with a standard laboratory diet and water ad libitum. All procedures involving animals were conducted in accordance with the NIH guidelines.
2.2. Experimental Groups
Mice were randomly divided into the following 3 groups ( = 6–10/group): (i) NC group (normal control, no treatments for 2 weeks, except for i.p. injection of saline on the 1st and 8th experimental day); (ii) IMH group (immunization with 2,4-dinitrophenylated keyhole limpet protein (DNP-KLH) on the 1st and 8th experimental day and being held in a holder restraint for 20 min a day for 2 weeks; (iii) IMEA group (immunization on the 1st and 8th experimental day + daily EA stimulation at the ST36 acupoint while being held in a holder restraint for 20 min a day for 2 weeks).
2.3. EA Stimulation and DNP-KLH Immunization
EA stimulation was performed as previously described . Briefly, a pair of stainless steel needles (0.2 mm in diameter and 3 cm long) was inserted into the ST36, which is located at the anterior tibial muscle, 5 mm laterally and lower from the anterior tubercle of the tibia, and at a point 5 mm distal from the first needle. The anode and cathode leads from an electrical stimulator were then connected to the two acupuncture needles, after which train-pulses (1 Hz, 0.25 ms pulse width, 3–5 V) were applied for 20 min. The stimulation intensity was determined to be the minimum voltage that induced a moderate muscle contraction.
All mice, except for the NC group, were immunized intraperitoneally with 4 g of DNP-KLH (Calbiochem, Gibbstown, NJ, USA) and 4 mg of aluminum hydroxide (Sigma, St Louis, MO, USA) on the 1st and 8th day during the 2-week experimental period. Immunization of BALB/c mice with DNP-KLH caused a significant increase in IgE and IL-4 production on the 14th day after first immunization; therefore, EA stimulation was conducted daily for 2 weeks .
2.4. Measurements of IgE in Serum Using ELISA
Serum was obtained from all groups of mice on the final experimental day and then stored at 20°C until analysis. The Serum level of IgE was measured using a quantitative sandwich enzyme-linked immunoassay kit (BD, Franklin Lakes, NJ, USA). Briefly, a 96-well microtiter plate (Coster, Cambridge, MA, USA) was incubated overnight at 4°C with anti-mouse IgE monoclonal antibody in coating buffer. The plate was then washed with PBS containing 0.05 tween 20 (Sigma, USA) and blocked with 5 FBS in PBS for 1 h at room temperature. Subsequently, 100 l of sample were loaded onto the plate, which was then incubated for 2 h at room temperature. Next, secondary peroxidase labeled biotinylated anti-mouse IgE monoclonal antibody in 5 FBS in PBS diluent was added and the plate was incubated for an additional 1 h. The plates were then treated with TMB substrate solution for 30 min, after which the reaction was stopped by the addition 50 l of 2N H2SO4 stop solution per well. Finally, the optical density at 450 nm was measured in a microplate reader (TECAN, Durham, NC, USA).
2.5. Oligonucleotide Chip Microarray
The hypothalamuses were obtained from all groups of mice on the final experimental day and then stored at 80°C until analysis. Total RNA was extracted from the hypothalamus by using TRIZOL reagent, and DNase I (Invitrogen Life Technologies, Rockville, MD, USA) was treated to eliminate genomic DNA contamination. Oligonucleotide chip microarray was performed using single round RNA amplification protocols, following the Affymetrix specifications (Affymetrix GeneChip Expression Analysis Technical Manual). Briefly, 3 g of total RNA was used to synthesize first-strand complementary DNA (cDNA) using oligonucleotide probes with 24 oligo-dT plus T7 promoter as primers (Proligo LLC, Boulder, CO, USA) and the Superscript Choice System (Life Technologies, Invitrogen, Milan, Italy). Following double-stranded cDNA synthesis, the products were purified by phenol-chloroform extraction, after which biotinylated antisense complimentary RNA (cRNA) was generated through in vitro transcription using a BioArray RNA High-Yield Transcript Labeling Kit (ENZO Life Sciences Inc., Farmingdale, NY, USA). Next, the biotinylated labeled cRNA was fragmented, and then 10 g of the total fragmented cRNA was hybridized to the GeneChip Mouse Gene 1.0 ST Array (P/N901178, Affymetrix Inc., Santa Clara, CA, USA). The Affymetrix Fluidics Station 400 was then used to wash and stain the chips, after which the nonhybridized target was removed. The samples were then incubated with a streptavidin-phycoerythrin conjugate to stain the biotinylated cRNA. Next, the staining was amplified using goat IgG as a blocking reagent and biotinylated antistreptavidin antibody (goat), followed by a second staining step using a streptavidin-phycoerythrin conjugate. The fluorescence was then detected using the Genechip System Confocal Scanner (Hewlett-Packard, Santa Clara, USA), and analysis of the data contained on each GeneChip was conducted using the GeneChip 3.1 software produced by Affymetrix using the default settings. Global scaling was then used to compare the chips, with all probe sets being scaled to a user-defined target intensity of 150.
2.6. Data Analysis
One-way ANOVA followed by Newman–Keuls multiple comparison test was used for statistical analysis of the IgE data. In all cases, P .05 was considered significant. Data are presented as the mean ± SEM.
For analysis of the microarray data, the MAS5 algorithm was used to evaluate the expression signals generated by the Affymetrix Mouse Gene 1.0 ST Array. Global scaling normalization was then performed, after which the normalized data were log-transformed using base 2. Next, a fold change and a Welch t-test were applied to select the differentially expressed genes (DEGs) using a fold change threshold of 1.2 and a P .05 to indicate significance. Each probe set used in the Affymetrix GeneChip produces a detection call, with P (present call) indicating good quality, M (marginal call) indicating intermediate quality and A (absent call) indicating relatively low reliability. Therefore, probe sets that resulted in A calls were removed to filter false positives. A volcano plot was used to better visualize and compare the two DEG methods. The 1.2-fold DEGs were clustered using the GenPle v2.3 software (ISTECH Inc., Ilsan, Korea) using hierarchical clustering with Pearson correlation as a similarity measure and complete linkage as the linkage method. In addition, gene ontology significance analysis was conducted to investigate the functional relationships among the 1.2-fold DEGs using high-throughput GoMiner. The 1.2-fold DEGs were then mapped to the relevant pathways using the GenPlex v2.4 software (ISTECH Inc., Ilsan, Korea). The pathway resources were provided by the KEGG database.
3.1. IgE Levels in DNP-KLH Immunized Mice Following EA Stimulation
BALB/c mice were immunized with DNP-KLH and aluminum hydroxide once a week. In the IMEA group, 20-min EA stimulation at the ST36 acupoint under restraint with an acryl holder was performed daily for 2 weeks, whereas mice in the IMH group were placed in the holder restraint for 20 min a day for 2 weeks, but not treated. After the 2-week experimental period, the serum levels of IgE were measured. Figure 1 shows serum levels of IgE in mice. The IgE levels in the IMH and the IMEA group were significantly increased when compared with the IgE levels of the NC group (P .001). However, EA treatment significantly reduced the elevated IgE levels (P .01, IMEA versus IMH).
3.2. DEGs in the Hypothalamus of Mice in the NC, IMH and IMEA Groups
The genes that were differentially expressed among the groups are listed in Table 1 In addition, the genes are grouped into functional categories and pathways based on the KEGG database in Table 2. There were 10 genes that were down-regulated in the IMH group (versus NC) and up-regulated in the IMEA group (versus IMH). Conversely, 17 genes were up-regulated in the IMH group (versus NC) and down-regulated in the IMEA group (versus IMH). Additionally, 12 genes were not significantly changed in the IMH group (versus NC), but were down-regulated in the IMEA group (versus IMH).
3.3. Clustering Analysis and Gene Ontology Analysis
As shown in Figure 2, microarray data from the normal control (NC) group and the experimental (IMH and IMEA) groups were combined and clustered using hierarchical clustering. To examine the functional relationships among the 1.2-fold DEGs, the difference in gene ontology was determined using significance analysis (Figure 3). The results revealed that genes that were up- and down-regulated in the experimental groups were associated with processes including signal transduction, cellular defense response, transport, transcription, biological processes, translation, the cell cycle and T cell proliferation.
DNP-KLH immunization is known to induce the elevation of serum IgE levels and Th2 cytokines, which play important roles in the development of various allergic disorders [7, 21]. In the present study, BALB/c mice in the IMH and IMEA groups were immunized with DNP-KLH and aluminum hydroxide. Although both groups showed a significant increase in IgE levels when compared with the NC group, IgE levels in the IMEA group were significantly lower than those in the IMH group (Figure 1). This result is consistent with the results of previous studies [7–9] and strongly suggests that EA treatment can modulate the immune response in the Th2 dominant condition.
Microarray analysis has been shown to be useful for the simultaneous profiling of global gene expression and identification of new genes or new functions of known genes [19, 20]. As described in the introduction, the hypothalamus is a primary center in the central nervous system responsible for neural-immune modulation  and has been shown to be affected by EA treatment [11, 12, 15]. Moreover, our recent study  showed that EA stimulation at ST36 can alter the transcriptional activity of the hypothalamus in normal rats. Therefore, we conducted microarray analysis of the hypothalamus of DNP-KLH immunized mice to determine if EA treatment of mice in a Th2-skewed condition can regulate gene expression in the hypothalamus.
Ten genes, including T-cell receptor alpha variable region family 13 subfamily 1 (Tcra-V13.1), trefoil factor 2 (spasmolytic protein 1) (Tff2), heat shock protein 1B (Hspa1b) and 2′–5′ oligoadenylate synthetase 1F (Oas1f), were down-regulated in the IMH group when compared with the NC group and up-regulated in the IMEA group when compared with the IMH group. Tcra-V13.1 is a subfamily of the T-cell receptor V chain and contributes to the recognition of self-MHC molecules or peptides that has polymorphisms associated with diabetes [22, 23]. The expression pattern of Tff2, which is a member of the newly discovered protective factor of gastrointestinal mucous, may be associated with gastric cancer . Hspa1b is involved in innate immunity and diverse biological processes including anti-apoptosis, response to stress and DNA repair [25, 26]. Oas1f is a transcription unit of a family of interferon-induced antiviral proteins that may play a role in the host's innate defense mechanism against viral infection [27, 28]. The differential expression of these genes indicates that EA treatment may increase some activities of innate immune and cellular defense at the transcriptional level.
Conversely, 17 genes including decay accelerating factor 2 (Daf2), cytochrome P450, family 2, subfamily c, polypeptide 38 (Cyp2c38), myosin light chain 2, precursor lymphocyte-specific (Mylc2pl), NAD(P)H dehydrogenase, and quinone 1 (Nqo1) were up-regulated in the IMH group when compared with the NC group and downregulated in the IMEA group when compared with the IMH group. In addition, 12 genes, including programmed cell death 1 ligand 2 (Pdcd1lg2), were downregulated in the IMEA group when compared with the IMH group, but not significantly changed in the IMH group when compared with the NC group. Daf2 is a glycosylphosphatidylinositol (GPI)-anchored membrane regulator that inhibits both the classical and alternative pathways of complement activation . Cyp2c38 is a member of a superfamily of ubiquitous hemoproteins that metabolize a vast array of foreign chemicals as well as endogenous compounds . Mylc2pl is a regulatory myosin light chain gene that is expressed specifically in precursor B and T lymphocytes and may play a role in lymphocyte development . Nqo1 is a cytosolic enzyme that catalyzes the metabolic reduction of quinines and derivatives and is an endogenous factor in the regulation of immune response and autoimmunity . Pdcd1lg2 is a ligand for programmed cell death-1, which is a receptor that plays an inhibitory role in T cell activation and has a potential role in Th2-mediated diseases . As shown in Table 2 and Figure 3, it appears that no major or specific biochemical pathways are regulated by EA treatment. However, the present data suggest that EA treatment widely regulates immune response in the hypothalamus at the transcriptional level.
In conclusion, EA treatment was found to reduce the serum levels of IgE that were elevated by DNP-KLH immunization in BALB/c mice. Furthermore, EA treatment altered the expression levels of 39 genes in the hypothalamus, many of which are involved in immune response (Figure 4). Taken together, these results suggest that EA treatment may induce immunomodulation by regulating gene expressions in the hypothalamus of DNP-KLH immunized mice.
The Kyung Hee University Research Fund in 2008 (KHU-20080600).
S. K. Kim and J. Kim contribute equally to this work.
- W. Biernacki and M. D. Peake, “Acupuncture in treatment of stable asthma,” Respiratory Medicine, vol. 92, no. 9, pp. 1143–1145, 1998.
- C.-J. Chen and H.-S. Yu, “Acupuncture treatment of urticaria,” Archives of Dermatology, vol. 134, no. 11, pp. 1397–1399, 1998.
- J. Che, “The effect of acupuncture on serum IgE levels in patients with chronic urticaria,” Journal of Traditional Chinese Medicine, vol. 26, no. 3, pp. 189–190, 2006.
- T. Takishima, S. Mue, G. Tamura, T. Ishihara, and K. Watanabe, “The bronchodilating effect of acupuncture in patients with acute asthma,” Annals of Allergy, vol. 48, pp. 44–49, 1982.
- E. Maggie, “The Th1/Th2 paradigm in allergy,” Immunotechnology, vol. 3, pp. 233–244, 1998.
- J. A. Woodfolk, “Cytokines as a therapeutic target for allergic diseases: a complex picture,” Current Pharmaceutical Design, vol. 12, no. 19, pp. 2349–2363, 2006.
- M.-B. Park, E. Ko, C. Ahn et al., “Suppression of IgE production and modulation of Th1/Th2 cell response by electroacupuncture in DNP-KLH immunized mice,” Journal of Neuroimmunology, vol. 151, no. 1-2, pp. 40–44, 2004.
- Y. Lee, S. K. Kim, Y. Kim et al., “The α-adrenoceptor mediation of the immunomodulatory effects of electroacupuncture in DNP-KLH immunized mice,” Neuroscience Letters, vol. 423, no. 2, pp. 149–152, 2007.
- S. K. Kim, Y. Lee, H. Cho, et al., “A parametric study on the immunomodulatory effects of electroacupuncture in DNP-KLH immunized mice,” Evidence-Based Complementary and Alternative Medicine, 2009.
- D. Wrona, “Neural-immune interactions: an integrative view of the bidirectional relationship between the brain and immune systems,” Journal of Neuroimmunology, vol. 172, no. 1-2, pp. 38–58, 2006.
- J. H. Chiu, H. C. Cheng, C. H. Tai, et al., “Electroacupuncture-induced neural activation detected by use of manganese-enhanced functional magnetic resonance imaging in rabbits,” American Journal of Veterinary Research, vol. 62, pp. 178–182, 2001.
- J. H. Chiu, M. S. Chung, H. C. Cheng, et al., “Different central manifestations in response to electroacupuncture at analgesic and nonanalgesic acupoints in rats: a manganese-enhanced functional magnetic resonance imaging study,” Canadian Journal of Veterinary Research, vol. 67, pp. 94–101, 2003.
- J. C. Hsieh, C. H. Tu, F. P. Chen, et al., “Activation of the hypothalamus characterizes the acupuncture stimulation at the analgesic point in human: a positron emission tomography study,” Neuroscience Letters, vol. 307, pp. 105–108, 2001.
- K. K. S. Hui, J. Liu, N. Makris et al., “Acupuncture modulates the limbic system and subcortical gray structures of the human brain: evidence from fMRI studies in normal subjects,” Human Brain Mapping, vol. 9, no. 1, pp. 13–25, 2000.
- V. Napadow, N. Kettner, J. Liu et al., “Hypothalamus and amygdala response to acupuncture stimuli in carpal tunnel syndrome,” Pain, vol. 130, no. 3, pp. 254–266, 2007.
- M.-T. Wu, J.-C. Hsieh, J. Xiong et al., “Central nervous pathway for acupunture stimulation: localization of processing with functional MR imaging of the brain—preliminary experience,” Radiology, vol. 212, no. 1, pp. 133–141, 1999.
- G. S. Choi, S. D. Oha, J. B. Han et al., “Modulation of natural killer cell activity affected by electroacupuncture through lateral hypothalamic area in rats,” Neuroscience Letters, vol. 329, no. 1, pp. 1–4, 2002.
- S.-W. Rho, G.-S. Choi, E.-J. Ko et al., “Molecular changes in remote tissues induced by electro-acupuncture stimulation at acupoint ST36,” Molecules and Cells, vol. 25, no. 2, pp. 178–183, 2008.
- K. Izuhara and H. Saito, “Microarray-based identification of novel biomarkers in asthma,” Allergology International, vol. 55, no. 4, pp. 361–367, 2006.
- S. Thornton, D. Sowders, B. Aronow et al., “DNA microarray analysis reveals novel gene expression profiles in collagen-induced arthritis,” Clinical Immunology, vol. 105, no. 2, pp. 155–168, 2002.
- H. Oka, Y. Emori, N. Kobayashi, Y. Hayashi, and K. Nomoto, “Suppression of allergic reactions by royal jelly in association with the restoration of macrophage function and the improvement of Th1/Th2 cell responses,” International Immunopharmacology, vol. 1, pp. 521–532, 2001.
- E. A. Simone, Y. Liping, D. R. Wegmann, and G. S. Eisenbarth, “T cell receptor gene polymorphisms associated with anti-insulin, autoimmune T cells in diabetes-prone NOD mice,” The Journal of Immunology, vol. 10, pp. 317–321, 1997.
- C. G. DiRienzo, G. F. Murphy, T. M. Friedman, and R. Korngold, “T-cell receptor Vα usage by effector CD4+Vβ11+ T cells mediating graft-versus-host disease directed to minor histocompatibility antigens,” Biology of Blood and Marrow Transplantation, vol. 13, no. 3, pp. 265–276, 2007.
- S.-Q. Shi, J.-T. Cai, and J.-M. Yang, “Expression of trefoil factors 1 and 2 in precancerous condition and gastric cancer,” World Journal of Gastroenterology, vol. 12, no. 18, pp. 3119–3122, 2006.
- S. D. Klein and B. Brüne, “Heat-shock protein 70 attenuates nitric oxide-induced apoptosis in RAW macrophages by preventing cytochrome c release,” Biochemical Journal, vol. 362, no. 3, pp. 635–641, 2002.
- E. Karoly, A. Fekete, N. F. Banki et al., “Heat shock protein 72 (HSPA1B) gene polymorphism and toll-like receptor (TLR) 4 mutation are associated with increased risk of urinary tract infection in children,” Pediatric Research, vol. 61, no. 3, pp. 371–374, 2007.
- S. Eskildsen, R. Hartmann, N. O. Kjeldgaard, and J. Justesen, “Gene structure of the murine 2′-5′-oligoadenylate synthetase family,” Cellular and Molecular Life Sciences, vol. 59, no. 7, pp. 1212–1222, 2002.
- T. Mashimo, P. Glaser, M. Lucas et al., “Structural and functional genomics and evolutionary relationships in the cluster of genes encoding murine 2′,5′-oligoadenylate synthetases,” Genomics, vol. 82, no. 5, pp. 537–552, 2003.
- X. Sun, C. D. Funk, C. Deng, A. Sahu, J. D. Lambris, and W. C. Song, “Role of decay-accelerating factor in regulating complement activation on the erythrocyte surface as revealed by gene targeting,” Proceedings of the National Academy of Sciences of the United States of America, vol. 96, pp. 628–633, 1999.
- G. Luo, D. C. Zeldin, J. A. Blaisdell, E. Hodgson, and J. A. Goldstein, “Cloning, and expression of murine CYP2Cs and their ability to metabolize arachidonic acid,” Archives of Biochemistry and Biophysics, vol. 357, pp. 45–57, 1998.
- E. M. Oltz, G. D. Yancopoulos, M. A. Morrow, et al., “A novel regulatory myosin light chain gene distinguishes pre-B cell subsets and is IL-7 inducible,” The EMBO Journal, vol. 11, pp. 2759–2767, 1992.
- K. Iskander, J. Li, S. Han, B. Zheng, and A. K. Jaiswal, “NQO1 and NQO2 regulation of humoral immunity and autoimmunity,” Journal of Biological Chemistry, vol. 281, no. 41, pp. 30917–30924, 2006.
- E. Oflazoglu, D. A. Swart, P. Anders-Bartholo et al., “Paradoxical role of programmed death-1 ligand 2 in Th2 immune responses in vitro and in a mouse asthma model in vivo,” European Journal of Immunology, vol. 34, no. 12, pp. 3326–3336, 2004.