High-throughput sequencing platforms are generating massive amounts of genetic variation data for diverse genomes, but it remains a challenge to pinpoint a small subset of functionally important ...variants. To fill these unmet needs, we developed the ANNOVAR tool to annotate single nucleotide variants (SNVs) and insertions/deletions, such as examining their functional consequence on genes, inferring cytogenetic bands, reporting functional importance scores, finding variants in conserved regions, or identifying variants reported in the 1000 Genomes Project and dbSNP. ANNOVAR can utilize annotation databases from the UCSC Genome Browser or any annotation data set conforming to Generic Feature Format version 3 (GFF3). We also illustrate a 'variants reduction' protocol on 4.7 million SNVs and indels from a human genome, including two causal mutations for Miller syndrome, a rare recessive disease. Through a stepwise procedure, we excluded variants that are unlikely to be causal, and identified 20 candidate genes including the causal gene. Using a desktop computer, ANNOVAR requires ~4 min to perform gene-based annotation and ~15 min to perform variants reduction on 4.7 million variants, making it practical to handle hundreds of human genomes in a day. ANNOVAR is freely available at http://www.openbioinformatics.org/annovar/.
Opioids are very effective analgesics, but they are also highly addictive. Methadone is used to treat opioid dependence (OD), acting as a selective agonist at the μ-opioid receptor encoded by the ...gene OPRM1. Determining the optimal methadone maintenance dose is time consuming; currently, no biomarkers are available to guide treatment. In methadone-treated OD subjects drawn from a case and control sample, we conducted a genome-wide association study of usual daily methadone dose. In African-American (AA) OD subjects (n=383), we identified a genome-wide significant association between therapeutic methadone dose (mean=68.0 mg, s.d.=30.1 mg) and rs73568641 (P=2.8 × 10
), the nearest gene (306 kilobases) being OPRM1. Each minor (C) allele corresponded to an additional ~20 mg day
of oral methadone, an effect specific to AAs. In European-Americans (EAs) (n=1027), no genome-wide significant associations with methadone dose (mean=77.8 mg, s.d.=33.9 mg) were observed. In an independent set of opioid-naive AA children being treated for surgical pain, rs73568641-C was associated with a higher required dose of morphine (n=241, P=3.9 × 10
). Similarly, independent genomic loci previously shown to associate with higher opioid analgesic dose were associated with higher methadone dose in the OD sample (AA and EA: n=1410, genetic score P=1.3 × 10
). The present results in AAs indicate that genetic variants influencing opioid sensitivity across different clinical settings could contribute to precision pharmacotherapy for pain and addiction.
Deposition of amyloid-β (Aβ) in the cerebral cortex is thought to be a pivotal event in Alzheimer's disease (AD) pathogenesis with a significant genetic contribution. Molecular imaging can provide an ...early noninvasive phenotype, but small samples have prohibited genome-wide association studies (GWAS) of cortical Aβ load until now. We employed florbetapir ((18)F) positron emission tomography (PET) imaging to assess brain Aβ levels in vivo for 555 participants from the Alzheimer's Disease Neuroimaging Initiative (ADNI). More than six million common genetic variants were tested for association to quantitative global cortical Aβ load controlling for age, gender and diagnosis. Independent genome-wide significant associations were identified on chromosome 19 within APOE (apolipoprotein E) (rs429358, P=5.5 × 10(-14)) and on chromosome 3 upstream of BCHE (butyrylcholinesterase) (rs509208, P=2.7 × 10(-8)) in a region previously associated with serum BCHE activity. Together, these loci explained 15% of the variance in cortical Aβ levels in this sample (APOE 10.7%, BCHE 4.3%). Suggestive associations were identified within ITGA6, near EFNA5, EDIL3, ITGA1, PIK3R1, NFIB and ARID1B, and between NUAK1 and C12orf75. These results confirm the association of APOE with Aβ deposition and represent the largest known effect of BCHE on an AD-related phenotype. BCHE has been found in senile plaques and this new association of genetic variation at the BCHE locus with Aβ burden in humans may have implications for potential disease-modifying effects of BCHE-modulating agents in the AD spectrum.
Anorexia nervosa (AN) is a mental illness with high mortality that most commonly afflicts adolescent female individuals. Clinical symptoms include chronic food refusal, weight loss and body image ...distortions. We carried out a genome-wide association study on 1033 AN cases and 3733 pediatric control subjects, all of whom were of European ancestry and were genotyped on the Illumina HumanHap610 platform (Illumina, San Diego, CA, USA). We confirmed that common single-nucleotide polymorphisms (SNPs) within OPRD1 (rs533123, P=0.0015) confer risk for AN, and obtained suggestive evidence that common SNPs near HTR1D (rs7532266, P=0.04) confer risk for restricting-type AN specifically. However, no SNPs reached genome-wide significance in our data, whereas top association signals were detected near ZNF804B, CSRP2BP, NTNG1, AKAP6 and CDH9. In parallel, we performed genome-wide analysis on copy number variations (CNVs) using the signal intensity data from the SNP arrays. We did not find evidence that AN cases have more CNVs than control subjects, nor do they have over-representation of rare or large CNVs. However, we identified several regions with rare CNVs that were only observed in AN cases, including a recurrent 13q12 deletion (1.5 Mb) disrupting SCAS in two cases, and CNVs disrupting the CNTN6/CNTN4 region in several AN cases. In conclusion, our study suggests that both common SNPs and rare CNVs may confer genetic risk to AN. These results point to intriguing genes that await further validation in independent cohorts for confirmatory roles in AN.
Although considerable evidence suggests that the chemical synapse is a lynchpin underlying affective disorders, how molecular insults differentially affect specific synaptic connections remains ...poorly understood. For instance, Neurexin 1a and 2 (NRXN1 and NRXN2) and CNTNAP2 (also known as CASPR2), all members of the neurexin superfamily of transmembrane molecules, have been implicated in neuropsychiatric disorders. However, their loss leads to deficits that have been best characterized with regard to their effect on excitatory cells. Notably, other disease-associated genes such as BDNF and ERBB4 implicate specific interneuron synapses in psychiatric disorders. Consistent with this, cortical interneuron dysfunction has been linked to epilepsy, schizophrenia and autism. Using a microarray screen that focused upon synapse-associated molecules, we identified Cntnap4 (contactin associated protein-like 4, also known as Caspr4) as highly enriched in developing murine interneurons. In this study we show that Cntnap4 is localized presynaptically and its loss leads to a reduction in the output of cortical parvalbumin (PV)-positive GABAergic (γ-aminobutyric acid producing) basket cells. Paradoxically, the loss of Cntnap4 augments midbrain dopaminergic release in the nucleus accumbens. In Cntnap4 mutant mice, synaptic defects in these disease-relevant neuronal populations are mirrored by sensory-motor gating and grooming endophenotypes; these symptoms could be pharmacologically reversed, providing promise for therapeutic intervention in psychiatric disorders.
Autism spectrum disorders (ASDs) comprise a constellation of highly heritable neuropsychiatric disorders. Genome-wide studies of autistic individuals have implicated numerous minor risk alleles but ...few common variants, suggesting a complex genetic model with many contributing loci. To assess commonality of biological function among rare risk alleles, we compared functional knowledge of genes overlapping inherited structural variants in idiopathic ASD subjects relative to healthy controls. In this study we show that biological processes associated with synapse function and neurotransmission are significantly enriched, with replication, in ASD subjects versus controls. Analysis of phenotypes observed for mouse models of copy-variant genes established significant and replicated enrichment of observable phenotypes consistent with ASD behaviors. Most functional terms retained significance after excluding previously reported ASD loci. These results implicate several new variants that involve synaptic function and glutamatergic signaling processes as important contributors of ASD pathophysiology and suggest a sizable pool of additional potential ASD risk loci.
Numerous investigations support decreased glutamatergic signaling as a pathogenic mechanism of schizophrenia, yet the molecular underpinnings for such dysregulation are largely unknown. In the ...post-mortem dorsolateral prefrontal cortex (DLPFC), we found striking decreases in tyrosine phosphorylation of N-methyl-D aspartate (NMDA) receptor subunit 2 (GluN2) that is critical for neuroplasticity. The decreased GluN2 activity in schizophrenia may not be because of downregulation of NMDA receptors as MK-801 binding and NMDA receptor complexes in postsynaptic density (PSD) were in fact increased in schizophrenia cases. At the postreceptor level, however, we found striking reductions in the protein kinase C, Pyk 2 and Src kinase activity that in tandem can decrease GluN2 activation. Given that Src serves as a hub of various signaling mechanisms affecting GluN2 phosphorylation, we postulated that Src hypoactivity may result from convergent alterations of various schizophrenia susceptibility pathways and thus mediate their effects on NMDA receptor signaling. Indeed, the DLPFC of schizophrenia cases exhibit increased PSD-95 and erbB4 and decreased receptor-type tyrosine-protein phosphatase-α (RPTPα) and dysbindin-1, each of which reduces Src activity via protein interaction with Src. To test genomic underpinnings for Src hypoactivity, we examined genome-wide association study results, incorporating 13 394 cases and 34 676 controls. We found no significant association of individual variants of Src and its direct regulators with schizophrenia. However, a protein-protein interaction-based network centered on Src showed significant enrichment of gene-level associations with schizophrenia compared with other psychiatric illnesses. Our results together demonstrate striking decreases in NMDA receptor signaling at the postreceptor level and propose Src as a nodal point of convergent dysregulations affecting NMDA receptor pathway via protein-protein associations.
Adults with psychotic disorders have dysconnectivity in critical brain networks, including the default mode (DM) and the cingulo-opercular (CO) networks. However, it is unknown whether such deficits ...are present in youth with less severe symptoms. We conducted a multivariate connectome-wide association study examining dysconnectivity with resting state functional magnetic resonance imaging in a population-based cohort of 188 youths aged 8-22 years with psychosis-spectrum (PS) symptoms and 204 typically developing (TD) comparators. We found evidence for multi-focal dysconnectivity in PS youths, implicating the bilateral anterior cingulate, frontal pole, medial temporal lobe, opercular cortex and right orbitofrontal cortex. Follow-up seed-based and network-level analyses demonstrated that these results were driven by hyper-connectivity among DM regions and diminished connectivity among CO regions, as well as diminished coupling between frontal and DM regions. Collectively, these results provide novel evidence for functional dysconnectivity in PS youths, which show marked correspondence to abnormalities reported in adults with established psychotic disorders.
Attention-deficit/hyperactivity disorder (ADHD) is a common and highly heritable disorder, but specific genetic factors underlying risk remain elusive. To assess the role of structural variation in ...ADHD, we identified 222 inherited copy number variations (CNVs) within 335 ADHD patients and their parents that were not detected in 2026 unrelated healthy individuals. Although no excess CNVs, either deletions or duplications, were found in the ADHD cohort relative to controls, the inherited rare CNV-associated gene set was significantly enriched for genes reported as candidates in studies of autism, schizophrenia and Tourette syndrome, including A2BP1, AUTS2, CNTNAP2 and IMMP2L. The ADHD CNV gene set was also significantly enriched for genes known to be important for psychological and neurological functions, including learning, behavior, synaptic transmission and central nervous system development. Four independent deletions were located within the protein tyrosine phosphatase gene, PTPRD, recently implicated as a candidate gene for restless legs syndrome, which frequently presents with ADHD. A deletion within the glutamate receptor gene, GRM5, was found in an affected parent and all three affected offspring whose ADHD phenotypes closely resembled those of the GRM5 null mouse. Together, these results suggest that rare inherited structural variations play an important role in ADHD development and indicate a set of putative candidate genes for further study in the etiology of ADHD.
The amygdala brain region has been implicated in the pathophysiology of schizophrenia through emotion processing. However, transcriptome messages in the amygdala of schizophrenia patients have not ...been well studied. We used RNA sequencing to investigate gene-expression profiling in the amygdala tissues, and identified 569 upregulated and 192 downregulated genes from 22 schizophrenia patients and 24 non-psychiatric controls. Gene functional enrichment analysis demonstrated that the downregulated genes were enriched in pathways such as 'synaptic transmission' and 'behavior', whereas the upregulated genes were significantly over-represented in gene ontology pathways such as 'immune response' and 'blood vessel development'. Co-expression-based gene network analysis identified seven modules including four modules significantly associated with 'synaptic transmission', 'blood vessel development' or 'immune responses'. Taken together, our study provides novel insights into the molecular mechanism of schizophrenia, suggesting that precision-tailored therapeutic approaches aimed at normalizing the expression/function of specific gene networks could be a promising option in schizophrenia.