Epigenetic profiling of Italian patients identified methylation sites associated with hereditary transthyretin amyloidosis
Clinical Epigenetics volume 12, Article number: 176 (2020)
Hereditary transthyretin (TTR) amyloidosis (hATTR) is a rare life-threatening disorder caused by amyloidogenic coding mutations located in TTR gene. To understand the high phenotypic variability observed among carriers of TTR disease-causing mutations, we conducted an epigenome-wide association study (EWAS) assessing more than 700,000 methylation sites and testing epigenetic difference of TTR coding mutation carriers vs. non-carriers. We observed a significant methylation change at cg09097335 site located in Beta-secretase 2 (BACE2) gene (standardized regression coefficient = −0.60, p = 6.26 × 10–8). This gene is involved in a protein interaction network enriched for biological processes and molecular pathways related to amyloid-beta metabolism (Gene Ontology: 0050435, q = 0.007), amyloid fiber formation (Reactome HSA-977225, q = 0.008), and Alzheimer’s disease (KEGG hsa05010, q = 2.2 × 10–4). Additionally, TTR and BACE2 share APP (amyloid-beta precursor protein) as a validated protein interactor. Within TTR gene region, we observed that Val30Met disrupts a methylation site, cg13139646, causing a drastic hypomethylation in carriers of this amyloidogenic mutation (standardized regression coefficient = −2.18, p = 3.34 × 10–11). Cg13139646 showed co-methylation with cg19203115 (Pearson’s r2 = 0.32), which showed significant epigenetic differences between symptomatic and asymptomatic carriers of amyloidogenic mutations (standardized regression coefficient = −0.56, p = 8.6 × 10–4). In conclusion, we provide novel insights related to the molecular mechanisms involved in the complex heterogeneity of hATTR, highlighting the role of epigenetic regulation in this rare disorder.
Hereditary transthyretin amyloidosis (hATTR; OMIM#105210) is a life-threatening disorder caused by transthyretin (TTR) protein misfolding. This causes amyloid fibril deposition in several tissues (e.g., peripheral nerves, heart, and gastrointestinal tract) [1, 2]. hATTR is characterized by extreme clinical heterogeneity including age of onset, penetrance, and clinical display [3,4,5]. To date, more than 130 amyloidogenic mutations have been identified in the coding regions of the TTR gene, which are the cause of hATTR . The prevalence of hATTR is estimated to be approximately 1/100,000 . hATTR endemic areas have been identified in Portugal and Sweden [4, 5]. Although both of these regions are affected by the same amyloidogenic mutation, Val30Met (rs28933979), the penetrance and age of onset are different: early age of onset and high penetrance in Portugal [4, 5, 8, 9] vs. late age of onset and low penetrance in Sweden and in non-endemic countries [3, 10, 11]. hATTR phenotypic heterogeneity is likely due to the contribution of genetic and non-genetic factors involved in the complex genotype–phenotype correlation observed [12,13,14,15,16,17,18]. Recent data strongly support the role of non-coding regulatory variation on TTR gene expression, as one of the mechanisms affecting the phenotypic manifestations observed in carriers of TTR amyloidogenic mutations [19,20,21,22]. Among genomic regulatory features, epigenetic modifications are key mechanisms in modulating a wide range of molecular functions and potential targets to develop novel treatments [23,24,25]. Of several epigenetic modifications, DNA methylation is the most studied with respect to human traits and diseases . In the context of monogenic disorders, methylation studies investigate the role of epigenetic changes involved in the phenotypic expression observed among carriers of disease-causing mutations [26,27,28]. While epigenetic modifications have the potential to be involved in hATTR pathogenic mechanisms, to our knowledge no study has explored methylation changes of patients affected by this life-threatening disease. In the present study, we conducted an epigenome-wide association study (EWAS) to identify DNA methylation sites associated with hATTR, investigating 48 carriers of TTR amyloidogenic mutations and 32 controls, comparing hATTR affected patients, asymptomatic carriers, and non-carriers.
We compared the methylation changes (measured as M values, i.e., the log2 ratio of the intensities of methylated probe versus unmethylated probe; beta-values are also graphically presented to provide more biologically interpretable data) in the peripheral blood of 48 carriers of TTR amyloidogenic mutations and 32 controls (non-carriers). Testing more than 700,000 methylation sites, an association survived epigenome-wide false discovery rate correction (FDR q < 0.05) at the cg09097335 site located in Beta-secretase 2 (BACE2) gene body (standardized regression coefficient = -0.60, p = 6.26 × 10–8, FDR q = 0.044). Carriers of TTR amyloidogenic mutations showed a significant hypomethylation when compared to controls (beta-value plot: Fig. 1; M value plot: Additional file 1). To understand whether methylation at this CpG site is associated with disease-associated genetic differences or post-disease processes, we compared hATTR patients, asymptomatic carriers of TTR mutations, and controls. Significant differences were observed for (i) hATTR patients vs. controls (standardized regression coefficient = −0.402, p = 5.7 × 10–4; Additional file 2, beta-value and M value plots) and ii) asymptomatic carriers vs. controls (standardized regression coefficient = −0.716, p = 3.21 × 10–5; Additional file 2: beta-value and M value plots), but no difference was present between hATTR patients vs. asymptomatic carriers (standardized regression coefficient = 0.137, p = 0.332). Similarly, considering different TTR mutations, we observed significant differences in cg09097335 methylation (i) between Val30Met carriers vs. controls (standardized regression coefficient = −0.587, p = 1.85 × 10–6; Additional file 2: beta-value and M value plots) and (ii) carriers of other TTR mutations vs. controls (standardized regression coefficient = −0.479, p = 1.27 × 10–3; Additional file 2: beta-value and M value plots), but not between Val30Met carriers vs. other TTR mutation carriers (standardized regression coefficient = 0.093, p = 0.424). Leveraging GTEx data , we observed a complementary transcriptomic regulation between TTR and BACE2 genes where the first is mainly expressed in its source organ (i.e., liver), while the second is expressed in target organs of TTR amyloid deposits (Fig. 2, Upper Panel). We investigated interactive proteins related to TTR and BACE2 loci based on multiple experimental and computational evidence, identifying five candidates with medium-to-highest interaction confidence (Fig. 3). These include FYN (FYN proto-oncogene, Src family tyrosine kinase; interaction score = 0.809), BACE1 (Beta-secretase 1; interaction score = 0.804), APP (amyloid-beta precursor protein; interaction score = 0.430), IGHV3-11 (immunoglobulin heavy variable 3–11; interaction score = 0.412), and ENSG00000259680 (uncharacterized protein similar to an immunoglobulin heavy variable 3/OR16 gene; interaction score = 0.412). Among them, TTR showed the highest interaction with APP protein (interaction score = 0.936). BACE2 protein interactive network (Fig. 3) showed functional enrichments for several biological processes and molecular pathways (Table 1). Among BACE2-related enrichments surviving FDR multiple testing correction, we observed: Alzheimer’s disease (AD, KEGG hsa05010, FDR q = 2.2 × 10–4) related to the interaction of BACE2 with APP and BACE1; membrane protein ectodomain proteolysis (GO: 0006509, FDR q = 0.007) and amyloid-beta metabolic process (GO: 0050435, FDR q = 0.007) related to BACE2-BACE1 interaction; protein metabolic process (GO: 0019538, FDR q = 0.043) related to the interaction of BACE with APP, BACE1, FYN, and IGHV3-11. The interactions of other proteins within BACE2 interactive network also highlighted amyloid-related functional enrichments: amyloid fiber formation (Reactome HSA-977225, FDR q = 0.008) and response to amyloid-beta (GO: 1904645, FDR q = 0.009) related to the interaction of APP with BACE1 and FYN, respectively.
Within TTR gene region, we observed that Val30Met mutation disrupts a methylation site, cg13139646, causing a drastic hypomethylation in Val30Met carriers when compared with carriers of other TTR mutations (standardized regression coefficient = −2.18, p = 3.34 × 10–11; Additional file 3: beta-value and M value plots). Since Val30Met mutation disrupts the probe extension for the cg13139646 methylation site (the Met allele removes the G nucleotide from the CpG site), its impact on the methylation changes in TTR gene should be confirmed using alternative typing methods. However, to conduct an initial exploratory analysis of the potential functional implications of cg13139646 disruption, we performed a co-methylation analysis with respect to cg13139646 in the non-carrier sample. We identified 34 methylation sites that are correlated with cg13139646 site (Pearson’s r2 > 0.20; Additional file 4). Considering these co-methylated CpG sites, we investigated epigenetic differences among hATTR patients, asymptomatic carriers of TTR mutations, and controls (Additional file 5). Applying a Bonferroni correction accounting for the number of CpG sites tested, we observed a significant methylation difference between hATTR patients and asymptomatic carriers of TTR amyloidogenic mutations at cg19203115 (standardized regression coefficient = −0.555, p = 8.6 × 10–4; Additional file 6: beta-value and M value plots). Nominally significant methylation differences were observed (i) between carriers vs. non-carriers at cg11481443 (standardized regression coefficient = −0.306, p = 3.4 × 10–3; Additional file 6: beta-value and M value plots) and cg02936398 (standardized regression coefficient = 0.177, p = 4.9 × 10–2; Additional file 6: beta-value and M value plots); (ii) hATTR patients and asymptomatic carriers at cg14311811 (standardized regression coefficient = −0.273, p = 3.8 × 10–2; Additional file 6: beta-value and M value plots). Considering symptoms reported by hATTR patients, we identified CpG sites co-methylated with cg13139646 (i.e., the site disrupted by Val30Met mutation) nominally associated with cardiac involvement (cg27392998, standardized regression coefficient = −0.235, p = 7.5 × 10–3; cg18038361, standardized regression coefficient = 0.227, p = 5 × 10–2; Additional file 6: beta-value and M value plots), carpal tunnel syndrome (cg16492377, standardized regression coefficient = −0.229, p = 1.8 × 10–2; Additional file 6: beta-value and M value plots), and peripheral nervous system involvement (cg14719951, standardized regression coefficient = 0.249, p = 3.5 × 10–2; Additional file 6: beta-value and M value plots). Since some of these CpG sites were mapped to loci located near TTR gene (Additional file 5), we analyzed TTR transcriptomic profile in hATTR target organs, observing a different pattern when compared to the expression of the surrounding genes (Fig. 2, Bottom Panel).
hATTR is a rare multi-organ disorder caused by TTR misfolding and consequently amyloid deposition in several tissues . This life-threatening condition is characterized by high clinical heterogeneity with respect to age of onset, penetrance, and phenotypic manifestation [1,2,3,4,5,6,7,8,9,10, 30]. Although TTR amyloidogenic mutations are the cause of TTR misfolding, non-coding variation and modifier genes are hypothesized to be involved in the wide variability of phenotypic manifestations observed in carriers of TTR disease-causing mutations [12, 15, 17,18,19,20,21,22]. Epigenetic modifications (e.g., DNA methylation changes) could also play an important role in the molecular network regulating the hATTR amyloidogenic process . To explore this hypothesis, we conducted an EWAS investigating more than 700,000 methylation sites in 48 carriers of TTR amyloidogenic mutations and 32 non-carriers. A CpG site (cg09097335) located in BACE2 gene was significantly hypomethylated in carriers when compared to non-carriers. This gene encodes Beta-secretase 2, a protein mainly known for its role in cleaving APP protein in amyloid-beta, which is a key factor involved in AD pathogenesis [31,32,33]. Differently from BACE1, which is the primary β-secretase protein cleaving APP to amyloid-beta, BACE2 is poorly expressed in the brain and its cleaving ability increases following an inflammatory response . APP processing occurs via three proteolytic cleavages caused by α- β- and γ-secretase . In non-amyloidogenic processes, α- and γ-secretases lead to the production of a smaller P3 fragment and APP intracellular domain, while, in the amyloidogenic pathway, β-secretase and γ-secretase produce amyloid-beta [35,36,37,38,39]. Our results also showed a high-confidence interaction between APP and TTR. Numerous studies explored the interactions between these two amyloidogenic proteins, displaying a relevant biological role of TTR in amyloid-beta aggregation and clearance in AD patients [40,41,42,43,44]. Specifically, TTR instability reduces the clearance of amyloid-beta, increasing amyloid toxicity in the brain [40,41,42]. Metal ions and interaction with other proteins could also affect TTR stability . Although TTR genetic reduction did not alter APP processing, immunohistochemical and biochemical studies showed that genetic reduction of TTR elevates Aβ deposition in the brains of transgenic mice harboring APPswe/PS1ΔE9TTR + / − transgenes . Under physiological conditions, the APP intracellular domain appears to be involved in epigenetically up-regulation of TTR to increase its amyloid-beta clearance activity . In AD patients with TTR Val30Met, a significant association between amyloid-beta levels and AD was identified [40, 44]. A putative amyloidogenic role of amyloid-beta in hATTR was also identified in a post-mortem analysis of a Val30Met carrier where both TTR and amyloid-beta were deposited in the cerebral leptomeningeal and cortical blood vessel walls with a part of the vessel wall occupied by a combination of TTR and amyloid-beta aberrant proteins . These previous findings strongly indicate an interplay between the pathogenic mechanisms involved in hATTR and AD. Our epigenome-wide study identified BACE2 as a potential key factor in this interaction. As previously discussed, BACE2 protein plays a minor role in APP cleaving in the brain [33, 34], while its activity increases in peripheral tissues under inflammatory response . Our transcriptomic analysis showed that BACE2 is expressed in tissues affected by TTR amyloid deposits (i.e., heart, nerves, colon, small intestine, and adipose tissues). Accordingly, the methylation change observed in the TTR-mutation carriers is possibly due to the role of BACE2 in response to the inflammation induced by TTR amyloidogenic process in peripheral tissues .
Within TTR gene region, we observed that Val30Met disrupts a CpG site, causing a drastic hypomethylation in the carriers of this mutation. SNPs at CpG sites disrupting the methylation reactions can be associated with changes in regulatory function.[48, 49]. Although the impact of Val30Met on cg13139646 has to be confirmed by alternative methods because of the cg13139646 prove extension disruption, we conducted an initial exploratory analysis, testing the co-methylation of this site with CpG sites in the surrounding regions (NC_000018.9: 28,171,000–30,171,500). Indeed, co-methylation patterns reflect specific molecular mechanisms responsible for the regulation of multiple genes located in the same region . In our analysis, some of the CpG sites identified map to TTR gene region, while others map in nearby loci. Considering hATTR target organs, these surrounding loci with co-methylated CpG sites have higher gene expression than TTR. We speculate that these genes may be involved in the formation of TTR amyloid deposits. This hypothesis is supported by the fact that co-methylated CpG sites are associated with hATTR traits. The strongest evidence was observed with respect to cg19203115 mapped in B4GALT6 gene. Considering a Bonferroni correction accounting for the number of co-methylated CpG sites tested, cg19203115 showed a significant difference in methylation levels between hATTR patients and asymptomatic carriers of TTR mutations. B4GALT6 gene encodes beta-1,4-galactosyltransferase 6, a type II membrane-bound glycoprotein that has exclusive specificity for the donor substrate UDP-galactose. B4GALT6 enzyme activity changes in response to inflammatory processes . B4GALT6 stimulates astrocyte activation through the catalyzation of lactosylceramide synthesis, which in turn controls the production of pro-inflammatory cytokines and chemokines . Hence, observing methylation changes in B4GALT6 may be associated with the inflammatory response to TTR amyloid deposits. Nominally significant differences were observed for CpG sites mapped in other surrounding loci: DSC2 (cg02936398, Carriers vs. Controls); DSG2 (cg14311811, hATTR patients vs. asymptomatic carriers); DSC3 (cg16492377, carpal tunnel syndrome in hATTR patients). DSC2 and DSG2 encode components of the desmosome. This protein complex is specialized for cell-to-cell adhesion in myocardial tissue and mutations in DSC2 and DSG2 genes are associated with arrhythmogenic right ventricular cardiomyopathy . In an in vivo study of myocardial inflammation, DSC2 overexpression was observed to lead to tissue necrosis, fibrosis, and calcification of ventricles . This process alters homeostasis among desmosomal proteins, inducing a cascade of different cell–cell interactions leading to cardiac remodeling . In hATTR, cardiac amyloid fibril depositions also led to tissue dysfunctions . Heart failure, restrictive cardiomyopathy, and rhythm disturbances (i.e., conduction system diseases, atrial fibrillation, and ventricular tachycardia) are the main clinical signs that occur after the accumulation of misfolded TTR protein [54,55,56]. Furthermore, transcriptomic interaction is observed between TTR and DSG2 to induce hypertrophic cardiomyopathy in animal models . In this context, methylation changes in DSC2 and DSG2 genes could reflect pathogenic processes in hATTR target organs. We also identified two CpG sites co-methylated with cg13139646 (i.e., the methylation site disrupted by Val30Met mutation) that are nominally associated with hATTR symptoms. Cg18038361 is located in TTR gene promoter region and is associated with cardiac involvement in hATTR patients. DNA methylation changes in promoter regions are well known to play an important role in gene expression regulation [58, 59]. Cg18038361 association may be linked to regulatory changes in TTR gene expression. Lastly, two CpG sites—cg16492377 and cg14719951, map to DSC3 transcription start site and gene body, respectively. In our analysis, methylation changes in these sites were associated with carpal tunnel syndrome and peripheral nervous system involvement in hATTR patients, respectively. DSC3 gene encodes the desmocollin-3 a calcium-dependent glycoprotein. Low DSC3 expression in human epidermis leads to a loss of tissue integrity . We speculate that cg16492377 methylation association with carpal tunnel syndrome may be related to changes in DSC3 transcriptomic regulation.
Although we provide novel findings regarding the role of methylation changes in hATTR, our study presents several limitations. Since hATTR is a rare disease, we investigated a relatively small sample size. Our calculation showed that the sample size investigated in our main analysis should provide > 80% statistical power to detect medium effect sizes (Δβ = 0.2; Additional file 7). However, large samples will be needed to investigate how epigenetic changes affect hATTR symptoms and differences across TTR amyloidogenic mutations. Our cohort showed age and sex differences between carriers of TTR amyloidogenic mutations and controls. Our analysis was adjusted for these confounding variables together with blood cell types, genetic principal components, and epigenetically determined smoking status. More balanced case–control groups are needed to investigate how epigenetic differences are associated differently between sexes and across age groups. We used transcriptomic data from GTEx project to explore the potential mechanisms related to the epigenetic associations identified. Further studies generating transcriptomic and epigenomic information across multiple informative tissues will provide a more comprehensive understanding of the molecular processes involved in hATTR. Although we provide some preliminary evidence showing that Val30Met mutation disrupts cg13139646 which appears co-methylated with CpG sites potentially associated with hATTR, our findings should be considered exploratory. Indeed, since Val30Met disrupts the probe extension for the cg13139646 methylation site, the impact of this mutation on the methylation changes in TTR gene should be confirmed using alternative methods.
Our study provided novel insights regarding hATTR pathogenesis, supporting the involvement of methylation changes in the amyloidogenic process induced by TTR disease-causing mutations. Further studies will be needed to characterize specific mechanisms underlying the epigenetic associations with particular attention to the potential role of amyloid-beta metabolic process and inflammatory response. The understanding of how methylation changes modulate the penetrance and the severity of TTR mutations could lead to the identification of novel targets to develop treatments and screening tools for the carriers. Additionally, similarly to what was observed with respect to genetic variation , it will be important to estimate the epigenetic similarity between hATTR and wild-type transthyretin amyloidosis.
Thirty-eight symptomatic patients and 10 asymptomatic TTR mutations carriers were recruited from three Italian centers for the treatment of systemic amyloidosis: “San Giovanni Calibita” Fatebenefratelli Hospital, Isola Tiberina—Rome, Fondazione Policlinico Universitario “A. Gemelli”—Rome, and Careggi University Hospital—Florence [16,17,18,19,20]. Thirty-two controls were recruited by the Department of Biology—University of Rome “Tor Vergata” (Table 2). hATTR diagnosis was based on the presence of clinical signs and symptoms and the presence of an amyloidogenic mutation on TTR gene. One hATTR patient is a carrier of a mutation (rs36204272) in an intronic region with a putative clinical impact . Carpal tunnel syndrome and cardiac involvement with confirmed TTR amyloid deposits are present in rs36204272 carrier. Information regarding the organ involvements was collected for each patient: peripheral and nerve involvement (nerve conduction study); cardiac involvement (electrocardiographic and echocardiography anomalies); gastrointestinal involvement (gastric paresis, stypsis, or diarrhea); autonomic neurological involvement (orthostatic hypotension and urinary incontinence); ocular involvement (vitreous opacities): and carpal tunnel syndrome (median nerve decompression) [11, 62,63,64]. In our previous analysis in this cohort , we observed that certain mutations were associated with specific clinical manifestations: Ile68Leu and Val122Ile with cardiac involvement; Val30Met with ocular symptoms; Phe64Leu with renal involvement. Conversely, autonomic neurological, peripheral neurological, and gastrointestinal symptoms were observed in most of the amyloidogenic mutations . The present study was approved under the protocol 39/18 by the Comitato Etico Indipendente, Fondazione Policlinico Tor Vergata—Rome, Italy.
DNA methylation analysis
DNA was extracted using the phenol/chloroform protocol  and purified through Amicon Ultra-0.5 mL Centrifugal Filters (EMD Millipore) to achieve a DNA concentration of 100 ng/µL. DNA concentration was checked via NanoDrop technology (ND-1000, Thermo Fisher Scientific) and Qubit Quantitation technology (High Accuracy & Sensitivity, Thermo Fisher). DNA methylation analysis was executed in two phases: the EZ DNA Methylation kit (Zymo Research) was used to perform sodium bisulfite conversion; the Illumina Infinium Methylation EPIC Chip (with over 850,000 methylation sites; Illumina Inc.) was used to quantify DNA methylation according to the standard Illumina protocol. The methylation array analysis was performed at the Connecting bio-research and Industry Center, Trieste—Italy.
Preprocessing, quality control, and normalization
The raw signal intensity files were processed and cleaned using R 3.6 and ChAMP package . The ratio of methylated and unmethylated intensities from idat files was converted into beta-values for further processing. The probes failing thresholds on detection p value (< 0.01) bead count, sites near SNPs—this method performed extensive characterization of probes on the EPIC and HM450 microarrays, including mappability to the latest genome build, genomic copy number of the 3΄ nested subsequence, and influence of polymorphisms including a previously unrecognized color channel switch for Type I probes, probes that align to multiple positions, sex chromosomes and outliers were removed. None of the samples failed quality control. The remaining 718,509 probes for 80 individuals were normalized with BMIQ. Batch effects were assessed using singular vector decomposition and corrected with ComBat method . The genomic lambda of the case–control association was 1.03, calculated using QQPerm package (https://cran.r-project.org/web/packages/QQperm/index.html).
Blood cell type composition, genetic variability estimation, and smoking prediction
A reference-based method was employed to adjust for the heterogeneity due to the cell type composition of the whole blood samples investigated . This method uses specific DNA methylation signatures derived from purified whole blood cell type as biomarkers of cell identity, to correct the beta-values. Cell proportions for five cell types (B cells, granulocytes, monocytes, natural killer cells, and T cells) were detected, and a linear regression was applied [66, 68]. To account for the genetic variability among the samples investigated, principal components (PCs) were calculate using the method proposed by Barfield, Almli . This approach allowed us to compute PCs based on CpGs selected for their proximity to SNPs. The data obtained can be used to adjust for population stratification in DNA methylation studies when genome-wide SNP data are unavailable . Cigarette smoke has a very large effect on DNA methylation profile, triggering alteration at multiple CpGs . Consequently, smoking status needs to be considered as a potential confounder in epigenetic association studies. EpiSmokEr package was used to classify the smoking status of each participant on the basis of their epigenetic profile . Briefly, EpiSmokEr is a prediction tool that provides smoking probabilities for each individual (never smoker, former smoker, and current smoker) using a set of 121 informative CpG sites .
We conducted an epigenome-wide analysis testing 718,509 methylation sites. First, we investigated the methylation changes (measured as M values, i.e., the log2 ratio of the intensities of methylated probe versus unmethylated probe; beta-values are plotted to provide more biologically interpretable data) between 48 carriers of a TTR amyloidogenic mutation) and 32 controls (i.e., non-carriers). Considering CpG sites that survive epigenome-wide multiple testing correction, we verified whether the associations observed were due to disease-associated genetic differences or post-disease processes, comparing (i) patients affected by hATTR vs. controls, (ii) asymptomatic carriers of TTR mutations vs. controls, and iii) patients affected by hATTR vs. asymptomatic carriers of TTR mutations. Additionally, we also tested whether the methylation changes observed in the case–control analysis were different between (i) Val30Met carriers vs. controls, (ii) carriers of other TTR mutations vs. controls and, iii) Val30Met patients vs. carriers of other TTR mutations.
To investigate the functionality of the CpG disrupted by Val30Met mutation, we analyzed its co-methylation with CpG sites in the surrounding region (NC_000018.9: 28,171,000–30,171,500). This region was selected based on the TTR regulatory mechanisms observed in previous studies [20,21,22]. We used cor() R function to calculate Pearson’s correlation coefficient testing 367 sites and considering the methylation levels (M values) using as reference the non-carriers. CpG sites with high co-methylation (Person’s r2 > 0.2) were investigated with respect to hATTR-related traits (i.e., carrier status, disease status, and symptoms). In all association analyses, we implemented a linear regression analysis including cell composition proportions, top three genetic PCs, epigenetically determined smoking status, age, and sex as covariates. The results of the regression analyses were reported as standardized regression coefficients and p values. FDR method  was applied to adjust the results for epigenome-wide testing and the q value < 0.05 was considered as the significance threshold. Co-expression analysis was conducted using GTEx v8  via the Multi-Gene Query available at https://www.gtexportal.org/. Ggplot2 R package  was employed to plot co-methylation pattern results. STRING v.11.0  was used to identify protein interaction with the loci identified, considering experiments, co-expression, co-occurrence, gene fusion, and neighborhood as active sources and an interaction score higher than 0.4 (medium confidence). The protein interaction network was investigated further conducting functional enrichments association related to the protein–protein interactions identified considering Gene Ontologies  for biological processes and molecular pathways available from Reactome Database  and Kyoto Encyclopedia of Genes and Genomes (KEGG) . FDR (q value < 0.05) was applied to account for multiple testing assuming the whole genome as the statistical background. A statistical power calculation was done using pwrEWAS tool  considering medium and small effect sizes (Δβ = 0.5 and 0.2, respectively) and multiple sample sizes.
Availability of data and materials
Data supporting the findings of this study are available within this article and its additional files.
Amyloid-beta precursor protein
Hereditary transthyretin amyloidosis
Epigenome-wide association studies
False discovery rate
FYN proto-oncogene, Src family tyrosine kinase
Immunoglobulin heavy variable 3–11
Palaninathan SK. Nearly 200 X-ray crystal structures of transthyretin: what do they tell us about this protein and the design of drugs for TTR amyloidoses? Curr Med Chem. 2012;19(15):2324–42.
Ueda M, Ando Y. Recent advances in transthyretin amyloidosis therapy. Transl Neurodegener. 2014;3:19.
Hellman U, Alarcon F, Lundgren HE, Suhr OB, Bonaiti-Pellie C, Plante-Bordeneuve V. Heterogeneity of penetrance in familial amyloid polyneuropathy, ATTR Val30Met, in the Swedish population. Amyloid. 2008;15(3):181–6.
Conceicao I. Clinical features of TTR-FAP in Portugal. Amyloid. 2012;19(Suppl 1):71–2.
Parman Y, Adams D, Obici L, Galan L, Guergueltcheva V, Suhr OB, et al. Sixty years of transthyretin familial amyloid polyneuropathy (TTR-FAP) in Europe: where are we now? A European network approach to defining the epidemiology and management patterns for TTR-FAP. Curr Opin Neurol. 2016;29(Suppl 1):S3–13.
Conceicao I, Damy T, Romero M, Galan L, Attarian S, Luigetti M, et al. Early diagnosis of ATTR amyloidosis through targeted follow-up of identified carriers of TTR gene mutations. Amyloid. 2019;26(1):3–9.
Ando Y, Coelho T, Berk JL, Cruz MW, Ericzon BG, Ikeda S, et al. Guideline of transthyretin-related hereditary amyloidosis for clinicians. Orphanet J Rare Dis. 2013;8:31.
Sousa A, Coelho T, Barros J, Sequeiros J. Genetic epidemiology of familial amyloidotic polyneuropathy (FAP)-type I in Povoa do Varzim and Vila do Conde (north of Portugal). Am J Med Genet. 1995;60(6):512–21.
Plante-Bordeneuve V, Carayol J, Ferreira A, Adams D, Clerget-Darpoux F, Misrahi M, et al. Genetic study of transthyretin amyloid neuropathies: carrier risks among French and Portuguese families. J Med Genet. 2003;40(11):e120.
Sousa A, Andersson R, Drugge U, Holmgren G, Sandgren O. Familial amyloidotic polyneuropathy in Sweden: geographical distribution, age of onset, and prevalence. Hum Hered. 1993;43(5):288–94.
Luigetti M, Conte A, Del Grande A, Bisogni G, Madia F, Lo Monaco M, et al. TTR-related amyloid neuropathy: clinical, electrophysiological and pathological findings in 15 unrelated patients. Neurol Sci. 2013;34(7):1057–63.
Alves-Ferreira M, Coelho T, Santos D, Sequeiros J, Alonso I, Sousa A, et al. A Trans-acting factor may modify age at onset in familial amyloid polyneuropathy ATTRV30M in Portugal. Mol Neurobiol. 2018;55(5):3676–83.
Bonaiti B, Olsson M, Hellman U, Suhr O, Bonaiti-Pellie C, Plante-Bordeneuve V. TTR familial amyloid polyneuropathy: does a mitochondrial polymorphism entirely explain the parent-of-origin difference in penetrance? Eur J Hum Genet. 2010;18(8):948–52.
Santos D, Coelho T, Alves-Ferreira M, Sequeiros J, Mendonca D, Alonso I, et al. Variants in RBP4 and AR genes modulate age at onset in familial amyloid polyneuropathy (FAP ATTRV30M). Eur J Hum Genet. 2016;24(5):756–60.
Soares ML, Coelho T, Sousa A, Batalov S, Conceicao I, Sales-Luis ML, et al. Susceptibility and modifier genes in Portuguese transthyretin V30M amyloid polyneuropathy: complexity in a single-gene disease. Hum Mol Genet. 2005;14(4):543–53.
Iorio A, De Angelis F, Di Girolamo M, Luigetti M, Pradotto L, Mauro A, et al. Most recent common ancestor of TTR Val30Met mutation in Italian population and its potential role in genotype-phenotype correlation. Amyloid. 2015;22(2):73–8.
Polimanti R, Di Girolamo M, Manfellotto D, Fuciarelli M. Functional variation of the transthyretin gene among human populations and its correlation with amyloidosis phenotypes. Amyloid. 2013;20(4):256–62.
Polimanti R, Di Girolamo M, Manfellotto D, Fuciarelli M. In silico analysis of TTR gene (coding and non-coding regions, and interactive network) and its implications in transthyretin-related amyloidosis. Amyloid. 2014;21(3):154–62.
Iorio A, De Angelis F, Di Girolamo M, Luigetti M, Pradotto LG, Mazzeo A, et al. Population diversity of the genetically determined TTR expression in human tissues and its implications in TTR amyloidosis. BMC Genomics. 2017;18(1):254.
Iorio A, De Lillo A, De Angelis F, Di Girolamo M, Luigetti M, Sabatelli M, et al. Non-coding variants contribute to the clinical heterogeneity of TTR amyloidosis. Eur J Hum Genet. 2017;25(9):1055–60.
De Lillo A, De Angelis F, Di Girolamo M, Luigetti M, Frusconi S, Manfellotto D, et al. Phenome-wide association study of TTR and RBP4 genes in 361,194 individuals reveals novel insights in the genetics of hereditary and wildtype transthyretin amyloidoses. Hum Genet. 2019;138(11–12):1331–40.
Polimanti R, Nunez YZ, Gelernter J. Increased risk of multiple outpatient surgeries in African-American Carriers of Transthyretin Val122Ile mutation is modulated by non-coding variants. J Clin Med. 2019;8(2).
Feinberg AP. The key role of epigenetics in human disease prevention and mitigation. N Engl J Med. 2018;378(14):1323–34.
Yan H, Tian S, Slager SL, Sun Z, Ordog T. Genome-wide epigenetic studies in human disease: a primer on-omic technologies. Am J Epidemiol. 2016;183(2):96–109.
Szymczak S, Dose J, Torres GG, Heinsen FA, Venkatesh G, Datlinger P, et al. DNA methylation QTL analysis identifies new regulators of human longevity. Hum Mol Genet. 2020.
Clissold RL, Ashfield B, Burrage J, Hannon E, Bingham C, Mill J, et al. Genome-wide methylomic analysis in individuals with HNF1B intragenic mutation and 17q12 microdeletion. Clin Epigenet. 2018;10(1):97.
Guay SP, Brisson D, Munger J, Lamarche B, Gaudet D, Bouchard L. ABCA1 gene promoter DNA methylation is associated with HDL particle profile and coronary artery disease in familial hypercholesterolemia. Epigenetics. 2012;7(5):464–72.
Magalhaes M, Rivals I, Claustres M, Varilh J, Thomasset M, Bergougnoux A, et al. DNA methylation at modifier genes of lung disease severity is altered in cystic fibrosis. Clin Epigenet. 2017;9:19.
Aguet F, Barbeira AN, Bonazzola R, Brown A, Castel SE, Jo B, et al. The GTEx Consortium atlas of genetic regulatory effects across human tissues. bioRxiv. 2019:787903.
Luigetti M, Romano A, Di Paolantonio A, Bisogni G, Sabatelli M. Diagnosis and treatment of hereditary transthyretin amyloidosis (hATTR) polyneuropathy: current perspectives on improving patient care. Ther Clin Risk Manag. 2020;16:109–23.
Zheng H, Koo EH. Biology and pathophysiology of the amyloid precursor protein. Mol Neurodegener. 2011;6(1):27.
Sassi C, Ridge PG, Nalls MA, Gibbs R, Ding J, Lupton MK, et al. Influence of coding variability in APP-abeta metabolism genes in sporadic Alzheimer’s disease. PLoS ONE. 2016;11(6):e0150079.
Wang Z, Xu Q, Cai F, Liu X, Wu Y, Song W. BACE2, a conditional beta-secretase, contributes to Alzheimer's disease pathogenesis. JCI Insight. 2019;4(1).
Voytyuk I, Mueller SA, Herber J, Snellinx A, Moechars D, van Loo G, et al. BACE2 distribution in major brain cell types and identification of novel substrates. Life Sci Alliance. 2018;1(1):e201800026.
Zhang YW, Thompson R, Zhang H, Xu H. APP processing in Alzheimer’s disease. Mol Brain. 2011;4:3.
Hartl D, May P, Gu W, Mayhaus M, Pichler S, Spaniol C, et al. A rare loss-of-function variant of ADAM17 is associated with late-onset familial Alzheimer disease. Mol Psychiatry. 2020;25(3):629–39.
Suh J, Choi SH, Romano DM, Gannon MA, Lesinski AN, Kim DY, et al. ADAM10 missense mutations potentiate beta-amyloid accumulation by impairing prodomain chaperone function. Neuron. 2013;80(2):385–401.
Tambini MD, Norris KA, D'Adamio L. Opposite changes in APP processing and human Abeta levels in rats carrying either a protective or a pathogenic APP mutation. Elife. 2020;9.
Durrant CS, Ruscher K, Sheppard O, Coleman MP, Ozen I. Beta secretase 1-dependent amyloid precursor protein processing promotes excessive vascular sprouting through NOTCH3 signalling. Cell Death Dis. 2020;11(2):98.
Alemi M, Silva SC, Santana I, Cardoso I. Transthyretin stability is critical in assisting beta amyloid clearance: relevance of transthyretin stabilization in Alzheimer’s disease. CNS Neurosci Ther. 2017;23(7):605–19.
Costa R, Goncalves A, Saraiva MJ, Cardoso I. Transthyretin binding to A-Beta peptide–impact on A-Beta fibrillogenesis and toxicity. FEBS Lett. 2008;582(6):936–42.
Du J, Murphy RM. Characterization of the interaction of beta-amyloid with transthyretin monomers and tetramers. Biochemistry. 2010;49(38):8276–89.
Sakai K, Asakawa M, Takahashi R, Ishida C, Nakamura R, Hamaguchi T, et al. Coexistence of transthyretin- and Abeta-type cerebral amyloid angiopathy in a patient with hereditary transthyretin V30M amyloidosis. J Neurol Sci. 2017;381:144–6.
Xiang Q, Bi R, Xu M, Zhang DF, Tan L, Zhang C, et al. Rare genetic variants of the transthyretin gene are associated with Alzheimer’s disease in Han Chinese. Mol Neurobiol. 2017;54(7):5192–200.
Choi SH, Leight SN, Lee VM, Li T, Wong PC, Johnson JA, et al. Accelerated Abeta deposition in APPswe/PS1deltaE9 mice with hemizygous deletions of TTR (transthyretin). J Neurosci. 2007;27(26):7006–10.
Kerridge C, Belyaev ND, Nalivaeva NN, Turner AJ. The Abeta-clearance protein transthyretin, like neprilysin, is epigenetically regulated by the amyloid precursor protein intracellular domain. J Neurochem. 2014;130(3):419–31.
Azevedo EP, Guimaraes-Costa AB, Bandeira-Melo C, Chimelli L, Waddington-Cruz M, Saraiva EM, et al. Inflammatory profiling of patients with familial amyloid polyneuropathy. BMC Neurol. 2019;19(1):146.
Gertz J, Varley KE, Reddy TE, Bowling KM, Pauli F, Parker SL, et al. Analysis of DNA methylation in a three-generation family reveals widespread genetic influence on epigenetic regulation. PLoS Genet. 2011;7(8):e1002228.
Hellman A, Chess A. Extensive sequence-influenced DNA methylation polymorphism in the human genome. Epigenet Chromatin. 2010;3(1):11.
Affinito O, Palumbo D, Fierro A, Cuomo M, De Riso G, Monticelli A, et al. Nucleotide distance influences co-methylation between nearby CpG sites. Genomics. 2020;112(1):144–50.
Mayo L, Trauger SA, Blain M, Nadeau M, Patel B, Alvarez JI, et al. Regulation of astrocyte activation by glycolipids drives chronic CNS inflammation. Nat Med. 2014;20(10):1147–56.
Ohno S. The genetic background of arrhythmogenic right ventricular cardiomyopathy. J Arrhythm. 2016;32(5):398–403.
Brodehl A, Belke DD, Garnett L, Martens K, Abdelfatah N, Rodriguez M, et al. Transgenic mice overexpressing desmocollin-2 (DSC2) develop cardiomyopathy associated with myocardial inflammation and fibrotic remodeling. PLoS ONE. 2017;12(3):e0174019.
Cheung CC, Roston TM, Andrade JG, Bennett MT, Davis MK. Arrhythmias in cardiac amyloidosis: challenges in risk stratification and treatment. Can J Cardiol. 2020;36(3):416–23.
Macedo AVS, Schwartzmann PV, de Gusmao BM, Melo MDT, Coelho-Filho OR. Advances in the treatment of cardiac amyloidosis. Curr Treat Options Oncol. 2020;21(5):36.
Henein MY, Suhr OB, Arvidsson S, Pilebro B, Westermark P, Hornsten R, et al. Reduced left atrial myocardial deformation irrespective of cavity size: a potential cause for atrial arrhythmia in hereditary transthyretin amyloidosis. Amyloid. 2018;25(1):46–53.
Prestes PR, Marques FZ, Lopez-Campos G, Lewandowski P, Delbridge LMD, Charchar FJ, et al. Involvement of human monogenic cardiomyopathy genes in experimental polygenic cardiac hypertrophy. Physiol Genomics. 2018;50(9):680–7.
Jjingo D, Conley AB, Yi SV, Lunyak VV, Jordan IK. On the presence and role of human gene-body DNA methylation. Oncotarget. 2012;3(4):462–74.
Meng H, Cao Y, Qin J, Song X, Zhang Q, Shi Y, et al. DNA methylation, its mediators and genome integrity. Int J Biol Sci. 2015;11(5):604–17.
Ayub M, Basit S, Jelani M, Ur Rehman F, Iqbal M, Yasinzai M, et al. A homozygous nonsense mutation in the human desmocollin-3 (DSC3) gene underlies hereditary hypotrichosis and recurrent skin vesicles. Am J Hum Genet. 2009;85(4):515–20.
Levine TD, Bland RJ. Incidence of nonamyloidogenic mutations in the transthyretin gene in patients with autonomic and small fiber neuropathy. Muscle Nerve. 2018;57(1):140–2.
Russo M, Mazzeo A, Stancanelli C, Di Leo R, Gentile L, Di Bella G, et al. Transthyretin-related familial amyloidotic polyneuropathy: description of a cohort of patients with Leu64 mutation and late onset. J Peripher Nerv Syst. 2012;17(4):385–90.
Mazzeo A, Russo M, Di Bella G, Minutoli F, Stancanelli C, Gentile L, et al. Transthyretin-related familial amyloid polyneuropathy (TTR-FAP): a single-center experience in Sicily, an Italian endemic area. J Neuromuscul Dis. 2015;2(s2):S39–48.
Gagliardi C, Perfetto F, Lorenzini M, Ferlini A, Salvi F, Milandri A, et al. Phenotypic profile of Ile68Leu transthyretin amyloidosis: an underdiagnosed cause of heart failure. Eur J Heart Fail. 2018;20(10):1417–25.
Budowle B, Moretti, T., Smith, J., and Dizinno, J. DNA typing protocols: molecular biology and forensic analysis. A Bio Techniques® Books Publications (Natick, MA: Eaton Publishing). 2000.
Tian Y, Morris TJ, Webster AP, Yang Z, Beck S, Feber A, et al. ChAMP: updated methylation analysis pipeline for Illumina BeadChips. Bioinformatics. 2017;33(24):3982–4.
Johnson WE, Li C, Rabinovic A. Adjusting batch effects in microarray expression data using empirical Bayes methods. Biostatistics. 2007;8(1):118–27.
Houseman EA, Accomando WP, Koestler DC, Christensen BC, Marsit CJ, Nelson HH, et al. DNA methylation arrays as surrogate measures of cell mixture distribution. BMC Bioinform. 2012;13:86.
Barfield RT, Almli LM, Kilaru V, Smith AK, Mercer KB, Duncan R, et al. Accounting for population stratification in DNA methylation studies. Genet Epidemiol. 2014;38(3):231–41.
Zeilinger S, Kuhnel B, Klopp N, Baurecht H, Kleinschmidt A, Gieger C, et al. Tobacco smoking leads to extensive genome-wide changes in DNA methylation. PLoS ONE. 2013;8(5):e63812.
Bollepalli S, Korhonen T, Kaprio J, Anders S, Ollikainen M. EpiSmokEr: a robust classifier to determine smoking status from DNA methylation data. Epigenomics. 2019;11(13):1469–86.
Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J Roy Stat Soc: Ser B (Methodol). 1995;57:289–300.
Wickham H. ggplot2: Elegant graphics for data analysis. Springer, New York. 2016.
Szklarczyk D, Gable AL, Lyon D, Junge A, Wyder S, Huerta-Cepas J, et al. STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 2019;47(D1):D607–13.
The Gene Ontology C. The gene ontology resource: 20 years and still GOing strong. Nucleic Acids Res. 2019;47(D1):D330-D8.
Jassal B, Matthews L, Viteri G, Gong C, Lorente P, Fabregat A, et al. The reactome pathway knowledgebase. Nucleic Acids Res. 2020;48(D1):D498–503.
Kanehisa M, Furumichi M, Tanabe M, Sato Y, Morishima K. KEGG: new perspectives on genomes, pathways, diseases and drugs. Nucleic Acids Res. 2017;45(D1):D353–61.
Graw S, Henn R, Thompson JA, Koestler DC. pwrEWAS: a user-friendly tool for comprehensive power estimation for epigenome wide association studies (EWAS). BMC Bioinform. 2019;20(1):218.
We thank the participants involved in this study and their caregivers.
This study was supported by an Investigator-Initiated Research from Pfizer Inc. to the University of Rome Tor Vergata. Pfizer Inc. had no role in the study design, data analysis, and data interpretation of the present study.
Ethics approval and consent to participate
This study was approved under the protocol 39/18 by the Comitato Etico Indipendente, Fondazione Policlinico Tor Vergata—Rome, Italy. Informed consent was obtained from each participant involved.
Drs. Fuciarelli and Polimanti received research grants from Pfizer Inc. to conduct epigenetic studies of hATTR. The other authors reported no biomedical financial interests or potential conflicts of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
. Methylation levels (M values) of cg09097335 site in carriers (cases) vs. non-carriers (controls) of amyloidogenic mutations. The regression line is shown in blue.
. Methylation change of cg09097335 site (upper panel: beta values; lower panel: M values) between i) hATTR patients (carriers of TTR amyloidogenic mutations with hATTR diagnosis) vs. controls, ii) asymptomatic carriers vs. controls, iii) V30M carriers vs. controls, and iv) carriers of other TTR mutation vs. controls. Standardized regression coefficient and p value reported for each comparison are derived from the analysis conducted on the M values.
. Methylation change of cg13139646 site between V30M carriers vs. controls (upper panel: beta values; bottom panel: M values). Standardized regression coefficient and p value reported for each comparison are derived from the analysis conducted on the M values.
. Co-methylation analysis (Pearson’s correlation based on M values) with respect to cg13139646 (red text).
. Association of cg13139646 co-methylated CpG sites (Pearson’s r2 > 0.2) with i) carrier status (cases vs. controls), ii) disease status (cases vs. asymptomatic) and, iii) symptoms (AGE, ANS, CTS, EYE, GI, HEA, PNS). Information about cg probe (cgID), chromosome localization (CHR), position (POS), determination coefficient (Pearson’s r2), standardized regression coefficients (Reg. Coef.), p value (Pval), mapped gene (GENE), gene region (Genomic Region), are reported. AGE = age of onset; ANS = autonomic nervous system; CTS = carpal tunnel syndrome; EYE = ocular involvement; GI = gastrointestinal involvement; HEA = cardiac involvement; PNS = peripheral nervous system. Red and underlined text = significant results surviving Bonferroni correction. Red text = nominally significant results.
. Methylation changes of co-methylated cg13139646-correlated CpG sites with respect to hATTR-related phenotypes: cg19203115, hATTR patients vs. asymptomatic carriers; cg11481443, cases vs. controls; cg02936398, cases vs. controls; cg14311811, hATTR patients vs. asymptomatic carriers; cg27392998, cardiac involvement in hATTR patients; cg18038361, cardiac involvement in hATTR patients; cg16492377, carpal tunnel syndrome in hATTR patients; cg14719951, peripheral nervous system involvement in hATTR patients.
. Statistical power calculations based on medium and small effect sizes (Δβ = 0.5 and 0.2, respectively) and multiple sample sizes.
About this article
Cite this article
De Lillo, A., Pathak, G.A., De Angelis, F. et al. Epigenetic profiling of Italian patients identified methylation sites associated with hereditary transthyretin amyloidosis. Clin Epigenet 12, 176 (2020). https://doi.org/10.1186/s13148-020-00967-6