Increased PTCHD4 expression via m6A modification of PTCHD4 mRNA promotes senescent cell survival.

RNA modifications, including N6-methyladenosine (m6A), critically modulate protein expression programs in a range of cellular processes. Although the transcriptomes of cells undergoing senescence are strongly regulated, the landscape and impact of m6A modifications during senescence are poorly understood. Here, we report a robust m6A modification of PTCHD4 mRNA, encoding Patched Domain-Containing Protein 4, in senescent cells. The METTL3/METTL14 complex was found to incorporate the m6A modification on PTCHD4 mRNA; addition of m6A rendered PTCHD4 mRNA more stable and increased PTCHD4 production. MeRIP RT-qPCR and eCLIP analyses were used to map this m6A modification to the last exon of PTCHD4 mRNA. Further investigation identified IGF2BP1, but not other m6A readers, as responsible for the stabilization and increased abundance of m6A-modified PTCHD4 mRNA. Silencing PTCHD4, a transmembrane protein, enhanced growth arrest and DNA damage in pre-senescent cells and sensitized them to senolysis and apoptosis. Our results indicate that m6A modification of PTCHD4 mRNA increases the production of PTCHD4, a protein associated with senescent cell survival, supporting the notion that regulating m6A modification on specific mRNAs could be exploited to eliminate senescent cells for therapeutic benefit.

Senescence suppresses the integrated stress response and activates a stress-enhanced secretory phenotype.

Senescence is a state of indefinite cell cycle arrest associated with aging, cancer, and age-related diseases. Here, using label-based mass spectrometry, ribosome profiling and nanopore direct RNA sequencing, we explore the coordinated interaction of translational and transcriptional programs of human cellular senescence. We find that translational deregulation and a corresponding maladaptive integrated stress response (ISR) is a hallmark of senescence that desensitizes senescent cells to stress. We present evidence that senescent cells maintain high levels of eIF2α phosphorylation, typical of ISR activation, but translationally repress production of the stress response transcription factor 4 (ATF4) by ineffective bypass of the inhibitory upstream open reading frames. Surprisingly, ATF4 translation remains inhibited even after acute proteotoxic and amino acid starvation stressors, resulting in a highly diminished stress response. Furthermore, absent a response, stress augments the senescence secretory phenotype, thus intensifying a proinflammatory state that exacerbates disease. Our results reveal a novel mechanism that senescent cells exploit to evade an adaptive stress response and remain viable.

DPP4 inhibition impairs senohemostasis to improve plaque stability in atherosclerotic mice.

Senescent vascular smooth muscle cells (VSMCs) accumulate in the vasculature with age and tissue damage and secrete factors that promote atherosclerotic plaque vulnerability and disease. Here, we report increased levels and activity of dipeptidyl peptidase 4 (DPP4), a serine protease, in senescent VSMCs. Analysis of the conditioned media from senescent VSMCs revealed a unique senescence-associated secretory phenotype (SASP) signature comprising many complement and coagulation factors; silencing or inhibiting DPP4 reduced these factors and increased cell death. Serum samples from persons with high risk for cardiovascular disease contained high levels of DPP4-regulated complement and coagulation factors. Importantly, DPP4 inhibition reduced senescent cell burden and coagulation and improved plaque stability, while single-cell resolution of senescent VSMCs reflected the senomorphic and senolytic effects of DPP4 inhibition in murine atherosclerosis. We propose that DPP4-regulated factors could be exploited therapeutically to reduce senescent cell function, reverse senohemostasis, and improve vascular disease.

Pleiotropic effects of BAFF on the senescence-associated secretome and growth arrest.

Senescent cells release a variety of cytokines, proteases, and growth factors collectively known as the senescence-associated secretory phenotype (SASP). Sustained SASP contributes to a pattern of chronic inflammation associated with aging and implicated in many age-related diseases. Here, we investigated the expression and function of the immunomodulatory cytokine BAFF (B-cell activating factor; encoded by the TNFSF13B gene), a SASP protein, in multiple senescence models. We first characterized BAFF production across different senescence paradigms, including senescent human diploid fibroblasts (WI-38, IMR-90) and monocytic leukemia cells (THP-1), and tissues of mice induced to undergo senescence. We then identified IRF1 (interferon regulatory factor 1) as a transcription factor required for promoting TNFSF13B mRNA transcription in senescence. We discovered that suppressing BAFF production decreased the senescent phenotype of both fibroblasts and monocyte-like cells, reducing IL6 secretion and SA-ß-Gal staining. Importantly, however, the influence of BAFF on the senescence program was cell type-specific: in monocytes, BAFF promoted the early activation of NF-κB and general SASP secretion, while in fibroblasts, BAFF contributed to the production and function of TP53 (p53). We propose that BAFF is elevated across senescence models and is a potential target for senotherapy.

Single-cell analysis of skeletal muscle macrophages reveals age-associated functional subpopulations.

Tissue-resident macrophages represent a group of highly responsive innate immune cells that acquire diverse functions by polarizing toward distinct subpopulations. The subpopulations of macrophages that reside in skeletal muscle (SKM) and their changes during aging are poorly characterized. By single-cell transcriptomic analysis with unsupervised clustering, we found 11 distinct macrophage clusters in male mouse SKM with enriched gene expression programs linked to reparative, proinflammatory, phagocytic, proliferative, and senescence-associated functions. Using a complementary classification, membrane markers LYVE1 and MHCII identified four macrophage subgroups: LYVE1-/MHCIIhi (M1-like, classically activated), LYVE1+/MHCIIlo (M2-like, alternatively activated), and two new subgroups, LYVE1+/MHCIIhi and LYVE1-/MHCIIlo. Notably, one new subgroup, LYVE1+/MHCIIhi, had traits of both M2 and M1 macrophages, while the other new subgroup, LYVE1-/MHCIIlo, displayed strong phagocytic capacity. Flow cytometric analysis validated the presence of the four macrophage subgroups in SKM and found that LYVE1- macrophages were more abundant than LYVE1+ macrophages in old SKM. A striking increase in proinflammatory markers (S100a8 and S100a9 mRNAs) and senescence-related markers (Gpnmb and Spp1 mRNAs) was evident in macrophage clusters from older mice. In sum, we have identified dynamically polarized SKM macrophages and propose that specific macrophage subpopulations contribute to the proinflammatory and senescent traits of old SKM.

Early SRC activation skews cell fate from apoptosis to senescence.

Cells responding to DNA damage implement complex adaptive programs that often culminate in one of two distinct outcomes: apoptosis or senescence. To systematically identify factors driving each response, we analyzed human IMR-90 fibroblasts exposed to increasing doses of the genotoxin etoposide and identified SRC as a key kinase contributing early to this dichotomous decision. SRC was activated by low but not high levels of etoposide. With low DNA damage, SRC-mediated activation of p38 critically promoted expression of cell survival and senescence proteins, while SRC-mediated repression of p53 prevented a rise in proapoptotic proteins. With high DNA damage, failure to activate SRC led to elevation of p53, inhibition of p38, and apoptosis. In mice exposed to DNA damage, pharmacologic inhibition of SRC prevented the accumulation of senescent cells in tissues. We propose that inhibiting SRC could be exploited to favor apoptosis over senescence in tissues to improve health outcomes.

AUF1 ligand circPCNX reduces cell proliferation by competing with p21 mRNA to increase p21 production.

Mammalian circRNAs can influence different cellular processes by interacting with proteins and other nucleic acids. Here, we used ribonucleoprotein immunoprecipitation (RIP) analysis to identify systematically the circRNAs associated with the cancer-related protein AUF1. Among the circRNAs interacting with AUF1 in HeLa (human cervical carcinoma) cells, we focused on hsa_circ_0032434 (circPCNX), an abundant target of AUF1. Overexpression of circPCNX specifically interfered with the binding of AUF1 to p21 (CDKN1A) mRNA, thereby promoting p21 mRNA stability and elevating the production of p21, a major inhibitor of cell proliferation. Conversely, silencing circPCNX increased AUF1 binding to p21 mRNA, reducing p21 production and promoting cell division. Importantly, eliminating the AUF1-binding region of circPCNX abrogated the rise in p21 levels and rescued proliferation. Therefore, we propose that the interaction of circPCNX with AUF1 selectively prevents AUF1 binding to p21 mRNA, leading to enhanced p21 mRNA stability and p21 protein production, thereby suppressing cell growth.

Mitochondrial RNA in Alzheimer's Disease Circulating Extracellular Vesicles.

Alzheimer's disease (AD) is the most common type of dementia. Amyloid ß (Aß) plaques, tau-containing neurofibrillary tangles, and neuronal loss leading to brain atrophy are pathologic hallmarks of AD. Given the importance of early diagnosis, extensive efforts have been undertaken to identify diagnostic and prognostic biomarkers for AD. Circulating extracellular vesicles (EVs) provide a platform for "liquid biopsy" biomarkers for AD. Here, we characterized the RNA contents of plasma EVs of age-matched individuals who were cognitively normal (healthy controls (HC)) or had mild cognitive impairment (MCI) due to AD or had mild AD dementia (AD). Using RNA sequencing analysis, we found that mitochondrial (mt)-RNAs, including MT-ND1-6 mRNAs and other protein-coding and non-coding mt-RNAs, were strikingly elevated in plasma EVs of MCI and AD individuals compared with HC. EVs secreted from cultured astrocytes, microglia, and neurons after exposure to toxic conditions relevant to AD pathogenesis (Aß aggregates and H2O2), contained mitochondrial structures (detected by electron microscopy) and mitochondrial RNA and protein. We propose that in the AD brain, toxicity-causing mitochondrial damage results in the packaging of mitochondrial components for export in EVs and further propose that mt-RNAs in plasma EVs can be diagnostic and prognostic biomarkers for MCI and AD.

Evolutionarily Selected Overexpression of the Cytokine BAFF Enhances Mucosal Immune Response Against P. falciparum.

We have previously shown that a variant of the TNFSF13B gene that we called BAFF-var increases the production of the cytokine BAFF, upregulating humoral immunity and increasing the risk for certain autoimmune diseases. In addition, genetic population signatures revealed that BAFF-var was evolutionarily advantageous, most likely by increasing resistance to malaria infection, which is a prime candidate for selective pressure. To evaluate whether the increased soluble BAFF (sBAFF) production confers protection, we experimentally assessed the role of BAFF-var in response to malaria antigens. Lysates of erythrocytes infected with Plasmodium falciparum (iRBCs) or left uninfected (uRBCs, control) were used to treat peripheral blood mononuclear cells (PBMCs) with distinct BAFF genotypes. The PBMCs purified from BAFF-var donors and treated with iRBCs showed different levels of specific cells, immunoglobulins, and cytokines as compared with BAFF-WT. In particular, a relevant differential effect on mucosal immunity B subpopulations have been observed. These findings point to specific immune cells and molecules through which the evolutionary selected BAFF-var may have improved fitness during P. falciparum infection.

Survey of senescent cell markers with age in human tissues.

Cellular senescence, triggered by sublethal damage, is characterized by indefinite growth arrest, altered gene expression patterns, and a senescence-associated secretory phenotype. While the accumulation of senescent cells during aging decreases tissue function and promotes many age-related diseases, at present there is no universal marker to detect senescent cells in tissues. Cyclin-dependent kinase inhibitors 2A (p16/CDKN2A) and 1A (p21/CDKN1A) can identify senescent cells, but few studies have examined the numbers of cells expressing these markers in different organs as a function of age. Here, we investigated systematically p16- and p21-positive cells in tissue arrays designed to include normal organs from persons across a broad spectrum of ages. Increased numbers of p21-positive and p16-positive cells with donor age were found in skin (epidermis), pancreas, and kidney, while p16-expressing cells increased in brain cortex, liver, spleen and intestine (colon), and p21-expressing cells increased in skin (dermis). The numbers of cells expressing p16 or p21 in lung did not change with age, and muscle did not appear to have p21- or p16-positive cells. In summary, different organs display different levels of the senescent proteins p16 and p21 as a function of age across the human life span.

Transcriptome signature of cellular senescence.

Cellular senescence, an integral component of aging and cancer, arises in response to diverse triggers, including telomere attrition, macromolecular damage and signaling from activated oncogenes. At present, senescent cells are identified by the combined presence of multiple traits, such as senescence-associated protein expression and secretion, DNA damage and ß-galactosidase activity; unfortunately, these traits are neither exclusively nor universally present in senescent cells. To identify robust shared markers of senescence, we have performed RNA-sequencing analysis across eight diverse models of senescence triggered in human diploid fibroblasts (WI-38, IMR-90) and endothelial cells (HUVEC, HAEC) by replicative exhaustion, exposure to ionizing radiation or doxorubicin, and expression of the oncogene HRASG12V. The intersection of the altered transcriptomes revealed 50 RNAs consistently elevated and 18 RNAs consistently reduced across all senescence models, including many protein-coding mRNAs and some non-coding RNAs. We propose that these shared transcriptome profiles will enable the identification of senescent cells in vivo, the investigation of their roles in aging and malignancy and the development of strategies to target senescent cells therapeutically.

Loss of miR-451a enhances SPARC production during myogenesis.

MicroRNAs (miRNAs) are small noncoding RNAs that critically regulate gene expression. Their abundance and function have been linked to a range of physiologic and pathologic processes. In aged monkey muscle, miR-451a and miR-144-3p were far more abundant than in young monkey muscle. This observation led us to hypothesize that miR-451a and miR-144-3p may influence muscle homeostasis. To test if these conserved microRNAs were implicated in myogenesis, we investigated their function in the mouse myoblast line C2C12. The levels of both microRNAs declined with myogenesis; however, only overexpression of miR-451a, but not miR-144-3p, robustly impeded C2C12 differentiation, suggesting an inhibitory role for miR-451a in myogenesis. Further investigation of the regulatory influence of miR-451a identified as one of the major targets Sparc mRNA, which encodes a secreted protein acidic and rich in cysteine (SPARC) that functions in wound healing and cellular differentiation. In mouse myoblasts, miR-451a suppressed Sparc mRNA translation. Together, our findings indicate that miR-451a is downregulated in differentiated myoblasts and suggest that it decreases C2C12 differentiation at least in part by suppressing SPARC biosynthesis.

NF90 regulation of immune factor expression in response to malaria antigens.

Nuclear factor 90 (NF90) is a dual DNA- and RNA-binding protein expressed ubiquitously in mammalian cells, including monocytes. Here, to elucidate the function of NF90 in the immune response, we analyzed systematically its influence on gene expression programs in the human monocytic cell line THP-1 expressing normal or reduced NF90 levels. RNA sequencing analysis revealed many mRNAs showing differential abundance in NF90-silenced cells, many of them encoding proteins implicated in the response to immune stimuli and malaria infection. The transcription of some of them (e.g. TNF, LILRB1, and CCL2 mRNAs) was modulated by silencing NF90. Ribonucleoprotein immunoprecipitation (RIP) analysis further revealed that a subset of these mRNAs associated directly with NF90. To understand how NF90 influenced globally the immune response to malaria infection, lysates of red blood cells infected with Plasmodium falciparum (iRBC lysates) or uninfected/mock-infected (uRBC lysates) were used to treat THP-1 cells as a surrogate of malaria infection. NF90 affected the stability of a few target mRNAs, but influenced more generally the translation and secretion of the encoded cytokines after treatment with either uRBC or iRBC lysates. Taken together, these results indicate that NF90 contributes to repressing the immune response in cells responding to P. falciparum infection and suggest that NF90 can be a therapeutic target in malaria.

Loss of RNA-binding protein GRSF1 activates mTOR to elicit a proinflammatory transcriptional program.

The RNA-binding protein GRSF1 (G-rich RNA sequence-binding factor 1) critically maintains mitochondrial homeostasis. Accordingly, loss of GRSF1 impaired mitochondrial respiration and increased the levels of reactive oxygen species (ROS), triggering DNA damage, growth suppression, and a senescent phenotype characterized by elevated production and secretion of interleukin (IL)6. Here, we characterize the pathways that govern IL6 production in response to mitochondrial dysfunction in GRSF1-depleted cells. We report that loss of GRSF1 broadly altered protein expression programs, impairing the function of respiratory complexes I and IV. The rise in oxidative stress led to increased DNA damage and activation of mTOR, which in turn activated NF-κB to induce IL6 gene transcription and orchestrate a pro-inflammatory program. Collectively, our results indicate that GRSF1 helps preserve mitochondrial homeostasis, in turn preventing oxidative DNA damage and the activation of mTOR and NF-κB, and suppressing a transcriptional pro-inflammatory program leading to increased IL6 production.

Cooperative translational control of polymorphic BAFF by NF90 and miR-15a.

Polymorphisms in untranslated regions (UTRs) of disease-associated mRNAs can alter protein production. We recently identified a genetic variant in the 3'UTR of the TNFSF13B gene, encoding the cytokine BAFF (B-cell-activating factor), that generates an alternative polyadenylation site yielding a shorter, more actively translated variant, BAFF-var mRNA. Accordingly, individuals bearing the TNFSF13B variant had higher circulating BAFF and elevated risk of developing autoimmune diseases. Here, we investigated the molecular mechanisms controlling the enhanced translation of BAFF-var mRNA. We identified nuclear factor 90 (NF90, also known as ILF3) as an RNA-binding protein that bound preferentially the wild-type (BAFF-WT mRNA) but not BAFF-var mRNA in human monocytic leukemia THP-1 cells. NF90 selectively suppressed BAFF translation by recruiting miR-15a to the 3'UTR of BAFF-WT mRNA. Our results uncover a paradigm whereby an autoimmunity-causing BAFF polymorphism prevents NF90-mediated recruitment of microRNAs to suppress BAFF translation, raising the levels of disease-associated BAFF.

Regulation of HuR structure and function by dihydrotanshinone-I.

The Human antigen R protein (HuR) is an RNA-binding protein that recognizes U/AU-rich elements in diverse RNAs through two RNA-recognition motifs, RRM1 and RRM2, and post-transcriptionally regulates the fate of target RNAs. The natural product dihydrotanshinone-I (DHTS) prevents the association of HuR and target RNAs in vitro and in cultured cells by interfering with the binding of HuR to RNA. Here, we report the structural determinants of the interaction between DHTS and HuR and the impact of DHTS on HuR binding to target mRNAs transcriptome-wide. NMR titration and Molecular Dynamics simulation identified the residues within RRM1 and RRM2 responsible for the interaction between DHTS and HuR. RNA Electromobility Shifts and Alpha Screen Assays showed that DHTS interacts with HuR through the same binding regions as target RNAs, stabilizing HuR in a locked conformation that hampers RNA binding competitively. HuR ribonucleoprotein immunoprecipitation followed by microarray (RIP-chip) analysis showed that DHTS treatment of HeLa cells paradoxically enriched HuR binding to mRNAs with longer 3'UTR and with higher density of U/AU-rich elements, suggesting that DHTS inhibits the association of HuR to weaker target mRNAs. In vivo, DHTS potently inhibited xenograft tumor growth in a HuR-dependent model without systemic toxicity.

LncRNA OIP5-AS1/cyrano suppresses GAK expression to control mitosis.

Some long noncoding RNAs (lncRNAs) can regulate gene expression programs, in turn affecting specific cellular processes. We sought to identify the mechanism through which the lncRNA OIP5-AS1, which is abundant in the cytoplasm, suppressed cell proliferation. Silencing of OIP5-AS1 in human cervical carcinoma HeLa cells triggered the appearance of many aberrant (monopolar, multipolar, misaligned) mitotic spindles. Through a combination of approaches to pull down mRNAs bound to OIP5-AS1 and identify proteins differentially expressed when OIP5-AS1 was silenced, we identified a subset of human cell cycle regulatory proteins encoded by mRNAs that interacted with OIP5-AS1 in HeLa cells. Further analysis revealed that GAK mRNA, which encodes a cyclin G-associated kinase important for mitotic progression, associated prominently with OIP5-AS1. The interaction between these two transcripts led to a reduction in GAK mRNA stability and GAK protein abundance, as determined in cells in which OIP5-AS1 levels were increased or decreased. Importantly, the aberrant mitotic cell division seen after silencing OIP5-AS1 was partly rescued if GAK was simultaneously silenced. These findings indicate that the abnormal mitoses seen after silencing OIP5-AS1 were caused by an untimely rise in GAK levels and suggest that OIP5-AS1 suppresses cell proliferation at least in part by reducing GAK levels.

High-purity circular RNA isolation method (RPAD) reveals vast collection of intronic circRNAs.

High-throughput RNA sequencing methods coupled with specialized bioinformatic analyses have recently uncovered tens of thousands of unique circular (circ)RNAs, but their complete sequences, genes of origin and functions are largely unknown. Given that circRNAs lack free ends and are thus relatively stable, their association with microRNAs (miRNAs) and RNA-binding proteins (RBPs) can influence gene expression programs. While exoribonuclease treatment is widely used to degrade linear RNAs and enrich circRNAs in RNA samples, it does not efficiently eliminate all linear RNAs. Here, we describe a novel method for the isolation of highly pure circRNA populations involving RNase R treatment followed by Polyadenylation and poly(A)+ RNA Depletion (RPAD), which removes linear RNA to near completion. High-throughput sequencing of RNA prepared using RPAD from human cervical carcinoma HeLa cells and mouse C2C12 myoblasts led to two surprising discoveries: (i) many exonic circRNA (EcircRNA) isoforms share an identical backsplice sequence but have different body sizes and sequences, and (ii) thousands of novel intronic circular RNAs (IcircRNAs) are expressed in cells. In sum, isolating high-purity circRNAs using the RPAD method can enable quantitative and qualitative analyses of circRNA types and sequence composition, paving the way for the elucidation of circRNA functions.

Identification of HuR target circular RNAs uncovers suppression of PABPN1 translation by CircPABPN1.

HuR influences gene expression programs and hence cellular phenotypes by binding to hundreds of coding and noncoding linear RNAs. However, whether HuR binds to circular RNAs (circRNAs) and impacts on their function is unknown. Here, we have identified en masse circRNAs binding HuR in human cervical carcinoma HeLa cells. One of the most prominent HuR target circRNAs was hsa_circ_0031288, renamed CircPABPN1 as it arises from the PABPN1 pre-mRNA. Further analysis revealed that HuR did not influence CircPABPN1 abundance; interestingly, however, high levels of CircPABPN1 suppressed HuR binding to PABPN1 mRNA. Evaluation of PABPN1 mRNA polysomes indicated that PABPN1 translation was modulated positively by HuR and hence negatively by CircPABPN1. We propose that the extensive binding of CircPABPN1 to HuR prevents HuR binding to PABPN1 mRNA and lowers PABPN1 translation, providing the first example of competition between a circRNA and its cognate mRNA for an RBP that affects translation.

HuR and GRSF1 modulate the nuclear export and mitochondrial localization of the lncRNA RMRP.

Some mitochondrial long noncoding RNAs (lncRNAs) are encoded by nuclear DNA, but the mechanisms that mediate their transport to mitochondria are poorly characterized. Using affinity RNA pull-down followed by mass spectrometry analysis, we found two RNA-binding proteins (RBPs), HuR (human antigen R) and GRSF1 (G-rich RNA sequence-binding factor 1), that associated with the nuclear DNA-encoded lncRNA RMRP and mobilized it to mitochondria. In cultured human cells, HuR bound RMRP in the nucleus and mediated its CRM1 (chromosome region maintenance 1)-dependent export to the cytosol. After RMRP was imported into mitochondria, GRSF1 bound RMRP and increased its abundance in the matrix. Loss of GRSF1 lowered the mitochondrial levels of RMRP, in turn suppressing oxygen consumption rates and modestly reducing mitochondrial DNA replication priming. Our findings indicate that RBPs HuR and GRSF1 govern the cytoplasmic and mitochondrial localization of the lncRNA RMRP, which is encoded by nuclear DNA but has key functions in mitochondria.