Generation of Monogenic Candidate Genes for Human Nephrotic Syndrome Using 3 Independent Approaches.

INTRODUCTION: Steroid-resistant nephrotic syndrome (SRNS) is the second most common cause of chronic kidney disease during childhood. Identification of 63 monogenic human genes has delineated 12 distinct pathogenic pathways.
METHODS: Here, we generated 2 independent sets of nephrotic syndrome (NS) candidate genes to augment the discovery of additional monogenic causes based on whole-exome sequencing (WES) data from 1382 families with NS.
RESULTS: We first identified 63 known monogenic causes of NS in mice from public databases and scientific publications, and 12 of these genes overlapped with the 63 known human monogenic SRNS genes. Second, we used a set of 64 genes that are regulated by the transcription factor Wilms tumor 1 (WT1), which causes SRNS if mutated. Thirteen of these WT1-regulated genes overlapped with human or murine NS genes. Finally, we overlapped these lists of murine and WT1 candidate genes with our list of 120 candidate genes generated from WES in 1382 NS families, to identify novel candidate genes for monogenic human SRNS. Using this approach, we identified 7 overlapping genes, of which 3 genes were shared by all datasets, including SYNPO. We show that loss-of-function of SYNPO leads to decreased CDC42 activity and reduced podocyte migration rate, both of which are rescued by overexpression of wild-type complementary DNA (cDNA), but not by cDNA representing the patient mutation.
CONCLUSION: Thus, we identified 3 novel candidate genes for human SRNS using 3 independent, nonoverlapping hypotheses, and generated functional evidence for SYNPO as a novel potential monogenic cause of NS.

NS is a kidney disease defined by proteinuria with resulting hypoalbuminemia, frequently causing edema and hyperlipidemia., In the clinical setting, “steroid-sensitive” versus “steroid-resistant” NS (SSNS vs. SRNS) can be distinguished based on the patient’s response to standard steroid therapy. SRNS with the histological correlate of focal segmental glomerulosclerosis (FSGS) invariably leads to end-stage renal failure. SRNS is the second most common cause of chronic kidney disease in childhood, with a reported incidence of 1.15 to 16.9 per 100,000 children., Over the past years, 63 genes have been identified as causing monogenic forms of SRNS (recessive or dominant) in humans with an onset <25 years.6, 7, 8, 9 Characterization of the cellular function of these 63 genes helped delineate 12 distinct pathogenic pathways of SRNS/FSGS.6, 7, 8, 9, 10, 11 As these genes are mainly expressed in podocytes, gene identification has helped to reveal a central role of the podocyte in the pathogenesis of SRNS.6, 7, 8, 9, 10 Monogenic causation accounts for 11.0% to 29.5% of patients with SRNS with onset before 25 years,,,, leaving up to 70% of cases in that age group genetically unsolved. Thus, we surmised that novel genes still remain to be discovered.

Most SRNS genes identified in recent years were discovered by WES. Thus, WES, as an unbiased approach, still represents the most efficient method to identify potential novel causes of SRNS in humans. However, as childhood-onset SRNS represents a rare disease with a limited number of patients worldwide, and as WES analysis often reveals multiple variants, an independent or “orthogonal” dataset of candidate genes may facilitate the WES filtering process. Here, we generated 2 datasets of candidate genes for SRNS. After validating both sets against a list of known human NS genes, we overlapped these with our list of unique candidate genes from WES in 1382 families with NS. Specifically, we generated the following lists of known or candidate genes: (i) 63 known human NS genes were used for validation, (ii) 63 known monogenic causes of murine NS were retrieved from the literature, (iii) 64 genes regulated by the transcription factor WT1 were used.,

By validating the mouse and WT1-derived candidate datasets, we found that 12 of the murine NS genes and 5 of the WT1 downstream targets overlapped with the validation set of known human NS genes. We then overlapped the mouse and WT1-derived candidate datasets with our unique list of 120 NS candidate genes (iv).,, This approach allowed us to identify 3 potential novel candidate genes (ITGB8, SEMA3G, Buerger, 2020, personal communication, functional studies in progress, synaptopodin [SYNPO]) for human NS.

Within this set of 3 novel candidate genes, we identified a homozygous mutation in the gene SYNPO in a family with steroid-dependent NS. We show that loss-of-function of SYNPO leads to decreased CDC42 activity and reduced podocyte migration rate, which were each rescued by overexpression of wild-type cDNA, but not by cDNA representing the patient mutation.

Thus, we identified 3 novel candidate genes for human SRNS using 3 independent nonoverlapping hypotheses, generating functional evidence for SYNPO as a potential novel candidate gene for human NS.

Methods

Research Subjects

This study was approved by the institutional review board of Boston Children’s Hospital. We obtained blood samples and pedigrees following informed consent from families with NS. The diagnoses of NS were based on published clinical criteria and renal biopsies criteria evaluated by renal pathologists. Clinical data were obtained using a standardized questionnaire (http://www.renalgenes.org).

WES and Mutation Calling

WES was performed in 1382 NS families from worldwide sources, as previously described, using Agilent SureSelect or Illumina Nextera human exome capture arrays (Thermo Fisher Scientific, Waltham, MA) with next generation sequencing on an Illumina (San Diego, CA) platform. In total we examined 1382 families (1624 individuals) with at least 1 affected individual with NS (1198 families had only 1 affected individual, 142 families had 2 affected individuals, 33 families had 3 affected individuals, 6 families had 4 affected individuals, 1 family had each 5 or 6 affected individuals, and 1 family had 8 affected individuals; families with multiple affected individuals mostly represented sibling cases, some of them also included affected parents, uncles, aunts, or cousins). Trio/duo evaluation was performed in case DNA of unaffected parents was available (trios: 320, duos [maternal or paternal DNA available]: 138). Sequence reads were mapped against the human reference genome (NCBI build 37/hg19) using CLC Genomics Workbench (version 6.5.1) (CLC bio, Aarhus, Denmark). Genetic location information is according to the February 2009 Human Genome Browser data, hg19 assembly (http://www.genome.ucsc.edu). Downstream processing of aligned BAM files were done using Picard and samtools60, and SNV calling was done using Genome Analysis Tool Kit 5. Mutation calling was performed in line with proposed guidelines, and the criteria were used as previously described.,,19, 20, 21, 22, 23, 24 Briefly, after alignment to the human reference genome, variants were filtered for most likely deleterious variants. The variants included were rare in the population with mean allele frequency <1% in dbSNP147 and with only 0–1 homozygotes in the adult genome database gnomAD. Synonymous and intronic variants that were not located within splice-site regions were excluded.

First, WES data were analyzed for 1 of the 63 known NS genes. In a subcohort of families we have identified a known monogenic NS gene in up to 29.5% of the cases; this information is not included in this study but is published elsewhere.,,, Unsolved cases were then subjected to an unbiased WES analysis. Subsequently, variant severity was stratified based on protein impact (truncating frameshift or nonsense mutations, essential or extended splice-site mutations, and missense mutations). Splice-site mutations were assessed by in silico tools MaxEnt and NNSPLICE splice-site mutation prediction scores,,, as well as conservation across human splice sites as described previously. Missense mutations were assessed based on SIFT, MutationTaster, and PolyPhen 2.0 conservation prediction scores22, 23, 24 and evolutionary conservation based on manually derived multiple sequence alignments.

Homozygosity Mapping

Homozygosity mapping was performed based on WES data. In brief, aligned BAM files were processed using Picard and SAMtools as described by other groups. Single nucleotide variant calling was performed using Genome Analysis Tool Kit. The resulting VCF files were used to generate homozygosity mapping data and visual outputs using the program Homozygosity Mapper. In case an individual had a relevant region of homozygosity (>60 MBp), we prioritized regions of homozygosity mapping in our evaluation of WES by initially analyzing variants identified in regions of homozygosity as published recently. Approximately 35% of the families, in which we identified 1 of the 120 candidate genes, had relevant homozygosity (>60 MBp). The 120 candidate genes each represented the strongest variants identified by this unbiased approach.

Generating a Candidate Gene Set of Monogenic Mouse Models of NS

Stringent criteria were used to generate a set of known monogenic mouse models of NS (Supplementary Figure S1) Two approaches were combined: (i) An “MGI Human-Mouse” Disease Connection search (http://www.informatics.jax.org) was performed by filtering for the terms: “nephrotic syndrome,” “proteinuria,” “foot process effacement,” or “glomerulonephritis.” Mouse models with nephritis rather than NS, hypertension models or autoimmune disease–related causes were excluded. (ii) A PubMed search (https://www.ncbi.nlm.nih.gov/pubmed/) was performed by searching for the following terms: “mouse model” AND “nephrotic”, “mice” AND “nephrotic”, “mice” AND “foot process effacement”, “mice” AND “proteinuria” AND “gene”. The same exclusion criteria were applied. In this way, 111 monogenic mouse models of nephrotic syndrome/proteinuria were identified (Supplementary Figure S1). Of these, only 75 models were based on global knockout, 36 tissue-specific mouse models were excluded, as they did not generate a hypothesis of a monogenic NS phenotype. Twelve of the 75 mouse NS genes were excluded from the study, as the murine NS phenotype was reported after the human NS gene was published in the literature, leaving us with 63 genes (Supplementary Tables S1 S2, and S3).

cDNA Cloning

Human SYNPO cDNA clones were purchased (Genscript, Piscataway, NJ; OHu26207). Human SYNPO full-length protein (GenBank accession NP_ 001159680) is encoded by GenBank accession NM_001166208. Mutagenesis was performed using the QuikChange II XL Site-Directed Mutagenesis Kit (Agilent Technologies, Santa Clara, CA). Expression constructs (pRK5-N-Myc, pCDNA6.2-N-GFP) were produced using LR Clonase (Invitrogen, Thermo Fisher Scientific) following the manufacturer’s instructions. Clones were validated prior to usage (Supplementary Figures S2 and S3).

Cell Lines

Experiments were done in HEK293T cells (ATCC, Manassas, VA) or immortalized human podocytes (gift from Moin Saleem, University of Bristol, Bristol, UK). Podocytes were cultured as previously described.

Short Hairpin RNA Transduction

Short hairpin RNA (shRNA) targeting human SYNPO was subcloned into pSIREN RetroQ for retroviral transduction using HEK293T cells 48 hours after transduction puromycin was added to the medium for selection of transduced cells at a final concentration of 2.4 μg/ml.

Reagents and Antibodies

The following primary antibodies were used: mouse monoclonal anti-SYNPO (G1D4; PROGEN, Heidelberg, Germany), rabbit anti-SYNPO (PA5–56997; Thermo Fisher Scientific), mouse horseradish peroxidase–linked anti-beta actin (ab20272; Abcam, Cambridge, UK). Donkey anti-mouse, anti-rabbit Alexa 488– and Alexa 594–conjugated secondary antibodies; 4′,6-diamidino-2-phenylindole (DAPI) staining reagents; Phalloidin-Alexa 594 were obtained from Invitrogen (Thermo Fisher Scientific). Horseradish peroxidase–labeled secondary antibodies were purchased from Santa Cruz Biotechnology (Dallas, TX).

Immunofluorescence and Confocal Laser Scanning Microscopy

For overexpression studies, human podocytes were transiently transfected using Lipofectamine 2000 (Thermo Fisher Scientific) according to the manufacturer’s instructions. Experiments were performed 24 to 48 hours after transfection. Cells were fixed for 15 minutes using 4% paraformaldehyde and permeabilized with 0.5% Triton-X 100. After blocking, sections were incubated overnight at 4°C with primary antibody. The cells were incubated in secondary antibodies for 60 minutes at room temperature followed by mounting in hardening medium with 4′,6-diamidino-2-phenylindole (DAPI). Confocal imaging was performed using the Leica (Wetzlar, Germany) SP5X system with an upright DM6000 microscope, images were processed with the Leica AF software suite.

G-LISA Activation Assays

Cells were transfected with wild-type or Mock constructs using Lipofectamine-2000. Transfected cells were incubated in RPMI medium with 10% fetal bovine serum for 24 hours and then in serum-free medium for 24 hours. RAC1 or CDC42 activity was determined using a G-LISA Activation Assay Biochem Kit (Cytoskeleton, Denver, CO), according to the manufacturer’s instructions.

Podocyte Migration Assay

Podocyte migration assay was performed using Incucyte ZOOM System (Essen Bioscience, Ann Arbor, MI) according to the manufacturer’s instructions. SEM is presented.

Statistical Analysis

Results are presented as SEM or mean ± SD for the indicated number of experiments. Statistical analysis was performed with 1-way analysis of variance using GraphPad (La Jolla, CA) Prism. P < 0.05 was considered statistically significant.

Results

Three Orthogonal Approaches to Generate Human NS Candidate Genes

To identify novel candidate genes for human monogenic NS, we generated 3 candidate gene lists that were considered “orthogonal” (i.e., independent) with regard to their biological and functional candidate hypothesis: (i) the first approach comprised a list of 63 published monogenic mouse models of NS or proteinuria (Figure 1, part B, Supplementary Figure S1, Supplementary Tables S1 and S2). Stringent search criteria were applied and 63 genes were identified (Supplementary Figure S1). To generate this list, 36 proteinuric mouse models based on conditional gene knockout were excluded. Furthermore, 12 mouse models with an equivalent human phenotype that were generated after the human NS gene was found were excluded (Supplementary Table S3). (ii) Second, we used a published list of 64 candidate genes based on the hypothesis that podocytic genes, downstream of the master regulator WT1, are potential causes of NS (Figure 1, part C, Supplementary Figure S4). (iii) Finally, we overlapped both candidate gene sets (mouse and WT1) with a list of 120 candidate genes that we generated by unbiased WES evaluation in a cohort of 1,382 families with NS (Figure 1, part D).

Figure 1

Venn diagram of 3 independent approaches to generate candidate genes for monogenic nephrotic syndrome (NS) in humans. (a) Overview of 63 known monogenic causes of NS in humans (blue oval) that serve as a positive control. Twelve of the 63 genes (19%) overlap with 63 genes of known mouse model for NS/proteinuria (green oval) and 5 of the 63 genes (8%) overlap with 64 Wilms Tumor 1 (WT1) downstream targets., (b) Overview of 63 known monogenic mouse models of NS/proteinuria (green oval) (Supplementary Table S1, Supplementary Figure S1). Twelve of the 63 mouse NS genes (19%) overlap with 63 known monogenic causes of human NS. Ten of the 63 mouse NS genes (16%) overlap with WT1 downstream targets (yellow oval), and 6 of the 63 mouse NS genes (10%) genes overlap with unique NS candidate genes resulting from unbiased whole exome sequencing (WES) (red oval). (c) Overview of 64 NS candidate genes (yellow oval) based on WT1 chromatin immunoprecipitation sequencing (ChIP-Seq) data on E18.5 mice published by Lefebvre et al. (Supplementary Figure S4). Five of the 64 genes (8%, KANK1, NPHS1, NPHS2, MYO1E, PTPRO) overlap with 63 known monogenic causes of NS in humans that served as a positive control. Ten of 64 genes (15.6%, MYO1E, PTPRO, CLIC5, CRIM1, MAFB, PODXL, CD55, ITGB8, SEMA3G, SYNPO) overlap with known monogenic mouse models of NS (green oval). Four of 64 genes (7.8%, ITGB8, SEMA3G, SYNPO, DAAM2) overlap with 120 unique NS candidate genes generated by us from unbiased WES (red oval). (d) Overview of 120 unique NS candidate genes generated from unbiased WES in 1382 NS families from a worldwide cohort (red oval). Four of the 120 genes (3.3%, ITGB8, SEMA3G, SYNPO, DAAM2) overlap with 64 WT1 downstream genes (yellow oval) and 6 of the 120 genes (5%, ITGB8, SEMA3G, SYNPO, NOS1AP, PIK3C2A, SEMA3A) overlap with mouse NS genes (green). Three of the 120 genes (2.5%, ITGB8, SEMA3G, SYNPO) overlap with both WT1 and mouse NS candidate genes. Bold font indicates that a family with a potentially causative mutation in this gene was detected. Double frame boxes indicate the 7 strongest candidate genes. ∗Buerger, personal communication, 2020; ∗∗Schneider; ∗∗∗Majmundar and Buerger.

Murine NS Genes

To test the validity of the mouse list, we examined it for overlap with the published 63 known human NS genes (Figure 1, parts A and B, blue: human NS, green: mouse NS).6, 7, 8, 9, 10 Twelve of the 63 NS mouse genes (19%) overlapped with the 63 known human NS genes (Figure 1, Supplementary Table S2) (MYO1E, PTPRO, ARHGDIA, CD2AP, ITGA3, ITSN2, KIRREL1, LAMA5, LAMB2, MAGI2, SCARB2, TNS2). This finding validated the mouse candidate list as suitable for the identification of novel potential monogenic human NS genes.

WT1 Downstream Targets in Podocytes Based on Chromatin Immunoprecipitation Sequencing in E18.5 Mouse Kidneys

WT1 is a master regulator for podocyte development and maintenance and is highly expressed in mature podocytes. Mutations in WT1 have been found to lead to different glomerular diseases, including isolated NS,, Denys-Drash Syndrome, or Frasier Syndrome., WT1 encodes a tumor suppressor protein and putative splicing cofactor, based on a zinc finger structure that has multiple key functions in kidney development, although exact mechanisms remain elusive.40, 41, 42, 43 Recently, different WT1-related studies using chromatin immunoprecipitation sequencing, exon array, and bioinformatical analyses have shown a specific role for WT1 in regulating the podocyte-specific transcriptome.,,, We therefore hypothesized that genes downstream of WT1 represent promising NS candidate genes and generated a second candidate list for this study (Figure 1, part C, WT1 genes: yellow Supplementary Figure S4) based on the following publications: (i) Motamedi et al. performed in vivo chromatin immunoprecipitation sequencing analyses on E18.5 kidneys. In this work, 36,512 WT1-associated regions were identified, which were reduced to 5547 or 2940 peaks, by setting an irreproducibility discovery rate cutoff of ≤0.1 or 0.01, respectively (Supplementary Figure S4). (ii) To identify podocyte-specific genes that are potentially regulated by WT1, Lefebvre et al. then evaluated these WT1 chromatin immunoprecipitation sequencing–associated genes, with a list of genes that represented the “Top200” genes highly expressed in podocytes (Supplementary Figure S4)., Sixty-four of the 192 podocytic genes were associated with WT1-bound regions (irreproducibility discovery rate ≤ 0.01) and harbored a WT1 motif, within 2.5 kb of the genes’ transcriptional start site.,

Thus, 64 podocytic genes were found to be potentially regulated by WT1 according to the criteria employed by Lefebvre et al. We therefore tested the hypothesis that some of these genes represent potential WT1-derived candidates for novel human NS genes (Figure 1, part C, blue: human NS, yellow: WT1 genes). To validate this candidate gene approach “a posteriori,” we overlapped the WT1 genes with human NS genes. Interestingly, 5 of the 64 WT1-derived genes (8%) overlapped with the 63 known human NS genes (Figure 1) (KANK1, NPHS1, NPHS2, MYO1E, PTPRO), of which 2 genes (MYO1E, PTPRO) overlapped with both the human and murine NS lists (Figure 1). This overlap demonstrated that the known human NS genes are enriched 26-fold in the WT1-derived list compared to the approximately 20,000 genes of the entire human genome (Supplementary Figure S5A). After validating the WT1 list against the human NS genes, we then compared the WT1 genes with the candidate list of murine monogenic NS genes (Figure 1). Here we found that 10 of the 64 WT1 candidate genes (15.6%) overlapped with the dataset of 63 mouse NS genes (MYO1E, PTPRO, CLIC5, CRIM1, MAFB, PODXL, CD55, ITGB8, SEMA3G, SYNPO). Hence, monogenic murine NS genes were enriched 50-fold within the WT1 candidate dataset compared to all genes of the entire human genome (Supplementary Figure S5B). We therefore concluded that due to this enrichment, the remaining 59 of the 64 WT1-derived genes, that have not been reported to be a human NS gene, may represent novel candidate genes for monogenic NS in humans.

Unique NS Candidate Genes Resulting From Unbiased WES in 1382 NS Families

Over the past 12 years, our laboratory has performed WES analyses in 1382 families in a worldwide cohort with NS. In a subcohort of families, we have identified a known monogenic NS gene in up to 29.5% of cases.,,, Cases unsolved for known NS genes were then subjected to unbiased WES analysis. Using this approach, we have discovered, functionally characterized, and published 25 novel NS genes.47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60 Mutation calling was performed in line with proposed unbiased guidelines. In addition to identifying known NS genes, we generated potential novel candidate genes for human monogenic NS in 120 families in which only one strong candidate gene resulted from WES evaluation (Figure 1d, red: WES candidates; also see Figure 1 of Waireko et al.). The following criteria were used,: Only rare, biallelic (recessive) variants with a mean allele frequency <1% and 0–1 homozygote allele carriers in the control genome/exome database gnomAD were included. Nonsynonymous variants and/or variants located within splice sites were analyzed. Subsequently, variant severity was stratified by the following: splice-site mutations were assessed by in silico prediction scores.,, Missense mutations were assessed based on SIFT, MutationTaster, and PolyPhen 2.0 prediction scores,22, 23, 24 and evolutionary conservation was evaluated based on manually derived multiple sequence alignments.

Having validated the mouse (Figure 1b, green) and WT1 (Figure 1c, yellow) candidate gene lists against the human (Figure 1a, blue) NS known genes, we then hypothesized that any genes from the mouse or WT1 candidate list that overlap with a WES-derived candidate gene may generate a strong candidate of our 120 WES candidates. We found that 6 of the WES genes overlapped with the mouse candidate list (ITGB8, 1 family; SEMA3G, 3 families, Buerger, 2020, personal communication, functional studies in progress; SYNPO, 1 family; NOS1AP, 2 families, Majmundar and Buerger; PIK3C2A, 2 families; SEMA3A, 1 family), whereas 4 genes overlapped with the WT1-derived candidate dataset (ITGB8, 1 family; SEMA3G, 3 families, Buerger, 2020, personal communication, functional studies in progress; SYNPO, DAAM2, 4 families, Schneider). Interestingly, 3 genes were included in all 3 of the candidate gene lists (ITGB8, 1 family; SEMA3G, 3 families, Buerger, 2020, personal communication, functional studies in progress; SYNPO, 1 family) (Figure 1, Supplementary Table S4).

Overlap Between Candidate Gene Sets

Seven genes from our WES-derived candidate dataset overlapped with either the mouse or the WT1 candidate gene lists (Figure 1, Supplementary Table S5). Interestingly, for 3 of these genes, we have identified different families with multiple, likely causative mutations (NOS1AP, Majmundar and Buerger; DAAM2, Schneider; SEMA3G, Buerger, 2020, personal communication, functional studies in progress) (Figure 1). Furthermore, we discovered families with NS with mutations in the following genes: PIK3C2A, SEMA3A, and ITGB8 (Figure 1, Supplementary Table S5).

In addition, we identified a consanguineous family (FA) harboring a homozygous missense mutation in the gene SYNPO (c.2540C>T, p.P847L) (Figure 2a, Supplementary Table S4). The patient was from Saudi Arabian ethnicity, consanguineous descent, and had an onset of NS at the age of 4 years (Supplementary Table S4). The here-identified variant was absent from the control database gnomAD, and the mutation yielded strong in silico prediction scores (Supplementary Table S4). The proline at position 847 has been evolutionarily conserved since Danio rerio (Figure 2a). The SYNPO mutation was present homozygously in the affected patient FA-24, whereas the healthy parents and siblings were heterozygous carriers of the mutation (Supplementary Figure S6). The patient was reported to initially respond to glucocorticoid treatment but remained steroid-dependent. No renal biopsy was performed. The patient continued to have proteinuria, with slow progression of renal disease despite ongoing steroid therapy. Extrarenal symptoms included muscular dystrophy caused by a homozygous LAMA2 mutation (c.2096G>T, p.R699M). The patient’s brother also suffered from this disease, but not from NS. Hypotonia developed at the age of 40 days. The patient has dysmorphic features with plagiocephaly and scoliosis and is reported to be wheelchair-bound.

Figure 2

Biallelic SYNPO mutation identified in family FA with nephrotic syndrome. (a) Exon structure of SYNPO is shown with arrow indicating position of mutation of patient FA-24 with nephrotic syndrome (NS). Position of start codon (ATG) and stop codon (TAA) are indicated. Exon numbers are marked on a black or white background. Evolutionary conservation is shown for sequence surrounding amino acid position P847 in SYNPO protein. (b) Immunofluorescence staining for Synpo and colocalization in rat glomeruli is shown for different cell type marker proteins: costaining with antibodies against CD31 (endothelial cells), αSMA (mesangial cells), nephrin (podocytes), and Wilms Tumor 1 (WT1) (podocyte nuclei). Insets are shown enlarged in the lower row. Note that there is a strong colocalization of Synpo staining only with nephrin, which labels podocyte slit membrane structures, thereby demonstrating Synpo localization to podocytes but not endothelial or mesangial cells. Bars = 5, 2.5 μm. (c) Single-cell type-specific average expression of Synpo. Data were modified from Karaiskos et al. The heat map is based on Z-scores. Endo, endothelial cell; Mes, mesangial cell; Podo, podocyte; Tub, tubule cell; Imm, immune cell. (d) Upper row: Transfection of a human podocyte cell line with negative control GFP_Mock (green). Cells were almost devoid of large actin fibers in perinuclear cytoplasm but displayed strong F-actin staining (red) within the cell periphery. Lower row: upon transfection with GFP_SYNPO wild-type cDNA (green), SYNPO colocalizes (yellow) with F-actin (red) in thick irregular cytoplasmic perinuclear actin clusters. Bar = 10 μm.

SYNPO Expression in Podocytes

To further characterize the deleteriousness of the patient-derived recessive SYNPO mutation, we performed functional in vitro studies. SYNPO encodes a linear proline-rich protein, highly present in neuronal dendrites and in foot processes of podocytes., Here, we confirm that Synpo colocalizes with nephrin, which, if mutated, causes NS, and which is a marker for podocytes and the glomerular slit membrane (Figure 2b). There was no colocalization with an endothelial cell marker (CD31) or with a mesangial cell marker (αSMA) in rat glomeruli (Figure 2b). Single-cell RNA sequencing databases also show a strong podocytic expression (Figure 2c). Thus, we confirm that Synpo is predominantly expressed in podocytes. SYNPO has been shown to be important for podocyte function and has been extensively studied as an actin-related protein in human podocytes., To determine its subcellular localization on overexpression in cultured human podocytes, we expressed GFP_SYNPO cDNA in a human podocyte cell line and tested for colocalization with phalloidin. On overexpression, SYNPO induced the formation of F-actin networks in podocytes and colocalized with F-actin in thick irregular or parallel arrangements of actin bundles in the perinuclear cytoplasm (Figure 2d, lower row), whereas cells that were transfected with an empty GFP vector did not harbor large actin fiber aggregates in the perinuclear region, but displayed strong F-actin staining within the cell periphery (Figure 2d, upper row).

Wild-Type SYNPO, But Not Mutant, Rescues CDC42 Activation

The small Rho-like GTPases (RHOA, RAC1, and CDC42) play a pivotal role in podocyte morphogenesis, migration, and the pathogenesis of monogenic forms of NS. SYNPO has been shown to induce stress fibers by competitively blocking Smurf1-mediated ubiquitination of RHOA. To examine whether the potential pathogenic mechanism of SYNPO loss-of-function in patient FA-24 is connected to altered RAC1/CDC42 activation, we performed G-LISA assays in cultured human podocytes. In shRNA-mediated SYNPO knockdown podocytes (Supplementary Figure S2), CDC42 activation was significantly reduced compared with scrambled shRNA cells (Figure 3a). The reduction of CDC42 activity was rescued by overexpressing murine Synpo wild-type cDNA. The Synpo cDNA construct reflecting the mutation of the NS patient FA-24 failed to rescue CDC42 activity (Figure 3a). There was no difference in the activation of RAC1 (Supplementary Figure S7).

Figure 3

CDC42 activity and podocyte migration rate (PMR) is reduced by SYNPO knockdown and rescued by transfection with wild-type Synpo, but not by complementary DNA (cDNA), reflecting the mutation of patient with nephrotic syndrome (NS) FA-24. (a) Active levels of CDC42 were measured by CDC42 G-LISA assay in a human podocyte cell line. Short hairpin RNA (shRNA)-mediated knockdown of SYNPO and overexpression of empty vector negative control (+ Mock) reduced active CDC42. Overexpression of mouse wild-type Synpo cDNA (+ WT) rescued this effect. A Synpo cDNA construct reflecting the mutation from NS patient FA-24 (+ P847L) failed to rescue the phenotype. P-values calculated by 1-way analysis of variance. ∗P < 0.01, ∗∗P < 0.01. NS, nonsignificant. (b) In a human podocyte cell line expressing scrambled shRNA, PMR is increased following serum addition (black) relative to serum-deplete conditions (gray). All subsequent experiments were performed in presence of serum. Knockdown of SYNPO (red) showed reduced PMR compared with scrambled shRNA (black). Migration was rescued by overexpression of wild-type Synpo construct (green). However, mouse Synpo cDNA constructs reflecting the mutation (p.P847L) of NS patient FA-24 only partially rescued PMR (pink).

Wild-Type Synpo, But Not the Patient Variant, Rescues Podocyte Migration

Because wild-type Synpo, but not the patient variant, rescued CDC42 activity upon SYNPO knockdown in cultured human podocytes (Figure 3a), we hypothesized that SYNPO may also regulate the podocyte migration rate (PMR), a well-established intermediate phenotype for SRNS disease genes.,, Using the IncuCyte ZOOM live cell imaging system, we found that knockdown of SYNPO reduced PMR compared with scrambled shRNA (Figure 3b). Migration phenotype was rescued by transient overexpression of murine wild-type Synpo cDNA. In contrast, a Synpo cDNA construct reflecting the mutation found in the NS patient failed to rescue the PMR (Figure 3b), suggesting deleteriousness of the mutation of NS family FA.

Discussion

In this study, we generated 2 independent lists of NS candidate genes to augment the discovery of additional monogenic causes of human NS from WES data generated in 1382 families with SRNS. A set of 63 known monogenic causes of NS in mice and another set of 64 genes described to be regulated by the transcription factor WT1, were identified. We then validated the murine and WT1 candidate gene sets by overlapping them with known human NS genes, used as a positive control toward our hypothesis: that the 2 sets of candidate genes may be relevant for the discovery of novel human SRNS genes. Finally, we compared the mouse and WT1 candidate genes (Figure 1) with our unique NS WES-derived candidate list. Hereby, we identified novel candidate genes as potential causes for monogenic human NS, including a SYNPO mutation in a family with steroid-dependent NS. We demonstrate that lack of SYNPO reduced CDC42 activity and PMR in cultured human podocytes, and show that wild-type, but not cDNA representing the patient mutation, was able to rescue the phenotype. This demonstrated the deleteriousness of this mutation in assays relevant for NS pathogenesis. Thus, we conclude that SYNPO may be a novel candidate gene for human NS.

SYNPO mutations have previously been discussed as potential monogenic causes of NS in humans. That assumption is supported by the following publications: (i) Dai et al. published 2 cases with heterozygous mutations in the SYNPO promoter region in patients with NS; however, no homozygous mutations in SYNPO in patients with NS have been reported. (ii) It has been shown that patients with FSGS,67, 68, 69 HIV-associated nephropathy, and IgA nephropathy have reduced expression of SYNPO in the glomerulus. (iii) Srivastava et al. observed a statistical difference of SYNPO expression between patients with FSGS and minimal change disease. (iv) Synpo-deficient mice display impaired recovery from protamine sulfate–induced foot process effacement and lipopolysaccharide-induced NS. Just recently, it was also shown that Synpo-deficient mice demonstrated increased susceptibility to Adriamycin nephropathy. (v) Digenic heterozygosity for Cd2ap and Synpo was sufficient to induce FSGS in mice. (vi) SYNPO was identified as an inhibitor of the Cdc42:IRSp53:Mena signaling, and it was suggested to play a role in a potential antiproteinuric signaling pathway, thus confirming SYNPO as a strong candidate gene for monogenic NS, if mutated.

SYNPO has also been well-studied in podocytes,, and has an established role in actin cytoskeleton regulation., In this context, gene silencing of SYNPO reduced the PMR, whereas the expression of SYNPO led to the induction of stress fibers. We here showed that knockdown of SYNPO in human podocytes leads to an impaired PMR and reduced CDC42 activity (Figure 3a and b). Of note, only wild-type Synpo, but not cDNA representing the patient’s mutation, was able to rescue the phenotype, demonstrating the deleteriousness of the patient’s homozygous allele. This is consistent with findings that implicated dysregulation of Rho-like small GTPases as a role in disease mechanism in several monogenic forms of human and murine SRNS.,,, We therefore suggest that the SYNPO mutation may partly exert its pathogenic effect via dysregulation of CDC42 signaling.

However, it should be emphasized that SYNPO, as of now, can be considered a candidate gene only for monogenic human NS. Identification of additional mutations in SYNPO in other families with NS and functional characterization is essential to substantiate its candidate status. Therefore, the validity of the here-identified mutation must be considered with caution, despite the loss-of-function data we generated.

Monogenic mouse disease genes have successfully been used in the past to identify novel human disease genes.76, 77, 78 Here, we show that 12 of the 63 known mouse NS candidate genes overlap with the 63 known human genes that cause NS if mutated (Figure 1). Before generating the 63 mouse candidate gene list, we had excluded 12 mouse NS genes that were published after the human NS phenotype was discovered, to maintain an unbiased approach (Supplementary Table S3). This emphasizes that the remaining mouse genes, without an associated known human NS phenotype, might be promising candidates for this disease. The fact that equivalent human mutations have not been discovered yet may be partially explained by the fact that childhood-onset SRNS is a rare disease, with a limited number of patients sequenced to date. Furthermore, human and mouse pathophysiology can differ, and some genes may lead to variable disease phenotypes depending on the organism.

In summary, we identified 3 novel candidate genes for human SRNS using 3 independent nonoverlapping hypotheses, based on monogenic mouse models of NS, WT1 downstream targets and unbiased WES, thereby generating further evidence for SYNPO as a potential novel cause of human NS. We demonstrate loss-of-function for the patient mutation and suggest that SYNPO mutations may lead to disease through CDC42 dysregulation.

Disclosure

FH is a co-founder and Scientific Advisory Board member of Goldfinch Biopharma Inc. All the other authors declared no competing interests.