1 Juvenile Diabetes Research Foundation Centre for Prevention of Type 1 Diabetes in Finland, Turku, Finland
2 Department of Virology, University of Turku, Turku, Finland
3 Department of Pediatrics, University of Oulu, Oulu, Finland
4 Department of Pediatrics, University of Pecs, Pecs, Hungary
5 Program of Developmental and Reproductive Biology, Biomedicum, University of Helsinki, Helsinki, Finland
6 Hospital for Children and Adolescents, University of Helsinki, Helsinki, Finland
7 Department of Pediatrics, University of Turku, Turku, Finland
8 Department of Pediatrics, Tampere University Hospital, Tampere, Finland
![]() |
ABSTRACT |
---|
![]() ![]() ![]() ![]() ![]() ![]() |
---|
Type 1 diabetes is caused by immune-mediated destruction of the insulin-producing pancreatic ß-cells that results from a complex interplay between a polygenic background and environmental factors. In addition to the HLA complex, the insulin gene (INS), the cytotoxic T-cellassociated protein 4 gene region (CTLA4), and the protein tyrosine phosphatase gene (PTPN22) are confirmed disease susceptibility loci (1,2). A number of other gene regions have also been found to show evidence of linkage or association with type 1 diabetes. The confirmation of true genetic effects in this complex disease remains a challenge for several reasons. One of the main hurdles is the modest disease predisposing effect of these loci (odds ratio [OR] < 2) that is further complicated by the genetic heterogeneity of the extensive patient cohorts that are required to ensure statistical power. As an example, recently two groups reported association of the small ubiquitin-like modifier 4 gene (SUMO4; OMIM 608829; 6q25) with type 1 diabetes (3,4). However, these findings were not confirmed in two large multiethnic family sets, indicating the complexity of identifying real underlying gene effects for type 1 diabetes (5,6).
The pancreatic ß-cell mass of an individual is tightly regulated according to insulin demand and reflects the balance between the rate of ß-cell replication and neogenesis and the rate of ß-cell apoptosis (7). During an autoimmune attack, the rate of ß-cell death and the capacity of ß-cell regeneration could be crucial factors in determining the progression of the imbalance in glucose homeostasis (8). It is known that a network of transcription factors provide genetic instructions for the differentiation and development of pancreatic ß-cells (9). Any alteration of this sophisticated mechanism could lead to reduced ß-cell mass or function and may even determine whether clinical diabetes develops.
Two recent reports have suggested that the gene encoding the pancreatic transcription factor PAX4 (OMIM 167413; 7q32) is associated with type 1 diabetes (10,11). In the current study, we aimed to assess this association in two populations with high and low incidence rates of type 1 diabetes. These populations also show conspicuous differences in terms of socioeconomic conditions, climatic, and other environmental factors.
![]() |
RESEARCH DESIGN AND METHODS |
---|
![]() ![]() ![]() ![]() ![]() ![]() |
---|
Finnish population series.
Children with type 1 diabetes were derived mainly from the Departments of Pediatrics at the universities in Turku, Helsinki (southern Finland), and Oulu (northern Finland) (n = 309, 50.5% males; mean age at diagnosis 8.4 ± 4.2 years). The control group comprised consecutive healthy infants born in the University Hospitals of Turku and Oulu who were randomly selected for this study (n = 535, 47.5% males).
Hungarian population series.
Patients were consecutively diagnosed in Baranya County (southern Hungary) over a 17-year period from 1980 to 1996 and recruited through the Hungarian Childhood Diabetes Registry (n = 189, 49.7% males; mean age at diagnosis 8.5 ± 4.8 years). Healthy schoolchildren born in the same geographical region and matched for sex and age were used as control subjects (n = 290, 50.7% males; mean age at sampling 10.6 ± 2.7 years).
Four single nucleotide polymorphisms (SNPs) in the PAX4 gene (OMIM 167413, 7q32) were analyzed using DNA sequencing (MegaBace1000 DNA Sequencer; GE Healthcare). Three nonsynonymous SNPs were located in exon 9 [rs712701 (P321H), rs2233583 (V330A), and rs2233584 (Ter333W)], whereas the rs2233578 SNP (R133W) was in exon 3. HLA typing was performed in all Finnish subjects using methods described earlier (12,13). The INS 23 HphI SNP was typed with a similar oligonucleotide hybridization format (T probe 5'-Europium-CTGTCTCCCAGA-3'; A probe 5'-Terbium-CTGTCACCCAGA-3'). The CTLA4 CT60 SNP was genotyped using minisequencing (SnuPe kit; GE Healthcare). Primer sequences are available on request.
Statistical analysis.
For comparisons of genotype frequencies between the groups, 2 statistics with Yates correction were used. P values <0.05 were considered significant. The effect of PAX4 and the other studied genes was also examined by logistic regression analysis (SPSS for Windows, version 11.0.1; SPSS, Chicago, IL). Power calculations were done using standard methods.
![]() |
RESULTS |
---|
![]() ![]() ![]() ![]() ![]() ![]() |
---|
The C/A variant at coding sequence position +1,168 is a nonsynonymous SNP (P321H; rs712701) that was found polymorphic. No differences in the genotype frequencies between cases and control subjects were observed in the Finnish and Hungarian series (Table 1). Accordingly, the comparison on the combined datasets showed no difference between patients and control subjects. The genotype distributions were in Hardy-Weinberg equilibrium both among cases and control subjects in the two series (P = 0.82 and 0.69, respectively, for the Finnish set; P = 0.50 and 0.96, respectively, for the Hungarian set). We stratified the Finnish patients according to age at diagnosis and sex and compared the PAX4 genotype frequencies with control subjects, but no differences were seen (Table 1). Moreover, to elucidate gene-gene interactions, we analyzed the distribution of the PAX4 genotypes in patients stratified according to the HLA DRB1-DQA1-DQB1, insulin gene (23) HphI, and CTLA4 CT60 genotypes and observed no differences (Table 1). The lack of association of PAX4 with type 1 diabetes was confirmed also by logistic regression analysis.
|
![]() |
DISCUSSION |
---|
![]() ![]() ![]() ![]() ![]() ![]() |
---|
In this study, we attempted to reproduce the findings of Biason-Lauber et al. (10) on the association of the gene encoding the PAX4 transcription factor and type 1 diabetes. However, we were unable to confirm a role for the +1,168 polymorphism of PAX4 in disease pathogenesis in the two populations studied with profoundly different socioeconomic conditions, climate, and environment. Similarly, no effect of PAX4 could be seen in any patient subgroups analyzed according to age at diagnosis, sex, HLA, and CTLA4 CT60 genotypes. Since PAX4 is a pancreatic transcriptional regulator, we hypothesized that its effect could be restricted to individuals with specific INS genotypes. We found, however, no evidence for that.
The Finnish population has the highest disease incidence in the world (55.3 · 100,0001 · year1 in 2004; A. Reunanen, personal communication), whereas the Hungarian population is among the low incidence countries in Europe (8.9 · 100,0001 · year1) (17). In addition to this difference, the two populations also show divergence in genetic background, socioeconomic characteristics, and environment. Although the allele frequencies were somewhat different in the two populations, the disease association was consistently lacking in both series. It is unlikely that the findings in our study are explained by interpopulation differences. The statistical power of our study to detect a difference of the magnitude reported by Biason-Lauber et al. (10) at = 0.05 was more than 99%, and it was more than 80% to detect a significant OR of 1.2. To explore the striking difference between our data and those reported by Biason-Lauber et al. (10), we applied the classical
2 goodness-of-fit test for Hardy-Weinberg equilibrium (HWE) on the PAX4+1,168 genotype data from the Swiss and German population. We found that the +1,168 genotype frequencies showed severe departure from the expected values, both among the cases and control subjects (
2 = 20.2, P = 0.00004;
2 = 73.9, P < 0.000001, respectively). Deviation from HWE might be seen in multiple SNPs in a strong linkage disequilibrium block, and it may also be caused by factors such as recent population migration or genetic drift. However, the most common reason for this phenomenon is genotyping error, as shown by the study of Xu et al. (18), who found that as much as 12% of the genotype data on 133 SNPs published in 75 studies were inconsistent with HWE in control subjects. Our comparison of the Swiss and German data on the +1,168 PAX4 genotype frequencies also showed that particularly the distribution of the C/A heterozygous genotype was affected by this bias, being 10.8% lower in cases and 16.2% higher among control subjects compared with the expected values. In large SNP genotyping efforts, indicators of assay quality other than HWE such as genotype call rate or Mendelian error rate may also help to reduce false-positive findings.
Holm et al. (10) observed in the 7q32 region an increase in the logarithm of odds from 2.2 to 5.3 in a pathway restricted linkage analysis using 408 Scandinavian multiplex families. They analyzed three intronic SNPs in the PAX4 gene and observed marginally significant P values in the haplotype transmission analysis. It has to be noted that these do not remain significant after correction for multiple comparisons. The T-C-C PAX4 haplotype, however, showed a P value of 0.001 in the transmission distortion test, but it was extremely rare.
In summary, we suggest that the PAX4+1,168 C/A variant does not play any essential role in genetic susceptibility to type 1 diabetes. The strong coherence between our datasets on two ethnic groups with conspicuously different disease incidences, socioeconomic characteristics, and profoundly diverse environment strengthens this finding. The striking observations made on the German and Swiss population cohorts require replication and a more detailed analysis of this gene to assess whether the difference between the two studies could be explained by heterogeneity in linkage disequilibrium patterns between ethnic groups or might reflect profound variation in the effect of PAX4 on type 1 diabetes susceptibility in particular patient groups.
![]() |
ACKNOWLEDGMENTS |
---|
We thank Eija Nirhamo, Piia Nurmi, Mia Karlson, Ritva Suominen, and Terhi Laakso for their skillful technical assistance.
Address correspondence and reprint requests to Dr. Robert Hermann, Department of Virology, University of Turku, Kiinamyllynkatu 13, FIN-20520 Turku, Finland. E-mail: robert.hermann{at}utu.fi
Received for publication May 2, 2005 and accepted in revised form June 10, 2005
SNP, single nucleotide polymorphism
![]() |
REFERENCES |
---|
![]() ![]() ![]() ![]() ![]() ![]() |
---|
HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Diabetes | Diabetes Care | Clinical Diabetes | Diabetes Spectrum | DOC News |