Autworks: a cross-disease network biology application for Autism and related disorders
BMC Medical Genomics volume 5, Article number: 56 (2012)
The genetic etiology of autism is heterogeneous. Multiple disorders share genotypic and phenotypic traits with autism. Network based cross-disorder analysis can aid in the understanding and characterization of the molecular pathology of autism, but there are few tools that enable us to conduct cross-disorder analysis and to visualize the results.
We have designed Autworks as a web portal to bring together gene interaction and gene-disease association data on autism to enable network construction, visualization, network comparisons with numerous other related neurological conditions and disorders. Users may examine the structure of gene interactions within a set of disorder-associated genes, compare networks of disorder/disease genes with those of other disorders/diseases, and upload their own sets for comparative analysis.
Autworks is a web application that provides an easy-to-use resource for researchers of varied backgrounds to analyze the autism gene network structure within and between disorders.
Autism spectrum disorder (ASD) is one of the most common [1, 2] and highly heritable [3, 4] neurodevelopmental disorders. Despite, and perhaps because of, its commonality, ASD is known to have a highly complex genetic etiology [5–8].
More than 100 genes have been reported to be perturbed in individuals with ASD , while only a few such genes have been verified by genome-wide association studies or gene expression analyses. There are a number of Mendelian disorders in which patients may manifest ASD traits, including fragile X syndrome , tuberous sclerosis complex , Rett syndrome , and Angelman syndrome  to name a few. Finally, common copy number variations have been found between ASD and other complex mental disorders, including schizophrenia (e.g., 1q21.1 duplication), bipolar disorder (e.g., 16p11.2 duplication), and epilepsy (e.g., 17q12 deletion) [9, 14].
The complexity of the autism genetic landscape suggests that networks of candidate genes will be of value to identify biological themes that contribute to the spectrum of autism phenotypes. Innovative work by Rzhetsky et al.  demonstrated that significant correlations of comorbidity were found across a wide range of disorders, ASD in particular. Related research demonstrated that combining the known genetic etiology of behaviorally related disorders with ASD through biological networks could provide particular insight into genes with great significance in ASD, as well as predictive power for novel autism-associated variants .
To capitalize on the value of network biology and a cross-disease understanding of autism, we built a web application called Autworks to allow for network scale cross-disorder analysis. We aimed to integrate pathogenic and genomic information spread across different data repositories and to augment this with novel features including network analysis and visualization tools lacking in existing applications. To this end, Autworks incorporates disease-specific, manually curated databases that provide gene associations for a few major mental disorders (e.g., [17–20]), pathogenic information from general-purpose genetic databases, including GeneCards , HuGE Navigator , and PharnGKB , as well as our own PubMed-based gene-disorder search results. Finally, we examined gene/protein interaction databases [24–28] and incorporated their functionality into Autworks, as these repositories mainly focus on general cellular interactions without disease-context.
Construction and content
Genes and disorders
Genes constitute nodes in our network analysis. 19,191 protein-coding human genes and their descriptive data are derived from the official list of human genes provided by the HUGO Gene Nomenclature Committee (HGNC) . Diseases are derived from the 2012 version of Medical Subject Headings (MeSH). Specifically, disorder names and their entry terms were gathered for all terms in sub-tree categories from C04 (Neoplasms) to C20 (Immune System Diseases) as well as F03 (Mental Disorders), in order to capture disorders with a possible genetic etiology.
For each disorder, gene-given disorder associations are derived from a PubMed search linking the given disorder name with gene names. Disorder names were derived as described above. Each MeSH term as well as all associated entry terms was queried for gene associations--the results were combined for a list of all gene-disorder associations. This produces 2,711 disorders with 660,090 associations with 16,685 genes.
Gene interaction network data is based on the Search Tool for Retrieval of Interacting Genes/Proteins (STRING) . In Autworks, the evidence for a gene interaction is broken down along five different lines of evidence, including protein-protein interactions, pathway information, and co-expressed arrays. Networks can be constructed using specific lines of evidence, omitting evidence not considered useful for a particular inquiry. By default, interactions with a medium score as defined by STRING (400 or greater) constitute valid interactions. Ensembl  was used to map proteins from STRING to genes in Autworks. This produces 367,308 interactions between genes in Autworks.
Autworks contains a measure of association between disorders based on the hypergeometric distribution. By considering the common genes shared by two disorders, this method tests the likelihood of seeing this amount of overlap (or more) by random chance. Enrichment values were calculated using the GNU Scientific Library. For each pair of disease gene sets [D1;D2] and the set of all human protein-coding genes associated with any disorder [G] we calculated the disease enrichment score as a p-value by applying the function gsl_cdf_hypergeometric Q(k; n1; n2; t) with the following parameters: k = |D1 ∩ D2| - 1 , n1 = |D2|, n2 = |G|, t = |D1|. Sets containing fewer than 5 genes were not considered for this analysis.
For each disease Autworks provides four distinct views: a cross-disorder view that visualizes a network of disorders against which the current disorder has significant genetic overlap, a list view of genes in the set ordered by significance, a gene network view that interactively shows the relationship of genes within a given disorder, as shown in Figure 1, and a genetic variant view listing known variants associated with the disease via genome wide association studies.
The gene network visualization component of Autworks contains several features that allow researchers to use only the genes and interactions they consider interesting. Researchers can choose to display interactions based on pathway databases, experimentally proven interactions, co-expression, or some combination of the above. Details of the data behind each interaction are linked to from the network visualization. In addition to the evidence behind the interaction, this view shows any diseases the two genes have in common, as well as places in the literature in which they occur together. Autworks also allows researchers to highlight genes that belong to other sets, are implicated in a large number of other disorders. For each interaction, specific information about the type of interaction as well as papers that mention both genes together are available.
As an demonstration of a single gene investigation in Autworks, here we show a Mendelian syndrome which shares phenotypic traits with autism spectrum disorder. Rubinstein-Taybi syndrome can be caused by a mutation in the gene which codes for CREB binding protein (CREBBP ). The main symptoms of this disorder include slow development of cognitive and motor skills, and researchers have speculated about classifying the syndrome among the autisms . A search for this gene in Autworks reveals that this gene occupies a place in the largest network of genes implicated in autism, mediating inter-actions between genes responsible for the calcium channel receptors (CACNA1C, CACNA1G, CACNA1H), mutations in which are responsible for Timothy syndrome, another syndromic variant of autism , and a larger network of genes associated with autism through genome wide association studies (GWAS) and gene expression experiments (AR, BCL2, GRIN2A, NFKB1, among others). These interactions are based on annotated gene pathways, including the MAPK signaling pathways  and the JAK-STAT signaling pathways . The position of this gene in the autism interaction network and the association of this gene with a disorder that has an autism-like phenotype make it an interesting candidate for further inquiry. This gene interaction based evidence for the association of CREBBP with autism genes is not known to have been reported in the literature before. Autworks makes this type of inquiry remarkably easy to perform.
Researchers may also approach the site with interest in a set of genes, possibly ones derived from one's own experiments, experiments reported in the literature, or biological processes. As a particular example of this, researchers have expressed interest in disruption of the glutamatergic synapse as a possible causative factor of autism. The pathway is well described in the Kyoto Encyclopedia of Genes and Genomes (KEGG), with a list of genes and interactions involved in the pathway. Importing the set of genes from this pathway into Autworks provides an enrichment analysis against the set of disorders in Autworks. As this set describes a neurologic function, it is significantly enriched for many different neurological disorders, including autism. To get a sense of how this relates to the network of interactions in autism, one can mark the genes occurring in this set on the autism network using Autworks as shown in Figure 1. One can clearly see a significant number of the genes in this pathway in the largest network of interacting autism genes, suggesting that this pathway may indeed play a significant role in autism.
Autworks provides a dynamic picture of the network of autism candidate genes and a set of network biology methods for researchers of varied backgrounds. These methods should assist in the analysis of the autism network and its relationship with the networks of other, related human disorders and in advancements towards a clearer understanding of the genetic landscape of autism.
Availability and requirements
Project name: autworks
Project home page: http://autworks.hms.harvard.edu
Operating system(s): platform independent
Programming Language: Ruby
Other requirements: HTML5 compliant web browser
Jon B, et al: Prevalence of autism spectrum disorders – autism and developmental disabilities monitoring network, 14 sites, United States, 2008. MMWR Morb Mortal Wkly Rep. 2012, 61 (SS03): 1-19. http://www.cdc.gov/mmwr/preview/mmwrhtml/ss6103a1.htm?s_cid=ss6103a1_w.
Fombonne E: Epidemiology of pervasive developmental disorders. Pediatr Res. 2009, 65 (6): 591-598. 10.1203/PDR.0b013e31819e7203.
Hallmayer J, Cleveland S, Torres A, Phillips J, Cohen B, Torigoe T, Miller J, Fedele A, Collins J, Smith K, Lotspeich L, Croen LA, Ozonoff S, Lajonchere C, Grether JK, Risch N: Genetic heritability and shared environmental factors among twin pairs with autism. Arch Gen Psychiatry. 2011, 68 (11): 1095-1102. 10.1001/archgenpsychiatry.2011.76.
Bailey A, Le Couteur A, Gottesman I, Bolton P, Simonoff E, Yuzda E, Rutter M: Autism as a strongly genetic disorder: evidence from a British twin study. Psychol Med. 1995, 25 (1): 63-77. 10.1017/S0033291700028099.
Sullivan PF, Daly MJ, O’Donovan M: Genetic architectures of psychiatric disorders: the emerging picture and its implications. Nat Rev Genet. 2012, 13 (8): 537-551. 10.1038/nrg3240.
Mefford HC, Batshaw ML, Hoffman EP: Genomics, intellectual disability, and autism. N Engl J Med. 2012, 366 (8): 733-743. 10.1056/NEJMra1114194.
Geschwind DH: Genetics of autism spectrum disorders. Trends Cogn Sci. 2011, 15 (9): 409-416. 10.1016/j.tics.2011.07.003.
State MW, Levitt P: The conundrums of understanding genetic risks for autism spectrum disorders. Nat Neurosci. 2011, 14 (12): 1499-1506. 10.1038/nn.2924.
Betancur C: Etiological heterogeneity in autism spectrum disorders: more than 100 genetic and genomic disorders and still counting. Brain Res. 2011, 1380: 42-77.
Wang LW, Berry-Kravis E, Hagerman RJ: Fragile X: leading the way for targeted treatments in autism. Neurotherapeutics. 2010, 7 (3): 264-274. 10.1016/j.nurt.2010.05.005.
Muzykewicz DA, Newberry P, Danforth N, Halpern EF, Thiele EA: Psychiatric comorbid conditions in a clinic population of 241 patients with tuberous sclerosis complex. Epilepsy Behav. 2007, 11 (4): 506-513. 10.1016/j.yebeh.2007.07.010.
Percy AK: Rett syndrome: exploring the autism link. Arch Neurol. 2011, 68 (8): 985-989. 10.1001/archneurol.2011.149.
Veltman MW, Craig EE, Bolton PF: Autism spectrum disorders in Prader-Willi and Angelman syndromes: a systematic review. Psychiatr Genet. 2005, 15 (4): 243-254. 10.1097/00041444-200512000-00006.
Malhotra D, Sebat J: CNVs: harbingers of a rare variant revolution in psychiatric genetics. Cell. 2012, 148 (6): 1223-1241. 10.1016/j.cell.2012.02.039.
Rzhetsky A, Wajngurt D, Park N, Zheng T: Probing genetic overlap among complex human phenotypes. Proc Natl Acad Sci U S A. 2007, 104: 11694-11699. 10.1073/pnas.0704820104.
Wall DP, et al: Comparative analysis of neurological disorders focuses genome-wide search for autism genes. Genomics. 2009, 93 (2): 120-129. 10.1016/j.ygeno.2008.09.015.
Xu LM, Li JR, Huang Y, Zhao M, Tang X, Wei L: AutismKB: an evidence-based knowledgebase of autism genetics. Nucleic Acids Res. 2012, 40 (Database issue): D1016-D1022. http://autismkb.cbi.pku.edu.cn.
Lill CM, Roehr JT, McQueen MB, Kawoura FK, Bagade S, et al: Comprehensive research synopsis and systematic meta-analyses in Parkinson’s disease genetics: the PDGene database. PLoS Genet. 2012, 8 (3): e1002548-10.1371/journal.pgen.1002548. http://www.pdgene.org.
Banerjee-Basu S, Packer A: SFARI gene: an evolving database for the autism research community. Dis Model Mech. 2010, 3 (3–4): 133-135. http://gene.sfari.org.
Allen NC, Bagade S, McQueen MB, Ioannidis JP, Kawoura FK, Khoury MJ, Tanzi RE, Bertram L: Systematic meta-analysis and field synopsis of genetic association studies in schizophrenia: the SzGene database. Nat Genet. 2008, 40 (7): 827-834. 10.1038/ng.171. http://www.szgene.org.
Stelzer G, Dalah I, Iny Stein T, Satanower Y, Rosen N, Nativ N, Oz-Levi D, Olender T, Belinky F, Bahir I, Krug H, Perco P, Mayer B, Kolker E, Safran M, Lancet D: In-silico human genomics with GeneCards. Human Genomics. 2011, 5 (6): 709-717. 10.1186/1479-7364-5-6-709. http://www.genecards.org.
Yu W, Gwinn M, Clyne M, Yesupriya A, Khoury MJ: A navigator for human genome epidemiology. Nat Genet. 2008, 40 (2): 124-125. 10.1038/ng0208-124. http://www.hugenavigator.net.
Hewett M, Oliver DE, Rubin DL, Easton KL, Stuart JM, Altman RB, Klein TE, Pharm GKB: The pharmacogenetics knowledge base. Nucleic Acids Res. 2002, 30 (1): 163-165. 10.1093/nar/30.1.163. http://www.pharmgkb.org.
Dahlquist KD, Salomonis N, Vranizan K, Lawlor SC, Conklin BR: GenMAPP, a new tool for viewing and analyzing data on biological pathways. Nat Genet. 2002, 31 (1): 19-20. 10.1038/ng0502-19. http://www.genmapp.org.
Mostafavi S, Ray D, Warde-Farley D, Grouios C, Morris Q: GeneMANIA: a real-time multiple association network integration algorithm for predicting gene function. Genome Biol. 2008, 9 (Suppl 1): S4-10.1186/gb-2008-9-s1-s4. http://www.genemania.org.
Szklarczyk D, Franceschini A, Kuhn M, Simonovic M, Roth A, Minguez P, Doerks T, Stark M, Muller J, Bork P, Jensen LJ, von Mering C: The STRING database in 2011: functional interaction networks of proteins, globally integrated and scored. Nucleic Acids Res. 2011, 39 (Database issue): D561-D568. http://string.embl.de.
Huttenhower C, Haley EM, Hibbs MA, Dumeaux V, Barrett DR, Coller HA: Troyanskaya OG: Exploring the human genome with functional maps. Genome Res. 2009, 19 (6): 1093-1106. 10.1101/gr.082214.108. http://hefalmap.princeton.edu.
Keshava Prasad TS, Goel R, Kandasamy K, Keerthikumar S, Kumar S, et al: Human Protein Reference Database -- 2009 update. Nucleic Acids Res. 2009, 37 (Database issue): D767-D772. http://www.hprd.org.
Seal RL, Gordon SM, Lush MJ, Wright MW, Bruford EA: genenames.org: the HGNC resources in 2011. Nucleic Acids Res. 2011, 39 (Database issue): D514-D519. http://genenames.org.
Flicek P, Amode MR, Barrell D, Beal K, Brent S, et al: Ensembl 2012. Nucleic Acids Res. 2011, 40 (Database issue): D84-D90. http://www.ensembl.org.
Stef M, et al: Spectrum of CREBBP gene dosage anomalies in Rubinstein-Taybi syndrome patients. Eur J Hum Genet. 2007, 15 (8): 843-847. 10.1038/sj.ejhg.5201847.
Galéra C, et al: Socio-behavioral characteristics of children with Rubeinstein-Taybi syndrome. J Autism Dev Discord. 2009, 39 (9): 1252-1260. 10.1007/s10803-009-0733-4.
Bader PL, et al: Mouse model of Timothy syndrome recapitulates triad of autistic traits. Proc Natl Acad Sci USA. 2011, 108 (37): 15432-15437. 10.1073/pnas.1112667108.
Cargnello M, Roux PP: Activation and function of the MAPKs and their substrates, the MAPK-activated protein kinases. Microbiol Mol Biol Rev. 2011, 75 (1): 50-83. 10.1128/MMBR.00031-10.
Vera J, Rateitschak K, Lange F, Kossow C, Wolkenhauer O, Jaster R: Systems biology of JAK-STAT signaling in human malignancies. Prog Biophys Mol Biol. 2011, 106 (2): 426-434. 10.1016/j.pbiomolbio.2011.06.013.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1755-8794/5/56/prepub
Funding: This work was supported by the National Science Foundation [0543480 to D.P.W, 0640809 to D.P.W]; and the National Institutes of Health [LM009261 to D.P.W].
The authors declare they have no competing interests.
DPW conceived the study and participated in application evaluation and design. THN and JYJ developed the application and wrote the manuscript and participated in application evaluation and design. TFD edited the manuscript and participated in application evaluation and design. All authors read and approved the final document.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Nelson, T.H., Jung, JY., DeLuca, T.F. et al. Autworks: a cross-disease network biology application for Autism and related disorders. BMC Med Genomics 5, 56 (2012). https://doi.org/10.1186/1755-8794-5-56