Genetic and Environmental Interaction in Type 1 Diabetes: a Relationship Between Genetic Risk Alleles and Molecular Traits of Enterovirus Infection?

PURPOSE OF REVIEW: We provide an overview of the current knowledge regarding the natural history of human type 1 diabetes (T1D) and the documented associations between virus infections (in particular the enteroviruses) and disease development. We review studies that examine whether T1D-specific risk alleles in genes involved in the function of the immune system can alter susceptibility to virus infections or affect the magnitude of the host antiviral response. We also highlight where the major gaps in our knowledge exist and consider possible implications that new insights gained from the discussed gene-environment interaction studies may bring. RECENT FINDINGS: A commonality between several of the studied T1D risk variants studied is their role in modulating the host immune response to viral infection. Generally, little support exists indicating that the risk variants increase susceptibility to infection and moreover, they usually appear to predispose the immune system towards a hyper-reactive state, decrease the risk of infection, and/or favor the establishment of viral persistence. In conclusion, although the current number of studies is limited, this type of research can provide important insights into the mechanisms that are central to disease pathogenesis and further describe how genetic and environmental factors jointly influence the risk of T1D development. The latter may provide genetic markers that could be used for patient stratification and for the selection of method(s) for disease prevention.

Genome-wide profiling of interleukin-4 and STAT6 transcription factor regulation of human Th2 cell programming.

Dissecting the molecular mechanisms by which T helper (Th) cells differentiate to effector Th2 cells is important for understanding the pathogenesis of immune-mediated diseases, such as asthma and allergy. Because the STAT6 transcription factor is an upstream mediator required for interleukin-4 (IL-4)-induced Th2 cell differentiation, its targets include genes important for this process. Using primary human CD4(+) T cells, and by blocking STAT6 with RNAi, we identified a number of direct and indirect targets of STAT6 with ChIP sequencing. The integration of these data sets with detailed kinetics of IL-4-driven transcriptional changes showed that STAT6 was predominantly needed for the activation of transcription leading to the Th2 cell phenotype. This integrated genome-wide data on IL-4- and STAT6-mediated transcription provide a unique resource for studies on Th cell differentiation and, in particular, for designing interventions of human Th2 cell responses.

Early T helper cell programming of gene expression in human.

Molecular mechanisms guiding naïve T helper cell differentiation into functionally specified effector cells are intensively studied. The rapidly growing knowledge is mainly achieved by using mouse cells or disease models. Comparatively exiguous data is gathered from human primary cells although they provide the "ultimate model" for immunology in man, have been exploited in many original studies paving the way for the field, and can be analyzed more easily than ever with the help of modern technology and methods. As usage of mouse models is unavoidable in translational research, parallel human and mouse studies should be performed to assure the relevancy of the hypothesis created during the basic research. In this review, we give an overview on the status of the studies conducted with human primary cells aiming at elucidating the mechanisms instructing the priming of T helper cell subtypes. The special emphasis is given to the recent high-throughput studies. In addition, by comparing the human and mouse studies we intend to point out the regulatory mechanisms and questions which are lacking examination with human primary cells.

Gut-adipose tissue axis in hepatic fat accumulation in humans.

BACKGROUND & AIMS: Recent evidence suggests that in animals gut microbiota composition (GMC) affects the onset and progression of hepatic fat accumulation. The aim of this study was to investigate in humans whether subjects with high hepatic fat content (HHFC) differ in their GMC from those with low hepatic fat content (LHFC), and whether these differences are associated with body composition, biomarkers and abdominal adipose tissue inflammation. METHODS: Hepatic fat content (HFC) was measured using proton magnetic resonance spectroscopy ((1)H MRS). Fecal GMC was profiled by 16S rRNA fluorescence in situ hybridization and flow cytometry. Adipose tissue gene expression was analyzed using Affymetrix microarrays and quantitative PCR. RESULTS: The HHFC group had unfavorable GMC described by lower amount of Faecalibacterium prausnitzii (FPrau) (p<0.05) and relatively higher Enterobacteria than the LHFC group. Metabolically dysbiotic GMC associated with HOMA-IR and triglycerides (p<0.05 for both). Several inflammation-related adipose tissue genes were differentially expressed and correlated with HFC (p<0.05). In addition, the expression of certain genes correlated with GMC dysbiosis, i.e., low FPrau-to-Bacteroides ratio. CONCLUSIONS: HHFC subjects differ unfavorably in their GMC from LHFC subjects. Adipose tissue inflammation may be an important link between GMC, metabolic disturbances, and hepatic fat accumulation.

Identification of early gene expression changes during human Th17 cell differentiation.

Th17 cells play an essential role in the pathogenesis of autoimmune and inflammatory diseases. Most of our current understanding on Th17 cell differentiation relies on studies carried out in mice, whereas the molecular mechanisms controlling human Th17 cell differentiation are less well defined. In this study, we identified gene expression changes characterizing early stages of human Th17 cell differentiation through genome-wide gene expression profiling. CD4(+) cells isolated from umbilical cord blood were used to determine detailed kinetics of gene expression after initiation of Th17 differentiation with IL1ß, IL6, and TGFß. The differential expression of selected candidate genes was further validated at protein level and analyzed for specificity in initiation of Th17 compared with initiation of other Th subsets, namely Th1, Th2, and iTreg. This first genome-wide profiling of transcriptomics during the induction of human Th17 differentiation provides a starting point for defining gene regulatory networks and identifying new candidates regulating Th17 differentiation in humans.

Coxsackievirus infection induces direct pancreatic ß cell killing but poor antiviral CD8+ T cell responses.

Coxsackievirus B (CVB) infection of pancreatic ß cells is associated with ß cell autoimmunity and type 1 diabetes. We investigated how CVB affects human ß cells and anti-CVB T cell responses. ß cells were efficiently infected by CVB in vitro, down-regulated human leukocyte antigen (HLA) class I, and presented few, selected HLA-bound viral peptides. Circulating CD8+ T cells from CVB-seropositive individuals recognized a fraction of these peptides; only another subfraction was targeted by effector/memory T cells that expressed exhaustion marker PD-1. T cells recognizing a CVB epitope cross-reacted with ß cell antigen GAD. Infected ß cells, which formed filopodia to propagate infection, were more efficiently killed by CVB than by CVB-reactive T cells. Our in vitro and ex vivo data highlight limited CD8+ T cell responses to CVB, supporting the rationale for CVB vaccination trials for type 1 diabetes prevention. CD8+ T cells recognizing structural and nonstructural CVB epitopes provide biomarkers to differentially follow response to infection and vaccination.

Coxsackievirus infection induces direct pancreatic ß-cell killing but poor anti-viral CD8+ T-cell responses.

Coxsackievirus B (CVB) infection of pancreatic ß cells is associated with ß-cell autoimmunity. We investigated how CVB impacts human ß cells and anti-CVB T-cell responses. ß cells were efficiently infected by CVB in vitro, downregulated HLA Class I and presented few, selected HLA-bound viral peptides. Circulating CD8+ T cells from CVB-seropositive individuals recognized only a fraction of these peptides, and only another sub-fraction was targeted by effector/memory T cells that expressed the exhaustion marker PD-1. T cells recognizing a CVB epitope cross-reacted with the ß-cell antigen GAD. Infected ß cells, which formed filopodia to propagate infection, were more efficiently killed by CVB than by CVB-reactive T cells. Thus, our in-vitro and ex-vivo data highlight limited T-cell responses to CVB, supporting the rationale for CVB vaccination trials for type 1 diabetes prevention. CD8+ T cells recognizing structural and non-structural CVB epitopes provide biomarkers to differentially follow response to infection and vaccination.

An integrative computational systems biology approach identifies differentially regulated dynamic transcriptome signatures which drive the initiation of human T helper cell differentiation.

BACKGROUND: A proper balance between different T helper (Th) cell subsets is necessary for normal functioning of the adaptive immune system. Revealing key genes and pathways driving the differentiation to distinct Th cell lineages provides important insight into underlying molecular mechanisms and new opportunities for modulating the immune response. Previous computational methods to quantify and visualize kinetic differential expression data of three or more lineages to identify reciprocally regulated genes have relied on clustering approaches and regression methods which have time as a factor, but have lacked methods which explicitly model temporal behavior. RESULTS: We studied transcriptional dynamics of human umbilical cord blood T helper cells cultured in absence and presence of cytokines promoting Th1 or Th2 differentiation. To identify genes that exhibit distinct lineage commitment dynamics and are specific for initiating differentiation to different Th cell subsets, we developed a novel computational methodology (LIGAP) allowing integrative analysis and visualization of multiple lineages over whole time-course profiles. Applying LIGAP to time-course data from multiple Th cell lineages, we identified and experimentally validated several differentially regulated Th cell subset specific genes as well as reciprocally regulated genes. Combining differentially regulated transcriptional profiles with transcription factor binding site and pathway information, we identified previously known and new putative transcriptional mechanisms involved in Th cell subset differentiation. All differentially regulated genes among the lineages together with an implementation of LIGAP are provided as an open-source resource. CONCLUSIONS: The LIGAP method is widely applicable to quantify differential time-course dynamics of many types of datasets and generalizes to any number of conditions. It summarizes all the time-course measurements together with the associated uncertainty for visualization and manual assessment purposes. Here we identified novel human Th subset specific transcripts as well as regulatory mechanisms important for the initiation of the Th cell subset differentiation.

SATB1 dictates expression of multiple genes including IL-5 involved in human T helper cell differentiation.

Special AT-rich binding protein 1 (SATB1) is a global chromatin organizer and a transcription factor regulated by interleukin-4 (IL-4) during the early T helper 2 (Th2) cell differentiation. Here we show that SATB1 controls multiple IL-4 target genes involved in human Th cell polarization or function. Among the genes regulated by SATB1 is that encoding the cytokine IL-5, which is predominantly produced by Th2 cells and plays a key role in the development of eosinophilia in asthma. We demonstrate that, during the early Th2 cell differentiation, IL-5 expression is repressed through direct binding of SATB1 to the IL-5 promoter. Furthermore, SATB1 knockdown-induced up-regulation of IL-5 is partly counteracted by down-regulating GATA3 expression using RNAi in polarizing Th2 cells. Our results suggest that a competitive mechanism involving SATB1 and GATA3 regulates IL-5 transcription, and provide new mechanistic insights into the stringent regulation of IL-5 expression during human Th2 cell differentiation.

Pancreas Whole Tissue Transcriptomics Highlights the Role of the Exocrine Pancreas in Patients With Recently Diagnosed Type 1 Diabetes.

Although type 1 diabetes (T1D) is primarily a disease of the pancreatic beta-cells, understanding of the disease-associated alterations in the whole pancreas could be important for the improved treatment or the prevention of the disease. We have characterized the whole-pancreas gene expression of patients with recently diagnosed T1D from the Diabetes Virus Detection (DiViD) study and non-diabetic controls. Furthermore, another parallel dataset of the whole pancreas and an additional dataset from the laser-captured pancreatic islets of the DiViD patients and non-diabetic organ donors were analyzed together with the original dataset to confirm the results and to get further insights into the potential disease-associated differences between the exocrine and the endocrine pancreas. First, higher expression of the core acinar cell genes, encoding for digestive enzymes, was detected in the whole pancreas of the DiViD patients when compared to non-diabetic controls. Second, In the pancreatic islets, upregulation of immune and inflammation related genes was observed in the DiViD patients when compared to non-diabetic controls, in line with earlier publications, while an opposite trend was observed for several immune and inflammation related genes at the whole pancreas tissue level. Third, strong downregulation of the regenerating gene family (REG) genes, linked to pancreatic islet growth and regeneration, was observed in the exocrine acinar cell dominated whole-pancreas data of the DiViD patients when compared with the non-diabetic controls. Fourth, analysis of unique features in the transcriptomes of each DiViD patient compared with the other DiViD patients, revealed elevated expression of central antiviral immune response genes in the whole-pancreas samples, but not in the pancreatic islets, of one DiViD patient. This difference in the extent of antiviral gene expression suggests different statuses of infection in the pancreas at the time of sampling between the DiViD patients, who were all enterovirus VP1+ in the islets by immunohistochemistry based on earlier studies. The observed features, indicating differences in the function, status and interplay between the exocrine and the endocrine pancreas of recent onset T1D patients, highlight the importance of studying both compartments for better understanding of the molecular mechanisms of T1D.

Quantitative proteomics reveals GIMAP family proteins 1 and 4 to be differentially regulated during human T helper cell differentiation.

T helper (Th) cells differentiate into functionally distinct effector cell subsets of which Th1 and Th2 cells are best characterized. Besides T cell receptor signaling, IL-12-induced STAT4 and T-bet- and IL-4-induced STAT6 and GATA3 signaling pathways are the major players regulating the Th1 and Th2 differentiation process, respectively. However, there are likely to be other yet unknown factors or pathways involved. In this study we used quantitative proteomics exploiting cleavable ICAT labeling and LC-MS/MS to identify IL-4-regulated proteins from the microsomal fractions of CD4(+) cells extracted from umbilical cord blood. We were able to identify 557 proteins of which 304 were also quantified. This study resulted in the identification of the down-regulation of small GTPases GIMAP1 and GIMAP4 by IL-4 during Th2 differentiation. We also showed that both GIMAP1 and GIMAP4 genes are up-regulated by IL-12 and other Th1 differentiation-inducing cytokines in cells induced to differentiate toward Th1 lineage and down-regulated by IL-4 in cells induced to Th2. Our results indicate that the GIMAP (GTPase of the immunity-associated protein) family of proteins is differentially regulated during Th cell differentiation.

Identification of novel Stat6 regulated proteins in IL-4-treated mouse lymphocytes.

Interleukin 4 (IL-4) has an indispensable role in the differentiation of naive T helper (Th) cells toward the Th2 phenotype and induction of B cells to produce the IgE class of Igs. By regulating these two cell types, IL-4 has a pre-eminent role in regulation of allergic inflammation. IL-4-mediated regulation of T and B cell functions is largely transmitted through signal transducer and activator of transcription 6 (Stat6). In this study, we have used metabolic labeling and 2-D electrophoresis to detect differences in the proteomes of IL-4 stimulated spleen mononuclear cells of Stat6-/- and wild type mice and MS/MS for protein identification. With this methodology, we identified 49 unique proteins from 21 protein spots to be differentially expressed. Interestingly, in Stat6-/- CD4(+) cells the expression of isoform 2 of core binding factor b (CBFb2) was enhanced. CBFb is a non-DNA binding cofactor for the Runx family of transcription factors, which have been implicated in regulation of Th cell differentiation. We also found cellular nucleic acid protein (CNBP) to be downregulated in Stat6-/- cells. None of the proteins identified in this study have previously been reported to be regulated via Stat6. The results highlight the importance of exploiting proteomics tools to complement the studies on Stat6 target genes identified through transcriptional profiling.

A practical comparison of methods for detecting transcription factor binding sites in ChIP-seq experiments.

BACKGROUND: Chromatin immunoprecipitation coupled with massively parallel sequencing (ChIP-seq) is increasingly being applied to study transcriptional regulation on a genome-wide scale. While numerous algorithms have recently been proposed for analysing the large ChIP-seq datasets, their relative merits and potential limitations remain unclear in practical applications. RESULTS: The present study compares the state-of-the-art algorithms for detecting transcription factor binding sites in four diverse ChIP-seq datasets under a variety of practical research settings. First, we demonstrate how the biological conclusions may change dramatically when the different algorithms are applied. The reproducibility across biological replicates is then investigated as an internal validation of the detections. Finally, the predicted binding sites with each method are compared to high-scoring binding motifs as well as binding regions confirmed in independent qPCR experiments. CONCLUSIONS: In general, our results indicate that the optimal choice of the computational approach depends heavily on the dataset under analysis. In addition to revealing valuable information to the users of this technology about the characteristics of the binding site detection approaches, the systematic evaluation framework provides also a useful reference to the developers of improved algorithms for ChIP-seq data.

Early Detection of Peripheral Blood Cell Signature in Children Developing ß-Cell Autoimmunity at a Young Age.

The appearance of type 1 diabetes (T1D)-associated autoantibodies is the first and only measurable parameter to predict progression toward T1D in genetically susceptible individuals. However, autoantibodies indicate an active autoimmune reaction, wherein the immune tolerance is already broken. Therefore, there is a clear and urgent need for new biomarkers that predict the onset of the autoimmune reaction preceding autoantibody positivity or reflect progressive ß-cell destruction. Here we report the mRNA sequencing-based analysis of 306 samples including fractionated samples of CD4+ and CD8+ T cells as well as CD4-CD8- cell fractions and unfractionated peripheral blood mononuclear cell samples longitudinally collected from seven children who developed ß-cell autoimmunity (case subjects) at a young age and matched control subjects. We identified transcripts, including interleukin 32 (IL32), that were upregulated before T1D-associated autoantibodies appeared. Single-cell RNA sequencing studies revealed that high IL32 in case samples was contributed mainly by activated T cells and NK cells. Further, we showed that IL32 expression can be induced by a virus and cytokines in pancreatic islets and ß-cells, respectively. The results provide a basis for early detection of aberrations in the immune system function before T1D and suggest a potential role for IL32 in the pathogenesis of T1D.

Th1 response and cytotoxicity genes are down-regulated in cutaneous T-cell lymphoma.

PURPOSE: Increased production of Th2 cytokines characterizes Sezary syndrome, the leukemic form of cutaneous T-cell lymphomas (CTCL). To identify the molecular background and to study whether shared by the most common CTCL subtype, mycosis fungoides, we analyzed the gene expression profiles in both subtypes. EXPERIMENTAL DESIGN: Freshly isolated cells from 30 samples, representing skin, blood, and enriched CD4(+) cell populations of mycosis fungoides and Sezary syndrome, were analyzed with Affymetrix (Santa Clara, CA) oligonucleotide microarrays, quantitative PCR, or immunohistochemistry. The gene expression profiles were combined with findings of comparative genomic hybridization of the same samples to identify chromosomal changes affecting the aberrant gene expression. RESULTS: We identified a set of Th1-specific genes [e.g., TBX21 (T-bet), NKG7, and SCYA5 (RANTES)] to be down-regulated in Sezary syndrome as well as in a proportion of mycosis fungoides samples. In both Sezary syndrome and mycosis fungoides blood samples, the S100P and LIR9 gene expression was up-regulated. In lesional skin, IL7R and CD52 were up-regulated. Integration of comparative genomic hybridization and transcriptomic data identified chromosome arms 1q, 3p, 3q, 4q, 12q, 16p, and 16q as likely targets for new CTCL-associated gene aberrations. CONCLUSIONS: Our findings revealed several new genes involved in CTCL pathogenesis and potential therapeutic targets. Down-regulation of a set of genes involved in Th1 polarization, including the major Th1-polarizing factor, TBX21, was for the first time associated with CTCL. In addition, a plausible explanation for the proliferative response of CTCL cells to locally produced interleukin-7 was revealed.

Genome-wide Analysis of STAT3-Mediated Transcription during Early Human Th17 Cell Differentiation.

The development of therapeutic strategies to combat immune-associated diseases requires the molecular mechanisms of human Th17 cell differentiation to be fully identified and understood. To investigate transcriptional control of Th17 cell differentiation, we used primary human CD4+ T cells in small interfering RNA (siRNA)-mediated gene silencing and chromatin immunoprecipitation followed by massive parallel sequencing (ChIP-seq) to identify both the early direct and indirect targets of STAT3. The integrated dataset presented in this study confirms that STAT3 is critical for transcriptional regulation of early human Th17 cell differentiation. Additionally, we found that a number of SNPs from loci associated with immune-mediated disorders were located at sites where STAT3 binds to induce Th17 cell specification. Importantly, introduction of such SNPs alters STAT3 binding in DNA affinity precipitation assays. Overall, our study provides important insights for modulating Th17-mediated pathogenic immune responses in humans.

Comparative analysis of human and mouse transcriptomes of Th17 cell priming.

Uncontrolled Th17 cell activity is associated with cancer and autoimmune and inflammatory diseases. To validate the potential relevance of mouse models of targeting the Th17 pathway in human diseases we used RNA sequencing to compare the expression of coding and non-coding transcripts during the priming of Th17 cell differentiation in both human and mouse. In addition to already known targets, several transcripts not previously linked to Th17 cell polarization were found in both species. Moreover, a considerable number of human-specific long non-coding RNAs were identified that responded to cytokines stimulating Th17 cell differentiation. We integrated our transcriptomics data with known disease-associated polymorphisms and show that conserved regulation pinpoints genes that are relevant to Th17 cell-mediated human diseases and that can be modelled in mouse. Substantial differences observed in non-coding transcriptomes between the two species as well as increased overlap between Th17 cell-specific gene expression and disease-associated polymorphisms underline the need of parallel analysis of human and mouse models. Comprehensive analysis of genes regulated during Th17 cell priming and their classification to conserved and non-conserved between human and mouse facilitates translational research, pointing out which candidate targets identified in human are worth studying by using in vivo mouse models.

NanoMiner - integrative human transcriptomics data resource for nanoparticle research.

The potential impact of nanoparticles on the environment and on human health has attracted considerable interest worldwide. The amount of transcriptomics data, in which tissues and cell lines are exposed to nanoparticles, increases year by year. In addition to the importance of the original findings, this data can have value in broader context when combined with other previously acquired and published results. In order to facilitate the efficient usage of the data, we have developed the NanoMiner web resource (http://nanominer.cs.tut.fi/), which contains 404 human transcriptome samples exposed to various types of nanoparticles. All the samples in NanoMiner have been annotated, preprocessed and normalized using standard methods that ensure the quality of the data analyses and enable the users to utilize the database systematically across the different experimental setups and platforms. With NanoMiner it is possible to 1) search and plot the expression profiles of one or several genes of interest, 2) cluster the samples within the datasets, 3) find differentially expressed genes in various nanoparticle studies, 4) detect the nanoparticles causing differential expression of selected genes, 5) analyze enriched Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways and Gene Ontology (GO) terms for the detected genes and 6) search the expression values and differential expressions of the genes belonging to a specific KEGG pathway or Gene Ontology. In sum, NanoMiner database is a valuable collection of microarray data which can be also used as a data repository for future analyses.

Gene expression profiling of immune-competent human cells exposed to engineered zinc oxide or titanium dioxide nanoparticles.

A comprehensive in vitro assessment of two commercial metal oxide nanoparticles, TiO2 and ZnO, was performed using human monocyte-derived macrophages (HMDM), monocyte-derived dendritic cells (MDDC), and Jurkat T cell leukemia-derived cell line. TiO2 nanoparticles were found to be non-toxic whereas ZnO nanoparticles caused dose-dependent cell death. Subsequently, global gene expression profiling was performed to identify transcriptional response underlying the cytotoxicity caused by ZnO nanoparticles. Analysis was done with doses 1 µg/ml and 10 µg/ml after 6 and 24 h of exposure. Interestingly, 2703 genes were significantly differentially expressed in HMDM upon exposure to 10 µg/ml ZnO nanoparticles, while in MDDCs only 12 genes were affected. In Jurkat cells, 980 genes were differentially expressed. It is noteworthy that only the gene expression of metallothioneins was upregulated in all the three cell types and a notable proportion of the genes were regulated in a cell type-specific manner. Gene ontology analysis revealed that the top biological processes disturbed in HMDM and Jurkat cells were regulating cell death and growth. In addition, genes controlling immune system development were affected. Using a panel of modified ZnO nanoparticles, we obtained an additional support that the cellular response to ZnO nanoparticles is largely dependent on particle dissolution and show that the ligand used to modify ZnO nanoparticles modulates Zn(2+) leaching. Overall, the study provides an extensive resource of transcriptional markers for mediating ZnO nanoparticle-induced toxicity for further mechanistic studies, and demonstrates the value of assessing nanoparticle responses through a combined transcriptomics and bioinformatics approach.

Genome-wide identification of novel genes involved in early Th1 and Th2 cell differentiation.

Th cell subtypes, Th1 and Th2, are involved in the pathogenesis or progression of many immune-mediated diseases, such as type 1 diabetes and asthma, respectively. Defining the molecular networks and factors that direct Th1 and Th2 cell differentiation will help to understand the pathogenic mechanisms causing these diseases. Some of the key factors regulating this differentiation have been identified, however, they alone do not explain the process in detail. To identify novel factors directing the early differentiation, we have studied the transcriptomes of human Th1 and Th2 cells after 2, 6, and 48 h of polarization at the genome scale. Based on our current and previous studies, 288 genes or expressed sequence tags, representing approximately 1-1.5% of the human genome, are regulated in the process during the first 2 days. These transcriptional profiles revealed genes coding for components of certain pathways, such as RAS oncogene family and G protein-coupled receptor signaling, to be differentially regulated during the early Th1 and Th2 cell differentiation. Importantly, numerous novel genes with unknown functions were identified. By using short-hairpin RNA knockdown, we show that a subset of these genes is regulated by IL-4 through STAT6 signaling. Furthermore, we demonstrate that one of the IL-4 regulated genes, NDFIP2, promotes IFN-gamma production by the polarized human Th1 lymphocytes. Among the novel genes identified, there may be many factors that play a crucial role in the regulation of the differentiation process together with the previously known factors and are potential targets for developing therapeutics to modulate Th1 and Th2 responses.