Human genetic analyses of organelles highlight the nucleus in age-related trait heritability
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).
-
GWAS Catalog, all associations v1.0.2NHGRI-EBI GWAS Catalog.
-
Neale Lab UKB Round 2 GWAS Summary StatisticsNeale lab, Broad Institute.
-
Genome-wide meta-analysis identifies 56 bone mineral density loci and reveals 14 loci associated with risk of fractureGEFOS Pooled Femoral Neck Summary Statistics.
-
Genetic associations at 53 loci highlight cell types and biological pathways relevant for kidney functionCKDGen Data at Medical Center - University of Freiburg, eGFRcrea and CKD.
-
Genomics of 1 million parent lifespans implicates novel pathways and common diseases and distinguishes survival chancesUniversity of Edinburgh DataShare, doi:10.7488/ds/2463.
-
Identification of 12 genetic loci associated with human healthspanZenodo, doi:10.5281/zenodo.1302861.
Article and author information
Author details
Funding
National Institutes of Health (T32AG000222)
- Rahul Gupta
National Institutes of Health (R35GM122455)
- Vamsi Mootha, MD MD
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
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,588
- views
-
- 260
- downloads
-
- 9
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
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)
Further reading
-
- Genetics and Genomics
Becker muscular dystrophy (BMD), an X-linked muscular dystrophy, is mostly caused by an in-frame deletion of Duchenne muscular dystrophy (DMD). BMD severity varies from asymptomatic to severe, associated with the genotype of DMD. However, the underlying mechanisms remain unclear. We established BMD mice carrying three representative exon deletions: ex45–48 del., ex45–47 del., and ex45–49 del. (d45–48, d45–47, and d45–49), with high frequencies and different severities in the human BMD hotspot. All three BMD mice showed muscle weakness, muscle degeneration, and fibrosis, but these changes appeared at different times for each exon deletion, consistent with the severities obtained by the natural history study of BMD. BMD mice showed site-specific muscle changes, unlike mdx mice, which showed diffuse muscle changes, and we demonstrated selective type IIa fiber reduction in BMD mice. Furthermore, BMD mice showed sarcolemmal neuronal nitric oxide synthase (nNOS) reduction and morphological capillary changes around type IIa fibers. These results suggest that capillary changes caused by nNOS reduction may be associated with the mechanism of skeletal muscle degeneration and type IIa fiber reduction in BMD mice. BMD mice may be useful in elucidating the pathomechanisms and developing vascular targeted therapies for human BMD.
-
- Genetics and Genomics
- Microbiology and Infectious Disease
Evolution of gene expression frequently drives antibiotic resistance in bacteria. We had previously (Patel and Matange, eLife, 2021) shown that, in Escherichia coli, mutations at the mgrB locus were beneficial under trimethoprim exposure and led to overexpression of dihydrofolate reductase (DHFR), encoded by the folA gene. Here, we show that DHFR levels are further enhanced by spontaneous duplication of a genomic segment encompassing folA and spanning hundreds of kilobases. This duplication was rare in wild-type E. coli. However, its frequency was elevated in a lon-knockout strain, altering the mutational landscape early during trimethoprim adaptation. We then exploit this system to investigate the relationship between trimethoprim pressure and folA copy number. During long-term evolution, folA duplications were frequently reversed. Reversal was slower under antibiotic pressure, first requiring the acquisition of point mutations in DHFR or its promoter. Unexpectedly, despite resistance-conferring point mutations, some populations under high trimethoprim pressure maintained folA duplication to compensate for low abundance DHFR mutants. We find that evolution of gene dosage depends on expression demand, which is generated by antibiotic and exacerbated by proteolysis of drug-resistant mutants of DHFR. We propose a novel role for proteostasis as a determinant of copy number evolution in antibiotic-resistant bacteria.