Alexandra Dumitriu

Genomewide association study for onset age in Parkinson disease

BackgroundAge at onset in Parkinson disease (PD) is a highly heritable quantitative trait for which a significant genetic influence is supported by multiple segregation analyses. Because genes associated with onset age may represent invaluable therapeutic targets to delay the disease, we sought to identify such genetic modifiers using a genomewide association study in familial PD. There have been previous genomewide association studies (GWAS) to identify genes influencing PD susceptibility, but this is the first to identify genes contributing to the variation in onset age.MethodsInitial analyses were performed using genotypes generated with the Illumina HumanCNV370Duo array in a sample of 857 unrelated, familial PD cases. Subsequently, a meta-analysis of imputed SNPs was performed combining the familial PD data with that from a previous GWAS of 440 idiopathic PD cases. The SNPs from the meta-analysis with the lowest p-values and consistency in the direction of effect for onset age were then genotyped in a replication sample of 747 idiopathic PD cases from the Parkinson Institute Biobank of Milan, Italy.ResultsMeta-analysis across the three studies detected consistent association (p < 1 × 10-5) with five SNPs, none of which reached genomewide significance. On chromosome 11, the SNP with the lowest p-value (rs10767971; p = 5.4 × 10-7) lies between the genes QSER1 and PRRG4. Near the PARK3 linkage region on chromosome 2p13, association was observed with a SNP (rs7577851; p = 8.7 × 10-6) which lies in an intron of the AAK1 gene. This gene is closely related to GAK, identified as a possible PD susceptibility gene in the GWAS of the familial PD cases.ConclusionTaken together, these results suggest an influence of genes involved in endocytosis and lysosomal sorting in PD pathogenesis.

MicroRNAs Located in the Hox Gene Clusters Are Implicated in Huntington's Disease Pathogenesis

Transcriptional dysregulation has long been recognized as central to the pathogenesis of Huntingtons disease (HD). MicroRNAs (miRNAs) represent a major system of post-transcriptional regulation, by either preventing translational initiation or by targeting transcripts for storage or for degradation. Using next-generation miRNA sequencing in prefrontal cortex (Brodmann Area 9) of twelve HD and nine controls, we identified five miRNAs (miR-10b-5p, miR-196a-5p, miR-196b-5p, miR-615-3p and miR-1247-5p) up-regulated in HD at genome-wide significance (FDR q-value<0.05). Three of these, miR-196a-5p, miR-196b-5p and miR-615-3p, were expressed at near zero levels in control brains. Expression was verified for all five miRNAs using reverse transcription quantitative PCR and all but miR-1247-5p were replicated in an independent sample (8HD/8C). Ectopic miR-10b-5p expression in PC12 HTT-Q73 cells increased survival by MTT assay and cell viability staining suggesting increased expression may be a protective response. All of the miRNAs but miR-1247-5p are located in intergenic regions of Hox clusters. Total mRNA sequencing in the same samples identified fifteen of 55 genes within the Hox cluster gene regions as differentially expressed in HD, and the Hox genes immediately adjacent to the four Hox cluster miRNAs as up-regulated. Pathway analysis of mRNA targets of these miRNAs implicated functions for neuronal differentiation, neurite outgrowth, cell death and survival. In regression models among the HD brains, huntingtin CAG repeat size, onset age and age at death were independently found to be inversely related to miR-10b-5p levels. CAG repeat size and onset age were independently inversely related to miR-196a-5p, onset age was inversely related to miR-196b-5p and age at death was inversely related to miR-615-3p expression. These results suggest these Hox-related miRNAs may be involved in neuroprotective response in HD. Recently, miRNAs have shown promise as biomarkers for human diseases and given their relationship to disease expression, these miRNAs are biomarker candidates in HD.

Gene Expression Profiles in Parkinson Disease Prefrontal Cortex Implicate FOXO1 and Genes under Its Transcriptional Regulation

Parkinson disease (PD) is a complex neurodegenerative disorder with largely unknown genetic mechanisms. While the degeneration of dopaminergic neurons in PD mainly takes place in the substantia nigra pars compacta (SN) region, other brain areas, including the prefrontal cortex, develop Lewy bodies, the neuropathological hallmark of PD. We generated and analyzed expression data from the prefrontal cortex Brodmann Area 9 (BA9) of 27 PD and 26 control samples using the 44K One-Color Agilent 60-mer Whole Human Genome Microarray. All samples were male, without significant Alzheimer disease pathology and with extensive pathological annotation available. 507 of the 39,122 analyzed expression probes were different between PD and control samples at false discovery rate (FDR) of 5%. One of the genes with significantly increased expression in PD was the forkhead box O1 (FOXO1) transcription factor. Notably, genes carrying the FoxO1 binding site were significantly enriched in the FDR–significant group of genes (177 genes covered by 189 probes), suggesting a role for FoxO1 upstream of the observed expression changes. Single-nucleotide polymorphisms (SNPs) selected from a recent meta-analysis of PD genome-wide association studies (GWAS) were successfully genotyped in 50 out of the 53 microarray brains, allowing a targeted expression–SNP (eSNP) analysis for 52 SNPs associated with PD affection at genome-wide significance and the 189 probes from FoxO1 regulated genes. A significant association was observed between a SNP in the cyclin G associated kinase (GAK) gene and a probe in the spermine oxidase (SMOX) gene. Further examination of the FOXO1 region in a meta-analysis of six available GWAS showed two SNPs significantly associated with age at onset of PD. These results implicate FOXO1 as a PD–relevant gene and warrant further functional analyses of its transcriptional regulatory mechanisms.

Copy Number Variation in Familial Parkinson Disease

Copy number variants (CNVs) are known to cause Mendelian forms of Parkinson disease (PD), most notably in SNCA and PARK2. PARK2 has a recessive mode of inheritance; however, recent evidence demonstrates that a single CNV in PARK2 (but not a single missense mutation) may increase risk for PD. We recently performed a genome-wide association study for PD that excluded individuals known to have either a LRRK2 mutation or two PARK2 mutations. Data from the Illumina370Duo arrays were re-clustered using only white individuals with high quality intensity data, and CNV calls were made using two algorithms, PennCNV and QuantiSNP. After quality assessment, the final sample included 816 cases and 856 controls. Results varied between the two CNV calling algorithms for many regions, including the PARK2 locus (genome-wide p = 0.04 for PennCNV and p = 0.13 for QuantiSNP). However, there was consistent evidence with both algorithms for two novel genes, USP32 and DOCK5 (empirical, genome-wide p-values<0.001). PARK2 CNVs tended to be larger, and all instances that were molecularly tested were validated. In contrast, the CNVs in both novel loci were smaller and failed to replicate using real-time PCR, MLPA, and gel electrophoresis. The DOCK5 variation is more akin to a VNTR than a typical CNV and the association is likely caused by artifact due to DNA source. DNA for all the cases was derived from whole blood, while the DNA for all controls was derived from lymphoblast cell lines. The USP32 locus contains many SNPs with low minor allele frequency leading to a loss of heterozygosity that may have been spuriously interpreted by the CNV calling algorithms as support for a deletion. Thus, only the CNVs within the PARK2 locus could be molecularly validated and associated with PD susceptibility.

Cyclin-G-associated kinase modifies α-synuclein expression levels and toxicity in Parkinson's disease: results from the GenePD Study

Although family history is a well-established risk factor for Parkinsons disease (PD), fewer than 5% of PD cases can be attributed to known genetic mutations. The etiology for the remainder of PD cases is unclear; however, neuronal accumulation of the protein α-synuclein is common to nearly all patients, implicating pathways that influence α-synuclein in PD pathogenesis. We report a genome-wide significant association (P = 3.97 × 10(-8)) between a polymorphism, rs1564282, in the cyclin-G-associated kinase (GAK) gene and increased PD risk, with a meta-analysis odds ratio of 1.48. This association result is based on the meta-analysis of three publicly available PD case-control genome-wide association study and genotyping from a new, independent Italian cohort. Microarray expression analysis of post-mortem frontal cortex from PD and control brains demonstrates a significant association between rs1564282 and higher α-synuclein expression, a known cause of early onset PD. Functional knockdown of GAK in cell culture causes a significant increase in toxicity when α-synuclein is over-expressed. Furthermore, knockdown of GAK in rat primary neurons expressing the A53T mutation of α-synuclein, a well-established model for PD, decreases cell viability. These observations provide evidence that GAK is associated with PD risk and suggest that GAK and α-synuclein interact in a pathway involved in PD pathogenesis. The GAK protein, a serine/threonine kinase, belongs to a family of proteins commonly targeted for drug development. This, combined with GAKs observed relationship to the levels of α-synuclein expression and toxicity, suggests that the protein is an attractive therapeutic target for the treatment of PD.

Decreased glutamic acid decarboxylase mRNA expression in prefrontal cortex in Parkinson's disease.

Parkinsons disease (PD) patients typically suffer from motor disorders but mild to severe cognitive deficits can also be present. Neuropathology of PD primarily involves loss of dopaminergic neurons in the substantia nigra, pars compacta, although more widespread pathology from the brainstem to the cerebral cortex occurs at different stages of the disease. Cognitive deficits in PD are thought to involve the cerebral cortex, and imaging studies have identified the dorsolateral prefrontal cortex (DLPFC) as a possible site for some of the symptoms. GABAergic neurons in the cerebral cortex play a key role in the modulation of pyramidal neurons and alterations in muscimol binding to GABA(A) receptors have been reported in Brodmann area 9 (BA9) of the prefrontal cortex in PD patients (Nishino et al., 1988). In order to further assess the likelihood that GABAergic activity is altered in the prefrontal cortex in PD, gene expression of the 67 kilodalton isoform of the GABA-synthesizing enzyme, glutamic acid decarboxylase (GAD67 encoded by the GAD1 gene), was examined in BA9 of post-mortem brains from 19 patients and 20 controls using isotopic in situ hybridization histochemistry. GAD67 mRNA labeling was examined and quantified on X-ray films and emulsion radioautographs. We show that GAD67 mRNA labeling is significantly lower in PD compared to control cases. Analysis of emulsion radioautographs indicates that GAD67 mRNA labeling is decreased in individual neurons and is not paralleled by a decrease in the number of GAD67 mRNA-labeled neurons. Analysis of expression data from a microarray study performed in 29 control and 33 PD samples from BA9 confirms that GAD67 expression is decreased in PD. Another finding from the microarray study is a negative relationship between GAD67 mRNA expression and age at death. Altogether, the results support the possibility that GABAergic neurotransmission is impaired in the DLPFC in PD, an effect that may be involved in some of the behavioral deficits associated with the disease.

Evaluation of Parkinson Disease Risk Variants as Expression-QTLs

The recent Parkinson Disease GWAS Consortium meta-analysis and replication study reports association at several previously confirmed risk loci SNCA, MAPT, GAK/DGKQ, and HLA and identified a novel risk locus at RIT2. To further explore functional consequences of these associations, we investigated modification of gene expression in prefrontal cortex brain samples of pathologically confirmed PD cases (N = 26) and controls (N = 24) by 67 associated SNPs in these 5 loci. Association between the eSNPs and expression was evaluated using a 2-degrees of freedom test of both association and difference in association between cases and controls, adjusted for relevant covariates. SNPs at each of the 5 loci were tested for cis-acting effects on all probes within 250 kb of each locus. Trans-effects of the SNPs on the 39,122 probes passing all QC on the microarray were also examined. From the analysis of cis-acting SNP effects, several SNPs in the MAPT region show significant association to multiple nearby probes, including two strongly correlated probes targeting the gene LOC644246 and the duplicated genes LRRC37A and LRRC37A2, and a third uncorrelated probe targeting the gene DCAKD. Significant cis-associations were also observed between SNPs and two probes targeting genes in the HLA region on chromosome 6. Expanding the association study to examine trans effects revealed an additional 23 SNP-probe associations reaching statistical significance (p<2.8×10−8) including SNPs from the SNCA, MAPT and RIT2 regions. These findings provide additional context for the interpretation of PD associated SNPs identified in recent GWAS as well as potential insight into the mechanisms underlying the observed SNP associations.

Postmortem Interval Influences α-Synuclein Expression in Parkinson Disease Brain.

Duplications and triplications of the α-synuclein (SNCA) gene increase risk for PD, suggesting increased expression levels of the gene to be associated with increased PD risk. However, past SNCA expression studies in brain tissue report inconsistent results. We examined expression of the full-length SNCA transcript (140 amino acid protein isoform), as well as total SNCA mRNA levels in 165 frontal cortex samples (101 PD, 64 control) using quantitative real-time polymerase chain reaction. Additionally, we evaluated the relationship of eight SNPs in both 5′ and 3′ regions of SNCA with the gene expression levels. The association between postmortem interval (PMI) and SNCA expression was different for PD and control samples: SNCA expression decreased with increasing PMI in cases, while staying relatively constant in controls. For short PMI, SNCA expression was increased in PD relative to control samples, whereas for long PMI, SNCA expression in PD was decreased relative to control samples.

Comparative multi-goal tradeoffs in systems engineering of microbial metabolism

BackgroundMetabolic engineering design methodology has evolved from using pathway-centric, random and empirical-based methods to using systems-wide, rational and integrated computational and experimental approaches. Persistent during these advances has been the desire to develop design strategies that address multiple simultaneous engineering goals, such as maximizing productivity, while minimizing raw material costs.ResultsHere, we use constraint-based modeling to systematically design multiple combinations of medium compositions and gene-deletion strains for three microorganisms (Escherichia coli, Saccharomyces cerevisiae, and Shewanella oneidensis) and six industrially important byproducts (acetate, D-lactate, hydrogen, ethanol, formate, and succinate). We evaluated over 435 million simulated conditions and 36 engineering metabolic traits, including product rates, costs, yields and purity.ConclusionsThe resulting metabolic phenotypes can be classified into dominant clusters (meta-phenotypes) for each organism. These meta-phenotypes illustrate global phenotypic variation and sensitivities, trade-offs associated with multiple engineering goals, and fundamental differences in organism-specific capabilities. Given the increasing number of sequenced genomes and corresponding stoichiometric models, we envisage that the proposed strategy could be extended to address a growing range of biological questions and engineering applications.

The 4p16.3 Parkinson Disease Risk Locus Is Associated with GAK Expression and Genes Involved with the Synaptic Vesicle Membrane

Genome-wide association studies (GWAS) have identified the GAK/DGKQ/IDUA region on 4p16.3 among the top three risk loci for Parkinson’s disease (PD), but the specific gene and risk mechanism are unclear. Here, we report transcripts containing the 3’ clathrin-binding domain of GAK identified by RNA deep-sequencing in post-mortem human brain tissue as having increased expression in PD. Furthermore, carriers of 4p16.3 PD GWAS risk SNPs show decreased expression of one of these transcripts, GAK25 (Gencode Transcript 009), which correlates with the expression of genes functioning in the synaptic vesicle membrane. Together, these findings provide strong evidence for GAK clathrin-binding- and J-domain transcripts’ influence on PD pathogenicity, and for a role for GAK in regulating synaptic function in PD.