Original Research ARTICLE
In Vivo and In Vitro Analysis in Coronary Artery Disease Related to Type 2 Diabetes
- IRCCS SDN, Naples, Italy
Aim: The leading cause of morbidity and mortality in patients with type 2 diabetes mellitus (DM) is coronary artery disease (CAD), a condition often asymptomatic but severe in these patients. Although glucose metabolism impairment and oxidative stress are known actors in the endothelial dysfunction/remodeling that occurs in diabetic patients, the relationship between cardiovascular disorders and DM is not fully understood. We have performed both an in vivo imaging and in vitro molecular analysis to investigate diabetic-specific CAD alterations.
Methods: Computed tomography coronary angiography (CTCA) was performed in a group of 20 diabetic patients with CAD (DM+CAD+), 20 non-diabetic with CAD (DM−CAD+), 10 diabetic non-CAD patients (DM+CAD−), and 20 non-diabetic healthy subjects (HS). Imaging quantitative parameters such as calcium score (Cascore), calcified plaque volume (CPV), non-calcified plaque volume (NCPV), total plaque volume (TPV), remodeling index (RI), and plaque burden were extracted for each CAD subject. Moreover, the expression levels of superoxide dismutase 2 (SOD2) and liver X receptor alpha (LXRα) genes were analyzed in the peripheral blood mononuclear cells, whereas hyaluronan (HA) concentrations were evaluated in the plasma of each subject.
Results: Imaging parameters, such as Cascore, CPV, RI, and plaque burden, were significantly higher in DM+CAD+ group, compared to DM−CAD+ (P = 0.019; P = 0.014; P < 0.001, P < 0.001, respectively). SOD2 mRNA was downregulated, while LXRα gene expression was upregulated in DM+CAD−, DM+CAD+, and DM−CAD+ groups compared to HS (P = 0.001, P = 0.03, and P = 0.001 for SOD2 and P = 0.006, P = 0.008, and P < 0.001 for LXRα, respectively). Plasmatic levels of HA were higher in DM−CAD+, DM+CAD−, and DM+CAD+ groups, compared to HS (P = 0.001 for the three groups). When compared to DM−CAD+, HA concentration was higher in DM+CAD− (P = 0.008) and DM+CAD+ (P < 0.001) with a significant difference between the two diabetic groups (P = 0.003). Moreover, HA showed a significant association with diabetes (P = 0.01) in the study population, and the correlation between HA levels and glycemia was statistically significant (ρ = 0.73, P < 0.001).
Conclusion: In our population, imaging parameters highlight a greater severity of CAD in diabetic patients. Among molecular parameters, HA is modulated by diabetic CAD-related alterations while SOD2 and LXRα are found to be more associated with CAD but do not discriminate between diabetic and non-diabetic subgroups.
Type 2 diabetes mellitus (DM) is the most important risk factor for the onset of coronary artery disease (CAD), causing glucose metabolism impairment and endothelial dysfunction mediated by oxidative stress and inflammation (1). A complex network of signaling pathways is involved in these pathological processes leading to the development and progression of cardiac dysfunction. In response to myocardial damage, the heart undergoes a progressive anatomical and functional transformation known as “remodeling” (2).
Several imaging modalities have been used to detect CAD in diabetic patients including invasive coronary angiography, myocardial scintigraphy and dobutamine stress echocardiography (3). Even if invasive coronary angiography is the gold standard for identifying obstructive lesions, it only depicts the lumen of the vessel, greatly underestimating the burden of atherosclerosis (4). Myocardial scintigraphy and dobutamine stress echocardiography highlight perfusion defects (inducible ischemia and necrosis), but they lack a direct visualization of coronary arteries (5). Unlike these, computed tomography coronary angiography (CTCA) is a powerful diagnostic tool to rule out CAD thanks to its high negative predictive value (6). It allows quantification of atherosclerotic burden providing comprehensive information about the location, severity, and features of coronary atherosclerotic plaques and can be useful for risk stratification (7, 8).
Atherosclerosis is a multistage pathological condition involving an imbalanced lipid metabolism and immune response leading to a chronic inflammation of the arterial wall with the formation of the atherosclerotic plaque and consequent thickening of vessel wall and lumen stenosis (9, 10).
The first step of atherosclerosis is endothelial dysfunction; atherosclerotic lesions initiate in regions characterized by low shear stress resulting in an increase of adhesiveness of circulating monocytes to the vessel wall and subendothelial accumulation of low-density lipoprotein (LDL) (9). Common cardiovascular risk factors, such as smoking, diabetes, hypertension, and hypercholesterolemia, are causes of dysfunction endothelial (10). The LDL particles in the intima are susceptible to oxidation by reactive oxygen species or other enzymes released from inflammatory cells. Oxidized LDL triggers the expression of adhesion molecules and the secretion of chemokines by endothelial cells that drive the intimal infiltration by immune cells forming the so-called “fatty streaks” especially consisting of monocyte-derived macrophage-like foam cells. Subsequently, vascular smooth muscle cells migrate and proliferate into the site of lesion producing an excessive amount of connective tissue with the consequent formation of the fibroatheromatous plaque leading to thickening of vessel wall and stenosis of coronary lumen (9–11). One of the major issues in CAD diagnosis and management is that symptoms onset in the advanced state of disease. Indeed, most individuals show no manifestations for long time before the first onset of symptoms, often with a fatal event.
Oxidative stress is a key component in the development and progression of DM and its vascular complications such as CAD (12, 13). The onset and progression of CAD involves multiple cell types, and whole-blood gene expression profiling has the potential to provide information about dynamic changes in disease states and on underlying disease mechanisms (14).
Superoxide dismutase 2 (SOD2) is one of the major antioxidant defense systems against free radicals (15). Mutations or polymorphisms of SOD2 gene are associated with DM progression and complications, where the reduction of total antioxidant capacity and depletion of plasma antioxidants could be related to induced-oxidative stress damage (16–21).
Nuclear liver X receptors (LXR) comprise liver X receptor alpha (LXRα) and liver X receptor beta (LXRβ), which are key regulators of macrophage function, controlling transcriptional programs involved in lipid homeostasis and inflammation. The inducible LXRα is highly expressed in macrophages, liver, adrenal gland, intestine, adipose tissue, lung, and kidney, whereas LXRβ is ubiquitously expressed (22). LXRs are involved in the regulation of cholesterol metabolism fundamental for the pathogenesis of CAD and inhibit atherogenesis, inflammation and autoimmune reactions (22). Furthermore, an additional role of LXRs is to contribute to glucose homeostasis, demonstrating potent glucose-lowering and insulin-sensitizing effects (23, 24). Despite extensive research in the field of LXR biology, however, very little is known about the regulation of expression and activity of these receptors.
Hyaluronan (HA) is present in low amount in normal blood vessels but increases in vascular diseases as well as in DM (25). It seems to have an important role in diabetic angiopathy (26–28) and is associated with an increased risk for developing CAD also in non-diabetic patients (29). HA is increased in vascular plaques, and its high metabolism causes their destabilization (30). Furthermore, the fragmentation of HA triggers inflammatory processes and activates leukocytes to produce superoxide radical causing oxidative stress (31).
To date, studies integrating parameters calculated by CTCA and biological markers in DM patients have been focused on the association between CTCA findings (mostly coronary artery calcium) and biological markers of inflammation (IL-6, IL-1β, TNF-α, hs-CRP, and YKL-40) and endothelial dysfunction (sVCAM-1, sICAM-1, and sICAM-3) (32–36). There are no data about the association between imaging parameters and gene expression profiling in DM.
In this study, we have analyzed the three above mentioned molecular markers that underlie important steps of the atherosclerotic process: endothelial dysfunction, oxidative stress, lipid homeostasis, and inflammation. In this regard, we have analyzed SOD2 and LXRα gene expression and HA plasmatic concentrations in a group of 20 diabetic patients with known CAD (DM+CAD+), 20 non-diabetic patients with CAD (DM−CAD+), 10 diabetic non-CAD patients (DM+CAD−), and 20 non-diabetic healthy subjects (HS). Furthermore, the purpose of our study was to investigate diabetic-specific CAD alterations using both quantitative imaging parameters derived from CTCA and molecular biomarkers.
Materials and Methods
Computed tomography coronary angiography was performed in 20 DM+CAD+ patients, 20 DM−CAD+ patients, 10 DM+CAD− patients, and 20 HS referred to our institution for suspected CAD. All clinical characteristics such as laboratory parameters, presence of cardiovascular risk factors, and medical history were accurately recorded.
Diabetes was defined as treatment with drugs or fasting blood glucose ≥126 mg/dL. Dyslipidemia was defined as treatment with drugs or fasting serum total cholesterol ≥240 mg/dL, or LDL cholesterol ≥140 mg/dL, or high-density lipoprotein cholesterol <40 mg/dL, or triglyceride ≥150 mg/dL. Hypertension was defined as treatment with drugs or systolic blood pressure (SBP) ≥140 mmHg or diastolic blood pressure (DBP)≥ 90 mmHg. Anthropometrical measurements including body weight and height were recorded and body mass index (BMI) was calculated. Blood pressure and resting heart rate were measured after ≥5 min rest with a sphygmomanometer. Physical activity in HS and patients was evaluated according to the WHO guidelines for adults 18–65 aged; specifically the performance of at least 150 min of moderate-intensity aerobical physical activity per week [50–70% of maximum heart rate (MHR)] or at least 75 min of vigorous-intensity aerobic physical activity throughout the week (70–80% MHR) (37). None of the recruited subjects had physical disabilities.
Patients with known history of cancer, cardiomyopathy, active infections, chronic or immune-mediated diseases, renal failure, hepatic failure, and not suitable for cardiac imaging (atrial fibrillation, arrhythmia, or pre-scan heart rate greater than 65 bpm) were excluded from the study to avoid confounding effects due to other variables.
Peripheral venous blood samples were collected after a 12 h overnight fasting immediately before i.v. cannulation for CTCA examination. All tubes were centrifuged within 30 min of collection at 1,900 g for 10 min at 4°C to separate plasma and cellular components. Aliquots of plasma were transferred into cryostat tubes and stored at −80°C until analysis. PBMNCs were isolated by Ficoll gradient using HISTOPAQUE-1077 (Sigma Diagnostics, MO, USA) and frozen at −80°C until total RNA extraction. All biological samples were stored at the IRCCS SDN Biobank (38). The study and the protocol were approved and reviewed by the institutional ethics committee (IRCCS Fondazione SDN, protocol no. 7-13). The study was performed in accordance with the ethical standards of the institutional ethics committee and with the Helsinki Declaration. A written informed consent was obtained from all subjects enrolled.
CT Angiography Protocol and Image Analysis
Computed tomography coronary angiographies were performed on a CT scanner (Discovery CT750 HD, GE Healthcare), with a 64 mm × 0.625 mm collimation, 350 ms rotation time, and 228 ms temporal resolution. A prospective ECG-triggered scan without contrast medium was used for calcium score (Cascore) evaluation followed by a retrospective scan with ECG tube current modulation. Contrast enhancement was obtained by a bolus tracking technique with scan starting when a region of interest placed in the ascending aorta at the pulmonary bifurcation reached a threshold of 150 Hounsfield unit (HU). Contrast material (iomeprol 400 mg I/mL, Iomeron 400, Bracco, Milan, Italy) was injected at 5–6 mL/s through an 18-gauge intravenous antecubital catheter and was followed by saline solution at the same flow. Tube voltage and contrast agent volume were adapted to patient anatomical features such as BMI, calcifications, or stents. Images were reconstructed with a section thickness of 0.625 mm and an increment of 0.4 mm; standard and sharp reconstruction filter kernels were used; an additional sharper convolution kernel was used in patients with stents or calcification. The best data set was chosen according to the phase of the cardiac cycle with lower artifacts and coronary motions. Images were sent to a dedicated offline workstation (GE Advantage workstation 4.6, GE Healthcare) where MIP, cMPR, and 3D volume rendering were generated. Cascore was calculated by the SmartScore tool to obtain the Agatston score. Total plaque volume (TPV), non-calcified plaque volume (NCPV), calcified plaque volume (CPV), and total lumen volume were measured for the major coronaries using the HU cutoff values reported in Ref. (39). The resulting values were summed to determine a per-patient plaque volume. Total vessel volume was determined summing TPV and total lumen volume. Plaque burden was obtained dividing TPV by total vessel volume (40). The remodeling index (RI) was calculated by dividing the cross-sectional vessel area at the site of maximum luminal narrowing including plaque by the cross-sectional vessel area in the most proximal atherosclerotic free segment chosen as reference (41). The total number of coronary artery segments exhibiting plaque (NSP) was determined according to the modified American Heart Association 16-segment classification (42) for each patient (less or more than 8 segments affected). Significant coronary stenosis was defined as a decrease in the luminal diameter of >50% in one or more of the major coronary arteries; the total number of coronaries (NCS) with significant stenosis was calculated for each patient (less or more than one stenotic vessel). All scans were analyzed by two experienced, independent radiologists; therefore, a consensus interpretation was arrived to obtain a final coronary CT diagnosis according to the international SCCT guidelines (43).
RNA Extraction and Reverse Transcription
Total RNA was isolated from PBMCs using TRIzol Reagent (Thermo Fischer Scientific, USA) as previously described (44). The quantity and quality of RNA were measured using the NanoDrop 1000 (Thermo Fischer Scientific, USA). Total RNA (0.5 µg) was reversed transcribed (RT) using the SuperScript® III First-Strand Synthesis SuperMix for qRT-PCR (Thermo Fischer Scientific, USA) according to the protocol of the manufacturer. The RT was performed using the Bio-Rad iCycler Thermal Cycler with the following protocol: incubation at 25°C for 10 min (primer annealing), 42°C for 30 min (cDNA synthesis), and 85°C for 5 min (termination of cDNA synthesis). Immediately after, the samples were cooled down and stored at −20°C.
Quantitative Real-time PCR
The optimal reference genes for the study were selected as previously reported (45). Gene expression was quantified on the MyiQ™ Single-Color Real-Time PCR Detection System (Bio-Rad Laboratories, USA). Primers pairs were designed through OLIGO 6.7 program, and their specificity was verified with the BLAST program for test of sequence homology, a test for secondary structures and optimization of multiplex setup. All primers were purchased from Life Technologies. All samples were run in triplicate for genes of interest and reference genes using 1 µL of cDNA and iQ™ SYBR® Green Supermix (Bio-Rad Laboratories, USA) in a 25 µL final volume reaction. The thermal profile employed was 3 min of initial step of denaturation at 95°C followed with denaturation for 15 s at 95°C, annealing at 60°C for 30 s, and elongation at 72°C for 30 s for 40 cycles. Melt curve analysis was performed to verify a single product species. Relative expression (fold change) was calculated by the 2−ΔΔCT method (46). Mean and SE were determined by averaging relative expression levels across three independent experiments, each determined in triplicate.
Plasmatic levels of HA were determined by enzyme-linked immunosorbent assay (ELISA) using Quantikine Hyaluronan Immunoassay kit (DHYAL0) (R&D Systems, Abingdon, UK), in accordance with the protocol supplied by the manufacturer. Briefly, samples were incubated with HA binding protein coated on microplates for 2 h at room temperature. After incubation, the microplates were washed five times with wash buffer, and further incubated with 100 µL of peroxidase labeled HA binding protein for 2 h at room temperature. After incubation, the microplates were again washed five times, and further incubated with 100 µL of peroxidase substrate for 30 min at room temperature in a dark room. The reaction was stopped by the addition of 100 µL of stop solution. The optical density of each well was determined using a microplate reader set to 450 nm within 30 min. HA concentration in each sample was calculated using the standard curve obtained with the purified HA solutions, included in the kit as references.
Statistical analysis was performed using R Core Team (version 3.03 Austria, Vienna). Continuous variables were expressed as mean ± SD or as median (1 quartile and 3 quartile). Data were tested for normality through the Shapiro–Wilk test and for homoscedasticity through the Levene test. For comparison between two groups, t-test was used if gaussianity was met; otherwise the Mann–Whitney U test was chosen. For comparison among four groups, the one-way analysis of variance was used if both gaussianity and homoscedasticity were met; otherwise the Kruskal–Wallis test was chosen. In case of statistical significance, the Tukey–Kramer test and the Conover test were used for multiple comparisons as parametric and non-parametric test, respectively. Categorical variables were expressed as percentage and were compared using the Fisher’s exact test. The Spearman correlation test was performed to assess linear relationship between variables; in case of binary variables, the association was tested by the Wilcoxon rank sum test. A P < 0.05 was considered for statistical significance (rounded to the third decimal place).
Clinical Characteristics of Study Groups
The baseline demographic and clinical characteristics of the study population are summarized in Table 1. Heart rate was significantly different between HS and DM+CAD+ (P < 0.01) and HS and DM−CAD+ (P < 0.01) since only 10% of HS was in treatment with beta blocker agents, while no statistical significance was found between both CAD groups and DM+CAD−. Considering the metabolic markers, glycemia was significantly higher in DM+CAD+ and DM+CAD− patients compared to HS (P < 0.01 and P < 0.001, respectively) and DM−CAD+ subjects (P < 0.01 and P < 0.001, respectively). Of diabetic patients, in DM+CAD+ group, 16% were insulin users, 64% were in treatment with antihyperglycemic agents, and 20% were not in treatment; in DM+CAD− group, the percentage of treatments were, respectively, 10% for insulin, 80% for antihyperglycemic drugs, and 10% were not treated. Total cholesterol, LDL- and HDL-cholesterol plasmatic concentrations did not significantly differ among the four groups, reflecting the effects of statin therapy to which 73.68% of DM+CAD+, 50% of DM−CAD+, 40% of DM+CAD−, and 5% of HS were subjected. Furthermore, SBP and DBP were not statistical different among the groups (P = 0.50 and P = 0.52, respectively). In this regard, hypertensive patients were in treatment with beta blocker agents (P = 0.008), calcium channel blockers (P = 0.68), and ACE inhibitors (P = 0.06).
There was significant difference between DM−CAD+ and DM+CAD+ according to NCS and NSP (P = 0.026, P = 0.04, respectively). Cascore was significantly higher in DM+CAD+ compared to DM−CAD+ (Figures 1 and 2): 1,068.7 (517.2–2,086.85) vs 214.05 (72.98–970.15) P = 0.019. As regards plaque characterization, CPV was significantly higher in DM+CAD+ [105.85 (51.2–341.73) mm3] compared to DM−CAD+ [42 (7.2–105.9) mm3] P = 0.014, but there was no significant difference according to NCPV and TPV: 519.85 (411.93–1,064.85) mm3 for DM+CAD+ and 421.85 (240.10–689.58) mm3 for DM−CAD+ P = 0.37 and 688.95 (470.05–1,436) mm3 for DM+CAD+ vs 454.45 (257.78–820.83) mm3 for DM−CAD+ P = 0.16, respectively. RI was 1.40 ± 0.24 for DM+CAD+ and 1 ± 0.19 for DM−CAD+ P < 0.001, and plaque burden was 0.45 ± 0.14 for DM+CAD+ and 0.27 ± 0.15 for DM−CAD+ P < 0.001. Results are summarized in Table 2.
Figure 1. (A,D) Non-contrast enhanced images showing calcium deposits (yellow) on the left descending coronary artery (LAD) in a non-diabetic CAD patient (DM−CAD+) and in a diabetic CAD patient (DM+CAD+), respectively. (B,E) cMPR of LAD is provided for DM−CAD+ and DM+CAD+. (C,F) Plaque characterization: the calcific (yellow) and non-calcific (pink) components of the plaque are highlighted; the vessel lumen is represented in green. DM+CAD+ displayed significantly higher coronary calcium values compared to DM−CAD+.
Figure 2. (A,B) Cross-sectional view and cMPR of the left descending coronary artery (LAD) in a non-diabetic CAD patient (DM−CAD+) and in a diabetic CAD patient (DM+CAD+). In cross-sectional images, the vessel lumen is represented in green whereas the calcific component of the plaque is red.
In our population, RI highly correlated with plaque burden (ρ = 0.65, P < 0.001). Cascore showed a positive correlation with NCPV (ρ = 0.83, P < 0.001), CPV (ρ = 0.96, P < 0.001), TPV (ρ = 0.88, P < 0.001), and plaque burden (ρ = 0.60, P < 0.001). Moreover, a significant correlation was found between plaque burden and NCPV (ρ = 0.57, P < 0.001), CPV (ρ = 0.64, P < 0.001), and TPV (ρ = 0.60, P = 0.001).
Gene Expression Profiling
We evaluated, by quantitative real-time PCR, SOD2 and LXRα gene expression in PBMCs from our population (Table 3). For both genes, ΔCT was computed and compared between the four groups. Molecular analysis showed that SOD2 mRNA was downregulated in DM+CAD− (ΔCT = 5.70 ± 3.28; fold change = 0.10 ± 0.03; P = 0.001), DM+CAD+ (ΔCT = 4.57 ± 3.56; fold change = 0.22 ± 0.08; P = 0.03), and DM−CAD+ (ΔCT = 5.75 ± 3.04; fold change = 0.10 ± 0.05; P = 0.001) compared to HS (ΔCT = 2.36 ± 2.61), with no statistically significant difference between the two CAD groups (with and without DM) (Figure 3A). LXRα gene expression was significantly upregulated in DM+CAD− (ΔCT = 2.77 ± 1.36; fold change = 4.51 ± 0.20 P = 0.006), DM+CAD+ (ΔCT = 3.27 ± 1.79; fold change = 3.19 ± 0.42; P = 0.008), and DM−CAD+ (ΔCT = 2.15 ± 1.16; fold change = 6.93 ± 0.70; P < 0.001) compared to the HS (ΔCT = 4.94 ± 2.14), with a significant difference between the two CAD groups (P = 0.03) (Figure 3B). No statistically significant correlation was found between SOD2 and Cascore (ρ = −0.04, P = 0.81), NCPV (ρ = −0.03, P = 0.85), and TPV (ρ = −0.02, P = 0.89) as well as between LXRα and Cascore (ρ = 0.13, P = 0.44), NCPV (ρ = 0.10, P = 0.55), and TPV (ρ = 0.10, P = 0.57).
Figure 3. (A) Superoxide dismutase 2 (SOD2) mRNA relative expression in PBMNCs of healthy subjects (HS), diabetic non-CAD patients (DM+CAD−), non-diabetic CAD patients (DM−CAD+), and diabetic CAD patients (DM+CAD+). (B) Liver X receptor alpha (LXRα) gene expression in PBMNCs of HS, DM+CAD−, DM−CAD+, and DM+CAD+. (C) Median plasma hyaluronan (HA) concentrations (ng/mL) in HS, DM+CAD−, DM−CAD+, and DM+CAD+ (*P < 0.05; **P < 0.01; ***P < 0.001). When not specified, significance is referred to the comparison vs HS. (D) Correlation between HA levels and glycemia in HS, DM+CAD−, DM−CAD+, and DM+CAD+ patients (ρ = 0.73, P < 0.001).
Comparison of HA Levels
In the HS group, mean concentration of plasma HA was 46.90 ± 23.79 ng/mL. Compared with HS, HA concentrations were higher in DM−CAD+ (90.05 ± 35.11 ng/mL; P = 0.001), DM+CAD− (105.56 ± 18.13 ng/mL; P = 0.001), and DM+CAD+ (120.74 ± 21.17 ng/mL; P = 0.001) (Figure 3C). When compared to DM−CAD+, HA concentration was significantly higher in DM+CAD− (P = 0.008) and DM+CAD+ (P < 0.001) with a significant difference between two diabetic groups (P = 0.003). Correlation of HA levels with Cascore, NCPV, and TPV revealed ρ = 0.29, P = 0.076, ρ = 0.30, P = 0.073, and ρ = 0.31, P = 0.063, respectively.
Risk Factors and Molecular Markers Analysis
Analysis on risk factors and molecular data showed no significant association between sex and SOD2, LXRα, and HA (P = 0.85; P = 0.21; P = 0.75, respectively) as well as regarding familiarity (P = 0.83; P = 0.64; P = 0.55, respectively) and smoke (P = 0.73; P = 0.17; P = 0.49, respectively). Furthermore, no significant correlations were found between BMI and the three molecular parameters (SOD2 ρ = 0.12, P = 0.35; LXRα ρ = −0.11, P = 0.35; HA ρ = 0.19, P = 0.15). HA was significantly correlated with age (ρ = 0.46, P < 0.001), unlike SOD2 (ρ = 0.03, P = 0.97) and LXRα (ρ = 0.03, P = 0.95). HA showed also a significant association with dyslipidemia (P = 0.01) and diabetes (P = 0.01) in the study population, while statistical analysis on hypertension revealed a P = 0.08. On the other hand, SOD2 and LXRα were not significantly associated with the previously mentioned risk factors (SOD2 vs dyslipidemia P = 0.45, SOD2 vs hypertension P = 0.63, SOD2 vs diabetes P = 0.47, LXRα vs dyslipidemia P = 0.85, LXRα vs hypertension P = 0.53, and LXRα vs diabetes P = 0.61). Correlation between HA levels and glycemia was statistically significant (ρ = 0.73, P < 0.001) (Figure 3D), while no significant correlation was found between SOD2 (P = 0.62) and LXRα (P = 0.55) gene expression and glycemia.
In this study, we have exploited an imaging and molecular based analysis to investigate diabetic-specific CAD alterations in selected groups of patients.
Calcium score, CPV, plaque burden, and RI were significantly higher in DM+CAD+ compared to DM−CAD+. Previous studied have examined CAD and plaque features in diabetic patients by CTCA. Diabetics showed extensive coronary artery calcium deposits and, therefore, a larger atherosclerotic plaque burden with a consequent higher risk for all-cause mortality than in non-diabetic patients (5, 47–54). Gao et al. (47) found that diabetics compared to non-diabetics have higher total coronary artery calcium, a higher proportion of coronary segments with plaque and multivessel obstructive disease. In a study by Van Werkhoven et al., obstructive CAD and the number of diseased segments, with obstructive and non-obstructive plaques, were higher in diabetics than non-diabetics. Total Agatston score was higher in diabetic patients (440 ± 786 vs 195 ± 404, P < 0.001) (5). Khazai et al. found that segment involvement score, segment stenosis score, and total plaque score were higher in diabetics but there was no significant difference in the number of non-calcified plaque between the two groups (50). In one study by Pundziute et al., diabetics showed more diseased segments and more segments with non-obstructive CAD, but Agatston score was similar between the two groups (54). Furthermore, Chu et al. detected more calcified plaques than mixed or non-calcified plaques in diabetics. Among the different degrees of stenosis, mild narrowing was most common, and no significant difference between non-obstructive stenosis and obstructive stenosis was observed (48). In agreement with the aforementioned works, in our study, DM+CAD+ presented more diseased coronaries in terms of coronary calcium, significant stenoses, atherosclerotic burden, and extent of disease. Furthermore, we have quantified RI in diabetic patients by CTCA providing an additional prognostic value comparable only to invasive procedures such as intravascular ultrasound (55, 56). A recent study analyzed CAD features comparing hypertensive, dyslipidemic, and diabetic patients by CTCA reporting the prevalence of positive remodeling as a qualitative parameter (57).
Positive coronary arterial remodeling is a compensatory enlargement of coronary arterial lumen in response to atherosclerotic plaque formation. Histopathological studies proved that positive remodeling is associated with infiltration of inflammatory cells, expression of pro-inflammatory cytokines, and increased protease activity (58, 59). Positive remodeling is associated with vulnerable plaque and progression of atherosclerosis. High plaque burden, together with positive remodeling, means more prone to rupture plaques in diabetic patients and, therefore, a worse prognosis and a major likelihood of cardiac event occurrence.
In the last decade, a great amount of data demonstrated a complex interaction between blood cells and the arterial wall with the consequent activation of oxidative and inflammatory pathways, leading to the development of CAD.
Our results showed that the expression levels of SOD2 gene were reduced in CAD patients compared to HS, while no significant difference was found between diabetic and non-diabetic CAD subjects. Previous studies reported controversial findings for the effect of SOD2 activity relative to CAD. A recent study by Peng et al. (60) showed that plasmatic concentration of SOD1 and SOD2 was higher in CAD than in healthy control. Our findings were in line with a gene expression study performed by Abdullah et al. (61) showing a downregulation of this gene in PBMCs of angiographically confirmed CAD patients (≥50% stenosis). These data indicate that SOD2 might serve as surrogate biomarker for CAD.
Data from in vitro and in vivo models have demonstrated a key role of LXRα in the regulation of processes involved in CAD and DM such as inflammation and glucose homeostasis (62, 63). Our findings reported that LXRα gene expression was significantly upregulated in DM+CAD+ and DM−CAD+ compared to HS. Although previous study by Dahlman et al. (64) investigated the association of LXRα and DM, we demonstrated also a differential expression of this gene between DM+CAD+ and DM−CAD+ groups suggesting this parameter as a possible biological hallmark for diabetic condition. HA plasmatic concentrations showed significant difference between diabetic and non-diabetic patients with higher values in patients affected by both DM and CAD suggesting a possible additive detrimental effect on endothelial dysfunction. A significant positive correlation was found between HA levels and glycemia in our study population. Our findings were in line with previous studies, also reporting a critical role for HA in DM-related atherosclerosis (26–29, 65). In vascular dysfunction, HA triggers smooth muscle cells’ dedifferentiation, which contributes to vessel wall thickening. Furthermore, HA is able to modulate inflammation by altering the adhesive properties of endothelial cells. In hyperglycemic conditions, HA accumulates in vessels and can contribute to the diabetic complications in macro- and microvasculature (25).
However, no study has yet examined the relationship between HA levels and vascular function assessed by CTCA. Our data suggested that serum HA levels positively correlated with poor glycemic control and angiopathy and, due to the pivotal role in favoring atherogenesis, this molecule could be used as a surrogate marker of vascular function.
In vitro molecular analysis represents a promising tool to stratify patients for CAD risk, while in vivo CTCA analysis is able to identify and characterize selective diabetic coronary features. These results gain clinical relevance, considering that most patients referred to elective invasive coronary angiography for CAD suspicious are not found to have obstructive CAD (66, 67). In patients with molecular alterations suggestive for CAD, we demonstrated by CTCA specific changes of coronary plaques in diabetic patients. Moreover, recently, this imaging technique has been used to evaluate its long-term prognostic value among patients with diabetes mellitus compared with non-diabetic subjects (68).
Nevertheless, our study has some limitations: the reduced sample size has influenced the statistical power; therapeutic treatments could have affected our results; a more accurate analysis with different genomic/proteomic techniques, on a wider pool of in vitro markers is needed to deeply investigate molecular and imaging phenotypic interplay in diabetic CAD patients. The analyzed biomarkers are not myocardial specific CAD molecules but can be downregulated or upregulated in blood also in presence of atherosclerotic processes involving peripheral arteries and/or supra aortic vessels. Moreover, recent studies have demonstrated that diabetes can be considered a CAD equivalent condition, independently from the clinical/imaging evidences of pathology (69, 70), determining the choice of specific diabetic-related CAD biomarkers attractive.
This study suggests an imaging and molecular based analysis to investigate cardiovascular alterations in diabetic patients. CTCA imaging parameters highlight a greater severity of CAD in diabetic patients. Among molecular parameters, HA is modulated by diabetic CAD-related alterations while SOD2 and LXRα are found to be more associated with CAD rather than to diabetes. Further studies are needed to better characterize the pathology and identify more specific biomarkers, also considering the complex multifactorial pathophysiology of CAD in diabetic patients.
The study was approved and performed in accordance with the ethical standards of the institutional ethics committee (IRCCS Fondazione SDN, protocol no. 7-13) and with the Helsinki Declaration. Written informed consent was obtained from all subjects for being included in the study.
CC designed and supervised the study. TI and EF recruited subjects, performed the experiments and data analysis, and wrote the manuscript. MA, MS, and CC reviewed the manuscript. All the authors read and approved the final manuscript and agreed to its submission.
Conflict of Interest Statement
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.
This research was supported by Italian Ministry of Health, “progetto Giovani Ricercatori 2011–2012” (Project code: GR-2011-02349436).
DM, diabetes mellitus; CAD, coronary artery disease; HS, healthy subjects; DM+CAD−, diabetic non-CAD patients; DM+CAD+, diabetic patients with CAD; DM−CAD+, non-diabetic patients with CAD; CTCA, computed tomography coronary angiography; HU, Hounsfield unit; Cascore, calcium score; CPV, calcified plaque volume; NCPV, non-calcified plaque volume; TPV, total plaque volume; RI, remodeling index; SOD2, superoxide dismutase 2; LXRα, liver X receptor alpha; HA, hyaluronan; NSP, number of coronary artery segments with plaque; NCS, number of coronaries with significant stenosis.
2. Schoenhagen P, Ziada KM, Vince DG, Nissen SE, Tuzcu EM. Arterial remodeling and coronary artery disease: the concept of “dilated” versus “obstructive” coronary atherosclerosis. J Am Coll Cardiol (2001) 38:297–306. doi:10.1016/S0735-1097(01)01374-2
3. Chopra S, Peter S. Screening for coronary artery disease in patients with type 2 diabetes mellitus: an evidence-based review. Indian J Endocrinol Metab (2012) 16:94–101. doi:10.4103/2230-8210.91202
4. Forte E, Aiello M, Inglese M, Infante T, Soricelli A, Tedeschi C, et al. Coronary artery aneurysms detected by computed tomography coronary angiography. Eur Heart J Cardiovasc Imaging (2016). doi:10.1093/ehjci/jew218
5. Van Werkhoven JM, Cademartiri F, Seitun S, Maffei E, Palumbo A, Martini C, et al. Diabetes: prognostic value of CT coronary angiography – comparison with a nondiabetic population. Radiology (2010) 256:83–92. doi:10.1148/radiol.1090600
6. Forte E, Inglese M, Infante T, Schiano C, Napoli C, Soricelli A, et al. Anomalous left main coronary artery detected by CT angiography. Surg Radiol Anat (2016) 38:987–90. doi:10.1007/s00276-016-1634-9
8. Rizvi A, Hartaigh BÓ, Danad I, Han D, Lee JH, Gransar H, et al. Diffuse coronary artery disease among other atherosclerotic plaque characteristics by coronary computed tomography angiography for predicting coronary vessel-specific ischemia by fractional flow reserve. Atherosclerosis (2017) 258:145–51. doi:10.1016/j.atherosclerosis.2017.01.018
12. Folli F, Corradi D, Fanti P, Davalli A, Paez A, Giaccari A, et al. The role of oxidative stress in the pathogenesis of type 2 diabetes mellitus micro- and macrovascular complications: avenues for a mechanistic-based therapeutic approach. Curr Diabetes Rev (2011) 7:313–24. doi:10.2174/157339911797415585
17. Ascencio-Montiel Ide J, Parra EJ, Valladares-Salgado A, Gómez-Zamudio JH, Kumate-Rodriguez J, Escobedo-de-la-Peña J, et al. SOD2 gene Val16Ala polymorphism is associated with macroalbuminuria in Mexican type 2 diabetes patients: a comparative study and meta-analysis. BMC Med Genet (2013) 14:110. doi:10.1186/1471-2350-14-110
19. Katakami N, Kaneto H, Matsuoka TA, Takahara M, Osonoi T, Saitou M, et al. Accumulation of oxidative stress-related gene polymorphisms and the risk of coronary heart disease events in patients with type 2 diabetes-an 8-year prospective study. Atherosclerosis (2014) 235:408–14. doi:10.1016/j.atherosclerosis.2014.05.936
20. Vats P, Sagar N, Singh TP, Banerjee M. Association of superoxide dismutases (SOD1 and SOD2) and glutathione peroxidase 1 (GPx1) gene polymorphisms with type 2 diabetes mellitus. Free Radic Res (2015) 49:17–24. doi:10.3109/10715762.2014.971782
21. Pourvali K, Abbasi M, Mottaghi A. Role of superoxide dismutase 2 gene Ala16Val polymorphism and total antioxidant capacity in diabetes and its complications. Avicenna J Med Biotechnol (2016) 8:48–56.
23. Liu X, Li G, Zhu H, Huang L, Liu Y, Ma C, et al. Beneficial effect of berberine on hepatic insulin resistance in diabetic hamsters possibly involves in SREBPs, LXRα and PPARα transcriptional programs. Endocr J (2010) 57:881–93. doi:10.1507/endocrj.K10E-043
24. Liu Y, Yan C, Wang Y, Nakagawa Y, Nerio N, Anghel A, et al. Liver X receptor agonist T0901317 inhibition of glucocorticoid receptor expression in hepatocytes may contribute to the amelioration of diabetic syndrome in db/db mice. Endocrinology (2006) 147:5061–8. doi:10.1210/en.2006-0243
25. Moretto P, Karousou E, Viola M, Caon I, D’Angelo ML, De Luca G, et al. Regulation of hyaluronan synthesis in vascular diseases and diabetes. J Diabetes Res (2015) 2015:167283. doi:10.1155/2015/167283
26. Heickendorff L, Ledet T, Rasmussen LM. Glycosaminoglycans in the human aorta in diabetes mellitus: a study of tunica media from areas with and without atherosclerotic plaque. Diabetologia (1994) 37:286–92. doi:10.1007/BF00398056
28. Morita M, Yano S, Ishibashi Y, Nakata N, Kurioka S, Sugimoto T. Close relationship between serum hyaluronan levels and vascular function in patients with type 2 diabetes. Biomarkers (2014) 19:493–7. doi:10.3109/1354750X.2014.940502
29. Papanastasopoulou C, Papastamataki M, Karampatsis P, Anagnostopoulou E, Papassotiriou I, Sitaras N. Cardiovascular risk and serum hyaluronic acid: a preliminary study in a healthy population of low/intermediate risk. J Clin Lab Anal (2017) 31(1). doi:10.1002/jcla.22010
30. Bot PT, Pasterkamp G, Goumans MJ, Strijder C, Moll FL, de Vries JP, et al. Hyaluronic acid metabolism is increased in unstable plaques. Eur J Clin Invest (2010) 40:818–27. doi:10.1111/j.1365-2362.2010.02326.x
32. Dayan A, Narin B, Biteker M, Aksoy S, Fotbolcu H, Duman D. Coronary calcium score, albuminuria and inflammatory markers in type 2 diabetic patients: associations and prognostic implications. Diabetes Res Clin Pract (2012) 98:98–103. doi:10.1016/j.diabres.2012.04.012
33. Harada K, Amano T, Uetani T, Yoshida T, Kato B, Kato M, et al. Association of inflammatory markers with the morphology and extent of coronary plaque as evaluated by 64-slice multidetector computed tomography in patients with stable coronary artery disease. Int J Cardiovasc Imaging (2013) 29:1149–58. doi:10.1007/s10554-013-0181-2
34. Kim HM, Lee BW, Song YM, Kim WJ, Chang HJ, Choi DH, et al. Potential association between coronary artery disease and the inflammatory biomarker YKL-40 in asymptomatic patients with type 2 diabetes mellitus. Cardiovasc Diabetol (2012) 11:84. doi:10.1186/1475-2840-11-84
35. von Scholten BJ, Reinhard H, Hansen TW, Schalkwijk CG, Stehouwer C, Parving HH, et al. Markers of inflammation and endothelial dysfunction are associated with incident cardiovascular disease, all-cause mortality, and progression of coronary calcification in type 2 diabetic patients with microalbuminuria. J Diabetes Complications (2016) 30:248–55. doi:10.1016/j.jdiacomp.2015.11.005
36. Zhang J, Lv Z, Zhao D, Liu L, Wan Y, Fan T, et al. Coronary plaque characteristics assessed by 256-slice coronary CT angiography and association with high-sensitivity C-reactive protein in symptomatic patients with type 2 diabetes. J Diabetes Res (2016) 2016:4365156. doi:10.1155/2016/4365156
38. Mirabelli P, Incoronato M, Coppola L, Infante T, Parente CA, Nicolai E, et al. SDN biobank: bioresource of human samples associated with functional and/or morphological bioimaging results for the study of oncological, cardiological, neurological, and metabolic diseases. Open J Biores (2017) 4:2. doi:10.5334/ojb.26
39. Tesche C, Plank F, De Cecco CN, Duguay TM, Albrecht MH, Varga-SzemesA Bayer RR, et al. Prognostic implications of coronary CT angiography-derived quantitative markers for the prediction of major adverse cardiac events. J Cardiovasc Comput Tomogr (2016) 10:458–65. doi:10.1016/j.jcct.2016.08.003
40. Joshi PH, Rinehart S, Vazquez G, Qian Z, Sharma A, Anderson H, et al. A peripheral blood gene expression score is associated with plaque volume and phenotype by intravascular ultrasound with radiofrequency backscatter analysis: results from the ATLANTA study. Cardiovasc Diagn Ther (2013) 3:5–14. doi:10.3978/j.issn.2223-3652.2013.01.02
41. Achenbach S, Ropers D, Hoffmann U, MacNeill B, Baum U, Pohle K, et al. Assessment of coronary remodeling in stenotic and nonstenotic coronary atherosclerotic lesions by multidetector spiral computed tomography. J Am Coll Cardiol (2004) 43:842–7. doi:10.1016/j.jacc.2003.09.053
42. Austen WG, Edwards JE, Frye RL, Gensini GG, Gott VL, Griffith LS, et al. A reporting system on patients evaluated for coronary artery disease. Report of the ad hoc committee for grading of coronary artery disease, council on cardiovascular surgery, American Heart Association. Circulation (1975) 51:5–40. doi:10.1161/01.CIR.51.4.5
43. Leipsic J, Abbara S, Achenbach S, Cury R, Earls JP, Mancini GBJ, et al. SCCT guidelines for the interpretation and reporting of coronary CT angiography: a report of the society of cardiovascular computed tomography guidelines committee. J Cardiovasc Comput Tomogr (2014) 8:342–58. doi:10.1016/j.jcct.2014.07.003
44. Schiano C, Rienzo M, Casamassimi A, Napoli C. Gene expression profile of the whole mediator complex in human osteosarcoma and normal osteoblasts. Med Oncol (2013) 30:739. doi:10.1007/s12032-013-0739-9
45. Rienzo M, Schiano C, Casamassimi A, Grimaldi V, Infante T, Napoli C. Identification of valid reference housekeeping genes for gene expression analysis in tumor neovascularization studies. Clin Transl Oncol (2013) 15:211–8. doi:10.1007/s12094-012-0904-1
47. Gao Y, Lu B, Sun ML, Hou ZH, Yu FF, Cao HL, et al. Comparison of atherosclerotic plaque by computed tomography angiography in patients with and without diabetes mellitus and with known or suspected coronary artery disease. Am J Cardiol (2011) 108:809–13. doi:10.1016/j.amjcard.2011.04.032
48. Chu ZG, Yang ZG, Dong ZH, Zhu ZY, Peng LQ, Shao H, et al. Characteristics of coronary artery disease in symptomatic type 2 diabetic patients: evaluation with CT angiography. Cardiovasc Diabetol (2010) 9:74. doi:10.1186/1475-2840-9-74
49. de AraújoGonçalves P, Garcia-Garcia HM, Carvalho MS, Dores H, Sousa PJ, Marques H, et al. Diabetes as an independent predictor of high atherosclerotic burden assessed by coronary computed tomography angiography: the coronary artery disease equivalent revisited. Int J Cardiovasc Imaging (2013) 29:1105–14. doi:10.1007/s10554-012-0168-4
50. Khazai B, Luo Y, Rosenberg S, Wingrove J, Budoff MJ. Coronary atherosclerotic plaque detected by computed tomographic angiography in subjects with diabetes compared to those without diabetes. PLoS One (2015) 10:e0143187. doi:10.1371/journal.pone.0143187
51. Palmieri V, Gravino E, Russo C, Salvati A, Lombardi C, Sauro R, et al. Coronary atherosclerosis burden by coronary computed tomography in type II diabetes with preclinical non-obstructive carotid atherosclerosis and without inducible myocardial ischemia. Diabetes Res Clin Pract (2017) 123:112–9. doi:10.1016/j.diabres.2016.11.024
52. Raggi P, Shaw LJ, Berman DS, Callister TQ. Prognostic value of coronary artery calcium screening in subjects with and without diabetes. J Am Coll Cardiol (2004) 43:1663–9. doi:10.1016/j.jacc.2003.09.068
53. Scholte AJ, Schuijf JD, Kharagjitsingh AV, Jukema JW, Pundziute G, van der Wall EE, et al. Prevalence of coronary artery disease and plaque morphology assessed by multi-slice computed tomography coronary angiography and calcium scoring in asymptomatic patients with type 2 diabetes. Heart (2008) 94:290–5. doi:10.1136/hrt.2007.121921
54. Pundziute G, Schuijf JD, Jukema JW, Boersma E, Scholte AJ, Kroft LJ, et al. Noninvasive assessment of plaque characteristics with multislice computed tomography coronary angiography in symptomatic diabetic patients. Diabetes Care (2007) 30:1113–9. doi:10.2337/dc06-2104
55. Reddy HK, Koshy SK, Foerst J, Sturek M. Remodeling of coronary arteries in diabetic patients-an intravascular ultrasound study. Echocardiography (2004) 21:139–44. doi:10.1111/j.0742-2822.2004.03014.x
56. Kim SH, Moon JY, Lim YM, Kim KH, Yang WI, Sung JH, et al. Association of insulin resistance and coronary artery remodeling: an intravascular ultrasound study. Cardiovasc Diabetol (2015) 14:74. doi:10.1186/s12933-015-0238-8
57. Tomizawa N, Nojo T, Inoh S, Nakamura S. Difference of coronary artery disease severity, extent and plaque characteristics between patients with hypertension, diabetes mellitus or dyslipidemia. Int J Cardiovasc Imaging (2015) 31:205–12. doi:10.1007/s10554-014-0542-5
59. Pasterkamp G, Schoneveld AH, van der Wal AC, Haudenschild CC, Clarijs RJ, Becker AE, et al. Relation of arterial geometry to luminal narrowing and histologic markers for plaque vulnerability: the remodeling paradox. J Am Coll Cardiol (1998) 32:655–62. doi:10.1016/S0735-1097(98)00304-0
60. Peng JR, Lu TT, Chang HT, Ge X, Huang B, Li WM. Elevated levels of plasma superoxide dismutases 1 and 2 in patients with coronary artery disease. Biomed Res Int (2016) 2016:3708905. doi:10.1155/2016/3708905
61. Abdullah MH, Othman Z, Noor HM, Arshad SS, Yusof AK, Jamal R, et al. Peripheral blood gene expression profile of atherosclerotic coronary artery disease in patients of different ethnicity in Malaysia. J Cardiol (2012) 60:192–203. doi:10.1016/j.jjcc.2012.05.009
63. Wójcicka G, Jamroz-Wiśniewska A, Horoszewicz K, Bełtowski J. Liver X receptors (LXRs). Part I: structure, function, regulation of activity, and role in lipid metabolism. Postepy Hig Med Dosw (2007) 61:736–59.
64. Dahlman I, Nilsson M, Gu HF, Lecoeur C, Efendic S, Ostenson CG, et al. Functional and genetic analysis in type 2 diabetes of liver X receptor alleles – a cohort study. BMC Med Genet (2009) 17(10):27. doi:10.1186/1471-2350-10-27
65. Xi W, Zhou Y, Lv S, Gao Q, Bu G, Wang Y, et al. Plasma hyaluronan and collateral development in patients with coronary artery disease. Coron Artery Dis (2010) 21:228–32. doi:10.1097/MCA.0b013e328338ccf3
67. Douglas PS, De Bruyne B, Pontone G, Patel MR, Norgaard BL, Byrne RA, et al. 1-year outcomes of FFRCT-guided care in patients with suspected coronary disease: the PLATFORM study. J Am Coll Cardiol (2016) 68:435–45. doi:10.1016/j.jacc.2016.05.057
68. Blanke P, Naoum C, Ahmadi A, Cheruvu C, Soon J, Arepalli C, et al. Long-term prognostic utility of coronary CT angiography in stable patients with diabetes mellitus. JACC Cardiovasc Imaging (2016) 9:1280–8. doi:10.1016/j.jcmg.2015.12.027
69. Kerkmeijer LS, Farhan S, Mehran R, Dangas GD. Diabetes mellitus and multivessel coronary artery disease: an ongoing battle for an ideal treatment strategy. Ann Transl Med (2017) 5:261. doi:10.21037/atm.2017.03.92
70. Haffner SM, Lehto S, Rönnemaa T, Pyörälä K, Laakso M. Mortality from coronary heart disease in subjects with type 2 diabetes and in nondiabetic subjects with and without prior myocardial infarction. N Engl J Med (1998) 339:229–34. doi:10.1056/NEJM199807233390404
Keywords: type 2 diabetes, coronary artery disease, computed tomography coronary angiography, biomarkers, atherosclerosis
Citation: Infante T, Forte E, Aiello M, Salvatore M and Cavaliere C (2017) In Vivo and In Vitro Analysis in Coronary Artery Disease Related to Type 2 Diabetes. Front. Endocrinol. 8:209. doi: 10.3389/fendo.2017.00209
Received: 23 February 2017; Accepted: 08 August 2017;
Published: 21 August 2017
Edited by:Gaetano Santulli, Columbia University, United States
Reviewed by:Carla Contaldi, Northwestern University, Italy
Jessica Gambardella, University of Salerno, Italy
Copyright: © 2017 Infante, Forte, Aiello, Salvatore and Cavaliere. 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) or licensor 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: Ernesto Forte, firstname.lastname@example.org
†These authors share first author position.