Identification of tubular injury microRNA biomarkers in urine: comparison of next-generation sequencing and qPCR-based profiling platforms
© Nassirpour et al.; licensee BioMed Central Ltd. 2014
Received: 1 December 2013
Accepted: 9 June 2014
Published: 18 June 2014
MicroRNAs (miRNAs) are small, non-coding RNAs that regulate protein levels post-transcriptionally. miRNAs play important regulatory roles in many cellular processes and have been implicated in several diseases. Recent studies have reported significant levels of miRNAs in a variety of body fluids, raising the possibility that miRNAs could serve as useful biomarkers. Next-generation sequencing (NGS) is increasingly employed in biomedical investigations. Although concordance between this platform and qRT-PCR based assays has been reported in high quality specimens, information is lacking on comparisons in biofluids especially urine. Here we describe the changes in miRNA expression patterns in a rodent model of renal tubular injury (gentamicin). Our aim is to compare RNA sequencing and qPCR based miRNA profiling in urine specimen from control and rats with confirmed tubular injury.
Our preliminary examination of the concordance between miRNA-seq and qRT-PCR in urine specimen suggests minimal agreement between platforms probably due to the differences in sensitivity. Our results suggest that although miRNA-seq has superior specificity, it may not detect low abundant miRNAs in urine samples. Specifically, miRNA-seq did not detect some sequences which were identified by qRT-PCR. On the other hand, the qRT-PCR analysis was not able to detect the miRNA isoforms, which made up the majority of miRNA changes detected by NGS.
To our knowledge, this is the first time that miRNA profiling platforms including NGS have been compared in urine specimen. miRNAs identified by both platforms, let-7d, miR-203, and miR-320, may potentially serve as promising novel urinary biomarkers for drug induced renal tubular epithelial injury.
KeywordsNext-generation sequencing NGS Small RNA TruSeq qRT-PCR Urine miRNA Gentamicin Kidney injury
There currently are no biomarkers that can accurately track, correlate, or predict the occurrence of renal injury, particularly when related to drug-related toxicity [1–3]. There is a growing interest in the role of micro Ribonucleic Acids (miRNAs) in the pathogenesis of renal diseases [4, 5] and a growing number of investigations are directed to characterizing their potential as biomarkers [6–9]. miRNAs are highly conserved, endogenous, small (19–25 nucleotides), non-coding RNAs, with the primary role implicated in post-transcriptional silencing. Since a single miRNA can regulate multiple transcripts, miRNAs can potentially impact multiple signaling pathways and their dysregulation may lead to a variety of different disorders including renal disease . In fact, several miRNAs have been implicated in various renal diseases, including diabetic nephropathy, hypertension, glomerulonephritis, renal cancer, and polycystic kidney disease [11–17]. miRNAs are secreted into the circulation by many different mechanisms including various lipid-containing vesicles, such as exosomes, microvesicles, and apoptotic bodies, but can also be protected from degradation by binding to RNA-binding proteins . Due to their stability, miRNAs are readily quantified in serum, plasma and other body fluids, and thus more investigations are being directed to characterizing their potential as biomarkers. In fact, study results from our laboratory and other participating members of the Health and Environmental Sciences Institute’s (HESI) Nephrotoxicity Committee have identified changes in miRNA levels in various kidney disease animal models. Although quantitative real-time PCR (qRT-PCR) has most commonly been used to profile annotated miRNome, newer technologies which promise increased sensitivity and specificity are becoming available at astonishing rates. For example, next-generation sequencing (NGS) is gaining popularity and has successfully been used to characterize miRNA profiles in various tissues  as well as bio-fluids including blood and, most recently, cerebral spinal fluid [19, 20]. Although NGS platform has been shown to generate highly reproducible, accurate data with high correlation to other platforms for high quality RNA samples [21–24], comparative data on media with low quality and quantity RNA are lacking.
The objective of this study was to characterize the miRNA expression changes using two different profiling platforms in urine specimens following induction of toxic renal tubular injury in rodents with gentamicin. Gentamicin is an aminoglycoside antibiotics which like other nephrotoxicants such as mercuric chloride and chromium induces injury of renal tubular epithelial cells . Gentamicin has also been reported to induce more severe nephrotoxicity than other aminoglycoside antibiotics when administered at high doses [26, 27] and has already been shown to induce miRNA changes in rats, albeit at a different dose and treatment regimen . Accordingly, we performed a global urinary miRNA expression profiling using both the more conventional qRT-PCR TaqMan Low Density Array (TLDA-A, Life Technologies) method as well as NGS in control and treated rats. The presence of tubular epithelial injury was confirmed by histopathology and assessment of a panel of urine protein biomarkers that have been qualified for use in non-clinical drug development : kidney injury molecule-1 (KIM-1), total protein, and beta 2- microglobulin.
Although a recent paper has described NGS characterization of exosomal and non-exosomal miRNAs in human urine , to our knowledge, this is the first time that the NGS and qRT-PCR based platforms have been compared in urine. We report that although both platforms were able to identify drug induced changes in miRNA expression, that there is discordance in the level of sensitivity of detection. Specifically, with exception of three miRNAs, miRNA-seq did not detect some miRNAs which were identified by qRT-PCR with threshold cycle (Ct) values lower than 32 in every urine specimen analyzed. These differences could be due to the level of sensitivity posed by NGS analysis for low abundant miRNAs which, are easily detected by the qRT-PCR method (performed with pre-amplification). On the other hand, the vast majority of miRNAs detected by NGS were in their isomeric forms, and thus could not have been detected by the pre-designed qRT-PCR primers available on the TLDA-A. Our results suggest that although miRNA-seq is specific, with few false positive calls, it may not detect certain low abundant miRNAs in low quality RNA specimen like those currently obtained from urine. However, the qRT-PCR analysis will miss a wealth of information, most importantly in the form of isoforms. In summary, we have compared the detection of urinary miRNA alterations in tubular injury with qRT-PCR and miRNA-seq. Each platform has its strengths and weaknesses and therefore, careful considerations should be made when selecting a platform of choice for urinary miRNA profiling.
Histopathology and urinary protein biomarker data indicate gentamicin induced tubular injury at D7 with 50 mg/kg
Levels of miRNAs in urine significantly changed following gentamicin induced tubule injury
Of the 173 miRNAs identified in urine specimens by qRT-PCR, 32 miRNAs changed significantly due to treatment with p-values < 0.05 and ≥ ±1.5 fold-change
Detector name (ABI)
miRNA ID v. 20
Small RNA Next Generation Sequencing can be used to evaluate miRNA expression in urinary samples
Urine samples from control (samples 1–7) and gentamicin treated rats (samples 32, 33, 39) were used for RNA isolation using the miRNeasy (Qiagen) extraction kit and analyzed for quality and quantity before NGS analysis
Percentages of miRNAs detected in control (7) and gentamicin-treated (3) specimen
Since a large proportion of reads (14-60%) remained unaligned to rat (Figure 4A), we investigated the origin of highly represented sequences in the reads that remained unaligned. In 2 of the 3 gentamicin treated specimen, a single sequence comprised over 40% of the unaligned reads. One sample was highly contaminated by a PCR primer and the other sample was contaminated with a sequence from the zebrafish. Similarly, in control samples we observed contamination from a variety of bacterial and fungal genomes (e.g. from Saccharomyces cerevisiae, Paenibacillus polymyxa, Ethanoligenens harbinense, Roseburia hominis, to name a few). This may indicate fungal and bacterial contamination during sample collection or alternatively these contaminations may have been inadvertently introduced during sample processing procedures. To what extent these contaminations affected the NGS results needs to be determined although among control samples, sample 4 had the least amount of contamination and detected the maximum number (111) of distinct miRNA species, whereas sample 2 that had the most contamination contributed 76 miRNAs.
We also investigated the possibility that some of these unaligned reads may have a dietary origin and thus performed additional mapping to the wheat, corn and soybean genomes. We observed that a proportion of our unaligned reads did originate from dietary genomes from wheat (2%-52%, mean = 19%), corn (0.1%-3.8%, mean = 1.2%) and soybean (0.1%-15%, mean = 3.8%).
Next we investigated various normalization strategies. Recent literature suggest that miRNA-Seq has its own challenges and applying the normalization strategies commonly used in mRNA-seq may need further refinement for use on miRNA sequencing data [31, 32]. We thus compared three different normalization strategies and three differential analysis methods. DESeq and EdgeR were applied to their respective normalizations and Limma was applied to RPM normalized data. EdgeR normalized values gave the best segregation and gave the maximum number of overlap with the differentially regulated miRNAs identified by qRT-PCR and hence are reported here.
Of the 146 miRNAs identified in urine specimens using NGS analysis, 14 had p-values < 0.05 and ≥ ±1.5 fold-change
Prediction of miRNA targets and pathway analysis relevant to renal injury
This study builds on a growing data set indicating that miRNAs in body fluids can be used as non-invasive biomarkers for detecting and monitoring various pathophysiological conditions. Gentamicin, an aminoglycoside antibiotic used to treat many types of gram negative bacterial infections, is a known nephrotoxicant. The underlying mechanism of gentamicin renal toxicity is assumed to be mainly related to the formation of free hydroxyl radicals and the inhibition of mitochondrial respiration, although other mechanisms have also been suggested [33–36]. The microscopic findings in the kidneys observed in this study were similar to those previously documented following administration of gentamicin to rats [37–39]. Therefore, in this study, gentamicin at a dose of 50 mg/kg/day successfully induced tubular injury 7 days post treatment (Figure 1). The objective of this study was to induce nephrotoxicity in rats with gentamicin and evaluate two miRNA profiling platforms in detecting significant changes in urine samples.
Sequencing allows for complete evaluation of all small RNAs, without primer bias or a need for species annotation and as NGS analysis becomes more affordable, more scientists are choosing this platform in their investigations. Several studies have examined the concordance between various profiling platforms including quantitative RT-PCR and miRNA-seq and although it is generally believed that for high quality specimens intra-platform variability is usually minimal, cross-platform comparisons do not necessarily correlate well [21–24]. Even less is known on comparisons in biofluids, especially urine. To address this issue, we have analyzed RNA derived from urine specimen using Life technologies’ TLDA-A and Illumina’s NGS platforms. We thus prepared small RNA sequencing libraries using Illumina’s TruSeq sample preparation kit from urinary total RNA extracted and used on the qRT-PCR analysis, and sequenced the samples on the HiSeq 2000 (Figure 3). Utilizing NGS allowed us to identify several kidney-related miRNA which were significantly changed in rats exposed to the kidney toxicant gentamicin. However, the majority of miRNAs identified by NGS were not the same as those identified on TLDA-A cards (Tables 1 and 4). Thus, we report for the first time, that there is a minimal agreement between NGS and qRT-PCR for low yield urinary miRNA analysis.
The observed minimal agreement between the two platforms could be due to a variety of reasons. One reason could be the fact that we only analyzed miRNAs available on the TLDA-A, which does not include most of the miRNAs detected by NGS. TLDA allows for quantitation of 335 mouse and 226 rat miRNAs based on Sanger mirBase v. 10 annotations. Out of 375 miRNA sequences on the TLDA-A, 218 had a corresponding identifier in Rat mirBase v. 20; 147 in mouse and 1 in human. 5 TLDA-A miRNAs could not be found in mirBase v. 20 as they were deprecated. Out of the 14 miRNAs that were detected to be differentially expressed using NGS (Table 4), five (rno-let-7d-5p, rno-miR-100-5p, rno-miR-203a-3p, rno-miR-21-5p, rno-miR-320-3p) were represented on the TLDA-A, while nine (rno-miR-378a-3p, mmu-miR-5100, rno-miR-30e-3p, rno-miR-125b-2-3p, rno-miR-320-5p, rno-miR-3473, rno-miR-21-3p, rno-miR-455-5p, and hsa-miR-7641) were not. In fact, these 9 miRNAs are in their isomeric forms (Figure 4B). This finding highlights a key limitation in the qRT-PCR analysis as the TLDA primers may only detect some isomiRs with extra bases on either ends, but will miss the isomers with deletions/substitutions. Hence, only a subset of miRNAs are queried on the TLDA-A, and as our results indicate, nine miRNAs differentially regulated in NGS could not be detected by the TLDA-A due to its dependency of primer design. Employing the full TLDA array capabilities to include both A and B cards may improve concordance. The NGS platform, on the other hand, seems to be highly sensitive to the presence of contaminants. As observed in this study, contaminations ranging from PCR primers to bacterial and fungal genomes could negatively affect detectability of low expressing miRNAs. Hence the sample with the least amount of contamination contributed the maximum species of miRNAs compared to ones with higher contamination. Additionally, as depicted in Table 3, out of the total 132 and 168 rat-miRNAs detected in TLDA-A and NGS control samples respectively, miRNAs in TLDA-A seem to be far more consistently detected than in NGS.Although the absence of primer-dimer sequences indicated that there was enough genetic material to sequence, a large proportion of the reads aligned to non-coding RNAs (introns, intergenic regions, and exons). Furthermore, a large proportion (20-60%) of reads could not be aligned to the rat genome, which may indicate the presence of exogenous RNAs, mRNA degradation, genetic material expelled through urine, and other contaminants mainly from dietary, bacterial and fungal origins (Figure 4). For example, most of the miRNAs that were mapped to mouse (25) and human (2) can be traced to 3’ regions of protein coding regions of rat. Furthermore, about 0.6-16% of detected miRNAs belonged to human and mouse. 24 out of 35 mouse miRNAs and 8 out of 9 human miRNAs detected in the samples were found in unannotated regions of rat genome. These findings, together with the fact that the number of miRNAs annotated to the rat is only about a third compared to human and half as many as in mouse strongly suggests that the rat miRnome may be under-annotated.
Another reason for the low concordance between the 2 platforms may be their differences in sensitivity. Specifically, following reverse transcription of miRNA targets, we performed a pre-amplification step to increase our detection capabilities for the low expressing miRNAs, which is a step not currently available for the miRNA NGS analysis. The use of cDNA pre-amplification has been proposed to increase the sensitivity of miRNA detection, although with increased variability for low abundant miRNAs [40, 41]. However, all three miRNAs: let-7d, miR-203, and miR-320 were reliably detected by qRT-PCR in every sample with low mean Ct values (miR-let-7d: 27.52 ± 0.77; miR-203: 24.66 ± 0.3; miR-320: 24.69 ± 1.1) which suggests relatively high levels of urinary expression. Therefore, it is likely that although the qRT-PCR method will yield a higher level of sensitivity (especially when combined with the pre-amplification step as was performed in this study), the NGS may actually provide higher specificity. Specifically, miRNA-seq did not detect some miRNAs which qRT-PCR detected (<32 Ct). In our study we were able to evaluate the NGS samples to a depth of 10 million reads. It would thus be interesting to examine if increasing the sequencing depth could improve detection sensitivity and concordance between the technologies.
These results suggest that although miRNA-seq is specific, due to its currently stringent RNA input requirements, it may not be able to detect certain low quantity (50 ng) and low quality (OD 260/280 < 1.8; RIN ~ 2) RNAs in urine at the conditions used in this study (50 base single end sequencing totaling 10 million reads). Additionally, although total RNA was extracted for each sample and divided between the TLDA and NGS analysis, the RNA was further enriched for small RNAs during the Illumina Truseq small RNA library preparation. Just as the pre-amplification step in the TLDA platform may introduce bias, steps such as RNA ligation may introduce technical bias during the NGS library construction . None the less, as reported previously, urine in general contains very little RNA , and this low quality and quantity of the RNA (as shown in Table 2) might be key to the low observed concordance between these platforms. Efforts to improve the quality of RNA to be extracted from urine specimen may decrease the inter-platform variability observed. Recently, as more commercially available kits are designed for this purpose, more investigators are evaluating the quality as well as quantity of miRNAs extracted from urine . These investigations together with improved methodology for enriching of small RNAs to be used in profiling platforms hold promise for improved concordance. Furthermore, although in this study we used cell-free urine probably containing both exosomal as well as non-exosomal miRNAs, enriching for exosomal RNAs may also shed light on differences in analysis between profiling platforms .Therefore, our data suggest that each platform has its own advantages and disadvantages and their utility should be decided based on the experimental objectives and study design limitations. Regardless of platform differences, the three miRNAs identified by both platforms, let-7d, miR-203 and miR-320, hold promise as tubular injury biomarkers due to their strong functional associations with cellular process or pathways relevant to renal disease. In fact, the main toxicology related function predicted by IPA was ‘Renal Nephritis’ (Figure 5A). Looking at each platform independently, it seems like the differentially regulated miRNAs identified by TLDA-A, namely miR-16-5p, 185-5p, 200b-3p, 208a-3p, 23a-3p, 30c-3p, 3118, 494-3p and 92a-3p which are associated with nephritis (from Ingenuity Pathway Analysis) may be promising biomarkers of tubule injury. Similarly, miR-100-5p and let-7a-5p implicated in renal inflammation and nephritis may be tubular injury biomarkers specifically detected by NGS. Taken together, these results suggest that miRNAs isolated from urine and profiled by either NGS or qRT-PCR could potentially serve as translational biomarkers for detection of some forms of renal injury.
Drug induced kidney injury is not only an unfortunate clinical manifestation but it also hinders pharmaceutical progress. Urine samples have the potential to provide a wealth of information both for clinical applications as well as translational investigators, and provide a perfect medium for biomarkers of injury. To our knowledge, this is the first study comparing urinary miRNA expression data generated by qRT-PCR to those from sequencing. Although, the high cost of these experiments prohibited us from performing technical replicates and limited our sample size (7 control and 3 treated), we were able to evaluate miRNA expression analysis employing two very different technologies and show that, for at least three miRNAs, that there was high inter-platform concordance.
In summary, we have shown that differentially expressed urinary miRNAs can be detected by two distinct profiling platforms, qRT-PCR and miRNA-seq. Our data indicate that the three miRNAs detected by both technologies may be useful urinary biomarkers for tubular injury.
All animal experiments were conducted in compliance with the UK Animals (Scientific Procedures) Act 1986 and its associated Guidelines and Codes of Practice. Studies were reviewed internally by the Pfizer, Sandwich Ethical Review Process and subsequently authorized by the UK Home Office as Project Licence 70/6903. Male Sprague Dawley rats (225-270 g, Charles River Laboratories) were maintained in a central animal facility housed individually in polycarbonate cages with autoclaved woodchip bedding (Lillico) enriched with Lillico paperwool (nesting material). The room environment was maintained at 21°C ± 2°C and 55 ± 10% relative humidity at all times in an alternating 12-h light–dark cycle. Rats were removed from cages during surgery, toxicant injection, and blood and urine collection. Animals were acclimated to the laboratory environment for a minimum of 5 days prior to initiation of dosing. Water (purified by reverse osmosis) and Certified Rodent Diet 5002 (PMI Feeds, Inc.) were provided ad libitum. Animals were euthanized by isoflurane gas anesthesia followed by exsanguinations.
Gentamicin Sulphate (Sigma) , in a sterile solution of 0.9% saline, pH to 7 (adjusted with 1 M NaOH), was administered by subcutaneous injection once daily to male rats at doses of 0, 10, 25 and 50 mg/kg/day for 7 days (maximum of 1 mL per site and no more than 4 sites in any 24 hour period).
Animals were placed into chilled metabolic cages for urine collection overnight. All animals were fasted, but provided free access to water. Total urine volume for each animal was recorded. Urine was briefly centrifuged (1000 g for 5 min at 4°C) to sediment any contaminants and cellular debris and was aliquoted and frozen at -80°C. Each urine protein biomarker was normalized to the total amount excreted in the volume of urine collected, and to the urine creatinine concentration (mg/ml) to normalize for dilutional effects. There was no change in the amount of creatinine excreted in the urine of control or treated animals over the time course of the studies (data not shown). KIM-1 and Beta-2-Microglobulin were measured using MSD technology, and urine protein was analyzed using Advia 1650 automated technology. Data are expressed as mean ± standard error. Statistical difference (p < 0.05) was calculated by one-way ANOVA. p < 0.05 is considered significant and is represented by a single asterisk (*) where applicable. All graphs were generated by Prism software (GraphPad Software).
Light microscopic examination
Kidneys were weighed at scheduled necropsy, preserved in 10% neutral buffered formalin, then processed through graded alcohol and xylene, infiltrated and embedded in paraffin, sectioned, and stained with hematoxylin and eosin for histological assessment. Sections for microscopic evaluation included one coronal section through the middle third of one kidney and one mid-sagittal section through the middle third of the other kidney. Histological sections for both kidneys included cortex, medulla and pelvis. The agreed upon standardized HESI/PSTC histology lexicon were used. Histopathology evaluation was performed by a board certified veterinary pathologist who had knowledge of the treatment groups and necropsy data (organ weights and macroscopic observations) but who was blinded as to clinical pathology datasets, including results from biomarker evaluations.
RNA extraction and quantitative Polymerase Chain Reaction (qPCR)
Total RNA was isolated from 200 μL of rat urine using a modified protocol (miRNeasy RNA extraction kit, Qiagen) according to the manufacturer’s instructions. Briefly, 700 μL of QIAzol reagent was added to 200 μL of urine. After vortexing vigorously with chloroform, the samples were then centrifuged at 12,000 g for 15 min at 4°C. The upper aqueous phase was transferred to a new tube and 1.5 volume of ethanol was added. The sample was then applied to the column and washed. The immobilized RNA was then collected from the membrane with 30 μL of RNase free water.
Total RNA concentration was measured at 260 nm using a NanoDrop 2000c spectrophotometer (Thermo Scientific). RNA concentration was further measured for urine RNA samples on a RNA 6000 Nano chip using Agilent 2100 Bioanalyzer (Agilent Technologies). This RNA was used for both the NGS as well the PCR assays. Quantitative miRNA analysis was performed using TaqMan miRNA assays from Applied Biosystems. RNA was reverse transcribed into complementary DNA (cDNA) using megaplex primers (16°C for 30 min; 42°C for 30 min; 85°C for 5 min). The cDNA product was then used in a pre-amplification step. Quantitative PCR was performed using TaqMan Low Density Array (TLDA A, Life Technologies) in a Viia7 instrument with the following temperature profile: 95°C for 10 min followed by 40 cycles of 95°C for 15 s and 60°C for 1 min. Only miRNAs that were detected with a Ct value of 32 or less and in at least 20% of all the samples were included in the analysis. We investigated several normalization methods, including invariant miRNAs, global mean, and lowess using the Expressionist version 7.6 (Genedata). Welch test was applied to each normalized data set to assess differential regulation between control and treated samples. All three normalization methods yielded similar results in terms of statistical significance (P-value < 0.05) and fold changes (FC magnitude > 1.5). However, we report the results derived from the lowess normalization as it produced the best match to the NGS data. All graphs were made by Spotfire DecisionSite (Spotfire Inc.). The miRNAs from TLDA-A, which are based on mirBase v.10, were mapped to mirBase v. 20 using their corresponding v.10 sequences. This allowed us to focus our efforts on the v. 20 miRNA sequences in our comparative analysis between TLDA-A and NGS.
miRNA reads were obtained from BGI. The 3’ and 5’ adaptors were trimmed off the reads and sequences between 16 and 35 nucleotides of average phred score of ≥ 30 were retained for further analysis. The reads were first mapped to rat miRNAs using bowtie  without allowing for mismatches. Instead of discarding the remaining reads, we explored the possibilities of RNA editing to include potential isomiRs . The reads were thus mapped to miRNAs that have at most 1 base mismatch, a substitution at the 5’ or 3’ ends, or have 1 or 2 extra bases at either ends. Although 1 base mismatch would include 1 base substitution/deletion/addition at the ends, we made sure that 1 base mismatch was inside the miRNA and that all categories were mutually exclusive. To account for heterogeneity, reads were further mapped to the precursor miRNAs. Reads that uniquely mapped to miRNAs were considered and were profiled as: 1) identical to miRNAs, 2) perfectly match to miRNAs but are of shorter length (no more than 5 bases short), 3) mapping to miRNAs with 1 substitution, 1 or 2 extra bases at 3’ and 5’ ends and with those with at most 1 mismatch (occurring inside the sequence and not on ends), 4) mapping to the precursor miRNA. Thus, steps 2–4 constitute “IsomiRs”-alternative forms of the mature miRNAs. We merged the counts of reads identical to a miR and those mapping to Isomeric forms. Additionally, since the total number of annotated mature miRNAs in rat is 728 compared to 1908 in mouse and 2578 in human, we accounted for the possibility of miRNA under-annotation in rat by expanding our mapping to include miRNAs from mouse and human. Thus, after aligning to rat, the remaining reads were further mapped to mouse and eventually to human miRNAs. MirBase v. 20 was used for miRNA sequence mapping. Since urine samples may contain exogenous RNAs and other small RNAs, reads were also mapped to the RFAM (RNA Families) database  and exon sequences in the rat. The other common contaminant in miRNA-Seq is primer-dimers which can adversely affect the quantification of miRNA abundance. Specifically, the Illumina protocol consists of manual gel band excision, which enhances the occurrence of primer-dimers due to incorrect cutting of the gel or insufficient sample to sequence. We thus included extensive primer-dimer analysis in our study. Furthermore, since normalization and differential expression analysis in mRNA-Seq are still being finalized by the scientific community, we employed three different widely used methods in our analysis. We first obtained a Reads per Million (RPM) count for a miRNA by dividing the number of reads mapping to that miRNA by the total number of reads mapping to all miRNAs and multiplying that by 1,000,000. This gave the abundance of a miRNA in a million reads. Additionally, the counts were normalized using scaling factors assuming a negative binomial distribution using DESeq (Differential Expression for Sequence Count Data)  and trimmed mean normalization using EdgeR (Empirical analysis of Digital Gene Expression in R)  packages in R-Bioconductor. We next used various filtering criteria and hierarchical clustering to achieve a clear segregation between the control and treated samples. Spearman correlation coefficient was used on the normalized values and a “complete” hierarchical clustering was then performed. To find the optimum differential analysis we investigated three different methods : Limma, variants of Fisher’s exact test from DESeq, and EdgeR packages from the Bioconductor software. A p-value cutoff of 0.05 was used for significance detection. In addition, to check for sequences that are of dietary origin, all unaligned sequences were aligned against wheat (Triticum aestivum), corn (Zea mays) and soybean (Glycine max). Sequences were first mapped to wheat, then the unaligned sequences were mapped to corn and then to soybean. The miRNAs from human and mouse that were detected in the rat samples were blasted against the rat database in NCBI. Since it also consisted of miRNAs in only isomeric forms, blast results that had >85% of the sequence aligned were considered for further analysis.
Prediction of toxicity effects for NGS/TaqMan significant miRNAs
The 14 miRNAs identified as significant by NGS and the 32 miRNAs identified by TaqMan were analyzed separately in Ingenuity System’s Ingenuity Pathway Analysis IPA tool. IPA Core Analysis was performed on each set and the most significant ‘Tox Functions’ along with relevant biological functions are reported. Experimentally verified gene targets were then further evaluated using IPA’s ‘Core Analysis to extract over-represented biological pathways with p-value < 0.01. Network analysis was performed on the three conserved miRNAs (let7a-5p, miR-203-3p, and miR-320b) using IPA’s ‘Grow’ function to add genes experimentally linked downstream of the miRNAs via ‘RNA/RNA interactions analysis: miRNA targeting’ relationship. The resulting network was expanded using the ‘Overlay’ function with IPA’s ‘Ingenuity Tox List’. The most overlapping list ‘Renal Necrosis/Cell Death’ was added and the unlinked genes trimmed off for better visualization.
We would like to thank Dr. Dale Morris, Dr. Denise Robinson-Gravatt and Pfizer’s Science and Technology Board, Dr. Magali Guffroy, Steven Bailey, Dr. Michael Agostino and Donald Fish for their generous help, support, and commitment. This study was conducted within the HESI Biomarkers of Nephrotoxicity Committee. HESI is a collaborative consortium consisting of pharmaceutical companies, academic institutions, and regulatory agencies with a common goal of attaining a deeper insight into the molecular mechanisms of kidney toxicity. The consortium aims to identify early preclinical safety biomarkers and we would like to especially acknowledge Dr. Jean-Charles Gautier, Jean-Francois Leonard and Dr. Raegan O'Lone for their critical reading of our manuscript.
- Vaidya VS, Ferguson MA, Bonventre JV: Biomarkers of acute kidney injury. Annu Rev Pharmacol Toxicol. 2008, 48: 463-493.PubMed CentralPubMedView ArticleGoogle Scholar
- Vaidya VS, Waikar SS, Ferguson MA, Collings FB, Sunderland K, Gioules C, Bradwin G, Matsouaka R, Betensky RA, Curhan GC, Bonventre JV: Urinary biomarkers for sensitive and specific detection of acute kidney injury in humans. Clin Translat Sci. 2008, 1 (3): 200-208.View ArticleGoogle Scholar
- Sasaki D, Yamada A, Umeno H, Kurihara H, Nakatsuji S, Fujihira S, Tsubota K, Ono M, Moriguch A, Watanabe K, Seki J: Comparison of the course of biomarker changes and kidney injury in a rat model of drug-induced acute kidney injury. Biomarkers. 2011, 16 (7): 553-566.PubMedView ArticleGoogle Scholar
- Khella HWZ, Bakhet M, Lichner Z, Romaschin AD, Jewett MAS, Yousef GM: MicroRNAs in kidney disease: an emerging understanding. Am J Kidney Dis. 2013, 61 (5): 798-808.PubMedView ArticleGoogle Scholar
- Saal S, Harvey SJ: MicroRNAs and the kidney: coming of age. Curr Opin Nephrol Hyperten. 2009, 18 (4): 317-323.View ArticleGoogle Scholar
- Akkina S, Becker BN: MicroRNAs in kidney function and disease. Translat Res. 2011, 157 (4): 236-240.View ArticleGoogle Scholar
- Neal CS, Michael MZ, Pimlott LK, Yong TY, Li JYZ, Gleadle JM: Circulating microRNA expression is reduced in chronic kidney disease. Nephrol Dialysis Transplant. 2011, 26 (11): 3794-3802.View ArticleGoogle Scholar
- White NM, Bao TT, Grigull J: miRNA profiling for clear cell renal cell carcinoma: biomarker discovery and identification of potential controls and consequences of miRNA dysregulation. J Urol. 2011, 186: 1077-1083.PubMedView ArticleGoogle Scholar
- Lorenzen JM, Thum T: Circulating and Urinary microRNAs in Kidney Disease. Clin J Am Soc Nephrol. 2012, 7 (9): 1528-1533.PubMedView ArticleGoogle Scholar
- Kasinath BS, Feliers D: The complex world of kidney microRNAs. Kidney Int. 2011, 80: 334-337.PubMedView ArticleGoogle Scholar
- White NM, Yousef GM: MicroRNAs: exploring a new dimension in the pathogenesis of kidney cancer. BMC Med. 2010, 8: 65-PubMed CentralPubMedView ArticleGoogle Scholar
- Yi Z, Fu Y, Zhao S: Differential expression of miRNA patterns in renal cell carcinoma and nontumorous tissues. J Cancer Res Clin Oncol. 2010, 136: 855-862.PubMedView ArticleGoogle Scholar
- Zampetaki A, Kiechl S, Drozdov I, Willeit P, Mayr U, Prokopi M, Mayr A, Weger S, Oberhollenzer F, Bonora E, Shah A, Willeit J, Mayr M: Plasma MicroRNA Profiling Reveals Loss of Endothelial MiR-126 and Other MicroRNAs in Type 2 Diabetes. Circulation Research. 2010, 107 (6): 810-817.PubMedView ArticleGoogle Scholar
- Pandey P, Qin S, Ho J: Systems biology approach to identify transcriptome reprogramming and candidate microRNA targets during the progression of polycystic kidney disease. BMC Syst Biol. 2011, 5: 56-PubMed CentralPubMedView ArticleGoogle Scholar
- Redova M, Svoboda M, Slaby O: MicroRNAs and their target gene networks in renal cell carcinoma. Biochem Biophys Res Comm. 2011, 405 (2): 153-156.PubMedView ArticleGoogle Scholar
- Heggermont WA, Heymans S: MicroRNAs Are Involved in End-Organ Damage During Hypertension. Hypertension. 2012, 60 (5): 1088-1093.PubMedView ArticleGoogle Scholar
- Alvarez ML, DiStefano JK: The role of non-coding RNAs in diabetic nephropathy: potential applications as biomarkers for disease development and progression. Diabetes Res Clin Pract. 2013, 99 (1): 1-11.PubMedView ArticleGoogle Scholar
- Meng F, Hackenberg M, Li Z, Yan J, Chen T: Discovery of Novel MicroRNAs in Rat Kidney Using Next Generation Sequencing and Microarray Validation. PLoS One. 2012, 7 (3): e34394-PubMed CentralPubMedView ArticleGoogle Scholar
- Wu Q, Wang C, Lu Z, Guo L, Ge Q: Analysis of serum genome-wide microRNAs for breast cancer detection. Clinica Chimica Acta. 2012, 413 (13–14): 1058-1065.View ArticleGoogle Scholar
- Burgos KL, Javaherian A, Bomprezzi R, Ghaffari L, Rhodes S, Courtright A, Tembe W, Kim S, Metpally R, Van Keuren-Jensen K: Identification of extracellular miRNA in human cerebrospinal fluid by next-generation sequencing. RNA. 2013, 19 (5): 712-722.PubMed CentralPubMedView ArticleGoogle Scholar
- Tam S, de Borja R, Tsao M-S, McPherson JD: Robust global microRNA expression profiling using next-generation sequencing technologies. Lab Invest. 2014, 94 (1): 350-PubMedView ArticleGoogle Scholar
- Llorens F, Hummel M, Pantano L, Pastor X, Vivancos A, Castillo E, Mattlin H, Ferrer A, Ingham M, Noguera M, Kofler R, Dohm J, Pluvinet R, Bayes M, Himmelbauer H, del Rio J, Marti E, Sumoy L: Microarray and deep sequencing cross-platform analysis of the mirRNome and isomiR variation in response to epidermal growth factor. BMC Genomics. 2013, 14 (1): 371-PubMed CentralPubMedView ArticleGoogle Scholar
- Kolbert CP, Feddersen RM, Rakhshan F, Grill DE, Simon G, Middha S, Jang JS, Simon V, Schultz DA, Zschunke M, Lingle W, Carr JM, Thompson EA, Oberg AL, Eckloff BW, Wieben ED, Li P, Yang P, Jen J: Multi-Platform Analysis of MicroRNA Expression Measurements in RNA from Fresh Frozen and FFPE Tissues. PLoS One. 2013, 8 (1): e52517-PubMed CentralPubMedView ArticleGoogle Scholar
- Git A, Dvinge H, Salmon-Divon M, Osborne M, Kutter C, Hadfield J, Bertone P, Caldas C: Systematic comparison of microarray profiling, real-time PCR, and next-generation sequencing technologies for measuring differential microRNA expression. RNA. 2010, 16 (5): 991-1006.PubMed CentralPubMedView ArticleGoogle Scholar
- Sieber M, Hoffmann D, Adler M: Comparative analysis of novel noninvasive renal biomarkers and metabonomic changes in a rat model of gentamicin nephrotoxicity. Toxicol Sci. 2009, 109: 336-349.PubMedView ArticleGoogle Scholar
- Zhou Y, Vaidya VS, Brown RP, Zhang J, Rosenzweig BA, Thompson KL, Miller TJ, Bonventre JV, Goering PL: Comparison of Kidney Injury Molecule-1 and Other Nephrotoxicity Biomarkers in Urine and Kidney Following Acute Exposure to Gentamicin, Mercury, and Chromium. Toxicol Sci. 2008, 101 (1): 159-170.PubMed CentralPubMedView ArticleGoogle Scholar
- Hottendorf GH, Barnett D, Gordon LL, Christensen EF, Madissoo H: Nonparallel nephrotoxicity dose–response curves of aminoglycosides. Antimicrobial Agents and Chemotherapy. 1981, 19 (6): 1024-1028.PubMed CentralPubMedView ArticleGoogle Scholar
- Saikumar J, Hoffmann D, Kim T-M, Gonzalez VR, Zhang Q, Goering PL, Brown RP, Bijol V, Park PJ, Waikar SS, Vaidya VS: Expression, Circulation, and Excretion Profile of MicroRNA-21, -155, and -18a Following Acute Kidney Injury. Toxicological Sciences. 2012, 129 (2): 256-267.PubMed CentralPubMedView ArticleGoogle Scholar
- Warnock DG, Peck CC: A roadmap for biomarker qualification. Nat Biotech. 2010, 28 (5): 444-445.View ArticleGoogle Scholar
- Cheng L, Sun X, Scicluna BJ, Coleman BM, Hill AF: Characterization and deep sequencing analysis of exosomal and non-exosomal miRNA in human urine. Kidney Int. 2013, doi: 10.3402/jev.v2i0.19671Google Scholar
- Garmire LX, Subramaniam S: Evaluation of normalization methods in mammalian microRNA-Seq data. RNA. 2012, 18 (6): 1279-1288.PubMed CentralPubMedView ArticleGoogle Scholar
- Zhou X, Oshlack A, Robinson MD: miRNA-Seq normalization comparisons need improvement. RNA. 2013, 19 (6): 733-734.PubMed CentralPubMedView ArticleGoogle Scholar
- Sue YM, Cheng CF, Chang CC: Antioxidation and anti-inflammation by haem oxygenase-1 contribute to protection by tetramethylpyrazine against gentamicin-induced apoptosis in murine renal tubular cells. Nephrol Dial Transplant. 2009, 24: 769-777.PubMedView ArticleGoogle Scholar
- Walker PD, Shah SV: Evidence suggesting a role for hydroxyl radical in gentamicin-induced acute renal failure in rats. J Clin Invest. 1988, 81: 334-341.PubMed CentralPubMedView ArticleGoogle Scholar
- Walker PD, Barri Y, Shah SV: Oxidant mechanisms in gentamicin nephrotoxicity. Ren Fail. 1999, 21 (3–4): 433-442.PubMedView ArticleGoogle Scholar
- Lopez-Novoa JM, Quiros Y, Vicente L, Morales AI, Lopez-Hernandez FJ: New insights into the mechanism of aminoglycoside nephrotoxicity: an integrative point of view. Kidney Int. 2011, 79 (1): 33-45.PubMedView ArticleGoogle Scholar
- Houghton DC, Hartnett M, Campbell-Boswell M, Porter G, Bennett W: A light and electron microscopic analysis of gentamicin nephrotoxicity in rats. Am J Pathol. 1976, 82 (3): 589-612.PubMed CentralPubMedGoogle Scholar
- Mingeot-Leclercq MP, Brasseur R, Schanck A: Molecular parameters involved in aminoglycoside nephrotoxicity. J Toxicol Environ Health. 1995, 44: 263-300.PubMedView ArticleGoogle Scholar
- Mingeot-Leclercq M-P, Tulkens PM: Aminoglycosides: Nephrotoxicity. Antimicrob Agents Chemother. 1999, 43 (5): 1003-1012.PubMed CentralPubMedGoogle Scholar
- Mestdagh P, Feys T, Bernard N, Guenther S, Chen C, Speleman F, Vandesompele J: High-throughput stem-loop RT-qPCR miRNA expression profiling using minute amounts of input RNA. Nucleic Acids Res. 2008, 36 (21): e143-PubMed CentralPubMedView ArticleGoogle Scholar
- Chen Y, Gelfond J, McManus L, Shireman P: Reproducibility of quantitative RT-PCR array in miRNA expression profiling and comparison with microarray analysis. BMC Genomics. 2009, 10 (1): 1-10.View ArticleGoogle Scholar
- Karim S, Helio P, HA E, Ana K, Sam G-J, Vincent M, Tamas D: Reducing ligation bias of small RNAs in libraries for next generation sequencing. Silence. 2012, 3 (1): 4-View ArticleGoogle Scholar
- Weber JA, Baxter DH, Zhang S, Huang DY, How Huang K, Jen Lee M, Galas DJ, Wang K: The MicroRNA Spectrum in 12 Body Fluids. Clin Chem. 2010, 56 (11): 1733-1741.PubMedView ArticleGoogle Scholar
- Chertow GM, Burdick E, Honour M, Bonventre JV, Bates DW: Acute Kidney Injury, Mortality, Length of Stay, and Costs in Hospitalized Patients. J Am Soc Nephrol. 2005, 16 (11): 3365-3370.PubMedView ArticleGoogle Scholar
- Langmead B, Trapnell C, Pop M, Salzberg S: Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009, 10 (3): R25-PubMed CentralPubMedView ArticleGoogle Scholar
- Lee LW, Zhang S, Etheridge A, Ma L, Martin D, Galas D, Wang K: Complexity of the microRNA repertoire revealed by next-generation sequencing. RNA. 2010, 16 (11): 2170-2180.PubMed CentralPubMedView ArticleGoogle Scholar
- Griffiths-Jones S, Bateman A, Marshall M, Khanna A, Eddy SR: Rfam: an RNA family database. Nucleic Acids Res. 2003, 31 (1): 439-441.PubMed CentralPubMedView ArticleGoogle Scholar
- Anders S, Huber W: Differential expression analysis for sequence count data. Genome Biol. 2010, 11 (10): R106-PubMed CentralPubMedView ArticleGoogle Scholar
- Robinson MD, McCarthy DJ, Smyth GK: edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010, 26 (1): 139-140.PubMed CentralPubMedView ArticleGoogle Scholar
- Soneson C, Delorenzi M: A comparison of methods for differential expression analysis of RNA-seq data. BMC Bioinformatics. 2013, 14 (1): 91-PubMed CentralPubMedView ArticleGoogle Scholar
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 credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.