Lack of evidence for increased transcriptional noise in aged tissues
Abstract
Aging is often associated with a loss of cell type identity that results in an increase in transcriptional noise in aged tissues. If this phenomenon reflects a fundamental property of aging remains an open question. Transcriptional changes at the cellular level are best detected by single-cell RNA sequencing (scRNAseq). However, the diverse computational methods used for the quantification of age-related loss of cellular identity have prevented reaching meaningful conclusions by direct comparison of existing scRNAseq datasets. To address these issues we created Decibel, a Python toolkit that implements side-to-side four commonly used methods for the quantification of age-related transcriptional noise in scRNAseq data. Additionally, we developed Scallop, a novel computational method for the quantification of membership of single cells to their assigned cell type cluster. Cells with a greater Scallop membership score are transcriptionally more stable. Application of these computational tools to seven aging datasets showed large variability between tissues and datasets, suggesting that increased transcriptional noise is not a universal hallmark of aging. To understand the source of apparent loss of cell type identity associated with aging, we analyzed cell type-specific changes in transcriptional noise and the changes in cell type composition of the mammalian lung. No robust pattern of cell type-specific transcriptional noise alteration was found across aging lung datasets. In contrast, age-associated changes in cell type composition of the lung were consistently found, particularly of immune cells. These results suggest that claims of increased transcriptional noise of aged tissues should be reformulated.
Data availability
Code availabilityThe Decibel and Scallop repositories can be found at https://gitlab.com/olgaibanez/decibel and https: //gitlab.com/olgaibanez/scallop, respectively. The reproducible Jupyter notebooks with the analyses carried out in this study can be found in figshare (https://doi.org/10.6084/m9.figshare.20402817.v1).
-
scRNAseq dataset of murine aging lungGene Expression Omnibus, GSE124872.
-
scRNAseq dataset of murine aging lung, spleen and kidneyGene Expression Omnibus, GSE132901.
-
Human Lung Cell AtlasSynapse, syn21041850.
-
scRNAseq datasets of adult mammalian lungsGene Expression Omnibus, GSE133747.
-
scRNAseq dataset of human aging pancreasGene Expression Omnibus, GSE81547.
-
scRNAseq dataset of human aging skinGene Expression Omnibus, GSE130973.
-
scRNAseq dataset of murine aging brainGene Expression Omnibus, GSE129788.
-
scRNAseq dataset of murine aging dermal fibroblastsGene Expression Omnibus, GSE111136.
Article and author information
Author details
Funding
la Caixa" Foundation " (LCF/BQ/IN18/11660065)
- Olga Ibañez-Solé
Instituto de Salud Carlos III (AC17/00012)
- Olga Ibañez-Solé
- Alex M Ascensión
- Ander Izeta
Instituto de Salud Carlos III (PI19/01621)
- Olga Ibañez-Solé
- Alex M Ascensión
- Ander Izeta
Ministerio de Ciencia e Innovación (PID2020-119715GB-I00)
- Olga Ibañez-Solé
- Alex M Ascensión
- Ander Izeta
European Regional Development Fund (MCIN/AEI/10.13039/501100011033)
- Olga Ibañez-Solé
- Alex M Ascensión
- Ander Izeta
H2020 Marie Skłodowska-Curie Actions (713673)
- Olga Ibañez-Solé
Eusko Jaurlaritza (PRE_2020_2_0081)
- Alex M Ascensión
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Copyright
© 2022, Ibañez-Solé 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
-
- 3,808
- views
-
- 386
- downloads
-
- 22
- 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
-
- Computational and Systems Biology
Protein engineering is a pivotal aspect of synthetic biology, involving the modification of amino acids within existing protein sequences to achieve novel or enhanced functionalities and physical properties. Accurate prediction of protein variant effects requires a thorough understanding of protein sequence, structure, and function. Deep learning methods have demonstrated remarkable performance in guiding protein modification for improved functionality. However, existing approaches predominantly rely on protein sequences, which face challenges in efficiently encoding the geometric aspects of amino acids’ local environment and often fall short in capturing crucial details related to protein folding stability, internal molecular interactions, and bio-functions. Furthermore, there lacks a fundamental evaluation for developed methods in predicting protein thermostability, although it is a key physical property that is frequently investigated in practice. To address these challenges, this article introduces a novel pre-training framework that integrates sequential and geometric encoders for protein primary and tertiary structures. This framework guides mutation directions toward desired traits by simulating natural selection on wild-type proteins and evaluates variant effects based on their fitness to perform specific functions. We assess the proposed approach using three benchmarks comprising over 300 deep mutational scanning assays. The prediction results showcase exceptional performance across extensive experiments compared to other zero-shot learning methods, all while maintaining a minimal cost in terms of trainable parameters. This study not only proposes an effective framework for more accurate and comprehensive predictions to facilitate efficient protein engineering, but also enhances the in silico assessment system for future deep learning models to better align with empirical requirements. The PyTorch implementation is available at https://github.com/ai4protein/ProtSSN.
-
- Computational and Systems Biology
- Microbiology and Infectious Disease
The Staphylococcus aureus clonal complex 8 (CC8) is made up of several subtypes with varying levels of clinical burden; from community-associated methicillin-resistant S. aureus USA300 strains to hospital-associated (HA-MRSA) USA500 strains and ancestral methicillin-susceptible (MSSA) strains. This phenotypic distribution within a single clonal complex makes CC8 an ideal clade to study the emergence of mutations important for antibiotic resistance and community spread. Gene-level analysis comparing USA300 against MSSA and HA-MRSA strains have revealed key horizontally acquired genes important for its rapid spread in the community. However, efforts to define the contributions of point mutations and indels have been confounded by strong linkage disequilibrium resulting from clonal propagation. To break down this confounding effect, we combined genetic association testing with a model of the transcriptional regulatory network (TRN) to find candidate mutations that may have led to changes in gene regulation. First, we used a De Bruijn graph genome-wide association study to enrich mutations unique to the USA300 lineages within CC8. Next, we reconstructed the TRN by using independent component analysis on 670 RNA-sequencing samples from USA300 and non-USA300 CC8 strains which predicted several genes with strain-specific altered expression patterns. Examination of the regulatory region of one of the genes enriched by both approaches, isdH, revealed a 38-bp deletion containing a Fur-binding site and a conserved single-nucleotide polymorphism which likely led to the altered expression levels in USA300 strains. Taken together, our results demonstrate the utility of reconstructed TRNs to address the limits of genetic approaches when studying emerging pathogenic strains.