- Research article
- Open Access
The beagle dog MicroRNA tissue atlas: identifying translatable biomarkers of organ toxicity
BMC Genomicsvolume 17, Article number: 649 (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.
Data from five samples per tissue type were included in miR-seq analysis to generate the dog miRNA tissue atlas with the following exceptions: only 4 of 5 colon samples and 3 sciatic nerve miRNA samples were included in downstream enrichment analysis. The 3 sciatic nerve samples comprised 2 individual samples and 3 samples were pooled due to low RNA yields. Quality assessments used for sample analysis included Spearman’s (Rank r 2) correlation of miRNA expression across all tissue samples (Fig. 1a), total tissue RNA (RNA integrity number [RIN]) the number of miRNAs, and the percent mapped reads in each tissue sample (Fig. 1b). Based on Spearman’s correlation, the majority of tissue expression profiles demonstrated good reproducibility between tissue samples of the same organ (average r 2 > 0.85) with the exception of muscle, duodenum, testis, and thymus (Fig. 1a). Tissues with similar miRNA expression patterns between organs included heart and muscle (average r 2 = 0.8) and tissues of the upper and lower gastrointestinal (GI) tract (duodenum, jejunum, ileum, and colon; average r 2 = 0.9). In contrast, the brain demonstrated the most distinct miRNA expression profile (average r 2 = 0.3) when compared to all other tissues examined.
miRNA sequencing: tissue enrichment analysis
miRNA alignment, annotation, and quantification
There are 453 annotated mature miRNA sequences for the dog in miRBase (v.21), far less than for rat (765) or human (2588). To increase the numbers of known miRNAs, a de novo annotation of dog miRNAs was performed using known rat and human sequences as guides (Fig. 2a). Sequence alignment with the combined rat, dog, and human miRNA annotation resulted in the identification of 857 putative homologous miRNAs within the dog genome. Of these, 650 miRNAs (355 from dog and 295 from rat and/or human) were detected by miR-seq and aligned to the genome of dogs, and 106 potential miRNA biomarkers were identified (Fig. 2b). A comprehensive list of the 106 potential miRNA biomarkers and the tissues in which they were enriched is presented in Fig. 2c. The similarity in sequences between the dog miR-seq data and the 295 miRNAs annotated using rat/human miRNA homologs strongly suggests that these miRNAs, while not currently annotated in the dog genome, are likely well conserved in the dog as well as in rat and human. Additionally, 20 of the 66 highly tissue expressed miRNAs were not annotated in the dog and instead were identified through rat/human homologs.
This study identified 298 miRNAs with normalized reads per million (RPM) greater than 100 in at least one tissue. Of these, 214 miRNAs were annotated dog miRNAs and 84 were identified using rat and/or human sequences. A total of 106 dog miRNA sequences were TE, 66 miRNAs were HTE and of these, 60 were enriched in only one tissue while the remaining 6 were enriched in two tissues (Fig. 3). Additional information about specific mature miRNA homologs in rats and humans is presented in Additional file 2: Figure S1. The top 20 TE/HTE miRNAs from the dog atlas (Additional file 3: Figure S2) and top 10 expressed miRNAs per tissue (Additional file 4: Figure S3) are also presented.
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.
Eleven of 15 miRNAs examined by miR-seq and Q-RT-PCR miRNA platforms demonstrated comparable expression levels of individual miRNAs selected for verification with r 2 values > 0.6 (Fig. 4). The exceptions (cfa-miR-432, -499, and -212) could not be differentiated between tissues analyzed by Q-RT-PCR. The 3 highly expressed ubiquitous miRNAs demonstrated good precision, as measured by the coefficient of variance (%CV = standard deviation/mean), for all tissues in both detection platforms; miRNA-SEQ (cfa-miR-16 [4.9 %], -29a [4.7 %], and -186 [1.8 %]) and Q-RT-PCR (cfa-miR-16 [4.1 %], -29a [4.6 %], and -186 [2.7 %]).
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
Three dogs were dosed twice daily for seven days with 500 mg/kg of Compound X or vehicle. At 24 h post-Day 7 dosing, two dogs demonstrated mild to moderate liver necrosis (Fig. 5a) and had an approximately 5-fold elevation in serum alanine aminotransferase (ALT) and aspartate aminotransferase (AST) which correlated with mild to moderate liver necrosis (Fig. 5b). Analysis of serum miRNA from these two dogs showed the highest elevation (ΔCt) values in cfa-miR-122 and cfa-miR-855 compared to control animals (Fig. 5c). Histopathological analysis did not identify microscopic changes in heart, muscle and pancreas in dogs treated with Compound X, and the levels of non-liver enriched miRNAs (cfa-miR-1, -133 and -216) did not significantly change when compared to control dogs (data not shown). Taken together, these results highlight the specificity of miR-122 and miR-885 for compound-induced liver injury.
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.
Elevations in liver enriched miR-122 and miR-885 correlated with increases in ALT and AST (Fig. 6a and b) and mild to moderate hepatocellular necrosis was observed in 2 of 6 animals (Dogs 1 and 3 on Day 14 [data not shown]). Dogs 1 and 3 had high correlation between elevations in serum miR-122 and ALT (r 2 = 0.79) and AST (r 2 = 0.90) levels while elevations in serum miR-885 were higher ALT (r 2 = 0.87) and AST (r 2 = 0.97). Additionally, Dog 3 demonstrated elevated miR-122 on Day 7 at the 24 h time point while ALT remained at baseline levels. The four other animals dosed with Compound Y did not demonstrate test article-related changes in serum chemistry, miRNA or liver histopathology. Elevation of serum cfa-miR-193b in the two dogs with liver injury (Dogs 1 and 3) was likely due to the ubiquitous expression of the miRNA in the dog, including the liver. The transient serum elevations of heart/muscle TE miRNAs (cfa-miR-1 and -133) and muscle HTE miRNA (cfa-miR-206) were not correlated with microscopic findings (data not shown) and may be due to injury during animal handling. The remaining non-liver miRNAs were not detected in any serum samples.
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
Ambros V. The functions of animal microRNAs. Nature. 2004;431(7006):350–5.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Guo X, Su B, Zhou Z, Sha J. Rapid evolution of mammalian X-linked testis microRNAs. BMC Genomics. 2009;10:97.
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.
Kozomara A, Griffiths-Jones S. miRBase: annotating high confidence microRNAs using deep sequencing data. Nucleic Acids Res. 2014;42(Database issue):D68–73.
Chen DB, Wang W. Human placental microRNAs and preeclampsia. Biol Reprod. 2013;88(5):130.
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.
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.
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.
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.
Liang Y, Ridzon D, Wong L, Chen C. Characterization of microRNA expression profiles in normal human tissues. BMC Genomics. 2007;8:166.
Guth BD. Preclinical cardiovascular risk assessment in modern drug development. Toxicol Sci. 2007;97(1):4–20.
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.
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.
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.
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.
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.
Adlakha YK, Saini N. Brain microRNAs and insights into biological functions and therapeutic potential of brain enriched miRNA-128. Mol Cancer. 2014;13:33.
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.
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.
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.
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.
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.
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.
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.
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.
McCarthy JJ. MicroRNA-206: the skeletal muscle-specific myomiR. Biochim Biophys Acta. 2008;1779(11):682–91.
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.
Malizia AP, Wang DZ. MicroRNAs in cardiomyocyte development. Wiley Interdiscip Rev Syst Biol Med. 2011;3(2):183–90.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Shi R, Chiang VL. Facile means for quantifying microRNA expression by real-time PCR. Biotechniques. 2005;39(4):519–25.
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.
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.
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.
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.
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.
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.
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.
Supplementary Methods. (DOCX 18 kb)
List of 106 enriched miRNAs by tissue. miRNAs are annotated as tissue enriched (TE) or highly tissue enriched (HTE). The family of miRNA and mature miRNAs in dog, human and rat are given. (PDF 55 kb)
Top 20 enriched miRNAs. 10 of top 20 are testis. (PDF 38 kb)
Top Ten Expressed (log10 RPM values) miRNA per tissue type. The list is sorted by dog miRNA tissue atlas enrichment and prevalence of top 10 expressed miRNAs across the atlas tissues. (PDF 45 kb)
Benefit and need for enhanced annotations for miRNA studies. (A) Pancreas enriched miR-217-3p (rno-miR-217-3p AUCAGUUCCUAAUGCAUUGCCU) identified in the dog miRNA tissue atlas was found as a novel un-annotated human miRNA. It is conserved in rats, dogs, and humans, but is only annotated in miRBase as rat. (B) There is evidence of un-annotated human reads aligning and expressed in human expression profile experiments on the miRBase website. This study suggests human miRBase annotation should reflect 2 mature microRNA sequences, miR-3p and -5p. (PDF 64 kb)
Q-RT-PCR assay list. (PDF 39 kb)
Top 20 expressed kidney miRNA. Top 20 expressed kidney miRNA found in the dog miRNA tissue atlas compared to previously published miRNA expression data from macro-dissection studies of the cortex and medulla of the cat and dog (Ichii et al). Shaded boxes indicate that the dog miRNA was among the top 10 expressed kidney miRNAs in the dog atlas; bolded miRNAs indicate that the miRNA was among the top 20. Comparison between dog whole kidney (dog altas) and data from dog cortex and medulla (Ichii et al) show good correlation with 7/10 miRNAs and 8/10 miRNAs expressed in the top 10 miRNAs, respectively, when compared to kidney expression observed in the dog atlas. Similarly, a comparison of dog whole kidney data (dog atlas) to cat cortex and medulla (Ichii et al) show good correlation as well with 7/10 and 6/10 miRNAs expressed in the top 10 miRNAs, respectively, when compared to dog whole kidney data (dog atlas). (PDF 41 kb)