Skip to main content
Download PDF

Comprehensive analysis of aberrantly expressed profiles of mRNA and its relationship with serum galactose-deficient IgA1 level in IgA nephropathy

Journal of Translational Medicine volume 17, Article number: 320 (2019) Cite this article

Abstract

Background

Immunoglobulin A nephropathy (IgAN) is the leading cause of end-stage kidney disease. Previous mRNA microarray profiling studies of IgAN revealed inconsistent data. We sought to identify the aberrantly expressed genes and biological pathways by integrating IgAN gene expression datasets in blood cells and performing systematically experimental validation. We also explored the relationship between target genes and galactose-deficient IgA1 (Gd-IgA1) in IgAN.

Methods

We retrieved Gene Expression Omnibus (GEO) datasets of IgAN. Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were used for functional analysis. Deep sequencing on RNA isolated from B cells was used for microarray validation. The relationship between target mRNA expressions and Gd-IgA1 levels in serum were also studied.

Results

Three studies with microarray expression profiling datasets met our inclusion criteria. We identified 655 dyregulated genes, including 319 up-regulated and 336 down-regulated genes in three GEO datasets with a total of 35 patients of IgAN and 19 healthy controls. Based on biological process in GO term, these dyregulated genes are mainly related to pentose-phosphate shunt, non-oxidative branch, post-embryonic camera-type eye development and leukocyte activation. KEGG pathway analysis of microarray data revealed that these aberrantly expressed genes were enriched in human T-cell leukemia virus 1 infection, proteoglycans in cancer, intestinal immune network for IgA production and autophagy. We further performed deep sequencing on mRNAs isolated from B cells of an independent set of five patients with IgAN and three healthy persons with the same clinical and demographic characteristics. Seventy-seven genes overlapped with 655 differentially regulated genes mentioned above, including 43 up-regulated and thirty-four down-regulated genes. We next investigated whether these genes expression correlated with Gd-IgA1 levels in IgAN patients. Pearson correlation analyses showed PTEN (phosphatase and tensin homolog) was the most powerful gene negatively correlated with Gd-IgA1 levels.

Conclusions

These results demonstrated that dyregulated genes in patients with IgAN were enriched in intestinal immune network for IgA production and autophagy process, and PTEN in B cells might be involved in the mechanism of Gd-IgA1 production.

Background

Immunoglobulin A nephropathy (IgAN), also known as Berger’s disease, is the most common glomerulonephritis in the world [1], approximately 20% of such patients will progress to end-stage kidney disease (ESKD) within 10 years of diagnosis [2]. The initial event in the pathogenesis of IgAN is the production of galactose-deficient IgA1 (Gd-IgA1) [3], which lead to IgA1 self-aggregation, IgA1–IgG immune complex formation, and sequential deposition in the kidney [4, 5]. To date, the cause of IgAN remains to be determined.

Array-based gene-expression analysis is a useful tool to evaluate and identify differently expressed genes between pathological and control samples and is able to provide an unbiased screening and assessment of gene expression. Through microarray profiles, IgAN has been reported to be interrelated with multiple genes, including IFI27, CXCL1 and STAT3 [6,7,8]. However, recent studies showed that the data from one microarray analysis were inconsistent with other studies. The discrepancy may be due to small sample sizes and varying results obtained from different microarray platforms. Comprehensive gene set enrichment analysis could help us overcome these limitations and control such confounding factors by increasing statistical power, leading to more robust, reproducible and accurate predictions for identifying differently expressed genes [9, 10]. Indeed, such studies have been successfully used to identify gene signatures in many diseases [11,12,13].

In this study, we first identified differentially expressed genes and biological pathways in blood cells involved in IgAN using the chip data in Gene Expression Omnibus (GEO) database. Then deep sequencing on RNA isolated from B cells was used for microarray validation. We also explored the relationship between target genes and Gd-IgA1 in IgAN.

Materials and methods

Sample collection

A total of 25 IgAN patients diagnosed in Tianjin Medical University General Hospital from May to August 2017, were enrolled in this study. The diagnosis was based on the deposition of IgA in the glomerular mesangium by immunofluorescence detection, as well as the lack of clinical or serological evidence of other inflammatory conditions, such as Henoch Schoenlein purpura. At the same time, 21 healthy volunteers whose age and gender matched with patients were recruited. 5 IgAN patients and three healthy controls were included in the RNA deep sequencing experiment, 20 IgAN patients and 18 healthy controls were used for RNA deep sequencing validation. Written informed consent was obtained from each patient and healthy participant.

Identification of individual studies and data processing

We performed an electronic search using the keywords ‘‘IgA nephropathy, IgA nephritis, or Berger’s disease’’ to identify studies that used microarray techniques through October 2018. We obtained mRNA expression data in GEO without language or publication date restrictions. Studies were included if they met the following criteria: (1) patients were diagnosed with IgAN based on IgAN diagnostic criteria [14]; (2) case–control studies; (3) all datasets were genome-wide; (4) all samples were from blood cells; (5) the number of cases and controls in each dataset must be ≥ 2; and (6) complete microarray raw data were available. We excluded any animal or duplicated studies. The whole genome raw expression data of included studies were downloaded from the National Center for Biotechnology Information GEO (http://www.ncbi.nlm.nih.gov/geo/). We used the Robust Multichip Average (RMA) algorithm in oligo package to normalize the raw expression data and generate normalized gene expression intensity [15]. After merging the data from different series, ComBat was used to adjust the known variances from different batches using an empirical Bayesian framework [16].

Differential expression genes analysis

Differential expression gene analysis was performed using R v3.2.2. The tool based on a t-test was used to detect differentially expressed genes between IgAN patients and controls. The log2 transformation was used to obtain the standardized expression values [17]. Significantly up-regulated genes were defined by as a logarithmic transformed fold-change (FC) > 0.26 and P value ≤ 0.05. Significantly down-regulated genes were defined by a logFC ≤ 0.26 and P value ≤ 0.05.

Functional analysis of differentially expressed genes

We next used Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses to interpret biological significance of the differentially expressed genes. The gene-associated GO terms were used according to http://www.geneontology.org/ [18], while KEGG categories organized the differentially expressed genes into gene pathways (http://www.kegg.jp/) [19]. The functional analysis was performed using the using R v3.2.2.

Assay for Gd-IgA1

Galactose-deficient IgA1 levels in serum were determined by ELISA. Serum Gd-IgA1 levels were detected using KM55 ELISA kit according to the manufacturer’s specifications (IBL, Japan) [20]. Serum samples were diluted in proportions of 1:50 with EIA buffer and incubation for 60 min with plate lid. Then washing four times with wash buffer, prepared labeled antibody was incubated for 30 min. Plates were washed and added 50 μL TMA solution incubation for 30 min in dark. At last, the color reaction was stopped and the absorbance was measured at 450/625 nm with an EL312 Bio-Kinetics microplate reader (Bio-Tek Instruments, Winooski, VT).

B lymphocytes isolation

About 5 mL venous blood sample was taken into ethylenediaminetetraacetic acid (EDTA) anticoagulated tubes. Peripheral blood mononuclear cells (PBMCs) were separated by density-gradient centrifugation on Ficoll (TBD, China), then washed three times with phosphate buffered saline (PBS) and resuspended in PBS + 1% bovine serum albumin (BSA). Peripheral B lymphocytes were isolated using CD19 positive magnetic beads (Miltenyi Biotec, USA) according to the manufacturer’s instructions.

RNA extraction and real time reverse transcription PCR (RT-PCR)

Total cellular RNA was extracted from CD19 positive B lymphocytes using the TRIZOL Reagent (Invitrogen, USA). RNA quantity was determined using NanoDrop ND-1000 spectrophotometer. cDNA was synthesized using 300 ng total RNA with revert first-strand cDNA kit according to manufacturer’s protocol (Promega, USA). Resulting cDNA was amplified with a 20 µL reaction mixture using SYBR Green PCR Master Mix (Roche, USA) in an Applied Biosystem 7500 Real-Time PCR System. And the primer pairs of PTEN and GAPDH were listed in Additional file 1: Table S1. The fold change between patients and controls was expressed by the 2−△△CT method. The GAPDH gene amplification was used as a reference standard to normalize the target signal.

RNA deep sequencing

RNA samples were sent to Shanghai Gminix, Biotechnology Co, Ltd. (Shanghai, China) for RNA deep sequencing on an Illumina HiSeq X10 sequencing platform with 15G PF data (Illumina, San Diego, CA). Procedures were performed as described in detail on the website of Gminix (http://www.gminix.com). The data discussed in this publication have been deposited in NCBI sequence read archive (SRA), and are accessible through SRA Series Accession number SUB6070993.

Statistical analysis

The differences in gene expression levels between groups were compared using ANOVA and P value < 0.05 was considered statistically significant. Statistical analysis was performed using SPSS 17.0 software.

Results

Characteristics of included studies

We obtained 16 potential relevant publications and retrieved eleven studies for full text view following the initial database search. Three studies with microarray expression profiling data sets met our inclusion criteria (Fig. 1) [7, 21, 22]. The final analysis input consisted of 35 patients of IgAN and 19 healthy controls. Details of these selected studies are summarized in Table 1, including GEO accession, microarray platform, patients’ age, gender, and baseline levels of serum creatinine and proteinuria. These studies were published between 2010 and 2016 and sample size ranged from 17 to 20.

Fig. 1
figure 1

Flowchart for selection of eligible studies

Full size image
Table 1 Characteristics of the individual studies included in the meta-analysis
Full size table

Differential gene expressions in IgAN patients

A total of 655 genes were identified differentially expressed in patients of IgAN group compared with controls across three microarray datasets (P < 0.05). Among these candidate implantation-associated genes, 319 were up-regulated and 336 were down-regulated in blood cells. A heat map visualization of the mRNA expression profile in IgAN was displayed in Additional file 2: Figure S1. We listed the top 10 most significantly up- or down-regulated genes in Table 2. The top 100 differentially expressed genes in GEO microarray were listed in Additional file 1: Table S2.

Table 2 The ten differentially expressed genes in IgAN patients and healthy control
Full size table

Function analysis

To analyze the related functions of the dyregulated genes, GO enrichment and KEGG pathway analysis were performed. After GO enrichment analysis, the top 10 most significant GO terms in biological process (BP), cellular component (CC), and molecular function (MF) were depicted in Fig. 2. Based on the BP, these dyregulated genes are mainly related to pentose-phosphate shunt, non-oxidative branch (GO: 0009052, four genes were enriched with P = 2.22 × 10−5), post-embryonic camera-type eye development (GO: 0031077, four genes were enriched with P = 5.05 × 10−5) and leukocyte activation (GO: 0045321, 69 genes were enriched with P = 6.60 × 10−5). For MF, the genes were mainly associated with protein kinase binding (GO: 0019901, 38 genes were enriched with P = 1.43 × 10−4), kinase binding (GO: 0019900, 41 genes were enriched with P = 2.75 × 10−4) and C-X3-C chemokine binding (GO: 0019960, three genes were enriched with P = 4.19 × 10−4). For CC, the genes were associated with mitochondrial matrix (GO: 0005759, 32 genes were enriched with P = 7.03 × 10−5), membrane protein complex (GO: 0098796, 42 genes were enriched with P = 7.29 × 10−5), and focal adhesion (GO: 0005925, 29 genes were enriched with P = 7.73 × 10−5).

Fig. 2
figure 2

The top 10 most significantly enriched GO terms identified in IgAN

Full size image

Figure 3 shows the top 11 significantly enriched pathways in blood cells. The KEGG enrichment pathway analysis of our microarray data revealed that these dyregulated genes were enriched in human T-cell leukemia virus 1 infection (hsa05166, 20 genes were enriched with P = 0.005), proteoglycans in cancer (hsa05205, 16 genes were enriched with P = 0.01), intestinal immune network for IgA production (hsa04672, six genes were enriched with P = 0.01), autophagy-animal (hsa04140, 11 genes were enriched with P = 0.01) and mitophagy-animal (hsa04137, 7 genes were enriched with P = 0.01).

Fig. 3
figure 3

The top 10 most significant KEGG pathways identified in IgAN

Full size image

Validation of selected targets

To further validate some of the results obtained, we performed deep sequencing on RNAs isolated from B cells of an independent set of five patients with IgAN and three healthy persons with the same clinical and demographic characteristics as those in the population used for microarray experiments (Table 3). A heat map visualization of the mRNA expression profile in IgAN was displayed in Fig. 4. The top 100 differential expression genes could be seen in Additional file 1: Table S3. We confirmed 77 differentially regulated genes, including 43 up-regulated and 34 down-regulated genes were overlapped with just identified 655 genes by microarray analysis. We listed 10 significantly up- or down-regulated genes in Table 4.

Table 3 Characteristics of individuals included for deep sequencing
Full size table
Fig. 4
figure 4

Heatmap illustration of the differentially expressed genes in IgAN. Expression levels are represented by red (high) and green (low expression). Samples are from five patients with IgAN (pink) and three healthy subjects (blue)

Full size image
Table 4 Ten differentially expressed genes confirmed by deep sequencing
Full size table

Correlation between dyregulated genes expression and Gd-IgA1 levels

After the confirmation that IgAN patients had abnormal expression of above 77 genes, we next investigated whether these genes expression was correlated with Gd-IgA1 levels in IgAN patients. Gd-IgA1 levels, evaluated by KM55 antibody, were significantly higher in serum of patients with IgAN than in those of healthy persons (P < 0.001; Fig. 5). The results identified eight genes whose expression was correlated with Gd-IgA1 levels (Additional file 1: Table 4). Pearson correlation analyses showed PTEN was the most powerful gene negatively correlated with Gd-IgA1 levels (r = − 0.83, P = 0.01, Fig. 6), suggesting that downregulated expression of PTEN might be involved in the mechanism of increased Gd-IgA1 production in IgAN.

Fig. 5
figure 5

Serum Gd-IgA1 levels in IgAN patients and healthy control

Full size image
Fig. 6
figure 6

The correlation between PTEN and Gd-IgA1

Full size image

PTEN mRNA levels in IgAN patients

PTEN takes the 21th place in the top hits overall in microarray and 320th place in B-cell RNA-seq datasets. We next isolated CD19 + B lymphocytes from 20 IgAN patients to verify the differently expression of PTEN mRNA. As shown in Fig. 7, the PTEN mRNA level of patients with IgAN is significantly lower than that of healthy control (1.12 ± 0.47 vs. 1.54 ± 0.68, P = 0.03).

Fig. 7
figure 7

The PTEN mRNA level in 20 patients with IgAN and 18 healthy control

Full size image

Discussion

It is the first study that provides an overview of three gene expression arrays in blood cells in IgAN and eventually discovers genes and pathways potentially involved in IgAN pathogenesis. Deep sequencing on RNA further verified these dyregulated genes in B cells in IgAN patients compared to healthy control samples, supporting the validity of our results.

In the era of high-throughput techniques, integration of heterogeneous data of microarray holds the promise of providing new insights of pathogenesis of IgAN. Functional analysis helps to understand differentially expressed genes from the point of view of canonical prior knowledge structured in the form of pathways diagrams. Our study provided more insightful biological interpretation by performing GO enrichment and KEGG pathway analyses in IgAN patients and healthy persons. We found that human T-cell leukemia virus 1 infection, proteoglycans in cancer, intestinal immune network for IgA production and autophagy were involved in IgAN. Large international GWASs have identified several susceptibility loci involved in dynamic environment of the intestinal mucosa and abnormal mucosal production of IgA against microbial antigens associated with susceptibility to IgAN [23, 24]. When the significant loci were tested for enrichment in KEGG pathways, Krzysztof showed the top overrepresented pathway was “Intestinal Immune Network for IgA Production”, which played a central role in the disease pathogenesis [25]. These data are consistent with the clinical observation that intestinal mucosal infections frequently trigger episodes of IgAN [26]. The recent random controlled trial, the NEFIGAN study, was designed to evaluate the effectiveness of Nefecon, an oral formulation that releases glucocorticosteroid budesonide in the lower ileum and ascending colon [27]. The result showed that Nefecon significantly decreased the level of proteinuria and maintained stable kidney function in patients with IgAN. These results clearly linked intestinal mucosal inflammatory disorders with the risk of IgAN. Evidences suggested that the agents, including blisibimod and atacicept, which target increasing levels of circulating IgA, may provide more specific therapy than what is currently available.

Autophagy is a highly conserved process that degrades cellular long-lived proteins and organelles. Accumulating evidence indicates that autophagy plays a critical role in chronic kidney disease. Recent GWAS have suggested that MTMR3 (encoding myotubularin-related phosphatase 3), participated in autophagy, is a susceptibility genes for IgAN [24]. And lower MTMR3 level was observed in IgAN patients in genome-wide gene expression analyses [28]. These data suggested autophagy may play an important role in the pathogenesis of IgAN.

Many genes tend to be expressed differentially in IgAN and the challenge is to identify the important genes and pathways associated with the disease. Characterizing the molecular and cellular events during the pathogenesis is an important endeavor. Recent studies indicated that Gd-IgA1 was the trigger factor of IgAN. After the confirmation 77 genes aberrantly expressed in IgAN by RNA sequence, we next investigated whether these genes expression was correlated with Gd-IgA1 levels in IgAN patients. The results identified eight genes whose expression was correlated with Gd-IgA1 levels. Pearson correlation analyses showed PTEN was the most powerful gene negatively correlated with Gd-IgA1 levels. PTEN (phosphatase and tensin homolog), has lipid/tyrosine phosphatase-dependent and -independent activities that alter cell behavior. Cox et al. [21] suggested that PTEN participated in the pathogenesis of IgAN. In line with our results, a previous study conducted on PBMCs of IgAN patients showed patients with IgAN exhibited lower PTEN expression [21, 29]. The precise function of this gene in production of Gd-IgA1 is still unclear. Grazia et al. demonstrated miR-148b significantly upregulated in patients with IgAN, and its putative target genes C1GALT1, INVS, and PTEN were down-regulated in these patients [30]. C1GALT1, encoding β1,3-galactosyltransferase, plays a critical role in the glycosylation process of IgA1 in IgAN patients [31, 32]. Whether there is a protein–protein interaction between C1GALT1 and PTEN or an abnormal upstream regulation mechanism is still unknown, future studies are needed to investigate the underlying molecular mechanisms of the gene alterations and develop novel strategies by targeting some of the identified genes/pathways as therapeutic tools to control IgAN.

Although many common susceptibility genes with IgAN were found by two methods, there were still several differential results between those previous microarray-analyses (discussed above) and our current RNA-seq analysis. Many top differentially expressed genes in IgAN patients and healthy control were not validated. The prevalence of IgAN shows differences among geographic and racial populations. The reported prevalence of IgAN is higher in Asia than in Europe and America. Association with the angiotensin-converting enzyme genes was observed in Asian but not in Western patients. Actually, most genetic association studies based on a candidate-gene approach have not been replicated between different races. Besides, genetic heterogeneity between ethnic groups, methodological problems (e.g., small sample size, different detection method) might be another reason contributing to this difference.

Despite being comprehensive, our analysis has several limitations. First, given that the extracted GEO datasets were based on comparisons of IgAN versus healthy controls, the predicted targets might not be specific to IgAN, but reflect general kidney disease. Second, the experiments took in whole-blood samples for gene-expression analyses were included, while those conducted in kidney biopsy samples were excluded. As we all know, Gd-IgA1 is mainly produced in B cells, which represents the first hit of the pathogenesis of IgAN. And in consideration of different transcriptase across tissues, we seek to identify the most robust different expressed genes in a blood cell transcriptome. Although we restricted input to blood cells datasets, some confounding factors (for example differences in blood cells types, chip platform, cohort sizes or patients’ ethnicity) are not taken into account, which collectively may affect the analysis output. Third, it is unclear how far the results of our analysis may be generalized to the world population of patients with IgAN, there were only eight Chinese subjects included for experimental validation. Further studies with larger sample sizes are needed to identify target mRNA and assess its relationship with serum galactose-deficient IgA1 level in IgAN.

Conclusion

In conclusion, these results demonstrated that dyregulated genes in blood cells of patients with IgAN are enriched in intestinal immune network for IgA production and autophagy. PTEN might be involved in the mechanism of increased Gd-IgA1 production in IgAN. Such altered genes and pathways could helpfully expand our knowledge of IgAN pathogenesis.

Availability of data and materials

Raw data used during the current study are available from the corresponding author on reasonable request for non-commercial use.

Abbreviations

IgAN:

immunoglobulin A nephropathy

GEO:

Gene Expression Omnibus

GO:

Gene Ontology

KEGG:

Kyoto Encyclopedia of Genes and Genomes

Gd-IgA1:

galactose-deficient IgA1

ESKD:

end-stage kidney disease

PTEN:

phosphatase and tensin homolog

References

  1. D’Amico G. The commonest glomerulonephritis in the world: IgA nephropathy. Q J Med. 1987;64:709–27.

    PubMed  Google Scholar 

  2. Li X, Liu Y, Lv J, Shi S, Liu L, Chen Y, Zhang H. Progression of iga nephropathy under current therapy regimen in a Chinese population. Clin J Am Soc Nephrol. 2014;9:484–9.

    Article  PubMed Central  PubMed  Google Scholar 

  3. Coppo R, Amore A. Aberrant glycosylation in iga nephropathy (IgAN). Kidney Int. 2004;65:1544–7.

    Article  CAS  PubMed  Google Scholar 

  4. Novak J, Rizk D, Takahashi K, Zhang X, Bian Q, Ueda H, Ueda Y, Reily C, Lai LY, Hao C, Novak L, Huang ZQ, Renfrow MB, Suzuki H, Julian BA. New insights into the pathogenesis of IgA nephropathy. Kidney Dis. 2015;1:8–18.

    Article  Google Scholar 

  5. Rodrigues JC, Haas M, Reich HN. IgA nephropathy. Clin J Am Soc Nephrol. 2017;12:677–86.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  6. Zhao Y, Zhu L, Zhou T, Zhang Q, Shi S, Liu L, Lv J, Zhang H. Urinary cxcl1: a novel predictor of IgA nephropathy progression. PLoS ONE. 2015;10:e0119033.

    Article  PubMed Central  PubMed  Google Scholar 

  7. Nagasawa Y, Okuzaki D, Muso E, Yamamoto R, Shinzawa M, Iwasaki Y, Iwatani H, Nakanishi T, Isaka Y, Nojima H. Ifi27 is a useful genetic marker for diagnosis of immunoglobulin a nephropathy and membranous nephropathy using peripheral blood. PLoS ONE. 2016;11:e0153252.

    Article  PubMed Central  PubMed  Google Scholar 

  8. Yamada K, Huang ZQ, Raska M, Reily C, Anderson JC, Suzuki H, Ueda H, Moldoveanu Z, Kiryluk K, Suzuki Y, Wyatt RJ, Tomino Y, Gharavi AG, Weinmann A, Julian BA, Willey CD, Novak J. Inhibition of stat3 signaling reduces IgA1 autoantigen production in IgA nephropathy. Kidney Int Rep. 2017;2:1194–207.

    Article  PubMed Central  PubMed  Google Scholar 

  9. Tseng GC, Ghosh D, Feingold E. Comprehensive literature review and statistical considerations for microarray meta-analysis. Nucleic Acids Res. 2012;40:3785–99.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  10. Rung J, Brazma A. Reuse of public genome-wide gene expression data. Nat Rev Genet. 2013;14:89–99.

    Article  CAS  PubMed  Google Scholar 

  11. Kohane IS, Valtchinov VI. Quantifying the white blood cell transcriptome as an accessible window to the multiorgan transcriptome. Bioinformatics. 2012;28:538–45.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  12. Wang Z, Chen X, Zhang D, Cao Y, Zhang L, Tang W. Pycard gene plays a key role in rapidly progressive glomerulonephritis: results of a weighted gene co-expression network analysis. Am J Nephrol. 2018;48:193–204.

    Article  CAS  PubMed  Google Scholar 

  13. Wang YZ, Zhu DY, Xie XM, Ding M, Wang YL, Sun LL, Zhang N, Shen E, Wang XX. Ea15, mir22, linc00472 as diagnostic markers for diabetic kidney disease. J Cell Physiol. 2018;234:8797–803.

    Article  PubMed  Google Scholar 

  14. Li GS, Zhang H, Lv JC, Shen Y, Wang HY. Variants of C1GALT1 gene are associated with the genetic susceptibility to iga nephropathy. Kidney Int. 2007;71:448–53.

    Article  CAS  PubMed  Google Scholar 

  15. Carvalho BS, Irizarry RA. A framework for oligonucleotide microarray preprocessing. Bioinformatics. 2010;26:2363–7.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  16. Johnson WE, Li C, Rabinovic A. Adjusting batch effects in microarray expression data using empirical Bayes methods. Biostatistics. 2007;8:118–27.

    Article  PubMed  Google Scholar 

  17. Gautier L, Cope L, Bolstad BM, Irizarry RA. Affy—analysis of affymetrix genechip data at the probe level. Bioinformatics. 2004;20:307–15.

    Article  CAS  PubMed  Google Scholar 

  18. Falcon S, Gentleman R. Using GOstats to test gene lists for GO term association. Bioinformatics. 2007;23:257–8.

    Article  CAS  PubMed  Google Scholar 

  19. Kanehisa M, Furumichi M, Tanabe M, Sato Y, Morishima K. KEGG: new perspectives on genomes, pathways, diseases and drugs. Nucleic Acids Res. 2017;45:D353–61.

    Article  CAS  PubMed  Google Scholar 

  20. Yasutake J, Suzuki Y, Suzuki H, Hiura N, Yanagawa H, Makita Y, Kaneko E, Tomino Y. Novel lectin-independent approach to detect galactose-deficient IgA1 in IgA nephropathy. Nephrol Dial Transplant. 2015;30:1315–21.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  21. Cox SN, Sallustio F, Serino G, Pontrelli P, Verrienti R, Pesce F, Torres DD, Ancona N, Stifanelli P, Zaza G, Schena FP. Altered modulation of WNT-beta-catenin and PI3K/AKT pathways in IgA nephropathy. Kidney Int. 2010;78:396–407.

    Article  CAS  PubMed  Google Scholar 

  22. Cox SN, Serino G, Sallustio F, Blasi A, Rossini M, Pesce F, Schena FP. Altered monocyte expression and expansion of non-classical monocyte subset in IgA nephropathy patients. Nephrol Dial Transplant. 2015;30:1122–232.

    Article  CAS  PubMed  Google Scholar 

  23. Yu XQ, Li M, Zhang H, Low HQ, Wei X, Wang JQ, Sun LD, Sim KS, Li Y, Foo JN, Wang W, Li ZJ, Yin XY, Tang XQ, Fan L, Chen J, Li RS, Wan JX, Liu ZS, Lou TQ, Zhu L, Huang XJ, Zhang XJ, Liu ZH, Liu JJ. A genome-wide association study in han chinese identifies multiple susceptibility loci for IgA nephropathy. Nat Genet. 2011;44:178–82.

    Article  PubMed  Google Scholar 

  24. Gharavi AG, Kiryluk K, Choi M, Li Y, Hou P, Xie J, Sanna-Cherchi S, Men CJ, Julian BA, Wyatt RJ, Novak J, He JC, Wang H, Lv J, Zhu L, Wang W, Wang Z, Yasuno K, Gunel M, Mane S, Umlauf S, Tikhonova I, Beerman I, Savoldi S, Magistroni R, Ghiggeri GM, Bodria M, Lugani F, Ravani P, Ponticelli C, Allegri L, Boscutti G, Frasca G, Amore A, Peruzzi L, Coppo R, Izzi C, Viola BF, Prati E, Salvadori M, Mignani R, Gesualdo L, Bertinetto F, Mesiano P, Amoroso A, Scolari F, Chen N, Zhang H, Lifton RP. Genome-wide association study identifies susceptibility loci for IgA nephropathy. Nat Genet. 2011;43:321–7.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  25. Kiryluk K, Li Y, Scolari F, Sanna-Cherchi S, Choi M, Verbitsky M, Fasel D, Lata S, Prakash S, Shapiro S, Fischman C, Snyder HJ, Appel G, Izzi C, Viola BF, Dallera N, Del Vecchio L, Barlassina C, Salvi E, Bertinetto FE, Amoroso A, Savoldi S, Rocchietti M, Amore A, Peruzzi L, Coppo R, Salvadori M, Ravani P, Magistroni R, Ghiggeri GM, Caridi G, Bodria M, Lugani F, Allegri L, Delsante M, Maiorana M, Magnano A, Frasca G, Boer E, Boscutti G, Ponticelli C, Mignani R, Marcantoni C, Di Landro D, Santoro D, Pani A, Polci R, Feriozzi S, Chicca S, Galliani M, Gigante M, Gesualdo L, Zamboli P, Battaglia GG, Garozzo M, Maixnerova D, Tesar V, Eitner F, Rauen T, Floege J, Kovacs T, Nagy J, Mucha K, Paczek L, Zaniew M, Mizerska-Wasiak M, Roszkowska-Blaim M, Pawlaczyk K, Gale D, Barratt J, Thibaudin L, Berthoux F, Canaud G, Boland A, Metzger M, Panzer U, Suzuki H, Goto S, Narita I, Caliskan Y, Xie J, Hou P, Chen N, Zhang H, Wyatt RJ, Novak J, Julian BA, Feehally J, Stengel B, Cusi D, Lifton RP, Gharavi AG. Discovery of new risk loci for IgA nephropathy implicates genes involved in immunity against intestinal pathogens. Nat Genet. 2014;46:1187–96.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  26. Wyatt RJ, Julian BA. IgA nephropathy. N Engl J Med. 2013;368:2402–14.

    Article  CAS  PubMed  Google Scholar 

  27. Fellstrom BC, Barratt J, Cook H, Coppo R, Feehally J, de Fijter JW, Floege J, Hetzel G, Jardine AG, Locatelli F, Maes BD, Mercer A, Ortiz F, Praga M, Sorensen SS, Tesar V, Del Vecchio L, Investigators NT. Targeted-release budesonide versus placebo in patients with IgA nephropathy (NEFIGAN): a double-blind, randomised, placebo-controlled phase 2b trial. Lancet. 2017;389:2117–27.

    Article  PubMed  Google Scholar 

  28. Zhou XJ, Nath SK, Qi YY, Cheng FJ, Yang HZ, Zhang Y, Yang W, Ma JY, Zhao MH, Shen N, Zhang H. Brief report: identification of MTMR3 as a novel susceptibility gene for lupus nephritis in northern han chinese by shared-gene analysis with IgA nephropathy. Arthritis Rheumatol. 2014;66:2842–8.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  29. Hu S, Bao H, Xu X, Zhou X, Qin W, Zeng C, Liu Z. Increased miR-374b promotes cell proliferation and the production of aberrant glycosylated IgA1 in B cells of IgA nephropathy. FEBS Lett. 2015;589:4019–25.

    Article  CAS  PubMed  Google Scholar 

  30. Serino G, Sallustio F, Cox SN, Pesce F, Schena FP. Abnormal miR-148b expression promotes aberrant glycosylation of IgA1 in IgA nephropathy. J Am Soc Nephrol. 2012;23:814–24.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  31. Gale DP, Molyneux K, Wimbury D, Higgins P, Levine AP, Caplin B, Ferlin A, Yin P, Nelson CP, Stanescu H, Samani NJ, Kleta R, Yu X, Barratt J. Galactosylation of IgA1 is associated with common variation in C1GALT1. J Am Soc Nephrol. 2017;28:2158–66.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  32. Pirulli D, Crovella S, Ulivi S, Zadro C, Bertok S, Rendine S, Scolari F, Foramitti M, Ravani P, Roccatello D, Savoldi S, Cerullo G, Lanzilotta SG, Bisceglia L, Zelante L, Floege J, Alexopoulos E, Kirmizis D, Ghiggeri GM, Frasca G, Schena FP, Amoroso A, European Ig ANC. Genetic variant of C1GALT1 contributes to the susceptibility to IgA nephropathy. J Nephrol. 2009;22:152–9.

    CAS  PubMed  Google Scholar 

Download references

Acknowledgements

The authors thank all the study subjects for their participation.

Funding

This study is supported by the National Natural Science Foundation (81600553, 81600554).

Author information

Authors and Affiliations

  1. Department of Nephrology, Tianjin Medical University General Hospital, NO. 154, Anshan Road, Heping District, Tianjin, People’s Republic of China

    Youxia Liu, Junya Jia & Tiekun Yan

  2. Department of Nephrology, The Second Hospital of Shandong University, Jinan, People’s Republic of China

    Xiangchun Liu

  3. Radiology Department, Tianjin Medical University General Hospital, Tianjin, People’s Republic of China

    Jie Zheng

Authors
  1. Youxia Liu
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Xiangchun Liu
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Junya Jia
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Jie Zheng
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Tiekun Yan
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

JZ prepared and collected samples. YL and XL performed all experiments and analyzed data. YL drafted the manuscript. JJ and TY conceived of the overall study. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Youxia Liu or Tiekun Yan.

Ethics declarations

Ethics approval and consent to participate

All subjects provided written informed consents. The study protocol was approved by the Institutional Ethical Committee of Tianjin medical university general hospital.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

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

Supplementary information

Additional file 1: Table S1.

Primers used to amplify the PTEN and GAPDH genes. Table S2. The top 100 differentially expressed genes in GEO microarray. Table S3. The top 100 differentially expressed genes by RNA-deep sequencing. Table S4. Correlation between 8 dyregulated genes expression and Gd-IgA1 levels.

Additional file 2: Figure S1.

Heatmap illustration of the patterns of change in a particular gene across different datasets. Expression levels are represented by red (high) and green (low expression). Samples are from 35 patients with IgAN (red) and 19 healthy subjects (green).

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

Liu, Y., Liu, X., Jia, J. et al. Comprehensive analysis of aberrantly expressed profiles of mRNA and its relationship with serum galactose-deficient IgA1 level in IgA nephropathy. J Transl Med 17, 320 (2019). https://doi.org/10.1186/s12967-019-2064-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12967-019-2064-3

Keywords

Journal of Translational Medicine

ISSN: 1479-5876

Contact us