Supplementary Figures
Epigenetic clock and methylation studies in the rhesus macaque
Steve Horvath1,2#, Joseph A. Zoller 2#, Amin Haghani1, Anna J. Jasinska 3, Ken Raj4, Charles E. Breeze5, Jason Ernst6, Kelli L. Vaughan7, Julie A. Mattison7
Correspondence:
1) Steve Horvath; E-mail: shorvath@mednet.ucla.edu; Address: Gonda Building, 695 Charles Young Drive South, Los Angeles, CA 90095
2) Julie A. Mattison; Email: Julie.mattison@nih.gov; Translational Gerontology Branch, National Institute on Aging, NIH.
Supplementary Figures
179185172214102104 13 46 7 37 70 58137 2
135206231123125224286107109 67 68 49 48 47 66 69124203283 44
194121199221 5 30 19 24 42 52 78 53 17 34 72 73 20
118213 12 39 64 86 71 87 18 40 26 75 27 82 90
205 14
136126225 25 36 62 89149 88 3 74 15 59 8222230119220 91 31 60
100202 4
232130132235146148134196 98209110228128201 1
15079 22122 77
120208105212233
43 97197198103210284 16147227 38 81
139 50 8011651 92 28 23 33 41193143145195115117 29 6144 32 94
113229106236101112237 11 55
129218 57 96 21 93 95142223127234287226200138140 61 65 83
114215211131133217111216285108
99 204219141207 9 84 54 63 45 56 76 85 166 10 35
156151152153154155175169187160181174180162186239258270191259248266244250275265267280241247253271272246277278256262269238240281263279264282257274252260276242249261251245268243255254273288 189161173167168 171165177183158190 192178159184 176157182163170164188
0.000.100.200.30
Cluster Dendrogram
hclust (*, "average") as.dist(1 - corSample)
Height
branch Tissue CanBeUsedForAging Age Female
Supplementary Figure 1. Unsupervised hierarchical clustering of tissue samples.
Average linkage hierarchical clustering based on the interarray correlation coefficient (Pearson correlation). A height cut-off of 0.05 led to branch colors that largely correspond to Tissue type (second panel).
The second color band encodes tissue type: green=liver, red=lung, blue=blood, adipose=turquoise, brown=cerebral cortex, yellow=kidney, black=muscle, pink=skin. A handful of arrays are severe outliers (left- and right most clusters) as indicated by the third color band (turquoise samples are severe outliers).
These technical outliers probably result from insufficient amounts of DNA. These outlying samples were removed from the analysis.
1
Supplementary Figure 2: Pan-tissue clock for rhesus macaque applied to different tissues. Each panel correlates the chronological age (x-axis) at the time of tissue collection with the leave-one-out (LOO) estimate of DNA methylation age. A) All tissues. Dots are colored by tissue type as indicated in the remaining panels. Results for B) adipose, C) blood, D) brain cortex, E) kidney, F) liver, G) lung, H) muscle, I) skin. Each panel reports the sample size (N), Pearson correlation, and the median absolute error, i.e., the median of the absolute difference between DNAmAgeLOO and chronological age. The dashed line is the diagonal y=x. The solid line corresponds to the least-squares regression line.
2
Supplementary Figure 3. EWAS was performed in each macaque tissue separately using the R function
"standardScreeningNumericTrait" from the "WGCNA" R package. Z statistics from correlation tests in different tissues. Upper panels report scatter plots for Z statistics in different tissues. Lower panels:
corresponding Pearson correlation coefficients.
3