1. Genetics and Genomics

Human genetic analyses of organelles highlight the nucleus in age-related trait heritability

  1. Rahul Gupta  Is a corresponding author
  2. Konrad J Karczewski
  3. Daniel Howrigan
  4. Benjamin M Neale  Is a corresponding author
  5. Vamsi K Mootha MD  Is a corresponding author
  1. Howard Hughes Medical Institute, Massachusetts General Hospital, United States
  2. Broad Institute of MIT and Harvard, United States
  3. Massachusetts General Hospital, United States
https://doi.org/10.7554/eLife.68610
Share this article
Cite this article
  1. Rahul Gupta
  2. Konrad J Karczewski
  3. Daniel Howrigan
  4. Benjamin M Neale
  5. Vamsi K Mootha MD
(2021)
Human genetic analyses of organelles highlight the nucleus in age-related trait heritability
eLife 10:e68610.
https://doi.org/10.7554/eLife.68610
  2. Download BibTeX
  3. Download .RIS

Abstract

Most age-related human diseases are accompanied by a decline in cellular organelle integrity, including impaired lysosomal proteostasis and defective mitochondrial oxidative phosphorylation. An open question, however, is the degree to which inherited variation in or near genes encoding each organelle contributes to age-related disease pathogenesis. Here, we evaluate if genetic loci encoding organelle proteomes confer greater-than-expected age-related disease risk. As mitochondrial dysfunction is a 'hallmark' of aging, we begin by assessing nuclear and mitochondrial DNA loci near genes encoding the mitochondrial proteome and surprisingly observe a lack of enrichment across 24 age-related traits. Within nine other organelles, we find no enrichment with one exception: the nucleus, where enrichment emanates from nuclear transcription factors. In agreement, we find that genes encoding several organelles tend to be 'haplosufficient', while we observe strong purifying selection against heterozygous protein-truncating variants impacting the nucleus. Our work identifies common variation near transcription factors as having outsize influence on age-related trait risk, motivating future efforts to determine if and how this inherited variation then contributes to observed age-related organelle deterioration.

Data availability

Heritability point estimates and standard errors for age-related traits are listed in Supplementary File 1. Genetic and phenotypic correlation point estimates and standard errors/p-values plotted in Figure 1B are available in Figure 1-Source data 1. Summary statistics from mtDNA-GWAS (plotted in Figure 2 and Figure 2-Figure supplement 9) are available in Source data 2. All gene-based enrichment analysis p-values and point estimates are available in Source data 1 and Source data 3. Period prevalence data for diseases in the UK can be obtained from Kuan et al. 2019. Gene-sets can be found using COMPARTMENTS (https://compartments.jensenlab.org), MitoCarta 2.0 (https://www.broadinstitute.org/files/shared/metabolism/mitocarta/human.mitocarta2.0.html), Lambert et al. 2018 (DOI: 10.1016/j.cell.2018.01.029), Frazier et al. 2019 (DOI: 10.1074/jbc.R117.809194), Finucane et al. 2018 (https://alkesgroup.broadinstitute.org/LDSCORE/), Kapopoulou et al. 2015 (DOI: 10.1111/evo.12819), and the Macarthur laboratory (https://github.com/macarthur-lab/gene_lists). Gene age estimates were obtained from Litman, Stein 2019 (DOI: 10.1053/j.seminoncol.2018.11.002). GWAS catalog annotations can be obtained from: https://www.ebi.ac.uk/gwas. Heritability estimates across UKB can be obtained at: https://nealelab.github.io/UKBB_ldsc/. UKB summary statistics can be obtained from Neale lab GWAS round 2: https://github.com/Nealelab/UK_Biobank_GWAS. Annotations for the Baseline v1.1 and BaselineLD v2.2 models as well as other relevant reference data, including the 1000G EUR reference panel, can be obtained from https://alkesgroup.broadinstitute.org/LDSCORE/. eQTL and expression data in human tissues can be obtained from GTEx (https://www.gtexportal.org). Constraint estimates can be found via gnomAD: https://gnomad.broadinstitute.org. See citations for publicly available GWAS meta-analysis summary statistics (28,29,51,52,30-37).

The following previously published data sets were used

Article and author information

Author details

  1. Rahul Gupta

    Department of Molecular Biology, Howard Hughes Medical Institute, Massachusetts General Hospital, Boston, United States
    For correspondence
    rahul@broadinstitute.org
    Competing interests
    No competing interests declared.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-8263-2455

  2. Konrad J Karczewski

    Broad Institute of MIT and Harvard, Cambridge, United States
    Competing interests
    Konrad J Karczewski, is a consultant for Vor Biopharma..

  3. Daniel Howrigan

    Analytic & Translational Genetics Unit, Massachusetts General Hospital, Boston, United States
    Competing interests
    No competing interests declared.

  4. Benjamin M Neale

    Broad Institute of MIT and Harvard, Cambridge, United States
    For correspondence
    bneale@broadinstitute.org
    Competing interests
    Benjamin M Neale, BMN is a member of the scientific advisory board at Deep Genomics and RBNC Therapeutics. BMN is a consultant for Camp4 Therapeutics, Takeda Pharmaceutical and Biogen..

  5. Vamsi K Mootha MD

    Molecular Biology, Howard Hughes Medical Institute, Massachusetts General Hospital, Boston, United States
    For correspondence
    vamsi@hms.harvard.edu
    Competing interests
    Vamsi K Mootha, is an advisor to and receives compensation or equity from Janssen Pharmaceuticals, 5am Ventures, and Raze Therapeutics..
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-9924-642X

Funding

National Institutes of Health (T32AG000222)

  • Rahul Gupta

National Institutes of Health (R35GM122455)

  • Vamsi K Mootha MD

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

Reviewing Editor

  1. Sara Hägg, Karolinska Institutet, Sweden

Version history

  1. Preprint posted: January 22, 2021 (view preprint)
  2. Received: March 21, 2021
  3. Accepted: August 30, 2021
  4. Accepted Manuscript published: September 1, 2021 (version 1)
  5. Accepted Manuscript updated: September 7, 2021 (version 2)
  6. Version of Record published: September 27, 2021 (version 3)

Copyright

© 2021, Gupta 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

  • 2,412
    views
  • 250
    downloads
  • 7
    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. Rahul Gupta
  2. Konrad J Karczewski
  3. Daniel Howrigan
  4. Benjamin M Neale
  5. Vamsi K Mootha MD
(2021)
Human genetic analyses of organelles highlight the nucleus in age-related trait heritability
eLife 10:e68610.
https://doi.org/10.7554/eLife.68610

Share this article

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

Categories and tags

Research organism

Further reading

    1. Computational and Systems Biology
    2. Genetics and Genomics

    Haplotype function score improves biological interpretation and cross-ancestry polygenic prediction of human complex traits

    Weichen Song, Yongyong Shi, Guan ning Lin
    Tools and Resources

    We propose a new framework for human genetic association studies: at each locus, a deep learning model (in this study, Sei) is used to calculate the functional genomic activity score for two haplotypes per individual. This score, defined as the Haplotype Function Score (HFS), replaces the original genotype in association studies. Applying the HFS framework to 14 complex traits in the UK Biobank, we identified 3619 independent HFS–trait associations with a significance of p < 5 × 10−8. Fine-mapping revealed 2699 causal associations, corresponding to a median increase of 63 causal findings per trait compared with single-nucleotide polymorphism (SNP)-based analysis. HFS-based enrichment analysis uncovered 727 pathway–trait associations and 153 tissue–trait associations with strong biological interpretability, including ‘circadian pathway-chronotype’ and ‘arachidonic acid-intelligence’. Lastly, we applied least absolute shrinkage and selection operator (LASSO) regression to integrate HFS prediction score with SNP-based polygenic risk scores, which showed an improvement of 16.1–39.8% in cross-ancestry polygenic prediction. We concluded that HFS is a promising strategy for understanding the genetic basis of human complex traits.

    1. Genetics and Genomics
    2. Immunology and Inflammation

    Transposable elements regulate thymus development and function

    Jean-David Larouche, Céline M Laumont ... Claude Perreault
    Research Article

    Transposable elements (TEs) are repetitive sequences representing ~45% of the human and mouse genomes and are highly expressed by medullary thymic epithelial cells (mTECs). In this study, we investigated the role of TEs on T-cell development in the thymus. We performed multiomic analyses of TEs in human and mouse thymic cells to elucidate their role in T-cell development. We report that TE expression in the human thymus is high and shows extensive age- and cell lineage-related variations. TE expression correlates with multiple transcription factors in all cell types of the human thymus. Two cell types express particularly broad TE repertoires: mTECs and plasmacytoid dendritic cells (pDCs). In mTECs, transcriptomic data suggest that TEs interact with transcription factors essential for mTEC development and function (e.g., PAX1 and REL), and immunopeptidomic data showed that TEs generate MHC-I-associated peptides implicated in thymocyte education. Notably, AIRE, FEZF2, and CHD4 regulate small yet non-redundant sets of TEs in murine mTECs. Human thymic pDCs homogenously express large numbers of TEs that likely form dsRNA, which can activate innate immune receptors, potentially explaining why thymic pDCs constitutively secrete IFN ɑ/β. This study highlights the diversity of interactions between TEs and the adaptive immune system. TEs are genetic parasites, and the two thymic cell types most affected by TEs (mTEcs and pDCs) are essential to establishing central T-cell tolerance. Therefore, we propose that orchestrating TE expression in thymic cells is critical to prevent autoimmunity in vertebrates.

    1. Genetics and Genomics

    Deciphering the complex relationship between type 2 diabetes and fracture risk with both genetic and observational evidence

    Pianpian Zhao, Zhifeng Sheng ... Hou-Feng Zheng
    Research Article

    The ‘diabetic bone paradox’ suggested that type 2 diabetes (T2D) patients would have higher areal bone mineral density (BMD) but higher fracture risk than individuals without T2D. In this study, we found that the genetically predicted T2D was associated with higher BMD and lower risk of fracture in both weighted genetic risk score (wGRS) and two-sample Mendelian randomization (MR) analyses. We also identified ten genomic loci shared between T2D and fracture, with the top signal at SNP rs4580892 in the intron of gene RSPO3. And the higher expression in adipose subcutaneous and higher protein level in plasma of RSPO3 were associated with increased risk of T2D, but decreased risk of fracture. In the prospective study, T2D was observed to be associated with higher risk of fracture, but BMI mediated 30.2% of the protective effect. However, when stratified by the T2D-related risk factors for fracture, we observed that the effect of T2D on the risk of fracture decreased when the number of T2D-related risk factors decreased, and the association became non-significant if the T2D patients carried none of the risk factors. In conclusion, the genetically determined T2D might not be associated with higher risk of fracture. And the shared genetic architecture between T2D and fracture suggested a top signal around RSPO3 gene. The observed effect size of T2D on fracture risk decreased if the T2D-related risk factors could be eliminated. Therefore, it is important to manage the complications of T2D to prevent the risk of fracture.