The pathogenic process responsible for the loss of dopaminergic neurons within the substantia nigra of patients with Parkinson disease (PD) is poorly understood. Current research supports the involvement of fibroblast growth factor (FGF20) in the survival of dopaminergic cells. FGF20 is a neurotrophic factor that is preferentially expressed within the substantia nigra of rat brain. The human homologue has been mapped to 8p21.3-8p22, which is within an area of PD linkage revealed through our published genomic screen. To test whether FGF20 influences risk of PD, we genotyped five single-nucleotide polymorphisms (SNPs) lying within the FGF20 gene, in a large family study. We analyzed our sample (644 families) through use of the pedigree disequilibrium test (PDT), the genotype PDT, the multilocus-genotype PDT, and the family-based association test to assess association between risk of PD and alleles, genotypes, multilocus genotypes, and haplotypes. We discovered a highly significant association of PD with one intronic SNP, rs1989754 (P=.0006), and two SNPs, rs1721100 (P=.02) and ss20399075 (P=.0008), located in the 3' regulatory region in our overall sample. Furthermore, we detected a haplotype (A-G-C-C-T) that is positively associated with risk of PD (P=.0003), whereas a second haplotype (A-G-G-G-C) was found to be negatively associated with risk of PD (P=.0009). Our results strongly support FGF20 as a risk factor for PD.

The functional structure of all biologically active molecules is dependent on intra- and inter-molecular interactions. This is especially evident for RNA molecules whose functionality, maturation, and regulation require formation of correct secondary structure through encoded base-pairing interactions. Unfortunately, intra- and inter-molecular base-pairing information is lacking for most RNAs. Here, we marry classical nuclease-based structure mapping techniques with high-throughput sequencing technology to interrogate all base-paired RNA in Arabidopsis thaliana and identify ∼200 new small (sm)RNA-producing substrates of RNA-DEPENDENT RNA POLYMERASE6. Our comprehensive analysis of paired RNAs reveals conserved functionality within introns and both 5' and 3' untranslated regions (UTRs) of mRNAs, as well as a novel population of functional RNAs, many of which are the precursors of smRNAs. Finally, we identify intra-molecular base-pairing interactions to produce a genome-wide collection of RNA secondary structure models. Although our methodology reveals the pairing status of RNA molecules in the absence of cellular proteins, previous studies have demonstrated that structural information obtained for RNAs in solution accurately reflects their structure in ribonucleoprotein complexes. Furthermore, our identification of RNA-DEPENDENT RNA POLYMERASE6 substrates and conserved functional RNA domains within introns and both 5' and 3' untranslated regions (UTRs) of mRNAs using this approach strongly suggests that RNA molecules are correctly folded into their secondary structure in solution. Overall, our findings highlight the importance of base-paired RNAs in eukaryotes and present an approach that should be widely applicable for the analysis of this key structural feature of RNA.

Tyrosine hydroxylase (TH) enzyme is a rate limiting enzyme in dopamine biosynthesis. Missense mutation in both alleles of the TH gene is known to cause dopamine-related phenotypes, including dystonia and infantile Parkinsonism. However, it is not clear if single allele mutation in TH modifies the susceptibility to the adult form of Parkinson disease (PD). We reported a novel deletion of entire TH gene in an adult with PD. The deletion was first identified by copy number variation (CNV) analysis in a genome-wide association study using Illumina Infinium BeadChips. After screening 635 cases and 642 controls, the deletion was found in one PD case but not in any control. The deletion was confirmed by multiple quantitative PCR (qPCR) assays. There is no additional exonic single nucleotide variant in the one copy of TH gene of the patient. The patient has an age-at-onset of 54 years, no evidence for dystonia, and was responsive to L-DOPA. This case supports the importance of the TH gene in PD pathogenesis and raises more attention to rare variants in candidate genes being a risk factor for Parkinson disease.

Alzheimer's disease (AD) and related dementias are a major public health challenge and present a therapeutic imperative for which we need additional insight into molecular pathogenesis. We performed a genome-wide association study and analysis of known genetic risk loci for AD dementia using neuropathologic data from 4,914 brain autopsies. Neuropathologic data were used to define clinico-pathologic AD dementia or controls, assess core neuropathologic features of AD (neuritic plaques, NPs; neurofibrillary tangles, NFTs), and evaluate commonly co-morbid neuropathologic changes: cerebral amyloid angiopathy (CAA), Lewy body disease (LBD), hippocampal sclerosis of the elderly (HS), and vascular brain injury (VBI). Genome-wide significance was observed for clinico-pathologic AD dementia, NPs, NFTs, CAA, and LBD with a number of variants in and around the apolipoprotein E gene (APOE). GalNAc transferase 7 (GALNT7), ATP-Binding Cassette, Sub-Family G (WHITE), Member 1 (ABCG1), and an intergenic region on chromosome 9 were associated with NP score; and Potassium Large Conductance Calcium-Activated Channel, Subfamily M, Beta Member 2 (KCNMB2) was strongly associated with HS. Twelve of the 21 non-APOE genetic risk loci for clinically-defined AD dementia were confirmed in our clinico-pathologic sample: CR1, BIN1, CLU, MS4A6A, PICALM, ABCA7, CD33, PTK2B, SORL1, MEF2C, ZCWPW1, and CASS4 with 9 of these 12 loci showing larger odds ratio in the clinico-pathologic sample. Correlation of effect sizes for risk of AD dementia with effect size for NFTs or NPs showed positive correlation, while those for risk of VBI showed a moderate negative correlation. The other co-morbid neuropathologic features showed only nominal association with the known AD loci. Our results discovered new genetic associations with specific neuropathologic features and aligned known genetic risk for AD dementia with specific neuropathologic changes in the largest brain autopsy study of AD and related dementias.

Age-related macular degeneration (AMD) is a common cause of blindness in older individuals. To accelerate the understanding of AMD biology and help design new therapies, we executed a collaborative genome-wide association study, including >17,100 advanced AMD cases and >60,000 controls of European and Asian ancestry. We identified 19 loci associated at P < 5 × 10(-8). These loci show enrichment for genes involved in the regulation of complement activity, lipid metabolism, extracellular matrix remodeling and angiogenesis. Our results include seven loci with associations reaching P < 5 × 10(-8) for the first time, near the genes COL8A1-FILIP1L, IER3-DDR1, SLC16A8, TGFBR1, RAD51B, ADAMTS9 and B3GALTL. A genetic risk score combining SNP genotypes from all loci showed similar ability to distinguish cases and controls in all samples examined. Our findings provide new directions for biological, genetic and therapeutic studies of AMD.

The dynamic range of cerebrospinal fluid (CSF) amyloid β (Aβ1-42) measurement does not parallel to cognitive changes in Alzheimer's disease (AD) and cognitively normal (CN) subjects across different studies. Therefore, identifying novel proteins to characterize symptomatic AD samples is important.

The human Werner and Bloom syndromes (WS and BS) are caused by deficiencies in the WRN and BLM RecQ helicases, respectively. WRN, BLM and their Saccharomyces cerevisiae homologue Sgs1, are particularly active in vitro in unwinding G-quadruplex DNA (G4-DNA), a family of non-canonical nucleic acid structures formed by certain G-rich sequences. Recently, mRNA levels from loci containing potential G-quadruplex-forming sequences (PQS) were found to be preferentially altered in sgs1Delta mutants, suggesting that G4-DNA targeting by Sgs1 directly affects gene expression. Here, we extend these findings to human cells. Using microarrays to measure mRNAs obtained from human fibroblasts deficient for various RecQ family helicases, we observe significant associations between loci that are upregulated in WS or BS cells and loci that have PQS. No such PQS associations were observed for control expression datasets, however. Furthermore, upregulated genes in WS and BS showed no or dramatically reduced associations with sequences similar to PQS but that have considerably reduced potential to form intramolecular G4-DNA. These findings indicate that, like Sgs1, WRN and BLM can regulate transcription globally by targeting G4-DNA.

Autism is a neurodevelopmental disorder of complex genetics, characterized by impairment in social interaction and communication, as well as repetitive behavior. Multiple lines of evidence, including alterations in levels of GABA and GABA receptors in autistic patients, indicate that the GABAergic system, which is responsible for synaptic inhibition in the adult brain, may be involved in autism. Previous studies in our lab indicated association of noncoding single nucleotide polymorphisms (SNPs) within a GABA receptor subunit gene on chromosome 4, GABRA4, and interaction between SNPs in GABRA4 and GABRB1 (also on chromosome 4), within Caucasian autism patients. Studies of genetic variation in African-American autism families are rare. Analysis of 557 Caucasian and an independent population of 54 African-American families with 35 SNPs within GABRB1 and GABRA4 strengthened the evidence for involvement of GABRA4 in autism risk in Caucasians (rs17599165, p=0.0015; rs1912960, p=0.0073; and rs17599416, p=0.0040) and gave evidence of significant association in African-Americans (rs2280073, p=0.0287 and rs16859788, p=0.0253). The GABRA4 and GABRB1 interaction was also confirmed in the Caucasian dataset (most significant pair, rs1912960 and rs2351299; p=0.004). Analysis of the subset of families with a positive history of seizure activity in at least one autism patient revealed no association to GABRA4; however, three SNPs within GABRB1 showed significant allelic association; rs2351299 (p=0.0163), rs4482737 (p=0.0339), and rs3832300 (p=0.0253). These results confirmed our earlier findings, indicating GABRA4 and GABRB1 as genes contributing to autism susceptibility, extending the effect to multiple ethnic groups and suggesting seizures as a stratifying phenotype.

Small non-coding RNAs (sncRNAs, <100 nts) are highly abundant RNAs that regulate diverse and often tissue-specific cellular processes by associating with transcription factor complexes or binding to mRNAs. While thousands of sncRNA genes exist in the human genome, no single resource provides searchable, unified annotation, expression and processing information for full sncRNA transcripts and mature RNA products derived from these larger RNAs.

Given the known strong relationship of DNA methylation with environmental exposure, we investigated whether brain regions affected in Parkinson disease (PD) were differentially methylated between PD cases and controls.

To identify rare causal variants underlying known loci that segregate with late-onset Alzheimer's disease (LOAD) in multiplex families.

Carbonic anhydrase-8 (CA8) is an intracellular protein that functions as an allosteric inhibitor of inositol trisphosphate receptor-1 (ITPR1) critical to intracellular Ca++ release, synaptic functions and neuronal excitability. We showed previously that murine nociception and analgesic responses are regulated by the expression of this gene in dorsal root ganglion (DRG) associated with a cis-eQTL. In this report, we identify an exon-level cis-eQTL (rs6471859) that regulates human DRG CA8 alternative splicing, producing a truncated 1,697bp transcript (e.g., CA8-204). Our functional genomic studies show the "G" allele at rs6471859 produces a cryptic 3'UTR splice site regulating expression of CA8-204. We developed constructs to study the expression and function of the naturally occurring CA8-204G transcript (G allele at rs6471859), CA8-204C (C allele at rs6471859 reversion mutation) and CA8-201 (full length transcript). CA8-204G transcript expression occurred predominantly in non-neuronal cells (HEK293), while CA8-204C expression was restricted to neuronal derived cells (NBL) in vitro. CA8-204G produced a stable truncated transcript in HEK293 cells that was barely detectable in NBL cells. We also show CA8-204 produces a stable peptide that inhibits pITPR1 and Ca++ release in HEK293 cells. These results imply homozygous G/G individuals at rs6471859, which are common in the general population, produce exclusively CA8-204G that is barely detectable in neuronal cells. CA8 null mutations that greatly impact neuronal functions are associated with severe forms of spinal cerebellar ataxia, and our data suggest G/G homozygotes should display a similar phenotype. To address this question, we show in vivo using AAV8-FLAG-CA8-204G and AAV8-V5-CA8-201 gene transfer delivered via intra-neural sciatic nerve injection (SN), that these viral constructs are able to transduce DRG cells and produce similar analgesic and anti-hyperalgesic responses to inflammatory pain. Immunohistochemistry (IHC) examinations of DRG tissues further show CA8-204G peptide is expressed in advillin expressing neuronal cells, but to a lesser extent compared to glial cells. These findings explain why G/G homozygotes that exclusively produce this truncated functional peptide in DRG evade a severe phenotype. These genomic studies significantly advance the literature regarding structure-function studies on CA8-ITPR1 critical to calcium signaling pathways, synaptic functioning, neuronal excitability and analgesic responses.

The poor outcomes in infant acute lymphoblastic leukemia (ALL) necessitate new treatments. Here we discover that EIF4E protein is elevated in most cases of infant ALL and test EIF4E targeting by the repurposed antiviral agent ribavirin, which has anticancer properties through EIF4E inhibition, as a potential treatment. We find that ribavirin treatment of actively dividing infant ALL cells on bone marrow stromal cells (BMSCs) at clinically achievable concentrations causes robust proliferation inhibition in proportion with EIF4E expression. Further, we find that ribavirin treatment of KMT2A-rearranged (KMT2A-R) infant ALL cells and the KMT2A-AFF1 cell line RS4:11 inhibits EIF4E, leading to decreases in oncogenic EIF4E-regulated cell growth and survival proteins. In ribavirin-sensitive KMT2A-R infant ALL cells and RS4:11 cells, EIF4E-regulated proteins with reduced levels of expression following ribavirin treatment include MYC, MCL1, NBN, BCL2 and BIRC5. Ribavirin-treated RS4:11 cells exhibit impaired EIF4E-dependent nuclear to cytoplasmic export and/or translation of the corresponding mRNAs, as well as reduced phosphorylation of the p-AKT1, p-EIF4EBP1, p-RPS6 and p-EIF4E signaling proteins. This leads to an S-phase cell cycle arrest in RS4:11 cells corresponding to the decreased proliferation. Ribavirin causes nuclear EIF4E to re-localize to the cytoplasm in KMT2A-AFF1 infant ALL and RS4:11 cells, providing further evidence for EIF4E inhibition. Ribavirin slows increases in peripheral blasts in KMT2A-R infant ALL xenograft-bearing mice. Ribavirin cooperates with chemotherapy, particularly L-asparaginase, in reducing live KMT2A-AFF1 infant ALL cells in BMSC co-cultures. This work establishes that EIF4E is broadly elevated across infant ALL and that clinically relevant ribavirin exposures have preclinical activity and effectively inhibit EIF4E in KMT2A-R cases, suggesting promise in EIF4E targeting using ribavirin as a means of treatment.

Recent high-throughput structure-sensitive genome-wide sequencing-based assays have enabled large-scale studies of RNA structure, and robust transcriptome-wide computational prediction of individual RNA structures across RNA classes from these assays has potential to further improve the prediction accuracy. Here, we describe HiPR, a novel method for RNA structure prediction at single-nucleotide resolution that combines high-throughput structure probing data (DMS-seq, DMS-MaPseq) with a novel probabilistic folding algorithm. On validation data spanning a variety of RNA classes, HiPR often increases accuracy for predicting RNA structures, giving researchers new tools to study RNA structure.

Findings regarding the association between mitochondrial DNA (mtDNA) variants and Alzheimer's disease (AD) are inconsistent.

The explosion of biobank data offers unprecedented opportunities for gene-environment interaction (GxE) studies of complex diseases because of the large sample sizes and the rich collection in genetic and non-genetic information. However, the extremely large sample size also introduces new computational challenges in G×E assessment, especially for set-based G×E variance component (VC) tests, which are a widely used strategy to boost overall G×E signals and to evaluate the joint G×E effect of multiple variants from a biologically meaningful unit (e.g., gene). In this work, we focus on continuous traits and present SEAGLE, a Scalable Exact AlGorithm for Large-scale set-based G×E tests, to permit G×E VC tests for biobank-scale data. SEAGLE employs modern matrix computations to calculate the test statistic and p-value of the GxE VC test in a computationally efficient fashion, without imposing additional assumptions or relying on approximations. SEAGLE can easily accommodate sample sizes in the order of 105, is implementable on standard laptops, and does not require specialized computing equipment. We demonstrate the performance of SEAGLE using extensive simulations. We illustrate its utility by conducting genome-wide gene-based G×E analysis on the Taiwan Biobank data to explore the interaction of gene and physical activity status on body mass index.

Recent Alzheimer's disease (AD) genetics findings from genome-wide association studies (GWAS) span progressively larger and more diverse populations and outcomes. Currently, there is no up-to-date resource providing harmonized and searchable information on all AD genetic associations found by GWAS, nor linking the reported genetic variants and genes with functional and genomic annotations.

Characterization of the genetic landscape of Alzheimer's disease (AD) and related dementias (ADD) provides a unique opportunity for a better understanding of the associated pathophysiological processes. We performed a two-stage genome-wide association study totaling 111,326 clinically diagnosed/'proxy' AD cases and 677,663 controls. We found 75 risk loci, of which 42 were new at the time of analysis. Pathway enrichment analyses confirmed the involvement of amyloid/tau pathways and highlighted microglia implication. Gene prioritization in the new loci identified 31 genes that were suggestive of new genetically associated processes, including the tumor necrosis factor alpha pathway through the linear ubiquitin chain assembly complex. We also built a new genetic risk score associated with the risk of future AD/dementia or progression from mild cognitive impairment to AD/dementia. The improvement in prediction led to a 1.6- to 1.9-fold increase in AD risk from the lowest to the highest decile, in addition to effects of age and the APOE ε4 allele.

Non-coding genetic variants outside of protein-coding genome regions play an important role in genetic and epigenetic regulation. It has become increasingly important to understand their roles, as non-coding variants often make up the majority of top findings of genome-wide association studies (GWAS). In addition, the growing popularity of disease-specific whole-genome sequencing (WGS) efforts expands the library of and offers unique opportunities for investigating both common and rare non-coding variants, which are typically not detected in more limited GWAS approaches. However, the sheer size and breadth of WGS data introduce additional challenges to predicting functional impacts in terms of data analysis and interpretation. This review focuses on the recent approaches developed for efficient, at-scale annotation and prioritization of non-coding variants uncovered in WGS analyses. In particular, we review the latest scalable annotation tools, databases and functional genomic resources for interpreting the variant findings from WGS based on both experimental data and in silico predictive annotations. We also review machine learning-based predictive models for variant scoring and prioritization. We conclude with a discussion of future research directions which will enhance the data and tools necessary for the effective functional analyses of variants identified by WGS to improve our understanding of disease etiology.

Structural variations (SVs) are important contributors to the genetics of numerous human diseases. However, their role in Alzheimer's disease (AD) remains largely unstudied due to challenges in accurately detecting SVs. Here, we analyzed whole-genome sequencing data from the Alzheimer's Disease Sequencing Project (ADSP, N=16,905 subjects) and identified 400,234 (168,223 high-quality) SVs. We found a significant burden of deletions and duplications in AD cases (OR=1.05, P=0.03), particularly for singletons (OR=1.12, P=0.0002) and homozygous events (OR=1.10, P<0.0004). On AD genes, the ultra-rare SVs, including protein-altering SVs in ABCA7, APP, PLCG2, and SORL1, were associated with AD (SKAT-O P=0.004). Twenty-one SVs are in linkage disequilibrium (LD) with known AD-risk variants, e.g., a deletion (chr2:105731359-105736864) in complete LD (R2=0.99) with rs143080277 (chr2:105749599) in NCK2. We also identified 16 SVs associated with AD and 13 SVs associated with AD-related pathological/cognitive endophenotypes. Our findings demonstrate the broad impact of SVs on AD genetics.