Healthy elderly individuals tend to have lower heart rate, increased heart rate variability, and phase advanced activity patterns relative to young subjects.

Experimental protocol for enrolment and monitoring of human subjects (A). As assessed by Zyphyr BioPatch, EKG measurements suggest that MESOR of average heart rate decreases with age (p=0.056) and that MESOR of heart rate variability increases with age (p=0.056). N=5 per group. MESOR and amplitude were tested using two-sided Wilcoxon rank sum exact test, while phase was tested by Kuiper’s two-sample test (B). Average activity counts across three axes (vector magnitude), as recorded by the actigraph device plotted throughout the day (top) and analyzed for circadian rhythm (bottom) (C). While the amplitudes of activity did not differ between the age groups, the older individuals trended towards an early phase (p∼0.055, Kupier’s two-sample test) compared to young individuals. N=7 for young and N=8 for old. Lines in the top panel of C are smoothed means (fit with penalized cubic regression splines) for data from each age group. Dots in the bottom panel of C are subject-level cosinor parameter estimates derived from cosinor fits to the actigaraphy data. Summary Figure (A) was made with BioRender.com. Boxplot midlines correspond to median values, while the lower and upper hinges correspond to the first and their quartiles, respectively. Boxplot whiskers extend to the smallest/largest points within 1.5 * IQR (Inter Quartile Range) of the lower/upper hinge.

© 2024, BioRender Inc. Any parts of this image created with BioRender are not made available under the same license as the Reviewed Preprint, and are © 2024, BioRender Inc.

The circadian transcriptome is differentially affected by entrainment with sera from young and old human subjects.

A visual representation of the serum starvation-serum addition protocol to synchronize the BMAL1-luciferase BJ-5TA fibroblasts (A, made with BioRender.com). When comparing the circadian transcriptome entrained by either young or old sera (n=4 sera per group), 1519 genes lost rhythmicity with age (B, top) while only 637 genes gained rhythmicity with age (B, bottom). Weighted BIC criterion with a threshold of 0.75 was used to assess rhythmicity (B). The number of genes rhythmic in young and old (according to weighted BIC > 0.75 criterion) that show a detectable change in MESOR (q < 0.05 criterion for MESOR difference) is 568, with MESOR increasing in 163 and decreasing in 405 genes. Out of these, 40 genes with increased MESOR (red) and 98 genes with decreased MESOR (blue) also satisfy the condition |log₂ FC| > 0.25 (C). We were only able to detect change in amplitude for a small number of genes, 39 genes had decreased amplitude in old and 2 had increased amplitude using CircaCompare. Only 30 genes with decreased amplitude also satisfy the condition |log₂ FC| > 0.1 (blue) (D). For phase, using a test provided by CircaCompare, under q < 0.05 cutoff for age-related phase differences in genes rhythmic in young and old (with BIC>0.75 cycling criteria), we detected 20 genes with advanced phase, and 34 with delayed phase (E).

© 2024, BioRender Inc. Any parts of this image created with BioRender are not made available under the same license as the Reviewed Preprint, and are © 2024, BioRender Inc.

Each type of circadian change is associated with different KEGG pathways by STRING analysis

Entrainment of the fibroblasts in culture with old serum significantly altered the circadian transcriptome despite use of the same cells in both young and old conditions, suggesting the aged serum itself affects the regulation of specific pathways in the cell. Significant KEGG Pathways (FDR <1x10-4 unless specified) are indicated next to each category of genes. Age related pathways such as Alzheimer’s Disease/oxidative phosphorylation are associated with a loss of rhythms in the old condition. Cell cycle and DNA replication pathways remain rhythmic in the old serum condition, but cycle with a decreased MESOR. Rhythmic genes were determined using CircaCompare. Rhythmicity was determined using the weighted BIC > 0.75 criterion while p-values for difference in MESOR and phase were determined using CircaCompare.

Synchronizing with old serum phase delays the expression profile of several core clock genes

Traces of molecular clock mRNA transcripts, the average curve (bold) of results of individual sera (faded). P-values for the difference in MESOR (P_#x1D759;M), amplitude (P_#x1D759;A), and phase (P_#x1D759;P) are shown for the comparison of young and old conditions. Several clock genes are significantly phase delayed (CRY1, CRY2, NR1D1, NR1D2, PER3, TEF) in response to synchronization with old serum. While BMAL1 is not phase delayed, the MESOR significantly increases with age. N=4 subjects per timepoint for both young and old groups.

There are no significant differences in grip strength, sympathetic, parapsympathetic nervous system indices, or cortisol levels between young and old individuals

Sympathetic nervous system (SNS, A, top) and parasympathetic nervous system (PNS, A, bottom) indices were measured using EKG. Amplitude, MESOR, and Phase were calculated from the resulting oscillations and the Median, quartiles, and SEMs are shown. Dots represent individuals. P>0.05. N=5 per group. MESOR and amplitude were tested using Wilcoxon rank sum exact test, while phase was tested by Kuiper’s two-sample test. Boxplot midlines correspond to median values, while the lower and upper hinges correspond to the first and third quartiles, respectively. Boxplot whiskers extend to the smallest/largest points within 1.5 * IQR (inter Quartile Range) of the lower/upper hinge (A). Serum cortisol levels are not statistically different by unpaired t-test (B, N=7 young, 8 old).

Demographic information for subjects in the study.

Young and old serum contain equivalent levels of cortisol and are equally effective at entraining cells in culture

(A) Cells synchronized with human serum from either young (N=7 subjects) or old (N=8 subjects) individuals did not show differences in amplitude (as measured in luminescence counts per second, bottom, left), period (bottom, middle), or phase (bottom, right). Data points represent average values for a specific patient’s serum over 2 experiments with 2-3 replicates per experiment. A visual representation of the individual replicates averaged in A (B) Summary statistics are displayed as mean +/- SEM. Means compared by unpaired t-test. Despite individual and run to run variability, cells synchronized with human serum from either young (N=7 subjects) or old (N=8 subjects) individuals showed similar amplitude (top, left), period (top, right), and phase (bottom) of BMAL1-luciferase rhythms. Data points represent individual replicates within an experiment. Wells run in the same experiment are displayed in the same color. Each subject’s sample was run in 2-3 replicates.

The circadian transcriptomes of cells synchronized with young or old sera deviate significantly on day two.

On day two of serum entrainment, young and old transcriptional rhythms differentiated. CircaCompare analysis of RNA sequencing revealed that on Day 2 of serum entrainment analysis (CT 36-58) the MESOR differences of cycling genes are larger between the cells entrained with young or old serum (A, top). Additionally, more genes were phase shifted in the old serum condition compared to young serum on Day 2 (B, top). The distribution of p-values shows an enrichment of low p-values on Day 2 for both MESOR and phase differences (A, bottom, B, bottom). The size of each circle in the scatter plots is proportional to -log10q, so the bigger the circle, the smaller the q-value for the difference between young and old for each metric.

Oxidative phosphorylation genes and chromosomal replication genes are affected by entrainment with old serum

Additional pathway analysis employing IPA further emphasizes finding with STRING analysis. Genes identified by IPA analysis to be associated with oxidative phosphorylation/mitochondrial dysfunction and eIF2 genes lose rhythmicity with age (BIC >0.75) (A, B, C). Cell cycle checkpoint control (CHK proteins) (D) and chromosomal replication genes (E) maintain their cycling in the aged condition, which is consistent with the continued division of cells.

eIF2 signaling pathway genes that maintain cycling with age phase advance in the old serum condition

From the IPA analyses, the eIF2 signaling pathway was enriched (BIC >0.75) in two datasets: 1) genes cycling in the young sera and not the old and 2) cycling genes phase-advanced in the old sera (q<0.05). A significant number of genes that maintain rhythmicity in the old sera condition are phase advanced as visualized in the heat map (A) or plot of peak phase (B). As the eIF2 pathway is involved in protein translation, this suggests that the timing of protein translation is phase advanced in the older serum and/or less rhythmic.

Genes with decreased MESOR in the IPA cell cycle control of chromosomal replication pathway

Table of genes that are rhythmic in both conditions with a decreased MESOR with old serum treatment compared to young serum treatment by CircaCompare. All genes are involved in the cell cycle/DNA replication pathway.

Genes with increased MESOR in the IPA cholesterol biosynthesis pathway

Table of genes that are rhythmic in both conditions with a decreased MESOR with old serum treatment compared to young serum treatment by CircaCompare. All genes are involved in the cholesterol biosynthesis pathway.