New Components of Drosophila Leg Development Identified through Genome Wide Association Studies

The adult Drosophila melanogaster body develops from imaginal discs, groups of cells set-aside during embryogenesis and expanded in number during larval stages. Specification and development of Drosophila imaginal discs have been studied for many years as models of morphogenesis. These studies are often based on mutations with large developmental effects, mutations that are often lethal in embryos when homozygous. Such forward genetic screens can be limited by factors such as early lethality and genetic redundancy. To identify additional genes and genetic pathways involved in leg imaginal disc development, we employed a Genome Wide Association Study utilizing the natural genetic variation in leg proportionality found in the Drosophila Genetic Reference Panel fly lines. In addition to identifying genes already known to be involved in leg development, we identified several genes involved in pathways that had not previously been linked with leg development. Several of the genes appear to be involved in signaling activities, while others have no known roles at this time. Many of these uncharacterized genes are conserved in mammals, so we can now begin to place these genes into developmental contexts. Interestingly, we identified five genes which, when their function is reduced by RNAi, cause an antenna-to-leg transformation. Our results demonstrate the utility of this approach, integrating the tools of quantitative and molecular genetics to study developmental processes, and provide new insights into the pathways and networks involved in Drosophila leg development.     

The Effect of Single Nucleotide Polymorphisms from Genome Wide Association Studies in Multiple Sclerosis on Gene Expression

Background

Multiple sclerosis (MS) is a complex neurological disorder. Its aetiology involves both environmental and genetic factors. Recent genome-wide association studies have identified a number of single nucleotide polymorphisms (SNPs) associated with susceptibility to (MS). We investigated whether these genetic variations were associated with alteration in gene expression.

Methods/Principal Findings

We used a database of mRNA expression and genetic variation derived from immortalised peripheral lymphocytes to investigate polymorphisms associated with MS for correlation with gene expression. Several SNPs were found to be associated with changes in expression: in particular two with HLA-DQA1, HLA-DQA2, HLA-DQB1, HLA-DRB1, HLA-DRB4 and HLA-DRB5, one with ZFP57, one with CD58, two with IL7 and FAM164A, and one with FAM119B, TSFM and KUB3. We found minimal cross-over with a recent whole genome expression study in MS patients.

Discussion

We have shown that many susceptibility loci in MS are associated with changes in gene expression using an unbiased expression database. Several of these findings suggest novel gene candidates underlying the effects of MS-associated genetic variation.     

Candidate Luminal B Breast Cancer Genes Identified by Genome,Gene Expression and DNA Methylation Profiling

Breast cancers (BCs) of the luminal B subtype are estrogen receptor-positive (ER+), highly proliferative, resistant to standard therapies and have a poor prognosis. To better understand this subtype we compared DNA copy number aberrations (CNAs), DNA promoter methylation, gene expression profiles, and somatic mutations in nine selected genes, in 32 luminal B tumors with those observed in 156 BCs of the other molecular subtypes. Frequent CNAs included 8p11-p12 and 11q13.1-q13.2 amplifications, 7q11.22-q34, 8q21.12-q24.23, 12p12.3-p13.1, 12q13.11-q24.11, 14q21.1-q23.1, 17q11.1-q25.1, 20q11.23-q13.33 gains and 6q14.1-q24.2, 9p21.3-p24,3, 9q21.2, 18p11.31-p11.32 losses. A total of 237 and 101 luminal B-specific candidate oncogenes and tumor suppressor genes (TSGs) presented a deregulated expression in relation with their CNAs, including 11 genes previously reported associated with endocrine resistance. Interestingly, 88% of the potential TSGs are located within chromosome arm 6q, and seven candidate oncogenes are potential therapeutic targets. A total of 100 candidate oncogenes were validated in a public series of 5,765 BCs and the overexpression of 67 of these was associated with poor survival in luminal tumors. Twenty-four genes presented a deregulated expression in relation with a high DNA methylation level. FOXO3, PIK3CA and TP53 were the most frequent mutated genes among the nine tested. In a meta-analysis of next-generation sequencing data in 875 BCs, KCNB2 mutations were associated with luminal B cases while candidate TSGs MDN1 (6q15) and UTRN (6q24), were mutated in this subtype. In conclusion, we have reported luminal B candidate genes that may play a role in the development and/or hormone resistance of this aggressive subtype.     

Robust Gene-Gene Interaction Analysis in Genome Wide Association Studies

Genome-wide association studies (GWAS) have successfully discovered hundreds of associations between genetic variants and complex traits. Most GWAS have focused on the identification of single variants. It has been shown that most of the variants that were discovered by GWAS could only partially explain disease heritability. The explanation for this missing heritability is generally believed to be gene-gene (GG) or gene-environment (GE) interactions and other structural variants. Generalized multifactor dimensionality reduction (GMDR) has been proven to be reasonably powerful in detecting GG and GE interactions; however, its performance has been found to decline when outlying quantitative traits are present. This paper proposes a robust GMDR estimation method (based on the L-estimator and M-estimator estimation methods) in an attempt to reduce the effects caused by outlying traits. A comparison of robust GMDR with the original MDR based on simulation studies showed the former method to outperform the latter. The performance of robust GMDR is illustrated through a real GWA example consisting of 8,577 samples from the Korean population using the Homeostasis Model Assessment of Insulin Resistance (HOMA-IR) level as a phenotype. Robust GMDR identified the KCNH1 gene to have strong interaction effects with other genes on the function of insulin secretion.     

Birth Cohort,Age, and Sex Strongly Modulate Effects of Lipid Risk Alleles Identified in Genome-Wide Association Studies

Insights into genetic origin of diseases and related traits could substantially impact strategies for improving human health. The results of genome-wide association studies (GWAS) are often positioned as discoveries of unconditional risk alleles of complex health traits. We re-analyzed the associations of single nucleotide polymorphisms (SNPs) associated with total cholesterol (TC) in a large-scale GWAS meta-analysis. We focused on three generations of genotyped participants of the Framingham Heart Study (FHS). We show that the effects of all ten directly-genotyped SNPs were clustered in different FHS generations and/or birth cohorts in a sex-specific or sex-unspecific manner. The sample size and procedure-therapeutic issues play, at most, a minor role in this clustering. An important result was clustering of significant associations with the strongest effects in the youngest, or 3^rd Generation, cohort. These results imply that an assumption of unconditional connections of these SNPs with TC is generally implausible and that a demographic perspective can substantially improve GWAS efficiency. The analyses of genetic effects in age-matched samples suggest a role of environmental and age-related mechanisms in the associations of different SNPs with TC. Analysis of the literature supports systemic roles for genes for these SNPs beyond those related to lipid metabolism. Our analyses reveal strong antagonistic effects of rs2479409 (the PCSK9 gene) that cautions strategies aimed at targeting this gene in the next generation of lipid drugs. Our results suggest that standard GWAS strategies need to be advanced in order to appropriately address the problem of genetic susceptibility to complex traits that is imperative for translation to health care.     

Comparison of 6q25 Breast Cancer Hits from Asian and European Genome Wide Association Studies in the Breast Cancer Association Consortium (BCAC)

The 6q25.1 locus was first identified via a genome-wide association study (GWAS) in Chinese women and marked by single nucleotide polymorphism (SNP) rs2046210, approximately 180 Kb upstream of ESR1. There have been conflicting reports about the association of this locus with breast cancer in Europeans, and a GWAS in Europeans identified a different SNP, tagged here by rs12662670. We examined the associations of both SNPs in up to 61,689 cases and 58,822 controls from forty-four studies collaborating in the Breast Cancer Association Consortium, of which four studies were of Asian and 39 of European descent. Logistic regression was used to estimate odds ratios (OR) and 95% confidence intervals (CI). Case-only analyses were used to compare SNP effects in Estrogen Receptor positive (ER+) versus negative (ER−) tumours. Models including both SNPs were fitted to investigate whether the SNP effects were independent. Both SNPs are significantly associated with breast cancer risk in both ethnic groups. Per-allele ORs are higher in Asian than in European studies [rs2046210: OR (A/G) = 1.36 (95% CI 1.26–1.48), p = 7.6×10⁻¹⁴ in Asians and 1.09 (95% CI 1.07–1.11), p = 6.8×10⁻¹⁸ in Europeans. rs12662670: OR (G/T) = 1.29 (95% CI 1.19–1.41), p = 1.2×10⁻⁹ in Asians and 1.12 (95% CI 1.08–1.17), p = 3.8×10⁻⁹ in Europeans]. SNP rs2046210 is associated with a significantly greater risk of ER− than ER+ tumours in Europeans [OR (ER−) = 1.20 (95% CI 1.15–1.25), p = 1.8×10⁻¹⁷ versus OR (ER+) = 1.07 (95% CI 1.04–1.1), p = 1.3×10⁻⁷, p_{heterogeneity} = 5.1×10⁻⁶]. In these Asian studies, by contrast, there is no clear evidence of a differential association by tumour receptor status. Each SNP is associated with risk after adjustment for the other SNP. These results suggest the presence of two variants at 6q25.1 each independently associated with breast cancer risk in Asians and in Europeans. Of these two, the one tagged by rs2046210 is associated with a greater risk of ER− tumours.     

Nanog 基因在卵巢癌和卵巢肿瘤干细胞中的表达及意义

目的:探讨胚胎干细胞关键因子Nanog基因mRNA及其蛋白在卵巢癌和卵巢癌肿瘤干细胞中的表达及意义。方法:选取10例正常卵巢上皮组织、10例卵巢良性肿瘤及60例卵巢癌组织,采用逆转录酶-聚合酶链反应(RT-PCR)方法和免疫组织化学PV-6000两步法检测Nanog mRNA和蛋白表达水平;采用无血清悬浮培养法从SKOV-3卵巢癌细胞株中分离培养肿瘤干细胞,流式细胞术鉴定肿瘤干细胞CD117表达,采用RT-PCR和Western Blot方法检测SKOV-3卵巢癌细胞及肿瘤干细胞中NanogmRNA及其蛋白的表达水平。结果:Nanog mRNA在卵巢癌组织中的表达水平均高于正常卵巢组织和卵巢良性肿瘤组织(P<0.05);Nanog mRNA在不同分化程度及临床分期的卵巢癌组织中表达水平不同,低分化组高于高分化组(P<0.05);III-IV期高于I-II期(P<0.05);免疫组化结果同RT-PCR。从SKOV-3卵巢癌细胞株中成功分离出肿瘤干细胞,SKOV-3卵巢癌细胞和肿瘤干细胞Nanog mRNA相对含量分别为0.6044±0.0368,0.8736±0.0537,差异具有统计学意义(P<0.05),两种细胞Nanog蛋白相对含量分别为0.6364±0.0169 1.2788±0.0314,差别具有统计学意义(P<0.05)。结论:Nanog基因在卵巢癌组织和SKOV-3细胞系中均高表达,其在组织中的表达强度与临床分期及病理分级关系密切,且在肿瘤干细胞中表达高于一般卵巢癌细胞,其与卵巢癌的发生发展关系密切,可能是卵巢癌干细胞的表面标志物,有望成为新的标志物。     

A Ten-MicroRNA Signature Identified from a Genome-Wide MicroRNA Expression Profiling in Human Epithelial Ovarian Cancer

Epithelial ovarian cancer (EOC) is the most common gynecologic malignancy. To identify the micro-ribonucleic acids (miRNAs) expression profile in EOC tissues that may serve as a novel diagnostic biomarker for EOC detection, the expression of 1722 miRNAs from 15 normal ovarian tissue samples and 48 ovarian cancer samples was profiled by using a quantitative real-time polymerase chain reaction (qRT-PCR) assay. A ten-microRNA signature (hsa-miR-1271-5p, hsa-miR-574-3p, hsa-miR-182-5p, hsa-miR-183-5p, hsa-miR-96-5p, hsa-miR-15b-5p, hsa-miR-182-3p, hsa-miR-141-5p, hsa-miR-130b-5p, and hsa-miR-135b-3p) was identified to be able to distinguish human ovarian cancer tissues from normal tissues with 97% sensitivity and 92% specificity. Two miRNA clusters of miR183-96-183 (miR-96-5p, and miR-182, miR183) and miR200 (miR-141-5p, miR200a, b, c and miR429) are significantly up-regulated in ovarian cancer tissue samples compared to those of normal tissue samples, suggesting theses miRNAs may be involved in ovarian cancer development.     

Association between Vitamin D Receptor Gene Polymorphisms and Breast Cancer Risk: A Meta-Analysis of 39 Studies

Background

The associations between vitamin D receptor (VDR) gene polymorphisms and breast cancer risk were comprehensively investigated to clarify issues that remain controversial.

Methodology/Principal Findings

An electronic search was conducted of several databases, including PubMed, the Cochrane library, Web of Science, EMBASE, CBM and CNKI, for papers that describe the association between Fok1, poly-A repeat, Bsm1, Taq1 or Apa1 polymorphisms of the VDR gene and breast cancer risk. Summary odds ratios and 95% confidence intervals (CI) were estimated based on a fixed-effect model (FEM) or random-effect model (REM), depending on the absence or presence of significant heterogeneity. A total of 39 studies met the inclusion criteria. A meta-analysis of high-quality studies showed that the Fok1 polymorphism of the VDR gene was associated with an increased risk of breast cancer (ff vs. Ff+FF, OR: 1.09, 95%CI: 1.02 to 1.16, p = 0.007). No significant associations were observed between the other polymorphisms and breast cancer risk. No positive results were detected by pooling the results of all relevant studies.

Conclusion

A meta-analysis of high-quality studies demonstrated that the Fok1 polymorphism of the VDR gene was closely associated with breast cancer risk.     

Association between Gene Polymorphism of Manganese Superoxide Dismutase and Prostate Cancer Risk

Manganese superoxide dismutase (MnSOD) is the most effective antioxidant enzyme in mitochondria and protects cells from reactive oxygen species‐induced oxidative damage. The aim of this study was to investigate the association between MnSOD Ala‐9Val gene polymorphism and prostate cancer (PCa) risk in Turkish men with prostate cancer. 33 patients with PCa and 81 control individuals were included in the study. We observed an association between MnSOD Ala/Ala frequency and a higher PCa risk. In addition, we found that the increased risk of early‐onset PCa (under age of 65) in the men homozygous for Ala allele was higher than the men homozygous for Val allele. However, we determined that MnSOD Ala‐9Val genotype was not associated with the aggressiveness of the disease. The results of our study suggest that MnSOD Ala/Ala genotype may influence on early‐onset of PCa patients, but no effect on subsequent development of the disease in Turkish men. However, our study has a limitation that is small numbers of individuals for cases and controls. Therefore, the presented study limited our statistical power to fully investigate the gene polymorphism on cancer risk. © 2013 Wiley Periodicals, Inc. J BiochemMol Toxicol 27:213‐218, 2013; View this article online at wileyonlinelibrary.com . DOI 10.1002/jbt.21472     

Association between DNA Damage Response and Repair Genes and Risk of Invasive Serous Ovarian Cancer

Background

We analyzed the association between 53 genes related to DNA repair and p53-mediated damage response and serous ovarian cancer risk using case-control data from the North Carolina Ovarian Cancer Study (NCOCS), a population-based, case-control study.

Methods/Principal Findings

The analysis was restricted to 364 invasive serous ovarian cancer cases and 761 controls of white, non-Hispanic race. Statistical analysis was two staged: a screen using marginal Bayes factors (BFs) for 484 SNPs and a modeling stage in which we calculated multivariate adjusted posterior probabilities of association for 77 SNPs that passed the screen. These probabilities were conditional on subject age at diagnosis/interview, batch, a DNA quality metric and genotypes of other SNPs and allowed for uncertainty in the genetic parameterizations of the SNPs and number of associated SNPs. Six SNPs had Bayes factors greater than 10 in favor of an association with invasive serous ovarian cancer. These included rs5762746 (median OR(odds ratio)_{per allele} = 0.66; 95% credible interval (CI) = 0.44–1.00) and rs6005835 (median OR_{per allele}  = 0.69; 95% CI  = 0.53–0.91) in CHEK2, rs2078486 (median OR_{per allele}  = 1.65; 95% CI = 1.21–2.25) and rs12951053 (median OR_{per allele}  = 1.65; 95% CI = 1.20–2.26) in TP53, rs411697 (median OR _{rare homozygote}  = 0.53; 95% CI  = 0.35 – 0.79) in BACH1 and rs10131 (median OR _{rare homozygote}  =  not estimable) in LIG4. The six most highly associated SNPs are either predicted to be functionally significant or are in LD with such a variant. The variants in TP53 were confirmed to be associated in a large follow-up study.

Conclusions/Significance

Based on our findings, further follow-up of the DNA repair and response pathways in a larger dataset is warranted to confirm these results.     

Magnitude of Stratification in Human Populations and Impacts on Genome Wide Association Studies

Genome-wide association studies (GWAS) may be biased by population stratification (PS). We conducted empirical quantification of the magnitude of PS among human populations and its impact on GWAS. Liver tissues were collected from 979, 59 and 49 Caucasian Americans (CA), African Americans (AA) and Hispanic Americans (HA), respectively, and genotyped using Illumina650Y (Ilmn650Y) arrays. RNA was also isolated and hybridized to Agilent whole-genome gene expression arrays. We propose a new method (i.e., hgdp-eigen) for detecting PS by projecting genotype vectors for each sample to the eigenvector space defined by the Human Genetic Diversity Panel (HGDP). Further, we conducted GWAS to map expression quantitative trait loci (eQTL) for the ∼40,000 liver gene expression traits monitored by the Agilent arrays. HGDP-eigen performed similarly to the conventional self-eigen methods in capturing PS. However, leveraging the HGDP offered a significant advantage in revealing the origins, directions and magnitude of PS. Adjusting for eigenvectors had minor impacts on eQTL detection rates in CA. In contrast, for AA and HA, adjustment dramatically reduced association findings. At an FDR = 10%, we identified 65 eQTLs in AA with the unadjusted analysis, but only 18 eQTLs after the eigenvector adjustment. Strikingly, 55 out of the 65 unadjusted AA eQTLs were validated in CA, indicating that the adjustment procedure significantly reduced GWAS power. A number of the 55 AA eQTLs validated in CA overlapped with published disease associated SNPs. For example, rs646776 and rs10903129 have previously been associated with lipid levels and coronary heart disease risk, however, the rs10903129 eQTL was missed in the eigenvector adjusted analysis.     

葡萄NIP基因家族的鉴定与表达分析

类根瘤菌26膜内在蛋白(nodulin 26 like intrinsic proteins,NIPs)是水通道蛋白的亚类,在植物营养获取和胁迫应答过程中发挥着重要作用。该研究利用多种生物信息学软件,对葡萄NIP家族基因进行分析,并采用RT PCR方法克隆得到4个NIP家族基因,利用qRT PCR方法分析非生物胁迫下NIP基因的表达特征。结果显示：(1)在葡萄基因组中,共鉴定到8个NIP基因,分布于葡萄4条染色体上,主要定位在质膜中;结构上含有6个跨膜结构域和两个典型的保守结构域NPA;氨基酸序列中存在很多个可能的磷酸化位点。(2)进化分析表明葡萄和拟南芥NIP基因具有较高的同源性,基因结构包含外显子数4~6个,保守基序种类和数量相似;基因启动子上游2 kb包含多种应答逆境和激素的顺式调控元件,其数量差异可能与基因本身功能相关。(3)NIP家族基因在不同组织中表达水平差异较大,多数成员在叶中表达水平较高,在茎中较低;成功克隆得到4个葡萄VvNIP基因,其长度分别为789 bp、606 bp、897 bp、789 bp,分别编码262、201、298、293个氨基酸。(4)qRT PCR结果显示,不同胁迫处理下NIP基因在葡萄叶片中的表达水平不同：低温处理下葡萄NIP基因大多呈显著下调表达;盐胁迫下,除VvNIP2 1、VvNIP4 2外其余家族基因均呈下调表达;干旱胁迫下VvNIP4 2显著上调。研究表明,VvNIP基因对多种胁迫均有响应,为葡萄逆境胁迫机制研究提供了参考。     

Utilizing Twins as Controls for Non-Twin Case-Materials in Genome Wide Association Studies

Twin registries around the globe have collected DNA samples from large numbers of monozygotic and dizygotic twins. The twin sample collections are frequently used as controls in disease-specific studies together with non-twins. This approach is unbiased under the hypothesis that twins and singletons are comparable in terms of allele frequencies; i.e. there are no genetic variants associated with being a twin per se. To test this hypothesis we performed a genome-wide association study comparing the allele frequency of 572,352 single nucleotide polymorphisms (SNPs) in 1,413 monozygotic (MZ) and 5,451 dizygotic (DZ) twins with 3,720 healthy singletons. Twins and singletons have been genotyped using the same platform. SNPs showing association with being a twin at P-value < 1 × 10^-5 were selected for replication analysis in 1,492 twins (463 MZ and 1,029 DZ) and 1,880 singletons from Finland. No SNPs reached genome-wide significance (P-value < 5 × 10^-8) in the main analysis combining MZ and DZ twins. In a secondary analysis including only DZ twins two SNPs (rs2033541 close to ADAMTSL1 and rs4149283 close to ABCA1) were genome-wide significant after meta-analysis with the Finnish population. The estimated proportion of variance on the liability scale explained by all SNPs was 0.08 (P-value=0.003) when MZ and DZ were considered together and smaller for MZ (0.06, P-value=0.10) compared to DZ (0.09, P-value=0.003) when analyzed separately. In conclusion, twins and singletons can be used in genetic studies together with general population samples without introducing large bias. Further research is needed to explore genetic variances associated with DZ twinning.     

Association between the XRCC1 Polymorphisms and Thyroid Cancer Risk: A Meta-Analysis from Case-Control Studies

Background

The previous published data on the association between the X-ray repair cross-conplementation group 1 (XRCC1) polymorphisms and thyroid cancer risk remained controversial. Hence, we performed a meta-analysis on all available studies that provided 1729 cases and 3774 controls (from 11 studies) for XRCC1 Arg399Gln, 1040 cases and 2487 controls for Arg194Trp (from 7 studies), and 1432 cases and 3356 controls for Arg280His (from 8 studies).

Methodology/Principal Findings

PubMed, CNKI, and EMBASE database were searched to identify relevant studies. Overall, no significant association was found between XRCC1 Arg399Gln (recessive model: OR = 0.95, 95% CI = 0.77–1.15; dominant model: OR = 0.89, 95% CI = 0.75–1.05; homozygote model: OR = 0.92, 95% CI = 0.69–1.23; Heterozygote model: OR = 0.91, 95% CI = 0.80–1.03; additive model: OR = 0.93, 95% CI = 0.81–1.07), Arg194Trp (recessive model: OR = 1.41, 95% CI = 0.62–3.23; dominant model: OR = 1.01, 95% CI = 0.77–1.34; homozygote model: OR = 1.42, 95% CI = 0.55–3.67; Heterozygote model: OR = 1.03, 95% CI = 0.85–1.26; additive model: OR = 1.08, 95% CI = 0.81–1.42), and Arg280His (recessive model: OR = 1.08, 95% CI = 0.56–2.10; dominant model: OR = 1.01, 95% CI = 0.84–1.22; homozygote model: OR = 1.00, 95% CI = 0.51–1.96; Heterozygote model: OR = 1.04, 95% CI = 0.75–1.42; additive model: OR = 1.03, 95% CI = 0.86–1.23) and thyroid cancer risk when all the eligible studies were pooled into the meta-analysis. In the further stratified and sensitivity analyses, significant association was still not found in these three genetic polymorphisms.

Conclusions/Significance

In summary, this meta-analysis indicates that XRCC1 Arg399Gln, Arg280His, and Arg194Trp are not associated with thyroid cancer.     

MultiGWAS: An integrative tool for Genome Wide Association Studies in tetraploid organisms

The genome‐wide association studies (GWASs) are essential to determine the genetic bases of either ecological or economic phenotypic variation across individuals within populations of the model and nonmodel organisms. For this research question, the GWAS replication testing different parameters and models to validate the results'' reproducibility is common. However, straightforward methodologies that manage both replication and tetraploid data are still missing. To solve this problem, we designed the MultiGWAS, a tool that does GWAS for diploid and tetraploid organisms by executing in parallel four software packages, two designed for polyploid data (GWASpoly and SHEsis) and two designed for diploid data (GAPIT and TASSEL). MultiGWAS has several advantages. It runs either in the command line or in a graphical interface; it manages different genotype formats, including VCF. Moreover, it allows control for population structure, relatedness, and several quality control checks on genotype data. Besides, MultiGWAS can test for additive and dominant gene action models, and, through a proprietary scoring function, select the best model to report its associations. Finally, it generates several reports that facilitate identifying false associations from both the significant and the best‐ranked association Single Nucleotide Polymorphisms (SNPs) among the four software packages. We tested MultiGWAS with public tetraploid potato data for tuber shape and several simulated data under both additive and dominant models. These tests demonstrated that MultiGWAS is better at detecting reliable associations than using each of the four software packages individually. Moreover, the parallel analysis of polyploid and diploid software that only offers MultiGWAS demonstrates its utility in understanding the best genetic model behind the SNP association in tetraploid organisms. Therefore, MultiGWAS probed to be an excellent alternative for wrapping GWAS replication in diploid and tetraploid organisms in a single analysis environment.     

A Genome Wide Association Study Links Glutamate Receptor Pathway to Sporadic Creutzfeldt-Jakob Disease Risk

We performed a genome-wide association (GWA) study in 434 sporadic Creutzfeldt-Jakob disease (sCJD) patients and 1939 controls from the United Kingdom, Germany and The Netherlands. The findings were replicated in an independent sample of 1109 sCJD and 2264 controls provided by a multinational consortium. From the initial GWA analysis we selected 23 SNPs for further genotyping in 1109 sCJD cases from seven different countries. Five SNPs were significantly associated with sCJD after correction for multiple testing. Subsequently these five SNPs were genotyped in 2264 controls. The pooled analysis, including 1543 sCJD cases and 4203 controls, yielded two genome wide significant results: rs6107516 (p-value=7.62x10^-9) a variant tagging the prion protein gene (PRNP); and rs6951643 (p-value=1.66x10^-8) tagging the Glutamate Receptor Metabotropic 8 gene (GRM8). Next we analysed the data stratifying by country of origin combining samples from the pooled analysis with genotypes from the 1000 Genomes Project and imputed genotypes from the Rotterdam Study (Total n=12967). The meta-analysis of the results showed that rs6107516 (p-value=3.00x10^-8) and rs6951643 (p-value=3.91x10^-5) remained as the two most significantly associated SNPs. Rs6951643 is located in an intronic region of GRM8, a gene that was additionally tagged by a cluster of 12 SNPs within our top100 ranked results. GRM8 encodes for mGluR8, a protein which belongs to the metabotropic glutamate receptor family, recently shown to be involved in the transduction of cellular signals triggered by the prion protein. Pathway enrichment analyses performed with both Ingenuity Pathway Analysis and ALIGATOR postulates glutamate receptor signalling as one of the main pathways associated with sCJD. In summary, we have detected GRM8 as a novel, non-PRNP, genome-wide significant marker associated with heightened disease risk, providing additional evidence supporting a role of glutamate receptors in sCJD pathogenesis.     

Association between the PARP1 Val762Ala Polymorphism and Cancer Risk: Evidence from 43 Studies

Background

Poly (ADP-ribose) polymerase-1 (PARP-1) plays critical roles in the detection and repair of damaged DNA, as well as cell proliferation and death. Numerous studies have examined the associations between PARP1 Val762Ala (rs1136410 T>C) polymorphism and cancer susceptibility; nevertheless, the findings from different research groups remain controversial.

Methods

We searched literatures from MEDLINE, EMBASE and CBM pertaining to such associations, and then calculated pooled odds ratio (OR) and 95% confidence interval (CI) by using random-effects model. The false-positive report probability (FPRP) analysis was used to confirm the validity of significant findings. Moreover, potential effects of rs1136410 variants on PARP1 mRNA expression were analyzed for three ethnicities by combining data from HapMap (genotype) and SNPexp (mRNA expression).

Results

The final meta-analysis incorporated 43 studies, consisting of 17,351 cases and 22,401 controls. Overall, our results did not suggest significant associations between Ala variant (Ala/Ala or Ala/Val genotype) and cancer risk. However, further stratification analysis showed significantly increased risk for gastric cancer (Ala/Ala vs. Val/Val: OR = 1.56, 95% CI = 1.01–2.42, Ala/Val vs. Val/Val: OR = 1.34, 95% CI = 1.14–1.58, dominant model: OR = 1.41, 95% CI = 1.21–1.65 and Ala vs. Val: OR = 1.29, 95% CI = 1.07–1.55). On the contrary, decreased risk for brain tumor (Ala/Val vs. Val/Val: OR = 0.77, 95% CI = 0.68–0.87, dominant model: OR = 0.77, 95% CI = 0.68–0.87 and Ala vs. Val: OR = 0.82, 95% CI = 0.74–0.91). Additionally, we found that the Ala carriers had a significantly increased risk in all models for Asians. Our mRNA expression data provided further biological evidence to consolidate this finding.

Conclusions

Despite some limitations, this meta-analysis found evidence for an association between the PAPR1 Val762Ala and cancer susceptibility within gastric cancer, brain tumor and Asian subgroups.