Association of endothelin genetic variants and hospitalized infection complications in end-stage renal disease (ESRD) patients

Background Infection is the second most common cause of mortality for patients with end-stage renal disease (ESRD), accompanying with immune dysfunction. Endothelin (EDN) is known to be related to inflammation; however, it is unknown whether genetic variants of the EDN gene family are associated with increased risk of hospitalized infection events. Methods Nineteen tagging single-nucleotide polymorphisms (tSNPs) of the EDN gene family were selected for genotyping a cohort of 190 ESRD patients. Patient demographics were recorded, the subtypes of infection events were identified, and association analysis between the EDN genetic variants and hospitalized infection events was performed. Results In this study, 106 patients were hospitalized for infection events. The leading events were pneumonia, bacteremia, and cellulitis. The minor allele of rs260741, rs197173, and rs926632 SNPs of EDN3 were found to be associated with reduced risk of hospitalized bacteremia events. Conclusions The minor allele of rs260741, rs197173, and rs926632 in EDN3 were associated with reduced risk of hospitalized bacteremia events in ESRD patients. Electronic supplementary material The online version of this article (10.1186/s12882-019-1349-3) contains supplementary material, which is available to authorized users.


Background
Infection is known as a common cause of morbidity and mortality in patients with end-stage renal disease (ESRD), accounting for 20% of total deaths in these patients [1]. More than half of patients hospitalized for infection events developed an unfavorable outcome, including prolonged hospitalization, intensive care unit (ICU) stay, or even death [2,3]. Impaired innate and adaptive immunity are associated with increased risk of infection in patients with renal failure [4]. A handful of cytokines was dysregulated in ESRD patients; for example, the production of interleukin (IL)-1, tumor necrosis factor (TNF)-α, and IL-6 was increased, whereas the bioavailability of IL-2 was reduced [5].
Genetic variants of immune-related genes have been shown to reflect risk of infection. Ferwerda et al. reported that the genetic variants of Toll-like receptor 4, which underlies the differential production of cytokines may affect the immune system, which in turn influenced the susceptibility of infection events and the risk of gram-negative bacterial infection [6]. Moreover, another study demonstrated that an IL-9 variant may influence the susceptibility of respiratory syncytial virus infection [7].
Endothelin-1 (ET-1) is the major isoform of the endothelin (EDN) family and is a potent vasoconstrictor. It is associated with proinflammatory cytokine production, fibrosis, and angiogenesis [8]. The release of lipopolysaccharide from bacteria impairs the integrity of the endothelial cell, resulting in endothelial cell injury and cytokine release [9]. The plasma ET-1 level is increased during sepsis and is correlated with the severity of sepsis [10]. ET-1 is known to increase reactive oxygen species (ROS) production, and associates with activation of nuclear factor-kappaB and inflammatory cytokines such as TNF-α, IL-1, and IL-6 [11]. In addition, Lin et al. reported that ET-1 is able to increase cyclooxygenase-2 expression and prostaglandin E2 release [12]. These findings suggest that ET-1 is involved in the inflammatory reactions and the severity of sepsis. Regarding ET-3, Sato et al. demonstrated that low level of ET-3 reduces inflammatory responses [13]. However, it is still unclear whether EDN genetic variants are associated with hospitalized infection events in ESRD patients. In this study, genetic association study was applied to investigate the correlations between EDN genetic variants and the hospitalized infection events in ESRD patients.

Study population
Patients who received dialysis for more than 3 months at Taipei Medical University Hospital between September 2013 and June 2014 were enrolled. Demographic and laboratory data of these patients was shown in Table 1. The erythropoietin resistance index (ERI) was determined to characterize patients' response to erythropoietin. ERI was calculated by dividing the weekly body-weight-adjusted epoetin dose by the hemoglobin concentration. Kt/V was used to measure dialysis adequacy. Outpatient and discharge medical records were used to determine the cause of ESRD. We prospectively followed up with the patients who were hospitalized for infection events. This study was approved by the Institutional Review Board of Taipei Medical University (Approval No. 201309026). A written informed consent form was obtained from all patients.

Infection events
We defined hospitalized "infection events" as bacteremia, pneumonia, cellulitis, urinary tract infection (UTI), peritonitis, and intra-abdominal infection (IAI). Bacteremia was confirmed by a positive blood culture result, and contamination was excluded by an infection specialist. Pneumonia was defined as the presence of clinical respiratory symptoms and the findings of increased infiltration on chest radiography. Cellulitis was determined as inflammation of the skin and subcutaneous tissues. UTI was defined as the presence of clinical symptoms and the detection of a pathogen in the urine. Peritonitis was diagnosed through clinical symptom examination and peritoneal fluid analysis. IAI was confirmed by the finding of intramural inflammation of the gastrointestinal (GI) tract without anatomic disruption. These diagnoses were ascertained by physicians and recorded in the discharge medical records.

Associations between EDN genetic variants and hospitalized infection events
The association analysis was performed for detecting the relationship of EDN genetic variants and hospitalized infection events. We haven't observed a significant association between the frequencies of both EDN1 and EDN2 genotypes and clinical outcomes (Additional file 1: Table  S2 and S3), but the rs882345, rs926632 and rs3026575 The P value was adjusted for age, sex, diabetes, hemoglobin, albumin, and the cause of ESRD. The P values of < 0.05 are shown in bold of EDN3 showed a modest association with hospitalized infection events (Table 2).

Associations between EDN genetic variants and hospitalized bacteremia events
The rs260741, rs6064764, rs197173, and rs926632 of EDN3 genotypes showed a significant association with hospitalized bacteremia events (Table 3). Patients carrying the minor allele of EDN3 rs260741 (AA/AG vs. GG), rs197173 (TT/TG vs. GG), and rs926632 (CC/CT vs. TT) had a lower risk of hospitalized bacteremia events (genotype model, dominant model, and allelic model P value < 0.05). The minor allele of rs6064764 (CC/CT vs. TT) was associated with increased risk of bacteremia according to the dominant model and allelic model P values (Table 3). No association of EDN1 or EDN2 SNPs with hospitalized bacteremia events were detected (Additional file 1: Table S4 and S5).

Haplotypes analysis of EDN3 with infection and bacteremia events
Next, we further validated the effects of haplotypes of EDN3 gene by pairwise linkage disequilibrium (LD) (Fig. 2). The haplotype frequency of EDN3 rs926632/ rs882345 variants among patients with hospitalized infection and bacteremia events were shown in Tables 4 and 5. The EDN3 haplotype had no significant association with the risk of infection or bacteremia events. To test the possible functional roles of the polymorphisms in EDN3, we queried expression quantitative trait loci (eQTL) of EDN3 via Genotype-Tissue Expression (GTEx) database in different types of tissues [14]. Of note, the low expression level of EDN3 was found in immune cells (Additional file 1: Figure S4). In addition, we analyzed the SNPs to understand the functional roles through using HaploReg V4.1. Interestingly, results indicated that the non-coding SNP rs260741 and rs6064764 of EDN3 were potentially related to the regulation of T cells activation by epigenetic modifications [15].

Discussion
Dialysis patients are exposed to high risk of infection, and related mortality [16,17]. Ishigami et al. reported that a low estimated glomerular filtration rate and high albuminuria were associated with increased risk of hospitalization for infection, including bloodstream infections, pneumonia, UTI, and cellulitis [18]. Among these infections, bacteremia was the most critical and responsible for three-quarters of infection-related mortality [19]. The risk factors for infection in dialysis patients include age, immunosuppressive therapy, poor hygiene, and low performance status [17,[20][21][22]. Age, albumin level, and category of infection have been shown to associate with poor outcomes [23]. Old age and lower albumin level were the risk factors for infection events. Indeed, both innate and adaptive immunity are known to change with age, which may result in a persistent low-grade inflammation and tissue damage [24]. Furthermore, T-cell repertoire (TCR) complexity was revealed to predict the EPO responsiveness. Thus, TCR repertoire diversity may indicate the immune responses to infection [25]. Malnutrition, inflammation, and atherosclerosis syndrome, characterized by patients with renal failure, are predictive of poor outcomes [26].   Importantly, malnutrition has also been reported to associate with chronic inflammation, followed by increased risk of infection [27]. Correa et al. reported that EDN signaling pathway is involved in the pathogenesis of mycobacteria tuberculosis infection, and that EDN receptor A or B signaling is critical for the host responses to the infection [28]. Wilson et al. indicated the functional roles of EDN genes in intestinal hypoperfusion during bacteremia [29]. Additionally, the precursor peptide ET-1 of EDN1 gene is correlated with the severity of pneumonia, ICU admission, and mortality [30]. Also, ET-1 has been shown to stimulate ROS production [9] and correlate with the risk of inflammatory diseases such as atherosclerosis. Our previous study showed genetic variants of EDN1 gene are associated with an increased risk of hospitalization for a cardiovascular event [31]. Pittet et al. demonstrated that the serum level of ET-1 was strongly correlated with reduced cardiac output in sepsis patients [32]. In this study, we found the correlation between EDN genetic variants and hospitalized infection events..
A previous study on EDN3 and EDN receptor type B (EDNRB) knockout mouse model showed that EDN3 and EDNRB were associated with severe enterocolitis [33]. The ET-3 peptide of EDN3 gene was reported to evoke an attenuated inflammation by inducing the EDN B2 receptor and nitric oxide production [13]. Another study also proposed that ET-3 reduces platelet-activating factor (PAF)-induced inflammation by directly binding to PAF [34]. Interestingly, the Human Protein Atlas showed that ET-3 is highly expressed in the GI tract and endocrine tissues [35]. During the development of enteric nervous system, ET-3 and its receptor, EDNRB, orchestrate the signaling cascades for enteric ganglion formation [36][37][38]. The enteric nervous system, composing of neurons and glial cells, have been shown to interact with the outer (microbiota) and inner environments (immune cells) [39]; the enteric microbiota is known to interact with the immune system [40]. These studies indicated that ET-3 may affect the inflammation process through the modification of enteric nervous system and interaction with the microbiota. The results provided indirect evidence for supporting our findings in this study.
This study has limitations. First, our small sample size may affect the statistical power in the analysis. Second, the gene expression and protein level of EDN were not measured from patients' samples. Further validation studies are needed in order to evaluate the associations of genetic variants with gene/protein expression. Third, a replication study should be conducted in the second population.

Conclusions
The minor allele of rs260741, rs197173, and rs926632 in EDN3 was associated with a reduced risk of hospitalized bacteremia events in patients with end-stage renal disease (ESRD).