Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Immunol., 14 January 2026

Sec. Mucosal Immunity

Volume 16 - 2025 | https://doi.org/10.3389/fimmu.2025.1665585

This article is part of the Research TopicMucosal Immune Dysregulation in IgA NephropathyView all articles

Altered acetate metabolism and signaling in IgA nephropathy: an integrated gut microbiome and glomerular spatial transcriptome analysis

Jung Hun Koh&#x;Jung Hun Koh1†Sehoon Park,&#x;Sehoon Park1,2†Minji KangMinji Kang3Ji In ParkJi In Park4Jangwook LeeJangwook Lee5Hyunjeong ChoHyunjeong Cho6Ji Eun KimJi Eun Kim7Hoonsik NamHoonsik Nam8Doyeon KimDoyeon Kim8Minshu LiMinshu Li8Sunghyouk ParkSunghyouk Park8Kyung Chul MoonKyung Chul Moon9Hyun Je KimHyun Je Kim3Yon Su Kim,Yon Su Kim1,2Dong Ki Kim,Dong Ki Kim1,2Hajeong Lee,*Hajeong Lee1,2*
  • 1Department of Internal Medicine, Seoul National University Hospital, Seoul, Republic of Korea
  • 2Department of Internal Medicine, Seoul National University College of Medicine, Seoul, Republic of Korea
  • 3Department of Biomedical Sciences, Seoul National University Graduate School, Seoul, Republic of Korea
  • 4Department of Internal Medicine, Kangwon National University Hospital, Chuncheon, Republic of Korea
  • 5Department of Internal Medicine, Dongguk University Ilsan Hospital, Ilsan, Republic of Korea
  • 6Department of Internal Medicine, Chungbuk National University Hospital, Cheongju, Republic of Korea
  • 7Department of Internal Medicine, Korea University Guro Hospital, Seoul, Republic of Korea
  • 8College of Pharmacy, Natural Products Research Institute, Seoul National University, Seoul, Republic of Korea
  • 9Department of Pathology, Seoul National University College of Medicine, Seoul, Republic of Korea

Introduction: IgA nephropathy (IgAN) is the most common primary glomerulonephritis, and emerging evidence implicates the gut microbiome in its pathogenesis. Additional studies focusing on the molecular mechanisms linking gut microbial signals to intraglomerular changes are warranted.

Methods: We performed 16S rRNA-based microbial profiling of fecal samples of 172 IgAN patients, 51 healthy controls, and other glomerular disease controls including 15 diabetic nephropathy, 35 minimal change disease, and 63 membranous nephropathy cases. Serum and fecal acetate levels were measured by liquid chromatography–mass spectrometry. Glomerular spatial transcriptomic profiling was performed with the GeoMx Digital Spatial Profiler. DESeq2 analysis was performed to identify differentially expressed genes, followed by gene ontology annotations.

Results: Beta diversity differed significantly between IgAN and healthy controls (p = 0.001). While no single taxon showed consistent differences in abundance, the methanogenesis from acetate pathway was significantly enriched in IgAN, accompanied by an increased proportion of major acetate-producing gut microbial genera. Serum acetate levels were elevated in IgAN (p = 0.03), while fecal acetate levels were comparable to those in healthy controls. In glomerular transcriptomes, functional annotations of 1,227 upregulated and 1,078 downregulated genes in IgAN indicated decreased activities of G protein-coupled receptors, short-chain fatty acid transporters, and beta-1,3-galactosyltransferases.

Discussion: IgAN is characterized by gut microbial enrichment in acetate metabolism and increased systemic acetate levels, along with altered intraglomerular expression of metabolic and signaling genes. These findings suggest a gut microbiome–glomerular signaling axis contributing to disease pathogenesis.

1 Introduction

Kidney diseases, spanning a broad spectrum of etiologies, are a growing health concern worldwide (1). Beyond traditional risk factors, gut microbial dysbiosis has been increasingly implicated in their pathophysiology (2, 3). Alterations in gut microbial composition and related metabolite profiles have been reported in chronic kidney disease (CKD) and in specific glomerular diseases (46), with associations to systemic consequences such as cardiorenal syndrome and mineral bone disorder (7, 8). The crosstalk between the gut microbiota and the kidney, dubbed the gut-kidney axis, involves multiple mediators including microbe-derived uremic toxins and other amino acid derivatives (911), whose causal associations to glomerular diseases were demonstrated by recent Mendelian randomization studies (11, 12). Together, a growing body of evidence indicates a key role of intestinal dysbiosis and metabolite imbalance across kidney diseases.

Immunoglobulin A (IgA) nephropathy (IgAN), the most prevalent primary glomerulonephritis worldwide, is closely connected to the gut microbiota in its immune-mediated pathophysiology. Central to its pathogenesis is the overproduction of galactose deficient-IgA1 (Gd-IgA1) and associated immune complexes (13, 14). Given the established role of IgA in mucosal immunity, host-microbiota interactions in mucosa-associated lymphoid tissues have been proposed as a key driver of Gd-IgA1 production (13). Supporting this concept, genome-wide association studies have also suggested genetic associations between IgAN and mucosal immune pathways including those related to microbial sensing and response (15, 16). These findings have motivated the introduction of tonsillectomy and enteric budesonide as therapeutic options (17, 18), and efforts are underway to utilize microbiota-modulating therapies including probiotics and fecal microbiota transplantation (19).

Nevertheless, the development of IgAN requires additional components, as described by the “multi-hit” model (14, 20). Mesangial cells actively contribute by binding IgA1 and driving proliferative signals (21). B cells play a pivotal role in producing Gd-IgA1 as well as autoantibodies against Gd-IgA1, although relatively little is known about their activation and potential interactions with the kidney glomeruli in IgAN (22). Unraveling the crosstalk among the gut microbiota, mucosal and systemic immune systems, and the kidney as a target organ is therefore critical for developing additional therapeutic strategies.

Data-driven, multi-omics approaches have the potential to yield novel insights into the complex pathophysiology of IgAN. Advances in high-throughput sequencing techniques have enabled comparative analyses of tonsillar or gut microbiome in patients with IgAN versus healthy or disease controls (21, 2325). However, reported taxonomic differences varied among studies, and their functional relevance remains largely unclear. Spatial transcriptomics represent another key advancement, allowing for substructure-specific profiling of gene expression in elaborately organized structures like the nephron (26). We previously explored human IgAN-specific glomerular transcriptomic changes to uncover proinflammatory signals associated with mesangial proliferation preceding overt morphologic changes (27).

In this study, we performed 16S rRNA-based gut microbiome profiling integrated with Nanostring GeoMx-based glomerular spatial transcriptomics to characterize the molecular signature of IgAN. Comparisons with healthy controls and other glomerular disease patients enabled identification of specific changes in both the gut microbiome and the glomerular transcriptome. Specifically, IgAN showed gut microbial alterations in short-chain fatty acid (SCFA) metabolism, particularly acetate, along with glomerular transcriptional changes in SCFA transporters and G protein-coupled receptors (GPCRs). Together, these findings suggest a previously underappreciated axis of gut microbial metabolites and impaired renal sensing via GPCR downregulation, offering a potential channel for the gut-kidney interaction in IgAN.

2 Materials and methods

2.1 Ethics approval

This study was conducted in accordance with the Declaration of Helsinki, and ethical approval was granted by the Institutional Review Board (IRB) of Seoul National University Hospital (SNUH) for gut microbiome analysis (IRB No. 2205-104-1325) and spatial transcriptomic profiling (IRB No. 2205-085-1324). Written informed consent was obtained from all patients included in the study prior to collection of any samples or clinical information.

2.2 Participant cohorts and biospecimen acquisition

For gut microbial analysis, stool samples were sourced from the KOrea Renal biobank NEtwoRk System TOward Next-generation analysis (KORNERSTONE) repository (28). Biopsy-proven cases of IgAN and other common glomerular diseases, namely diabetic nephropathy (DN), minimal change disease (MCD), and membranous nephropathy (MN), were included. In addition, stool samples from consenting live kidney donors at SNUH were included as healthy controls. The following inclusion criteria were applied: (1) age ≥18 years, (2) estimated glomerular filtration rate (eGFR) ≥60 mL/min/1.73 m2 to minimize the potential confounding from uremic metabolites we previously observed (29), (3) no exposure to antibiotics or immunosuppressant treatment within one month of sample collection, (4) no history of bowel resection or inflammatory bowel disease.

For spatial transcriptomics analysis, IgAN samples were chosen from archived formalin-fixed paraffin-embedded (FFPE) slides, with biopsies performed between 2018 and 2022. Time-zero allograft biopsies after living donor kidney transplantation served as healthy controls, while DN, MCD, and MN samples archived between 2009 and 2021 were selected as disease controls. The selected participants had age 18–70 years, eGFR ≥30 mL/min/1.73 m2, and no immunosuppressant admission prior to biopsy to minimize the confounding from changes related to advanced aging, late-stage chronic kidney disease, or treatment decisions. Cases with fewer than 10 sampled glomeruli and non-DN cases with DN involvement on pathology were excluded.

2.3 Gut microbiome profiling

Stool DNA extraction and sequencing protocols were performed as described previously (29). Stool samples were collected on the day of kidney biopsy, immediately stored at -20 °C, then moved to -80 °C storage within 24 hours. DNA was extracted with the QIAamp Fast DNA Stool Mini Kit (Qiagen, Hilden, Germany). The V3-V4 hypervariable region of the 16S rRNA gene was amplified, followed by library preparation and sequencing through the Illumina MiSeq system according to the manufacturer’s protocol (Illumina, CA, USA).

For data analysis, raw sequence reads in FASTQ format were imported into the Quantitative Insights Into Microbial Ecology (QIIME2) microbiome analysis platform. The sequences were processed through DADA2 for filtering, trimming, and correction for low-quality reads and chimeric sequences (30). The denoised amplicon sequence variants were then annotated up to the genus level using the SILVA 138 database as reference (31). For diversity analysis, the samples were rarefied to a sampling depth of 10,000, and Shannon diversity indices and Bray-Curtis distances were computed through the QIIME2 q2-diversity plugin. Prediction of metagenomic functions was performed through Phylogenetic Investigation of Communities by Reconstruction of Unobserved States 2 (PICRUSt2) (32), and the predicted functions for each sample were expressed in terms of MetaCyc pathway abundances (33). Microbial genera and functional pathway abundances were compared between IgAN and each control group using ANOVA-like differential expression 2 (ALDEx2) (34), where false-discovery rates (FDR) < 0.05 were considered statistically significant.

2.4 Spatial transcriptomic profiling of glomeruli

Slide preparation and processing procedures for spatial transcriptomics profiling were described in detail previously (27). Briefly, 5 µm FFPE kidney biopsy sections underwent deparaffinization and epitope retrieval. A total of 18,677 target genes were labeled via in situ hybridization with the GeoMx Whole Transcriptome Atlas, which contains ultraviolet-photocleavable oligonucleotide identifiers. Three representative glomeruli per sample were selected as regions of interest. The oligonucleotide identifiers were collected, amplified through polymerase chain reactions, and sequenced on an Illumina NovaSeq 6000.

Data preprocessing and quality control were performed with the GeoMx Digital Spatial Profiler Data Analysis Suite (v2.4), removing low-performing genes, defined as those expressed in fewer than 50% of samples or below the limit of quantitation (LOQ). The LOQ was set as 2.0 standard deviations above the geometric mean of negative probes. The DESeq2 R package (v3.6.2), which employs negative binomial generalized linear models with median-of-ratios normalization, was used to identify differentially expressed genes (DEGs) between IgAN and each control group (35). Significant DEGs, defined as those with FDR < 0.10, were extracted for each comparison. Finally, the significant DEGs that show consistently high or low expression in IgAN across all comparisons were processed through the ToppGene Suite for Gene Ontology-based functional enrichment analysis, where a FDR < 0.05 was considered significant (36).

2.5 Serum and fecal acetate quantification by LC-MS

Acetate concentrations in serum and fecal samples were quantified using liquid chromatography–mass spectrometry (LC-MS) after derivatization with 3-nitrophenylhydrazine (3-NPH), as previously described (37). Fecal samples (~20 mg) were extracted in assay diluent buffer (water:acetonitrile, 1:1), centrifuged, and the supernatants were used for derivatization. For serum, 20 µL of sample was directly derivatized without prior extraction. Derivatized samples were incubated at 40 °C for 30 minutes and then diluted for analysis. LC-MS analysis was performed using a Q-Exactive Focus Orbitrap mass spectrometer (Thermo Fisher Scientific) coupled to a BEH C18 column (2.1 × 100 mm, 1.7 µm, Waters) under electrospray ionization in negative mode. Chromatographic separation was achieved using a water/acetonitrile gradient containing 0.1% formic acid at a flow rate of 0.35 mL/min. The mass spectrometer was operated in full scan mode (m/z 70–900), and acetate was identified and quantified by the retention time and m/z (194.1) of its derivatized form, which matched a derivatized acetate standard.

2.6 Statistical analysis and data visualization

All statistical analysis was performed in R (v4.4.0) except for gut microbial diversity analysis, which was performed on the QIIME2 platform using the q2-diversity plugin. For clinical characteristics, normality of continuous variables was assessed using histograms and skewness evaluations. For statistical comparisons between groups, Mann-Whitney U tests were used for two groups and Kruskal-Wallis was used for three or more groups. Linear associations between continuous variables were evaluated using the Pearson correlation coefficient. Permutational multivariate analysis of variance (PERMANOVA) with 999 permutations was used to evaluate differences in beta diversities among groups using the vegan R package (v2.7-1) (38). The Benjamini–Hochberg method was applied in calculations of false-discovery rates in cases of multiple testing. All boxplots and scatterplots were created with the ggpubr R package (v0.6.1) (39). Dotplot of the enriched Gene Ontology terms was visualized with the scToppR R package (v0.99.0) (40).

3 Results

3.1 Participant demographics and clinical profiles

A total of 336 participants were included in the gut microbiome analysis comprising 172 IgAN patients, 51 healthy donor controls, and 113 disease controls (Table 1). The IgAN patients had a mean age of 41.8, and 52% were female. The mean body mass index was 24.6 kg/m2 for IgAN, comparable to control groups, which ranged between 23.7 and 25.5 kg/m2. The mean eGFR ranged from 85.4 to 104.9 mL/min/1.73 m2 for each group.

Table 1
www.frontiersin.org

Table 1. Baseline clinical characteristics of study participants subject to gut microbial profiling.

For spatial transcriptomic analysis, 8 biopsy-proven IgAN cases were subject to spatial transcriptomics profiling along with 10 healthy controls and 35 disease controls (Table 2). The IgAN patients had a mean age of 36.6, younger than control groups. They were more female, while control groups had a male majority. Lower eGFRs were observed among other disease controls, with the DN group having the lowest mean eGFR at 56.3 mL/min/1.73 m2.

Table 2
www.frontiersin.org

Table 2. Baseline clinical characteristics of study participants subject to glomerular spatial transcriptomics profiling.

3.2 IgAN-specific functional enrichment of gut microbiota

A total of 466 microbial genera and 332 predicted metagenomic functional pathways were identified. While alpha diversities did not differ among groups (Figure 1A), beta diversities were significantly differed among diagnostic groups (PERMANOVA, p = 0.001) (Figure 1B). IgAN microbiota composition was distinct from that of healthy controls (p = 0.001), MCD (p = 0.034), and MN (p = 0.013) but not from DN (p = 0.34). Significant differences between IgAN and healthy controls were consistently observed across different beta diversity metrics including Jaccard indices (p = 0.001), unweighted UniFrac (p = 0.006), and weighted UniFrac (p = 0.001), as were the differences among the five diagnostic groups (Supplementary Figure 1). Beta diversities based on the Bray-Curtis metric also significantly differed among groups (p < 0.001) when major covariates (age, sex, and eGFR) were included in the model.

Figure 1
Three charts compare five groups: IgAN, HC, DN, MCD, and MN. Chart A shows box plots of the Shannon index with a p-value of 0.72. Chart B is a PCA plot with a p-value of 0.001, showing clusters for each group. Chart C shows box plots of pathway abundance. Each group is represented by different colors.

Figure 1. Characteristics of the gut microbiome in IgA nephropathy (IgAN). (A) Alpha diversity based on Shannon indices. The Kruskal-Wallis test was used to evaluate statistical significance. (HC: healthy controls, DN: diabetic nephropathy, MCD: minimal change disease, MN: membranous nephropathy) (B) Beta diversity based on Bray-Curtis dissimilarity indices represented as a principal coordinate analysis plot. Statistical significance was evaluated by permutational multivariate analysis of variance (PERMANOVA). (PC: principal coordinate) (C) Predicted abundances of the methanogenesis from acetate MetaCyc pathway grouped by diagnosis. Absolute pathways abundances were derived from PICRUSt2-predicted metagenome functions.

In taxonomic differential abundance analysis, eight genera (Bacteroides, Parabacteroides, Alistipes, Blautia, Romboutsia, Dorea, Lachnospira, and Butyricicoccus) differed between IgAN and healthy controls, but none were consistently different from all disease controls (Table 3). As for functional analysis, Methanogenesis from Acetate (METH-ACETATE-PWY) was the only pathway significantly enriched in IgAN relative to both healthy controls and to each of the different glomerular disease controls, as visualized in a plot of the predicted pathway abundances (Figure 1C).

Table 3
www.frontiersin.org

Table 3. Microbial genera with statistically significant differential abundance for IgA nephropathy compared to healthy and disease controls.

3.3 Altered acetate metabolism in IgAN

To assess the generation of acetate by the gut microbiota, we next compared relative acetate-producing potential based on the abundances of eight common gut microbial genera that contain species known to contribute to acetate production (Figure 2A). IgAN showed the highest mean total relative abundances of the eight genera at 39.8%, followed by disease controls at 37.3% (p = 0.037, FDR = 0.037) and healthy control at 30.6% (p < 0.001, FDR < 0.001). Among the 8 genera, Bacteroides and Blautia revealed higher abundance in IgAN than in healthy controls (Table 3). MCD (p = 0.008, FDR = 0.041) and MN (p = 0.014, FDR = 0.046) also showed elevated proportions of major acetate producers. Meanwhile, the acetotrophic methanogens recognized to possess the methanogenesis from acetate pathway, Methanosarcina and Methanothrix, were not detected in any sample. In addition to this, we found that IgAN patients had higher serum acetate levels than healthy controls (p = 0.03) but similar fecal acetate levels (Figure 2B). There was negligible correlation between serum and fecal acetate levels (R = -0.17, p = 0.15) (Figure 2C). The identified IgAN-associated gut microbial signatures, namely the relative acetate-producing potential and the methanogenesis from acetate pathway, did not show significant correlations with baseline clinical parameters or with serum and fecal acetate levels (Supplementary Table 1).

Figure 2
Three-panel figure showing microbiome and metabolite data. Panel A presents a stacked bar graph and bar chart of microbiome genus relative abundance across five groups: IgAN, HC, DN, MCD, and MN. Panel B includes box plots comparing serum and fecal acetate levels between IgAN and HC groups, with significant p-value for serum acetate. Panel C displays a scatter plot correlating serum and fecal acetate, with IgAN and HC differentiation, showing a weak correlation. Legends indicate genus and group color codes.

Figure 2. Acetate metabolism in IgA nephropathy (IgAN). (A) Relative abundances of major acetate-producing microbial taxa. The bar plot depicts the relative abundances of major acetate-producing genera for each sample, compiled into mean relative abundances for each diagnosis group. Pairwise Mann-Whitney U tests were performed for the combined mean relative abundances, and asterisks indicate statistically significant differences based on false discovery rates (FDR): *FDR < 0.05, ***FDR < 0.001. (HC: healthy controls, DN: diabetic nephropathy, MCD: minimal change disease, MN: membranous nephropathy) (B) Serum and fecal acetate levels in IgAN (N = 55) and HC (N = 23). Measurements from liquid chromatography–mass spectrometry are shown as intensities in arbitrary units (a.u.), and statistical significance was evaluated with Mann-Whitney U tests. (C) Correlation between measured serum and fecal acetate levels. Pearson correlation coefficient is shown, along with the corresponding p-value.

3.4 Glomerular spatial transcriptomics reveals altered signaling pathways in IgAN

Glomerular expression profiles of 17,834 genes were subject to analysis. Relative log expression plots showed no extreme outliers in count distributions (Supplementary Figure 2A), and compartment-specific gene expressions supported appropriate region selection (Supplementary Figure 2B). A total of 1,227 and 1,078 genes were consistently upregulated and downregulated, respectively, in IgAN relative to healthy controls and to each disease control group (Figure 3A, Supplementary Tables 2, 3). The highest fold increase was found in COL3A1 (collagen type III alpha 1 chain), a key component of the extracellular matrix. The highly expressed DEGs were associated with 1,209 Gene Ontology terms (Supplementary Table 4), where the top categories included cell adhesion, extracellular matrix organization, protein complex assembly, and mitochondrial function (Table 4). Meanwhile, the lowly expressed DEGs were enriched in 18 functional terms (Table 4, Supplementary Table 5). GPCR activity was the top enriched function, alongside terms pertaining to olfactory receptors, which belong to the GPCR superfamily. Beta-1,3-galactosyltransferase activity, involved in O-galactosylation of glycoproteins including immunoglobulins, also showed significant enrichment. Moreover, the majority of SCFA transmembrane transporter and calcium, potassium:sodium antiporter genes were downregulated in IgAN (Figure 3B). As for the acetate-sensing GPCRs, IgAN displayed undetectable GPR41 (FFAR3) expression as well as lower GPR43 (FFAR2) and Olfr78 (OR51E2) expressions, although the latter was only significant relative to MCD and MN (Figure 4).

Figure 3
A) Volcano plot showing gene expression data with log2 fold change on the x-axis and negative log10 adjusted p-value on the y-axis. Genes are indicated by dots, with blue and red highlights for significant genes. B) Dot plot displaying gene ontology enrichment with gene ratio on the x-axis and terms on the y-axis. Dot colors represent the false discovery rate, and dot sizes show the number of genes in each gene set. Biological process terms are in orange, and molecular function terms are in blue.

Figure 3. Differential gene expression analysis in the glomerular transcriptome of IgA nephropathy (IgAN). (A) Volcano plot of differentially expressed genes (DEGs) for the glomerulus of IgAN compared to healthy controls. Red and blue colors indicate DEGs with log2 fold change above 1.5 or below -1.5, respectively. (B) Gene Ontology (GO) terms significantly enriched among DEGs with decreased expression in IgAN. The x-axis, gene ratio, shows the ratio between the count of contributing DEGs and the total gene count for each GO term. Dot sizes show the absolute number of DEGs contributing to each GO term, while the dot colors correspond to the false discovery rate (FDR).

Table 4
www.frontiersin.org

Table 4. Top gene ontologies for differentially expressed genes (DEGs) consistently upregulated or downregulated in IgA nephropathy.

Figure 4
Box plots displaying expression levels of GPR41 (FFAR3), GPR43 (FFAR2), and Olfr78 (OR51E2) genes across different groups: IgAN, HC, DN, MCD, and MN. GPR41 shows multiple significant differences, marked by asterisks, indicating varying expression levels. GPR43 and Olfr78 also show statistical significance across the groups. Counts per million are plotted on the y-axis, with significant differences denoted by asterisks (*, **, ***) indicating p-values.

Figure 4. Glomerular expression of short-chain fatty acid-sensing G protein-coupled receptors. Statistically significant differences based on pairwise Mann-Whitney U tests are indicated with asterisks: *p < 0.05, **p < 0.01, ***p < 0.001, ***p < 0.0001.

4 Discussion

In this study, we integrated gut microbiome profiling and glomerular spatial transcriptomics to characterize the microbial and transcriptomic changes in IgAN that potentially modulate the gut-kidney axis. The gut microbiome in IgAN showed an enrichment of major acetate-producing taxa and the methanogenesis from acetate pathway, accompanied by elevated serum acetate levels. At the same time, the glomerular transcriptome of IgAN demonstrated a downregulation of SCFA transporters and SCFA-sensing GPCRs, suggesting a mismatch between systemic metabolite availability and kidney sensing capacity.

Mounting evidence suggests acetate plays a key role in regulating immunity and inflammation, including in kidney diseases. Primarily generated through gut bacterial fermentation, acetate is the predominant SCFA, constituting 60-75% of total gut SCFAs (41, 42) and over 90% of those in serum (43). As both an energy source and a signaling molecule, it supports maintenance of the intestinal epithelial barrier while suppressing proinflammatory mediators (41). In vivo, acetate prevented acute kidney injury in an ischemia-reperfusion injury model and reduced kidney fibrosis in an unilateral ureteral obstruction model (44, 45), and it also attenuated clinical manifestations in a mouse model of IgAN (46). Hence, acetate may have renoprotective effects specific to IgAN as well as across kidney disease more broadly.

Our gut microbial analysis suggests that acetate production and consumption in the gut are both modulated in IgAN. Functional analysis indicated an increase in predicted enzymatic activities attributable to acetate consumption by methanogens. While methanogens have been associated with several diseases including inflammatory bowel disease (47), research on their role in the human gut is relatively sparse. Since methanogens were not identified in our taxonomic annotations, possibly because archaea comprise only 0.1% of reads in non-targeted sequencing (48), the observed metabolic activities may also reflect contributions from other acetate utilizers. As for acetate production, a previous study reported significantly decreased fecal SCFA levels in IgAN, with correlations to gut microbiota composition and clinical parameters, suggesting that disruption of SCFA levels could be relevant to IgAN (49). However, we did not observe differences in fecal acetate levels, and the relative increase in major acetate-producing genera in IgAN suggests that acetate production may not be compromised in IgAN. Our results are therefore most consistent with a simultaneous increase in the consumption and production of gut acetate in IgAN.

In addition, the glomerular transcriptome of IgAN shows potential evidence of altered acetate signaling in the kidney. Acetate in circulation can signal across diverse tissues through specific GPCRs (GPR41, GPR43, Olfr78) (50), which were downregulated in the glomerular transcriptome of IgAN: GPR41 was undetectable, and GPR43 and Olfr78 had significantly lower expression compared to MCD and MN. As GPR41 expression was also undetected in some of the healthy controls, the observed lack of GPR41 expression in IgAN likely reflects its low baseline expression below the detection threshold. We also observed a downregulation of SCFA transmembrane transporters, namely SMCT (SLC5A8), OAT4 (SLC22A9), OAT5 (SLC22A10), and UST5 (SLC22A25). These transporters are mainly associated with renal tubules rather than the glomerulus, but they are also part of the organic anion transporter family, whose members can mediate interorgan communication by modulating metabolite levels (51). The downregulation of olfactory receptor activity in IgAN could also reflect decreased SCFA signaling, as some renal olfactory receptors have physiological functions. Olfr78, as a prime example, localizes to the juxtaglomerular apparatus and modulates renin homeostasis through acetate and propionate signaling (52). Other downregulated orphan olfactory receptors may similarly localize to the glomerulus and contribute to acetate signaling. While these findings may indicate compensatory responses to elevated acetate levels, negative feedback is not universal. Indeed, a recent study showed that acetate exposure upregulated, rather than downregulated, GPR43 expression in kidney tubular cells (53). Hence, the glomerular downregulation of SCFA-related genes in IgAN may reflect pathologic suppression of acetate signaling, which would leave the kidney vulnerable to inflammation and oxidative stress.

Collectively, our analyses present a nuanced picture of acetate biology in IgAN that involves potential changes from synthesis in the gut to signaling in the glomerulus. While serum acetate levels were elevated in IgAN, this may be a reactive change rather than a pathological mechanism. Elevations in circulating acetate levels have been observed in acute stress conditions such as bacterial infections, where acetate facilitates memory T cell activity (54). Thus, the absolute increase in serum acetate levels in IgAN may not meet physiologic needs in acute inflammatory conditions, exacerbated by the downregulation of potential SCFA sensors in the glomerulus. The lack of correlation between fecal and serum acetate levels further demonstrates the complexity of acetate dynamics in IgAN. Previous studies in healthy individuals also did not report clear correlations between the two measurements despite rapid absorption of fecal acetate (43, 55).

In the kidney, acetate signaling may modulate diverse molecular pathways, with potential links to pathophysiologic processes underlying glomerular diseases. One possibly relevant mechanism is the nuclear transcription factor kappa B (NF-κB) pathway, a key regulator of inflammation across glomerular diseases and CKD (5658). SCFA-activated GPR43 inhibits the NF-κB pathway and oxidative stress in DN, suggesting protective actions of acetate (59). At the same time, SCFAs can also amplify signaling from aryl hydrocarbon receptors, which can stimulate NF-κB-driven inflammation and suppress Nrf2 antioxidant pathways (60, 61). Acetate may also affect the kidney by modulating the renin-angiotensin system. Given the increased circulating acetate levels in IgAN, Olfr78 in the juxtaglomerular apparatus could be highly activated, promoting renin secretion and leading to glomerular hypertension. Indeed, a study of early DN identified elevated serum acetate levels in a rat model and suggested that excess acetate may promote kidney injury through the renin-angiotensin system (62). Also downstream of the renin-angiotensin system is the Wnt/β-catenin pathway, whose persistent activation can lead to kidney fibrosis in glomerular diseases and CKD (63, 64). Conversely, the counterbalancing effects of GPR41 on blood pressure, together with the downregulation of Olfr78, could blunt or even reverse these effects on the regulation of renin (65). As such, renal acetate signaling in IgAN may affect multiple molecular mechanisms, whose contributions to pathophysiology are likely context dependent.

Beyond its effects on glomerular inflammation, acetate may also contribute to the pathophysiology of IgAN through signaling in other cell populations. Although not directly addressed by the current analysis, earlier in vivo studies have suggested that acetate can regulate gut mucosal immune cells. Specifically, acetate supplementation induced gut IgA production and regulated IgA reactivity to microbes through stimulation of dendritic cells, which promote IgA class switching and gut homing in B cells (66). In certain contexts, acetate can also induce dendritic cells to produce B-cell activating factor (BAFF) (66, 67), an emerging therapeutic target whose levels are elevated in IgAN (14). In addition, acetate stimulates CD4+CCR6+ T follicular helper-like cells by promoting CCL20 production in colonic epithelial cells, regulating T cell-dependent IgA production in germinal centers by modulating Toll-like receptor (TLR) signaling (68). This may also be relevant to IgAN pathogenesis, as a recent study found that transplanting gut microbiota from IgAN patients into mice activates TLR4 signaling and induces the disease phenotype, while TLR4 inhibition suppresses Gd-IgA1 production (69). In all, acetate can modulate IgA biology through its interactions with gut immune cells, such as dendritic cells and helper T cells, which may have potential relevance to the pathogenesis of IgAN.

Recent studies associate gut dysbiosis and metabolite dysregulation with kidney injury across glomerular diseases. Uremic toxin generation and impaired SCFA production, along with intestinal barrier dysfunction, have been identified as major pathways whereby gut dysbiosis contribute to CKD (70). The gut microbiome in MN shows reductions in Lactobacillus, whose tryptophan-derived metabolites can attenuate kidney damage by inhibition of aryl hydrocarbon receptor signaling (71). Significant gut dysbiosis was also observed in MCD that includes decreased abundance of the butyrate producer Faecalibacterium (72). As for DN, the gut microbiome is associated with multiple potentially pathogenic mechanisms including increased endotoxin levels, decreased SCFA levels, and dysregulated bile acid metabolism (73). As such, SCFAs are associated with multiple kidney diseases, typically showing reduced production, but their mechanistic roles remain incompletely defined.

SCFAs, in particular acetate, may therefore play a disease-specific role in IgAN. Our analysis suggests an overall increased activity of acetate metabolism in the gut, where it can boost IgA production and reactivity as previously discussed. At the same time, downregulation of genes responsible for SCFA transport and sensing in the glomerulus would limit acetate uptake and signaling, which can contribute to glomerular injury through impaired regulation of anti-inflammatory pathways, disruption of renin signaling, and promotion of local metabolic stress. The observed elevation in serum acetate may represent both a systemic response to inflammation and a compensatory adaptation to defective acetate signaling in the kidney.

Human studies demonstrate the clinical relevance of gut dysbiosis in IgAN, but the contribution of SCFAs such as acetate remains to be established. Gut microbiota profile could predict treatment response in IgAN (74), and a pilot trial suggested fecal microbiota transplantation (FMT) may serve as an adjunctive treatment (75). Nonetheless, despite several reports of decreased serum and fecal acetate levels in IgAN and in broader CKD (49, 76), SCFAs were not among the intestinal metabolites altered by FMT in IgAN in the pilot study. Adding to this complexity, our glomerular transcriptomic data indicated that acetate signaling in the kidney may depend on systemic availability as well as renal sensing capacity, highlighting the need to assess associated receptor functions when interpreting SCFA biology in kidney disease.

The current study has several limitations. First, the 16S rRNA data had limited resolution. Taxonomic annotation was restricted to the genus level, and certain taxa may not have been detected, particularly archaeal species for which the primer designs and taxonomic databases are less optimized. Functional annotation was also performed indirectly through PICRUSt2, which first derives the predicted metagenome from the 16S data and then algorithmically infers pathway composition. Further validation through metagenomic sequencing and functional assays is needed to verify the observed genus- and pathway-level associations. Second, due to the relatively low number of reads in our spatial transcriptomics data, we could not reliably localize gene expression patterns within the glomerulus, and some genes with low expression levels may have not been detected, as was likely the case for GPR41. Third, while efforts were made to select representative cases, our study may not reflect the entire spectrum of IgAN. The cases included in this study are all relatively young Koreans, which may not fully capture the nature of the disease in other ethnicities, especially given the known differences in clinical severity and gender distributions in Asian-Pacific populations (14). The low number of samples subject to glomerular transcriptomics also limits the generalizability of our findings. Finally, because this is a data-driven, observational study, the contributions of metabolites, glomerular genes, and gut microbial functions to IgAN pathophysiology remain to be established by experimental studies. Other SCFAs may also signal through the involved GPCRs and contribute to the postulated mechanisms, and additional confounding factors, such as dietary patterns, could influence gut microbial composition as well as SCFA levels. There may also be important transcriptional and microbial characteristics not captured by our analysis due to the limited sample size and the study design, which compared IgAN with four different control groups separately to minimize false positive results.

In conclusion, we have identified elevations in gut microbial acetate production and consumption activities as well as glomerular downregulation of GPCR, SCFA transport, and galactosylation functions in IgAN. Our results suggest that acetate may contribute to the pathophysiology of IgAN by mediating the gut-kidney interaction through GPCR signaling.

Data availability statement

The data presented in the study are deposited in the NCBI SRA repository, accession number PRJNA1399494.

Ethics statement

The studies involving humans were approved by Seoul National University Hospital. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.

Author contributions

JHK: Writing – original draft, Writing – review & editing, Data curation, Formal Analysis, Investigation, Methodology, Software, Visualization. SP: Writing – original draft, Writing – review & editing, Conceptualization, Data curation, Formal Analysis, Funding acquisition, Investigation, Methodology, Project administration, Software, Visualization. MK: Writing – review & editing, Formal Analysis, Investigation, Methodology, Resources, Software. JP: Writing – review & editing, Conceptualization, Data curation, Investigation, Resources. JL: Writing – review & editing, Conceptualization, Data curation, Investigation, Resources. HC: Writing – review & editing, Conceptualization, Data curation, Investigation, Resources. JEK: Writing – review & editing, Conceptualization, Data curation, Investigation, Resources. HN: Writing – review & editing, Formal Analysis, Investigation, Methodology, Resources, Software. SP: Writing – review & editing, Conceptualization, Formal Analysis, Investigation, Methodology, Resources. KM: Writing – review & editing, Conceptualization, Data curation, Investigation, Methodology, Resources. HK: Writing – review & editing, Conceptualization, Investigation, Methodology, Resources, Software. YK: Writing – review & editing, Conceptualization, Data curation, Investigation, Methodology, Resources. DKK: Writing – review & editing, Conceptualization, Data curation, Investigation, Methodology, Resources. HL: Writing – review & editing, Writing – original draft, Conceptualization, Data curation, Formal Analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Supervision. DK: Writing – review & editing, Formal Analysis, Investigation, Methodology, Resources. ML: Writing – review & editing, Formal Analysis, Investigation, Methodology, Resources.

Funding

The author(s) declared that financial support was received for work and/or its publication. This research was supported by the National Research Foundation of Korea (NRF) funded by the Ministry of Science, ICT, and Future Planning (grant numbers RS-2019-NR039661, RS-2022-NR070373, RS-2023-00219548, and RS-2024-00345867). This research was also supported by a grant of the Korea Health Technology R&D Project through the Korea Health Industry Development Institute (KHIDI), funded by the Ministry of Health & Welfare, Republic of Korea (grant number: RS-2024-00403375).

Acknowledgments

The Biospecimens and data used for this study were provided by the Biobank of Seoul National University Hospital, a member of the Korea Biobank Network (project No. 2024ER050800).

Conflict of interest

The authors declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declared that generative AI was not used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2025.1665585/full#supplementary-material

References

1. Li Z, He R, Wang Y, Qu Z, Liu J, Yu R, et al. Global trends of chronic kidney disease from 1990 to 2021: A systematic analysis for the global burden of disease study 2021. BMC Nephrol. (2025) 26:385. doi: 10.1186/s12882-025-04309-7

PubMed Abstract | Crossref Full Text | Google Scholar

2. Zhu J, Fu Y, Olovo CV, Xu J, Wu Q, Wei W, et al. The influence of gut microbiota on the gut-brain-kidney axis and its implications for chronic kidney disease. Front Microbiol. (2025) 16:1535356. doi: 10.3389/fmicb.2025.1535356

PubMed Abstract | Crossref Full Text | Google Scholar

3. Tao P, Huo J, and Chen L. Bibliometric analysis of the relationship between gut microbiota and chronic kidney disease from 2001–2022. Integr Med Nephrol Andrology. (2024) 11:e00017. doi: 10.1097/imna-d-23-00017

Crossref Full Text | Google Scholar

4. Miao H, Liu F, Wang YN, Yu XY, Zhuang S, Guo Y, et al. Targeting lactobacillus johnsonii to reverse chronic kidney disease. Signal Transduct Target Ther. (2024) 9:195. doi: 10.1038/s41392-024-01913-1

PubMed Abstract | Crossref Full Text | Google Scholar

5. Zhang L, Hu L, Tan L, Zhang Z, Chen M, Gan W, et al. The dysbiosis of gut microbiota and dysregulation of metabolites in iga nephropathy and membranous nephropathy. Front Med (Lausanne). (2025) 12:1618947. doi: 10.3389/fmed.2025.1618947

PubMed Abstract | Crossref Full Text | Google Scholar

6. Wang X, Liu X, Gong F, Jiang Y, Zhang C, Zhou W, et al. Targeting gut microbiota for diabetic nephropathy treatment: probiotics, dietary interventions, and fecal microbiota transplantation. Front Endocrinol (Lausanne). (2025) 16:1621968. doi: 10.3389/fendo.2025.1621968

PubMed Abstract | Crossref Full Text | Google Scholar

7. Lai Y, Zhu Y, Zhang X, Ding S, Wang F, Hao J, et al. Gut microbiota-derived metabolites: potential targets for cardiorenal syndrome. Pharmacol Res. (2025) 214:107672. doi: 10.1016/j.phrs.2025.107672

PubMed Abstract | Crossref Full Text | Google Scholar

8. Evenepoel P, Stenvinkel P, Shanahan C, and Pacifici R. Inflammation and gut dysbiosis as drivers of ckd-mbd. Nat Rev Nephrol. (2023) 19:646–57. doi: 10.1038/s41581-023-00736-7

PubMed Abstract | Crossref Full Text | Google Scholar

9. Corradi V, Caprara C, Barzon E, Mattarollo C, Zanetti F, Ferrari F, et al. A possible role of P-cresyl sulfate and indoxyl sulfate as biomarkers in the prediction of renal function according to the gfr (G) categories. Integr Med Nephrol Andrology. (2024) 11:e24–00002. doi: 10.1097/imna-d-24-00002

Crossref Full Text | Google Scholar

10. Miao H, Zhang SJ, Wu X, Li P, and Zhao YY. Tryptophan metabolism as a target in gut microbiota, ageing and kidney disease. Int J Biol Sci. (2025) 21:4374–87. doi: 10.7150/ijbs.115359

PubMed Abstract | Crossref Full Text | Google Scholar

11. Cao BN, Zhang CY, Wang Z, and Wang YX. Causal relationship between 412 gut microbiota, 1,400 blood metabolites, and diabetic nephropathy: A randomized mendelian study. Front Endocrinol (Lausanne). (2024) 15:1450428. doi: 10.3389/fendo.2024.1450428

PubMed Abstract | Crossref Full Text | Google Scholar

12. Song S, Ning L, and Yu J. Elucidating the causal relationship between gut microbiota, metabolites, and diabetic nephropathy in european patients: revelations from genome-wide bidirectional mendelian randomization analysis. Front Endocrinol (Lausanne). (2024) 15:1391891. doi: 10.3389/fendo.2024.1391891

PubMed Abstract | Crossref Full Text | Google Scholar

13. Gesualdo L, Di Leo V, and Coppo R. The mucosal immune system and iga nephropathy. Semin Immunopathol. (2021) 43:657–68. doi: 10.1007/s00281-021-00871-y

PubMed Abstract | Crossref Full Text | Google Scholar

14. Cheung CK, Alexander S, Reich HN, Selvaskandan H, Zhang H, and Barratt J. The pathogenesis of iga nephropathy and implications for treatment. Nat Rev Nephrol. (2025) 21:9–23. doi: 10.1038/s41581-024-00885-3

PubMed Abstract | Crossref Full Text | Google Scholar

15. Kiryluk K, Sanchez-Rodriguez E, Zhou XJ, Zanoni F, Liu L, Mladkova N, et al. Genome-wide association analyses define pathogenic signaling pathways and prioritize drug targets for iga nephropathy. Nat Genet. (2023) 55:1091–105. doi: 10.1038/s41588-023-01422-x

PubMed Abstract | Crossref Full Text | Google Scholar

16. Kiryluk K, Li Y, Scolari F, Sanna-Cherchi S, Choi M, Verbitsky M, et al. Discovery of new risk loci for iga nephropathy implicates genes involved in immunity against intestinal pathogens. Nat Genet. (2014) 46:1187–96. doi: 10.1038/ng.3118

PubMed Abstract | Crossref Full Text | Google Scholar

17. Lafayette R, Kristensen J, Stone A, Floege J, Tesar V, Trimarchi H, et al. Efficacy and safety of a targeted-release formulation of budesonide in patients with primary iga nephropathy (Nefigard): 2-year results from a randomised phase 3 trial. Lancet. (2023) 402:859–70. doi: 10.1016/S0140-6736(23)01554-4

PubMed Abstract | Crossref Full Text | Google Scholar

18. Feriozzi S and Polci R. The role of tonsillectomy in iga nephropathy. J Nephrol. (2016) 29:13–9. doi: 10.1007/s40620-015-0247-4

PubMed Abstract | Crossref Full Text | Google Scholar

19. Dong Z, Zhang R, Shen L, Ji HF, He H, Ji X, et al. Gut microbiota and immunoglobulin a nephropathy: exploration of dietary intervention and treatment strategies. Food Sci Nutr. (2025) 13:e70218. doi: 10.1002/fsn3.70218

PubMed Abstract | Crossref Full Text | Google Scholar

20. Pattrapornpisut P, Avila-Casado C, and Reich HN. Iga nephropathy: core curriculum 2021. Am J Kidney Dis. (2021) 78:429–41. doi: 10.1053/j.ajkd.2021.01.024

PubMed Abstract | Crossref Full Text | Google Scholar

21. Luvizotto MJ, Menezes-Silva L, Woronik V, Monteiro RC, and Camara NOS. Gut-kidney axis in iga nephropathy: role on mesangial cell metabolism and inflammation. Front Cell Dev Biol. (2022) 10:993716. doi: 10.3389/fcell.2022.993716

PubMed Abstract | Crossref Full Text | Google Scholar

22. Popova A, Slisere B, Racenis K, Kuzema V, Karklins R, Saulite M, et al. Iga class-switched cd27-cd21+ B cells in iga nephropathy. Nephrol Dial Transplant. (2025) 40:505–15. doi: 10.1093/ndt/gfae173

PubMed Abstract | Crossref Full Text | Google Scholar

23. Sanchez-Russo L, Rajasekaran A, Bin S, Faith J, and Cravedi P. The gut and kidney crosstalk in immunoglobulin a nephropathy. Kidney360. (2022) 3:1630–9. doi: 10.34067/KID.0002382022

PubMed Abstract | Crossref Full Text | Google Scholar

24. Zhao J, Bai M, Ning X, Qin Y, Wang Y, Yu Z, et al. Expansion of escherichia-shigella in gut is associated with the onset and response to immunosuppressive therapy of iga nephropathy. J Am Soc Nephrol. (2022) 33:2276–92. doi: 10.1681/ASN.2022020189

PubMed Abstract | Crossref Full Text | Google Scholar

25. Park JI, Kim TY, Oh B, Cho H, Kim JE, Yoo SH, et al. Comparative analysis of the tonsillar microbiota in iga nephropathy and other glomerular diseases. Sci Rep. (2020) 10:16206. doi: 10.1038/s41598-020-73035-x

PubMed Abstract | Crossref Full Text | Google Scholar

26. Jain S and Eadon MT. Spatial transcriptomics in health and disease. Nat Rev Nephrol. (2024) 20:659–71. doi: 10.1038/s41581-024-00841-1

PubMed Abstract | Crossref Full Text | Google Scholar

27. Park S, Kang M, Kim YC, Kim DK, Oh KH, Joo KW, et al. Glomerular spatial transcriptomics of iga nephropathy according to the presence of mesangial proliferation. Sci Rep. (2024) 14:2211. doi: 10.1038/s41598-024-52581-8

PubMed Abstract | Crossref Full Text | Google Scholar

28. Kang E, Kim Y, Kim YC, Kim E, Lee N, Kim Y, et al. Biobanking for glomerular diseases: A study design and protocol for korea renal biobank network system toward next-generation analysis (Kornerstone). BMC Nephrol. (2020) 21:367. doi: 10.1186/s12882-020-02016-z

PubMed Abstract | Crossref Full Text | Google Scholar

29. Kim JE, Kim HE, Park JI, Cho H, Kwak MJ, Kim BY, et al. The association between gut microbiota and uremia of chronic kidney disease. Microorganisms. (2020) 8:907. doi: 10.3390/microorganisms8060907

PubMed Abstract | Crossref Full Text | Google Scholar

30. Callahan BJ, McMurdie PJ, Rosen MJ, Han AW, Johnson AJ, and Holmes SP. Dada2: high-resolution sample inference from illumina amplicon data. Nat Methods. (2016) 13:581–3. doi: 10.1038/nmeth.3869

PubMed Abstract | Crossref Full Text | Google Scholar

31. Glockner FO, Yilmaz P, Quast C, Gerken J, Beccati A, Ciuprina A, et al. 25 years of serving the community with ribosomal rna gene reference databases and tools. J Biotechnol. (2017) 261:169–76. doi: 10.1016/j.jbiotec.2017.06.1198

PubMed Abstract | Crossref Full Text | Google Scholar

32. Douglas GM, Maffei VJ, Zaneveld JR, Yurgel SN, Brown JR, Taylor CM, et al. Picrust2 for prediction of metagenome functions. Nat Biotechnol. (2020) 38:685–8. doi: 10.1038/s41587-020-0548-6

PubMed Abstract | Crossref Full Text | Google Scholar

33. Caspi R, Billington R, Ferrer L, Foerster H, Fulcher CA, Keseler IM, et al. The metacyc database of metabolic pathways and enzymes and the biocyc collection of pathway/genome databases. Nucleic Acids Res. (2016) 44:D471–80. doi: 10.1093/nar/gkv1164

PubMed Abstract | Crossref Full Text | Google Scholar

34. Fernandes AD, Macklaim JM, Linn TG, Reid G, and Gloor GB. Anova-like differential expression (Aldex) analysis for mixed population rna-seq. PloS One. (2013) 8:e67019. doi: 10.1371/journal.pone.0067019

PubMed Abstract | Crossref Full Text | Google Scholar

35. Love MI, Huber W, and Anders S. Moderated estimation of fold change and dispersion for rna-seq data with deseq2. Genome Biol. (2014) 15:550. doi: 10.1186/s13059-014-0550-8

PubMed Abstract | Crossref Full Text | Google Scholar

36. Chen J, Bardes EE, Aronow BJ, and Jegga AG. Toppgene suite for gene list enrichment analysis and candidate gene prioritization. Nucleic Acids Res. (2009) 37:W305–11. doi: 10.1093/nar/gkp427

PubMed Abstract | Crossref Full Text | Google Scholar

37. Meng X, Pang H, Sun F, Jin X, Wang B, Yao K, et al. Simultaneous 3-nitrophenylhydrazine derivatization strategy of carbonyl, carboxyl and phosphoryl submetabolome for lc-ms/ms-based targeted metabolomics with improved sensitivity and coverage. Anal Chem. (2021) 93:10075–83. doi: 10.1021/acs.analchem.1c00767

PubMed Abstract | Crossref Full Text | Google Scholar

38. Oksanen J, Simpson G, Blanchet F, Kindt R, Legendre P, Minchin P, et al. Vegan: community ecology package. (2025).

Google Scholar

39. Kassambara A. Ggpubr: 'Ggplot2' Based publication ready plots. (2025).

Google Scholar

40. Granger B and Berto S. Sctoppr: A coding-friendly R interface to toppgene. Bioinformatics. (2024) 40. doi: 10.1093/bioinformatics/btae582

PubMed Abstract | Crossref Full Text | Google Scholar

41. Hosmer J, McEwan AG, and Kappler U. Bacterial acetate metabolism and its influence on human epithelia. Emerg Top Life Sci. (2024) 8:1–13. doi: 10.1042/ETLS20220092

PubMed Abstract | Crossref Full Text | Google Scholar

42. Zhang D, Jian YP, Zhang YN, Li Y, Gu LT, Sun HH, et al. Short-chain fatty acids in diseases. Cell Commun Signal. (2023) 21:212. doi: 10.1186/s12964-023-01219-9

PubMed Abstract | Crossref Full Text | Google Scholar

43. Manokasemsan W, Jariyasopit N, Poungsombat P, Kaewnarin K, Wanichthanarak K, Kurilung A, et al. Quantifying fecal and plasma short-chain fatty acids in healthy thai individuals. Comput Struct Biotechnol J. (2024) 23:2163–72. doi: 10.1016/j.csbj.2024.05.007

PubMed Abstract | Crossref Full Text | Google Scholar

44. Andrade-Oliveira V, Amano MT, Correa-Costa M, Castoldi A, Felizardo RJ, de Almeida DC, et al. Gut bacteria products prevent aki induced by ischemia-reperfusion. J Am Soc Nephrol. (2015) 26:1877–88. doi: 10.1681/ASN.2014030288

PubMed Abstract | Crossref Full Text | Google Scholar

45. Kawabata C, Hirakawa Y, Inagi R, and Nangaku M. Acetate attenuates kidney fibrosis in an oxidative stress-dependent manner. Physiol Rep. (2023) 11:e15774. doi: 10.14814/phy2.15774

PubMed Abstract | Crossref Full Text | Google Scholar

46. Tan J, Dong L, Jiang Z, Tan L, Luo X, Pei G, et al. Probiotics ameliorate iga nephropathy by improving gut dysbiosis and blunting nlrp3 signaling. J Transl Med. (2022) 20:382. doi: 10.1186/s12967-022-03585-3

PubMed Abstract | Crossref Full Text | Google Scholar

47. Chaudhary PP, Conway PL, and Schlundt J. Methanogens in humans: potentially beneficial or harmful for health. Appl Microbiol Biotechnol. (2018) 102:3095–104. doi: 10.1007/s00253-018-8871-2

PubMed Abstract | Crossref Full Text | Google Scholar

48. Borrel G, Brugere JF, Gribaldo S, Schmitz RA, and Moissl-Eichinger C. The host-associated archaeome. Nat Rev Microbiol. (2020) 18:622–36. doi: 10.1038/s41579-020-0407-y

PubMed Abstract | Crossref Full Text | Google Scholar

49. Chai L, Luo Q, Cai K, Wang K, and Xu B. Reduced fecal short-chain fatty acids levels and the relationship with gut microbiota in iga nephropathy. BMC Nephrol. (2021) 22:209. doi: 10.1186/s12882-021-02414-x

PubMed Abstract | Crossref Full Text | Google Scholar

50. Foresto-Neto O, Ghirotto B, and Camara NOS. Renal sensing of bacterial metabolites in the gut-kidney axis. Kidney360. (2021) 2:1501–9. doi: 10.34067/KID.0000292021

PubMed Abstract | Crossref Full Text | Google Scholar

51. Nigam SK, Bush KT, Martovetsky G, Ahn SY, Liu HC, Richard E, et al. The organic anion transporter (Oat) family: A systems biology perspective. Physiol Rev. (2015) 95:83–123. doi: 10.1152/physrev.00025.2013

PubMed Abstract | Crossref Full Text | Google Scholar

52. Shepard BD and Pluznick JL. How does your kidney smell? Emerging roles for olfactory receptors in renal function. Pediatr Nephrol. (2016) 31:715–23. doi: 10.1007/s00467-015-3181-8

PubMed Abstract | Crossref Full Text | Google Scholar

53. Wang H, Kang T, and Li W. Reduced maternal scfas in gdm diminish gpr43 signaling and induce offspring cakut. Commun Biol. (2025) 8:1063. doi: 10.1038/s42003-025-08469-y

PubMed Abstract | Crossref Full Text | Google Scholar

54. Balmer ML, Ma EH, Bantug GR, Grahlert J, Pfister S, Glatter T, et al. Memory cd8(+) T cells require increased concentrations of acetate induced by stress for optimal function. Immunity. (2016) 44:1312–24. doi: 10.1016/j.immuni.2016.03.016

PubMed Abstract | Crossref Full Text | Google Scholar

55. Yamamura R, Nakamura K, Kitada N, Aizawa T, Shimizu Y, Nakamura K, et al. Associations of gut microbiota, dietary intake, and serum short-chain fatty acids with fecal short-chain fatty acids. Biosci Microbiota Food Health. (2020) 39:11–7. doi: 10.12938/bmfh.19-010

PubMed Abstract | Crossref Full Text | Google Scholar

56. Chang M, Shi X, Yang B, Li P, Zhang Y, Zhang Q, et al. Modified huangqi chifeng decoction alleviates podocyte injury on rat with experimental membranous nephropathy. Ren Fail. (2025) 47:2459896. doi: 10.1080/0886022X.2025.2459896

PubMed Abstract | Crossref Full Text | Google Scholar

57. Tian R, Wang L, Chen A, Huang L, Liang X, Wang R, et al. Sanqi oral solution ameliorates renal damage and restores podocyte injury in experimental membranous nephropathy via suppression of nfkappab. BioMed Pharmacother. (2019) 115:108904. doi: 10.1016/j.biopha.2019.108904

PubMed Abstract | Crossref Full Text | Google Scholar

58. Yuan Q, Tang B, and Zhang C. Signaling pathways of chronic kidney diseases, implications for therapeutics. Signal Transduct Target Ther. (2022) 7:182. doi: 10.1038/s41392-022-01036-5

PubMed Abstract | Crossref Full Text | Google Scholar

59. Huang W, Man Y, Gao C, Zhou L, Gu J, Xu H, et al. Short-chain fatty acids ameliorate diabetic nephropathy via gpr43-mediated inhibition of oxidative stress and nf-kappab signaling. Oxid Med Cell Longev. (2020) 2020:4074832. doi: 10.1155/2020/4074832

PubMed Abstract | Crossref Full Text | Google Scholar

60. Zheng W, Liu M, Lv X, He C, Yin J, and Ma J. Ahr governs lipid metabolism: the role of gut microbiota. Front Microbiol. (2025) 16:1442282. doi: 10.3389/fmicb.2025.1442282

PubMed Abstract | Crossref Full Text | Google Scholar

61. Wang YN, Miao H, Yu XY, Guo Y, Su W, Liu F, et al. Oxidative stress and inflammation are mediated via aryl hydrocarbon receptor signalling in idiopathic membranous nephropathy. Free Radic Biol Med. (2023) 207:89–106. doi: 10.1016/j.freeradbiomed.2023

PubMed Abstract | Crossref Full Text | Google Scholar

62. Lu CC, Hu ZB, Wang R, Hong ZH, Lu J, Chen PP, et al. Gut microbiota dysbiosis-induced activation of the intrarenal renin-angiotensin system is involved in kidney injuries in rat diabetic nephropathy. Acta Pharmacol Sin. (2020) 41:1111–8. doi: 10.1038/s41401-019-0326-5

PubMed Abstract | Crossref Full Text | Google Scholar

63. Miao H, Wang YN, Su W, Zou L, Zhuang SG, Yu XY, et al. Sirtuin 6 protects against podocyte injury by blocking the renin-angiotensin system by inhibiting the wnt1/beta-catenin pathway. Acta Pharmacol Sin. (2024) 45:137–49. doi: 10.1038/s41401-023-01148-w

PubMed Abstract | Crossref Full Text | Google Scholar

64. Zhou L and Liu Y. Wnt/beta-catenin signaling and renin-angiotensin system in chronic kidney disease. Curr Opin Nephrol Hypertens. (2016) 25:100–6. doi: 10.1097/MNH.0000000000000205

PubMed Abstract | Crossref Full Text | Google Scholar

65. Pluznick JL. Microbial short-chain fatty acids and blood pressure regulation. Curr Hypertens Rep. (2017) 19:25. doi: 10.1007/s11906-017-0722-5

PubMed Abstract | Crossref Full Text | Google Scholar

66. Wu W, Sun M, Chen F, Cao AT, Liu H, Zhao Y, et al. Microbiota metabolite short-chain fatty acid acetate promotes intestinal iga response to microbiota which is mediated by gpr43. Mucosal Immunol. (2017) 10:946–56. doi: 10.1038/mi.2016.114

PubMed Abstract | Crossref Full Text | Google Scholar

67. Yang W, Xiao Y, Huang X, Chen F, Sun M, Bilotta AJ, et al. Microbiota metabolite short-chain fatty acids facilitate mucosal adjuvant activity of cholera toxin through gpr43. J Immunol. (2019) 203:282–92. doi: 10.4049/jimmunol.1801068

PubMed Abstract | Crossref Full Text | Google Scholar

68. Takeuchi T, Miyauchi E, Kanaya T, Kato T, Nakanishi Y, Watanabe T, et al. Acetate differentially regulates iga reactivity to commensal bacteria. Nature. (2021) 595:560–4. doi: 10.1038/s41586-021-03727-5

PubMed Abstract | Crossref Full Text | Google Scholar

69. Zhu Y, He H, Sun W, Wu J, Xiao Y, Peng Y, et al. Iga nephropathy: gut microbiome regulates the production of hypoglycosilated iga1 via the tlr4 signaling pathway. Nephrol Dial Transplant. (2024) 39:1624–41. doi: 10.1093/ndt/gfae052

PubMed Abstract | Crossref Full Text | Google Scholar

70. Li XJ, Shan QY, Wu X, Miao H, and Zhao YY. Gut microbiota regulates oxidative stress and inflammation: A double-edged sword in renal fibrosis. Cell Mol Life Sci. (2024) 81:480. doi: 10.1007/s00018-024-05532-5

PubMed Abstract | Crossref Full Text | Google Scholar

71. Miao H, Wang YN, Yu XY, Zou L, Guo Y, Su W, et al. Lactobacillus species ameliorate membranous nephropathy through inhibiting the aryl hydrocarbon receptor pathway via tryptophan-produced indole metabolites. Br J Pharmacol. (2024) 181:162–79. doi: 10.1111/bph.16219

PubMed Abstract | Crossref Full Text | Google Scholar

72. Zhang Y, Zhou Y, Cui W, Wang Z, Wang X, Wu F, et al. Characterization and diagnostic value of the gut microbial composition in patients with minimal change disease. Front Physiol. (2022) 13:1070569. doi: 10.3389/fphys.2022.1070569

PubMed Abstract | Crossref Full Text | Google Scholar

73. Chu C, Behera TR, Huang Y, Qiu W, Chen J, and Shen Q. Research progress of gut microbiome and diabetic nephropathy. Front Med (Lausanne). (2024) 11:1490314. doi: 10.3389/fmed.2024.1490314

PubMed Abstract | Crossref Full Text | Google Scholar

74. Dong Y, Yan G, Zhang Y, Zhou Y, and Shang J. Gut microbiota as a predictive tool for outcomes in iga nephropathy. Ren Fail. (2025) 47:2514184. doi: 10.1080/0886022X.2025.2514184

PubMed Abstract | Crossref Full Text | Google Scholar

75. Zhi W, Li A, Wang Q, Yuan X, Qing J, Zhang C, et al. Safety and efficacy assessment of fecal microbiota transplantation as an adjunctive treatment for iga nephropathy: an exploratory clinical trial. Sci Rep. (2024) 14:22935. doi: 10.1038/s41598-024-74171-4

PubMed Abstract | Crossref Full Text | Google Scholar

76. Wang S, Lv D, Jiang S, Jiang J, Liang M, Hou F, et al. Quantitative reduction in short-chain fatty acids, especially butyrate, contributes to the progression of chronic kidney disease. Clin Sci (Lond). (2019) 133:1857–70. doi: 10.1042/CS20190171

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: glomerulonephritis, gut microbiome, IgA nephropathy, short-chain fatty acid, spatial transcriptomics

Citation: Koh JH, Park S, Kang M, Park JI, Lee J, Cho H, Kim JE, Nam H, Kim D, Li M, Park S, Moon KC, Kim HJ, Kim YS, Kim DK and Lee H (2026) Altered acetate metabolism and signaling in IgA nephropathy: an integrated gut microbiome and glomerular spatial transcriptome analysis. Front. Immunol. 16:1665585. doi: 10.3389/fimmu.2025.1665585

Received: 14 July 2025; Accepted: 18 December 2025; Revised: 04 November 2025;
Published: 14 January 2026.

Edited by:

Gabriela Barrientos, National Scientific and Technical Research Council (CONICET), Argentina

Reviewed by:

Xiaoyong Yu, Shaanxi Provincial Hospital of Traditional Chinese Medicine, China
Yijun Dong, First Affiliated Hospital of Zhengzhou University, China
Sho Hamaguchi, Juntendo University, Japan

Copyright © 2026 Koh, Park, Kang, Park, Lee, Cho, Kim, Nam, Kim, Li, Park, Moon, Kim, Kim, Kim and Lee. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Hajeong Lee, bWRoamxlZTlAc251LmFjLmty

†These authors share first authorship

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.