Skip to main content
Download PDF
  • Research article
  • Open access
  • Published:

Changes in DNA methylation in naïve T helper cells regulate the pathophysiological state in minimal-change nephrotic syndrome

BMC Research Notes volume 10, Article number: 480 (2017) Cite this article

Abstract

Background

DNA methylation plays a crucial role in regulating transcription, and changes in DNA methylation affect gene expression and disease development. Minimal change nephrotic syndrome (MCNS) has been reported to involve immunological disturbances. Since the characteristic features of the disease include recurrent relapse and sex and age preference, the disease pathogenesis may be partly related to epigenetic changes. However, little is known about these changes.

Methods

We analyzed genome-wide DNA methylation using the microarray-based integrated analysis of methylation by isoschizomers method. This method was used to evaluate methylation in monocytes (patient number; n = 6) and naïve T helper cells (n = 4) from the peripheral blood of MCNS patients both in relapse and following remission and that of healthy controls (n = 5).

Results

In total, 85 co-occurring genes were identified in naïve T helper cells, while 4 such genes were identified in monocytes, which were common among the 3 following comparisons for changes in DNA methylation using sample pairs: (1) relapse versus remission, (2) relapse versus controls, and (3) remission versus controls. In 82 of 85 co-occurring genes (96.5%) in naïve T helper cells, the level of DNA methylation was altered according to disease activity, but was not related to disease activity in the 4 genes detected in monocytes.

Conclusions

Therefore, in 82 co-occurring genes in naïve T helper cells, the regulation of DNA methylation was well correlated with the clinical and pathophysiological state. Our genome-wide approach to analyze DNA methylation provides further insight into the pathogenesis of MCNS and indicates potential prediction and diagnostic tool for the disease.

Background

Minimal-change nephrotic syndrome (MCNS) is the most common cause of nephrotic syndrome in children. It is characterized by massive proteinuria and hypoalbuminemia in a relapse/remission course without histological evidence of immune-mediated inflammatory damage. These manifestations are typically reversible with corticosteroid therapy. Although the pathogenesis of MCNS remains unknown, immunological disruption has been implicated in this disease [1], and T cell-derived vascular permeability factors have been shown to be responsible for alterations in glomerular permeability [2,3,4]. During childhood, the incidence of MCNS has been reported to be twofold higher in male children, with a prevalence that is inversely proportional to age. Recurrent relapse tends to decrease after adolescence [5, 6]. The characteristic features of MCNS include a recurrent relapse/remission course, gender preference, age preference of onset and relapse, and steroid response in most patients, therefore, a single genetic defect is likely not responsible for the disease. However, epigenetic alterations may occur without a direct change in the genetic sequence and regulate the phenotype in relapsed MCNS patients.

Epigenetics is the study of mitotically heritable changes, which alter gene expression without direct DNA sequence alterations. DNA methylation, one of the principal epigenetic mechanisms in mammals, involves the covalent addition of a methyl group to a cytosine residue followed by a guanine residue [7]. DNA methylation regulates gene expression and is essential for differentiation, embryonic development [8], genomic imprinting [9], and X-chromosome inactivation [10]. DNA methylation within the promoter region of a gene is commonly associated with transcriptional inactivation, whereas demethylation contributes to transcriptional activation. Changes in the DNA methylation profile can also lead to differences in gene expression programs and thus influence disease development [11].

Audard et al. reported that nuclear factor related to kappa B binding protein (NFRKB) is highly expressed in the nuclear compartment during relapse and that NFRKB promotes the hypomethylation of genomic DNA, suggesting epigenetic involvement in the pathogenesis of MCNS [12]. Elie et al. suggested that alterations in epigenetic modifications induced by external cues such as viral infection affect MCNS progression; this is supported by the fact that MCNS relapse is frequently triggered by external or internal environmental factors. Such environmental changes affect the epigenotype and alter gene expression [13]. We previously reported that the DNA methylation states of 3 genes, GATA binding protein 2 (GATA2), pre-B cell leukemia homeobox 4 (PBX4), and nyctalopin (NYX) in naïve T helper cells (Th0s), but not in monocytes, significantly differed between relapse and remission in affected patients [14]. Additionally, epigenetic regulation in Th0s underlies the pathogenesis of MCNS, the disturbance of which has been implicated in MCNS development [14].

In this study, we examined whether DNA methylation states change according to the clinical states of MCNS. We identified candidate genes whose DNA methylation states change during relapse and remission and compared them to those in healthy controls to gain insight into MCNS pathogenesis.

Results

Study population characteristics

The mean total follow-up period for patients until December 2012 was 122.2 months (range 53–204 months). All patients were steroid-responsive and experienced frequent relapses. The mean total relapse time was 8.3 (range 4–15 times) during the follow-up period. Although it has been already known that the likelihood that a child with primary nephrotic syndrome who responds during eight weeks of intensive initial steroid treatment has MCNS is quite high and that renal biopsy need not be performed as prognosis in these patients is considered to be very favorable [15], renal biopsy was performed for histological evaluation of the side effects of cyclosporine in 4 patients who were administered with cyclosporine for frequent relapse during the follow-up periods. The histological findings were consistent with those of MCNS. The mean age at sampling in relapse was 13 years 5 months (range 6–19 years and 10 months) and the mean number of relapses until sampling was 6 (range 3–13 relapses) (Table 1). The mean sampling interval from relapse to subsequent remission was 25 weeks (range 5–71 weeks). The therapeutic conditions were similar at sampling for relapse and remission in all subjects except for patients 2 and 6, who received a corticosteroid (patients 2 and 6) and an immunosuppressant (patient 2) at remission, but not at relapse sampling.

Table 1 MCNS patient characteristics at sampling
Full size table

Microarray-based integrated analysis of methylation by isoschizomers (MIAMI)

We first investigated differences in genome-wide DNA methylation between relapse and remission of MCNS patients (Fig. 1a–c), between relapse and controls (Fig. 1d–f), and between remission and controls (Fig. 1g–i), both in monocytes (Fig. 1b, e, h) and in Th0s (Fig. 1c, f, i). For each comparison, the numbers of probes with significant changes in DNA methylation have been shown in Table 2. The MIAMI results were confirmed using bisulfite-pyrosequencing analysis as previously reported [14].

Fig. 1
figure 1

Three sets of comparisons with microarray-based integrated analysis of methylation by isoschizomers (MIAMI). The 3 sets of comparisons are indicated schematically: comparison between a relapse and remission, d relapse and controls, and g remission and controls. Scatter plots of the signals obtained with each probe in MIAMI analysis are represented; comparisons between relapse and remission in monocytes (b) and in naïve T helper cells (Th0s) (c), between relapse and controls in monocytes (e) and in Th0s (f), and between remission and controls in monocytes (h) and Th0s (i) are shown. The y-axis represents the signal intensity ratios of each probe between 2 samples digested with HpaII, a methyl-sensitive restriction enzyme that indicates the differences in DNA methylation between the 2 samples. The x-axis represents the ratios between the same set of samples digested with MspI, which is a methyl-insensitive isoschizomer of HpaII. Thus, the x-axis indicates the digestion efficiency or accessibility of restriction enzymes. For statistical analysis, we used the Mahalanobis distance, which is a simplistic approach for estimating the standard deviation of distance between each sample point and the center of mass. The threshold values were determined at a 99% CI of the Mahalanobis distance from the center of the mass (solid line) and 99% CI of the MspI treatment signal plotted on the x-axis (dash line). The probes in the upper area of the 99% CI of Mahalanobis distance and concomitantly within 99% CI of MspI digestion corresponded to significantly less methylated in relapse (b, c, e, and f) and in remission (h and i) than each comparative. In contrast, the probes in the lower area indicated significantly more methylated in the relapse (b, c, e, and f) or in remission (h and i). The numbers in each figure indicate detected probe numbers in each area; (minus) indicates less methylated probes, while (plus) indicates more methylated probes. The number of probes detected is shown in Table 2

Full size image
Table 2 Number of probes in which DNA methylation changed significantly
Full size table

For all the probes, the Mahalanobis distance distributions were significantly different (P < 0.0001) between monocytes and Th0s in all comparisons (Fig. 2a–c). These results indicate that the regulation of DNA methylation differs between monocytes and Th0s in all comparisons. Specifically, regulation differed between relapse and remission, even in the same individuals (Fig. 2a), as we previously reported [14], as well as between patients both in relapse and in remission and healthy controls (Fig. 2b, c).

Fig. 2
figure 2

Distance distributions between monocytes and naïve T helper cells (Th0s) in the 3 comparisons. The distance distributions of all probes from the center of the mass were significantly different (P < 0.0001) between monocytes and Th0s in each of the 3 sets of comparisons: comparison between a relapse and remission, b relapse and controls, and c remission and controls

Full size image

Genes co-occurring among three comparisons in monocytes and Th0s derived from patients in relapse and remission, and healthy controls

We investigated whether there are probes in monocytes (Fig. 3a) and Th0s (Fig. 3b) that could be detected to be common among the 3 comparisons. The number of genes in Th0s that were identified in all 3 comparisons was 85 (Fig. 3b); however, only 4 genes in monocytes were identified to be common among all 3 comparisons (Fig. 3a). In this paper, we refer to these genes as ‘co-occurring’ genes. The distribution of the 85 probes was not skewed, indicating that the results were not affected by technical bias (Fig. 3c). Further analysis of the DNA methylation status in these genes revealed that 82 of 85 co-occurring genes in Th0s (Table 3) showed changes in methylation according to disease activity.

Fig. 3
figure 3

Co-occurring genes* among the 3 comparisons in monocytes and Th0s. The number of probes plotted outside the 99% CI of the Mahalanobis distance within MspI 99% CI among the 3 comparisons are shown in a Venn diagram for monocytes (a) and Th0s (b). There were 85 co-occurring genes in Th0s in the 3 comparisons and 4 in monocytes. c The probes of the 85 co-occurring genes are indicated in the scatter plots comparing relapse and remission as a representative of 3 comparisons. The distributions of the 85 probes were equitable and appropriate. *“Co-occuring” genes are those identified to be common among all 3 comparisons

Full size image
Table 3 List of the 82 genes co-occurring among the three sets of comparisons in Th0 cells
Full size table

Among the 82 genes, the degree of methylation in 40 genes was higher in relapse than in remission (Fig. 4a) and controls (Fig. 4b), and higher in remission than in controls (Fig. 4a), i.e. relapse > remission > controls. The difference between the control and relapse (Fig. 4b) was comparable to the sum of the differences between the control to remission and remission to relapse (Fig. 4a, see c for the concept). The remaining 42 genes were less methylated in relapse than in remission (Fig. 4d) and controls (Fig. 4e) and less methylated in remission than in controls (Fig. 4d), i.e. relapse < remission < controls. Signal intensity ratios between relapse and controls (Fig. 4e) showed similar cumulative differences between relapse to remission and remission to controls (Fig. 4d, see f for the concept), indicating the accuracy of the 3 different experiments comparing relapse versus remission, relapse versus control, and remission versus control (Additional file 1: Table S1). The remaining 3 genes from the 85 genes detected in Th0s (Additional file 1: Table S1) and the 4 co-occurring genes in monocytes showed no specific changes in DNA methylation in association with disease activity (data not shown).

Fig. 4
figure 4

DNA methylation changes in 82 of 85 genes detected in Th0s were in accordance with disease activity. Among the 85 genes co-occurring in Th0s, 82 showed changes in methylation status according to disease activity. Forty of the 82 genes showed significantly higher methylation at relapse compared to remission (a) and controls (b), or at remission compared to controls (a) in each comparison. Additionally, the remaining 42 of the 82 genes were always less methylated at relapse compared to remission (d) and controls (e), or at remission compared to controls (d) in each comparison. The difference between control and relapse (b and e) was comparable to the sums of differences between control to remission and remission to relapse (a and d) at every probe of 82 genes, both in increasing and decreasing methylation towards relapse, as shown in c and f. The line drawn in the same color indicates the methylation ratio obtained from the same probe spotted in each array for 3 comparisons (a, b, d, and e)

Full size image

Ingenuity pathway analysis (IPA) of 85 co-occurring genes in Th0s

The 85 co-occurring genes in Th0s were analyzed by IPA [16, 17], which generated 2-network diagrams. The associated functions of the genes in network 1, which contained 19 of the 85 genes, included cell-mediated immune response, cellular development, and cellular function and maintenance (Fig. 5a), with a score of 47. Genes in network 2, which contained 15 of the 85 genes, had roles in gene expression, cellular development, and cellular growth and proliferation (Fig. 5b), with a score of 31. Top molecules and top canonical pathways detected by IPA are shown in Table 4.

Fig. 5
figure 5

Two networks diagrams generated by ingenuity pathway analysis with the 85 genes. The associated network functions of network 1 (a) were cell-mediated immune response, cellular development, and cellular function and maintenance and included 19 of the 85 genes, with a score of 47, and those of network 2 (b) were gene expression, cellular development, and cellular growth and proliferation, and included 15 of the 85 genes, with a score of 31. The networks are based on known protein–protein interactions and functional relationships. The intensity of gene (node) color in the networks indicates the degree of hypermethylation (red) or hypomethylation (gray). Nodes are displayed using various shapes that represent the functional class of gene products

Full size image
Table 4 Ingenuity pathway analysis of the 85 co-occurring genes
Full size table

Co-occurring genes between Th0s isolated from patients in relapse and remission, and healthy controls

We analyzed 206 genes that were common in the 2 comparisons, including in relapse versus remission and relapse versus controls (Additional file 2: Table S2), and 58 genes in relapse remission versus controls (Additional file 3: Table S3) in Th0 cells, shown in the Venn diagram in Fig. 3b. The former 206 genes are proteinuria-related genes, while the latter 58 genes are disease-related genes. Among the 206 co-occurring genes, 186 genes were less methylated in the relapse (proteinuric) state than in the remission and control groups (non-proteinuric state), while the remaining 20 were more methylated in the relapse (proteinuric) state than in other samples (non-proteinuric) (Additional file 2: Table S2). Furthermore, among the latter 58 co-occurring genes, 20 were less methylated in both relapse and remission (disease) than in controls (non-disease), while the remaining 37 genes were more methylated in relapse and remission (disease) than in controls (non-disease) (Additional file 3: Table S3). The consistency of the results from the 2 different experiments using each probe indicates the accuracy of the experiment and the high probability of detection. We further examined these genes using DAVID bioinformatics resources 6.7 [18]. The 5 KEGG pathways identified the 206 proteinuria-associated genes, cytokine-cytokine receptor interaction (hsa04060), and hematopoietic cell lineage (hsa04640), glycine, serine, and threonine metabolism (hsa00260), and 56 disease-associated genes, including 1 KEGG pathway of focal adhesion (hsa04510) (Additional file 4: Table S4, Table S5).

Discussion

By analyzing DNA methylation in Th0s and monocytes sorted from peripheral blood mononuclear cells (PBMCs) isolated from MCNS patients and healthy controls, we successfully identified 82 genes in Th0s whose DNA methylation levels varied according to disease activity, i.e., relapse, remission, or healthy (Fig. 6). However, in monocytes, only 4 genes were co-occurring. It has long been implicated that helper T cells rather than monocytes are involved in MCNS pathogenesis [1,2,3,4]. The current results suggest that the epigenetic status of helper T cells is altered and is associated with MCNS pathogenesis. DNA methylation analysis possibly enabled the identification of relevant genes associated with disease activity, particularly between remission and healthy controls, despite situations in which patients showed no symptoms and no proteinuria, such that their physical status was similar to that of healthy controls. Genes were less likely to be detected in the expression analysis when they caused similar phenotypes and clinical symptoms.

Fig. 6
figure 6

Scheme of correlation between DNA methylation change in 82 genes and disease activity. DNA methylation of 82 genes identified in Th0s varied according to disease activity, i.e., relapse, remission, or healthy, in 42 genes increasing and in 40 genes decreasing DNA methylation in relapse

Full size image

We previously identified 3 genes (GATA2, PBX4, and NYX) showing significantly lower methylation in Th0s in the relapse state than in the remission state by using regression analysis [14], which was originally used in the MIAMI method [19, 20]. In the current study, we used the Mahalanobis distance rather than regression analysis and determined that probes showing 99% CI were significant. Mahalanobis distance is the distance of the sample points from the center of the mass. Therefore, this value shows greater correlation with the differences in DNA methylation ratios between 2 samples in paired comparisons compared to regression analysis. The distributions of distance for all probes from the center of the mass were significantly different (P < 0.0001) between monocytes and Th0s, even though they were from the same individuals, and were significantly greater in Th0s than in monocytes in all 3 comparisons (Fig. 2a–c). This indicates that the steady-state regulation of DNA methylation differs between monocytes and Th0s and that more significant changes in DNA methylation were observed in Th0s than in monocytes in each comparison.

MIAMI is used for semi-quantitative and comparative analysis of DNA methylation. In 82 (95.6%) of the 85 co-occurring genes in 3 permutations of 3 Th0s samples, DNA methylation changes occurred in the same direction and was consistent in all 3 comparisons, and these changes were in accordance with the level of disease activity. This finding indicates that DNA methylation of the 82 genes is very relevant to MCNS development. Furthermore, 2 associated network functions were identified by IPA, involving 34 of the 85 co-occurring genes. The network functions contained nuclear factor kappa B (NF-κB) (Fig. 5a) and tumor necrosis factor alpha (TNF-α) (Fig. 5b) networks, both of which have been previously implicated in MCNS pathogenesis [21,22,23,24,25]. Particularly, TNF was less methylated in relapse than in remission, which is consistent with previous reports of elevated TNF-α levels in relapse [25, 26]. Moreover, the proteinuric state was associated with 206 co-occurring genes, while the disease was associated with 58 co-occurring genes, demonstrating that the direction of changes in DNA methylation in Th0s was consistent with disease activity (Fig. 3b) (Additional file 2: Table S2, Additional file 3: Table S3).

In this study, we focused on immunocytes and their epigenetic states. We purified and used 2 different cell types from PBMCs isolated from the participants: monocytes, which are precursors of antigen-presenting cells derived from the myeloid cell series, and Th0s, derived from the lymphoid system, because different types of blood cells have different methylation patterns that reflect differences in gene expression during cell differentiation into specific lineages [14]. An increasing body of evidence suggests that epigenetic regulation is essential for the differentiation of blood cells such as T helper type 2 (Th2) cells [27, 28]. Therefore, we isolated and studied monocytes and Th0s rather than total PBMCs, which contain a mixture of various cell types in different developmental stages and in differing epigenetic states. This restricted the analysis to epigenetic changes related to disease development, avoiding the inclusion of changes caused by differences in cell composition in the analysis with isolated cell types. Although Th2 activation has been shown to be involved in MCNS [21, 29,30,31,32] we hypothesized that epigenetic changes occur in Th0s before differentiation into Th2 cells, predisposing individuals to the active state of the disease. Examination of Th0s revealed that relevant genes showed modified DNA methylation states associated with disease activity. This is consistent with the results of a previous study showing that immature lymphocytes play a role in MCNS pathogenesis [33].

We found that variations in DNA methylation are associated with disease activity, proteinuria, and disease development in Th0s rather than in monocytes. Furthermore, the functions of the 82 genes associated with disease activity were found to be related to the TNF-α or NF-κB pathways, which have been reported to be related to MCNS pathogenesis. These results may suggest that the DNA methylation changes in Th0s were not only concomitant with disease activity, but also functionally relevant.

The limitation of this study is that limited number of patients and heterogeneity in treatment states were included. Genes that are sensitive to prednisolone and cyclosporine exposure in DNA methylation may be included in the detected genes. Further studies that include a larger number of MCNS patients without any immunotherapies both at relapse and remission are required for the identification of candidate genes for MCNS.

Conclusions

Based on our data though it is preliminary, regulation of DNA methylation in Th0s may underlie the pathogenesis and pathophysiological states of MCNS. This suggests that DNA methylation analysis could be potentially a useful diagnostic tool that can be used to distinguish different disease states, particularly between remission and cured states.

Methods

Patients

Samples for MIAMI analysis were obtained from 6 male patients with MCNS (Table 1) at relapse and following complete remission. All patients had developed nephrotic syndrome before 16 years of age, and were diagnosed according to the criteria of the International Study of Kidney Disease in Children. In brief, the diagnostic criteria for MCNS are (1) heavy proteinuria of ≥0.96 g/day m2 BSA, determined using a urine sample collected overnight; (2) hypoalbuminemia of ≤2.5 g/dL; (3) no evidence of underlying systemic diseases nor exposure to agents that are known to be associated with NS at its onset; (4) response during 8 week of the initial treatment with 60 mg/24 h m2 prednisone in three divided daily doses for 4 week, followed by administration of 40 mg/24 h m2 in divided doses on 3 consecutive days per week for 4 week. Response was defined as the reduction in the rate of urinary excretion of protein to <4 mg/h m2 (Albustix, 0 to trace) for 3 consecutive days. Relapse involved the reappearance of proteinuria of 40 mg/h m2 (Albustix, + +or greater) for 3 consecutive days [15]. A frequent relapser is a patient who experiences 2 or more relapses during a 6-month period subsequent to treatment [15].

Informed consent was obtained from all parents for the children and from the patients who were older/adolescents. The healthy controls included five male voluntary participants who were 23 years old. This study was approved by the Ethics Committee of Gunma University Graduate School of Medicine, Japan (Receipt Number 89).

Cell separation

PBMCs were isolated from 20 mL samples of anti-coagulated blood obtained by gradient separation using the Lymphprep™ Tube (Axis-Shield PoC AS, Oslo, Norway). Monocytes and Th0s were separated using an autoMACS Pro Separator (Miltenyi Biotec, Bergisch Gladbach, Germany) from PBMCs that had been magnetically labeled with CD14, CD4, and CD45RO MicroBeads (Miltenyi Biotec). To obtain monocytes, the CD14-positive (CD14+) fraction was collected as monocytes, whereas CD14-negative, CD4-positive, and CD45RO-negative (CD14−CD4+CD45RO−) fractions were collected as Th0s. Flow cytometric analysis using a MACSQuant® Analyzer (Miltenyi Biotec) showed that the precision of cell separation was 96.2% for CD14+ cells and 94.46% for CD4-positive and CD45RA-positive cells as a CD45RO-negative fraction [14].

CD14+ monocytes were obtained from all 6 patients, while CD14−CD4+CD45RO− Th0 cells were obtained from 4 of these patients (Table 1a) for MIAMI analysis, both at relapse and following complete remission. Genomic DNA (gDNA) was extracted from these cells as previously described [34]. Consistent amounts of extracted gDNA (300 ng from monocytes and 250 ng from Th0s) from each selected cell type were pooled and used for subsequent DNA methylation analysis using MIAMI. This was done to exclude the influence of epigenetic differences among cell types, identify specific changes resulting from the clinical courses of MCNS, and reduce noise caused by individual differences.

MIAMI analysis

To determine methylation states, we used the MIAMI method, which provides high-throughput global analysis of DNA methylation [19, 20]. Since this method is based on the comparison of the digestion efficiency between 2 different samples by methyl-sensitive HpaII restriction enzymes, it indicates the differences in DNA methylation between the samples. The method was performed as described previously using 1.8 and 1.0 µg of pooled gDNA from the monocytes of 6 patients and Th0s of 4 patients, respectively [14, 19, 20]. Briefly, isoschizomers (HpaII and MspI) that recognize the same recognition site (CCGG) were used. Pooled gDNA was digested with HpaII, a methylation-sensitive restriction enzyme that cleaves only unmethylated DNA, and then the samples were adapter-ligated and amplified by PCR using primers designed against the adapter sequences. The samples were further digested with MspI, a methylation-insensitive enzyme that digests CCGG sites irrespective of their methylation status, and were amplified again with the same primer set (HpaII–MspI treatment). The second MspI treatment yielded amplicons from unmethylated DNA fragments only. Hence, only HpaII-cleavable unmethylated DNA fragments were amplified and were quantified based on their fluorescence intensities by microarray analysis. The amplified products were then labeled with Cy3 (remission or control samples) or Cy5 (relapse or remission samples) and co-hybridized (i.e. remission vs relapse, control vs relapse, or control vs remission) to a microarray spotted with 38,172 sixty-mer oligonucleotides covering the vicinity of the transcription start sites of 14,978 genes. Following hybridization, the microarray was scanned and the fluorescence intensities obtained were quantified and normalized. The same pooled gDNA samples were treated first with MspI rather than HpaII (MspI–MspI treatment) and analyzed on a duplicate array to correct for false-positives caused by single-nucleotide polymorphisms or incomplete digestion.

IPA

Gene networks and canonical pathways representing key genes were identified using the curated IPA database (http://www.ingenuity.com/). The data set containing gene identifiers and corresponding fold-changes was uploaded into the web-delivered application, and each gene identifier was mapped to its corresponding gene object in the Ingenuity Pathways Knowledge Base. IPA generates hypothetical protein–protein interactions, which are not limited to “binding”, but also include “activation”, “inhibition”, “expression”, and other interactions described in the literature. Functional analysis identified the biological functions and/or diseases that were most significant to the data sets. Fisher’s exact test was performed to calculate a P value to indicate the probability that each biological function and/or disease assigned to the data set was because of chance alone. The data set was mined for significant pathways using the IPA library of canonical pathways, using either (1) a ratio of the number of genes from the data set that mapped to the pathway divided by the total number of genes that mapped to the canonical pathway or (2) a Fisher’s exact test to calculate a P value determining the probability that the association between the genes in the data set and the canonical pathway was explained by chance alone [19, 20].

DAVID bioinformatics resources

DAVID bioinformatics resources (http://david.abcc.ncifcrf.gov) consist of an integrated biological knowledgebase and analytic tools aimed at systematically extracting biological meaning from large gene/protein lists. DAVID is a high-throughput and integrated data-mining environment that can analyze gene lists derived from high-throughput genomic experiments. The procedure first requires uploading a gene list containing any number of common gene identifiers, followed by analysis using one or more text and pathway-mining tools such as gene functional classification, functional annotation chart or clustering, and a functional annotation table. The major statistical methods have been described previously [18], and Additional file 2: Table S2 shows the associated parameters used in DAVID.

Statistical analysis

Mahalanobis distance was calculated using R version 2.0. The threshold values were determined based on 99% confidence intervals (Fig. 1). Statistical analyses of the distance distributions from the center of mass (Mahalanobis distance) for each probe used in the assays were performed using a non-parametric Mann–Whitney U test (Fig. 2). Statistical analyses were performed using GraphPad PRISM5 software (GraphPad Software, La Jolla, CA, USA) with the significance level set at P < 0.05.

References

  1. Shalhoub RJ. Pathogenesis of lipoid nephrosis: a disorder of T-cell function. Lancet. 1974;2:556–9.

    Article  CAS  PubMed  Google Scholar 

  2. Heslan JM, Branellec A, Laurent J, et al. The vascular permeability factor is a T lymphocyte product. Nephron. 1986;42:187–8.

    Article  CAS  PubMed  Google Scholar 

  3. Maruyama K, Tomizawa S, Shimabukuro N, et al. Studies of vascular permeability factor and inhibitory effect of supernatants derived from T lymphocytes culture in minimal change nephrotic syndrome on rat kidney capillaries. Nephron. 1989;51:73–6.

    Article  CAS  PubMed  Google Scholar 

  4. Koyama A, Fujisaki M, Kobayashi M, et al. A glomerular permeability factor produced by human T cell hybridomas. Kidney Int. 1991;40:453–60.

    Article  CAS  PubMed  Google Scholar 

  5. Nephrotic syndrome in children. Prediction of histopathology from clinical and laboratory characteristics at time of diagnosis. A report of the International study of kidney disease in children. Kidney Int. 1978;13:159–65.

    Article  Google Scholar 

  6. Tarshish P, Tobin JN, Bernstein J, et al. Prognostic significance of the early course of minimal change nephrotic syndrome: report of the international study of kidney disease in children. J Am Soc Nephrol. 1997;8:769–76.

    CAS  PubMed  Google Scholar 

  7. Jeltsch A. Beyond Watson and Crick: DNA methylation and molecular enzymology of DNA methyltransferases. ChemBioChem. 2002;3:274–93.

    Article  CAS  PubMed  Google Scholar 

  8. Li E, Bestor TH, Jaenisch R. Targeted mutation of the DNA methyltransferase gene results in embryonic lethality. Cell. 1992;69:915–26.

    Article  CAS  PubMed  Google Scholar 

  9. Li E, Beard C, Jaenisch R. Role for DNA methylation in genomic imprinting. Nature. 1993;366:362–5.

    Article  CAS  PubMed  Google Scholar 

  10. Heard E, Clerc P, Avner P. X-chromosome inactivation in mammals. Annu Rev Genet. 1997;31:571–610.

    Article  CAS  PubMed  Google Scholar 

  11. Egger G, Liang G, Aparicio A, et al. Epigenetics in human disease and prospects for epigenetic therapy. Nature. 2004;429:457–63.

    Article  CAS  PubMed  Google Scholar 

  12. Audard V, Pawlak A, Candelier M, et al. Upregulation of nuclear factor-related kappa B suggests a disorder of transcriptional regulation in minimal change nephrotic syndrome. PLoS One. 2012;27:30523.

    Article  Google Scholar 

  13. Elie V, Fakhoury M, Deschenes G, et al. Physiopathology of idiopathic nephrotic syndrome: lessons from glucocorticoids and epigenetic perspectives. Pediatr Nephrol. 2012;27:1249–56.

    Article  PubMed  Google Scholar 

  14. Kobayashi Y, Aizawa A, Takizawa T, et al. DNA methylation changes between relapse and remission of minimal change nephrotic syndrome. Pediatr Nephrol. 2012;27:2233–41.

    Article  PubMed  PubMed Central  Google Scholar 

  15. The primary nephrotic syndrome in children. Identification of patients with minimal change nephrotic syndrome from initial response to prednisone. A report of the international study of kidney disease in children. J Pediatr. 1981;98:561–4.

    Article  Google Scholar 

  16. Rodenhiser DI, Andrews J, Kennette W, et al. Epigenetic mapping and functional analysis in a breast cancer metastasis model using whole-genome promoter tiling microarrays. Breast Cancer Res. 2008;10:R62.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Müller T, Schrötter A, Looße C, et al. Sense and nonsense of pathway analysis software in proteomics. J Proteome Res. 2011;10:5398–408.

    Article  PubMed  Google Scholar 

  18. Huang DW, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4:44–57.

    Article  CAS  Google Scholar 

  19. Hatada I, Fukasawa M, Kimura M, et al. Genome-wide profiling of promoter methylation in human. Oncogene. 2006;25:3059–64.

    Article  CAS  PubMed  Google Scholar 

  20. Hatada I, Morita S, Kimura M, et al. Genome-wide demethylation during neural differentiation of P19 embryonal carcinoma cells. J Hum Genet. 2008;53:185–91.

    Article  CAS  PubMed  Google Scholar 

  21. Sahali D, Pawlak A, Valanciute A, et al. A novel approach to investigation of the pathogenesis of active minimal-change nephrotic syndrome using subtracted cDNA library screening. J Am Soc Nephrol. 2002;13:1238–47.

    CAS  PubMed  Google Scholar 

  22. Valanciuté A, le Gouvello S, Solhonne B, et al. NF-kappa B p65 antagonizes IL-4 induction by c-maf in minimal change nephrotic syndrome. J Immunol. 2004;172:688–98.

    Article  PubMed  Google Scholar 

  23. Coppo R, Camilla R, Porcellini MG, et al. Saquinavir in steroid-dependent and -resistant nephrotic syndrome: a pilot study. Nephrol Dial Transplant. 2012;27:1902–10.

    Article  CAS  PubMed  Google Scholar 

  24. Okuyama S, Komatsuda A, Wakui H, et al. Up-regulation of TRAIL mRNA expression in peripheral blood mononuclear cells from patients with minimal-change nephrotic syndrome. Nephrol Dial Transpl. 2005;20:539–44.

    Article  CAS  Google Scholar 

  25. Bakr A, Shokeir M, El-Chenawi F, et al. Tumor necrosis factor-alpha production from mononuclear cells in nephrotic syndrome. Pediatr Nephrol. 2003;18:516–20.

    PubMed  Google Scholar 

  26. Suranyi MG, Guasch A, Hall BM, et al. Elevated levels of tumor necrosis factor-alpha in the nephrotic syndrome in humans. Am J Kidney Dis. 1993;21:251–9.

    Article  CAS  PubMed  Google Scholar 

  27. Santangelo S, Cousins DJ, Winkelmann NE, et al. DNA methylation changes at human Th2 cytokine genes coincide with DNase I hypersensitive site formation during CD4(+) T cell differentiation. J Immunol. 2002;169:1893–903.

    Article  CAS  PubMed  Google Scholar 

  28. Yamashita M, Ukai-Tadenuma M, Kimura M, et al. Identification of a conserved GATA3 response element upstream proximal from the interleukin-13 gene locus. J Biol Chem. 2002;277:42399–408.

    Article  CAS  PubMed  Google Scholar 

  29. Yap HK, Cheung W, Murugasu B, et al. Th1 and Th2 cytokine mRNA profiles in childhood nephrotic syndrome: evidence for increased IL-13 mRNA expression in relapse. J Am Soc Nephrol. 1999;10:529–37.

    CAS  PubMed  Google Scholar 

  30. Cho BS, Yoon SR, Jang JY, et al. Up-regulation of interleukin-4 and CD23/FcepsilonRII in minimal change nephrotic syndrome. Pediatr Nephrol. 1999;13:199–204.

    Article  CAS  PubMed  Google Scholar 

  31. Komatsuda A, Wakui H, Iwamoto K, et al. GATA-3 is upregulated in peripheral blood mononuclear cells from patients with minimal change nephrotic syndrome. Clin Nephrol. 2009;71:608–16.

    Article  CAS  PubMed  Google Scholar 

  32. Ikeuchi Y, Kobayashi Y, Arakawa H, et al. Polymorphisms in interleukin-4-related genes in patients with minimal change nephrotic syndrome. Pediatr Nephrol. 2009;24:489–95.

    Article  PubMed  Google Scholar 

  33. Sellier-Leclerc AL, Duval A, Riveron S, et al. A humanized mouse model of idiopathic nephrotic syndrome suggests a pathogenic role for immature cells. J Am Soc Nephrol. 2007;18:2732–9.

    Article  PubMed  Google Scholar 

  34. Kobayashi Y, Arakawa H, Suzuki M, et al. Polymorphisms of interleukin-4–related genes in Japanese children with minimal change nephrotic syndrome. Am J Kidney Dis. 2003;42:271–6.

    Article  CAS  PubMed  Google Scholar 

Download references

Authors’ contributions

YK, AA and TT designed the research; YK and AA performed the research; YK, AA and KI analyzed the data; YK and AA wrote the manuscript; TT, IH and HA gave comments on the manuscript and revised it. All authors read and approved the final manuscript.

Acknowledgements

The authors thank Ms. Tomoko Endo, Ms. Chinori Iijima, Ms. Sachiko Hayashi, and Ms. Kiyoe Ishii for their technical assistance, Dr. Chikage Yoshizawa, Dr. Hiromi Horiguchi, Dr. Yuka Ikeuchi, Dr. Satoko Kakegawa and Dr. Toshio Watanabe for their contribution in collecting patients’ samples.

Competing interests

The authors declare that they have no competing interests.

Availability of data and materials

All data generated or analysed during this study are included in this published article and its Additional files.

Consent for publication

Consent for publication has been obtained from all the participants/parents or legal guardian.

Ethics approval and consent to participate

This study was approved by the Ethics Committee of Gunma University Graduate School of Medicine, Japan (Receipt Number 89). Consent and permission to participate have been obtained from all the participants. Parental consent for participation of children under 16 years of age have been obtained. Participants and parents of children provided written consent.

Funding

This study was supported by JSPS KAKENHI Grant Numbers 19591238, 25430177, and 25670265; Kawano Masanori Memorial Foundation for Promotion of Pediatrics; the Joint Research Program of the Institute for Molecular and Cellular Regulation, Gunma University (Research Project Number: 11027); and Takeda Science Foundation.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Author information

Authors and Affiliations

  1. Department of Pediatrics, Gunma University Graduate School of Medicine, 3-39-22 Showa-machi, Maebashi, Gunma, 371-8511, Japan

    Yasuko Kobayashi, Akira Aizawa, Takumi Takizawa & Hirokazu Arakawa

  2. Life Science Tokyo Advanced Research Center (L-StaR), Hoshi University School of Pharmacy and Pharmaceutical Sciences, Tokyo, Japan

    Katsuhide Igarashi

  3. Laboratory of Genome Science, Biosignal Genome Resource Center, Gunma University Institute for Molecular and Cellular Regulation, Maebashi, Japan

    Izuho Hatada

  4. Center for Fiber and Textile Science, Kyoto Institute of Technology, Matsugasaki, Kyoto, Japan

    Akira Aizawa

Authors
  1. Yasuko Kobayashi
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Akira Aizawa
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Takumi Takizawa
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Katsuhide Igarashi
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Izuho Hatada
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Hirokazu Arakawa
    View author publications

    You can also search for this author in PubMed Google Scholar

Corresponding author

Correspondence to Yasuko Kobayashi.

Additional files

Additional file 1. MIAMI results of 85 co-occurring genes in Th0s identified in all 3 comparisons.

13104_2017_2719_MOESM2_ESM.xlsx

Additional file 2. MIAMI results of 206 co-occurring genes in Th0s identified in the 2 comparisons; relapse versus remission and relapse versus controls (proteinuria-related genes).

13104_2017_2719_MOESM3_ESM.xlsx

Additional file 3. MIAMI results of 58 co-occurring genes in Th0s identified in the 2 comparisons; relapse versus controls and remission versus controls (disease-related genes).

13104_2017_2719_MOESM4_ESM.docx

Additional file 4. Table S4: Analysis with DAVID bioinformatics resources: Analysis of the 206 proteinuria oriented co-occurring genes. Table S5: Analysis with DAVID bioinformatics resources: Analysis of the 58 disease-related co-occurring genes.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Kobayashi, Y., Aizawa, A., Takizawa, T. et al. Changes in DNA methylation in naïve T helper cells regulate the pathophysiological state in minimal-change nephrotic syndrome. BMC Res Notes 10, 480 (2017). https://doi.org/10.1186/s13104-017-2719-1

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13104-017-2719-1

Keywords

BMC Research Notes

ISSN: 1756-0500

Contact us