- Research article
- Open Access
The beagle dog MicroRNA tissue atlas: identifying translatable biomarkers of organ toxicity
- Erik M. Koenig†1Email author,
- Craig Fisher†1,
- Hugues Bernard1,
- Francis S. Wolenski1,
- Joseph Gerrein1,
- Mary Carsillo1,
- Matt Gallacher1,
- Aimy Tse1,
- Rachel Peters1,
- Aaron Smith2,
- Alexa Meehan1,
- Stephen Tirrell1 and
- Patrick Kirby1
© The Author(s). 2016
Received: 6 November 2015
Accepted: 21 July 2016
Published: 17 August 2016
MicroRNAs (miRNA) are varied in length, under 25 nucleotides, single-stranded noncoding RNA that regulate post-transcriptional gene expression via translational repression or mRNA degradation. Elevated levels of miRNAs can be detected in systemic circulation after tissue injury, suggesting that miRNAs are released following cellular damage. Because of their remarkable stability, ease of detection in biofluids, and tissue specific expression patterns, miRNAs have the potential to be specific biomarkers of organ injury. The identification of miRNA biomarkers requires a systematic approach: 1) determine the miRNA tissue expression profiles within a mammalian species via next generation sequencing; 2) identify enriched and/or specific miRNA expression within organs of toxicologic interest, and 3) in vivo validation with tissue-specific toxicants. While miRNA tissue expression has been reported in rodents and humans, little data exists on miRNA tissue expression in the dog, a relevant toxicology species. The generation and evaluation of the first dog miRNA tissue atlas is described here.
Analysis of 16 tissues from five male beagle dogs identified 106 tissue enriched miRNAs, 60 of which were highly enriched in a single organ, and thus may serve as biomarkers of organ injury. A proof of concept study in dogs dosed with hepatotoxicants evaluated a qPCR panel of 15 tissue enriched miRNAs specific to liver, heart, skeletal muscle, pancreas, testes, and brain. Dogs with elevated serum levels of miR-122 and miR-885 had a correlative increase of alanine aminotransferase, and microscopic analysis confirmed liver damage. Other non-liver enriched miRNAs included in the screening panel were unaffected. Eli Lilly authors created a complimentary Sprague Dawely rat miRNA tissue atlas and demonstrated increased pancreas enriched miRNA levels in circulation, following caerulein administration in rat and dog.
The dog miRNA tissue atlas provides a resource for biomarker discovery and can be further mined with refinement of dog genome annotation. The 60 highly enriched tissue miRNAs identified within the dog miRNA tissue atlas could serve as diagnostic biomarkers and will require further validation by in vivo correlation to histopathology. Once validated, these tissue enriched miRNAs could be combined into a powerful qPCR screening panel to identify organ toxicity during early drug development.
MicroRNAs (miRNA) are varied in length, under 25 nucleotides, single-stranded noncoding RNA that regulate post-transcriptional gene expression via translational repression or mRNA degradation . miRNAs have properties which makes them ideal candidates for bio-fluid based biomarkers. miRNAs are detectable in a wide variety of biofluids [2–5], are stable in serum and can be quantified using sensitive and specific qPCR assays. miRNAs have been associated with and put forth as putative biomarkers of human disease, including hepatitis C (miR-122) , cardiovascular diseases (miR-192)  and various types of cancers .
Many miRNAs are highly conserved among mammalian species and demonstrate tissue-specific expression [1, 9–14]. In humans, there are approximately 2600 annotated human miRNAs in miRBase (v. 21) . Studies comparing expression among vertebrates (zebrafish, chickens, and mice) have demonstrated that miRNA expression is not strictly conserved . Not all human miRNAs are detected in other mammalian species [10, 14, 16]. However, there are examples of validated biomarkers such as miR-122, which is detectable in the blood after liver injury in mice, rats, dogs and humans [9, 17–20]. Thus, miRNA biomarkers have great potential as translatable tools to monitor for organ specific injury throughout all phases of drug safety assessment.
The identification of novel miRNA biomarkers requires a systematic approach to: 1) determine the miRNA tissue expression profiles within a mammalian species via next generation sequencing (miR-seq); 2) identify enriched miRNA expression within organs of toxicologic interest; and 3) validate in vivo using select toxicants. miRNA expression atlases and candidate biomarkers of organ specific injury have been investigated to varying degrees in rodents and humans [5, 12, 21]. Dogs are a relevant preclinical species that are used in nonclinical safety assessments [22, 23]. However, little is known about dog miRNA expression or annotation of sequences.
This work describes an miRNA atlas of 16 tissues from male Marshall Beagle dogs that was compiled using miR-seq. In the miRNA dog tissue atlas presented, 106 tissue enriched miRNAs were identified, 66 of which demonstrated a high level of enrichment in different tissues. Comparison of the dog atlas to currently available literature suggests that most highly enriched miRNAs are conserved across the target organs. The most unique miRNAs identified in the dog tissue atlas were found in the central nervous system, which is consistent with other species [5, 24]. Quantitative reverse transcription PCR (Q-RT-PCR) of tissue extracts confirmed the enrichment of most miRNAs identified by miR-seq analysis. Finally, the specificity of the candidate miRNA biomarkers was demonstrated through proof of concept (POC) studies with established liver toxicants.
In parallel to this effort, Eli Lilly authors have constructed a rat miRNA tissue atlas using miR-seq in 5 male and 5 female Sprague Dawley rats and 21 and 23 organs, respectively. Publication of both miRNA tissue atlases will provide an invaluable resource for the characterization of potential biomarkers of organ injury in the two most common nonclinical species used in drug development safety assessment.
The ultimate objective of this collaboration will be to assemble a serum based Q-RT-PCR screening panel of validated miRNA biomarkers of organ injury. The basic miRNA panels described herein and in the complimentary rat miRNA tissue atlas manuscript may be valuable tools in the early stages of drug development to understand test article-related target organ profiles.
Tissue sample collection
Sixteen different tissues were harvested from five 10-month old male Marshall Beagle dogs at Bioreclamation, LLC including: liver, heart, testis, lung, skeletal muscle (quadriceps), kidney, thymus, brain, sciatic nerve, pancreas, small intestine (duodenum, jejunum, and ileum), colon, bone marrow, and plasma. In order to minimize the contribution of blood cells to miRNA expression detected in individual organs, dogs were perfused with phosphate buffered saline prior to collection of organs analyzed in the dog tissue atlas. Tissue samples were placed in RNAlater® solution (Life Technologies), snap frozen in liquid nitrogen, and then stored at -80 °C until RNA isolation and miR-seq analysis.
For tissues, samples were homogenized with lysis buffer (RLT [Qiagen] + 1 %-β mercaptoethanol) using a rotor-stator homogenizer. RNA was extraced from tissue homogenate with KingFisher™ Pure RNA Tissue Kit (Thermo Scientific). Briefly, tissue homogenates were combined with magnetic beads and ethanol, and loaded onto a KingFisher Magnetic Particle Processor (Thermo Scientific). Samples were DNase-treated, washed, and eluted in RNase-free water. Due to the low level of total RNA isolated from dog sciatic nerve, the five individual samples were pooled and split into three sample replicates for analysis. RNA isolation of the bone marrow was performed with 3-fold volume of TriReagent LS (Life Technologies). Additional details are in the Additional file 1.
Assessment of RNA purity, quantity, and integrity
The purity and quantity of total RNA samples were determined by absorbance readings at 260 and 280 nm using a NanoDrop ND-1000 UV-Vis spectrophotometer (Thermo Scientific). The integrity of total RNA was assessed by capillary electrophoresis using an Agilent Bioanalyzer 2100 (Agilent).
miR-seq annotation and sample analysis
Dog miRNAs were aligned to the canine genome identified by the annotation prefix canis familiaris annotation (cfa). The annotation of mature miRNAs within miRBase was much smaller for dog (453 miRNAs) than for rat (765) or human (2588). Given the hypothesis that most miRNAs are conserved across mammalian species, a more comprehensive annotation of the dog miRNA tissue atlas was conducted using annotation from dog, rat, and human mature miRNA sequences. This was accomplished in four sequential steps: 1) identification of mature miRNA sequences for dog, rat, and human (n = 3896) using miRBase v.21; 2) consolidation of mature miRNA by eliminating sequences that were conserved across all three species (n = 3355); 3) alignment of sequences from the dog genome using Omicsoft Sequence Aligner (OSA) v4, and retained sequences that aligned with zero mismatches (n = 1087); and 4) merging of miRNAs that aligned to approximately the same location. This process resulted in the annotation of 857 miRNAs in the dog genome, including 314 putative dog miRNAs that were annotated using rat and/or human miRNA sequences which were not present in the current dog miRNA genome.
Total RNA from dog tissue samples were analyzed using miR-seq. A TruSeq Small RNA Library Kit (Illumina) was used for library construction of dog tissue RNA. Sequencing was performed on the GAIIx (Illumina) at 36 base pair read length and targeting 12 million reads per sample. Adaptor sequences were clipped and OSA v4 (http://omicsoft.com/osa) was used to align the reads to the dog genome (CanFam3.1)  allowing zero mismatches and excluding any reads that aligned to greater than 10 genomic locations. Expression levels were quantified using bedtools (https://github.com/arq5x/bedtools2). Aligned reads were counted in miRNA expression only if they overlapped with at least 70 % of known or putative miRNA loci. Additional details of analyses are described in the Additional file 1.
Hereafter, we define two terms to describe the enrichment status of miRNA: tissue enriched (TE) and highly tissue enriched (HTE). Briefly, an miRNA was considered to be TE if the fold-change level within the tissue of interest was > 5, and the FDR-corrected Wilcoxon Rank-Sum p-value was < 0.05 . To be classified as HTE, the miRNA median expression required > 5 fold higher expressed compared to the maximum expression value observed in any other tissue analyzed and was only detected in ≤ 2 tissues within the atlas. Details of the criteria used for these terms are described in the Additional file 1.
miRNA atlas verification
Q-RT-PCR analysis was conducted to verify the expression levels of candidate miRNAs identified via miR-seq. Briefly, total RNA was reverse transcribed using the Applied Biosystems TaqMan MicroRNA RT Kit (Life Technologies) according to manufacturer instructions. For each PCR reaction, cDNA was mixed with Applied Biosystems TaqMan Universal PCR Master Mix, No AmpErase® UNG (Life Technologies). All amplifications were performed in triplicate using a 7900HT Fast Real-Time PCR System (Applied Biosystems). Technical replicate threshold cycle (Ct) values were averaged for each sample. Ct values greater than 38 were removed from the analysis because they were within a 10-fold level of the final cycle and were deemed less reliable or not expressed. Relative miRNA expression quantification was derived using the ΔCt method . Additional experimental details are described in the Additional file 1.
Atlas utilization in dog toxicology studies
Naïve male beagle dogs (approximately 10 months old) were acquired from Marshall BioResources. Animals were housed in a facility accredited by the Association for Assessment and Accreditation of Laboratory Animal Care International. At the end of study, liver samples were harvested and placed in neutral buffered formalin for a minimum of 24 h before trimming, embedding, and histopathologic processing. Slides for histopathologic assessment were prepared using hematoxylin and eosin staining. Histopathologic assessment of liver sections was conducted by a board certified veterinary pathologist.
Assessment of a 5 miRNA biomarker screening panel in a 7-day repeat-dose study with Compound X
Three male dogs per group were administered either vehicle (10 % hydroxypropyl-beta-cyclodextrin [HP-β-CD] + 3.5 % NaHCO3) or Compound X at 500 mg/kg via oral gavage twice daily for seven consecutive days. Animals were fasted overnight prior to the end of study on Day 8 (24 h after the last dose). Liver tissues were harvested for histopathologic assessment and serum was collected for serum chemistry and Q-RT-PCR analysis using a 5 miRNA biomarker screening panel. The miRNAs evaluated were cfa-miR-122 and cfa-miR-885 (both considered to be highly liver enriched) as well as putative negative controls for heart (cfa-miR-1), pancreas (cfa-miR-216a), and muscle (cfa-miR-133a).
Assessment of a 22 miRNA biomarker screening panel in a 14-day repeat-dose study with Compound Y
Six male dogs per group were administered either vehicle (0.5 % methylcellulose) or 600 mg/kg of Compound Y via oral gavage daily for 14 days. Blood samples were collected on Days 1, 7, and 14 (at predose, and 24 h postdose) for serum chemistry and Q-RT-PCR analysis of 22 miRNAs. The miRNAs evaluated were cfa-miR-122 and cfa-miR-885, and an extended panel of non-liver enriched miRNAs for heart/muscle (cfa-miR-1, -133a/b, and -208b), heart (cfa-miR-499), muscle (cfa-miR-206), brain and sciatic nerve (cfa-miR-212 and -432), testis (cfa-miR-34b/c), pancreas (cfa-miR-216a/b), putative liver (cfa-miR-21 and -192) and ubiquitous controls (cfa-miR-16, -29a, and -186). Dogs were euthanized on Day 15 (24 h after the last dose) and liver tissue was harvested for microscopic assessment.
Results and discussion
Atlas tissue, RNA, and sequencing quality
This work describes the first comprehensive miRNA atlas of 16 tissues from male beagle dogs that was generated using miR-seq. Similar miRNA tissue atlases exist for other toxicologically relevant nonclinical species, such as the rat and monkey [5, 28]. Collectively, these atlases may contribute to the identification and monitorability of organ specific toxicities and provide an additional tool in the drug development process.
miRNA sequencing: tissue enrichment analysis
miRNA alignment, annotation, and quantification
Brain tissue had the greatest number of enriched miRNAs (55), including 30 HTE and 25 TE miRNAs (Fig. 3). This is consistent with previous findings that approximately 70 % of miRNAs are expressed in the brain . Rat studies identified numerous enriched brain miRNA, including miR-9, -124, -128, -184, and -219 [5, 9, 30]. There are 14 known enriched miRNAs in the human brain . Four dog brain HTE miRNAs (cfa-miR-9, -124, -128, and -219) share sequence homology with rat and human, and are cross-species biomarker candidates for brain injury. Of these, miR-9 was among the highest expressed miRNAs in the dog atlas (Additional file 3: Figure S2).
The testes had the second largest number of tissue enriched miRNAs (23), including 16 HTE and 7 TE (Fig. 3). Interestingly, 10 testes HTE miRNAs were among the most highly expressed miRNAs across all tissues, including the two highest level of expression miRNAs (cfa-miR-508b and -202) observed within the entire dog miRNA tissue atlas (Additional file 3: Figure S2). Except for miR-34b/c, which has a role germ cell maturation in murine testes [31, 32], none of these miRNAs have reported associations with testis. Testes enriched cfa-miR-34b/c, and the other testes HTE dog miRNA identified here may be invaluable as biomarkers for an organ that presently lacks a definitive diagnostic test for tissue injury.
There were 12 enriched miRNAs observed in the pancreas (7 HTE and 5 TE) (Fig. 3). Of these, only cfa-miR-216a correlated with observed expression in both mouse and rat pancreas . A recent rat study demonstrated that plasma levels of miR-216a and miR-216b increased after pancreatic injury . The pancreatic HTE cfa-miR-217-3p had an annotated homolog in rat but not human (Additional file 5: Figure S4). This novel finding suggests that human miR-217 should in fact be two separate mature miRNA sequences: miR-217-5p and miR-217-3p.
The two liver HTE miRNAs in dogs were miR-122 and the reverse complement miR-3591 (Fig. 3). In the literature, miR-122 is liver specific for mice, rats, and humans and has been used as a biomarker for liver injury [5, 10, 12]. miR-855 is a liver TE miRNA reported to be significantly elevated in sera from patients with liver disease . Others have reported miR-192 as a liver specific biomarker in multiple species [9, 20]. In contrast, the findings presented in this dog atlas show that cfa-miR-192 was expressed at relatively high levels (>650 RPM) in all tissues. Because this miRNA did not meet the TE miRNA criteria, it does not appear to be a suitable biomarker of liver injury in the dog.
A single HTE miRNA was identified in the dog heart (cfa-miR-499) and in skeletal muscle (cfa-miR-206), and both tissues had the same 5 TE miRNAs (cfa-miR-1, -133a-5p, -133a-3p, -133b, and -208) (Fig. 3). Enrichment of miR-499 in the heart has been demonstrated in rat, monkey, and human [36, 37]. Likewise, miR-206 is enriched in the skeletal muscle from rodents and humans . The overlap of heart and muscle miRNA expression is supported by the Spearman’s correlation (Fig. 1) and is consistent with the literature . The family of ‘myomiRs’ (miR-1, -133a, -133b, and -208) is known to be expressed in both cardiac and muscle tissues [39, 40], and was enriched in the myocardium of rats, dogs, and monkeys .
Analysis of dog plasma identified 3 HTE (cfa-miR-423a, -144-5p, and -16-3p) and 8 TE miRNAs (Fig. 3). Of the plasma HTE miRNAs, only miR-144 is known to be blood cell specific  and is involved in erythrocyte homeostasis . Analysis of dog bone marrow identified 1 HTE miRNA (miR-8865) and 11 TE miRNAs. Only TE miRNAs were identified in the thymus (3) and GI tract (7). There were no HTE or TE miRNAs found in either the dog lung, kidney or duodenum.
There were 3 dog miRNAs (cfa-miR-16, -29a, and -186) with ubiquitous expression across all tissues. miR-16 is regarded as a stably expressed miRNA in human and has been used as a reference miRNA for normalization of serum miRNA biomarkers for breast, prostate, and colorectal cancers [43–45]. Both miR-29a and miR-186 are some of the most frequently detected miRNAs in plasma and serum . Likewise, all three miRNAs were abundantly detected, have low CV < 10 %, and standard deviation of 1 or below when compared across all tissues analyzed in the dog miRNA atlas. While these miRNAs were used to normalize miR-seq validation by qPCR data, they could not be used in POC studies because their serum levels unexpectedly differed between control animals and dosed animals.
Tissue samples were analyzed by Q-RT-PCR to confirm the expression levels of miRNAs identified as HTE or TE by miR-seq analysis. A total of 22 miRNAs (Additional file 6: Figure S5) were selected for qPCR validation including the following 14 biomarker candidates of organ toxicity: liver (cfa-miR-122 and -885), pancreas (cfa-miR-216a/b); heart (cfa-miR-499); muscle (cfa-miR-206); heart/muscle (cfa-miR-1, -133a/b, and -208); testis (cfa-miR-34b/c); and brain and sciatic nervous tissues (cfa-miR-212, -432, and -885), and 5 miRNAs reported in the literature (cfa-miR-21, -192, -193a/b, and -200). For Q-RT-PCR normalization controls, 3 miRNAs identified in the dog atlas (cfa-miR-16, -29a, and -186) were selected on the basis that they are ubiquitously expressed at similar levels in all tissues examined.
Highest expressed miRNA in the kidney
Although enriched miRNAs were not detected in the dog kidney, this is a unique organ because urine can be evaluated for biomarkers. Thus, the highest expressed kidney miRNAs may possess the greatest potential as urinary biomarkers of kidney injury. Circulating blood miRNAs are not thought to pass through the kidney into the urine intact . High expression in dogs was observed for cfa-miR-10a, -10b, -22, -181a, -191, -192, and -378 (Additional file 4: Figure S3). There was a strong overlap between the miRNA identified here in the dog whole kidney and the highest expressed miRNA identified by Ichii et al. in the dog cortex (7/10) and medulla (8/10) (Additional file 7: Figure S6) . Kidney enriched miR-10a and miR-192 were previously identified as potential circulating biomarkers of renal injury in rats , however miR-192 had high expression in multiple dog tissues and was not identified as kidney enriched in the current study. Other miRNAs have been implicated in various kidney diseases, including in renal fibrosis (miR-22) , nephritic syndrome (miR-181a) , and renal cell carcinoma (miR-378) . Levels of miR-191 and miR-378 were elevated in the urine of rats dosed with nephrotoxicants [52, 53]. It remains to be determined whether any of these highly expressed renal miRNAs can be detected in the urine and outperform standard protein-based biomarkers of renal injury.
Liver safety biomarkers POC studies
One concern for using miRNAs as circulating biomarkers is how to normalize across samples . While studies have investigated the use of spiking samples with miRNA from exogenous sources such as Arabidopsis or C. elegans [55, 56] or normalization by serum or plasma volume, miRNA in the current study were normalized using total counts. The expression patterns of miRNA in plasma versus serum are similar [46, 56, 57]. Due to sample availability only serum was analyzed in the dog POC studies.
Assessment of a 5 miRNA biomarker screening panel after 7 days of dosing with Compound X
Assessment of a 20 miRNA biomarker screening panel after 14 days of dosing with Compound Y
Six dogs were dosed daily with 600 mg/kg of Compound Y or vehicle control for 14 days. Serum samples collected on Days 1, 7, and 14 (predose and 24 h postdose) were analyzed for ALT and AST levels. Serum was also analyzed by Q-RT-PCR for a panel of 20 tissue enriched and potential miRNA biomarkers, including those identified for liver (cfa-miR-122 and -885), heart/muscle (cfa-miR-1, -133, and -206), testis (miR-34b/c), pancreas (cfa-miR-216), brain (cfa-miR-212), and ubiquitously expressed cfa-miR-193b.
Elevation of miR-122 in the plasma/serum after administration of hepatotoxic compounds in rats [58, 59] or in various human disease conditions [19, 60] is well established. The results from these two POC studies represent the first published data demonstrating the utility of miR-122 and miR-885 as potential biomarkers of liver injury in the dog. Levels of miR-122 increased in correlation with minimal liver injury in the absence of ALT elevations, and thus may have superior sensitivity than standard liver function tests. This is consistent with a previous rat study that detected an increase in miR-122 levels prior to ALT elevation . While the sensitivity of miR-122 versus ALT has not been definitively determined in any species, these collective studies indicate that miR-122 is a translatable diagnostic to detect liver injury .
miRNAs hold promise as circulating biomarkers of organ specific injury. While recent studies focused on miRNA tissue expression in rodent species and humans, there is a significant lack of miRNA data available for dogs. The beagle dog is an important nonclinical species for assessment of human drug safety. Identification of miRNA biomarkers of tissue specific injury in the dog may improve prediction and monitorability of potentially adverse effects of compounds intended for human use. This miRNA tissue atlas may also serve as a reference for the identification of novel biomarkers of organ damage in nonclinical species which are translatable to patients in the clinic. Analysis of the dog miRNA tissue atlas identified 66 HTE miRNAs which offer a starting point for exploratory pre-clinical safety monitoring, further proof of concept and biomarker validation studies beyond the liver biomarker potential described herein.
The dog miRNA annotation has been expanded through the identification of homologs to rat and human miRNAs in the dog genome which demonstrate a high level of sequence homology and similar levels of tissue expression among rats and humans. There are likely to be additional miRNAs in dogs that might have been tissue enriched but were not annotated. As the quality of annotation for dog miRNA improves, so will the resolution and quality of the dog miRNA database. The expression data may be re-evaluated at a later date to see if previously unannotated miRNA sequences can be identified. Conversely, information gathered from the dog miRNA atlas can be used to fill in gaps in the human annotation, as in the case with miR-217. The data generated for the dog miRNA tissue atlas is publically available for further analysis (The accession number is: GSE83278).
In collaboration with Eli Lilly, co-analyses of the dog and rat miRNA sequencing data determined which tissue enriched miRNAs were conserved between the two species. The next phase of the atlas will be to further demonstrate in vivo the correlation of tissue enriched miRNA to tissue specific injury. Subsequent animal studies are necessary to determine the sensitivity and specificity of the miRNA qPCR panel. Candidate biomarkers may eventually be used as an early readout of tissue injury in drug development.
ΔCt, change in threshold cycle; ALT, alanine aminotransferase; AST, aspartate aminotransferase; cDNA, complementary deoxyribonucleic acid; cfa, canis familiaris annotation; Ct, threshold cycle; CV, coefficient of variation; FDR, false discovery rate; GI, gastrointestinal; HP-β-CD, hydroxypropyl-beta-cyclodextrin; HTE, highly tissue enriched; miRNA was considered HTE if the miRNA median expression required > 5 fold higher expressed compared to the maximum expression value observed in any other tissue analyzed and was only detected in ≤ 2 tissues within the atlas; miRNA, microRNA; mRNA, messenger ribonucleic acid; OSA, Omicsoft Sequence Aligner; PCR, polymerase chain reaction; POC, proof of concept; qPCR, see Q-RT-PCR; Q-RT-PCR, quantitative reverse transcription polymerase chain reaction; RIN, ribonucleic acid integrity number; RNA, ribonucleic acid; RPM, reads per million; TE, tissue enriched; miRNA was considered to be TE if the fold-change levels within the tissue of interest was > 5, and the FDR-corrected Wilcoxon Rank-Sum p-value was < 0.05
Bioreclamation; Takeda Molecular Pathology and Comparative Medicine Groups.
John Cally, Mark Farmen and Han Wu (Eli Lilly non-author contributors).
Experimental concept and design, sample selection: EK, CF, PK, MC. Tissue collection, trimming, and pathology: MG, AT, MC, RP, ST. Molecular assays supporting Atlas and POC verification studies (RT-qPCR, miR-seq): EK, HB, ST. Sequencing alignment and miRNA enrichment analysis: EK, JG. POC study and biomarker interpretation: EK, CF, HB, AS, PK. Manuscript writing and figures: EK, CF, FW, HB, AM, PK. Manuscript review: All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
This study complies with all applicable sections of the final rules of the Animal Welfare Act (Code of Federal Regulations, Title 9), The Public Health Service Policy on Humane Care and Use of Laboratory Animal Welfare, and The Guide for the Care and Use of Laboratory Animals from the National Research Council.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.
- Ambros V. The functions of animal microRNAs. Nature. 2004;431(7006):350–5.View ArticlePubMedGoogle Scholar
- Mitchell PS, Parkin RK, Kroh EM, Fritz BR, Wyman SK, Pogosova-Agadjanyan EL, Peterson A, Noteboom J, O’Briant KC, Allen A, et al. Circulating microRNAs as stable blood-based markers for cancer detection. Proc Natl Acad Sci U S A. 2008;105(30):10513–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Gilad S, Meiri E, Yogev Y, Benjamin S, Lebanony D, Yerushalmi N, et al. Serum MicroRNAs Are Promising Novel Biomarkers. PLoS One. 2008;3(9):e3148. doi:10.1371/journal.pone.0003148.View ArticlePubMedPubMed CentralGoogle Scholar
- Etheridge A, Lee I, Hood L, Galas D, Wang K. Extracellular microRNA: a new source of biomarkers. Mutat Res. 2011;717(1-2):85–90.View ArticlePubMedPubMed CentralGoogle Scholar
- Minami K, Uehara T, Morikawa Y, Omura K, Kanki M, Horinouchi A, Ono A, Yamada H, Ohno Y, Urushidani T. miRNA expression atlas in male rat. Scientific data. 2014;1:140005.View ArticlePubMedPubMed CentralGoogle Scholar
- Koberle V, Waidmann O, Kronenberger B, Andrei A, Susser S, Fuller C, Perner D, Zeuzem S, Sarrazin C, Piiper A. Serum microRNA-122 kinetics in patients with chronic hepatitis C virus infection during antiviral therapy. J Viral Hepat. 2013;20(8):530–5.View ArticlePubMedGoogle Scholar
- Matsumoto S, Sakata Y, Suna S, Nakatani D, Usami M, Hara M, Kitamura T, Hamasaki T, Nanto S, Kawahara Y, et al. Circulating p53-responsive microRNAs are predictive indicators of heart failure after acute myocardial infarction. Circ Res. 2013;113(3):322–6.View ArticlePubMedGoogle Scholar
- Schwarzenbach H, Nishida N, Calin GA, Pantel K. Clinical relevance of circulating cell-free microRNAs in cancer. Nat Rev Clin Oncol. 2014;11(3):145–56.View ArticlePubMedGoogle Scholar
- Laterza OF, Lim L, Garrett-Engele PW, Vlasakova K, Muniappa N, Tanaka WK, Johnson JM, Sina JF, Fare TL, Sistare FD, et al. Plasma MicroRNAs as sensitive and specific biomarkers of tissue injury. Clin Chem. 2009;55(11):1977–83.View ArticlePubMedGoogle Scholar
- Bentwich I, Avniel A, Karov Y, Aharonov R, Gilad S, Barad O, Barzilai A, Einat P, Einav U, Meiri E, et al. Identification of hundreds of conserved and nonconserved human microRNAs. Nat Genet. 2005;37(7):766–70.View ArticlePubMedGoogle Scholar
- Ason B, Darnell DK, Wittbrodt B, Berezikov E, Kloosterman WP, Wittbrodt J, Antin PB, Plasterk RH. Differences in vertebrate microRNA expression. Proc Natl Acad Sci U S A. 2006;103(39):14385–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Landgraf P, Rusu M, Sheridan R, Sewer A, Iovino N, Aravin A, Pfeffer S, Rice A, Kamphorst AO, Landthaler M, et al. A mammalian microRNA expression atlas based on small RNA library sequencing. Cell. 2007;129(7):1401–14.View ArticlePubMedPubMed CentralGoogle Scholar
- Guo X, Su B, Zhou Z, Sha J. Rapid evolution of mammalian X-linked testis microRNAs. BMC Genomics. 2009;10:97.View ArticlePubMedPubMed CentralGoogle Scholar
- Iwama H, Kato K, Imachi H, Murao K, Masaki T. Human microRNAs originated from two periods at accelerated rates in mammalian evolution. Mol Biol Evol. 2013;30(3):613–26.View ArticlePubMedGoogle Scholar
- Kozomara A, Griffiths-Jones S. miRBase: annotating high confidence microRNAs using deep sequencing data. Nucleic Acids Res. 2014;42(Database issue):D68–73.View ArticlePubMedGoogle Scholar
- Chen DB, Wang W. Human placental microRNAs and preeclampsia. Biol Reprod. 2013;88(5):130.View ArticlePubMedPubMed CentralGoogle Scholar
- Harrill AH, Eaddy JS, Rose K, Cullen JM, Ramanathan L, Wanaski S, Collins S, Ho Y, Watkins PB, Lecluyse EL. Liver biomarker and in vitro assessment confirm the hepatic origin of aminotransferase elevations lacking histopathological correlate in beagle dogs treated with GABAA receptor antagonist NP260. Toxicol Appl Pharmacol. 2014;277(2):131–7.View ArticlePubMedGoogle Scholar
- Wang K, Zhang S, Marzolf B, Troisch P, Brightman A, Hu Z, Hood LE, Galas DJ. Circulating microRNAs, potential biomarkers for drug-induced liver injury. Proc Natl Acad Sci U S A. 2009;106(11):4402–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhang Y, Jia Y, Zheng R, Guo Y, Wang Y, Guo H, Fei M, Sun S. Plasma microRNA-122 as a biomarker for viral-, alcohol-, and chemical-related hepatic diseases. Clin Chem. 2010;56(12):1830–8.View ArticlePubMedGoogle Scholar
- Starkey Lewis PJ, Dear J, Platt V, Simpson KJ, Craig DG, Antoine DJ, French NS, Dhaun N, Webb DJ, Costello EM, et al. Circulating microRNAs as potential markers of human drug-induced liver injury. Hepatology (Baltimore, Md). 2011;54(5):1767–76.View ArticleGoogle Scholar
- Liang Y, Ridzon D, Wong L, Chen C. Characterization of microRNA expression profiles in normal human tissues. BMC Genomics. 2007;8:166.View ArticlePubMedPubMed CentralGoogle Scholar
- Guth BD. Preclinical cardiovascular risk assessment in modern drug development. Toxicol Sci. 2007;97(1):4–20.View ArticlePubMedGoogle Scholar
- Schultze AE, Walker DB, Turk JR, Tarrant JM, Brooks MB, Pettit SD. Current practices in preclinical drug development: gaps in hemostasis testing to assess risk of thromboembolic injury. Toxicol Pathol. 2013;41(3):445–53.View ArticlePubMedGoogle Scholar
- Shao NY, Hu HY, Yan Z, Xu Y, Hu H, Menzel C, Li N, Chen W, Khaitovich P. Comprehensive survey of human brain microRNA by deep sequencing. BMC Genomics. 2010;11:409.View ArticlePubMedPubMed CentralGoogle Scholar
- Hoeppner MP, Lundquist A, Pirun M, Meadows JR, Zamani N, Johnson J, Sundstrom G, Cook A, FitzGerald MG, Swofford R, et al. An improved canine genome and a comprehensive catalogue of coding genes and non-coding transcripts. PLoS One. 2014;9(3):e91172.View ArticlePubMedPubMed CentralGoogle Scholar
- Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J. R. Stat. Soc. Ser. B Methodol. 1995;1:289–300. http://www.stat.purdue.edu/~doerge/BIOINFORM.D/FALL06/Benjamini%20and%20Y%20FDR.pdf.
- Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2001;25(4):402–8.View ArticlePubMedGoogle Scholar
- Peng X, Thierry-Mieg J, Thierry-Mieg D, Nishida A, Pipes L, Bozinoski M, Thomas MJ, Kelly S, Weiss JM, Raveendran M, et al. Tissue-specific transcriptome sequencing analysis expands the non-human primate reference transcriptome resource (NHPRTR). Nucleic Acids Res. 2015;43(Database issue):D737–42.View ArticlePubMedGoogle Scholar
- Adlakha YK, Saini N. Brain microRNAs and insights into biological functions and therapeutic potential of brain enriched miRNA-128. Mol Cancer. 2014;13:33.View ArticlePubMedPubMed CentralGoogle Scholar
- Ogata K, Sumida K, Miyata K, Kushida M, Kuwamura M, Yamate J. Circulating miR-9* and miR-384-5p as potential indicators for trimethyltin-induced neurotoxicity. Toxicol Pathol. 2015;43(2):198–208.View ArticlePubMedGoogle Scholar
- Bouhallier F, Allioli N, Lavial F, Chalmel F, Perrard MH, Durand P, Samarut J, Pain B, Rouault JP. Role of miR-34c microRNA in the late steps of spermatogenesis. RNA. 2010;16(4):720–31.View ArticlePubMedPubMed CentralGoogle Scholar
- Bao J, Li D, Wang L, Wu J, Hu Y, Wang Z, Chen Y, Cao X, Jiang C, Yan W, et al. MicroRNA-449 and microRNA-34b/c function redundantly in murine testes by targeting E2F transcription factor-retinoblastoma protein (E2F-pRb) pathway. J Biol Chem. 2012;287(26):21686–98.View ArticlePubMedPubMed CentralGoogle Scholar
- Goodwin D, Rosenzweig B, Zhang J, Xu L, Stewart S, Thompson K, Rouse R. Evaluation of miR-216a and miR-217 as potential biomarkers of acute pancreatic injury in rats and mice. Biomarkers. 2014;19(6):517–29.View ArticlePubMedGoogle Scholar
- Endo K, Weng H, Kito N, Fukushima Y, Iwai N. MiR-216a and miR-216b as markers for acute phased pancreatic injury. Biomed Res. 2013;34(4):179–88.View ArticlePubMedGoogle Scholar
- Gui J, Tian Y, Wen X, Zhang W, Zhang P, Gao J, Run W, Tian L, Jia X, Gao Y. Serum microRNA characterization identifies miR-885-5p as a potential marker for detecting liver pathologies. Clin Sci (Lond). 2011;120(5):183–93.View ArticleGoogle Scholar
- Xu J, Zhao J, Evan G, Xiao C, Cheng Y, Xiao J. Circulating microRNAs: novel biomarkers for cardiovascular diseases. J Mol Med (Berl). 2012;90(8):865–75.View ArticleGoogle Scholar
- Qin H, Chen GX, Liang MY, Rong J, Yao JP, Liu H, Wu ZK. The altered expression profile of microRNAs in cardiopulmonary bypass canine models and the effects of mir-499 on myocardial ischemic reperfusion injury. J Transl Med. 2013;11(1):154.View ArticlePubMedPubMed CentralGoogle Scholar
- McCarthy JJ. MicroRNA-206: the skeletal muscle-specific myomiR. Biochim Biophys Acta. 2008;1779(11):682–91.View ArticlePubMedPubMed CentralGoogle Scholar
- Ma G, Wang Y, Li Y, Cui L, Zhao Y, Zhao B, Li K. MiR-206, a key modulator of skeletal muscle development and disease. Int J Biol Sci. 2015;11(3):345–52.View ArticlePubMedPubMed CentralGoogle Scholar
- Malizia AP, Wang DZ. MicroRNAs in cardiomyocyte development. Wiley Interdiscip Rev Syst Biol Med. 2011;3(2):183–90.View ArticlePubMedPubMed CentralGoogle Scholar
- Vacchi-Suzzi C, Hahne F, Scheubel P, Marcellin M, Dubost V, Westphal M, Boeglen C, Buchmann-Moller S, Cheung MS, Cordier A, et al. Heart structure-specific transcriptomic atlas reveals conserved microRNA-mRNA interactions. PLoS One. 2013;8(1):e52442.View ArticlePubMedPubMed CentralGoogle Scholar
- Rasmussen KD, Simmini S, Abreu-Goodger C, Bartonicek N, Di Giacomo M, Bilbao-Cortes D, Horos R, Von Lindern M, Enright AJ, O’Carroll D. The miR-144/451 locus is required for erythroid homeostasis. J Exp Med. 2010;207(7):1351–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Davoren PA, McNeill RE, Lowery AJ, Kerin MJ, Miller N. Identification of suitable endogenous control genes for microRNA gene expression analysis in human breast cancer. BMC Mol Biol. 2008;9:76.View ArticlePubMedPubMed CentralGoogle Scholar
- Chang KH, Mestdagh P, Vandesompele J, Kerin MJ, Miller N. MicroRNA expression profiling to identify and validate reference genes for relative quantification in colorectal cancer. BMC Cancer. 2010;10:173.View ArticlePubMedPubMed CentralGoogle Scholar
- Schaefer A, Jung M, Miller K, Lein M, Kristiansen G, Erbersdobler A, Jung K. Suitable reference genes for relative quantification of miRNA expression in prostate cancer. Exp Mol Med. 2010;42(11):749–58.View ArticlePubMedPubMed CentralGoogle Scholar
- Blondal T, Jensby Nielsen S, Baker A, Andreasen D, Mouritzen P, Wrang Teilum M, Dahlsveen IK. Assessing sample and miRNA profile quality in serum and plasma or other biofluids. Methods. 2013;59(1):S1–6.View ArticlePubMedGoogle Scholar
- Ichii O, Otsuka S, Ohta H, Yabuki A, Horino T, Kon Y. MicroRNA expression profiling of cat and dog kidneys. Res Vet Sci. 2014;96(2):299–303.View ArticlePubMedGoogle Scholar
- Wang JF, Zha YF, Li HW, Wang F, Bian Q, Lai XL, Yu G. Screening plasma miRNAs as biomarkers for renal ischemia-reperfusion injury in rats. Med Sci Monit. 2014;20:283–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Long J, Badal SS, Wang Y, Chang BH, Rodriguez A, Danesh FR. MicroRNA-22 is a master regulator of bone morphogenetic protein-7/6 homeostasis in the kidney. J Biol Chem. 2013;288(51):36202–14.View ArticlePubMedPubMed CentralGoogle Scholar
- Sui W, Lin H, Li H, Yan Q, Chen J, Dai Y. Circulating microRNAs as potential biomarkers for nephrotic syndrome. Iran J Kidney Dis. 2014;8(5):371–6.PubMedGoogle Scholar
- Redova M, Poprach A, Nekvindova J, Iliev R, Radova L, Lakomy R, Svoboda M, Vyzula R, Slaby O. Circulating miR-378 and miR-451 in serum are potential biomarkers for renal cell carcinoma. J Transl Med. 2012;10:55.View ArticlePubMedPubMed CentralGoogle Scholar
- Nassirpour R, Mathur S, Gosink MM, Li Y, Shoieb AM, Wood J, O’Neil SP, Homer BL, Whiteley LO. Identification of tubular injury microRNA biomarkers in urine: comparison of next-generation sequencing and qPCR-based profiling platforms. BMC Genomics. 2014;15:485.View ArticlePubMedPubMed CentralGoogle Scholar
- Kanki M, Moriguchi A, Sasaki D, Mitori H, Yamada A, Unami A, Miyamae Y. Identification of urinary miRNA biomarkers for detecting cisplatin-induced proximal tubular injury in rats. Toxicology. 2014;324:158–68.View ArticlePubMedGoogle Scholar
- Meyer SU, Pfaffl MW, Ulbrich SE. Normalization strategies for microRNA profiling experiments: a ‘normal’ way to a hidden layer of complexity? Biotechnol Lett. 2010;32(12):1777–88.View ArticlePubMedGoogle Scholar
- Shi R, Chiang VL. Facile means for quantifying microRNA expression by real-time PCR. Biotechniques. 2005;39(4):519–25.View ArticlePubMedGoogle Scholar
- Moldovan L, Batte KE, Trgovcich J, Wisler J, Marsh CB, Piper M. Methodological challenges in utilizing miRNAs as circulating biomarkers. J Cell Mol Med. 2014;18(3):371–90.View ArticlePubMedPubMed CentralGoogle Scholar
- Wang K, Yuan Y, Cho JH, McClarty S, Baxter D, Galas DJ. Comparing the MicroRNA spectrum between serum and plasma. PLoS One. 2012;7(7):e41561.View ArticlePubMedPubMed CentralGoogle Scholar
- Starckx S, Batheja A, Verheyen GR, Jonghe SD, Steemans K, Dijck BV, Singer M, Bogdan N, Snoeys J, Vinken P, et al. Evaluation of miR-122 and other biomarkers in distinct acute liver injury in rats. Toxicol Pathol. 2013;41(5):795–804.View ArticlePubMedGoogle Scholar
- Lardizabal MN, Rodriguez RE, Nocito AL, Daniele SM, Palatnik JF, Veggi LM. Alteration of the microRNA-122 regulatory network in rat models of hepatotoxicity. Environ Toxicol Pharmacol. 2014;37(1):354–64.View ArticlePubMedGoogle Scholar
- Yang X, Salminen WF, Shi Q, Greenhaw J, Gill PS, Bhattacharyya S, Beger RD, Mendrick DL, Mattes WB, James LP. Potential of extracellular microRNAs as biomarkers of acetaminophen toxicity in children. Toxicol Appl Pharmacol. 2015;284(2):180–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Antoine DJ, Dear JW, Lewis PS, Platt V, Coyle J, Masson M, Thanacoody RH, Gray AJ, Webb DJ, Moggs JG, et al. Mechanistic biomarkers provide early and sensitive detection of acetaminophen-induced acute liver injury at first presentation to hospital. Hepatology. 2013;58(2):777–87.View ArticlePubMedPubMed CentralGoogle Scholar