- Research article
- Open Access
Systems biology approach to identify transcriptome reprogramming and candidate microRNA targets during the progression of polycystic kidney disease
BMC Systems Biology volume 5, Article number: 56 (2011)
Autosomal dominant polycystic kidney disease (ADPKD) is characterized by cyst formation throughout the kidney parenchyma. It is caused by mutations in either of two genes, PKD1 and PKD2. Mice that lack functional Pkd1 (Pkd1-/-), develop rapidly progressive cystic disease during embryogenesis, and serve as a model to study human ADPKD. Genome wide transcriptome reprogramming and the possible roles of micro-RNAs (miRNAs) that affect the initiation and progression of cyst formation in the Pkd1-/-have yet to be studied. miRNAs are small, regulatory non-coding RNAs, implicated in a wide spectrum of biological processes. Their expression levels are altered in several diseases including kidney cancer, diabetic nephropathy and PKD.
We examined the molecular pathways that modulate renal cyst formation and growth in the Pkd1-/-model by performing global gene-expression profiling in embryonic kidneys at days 14.5 and 17.5. Gene Ontology and gene set enrichment analysis were used to identify overrepresented signaling pathways in Pkd1-/-kidneys. We found dysregulation of developmental, metabolic, and signaling pathways (e.g. Wnt, calcium, TGF-β and MAPK) in Pkd1-/-kidneys. Using a comparative transcriptomics approach, we determined similarities and differences with human ADPKD: ~50% overlap at the pathway level among the mis-regulated pathways was observed. By using computational approaches (TargetScan, miRanda, microT and miRDB), we then predicted miRNAs that were suggested to target the differentially expressed mRNAs. Differential expressions of 9 candidate miRNAs, miRs-10a, -30a-5p, -96, -126-5p, -182, -200a, -204, -429 and -488, and 16 genes were confirmed by qPCR. In addition, 14 candidate miRNA:mRNA reciprocal interactions were predicted. Several of the highly regulated genes and pathways were predicted as targets of miRNAs.
We have described global transcriptional reprogramming during the progression of PKD in the Pkd1-/-model. We propose a model for the cascade of signaling events involved in cyst formation and growth. Our results suggest that several miRNAs may be involved in regulating signaling pathways in ADPKD. We further describe novel putative miRNA:mRNA signatures in ADPKD, which will provide additional insights into the pathogenesis of this common genetic disease in humans.
Autosomal dominant polycystic kidney disease (ADPKD) is characterized by fluid-filled cysts that are thought to result from abnormal cell proliferation and deregulated apoptosis, increased secretion of fluids into the tubular lumen, irregular cell-matrix interactions, and defective cellular polarity [1, 2]. Thus, normal parenchyma is replaced by a cystic epithelium and fibrotic tissue . Genetic mutations in PKD1 (encoding polycystin 1; PC-1) are responsible for majority of cases of ADPKD, the remainder are due to loss of PKD2 (encoding polycystin 2; PC-2). Loss of PC-1 or PC2 expression results in disruption of intracellular Ca2+ levels, which may lead to abnormal proliferation of tubule epithelial cells [4–7]. Additionally, the involvement of PC-1 in various pathways related to proliferation, such as G-protein signaling, Wnt signaling, AP-1, and cell cycle arrest has been reported [8–12]. However, the manner in which these diverse pathways are integrated into cellular circuitry and regulated during progression of ADPKD is not well studied.
MicroRNAs (miRNAs) are small endogenous non-protein encoding RNAs that post-transcriptionally modulate gene expression by binding to the 3'UTR of target mRNAs . They are involved in many biological processes including cell differentiation, cell proliferation, cell mobility and apoptosis  and are associated with many diseases including cancer, hypertension, diabetes, and kidney dysfunction . For example, Kato et al reported a role for miR-192 in diabetic nephropathy . Also, over-expression of the miR-17-92 cluster may play a role in renal cell carcinoma . Further, the importance of miRNAs that are expressed in kidney is supported by mouse knockout studies [18–21]. For example, eliminating Dicer, a key enzyme in miRNA biogenesis, from podocytes, a cell type required for the formation of the size exclusion barrier in the glomerulus, results in progressive loss of podocyte function [20, 21]. Some studies have also suggested a role for miRNAs in ADPKD [22–24].
The Pkd1-/-mouse model develops cystic disease caused by mutation of the same gene responsible for the majority of human ADPKD, and provides a system to study the pathogenesis of ADPKD . However, a systematic, large-scale study, elucidating global changes in gene expression during disease progression in the Pkd1-/-mouse model has yet to be reported. Using the Pkd1-/-model to study the pathogenesis of ADPKD offers the ability to compare gene expression in pre-diseased and diseased kidneys. In the current investigation, we (1) use the Pkd1-/-model to explore the transcriptional changes that occur in ADPKD on a whole genome scale, (2) undertake a comparative transcriptomics approach to determine similarities and differences with human ADPKD, and (3) investigate whether these changes might be related to changes in miRNA expression. We systemically predicted the possible miRNAs that may be associated with the changes in mRNA expression levels during disease progression that were determined by gene expression microarray analysis. We predicted miRNAs that could target signaling pathways in ADPKD. Our results suggest that several miRNAs may be involved in regulating the genetic switches in ADPKD. We further describe several miRNAs and putative miRNA-mRNA signatures, which were previously not reported in ADPKD.
Kidneys from wild type (WT; Pkd1+/+) and Pkd1-null, designated as Pkd1-/-, mutant littermates were investigated at the embryonic ages 14 (E14.5) and 17 (E17.5). The kidneys were fixed in 4% paraformaldehyde overnight, paraffin-embedded, and stained with hematoxylin and eosin. The sections were visualized with a Nikon Eclipse 80i microscope and photographed with a Qimaging Retiga 2000R Fast 1394 camera using NIS-Elements Basic Research 2.34 software (Micro Video Instruments, Avon, MA).
Total RNA extraction and Microarrays
Total RNA was extracted from embryonic kidneys at, E14.5 and E17.5, using the Qiagen miRNeasy Mini kit (Qiagen, Valencia, CA). Three pairs of kidneys at E14.5 and E17.5 were processed. Each pair was analyzed separately. The integrity and purity of the mRNA samples were assessed prior to hybridization using Bioanalyser 2100 with mRNA Nanochips (Agilent Technologies). RNA hybridization was performed, and gene expression profiles were determined using Illumina Mouse Sentrix 6 version 2 Beadchips. The data extraction was performed by using an Illumina Bead Studio V188.8.131.52 software with the output being raw, non-normalized bead summary values.
The raw data matrix extracted from Beadstudio was uploaded in Bioconductor R version 2.9.1 for downstream analysis. The raw data were read using the beadarray package available through the Bioconductor project. The raw intensities were background adjusted by the Subtract method. The log2 summarized data were quantile normalized. The limma package  with empirical Bayes method was used to assess the differentially expressed genes.
The data discussed in this publication have been deposited in NCBI's Gene Expression Omnibus (GEO) and are accessible through GEO Series accession number GSE24352 (http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE24352).
Gene set enrichment analysis and Gene Ontology enrichment analysis
Gene set enrichment analysis (GSEA) (http://www.broad.mit.edu/gsea/) was used to identify potential gene pathways and key transcription factors (TFs) that may modulate cystogenesis (Figure 1). The GSEA C2 database (curated gene sets from BioCarta, KEGG, Signaling Pathway database, Signaling Gateway, Signal Transduction Knowledge Environment, Human Protein Reference database, GenMAPP, Sigma-Aldrich Pathways, Gene Arrays, BioScience Corp., Human Cancer Genome Anatomy Consortium database) that includes well-studied metabolic and signaling pathways and published microarray data sets was used for pathway analysis. The GSEA C3 database containing shared and evolutionarily conserved TF binding motifs defined by the TRANSFAC database was used for TF analysis. The description of each gene set can be found on the GSEA Molecular Signatures Database website: http://www.broad.mit.edu/gsea/msigdb/index.jsp Enrichment of Gene Ontology (GO) categories (Figure 1) in the differentially expressed genes was determined using Bioconductor's GOstats package  for both up- and down-regulated genes.
Comparison with other datasets on ADPKD
Two datasets on ADPKD with mutation in Pkd1, Pkd1L3/L3 and human ADPKD  were obtained from the GEO database (http://www.ncbi.nlm.nih.gov/geo/; ). Raw data were log2 summarized and quantile normalized. The limma package  with empirical Bayes method was used to assess the differentially expressed genes between mutants and controls. These were compared with the differentially expressed genes obtained in comparisons using custom written Perl scripts.
As shown in Figure 1, a combinatorial strategy was used where target miRNAs were predicted for the differentially expressed genes using four algorithms, TargetScan (http://www.targetscan.org/; ), miRanda (http://www.microrna.org/; ), miRDB (http://mirdb.org/miRDB/; ) and microT (http://diana.cslab.ece.ntua.gr/microT/; ). To identify the miRNAs commonly predicted by two or more algorithms, results were intersected (Additional file 1) using custom written Perl scripts (Additional file 2).
Quantitative real-time PCR
Total RNA was extracted from Pkd1-/-and WT kidneys. Three independent biological replicates (different from those used for microarrays) were used and all reactions were run in duplicates for all the expression analysis of genes and miRNAs. For all the quantitative real-time-PCR (qPCR) assays, the ABI PRISM 7300 Sequence Detection System was used. The comparative CT method was used to obtain relative quantitation of genes and miRNAs as per the manufacturer's protocol (Qiagen). For gene expression analysis, cDNA was made using Superscript III reverse transcriptase (Invitrogen). qPCRs were performed, using Power SYBR® Green PCR Master Mix (Invitrogen). Primer sets for P2rx7, Cpeb3, Hdac9, Sox6, Ltbp1, Calcr, Pitx2, Fgfr3, Fgf10, Adam22, Ddx3y, F2rl2, Grap2, Edil3, Mysm1, Alg6 and Alg8 are available in Additional file 3. 18S rRNA was used as the endogenous control.
miRNA reverse transcription was performed using the TaqMan microRNA Reverse Transcription Kit and miRNA-specific primers (Applied Biosystems, Foster City, CA). Real-time TaqMan miRNA-assays were used to quantify miRNA expression using the TaqMan Master Mix and validated primer/probe sets (Applied Biosystems). miRNA levels were normalized to snoRNA-202. Cycling conditions for TaqMan PCR consisted of an initial incubation at 50°C for 2 min and 95°C for 15 sec and 60°C for 1 min.
Pkd1-/- mice model and design of experiment
To study the detailed changes in molecular profiles during the progression of ADPKD, we generated and compared gene expression profiles of Pkd1-/-embryonic kidneys and age-matched WT kidneys at two stages, E14.5 and E17.5 (Figure 1; Additional file 4). Pkd1-/-embryos develop rapidly progressive kidney cysts during embryogenesis, with null mutant kidneys showing no cysts at E14.5 and marked cystic changes by E17.5 (Figure 2). Unsupervised hierarchical clustering was able to discriminate WT and mutant kidney samples at both time points (Figure 3). At the same time, the cluster analysis was also able to distinguish changes between E14.5 and E17.5 kidneys (Figure 3). Table 1 shows a summary of number of differentially expressed genes and pathways for all of the comparisons. Genes showing a greater than 2-fold difference in expression between WT and mutant kidneys at E14.5 as well as at E17.5 (empirical Bayes moderated t-statistic, unequal variance, uncorrected p-value ≤0.05), were considered to be differentially expressed whereas in the comparisons of WT kidneys at E14.5 vs E17.5 and Pkd1-/-kidneys at E14.5 vs E17.5, genes with ≥ 2-fold difference in expression at p-value ≤0.05 (corrected for multiple testing by Benjamini-Hochberg method) were considered as differentially expressed. The expression of 454 genes was significantly changed between WT and mutant kidneys at E14.5 (Table 1, Additional file 5) whereas 884 genes were significantly changed at E17.5 between WT and mutant kidneys (Table 1, Additional file 6). The comparison of E14.5 and E17.5 WT kidneys yielded 1189 differentially expressed genes (Additional file 7) whereas 2287 genes were differentially expressed in comparison of Pkd1-/-kidneys at E14.5 vs E17.5 (Table 1, Additional file 8). Comparing differentially expressed genes observed in Pkd1-/-at E14.5 vs E17.5 (2287 genes) and in WT at E14.5 vs E17.5 (1189 genes) identified genes that were specifically changing during development in diseased (Additional file 9) or healthy conditions (Additional file 10) as shown in Figure 4. This comparison also indicated genes that were regulated during aging from E14.5 to E17.5 regardless of the genotype (Figure 4; Additional file 11). These results indicate that maturation accounted for the greatest number of changes in gene expression, more so than the cyst formation. Nevertheless, 1397 genes (Additional file 9) could be identified for which changes in gene expression levels were specific for Pkd1-/-kidneys. Thus, this analysis served to identify a set of genes specifically changing in PKD that can yield insights about the regulation of gene expression during cystogenesis.
Profiling gene expression changes in Pkd1-/- mutants at E14.5 and E17.5
At E14.5, Pkd1-/-mutants do not exhibit any cyst formation (Figure 2). Therefore gene expression changes at E14.5 may provide an indication of signaling pathways that cause cysts rather than being a consequence of cyst formation. Among the most interesting genes identified by microarrays were P2rx7, Cer1, Frzb, Wnt7b, Dvl3, Gpr34, and Gpr116. These genes are representative of calcium, Wnt and GPCR signaling, and their roles are discussed later.
Pkd1-/-mutant embryos show numerous large renal cysts at E17.5. Microarray analysis showed up-regulation of multiple developmental genes such as Bmp8b, Grem1, Mmp12, Sema3c, Sema4c and Sema6c and transcription factors (TFs)- Gli2, Foxo1, Foxp2, Hoxb8, Pou2f1 and Zbtb32 at E17.5 in Pkd1-/-mutant kidneys compared to WT. Many of these genes are essential for ureteric bud formation, outgrowth, and branching during kidney development. On the other hand, the expression of genes associated with the differentiation of specific nephron segments, such as Umod, Pck1, Fmo2, Slc12a3, Pvalb, and Angpt2 were down-regulated in the null mutants. These data suggest that Pkd1-/-mutants have failed to maintain nephron segment-specific transcriptional signatures. Additionally, we found changes in expression of genes previously related to cyst growth and progression such as components of MAPK and JAK/STAT signaling (these are detailed later in the following sections).
Gene pathway analysis using GO term enrichment and Gene Set Enrichment Analysis
Gene Ontology (GO) and Gene Set Enrichment Analysis (GSEA; http://www.broad.mit.edu/gsea/) are two bioinformatics approaches to identify pathways involved in specific biological processes such as development and disease. GO term enrichment uses association of Gene Ontology terms to genes in a selected gene list. This method discovers what a set of genes may have in common by examining annotations and finding significant shared GO terms whereas GSEA determines whether a known set of genes shows statistically significant, concordant differences between two biological conditions, e.g. WT and Pkd1-/-.
Table 2 and Additional file 12 show the number of enriched GO terms associated with up- and down-regulated genes for each comparison. Of the 1892 gene sets tested in GSEA, we found that 117 (46 up- and 71 down-regulated) pathways at E14.5 and 75 (52 up- and 23 down-regulated) pathways at E17.5 were dysregulated in the mutants. We defined overrepresented pathways by a nominal (NOM) p-value ≤ 0.05. Table 3 and Additional file 13 show enriched gene sets in mutants for up- and down-regulated genes at E14.5 and E17.5. Some pathways may be represented by multiple independent gene sets.
We found that most of the down-regulated gene sets in Pkd1-/-mutants compared to WT at E14.5 and E17.5 represent metabolic pathways (Table 3 and Additional file 13), as suggested by GSEA and GO analysis. In contrast, both the GSEA and GO analysis of up-regulated gene sets in Pkd1-/-mutants compared to WT at E14.5 and E17.5 suggested that Pkd1-/-mutants displayed a rich gene transcriptional profile for kidney development and regeneration, including Wnt, calcium, MAPK and TGFβ signaling (Additional file 13). Among the pathways identified by GO and GSEA were the following:
Activation of mitogenic signaling pathways
We identified 11 up-regulated gene sets in Pkd1-/-mutants compared to WT associated with mitogenic signaling at E17.5, including those associated with growth factor/receptor tyrosine kinase (RTK) signaling (e.g. FGFs, EGF) and extracellular cellular matrix (ECM)/integrin signaling. Previous studies have suggested that activation of EGFs, FGFs and their receptors could promote tubular epithelial cell proliferation and cyst formation in PKD [35–37]. We found Fgfr1 (3.10 fold), Fgfr3 (10.16 fold), Fgf10 (4.18 fold), Nrg1 (3.37 fold) and Prkcb (2.56 fold) up-regulated in Pkd1-/-mutants compared to WT at E17.5. Although gene sets associated with G protein-coupled receptor (GPCR) were not definitively enriched (p = 0.07) in the cystic kidneys at E17.5 (two gene sets associated with GPCR were enriched at E14.5 in Pkd1-/-mutants compared to WT) as shown in Additional file 13, we identified 17 up-regulated individual genes associated with GPCR. Of interest, genes involved in cAMP mediated signaling and calcium regulation such as Pde4b (7.22 fold), Calcr (6.14 fold), and Sstr2 (11.4 fold) were all up-regulated, while negative regulator of GPCR signaling such as Gprasp1 (59.82 fold), and Rgs3 (5.14 fold)  were down-regulated. These results suggest changes in GPCR signaling that might be associated with intracellular cAMP and calcium regulation in the kidneys of Pkd1-/-animals.
Increased cAMP has been shown to promote renal cystic epithelial proliferation in PKD . We observed up-regulation of adenylyl cyclase, Adcy7, which may increase cAMP production in the renal cysts, which in turn promotes renal cystic epithelial proliferation in PKD . The change in Ca2+ homeostasis is an important feature of ADPKD and may lead to increased levels of cAMP [39–41]. Although the gene sets for calcium/calcineurin/NFAT signaling were not significantly enriched (p = 0.14), we found that Crebbp (2.08 fold), P2rx7 (11.66 fold), Prkcb (2.55 fold) and Traf2 (2.14 fold) were up-regulated in the Pkd1-/-cystic kidneys, suggesting a situation where limitation of intracellular Ca2+ could promote cAMP mediated activation of B-Raf/MEK/ERK pathway and cellular proliferation [39–41].
Previous studies suggested that aberrant activation of ERK/MAPK signaling might modulate cyst growth in ADPKD. We found three gene sets for ERK/MAPK signaling cascades up-regulated in Pkd1-/-kidneys compared to WT at E17.5 (Table 3 and Additional file 13) including up-regulation of Fgfr3 (10.16-fold), Map3k12 (2.32-fold), Mink1 (23.48-fold), Rps6ka5 (2.14-fold), and Fosb (2.25-fold). Although we found slight up-regulation of Akt1 (1.71-fold) and Eif4e (2.33-fold), no definitive enrichment of the mTOR pathway could be determined. This may be important from a therapeutic perspective: the mTOR pathway has been considered as a drug target in arresting ADPKD progression, despite recent studies in humans with ADPKD that did not demonstrate a therapeutic effect of Sirolimus and Everolimus, inhibitors of the mTOR pathway, on cyst progression [42, 43].
Activation of angiogenic and immune/inflammatory pathways
Activation of immunoregulatory and inflammatory pathways may be involved in cyst growth in ADPKD. Multiple up-regulated gene sets in the Pkd1-/-kidneys compared to WT at E17.5 were associated with immune/inflammatory (JAK-STAT (n = 3), immunoregulation and inflammation (n = 2)) responses including up-regulation of Il1rl1 (3.79-fold), Il28ra (2.26-fold), Il4i1 (7.35-fold), Jakmip1 (2.6-fold) and Jakmip2 (2.57-fold). We found increased expression of cytokines such as Cxcl1 (2.19-fold), Cxcl16 (2.26-fold), and Clcf1 (8.57-fold). We also found two dysregulated pathways for aging in Pkd1-/-kidneys compared to WT at E17.5 (Additional file 13).
Activation of HDAC inhibitor pathways
Recently it has been shown that inhibition of class I HDACs is able to suppress Pkd2 related phenotypes . Our studies have not implicated any class I HDACs in the pathogenesis of PKD1. On the other hand, we found down-regulation of a member of the class II HDACs, Hdac9 (22.29-fold) in Pkd1-/-mutants compared to WT at E17.5. Additionally, we found multiple up-regulated gene sets belonging to the class of histone deacetylase inhibitors (n = 7; Additional file 13) in Pkd1-/-kidneys compared to WT at E17.5. HDAC9 plays a role in hematopoiesis, and its deregulated expression, along with altered expression of TGF-β2, may be associated with human cancer and Peters' anomaly [44, 45]. As reported in the KEGG pathway database (http://www.genome.jp/kegg/pathway.html), histone deacetylases are involved in many pathways including signal transduction, the notch signaling pathway, cell growth, and death/cell cycle. Their precise role in Pkd1 dependent ADPKD remains obscure. However, strong differential regulation of genes involved in the class II HDAC regulatory cascade suggests that histone modification may be an important event in transcriptome reprogramming during PKD progression.
Comparison with other ADPKD data sets in mouse and human
We aimed to derive correlations between the gene expression changes in PKD and a set of pathways that modulate disease severity in our mouse model and human ADPKD. We compared our study with another published study using the Pkd1L3/L3mouse model  and with the data set available on human ADPKD . Chen et al. described members of three pathways, Wnt, Notch, and BMP, as differentially regulated in the Pkd1L3/L3model compared to the control WT . Our data from the Pkd1-/-model confirmed the differential regulation of members of these pathways (Additional files 14 and 15). A total of 102 dysregulated genes were common between the Pkd1L3/L3model data and our set of differentially expressed genes at E17.5 (comparison: Pkd1-/-mutants vs. WT), whereas 36 dysregulated genes were common between our Pkd1-/-model at E14.5 (comparison: Pkd1-/-mutants vs. WT) and the Pkd1L3/L3model. A ~20% commonality was observed at the pathway level between our Pkd1-/-model and the Pkd1L3/L3model. The comparison between our data set (differentially expressed genes at E17.5) with human ADPKD data  showed a 50% overlap between significantly enriched pathways (Table 4, Additional file 16), which included several important pathways such as calcium signaling, Wnt, MAPK, TGFβ signaling pathways, immune/inflammatory responses, and Notch signaling pathways. On the individual gene level, 20% (n = 314) of genes were similarly changed in both our mouse model and human ADPKD (Table 4, Additional file 16).
Transcription factor analysis of Pkd1-/- kidneys
Computational analysis of genetic regulatory regions of co-expressed genes can furnish additional information on the regulation of a gene set by specific transcription factors (TFs) . Using GSEA, we searched for overrepresented TF promoter binding motifs among differentially expressed genes (Table 1). We defined overrepresented TF binding motifs by a NOM p-value ≤ 0.05. Additional file 17 lists all the dysregulated gene sets with shared TF binding sites for renal development, mitogen-mediated proliferation, cell cycle, epithelial-mesenchymal transition, angiogenesis, and immune/inflammatory response. Several of these enriched transcription factors, such as Gli2 (4.39 fold), Foxo1 (2.51 fold) and Pou2F1 (2.88 fold) were also differentially regulated in our microarray analyses, providing additional evidence for their functional relevance in PKD.
Prediction of miRNA and miRNA-target interactions, and analysis of miRNA expression
Large scale transcriptional reprogramming during the progression of PKD is consistent with a complex, multi-layered regulatory process, one of the possible regulators of which are miRNAs. Insights into the biological pathways potentially regulated by miRNAs in PKD were obtained by predicting target miRNAs (Additional file 1) using four prediction tools, TargetScan , miRanda , miRDB  and microT  for the differentially expressed genes in cystic kidneys at E14.5 and E17.5 (Table 5). Results were overlapped using custom written Perl scripts (Additional file 2). Integrating the results from several prediction algorithms helps in reducing false positives . More than two tools predicted over 5000 mRNA-miRNA interactions, which represented approximately 500 dysregulated genes in cystic kidneys were predicted targets of at least one miRNA. A total of 344 miRNAs were predicted to target 179 out of 454 dysregulated genes at E14.5 by different algorithms (Table 5a and Additional file 18), whereas 372 out of 884 dysregulated genes were predicted targets of 424 miRNAs at E17.5 (Table 5b and Additional file 19). Of these 372 and 179 dysregulated genes at E17.5 and E14.5 respectively, about 82% were predicted targets of multiple miRNAs (Table 5, Additional files 18 and 19), thus suggesting a role of miRNAs in PKD transcriptional reprogramming. These results are consistent with the previous studies estimating that nearly all mammalian genes may be regulated by miRNAs .
Validation of expression by qPCR
Using qPCR assays for the Pkd1-/-and WT samples, we confirmed the microarray results on a subset of genes that were- (i) differentially expressed on microarrays, (ii) known/suggested to be involved in ADPKD, and (iii) targets of miRNAs predicted by two or more tools. These genes were P2rx7, Cpeb3, Hdac9, Sox6, Calcr, Pitx2, Fgfr3, Fgf10, Adam22, Ddx3y, F2rl2, Grap2, Edil3, Mysm1 and Alg6 (Figures 5 and6). Qualitatively, the qPCR validation data agreed with the microarray data (Additional file 20). On the other hand, we found that for certain genes such as Adam22, Grap2 and F2rl2 the fold changes observed on microarrays varied from those obtained by qPCR (Additional file 20). Several factors may be responsible for these variations, among which include the distinct platforms used for microarrays and qPCR, differences in qPCR amplicon/primers and hybridization probes, non-specific and cross-hybridization and issues associated with amplification of mRNA samples [46, 47].
Computational identification of a large number of miRNA-mRNA target interactions suggested that the expression of miRNAs might also change during the progression of PKD. We tested this hypothesis by determining the differential expression of 9 miRNAs (mmu-miR-10a, mmu-miR-30a-5p, mmu-miR-96, mmu-miR-126-5p, mmu-miR-182, mmu-miR-200a, mmu-miR-204, mmu-miR-429, and mmu-miR-488) between WT and Pkd1-/-genotypes at E14.5 and E17.5 (Figures 7 and8). These nine miRNAs are expressed in kidney , but have not been previously associated with ADPKD. In our computational analysis, two or more tools predicted them as targets for differentially expressed genes at E14.5 and E17.5. The significance of changes in relative expression of miRNAs among the two groups was tested by Student's t-test, and a cutoff of 1.2 fold  and P ≤ 0.05 was used to determine the differential miRNA accumulation.
Among the 9 miRNAs tested at E14.5, only two miRNAs, miR-488 and miR-204, were down-regulated in Pkd1-/-kidneys compared to WT, miR-96 did not change, and the remaining 6 miRNAs were up-regulated in Pkd1-/-kidneys compared to WT (Figure 7). In contrast, at E17.5 miR-96, miR-182 and miR-30a were up-regulated in Pkd1-/-kidneys compared to WT. miR-200a did not show significant variation between the WT and Pkd1-/-kidneys at E17.5 (Figure 8). This indicates that miRNA expression also undergoes reprogramming during PKD. Given that renal cystic- and control- kidneys contained a mixture of cell types, the moderate effects we observed in changes in miRNA levels may in reality be far stronger in specific cell/tissue types.
Most miRNAs negatively regulate the levels of expression of their targets: therefore, miRNA-mRNA target pair should be inversely correlated for their expression levels i.e. if the expression of a given miRNA is up-regulated, levels of its target genes should be down-regulated and vice-versa. Indeed, this was observed for a large number of predicted miRNA:mRNA interactions. At E14.5, expression of 9 miRNAs was inversely related to expression of 46 potential target genes, indicating that these relationships may be functional miRNA-target combinations in ADPKD (Additional file 21). Similarly, 78 genes that were potential targets of 9 miRNAs were inversely related by their expression at E17.5 (Additional file 21). Several of the genes verified with qPCRs (Figures 56), and predicted as targets of the miRNAs, showed a reciprocal relationship (Table 6). For example, down-regulation of miR-200a at E17.5 was correlated with up-regulation of Pitx2 in the Pkd1-/-mutants (Additional file 21). Similarly, down-regulation of miR-10a, miR-126-5p, miR-204, and miR-488 at E17.5 were inversely correlated with up-regulation of Ltbp1, Edil3, P2rx7, and Fgfr3 respectively (Additional file 21). Conversely, we found that up-regulation of miR-30a-5p, miR-96 and miR-182 at E17.5, and miR-429 at E14.5 were reciprocally correlated with down-regulation of Cpeb3, Sox6, Hdac9, and Ddx3y respectively in Pkd1-/-mutants (Additional file 21). Moreover, miRs-10a, -30a-5p, -96, -126-5p, -182, -200a, -204, -429, and -488 have not been previously linked to PKD. Thus, our analysis suggests a potential important role for miRNAs in PKD, though we emphasize that these predicted miRNA:mRNA interaction pairs are subject to further experimental validation.
We analyzed the functional enrichment of predicted and validated (by qPCR) miRNAs for differentially expressed genes in each comparison in an attempt to uncover the functional meaning among these miRNAs (Figure 9). Genes that were commonly targeted by the differentially expressed miRNAs in Pkd1-/-samples were clustered in 18 biochemical pathways. The overall analysis highlighted that signal transduction pathways such as calcium, VEGF, Notch, and MAPK signaling (Figure 9) may be regulated by miRNA. For example, miR-30a-5p may be involved in histone deactylase inhibitor pathways, apoptosis, calcium and Wnt signaling (Figure 9); miR-10a may be involved in TGF-β and hedgehog signaling; miR-204 may be involved in calcium signaling while miR-488 may be involved in MAPK signaling by targeting Fgfr3 (Figure 9). Thus, we generated a comprehensive atlas of miRNA-target genes and pathways during the progression of PKD. This dataset will serve as a useful resource for future investigations.
In the current investigation, we have attempted to uncover the gene expression signatures associated with the initiation and progression of ADPKD, to decipher the signal transduction pathways associated with these changes in gene expression, and to explore the potential role of miRNAs in effecting these gene expression and signal transduction differences between normal and PKD kidneys. We used a combinatorial approach involving microarrays, data mining and prediction of target miRNAs to profile changes in gene expression during progression of ADPKD, integrate various signaling pathways and present a possible cellular signaling circuitry in PKD (Figure 10). In parallel to the mRNA expression profiling, we also determined the changes in expression of several miRNAs that were computationally identified for their probable roles in ADPKD, thus allowing us to predict several miRNA-mRNA reciprocal interactions. By comparing our datasets with those available for human ADPKD, we found common mis-regulated pathways in the developmental and signaling pathways. Identification of commonly regulated pathways and genes between mice and humans is essential for designing effective therapeutics, as these genes/pathways would be targets for pharmacological intervention in humans with PKD.
Transcriptome changes in PKD
A high level of complexity in the molecular pathobiology was found in Pkd1-/-cystic kidneys. We identified genes that were changed only in Pkd1-/-animals such as Aqp7, Bmp6, Bmp8b, Ddx3y, Kl, Psen1 and Pitx2. This shows changes in expression of developmental and nephron-segment specific genes. Further we identified changes in expression of genes that may be associated with pre-cystic stages leading to cyst formation such as changes in the components of Wnt signaling. We identified changes in gene expression related to cyst growth such as components of MAPK (e.g. Map3k12, Fgf10, Prkcb, Traf2) and TGF-β (e.g. Pitx2) signaling that were up-regulated in the Pkd1-/-animals. The GO term enrichment and GSEA revealed changes in calcium, MAPK, Wnt, TGF-β and JAK/STAT signaling pathways in the Pkd1-/-animals. Additionally, this analysis showed that the class II HDAC, Hdac9 was strongly down-regulated in Pkd1-/-animals. This could effect more global changes in gene expression.
Mitogen-Activated Protein Kinases play central roles in the signal transduction pathways that affect gene expression, cell proliferation, differentiation and apoptosis . Although some studies have suggested that cellular proliferation may be involved in initial cyst formation in PKD, other evidence suggests that abnormal proliferation contributes to cyst growth rather than cyst initiation [51–54]. Our data (Table 3 and Additional file 6), showing a greater increase in MAPK-associated pathways at E17.5 as compared with E14.5, are more consistent with the latter hypothesis that proliferation is more important in cyst growth than in initiation.
Wnt signaling in ADPKD
Previous studies have shown that the activation of canonical Wnt signaling  and inhibition of non-canonical Wnt signaling (such as decreased activation of c-jun N-terminal kinase, JNK, and transient changes in intracellular Ca2+ concentrations) may play roles in cyst formation in PKD . We found Wnt7b (1.6 fold), Fzd6 (1.3 fold), Dvl3 (1.6 fold) and Lef1 (1.7 fold) were up-regulated while Apc (2.4 fold), an inhibitor of canonical Wnt signaling was down-regulated at E14.5. Although microarray probes failed to detect increased Wnt7a expression, our independent study confirms up-regulation of Wnt7a in Pkd1-/-animals (Qin et al, submitted). Moreover, in a study of Pkd1L3/L3mouse model Wnt7a and Wnt7b were also found up-regulated in mutant kidneys . The GSEA results also showed enrichment of canonical Wnt signaling for up-regulated genes in Pkd1-/-mutants at E14.5 (Table 3). Canonical Wnt signaling was also up-regulated at E17.5, as determined by both GO analysis and GSEA (Tables 2 and 3, Additional files 12 and 13). Together, these results suggested that there may be an increase in canonical Wnt signaling in Pkd1-/-animals. Additionally, the GO term enrichment analysis showed down-regulation of non-canonical Wnt signaling in Pkd1-/-animals at E14.5 (Table 2, Additional file 12).
In contradistinction to the GO and GSEA analyses suggesting increased Wnt signaling in Pkd1-/-kidneys, we found that two potential inhibitors of Wnt signaling, Cer1 and Frzb (respectively increased 12 fold and 9.7 fold at E14.5), were up-regulated in Pkd1-/-animals at E14.5. Cer1 and Frzb belong to the family of sFRPs (secreted Frizzled Related Proteins) that act as inhibitors of Wnt signaling in animals such as Xenopus and chick by interfering with the binding of Wnt proteins to Frizzled transmembrane receptors . Frzb is a known Wnt antagonist in mammals [57, 58]. However the mouse Cer1 may not inhibit Wnt signaling . At present we are unable to explain the discrepancy between the overall informatic result that canonical Wnt signaling is increased in Pkd1-/-kidneys, and the observed increased expression of Wnt antagonists Cer1 and Frzb. It is possible that Cer1 and/or Frzb have a greater ability to interfere with non-canonical vs. canonical Wnt signaling, and that this will be borne out by future studies.
Calcium signaling in ADPKD
The polycystin-1/2 complex (PC1/PC2) mediates Ca2+ entry into the cell in response to mechanical stimulation of the primary cilium . This in turn triggers additional release of Ca2+ from intracellular stores. Thus, disruption of the polycystin pathway alters intracellular Ca2+ homeostasis [1, 56]. It has been suggested that altered intracellular Ca2+ levels increases cAMP levels by stimulating adenylyl cyclases . In turn, increased cAMP paradoxically stimulates MAPK/ERK signaling and promotes renal cystic epithelial proliferation in Pkd1-/-cells . Our results showed that in Pkd1-/-animals at E17.5 an adenylyl cyclase, Adcy7 was up-regulated. Thus, increased gene expression of an adenylyl cyclase may be an additional contributing factor to increased cAMP levels in ADPKD.
Novel regulatory mechanisms in ADPKD
Our study identifies several new nodes of regulation in ADPKD. For example, we found a strong down-regulation of Alg6 (asparagine-linked glycosylation 6 homolog) in the Pkd1-/-genotype at E17.5. ALGs are important components of the glycosylation pathway involving the endoplasmic reticulum (ER) and the Golgi apparatus. Alg6 encodes a member of the glucosyltransferase family that catalyzes the addition of the first glucose residue to the growing lipid-linked oligosaccharide precursor of N-linked glycosylation. Mutations in this gene are associated with congenital disorders of glycosylation type Ic . Indeed, in Pkd1-/-mice, defective glycosylation of protein has been observed, that may have important implications for ADPKD [62, 63]. As Alg6 encodes an ER glucosyltransferase , its specific downregulation in the Pkd1-/-genotype indicates metabolic reprogramming in the ER and Golgi complex. Thus, exploring Alg6 function in knock-down studies in kidneys may provide new insights regarding the disease processes underlying ADPKD.
The comparison of the Pkd1-/-mouse model and human ADPKD data sets revealed significant overlaps. Some representative of commonly up-regulated genes include Adcy7 (Calcium signaling), Crebbp (Wnt, Jak-STAT, calcium, TGFβ signaling), Cxcl1 (Jak-STAT signaling), Vegfb (VEGF signaling), Dvl2 (Wnt signaling), Tgif1 and Smurf1 (TGFβ signaling). Representative down-regulated genes that were commonly changed between the two data sets include Umod (loop of Henle marker), Rgs3 (GPCR signaling), Pvalb (distal tubule marker), Pck1 (proximal tubule marker), and Lamc3 (Focal adhesion). Thus, these genes may be used as prognostic markers for ADPKD.
Possible regulatory roles of miRNAs in ADPKD
We examined the possible involvement of miRNAs in ADPKD in the Pkd1-/-mouse model. While several miRNAs can target a single transcript, each miRNA can also have multiple targets [13, 31]. Thus, moderate alterations in miRNA levels can have profound effects on the accumulation of their targets [65, 66]. However, their role in ADPKD, both at the individual as well as the systems levels, is still under investigation, as only three studies have directly demonstrated changes in miRNA expression in ADPKD [18, 22, 23]. We have addressed this problem by systemically estimating the possible miRNAs that may target genes deregulated during disease progression (Additional file 1). We have predicted and verified (by qRT-PCR) reciprocal expression of several miRNAs and their predicted targets in Pkd1-/-animals (Table 6). We observed that miRNAs: miRs-10a, -30a-5p, -96, -126-5p, -182, -200a, -204, -429, and -488; and the miRNA-mRNA interactions such as miR-126-5p-Fgf10, miR-488-Fgfr3, miR-182-Hdac9, miR-204-P2rx7 and miR-96-Sox6 (as shown in Table 6) have not been previously reported in ADPKD.
Signaling pathways are ideal candidates for miRNA-mediated regulation, as their components require finely tuned temporal changes in their expression [65, 67]. miRNAs may affect the responsiveness of cells to signaling molecules such as TGFβ, Wnt, Notch, and epidermal growth factor . Once a miRNA targets an inhibitor of a signaling cascade, it serves as a positive regulator by either amplifying signal strength or duration, or empowering cell responsiveness to otherwise sub-threshold stimuli. For example, miR-126 promotes angiogenesis and vascular integrity  by inhibiting the production of natural repressor (SPRED1 and PIK3R2) of VEGF signaling, suggesting that it may serve as an effective target for anti-angiogenic therapies. We found that miR-126-5p may be involved in calcium, EGF, MAPK signaling, and neuroactive ligand-receptor interaction (Figure 9; Additional file 21). Previous studies [13, 65, 67] have shown that a single miRNA can act simultaneously on multiple signaling pathways to coordinate their biological effects and concurrently several miRNAs may regulate single pathway (Figure 9). Moreover, the predicted miRNA:mRNA interactions suggested that the involvement of miRNAs in the signaling pathways through their target mRNAs may lead to the dysregulation of these pathways. For example, the interaction of miR-204:P2rx7 (Figure 10) suggested a possibility that the down-regulation of miR-204 may lead to up-regulation of its target gene, P2rx7 at E17.5 in Pkd1-/-kidneys. P2rx7 is a cell-surface, ligand-gated cation channel and its stimulation leads to alteration in the intracellular Ca2+ levels that may further result into dysregulation of calcium signaling. Therefore, by targeting P2rx7, miR-204 may disrupt calcium signaling. Similarly, at E17.5 in Pkd1-/-animals, the up-regulation of Fgfr3 and Fgf10 (components of MAPK signaling) and down-regulation of their target miRNAs- miR-488 and miR-126-5p respectively, may stimulate MAPK signaling and cell proliferation in Pkd1-/-samples.
The functional correlation between the differentially expressed mRNAs and miRNAs as a module in ADPKD revealed a tight post-transcriptional regulatory network at the mRNA level whose alteration might contribute to increased immune response, by either direct miRNA targeting or through secondary proteins. Further progress in the understanding of miRNA activity, the identification of miRNA signatures in different states, and the advancement of miRNA manipulation techniques will be valuable for deciphering the roles of individual miRNAs in PKD. The overall discovery of differentially regulated miRNAs in the different diseases is expected not only to broaden our biological understanding of these diseases, but more importantly, to identify candidate miRNAs as potential targets for future clinical applications.
We have attempted to reveal the molecular players of PKD in a Pkd1-/-mouse model. Taken together, our study suggests the presence of complex layers of regulation in ADPKD. Our microarray data revealed genes that specifically changed during disease condition. Further we identified genes related to pre-cystic stages and cyst progression. Our model suggests a cascade of signaling events involving up-regulation of canonical- and down-regulation of non-canonical- Wnt signaling that may result in decreased intracellular Ca2+ concentration, resulting in an increase of intracellular cAMP levels that in turn stimulates MAPK/ERK signaling leading to proliferation, followed by increased Jak-STAT signaling and inflammation (Figure 10), ultimately leading to renal failure. Further, we determined gene expression signatures common between the Pkd1-/-mouse model and human ADPKD. These could be used as prognostic markers of disease progression in PKD. Moreover, we add various new components including Alg6, Hdac9 and several miRNAs to the regulatory layers of ADPKD. We predicted that several of the differentially regulated genes are miRNA targets and miRs-10a, -30a-5p, -96, -126-5p, -182, -200a, -204, -429, and -488 may be important players in cellular signalling events leading to PKD. Furthermore, it is interesting to note that these miRNAs have not been previously reported in PKD. It has been proposed that a single miRNA can target more than hundred genes and one gene can be the target of several miRNAs [69, 70]. A future challenge will be to systemically identify all of the miRNAs affecting, and regulated by, the dysregulated cell signalling pathways in ADPKD. Extensive functional analyses of these miRNAs and their target genes by performing knockout and over-expression studies, individually and in combination, are likely to open up new avenues for PKD research.
Torres VE, Harris PC: Polycystic kidney disease: genes, proteins, animal models, disease mechanisms and therapeutic opportunities. J Intern Med. 2007, 261 (1): 17-31. 10.1111/j.1365-2796.2006.01743.x
Wilson PD, Goilav B: Cystic disease of the kidney. Annu Rev Pathol. 2007, 2: 341-368. 10.1146/annurev.pathol.2.010506.091850
Wilson PD: Polycystic kidney disease. N Engl J Med. 2004, 350 (2): 151-164. 10.1056/NEJMra022161
Lantinga-van Leeuwen IS, Dauwerse JG, Baelde HJ, Leonhard WN, van de Wal A, Ward CJ, Verbeek S, DeRuiter MC, Breuning MH, de Heer E, et al.: Lowering of Pkd1 expression is sufficient to cause polycystic kidney disease. Human Molecular Genetics. 2004, 13 (24): 3069-3077. 10.1093/hmg/ddh336
Chang MY, Parker E, Ibrahim S, Shortland JR, Nahas ME, Haylor JL, Ong AC: Haploinsufficiency of Pkd2 is associated with increased tubular cell proliferation and interstitial fibrosis in two murine Pkd2 models. Nephrol Dial Transplant. 2006, 21 (8): 2078-2084. 10.1093/ndt/gfl150
Hanaoka K, Guggino WB: cAMP regulates cell proliferation and cyst formation in autosomal polycystic kidney disease cells. J Am Soc Nephrol. 2000, 11 (7): 1179-1187.
Yamaguchi T, Nagao S, Wallace DP, Belibi FA, Cowley BD, Pelling JC, Grantham JJ: Cyclic AMP activates B-Raf and ERK in cyst epithelial cells from autosomal-dominant polycystic kidneys. Kidney Int. 2003, 63 (6): 1983-1994. 10.1046/j.1523-1755.2003.00023.x
Kim E, Arnould T, Sellin L, Benzing T, Comella N, Kocher O, Tsiokas L, Sukhatme VP, Walz G: Interaction between RGS7 and polycystin. Proc Natl Acad Sci USA. 1999, 96 (11): 6371-6376. 10.1073/pnas.96.11.6371
Kim E, Arnould T, Sellin LK, Benzing T, Fan MJ, Gruning W, Sokol SY, Drummond I, Walz G: The polycystic kidney disease 1 gene product modulates Wnt signaling. J Biol Chem. 1999, 274 (8): 4947-4953. 10.1074/jbc.274.8.4947
Delmas P, Nomura H, Li X, Lakkis M, Luo Y, Segal Y, Fernandez-Fernandez JM, Harris P, Frischauf AM, Brown DA, et al.: Constitutive activation of G-proteins by polycystin-1 is antagonized by polycystin-2. J Biol Chem. 2002, 277 (13): 11276-11283. 10.1074/jbc.M110483200
Parnell SC, Magenheimer BS, Maser RL, Zien CA, Frischauf AM, Calvet JP: Polycystin-1 activation of c-Jun N-terminal kinase and AP-1 is mediated by heterotrimeric G proteins. J Biol Chem. 2002, 277 (22): 19566-19572. 10.1074/jbc.M201875200
Bhunia AK, Piontek K, Boletta A, Liu L, Qian F, Xu PN, Germino FJ, Germino GG: PKD1 induces p21(waf1) and regulation of the cell cycle via direct activation of the JAK-STAT signaling pathway in a process requiring PKD2. Cell. 2002, 109 (2): 157-168. 10.1016/S0092-8674(02)00716-X
Bartel DP: MicroRNAs: Target Recognition and Regulatory Functions. Cell. 2009, 136 (2): 215-233. 10.1016/j.cell.2009.01.002
Shivdasani RA: MicroRNAs: regulators of gene expression and cell differentiation. Blood. 2006, 108 (12): 3646-3653. 10.1182/blood-2006-01-030015
Zhang BH, Pan XP, Cobb GP, Anderson TA: microRNAs as oncogenes and tumor suppressors. Developmental Biology. 2007, 302 (1): 1-12. 10.1016/j.ydbio.2006.08.028
Kato M, Zhang J, Wang M, Lanting L, Yuan H, Rossi J, Natarajan R: MicroRNA-192 in diabetic kidney glomeruli and its function in TGF-beta-induced collagen expression via inhibition of E-box repressors. Proc Natl Acad Sci USA. 2007, 104 (9): 3432-3437. 10.1073/pnas.0611192104
Chow TF, Mankaruos M, Scorilas A, Youssef Y, Girgis A, Mossad S, Metias S, Rofael Y, Honey RJ, Stewart R, et al.: The miR-17-92 cluster is over expressed in and has an oncogenic effect on renal cell carcinoma. J Urol. 2010, 183 (2): 743-751. 10.1016/j.juro.2009.09.086
Agrawal R, Tran U, Wessely O: The miR-30 miRNA family regulates Xenopus pronephros development and targets the transcription factor Xlim1/Lhx1. Development. 2009, 136 (23): 3927-3936. 10.1242/dev.037432
Shi S, Yu L, Chiu C, Sun Y, Chen J, Khitrov G, Merkenschlager M, Holzman LB, Zhang W, Mundel P, et al.: Podocyte-selective deletion of dicer induces proteinuria and glomerulosclerosis. J Am Soc Nephrol. 2008, 19 (11): 2159-2169. 10.1681/ASN.2008030312
Ho J, Ng KH, Rosen S, Dostal A, Gregory RI, Kreidberg JA: Podocyte-specific loss of functional microRNAs leads to rapid glomerular and tubular injury. J Am Soc Nephrol. 2008, 19 (11): 2069-2075. 10.1681/ASN.2008020162
Harvey SJ, Jarad G, Cunningham J, Goldberg S, Schermer B, Harfe BD, McManus MT, Benzing T, Miner JH: Podocyte-specific deletion of dicer alters cytoskeletal dynamics and causes glomerular disease. J Am Soc Nephrol. 2008, 19 (11): 2150-2158. 10.1681/ASN.2008020233
Lee SO, Masyuk T, Splinter P, Banales JM, Masyuk A, Stroope A, Larusso N: MicroRNA15a modulates expression of the cell-cycle regulator Cdc25A and affects hepatic cystogenesis in a rat model of polycystic kidney disease. J Clin Invest. 2008, 118 (11): 3714-3724. 10.1172/JCI34922
Pandey P, Brors B, Srivastava PK, Bott A, Boehn SN, Groene HJ, Gretz N: Microarray-based approach identifies microRNAs and their target functional patterns in polycystic kidney disease. BMC Genomics. 2008, 9: 624- 10.1186/1471-2164-9-624
Tran U, Zakin L, Schweickert A, Agrawal R, Doger R, Blum M, De Robertis EM, Wessely O: The RNA-binding protein bicaudal C regulates polycystin 2 in the kidney by antagonizing miR-17 activity. Development. 2010, 137 (7): 1107-1116. 10.1242/dev.046045
Lu W, Shen X, Pavlova A, Lakkis M, Ward CJ, Pritchard L, Harris PC, Genest DR, Perez-Atayde AR, Zhou J: Comparison of Pkd1-targeted mutants reveals that loss of polycystin-1 causes cystogenesis and bone defects. Hum Mol Genet. 2001, 10 (21): 2385-2396. 10.1093/hmg/10.21.2385
Wettenhall JM, Smyth GK: limmaGUI: a graphical user interface for linear modeling of microarray data. Bioinformatics. 2004, 20 (18): 3705-3706. 10.1093/bioinformatics/bth449
Falcon S, Gentleman R: Using GOstats to test gene lists for GO term association. Bioinformatics. 2007, 23 (2): 257-258. 10.1093/bioinformatics/btl567
Chen WC, Tzeng YS, Li H: Gene expression in early and progression phases of autosomal dominant polycystic kidney disease. BMC Res Notes. 2008, 1: 131- 10.1186/1756-0500-1-131
Song X, Di Giovanni V, He N, Wang K, Ingram A, Rosenblum ND, Pei Y: Systems biology of autosomal dominant polycystic kidney disease (ADPKD): computational identification of gene expression pathways and integrated regulatory networks. Hum Mol Genet. 2009, 18 (13): 2328-2343. 10.1093/hmg/ddp165
Barrett T, Troup DB, Wilhite SE, Ledoux P, Rudnev D, Evangelista C, Kim IF, Soboleva A, Tomashevsky M, Marshall KA, et al.: NCBI GEO: archive for high-throughput functional genomic data. Nucleic Acids Res. 2009, 37 (Database issue): D885-890.
Lewis BP, Burge CB, Bartel DP: Conserved seed pairing, often flanked by adenosines, indicates that thousands of human genes are microRNA targets. Cell. 2005, 120 (1): 15-20. 10.1016/j.cell.2004.12.035
John B, Enright AJ, Aravin A, Tuschl T, Sander C, Marks DS: Human microRNA targets (vol 2, pg 1862, 2005). Plos Biology. 2005, 3 (7): 1328-1328.
Wang XW, El Naqa IM: Prediction of both conserved and nonconserved microRNA targets in animals. Bioinformatics. 2008, 24 (3): 325-332. 10.1093/bioinformatics/btm595
Maragkakis M, Reczko M, Simossis VA, Alexiou P, Papadopoulos GL, Dalamagas T, Giannopoulos G, Goumas G, Koukis E, Kourtis K, et al.: DIANA-microT web server: elucidating microRNA functions through target prediction. Nucleic Acids Research. 2009, 37: W273-W276. 10.1093/nar/gkp292
Sweeney W, Avner E: Functional activity of epidermal growth factor receptors in autosomal recessive polycystic kidney disease. Am J Physiol. 1998, 275: F387-F394.
Sweeney WE, Chen Y, Nakanishi K, Frost P, ED A: Treatment of polycystic kidney disease with a novel tyrosine kinase inhibitor. Kidney International. 2000, 57 (1): 33-40. 10.1046/j.1523-1755.2000.00829.x
Kuo NT, Norman JT, PD W: Acidic FGF regulation of hyperproliferation of fibroblasts in human autosomal dominant polycystic kidney disease. Biochem Mol Med. 1997, 61 (2): 178-191. 10.1006/bmme.1997.2583
Hollinger S, Hepler JR: Cellular regulation of RGS proteins: modulators and integrators of G protein signaling. Pharmacol Rev. 2002, 54 (3): 527-559. 10.1124/pr.54.3.527
Yamaguchi T, Hempson SJ, Reif GA, Hedge AM, Wallace DP: Calcium restores a normal proliferation phenotype in human polycystic kidney disease epithelial cells. J Am Soc Nephrol. 2006, 17 (1): 178-187. 10.1681/ASN.2006080912
Harris PC, Torres VE: Polycystic kidney disease. Annu Rev Med. 2009, 60: 321-337. 10.1146/annurev.med.60.101707.125712
Ong AC, Harris PC: Molecular pathogenesis of ADPKD: the polycystin complex gets complex. Kidney Int. 2005, 67 (4): 1234-1247. 10.1111/j.1523-1755.2005.00201.x
Serra AL, Poster D, Kistler AD, Krauer F, Raina S, Young J, Rentsch KM, Spanaus KS, Senn O, Kristanto P, et al.: Sirolimus and Kidney Growth in Autosomal Dominant Polycystic Kidney Disease. New England Journal of Medicine. 2010, 26 (9): 820-9. 363,
Walz G, Budde K, Mannaa M, Nürnberger J, Wanner C, Sommerer C, Kunzendorf U, Banas B, Hörl WH, Obermüller N, et al.: Everolimus in Patients with Autosomal Dominant Polycystic Kidney Disease. New England Journal of Medicine. 2010, 363 (9): 830-840. 10.1056/NEJMoa1003491
Petrie K, Guidez F, Howell L, Healy L, Waxman S, Greaves M, Zelent A: The histone deacetylase 9 gene encodes multiple protein isoforms. J Biol Chem. 2003, 278 (18): 16059-16072. 10.1074/jbc.M212935200
David D, Cardoso J, Marques B, Marques R, Silva ED, Santos H, Boavida MG: Molecular characterization of a familial translocation implicates disruption of HDAC9 and possible position effect on TGFbeta2 in the pathogenesis of Peters' anomaly. Genomics. 2003, 81 (5): 489-503. 10.1016/S0888-7543(03)00046-6
Chuaqui RF, Bonner RF, Best CJM, Gillespie JW, Flaig MJ, Hewitt SM, Phillips JL, Krizman DB, Tangrea MA, Ahram M, et al.: Post-analysis follow-up and validation of microarray experiments. Nature Genetics. 2002, 32: 509-514. 10.1038/ng1034
Morey JS, Ryan JC, Dolah FMV: Microarray validation: factors influencing correlation between oligonucleotide microarrays and real-time PCR. Biol Proced Online. 2006, 8 (1): 175-193. 10.1251/bpo126
Saal S, Harvey SJ: MicroRNAs and the kidney: coming of age. Curr Opin Nephrol Hypertens. 2009, 18 (4): 317-323. 10.1097/MNH.0b013e32832c9da2
Gaj S, Eijssen L, Mensink RP, Evelo CTA: Validating nutrient-related gene expression changes from microarrays using RT2 PCR-arrays. Genes and Nutrition. 2008, 3 (3-4): 153-157. 10.1007/s12263-008-0094-1
Marshall C: Specificity of receptor tyrosine kinase signaling: transient versus sustained extracellular signal-regulated kinase activation. Cell. 1995, 80: 179-185. 10.1016/0092-8674(95)90401-8
Piontek K, Menezes LF, Garcia-Gonzalez MA, Huso DL, Germino GG: A critical development switch defines the kinetics of kidney cyst formation after loss of Pkd1. Nature Medicine. 2007, 13 (12): 1490-1495. 10.1038/nm1675
Patel V, Li L, Cobo-Stark P, Shao X, Somlo S, Lin F, P I: Acute kidney injury and aberrant planar cell polarity induce cyst formation in mice lacking renal cilia. Hum Mol Genet. 2008, 17 (11): 1578-1590. 10.1093/hmg/ddn045
Happe H, Leonhard WN, Van der Wal A, Van de Water B, Lantinga-van Leeuwen IS, Breuning MH, De Heer E, DJ P: Toxic tubular injury in kidneys from Pkd1-deletion mice accelerates cystogenesis accompanied by dysregulated planar cell polarity and canonical Wnt signaling pathways. Hum Mol genet. 2009, 18 (14): 2532-2542. 10.1093/hmg/ddp190
Koupepidou P, Felekkis KN, Kranzlin B, Sticht C, Gretz N, C D: Cyst formation in the PKD2 (1-703) transgenic rat precedes deregulation of proliferation-related pathways. BMC Nephrology. 2010, 11 (23):
Saadi-Kheddouci S, Berrebi D, Romagnolo B, Cluzeaud F, Peuchmaur M, Kahn A, Vandewalle A, Perret C: Early development of polycystic kidney disease in transgenic mice expressing an activated mutant of the beta-catenin gene. Oncogene. 2001, 20 (42): 5972-5981. 10.1038/sj.onc.1204825
Torres VE, Harris PC: Autosomal dominant polycystic kidney disease: the last 3 years. Kidney International. 2009, 76: 149-168. 10.1038/ki.2009.128
Kawano Y, Kypta R: Secreted antagonists of the Wnt signalling pathway. Journal of Cell Sceince. 2003, 116: 2627-2634. 10.1242/jcs.00623.
Bafico A, Gazit A, Pramila T, Finch PW, Yaniv A, Aaronson SA: Interaction of Frizzled Related Protein (FRP) with Wnt Ligands and the Frizzled Receptor Suggests Alternative Mechanisms for FRP Inhibition of Wnt Signaling. Journal of Biological Chemistry. 1999, 274 (23): 16180-16187. 10.1074/jbc.274.23.16180
Belo JA, Bachiller D, Agius E, Kemp C, Borges AC, Marques S, Piccolo S, de Robertis EM: Cerberus-like is a secreted BMP and nodal antagonist not essential for mouse development. Genesis. 2000, 26: 265-270. 10.1002/(SICI)1526-968X(200004)26:4<265::AID-GENE80>3.0.CO;2-4
Torres VE, Harris PC: Mechanisms of Disease: autosomal dominant and recessive polycystic kidney diseases. Nature Clinical Practice Nephrology. 2007, 2 (1): 40-55.
Imbach T, Grunewald S, Schenk B, Burda P, Schollen E, Wevers RA, Jaeken J, de Klerk JB, Berger EG, Matthijs G, et al.: Multi-allelic origin of congenital disorder of glycosylation (CDG)-Ic. Hum Genet. 2000, 106 (5): 538-545. 10.1007/s004390050022
Charron AJ, Nakamura S, Bacallao R, Wandinger-Ness A: Compromised Cytoarchitecture and Polarized Trafficking in Autosomal Dominant Polycystic Kidney Disease Cells. Journal of Cell Biology. 2000, 149 (1): 111-124. 10.1083/jcb.149.1.111
Qin S, Taglienti M, Nauli SM, Contrino L, Takakura A, Zhou J, JA K: Failure to ubiquitinate c-Met leads to hyperactivation of mTOR signaling in a mouse model of autosomal dominant polycystic kidney disease. Journal of Clinical Investigation. 2010, 120 (10): 3617-3628. 10.1172/JCI41531
Reiss G, te Heesen S, Zimmerman J, Robbins PW, Aebi M: Isolation of the ALG6 locus of Saccharomyces cerevisiae required for glucosylation in the N-linked glycosylation pathway. Glycobiology. 1996, 6 (5): 493-498. 10.1093/glycob/6.5.493
Flynt AS, Lai EC: Biological principles of microRNA-mediated regulation: shared themes amid diversity. Nature Reviews Genetics. 2008, 9 (11): 831-842. 10.1038/nrg2455
Hooten NN, Abdelmohsen K, Gorospe M, Ejiogu N, Zonderman AB, Evans MK: microRNA Expression Patterns Reveal Differential Expression of Target Genes with Age. Plos One. 2010, 5 (5).
Inui M, Martello G, Piccolo S: MicroRNA control of signal transduction. Nature Reviews Molecular Cell Biology. 2010, 11 (4): 252-263.
Fish JE, Santoro MM, Morton SU, Yu SH, Yeh RF, Wythe JD, Lvey KN, Bruneau BG, Stainier DYR, Srivastava D: MiR-126 regulates angiogenic signaling and vascular integrity. Developmental Cell. 2008, 15 (2): 272-284. 10.1016/j.devcel.2008.07.008
Jackson RJ, N S: How do microRNAs regulate gene expression?. Sci STKE 2007. 2007, 2007 (367): re1-
Neilson JR, Zheng GX, Burge CB, Sharp PA: Dynamic regulation of miRNA expression in ordered stages of cellular development. Genes Dev. 2007, 21 (5): 578-589. 10.1101/gad.1522907
This work was supported by a National Institutes of Health/National Institutes of Diabetes and Digestive and Kidney Diseases grant to JAK. Drs. Pandey and Kreidberg acknowledge the support from P50 DK074030-01. SQ acknowledges support form the American Heart Association Founders Affiliate. JH was supported by Kidney Research Scientist Core Education and National Training Program (KRESCENT) postdoctoral fellowships and was a Fellow of the Pediatric Scientist Development Program (supported by the March of Dimes, Pediatric Chairs of Canada, and the SickKids Foundation). We thank Dr. Meredith Seamon for critically reading the manuscript. There is no conflict of interest for any of the authors of this paper.
PP designed and carried out the study, analysed and interpreted the data and drafted the manuscript. SQ participated in the experimental validation of the data. JH helped in animal experiments. JAK helped in designing the study, and participated in drafting the manuscript. All authors have given final approval of the version to be published.
Electronic supplementary material
Additional file 1:. Step-wise approaches used to predict target miRNAs for the significantly regulated genes obtained at three comparisons, shown in Additional file 1. We used four prediction tools- TargetScan, miRanda, miRDB and microT and the results obtained from individual tool were overlapped by custom written Perl scripts. Only those miRNA-mRNA pairs were considered for further study that were predicted by atleast two tools. (PDF 213 KB)
Additional file 2:Perl script used for integrating prediction tools results. Perl script used for integrating prediction tools results (PDF 45 KB)
Additional file 3:Gene-primers used for qPCR. This table shows list of primers used for qPCR analysis of genes (XLS 34 KB)
Additional file 4:. This figure shows the experimental design implied in this study. It shows four comparisons including two stages-E14.5 and E17.5 and two genotypes- Pkd1-/-and wild-type. (PDF 152 KB)
Additional file 5:Significantly regulated genes at E14.5. This table shows the significantly regulated genes obtained for mutant versus wild-type comparison at E14.5 (XLS 199 KB)
Additional file 6:Significantly regulated genes at E17.5. This table shows the significantly regulated genes obtained for mutant versus wild-type comparison at E17.5 (XLS 402 KB)
Additional file 7:Significantly regulated genes in wild-types across time points E14.5 to E17.5. This table shows the significantly regulated genes obtained for wild-types when compared at E14.5 and E17.5 (XLS 411 KB)
Additional file 8:Significantly regulated genes in mutants across time points E14.5 to E17.5. This table shows the significantly regulated genes obtained for mutant when compared at E14.5 and E17.5 (XLS 782 KB)
Additional file 9:Significantly regulated genes specific for disease transition. This table shows the significantly regulated genes that are specific for transition under disease condition. (XLS 660 KB)
Additional file 10:Significantly regulated genes specific for healthy transition. This table shows the significantly regulated genes that are specific for transition under healthy condition. (XLS 160 KB)
Additional file 11:Significantly regulated genes specific to aging. This table shows the significantly regulated genes that are specific for aging and common to both diseased and healthy conditions. (XLS 515 KB)
Additional file 12:Gene Ontology terms- overrepresented biological processes. This table shows the overrepresented Biological Processes from Gene ontology for significantly regulated genes. Sheets 1 through 3 show overrepresented biological processes for up- and down-regulated genes obtained for E14.5 mutant versus wild-type comparison, E17.5 mutant versus wild-type, and mutants at E14.5 versus E17.5. (XLS 200 KB)
Additional file 13:Gene Set Enrichment Analysis- overrepresented gene sets or pathways. This table shows the overrepresented gene-sets/pathways obtained from Gene Set Enrichment Analysis for significantly regulated genes. Sheets 1 through 3 show overrepresented gene-sets for up- and down-regulated genes obtained for E14.5 mutant versus wild-type comparison, E17.5 mutant versus wild-type, and mutants at E14.5 versus E17.5. (XLS 60 KB)
Additional file 14: Pkd1L3/L3and Pkd1-/-mouse models. This figure shows some of the early and late responsive genes we obtained in comparison between Pkd1L3/L3and Pkd1-/-mouse models. (PDF 292 KB)
Additional file 15: Pkd1L3/L3and Pkd1-/-mouse models. This table shows the common genes and pathways obtained from comparing the Pkd1L3/L3and Pkd1-/-mouse models. (XLS 165 KB)
Additional file 16: Pkd1-/-mouse model and human ADPKD. This table shows the common genes and pathways obtained from comparing the Pkd1-/-mouse model and human ADPKD. (XLS 444 KB)
Additional file 17:Gene Set Enrichment Analysis- overrepresented gene sets with shared Transcription Factors. This table shows the overrepresented overrepresented gene sets with shared Transcription Factors obtained from Gene Set Enrichment Analysis for significantly regulated genes. Sheets 1 through 3 show overrepresented gene-sets with shared Transcription Factors for up- and down-regulated genes obtained for E14.5 mutant versus wild-type comparison, E17.5 mutant versus wild-type, and mutants at E14.5 versus E17.5. (XLS 52 KB)
Additional file 18:Predicted target miRNAs for significantly regulated genes at E14.5. This table shows the predicted target miRNAs for significantly regulated genes at E14.5. Four tools-TargetScan, miRanda, miRDB and microT were used for prediction and results were overlapped. Only those miRNA:gene pairs are shown that are predicted by atleast two tools. (XLS 601 KB)
Additional file 19:Predicted target miRNAs for significantly regulated genes at E17.5. This table shows the predicted target miRNAs for significantly regulated genes at E17.5. Four tools-TargetScan, miRanda, miRDB and microT were used for prediction and results were overlapped. Only those miRNA:gene pairs are shown that are predicted by atleast two tools. (XLS 1 MB)
Additional file 20:Comparison of microarrays and qPCR results. This table shows the fold change obtained from both microarrays and qPCR for significantly regulated genes (that were validated by qPCR). (XLS 47 KB)
Additional file 21:miRNA:target inverse relationship. This table shows the predicted miRNA:mRNA inverse relationship for significantly regulated genes and 9 selected miRNAs. The expression of these nine miRNAs was validated by qPCR. Only those miRNA:mRNA pairs are shown that are predicted by atleast two tools and show inverse correlation to each other. (XLS 48 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Pandey, P., Qin, S., Ho, J. et al. Systems biology approach to identify transcriptome reprogramming and candidate microRNA targets during the progression of polycystic kidney disease. BMC Syst Biol 5, 56 (2011). https://doi.org/10.1186/1752-0509-5-56
- Gene Ontology
- Cyst Formation
- Autosomal Dominant Polycystic Kidney Disease
- Cystic Kidney
- GPCR Signaling