Skip to main content
Advertisement
  • Other Publications
    • EMBO Press
    • Molecular Systems Biology (Home)
    • The EMBO Journal
    • EMBO reports
    • EMBO Molecular Medicine
    • Life Science Alliance
Login

   

Search

Advanced Search

Journal

  • Home
  • Current Issue
  • Archive
  • Subject Collections
  • Review Series & Focuses
  • Podcasts & Videos

Authors & Referees

  • Submit
  • Author Guidelines
  • Aims & Scope
  • Editors & Board
  • Transparent Process
  • Bibliometrics
  • Open Access
  • Referee Guidelines

Info

  • E-Mail Editorial Office
  • Alerts
  • RSS Feeds
  • Reprints & Permissions
  • Advertise & Sponsor
  • Media Partners
  • News & Press
  • Customer Service
  • Home
  • Molecular Systems Biology: 10 (12)

Open Access

Transparent Process

Article

Integrated systems analysis reveals a molecular network underlying autism spectrum disorders

Jingjing Li, Minyi Shi, Zhihai Ma, Shuchun Zhao, Ghia Euskirchen, Jennifer Ziskin, Alexander Urban, Joachim Hallmayer, Michael Snyder
DOI 10.15252/msb.20145487 | Published online 30.12.2014
Molecular Systems Biology (2014) 10, 774
Jingjing Li
Department of Genetics, Stanford Center for Genomics and Personalized Medicine Stanford University School of Medicine, Stanford, CA, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Minyi Shi
Department of Genetics, Stanford Center for Genomics and Personalized Medicine Stanford University School of Medicine, Stanford, CA, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Zhihai Ma
Department of Genetics, Stanford Center for Genomics and Personalized Medicine Stanford University School of Medicine, Stanford, CA, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Shuchun Zhao
Department of Pathology, Stanford University School of Medicine, Stanford, CA, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Ghia Euskirchen
Department of Genetics, Stanford Center for Genomics and Personalized Medicine Stanford University School of Medicine, Stanford, CA, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Jennifer Ziskin
Department of Pathology, Stanford University School of Medicine, Stanford, CA, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Alexander Urban
Department of Psychiatry & Behavioral Sciences, Stanford University School of Medicine, Stanford, CA, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Joachim Hallmayer
Department of Psychiatry & Behavioral Sciences, Stanford University School of Medicine, Stanford, CA, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Michael Snyder
Department of Genetics, Stanford Center for Genomics and Personalized Medicine Stanford University School of Medicine, Stanford, CA, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site

Author Affiliations

  1. Jingjing Li1,†,
  2. Minyi Shi1,†,
  3. Zhihai Ma1,†,
  4. Shuchun Zhao2,
  5. Ghia Euskirchen1,
  6. Jennifer Ziskin2,
  7. Alexander Urban3,
  8. Joachim Hallmayer3 and
  9. Michael Snyder*,1
  1. 1Department of Genetics, Stanford Center for Genomics and Personalized Medicine Stanford University School of Medicine, Stanford, CA, USA
  2. 2Department of Pathology, Stanford University School of Medicine, Stanford, CA, USA
  3. 3Department of Psychiatry & Behavioral Sciences, Stanford University School of Medicine, Stanford, CA, USA
  1. ↵*Corresponding author. Tel: +1 650 736 8099; E‐mail: mpsnyder{at}stanford.edu
  1. ↵† These authors contributed equally to this work

View Abstract
  • Article
  • Figures & Data
  • Transparent Process
Loading

Article Figures & Data

Figures

  • Supplementary Materials
  • Figure 1.
    • Download figure
    • Open in new tab
    • Download powerpoint
    Figure 1. A modular protein interaction network with modules containing enrichment of autism‐associated genes

    Two topological modules (#2 and #13) on human protein interaction network showed significant enrichment for autism genes (in red). The topological modules are physical clusters on the network where their member genes intensively interact with each other but sparsely interact with non‐member genes on the network. A zoom‐in view of module #13 is also shown, where known autism genes and genes affected by ASD‐associated de novo CNVs were colorized in red and green, respectively. Genes annotated by both were in blue. The false discovery rate indicates its significant enrichment for the known autism genes.

  • Figure 2.
    • Download figure
    • Open in new tab
    • Download powerpoint
    Figure 2. Candidate genes from sequencing screens

    1. An overview of the identified loci from whole‐genome and exome sequencing. Evolutionary conservation is quantified by GERP++ score, where the higher scores indicate greater selective pressure on the genomic loci. For genes with multiple significant loci, the most conserved residue is considered. Variants absent in the 1,000 Genome dataset are considered rare variants. The genes were colorized based on the fraction of deleterious mutations predicted by MutationTaster among all the identified mutations in the gene (MRI image of the corpus callosum: Allen Institute of Brain Science).

    2. Validation using another larger patient cohort. In this dataset, variants with allele frequencies with increased absolute differences between cases and controls are more likely to affect genes that were also detected in our study (red line). The allele frequency difference is the absolute difference between cases and controls. This trend cannot be observed by 10,000 simulations (blue line for one randomized dataset).

  • Figure 3.
    • Download figure
    • Open in new tab
    • Download powerpoint
    Figure 3. Expression analysis of the synaptic module

    1. Dichotomized expression of the genes in module #13 across 295 brain sections. Relative abundance of each gene across the 295 brain sections was hierarchically clustered to reveal gene groups exhibiting similar expression patterns across tissues. Group 1 genes showed elevated expression in 175 regions (T1, e.g., corpus callosum) relative to other brain sections, and Group 2 genes showed high expression in 120 regions (T2, e.g., hippocampal regions) relative to other brain sections.

    2. RNA‐sequencing of four different brain regions from a healthy subject. The brain regions include the Brodmann areas 9 (BA9), 40 (BA40), the amygdala (AMY) and the corpus callosum (CC), which revealed the same observation as from the microarray analyses. Group 1 (red) and 2 (blue) genes were compared with 1,000 randomly sampled genes (gray) from the transcriptome in each brain region. The raw FPKM values were normalized into the cumulative density functions based on kernel density estimation. The elevation of Group 2 genes across all brain regions and the greatest increase of Group 1 genes in the corpus callosum were all statistically significant (P < 1e‐5, Wilcoxon rank‐sum test).

    3. RNA‐sequencing of the corpus callosum transcriptomes from six non‐autistic individuals. FPKM quantifies the absolute expression of genes in each group. The two groups have similar expression in the corpus callosum (P > 0.5, Wilcoxon rank‐sum test), which are all above the transcriptome background (P < 4.87e‐6, Wilcoxon rank‐sum test), suggesting that both sub‐components are active in this tissue.

  • Figure 4.
    • Download figure
    • Open in new tab
    • Download powerpoint
    Figure 4. Cell‐type expression of module #13 in oligodendrocytes

    1. Immunohistochemistry analysis in the corpus callosum. Staining of LRP2 in the human corpus callosum reveals that the major cell population in the corpus callosum is the oligodendrocytes (the blue round nuclei), which express LRP2 stained in brown. A zoom‐in view is shown in the inset.

    2. Neural cell‐type expression of the orthologous module #13 in mouse brain. Gene expression in different neural cell types was hierarchically clustered into the three major cell types in brain (neurons, oligodendrocytes and astrocytes). The clustering grouped genes in module #13 into a neuron cluster and a glial cluster, enriched for Group 1 and 2 genes, respectively. The fraction of Group 1 (red) and 2 (blue) genes in the glial and neuronal clusters were represented by the pie charts, with statistical significance determined by a chi‐square test.

    3. Overall expression of module #13 in cultured oligodendrocyte precursor cells (OPCs). Group 1 and 2 were expressed at a similar level as the transcriptome background in OPCs. The statistical significance was determined by Wilcoxon rank‐sum test, and the error bars represent one standard error.

    4. The role of the module in oligodendrocyte (OL) development. Differentiation of OPCs into mature myelinating OLs (MOG+) led to a significant up‐regulation of Group 1 genes (left, OPCs → mature OLs). On the other hand, conditional knockout (CKO) of the master myelination factor MRF from mature OLs led to a significant up‐regulation of Group 2 genes (right, mature OLs → MRF CKO). The statistical significance was determined by Wilcoxon rank‐sum test.

    5. A proposed model. Up‐regulation is associated with, or likely to contribute to, the differentiation of OPCs into mature myelinating OLs. The mature OLs acquire their myelination capacity by activating the MRF‐mediated regulatory network, which also serves to repress expression of Group 2 genes.

  • Figure 5.
    • Download figure
    • Open in new tab
    • Download powerpoint
    Figure 5. Integrative analysis of the genetic alteration in this study

    1. Enrichment of the differentially expressed genes in module #13. RNA‐sequencing of the corpus callosum of autism patients and their matched controls. Enrichment was not observed for the genes in the human synaptome or the collection of known autism genes (excluding genes in this module). Statistical significance was determined by hypergeometric test.

    2. The mutation pattern of the genes from the innermost layers of the interaction network (K ≥ 10) to the periphery layer (K = 1). Genes in the central and periphery layers in this module are more likely to be affected, while the trend cannot be observed in 10,000 random simulations. For individual bins, significant enrichment and depletion were observed in the central layers (K ≥ 10) and the intermediate layers (3 ≤ K < 6), respectively. Statistical significance of the enrichment was determined by hypergeometric test. 10,000 random permutations were performed to determine the statistical significance of the curve.

    3. Compositional bias of the mutated genes in central layers. The mutated genes in central layers are more biased toward the corpus callosum‐specific subcomponent; this trend is not observed in background or other mutated genes with varying degree of K. Statistical significance of the enrichment was determined by hypergeometric test.

    4. Positive correlation between network coreness and gene expression in the corpus callosum. RNA‐sequencing of the corpus callosum of six non‐autistic individuals revealed a positive correlation, suggesting the central layers may play critical roles in the corpus callosum. Two outlier genes, DYNLL1 and BCAS1, are separately labeled due to their extreme expression in this tissue. The correlation coefficient r and its statistical significance were computed using Spearman's correlation.

    5. Predicted sub‐complexes within this module. Genes in this module are topologically clustered to form sub‐complexes, with the significantly mutated genes labeled in blue, mis‐expressed genes in the corpus callosum labeled in green, and both in red. Two clusters, #6 for SHANK‐DLGAP complexes and #6 for LRP2, and its binding partners, are enriched for the mis‐expressed or mutated genes, respectively. Statistical significance of the enrichment was determined by hypergeometric test.

Supplementary Materials

  • Figures
  • Supplementary Figure S1 [msb145487-sup-0001-FigS1.pdf]

  • Supplementary Figure S2 [msb145487-sup-0002-FigS2.pdf]

  • Supplementary Figure S3 [msb145487-sup-0003-FigS3.pdf]

  • Supplementary Figure S4 [msb145487-sup-0004-FigS4.pdf]

  • Supplementary Figure S5 [msb145487-sup-0005-FigS5.pdf]

  • Supplementary Figure S6 [msb145487-sup-0006-FigS6.pdf]

  • Supplementary Figure S7 [msb145487-sup-0007-FigS7.pdf]

  • Supplementary Figure S8 [msb145487-sup-0008-FigS8.pdf]

  • Supplementary Figure S9 [msb145487-sup-0009-FigS9.pdf]

  • Supplementary Figure S10 [msb145487-sup-0010-FigS10.pdf]

  • Supplementary Figure S11 [msb145487-sup-0011-FigS11.pdf]

  • Supplementary Figure S12 [msb145487-sup-0012-FigS12.pdf]

  • Supplementary Figure S13 [msb145487-sup-0013-FigS13.pdf]

  • Supplementary Figure S14 [msb145487-sup-0014-FigS14.pdf]

  • Supplementary Table S1 [msb145487-sup-0015-TableS1.pdf]

  • Supplementary Table S2 [msb145487-sup-0016-TableS2.pdf]

  • Supplementary Table S3 [msb145487-sup-0017-TableS3.pdf]

  • Supplementary Table S4 [msb145487-sup-0018-TableS4.pdf]

  • Supplementary Table S5 [msb145487-sup-0019-TableS5.pdf]

  • Supplementary Table S6 [msb145487-sup-0020-TableS6.pdf]

  • Supplementary Legends [msb145487-sup-0021-Legends.docx]

  • Supplementary Dataset S1 [msb145487-sup-0022-DatasetS1.xlsx]

  • Supplementary Dataset S2 [msb145487-sup-0023-DatasetS2.xlsx]

  • Supplementary Dataset S3 [msb145487-sup-0024-DatasetS3.xlsx]

  • Supplementary Dataset S4 [msb145487-sup-0025-DatasetS4.xlsx]

Previous Article in this IssueNext Article in this Issue
Back to top

  • PDF
  • Share
  • Export
  • Print
Loading

PDF

Review Process

In this Issue
Volume 10, Issue 12
01 December 2014 | pp -
Molecular Systems Biology: 10 (12)
Alert me when this article is cited
Alert me if a correction is posted

Article

  • Article
    • Abstract
    • Synopsis
    • Introduction
    • Results
    • Discussion
    • Materials and Methods
    • Acknowledgements
    • Author contributions
    • Conflict of interest
    • Supplementary Information
    • References
  • Figures & Data
  • Transparent Process

News & Press

Molecular network identified underlying autism spectrum disorders
EMBO

Request Permissions

Subject Areas

  • Chromatin, Epigenetics, Genomics & Functional Genomics
  • Genome-Scale & Integrative Biology
  • Network Biology

Journal

  • Home
  • Latest Content
  • Archive
  • Bibliometrics
  • E-Mail Editorial Office

Authors & References

  • Aims & Scope
  • Editors & Board
  • Transparent Process
  • Author Guidelines
  • Referee Guidelines
  • Open Access
  • Submit

Info

  • Alerts
  • RSS Feeds
  • Reprints & Permissions
  • Advertise & Sponsor
  • News & Press
  • Customer Service

EMBO

  • Funding & Awards
  • Events
  • Science Policy
  • Members
  • About EMBO

Online ISSN  1744-4292

Copyright© 2018 EMBO

This website is best viewed using the latest versions of all modern web browsers. Older browsers may not display correctly.