Abstract
Atrial fibrillation (AF) is one of the most prevalent heart rhythm disorder. The causes of AF include age, male sex, diabetes, hypertension, valve disease, and systolic/diastolic dysfunction. But on molecular level, its mechanisms are largely unknown. In this study, we collected 10 patients with persistent atrial fibrillation, 10 patients with paroxymal atrial fibrillation and 10 healthy individuals and did Methylation EPICBead Chip and RNA sequencing. By analyzing the methylation and gene expression data using machine learning based feature selection method Boruta, we identified the key genes that were strongly associated with AF and found their interconnections. The results suggested that the methylation of KIF15 may regulate the expression of PSMC3, TINAG, and NUDT6. The identified AF associated methylation-expression regulations may help understand the molecular mechanisms of AF from a multi-omics perspective.
Introduction
Atrial fibrillation (AF), one of the most prevalent heart rhythm disorders, is a potential cause of heart failure and ischemic stroke with high morbidity and mortality (; ). The cause of AF is multifactorial which include age, male sex, diabetes, hypertension, valve disease, and systolic/diastolic dysfunction (; ; ). Depends on how often atrial fibrillation occurs and how it responds to treatment, AF is roughly divided into two major subtypes-paroxysmal atrial fibrillation (PAF) and persistent atrial fibrillation (PeAF). In the treatment of AF, drugs were the first choice, non-drug therapies were used only when drug therapy failed or patients could not tolerate the medication. In contrast to the extensive knowledge of etiology, the underlying mechanism of AF remains elusive. Further study of the potential mechanisms of AF could provide novel strategies for the treatment and management to increase quality of life and reduce economic burden of social ().
With the development of next generation sequencing (NGS) technologies, growing evidence have demonstrated that AF is a disease with a significant genetic contribution. Previous studies have filtered novel genetic variants and candidate genes including transcriptional factor genes (PITX2, PRRX1, ZFHX3, NKX2.5, TBX5), ion channel genes (KCNN3, HCN4, CACNA1C, SCN5A, KCNQ1, KCNH2), and caveolin genes (CAV1 and CAV2) (; ; ; ; ). However, these genes explain only a small fraction of the biology and genetic underpinnings of AF.
Epidemiological studies have demonstrated that genetic, environmental, behavioral, and clinical factors contribute to AF pathogenic mechanism (). Emelia J. B performed genome-wide methylation using whole blood samples from 183 prevalent AF and 220 incident AF cases. They examined the association between DNA methylation and GWAS loci, suggesting DNA methylation might be a possible mechanism through which AF-specific genetic variations affect gene regulation (). To date, only a few studies have investigated differential DNA methylation as a predictor biomarker at specific candidate loci that were previously associated with AF.
Therefore, we applied DNA methylation profiling study to identify the likely rare damaging variants and putative candidate genes from 10 patients with persistent atrial fibrillation (PeAF), 10 patients with paroxymal atrial fibrillation (PAF) and 10 healthy individuals. Interestingly, we identified top 10 genes (KIF15, ABCA3, FOXG1, VGF, PDE4D, EIF3C, CNTNAP5, SHOX2, VGF, TRIM59) as functional candidate genes and the expression level are significantly increased in PeAF and PAF patients than control. Given the importance of DNA methylation to gene expression, we investigated the gene expression of the same participants using RNA sequencing. We also defined top 10 genes (EPN3, EMD, SMCO4, F2RL2, TMED1, PSMC3, PDZD11, NUDT6, TINAG, GALNT5) in gene expression data and the expression pattern of these genes was significantly different between PeAF and PAF. These results have improved our understanding of the underlying mechanism and offer new insights into the potential pathway of AF, which could provide novel therapeutic option for this disease.
Materials and Methods
Atrial Fibrillation Patients
Ten patients with paroxymal atrial fibrillation (g1), 10 patients with persistent atrial fibrillation (continuous atrial fibrillation lasting more than 12 months) (g2) and 10 healthy individuals (g3) were enrolled in this study (Table 1). All patients were subjected to detailed medical evaluation, which included medical history, physical examination, electrocardiography (ECG), and echocardiography. Patients with chronic heart failure, coronary heart disease, cardiomyopathy, hyperthyroidism or chronic pulmonary heart disease were excluded.
TABLE 1
| Coronary | |||||||||||
| Age | Weight | Height | Diabetes | angiography | LVEF | Left atrial | Duration of | ||||
| No. | (years) | Gender | (Kg) | (cm) | Smoking | Hypertension | mellitus | or CTA | (%) | diameter (mm) | AF (years) |
| 1 | 69 | Male | 76 | 169 | No | Yes | No | Negative | 70 | 40 | – |
| 2 | 63 | Male | 64 | 170 | No | No | No | Negative | 59 | 46 | – |
| 3 | 63 | Male | 70 | 170 | No | No | No | Negative | 66 | 39 | – |
| 4 | 69 | Male | 67 | 173 | No | Yes | No | Negative | 67 | 46 | – |
| 5 | 69 | Male | 75 | 165 | No | No | No | Negative | 70 | 36 | – |
| 6 | 61 | Male | 76 | 176 | No | Yes | Yes | Negative | 60 | 42 | – |
| 7 | 64 | Male | 52 | 168 | Yes | No | No | Negative | 64 | 40 | – |
| 8 | 64 | Male | 71 | 181 | No | Yes | Yes | Negative | 63 | 39 | – |
| 9 | 61 | Male | 87 | 167 | Yes | Yes | No | Negative | 62 | 37 | – |
| 10 | 66 | Male | 82 | 173 | No | Yes | No | Negative | 63 | 42 | – |
| 11 | 63 | Male | 86 | 176 | No | Yes | No | Negative | 57 | 46 | 2.5 |
| 12 | 63 | Male | 80 | 178 | No | No | No | Negative | 68 | 55 | 3 |
| 13 | 64 | Male | 70 | 170 | No | No | No | Negative | 67 | 41 | 4 |
| 14 | 64 | Male | 84 | 164 | No | Yes | No | Negative | 55 | 48 | 2 |
| 15 | 65 | Male | 73 | 169 | No | Yes | No | Negative | 69 | 55 | 3.5 |
| 16 | 66 | Male | 66 | 168 | No | Yes | No | Negative | 64 | 45 | 4 |
| 17 | 67 | Male | 80 | 175 | No | Yes | Yes | Negative | 59 | 47 | 2.5 |
| 18 | 67 | Male | 73 | 165 | Yes | Yes | No | Negative | 59 | 47 | 3 |
| 19 | 63 | Male | 61 | 164 | No | No | No | Negative | 73 | 49 | 2 |
| 20 | 67 | Male | 90 | 178 | No | Yes | No | Negative | 70 | 58 | 2.5 |
Demographic characteristics of AF patients.
LVEF, left ventricular ejection fraction; AF, atrial fibrillation; CTA, CT angiography.
The study was conducted in accordance with the Declaration of Helsinki, and the protocol used to collect human heart tissue was approved by the Ethics Committee of Shanghai East Hospital (DI: 0402017).
Written informed consents to participate in this study were provided by all the enrolled patients before operation of fibrillation ablation. The left atrial appendage tissues which were abandoned during isolated surgical ablation were collected. Normal left atrial appendages were collected from healthy male donors.
The Methylation Profiles
The DNA methylation status of 850K probes in the 30 samples was measured using Methylation EPICBead Chip. The raw data was quality controlled and preprocessed using R/Bioconductor package minfi1 (). The beta value ranged from 0 to 1 was calculated to represent how each position was methylated. 1 meant high methylation and 0 meant low methylation.
The RNA Sequencing Profiles
The total RNAs were extracted using RNeasy Mini Kit (Cat#74106, Qiagen) and the RNA integrity was checked using Agilent Bioanalyzer 2100 (Agilent technologies, Santa Clara, CA, United States). Qualified total RNA was further purified by RNAClean XP Kit (Cat A63987, Beckman Coulter Inc., Kraemer Boulevard, Brea, CA, United States) and RNase-Free DNase Set (Cat#79254, QIAGEN, GmBH, Germany). Pair-end sequencing reads were generated using Illumina data collection software. First, the reads were mapped onto human reference genome GRCh38 using Hisat2 (version:2.0.42) (). Then, Stringtie (version:1.3.03) () was used to calculate the FPKM (Fragments Per Kilobase of exon model per Million mapped reads).
Feature Selection Algorithm
There were 866,091 methylation probes and 50,868 RNAs. The number of features were extremely large. It was difficult to select key features using traditional statistical methods. Therefore, we adopted the latest machine learning based feature selection method Boruta to get the key methylation probes and RNAs.
Boruta () is a feature selection method based on random forest. It can select sample group relevant features effectively. First, it will shuffle the features to create many permuted datasets. Then, it will evaluate the importance score of each feature in the original actual dataset and permuted datasets. Then, it will compare the actual importance score with permuted scores and find the features with significantly higher actual importance scores than permuted scores. After multiple iterations, it will select all the sample group relevant features. The python code from https://github.com/scikit-learn-contrib/boruta_py was used to apply the Boruta feature selection algorithm.
Classification Predictor
To evaluate how well the selected features can classify the samples, we built an SVM (Support Vector Machine) classifier using the methylation data and another RNA-Seq data based SVM classifier. The svm function in R package e101714 () was used to apply the SVM classification algorithm.
LOOCV (leave-one-out cross validation) was used to objectively evaluate the classification performance. Each time, one sample was treated as test sample while all the other samples were used to train the model. After 30 rounds, all samples had been tested once and the overall accuracy was calculated based on the confusion matrix. In confusion matrix, the actual sample groups were compared with predicted sample groups.
Results
The Key Methylation Features Identified With Boruta
We ran Boruta feature selection algorithm on the methylation data and got 10 key methylation features. These 10 key methylation features were listed in Table 2. The probes were annotated to genome positions (Genome Build 37) and genes using the official annotation file from Illumina. Sometime, one position may be associated with multiple genes. Therefore, the 10 methylation probes can be mapped onto 15 gene symbols. We checked the GO annotation of these genes and found that cg16703882 (SHOX2) was associated with GO:0007507: heart development which was closely relevant to AF.
TABLE 2
| ILMNID | Chromosome | Position | Strand | UCSC Ref gene name |
| cg00702638 | 3 | 44803293 | R | KIF15; KIAA1143 |
| cg02331561 | 16 | 2391081 | F | ABCA17P; ABCA3 |
| cg02991338 | 14 | 29236017 | R | FOXG1 |
| cg04084157 | 7 | 100809049 | F | VGF |
| cg05995159 | 5 | 59325256 | R | PDE4D |
| cg06357615 | 16 | 28403195 | R | MIR6862-2; MIR6862-1; EIF3CL; EIF3C |
| cg11344566 | 2 | 124782885 | F | CNTNAP5 |
| cg16703882 | 3 | 157823479 | R | SHOX2 |
| cg21186299 | 7 | 100808810 | R | VGF |
| cg26856080 | 3 | 160167746 | R | TRIM59 |
The 10 key methylation features identified by Boruta.
We plotted the heatmap of these 10 key methylation features in Figure 1. It can be seen that most of the 10 patients with paroxymal atrial fibrillation (g1), 10 patients with persistent atrial fibrillation (g2) and 10 healthy individuals (g3) were cluster into the right groups.
FIGURE 1
The Key Gene Expression Features Identified With Boruta
Similarly, we ran Boruta feature selection algorithm on the RNA-Seq gene expression data and got 10 key gene expression features. The 10 key gene expression features were given in Table 3. We also plotted their heatmap in Figure 2. The clusters were also largely correct.
TABLE 3
| Gene ID | Name | Description | GRCh38 locus |
| ENSG00000049283 | EPN3 | Epsin 3 | 17:50532543-50543750 |
| ENSG00000102119 | EMD | Emerin | X:154379197-154381523 |
| ENSG00000166002 | SMCO4 | Single-pass membrane protein with coiled-coil domains 4 | 11:93478472-93543508 |
| ENSG00000164220 | F2RL2 | Coagulation factor II (thrombin) receptor-like 2 | 5:76615482-76623434 |
| ENSG00000099203 | TMED1 | Transmembrane p24 trafficking protein 1 | 19:10832438-10836318 |
| ENSG00000165916 | PSMC3 | Proteasome 26S subunit, atpase 3 | 11:47418769-47426473 |
| ENSG00000120509 | PDZD11 | PDZ domain containing 11 | X:70286595-70290514 |
| ENSG00000170917 | NUDT6 | Nudix hydrolase 6 | 4:122888697-122922968 |
| ENSG00000137251 | TINAG | Tubulointerstitial nephritis antigen | 6:54307859-54390152 |
| ENSG00000136542 | GALNT5 | Polypeptide N-acetylgalactosaminyltransferase 5 | 2:157257598-157314211 |
The 10 key gene expression features identified by Boruta.
FIGURE 2
The Classification Performance of Methylation and Gene Expression
From Figures 1, 2, we can see that both methylation and gene expression features can correctly cluster most samples. But we would like to evaluate their performance objectively and quantitatively. Therefore, we applied LOOCV to test the SVM classifiers of methylation and gene expression features. The confusion matrixes of methylation features and gene expression features were listed in Tables 4, 5, respectively.
TABLE 4
| Predicted g1 | Predicted g2 | Predicted g3 | |
| Actual g1 | 8 | 2 | 0 |
| Actual g2 | 3 | 7 | 0 |
| Actual g3 | 0 | 0 | 10 |
The confusion matrix of key methylation features.
TABLE 5
| Predicted g1 | Predicted g2 | Predicted g3 | |
| Actual g1 | 9 | 0 | 1 |
| Actual g2 | 0 | 9 | 1 |
| Actual g3 | 1 | 0 | 9 |
The confusion matrix of key gene expression features.
From Table 3, we can see that the AF patients (g1 + g2) and healthy individuals (g3) were perfectly classified using the methylation features, but the methylation data did not have great performance on classifying the subtype of AF (g1 and g2). From Table 4, we can see that the two subtype of AF, paroxysmal atrial fibrillation (g1) and persistent atrial fibrillation (g2), had very different gene expression pattern. In other words, the methylation data and gene expression data complement each other. The methylation data can be used to predict the AF and the gene expression data can be used to classify the subtypes of AF.
We checked the wrongly predicted samples in Tables 3, 4. They were different. Within the 30 samples, 22 samples had the same predicted labels by expression and methylation. All these 22 samples were correctly predicted. For the 8 inconsistent samples between expression and methylation predictions, at least one of the two predictions (the expression-based prediction and the methylation-based prediction) was correct. In other words, all the samples can be corrected classified based on either expression or methylation. The expression-based prediction and the methylation-based prediction were complementary.
The Methylation-Expression Regulation Network
We mapped the genes of methylation features and gene expression features to the STRING network (Version 11.05) () and visualized the network using R package igraph ()6 to identify the potential relationship between two candidate genes sets. The methylation-expression regulation network was shown in Figure 3. In the network, the methylation and expression genes were marked in red and green. The methylation genes located in three clusters: EIF3CL-EIF3C, KIF15-TRIM59, and SHOX2-FOXG1-CNTNAP5. The three expression genes (PSMC3, TINAG, and NUDT6) were connected with methylation gene KIF15. Even EPN3 can be indirectly connected to KIF15. That made KIF15 at the center of the network. These results suggested that KIF5 may play important roles in the pathogenesis of AF.
FIGURE 3
Discussion
DNA methylation, a pre-transcriptional modification characterized by the addition of methyl groups to specific nucleotides, regulates the stability of gene expression states and maintains genome integrity by collaborating with proteins that modify nucleosomes (; ; ). Previous studies considered that changes in DNA methylation states contribute to the regulation of biological processes underlying AF, such as fibrosis, atrial dilatation, atrial fibroblast proliferation and differentiation from fibroblasts into myofibroblasts (). To further enhance the biological understanding of the atrial fibrillation, our study focused on DNA methylation, particularly with respect to how it relates to mRNA expression. Among our two gene sets of top 10 genes, we found PDED4, SHOX2, and EMD were the most important genes for AF which have been reported associated with AF in previous reference.
Atrial fibrillation is reported to be associated with a profound remodeling of membrane receptors and alterations in cAMP dependent regulation of Ca2+ handling. PDE4 is expressed in human atrial myocytes and accounts for approximately 15% of PDE (phosphodiesterase) activity (). PDE4D encoded protein has 3’,5’-cyclic-AMP phosphodiesterase activity and degrades cAMP, which acts as a signal transduction molecule in multiple cell types and represents the major PDE4 subtype (). The activity of PDEs decreased with age, and the relative PDED4 activity was lower in patients with permanent atrial fibrillation than in age-matched sinus rhythm controls (). Previous study provided evidence that patients with pAF were found to have a decreased PDE4 activity as compared with patients in sinus rhythm ().
Short Stature Homeobox 2 (SHOX2) is a member of the homeobox family of genes in which mutations associated with early-onset and familial AF (). SHOX2 is considered as a key regulator of sinus node development of which deficiency could lead to bradycardia in animal models (). Previous study demonstrated SHOX2 was susceptible for SND and AF by screening 98 SND patients and 450 individuals with AF. In the heart development of mouse and zebrafish, they also proved SHOX2 plays an important role, the mutation of SHOX2 could lead to severe bradycardia (; ).
Emerin (EMD) encodes a serine-rich nuclear membrane protein which located on the cytoplasmic surface of the inner nuclear membrane and related to X-linked Emery-Dreifuss muscular dystrophy (EDMD) (). Previous study found a nonsense mutation in EMD from two EDMD families which is associated with X-linked recessive inheritance, result in serious cardiac complication, including AF (). Cardiologic assessment revealed slow atrial fibrillation in a recent case of a 65-year-old male patient with a hemizygous duplication of 5 bases in exon 6 of the EMD, gene on the X chromosome (; ; ).
Statements
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation, to any qualified researcher.
Ethics statement
The studies involving human participants were reviewed and approved by the Ethics Committee of Shanghai East Hospital. The patients/participants provided their written informed consent to participate in this study. Written informed consent was obtained from the individual(s) for the publication of any potentially identifiable images or data included in this manuscript.
Author contributions
YZ and CG conceived and designed the experiments. BL, XS, and KD performed the experiments. ML, YQ, and SZ analyzed the data. BL, XS, and YZ wrote the manuscript. All authors read and approved the final version of the manuscript.
Funding
The present study was funded by the National Key Research and Development Program (Grant No. 2018-YFC-1312505 to YZ), the National Natural Science Foundation of China (Grant No. 81770408 to CG), and (Grant No. 81900280 to XS).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Footnotes
1.^https://bioconductor.org/packages/minfi/
2.^http://ccb.jhu.edu/software/hisat2/
3.^http://ccb.jhu.edu/software/stringtie/
References
1
AryeeM. J.JaffeA. E.Corrada-BravoH.Ladd-AcostaC.FeinbergA. P.HansenK. D.et al (2014). Minfi: a flexible and comprehensive bioconductor package for the analysis of infinium DNA methylation microarrays.Bioinformatics301363–1369. 10.1093/bioinformatics/btu049
2
AsmaratsL.Cruz-GonzalezI.Nombela-FrancoL.ArzamendiD.PeralV.NietlispachF.et al (2019). Recurrence of device-related thrombus after percutaneous left atrial appendage closure.Circulation1401441–1443. 10.1161/CIRCULATIONAHA.119.040860
3
BerkE.ChristT.SchwarzS.RavensU.KnautM.KaumannA. J. (2016). In permanent atrial fibrillation, PDE3 reduces force responses to 5-HT, but PDE3 and PDE4 do not cause the blunting of atrial arrhythmias.Br. J. Pharmacol.1732478–2489. 10.1111/bph.13525
4
BlaschkeR. J.HahurijN. D.KuijperS.JustS.WisseL. J.DeisslerK.et al (2007). Targeted mutation reveals essential functions of the homeodomain transcription factor Shox2 in sinoatrial and pacemaking development.Circulation1151830–1838. 10.1161/CIRCULATIONAHA.106.637819
5
BrissetM.Ben YaouR.CarlierR. Y.ChanutA.NicolasG.RomeroN. B.et al (2019). X-linked Emery-Dreifuss muscular dystrophy manifesting with adult onset axial weakness, camptocormia, and minimal joint contractures.Neuromuscul. Disord.29678–683. 10.1016/j.nmd.2019.06.009
6
CapanniC.Del CocoR.MattioliE.CamozziD.ColumbaroM.SchenaE.et al (2009). Emerin-prelamin a interplay in human fibroblasts.Biol. Cell101541–554. 10.1042/BC20080175
7
ChangC.-C.LinC.-J. (2011). LIBSVM: a library for support vector machines.ACM Trans. on Intell. Syst. Technol.21–27. 10.1145/1961189.1961199
8
ChughS. S.BlackshearJ. L.ShenW. K.HammillS. C.GershB. J. (2001). Epidemiology and natural history of atrial fibrillation: clinical implications.J. Am. Coll. Cardiol.37371–378. 10.1016/s0735-1097(00)01107-4
9
CsardiG.NepuszT. (2006). The igraph software package for complex network research.Inter J. Complex Syst.16951–9.
10
EllinorP. T.LunettaK. L.GlazerN. L.PfeuferA.AlonsoA.ChungM. K.et al (2010). Common variants in KCNN3 are associated with lone atrial fibrillation.Nat. Genet.42240–244. 10.1038/ng.537
11
HoffmannS.PaoneC.SumerS. A.DieboldS.WeissB.RoethR.et al (2019). Functional characterization of rare variants in the shox2 gene identified in sinus node dysfunction and atrial fibrillation.Front. Genet.10:648. 10.3389/fgene.2019.00648
12
KimD.LangmeadB.SalzbergS. L. (2015). HISAT: a fast spliced aligner with low memory requirements.Nat. Methods12357–360. 10.1038/nmeth.3317
13
KisselJ. T.DimbergE. L.Emslie-SmithA. M.SelcenD.KeeganB. M. (2009). A 49-year-old man with contractures, weakness, and cardiac arrhythmia.Neurology722036–2043. 10.1212/01.wnl.0000350522.42635.09
14
KursaM.RudnickiW. (2010). Feature selection with the boruta package.J. Statist. Softw. Artic.361–13. 10.18637/jss.v036.i11
15
LinH.YinX.XieZ.LunettaK. L.LubitzS. A.LarsonM. G.et al (2017). Methylome-wide association study of atrial fibrillation in framingham heart study.Sci. Rep.7:40377. 10.1038/srep40377
16
LipG. Y. H.LaneD. A.BullerH.ApostolakisS. (2013). Development of a novel composite stroke and bleeding risk score in patients with atrial fibrillation: the AMADEUS study.Chest1441839–1847. 10.1378/chest.13-1635
17
MaB.WilkerE. H.Willis-OwenS. A.ByunH. M.WongK. C.MottaV.et al (2014). Predicting DNA methylation level across human tissues.Nucleic Acids Res.423515–3528. 10.1093/nar/gkt1380
18
MaJ. F.YangF.MahidaS. N.ZhaoL.ChenX.ZhangM. L.et al (2016). TBX5 mutations contribute to early-onset atrial fibrillation in chinese and caucasians.Cardiovasc. Res.109442–450. 10.1093/cvr/cvw003
19
MiltonA. G.AykanatV. M.Hamilton-BruceM. A.NezicM.JannesJ.KoblarS. A. (2011). Association of the phosphodiesterase 4D (PDE4D) gene and cardioembolic stroke in an Australian cohort.Int. J. Stroke6480–486. 10.1111/j.1747-4949.2011.00616.x
20
MolinaC. E.LeroyJ.RichterW.XieM.ScheitrumC.LeeI. O.et al (2012). Cyclic adenosine monophosphate phosphodiesterase type 4 protects against atrial arrhythmias.J. Am. Coll. Cardiol.592182–2190. 10.1016/j.jacc.2012.01.060
21
NielsenJ. B.ThorolfsdottirR. B.FritscheL. G.ZhouW.SkovM. W.GrahamS. E.et al (2018). Biobank-driven genomic discovery yields new insight into atrial fibrillation biology.Nat. Genet.501234–1239. 10.1038/s41588-018-0171-3
22
OgawaH.SenooK.AnY.ShantsilaA.ShantsilaE.LaneD. A.et al (2017). Clinical features and prognosis in patients with atrial fibrillation and prior stroke: comparing the fushimi and darlington AF registries.eBio Med.18199–203. 10.1016/j.ebiom.2017.03.022
23
OlesenM. S.NielsenM. W.HaunsoS.SvendsenJ. H. (2014). Atrial fibrillation: the role of common and rare genetic variants.Eur. J. Hum. Genet.22297–306. 10.1038/ejhg.2013.139
24
PerteaM.PerteaG. M.AntonescuC. M.ChangT.-C.MendellJ. T.SalzbergS. L. (2015). StringTie enables improved reconstruction of a transcriptome from RNA-seq reads.Nat. Biotechnol.33290–295. 10.1038/nbt.3122
25
SakataK.ShimizuM.InoH.YamaguchiM.TeraiH.FujinoN.et al (2005). High incidence of sudden cardiac death with conduction disturbances and atrial cardiomyopathy caused by a nonsense mutation in the STA gene.Circulation1113352–3358. 10.1161/CIRCULATIONAHA.104.527184
26
SchnabelR. B.SullivanL. M.LevyD.PencinaM. J.MassaroJ. M.D’AgostinoR. B.et al (2009). Development of a risk score for atrial fibrillation (Framingham Heart Study): a community-based cohort study.Lancet373739–745. 10.1016/S0140-6736(09)60443-8
27
ShenK.TuT.YuanZ.YiJ.ZhouY.LiaoX.et al (2017). DNA methylation dysregulations in valvular atrial fibrillation.Clin. Cardiol.40686–691. 10.1002/clc.22715
28
SinnerM. F.TuckerN. R.LunettaK. L.OzakiK.SmithJ. G.TrompetS.et al (2014). Integrating genetic, transcriptional, and functional analyses to identify 5 novel genes for atrial fibrillation.Circulation1301225–1235. 10.1161/CIRCULATIONAHA.114.009892
29
SzklarczykD.GableA. L.LyonD.JungeA.WyderS.Huerta-CepasJ.et al (2018). STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets.Nucleic Acids Res.47D607–D613. 10.1093/nar/gky1131
30
TaoH.ShiK. H.YangJ. J.LiJ. (2016). Epigenetic mechanisms in atrial fibrillation: new insights and future directions.Trends Cardiovasc. Med.26306–318. 10.1016/j.tcm.2015.08.006
31
Vicente-SteijnR.KelderT. P.TertoolenL. G.WisseL. J.PijnappelsD. A.PoelmannR. E.et al (2017). RHOA-ROCK signalling is necessary for lateralization and differentiation of the developing sinoatrial node.Cardiovasc. Res.1131186–1197. 10.1093/cvr/cvx104
32
VoukalisC.LipG. Y.ShantsilaE. (2016). Emerging Tools for Stroke Prevention in Atrial Fibrillation.eBio Med.426–39. 10.1016/j.ebiom.2016.01.017
33
YeW.WangJ.SongY.YuD.SunC.LiuC.et al (2015). A common Shox2-Nkx2-5 antagonistic mechanism primes the pacemaker cell fate in the pulmonary vein myocardium and sinoatrial node.Development1422521–2532. 10.1242/dev.120220
34
YehY. H.EhrlichJ. R.QiX.HebertT. E.ChartierD.NattelS. (2007). Adrenergic control of a constitutively active acetylcholine-regulated potassium current in canine atrial cardiomyocytes.Cardiovasc. Res.74406–415. 10.1016/j.cardiores.2007.01.020
35
ZhaoG.ZhouJ.GaoJ.LiuY.GuS.ZhangX.et al (2017). Genome-wide DNA methylation analysis in permanent atrial fibrillation.Mol. Med. Rep.165505–5514. 10.3892/mmr.2017.7221
36
ZhaoJ.LiuT.LiG. (2014). Relationship between two arrhythmias: sinus node dysfunction and atrial fibrillation.Arch. Med. Res.45351–355. 10.1016/j.arcmed.2014.04.005
37
ZhongJ.AghaG.BaccarelliA. A. (2016). The role of DNA Methylation in cardiovascular risk and disease: methodological aspects, study design, and data analysis for epidemiological studies.Circ. Res.118119–131. 10.1161/CIRCRESAHA.115.305206
Summary
Keywords
atrial fibrillation, methylation, multi-omics, feature selection, classification
Citation
Liu B, Shi X, Ding K, Lv M, Qian Y, Zhu S, Guo C and Zhang Y (2020) The Joint Analysis of Multi-Omics Data Revealed the Methylation-Expression Regulations in Atrial Fibrillation. Front. Bioeng. Biotechnol. 8:187. doi: 10.3389/fbioe.2020.00187
Received
19 December 2019
Accepted
26 February 2020
Published
12 March 2020
Volume
8 - 2020
Edited by
Peilin Jia, The University of Texas Health Science Center at Houston (UTHealth), United States
Reviewed by
Yan Yao, Chinese Academy of Medical Sciences, China; Jian Zhong Zhou, The First Affiliated Hospital of Chongqing Medical University, China
Updates
Copyright
© 2020 Liu, Shi, Ding, Lv, Qian, Zhu, Guo and Zhang.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Changfa Guo, guo.changfa@zs-hospital.sh.cnYangyang Zhang, zhangyangyang_wy@vip.sina.com; zhangyangyang_wy@tongji.edu.cn
†These authors have contributed equally to this work
This article was submitted to Bioinformatics and Computational Biology, a section of the journal Frontiers in Bioengineering and Biotechnology
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.