Next Article in Journal
Uncovering the Metabolic and Stress Responses of Human Embryonic Stem Cells to FTH1 Gene Silencing
Next Article in Special Issue
Electrophysiological Remodeling: Cardiac T-Tubules and ß-Adrenoceptors
Previous Article in Journal
Sustainable Agriculture through Multidisciplinary Seed Nanopriming: Prospects of Opportunities and Challenges
Previous Article in Special Issue
Remodeling of Cardiac Gap Junctional Cell–Cell Coupling
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Artificial Intelligence-Assisted Identification of Genetic Factors Predisposing High-Risk Individuals to Asymptomatic Heart Failure

1
Division of Cardiology, Department of Internal Medicine, Chang Gung Memorial Hospital, Keelung 204, Taiwan
2
Community Medicine Research Center, Chang Gung Memorial Hospital, Keelung 204, Taiwan
3
College of Medicine, Chang Gung University, Taoyuan 333, Taiwan
4
Department of Thoracic and Cardiovascular Surgery, Chang Gung Memorial Hospital, New Taipei City 333, Taiwan
5
Advanced Tech BU, Acer Inc., New Taipei City 221, Taiwan
6
Institute of Molecular and Genomic Medicine, National Health Research Institute, Zhunan 350, Taiwan
7
Department of Nursing, Chang Gung University of Science and Technology, Taoyuan 333, Taiwan
8
Department of Ophthalmology, Chang Gung Memorial Hospital, Keelung 204, Taiwan
9
National Institute of Infectious Diseases and Vaccinology, National Health Research Institutes, Zhunan 350, Taiwan
10
Department & Graduate Institute of Microbiology and Immunology, National Defense Medical Center, Taipei 114, Taiwan
11
Department of Life Sciences and Institute of Genome Sciences, National Yang Ming Chiao Tung University, Taipei 112, Taiwan
12
Aging and Health Research Center, National Yang Ming Chiao Tung University, Taipei 112, Taiwan
*
Authors to whom correspondence should be addressed.
Submission received: 31 July 2021 / Revised: 7 September 2021 / Accepted: 13 September 2021 / Published: 15 September 2021
(This article belongs to the Special Issue Electrical Remodeling in Cardiac Disease)

Abstract

:
Heart failure (HF) is a global pandemic public health burden affecting one in five of the general population in their lifetime. For high-risk individuals, early detection and prediction of HF progression reduces hospitalizations, reduces mortality, improves the individual’s quality of life, and reduces associated medical costs. In using an artificial intelligence (AI)-assisted genome-wide association study of a single nucleotide polymorphism (SNP) database from 117 asymptomatic high-risk individuals, we identified a SNP signature composed of 13 SNPs. These were annotated and mapped into six protein-coding genes (GAD2, APP, RASGEF1C, MACROD2, DMD, and DOCK1), a pseudogene (PGAM1P5), and various non-coding RNA genes (LINC01968, LINC00687, LOC105372209, LOC101928047, LOC105372208, and LOC105371356). The SNP signature was found to have a good performance when predicting HF progression, namely with an accuracy rate of 0.857 and an area under the curve of 0.912. Intriguingly, analysis of the protein connectivity map revealed that DMD, RASGEF1C, MACROD2, DOCK1, and PGAM1P5 appear to form a protein interaction network in the heart. This suggests that, together, they may contribute to the pathogenesis of HF. Our findings demonstrate that a combination of AI-assisted identifications of SNP signatures and clinical parameters are able to effectively identify asymptomatic high-risk subjects that are predisposed to HF.

1. Introduction

Heart failure (HF) is an important public health problem that is associated with high morbidity, high mortality, and a burden to healthcare [1]. It is considered to be a progressive disorder and is caused by a range of different risk factors; this leads to a heterogenic pathophysiology. The symptoms of HF include effort intolerance and/or fluid retention, as well as dyspnea, fatigue, and pulmonary congestion. The American College of Cardiology and American Heart Association (ACC/AHA) have categorized HF into four stages: Stage A is defined as having risk factors for HF only; Stage B is defined as having structural heart disease without any current or prior symptoms of HF; Stage C is defined as symptomatic HF; and Stage D is HF refractory to treatment [2].
Stage B HF, which is asymptomatic, is a risk factor for mortality. Stage B HF is defined by the ACC/AHA as best characterized by an increase in the left ventricular (LV) mass (LVM), an increase in the left atrial dimensions, the presence of LV geometric patterns indicative of adverse remodeling (i.e., concentric remodeling and/or both concentric and eccentric hypertrophy), and a reduction in the LV ejection fraction (LVEF) or diastolic dysfunction. However, Stage B HF subjects are asymptomatic, thus it is highly unlikely this clinically silent population will routinely receive a detailed examination of their heart. This leads to a five-fold increase in mortality risk among such individuals, as well as the transition from Stage B HF to Stage C HF, which is associated with elevated rates of hospitalization and death [3]. The ACC/AHA guidelines have emphasized the importance of the appropriate treatment of Stage B HF subjects to prevent the development of symptomatic HF [2] and have highlighted the need for early detection strategies to identify these clinically asymptomatic Stage B HF subjects.
The risk prediction of a complex disease such as the asymptomatic Stage B HF is currently a challenge as well as is an unmet clinical need. It is well documented that the risk factors for HF include environmental factors, metabolic derangements, and genetic factors [2], and there is also an increasing appreciation that there is an underlying strong heritable component. This strengthens the importance of discovering genetic factors that contribute to the underlying mechanism in an attempt to reveal novel targets for the prevention and treatment of HF. Notably, recent studies of single nucleotide polymorphisms (SNPs) have suggested that there is an association between certain genetic factors and an increased risk of HF [4,5]. However, the role of those factors and the molecular mechanism(s) behind the pathogenesis of HF remain incompletely understood.
Artificial intelligence (AI) techniques have been applied previously to cardiovascular diseases, namely model prediction of the presence of HF, estimation of the HF subtype, and assessment of the severity of HF [6]. Currently, most studies of AI-assisted HF prediction have used clinical features and focused on the recognition of subtypes for prognosis [7,8], such as destabilizations, re-hospitalizations, and mortality [9,10]. In addition, the techniques of AI machine-learning also have great potential for delineating complex biological processes, in particular those involving interactions between the multiple genetic factors and biochemical pathways that accelerate the development of HF. Accordingly, in this study, we applied an AI-assisted methodology to identify the genetic factors in a high-risk population that are potentially associated with asymptomatic Stage B HF. This involved carrying out genome-wide SNPs screening. Furthermore, we also performed protein connectivity mapping of the genes in which the SNPs are located in order to pinpoint their potential role in the molecular pathogenesis of HF in terms of functional connectivity and protein–protein interaction networks.

2. Materials and Methods

2.1. Study Subjects

Between February 2019 and November 2019, 162 prospectively recruited participants from the Northeastern Taiwan Community Medicine Research Cohort (NTCMRC, ClinicalTrials.gov Identifier: NCT04839796) were enrolled and examined in the cardiology outpatient clinic of Chang Gung Memorial Hospital, Keelung. All subjects received a clinical examination, blood tests, electrocardiogram (ECG) and echocardiography evaluation, as well as having a complete personal and past medical history takenrecording. The 10-year and lifetime risk of atherosclerotic cardiovascular disease (ASCVD) was calculated for each subject [11]. The inclusion criteria were as follows: subjects were over 30 years old with a 10-year ASCVD risk ≥20%. The exclusion criteria were that the subject was already known to have clinical HF at Stage C or D; was suffering from atrial fibrillation (as identified by having a previous diagnosis of atrial fibrillation or paroxysmal atrial fibrillation; found to have documented atrial fibrillation on ECG or during the echocardiography exam); and/or was pregnant. Informed signed consent was obtained from all participants. This study protocol conforms to the ethical guidelines of the 1975 Declaration of Helsinki and was approved by the Institutional Review Board of Chang Gung Medical Foundation (IRB No: 201800802B0 and 202000077B0A3).

2.2. Clinical Assessment

At recruitment, all participants provided a detailed personal history and received a full physical examination, completed a questionnaire, and underwent various biochemical tests, including assays for n-terminal pro-brain natriuretic peptide (NT-proBNP), high sensitivity-Troponin T (hs-Tnt), and high sensitivity c-reactive protein (hs-CRP), as well as a two-dimensional (2D) echocardiography examination. Blood pressure was measured using the average of two seated measurements. Heart rate was measured via a resting 12-lead ECG. Body mass index was calculated as weight divided by height2 and expressed as kg/m2. Diabetes mellitus was defined as a fasting glucose of ≥126 mg/dL, a random glucose of ≥200 mg/dL, or the use of hypoglycemic medications. Previous history of coronary heart disease was used to identify if subjects had angina pectoris with a positive exercise test result; a history of myocardial infarction; angiographic evidence of significant (>75%) coronary artery stenosis after intra-coronary nitroglycerine 50–200 μg administration; a history of percutaneous coronary revascularization; or coronary artery bypass grafting. Current smoking status was defined as having smoked more than 100 cigarettes in their lifetime and having smoked within 1 month before enrollment.

2.3. Biochemical Analysis

Blood specimens were collected in citrate-treated tubes at recruitment. After centrifuging for at least 15 min, the plasma component was frozen and shipped on dry ice to the core laboratory center of our hospital, at which the samples were stored at −80 °C for the subsequent measurement of cytokines and inflammatory markers. Plasma hs-CRP was measured in duplicates by an enzyme-linked immunosorbent assay on the basis of purified protein and polyclonal anti-C-reactive protein antibodies (IMMULITE hs-CRP, Diagnostic Products Corporation, Los Angeles, CA, USA). The lower limit of this assay was 0.10 mg/L and the coefficient of variation was ≤5% at the 0.20 mg/L C-reactive protein level. The plasma concentrations of hs-Tnt and NT-proBNP were measured using appropriate sandwich enzyme-linked immunosorbent assay kits and the monoclonal antibody targeting the relevant cytokine (R&D Systems, Inc., Minneapolis, MN, USA). Other related biomarkers, namely leptin and adiponectin, were measured using commercially available enzyme-linked immunosorbent assays (Boster Biological Technology, Pleasanton, CA, USA).

2.4. Echocardiography

Echocardiography was performed within one month of recruitment. A comprehensive transthoracic Doppler echocardiography was performed using a commercially available machine (Vivid E9 system, General Electrics, Boston, MA, USA) with a M5S probe. LV end-diastolic and end-systolic volume were measured from the apical two-chamber and four-chamber view, and the LVEF was calculated using the modified biplane Simpson’s rule. The LVM index was measured according to the American Society of Echocardiography formula. Conventional Doppler parameters were measured according to a standardized examination procedure and these were early (E) and late diastolic transmitral flow velocity; deceleration time of E; average of the septal annular mitral early diastolic, late diastolic, and systolic tissue velocities (E’); and the ratio of E/E’. The pulmonary artery systolic pressure was calculated using the modified Bernoulli equation from the tricuspid regurgitation peak jet velocity and estimated right atrial pressure (from the respiratory variation of the inferior vena cava diameter). Two-dimensional strain analysis was performed using custom 2D strain-imaging software (EchoPac, GE Ultrasound, Boston, MA, USA). The endocardial borders were traced from the end-systolic frame of the 2D images. Interactive software then automatically tracked myocardial motion and divided each image into six segments. Numerical and graphical displays of the deformation parameters (reflecting the average value for the tracking of all acoustic markers in each segment) were then generated for all six segments from each view. Longitudinal peak systolic strain was acquired for the apical two-chamber, three chamber, and four-chamber views. If further shortening occurred after the end of the systole, this was measured as the peak strain. Global longitudinal strain was calculated as the average longitudinal strain of the segments of two-chamber, three-chamber, and four-chamber views. All examinations of echocardiography were performed and analyzed by one experienced physician (Dr. Ning-I Yang) of excellent reproducibility [12] who was blinded to the subject data.

2.5. Definition of Stages A and Stage B Heart Failure

HF preclinical stages were assessed based on the clinical history and echocardiographic data. Stage A HF was defined as the presence of risk factors such as arterial hypertension; type 2 diabetes; obesity; metabolic syndrome; a documented clinical history of atherosclerotic disease or the use of cardiotoxins when there was no evidence of structural heart disease; or the signs/symptoms of HF. Stage B HF was defined as the presence of structural heart disease or the detection of diastolic dysfunction on the echocardiographic examination. The latter subject needed to fulfill one of the following criteria: (1) LV hypertrophy as defined by an LV mass index of >95 g/m2 in women or of >115 g/m2 in men; (2) LV dilatation as defined by an LV end diastolic volume index of >95 mL/m2; (3) concentric remodeling defined by a relative wall thickness of >0.42; (4) asymptomatic LV dysfunction, including an LVEF of < 50% and/or diastolic dysfunction (≥grade II) without clinical signs and/or symptoms of HF; and (5) more than mild mitral or aortic regurgitation.

2.6. DNA Extraction of White Blood Cells

Peripheral venous blood was collected from the subjects and processed on the same day. Each blood sample was centrifuged at 3000 rpm for 10 min at 4 °C to separate the serum from cells. Genomic DNA was then isolated from peripheral white blood cells using the phenol/chloroform DNA extraction method after lysis of red blood cells. Finally, precipitation and washing using 95% isopropanol followed by 80% alcohol were used to obtain the total genomic DNA.

2.7. Whole-Genome SNP Analysis

To identify single nucleotide polymorphisms (SNPs), we genotyped the genomic DNA using AxiomTM Genome-Wide TWB 2.0 array plates, which included 686,463 SNPs. Genotyping analyses were performed on 117 high-risk (ASCVD risk ≥ 20%) subjects, including 83 Stage A and 34 Stage B asymptomatic HF subjects. The stages of HF were defined by the echocardiographic examination results. Among the SNPs identified, those with a minor allele frequency rate of 0 or those SNPs with a missing rate of more than 10% were excluded from further analysis. In total, 392,885 SNPs were available for further analyses.

2.8. AI-Assisted Discovery of Candidate SNPs

For the model training and testing, all machine-learning analyses were performed using R Version 3.5.3 (using the random Forest, e1071, glmnet, rpart, caret, and cvAUC packages). We used three supervised algorithms to select important features, namely the random forest (RF), support vector machine (SVM), and least absolute shrinkage and selection operator (LASSO) methods, with the input dataset having a train-to-validation split ratio of 80:20. The SNPs were ranked by the summation of the selected counts using 100-time bootstrapped random samples and the three machine-learning methods. The minimum features needed were the highest performance for the area under the curve (AUC) and the accuracy rate; these were calculated for the four different machine-learning models (RF, SVM, LASSO, and the decision tree).

2.9. Protein Interaction Network

The locations of AI-identified SNP signatures were mapped to their relevant genes. To explore the regulatory mechanisms and potential pathways that the genes may be involved in, the six protein-coding genes from the twelve signature genes were subjected to protein–protein interaction analyses using the BioGRID database [13]. The results are displayed as a graphical network using the open-source software Cytoscape [14].

2.10. Statistics

Two independent sample t-tests were used to compare differences between the continuous variables derived from the groups. The results are presented as means ± standard deviations (SD). The χ2 test was used to examine the distribution of the categorical variables and results are expressed as frequencies and percentages between the groups. A multiple logistic regression model was used to determine the strength of association between the selected parameters and the presence of Stage B HF. The statistical software used for this study was SPSS 25.0 (IBM Corporation, Armonk, NY, USA).

3. Results

3.1. Clinical, Biochemical, and Echocardiography Data

A total of 162 subjects were recruited into the study between February 2019 and November 2019. The demographic and clinical baseline characteristics of the subjects are shown in Table 1. Of these, 113 subjects (70%) had Stage A HF and 49 subjects (30%) had Stage B (asymptomatic) HF. Both groups had a male predominance and the ASCVD risk scores in both groups were in the high category (>20%). The risk factors for HF included coronary artery disease as well as hypertension and diabetes, and these were similar for both groups. Biochemical analysis showed that the Stage B HF group had higher levels of NT-proBNP (Figure 1A, 117.35 ± 114.93 vs. 78.30 ± 87.55, p = 0.040) and adiponectin (13.13 ± 10.32 vs. 8.74 ± 8.86, p = 0.011). Our results demonstrated that elevated serum levels of NT-proBNP and adiponectin correlated with the progression of HF in asymptomatic patients. Despite the beneficial effects of adiponectin on cardiometabolic traits, increased adiponectin has also been found during HF progression, the so-called adiponectin paradox [15]. Echocardiography analysis showed that both groups had adequate LV systolic function and right ventricle systolic function, as shown by the LVEF, global longitudinal strain analysis, and tricuspid annular plane systolic excursion parameters. Since Stage B HF is defined as having LV hypertrophy, dilatation, increased relative wall thickness, and/or LVdysfunction, the corresponding echocardiographic parameters were greater in the Stage B group as would be expected (Table 2; Figure 1B). Following multiple logistic regression analyses, it was found that an elevated NT-proBNP level was associated with a prediction of Stage B HF (odds ratio: 1.005, 95% confidence interval 1.000–1.010, p = 0.032; Table 3).

3.2. AI-Assisted Identification of SNP Signature: Model Selection, Performance, and Validation

Genome-wide association studies (GWAS) have been able to identify thousands of SNPs that are linked to complex human diseases. To categorize loci associated with asymptomatic HF, namely Stage B, we combined AI-assisted analysis with GWAS using the data from the 117 high-risk (ASCVD risk ≥ 20%) subjects. This study comprised three parts, which were as follows: (1) the performance of GWAS on 83 subjects with Stage A HF and 34 subjects with asymptomatic Stage B HF using the Axiom Genome-Wide TWB 2.0 Array; (2) feature selection; and (3) model derivation and validation. During the feature selection process, only SNPs with adjusted p-values < 1 × 10−7 (0.05/392885) and an odds ratio of > 1 in 392,885 SNPs, whose p-values were calculated by χ2 test and adjusted with Bonferroni correction, were used for the correlation analysis of the HF stage. Multiple feature importance methods within the various different machine-learning algorithms consisted of feature importance in the random forest (RF) approach, weighted support vector in the support vector machine (SVM) approach, and a shrinkage coefficient >0 in the least absolute shrinkage and selection operator (LASSO) approach; these were integrated into the ranking features importance. For model derivation and validation, subjects were randomly assigned into either a training (80%) or a validation (20%) set after selection of the significant SNPs by these three machine-learning (RF, SVM, and LASSO) methods. The above process was repeated 100 times using different combinations of subject to form the training and validation sets. The SNPs were ranked using their summarized counts from the 100-times random samples obtained from the three machine-learning methods (Figure 2A). Although we didn’t have a large sample size, we used three screening processes to obtain reliable results. Firstly, the adjusted p-value (Bonferroni correction) was used to identify the SNPs with strong evidence. Secondly, the positive odds ratio was used to provide reasonable explanations on the risk alleles of the disease. Thirdly, the bootstrap (resampling technique) method with multiple replication times was used to approximate the true population distribution. In addition, different variable selection methods, including machine-learning and statistical methods, were used to rank the importance of SNPs. Notably, our data revealed that the top-ranked SNPs were frequently picked up by the machine using the variable selection and bootstrapped random sample (reproducibility) methods. This means that, although the training and validation samples were different in every bootstrap replication, the top-ranked SNPs consistently showed their importance in distinguishing between Stage A and B patients. Of the four machine-learning models tested (RF, SVM, LASSO, and the decision tree), the SVM model gave the best performance in terms of AUC and accuracy rate. The top 20 SNPs selected by AI had the best prediction performance and the accuracy rate was 0.899 while the AUC was 0.931 when differentiating between Stage A and Stage B (Table 4; Figure 2B–D). Among the top 20 SNPs, 13 SNPs had previously been annotated and mapped either to protein-coding genes, pseudogenes, or non-coding RNA (ncRNA) genes. Furthermore, using these annotated 13 SNPs, the prediction still performed well with an accuracy rate of 0.857 and AUC of 0.912 when differentiating between Stage A and B subjects (Figure 2E).

3.3. The AI-Selected SNP Signature Identifies Asymptomatic High-Risk Subjects That Are Predisposed towards Progression to Heart Failure

The haplotypes distribution, expression patterns, and functions of the genes containing the 13 SNPs that were identified using the Stage A and B subjects are summarized in Table 5. These thirteen genes include six protein-coding genes, one pseudogene, and six ncRNA genes. The SNPs in the five protein-coding genes (GAD2, APP, RASGEF1C, MACROD2, and DMD) and one pseudogene (PGAM1P5) are intron variants. The SNP located in the protein-coding gene DOCK1 is an upstream variant. When the protein-coding genes were examined, we were able to make the following observations (Figure 3A): (1) in the GAD2, APP, and RASGEF1C genes, the incidences of the G/G genotypes were significantly increased in the Stage B subjects (p < 0.001); (2) in the MACROD2 gene, the incidence of the A/A genotype was significantly increased in the Stage B subjects (p < 0.001); and (3) in the DMD and DOCK1 genes, the incidences of the C/T genotypes were significantly increased in the Stage B subjects (p < 0.001). Notably, expression of these four protein-coding genes, namely RASGEF1C, MACROD2, DMD and DOCK1, as well as the pseudogene PGAM1P5, can be detected in the heart [16], suggesting that these candidate genes might play a role in the pathogenesis of HF.
When the non-coding RNA (ncRNA) genes were examined, two SNPs were intron variants of long intergenic ncRNA genes (LINC01968 and LINC00687), while the remain-ing four SNPs were intron variants of ncRNA genes that have uncharacterized functions (LOC105372209, LOC101928047, LOC105372208, and LOC105371356). When the ncRNA genes were examined, we were able to make the following observations (Figure 3B): (1) in LINC01968 and LOC105372208, the C/T genotypes are significantly increased in Stage B subjects (p < 0.001); (2) in LINC00687, the incidence of the C/C genotype is significantly increased in Stage B subjects (p < 0.001); (3) in LOC105372209, LOC105372210, and LOC105371356, the A/G genotypes are significantly increased in Stage B subjects (p < 0.001); and (4) in the LOC101928047 and LOC101928004, the incidences of the T/T geno-types are significantly increased in Stage B subjects (p < 0.001). It should also be noted that the expression of three of these ncRNA, namely LOC105372209, LOC101928047, and LOC105372208, are detectable in the heart [17], which suggests that they might also be involved in HF progression.

3.4. The Protein Interaction Network of the Genes for which the AI-Assistance Identified where the SNPs Are Located

To explore the potential role in the pathogenesis of the high-risk Stage B subjects of the identified protein-coding genes where the SNPs are located, we conducted a protein–protein interaction network analysis in order to identify any potential pathways that might affect cardiac function. Interestingly, five hub proteins, namely TERF1, TERF2, TRIM25, KIAA1429, and PRKACA, were identified in the protein–protein interaction network (marked in blue). This indicates that these proteins are likely to serve as hubs to connect DMD, RASGEF1C, MACROD2, DOCK1, and PGAM1P5 in the heart (Figure 4). In addition, within the protein interaction network, DMD seems to connect with APP, which is expressed in the brain, whereas DOCK1 seems to connect with GAD2, which is expressed in the pancreas (Figure 4). Currently, we cannot rule out the possibility that these two interactions might occur in vivo in cardiomyocytes via a novel mechanism that involves the transportation of the proteins or peptides outside the heart. For example, Aβ, the product of APP, might be released from the brain and delivered to the heart via the circulation where it can interact with DMD, thereby affecting cardiac function.

4. Discussion

The need to halt the progression of asymptomatic pre-clinical HF cannot be over emphasized and thus there is an urgent need for new diagnostic and management tools. The early identification of Stage B HF subjects can be challenging, despite the fact that there is a clear association between traditional risk factors and the development of HF [18]. Notwithstanding the above, the majority of individuals with hypertension and prior myocardial infarction do not eventually develop new-onset HF, a concept that is often referred to as the “prevention paradox”. In this study, we have shown that AI-assisted machine-learning is able to identify SNPs that are potentially associated with the risk of Stage B preclinical HF in high-risk individuals.

4.1. The Genomics of Heart Failure

The presence of heritable, polygenic components related to symptomatic and asymptomatic cardiovascular disease have been long recognized [19]. Many loci are associated with cardiovascular risk factors and diseases, and these have provided insights into the possible mechanisms that underlie the disease. However, there remains a great need for efforts to translate these genomic mechanisms into clinical practice. The application of GWAS has made the identification of loci possible, which is possibly related to the occurrence of HF and the mortality associated with HF. However, due to the heterogeneous pattern of HF, very few GWAS studies have been able to be replicated. Based on the results of GWAS studies, most cardiovascular diseases seem to be influenced by a large number of loci and these variants themselves are seldom the causal variants of the disease. Interestingly, the aggregation of these minor loci account for 10% to 36% of the inherited variation in hyperlipidemia [20], type 2 diabetes [21], myocardial infarction [22], and HF [23,24]. In this study, we hypothesized that genetic factors have an important role in the progression of asymptomatic HF in high-risk subjects. Using a traditional approach, our study subjects were sieved and we confirmed that the high-risk subjects were without any identifiable clinical manifestations of HF. The subjects included in this study were homogenous without other obvious cardiac diseases. In previous studies, the lack of homogeneity has been a potential bias. Our results provide a possible precision medicine approach to the early identification of individuals who are asymptomatic but at high risk of HF. Current guidelines are lacking on how to precisely predict the progression to HF and there is a lack of specific preventive measures or treatments that can used to help these asymptomatic high-risk patients.

4.2. Integrate Artificial Intelligence into the Traditional Prediction Model

In the era of precision medicine, AI or machine-learning algorithms have a number of advantages over the traditional regression model approach [25,26]. Currently, machine-learning or AI is utilized for predicting the prognosis for HF and this has used large or multifaceted datasets such as electronic health records [27] or multi-omics data [28]. However, these datasets are biased due to heterogeneous causes of HF. These prediction models created a “black-box” algorithm that shows only limited improvement over traditional logistic regression prediction models [29]. The application of AI (machine-learning) needs to be incorporated into the traditional approach rather than there being a case of mutual exclusion. In this study, we screened high-risk asymptomatic subjects using a traditional scoring system and this was then followed by AI-assisted prediction using whole genome SNPs. In previous studies, the AUCs of the prognosis prediction, generated by their identified SNPs, combined with traditional risk factors ranged from 0.56 to 0.77 [30,31,32]. Our approach using traditional risk factors, AI-assisted analysis, and a combination of 13 novel SNPs gave an AUC of 0.91. Our findings demonstrate that using an integrated application of traditional and AI-assisted approaches dramatically improves HF prediction.

4.3. Individual SNPs and Heart Failure Progression

Although prior studies using subjects with European ancestry have identified associations between candidate SNPs in the introns of PITX2, ABO, ACTN2, MYOZ1, SYNPO2L, BAG3, and CDKN1A [23,33], such an analysis of a combined database constructed from a number of different cohorts results in an increase in the heterogeneity of the etiology and clinical manifestation of HF; this then leads to a reduction in the statistical power [33]. Other studies focusing on a Han population found an association between a prognosis of heterogeneous cardiomyopathy with SNPs associated with LGALS3 [34], AGCT, SLC25A13, HRG, APOB, SOD3, SYNM, and TLN2 [30]. In the present study, three organ-specific clusters of SNPs were identified (Figure 4). Specifically, APP and GAD2 are mainly expressed in the brain and pancreas, respectively, while the other SNPs, namely DMD, MACROD2, PGAM1P5, RASGEF1C, and DOCK1, are mainly expressed in the heart. The GAD2 polymorphism has been reported to be associated with eating behaviors among women [35] and the risk of obesity [36,37]. Obesity is a well-known risk factor for a HF progression. Another gene associated with obesity found in our signature is the APP gene, which is upregulated in mitochondria and regulates mitochondrial function [38]. Another gene, MACROD2, which is one of three mono ADP-ribosylases in humans, has been reported to act as a transcriptional regulator of adipogenesis and obesity in a Han population [39]. The DOCK1 gene, an atypical Rac activator, has been associated with obesity in a Yup’ik population [40] and is required for cardiovascular development [41]. Mutation of DMD results in muscular dystrophy that can be complicated by the presence of HF and irregular heart rhythms [42]. Finally, there is increasing evidence that lncRNAs are able to affect the expression of protein-coding genes by competitively binding to shared miRNAs, which then reduces the degradation of protein-coding genes. Several studies have found that lncRNA-associated competing endogenous RNA cross-talks with cardiovascular disease pathogenetic processes [43,44].

4.4. Natriuretic Peptides and the Brain

Cardiovascular disorders share many risk factors with Alzheimer’s disease and other memory disorders. NT-proBNP has been found to be an independent risk marker for the incident of dementia and Alzheimer’s disease [45], with higher levels of NT-proBNP being significantly associated with a smaller total grey matter volume [46]. One possible explanation for the relationship between NT-proBNP and dementia may be that individuals with an elevated level of NT-proBNP are more likely to suffer from clinically identified and silent brain ischemic events. It is interesting to note that changes in NT-proBNP are still associated with dementia even after adjusting for CVD risk factors and stroke [47]. Another plausible explanation is that NT-proBNP has a role as a marker of myocardial stress, whereby it reflects the mechanisms leading to progressive subclinical cardiac dysfunction with concomitant myocardial [48] and retinal microvascular damage [49]. Finally, it should be noted that the expression of NT-proBNP is also elevated during stroke and has been found to be associated with increased mortality from stroke [50].

4.5. The Cardiac Natriuretic Peptide System

In our study population, we found that Stage B subjects had a higher level of NT-proBNP. It is well known that cardiac hormones and their prohormones are involved in cardiovascular hemostasis via the regulation of natriuresis, diuresis, vasodilatation, and the inhibition of the renin–angiotensin–aldosterone system (RAAS). BNP is a natriuretic hormone that was initially identified in the brain but is released mainly from the heart, particularly as a response when a ventricle is subject to high ventricular filling pressure [51]. Cleavage of the prohormone pro-BNP produces two forms; these are the biologically active 32 amino acid BNP and the biologically inactive 76 amino acid NT-proBNP. These natriuretic peptides play an important role in the diagnosis of patients who are suspected to have HF [52]. It has been shown previously that when an individual is at risk of HF, BNP-based screening and collaborative care is able to reduce the combined rates of LV systolic-diastolic dysfunction and clinical HF [53]. In addition, NT-proBNP-guided RAAS antagonists and beta-blocker therapy in diabetic subjects have been shown to be beneficial and to help prevent cardiac events [54].

4.6. The Heart-to-Brain Connection

Alzheimer’s disease (AD) and HF with a preserved fraction are age-related disorders that can coexist; they also have common risk factors, a similar epidemiological stratification, and involve common triggers, including oxidative stress, inflammation, and hypoxia. An examination of elderly AD patients has identified the presence of subclinical heart disease, including LV hypertrophy, aortic valve thickening, and aortic regurgitation. The hallmark of AD is the deposition of amyloid plaques, which consist primarily of a 40–42 amino acid peptide called amyloid-β (Aβ). These peptides aggregate into fibrils that then form an ordered β-sheet structure. The amyloid precursor protein is known as the precursor protein for AD-related amyloid Aβ [55]. Aβ deposition in the walls of the cerebral blood vessels is a hallmark lesion of cerebral amyloid angiopathy. In addition, APP and amyloid beta precursor such as protein 2 (APLP2) have also been found to be expressed in cardiomyocytes when heart pathology is present [56]; thus, Aβ may play a role in cardiomyocyte degeneration during HF [57]. Inclusions in the cardiomyocytes of an aging heart are described as being basophilic degenerations of the heart; this has been found to be correlated with age, the degree of myocardial fibrosis in individuals with arterial hypertension, and the severity of cerebral amyloid angiopathy. The fragments detected as part of cardiac basophilic degeneration indicate the presence of specific inclusion body pathology that is related to amyloid precursor protein metabolism. The severity of cerebral amyloid angiopathy has been found to be related to the amyloid precursor protein-derived amyloid β-protein, which suggests a possible link between myocardial and cerebrovascular amyloid precursor protein-related lesions [58].

4.7. Limitations and Future Perspectives

This study has several limitations. First, this is a cross-sectional study involving a limited number of subjects from north-east Taiwan and therefore the results need to be validated using both larger scale cohorts and long-term follow ups; furthermore, the findings may not be applicable to western populations. Second, the 13 SNPs’ signatures, identified by AI-assisted whole-genome SNP analysis, are associated with HF progression in high-risk subjects. Thus, more studies are needed to clarify how these SNPs affect the functions of these genes and how any relevant changes are involved in the underlying mechanisms behind HF progression.

5. Conclusions

This study demonstrates the potential of employing AI machine-learning models to augment traditional methods when predicting genetic predispositions to HF in a high-risk population. Knowledge of the major SNPs associated with preclinical HF provides insights into the relationships between complex pathways and also highlights various key genes that potentially are targets for risk stratification, therapy, and drug development. Our results demonstrate that the application of traditional risk stratification, followed by AI-assisted analysis, is able to raise prediction performance when used on well-characterized, homogenous, and phenotypically identified subjects.

Author Contributions

All the authors contributed to the manuscript preparation. N.-I.Y., C.-H.Y., Y.-J.C., T.-F.T., C.-C.L., and H.-K.S. co-designed the topics and manuscript structure, prepared Figures, and drafted the manuscript. L.-T.K., C.-T.M., Y.-C.S., and M.-J.H. recruited subjects and defined the clinical stage. T.-H.T., C.-H.L., and Y.-H.C. contributed to AI analyses. P.W.-C.H. contributed to the SNPs and protein networking. All authors have read and agreed to the published version of the manuscript.

Funding

We acknowledge the support provided by grants from the Ministry of Health and Welfare (Smart Healthcare for Obesity Therapeutics; PD-109-GP-02 and MG-110-GP-03 to HKS), and from Chang Gung Memorial Hospital (CRRPG2H0181–183 to NIY; CORPG2H0041–0043 and CMRPG2H000091–0093 to CHY; CMRPG2K0141-142 to CCL; and CLRPG2L0051 to YCS).

Institutional Review Board Statement

This study protocol conforms to the ethical guidelines of the 1975 Declaration of Helsinki and was approved by the Institutional Review Board of Chang Gung Medical Foundation (IRB No: 201800802B0 and 202000077B0A3).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

Due to ethical restrictions, the data presented in this study are available on request from the corresponding author.

Acknowledgments

We thank the following core facilities. Firstly, the National Center for Genome Medicine for the technical support in the genotyping. Secondly, we would also like to acknowledge the participant recruitment and sample preservation carried out as part of the Northeastern Taiwan Community Medicine Research Cohort (NTCMRC, ClinicalTrials.gov Identifier: NCT04839796).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Schocken, D.D.; Benjamin, E.J.; Fonarow, G.C.; Krumholz, H.M.; Levy, D.; Mensah, G.A.; Narula, J.; Shor, E.S.; Young, J.B.; Hong, Y.; et al. Prevention of heart failure: A scientific statement from the American Heart Association Councils on Epidemiology and Prevention, Clinical Cardiology, Cardiovascular Nursing, and High Blood Pressure Research; Quality of Care and Outcomes Research Interdisciplinary Working Group; and Functional Genomics and Translational Biology Interdisciplinary Working Group. Circulation 2008, 117, 2544–2565. [Google Scholar] [CrossRef] [Green Version]
  2. Yancy, C.W.; Jessup, M.; Bozkurt, B.; Butler, J.; Casey, D.E., Jr.; Drazner, M.H.; Fonarow, G.C.; Geraci, S.A.; Horwich, T.; Januzzi, J.L.; et al. 2013 ACCF/AHA guideline for the management of heart failure: A report of the American College of Cardiology Foundation/American Heart Association Task Force on Practice Guidelines. J. Am. Coll. Cardiol. 2013, 62, e147–e239. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Ammar, K.A.; Jacobsen, S.J.; Mahoney, D.W.; Kors, J.A.; Redfield, M.M.; Burnett, J.C., Jr.; Rodeheffer, R.J. Prevalence and prognostic significance of heart failure stages: Application of the American College of Cardiology/American Heart Association heart failure staging criteria in the community. Circulation 2007, 115, 1563–1570. [Google Scholar] [CrossRef] [PubMed]
  4. Tavira, B.; Gomez, J.; Ortega, F.; Tranche, S.; Diaz-Corte, C.; Alvarez, F.; Ortiz, A.; Santos, F.; Sanchez-Nino, M.D.; Coto, E. A CLCNKA polymorphism (rs10927887; p.Arg83Gly) previously linked to heart failure is associated with the estimated glomerular filtration rate in the RENASTUR cohort. Gene 2013, 527, 670–672. [Google Scholar] [CrossRef]
  5. Fatini, C.; Sticchi, E.; Marcucci, R.; Verdiani, V.; Nozzoli, C.; Vassallo, C.; Emdin, M.; Abbate, R.; Gensini, G.F. S38G single-nucleotide polymorphism at the KCNE1 locus is associated with heart failure. Heart Rhythm 2010, 7, 363–367. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Tripoliti, E.E.; Papadopoulos, T.G.; Karanasiou, G.S.; Naka, K.K.; Fotiadis, D.I. Heart Failure: Diagnosis, Severity Estimation and Prediction of Adverse Events Through Machine Learning Techniques. Comput. Struct. Biotechnol. J. 2017, 15, 26–47. [Google Scholar] [CrossRef] [Green Version]
  7. Yu, S.-N.; Lee, M.-Y. Bispectral analysis and genetic algorithm for congestive heart failure recognition based on heart rate variability. Comput. Biol. Med. 2012, 42, 816–825. [Google Scholar] [CrossRef] [PubMed]
  8. Olsen, C.R.; Mentz, R.J.; Anstrom, K.J.; Page, D.; Patel, P.A. Clinical applications of machine learning in the diagnosis, classification, and prediction of heart failure. Am. Heart J. 2020, 229, 1–17. [Google Scholar] [CrossRef]
  9. Zheng, Y.; Guo, X.; Qin, J.; Xiao, S. Computer-assisted diagnosis for chronic heart failure by the analysis of their cardiac reserve and heart sound characteristics. Comput. Methods Programs Biomed. 2015, 122, 372–383. [Google Scholar] [CrossRef]
  10. Angraal, S.; Mortazavi, B.J.; Gupta, A.; Khera, R.; Ahmad, T.; Desai, N.R.; Jacoby, D.L.; Masoudi, F.A.; Spertus, J.A.; Krumholz, H.M. Machine Learning Prediction of Mortality and Hospitalization in Heart Failure With Preserved Ejection Fraction. JACC Heart Fail. 2020, 8, 12–21. [Google Scholar] [CrossRef]
  11. Goff, D.C., Jr.; Lloyd-Jones, D.M.; Bennett, G.; Coady, S.; D’Agostino, R.B.; Gibbons, R.; Greenland, P.; Lackland, D.T.; Levy, D.; O’Donnell, C.J.; et al. 2013 ACC/AHA guideline on the assessment of cardiovascular risk: A report of the American College of Cardiology/American Heart Association Task Force on Practice Guidelines. Circulation 2014, 129, S49–S73. [Google Scholar] [CrossRef] [Green Version]
  12. Yang, N.I.; Wang, C.H.; Hung, M.J.; Chen, Y.C.; Wu, I.W.; Lee, C.C.; Wu, M.S.; Kuo, L.T.; Cheng, C.W.; Cherng, W.J. Real-time three-dimensional echocardiography provides advanced haemodynamic information associated with intra-dialytic hypotension in patients with autonomic dysfunction. Nephrol Dial Transpl. 2010, 25, 249–254. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Oughtred, R.; Rust, J.; Chang, C.; Breitkreutz, B.J.; Stark, C.; Willems, A.; Boucher, L.; Leung, G.; Kolas, N.; Zhang, F.; et al. The BioGRID database: A comprehensive biomedical resource of curated protein, genetic, and chemical interactions. Protein Sci. 2021, 30, 187–200. [Google Scholar] [CrossRef] [PubMed]
  14. Shannon, P.; Markiel, A.; Ozier, O.; Baliga, N.S.; Wang, J.T.; Ramage, D.; Amin, N.; Schwikowski, B.; Ideker, T. Cytoscape: A software environment for integrated models of biomolecular interaction networks. Genome Res. 2003, 13, 2498–2504. [Google Scholar] [CrossRef]
  15. Menzaghi, C.; Trischitta, V. The Adiponectin Paradox for All-Cause and Cardiovascular Mortality. Diabetes 2018, 67, 12–22. [Google Scholar] [CrossRef] [Green Version]
  16. Uhlen, M.; Fagerberg, L.; Hallstrom, B.M.; Lindskog, C.; Oksvold, P.; Mardinoglu, A.; Sivertsson, A.; Kampf, C.; Sjostedt, E.; Asplund, A.; et al. Proteomics. Tissue-based map of the human proteome. Science 2015, 347, 1260419. [Google Scholar] [CrossRef]
  17. Fagerberg, L.; Hallstrom, B.M.; Oksvold, P.; Kampf, C.; Djureinovic, D.; Odeberg, J.; Habuka, M.; Tahmasebpoor, S.; Danielsson, A.; Edlund, K.; et al. Analysis of the human tissue-specific expression by genome-wide integration of transcriptomics and antibody-based proteomics. Mol. Cell. Proteom. 2014, 13, 397–406. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Echouffo-Tcheugui, J.B.; Greene, S.J.; Papadimitriou, L.; Zannad, F.; Yancy, C.W.; Gheorghiade, M.; Butler, J. Population risk prediction models for incident heart failure: A systematic review. Circ. Heart Fail. 2015, 8, 438–447. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. O’Donnell, C.J.; Nabel, E.G. Genomics of cardiovascular disease. N. Engl. J. Med. 2011, 365, 2098–2109. [Google Scholar] [CrossRef]
  20. Teslovich, T.M.; Musunuru, K.; Smith, A.V.; Edmondson, A.C.; Stylianou, I.M.; Koseki, M.; Pirruccello, J.P.; Ripatti, S.; Chasman, D.I.; Willer, C.J.; et al. Biological, clinical and population relevance of 95 loci for blood lipids. Nature 2010, 466, 707–713. [Google Scholar] [CrossRef]
  21. Voight, B.F.; Scott, L.J.; Steinthorsdottir, V.; Morris, A.P.; Dina, C.; Welch, R.P.; Zeggini, E.; Huth, C.; Aulchenko, Y.S.; Thorleifsson, G.; et al. Twelve type 2 diabetes susceptibility loci identified through large-scale association analysis. Nat. Genet. 2010, 42, 579–589. [Google Scholar] [CrossRef]
  22. Schunkert, H.; Konig, I.R.; Kathiresan, S.; Reilly, M.P.; Assimes, T.L.; Holm, H.; Preuss, M.; Stewart, A.F.; Barbalic, M.; Gieger, C.; et al. Large-scale association analysis identifies 13 new susceptibility loci for coronary artery disease. Nat. Genet. 2011, 43, 333–338. [Google Scholar] [CrossRef]
  23. Arvanitis, M.; Tampakakis, E.; Zhang, Y.; Wang, W.; Auton, A.; 23andMe Research Team; Dutta, D.; Glavaris, S.; Keramati, A.; Chatterjee, N.; et al. Genome-wide association and multi-omic analyses reveal ACTN2 as a gene linked to heart failure. Nat. Commun. 2020, 11, 1122. [Google Scholar] [CrossRef] [Green Version]
  24. Lindgren, M.P.; PirouziFard, M.; Smith, J.G.; Sundquist, J.; Sundquist, K.; Zoller, B. A Swedish Nationwide Adoption Study of the Heritability of Heart Failure. JAMA Cardiol. 2018, 3, 703–710. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Goldstein, B.A.; Navar, A.M.; Carter, R.E. Moving beyond regression techniques in cardiovascular risk prediction: Applying machine learning to address analytic challenges. Eur. Heart J. 2017, 38, 1805–1814. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Sajeev, S.; Maeder, A.; Champion, S.; Beleigoli, A.; Ton, C.; Kong, X.; Shu, M. Deep Learning to Improve Heart Disease Risk Prediction; Springer: Cham, Switzerland, 2019; pp. 96–103. [Google Scholar]
  27. Bots, S.H.; Onland-Moret, N.C.; Tulevski, II; van der Harst, P.; Cramer, M.J.M.; Asselbergs, F.W.; Somsen, G.A.; den Ruijter, H.M. Heart failure medication dosage and survival in women and men seen at outpatient clinics. Heart 2021. [Google Scholar] [CrossRef] [PubMed]
  28. Goldstein, B.A.; Hubbard, A.E.; Cutler, A.; Barcellos, L.F. An application of Random Forests to a genome-wide association dataset: Methodological considerations & new findings. BMC Genet. 2010, 11, 49. [Google Scholar] [CrossRef] [Green Version]
  29. Desai, R.J.; Wang, S.V.; Vaduganathan, M.; Evers, T.; Schneeweiss, S. Comparison of Machine Learning Methods With Traditional Models for Use of Administrative Claims with Electronic Medical Records to Predict Heart Failure Outcomes. JAMA Netw. Open 2020, 3, e1918962. [Google Scholar] [CrossRef]
  30. Li, S.; Sun, Y.; Hu, S.; Hu, D.; Li, C.; Xiao, L.; Chen, Y.; Li, H.; Cui, G.; Wang, D.W. Genetic risk scores to predict the prognosis of chronic heart failure patients in Chinese Han. J. Cell. Mol. Med. 2020, 24, 285–293. [Google Scholar] [CrossRef] [Green Version]
  31. Davies, R.W.; Dandona, S.; Stewart, A.F.; Chen, L.; Ellis, S.G.; Tang, W.H.; Hazen, S.L.; Roberts, R.; McPherson, R.; Wells, G.A. Improved prediction of cardiovascular disease based on a panel of single nucleotide polymorphisms identified through genome-wide association studies. Circ. Cardiovasc. Genet. 2010, 3, 468–474. [Google Scholar] [CrossRef] [Green Version]
  32. Lello, L.; Raben, T.G.; Yong, S.Y.; Tellier, L.; Hsu, S.D.H. Genomic Prediction of 16 Complex Disease Risks Including Heart Attack, Diabetes, Breast and Prostate Cancer. Sci. Rep. 2019, 9, 15286. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Shah, S.; Henry, A.; Roselli, C.; Lin, H.; Sveinbjornsson, G.; Fatemifar, G.; Hedman, A.K.; Wilk, J.B.; Morley, M.P.; Chaffin, M.D.; et al. Genome-wide association and Mendelian randomisation analysis provide insights into the pathogenesis of heart failure. Nat. Commun. 2020, 11, 163. [Google Scholar] [CrossRef]
  34. Zhang, Y.; Wang, Y.; Zhai, M.; Gan, T.; Zhao, X.; Zhang, R.; An, T.; Huang, Y.; Zhou, Q.; Zhang, J. Influence of LGALS3 gene polymorphisms on susceptibility and prognosis of dilated cardiomyopathy in a Northern Han Chinese population. Gene 2018, 642, 293–298. [Google Scholar] [CrossRef]
  35. Choquette, A.C.; Lemieux, S.; Tremblay, A.; Drapeau, V.; Bouchard, C.; Vohl, M.C.; Perusse, L. GAD2 gene sequence variations are associated with eating behaviors and weight gain in women from the Quebec family study. Physiol. Behav. 2009, 98, 505–510. [Google Scholar] [CrossRef]
  36. Chen, K.C.; Lin, Y.C.; Chao, W.C.; Chung, H.K.; Chi, S.S.; Liu, W.S.; Wu, W.T. Association of genetic polymorphisms of glutamate decarboxylase 2 and the dopamine D2 receptor with obesity in Taiwanese subjects. Ann. Saudi Med. 2012, 32, 121–126. [Google Scholar] [CrossRef]
  37. Prakash, J.; Mittal, B.; Awasthi, S.; Srivastava, N. Association of the −243A>G, +61450C>A Polymorphisms of the Glutamate Decarboxylase 2 (GAD2) Gene with Obesity and Insulin Level in North Indian Population. Iran J. Public Health 2016, 45, 460–468. [Google Scholar] [PubMed]
  38. An, Y.A.; Crewe, C.; Asterholm, I.W.; Sun, K.; Chen, S.; Zhang, F.; Shao, M.; Funcke, J.B.; Zhang, Z.; Straub, L.; et al. Dysregulation of Amyloid Precursor Protein Impairs Adipose Tissue Mitochondrial Function and Promotes Obesity. Nat. Metab. 2019, 1, 1243–1257. [Google Scholar] [CrossRef] [PubMed]
  39. Chang, Y.C.; Hee, S.W.; Lee, W.J.; Li, H.Y.; Chang, T.J.; Lin, M.W.; Hung, Y.J.; Lee, I.T.; Hung, K.Y.; Assimes, T.; et al. Genome-wide scan for circulating vascular adhesion protein-1 levels: MACROD2 as a potential transcriptional regulator of adipogenesis. J. Diabetes Investig. 2018, 9, 1067–1074. [Google Scholar] [CrossRef] [Green Version]
  40. Vaughan, L.K.; Wiener, H.W.; Aslibekyan, S.; Allison, D.B.; Havel, P.J.; Stanhope, K.L.; O’Brien, D.M.; Hopkins, S.E.; Lemas, D.J.; Boyer, B.B.; et al. Linkage and association analysis of obesity traits reveals novel loci and interactions with dietary n-3 fatty acids in an Alaska Native (Yup’ik) population. Metabolism 2015, 64, 689–697. [Google Scholar] [CrossRef] [Green Version]
  41. Sanematsu, F.; Hirashima, M.; Laurin, M.; Takii, R.; Nishikimi, A.; Kitajima, K.; Ding, G.; Noda, M.; Murata, Y.; Tanaka, Y.; et al. DOCK180 is a Rac activator that regulates cardiovascular development by acting downstream of CXCR4. Circ. Res. 2010, 107, 1102–1105. [Google Scholar] [CrossRef]
  42. Gao, Q.Q.; McNally, E.M. The Dystrophin Complex: Structure, Function, and Implications for Therapy. Compr. Physiol. 2015, 5, 1223–1239. [Google Scholar] [CrossRef] [Green Version]
  43. Xue, X.; Luo, L. LncRNA HIF1A-AS1 contributes to ventricular remodeling after myocardial ischemia/reperfusion injury by adsorption of microRNA-204 to regulating SOCS2 expression. Cell Cycle 2019, 18, 2465–2480. [Google Scholar] [CrossRef]
  44. Wang, Q.C.; Wang, Z.Y.; Xu, Q.; Chen, X.L.; Shi, R.Z. lncRNA expression profiles and associated ceRNA network analyses in epicardial adipose tissue of patients with coronary artery disease. Sci. Rep. 2021, 11, 1567. [Google Scholar] [CrossRef]
  45. Tynkkynen, J.; Hernesniemi, J.A.; Laatikainen, T.; Havulinna, A.S.; Salo, P.; Blankenberg, S.; Zeller, T.; Salomaa, V. High-sensitivity cardiac troponin I and NT-proBNP as predictors of incident dementia and Alzheimer’s disease: The FINRISK Study. J. Neurol. 2017, 264, 503–511. [Google Scholar] [CrossRef]
  46. Ferguson, I.T.; Elbejjani, M.; Sabayan, B.; Jacobs, D.R., Jr.; Meirelles, O.; Sanchez, O.A.; Tracy, R.; Bryan, N.; Launer, L.J. N-Terminal pro-Brain Natriuretic Peptide and Associations With Brain Magnetic Resonance Imaging (MRI) Features in Middle Age: The CARDIA Brain MRI Study. Front. Neurol. 2018, 9, 307. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Mirza, S.S.; de Bruijn, R.F.; Koudstaal, P.J.; van den Meiracker, A.H.; Franco, O.H.; Hofman, A.; Tiemeier, H.; Ikram, M.A. The N-terminal pro B-type natriuretic peptide, and risk of dementia and cognitive decline: A 10-year follow-up study in the general population. J. Neurol. Neurosurg. Psychiatry 2016, 87, 356–362. [Google Scholar] [CrossRef]
  48. Mitchell, A.; Misialek, J.R.; Folsom, A.R.; Duprez, D.; Alonso, A.; Jerosch-Herold, M.; Sanchez, O.A.; Watson, K.E.; Sallam, T.; Konety, S.H. Usefulness of N-terminal Pro-brain Natriuretic Peptide and Myocardial Perfusion in Asymptomatic Adults (from the Multi-Ethnic Study of Atherosclerosis). Am. J. Cardiol. 2015, 115, 1341–1345. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Mutlu, U.; Ikram, M.A.; Hofman, A.; de Jong, P.T.; Klaver, C.C.; Ikram, M.K. N-Terminal Pro-B-Type Natriuretic Peptide Is Related to Retinal Microvascular Damage: The Rotterdam Study. Arterioscler. Thromb. Vasc. Biol. 2016, 36, 1698–1702. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Tu, W.J.; Ma, G.Z.; Ni, Y.; Hu, X.S.; Luo, D.Z.; Zeng, X.W.; Liu, Q.; Xu, T.; Yu, L.; Wu, B. Copeptin and NT-proBNP for prediction of all-cause and cardiovascular death in ischemic stroke. Neurology 2017, 88, 1899–1905. [Google Scholar] [CrossRef] [PubMed]
  51. Kinnunen, P.; Vuolteenaho, O.; Ruskoaho, H. Mechanisms of atrial and brain natriuretic peptide release from rat ventricular myocardium: Effect of stretching. Endocrinology 1993, 132, 1961–1970. [Google Scholar] [CrossRef] [PubMed]
  52. Ponikowski, P.; Voors, A.A.; Anker, S.D.; Bueno, H.; Cleland, J.G.F.; Coats, A.J.S.; Falk, V.; Gonzalez-Juanatey, J.R.; Harjola, V.P.; Jankowska, E.A.; et al. 2016 ESC Guidelines for the diagnosis and treatment of acute and chronic heart failure: The Task Force for the diagnosis and treatment of acute and chronic heart failure of the European Society of Cardiology (ESC)Developed with the special contribution of the Heart Failure Association (HFA) of the ESC. Eur. Heart J. 2016, 37, 2129–2200. [Google Scholar] [CrossRef] [PubMed]
  53. Ledwidge, M.; Gallagher, J.; Conlon, C.; Tallon, E.; O’Connell, E.; Dawkins, I.; Watson, C.; O’Hanlon, R.; Bermingham, M.; Patle, A.; et al. Natriuretic peptide-based screening and collaborative care for heart failure: The STOP-HF randomized trial. JAMA 2013, 310, 66–74. [Google Scholar] [CrossRef] [PubMed]
  54. Huelsmann, M.; Neuhold, S.; Resl, M.; Strunk, G.; Brath, H.; Francesconi, C.; Adlbrecht, C.; Prager, R.; Luger, A.; Pacher, R.; et al. PONTIAC (NT-proBNP selected prevention of cardiac events in a population of diabetic patients without a history of cardiac disease): A prospective randomized controlled trial. J. Am. Coll. Cardiol. 2013, 62, 1365–1372. [Google Scholar] [CrossRef] [Green Version]
  55. Kang, J.; Lemaire, H.G.; Unterbeck, A.; Salbaum, J.M.; Masters, C.L.; Grzeschik, K.H.; Multhaup, G.; Beyreuther, K.; Muller-Hill, B. The precursor of Alzheimer’s disease amyloid A4 protein resembles a cell-surface receptor. Nature 1987, 325, 733–736. [Google Scholar] [CrossRef] [PubMed]
  56. Slunt, H.H.; Thinakaran, G.; Von Koch, C.; Lo, A.C.; Tanzi, R.E.; Sisodia, S.S. Expression of a ubiquitous, cross-reactive homologue of the mouse beta-amyloid precursor protein (APP). J. Biol. Chem. 1994, 269, 2637–2644. [Google Scholar] [CrossRef]
  57. Greco, S.; Zaccagnini, G.; Fuschi, P.; Voellenkle, C.; Carrara, M.; Sadeghi, I.; Bearzi, C.; Maimone, B.; Castelvecchio, S.; Stellos, K.; et al. Increased BACE1-AS long noncoding RNA and beta-amyloid levels in heart failure. Cardiovasc. Res. 2017, 113, 453–463. [Google Scholar] [CrossRef]
  58. Kramer, L.M.; Brettschneider, J.; Lennerz, J.K.; Walcher, D.; Fang, L.; Rosenbohm, A.; Balakrishnan, K.; Benckendorff, J.; Moller, P.; Just, S.; et al. Amyloid precursor protein-fragments-containing inclusions in cardiomyocytes with basophilic degeneration and its association with cerebral amyloid angiopathy and myocardial fibrosis. Sci. Rep. 2018, 8, 16594. [Google Scholar] [CrossRef]
Figure 1. Representative echocardiographic images and the NT-proBNP levels of subjects in Stage A and B heart failure (HF). (A) The NT-proBNP levels of the subjects in Stage A and B HF. * indicated p < 0.05 by independent t-test. (B) Upper panel: a 68-year-old male subject with Stage A HF (body surface area (BSA) of 1.8 m2). His echocardiographic parameters, including interventricular septum thickness (IVS; A1, 7 mm, normal range 6–10 mm), left ventricle end diastolic volume index (LVEDVI; A2, 36 mL/m2; normal range < 95 mL/m2), and left atrium volume index (LAVI; A3, 31 ml/m2; normal range < 34 mL/m2) were within normal ranges. Bottom panel: a 70-year-old male subject with asymptomatic Stage B HF (BSA 1.9 m2). The echocardiographic examination results demonstrated myocardial remodeling, including thickened IVS (B1, 13 mm) and enlarged LVEDI (B2, 102 mL/m2) and LAVI (B3, 51 mL/m2). The red arrow indicates that the echocardiographic pictures were taken at the end-diastolic phase in A2 and B2, and at the end-systolic phase in A3 and B3.
Figure 1. Representative echocardiographic images and the NT-proBNP levels of subjects in Stage A and B heart failure (HF). (A) The NT-proBNP levels of the subjects in Stage A and B HF. * indicated p < 0.05 by independent t-test. (B) Upper panel: a 68-year-old male subject with Stage A HF (body surface area (BSA) of 1.8 m2). His echocardiographic parameters, including interventricular septum thickness (IVS; A1, 7 mm, normal range 6–10 mm), left ventricle end diastolic volume index (LVEDVI; A2, 36 mL/m2; normal range < 95 mL/m2), and left atrium volume index (LAVI; A3, 31 ml/m2; normal range < 34 mL/m2) were within normal ranges. Bottom panel: a 70-year-old male subject with asymptomatic Stage B HF (BSA 1.9 m2). The echocardiographic examination results demonstrated myocardial remodeling, including thickened IVS (B1, 13 mm) and enlarged LVEDI (B2, 102 mL/m2) and LAVI (B3, 51 mL/m2). The red arrow indicates that the echocardiographic pictures were taken at the end-diastolic phase in A2 and B2, and at the end-systolic phase in A3 and B3.
Cells 10 02430 g001aCells 10 02430 g001b
Figure 2. Prediction models built based on the clinical features and genome-wide SNPs were used to classify Stage A and Stage B using machine-learning methods. (A) Machine-learning workflow to select reliable and predictable features to build prediction models for the group classification. (B) The accuracy and AUC of the SVM model with different feature selection numbers. (C) Classification performance of the four machine-learning algorithms with selected feature subsets from genome-wide SNPs. (D) The SVM confusion matrix used for prediction accuracy. (E) Classification of selected SNPs.
Figure 2. Prediction models built based on the clinical features and genome-wide SNPs were used to classify Stage A and Stage B using machine-learning methods. (A) Machine-learning workflow to select reliable and predictable features to build prediction models for the group classification. (B) The accuracy and AUC of the SVM model with different feature selection numbers. (C) Classification performance of the four machine-learning algorithms with selected feature subsets from genome-wide SNPs. (D) The SVM confusion matrix used for prediction accuracy. (E) Classification of selected SNPs.
Cells 10 02430 g002
Figure 3. Pie charts indicating the genotype frequencies of the SNPs identified by AI-assisted analysis using the SNP datasets obtained from Stage A and B subjects. (A) Pie charts representing the AI-assisted SNPs within the protein-coding genes, excluding the pseudogene (PGAM1P5). Abbreviations include GAD2: glutamate decarboxylase 2; APP: amyloid beta precursor protein; RASGEF1C: RasGEF domain family member 1C; MACROD2: mono-ADP ribosylhydrolase 2; DMD: dystrophin; and DOCK1: dedicator of cytokinesis 1. (B) Pie charts representing the AI-assisted SNPs within the non-coding genes. rs8084397 is the SNP of LOC105372209 and LOC105372210, and rs201036 is the SNP of LOC101928047 and LOC101928004. # indicates that the signaling by the SNP array was lower than the calling rate. ** p < 0.005 by Chi-square test.
Figure 3. Pie charts indicating the genotype frequencies of the SNPs identified by AI-assisted analysis using the SNP datasets obtained from Stage A and B subjects. (A) Pie charts representing the AI-assisted SNPs within the protein-coding genes, excluding the pseudogene (PGAM1P5). Abbreviations include GAD2: glutamate decarboxylase 2; APP: amyloid beta precursor protein; RASGEF1C: RasGEF domain family member 1C; MACROD2: mono-ADP ribosylhydrolase 2; DMD: dystrophin; and DOCK1: dedicator of cytokinesis 1. (B) Pie charts representing the AI-assisted SNPs within the non-coding genes. rs8084397 is the SNP of LOC105372209 and LOC105372210, and rs201036 is the SNP of LOC101928047 and LOC101928004. # indicates that the signaling by the SNP array was lower than the calling rate. ** p < 0.005 by Chi-square test.
Cells 10 02430 g003
Figure 4. Protein–protein interaction network of the genes containing the AI-assisted identification of SNPs. The protein–protein interaction network was established using the protein-coding genes that contain a SNP. The red dots indicate the genes containing the AI-assisted SNPs and the blue dots indicate the proteins that are connected to the proteins containing SNPs within the network. PGAM1P5 is annotated as a pseudogene. Abbreviations include VIRMA: Vir-like M6A Methyltransferase-associated; PRKACA: protein kinase CAMP-activated catalytic subunit alpha; TERF1: telomeric repeat binding factor 1; TERF2: telomeric repeat binding factor 2; and TRIM25: tripartite motif containing 25.
Figure 4. Protein–protein interaction network of the genes containing the AI-assisted identification of SNPs. The protein–protein interaction network was established using the protein-coding genes that contain a SNP. The red dots indicate the genes containing the AI-assisted SNPs and the blue dots indicate the proteins that are connected to the proteins containing SNPs within the network. PGAM1P5 is annotated as a pseudogene. Abbreviations include VIRMA: Vir-like M6A Methyltransferase-associated; PRKACA: protein kinase CAMP-activated catalytic subunit alpha; TERF1: telomeric repeat binding factor 1; TERF2: telomeric repeat binding factor 2; and TRIM25: tripartite motif containing 25.
Cells 10 02430 g004
Table 1. Baseline characteristics of both groups.
Table 1. Baseline characteristics of both groups.
Stage A HF 1
N = 113
Stage B HF
N = 49
p-Value
Age, years68.35 ± 6.1969.88 ± 6.530.156
Male, n (%)72 (63.7)32 (65.3)1.00
BMI 2, kg/m 226.60 ± 3.7627.21 ± 3.800.342
Smoking, n (%)22 (19.5)12 (24.5)0.53
ASCVD 3 risk score, %24 ± 11.2527.46 ± 11.830.087
Coronary artery disease, n (%)34 (30.1)11 (22.4)0.35
Hypertension, n (%)93 (82.3)44 (89.9)0.343
Diabetes mellitus, n (%)40 (35.4)16 (32.7)0.858
Urine Alb/Cre 4 (mg/g)67.21 ± 173.22175.26 ± 379.360.084
NT-proBNP 5 (pg/mL)78.30 ± 87.55117.35 ± 114.930.040
Hs-Tnt 6 (ng/L)10.37 ± 5.47 12.70 ± 9.67 0.128
HbA1c 7 (%)6.44 ± 9.856.58 ± 1.560.549
Creatinine (mg/dL)1.00 ± 0.271.06 ± 0.400.287
AST 8 (U/L)24.79 ± 11.4430.80 ± 37.520.277
ALT 9 (U/L)29.05 ± 17.1338.19 ± 66.930.351
Sodium (mEq/L)141.41 ± 1.92141.65 ± 2.560.501
Potassium (mEq/L)4.33 ± 0.404.27 ± 0.410.393
Albumin (g/dL)4.58 ± 0.234.52 ± 0.290.194
Total bilirubin (mg/dL)0.68 ± 0.280.66 ± 0.250.561
Total protein (g/dL)7.36 ± 0.407.31 ± 0.360.494
Uric acid (mg/dL)6.03 ± 1.485.73 ± 1.550.252
HDL 10 (mg/dL)51.83 ± 13.4 52.24 ± 15.280.863
LDL 11 (mg/dL)111.71 ± 26.09111.02 ± 26.660.879
Total cholesterol (mg/dL)178.60 ± 28.71180 ± 35.05 0.758
Insulin (μIU/mL)14.29 ± 7.5013.79 ± 8.670.715
Hs-CRP (mg/L) 122.69 ± 11.002.21 ± 2.300.772
Ferritin (ng/mL)303.39 ± 240.80280.78 ± 186.700.559
Adiponectin (μg/mL)8.74 ± 8.8613.13 ± 10.320.011
Leptin (ng/mL)11.47 ± 11.4110.77 ± 9.650.708
1. HF: heart failure; 2 BMI: body mass index; 3 ASCVD: atherosclerotic cardiovascular disease; 4 Alb/Cre: albumin creatinine ratio; 5 NT-proBNP: N-terminal pro b-type natriuretic peptide; 6 Hs-Tnt: high-sensitivity Troponin T; 7 HbA1c: hemoglobin A1c; 8 AST: aspartate aminotransferase; 9 ALT: alanine aminotransferase; 10 HDL: high-density lipoprotein; 11 LDL: low-density lipoprotein; and 12 Hs-CRP: high-sensitivity C-reactive protein.
Table 2. Echocardiographic parameters of both groups.
Table 2. Echocardiographic parameters of both groups.
Stage A HF 1
N = 113
Stage B HF
N = 49
p-Value
Left atrial volume index (mL/m2)34.61 ± 9.6242.92 ± 12.32<0.001
Left ventricle EDVI 2 (mL/m2)54.22 ± 12.2162.39 ± 19.790.009
Left ventricle ESVI 3 (mL/m2)21.10 ± 5.6923.64 ± 9.100.075
Left ventricle EF 4, 2D (%)63.54 ± 6.5863.92 ± 6.180.736
Left ventricle mass index (g/m2)75.76 ± 17.56104.45 ± 28.73<0.001
Relative wall thickness0.28 ± 0.060.29 ± 0.070.201
Peak GLS 5 (%)−17.22 ± 3.06−18.04 ± 3.310.134
Mitral valve E 6 velocity (cm/s)0.69 ± 0.190.74 ± 0.200.207
Mitral valve A 7 velocity (cm/s)0.84 ± 0.230.96 ± 0.200.003
Mitral valve deceleration time (ms)204.45 ± 64.00201.50 ± 47.240.775
Mitral valve E/A ratio1.28 ± 3.360.78 ± 0.210.306
Tissue doppler septal S’ 8 (cm/s)6.28 ± 1.305.98 ± 1.120.177
Tissue doppler septal E’ 9 (cm/s)5.98 ± 1.124.95 ± 1.37<0.001
Septal E/E’11.97 ± 3.7415.26 ± 3.66<0.001
TAPSE 10 (mm)22.46 ± 3.2227.37 ± 34.280.307
1. HF: heart failure; 2 EDVI: end-diastolic volume index; 3 ESVI: end-systolic volume Index; 4 EF: ejection fraction; 5 GLS: global longitudinal strain; 6 E: early diastolic transmitral flow; 7 A: late diastolic transmitral flow; 8 E’: septal annulus mitral early diastolic tissue velocity; 9 S’: septal annulus mitral systolic tissue velocity; and 10 TAPSE: tricuspid annular plane systolic excursion.
Table 3. Multiple logistic regression analysis for Stage B heart failure prediction.
Table 3. Multiple logistic regression analysis for Stage B heart failure prediction.
Odds Ratio95% Confidence Intervalp Value
Male1.4560.543–3.9020.455
Hypertension1.9760.574–6.8050.280
Diabetes Mellitus0.8510.329–2.2010.740
Smoking2.396 0.836–6.8620.104
Urine Alb/Cre 11.0021.000–1.0030.065
NT-proBNP 21.0051.000–1.0100.032
Adiponectin1.0310.986–1.0780.182
GLS 30.9290.809–1.0670.298
1 Alb/Cre: albumin creatinine ratio; 2 NT-proBNP: N-terminal pro b-type natriuretic peptide; and 3 GLS: global longitudinal strain.
Table 4. Top 20 SNPs of both groups.
Table 4. Top 20 SNPs of both groups.
SNPsPositionGeneAlleleStage A HF
N = 113
Stage B HF
N = 49
p-Value
rs1999241chr10:26277498GAD2 A A/G A/G G/0 034/44/4/117/6/11/0<0.001
rs7645985chr3:194777389LINC01968C C/T C/T T78/4/121/13/0<0.001
rs6516709chr21:25893395APPA A/G A/G G37/43/312/11/11<0.001
rs6078354chr20:11820605LINC00687 C C/C T/T T3/22/587/17/10<0.001
rs7725201chr5:180102124RASGEF1CA A/G A/G G45/38/019/9/6<0.001
rs10859918chr12:95652910PGAM1P5G G/T G/T T/0 045/26/11/18/23/3/0<0.001
rs6110516chr20:15097697MACROD2A G/G G/A A13/70/013/18/3<0.001
rs4693641chr4:83755117NoneT T/C C/C T83/0/027/1/6<0.001
rs5928104chrX:32899133DMD C C/T T/C T3/80/00/27/7<0.001
rs8084397chr18:76251014LOC105372209; LOC105372210A A/G A/G G/0 058/19/5/110/21/3/0<0.001
rs4715127chr6:49320956NoneC C/C T/T T6/22/550/22/12<0.001
rs2496369chr6:49156247NoneG G/G T/T T1/25/574/19/11<0.001
Affx-2716217chr10:126685478DOCK1C C/T C/T T53/23/712/21/1<0.001
rs56352414chr22:48917877NoneC C/T C/T T50/32/116/10/8<0.001
rs2806810chr13:103997325NoneA A/C A/C C38/42/39/16/9<0.001
rs4934985chr10:33616333NoneC C/C T/T T/0 05/30/48/012/13/8/1<0.001
rs6912291chr6:110065139NoneC C/C T/T T1/34/486/3/25<0.001
rs201036chr6:6708885LOC101928047; LOC101928004C C/T C/T T32/39/122/17/15<0.001
rs9965164chr18:76220796LOC105372208C C/T C/T T59/19/510/21/3<0.001
rs3813579chr16:79715379LOC105371356A A/G A/G G43/31/95/21/8<0.001
Allele: ‘0’ is no allele appearance.
Table 5. Stage A vs. Stage B: feature importance.
Table 5. Stage A vs. Stage B: feature importance.
RankSNPsTypeGeneGene TypeFunctionExpression
Pattern
1rs1999241Intron variantGAD2Protein-codingMajor autoantigen in insulin-dependent diabetesPancreas and brain
2rs7645985Intron variantLINC01968lcRNA 1N/ATestis and placenta
3rs6516709Intron variantAPPProtein-codingNeurite growth, neuronal adhesion, and axonogenesisBrain
4rs6078354Intron variantLINC00687lcRNA 1N/ATestis
5rs7725201Stop gainedRASGEF1CProtein-codingGuanine nucleotide exchange factor Ubiquitous
6rs10859918Intron variantPGAM1P5PseudoPseudogeneUbiquitous
7rs6110516Intron variantMACROD2Protein-codingRemoving ADP-ribose from mono-ADP-ribosylated proteinsUbiquitous
8rs5928104Intron variantDMDProtein-coding
  • The DGC 3 bridges the inner cytoskeleton and ECM 4
  • Involved in DMD 5, BMD 6, or cardiomyopathy
Ubiquitous
9rs8084397Intron variantLOC105372209 and LOC105372210ncRNA 2N/AHeart, testis, placenta, and brain
10Affx-2716217Upstream variantDOCK1Protein-coding
  • Dedicator of the cytokinesis protein
  • Phagocytosis and cell migration
Ubiquitous
11rs201036Non-coding transcript; intron variantLOC101928047 and LOC101928004ncRNA 2N/AHeart, kidney, bone marrow, and fat
12rs9965164Non-coding transcriptLOC105372208ncRNA 2N/AHeart, testis, and placenta
13rs3813579Intron variantLOC105371356ncRNA 2N/ALiver, testis, kidney, and skin
1 lcRNA: long intergenic non-protein-coding RNA; 2 ncRNA: uncharacterized model non-protein-coding RNA; 3 DGC: dystrophin-glycoprotein complex; 4 ECM: extracellular matrix; 5 DMD: Duchenne muscular dystrophy; and 6 BMD: Becker muscular dystrophy.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Yang, N.-I.; Yeh, C.-H.; Tsai, T.-H.; Chou, Y.-J.; Hsu, P.W.-C.; Li, C.-H.; Chan, Y.-H.; Kuo, L.-T.; Mao, C.-T.; Shyu, Y.-C.; et al. Artificial Intelligence-Assisted Identification of Genetic Factors Predisposing High-Risk Individuals to Asymptomatic Heart Failure. Cells 2021, 10, 2430. https://0-doi-org.brum.beds.ac.uk/10.3390/cells10092430

AMA Style

Yang N-I, Yeh C-H, Tsai T-H, Chou Y-J, Hsu PW-C, Li C-H, Chan Y-H, Kuo L-T, Mao C-T, Shyu Y-C, et al. Artificial Intelligence-Assisted Identification of Genetic Factors Predisposing High-Risk Individuals to Asymptomatic Heart Failure. Cells. 2021; 10(9):2430. https://0-doi-org.brum.beds.ac.uk/10.3390/cells10092430

Chicago/Turabian Style

Yang, Ning-I, Chi-Hsiao Yeh, Tsung-Hsien Tsai, Yi-Ju Chou, Paul Wei-Che Hsu, Chun-Hsien Li, Yun-Hsuan Chan, Li-Tang Kuo, Chun-Tai Mao, Yu-Chiau Shyu, and et al. 2021. "Artificial Intelligence-Assisted Identification of Genetic Factors Predisposing High-Risk Individuals to Asymptomatic Heart Failure" Cells 10, no. 9: 2430. https://0-doi-org.brum.beds.ac.uk/10.3390/cells10092430

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