Association of Polymorphisms and Haplotypes in the Insulin-Like Growth Factor 1 Receptor (IGF1R) Gene with the Risk of Breast Cancer in Korean Women

The insulin-like growth factor (IGF) signaling pathway plays an important role in cancer biology. The IGF 1 receptor (IGF1R) overexpression has been associated with a number of hematological neoplasias and solid tumors including breast cancer. However, molecular mechanism involving IGF1R in carcinogenic developments is clearly not known. We investigated the genetic variations across the IGF1R polymorphism and the risk of breast cancer risk in Korean women. A total of 1418 individuals comprising 1026 breast cancer cases and 392 age-matched controls of Korean were included for the analysis. Genomic DNA was extracted from whole blood and single nucleotide polymorphisms (SNPs) were analyzed on the GoldenGate Assay system by Illumina’s Custom Genetic Analysis service. SNPs were selected for linkage disequilibrium (LD) analysis by Haploview. We genotyped total 51 SNPs in the IGF1R gene and examined for association with breast cancer. All the SNPs investigated were in Hardy-Weinberg equilibrium. These SNPs tested were significantly associated with breast cancer risk, after correction for multiple comparisons by adjusting for age at diagnosis, BMI, age at menarche, and age at first parturition. Among 51 IGF1R SNPs, five intron located SNPs (rs8032477, rs7175052, rs12439557, rs11635251 and rs12916884) with homozygous genotype (variant genotype) were associated with decreased risk of breast cancer. Fisher’s combined p-value for the five SNPs was 0.00032. Three intron located SNPs with heterozygous genotypes also had decreased risk of breast cancer. Seven of the 51 IGF1R SNPs were in LD and in one haplotype block, and were likely to be associated with breast cancer risk. Overall, this case-control study demonstrates statistically significant associations between breast cancer risk and polymorphisms in IGF1R gene.     

SNP identification, linkage disequilibrium, and haplotype analysis for a 200-kb genomic region in a Korean population

Understanding patterns of linkage disequilibrium (LD) across genomes may facilitate association mapping studies to localize genetic variants influencing complex diseases, a recognition that led to the International Haplotype Mapping Project (HapMap). Divergent patterns of haplotype frequency and LD across global populations require that the HapMap database be supplemented with haplotype and LD data from additional populations. We conducted a pilot study of the LD and haplotype structure of a genomic region in a Korean population. A total of 165 SNPs were identified in a 200-kb region of 22q13.2 by direct sequencing. Unphased genotype data were generated for 76 SNPs in 90 unrelated Korean individuals. LD, haplotype diversity, and recombination rates were assessed in this region and compared with the HapMap database. The pattern of LD and haplotype frequencies of Korean samples showed a high degree of similarity with Japanese data. There was a strong correlation between high LD and low recombination frequency in this region. We found considerable similarities in local LD patterns between three Asian populations (Han Chinese, Japanese, and Korean) and the CEPH population. Haplotype frequencies were, however, significantly different between them. Our results should further the understanding of distinctive Korean genomic features and assist in designing appropriate association studies.     

Interpopulation linkage disequilibrium patterns of GABRA2 and GABRG1 genes at the GABA cluster locus on human chromosome 4

GABRA2 and GABRG1, which encode the alpha-2 and gamma-1 subunits, respectively, of the GABA(A) receptor, are located in a cluster on chromosome 4p. The GABRA2 locus has been found to be associated with alcohol dependence in several studies, but no functional variant that can account for this association has been identified. To understand the reported associations, we sought to understand the linkage disequilibrium (LD) patterns and haplotype structures of these genes. With close intergenic distance, approximately 90 kb, it was anticipated that some markers might show intergenic LD. Variation in 13-SNP haplotype block structure was observed in five different populations: European American, African American, Chinese (Han and Thai), Thai, and Hmong. In the Hmong, a 280-kb region of considerably higher LD spans the intergenic region, whereas in other populations, there were two or more LD blocks that cross this region. These findings may aid in understanding the genetic association of this locus with alcohol dependence in several populations.     

Comparative study of the linkage disequilibrium of an ENCODE region, chromosome 7p15, in Korean, Japanese, and Han Chinese samples

The extent and pattern of linkage disequilibrium (LD) in the human genome provide important information for disease gene mapping. Previous studies have shown that LDs vary depending on chromosomal regions and populations. As the Asian samples of the International HapMap Project consisted of Japanese and Chinese populations, it was of interest whether we could use the HapMap data as a reference to carry out association studies of common complex diseases in a closely related population, such as Koreans. We have compared the LD and recombination patterns defined by single-nucleotide polymorphisms (SNPs) in ENCODE region ENm010, chromosome 7p15.2, in Korean, Japanese, and Chinese samples and further tested the robustness of tagSNPs among the Asian samples. We genotyped 792 SNPs in 500 kb (chromosome 7: 26699793-27199792, NCBI build 34) from 90 unrelated Koreans by fluorescence polarization detection and compared the data with Asian data from the HapMap project. Despite some differences in the position of high LD region boundaries, the overall patterns of LD were remarkably similar across the three samples, reflecting strong genetic affinities among them. Furthermore, the haplotype tag SNP transferability across the three samples was greater than 90%. Our results support the initial suggestion that the populations genotyped in the HapMap project might serve as reference populations for the selection of tagSNPs in association studies.     

与宁夏人群强直性脊柱炎关联的新基因淋巴毒素-α(LTA)的识别

淋巴毒素-α(Lymphotoxin-alpha,LTA)基因与系统性红斑狼疮、银屑病及类风湿性关节炎的遗传性有关,但目前还未有关于LTA基因与强直性脊柱炎(Ankylosing spondylitis,AS)关联的报道。文章采用病例-对照设计,在宁夏人群中对人类白细胞抗原(Human leukocyte antigen,HLA)的Ⅲ类基因约58 kb区域进行了高密度标志的基因组扫描,在33个SNPs及其单倍型中,仅定位于LTA基因中的SNPs组成的TCC单倍型的分布在病例-对照间比较有统计学意义(P=0.0005)。在宁夏群体(病例组:300,对照组:385)中发现,LTA基因中的rs909253 T/C多态性在AS患者中出现的频率显著高于正常对照(28.5%vs 19.7%,P=2×10-6)。结果表明LTA基因变异和AS易感性之间存在相关性,由此识别LTA基因可能与宁夏人群AS关联。     

The extent of linkage disequilibrium in rice (Oryza sativa L.)

Despite its status as one of the world's major crops, linkage disequilibrium (LD) patterns have not been systematically characterized across the genome of Asian rice (Oryza sativa). Such information is critical to fully exploit the genome sequence for mapping complex traits using association techniques. Here we characterize LD in five 500-kb regions of the rice genome in three major cultivated rice varieties (indica, tropical japonica, and temperate japonica) and in the wild ancestor of Asian rice, Oryza rufipogon. Using unlinked SNPs to determine the amount of background linkage disequilibrium in each population, we find that the extent of LD is greatest in temperate japonica (probably >500 kb), followed by tropical japonica (approximately 150 kb) and indica (approximately 75 kb). LD extends over a shorter distance in O. rufipogon (<40 kb) than in any of the O. sativa groups assayed here. The differences in the extent of LD among these groups are consistent with differences in outcrossing and recombination rate estimates. As well as heterogeneity between groups, our results suggest variation in LD patterns among genomic regions. We demonstrate the feasibility of genomewide association mapping in cultivated Asian rice using a modest number of SNPs.     

Association of the ABCA1 gene polymorphisms with type 2 DM in a Japanese population

To examine the association of the ATP-binding cassette transporter 1 (ABCA1) gene with type 2 diabetes (DM), we studied genetic polymorphisms of the ABCA1 gene including its linkage disequilibrium (LD) and haplotype analyses using a Japanese population. A sample set (DM:72, IGT:75, and NGT:227) was genotyped with 34 SNPs distributed from the promoter region to the last exon of the ABCA1 gene. LD between SNPs was assessed in pairwise manner. Among 13 LD blocks constructed, an LD block at the 5'-region showed a significant difference in the haplotype distribution between the study groups (NGT vs. IGT + DM: overall p = 0.0180; NGT vs. DM: 0.0001). Fisher's exact probability test (NGT vs. DM) showed a significant association of the haplotype 2 of the LD block (p = 0.0001), with an odds ratio (OR) of 2.53 (95%CI:1.62-4.12). Diplotype analysis also showed a significant association of the diplotypes with the haplotype 2 (OR:2.59, 95%CI:1.48-4.54, p = 0.0013).     

Linkage disequilibrium mapping of Arabidopsis CRY2 flowering time alleles

The selfing plant Arabidopsis thaliana has been proposed to be well suited for linkage disequilibrium (LD) mapping as a means of identifying genes underlying natural trait variation. Here we apply LD mapping to examine haplotype variation in the genomic region of the photoperiod receptor CRYPTOCHROME2 and associated flowering time variation. CRY2 DNA sequences reveal strong LD and the existence of two highly differentiated haplogroups (A and B) across the gene; in addition, a haplotype possessing a radical glutamine-to-serine replacement (AS) occurs within the more common haplogroup. Growth chamber and field experiments using an unstratified population of 95 ecotypes indicate that under short-day photoperiod, the AS and B haplogroups are both highly significantly associated with early flowering. Data from six genes flanking CRY2 indicate that these haplogroups are limited to an approximately 65-kb genomic region around CRY2. Whereas the B haplogroup cannot be delimited to <16 kb around CRY2, the AS haplogroup is characterized almost exclusively by the nucleotide polymorphisms directly associated with the serine replacement in CRY2; this finding strongly suggests that the serine substitution is directly responsible for the AS early flowering phenotype. This study demonstrates the utility of LD mapping for elucidating the genetic basis of natural, ecologically relevant variation in Arabidopsis.     

Family-based analysis using a dense single-nucleotide polymorphism-based map defines genetic variation at PSORS1, the major psoriasis-susceptibility locus

Psoriasis is a common skin disorder of multifactorial origin. Genomewide scans for disease susceptibility have repeatedly demonstrated the existence of a major locus, PSORS1 (psoriasis susceptibility 1), contained within the major histocompatibility complex (MHC), on chromosome 6p21. Subsequent refinement studies have highlighted linkage disequilibrium (LD) with psoriasis, along a 150-kb segment that includes at least three candidate genes (encoding human leukocyte antigen-C [HLA-C], alpha-helix-coiled-coil-rod homologue, and corneodesmosin), each of which has been shown to harbor disease-associated alleles. However, the boundaries of the minimal PSORS1 region remain poorly defined. Moreover, interpretations of allelic association with psoriasis are compounded by limited insight of LD conservation within MHC class I interval. To address these issues, we have pursued a high-resolution genetic characterization of the PSORS1 locus. We resequenced genomic segments along a 220-kb region at chromosome 6p21 and identified a total of 119 high-frequency SNPs. Using 59 SNPs (18 coding and 41 noncoding SNPs) whose position was representative of the overall marker distribution, we genotyped a data set of 171 independently ascertained parent-affected offspring trios. Family-based association analysis of this cohort highlighted two SNPs (n.7 and n.9) respectively lying 7 and 4 kb proximal to HLA-C. These markers generated highly significant evidence of disease association (P<10-9), several orders of magnitude greater than the observed significance displayed by any other SNP that has previously been associated with disease susceptibility. This observation was replicated in a Gujarati Indian case/control data set. Haplotype-based analysis detected overtransmission of a cluster of chromosomes, which probably originated by ancestral mutation of a common disease-bearing haplotype. The only markers exclusive to the overtransmitted chromosomes are SNPs n.7 and n.9, which define a 10-kb PSORS1 core risk haplotype. These data demonstrate the power of SNP haplotype-based association analyses and provide high-resolution dissection of genetic variation across the PSORS1 interval, the major susceptibility locus for psoriasis.     

Genetic determinants of hair and eye colours in the Scottish and Danish populations

Background

The Bovine HapMap Consortium has generated assay panels to genotype ~30,000 single nucleotide polymorphisms (SNPs) from 501 animals sampled from 19 worldwide taurine and indicine breeds, plus two outgroup species (Anoa and Water Buffalo). Within the larger set of SNPs we targeted 101 high density regions spanning up to 7.6 Mb with an average density of approximately one SNP per 4 kb, and characterized the linkage disequilibrium (LD) and haplotype block structure within individual breeds and groups of breeds in relation to their geographic origin and use.

Results

From the 101 targeted high-density regions on bovine chromosomes 6, 14, and 25, between 57 and 95% of the SNPs were informative in the individual breeds. The regions of high LD extend up to ~100 kb and the size of haplotype blocks ranges between 30 bases and 75 kb (10.3 kb average). On the scale from 1–100 kb the extent of LD and haplotype block structure in cattle has high similarity to humans. The estimation of effective population sizes over the previous 10,000 generations conforms to two main events in cattle history: the initiation of cattle domestication (~12,000 years ago), and the intensification of population isolation and current population bottleneck that breeds have experienced worldwide within the last ~700 years. Haplotype block density correlation, block boundary discordances, and haplotype sharing analyses were consistent in revealing unexpected similarities between some beef and dairy breeds, making them non-differentiable. Clustering techniques permitted grouping of breeds into different clades given their similarities and dissimilarities in genetic structure.

Conclusion

This work presents the first high-resolution analysis of haplotype block structure in worldwide cattle samples. Several novel results were obtained. First, cattle and human share a high similarity in LD and haplotype block structure on the scale of 1–100 kb. Second, unexpected similarities in haplotype block structure between dairy and beef breeds make them non-differentiable. Finally, our findings suggest that ~30,000 uniformly distributed SNPs would be necessary to construct a complete genome LD map in Bos taurus breeds, and ~580,000 SNPs would be necessary to characterize the haplotype block structure across the complete cattle genome.     

Identification of functional genetic variations underlying drought tolerance in maize using SNP markers

Single nucleotide polymorphism (SNP) is a common form of genetic variation and popularly exists in maize genome. An Illumina GoldenGate assay with 1 536 SNP markers was used to genotype maize inbred lines and identified the functional genetic variations underlying drought tolerance by association analysis. Across 80 lines, 1 006 polymorphic SNPs (65.5% of the total) in the assay with good call quality were used to estimate the pattern of genetic diversity, population structure, and familial relatedness. The analysis showed the best number of fixed subgroups was six, which was consistent with their original sources and results using only simple sequence repeat markers. Pairwise linkage disequilibrium (LD) and association mapping with phenotypic traits investigated under water-stressed and well-watered regimes showed rapid LD decline within 100-500 kb along the physical distance of each chromosome, and that 29 SNPs were associated with at least two phenotypic traits in one or more environments, which were related to drought-tolerant or drought-responsive genes. These drought-tolerant SNPs could be converted into functional markers and then used for maize improvement by marker-assisted selection.     

Single nucleotide polymorphisms and the linkage disequilibrium at the LDL receptor gene in Koreans

The low-density lipoprotein (LDL) receptor gene is one of the important genes in cardiovascular biology, but gives difficulty in single nucleotide polymorphism (SNP) marker choice. Twelve SNPs at the LDL receptor locus were genotyped in 117 Korean individuals. Six SNPs (rs8111982, rs12983082, rs1003723, rs5925, rs6413504, and rs2738464) were polymorphic in the Korean sample. And 8 SNPs (rs5931, rs811 1982, rs12983082, rs11669576, rs1003723, rs5925, rs641 3504, and rs2738464) were polymorphic in at least one population among the five ethnic groups (European: CEU; Chinese: CHB; Japanese: JPT; African:YRI; and Korean: KOR). Phylogenetic trees constructed based on the twelve SNPs showed a pattern clustering KOR with CHB and JPT as more closely related than either YRI or CEU. The plot of the pairwise linkage disequilibrium (LD) values for the exon 13 SNP (rs5925) against the physical distances between pairs of SNPs showed a tendency of LD decrease in the upstream of rs5925, while the low r² values are present in the downstream of rs5925 (≥0.20). The three SNPs (rs12983082, rs1003723, and rs5925) which are ethnically polymorphic are likely to be useful SNP markers for the LDL receptor locus. The interblock region between exon 13 and exon 18 delimited by the less stringent block definition (four gamete rule) contains a feasible recombination hotspot in the LDL receptor gene in Koreans.     

Genetic Association Study of Adiposity and Melanocortin-4 Receptor (MC4R) Common Variants: Replication and Functional Characterization of Non-Coding Regions

Common genetic variants 3′ of MC4R within two large linkage disequilibrium (LD) blocks spanning 288 kb have been associated with common and rare forms of obesity. This large association region has not been refined and the relevant DNA segments within the association region have not been identified. In this study, we investigated whether common variants in the MC4R gene region were associated with adiposity-related traits in a biracial population-based study. Single nucleotide polymorphisms (SNPs) in the MC4R region were genotyped with a custom array and a genome-wide array and associations between SNPs and five adiposity-related traits were determined using race-stratified linear regression. Previously reported associations between lower BMI and the minor alleles of rs2229616/Val103Ile and rs52820871/Ile251Leu were replicated in white female participants. Among white participants, rs11152221 in a proximal 3′ LD block (closer to MC4R) was significantly associated with multiple adiposity traits, but SNPs in a distal 3′ LD block (farther from MC4R) were not. In a case-control study of severe obesity, rs11152221 was significantly associated. The association results directed our follow-up studies to the proximal LD block downstream of MC4R. By considering nucleotide conservation, the significance of association, and proximity to the MC4R gene, we identified a candidate MC4R regulatory region. This candidate region was sequenced in 20 individuals from a study of severe obesity in an attempt to identify additional variants, and the candidate region was tested for enhancer activity using in vivo enhancer assays in zebrafish and mice. Novel variants were not identified by sequencing and the candidate region did not drive reporter gene expression in zebrafish or mice. The identification of a putative insulator in this region could help to explain the challenges faced in this study and others to link SNPs associated with adiposity to altered MC4R expression.     

Genome-wide association studies demonstrate that TASP1 contributes to increased muscle fiber diameter

Muscle fiber diameter is an economically important trait because it affects meat yield and quality. However, the genetic basis underlying muscle fiber diameter has not been determined. In this study, we collected THREE muscular histological phenotypes in 479 ducks from an F₂ segregating population generated by mallard × Pekin duck crosses. We performed genome-wide association studies (GWAS) and identified a quantitative trait locus (QTL) significantly associated with muscle fiber diameter on chromosome 3. Then, we discovered the selection signatures using the fixation index among 40 mallards and 30 Pekin ducks in this QTL region. Furthermore, we characterized the recombination event in this QTL region and identified a 6-kb block located on TASP1 that was significantly associated with muscle fiber diameter. Finally, five SNPs were screened as potential causative mutations within the 6-kb block. In conclusion, we demonstrated that TASP1 contributes to an increase in muscle fiber diameter, which helps to characterize muscle development and contributes to the genetic improvement of meat yield and quality in livestock.Subject terms: Genome-wide association studies, Gene expression     

Prenatal famine and genetic variation are independently and additively associated with DNA methylation at regulatory loci within IGF2/H19

Both the early environment and genetic variation may affect DNA methylation, which is one of the major molecular marks of the epigenome. The combined effect of these factors on a well-defined locus has not been studied to date. We evaluated the association of periconceptional exposure to the Dutch Famine of 1944-45, as an example of an early environmental exposure, and single nucleotide polymorphisms covering the genetic variation (tagging SNPs) with DNA methylation at the imprinted IGF2/H19 region, a model for an epigenetically regulated genomic region. DNA methylation was measured at five differentially methylated regions (DMRs) that regulate the imprinted status of the IGF2/H19 region. Small but consistent differences in DNA methylation were observed comparing 60 individuals with periconceptional famine exposure with unexposed same-sex siblings at all IGF2 DMRs (P(BH)<0.05 after adjustment for multiple testing), but not at the H19 DMR. IGF2 DMR0 methylation was associated with IGF2 SNP rs2239681 (P(BH) = 0.027) and INS promoter methylation with INS SNPs, including rs689, which tags the INS VNTR, suggesting a mechanism for the reported effect of the VNTR on INS expression (P(BH) = 3.4 × 10(-3)). Prenatal famine and genetic variation showed similar associations with IGF2/H19 methylation and their contributions were additive. They were small in absolute terms (<3%), but on average 0.5 standard deviations relative to the variation in the population. Our analyses suggest that environmental and genetic factors could have independent and additive similarly sized effects on DNA methylation at the same regulatory site.     

Mucin variable number tandem repeat polymorphisms and severity of cystic fibrosis lung disease: significant association with MUC5AC

Variability in cystic fibrosis (CF) lung disease is partially due to non-CFTR genetic modifiers. Mucin genes are very polymorphic, and mucins play a key role in the pathogenesis of CF lung disease; therefore, mucin genes are strong candidates as genetic modifiers. DNA from CF patients recruited for extremes of lung phenotype was analyzed by Southern blot or PCR to define variable number tandem repeat (VNTR) length polymorphisms for MUC1, MUC2, MUC5AC, and MUC7. VNTR length polymorphisms were tested for association with lung disease severity and for linkage disequilibrium (LD) with flanking single nucleotide polymorphisms (SNPs). No strong associations were found for MUC1, MUC2, or MUC7. A significant association was found between the overall distribution of MUC5AC VNTR length and CF lung disease severity (p = 0.025; n = 468 patients); plus, there was robust association of the specific 6.4 kb HinfI VNTR fragment with severity of lung disease (p = 6.2×10(-4) after Bonferroni correction). There was strong LD between MUC5AC VNTR length modes and flanking SNPs. The severity-associated 6.4 kb VNTR allele of MUC5AC was confirmed to be genetically distinct from the 6.3 kb allele, as it showed significantly stronger association with nearby SNPs. These data provide detailed respiratory mucin gene VNTR allele distributions in CF patients. Our data also show a novel link between the MUC5AC 6.4 kb VNTR allele and severity of CF lung disease. The LD pattern with surrounding SNPs suggests that the 6.4 kb allele contains, or is linked to, important functional genetic variation.     

Association of SERPINE2 gene with the risk of chronic obstructive pulmonary disease and spirometric phenotypes in northern Han Chinese population

Chronic obstructive pulmonary disease (COPD) is a complex human disease influenced by multiple genes and environmental factors. The SERPINE2 gene has recently been demonstrated to be associated with COPD onset in a non-East Asian population. In this study, we genotyped 20 single nucleotide polymorphisms (SNPs) in SERPINE2 from 310 cases and 203 controls, all of which belong to the Han from North China. Genotype frequencies were compared between the cases and the controls and analyzed for statistical significance. Two SNPs (rs729631 and rs975278), which are in strong linkage disequilibrium (LD) and locate in block 1 on the LD map of our samples, showed significant association both with the risk of COPD and decline in baseline lung function after Bonferroni correction (P < 0.05). This study provides further evidences for SERPINE2 gene as a COPD susceptible gene, and block 1 of SERPINE2 appears to be the genetic variant region that affects the Han Chinese.     

The variable number of tandem repeats element in DAT1 regulates in vitro dopamine transporter density

Background

The selection of markers in association studies can be informed through the use of haplotype blocks. Recent reports have determined the genomic architecture of chromosomal segments through different haplotype block definitions based on linkage disequilibrium (LD) measures or haplotype diversity criteria. The relative applicability of distinct block definitions to association studies, however, remains unclear. We compared different block definitions in 6.1 Mb of chromosome 17q in 189 unrelated healthy individuals. Using 137 single nucleotide polymorphisms (SNPs), at a median spacing of 15.5 kb, we constructed haplotype block maps using published methods and additional methods we have developed. Haplotype tagging SNPs (htSNPs) were identified for each map.

Results

Blocks were found to be shorter and coverage of the region limited with methods based on LD measures, compared to the method based on haplotype diversity. Although the distribution of blocks was highly variable, the number of SNPs that needed to be typed in order to capture the maximum number of haplotypes was consistent.

Conclusion

For the marker spacing used in this study, choice of block definition is not important when used as an initial screen of the region to identify htSNPs. However, choice of block definition has consequences for the downstream interpretation of association study results.     

Association of the programmed cell death 1 (PDCD1) gene polymorphism with ankylosing spondylitis in the Korean population

The PD-1 (programmed death 1) molecule is a negative regulator of T cells. PDCD1 (programmed cell death 1) has been reported to have a genetic association in systemic lupus erythematosus and rheumatoid arthritis in Caucasians. However, there are no reports on the association between this gene and ankylosing spondylitis (AS). The present study investigated the association of the PD-1 polymorphisms and the haplotypes with AS in a Korean population sample. In a case-control association study, two single-nucleotide polymorphisms, PD-1.5 C/T and PD-1.9 T/C, were genotyped in 95 AS patients and 130 healthy controls. The T allele of the PD-1.9 polymorphism was more frequent in the Korean male population with AS than in the Korean male controls (21.0% versus 6.9%, odds ratio 1.89, 95% confidence interval 1.483 to 2.408). The frequency of the CT haplotype (PD-1.5 C/T and PD-1.9 T/C) was higher in the AS patients (19%) than the controls (5.4%) (odds ratio 1.83, 95% confidence interval 1.559 to 2.521). The PD-1 polymorphism was demonstrated in Korean AS patients. The results suggest a genetic association between the PD-1 polymorphism and susceptibility to AS.