circRNAs expressed in human peripheral blood are associated with human aging phenotypes, cellular senescence and mouse lifespan

Circular RNAs (circRNAs) are an emerging class of non-coding RNA molecules that are thought to regulate gene expression and human disease. Despite the observation that circRNAs are known to accumulate in older organisms and have been reported in cellular senescence, their role in aging remains relatively unexplored. Here, we have assessed circRNA expression in aging human blood and followed up age-associated circRNA in relation to human aging phenotypes, mammalian longevity as measured by mouse median strain lifespan and cellular senescence in four different primary human cell types. We found that circRNAs circDEF6, circEP300, circFOXO3 and circFNDC3B demonstrate associations with parental longevity or hand grip strength in 306 subjects from the InCHIANTI study of aging, and furthermore, circFOXO3 and circEP300 also demonstrate differential expression in one or more human senescent cell types. Finally, four circRNAs tested showed evidence of conservation in mouse. Expression levels of one of these, circPlekhm1, was nominally associated with lifespan. These data suggest that circRNA may represent a novel class of regulatory RNA involved in the determination of aging phenotypes, which may show future promise as both biomarkers and future therapeutic targets for age-related disease. Electronic supplementary material The online version of this article (10.1007/s11357-019-00120-z) contains supplementary material, which is available to authorized users.


Introduction
Aging is a multifactorial process leading to gradual deterioration of physical and physiological functionality at the cellular, tissue and organ levels. It is the primary risk factor for chronic aging pathologies more exons. They have been proposed to be key regulators of gene expression by various mechanisms including sequestration of RNA-binding proteins and miRNAs or by acting as a competitor of linear splicing of their cognate genes (Memczak et al. 2013). The possibility that a single circRNA could sequester several such RNA regulators suggests that this class of non-coding RNAs could modulate many cellular and physiological processes through multiple pathways. circRNAs are known to accumulate in older organisms (Gruner et al. 2016), and some have been reported to be implicated in cellular senescence (Du et al. 2017;Du et al. 2016). Despite these promising findings, their role in aging remains relatively unexplored.
We hypothesized that expression of some circRNAs may be associated with advancing age, aging phenotypes, lifespan or cellular senescence. Changes in circRNA expression over a 5-year period were assessed in relation to age, combined parental longevity score (PLS) and hand grip strength. We then assessed expression levels of 15 circRNAs in early passage and late passage primary human dermal fibroblasts, cardiomyocytes, astrocytes and vascular endothelial cells. Finally, the junction sequences of relevant exons were examined for conservation between mouse and humans and where evidence was present that the back-spliced junction, and thus, the circular RNA were conserved; we assessed expression in relation to longevity in six strains of mice with differential median strain longevities.
We present here evidence that although effects on age itself did not replicate in the wider sample set, the expression levels of circEP300 (β = − 0.065, P = 0.001) and circFOXO3 (β = − 0.060, P = 0.002) were negatively associated with parental longevity score. circDEF6 was positively associated with parental longevity score (β = 0.070, P = 0.024) although this did not reach multiple testing thresholds. circFNDC3B was also nominally associated with hand grip strength (β = 0.004, P = 0.039). circRNAs (7/12 (58%)) expressed in senescent human primary astrocytes, endothelial cells, fibroblasts or cardiomyocytes also demonstrated dysregulated expression in one or more cell types. Comparative sequence analysis suggested that four circRNAs may be conserved in mice. When assessed, circPlekhm1 transcript level in spleen was also demonstrated to be positively associated with mouse median strain lifespan (β = 0.0025; P = 0.017). These results suggest that some age-related circRNAs may play roles in molecular drivers of aging such as cellular senescence, and hence may represent potential contributors to lifespan or other human aging phenotypes.

InCHIANTI cohort and selection of participants
The InCHIANTI study of Aging is a population study of aging (Ferrucci et al. 2000). Participants undertook detailed assessment of health and lifestyle parameters at baseline, and again at three subsequent follow-ups (FU2 2004-2006, FU3 2007-2010and FU4 2012-2014. The present study used participants from the third and fourth follow-up visits (FU3 and FU4). RNA samples and clinical/phenotypic data were already available for 698 participants at FU3. The collection of the FU4 samples and data comprise part of this study. During the FU4 interviews in 2012/2013, blood and clinical/phenotypic data were collected from 455 study participants. These data were cross-checked against RNA samples and clinical/phenotypic data already held from FU3, to ensure that sample and phenotypic data was available from both collections. Sample-associated data included measures of potential confounding factors such as BMI, sex, level of education (none, elementary, secondary, high school and university), study site, smoking and white blood counts (neutrophil, lymphocyte, monocyte, eosinophil percentages). Characteristics of the study population are given in Table 1. Informed consent was obtained from all participants. Ethical approval was obtained from the Instituto Nazionale Riposo e Cura Anziani institutional review board, Italy.
Generation of circRNA profiles from old and young human peripheral blood Circular RNA profiles were initially generated in parallel from two sets of pooled peripheral blood total RNA samples using a modified 'CircleSeq' procedure (Lopez-Jimenez et al. 2018). 2 μg RNA (RNA integrity number (RIN) = 6.4) was assessed in two separate pools from 20 'young' samples (median age = 33 years, range 30-36 years, 55% female, 45% male; RIN 5.6) and 20 'old' samples (median age 87 years, range 86-95 years, 90% female 10% male, RIN 7.7). Each pooled sample was divided into two aliquots, one of which was treated with 20 units RNAse R (Epicentre, Madison, USA) at 30°C for 30 min to remove linear RNA, the other sample being mock-treated using 1 μL RNase-free water in place of the enzyme. Both aliquots were cleaned and concentrated using 2 volumes of RNA clean beads (Beckman Coulter, Indianapolis, USA) to remove the enzyme. The results of the RNase R treatment were confirmed on a highsensitivity RNA screentape (Agilent, Santa Clara, USA). Ribosomal RNA was removed, and indexed sequencing libraries made using the libraries were determined by qPCR and adjusted for size using Tapestation D1000 analysis (Agilent, Santa Clara, USA). Ribosomal RNA was removed, and indexed sequencing libraries made using the Illumina RNASeq protocol. The library concentrations were determined by qPCR and adjusted for size using the data from the Tapestation D1000 analysis. Libraries were pooled in equimolar quantities, denatured and diluted to 12.0 pM + 1% PhiX for clustering and then underwent 125 paired-end Illumina sequencing in four lanes using TruSeq SBS reagents (V3).

Analysis of circRNA profiles
RNase R and mock-treated sequence data were assembled, and putative circular RNAs were identified using PTESFinder (Izuogu et al. 2016) with the human genome (hg19) reference files provided with the software, a segment size of 65 and a uniqueness score of 7. The remaining parameters were left to default settings. To calculate a comparable measure of circular RNA abundance between samples, we used a measure termed back-spliced reads per million mapped reads (bpm) for each circular RNA defined as where J i is the number of reads mapped to the backspliced junction of the circular RNA, c is the number of reads mapped to canonical sites of the gene with the circular RNA and n is the number of circular RNAs identified. This measure is designed to be similar to the commonly used reads per kilobase per million mapped reads (RPKM) metric used regularly to estimate gene expression from RNA-Seq data. In addition to circular RNA detection using PTESFinder, reads from all samples were also mapped to the human genome reference (hg19) obtained from iGenomes using Tophat v2.1.0 with the pre-set sensitive alignment parameters in paired-end mode (Trapnell et al. 2009). The number of reads mapping to each exon of each gene was then calculated using FeatureCounts v2.0.0 with parameters for unstranded alignment, paired reads, count multimapping reads and assigning reads to overlapping features (Liao et al. 2013;Liao et al. 2014). Counts were used to calculate RPKM per exon using the standard method to compare the expression of each exon across samples.
Pathway analysis of differentially regulated circRNA host genes circRNAs showing expression differences between the pooled old and the pooled young samples were ranked by RPKM and fold change. To assess whether circRNAs demonstrating expression differences between young and old pools were enriched in genes derived from specific molecular or biochemical function groups, we carried out a Cytoscape version 2.5.2 plug-in ClueGO analysis. This platform queries over-representation of query genes in specific KEGG, REACTOME and WikiPathways (Bindea et al. 2009). The linear genes hosting the top 10% most abundantly expressed circRNAs in young and old pools for the circRNA profile were queried against KEGG_20.11.2017, REACTOME_Pathways_20.11.2017 and WikiPathways_20.11.2017. Outputs were selected based on 'enrichment/depletion' through a two-sided hypergeometric test with Bonferroni step down for P value correction with the selected ontology reference set of chosen genes. The GO terms were used to group functional pathways, and the leading functional grouping was based on highest significant kappa score.
Design of qPCR assays for circRNA validation Levels of individual circRNA in young and old pools were ranked by abundance. circRNAs demonstrating evidence of altered expression with age fell into three classes: those expressed exclusively in old, those expressed exclusively in young, and those expressed in both young and old, but with evidence that levels were different between the pools. We selected five circRNAs exclusively expressed in young (circITGAX, circPLEKHM1, circDEF6, circATP6V0A1 and circASAP1), five exclusively expressed in the old (circFOXO3, circFNDC3B, circAFF1, circCDYL and circXPO7), as well as five expressed in both pools but demonstrating evidence of altered expression (circMIB1, circMETTL3, circBCL11B, circZC3H18 and circEP300), where sequence and assay design constraints allowed for to design specific assays to unique back-spliced junction for qRTPCR follow-up.
circRNA probe design Custom-designed qRTPCR assays for quantification of relative expression were designed to unique backspliced circRNA junctions (Thermo Fisher, Foster City, USA), the sequences of which are given in Online Resource 1. Each target sequence was checked for the presence of single nucleotide polymorphisms in potential primer or probe binding regions prior to ordering. Assays were ordered as custom single-tube assays from Thermo Fisher (Foster City, USA). Each circRNA probe was validated using standard curve analysis using 1:10 serial dilutions of synthetic oligonucleotides homologous to the back-spliced junctions.
Assessment of associations between circRNA expression and aging phenotypes in the InCHIANTI cohort RNA samples and phenotypic data were available from 306 individuals at both follow-up 3 (FU3) and followup 4 (FU4) of the InCHIANTI study of aging. Characteristics of participants are given in Table 1. We assessed the expression of 15 age-associated circRNAs demonstrating the most marked differential expression with age between young and old pools as described above. Aging parameters assessed were age itself, parental longevity score (PLS) and hand grip strength. Participants aged 65 + years were categorised for PLS based on the age at death of their parents. Short, intermediate and long-lived cut-offs were calculated separately for mothers and fathers based on the normal distribution of age at death in the cohort, as described in Dutta et al. (2013a). Mothers and fathers aged < 49 years or < 52 years at death respectively were classed as premature and excluded. To standardize parental age of death, a Z score was generated for combined maternal and paternal measures of parental longevity. Hand-grip strength was measured in kilograms using a dynamometer, with repeated measurements at both FU3 and FU4.
Reverse transcription and pre-amplification of circRNAs in human peripheral blood RNA cDNA synthesis was carried out using 100 ng total RNA using the High-Capacity cDNA Reverse Transcription Kit (Thermo Fisher, Foster City, USA) according to manufacturer's instructions (Fisher Scientific, New Hampshire, USA) in a final reaction volume of 10.0 μL per sample. Reactions (samples in 96-well plates) were run at 25°C for 10 min, 37°C for 120 min, 85°C for 5 min followed by an inactivation period for 95°C for 10 min. Pre-amplification of circRNA expression was carried out using 5 μL TaqMan PreAmp master mix (Thermo Fisher, Foster City, USA), 2.5 μL pooled assay mix and 2.5 μL cDNA in a final reaction volume of 10 μL per sample. Cycling conditions were one cycle of 95°C for 10 min followed by 14 cycles of 95°C for 15 s with 60°C for 4 min followed by 95°C for 10 min. Pre-amplified samples were then diluted 1:10 and maintained on ice prior to analysis.

Assessment of associations between circRNA expression in peripheral blood RNA and human aging phenotypes
The expression profiles of selected circRNAs were then measured in total peripheral blood mRNA using customdesigned OpenArray plates on the Thermo Fisher 12K Flex platform (Thermo Fisher, Foster City, USA

Assessment of circRNA conservation between mouse and human
We assessed whether the 15 circRNAs identified in our human study were likely to be conserved in mouse by aligning the mouse and human exon junction sequences using the Blat tool in the UCSC genome browser (https://genome.ucsc.edu). Quantitative real-time PCR assays were developed to unique back-spliced junctions of conserved circRNAs. Probe and primer sequences are given in Online Resource 2. circRNA expression was then measured in mouse spleen and muscle tissue and assessed in relation to lifespan by analysis of levels in six strains of male mice (A/J, NOD.B10Sn-H2 b /J, PWD/PhJ, 129S1 /SvlmJ, C57BL/6J and WSB/EiJ) selected on the basis of divergent median strain longevity (Yuan et al. 2009). Animal husbandry, handling, animal characteristics and sample preparation protocols have been previously described (Lee et al. 2016). Tissue samples were obtained from cross-sectional study conducted in the same compartment and in the same period of time as described in Yuan et al. (2009). Spleen and quadricep muscle tissues were excised immediately after sacrifice and shipped from the Jackson Laboratory using RNAlater-ICE Collection protocol (Life Technologies, Carlsbad, CA). In this method, tissues are submerged in RNAlater stabilization solution; an aqueous tissue storage reagent used to rapidly permeate tissues and stabilize RNA from fresh specimens and stored at -20°C or below for later use.

RNA extraction and reverse transcription from mouse tissues
Total RNA was extracted using the TRI Reagent/ chloroform phase separation according to manufacturer's instructions. Briefly, tissues stored in RNA later were drained, and then placed in 1 mL TRI Reagent solution containing 10 mM MgCl 2 . Samples were homogenized for 15 min (spleen) or 30 min (muscle) using bead mills (Retsch Technology GmbH, Haan, Germany). This was followed by a phase separation using chloroform. Total RNAs in the separated RNAs were precipitated from the aqueous phase through overnight incubation with isopropanol at − 20°C. The following morning, RNA pellets were washed twice with ethanol and resuspended in RNase-free dH 2 O. Complementary DNA (cDNA) was generated from 100 ng RNA using the Evocript Universal cDNA Master Synthesis kit according to the manufacturer's instructions (Roche, Switzerland).
Assessment of circRNA expression in mouse spleen and muscle circRNAs selected on the basis of interspecies sequence conservation were validated in mouse spleen and muscle tissue. Expression levels of conserved circRNAs were assessed in relation to median strain lifespan by relative quantification. Quantitative qRTPCR was carried out for circRNAs (circFoxo3, circMib1, circPlekhm1 and circXpo7) in relation to the Pol2ra, Trfc and Ipo8 endogenous control genes, selected on the basis of lack of age association in a previous study (Harries et al. 2011). Reaction mixes contained cDNA (0.5 μL), TaqMan Universal PCR mastermix II (2.5 μL, no AmpErase UNG (Thermo Fisher, Foster City, USA), dH 2 O (1.75 μL, Fisher Scientific, USA), and TaqMan gene assay (0.25 μL, Thermo Fisher, Foster City, USA) in a 5 μL final reaction volume. The reaction mixes were centrifuged at 3000 rpm, vortexed and centrifuged again at 3000 rpm and transferred to 384-well qRTPCR plates. qRTPCR was run at 50°C for 2 min, 95°C for 10 min and 50 cycles of 15 s at 95°C for 30 s and 1 min at 60°C. Each sample assay was conducted in three technical triplicates. Expression levels of circRNAs in young and old mouse tissues were measured relative to the geometric mean of the entire set of transcripts, with the expression of each individual circRNA normalised to the global mean of expression of each circRNA across the samples. Linear regression analysis was carried out to assess the association of expression of circRNA using StataSE15 (StataCorp, TX, USA).

circRNA profile in peripheral blood of aging humans
One hundred sixty-six to 167M reads were obtained from the RNAse R-treated pools and 157-163M reads from the mock-treated pools with a mean Q score of 34.6-35.1 and total error rate of 0.47-0.53%. A total of 2207 circRNAs were expressed in human peripheral blood. Of these, 184 circRNAs were found in both the young and old samples, 431 were exclusively expressed in the young sample pool and 1592 were exclusively expressed in the old sample pool (Online Resource 3). We selected 15 circRNAs for further analysis: 5 expressed exclusively in the young pool, 5 expressed exclusively in the old pool and 5 expressed in both pools but showing the most discrepant expression for further study. These were circITGAX, circPLEKHM1, circDEF6, circATP6V0A1 and circASAP1 which showed exclusive expression in the young; circFOXO3, circFNDC3B, circAFF1, circCDYL and circXPO7 which showed exclusive expression in the old; and circMIB1, circMETTL3, circEP300, circZC3H18 and circBCL11B that were expressed, but differentially so in both sample pools.

Pathway analysis of circRNA expressed in aging humans
Pathway enrichment for the genes hosting the top 10% most abundant circRNAs in each of young and old pooled peripheral blood samples was performed using ClueGO cytoscape (Bindea et al. 2009). In the young peripheral blood, the top 10% most abundant circRNAs derived from genes associated with negative regulation of ATP metabolic processes and in transmission of synaptic signals. The leading edge genes hosting circRNAs for negative regulation of ATP processes were SNCA, STAT3 and UFSP2, whilst those associated with synaptic vesicle endocytosis were FCH02, PICALM, PIP5K1C and SNCA. Genes hosting circRNAs were primarily localised in pathways involved in phagocytosis, circadian regulation, cancer pathways and golgi-associated vesicle budding in the blood from aged donors (Table 2).
circPLEKHM1, circMETTL and circFNDC3B expression levels are associated with aging phenotypes in humans The structures of the 15 circRNAs selected for follow-up were predicted based on the sequencing read depth for each exon and are presented in Fig.  1. Exon structures presented as read depth plots are given in Online Resource 4. Although we demonstrated no associations with age itself, we did identify associations between some circRNAs and human aging phenotypes. circEP300 and circFOXO3 both demonstrated negative associations with combined parental longevity score (β = − 0.065 and − 0.060; P = 0.001 and 0.002 respectively), after adjustment for multiple testing. circDEF6 was positively correlated with parental longevity scores but demonstrated nominal significance only (β = 0.070, P = 0.024) ( Table 3, Fig. 2). A positive association was also identified both cross-sectionally (β = 0.004, P = 0.039) and longitudinally (β = 0.004, P = 0.038) between circFNDC3B expression and hand grip strength (Table 4, Fig. 3), although these were nominal only.
circRNAs are differentially expressed in early passage and late passage cells Twelve of 15 circRNAs tested were expressed in astrocytes, endothelial cells, fibroblasts or astrocytes. Seven (58%) of these demonstrated differential expression between early and late passage cells of one or more cell type (Table 5). circAFF1 and circFOXO3 demonstrated associations in more than one cell type although direction of effect was concordant only for circFOXO3 (in cardiomyocytes and fibroblasts). circCDYL, circEP300, circMIB1, circZC3H18 and circMETTL3 were differentially expressed in only one cell type. circBCL11B, circDEF6 and circITGAX were not expressed in any cell type tested.

Differential expression of circRNAs between mice of different median strain longevities
In silico analyses suggested that four circRNAs (circFoxo3, circMib1, circPlekhm1 and circXpo7) may have conserved back-spliced junction in the mouse. Associations with longevity were then assessed in spleen and muscle tissue from young (6 months) and old (20-22 months) mouse strains of six different median strain longevities. circMib1 and circXpo7 were expressed only in spleen, whereas circFoxo3 and circPlekhm1 were expressed in both tissues ( Table 6). The expression of circPlekhm1 demonstrated a nominal positive correlation with median lifespan in young and old (β = 0.0013, P = 0.016) as well as in spleen of young mice (β = 0.0025, P = 0.017), although these were not significant after adjustment

Discussion
Circular RNAs (circRNAs) are an emerging class of regulatory RNA molecule thought to play a role in human disease (Haque and Harries 2017). These molecules have no free ends, and as such are exonuclease resistant. circRNAs accumulate in aged organisms (Gruner et al. 2016) and have been suggested to play a role in cellular senescence (Du et al. 2017;Du et al. 2016). We hypothesised that the human circRNAome may differ in aged humans compared with younger subjects and that these changes may also be associated with cellular senescence or with longevity in animal models. We identified > 2000 circRNAs in total RNA from human blood, some of which were expressed exclusively in samples from older donors. GSEA pathways enrichment analysis of genes hosting the top 10% most abundant circRNAs in elderly donors suggested that pathways involved in phagocytosis, circadian regulation, cancer pathways and golgi-associated vesicles were the most enriched in these genes. We demonstrated that three circRNAs (circDEF6, circFOXO3 and circEP300) were associated with measures of parental longevity, and one (circFNDC3B) was associated with hand grip strength both longitudinally and cross-sectionally. Furthermore, 7 of 12 circRNAs expressed in human senescent cells of different cell types demonstrated dysregulated expression in one or more cell type and 1 of 4 circRNAs demonstrating conserved expression were associated with median strain longevity in spleen tissue from young mice. These findings are consistent with the hypothesis that some circRNAs have roles in molecular aging and the determination of mammalian aging phenotypes. circRNAs generated from the FOXO3 and EP300 genes were negatively associated with measures of human parental longevity and also demonstrated dysregulated expression in human senescent cells. circRNAs deriving from the FOXO3 gene have previously been demonstrated to regulate cell cycle when manipulated by gene knockdown in mouse embryonic fibroblasts, cardiac fibroblasts or mammary cancer cell lines (Du Table 3  Beta coefficients, p values and 95% confidence intervals (95% CI) are given for associations between circRNAs expression and combined parental longevity (PLS) score. Two hundred ninety-one samples were assessed. Genes demonstrating statistically significant results below the multiple testing limit of 0.003 are indicated in italics, whilst those demonstrating nominal associations only are given in bold type Fig. 2 circRNA expression is associated with combined parental longevity. Forest plot illustrating the association between peripheral blood circRNA expression and combined human parental longevity score (PLS) in participants from the InCHIANTI study of aging. N = 306 individuals. The beta-coefficient of the association is given on the X-axis, and the identity of the gene is given on the Y-axis. Lines attached to each data point represent 95% confidence intervals (95% CI). Statistical significance is indicated by stars, *< 0.05, **< 0.005 et al. 2016). Furthermore, FOXO3 circular RNAs also demonstrate elevated expression and association with cellular senescence in the heart tissue of mice and humans (Du et al. 2017). It is not clear whether the previously reported circular FOXO3 transcripts have the same structure as the one we have identified, since previous studies do not give its exon structure. A circRNA from the FOXO3 gene identical to the one we have identified has also previously been demonstrated to inhibit myoblast differentiation in mouse cells (Li et al. 2019). Genetic variation in the FOXO3 gene itself has previously been associated with extreme longevity (Flachsbart et al. 2017;Fuku et al. 2016) and has also been associated with maintenance of telomere length (Davy et al. 2018). circRNAs deriving from the EP300 gene have not been previously reported. EP300 encodes the repressor histone acetyltransferase protein p300, which also has roles as a transcriptional corepressor protein. EP300 has been implicated in modulation of FOXO3 activity (Mahmud et al. 2019) (Jeung et al. 2016). Inhibition of EP300 has been shown to mimic calorific restriction in human and mouse cells (Pietrocola et al. 2018); calorific restriction is of course a well-known modifier of lifespan in many species (Austad 1989;Hansen et al. 2008;Kapahi et al. 2004;Mitchell et al. 2010). This protein is also a master regulator of autophagy, which is a pivotal factor in stem cell maintenance and evasion of cellular senescence (Vijayakumar and Cho 2019). circFNDC3B was positively associated with hand grip strength. Although these associations were nominal only, they were present both cross-sectionally and longitudinally. An average person may lose~20-40% of skeletal muscle mass as well as muscle strength from by the time they reach 80 years of age (Carmeli et al. 2002;Doherty 2003) and decline in skeletal muscle strength is predictive of disability and mortality in humans (Giampaoli et al. 1999;Rantanen et al. 1999;Rantanen et al. 2012). Circular RNAs originating from this gene have been reported previously, and suggested to possess tumour suppressor activity (Liu et al. 2018).

and in
The results generated from our mouse data suggest that circPlekhm1, which was associated with median strain longevity, may drive longevity, rather than being consequential to it, since the associations are present in the spleen RNA of young mice alone. The Plekhm1 gene encodes a multivalent adaptor protein that integrates endocytic and autophagic pathways at the lysosome (McEwan and Dikic 2015). Its role in lifespan may therefore stem from moderation of lysosomal trafficking since lysosomes play a critical part in successful aging and longevity (Carmona-Gutierrez et al. 2016;Simonsen et al. 2007).
Our study has both strengths and weaknesses. It represents one of the first circRNA profiles in aging human peripheral blood and provides data not only population-level epidemiological evidence for a role in human aging phenotypes, or mammalian lifespan, but also in vitro evidence that some circRNA may influence cell senescence phenotypes. Weaknesses individuals. The beta-coefficient of the association is given on the X-axis, and the identity of the gene is given on the Y-axis. Lines attached to each data point represent 95% confidence intervals (95% CI). Statistical significance is indicated by stars, *< 0.05, **< 0.005   circRNA expression is reported here in relation to median strain longevity. Data are assessed separately for young and old animals of each strain. N = 67 (muscle); 90 (spleen). Results reaching statistical significance are indicated in bold typeface IQR interquartile range, ND not detected include a relatively low power to detect effects of in the population study, which might be attributed to the biological variation in circRNA levels and limitations in samples size and power. Nevertheless, we were able to identify some interesting associations, which likely represent the largest effects. Future work could include validation of epidemiological data in larger sample sets and also functional delineation of the molecular effects of the circRNA in question. Our data provide evidence that circRNAs may play an important role in the determination of mammalian aging phenotypes. circRNAs are inherently stable, due to their exonuclease resistance, and are found not only in tissues relevant to human diseases, but also in the circulation, raising the possibility that they may prove useful as biomarkers of disease or targets for molecular therapies in the future.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http:// creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.