A concerted increase in readthrough and intron retention drives transposon expression during aging and senescence
eLife assessment
This study presents fundamental findings on the role of transcription readout and intron retention in transposon expression during aging in mammals. The evidence supporting the claims of the authors is compelling, strongly supporting the authors' claims. The work will be of interest to scientists studying aging, transcription regulation, and epigenetics.
https://doi.org/10.7554/eLife.87811.3.sa0Fundamental: Findings that substantially advance our understanding of major research questions
- Landmark
- Fundamental
- Important
- Valuable
- Useful
Compelling: Evidence that features methods, data and analyses more rigorous than the current state-of-the-art
- Exceptional
- Compelling
- Convincing
- Solid
- Incomplete
- Inadequate
During the peer-review process the editor and reviewers write an eLife Assessment that summarises the significance of the findings reported in the article (on a scale ranging from landmark to useful) and the strength of the evidence (on a scale ranging from exceptional to inadequate). Learn more about eLife Assessments
Abstract
Aging and senescence are characterized by pervasive transcriptional dysfunction, including increased expression of transposons and introns. Our aim was to elucidate mechanisms behind this increased expression. Most transposons are found within genes and introns, with a large minority being close to genes. This raises the possibility that transcriptional readthrough and intron retention are responsible for age-related changes in transposon expression rather than expression of autonomous transposons. To test this, we compiled public RNA-seq datasets from aged human fibroblasts, replicative and drug-induced senescence in human cells, and RNA-seq from aging mice and senescent mouse cells. Indeed, our reanalysis revealed a correlation between transposons expression, intron retention, and transcriptional readthrough across samples and within samples. Both intron retention and readthrough increased with aging or cellular senescence and these transcriptional defects were more pronounced in human samples as compared to those of mice. In support of a causal connection between readthrough and transposon expression, analysis of models showing induced transcriptional readthrough confirmed that they also show elevated transposon expression. Taken together, our data suggest that elevated transposon reads during aging seen in various RNA-seq dataset are concomitant with multiple transcriptional defects. Intron retention and transcriptional readthrough are the most likely explanation for the expression of transposable elements that lack a functional promoter.
Introduction
Repetitive DNA makes up at least 50% of the human genome (de Koning et al., 2011) and has been linked to genomic instability, cancer and – somewhat less consistently – species longevity (Pabis, 2021; Khristich and Mirkin, 2020). Transposons are one particularly abundant class of repetitive sequences found in the nuclear genome of eukaryotes and are strongly associated with cancer-causing structural variants like DNA deletions, insertions, or inversions (Rodriguez-Martin et al., 2020). Although transposons are often considered parasitic, harmful or, at best, neutral recent evidence suggests they may contribute to genomic diversification on evolutionary timescales and to the cellular stress-response by providing transcription factor-binding sites (Villanueva-Cañas et al., 2019).
The major transposon families in the human genome are long interspersed nuclear element (LINE), short interspersed nuclear element (SINE), long terminal repeat (LTR) transposons, and DNA transposons comprising 21%, 11%, 8%, and 3% of the human genome, respectively (Kazazian and Moran, 2017).
In recent years, transposon expression, particularly of elements belonging to the LINE-1 family, has been hypothesized to be a near universal marker of aging. Several lines of evidence support a link between transposon reactivation and aging (Gorbunova et al., 2021).
Using RNA-seq, for example, to quantify reads mapping to transposable elements it was shown that expression of these was higher in 10- vs 1-day-old nematodes (LaRocca et al., 2020), 40- vs 10-day-old flies (Wood et al., 2016) and in muscle and liver of aged mice (De Cecco et al., 2019). Transposon expression also increases during in vitro senescence (Colombo et al., 2018) and in fibroblasts isolated from people between the ages of 1 and 96 years (LaRocca et al., 2020), whereas lifespan extending mutations and interventions in mice reduce transposon expression (Wahl et al., 2021).
However, not all studies report such consistent age-related increases across all transposon classes (Ghanam et al., 2019). Moreover, it would be a mistake to conflate changes in RNA-seq-, RNA-, DNA-, and protein-based measurements of transposon expression as evidence for one and the same phenomenon.
Quantification of transposon expression using RNA-seq or PCR techniques remains challenging given their repetitive nature and the fact that 98–99% of all transposons are co-expressed with neighboring transcriptional units (Deininger et al., 2017; Stow et al., 2021). In fact, only a small fraction of LINE-1 and SINE elements, the latter relying on co-transposition by LINEs, are expressed from a functional promoter. It is believed that such ‘hot’ LINE-1 loci drive most transposition events and genomic instability (Deininger et al., 2017; Rodriguez-Martin et al., 2020).
The large number of co-expressed transposons would likely mask any signal from autonomous elements in standard RNA-seq experiments. Therefore, alternative explanations are needed for increased transposon expression in such datasets and these could include age-related changes in transposon adjacent genes for example intron retention and transcriptional readthrough.
Most eukaryotic genes contain introns which need to be removed during a complex co-transcriptional process called splicing. Defects in splicing are the cause of many, often neurologic, hereditary conditions (Scotti and Swanson, 2016). One such type of splicing defect is intron retention. Although basal levels of intron retention may be benign or even physiologic, high levels are considered harmful (Zheng et al., 2020). Aging is accompanied by many changes to splicing patterns (Wang et al., 2018), including intron retention, which is increased during cellular senescence (Yao et al., 2020) and in aging Drosophila, mouse hippocampus and human prefrontal cortex. In addition, levels in human AD brain tissues are elevated even further than in the aged brain (Adusumalli et al., 2019).
Transcription is terminated by the cleavage and polyadenylation machinery when the polymerase transcribes through the polyadenylation signal leading to pre-mRNA cleavage, polymerase pausing, conformational change, and eventually polymerase detachment (Rosa-Mercado et al., 2021). However, this termination may fail leading to so-called readthrough transcription which increases during stress conditions (Rosa-Mercado and Steitz, 2022) and during cellular senescence, at least for a subset of genes (Muniz et al., 2017).
Emerging evidence suggests that intron retention, readthrough, and transposon expression are linked (Rosa-Mercado et al., 2021; Hadar et al., 2022), but this has not been studied in the context of aging. In this manuscript, we show that intron retention and readthrough taken together can explain most of the apparent increase in transposon expression seen in RNA-seq datasets of aging.
Results
Genomic location of transposons implicates them as a marker of multiple transcriptional defects
During aging and senescence the expression of all four major transposon families (LTR, DNA, LINE, and SINE) increases. However, such a global increase in transposon expression cannot be explained by transposon biology alone since most transposons are ancient, and, having accumulated many mutations, lack a functional promoter (Deniz et al., 2019). Instead, we reasoned that transposon expression is associated with the expression of nearby genes.
Consistently, out of 4.5 million annotated transposons 60% were intragenic (most of them intronic), 10% within 5 kb of a gene and the remaining 30% were intergenic, albeit usually still within 100 kb (Table 1). For transposons whose expression significantly changed with age or senescence, the number of intronic and downstream-of-gene transposons is even higher.
Based on these findings, we studied intronic transposons as a potential marker for age-related intron retention and downstream-of-gene transposons as a marker of age-related readthrough transcription.
Elevated transposon expression and intron retention with age
To confirmed that age-related transposon expression is indeed elevated as shown by LaRocca et al., 2020, we analyzed the transcriptomic dataset originally generated by Fleischer et al., 2018, which included 143 fibroblasts isolated from donors between the ages of 0–96. This dataset will be hereafter referred to as the ‘aging dataset’. We find that transposon expression increases slowly and non-linearly during aging, with a sudden and sharp increase around 80 years of age (Figure 1A, Figure 1—figure supplement 1) and this increase was highly strand specific (Figure 1—figure supplement 2).
Next, we found that intron retention is elevated with aging (Figure 1B), consistent with the study by Yao et al., 2020, and that samples with high transposon expression also showed elevated intron retention (Figure 1C).
In our analysis, we corrected for gene co-expression by normalizing the expression levels of introns and transposons to the expression levels of the closest gene. Another complementary approach would be to plot and compare the raw, unadjusted log-fold changes for genes, introns, transposons, and other elements. In this analysis, we found that, on average, age-related upregulation of transposons and introns is more pronounced than upregulation of genes (Figure 1—figure supplement 3A), consistent with Figure 1A, B.
Next, we reanalyzed a dataset originally published by Colombo et al., 2018 which will be referred to as the ‘senescence dataset’. Here, we tested whether MDAH041 cells induced into senescence with H2O2, 5-azacytidine, adriamycin, or serial passaging show elevated transposon expression and intron retention. Consistent with the original findings, all these treatments increased the expression of transposons (Figure 1D, Figure 1—figure supplement 3B). The retention of introns was also increased (Figure 1E) and expression levels of introns and transposons correlated with each other (Figure 1F).
Having established that transposons reads and retained introns are elevated with aging and senescence, we tested whether this is conserved between mouse and human. Transposon expression (p < 0.0001) was elevated, and intron retention showed a trend toward elevation (p = 0.14), in the liver of 26-month-old mice, suggesting these age-related transcriptional defects might be conserved between species (Figure 1—figure supplement 4A, B; Figure 1—figure supplement 5). In contrast, we found that increased expression of transposons and intron retention in models of cellular senescence is not conserved between species. To the contrary, early passage mouse fibroblasts showed elevated transposon expression and intron retention compared to senescent fibroblasts (Figure 1—figure supplement 4C, D), which contrasts with the mouse liver and human senescence data.
Transposon expression as a marker of intron retention
So far these age-related events have been considered independently, but given the localization of transposons being mostly intragenic we reasoned that intron retention could bias the transposon findings. Indeed, we found that samples with high intron retention also show high transposon expression. This is true both in the aging dataset (Figure 1C) and in the cellular senescence dataset (Figure 1F), where intronic reads accounted for most of the observed variability in transposon expression (or vice versa). Looking at individual genomic loci, we also found a strong correlation between the log-fold change of a transposon and the log-fold change of the intron it is located in (Figure 1—figure supplement 6).
As a more direct test of this hypothesis, we determined the correlation between the counts of retained introns and the counts of intronic transposons for each sample individually. Non-intronic transposons and randomized transposons served as a negative control. Consistent with the log-fold change data, we found a specific correlation between intronic transposon and intron counts in both datasets, with stronger correlations in the senescence dataset (Figure 1—figure supplement 7).
Although causality is hard to establish from our datasets, we also tested whether transposon expression could affect intron retention. If transposon expression was causally linked to intron retention, the most likely mechanism would be via an impairment of the splice junction. In this case, we would expect more transposons near the splice junctions of retained (differentially expressed) introns. However, we observed a depletion of transposons at the splice junctions of all introns (Figure 1—figure supplement 8A) and of differentially expressed introns (Figure 1—figure supplement 8B, E). Transposon densities were also comparable between these two classes of introns (Figure 1—figure supplement 8C, F). Moreover, introns with and without expressed transposons showed similar levels of age-related upregulation (Figure 1—figure supplement 9). Based on our data, we propose a model (Figure 2) where intron retention contributes to age-related upregulation of transposons, rather than transposons contributing to intron retention.
Expression of extragenic transposons as a marker of readthrough transcription
Most transposons showing increased expression with aging or cellular senescence mapped within genes or introns and their expression changes may be attributed to intron retention. How can we explain expression of the remaining 14–30% of transposons (Table 1)?
To answer this question, we analyzed the distribution of transposons mapping to extragenic regions. Out of all the transposons, showing significantly changed expression during aging or cellular senescence, most were located within 5 kb of genes, primarily downstream of genes (Figure 3A). The distribution of transposons in the aging and senescence dataset was very similar except for a wider distribution of senescence-associated transposons within a 100-kb window around genes (Figure 3B). This contrasts with the almost completely even distribution of randomly permuted transposons.
Given the proximity of significantly changed transposons and genes, this suggested that an age-related increase in transcriptional readthrough could explain elevated transposon expression. To test this, we next quantified readthrough in a 10- to 20-kb region downstream of genes. We chose a region not immediately adjacent to genes to avoid biases through basal levels of readthrough, although the data were comparable for closer distances (Figure 3—figure supplement 1).
Indeed, readthrough transcription was increased with aging (Figure 4A, B), and in the senescence dataset, levels of readthrough transcription increased with replicative or chemically induced senescence (Figure 4D, E). Transposon expression correlated with readthrough levels on a per sample basis in both the aging and the senescence dataset (Figure 4C, F).
As shown before for intron retention, these findings may not translate to other species. We found that readthrough was unchanged in the liver of 26-month-old mice (p = 0.22, Figure 4—figure supplement 1A). In addition, readthrough decreased with replicative senescence of mouse fibroblasts, contrary to our expectations (Figure 4—figure supplement 1C, D).
Next, we quantified readthrough across a 100-kb region downstream of genes. We found that both aging and senescence lead to readthrough across the whole region, although senescence was associated with more extensive readthrough further downstream of genes (Figure 4—figure supplement 2, Figure 3B). As another line of evidence we quantified transcription upstream of genes, which is called read-in. This kind of transcription is often attributed to widespread readthrough extending from the end of a gene to the start of another (Roth et al., 2020). Consistent with the extensive readthrough levels observed we found elevated read-in transcription during cellular senescence (Figure 4—figure supplement 3A) but not aging (Figure 4—figure supplement 3B).
Transposon counts followed a similar decaying pattern to readthrough counts downstream of genes (Figure 4—figure supplement 4) and genes with high age-related readthrough were associated with strong age-related upregulation of gene proximal transposons (Figure 4—figure supplement 5). Intergenic transposons showed such low expression levels that their contribution to total counts was negligible (Figure 4—figure supplement 6).
Having shown that samples with higher levels of readthrough also have higher transposon expression (per sample correlation), we asked whether transposons are more often than expected located in readthrough regions (per locus correlation). Indeed as expected, in both datasets, transposons significantly changed with age or senescence were more often found in significant readthrough regions compared to all transposons. Moreover, there was a significant correlation between the log-fold changes of transposons downstream of genes with the log-fold changes of readthrough downstream of the same genes (Figure 4—figure supplement 7). Transposons located elsewhere showed a weaker correlation with readthrough.
To further test whether these correlations between log-fold changes are also reflected on the count level, we determined the correlation between readthrough counts and the counts of transposons downstream of genes for each sample individually. Transposons located upstream of genes and randomized transposons served as a negative control. Consistent with the above data, we found a specific correlation between transposon and readthrough counts in both datasets (Figure 4—figure supplement 8).
Finally, we asked whether readthrough levels are predictive of transposon expression on a per sample basis independently of intron retention. Using the aging dataset with 143 samples we compared a linear model just including intron retention with one also including readthrough and found the second model to perform significantly better (p < 0.0001). In addition, we found that the correlation between readthrough levels and transposon expression remained significant even when correcting for intron retention using partial correlation (R = 0.53, p < 0.0001). We did not test this in the senescence dataset due to its limited size.
Readthrough is a causal inducer of transposon expression
Based on this correlative evidence (Figure 4), we hypothesized that specific inducers of readthrough would also increase transposon expression. To test this, we reanalyzed literature data from studies involving multiple treatments that promote readthrough. These included viral infection, hyperosmotic stress (KCl treatment), heat shock, and oxidative stress.
We found that KCl treatment of HEK293 cells for 1 hr (Rosa-Mercado et al., 2021) indeed strongly induced both readthrough (Figure 5A) and transposon expression (Figure 5B) as determined by RNA-seq and TTLseq. Consistent with this, in our reanalysis of a study using mNET-seq profiling of nascent mRNA (Bauer et al., 2018), KCl treatment of HEK293 cells also induced readthrough (Figure 5D) and transposon expression (Figure 5E), whereas influenza infection caused an intermediate phenotype. Finally, we analyzed RNA-seq data from NIH-3T3 mouse fibroblasts subjected to three stressors over 2 hr (Vilborg et al., 2017). Heatshock caused extensive readthrough (Figure 5—figure supplement 1A) and transposon expression (Figure 5—figure supplement 1B), whereas hydrogen peroxide treatment and KCl caused moderate expression of these.
We also reasoned that readthrough should specifically increase the expression of transposons that are close to genes. Consistent with this notion, induced readthrough led to preferential expression of gene proximal transposons (i.e. those within 25 kb of genes), when compared with senescence or aging (Figure 5—figure supplement 2A). Moreover, the fraction of significantly upregulated transposons found outside of genes was higher with induced readthrough than with senescence or aging (Figure 5—figure supplement 2B).
Increased expression of autonomous LINE-1 elements with senescence but not aging
Many expressed LINE-1 elements are found within genes, which poses a challenge for the quantification of active LINE-1s. Given this proximity it is unsurprising that a substantial correlation between the expression of transposons and neighboring genes (Figure 6) is seen in RNA-seq.
To mitigate this problem, we selected 250 transposons that showed higher than expected expression levels compared to the nearest gene. We reasoned that LINE-1s with a functional promoter (‘active’ LINE-1s) would be more likely autonomous and should be enriched among transposons that show a divergent expression profile from their neighboring genes.
We found no enrichment for active LINE-1 elements in the aging dataset (Figure 6A). In contrast, transposons significantly upregulated with senescence and, among those, especially the ones with higher-than-expected expression, showed significant enrichment for active LINE-1s (Figure 6B). Moreover, we found that almost 100% of active LINE-1s were increased with senescence which was higher than for all LINE-1s and other transposons.
Finally, we performed the same analysis with a different definition of active LINE-1s as those potentially encoding functional ORF1p and ORF2p. In this analysis, we included all elements encoding an open reading frame >900 nt with high sequence similarity to ORF1p and >3400 nt with high sequence similarity to ORF2p (p < 0.05 for pairwise alignment). Out of 122,626 annotated LINE-1s 585 encoded a putative ORF1p and 133 encoded an ORF2p.
ORF1p-encoding LINE-1s were significantly enriched among transposons with higher-than-expected expression and among transposons significantly changed with senescence or aging. These findings were more robust for senescence than for aging. ORF2p-encoding LINE-1s, in contrast, were too rare to be studied reliably as to their enrichment.
Discussion
We show that transposons would serve as a useful biomarker of aging, even if most of the transposon signal in RNA-seq was due to co-expression with nearby transcriptional units, because they capture different kinds of transcriptional dysfunction.
Introns
Intron retention is one prominent example of dysfunctional splicing during aging, but the association with age-related transposon expression has not been evaluated. A recent study by Gualandi et al., 2022 using RNA-seq data from the 1000 Genomes project found a strong per sample correlation between intron and transposon expression in the general population. In line with these results, we showed that intron retention may explain a large fraction of transposon reads in aging and senescence datasets (Figure 1).
Readthrough, read-in, and transposons
In a study of oncogene-induced senescence, Muniz et al., 2017 found increased readthrough transcription at 91 so-called START loci, which was associated with suppression of neighboring genes through the production of antisense RNAs and several of these antisense RNAs were necessary for the maintenance of the senescent state. Our data are consistent with their model, although we find even more extensive readthrough transcription. Analyzing hundreds of loci we show that senescence promotes transcriptional readthrough leading to the overexpression of intergenic transposons (>25 kb from the nearest gene) and read-in at neighboring genes (Figure 3B, Figure 4—figure supplement 3B). In contrast, aging-related readthrough was qualitatively different, producing shorter readthrough and no read-in (Figure 3B, Figure 4—figure supplement 3A).
Our correlative data suggested that readthrough transcription could account for the expression of many transposons that are outside of genes and thus cannot be explained by intron retention (Figure 4C, F). To test this hypothesis, we reanalyzed data from models of induced readthrough showing that readthrough is indeed a causal inducer for the expression of gene proximal transposons after hyperosmotic stress and viral infection (Figure 5). Given the technical challenges of measuring readthrough and the unreliability of intron retention detection tools (David et al., 2022), our results are likely underestimating the true magnitude of correlation between readthrough and transposon expression.
In a previous study, it was observed that senescence leads to increased expression of transposons and of genes with low abundances (Zhang et al., 2021). This ‘leakage’ was attributed to loss of H3K9me3 histone marks and increased heterochromatin accessibility, but is also consistent with readthrough transcription, as we observed. Indeed, it has been suggested that readthrough is the most likely cause for apparent intergenic transcription (Agostini et al., 2021).
Is there genuine expression of autonomous transposons during aging?
The distribution of transposons showing apparent age-related expression within the genome is an important question since intragenic transposons are more likely to be co-expressed with genes. Using Sanger sequencing of long-range amplicons to study replicative senescence, De Cecco et al., 2019 found a large fraction of transposons mapping to intergenic (69%) regions and few to intragenic/gene proximal regions (31%). In contrast using RNA-seq data, we find that most transposons are located within genes (≥70%; Table 1), although, there are more intergenic transposons significantly changed with senescence than with aging (Figure 3). One explanation for this discrepancy is the mapping approach used. Since their analysis did not use an expectation–maximization algorithm to map back transposons onto the genome, it may underestimate the number of intragenic and gene-adjacent transposons in contrast to our analysis using TELocal.
Given that inactive LINE-1 loci greatly outnumber active loci (Deininger et al., 2017; Stow et al., 2021), and we saw only limited evidence for active LINE-1 transcription in fibroblasts isolated from donors of varying ages (Figure 6A), our data appear to question the role of transposons in organismal aging. An intriguing explanation is provided by De Cecco et al., 2019, who found that LINE-1 ORF1 protein-positive cells colocalize with a small subset of senescent cells. Most likely, such changes would not be obvious from bulk sequencing, perhaps explaining why we did not see elevated LINE-1 expression in the aging dataset. In contrast, the signal would be readily detectable in cultures induced into senescence where most cells do enter senescence, which is consistent with our findings (Figure 6B).
As for another explanation, we speculate that LINE-1 elements expressed through intron retention and readthrough are compatible with translation via leaky scanning or re-initiation (Mouilleron et al., 2016), accounting for increased LINE-1 ORF1 protein levels detected in senescent cells, even if leaky scanning is orders of magnitude less efficient than normal initiation.
We did not study the expression of LTR or SINE elements in detail, which have been implicated in organismal aging (Liu et al., 2023) and age-related diseases such as macular degeneration (de Koning et al., 2011), because we lack consensus annotations for active and inactive members of these families.
Summary and limitations
Taken together, our work and existing literature suggest that RNA-seq, and by extension qPCR, of transposable elements overestimates their true expression and their age-related changes. Such data alone are insufficient to support a ‘transposon hypothesis of aging’ that postulates an increase in autonomous expression of such elements. Instead, a concerted increase in transposon expression, intron retention, and transcriptional readthrough points to a global increase in transcriptional dysfunction during aging and senescence, at least in human cells. Nevertheless, there are several other techniques such as immunoblotting for ORF1 protein, rapid amplification of cDNA ends (RACE), or genomic PCR that may provide evidence in favor of age-related reactivation of potentially autonomous transposons, specifically LINE-1 elements (De Cecco et al., 2019).
Furthermore, since most of our analysis focused on cellular models of aging and senescence, we cannot exclude that RNA-seq data performed on tissues samples would more faithfully and specifically capture the expression of autonomous transposable elements.
Although we have shown increased transcriptional dysfunction with aging and cellular senescence, so far we lack a model to explain these findings. In this context, it is noteworthy that the elongation rate of Pol II is correlated with readthrough transcription (Fong et al., 2015), intron processing, and may increase with aging (Debès et al., 2023), potentially providing a unifying explanation. Another potential explanation that could account for both elevated readthrough and expression of intergenic transposon is loss of repressive heterochromatin (Zhang et al., 2021) during aging or senescence. We were unable to test this hypothesis directly since we only had access to RNA-seq data for these samples.
Methods
Dataset selection, read preprocessing and alignment
We searched the SRA for suitable datasets that used total RNA extraction because this improves the detection of non-canonical transcripts that may be degraded before successful polyadenylation. Publicly available human RNA-seq datasets of aging and senescence, respectively, were obtained from Colombo et al., 2018 (GSE60340) and Fleischer et al., 2018 (GSE113957). Human RNA-seq and TT-TimeLapse-seq (TT-TL-seq) data were obtained from Rosa-Mercado et al., 2021 (GSE152059), human native elongating transcript sequencing (mNET-seq) data from Bauer et al., 2018 (PRJNA432639), and mouse RNA-seq data from Vilborg et al., 2017 (GSE98906).
Mouse liver RNA-seq data from ad libitum and dietary restricted mice were obtained from Hahn et al., 2017 (GSE92486) and mouse replicative senescence data from Wang et al., 2022 (GSE179880).
Reads were filtered, adaptor-trimmed and repaired using fastp (Chen et al., 2018) with a min Phred score of 25. Afterwards reads passing these filters were aligned to the human genome (GRCh38) or mouse genome (GRCm39) using STAR 2.7.8a (Dobin et al., 2013) in single-pass mode. To optimize the quantification of transposons the STAR parameters --outFilterMultimapNmax and --winAnchorMultimapNmax were set to 100. The full pipeline with all the configuration files is available on github (https://github.com/pabisk/aging_transposons copy archived at Pabis, 2024).
Detection of readthrough and intron retention
We employed a modified Automatic Readthrough Transcription Detection (ARTDeco) pipeline, using their provided annotations (Roth et al., 2020) but with a different quantification approach. Readthrough candidate regions in the ARTDeco pipeline were defined based on the longest transcript to exclude the effects of differential 3′-UTR usage. Genic and readthrough reads were counted using featureCounts from the Rsubread package in 10 kb bins downstream of genes, in a window between 0 and 100 kb. Readthrough regions with low coverage, in close proximity to, or overlapping, nearby protein-coding genes were excluded from the analysis. Intron retention was detected using the intron REtention Analysis and Detector (iREAD), which is a tool that detects retention events based on both splice junction reads and intron expression levels (Li et al., 2020). All samples were treated as unstranded.
Locus-specific detection of transposon reads and transcripts
We used TELocal v1.1.1 (mhammell-laboratory/TElocal) to quantify transposon reads, which is functionally similar to TETranscripts (Jin et al., 2015) except it provides locus-specific expression levels. Both tools use an expectation maximization algorithm to assign ambiguous reads. Annotation files were downloaded from the mhammell lab and from genecode (v39 for human and M27 for mouse). TELocal was also used to quantify transcriptomic counts. Again, all samples were treated as unstranded.
Expression and normalization of genes, introns, transposons, and readthrough transcripts
Transcripts and elements with low expression counts were filtered. After filtering, intron, readthrough, read-in, or transposon elements were combined with transcript counts for the analysis in DESeq2. For the aging fibroblast dataset age was normalized and centered and the design model considered age + gender + ethnicity. For the induced senescence dataset all four inducers were included together and the remaining treatment groups were considered non-senescent.
In order to study expression of these elements unbiased by gene co-expression, we selected the top 1000 overexpressed elements and divided their read counts by the read counts of the nearest gene. When not specified in the annotation the nearest genes or elements were identified with the bedtools closest or intersect function (Quinlan and Hall, 2010), as appropriate. Finally, to provide a single expression value for each sample we used the mean log of the normalized expression values.
Expression of potentially active and protein-coding LINE-1 elements
We defined potentially active LINE-1 elements as those that have a functional promoter (McKerrow and Fenyö, 2020). Whereas to define potentially protein-coding LINE-1 elements we predicted open reading frames for each annotated element using the predORF function from the systemPipeR package. Global pairwise alignment to ORF1p (uniprot: Q9UN81) and ORF2p (uniprot: O00370) was performed and only elements with a significant alignment were included for further analyses (p < 0.05).
Data availability
The current manuscript is a computational study, so no data have been generated for this manuscript.
-
NCBI Gene Expression OmnibusID GSE60340. Pathway Profiling of Replicative and Induced Senescence.
-
NCBI Gene Expression OmnibusID GSE113957. Predicting age from the transcriptome of human dermal fibroblasts.
-
NCBI Gene Expression OmnibusID GSE152059. TT-TL-seq reveals transcriptional profiles that accompany DoG induction after hyperosmotic stress.
-
NCBI BioProjectID PRJNA257655. Meromictic Lake Metagenome.
-
NCBI Gene Expression OmnibusID GSE98906. Comparative analysis reveals genomic features of stress-induced transcriptional readthrough.
-
NCBI Gene Expression OmnibusID GSE92486. Dietary restriction protects from age-associated DNA methylation and induces epigenetic reprogramming of lipid metabolism.
-
NCBI Gene Expression OmnibusID GSE179880. Unveiling E2F4, TEAD1 and AP-1 as regulatory transcription factors of the replicative senescence program by multi-omics analysis.
References
-
fastp: an ultra-fast all-in-one FASTQ preprocessorBioinformatics 34:i884–i890.https://doi.org/10.1093/bioinformatics/bty560
-
A comprehensive approach to expression of L1 lociNucleic Acids Research 45:e31.https://doi.org/10.1093/nar/gkw1067
-
Regulation of transposable elements by DNA modificationsNature Reviews. Genetics 20:417–431.https://doi.org/10.1038/s41576-019-0106-6
-
STAR: ultrafast universal RNA-seq alignerBioinformatics 29:15–21.https://doi.org/10.1093/bioinformatics/bts635
-
New Insights into Chronological Mobility of Retrotransposons In VivoOxidative Medicine and Cellular Longevity 2019:2818415.https://doi.org/10.1155/2019/2818415
-
Mobile DNA in health and diseaseThe New England Journal of Medicine 377:361–370.https://doi.org/10.1056/NEJMra1510092
-
On the wrong DNA track: Molecular mechanisms of repeat-mediated genome instabilityThe Journal of Biological Chemistry 295:4134–4170.https://doi.org/10.1074/jbc.REV119.007678
-
Death of a dogma: eukaryotic mRNAs can code for more than one proteinNucleic Acids Research 44:14–23.https://doi.org/10.1093/nar/gkv1218
-
Triplex and other DNA motifs show motif-specific associations with mitochondrial DNA deletions and species lifespanMechanisms of Ageing and Development 194:111429.https://doi.org/10.1016/j.mad.2021.111429
-
SoftwareAging_Transposons, version swh:1:rev:408f6c28826dee3dcb6895a2ce562fe036de372dSoftware Heritage.
-
Who let the DoGs out? - biogenesis of stress-induced readthrough transcriptsTrends in Biochemical Sciences 47:206–217.https://doi.org/10.1016/j.tibs.2021.08.003
-
ARTDeco: automatic readthrough transcription detectionBMC Bioinformatics 21:214.https://doi.org/10.1186/s12859-020-03551-0
-
Organ-, sex- and age-dependent patterns of endogenous L1 mRNA expression at a single locus resolutionNucleic Acids Research 49:5813–5831.https://doi.org/10.1093/nar/gkab369
-
Healthy aging interventions reduce repetitive element transcriptsThe Journals of Gerontology. Series A, Biological Sciences and Medical Sciences 76:805–810.https://doi.org/10.1093/gerona/glaa302
-
Intron Retention as a Mode for RNA-Seq data analysisFrontiers in Genetics 11:586.https://doi.org/10.3389/fgene.2020.00586
Article and author information
Author details
Funding
National University of Singapore (NUHSRO/2020/114)
- Brian K Kennedy
The funders had no role in study design, data collection, and interpretation, or the decision to submit the work for publication.
Acknowledgements
We would like to thank Juliane Liepe, Marco Malavolta, and Chin-Tong Ong for support and/or their help in improving this manuscript and VitaDAO for financial support.
Version history
- Preprint posted:
- Sent for peer review:
- Reviewed Preprint version 1:
- Reviewed Preprint version 2:
- Version of Record published:
Cite all versions
You can cite all versions using the DOI https://doi.org/10.7554/eLife.87811. This DOI represents all versions, and will always resolve to the latest one.
Copyright
© 2023, Pabis et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 1,616
- views
-
- 139
- downloads
-
- 2
- 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
-
- Chromosomes and Gene Expression
- Computational and Systems Biology
Genes are often regulated by multiple enhancers. It is poorly understood how the individual enhancer activities are combined to control promoter activity. Anecdotal evidence has shown that enhancers can combine sub-additively, additively, synergistically, or redundantly. However, it is not clear which of these modes are more frequent in mammalian genomes. Here, we systematically tested how pairs of enhancers activate promoters using a three-way combinatorial reporter assay in mouse embryonic stem cells. By assaying about 69,000 enhancer-enhancer-promoter combinations we found that enhancer pairs generally combine near-additively. This behaviour was conserved across seven developmental promoters tested. Surprisingly, these promoters scale the enhancer signals in a non-linear manner that depends on promoter strength. A housekeeping promoter showed an overall different response to enhancer pairs, and a smaller dynamic range. Thus, our data indicate that enhancers mostly act additively, but promoters transform their collective effect non-linearly.
-
- Chromosomes and Gene Expression
- Medicine
LncRNAs are involved in modulating the individual risk and the severity of progression in metabolic dysfunction-associated fatty liver disease (MASLD), but their precise roles remain largely unknown. This study aimed to investigate the role of lncRNA Snhg3 in the development and progression of MASLD, along with the underlying mechanisms. The result showed that Snhg3 was significantly downregulated in the liver of high-fat diet-induced obesity (DIO) mice. Notably, palmitic acid promoted the expression of Snhg3 and overexpression of Snhg3 increased lipid accumulation in primary hepatocytes. Furthermore, hepatocyte-specific Snhg3 deficiency decreased body and liver weight, alleviated hepatic steatosis and promoted hepatic fatty acid metabolism in DIO mice, whereas overexpression induced the opposite effect. Mechanistically, Snhg3 promoted the expression, stability and nuclear localization of SND1 protein via interacting with SND1, thereby inducing K63-linked ubiquitination modification of SND1. Moreover, Snhg3 decreased the H3K27me3 level and induced SND1-mediated chromatin loose remodeling, thus reducing H3K27me3 enrichment at the Pparg promoter and enhancing PPARγ expression. The administration of PPARγ antagonist T0070907 improved Snhg3-aggravated hepatic steatosis. Our study revealed a new signaling pathway, Snhg3/SND1/H3K27me3/PPARγ, responsible for mice MASLD and indicates that lncRNA-mediated epigenetic modification has a crucial role in the pathology of MASLD.