Big flies, small repeats: the "Thr-Gly" region of the period gene in Diptera

The region of the clock gene period (per) that encodes a repetitive tract of threonine-glycine (Thr-Gly) pairs has been compared between Dipteran species both within and outside the Drosophilidae. All the non- Drosophilidae sequences in this region are short and present a remarkably stable picture compared to the Drosophilidae, in which the region is much larger and extremely variable, both in size and composition. The accelerated evolution in the repetitive region of the Drosophilidae appears to be mainly due to an expansion of two ancestral repeats, one encoding a Thr-Gly dipeptide and the other a pentapeptide rich in serine, glycine, and asparagine or threonine. In some drosophilids the expansion involves a duplication of the pentapeptide sequence, but in Drosophila pseudoobscura both the dipeptide and the pentapeptide repeats are present in larger numbers. In the nondrosophilids, however, the pentapeptide sequence is represented by one copy and the dipeptide by two copies. These observations fulfill some of the predictions of recent theoretical models that have simulated the evolution of repetitive sequences.      

Genome-wide gene-environment study identifies glutamate receptor gene GRIN2A as a Parkinson's disease modifier gene via interaction with coffee

Our aim was to identify genes that influence the inverse association of coffee with the risk of developing Parkinson''s disease (PD). We used genome-wide genotype data and lifetime caffeinated-coffee-consumption data on 1,458 persons with PD and 931 without PD from the NeuroGenetics Research Consortium (NGRC), and we performed a genome-wide association and interaction study (GWAIS), testing each SNP''s main-effect plus its interaction with coffee, adjusting for sex, age, and two principal components. We then stratified subjects as heavy or light coffee-drinkers and performed genome-wide association study (GWAS) in each group. We replicated the most significant SNP. Finally, we imputed the NGRC dataset, increasing genomic coverage to examine the region of interest in detail. The primary analyses (GWAIS, GWAS, Replication) were performed using genotyped data. In GWAIS, the most significant signal came from rs4998386 and the neighboring SNPs in GRIN2A. GRIN2A encodes an NMDA-glutamate-receptor subunit and regulates excitatory neurotransmission in the brain. Achieving P_2df = 10⁻⁶, GRIN2A surpassed all known PD susceptibility genes in significance in the GWAIS. In stratified GWAS, the GRIN2A signal was present in heavy coffee-drinkers (OR = 0.43; P = 6×10⁻⁷) but not in light coffee-drinkers. The a priori Replication hypothesis that “Among heavy coffee-drinkers, rs4998386_T carriers have lower PD risk than rs4998386_CC carriers” was confirmed: OR_Replication = 0.59, P_Replication = 10⁻³; OR_Pooled = 0.51, P_Pooled = 7×10⁻⁸. Compared to light coffee-drinkers with rs4998386_CC genotype, heavy coffee-drinkers with rs4998386_CC genotype had 18% lower risk (P = 3×10⁻³), whereas heavy coffee-drinkers with rs4998386_TC genotype had 59% lower risk (P = 6×10⁻¹³). Imputation revealed a block of SNPs that achieved P_2df<5×10⁻⁸ in GWAIS, and OR = 0.41, P = 3×10⁻⁸ in heavy coffee-drinkers. This study is proof of concept that inclusion of environmental factors can help identify genes that are missed in GWAS. Both adenosine antagonists (caffeine-like) and glutamate antagonists (GRIN2A-related) are being tested in clinical trials for treatment of PD. GRIN2A may be a useful pharmacogenetic marker for subdividing individuals in clinical trials to determine which medications might work best for which patients.     

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.     

Dissection of a metastatic gene expression signature into distinct components

Background

Metastasis, the process whereby cancer cells spread, is in part caused by an incompletely understood interplay between cancer cells and the surrounding stroma. Gene expression studies typically analyze samples containing tumor cells and stroma. Samples with less than 50% tumor cells are generally excluded, thereby reducing the number of patients that can benefit from clinically relevant signatures.

Results

For a head-neck squamous cell carcinoma (HNSCC) primary tumor expression signature that predicts the presence of lymph node metastasis, we first show that reduced proportions of tumor cells results in decreased predictive accuracy. To determine the influence of stroma on the predictive signature and to investigate the interaction between tumor cells and the surrounding microenvironment, we used laser capture microdissection to divide the metastatic signature into six distinct components based on tumor versus stroma expression and on association with the metastatic phenotype. A strikingly skewed distribution of metastasis associated genes is revealed.

Conclusion

Dissection of predictive signatures into different components has implications for design of expression signatures and for our understanding of the metastatic process. Compared to primary tumors that have not formed metastases, primary HNSCC tumors that have metastasized are characterized by predominant down-regulation of tumor cell specific genes and exclusive up-regulation of stromal cell specific genes. The skewed distribution agrees with poor signature performance on samples that contain less than 50% tumor cells. Methods for reducing tumor composition bias that lead to greater predictive accuracy and an increase in the types of samples that can be included are presented.     

Exploring gene-environment interactions in Parkinson’s disease

The objective of this study was to explore combined effects of four candidate susceptibility genes and two exposures on Parkinson’s disease (PD) risk; namely, α-synuclein (SNCA) promoter polymorphism REP1, microtubule-associated protein tau (MAPT) H1/H2 haplotypes, apolipoprotein E (APOE) ε2/ε3/ε4 polymorphism, ubiquitin carboxy-terminal esterase L1 (UCHL1) S18Y variant, cigarette smoking and caffeinated coffee consumption. 932 PD patients and 664 control subjects from the NeuroGenetics Research Consortium, with complete data on all six factors, were studied. Uniform protocols were used for diagnosis, recruitment, data collection and genotyping. A logistic regression model which included gene-exposure interactions was applied. Likelihood ratio tests (LRTs) were used for significance testing and Bayesian inference was used to estimate odds ratios (ORs). MAPT (P = 0.007), SNCA REP1 (P = 0.012), smoking (P = 0.001), and coffee (P = 0.011) were associated with PD risk. Two novel interactions were detected: APOE with coffee (P = 0.005), and REP1 with smoking (P = 0.021). While the individual main effects were modest, each yielding OR < 1.6, the effects were cumulative, with some combinations reaching OR = 12.6 (95% CI: 5.9–26.8). This study provides evidence for the long-held notion that PD risk is modulated by cumulative and interactive effects of genes and exposures. Furthermore, the study demonstrates that while interaction studies are useful for exploring risk relationships that might otherwise go undetected, results should be interpreted with caution because of the inherent loss of power due to multiple testing. The novel findings of this study that warrant replication are the evidence for interaction of coffee with APOE, and of smoking with REP1 on PD risk. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.