- Research article
- Open Access
Identification of nephropathy candidate genes by comparing sclerosis-prone and sclerosis-resistant mouse strain kidney transcriptomes
BMC Nephrology volume 13, Article number: 61 (2012)
The genetic architecture responsible for chronic kidney disease (CKD) remains incompletely described. The Oligosyndactyly (Os) mouse models focal and segmental glomerulosclerosis (FSGS), which is associated with reduced nephron number caused by the Os mutation. The Os mutation leads to FSGS in multiple strains including the ROP-Os/+. However, on the C57Bl/6J background the mutation does not cause FSGS, although nephron number in these mice are equivalent to those in ROP-Os/+ mice. We exploited this phenotypic variation to identify genes that potentially contribute to glomerulosclerosis.
To identify such novel genes, which regulate susceptibility or resistance to renal disease progression, we generated and compared the renal transcriptomes using serial analysis of gene expression (SAGE) from the sclerosis-prone ROP-Os/+ and sclerosis resistant C57-Os/+ mouse kidneys. We confirmed the validity of the differential gene expression using multiple approaches. We also used an Ingenuity Pathway Analysis engine to assemble differentially regulated molecular networks. Cell culture techniques were employed to confirm functional relevance of selected genes.
A comparative analysis of the kidney transcriptomes revealed multiple genes, with expression levels that were statistically different. These novel, candidate, renal disease susceptibility/resistance genes included neuropilin2 (Nrp2), glutathione-S-transferase theta (Gstt1) and itchy (Itch). Of 34 genes with the most robust statistical difference in expression levels between ROP-Os/+ and C57-Os/+ mice, 13 and 3 transcripts localized to glomerular and tubulointerstitial compartments, respectively, from micro-dissected human FSGS biopsies. Network analysis of all significantly differentially expressed genes identified 13 connectivity networks. The most highly scored network highlighted the roles for oxidative stress and mitochondrial dysfunction pathways. Functional analyses of these networks provided evidence for activation of transforming growth factor beta (TGFβ) signaling in ROP-Os/+ kidneys despite similar expression of the TGFβ ligand between the tested strains.
These data demonstrate the complex dysregulation of normal cellular functions in this animal model of FSGS and suggest that therapies directed at multiple levels will be needed to effectively treat human kidney diseases.
Focal segmental glomerulosclerosis (FSGS) is a leading cause of nephrotic syndrome in adults and is particularly prevalent in the African American population [1–3]. The incidence of FSGS is increasing in young children [4, 5]. Immunosuppressive therapies are the mainstay of treatment, though most patients with FSGS are resistant to treatment . Unfortunately, a large proportion of FSGS patients, who do not respond to treatment, frequently continue to lose kidney function and progress to end stage renal disease (ESRD), requiring dialysis or kidney transplantation. Moreover, recurrent FSGS is a leading cause of proteinuria and allograft loss after kidney transplantation [7–9].
In recent years, genes have been linked or associated with both familial and sporadic FSGS [10–16], which supports the premise that genetic factors play a major pathophysiologic role . To identify novel, candidate FSGS susceptibility genes, we used Serial Analysis of Gene Expression (SAGE) to characterize the kidney transcriptome in the Oligosyndactyly (Os) mouse, a genetic model of oligomeganephronia and progressive FSGS. Homozygous Os/Os mice die in utero at the 64-cells stage of development , but heterozygote Os/+ mice have a skeletal phenotype mainly in the form of fused digits (oligosyndactyly), and a renal phenotype in the form of 50% reduction in nephron number [19, 20]. The ROP/GnLe strain (ROP-Os/+), was established at the Jackson Lab and has been maintained by sibling mating. As the ROP-Os/+ mice age, they develop a histopathologic lesion similar to human FSGS. The Os mutation causes on average 50% reduction of nephron number, an effect which is independent of genetic background. The ROP-Os/+ mice model progressive kidney disease, as demonstrated by glomerulosclerosis, proteinuria, and increase in creatinine which reflects a diminished glomerular filtration rate (, and El-Meanawy, unpublished data). The Os mutation produces a similar renal phenotype when bred on other genetic backgrounds, such as FvB and C3H [20, 22–24]. However, a congenic mouse strain, the C57Bl/6J-Os/+ (C57-Os/+), which has the Os mutation on a C57Bl/6J background, develops oligosyndactyly and oligomeganephronia but not FSGS , thus providing an ideal model for transcriptome comparison with the ROP-Os/+ mouse.
SAGE captures an unbiased, quantitative snapshot of gene expression patterns [25–27]. SAGE libraries are comprised of cDNA sequence tags derived from the 3’ end of the cDNA pool, and a count of SAGE-derived cDNA tags provides a quantitative representation of the corresponding mRNAs in the sample. We previously used SAGE to generate a tag library from normal mouse kidney, which yielded a potassium channel not known to be expressed in the kidney [28, 29]. In the current report, we exploited phenotypic differences between the sclerosis-prone ROP-Os/+ and the sclerosis-resistant C57-Os/+ mice to discover candidate FSGS susceptibility and/or protective genes by comparing the kidney transcriptomes between the two mouse strains using SAGE. To identify pathophysiologically relevant genes, we specifically examined transcriptome profiles in young mice after completion of nephrogenesis, but prior to the development of significant FSGS. Statistical analysis of the SAGE libraries from ROP-Os/+ and C57-Os/+ mice kidneys, followed by assembly of differentially expressed transcripts into networks, allowed identification of genes which potentially play a role in the pathogenesis of glomerulosclerosis in the Os mouse model.
Animals and tissue collection
Animal experiments were conducted in accordance with the National Institute of Health guidelines and approved by the Institutional Animal Care and Utilization committee at Case Western Reserve University and the Medical College of Wisconsin. ROP-Os/+ and C57-Os/+ mice were purchased from Jackson Laboratories (Bar Harbor, ME) and housed in the animal facilities at the above institutions. We phenotyped mice at 6,10,12, and 16 weeks by comparing urine protein to creatinine ratio,glomerular surface area, and sclerosis score in Periodic Acid Schiff (PAS) stained sections as previously described . Glomerular number was determined by the maceration method as previously described . The kidneys were collected for SAGE library construction at 6 weeks of age to avoid identification of developmental genes (because mouse nephrogenesis continues up to 4 weeks of age), and prior to the onset of significant renal disease. After sacrifice, the kidneys were promptly removed and the renal medulla was dissected away and discarded. A small portion of each kidney was fixed in formaldehyde, sectioned, and stained with PAS for histopathological analysis.
Generation of SAGE libraries
Poly(A)+ RNA was isolated from renal cortices of ROP-Os/+ and C57Bl/6-Os/+ mice using ion-exchange column kits from Qiagen (Valencia, CA). The quality of the isolated RNA was analyzed by agarose gel electrophoresis and RNA concentration was determined by UV spectrophotometry . RNA from three mice per strain was pooled and cDNA was synthesized from a starting amount of 5 μg poly(A)+ RNA using cDNA synthesis kit, according to manufacturer recommendation (Invitrogen). A fraction of first and second strand reactions was labeled with 32P] to evaluate synthesis efficiency. cDNA quality was assessed using agarose gel electrophoresis and by measuring radioisotope incorporation. SAGE steps were performed as described previously  using the restriction endonuclease Nla-III as the anchoring enzyme . We used SAGE software to analyze the SAGE produced concatemer sequences, which identifies the cDNA tags flanked by Nla-III restriction enzyme recognition sequence. In addition, the software identifies and excludes duplicate ditags, which are likely to be PCR amplification artifacts. Confirmation of SAGE tag gene identity, as well as relative expression levels were performed using northern blots, RNase protection assay or real-time RT-PCR. The confirmatory methods and data are in a Additional file 1.
Statistical analysis of SAGE tags
We applied several statistical models, which employ Bayesian methods or Poisson distributions, to identify tags that were significantly different between the ROP-Os/+ and C57-Os/+ kidney libraries [33–36]. The method allows comparison of unequal libraries based on the assumption that tag counts are binomially distributed. Regardless of the method used to identify statistically significant differentially expressed tags, rank orders and p-values were comparable.
Annotation of SAGE tags
Accurate annotation of SAGE tag libraries requires a combination of methods. Accordingly, we used multiple approaches to annotate SAGE tags including the online engine “SAGE genie” (http://cgap.nci.nih.gov/SAGE/mSEM), and manual nucleotide blast (blastn) searches. We also constructed genome-wide mouse cDNA data sets from NCBI-Ref-Seq and Riken-Fantom collections and designed a software program, which identifies and catalogues the 3’-most “CATG” sequence string (Nla-III recognition sequence) with the 3’ flanking 10 nucleotides. These tag-to-gene database tables were used to annotate the SAGE tags. Ambiguous tags, which mapped to multiple genes, and single copy tags were excluded. Tags were assigned gene identifiers if the tag mapped to a sequence that contained the most 3’-CATG sequence. Sequence databases were searched in the following priority: 1) our own libraries containing tags extracted from Ref-Seq and Fantom cDNA databases; 2) the mouse mSAGE expression matrix using SAGE Genie; or 3) the mouse mRNA (non-redundant [nr] and expressed sequence tags [dbest]) databases using blastn.
Assembly of transcriptome networks and pathways
We utilized the Ingenuity Pathway Analysis (IPA) system (Ingenuity Systems, Mountain View, CA, USA, http://www.ingenuity.com) to visualize the SAGE data as biological networks. To perform a “connectivity” analysis, a table containing selected gene identifiers, defined as Focus Genes, and their corresponding relative mRNA expression ratios between ROP-Os/+ and C57-Os/+ kidneys was uploaded into the application. Criteria for selection as Focus Genes included: 1) a greater than 2-fold difference in expression levels between the ROP-Os/+ and C57-Os/+ kidney libraries; and/or 2) a p-value < 0.05 for differential expression after correction for multiple testing. Three hundred and nine genes met these criteria, were mapped to their corresponding gene objects in the Ingenuity Pathway Knowledge Base and subsequently overlaid onto a global molecular network developed from information contained in Ingenuity Pathways Knowledge Base. Networks were then algorithmically generated based on the connectivity of the eligible genes. Twenty-eight Focus Genes failed to generate a shared molecular network or pathway. Function of the connectivity networks was inferred from the functions of the genes in the network as annotated in the Knowledge Base.
An additional set of networks was generated using the functional information in the Ingenuity Pathways Knowledge Base (“functional” analysis). Genes from the SAGE kidney libraries, which met the cutoffs defined above and that were associated with biological functions in the Ingenuity Pathways Knowledge Base, were considered for this analysis. Right-tailed Fisher’s exact test was used to calculate a p-value reflecting the probability that each biological function assigned to that network is due to chance alone. Networks were ranked according to the combinatorial p-value of differentially expressed genes represented in the network.
TGFβ Western blot analysis
Western blot of kidney protein extracts were performed as previously described  using TGFβ specific antibody. Equal loading was verified by stripping the membrane and probing for actin expression.
TGFβ activity was determined using mink lung epithelial cells expressing the firefly luciferase reporter gene under the control of the minimal, TGFβ-responsive, plasminogen activator inhibitor-1 (PAI-1) promoter  (gift from Dr. Daniel Rifkin, New York University). The culture conditions were as previously described . The effect of Itch overexpression on TGFβ was determined from two nearly confluent six-well tissue culture plates, which were transiently transfected with either mouse Itch cDNA which was cloned in pCDNA3.1 or empty vector control, using lipofectamine (Invitrogen) according to manufacturer recommendations. Twenty-four hours after transfection, TGFβ (R&D system) was added to the culture media at 200 ng/ml. The cultures were maintained for 8 h and luciferase activity was assayed using a Promega kit according to manufacturer instructions. Relative luminescence units were determined compared to control conditions (empty pCDNA3.1 plasmid transfected cells).
Comparing differentially expressed Os mice genes to human FSGS expression data
We compared 40 SAGE-derived transcripts with the most divergent expression levels between ROP-Os/+ versus C57-Os/+ kidneys, with the human microarray expression data, comprised of samples obtained from micro-dissected human kidney biopsies from patients with FSGS . The fold change and significance (assessed by q-value, corrected for multiple testing) was calculated by the Significance Analysis of Microarrays method (SAM)  using the MeV software from the Institute of Genome Research (http://www.jcvi.org).
Identification of glomerular enriched genes in the Os model of FSGS
To identify genes which are differentially expressed in our model system, and are known to be expressed in the glomerulus, we searched for common genes between the glombase  gene list, SAGE glomerular tags  and annotated transcripts with the most significant differential expression between ROP-Os/+ versus C57-Os/+ kidneys. We then used Pubmed and Ingenuity pathway analysis to mine the literature for possible links between these genes and nephropathy.
We found that the urine protein to creatinine ratio was significantly higher in ROP-Os/+ mice compared to both C57-Os/+ and ROP-+/+ at all time points. The proteinuria in ROP-Os/+ mice reached a peak at 12 weeks. The glomerular surface area at 6 weeks was similar between the ROP-Os/+ and C57-Os/+ (5646 ± 286 and 5708 ± 460 μm2, respectively) and significantly larger than glomerular surface area of ROP-+/+ (3984 ± 456 μm2). However, by 16 weeks of age, the surface area had increased significantly in the ROP-Os/+ compared to C57-Os/+glomeruli (8990 ± 628 vs 6458 ± 1259 μm2). Furthermore, the ROP-Os/+ mice showed progressive increase in sclerosis score over time (Additional file 1).
Os/+ kidney SAGE libraries
We constructed SAGE libraries from renal cortical RNA extracted from ROP-Os/+ and C57-Os/+ kidneys at 6 weeks of age (Table 1), a time point which permitted completion of differentiation, but prior to observation of significant glomerulosclerosis. The combined ROP-Os/+ and C57-Os/+ kidney SAGE libraries contained 49,588 cDNA sequence tags, of which 20,594 were unique. We did not observe “loss of AT-rich tags” bias  and the libraries demonstrated 54.32% A+T and 45.68% C+G nucleotides. Table 2 shows the 50 most statistically significant, differentially expressed and annotated tags. We compared the annotated genes in our SAGE libraries to publically available wild type mouse Nla-III-anchored SAGE libraries. Of the genes not previously described in SAGE kidney libraries, we observed significant expression of arginine-glutamic acid dipeptide (RE) repeats (Rere or Atrophin-2), insulin-like growth factor binding protein (Igfbp7), aminolevulinic acid synthase (Alas1), glutathione S-transferase, theta 1 (Gstt1), and propionyl coenzyme A carboxylase (Pccb). Interestingly, all of these transcripts were induced to a detectable level in ROP-Os/+, suggesting that deep SAGE coverage can unmask changes in kidney gene expression patterns induced in a disease model. To identify differentially expressed glomerular genes, we searched for common genes between glombase , SAGE glomerular tags  and 220 annotated transcripts with the most significant differential expression between ROP-Os/+ versus C57-Os/+ kidneys. This comparison yielded 67 genes, of which 53 were upregulated, and only 14 were downregulated in the ROP-Os/+ mouse (Table 3). We confirmed a reduction in glomerular number in 6 week old ROP-Os/+ mice (data not shown), suggesting that upregulated gene expression is not due to altered glomerular mass.
Assembly of differentially expressed networks and pathways
Using the Ingenuity Pathway Analysis software to perform a connectivity analysis (see Methods), we identified multiple molecular networks of genes differentially expressed in ROP-Os/+ and C57-Os/+ kidneys. These interactions are developed from published literature describing either physical or functional interactions between the molecules. The 13 most significant gene networks, based on connectivity and their imputed molecular and cellular functions, are shown in Table 4; corresponding diagrams of the actual connectivity networks are shown in Figure 1 and the supplemental data (Additional file 1).
Because TGFβ has been implicated in the pathogenesis of glomerulosclerosis, we next chose to overlay the canonical TGFβ signaling network onto our connectivity networks. TGFβ1 signaling pathway target genes, Hnf4a, Itch, and Map3K7ip1, were upregulated in the ROP-Os/+ kidneys, in agreement with published microarray analysis of laser-captured/microdissected human glomeruli from FSGS biopsies . Overlaying TGFβ1-regulated signaling pathways on to the connectivity networks of differentially expressed genes demonstrated that directional changes in gene expression patterns (up or down) in multiple networks (1, 2, 6, 8, and 13) could be a result of TGFβ1 signaling. Figure 2 shows molecular network 1 (Table 4) with the TGFβ signature identified.
We also performed a functional analysis using the Ingenuity Pathway Analysis system (see Methods). These data demonstrated functional clustering of differentially expressed genes within known canonical pathways, and highlighted mitochondrial dysfunction and oxidative stress responses as mechanisms for disease (Figure 2). Interestingly, a Nrf2-dependent stress response, which is the target for a new class of diabetic nephropathy drug , appears to also be activated in ROP-Os/+ kidneys.
Surprisingly, TGFβ Tag was not differentially expressed between the sclerotic and non-sclerotic kidneys. Because the pathway analyses identified a prominent role for TGFβ in ROP-Os/+−dependent FSGS pathogenesis, we there fore compared TGFβ proteinbetween ROP-Os/+ and C57-Os/+ kidneys using Western blot analysis, and found no significant difference between the two strains (Figure 3).
Itch overexpression increases sensitivity to TGFβ invitro
To test the involvement of TGFβ-related candidate genes in glomerular fibrogenesis, we evaluated the effect of Itch, which was upregulated in ROP-Os/+ kidneys, on TGFβ signaling. This gene was selected because Itch, an E3 ubiquitin ligase, amplifies TGFβ signaling by facilitating the interaction between the TGFβ receptor and SMAD2 . For these experiments we used mink lung epithelial cells which express the firefly luciferase under the control of minimal TGFβ-responsive PAI-1 promoter . As seen in Figure 4 cells transfected with Itch cDNA showed higher luciferase activity in response to TGFβ, compared to control cells transfected with empty vector, suggesting that enhanced Itch expression may be relevant to the ROP-Os/+ phenotype by stimulating TGFβ signaling.
Differentially expressed genes in murine FSGS are similarly regulated in human FSGS
To determine if the comparative transcriptome analysis illuminated changes in gene expression in human FSGS, we compared the 40 most significant, differentially expressed genes derived from the ROP-Os/+ vs. C57-Os/+ analysis, with microarray data from kidney biopsy samples derived from patients with FSGS . Of the 40 queried genes, corresponding oligonucleotides from six (ACSM2A, TSPAN33, ERRFI1, NDUFA12, TMEM27, and RHPN2) were not spotted on the array chip. Of the remaining 34 genes, 17 were concordantly regulated in the human glomerulus and ROP-Os/+ mouse kidney, and altered mRNA expression in the human samples was statistically significant for 13 of the 17 genes. Two of the 13 genes were also concordantly regulated in the human tubulointerstitium and ROP-Os/+ mouse kidney (Table 5). An additional concordantly regulated gene, ID2, which encodes inhibitor of DNA binding 2, localized exclusively to the human tubulointerstitium.
Focal and segmental glomerulosclerosis is the histopathologic pattern of a spectrum of renal diseases which start with the glomerulus. There exists the likelihood that multiple insults converge on a common pathogenic pathway. Monogenic forms of familial FSGS due to mutations in TRPC6, α-actinin-4, podocin or APOL1 have been previously described [15, 16, 49]. To identify novel, candidate FSGS susceptibility genes, we compared SAGE-generated kidney transcriptomes from ROP-Os/+ and C57-Os/+mouse strains that are glomerulosclerosis-susceptible and resistant, respectively. Both strains carry the Os mutation which causes reduction in nephron number. The quantitative decrease in nephron number is similar between the two strains, however, the C57-Os/+ strain is protected against the development of renal disease [22, 50]. The glomerulosclerosis phenotype that is associated with the Os-induced reduction in nephron number has been shown in other strains, but it appears the C57-BL/6 genetic background is unique in its resistance to glomerulosclerosis [23, 24]. The increase in glomerular surface area in the ROP-Os/+ strain reflects a hyperfiltration physiology which appear to be not as prominent in the C57-Os/+ strain. Comparison between the ROP-Os/+ and C57-Os/+ transcriptome is therefore a useful approach to identify potential nephropathy genes.
Potentialcandidate nephropathy genes and pathways
There are >100 potentially different expressed genes with a p-value of <0.05 between the mouse strains studied in this report. It is likely that the sclerosis phenotype of the ROP-Os/+ mice (or the resistance of the C57-Os/+ to sclerosis) can be attributed to the combined effect of multiple genes. The renal pathology is due to the expression pattern of these genes on top of the reduced nephron number with its consequent hyperfiltration. The FSGS phenotype is not merely due to difference in the expression of background genes, as we contrasted the wild type ROP versus C57BL/6 kidney transcriptome using microarray, and did not detect differential expression of the SAGE identified genes (data not shown).
In addition to the TGFβ-related genes, which were evaluated in detail in Results, SAGE data showed that neuropilin-2 (Nrp2), which is expressed in the glomerulus and regulates vascular endothelial growth factor (Vegf) signaling [51, 52], is downregulated in the sclerosis prone ROP-Os/+ mouse. Nrp2 knockout mice show progressive glomerular damage after injection of the podocyte toxin, adriamycin, in contrast to their wild type littermates, who recover after the initial injury (J. Sedor, unpublished data).
Pea15 gene, which is overexpressed in the sclerotic mouse kidney, is a death effector domain-containing protein and promoter of autophagy. PEA15 inhibits caspase activation and increases ERK activity . Transgenic mice overexpressing PEA15 have glomerular mesangial expansion and a histological pattern similar to diabetic nephropathy . PEA15 induces the expression of the glucose transporter Glut1 in skeletal muscle cells . Interestingly, Glut1 is also overexpressed in glomeruli from the FvB-Os/+ mice (Os mutation bred on FvB genetic background) , suggesting that PEA15 regulation of Glut1 may play a role in non-diabetic glomerular scarring.
In the ROP-Os/+ strain a number of ubiquitin-proteasome pathway genes are differentially expressed; Psmb5, Psma3, Fbxl12, and Itch gene expression was enhanced, while Siah1a was downregulated. The differential expression of these genes suggests a role for aberrant protein degradation in glomerulosclerosis phenotype of the Os mice. Itch has been linked to the regulation of a multitude of signaling cascades, including TGFβ and EGF through ubiquitin and non-ubiquitin mediated mechanisms [47, 56]. We tested the effect of Itch overexpression in vitro as a proof of concept that it regulates TGFβ signaling, which has been widely implicated in the pathogenesis of glomerulosclerosis. The experimental data support the hypothesis that Itch overexpression increases TGFβ signaling, independent of TGFβ ligand level (Figures 3 and 4). Itch expression is regulated by the Src kinase Fyn [57–60], which regulates podocyte function through phosphorylation of nephrin [61, 62]. Bigenic Fyn/Cd2ap heterozygotes demonstrated an FSGS phenotype . Because of the links with Fyn, TGFβ and other TGFβ-related signaling molecules, such as Pcbp1 , we speculate that Itch may function as a molecular rheostat, by regulating downstream TGFβ signaling (and FSGS pathophysiology) independent of ligand concentration.
Laser capture microdissection-microarray analysis of FSGS glomeruli demonstrated multiple changes consistent with activity of TGFβ signaling . Although we did not identify differences in TGFβ ligand SAGE tag expression in our libraries, and Western blot analysis showed similar TGFβ protein expression between ROP-Os/+ and C57-Os/+ mouse kidneys, we postulate that the TGFβ pathway in the ROP-Os/+ kidney is upregulated by downstream molecules, such as Itch and SnoN.
Using the Ingenuity Pathway Analysis engine, the clustering of genes involved in oxidative stress response like the Nrf2 response genes and GPX suggests a role for an electrophile or oxidative stress in the mechanisms promoting renal injury in the ROP-Os/+ model . Recently, bardoxolone methyl, which activates the Keap1-Nrf2 anti-oxidant pathway, was shown to protect kidney function in patients with type 2 diabetes .
Genes overexpressed in the C57-Os/+ mouse may be protective. For example, glutathione-S-transferase theta (Gstt1) tag counts were 18.6 fold higher in the C57-Os/+ compared to ROP-Os/+ mice, and it was shown that the deletion of this gene increases the likelihood of ESRD in diabetic patients . Also, upregulation of Hsp90ab1, Vcp, Prdx1, and Prdx2 genes in the sclerosis resistant C57-Os/+ strain could indicate a robust protective response to the increased oxidative stress induced by the Os mutation.
Differentially regulated genes known to be involved in renal pathology
Connective tissue growth factor (Ctgf), integrin-β1 (Itgb1), and secreted phosphoprotein 1 (Spp1) were all upregulated in the sclerosis-prone ROP-Os/+ mouse kidney. These genes are also upregulated in a other renal diseases characterized by fibrosis [67–72]. On the other hand, complement factor H (Ctf) and prosaposin (Psap) are relatively downregulated in ROP-Os/+ kidneys, similar to observations in membranoproliferative glomerulonephritis (MPGN) and tubular damage [73, 74].
The ROP-Os/+ mouse and human FSGS
By comparing our ROP-Os/+SAGE libraries to previously published kidney SAGE libraries, we identified 13 genes, which were concordantly regulated in ROP-Os/+ kidneys and human FSGS kidney biopsies, all of which localize to the glomerulus. None of these have been previously detected using either SAGE or microarray libraries from normal kidneys, (http://www.ncbi.nlm.nih.gov/geo/) [43, 75], suggesting that novel glomerular gene sets are induced as part of the pathophysiology of FSGS. Alternatively, the regulation of gene expression in fibrotic glomeruli could be a manifestation (rather than a cause) of fibrogenesis, though we attempted to minimize this possibility by selecting mice for study at an age prior to the appearance of glomerulosclerosis.
Limitations of SAGE
SAGE is a powerful tool for gene and pathway discovery. However, like other global gene expression analysis methods, it has inherent limitations. The lack of a uniform annotation tool for SAGE tags makes data analysis burdensome. Although SAGE permits assembly of a comprehensive transcriptome, the sequencing costs can be significant. As a result, most labs have resorted to microarray technology to achieve adequate, but less comprehensive coverage.
In this paper we have identified multiple candidate nephropathy genes which might regulate FSGS pathogenesis in a mouse model. The identified candidate genes and pathways need to be validated in knockout mice and/or other renal disease models, in order to discover targets for rational drug design or novel renal disease susceptibility markers. The number of networks derived from the differentially expressed gene set underscores the complexity of renal disease pathogenesis, and suggests that therapies directed at multiple pathways will be needed to effectively treat human kidney diseases.
Chun MJ, Korbet SM, Schwartz MM, Lewis EJ: Focal segmental glomerulosclerosis in nephrotic adults: Presentation, prognosis, and response to therapy of the histologic variants. J Am Soc Nephrol. 2004, 15: 2169-2177. 10.1097/01.ASN.0000135051.62500.97.
Bakir AA, Bazilinski NG, Rhee HL, Ainis H, Dunea G: Focal segmental glomerulosclerosis. A common entity in nephrotic black adults. Arch Intern Med. 1989, 149: 1802-1804. 10.1001/archinte.1989.00390080072016.
Bakir AA, Share DS, Levy PS, Arruda JAL, Dunea G: Focal segmental glomerulosclerosis in adult African Americans. Clin Nephrol. 1996, 46: 306-311.
Bonilla-Felix M, Parra C, Dajani T, Ferris M, Swinford RD, Portman RJ, et al: Changing patterns in the histopathology of idiopathic nephrotic syndrome in children. Kidney Int. 1999, 55: 1885-1890. 10.1046/j.1523-1755.1999.00408.x.
Gulati S, Sharma AP, Sharma RK, Gupta A: Changing trends of histopathology in childhood nephrotic syndrome. Amer J Kidney Dis. 1999, 34: 646-650. 10.1016/S0272-6386(99)70388-4.
Gipson DS, Trachtman H, Kaskel FJ, Greene TH, Radeva MK, Gassman JJ, et al: Clinical trial of focal segmental glomerulosclerosis in children and young adults. Kidney Int. 2011, 80: 868-878. 10.1038/ki.2011.195.
Fencl F, Simkova E, Vondrak K, Janda J, Chadimova M, Stejskal J, et al: Recurrence of nephrotic proteinuria in children with focal segmental glomerulosclerosis after renal transplantation treated with plasmapheresis and immunoadsorption: case reports. Transplant Proc. 2007, 39: 3488-3490. 10.1016/j.transproceed.2007.09.045.
Fine RN: Recurrence of nephrotic syndrome/focal segmental glomerulosclerosis following renal transplantation in children. Pediatr Nephrol. 2007, 22: 496-502. 10.1007/s00467-006-0361-6.
Itami N, Akutus Y, Hattori M, Itoh K: Recurrence of nephrotic syndrome in focal segmental glomerulosclerosis after transplantation. Pediatr Nephrol. 1990, 4: 298-299.
Kaplan JM, Kim SH, North KN, Rennke H, Correia LA, Tong HQ, et al: Mutations in ACTN4, encoding à-actinin-4, cause familial focal segmental glomerulosclerosis. Nat Genet. 2000, 24: 251-256. 10.1038/73456.
Winn MP, Conlon PJ, Lynn KL, Farrington MK, Creazzo T, Hawkins AF, et al: A Mutation in the TRPC6 Cation Channel Causes Familial Focal Segmental Glomerulosclerosis. Science. 2005, 308: 1801-1804. 10.1126/science.1106215.
Reiser J, Polu KR, Moller CC, Kenlan P, Altintas MM, Wei C, et al: TRPC6 is a glomerular slit diaphragm-associated channel required for normal renal function. Nat Genet. 2005, 37: 739-744. 10.1038/ng1592.
Kopp JB: Gene expression in kidney using transgenic approaches. Exp Nephrol. 1997, 5: 157-167.
Kao WH, Klag MJ, Meoni LA, Reich D, Berthier-Schaad Y, Li M, et al: MYH9 is associated with nondiabetic end-stage renal disease in African Americans. Nat Genet. 2008, 40: 1185-1192. 10.1038/ng.232.
Genovese G, Friedman DJ, Ross MD, Lecordier L, Uzureau P, Freedman BI, et al: Association of trypanolytic ApoL1 variants with kidney disease in African Americans. Science. 2010, 329: 841-845. 10.1126/science.1193032.
Tzur S, Rosset S, Shemer R, Yudkovsky G, Selig S, Tarekegn A, et al: Missense mutations in the APOL1 gene are highly associated with end stage kidney disease risk previously attributed to the MYH9 gene. Hum Genet. 2010, 128: 345-350. 10.1007/s00439-010-0861-0.
Butkus DE: Familial clustering of end-stage renal disease in Mississippi. J Miss State Med Assoc. 2002, 43: 71-77.
Magnuson T, Epstein CJ: Oligosyndactyly: a lethal mutation in the mouse that results in mitotic arrest very early in development. Cell. 1984, 38: 823-833. 10.1016/0092-8674(84)90277-0.
Stewart AD, Stewart J: Studies on syndrome of diabetes insipidus associated with oligosyndactyly in mice. Am J Physiol. 1969, 217: 1191-1198.
He C, Zalups RK, Henderson DA, Striker GE, Striker LJ: Molecular analysis of spontaneous glomerulosclerosis in Os/+ mice, a model with reduced nephron mass. Am J Physiol. 1995, 269: F266-F273.
Zalups RK: The Os/+ mouse: a genetic animal model of reduced renal mass. Am J Physiol. 1993, 264: F53-F60.
Esposito C, He CJ, Striker GE, Zalups RK, Striker LJ: Nature and severity of the glomerular response to nephron reduction is strain-dependent in mice. Am J Pathol. 1999, 154: 891-897. 10.1016/S0002-9440(10)65336-9.
Lenz O, Zheng F, Vilar J, Doublier S, Lupia E, Schwedler S, et al: The inheritance of glomerulosclerosis in mice is controlled by multiple quantitative trait loci. Nephrol Dial Transplant. 1998, 13: 3074-3078. 10.1093/ndt/13.12.3074.
Wang Y, Heilig KO, Minto AW, Chen S, Xiang M, Dean DA, et al: Nephron-deficient Fvb mice develop rapidly progressive renal failure and heavy albuminuria involving excess glomerular GLUT1 and VEGF. Lab Invest. 2009, 90: 83-97.
Velculescu VE, Zhang L, Vogelstein B, Kinzler KW: Serial analysis of gene expression. Science. 1995, 270: 484-487. 10.1126/science.270.5235.484.
Velculescu VE, Madden SL, Zhang L, Lash AE, Yu J, Rago C, et al: Analysis of human transcriptomes. Nat Genet. 1999, 23: 387-388.
Velculescu VE, Zhang L, Zhou W, Vogelstein J, Basrai MA, Bassett DE, et al: Characterization of the yeast transcriptome. Cell. 1997, 88: 243-251. 10.1016/S0092-8674(00)81845-0.
El-Meanawy MA, Schelling JR, Barathan S, Iyengar SK, Sedor JR: Serial analysis of gene expression (SAGE). Methods in Molecular Medicine: Renal Disease. Guide to Technical Approaches. Volume 86. Edited by: Goligorsky MS. 2003, The Humana Press, Totowa, 1
El-Meanawy MA, Schelling JR, Pozuelo F, Churpek MM, Ficker EK, Iyengar S, et al: Use of serial analysis of gene expression to generate kidney expression libraries. Am J Physiol Renal Physiol. 2000, 279: F383-F392.
He CJ, Esposito C, Phillips C, Zalups RK, Henderson DA, Striker GE, et al: Dissociation of glomerular hypertrophy, cell proliferation, and glomerulosclerosis in mouse strains heterozygous for a mutation (Os) which induces a 50% reduction in nephron number. J Clin Invest. 1996, 97: 1242-1249. 10.1172/JCI118539.
Kaufman JM, Hardy R, Hayslett JP: Age-dependent characteristics of compensatory renal growth. Kidney Int. 1975, 8: 21-26. 10.1038/ki.1975.72.
Sambrook J, Fritsch EF, Maniatis T: Molecular cloning: a laboratory manual. 1989, Cold Springs Harbor, Cold Springs Harbor, 2
Cooper GF, Herskovits E: A Bayesian method for the induction of probabilistic networks from data. Mach Learn. 1992, 9: 309-347.
Schoenberg BS: Calculating Confidence Intervals for Rates and Ratios. Neuroepidemiology. 1983, 2: 257-265. 10.1159/000110529.
van Kampen AH, van Schaik BD, Pauws E, Michiels EM, Ruijter JM, Caron HN, et al: USAGE: a web-based approach towards the analysis of SAGE data. Serial Analysis of Gene Expression. Bioinformatics. 2000, 16: 899-905. 10.1093/bioinformatics/16.10.899.
Kal AJ, van Zonneveld AJ, Benes V, van den Berg M, Koerkamp MG, Albermann K, et al: Dynamics of gene expression revealed by comparison of serial analysis of gene expression transcript profiles from yeast grown on two different carbon sources [In Process Citation]. MolBiol Cell. 1999, 10: 1859-1872.
Laemmli EK: Cleavage of structural proteins during the assembly of the head of bacteriophage T4. Nature. 1970, 227: 680-685. 10.1038/227680a0.
Abe M, Harpel JG, Metz CN, Nunes I, Loskutoff DJ, Rifkin DB: An Assay for Transforming Growth Factor-β Using Cells Transfected with a Plasminogen Activator Inhibitor-1 Promoter-Luciferase Construct. Anal Biochem. 1994, 216: 276-284. 10.1006/abio.1994.1042.
Khan S, Lakhe-Reddy S, McCarty JH, Sorenson CM, Sheibani N, Reichardt LF, et al: Mesangial Cell Integrin [alpha]v[beta]8 Provides Glomerular Endothelial Cell Cytoprotection by Sequestering TGF-[beta] and Regulating PECAM-1. Am J Pathol. 2011, 178: 609-620. 10.1016/j.ajpath.2010.10.031.
Berthier CC, Zhang H, Schin M, Henger A, Nelson RG, Yee B, et al: Enhanced expression of Janus kinase-signal transducer and activator of transcription pathway members in human diabetic nephropathy. Diabetes. 2009, 58: 469-477.
Tusher VG, Tibshirani R, Chu G: Significance analysis of microarrays applied to the ionizing radiation response. Proc Natl Acad Sci USA. 2001, 98: 5116-5121. 10.1073/pnas.091062498.
He L, Sun Y, Patrakka J, Mostad P, Norlin J, Xiao Z, et al: Glomerulus-specific mRNA transcripts and proteins identified through kidney expressed sequence tag database analysis. Kidney Int. 2007, 71: 889-900. 10.1038/sj.ki.5002158.
Nystrom J, Fierlbeck W, Granqvist A, Kulak SC, Ballermann BJ: A human glomerular SAGE transcriptome database. BMC Nephrol. 2009, 10: 13-27. 10.1186/1471-2369-10-13.
Margulies EH, Kardia SL, Innis JW: Identification and prevention of a GC content bias in SAGE libraries. Nucleic Acids Res. 2001, 29: E60-10.1093/nar/29.12.e60.
Bennett MR, Czech KA, Arend LJ, Witte DP, Devarajan P, Potter SS: Laser Capture Microdissection-Microarray Analysis of Focal Segmental Glomerulosclerosis Glomeruli. Nephron Exp Nephrol. 2007, 107: e30-e40. 10.1159/000106775.
Pergola PE, Krauth M, Huff JW, Ferguson DA, Ruiz S, Meyer CJ, et al: Effect of bardoxolone methyl on kidney function in patients with T2D and Stage 3b-4 CKD. Am J Nephrol. 2011, 33: 469-476. 10.1159/000327599.
Bai Y, Yang C, Hu K, Elly C, Liu YC: Itch E3 ligase-mediated regulation of TGF-beta signaling by modulating smad2 phosphorylation. Mol Cell. 2004, 15: 825-831. 10.1016/j.molcel.2004.07.021.
Reich HN, Tritchler D, Cattran DC, Herzenberg AM, Eichinger F, Boucherot A, et al: A Molecular Signature of Proteinuria in Glomerulonephritis. PLoS One. 2010, 5: e13451-10.1371/journal.pone.0013451.
Kaplan J, Pollak MR: Familial focal segmental glomerulosclerosis. Curr Opin Nephrol Hypertens. 2001, 10: 183-187. 10.1097/00041552-200103000-00005.
Striker LJ: Nephron reduction in man–lessons from the Os mouse. Nephrol Dial Transplant. 1998, 13: 543-545. 10.1093/ndt/13.3.543.
Barathan S, Konieczkowski M, El-Meanawy MA, Khan S, Wang B, Dodig T, et al: Expression analysis demonstrates the cell guidance signal neuropilin 2b is reduced in sclerosis-prone ROP-Os/+ but not sclerosis-resistant C57-Os/+ mouse kidneys. J Am Soc Nephrol. 2001, 12: 809A-Ref Type: Abstract
Gluzman-Poltorak Z, Cohen T, Herzog Y, Neufeld G: Neuropilin-2 is a receptor for the vascular endothelial growth factor (VEGF) forms VEGF-145 and VEGF-165. J BiolChem. 2000, 275: 29922-
Bock BC, Tagscherer KE, Fassl A, Kr+ñmer A, Oehme I, Zentgraf HW, et al: The PEA-15 Protein Regulates Autophagy via Activation of JNK. J BiolChem. 2010, 285: 21644-21654.
Oriente F, Iovino S, Cassese A, Romano C, Miele C, Troncone G, et al: Overproduction of phosphoprotein enriched in diabetes (PED) induces mesangial expansion and upregulates protein kinase C-beta activity and TGF-beta1 expression. Diabetologia. 2009, 52: 2642-2652. 10.1007/s00125-009-1528-z.
Condorelli G, Vigliotta G, Iavarone C, Caruso M, Tocchetti CG, Andreozzi F, et al: PED/PEA-15 gene controls glucose transport and is overexpressed in type 2 diabetes mellitus. EMBO J. 1998, 17: 3858-3866. 10.1093/emboj/17.14.3858.
Omerovic J, Santangelo L, Puggioni EM-R, Marrocco J: Dall'Armi C, Palumbo C et al.: The E3 ligase Aip4/Itch ubiquitinates and targets ErbB-4 for degradation. The. FASEB J. 2007, 21: 2849-2862. 10.1096/fj.06-7925com.
Azakir BA, Angers A: Reciprocal regulation of the ubiquitin ligase Itch and the epidermal growth factor receptor signaling. Cell Signal. 2009, 21: 1326-1336. 10.1016/j.cellsig.2009.03.020.
Lallemand F, Seo SR, Ferrand N, Pessah M, L'Hoste S, Rawadi G, et al: AIP4 restricts transforming growth factor-beta signaling through a ubiquitination-independent mechanism. J BiolChem. 2005, 280: 27645-27653.
Yang C, Zhou W, Jeon MS, Demydenko D, Harada Y, Zhou H, et al: Negative regulation of the E3 ubiquitin ligase itch via Fyn-mediated tyrosine phosphorylation. Mol Cell. 2006, 21: 135-141. 10.1016/j.molcel.2005.11.014.
Magnifico A, Ettenberg S, Yang C, Mariano J, Tiwari S, Fang S, et al: WW domain HECT E3s target Cbl RING finger E3s for proteasomal degradation. J BiolChem. 2003, 278: 43169-43177.
Verma R, Wharram B, Kovari I, Kunkel R, Nihalani D, Wary KK, et al: Fyn binds to and phosphorylates the kidney slit diaphragm component nephrin. J BiolChem. 2003, 278: 20716-20723.
Azhibekov TA, Wu Z, Padiyar A, Bruggeman LA, Simske JS: TM4SF10 and ADAP interaction in podocytes: role in Fyn activity and nephrin phosphorylation. Am J Physiol Cell Physiol. 2011, 301: C1351-C1359. 10.1152/ajpcell.00166.2011.
Huber TB, Kwoh C, Wu H, Asanuma K, Godel M, Hartleben B, et al: Bigenic mouse models of focal segmental glomerulosclerosis involving pairwise interaction of CD2AP, Fyn, and synaptopodin. J Clin Invest. 2006, 116: 1337-1345. 10.1172/JCI27400.
Thiele BJ, Doller A, Kahne T, Pregla R, Hetzer R, Regitz-Zagrosek V: RNA-binding proteins heterogeneous nuclear ribonucleoprotein A1, E1, and K are involved in post-transcriptional control of collagen I and III synthesis. Circ Res. 2004, 95: 1058-1066. 10.1161/01.RES.0000149166.33833.08.
Nkova-Kostova AT, Liby KT, Stephenson KK, Holtzclaw WD, Gao X, Suh N, et al: Extremely potent triterpenoid inducers of the phase 2 response: correlations of protection against oxidant and inflammatory stress. Proc Natl Acad Sci USA. 2005, 102: 4584-4589. 10.1073/pnas.0500815102.
Yang Y, Kao MT, Chang CC, Chung SY, Chen CM, Tsai JJ, et al: Glutathione S-transferase T1 deletion is a risk factor for developing end-stage renal disease in diabetic patients. Int J Mol Med. 2004, 14: 855-859.
Turk T, Leeuwis JW, Gray J, Torti SV, Lyons KM, Nguyen TQ, et al: BMP Signaling and Podocyte Markers Are Decreased in Human Diabetic Nephropathy in Association With CTGF Overexpression. J HistochemCytochem. 2009, 57: 623-631.
Mason RM: Connective tissue growth factor(CCN2), a pathogenic factor in diabetic nephropathy. What does it do? How does it do it?. J Cell Commun Signal. 2009, 3: 95-104. 10.1007/s12079-009-0038-6.
Umezono T, Toyoda M, Kato M, Miyauchi M, Kimura M, Maruyama M, et al: Glomerular expression of CTGF, TGF-beta 1 and type IV collagen in diabetic nephropathy. J Nephrol. 2006, 19: 751-757.
Hudkins KL, Giachelli CM, Eitner F, Couser WG, Johnson RJ, Alpers CE: Osteopontin expression in human crescentic glomerulonephritis. Kidney Int. 2000, 57: 105-116.
Pichler R, Giachelli CM, Lombardi D, Pippin J, Gordon K, Alpers CE, et al: Tubulointerstitial disease in glomerulonephritis: Potential role of osteopontin (uropontin). Am J Pathol. 1994, 144: 915-926.
Mooney A, Jackson K, Bacon R, Streuli C, Edwards G, Bassuk J, et al: Type IV collagen and laminin regulate glomerular mesangial cell susceptibility to apoptosis via beta(1) integrin-mediated survival signals. Am J Pathol. 1999, 155: 599-606. 10.1016/S0002-9440(10)65155-3.
Brera-Abeleda MA, Nishimura C, Smith JL, Sethi S, McRae JL, Murphy BF, et al: Variations in the complement regulatory genes factor H (CFH) and factor H related 5 (CFHR5) are associated with membranoproliferative glomerulonephritis type II (dense deposit disease). J Med Genet. 2006, 43: 582-589.
Matsuda J, Kido M, Tadano-Aritomi K, Ishizuka I, Tominaga K, Toida K, et al: Mutation in saposin D domain of sphingolipid activator protein gene causes urinary system defects and cerebellar Purkinje cell degeneration with accumulation of hydroxy fatty acid-containing ceramide in mouse. Hum Mol Genet. 2004, 13: 2709-2723. 10.1093/hmg/ddh281.
Lindenmeyer MT, Eichinger F, Sen K, Anders HJ, Edenhofer I, Mattinzoli D, et al: Systematic analysis of a novel human renal glomerulus-enriched gene expression dataset. PLoS One. 2010, 5: e11545-10.1371/journal.pone.0011545.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2369/13/61/prepub
This work was supported in part by NIDDK grants 5K08DK002699 and 5T32DK007470. We acknowledge the support of the Applied Systems Biology Core of O'Brien Kidney Research Core (P30 DK081943).We thank all participating centers of the European Renal cDNA Bank - Kröner-Fresenius biopsy bank (ERCB-KFB) and their patients for their cooperation. We thank Mr. Michael Linnen (currently, at Morehead Associates) for his generous effort in developing the TAG sequence extraction software. We thank Dr. Frank Park (Medical College of Wisconsin) for his critical reading of this manuscript.
The authors declare that they have no competing interests.
AE-M carried out phenotyping of the mice, construction and analysis of SAGE libraries, confirmation of SAGE tags, IPA analysis, TGF-β assays, and manuscript preparation. JS participated in phenotyping, experimental design, and manuscript preparation. JS participated in phenotyping, experimental design, and manuscript preparation. SI assisted with database mining and experimental design. KG conducted statistical analysis of SAGE tags. VN and MK performed the comparative analysis of SAGE tags to human kidney disease genes. PH, SB, FP and EE participated in tag expression confirmation and SAGE tag annotation. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Supplemental methods, real-time PCR primers sequences, and full table of statistically significant, differentially expressed SAGE Tags (Table 2S).(PDF 205 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
About this article
Cite this article
El-Meanawy, A., Schelling, J.R., Iyengar, S.K. et al. Identification of nephropathy candidate genes by comparing sclerosis-prone and sclerosis-resistant mouse strain kidney transcriptomes. BMC Nephrol 13, 61 (2012). https://doi.org/10.1186/1471-2369-13-61
- Ingenuity Pathway Analysis
- Mouse Kidney
- Nephron Number
- Ingenuity Pathway Knowledge
- Mink Lung Epithelial Cell