Abstract
The disproportionate risk for idiopathic proteinuric podocytopathies in Black people is explained, in part, by the presence of two risk alleles (G1 or G2) in the APOL1 gene. The pathogenic mechanisms responsible for this genetic association remain incompletely understood. We analyzed glomerular RNASeq transcriptomes from patients with idiopathic nephrotic syndrome of which 72 had inferred African ancestry (AA) and 152 did not (noAA). Using gene coexpression networks we found a significant association between APOL1 risk allele number and the coexpression metamodule 2 (MM2), even after adjustment for eGFR and proteinuria at biopsy. Unadjusted Kaplan-Meier curves showed that unlike noAA, AA with the highest tertile of MM2 gene activation scores were less likely to achieve complete remission (p≤0.014). Characteristic direction (ChDir) identified a signature of 1481 genes, which separated patients with APOL1 risk alleles from those homozygous for reference APOL1 . Only in AA, the tertile with the highest activation scores of these 1481 genes was less likely to achieve complete remission (p≤0.022) and showed a trend to faster progression to the composite event of kidney failure or loss of 40% eGFR (p≤0.099). The MM2 and ChDir genes significantly overlapped and were both enriched for Epithelial Mesenchymal Transition and inflammation terms. Finally, MM2 significantly overlapped with a parietal epithelial cell (PEC)-identity gene signature but not with a podocyte identity signature. Podocytes expressing variant APOL1s may generate inflammatory signals that activate PECs by paracrine mechanisms contributing to APOL1 nephropathy.