Microarray-based approach identifies microRNAs and their target functional patterns in polycystic kidney disease
© Pandey et al; licensee BioMed Central Ltd. 2008
Received: 11 September 2008
Accepted: 23 December 2008
Published: 23 December 2008
MicroRNAs (miRNAs) play key roles in mammalian gene expression and several cellular processes, including differentiation, development, apoptosis and cancer pathomechanisms. Recently the biological importance of primary cilia has been recognized in a number of human genetic diseases. Numerous disorders are related to cilia dysfunction, including polycystic kidney disease (PKD). Although involvement of certain genes and transcriptional networks in PKD development has been shown, not much is known how they are regulated molecularly.
Given the emerging role of miRNAs in gene expression, we explored the possibilities of miRNA-based regulations in PKD. Here, we analyzed the simultaneous expression changes of miRNAs and mRNAs by microarrays. 935 genes, classified into 24 functional categories, were differentially regulated between PKD and control animals. In parallel, 30 miRNAs were differentially regulated in PKD rats: our results suggest that several miRNAs might be involved in regulating genetic switches in PKD. Furthermore, we describe some newly detected miRNAs, miR-31 and miR-217, in the kidney which have not been reported previously. We determine functionally related gene sets, or pathways to reveal the functional correlation between differentially expressed mRNAs and miRNAs.
We find that the functional patterns of predicted miRNA targets and differentially expressed mRNAs are similar. Our results suggest an important role of miRNAs in specific pathways underlying PKD.
MicroRNAs (miRNAs) are known to regulate the expression of key genes relevant to cancer and potentially other diseases [1–3]. MiRNAs are short noncoding RNAs of about 22 nt length that have recently been shown to play important roles in mammalian gene expression [1, 3–5]. They induce posttranscriptional gene repression by blocking protein translation (by binding to the 3' UTR of their target genes) or by inducing mRNA degradation, and have the potential to play central roles in physiological and pathological conditions. Recently, it has been shown that miRNAs can also increase translation [6, 7]. The physiological conditions of a cell seem to affect the recruitment of regulatory proteins, which can alter the effect of a miRNA.
MiRNAs are transcribed as long primary transcripts (pri-miRNAs), some of them being polycistronic, which are processed in the cell nucleus by an enzyme called Drosha, yielding precursor miRNAs (pre-miRNAs) that exhibit a characteristic stem-loop sequence . These are exported into the cytosol where mature miRNAs are generated by the RNAse Dicer, producing the small single-stranded miRNA . Translational inhibition, which seems to be the major mode of action in animals, is performed by a riboprotein complex called RNA-induced silencing complex (RISC) consisting of the miRNA and proteins of the argonaute family [9, 10].
MiRNAs are involved in several cellular processes, including cellular differentiation [11, 12], organism development [13, 14], and apoptosis [15, 16]. While all of these are conserved in metazoans, the number of conserved miRNAs between mammals suggests that there are additional functions only found in vertebrates , e.g. controlling hematopoietic differentiation . Recent studies provide growing evidence for the involvement of miRNAs in cancer pathomechanisms [18–21]. However, to date nothing is known regarding miRNAs in the context of Polycystic Kidney Diseases.
Cilia and flagella are ancient, evolutionary conserved organelles that project from cell surfaces to perform diverse biological roles, including whole-cell locomotion, movement of fluid, chemo-, mechano-, and photosensation, and sexual reproduction. The concept of ciliopathies has helped in advancing a unifying theory of cystic kidney diseases . This theory states that the products of all genes that are mutated in cystic kidney diseases in humans, mice, or zebrafish are expressed in primary cilia or centrosomes of renal epithelial cells .
There are numerous disorders linked to basal body and/or cilia dysfunction, including polycystic kidney disease (PKD), primary ciliary dyskinesia (PCD), nephronophthisis (NPHP1–9) , Senior-Loken syndrome, Joubert syndrome, Meckel syndrome, oral-facial-digital syndrome, Alström syndrome and Bardet-Biedl syndrome . These syndromes are typically associated with one or more of the symptoms like cystic kidneys , retinal degeneration and retinitis pigmentosa , situs inversus [26, 27], anosmia , respiratory problems , infertility , hydrocephalus , other ailments like obesity, diabetes, liver fibrosis, hypertension, heart malformations, skeletal anomalies (e.g. polydactyly), cognitive impairment and developmental defects such as exencephaly [16, 26].
Although the mechanisms of the cyst formation are not clearly understood, they are postulated to involve improper functioning of several pathways including cell proliferation, apoptosis, cell polarity, and fluid secretion . Woo  has described apoptotic cells in glomeruli, cyst walls, and in both cystic and non-cystic tubules of the polycystic kidneys. Apoptotic loss of renal tissue may be associated with the progressive deterioration of renal function that occurs in patients with autosomal dominant polycystic kidney disease (ADPKD). There is evidence that genes involved in the regulation of cell proliferation, such as p53, c-fos, cyclin D1, and c-myc may be involved in the control of apoptosis . Veis and colleagues have shown overexpression of c-myc in human ADPKD in association with increased levels of apoptosis and cell proliferation . It is clear that pathogenesis of PKD is very complicated and involves multiple molecular pathways with overlapping, complementary, or opposing effects. There are several signalling pathways that have been implicated in ciliary function . The dysregulation of mitogen-activated protein kinases (MAPKs) in the cyst epithelium of pcy (polycystic kidney disease) mice, carrying a missense mutation in NPHP3 , is a downstream consequence of disturbed renal monocilia function . The proteins, implicated in the formation of renal cysts in tuberous sclerosis, have been found at the ciliary base. They form a complex that inhibits mTOR resulting in retarded cyst formation in rats with PKD . Cilia-mediated signalling acts as a switch between canonical and non-canonical Wnt signalling pathways .
Rats or mice have been used as common model systems for the study of PKD. The Hannover rat, Han:Sprague-Dawley (SPRD)-cy rat is an accepted model for human PKD [38, 39] and has been efficiently used since more than a decade. It is an autosomal dominant model for PKD resulting in cyst formation and slowly progressive chronic renal failure . In the current investigation, we explore the transcriptional changes that occur in PKD and investigate if these changes could be related to miRNAs. We use a microarray-based approach to profile the transcriptional (mRNA) changes as well as changes in the miRNA expression patterns in PKD. We use the Han:SPRD cy/+ rat model  for our current investigation. Our results suggest several miRNAs may be involved in regulating the genetic switches in PKD. Furthermore we describe some newly detected miRNAs in the kidney.
Animals and physiological state
Inbred homozygous unaffected (+/+) and heterozygous affected (cy/+) Han:SPRD rats exhibiting PKD were investigated. The Han:SPRD-cy rats used in this study carry a dominant mutation that causes cystic kidneys and is an accepted model for human PKD . After approximately 40 generations of inbreeding this substrain was registered as PKD/Mhm-cy inbred strain of rats: polycystic kidney diseases, Mannheim, Germany hereafter designated as PKD/Mhm. From each of the aforementioned PKD/Mhm (cy/+) and PKD/Mhm (+/+) animals, littermates were investigated. The animals were sacrificed by cervical dislocation at day 36. On the day of sacrifice, body weight was determined. Following dislocation, the left and right kidneys were immediately removed and weights were determined. Kidneys were preserved for histological analyses and the genotypes were confirmed. Histological analysis of pathogenesis was performed using hematoxylin-esosin (HE) staining of the kidney section (3 μm) followed by cyst grading under light microscope [41, 42]; kidneys were graded on 1–5 scale as previously described [41, 42].
We had ethical approval to carry out this work on animals by Regierungspraesidium Nordbaden and Internal Review Board.
RNA isolation and Affymetrix microarray
Total RNA was extracted using TRIzol method according to manufacturer's protocol (Invitrogen Life Technologies). cDNA synthesis was performed using the SuperScript Choice System (Invitrogen Life Technologies, Invitrogen Corporation) according to manufacturer's protocol. Biotin-labelled cRNA was produced using ENZO BioArray High Yield RNA Transcript Labelling Kit. The standard protocol from Affymetrix (Santa Clara, CA) with 3.3 μL of cDNA was used for the in vitro transcription (IVT). Cleanup of the IVT product was done using CHROMA SPIN-100 columns (Clontech, USA). Spectrophotometric analysis was used for quantification of cRNA with acceptable A260/A280 ratio of 1.9 to 2.1. After that, the cRNA was fragmented using Affymetrix defined protocol. Labelled and fragmented cRNA was hybridized to Affymetrix Rat 230_2 microarrays for 16 hrs at 45°C using Affymetrix defined protocol. cRNA in the range of 15 μg was used for all the 12 microarrays. Microarrays were washed using an Affymetrix fluidics station 450 and stained initially with streptavidin/phycoerytherin. For each sample the signal was further enhanced by incubation with biotinylated goat anti-streptavidin followed by a second incubation with streptavidin/phycoerytherin, and a second round of intensities were measured. Microarrays were scanned with an Affymetrix scanner controlled by the Affymetrix Microarray Suite software. A total of 12 Affymetrix whole genome arrays (from 6 healthy and 6 diseased biological replicates) were performed.
Locked Nucleic Acid (LNA) based miRNA hybridization assays were done as described by Castoldi et al. . Briefly, 5 μg RNA was hybridized to a miRNA microarray (miChip v6.0), containing probes for ~300 miRNAs, based on LNA-modified and Tm-normalized oligonucleotide capture probes (miRCURY Array probes, Exiqon) spotted onto Codelink (GE Healthcare) in multiple replicates, and scanned using an Axon Scanner 4000B. Microarray images were analyzed using the Genepix Pro 4.0 software . As for the Affymetrix whole genome arrays, 6 biologically replicated arrays, each for healthy and diseased animals (total of 12 chips) were performed.
Statistical and bioinformatics analysis of microarray data
Microarray data obtained from Affymetrix chips were analyzed using SAS Micro-Array Solution version 1.3 (SAS, USA). Normalization of raw data was performed by fitting a mixed linear model across all arrays in the experiment. Log2-transformed scores of all spot measures were subjected to normalization. Identification of differentially expressed genes was carried out by Mixed Model Analysis (MMA) using ANOVA approach. A Bonferroni adjustment for multiple testing with α ≤ 0.05 was used to calculate the statistical significance threshold/cut-off [negative log10 (α)]. Custom CDFs version 8 (UniGene CDF files)  were used for annotations. Significantly up-regulated genes in PKD/Mhm (cy/+) animals were defined as those with >0 log-fold-changes in expression and down-regulated as those having <0 log-fold-changes in expression at p-value < 0.05 (adjusted for multiple testing).
The raw data from miRNA chips were also normalized and analyzed using the same modules in SAS Micro-Array Solution version 1.3, as described above. Significantly up-regulated miRNAs in PKD/Mhm (cy/+) animals were defined as those with >0 log-fold-changes in expression and down-regulated as those having <0 log-fold-changes in expression at p-value < 0.05 (adjusted for multiple testing).
Genes previously reported as targets for the significantly regulated miRNAs, and differentially regulated in PKD, were obtained from Argonaute . In addition, targets for the significantly regulated miRNAs, in the genes differentially regulated in PKD, were predicted using two algorithms: TargetScan  and miRanda . To identify the genes commonly predicted by both algorithms, results were intersected using a Perl script; the intersection of TargetScan and miRanda algorithms was used for further analysis.
Identification of pathways possibly affected by miRNAs and their target genes in PKD
The pathways affected in PKD were determined for the significantly regulated genes from KEGG , GeneOntology (GO) , Biocarta http://www.biocarta.com/ and the Molecular Signature database . To check whether significantly regulated genes are overrepresented in a pathway in PKD or not, Fisher's exact test  without correction for multiple testing was used. Pathways showing p-value less than 0.05 were considered as significantly enriched.
Reverse transcription and quantitative real-time PCR (qPCR) analysis of miRNAs
qPCR assays for miRNAs were performed using "TaqMan MicroRNA Assays" (Applied Biosystems, USA). The reverse transcriptase reactions were performed as per the manufacturer's protocol. Real-time PCR was also performed using a standard TaqMan PCR kit protocol on the Applied Biosystems 7000 Sequence Detection System. 6 biological replicates were used for analysis and all the reactions were run in quadruplicates. The comparative CT method for relative quantitation of gene expression was used to determine miRNA expression levels among the normal animals and PKD/Mhm (cy/+) animals according to manufacturer's protocol. miR-193a was used as an endogenous control to normalize the expression levels of targets. The miR-193a served as a good choice for endogenous control because its expression was almost uniform in all the samples on the miRNA-chips (0.99 ± 0.026) and it was not differentially regulated among the tested samples. miRNAs with ubiquitous and stable expression values are superior normalizers to other RNAs such as 5S rRNA, U6 snRNA or total RNA [52, 53]. Consistency in expression of miR-193a across all the samples was evident in qPCR assays (mean CT-value 29.51 ± 0.19). The significance of differences in relative expression of miRNAs among the two groups was tested by One-way ANOVA method in SAS version 9.1.
Pathological status of PKD/Mhm(cy/+) and healthy PKD/Mhm (+/+) animals
Profiling miRNA expression change during PKD
Significantly regulated miRNAs on Exiqon chips from SAS analysis
Name of miRNAs on chip
miRNA name in rat*
Lsmean animal PKD
Lsmean animal Control
Log Fold Change PKD_Ctrl
-log10(p-value) for log fold change of PKD_Ctrl
Microarray analysis of genes involved in PKD
Top 25 up-regulated and 25 down-regulated genes on Affymetrix chips from SAS analysis
Log_fold_change (>0 = up-regulated; <0 = down-regulated)
-log10(p-value) for Estimate of diseased_healthy
1368821_at, 1368822_at, 1371331_at, 1394119_at
1389189_at, 1396539_at, 1398294_at
1367823_at, 1375144_at, 1386940_at, 1388312_at
1369400_a_at, 1388044_at, 1388063_a_at, 1398320_at
1367571_a_at, 1371206_a_at, 1398322_at
Identifying miRNA-target interaction
miRNAs and their targets
a: miRNAs and their previously reported targets (from Argonaute database). miRNAs and their corresponding targets are both differentially regulated during PKD.
ABCB9, BTG1, DSCR1, RASD1
ABC transporters General
Adherens junction, Maturity onset diabetes of the, Focal adhesion, **Long term depression
ABCB9, CPB2, IRS1, MAPK10
ABC transporters General, Complement and coagulation cas, Adipocytokine signaling pathwa, Insulin signaling pathway, Type II diabetes mellitus, Fc epsilon RI signaling pathwa, Focal adhesion, **GnRH signaling pathway, **MAPK signaling pathway, Toll like receptor signaling p, Wnt signaling pathway
BTG1, SARA1, YWHAB
**Long term depression
Adherens junction, Axon guidance, Focal adhesion, Leukocyte transendothelial mig, Regulation of actin cytoskelet, TGF beta signaling pathway, T cell receptor signaling path, Tight junction, Wnt signaling pathway
ATP2B2, DNM1L, EGR3, PPP1R9B, YWHAB
**Calcium signaling pathway, Cell cycle
HRH3, NCDN, SLC23A2
**Neuroactive ligand receptor in
b: miRNAs and their targets (from TargetScan and miRanda). miRNAs and their corresponding targets are both differentially regulated during PKD.
ABCB9, BTG2, CACNB2, CNR1, COL3A1, GLRA2, IRS1, NEK6, PDE7B, PTPN5, SV2A, SYT4, YWHAB, NR5A2, NTRK3
ABC transporters General, **MAPK signaling pathway, **Neuroactive ligand receptor in, **Cell Communication, **ECM receptor interaction, Focal adhesion, Adipocytokine signaling pathwa, Insulin signaling pathway, Type II diabetes mellitus, Purine metabolism, Cell cycle, Maturity onset diabetes of the
ADAMTS1, ATP2B2, CACNB2, CDH13, CNR1, DUSP6, EGR3, EPHA7, GRIK2, GRM5, GRM7, HOXA1, MMP14
**Calcium signaling pathway, **MAPK signaling pathway, **Neuroactive ligand receptor in, Axon guidance, Gap junction, **Long term depression, Long term potentiation, **GnRH signaling pathway
ABCB9, COL3A1, EPHA7, GAS7, OTX1
ABC transporters General, **Cell Communication, **ECM receptor interaction, Focal adhesion, Axon guidance
Cell cycle, **MAPK signaling pathway, TGF beta signaling pathway
CLU, FN1, GRIK2, KCNH5, NR4A2, RAP1B
**Cell Communication, **ECM receptor interaction, Focal adhesion, Regulation of actin cytoskelet, **Neuroactive ligand receptor in, Leukocyte transendothelial mig, Long term potentiation, **MAPK signaling pathway
**Calcium signaling pathway
ATP2B2, CNR1, SNF1LK
**Calcium signaling pathway, **Neuroactive ligand receptor in
**Calcium signaling pathway
CCKBR, COLEC12, EPHA7, GRM7, JUN, KCNAB1, LIN7A, MAMDC1, PRRX1
**Calcium signaling pathway, **Neuroactive ligand receptor in, Axon guidance, B cell receptor signaling path, Focal adhesion, **GnRH signaling pathway, **MAPK signaling pathway, T cell receptor signaling path, Toll like receptor signaling p, Wnt signaling pathway
**Neuroactive ligand receptor in
ATP2B2, NR4A2, SNF1LK, RET, SLC11A2
**Calcium signaling pathway
COLEC12, DNM3, GRIK2, IGF1R, IL12B, KCNH7, KCNIP4, SCN3A, SCN5A
**Neuroactive ligand receptor in, Adherens junction, Focal adhesion, **Long term depression, **Cytokine cytokine receptor int, **Jak STAT signaling pathway, Toll like receptor signaling p, **Type I diabetes mellitus
BTG2, CALCR, DLL1, MYRIP
**Neuroactive ligand receptor in, Notch signaling pathway
Amyotrophic lateral sclerosis, Apoptosis, **Jak STAT signaling pathway, Neurodegenerative Disorders
Secondly, we mapped the differentially regulated genes for miRNA targets for the differentially regulated miRNAs using two tools, TargetScan  and miRanda . Only commonly predicted targets for different miRNAs were taken for further analysis. A number of 65 genes were identified as miRNA targets, common from TargetScan as well as miRanda (Table 3b), out of which 31 genes participate in regulating 33 pathways, which could further be grouped into 13 functional categories (Figure 6b). The intersection of TargetScan and miRanda results gave target genes for 20 miRNAs out of 30 differentially regulated miRNAs. This shows that available prediction tools are not yet fully optimal. Altogether, we obtained 89 target genes differentially expressed on Affymetrix chips (Figure 2). These were associated with 36 pathways, out of which 10 were significantly enriched by ORA (Table 3a and 3b; Figure 2).
Selected inverse miRNA-target relation identified
Log fold change miRNA
Log fold change gene
**Cell Communication, **ECM receptor interaction, Focal adhesion
**Cell Communication, **ECM receptor interaction, Focal adhesion
Adherens junction, Focal adhesion, **Long term depression
**GnRH signaling pathway
**MAPK signaling pathway
**Neuroactive ligand receptor in
**Cell Communication, **ECM receptor interaction, Focal adhesion
Cell cycle, **MAPK signaling pathway, TGF beta signaling pathway
Adherens junction, Axon guidance, Focal adhesion, Leukocyte transendothelial mig, Regulation of actin cytoskelet, T cell receptor signaling path, TGF beta signaling pathway, Tight junction, Wnt signaling pathway
**Cell Communication, **ECM receptor interaction, Focal adhesion, Regulation of actin cytoskelet
Focal adhesion, Leukocyte transendothelial mig, Long term potentiation, **MAPK signaling pathway
**Neuroactive ligand receptor in
B cell receptor signaling path, Focal adhesion, **GnRH signaling pathway, **MAPK signaling pathway, T cell receptor signaling path, Toll like receptor signaling p, Wnt signaling pathway
Adherens junction, Focal adhesion, **Long term depression
**Neuroactive ligand receptor in
In the current investigation, we explore the possible involvement of miRNAs in PKD. In a well characterized rat model system (Han:SPRD cy/+ rat model; ), we used a combinatorial approach involving data mining and microarray analysis, to profile the miRNAs involved in PKD. In parallel to the miRNA expression profiling, we also describe the changes in mRNA transcript patterns in PKD using genome-wide Affymetrix arrays. Whereas, the gene expression arrays e.g. Affymetrix arrays have already established themselves for measuring large-scale differential regulation of mRNAs, microarrays have recently been developed to measure miRNA expression changes accurately [43, 62]. We used arrays based on LNA technology . The cluster analysis of the control and PKD miRNA-microarrays showed tight clustering PKD samples (Figure 3). Few studies have appeared showing the use of proteomics based approaches for identifying specific miRNA targets [63–66]. Several of these studies also applied mRNA-expression arrays to derive correlations between the protein- and messenger- turndown. Over-expressing miRNAs resulted in repressed targets, both at protein and messenger levels: how much both processes contribute to down-regulation of targets depends on individual miRNA-mRNA pairs . These studies show that mRNA profiling are indeed valuable tools for studying mRNA-miRNA interaction, although additional information on fine-tuning of proteome-expression may be obtained by including proteomics approaches.
Our study, profiling the changes in the mRNA transcripts in PKD, revealed large-scale changes. Although some transcription factors, such as SP1 , JUN and c-myc  have been implicated to play a role in PKD, the magnitude of transcriptional changes (>900 genes) suggest involvement of other regulators. Moreover, the regulatory basis of changes in the expression of transcription factors involved in PKD remains poorly understood. On the other hand, miRNAs have emerged recently as key regulators of gene expression in many developmental [69, 70] and disease events [71, 72]. Our miRNA study shows that changes in gene expression during PKD also involve miRNAs.
When we compared the mRNA and miRNA profiles, differentially regulated in PKD, with Argonaute (a comprehensive database on miRNAs; [45, 71]), there were few genes reported as miRNA target like tropomyosin 1, alpha (TPM1) as a target of miR-21, the beta polypeptide of tyrosine 3-monooxygenase/tryptophan 5-monooxygenase activation protein (YWHAB), regulatory subunit 9B of protein phosphatase 1 (PPP1R9B), early growth response 3 (EGR3) and dynamin 1-like (DNM1L) as targets of miR-31, plysia ras-related homolog A2 (RHOA) as targets of miR-217, etc. (Table 3a). But none of these genes had previously been associated with PKD. Over-representation analysis showed the enrichment of five pathways namely, GnRH signaling pathway, MAPK signaling pathway, Long term depression, Calcium signaling pathway and Neuroactive ligand receptor interaction. The comparison of differentially regulated genes/miRNAs to Argonaute database indicates that previously known miRNA-target interactions (Table 3a) could play important roles in PKD.
Parallel profiling of the transcripts as well as miRNAs on the same set of samples gives us insight into potential interactions between miRNAs and differentially expressed targets in PKD. This parallel profiling allowed us to scrutinize for the negative expression patterns for the miRNA-target relation at a given time point (Table 4). A negative relationship of the expression patterns between miRNA and its target mRNA is an important parameter for determining their interactions because miRNA, in general, are regarded as negative regulators of their targets.
Microarray analysis of miRNAs revealed down regulation of 29 out of 30 differentially regulated miRNAs, which corresponded to increased expression of several genes related to pathways upregulated during PKD. To further uncover the functional correlation between differentially expressed mRNAs and miRNAs, we determined functionally related gene sets, or pathways. The differentially regulated mRNAs were associated with 24 functional categories, which included several enriched pathways important to renal diseases. Previous studies of cystic kidneys implicated several pathways thought to contribute to the pathogenesis of renal cysts' formation like mTOR signalling, MAPK (Mitogen-Activated Protein Kinase) signalling, Wnt signalling, and the TGF-β (Transforming Growth Factor-β) pathway.
MAPKs play important roles in the cell by transmitting extracellular signals from the cell membrane to the nucleus . MAPKs are activated by various stimuli, influencing cell proliferation, differentiation, and apoptosis. Aberrant regulation of MAPKs and other signalling pathways has been reported to be consistent with altered regulation of cell proliferation and differentiation observed in renal cystic disease . Sustained activation of MAPKs in kidney epithelial cells inhibits normal epithelial phenotype and formation of adherens junctions . The expression of genes involved in MAPK signalling pathway like MAP3K1, JUN, and MYC was significantly higher in the PKD animals, indicating the activation of MAPK signalling pathway.
Wnt signalling is essential for renal development. Recent research has revealed an unexpected intersection between Wnt signalling and PKD . It has been reported that canonical Wnt signalling seemed mandatory for early renal development, but persistent β-catenin signalling seemed to trigger cyst formation at later developmental stages . Components of the Wnt signalling pathway like FZD2, JUN, MYC, and RHOA were significantly upregulated in our PKD animals.
TGF-β promotes renal cell hypertrophy and stimulates extracellular matrix accumulation in several renal diseases, including diabetic nephropathy  and PKD . It activates the inhibitors of the proteases e.g. tissue inhibitors of metalloproteinases and plasminogen activator inhibitor 1 . TGF-β 1, TGF-β 2 and TGF-β 3 were consistently upregulated in PKD-affected animals. The expression of extracellular matrix components, such as collagen triple helix repeat containing 1 (CTHRC1), and fibronectin 1 (FN1) was significantly higher in animals with PKD. One of the integrins, integrin beta 1 (ITGB1) was upregulated. Therefore, there was a clear activation of the TGF-β signaling pathway in the PKD animals with a significant increase of the synthesis of the extracellular matrix components and inhibition of the proteases that digest this matrix. Recently Kato et al. showed that miRNA miR-192 has been associated with TGF-β pathway in diabetic kidney glomeruli . Its implication in PKD has not been shown, though.
We could not find negative co-relations between the predicted targets of the only up-regulated miRNA, miR-21, whose expression was verified by qPCR analysis too. All the predicted targets were up-regulated. It should be noted that alternate action mechanisms for miRNAs also exist. The action of miRNAs may not to be reflected on the level of their target mRNAs as they are believed to block or attenuate further translation of mRNAs to protein. In such conditions, miRNAs will exert their regulatory role on the level of translation. Moreover, a new mode of action where the miRNAs act as positive regulators has also been defined recently . Depending upon the state of a cell, a miRNA can act as positive or negative regulator .
The specific cellular pathways that were found to be associated with dysregulated miRNAs or differentially expressed transcripts in PKD can be used to shape some initial hypotheses on how alteration of miRNA expression may be directly involved in the disease. Furthermore, the functional correlation between the differentially expressed mRNAs and miRNAs in PKD revealed a tight posttranscriptional regulation network at both mRNA and protein level.
While this work establishes a multi-tiered approach for investigation of miRNA mode of genetic regulations of PKD, the results await further investigation. PKD is an important disease as 8–10% of all patients on renal replacement therapy including haemodialysis or transplantation, suffer from PKD . The molecular components have just started to become apparent, and we add another layer of regulators, namely miRNAs. We predict that several of the differentially regulated genes are miRNA targets and miR-21, miR-31, miR-128, miR-147 and miR-217 may be the important players in such interaction. It is interesting to note that miR-31 and miR-217 have not been previously reported in kidney. miR-21 has been reported in kidney  but its function has yet not been established. It has been speculated in previous studies that a single miRNA could target more than hundred genes and one gene could be the target of several miRNAs [61, 80]. Knockout and over-expression studies will provide further insight into the regulatory interaction between these miRNAs and their targets in order to properly understand PKD development and design new therapeutic measures.
We thank Dr. Xiaolei Yu for curated pathway lists, Dr. Mirco Castoldi for performing Exiqon miRNA arrays at EMBL, Heidelberg, Germany and Maria Saile for the technical support. This work was supported by the DFG research Training Group 886 'Molecular Imaging Methods for the analysis of Gene and Protein Expression'. There is no kind of conflict of interest for any of the authors of this paper.
- He L, Hannon GJ: MicroRNAs: small RNAs with a big role in gene regulation. Nat Rev Genet. 2004, 5: 522-531. 10.1038/nrg1379.PubMedView Article
- Pasquinelli AE, Hunter S, Bracht J: MicroRNAs: a developing story. Curr Opin Genet Dev. 2005, 15: 200-205. 10.1016/j.gde.2005.01.002.PubMedView Article
- Esquela-Kerscher A, Slack FJ: Oncomirs – microRNAs with a role in cancer. Nat Rev Cancer. 2006, 6: 259-269. 10.1038/nrc1840.PubMedView Article
- Ambros V: The functions of animal microRNAs. Nature. 2004, 431: 350-355. 10.1038/nature02871.PubMedView Article
- Pasquinelli AE: Demystifying small RNA pathways. Dev Cell. 2006, 10: 419-424. 10.1016/j.devcel.2006.03.005.PubMedView Article
- Buchan JR, Parker R: Molecular biology. The two faces of miRNA. Science. 2007, 318: 1877-1878. 10.1126/science.1152623.PubMedView Article
- Vasudevan S, Tong Y, Steitz JA: Switching from repression to activation: microRNAs can up-regulate translation. Science. 2007, 318: 1931-1934. 10.1126/science.1149460.PubMedView Article
- Cullen BR: Transcription and processing of human microRNA precursors. Mol Cell. 2004, 16: 861-865. 10.1016/j.molcel.2004.12.002.PubMedView Article
- Tang G: siRNA and miRNA: an insight into RISCs. Trends Biochem Sci. 2005, 30: 106-114. 10.1016/j.tibs.2004.12.007.PubMedView Article
- Tijsterman M, Plasterk RH: Dicers at RISC; the mechanism of RNAi. Cell. 2004, 117: 1-3. 10.1016/S0092-8674(04)00293-4.PubMedView Article
- Chang S, Johnston RJ, Frokjaer-Jensen C, Lockery S, Hobert O: MicroRNAs act sequentially and asymmetrically to control chemosensory laterality in the nematode. Nature. 2004, 430: 785-789. 10.1038/nature02752.PubMedView Article
- Johnston RJ, Hobert O: A microRNA controlling left/right neuronal asymmetry in Caenorhabditis elegans. Nature. 2003, 426: 845-849. 10.1038/nature02255.PubMedView Article
- Lee RC, Feinbaum RL, Ambros V: The C. elegans heterochronic gene lin-4 encodes small RNAs with antisense complementarity to lin-14. Cell. 1993, 75: 843-854. 10.1016/0092-8674(93)90529-Y.PubMedView Article
- Reinhart BJ, Slack FJ, Basson M, Pasquinelli AE, Bettinger JC, Rougvie AE, Horvitz HR, Ruvkun G: The 21-nucleotide let-7 RNA regulates developmental timing in Caenorhabditis elegans. Nature. 2000, 403: 901-906. 10.1038/35002607.PubMedView Article
- Brennecke J, Hipfner DR, Stark A, Russell RB, Cohen SM: bantam encodes a developmentally regulated microRNA that controls cell proliferation and regulates the proapoptotic gene hid in Drosophila. Cell. 2003, 113: 25-36. 10.1016/S0092-8674(03)00231-9.PubMedView Article
- Xu P, Vernooy SY, Guo M, Hay BA: The Drosophila microRNA Mir-14 suppresses cell death and is required for normal fat metabolism. Curr Biol. 2003, 13: 790-795. 10.1016/S0960-9822(03)00250-1.PubMedView Article
- Chen CZ, Li L, Lodish HF, Bartel DP: MicroRNAs modulate hematopoietic lineage differentiation. Science. 2004, 303: 83-86. 10.1126/science.1091903.PubMedView Article
- Gregory RI, Shiekhattar R: MicroRNA biogenesis and cancer. Cancer Res. 2005, 65: 3509-3512. 10.1158/0008-5472.CAN-05-0298.PubMedView Article
- He L, Thomson JM, Hemann MT, Hernando-Monge E, Mu D, Goodson S, Powers S, Cordon-Cardo C, Lowe SW, Hannon GJ, Hammond SM: A microRNA polycistron as a potential human oncogene. Nature. 2005, 435: 828-833. 10.1038/nature03552.PubMedView Article
- Lu J, Getz G, Miska EA, varez-Saavedra E, Lamb J, Peck D, Sweet-Cordero A, Ebert BL, Mak RH, Ferrando AA, Downing JR, Jacks T, Horvitz HR, Golub TR: MicroRNA expression profiles classify human cancers. Nature. 2005, 435: 834-838. 10.1038/nature03702.PubMedView Article
- McManus MT: MicroRNAs and cancer. Semin Cancer Biol. 2003, 13: 253-258. 10.1016/S1044-579X(03)00038-5.PubMedView Article
- Hildebrandt F, Zhou W: Nephronophthisis-associated ciliopathies. J Am Soc Nephrol. 2007, 18: 1855-1871. 10.1681/ASN.2006121344.PubMedView Article
- Otto EA, Trapp ML, Schultheiss UT, Helou J, Quarmby LM, Hildebrandt F: NEK8 mutations affect ciliary and centrosomal localization and may cause nephronophthisis. J Am Soc Nephrol. 2008, 19: 587-592. 10.1681/ASN.2007040490.PubMedView Article
- Inglis PN, Boroevich KA, Leroux MR: Piecing together a ciliome. Trends Genet. 2006, 22: 491-500. 10.1016/j.tig.2006.07.006.PubMedView Article
- Hildebrandt F, Otto E: Cilia and centrosomes: a unifying pathogenic concept for cystic kidney disease?. Nat Rev Genet. 2005, 6: 928-940.PubMedView Article
- Afzelius BA: Cilia-related diseases. J Pathol. 2004, 204: 470-477. 10.1002/path.1652.PubMedView Article
- Davenport JR, Yoder BK: An incredible decade for the primary cilium: a look at a once-forgotten organelle. Am J Physiol Renal Physiol. 2005, 289: F1159-F1169. 10.1152/ajprenal.00118.2005.PubMedView Article
- Kulaga HM, Leitch CC, Eichers ER, Badano JL, Lesemann A, Hoskins BE, Lupski JR, Beales PL, Reed RR, Katsanis N: Loss of BBS proteins causes anosmia in humans and defects in olfactory cilia structure and function in the mouse. Nat Genet. 2004, 36: 994-998. 10.1038/ng1418.PubMedView Article
- Bouillet P, Robati M, Bath M, Strasser A: Polycystic kidney disease prevented by transgenic RNA interference. Cell Death Differ. 2005, 12: 831-833. 10.1038/sj.cdd.4401603.PubMedView Article
- Igarashi P, Somlo S: Genetics and pathogenesis of polycystic kidney disease. J Am Soc Nephrol. 2002, 13: 2384-2398. 10.1097/01.ASN.0000028643.17901.42.PubMedView Article
- Woo D: Apoptosis and loss of renal tissue in polycystic kidney diseases. N Engl J Med. 1995, 333: 18-25. 10.1056/NEJM199507063330104.PubMedView Article
- Steller H: Mechanisms and genes of cellular suicide. Science. 1995, 267: 1445-1449. 10.1126/science.7878463.PubMedView Article
- Veis DJ, Sorenson CM, Shutter JR, Korsmeyer SJ: Bcl-2-deficient mice demonstrate fulminant lymphoid apoptosis, polycystic kidneys, and hypopigmented hair. Cell. 1993, 75: 229-240. 10.1016/0092-8674(93)80065-M.PubMedView Article
- Fliegauf M, Benzing T, Omran H: When cilia go bad: cilia defects and ciliopathies. Nat Rev Mol Cell Biol. 2007, 8: 880-893. 10.1038/nrm2278.PubMedView Article
- Olbrich H, Fliegauf M, Hoefele J, Kispert A, Otto E, Volz A, Wolf MT, Sasmaz G, Trauer U, Reinhardt R, Sudbrak R, Antiqnac C, Gretz N, Walz G, Schermer B, Benzing T, Hildebrandt F, Omran H: Mutations in a novel gene, NPHP3, cause adolescent nephronophthisis, tapeto-retinal degeneration and hepatic fibrosis. Nat Genet. 2003, 34: 455-459. 10.1038/ng1216.PubMedView Article
- Omori S, Hida M, Fujita H, Takahashi H, Tanimura S, Kohno M, Awazu M: Extracellular signal-regulated kinase inhibition slows disease progression in mice with polycystic kidney disease. J Am Soc Nephrol. 2006, 17: 1604-1614. 10.1681/ASN.2004090800.PubMedView Article
- Tao Y, Kim J, Schrier RW, Edelstein CL: Rapamycin markedly slows disease progression in a rat model of polycystic kidney disease. J Am Soc Nephrol. 2005, 16: 46-51. 10.1681/ASN.2004080660.PubMedView Article
- Kaspareit-Rittinghausen J, Deerberg F, Rapp KG, Wcislo A: A new rat model for polycystic kidney disease of humans. Transplant Proc. 1990, 22: 2582-2583.PubMed
- Gretz N, Kranzlin B, Pey R, Schieren G, Bach J, Obermuller N, Ceccherini I, Kloting I, Rohmeiss P, Bachmann S, Hafner M: Rat models of autosomal dominant polycystic kidney disease. Nephrol Dial Transplant. 1996, 11 (Suppl 6): 46-51.PubMedView Article
- Brown JH, Bihoreau MT, Hoffmann S, Kranzlin B, Tychinskaya I, Obermuller N, Podlich D, Boehn SN, Kaisaki PJ, Megel N, Danoy P, Copley RR, Broxholme J, Witzgall R, Lathrop M, Gretz N, Gaugier D: Missense mutation in sterile alpha motif of novel protein SamCystin is associated with polycystic kidney disease in (cy/+) rat. J Am Soc Nephrol. 2005, 16: 3517-3526. 10.1681/ASN.2005060601.PubMedView Article
- Bihoreau MT, Megel N, Brown JH, Kranzlin B, Crombez L, Tychinskaya Y, Broxholme J, Kratz S, Bergmann V, Hoffmann S, Gauguier D, Gretz N: Characterization of a major modifier locus for polycystic kidney disease (Modpkdr1) in the Han:SPRD(cy/+) rat in a region conserved with a mouse modifier locus for Alport syndrome. Hum Mol Genet. 2002, 11: 2165-2173. 10.1093/hmg/11.18.2165.PubMedView Article
- Gretz N, Ceccherini I, Kranzlin B, Kloting I, Devoto M, Rohmeiss P, Hocher B, Waldherr R, Romeo G: Gender-dependent disease severity in autosomal polycystic kidney disease of rats. Kidney Int. 1995, 48: 496-500. 10.1038/ki.1995.319.PubMedView Article
- Castoldi M, Schmidt S, Benes V, Noerholm M, Kulozik AE, Hentze MW, Muckenthaler MU: A sensitive array for microRNA expression profiling (miChip) based on locked nucleic acids (LNA). RNA. 2006, 12: 913-920. 10.1261/rna.2332406.PubMedView Article
- Dai M, Wang P, Boyd AD, Kostov G, Athey B, Jones EG, Bunney WE, Myers RM, Speed TP, Akil H, Watson SJ, Meng F: Evolving gene/transcript definitions significantly alter the interpretation of GeneChip data. Nucleic Acids Res. 2005, 33: e175-10.1093/nar/gni179.PubMedView Article
- Shahi P, Loukianiouk S, Bohne-Lang A, Kenzelmann M, Kuffer S, Maertens S, Eils R, Grone HJ, Gretz N, Brors B: Argonaute – a database for gene regulation by mammalian microRNAs. Nucleic Acids Res. 2006, 34: D115-D118. 10.1093/nar/gkj093.PubMedView Article
- Lewis BP, Shih IH, Jones-Rhoades MW, Bartel DP, Burge CB: Prediction of mammalian microRNA targets. Cell. 2003, 115: 787-798. 10.1016/S0092-8674(03)01018-3.PubMedView Article
- John B, Enright AJ, Aravin A, Tuschl T, Sander C, Marks DS: Human MicroRNA targets. PLoS Biol. 2004, 2: e363-10.1371/journal.pbio.0020363.PubMedView Article
- Aoki KF, Kanehisa M: Using the KEGG database resource. Curr Protoc Bioinformatics. 2005, Chapter 1: Unit
- Harris MA, Clark J, Ireland A, Lomax J, Ashburner M, Foulger R, Eilbeck K, Lewis S, Marshall B, Mungall C, Richter J, Rubin GM, Blake JA, Bult C, Dolan M, Drabkin H, Eppig JT, Hill DP, Ni L, Ringwald M, Balakrishnan R, Cherry JM, Christie KR, Costanzo MC, Dwight SS, Engel S, Fisk DG, Hirschman JE, Hong EL, Nash RS, Sethuraman A, Theesfeld CL, Botstein D, Dolinski K, Feierbach B, Berardini T, Mundodi S, Rhee SY, Apweiler R, Barrell D, Camon E, Dimmer E, Lee V, Chisholm R, Gaudet P, Kibbe W, Kishore R, Schwarz EM, Sternberg P, Gwinn M, Hannick L, Wortman J, Berriman M, Wood V, de la Cruz N, Tonellato P, Jaiswal P, Seigfried T, White R: The Gene Ontology (GO) database and informatics resource. Nucleic Acids Res. 2004, 32: D258-D261. 10.1093/nar/gkh066.PubMedView Article
- Subramanian A, Kuehn H, Gould J, Tamayo P, Mesirov JP: GSEA-P: a desktop application for Gene Set Enrichment Analysis. Bioinformatics. 2007, 23: 3251-3253. 10.1093/bioinformatics/btm369.PubMedView Article
- Manoli T, Gretz N, Grone HJ, Kenzelmann M, Eils R, Brors B: Group testing for pathway analysis improves comparability of different microarray datasets. Bioinformatics. 2006, 22: 2500-2506. 10.1093/bioinformatics/btl424.PubMedView Article
- Bottoni A, Zatelli MC, Ferracin M, Tagliati F, Piccin D, Vignali C, Calin GA, Negrini M, Croce CM, Degli Uberti EC: Identification of differentially expressed microRNAs by microarray: a possible role for microRNA genes in pituitary adenomas. J Cell Physiol. 2007, 210: 370-377. 10.1002/jcp.20832.PubMedView Article
- Peltier HJ, Latham GJ: Normalization of microRNA expression levels in quantitative RT-PCR assays: identification of suitable reference RNA targets in normal and cancerous human solid tissues. RNA. 2008, 14: 844-852. 10.1261/rna.939908.PubMedView Article
- Berezikov E, Guryev V, van de BJ, Wienholds E, Plasterk RH, Cuppen E: Phylogenetic shadowing and computational identification of human microRNA genes. Cell. 2005, 120: 21-24. 10.1016/j.cell.2004.12.031.PubMedView Article
- Murcia NS, Sweeney WE, Avner ED: New insights into the molecular pathophysiology of polycystic kidney disease. Kidney Int. 1999, 55: 1187-1197. 10.1046/j.1523-1755.1999.00370.x.PubMedView Article
- Grantham JJ: The etiology, pathogenesis, and treatment of autosomal dominant polycystic kidney disease: recent advances. Am J Kidney Dis. 1996, 28: 788-803. 10.1016/S0272-6386(96)90378-9.PubMedView Article
- Lanoix J, D'Agati V, Szabolcs M, Trudel M: Dysregulation of cellular proliferation and apoptosis mediates human autosomal dominant polycystic kidney disease (ADPKD). Oncogene. 1996, 13: 1153-1160.PubMed
- Sorenson CM, Sheibani N: Sustained activation of MAPK/ERKs signaling pathway in cystic kidneys from bcl-2 -/- mice. Am J Physiol Renal Physiol. 2002, 283: F1085-F1090.PubMedView Article
- Benzing T, Simons M, Walz G: Wnt signaling in polycystic kidney disease. J Am Soc Nephrol. 2007, 18: 1389-1398. 10.1681/ASN.2006121355.PubMedView Article
- Shillingford JM, Murcia NS, Larson CH, Low SH, Hedgepeth R, Brown N, Flask CA, Novick AC, Goldfarb DA, Kramer-Zucker A, Walz G, Piontek KB, Germino GG, Weimbs T: The mTOR pathway is regulated by polycystin-1, and its inhibition reverses renal cystogenesis in polycystic kidney disease. Proc Natl Acad Sci USA. 2006, 103: 5466-5471. 10.1073/pnas.0509694103.PubMedView Article
- Jackson RJ, Standart N: How do microRNAs regulate gene expression?. Sci STKE. 2007, 2007: re1-10.1126/stke.3672007re1.PubMedView Article
- O'Connell RM, Taganov KD, Boldin MP, Cheng G, Baltimore D: MicroRNA-155 is induced during the macrophage inflammatory response. Proc Natl Acad Sci USA. 2007, 104: 1604-1609. 10.1073/pnas.0610731104.PubMedView Article
- Vinther J, Hedegaard MM, Gardner PP, Andersen JS, Arctander P: Identification of miRNA targets with stable isotope labeling by amino acids in cell culture. Nucleic Acids Res. 2006, 34: e107-10.1093/nar/gkl590.PubMedView Article
- Tian Z, Greene AS, Pietrusz JL, Matus IR, Liang M: MicroRNA-target pairs in the rat kidney identified by microRNA microarray, proteomic, and bioinformatic analysis. Genome Res. 2008, 18: 404-411. 10.1101/gr.6587008.PubMedView Article
- Selbach M, Schwanhausser B, Thierfelder N, Fang Z, Khanin R, Rajewsky N: Widespread changes in protein synthesis induced by microRNAs. Nature. 2008, 455: 58-63. 10.1038/nature07228.PubMedView Article
- Baek D, Villen J, Shin C, Camargo FD, Gygi SP, Bartel DP: The impact of microRNAs on protein output. Nature. 2008, 455: 64-71. 10.1038/nature07242.PubMedView Article
- Lantinga-van LI, Leonhard WN, Dauwerse H, Baelde HJ, van Oost BA, Breuning MH, Peters DJ: Common regulatory elements in the polycystic kidney disease 1 and 2 promoter regions. Eur J Hum Genet. 2005, 13: 649-659. 10.1038/sj.ejhg.5201392.View Article
- Lu W, Peissel B, Babakhanlou H, Pavlova A, Geng L, Fan X, Larson C, Brent G, Zhou J: Perinatal lethality with kidney and pancreas defects in mice with a targetted Pkd1 mutation. Nat Genet. 1997, 17: 179-181. 10.1038/ng1097-179.PubMedView Article
- Perera RJ, Ray A: MicroRNAs in the search for understanding human diseases. BioDrugs. 2007, 21: 97-104. 10.2165/00063030-200721020-00004.PubMedView Article
- Zhou B, Wang S, Mayr C, Bartel DP, Lodish HF: miR-150, a microRNA expressed in mature B and T cells, blocks early B cell development when expressed prematurely. Proc Natl Acad Sci USA. 2007, 104: 7080-7085. 10.1073/pnas.0702409104.PubMedView Article
- Bagnyukova TV, Pogribny IP, Chekhun VF: MicroRNAs in normal and cancer cells: a new class of gene expression regulators. Exp Oncol. 2006, 28: 263-269.PubMed
- Kumar MS, Lu J, Mercer KL, Golub TR, Jacks T: Impaired microRNA processing enhances cellular transformation and tumorigenesis. Nat Genet. 2007, 39: 673-677. 10.1038/ng2003.PubMedView Article
- Marshall CJ: 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.PubMedView Article
- Sheibani N, Sorenson CM, Frazier WA: Differential modulation of cadherin-mediated cell-cell adhesion by platelet endothelial cell adhesion molecule-1 isoforms through activation of extracellular regulated kinases. Mol Biol Cell. 2000, 11: 2793-2802.PubMedView Article
- Chen S, Jim B, Ziyadeh FN: Diabetic nephropathy and transforming growth factor-beta: transforming our view of glomerulosclerosis and fibrosis build-up. Semin Nephrol. 2003, 23: 532-543. 10.1053/S0270-9295(03)00132-3.PubMedView Article
- Obermuller N, Morente N, Kranzlin B, Gretz N, Witzgall R: A possible role for metalloproteinases in renal cyst development. Am J Physiol Renal Physiol. 2001, 280: F540-F550.PubMed
- Ziyadeh FN, Sharma K: Role of transforming growth factor-beta in diabetic glomerulosclerosis and renal hypertrophy. Kidney Int Suppl. 1995, 51: S34-S36.PubMed
- Kato M, Zhang J, Wang M, Lanting L, Yuan H, Rossi JJ, 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: 3432-3437. 10.1073/pnas.0611192104.PubMedView Article
- Calvet JP: Molecular genetics of polycystic kidney disease. J Nephrol. 1998, 11: 24-34.PubMed
- Neilson JR, Zheng GX, Burge CB, Sharp PA: Dynamic regulation of miRNA expression in ordered stages of cellular development. Genes Dev. 2007, 21: 578-589. 10.1101/gad.1522907.PubMedView Article
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.