An investigation of genome-wide studies reported susceptibility loci for ulcerative colitis shows limited replication in north Indians

Genome-Wide Association studies (GWAS) of both Crohn's Disease (CD) and Ulcerative Colitis (UC) have unearthed over 40 risk conferring variants. Recently, a meta-analysis on UC revealed several loci, most of which were either previously associated with UC or CD susceptibility in populations of European origin. In this study, we attempted to replicate these findings in an ethnically distinct north Indian UC cohort. 648 UC cases and 850 controls were genotyped using Infinium Human 660W-quad. Out of 59 meta-analysis index SNPs, six were not in the SNP array used in the study. Of the remaining 53 SNPs, four were found monomorphic. Association (p<0.05) at 25 SNPs was observed, of which 15 were CD specific. Only five SNPs namely rs2395185 (HLA-DRA), rs3024505 (IL10), rs6426833 (RNF186), rs3763313 (BTNL2) and rs2066843 (NOD2) retained significance after Bonferroni correction. These results (i) reveal limited replication of Caucasian based meta-analysis results; (ii) reiterate overlapping molecular mechanism(s) in UC and CD; (iii) indicate differences in genetic architecture between populations; and (iv) suggest that resources such as HapMap need to be extended to cover diverse ethnic populations. They also suggest a systematic GWAS in this terrain may be insightful for identifying population specific IBD risk conferring loci and thus enable cross-ethnicity fine mapping of disease loci.     

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

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

Vitamin D Insufficiency in Arabs and South Asians Positively Associates with Polymorphisms in GC and CYP2R1 Genes

Background

A number of genetic studies have reported an association between vitamin D related genes such as group-specific component gene (GC), Cytochrome P450, family 2, subfamily R, polypeptide 1 (CYP2R1) and 7-dehydrocholesterol reductase/nicotinamide-adenine dinucleotide synthetase 1 (DHCR7/NADSYN1) and serum levels of the active form of Vitamin D, 25 (OH) D among African Americans, Caucasians, and Chinese. Little is known about how genetic variations associate with, or contribute to, 25(OH)D levels in Arabs populations.

Methods

Allele frequencies of 18 SNPs derived from CYP2R1, GC, and DHCR7/NADSYN1 genes in 1549 individuals (Arabs, South Asians, and Southeast Asians living in Kuwait) were determined using real time genotyping assays. Serum levels of 25(OH)D were measured using chemiluminescence immunoassay.

Results

GC gene polymorphisms (rs17467825, rs3755967, rs2282679, rs7041 and rs2298850) were found to be associated with 25(OH)D serum levels in Arabs and South Asians. Two of the CYP2R1 SNPs (rs10500804 and rs12794714) and one of GC SNPs (rs1155563) were found to be significantly associated with 25(OH)D serum levels only in people of Arab origin. Across all three ethnicities none of the SNPs of DHCR7/NADSYN1 were associated with serum 25(OH)D levels and none of the 18 SNPs were significantly associated with serum 25(OH)D levels in people from South East Asia.

Conclusion

Our data show for the first time significant association between the GC (rs2282679 and rs7041), CYP2R1 (rs10741657) SNPs and 25(OH)D levels. This supports their roles in vitamin D Insufficiency in Arab and South Asian populations respectively. Interestingly, two of the CYP2R1 SNPs (rs10500804 and rs12794714) and one GC SNP (rs1155563) were found to correlate with vitamin D in Arab population exclusively signifying their importance in this population.     

Comprehensive association analysis of nine candidate genes with serum 25-hydroxy vitamin D levels among healthy Caucasian subjects

Vitamin D deficiency is a common public health problem in the US. It is related to the high risk of rickets, osteoporosis and other diseases. Currently, serum 25-hydroxy vitamin D [25(OH)D] concentration is the best indicator of vitamin D status, and determination of its deficiency or sufficiency. This level has high heritability (28–80%). However, genes contributing to the wide variation in serum 25(OH)D are generally unknown. In this study, we screened nine important genes in vitamin D metabolic pathways using 49 single nucleotide polymorphism (SNP) markers in a group of 156 unrelated healthy Caucasian subjects. Significant confounding factors that may affect serum 25(OH)D variations were used as covariates for the association analyses. An association test for quantitative trait was performed to evaluate the association between candidate genes and serum 25(OH)D levels. Permutation was conducted for correcting multiple testing problems. Evidence of association was observed at SNPs in the CYP2R1 (cytochrome P450, family 2, subfamily R, polypeptide 1) and the GC (vitamin D binding protein) gene. Next, we performed a replication study for six promising SNPs in the gene CYP2R1 and GC, using another group of 340 unrelated healthy Caucasian subjects. Association analyses were conducted in the replication cohort (n = 340) and the pooled cohort (n = 496). The CYP2R1 gene and the GC gene remain significant in the pooled cohort. The results suggest that the CYP2R1 and GC genes may contribute to the variation of serum 25(OH)D levels in healthy populations.     

A replication study of GWAS-derived lipid genes in Asian Indians: the chromosomal region 11q23.3 harbors loci contributing to triglycerides

Recent genome-wide association scans (GWAS) and meta-analysis studies on European populations have identified many genes previously implicated in lipid regulation. Validation of these loci on different global populations is important in determining their clinical relevance, particularly for development of novel drug targets for treating and preventing diabetic dyslipidemia and coronary artery disease (CAD). In an attempt to replicate GWAS findings on a non-European sample, we examined the role of six of these loci (CELSR2-PSRC1-SORT1 rs599839; CDKN2A-2B rs1333049; BUD13-ZNF259 rs964184; ZNF259 rs12286037; CETP rs3764261; APOE-C1-C4-C2 rs4420638) in our Asian Indian cohort from the Sikh Diabetes Study (SDS) comprising 3,781 individuals (2,902 from Punjab and 879 from the US). Two of the six SNPs examined showed convincing replication in these populations of Asian Indian origin. Our study confirmed a strong association of CETP rs3764261 with high-density lipoprotein cholesterol (HDL-C) (p?=?2.03×10(-26)). Our results also showed significant associations of two GWAS SNPs (rs964184 and rs12286037) from BUD13-ZNF259 near the APOA5-A4-C3-A1 genes with triglyceride (TG) levels in this Asian Indian cohort (rs964184: p?=?1.74×10(-17); rs12286037: p?=?1.58×10(-2)). We further explored 45 SNPs in a ～195 kb region within the chromosomal region 11q23.3 (encompassing the BUD13-ZNF259, APOA5-A4-C3-A1, and SIK3 genes) in 8,530 Asian Indians from the London Life Sciences Population (LOLIPOP) (UK) and SDS cohorts. Five more SNPs revealed significant associations with TG in both cohorts individually as well as in a joint meta-analysis. However, the strongest signal for TG remained with BUD13-ZNF259 (rs964184: p?=?1.06×10(-39)). Future targeted deep sequencing and functional studies should enhance our understanding of the clinical relevance of these genes in dyslipidemia and hypertriglyceridemia (HTG) and, consequently, diabetes and CAD.     

Genome-Wide Association Study Identifies Novel Pharmacogenomic Loci For Therapeutic Response to Montelukast in Asthma

Background

Genome-wide association study (GWAS) is a powerful tool to identify novel pharmacogenetic single nucleotide polymorphisms (SNPs). Leukotriene receptor antagonists (LTRAs) are a major class of asthma medications, and genetic factors contribute to variable responses to these drugs. We used GWAS to identify novel SNPs associated with the response to the LTRA, montelukast, in asthmatics.

Methods

Using genome-wide genotype and phenotypic data available from American Lung Association - Asthma Clinical Research Center (ALA-ACRC) cohorts, we evaluated 8-week change in FEV₁ related to montelukast administration in a discovery population of 133 asthmatics. The top 200 SNPs from the discovery GWAS were then tested in 184 additional samples from two independent cohorts.

Results

Twenty-eight SNP associations from the discovery GWAS were replicated. Of these, rs6475448 achieved genome-wide significance (combined P = 1.97 x 10^-09), and subjects from all four studies who were homozygous for rs6475448 showed increased ΔFEV₁ from baseline in response to montelukast.

Conclusions

Through GWAS, we identified a novel pharmacogenomic locus related to improved montelukast response in asthmatics.     

A GC polymorphism associated with serum 25-hydroxyvitamin D level is a risk factor for hip fracture in Japanese patients with rheumatoid arthritis: 10-year follow-up of the Institute of Rheumatology,Rheumatoid Arthritis cohort study

Introduction

Vitamin D deficiency has been reported to be common in patients with rheumatoid arthritis (RA) who have a higher prevalence of osteoporosis and hip fracture than healthy individuals. Genetic variants affecting serum 25-hydroxyvitamin D (25(OH)D) concentration, an indicator of vitamin D status, were recently identified by genome-wide association studies of Caucasian populations. The purpose of this study was to validate the association and to test whether the serum 25(OH)D-linked genetic variants were associated with the occurrence of hip fracture in Japanese RA patients.

Methods

DNA samples of 1,957 Japanese RA patients were obtained from the Institute of Rheumatology, Rheumatoid Arthritis (IORRA) cohort DNA collection. First, five single nucleotide polymorphisms (SNPs) that were reported to be associated with serum 25(OH)D concentration by genome-wide association studies were genotyped. The SNPs that showed a significant association with serum 25(OH)D level in the cross-sectional study were used in the longitudinal analysis of hip fracture risk. The genetic risk for hip fracture was determined by a multivariate Cox proportional hazards model in 1,957 patients with a maximum follow-up of 10 years (median, 8 years).

Results

Multivariate linear regression analyses showed that rs2282679 in GC (the gene encoding group-specific component (vitamin D binding protein)) locus was significantly associated with lower serum 25(OH)D concentration (P = 8.1 × 10^-5). A Cox proportional hazards model indicated that rs2282679 in GC was significantly associated with the occurrence of hip fracture in a recessive model (hazard ratio (95% confidence interval) = 2.52 (1.05-6.05), P = 0.039).

Conclusions

A two-staged analysis demonstrated that rs2282679 in GC was associated with serum 25(OH)D concentration and could be a risk factor for hip fracture in Japanese RA patients.     

Common variation in vitamin D pathway genes predicts circulating 25-hydroxyvitamin D Levels among African Americans

Vitamin D is implicated in a wide range of health outcomes, and although environmental predictors of vitamin D levels are known, the genetic drivers of vitamin D status remain to be clarified. African Americans are a group at particularly high risk for vitamin D insufficiency but to date have been virtually absent from studies of genetic predictors of circulating vitamin D levels. Within the Southern Community Cohort Study, we investigated the association between 94 single nucleotide polymorphisms (SNPs) in five vitamin D pathway genes (GC, VDR, CYP2R1, CYP24A1, CYP27B1) and serum 25-hydroxyvitamin D (25(OH)D) levels among 379 African American and 379 Caucasian participants. We found statistically significant associations with three SNPs (rs2298849 and rs2282679 in the GC gene, and rs10877012 in the CYP27B1 gene), although only for African Americans. A genotype score, representing the number of risk alleles across the three SNPs, alone accounted for 4.6% of the variation in serum vitamin D among African Americans. A genotype score of 5 (vs. 1) was also associated with a 7.1 ng/mL reduction in serum 25(OH)D levels and a six-fold risk of vitamin D insufficiency (<20 ng/mL) (odds ratio 6.0, p = 0.01) among African Americans. With African ancestry determined from a panel of 276 ancestry informative SNPs, we found that high risk genotypes did not cluster among those with higher African ancestry. This study is one of the first to investigate common genetic variation in relation to vitamin D levels in African Americans, and the first to evaluate how vitamin D-associated genotypes vary in relation to African ancestry. These results suggest that further evaluation of genetic contributors to vitamin D status among African Americans may help provide insights regarding racial health disparities or enable the identification of subgroups especially in need of vitamin D-related interventions.     

Identification of PLCL1 gene for hip bone size variation in females in a genome-wide association study

Osteoporosis, the most prevalent metabolic bone disease among older people, increases risk for low trauma hip fractures (HF) that are associated with high morbidity and mortality. Hip bone size (BS) has been identified as one of the key measurable risk factors for HF. Although hip BS is highly genetically determined, genetic factors underlying the trait are still poorly defined. Here, we performed the first genome-wide association study (GWAS) of hip BS interrogating approximately 380,000 SNPs on the Affymetrix platform in 1,000 homogeneous unrelated Caucasian subjects, including 501 females and 499 males. We identified a gene, PLCL1 (phospholipase c-like 1), that had four SNPs associated with hip BS at, or approaching, a genome-wide significance level in our female subjects; the most significant SNP, rs7595412, achieved a p value of 3.72x10(-7). The gene's importance to hip BS was replicated using the Illumina genotyping platform in an independent UK cohort containing 1,216 Caucasian females. Two SNPs of the PLCL1 gene, rs892515 and rs9789480, surrounded by the four SNPs identified in our GWAS, achieved p values of 8.62x10(-3) and 2.44x10(-3), respectively, for association with hip BS. Imputation analyses on our GWAS and the UK samples further confirmed the replication signals; eight SNPs of the gene achieved combined imputed p values<10(-5) in the two samples. The PLCL1 gene's relevance to HF was also observed in a Chinese sample containing 403 females, including 266 with HF and 177 control subjects. A SNP of the PLCL1 gene, rs3771362 that is only approximately 0.6 kb apart from the most significant SNP detected in our GWAS (rs7595412), achieved a p value of 7.66x10(-3) (odds ratio = 0.26) for association with HF. Additional biological support for the role of PLCL1 in BS comes from previous demonstrations that the PLCL1 protein inhibits IP3 (inositol 1,4,5-trisphosphate)-mediated calcium signaling, an important pathway regulating mechanical sensing of bone cells. Our findings suggest that PLCL1 is a novel gene associated with variation in hip BS, and provide new insights into the pathogenesis of HF.     

Genetic Polymorphisms in Vitamin D Metabolism and Signaling Genes and Risk of Breast Cancer: A Nested Case-Control Study

Genetic polymorphisms in vitamin D metabolism and signaling genes have been inconsistently associated with risk of breast cancer, though few studies have examined SNPs in vitamin D-related genes other than the vitamin D receptor (VDR) gene and particularly have not examined the association with the retinoid X receptor alpha (RXRA) gene which may be a key vitamin D pathway gene. We conducted a nested case-control study of 734 cases and 1435 individually matched controls from a population-based prospective cohort study, the Northern Sweden Mammary Screening Cohort. Tag and functional SNPs were genotyped for the VDR, cytochrome p450 24A1 (CYP24A1), and RXRA genes. We also genotyped specific SNPs in four other genes related to vitamin D metabolism and signaling (GC/VDBP, CYP2R1, DHCR7, and CYP27B1). SNPs in the CYP2R1, DHCR7, and VDBP gene regions that were associated with circulating 25(OH)D concentration in GWAS were also associated with plasma 25(OH)D in our study (p-trend <0.005). After taking into account the false discovery rate, these SNPs were not significantly associated with breast cancer risk, nor were any of the other SNPs or haplotypes in VDR, RXRA, and CYP24A1. We observed no statistically significant associations between polymorphisms or haplotypes in key vitamin D-related genes and risk of breast cancer. These results, combined with the observation in this cohort and most other prospective studies of no association of circulating 25(OH)D with breast cancer risk, do not support an association between vitamin D and breast cancer risk.     

Large scale replication study of the association between HLA class II/BTNL2 variants and osteoarthritis of the knee in European-descent populations

Osteoarthritis (OA) is the most common form of arthritis and a major cause of disability. This study evaluates the association in Caucasian populations of two single nucleotide polymorphisms (SNPs) mapping to the Human Leukocyte Antigen (HLA) region and deriving from a genome wide association scan (GWAS) of knee OA in Japanese populations. The frequencies for rs10947262 were compared in 36,408 controls and 5,749 knee OA cases from European-descent populations. rs7775228 was tested in 32,823 controls and 1,837 knee OA cases of European descent. The risk (major) allele at rs10947262 in Caucasian samples was not significantly associated with an odds ratio (OR) =?1.07 (95%CI 0.94 -1.21; p?=?0.28). For rs7775228 the meta-analysis resulted in OR?=?0.94 (95%CI 0.81-1.09; p?=?0.42) for the allele associated with risk in the Japanese GWAS. In Japanese individuals these two SNPs are in strong linkage disequilibrium (LD) (r(2)?=?0.86) with the HLA class II haplotype DRB1*1502 DQA1*0103 DQB1*0601 (frequency 8%). In Caucasian and Chinese samples, using imputed data, these SNPs appear not to be in LD with that haplotype (r(2)<0.07). The rs10947262 and rs7775228 variants are not associated with risk of knee OA in European descent populations and they do not appear tag the same HLA class II haplotype as they do in Japanese individuals.     

Population structure of Hispanics in the United States: the multi-ethnic study of atherosclerosis

Using ∼60,000 SNPs selected for minimal linkage disequilibrium, we perform population structure analysis of 1,374 unrelated Hispanic individuals from the Multi-Ethnic Study of Atherosclerosis (MESA), with self-identification corresponding to Central America (n = 93), Cuba (n = 50), the Dominican Republic (n = 203), Mexico (n = 708), Puerto Rico (n = 192), and South America (n = 111). By projection of principal components (PCs) of ancestry to samples from the HapMap phase III and the Human Genome Diversity Panel (HGDP), we show the first two PCs quantify the Caucasian, African, and Native American origins, while the third and fourth PCs bring out an axis that aligns with known South-to-North geographic location of HGDP Native American samples and further separates MESA Mexican versus Central/South American samples along the same axis. Using k-means clustering computed from the first four PCs, we define four subgroups of the MESA Hispanic cohort that show close agreement with self-identification, labeling the clusters as primarily Dominican/Cuban, Mexican, Central/South American, and Puerto Rican. To demonstrate our recommendations for genetic analysis in the MESA Hispanic cohort, we present pooled and stratified association analysis of triglycerides for selected SNPs in the LPL and TRIB1 gene regions, previously reported in GWAS of triglycerides in Caucasians but as yet unconfirmed in Hispanic populations. We report statistically significant evidence for genetic association in both genes, and we further demonstrate the importance of considering population substructure and genetic heterogeneity in genetic association studies performed in the United States Hispanic population.     

A Genome-Wide Association Study of Total Serum and Mite-Specific IgEs in Asthma Patients

Immunoglobulin E (IgE) is one of the central players in asthma and allergic diseases. Although the serum IgE level, a useful endophenotype, is generally increased in patients with asthma, genetic factors influencing IgE regulation in asthma are still not fully understood. To identify the genetic variations associated with total serum and mite-specific IgEs in asthmatics, a genome-wide association study (GWAS) of 657,366 single nucleotide polymorphisms (SNPs) was performed in 877 Korean asthmatics. This study found that several new genes might be associated with total IgE in asthmatics, such as CRIM1 (rs848512, P = 1.18×10⁻⁶; rs711254, P = 6.73×10⁻⁶), ZNF71 (rs10404342, P = 7.60×10⁻⁶), TLN1 (rs4879926, P = 7.74×10⁻⁶), and SYNPO2 (rs1472066, P = 8.36×10⁻⁶; rs1038770, P = 8.66×10⁻⁶). Regarding the association of specific IgE to house dust mites, it was observed that intergenic SNPs nearby to OPRK1 and LOC730217 might be associated with Dermatophagoides pteronyssinus (D.p.) and Dermatophagoides farinae (D.f.) in asthmatics, respectively. In further pathway analysis, the phosphatidylinositol signaling system and adherens junction pathways were estimated to play a role in the regulation of total IgE levels in asthma. Although functional evaluations and replications of these results in other populations are needed, this GWAS of serum IgE in asthmatics could facilitate improved understanding of the role of the newly identified genetic variants in asthma and its related phenotypes.     

Vitamin D binding protein and monocyte response to 25-hydroxyvitamin D and 1,25-dihydroxyvitamin D: analysis by mathematical modeling

Vitamin D binding protein (DBP) plays a key role in the bioavailability of active 1,25-dihydroxyvitamin D (1,25(OH)(2)D) and its precursor 25-hydroxyvitamin D (25OHD), but accurate analysis of DBP-bound and free 25OHD and 1,25(OH)(2)D is difficult. To address this, two new mathematical models were developed to estimate: 1) serum levels of free 25OHD/1,25(OH)(2)D based on DBP concentration and genotype; 2) the impact of DBP on the biological activity of 25OHD/1,25(OH)(2)D in vivo. The initial extracellular steady state (eSS) model predicted that 50 nM 25OHD and 100 pM 1,25(OH)(2)D), <0.1% 25OHD and <1.5% 1,25(OH)(2)D are 'free' in vivo. However, for any given concentration of total 25OHD, levels of free 25OHD are higher for low affinity versus high affinity forms of DBP. The eSS model was then combined with an intracellular (iSS) model that incorporated conversion of 25OHD to 1,25(OH)(2)D via the enzyme CYP27B1, as well as binding of 1,25(OH)(2)D to the vitamin D receptor (VDR). The iSS model was optimized to 25OHD/1,25(OH)(2)D-mediated in vitro dose-responsive induction of the vitamin D target gene cathelicidin (CAMP) in human monocytes. The iSS model was then used to predict vitamin D activity in vivo (100% serum). The predicted induction of CAMP in vivo was minimal at basal settings but increased with enhanced expression of VDR (5-fold) and CYP27B1 (10-fold). Consistent with the eSS model, the iSS model predicted stronger responses to 25OHD for low affinity forms of DBP. Finally, the iSS model was used to compare the efficiency of endogenously synthesized versus exogenously added 1,25(OH)(2)D. Data strongly support the endogenous model as the most viable mode for CAMP induction by vitamin D in vivo. These novel mathematical models underline the importance of DBP as a determinant of vitamin D 'status' in vivo, with future implications for clinical studies of vitamin D status and supplementation.     

Potential novel candidate polymorphisms identified in genome-wide association study for breast cancer susceptibility

Previous genome-wide association studies (GWAS) have shown several risk alleles to be associated with breast cancer. However, the variants identified so far contribute to only a small proportion of disease risk. The objective of our GWAS was to identify additional novel breast cancer susceptibility variants and to replicate these findings in an independent cohort. We performed a two-stage association study in a cohort of 3,064 women from Alberta, Canada. In Stage I, we interrogated 906,600 single nucleotide polymorphisms (SNPs) on Affymetrix SNP 6.0 arrays using 348 breast cancer cases and 348 controls. We used single-locus association tests to determine statistical significance for the observed differences in allele frequencies between cases and controls. In Stage II, we attempted to replicate 35 significant markers identified in Stage I in an independent study of 1,153 cases and 1,215 controls. Genotyping of Stage II samples was done using Sequenom Mass-ARRAY iPlex platform. Six loci from four different gene regions (chromosomes 4, 5, 16 and 19) showed statistically significant differences between cases and controls in both Stage I and Stage II testing, and also in joint analysis. The identified variants were from EDNRA, ROPN1L, C16orf61 and ZNF577 gene regions. The presented joint analyses from the two-stage study design were not significant after genome-wide correction. The SNPs identified in this study may serve as potential candidate loci for breast cancer risk in a further replication study in Stage III from Alberta population or independent validation in Caucasian cohorts elsewhere.     

Genetic Variants in CASP3, BMP5, and IRS2 Genes May Influence Survival in Prostate Cancer Patients Receiving Androgen-Deprivation Therapy

Several genome-wide association studies (GWAS) have been conducted to identify the common single nucleotide polymorphisms (SNPs) that influence the risk of prostate cancer. It was hypothesized that some prostate cancer-associated SNPs might relate to the clinical outcomes in patients treated for prostate cancer using androgen-deprivation therapy (ADT). A cohort of 601 patients who have received ADT for prostate cancer was genotyped for 29 SNPs that have been associated with prostate cancer in Cancer Genetic Markers of Susceptibility GWAS, and within the genes that have been implicated in cancer. Prognostic significance of these SNPs on the disease progression, prostate cancer-specific mortality (PCSM) and all-cause mortality (ACM) after ADT were assessed by Kaplan-Meier analysis and Cox regression model. Three SNPs, namely CASP3 rs4862396, BMP5 rs3734444 and IRS2 rs7986346, were found to be closely associated with the ACM (P≤0.042), and BMP5 rs3734444 and IRS2 rs7986346 were also noted to be significantly related to the PCSM (P≤0.032) after adjusting for the known clinicopathologic predictors. Moreover, patients carrying a greater number of unfavorable genotypes at the loci of interest had a shorter time to ACM and PCSM during ADT (P for trend <0.001). Our results suggest that CASP3 rs4862396, BMP5 rs3734444 and IRS2 rs7986346 may affect the survival in patients after ADT for prostate cancer, and the analysis of these SNPs can help identify patients at higher risk of poor outcome.     

The Challenges of Genome-Wide Interaction Studies: Lessons to Learn from the Analysis of HDL Blood Levels

Genome-wide association studies (GWAS) have revealed 74 single nucleotide polymorphisms (SNPs) associated with high-density lipoprotein cholesterol (HDL) blood levels. This study is, to our knowledge, the first genome-wide interaction study (GWIS) to identify SNP×SNP interactions associated with HDL levels. We performed a GWIS in the Rotterdam Study (RS) cohort I (RS-I) using the GLIDE tool which leverages the massively parallel computing power of Graphics Processing Units (GPUs) to perform linear regression on all genome-wide pairs of SNPs. By performing a meta-analysis together with Rotterdam Study cohorts II and III (RS-II and RS-III), we were able to filter 181 interaction terms with a p-value<1 · 10⁻⁸ that replicated in the two independent cohorts. We were not able to replicate any of these interaction term in the AGES, ARIC, CHS, ERF, FHS and NFBC-66 cohorts (N _total = 30,011) when adjusting for multiple testing. Our GWIS resulted in the consistent finding of a possible interaction between rs774801 in ARMC8 (ENSG00000114098) and rs12442098 in SPATA8 (ENSG00000185594) being associated with HDL levels. However, p-values do not reach the preset Bonferroni correction of the p-values. Our study suggest that even for highly genetically determined traits such as HDL the sample sizes needed to detect SNP×SNP interactions are large and the 2-step filtering approaches do not yield a solution. Here we present our analysis plan and our reservations concerning GWIS.     

Development of GMDR-GPU for Gene-Gene Interaction Analysis and Its Application to WTCCC GWAS Data for Type 2 Diabetes

Although genome-wide association studies (GWAS) have identified a significant number of single-nucleotide polymorphisms (SNPs) associated with many complex human traits, the susceptibility loci identified so far can explain only a small fraction of the genetic risk. Among other possible explanations, the lack of a comprehensive examination of gene–gene interaction (G×G) is often considered a source of the missing heritability. Previously, we reported a model-free Generalized Multifactor Dimensionality Reduction (GMDR) approach for detecting G×G in both dichotomous and quantitative phenotypes. However, the computational burden and less efficient implementation of the original programs make them impossible to use for GWAS. In this study, we developed a graphics processing unit (GPU)-based GMDR program (named GWAS-GPU), which is able not only to analyze GWAS data but also to run much faster than the earlier version of the GMDR program. As a demonstration of the program, we used the GMDR-GPU software to analyze a publicly available GWAS dataset on type 2 diabetes (T2D) from the Wellcome Trust Case Control Consortium. Through an exhaustive search of pair-wise interactions and a selected search of three- to five-way interactions conditioned on significant pair-wise results, we identified 24 core SNPs in six genes (FTO: rs9939973, rs9940128, rs9922047, rs1121980, rs9939609, rs9930506; TSPAN8: rs1495377; TCF7L2: rs4074720, rs7901695, rs4506565, rs4132670, rs10787472, rs11196205, rs10885409, rs11196208; L3MBTL3: rs10485400, rs4897366; CELF4: rs2852373, rs608489; RUNX1: rs445984, rs1040328, rs990074, rs2223046, rs2834970) that appear to be important for T2D. Of these core SNPs, 11 in FTO, TSPAN8, and TCF7L2 have been reported to be associated with T2D, obesity, or both, providing an independent replication of previously reported SNPs. Importantly, we identified three new susceptibility genes; i.e., L3MBTL3, CELF4, and RUNX1, for T2D, a finding that warrants further investigation with independent samples.     

Seed predation,disease and germination on landslides in Neotropical lower montane wet forest

Abstract. Seed mortality (caused by predators and pathogens) and germination were compared between Puerto Rico and Costa Rica on landslides in lower montane wet forest. Seeds of six common species on five Puerto Rican landslides and four common species on two Costa Rican landslides were used with a Cecropia species and a Gonzalagunia species included at both sites. In the Puerto Rican experiments Cecropia schreberiana was the only species to show significant seed predation (which was due to insects), pathogens grew from all species and fewer seeds were lost to predators than pathogens. Also in Puerto Rico mean germination across all species was 57 % before dispersal (filled seeds collected while still on the tree) and 71 % after, with Phytolacca rivinoides seeds germinating most abundantly, followed in descending order by Ocotea leucoxylon, Cecropia spec, Miconia racemosa, Palicourea riparia and Gonzalagunia spicata. In the Costa Rican experiments three species had significant predation: Cecropia polyphlebia and Urera caracasana (both due to insects) and Witheringia coccoloboides (due to mammals); pathogenic disease caused more seed loss than predation, and germination was high (61 % pre-dispersal, 69 % post-dispersal). Similarities between these island and mainland sites included (1) percentage of seeds lost to predation and percentage lost to pathogens (all in the 5–15 % range), (2) generalist pathogens which claimed more seeds than predators and (3) majority germination with a general increase after dispersal. Finally sites were dissimilar only in the number of species with significant predation loss and whether it was by insects or mammals, casting doubt on the traditional island/mainland dichotomy.