The Journal of Allergy and Clinical Immunology
Volume 115, Issue 1 , Pages 80-87, January 2005

Association analysis of common variants of STAT6, GATA3, and STAT4 to asthma and high serum IgE phenotypes

  • Maritta Pykäläinen, MSc

      Affiliations

    • Turku Centre for Biotechnology, University of Turku and Åbo Akademi University
    • Turku Graduate School of Biomedical Sciences
    • Corresponding Author InformationReprint requests: Maritta Pykäläinen, MSc, Turku Centre for Biotechnology, University of Turku and Åbo Akademi University, PO Box 124, FIN-20520, Turku, Finland.
  • ,
  • Riikka Kinos, PhD

      Affiliations

    • Turku Centre for Biotechnology, University of Turku and Åbo Akademi University
  • ,
  • Sanna Valkonen, MSc

      Affiliations

    • Turku Centre for Biotechnology, University of Turku and Åbo Akademi University
  • ,
  • Pia Rydman, PhD

      Affiliations

    • GeneOS Ltd, Helsinki
  • ,
  • Maritta Kilpeläinen, MD, PhD

      Affiliations

    • Department of Pulmonary Diseases and Clinical Allergology, Turku
  • ,
  • Lauri A. Laitinen, MD, PhD

      Affiliations

    • Department of Pulmonary Medicine, Helsinki University Central Hospital
  • ,
  • Jussi Karjalainen, MD, PhD

      Affiliations

    • Department of Respiratory Medicine
  • ,
  • Markku Nieminen, MD, PhD

      Affiliations

    • Department of Respiratory Medicine
  • ,
  • Mikko Hurme, MD, PhD

      Affiliations

    • Department of Microbiology and Immunology, Tampere University Hospital
  • ,
  • Juha Kere, MD, PhD

      Affiliations

    • Department of Medical Genetics, Biomedicum Helsinki, University of Helsinki
    • Department of Biosciences at Novum, Karolinska Institute, Stockholm
  • ,
  • Tarja Laitinen, MD, PhD

      Affiliations

    • GeneOS Ltd, Helsinki
    • Department of Medical Genetics, Biomedicum Helsinki, University of Helsinki
  • ,
  • Riitta Lahesmaa, MD, PhD

      Affiliations

    • Turku Centre for Biotechnology, University of Turku and Åbo Akademi University

Received 15 January 2004; received in revised form 27 September 2004; accepted 4 October 2004.

Turku, Helsinki, and Tampere, Finland, and Stockholm, Sweden

Article Outline

Background

Immune responses characterized by TH2 type cells and IgE are important for the development of asthma and atopy. The transcription factors STAT6, GATA3, and STAT4 mediate the cytokine-induced development of naïve CD4+ T cells into either TH1 or TH2 type.

Objective

We studied genetic variation of the STAT6, GATA3, and STAT4 genes and examined whether single nucleotide polymorphisms (SNPs) in these loci were associated with asthma or serum high IgE levels in the Finnish asthmatic families.

Methods

With denaturing high-performance liquid chromatography we screened all exons and exon-intron boundaries of the genes in 14 to 22 patients. All identified SNPs were genotyped in 120 nuclear families, and the haplotypes were analyzed by Haplotype Pattern Mining based statistical analysis. When potential association was observed, the analysis was replicated among 245 asthmatic patients and 405 population-based control subjects.

Results

A total of 23 SNPs were identified, of which 8 were not previously listed in the SNP database. Interestingly, a haplotype analysis of GATA3 showed 3 related haplotypes that associated with different asthma and atopy related phenotypes among both the family and case-control data sets. For STAT6 and STAT4, no significant association to asthma or serum total IgE levels was observed.

Conclusions

We identified a panel of novel SNPs in genes coding for proteins important in the TH1/TH2 cell differentiation. SNPs of the GATA3 gene showed an initial association to asthma-related phenotypes. Elucidation of the importance of the identified panel of SNPs in other TH1/TH2 mediated diseases will be of great interest.

Key words: Association analysis, common variant, STAT6, GATA3, STAT4, asthma, serum IgE, haplotype

Abbreviations used: DHPLC, Denaturing high-performance liquid chromatography, GATA3, GATA-binding protein 3, HPM, Haplotype Pattern Mining, ILR, Interleukin receptor, SNP, Single nucleotide polymorphism, STAT4, Signal transductor and activator of transcription 4, STAT6, Signal transductor and activator of transcription 6, UTR, Untranslated region

 

Specific activation of STAT6 is critical for a number of IL-4 mediated effects, including induction of TH2 and IgE responses, and antigen-induced airway inflammation and hyperresponsiveness.1, 2, 3, 4, 5 Human STAT6 (12q13.3-14.1) spans 19 kb of genomic DNA containing 23 exons.6 Several genome-wide linkage studies in both asthma and atopy have shown linkage to 12q, suggesting that STAT6 is a functional and positional candidate gene for asthma-related traits.7, 8, 9, 10, 11, 12 The most profound association study was performed among the white asthmatic sibpairs. In the study, 14 markers were screened, of which the one in the intron 18 showed an association (P=.007) to serum IgE level and the GT repeat polymorphism in the exon 1 (P=.001) to elevated eosinophils.13 The significance of STAT6 variations was investigated also in Japanese studies; one showed an association between the GT repeat and allergic phenotypes,14 and another reported the same repeat to be associated with Japanese childhood asthma.15 In another Japanese asthma study, an SNP in the 3′ untranslated region (UTR) of the STAT6 gene was shown to be associated with mild atopic asthma (P=.004); this association was not, however, replicated among the British patients16 or 2 other studies among Japanese patients.14, 15

Human GATA3 (10p15) consists of 6 exons spanning 19 kb of genomic DNA.17 GATA3 belongs to the highly conserved family of transcription activators of which only GATA3 is expressed in T cells.18, 19, 20 GATA3 is a key factor of TH2 cell development by regulating the production of TH2-specific cytokines IL-4, IL-5, and IL-13.21, 22, 23 GATA3 also has a STAT6-independent autoactivation pathway.24 Recently, an alternative form of GATA3 promoter and an alternative first exon (exon 1a) were identified.25

STAT4 deficient mice have markedly reduced IL-12–induced TH1 responses.26 While promoting TH1-type cell development STAT4 is a negative regulator for TH2 cell differentiation.27, 28 STAT4 spans 120 kb on chromosome 2q32.2-q32.3 across one large intron.29 The gene consists of 24 exons. Also 2q has been reported to be linked to asthma-related phenotypes.9, 30

In this study, we investigated the role of genetic variations in the 3 genes, STAT6, STAT4, and GATA3, central for T-helper cell differentiation in the predisposition to asthma or high levels of serum IgE in the Finnish asthmatic families.

Back to Article Outline

Methods 

Study subjects in SNP identification 

Fourteen Finnish unrelated study subjects were screened to detect polymorphisms in the STAT6 gene and 22 subjects for STAT4 and GATA3 loci. Seven were non-atopic control individuals (average IgE, 143 kU/L; range, 50-310 kU/L). They had no history of asthma or atopy, had normal pulmonary function, showed no bronchial hyperreactivity in histamine challenge test, and had negative RAST and skin prick test results. The rest were young atopic adults (average IgE, 4697 kU/L; range, 1,143 to 27,716 kU/L) of which 10 had asthmatic symptoms. These patients were tested as the non-atopic group. They had a history of asthma and atopy, a clinically diagnosed significant reversible bronchus obstruction, and/or bronchial hyperreactivity during histamine challenge test, and at least one positive RAST and/or skin prick test result. Genomic DNA was extracted from peripheral white blood cells by propanol/ethanol-extraction. Serum total IgE level was measured by Diagnostics CAP FEIA (Kabi Pharmacia, Uppsala, Sweden).

Study subjects in the association analyses 

For association study, our DNA material consisted of Finnish study families recruited on the basis of asthma. The details of the recruitment and diagnostic testing have been described previously.31, 32 Phenotypes available for this cohort (average age, 43 years; 37 years among individuals with high and 45 years among individuals with low serum IgE) were asthma status, serum total IgE level (Diagnostics CAP FEIA), and blood eosinophil level. Among the genotyped samples there were no non-atopic individuals with high IgE values. For STAT6, STAT4, and GATA3, an association analysis was performed with a subset of the cohort including all the trios with full phenotypic and genotypic information (a total of 120 trios and 360 study subjects). For GATA3 the initial association was replicated in a population-based case-control sample including 245 Finnish asthmatic patients and 405 control subjects. Details of the recruitment and clinical evaluation of the patients have been described.33 Studied phenotypes for this cohort were asthma status, serum total IgE level, skin prick test results, and blood eosinophil level. Total serum IgE was measured by the immunoluminometric method (Ciba Corning Diagnostics, Halsted, United Kingdom). Genomic DNA of all study individuals was extracted from peripheral white blood cells by using a standard nonenzymatic method.

SNP screening 

For STAT6, GATA3, and STAT4, all the exons, 1 to 2 kb of the 5′ and 3′ flanking sequences of the genes, and intron-exon boundaries were screened for polymorphisms by using temperature-modulated heteroduplex analysis based denaturing high-performance liquid chromatography (DHPLC) method and carried out on the automatic Wave DNA Fragment Analysis System (Transgenomic, San Jose, Calif).34 The overlapping PCR fragments covered 9.8 kb (~51.6%) of the STAT6, 8 kb (~42.1 %) of the GATA3, and 11 kb (~9.1 %) of STAT4 loci. Screening covered Gene bank sequences NT_029419 for STAT6, NT_077569 for GATA3, and NT_005403 for STAT4. PCR primers (given as requested) were designed by the Primer Express software (Applied Biosystems, Foster City, Calif); each amplicon was 300 to 800 bp. When possible, the amplicons were designed to overlap 50 to 200 bp. PCR reactions were carried out in 50- to 100-μL volumes containing 100 to 200 ng of template DNA, 0.2 mmol/L deoxyribonucleoside triphosphate (MBI Fermentas, Vilnius, Lithuania), 400 to 800 nmol/L of each PCR primer (Medprobe, Oslo, Norway), 2 to 3.5 U of Dynazyme EXT polymerase, 10 X Dynazyme Tritonfree polymerase buffer (Finnzymes, Espoo, Finland), and 5% dimethyl sulfoxide (Finnzymes). For DHPLC analysis the melting temperatures for each fragment were predicted with WAVEMAKER software (Transgenomic) and Stanford's DHPLC Melt Program (http://insertion.stanford.edu/melt.html). Observed heteroduplexes were confirmed with direct sequencing, both the group of samples giving the heteroduplex as well as those giving homoduplex profile (ABI PRISM 3100 Genetic Analyzer using BigDye Terminator Cycle Sequencing Kit; Applied Biosystems).

SNP genotyping 

Genotyping was performed by SNaPshot (Applied Biosystems)35 according to the manufacturer's directions. Primer sequences are available as requested. PCR samples were purified with Shrimp Alkaline Phosphatase (SAP) and Exonuclease I (ExoI) enzymes (UBS, Cleveland, Ohio). The SnaPshot reactions were carried out in multiplex reactions in a total volume of 10 μL containing 5 μL of ready reaction mix (diluted 2:3 with 20 mmol/L Mg2+ and 800 mmol/L Tris-HCl, pH 9.0), 0.1 to 2.8 μmol/L of each genotyping primer, and 1:10 to 1:20 dilutions of the purified PCR products. The thermal cycling step consisted of the denaturing step at 96°C for 10 seconds, the annealing step at 50°C for 5 seconds, and 25 cycles of the extension step at 60°C for 30 seconds. For post-extension treatment the samples were purified with 1 U SAP. The samples were run with the ABI PRISM 3100, and obtained genotype data were analyzed with GeneScan 2.1 software (Applied Biosystems).

The genotyping success rate varied from 94.8% to 99.1%. All the markers were in Hardy-Weinberg equilibrium, and Mendel errors were found <0.01%.36

Statistical analysis 

Haplotype analysis with Haplotype Pattern Mining (HPM) algorithm37 was performed for dichotomized traits: asthma (affected and control) and high serum total IgE level (serum IgE level>100 kU/L as high responder and serum IgE level ≤ 100 kU/L as low responder), or the combinations of asthma with different serum total IgE levels and associated blood eosinophil accounts. Haplotyping was done within each trio, and 4 independent chromosomes were obtained. In case of ambiguities resulting from missing genotype data, identical heterozygotic genotypes in all of the family members, or mendelian incompatibilities, the alleles were discarded. If the child was affected, the transmitted chromosomes were considered disease-associated and the nontransmitted chromosomes as controls. If one of the parents was affected, his/her chromosomes were considered disease-associated and the spouse's chromosomes as controls. When both the parent and the child were affected, only the nontransmitted chromosome of the unaffected parent was considered as the control and the remaining 3 as disease-associated ones. The input for HPM consists of these haplotypes.37 The maximum length of the haplotype patterns and the number of gaps allowed were given as parameters for the algorithm. For population-based case-control replication data set, the haplotypes were constructed of the data from unrelated individuals by using SNPhap algorithm (http://www-gene.cimr.cam.ac.uk/clayton/software/), and those haplotypes were used in the similar way as input data for HPM. The HPM finds all the phenotype-associated haplotype patterns that exceed a given threshold value by a χ2 test. Each marker is given a marker-wise score based on the number of qualified patterns spanning across the marker. The significance of the observed P value can be tested by permutation where the grouping of chromosomes is done randomly, and the scores are recalculated several times to obtain marker-wise P values that are comparable with each other.37 Quantitative statistical analysis methods were used to evaluate associations between specific haplotype carrier status and elevated total serum IgE levels detected by HPM. For this purpose, total serum IgE measurements were transformed into a normally distributed form by loge transformation and adjusted for age and gender. Resulting new variable was analyzed with independent samples t test and linear regression analysis.

Back to Article Outline

Results 

Identification of polymorphisms 

We screened all the exons and the exon-intron boundaries of STAT6, GATA3, and STAT4 for polymorphisms in 14 to 22 Finnish study subjects. To improve our chances to find functionally important SNPs in IgE-mediated disorders, the pool of screened individuals was enriched with atopic patients with exceptionally high serum total IgE values. All identified heteroduplexes were verified with direct sequencing. Among these individuals, a total of 5 SNPs in STAT6, 9 SNPs in GATA3, and 9 SNPs in STAT4 were identified (Fig 1), of which 8 were not previously reported into the SNP database (dbSNP/NCBI, http://www.ncbi.nlm.nih.gov/SNP/). Two SNPs identified in STAT6 were intronic, and 3 were located in the 3′UTR region of the gene. Nine non-coding variants of GATA3 were observed. One STAT4 polymorphism causing an amino acid substitution (C42113438T, located in the DNA-binding domain of the gene) in exon 16 was identified in one of the study subjects. This SNP was identified as a heteroduplex in DHPLC screening and verified by sequencing.

  • View full-size image.
  • Fig 1. 

    Genomic structure of the (A) STAT6, (B) GATA3, and (C) STAT4 and identified polymorphisms (not in scale). The SNPs are named according to their position in NT_029419 (STAT6), NT_077569 (GATA3), and NT_005403 (STAT4), respectively. Novel SNPs not previously reported in the database are marked with an asterisk.

These 23 SNPs identified were chosen for the haplotype association analysis.

Association analysis 

To analyze the association of 23 SNPs with either asthma or high (>100 kU/L) serum total IgE level, the allele/haplotype association analysis was performed among Finnish nuclear families recruited on the basis of asthma. The data set revealed 199 affected and 161 control chromosomes for asthma; moreover, there were 273 high serum IgE response and 199 low IgE response associated chromosomes.

We found no significant allele association for any of the loci studied (Table I). All the SNPs detected in STAT6 gene were relatively common among the genotyped subjects (28% to 50%). In GATA3 three of 9 SNPs were rare (<5%), whereas in STAT4 all except 1 SNP were rare (<6%). For haplotype analysis we allowed 1 wild card (1 error or missing allele) in the haplotype patterns. For either STAT6 or STAT4, no haplotype association was found (χ2 values in all 2 × 2 contingency tables < 4.9, data not shown). The most common haplotypes explaining the majority of haplotype diversity of both loci are shown in Tables E1 through E3 (in the Journal's Online repository at www.mosby.com/jaci). Interestingly, for GATA3 the haplotype analysis identified 3 potentially associating haplotypes (Table II). One of the haplotypes (CGTAGCTCC) was found to be associated with high total serum IgE levels and was studied further. Subjects were divided into haplotype carriers and noncarriers; their total serum IgE values were transformed into a normally distributed form and analyzed with quantitative analysis methods. Haplotype (CGTAGCTCC) carriers showed significant association to elevated total serum IgE levels in t test (Table III) and in linear regression analysis (P=.01) when age, gender, and asthma status were used as covariates.

Table I. Comparisons of allele frequencies of the SNPs in STAT6, GATA3, and STAT4 in high IgE (273 high IgE response associated vs 161 low response chromosomes) and asthma associated chromosome (199 asthma associated vs 161 control chromosomes), respectively
High IgEAsthma
AlleleAssociated(%)Not associated(%)Associated(%)Not associated(%)
STAT6 specific markers named on basis of their position on NT_029419
19647081G55.865.556.961.3
A44.234.543.138.7
19645488G55.763.356.260.3
A44.336.743.839.7
19633406G77.372.876.074.2
A22.727.224.025.8
19633134A51.851.550.251.2
G48.248.549.848.8
19632954A54.852.457.351.5
G45.247.642.748.5
GATA3 specific markers named on basis of their position on NT_077569
2450477C92.186.290.488.4
T7.913.89.611.6
2458625G99.399.5100.098.8
A0.70.50.01.2
2458701T98.597.598.598.1
G1.52.51.51.9
2459095A21.516.020.919.3
G78.584.079.180.7
2460264G60.659.163.858.3
A39.440.936.241.7
2460466C98.9100.099.099.4
T1.10.01.00.6
2463543C91.691.391.890.5
T8.48.78.29.5
2474305C15.915.517.214.6
T84.184.582.885.4
2478975C89.388.485.991.4
T10.711.614.18.6
STAT4 specific markers named on basis of their position on NT_005403
42225344G98.995.598.596.3
A1.14.51.53.7
42132257C34.933.839.935.0
T65.166.260.165.0
42128770T97.896.596.596.9
C2.23.53.53.1
42113438C100.0100.0100.0100.0
T0.00.00.00.0
42110400T98.998.098.597.6
C1.12.01.52.4
42105770C99.697.599.098.2
T0.42.51.01.8
42105023C97.094.597.494.4
T3.05.52.65.6
42103557T98.195.497.995.6
C1.94.62.14.4
42103821T100.0100.0100.0100.0
C0.00.00.00.0

A rare SNP found in ∼4.7% of screened individuals, but not present among the genotyped families.

Table II. The best associating haplotypes of GATA3 gene detected by HPM in different phenotypic comparisons in both the primary (family collection) and replication (case-control) data sets
Associating haplotypesHaplotype frequencies
245047724586252458701245909524602642460466246354324743052478975I data setII data setPhenotypic comparisons showing associationχ2P value
CGTGGCCTC0.170.11I data set: asthma patients vs controls5.8.2
CGTAGCTCC0.110.08I data set: asthma patients with serum IgE >100 kU/L vs ≤100 kU/L13.1.004
I data set: asthma patients with serum IgE >200 kU/L vs ≤200 kU/L11.7.004
CGTAGCCTC0.090.10I data set: asthma patients with blood eosinophil count >0.44 × 109/L vs ≤0.44 ×109/L3.3.4
II data set: asthma patients with blood eosinophil count >0.2 × 109/L vs ≤0.2 × 109/L8.9.01
II data set: asthma patients with serum IgE >100 IU/L vs ≤100 IU/L4.6.2

The alleles shared between the haplotypes are shown in bold.

I data set, Primary data set; II data set, replication data set.

The best nominal χ2 value for the associating haplotype pattern in 2 × 2 contingency table.

Experimental P value observed by 10,000 permutations.

Table III. Independent 2-sample t test of serum total IgE levels among haplotype CGTAGCTCC carriers in data set I
Associating haplotype Mean total serum IgE (kU/L)
245047724586252458701245909524602642460466246354324743052478975Independent samples t test showing associationCarriersNon-carrierst StatP value
CGTAGCTCCI data set: carriers vs non-carriers2941512.4.02
I data set, asthmatic patients: carriers vs non-carriers5642382.2.04

The alleles shared between the haplotypes are shown in bold.

I data set, Primary data set.

The comparisons with related phenotypes identified partially different GATA3 haplotypes. Six of 9 alleles were shared between the haplotypes. Most of the results remained statistically significant after 10,000 permutations. Obviously we have tested multiple phenotypes, genes, and markers simultaneously. Therefore, the marker-wise P values should be corrected with number of independent tests performed. However, because the phenotypes are closely related and the gene-specific markers are inherited in linkage disequilibrium, it is difficult to determine the appropriate statistical correction. To evaluate the biologic significance of our findings, we decided to replicate the genotyping of GATA3 in another, independent data set. The Finnish case-control data set (245 asthmatic patients and 405 age- and gender-matched control subjects) expresses very similar clinical characteristics (onset of chronic asthma in adulthood and demand for daily medication) as the previously analyzed family collection. The frequencies of the observed haplotypes in the family data set and predicted haplotypes in the case-control data set were comparable. Interestingly, 1 of the 3 haplotypes showing association already in the family data was also observed to associate in the replication data (Table II).

Back to Article Outline

Discussion 

On the basis of chromosomal linkage studies and their function, STAT6 and STAT4 are positional candidates for asthma and atopy.8, 9, 10, 11, 12, 30 GATA3 and T-BET are essential transcription factors in the differentiation of T helper lymphocytes. T-BET was not associated with asthma-related traits on the basis of 15 SNPs genotyped in this same patient cohort used in the present study.38

We found 23 SNPs in the STAT6, GATA3, and STAT4 genes. No coding variants leading to amino acid changes were found in the STAT6 or GATA3 genes, suggesting a highly conserved exon structure of the genes. This is consistent with previous SNP screening studies.13, 14, 16, 39, 40 We concentrated our genotyping efforts into the SNPs (5 to 9 markers per locus) found among the Finnish individuals. Therefore, our study does not include all the SNPs described in the public database in the intronic regions that were not covered by our screening process. The genotyped markers, however, were inherited in linkage disequilibrium (Tables E1 through E3 in the Online Repository), and we should have seen a bias in haplotype distribution among affected individuals also in case the causative SNP is missing. All the SNPs found in STAT6 had been previously reported13, 16 and are relatively common (frequency>22 %) (Table I). In the present study we did not genotype the GT repeat in the first exon of STAT6. This GT repeat has been shown to be associated with elevated eosinophil levels among whites13 and with allergic phenotypes14 and childhood asthma15 among Japanese individuals. SNP in intron 18, shown to associate with total IgE levels (P=.0070) in one study,13 did not behave that way in our study. Genotyping of the 5 SNPs identified among the screened Finnish individuals (28 chromosomes) did not support previous reports.13, 16 There was no association of the 3′UTR SNP (G19633406A) described by Gao et al16 among Finnish families, which, together with previously reported negative associations,13, 14, 15, 16 suggests that the 3′UTR variant does not have an important role in the pathogenesis of asthma and allergy. Thus, the biologic importance of described polymorphisms in STAT6 remains unclear. Apparently both ethnic and clinical differences in study populations may add to the complexity of the results.

GATA3 haplotype analysis in both the family and case-control data sets suggested an association with asthma-related phenotypes (Table II). Three partially different (6/9 alleles shared) associating haplotypes were identified. The allelic differences between these haplotypes may be due to the ethnic background of the 2 data sets, family collection and population-based, analyzed. When studying complex disorders such as asthma and atopy, susceptibility alleles predisposing to disease can differ between founder and outbred populations. Because none of the SNPs of GATA3 altered the amino acid sequence, it may be difficult to directly predict whether the SNPs have a functional role in the disease susceptibility. None of the SNPs seemed to locate in the functional domains of GATA3 (TRANSFAC professional 7.4.1; BIOBASE GmbH, Wolfenbüttel, Germany). Still, it is possible that the SNPs relatively near to intron/exon borders have an effect on mRNA splicing, whereas the SNPs at 5′and 3′UTR sequences affect the efficiency of translation or RNA processing. There are studies indicating that this type of noncoding variants causes changes in the stability and splicing of the encoded mRNAs and, hence, has biologic importance.41, 42, 43, 44, 45, 46, 47, 48 Therefore, the haplotype associations in GATA3 we describe here may be important. The SNPs G2458625A, T2458701G, and A2459095G are all located 1000 to 1500 bp upstream of the coding sequence of GATA3, whereas 2 SNPs in the 5′UTR, G2460264A and C2460466T, are located at position −250 to −50 bp. Being at promoter regions, these may have an unknown effect to transcription. An SNP in the 3′UTR, C2478975T, is located 92 bp downstream of STOP codon. Only 1 SNP of GATA3, C2474305T, is located near the intron/exon border, being 27 bp upstream of the fifth exon. However, further genetic and functional studies are needed for the final conclusions of the functional impact of these SNPs.

In conclusion, we have identified a panel of novel SNPs in genes encoding key proteins that orchestrate immune responses. For STAT6, GATA3, and STAT4 genes this is the first asthma association study among the Finnish asthmatic families. To our knowledge, for STAT4 and GATA3 there are also no previously published association analyses. In addition, a putative association between GATA3 haplotypes and serum high IgE and asthma traits was observed. Elucidation of the importance of this panel of SNPs in a variety of diseases, also other than asthma, characterized by an imbalance between TH1/TH2 responses will be of great interest.

Back to Article Outline

 

We thank Tia Elo, Heidi Jalonen, Tuija Kyrölä, Tuuli Lappalainen, and Eveliina Virtanen for technical assistance and Emmi Ylikoski for comments and discussions. We also acknowledge Dr Johannes Savolainen for collecting patient samples for SNP screening.

Back to Article Outline

Appendix. Supplementary data 

Table E1a-c The most common haplotypes spanning across the STAT6 loci

Back to Article Outline

References 

  1. Shimoda K, van Deursen J, Sangster MY, Sarawar SR, Carson RT, Tripp RA, et al. Lack of IL-4-induced Th2 response and IgE class switching in mice with disrupted Stat6 gene. Nature. 1996;380:630–633
  2. Takeda K, Kamanaka M, Tanaka T, Kishimoto T, Akira S. Impaired IL-13-mediated functions of macrophages in STAT6-deficient mice. J Immunol. 1996;157:3220–3222
  3. Kaplan MH, Schindler U, Smiley ST, Grusby MJ. Stat6 is required for mediating responses to IL-4 and for development of Th2 cells. Immunity. 1996;4:313–319
  4. Akimoto T, Numata F, Tamura M, Takata Y, Higashida N, Takashi T, et al. Abrogation of bronchial eosinophilic inflammation and airway hyperreactivity in signal transducers and activators of transcription (STAT)6-deficient mice. J Exp Med. 1998;187:1537–1542
  5. Kuperman D, Schofield B, Wills-Karp M, Grusby MJ. Signal transducer and activator of transcription factor 6 (Stat6)- deficient mice are protected from antigen-induced airway hyperresponsiveness and mucus production. J Exp Med. 1998;187:939–948
  6. Patel BK, Keck CL, O'Leary RS, Popescu NC, LaRochelle WJ. Localization of the human stat6 gene to chromosome 12q13.3-q14.1, a region implicated in multiple solid tumors. Genomics. 1998;52:192–200
  7. Hoffjan S, Ober C. Present status on the genetic studies of asthma. Curr Opin Immunol. 2002;14:709–717
  8. Barnes KC, Neely JD, Duffy DL, Freidhoff LR, Breazeale DR, Schou C, et al. Linkage of asthma and total serum IgE concentration to markers on chromosome 12q: evidence from Afro-Caribbean and Caucasian populations. Genomics. 1996;37:41–50
  9. CSGA . A genome-wide search for asthma susceptibility loci in ethnically diverse populations: the Collaborative Study on the Genetics of Asthma (CSGA). Nat Genet. 1997;15:389–392
  10. Barnes KC, Freidhoff LR, Nickel R, Chiu YF, Juo SH, Hizawa N, et al. Dense mapping of chromosome 12q13.12-q23.3 and linkage to asthma and atopy. J Allergy Clin Immunol. 1999;104(pt 1):485–491
  11. Nickel R, Wahn U, Hizawa N, Maestri N, Duffy DL, Barnes KC, et al. Evidence for linkage of chromosome 12q15-q24.1 markers to high total serum IgE concentrations in children of the German Multicenter Allergy Study. Genomics. 1997;46:159–162
  12. Wilkinson J, Grimley S, Collins A, Thomas NS, Holgate ST, Morton N. Linkage of asthma to markers on chromosome 12 in a sample of 240 families using quantitative phenotype scores. Genomics. 1998;53:251–259
  13. Duetsch G, Illig T, Loesgen S, Rohde K, Klopp N, Herbon N, et al. STAT6 as an asthma candidate gene: polymorphism-screening, association and haplotype analysis in a Caucasian sib-pair study. Hum Mol Genet. 2002;11:613–621
  14. Tamura K, Arakawa H, Suzuki M, Kobayashi Y, Mochizuki H, Kato M, et al. Novel dinucleotide repeat polymorphism in the first exon of the STAT-6 gene is associated with allergic diseases. Clin Exp Allergy. 2001;31:1509–1514
  15. Shao C, Suzuki Y, Kamada F, Kanno K, Tamari M, Hasegawa K, et al. Linkage and association of childhood asthma with the chromosome 12 genes. J Hum Genet. 2004;49:115–122
  16. Gao PS, Mao XQ, Roberts MH, Arinobu Y, Akaiwa M, Enomoto T, et al. Variants of STAT6 (signal transducer and activator of transcription 6) in atopic asthma (letter). J Med Genet. 2000;37:380–382
  17. Joulin V, Bories D, Eleouet JF, Labastie MC, Chretien S, Mattei MG, et al. A T-cell specific TCR delta DNA binding protein is a member of the human GATA family. Embo J. 1991;10:1809–1816
  18. Ho IC, Vorhees P, Marin N, Oakley BK, Tsai SF, Orkin SH, et al. Human GATA-3: a lineage-restricted transcription factor that regulates the expression of the T cell receptor alpha gene. Embo J. 1991;10:1187–1192
  19. Oosterwegel M, Timmerman J, Leiden J, Clevers H. Expression of GATA-3 during lymphocyte differentiation and mouse embryogenesis. Dev Immunol. 1992;3:1–11
  20. Ting CN, Olson MC, Barton KP, Leiden JM. Transcription factor GATA-3 is required for development of the T-cell lineage. Nature. 1996;384:474–478
  21. Zheng W, Flavell RA. The transcription factor GATA-3 is necessary and sufficient for Th2 cytokine gene expression in CD4 T cells. Cell. 1997;89:587–596
  22. Zhang DH, Cohn L, Ray P, Bottomly K, Ray A. Transcription factor GATA-3 is differentially expressed in murine Th1 and Th2 cells and controls Th2-specific expression of the interleukin-5 gene. J Biol Chem. 1997;272:21597–21603
  23. Zhang DH, Yang L, Cohn L, Parkyn L, Homer R, Ray P, et al. Inhibition of allergic inflammation in a murine model of asthma by expression of a dominant-negative mutant of GATA-3. Immunity. 1999;11:473–482
  24. Ouyang W, Lohning M, Gao Z, Assenmacher M, Ranganath S, Radbruch A, et al. Stat6-independent GATA-3 autoactivation directs IL-4-independent Th2 development and commitment. Immunity. 2000;12:27–37
  25. Asnagli H, Afkarian M, Murphy KM. Cutting edge: identification of an alternative GATA-3 promoter directing tissue-specific gene expression in mouse and human. J Immunol. 2002;168:4268–4271
  26. Kaplan MH, Sun YL, Hoey T, Grusby MJ. Impaired IL-12 responses and enhanced development of Th2 cells in Stat4- deficient mice. Nature. 1996;382:174–177
  27. Ouyang W, Ranganath SH, Weindel K, Bhattacharya D, Murphy TL, Sha WC, et al. Inhibition of Th1 development mediated by GATA-3 through an IL-4- independent mechanism. Immunity. 1998;9:745–755
  28. Murphy KM, Ouyang W, Farrar JD, Yang J, Ranganath S, Asnagli H, et al. Signaling and transcription in T helper development. Annu Rev Immunol. 2000;18:451–494
  29. Yamamoto K, Kobayashi H, Arai A, Miura O, Hirosawa S, Miyasaka N. cDNA cloning, expression and chromosome mapping of the human STAT4 gene: both STAT4 and STAT1 genes are mapped to 2q32.2->q32.3. Cytogenet Cell Genet. 1997;77:207–210
  30. Wjst M, Fischer G, Immervoll T, Jung M, Saar K, Rueschendorf F, et al. A genome-wide search for linkage to asthma. German Asthma Genetics Group. Genomics. 1999;58:1–8
  31. Laitinen T, Kauppi P, Ignatius J, Ruotsalainen T, Daly MJ, Kaariainen H, et al. Genetic control of serum IgE levels and asthma: linkage and linkage disequilibrium studies in an isolated population. Hum Mol Genet. 1997;6:2069–2076
  32. Laitinen T, Ollikainen V, Lazaro C, Kauppi P, de Cid R, Anto JM, et al. Association study of the chromosomal region containing the FCER2 gene suggests it has a regulatory role in atopic disorders. Am J Respir Crit Care Med. 2000;161(pt 1):700–706
  33. Karjalainen J, Hulkkonen J, Pessi T, Huhtala H, Nieminen MM, Aromaa A, et al. The IL1A genotype associates with atopy in nonasthmatic adults. J Allergy Clin Immunol. 2002;110:429–434
  34. Kuklin A, Munson K, Gjerde D, Haefele R, Taylor P. Detection of single-nucleotide polymorphisms with the WAVE DNA fragment analysis system. Genet Test. 1997;1:201–206
  35. Turner D, Choudhury F, Reynard M, Railton D, Navarrete C. Typing of multiple single nucleotide polymorphisms in cytokine and receptor genes using SNaPshot. Hum Immunol. 2002;63:508–513
  36. O'Connell JR, Weeks DE. PedCheck: a program for identification of genotype incompatibilities in linkage analysis. Am J Hum Genet. 1998;63:259–266
  37. Toivonen HT, Onkamo P, Vasko K, Ollikainen V, Sevon P, Mannila H, et al. Data mining applied to linkage disequilibrium mapping. Am J Hum Genet. 2000;67:133–145
  38. Ylikoski E, Kinos R, Sirkkanen N, Pykäläinen M, Savolainen J, Laitinen LA, et al. Association study of 15 novel single nucleotide polymorphisms of the T-bet locus among Finnish asthma families. Clin Exp Allergy. 2004;34:1049–1055
  39. Nagarkatti R, Ghosh B. Identification of single-nucleotide and repeat polymorphisms in two candidate genes, interleukin 4 receptor (IL4RA) and signal transducer and activator of transcription protein 6 (STAT6), for Th2-mediated diseases. J Hum Genet. 2002;47:684–687
  40. Heinzmann A, Grotherr P, Jerkic SP, Lichtenberg A, Braun S, Kruse S, et al. Studies on linkage and association of atopy with the chromosomal region 12q13-24. Clin Exp Allergy. 2000;30:1555–1561
  41. Shen LX, Basilion JP, Stanton VP. Single-nucleotide polymorphisms can cause different structural folds of mRNA. Proc Natl Acad Sci U S A. 1999;96:7871–7876
  42. Gibson AW, Edberg JC, Wu J, Westendorp RG, Huizinga TW, Kimberly RP. Novel single nucleotide polymorphisms in the distal IL-10 promoter affect IL-10 production and enhance the risk of systemic lupus erythematosus. J Immunol. 2001;166:3915–3922
  43. Nickel RG, Casolaro V, Wahn U, Beyer K, Barnes KC, Plunkett BS, et al. Atopic dermatitis is associated with a functional mutation in the promoter of the C-C chemokine RANTES. J Immunol. 2000;164:1612–1616
  44. Rosenwasser LJ, Klemm DJ, Dresback JK, Inamura H, Mascali JJ, Klinnert M, et al. Promoter polymorphisms in the chromosome 5 gene cluster in asthma and atopy. Clin Exp Allergy. 1995;25(suppl 2):74–78discussion 95-6
  45. Von Ahsen N, Oellerich M. The intronic prothrombin 19911A>G polymorphism influences splicing efficiency and modulates effects of the 20210G>A polymorphism on mRNA amount and expression in a stable reporter gene assay system. Blood. 2004;103:586–593
  46. Tokuhiro S, Yamada R, Chang X, Suzuki A, Kochi Y, Sawada T, et al. An intronic SNP in a RUNX1 binding site of SLC22A4, encoding an organic cation transporter, is associated with rheumatoid arthritis. Nat Genet. 2003;35:341–348
  47. Mann V, Hobson EE, Li B, Stewart TL, Grant SFA, Robins SP, et al. A COL1A1 SP1 binding site polymorphism predisposes to osteoporotic fracture by affecting bone density and quality. J Clin Invest. 2001;107:899–907
  48. Fishman D, Faulds G, Jeffery R, Mohamed-Ali V, Yudkin JS, Humpries S, et al. The effect of novel polymorphisms in the interleukin-6 (IL-6) gene on IL-6 transcription and plasma IL-6 levels, and an association with the systemic-onset juvenile chronic arthritis. J Clin Invest. 1998;102:1369–1376

 Supported by grants from the Academy of Finland, National Technology Agency (TEKES), Väinö and Laina Kivi's Foundation, The Finnish Society of Allergology and Immunology, Ida Montin Foundation, and Turku University Hospital Fund.

PII: S0091-6749(04)02625-9

doi:10.1016/j.jaci.2004.10.006

The Journal of Allergy and Clinical Immunology
Volume 115, Issue 1 , Pages 80-87, January 2005