Investigation of the variants at the binding site of inflammatory transcription factor NF-κB in patients with end-stage renal disease
BMC Nephrology volume 20, Article number: 300 (2019)
A chronic inflammatory state is a prominent feature in patients with end-stage renal disease (ESRD). Nuclear factor-kappa B (NF-κB) is a transcription factor that regulates the expression of genes involved in inflammation. Some genetic studies have demonstrated that the NF-κB genetic mutation could cause kidney injury and kidney disease progression. However, the association of a gene polymorphism in the transcription factor binding site of NF-κB with kidney disease is not clear.
We used the Taiwan Biobank database, the University of California, Santa Cruz, reference genome, and a chromatin immunoprecipitation sequencing database to find single nucleotide polymorphisms (SNPs) at potential binding sites of NF-κB. In addition, we performed a case–control study and genotyped 847 patients with ESRD and 846 healthy controls at Tri-Service General Hospital from 2015 to 2016. Furthermore, we used the ChIP assay to identify the binding activity of different genotypes and used Luciferase reporter assay to examine the function of the rs9395890 polymorphism.
The results of biometric screening in the databases revealed 15 SNPs with the potential binding site of NF-κB. Genotype distributions of rs9395890 were significantly different in ESRD cases and healthy controls (P = 0.049). The ChIP assay revealed an approximately 1.49-fold enrichment of NF-κB of the variant type TT when compared to that of the wild-type GG in rs9395890 (P = 0.027; TT = 3.20 ± 0.16, GT = 2.81 ± 0.20, GG = 1.71 ± 0.18). The luciferase reporter assay showed that the NF-κB binding site activity in T allele was slightly higher than that in G allele, though it is not significant.
Our findings indicate that rs9395890 is associated with susceptibility to ESRD in Taiwan population.
According to the United States Renal Data System annual report, the incidence and prevalence of end-stage renal disease (ESRD) in Taiwan are among the highest in the world [1, 2]. Taiwan has a 9.8% prevalence of chronic kidney disease (CKD) . The incidence of ESRD in Taiwanese individuals is 450 per million people [2, 4]. CKD-related and ESRD-related costs in Taiwan are US$25,576 per patient-year . CKD development and progression to ESRD involve complex interactions between multiple genetic and environmental factors . Chronic inflammation is an important component of CKD and ESRD. It has a unique role in their pathophysiology and contributes to cardiovascular and all-cause mortality, as well as the development of protein-energy wasting [7, 8]. Genetic factors are important risk factors in the pathogenesis of CKD. The heritability of ESRD is 31.1% in the Taiwanese population .
Nuclear factor-kappa B (NF-κB) is an important transcription factor in inflammation and promotes the expression of genes involved in inflammation, such as cytokines and adhesion molecules. NF-κB comprises a family of dimeric transcription factors that regulate the expression of numerous genes involved in inflammation and cell proliferation .
NF-κB pathways are activated after a potent stimulus from members of the interleukin-1 and tumor necrosis factor superfamilies or lipopolysaccharides, which rapidly degrade IκB within minutes. Degradation of IκB releases NF-κB. After NF-κB is activated, it moves into the nucleus and induces transcription and expression of specific genes, resulting in inflammation, apoptosis, cell proliferation and differentiation, possibly leading to CKD [11,12,13,14,15,16,17,18].
Several studies have shown that inhibitors of NF-κB activation can regulate the inflammatory response of glomerular mesangial cells . The pathogenesis of glomerular mesangial cell inflammation in patients with kidney disease has been associated with NF-κB activation . A recent study showed that when patients with kidney disease have proteinuria, the NF-κB inflammatory reaction and expression of proinflammatory genes are accelerated [21,22,23].
Some genetic studies have shown the association of NF-κB genetic mutations with kidney failure and kidney disease progression [10, 24, 25]. NF-κB is an important transcription factor in inflammation. The polymorphisms in NF-κB transcription binding sites have yet to be identified. Therefore, we used bioinformatics technology and the Taiwan Biobank and chromatin immunoprecipitation sequencing (ChIP-Seq) databases to find NF-κB transcription binding site polymorphisms in a Han population. We then performed a case–control study to investigate the association between the polymorphisms and ESRD.
Bioinformatics analysis in the screening of gene processes
We performed a three-step process for screening of genes (Fig. 1).
Screening of genetic variation in a Taiwanese population through a quality control program
First, we used the single nucleotide polymorphism (SNP) database from the Taiwan Biobank. This database includes 58,917,994 SNPs and 997 next-generation sequencing (NGS) samples. Then, we used a human reference genome downloaded from the University of California, Santa Cruz (UCSC; GRCh37/hg19), and all SNPs within 500 kb upstream and downstream of each candidate SNP from the UCSC genome browser (https://genome.ucsc.edu/). We deleted structural variants (insertion/deletion, deletion) because there was no way to use the multifunctional mass spectrometer (mass array) for genotyping. We kept the remaining variants for study and deleted variants with a call rate of less than 90% at the position. Finally, the remaining SNPs were used for further alignment.
Sequence alignment techniques using bioinformatics analysis of genetic variations that may affect NF-κB binding
Second, we analyzed genetic variants that may affect NF-κB binding by using bioinformatics sequence alignment techniques and identified the variants located in the transcription factor binding site (TFBS). Prior studies have confirmed that the structure of NF-κB is a dimer consisting of five different related structural proteins: p50, p52, p65 (RelA), RelB, and c-Rel. The combination of the p50 protein and the p65 protein is found [26,27,28] in almost all cells; as a result, this study explored only the heterodimer of p50/p65. In the past, the TFBS sequence of the identified transcription factor was 5′-GGGRNNYYCC-3′ (R = A or G; N = A, C, G, or T; Y = C or T). We aligned this motif in all 36,041,790 SNPs and in its nearby sequences within 500 kb that included this motif and found 40,137 SNPs that may affect the binding activity.
Confirmation by ChIP-Seq that these mutations bind to these positions
Third, we further confirmed that these variations do combine with these locations through ChIP-Seq. After the above screening, we used the method of a previous study on the human genome on the NF-κB ChIP-Seq for analysis of the results of further screening . The study was performed using B cells for ChIP-Seq analysis and analysis of the NF-κB five structural proteins of TFBS and was published in the online Gene Expression Omnibus (GEO) database (GSE55105). We extracted the results of the p50–p65 dimer follow-up screening and found a total of 5112 sequences with alignment to 15 SNPs.
Next, we conducted a case–control study to identify SNPs related to the NF-κB binding site associated with ESRD. In this study, we collected blood samples and social demographic data from patients admitted to the Tri-Service General Hospital in Taipei, Taiwan, between 2015 and 2016 and then performed real-time polymerase chain reaction (PCR) and genotyping.
We collected data from 847 hemodialysis (HD) patients (male, 50.8%; female, 49.2%; age, 71.84 ± 12.93) from the Tri-Service General Hospital in Taipei, Taiwan. CKD was defined according to the Kidney Disease Outcomes Quality Initiative definitions, and the estimated glomerular filtration rate (eGFR) was calculated using the Modification of Diet in Renal Disease study equation [30, 31]. Study patients were defined as having an eGFR ≤15 mL/min/1.73 m2 with clinical signs of uremic syndrome requiring HD. All patients were over 20 years old and had been on HD for more than 6 months. Patients were excluded if they had autoimmune disease, malignancy, or acute or chronic infection. Demographic data of the HD patients included age, sex, diabetes, HD duration, hypertension, education level, and blood biochemical values (white and red blood cells, hemoglobin, blood urea nitrogen [BUN], creatinine, albumin, proteinuria, blood sugar AC, triglycerides, cholesterol, sodium, potassium, calcium, phosphorus, and eGFR). The 846 healthy controls (male, 44.6%; female, 55.4%; age, 73.50 ± 7.21) had no history of renal disease, and their eGFR was ≥60 mL/min/1.73 m2. The control group was composed from those undergoing a physical examination at Tri-Service General Hospital. The healthy controls had no microalbuminuria, proteinuria, or hematuria and had normal abdominal/renal ultrasonography findings.
The study was reviewed and approved by the institutional ethics committee of the Tri-Service General Hospital (TSGH-1-104-05-006, TSGH-2-106-05-127). After full explanation of the study, written informed consent was obtained from all participants. All clinical and biological samples were collected, and DNA was genotyped following patient consent.
Genomic DNA extraction and genotyping
The blood samples were extracted from the laboratory by phenol chloroform and stored in a − 20 °C refrigerator for subsequent genotyping and experimental use. Genomic DNA used standard procedures for proteinase K (Invitrogen, Carlsbad, CA, USA) digestion and phenol/chloroform  peripheral blood sample separation; then, the samples were genotyped by iPLEX Gold SNP . We assessed the genotyping experiment quality by intrareplication validation. The concordance rate of interreplication validation of 78 samples (approximately 5%) was 100%. Secondary genotyping was performed on 10 random blood samples with PCR according to a previously described protocol for intrareplication validation . After genotyping replication was conducted twice, the concordance rate was 100% between the two genotyping methods.
Chromatin immunoprecipitation assay and qPCR
We included 9 ESRD patients, GG (n = 3), GT (n = 3), and TT (n = 3), for the ChIP assay. ChIP assays were conducted by using the ChIP Kit ab500 (Abcam, USA) and NF-κB antibody (Proteintech) according to the manufacturer’s instructions. The immunoprecipitate was eluted with 100 μl DNA purifying slurry, and 2 μl of DNA was used in qPCR. Input DNA and NF-κB-enriched DNA fragments were amplified by using qPCR in a 7500 Fast Real-Time PCR System (Applied Biosystems) with primers 5′-ATTCTCACCATGGGAATGG-3′ and 5’GAGGACAGCAAGGTAATAG-3′. The results are shown as percentage input.
Transient transfection and luciferase assay
NF-κB binding site SNP rs9395890 reporter (from 53820675 to 53821295, 620 bp) was amplified by polymerase chain reaction from one home-made genomic DNA library with the primer pair: 5′: 5′-GGGGTACCGCATCTACGTTCTTAAATGGCC-3′ and 3′: 5′-GGAAGATCTCCTACAGAACCATTACACTCTC-3′ and subcloned into a pGL3 basal reporter (Promega, USA) cut at KpnI and BglII sites. After the sequence verification, we further changed the current T allele into G allele using the QuickChange Lightening Site-directed mutagenesis kit (Agilent Technology). HEK293 cells were grown in Dulbecco’s modified Eagle’s medium supplemented with 10% charcoal/dextran-treated fetal bovine serum. The cells in each well (24-well plate) were transfected with total 1 μg DNA and jetPEI (PolyPlus-transfection, Illkirch, France) according to the manufacturer’s protocol. Luciferase activity was assessed after 24 h post transfection using the Promega Luciferase Assay kit and expressed as mean relative light units (RLU) of two transfected sets. Results shown are representative of at least three independent experiments.
We assessed the correlation between genetic variants and mRNA expression of the corresponding genes. Expression quantitative trait loci (eQTL) analysis was also performed using data from the GTEx portal database (https://www.gtexportal.org/home/) and the HapMap Project by a general linear regression model in an additive genetic model .
Statistical analysis was performed with R software, version 3.3.1 (R Project for Statistical Computing, Vienna, Austria). Demographic and clinical data between the groups were compared with Student’s t-test, and the results for continuous variables were given as the mean ± SD. The allele and genotype frequencies between the different groups were compared with the χ2 test when appropriate. The results of ChIP assay qPCR cycles were compared with ANOVA. The genetic polymorphism of ESRD risk was calculated using dominant/recessive models. The odds ratios (ORs) and corresponding 95% confidence intervals (CIs) for assessing the effect of the genotype distribution, allele frequencies and binding site activity on ESRD were calculated by logistic regression analysis with adjustment for relevant significant variables. Statistical significance was defined at the 95% level (P < 0.05).
Screening of genes
We screened for genetic variations in 997 samples from the NGS database in the Taiwan Biobank to determine the total number (58,917,994) of genetic variants in Taiwanese genomes: 11,423,191 were structural variants (insertion/deletion, deletion). There was no way to use the multifunctional mass spectrometer (mass array) for subsequent analysis, and thus we kept only the remaining variants for further study. Therefore, a total of 47,494,803 SNPs were analyzed in detail. Following a quality control program that involved deleting variants with a call rate of less than 90% at the position, 36,041,790 SNPs remained; we then performed sequence alignment analysis.
National Center for biotechnology information
We downloaded the human reference gene sequence of GRCh37/hg19 from the National Center for Biotechnology Information (NCBI) in combination with the human biological database in Taiwan and found that NF-κB (p50–p65) contained 271,063 potential variants in the human genome based on the sequence of the above binding sites. Of the TFBSs, we compared the remaining 36,041,790 SNPs in the first step and found that 3,121,467 SNPs were near the 271,063 potential TFBSs, of which 40,137 SNPs were even in the TFBS of NF-κB. Additionally, mutation of this site will likely result in NF-κB (p50–p65) being unable to bind. Finally, a total of 5766 SNPs with a minor allele frequency > 5% were screened for further follow-up by ChIP-Seq analysis due to the limited number of samples subject to subsequent analysis in this study .
Gene expression omnibus
In the GEO database, there were 5112 positions in the TFBS associated with the p50–p65 dimer. After validating these results with the results of the second stage, the remaining 15 SNPs are shown in Table 1. For SNPs near the DNA sequence, the SNP position as the center ±9 base pairs (the bold font indicates NF-κB) was the expected TFBS. The 15 SNP variations may affect NF-κB binding activity. Finally, we used 15 SNPs obtained from the bioinformatics technology results and the ChIP-Seq database to confirm the relationship with ESRD in this study.
The characteristics of the 846 ESRD and 847 control group subjects are presented in Table 2. The causes of ESRD were diabetes mellitus (DM) in 215 patients (25%), hypertensive nephropathy in 164 (19%), systemic nephropathy in 252 (29%), and other and unknown causes in 136 (16%). There was no significant difference in body mass index. Significant differences in sex, age, DM, hypertension, BUN, serum creatinine, GFR, blood sugar AC, total cholesterol, and triglycerides were observed between patients with ESRD and controls (P < 0.001).
Association analyses of NF-κB binding site gene polymorphisms with susceptibility to ESRD
In the gene screening process, the call rate of all 15 SNPs was > 90%, and the genotypes of these SNPs were in Hardy–Weinberg equilibrium (P > 0.05). When we calculated our sample size, the power was > 50% and the OR was set at 1.5 to detect the real effects of expected NF-κB binding site SNPs. Two SNPs were nonfrequency SNPs under the allele model (rs2851583, rs76552560), and a suitable primer could not be found for three SNPs (rs11234413, rs3826454, rs67087171). Finally, genotyping results were obtained for 10 SNPs. Our results showed that SNP rs9395890 had a significant association with ESRD risk according to genotype (P = 0.041; Table 3).
Allele frequencies for the NF-κB binding site gene polymorphisms with susceptibility to ESRD
There was a significant association (P = 0.049; Table 3) between rs9395890 and ESRD. The SNP rs9395890 with the T allele is associated with ESRD risk in the allele model (P = 0.013; odds ratio [OR] = 1.31, 95% confidence interval [CI] = 1.06 to 1.62; Table 4) and recessive model (P = 0.008; odds ratio [OR] = 1.31, 95% confidence interval [CI] = 1.07 to 1.60; Table 4). There were no significant differences in genotype or allele frequencies in the other nine SNPs between patients with ESRD and controls (Additional file 1).
ChIP assay identified NF-κB binding site rs9395890 enrichment
We included 9 ESRD patients, GG (n = 3), GT (n = 3), TT (n = 3), in the ChIP assay experiments. Real-time qPCR was performed to measure the amount of NF-κB-enriched DNA fragments. The ChIP assay revealed an approximately 1.49-fold enrichment of NF-κB of the variant type TT when compared to that of the wild-type GG in rs9395890 (P = 0.027; TT = 3.20 ± 0.16, GT = 2.81 ± 0.20, GG = 1.71 ± 0.18). When the SNP rs9395890 was type TT, the NF-κB transcription binding activity was higher than that of the GG type (TT: P < 0.001; odds ratio [OR] = 1.50, 95% confidence interval [CI] = 1.20–1.77; GT: P < 0.001; odds ratio [OR] = 1.12, 95% confidence interval [CI] = 1.01–1.38; Table 5, Fig. 2).
Comparison of NF-κB binding activity of T and G allele of the rs9395890 T/G
To establish whether the SNP rs9395890 were functional. We investigated influenced NF-κB binding site activity using a luciferase reporter assay in HEK293 cells. The results are shown in Fig. 3. First, we examined the background of reporter activity in T and G allele by a dose dependent of the amount of pGL3.MLIP-IT1-LUC. It showed that the background activity were higher upon the increasing amount of pGL3.MLIP-IT1-LUC, however, there were no difference in that between T and G allele (Fig. 3a). Further, the functionality of T allele and G allele were observed by upon overexpression of p65 or not in the reporter assay, respectively. Data showed that the luciferase activity in T allele (2.7X) was slightly higher than that in G allele (2.5X), though there are no significant difference (P = 0.589).
The SNP rs9395890 in silico functional validation
When we mapped SNPs using the NCBI database (https://www.ncbi.nlm.nih.gov/projects/SNP/snp_ref.cgi?rs=9395890), rs9395890 was located in intron, chromosome 6, LOC101927189, and the nearest downstream gene was MLIP-IT1. We further conducted MLIP-IT1 mRNA expression quantitative trait loci (eQTL) analysis by searching the GTEx portal (https://www.gtexportal.org/home/). We found that for MLIP-IT1, the T allele was significantly associated with increased expression levels in the following categories: Thyroid (P = 2.6*10− 9; Fig. 4b), Skin - Not Sun Exposed (Suprapubic) (P = 4.8*10− 7; Fig. 4c), Skin - Sun Exposed (Lower leg) (P = 5.5*10− 7; Fig. 4c), Esophagus - Mucosa (P = 1.0*10− 4; Fig. 4a), Minor Salivary Gland (P = 0.1; Fig. 3a), Breast - Mammary Tissue (P = 0.06; Fig. 4a), Adipose - Subcutaneous (P = 0.02; Fig. 4a), Adipose - Visceral (Omentum) (P = 0.2; Fig. 4a), Nerve - Tibial (P = 0.5; Fig. 4a), Vagina (P = 0.7; Fig. 4a), and Small Intestine - Terminal Ileum (P = 0.4; Fig. 4a).
Our results suggested that there is a significant correlation between rs9395890 and ESRD risk. This genetic association study employed bioinformatics technology and epidemiological approaches that make it different from other studies. Previous reports included more genetic and molecular epidemiological studies of ESRD in genome-wide association studies (GWAS). GWAS can explore the etiological contribution of genetic variants throughout the whole genome without applying previously hypothesis. However there are very few detected causal variants . Therefore, we provided an approach to use a hybrid method consisting of candidate gene and epidemiologic approaches.
The research of inflammatory transcription factor (NF-κB) associated SNPs has been investigated in a few previous studies [38,39,40,41]. However, we addressed the importance of genetic polymorphisms in determining ESRD in this study. We were able to identify loci and information about which genes were associated with complex diseases [42, 43]. In our study, we used methodological approach that combined the NGS, NCBI, and GEO online databases to find target SNPs and used epidemiological methods to confirm the findings in a case–control study. A previous study in 2014 also used publicly available genomic data and bioinformatics platforms to provide additional evidence for the TFBSs of SNPs of the ERα-regulating sequence at 21q22.3, which are important in determining breast cancer progression .
Immune and inflammatory factors have important roles in the pathogenesis of kidney diseases [44, 45]. Based on previous studies, the transcription factor NF-κB regulates the expression of various genes that have an important role in the regulation of immunity and inflammation in disease . NF-κB regulates T cells, particularly the T helper 17 cells, which mainly affect the pathogenesis of autoimmunity and inflammation . Several studies have shown the cell-intrinsic role of NF-κB in T cell generation [47, 48]. In the NF-κB pathway, when cells are unstimulated, NF-κB is bound to IkBa and IkBb in the cytoplasm, which prevents NF-κB from entering the nucleus . When these cells are stimulated, specific kinases phosphorylate IkB, allowing degradation by proteasomes [50, 51]. The NF-κB released from IkB results in the passage of NF-κB into the nucleus, and NF-κB binds to target sequences in the promoter regions of target genes, leading to the expression of many genes involved in immune and inflammatory responses .
The NF-κB signaling pathway regulated renal inflammation and the progression of ESRD. Histological evidence of NF-κB activation has been associated with human renal disease with diabetes, glomerular disease, and acute kidney injury . The NF-κB transcription of multiple proinflammatory molecules, such as cytokines, chemokines, allograft antigens, adhesion molecules, and reactive oxygen, in response to renal injury . The SNPs at NF-κB transcription binding site are functional polymorphisms that might regulatory polymorphisms situated in the noncoding regions of the genes which may affect gene product protein due to the transcriptional alterations .
In the past, we knew that the NF-κB transcription factor binding site was involved in the regulation of downstream inflammatory genes, which in turn affected the progression of disease and the deterioration of inflammation. However, the results of this study found an association between ESRD risk and the NF-κB binding site SNP rs9395890. Furthermore, we used a ChIP assay to identify NF-KB binding activity with different genotypes. We found that the NF-KB binding activity at SNP rs9395890 with the TT type was higher than that of the GG type. And we assessed the functionality of the NF-κB binding site rs9395890 T/G polymorphism for effects activity by luciferase reporter assay. Our experimental Fig. 3 showed that the transcriptional activity of the T allele was higher than G allele, but the relative light units (RLU) data was no significant difference between with T and G allele. So far there were no study about the rs9395890 and MLIP-IT1. However it might be the distance between rs9395890 and MLIP-IT1 is too far away.
Furthermore, the results from GTEx portal demonstrated that the T allele was significantly associated with increasing expression levels of rs9395890 in multiple tissues, suggesting that rs9395890 may modulate the risk of ESRD, possibly through a mechanism of modulating gene expression .
SNP rs9395890 is an intron variant located on chromosome 6: 53820994 in front of the MLIP-IT1 gene − 42694 bp. MLIP-IT1 is a noncoding RNA gene, and MLIP-IT1 is a responding gene of rs9395890. Noncoding RNA is not translated into protein but causes transcription factor binding protein and expression of downstream genetics. We suspect that a mutation in this site will affect the function of this gene in MLIP-IT1, which increases the risk of ESRD. To our knowledge, few studies have reported MLIP-IT1 and rs9395890 . DNA is transcribed to mRNA by transcription factors, which then initiate their function. Noncoding RNA occurs during DNA transcription to RNA, when a portion of RNA cannot become mRNA. Noncoding RNA regulates gene transcription function and protein transport. More studies have focused on noncoding RNA and its association with chromatin remodeling, gene transcription, protein transport, and trafficking. Noncoding RNA also has important roles in most human diseases, including coronary artery diseases, autoimmune diseases, neurological disorders, and various cancers [37, 43, 57]. Specifically, we found that the rs9395890 T allele was associated with the risk of ESRD. The T allele mRNA expression levels were higher than those of the G allele in thyroid, skin and mucosa inflammation disease according to data from the GTEx portal. These results are consistent with our ChIP assay data (TT binding activity higher than GG; Fig. 2) .
However, we did not confirm that the NF-κB transcription binding site SNP rs9395890 and the responding gene MLIP-IT1 regulated the mechanism of ESRD risk. Therefore, an experiment to identify the association between rs9395890 MLIP-IT1 RNA expression and ESRD risk is necessary in the future.
Our study has some limitations. First, to our knowledge, no studies have related SNPs of NF-κB transcription binding sites to disease. Our study used bioinformatics technology, that is, the NGS, NCBI, and GEO online databases, to screen transcription binding site genetics. Furthermore, we used our case–control groups for genotyping to confirm that rs9395890 was associated with ESRD. We used GEO database-involved B cells for ChIP-Seq analysis, but it was difficult to obtain renal cells to repeat verification. Second, the odds ratio of rs9395890 was very low, but this is a limitation of an observational study. Third, our study sample size was not large enough. Bonferroni correction could not be performed. However, only one SNP was significantly correlated in our study, and the results of functional analysis were indeed related to ESRD risk. Our study included both a genetic association test and a functional analysis, and the results were consistent (p < 0.05 in both tests). Because of the double statistical test setting, we consider that the type 1 error rate in our setting is less than that in general genetic association studies using Bonferroni correction, and thus the evidence level provided by our study is sufficient even though we cannot conduct Bonferroni correction. In summary, we conclude that SNP rs9395890 plays a key role in the incidence of ESRD.
Our study demonstrated that SNP rs9395890 might contribute to NF-κB transcription binding site ability and might exert an effect on MLIP-IT1 activity. The function of MLIP-IT1 with regard to ESRD progression risk and survival should be explored further.
Availability of data and materials
The data analyzed in this study can be accessed by sending a request to the corresponding author.
Blood urea nitrogen
- ChIP-Seq :
Chromatin immunoprecipitation sequencing
Chronic kidney disease
Estimated glomerular filtration rate
End-stage renal disease
Gene Expression Omnibus
- HD :
National Center for Biotechnology Information
Single nucleotide polymorphisms
Transcription factor binding site
Tri-Service General Hospital
University of California, Santa Cruz
Genome-wide association studies
MD., B., Annual Data Report: Atlas of Chronic Kidney Disease and End-Stage Renal Disease in the United States. National Institutes of Health, National Institute of Diabetes and Digestive and Kidney Diseases: U.S. Renal Data System, USRDS., in U.S. Renal Data System, USRDS. 2013; U.S.
MD., B., Annual Data Report: Atlas of Chronic Kidney Disease and End-Stage Renal Disease in the United States. National Institutes of Health, National Institute of Diabetes and Digestive and Kidney Diseases: U.S. Renal Data System, USRDS. 2015.
Kuo HW, et al. Epidemiological features of CKD in Taiwan. Am J Kidney Dis. 2007;49(1):46–55.
Lin YC, Hsu CY, Kao CC, Chen TW, Chen HH, Hsu CC, Wu MS. Incidence and prevalence of ESRD in Taiwan renal registry data system (TWRDS): 2005-2012. Acta Nephrologica 2014; 28(2): 65–68.
Yang WC, Hwang SJ, Chiang SS, Chen HF, Tsai ST. The impact of diabetes on economic costs in dialysis patients: experiences in Taiwan. Diabetes Res Clin Pract. 2001;54:S47–54.
Satko SG, Sedor JR, Iyengar SK, Freedman BI. Familial clustering of chronic kidney disease. Semin Dial. 2007;20(3):229–36.
Akchurin M, Kaskel F. Update on inflammation in chronic kidney disease. Blood Purif. 2015;39(1–3):84–92.
S., A., Renal disease: environment, race, or genes. Ethn Dis. 2006; 16(S2): 35–39.
Wu HH, Kuo CF, Li IJ, Weng CH, Lee CC, Tu KH, Liu SH, Chen YC, Yang CW, Luo SF, See LC, Yu KH, Huang LH, Zhang W, Doherty M, Tian YC. Family aggregation and heritability of ESRD in Taiwan: a population-based study. Am J Kidney Dis. 2017;70(5):619–26.
Guijarro C, Egido J. Transcription factor-_B (NF-_B) and renal disease. Kidney Int. 2001;59(2):415–24.
Chapter6 Cell Signaling and Apoptosis. Figure 6-F-2. Molecular Biology Web Book. Available from: http://www.web-books.com/MoBio/Free/Ch6F2.htm. Accessed 19 Mar 2018.
Chen F, Castranova V, Shi X. New insights into the role of nuclear factor-κB in cell growth regulation. Am J Pathol. 2001;159(2):387–97.
Ebrahim Zandi MK. Bridging the gap: composition, regulation, and physiological function of the IκB kinase complex. Mol Cell Biol. 1999;19(7):4547–51.
Fu H, Sadis S, Rubin DM, Glickman M, van Nocker S, Finley D, Vierstra RD. Multiubiquitin chain binding and protein degradation are mediated by distinct domains within the 26S proteasome subunit. MCBI J Biol Chem. 1998;273(4):1970–81.
Moscat J, Diaz‐Meco MT, Rennert P. NF-κB activation by protein kinase C isoforms and B-cell function. EMBO Rep. 2003;4(1):31–6.
Kaufman CK, Fuchs E. It's got you covered: NF-κB in the epidermis. J Cell Biol. 2000;149(5):999–1004.
Rothwarf DM, K.M., The NF-_B activation pathway: a paradigm in information transfer from membrane to nucleus. Sci Signal Transduction Knowledge Environ 1999; 5: RE1.
Silverman N, Maniatis T. NF-kappaB signaling pathways in mammalian and insect innate immunity. Genes Dev. 2001;15(18):2321–42.
Massy ZA, Guijarro C, O’Donnell MP, Kim Y, Kashtan CE, Egido J, Kasiske BL, Keane WF. The central role of nuclear factor-kB in mesangial cell activation. Kidney Int Suppl. 1999;71:S76–9.
Li N, Karin M. Ionizing radiation and short wavelength UV activate NF-_B through two distinct mechanisms. Proc Natl Acad Sci U S A. 1998;95(22):13012–7.
Elliott PJ, Pien CS, McCormack TA, Chapman ID, Adams J. Proteasome inhibi-tion: a novel mechanism to combat asthma. J Allergy Clin Immunol. 1999;104(2 Pt 1):294–300.
Hughes K, Antonsson Å, Grundström T. Calmodulin dence of NF_B activation. FEBS Lett. 1998;441:132–6.
Meyer S, Kohler NG, Joly A. Cyclosporine a is an uncompetitive inhibitor of proteasome activity and prevents NF-kappaB activation. FEBS Lett. 1997;413(2):354–8.
Haisong Zhang S-CS. NF-κB in inflammation and renal diseases. Cell Biosci. 2015;5:63.
Misra MK, Mishra A, Pandey SK, Kapoor R, Sharma RK, Agrawal S. Association of functional genetic variants of transcription factor Forkhead box P3 and nuclear factor-κB with end-stage renal disease and renal allograft outcome. Gene. 2016;581(1):57–65.
Wulczyn FG, Naumann M, Scheidereit C. Candidate proto-oncogene bcl-3 encodes a subunit-specific inhibitor of transcription factor NF-kappa B. Nature. 1992;358(6387):597–9.
Hayden MS, Ghosh S. Signaling to NF-kappaB. Genes Dev. 2004;18(18):195–224.
Perkins ND. Post-translational modifications regulating the activity and function of the nuclear factor kappa B pathway. Oncogene. 2006;25(51):6717–30.
Zhao B, Barrera LA, Ersing I, Willox B, Schmidt SC, Greenfeld H, Zhou H, Mollo SB, Shi TT, Takasaki K, Jiang S, Cahir-McFarland E, Kellis M, Bulyk ML, Kieff E, Gewurz BE. The NF-κB genomic landscape in Lymphoblastoid B-cells. Cell Rep. 2014;8(5):1595–606.
Basseres DS, Baldwin AS. Nuclear factor-kappaB and inhibitor of kappaB kinase pathways in oncogenic initiation and progression. Oncogene. 2006;25(51):6817–30.
Levey AS, Bosch JP, Lewis JB, Greene T, Rogers N, Roth D. A more accurate method to estimate glomerular filtration rate from serum creatinine: a new prediction equation. Modification of diet in renal disease study group. Ann Intern Med. 1999;130(6):461–70.
Tan SC, Yiap BC. DNA, RNA, and protein extraction: the past and the present. J Biomed Biotechnol. 2009:574398. https://doi.org/10.1155/2009/574398.
Perkel J. SNP genotyping: six technologies that keyed a revolution. Nat Methods. 2008;5(5):447–54.
Yang HY, Lu KC, Fang WH, Lee HS, Wu CC, Huang YH, Lin YF, Kao SY, Lai CH, Chu CM, Su SL. Impact of interaction of cigarette smoking with angiotensin-converting enzyme polymorphisms on end-stage renal disease risk in a Han Chinese population. J Renin-Angiotensin-Aldosterone Syst. 2015;16(1):203–10.
GTEx-Portal online database. 2018. Available from: https://www.gtexportal.org/home/. Accessed 19 Feb 2019.
Mathelier A, Zhao X., Zhang AW, Parcy F, Worsley-Hunt R, Arenillas DJ, Buchman S, Chen CY, Chou A, Ienasescu H, Lim J, Shyr C, Tan G, Zhou M, Lenhard B, Sandelin A, Wasserman WW., JASPAR 2014 an extensively expanded and updated open-access database of transcription factor binding profiles. Nucleic Acids Res. 2014; 42: D142–D147.
Fugger L, McVean G, Bell JI. Genomewide Association Studies and Common Disease — Realizing Clinical Utility. N Engl J Med. 2012;367:2370–1.
Udalova IA, Mott R, Field D, Kwiatkowski D. Quantitative prediction of NF-κB DNA– protein interactions. PNAS. 2002;99(12):8167–72.
Wong D, Teixeira A, Oikonomopoulos S, Humburg P, Lone IN, Saliba D, Siggers T, Bulyk M, Angelov D, Dimitrov S, Udalova IA, Ragoussiscorresponding J. Extensive characterization of NF-κB binding uncovers non-canonical motifs and advances the interpretation of genetic functional traits. Genome Biol. 2011;12(7):R70.
Jamshidi M, Fagerholm R, Khan S, Aittomäki K, Czene K, Darabi H, Li J, et al. SNP-SNP interaction analysis of NF-κB signaling pathway on breast cancer survival. Oncotarget. 2015;6(35):37979–94.
Bečanović K, Nørremølle A, Neal SJ, Kay C, Collins JA, Arenillas D, Lilja T, Gaudenzi G, Manoharan S, et al. A SNP in the HTT promoter alters NF-κB binding and is a bidirectional genetic modifier of Huntington disease. Nat Neurosci. 2015;18(6):807–16.
Ortiz-Barahona A, Villar D, Pescador N, Amigo J, del Peso L. Genome-wide identification of hypoxia-inducible factor binding sites and target genes by a probabilistic model integrating transcription-profiling data and in silico binding site prediction. Nucleic Acids Res. 2010;38(7):2332–45.
Hsiung CN, Chu HW, Huang YL, Chou WC, Hu LY, Hsu HM, Wu PE, Hou MF, Yu JC, Shen CY. Functional variants at the 21q22.3 locus involved in breast cancer progression identified by screening of genome-wide estrogen response elements. Breast Cancer Res. 2014;16(5):455.
Kurts C, Panzer U, Anders HJ, Rees AJ. The immune system and kidney disease: basic concepts and clinical implications. Nature Reviews Immunol. 2013;13:738–53.
Imig JD, Ryan MJ. Immune and inflammatory role in renal disease. Compr Physiol. 2013;3(2):957–76.
Dong C. TH17 cells in development: an updated view of their molecular identity and genetic programming. Nat Rev Immunol. 2008;8(5):337–48.
Isomura I, Palmer S, Grumont RJ, Bunting K, Hoyne G, Wilkinson N, Banerjee A, Proietto A, Gugasyan R, Wu L, McNally A, Steptoe RJ, Thomas R, Shannon MF, Gerondakis S. c-Rel is required for the development of thymic Foxp3+ CD4 regulatory T cells. J Exp Med. 2009;206(13):3001–14.
Long M, P.S., Strickland I, Hayden MS, Ghosh S., nuclear factor-kappaB modulates regulatory T cell development by directly regulating expression of Foxp3 transcription factor. Immunity. 2009; 31(6): 921–931.
Baldwin AS Jr. The NF-kappa B and I kappa B proteins: new discoveries and insights. Annu Rev Immunol. 1996;14:649–83.
Chen ZJ, Parent L, Maniatis T. Site-specific phosphorylation of IkappaBalpha by a novel ubiquitination-dependent protein kinase activity. Cell. 1996;84:853–62.
DiDonato J, Mercurio F, Rosette C, Wu-Li J, Suyang H, Ghosh S, Karin M. Mapping of the inducible IkappaB phosphorylation sites that signal its ubiquitination and degradation. Mol Cell Biol. 1996;16(4):1295–304.
Barnes PJ, Karin M. Nuclear factor-kB - a pivotal transcription factor in chronic inflammatory disease. N Engl J Med. 1997;336(15):1066–71.
Sanz AB, Sanchez-Niño MD, Ramos AM, Moreno JA, Santamaria B, Ruiz-Ortega M, Egido J, Ortiz A. NF-kappaB in renal inflammation. J Am Soc Nephrol. 2010;21(8):1254–62.
Li Q, Verma IM. NF-kappaB regulation in the immune system. Nat Rev Immunol. 2002;2(10):725–34.
Chorley BN, Wang X, Campbell MR, Pittman GS, Noureddine MA, Bell DA. Discovery and verification of functional single nucleotide polymorphisms in regulatory genomic regions: current and developing technologies. Mutat Res. 2008;659(1–2):147–57.
GRCh37. Human. 2018. https://www.ncbi.nlm.nih.gov/assembly/GCF_000001405.13/.
Li J, Xuan Z, Liu C. Long non-coding RNAs and complex human diseases. Int J Mol Sci. 2013;14(9):18790–808.
A particular acknowledgment to the medical staff of the dialysis centers involved in the study.
Ethical approval and consent to participate
The study was reviewed and approved by the institutional ethics committee of the Tri-Service General Hospital (TSGH-1-104-05-006, TSGH-2-106-05-127). All the patients were informed about the nature and the aim of the study and gave their consent to participate in written format.
This study was supported by grants from the Taoyuan Armed Forces General Hospital, Cheng Hsin General Hospital, National Defense Medical Center, Taiwan, ROC, and Ministry of Science and Technology, Taiwan, ROC (TAFGH-10429, TAFGH-10801, CH-NDMC-105-5, CH-NDMC-106-4, CH-NDMC-107-13, CH-NDMC-108-6, MAB-1050113, MAB-106-105, MAB-107-074, MAB-108-047, MOST107–2314-B016–052-MY3). These sponsorships are used for stationery and printing, research staff salary, English editing.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Yang, JH., Chen, WT., Lee, MC. et al. Investigation of the variants at the binding site of inflammatory transcription factor NF-κB in patients with end-stage renal disease. BMC Nephrol 20, 300 (2019). https://doi.org/10.1186/s12882-019-1471-2