BioMed Research International

BioMed Research International / 2017 / Article

Research Article | Open Access

Volume 2017 |Article ID 8306736 | 9 pages | https://doi.org/10.1155/2017/8306736

Integrative Analysis of miRNA and mRNA Expression Profiles in Calcium Oxalate Nephrolithiasis Rat Model

Academic Editor: Sergio Claudio Saccà
Received01 Sep 2017
Accepted06 Nov 2017
Published17 Dec 2017

Abstract

The microRNA (miRNA) expression profiles and their biological functions in calcium oxalate nephrolithiasis remain unclear. In this study, we investigate the miRNA and mRNA expression profiles of kidney tissues in calcium oxalate stone rats. 16 Sprague Dawley rats were divided into control group and stone-forming group. 24-hour urine samples and kidney tissues were collected for biochemical and histological determination after 4 weeks. MiRNA and mRNA microarray were applied to evaluate the miRNA and mRNA expression profiles. To validate the microarray results, the quantitative real-time PCR (qRT-PCR) was performed. A total of 38 miRNAs and 2728 mRNAs were significantly and differentially expressed in kidney tissues of stone-forming group versus control group. Gene Ontology (GO) analysis revealed that most of the target genes were enriched in terms of oxidation reduction, ion transport, inflammatory response, and response to wounding. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis of these targets highlights their critical role in cytokine-cytokine receptor interaction, gap junction, and chemokine signaling pathway. Furthermore, the reliability of the microarray-based results was confirmed by using qRT-PCR determination. The miRNA and mRNA expressions in calcium oxalate stone rat kidneys might provide a basis for further research on urolithiasis mechanism.

1. Introduction

Urolithiasis is a worldwide disease with high morbidity and economic costs. The kidney stone prevalence in US showed a linear increase over the past four decades, that is, from 3.8% during 1976–1980 to 8.8% during 2007–2010 [1, 2]. It has been reported that the prevalence in China was 6.50% [3]. More than half of the stone patients experience a recurrence within 5–10 years of a first episode [4]. Calcium oxalate, the primary constituent of urolithiasis, accounts for approximately 80% of the total urinary tract stone [5], whose pathogenesis remains a mystery though existing theories showed that it was a multifactorial disorder including genetic, metabolic, and dietary factors [6].

MicroRNAs (miRNAs), a class of endogenous noncoding small RNA molecules of 21–23 nucleotides, play a crucial regulatory role in degrading target mRNAs or repressing their translation by binding to complementary sites on target mRNAs in the 3′UTR [7]. In eukaryotic cells, miRNAs can regulate 31% of the target genes and participate in multiple processes of biological development and diseases [8]. Recent studies implicated that up- or downregulation of certain miRNAs might play a major role in the cause of hypercalciuria, which is an important risk factor for calcium urolithiasis [9]. However, there is a paucity of research about miRNA expression patterns and their biological functions in urolithiasis, which seems important to further understand the relationship between genes and pathophysiological process of urolithiasis.

In this study, we applied miRNA, mRNA microarray technique, and bioinformatics databases to integrally analyze the miRNA and mRNA expression profiles of kidney tissues in calcium oxalate stone rats, aiming at exploring more specific and effective molecular targets for urolithiasis.

2. Materials and Methods

2.1. Animal Model of Urolithiasis

Hyperoxaluria is one of the most common risk factors of calcium oxalate nephrolithiasis, which can be induced by sustained administration of oxalate or its precursors such as ethylene glycol (EG). During administration of EG with ammonium chloride (AC), the urinary excretion of oxalate increased rapidly and calcium oxalate supersaturation increased accordingly. Then crystalluria was produced and was followed by calcium oxalate nephrolithiasis [10]. In the present study, a total of 16 male Sprague Dawley rats ( g) were utilized and purchased from Guangdong Medical Laboratory Animal Center, China. They were kept in standard laboratory conditions (under natural light and dark cycles) and received food and water ad libitum before the experiment. 16 rats were randomly divided into two groups of 8 rats each. In the control group, rats were fed with free access to food and normal drinking water, while, in the stone-forming group, rats were fed with free access to food and drinking water containing 1% EG for 4 weeks. Besides, 1% AC was given by gavage for 5 days at the early stage of experiment. The experimental protocol was approved by the Ethics Review Board of Guangzhou Medical University, Guangzhou, China.

2.2. Collection of Urine and Kidney Specimen

The duration of the experiment was 4 weeks. 24-hour urine samples were collected on the day before the end of the experiment. At the end of the experimental period, all the rats were sacrificed by overdose of pentobarbital sodium (50 mg/kg, i.p.) and then both sides of the kidneys were removed after sufficient transcardial perfusion. One kidney was stored in liquid nitrogen for molecular biological studies and the other was fixed in 4% paraformaldehyde for histological studies.

2.3. Biochemical Determinations of Urine and Histological Studies

The urine oxalate and citrate were measured using an ion chromatograph (883 Basic IC plus, Metrohm AG, Switzerland) while the urine calcium was determined by the automatic biochemical analyzer. Rat kidney tissues were embedded in paraffin. Cross-sections of 5 μm were prepared and stained with hematoxylin/eosin (HE). Then the tissues and crystals were examined using polarizing microscope.

2.4. RNA Extraction

Total RNA was extracted from kidney cortex tissue adjacent to stones using Trizol Reagent (Ambion, USA) according to the manufacturer’s protocol. The purity and concentration were assessed by NanoPhotometer (Implen, Germany).

2.5. Microarray Analysis

Three stone-forming rats and three control rats were randomly selected for microarray analysis. The Agilent Rat miRNA Microarray Release 19.0 and Agilent Rat Microarray  K were used to compare miRNA and mRNA expression profiles, respectively, in both control and stone-forming rats. Microarray analyses were performed by CapitalBio Corporation (Beijing, China) and Biotechnology Corporation (Shanghai, China), respectively.

2.6. Real-Time RT-PCR

For miRNA analysis, total RNA was used for cDNA synthesis and quantitative detection according to the protocol of the All-in-One™ miRNA qRT-PCR Detection Kit (GeneCopoeia, Guangzhou, China). For mRNA analysis, total RNA was reverse-transcribed to cDNA using Takara reverse transcription kit (Takara, Dalian, China). PCR was conducted according to standard protocol of SYBR Premix Ex Taq (Takara, Dalian, China). The primers were chemically synthesized by Shanghai Generay Biotech Co., Ltd. U6 and GAPDH were used as the internal controls of miRNA and mRNA, respectively. All reactions were performed in triplicate.

2.7. Target Prediction and Function Analysis

Five bioinformatics databases, including TargetScan, miRanda, miRDB, miRWalk, and PITA, were utilized to predict the potential target genes. Genes present in three of the five databases were considered as potential target mRNAs. Then these genes were compared with those from mRNA microarray, and the overlaps of them were determined. As miRNA exerts its biological functions through the suppression of target genes, that is, miRNA negatively modulates gene expression, we intersected the predicted target mRNAs of upregulated miRNAs with the downregulated mRNAs from microarray, and vice versa. In order to determine the biological functions and pathways that were potentially functional in the formation of kidney stone, we performed Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) annotations of the predicted target mRNAs by DAVID Bioinformatics Resources.

2.8. Statistical Analysis

Continuous variables were expressed as means ± standard deviation; then Student’s -test was applied for statistical analysis by using SPSS 13.0. The raw data of microarray were analyzed by using GeneSpring and the RVM (random variance model) -test was performed to filter the differentially expressed miRNAs/mRNAs between the two groups [11]. A value of <0.05 was considered to be statistically significant.

3. Results

3.1. Urine Biochemical Analysis and Histological Observation

Comparisons of biochemical analysis of urine between two groups of rats are illustrated in Figure 1. After a four-week ethylene glycol administration, the average urine oxalate excretion in stone-forming rats was significantly enhanced compared to the control rats ( versus  mmol/24 h), whereas urinary calcium was significantly reduced ( versus  mmol/24 h). Urine citrate and volume, however, were not significantly altered in the two groups.

The histological changes of kidney sections in the two groups of rats are shown in Figure 2. In stone-forming group, the crystal deposits were displayed in the renal tubular lumen of cortex and medulla region and were irregular in shape and size. Moreover, most of the renal tubules were distended. By contrast, the kidney of control group showed normal morphology and no crystals were observed in them at all.

3.2. Alteration of miRNA and mRNA Expression Profiles

The miRNA expression profiles were significantly different in rat kidney tissues between the two groups. Of all the miRNAs measured, 38 mature miRNAs were identified. In the rats of stone-forming group, the deferential expressions revealed 2-fold changes compared to that in the rats of control group (), including 18 upregulated miRNAs and 20 downregulated ones. Among these differentially expressed miRNAs, rno-miR-130b-3p, rno-miR-132-3p, rno-miR-181a-1-3p, rno-miR-222-3p, rno-miR-351-5p, and rno-miR-21-3p had the largest positive fold changes, while rno-miR-335, rno-miR-192-3p, rno-miR-194-5p, rno-miR-192-5p, rno-miR-499-5p, and rno-miR-210-3p had the largest negative fold changes (Table 1).


Upregulated miRNAsFold change valueDownregulated miRNAsFold change value

rno-miR-130b-3p153.930.0033rno-miR-33592.070.038
rno-miR-132-3p135.680.016rno-miR-192-3p56.460.036
rno-miR-181a-1-3p88.120.0033rno-miR-194-5p3.750.011
rno-miR-222-3p44.970.012rno-miR-192-5p3.680.012
rno-miR-351-5p33.890.0033rno-miR-499-5p2.980.023
rno-miR-21-3p28.020.043rno-miR-210-3p2.890.041
rno-miR-146b-5p10.730.0068rno-miR-200b-5p2.880.030
rno-miR-503-5p4.710.0011rno-miR-3472.880.048
rno-miR-21-5p4.680.011rno-miR-29c-3p2.640.00048
rno-miR-34b-5p4.030.016rno-miR-29c-5p2.560.0032
rno-miR-542-5p3.620.039rno-miR-185-5p2.390.013
rno-miR-223-3p3.210.012rno-miR-218a-5p2.380.012
rno-miR-18a-5p2.980.038rno-miR-378a-5p2.340.0033
rno-miR-214-3p2.830.024rno-miR-378b2.270.011
rno-miR-542-3p2.650.015rno-miR-30e-3p2.260.0033
rno-miR-199a-5p2.580.016rno-miR-30c-5p2.240.0042
rno-miR-322-5p2.420.014rno-miR-30e-5p2.230.014
rno-miR-142-3p2.290.011rno-miR-203a-3p2.170.036
rno-miR-200a-3p2.100.0033
rno-miR-4292.010.0042

In order to investigate potential targets of the altered miRNAs, we also determined the mRNA expression profiles of the same kidney tissues by using Agilent microarray. The results showed that, among the 2728 differentially expressed mRNAs, 1350 mRNAs were significantly increased whereas 1378 mRNAs were significantly decreased in stone-forming group. The top 10 upregulated genes were Il19, Reg3b, Havcr1, C4bpa, Mmp7, Il24, Fgg, Fgb, Gpnmb, and Cxcl2. Among the downregulated genes, the top 10 were Cyp2c11, LOC102550988, RGD1306462, Serpinb12, LOC361914, Cyp1a1, Cacng5, Iqch, Pkd2l2, and LOC691551 (Table 2).


Upregulated mRNAsFold change valueDownregulated mRNAsFold change value

Il191212.010.00028Cyp2c11307.340.0018
Reg3b739.23LOC102550988118.350.0065
Havcr1404.69RGD130646269.720.0083
C4bpa402.770.0065Serpinb1267.620.0063
Mmp7354.620.0047LOC36191454.30.0065
Il24242.970.0036Cyp1a145.190.00095
Fgg187.730.00050Cacng544.060.00011
Fgb179.11Iqch37.510.00064
Gpnmb139.070.00035Pkd2l229.500.0063
Cxcl2105.480.0067LOC69155129.270.0074

3.3. miRNA Target Gene Prediction

As miRNAs negatively regulate mRNAs expression via degradation or translation inhibition, we intersect the predicted target mRNAs of upregulated miRNAs with the downregulated mRNAs from microarray, and vice versa. Among the 1378 significantly decreased mRNAs, 141 genes were the predicted targets of the 18 upregulated miRNAs. Similarly, 187 in 1350 overexpressed mRNAs and 20 downregulated miRNAs were predicted as targets. According to the above method, a total of 797 significant miRNA-mRNA pairs were predicted, consisting of the 38 differentially expressed miRNAs and 328 mRNAs. The predicted targets of these miRNAs were then used for downstream GO analysis and KEGG pathway analysis.

3.4. GO and KEGG Pathway Analysis

Functional enrichment analysis of the target genes of the differentially expressed miRNAs was performed by means of Gene Ontology system, and a few representative biological processes were revealed in Table 3. According to the GO analysis, there were 362 GO terms significantly enriched (). The most highly enriched biological function of the upregulated miRNAs included oxidation reduction (GO:0055114), response to hormone stimulus (GO:0009725), and ion transport (GO:0006811), while response to wounding (GO:0009611), inflammatory response (GO:0006954), and defense response (GO:0006952) were associated with the downregulated miRNAs.


GO IDGO termTarget genes value

Upregulated miRNAs
GO:0055114Oxidation reductionME1, BCKDHA, LDHC, SUOX, GPD1, ADHFE1, CYP1A1, HSD3B6, GLUD1, CYCS, HGD, GCLM, HIBADH, ACOX3, FDFT1, IVD, DHCR7, HSD11B1, AKR7A3, DMGDH, MECR, ALDH9A1, NQO2
GO:0009725Response to hormone stimulusME1, BCKDHA, AR, CRYAB, STAT5A, GGH, AQP1, PPARGC1B, GHRHR, PRSS8, EIF4EBP2, WFDC1, GHRL, TGFBR3, GHR, GNG7
GO:0006811Ion transportSLC36A1, SLC38A3, CASR, FXYD4, AQP1, FXYD6, BSND, SLC23A1, ATP6V0E2, P2RY4, ATP6V1E1, SLC39A8, SCN4B, SLC13A2, PLLP, SLC5A6, SLC4A1, ATP6V0D10.00018

Downregulated miRNAs
GO:0009611Response to woundingCXCL1, MASP1, SCN3A, TACR1, CCR1, CXCL2, GJA1, TLR4, CCL7, TGFB2, TIMP1, VCAM1, PCSK1, CASP3, FGG, CD44, SERPINE1, IL1B, REG3G, KLF6, TNFSF4, LYN, MAP1B, IL1RN, C4BPA, IL24, P2RY12, FCGR2B, F3, CCR2, NPPB, IGFBP1, PTAFR, IGFBP4
GO:0006954Inflammatory responseCXCL1, TNFSF4, LYN, MASP1, TACR1, CCR1, IL1RN, CXCL2, TLR4, C4BPA, CCL7, VCAM1, FGG, CD44, FCGR2B, F3, CCR2, IL1B, NPPB, REG3G, PTAFR, IGFBP4
GO:0006952Defense responseCXCL1, MASP1, FGR, GRIK2, TACR1, CCR1, CXCL2, TLR4, CCL7, VCAM1, FGG, CD44, IL1B, REG3G, PTPRC, TNFSF4, LYN, IL1RN, C4BPA, FCGR2B, F3, CCR2, NPPB, PTAFR, IGFBP4

Potential functional pathways for the targets of differentially expressed miRNAs were then predicted by KEGG pathway system. Based on the KEGG pathway analysis, 19 pathways were predicted to be significantly related to stone formation (). Cytokine-cytokine receptor interaction, gap junction, and chemokine signaling pathway were the main pathways involved in these predicted pathways (Table 4).


KEGG pathway termTarget genes value

rno00280: valine, leucine, isoleucine degradationBCKDHA, MCCC2, IVD, ECHS1, ACAT1, HIBADH, ALDH9A10.0018
rno04060: cytokine-cytokine receptor interactionCCL22, CNTF, TNFSF4, IL9R, TNFRSF12A, CLCF1, CCR1, CXCL2, CCR2, IL21R, IL1B, FAS, GHR0.010
rno04115: p53 signaling pathwayCASP3, CDKN1A, RRM2, CASP8, CYCS, SERPINE1, FAS0.011
rno04540: gap junctionTUBB2B, TUBB5, GJA1, GUCY1A3, GUCY1B3, HTR2B, TUBB30.027
rno00380: tryptophan metabolismCCBL1, CYP1A1, ECHS1, ACAT1, ALDH9A10.033
rno04080: neuroactive ligand-receptor interactionADORA2B, GRIK2, TACR1, GRIA3, GHRHR, P2RY4, SSTR1, CNR1, P2RY2, CCR2, HTR2B, PTAFR, ADRA1D, GHR0.035

3.5. Validation of the Microarray Expression Data by Real-Time PCR

To validate the reliability of our microarray-based results, qRT-PCR was performed in 5 pairs of matched stone-forming/control rat kidneys. From the differentially expressed miRNAs and mRNAs observed by microarray, 15 miRNAs and 10 mRNAs were selected. Relative expression levels of the selected miRNAs and mRNAs were depicted in Figure 3. Consistent with the microarray data, real-time PCR confirmed that, compared with controls, rno-miR-132-3p, rno-miR-181a-1-3p, rno-miR-222-3p, and rno-miR-351-5p were significantly increased, while rno-miR-192-3p, rno-miR-194-5p, rno-miR-29c-3p, rno-miR-185-5p, and rno-miR-30c-5p were significantly decreased in stone-forming rat kidneys. Besides, the expressions of Il19, Reg3b, Havcr1, C4bpa, and Mmp7 were significantly increased, and those of Serpinb12, Cacng5, Pkd2l2, and Casr were significantly decreased in stone-forming group, which is in line with the microarray data, too. Although the magnitude of changes differed between the two methods, these results demonstrated a high consistency between the microarray and RT-qPCR based results except for rno-miR-130b-3p and Cyp2c11.

4. Discussion

In the present study, we found 38 miRNAs and 2728 mRNAs differentially expressed in kidney tissues between calcium oxalate stone-forming rats and the controls. 328 potential target genes for the 38 miRNAs were predicted by intersecting five bioinformatics databases and differentially expressed mRNAs, resulting in 797 miRNA-mRNA pairs. GO term and KEGG pathway analysis further revealed that the miRNA-mRNA contributed a lot to the calcium oxalate stone formation. In addition, the expression of several miRNAs and mRNAs was confirmed via qRT-PCR and the outcomes were highly consistent with the microarray data, indicating that our microarray-based determinations were reliable.

As posttranscriptional gene regulators, miRNAs are involved in diverse pathophysiological processes of various diseases [8]. There are growing evidences that miRNAs play a vital role in the formation of kidney stone. Hu et al. [12] reported that the miR-155 level in serum and urine was significantly increased in patients with nephrolithiasis and proposed that miR-155 might influence the nephrolithiasis pathophysiological process by regulating inflammatory cytokines expression. Though there was no difference of miR-155 in the expression between the two groups of our study, we found that the targets of miRNAs were enriched in the biological function of inflammatory response. Another study also reported that miRNAs were significantly differently expressed in proximal renal tubular cells due to exposure to calcium oxalate monohydrate crystal and some of them might be involved in the insulin signaling pathway and type 2 diabetes mellitus [13].

Recently, Liu and his colleagues [14] established a hyperoxaluric rat model induced by feeding of ethylene glycol and ammonium chloride and demonstrated that the target genes and signaling pathways were involved in the pathogenesis of hyperoxaluria. They showed that 28 miRNAs were differentially expressed between the experimental and control groups. In addition, the GO and KEGG analyses revealed that these miRNAs were potentially associated with the insulin resistance and phosphatidylinositol-bisphosphonate 3-kinase/AKT serine threonine kinase signaling pathways. Among the miRNAs differentially expressed in Liu’s study, rno-miR-132-3p, rno-miR-146b-5p, rno-miR-223-3p, rno-miR-21-5p, and rno-miR-214-3p were upregulated, which was just the same as our findings. However, their rat models were created by free access to water containing 1% EG and 0.5% AC for two weeks. Since kidney stone is a chronic disease, 2 weeks of modeling time might be not conformable with regard to the natural course of stone formation. In addition, AC was reported to have nephrotoxicity, which might cause kidney injury and in consequence affect the outcomes of the study when used for a longer time. Therefore, our rat models were induced by drinking water with 1% EG for 4 weeks, and 1% AC was given as gavage during the early 5 days. Theoretically, this model would be more convincing and the obtained data more reliable.

Lu et al. [9] reported that there were 19 miRNAs significantly expressed in genetic hypercalciuric stone-forming rat kidneys compared to normal kidneys, including 10 upregulated miRNAs and 9 downregulated ones. Among them, rno-miR-674-5p, rno-miR-672-5p, rno-miR-138-1-3p, and rno-miR-21-3p were significantly expressed and might play important roles in the regulatory network. KEGG pathway analysis showed that NF-κB signaling pathway was significantly enriched in target genes of differentially expressed miRNAs, suggesting that it might be associated with hypercalciuria urolithiasis. Out of the 19 reported miRNAs that were dysregulated, miR-21-3p was the only one miRNA identified as the upregulated gene in the above-mentioned two studies.

With respect to oncogenic miRNA, miR-21 was traditionally regarded as highly associated with variety of tumors, such as breast cancer, ovarian cancer, or urothelial carcinoma [1517]. Recent study suggested that miR-21 might participate in the development of fibrosis via promoting the proliferation of interstitial fibroblasts and increasing the abnormal deposition of the extracellular matrix [18]. Other investigators showed that miR-21 could inhibit CD40 expression through SIRT1-NF-κB signaling, which might negatively regulate inflammatory processes in renal tubule epithelial cell [1921].

Moreover, miR-21 was regarded as modulatory reaction of oxidative stress since the SOD expression was decreased and ROS formation increased [22]. It prompts us to believe that, no matter whether mediated by hypercalciuria or hyperoxaluria, miR-21 was the common miRNA that was overexpressed in the process of stone formation, which might subsequently cause oxidative stress or stimulate renal tubular epithelial cells apoptosis and tubulointerstitial fibroblast proliferation, thus resulting in the fibrotic dysfunction of kidney by certain downstream signaling pathways such as PTEN/Akt pathway [23]. Therefore, miR-21 might be crucial in the processes of stone formation, and the mechanism needs to be further investigated.

Based on the GO analysis, the most highly enriched biological functions of dysregulated miRNAs targets are oxidation reduction, ion transport, response to wounding and inflammation, and so on. Similar study showed that oxidative stress was intimately involved in signaling molecule changes, renal tubular epithelial cell injury, and inflammation during the renal stone formation [6, 24, 25]. The oxidants such as ROS could severely affect the cells structure and function in consequence of reacting with all the cell basic constituents [26]. Both in vivo and in vitro animal experiments demonstrated that ROS could generate the formation of kidney stone, and treatment with antioxidants significantly reduces the crystal deposition [2729]. Obviously, the crystal-associated pathological changes in kidney resulted from the damage effects of oxidants and from the oxidants-mediated changes in gene expression and signal transduction.

Inflammation is a complex biological response to various irritants and inflammatory cells. The inflammatory cells, such as monocyte and macrophages, are usually found to migrate to the sites of crystal deposition, which is probably due to the effects of inflammatory cytokines, chemokines, and adhesion molecules released by renal tubular epithelial cells [21, 30]. Normally, once crystals have been formed, they would be soon cleared away by the protective response, that is, producing macromolecules and phagocytose to eliminate them, rendering them harmless [31]. In agreement, we showed that the expression of rno-miR-192, rno-miR-194, and rno-miR-499 was higher in the stone-forming group while their potential target gene chemokine receptor 2 (CCR2) was lower. Another study demonstrated that CCR2, the receptor of chemoattractant protein-1 (MCP-1), could stimulate macrophages infiltration, tubulointerstitial injury, and fibrosis, which was closely associated with stone formation [32]. Therefore, we proposed that, in the process of kidney stone formation, the overexpression of CCR2 mediated by relevant miRNAs, such as rno-miR-192, rno-miR-194, and rno-miR-499, would induce the inflammation and damage to the renal tubular epithelial cells and promote nephrolithiasis. Ion transport is another major pathogenesis for the formation of nephrolithiasis. In this study, we showed that rno-miR-223-3p was upregulated whereas CaSR was downregulated in the kidney of stone-forming group, indicating that CaSR was a potential target gene of rno-miR-223-3p. Recent research showed that the expression of CLDN14 was increased in mice with high dietary calcium intake. It was assumed that this regulation might be via the activation of CaSR, blockage of the paracellular reabsorption of Ca2+, and facilitation of the renal Ca2+ losses [3335]. Therefore, the disturbance of this signaling pathway might contribute to calcium changes by renal calcium excretion, which is closely associated with stone formation.

In addition, as predicted by the bioinformatics databases, SLC4A1 was the potential target gene of rno-miR-34b, rno-miR-146b, rno-miR-214, rno-miR-223, and rno-miR-351. SLC4A1, which encoded anion exchanger-1 (AE1), was manifested as bicarbonate reabsorption and urinary acidification [36]. SLC4A1 was differentially downregulated in stone-forming group in our study. Previous reports suggested that the abnormal function of SLC4A1 was a risk factor for idiopathic calcium nephrolithiasis, and the mutations of it might cause distal renal tubular acidosis, which also promoted recurrent renal stone formation [37, 38]. Therefore, in the process of stone formation, the upregulation of these miRNAs might downregulate the expression of SLC4A1, which impacted the acid-base homeostasis and led to stone formation.

In previous studies, pathways such as p38 MAPK pathway and TGF pathway had been reported to be related to the process of stone formation [28, 39]. However, our results did not draw the same conclusion. The KEGG pathway analysis indicated that cytokine-cytokine receptor interaction and chemokine signaling pathway were significantly enriched in stone-forming rats. It seems that KEGG pathways are consistent with the GO terms because both of them were involved in the development of inflammation.

The findings of the present study may provide novel biomarkers for diagnosis, prognosis, and clinical trials of calcium oxalate nephrolithiasis. However, as the urine calcium excretion was lower in the stone-forming group, we are not so sure whether the miRNAs and mRNAs observed in this animal stone model are consistent with the natural pathogenesis of human stone disease. Other limitations of this study are whether the altered mRNAs could be translated into the similar proteins and what the exact target gene and biological function of these dysregulated miRNAs are. All these issues need further study for clarification.

5. Conclusion

In summary, the identification of differentially expressed miRNAs and mRNAs in the kidneys between calcium oxalate stone-formation rats and normal ones was determined. Some pathway analysis in this study, such as cytokine-cytokine receptor interaction, gap junction, and chemokine signaling pathway, provides a basis and new direction for further research on urolithiasis mechanism.

Conflicts of Interest

All authors have no conflicts of interest.

Acknowledgments

This study was supported in part by research grants from Science and Technology Project of Guangzhou (201604020001 and 201607010359), General Project of Guangzhou Education Bureau (1201410963), National Natural Science Foundation of China (81570633 and 81400708), Educational Commission of Guangdong Province of China (2015KTSCX115), and Foundation of Distinguished Young Talents in Higher Education of Guangdong (2014KQNCX121).

References

  1. C. D. Scales Jr., A. C. Smith, J. M. Hanley, and C. S. Saigal, “Urologic diseases in America Project. Prevalence of kidney stones in the United States,” European Urology, vol. 62, no. 1, pp. 160–165, 2012. View at: Publisher Site | Google Scholar
  2. K. K. Stamatelou, M. E. Francis, C. A. Jones, L. M. Nyberg Jr., and G. C. Curhan, “Time trends in reported prevalence of kidney stones in the United States: 1976–1994,” Kidney International, vol. 63, no. 5, pp. 1817–1823, 2003. View at: Publisher Site | Google Scholar
  3. G. Zeng, Z. Mai, S. Xia et al., “A cross-sectional survey of urolithiasis prevalence in China,” Chinese Journal of Urology, vol. 36, no. 7, pp. 528–532, 2015. View at: Google Scholar
  4. J. Uribarri, M. S. Oh, and H. J. Carroll, “The first kidney stone,” Annals of Internal Medicine, vol. 111, no. 12, pp. 1006–1009, 1989. View at: Publisher Site | Google Scholar
  5. E. M. Worcester and F. L. Coe, “Clinical practice. Calcium kidney stones,” The New England Journal of Medicine, vol. 363, no. 10, pp. 954–963, 2010. View at: Publisher Site | Google Scholar
  6. O. W. Moe, “Kidney stones: pathophysiology and medical management,” The Lancet, vol. 367, no. 9507, pp. 333–344, 2006. View at: Publisher Site | Google Scholar
  7. V. Ambros, “microRNAs: tiny regulators with great potential,” Cell, vol. 107, no. 7, pp. 823–826, 2001. View at: Publisher Site | Google Scholar
  8. B. P. Lewis, I.-H. Shih, M. W. Jones-Rhoades, D. P. Bartel, and C. B. Burge, “Prediction of mammalian microRNA targets,” Cell, vol. 115, no. 7, pp. 787–798, 2003. View at: Publisher Site | Google Scholar
  9. Y. Lu, B. Qin, H. Hu et al., “Integrative microRNA-gene expression network analysis in genetic hypercalciuric stone-forming rat kidney,” PeerJ, vol. 2016, no. 3, article e1884, 2016. View at: Publisher Site | Google Scholar
  10. S. R. Khan, “Animal models of calcium oxalate kidney stone formation,” in Animal Models for the Study of Human Disease, pp. 483–498, Elsevier, 2013. View at: Google Scholar
  11. G. W. Wright and R. M. Simon, “A random variance model for detection of differential gene expression in small microarray experiments,” Bioinformatics, vol. 19, no. 18, pp. 2448–2455, 2003. View at: Publisher Site | Google Scholar
  12. Y. Y. Hu, W. D. Dong, Y. F. Xu et al., “Elevated levels of miR-155 in blood and urine from patients with nephrolithiasis,” BioMed Research International, vol. 2014, Article ID 295651, 7 pages, 2014. View at: Publisher Site | Google Scholar
  13. B. Wang, B. Wu, J. Liu et al., “Analysis of altered MicroRNA expression profiles in proximal renal tubular cells in response to calcium oxalate monohydrate crystal adhesion: Implications for kidney stone disease,” PLoS ONE, vol. 9, no. 7, Article ID e101306, 2014. View at: Publisher Site | Google Scholar
  14. Z. Liu, H. Jiang, J. Yang et al., “Analysis of altered microRNA expression profiles in the kidney tissues of ethylene glycol-induced hyperoxaluric rats,” Molecular Medicine Reports, vol. 14, no. 5, pp. 4650–4658, 2016. View at: Publisher Site | Google Scholar
  15. P. M. Báez-Vega, I. M. E. Vargas, F. Valiyeva et al., “Targeting miR-21-3p inhibits proliferation and invasion of ovarian cancer cells,” Oncotarget , vol. 7, no. 24, pp. 36321–36337, 2016. View at: Publisher Site | Google Scholar
  16. I. Casaburi, M. G. Cesario, A. Donà et al., “Androgens downregulate miR-21 expression in breast cancer cells underlining the protective role of androgen receptor,” Oncotarget , vol. 7, no. 11, pp. 12651–12661, 2016. View at: Publisher Site | Google Scholar
  17. R. Ohno, H. Uozaki, Y. Kikuchi et al., “Both cancerous miR-21 and stromal miR-21 in urothelial carcinoma are related to tumour progression,” Histopathology, vol. 69, no. 6, pp. 993–999, 2016. View at: Publisher Site | Google Scholar
  18. X. Liu, Q. Hong, Z. Wang, Y. Yu, X. Zou, and L. Xu, “Transforming growth factor-β-sphingosine kinase 1/S1P signaling upregulates microRNA-21 to promote fibrosis in renal tubular epithelial cells,” Experimental Biology and Medicine, vol. 241, no. 3, pp. 265–272, 2016. View at: Publisher Site | Google Scholar
  19. H. Li and E. P. Nord, “IL-8 amplifies CD40/CD154-mediated ICAM-1 production via the CXCR-1 receptor and p38-MAPK pathway in human renal proximal tubule cells,” American Journal of Physiology—Renal Physiology, vol. 296, no. 2, pp. F438–F445, 2009. View at: Publisher Site | Google Scholar
  20. Q. Lin, Y. Geng, M. Zhao, S. Lin, Q. Zhu, and Z. Tian, “MiR-21 regulates TNF-alpha-induced CD40 expression via the SIRT1-NF-kappaB pathway in renal inner medullary collecting duct cells,” Cellular Physiology and Biochemistry, vol. 41, no. 1, pp. 124–136, 2017. View at: Publisher Site | Google Scholar
  21. P. Pontrelli, M. Ursi, E. Ranieri et al., “CD40L proinflammatory and profibrotic effects on proximal tubular epithelial cells: Role of NF-κB and Lyn,” Journal of the American Society of Nephrology, vol. 17, no. 3, pp. 627–636, 2006. View at: Publisher Site | Google Scholar
  22. F. Fleissner, V. Jazbutyte, J. Fiedler et al., “Short communication: asymmetric dimethylarginine impairs angiogenic progenitor cell function in patients with coronary artery disease through a MicroRNA-21-dependent mechanism,” Circulation Research, vol. 107, no. 1, pp. 138–143, 2010. View at: Publisher Site | Google Scholar
  23. A. D. McClelland, M. Herman-Edelstein, R. Komers et al., “miR-21 promotes renal fibrosis in diabetic nephropathy by targeting PTEN and SMAD7,” Clinical Science, vol. 129, no. 12, pp. 1237–1249, 2016. View at: Publisher Site | Google Scholar
  24. S. R. Khan, “Hyperoxaluria-induced oxidative stress and antioxidants for renal protection,” Urolithiasis, vol. 33, no. 5, pp. 349–357, 2005. View at: Publisher Site | Google Scholar
  25. S. R. Khan, “Reactive oxygen species, inflammation and calcium oxalate nephrolithiasis,” Translational Andrology and Urology, vol. 3, no. 3, pp. 256–276, 2014. View at: Google Scholar
  26. W. Dröge, “Free radicals in the physiological control of cell function,” Physiological Reviews, vol. 82, no. 1, pp. 47–95, 2002. View at: Publisher Site | Google Scholar
  27. H. S. Huang, C. F. Chen, C. T. Chien, and J. Chen, “Possible biphasic changes of free radicals in ethylene glycol-induced nephrolithiasis in rats,” BJU International, vol. 85, no. 9, pp. 1143–1149, 2000. View at: Publisher Site | Google Scholar
  28. S. R. Khan, “Reactive oxygen species as the molecular modulators of calcium oxalate kidney stone formation: evidence from clinical and experimental investigations,” The Journal of Urology, vol. 189, no. 3, pp. 803–811, 2013. View at: Publisher Site | Google Scholar
  29. T. Umekawa, H. Tsuji, H. Uemura, and S. R. Khan, “Superoxide from NADPH oxidase as second messenger for the expression of osteopontin and monocyte chemoattractant protein-1 in renal epithelial cells exposed to calcium oxalate crystals,” BJU International, vol. 104, no. 1, pp. 115–120, 2009. View at: Publisher Site | Google Scholar
  30. M. Asselman, A. Verhulst, M. E. De Broe, and C. F. Verkoelen, “Calcium oxalate crystal adherence to hyaluronan-, osteopontin-, and CD44-expressing injured/regenerating tubular epithelial cells in rat kidneys,” Journal of the American Society of Nephrology, vol. 14, no. 12, pp. 3155–3166, 2003. View at: Publisher Site | Google Scholar
  31. S. Kusmartsev, P. R. Dominguez-Gutierrez, B. K. Canales, V. G. Bird, J. Vieweg, and S. R. Khan, “Calcium oxalate stone fragment and crystal phagocytosis by human macrophages,” The Journal of Urology, vol. 195, no. 4, part 1, pp. 1143–1151, 2016. View at: Google Scholar
  32. A. S. Awad, G. R. Kinsey, K. Khutsishvili, T. Gao, W. K. Bolton, and M. D. Okusa, “Monocyte/macrophage chemokine receptor CCR2 mediates diabetic renal injury,” American Journal of Physiology—Renal Physiology, vol. 301, no. 6, pp. F1358–F1366, 2011. View at: Publisher Site | Google Scholar
  33. H. Dimke, P. Desai, J. Borovac, A. Lau, W. Pan, and R. T. Alexander, “Activation of the Ca2+-sensing receptor increases renal claudin-14 expression and urinary Ca2+ excretion,” American Journal of Physiology—Renal Physiology, vol. 304, no. 6, pp. F761–F769, 2013. View at: Publisher Site | Google Scholar
  34. Y. Gong, N. Himmerkus, A. Plain, M. Bleich, and J. Hou, “Epigenetic regulation of microRNAs controlling CLDN14 expression as a mechanism for renal calcium handling,” Journal of the American Society of Nephrology, vol. 26, no. 3, pp. 663–676, 2015. View at: Publisher Site | Google Scholar
  35. Y. Gong, V. Renigunta, N. Himmerkus et al., “Claudin-14 regulates renal Ca++ transport in response to CaSR signalling via a novel microRNA pathway,” EMBO Journal, vol. 31, no. 8, pp. 1999–2012, 2012. View at: Publisher Site | Google Scholar
  36. A. Pushkin and I. Kurtz, “SLC4 base (HCO3, CO32−) transporters: classification, function, structure, genetic diseases, and knockout models,” American Journal of Physiology—Renal Physiology, vol. 290, no. 3, pp. F580–F599, 2006. View at: Publisher Site | Google Scholar
  37. S. L. Alper, “Molecular physiology of SLC4 anion exchangers,” Experimental Physiology, vol. 91, no. 1, pp. 153–161, 2006. View at: Publisher Site | Google Scholar
  38. G. Gambaro, F. Marchini, A. Piccoli, M. A. Nassuato, F. Bilora, and B. Baggio, “The abnormal red-cell oxalate transport is a risk factor for idiopathic calcium nephrolithiasis: a prospective study,” Journal of the American Society of Nephrology, vol. 7, no. 4, pp. 608–612, 1996. View at: Google Scholar
  39. D. He, Y. Lu, H. Hu et al., “The wnt11 signaling pathway in potential cellular EMT and osteochondral differentiation progression in nephrolithiasis formation,” International Journal of Molecular Sciences, vol. 16, no. 7, pp. 16313–16329, 2015. View at: Publisher Site | Google Scholar

Copyright © 2017 Chuangxin Lan et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


More related articles

763 Views | 347 Downloads | 2 Citations
 PDF  Download Citation  Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

We are committed to sharing findings related to COVID-19 as quickly and safely as possible. Any author submitting a COVID-19 paper should notify us at help@hindawi.com to ensure their research is fast-tracked and made available on a preprint server as soon as possible. We will be providing unlimited waivers of publication charges for accepted articles related to COVID-19. Sign up here as a reviewer to help fast-track new submissions.