The interleukin-6 receptor Asp358Ala single nucleotide polymorphism rs2228145 confers increased proteolytic conversion rates by ADAM proteases

The pleiotropic activities of Interleukin (IL-)6 are controlled by membrane-bound and soluble forms of the IL-6 receptor (IL-6R) in processes called classic and trans-signaling, respectively. The coding single nucleotide polymorphism (SNP) rs2228145 of the Interleukin 6 receptor (IL-6R Asp358Ala variant) is associated with a 2-fold increase in soluble IL-6R (sIL-6R) serum levels resulting in reduced IL-6-induced C-reactive protein (CRP) production and a reduced risk for coronary heart disease. It was suggested that the increased sIL-6R level leads to decreased IL-6 classic or increased IL-6 trans-signaling. Irrespective of the functional outcome of increased sIL-6R serum level, it is still under debate, whether the increased sIL-6R serum levels emerged from differential splicing or ectodomain shedding. Here we show that increased proteolytic ectodomain shedding mediated by the A Disintegrin and metalloproteinase domain (ADAM) proteases ADAM10 and ADAM17 caused increased sIL-6R serum level in vitro as well as in healthy volunteers homozygous for the IL-6R Asp358Ala allele. Differential splicing of the IL-6R appears to have only a minor effect on sIL-6R level. Increased ectodomain shedding resulted in reduced cell-surface expression of the IL-6R Asp358Ala variant compared to the common IL-6R variant. In conclusion, increased IL-6R ectodomain shedding is a mechanistic explanation for the increased serum IL-6R levels found in persons homozygous for the rs2228145 IL-6R Asp358Ala variant.     

The role of osteopontin (OPN/SPP1) haplotypes in the susceptibility to Crohn's disease

Background

Osteopontin represents a multifunctional molecule playing a pivotal role in chronic inflammatory and autoimmune diseases. Its expression is increased in inflammatory bowel disease (IBD). The aim of our study was to analyze the association of osteopontin (OPN/SPP1) gene variants in a large cohort of IBD patients.

Methodology/Principal Findings

Genomic DNA from 2819 Caucasian individuals (n = 841 patients with Crohn''s disease (CD), n = 473 patients with ulcerative colitis (UC), and n = 1505 healthy unrelated controls) was analyzed for nine OPN SNPs (rs2728127, rs2853744, rs11730582, rs11739060, rs28357094, rs4754 = p.Asp80Asp, rs1126616 = p.Ala236Ala, rs1126772 and rs9138). Considering the important role of osteopontin in Th17-mediated diseases, we performed analysis for epistasis with IBD-associated IL23R variants and analyzed serum levels of the Th17 cytokine IL-22. For four OPN SNPs (rs4754, rs1126616, rs1126772 and rs9138), we observed significantly different distributions between male and female CD patients. rs4754 was protective in male CD patients (p = 0.0004, OR = 0.69). None of the other investigated OPN SNPs was associated with CD or UC susceptibility. However, several OPN haplotypes showed significant associations with CD susceptibility. The strongest association was found for a haplotype consisting of the 8 OPN SNPs rs2728127-rs2853744-rs11730582-rs11439060-rs28357094-rs112661-rs1126772-rs9138 (omnibus p-value = 2.07×10⁻⁸). Overall, the mean IL-22 secretion in the combined group of OPN minor allele carriers with CD was significantly lower than that of CD patients with OPN wildtype alleles (p = 3.66×10⁻⁵). There was evidence for weak epistasis between the OPN SNP rs28357094 with the IL23R SNP rs10489629 (p = 4.18×10⁻²) and between OPN SNP rs1126616 and IL23R SNP rs2201841 (p = 4.18×10⁻²) but none of these associations remained significant after Bonferroni correction.

Conclusions/Significance

Our study identified OPN haplotypes as modifiers of CD susceptibility, while the combined effects of certain OPN variants may modulate IL-22 secretion.     

The role of inflammatory pathway genetic variation on maternal metabolic phenotypes during pregnancy

Background

Since mediators of inflammation are associated with insulin resistance, and the risk of developing diabetes mellitus and gestational diabetes, we hypothesized that genetic variation in members of the inflammatory gene pathway impact glucose levels and related phenotypes in pregnancy. We evaluated this hypothesis by testing for association between genetic variants in 31 inflammatory pathway genes in the Hyperglycemia and Adverse Pregnancy Outcome (HAPO) cohort, a large multiethnic multicenter study designed to address the impact of glycemia less than overt diabetes on pregnancy outcome.

Results

Fasting, 1-hour, and 2-hour glucose, fasting and 1-hour C-peptide, and HbA1c levels were measured in blood samples obtained from HAPO participants during an oral glucose tolerance test at 24-32 weeks gestation. We tested for association between 458 SNPs mapping to 31 genes in the inflammatory pathway and metabolic phenotypes in 3836 European ancestry and 1713 Thai pregnant women. The strongest evidence for association was observed with TNF alpha and HbA1c (rs1052248; 0.04% increase per allele C; p-value = 4.4×10⁻⁵), RETN and fasting plasma glucose (rs1423096; 0.7 mg/dl decrease per allele A; p-value = 1.1×10⁻⁴), IL8 and 1 hr plasma glucose (rs2886920; 2.6 mg/dl decrease per allele T; p-value = 1.3×10⁻⁴), ADIPOR2 and fasting C-peptide (rs2041139; 0.55 ug/L decrease per allele A; p-value = 1.4×10⁻⁴), LEPR and 1-hour C-peptide (rs1171278; 0.62 ug/L decrease per allele T; p-value = 2.4×10⁻⁴), and IL6 and 1-hour plasma glucose (rs6954897; −2.29 mg/dl decrease per allele G, p-value = 4.3×10⁻⁴).

Conclusions

Based on the genes surveyed in this study the inflammatory pathway is unlikely to have a strong impact on maternal metabolic phenotypes in pregnancy although variation in individual members of the pathway (e.g. RETN, IL8, ADIPOR2, LEPR, IL6, and TNF alpha,) may contribute to metabolic phenotypes in pregnant women.     

The IL-33-ST2L Pathway Is Associated with Coronary Artery Disease in a Chinese Han Population

The effects of interleukin-33 (IL-33) on the immune system have been clearly demonstrated; however, in cardiovascular diseases, especially in coronary artery disease (CAD), these effects have not yet been clarified. In this study, we investigate the genetic role of the IL-33-ST2L pathway in CAD. We performed three-stage case-control association analyses on a total of 4,521 individuals with CAD and 4,809 controls via tag SNPs in the genes encoding IL-33 and ST2L-IL-1RL1. One tag SNP in each gene was significantly associated with CAD (rs7025417^T in IL33, p_adj = 1.19 × 10⁻²⁸, OR = 1.39, 95% CI: 1.31–1.47; rs11685424^G in IL1RL1, p_adj = 6.93 × 10⁻³⁰, OR = 1.40, 95% CI: 1.32–1.48). Combining significant variants in two genes, the risk for CAD increased nearly 5-fold (p_adj = 8.90 × 10⁻²¹, OR = 4.98, 95% CI: 3.56–6.97). Traditional risk factors for CAD were adjusted for the association studies by SPSS with logistic regression analysis. With the two variants above, both located within the gene promoter regions, reporter gene analysis indicated that the rs7025417 C>T and rs11685424 A>G changes resulted in altered regulation of IL33 and IL1RL1 gene expression, respectively (p < 0.005). Further studies revealed that the rs7025417 genotype was significantly associated with plasma IL-33 levels in the detectable subjects (n = 227, R² = 0.276, p = 1.77 × 10⁻¹⁷): the level of IL-33 protein increased with the number of rs7025417 risk (T) alleles. Based on genetic evidence in humans, the IL-33-ST2L pathway appears to have a causal role in the development of CAD, highlighting this pathway as a valuable target for the prevention and treatment of CAD.     

染色体1q12 区段IL6R 基因多态性与儿童哮喘易感性的研究

目的:检测染色体1q12区段IL6R基因多态性与儿童发生支气管哮喘易感性的关系。方法:150名支气管哮喘患儿为支气管哮喘组。150名健康儿童为对照组。采用质谱单核苷酸多态性(SNP)技术,对两组儿童的IL6R基因进行分析。结果:两组IL6基因位点的分布符合Hardy-Weinburg平衡定律。两组IL6R基因rs4845374位点的基因型与等位基因相比较,无明显差异(X2值分别为3.442和3.701;P值分别为0.179和0.088)。两组IL6R基因rs2228145位点的基因型与等位基因相比较,差异均有统计学意义(X2值分别为6.635和9.200;P值分别为0.036和0.003)。IL6R基因rs2228145位点基因型的变异等位基因T、C与支气管哮喘存在密切联系,CC、TT型出现支气管哮喘的风险均比CT型高。结论:IL6R基因rs4845374位点与儿童支气管哮喘无相关性,而rs2228145位点多态性与儿童支气管哮喘有相关性。     

The Soluble Interleukin-6 Receptor Is a Mediator of Hematopoietic and Skeletal Actions of Parathyroid Hormone

Both PTH and IL-6 signaling play pivotal roles in hematopoiesis and skeletal biology, but their interdependence is unclear. The purpose of this study was to evaluate the effect of IL-6 and soluble IL-6 receptor (sIL-6R) on hematopoietic and skeletal actions of PTH. In the bone microenvironment, PTH stimulated sIL-6R protein levels in primary osteoblast cultures in vitro and bone marrow in vivo in both IL-6^+/+ and IL-6^−/− mice. PTH-mediated hematopoietic cell expansion was attenuated in IL-6^−/− compared with IL-6^+/+ bone marrow, whereas sIL-6R treatment amplified PTH actions in IL-6^−/− earlier than IL-6^+/+ marrow cultures. Blocking sIL-6R signaling with sgp130 (soluble glycoprotein 130 receptor) inhibited PTH-dependent hematopoietic cell expansion in IL-6^−/− marrow. In the skeletal system, although intermittent PTH administration to IL-6^+/+ and IL-6^−/− mice resulted in similar anabolic actions, blocking sIL-6R significantly attenuated PTH anabolic actions. sIL-6R showed no direct effects on osteoblast proliferation or differentiation in vitro; however, it up-regulated myeloid cell expansion and production of the mesenchymal stem cell recruiting agent, TGF-β1 in the bone marrow microenvironment. Collectively, sIL-6R demonstrated orphan function and mediated PTH anabolic actions in bone in association with support of myeloid lineage cells in the hematopoietic system.     

Common variants in CRP and LEPR influence high sensitivity C-reactive protein levels in North Indians

Background

High sensitivity C-reactive protein (hsCRP) levels are shown to be influenced by genetic variants in Europeans; however, little is explored in Indian population.

Methods

Herein, we comprehensively evaluated association of all previously reported genetic determinants of hsCRP levels, including 18 cis (proximal to CRP gene) and 73 trans-acting (distal to CRP gene) variants in 4,200 North Indians of Indo-European ethnicity. First, we evaluated association of 91 variants from 12 candidate loci with hsCRP levels in 2,115 North Indians (1,042 non-diabetic subjects and 1,073 patients with type 2 diabetes). Then, cis and trans-acting variants contributing maximally to hsCRP level variation were further replicated in an independent 2,085 North Indians (1,047 patients with type 2 diabetes and 1,038 non-diabetic subjects).

Results

We found association of 12 variants from CRP, LEPR, IL1A, IL6, and IL6R with hsCRP levels in non-diabetic subjects. However, only rs3093059-CRP [β = 0.33, P = 9.6×10⁻⁵] and the haplotype harboring rs3093059 risk allele [β = 0.32 µg/mL, P = 1.4×10⁻⁴/P_perm = 9.0×10⁻⁴] retained significance after correcting for multiple testing. The cis-acting variant rs3093059-CRP had maximum contribution to the variance in hsCRP levels (1.14%). Among, trans-acting variants, rs1892534-LEPR was observed to contribute maximally to hsCRP level variance (0.59%). Associations of rs3093059-CRP and rs1892534-LEPR were confirmed by replication and attained higher significance after meta-analysis [β_meta = 0.26/0.22; P_meta = 4.3×10⁻⁷/7.4×10⁻³ and β_meta = −0.15/−0.12; P_meta = 2.0×10⁻⁶/1.6×10⁻⁶ for rs3093059 and rs1892534, respectively in non-diabetic subjects and all subjects taken together].

Conclusion

In conclusion, we identified rs3093059 in CRP and rs1892534 in LEPR as major cis and trans-acting contributor respectively, to the variance in hsCRP levels in North Indian population.     

A genome-wide association scan on the levels of markers of inflammation in Sardinians reveals associations that underpin its complex regulation

Identifying the genes that influence levels of pro-inflammatory molecules can help to elucidate the mechanisms underlying this process. We first conducted a two-stage genome-wide association scan (GWAS) for the key inflammatory biomarkers Interleukin-6 (IL-6), the general measure of inflammation erythrocyte sedimentation rate (ESR), monocyte chemotactic protein-1 (MCP-1), and high-sensitivity C-reactive protein (hsCRP) in a large cohort of individuals from the founder population of Sardinia. By analysing 731,213 autosomal or X chromosome SNPs and an additional ∼1.9 million imputed variants in 4,694 individuals, we identified several SNPs associated with the selected quantitative trait loci (QTLs) and replicated all the top signals in an independent sample of 1,392 individuals from the same population. Next, to increase power to detect and resolve associations, we further genotyped the whole cohort (6,145 individuals) for 293,875 variants included on the ImmunoChip and MetaboChip custom arrays. Overall, our combined approach led to the identification of 9 genome-wide significant novel independent signals—5 of which were identified only with the custom arrays—and provided confirmatory evidence for an additional 7. Novel signals include: for IL-6, in the ABO gene (rs657152, p = 2.13×10⁻²⁹); for ESR, at the HBB (rs4910472, p = 2.31×10⁻¹¹) and UCN119B/SPPL3 (rs11829037, p = 8.91×10⁻¹⁰) loci; for MCP-1, near its receptor CCR2 (rs17141006, p = 7.53×10⁻¹³) and in CADM3 (rs3026968, p = 7.63×10⁻¹³); for hsCRP, within the CRP gene (rs3093077, p = 5.73×10⁻²¹), near DARC (rs3845624, p = 1.43×10⁻¹⁰), UNC119B/SPPL3 (rs11829037, p = 1.50×10⁻¹⁴), and ICOSLG/AIRE (rs113459440, p = 1.54×10⁻⁰⁸) loci. Confirmatory evidence was found for IL-6 in the IL-6R gene (rs4129267); for ESR at CR1 (rs12567990) and TMEM57 (rs10903129); for MCP-1 at DARC (rs12075); and for hsCRP at CRP (rs1205), HNF1A (rs225918), and APOC-I (rs4420638). Our results improve the current knowledge of genetic variants underlying inflammation and provide novel clues for the understanding of the molecular mechanisms regulating this complex process.     

INSIG1 influences obesity-related hypertriglyceridemia in humans

In our analysis of a quantitative trait locus (QTL) for plasma triglyceride (TG) levels [logarithm of odds (LOD) = 3.7] on human chromosome 7q36, we examined 29 single nucleotide polymorphisms (SNPs) across INSIG1, a biological candidate gene in the region. Insulin-induced genes (INSIGs) are feedback mediators of cholesterol and fatty acid synthesis in animals, but their role in human lipid regulation is unclear. In our cohort, the INSIG1 promoter SNP rs2721 was associated with TG levels (P = 2 × 10⁻³ in 1,560 individuals of the original linkage cohort, P = 8 × 10⁻⁴ in 920 unrelated individuals of the replication cohort, combined P = 9.9 × 10⁻⁶). Individuals homozygous for the T allele had 9% higher TG levels and 2-fold lower expression of INSIG1 in surgical liver biopsy samples when compared with individuals homozygous for the G allele. Also, the T allele showed additional binding of nuclear proteins from HepG2 liver cells in gel shift assays. Finally, the variant rs7566605 in INSIG2, the only homolog of INSIG1, enhances the effect of rs2721 (P = 0.00117). The variant rs2721 alone explains 5.4% of the observed linkage in our cohort, suggesting that additional, yet-undiscovered genes and sequence variants in the QTL interval also contribute to alterations in TG levels in humans.     

PAX6 haplotypes are associated with high myopia in Han chinese

Background

The paired box 6 (PAX6) gene is considered as a master gene for eye development. Linkage of myopia to the PAX6 region on chromosome 11p13 was shown in several studies, but the results for association between myopia and PAX6 were inconsistent so far.

Methodology/Principal Findings

We genotyped 16 single nucleotide polymorphisms (SNPs) in the PAX6 gene and its regulatory regions in an initial study for 300 high myopia cases and 300 controls (Group 1), and successfully replicated the positive results with another independent group of 299 high myopia cases and 299 controls (Group 2). Five SNPs were genotyped in the replication study. The spherical equivalent of subjects with high myopia was ≤−8.0 dioptres. The PLINK package was used for genetic data analysis. No association was found between each of the SNPs and high myopia. However, exhaustive sliding-window haplotype analysis highlighted an important role for rs12421026 because haplotypes containing this SNP were found to be associated with high myopia. The most significant results were given by the 4-SNP haplotype window consisting of rs2071754, rs3026393, rs1506 and rs12421026 (P = 3.54×10⁻¹⁰, 4.06×10⁻¹¹ and 1.56×10⁻¹⁸ for Group 1, Group 2 and Combined Group, respectively) and the 3-SNP haplotype window composed of rs3026393, rs1506 and rs12421026 (P = 5.48×10⁻¹⁰, 7.93×10⁻¹² and 6.28×10⁻²³ for the three respective groups). The results remained significant after correction for multiple comparisons by permutations. The associated haplotyes found in a previous study were also successfully replicated in this study.

Conclusions/Significance

PAX6 haplotypes are associated with susceptibility to the development of high myopia in Chinese. The PAX6 locus plays a role in high myopia.     

Genome-wide association meta-analysis of cortical bone mineral density unravels allelic heterogeneity at the RANKL locus and potential pleiotropic effects on bone

Previous genome-wide association (GWA) studies have identified SNPs associated with areal bone mineral density (aBMD). However, this measure is influenced by several different skeletal parameters, such as periosteal expansion, cortical bone mineral density (BMD_C) cortical thickness, trabecular number, and trabecular thickness, which may be under distinct biological and genetic control. We have carried out a GWA and replication study of BMD_C, as measured by peripheral quantitative computed tomography (pQCT), a more homogenous and valid measure of actual volumetric bone density. After initial GWA meta-analysis of two cohorts (ALSPAC n = 999, aged ∼15 years and GOOD n = 935, aged ∼19 years), we attempted to replicate the BMD_C associations that had p<1×10⁻⁵ in an independent sample of ALSPAC children (n = 2803) and in a cohort of elderly men (MrOS Sweden, n = 1052). The rs1021188 SNP (near RANKL) was associated with BMD_C in all cohorts (overall p = 2×10⁻¹⁴, n = 5739). Each minor allele was associated with a decrease in BMD_C of ∼0.14SD. There was also evidence for an interaction between this variant and sex (p = 0.01), with a stronger effect in males than females (at age 15, males −6.77mg/cm³ per C allele, p = 2×10⁻⁶; females −2.79 mg/cm³ per C allele, p = 0.004). Furthermore, in a preliminary analysis, the rs1021188 minor C allele was associated with higher circulating levels of sRANKL (p<0.005). We show this variant to be independent from the previously aBMD associated SNP (rs9594738) and possibly from a third variant in the same RANKL region, which demonstrates important allelic heterogeneity at this locus. Associations with skeletal parameters reflecting bone dimensions were either not found or were much less pronounced. This finding implicates RANKL as a locus containing variation associated with volumetric bone density and provides further insight into the mechanism by which the RANK/RANKL/OPG pathway may be involved in skeletal development.     

Caucasian and Asian specific rheumatoid arthritis risk loci reveal limited replication and apparent allelic heterogeneity in north Indians

Genome-wide association studies and meta-analysis indicate that several genes/loci are consistently associated with rheumatoid arthritis (RA) in European and Asian populations. To evaluate the transferability status of these findings to an ethnically diverse north Indian population, we performed a replication analysis. We investigated the association of 47 single-nucleotide polymorphisms (SNPs) at 43 of these genes/loci with RA in a north Indian cohort comprising 983 RA cases and 1007 age and gender matched controls. Genotyping was done using Infinium human 660w-quad. Association analysis by chi-square test implemented in plink was carried out in two steps. Firstly, association of the index or surrogate SNP (r2>0.8, calculated from reference GIH Hap-Map population) was tested. In the second step, evidence for allelic/locus heterogeneity at aforementioned genes/loci was assessed for by testing additional flanking SNPs in linkage equilibrium with index/surrogate marker.Of the 44 European specific index SNPs, neither index nor surrogate SNPs were present for nine SNPs in the genotyping array. Of the remaining 35, associations were replicated at seven genes namely PTPN22 (rs1217407, p = 3×10⁻³); IL2–21 (rs13119723, p = 0.008); HLA-DRB1 (rs660895, p = 2.56×10⁻⁵; rs6457617, p = 1.6×10⁻⁰⁹; rs13192471, p = 6.7×10⁻¹⁶); TNFA1P3 (rs9321637, p = 0.03); CCL21 (rs13293020, p = 0.01); IL2RA (rs2104286, p = 1.9×10⁻⁴) and ZEB1 (rs2793108, p = 0.006). Of the three Asian specific loci tested, rs2977227 in PADI4 showed modest association (p<0.02). Further, of the 140 SNPs (in LE with index/surrogate variant) tested, association was observed at 11 additional genes: PTPRC, AFF3, CD28, CTLA4, PXK, ANKRD55, TAGAP, CCR6, BLK, CD40 and IL2RB. This study indicates limited replication of European and Asian index SNPs and apparent allelic heterogeneity in RA etiology among north Indians warranting independent GWAS in this population. However, replicated associations of HLA-DRB1, PTPN22 (which confer ∼50% of the heritable risk to RA) and IL2RA suggest that cross-ethnicity fine mapping of such loci is apposite for identification of causal variants.     

Non-synonymous polymorphisms in the P2RX 4 are related to bone mineral density and osteoporosis risk in a cohort of Dutch fracture patients

In the present study we investigated whether single nucleotide polymorphisms (SNPs) in the P2RX₄, which alter the P2X₄R function, are associated with the development of osteoporosis and whether an interaction between the P2X₄R and P2X₇R confer a synergistic effect of these two receptors on osteoporosis risk. Patients with fracture (690 females and 231 males, aged ≥50 years) were genotyped for three non-synonymous P2X₄R SNPs. Bone mineral density (BMD) was measured at the total hip, lumbar spine, and femoral neck. Subject carrying the variant allele of the Tyr315Cys polymorphism showed a 2.68-fold (95 % CI, 1.20–6.02) higher risk of osteoporosis compared with wild-type subject. Furthermore, significant lower lumbar spine BMD values were observed in subjects carrying the Cys315 allele as compared with wild-type (0.85 ± 0.17 and 0.93 ± 0.17 g/cm², respectively; p < 0.001). Assuming a recessive model, carriers of the variant allele of the Ser242Gly polymorphism showed increased BMD values at the lumbar spine compare to wild-type subject (1.11 ± 0.35 and 0.92 ± 0.17 g/cm², respectively; p = 0.0045). This is the first study demonstrating an association of non-synonymous polymorphisms in the P2RX₄ and the risk of osteoporosis, suggesting a role of the P2X₄R in the regulation of bone mass.

Electronic supplementary material

The online version of this article (doi:10.1007/s11302-012-9337-0) contains supplementary material, which is available to authorized users.     

Investigation of relationship between IL-6 gene variants and hypertension in Turkish population

Hypertension (HT) is a common and life threating health problem worldwide leading to stroke, heart attack and renal failure. It is characterized by elevated blood pressure forced heart load. Human interleukin-6 (IL-6) and C- reactive protein (CRP) are known to be involved in inflammatory processes. IL-6 gene is a polymorphic gene which −174 G/C is a common and −572 G/C is a rare polymorphisms identified in promoter region. Publications on IL-6 gene polymorphisms raised the question whether this gene polymorphisms lead to susceptibility to HT or not. To investigate the effects of IL-6 gene −174 G/C (rs 1800795) and −572 G/C (rs1800796) polymorphisms on plasma IL-6 and CRP levels and their associations with hypertension disease in Turkish population we analyzed −174 G/C and −572 G/C polymorphisms and plasma IL-6 and CRP levels in 111 healthy controls and 108 hypertension patients from Adıyaman, Turkey. We determined the genotypes using polymerase chain reaction-restriction fragment length polymorphism and analyzed plasma levels of IL-6 by ELISA and CRP by automated standard biochemical methods. We have found no statistically significant differences between IL-6 gene −174 G/C and −572 G/C genotypes and allelic frequencies and IL-6 and CRP plasma levels and HT (p > 0.05). No CC genotype was found in control subjects for −572 G/C polymorphism. In conclusion, we found relation to −174 G/C and −572 G/C gene variants between neither IL-6 and CRP levels nor hypertension. The −572 G allele and GG genotype are predominant in Turkish population in Adıyaman, Turkey whereas the CC genotype is very rare.     

Association of variants at 1q32 and STAT3 with ankylosing spondylitis suggests genetic overlap with Crohn's disease

Ankylosing spondylitis (AS) is a common inflammatory arthritic condition. Overt inflammatory bowel disease (IBD) occurs in about 10% of AS patients, and in addition 70% of AS cases may have subclinical terminal ileitis. Spondyloarthritis is also common in IBD patients. We therefore tested Crohn''s disease susceptibility genes for association with AS, aiming to identify pleiotropic genetic associations with both diseases. Genotyping was carried out using Sequenom and Applied Biosystems TaqMan and OpenArray technologies on 53 markers selected from 30 Crohn''s disease associated genomic regions. We tested genotypes in a population of unrelated individual cases (n = 2,773) and controls (n = 2,215) of white European ancestry for association with AS. Statistical analysis was carried out using a Cochran-Armitage test for trend in PLINK. Strong association was detected at chr1q32 near KIF21B (rs11584383, P = 1.6×10⁻¹⁰, odds ratio (OR) = 0.74, 95% CI:0.68–0.82). Association with disease was also detected for 2 variants within STAT3 (rs6503695, P = 4.6×10⁻⁴. OR = 0.86 (95% CI:0.79–0.93); rs744166, P = 2.6×10⁻⁵, OR = 0.84 (95% CI:0.77–0.91)). Association was confirmed for IL23R (rs11465804, P = 1.2×10⁻⁵, OR = 0.65 (95% CI:0.54–0.79)), and further associations were detected for IL12B (rs10045431, P = 5.2×10⁻⁵, OR = 0.83 (95% CI:0.76–0.91)), CDKAL1 (rs6908425, P = 1.1×10⁻⁴, OR = 0.82 (95% CI:0.74–0.91)), LRRK2/MUC19 (rs11175593, P = 9.9×10⁻⁵, OR = 1.92 (95% CI: 1.38–2.67)), and chr13q14 (rs3764147, P = 5.9×10⁻⁴, OR = 1.19 (95% CI: 1.08–1.31)). Excluding cases with clinical IBD did not significantly affect these findings. This study identifies chr1q32 and STAT3 as ankylosing spondylitis susceptibility loci. It also further confirms association for IL23R and detects suggestive association with another 4 loci. STAT3 is a key signaling molecule within the Th17 lymphocyte differentiation pathway and further enhances the case for a major role of this T-lymphocyte subset in ankylosing spondylitis. Finally these findings suggest common aetiopathogenic pathways for AS and Crohn''s disease and further highlight the involvement of common risk variants across multiple diseases.     

Genome-wide Association Study Reveals Multiple Nasopharyngeal Carcinoma-Associated Loci within the HLA Region at Chromosome 6p21.3

Nasopharyngeal carcinoma (NPC) is a multifactorial malignancy closely associated with genetic factors and Epstein-Barr virus infection. To identify the common genetic variants linked to NPC susceptibility, we conducted a genome-wide association study (GWAS) in 277 NPC patients and 285 healthy controls within the Taiwanese population, analyzing 480,365 single-nucleotide polymorphisms (SNPs). Twelve statistically significant SNPs were identified and mapped to chromosome 6p21.3. Associations were replicated in two independent sets of case-control samples. Two of the most significant SNPs (rs2517713 and rs2975042; p_combined = 3.9 × 10⁻²⁰ and 1.6 × 10⁻¹⁹, respectively) were located in the HLA-A gene. Moreover, we detected significant associations between NPC and two genes: specifically, gamma aminobutyric acid b receptor 1 (GABBR1) (rs29232; p_combined = 8.97 × 10⁻¹⁷) and HLA-F (rs3129055 and rs9258122; p_combined = 7.36 × 10⁻¹¹ and 3.33 × 10⁻¹⁰, respectively). Notably, the association of rs29232 remained significant (residual p < 5 × 10⁻⁴) after adjustment for age, gender, and HLA-related SNPs. Furthermore, higher GABA_B receptor 1 expression levels can be found in the tumor cells in comparison to the adjacent epithelial cells (p < 0.001) in NPC biopsies, implying a biological role of GABBR1 in NPC carcinogenesis. To our knowledge, it is the first GWAS report of NPC showing that multiple loci (HLA-A, HLA-F, and GABBR1) within chromosome 6p21.3 are associated with NPC. Although some of these relationships may be attributed to linkage disequilibrium between the loci, the findings clearly provide a fresh direction for the study of NPC development.     

CFH, VEGF, and PEDF genotypes and the response to intravitreous injection of bevacizumab for the treatment of age-related macular degeneration

We determined whether there is an association between complement factor H (CFH), high-temperature requirement A-1 (HTRA1), vascular endothelial growth factor (VEGF), and pigment epithelium-derived factor (PEDF) genotypes and the response to treatment with a single intravitreous injection of bevacizumab for age-related macular degeneration (AMD). Eighty-three patients with exudative AMD treated by bevacizumab injection were genotyped for three single nucleotide polymorphisms (SNPs; rs800292, rs1061170, rs1410996) in the CFH gene, a rs11200638-SNP in the HTRA1 gene, three SNPs (rs699947, rs1570360, rs2010963) in the VEGF gene, and four SNPs (rs12150053, rs12948385, rs9913583, rs1136287) in the PEDF gene using a TaqMan assay. The CT genotype (heterozygous) of CFH-rs1061170 was more frequently represented in nonresponders in vision than TT genotypes (nonrisk allele homozygous) at the time points of 1 and 3 months, while there was no CC genotype (risk allele homozygous) in our study cohort (p = 7.66 × 10⁻³, 7.83 × 10⁻³, respectively). VEGF-rs699947 was also associated with vision changes at 1 month and PEDF-rs1136287 at 3 months (p = 5.11 × 10⁻³, 2.05 × 10⁻², respectively). These variants may be utilized for genetic biomarkers to estimate visual outcomes in the response to intravitreal bevacizumab treatment for AMD.     

A Weighted False Discovery Rate Control Procedure Reveals Alleles at FOXA2 that Influence Fasting Glucose Levels

Association signals in GWAS are usually prioritized solely by p values. Here, we attempt to improve the power of GWAS by using a weighted false discovery rate control procedure to detect associations of low-frequency variants with effect sizes similar to or even larger than those of common variants. We used the Affymetrix Genome-Wide Human SNP Array 6.0 to test for association with fasting glucose levels in the Atherosclerosis Risk in Communities Study (ARIC) population. In addition to finding several previously identified sequence variations, we identified a low-frequency variant (rs1209523; minor allele frequency = 0.043) near FOXA2 that was associated with fasting glucose levels in European Americans (EAs) (n = 7428, p value = 1.3 × 10⁻⁵). The association between rs1209523 and glucose levels was also significant in African Americans (AAs) (n = 2029, p value = 6.7 × 10⁻³) of the ARIC and was confirmed by replication in both EAs and AAs of the Dallas Heart Study (n = 963 and 1571, respectively; p values = 5.3 × 10⁻³ and 5.8 × 10⁻⁴, respectively) and in EAs of the Cooper Center Longitudinal Study (n = 2862; p value = 1.6 × 10⁻²). A meta-analysis of these five populations yielded an estimated effect size of −1.31 mg/dl per minor allele (p value = 2.2 × 10⁻¹¹). This study reveals that there is a cache of less-frequent variants in GWAS arrays that can be identified via analytical approaches accounting for allele frequencies.