Heritability enrichment in context-specific regulatory networks improves phenotype-relevant tissue identification

  1. Zhanying Feng
  2. Zhana Duren
  3. Jingxue Xin
  4. Qiuyue Yuan
  5. Yaoxi He
  6. Bing Su
  7. Wing Hung Wong  Is a corresponding author
  8. Yong Wang  Is a corresponding author
  1. Chinese Academy of Sciences, China
  2. Clemson University, United States
  3. Stanford University, United States

Abstract

Systems genetics holds the promise to decipher complex traits by interpreting their associated SNPs through gene regulatory networks derived from comprehensive multi-omics data of cell types, tissues, and organs. Here, we propose SpecVar to integrate paired chromatin accessibility and gene expression data into context-specific regulatory network atlas and regulatory categories, conduct heritability enrichment analysis with GWAS summary statistics, identify relevant tissues, and depict common genetic factors acting in the shared regulatory networks between traits by relevance correlation. Our method improves power upon existing approaches by associating SNPs with context-specific regulatory elements to assess heritability enrichments and by explicitly prioritizing gene regulations underlying relevant tissues. Ablation studies, independent data validation, and comparison experiments with existing methods on GWAS of six phenotypes show that SpecVar can improve heritability enrichment, accurately detect relevant tissues, and reveal causal regulations. Furthermore, SpecVar correlates the relevance patterns for pairs of phenotypes and better reveals shared SNP associated regulations of phenotypes than existing methods. Studying GWAS of 206 phenotypes in UK-Biobank demonstrates that SpecVar leverages the context-specific regulatory network atlas to prioritize phenotypes' relevant tissues and shared heritability for biological and therapeutic insights. SpecVar provides a powerful way to interpret SNPs via context-specific regulatory networks and is available at https://github.com/AMSSwanglab/SpecVar.

Data availability

Codes and regulatory network resources are available at https://github.com/AMSSwanglab/SpecVar. Expression and chromatin accessibility data were summarized in Table S1. GWAS data used: GWAS summary statistics of LDL and TC were downloaded at http://csg.sph.umich.edu/willer/public/lipids2013/; GWAS summary statistics of EA (GCST006442), CP (GCST006572), BrainShape (GCST90012880-GCST90013164), and Face (GCST009464) were downloaded at GWAS catalog https://www.ebi.ac.uk/gwas/summary-statistics; GWAS summary statistics of UK-Biobank were downloaded at http://www.nealelab.is/uk-biobank. The LDSC genetic correlation and phenotypic correlation computed from individual phenotypic data were downloaded at https://ukbb-rg.hail.is/.

Article and author information

Author details

  1. Zhanying Feng

    Academy of Mathematics and Systems Science, Chinese Academy of Sciences, Beijing, China
    Competing interests
    The authors declare that no competing interests exist.
  2. Zhana Duren

    Department of Genetics and Biochemistry, Clemson University, Greenwood, United States
    Competing interests
    The authors declare that no competing interests exist.
  3. Jingxue Xin

    Department of Statistics, Stanford University, Palo Alto, United States
    Competing interests
    The authors declare that no competing interests exist.
  4. Qiuyue Yuan

    Department of Genetics and Biochemistry, Clemson University, Greenwood, United States
    Competing interests
    The authors declare that no competing interests exist.
  5. Yaoxi He

    Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming, China
    Competing interests
    The authors declare that no competing interests exist.
  6. Bing Su

    Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming, China
    Competing interests
    The authors declare that no competing interests exist.
  7. Wing Hung Wong

    Department of Statistics, Stanford University, Palo Alto, United States
    For correspondence
    whwong@stanford.edu
    Competing interests
    The authors declare that no competing interests exist.
  8. Yong Wang

    Academy of Mathematics and Systems Science, Chinese Academy of Sciences, Beijing, China
    For correspondence
    zyfeng@amss.ac.cn
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-0695-5273

Funding

National Key Research and Development Program of China (2022YFA1004800,2020YFA0712402)

  • Yong Wang

Strategic Priority Research Program of the Chinese Academy of Science (XDPB17)

  • Yong Wang

CAS Young Scientists in Basic esearch (YSBR-077)

  • Yong Wang

National Natural Science Foundation of China (12025107,11871463,11688101)

  • Yong Wang

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Copyright

© 2022, Feng et al.

This article is distributed under the terms of the Creative Commons Attribution License permitting unrestricted use and redistribution provided that the original author and source are credited.

Metrics

  • 1,199
    views
  • 128
    downloads
  • 2
    citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Zhanying Feng
  2. Zhana Duren
  3. Jingxue Xin
  4. Qiuyue Yuan
  5. Yaoxi He
  6. Bing Su
  7. Wing Hung Wong
  8. Yong Wang
(2022)
Heritability enrichment in context-specific regulatory networks improves phenotype-relevant tissue identification
eLife 11:e82535.
https://doi.org/10.7554/eLife.82535

Share this article

https://doi.org/10.7554/eLife.82535

Further reading

    1. Computational and Systems Biology
    Masaaki Uematsu, Jeremy M Baskin
    Tools and Resources

    Plasmid construction is central to life science research, and sequence verification is arguably its costliest step. Long-read sequencing has emerged as a competitor to Sanger sequencing, with the principal benefit that whole plasmids can be sequenced in a single run. Nevertheless, the current cost of nanopore sequencing is still prohibitive for routine sequencing during plasmid construction. We develop a computational approach termed Simple Algorithm for Very Efficient Multiplexing of Oxford Nanopore Experiments for You (SAVEMONEY) that guides researchers to mix multiple plasmids and subsequently computationally de-mixes the resultant sequences. SAVEMONEY defines optimal mixtures in a pre-survey step, and following sequencing, executes a post-analysis workflow involving sequence classification, alignment, and consensus determination. By using Bayesian analysis with prior probability of expected plasmid construction error rate, high-confidence sequences can be obtained for each plasmid in the mixture. Plasmids differing by as little as two bases can be mixed as a single sample for nanopore sequencing, and routine multiplexing of even six plasmids per 180 reads can still maintain high accuracy of consensus sequencing. SAVEMONEY should further democratize whole-plasmid sequencing by nanopore and related technologies, driving down the effective cost of whole-plasmid sequencing to lower than that of a single Sanger sequencing run.

    1. Biochemistry and Chemical Biology
    2. Computational and Systems Biology
    Shinichi Kawaguchi, Xin Xu ... Toshie Kai
    Research Article

    Protein–protein interactions are fundamental to understanding the molecular functions and regulation of proteins. Despite the availability of extensive databases, many interactions remain uncharacterized due to the labor-intensive nature of experimental validation. In this study, we utilized the AlphaFold2 program to predict interactions among proteins localized in the nuage, a germline-specific non-membrane organelle essential for piRNA biogenesis in Drosophila. We screened 20 nuage proteins for 1:1 interactions and predicted dimer structures. Among these, five represented novel interaction candidates. Three pairs, including Spn-E_Squ, were verified by co-immunoprecipitation. Disruption of the salt bridges at the Spn-E_Squ interface confirmed their functional importance, underscoring the predictive model’s accuracy. We extended our analysis to include interactions between three representative nuage components—Vas, Squ, and Tej—and approximately 430 oogenesis-related proteins. Co-immunoprecipitation verified interactions for three pairs: Mei-W68_Squ, CSN3_Squ, and Pka-C1_Tej. Furthermore, we screened the majority of Drosophila proteins (~12,000) for potential interaction with the Piwi protein, a central player in the piRNA pathway, identifying 164 pairs as potential binding partners. This in silico approach not only efficiently identifies potential interaction partners but also significantly bridges the gap by facilitating the integration of bioinformatics and experimental biology.