PKR activation-induced mitochondrial dysfunction in HIV-transgenic mice with nephropathy

  1. Teruhiko Yoshida  Is a corresponding author
  2. Khun Zaw Latt
  3. Avi Z Rosenberg
  4. Briana A Santo
  5. Komuraiah Myakala
  6. Yu Ishimoto
  7. Yongmei Zhao
  8. Shashi Shrivastav
  9. Bryce A Jones
  10. Xiaoping Yang
  11. Xiaoxin X Wang
  12. Vincent M Tutino
  13. Pinaki Sarder
  14. Moshe Levi
  15. Koji Okamoto
  16. Cheryl A Winkler
  17. Jeffrey B Kopp
  1. Kidney Disease Section, Kidney Diseases Branch, NIDDK, NIH, United States
  2. Department of Pathology, Johns Hopkins Medical Institutions, United States
  3. Department of Pathology and Anatomical Sciences, Jacobs School of Medicine & Biomedical Sciences, University at Buffalo, United States
  4. Department of Biochemistry and Molecular & Cellular Biology, Georgetown University, United States
  5. Polycystic Kidney Disease Section, Kidney Diseases Branch, NIDDK, NIH, United States
  6. Frederick National Laboratory for Cancer Research, NCI, NIH, United States
  7. College of Medicine, University of Florida, United States
  8. Nephrology Endocrinology and Vascular Medicine, Tohoku University Hospital, Japan
8 figures, 1 table and 1 additional file

Figures

PKR inhibition by C16 ameliorates Tg26 mice kidney phenotype.

(A–C) Shown are the following: plasma creatinine (mg/dl, mean with SD), urinary albumin-to-creatinine ratio (mg/g creatinine, mean with SEM), urinary albumin-to-creatinine ratio (mg/g creatinine) of …

Figure 2 with 2 supplements
Overview of bulk RNA-seq and single-nucleus RNA-seq experiments.

(A) Shown is the workflow of the bulk RNA-seq and single-nucleus RNA-seq experiments. (B) Principal component analysis plot of bulk RNA-seq results. (C) Uniform manifold approximation and projection …

Figure 2—figure supplement 1
Pre-processing of single-nuclear RNA-seq data.

(A) Uniform manifold approximation and projection (UMAP) plot of single-nuclear RNA-seq data before doublet removal from 8 samples, 57,061 cells, showing 27 clusters. (B) UMAP plot showing detected …

Figure 2—figure supplement 2
Additional plots of single-nuclear RNA-seq data.

(A) Uniform manifold approximation and projection (UMAP) plot of single-nuclear RNA-seq data from each 8 samples. (B) Feature plot showing Eif2ak2 expression.

Figure 3 with 1 supplement
Oxidative phosphorylation genes are downregulated in Tg26 mice and downregulation is reversed by PKR inhibition using C16.

(A) Shown is the enrichment plot of oxidative phosphorylation pathway based on bulk mRNA-seq comparing Tg26 and wild-type (WT). (B) Shown is the enrichment plot of oxidative phosphorylation pathway …

Figure 3—figure supplement 1
Additional data about mitochondria and S phase score of podocytes.

(A–D) Quantification of western blot probed against mitochondrial subunits normalized by VDAC (mean with SD). (E and F) Relative expression ratio of Nd1 and 16S divided by Hk2 (mean with SEM). (G) …

Figure 4 with 3 supplements
PT-Mito and PT-Inj cluster characterization.

(A) Shown are dot plots showing the top five marker genes in each of the PT clusters (PT-S1, PT-S1/S2, PT-S2, PT-S2/S3, PT-S3, PT-Inj, PT-Mito). (B, C) In situ hybridization of mt-Co1 and mt-Atp6

Figure 4—figure supplement 1
Additional in situ hybridization images.

(A) In situ hybridization images probing dapB (negative control probe) showed no signals. (B) In situ hybridization images probing Ppib (positive control probe) showed strong signals.

Figure 4—figure supplement 2
PT-Mito cluster detection of publicly available human kidney single-nuclear RNA-seq data (GSE131882).

(A) Uniform manifold approximation and projection (UMAP) plot of human kidney single-nuclear RNA-seq data shows 16 clusters. Clusters 1, 4 are proximal tubule (PT) clusters, and cluster 7 is PT-Mito …

Figure 4—figure supplement 3
Single-nuclear RNA-seq data comparison with ischemic reperfusion injury model.

(A) Dot plot showing expression of marker genes reported in new proximal tubular clusters by mouse ischemic reperfusion acute kidney injury model (Kirita.PNAS.2020). (B) Dot plot showing expression …

STAT3 activation downstream of PKR.

(A, B) Shown is mapping of STAT3 regulating genes comparing Tg26 vs WT (Z-score 6.286), Tg26 C16 vs Tg26 (Z-score –3.182) by bulk mRNA-seq. Red color indicates upregulation and green color indicates …

PKR inhibition with C16 reverses mitochondrial dysfunction in Tg26 glomeruli and proximal tubules.

(A) Shown are oxygen consumption rate (OCR) measurements during cell mitochondrial stress testing, using extracted glomeruli from wild-type (WT), Tg26, and C16 Tg26 kidney tissue. (B) Maximum …

HIV-1 gene expression causes Tg26 podocytes dedifferentiation.

(A) Shown is dot plot demonstrating HIV-1 gene expression levels in each cluster detected by snRNA-seq. (B–D) p57 staining of kidney showing podocyte loss and dedifferentiation (scale bars are 50 …

Cell-cell interaction analysis shows activated ligand-receptor interaction: platelet-derived growth factor (PDGF)-D-PDGFR-B pathway in Tg26.

(A–C) Circles plot showing cell-cell interaction weight in each sample at the resolution of each identified cell cluster. (D–F) Shown are dot plots depicting results from cell-cell interaction …

Tables

Table 1
Activation Z-scores from bulk RNA-seq data.

Activation Z-scores of transcription factors by upstream regulator analysis of Ingenuity Pathway Analysis. Triple Z-score was calculated by multiplying three Z-scores in each comparison.

Upstream regulatorActivation Z-scoreWT C16 vs WTActivation Z-scoreTg26 vs WTActivation Z-scoreTg26 C16 vs Tg26Triple Z-score
STAT3–2.1836.286–3.18243.66
FOXO1–3.6183.991–2.69738.94
GLI1–2.3355.517–2.93637.82
EPAS1–2.6653.554–3.24230.71
SREBF1–3.0133.672–2.6329.10
CREB1–2.3254.181–2.91628.35
ATF4–2.6064.445–2.02823.49
SMAD2–2.2072.791–2.33514.38

Additional files

Download links