Gene-Gene Interactions Lead to Higher Risk for Development of Type 2 Diabetes in an Ashkenazi Jewish Population

Background

Evidence has accumulated that multiple genetic and environmental factors play important roles in determining susceptibility to type 2 diabetes (T2D). Although variants from candidate genes have become prime targets for genetic analysis, few studies have considered their interplay. Our goal was to evaluate interactions among SNPs within genes frequently identified as associated with T2D.

Methods/Principal Findings

Logistic regression was used to study interactions among 4 SNPs, one each from HNF4A[rs1884613], TCF7L2[rs12255372], WFS1[rs10010131], and KCNJ11[rs5219] in a case-control Ashkenazi sample of 974 diabetic subjects and 896 controls. Nonparametric multifactor dimensionality reduction (MDR) and generalized MDR (GMDR) were used to confirm findings from the logistic regression analysis. HNF4A and WFS1 SNPs were associated with T2D in logistic regression analyses [P<0.0001, P<0.0002, respectively]. Interaction between these SNPs were also strong using parametric or nonparametric methods: the unadjusted odds of being affected with T2D was 3 times greater in subjects with the HNF4A and WFS1 risk alleles than those without either (95% CI = [1.7–5.3]; P≤0.0001). Although the univariate association between the TCF7L2 SNP and T2D was relatively modest [P = 0.02], when paired with the HNF4A SNP, the OR for subjects with risk alleles in both SNPs was 2.4 [95% CI = 1.7–3.4; P≤0.0001]. The KCNJ11 variant reached significance only when paired with either the HNF4A or WFSI SNPs: unadjusted ORs were 2.0 [95% CI = 1.4–2.8; P≤0.0001] and 2.3 [95% CI = 1.2-4.4; P≤0.0001], respectively. MDR and GMDR results were consistent with the parametric findings.

Conclusions

These results provide evidence of strong independent associations between T2D and SNPs in HNF4A and WFS1 and their interaction in our Ashkenazi sample. We also observed an interaction in the nonparametric analysis between the HNF4A and KCNJ11 SNPs (P≤0.001), demonstrating that an independently non-significant variant may interact with another variant resulting in an increased disease risk.     

A Genome-Wide Association Study Reveals Variants in ARL15 that Influence Adiponectin Levels

The adipocyte-derived protein adiponectin is highly heritable and inversely associated with risk of type 2 diabetes mellitus (T2D) and coronary heart disease (CHD). We meta-analyzed 3 genome-wide association studies for circulating adiponectin levels (n = 8,531) and sought validation of the lead single nucleotide polymorphisms (SNPs) in 5 additional cohorts (n = 6,202). Five SNPs were genome-wide significant in their relationship with adiponectin (P≤5×10⁻⁸). We then tested whether these 5 SNPs were associated with risk of T2D and CHD using a Bonferroni-corrected threshold of P≤0.011 to declare statistical significance for these disease associations. SNPs at the adiponectin-encoding ADIPOQ locus demonstrated the strongest associations with adiponectin levels (P-combined = 9.2×10⁻¹⁹ for lead SNP, rs266717, n = 14,733). A novel variant in the ARL15 (ADP-ribosylation factor-like 15) gene was associated with lower circulating levels of adiponectin (rs4311394-G, P-combined = 2.9×10⁻⁸, n = 14,733). This same risk allele at ARL15 was also associated with a higher risk of CHD (odds ratio [OR] = 1.12, P = 8.5×10⁻⁶, n = 22,421) more nominally, an increased risk of T2D (OR = 1.11, P = 3.2×10⁻³, n = 10,128), and several metabolic traits. Expression studies in humans indicated that ARL15 is well-expressed in skeletal muscle. These findings identify a novel protein, ARL15, which influences circulating adiponectin levels and may impact upon CHD risk.     

Differential effects of MYH9 and APOL1 risk variants on FRMD3 Association with Diabetic ESRD in African Americans

Single nucleotide polymorphisms (SNPs) in MYH9 and APOL1 on chromosome 22 (c22) are powerfully associated with non-diabetic end-stage renal disease (ESRD) in African Americans (AAs). Many AAs diagnosed with type 2 diabetic nephropathy (T2DN) have non-diabetic kidney disease, potentially masking detection of DN genes. Therefore, genome-wide association analyses were performed using the Affymetrix SNP Array 6.0 in 966 AA with T2DN and 1,032 non-diabetic, non-nephropathy (NDNN) controls, with and without adjustment for c22 nephropathy risk variants. No associations were seen between FRMD3 SNPs and T2DN before adjusting for c22 variants. However, logistic regression analysis revealed seven FRMD3 SNPs significantly interacting with MYH9—a finding replicated in 640 additional AA T2DN cases and 683 NDNN controls. Contrasting all 1,592 T2DN cases with all 1,671 NDNN controls, FRMD3 SNPs appeared to interact with the MYH9 E1 haplotype (e.g., rs942280 interaction p-value = 9.3E⁻⁷ additive; odds ratio [OR] 0.67). FRMD3 alleles were associated with increased risk of T2DN only in subjects lacking two MYH9 E1 risk haplotypes (rs942280 OR = 1.28), not in MYH9 E1 risk allele homozygotes (rs942280 OR = 0.80; homogeneity p-value = 4.3E⁻⁴). Effects were weaker stratifying on APOL1. FRMD3 SNPS were associated with T2DN, not type 2 diabetes per se, comparing AAs with T2DN to those with diabetes lacking nephropathy. T2DN-associated FRMD3 SNPs were detectable in AAs only after accounting for MYH9, with differential effects for APOL1. These analyses reveal a role for FRMD3 in AA T2DN susceptibility and accounting for c22 nephropathy risk variants can assist in detecting DN susceptibility genes.     

Association of Levels of Mannose-Binding Lectin and the MBL2 Gene with Type 2 Diabetes and Diabetic Nephropathy

Objective

To investigate the association of Mannose-binding lectin (MBL) and the MBL2 gene with type 2 diabetes and diabetic nephropathy and the influence of MBL2 polymorphisms on serum MBL levels.

Methods

The study population included 675 type 2 diabetic patients with or without nephropathy and 855 normoglycemic controls. The single nucleotide polymorphisms (SNPs) of rs1800450, rs1800451, and rs11003125 of the MBL2 gene were determined by the Multiplex Snapshot method. Serum MBL levels were measured by enzyme-linked immune sorbent assay.

Results

Rs1800450 and rs11003125 SNPs demonstrated strong linkage disequilibrium in the study population (r² = 0.97). The haplotypes constructed from the G allele of rs1800450 and the C allele of rs11003125 increased the risk for type 2 diabetes (OR = 1.2, 95% CI = 1.1–1.4, P = 0.01). For rs1800450, GG and GA genotypes were associated with type 2 diabetes (P = 0.02, 0.01, respectively). For rs11003125, the GC genotype frequency was significantly different between patients and controls (18.1% vs. 24.9%, P = 0.001). Analyses of genotypes and allele frequency distributions among patients with normal UAE, microalbuminuria, and macroalbuminuria showed that there was no obvious evidence of association between the MBL2 gene and diabetic nephropathy. Subjects with the GG genotype of rs1800450 and the CC genotype of rs11003125 had much higher serum MBL levels.

Conclusions

The rs1800450 and rs11003125 SNPs of the MBL2 gene have strong linkage disequilibrium and are associated with type 2 diabetes in the North Chinese Han population. No association was observed between the MBL2 gene and diabetic nephropathy. Subjects with the GG genotype of rs1800450 and the CC genotype of rs11003125 had much higher serum MBL levels. An association between elevated serum MBL and diabetic nephropathy was also observed.     

Genetic Associations with Diabetes: Meta-Analyses of 10 Candidate Polymorphisms

Aims

The goal of our study is to investigate the combined contribution of 10 genetic variants to diabetes susceptibility.

Methods

Bibliographic databases were searched from 1970 to Dec 2012 for studies that reported on genetic association study of diabetes. After a comprehensive filtering procedure, 10 candidate gene variants with informative genotype information were collected for the current meta-anlayses. Using the REVMAN software, odds ratios (ORs) with 95% confidence intervals (CIs) were calculated to evaluate the combined contribution of the selected genetic variants to diabetes.

Results

A total of 37 articles among 37,033 cases and 54,716 controls were involved in the present meta-analyses of 10 genetic variants. Three variants were found to be significantly associated with type 1 diabetes (T1D): NLRP1 rs12150220 (OR = 0.71, 95% CI = 0.55–0.92, P = 0.01), IL2RA rs11594656 (OR = 0.86, 95% CI = 0.82–0.91, P<0.00001), and CLEC16A rs725613 (OR = 0.71, 95% CI = 0.55–0.92, P = 0.01). APOA5 −1131T/C polymorphism was shown to be significantly associated with of type 2 diabetes (T2D, OR = 1.27, 95% CI = 1.03–1.57, P = 0.03). No association with diabetes was showed in the meta-analyses of other six genetic variants, including SLC2A10 rs2335491, ATF6 rs2070150, KLF11 rs35927125, CASQ1 rs2275703, GNB3 C825T, and IL12B 1188A/C.

Conclusion

Our results demonstrated that IL2RA rs11594656 and CLEC16A rs725613 are protective factors of T1D, while NLRP1 rs12150220 and APOA5 −1131T/C are risky factors of T1D and T2D, respectively.     

An Exploratory Study of the Association between KCNB1 rs1051295 and Type 2 Diabetes and Its Related Traits in Chinese Han Population

Since the KCNB1 encoding Kv2.1 channel accounts for the majority of Kv currents modulating insulin secretion by pancreatic islet beta-cells, we postulated that KCNB1 is a plausible candidate gene for genetic variation contributing to the variable compensatory secretory function of beta-cells in type-2 diabetes (T2D). We conducted two studies, a case-control study and a cross-section study, to investigate the association of common single-nucleotide polymorphisms (SNPs) in KCNB1 with T2D and its linking traits. In the case-control study, we first examined the association of 20 tag SNPs of KCNB1 with T2D in a population with 226 T2D patients and non-diabetic subjects (screening study). We then identified the association in an enlarged population of 412 T2D patients and non-diabetic subjects (replication study). In the cross-sectional study, we investigated the linkage between the candidate SNP rs1051295 and T2D by comparing beta-cell function and insulin sensitivity among rs1051295 genotypes in a general population of 1051 subjects at fasting and after glucose loading (oral glucose tolerance tests, OGTT) in 84 fasting glucose impaired subjects, and several T2D-related traits. We found that among the 19 available tag SNPs, only the KCNB1 rs1051295 was associated with T2D (P = 0.027), with the rs1051295 TT genotype associated with an increased risk of T2D compared with genotypes CC (P = 0.009). At fasting, rs1051295 genotype TT was associated with a 9.8% reduction in insulin sensitivity compared to CC (P = 0.008); along with increased plasma triglycerides (TG) levels (TT/CC: P = 0.046) and increased waist/hip (W/H) ratio (TT/CC: P = 0.013; TT/TC: P = 0.002). OGTT confirmed that genotype TT exhibited reduced insulin sensitivity by 16.3% (P = 0.030) compared with genotype TC+CC in a fasting glucose impaired population. The KCNB1 rs1051295 genotype TT in the Chinese Han population is associated with decreased insulin sensitivity and increased plasma TG and W/H ratio, which together contribute to an increased risk for T2D.     

PTPN22 Association in Systemic Lupus Erythematosus (SLE) with Respect to Individual Ancestry and Clinical Sub-Phenotypes

Protein tyrosine phosphatase non-receptor type 22 (PTPN22) is a negative regulator of T-cell activation associated with several autoimmune diseases, including systemic lupus erythematosus (SLE). Missense rs2476601 is associated with SLE in individuals with European ancestry. Since the rs2476601 risk allele frequency differs dramatically across ethnicities, we assessed robustness of PTPN22 association with SLE and its clinical sub-phenotypes across four ethnically diverse populations. Ten SNPs were genotyped in 8220 SLE cases and 7369 controls from in European-Americans (EA), African-Americans (AA), Asians (AS), and Hispanics (HS). We performed imputation-based association followed by conditional analysis to identify independent associations. Significantly associated SNPs were tested for association with SLE clinical sub-phenotypes, including autoantibody profiles. Multiple testing was accounted for by using false discovery rate. We successfully imputed and tested allelic association for 107 SNPs within the PTPN22 region and detected evidence of ethnic-specific associations from EA and HS. In EA, the strongest association was at rs2476601 (P = 4.7×10⁻⁹, OR = 1.40 (95% CI = 1.25–1.56)). Independent association with rs1217414 was also observed in EA, and both SNPs are correlated with increased European ancestry. For HS imputed intronic SNP, rs3765598, predicted to be a cis-eQTL, was associated (P = 0.007, OR = 0.79 and 95% CI = 0.67–0.94). No significant associations were observed in AA or AS. Case-only analysis using lupus-related clinical criteria revealed differences between EA SLE patients positive for moderate to high titers of IgG anti-cardiolipin (aCL IgG >20) versus negative aCL IgG at rs2476601 (P = 0.012, OR = 1.65). Association was reinforced when these cases were compared to controls (P = 2.7×10⁻⁵, OR = 2.11). Our results validate that rs2476601 is the most significantly associated SNP in individuals with European ancestry. Additionally, rs1217414 and rs3765598 may be associated with SLE. Further studies are required to confirm the involvement of rs2476601 with aCL IgG.     

SORL1 Is Genetically Associated with Late-Onset Alzheimer’s Disease in Japanese,Koreans and Caucasians

To discover susceptibility genes of late-onset Alzheimer’s disease (LOAD), we conducted a 3-stage genome-wide association study (GWAS) using three populations: Japanese from the Japanese Genetic Consortium for Alzheimer Disease (JGSCAD), Koreans, and Caucasians from the Alzheimer Disease Genetic Consortium (ADGC). In Stage 1, we evaluated data for 5,877,918 genotyped and imputed SNPs in Japanese cases (n = 1,008) and controls (n = 1,016). Genome-wide significance was observed with 12 SNPs in the APOE region. Seven SNPs from other distinct regions with p-values <2×10⁻⁵ were genotyped in a second Japanese sample (885 cases, 985 controls), and evidence of association was confirmed for one SORL1 SNP (rs3781834, P = 7.33×10⁻⁷ in the combined sample). Subsequent analysis combining results for several SORL1 SNPs in the Japanese, Korean (339 cases, 1,129 controls) and Caucasians (11,840 AD cases, 10,931 controls) revealed genome wide significance with rs11218343 (P = 1.77×10⁻⁹) and rs3781834 (P = 1.04×10⁻⁸). SNPs in previously established AD loci in Caucasians showed strong evidence of association in Japanese including rs3851179 near PICALM (P = 1.71×10⁻⁵) and rs744373 near BIN1 (P = 1.39×10⁻⁴). The associated allele for each of these SNPs was the same as in Caucasians. These data demonstrate for the first time genome-wide significance of LOAD with SORL1 and confirm the role of other known loci for LOAD in Japanese. Our study highlights the importance of examining associations in multiple ethnic populations.     

Familial Young-Onset Diabetes,Pre-Diabetes and Cardiovascular Disease Are Associated with Genetic Variants of DACH1 in Chinese

In Asia, young-onset type 2 diabetes (YOD) is characterized by obesity and increased risk for cardiovascular disease (CVD). In a genome-wide association study (GWAS) of 99 Chinese obese subjects with familial YOD diagnosed before 40-year-old and 101 controls, the T allele of rs1408888 in intron 1 of DACH1(Dachshund homolog 1) was associated with an odds ratio (OR) of 2.49(95% confidence intervals:1.57–3.96, P = 8.4×10⁻⁵). Amongst these subjects, we found reduced expression of DACH1 in peripheral blood mononuclear cells (PBMC) from 63 cases compared to 65 controls (P = 0.02). In a random cohort of 1468 cases and 1485 controls, amongst top 19 SNPs from GWAS, rs1408888 was associated with type 2 diabetes with a global P value of 0.0176 and confirmation in a multiethnic Asian case-control cohort (7370/7802) with an OR of 1.07(1.02–1.12, P_meta = 0.012). In 599 Chinese non-diabetic subjects, rs1408888 was linearly associated with systolic blood pressure and insulin resistance. In a case-control cohort (n = 953/953), rs1408888 was associated with an OR of 1.54(1.07–2.22, P = 0.019) for CVD in type 2 diabetes. In an autopsy series of 173 non-diabetic cases, TT genotype of rs1408888 was associated with an OR of 3.31(1.19–9.19, P = 0.0214) and 3.27(1.25–11.07, P = 0.0184) for coronary heart disease (CHD) and coronary arteriosclerosis. Bioinformatics analysis revealed that rs1408888 lies within regulatory elements of DACH1 implicated in islet development and insulin secretion. The T allele of rs1408888 of DACH1 was associated with YOD, prediabetes and CVD in Chinese.     

The A Allele of the rs1990760 Polymorphism in the IFIH1 Gene Is Associated with Protection for Arterial Hypertension in Type 1 Diabetic Patients and with Expression of This Gene in Human Mononuclear Cells

Background

The rs1990760 polymorphism of interferon induced with helicase C domain 1 (IFIH1) has been associated with type 1 diabetes mellitus (T1DM). Here, we investigated whether this polymorphism is associated with T1DM or its clinical characteristics in a Brazilian population, and if IFIH1 gene expression in mononuclear cells from T1DM patients differs according to the genotypes of this polymorphism. A meta-analysis was also conducted to evaluate if the rs1990760 polymorphism is associated with T1DM.

Methods

Frequencies of the rs1990760 polymorphism were analyzed in 527 T1DM patients and in 517 healthy subjects. IFIH1 gene expressions according to genotypes were measured in a sub-sample of 26 T1DM patients by quantitative real-time PCR.

Results

Our data show the association of the A allele with risk to T1DM under a dominant model of inheritance [odds ratio (OR) = 1.421, P = 0.037], adjusting for ethnicity. The meta-analysis revealed significant association between the rs199760A allele and risk for T1DM for all analyzed inheritance models. Surprisingly, T1DM patients carrying the A allele showed lower levels of systolic (P = 0.001) and diastolic (P = 1×10⁻¹⁰) blood pressures as compared to G/G carriers. Furthermore, the A/A genotype seems to be associated with protection to arterial hypertension (AH) after adjustment for covariates (OR = 0.339, P = 0.019). IFIH1 gene expression in mononuclear cells from 26 T1DM patients did not differ among genotypes (P = 0.274). Nevertheless, IFIH1 gene expression was increased in mononuclear cells from T1DM patients with AH as compared with T1DM patients without AH [6.7 (1.7–2.0) vs. 1.8 (1.3–7.1) arbitrary units; P = 0.036]. The association with blood pressures and AH was not observed in patients with type 2 diabetes mellitus.

Conclusions

Our results indicate that the rs1990760 polymorphism is associated with T1DM. Interestingly, the rs1990760 A allele seems to be associated with protection for AH in T1DM patients. Further studies are needed to confirm the association with AH.     

DNA Methyl Transferase (DNMT) Gene Polymorphisms Could Be a Primary Event in Epigenetic Susceptibility to Schizophrenia

DNA methylation has been implicated in the etiopathology of various complex disorders. DNA methyltransferases are involved in maintaining and establishing new methylation patterns. The aim of the present study was to investigate the inherent genetic variations within DNA methyltransferase genes in predisposing to susceptibility to schizophrenia. We screened for polymorphisms in DNA methyltransferases, DNMT1, DNMT3A, DNMT3B and DNMT3L in 330 schizophrenia patients and 302 healthy controls for association with Schizophrenia in south Indian population. These polymorphisms were also tested for subgroup analysis with patient''s gender, age of onset and family history. DNMT1 rs2114724 (genotype P = .004, allele P = 0.022) and rs2228611 (genotype P = 0.004, allele P = 0.022) were found to be significantly associated at genotypic and allelic level with Schizophrenia in South Indian population. DNMT3B rs2424932 genotype (P = 0.023) and allele (P = 0.0063) increased the risk of developing schizophrenia in males but not in females. DNMT3B rs1569686 (genotype P = 0.027, allele P = 0.033) was found to be associated with early onset of schizophrenia and also with family history and early onset (genotype P = 0.009). DNMT3L rs2070565 (genotype P = 0.007, allele P = 0.0026) confers an increased risk of developing schizophrenia at an early age in individuals with family history. In-silico prediction indicated functional relevance of these SNPs in regulating the gene. These observations might be crucial in addressing and understanding the genetic control of methylation level differences from ethnic viewpoint. Functional significance of genotype variations within the DNMTs indeed suggest that the genetic nature of methyltransferases should be considered while addressing epigenetic events mediated by methylation in Schizophrenia.     

IL2RA/CD25 Gene Polymorphisms: Uneven Association with Multiple Sclerosis (MS) and Type 1 Diabetes (T1D)

Background

IL-2 receptor (IL2R) alpha is the specific component of the high affinity IL2R system involved in the immune response and in the control of autoimmunity.

Methods and Results

Here we perform a replication and fine mapping of the IL2RA gene region analyzing 3 SNPs previously associated with multiple sclerosis (MS) and 5 SNPs associated with type 1 diabetes (T1D) in a collection of 798 MS patients and 927 matched Caucasian controls from the south of Spain. We observed association with MS in 6 of 8 SNPs. The rs1570538, at the 3′- UTR extreme of the gene, previously reported to have a weak association with MS, is replicated here (P = 0.032). The most associated T1D SNP (rs41295061) was not associated with MS in the present study. However, the rs35285258, belonging to another independent group of SNPs associated with T1D, showed the maximal association in this study but different risk allele. We replicated the association of only one (rs2104286) of the two IL2RA SNPs identified in the recently performed genome-wide association study of MS.

Conclusions

These findings confirm and extend the association of this gene with MS and reveal a genetic heterogeneity of the associated polymorphisms and risk alleles between MS and T1D suggesting different immunopathological roles of IL2RA in these two diseases.     

A Large Scale Gene-Centric Association Study of Lung Function in Newly-Hired Female Cotton Textile Workers with Endotoxin Exposure

Background

Occupational exposure to endotoxin is associated with decrements in pulmonary function, but how much variation in this association is explained by genetic variants is not well understood.

Objective

We aimed to identify single nucleotide polymorphisms (SNPs) that are associated with the rate of forced expiratory volume in one second (FEV₁) decline by a large scale genetic association study in newly-hired healthy young female cotton textile workers.

Methods

DNA samples were genotyped using the Illumina Human CVD BeadChip. Change rate in FEV₁ was modeled as a function of each SNP genotype in linear regression model with covariate adjustment. We controlled the type 1 error in study-wide level by permutation method. The false discovery rate (FDR) and the family-wise error rate (FWER) were set to be 0.10 and 0.15 respectively.

Results

Two SNPs were found to be significant (P<6.29×10⁻⁵), including rs1910047 (P = 3.07×10⁻⁵, FDR = 0.0778) and rs9469089 (P = 6.19×10⁻⁵, FDR = 0.0967), as well as other eight suggestive (P<5×10⁻⁴) associated SNPs. Gene-gene and gene-environment interactions were also observed, such as rs1910047 and rs1049970 (P = 0.0418, FDR = 0.0895); rs9469089 and age (P = 0.0161, FDR = 0.0264). Genetic risk score analysis showed that the more risk loci the subjects carried, the larger the rate of FEV₁ decline occurred (P _trend = 3.01×10⁻¹⁸). However, the association was different among age subgroups (P = 7.11×10⁻⁶) and endotoxin subgroups (P = 1.08×10⁻²). Functional network analysis illustrates potential biological connections of all interacted genes.

Conclusions

Genetic variants together with environmental factors interact to affect the rate of FEV₁ decline in cotton textile workers.     

Replication of Association between Schizophrenia and Chromosome 6p21-6p22.1 Polymorphisms in Chinese Han Population

Chromosome 6p21-p22.1, spanning the extended major histocompatibility complex (MHC) region, is a highly polymorphic, gene-dense region. It has been identified as a susceptibility locus of schizophrenia in Europeans, Japanese, and Chinese. In our previous two-stage genome-wide association study (GWAS), polymorphisms of zinc finger with KRAB and SCAN domains 4 (ZKSCAN4), nuclear factor-κB-activating protein-like (NKAPL), and piggyBac transposable element derived 1 (PGBD1), localized to chromosome 6p21-p22.1, were strongly associated with schizophrenia. To further investigate the association between polymorphisms at this locus and schizophrenia in the Chinese Han population, we selected eight other single-nucleotide polymorphisms (SNPs) distributed in or near these genes for a case-control association study in an independent sample of 902 cases and 1,091 healthy controls in an attempt to replicate the GWAS results. Four of these eight SNPs (rs12214383, rs1150724, rs3800324, and rs1997660) displayed a nominal difference in allele frequencies between the case and control groups. The association between two of these SNPs and schizophrenia were significant even after Bonferroni correction (rs12000: allele A>G, P = 2.50E-04, odds ratio [OR] = 1.27, 95% confidence interval [CI] = 1.12–1.45; rs1150722: allele C>T, P = 4.28E-05, OR = 0.55, 95% CI = 0.41–0.73). Haplotype ATTGACGC, comprising these eight SNPs (rs2235359, rs2185955, rs12214383, rs12000, rs1150724, rs1150722, rs3800324, and rs1997660), was significantly associated with schizophrenia (P = 6.60E-05). We also performed a combined study of this replication sample and the first-stage GWAS sample. The combined study revealed that rs12000 and rs1150722 were still strongly associated with schizophrenia (rs12000: allele G>A, P _combined  = 0.0019, OR = 0.81; rs1150722: allele G>A, P _combined  = 3.00E-04, OR = 0.61). These results support our findings that locus 6p21-p22.1 is significantly associated with schizophrenia in the Chinese Han population and encourage further studies of the functions of these genetic factors.     

An ACACB Variant Implicated in Diabetic Nephropathy Associates with Body Mass Index and Gene Expression in Obese Subjects

Acetyl coenzyme A carboxylase B gene (ACACB) single nucleotide polymorphism (SNP) rs2268388 is reproducibly associated with type 2 diabetes (T2DM)-associated nephropathy (DN). ACACB knock-out mice are also protected from obesity. This study assessed relationships between rs2268388, body mass index (BMI) and gene expression in multiple populations, with and without T2DM. Among subjects without T2DM, rs2268388 DN risk allele (T) associated with higher BMI in Pima Indian children (n = 2021; p-additive = 0.029) and African Americans (AAs) (n = 177; p-additive = 0.05), with a trend in European Americans (EAs) (n = 512; p-additive = 0.09), but not Germans (n = 858; p-additive = 0.765). Association with BMI was seen in a meta-analysis including all non-T2DM subjects (n = 3568; p-additive = 0.02). Among subjects with T2DM, rs2268388 was not associated with BMI in Japanese (n = 2912) or EAs (n = 1149); however, the T allele associated with higher BMI in the subset with BMI≥30 kg/m² (n = 568 EAs; p-additive = 0.049, n = 196 Japanese; p-additive = 0.049). Association with BMI was strengthened in a T2DM meta-analysis that included an additional 756 AAs (p-additive = 0.080) and 48 Hong Kong Chinese (p-additive = 0.81) with BMI≥30 kg/m² (n = 1575; p-additive = 0.0033). The effect of rs2268388 on gene expression revealed that the T risk allele associated with higher ACACB messenger levels in adipose tissue (41 EAs and 20 AAs with BMI>30 kg/m²; p-additive = 0.018) and ACACB protein levels in the liver tissue (mixed model p-additive = 0.03, in 25 EA bariatric surgery patients with BMI>30 kg/m² for 75 exams). The T allele also associated with higher hepatic triglyceride levels. These data support a role for ACACB in obesity and potential roles for altered lipid metabolism in susceptibility to DN.     

Association between Human Prothrombin Variant (T165M) and Kidney Stone Disease

We previously reported the association between prothrombin (F2), encoding a stone inhibitor protein - urinary prothrombin fragment 1 (UPTF1), and the risk of kidney stone disease in Northeastern Thai patients. To identify specific F2 variation responsible for the kidney stone risk, we conducted sequencing analysis of this gene in a group of the patients with kidney stone disease. Five intronic SNPs (rs2070850, rs2070852, rs1799867, rs2282687, and rs3136516) and one exonic non-synonymous single nucleotide polymorphism (nsSNP; rs5896) were found. The five intronic SNPs have no functional change as predicted by computer programs while the nsSNP rs5896 (c.494 C>T) located in exon 6 results in a substitution of threonine (T) by methionine (M) at the position 165 (T165M). The nsSNP rs5896 was subsequently genotyped in 209 patients and 216 control subjects. Genotypic and allelic frequencies of this nsSNP were analyzed for their association with kidney stone disease. The frequency of CC genotype of rs5896 was significantly lower in the patient group (13.4%) than that in the control group (22.2%) (P = 0.017, OR 0.54, 95% CI 0.32–0.90), and the frequency of C allele was significantly lower in the patient group (36.1%) than that in the control group (45.6%) (P = 0.005, OR 0.68, 95% CI 0.51–0.89). The significant differences of genotype and allele frequencies were maintained only in the female group (P = 0.033 and 0.003, respectively). The effect of amino-acid change on UPTF1 structure was also examined by homologous modeling and in silico mutagenesis. T165 is conserved and T165M substitution will affect hydrogen bond formation with E180. In conclusion, our results indicate that prothrombin variant (T165M) is associated with kidney stone risk in the Northeastern Thai female patients.     

Case-Control Analysis of SNPs in GLUT4, RBP4 and STRA6: Association of SNPs in STRA6 with Type 2 Diabetes in a South Indian Population

Background

The inverse relationship between GLUT4 and RBP4 expression is known to play a role in the pathogenesis of type 2 diabetes. Elevated levels of RBP4 were shown to cause insulin resistance in muscles and liver. Identification of STRA6 as a cell surface receptor for RBP4 provides further link in this axis and hence we analyzed SNPs in these three genes for association with type 2 diabetes in a South Indian population.

Methodology/Principal Findings

Selected SNPs in the three genes were analyzed in a total of 2002 individuals belonging to Dravidian ethnicity, South India, by Tetra Primer ARMS PCR or RFLP PCR. Allele frequencies and genotype distribution were calculated in cases and controls and were analyzed for association by Chi-squared test and Logistic regression. Haplotype analysis was carried out for each gene by including all the markers in a single block. We observed a significant association of three SNPs, rs974456, rs736118, and rs4886578 in STRA6 with type 2 diabetes (P = 0.001, OR 0.79[0.69–0.91], P = 0.003, OR 0.81[0.71–0.93], and P = 0.001, OR 0.74[0.62–0.89] respectively). None of the SNPs in RBP4 and GLUT4 showed any association with type 2 diabetes. Haplotype analysis revealed that two common haplotypes H1 (111, P = 0.001, OR 1.23[1.08–1.40]) and H2 (222, P = 0.002 OR 0.73[0.59–0.89]) in STRA6, H6 (2121, P = 0.006, OR 1.69[1.51–2.48]) in RBP4 and H4 (2121, P = 0.01 OR 1.41[1.07–1.85]) in GLUT4 were associated with type 2 diabetes.

Conclusion

SNPs in STRA6, gene coding the cell surface receptor for RBP4, were significantly associated with type 2 diabetes and further genetic and functional studies are required to understand and ascertain its role in the manifestation of type 2 diabetes.     

Genetic Determinants of Long-Term Changes in Blood Lipid Concentrations: 10-Year Follow-Up of the GLACIER Study

Recent genome-wide meta-analyses identified 157 loci associated with cross-sectional lipid traits. Here we tested whether these loci associate (singly and in trait-specific genetic risk scores [GRS]) with longitudinal changes in total cholesterol (TC) and triglyceride (TG) levels in a population-based prospective cohort from Northern Sweden (the GLACIER Study). We sought replication in a southern Swedish cohort (the MDC Study; N = 2,943). GLACIER Study participants (N = 6,064) were genotyped with the MetaboChip array. Up to 3,495 participants had 10-yr follow-up data available in the GLACIER Study. The TC- and TG-specific GRSs were strongly associated with change in lipid levels (β = 0.02 mmol/l per effect allele per decade follow-up, P = 2.0×10⁻¹¹ for TC; β = 0.02 mmol/l per effect allele per decade follow-up, P = 5.0×10⁻⁵ for TG). In individual SNP analysis, one TC locus, apolipoprotein E (APOE) rs4420638 (β = 0.12 mmol/l per effect allele per decade follow-up, P = 2.0×10⁻⁵), and two TG loci, tribbles pseudokinase 1 (TRIB1) rs2954029 (β = 0.09 mmol/l per effect allele per decade follow-up, P = 5.1×10⁻⁴) and apolipoprotein A-I (APOA1) rs6589564 (β = 0.31 mmol/l per effect allele per decade follow-up, P = 1.4×10⁻⁸), remained significantly associated with longitudinal changes for the respective traits after correction for multiple testing. An additional 12 loci were nominally associated with TC or TG changes. In replication analyses, the APOE rs4420638, TRIB1 rs2954029, and APOA1 rs6589564 associations were confirmed (P≤0.001). In summary, trait-specific GRSs are robustly associated with 10-yr changes in lipid levels and three individual SNPs were strongly associated with 10-yr changes in lipid levels.