Epigenomic signature of accelerated ageing in progeroid Cockayne syndrome

Abstract Cockayne syndrome (CS) and UV‐sensitive syndrome (UVSS) are rare genetic disorders caused by mutation of the DNA repair and multifunctional CSA or CSB protein, but only CS patients display a progeroid and neurodegenerative phenotype, providing a unique conceptual and experimental paradigm. As DNA methylation (DNAm) remodelling is a major ageing marker, we performed genome‐wide analysis of DNAm of fibroblasts from healthy, UVSS and CS individuals. Differential analysis highlighted a CS‐specific epigenomic signature (progeroid‐related; not present in UVSS) enriched in three categories: developmental transcription factors, ion/neurotransmitter membrane transporters and synaptic neuro‐developmental genes. A large fraction of CS‐specific DNAm changes were associated with expression changes in CS samples, including in previously reported post‐mortem cerebella. The progeroid phenotype of CS was further supported by epigenomic hallmarks of ageing: the prediction of DNAm of repetitive elements suggested an hypomethylation of Alu sequences in CS, and the epigenetic clock returned a marked increase in CS biological age respect to healthy and UVSS cells. The epigenomic remodelling of accelerated ageing in CS displayed both commonalities and differences with other progeroid diseases and regular ageing. CS shared DNAm changes with normal ageing more than other progeroid diseases do, and included genes functionally validated for regular ageing. Collectively, our results support the existence of an epigenomic basis of accelerated ageing in CS and unveil new genes and pathways that are potentially associated with the progeroid/degenerative phenotype.


| INTRODUC TI ON
Progeroid syndromes are a group of rare genetic disorders in which ageing appears to be greatly accelerated (Burtner & Kennedy, 2010).
These conditions are genetically, biochemically, and clinically heterogeneous, with the extent and combination of ageing-related defects (Aryee et al., 2014) being specific to each progeroid syndrome.
While the Hutchinson-Gilford progeria syndrome (HGPS) is due to a defective nuclear lamina as a consequence of mutations in LMNA gene, the Down syndrome (DS) is due to trisomy of chromosome 21.The Werner syndrome (WS) and Cockayne syndrome (CS) result from mutation in multifunctional proteins best known for their role in DNA repair (Rieckher et al., 2021).
Here, we focus on CS which is a devastating multisystemic disease characterized by neuro-developmental defects, precocious ageing, and in most cases cutaneous photosensitivity.It is a progressive disorder with most symptoms appearing relatively early in life and worsening with time (Natale, 2011).Multiple degrees of severity have been described, including a prenatal-onset form known as cerebro-oculo-facioskeletal (COFS) syndrome, an early-onset severe form (CS-II, average age of death 5-6 years), a classical or moderate form (CS-I, mean life expectancy around 16 years), and a mild or lateonset form (CS-III) (Laugel et al., 2010).
Cockayne syndrome is due to mutations in either CSB or CSA genes (approximately 70:30 ratio) (Calmels et al., 2018;Laugel et al., 2010), which code repair factors of UV-induced DNA damage through the transcription-coupled nucleotide excision repair (TC-NER) pathway (Mayne & Lehmann, 1982).CS defects, including the precocious ageing phenotype, have been essentially ascribed to impaired TC-NER.Nevertheless, three patients carrying mutations in CSB or CSA developed an independent syndrome, the UV-sensitive syndrome (UVSS) that is characterized by severe photosensitivity but no sign of premature ageing (Horibata et al., 2004;Nardo et al., 2009).Only a few UVSS patients, which display a common phenotype, are known worldwide (Nardo et al., 2009).Most UVSS cases are due to mutations in UVSSA, a protein that acts in TC-NER in concert with CSA and CSB (Sarasin, 2012;van der Weegen et al., 2020).
These case reports, at least those mutated in CSB or CSA, de facto uncouple the DNA repair defect from precocious ageing and suggest that impairment of other functions of multifunctional CSA and CSB proteins are responsible for progeroid defects.
Indeed, CSA and CSB are also transcription factors (Kristensen et al., 2013), CSB is implicated in chromatin remodelling (Citterio et al., 2000;Newman et al., 2006), and have been detected in mitochondria where they act in repairing oxidative damage (Kamenisch et al., 2010) and promoting transcription (Berquist et al., 2012) of mitochondrial DNA.Additionally, we showed that in CS (but not UVSS) cells, the CSA/CSB impairment results in overexpression of the HTRA3 protease leading to degradation of the mitochondrial DNA polymerase POLG1, and consequent mitochondrial dysfunction (Chatre et al., 2015).CSA and CSB mutations represent therefore an exceptional case for disentangling downstream genes and pathways involved in dramatically different phenotypes, since their impairment results (CS) or not (UVSS) in accelerated ageing.We recently associated the CS-specific HTRA3/POLG1/mitochondrial alterations with replicative senescence in normal cells, and identified CSB depletion as an early trigger of this process, mechanistically linking progeroid defects with a major ageing-related process (Crochemore et al., 2019).
Different molecular hallmarks of ageing (e.g., epigenetic alterations, senescence and mitochondrial dysfunction) have been identified (Lopez-Otin et al., 2013).DNA methylation (DNAm), an epigenetic mechanism corresponding to the addition of a methyl group in the fifth position of a cytosine (Moore et al., 2013), is of particular relevance in the study of ageing and age-related conditions, as it substantially contributes to genomic stability and can regulate agedependent gene expression (Ciccarone et al., 2018).DNAm globally decreases with ageing, largely as a consequence of demethylation of repetitive retroelements (LINEs, Alu sequences), driving genomic instability (Ciccarone et al., 2018).At the same time, age-dependent hypermethylation or hypomethylation events have been reported for specific CpG sites located within genic sequences, although their functional consequences for gene expression changes are not straightforward (Yuan et al., 2015).Subsets of these specific sets of CpGs have been used to build epigenetic clocks, that is, predictors epigenetic clock returned a marked increase in CS biological age respect to healthy and UVSS cells.The epigenomic remodelling of accelerated ageing in CS displayed both commonalities and differences with other progeroid diseases and regular ageing.
CS shared DNAm changes with normal ageing more than other progeroid diseases do, and included genes functionally validated for regular ageing.Collectively, our results support the existence of an epigenomic basis of accelerated ageing in CS and unveil new genes and pathways that are potentially associated with the progeroid/degenerative phenotype.

K E Y W O R D S
ageing, Cockayne syndrome, DNA methylation, epigenetic clock, progeroid diseases, UVsensitive syndrome of chronological age that are also informative of biological age and mortality risk (Horvath, 2013;Levine et al., 2018;Lu et al., 2019).
Previous studies reported alterations of the DNAm profile in WS-, HGPS-and DS-derived cells compared to healthy controls (Guastafierro et al., 2017;Heyn et al., 2013;Horvath et al., 2015;Kohler et al., 2020) and discussed their possible contribution to the progeroid phenotype based on the comparison of changes occurring in physiological ageing.Interestingly, mild increased epigenetic age calculated using epigenetic clocks was reported for all the three conditions (Horvath et al., 2015;Horvath & Raj, 2018;Maierhofer et al., 2017).On the contrary, the epigenetic landscape of CS has not been described so far, and differences with non-progeroid UVSS have not been addressed.
In the present study we performed genome-wide analysis of the methylome of CS, UVSS and WT fibroblasts.We showed that large genome-wide DNAm changes occur in mutated cells (UVSS and CS).
In order to focus on the DNAm alterations that are associated with the age-accelerated phenotype of CS, we took advantage of the availability of fibroblasts from two very rare UVSS patients, that are mutated but non-progeroid.We thus included UVSS in the group of 'Non-Progeroid' controls, together with WT cells, and compared them with CS cells.This analysis identified CpG sites mapping in genes enriched in three major functional categories, and, importantly, the most significant of them display a methylation/transcription correlation according to CS transcriptomic datasets.We also demonstrated a dramatic increase of the epigenetic age of CS fibroblasts, and comparison with available datasets revealed a CS epigenetic signature closer to physiological ageing than HGPS, WS or DS, underscoring CS as a relevant model of accelerated ageing.In support of this finding, we report that a majority of top differentially methylated genes have been functionally validated in age-related processes.

| Identification of differentially methylated positions/regions (DMPs/DMRs) associated with the CS progeroid phenotype
We analysed genome-wide DNAm in dermal fibroblasts isolated from seven CS patients, two UVSS patients (including the only case mutated in CSA) and three healthy subjects (WT), with no siblings, and at similar and early passage number (PN 14).These fibroblasts are available because of single biopsies obtained for diagnostic purposes of the DNA repair defect, and are not repeatable longitudinally.An overview of the study design is reported in Figure S1.

Samples were processed using the Infinium HumanMethylation450
BeadChip or Infinium MethylationEPIC BeadChip microarrays, as indicated in Figure 1a, and seven samples were assessed using both platforms.After data pre-processing (see 5. Experimental procedures), we selected the probes that were in common between the two arrays (452,567 probes).In a first exploratory analysis, principal component analysis (PCA) was applied to all the probes of the microarray to identify the main sources of variability in genome-wide DNA methylation profiles.PCA showed a separation of samples according to the type of platform on which they had been processed (Figure S2A); within each batch, samples tended to cluster according to the mutational state (Figure S2B).After applying the ComBat algorithm to adjust for the platform effect, no differences between the 450 k and EPIC samples were observed (Figure S2C), and samples still clustered together according to the mutational state (Figure S2D).Along the first component, mutated cells (UVSS and CS) were separated respect to WT cells, with the largest differences exhibited between WT and CS-I cells (classical form).The second component identified a clear separation between CS-I cells mutated for CSA and CSB, which, to date, are clinically indistinguishable (Laugel et al., 2010).However, CS-I CSA and CS-I CSB derived from male and female patients, respectively.We therefore removed the probes mapping on X and Y chromosomes (10,585 probes) and performed again the PCA analysis, which returned a result substantially overlapping the previous one (Figure S2E-H).Distinct CS-I CSA and CS-I CSB clusters were confirmed also upon removal of 1184 autosomic probes whose methylation level has been shown to be affected by sex in whole blood (Singmann et al., 2015;  PCA results reported in the previous paragraph suggested the existence of large genome-wide DNA methylation differences between mutated and WT cells.Among the CpG sites showing DNAm changes, we aimed focusing on those specifically related to the progeroid phenotype that characterizes CS samples (independently from the subtype) and distinguishes them not only from WT, but also from UVSS samples, and thereby identify CS-specific DNA methylation changes.For this, we compared CpG DNAm between CS-derived fibroblasts ('Progeroid' group) versus WT and UVSS cells, merged in the same 'Non-Progeroid' group, thus prioritizing CS-specific epigenetic changes (progeroid-related) and understating those associated with the impaired response to UV, common to CS and UVSS (Figure 1d).
We first searched for differentially methylated positions (DMPs) between the non-progeroid and progeroid samples (Figure S1).Out of 441,982 analysed probes, we unveiled 11,597 DMPs (ANOVA; BH-corrected p-value < 0.05; absolute difference in mean DNAm values >10%; 2.6% of total probes), 9872 of which were hypomethylated and 1725 hypermethylated in the progeroid group (top-ranking hits in Figure 1e and Tables S1, S2 columns A-X).About two third of DMPs (7277) mapped in a genic sequence.Results of the Fisher's exact test to assess whether DMPs are over-or under-represented in specific genomic regions revealed that both hypo-and hypermethylated DMPs were under-represented (p-value < 0.05) in Islands (DNA segments of >200 bp in length, with high CpG density, >50% GC percentage, >60% observed/expected CpG ratio).Hypomethylated DMPs were strongly enriched in Open Sea (isolated CpGs in the genome) and Shelves (regions located within 2-4 kbp from CpG islands), whereas hypermethylated DMPs were enriched in Shores (regions located within 0-2 kbp from CpG islands), Shelves and Open Sea regions (Figure 1f).Unsupervised hierarchical clustering on the list of 11,597 DMPs confirmed a clear separation between the progeroid and non-progeroid groups (Figure 1g).Furthermore, CS subtypes tended to cluster in separate subgroups, such as WT and UVSS cells.
To strengthen the analysis of genomic regions undergoing DNA methylation changes specifically in CS, we searched for differentially methylated regions (DMRs).We focused on probes located on CpG Islands and surrounding regions (Shores and Shelves) of geneassociated regions, and applied a pipeline based on a Multivariate Analysis of Variance (MANOVA) of adjacent probes, which has been validated for relevant age-associated changes (Bacalini et al., 2015).
We identified 1817 DMRs in progeroid versus non-progeroid groups, of which 1416 were hypomethylated and 401 hypermethylated (MANOVA BH-corrected p-value < 0.05; Table S3).The top-ranking hits of DMRs list are reported in Table S4.These 1817 DMRs map in 1498 genes (1140 of which are hypomethylated and 358 hypermethylated), since some genes are covered by multiple DMRs.
Importantly, in several DMRs, UVSS samples displayed an intermediary DNAm profile between WT and CS samples (Figure S3A) rather than being close to WT (Figure S3B).We reasoned that some DNAm changes in this list could still be linked to impaired DNA repair (a feature common to UVSS and CS cells) rather than precocious ageing (only in CS).We thus filtered the DMRs list in order to retain only those that unambiguously separated progeroid (CS) samples and non-progeroid (WT and UVSS) samples, thus having the highest probability of being implicated in the mechanism leading to precocious ageing (see Supplementary information, "Differential methylation analysis").This approach resulted in a shorter and more stringent list of 222 hits that we named StringentDMRs (Figure S1 and Table S5), and contained 141 hypo-and 81 hypermethylated regions in progeroid compared to non-progeroid.Hierarchical clustering on the StringentDMRs list confirmed a clear separation between the progeroid and non-progeroid group, a clear clustering of CS subtypes, UVSS and WT cells in separate subgroups (Figure S3C), and the maintenance of top-ranking position of individual genes compared to the non-filtered DMRs list.
The forthy top-ranking hits of this list are reported in Table 1.For each of these genes, a literature survey was done to identify those reported to play a functional role in ageing.Interestingly, about half (11/20) of the top hypomethylated genes have been functionally implicated in ageing or ageing-related processes (e.g., senescence) by gain (overexpression) or loss (knock-out/down) of function experiments in human cells or animal models (Table S6).Other genes have not been tested or have not shown a role (e.g., ESYT3) in ageing.Similarly, 43% (9/21) of the top hypermethylated genes have a functional role in ageing, some with a higher number of reports such as RPTOR, HDAC4 and ATP132A (Table S6).The remaining top 12 hypermethylated genes have not been challenged in an ageing context.

| Functional enrichment analysis of positions and regions with CS-specific DNA methylation
To characterize the potential biological impact of CS-specific DNAm changes, we performed a gene set enrichment analysis (GSEA) on the genes associated with methylated positions (MPs) and methylated regions (MRs) that resulted from the progeroid vs non-progeroid comparison, using the Genome Ontology (GO) database.GSEA is a method that determines whether an a priori F I G U R E 1 Global features of DNA methylation in CS skin fibroblasts and profile of the differentially methylated positions in progeroid versus non-progeroid groups.(a) Characteristics of primary skin fibroblasts derived from three healthy donors (WT), two UVSS patients and seven CS patients.A plus indicates the platform used for the DNAm analysis (seven samples analysed with both platforms), and the presence of a progeroid phenotype in individuals.For each cell type is shown the age at skin biopsy, age at death/age at latest follow up.The star (*) indicates that the individual was alive at the indicated age (A.Sarasin, personal communication).All age columns are expressed in years.F, female; M, male; N, no; S, suspected; Y, yes.More information on these patients (except UVSTA24VI) is available in Laugel et al. (2010)  defined set of genes shows statistically significant concordant differences between two biological states (here progeroid versus non-progeroid).Since both MPs and MRs results displayed high semantic redundancy of significantly enriched GO terms (i.e., Adjusted p-value < 0.01), we removed the dispensable terms using REVIGO (Supek et al., 2011).The treemap visualization of the remaining terms (Biological processes), which allows hierarchical grouping of semantically related terms, identified 'superclusters' (meta-terms), prevalently genes implicated in embryonic morphogenesis/development, and also transport of ions, synaptic transmission and cell-cell adhesion, that were shared among MPs (Figure S4A) and MRs (Figure S4B).This indicates that genes associated with MPs and MRs of CS, and thereby with the progeroid phenotype, are functionally comparable, and underscores the robustness of our analyses.To note, we observed an enrichment in neurotransmitters transport in MRs, and regulation of transcription in MPs.
We then focused on the functional analysis related to MRs, that we assume to be more prone to influence the phenotype than MPs.
We identified 80 GO terms significantly associated with methylation changes between progeroid and non-progeroid samples (Figure 2).
These GO terms include 59 Biological processes, 8 Molecular functions and 13 Cellular components that in total contained 6398 genes (Table S7, sheet 1), 544 of which were associated with DMRs also identified by our analysis (Figure S4C, listed in Table S7, sheet 2).
We found that DMR-associated genes are mainly enriched in these categories: activation of transcription, embryogenesis/development, nervous system/synapse and membrane transport (colour bars in the bottom of Figure 2, the individual genes are indicated in Table S7, sheet 3).Most of these GOs (64/80) have been associated Filtering and functional GSEA of methylated regions in progeroid versus non-progeroid groups.Plots reporting the correspondence between the enriched GO terms identified in the progeroid vs non-progeroid GSEA analysis (right y-axis) and the genes associated with DMRs, x-axis.This representation displays genes in the x-axis (listed in Table S7, sheet 3), and their corresponding GO term annotation on the y-axis.GO terms on the y-axis are grouped according to GO categories (Biological process, Molecular function, Cellular component).Genes may belong to more than one GO category as well as more than one GO term.Indeed, the same gene can be annotated with in multiple semantically similar terms, and the three GO categories provide complementary biological information.Numbered red frames indicate identified clusters of genes (related GO terms/GO categories); #1: transcription factors implicated in embryogenesis/ development, #2: ion/neurotransmitter membrane transporters, #3: synaptic neuro-developmental genes.These genes are more precisely identified because they belong to at least two independent GO categories (the colour code of the corresponding GO category/ies is/are reported below the x-axis, and their name is reported on the right).
To characterize the functional relevance of the 544 DMRassociated genes, we simultaneously classified all genes according to their GO annotation, which was subdivided in the main GO categories (Biological process, Molecular function, Cellular component), with GO terms clustered by their semantic proximity (right side of Figure 2).This representation allows a more precise definition of differentially methylated genes in CS by revealing genes that belong to at least two independent GO categories (for instance the biological process 'embryogenesis' and the molecular function 'transcription factor').We thus observed that DMRs-associated genes are principally enriched in transcription factors implicated in embryogenesis/development (red frame #1, Figure 2), ion/neurotransmitter membrane transporters (red frames #2) and synaptic neuro-developmental genes (red frames #3).Genes present in only one of the three categories (Biological process or Molecular function or Cellular component) were reported separately (Figure S5).These last clusters included genes involved in skeletal and neuronal development, as well as cytokine production, sensory perception and hormone regulation (red frames #6, #5, #7, #8 and #4, respectively, Figure S5, the individual genes are indicated in Table S7, sheet 4), suggesting that they have a more tissue-specific function than those reported in Figure 2.
We finally compared the 544 DMR-associated genes of our GSEA analysis for their possible implication in human ageing.We used two resources, the AgingAtlas (Aging Atlas Consortium, 2021) and the GenAge (de Magalhaes et al., 2009) databases, which provide a short list of genes (502 and 307 genes, respectively) selected for their role in ageing.Despite the small size of these lists, we found a certain overlap with our 544 DMR-associated genes (15 and 11 genes in common with AgingAtlas and GenAge, respectively; Table S7, sheet 2), six of which common to both lists.

| Epigenetic hallmarks of accelerated ageing in CS and comparison with other progeroid syndromes and normal ageing
To better characterize DNAm remodelling of progeroid CS, we evaluated possible hallmarks of epigenetic ageing in our dataset.
First, we exploited the DNAm dataset to predict the methylation status of repetitive elements, which are consistently hypomethylated in ageing, using the REMP algorithm (Zheng et al., 2017).
This analysis unveiled a clear hypomethylation of Alu sequences in CS samples, in particular in CS-I, compared to WT and UVSS (Figure 3a).Conversely the DNAm pattern of LINE-1 sequences did not appear to be correlated with the progeroid condition (Figure 3b).
We then evaluated our samples for their epigenetic age, which is considered a reliable marker of biological age (Horvath & Raj, 2018).
Using the 'Skin&Blood clock', which was explicitly designed for in vitro studies with fibroblasts (Horvath & Raj, 2018), we found that CS samples tended to have an epigenetic age substantially higher than their chronological age (Figure 3c  We first calculated the correlation matrix between DNAm changes observed in the four datasets (Materials and Methods).
When we considered the entire set of microarray probes (452,567 probes; Figure S7A), the highest correlation was observed between WS and DS datasets (DNAm measured in whole blood in both the cases).CS dataset showed the highest correlations with MSC dataset (positive correlation) and with HGPS dataset (negative correlation).When we focused on the 11,597 CS-DMPs (Table S2 and  We then confronted the 1817 CS-specific DMRs with DMRs identified in the other datasets.No DMRs were identified in WS using the same significance threshold applied to the other datasets.

As summarized in
Conversely, 25% (455 DMRs), 17% (302 DMRs) and 33% (605 DMRs) of CS-specific DMRs were also significant in HGPS, DS and normal ageing, respectively (Table 2).These overlaps were greater than expected by chance (p-value < 0.05, Fisher's exact test; odds ratio of 2.1, 1.5 and 2 for HGPS, DS and MSC datasets, respectively).Fortyeight DMRs were common to all datasets (except WS; Figure S7F and Table S3, sheet 2).Differently from DMPs, the classification of DMRs as hyper-or hypo-methylated is not a straightforward task, as the same DMR can include both probes that gain and probes that lose methylation compared to controls.We assigned the DMRs as hypermethylated or hypomethylated by considering the methylation changes of the most significant CpG probe within each DMR, that

| Correlation between CS-specific DNAm and gene expression
Transcriptomic data were not available from our samples; we therefore assessed by RT-qPCR the expression of several key genes identified as differentially methylated in the progeroid vs non-progeroid comparison.We selected DMRs according to one or more of the following criteria: the lowest p-values, the lowest deltas of methylation between WT and UVSS, presence also in related datasets (normal ageing and/or other progeroid diseases, Figure 4a and Figure S8).
Fourteen of the selected genes were part of the StringentDMRs list.
The correlation between DNAm and transcription for each cell sample and for each of the 31 selected genes is shown in Figure 4a and (a) For each cell type of our dataset, the expression of selected genes was plotted against the percentage of their DNAm.Out of 31 genes tested, 11 genes with a significant correlation between DNA methylation and gene expression are shown, eight of which were inversely correlated (EPB41, EPB49, ASAH1, HOXA11, VARS, SLC1A5, PRDM16, ZIC1), and three were directly correlated (SLC7A1, NDUFC2, CLDND1).
The linear regression curve is indicated in red.The Pearson's correlation coefficient was used to assess the correlation between methylation and transcription levels.The R squared (R 2 ) and p-value (p) are indicated in each graph.The correlation was considered significant when p < 0.05 (gene name underscored in yellow).The characteristics of selected genes are indicated below.The direction of methylation changes (hypermethylation and hypomethylation in the progeroid vs non-progeroid group) is indicated with gold and green highlight of the x-axis label (% Methylation), respectively.Genes that showed differential methylation also in normal ageing and/or other progeroid diseases (10/11) are framed in black (rectangles).Genes present in the top list of DMRs and/or StringentDMRs (6/11) are framed in dark blue (rounded rectangles).Genes selected specifically for their function (HOXA11, CLDND1) are framed in grey (ovals).dTFs are identified with a red square on the left upper corner.CLDND1 was selected for being specifically differentially methylated in CS.These characteristics/criteria are not mutually exclusive.The remaining 20 genes are shown in Finally, out of the 206 StringentDMRs-associated genes, 88 (43%) were differentially expressed in at least one transcriptome dataset (Figure 4e).Several StringentDMRs-associated genes were similarly up-or down-regulated in two different CS transcriptome datasets (e.g., PDE4B), and one, ZIC4, in three datasets (Figure 4e).
Although no StringentDMRs-associated gene was consistently differentially expressed in the four datasets, one gene associated with the DMR list, IRX3 was identified in the four CS datasets.We then assessed the expression of these last three genes in the fibroblasts used in the present study.RT-qPCR analysis revealed significantly higher levels of ZIC4 and IRX3 transcripts in CS (progeroid conditions) versus WT + UVSS cells (non-progeroid conditions), and no changes for PDE4B (Figure 4f), indicating a possible correlation also in this paradigm.

| DISCUSS ION
The cause of premature ageing features in Cockayne syndrome remains largely unknown, and increasing evidences suggest that the DNA repair defect alone does not explain the progeroid phenotype (Velez-Cruz & Egly, 2013).Our analysis relies on the methodological assumption that DNAm changes that are relevant for the progeroid phenotype of CS are not shared with UVSS (mutated but not progeroid).Using this approach, we identified a methylation signature that specifically characterizes CS and can be associated with its progeroid phenotype.Furthermore, using well established hallmarks of ageing, we confirmed the existence of an epigenetic basis of accelerated ageing in CS.
A constraint of this type of studies is that, due to the very low prevalence of CS and especially the UVSS conditions, this DNAm dataset has a limited sample size.This was also the case of recent studies on DNAm profiles in rare diseases that identified diseasespecific epigenetic signatures, and relied on a small number of samples (3-4 per group; Levy et al., 2022;Vento-Tormo et al., 2016).
Specific methods such as the limma R package have been developed to palliate the known limitation and power of differential expression analysis with few replicates, which we also used.Additionally, pointing to decrease the risk of false positives or biased associations inherent to the small sample size of the dataset, our methodological approach is based on the use of multiple and diverse statistical analyses of DNAm data, and the integration of complementary measurements (DNAm and expression), which collectively provide robustness to CS-specific epigenomic signatures (Kolodziej-Wojnar et al., 2020).

| A CS-specific epigenetic signature associated with the progeroid phenotype
To identify DNAm changes that are linked to the precocious ageing phenotype of CS patients, we compared 'Progeroid' (CS) vs 'Non-Progeroid' (WT + UVSS) samples in a differential analysis.We identified 11,597 DMPs in CS, 85% of which were hypomethylated, consistent with a study performed during regular ageing in cultured cells (Fernandez et al., 2015) but differently from DMPs identified in progeroid WS and in the ageing blood, where only 45% and 46% of DMP were hypomethylated, respectively (Guastafierro et al., 2017;Marttila et al., 2015).Moreover, in CS, hypomethylated DMPs were enriched in isolated CpGs (Open Sea) whereas hypermethylated DMPs were enriched in regulatory regions flanking gene promoters (Shore, Shelf; Heyn et al., 2012;Madrigano et al., 2012).Thus, despite they represent the minority of DNAm changes, hypermethylated DMPs are more prone to impact on gene transcription (Ciccarone et al., 2018) in CS.
To better evaluate the impact of DNAm changes on genes, we focused on DMRs encompassing genic sequences.This analysis identified 1817 DMRs in CS, 78% of which were hypomethylated.
These DMRs cover 1498 genes, since some DMRs are present more than once in the same gene.We finally applied a filtering to retain only the DMRs where CS samples clustered separately from WT + UVSS, and that are reasonably specific to the progeroid/ neurodegenerative phenotype (uncoupled from the DNA repair deficiency).Using this pipeline and thanks to the pivotal role of UVSS samples (mutated but non-progeroid), we identified a CSspecific epigenetic signature highlighting a restricted list of 222 DMRs (64% hypomethylated).Successive steps of selection (from DMPs to DMRs, and then to StringentDMRs) obviously restricted the number of hits but, importantly, the main gene categories remained highly represented along the successive analyses (see below).Indeed, the top 20 hypo-and hypermethylated associated genes in the 222 StringentDMRs that constitute the CS epigenetic signature, in large part encode transcription factors implicated in developmental processes and membrane transporters (ions, amines, lipids) that is, two categories also found in DMRs analyses of single genes as well as pathways (see below).Importantly, the relevance of this concentric strategy to identify ageing-related differentially methylated gene in CS was confirmed by the observation that half of the top hits correspond to genes functionally validated by experimental and genetic approaches in humans or animal models for playing a role in ageing or age-related processes.
Consequently, the other half of genes constitutes a pool of novel genes potentially involved in ageing or specific to accelerated ageing, revealed by our approach, and which will be necessary to validate in further studies.
We do not exclude that our analysis, in particular for StringentD-MRs, has been too stringent and relevant genes have been excluded, and/or these findings could be complemented if novel available UVSS fibroblasts at a compatible passage number become available or additional significant CS cases are studied.However, we consider that the major common DNAm differences result from samples derived from patients with different genetic and clinical conditions.
Analysis of pathways rather than single genes, using a GSEA and the GO annotation, revealed GO terms that show statistically significant concordant differences between progeroid and non-progeroid samples.Genes associated with DMRs were particularly enriched (high NES score) in the following biological processes/functions: (i) developmental transcription factors and regulators, (ii) ion/neurotransmitter transporters and (iii) synaptic neuro-developmental genes, in remarkable agreement with single gene analyses above.
i. Gene enrichment in embryogenesis/developmental transcription factors (dTFs) prevails over all other GO terms, suggesting a major implication in CS-specific defects, and to a larger extent than observed in other precocious ageing diseases (Horvath et al., 2015) and during regular ageing (Marttila et al., 2015).
These findings are consistent with the quasi-programmed theory of ageing that claims detrimental over-activity of developmental factors during ageing (Blagosklonny, 2013).According to this theory, if processes important early in life are not properly switched off in the adult, they progressively lead to ageing-related alterations.Accordingly, developmental genes are normally tightly controlled and fail to be so during ageing (Yang et al., 2018).We report modifications of DNAm in CS that are possibly responsible for deregulation of developmental genes, suggesting that this mechanism, whether or not of stochastic origin, plays a role in the ageing phenotype.
Consistently, one of the main 'rejuvenating' interventions proposed to delay ageing is partial cellular reprogramming (Ocampo et al., 2016), which, by transiently turning on master regulators (OCT4, KLF4, SOX2, cMYC) involves re-activation of developmental programmes (Wang et al., 2020).Interestingly, the amelioration of ageing signs during partial reprogramming was shown to be driven by epigenomic remodelling (Ocampo et al., 2016), which supports the notion of a link between DNAm changes and developmental TFs unveiled in our study.Indeed, partial reprogramming was described to lead to a reduction of epigenetic age in human cells (DNAm-based clock; Olova et al., 2019), and restoration of youthful whole DNA methylation patterns in old mice (Lu et al., 2020).
ii. Functional analysis of CS methylated regions reveals enrichment also in ion/neurotransmitter transporters, both passive and active carriers.Passive transporters were mainly voltage-gated potassium channels (KCNA3, KCNA4, KCNC2, KCND2, KCNE4, KCNG4, KCNH1, KCNH2, KCNIP4, KCNQ1), which are implicated in regulating the membrane potential in excitable cells like neuronal or muscle cells (Zironi et al., 2010), but are also expressed in fibroblasts (Li et al., 2009).The involvement of these transporters in neurodegeneration has been evoked in several neurodegenerative diseases (i.e., Alzheimer's and Parkinson's disease; Sesti, 2016), and they may play a role in neurodegeneration in CS.Interestingly, remodelling of K + channels in dermal fibroblasts has been associated with the successful ageing of centenarians (Zironi et al., 2010), whereas overexpression of K + channels was observed in HGPS cells (Zironi et al., 2018).
Genes for active transmembrane transporters mainly involve the solute carrier group of transport proteins (SLCs) that carry a wide variety of substrates including ions, nutrients and/or neurotransmitters through plasma/organelle membranes, also in fibroblasts (De Sanctis et al., 2016).Only a few members of this large family of transporters (>400 proteins) have been associated with ageing (Chekroud et al., 2012).SLC1A5 and SLC7A1 showed a strong direct or inverse correlation, respectively, with DNAm levels in CS patient and control samples.
iii.Gene enrichment in the functional analysis of CS methylated regions also identified synaptic neuro-developmental genes.
This finding is of relevance since progeroid CS is also characterized by brain and nervous system abnormalities (Karikkineth et al., 2017), providing a novel pool of targets to analyse in the context of age-related neuronal pathophysiology (e.g., neurodegeneration).
To note, additional genes with just one of these characteristics, namely developmental genes, membrane transporters and synapsis-related genes, were also enriched in our analysis.Further, we observed components of the cell surface and extracellular compartment, including different types of cell adhesion, namely adherens junctions (ACTN1, APC, CDH3,4,7,15,JUP,DSP), tight junctions (CLDN6, 10) and gap junctions (GJA4, GJA8, GJC1).The expression of CLDN1 strongly correlated with DNAm levels in our samples.
Cell adhesion factors play an important role in maintaining tissue architecture as well as tissue homeostasis due to their involvement in the formation of epithelial barriers that regulate exchanges, and is therefore not surprising that they are possibly involved in ageing (Parrish, 2017).Accordingly, higher expression of protein involved in intercellular (tight) junctions was recently found in healthy centenarians compared to controls (Santos-Lozano et al., 2020).The role of these junctions, particularly during ageing, remains poorly understood outside epithelia to which they seem restricted, despite fibroblasts express adhesion molecules (Ko et al., 2000).

| Epigenetic hallmarks of accelerated ageing in CS
We applied to our dataset a methylation age predictor specifically developed for human fibroblasts (Horvath & Raj, 2018), which outperforms the pan-tissue age estimator (Horvath, 2013) and other clocks for this cell type.Using the Skin&Blood clock, we showed acceleration of epigenetic age specifically in CS, but not in UVSS.Although we cannot exclude that this observation is related to the small sample size and/or uneven distribution of sexes across the groups, we think that these possibilities are unlikely.Indeed, epigenetic age calculated using the Skin&Blood clock was highly correlated with chronological age in an independent dataset, and these additional WT samples showed age acceleration values similar to WT and UVSS samples in our study.Moreover, whereas the pan-tissue epigenetic clock has been shown to return higher epigenetic ages in males compared to females, this has not been reported for the Skin&Blood clock; accordingly, epigenetic age acceleration was significantly larger in the progeroid compared to the non-progeroid group also when correcting for sex as possible confounding factor.Notably, the mean epigenetic age acceleration in CS was at least 15.5 years, that is, higher than those previously reported in blood from WS or DS (addition of 6.4 and 6.6 years, respectively, for 10-60 year old patients) (Horvath et al., 2015;Maierhofer et al., 2017) (Borelli et al., 2015;Horvath et al., 2015).
In addition, CS fibroblasts displayed other epigenetic changes that have been largely reported in studies on ageing (Ciccarone et al., 2018).CS cells tended to have lower DNAm values compared to UVSS or WT cells, resembling the global hypomethylation described in physiological ageing contexts (Madrigano et al., 2012), and ascribed to heterochromatin loss in non-genic repetitive sequences with ageing (Ciccarone et al., 2018).By inferring the methylation of locus-specific repetitive elements from microarray data, we found that CS cells displayed a marked hypomethylation of Alu elements, while LINE elements were not affected by the syndrome.Thus, the methylation patterns of CS cells exacerbate those that are characteristic of normal ageing (Cardelli, 2018).This result seems specific for our dataset, as in whole blood from WS patients no differences in DNAm of repetitive elements (experimentally measured) were found compared to healthy controls (Maierhofer et al., 2019), or this is due to differences between fibroblasts and blood cells.

| CS has DNAm characteristics of regular ageing and common marks with progeroid conditions
The CS-specific epigenetic signature was compared with physiological ageing (Fernandez et al., 2015) and other progeroid diseases [DS (Horvath et al., 2015), WS (Guastafierro et al., 2017) for DS and WS), a factor that may affect our results.
Remarkably, several CS characteristics were distinct from normal ageing.For instance, significant hypomethylated regions were more numerous in CS than normal ageing, which rather displayed hypermethylated regions (Marttila et al., 2015).Moreover, in normal ageing hypermethylation was associated with developmental processes and DNA binding/transcription, whereas hypomethylation did not enrich a specific set of genes.Genes involved in transcriptional regulation of developmental processes, as indicated above, were prevalently hypomethylated in CS.Similarly, the enriched functional GO categories identified by our GSEA analysis are prevalently composed (>95%) of genes unknown in databases of human age-related genes, but these databases are far from being exhaustive.
From the CS-specific pool, we identified 48 DMR-related genes that are common to the other datasets analysed (normal ageing, HGPS, DS).Among them, the expression of EPB49, KANK2, and to some extent also HOXA9, was inversely correlated with DNAm, suggesting that these methylation changes are relevant for CS.EPB49 has been recently associated, together with other components of the haem metabolism, with human ageing (Timmers et al., 2020).
Other genes in this list, like NAGS, coding for a N-acetyl glutamate synthase linked to glutamine/glutamate metabolism, and GNL1, coding for a little-known GTPase reported to modulate cell proliferation (Krishnan et al., 2018), have not been associated with ageing, and the function of the protein coded by FAM43A is still unknown.

| Top-ranked genes with DNAm changes display a good correlation with expression changes
A general trend in ageing-related DNAm studies is that changes (hypo-or hypermethylation) are not necessarily correlated with changes in gene expression, even in regions with a highly reproducible epigenetic remodelling during ageing that define the 'methylation clock' (Horvath, 2013).We first analysed transcription levels of key genes directly in patient-derived fibroblasts.Transcription of 31 genes that displayed the highest DNAm changes in CS (top 50 hypo/ hypermethylated genes) and/or were in common with other ageing datasets, unveiled a tendency to inverse correlation (high expression vs low methylation levels), although direct correlation was observed too, and one third of cases was statistically robust.The remaining cases showed no correlation or a tendency to direct correlation, at least within the limited number of patient samples available.The DNAm/transcription correlation appears remarkably larger than in reported cases of ageing-related datasets (for instance less than 4% in normal ageing; Marttila et al., 2015), perhaps as a result of the robustness of our data analysis (DMR instead of DMP, multiple filtering processes based on biological assumptions), or a characteristic of DNAm changes in CS.Moreover, we may have underestimated the correlation between DNAm and transcription, since our experimental procedures (i.e., the use of Infinium microarray on bisulfite-treated DNA), like the published databases that we used for comparisons, do not distinguish DNA methylation from DNA hydroxymethylation.The latter epigenetic modification reported a more direct correlation with transcription in ageing fibroblasts.
To further understand the consequence of DNAm changes on gene expression, as genome-wide transcriptomic data are not available in primary fibroblasts used in this study, DNAm/transcription correlations were assessed using four published CS transcriptomes performed in post-mortem cerebella (Okur et al., 2020;Wang et al., 2014), iPSC-derived neurons (Vessoni et al., 2016) and MSCs (Wang et al., 2020).These datasets were relevant because they were originated from patient-derived cells or native organs, We are aware that other DNAm changes may not be associated with transcription alterations, and further genome-wide studies will determine case by case when this occurs.Interestingly, this analysis identified expression changes mainly in genes enriched in transmembrane ion channels of the nervous system, in agreement with neurodegeneration observed in CS patients and also during regular ageing (Sesti, 2016).Among the genes identified in this analysis, the differentially methylated genes ZIC4 and IRX3 were of particular interest as they were consistently down-regulated or up-regulated, respectively, in multiple CS transcriptomic datasets.In patient-derived fibroblasts tested here, both ZIC4 and IRX3 were up-regulated.Consistent upregulation of IRX3 in most datasets from different patients, and distinct cells/native organs makes it an interesting candidate for future mechanistic studies.

| CON CLUS IONS
In summary, we identified CS-specific DNAm changes by filtering out changes not associated with the progeroid and neurodegenerative phenotype, a unique characteristic of this pathological paradigm.
Functional analyses revealed enrichment in developmental transcription factors, ion/neurotransmitter transporters and synaptic neuro-developmental genes, that appear thus critical for accelerated ageing/neurodegeneration.A large fraction of DNAm changes have been associated with expression changes in available datasets, in particular post-mortem cerebella.In addition to epigenetic signatures in common between CS and physiological ageing, DNAm changes observed in CS may help identifying transcriptional changes that would appear at a later time in regular ageing.Finally, our results demonstrate that CS is a pertinent model to study ageing and unveil pathways and genes that could be relevant for future mechanistic studies.
Figure S2I-L).The global pattern of methylation values was not affected by the ComBat batch correction (Figure S2M,N).We performed the subsequent analyses using the batchcorrected data, removing the probes on X and Y chromosomes (441,982 probes left) and, for the seven replicated samples, averaging methylation values from the 450 k and EPIC platforms.PCA and boxplots calculated on this merged dataset are reported in Figure 1b,c, respectively.Interestingly, boxplots of the methylation values showed a global hypomethylation in CS cells compared to WT and UVSS (Figure 1c).
. (b) Principal component analysis (PCA) of DNAm levels of probes assessed in the Infinium450k and InfiniumEPIC beadchips (merged dataset), after batch correction; dots are coloured according to the pathological (CS, UVSS) or control (WT) condition.(c) Boxplot of global percentage of DNAm in the merged dataset after batch correction.(d) Graphical representation of the approach used to identify CS-specific (progeroidrelated) epigenetic signature.'CS-specific' refers to epigenetic changes present in CS but not in UVSS samples that, although harbouring the mutation, do not show the progeroid phenotype.Progeroid-associated epigenetic changes result from comparing DNAm between CS ('Progeroid' group) versus WT and UVSS cells, merged in the same 'Non-Progeroid' group.A CS-specific CpG site will tend to have similar DNAm values in CS samples and differ from those in WT and UVSS cells (top panel).On the contrary, a CpG site whose DNAm is associated with the phenotype shared between CS and UVSS, like photosensitivity, will have similar values in CS and UVSS cells and differ from WT cells (bottom panel).Analysis of variance between the progeroid and the non-progeroid group will tend to return significant p-values in the first but not in the second case.(e) Volcano plot showing nominal p-values vs the difference in mean percentage of methylation between the progeroid and the non-progeroid group.Green: hypomethylation; gold: hypermethylation.(f) Enrichment of DMPs in different genomic regions.Odds ratio values are reported on the y-axis.Green: hypomethylation; gold: hypermethylation.The position of the various regions is schematized in the bottom.(g) Heatmap and hierarchical clustering on DMPs methylation values in the various patient and donor-derived cells.TA B L E 1 Top 20 hypo-and hypermethylated StringentDMRs of progeroid versus non-progeroid conditions.
and FigureS6A).Epigenetic age acceleration, calculated as the residuals from the linear regression of epigenetic age on chronological age, was significantly higher in progeroid compared to the non-progeroid (WT and UVSS) group (pvalue < 0.05; difference in mean residuals between the two groups: 15.5 years; Figure3dand FigureS6B) and was particularly evident in CS-I samples.Statistical significance was observed also upon correction for sex.Other available epigenetic clocks, like the pan-tissue clock(Horvath, 2013; Figure S6E,F), the Hannum's clock(Hannum et al., 2013; Figure S6I,J), the PhenoAge(Levine et al., 2018;  Figure S6M,N) and their derivates based on principal component analysis (Higgins-Chen et al., 2022; Figure S6C,D,H,K,L,O,P), were notable to detect a significant epigenetic age acceleration in CS.This was in part expected, as these clocks were reported to perform poorly in fibroblasts(Horvath, 2013) or were developed for whole blood samples(Hannum's clock and PhenoAge).The analysis of an independent dataset including 31 neonatal and 7 adult untreated fibroblasts (GSE197724) confirmed that only the Skin&Blood was accurate in predicting the age of healthy control (WT) samples, supporting its higher reliability in these experimental settings (Turquoise symbols; all panels of FigureS6).To assess whether CS-specific DNAm changes that emerged from the progeroid vs non-progeroid comparison are shared with normal ageing and/or other progeroid diseases, we assessed available Infinium datasets including: (1) other progeroid syndromes: HGPS (fibroblasts from the GSE196148 dataset), WS (whole blood;Guastafierro et al., 2017) and DS (whole blood;Horvath et al., 2015); and (2) mesenchymal stem cells (MSC) from individuals of different age, as a proxy of age-associated changes in human fibroblasts(Fernandez et al., 2015).In HGPS, WS and DS datasets, we compared progeroid samples versus healthy controls; in the MSC dataset, we compared cells from old donors versus cells from young donors.For each dataset, we retrieved the lists of DMPs and DMRs, using the same analytical pipeline applied for CS dataset (FigureS1).

Figure S7B )
FigureS7B), the positive correlation of the CS-MSC pair increased, suggesting that DNAm changes observed in CS partially resemble those related to normal ageing.On the other side, also the negative correlation of the CS-HGPS increased, suggesting that DNAm is differentially remodelled in the two progeroid conditions.These results were supported by an alternative analysis of the data, in which we evaluated the overlap between the DMPs and DMRs lists in the four comparisons (Table2 and FigureS7C,D).
is, the probe with the smallest BH-corrected p-value resulting from the ANOVA test.To note, also taking into account the complexity of defining the direction of methylation changes in DMRs, none of the 48 DMRs changed in the same direction in all categories (CS, HGPS, DS and normal ageing; Figure S7G,H).Similarly to what observed for DMPs, MSC ageing had the highest proportion of DMRs concordant with CS.

Figure S8 .
Figure S8.Fifteen genes showed a strong (0.4 ≤ R2 ≤ 0.9) and seven genes a moderate (0.4 ≤ R2 ≤ 0.2) correlation between DNAm and RNA expression.The majority of tested genes displayed an inverse methylation/transcription correlation whereas 12 DMRs were poorly FigureS7.(b) Heatmap of odds ratios representing the strength of association between genes associated with DMPs, DMRs or StringentDMRs and genes differentially expressed in different CS transcriptomic datasets (from iPSC-derived Neurons(Vessoni et al., 2016), iPSC-derived MSCs(Wang et al., 2020) and post-mortem Cerebellum [_1(Wang et al., 2014) and _2(Okur et al., 2020)].The statistical significance of overlaps was tested with Fisher's exact tests, and the corresponding pvalues are indicated in orange.(c) Venn diagram showing the total number of common genes between genes associated with StringentDMRs and genes differentially expressed inWang et al. (2014) (Post-mortem Cerebellum_1) (upper panel).Venn diagrams displaying among the Hypo-(blue) or Hypermethylated (magenta) StringentDMRs, number of genes in common with genes found Down-(lower-left panel) or Upregulated (lower-right panel) in the same dataset.(d) Clustering of semantically related GO terms (left panel; molecular function right panel; biological process) identified in GSEA (green circles) and in the enrichment analysis of Wang et al. (2014) (purple circles) using REVIGO.Common GO terms are represented by the orange circles.Semantically similar terms are grouped in clusters identified by a 'meta-term' (black).The significativity of GO term overlap and odds ratios are displayed in red.(e) Heatmap representing the direction of expression changes (green; up-regulated, red; down-regulated) of the 88 unique genes in common between genes associated with StringentDMRs and genes differentially expressed in the four previously mentioned CS transcriptomic datasets.Genes already tested for expression in Figure4aand FigureS7are underlined in orange.(f) Quantitative RT-qPCR of ZIC4, PDE4B and IRX3 transcripts in the progeroid vs non-progeroid groups.ZIC4 is common to StringentDMRs and three out of four transcriptomic datasets (TDs); PDE4B is common to StringentDMRs and 2 out 4 CS TDs, and IRX3 is common to DMRs and the 4 CS TDs.n = 3 independent experiments; mean ± SD; n.s, non-significant; p, p-value; unpaired t-test vs non-progeroid.a direct correlation (14 genes: hypomethylated/down-regulated; 7 genes: hypermethylated/up-regulated; Figure4c).Importantly, comparison of enrichment analyses demonstrated a strong association between the enriched GO terms identified [Odds ratio: 130.6 (Molecular function) and 69.8 (Biological process); Figure4d] in the present DNA methylation study and in the transcriptome of CS cerebella fromWang et al. (2014).Moreover, in this last comparison, the semantic clustering of common enriched GO term with REVIGO revealed that genes displaying a methylation/transcription correlation are enriched in transmembrane ion transporters of the nervous system (Figure4d).
allowing a comparison with patients-derived fibroblasts in our study, and in the absence of immortalization processes, which may interfere with pathways implicated in ageing.All these datasets displayed a strong association (Odd ratios > 1) with DNAm changes observed in our study (at the level of DMPs and DMRs, or StringentDMRs).It should be underscored that expression changes observed in these datasets may likely include changes occurring also in UVSS cells, and that are consequently linked to DNA repair deficiency rather than the precocious ageing defects.To reduce this possible bias, we focused on the CS transcriptome of cerebella fromWang et al. (2014) that displayed the strongest and most significant association with expression changes observed in StringentDMRs, where DNAm changes in common with UVSS have been removed.Results confirmed a robust DNAm/expression correlation, obviously at a larger scale than the one assessed in patient fibroblasts where only 31 top-ranked genes were tested.To note, this large-scale correlation was strong despite DNAm data originated from CS patient-derived fibroblasts in culture, and transcriptome data from post-mortem cerebella of other CS patients (i.e., carrying distinct clinical defects and diverse genomic backgrounds), underscoring the robustness of our findings.

CHR CpG island name Relation to the CpG island Gene(s) Adjusted p-value Description
Note: DMRs are ranked according to the adjusted p-value.Mapping information on DMRs are given by the chromosome number, the CpG Island coordinate, the relation respect to the CpG island and the associated gene.A description of each gene is also provided.Bold lines indicate StringentDMRs and corresponding genes assessed for expression.Abbreviations: /, unknown; dTF, developmental transcription factor; Ms, muscle; Mt, mitochondrial; Rb, ribosome; TF, transcription factor.

Table 2
, out of 11,597 CS-specific DMPs, 1822 (16%) were also detected in the HGPS dataset, 1 in the WS dataset, datasets showed the same direction of epigenetic changes (hyperor hypomethylation compared to the respective controls) while this percentage was lower for HGPS and DS datasets (51.21 and 38.83% of concordant DMPs, respectively; Table 2 and Figure S7D,E).
Epigenetic age acceleration in CS and commonalities with regular ageing.Density of DNAm of each sample for the repetitive Alu (a) and LINE-1 (b) elements.(c) Scatter plots of epigenetic age, calculated with the Skin&Blood clock, vs chronological age.The black line represents the regression of epigenetic age on chronological age.The age of samples (only CS) with a filled symbol was calculated on the last report/death since the age at the time of biopsy was not known, thereby their chronological age (xaxis) is likely overestimated.(d) Boxplots of epigenetic age acceleration values according to the Skin&Blood clock in WT, UVSS and CS groups.
E 3 TA B L E 2 Comparison of the CS-specific DMPs/DMRs with those emerged in other progeroid conditions and normal ageing.Comparison In common with CS Odd ratio Hypomethylated/Hypermethylated Concordant hypo−/ hypermethylation (%) Note: Number of common DMPs/DMRs between those identified in progeroid vs non-progeroid (11,597 DMPs, 1817 DMRs) conditions and DMPs/ DMRs identified in other studies performed in progeroid syndromes (HGPS, WS, and DS) or during normal ageing (MSC).For each, the odd ratios of Fisher's Exact Tests (p-value < 0.05), the number (out of total) and percentage of concordant hyper-/hypo-DMPs/DMRs (same direction) is also provided.