Next Article in Journal
Toward the Pathogenicity of the SLC26A4 p.C565Y Variant Using a Genetically Driven Mouse Model
Next Article in Special Issue
BET Protein-Mediated Transcriptional Regulation in Heart Failure
Previous Article in Journal
Evolution as a Guide to Designing xeno Amino Acid Alphabets
Previous Article in Special Issue
Epigenetic Changes Governing Scn5a Expression in Denervated Skeletal Muscle
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 22
Issue 6
10.3390/ijms22062790
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Differential Methylation in the GSTT1 Regulatory Region in Sudden Unexplained Death and Sudden Unexpected Death in Epilepsy

by
Steffan Noe Christiansen
1,†,
Stine Bøttcher Jacobsen
1,*,†,
Jeppe Dyrberg Andersen
1,
Marie-Louise Kampmann
1,
Linea Christine Trudsø
2,
Kristine Boisen Olsen
3,
Jacob Tfelt-Hansen
3,4,
Jytte Banner
3 and
Niels Morling
1
1
Section of Forensic Genetics, Department of Forensic Medicine, Faculty of Health and Medical Sciences, University of Copenhagen, DK-1353 Copenhagen, Denmark
2
Department of Biomedical Sciences, Faculty of Health and Medical Sciences, University of Copenhagen, DK-1353 Copenhagen, Denmark
3
Section of Forensic Pathology, Department of Forensic Medicine, Faculty of Health and Medical Sciences, University of Copenhagen, DK-1353 Copenhagen, Denmark
4
The Heart Centre, Department of Cardiology, Copenhagen University Hospital Rigshospitalet, DK-1353 Copenhagen, Denmark
*
Author to whom correspondence should be addressed.
†
These authors contributed equally.
Int. J. Mol. Sci. 2021, 22(6), 2790; https://0-doi-org.brum.beds.ac.uk/10.3390/ijms22062790
Submission received: 30 January 2021 / Revised: 27 February 2021 / Accepted: 4 March 2021 / Published: 10 March 2021
(This article belongs to the Special Issue Epigenetic Mechanisms of Cardiac Disease)
Browse Figures
Versions Notes

Abstract

:
Sudden cardiac death (SCD) is a diagnostic challenge in forensic medicine. In a relatively large proportion of the SCDs, the deaths remain unexplained after autopsy. This challenge is likely caused by unknown disease mechanisms. Changes in DNA methylation have been associated with several heart diseases, but the role of DNA methylation in SCD is unknown. In this study, we investigated DNA methylation in two SCD subtypes, sudden unexplained death (SUD) and sudden unexpected death in epilepsy (SUDEP). We assessed DNA methylation of more than 850,000 positions in cardiac tissue from nine SUD and 14 SUDEP cases using the Illumina Infinium MethylationEPIC BeadChip. In total, six differently methylated regions (DMRs) between the SUD and SUDEP cases were identified. The DMRs were located in proximity to or overlapping genes encoding proteins that are a part of the glutathione S-transferase (GST) superfamily. Whole genome sequencing (WGS) showed that the DNA methylation alterations were not caused by genetic changes, while whole transcriptome sequencing (WTS) showed that DNA methylation was associated with expression levels of the GSTT1 gene. In conclusion, our results indicate that cardiac DNA methylation is similar in SUD and SUDEP, but with regional differential methylation in proximity to GST genes.

1. Introduction

The sudden and unexpected death of a young individual is a devastating event in families and remains a key challenge to healthcare systems and forensic examinations worldwide. In some cases, death is the first disease manifestation, and this demonstrates that timely diagnosis is essential to prevent premature death [1,2]. The majority of sudden cardiac deaths (SCDs) in the young are suspected to result from undiagnosed structural cardiovascular disorders that can be identified at post-mortem investigations [3,4]. However, in cases of a structurally normal heart, routine autopsies often fail to identify a disease as the cause of death. These sudden unexplained deaths (SUDs) are presumed to be caused by undiagnosed cardiac arrhythmias and are, in the case of a negative histological and toxicological investigation, termed sudden arrhythmic death syndrome [5,6]. Diseases most commonly associated with SUD are inherited cardiac channelopathies such as long QT syndrome (LQTS), Brugada Syndrome (BrS), and catecholaminergic polymorphic ventricular tachycardia (CPVT) [7,8]. Extensive research on the genetics of inherited cardiac diseases have been performed and several studies report successful identification of causal genetic variants in relation to SUD [9,10,11]. Furthermore, an association between methylation of DNA and cardiac disease has been suggested [12,13,14,15]. However, the role of epigenetic patterns in SUD remains largely undetermined.
Epilepsy, diabetes, and psychiatric diseases have also been associated with a higher risk of SCD in young individuals [16,17,18]. Sudden unexpected death in epilepsy (SUDEP) accounts for the majority of premature mortality among patients with epilepsy [19,20]. SUD and SUDEP share common circumstances of death with both events occurring predominantly during sleep or at rest [21,22,23]. This may indicate similar pathophysiological mechanisms. Despite extensive research, the etiology of SUDEP remains largely unknown. Several risk factors have been suggested with the most predominant being the presence of generalised tonic-clonic seizures [24,25]. However, a combination of neurological, cardiac, and respiratory factors has also been proposed to contribute to death [22,26]. Cardiac channelopathies may predispose to SUDEP as channelopathies can present with both arrhythmia and seizure [27,28,29].
A genetic basis of SUDEP has been suggested, and several studies have identified pathogenic variants in genes co-expressed in the brain and the heart [30,31,32,33,34]. However, the potential contribution of epigenetic factors to the risk of SUDEP as well as the effect of gene expression levels remain largely unknown.
Epigenetic patterns are tissue specific and, hence, studies of epigenetic phenomena in various tissues require access to the tissue of interest. A major limitation of epigenetic studies of heart diseases is, therefore, the general unavailability of cardiac tissue. However, in post-mortem investigations, the accessibility of tissues from autopsied individuals offers a unique opportunity for epigenetic investigations of cardiac tissue. A key challenge of post-mortem genetic and epigenetic testing is, though, the post-mortem degradation of the tissue. Post-mortem decomposition may result in widespread fragmentation of nucleic acids [35]. The quality of the nucleic acids can, furthermore, be reduced by the tissue storage methods. Formalin fixed and paraffin embedded (FFPE) tissue has for many years been widely used as preparation for histological examinations and long-term storage. Unfortunately, FFPE tissue may induce sequence artefacts [36,37].
In this study, we investigated DNA methylation in FFPE cardiac tissues from SUD and SUDEP cases and identified differentially methylated regions (DMRs) between the two groups. The heart tissue from SUDEP victims was also examined for RNA expression. The results were compared to DNA variation investigated using whole genome sequencing (WGS).

2. Results

2.1. Characteristics of the Study Population

The study included ten males and four females diagnosed with SUDEP and six males and three females diagnosed with SUD (Table 1). The median age was 35 years in both groups with the ranges 0–50 years and 3–49 years among the SUD and SUDEP cases, respectively. Deaths in both groups were predominantly unwitnessed and occurred during sleep or rest. There were no significant differences between the two investigated groups except for the difference in the expected detected antiepileptic drugs (Table 1).

2.2. Allele Frequencies in Sudden Unexplained Death and Sudden Unexpected Death in Epilepsy

The allele frequencies of variants in SUD and SUDEP cases were assessed. No statistically significant difference between SUD and SUDEP was identified.

2.3. Principal Component Analysis

The principal component analysis (PCA) was conducted to investigate the main patterns in the DNA methylation data (Figure 1). The analysis revealed that PC2 tended to separate older individuals from younger individuals. However, the analysis revealed no overall difference between SUDEP and SUD.

2.4. Regulatory Association of the GSTT1 DMR

One DMR between the SUD and SUDEP cases was identified using the recommended default kernel-smoothed false discovery rate (FDR) threshold of DMRcate (version 2.2.3) [38]. The DMR was located in proximity to the GSTT1 gene (Figure 2A). To evaluate the regulatory relevance of the GSTT1 DMR, DNase I sensitivity data from the Encyclopedia of DNA Elements (ENCODE) project was used for the mapping of likely cis-regulatory elements. Three clusters overlapped the GSTT1 DMR with the second cluster scoring the maximal value for DNase I sensitivity. The expression data was also used to evaluate the regulatory relevance of the GSTT1 DMR. The Spearman’s correlation coefficients among the methylation of each single CpG position within the DMRs and the median Trimmed Mean of M-values (TMM) per kb values were calculated for all CpG positions within the DMR. Among the 13 CpG positions overlapping the DMR in the GSTT1 promoter region, five were negatively correlated (p-value ≤ 0.05) with the median TMM per kb for GSTT1 (Figure 2B).

2.5. Differentially Methylated Regions between Sudden Unexplained Death and Sudden Unexpected Death in Epilepsy

Since only one DMR was detected, we decided to increase the FDR threshold to 0.001 well aware that it would increase the risk of type I errors. Six DMRs were identified between the SUD and SUDEP using the less stringent thresholds for DMR identification (Table 2). The Δβ among the CpG positions ranged from 0.10–0.19 within the six DMRs. Five of the DMRs overlapped a gene (GSTT1, HLX, GSTM1, LOC391322, or TRIM5) while the remaining DMR (chr22:24,348,549-24,348,715) was located in proximity to the GSTT4 gene. Three of the DMRs (chr22:24,384,105-24,384,944, chr22:24,373,322-24,373,618, and chr22:24,348,549-24,348,715) were located within the same genomic region (<40 kb) as the genes, GSTT4, LOC391322, GSTT1, and the GSTT1 antisense RNA 1 (GSTT1-AS1).

3. Discussion

SUD and SUDEP have similar circumstances of death and remain key diagnostic challenges in autopsies. Both of them are diagnoses of exclusion of all other possible causes of death, and the phenotypes of the deceased individuals are often heterogeneous and vague [39]. In this study, we aimed to assess and compare genetic and epigenetic patterns in cardiac tissues from SUD and SUDEP cases to highlight differences between the two groups.
The methylation patterns of consecutive CpG sites in SUD and SUDEP cases were assessed to identify DMRs. We show that a region upstream the GSTT1 gene was differentially methylated between the SUD and SUDEP cases. By incorporating DNase I sensitivity and gene expression data, we found that the DMR is associated with GSTT1 gene expression. The DNA methylation was negatively correlated with gene expression suggesting a repressive function of DNA methylation on GSTT1 gene expression.
Only six DMRs were observed between the SUD and SUDEP cases. Five of these failed to pass the pre-determined FDR threshold. We performed no detailed analyses of these DMRs due to the relatively high Min. FDR for the regions (Table 2). However, two interesting patterns among the DMRs were observed. First, we observed that three DMRs were located in a genomic region from chr22:24,348,550-24,384,944 constituted by less than 40 kb. The region overlapped the glutathione transferases, GSTT4, GSTT1, and the GSTT1 antisense RNA 1 (GSTT1-AS1), and thus it could be speculated that these DMRs are a part of the same biological feature. However, additional research is required to confirm this. Second, among the remaining three DMRs outside the 40 kb region, one of the DMRs was located in proximity to GSTM1, which encodes a protein that is part of the glutathione S-transferase (GST) superfamily just as GSTT1 and GSTT4.
GSTs conjugate glutathione to toxic compounds such as products from oxidative metabolism, therapeutic drugs, and electrophilic molecules [40]. The normal function of GSTs is, therefore, essential for cellular detoxification and protection against oxidative stress. We cannot exclude that the observed differential methylation is a response to the intake of antiepileptic drugs in the SUDEP cases. However, a current study suggests an association between GSTT1 dysregulation and dilated cardiomyopathy [41]. Heritable cardiomyopathies with a mild phenotype may cause SUD in victims where the autopsy findings are insufficient to determine a death cause [42]. Haas and colleagues reported decreased expression of Gstt1 and Gstt2 in a mouse model following transverse aortic constriction-induced heart failure and suggest that these genes are involved in adaptive or maladaptive cascades of the failing heart [41]. To our knowledge, there are no other direct associations between GSTT1 and cardiac diseases associated with SUD. However, members of the GST superfamily have been associated with inhibition of the ryanodine receptor 2 [43]. The gene encoding this protein, RYR2, is one of the most commonly mutated genes in SUD [9,10], and variants within this gene are likely increasing the risk of SUDEP [30,44,45].
In addition to the association with cardiac disease, GSTT1 null variants have been associated with epilepsy [46,47]. Metabolic dysregulation and oxidative stress is recognised as both a cause and consequence of epileptic seizures and may explain the relevance of GSTT1 in SUDEP cases [48]. Based on the current knowledge about the GST superfamily, it cannot be excluded that the DMRs are associated with epilepsy and heart disease.
No statistically significant difference in allele frequencies was identified between the SUD and SUDEP cases. This may be explained by the heterogeneous nature of the diseases associated with SUD and SUDEP and the lack of statistical power due to the relatively small study population. Even though the aetiology of SUDEP is largely undetermined, the underlying pathophysiology is suspected to be multifactorial and complex [49,50,51,52]. SUD is likely caused by cardiac channelopathies such as LQTS, BrS, and CPVT, all of which are suspected to result from rare genetic variants [9,53]. It was, therefore, unlikely to identify a common causal genetic variant within the SUD and SUDEP cohorts. We did not identify genetic variation in proximity to the GSTT1 gene that could explain the differential methylation between SUD and SUDEP cases. Neither did we find a correlation between genotypes in the GSTT1 region and GSTT1 RNA levels.
By using a larger study population, the power to detect significant differences in allele frequencies and differential methylation between the two groups would be higher. However, since both SUD and especially SUDEP are rare events, it is very difficult to obtain large study populations. In a 7-year period (2000–2006), only 50 definite and probable SUDEPs in the age group 1–35 years were identified among all residents in Denmark [16]. We observed that a part of the variation in methylation patterns was likely explained by age. This shows the importance of incorporating age as covariable in the statistical model for the identification of DMRs. We also investigated if the degree of post-mortem decomposition of the cardiac tissue could explain a part of the variation in methylation patterns. We did not find an association between the overall methylation patterns and the degree of post-mortem decomposition (data not shown). When using a relatively small sample size, it is difficult to obtain statistical significance after correction for multiple testing. Small methylation differences of single positions are therefore difficult to detect, but if the changes are consistent among a region, the statistical power will be greater [54]. Thus, we examined DMRs, which are arguably functionally more important compared to differentially methylated positions [55]. The use of FFPE tissue limits the quality of several molecular analyses due to the degradation of nucleic acids. However, array-based DNA methylation measurements of frozen and FFPE tissue show a good correlation [56,57].
We did not have the opportunity to compare our results with cardiac tissues from deceased individuals that died due to a non-cardiac event. Thus, we cannot determine if both SUD and SUDEP are differentially methylated compared to other deceased individuals.
With the use of a unique collection of cardiac tissue, it was demonstrated that the DNA methylation profiles of SUD and SUDEP were similar overall, except for the DMRs identified in proximity to genes encoding GSTs. Future studies are, therefore, needed to elucidate the role of DNA methylation in assumed cardiac deaths.

4. Materials and Methods

4.1. Study Population

The tissues were collected from autopsies of SUD and SUDEP cases that were carried out in 2009–2011 at the Section of Forensic Pathology, Department of Forensic Medicine, Faculty of Health and Medical Sciences, University of Copenhagen, Denmark. The cardiac tissues were fixed in 4% buffered formaldehyde for 48 hours and embedded in paraffin. Whole blood was stored at −20 °C in tubes coated with EDTA. Differences in clinical characteristics, such as age, sex, and circumstances of death, between SUD and SUDEP were assessed using the Fisher’s exact test for categorical variables and the Mann–Whitney U test for continuous variables.

4.2. Whole Genome Sequencing

WGS was performed as previously described [58]. Briefly, spin column-based DNA extraction from whole blood was performed, followed by WGS library preparation, and subsequent paired-end sequencing (2 × 150 bp) on a NextSeq 500 using the High Output kit v2 (Illumina, San Diego, CA, USA). Quality control of sequencing data and variant detection was performed using an in-house pipeline described in detail in [58].
Differences in allele frequencies between SUD and SUDEP cases were assessed with a two-sided Fisher’s exact test using the fisher.test function in the statistical software R (version 3.6.0) [59]. Statistically significant variants were defined as variants with a Bonferroni corrected p-value < 0.05.

4.3. Whole Transcriptome Sequencing

RNA was extracted from FFPE tissue of the anterior wall of the left cardiac ventricle. Two 20 µm FFPE tissue slides were used for each extraction. Excess paraffin was removed using xylene followed by a wash with 96–100% ethanol. Total RNA extraction was performed using the RNeasy FFPE Kit (Qiagen, Hilden, Germany) following manufacturer’s recommendations. Whole transcriptome library preparation was performed using the SMARTer® Stranded Total RNA-Seq Kit v2-Pico Input Mammalian (Takara Bio Europe, St-Germain-en-Laye, France). Due to the fragmented nature of RNA extracted from FFPE tissue, no further fragmentation of the RNA was performed. cDNA synthesis was performed using random N6 primers followed by the addition of adapters and depletion of ribosomal cDNA. The quality and quantity of cDNA libraries were assessed using the Bioanalyzer 2100 instrument (Agilent Technologies, Santa Clara, CA, USA) and Qubit dsDNA HS Assay (Invitrogen, Waltham, MA, USA), respectively. Single-read sequencing (1 × 75 bp) on a NextSeq 500 using the High Output kit v2 (Illumina, San Diego, CA, USA) was performed. Biological triplicates were investigated for each of the 14 SUDEP cases. Whole transcriptome sequencing (WTS) data was only investigated in the SUDEP victims.

4.3.1. Preprocessing and Alignment of RNA Data

The BCL files were converted to FASTQ files using bcl2fastq version 2.17.1.14 (Illumina, San Diego, CA, USA), and the sequencing quality was assessed using FastQC version 0.11.2 [60]. Consecutive stretches of low-quality bases (Q < 30) at the 5′ and 3′ termini as well as adapters were removed using AdapterRemoval version 2.1.3 [61]. Sequencing reads shorter than 20 bp were discarded.
The reads were aligned using STAR version 2.5.3a [62] using default settings except for two parameters. Reads aligned to multiple loci were classified as unmapped (--outFilterMultimapNmax 1), and the unannotated splice junction prediction was disabled (--alignIntronMax 1). The gene counts were obtained with STAR based on the human genome assembly GRCh37 (hg19) FASTA and GTF files from the GENCODE consortium release 27 [63].

4.3.2. Normalisation of Read Counts

Reads in annotated non-coding, ribosomal, mitochondrial, or nuclear mitochondrial DNA segment (NUMT) sequences were excluded from analysis as previously described [64]. Normalisation of RNA data was performed according to the TMM method [65] using EdgeR [66] in R (version 3.6.0). Additionally, the read counts were normalised according to the gene length, and the normalised gene counts were presented as TMM per kilobase (kb).

4.4. Quantification of DNA Methylation

DNA from the FFPE tissues was extracted using five 20 µm slides from the tissue blocks using the QIAamp DNA FFPE Tissue Kit (Qiagen, Hilden, Germany) following the manufacturer’s recommendation with the only exception that the Proteinase K digestion step was conducted overnight until the tissue was completely dissolved.
The quality of the DNA extractions was evaluated using the Infinium FFPE QC Kit (Illumina, San Diego, CA, USA). Extracted DNA with an amplification cycle difference > 5 compared to the Infinium FFPE QC template DNA was excluded. The bisulfite conversion was performed to unveil methylated DNA bases using EZ DNA Methylation™ Kit (Zymo Research, Irvine, CA, USA) with 250–500 ng DNA. The elution of DNA was performed using 10 µL elution buffer. The eluate was subsequently restored using the Infinium HD FFPE DNA Restore Kit (Illumina, San Diego, CA, USA). The DNA methylation levels were measured using the Infinium MethylationEPIC BeadChip (Illumina, San Diego, CA, USA) following the manufacturer’s protocol (# 15,021,525 Rev. B). The slides were scanned using the iScan system (Illumina, San Diego, CA, USA).
All data analysis was conducted in R (version 3.6.0) using tidyverse (version 1.3.0) (Wickham et al., 2019), patchwork (version 1.1.1), and RColorBrewer (version 1.1-2).
The raw .idat files were imported into R using the minfi package (version 1.32.0) [67,68]. The samples with a mean detection p-value > 0.05 were excluded. The R package, quantro (v.1.20) [69], was used to perform a data-driven test to ensure that the assumption for conducting Quantile Normalisation of the data (no global differences in distributions between the SUD and SUDEP cases) was not violated. The signals were normalised using the preprocessQuantile function [68,70] and subsequently converted into methylation values. β- and M-values were calculated using the getBeta and getM functions, respectively. The three functions are all implemented in minfi. The PCA was conducted based on M-values and visualised using the ggbiplot package (version 0.55).

4.5. Differentially Methylated Regions

Regions of consecutive CpG positions with different methylation between SUD and SUDEP cases, i.e., DMRs, were identified using DMRcate (version 2.2.3) [38]. A CpG M-value matrix was annotated with chromosomal positions and test statistics using a design matrix with gender, age, and sample slide information as co-variables. The DMRs were identified using the recommended scaling factor = 2 and a minimum of three consecutive CpG positions in the DMR. The analysis of the DMRs was conducted using the plyranges (version 1.6.10) [38], GenomicRanges (version 1.38.0) [71], org.Hs.eg.db (version 3.10), annotate (version 1.64), and TxDb.Hsapiens.UCSC.hg19.knownGene (version 3.2.2) packages. Only regions with Δβ > 0.10 were considered as DMRs.
The GSTT1 DMR was visualised using the Gviz package (version 1.30.3) [72] from R (version 3.6.0). The Hg19 RefSeq gene track was constructed using the UcscTrack function within the Gviz package and org.Hs.eg.db. Information of the DNase sensitive regions (wgEncodeRegDnaseClusteredV3) were downloaded from the University of California Santa Cruz (UCSC) Genome Browser. The regions were originally identified in 125 cell types and tissues as a part of the ENCODE project [73,74]. The Spearman’s correlation coefficients between the degrees of DNA methylation of the CpG positions and the gene expression in SUDEP were calculated using the cor.test function in R.

Author Contributions

Conceptualisation: S.N.C., S.B.J., J.D.A., M.-L.K., J.B. and N.M.; Methodology: S.N.C., S.B.J., J.D.A., M.-L.K. and L.C.T.; Resources: K.B.O. and J.B.; Investigation: S.N.C., S.B.J., M.-L.K. and L.C.T.; Formal analysis: S.N.C., S.B.J. and J.D.A.; Supervision: J.T.-H., J.B. and N.M.; Writing—Original Draft Preparation: S.N.C. and S.B.J., Writing—Review and Editing: S.N.C., S.B.J., J.D.A., M.-L.K., J.T.-H. and N.M., All authors have read and agreed to the published version of the manuscript.

Funding

This work was kindly supported by Ellen and Aage Andersen’s foundation and the Novo Nordic foundation No. NFF18OC0031634.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki, and approved by the Committees on Health Research Ethics in the Capital Region of Denmark (H-2-2012-017) and the Danish Data Protection Agency (2011-54-1262).

Informed Consent Statement

Due to autopsies being performed prior to 1 January 2012, informed consent was waived in accordance with the Danish Health Act and the Consolidated Act no. 1338 of 1 September 2020 on Research Ethics Review of Health Research Projects and Health Data Research Projects.

Data Availability Statement

Data cannot be shared publicly because of data protection legislations (Danish Data Protection Agency no. 2011-54-1262). Data are available from the University of Copenhagen, Department of Forensic Medicine Institutional Data Access/Ethics Committee (contact via [email protected], University of Copenhagen, Department of Forensic Medicine) for researchers who meet the criteria for access to confidential data.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Skjelbred, T.; Lynge, T.H.; Nielsen, J.; Winkel, B.G.; Tfelt-Hansen, J. Symptoms and healthcare contact preceding sudden cardiac death in persons aged 1–49 years. Trends Cardiovasc. Med. 2020, 31, 119–124. [Google Scholar] [CrossRef] [PubMed]
  2. Mellor, G.; Raju, H.; De Noronha, S.V.; Papadakis, M.; Sharma, S.; Behr, E.R.; Sheppard, M.N. Clinical characteristics and circumstances of death in the sudden arrhythmic death syndrome. Circ. Arrhythmia Electrophysiol. 2014, 7, 1078–1083. [Google Scholar] [CrossRef] [Green Version]
  3. Winkel, B.G.; Holst, A.G.; Theilade, J.; Kristensen, I.B.; Thomsen, J.L.; Ottesen, G.L.; Bundgaard, H.; Svendsen, J.H.; Haunsø, S.; Tfelt-Hansen, J. Nationwide study of sudden cardiac death in persons aged 1–35 years. Eur. Heart J. 2011, 32, 983–990. [Google Scholar] [CrossRef] [Green Version]
  4. Bagnall, R.D.; Weintraub, R.G.; Ingles, J.; Duflou, J.; Yeates, L.; Lam, L.; Davis, A.M.; Thompson, T.; Connell, V.; Wallace, J.; et al. A Prospective Study of Sudden Cardiac Death among Children and Young Adults. N. Engl. J. Med. 2016, 374, 2441–2452. [Google Scholar] [CrossRef]
  5. Risgaard, B.; Winkel, B.G.; Jabbari, R.; Behr, E.R.; Ingemann-Hansen, O.; Thomsen, J.L.; Ottesen, G.L.; Gislason, G.H.; Bundgaard, H.; Haunsø, S.; et al. Burden of Sudden Cardiac Death in Persons Aged 1 to 49 Years. Circ. Arrhythmia Electrophysiol. 2014, 7, 205–211. [Google Scholar] [CrossRef] [Green Version]
  6. Semsarian, C.; Ingles, J.; Wilde, A.A.M. Sudden cardiac death in the young: The molecular autopsy and a practical approach to surviving relatives. Eur. Heart J. 2015, 36, 1290–1296. [Google Scholar] [CrossRef]
  7. Basso, C.; Carturan, E.; Pilichou, K.; Rizzo, S.; Corrado, D.; Thiene, G. Sudden cardiac death with normal heart: Molecular autopsy. Cardiovasc. Pathol. 2010, 19, 321–325. [Google Scholar] [CrossRef] [PubMed]
  8. Schwartz, P.J.; Ackerman, M.J.; Antzelevitch, C.; Bezzina, C.R.; Borggrefe, M.; Cuneo, B.F.; Wilde, A.A.M. Inherited cardiac arrhythmias. Nat. Rev. Dis. Prim. 2020, 6, 58. [Google Scholar] [CrossRef] [PubMed]
  9. Lahrouchi, N.; Raju, H.; Lodder, E.M.; Papatheodorou, E.; Ware, J.S.; Papadakis, M.; Tadros, R.; Cole, D.; Skinner, J.R.; Crawford, J.; et al. Utility of Post-Mortem Genetic Testing in Cases of Sudden Arrhythmic Death Syndrome. J. Am. Coll. Cardiol. 2017, 69, 2134–2145. [Google Scholar] [CrossRef]
  10. Christiansen, S.L.; Hertz, C.L.; Ferrero-Miliani, L.; Dahl, M.; Weeke, P.E.; LuCamp; Ottesen, G.L.; Frank-Hansen, R.; Bundgaard, H.; Morling, N. Genetic investigation of 100 heart genes in sudden unexplained death victims in a forensic setting. Eur. J. Hum. Genet. 2016, 24, 1797–1802. [Google Scholar] [CrossRef]
  11. Neubauer, J.; Lecca, M.R.; Russo, G.; Bartsch, C.; Medeiros-Domingo, A.; Berger, W.; Haas, C. Exome analysis in 34 sudden unexplained death (SUD) victims mainly identified variants in channelopathy-associated genes. Int. J. Legal Med. 2018, 132, 1057–1065. [Google Scholar] [CrossRef]
  12. Meder, B.; Haas, J.; Sedaghat-Hamedani, F.; Kayvanpour, E.; Frese, K.; Lai, A.; Nietsch, R.; Scheiner, C.; Mester, S.; Bordalo, D.M.; et al. Epigenome-Wide Association Study Identifies Cardiac Gene Patterning and a Novel Class of Biomarkers for Heart Failure. Circulation 2017, 136, 1528–1544. [Google Scholar] [CrossRef] [PubMed]
  13. Grunert, M.; Dorn, C.; Cui, H.; Dunkel, I.; Schulz, K.; Schoenhals, S.; Sun, W.; Berger, F.; Chen, W.; Sperling, S.R. Comparative DNA methylation and gene expression analysis identifies novel genes for structural congenital heart diseases. Cardiovasc. Res. 2016, 112, 464–477. [Google Scholar] [CrossRef]
  14. Haas, J.; Frese, K.S.; Park, Y.J.; Keller, A.; Vogel, B.; Lindroth, A.M.; Weichenhan, D.; Franke, J.; Fischer, S.; Bauer, A.; et al. Alterations in cardiac DNA methylation in human dilated cardiomyopathy. EMBO Mol. Med. 2013, 5, 413–429. [Google Scholar] [CrossRef]
  15. Movassagh, M.; Choy, M.K.; Knowles, D.A.; Cordeddu, L.; Haider, S.; Down, T.; Siggens, L.; Vujic, A.; Simeoni, I.; Penkett, C.; et al. Distinct epigenomic features in end-stage failing human hearts. Circulation 2011, 124, 2411–2422. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Holst, A.G.; Winkel, B.G.; Risgaard, B.; Nielsen, J.B.; Rasmussen, P.V.; Haunsø, S.; Sabers, A.; Uldall, P.; Tfelt-Hansen, J. Epilepsy and risk of death and sudden unexpected death in the young: A nationwide study. Epilepsia 2013, 54, 1613–1620. [Google Scholar] [CrossRef]
  17. Risgaard, B.; Waagstein, K.; Winkel, B.G.; Jabbari, R.; Lynge, T.H.; Glinge, C.; Albert, C.; Correll, C.U.; Hauns, S.; Fink-Jensen, A.; et al. Sudden cardiac death in young adults with previous hospital-based psychiatric inpatient and outpatient treatment: A nationwide cohort study from Denmark. J. Clin. Psychiatry 2015, 76, 1122–1129. [Google Scholar] [CrossRef] [PubMed]
  18. Lynge, T.H.; Svane, J.; Pedersen-Bjergaard, U.; Gislason, G.; Torp-Pedersen, C.; Banner, J.; Risgaard, B.; Winkel, B.G.; Tfelt-Hansen, J. Sudden cardiac death among persons with diabetes aged 1-49 years: A 10-year nationwide study of 14 294 deaths in Denmark. Eur. Heart J. 2020, 41, 2699–2706. [Google Scholar] [CrossRef]
  19. Thurman, D.J.; Logroscino, G.; Beghi, E.; Hauser, W.A.; Hesdorffer, D.C.; Newton, C.R.; Scorza, F.A.; Sander, J.W.; Tomson, T. The burden of premature mortality of epilepsy in high-income countries: A systematic review from the Mortality Task Force of the International League against Epilepsy. Epilepsia 2017, 58, 17–26. [Google Scholar] [CrossRef]
  20. Shorvon, S.; Tomson, T. Sudden unexpected death in epilepsy. Lancet 2011, 378, 2028–2038. [Google Scholar] [CrossRef]
  21. Ali, A.; Wu, S.; Issa, N.P.; Rose, S.; Towle, V.L.; Warnke, P.; Tao, J.X. Association of sleep with sudden unexpected death in epilepsy. Epilepsy Behav. 2017, 76, 1–6. [Google Scholar] [CrossRef]
  22. Chahal, C.A.A.; Salloum, M.N.; Alahdab, F.; Gottwald, J.A.; Tester, D.J.; Anwer, L.A.; So, E.L.; Murad, M.H.; St Louis, E.K.; Ackerman, M.J.; et al. Systematic Review of the Genetics of Sudden Unexpected Death in Epilepsy: Potential Overlap With Sudden Cardiac Death and Arrhythmia-Related Genes. J. Am. Heart Assoc. 2020, 9, e012264. [Google Scholar] [CrossRef] [PubMed]
  23. Bagnall, R.D.; Singer, E.S.; Tfelt-Hansen, J. Sudden Cardiac Death in the Young. Hear. Lung Circ. 2020, 29, 498–504. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. DeGiorgio, C.M.; Markovic, D.; Mazumder, R.; Moseley, B.D. Ranking the leading risk factors for sudden unexpected death in epilepsy. Front. Neurol. 2017, 8, 473. [Google Scholar] [CrossRef]
  25. Devinsky, O.; Hesdorffer, D.C.; Thurman, D.J.; Lhatoo, S.; Richerson, G. Sudden unexpected death in epilepsy: Epidemiology, mechanisms, and prevention. Lancet Neurol. 2016, 15, 1075–1088. [Google Scholar] [CrossRef]
  26. Manolis, T.A.; Manolis, A.A.; Melita, H.; Manolis, A.S. Sudden unexpected death in epilepsy: The neuro-cardio-respiratory connection. Seizure 2019, 64, 65–73. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Johnson, J.N.; Hofman, N.; Haglund, C.M.; Cascino, G.D.; Wilde, A.A.M.; Ackerman, M.J. Identification of a possible pathogenic link between congenital long QT syndrome and epilepsy. Neurology 2009, 72, 224–231. [Google Scholar] [CrossRef] [Green Version]
  28. Anderson, J.H.; Bos, J.M.; Cascino, G.D.; Ackerman, M.J. Prevalence and spectrum of electroencephalogram-identified epileptiform activity among patients with long QT syndrome. Hear. Rhythm 2014, 11, 53–57. [Google Scholar] [CrossRef] [PubMed]
  29. Goldman, A.M.; Glasscock, E.; Yoo, J.; Chen, T.T.; Klassen, T.L.; Noebels, J.L. Arrhythmia in heart and brain: KCNQ1 mutations link epilepsy and sudden unexplained. Sci. Transl. Med. 2009, 2, 2ra6. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Bagnall, R.D.; Crompton, D.E.; Petrovski, S.; Lam, L.; Cutmore, C.; Garry, S.I.; Sadleir, L.G.; Dibbens, L.M.; Cairns, A.; Kivity, S.; et al. Exome-based analysis of cardiac arrhythmia, respiratory control, and epilepsy genes in sudden unexpected death in epilepsy. Ann. Neurol. 2016, 79, 522–534. [Google Scholar] [CrossRef]
  31. Bagnall, R.D.; Ingles, J.; Yeates, L.; Berkovic, S.F.; Semsarian, C. Exome sequencing–based molecular autopsy of formalin-fixed paraffin-embedded tissue after sudden death. Genet. Med. 2017, 19, 1127–1133. [Google Scholar] [CrossRef] [Green Version]
  32. Tu, E.; Bagnall, R.D.; Duflou, J.; Semsarian, C. Post-mortem review and genetic analysis of sudden unexpected death in epilepsy (SUDEP) cases. Brain Pathol. 2011, 21, 201–208. [Google Scholar] [CrossRef] [PubMed]
  33. Coll, M.; Allegue, C.; Partemi, S.; Mates, J.; Del Olmo, B.; Campuzano, O.; Pascali, V.; Iglesias, A.; Striano, P.; Oliva, A.; et al. Genetic investigation of sudden unexpected death in epilepsy cohort by panel target resequencing. Int. J. Legal Med. 2016, 130, 331–339. [Google Scholar] [CrossRef]
  34. Aurlien, D.; Leren, T.P.; Taubøll, E.; Gjerstad, L. New SCN5A mutation in a SUDEP victim with idiopathic epilepsy. Seizure 2009, 18, 158–160. [Google Scholar] [CrossRef] [Green Version]
  35. Jones, W.; Greytak, S.; Odeh, H.; Guan, P.; Powers, J.; Bavarva, J.; Moore, H.M. Deleterious effects of formalin-fixation and delays to fixation on RNA and miRNA-Seq profiles. Sci. Rep. 2019, 9, 6980. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Gilbert, M.T.P.; Sanchez, J.J.; Haselkorn, T.; Jewell, L.D.; Lucas, S.B.; Van Marck, E.; Børsting, C.; Morling, N.; Worobey, M. Multiplex PCR with minisequencing as an effective high-throughput SNP typing method for formalin-fixed tissue. Electrophoresis 2007, 28, 2361–2367. [Google Scholar] [CrossRef] [PubMed]
  37. Do, H.; Dobrovic, A. Sequence artifacts in DNA from formalin-fixed tissues: Causes and strategies for minimization. Clin. Chem. 2015, 61, 64–71. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Peters, T.J.; Buckley, M.J.; Statham, A.L.; Pidsley, R.; Samaras, K.; V Lord, R.; Clark, S.J.; Molloy, P.L. De novo identification of differentially methylated regions in the human genome. Epigenetics Chromatin 2015, 8, 6. [Google Scholar] [CrossRef] [Green Version]
  39. Bagnall, R.D.; Semsarian, C. Role of the molecular autopsy in the investigation of sudden cardiac death. Prog. Pediatr. Cardiol. 2017, 45, 17–23. [Google Scholar] [CrossRef]
  40. Allocati, N.; Masulli, M.; Di Ilio, C.; Federici, L. Glutathione transferases: Substrates, inihibitors and pro-drugs in cancer and neurodegenerative diseases. Oncogenesis 2018, 7, 8. [Google Scholar] [CrossRef]
  41. Haas, J.; Mester, S.; Lai, A.; Frese, K.S.; Sedaghat-Hamedani, F.; Kayvanpour, E.; Rausch, T.; Nietsch, R.; Boeckel, J.; Carstensen, A.; et al. Genomic structural variations lead to dysregulation of important coding and non-coding RNA species in dilated cardiomyopathy. EMBO Mol. Med. 2018, 10, 107–120. [Google Scholar] [CrossRef] [PubMed]
  42. Shanks, G.W.; Tester, D.J.; Ackerman, J.P.; Simpson, M.A.; Behr, E.R.; White, S.M.; Ackerman, M.J. Importance of Variant Interpretation in Whole-Exome Molecular Autopsy. Circulation 2018, 137, 2705–2715. [Google Scholar] [CrossRef]
  43. Dulhunty, A.F.; Hewawasam, R.; Liu, D.; Casarotto, M.G.; Board, P.G. Regulation of the cardiac muscle ryanodine receptor by glutathione transferases. Drug Metab. Rev. 2011, 43, 236–252. [Google Scholar] [CrossRef]
  44. Glasscock, E. Genomic biomarkers of SUDEP in brain and heart. Epilepsy Behav. 2014, 38, 172–179. [Google Scholar] [CrossRef] [Green Version]
  45. Lehnart, S.E.; Mongillo, M.; Bellinger, A.; Lindegger, N.; Chen, B.-X.; Hsueh, W.; Reiken, S.; Wronska, A.; Drew, L.J.; Ward, C.W.; et al. Leaky Ca2+ release channel/ryanodine receptor 2 causes seizures and sudden cardiac death in mice. J. Clin. Invest. 2008, 118, 2230–2245. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Prabha, T.S.; Kumaraswami, K.; Kutala, V.K. Association of GSTT1 and GSTM1 polymorphisms in South Indian epilepsy patients. Indian J. Exp. Biol. 2016, 54, 783–787. [Google Scholar] [PubMed]
  47. Ercegovac, M.; Jovic, N.; Sokic, D.; Savic-Radojevic, A.; Coric, V.; Radic, T.; Nikolic, D.; Kecmanovic, M.; Matic, M.; Simic, T.; et al. GSTA1, GSTM1, GSTP1 and GSTT1 polymorphisms in progressive myoclonus epilepsy: A Serbian case-control study. Seizure 2015, 32, 30–36. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Pearson-Smith, J.N.; Patel, M. Metabolic dysfunction and oxidative stress in epilepsy. Int. J. Mol. Sci. 2017, 18, 2365. [Google Scholar] [CrossRef] [Green Version]
  49. Coll, M.; Oliva, A.; Grassi, S.; Brugada, R.; Campuzano, O. Update on the genetic basis of sudden unexpected death in epilepsy. Int. J. Mol. Sci. 2019, 20, 1979. [Google Scholar] [CrossRef] [Green Version]
  50. Nashef, L.; Walker, F.; Allen, P.; Sander, J.W.A.S.; Shorvon, S.D.; Fish, D.R. Apnoea and bradycardia during epileptic seizures: Relation to sudden death in epilepsy. J. Neurol. Neurosurg. Psychiatry 1996, 60, 297–300. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  51. Ryvlin, P.; Nashef, L.; Lhatoo, S.D.; Bateman, L.M.; Bird, J.; Bleasel, A.; Boon, P.; Crespel, A.; Dworetzky, B.A.; Høgenhaven, H.; et al. Incidence and mechanisms of cardiorespiratory arrests in epilepsy monitoring units (MORTEMUS): A retrospective study. Lancet Neurol. 2013, 12, 966–977. [Google Scholar] [CrossRef]
  52. Stewart, M.; Silverman, J.B.; Sundaram, K.; Kollmar, R. Causes and Effects Contributing to Sudden Death in Epilepsy and the Rationale for Prevention and Intervention. Front. Neurol. 2020, 11, 765. [Google Scholar] [CrossRef] [PubMed]
  53. Khera, A.V.; Mason-Suares, H.; Brockman, D.; Wang, M.; VanDenburgh, M.J.; Senol-Cosar, O.; Patterson, C.; Newton-Cheh, C.; Zekavat, S.M.; Pester, J.; et al. Rare Genetic Variants Associated With Sudden Cardiac Death in Adults. J. Am. Coll. Cardiol. 2019, 74, 2623–2634. [Google Scholar] [CrossRef]
  54. Robinson, M.D.; Kahraman, A.; Law, C.W.; Lindsay, H.; Nowicka, M.; Weber, L.M.; Zhou, X. Statistical methods for detecting differentially methylated loci and regions. Front. Genet. 2014, 5, 234. [Google Scholar] [CrossRef] [Green Version]
  55. Teschendorff, A.E.; Relton, C.L. Statistical and integrative system-level analysis of DNA methylation data. Nat. Rev. Genet. 2018, 19, 129–147. [Google Scholar] [CrossRef]
  56. Kling, T.; Wenger, A.; Beck, S.; Carén, H. Validation of the MethylationEPIC BeadChip for fresh-frozen and formalin-fixed paraffin-embedded tumours. Clin. Epigenetics 2017, 9, 33. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Andersen, G.B.; Hager, H.; Hansen, L.L.; Tost, J. Improved reproducibility in genome-wide DNA methylation analysis for PAXgene-fixed samples compared with restored formalin-fixed and paraffin-embedded DNA. Anal. Biochem. 2015, 468, 50–58. [Google Scholar] [CrossRef] [Green Version]
  58. Trudsø, L.C.; Andersen, J.D.; Jacobsen, S.B.; Christiansen, S.L.; Congost-Teixidor, C.; Kampmann, M.; Morling, N. A comparative study of single nucleotide variant detection performance using three massively parallel sequencing methods. PLoS ONE 2020, 15, e0239850. [Google Scholar] [CrossRef]
  59. R Core Team. R: A language and environment for statistical computing. Available online: https://www.r-project.org/ (accessed on 5 March 2021).
  60. Andrews, S. FastQC—A quality control tool for high throughput sequence data. Available online: https://www.bioinformatics.babraham.ac.uk/projects/fastqc/ (accessed on 5 March 2021).
  61. Lindgreen, S. AdapterRemoval: Easy cleaning of next-generation sequencing reads. BMC Res. Notes 2012, 17, 76. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Dobin, A.; Davis, C.A.; Schlesinger, F.; Drenkow, J.; Zaleski, C.; Jha, S.; Batut, P.; Chaisson, M.; Gingeras, T.R. STAR: Ultrafast universal RNA-seq aligner. Bioinformatics 2013, 29, 15–21. [Google Scholar] [CrossRef]
  63. Harrow, J.; Frankish, A.; Gonzalez, J.M.; Tapanari, E.; Diekhans, M.; Kokocinski, F.; Aken, B.L.; Barrell, D.; Zadissa, A.; Searle, S.; et al. GENCODE: The reference human genome annotation for the ENCODE project. Genome Res. 2012, 22, 1760–1774. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Andersen, J.D.; Jacobsen, S.B.; Trudsø, L.C.; Kampmann, M.-L.; Banner, J.; Morling, N. Whole genome and transcriptome sequencing of post-mortem cardiac tissues from sudden cardiac death victims identifies a gene regulatory variant in NEXN. Int. J. Legal Med. 2019, 133, 1699–1709. [Google Scholar] [CrossRef] [PubMed]
  65. Robinson, M.D.; Oshlack, A. A scaling normalization method for differential expression analysis of RNA-seq data. Genome Biol. 2010, 11, R25. [Google Scholar] [CrossRef] [Green Version]
  66. Robinson, M.D.; McCarthy, D.J.; Smyth, G.K. edgeR: A Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 2009, 26, 139–140. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  67. Fortin, J.P.; Triche, T.J.; Hansen, K.D. Preprocessing, normalization and integration of the Illumina HumanMethylationEPIC array with minfi. Bioinformatics 2017, 33, 558–560. [Google Scholar] [CrossRef]
  68. Aryee, M.J.; Jaffe, A.E.; Corrada-Bravo, H.; Ladd-Acosta, C.; Feinberg, A.P.; Hansen, K.D.; Irizarry, R.A. Minfi: A flexible and comprehensive Bioconductor package for the analysis of Infinium DNA methylation microarrays. Bioinformatics 2014, 30, 1363–1369. [Google Scholar] [CrossRef] [Green Version]
  69. Hicks, S.C.; Irizarry, R.A. quantro: A data-driven approach to guide the choice of an appropriate normalization method. Genome Biol. 2015, 16, 117. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  70. Touleimat, N.; Tost, J. Complete pipeline for Infinium® Human Methylation 450K BeadChip data processing using subset quantile normalization for accurate DNA methylation estimation. Epigenomics 2012, 4, 325–341. [Google Scholar] [CrossRef] [Green Version]
  71. Lee, S.; Cook, D.; Lawrence, M. Plyranges: A grammar of genomic data transformation. Genome Biol. 2019, 20, 4. [Google Scholar] [CrossRef] [PubMed]
  72. Lawrence, M.; Huber, W.; Pagès, H.; Aboyoun, P.; Carlson, M.; Gentleman, R.; Morgan, M.T.; Carey, V.J. Software for Computing and Annotating Genomic Ranges. PLoS Comput. Biol. 2013, 9, e1003118. [Google Scholar] [CrossRef] [PubMed]
  73. Hahne, F.; Ivanek, R. Visualizing genomic data using Gviz and bioconductor. In Methods in Molecular Biology; Humana Press Inc.: Totowa, NJ, USA, 2016; Volume 1418, pp. 335–351. ISBN 978-1-4939-3578-9. [Google Scholar]
  74. ENCODE Project Consortium. An integrated encyclopedia of DNA elements in the human genome. Nature 2012, 489, 57–74. [Google Scholar] [CrossRef] [PubMed]
Ijms 22 02790 g001 550
Figure 1. Principal component analysis plot. The principal component analysis was based on M-values. The samples are coloured according to age and the shape of the symbols shows diagnosis. PC: Principal component. SUD: sudden unexplained death. SUDEP: sudden unexpected death in epilepsy.
Figure 1. Principal component analysis plot. The principal component analysis was based on M-values. The samples are coloured according to age and the shape of the symbols shows diagnosis. PC: Principal component. SUD: sudden unexplained death. SUDEP: sudden unexpected death in epilepsy.
Ijms 22 02790 g001
Ijms 22 02790 g002 550
Figure 2. Regional differences of β-values from sudden unexplained death and sudden unexpected death in epilepsy cases overlapping the GSTT1 region. (A) DNA methylation is shown for CpG positions in the differentially methylated region (DMR) overlapping the GSTT1 region in ventricular tissues from sudden unexplained death (SUD) and sudden unexpected death in epilepsy (SUDEP). ENCODE DNase I sensitive regions are shown in blue with colour intensities corresponding to the signal strength. The lines connecting the points in the β-value track indicate the mean for each group. The DMR is located at chr22:24,384,105-24,384,944 covering 13 CpG positions. The colors in the chromosome structure correspond to chromosome staining bands and the vertical red line indicates the chromosomal location of the DMR. (B) The median of the three Trimmed Mean of M-values (TMM) per kb values from each individual was plotted against the β-value of each of the positions that were statistical significantly correlated (p ≤ 0.05) with GSTT1 expression. Each CpG position is labelled with the CpG loci IDs according to the CpG nomenclature of Illumina.
Figure 2. Regional differences of β-values from sudden unexplained death and sudden unexpected death in epilepsy cases overlapping the GSTT1 region. (A) DNA methylation is shown for CpG positions in the differentially methylated region (DMR) overlapping the GSTT1 region in ventricular tissues from sudden unexplained death (SUD) and sudden unexpected death in epilepsy (SUDEP). ENCODE DNase I sensitive regions are shown in blue with colour intensities corresponding to the signal strength. The lines connecting the points in the β-value track indicate the mean for each group. The DMR is located at chr22:24,384,105-24,384,944 covering 13 CpG positions. The colors in the chromosome structure correspond to chromosome staining bands and the vertical red line indicates the chromosomal location of the DMR. (B) The median of the three Trimmed Mean of M-values (TMM) per kb values from each individual was plotted against the β-value of each of the positions that were statistical significantly correlated (p ≤ 0.05) with GSTT1 expression. Each CpG position is labelled with the CpG loci IDs according to the CpG nomenclature of Illumina.
Ijms 22 02790 g002
Table
Table 1. Descriptive data of the study population.
Table 1. Descriptive data of the study population.
CharacteristicsSUD (n = 9)SUDEP (n = 14)p-Value
Males, n (%)6 (67)10 (71)1.000
Median age at time of death, years (range)35 (0–50)35 (3–49)0.850
Witnessed death, n (%)2 (22)1 (7)0.538
Performed toxicology screen, n (%)6 (67)11 (79)0.643
No compounds detected, n (%)2 (33)2 (18)0.584
One compound detected, n (%)2 (33)5 (46)1.000
More than one compound detected, n (%)2 (33)4 (36)1.000
Antiepileptic drug(s) detected, n (%)0 (0)8 (73)0.009
Activity at time of death
Sleep, n (%)2 (22)5 (36)0.657
At rest, n (%)3 (33)3 (21)0.643
Passenger in car, n (%)1 (11)0 (0)0.391
Unknown, n (%)3 (33)6 (43)1.000
SUD: Sudden unexplained death. SUDEP = Sudden unexpected death in epilepsy.
Table
Table 2. Chromosomal regions with different methylation levels in sudden unexplained death compared to sudden unexpected death in epilepsy sorted by minimum false discovery rate.
Table 2. Chromosomal regions with different methylation levels in sudden unexplained death compared to sudden unexpected death in epilepsy sorted by minimum false discovery rate.
ChrStartEndSize CpG Min. FDRΔβNearest GeneDistance
chr2224,384,10524,384,944840133.81×10–230.15GSTT10
chr1221,053,841221,054,825985106.04×10–06−0.12HLX0
chr1110,230,252110,230,63338271.12×10–05−0.14GSTM10
chr2224,373,32224,373,61829731.17×10–040.19LOC3913220
chr115,959,6585,959,94528834.61×10–040.11TRIM50
chr2224,348,54924,348,71516736.15×10–04−0.16GSTT41290
Chr: Chromosome. Size: Size of the differentially methylated region (DMR) in bases. CpG: Number of CpG sites in the region. Min. FDR: Minimum false discovery rate (FDR) of the kernel-smoothed estimate of the DMR. Δβ: Mean β-value coefficient across the DMR between sudden unexplained death (SUD) and sudden unexpected death in epilepsy (SUDEP). Nearest gene: Nearest or overlapping gene for the DMR. Distance: Distance to nearest gene in bases. The DMR overlapped a gene region for values equal to zero.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Christiansen, S.N.; Jacobsen, S.B.; Andersen, J.D.; Kampmann, M.-L.; Trudsø, L.C.; Olsen, K.B.; Tfelt-Hansen, J.; Banner, J.; Morling, N. Differential Methylation in the GSTT1 Regulatory Region in Sudden Unexplained Death and Sudden Unexpected Death in Epilepsy. Int. J. Mol. Sci. 2021, 22, 2790. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms22062790

AMA Style

Christiansen SN, Jacobsen SB, Andersen JD, Kampmann M-L, Trudsø LC, Olsen KB, Tfelt-Hansen J, Banner J, Morling N. Differential Methylation in the GSTT1 Regulatory Region in Sudden Unexplained Death and Sudden Unexpected Death in Epilepsy. International Journal of Molecular Sciences. 2021; 22(6):2790. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms22062790

Chicago/Turabian Style

Christiansen, Steffan Noe, Stine Bøttcher Jacobsen, Jeppe Dyrberg Andersen, Marie-Louise Kampmann, Linea Christine Trudsø, Kristine Boisen Olsen, Jacob Tfelt-Hansen, Jytte Banner, and Niels Morling. 2021. "Differential Methylation in the GSTT1 Regulatory Region in Sudden Unexplained Death and Sudden Unexpected Death in Epilepsy" International Journal of Molecular Sciences 22, no. 6: 2790. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms22062790

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop