Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
. 2024 Feb 23;383(6685):eadd6371.
doi: 10.1126/science.add6371. Epub 2024 Feb 23.

SOX9 switch links regeneration to fibrosis at the single-cell level in mammalian kidneys

Affiliations

SOX9 switch links regeneration to fibrosis at the single-cell level in mammalian kidneys

Shikhar Aggarwal et al. Science. .

Abstract

The steps governing healing with or without fibrosis within the same microenvironment are unclear. After acute kidney injury (AKI), injured proximal tubular epithelial cells activate SOX9 for self-restoration. Using a multimodal approach for a head-to-head comparison of injury-induced SOX9 lineages, we identified a dynamic SOX9 switch in repairing epithelia. Lineages that regenerated epithelia silenced SOX9 and healed without fibrosis (SOX9on-off). By contrast, lineages with unrestored apicobasal polarity maintained SOX9 activity in sustained efforts to regenerate, which were identified as a SOX9on-on Cadherin6pos cell state. These reprogrammed cells generated substantial single-cell WNT activity to provoke a fibroproliferative response in adjacent fibroblasts, driving AKI to chronic kidney disease. Transplanted human kidneys displayed similar SOX9/CDH6/WNT2B responses. Thus, we have uncovered a sensor of epithelial repair status, the activity of which determines regeneration with or without fibrosis.

PubMed Disclaimer

Figures

Fig. 1.
Fig. 1.. scRNA-seq reveals a dynamic SOX9/CDH6 switch within the repairing nephron epithelium.
(A) Identification of a dynamic Sox9 switch. Shown are the schema of lineage-tracing of injury-induced Sox9pos cells and co-immunoanalysis for SOX9 and the basolateral polarity marker ATP1A1, showing that cells with disrupted basolateral polarity activated SOX9. Note the co-localization of tdT with SOX9pos cells at 48 hours after IRI (IRI 48 hours: Sox9on cells, yellow arrows). By day 10 after IRI, the lineage of Sox9on cells that restored basolateral polarity had silenced SOX9 (Sox9on-off cell state, regenerated epithelia), whereas the lineage with unrestored ATP1A1 maintained SOX9 activity (Sox9on-on cells, yellow arrows). (B) Schema of isolating single cells to compare IRI 48 hours Sox9pos (n = 2) versus day 10 Sox9pos cells, with dot plot showing average gene expression values and percentage of cells expressing cell-type- specific markers of differentiation, injury (Havcr1 and Lcn2), and repair (Sox9) response. (C) Time-resolved UMAP projection of Sox9pos cells within the proximal tubule cluster. Arrows show distinct clusters at day 10 after IRI containing Sox9pos cells compared with 48 hours after IRI. (D) Topmost enriched GO cellular components and biological processes in IRI day 10 Sox9pos cells versus 48 hour Sox9pos cells. (E) Heatmap of genes driving the GO terms. Shown is the enrichment of epithelial cell-cell adhesion machinery, including Cdh6, within day 10 Sox9pos cells. (F) Time-resolved UMAP projection of Cdh6pos cells (arrows) within the PT cluster. (G and H) Co-immunostaining for CDH6 and SOX9 (G) and ATP1A1 (H) showing CDH6 expression restricted to day 10 SOX9pos cells (G) and to the Sox9 lineage with unrestored ATP1A1-based basolateral polarity versus the lineage that restored polarity (H), indicating that CDH6 demarcates Sox9on-on cells (SOX9posCDH6pos cell state). (I) Co-immunostaining for SOX9 and the epithelial tight junction marker ZO-1 and immunoblot for SOX9 and CDH6 showing that cells that lacked cell-cell contact (subconfluent state) activated SOX9 and CDH6, which waned upon restoration of tight junctions (confluent monolayer). Box inset highlights the region depicted at high magnification. All n = 3 animals per time point unless otherwise stated. Data are shown as mean ± SEM. Scale bars, 100 μm. For total cells counted, see table S1.
Fig. 2.
Fig. 2.. CDH6 status is tightly linked to repair with or without αSMApos myofibroblast response through single-cell Wnt activity.
(A) Schema of lineage-tracing of IRI-induced Sox9pos cells and co-immunostaining showing αSMApos myofibroblasts adjacent to SOX9postdTpos cells (Sox9on-on cells, arrows) but no detectable αSMA activity around the Sox9on-off cells (arrow). Left and right magnified panels correspond to the foci highlighted by corresponding left and right box insets. (B) Co-immunoanalysis showing αSMApos myofibroblasts encasing the CDH6pos Sox9 lineages at single-cell spatial distance; by contrast, the paucity of such response around the CDH6negSox9-lineage can be observed (arrow). Circle and box highlight such foci, with the magnified panel corresponding to the foci highlighted by the box inset. (C) Schemas for purifying Slc34a1pos normal PTECs (control tdTpos cells) and Sox9pos cells 48 hours and day 14 after IRI and principal components analysis (PCA) plot showing distinct transcriptomic profiles of purified tdTpos cells. Sox9 transcripts were enriched in purified tdTpos cells. (D) qPCR of purified tdTpos cells confirming Sox9 and Cdh6 enrichment. (E) UMAP projections for Sox11 at 48 hours and day 10 after IRI. Arrows highlight the distinct day 10 cluster composed of these cells. Single Sox9pos and Cdh6pos cells also contributed to this cluster (Fig. 1, B and E). (F) Volcano plot demonstrating enrichment of Wnt ligands within day 14 versus 48 hour tdTpos cells. (G) UMAP projections for Wnt4 and Wnt7b at 48 hours and day 10 after IRI. Arrows highlight the distinct day 10 cluster, which contained such single Wnt4pos and Wnt7bpos cells. Note that the same distinct cluster was composed of day 10 Sox9pos, Cdh6pos, and Sox11pos cells [see also Fig. 1, B and E, and this figure (E)]. (H) Upset plot analysis of scRNAseq datasets showing subset of single Cdh6pos cells enriched with Wnt4, including Cdh6posSox11pos subsets. (I) qPCR of purified tdTpos cells confirming Wnt4 enrichment in day 14 post-IRI tdTpos cells versus their 48-hour counterparts. (J) Schema for labeling early Axin2pos cells after IRI and immunoanalysis for CDH6 showing co-localization of tdTpos and GFPpos cells with tdTposGFPpos cells located adjacent to CDH6pos subsets. Circles and box highlight single-cell, biologically active Wnt-enriched niches robustly linked with subsets of CDH6pos cells. Magnified panels correspond to the foci highlighted by box inset. All n = 3 animals per time point unless otherwise stated. Data are shown as mean ± SEM. Scale bars, 100 μm. For total cells counted, see table S1.
Fig. 3.
Fig. 3.. Sox9on-on cell state is the main driver of myofibroblast formation and maintenance.
(A to C) Schema of Wls removal from Sox9on-on cells, outer medullary region representative image, and blinded co-immunoanalysis for SOX9 and αSMA (n = 5 animals/group). No difference can be seen in Sox9on-on cells between the two groups (B). The mosaic tissue–damaged microenvironment displayed significantly (P < 0.01) reduced αSMApos myofibroblasts adjacent to tdTpos demarcated Wls-KO tubules (arrows) versus tdTneg tubules with intact Wls [(C), arrowheads]. (D) Schema of isolating early Axin2pos cells after IRI with FACS plot showing that nearly all tdTpos cells co-expressed GFP. (E) PCA plot of purified tdTposGFPpos and tdTnegGFPpos cells. (F) GO analysis showing Wnt pathway among the top 10 terms, confirming that Axin2pos cells demarcated early WRCs. (G and H) Volcano plot showing the molecular signature of early WRCs, with qPCR confirming Axin2 enrichment within tdTposGFPpos cells, thus validating reporter animals. (I) UMAP projection of Axin2, Nkd1, and Nkd2 cells. Such cells were within the same PdgfrbposCola1posCol3a1pos cluster (see also fig. S9). (J) Experimental outline of removal of Sox9on-on activity during the AKI to CKD transition. (K and L) Trichrome staining (K) and quantitative scores of interstitial fibrosis (L) (blinded analysis, n = 5 Sox9-Ctrl and 6 Sox9-cKO animals). (M and N) Co-immunostaining for SOX9, KIM1, and αSMA (M) and blinded co-immunoanalysis showing head-head comparison between KIM1+ regions for αSMA activity (N) in Sox9-cKO versus Sox9-Ctrl animals. (O) qPCR analysis of genes associated with fibrosis in the kidneys from Sox9-cKO versus Sox9-Ctrl animals. (P and Q) Western blot confirming SOX9 knock-down in subconfluent primary TECs (P), with qPCR showing Sox9 and Wnt4 down-regulation (Q). (R and S) RNAscope image (R) and analysis (S) showing significant (P < 0.01) Wnt4 reduction. (T) qPCR analysis of Wnts under scrutiny and Axin2, Nkd1, and Nkd2 in the kidneys from Sox9-cKO versus Sox9-Ctrl. All images are representative images. Data are shown as mean ± SEM. *P < 0.05, **P < 0.01 unpaired two-sided Student’s t test; ***P < 0.01, paired Student’s t test. Scale bars in whole scanned images [(K) and (R)], 1000 μm; all others, 100 μm. For total cells counted, see table S1.
Fig. 4.
Fig. 4.. Epigenetic reprogramming of single Sox9on-on nuclei to a nephron progenitor–like state contrasts with the Sox9on-off counterpart, which reverts to normal PTEC.
(A) Schema of the workflow to obtain single nuclei of the FACS-enriched tdTpos Sox9-descendants at day 10 after injury. Rectangle (inset) highlights the dissected region showing the inner cortices and outer medulla, the site of relatively extensive PTEC loss. (B) Integration of scRNA-seq and snATAC-seq showing the cellular identity of single-tdTpos nuclei (n = 2 animals). (C) Identification of Sox9on-on and Sox9on-off single nuclei: Chromatin accessibility revealing clustering of tdTpos descendants based on relatively open and closed chromatin accessibility state of Sox9. Sox9on-off nuclei exhibited open chromatin accessibility for Hnf4a, a known marker of healthy, mature PTECs, suggesting that the Sox9-lineage that regenerated normal PTECs closed chromatin accessibility for Sox9. (D) UMAP representation of chromatin accessibility analysis showing the relatively open chromatin state of Havcr1, Wnt7b, Wnt4, and Cdh6 versus their accessibility state in Sox9on-off nuclei. (E and F) Sox9 trajectory analysis with scRNAseq-imputed gene expression (E) and pseudotime (F) scale showing the Sox9on-on → Sox9on-off transition. (G) Heatmaps of cross-platform linked genes involved in the transcriptional cascade during Sox9on-on → Sox9on-off transition. Note that the highlighted genes in Sox9on-on nuclei are linked with tissue development and/or nephrogenesis, and the Sox9on-off nuclei displayed open chromatin accessibility for Hnf4β, another known marker of healthy, mature PTECs in addition to Hnf4a. (H) Schema of EDU regime and co-immunostaining for SOX9 and EDU showing SOX9postdTposEDUpos cells. Two representative images are shown, with the lower panel demonstrating a cluster of such cells. n = 3 animals. (I and J) Co-immunostaining for SMARCC1 and CDH6 showing SMARCC1 expression confined to day 10 CDH6pos Sox9 lineage cells (I), and RUNX1 and SOX9 co-immunostaining showing RUNX1 restricted to Sox9on-on versus Sox9on-off lineage (J) as predicted by (G), thus validating Sox9 lineage-specific snATAC-seq datasets. n = 2 animals. Scale bars, 100 μm.
Fig. 5.
Fig. 5.. Human renal allografts display dynamic SOX9/CDH6/WNT2B activity, with CDH6pos cells demarcating SOX9on-on activity and fibrotic foci.
(A) Co-immunostaining for SOX9 and LTL in biopsies obtained before implantation (uninjured) and IRI-induced AKI after transplantation (same allograft) showing early SOX9 activation. (B) Box plot showing SOX9 levels at different time points within kidney transplant protocol biopsies (n = 163). (C and D) Dot plots (C) and box plots (D) showing categorization of patients and kidney function according to SOX9 levels at 1 year (comparison by Mann-Whitney U test; n = 35). (E) Histograms showing the number of patients with different degrees of kidney fibrosis, as estimated by ci-score according to Banff classification (comparison by chi-square test; n = 39). (F) Box plot showing SOX9 and CDH6 levels categorized according to a previously reported model discriminating the transcriptome of kidney transplant biopsies in successful repair (1), transition to chronic injury (2), or CKD (3). LD, biopsies obtained from living donors at the time of transplantation (49). (G) Heatmap showing expression correlation of genes of interest at 3 and 12 months after transplantation (Spearman r correlation coefficient; n = 72). Arrows highlight the identified molecular signatures of SOX9/CDH6/WNT, including AXIN2 and NKD1 in human kidneys. Note that CDH6 is the topmost correlated gene. (H) UMAP showing distinct PTEC clustering (29,180 genes × 24,070 cells). (I) Dot plot showing average gene expression values and percentage of cells expressing markers of differentiation and injury, SOX9 and CDH6, by each identified cluster of PTECs. Cluster 8 consisted of SOX9-, CDH6-, and WNT2B-expressing cells (see also figs. S22, C to E, and S23). (J) Feature plot displaying the normalized transcript expression for the respective genes. Magic mRNAassay of renal cells is depicted. (K) Analysis comparing the percentage of cluster 8 and the interstitial fibrosis and tubular atrophy (IFTA) grade for each patient. Percentage of cluster 8 is calculated with respect to the PTEC number in each patient. P < 0.05 based on the Kruskal-Wallis test and Pearson correlation analysis. (L) UMAP showing time-resolved PTECs clustering after transplantation (shown as days after transplantation). Note the emergence of cluster 8 with time after transplantation. (M) Cluster-type-specific genes analysis revealed CDH6 as being among the top two driver genes that underlie the dynamic activity in cluster 8. (N) Co-immunostaining showing that CDH6pos cells displayed a tight intimate association with ACTA2pos myofibroblasts within human kidney allograft, with CDH6neg foci showing no ACTA2pos myofibroblasts. (O) Immunoblot showing that subconfluent human primary PTECs activated SOX9, which waned upon confluency. (P and Q) Co-immunoanalysis (P) and qPCR (Q) showing reduction in CDH6pos cells and Cdh6 mRNA, respectively, upon removal of Sox9on-on activity (blinded analysis, unpaired two-sided Student’s t test; data are shown as mean ± SEM). Scale bars in whole scanned image (P), 1000 μm; all others, 100 μm.

Comment in

Similar articles

Cited by

References

    1. Rockey DC, Bell PD, Hill JA, Fibrosis—A common pathway to organ injury and failure. N. Engl. J. Med 372, 1138–1149 (2015). doi: 10.1056/NEJMra1300575; pmid: - DOI - PubMed
    1. Henderson NC, Rieder F, Wynn TA, Fibrosis: From mechanisms to medicines. Nature 587, 555–566 (2020). doi: 10.1038/s41586-020-2938-9 - DOI - PMC - PubMed
    1. Konieczny P, Naik S, Healing without scarring. Science 372, 346–347 (2021). doi: 10.1126/science.abi5770 - DOI - PubMed
    1. Kumar S, Cellular and molecular pathways of renal repair after acute kidney injury. Kidney Int 93, 27–40 (2018). doi: 10.1016/j.kint.2017.07.030; - DOI - PubMed
    1. Liu J et al., Molecular characterization of the transition from acute to chronic kidney injury following ischemia/reperfusion. JCI Insight 2, e94716 (2017). doi: 10.1172/jci.insight.94716; - DOI - PMC - PubMed

Substances