Genetic association for renal traits among participants of African ancestry reveals new loci for renal function

Chronic kidney disease (CKD) is an increasing global public health concern, particularly among populations of African ancestry. We performed an interrogation of known renal loci, genome-wide association (GWA), and IBC candidate-gene SNP association analyses in African Americans from the CARe Renal Consortium. In up to 8,110 participants, we performed meta-analyses of GWA and IBC array data for estimated glomerular filtration rate (eGFR), CKD (eGFR <60 mL/min/1.73 m²), urinary albumin-to-creatinine ratio (UACR), and microalbuminuria (UACR >30 mg/g) and interrogated the 250 kb flanking region around 24 SNPs previously identified in European Ancestry renal GWAS analyses. Findings were replicated in up to 4,358 African Americans. To assess function, individually identified genes were knocked down in zebrafish embryos by morpholino antisense oligonucleotides. Expression of kidney-specific genes was assessed by in situ hybridization, and glomerular filtration was evaluated by dextran clearance. Overall, 23 of 24 previously identified SNPs had direction-consistent associations with eGFR in African Americans, 2 of which achieved nominal significance (UMOD, PIP5K1B). Interrogation of the flanking regions uncovered 24 new index SNPs in African Americans, 12 of which were replicated (UMOD, ANXA9, GCKR, TFDP2, DAB2, VEGFA, ATXN2, GATM, SLC22A2, TMEM60, SLC6A13, and BCAS3). In addition, we identified 3 suggestive loci at DOK6 (p-value = 5.3×10⁻⁷) and FNDC1 (p-value = 3.0×10⁻⁷) for UACR, and KCNQ1 with eGFR (p = 3.6×10⁻⁶). Morpholino knockdown of kcnq1 in the zebrafish resulted in abnormal kidney development and filtration capacity. We identified several SNPs in association with eGFR in African Ancestry individuals, as well as 3 suggestive loci for UACR and eGFR. Functional genetic studies support a role for kcnq1 in glomerular development in zebrafish.     

Genetic variants at chromosomes 2q35, 5p12, 6q25.1, 10q26.13, and 16q12.1 influence the risk of breast cancer in men

Male breast cancer accounts for approximately 1% of all breast cancer. To date, risk factors for male breast cancer are poorly defined, but certain risk factors and genetic features appear common to both male and female breast cancer. Genome-wide association studies (GWAS) have recently identified common single nucleotide polymorphisms (SNPs) that influence female breast cancer risk; 12 of these have been independently replicated. To examine if these variants contribute to male breast cancer risk, we genotyped 433 male breast cancer cases and 1,569 controls. Five SNPs showed a statistically significant association with male breast cancer: rs13387042 (2q35) (odds ratio (OR)  = 1.30, p = 7.98×10⁻⁴), rs10941679 (5p12) (OR = 1.26, p = 0.007), rs9383938 (6q25.1) (OR = 1.39, p = 0.004), rs2981579 (FGFR2) (OR = 1.18, p = 0.03), and rs3803662 (TOX3) (OR = 1.48, p = 4.04×10⁻⁶). Comparing the ORs for male breast cancer with the published ORs for female breast cancer, three SNPs—rs13387042 (2q35), rs3803662 (TOX3), and rs6504950 (COX11)—showed significant differences in ORs (p<0.05) between sexes. Breast cancer is a heterogeneous disease; the relative risks associated with loci identified to date show subtype and, based on these data, gender specificity. Additional studies of well-defined patient subgroups could provide further insight into the biological basis of breast cancer development.     

Common variants at 9p21 and 8q22 are associated with increased susceptibility to optic nerve degeneration in glaucoma

Optic nerve degeneration caused by glaucoma is a leading cause of blindness worldwide. Patients affected by the normal-pressure form of glaucoma are more likely to harbor risk alleles for glaucoma-related optic nerve disease. We have performed a meta-analysis of two independent genome-wide association studies for primary open angle glaucoma (POAG) followed by a normal-pressure glaucoma (NPG, defined by intraocular pressure (IOP) less than 22 mmHg) subgroup analysis. The single-nucleotide polymorphisms that showed the most significant associations were tested for association with a second form of glaucoma, exfoliation-syndrome glaucoma. The overall meta-analysis of the GLAUGEN and NEIGHBOR dataset results (3,146 cases and 3,487 controls) identified significant associations between two loci and POAG: the CDKN2BAS region on 9p21 (rs2157719 [G], OR = 0.69 [95%CI 0.63–0.75], p = 1.86×10⁻¹⁸), and the SIX1/SIX6 region on chromosome 14q23 (rs10483727 [A], OR = 1.32 [95%CI 1.21–1.43], p = 3.87×10⁻¹¹). In sub-group analysis two loci were significantly associated with NPG: 9p21 containing the CDKN2BAS gene (rs2157719 [G], OR = 0.58 [95% CI 0.50–0.67], p = 1.17×10⁻¹²) and a probable regulatory region on 8q22 (rs284489 [G], OR = 0.62 [95% CI 0.53–0.72], p = 8.88×10⁻¹⁰). Both NPG loci were also nominally associated with a second type of glaucoma, exfoliation syndrome glaucoma (rs2157719 [G], OR = 0.59 [95% CI 0.41–0.87], p = 0.004 and rs284489 [G], OR = 0.76 [95% CI 0.54–1.06], p = 0.021), suggesting that these loci might contribute more generally to optic nerve degeneration in glaucoma. Because both loci influence transforming growth factor beta (TGF-beta) signaling, we performed a genomic pathway analysis that showed an association between the TGF-beta pathway and NPG (permuted p = 0.009). These results suggest that neuro-protective therapies targeting TGF-beta signaling could be effective for multiple forms of glaucoma.     

TRAIL/TRAIL receptor system and susceptibility to multiple sclerosis

The TNF-related apoptosis inducing ligand (TRAIL)/TRAIL receptor system participates in crucial steps in immune cell activation or differentiation. It is able to inhibit proliferation and activation of T cells and to induce apoptosis of neurons and oligodendrocytes, and seems to be implicated in autoimmune diseases. Thus, TRAIL and TRAIL receptor genes are potential candidates for involvement in susceptibility to multiple sclerosis (MS). To test whether single-nucleotide polymorphisms (SNPs) in the human genes encoding TRAIL, TRAILR-1, TRAILR-2, TRAILR-3 and TRAILR-4 are associated with MS susceptibility, we performed a candidate gene case-control study in the Spanish population. 59 SNPs in the TRAIL and TRAIL receptor genes were analysed in 628 MS patients and 660 controls, and validated in an additional cohort of 295 MS patients and 233 controls. Despite none of the SNPs withstood the highly conservative Bonferroni correction, three SNPs showing uncorrected p values<0.05 were successfully replicated: rs4894559 in TRAIL gene, p = 9.8×10⁻⁴, OR = 1.34; rs4872077, in TRAILR-1 gene, p = 0.005, OR = 1.72; and rs1001793 in TRAILR-2 gene, p = 0.012, OR = 0.84. The combination of the alleles G/T/A in these SNPs appears to be associated with a reduced risk of developing MS (p = 2.12×10⁻⁵, OR = 0.59). These results suggest that genes of the TRAIL/TRAIL receptor system exerts a genetic influence on MS.     

Genetic determinants of cardiovascular events among women with migraine: a genome-wide association study

Background

Migraine is associated with an increased risk for cardiovascular disease (CVD). Both migraine and CVD are highly heritable. However, the genetic liability for CVD among migraineurs is unclear.

Methods

We performed a genome-wide association study for incident CVD events during 12 years of follow-up among 5,122 migraineurs participating in the population-based Women''s Genome Health Study. Migraine was self-reported and CVD events were confirmed after medical records review. We calculated odds ratios (OR) and 95% confidence intervals (CI) and considered a genome-wide p-value <5×10⁻⁸ as significant.

Results

Among the 5,122 women with migraine 164 incident CVD events occurred during follow-up. No SNP was associated with major CVD, ischemic stroke, myocardial infarction, or CVD death at the genome-wide level; however, five SNPs showed association with p<5×10⁻⁶. Among migraineurs with aura rs7698623 in MEPE (OR = 6.37; 95% CI 3.15–12.90; p = 2.7×10⁻⁷) and rs4975709 in IRX4 (OR = 5.06; 95% CI 2.66–9.62; p = 7.7×10⁻⁷) appeared to be associated with ischemic stroke, rs2143678 located close to MDF1 with major CVD (OR = 3.05; 95% CI 1.98–4.69; p = 4.3×10⁻⁷), and the intergenic rs1406961 with CVD death (OR = 12.33; 95% CI 4.62–32.87; p = 5.2×10⁻⁷). Further, rs1047964 in BACE1 appeared to be associated with CVD death among women with any migraine (OR = 4.67; 95% CI 2.53–8.62; p = 8.0×10⁻⁷).

Conclusion

Our results provide some suggestion for an association of five SNPs with CVD events among women with migraine; none of the results was genome-wide significant. Four associations appeared among migraineurs with aura, two of those with ischemic stroke. Although our population is among the largest with migraine and incident CVD information, these results must be treated with caution, given the limited number of CVD events among women with migraine and the low minor allele frequencies for three of the SNPs. Our results await independent replication and should be considered hypothesis generating for future research.     

A genome-wide scan of Ashkenazi Jewish Crohn's disease suggests novel susceptibility loci

Crohn''s disease (CD) is a complex disorder resulting from the interaction of intestinal microbiota with the host immune system in genetically susceptible individuals. The largest meta-analysis of genome-wide association to date identified 71 CD–susceptibility loci in individuals of European ancestry. An important epidemiological feature of CD is that it is 2–4 times more prevalent among individuals of Ashkenazi Jewish (AJ) descent compared to non-Jewish Europeans (NJ). To explore genetic variation associated with CD in AJs, we conducted a genome-wide association study (GWAS) by combining raw genotype data across 10 AJ cohorts consisting of 907 cases and 2,345 controls in the discovery stage, followed up by a replication study in 971 cases and 2,124 controls. We confirmed genome-wide significant associations of 9 known CD loci in AJs and replicated 3 additional loci with strong signal (p<5×10⁻⁶). Novel signals detected among AJs were mapped to chromosomes 5q21.1 (rs7705924, combined p = 2×10⁻⁸; combined odds ratio OR = 1.48), 2p15 (rs6545946, p = 7×10⁻⁹; OR = 1.16), 8q21.11 (rs12677663, p = 2×10⁻⁸; OR = 1.15), 10q26.3 (rs10734105, p = 3×10⁻⁸; OR = 1.27), and 11q12.1 (rs11229030, p = 8×10⁻⁹; OR = 1.15), implicating biologically plausible candidate genes, including RPL7, CPAMD8, PRG2, and PRG3. In all, the 16 replicated and newly discovered loci, in addition to the three coding NOD2 variants, accounted for 11.2% of the total genetic variance for CD risk in the AJ population. This study demonstrates the complementary value of genetic studies in the Ashkenazim.     

Interleukin 12B (IL12B) genetic variation and pulmonary tuberculosis: a study of cohorts from The Gambia, Guinea-Bissau, United States and Argentina

We examined whether polymorphisms in interleukin-12B (IL12B) associate with susceptibility to pulmonary tuberculosis (PTB) in two West African populations (from The Gambia and Guinea-Bissau) and in two independent populations from North and South America. Nine polymorphisms (seven SNPs, one insertion/deletion, one microsatellite) were analyzed in 321 PTB cases and 346 controls from Guinea-Bissau and 280 PTB cases and 286 controls from The Gambia. For replication we studied 281 case and 179 control African-American samples and 221 cases and 144 controls of European ancestry from the US and Argentina. First-stage single locus analyses revealed signals of association at IL12B 3′ UTR SNP rs3212227 (unadjusted allelic p = 0.04; additive genotypic p = 0.05, OR = 0.78, 95% CI [0.61–0.99]) in Guinea-Bissau and rs11574790 (unadjusted allelic p = 0.05; additive genotypic p = 0.05, OR = 0.76, 95% CI [0.58–1.00]) in The Gambia. Association of rs3212227 was then replicated in African-Americans (rs3212227 allelic p = 0.002; additive genotypic p = 0.05, OR = 0.78, 95% CI [0.61–1.00]); most importantly, in the African-American cohort, multiple significant signals of association (seven of the nine polymorphisms tested) were detected throughout the gene. These data suggest that genetic variation in IL12B, a highly relevant candidate gene, is a risk factor for PTB in populations of African ancestry, although further studies will be required to confirm this association and identify the precise mechanism underlying it.     

Association of six single nucleotide polymorphisms with gestational diabetes mellitus in a Chinese population

Background

To investigate whether the candidate genes that confer susceptibility to type 2 diabetes mellitus are also correlated with gestational diabetes mellitus (GDM) in pregnant Chinese women.

Methodology/Principal Findings

In this study, 1764 unrelated pregnant women were recruited, of which 725 women had GDM and 1039 served as controls. Six single nucleotide polymorphisms (rs7754840 in CDKAL1, rs391300 in SRR, rs2383208 in CDKN2A/2B, rs4402960 in IGF2BP2, rs10830963 in MTNR1B, rs4607517 in GCK) were genotyped using TaqMan allelic discrimination assays. The genotype and allele distributions of each SNP between the GDM cases and controls and the combined effects of alleles for the risk of developing GDM were analyzed. We found that the rs4402960, rs2383208 and rs391300 were statistically associated with GDM (OR = 1.207, 95%CI = 1.029–1.417, p = 0.021; OR = 1.242, 95%CI = 1.077–1.432, p = 0.003; OR = 1.202, 95%CI = 1.020–1.416, P = 0.028, respectively). In addition, the effect was greater under a recessive model in rs391300 (OR = 1.820, 95%CI = 1.226–2.701, p = 0.003). Meanwhile, the joint effect of these three loci indicated an additive effect of multiple alleles on the risk of developing GDM with an OR of 1.196 per allele (p = 1.08×10⁻⁴). We also found that the risk alleles of rs2383208 (b = −0.085, p = 0.003), rs4402960 (b = −0.057, p = 0.046) and rs10830963 (b = −0.096, p = 0.001) were associated with HOMA-B, while rs7754840 was associated with decrease in insulin AUC during a 100 g OGTT given at the time of GDM diagnosis (b = −0.080, p = 0.007).

Conclusions/Significance

Several risk alleles of type 2 diabetes were associated with GDM in pregnant Chinese women. The effects of these SNPs on GDM might be through the impairment of beta cell function and these risk loci contributed additively to the disease.     

Association of CAPN10 SNPs and haplotypes with polycystic ovary syndrome among South Indian Women

Polycystic Ovary Syndrome (PCOS) is known to be characterized by metabolic disorder in which hyperinsulinemia and peripheral insulin resistance are central features. Given the physiological overlap between PCOS and type-2 diabetes (T2DM), and calpain 10 gene (CAPN10) being a strong candidate for T2DM, a number of studies have analyzed CAPN10 SNPs among PCOS women yielding contradictory results. Our study is first of its kind to investigate the association pattern of CAPN10 polymorphisms (UCSNP-44, 43, 56, 19 and 63) with PCOS among Indian women. 250 PCOS cases and 299 controls from Southern India were recruited for this study. Allele and genotype frequencies of the SNPs were determined and compared between the cases and controls. Results show significant association of UCSNP-44 genotype CC with PCOS (p = 0.007) with highly significant odds ratio when compared to TC (OR = 2.51, p = 0.003, 95% CI = 1.37–4.61) as well as TT (OR = 1.94, p = 0.016, 95% CI = 1.13–3.34). While the haplotype carrying the SNP-44 and SNP-19 variants (21121) exhibited a 2 fold increase in the risk for PCOS (OR = 2.37, p = 0.03), the haplotype containing SNP-56 and SNP-19 variants (11221) seems to have a protective role against PCOS (OR = 0.20, p = 0.004). Our results support the earlier evidence for a possible role of UCSNP-44 of the CAPN10 gene in the manifestation of PCOS.     

Genetic association study identifies HSPB7 as a risk gene for idiopathic dilated cardiomyopathy

Dilated cardiomyopathy (DCM) is a structural heart disease with strong genetic background. Monogenic forms of DCM are observed in families with mutations located mostly in genes encoding structural and sarcomeric proteins. However, strong evidence suggests that genetic factors also affect the susceptibility to idiopathic DCM. To identify risk alleles for non-familial forms of DCM, we carried out a case-control association study, genotyping 664 DCM cases and 1,874 population-based healthy controls from Germany using a 50K human cardiovascular disease bead chip covering more than 2,000 genes pre-selected for cardiovascular relevance. After quality control, 30,920 single nucleotide polymorphisms (SNP) were tested for association with the disease by logistic regression adjusted for gender, and results were genomic-control corrected. The analysis revealed a significant association between a SNP in HSPB7 gene (rs1739843, minor allele frequency 39%) and idiopathic DCM (p = 1.06×10⁻⁶, OR = 0.67 [95% CI 0.57–0.79] for the minor allele T). Three more SNPs showed p < 2.21×10⁻⁵. De novo genotyping of these four SNPs was done in three independent case-control studies of idiopathic DCM. Association between SNP rs1739843 and DCM was significant in all replication samples: Germany (n = 564, n = 981 controls, p = 2.07×10⁻³, OR = 0.79 [95% CI 0.67–0.92]), France 1 (n = 433 cases, n = 395 controls, p = 3.73×10⁻³, OR = 0.74 [95% CI 0.60–0.91]), and France 2 (n = 249 cases, n = 380 controls, p = 2.26×10⁻⁴, OR = 0.63 [95% CI 0.50–0.81]). The combined analysis of all four studies including a total of n = 1,910 cases and n = 3,630 controls showed highly significant evidence for association between rs1739843 and idiopathic DCM (p = 5.28×10⁻¹³, OR = 0.72 [95% CI 0.65–0.78]). None of the other three SNPs showed significant results in the replication stage.This finding of the HSPB7 gene from a genetic search for idiopathic DCM using a large SNP panel underscores the influence of common polymorphisms on DCM susceptibility.     

FTO gene associated fatness in relation to body fat distribution and metabolic traits throughout a broad range of fatness

Background

A common single nucleotide polymorphism (SNP) of FTO (rs9939609, T/A) is associated with total body fatness. We investigated the association of this SNP with abdominal and peripheral fatness and obesity-related metabolic traits in middle-aged men through a broad range of fatness present already in adolescence.

Methodology/Principal Findings

Obese young Danish men (n = 753, BMI≥31.0 kg/m²) and a randomly selected group (n = 879) from the same population were examined in three surveys (mean age 35, 46 and 49 years, respectively). The traits included anthropometrics, body composition, oral glucose tolerance test, blood lipids, blood pressure, fibrinogen and aspartate aminotransferase. Logistic regression analysis was used to assess the age-adjusted association between the phenotypes and the odds ratios for the FTO rs9939609 (TT and TA genotype versus the AA genotype), for anthropometrics and body composition estimated per unit z-score. BMI was strongly associated with the AA genotype in all three surveys: OR = 1.17, p = 1.1*10⁻⁶, OR = 1.20, p = 1.7*10⁻⁷, OR = 1.17, p = 3.4*10⁻³, respectively. Fat body mass index was also associated with the AA genotype (OR = 1.21, p = 4.6*10⁻⁷ and OR = 1.21, p = 1.0*10⁻³). Increased abdominal fatness was associated with the AA genotype when measured as waist circumference (OR = 1.21, p = 2.2*10⁻⁶ and OR = 1.19, p = 5.9*10⁻³), sagittal abdominal diameter (OR = 1.17, p = 1.3*10⁻⁴ and OR = 1.18, p = 0.011) and intra-abdominal adipose tissue (OR = 1.21, p = 0.005). Increased peripheral fatness measured as hip circumference (OR = 1.19, p = 1.3*10⁻⁵ and OR = 1.18, p = 0.004) and lower body fat mass (OR = 1.26, p = 0.002) was associated with the AA genotype. The AA genotype was significantly associated with decreased Stumvoll insulin sensitivity index (OR = 0.93, p = 0.02) and with decreased non-fasting plasma HDL-cholesterol (OR = 0.57, p = 0.037), but not with any other of the metabolic traits. However, all significant results for both body fat distribution and metabolic traits were explained by a mediating effect of total fat mass.

Conclusion

The association of the examined FTO SNP to general fatness throughout the range of fatness was confirmed, and this association explains the relation between the SNP and body fat distribution and decreased insulin sensitivity and HDL-cholesterol. The SNP was not significantly associated with other metabolic traits suggesting that they are not derived from the general accumulation of body fat.     

Genome-wide association study of blood pressure extremes identifies variant near UMOD associated with hypertension

Hypertension is a heritable and major contributor to the global burden of disease. The sum of rare and common genetic variants robustly identified so far explain only 1%–2% of the population variation in BP and hypertension. This suggests the existence of more undiscovered common variants. We conducted a genome-wide association study in 1,621 hypertensive cases and 1,699 controls and follow-up validation analyses in 19,845 cases and 16,541 controls using an extreme case-control design. We identified a locus on chromosome 16 in the 5′ region of Uromodulin (UMOD; rs13333226, combined P value of 3.6×10⁻¹¹). The minor G allele is associated with a lower risk of hypertension (OR [95%CI]: 0.87 [0.84–0.91]), reduced urinary uromodulin excretion, better renal function; and each copy of the G allele is associated with a 7.7% reduction in risk of CVD events after adjusting for age, sex, BMI, and smoking status (H.R. = 0.923, 95% CI 0.860–0.991; p = 0.027). In a subset of 13,446 individuals with estimated glomerular filtration rate (eGFR) measurements, we show that rs13333226 is independently associated with hypertension (unadjusted for eGFR: 0.89 [0.83–0.96], p = 0.004; after eGFR adjustment: 0.89 [0.83–0.96], p = 0.003). In clinical functional studies, we also consistently show the minor G allele is associated with lower urinary uromodulin excretion. The exclusive expression of uromodulin in the thick portion of the ascending limb of Henle suggests a putative role of this variant in hypertension through an effect on sodium homeostasis. The newly discovered UMOD locus for hypertension has the potential to give new insights into the role of uromodulin in BP regulation and to identify novel drugable targets for reducing cardiovascular risk.     

Circadian gene variants and susceptibility to type 2 diabetes: a pilot study

Background

Disruption of endogenous circadian rhythms has been shown to increase the risk of developing type 2 diabetes, suggesting that circadian genes might play a role in determining disease susceptibility. We present the results of a pilot study investigating the association between type 2 diabetes and selected single nucleotide polymorphisms (SNPs) in/near nine circadian genes. The variants were chosen based on their previously reported association with prostate cancer, a disease that has been suggested to have a genetic link with type 2 diabetes through a number of shared inherited risk determinants.

Methodology/Principal Findings

The pilot study was performed using two genetically homogeneous Punjabi cohorts, one resident in the United Kingdom and one indigenous to Pakistan. Subjects with (N = 1732) and without (N = 1780) type 2 diabetes were genotyped for thirteen circadian variants using a competitive allele-specific polymerase chain reaction method. Associations between the SNPs and type 2 diabetes were investigated using logistic regression. The results were also combined with in silico data from other South Asian datasets (SAT2D consortium) and white European cohorts (DIAGRAM+) using meta-analysis. The rs7602358G allele near PER2 was negatively associated with type 2 diabetes in our Punjabi cohorts (combined odds ratio [OR] = 0.75 [0.66–0.86], p = 3.18×10⁻⁵), while the BMAL1 rs11022775T allele was associated with an increased risk of the disease (combined OR = 1.22 [1.07–1.39], p = 0.003). Neither of these associations was replicated in the SAT2D or DIAGRAM+ datasets, however. Meta-analysis of all the cohorts identified disease associations with two variants, rs2292912 in CRY2 and rs12315175 near CRY1, although statistical significance was nominal (combined OR = 1.05 [1.01–1.08], p = 0.008 and OR = 0.95 [0.91–0.99], p = 0.015 respectively).

Conclusions/significance

None of the selected circadian gene variants was associated with type 2 diabetes with study-wide significance after meta-analysis. The nominal association observed with the CRY2 SNP, however, complements previous findings and confirms a role for this locus in disease susceptibility.     

Risk alleles for systemic lupus erythematosus in a large case-control collection and associations with clinical subphenotypes

Systemic lupus erythematosus (SLE) is a genetically complex disease with heterogeneous clinical manifestations. Recent studies have greatly expanded the number of established SLE risk alleles, but the distribution of multiple risk alleles in cases versus controls and their relationship to subphenotypes have not been studied. We studied 22 SLE susceptibility polymorphisms with previous genome-wide evidence of association (p<5×10⁻⁸) in 1919 SLE cases from 9 independent Caucasian SLE case series and 4813 independent controls. The mean number of risk alleles in cases was 15.1 (SD 3.1) while the mean in controls was 13.1 (SD 2.8), with trend p = 4×10⁻¹²⁸. We defined a genetic risk score (GRS) for SLE as the number of risk alleles with each weighted by the SLE risk odds ratio (OR). The OR for high-low GRS tertiles, adjusted for intra-European ancestry, sex, and parent study, was 4.4 (95% CI 3.8–5.1). We studied associations of individual SNPs and the GRS with clinical manifestations for the cases: age at diagnosis, the 11 American College of Rheumatology classification criteria, and double-stranded DNA antibody (anti-dsDNA) production. Six subphenotypes were significantly associated with the GRS, most notably anti-dsDNA (OR_high-low = 2.36, p = 9e−9), the immunologic criterion (OR_high-low = 2.23, p = 3e−7), and age at diagnosis (OR_high-low = 1.45, p = 0.0060). Finally, we developed a subphenotype-specific GRS (sub-GRS) for each phenotype with more power to detect cumulative genetic associations. The sub-GRS was more strongly associated than any single SNP effect for 5 subphenotypes (the above plus hematologic disorder and oral ulcers), while single loci are more significantly associated with renal disease (HLA-DRB1, OR = 1.37, 95% CI 1.14–1.64) and arthritis (ITGAM, OR = 0.72, 95% CI 0.59–0.88). We did not observe significant associations for other subphenotypes, for individual loci or the sub-GRS. Thus our analysis categorizes SLE subphenotypes into three groups: those having cumulative, single, and no known genetic association with respect to the currently established SLE risk loci.     

An exploratory study on CLU, CR1 and PICALM and Parkinson disease

Background

Recent GWAS and subsequent confirmation studies reported several single-nucleotide polymorphisms (SNPs) at the CLU, CR1 and PICALM loci in association with late-onset Alzheimer''s disease (AD). Parkinson disease (PD) shares several clinical and pathologic characteristics with AD; we therefore explored whether these SNPs were also associated with PD risk.

Methodology/Principal Findings

791 non-Hispanic Whites cases and 1,580 matched controls were included in the study. Odds ratios (OR) and 95% confidence intervals (CI) were obtained from logistic regression models. rs11136000 at the CLU locus was associated with PD risk under the recessive model (comparing TT versus CC+CT: OR = 0.71, 95% CI: 0.55-0.92, p = 0.008) after adjusting for year of birth, gender, smoking, and caffeine intake. Further adjustment for family history of PD and ApoE ε4 status did not change the result. In addition, we did not find evidence for effect modification by ApoE or known PD risk factors. The association, however, appeared to be stronger for PD with dementia (OR = 0.49, 95% CI: 0.27-0.91) than for PD without dementia (OR = 0.81, 95% CI: 0.61-1.06). The two other SNPs, rs6656401 from CR1, and rs3851179 from PICALM region were not associated with PD (p>0.05).

Conclusion

Our exploratory analysis suggests an association of CLU with PD. This exploratory finding and the role of dementia in explaining this finding needs further investigation.     

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.     

JC Polyomavirus Infection Is Strongly Controlled by Human Leucocyte Antigen Class II Variants

JC polyomavirus (JCV) carriers with a compromised immune system, such as in HIV, or subjects on immune-modulating therapies, such as anti VLA-4 therapy may develop progressive multifocal leukoencephalopathy (PML) which is a lytic infection of oligodendrocytes in the brain. Serum antibodies to JCV mark infection occur only in 50–60% of infected individuals, and high JCV-antibody titers seem to increase the risk of developing PML. We here investigated the role of human leukocyte antigen (HLA), instrumental in immune defense in JCV antibody response. Anti-JCV antibody status, as a surrogate for JCV infection, were compared to HLA class I and II alleles in 1621 Scandinavian persons with MS and 1064 population-based Swedish controls and associations were replicated in 718 German persons with MS. HLA-alleles were determined by SNP imputation, sequence specific (SSP) kits and a reverse PCR sequence-specific oligonucleotide (PCR-SSO) method. An initial GWAS screen displayed a strong HLA class II region signal. The HLA-DRB1*15 haplotype was strongly negatively associated to JCV sero-status in Scandinavian MS cases (OR = 0.42, p = 7×10⁻¹⁵) and controls (OR = 0.53, p = 2×10⁻⁵). In contrast, the DQB1*06:03 haplotype was positively associated with JCV sero-status, in Scandinavian MS cases (OR = 1.63, p = 0.006), and controls (OR = 2.69, p = 1×10⁻⁵). The German dataset confirmed these findings (OR = 0.54, p = 1×10⁻⁴ and OR = 1.58, p = 0.03 respectively for these haplotypes). HLA class II restricted immune responses, and hence CD4+ T cell immunity is pivotal for JCV infection control. Alleles within the HLA-DR1*15 haplotype are associated with a protective effect on JCV infection. Alleles within the DQB1*06:03 haplotype show an opposite association. These associations between JC virus antibody response and human leucocyte antigens supports the notion that CD4+ T cells are crucial in the immune defence to JCV and lays the ground for risk stratification for PML and development of therapy and prevention.     

Association between anaemia, iron deficiency anaemia, neglected parasitic infections and socioeconomic factors in rural children of West Malaysia

Background

Given that micronutrient deficiency, neglected intestinal parasitic infections (IPIs) and poor socioeconomic status are closely linked, we conducted a cross-sectional study to assess the relationship between IPIs and nutritional status of children living in remote and rural areas in West Malaysia.

Methods/Findings

A total of 550 children participated, comprising 520 (94.5%) school children aged 7 to 12 years old, 30 (5.5%) young children aged 1 to 6 years old, 254 (46.2%) boys and 296 (53.8%) girls. Of the 550 children, 26.2% were anaemic, 54.9% iron deficient and 16.9% had iron deficiency anaemia (IDA). The overall prevalence of helminths was 76.5% comprising Trichuris trichiura (71.5%), Ascaris lumbricoides (41.6%) and hookworm infection (13.5%). It was observed that iron deficiency was significantly higher in girls (p = 0.032) compared to boys. Univariate analysis demonstrated that low level of mother''s education (OR = 2.52; 95% CI = 1.38–4.60; p = 0.002), non working parents (OR = 2.18; 95% CI = 2.06–2.31; p = 0.013), low household income (OR = 2.02; 95% CI = 1.14–3.59; p = 0.015), T. trichiura (OR = 2.15; 95% CI = 1.21–3.81; p = 0.008) and A. lumbricoides infections (OR = 1.63; 95% CI = 1.04–2.55; p = 0.032) were significantly associated with the high prevalence of IDA. Multivariate analysis confirmed that low level of mother''s education (OR = 1.48; 95 CI% = 1.33–2.58; p<0.001) was a significant predictor for IDA in these children.

Conclusion

It is crucial that a comprehensive primary health care programme for these communities that includes periodic de-worming, nutrition supplement, improved household economy, education, sanitation status and personal hygiene are taken into consideration to improve the nutritional status of these children.