Skip to main content

Association of gut microbiota with idiopathic membranous nephropathy



The prevalence of idiopathic membranous nephropathy (IMN) is increasing worldwide and the gut microbiota is recognized to play a role in its pathology. The aim of this study was to understand the involvement of the gut–kidney axis in IMN by analyzing the composition of the gut microbiota of biopsy-proven IMN patients compared with healthy controls (HC).


Fecal samples from 30 patients with IMN diagnosed by renal biopsy and 30 healthy co-residents (control group) were collected for analysis in the Nephrology Department of the Second Affiliated Hospital of Harbin Medical University. The microbiota composition was analyzed by a 16S rRNA microbial profiling approach.


The results indicated that the α- and β-diversity of IMN patients differed significantly from those of the HC groups (P < 0.05). At the phylum level, IMN patients showed an increased abundance of Proteobacteria but a reduced abundance of Bacteroidota compared with the HC group. Actinobacteriota abundance showed a strong negative correlation with the estimated glomerular filtration rate. At the genus level, Faecalibacterium, Agathobacter, and Bacteroides were less abundant in the IMN group than in the HC group (LDA score > 2). Abundant bacterial functions related to lipid metabolism were observed among IMN group.


Patients with IMN appear to have an altered gut microbiome, which could provide reference for future research on the interaction mechanism between the intestinal flora and IMN.

Peer Review reports


Membranous nephropathy (MN), the main pathological type of adult nephrotic syndrome (NS), is an organ-specific autoimmune disease. The pathological characteristics of MN are diffuse deposition of glomerular basement membrane that forms subepithelial immune complexes, and diffuse thickening of the basement membrane [1]. According to the etiology, MN can be classified into idiopathic membranous nephropathy (IMN) and secondary membranous nephropathy (SMN). IMN is predominantly mediated by antibodies targeting the M-type phospholipase A2 receptor (anti-PLA2R) (85%) and thrombospondin type 1 domain containing 7A (anti-THSD7A) (3%–5%) [2]. Most of the autoantibodies from IMN patients belong to the IgG4 subclass [3]. IMN is characterized by heterogeneity of clinical outcomes, with 40% of patients undergoing spontaneous remission, and approximately 20% of patients progressing to severe kidney failure requiring renal replacement therapy [4]. The pathogenesis of MN is not yet clear, and may be related to environmental pollution [5], heredity, or autoimmunity.

The gut microbiome, residing in the intestine, is recognized as an important factor that contributes to both health and disease [6]. The composition of the intestinal flora is related to various diseases, such as inflammatory bowel disease [7,8], obesity [9], diabetes [10, 11], rheumatoid arthritis [12], atopic manifestations [13], liver sclerosis [14], cardiovascular disease [15], mental disease [16, 17], colorectal cancer [18, 19], and chronic kidney disease. The nonpodocyte circulating antigen, cationic bovine serum albumin, is believed to be the main antigen in early-childhood MN and may be affected by intestinal barrier formation and the use of infant formula [20]. Hence, we hypothesized that gut dysbiosis may play an important role in the development or exacerbation of IMN.

On the basis of next-generation sequencing technology, several human studies investigating the gut microbiota and IMN have been published recently [21,22,23,24]. Although all these studies indicated a compositional change in the fecal microbiota of IMN patients, their conclusions were not consistent. One possible explanation is that the participants were from different geographical locations and had different dietary habits, which may impact on gut microbial communities.

We thus conducted this study to investigate the fecal microbiota of IMN patients. We analyzed the bacterial community composition and diversity among patients and healthy individuals (cohabitants of the patients) using 16S rRNA gene sequencing. The role of important clinical parameters was also evaluated.


Study participants

A total of 30 hospitalized patients with IMN were enrolled in this study from May 2019 to December 2020 at the Second Affiliated Hospital of Harbin Medical University, along with 30 co-residents defined as healthy controls (HC) for whom urinary protein abnormalities had been ruled out. All participants, including patients and HC, were natives of Heilongjiang Province and were of Han ethnic group. The experimental group and the control group had similar dietary habits and living conditions.

The enrollment criteria for patients and their families were that they had not been administered antibiotics, probiotics, hormones, or immunosuppressants, and had not undergone any significant changes to diet or medication for at least 3 months. Exclusion criteria included diabetes, systemic lupus erythematosus, hepatitis B, hepatitis C, HIV, syphilis, and other diseases that can lead to SMN, as well as malignant tumors, gastrointestinal diseases, severe uncontrolled hypertension (diastolic blood pressure ≥ 120 mmHg and/or systolic blood pressure ≥ 220 mmHg), acute myocardial infarction or stroke within the last 6 months, and suspected/confirmed renal vascular disease. The demographic data and clinical features including sex, age, body mass index (BMI), estimated glomerular filtration rate (eGFR), total cholesterol (TC), triglyceride (TG), total protein (TP), alanine aminotransfease, aspartate aminotransferase, blood urea nitrogen, lactic dehydrogenase, high-density lipoprotein cholesterol (HDL-C), low-density lipoprotein (LDL-C), serum albumin were recorded.

This research received approval from the Ethics Committee from the Second Affiliated Hospital of Harbin Medical University in Harbin, Heilongjiang Province, China. All subjects provided written informed consent.

Stool sample collection

Fresh fecal samples were collected from the participants in the morning. After collection, the samples were immediately frozen and stored at − 80 ℃ prior to analyses.

DNA extraction and PCR amplification

Microbial community genomic DNA was extracted from fecal samples using the E.Z.N.A.® soil DNA Kit (Omega Bio-tek, Norcross, GA, USA) in accordance with the manufacturer’s instructions. The DNA extract was checked on a 1% agarose gel, and the DNA concentration and purity were determined using the NanoDrop 2000 UV–vis spectrophotometer (Thermo Fisher Scientific, Wilmington, MA, USA). The hypervariable region V3-V4 of the bacterial 16S rRNA gene was amplified with the primer pair 338F (ACTCCTACGGGAGGCAGCAG) and 806R (GGACTACHVGGGTATCTAAT) on an ABIGeneAmp®9700 PCR thermocycler (ABI, CA, USA). PCR amplification of the 16S rRNA gene was performed as follows: initial denaturation at 95 ℃ for 3 min, followed by 27 cycles of denaturing at 95 ℃ for 30 s, annealing at 55 ℃ for 30 s, and extension at 72 ℃ for 45 s, and a final extension step at 72 ℃ for 10 min and incubation at 10 ℃. The PCR mixture contained 4 µL of 5 × TransStartFastPfu buffer, 2 µL of 2.5 mM dNTPs, 0.8 µL of forward primer (5 µM), 0.8 µL of reverse primer (5 µM), 0.4 µL of TransStartFastPfu DNA polymerase, 10 ng of template DNA, and ddH2O up to a final volume of 20 µL. PCR reactions were performed in triplicate. The PCR product was then extracted from a 2% agarose gel and purified using the AxyPrep DNA Gel Extraction Kit (Axygen Biosciences, Union City, CA, USA) in accordance with the manufacturer’s instructions and quantified using the Quantus™ Fluorometer (Promega, USA).

Illumina MiSeq sequencing

Purified amplicons were pooled in equimolar amounts and subjected to paired-end sequencing on an Illumina MiSeq PE300 platform/NovaSeqPE250 platform (Illumina, San Diego, CA, USA) following the standard protocols of Majorbio Bio-Pharm Technology Co. Ltd. (Shanghai, China).

Processing of sequencing data

The raw 16S rRNA gene sequencing reads were demultiplexed, quality-filtered by fastp version 0.20.0, and merged by FLASH version 1.2.7 (, in accordance with the following criteria: (i) the 300 bp reads were truncated at any site with an average quality score of < 20 over a 50 bp sliding window, and truncated reads shorter than 50 bp and reads containing ambiguous characters were discarded; (ii) only overlapping sequences longer than 10 bp were assembled according to their overlapping sequence. The maximum mismatch ratio of overlapping regions was 0.2. Reads that could not be assembled were discarded; and (iii) samples were distinguished according to their barcode and primers. The sequence direction was adjusted, barcodes had to be an exact match, and two nucleotide mismatches were permitted in primer matching.

Operational taxonomic units (OTUs) with a 97% similarity cutoff were clustered using UPARSE version 7.1, and chimeric sequences were identified and removed. The taxonomy of each OTU representative sequence was analyzed by RDP Classifier version 2.2 against the 16S rRNA database using a confidence threshold of 0.7. The current species annotation results are from the silva138 version(

Statistical analysis

Continuous data with or without a normal distribution were expressed as the mean ± SD or the median. Normally distributed data were assessed using Levene’s test to judge the homogeneity of the variance, and t-tests were performed for homogeneous variance. If the variance was non-uniform, Mann–Whitney U tests were used and counts were analyzed using chi-squared tests.

A rarefaction curve was plotted to evaluate the sufficiency of the sample size and to estimate bacterial richness. α-Diversity parameters, Sobs and Chao indices were used to estimate species richness while Shannon and Simpson’s indices were used to estimate species diversity, were analyzed by software mothur (version v.1.30.1). β-Diversity provides a comparison of the taxonomic profiles between pairs of individual samples. β-Diversity was calculated based on the Bray–Curtis distance matrices and displayed using principal coordinates analysis (PCoA) by the R software (version 3.3.1), and the analysis of differences between groups was tested by analysis of similarities (ANOSIM). Both α-diversity and β-diversity were calculated at the OTU level. Bacterial taxonomic comparisons at the phylum and genus levels were performed between two groups using the Wilcoxon rank sum test. Linear discriminant analysis (LDA) effect size (LEfSe) was used to identify the characteristic microbiota and explain the differences between patients and HC. Only taxa with an LDA score > 2 and a significance of α < 0.05 were shown. The results were plotted in a cladogram based on their phylogenetic relationship. The correlation between biochemical indicators and various microbes were calculated by Spearman’s rank correlation coefficient and visualized by heatmap in R using the pheatmap package. Relative predictive performance analysis was performed by random forest analysis in R using the randomForest package and receiver operating characteristic (ROC) curve analysis. The gut microbiota was explored based on Kyoto Encyclopedia of Genes and Genomes (KEGG, using Phylogenetic Investigation of Communities by Reconstruction of Unobserved States (PICRUSt). Between-group differences in functional pathways indicated by taxa variations were assessed using the Mann–Whitney U test. P values < 0.05 were considered statistically significant. All statistical analyses were performed using the computer software Statistical package for the Social Sciences, version 25.0 (IBM Corp., Armonk, NY, USA).


General characteristics of all participants

The study included 30 IMN patients, confirmed by renal biopsy, with a mean age of 51.3 ± 10.25 years. All IMN patients had an estimated glomerular filtration rate(eGFR) ≥ 55 mL/min/1.73 m2. The baseline characteristics of the IMN and control groups are summarized. (Table 1, Supplemental Table 1).

Table 1 Baseline characteristics of participants

Intestinal flora diversity

A total of 2,075,757,404 reads were obtained for the 60 subjects by V4V5 16S rRNA pyrosequencing. After quality control and paired-end read merging, we obtained 1,418,592,102 high-quality reads, accounting for 68.34% of the total reads. A total of 3,448,102 sequences for all samples were used in downstream bioinformatics analysis, and the average sequence length of the merged sequences was 412 bp.

A total of 916 OTUs were obtained at a 97% homology level cutoff. The number of OTUs was larger for the HC group. As exhibited in the Venn diagram (Fig. 1A), the number of common OTUs between the HC and IMN groups was 750; the IMN group had 88 specific OTUs and the HC group had 78 specific OTUs not shared by the other group.

Fig. 1
figure 1

Gut microbiota compositions differed between IMN and HC groups. A Venn diagram of IMN and HC; B rarefaction analysis of both groups.The curves in all samples were near saturation, which means the sequencing depth was sufficient to capture most gene diversity

To determine whether the sequencing adequately captured the diversity of the gut microbiota, rarefaction analysis was performed. The curves in all samples were near saturation, which meant that the sequencing depth was sufficient to capture most of the gene diversity (Fig. 1B). Chao and Sobs indices on the OTU profile were used to estimate community richness. Shannon and Simpson indices were used to estimate community diversity. α-Diversity analysis revealed significant differences among groups; i.e., there was a significantly lower diversity in IMN patients compared with the HC group (Chao, P = 0.048; Sobs, P = 0.024; Shannon, P = 0.016; Simpson, P = 0.032) (Table 2, Fig. 2, Supplemental Table 2).

Table 2 α-Diversity between HC and IMN
Fig. 2
figure 2

Composition of gut microbiota is significantly altered in primary hypothyroidism patients. α-Diversity indices (Chao, Shannon, Simpson and Sobs index) of intestinal flora in healthy individuals and Idiopathic membranous nephropathy patients. The Wilcoxon rank sum test was used to detect significant changes. *P < 0.05; **P < 0.01. HC: healthy individuals, n = 30; IMN: Idiopathic membranous nephropathy, n = 30

Principal co-ordinates analysis (PCoA) searches for main coordinates based on a distance matrix; i.e., shows an higher similarity when OTUs decreased. PCoA based on Bray–Curtis dissimilarity at the OTU level revealed that the microbiota composition of IMN and HC patients was significantly different (P = 0.002, r = 0.066) (Fig. 3).

Fig. 3
figure 3

A principal component (PCoA) score plot based on Bray–Curtis distance matrix for all participants. Each point represents the composition of the intestinal microbiota of one participant. The ellipses do not represent statistical significance but rather serve as a visual guide to illustrate group differences. The ANOSIM test was used

Bacterial taxa differences between IMN patients and healthy controls

Both IMN patients and HC showed typical microbiome structures. Most bacteria identified belonged to the phyla Bacteroidota, Firmicutes, Proteobacteria, and Actinobacteriota. IMN patients showed an increased abundance of Proteobacteria and Actinobacteriota and a decreased abundance of Bacteroidota and Firmicutes (Fig. 4A).

Fig. 4
figure 4

Relative abundance and taxonomic differences of fecal microbiota in the IMN patients and healthy controls. A Microbiome composition in the two groups at the phylum level. B Microbiome composition in the two groups at the genus level. The figure shows the top 20 species in each group based on their relative abundances. C, D LDA and LEfSe of IMN and HC

Within the Proteobacteria, the abundance of Enterobacteriaceae was significantly higher in IMN patients than in the HC at the family level (P < 0.05). For the Bacteroidota family, a significant decrease was found in the abundance of Bacteroidaceae in the IMN group compared with the HC group (P < 0.01) (Table 3).

Table 3 Relative abundance of fecal microbiota in IMN patients and healthy controls

At the genus level, Bacteroides, Faecalibacterium, Escherichia-Shigella, Agathobacter, and Citrobacter were the most abundant genera in both groups. Some of these genera showed a significant difference between these two groups (Fig. 4B, Table 3).

To determine different taxa from the phylum to genus level among the two groups, the LEfSe algorithm was used. Bacterial species that differed significantly (P < 0.05) between the IMN patients and HC were screened using LEfSe to construct a cladogram (Fig. 4C, D).

The histogram obtained by LDA showed the increased abundances of 19 genera in the IMN group. Among the HC group, the abundances of 29 genera were significantly higher than for the IMN patients. (Supplemental Table 6).

Association between the fecal microbiota and IMN clinical characteristics

Correlations between the differences in gut microbiota and clinical biomarkers were further evaluated by Spearman’s rank correlation analysis. At the phylum level, Actinobacterioia showed a strong negative correlation with the eGFR (r =  − 0.414, P = 0.023), while Campylobacterota showed a positive correlation with LDL-C (r = 0.375, P = 0.041) and Bacteroidota showed a positive correlation with TP (r = 0.373, P = 0.043) (Fig. 5A, Supplemental Table 3).

Fig. 5
figure 5

Heatmaps showing correlations between fecal microbiota and IgAN clinical parameters. A Correlations between microbiota phyla and IMN clinical characteristics. B Correlations between microbiota genera and IMN clinical characteristics. The correlations were performed by heatmap in R version 3.3.1 ( using the pheatmap package. The figure shows the top 20 species in each group based The intensity of the color indicates the r value (correlation). The red color represent positive score and the blue color represent negative score. *P < 0.05 and **P < 0.01

At the genus level, Bacteroides, Bifidobacterium, and Veillonella showed strong negative correlations with TG (r =  − 0.420, P = 0.021; r =  − 0.447, P = 0.013; − 0.555, P = 0.001). Dorea also showed a negative correlation with HDL-C (r =  − 0.390, P = 0.033). Faecalibacterium and Roseburia exhibited a negative correlation with BMI (r =  − 0.409, P = 0.024; r =  − 0.390, P = 0.033). A negative correlation also existed between Blautia and age (r =  − 0.397, P = 0.030) (Fig. 5B, Supplemental Table 4).

Microbial markers for the potential detection of IMN

Four bacteria, determined by random forest analysis, were found to be characteristic of IMN (Fig. 6A). To illustrate the microbial signatures of IMN patients and HC and to further build a predictive model according to the fecal microbiota profiles of the significantly different taxa abundances at the genus level, ROC curves for classifying IMN from HC were developed. We could detect IMN individuals accurately based on the combination of the four genera (g__Moryella, g__Lachnospira, g__Gemella, g__Lachnospiraceae_NC200 4_group), as indicated by the area under the receiver operating curve (AUC) of up to 0.96 (Fig. 6B).

Fig. 6
figure 6

The predictive model based on the relative abundances at genus level by using RF model. A The importance of each genus in the predictive model was evaluated by the mean decreasing accuracy and the Gini coefficient. B ROC curve analysis generated by the RF using 4 genera in the fecal microbiota

Functional pathway predictions

From metagenome predictions based on PICRUSt analysis, abundant bacterial functions related to lipid metabolism were observed among IMN cases and the controls. At the level of the KEGG pathway, microbial gene functions included pathways involved in bacterial invasion of epithelial cells (ko05100), Alpha-linolenic acid metabolism (ko00592), Staphylococcus aureus infection (ko05150), and Arachidonic acid metabolism (ko00590), which were higher in the fecal microbiome of the IMN group compared with the control group (P < 0.05; Fig. 7, Supplemental Table 5). This indicated the differential function of the fecal microbiota between IMN patients and healthy individuals.

Fig. 7
figure 7

Predicted microbial functions using PICRUSt. Significant KEGG pathways ( at level 3 between the IMN patients and healthy controls. P < 0.05 is statistically significant


The intestinal microecosystem is continuously evolving, and its microbial diversity and abundance play important roles in maintaining a normal physiological state; however, the host also influences the intestinal microflora [25]. In this study, 16S rRNA sequencing was performed on stool samples from IMN patients and healthy co-residents. The relative abundances of each intestinal microbial species and their diversities and compositions were compared between IMN patients and HC. The results showed significantly decreased intestinal flora diversity and an altered intestinal flora composition in the IMN patients compared with those of the HC.

External factors (such as changes in diet or geographical region) and internal factors (such as physiological changes in response to a disease process) can both induce shifts in microbiota composition and function [26, 27]. We recruited the patients co-residents as the control group to reduce the influence of diet and living conditions on the results. Our findings demonstrated that the community richness of the fecal microbiota in the IMN patients was significantly lower than that in the HC, which was consistent with previous observations [21,22,23,24]. However, there were differences in the taxonomic profile of the fecal microbiota compared with previous research.

Bacteroidota, Firmicutes, Proteobacteria, and Actinobacteriota were the predominant phyla in both groups, which was consistent with reports from previous Chinese studies [21]. When compared with the HC, Firmicutes and Bacteroidota were depleted in patients with IMN. The phylum Bacteroidota produces acetate and propionate, whereas Firmicutes produces butyrate as its primary metabolic end product [28]. Proteobacteria not only possesses genes that encode the urease enzyme [29] but also promotes inflammatory uremic toxins [30].

Proteobacteria is regulated by a gamma-Proteobacteria-specific IgA response, thus a higher abundance of gamma-Proteobacteria promotes the production of inflammation [31]. Proteobacteria and Fusobacteria have also been found to be enriched in type 2 diabetic patients with chronic kidney disease [32].

The main genera that differed between the IMN and HC groups were Agathobacter, Faecalibacterium, Bacteroides, and Subdoligranulum. However, the physiological and pathological effects of these genera have not been investigated in detail. Bacteroides is a dominant intestinal organism that plays a pivotal role in improving nutrient utilization, accelerating intestinal mucosal angiogenesis, developing the immune system, improving immunity, and maintaining the balance of intestinal flora. Bacteroides can produce short-chain fatty acids in the colon, and significantly enhance the growth of colonic lamina propria Treg cells after colonization. Such enhancement further exerts immunomodulatory effects [22]. Faecalibacterium is abundant in the fecal microbiome of healthy humans. It is a member of Clostridium IV (Firmicutes) and produces short-chain fatty acids (including butyrate) [33].

Collectively, IMN correlates with changes in abundance of some bacterial taxa. Spearman correlation analysis confirmed that opportunistic pathogens, such as Actinobacteriota, negatively correlated with the eGFR, which is classical marker of renal damage. Patients with IMN often exhibit hyperlipidemia. Dorea, which has been shown to be closely related to obesity [34], was negatively correlated with HDL-C levels and was significantly increased in IMN patients. Bacteroides, Bifidobacterium, and Veillonella, associated with adiposity and lipid levels [35], were negatively correlated with TG in our research. PICRUSt analysis indicated that the gut microbiota that were altered in abundance between IMN patients and HC were related to lipid metabolism. In 1982, Moorhead [36] proposed the hypothesis that lipids exert renal toxicity. In recent years, research has shown that abnormal blood lipid levels can cause renal glomerular membrane hyperplasia, epidural matrix aggregation, renal glomerular epithelial cell damage, and inflammatory cell immersion, all of which can damage the kidney directly [37]. Hyperlipidemia, defined as one of the characteristics of IMN, increases blood viscosity and the probability of venous thrombosis, and accelerates the progression of renal damage.

The results of our study demonstrated that IMN could be distinguished relatively accurately from HC via the detection of four specific genera in feces. However, further studies are required to establish the exact role of the gut microbiota in IMN progression and to identify diagnostic and predictive markers.

Current treatments, including supportive therapy, corticosteroids and immunosuppressive agents, are not effective in all patients. In one reported case, fecal microbiota transplantation as used to treat a patient with MN and chronic diarrhea, whose symptoms ameliorated and renal function improved [38]. Research indicates that the gut microbiota may contribute to the pathogenesis of this disease.

Our study had certain limitations. First, 16S rRNA sequencing was performed among participants from the same area of residence (Han nationality). Second, our sample size was relatively small; therefore, a large-scale study involving different populations and strict dietary controls are needed to confirm our results. We minimized the influence of other factors, such as weight, BMI, sex, and diet in our study; however, the proportion of certain ingredients in the diet (such as salt or fatty acids) was not controlled. Third, to identify further novel microbial and metabolic biomarkers and to further characterize the biological mechanisms and pathways involved, bacterial metagenomic sequencing and metabonomic analysis are needed in future studies. Finally, a cross-sectional study cannot confirm a causal relationship because of variations in the gut microbiota over time.

Overall, the up- or down-regulation of these strains may be closely related to the occurrence of IMN. The molecular mechanisms involved in the relationship between IMN and intestinal bacterial changes will be studied in future analyses.


By comparing the intestinal flora of IMN and HC, we found that microbial diversity was significantly different in IMN patients. The variations identified in the gut microbiome will provide a preliminary reference for future research on the interaction mechanism between the intestinal flora and IMN.

Availability of data and materials

All data generated or analyzed during this study are included in this published article.



Analysis of similarities


Area under curve


Body mass index


Estimated glomerular filtration rate


End-stage renal disease


Fecal microbiota transplantation


High density liptein cholesterol


Immunoglobulin A nephropathy


Idiopathic membranous nephropathy


Kyoto Encyclopedia of Genes and Genomes


Low-Density Lipoprotein Cholesterol


Linear discriminant analysis effect size


Operational taxonomic units


Principal coordinates analysis


Phylogenetic Investigation of Communities by Reconstruction of Unobserved States


Random forest


Secondary membranous nephropathy


Serum total cholesterol




Total protein


  1. Couser WG. The pathogenesis of human membranous nephropathy: we are (almost) there. Kidney Int. 2020;97(5):849–52.

    Article  CAS  Google Scholar 

  2. Beck L, Bonegio R, Lambeau G, Beck D, Powell D, Cummins T, Klein J, Salant D. M-type phospholipase A2 receptor as target antigen in idiopathic membranous nephropathy. N Engl J Med. 2009;361(1):11–21.

    Article  CAS  Google Scholar 

  3. Vriese AD, Glassock RJ, Nath KA, Sethi S, Fervenza FC. A proposal for a serology-based approach to membranous nephropathy. J Am Soc Nephrol. 2017;28(2):421–30.

  4. Ronco P, Debiec H. Pathophysiological advances in membranous nephropathy: time for a shift in patient’s care. Lancet. 2015;385(9981):1983–92.

    Article  Google Scholar 

  5. Xu X, Wang G, Chen N, Lu T, Nie S, Xu G, Zhang P, Luo Y, Wang Y, Wang X. Long-term exposure to air pollution and increased risk of membranous nephropathy in China. J Am Soc Nephrol. 2016;27:3739–46.

    Article  CAS  Google Scholar 

  6. Meijers B, Evenepoel P, Anders HJ. Intestinal microbiome and fitness in kidney disease. Nat Rev Nephrol. 2019;15(9):531–45.

    Article  Google Scholar 

  7. Frank DN, St. Amand AL, Feldman RA. Molecular-phylogenetic characterization of microbial community imbalances in human inflammatory bowel diseases. Proc Natl Acad Sci. 2007;104(34):13780–5.

    Article  CAS  Google Scholar 

  8. Ott S, Musfeldt M, Wenderoth D, Hampe J, Brant O, Fölsch U, Timmis K, Schreiber S. Reduction in diversity of the colonic mucosa associated bacterial microflora in patients with active inflammatory bowel disease. Gut. 2004;53(5):685–93.

    Article  CAS  Google Scholar 

  9. Turnbaugh PJ, Hamady M, Yatsunenko T, Cantarel BL, Duncan A, Ley RE, Sogin ML, Jones WJ, Roe BA, Affourtit JP. A core gut microbiome in obese and lean twins. Nature. 2009;457(7228):480.

    Article  CAS  Google Scholar 

  10. Karlsson FH, Tremaroli V, Nookaew I, Bergstroem G, Behre CJ, Fagerberg B, Nielsen J, Baeckhed F. Gut metagenome in European women with normal, impaired and diabetic glucose control. Nature. 2013;498(7452):99–103.

    Article  CAS  Google Scholar 

  11. Brown CT, Davis-Richardson AG, Giongo A, Gano KA, Crabb DB, Mukherjee N, Casella G, Drew JC, Ilonen J, Knip M. Gut microbiome metagenomics analysis suggests a functional model for the development of autoimmunity for type 1 diabetes. PLoS One. 2011;6(10):e25792.

    Article  CAS  Google Scholar 

  12. Vaahtovuo J, Munukka E, Korkeamaki M, Luukkainen R, Toivanen P, KorkeamäKi M, Luukkainen R, Toivanen P. Fecal microbiota in early rheumatoid arthritis. J Rheumatol. 2008;35(8):1500–5.

    CAS  PubMed  Google Scholar 

  13. Penders J, Thijs C, Van D, Kummeling I, Snijders B, Stelma F, Adams H, Ree RV, Stobberingh EE. Gut microbiota composition and development of atopic manifestations in infancy: the KOALA Birth Cohort Study. Gut. 2007;56(5):661–7.

    Article  CAS  Google Scholar 

  14. Chen Y, Yang F, Lu H, Wang B, Li L. Characterization of fecal microbial communities in patients with liver cirrhosis. Hepatology. 2011;54(2):562–72.

    Article  Google Scholar 

  15. Wang Z, Klipfell E, Bennett BJ, Koeth R, Levison BS, Dugar B, Feldstein AE, Britt EB, Fu X, Chung YM. Gut flora metabolism of phosphatidylcholine promotes cardiovascular disease. Nature. 2011;472(7341):57–63.

    Article  CAS  Google Scholar 

  16. Song Y, Liu C, Finegold SM. Real-time PCR quantitation of clostridia in feces of autistic children. Appl Environ Microbiol. 2004;70(11):6459–65.

    Article  CAS  Google Scholar 

  17. Parracho HM. Differences between the gut microflora of children with autistic spectrum disorders and that of healthy children. J Med Microbiol. 2005;54(10):987–91.

    Article  Google Scholar 

  18. Scanlan P, Shanahan F, Clune Y, Collins J, O’Sullivan G, O’Riordan M, Holmes E, Wang Y, Marchesi J. Culture-independent analysis of the gut microbiota in colorectal cancer and polyposis. Environ Microbiol. 2008;10(3):789–98.

    Article  CAS  Google Scholar 

  19. Sobhani I, Tap J, Roudot-Thoraval F, Roperch J, Letulle S, Langella P, Corthier G, Tran Van Nhieu J, Furet J. Microbial dysbiosis in colorectal cancer (CRC) patients. PloS One. 2011;6(1):e16393.

    Article  CAS  Google Scholar 

  20. Debiec H, Lefeu F, Kemper M, Niaudet P, Deschênes G, Remuzzi G, Ulinski T, Ronco P. Early-childhood membranous nephropathy due to cationic bovine serum albumin. N Engl J Med. 2011;364(22):2101–10.

    Article  CAS  Google Scholar 

  21. Zhang J, Luo D, Lin Z, Zhou W, Rao J, Li Y, Wu J, Peng H, Lou T. Dysbiosis of gut microbiota in adult idiopathic membranous nephropathy with nephrotic syndrome. Microb Pathog. 2020;147:104359.

    Article  CAS  Google Scholar 

  22. Lang R, Wang X, Li A, Liang Y, Zhu B, Shi B, Zheng Y, Yu R. Effects of Jian Pi Qu Shi Formula on intestinal bacterial flora in patients with idiopathic membranous nephropathy: a prospective randomized controlled trial. Chronic Dis Transl Med. 2020;6(2):124–33.

    PubMed  PubMed Central  Google Scholar 

  23. Dong R, Bai M, Zhao J, Wang D, Ning X, Sun S. A comparative study of the gut microbiota associated with immunoglobulin a nephropathy and membranous nephropathy. Front Cell Infect Microbiol. 2020;10:557368.

    Article  CAS  Google Scholar 

  24. He H, Lin M, You L, Chen T, Liang Z, Li D, Xie C, Xiao G, Ye P, Kong Y, et al. Gut microbiota profile in adult patients with idiopathic nephrotic syndrome. Biomed Res Int. 2021;2021:8854969.

    PubMed  PubMed Central  Google Scholar 

  25. Huang X, Fan X, Ying J, Chen S. Emerging trends and research foci in gastrointestinal microbiome. J Transl Med. 2019;17(1):67.

    Article  Google Scholar 

  26. Dethlefsen L, McFall-Ngai M, Relman DA. An ecological and evolutionary perspective on human-microbe mutualism and disease. Nature. 2007;449(7164):811–8.

    Article  CAS  Google Scholar 

  27. Yatsunenko T, Rey F, Manary M, Trehan I, Dominguez-Bello M, Contreras M, Magris M, Hidalgo G, Baldassano R, Anokhin A, et al. Human gut microbiome viewed across age and geography. Nature. 2012;486(7402):222–7.

    Article  CAS  Google Scholar 

  28. Macfarlane S, Macfarlane G. Regulation of short-chain fatty acid production. Proc Nutr Soc. 2003;62(1):67–72.

    Article  CAS  Google Scholar 

  29. Wong J, Piceno Y, DeSantis T, Pahl M, Andersen G, Vaziri N. Expansion of urease- and uricase-containing, indole- and p-cresol-forming and contraction of short-chain fatty acid-producing intestinal microbiota in ESRD. Am J Nephrol. 2014;39(3):230–7.

    Article  CAS  Google Scholar 

  30. Onal E, Afsar B, Covic A, Vaziri N, Kanbay M. Gut microbiota and inflammation in chronic kidney disease and their roles in the development of cardiovascular disease. Hypertens Res. 2019;42(2):123–40.

    Article  Google Scholar 

  31. Mirpuri J, Raetz M, Sturge C, Wilhelm C, Benson A, Savani R, Hooper L, Yarovinsky F. Proteobacteria-specific IgA regulates maturation of the intestinal microbiota. Gut Microbes. 2014;5(1):28–39.

    Article  Google Scholar 

  32. Salguero M, Al-Obaide M, Singh R, Siepmann T, Vasylyeva T. Dysbiosis of Gram-negative gut microbiota and the associated serum lipopolysaccharide exacerbates inflammation in type 2 diabetic patients with chronic kidney disease. Exp Ther Med. 2019;18(5):3461–9.

    CAS  PubMed  PubMed Central  Google Scholar 

  33. de Paiva C, Jones D, Stern M, Bian F, Moore Q, Corbiere S, Streckfus C, Hutchinson D, Ajami N, Petrosino J, et al. Altered mucosal microbiome diversity and disease severity in Sjögren syndrome. Sci Rep. 2016;6:23561.

    Article  Google Scholar 

  34. Companys J, Gosalbes M, Pla-Pagà L, Calderón-Pérez L, Llauradó E, Pedret A, Valls R, Jiménez-Hernández N, Sandoval-Ramirez B, Del Bas J, et al. Gut microbiota profile and its association with clinical variables and dietary intake in overweight/obese and lean subjects: a cross-sectional study. Nutrients. 2021;13(6):2032.

    Article  CAS  Google Scholar 

  35. Moayyeri A, Hammond C, Valdes A, Spector T. Cohort Profile: TwinsUK and healthy ageing twin study. Int J Epidemiol. 2013;42(1):76–85.

    Article  Google Scholar 

  36. Moorhead J, Wheeler D, Varghese Z. Glomerular structures and lipids in progressive renal disease. Am J Med. 1989;87:12N-20N.

    CAS  PubMed  Google Scholar 

  37. Ma H, Sandor D, Beck L. The role of complement in membranous nephropathy. Semin Nephrol. 2013;33(6):531–42.

    Article  CAS  Google Scholar 

  38. Zhou G, Zeng J, Peng L, Wang L, Zheng W, Wu Di, Yang Y. Fecal microbiota transplantation for membranous nephropathy. CEN Case Rep. 2021;10(2):261–4.

    Article  Google Scholar 

Download references


We thank Liwen Bianji (Edanz) ( for editing the language of a draft of this manuscript.



Author information

Authors and Affiliations



Research idea and study design: MFL, LXW, XYD. Data acquisition:MFL, LXW, JS, QSZ, HY, YZ, LX, CZ, RZ, XYD. Data analysis/interpretation: MFL, LXW, XYD. Manuscript drafting: MFL, XYD. Editing and revising: MFL, XYD. All authors have approved the submitted version and agreed to be accountable for the author’s own contributions.

Corresponding author

Correspondence to Xuanyi Du.

Ethics declarations

Ethics approval and consent to participate

The study was conducted according to the Declaration of Helsinki. It was approved by the Medical Ethics Committee of he Second Affiliated Hospital of Harbin Medical University (No.KY2016-221). All participants provided written informed consent in this study.

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: Supplemental Table 1.

Baseline characteristics of participants.

Additional file 2: Supplemental Table 2.

α-Diversity was assessed with 4 different diversity analysis methods.

Additional file 3: Supplemental Table 3. 

Gut microbiota and clinical biomarkers evaluated by Spearson correlation analysis at the phylum level.

Additional file 4: Supplemental Table 4. 

Gut microbiota and clinical biomarkers evaluated by Spearson correlation analysis at the genus level.

Additional file 5: Supplemental Table 5. 

KEGG level 3 pathways of abundance changed intestinal flora.

Additional file 6: Supplemental Table 6.

Bacterial species that differed significantly between the IMN patients and HC by LEfSe.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Li, M., Wei, L., Sun, J. et al. Association of gut microbiota with idiopathic membranous nephropathy. BMC Nephrol 23, 164 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Idiopathic membranous nephropathy
  • Gut microbiota
  • 16S rRNA sequencing