Integrative analysis of gene expression, DNA methylation, physiological traits, and genetic variation in human skeletal muscle.

We integrate comeasured gene expression and DNA methylation (DNAme) in 265 human skeletal muscle biopsies from the FUSION study with >7 million genetic variants and eight physiological traits: height, waist, weight, waist-hip ratio, body mass index, fasting serum insulin, fasting plasma glucose, and type 2 diabetes. We find hundreds of genes and DNAme sites associated with fasting insulin, waist, and body mass index, as well as thousands of DNAme sites associated with gene expression (eQTM). We find that controlling for heterogeneity in tissue/muscle fiber type reduces the number of physiological trait associations, and that long-range eQTMs (>1 Mb) are reduced when controlling for tissue/muscle fiber type or latent factors. We map genetic regulators (quantitative trait loci; QTLs) of expression (eQTLs) and DNAme (mQTLs). Using Mendelian randomization (MR) and mediation techniques, we leverage these genetic maps to predict 213 causal relationships between expression and DNAme, approximately two-thirds of which predict methylation to causally influence expression. We use MR to integrate FUSION mQTLs, FUSION eQTLs, and GTEx eQTLs for 48 tissues with genetic associations for 534 diseases and quantitative traits. We identify hundreds of genes and thousands of DNAme sites that may drive the reported disease/quantitative trait genetic associations. We identify 300 gene expression MR associations that are present in both FUSION and GTEx skeletal muscle and that show stronger evidence of MR association in skeletal muscle than other tissues, which may partially reflect differences in power across tissues. As one example, we find that increased RXRA muscle expression may decrease lean tissue mass.

Interactions between genetic variation and cellular environment in skeletal muscle gene expression.

From whole organisms to individual cells, responses to environmental conditions are influenced by genetic makeup, where the effect of genetic variation on a trait depends on the environmental context. RNA-sequencing quantifies gene expression as a molecular trait, and is capable of capturing both genetic and environmental effects. In this study, we explore opportunities of using allele-specific expression (ASE) to discover cis-acting genotype-environment interactions (GxE)-genetic effects on gene expression that depend on an environmental condition. Treating 17 common, clinical traits as approximations of the cellular environment of 267 skeletal muscle biopsies, we identify 10 candidate environmental response expression quantitative trait loci (reQTLs) across 6 traits (12 unique gene-environment trait pairs; 10% FDR per trait) including sex, systolic blood pressure, and low-density lipoprotein cholesterol. Although using ASE is in principle a promising approach to detect GxE effects, replication of such signals can be challenging as validation requires harmonization of environmental traits across cohorts and a sufficient sampling of heterozygotes for a transcribed SNP. Comprehensive discovery and replication will require large human transcriptome datasets, or the integration of multiple transcribed SNPs, coupled with standardized clinical phenotyping.

The effect of age and gender on the genetic regulation of serum 25-hydroxyvitamin D - the FIN-D2D population-based study.

In addition to sunlight and dietary sources, several genes in the metabolic pathway of vitamin D affect serum 25-hydroxyvitamin D (25OHD) concentration. It is not known whether this genetic regulation is influenced by host characteristics. We investigated the effect of age and gender on the genetic regulation of serum 25OHD concentration. In total, 2868 Finnish men and women aged 45-74 years participated in FIN-D2D population-based health survey in 2007. Of the 2822 participants that had serum 25OHD concentration available, 2757 were successfully genotyped. Age and gender-dependent association of SNPs with serum 25OHD concentration was studied in 10 SNPs with previously found association with vitamin D metabolites. Associations of 3 SNPs with serum 25OHD concentration were dependent on age with greater effects on younger (≤60 y) than older (>60 y) adults (rs10783219 in VDR, rs12512631 in GC and rs3794060 in NADSYN1/DHCR7; pinteraction = 0.03, 0.02 and 0.01, respectively). The results suggested a novel association between serum 25OHD concentration and rs8082391 in STAT5B gene in men but not in women (pinteraction = 0.01). After multiple testing correction with false discovery rate method, two age-dependent interactions (rs3794060 in NADSYN1/DHCR7 gene and rs12512631 in GC gene) remained statistically significant. This is the first study to suggest that genetic regulation of serum 25OHD concentration is age-dependent. Our results also indicated a novel association between serum 25OHD concentration and SNP in STAT5B gene in men. The results need to be confirmed in future studies preferably in a larger sample.

Genetic determinants of serum 25-hydroxyvitamin D concentration during pregnancy and type 1 diabetes in the child.

OBJECTIVE: The in utero environment plays an important role in shaping development and later life health of the fetus. It has been shown that maternal genetic factors in the metabolic pathway of vitamin D associate with type 1 diabetes in the child. In this study we analyzed the genetic determinants of serum 25-hydroxyvitamin D (25OHD) concentration during pregnancy in mothers whose children later developed type 1 diabetes and in control mothers. STUDY DESIGN: 474 mothers of type 1 diabetic children and 348 mothers of non-diabetic children were included in the study. We previously selected 7 single nucleotide polymorphisms (SNPs) in four genes in the metabolic pathway of vitamin D vitamin based on our previously published data demonstrating an association between genotype and serum 25OHD concentration. In this re-analysis, possible differences in strength in the association between the SNPs and serum 25OHD concentration in mothers of type 1 diabetic and non-diabetic children were investigated. Serum 25OHD concentrations were previously shown to be similar between the mothers of type 1 diabetic and non-diabetic children and vitamin D deficiency prevalent in both groups. RESULTS: Associations between serum 25OHD concentration and 2 SNPs, one in the vitamin D receptor (VDR) gene (rs4516035) and one in the group-specific component (GC) gene (rs12512631), were stronger during pregnancy in mothers whose children later developed type 1 diabetes than in mothers whose children did not (pinteraction = 0.03, 0.02, respectively). CONCLUSIONS: We show for the first time that there are differences in the strength of genetic determinants of serum 25OHD concentration during pregnancy between the mothers of type 1 diabetic and non-diabetic children. Our results emphasize that the in utero environment including maternal vitamin D metabolism should be important lines of investigation when searching for factors that lead to early programming of type 1 diabetes.

Genome-wide meta-analysis of 241,258 adults accounting for smoking behaviour identifies novel loci for obesity traits.

Few genome-wide association studies (GWAS) account for environmental exposures, like smoking, potentially impacting the overall trait variance when investigating the genetic contribution to obesity-related traits. Here, we use GWAS data from 51,080 current smokers and 190,178 nonsmokers (87% European descent) to identify loci influencing BMI and central adiposity, measured as waist circumference and waist-to-hip ratio both adjusted for BMI. We identify 23 novel genetic loci, and 9 loci with convincing evidence of gene-smoking interaction (GxSMK) on obesity-related traits. We show consistent direction of effect for all identified loci and significance for 18 novel and for 5 interaction loci in an independent study sample. These loci highlight novel biological functions, including response to oxidative stress, addictive behaviour, and regulatory functions emphasizing the importance of accounting for environment in genetic analyses. Our results suggest that tobacco smoking may alter the genetic susceptibility to overall adiposity and body fat distribution.

The genetic regulatory signature of type 2 diabetes in human skeletal muscle.

Type 2 diabetes (T2D) results from the combined effects of genetic and environmental factors on multiple tissues over time. Of the >100 variants associated with T2D and related traits in genome-wide association studies (GWAS), >90% occur in non-coding regions, suggesting a strong regulatory component to T2D risk. Here to understand how T2D status, metabolic traits and genetic variation influence gene expression, we analyse skeletal muscle biopsies from 271 well-phenotyped Finnish participants with glucose tolerance ranging from normal to newly diagnosed T2D. We perform high-depth strand-specific mRNA-sequencing and dense genotyping. Computational integration of these data with epigenome data, including ATAC-seq on skeletal muscle, and transcriptome data across diverse tissues reveals that the tissue-specific genetic regulatory architecture of skeletal muscle is highly enriched in muscle stretch/super enhancers, including some that overlap T2D GWAS variants. In one such example, T2D risk alleles residing in a muscle stretch/super enhancer are linked to increased expression and alternative splicing of muscle-specific isoforms of ANK1.

Maternal VDR variants rather than 25-hydroxyvitamin D concentration during early pregnancy are associated with type 1 diabetes in the offspring.

AIMS/HYPOTHESIS: We investigated whether single nucleotide polymorphisms (SNPs) associated with 25-hydroxyvitamin D concentration in the metabolic pathway of vitamin D show different genotype distributions between Finnish families with an offspring with type 1 diabetes (cases) and families with a healthy offspring (controls). METHODS: A total of 31 SNPs in eight genes were studied in case and control mothers and family members (offspring with type 1 diabetes and healthy siblings, healthy control children and fathers) (n = 2,854). The 25-hydroxyvitamin D concentration was studied in 474 case and 348 matched control mothers during pregnancy. RESULTS: The genotype distributions of 13 SNPs (in the following genes: 7-dehydrocholesterol reductase NADSYN1/DHCR7, vitamin D receptor VDR, group-specific component GC and CYP27A1) that showed a nominal association with 25-hydroxyvitamin D concentration (p < 0.05) were compared between case and control families. SNPs in VDR had different genotype distributions between the case and control mothers (rs1544410, p = 0.007; rs731236, p = 0.003; rs4516035, p = 0.015), two SNPs (rs1544410 and rs731236) remaining significant after correction for multiple testing using a false discovery rate. The mean 25-hydroxyvitamin D concentrations during pregnancy did not differ between the case and control mothers. CONCLUSIONS/INTERPRETATION: Our preliminary results suggest that the maternal genotypes of SNPs in VDR may influence the in utero environment and thus contribute to the early programming of type 1 diabetes in the fetus. It is possible that the effects are only relevant in the presence of vitamin D insufficiency.

Age-related clonal hematopoiesis associated with adverse outcomes.

BACKGROUND: The incidence of hematologic cancers increases with age. These cancers are associated with recurrent somatic mutations in specific genes. We hypothesized that such mutations would be detectable in the blood of some persons who are not known to have hematologic disorders. METHODS: We analyzed whole-exome sequencing data from DNA in the peripheral-blood cells of 17,182 persons who were unselected for hematologic phenotypes. We looked for somatic mutations by identifying previously characterized single-nucleotide variants and small insertions or deletions in 160 genes that are recurrently mutated in hematologic cancers. The presence of mutations was analyzed for an association with hematologic phenotypes, survival, and cardiovascular events. RESULTS: Detectable somatic mutations were rare in persons younger than 40 years of age but rose appreciably in frequency with age. Among persons 70 to 79 years of age, 80 to 89 years of age, and 90 to 108 years of age, these clonal mutations were observed in 9.5% (219 of 2300 persons), 11.7% (37 of 317), and 18.4% (19 of 103), respectively. The majority of the variants occurred in three genes: DNMT3A, TET2, and ASXL1. The presence of a somatic mutation was associated with an increase in the risk of hematologic cancer (hazard ratio, 11.1; 95% confidence interval [CI], 3.9 to 32.6), an increase in all-cause mortality (hazard ratio, 1.4; 95% CI, 1.1 to 1.8), and increases in the risks of incident coronary heart disease (hazard ratio, 2.0; 95% CI, 1.2 to 3.4) and ischemic stroke (hazard ratio, 2.6; 95% CI, 1.4 to 4.8). CONCLUSIONS: Age-related clonal hematopoiesis is a common condition that is associated with increases in the risk of hematologic cancer and in all-cause mortality, with the latter possibly due to an increased risk of cardiovascular disease. (Funded by the National Institutes of Health and others.).

Association of serum 25-hydroxyvitamin D with lifestyle factors and metabolic and cardiovascular disease markers: population-based cross-sectional study (FIN-D2D).

OBJECTIVES: Low serum 25-hydroxyvitamin D (25OHD) level has been associated with an increased risk of several chronic diseases. Our aim was to determine lifestyle and clinical factors that are associated with 25OHD level and to investigate connection of 25OHD level with metabolic and cardiovascular disease markers. DESIGN: In total, 2868 Finnish men and women aged 45-74 years participated in FIN-D2D population-based health survey in 2007. Participants that had a serum sample available (98.4%; n = 2822) were included in this study. 25OHD was measured with chemiluminescent microparticle immunoassay method. RESULTS: The mean 25OHD level was 58.2 nmol/l in men (n = 1348) and 57.1 nmol/l in women (n = 1474). Mean 25OHD level was lower in the younger age groups than in the older ones (p<0.0001 both in men and women). This study confirmed that low physical activity (p<0.0001 both in men and women), smoking (p = 0.0002 in men and p = 0.03 in women) and high BMI (p<0.0001 in women) are factors that independently associate with low 25OHD level. Of the metabolic and cardiovascular disease markers high triglyceride concentration (p = 0.02 in men and p = 0.001 in women) and high apolipoprotein B/apolipoprotein A1 ratio (p = 0.04 in men and p = 0.03 in women) were independently associated with low 25OHD level. CONCLUSIONS: Higher age did not predict lower 25OHD level in this study population of aged 45-74 years which may derive from a healthy life-style of "active pensioners". Low physical activity and smoking came up as independent lifestyle factors associated with low 25OHD level. Defining the molecular mechanisms behind the associations of 25OHD with low physical activity and smoking are important objective in future studies. The association of 25OHD with BMI, high triglyceride concentration and apolipoprotein B/apolipoprotein A1 ratio may be related to the role of vitamin D in inflammation, but more detailed studies are needed.

Trans-ethnic fine-mapping of lipid loci identifies population-specific signals and allelic heterogeneity that increases the trait variance explained.

Genome-wide association studies (GWAS) have identified ~100 loci associated with blood lipid levels, but much of the trait heritability remains unexplained, and at most loci the identities of the trait-influencing variants remain unknown. We conducted a trans-ethnic fine-mapping study at 18, 22, and 18 GWAS loci on the Metabochip for their association with triglycerides (TG), high-density lipoprotein cholesterol (HDL-C), and low-density lipoprotein cholesterol (LDL-C), respectively, in individuals of African American (n = 6,832), East Asian (n = 9,449), and European (n = 10,829) ancestry. We aimed to identify the variants with strongest association at each locus, identify additional and population-specific signals, refine association signals, and assess the relative significance of previously described functional variants. Among the 58 loci, 33 exhibited evidence of association at P<1 × 10(-4) in at least one ancestry group. Sequential conditional analyses revealed that ten, nine, and four loci in African Americans, Europeans, and East Asians, respectively, exhibited two or more signals. At these loci, accounting for all signals led to a 1.3- to 1.8-fold increase in the explained phenotypic variance compared to the strongest signals. Distinct signals across ancestry groups were identified at PCSK9 and APOA5. Trans-ethnic analyses narrowed the signals to smaller sets of variants at GCKR, PPP1R3B, ABO, LCAT, and ABCA1. Of 27 variants reported previously to have functional effects, 74% exhibited the strongest association at the respective signal. In conclusion, trans-ethnic high-density genotyping and analysis confirm the presence of allelic heterogeneity, allow the identification of population-specific variants, and limit the number of candidate SNPs for functional studies.

Large-scale association analysis provides insights into the genetic architecture and pathophysiology of type 2 diabetes.

To extend understanding of the genetic architecture and molecular basis of type 2 diabetes (T2D), we conducted a meta-analysis of genetic variants on the Metabochip, including 34,840 cases and 114,981 controls, overwhelmingly of European descent. We identified ten previously unreported T2D susceptibility loci, including two showing sex-differentiated association. Genome-wide analyses of these data are consistent with a long tail of additional common variant loci explaining much of the variation in susceptibility to T2D. Exploration of the enlarged set of susceptibility loci implicates several processes, including CREBBP-related transcription, adipocytokine signaling and cell cycle regulation, in diabetes pathogenesis.

Single-nucleotide polymorphisms in the IL2RA gene are associated with age at diagnosis in late-onset Finnish type 1 diabetes subjects.

The onset of type 1 diabetes can occur at any age, with as many as half of all cases diagnosed after age 15. Despite this wide distribution in age at diagnosis, most genetic studies focus on cases diagnosed in childhood or during early adulthood. To better understand the genetics of late-onset type 1 diabetes, we collected a Finnish case/control cohort with all cases diagnosed between ages 15 and 40. We genotyped 591 probands and 1,538 control subjects at regions well established as susceptibility loci in early onset type 1 diabetes. These loci were then tested for disease association and age-at-diagnosis effects. Using logistic regression, we found that single-nucleotide polymorphisms (SNPs) at the INS, PTPN22, and IFIH1 loci were associated with late-onset disease (OR (95%CI) = 0.57(0.47-0.69), p = 2.77 x 10(-9); OR (95%CI) = 1.50 (1.27-1.78), p = 3.98 x 10(-6); and OR (95%CI) = 0.81(0.71-0.93), p = 0.0028, respectively). In contrast, a disease association was not detected for two SNPs at the IL2RA locus (rs11594656 and rs41295061). Despite this, we did find an independent age-at-diagnosis effect for each IL2RA SNP using a multivariate Cox proportional hazards model (p = 0.003, 0.002, respectively). Taken together, polymorphisms at the IL2RA locus were a major determinant of age at diagnosis in our cohort with an effect at par with the HLA-DQ2/DQ8 genotype as measured by hazard ratios. These findings suggest that the IL2RA locus controls both the susceptibility to disease and its time of occurrence. Thus, we believe the IL2/IL2R axis represents a potential therapeutic target for delaying the onset of disease.

A genome-wide scan for type 1 diabetes susceptibility genes in nuclear families with multiple affected siblings in Finland.

BACKGROUND: A genome-wide search for genes that predispose to type 1 diabetes using linkage analysis was performed using 900 microsatellite markers in 70 nuclear families with affected siblings from Finland, a population expected to be more genetically homogeneous than others, and having the highest incidence of type 1 diabetes in the world and, yet, the highest proportion in Europe of cases (10%) carrying neither of the highest risk HLA haplotypes that include DR3 or DR4 alleles. RESULTS: In addition to the evidence of linkage to the HLA region on 6p21 (nominal p = 4.0 x 10-6), significant evidence of linkage in other chromosome regions was not detected with a single-locus analysis. The two-locus analysis conditional on the HLA gave a maximum lod score (MLS) of 3.1 (nominal p = 2 x 10-4) on chromosome 9p13 under an additive model; MLS of 2.1 (nominal p = 6.1 x 10-3) on chromosome 17p12 and MLS of 2.5 (nominal p = 2.9 x 10-3) on chromosome 18p11 under a general model. CONCLUSION: Our genome scan data confirmed the primary contribution of the HLA genes also in the high-risk Finnish population, and suggest that non-HLA genes also contribute to the familial clustering of type 1 diabetes in Finland.

A genome-wide association study of type 2 diabetes in Finns detects multiple susceptibility variants.

Identifying the genetic variants that increase the risk of type 2 diabetes (T2D) in humans has been a formidable challenge. Adopting a genome-wide association strategy, we genotyped 1161 Finnish T2D cases and 1174 Finnish normal glucose-tolerant (NGT) controls with >315,000 single-nucleotide polymorphisms (SNPs) and imputed genotypes for an additional >2 million autosomal SNPs. We carried out association analysis with these SNPs to identify genetic variants that predispose to T2D, compared our T2D association results with the results of two similar studies, and genotyped 80 SNPs in an additional 1215 Finnish T2D cases and 1258 Finnish NGT controls. We identify T2D-associated variants in an intergenic region of chromosome 11p12, contribute to the identification of T2D-associated variants near the genes IGF2BP2 and CDKAL1 and the region of CDKN2A and CDKN2B, and confirm that variants near TCF7L2, SLC30A8, HHEX, FTO, PPARG, and KCNJ11 are associated with T2D risk. This brings the number of T2D loci now confidently identified to at least 10.

Nuclear and nucleolar targeting of influenza A virus NS1 protein: striking differences between different virus subtypes.

Influenza A virus nonstructural protein 1 (NS1A protein) is a virulence factor which is targeted into the nucleus. It is a multifunctional protein that inhibits host cell pre-mRNA processing and counteracts host cell antiviral responses. We show that the NS1A protein can interact with all six human importin alpha isoforms, indicating that the nuclear translocation of NS1A protein is mediated by the classical importin alpha/beta pathway. The NS1A protein of the H1N1 (WSN/33) virus has only one N-terminal arginine- or lysine-rich nuclear localization signal (NLS1), whereas the NS1A protein of the H3N2 subtype (Udorn/72) virus also has a second C-terminal NLS (NLS2). NLS1 is mapped to residues 35 to 41, which also function in the double-stranded RNA-binding activity of the NS1A protein. NLS2 was created by a 7-amino-acid C-terminal extension (residues 231 to 237) that became prevalent among human influenza A virus types isolated between the years 1950 to 1987. NLS2 includes basic amino acids at positions 219, 220, 224, 229, 231, and 232. Surprisingly, NLS2 also forms a functional nucleolar localization signal NoLS, a function that was retained in H3N2 type virus NS1A proteins even without the C-terminal extension. It is likely that the evolutionarily well-conserved nucleolar targeting function of NS1A protein plays a role in the pathogenesis of influenza A virus.

How well does the metabolic syndrome defined by five definitions predict incident diabetes and incident coronary heart disease in a Chinese population?

We evaluate the ability of the metabolic syndrome (MetS) defined by five definitions for predicting both incident CHD and diabetes combined, diabetes alone, and CHD alone in a Chinese population. The screening survey for type 2 diabetes was conducted in 1994. A follow-up study of 541 high-risk non-diabetic individuals who were free of CHD at baseline was carried out in 1999 in Beijing area. The MetS was defined by the World Health Organization (WHO), European Group for the Study of Insulin Resistance (EGIR), American College of Endocrinology (ACE), the International Diabetes Federation (IDF), and the National Cholesterol Education Program and the American Heart Association (AHA) (updated NCEP) criteria. From a multiple logistic regression adjusting for age, sex, education, occupation, smoking, family history of diabetes, and total cholesterol, the relative risk of the ACE-defined MetS for incident diabetes alone (67 cases) was 2.29 (95% CI, 1.20-4.34). The MetS defined by the five definitions was associated with a 1.8-3.9 times increased risk for both incident CHD and diabetes combined (59 cases), and with a 1.9-3.0 times for total incident diabetes (126 cases). None of the five definitions predicted either incident CHD alone (177 cases) or total incident CHD (236 cases). In conclusion, the MetS defined by the current definitions appears to be more effective at predicting incident diabetes.

Tag SNP selection for Finnish individuals based on the CEPH Utah HapMap database.

The pattern and nature of linkage disequilibrium in the human genome is being studied and catalogued as part of the International HapMap Project [:2003 Nature 426:789-796]. A key goal of the HapMap Project is to enable identification of tag single nucleotide polymorphisms (SNPs) that capture a substantial portion of common human genetic variability while requiring only a small fraction of SNPs to be genotyped [International HapMap Consortium, 2005: Nature 437:1299-1320]. In the current study, we examined the effectiveness of using the CEU HapMap database to select tag SNPs for a Finnish sample. We selected SNPs in a 17.9-Mb region of chromosome 14 based on pairwise linkage disequilibrium (r(2)) estimates from the HapMap CEU sample, and genotyped 956 of these SNPs in 1,425 Finnish individuals. An excess of SNPs showed significantly different allele frequencies between the HapMap CEU and the Finnish samples, consistent with population-specific differences. However, we observed strong correlations between the two samples for estimates of allele frequencies, r(2) values, and haplotype frequencies. Our results demonstrate that the HapMap CEU samples provide an adequate basis for tag SNP selection in Finnish individuals, without the need to create a map specifically for the Finnish population, and suggest that the four-population HapMap data will provide useful information for tag SNP selection beyond the specific populations from which they were sampled.

Differential expression of IFN regulatory factor 4 gene in human monocyte-derived dendritic cells and macrophages.

In vitro human monocyte differentiation to macrophages or dendritic cells (DCs) is driven by GM-CSF or GM-CSF and IL-4, respectively. IFN regulatory factors (IRFs), especially IRF1 and IRF8, are known to play essential roles in the development and functions of macrophages and DCs. In the present study, we performed cDNA microarray and Northern blot analyses to characterize changes in gene expression of selected genes during cytokine-stimulated differentiation of human monocytes to macrophages or DCs. The results show that the expression of IRF4 mRNA, but not of other IRFs, was specifically up-regulated during DC differentiation. No differences in IRF4 promoter histone acetylation could be found between macrophages and DCs, suggesting that the gene locus was accessible for transcription in both cell types. Computer analysis of the human IRF4 promoter revealed several putative STAT and NF-kappaB binding sites, as well as an IRF/Ets binding site. These sites were found to be functional in transcription factor-binding and chromatin immunoprecipitation experiments. Interestingly, Stat4 and NF-kappaB p50 and p65 mRNAs were expressed at higher levels in DCs as compared with macrophages, and enhanced binding of these factors to their respective IRF4 promoter elements was found in DCs. IRF4, together with PU.1, was also found to bind to the IRF/Ets response element in the IRF4 promoter, suggesting that IRF4 protein provides a positive feedback signal for its own gene expression in DCs. Our results suggest that IRF4 is likely to play an important role in myeloid DC differentiation and gene regulatory functions.

NF-{kappa}B is transported into the nucleus by importin {alpha}3 and importin {alpha}4.

NF-kappaB transcription factors are retained in the cytoplasm in an inactive form until they are activated and rapidly imported into the nucleus. We identified importin alpha3 and importin alpha4 as the main importin alpha isoforms mediating TNF-alpha-stimulated NF-kappaB p50/p65 heterodimer translocation into the nucleus. Importin alpha3 and alpha4 are close relatives in the human importin alpha family. We show that importin alpha3 isoform also mediates nuclear import of NF-kappaB p50 homodimer in nonstimulated cells. Importin alpha3 is shown to directly bind to previously characterized nuclear localization signals (NLSs) of NF-kappaB p50 and p65 proteins. Importin alpha molecules are known to have armadillo repeats that constitute the N-terminal and C-terminal NLS binding sites. We demonstrate by site-directed mutagenesis that NF-kappaB p50 binds to the N-terminal and p65 to the C-terminal NLS binding site of importin alpha3. In vitro competition experiments and analysis of cellular NF-kappaB suggest that NF-kappaB binds to importin alpha only when it is free of IkappaBalpha. The present study demonstrates that the nuclear import of NF-kappaB is a highly regulated process mediated by a subset of importin alpha molecules.