Next Article in Journal
Analyses of AUC(0–12) and C0 Compliances within Therapeutic Ranges in Kidney Recipients Receiving Cyclosporine or Tacrolimus
Next Article in Special Issue
Serum Levels of Soluble Triggering Receptor Expressed on Myeloid Cells-1 Associated with the Severity and Outcome of Acute Ischemic Stroke
Previous Article in Journal
Considerations for Cannabis Use to Treat Pain in Sickle Cell Disease
Previous Article in Special Issue
Anatomical Considerations for Endovascular Intervention for Extracranial Carotid Disease: A Review of the Literature and Recommended Guidelines
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Sequence Variation in the DDAH1 Gene Predisposes for Delayed Cerebral Ischemia in Subarachnoidal Hemorrhage

1
Institute of Clinical Pharmacology and Toxicology, University Medical Center Hamburg-Eppendorf, 20246 Hamburg, Germany
2
Department of Anesthesiology, University Medical Center Hamburg-Eppendorf, 20246 Hamburg, Germany
*
Author to whom correspondence should be addressed.
J. Clin. Med. 2020, 9(12), 3900; https://doi.org/10.3390/jcm9123900
Submission received: 30 October 2020 / Revised: 19 November 2020 / Accepted: 27 November 2020 / Published: 1 December 2020
(This article belongs to the Special Issue Novel Approaches for the Treatment of Cerebrovascular Disorders)

Abstract

:
Delayed cerebral ischemia (DCI) often causes poor long-term neurological outcome after subarachnoidal hemorrhage (SAH). Asymmetric dimethylarginine (ADMA) inhibits nitric oxide synthase (NOS) and is associated with DCI after SAH. We studied single nucleotide polymorphisms (SNPs) in the NOS3, DDAH1, DDAH2, PRMT1, and AGXT2 genes that are part of the L-arginine–ADMA–NO pathway, and their association with DCI. We measured L-arginine, ADMA and symmetric dimethylarginine (SDMA) in plasma and cerebrospinal fluid (CSF) of 51 SAH patients at admission; follow-up was until 30 days post-discharge. The primary outcome was the incidence of DCI, defined as new infarctions on cranial computed tomography, which occurred in 18 of 51 patients. Clinical scores did not significantly differ in patients with or without DCI. However, DCI patients had higher plasma ADMA and SDMA levels and higher CSF SDMA levels at admission. DDAH1 SNPs were associated with plasma ADMA, whilst AGXT2 SNPs were associated with plasma SDMA. Carriers of the minor allele of DDAH1 rs233112 had a significantly increased relative risk of DCI (Relative Risk = 2.61 (1.25–5.43), p = 0.002). We conclude that the DDAH1 gene is associated with ADMA concentration and the incidence of DCI in SAH patients, suggesting a pathophysiological link between gene, biomarker, and clinical outcome in patients with SAH.

1. Introduction

Subarachnoidal hemorrhage (SAH) is a major neurosurgical emergency. Vasospasm and cerebral ischemia often occur with some delay during intensive care unit (ICU) treatment of SAH; delayed cerebral ischemia (DCI) is associated with poor clinical outcome [1]. Despite clinical scores and advanced imaging techniques, identification of patients who will develop DCI or not has remained difficult [2].
The L-arginine–nitric oxide (NO) pathway is a major regulator of cerebral vascular tone [3]. NO-mediated cerebral vasodilation is impaired in SAH, thereby eliciting vasospasm and, possibly, DCI [4]. Multiple mechanisms contribute to this reduction in NO. Cell-free hemoglobin in the subarachnoidal space has been shown to scavenge NO [5]; removal of hemoglobin by haptoglobin prevents hemoglobin-induced cerebral vasospasm [6]. Various strategies to restore pulmonary NO levels, e.g., the use of inhaled NO [7], exogenous NO donors, or phosphodiesterase-5 (PDE5) inhibitors, have been shown to mitigate the risk of cerebral vasospasm and ischemia after SAH. One experimental study in haptoglobin 2-2 transgenic mice tested the effects of L-citrulline that is converted to L-arginine after absorption, and found improved cerebral perfusion and neurological outcome [8]. Restoring eNOS substrate concentration has been shown to be effective only when this enzyme is blocked by a competitive inhibitor [9]. Asymmetric dimethylarginine (ADMA) competitively inhibits NO synthesis by displacing the substrate, L-arginine, from the enzyme. We demonstrated in a previous study that higher levels of the endogenous inhibitor of NO synthase (NOS), ADMA, and its congener, symmetric dimethylarginine (SDMA), are markers of DCI and poor neurological outcome in SAH patients [10].
Both ADMA and SDMA are formed during posttranslational protein methylation by protein arginine N-methyltransferases (PRMTs) [11]. In many organs including the vasculature, ADMA is the main product of this process; ADMA is cleaved by either of two isoforms of dimethylarginine dimethylaminohydrolase (DDAH1 and DDAH2) [12]. SDMA is not degraded by DDAHs, but both SDMA and ADMA are substrates of alanine glyoxylate aminotransferase-2 (AGXT2) [13,14].
Figure 1 provides an overview of enzymatic pathways regulating L-arginine and dimethylarginine levels. We have described this pathway and its regulation in detail before [15].
Clinically, ADMA is associated with cardiovascular events and mortality [16,17], whilst SDMA has been associated with mortality after ischemic stroke [18,19]. ADMA is elevated in plasma and/or cerebrospinal fluid (CSF) of patients with SAH [20,21,22] as well as in animal models of SAH [23]. In our previous study [10], we found ADMA and SDMA levels in plasma and CSF of SAH patients to be associated with the incidence of DCI. Therefore, elevated levels of dimethylarginines may result in dysfunction of endothelial NO production, impaired cerebrovascular autoregulation, and vasospasm.
Genotyping studies showed that single nucleotide polymorphisms (SNPs) in the NOS3 gene are associated with an increased incidence of DCI after SAH [24,25] and with poor long-term recovery from SAH [26]. Several single nucleotide polymorphisms (SNPs) in the genes encoding for DDAH1, DDAH2, AGXT2, and PRMT1 have been shown to relate to the plasma concentrations of ADMA or SDMA, respectively, and to be associated with clinical findings of cardiovascular diseases [13,27,28,29]. Based on this data, we hypothesized that genetic variation caused by SNPs in the genes encoding enzymes relating to the L-arginine–dimethylarginine metabolic pathways may be associated with elevated baseline levels of ADMA and SDMA and may contribute to individual susceptibility to DCI after SAH. To test this hypothesis, we genotyped representative SNPs in each of the genes of interest for the L-arginine–dimethylarginine–NO pathway in a cohort of patients admitted with acute SAH.
Figure 1. Schematic representation of biosynthesis and metabolism of asymmetric and symmetric dimethylarginine. L-arginine residues within specific proteins are subject to methylation by protein arginine N-methyltransferases (PRMTs). After protein hydrolysis, asymmetric (ADMA) and symmetric dimethylarginine (SDMA) are released. ADMA is a competitive inhibitor of nitric oxide synthases (NOS). ADMA, but not SDMA, is degraded by dimethylarginine dimethylaminohydrolases (DDAH1 and DDAH2) into L-citrulline and dimethylamine (DMA). Both dimethylarginines may be cleaved by an alternative pathway through alanine glyoxylate aminotransferase 2 (AGXT2), resulting in the formation of symmetric or asymmetric dimethylguanidinovaleric acid (DMGV and DM’GV). Figure reproduced with permission from [30] under CC BY-NC (SAGE Publishers).
Figure 1. Schematic representation of biosynthesis and metabolism of asymmetric and symmetric dimethylarginine. L-arginine residues within specific proteins are subject to methylation by protein arginine N-methyltransferases (PRMTs). After protein hydrolysis, asymmetric (ADMA) and symmetric dimethylarginine (SDMA) are released. ADMA is a competitive inhibitor of nitric oxide synthases (NOS). ADMA, but not SDMA, is degraded by dimethylarginine dimethylaminohydrolases (DDAH1 and DDAH2) into L-citrulline and dimethylamine (DMA). Both dimethylarginines may be cleaved by an alternative pathway through alanine glyoxylate aminotransferase 2 (AGXT2), resulting in the formation of symmetric or asymmetric dimethylguanidinovaleric acid (DMGV and DM’GV). Figure reproduced with permission from [30] under CC BY-NC (SAGE Publishers).
Jcm 09 03900 g001

2. Materials and Methods

2.1. Study Participants and Protocol

We included 51 consecutive patients who had been admitted to the Intensive Care Unit of the University Medical Center Hamburg-Eppendorf with acute subarachnoidal hemorrhage. The aneurysmal pathogenesis was diagnosed using CT angiography or digital subtraction angiography; in 47 patients, an aneurysm was detected, and in four patients a perimesencephalic subarachnoidal hemorrhage was diagnosed. In 38 patients, external ventricular drain (EVD) insertion was performed on the day of admission because elevation of intracranial pressure was documented. The local ethics committee approved the study protocol (Ethics Committee of the Board of Physicians of Hamburg; PV3622), and written informed consent was obtained from all participants or their next of kin.
The aneurysms were secured by intravascular coiling or surgical clipping as early as possible, at least within 24 h after admission. In eight cases, a non-interventional treatment approach was chosen. Subsequent adjunctive therapies and CT scans were administered when appropriate in adherence with international guidelines defined by the Neurocritical Care Society’s Multidisciplinary Consensus Conference and the American Heart Association (AHA) [31,32]. In accordance with local protocols, oral or intravenous nimodipine was used to prevent cerebral vasospasm in 47 patients. During poor neurological condition or during mechanical ventilation, nimodipine was given intravenously, otherwise orally. Blood gas analysis was performed several times per day, and at least one blood sample per day for routine laboratory parameters was drawn during intensive care therapy, depending on clinical needs. Normoventilation (PaCO2 35–45 mmHg) and normovolemia were controlled closely, and blood pressure was monitored continuously by invasive monitoring. Norepinephrine was primarily used for therapy of episodes of hypotension.

2.2. Clinical Assessment and Follow-up

Patients were included on the day of bleeding and were reviewed on days three, six, eight, 12, and 15 after hemorrhage, including a standardized medical history, blood and CSF samples drawn from arterial catheters and external ventricular drains, respectively, and the Glasgow Coma Scale (GCS). CSF samples were only collected from patients with EVD inserted for clinical needs; for ethical reasons, no lumbar punctures were performed solely for study purposes.
At admission, clinical assessment of patients was performed according to standardized procedures including the World Federation of Neurological Surgeons (WFNS) score [33], the GCS [34], and the Hunt and Hess score [35] by board-certified neurosurgeons. Severe clinical condition was defined by a Hunt and Hess score of 3–5. For all patients, the initial CT scan was assessed according to the Fisher classification [36].
The primary outcome parameter was DCI, which was assessed by cranial CT or magnetic resonance imaging performed when signs of neurological impairment became visible to identify new cerebral infarction according to current standards [37]. Transcranial Doppler (TCD) measurements of cerebral arterial flow were taken daily by the same investigator to identify cerebral vasospasm. For statistical analyses, we calculated the maximum increase in mean flow velocity (MFV) related to the day of bleeding (baseline). TCD vasospasm was defined as MFV > 140 cm/s [38].
At discharge, neurological outcome was determined using the Glasgow Outcome Scale (GOS). Poor neurological outcome was defined as GOS < 4. In addition, outcome including the GOS was assessed at 3 months after discharge by patient visits or telephone interviews.

2.3. Selection of SNPs

We performed an extensive literature search using PubMed and the National Center for Biotechnology Information (NCBI) SNP database for publications linking SNPs in the DDAH1, DDAH2, ARG1, ARG2, AGXT2, NOS3, and PRMT1 genes to outcome parameters related to the pathophysiology of subarachnoidal hemorrhage, ischemic stroke, and vasospasm. These genes code for the key enzymes in the L-arginine–ADMA/SDMA–NO pathway (Figure 1) [30]. Subsequently, SNPs in the same gene were tested for linkage disequilibrium (LD) using the LDLink database of the National Cancer Institute, USA [39], and the list of SNPs was reduced to the best described SNP for this specific group when strong LD was present; only SNPs that are not in LD were included in the analyses (Table A1).

2.4. DNA Isolation and SNP Genotyping

DNA from blood samples was isolated using QIAamp DNA Blood Mini Kit according to the manufacturer’s protocol (Qiagen, Hilden, Germany). DNA was eluted from the columns in a volume of 20 µL. To further purify and concentrate the isolated DNA for SNP genotyping, DNA was precipitated by ethanol precipitation. The DNA pellet was washed once and the resuspended in 20–100 µL H20. DNA concentration was determined using a nanophotometer NP60 (Implen, Munich, Germany).
SNP genotyping was performed using single-tube human TaqMan SNP Genotyping Assays (Thermofisher, Waltham, MS, USA). Each reaction mix contained 10 ng template DNA, 0.5µL specific TaqMan SNP Genotyping Assay and 5 µL of 2xTaqPath ProAmp™ Mastermix in a total reaction volume of 10 µL. The polymerase chain reaction (PCR) was performed in a QuantStudio™ 5 Real-Time PCR System (Thermofisher, Waltham, MS, USA) using the following PCR program: pre-read of 30 s at 60 °C, enzyme activation at 95 °C for 5 min, 40 cycles of denaturation (5 s, 95 °C) and annealing (30 s, 60 °C), followed by a last post-read of 30 s at 60 °C. Allelic calls were identified by Quant Studio Design and Analysis Software (Thermofisher, Waltham, MS, USA).

2.5. Biochemical Analyses

At baseline, non-fasting blood and CSF samples were drawn from arterial catheters (N = 51) and external ventricular drains (when available, N = 38), respectively. For ethical reasons, no patient received any interventional measures for study-related sample collection. Samples were centrifuged immediately and stored at −20 °C until analysis. ADMA, SDMA, and L-arginine were determined by liquid chromatography-tandem mass spectrometry (LC–MS/MS) using a previously validated method [40]. In brief, 25 μL of plasma or CSF was diluted with stable isotope-labeled internal standards. Proteins were precipitated with methanol; the guanidine compounds were converted to their butyl esters and analyzed by LC–MS/MS (Varian 1200 MS, Agilent Technologies, Santa Clara, CA, USA). Quantification was performed by calculation of peak area ratios and calibration with known concentrations of analytes in dialyzed EDTA plasma. The analytical range of the method was validated from 0.05 to 4 μmol/L for ADMA and SDMA, respectively, and from 0.5 to 250 μmol/L for L-arginine, and mean coefficients of variation were ≤5% for all analytes. All other laboratory values were measured using routine clinical laboratory methods.

2.6. Statistical Analyses

All statistical analyses were performed using SPSS (version 21; IBM Corporation, Armonk, NY, USA). All variables were tested for normal distribution using the Kolmogorov–Smirnov test. Differences between groups were tested for significance by using either the nonparametric Mann–Whitney U test for two groups or the Kruskal–Wallis analysis of variance for more than two groups. Correlations were calculated using linear regression or the Spearman test. Allele frequencies of the genes of interest (GOI) were compared to the European reference population or between groups using contingency tables and Fisher’s exact test. Biomarker levels between genotypes were tested for statistically significant differences by ANOVA followed by the Scheffé f-test, and associations of genotypes with outcome (DCI) were tested for significance by the Chi2 test. Haplotypes were constructed for the DDAH1 gene based upon the frequency of homozygosity for the major and minor alleles, respectively. The associations of DDAH1 haplotypes with outcome (DCI) were tested for significance by the Chi2 test. Data are presented as median with 25th and 75th percentiles, or as mean with standard error of the mean. For all tests, p < 0.05 was considered significant.

3. Results

3.1. Baseline Characteristics and Clinical Outcome of the Patients

We included in this study 51 patients who were admitted to our Emergency Department with evidence for subarachnoidal hemorrhage. Patients had a mean age of 54.8 ± 13.1 years; 25% of them were men. A total of 43 patients were treated interventionally with either clipping or coiling of the aneurysmal location; 47 patients received nimodipine treatment. The Glasgow Coma Scale was 8.6 ± 5.6 on the day of admission, the WFNS score was 3.4 ± 1.7, and the Hunt and Hess index was 3.1 ± 1.4. Further demographic and anthropometric characteristics of the cohort are listed in Table 1, including risk factor profiles and long-term medical treatment. Thirty-eight patients had an EVD placed on the day of admission for clinical reasons. We found a significantly poorer survival rate in patients with EVD (Relative Risk (RR), 3.897 (95% CI 1.002–15.160), p = 0.039), but no significant associations of EVD with the incidence of DCI (RR, 1.271 (95% CI 0.858–1.883), p = 0.286) nor vasospasm (RR, 1.542 (95% CI 0.710–3.347), p = 0.305).
A total of 18 patients developed DCI during the course of in-patient treatment. At admission, there were no significant differences in any of the demographic and anthropometric parameters between patients who later developed DCI or not (Table 1). Initially, patients who later developed DCI scored even better in the Glasgow Coma Scale than patients who did not (10.9 ± 4.6 vs. 7.2 ± 5.7; p = 0.02; Table 2), whilst there were no significant differences in the Hunt and Hess Score, Fisher Score, and WFNS Score (all p = n.s.; Table 2). However, the clinical outcome of patients with DCI was more unfavorable than for patients without DCI: There was a trend towards a higher lethality, longer ICU treatment, and longer mechanical ventilation (all p = n.s.), as well as a significantly higher rate of re-intervention (p = 0.047; Table 2). The Glasgow Outcome Score at discharge and three months later tended to be lower for patients with DCI (2.8 ± 1.6 vs. 3.6 ± 1.5 at discharge, p = 0.056; 3.1 ± 1.1 vs. 4.0 ± 1.0 three months later; p = 0.071). Vasospasms occurred slightly, but not significantly more frequently in patients with DCI (10 vs. 8 patients) than without DCI (12 vs. 21 patients; p = n.s.).

3.2. Association of L-Arginine and Dimethylarginines with DCI

For L-arginine, there was no significant difference in plasma or CSF between patients with or without DCI. By contrast, the plasma concentrations of ADMA (0.58 ± 0.15 vs. 0.49 ± 0.17 µmol/L; p = 0.04) and SDMA at admission (0.73 ± 0.26 vs. 0.52 ± 0.17 µmol/L; p = 0.03) were significantly higher in patients who later developed DCI than in those who did not. SDMA concentration was also significantly higher in CSF of patients with DCI than without (0.31 ± 0.10 vs. 0.25 ± 0.08 µmol/L; p = 0.03; Figure 2).

3.3. Allele Frequencies of Genes of Interest in the Study Cohort Compared to the European Reference Population

By comparison with the European reference population of the 1000 Genomes project, there were significant deviations from the expected distribution of alleles in four of the SNPs that we studied: For NOS3 rs891512, the major allele was more frequent in our cohort, whilst for the DDAH1 rs1241321 and DDAH1 rs480414 SNPs and the DDAH2 rs805304 SNP, the minor allele was more frequent in our cohort (all p < 0.05). Table A2 shows allele frequencies of all analyzed SNPs.

3.4. Association of SNPs with Biomarkers in Plasma and in Cerebrospinal Fluid

There was a significant association of the DDAH1 rs233112 SNP with ADMA plasma concentration, whilst this SNP was not significantly associated with CSF ADMA concentration (Figure 3), nor were any of the other SNPs that we analyzed. By contrast, the DDAH1 SNPs rs480414 and rs1241321 were not significantly associated with ADMA plasma concentration (Figure 3c,d). The AGXT2 SNP rs16899974 was significantly associated with SDMA concentration in plasma but not in CSF (Figure 4a,b) and showed a weak, non-significant association with ADMA concentration in plasma, but not in CSF (Figure 4c,d).

3.5. Association of SNPs with the Incidence of Delayed Cerebral Ischemia

All three SNPs in the DDAH1 gene were significantly associated with the incidence of DCI. In each of the three SNPs, patients with DCI more frequently carried the minor allele than patients without DCI (Figure 5). Patients who were homozygous for the minor allele of DDAH1 rs233112 had a relative risk of 2.61 (95% CI, 1.25–5.43) of developing DCI (p = 0.002) as compared to patients who were homozygous for the major allele. By contrast, SNPs in the DDAH2, AGXT2, PRMT1, or NOS3 genes showed no significant associations with DCI.
When we combined genotypes for the three DDAH1 SNPs and analyzed the associations of DDAH1 haplotypes with outcome, we found that as compared to patients homozygous for all three major alleles, individuals homozygous for ≥1 minor allele had a significantly higher relative risk for DCI (RR, 1.316 (95% CI 1.072–1.614), p = 0.011) and for poor neurological outcome at ICU discharge (RR, 2.273 (95% CI 1.681–3.072), p < 0.0001) and after 3 months (RR, 1.299 (95% CI 1.110–1.519), p = 0.001). No significant different was observed for the relative risk of death and vasospasm, respectively (Table 3). Patients who were not homozygous for any of the DDAH1 SNPs had an intermediate risk (Table 3).

4. Discussion

Our study has three major findings: Firstly, the plasma concentrations of ADMA and SDMA and the CSF concentration of SDMA at the time of admission are positively associated with the incidence of DCI after subarachnoidal hemorrhage. Secondly, common genetic polymorphisms in the DDAH1 gene are responsible for differences in baseline ADMA plasma concentration, and common genetic polymorphisms in the AGXT2 gene are responsible for differences in plasma SDMA concentration; by contrast, the CSF concentrations of these metabolites were not significantly associated with these SNPs. Thirdly, and most importantly, patients with DCI more frequently carried the minor alleles of any of the three common DDAH1 SNPs that we analyzed; patients homozygous for the minor allele of DDAH1 rs233112 had a significantly higher probability of DCI than patients homozygous for the major allele. Patients carrying the DDAH1 haplotype with the most minor alleles had the highest relative risk for DCI and poor neurological outcome as compared to patients homozygous for the major allele of all three DDAH1 SNPs.
ADMA is an endogenous inhibitor of nitric oxide synthesis. Its plasma levels are positively associated with the incidence of major adverse cardiovascular events and mortality [16]. By inhibiting endothelial NO release [41], high ADMA concentration may cause vasospasm and contribute to ischemia. We have reported a significant association of plasma ADMA with magnetic resonance imaging markers of subclinical vascular brain injury in the Framingham Offspring cohort [42]. These observations are in line with our current results. We have reported before that ADMA plasma concentration increases significantly more in SAH patients developing DCI than in those who do not [10]. However, the initial L-arginine/ADMA ratio (marking the availability of NOS substrate (L-arginine) over NOS inhibitor (ADMA)) at admission was significantly lower in patients who later developed DCI than in those who did not. We therefore hypothesized that genetic predisposition may underlie variable ADMA and SDMA concentrations at admission and be related to outcomes of SAH patients.
In a previous genome-wide association study, we demonstrated that plasma ADMA concentration is associated with the DDAH1 gene, whilst SDMA plasma concentration is associated with the AGXT2 gene in the general population [13]. In our present study, we genotyped common SNPs in genes encoding crucial enzymes of the L-arginine–ADMA–NO pathway, and we confirm in this critically ill population a link between DDAH1 and ADMA plasma concentration and between AGXT2 and SDMA plasma concentration.
Another observation of our study is that the plasma concentrations, but not the CSF levels, of ADMA and SDMA were significantly associated with the DDAH1 and AGXT2 genes, respectively. DDAH1 is highly expressed in liver and kidneys [12] and AGXT2 is primarily located in the renal proximal tubular epithelial cells [43]. These results suggest that ADMA and SDMA circulating in plasma may be cleared by these two enzymatic pathways, whilst dimethylarginines in CSF may not be directly affected by these enzymes. This conclusion is in line with our previous observation that plasma and CSF ADMA concentrations are differentially regulated in Alzheimer’s disease patients [44]. Differential regulation of ADMA in plasma versus CSF was also reported after experimental middle cerebral artery occlusion in rats [45]; ADMA and SDMA concentrations in CSF were found to be elevated immediately after hyperacute ischemic stroke and associated with stroke severity [46].
SDMA does not interfere with NO synthesis directly, but it indirectly modulates NO generation by limiting the intracellular availability of L-arginine via inhibition of the cationic amino acid transporter [47]. Several studies have linked SDMA plasma concentration with ischemic cerebral stroke: In a study of 394 survivors of the first 30 days after acute ischemic stroke, we found that SDMA is associated with total mortality during 7.4 years of follow-up [19]. In another study, we found that cumulative survival during the first 30 days after ischemic stroke decreased significantly with increasing tertiles of plasma SDMA [18]. SDMA is formed during posttranslational methylation of proteins by type II PRMTs [11]; intriguingly, myelin-basic protein, an essential component of myelin sheaths in neuronal tissues, is subject to symmetric dimethylation [48]. Thus, SDMA might be a marker of neurological damage, and in our present study high CSF SDMA levels at admission may point to early neuronal damage in SAH. This observation clearly deserves further investigation.
Beyond its association with ADMA levels, there was a strong association of the DDAH1 gene with outcome after SAH: All three analyzed DDAH1 SNPs were associated with the incidence of DCI, and as compared to a combined genotype (haplotype) containing homozygosity for all three SNPs, haplotypes with an increasing number of minor alleles showed increasingly worse outcome. Thus, functionality of DDAH1 may be an important, in part genetically influenced determinant of the incidence of cerebral ischemia after SAH.
Whilst presenting some exciting observations, our study has limitations. Firstly, the small sample size limits our ability to extrapolate our findings to various levels of severity of SAH. Specifically, the significant deviation of allele distribution in our critically ill cohort from the expected allele frequencies may suggest that genes regulating vascular NO production may predispose for SAH. Clearly, larger cohorts need to be genotyped to confirm and extend our findings. Secondly, due to the strict limitations linked to a clinical study in critically ill patients, we had no opportunity to investigate the mechanisms that may underlie the possible differential regulation of plasma and CSF dimethylarginine levels. This might be better addressed in animal models.

5. Conclusions

In conclusion, our study provides evidence that the L-arginine–ADMA–NO metabolic pathway is involved in the pathophysiology of cerebrovascular ischemic complications after subarachnoidal hemorrhage. The DDAH1 gene, probably via determining the metabolism of ADMA, is prominently involved in mediating this risk. Carriers of genetic variants of the DDAH1 gene that lead to malfunction of the enzyme may carry a genetic risk of DCI after SAH.

Author Contributions

Conceptualization, J.H., D.A. and R.B.; Data curation, D.A., M.S.-S., T.B., J.Z. and R.B.; Formal analysis, J.H., D.A., M.S.-S., T.B., J.Z. and R.B.; Methodology, J.H. and R.B.; Supervision, J.H. and R.B.; Writing—original draft, J.H.; Writing—review and editing, D.A., M.S.-S., T.B., J.Z. and R.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

We gratefully thank Mariola Kastner for excellent technical assistance.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Appendix A.1. Selection and Clinical Significance of the Selected Single Nucleotide Polymorphisms

In this study, we aimed to investigate whether genetic variance in the genes encoding for enzymes relating to the L-arginine–dimethylarginine metabolic pathways may be associated with elevated baseline levels of ADMA and SDMA and may contribute to individual susceptibility to DCI after SAH. A number of enzymes metabolize the biomarkers associated with this pathway (see Figure 1 in the main manuscript): Endothelial nitric oxide synthase (eNOS) is the enzyme that generates nitric oxide (NO), the major vasodilator mediator, from L-arginine. Another enzyme that metabolizes L-arginine is arginase. By converting L-arginine into L-ornithine, arginase may reduce the availability of substrate for NO synthesis and thereby indirectly influence NO production. Arginase is present in two isoforms (arginase-1 and arginase-2). NOS is inhibited by methylarginines; NG-monomethyl-L-arginine (L-NMMA), asymmetric and symmetric dimethylarginine (ADMA and SDMA) are present endogenously. They are formed during the process of posttranslational protein methylation by a family of enzymes called protein arginine N-methyltransferases (PRMTs). PRMT1 is the most abundant isoform of these enzymes that is present in many cell types. L-NMMA and ADMA are metabolized by either of two isoforms of dimethylarginine dimethylaminohydrolase (DDAH1 and DDAH2), which are present in different tissues and apparently underlie different mechanisms of regulation. ADMA and SDMA, in turn, are both substrates of alanine glyoxylate aminotransferase-2 (AGXT2). We have described this pathway and its regulation in detail before [15].
Genetic variability studies on the corresponding genes have been performed previously by numerous groups of investigators. Table A1 gives an overview of the SNPs that were selected for this study, including references to previous studies linking variability to outcome parameters related to the pathophysiology of subarachnoidal hemorrhage, ischemic stroke, and vasospasm. For some of the SNPs (e.g., in the NOS3, DDAH1, DDAH2, and AGXT2 genes), we made a selection of references relating the SNP to clinical parameters linked to the topic of our present study; for some other SNPs, there was comparatively little published information available (e.g., in the ARG1, ARG2, and PRMT1 genes). For the NOS3 gene, a vast number of different SNPs have been described; however, many of these are in linkage disequilibrium (LD) based on information available in the LDLink database of the National Cancer Institute, USA [39], so we selected one tag SNP out of a specific group. For ARG1, ARG2, and PRMT1, we used the same strategy of testing for LD in order to identify SNPs spread over the whole gene even if no clinical associations have been published yet.

Appendix A.2. Construction of DDAH1 Haplotypes

To assess the association of the DDAH1 gene with clinical outcome after SAH beyond the effects of individual SNPs in this gene, we constructed haplotypes based upon the numbers of homozygosities in these three gene loci. Haplotype 1 was defined as homozygous presence of the major allele for all three loci; this haplotype was used as reference in the statistical analyses (n = 8). Haplotype 2 was an intermediate haplotype, comprising individuals who were homozygous for the minor allele in none of the three DDAH1 loci and not simultaneously homozygous for the major alleles of all three loci (n = 22). Haplotype 3 comprised individuals who carried at least one homozygosity for the minor allele in any of the loci (n = 21). This classification allowed us to study the gene dose effect of DDAH1 on outcome, with the frequency of the major alleles gradually decreasing from haplotype 1 to haplotype 3. The associations of haplotypes with the following clinical outcome measures were assessed by Chi2 test: incident DCI, death, incident cerebral vasospasm, GOS at discharge, GOS at three months after hospital discharge.
Table A1. Summary of SNPs selected for this study and their previously published associations with clinical parameters.
Table A1. Summary of SNPs selected for this study and their previously published associations with clinical parameters.
Gene/SNPPositionGenetic ConsequenceClinical RationaleRelated References
NOS3
rs1799983Exon 8Missense variantMinor allele associated with worse outcome after SAH and risk of ischemic stroke.[26,49,50,51,52]
rs2070744Intron 1Intron variantMinor allele associated with increased risk of DCI and cardiac instability in patients with SAH; proposed as a susceptibility locus for cerebral infarction.[53,54,55]
rs891512Intron 23Intron variantMinor allele associated with ankle-brachial index, with higher systolic and diastolic blood pressure, and risk of coronary artery disease.[56,57,58,59,60]
DDAH1
rs1241321Intron 1Intron variantMinor allele associated with higher risk of all-cause mortality and the combined endpoint of cardiovascular death, nonfatal myocardial infarction and stroke.[61]
rs233112Exon 63′ UTR variantMinor allele associated with arterial stiffness and pulse wave reflection.[62]
rs480414Intron 1Intron variantMinor allele associated with a lower incidence of pulmonary hypertension in patients with bronchopulmonary dysplasia.[63]
DDAH2
rs805304 2 kb upstream variantMinor allele associated with higher prevalence of hypertension, higher ADMA concentration in diabetic renal impairment, and decreased risk of myocardial infarction.[64,65,66]
rs2272592 2 kb upstream variantMinor allele associated with type 2 diabetes.[67]
ARG1
rs2246012Intron 2Intron variantNo published information available.
rs2781667Intron 1Intron variantMajor allele associated with higher severity of erectile dysfunction.[68]
ARG2
rs3742879Intron 7Intron variantMajor allele associated with lower exhaled nitric oxide in children, and more severe airway obstruction in patients with asthma.[69,70]
rs3759757 8 kb upstream variantNo published information available.
AGXT2
rs37369Exon 4Missense variantMinor allele associated with higher methylarginine levels, higher diastolic blood pressure, and higher risk of coronary artery disease.[71,72,73]
rs16899974Exon 14Missense variantMinor allele associated with higher levels of SDMA, with atrial fibrillation and ischemic stroke, and higher risk of coronary artery disease.[71,73,74]
PRMT1
rs10415880Intron 8Intron variantMinor allele associated with arteriovenous fistula malfunction in male hemodialysis patients.[75]
rs975484 Upstream transcript variantNo published information available.
Abbreviations: NOS3, endothelial nitric oxide synthase; DDAH1, dimethylarginine dimethylaminohydrolase-1; DDAH2, dimethylarginine dimethylaminohydrolase-2; ARG1, arginase-1; ARG2, arginase-2; AGXT2, alanine glyoxylate aminotransferase-2; PRMT1, protein-arginine N-methyltransferase-1.

Appendix B

Table A2. Allele frequencies of studied SNPs.
Table A2. Allele frequencies of studied SNPs.
Gene/SNPMajor/Minor AlleleAllele Frequency MeasuredAllele Frequency Expected #p
NOS3
rs1799983G/T0.706/0.2940.656/0.3440.223
rs2070744T/C0.598/0.4020.562/0.4380.530
rs891512G/A0.843/0.1570.751/0.2490.039
DDAH1
rs1241321A/G0.598/0.4020.712/0.2880.023
rs233112T/C0.667/0.3330.617/0.3830.338
rs480414G/A0.598/0.4020.698/0.3020.043
DDAH2
rs805304T/G0.529/0.4710.639/0.3610.032
rs2272592C/T0.892/0.1080.844/0.1560.246
ARG1
rs2246012C/T0.853/0.1470.680/0.3200.677
rs2781667T/C0.676/0.3240.829/0.1711.000
ARG2
rs3742879G/C0.696/0.3040.741/0.2590.911
rs3759757A/G0.755/0.2450.687/0.3130.813
AGXT2
rs37369C/A0.075/0.9250.773/0.2271.000
rs16899974C/T0.185/0.8150.913/0.0860.385
PRMT1
rs10415880G/A0.407/0.5930.738/0.2620.158
rs975484C/G0.685/0.3150.658/0.3420.197
# Expected allele frequencies were calculated using the allele distributions in the European population of the 1000 Genomes Project (N = 1006) [76]. The p values denote the level of statistical significance for the comparison of absolute frequencies of the minor and major allele in Fisher’s exact test (two-sided). Alleles for which the distribution significantly between the study population and the European reference population are printed in bold type. A p < 0.05 was considered significant. Abbreviations: NOS3, endothelial nitric oxide synthase; DDAH1, dimethylarginine dimethylaminohydrolase-1; DDAH2, dimethylarginine dimethylaminohydrolase-2; ARG1, arginase-1; ARG2, arginase-2; AGXT2, alanine glyoxylate aminotransferase-2; PRMT1, protein-arginine N-methyltransferase-1.

References

  1. Harrod, C.G.; Bendok, B.R.; Batjer, H.H. Prediction of cerebral vasospasm in patients presenting with aneurysmal subarachnoid hemorrhage: A review. Neurosurgery 2005, 56, 633–654. [Google Scholar] [CrossRef]
  2. Miller, C.M.; Palestrant, D.; Schievink, W.I.; Alexander, M.J. Prolonged transcranial Doppler monitoring after aneurysmal subarachnoid hemorrhage fails to adequately predict ischemic risk. Neurocrit. Care 2011, 15, 387–392. [Google Scholar] [CrossRef]
  3. Bauknight, G.C., Jr.; Faraci, F.M.; Heistad, D.D. Endothelium-derived relaxing factor modulates noradrenergic constriction of cerebral arterioles in rabbits. Stroke 1992, 23, 1522–1525. [Google Scholar] [CrossRef] [Green Version]
  4. Pluta, R.M. Delayed cerebral vasospasm and nitric oxide: Review, new hypothesis, and proposed treatment. Pharmacol. Ther. 2005, 105, 23–56. [Google Scholar] [CrossRef] [PubMed]
  5. Li, Q.; Chen, Y.; Li, B.; Luo, C.; Zuo, S.; Liu, X.; Zhang, J.H.; Ruan, H.; Feng, H. Hemoglobin induced NO/cGMP suppression Deteriorate Microcirculation via Pericyte Phenotype Transformation after Subarachnoid Hemorrhage in Rats. Sci. Rep. 2016, 6, 22070. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Hugelshofer, M.; Buzzi, R.M.; Schaer, C.A.; Richter, H.; Akeret, K.; Anagnostakou, V.; Mahmoudi, L.; Vaccani, R.; Vallelian, F.; Deuel, J.W.; et al. Haptoglobin administration into the subarachnoid space prevents hemoglobin-induced cerebral vasospasm. J. Clin. Investig. 2019, 129, 5219–5235. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Terpolilli, N.A.; Feiler, S.; Dienel, A.; Müller, F.; Heumos, N.; Friedrich, B.; Stover, J.; Schöller, K.; Plesnila, N. Nitric oxide inhalation reduces brain damage, prevents mortality, and improves neurological outcome after subarachnoid hemorrhage by resolving early pial microvasospasms. J. Cereb. Blood Flow Metab. 2016, 36, 2096–2107. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Pradilla, G.; Garzon-Muvdi, T.; Ruzevick, J.J.; Bender, M.; Edwards, L.; Momin, E.N.; Thompson, R.C.; Tamargo, R.J. Systemic L-citrulline prevents cerebral vasospasm in haptoglobin 2-2 transgenic mice after subarachnoid hemorrhage. Neurosurgery 2012, 70, 747–756, discussion 756–757. [Google Scholar] [CrossRef]
  9. Böger, R.H.; Ron, E.S. L-Arginine improves vascular function by overcoming deleterious effects of ADMA, a novel cardiovascular risk factor. Altern. Med. Rev. 2005, 10, 14–23. [Google Scholar]
  10. Appel, D.; Seeberger, M.; Schwedhelm, E.; Czorlich, P.; Goetz, A.E.; Böger, R.H.; Hannemann, J. Asymmetric and symmetric dimethylarginines are markers of delayed cerebral ischemia and neurological outcome in patients with subarachnoid hemorrhage. Neurocrit. Care 2018, 29, 84–93. [Google Scholar] [CrossRef]
  11. Bedford, M.T.; Clarke, S.G. Protein arginine methylation in mammals: Who, what, and why. Mol. Cell 2009, 33, 1–13. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Leiper, J.M.; Santa Maria, J.; Chubb, A.; MacAllister, R.J.; Charles, I.G.; Whitley, G.S.; Vallance, P. Identification of two human dimethylarginine dimethylaminohydrolases with distinct tissue distributions and homology with microbial arginine deiminases. Biochem. J. 1999, 343 Pt 1, 209–214. [Google Scholar] [CrossRef]
  13. Lüneburg, N.; Lieb, W.; Zeller, T.; Chen, M.H.; Maas, R.; Carter, A.M.; Xanthakis, V.; Glazer, N.L.; Schwedhelm, E.; Sehsadri, S.; et al. Genome-wide association study of L-arginine and dimethylarginines reveals novel metabolic pathway for symmetric dimethylarginine. Circ. Cardiovasc. Genet. 2014, 7, 864–872. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Ogawa, T.; Kimoto, M.; Sasaoka, K. Dimethylarginine:pyruvate aminotransferase in rats. Purification, properties, and identity with alanine: Glyoxylate aminotransferase 2. J. Biol. Chem. 1990, 265, 20938–20945. [Google Scholar] [PubMed]
  15. Böger, R.; Hannemann, J. Dual role of the L-arginine-ADMA-NO pathway in systemic hypoxic vasodilation and pulmonary hypoxic vasoconstriction. Pulm. Circ. 2020, 10. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Böger, R.H.; Maas, R.; Schulze, F.; Schwedhelm, E. Asymmetric dimethylarginine (ADMA) as a prospective marker of cardiovascular disease and mortality--An update on patient populations with a wide range of cardiovascular risk. Pharmacol. Res. 2009, 60, 481–487. [Google Scholar] [CrossRef]
  17. Böger, R.H.; Sullivan, L.M.; Schwedhelm, E.; Wang, T.J.; Maas, R.; Benjamin, E.J.; Schulze, F.; Xanthakis, V.; Benndorf, R.A.; Vasan, R.S. Plasma asymmetric dimethylarginine and incidence of cardiovascular disease and death in the community. Circulation 2009, 119, 1592–1600. [Google Scholar] [CrossRef]
  18. Lüneburg, N.; von Holten, R.A.; Töpper, R.F.; Schwedhelm, E.; Maas, R.; Böger, R.H. Symmetric dimethylarginine is a marker of detrimental outcome in the acute phase after ischaemic stroke: Role of renal function. Clin. Sci. 2012, 122, 105–111. [Google Scholar] [CrossRef] [Green Version]
  19. Schulze, F.; Carter, A.M.; Schwedhelm, E.; Ajjan, R.; Maas, R.; von Holten, R.A.; Atzler, D.; Grant, P.J.; Böger, R.H. Symmetric dimethylarginine predicts all-cause mortality following ischemic stroke. Atherosclerosis 2010, 208, 518–523. [Google Scholar] [CrossRef]
  20. Jung, C.S.; Oldfield, E.H.; Harvey-White, J.; Espey, M.G.; Zimmermann, M.; Seifert, V.; Pluta, R.M. Association of an endogenous inhibitor of nitric oxide synthase with cerebral vasospasm in patients with aneurysmal subarachnoid hemorrhage. J. Neurosurg. 2007, 107, 945–950. [Google Scholar] [CrossRef]
  21. Li, H.; Wu, W.; Liu, M.; Zhang, X.; Zhang, Q.R.; Ni, L.; Hang, C.H. Increased cerebrospinal fluid concentrations of asymmetric dimethylarginine correlate with adverse clinical outcome in subarachnoid hemorrhage patients. J. Clin. Neurosci. 2014, 21, 1404–1408. [Google Scholar] [CrossRef] [PubMed]
  22. Lindgren, C.; Hultin, M.; Koskinen, L.O.; Lindvall, P.; Borota, L.; Naredi, S. ADMA levels and arginine/ADMA ratios reflect severity of disease and extent of inflammation after subarachnoid hemorrhage. Neurocrit. Care 2014, 21, 91–101. [Google Scholar] [CrossRef] [PubMed]
  23. Jung, C.S.; Iuliano, B.A.; Harvey-White, J.; Espey, M.G.; Oldfield, E.H.; Pluta, R.M. Association between cerebrospinal fluid levels of asymmetric dimethyl-L-arginine, an endogenous inhibitor of endothelial nitric oxide synthase, and cerebral vasospasm in a primate model of subarachnoid hemorrhage. J. Neurosurg. 2004, 101, 836–842. [Google Scholar] [CrossRef] [Green Version]
  24. Ko, N.U.; Rajendran, P.; Kim, H.; Rutkowski, M.; Pawlikowska, L.; Kwok, P.Y.; Higashida, R.T.; Lawton, M.T.; Smith, W.S.; Zaroff, J.G.; et al. Endothelial nitric oxide synthase polymorphism (-786T->C) and increased risk of angiographic vasospasm after aneurysmal subarachnoid hemorrhage. Stroke 2008, 39, 1103–1108. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Starke, R.M.; Kim, G.H.; Komotar, R.J.; Hickman, Z.L.; Black, E.M.; Rosales, M.B.; Kellner, C.P.; Hahn, D.K.; Otten, M.L.; Edwards, J.; et al. Endothelial nitric oxide synthase gene single-nucleotide polymorphism predicts cerebral vasospasm after aneurysmal subarachnoid hemorrhage. J. Cereb. Blood Flow Metab. 2008, 28, 1204–1211. [Google Scholar] [CrossRef]
  26. Alexander, S.; Poloyac, S.; Hoffman, L.; Gallek, M.; Dianxu, R.; Balzer, J.; Kassam, A.; Conley, Y. Endothelial nitric oxide synthase tagging single nucleotide polymorphisms and recovery from aneurysmal subarachnoid hemorrhage. Biol. Res. Nurs. 2009, 11, 42–52. [Google Scholar] [CrossRef]
  27. Abhary, S.; Burdon, K.P.; Kuot, A.; Javadiyan, S.; Whiting, M.J.; Kasmeridis, N.; Petrovsky, N.; Craig, J.E. Sequence variation in DDAH1 and DDAH2 genes is strongly and additively associated with serum ADMA concentrations in individuals with type 2 diabetes. PLoS ONE 2010, 5, e9462. [Google Scholar] [CrossRef] [Green Version]
  28. Caplin, B.; Nitsch, D.; Gill, H.; Hoefield, R.; Blackwell, S.; MacKenzie, D.; Cooper, J.A.; Middleton, R.J.; Talmud, P.J.; Veitch, P.; et al. Circulating methylarginine levels and the decline in renal function in patients with chronic kidney disease are modulated by DDAH1 polymorphisms. Kidney Int. 2010, 77, 459–467. [Google Scholar] [CrossRef] [Green Version]
  29. Thaha, M.; Nilamsari, W.P.; Yusuf, M.; Amin, M. Profile of PRMT-1 gene polymorphism in hemodialysis patients with increased ADMA levels. Acta Med. Indones. 2014, 46, 97–103. [Google Scholar]
  30. Hannemann, J.; Zummack, J.; Hillig, J.; Böger, R. Metabolism of asymmetric dimethylarginine in hypoxia: From bench to bedside. Pulm. Circ. 2020, 10. [Google Scholar] [CrossRef] [Green Version]
  31. Connolly, E.S., Jr.; Rabinstein, A.A.; Carhuapoma, J.R.; Derdeyn, C.P.; Dion, J.; Higashida, R.T.; Hoh, B.L.; KIrkness, C.J.; Naidech, A.M.; Ogilvy, C.S.; et al. Guidelines for the management of aneurysmal subarachnoid hemorrhage: A guideline for healthcare professionals from the American Heart Association/American Stroke Association. Stroke 2012, 43, 1711–1737. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Diringer, M.N.; Bleck, T.P.; Claude Hemphill, J., 3rd; Menon, D.; Shutter, L.; Vespa, P.; Bruder, N.; Connolly, E.S., Jr.; Citerio, G.; Gress, D.; et al. Critical care management of patients following aneurysmal subarachnoid hemorrhage: Recommendations from the Neurocritical Care Society’s Multidisciplinary Consensus Conference. Neurocrit. Care 2011, 15, 211–240. [Google Scholar] [CrossRef] [PubMed]
  33. Drake, C.G.; Hunt, W.E.; Sano, K. Report of World Federation of Neurological Surgeons Committee on a Universal Subarachnoid Hemorrhage Grading Scale. J. Neurosurg. 1988, 68, 985–986. [Google Scholar] [CrossRef]
  34. Gotoh, O.; Tamura, A.; Yasui, N.; Suzuki, A.; Hadeishi, H.; Sano, K. Glasgow Coma Scale in the prediction of outcome after early aneurysm surgery. Neurosurgery 1996, 39, 19–24. [Google Scholar] [CrossRef] [PubMed]
  35. Hunt, W.E.; Hess, R.M. Surgical risk as related to time of intervention in the repair of intracranial aneurysms. J. Neurosurg. 1968, 28, 14–20. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Fisher, C.M.; Kistler, J.P.; Davis, J.M. Relation of cerebral vasospasm to subarachnoid hemorrhage visualized by computerized tomographic scanning. Neurosurgery 1980, 6, 1–9. [Google Scholar] [CrossRef]
  37. Vergouwen, M.D.; Vermeulen, M.; van Gijn, J.; Rinkel, G.J.; Wijdicks, E.F.; Muizelaar, J.P.; Mendelow, A.D.; Juvela, S.; Yonas, H.; Terbrugge, K.G.; et al. Definition of delayed cerebral ischemia after aneurysmal subarachnoid hemorrhage as an outcome event in clinical trials and observational studies: Proposal of a multidisciplinary research group. Stroke 2010, 41, 2391–2395. [Google Scholar] [CrossRef] [Green Version]
  38. Westermaier, T.; Pham, M.; Stetter, C.; Willner, N.; Solymosi, L.; Ernestus, R.I.; Vince, G.H.; Kunze, E. Value of transcranial Doppler, perfusion-CT and neurological evaluation to forecast secondary ischemia after aneurysmal SAH. Neurocrit. Care 2014, 20, 406–412. [Google Scholar] [CrossRef]
  39. Machiela, M.J.; Chanock, S.J. LDlink: A web-based application for exploring population-specific haplotype structure and linking correlated alleles of possible functional variants. Bioinformatics 2015, 31, 3555–3557. [Google Scholar] [CrossRef]
  40. Schwedhelm, E.; Tan-Andresen, J.; Maas, R.; Riederer, U.; Schulze, F.; Böger, R.H. Liquid chromatography-tandem mass spectrometry method for the analysis of asymmetric dimethylarginine in human plasma. Clin. Chem. 2005, 51, 1268–1271. [Google Scholar] [CrossRef] [Green Version]
  41. Böger, R.H.; Sydow, K.; Borlak, J.; Thum, T.; Lenzen, H.; Schubert, B.; Tsikas, D.; Bode-Böger, S.M. LDL cholesterol upregulates synthesis of asymmetrical dimethylarginine in human endothelial cells: Involvement of S-adenosylmethionine-dependent methyltransferases. Circ. Res. 2000, 87, 99–105. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Pikula, A.; Böger, R.H.; Beiser, A.S.; Maas, R.; DeCarli, C.; Schwedhelm, E.; Himali, J.J.; Schulze, F.; Au, R. Kelly-Hayes, M.; et al. Association of plasma ADMA levels with MRI markers of vascular brain injury: Framingham offspring study. Stroke 2009, 40, 2959–2964. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Rodionov, R.N.; Murry, D.J.; Vaulman, S.F.; Stevens, J.W.; Lentz, S.R. Human alanine-glyoxylate aminotransferase 2 lowers asymmetric dimethylarginine and protects from inhibition of nitric oxide production. J. Biol. Chem. 2010, 285, 5385–5391. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Arlt, S.; Schulze, F.; Eichenlaub, M.; Maas, R.; Lehmbeck, J.T.; Schwedhelm, E.; Jahn, H.; Böger, R.H. Asymmetrical dimethylarginine is increased in plasma and decreased in cerebrospinal fluid of patients with Alzheimer’s disease. Dement Geriatr. Cogn. Disord. 2008, 26, 58–64. [Google Scholar] [CrossRef]
  45. Xu, L.; Wang, B.; Kaur, K.; Kho, M.F.; Cooke, J.P.; Giffard, R.G. NOx and ADMA changes with focal ischemia, amelioration with the chaperonin GroEL. Neurosci. Lett. 2007, 418, 201–204. [Google Scholar] [CrossRef] [Green Version]
  46. Brouns, R.; Marescau, B.; Possemiers, I.; Sheorajpanday, R.; De Deyn, P.P. Dimethylarginine levels in cerebrospinal fluid of hyperacute ischemic stroke patients are associated with stroke severity. Neurochem. Res. 2009, 34, 1642–1649. [Google Scholar] [CrossRef]
  47. Closs, E.I.; Basha, F.Z.; Habermeier, A.; Förstermann, U. Interference of L-arginine analogues with L-arginine transport mediated by the y+ carrier hCAT-2B. Nitric Oxide 1997, 1, 65–73. [Google Scholar] [CrossRef]
  48. McBride, A.E.; Silver, P.A. State of the arg: Protein methylation at arginine comes of age. Cell 2001, 106, 5–8. [Google Scholar] [CrossRef] [Green Version]
  49. Malik, R.; Rannikmae, K.; Traylor, M.; Georgakis, M.K.; Sargurupremraj, M.; Markus, H.S.; Hopewell, J.C.; Debette, S.; Sudlow, C.L.M.; Dichgans, M. Genome-wide meta-analysis identifies 3 novel loci associated with stroke. Ann. Neurol. 2018, 84, 934–939. [Google Scholar] [CrossRef]
  50. Morris, C.M.; Ballard, C.G.; Allan, L.; Rowan, E.; Stephens, S.; Firbank, M.; Ford, G.A.; Kenny, R.A.; O’Brien, J.T.; Kalaria, R.N. NOS3 gene rs1799983 polymorphism and incident dementia in elderly stroke survivors. Neurobiol. Aging 2011, 32, e551–e556. [Google Scholar] [CrossRef]
  51. Wang, M.; Jiang, X.; Wu, W.; Zhang, D. Association of G894T polymorphism in endothelial nitric oxide synthase gene with the risk of ischemic stroke: A meta-analysis. Biomed. Rep. 2013, 1, 144–150. [Google Scholar] [CrossRef] [PubMed]
  52. Zhang, X.; Lynch, A.I.; Davis, B.R.; Ford, C.E.; Boerwinkle, E.; Eckfeldt, J.H.; Leiendecker-Foster, C.; Arnett, D.K. Pharmacogenetic association of NOS3 variants with cardiovascular disease in patients with hypertension: The GENHAT study. PLoS ONE 2012, 7, e34217. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Hendrix, P.; Foreman, P.M.; Harrigan, M.R.; Fisher, W.S., 3rd; Vyas, N.A.; Lipsky, R.H.; Lin, M.; Walters, B.C.; Tubbs, R.S.; Shoja, M.M.; et al. Endothelial nitric oxide synthase polymorphism is associated with delayed cerebral ischemia following aneurysmal subarachnoid hemorrhage. World Neurosurg. 2017, 101, 514–519. [Google Scholar] [CrossRef] [PubMed]
  54. Hendrix, P.; Foreman, P.M.; Harrigan, M.R.; Fisher, W.S.R.; Vyas, N.A.; Lipsky, R.H.; Lin, M.; Walters, B.C.; Tubbs, R.S.; Shoja, M.M.; et al. The role of endothelial nitric oxide synthase -786 T/C polymorphism in cardiac instability following aneurysmal subarachnoid hemorrhage. Nitric Oxide 2017, 71, 52–56. [Google Scholar] [CrossRef] [PubMed]
  55. Yamada, Y.; Kato, K.; Oguri, M.; Yoshida, T.; Yokoi, K.; Watanabe, S.; Metoki, N.; Yoshida, H.; Satoh, K.; Ichihara, S.; et al. Association of genetic variants with atherothrombotic cerebral infarction in Japanese individuals with metabolic syndrome. Int. J. Mol. Med. 2008, 21, 801–808. [Google Scholar] [CrossRef] [Green Version]
  56. Kardia, S.L.; Greene, M.T.; Boerwinkle, E.; Turner, S.T.; Kullo, I.J. Investigating the complex genetic architecture of ankle-brachial index, a measure of peripheral arterial disease, in non-hispanic whites. BMC Med. Genom. 2008, 1, 16. [Google Scholar] [CrossRef] [Green Version]
  57. Kullo, I.J.; Greene, M.T.; Boerwinkle, E.; Chu, J.; Turner, S.T.; Kardia, S.L. Association of polymorphisms in NOS3 with the ankle-brachial index in hypertensive adults. Atherosclerosis 2008, 196, 905–912. [Google Scholar] [CrossRef] [Green Version]
  58. Li, X.; Lin, Y.; Zhang, R. Associations between endothelial nitric oxide synthase gene polymorphisms and the risk of coronary artery disease: A systematic review and meta-analysis of 132 case-control studies. Eur. J. Prev. Cardiol. 2019, 26, 160–170. [Google Scholar] [CrossRef]
  59. Seelenfreund, D.; Lobos, S.; Quesada, A.; Saavedra, J.; Wolff, C.; López-Stewart, G.; Araya, A.; Durruty, P. Association of the intronic polymorphism rs891512 (G24943A) of the endothelial nitric oxide synthase gene with hypertension in Chilean type 2 diabetes patients. Diabetes Res. Clin. Pract. 2012, 96, e47–e49. [Google Scholar] [CrossRef]
  60. Vimaleswaran, K.S.; Franks, P.W.; Barroso, I.; Brage, S.; Ekelund, U.; Wareham, N.J.; Loos, R.J.F. Habitual energy expenditure modifies the association between NOS3 gene polymorphisms and blood pressure. Am. J. Hypertens 2008, 21, 297–302. [Google Scholar] [CrossRef] [Green Version]
  61. Lu, T.M.; Lin, S.J.; Lin, M.W.; Hsu, C.P.; Chung, M.Y. The association of dimethylarginine dimethylaminohydrolase 1 gene polymorphism with type 2 diabetes: A cohort study. Cardiovasc. Diabetol. 2011, 10, 16. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Schnabel, R.; Larson, M.G.; Dupuis, J.; Lunetta, K.L.; Lipinska, I.; Meigs, J.B.; Yin, X.; Rong, Y.; Vita, J.A.; Newton-Cheh, C.; et al. Relations of inflammatory biomarkers and common genetic variants with arterial stiffness and wave reflection. Hypertension 2008, 51, 1651–1657. [Google Scholar] [CrossRef] [PubMed]
  63. Trittmann, J.K.; Gastier-Foster, J.M.; Zmuda, E.J.; Frick, J.; Rogers, L.K.; Vieland, V.J.; Chicoine, L.G.; Nelin, L.D. A single nucleotide polymorphism in the dimethylarginine dimethylaminohydrolase gene is associated with lower risk of pulmonary hypertension in bronchopulmonary dysplasia. Acta Paediatr. 2016, 105, e170–e175. [Google Scholar] [CrossRef] [Green Version]
  64. Maas, R.; Erdmann, J.; Lüneburg, N.; Stritzke, J.; Schwedhelm, E.; Meisinger, C.; Peters, A.; Weil, J.; Schunkert, H.; Böger, R.H.; et al. Polymorphisms in the promoter region of the dimethylarginine dimethylaminohydrolase 2 gene are associated with prevalence of hypertension. Pharmacol Res. 2009, 60, 488–493. [Google Scholar] [CrossRef]
  65. Marra, M.; Marchegiani, F.; Ceriello, A.; Sirolla, C.; Boemi, M.; Franceschi, C.; Spazzafumo, L.; Testa, I.; Bonfigli, A.R.; Cucchi, M.; et al. Chronic renal impairment and DDAH2-1151 a/c polymorphism determine ADMA levels in type 2 diabetic subjects. Nephrol. Dial. Transpl. 2013, 28, 964–971. [Google Scholar] [CrossRef] [Green Version]
  66. Perez-Hernandez, N.; Vargas-Alarcon, G.; Arellano-Zapoteco, R.; Martinez-Rodriguez, N.; Fragoso, J.M.; Aptilon-Duque, G.; Posadas-Sanchez, R.; Posadas-Romero, C.; Juarez-Cedillo, T.; Dominguez-Lopez, M.L.; et al. Protective role of DDAH2 (rs805304) gene polymorphism in patients with myocardial infarction. Exp. Mol. Pathol. 2014, 97, 393–398. [Google Scholar] [CrossRef] [PubMed]
  67. Seo, H.A.; Kim, S.W.; Jeon, E.J.; Jeong, J.Y.; Moon, S.S.; Lee, W.K.; Kim, J.G.; Lee, I.K.; Park, K.G. Association of the DDAH2 gene polymorphism with type 2 diabetes and hypertension. Diabetes Res. Clin. Pract. 2012, 98, 125–131. [Google Scholar] [CrossRef] [PubMed]
  68. Lacchini, R.; Muniz, J.J.; Nobre, Y.T.; Cologna, A.J.; Martins, A.C.; Tanus-Santos, J.E. Relationship between arginase 1 and arginase 2 levels and genetic polymorphisms with erectile dysfunction. Nitric Oxide 2015, 51, 36–42. [Google Scholar] [CrossRef]
  69. Salam, M.T.; Bastain, T.M.; Rappaport, E.B.; Islam, T.; Berhane, K.; Gauderman, W.J.; Gilliland, F.D. Genetic variations in nitric oxide synthase and arginase influence exhaled nitric oxide levels in children. Allergy 2011, 66, 412–419. [Google Scholar] [CrossRef] [Green Version]
  70. Vonk, J.M.; Postma, D.S.; Maarsingh, H.; Bruinenberg, M.; Koppelman, G.H.; Meurs, H. Arginase 1 and arginase 2 variations associate with asthma, asthma severity and beta2 agonist and steroid response. Pharm. Genom. 2010, 20, 179–186. [Google Scholar] [CrossRef]
  71. Amir, M.; Hassanein, S.I.; Abdel Rahman, M.F.; Gad, M.Z. AGXT2 and DDAH-1 genetic variants are highly correlated with serum ADMA and SDMA levels and with incidence of coronary artery disease in Egyptians. Mol. Biol. Rep. 2018, 45, 2411–2419. [Google Scholar] [CrossRef] [PubMed]
  72. Caplin, B.; Wang, Z.; Slaviero, A.; Tomlinson, J.; Dowsett, L.; Delahaye, M.; Salama, A.; Wheeler, D.C.; Leiper, J. Alanine-glyoxylate aminotransferase-2 metabolizes endogenous methylarginines, regulates NO, and controls blood pressure. Arter. Thromb. Vasc. Biol. 2012, 32, 2892–2900. [Google Scholar] [CrossRef] [Green Version]
  73. Kittel, A.; Müller, F.; König, J.; Mieth, M.; Sticht, H.; Zolk, O.; Kralj, A.; Heinrich, M.R.; Fromm, M.F.; Maas, R. Alanine-glyoxylate aminotransferase 2 (AGXT2) polymorphisms have considerable impact on methylarginine and beta-aminoisobutyrate metabolism in healthy volunteers. PLoS ONE 2014, 9, e88544. [Google Scholar] [CrossRef] [PubMed]
  74. Seppälä, I.; Kleber, M.E.; Bevan, S.; Lyytikäinen, L.P.; Oksala, N.; Hernesniemi, J.A.; Mäkelä, K.M.; Rothwell, P.M.; Sudlow, C.; Dichgans, M.; et al. Associations of functional alanine-glyoxylate aminotransferase 2 gene variants with atrial fibrillation and ischemic stroke. Sci. Rep. 2016, 6, 23207. [Google Scholar] [CrossRef] [PubMed]
  75. Lee, K.H.; Tsai, W.J.; Chen, Y.W.; Yang, W.C.; Lee, C.Y.; Ou, S.M.; Chen, Y.T.; Chien, C.C.; Lee, P.C.; Chung, M.Y.; et al. Genotype polymorphisms of genes regulating nitric oxide synthesis determine long-term arteriovenous fistula patency in male hemodialysis patients. Ren. Fail. 2016, 38, 228–237. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  76. Auton, A.; Brooks, L.D.; Durbin, R.M.; Garrison, E.P.; Kang, H.M.; Korbel, J.O.; Marchini, J.L.; McCarthy, S.; McVean, G.A.; Abecasis, G.R. A global reference for human genetic variation. Nature 2015, 526, 68–74. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 2. Concentrations of ADMA, SDMA, and L-arginine in (a) plasma and (b) CSF at the time of admission. A total of 51 patients with SAH were categorized into those who developed delayed cerebral ischemia (N = 18) or not (N = 33). Data are mean ± S.E.M. * p < 0.05, ** p < 0.01 vs. no DCI group. Abbreviations: ADMA, asymmetric dimethylarginine; CSF, cerebrospinal fluid; DCI, delayed cerebral ischemia; SDMA, symmetric dimethylarginine; S.E.M., standard error of the mean
Figure 2. Concentrations of ADMA, SDMA, and L-arginine in (a) plasma and (b) CSF at the time of admission. A total of 51 patients with SAH were categorized into those who developed delayed cerebral ischemia (N = 18) or not (N = 33). Data are mean ± S.E.M. * p < 0.05, ** p < 0.01 vs. no DCI group. Abbreviations: ADMA, asymmetric dimethylarginine; CSF, cerebrospinal fluid; DCI, delayed cerebral ischemia; SDMA, symmetric dimethylarginine; S.E.M., standard error of the mean
Jcm 09 03900 g002
Figure 3. Concentrations of ADMA in plasma or CSF in relation to genotypes of the DDAH1 SNP rs233112 (a,b), the DDAH1 SNP rs480414 (c), and the DDAH1 SNP rs124321 (d). Data are mean ± S.E.M. Abbreviations: ADMA, asymmetric dimethylarginine; CSF, cerebrospinal fluid; n.s.; not significant; S.E.M., standard error of the mean; SNP, single nucleotide polymorphism.
Figure 3. Concentrations of ADMA in plasma or CSF in relation to genotypes of the DDAH1 SNP rs233112 (a,b), the DDAH1 SNP rs480414 (c), and the DDAH1 SNP rs124321 (d). Data are mean ± S.E.M. Abbreviations: ADMA, asymmetric dimethylarginine; CSF, cerebrospinal fluid; n.s.; not significant; S.E.M., standard error of the mean; SNP, single nucleotide polymorphism.
Jcm 09 03900 g003
Figure 4. Plasma and CSF concentrations of SDMA (a,c) and ADMA (b,d) in relation to genotypes of the AGXT2 rs168999974 SNP. Data are mean ± S.E.M. Abbreviations: ADMA, asymmetric dimethylarginine; CSF, cerebrospinal fluid; SDMA, symmetric dimethylarginine; S.E.M., standard error of the mean; SNP, single nucleotide polymorphism.
Figure 4. Plasma and CSF concentrations of SDMA (a,c) and ADMA (b,d) in relation to genotypes of the AGXT2 rs168999974 SNP. Data are mean ± S.E.M. Abbreviations: ADMA, asymmetric dimethylarginine; CSF, cerebrospinal fluid; SDMA, symmetric dimethylarginine; S.E.M., standard error of the mean; SNP, single nucleotide polymorphism.
Jcm 09 03900 g004
Figure 5. Frequency distribution of genotypes for (a) DDAH1 SNP rs233112, (b) DDAH1 SNP rs480414, and (c) DDAH1 SNP rs1241321 in SAH patients with or without DCI. For all three SNPs, the minor allele was more frequently observed in patients with DCI than in those without (p values as indicated in the graphs). Abbreviations: DDAH1, dimethylarginine dimethylaminohydrolase 1; DCI, delayed cerebral ischemia; SAH, subarachnoidal hemorrhage; SNP, single nucleotide polymorphism.
Figure 5. Frequency distribution of genotypes for (a) DDAH1 SNP rs233112, (b) DDAH1 SNP rs480414, and (c) DDAH1 SNP rs1241321 in SAH patients with or without DCI. For all three SNPs, the minor allele was more frequently observed in patients with DCI than in those without (p values as indicated in the graphs). Abbreviations: DDAH1, dimethylarginine dimethylaminohydrolase 1; DCI, delayed cerebral ischemia; SAH, subarachnoidal hemorrhage; SNP, single nucleotide polymorphism.
Jcm 09 03900 g005
Table 1. Baseline characteristics of the patients.
Table 1. Baseline characteristics of the patients.
AllPatients with DCIPatients without DCIp
Demographics
Number of patients511833
Age, years (mean ± SD)54.8 ± 13.158.6 ± 14.152.6 ± 12.20.140
Sex, male/female (%)13/38 (25.5/75.5)2/16 (11.2/88.8)11/22 (33.3/66.7)0.102
BMI, kg/m2 (mean ± SD)25.9 ± 3.926.8 ± 4.925.4 ± 3.20.333
Risk factors, n (%)
Hypertension21 (41.2)9 (50.0)12 (36.4)0.385
Diabetes mellitus3 (5.9)1 (5.6)2 (6.1)1.000
Renal failure1 (2.0)1 (5.6)0 (0.0)0.353
Smoker, current + ex11 (21.6)4 (22.2)7 (21.2)1.000
Long-term medication, n (%)
Angiotensin II blockers6 (11.8)2 (11.2)4 (12.1)0.586
Beta blockers8 (15.7)3 (16.7)5 (15.2)0.591
Diuretics3 (5.9)1 (5.6)2 (6.1)0.718
Calcium channel blockers5 (9.8)2 (11.2)3 (9.1)0.585
Statins6 (11.8)2 (11.2)4 (12.1)0.646
NSAR3 (5.9)1 (5.6)2 (6.1)0.718
Platelet inhibitors5 (9.8)2 (11.2)3 (9.1)0.585
Oral anticoagulants1 (2.0)0 (0.0)1 (3.0)0.647
Bleeding location, n (%)
Anterior cerebral circulation35 (64.8)12 (66.7)23 (69.7)0.532
Posterior cerebral circulation12 (22.2)4 (22.2)8 (24.2)0.579
Aneurysm treatment, n (%)
Clipping15 (29.4)6 (33.3)9 (27.3)0.442
Coiling28 (54.9)9 (50.0)18 (54.5)0.492
Non-interventional8 (15.7)2 (11.1)6 (18.2)0.409
Nimodipin treatment47 (92.2)17 (94.4)30 (90.9)0.557
Abbreviations: BMI, body mass index; DCI, delayed cerebral ischemia; NSAR, non-steroidal antirheumatic drugs; SD, standard deviation. p values are given for differences between patients with and without DCI.
Table 2. Clinical scores and outcome of SAH patients with or without DCI.
Table 2. Clinical scores and outcome of SAH patients with or without DCI.
AllPatients with DCIPatients without DCIp
Number of patients511833
Scores at admission (mean ± SD)
GCS8.6 ± 5.610.9 ± 4.67.2 ± 5.70.021
Hunt and Hess3.1 ± 1.42.9 ± 1.23.2 ± 1.50.622
Fisher Score3.3 ± 0.73.2 ± 0.83.4 ± 0.70.416
WFNS3.4 ± 1.72.9 ± 1.53.6 ± 1.80.134
Outcome, n (%)
Vasospasm22 (43.1)10 (55.6)12 (36.4)0.242
In-hospital lethality7 (13.7)3 (16.7)4 (12.1)0.686
Duration of ICU treatment (d)21.0 ± 11.522.7 ± 11.220.1 ± 11.70.449
Mechanical ventilation (h)237.5 ± 272.5313.3 ± 251.7196.1 ± 278.20.144
Re-intervention rate5 (9.8)4 (22.2)1 (3.0)0.047
GOS at discharge3.4 ± 1.62.8 ± 1.63.6 ± 1.50.056
GOS three months later3.6 ± 1.13.1 ± 1.14.0 ± 1.00.071
Abbreviations: GCS, Glasgow Coma Scale; WFNS, World Federation of Neurological Surgeons; GOS, Glasgow Outcome Scale. p values are given for differences between patients with and without DCI.
Table 3. Association of the combined DDAH1 haplotype with outcome parameters.
Table 3. Association of the combined DDAH1 haplotype with outcome parameters.
DDAH1 HaplotypeRelative Risk95% CIp
(3× homozygous major versus 0× homozygous minor)
DCI1.1030.925–1.3150.273
Vasospasm1.4831.013–2.1700.039
Death1.0000.898–1.1131.000
GOS < 4 at discharge1.3641.105–1.6830.003
GOS < 4 after 3 months1.4751.231–1.766<0.0001
(3× homozygous major versus ≥ 1× homozygous minor)
DCI1.3161.072–1.6140.011
Vasospasm0.9830.774–1.2471.000
Death0.9290.460–1.8741.000
GOS < 4 at discharge2.2731.681–3.072<0.0001
GOS < 4 after 3 months1.2991.110–1.5190.001
Details of the haplotype construction are given in Appendix A. Abbreviations: DCI, delayed cerebral ischemia; GOS, Glasgow outcome scale. p values presented were calculated by Chi2 test; significant p values are printed in bold type. Abbreviation: 95% CI, 95% confidence interval.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Hannemann, J.; Appel, D.; Seeberger-Steinmeister, M.; Brüning, T.; Zummack, J.; Böger, R. Sequence Variation in the DDAH1 Gene Predisposes for Delayed Cerebral Ischemia in Subarachnoidal Hemorrhage. J. Clin. Med. 2020, 9, 3900. https://doi.org/10.3390/jcm9123900

AMA Style

Hannemann J, Appel D, Seeberger-Steinmeister M, Brüning T, Zummack J, Böger R. Sequence Variation in the DDAH1 Gene Predisposes for Delayed Cerebral Ischemia in Subarachnoidal Hemorrhage. Journal of Clinical Medicine. 2020; 9(12):3900. https://doi.org/10.3390/jcm9123900

Chicago/Turabian Style

Hannemann, Juliane, Daniel Appel, Miriam Seeberger-Steinmeister, Tabea Brüning, Julia Zummack, and Rainer Böger. 2020. "Sequence Variation in the DDAH1 Gene Predisposes for Delayed Cerebral Ischemia in Subarachnoidal Hemorrhage" Journal of Clinical Medicine 9, no. 12: 3900. https://doi.org/10.3390/jcm9123900

APA Style

Hannemann, J., Appel, D., Seeberger-Steinmeister, M., Brüning, T., Zummack, J., & Böger, R. (2020). Sequence Variation in the DDAH1 Gene Predisposes for Delayed Cerebral Ischemia in Subarachnoidal Hemorrhage. Journal of Clinical Medicine, 9(12), 3900. https://doi.org/10.3390/jcm9123900

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