Pharmacogenetic associations of MMP9 and MMP12 variants with cardiovascular disease in patients with hypertension

Objectives

MMP-9 and -12 function in tissue remodeling and may play roles in cardiovascular disease (CVD). We assessed associations of four MMP polymorphisms and three antihypertensive drugs with cardiovascular outcomes.

Methods

Hypertensives (n = 42,418) from a double-blind, randomized, clinical trial were randomized to chlorthalidone, amlodipine, lisinopril, or doxazosin treatment (mean follow up, 4.9 years). The primary outcome was coronary heart disease (CHD). Secondary outcomes included combined CHD, all CVD outcomes combined, stroke, heart failure (HF), and mortality. Genotype-treatment interactions were tested.

Results

There were 38,698 participants genotyped for at least one of the polymorphisms included here. For MMP9 R668Q (rs2274756), lower hazard ratios (HRs) were found for AA subjects for most outcomes when treated with chlorthalidone versus amlodipine (eg., CCHD: GG = 1.00, GA = 1.01, AA = 0.64; P = 0.038). For MMP9 R279Q (rs17576), modest pharmacogenetic findings were observed for combined CHD and the composite CVD outcome. For MMP12 N122S (rs652438), lower HRs were observed for CHD in subjects carrying at least one G allele and being treated with chlorthalidone versus lisinopril (CHD: AA = 1.07, AG = 0.80, GG = 0.49; P = 0.005). In the lisinopril-amlodipine comparison, higher HRs were observed for participants having at least one G allele at the MMP12 N122S locus (CHD: AA = 0.94, AG = 1.19, GG = 1.93; P = 0.041). For MMP12 −82A>G (rs2276109), no pharmacogenetic effect was found for the primary outcome, although lower HRs were observed for AA homozygotes in the chlorthalidone-amlodipine comparison for HF (P = 0.015).

Conclusions

We observed interactions between antihypertensive drugs and MMP9 and MMP12 for CHD and composite CVD. The data suggest that these genes may provide useful clinical information with respect to treatment decisions.     

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.     

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 SNP rs6903956 on chromosome 6p24.1 with angiographical characteristics of coronary atherosclerosis in a chinese population

Objective

To explore the association between rs6903956 and severity of coronary artery disease (CAD) in a Chinese population.

Methods

A cohort of 1075 consecutive patients who underwent coronary arteriography for suspected or known coronary atherosclerosis was enrolled in our study. Coronary atherosclerosis severity was defined by Gensini''s Score System and counts of diseased vessels.

Results

Gensini score frequencies and counts of diseased vessels differed among GG, AG, AA genotype groups at the rs6903956 locus (p = 0.025 for Gensini score frequencies vs. p = 0.024 for counts of diseased vessels, respectively). A univariate logistic regression analysis revealed that the genotype distribution of this SNP was associated significantly with angiographical characteristics of coronary atherosclerosis risk (p = 0.030, odds ratio (OR) = 1.444, 95% confidence interval (CI) = 1.036∼2.013 for AG vs. GG; p = 0.021, OR = 5.896, 95% CI = 1.299∼26.750 for AA vs. GG and p = 0.007, OR = 1.564, 95% CI = 1.132∼2.162 for combined (AG+AA) vs. GG). A multivariate logistic regression analysis indicated that the genotype distribution of the rs6903956 polymorphism be associated significantly with the angiographical characteristics of coronary atherosclerosis risk (p = 0.004, OR = 1.578, 95% CI = 1.155∼2.154 for GG vs. AG vs. AA; p = 0.013, OR = 1.541, 95% CI = 1.097∼2.163 for GG vs. GA+ AA). A stratification analysis revealed that male subjects and smoking subjects had a higher frequency of the rs6903956 heterozygous mutant among higher Gensini score subjects than among lower Gensini score subjects (p = 0.023, OR = 1.579, 95% CI = 1.064∼2.344 for male subgroup; p = 0.005, OR = 2.075, 95% CI = 1.249∼3.448 for smoking subgroup).

Conclusions

Allele A is a risk factor for CAD and the G-to-A allele substitution may underlie the association between rs6903956 and CAD.     

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.     

Large-scale pathway-based analysis of bladder cancer genome-wide association data from five studies of European background

Pathway analysis of genome-wide association studies (GWAS) offer a unique opportunity to collectively evaluate genetic variants with effects that are too small to be detected individually. We applied a pathway analysis to a bladder cancer GWAS containing data from 3,532 cases and 5,120 controls of European background (n = 5 studies). Thirteen hundred and ninety-nine pathways were drawn from five publicly available resources (Biocarta, Kegg, NCI-PID, HumanCyc, and Reactome), and we constructed 22 additional candidate pathways previously hypothesized to be related to bladder cancer. In total, 1421 pathways, 5647 genes and ∼90,000 SNPs were included in our study. Logistic regression model adjusting for age, sex, study, DNA source, and smoking status was used to assess the marginal trend effect of SNPs on bladder cancer risk. Two complementary pathway-based methods (gene-set enrichment analysis [GSEA], and adapted rank-truncated product [ARTP]) were used to assess the enrichment of association signals within each pathway. Eighteen pathways were detected by either GSEA or ARTP at P≤0.01. To minimize false positives, we used the I² statistic to identify SNPs displaying heterogeneous effects across the five studies. After removing these SNPs, seven pathways (‘Aromatic amine metabolism’ [P_GSEA = 0.0100, P_ARTP = 0.0020], ‘NAD biosynthesis’ [P_GSEA = 0.0018, P_ARTP = 0.0086], ‘NAD salvage’ [P_ARTP = 0.0068], ‘Clathrin derived vesicle budding’ [P_ARTP = 0.0018], ‘Lysosome vesicle biogenesis’ [P_GSEA = 0.0023, P_ARTP<0.00012], ’Retrograde neurotrophin signaling’ [P_GSEA = 0.00840], and ‘Mitotic metaphase/anaphase transition’ [P_GSEA = 0.0040]) remained. These pathways seem to belong to three fundamental cellular processes (metabolic detoxification, mitosis, and clathrin-mediated vesicles). Identification of the aromatic amine metabolism pathway provides support for the ability of this approach to identify pathways with established relevance to bladder carcinogenesis.     

Pharmacogenetic association of NOS3 variants with cardiovascular disease in patients with hypertension: the GenHAT study

Nitric oxide synthase 3 (NOS3) catalyzes production of NO in the endothelium and may play a role in cardiovascular disease (CVD). We assessed the pharmacogenetic associations of three NOS3 polymorphisms and three antihypertensive drugs with CVD outcomes. Hypertensive subjects (n = 30,280) from a multi-center, double-blind clinical trial were randomized to chlorthalidone, amlodipine, or lisinopril treatment (mean follow up, 4.9 years). Outcomes included coronary heart disease (CHD: fatal CHD and nonfatal myocardial infarction); stroke; heart failure (fatal, requiring hospitalization, or outpatient treatment); all-cause mortality; and end-stage renal disease (ESRD). Main effects of NOS3 variants on outcome and genotype-treatment interactions were tested. For NOS3 −690 C>T (rs3918226), a higher hazard ratio (HR) was found in minor allele carriers for CHD (CC = 1.00, CT+TT = 1.12 (95% confidence interval (CI) = 1.00–1.26), P = 0.048). For NOS3 −922 A>G (rs1800779), a higher HR was found in minor allele carriers for heart failure (AA = 1.00, AG+GG = 1.10 (CI = 1.00–1.21), P = 0.046). Significant pharmacogenetic findings were observed for stroke and all-cause mortality. For −690 C>T, a lower HR was observed for stroke in minor allele carriers when treated with amlodipine versus lisinopril (CC = 0.85 (CI = 0.73–0.99), CT+TT = 0.49 (CI = 0.31–0.80), P = 0.04). For glu298asp G>T (rs1799983), a lower HR was observed for all-cause mortality in minor allele carriers when treated with amlodipine versus lisinopril (GG = 1.01 (CI = 0.91–1.13), GT+TT = 0.85 (CI = 0.75–0.97), P = 0.04). We observed significant associations with NOS3 variants and CHD and heart failure and significant pharmacogenetic effects for stroke and all cause mortality. This suggests that NOS3 variants may potentially provide useful clinical information with respect to treatment decisions in the future.     

Meta-analysis for genome-wide association study identifies multiple variants at the BIN1 locus associated with late-onset Alzheimer's disease

Recent GWAS studies focused on uncovering novel genetic loci related to AD have revealed associations with variants near CLU, CR1, PICALM and BIN1. In this study, we conducted a genome-wide association study in an independent set of 1034 cases and 1186 controls using the Illumina genotyping platforms. By coupling our data with available GWAS datasets from the ADNI and GenADA, we replicated the original associations in both PICALM (rs3851179) and CR1 (rs3818361). The PICALM variant seems to be non-significant after we adjusted for APOE e4 status. We further tested our top markers in 751 independent cases and 751 matched controls. Besides the markers close to the APOE locus, a marker (rs12989701) upstream of BIN1 locus was replicated and the combined analysis reached genome-wide significance level (p = 5E-08). We combined our data with the published Harold et al. study and meta-analysis with all available 6521 cases and 10360 controls at the BIN1 locus revealed two significant variants (rs12989701, p = 1.32E-10 and rs744373, p = 3.16E-10) in limited linkage disequilibrium (r² = 0.05) with each other. The independent contribution of both SNPs was supported by haplotype conditional analysis. We also conducted multivariate analysis in canonical pathways and identified a consistent signal in the downstream pathways targeted by Gleevec (P = 0.004 in Pfizer; P = 0.028 in ADNI and P = 0.04 in GenADA). We further tested variants in CLU, PICALM, BIN1 and CR1 for association with disease progression in 597 AD patients where longitudinal cognitive measures are sufficient. Both the PICALM and CLU variants showed nominal significant association with cognitive decline as measured by change in Clinical Dementia Rating-sum of boxes (CDR-SB) score from the baseline but did not pass multiple-test correction. Future experiments will help us better understand potential roles of these genetic loci in AD pathology.     

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.     

Multiple analytical approaches reveal distinct gene-environment interactions in smokers and non smokers in lung cancer

Complex disease such as cancer results from interactions of multiple genetic and environmental factors. Studying these factors singularly cannot explain the underlying pathogenetic mechanism of the disease. Multi-analytical approach, including logistic regression (LR), classification and regression tree (CART) and multifactor dimensionality reduction (MDR), was applied in 188 lung cancer cases and 290 controls to explore high order interactions among xenobiotic metabolizing genes and environmental risk factors. Smoking was identified as the predominant risk factor by all three analytical approaches. Individually, CYP1A1*2A polymorphism was significantly associated with increased lung cancer risk (OR = 1.69;95%CI = 1.11–2.59,p = 0.01), whereas EPHX1 Tyr113His and SULT1A1 Arg213His conferred reduced risk (OR = 0.40;95%CI = 0.25–0.65,p<0.001 and OR = 0.51;95%CI = 0.33–0.78,p = 0.002 respectively). In smokers, EPHX1 Tyr113His and SULT1A1 Arg213His polymorphisms reduced the risk of lung cancer, whereas CYP1A1*2A, CYP1A1*2C and GSTP1 Ile105Val imparted increased risk in non-smokers only. While exploring non-linear interactions through CART analysis, smokers carrying the combination of EPHX1 113TC (Tyr/His), SULT1A1 213GG (Arg/Arg) or AA (His/His) and GSTM1 null genotypes showed the highest risk for lung cancer (OR = 3.73;95%CI = 1.33–10.55,p = 0.006), whereas combined effect of CYP1A1*2A 6235CC or TC, SULT1A1 213GG (Arg/Arg) and betel quid chewing showed maximum risk in non-smokers (OR = 2.93;95%CI = 1.15–7.51,p = 0.01). MDR analysis identified two distinct predictor models for the risk of lung cancer in smokers (tobacco chewing, EPHX1 Tyr113His, and SULT1A1 Arg213His) and non-smokers (CYP1A1*2A, GSTP1 Ile105Val and SULT1A1 Arg213His) with testing balance accuracy (TBA) of 0.6436 and 0.6677 respectively. Interaction entropy interpretations of MDR results showed non-additive interactions of tobacco chewing with SULT1A1 Arg213His and EPHX1 Tyr113His in smokers and SULT1A1 Arg213His with GSTP1 Ile105Val and CYP1A1*2C in nonsmokers. These results identified distinct gene-gene and gene environment interactions in smokers and non-smokers, which confirms the importance of multifactorial interaction in risk assessment of lung cancer.     

Age and sex composition of seals killed by polar bears in the eastern Beaufort Sea

Background

Polar bears (Ursus maritimus) of the Beaufort Sea enter hyperphagia in spring and gain fat reserves to survive periods of low prey availability. We collected information on seals killed by polar bears (n = 650) and hunting attempts on ringed seal (Pusa hispida) lairs (n = 1396) observed from a helicopter during polar bear mark-recapture studies in the eastern Beaufort Sea in spring in 1985–2011. We investigated how temporal shifts in ringed seal reproduction affect kill composition and the intraspecific vulnerabilities of ringed seals to polar bear predation.

Principal Findings

Polar bears primarily preyed on ringed seals (90.2%) while bearded seals (Erignathus barbatus) only comprised 9.8% of the kills, but 33% of the biomass. Adults comprised 43.6% (150/344) of the ringed seals killed, while their pups comprised 38.4% (132/344). Juvenile ringed seals were killed at the lowest proportion, comprising 18.0% (62/344) of the ringed seal kills. The proportion of ringed seal pups was highest between 2007–2011, in association with high ringed seal productivity. Half of the adult ringed seal kills were ≥21 years (60/121), and kill rates of adults increased following the peak of parturition. Determination of sex from DNA revealed that polar bears killed adult male and adult female ringed seals equally (0.50, n = 78). The number of hunting attempts at ringed seal subnivean lair sites was positively correlated with the number of pup kills (r² = 0.30, P = 0.04), but was not correlated with the number of adult kills (P = 0.37).

Conclusions/Significance

Results are consistent with decadal trends in ringed seal productivity, with low numbers of pups killed by polar bears in spring in years of low pup productivity, and conversely when pup productivity was high. Vulnerability of adult ringed seals to predation increased in relation to reproductive activities and age, but not gender.     

Validation of a new test for Schistosoma haematobium based on detection of Dra1 DNA fragments in urine: evaluation through latent class analysis

Background

Diagnosis of urogenital schistosomiasis in chronically infected adults is challenging but important, especially because long term infection of the bladder and urinary tract can have dire consequences. We evaluated three tests for viable infection: detection of parasite specific DNA Dra1 fragments, haematuria and presence of parasite eggs for sensitivity (Se) and specificity (Sp).

Methods

Over 400 urine specimens collected from adult volunteers in an endemic area in Western Nigeria were assessed for haematuria then filtered in the field, the filter papers dried and later examined for eggs and DNA. The results were stratified according to sex and age and subjected to Latent Class analysis.

Conclusions

Presence of Dra1 in males (Se = 100%; Sp = 100%) exceeded haematuria (Se = 87.6%: Sp = 34.7%) and detection of eggs (Se = 70.1%; Sp = 100%). In females presence of Dra1 was Se = 100%: Sp = 100%, exceeding haematuria (Se = 86.7%: Sp = 77.0%) and eggs (Se = 70.1%; Sp = 100%). Dra1 became undetectable 2 weeks after praziquantel treatment. We conclude detection of Dra1 fragment is a definitive test for the presence of Schistosoma haematobium infection.     

Evidence for more than one Parkinson's disease-associated variant within the HLA region

Parkinson''s disease (PD) was recently found to be associated with HLA in a genome-wide association study (GWAS). Follow-up GWAS''s replicated the PD-HLA association but their top hits differ. Do the different hits tag the same locus or is there more than one PD-associated variant within HLA? We show that the top GWAS hits are not correlated with each other (0.00≤r²≤0.15). Using our GWAS (2000 cases, 1986 controls) we conducted step-wise conditional analysis on 107 SNPs with P<10⁻³ for PD-association; 103 dropped-out, four remained significant. Each SNP, when conditioned on the other three, yielded P_SNP1 = 5×10⁻⁴, P_SNP2 = 5×10⁻⁴, P_SNP3 = 4×10⁻³ and P_SNP4 = 0.025. The four SNPs were not correlated (0.01≤r²≤0.20). Haplotype analysis (excluding rare SNP2) revealed increasing PD risk with increasing risk alleles from OR = 1.27, P = 5×10⁻³ for one risk allele to OR = 1.65, P = 4×10⁻⁸ for three. Using additional 843 cases and 856 controls we replicated the independent effects of SNP1 (P_{conditioned-on-SNP4} = 0.04) and SNP4 (P_{conditioned-on-SNP1} = 0.04); SNP2 and SNP3 could not be replicated. In pooled GWAS and replication, SNP1 had OR_{conditioned-on-SNP4} = 1.23, P_{conditioned-on-SNP4} = 6×10⁻⁷; SNP4 had OR_{conditioned-on-SNP1} = 1.18, P_{conditioned-on-SNP1} = 3×10⁻³; and the haplotype with both risk alleles had OR = 1.48, P = 2×10⁻¹². Genotypic OR increased with the number of risk alleles an individual possessed up to OR = 1.94, P = 2×10⁻¹¹ for individuals who were homozygous for the risk allele at both SNP1 and SNP4. SNP1 is a variant in HLA-DRA and is associated with HLA-DRA, DRB5 and DQA2 gene expression. SNP4 is correlated (r² = 0.95) with variants that are associated with HLA-DQA2 expression, and with the top HLA SNP from the IPDGC GWAS (r² = 0.60). Our findings suggest more than one PD-HLA association; either different alleles of the same gene, or separate loci.     

Prediction of promiscuous p-glycoprotein inhibition using a novel machine learning scheme

Background

P-glycoprotein (P-gp) is an ATP-dependent membrane transporter that plays a pivotal role in eliminating xenobiotics by active extrusion of xenobiotics from the cell. Multidrug resistance (MDR) is highly associated with the over-expression of P-gp by cells, resulting in increased efflux of chemotherapeutical agents and reduction of intracellular drug accumulation. It is of clinical importance to develop a P-gp inhibition predictive model in the process of drug discovery and development.

Methodology/Principal Findings

An in silico model was derived to predict the inhibition of P-gp using the newly invented pharmacophore ensemble/support vector machine (PhE/SVM) scheme based on the data compiled from the literature. The predictions by the PhE/SVM model were found to be in good agreement with the observed values for those structurally diverse molecules in the training set (n = 31, r ² = 0.89, q ² = 0.86, RMSE = 0.40, s = 0.28), the test set (n = 88, r ² = 0.87, RMSE = 0.39, s = 0.25) and the outlier set (n = 11, r ² = 0.96, RMSE = 0.10, s = 0.05). The generated PhE/SVM model also showed high accuracy when subjected to those validation criteria generally adopted to gauge the predictivity of a theoretical model.

Conclusions/Significance

This accurate, fast and robust PhE/SVM model that can take into account the promiscuous nature of P-gp can be applied to predict the P-gp inhibition of structurally diverse compounds that otherwise cannot be done by any other methods in a high-throughput fashion to facilitate drug discovery and development by designing drug candidates with better metabolism profile.     

High-density SNP mapping of the HLA region identifies multiple independent susceptibility loci associated with selective IgA deficiency

Selective IgA deficiency (IgAD; serum IgA<0.07 g/l) is the most common form of human primary immune deficiency, affecting approximately 1∶600 individuals in populations of Northern European ancestry. The polygenic nature of IgAD is underscored by the recent identification of several new risk genes in a genome-wide association study. Among the characterized susceptibility loci, the association with specific HLA haplotypes represents the major genetic risk factor for IgAD. Despite the robust association, the nature and location of the causal variants in the HLA region remains unknown. To better characterize the association signal in this region, we performed a high-density SNP mapping of the HLA locus and imputed the genotypes of common HLA-B, -DRB1, and -DQB1 alleles in a combined sample of 772 IgAD patients and 1,976 matched controls from 3 independent European populations. We confirmed the complex nature of the association with the HLA locus, which is the result of multiple effects spanning the entire HLA region. The primary association signal mapped to the HLA-DQB1*02 allele in the HLA Class II region (combined P = 7.69×10⁻⁵⁷; OR = 2.80) resulting from the combined independent effects of the HLA-B*0801-DRB1*0301-DQB1*02 and -DRB1*0701-DQB1*02 haplotypes, while additional secondary signals were associated with the DRB1*0102 (combined P = 5.86×10⁻¹⁷; OR = 4.28) and the DRB1*1501 (combined P = 2.24×10⁻³⁵; OR = 0.13) alleles. Despite the strong population-specific frequencies of HLA alleles, we found a remarkable conservation of these effects regardless of the ethnic background, which supports the use of large multi-ethnic populations to characterize shared genetic association signals in the HLA region. We also provide evidence for the location of association signals within the specific extended haplotypes, which will guide future sequencing studies aimed at characterizing the precise functional variants contributing to disease pathogenesis.     

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

Background

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

Methodology/Principal Findings

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

Conclusions/Significance

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

Genome-wide association of familial late-onset Alzheimer's disease replicates BIN1 and CLU and nominates CUGBP2 in interaction with APOE

Late-onset Alzheimer''s disease (LOAD) is the most common form of dementia in the elderly. The National Institute of Aging-Late Onset Alzheimer''s Disease Family Study and the National Cell Repository for Alzheimer''s Disease conducted a joint genome-wide association study (GWAS) of multiplex LOAD families (3,839 affected and unaffected individuals from 992 families plus additional unrelated neurologically evaluated normal subjects) using the 610 IlluminaQuad panel. This cohort represents the largest family-based GWAS of LOAD to date, with analyses limited here to the European-American subjects. SNPs near APOE gave highly significant results (e.g., rs2075650, p = 3.2×10⁻⁸¹), but no other genome-wide significant evidence for association was obtained in the full sample. Analyses that stratified on APOE genotypes identified SNPs on chromosome 10p14 in CUGBP2 with genome-wide significant evidence for association within APOE ε4 homozygotes (e.g., rs201119, p = 1.5×10⁻⁸). Association in this gene was replicated in an independent sample consisting of three cohorts. There was evidence of association for recently-reported LOAD risk loci, including BIN1 (rs7561528, p = 0.009 with, and p = 0.03 without, APOE adjustment) and CLU (rs11136000, p = 0.023 with, and p = 0.008 without, APOE adjustment), with weaker support for CR1. However, our results provide strong evidence that association with PICALM (rs3851179, p = 0.69 with, and p = 0.039 without, APOE adjustment) and EXOC3L2 is affected by correlation with APOE, and thus may represent spurious association. Our results indicate that genetic structure coupled with ascertainment bias resulting from the strong APOE association affect genome-wide results and interpretation of some recently reported associations. We show that a locus such as APOE, with large effects and strong association with disease, can lead to samples that require appropriate adjustment for this locus to avoid both false positive and false negative evidence of association. We suggest that similar adjustments may also be needed for many other large multi-site studies.     

Patient's origin and lifestyle associated with CTX-M-producing Escherichia coli: a case-control-control study

Background

Global dissemination of Escherichia coli producing CTX-M extended-spectrum β-lactamases (ESBL) is a public health concern. The aim of the study was to determine factors associated with CTX-M- producing E. coli infections among patients hospitalised in the Assistance Publique-Hôpitaux de Paris, the largest hospital system in France (23 000 beds), through a prospective case-control-control study.

Methods/Principal Findings

From November 2008 to June 2009, 152 inpatients with a clinical sample positive for CTX-M-producing E. coli (cases), 152 inpatients with a clinical sample positive for non ESBL-producing E. coli on the day or within the three days following case detection (controls C1), and 152 inpatients with culture-negative clinical samples since the beginning of hospitalisation and until three days after case detection (controls C2) were included in ten hospitals of the Paris area. Factors studied were related to patient''s origin, lifestyle and medical history as well as care during hospitalisation. Those independently associated with CTX-M-producing E. coli were determined. Three independent factors were common to the two case-control comparisons: birth outside of Europe (cases vs C1: OR₁ = 2.4; 95%CI = [1.3–4.5] and cases vs C2: OR₂ = 3.1; 95%CI = [1.4–7.0]), chronic infections (OR₁ = 2.9; 95%CI = [1.3–6.9] and OR₂ = 8.7; 95%CI = [2.0–39.7]), and antibiotic treatment between hospital admission and inclusion (OR₁ = 2.0; 95%CI = [1.0–3.8] and OR₂ = 3.3; 95%CI = [1.5–7.2]). Cases were also more likely to be (i) functionally dependent before hospitalisation than C2 (OR₂ = 7.0; 95%CI = [2.1–23.5]) and (ii) living in collective housing before hospitalisation than C2 (OR₂ = 15.2; 95%CI = [1.8–130.7]) when CTX-M-producing E. coli was present at admission.

Conclusion

For the first time, patient''s origin and lifestyle were demonstrated to be independently associated with isolation of CTX-M-producing E. coli, in addition to health care-related factors.     

Genetic and non-genetic influences during pregnancy on infant global and site specific DNA methylation: role for folate gene variants and vitamin B12

Inter-individual variation in patterns of DNA methylation at birth can be explained by the influence of environmental, genetic and stochastic factors. This study investigates the genetic and non-genetic determinants of variation in DNA methylation in human infants. Given its central role in provision of methyl groups for DNA methylation, this study focuses on aspects of folate metabolism. Global (LUMA) and gene specific (IGF2, ZNT5, IGFBP3) DNA methylation were quantified in 430 infants by Pyrosequencing®. Seven polymorphisms in 6 genes (MTHFR, MTRR, FOLH1, CβS, RFC1, SHMT) involved in folate absorption and metabolism were analysed in DNA from both infants and mothers. Red blood cell folate and serum vitamin B₁₂ concentrations were measured as indices of vitamin status. Relationships between DNA methylation patterns and several covariates viz. sex, gestation length, maternal and infant red cell folate, maternal and infant serum vitamin B₁₂, maternal age, smoking and genotype were tested. Length of gestation correlated positively with IGF2 methylation (rho = 0.11, p = 0.032) and inversely with ZNT5 methylation (rho = −0.13, p = 0.017). Methylation of the IGFBP3 locus correlated inversely with infant vitamin B₁₂ concentration (rho = −0.16, p = 0.007), whilst global DNA methylation correlated inversely with maternal vitamin B₁₂ concentrations (rho = 0.18, p = 0.044). Analysis of common genetic variants in folate pathway genes highlighted several associations including infant MTRR 66G>A genotype with DNA methylation (χ² = 8.82, p = 0.003) and maternal MTHFR 677C>T genotype with IGF2 methylation (χ² = 2.77, p = 0.006). These data support the hypothesis that both environmental and genetic factors involved in one-carbon metabolism influence DNA methylation in infants. Specifically, the findings highlight the importance of vitamin B₁₂ status, infant MTRR genotype and maternal MTHFR genotype, all of which may influence the supply of methyl groups for DNA methylation. In addition, gestational length appears to be an important determinant of infant DNA methylation patterns.