Significant association between HLA-B*35:01 and onset of drug-induced liver injury caused by Kampo medicines in Japanese patients.

AIM: Drug-induced liver injury (DILI) is a severe and life-threatening immune-mediated adverse effect, occurring rarely among treated patients. We examined genomic biomarkers in the Japanese population that predict the onset of DILI after using a certain class of drugs, such as Kampo products (Japanese traditional medicines). METHODS: A total of 287 patients diagnosed as DILI by hepatology specialists were recruited after written informed consent was obtained. A genome-wide association analysis and human leukocyte antigen (HLA) typing in four digits were performed. RESULTS: We found a significant association (p = 9.41 × 10-10 ) of rs146644517 (G > A) with Kampo product-related DILI. As this polymorphism is located in the HLA region, we evaluated the association of HLA types and found that 12 (63.2%) of 19 Kampo-DILI patients contained HLA-B*35:01, whereas only 15.2% were positive for this HLA among healthy volunteers. The odds ratio was 9.56 (95% confidence interval 3.75-24.46; p = 2.98 × 10-6 , corrected p = 4.17 × 10-5 ), and it increased to 13.55 compared with the DILI patients not exposed to Kampo products. The individual crude drug components in the Kampo products, including Scutellaria root (ougon in Japanese), rhubarb (daiou), Gardenia fruit (sanshishi), and Glycyrrhiza (kanzou), were significantly associated with HLA-B*35:01. CONCLUSIONS: HLA-B*35:01 is a genetic risk factor and a potential predictive biomarker for Kampo-induced DILI in the Japanese population.

Association of an IGHV3-66 gene variant with Kawasaki disease.

In a meta-analysis of three GWAS for susceptibility to Kawasaki disease (KD) conducted in Japan, Korea, and Taiwan and follow-up studies with a total of 11,265 subjects (3428 cases and 7837 controls), a significantly associated SNV in the immunoglobulin heavy variable gene (IGHV) cluster in 14q33.32 was identified (rs4774175; OR = 1.20, P = 6.0 × 10-9). Investigation of nonsynonymous SNVs of the IGHV cluster in 9335 Japanese subjects identified the C allele of rs6423677, located in IGHV3-66, as the most significant reproducible association (OR = 1.25, P = 6.8 × 10-10 in 3603 cases and 5731 controls). We observed highly skewed allelic usage of IGHV3-66, wherein the rs6423677 A allele was nearly abolished in the transcripts in peripheral blood mononuclear cells of both KD patients and healthy adults. Association of the high-expression allele with KD strongly indicates some active roles of B-cells or endogenous immunoglobulins in the disease pathogenesis. Considering that significant association of SNVs in the IGHV region with disease susceptibility was previously known only for rheumatic heart disease (RHD), a complication of acute rheumatic fever (ARF), these observations suggest that common B-cell related mechanisms may mediate the symptomology of KD and ARF as well as RHD.

Association of HLA-DRB1*04:05 allele with drug-induced interstitial lung disease in Japanese population.

Drug-induced interstitial lung disease (DILD) is a life-threatening adverse reaction. The Japanese population is more susceptible to DILD as compared with other populations, suggesting its pathogenesis could vary depending on ethnic genetic background. We conducted case-control studies to elucidate the association between DILD and HLA alleles in the Japanese. The 177 clinically diagnosed DILD patients and 3002 healthy controls for exploration and 55 DILD patients and 201 healthy controls for validation were genotyped for four HLA genes. HLA-DRB1*04:05 was significantly associated with DILD (corrected p = 0.014); this was also validated in the other set of patients/controls. Chemical drugs other than protein therapeutics showed this association (p = 1.7 × 10-4) . The Japanese population showed a higher HLA-DRB1*04:05 frequency than most other populations. In conclusion, HLA-DRB1*04:05 could be associated with DILD susceptibility in Japanese individuals, and its high general frequency may explain the high reported incidence of DILD in Japanese.

Correction to: Investigation of novel variations of ORAI1 gene and their association with Kawasaki disease.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Investigation of novel variations of ORAI1 gene and their association with Kawasaki disease.

ORAI1 encodes a calcium channel essential in the store-operated calcium entry mechanism. A previous genetic association study identified a rare in-frame insertion variant of ORAI1 conferring Kawasaki disease (KD). To deepen our understanding of the involvement of rare variants of ORAI1 in KD pathogenesis, we investigated 3812 patients with KD and 2644 healthy individuals for variations in the protein-coding region of ORAI1. By re-sequencing the study participants' DNA, 27 variants with minor allele frequencies (MAFs) < 0.01 that had not been examined in the previous study were identified. Although no significant association with KD was observed either in single-variant analyses or in a collapsing method analysis of the 27 variants, stratification by MAFs, variant types, and predicted deleteriousness revealed that six rare, deleterious, missense variants (MAF < 0.001, CADD C-score ≥ 20) were exclusively present in KD patients, including three refractory cases (OR = ∞, P = 0.046). The six missense variants include p.Gly98Asp, which has been demonstrated to result in gain of function leading to constitutive Ca2+ entry. Conversely, five types of frameshift variants, all identified near the N terminus and assumed to disrupt ORAI1 function, showed an opposite trend of association (OR = 0.35, P = 0.24). These findings support our hypothesis that genetic variations causing the upregulation of the Ca2+/NFAT pathway confer susceptibility to KD. Our findings also provide insights into the usefulness of stratifying the variants based on their MAFs and on the direction of the effects on protein function when conducting association studies using the gene-based collapsing method.

Case Report: Identification of a CARD8 variant in all three patients with PFAPA syndrome complicated with Kawasaki disease.

Background: Periodic fever, aphthous stomatitis, pharyngitis, and cervical adenitis (PFAPA syndrome), and Kawasaki disease (KD) are both considered to be disorders of the innate immune system, and the potential role of inflammasome activation in the immunopathogenesis of both diseases has been previously described. Case presentation: Herein, we report the clinical courses of three patients who presented a rare combination of PFAPA syndrome and KD. Two patients who presented KD later developed the PFAPA syndrome, of whom one developed recurrent KD 2 years after the initial diagnosis. The third patient developed KD one year after the onset of PFAPA syndrome. The presence of both of these conditions within individual patients, combined with the knowledge that inflammasome activation is involved in both PFAPA syndrome and KD, suggests a shared background of inflammatory dysregulation. To elucidate the mechanism underlying shared inflammatory dysregulation, we investigated the roles of Nod-like receptors (NLRs) and their downstream inflammasome-related genes. All the patients had a frameshift variant in CARD8 (CARD8-FS). A previous study demonstrated a higher frequency of CARD8-FS, whose product loses CARD8 activity and activates the NLRP3 inflammasome, in patients with the PFAPA syndrome. Additionally, the NLRP3 inflammasome is known to be activated in patients with KD. Together, these results suggest that the CARD8-FS variant may also be essential in KD pathogenesis. As such, we analyzed the CARD8 variants among patients with KD. However, we found no difference in the variant frequency between patients with KD and the general Japanese population. Conclusions: We report the clinical courses of three patients with a rare combination of PFAPA syndrome and KD. All the patients had the CARD8-FS variant. However, we could not find a difference in the variant frequency between patients with KD and the general Japanese population. As the frequency of KD is much higher than that of PFAPA among Japanese patients, and the cause of KD is multifactorial, it is possible that only a small portion of patients with KD harbor CARD8-FS as a causative gene.

Association study of matrix metalloproteinase-12 gene polymorphisms and asthma in a Japanese population.

BACKGROUND: Matrix metalloproteinase 12 gene (MMP12) has been shown to be associated with asthma in a Caucasian population. In this study, we investigate whether single-nucleotide polymorphisms (SNPs) of MMP12 are associated with a risk for asthma in a Japanese population. METHODS: We tested for an association between SNPs in MMP12 and asthma, including its severity, in a Japanese population (630 pediatric and 417 adult patients with atopic asthma and 336 children and 632 adults as controls). The rs652438 A and G variants (N357S) were generated by site-directed mutagenesis and an assay with artificial peptide substrates was used to compare two types of MMP12 activity. The effect of MMP12 inhibition with MMP12-specific small interfering RNA (siRNA) on chemokine secretion from airway epithelial cells was also tested in vitro. RESULTS: N357S showed a p value <0.05 for childhood and combined (adult plus childhood) asthma in the dominant model [odds ratio (OR) 1.60, 95% confidence interval (CI) 1.00-2.56, p = 0.047; OR 1.40, 95% CI 1.04-1.89, p = 0.028, respectively]. This risk variant is associated with asthma severity in adult patients. In the functional assay, the minor-allele enzyme showed significantly lower activity than the major-allele enzyme. MMP12-specific siRNA suppressed IP-10 secretion from airway epithelial cells upon stimulation with IFN-ß. CONCLUSIONS: Our results suggest that MMP12 confers susceptibility to asthma and is associated with asthma severity in a Japanese population. MMP12 may be associated with asthma through inappropriate attraction of leukocytes to the inflamed tissue.

Association of the MMP9 gene with childhood cedar pollen sensitization and pollinosis.

Matrix metalloproteinase 9 (MMP9) gene has been shown to be involved in the pathogenesis of allergic rhinitis (AR) and asthma. Previous studies suggested that single-nucleotide polymorphisms (SNPs) of the MMP9 gene conferred a risk for childhood asthma. However, whether the SNPs confer a risk for AR has not been previously investigated. The objective of this study was to investigate whether SNPs of the MMP9 gene are associated with risk of seasonal AR (pollinosis), perennial AR and allergen sensitization. A total of 670 school children were recruited in Japan and genotyped for functional polymorphism in the promoter (-1590C/T: rs3918242) and three amino-acid substitutions (R297Q: rs17576; P574R: rs2250889; R668Q: rs17577). Serum levels of total and specific IgE were determined. Disease status and other clinical characteristics of the subjects were investigated using a questionnaire. Associations between the MMP9 SNPs and both AR and serum IgE levels were evaluated. -1590C/T showed significant association with cedar pollinosis (corrected P (Pcor)=0.039). R668Q was in strong linkage disequilibrium (LD) with -1590C/T and showed significant association with cedar pollinosis (Pcor=0.023) and serum cedar pollen-specific IgE level (Pcor=0.022). A haplotype associated with -1590T and 668Q showed a significant association with cedar pollinosis, orchard grass pollinosis and cedar pollen-specific IgE (Pcor=0.0012, Pcor=0.0059 and Pcor=0.0041, respectively). R297Q and P574R were in weak LD with the rest of the SNPs and did not show significant association with disease. Compared with wild-type MMP9 protein (279R-574P-668R), a variant enzyme (279R-574P-668Q) that showed association with pollinosis had lower activity. However, lower enzyme activity was not associated with disease risk because another variant (279Q-574R-668R) showed lower enzyme activity but was not associated with pollinosis. The -1590T allele and its corresponding haplotype was associated with higher promoter activity and with pollen-specific IgE levels and pollinosis, suggesting that -1590C/T may have more impact on sensitization and disease development than R668Q. Our results suggest that the MMP9 gene confers susceptibility to cedar pollinosis in Japanese children. The MMP9 gene may be associated with pollinosis through sensitization processes.

Markers of Memory CD8 T Cells Depicting the Effect of the BNT162b2 mRNA COVID-19 Vaccine in Japan.

BNT162b2, a nucleoside-modified mRNA vaccine for SARS-CoV-2 spike glycoprotein (S), provides approximately 95% efficacy for preventing COVID-19. However, it remains unclear how effectively memory CD8+ T cells are generated and which genetic and environmental factors affect the generation and function of memory CD8+ T cells elicited by this vaccine. Here, we investigated the frequency and functions of memory CD8+ T cells 3 weeks after the second vaccination in the Japanese population. Using a peptide-MHC pentamer, we detected an increased number of memory CD8+ T cells together with increased serum anti-S protein antibody in females compared with that in males, but the frequency of pentamer-positive cells was not positively correlated with antibody titers. Memory precursor effector cells (KLRG1-CD127+) among both CD8+ cells and pentamer+ cells and effector cells (CD38-HLA-DR+) among pentamer+ cells were more abundant in females than in males. Upon S protein-mediated stimulation of T cells, the intensity of CD107a and granzyme B expression was increased in females compared with that in males, indicating stronger memory CD8+ T cell responses in females than in males. Our studies showed that the BNT162b2 vaccine elicits increased memory CD8+ T cell proliferation and secondary CTL responses in females compared with those in males in the Japanese population. These findings provide an important basis for the distinct sex difference in cellular immune responses to mRNA vaccination and suggest that memory precursor effector cells can be one of markers to evaluate and boost cellular immunity induced by BNT162b2.

Increase in eosinophil-derived neurotoxin level in school children with allergic disease.

Background: Eosinophils are major effector cells of allergic disease and excellent markers of eosinophilic inflammation. Accurate and reliable biomarkers are helpful in the diagnosis, treatment, and control of allergic disease. Objective: This study aimed to investigate an alternate marker of eosinophilic inflammation, eosinophil-derived neurotoxin (EDN), in a number of allergic diseases. Methods: Three hundred ninety-six elementary school-age children with various allergic conditions were recruited for this study. Subgroups included food allergies (FAs), atopic dermatitis (AD), bronchial asthma (BA), and allergic rhinitis (AR). EDN levels in these groups were compared to those in 93 healthy controls (HC). Results: All subjects with allergic disease had elevated levels of serum EDN (median [interquartile range]: FA, 124.2 ng/mL [59.13-160.5 ng/mL]; AD, 110.8 ng/mL [57.52-167.9 ng/mL]; BA, 131.5 ng/mL [60.60-171.0 ng/mL]; AR, 91.32 ng/mL [46.16-145.0 ng/mL]) compared to HC (38.38 ng/mL [32.40-55.62 ng/mL]) (p < 0.0001). These elevated levels were consistent throughout the age range (6-12 years) of the healthy study subjects (p = 0.0679). EDN levels also correlated well with total immunoglobulin E (Rs = 0.5599, p < 0.0001). Looking at all individuals with an allergic disease, the area under the curve was 0.790. Conclusions: Direct measures of eosinophilic inflammation are needed for accurate diagnosis, treatment, and monitoring of allergic diseases. EDN may be a worthy biomarker of eosinophil activity and a useful screening tool for allergic diseases including FA, AD, BA, and AR.

A genome-wide association study identifies RNF213 as the first Moyamoya disease gene.

Moyamoya disease (MMD) shows progressive cerebral angiopathy characterized by bilateral internal carotid artery stenosis and abnormal collateral vessels. Although ∼ 15% of MMD cases are familial, the MMD gene(s) remain unknown. A genome-wide association study of 785,720 single-nucleotide polymorphisms (SNPs) was performed, comparing 72 Japanese MMD patients with 45 Japanese controls and resulting in a strong association of chromosome 17q25-ter with MMD risk. This result was further confirmed by a locus-specific association study using 335 SNPs in the 17q25-ter region. A single haplotype consisting of seven SNPs at the RNF213 locus was tightly associated with MMD (P = 5.3 × 10(-10)). RNF213 encodes a really interesting new gene finger protein with an AAA ATPase domain and is abundantly expressed in spleen and leukocytes. An RNA in situ hybridization analysis of mouse tissues indicated that mature lymphocytes express higher levels of Rnf213 mRNA than their immature counterparts. Mutational analysis of RNF213 revealed a founder mutation, p.R4859K, in 95% of MMD families, 73% of non-familial MMD cases and 1.4% of controls; this mutation greatly increases the risk of MMD (P = 1.2 × 10(-43), odds ratio = 190.8, 95% confidence interval = 71.7-507.9). Three additional missense mutations were identified in the p.R4859K-negative patients. These results indicate that RNF213 is the first identified susceptibility gene for MMD.

Association between Serum Biotin Levels and Cedar Pollinosis in Japanese Schoolchildren.

Biotin is a water-soluble B complex vitamin and coenzyme of five types of carboxylase and plays crucial roles in fatty acid, glucose, and amino acid metabolism. Nutritional biotin deficiency and defective enzymes essential for biotin metabolism cause inflammatory diseases such as eczema-like dermatitis and Crohn's disease; however, little is known about the pathophysiological roles of biotin. This study investigated the relationship between biotin metabolism and human allergic sensitization and diseases by measuring serum levels of biotin, total immunoglobulin E (IgE) and allergen-specific IgEs in more than 400 Japanese schoolchildren aged 6 to 12. The prevalence of allergic diseases, and environmental and life-style factors were also examined by a questionnaire. Like total IgE, serum biotin levels of children showed a log-normal distribution. Meanwhile, Spearman's rank correlation analysis showed weak but significant positive associations between serum biotin levels and total IgE (rho=0.147, p=0.0029) as well as allergen-specific IgEs against egg whites (rho=0.215, p=0.00013), cedar pollen (rho=0.176, p=0.00036), and cat dander (rho=0.130, p=0.0085). Furthermore, mean serum biotin levels in children with cedar pollinosis, but not with other allergic diseases such as asthma and allergic rhinitis, were significantly higher than in those without (p=0.0015). These results suggest a correlation between serum biotin levels and the development of cedar pollinosis. Further prospective studies are needed to evaluate the causal relationship between biotin metabolism and cedar pollen sensitization and pollinosis development.

Replication of genetic association studies in asthma and related phenotypes.

In asthma genetics, the association of highly replicated susceptibility genes lacks consistency across populations. To identify genuine associations, we investigated the reproducibility of the 23 most promising asthma and asthma-related candidate genes in a moderately sized sample from the Japanese population. We compared the frequency of 33 polymorphisms in unrelated cases and controls and tested for their association with asthma, atopy and serum total IgE levels using allele frequency, codominant, dominant and recessive genotype models. On the basis of the consistency of our findings with previous meta-analyses and large replication studies, IL13, TNF, ADAM33, IL4RA and TBXA2R might represent common major asthma and asthma-related trait genes. Individual gene assessment was extended to the interactions between two polymorphisms using our original method. Interactions between TBXA2R and ADAM33, and IL4RA and C3 were suggested to increase the risk for childhood and all asthma (adult and childhood asthma combined). The confirmation of previously reported associations between gene polymorphisms and phenotypes was problematic when as few as several hundred samples per group were used. Stratification of the subjects by environmental factors or other confounding factors may be necessary to improve the sensitivity and reliability of association results.

Single-molecular real-time deep sequencing reveals the dynamics of multi-drug resistant haplotypes and structural variations in the hepatitis C virus genome.

While direct-acting antivirals (DAAs) for hepatitis C virus (HCV) have dramatically progressed, patients still suffer from treatment failures. For the radical eradication of HCV, a deeper understanding of multiple resistance-associated substitutions (RASs) at the single-clone level is essential. To understand HCV quasispecies and their dynamics during DAA treatment, we applied single-molecule real-time (SMRT) deep sequencing on sera from 12 patients with genotype-1b HCV infections with DAA treatment failures, both pre- and post-treatment. We identified >3.2 kbp sequences between NS3 and NS5A genes of 187,539 clones in total, classifying into haplotype codes based on the linkage of seven RAS loci. The number of haplotype codes during the treatment, per sample, significantly decreased from 14.67 ± 9.12 to 6.58 ± 7.1, while the number of nonsynonymous codons on the seven RAS loci, per clone, significantly increased from 1.50 ± 0.92 to 3.64 ± 0.75. In five cases, the minority multi-drug resistant haplotypes at pre-treatment were identical to the major haplotypes at relapse. Moreover, various structural variations (SVs) were detected and their dynamics analysed. These results suggest that SMRT deep sequencing is useful for detecting minority haplotypes and SVs, and to evaluate the dynamics of viral genomes at the single-clone level.

Association of TNFRSF4 gene polymorphisms with essential hypertension.

BACKGROUND: Essential hypertension is a complex disorder that results from the interaction of a number of susceptibility genes and environmental factors. The TNFRSF4 (tumor necrosis factor receptor superfamily, member 4) gene was one of the genes that showed altered renal expression in long-term salt loading in mice. Moreover, association of the TNFRSF4 and TNFSF4 (tumor necrosis factor (ligand) superfamily, member 4) genes with myocardial infarction was recently reported. Since essential hypertension is a well-known risk factor for myocardial infarction, we hypothesized that TNFRSF4 could be a susceptibility gene for essential hypertension. METHODS: We performed a case-control study of TNFRSF4 in two independent population. RESULTS: Extensive investigation of single nucleotide polymorphisms of the entire gene suggested that it resided in one linkage disequilibrium block, and four single nucleotide polymorphisms in the 5' flanking region sufficiently represented major haplotypes. In the combined population, the frequency of the most frequent haplotype, C-C-A-A, was significantly lower (P = 8.07 x 10(-5)) and that of the second most frequent haplotype, C-T-G-A, was significantly higher (P = 6.07 x 10(-4)) in hypertensive subjects than in control subjects. This difference was observed only in female patients. The C-T-G-A haplotype showed a lower promoter activity than other haplotypes, suggesting a relationship with disease susceptibility. CONCLUSION: Our results suggest that TNFRSF4 is a female-specific susceptible gene for essential hypertension.

IL-10 gene polymorphism, but not TGF-beta1 gene polymorphisms, is associated with food allergy in a Japanese population.

The regulatory IL-10 and TGF-beta1 cytokine gene polymorphisms have been associated with allergic diseases in different populations, like Caucasian, Chinese and Indians. We investigated the association between the polymorphisms IL-10 A-1082G, C-819T, C-627A and TGF-beta1 T+869C, G+915C, C-509T and food allergy in Japanese children. One hundred and eleven children with food allergy and 115 atopic control children without food allergy were recruited. DNA samples from these subjects were genotyped by using PCR. The odds ratio of IL-10 -1082 AA genotype was 2.5 (95% CI, 1.0-6.4) for food allergy risk when compared with atopic control subjects (p = 0.03). There were no significative differences in the frequency of TGF-beta1 gene polymorphisms between both groups. Our results indicate that IL-10 A-1082G gene polymorphism is associated with food allergy susceptibility in atopic Japanese children.

Identification of functional polymorphisms in the promoter region of the human PICK1 gene and their association with methamphetamine psychosis.

OBJECTIVE: Protein interacting with C-kinase-1 (PICK1) plays a role in the targeting and clustering of dopamine transporter, which is the primary target site for the abused drug methamphetamine. Based on the interaction of PICK1 with dopamine transporter, it is of particular interest to investigate the association between the PICK1 gene and methamphetamine abusers. METHOD: The authors studied the association between PICK1 gene polymorphisms and methamphetamine abusers in a Japanese group. Two hundred and eight methamphetamine abusers and 218 healthy comparison subjects were enrolled in the study. Furthermore, the authors also examined the effects of single nucleotide polymorphisms (SNPs) in the promoter and 5'-untranslated region on transcription levels of PICK1. RESULTS: The authors identified four highly frequent SNPs, rs737622 (-332 C/G) and rs3026682 (-205 G/A) in the promoter region and rs713729 (T/A) in intron3 and rs2076369 (T/G) in intron4. Of these SNPs, rs713729 was significantly associated with methamphetamine abusers in general, and rs713729 and rs2076369 were significantly associated with those with spontaneous relapse of psychosis. Furthermore, haplotype analysis revealed that specific haplotypes of these SNPs were associated with methamphetamine abusers. A gene reporter assay revealed that the two SNPs in the promoter region significantly altered transcriptional activity. CONCLUSIONS: Our findings suggest that the PICK1 gene may be implicated in the susceptibility to spontaneous relapse of methamphetamine psychosis and that, as an intracellular adapter protein, PICK1 may play a role in the pathophysiology of methamphetamine psychosis.

The GSTP1 gene is a susceptibility gene for childhood asthma and the GSTM1 gene is a modifier of the GSTP1 gene.

BACKGROUND: Bronchial asthma is a chronic airway disorder characterized by bronchial inflammation. Oxidative stress is a key component of inflammation. Glutathione S-transferase P1 (GSTP1), the abundant isoform of glutathione S-transferases (GSTs) in lung epithelium, plays a key role in cellular protection against oxidative stress. Several studies have shown that the GSTP1 geneis involved in the pathogenesis of asthma and a gene-gene interaction may occur within the GST gene superfamily. METHODS: We screened single-nucleotide polymorphisms (SNPs) at the GSTP1 locus and performed an association study in the Japanese population using two independent case-control groups (group 1: 391 pediatric patients with asthma, 462 adult patients with asthma, and 639 controls, and group 2: 115 pediatric patients with asthma and 184 controls). The effect of GSTM1 null/present genotype on the association between GSTP1 Ile105Val and asthma was also investigated. RESULTS: We identified 20 SNPs at this locus and found this region consisted of one linkage disequilibrium block represented by four SNPs (tag SNPs). The association between the Ile105Val polymorphism in the GSTP1 gene and childhood asthma was significant in both groups (p = 0.047 in group 1, and p = 0.021 in group 2). This association was only significant in patients with GSTM1-positive genotype in both groups (group 1: GSTM1 present p = 0.013 and GSTM1 null p = 0.925, and group 2: GSTM1 present p = 0.015 and GSTM1 null p = 0.362). CONCLUSIONS: These findings suggest that the GSTP1 gene is a childhood asthma susceptible gene, and the GSTM1 gene is a modifier gene of GSTP1 for the risk of childhood asthma in the Japanese population.