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.     

PRKCA: A Positional Candidate Gene for Body Mass Index and Asthma

Asthma incidence and prevalence are higher in obese individuals. A potential mechanistic basis for this relationship is pleiotropy. We hypothesized that significant linkage and candidate-gene association would be found for body mass index (BMI) in a population ascertained on asthma affection status. Linkage analysis for BMI was performed on 657 subjects in eight Costa Rican families enrolled in a study of asthma. Family-based association studies were conducted for BMI with SNPs within a positional candidate gene, PRKCA. SNPs within PRKCA were also tested for association with asthma. Association studies were conducted in 415 Costa Rican parent-child trios and 493 trios participating in the Childhood Asthma Management Program (CAMP). Although only modest evidence of linkage for BMI was obtained for the whole cohort, significant linkage was noted for BMI in females on chromosome 17q (peak LOD = 3.39). Four SNPs in a candidate gene in this region (PRKCA) had unadjusted association p values < 0.05 for BMI in both cohorts, with the joint p value for two SNPs remaining significant after adjustment for multiple comparisons (rs228883 and rs1005651, joint p values = 9.5 × 10⁻⁵ and 5.6 × 10⁻⁵). Similarly, eight SNPs had unadjusted association p values < 0.05 for asthma in both populations, with one SNP remaining significant after adjustment for multiple comparisons (rs11079657, joint p value = 2.6 × 10⁻⁵). PRKCA is a pleiotropic locus that is associated with both BMI and asthma and that has been identified via linkage analysis of BMI in a population ascertained on asthma.     

Variants near FOXE1 are associated with hypothyroidism and other thyroid conditions: using electronic medical records for genome- and phenome-wide studies

We repurposed existing genotypes in DNA biobanks across the Electronic Medical Records and Genomics network to perform a genome-wide association study for primary hypothyroidism, the most common thyroid disease. Electronic selection algorithms incorporating billing codes, laboratory values, text queries, and medication records identified 1317 cases and 5053 controls of European ancestry within five electronic medical records (EMRs); the algorithms'' positive predictive values were 92.4% and 98.5% for cases and controls, respectively. Four single-nucleotide polymorphisms (SNPs) in linkage disequilibrium at 9q22 near FOXE1 were associated with hypothyroidism at genome-wide significance, the strongest being rs7850258 (odds ratio [OR] 0.74, p = 3.96 × 10⁻⁹). This association was replicated in a set of 263 cases and 1616 controls (OR = 0.60, p = 5.7 × 10⁻⁶). A phenome-wide association study (PheWAS) that was performed on this locus with 13,617 individuals and more than 200,000 patient-years of billing data identified associations with additional phenotypes: thyroiditis (OR = 0.58, p = 1.4 × 10⁻⁵), nodular (OR = 0.76, p = 3.1 × 10⁻⁵) and multinodular (OR = 0.69, p = 3.9 × 10⁻⁵) goiters, and thyrotoxicosis (OR = 0.76, p = 1.5 × 10⁻³), but not Graves disease (OR = 1.03, p = 0.82). Thyroid cancer, previously associated with this locus, was not significantly associated in the PheWAS (OR = 1.29, p = 0.09). The strongest association in the PheWAS was hypothyroidism (OR = 0.76, p = 2.7 × 10⁻¹³), which had an odds ratio that was nearly identical to that of the curated case-control population in the primary analysis, providing further validation of the PheWAS method. Our findings indicate that EMR-linked genomic data could allow discovery of genes associated with many diseases without additional genotyping cost.     

Genetic variants at 13q12.12 are associated with high myopia in the Han Chinese population

High myopia, which is extremely prevalent in the Chinese population, is one of the leading causes of blindness in the world. Genetic factors play a critical role in the development of the condition. To identify the genetic variants associated with high myopia in the Han Chinese, we conducted a genome-wide association study (GWAS) of 493,947 SNPs in 1088 individuals (419 cases and 669 controls) from a Han Chinese cohort and followed up on signals that were associated with p < 1.0 × 10⁻⁴ in three independent cohorts (combined, 2803 cases and 5642 controls). We identified a significant association between high myopia and a variant at 13q12.12 (rs9318086, combined p = 1.91 × 10⁻¹⁶, heterozygous odds ratio = 1.32, and homozygous odds ratio = 1.64). Furthermore, five additional SNPs (rs9510902, rs3794338, rs1886970, rs7325450, and rs7331047) in the same linkage disequilibrium (LD) block with rs9318086 also proved to be significantly associated with high myopia in the Han Chinese population; p values ranged from 5.46 × 10⁻¹¹ to 6.16 × 10⁻¹⁶. This associated locus contains three genes—MIPEP, C1QTNF9B-AS1, and C1QTNF9B. MIPEP and C1QTNF9B were found to be expressed in the retina and retinal pigment epithelium (RPE) and are more likely than C1QTNF9B-AS1 to be associated with high myopia given the evidence of retinal signaling that controls eye growth. Our results suggest that the variants at 13q12.12 are associated with high myopia.     

Genome-wide screen identifies rs646776 near sortilin as a regulator of progranulin levels in human plasma

Recent studies suggest progranulin (GRN) is a neurotrophic factor. Loss-of-function mutations in the progranulin gene (GRN) cause frontotemporal lobar degeneration (FTLD), a progressive neurodegenerative disease affecting ∼10% of early-onset dementia patients. Using an enzyme-linked immunosorbent assay, we previously showed that GRN is detectable in human plasma and can be used to predict GRN mutation status. This study also showed a wide range in plasma GRN levels in non-GRN mutation carriers, including controls. We have now performed a genome-wide association study of 313,504 single-nucleotide polymorphisms (SNPs) in 533 control samples and identified on chromosome 1p13.3 two SNPs with genome-wide significant association with plasma GRN levels (top SNP rs646776; p = 1.7 × 10⁻³⁰). The association of rs646776 with plasma GRN levels was replicated in two independent series of 508 controls (p = 1.9 × 10⁻¹⁹) and 197 FTLD patients (p = 6.4 × 10⁻¹²). Overall, each copy of the minor C allele decreased GRN levels by ∼15%. SNP rs646776 is located near sortilin (SORT1), and the minor C allele of rs646776 was previously associated with increased SORT1 mRNA levels. Supporting these findings, overexpression of SORT1 in cultured HeLa cells dramatically reduced GRN levels in the conditioned media, whereas knockdown of SORT1 increased extracellular GRN levels. In summary, we identified significant association of a locus on chromosome 1p13.3 with plasma GRN levels through an unbiased genome-wide screening approach and implicated SORT1 as an important regulator of GRN levels. This finding opens avenues for future research into GRN biology and the pathophysiology of neurodegenerative diseases.     

Variants in TF and HFE Explain ∼40% of Genetic Variation in Serum-Transferrin Levels

Only a small proportion of genetic variation in complex traits has been explained by SNPs from genome-wide association studies (GWASs). We report the results from two GWASs for serum markers of iron status (serum iron, serum transferrin, transferrin saturation with iron, and serum ferritin), which are important in iron overload (e.g., hemochromatosis) and deficiency (e.g., anemia) conditions. We performed two GWASs on samples of Australians of European descent. In the first GWAS, 411 adolescent twins and their siblings were genotyped with 100K SNPs. rs1830084, 10.8 kb 3′ of TF, was significantly associated with serum transferrin (p total association test = 1.0 × 10⁻⁹; p within-family test = 2.2 × 10⁻⁵). In the second GWAS on an independent sample of 459 female monozygotic (MZ) twin pairs genotyped with 300K SNPs, we found rs3811647 (within intron 11 of TF, HapMap CEU r² with rs1830084 = 0.86) was significantly associated with serum transferrin (p = 3.0 × 10⁻¹⁵). In the second GWAS, we found two additional and independent SNPs on TF (rs1799852 and rs2280673) and confirmed the known C282Y mutation in HFE to be independently associated with serum transferrin. The three variants in TF (rs3811647, rs1799852 and rs2280673) plus the HFE C282Y mutation explained ~40% of genetic variation in serum transferrin (p = 7.8 × 10⁻²⁵). These findings are potentially important for our understanding of iron metabolism and of regulation of hepatic protein secretion, and also strongly support the hypothesis that the genetic architecture of some endophenotypes may be simpler than that of disease.     

A genome-wide association search for type 2 diabetes genes in African Americans

African Americans are disproportionately affected by type 2 diabetes (T2DM) yet few studies have examined T2DM using genome-wide association approaches in this ethnicity. The aim of this study was to identify genes associated with T2DM in the African American population. We performed a Genome Wide Association Study (GWAS) using the Affymetrix 6.0 array in 965 African-American cases with T2DM and end-stage renal disease (T2DM-ESRD) and 1029 population-based controls. The most significant SNPs (n = 550 independent loci) were genotyped in a replication cohort and 122 SNPs (n = 98 independent loci) were further tested through genotyping three additional validation cohorts followed by meta-analysis in all five cohorts totaling 3,132 cases and 3,317 controls. Twelve SNPs had evidence of association in the GWAS (P<0.0071), were directionally consistent in the Replication cohort and were associated with T2DM in subjects without nephropathy (P<0.05). Meta-analysis in all cases and controls revealed a single SNP reaching genome-wide significance (P<2.5×10⁻⁸). SNP rs7560163 (P = 7.0×10⁻⁹, OR (95% CI) = 0.75 (0.67–0.84)) is located intergenically between RND3 and RBM43. Four additional loci (rs7542900, rs4659485, rs2722769 and rs7107217) were associated with T2DM (P<0.05) and reached more nominal levels of significance (P<2.5×10⁻⁵) in the overall analysis and may represent novel loci that contribute to T2DM. We have identified novel T2DM-susceptibility variants in the African-American population. Notably, T2DM risk was associated with the major allele and implies an interesting genetic architecture in this population. These results suggest that multiple loci underlie T2DM susceptibility in the African-American population and that these loci are distinct from those identified in other ethnic populations.     

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.     

Association of JAG1 with Bone Mineral Density and Osteoporotic Fractures: A Genome-wide Association Study and Follow-up Replication Studies

Bone mineral density (BMD), a diagnostic parameter for osteoporosis and a clinical predictor of fracture, is a polygenic trait with high heritability. To identify genetic variants that influence BMD in different ethnic groups, we performed a genome-wide association study (GWAS) on 800 unrelated Southern Chinese women with extreme BMD and carried out follow-up replication studies in six independent study populations of European descent and Asian populations including 18,098 subjects. In the meta-analysis, rs2273061 of the Jagged1 (JAG1) gene was associated with high BMD (p = 5.27 × 10⁻⁸ for lumbar spine [LS] and p = 4.15 × 10⁻⁵ for femoral neck [FN], n = 18,898). This SNP was further found to be associated with the low risk of osteoporotic fracture (p = 0.009, OR = 0.7, 95% CI 0.57–0.93, n = 1881). Region-wide and haplotype analysis showed that the strongest association evidence was from the linkage disequilibrium block 5, which included rs2273061 of the JAG1 gene (p = 8.52 × 10⁻⁹ for LS and 3.47 × 10⁻⁵ at FN). To assess the function of identified variants, an electrophoretic mobility shift assay demonstrated the binding of c-Myc to the “G” but not “A” allele of rs2273061. A mRNA expression study in both human bone-derived cells and peripheral blood mononuclear cells confirmed association of the high BMD-related allele G of rs2273061 with higher JAG1 expression. Our results identify the JAG1 gene as a candidate for BMD regulation in different ethnic groups, and it is a potential key factor for fracture pathogenesis.     

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.     

Integration of Genomic and Genetic Approaches Implicates IREB2 as a COPD Susceptibility Gene

Chronic obstructive pulmonary disease (COPD) is the fourth leading cause of death worldwide and is influenced by both genetic determinants and smoking. We identified genomic regions from 56 lung-tissue gene-expression microarrays and used them to select 889 SNPs to be tested for association with COPD. We genotyped SNPs in 389 severe COPD cases from the National Emphysema Treatment Trial and 424 cigarette-smoking controls from the Normative Aging Study. A total of 71 autosomal SNPs demonstrated at least nominal significance with COPD susceptibility (p = 3.4 × 10⁻⁶ to 0.05). These 71 SNPs were evaluated in a family-based study of 127 probands with severe, early-onset COPD and 822 of their family members in the Boston Early-Onset COPD Study. We combined p values from the case-control and family-based analyses, setting p = 5.60 × 10⁻⁵ as a conservative threshold for significance. Three SNPs in the iron regulatory protein 2 (IREB2) gene met this stringent threshold for significance, and four other IREB2 SNPs demonstrated combined p < 0.02. We demonstrated replication of association for these seven IREB2 SNPs (all p values ≤ 0.02) in a family-based study of 3117 subjects from the International COPD Genetics Network; combined p values across all cohorts for the main phenotype of interest ranged from 1.6 × 10⁻⁷ to 6.4 × 10⁻⁴. IREB2 protein and mRNA were increased in lung-tissue samples from COPD subjects in comparison to controls. In summary, gene-expression and genetic-association results have implicated IREB2 as a COPD susceptibility gene.     

Genome-wide scan identifies TNIP1, PSORS1C1, and RHOB as novel risk loci for systemic sclerosis

Systemic sclerosis (SSc) is an orphan, complex, inflammatory disease affecting the immune system and connective tissue. SSc stands out as a severely incapacitating and life-threatening inflammatory rheumatic disease, with a largely unknown pathogenesis. We have designed a two-stage genome-wide association study of SSc using case-control samples from France, Italy, Germany, and Northern Europe. The initial genome-wide scan was conducted in a French post quality-control sample of 564 cases and 1,776 controls, using almost 500 K SNPs. Two SNPs from the MHC region, together with the 6 loci outside MHC having at least one SNP with a P<10⁻⁵ were selected for follow-up analysis. These markers were genotyped in a post-QC replication sample of 1,682 SSc cases and 3,926 controls. The three top SNPs are in strong linkage disequilibrium and located on 6p21, in the HLA-DQB1 gene: rs9275224, P = 9.18×10⁻⁸, OR = 0.69, 95% CI [0.60–0.79]; rs6457617, P = 1.14×10⁻⁷ and rs9275245, P = 1.39×10⁻⁷. Within the MHC region, the next most associated SNP (rs3130573, P = 1.86×10⁻⁵, OR = 1.36 [1.18–1.56]) is located in the PSORS1C1 gene. Outside the MHC region, our GWAS analysis revealed 7 top SNPs (P<10⁻⁵) that spanned 6 independent genomic regions. Follow-up of the 17 top SNPs in an independent sample of 1,682 SSc and 3,926 controls showed associations at PSORS1C1 (overall P = 5.70×10⁻¹⁰, OR:1.25), TNIP1 (P = 4.68×10⁻⁹, OR:1.31), and RHOB loci (P = 3.17×10⁻⁶, OR:1.21). Because of its biological relevance, and previous reports of genetic association at this locus with connective tissue disorders, we investigated TNIP1 expression. A markedly reduced expression of the TNIP1 gene and also its protein product were observed both in lesional skin tissue and in cultured dermal fibroblasts from SSc patients. Furthermore, TNIP1 showed in vitro inhibitory effects on inflammatory cytokine-induced collagen production. The genetic signal of association with TNIP1 variants, together with tissular and cellular investigations, suggests that this pathway has a critical role in regulating autoimmunity and SSc pathogenesis.     

Abdominal aortic aneurysm is associated with a variant in low-density lipoprotein receptor-related protein 1

Abdominal aortic aneurysm (AAA) is a common cause of morbidity and mortality and has a significant heritability. We carried out a genome-wide association discovery study of 1866 patients with AAA and 5435 controls and replication of promising signals (lead SNP with a p value < 1 × 10⁻⁵) in 2871 additional cases and 32,687 controls and performed further follow-up in 1491 AAA and 11,060 controls. In the discovery study, nine loci demonstrated association with AAA (p < 1 × 10⁻⁵). In the replication sample, the lead SNP at one of these loci, rs1466535, located within intron 1 of low-density-lipoprotein receptor-related protein 1 (LRP1) demonstrated significant association (p = 0.0042). We confirmed the association of rs1466535 and AAA in our follow-up study (p = 0.035). In a combined analysis (6228 AAA and 49182 controls), rs1466535 had a consistent effect size and direction in all sample sets (combined p = 4.52 × 10⁻¹⁰, odds ratio 1.15 [1.10–1.21]). No associations were seen for either rs1466535 or the 12q13.3 locus in independent association studies of coronary artery disease, blood pressure, diabetes, or hyperlipidaemia, suggesting that this locus is specific to AAA. Gene-expression studies demonstrated a trend toward increased LRP1 expression for the rs1466535 CC genotype in arterial tissues; there was a significant (p = 0.029) 1.19-fold (1.04–1.36) increase in LRP1 expression in CC homozygotes compared to TT homozygotes in aortic adventitia. Functional studies demonstrated that rs1466535 might alter a SREBP-1 binding site and influence enhancer activity at the locus. In conclusion, this study has identified a biologically plausible genetic variant associated specifically with AAA, and we suggest that this variant has a possible functional role in LRP1 expression.     

Sequence Variants in Three Loci Influence Monocyte Counts and Erythrocyte Volume

Blood cells participate in vital physiological processes, and their numbers are tightly regulated so that homeostasis is maintained. Disruption of key regulatory mechanisms underlies many blood-related Mendelian diseases but also contributes to more common disorders, including atherosclerosis. We searched for quantitative trait loci (QTL) for hematology traits through a whole-genome association study, because these could provide new insights into both hemopoeitic and disease mechanisms. We tested 1.8 million variants for association with 13 hematology traits measured in 6015 individuals from the Australian and Dutch populations. These traits included hemoglobin composition, platelet counts, and red blood cell and white blood cell indices. We identified three regions of strong association that, to our knowledge, have not been previously reported in the literature. The first was located in an intergenic region of chromosome 9q31 near LPAR1, explaining 1.5% of the variation in monocyte counts (best SNP rs7023923, p = 8.9 × 10⁻¹⁴). The second locus was located on chromosome 6p21 and associated with mean cell erythrocyte volume (rs12661667, p = 1.2 × 10⁻⁹, 0.7% variance explained) in a region that spanned five genes, including CCND3, a member of the D-cyclin gene family that is involved in hematopoietic stem cell expansion. The third region was also associated with erythrocyte volume and was located in an intergenic region on chromosome 6q24 (rs592423, p = 5.3 × 10⁻⁹, 0.6% variance explained). All three loci replicated in an independent panel of 1543 individuals (p values = 0.001, 9.9 × 10⁻⁵, and 7 × 10⁻⁵, respectively). The identification of these QTL provides new opportunities for furthering our understanding of the mechanisms regulating hemopoietic cell fate.     

Identification of KCNJ15 as a Susceptibility Gene in Asian Patients with Type 2 Diabetes Mellitus

Recent advances in genome research have enabled the identification of new genomic variations that are associated with type 2 diabetes mellitus (T2DM). Via fine mapping of SNPs in a candidate region of chromosome 21q, the current study identifies potassium inwardly-rectifying channel, subfamily J, member 15 (KCNJ15) as a new T2DM susceptibility gene. KCNJ15 is expressed in the β cell of the pancreas, and a synonymous SNP, rs3746876, in exon 4 (C566T) of this gene, with T allele frequency among control subjects of 3.1%, showed a significant association with T2DM affecting lean individuals in three independent Japanese sample sets (p = 2.5 × 10⁻⁷, odds ratio [OR] = 2.54, 95% confidence interval [CI] = 1.76–3.67) and with unstratified T2DM (p = 6.7 × 10⁻⁶, OR = 1.76, 95% CI = 1.37–2.25). The diabetes risk allele frequency was, however, very low among Europeans in whom no association between this variant and T2DM could be shown. Functional analysis in human embryonic kidney 293 cells demonstrated that the risk allele of the synonymous SNP in exon 4 increased KCNJ15 expression via increased mRNA stability, which resulted in the higher expression of protein as compared to that of the nonrisk allele. We also showed that KCNJ15 is expressed in human pancreatic β cells. In conclusion, we demonstrated a significant association between a synonymous variant in KCNJ15 and T2DM in lean Japanese patients with T2DM, suggesting that KCNJ15 is a previously unreported susceptibility gene for T2DM among Asians.     

��-Catenin Is Genetically and Biologically Associated with Cortical Cataract and Future Alzheimer-Related Structural and Functional Brain Changes

Multiple lines of evidence suggest that specific subtypes of age-related cataract (ARC) and Alzheimer disease (AD) are related etiologically. To identify shared genetic factors for ARC and AD, we estimated co-heritability of quantitative measures of cataract subtypes with AD-related brain MRI traits among 1,249 members of the Framingham Eye Study who had a brain MRI scan approximately ten years after the eye exam. Cortical cataract (CC) was found to be co-heritable with future development of AD and with several MRI traits, especially temporal horn volume (THV, ρ = 0.24, P<10⁻⁴). A genome-wide association study using 187,657 single nucleotide polymorphisms (SNPs) for the bivariate outcome of CC and THV identified genome-wide significant association with CTNND2 SNPs rs17183619, rs13155993 and rs13170756 (P<2.6×10⁻⁷). These SNPs were also significantly associated with bivariate outcomes of CC and scores on several highly heritable neuropsychological tests (5.7×10⁻⁹≤P<3.7×10⁻⁶). Statistical interaction was demonstrated between rs17183619 and APP SNP rs2096488 on CC (P = 0.0015) and CC-THV (P = 0.038). A rare CTNND2 missense mutation (G810R) 249 base pairs from rs17183619 altered δ-catenin localization and increased secreted amyloid-β_1–42 in neuronal cell culture. Immunohistopathological analysis of lens tissue obtained from two autopsy-confirmed AD subjects and two non-AD controls revealed elevated expression of δ-catenin in epithelial and cortical regions of lenses from AD subjects compared to controls. Our findings suggest that genetic variation in delta catenin may underlie both cortical lens opacities in mid-life and subsequent MRI and cognitive changes that presage the development of AD.     

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.     

Genetic determinants of serum testosterone concentrations in men

Testosterone concentrations in men are associated with cardiovascular morbidity, osteoporosis, and mortality and are affected by age, smoking, and obesity. Because of serum testosterone''s high heritability, we performed a meta-analysis of genome-wide association data in 8,938 men from seven cohorts and followed up the genome-wide significant findings in one in silico (n = 871) and two de novo replication cohorts (n = 4,620) to identify genetic loci significantly associated with serum testosterone concentration in men. All these loci were also associated with low serum testosterone concentration defined as <300 ng/dl. Two single-nucleotide polymorphisms at the sex hormone-binding globulin (SHBG) locus (17p13-p12) were identified as independently associated with serum testosterone concentration (rs12150660, p = 1.2×10⁻⁴¹ and rs6258, p = 2.3×10⁻²²). Subjects with ≥3 risk alleles of these variants had 6.5-fold higher risk of having low serum testosterone than subjects with no risk allele. The rs5934505 polymorphism near FAM9B on the X chromosome was also associated with testosterone concentrations (p = 5.6×10⁻¹⁶). The rs6258 polymorphism in exon 4 of SHBG affected SHBG''s affinity for binding testosterone and the measured free testosterone fraction (p<0.01). Genetic variants in the SHBG locus and on the X chromosome are associated with a substantial variation in testosterone concentrations and increased risk of low testosterone. rs6258 is the first reported SHBG polymorphism, which affects testosterone binding to SHBG and the free testosterone fraction and could therefore influence the calculation of free testosterone using law-of-mass-action equation.     

Quantitative Trait Loci for CD4:CD8 Lymphocyte Ratio Are Associated with Risk of Type 1 Diabetes and HIV-1 Immune Control

Abnormal expansion or depletion of particular lymphocyte subsets is associated with clinical manifestations such as HIV progression to AIDS and autoimmune disease. We sought to identify genetic predictors of lymphocyte levels and reasoned that these may play a role in immune-related diseases. We tested 2.3 million variants for association with five lymphocyte subsets, measured in 2538 individuals from the general population, including CD4+ T cells, CD8+ T cells, CD56+ natural killer (NK) cells, and the derived measure CD4:CD8 ratio. We identified two regions of strong association. The first was located in the major histocompatibility complex (MHC), with multiple SNPs strongly associated with CD4:CD8 ratio (rs2524054, p = 2.1 × 10⁻²⁸). The second region was centered within a cluster of genes from the Schlafen family and was associated with NK cell levels (rs1838149, p = 6.1 × 10⁻¹⁴). The MHC association with CD4:CD8 replicated convincingly (p = 1.4 × 10⁻⁹) in an independent panel of 988 individuals. Conditional analyses indicate that there are two major independent quantitative trait loci (QTL) in the MHC region that regulate CD4:CD8 ratio: one is located in the class I cluster and influences CD8 levels, whereas the second is located in the class II cluster and regulates CD4 levels. Jointly, both QTL explained 8% of the variance in CD4:CD8 ratio. The class I variants are also strongly associated with durable host control of HIV, and class II variants are associated with type-1 diabetes, suggesting that genetic variation at the MHC may predispose one to immune-related diseases partly through disregulation of T cell homeostasis.