BACKGROUND
the drug discovery process is a long and expensive endeavor with many experimental molecules failing in late stage in vivo toxicology studies.
the high attrition rate of drugs in late stage clinical trials lead to the call for investigation of the utilization of species translatable biomarkers to select safer compounds early in the drug development process  <cit> .
in order to better understand toxicities in preclinical and clinical settings and to speed the development of new therapies, body fluid based biomarkers of organ damage have been sought to both identify and monitor toxicities.
serum or plasma based biomarkers capable of monitoring tissues affected by drug toxicity could allow for the detection of tissue injury from simple blood draws, thus permitting larger numbers of compounds to be tested more rapidly in studies that are less expensive and use fewer animals compared to traditional means of toxicity evaluation.
the properties of an ideal biomarker include tissue specificity, a strong correlation to organ dysfunction or toxicity, a strong association to histopathologic changes, robust sensitivity, predictability, accessibility through a noninvasive method, readily-available rapid and inexpensive assays for quantitative measurement and direct translatability across preclinical species and human.
serum based biomarkers of organ dysfunction have been useful in preclinical and clinical settings.
serum ast  and alt  were discovered in  <dig>  <cit>  and have been used in combination with total bilirubin, gamma-glutamyl transferase , alkaline phosphatase , albumin and bile acids as markers of liver injury to de-risk compounds in preclinical drug development  <cit>  as well as being utilized for monitoring liver toxicity in clinical drug development  <cit> .
additional markers of hepatotoxicity are currently under investigation  <cit> .
cardiac troponin i is a sensitive and specific serum based biomarker of cardiac injury in humans and mice and has been employed to mitigate cardiac toxicity in mouse toxicology studies with novel kinase inhibitors  <cit> .
urine based biomarkers for evaluation of renal toxicity include bun , creatinine, urinary albumin, kim- <dig> , ngal  and cystatin c  <cit> .
amylase and lipase have been used as markers for pancreatitis, but their diagnostic capability is limited by sensitivity and specificity  <cit> .
these preclinical and clinical biomarkers are carefully reviewed in dieterle,  <dig>  <cit> .
despite the success with protein based biomarkers to date, gaps in our ability to detect toxicities in a number of organs prevent the widespread utilization of biomarkers in toxicology studies.
current gaps in the minimally-invasive biomarker portfolio include gastrointestinal injury, cns  injury, reproductive toxicity and specific markers of glomerular and liver injury.
clinical chemistry parameters offer important information regarding the health of a given organ but can be limited by lack of specificity and sensitivity as in the cases of alt and ast for liver injury and amylase and lipase for injury in the exocrine pancreas  <cit> .
while protein based biomarkers provide useful information, they can be limited in their translation across species and to the clinical setting due to species differences in the biomarkers and the need for species specific immunoassay reagents.
traditional pre-clinical methods for assessing the potential toxicity of a compound include clinical chemistry, histopathology, hematology, physiological measurements and clinical observations, but the concordance rate to human toxicity was reported by olsen et al.
to be 71 % leaving room for significant improvement in preclinical toxicology studies  <cit> .
unfortunately, the concordance of clinical chemistry or other biomarkers was not directly compared to establish concordance between preclinical and clinical studies; several groups advised more use of biomarkers in preclinical and clinical studies  <cit> .
in a subsequent study, adverse drug reactions in each target organ from compounds approved in japan from 2001– <dig> were found to be 48 % concordant between preclinical toxicology studies and humans while changes in laboratory parameters displayed a 30-70 % correlation  <cit> .
it is evident that additional endpoints indicative of organ toxicities that are translatable and/or predictive from preclinical species to human are necessary to improve the safety of compounds in preclinical drug development and clinical trials.
with the end goal of identifying tissue specific, species translatable blood based biomarkers of tissue injury, we set out to characterize the rat microrna expression profile in various organs and create a rat microrna body atlas to identify specific and enriched mirnas in individual tissues.
micrornas are 21– <dig> nucleotide long rna molecules whose cellular function appears to be to inhibit translation or cause degradation of their mrna target or targets  <cit> .
some mirnas are expressed in high abundance in a tissue specific manner, are resistant to degradation in the serum, are well conserved in preclinical species including mouse, rat, macaque and human, and are detectable in a wide variety of biologic fluids .
in addition, sensitive, specific and readily available qpcr based assays exist to measure mirnas making them attractive candidates for serum based biomarkers for organ injury.
perhaps the most well characterized mirna for use as a toxicity biomarker is the liver enriched mir-122-5p.
this mirna was increased in the serum of rats that were given the liver toxicant trichlorobromomethane, but not increased in the serum of rats administered the skeletal muscle toxicants  <dig> , <dig>  tetramethyl-p-phenylenediamine  and an hmg-coa reductase inhibitor  <cit> .
ast and alt were increased when animals were given liver or skeletal muscle toxicants, demonstrating that mir- <dig> was indicative of liver injury only and not muscle injury, thus displaying an advantage in specificity over alt and ast.
mice treated with acetaminophen displayed increases of mir- <dig> in the serum 2 h before an increase of ast/alt  <cit> .
additionally, mir- <dig> and liver enriched mir- <dig> were increased in the serum of patients who had acetaminophen induced hepatotoxicity demonstrating the potential for these mirnas to be used as both preclinical and clinical biomarkers for liver injury  <cit> .
in a further study mir- <dig>  hmgb <dig>  cytokeratin  <dig> and gldh were able to predict patients who would develop liver failure as a result of acetaminophen overdose even when alt remained in the normal range  <cit> .
while several mirna tissue profiling experiments have been reported in mouse, rat and human, these rodent experiments do not encompass the range of tissues relevant to toxicology studies for drug candidates.
in addition, published studies relied upon microarray and/or qpcr based technologies which could suffer from a lack of specificity and sensitivity.
microarrays also do not allow for the identification of novel mirnas due to their dependence on known sequences and may not account for tissue expression of isomirs, which are mirnas that have variations from the reference mirna sequence  <cit> .
qpcr also lacks specificity in cases where  <dig> mirnas have a closely aligned sequence.
as an example, mir- <dig> expressed in skeletal muscle and heart and mir- <dig> expressed specifically in skeletal muscle differ by only  <dig> bases.
several of the current qpcr based assays for mir- <dig> also detect synthetic mir- <dig>  indicating that caution must be taken when using this approach.
next generation sequencing  is a preferred method for mirna profiling because it makes direct, unbiased and quantitative measures of known and unknown mirnas.
in order to realize the true potential of mirnas as biomarkers of organ damage it is critical to comprehensively understand the tissue expression pattern of mirs and isomirs, and this can be best achieved using ngs.
therefore we have constructed a rat mirna body atlas using illumina mirna sequencing from  <dig> male and  <dig> female sprague dawley rats and  <dig> and  <dig> organs of toxicologic interest from each sex, respectively .
in parallel to this effort, several of the authors have constructed a similar canine body atlas from  <dig> tissues from  <dig> male marshall beagle dogs.
publication of both mirna body atlases will provide an invaluable resource for the characterization of potential biomarkers of organ injury in  <dig> most commonly used species for preclinical safety assessment.
here we present the findings of the rat mirna body atlas in which we have identified  <dig>  previously unreported mirnas in the rat and detailed numerous tissue specific and enriched mirnas.
we have discovered several novel tissue specific mirnas in each species, some of which are conserved across both rat and dog, and some which are species specific .
comparison of the rat and dog atlases yielded candidate liver and pancreas specific biomarkers.
the potential pancreas biomarkers were tested in rat and dog toxicology studies using the pancreatic toxicant caerulein, while the potential liver biomarkers were tested in toxicology studies conducted by takeda.
in these proof-of-concept studies, mirnas performed as well as or better than traditional markers of pancreatic injury, demonstrated a much larger dynamic range and increased tissue specificity over traditional biomarkers  <dig>
methods
necropsy
five male and five female sprague dawley rats 12– <dig> weeks in age were obtained from charles rivers laboratories and euthanized by co <dig> asphyxiation.
the rats were necropsied and tissues were placed into rnalater® within 5 min of asphyxiation.
organs  were processed in batches according to approximate tissue size as follows and placed into rnalater® at 4 °c for 24 h and were transferred to-20c the following day.
in batch # <dig>  greater than 150 mg of liver, stomach, ileum, jejunum and duodenum were placed into rnalater® for subsequent total rna isolation.
the following portions of the intestine corresponding to the segment that is routinely collected for histology were collected: the stomach with approximately 5 cm of the duodenum attached, a 10 cm section from the midpoint of the jejunum and the distal 5 cm of ileum.
a syringe containing rnalater® was used to flush the contents from the lumen of the collected segments prior to placing them into the rnalater® collection tube.
in batch # <dig>  ~150 mg of kidney, pancreas, brain, testis, biceps and soleus were placed into rnalater® for subsequent total rna isolation.
in batch # <dig>  30-50 mg of adrenal, heart, ovary and uterus were placed into rnalater® for subsequent total rna isolation.
rna extraction
tissues were removed from rnalater®, frozen in liquid nitrogen and pulverized according to the covaris tissue pulverization protocol  <cit> .
tissues were placed in an appropriate volume of qiazol and homogenized using lysing d matrix tubes .
qiazol volume was scaled up according to the mass of tissue to be homogenized.
total rna was next isolated and eluted into 80ul of 95 °c molecular biology grade water.
rna quality was assessed using the agilent  <dig> bioanalyzer using the rna pico kit.
all rnas except for pancreas, 3/ <dig> hippocampus samples and  <dig> uterus sample had rin  values of  <dig> or higher with an average rin value of  <dig> .
for samples that displayed rin values lower than  <dig> the rna extraction protocol was repeated using additional qiazol/tissue homogenate.
most samples achieved rin values greater than  <dig> as shown in fig.
<dig>  whole blood was placed into rnalater® for blood and all rnas were isolated using the mirneasy kit from qiagen following the manufacturer’s instructions.
dorsal root ganglion rnas were provided by takeda pharmaceuticals.
rnas were quantitated using a nanodrop and, if necessary, were concentrated by speed vac to obtain an appropriate concentration for library preparation prior to sequencing.fig.
<dig> the tissues used for construction of the rat mirna body atlas and boxplots of the rin values associated with the tissues from  <dig> male and  <dig> female sprague dawley rats are illustrated above.
rnas were extracted from rat tissues and were tested for rna quality using the agilent bioanalyzer.
most rnas were of very high quality as judged by rin values.
pancreas rnas had the lowest rin values and also had the fewest sequences that mapped to mature mirnas
mirna seq library construction
library construction was performed according to illumina’s truseq small rna sample prep protocol .
library quality was assessed using the agilent bioanalyzer dna chip and quantified using ribo green.
libraries were subjected to miseq analysis to ensure library quality.
body atlas libraries were assigned to gels and lanes on the flow cells in a manner designed to mitigate batch effect.
subsequent hiseq analysis was performed and 50 bp single end reads were generated on the illumina hiseq  <dig> platform resulting in 4– <dig> million reads per sample.
data analysis by institution
read mapping and quantification
data was analyzed by eli lilly, maastricht university and niehs using separate analysis methods.
the output generated from each institution was then collated to create a venn diagram of the tissue enriched mirnas identified by each analysis method.
read mapping and quantification for the rat and dog mirna body atlases was performed in the following manner by eli lilly: fastq files were processed to remove adaptor sequences and discard reads < 17 bp in length with cutadapt v.  <dig> .
<dig>  <cit>  with options –a tggaattctcgggtgccaagg–quality–base 33–q 20–match-read-wildcards –m  <dig>  all trimmed reads that contained an ‘n’ were discarded.
identical sequences from the same sample were combined into a single sequence in the form expected by mirdeep <dig>  <cit> .
that is, with ids of the form “nnn_#_x#”, where nnn is a unique alphabetic code identifying the sample, the first ‘#’ is a unique integer, and ‘x#’ records the number of times this sequence occurred in the input for the sample.
quantifier.pl from the mirdeep <dig> package v.  <dig> .
<dig> was used to generate mir alignment files  against known mirs from mirbase  <dig> with options–p–p < org > _hairpin.fa –m < org > _mature.fa–r trimmed_reads.fa–d.
this was done separately for rat, mouse, human, and c. elegans known mirs.
the.mrd files were then parsed with a custom perl script as follows.
each isomir sequence in an alignment was associated with the corresponding mature mir identifier.
if it aligned to a mirna precursor, but not with an expected mature mir sequence, it was identified as < mir > −pre to indicate this.
usually these correspond with reverse strand mirs that have not yet been annotated as alternate mature forms.
sometimes these appear to be microrna offset rnas  <cit> , and sometimes they appear to represent incompletely processed sequence.
if a given sequence was identified as aligning with more than one precursor, it was associated with all potential names as a composite name.
that is, a sequence that aligned to both let-7a-5p and let-7f-5p was assigned to the composite mature mir let-7a-5p;let-7f-5p to indicate for subsequent analyses that this identification was ambiguous.
all sequences that had not yet been identified against known rat mirs were then looked for in the analysis with respect to known mouse mirs, and so on against human and finally c. elegans mirs.
this process allows for identification of conserved rodent mirs that have not yet been annotated in rat, and against conserved mammalian mirs that have not yet been annotated in rat and mouse.
in parallel to the above identification of known mirs, we used mirdeep2’s novel mir identification process.
additional details are described in additional file 1: supplementary methods.
tissue specific and enriched mirna identification
data analysis of the rat and dog sequencing data was performed by eli lilly as follows.
tissue specific mirnas are defined as mirnas that are expressed at a raw count of greater than  <dig> in  <dig> out of  <dig> rats in  <dig> organ/ tissue only.
for example, the brain may be the only organ with mir-x expression, but the hypothalamus may be the only tissue that has mir-x expressed.
therefore, mir-x would be specific to the hypothalamus tissue.
tissue enriched mirnas are defined as mirnas that are expressed at a raw count of greater than  <dig> in  <dig> out of  <dig> rats in a few organs/tissues, typically 3– <dig> tissues.
the definition of tissue enriched was purposely vague because a mirna may be enriched in several organs, but could still be a useful biomarker if toxicity is only present in  <dig> of the organs the mirna is enriched in.
sex specific mirnas are identified as mirnas that are expressed in  <dig> tissue/organ at a raw count of greater than  <dig> in 4/ <dig> males or females and in 0/ <dig> rats in the opposite sex.
raw and normalized data were analyzed for tissue specific and tissue enriched mirnas.
animal level data is obtained after preprocessing of the raw sequencing data.
for each microrna, there are measurements from  <dig> animals  and potentially  <dig> different tissues.
these  <dig> tissues were from  <dig> organs.
in order to find tissue specific mirnas, tissue level counts were obtained for each microrna by summing over all the animals.
the tissue level counts were aggregated to organ level by selecting the maximum of the tissue counts to represent the organ count.
for example, brain  has  <dig> different tissue types , the maximum of these  <dig> tissue level counts will be used to represent the expression levels of brain.
after these data processing steps, a count data matrix is generated, where each row corresponds to a mir and each column corresponds to an organ type.
non-negative matrix factorization   <cit>  is a commonly used unsupervised learning algorithm to decompose a nonnegative matrix x into  <dig> non-negative matrices w and h. in our case, x is the observed count matrix of dimension nx <dig>  where n is the number of mirs.
each column of w defines a mir group and each column of h represents the expression pattern of the mir group corresponding to the organ type.
based on this decomposition, we can find a representative set of mirs that corresponds to the specific high expression pattern for the organ.
to identify tissue specific mirs, it has to be organ specific .
thus we start from each organ specific mir, and examine if it is also highly expressed in a particular type of tissue or a few tissues, which belong to the current organ.
for example, a mirna may be specific to brain , but it may only be specific to the hippocampus .
to determine the expression level of each mir, we fitted a two-component mixture of poisson distribution to the animal level tissue counts data.
the component with larger intensity of the poisson distribution is naturally used to model high counts and the component with lower intensity is used to model background noise.
the larger component of the poisson mixture model naturally corresponds to high expression level, whereas the smaller component of the poisson mixture model naturally corresponds to the low expression level.
for an organ specific microrna to be tissue specific, it has to be expressed at a high level in  <dig> tissue of this organ.
after the mixture model is fitted and the intensities are estimated, we can determine which tissues the microrna is expressed in according to their estimated intensities and then declare tissue specific micrornas afterwards.
additionally, the raw sequencing data was manually analyzed in each tissue for each mirna and statistical analysis of the manually identified tissue specific and enriched mirnas was then conducted to verify or dispute the tissue specificity or enrichment.
read mapping and quantitation by maastricht university
read mapping and quantification was performed in the following manner by maastricht university: mirdeep <dig>  <cit>  was used to predict putative novel mirnas.
the raw fastq files from the  <dig> rat samples were mapped to the rat genome .
mirdeep <dig> output was then parsed to keep only the  <dig> predicted mirna precursors with a score of  <dig> or above.
raw reads were trimmed off the 3′ adapter sequence using the java algorithm “adaptor remover” from the freely available uea srna workbench .
the parameters were set to look for a perfect match with the first  <dig> bases of the 3′ adapter to be considered a hit.
the trimmed sequences with a size below  <dig> or above 35 bp were discarded.
during this trimming step, every identical sequence was pulled and each iteration was summed.
the trimmed reads were then mapped to the list of rat precursor mirna obtained either from mirbase, or generated de novo from the mirdeep <dig> prediction using a fast short read mapping software patman  <cit> , configured to allow no mismatches or gap.
patman output was then parsed to 1) assign a unique name to each unique sequence consisting of the pre-mirna name, the pre-mirna arms  and the mapping coordinates .
2) divide the total read count of each unique sequence by the number of assigned loci in the complete mapping.
3) create a mature output table by pulling the different isomir count for each mirna species .
the parsed patman output were normalized with the trimmed mean of m-values method  and filtered to remove all mirnas where the normalized read counts were fewer than  <dig> reads in all the  <dig> samples.
total read counts were then computed for each individual mirna.
tissue enriched mirna identification by maastricht university
a mirna was called tissue specific or tissue enriched when the percentage of reads for a single tissue  was respectively higher than  <dig> or 50 % of the total read count.
read mapping and quantitation by niehs
in order to identify isomirs the raw count number of each given isomir was first converted to the percentage of expression compared to the mature count observed for its corresponding mirna.
by comparing this percentage, we were able to report the mirna for which the most expressed isomirs differs between tissues.
analysis of the data performed at niehs was executed in the following manner: raw fastq files with the mirna-seq reads were checked for quality using a tqe standard based on trimming the adapters, quality filtering at q ≤  <dig> and the elimination of short reads <  <dig> bases long.
in some cases when minor adapter sequences  were observed after initial process, additional rounds of processing were applied until all the criteria were met.
in the end, reads longer than  <dig> base pairs were also eliminated, resulting in quality filtered reads between 14–24 bp.
reads passing the tqe filtering were then aligned to the rat mirbase database v. <dig> using bwa version  <dig> .1-r104q with “-n 2” and default parameters otherwise.
only perfect matches were further used in the quantification, which was performed by summarizing read count level measurements at each mature mirbase mirs.
datapoints with zeros were imputed with  <dig> .
total counts for each mir were normalized by transcripts per million reads .
prior to statistical analysis, a transformation from floats to integers was performed by ceiling the data.
tissue enriched mirna identification by niehs
to detect the tissue-enriched micrornas, we adopted the one-vs-rest strategy when applying statistic models.
“differentially” expressed mirs  in  <dig> tissue against the combination of all other tissues were identified using a quasi-poisson approach  with shrunken dispersion estimates  <cit> .
during the modeling process, the offset was computed on the 3rd quartile.
significant demirs were detected at a nominal p-value <  <dig> .
rat and dog toxicity studies
caerulein studies: studies were run at covance laboratories in greenfield, in and histopathology was performed by a board certified pathologist.
due to the short half-life of caerulein  <cit> ,  <dig> male sprague dawley rats per group  were treated with  <dig> intraperitoneal doses of vehicle or caerulein , 1 h apart, at  <dig>   <dig> and 50 μg/kg.
blood was collected at  <dig>   <dig>   <dig>   <dig> and 48 h after dosing and examined for clinical chemistry, hematology and microrna changes with terminal serum taken at  <dig>   <dig> and 48 h. urine was collected by cystocentesis at  <dig>   <dig>   <dig>   <dig> and 48 h and examined for total protein and specific gravity.
toxicokinetic analysis was performed at  <dig> ,  <dig>   <dig>   <dig>   <dig> and 48 h post dose.
histopathology of adrenal glands, bone , brain stem, cerebellum, cerebrum, duodenum, heart, ileum, jejunum, kidneys, liver, lung, muscle , pancreas, spleen, stomach and thymus was examined at  <dig>   <dig> and 48 h post dose.
for canine studies,  <dig> beagle dogs per group,  <dig> male and  <dig> female obtained from covance research products, were given  <dig> doses  <dig> h apart of caerulein at  <dig>   <dig> or  <dig> μg/kg.
blood was collected at pre-dose,  <dig>   <dig>   <dig>   <dig>   <dig> and  <dig> h after the third dose and examined for clinical chemistry, hematology and mirna changes.
urine was collected by cystocentesis and examined for changes in specific gravity, bun, creatinine and other urinary parameters at necropsy at  <dig> and 72 h. histopathology on organs including adrenal glands, bone , brain stem, cerebellum, cerebrum, colon, epididymis, eyes, heart, ileum, jejunum, kidneys, liver, lung, muscle , sciatic nerve, ovaries, pancreas, prostate, spleen, stomach, testis, thymus, thyroid, uterus and vagina was performed at  <dig> and 72 h post dose.
animals were euthanized by barbiturate overdose or anesthesia and exsanguination at necropsy.
qpcr of serum samples from caerulein studies
for each sample, 100 μl of serum from all studies was subjected to total rna isolation, including mirnas, using the mirneasy tm kit from qiagen .
upon the addition of qiazol 5 μl of 5 μm synthetic cel-mir-55-3p  was spiked into each sample in the rat caerulein and dog caerulein studies.
rnas were eluted in 80ul of 95 °c h2o and a water extraction control with spike-in was also included.
reverse transcription was executed with the mircury lna microrna universal cdna synthesis kit .
10 μl of the extracted rnas were used in 20 μl rt reactions.
a no template control  and two -rt controls were also included.
real-time pcr was conducted using exiqon  primers and exilent sybr green master mix .
samples were diluted 1: <dig> prior to qpcr and 5 μl of the diluted samples were used in 10 μl reactions.
data was collected on a 7900ht real-time pcr instrument .
statistical analysis of rat and dog caerulein qpcr data
for each study, the fold change data were log <dig> transformed and fitted to a two-way anova model with repeated measures.
because each animal has multiple measurements, we allow these measurements to be auto-correlated.
once we obtained the parameter estimates, for the dog study, all treatment means and 95 % confidence intervals at each time point are compared against the pre-dose group, which have fc fixed at  <dig>  for the rat study, low dose and high dose groups are compared to untreated group at each time point.
qpcr validation of tissue specific and enriched mirna
tissue enriched and tissue specific mirnas were chosen for qpcr validation.
the rnas from  <dig> rats and  <dig> tissues were normalized by mass  and subjected to cdna synthesis using exiqon universal cdna synthesis kits  and exilent sybr green mastermix .
data was collected on a 7900ht or viia <dig> real-time pcr instrument .
assays were designed to the most highly expressed mirna or isomir sequence present that mapped to the mature mirna.
for example, many isomirs mapped to mir-xyz.
for qpcr analysis we designed primers to the most highly expressed mirna sequence that mapped to mir-xyz as determined by deep sequencing.
ct values of each mirna were compared to one another to determine tissue specificity or enrichment.
a ct value was considered acceptable if the ct in a particular tissue was  <dig> cts lower than the signal in the no rt or no template  controls unless there was no signal in the no rt or ntc.
if an assay did not display a signal that was not  <dig> cts away from the no rt or ntc it was not included in the analysis.
additionally, if the melt curves for an assay displayed more than one peak the assay was not included in the analysis.
qpcr analysis of pancreas enriched mirnas
total rna from all male and female rats used in the study were normalized by mass using 250 ng of rna per sample and subjected to cdna synthesis using exiqon universal cdna synthesis kits  and exilent sybr green mastermix .
data was collected on a 7900ht or viia <dig> real-time pcr instrument .
RESULTS
our first objective was to identify mirnas expressed in rat tissues.
members of the hesi consortium were queried on which rat organs are most commonly injured during toxicology studies and therefore, would represent the most important organs to investigate for potential biomarkers of injury.
this resulted in the list of organs detailed in fig.
<dig>  dorsal root ganglion was added as an organ of investigative interest.
the brain lacks serum and cerebrospinal fluid  based biomarkers of injury, therefore the mirna content of the cerebrum, cerebellum, hippocampus and brainstem were interrogated.
additionally, whole blood was included in the study in order to determine what the contribution of blood mirnas may be to highly vascularized tissues.
total rna was isolated from the tissues of male and female sd rats and the quality was assessed using rin values  to be sure that small rnas would be sequenced as opposed to degraded rnas.
the rat microrna sequencing data, as analyzed by eli lilly, generated  <dig>  sequences that mapped to previously-identified mature micrornas in rat, mouse and human.
this resulted in an increase of  <dig>  mirnas in the rat as compared to the  <dig> mature rat mirnas in mirbase v <dig>  while eli lilly’s data analysis serves as the main source of data for discussion in this article we also collaborated with hesi members to analyze the data in an attempt to maximize our ability to discover and confirm the tissue enrichment profile of mirnas .fig.
<dig> workflow diagram.
the small rna content of rat organs was interrogated by deep sequencing and the data was given to  <dig> institutions for independent data analysis.
reads were mapped to rat, mouse and human mirnas from mirbase v <dig> and v <dig> using mirdeep <dig> and tissue specific and tissue enriched mirnas were identified
our second objective was to discover novel and confirm previously identified tissue specific and/or tissue enriched mirnas and/or isomirs that could be used as serum based biomarkers of organ injury in the rat.
a variety of data analysis software packages and statistical methods exist to identify and predict mirnas and perform differential expression on mirna-seq data  <cit> .
since these methods have different strengths and weaknesses and the most suitable data analysis method is dependent on the output required, different methods were employed by scientists at eli lilly, niehs and maastricht university to analyze the sequencing data in order to maximize discovery of tissue enriched mirnas.
the tissue enriched mirnas identified by each institution were included in a venn diagram in order to determine which tissue enriched mirnas were found in common between each data analysis method .
tissue specific mirnas identified by lilly are reported in table  <dig> and were not included in the venn diagram.
maastricht university identified  <dig> tissue enriched mirnas, eli lilly identified  <dig> tissue enriched mirnas and niehs identified  <dig> tissue enriched mirnas.
maastricht university and eli lilly had  <dig> mirnas in common, maastricht university and niehs had  <dig> mirnas in common, eli lilly and niehs had  <dig> mirnas in common and all  <dig> institutions had  <dig> tissue enriched mirnas in common .
differences in the results from each institution are due to the processing steps prior to alignment, the reference sequences the sequencing data were aligned to and the methods of identifying tissue specific and/or enriched mirnas.
despite which approach was taken each analysis identified many of the same tissue enriched mirnas and the data has been made publically available for researchers in industry, academia or government to explore and potentially discover additional tissue specific and/or enriched mirnas which may serve as novel biomarkers for safety, pharmacodynamic or other effects.fig.
<dig> venn diagram of data analysis methods.
analysis of the mirna sequencing data was independently performed by eli lilly, maastricht university and the niehs using different data analysis methods.
eli lilly identified  <dig> tissue enriched mirnas, maastricht university identified  <dig> tissue enriched mirnas and the niehs identified  <dig> tissue enriched mirnas.
the number of mirnas in common between maastricht university and eli lilly was  <dig>  the number of mirnas in common between maastricht university and the niehs was  <dig> and the number of mirnas in common between eli lilly and the niehs was  <dig>  the number of mirnas in common between all  <dig> analysis was 165
tissue specific mirnas were defined as mirnas there are expressed at a raw count of > <dig> in 8/ <dig> rats in  <dig> tissue only.
sequencing data was analyzed by non-negative matrix factorization which identified tissue specific mirnas.
organs not listed in the table did not have tissue specific mirnas
tissue specific and enriched mirnas identified by lilly are shown in  and  and the tissue enriched mirnas identified by niehs and maastricht university are shown in .
the tissue specificity/enrichment of many other tissue enriched mirnas including muscle and heart enriched mirs-133a-3p, 499-5p, mir-1a-3p as well as pancreas enriched mirs-216a-5p, 217-5p and mir-375-3p were confirmed  .
the brain, testis and liver have the most tissue specific and enriched mirnas and interestingly, tissue enriched mirnas in the brain are often expressed in the intestine.
the kidney has the fewest tissue specific and enriched mirnas of any organ examined, while skeletal muscle and heart share several of the same tissue enriched mirnas.
while previously identified tissue specific and enriched mirnas were verified, previously unidentified tissue specific/enriched mirnas were discovered including, but not limited to, pancreas enriched mirs-148a-3p, 141-3p and 200c-3p .
in order to verify that mirnas were, in fact, tissue specific or enriched we conducted qpcr verification of some of the tissue specific and/or enriched mirnas .
as shown in fig.
<dig>   <dig> of the  <dig> mirnas tested correlated very well to the sequencing data with a correlation of −  <dig>  to −  <dig> ,  <dig> mirnas displayed a correlation of −  <dig>  to −  <dig>  and  <dig> mirnas displayed poorer correlations from −  <dig>  to −  <dig>  .
mirs-196c-5p and 206-5p were not amplified by qpcr and, therefore, did not display tissue specificity or enrichment by qpcr.fig.
<dig> heat map of tissue specific and tissue enriched mirnas.
examples of tissue specific and enriched mirnas identified by this and previously published studies are displayed with tissues on the x-axis and mirnas listed on the y-axis.
mirna expression from least to highest expression is indicated as a shift from blue to red with blue representing the least expression and dark red indicating maximal expression
fig.
<dig> correlation of qpcr analysis of tissue specific and/or enriched mirnas.
tissue specific and/or enriched mirnas were chosen according to tissues that lack or require additional biomarkers of organ status and were measured in  <dig> tissues from  <dig> female rats by qpcr.
the qpcr data was compared to the sequencing raw counts for each mirna.
a perfect correlation is equal to −  <dig> because the ct value decreases as sequencing abundance increases.
twenty seven of the thirty nine mirnas tested correlated very well to the sequencing data with a correlation of −  <dig>  to −  <dig> ,  <dig> mirnas displayed a correlation of −  <dig>  to −  <dig>  and  <dig> mirnas displayed poorer correlations from −  <dig>  to −  <dig> .
mirs-196c-5p and 206-5p were not amplified by qpcr and, therefore, did not display tissue specificity or enrichment by qpcr
while each institution was able to identify tissue specific and enriched mirnas it is also imperative to understand whether the isomirs of the tissue specific and enriched mirnas are expressed in the same tissues as the parent mirna as isomir expression in other tissues could lead to inaccurate interpretation of serum mirna profiles with respect to toxicity.
eli lilly investigated the isomir expression of pancreas tissue specific and enriched mirnas and the analysis revealed that isomirs generally mirror their parent mirna expression , but some isomirs are more tissue specific than others as shown for mir- <dig> in the intestines  and mir-217-5p in the pancreas .
for example, the mature mir-217-5p sequence is 5′tactgcatcaggaactgactgg-3′ and the isomir-217-5p sequence 5′-tactgcatcaggaactgactggac-3′ are both highly enriched in the pancreas, but the mir-217-5p isomir lacks expression in the stomach, intestine and ovary while other isomirs of mir-217-5p are expressed in the stomach, intestine and ovary .
isomirs in other tissues such as isomirs of mir- <dig> in the intestine , mir-192-5p in the liver  and mir-3473f-pre in the testis  display similar isomir expression.
it appears that interpretation of serum mirna toxicity biomarker data should not be compromised by expression of isomirs from tissues lacking parent mirna expression.
rat and dog pancreas toxicity studies
conserved pancreas tissue enriched mirnas were examined in the rat and dog models of caerulein pancreatic toxicity in order to begin to characterize mirnas as species translatable serum based biomarkers of pancreatic injury .
rats were treated with vehicle,  <dig> and 50 μg/kg of caerulein and the serum was examined for clinical chemistry and mirna changes at  <dig>   <dig>   <dig>   <dig> and 48 h with histopathologic assessment of organs at  <dig>   <dig> and 48 h. minimal and slight pancreatic necrosis was noted in  <dig> vehicle treated rats at  <dig> and 48 h which had no clinical chemistry or mirna correlates and were considered to be consistent with common background while minimal, slight, moderate and marked pancreatic necrosis was noted in the  <dig> and 50 μg/kg treated rats that was interpreted to be dose related .
clinical chemistry parameters including amylase and lipase , mir-122-5p , mir-133a-3p , 148a-3p , 208a-3p  and pancreas mirnas conserved between rat and dog  were examined for changes in the serum.
mir-216a-5p, amylase and lipase were increased in the serum concurrently at 1 h, remained elevated until 8 h while mir-216a-5p remained elevated until 24 h. mir-217-5p displayed similar kinetics to mir-216a-5p except mir-217-5p generally had a larger dynamic range and remained elevated until 48 h. both mirs-216a-5p and 217-5p displayed much larger dynamic ranges than amylase or lipase and remained elevated longer.
mir-216b-5p displayed serum increases that were kinetically similar to amylase and lipase, but had a dynamic range that was less than lipase and greater than amylase.
mir-375-3p was increased from 4–24 h in the  <dig> and 50 μg/kg groups and returned to approximately vehicle level by 48 h. the pancreas enriched mirnas conserved between rat and dog  that were tested displayed increases in the serum similar to amylase , whereas mirs 320-3p,  <dig> and  <dig> were not increased .table  <dig> pancreas enriched mirnas conserved between rat and dog
tacagtactgtgataactgatc
tissue enriched mirnas conserved between rat and dog were identified and are shown above.
these mirnas were analyzed for changes in the serum in rat and dog caerulein toxicology studies.
bases underlined indicate bases that are not conserved between the rat and dog isomirs with the highest sequencing counts
fig.
<dig> histopathologic assessment of the pancreases of rats treated with  <dig> and 50ug/kg of caerulein.
histopathologic assessment of organs was performed at  <dig>   <dig> and 48 h with dose, time and animal number indicated on the x-axis and the grade of injury indicated on the y-axis .
rat pancreases displayed minimal and slight pancreatic necrosis in  <dig> vehicle treated rats at  <dig> and 48 h while dose responsive increases in the severity of necrosis was noted in the  <dig> and 50ug/kg treatred rats at  <dig>   <dig> and 48 h
fig.
<dig> clinical chemistry and serum mirna analysis of caerulein treated rats.
rats treated with  <dig>   <dig> and 50ug/kg of caerulein were examined for clinical chemistry and mirna changes .
time is indicated on the x-axis and the log  <dig> of the fold change, the mean as well as 95 % confidence interval are indicated on the y-axis.
mirnas were incereased in the serum of rats with  <dig> and  <dig> ug/kg of caerulein coincident with increases in amylase and lipase with maximal increases occurring at 8–24 h for mirs-216a-5p, 217-5p and 375-3p.
values are normalized to vehicle and to the cel-mir-55-3p spike in
marshall beagle dogs were treated with vehicle,  <dig>   <dig> and 45 μg/kg of caerulein and the serum was examined for clinical chemistry and mirna changes at  <dig>   <dig>   <dig>   <dig>   <dig> and 72 h with histopathologic assessment of organs at  <dig> and 72 h. pancreatic necrosis was noted in  <dig> dog at 24 h in the 15 μg/kg group and  <dig> and 72 h in the 45 μg/kg group .
amylase and lipase were not increased greater than two fold in the vehicle or 3 μg/kg treated dogs at any time point.
amylase and lipase were increased in the 15 μg/kg treated dogs with the largest increases at 1 h and decreasing by 8–24 h while both were increased maximally at 4-8 h in the 45 μg/kg treated group.
mir-216b-5p was increased in  <dig> vehicle treated animal and mir-217-5p was increased in  <dig> vehicle treated animals with no histopathologic or clinical chemistry correlates .
in 3 μg/kg treated dogs mir-216b-5p and mir-141-3p were increased at some time points, but mir-216b-5p was not increased beyond the level of the vehicle treated dog and did not display consistent time dependent increases.
mir-141-3p and mir-216a-5p were increased beyond vehicle treated dogs and did not display consistent time dependent increases.
mir-375-3p did display increases as compared to vehicle reaching its maximum increase at 1 h and remaining elevated until 24 h. in 15 μg/kg treated dogs mirs-216a-5p, 216b-5p, mir-375-3p and 141-3p were increased in a time dependent manner with maximal increases at 4–8 h .
mir-216a-5p displayed the greatest dynamic range reaching as high as  <dig> fold above vehicle.
in the 45 μg/kg treated dogs amylase, lipase and mirnas displayed dose and time dependent increases.
mir-216a-5p was increased maximally at 4–8 h and remained elevated until 72 h while levels as high as  <dig>  fold above vehicle were observed in individual animals.
mir-216b-5p reached levels as high as  <dig> fold above vehicle in individual animals and remained elevated until 24 h. the remaining mirnas did increase and remained elevated generally as long, but not to the degree of mirs-216a-5p and 216b-5p.
interestingly, liver enriched mir-122-5p was increased in a dose and time dependent manner while ast and alt were increased in the 45 μg/kg treated dogs while no histopatholgic correlate was observed .
the rat caerulein study did not generate similar mir- <dig> increases.
as demonstrated by these studies, mirnas can serve as species translatable serum based biomarkers of pancreatic injury.fig.
<dig> histopathologic assessment of the pancreases of dogs treated with caerulein.
marshall beagle dogs were treated with  <dig>   <dig> and  <dig> ug/kg of caerulein and histopathologic assessment of organs was performed at  <dig> and 72 h with dose, time and animal number indicated on the x-axis and the grade of injury indicated on the y-axis .
pancreatic degeneration/necrosis was not noted in dogs treated with 3ug/kg of caerulein.
pancreatic degeneration/necrosis was noted in  <dig> dog at 24 h in the 15ug/kg group and displayed dose dependent increases in severity at  <dig> and 72 h in the 45ug/kg group
fig.
<dig> clinical chemistry and serum mirna analysis of caerulein treated dogs.
dogs treated with  <dig>   <dig> and  <dig> ug/kg of caerulein were examined for clinical chemistry and serum mirna changes .
time is indicated on the x-axis and the log  <dig> of the fold change, the mean as well as 95 % confidence interval are indicated on the y-axis.. mirnas were increased in the serum of dogs trreated with 3ug/kg of caerulein prior to increases in amylase and lipase.
mirnas were dose dependently increased in the serum of dogs treated with  <dig> and  <dig> ug/kg of caerulein and correlated to amylase and lipase increases while mir-216a-5p remained elevated until 24 h. values are normalized to vehicle and to the cel-mir-55-3p spike in
discussion
we have generated a rat mirna body atlas and identified an additional ~ <dig> potential mirnas in the rat as compared to the current number of rat mirnas in mirbase v <dig>  we have also confirmed many previously identified tissue specific and enriched mirnas and have discovered many previously unknown rat tissue specific and enriched mirnas.
rat and dog mirna atlases were compared in order to determine which mirnas are conserved and may represent species translatable biomarkers of organ injury.
subsequently, the conserved pancreas enriched mirnas in rats and dogs were demonstrated in pancreas toxicology studies to be sensitive blood-based biomarkers of caerulein induced pancreatic injury in both species.
the aims of this study were  <dig> fold.
first we sought to comprehensively catalogue the mirna and isomir content of the major organs of toxicologic interest in the rat in order to discover potential mirna biomarkers of organ injury.
secondly, we wanted to identify tissue specific and tissue enriched mirnas that may serve as biomarkers of organ injury.
thirdly, we wanted to test tissue specific and/or enriched mirnas that are conserved in the rat and dog in toxicology studies to determine whether conserved mirnas could be used as serum based biomarkers of organ injury in both species.
while others have generated similar data using microarray or qpcr, next generation sequencing is advantageous in that it allows for the discovery of mirna sequences without the need for prior genome annotation  <cit> .
analysis of the sequencing data refined the tissue specific/enrichment of previously reported tissue specific mirnas and revealed many previously undiscovered tissue specific and enriched mirnas.
identification of tissue specific and enriched mirs and isomirs
independent data analysis, utilizing different techniques, was performed by members of hesi.
approximately  <dig> tissue enriched mirnas were found to be in common between the data analysis methods and additional tissue enriched mirnas were identified by each institution.
in order to ensure that interpretation of serum mirna profiles is reflective of organ specific toxicity we sought to determine the tissue specificity and enrichment of the isomirs of tissue specific and enriched mirnas identified by eli lilly.
eli lilly identified many sequences that mapped to mature mirnas which consisted of sequences that differed from the mature mirna sequence.
these sequences, termed isomirs  <cit>  contained additional bases on the 5′ and 3′ ends as compared to the reference mirna sequences and were also evaluated for tissue specificity and enrichment.
we determined that the isomirs of tissue specific and enriched mirnas, with some exceptions, did generally mirror the expression of the reference mirna sequence indicating that serum mirna profiles should accurately reflect organ toxicity.
some isomirs did display increased tissue specificity/enrichment as compared to the parent mirna and these could serve as more specific markers of organ injury when faced with multiple organ toxicity.
eli lilly’s definition of a tissue enriched mirna was purposefully vague because some mirnas may be expressed in multiple organs at a higher level than the remaining organs, but may still be reasonable candidate biomarkers of organ injury depending on the observed toxicity.
it is possible that a mirna that is enriched in several tissues could be a useful marker of tissue injury if toxicity is restricted to a single organ that the mirna is expressed in.
in the face of multiple organ toxicity this approach may not be reasonable unless additional biomarkers could improve interpretation of the injury.
consequently, some mirnas may not appear to be particularly tissue enriched as shown in additional file 5: figure s <dig>
rat and dog pancreas toxicity studies
in order to characterize pancreas enriched mirnas conserved between rat and dog as potential biomarkers of pancreatic injury we conducted rat and dog caerulein studies in which animals were treated with caerulein and concurrent clinical chemistry and histopathology were assessed along with mirna changes in time course studies.
amylase and lipase increases were noted from 1–8 h in rats in both  <dig> and 50 μg/kg dose groups while pancreatic necrosis was noted at  <dig>   <dig> and 48 h. mir-375-3p has been reported to be enriched in islets and the mirna with the highest intra-islet expression  <cit>  and in our study was increased from 4–24 h in the  <dig> and 50 μg/kg groups, returning to approximately vehicle level by 48 h. mir-216a-5p and mir-217-5p remained elevated in the serum of rats longer than amylase or lipase and had a much greater dynamic range which could be advantageous if detection of pancreatic injury is not able to be examined at earlier time points.
dogs treated with caerulein displayed serum amylase and lipase increases in the 15 μg/kg group from 1–8 h while showing increases in the 50 μg/kg group from 1–24 h. mir-216a-5p levels increased prior to amylase and lipase and correlated well to amylase and lipase in the  <dig> and 50 μg/kg groups, but remained elevated longer and to a much larger degree than amylase and lipase.
mir-216b displayed similar kinetics to mir-216a-5p, but with a reduced dynamic range while mir-217-5p displayed increases similar to amylase.
interestingly, mir-217-5p displayed very large increases for the longest amount of time in the rat and this result was not reflected in the dog.
this difference could reflect species differences in the toxicity of the compound, differences in the relative amount of the mirna in dog pancreas or differences in the stability or clearance of the mirna in dog serum.
mirs-141-3p and 200c-3p appeared to increase prior to amylase or lipase in the 3 μg/kg group and displayed increases similar to amylase in the  <dig> and 45 μg/kg groups while mir-148a-3p displayed increases similar to amylase.
mir-375-3p was increased prior to amylase and lipase at 3 μg/kg, correlated well with amylase and lipase in the  <dig> and 45 μg/kg dose groups and had a much larger dynamic range in the dog study as compared to the rat study.
the other conserved pancreas enriched micrornas tested were not correlated to clinical chemistry or histopathologic changes.
these observations are important and demonstrate that mirnas may display different kinetics in different species.
therefore, different mirnas may have more favorable characteristics as biomarkers in different species and should be thoroughly investigated in each.
interestingly, liver enriched mir-122-5p was increased in a dose and time dependent manner while ast and alt were increased in the 45ug/kg treated dogs with no observed histopatholgic correlate.
the rat caerulein study did not generate similar mir- <dig> increases.
it is possible that, due to species differences in injury and physical proximity, pancreatic damage induced by caerulein in dogs caused some perturbation of the membranes of hepatocytes resulting in ast, alt and mir-122-5p increases in the serum while not producing similar results in the rat.
additionally, individual rats and dogs displayed large differences in the amount of mirna increased in the serum.
this may be due to individual animal variation in the pancreatic damage, differing stability of the mirnas or lack of mirna normalization similar to tubulin or gapdh used for mrna comparisons as all samples were normalized by serum volume.
we tested several conserved pancreas enriched mirnas in rat and dog pancreas toxicity studies and determined that the mirnas perform as well or better than traditional markers of pancreas injury.
based on these studies mir-216a-5p displayed the largest dynamic range and correlated well with amylase and lipase indicating that it is the best translatable mirna biomarker that correlated to pancreatic injury from rat to dog.
it is possible that additional pancreas specific and/or enriched mirnas may supplement the interpretation of pancreatic injury or serve as better biomarkers than mir-216a-5p in pancreatic toxicity studies utilizing different pancreatic toxicants.
the data generated in this study verified our rationale for generating a library of rat tissue/organ specific and enriched mirnas.
these data can serve as a resource for future research by the scientific community for discovery and characterization of additional mirna biomarkers as well as a resource for more comprehensive isomir analysis.
additional considerations
while the work here has advanced the knowledge of the mirna and isomir content of rat organs it is important to note that the quantitative capacity of mirna sequencing may be inaccurate due to ligase bias as certain mirnas may be under or over-represented and potentially not detected .
the sequencing data and subsequent analysis did identify many tissue enriched mirnas previously discovered through other methods, therefore it is reasonable to infer that ligase bias may not be dependent on the rna content of a particular sample, but may be dependent on the base composition of an individual mirna.
ligase bias most certainly affected the quantitative capacity of mirna sequencing in this study  as the quantitation of mirs-216a-5p, 375-3p and 148a-3p did differ between the sequencing data and qpcr.
it is also possible that ligase bias could affect quantitation of isomirs and this could affect which mirna sequences are targeted for detection in the serum by qpcr or other technologies.
qpcr assays in the caerulein studies conducted here were chosen based upon the sequences that aligned to the mature mirna with the highest expression level.
this may have resulted in picking an isomir sequence that was expressed at a different level than the sequencing data indicated.
if a mirna is not expressed at a high level in a tissue then it would be reasonable to conclude that little would be detectable in the serum and this could be an additional reason why some mirnas did not appear to perform as well as others in the toxicology studies conducted here.
in addition, as methods for improving ligase bias improve, reanalysis of the mirna content of rat organs may need to be performed or analytical models may be developed that could account for ligase bias to gain a true understanding of mirna quantity in each organ and potentially detect additional mirnas.
isomirs generally displayed more variation on the 3′ ends which can affect the selection of the most appropriate qpcr assay to use for detection, therefore selection of the sequence present in the greatest amount should result in optimal detection of injury if the mirna is released upon cellular injury, stable in the blood, present in sufficient amount and ligase bias is not a confounding factor.
improvements in the quantitative capacity of mirna sequencing could also result in changing which mirna is defined as the mature mirna sequence in mirbase and other mirna cataloguing databases.
an advantage of mirnas as biomarkers of organ injury is the potential for species translatability of the biomarker for organ toxicity using the same assay to detect the mirna.
it is possible that expression of mirs and isomirs could change in the tissues of compound treated animals.
therefore, it would be worthwhile to determine if the tissue specificity/enrichment and expression levels of mirna sequences change in response to toxic insult.
ultimately, toxicology studies examining the kinetics of tissue specific and enriched mirnas utilizing a variety of compounds eliciting toxicity in a particular organ in different species will likely be needed to establish which micrornas are the best species translatable serum based biomarkers of organ injury.
studies may also need to be conducted with a wide variety of toxicants to determine which mirnas are subject to the least amount of change in expression level in serum in order to determine if a serum mirna could be used for normalization.
additionally, determining reference ranges for mirnas of interest will be a topic that will need to be investigated with additional toxicology studies in the future.
this could assist researchers in determining what degree of change correlates to what degree of organ injury.
clearly, mirnas will need to be rigorously evaluated in a manner similar to other biomarkers  <cit> .
while the major organs of toxicologic interest were examined, not all tissues in the rat or dog were examined.
therefore, tissues such as the thyroid, parathyroid, lung, spleen, fat, bone, cartilage, skin, spinal cord, peripheral neurons and other tissues could also be examined in future studies for their mirna content.
moreover, as shown in this study, different data analysis methods can generate different numbers of tissue enriched mirnas.
these data may yield additional candidate biomarkers of organ injury if different data analysis methods are employed.
as further characterization of serum mirna biomarkers is conducted, a qpcr based panel or a method utilizing multiplexing of mirnas representative of each organ with multiple mirnas per organ may be generated and used regularly in toxicology studies utilizing biomarkers as the major or a complimentary method for determining organ toxicity.
this method could ultimately detect toxicity in organs with validated mirna biomarkers using a small volume of serum and may significantly reduce the cost of preclinical and clinical drug development.
CONCLUSIONS
using a next generation sequencing approach, we have generated a rat mirna body atlas.
in that process we discovered an additional  <dig>  potential rat mirnas, verified previously identified tissue specific and enriched mirnas, identified many new tissue specific and enriched mirnas, and analyzed isomir expression of tissue specific and enriched mirnas.
in addition, we provided proof of concept that pancreas enriched mirnas conserved between rat and dog could be used as serum based markers of pancreatic injury in rat and dog caerulein toxicology studies.
in collaboration with takeda and hesi we have analyzed the dog mirna sequencing data and determined which tissue specific and enriched mirnas were conserved between the  <dig> species.
finally, we determined that mir-216a-5p appears to perform as well or better as a marker of pancreatic injury than amylase, lipase or other pancreas enriched mirnas in rat and dog caerulein toxicity studies.
additional files
additional file 1: supplementary methods.
eli lilly mirdeep  <dig> novel mirna identification
additional file 2: sequences of mirnas and exiqon qpcr assay identification.
additional file 3: table s <dig>  list of mirnas in common in hesi analysis.
tissue enriched mirnas mutually identified by independent data analysis techniques employed by eli lilly, niehs and maastricht university are listed including tissue  in which the mirna is most abundantly expressed.
sequences of mirnas and exiqon qpcr assay identification.
additional file 4: table s <dig>  tissue specific and enriched mirnas identified by lilly.
tissue specific and tissue enriched mirnas identified by eli lilly are listed under the organ in which they are most abundantly expressed.
tissue specific mirnas are in bold text and tissue enriched mirnas are in non-bold text.
additional file 5: figure s <dig>  tissue specific and enriched mirnas identified by lilly.
tissue specific and enriched mirnas identified by eli lilly are displayed in a heat map with tissues listed on the x-axis and mirnas listed on the y-axis.
expression levels are indicated by a range of colors with blue indicating low expression and red indicating high expression.
p-values are included next to each mirna indicating the degree of tissue specificity and/or enrichment.
some mirnas may be enriched in several tissues resulting in a high p-value.
additional file 6: table s <dig>  tissue enriched mirnas identified by niehs.
tissue enriched mirnas identified by niehs are listed under the organ in which they are most abundantly expressed.
additional file 7: table s <dig>  tissue enriched mirnas identified by maastricht university.
tissue enriched mirnas identified by maastricht university are listed under the organ in which they are most abundantly expressed.
additional file 8: table s <dig>  mirnas tested in qpcr validation.
tissue enriched mirnas and their sequences are displayed.
mirna abundance was measured through qpcr in order to correlate their tissue enrichment status to deep sequencing results.
additional file 9: table s <dig>  correlations of mirnas to sequencing data.
correlation of tissue enriched mirna measurements using deep sequencing and qpcr are depicted with a correlation of − <dig>  being highly correlated and − <dig>  representing poor correlation.
additional file 10: figure s <dig>  isomirs of pancreas enriched mirnas.
the isomirs of pancreas enriched mirnas and their sequences are displayed in a heat map with tissues listed on the x-axis and isomirs listed on the y-axis.
expression levels are indicated by a range of colors with blue indicating low expression and red indicating high expression.
isomirs of mir-217-5p display more expanded or restricted expression with respect to tissues.
additional file 11: figure s <dig>  isomirs of mir- <dig> in the intestine.
the isomirs of select intestine enriched mirnas and their sequences are displayed in a heat map with tissues listed on the x-axis and isomirs listed on the y-axis.
expression levels are indicated by a range of colors with blue indicating low expression and red indicating high expression.
isomirs of mir- <dig> display more expanded or restricted expression with respect to tissues.
additional file 12: figure s <dig>  isomirs of mir-192-5p in the liver.
the isomirs of select liver enriched mirnas and their sequences are displayed in a heat map with tissues listed on the x-axis and isomirs listed on the y-axis.
expression levels are indicated by a range of colors with blue indicating low expression and red indicating high expression.
isomirs of mir-192-5p display more expanded or restricted expression with respect to tissues.
additional file 13: figure s <dig>  mir-3473f-pre isomirs in the testis.
the isomirs of select testis enriched mirnas and their sequences are displayed in a heat map with tissues listed on the x-axis and isomirs listed on the y-axis.
expression levels are indicated by a range of colors with blue indicating low expression and red indicating high expression.
isomirs of mir-3473f-pre display more expanded or restricted expression with respect to tissues.
additional file 14: figure s <dig>  mir-3473f-pre isomirs in the testis.
the isomirs of select testis enriched mirnas and their sequences are displayed in a heat map with tissues listed on the x-axis and isomirs listed on the y-axis.
expression levels are indicated by a range of colors with blue indicating low expression and red indicating high expression.
isomirs of mir-3473f-pre display more expanded or restricted expression with respect to tissues.
additional file 15: table s <dig>  p-values of rat caerulein toxicity study.
statistical analysis of mirna and clinical chemistry biomarkers measured in the rat caerulein toxicity study are presented in this table with their associated fold changes and p-values indicated at 1 h , 4 h , 8 h , 24 h  and 48 h .
additional file 16: table s <dig>
p-values of dog caerulein toxicity study.
statistical analysis of mirna biomarkers measured in the dog caerulein toxicity study are displayed in this table with their associated p-values, lsmean, lower confidence intervals, upper confidence intervals and t ratios indicated at 1 h , 4 h , 8 h , 24 h  and 48 h  for the each dose group.
additional file 17: figure s <dig>  representation of rat mirna sequencing data.
rat mirna sequencing data is represented with the normalized reads per million on the y-axis and the tissues on the x-axis.
each color within the figure represents an individual mirna while the height of the colored bar represents the percentage of reads that a particular mirna makes up within that tissue.
mirs-216a-5p, 375-3p and 148a-3p are indicated in the yellow, green and pink bars in the pancreas tissue respectively.
additional file 18: figure s <dig>  qpcr of pancreas enriched mirnas.
total rna from the pancreas tissue of each individual rat was normalized by mass and examined by qpcr for expression of mirs-216a-5p, 375-3p and 148a-3p.
ct values are indicated on the y-axis and the rat pancreas samples are indicated on the x-axis.
each mirna appears to be expressed at approximately the same ct value indicating they are present at roughly the same quantity in the pancreas.
abbreviations
alpalkaline phosphatase
astaspartate aminotransferase
altalanine aminotransferase
bunblood urea nitrogen
cnscentral nervous system
co2carbon dioxoide
creatinineurinary albumin
kim-1kidney injury molecule-1
ctnicardiac troponin i
ggtgamma-glutamyl transferase
gldhglutamate dehydrogenase
hesihealth and environmental sciences institute
hmgb1high mobility group protein b1
mirmicrorna
ngalneutrophil gelatinase associated lipocalin
qpcrquantitative polymerase chain reaction
riscrna induced silencing complex
rnaribonucleic acid
sdsprague dawley
tmpd <dig> , <dig>  tetramethyl-p-phenylenediamine
<dig> this research was performed in collaboration with hesi genomics committee

 the drug discovery process is a long and expensive endeavor with many experimental molecules failing in late stage in vivo toxicology studies. the high attrition rate of drugs in late stage clinical trials lead to the call for investigation of the utilization of species translatable biomarkers to select safer compounds early in the drug development process. in order to better understand toxicities in preclinical and clinical settings and to speed the development of new therapies, body fluid based biomarkers of organ damage have been sought to both identify and monitor toxicities. serum or plasma based biomarkers capable of monitoring tissues affected by drug toxicity could allow for the detection of tissue injury from simple blood draws, thus permitting larger numbers of compounds to be tested more rapidly in studies that are less expensive and use fewer animals compared to traditional means of toxicity evaluation. the properties of an ideal biomarker include tissue specificity, a strong correlation to organ dysfunction or toxicity, a strong association to histopathologic changes, robust sensitivity, predictability, accessibility through a noninvasive method, readily-available rapid and inexpensive assays for quantitative measurement and direct translatability across preclinical species and human. serum based biomarkers of organ dysfunction have been useful in preclinical and clinical settings. serum ast and alt were discovered in and have been used in combination with total bilirubin, gamma-glutamyl transferase, alkaline phosphatase, albumin and bile acids as markers of liver injury to de-risk compounds in preclinical drug development as well as being utilized for monitoring liver toxicity in clinical drug development. additional markers of hepatotoxicity are currently under investigation. cardiac troponin i is a sensitive and specific serum based biomarker of cardiac injury in humans and mice and has been employed to mitigate cardiac toxicity in mouse toxicology studies with novel kinase inhibitors. urine based biomarkers for evaluation of renal toxicity include bun, creatinine, urinary albumin, kim-, ngal and cystatin c. amylase and lipase have been used as markers for pancreatitis, but their diagnostic capability is limited by sensitivity and specificity. these preclinical and clinical biomarkers are carefully reviewed in dieterle,. despite the success with protein based biomarkers to date, gaps in our ability to detect toxicities in a number of organs prevent the widespread utilization of biomarkers in toxicology studies. current gaps in the minimally-invasive biomarker portfolio include gastrointestinal injury, cns injury, reproductive toxicity and specific markers of glomerular and liver injury. clinical chemistry parameters offer important information regarding the health of a given organ but can be limited by lack of specificity and sensitivity as in the cases of alt and ast for liver injury and amylase and lipase for injury in the exocrine pancreas. while protein based biomarkers provide useful information, they can be limited in their translation across species and to the clinical setting due to species differences in the biomarkers and the need for species specific immunoassay reagents. traditional pre-clinical methods for assessing the potential toxicity of a compound include clinical chemistry, histopathology, hematology, physiological measurements and clinical observations, but the concordance rate to human toxicity was reported by olsen et al. to be 71 % leaving room for significant improvement in preclinical toxicology studies. unfortunately, the concordance of clinical chemistry or other biomarkers was not directly compared to establish concordance between preclinical and clinical studies; several groups advised more use of biomarkers in preclinical and clinical studies. in a subsequent study, adverse drug reactions in each target organ from compounds approved in japan from 2001– were found to be 48 % concordant between preclinical toxicology studies and humans while changes in laboratory parameters displayed a 30-70 % correlation. it is evident that additional endpoints indicative of organ toxicities that are translatable and/or predictive from preclinical species to human are necessary to improve the safety of compounds in preclinical drug development and clinical trials. with the end goal of identifying tissue specific, species translatable blood based biomarkers of tissue injury, we set out to characterize the rat microrna expression profile in various organs and create a rat microrna body atlas to identify specific and enriched mirnas in individual tissues. micrornas are 21– nucleotide long rna molecules whose cellular function appears to be to inhibit translation or cause degradation of their mrna target or targets. some mirnas are expressed in high abundance in a tissue specific manner, are resistant to degradation in the serum, are well conserved in preclinical species including mouse, rat, macaque and human, and are detectable in a wide variety of biologic fluids. in addition, sensitive, specific and readily available qpcr based assays exist to measure mirnas making them attractive candidates for serum based biomarkers for organ injury. perhaps the most well characterized mirna for use as a toxicity biomarker is the liver enriched mir-122-5p. this mirna was increased in the serum of rats that were given the liver toxicant trichlorobromomethane, but not increased in the serum of rats administered the skeletal muscle toxicants, tetramethyl-p-phenylenediamine and an hmg-coa reductase inhibitor. ast and alt were increased when animals were given liver or skeletal muscle toxicants, demonstrating that mir- was indicative of liver injury only and not muscle injury, thus displaying an advantage in specificity over alt and ast. mice treated with acetaminophen displayed increases of mir- in the serum 2 h before an increase of ast/alt. additionally, mir- and liver enriched mir- were increased in the serum of patients who had acetaminophen induced hepatotoxicity demonstrating the potential for these mirnas to be used as both preclinical and clinical biomarkers for liver injury. in a further study mir- hmgb cytokeratin and gldh were able to predict patients who would develop liver failure as a result of acetaminophen overdose even when alt remained in the normal range. while several mirna tissue profiling experiments have been reported in mouse, rat and human, these rodent experiments do not encompass the range of tissues relevant to toxicology studies for drug candidates. in addition, published studies relied upon microarray and/or qpcr based technologies which could suffer from a lack of specificity and sensitivity. microarrays also do not allow for the identification of novel mirnas due to their dependence on known sequences and may not account for tissue expression of isomirs, which are mirnas that have variations from the reference mirna sequence. qpcr also lacks specificity in cases where mirnas have a closely aligned sequence. as an example, mir- expressed in skeletal muscle and heart and mir- expressed specifically in skeletal muscle differ by only bases. several of the current qpcr based assays for mir- also detect synthetic mir- indicating that caution must be taken when using this approach. next generation sequencing is a preferred method for mirna profiling because it makes direct, unbiased and quantitative measures of known and unknown mirnas. in order to realize the true potential of mirnas as biomarkers of organ damage it is critical to comprehensively understand the tissue expression pattern of mirs and isomirs, and this can be best achieved using ngs. therefore we have constructed a rat mirna body atlas using illumina mirna sequencing from male and female sprague dawley rats and and organs of toxicologic interest from each sex, respectively. in parallel to this effort, several of the authors have constructed a similar canine body atlas from tissues from male marshall beagle dogs. publication of both mirna body atlases will provide an invaluable resource for the characterization of potential biomarkers of organ injury in most commonly used species for preclinical safety assessment. here we present the findings of the rat mirna body atlas in which we have identified previously unreported mirnas in the rat and detailed numerous tissue specific and enriched mirnas. we have discovered several novel tissue specific mirnas in each species, some of which are conserved across both rat and dog, and some which are species specific. comparison of the rat and dog atlases yielded candidate liver and pancreas specific biomarkers. the potential pancreas biomarkers were tested in rat and dog toxicology studies using the pancreatic toxicant caerulein, while the potential liver biomarkers were tested in toxicology studies conducted by takeda. in these proof-of-concept studies, mirnas performed as well as or better than traditional markers of pancreatic injury, demonstrated a much larger dynamic range and increased tissue specificity over traditional biomarkers methods necropsy five male and five female sprague dawley rats 12– weeks in age were obtained from charles rivers laboratories and euthanized by co asphyxiation. the rats were necropsied and tissues were placed into rnalater® within 5 min of asphyxiation. organs were processed in batches according to approximate tissue size as follows and placed into rnalater® at 4 °c for 24 h and were transferred to-20c the following day. in batch # greater than 150 mg of liver, stomach, ileum, jejunum and duodenum were placed into rnalater® for subsequent total rna isolation. the following portions of the intestine corresponding to the segment that is routinely collected for histology were collected: the stomach with approximately 5 cm of the duodenum attached, a 10 cm section from the midpoint of the jejunum and the distal 5 cm of ileum. a syringe containing rnalater® was used to flush the contents from the lumen of the collected segments prior to placing them into the rnalater® collection tube. in batch # ~150 mg of kidney, pancreas, brain, testis, biceps and soleus were placed into rnalater® for subsequent total rna isolation. in batch # 30-50 mg of adrenal, heart, ovary and uterus were placed into rnalater® for subsequent total rna isolation. rna extraction tissues were removed from rnalater®, frozen in liquid nitrogen and pulverized according to the covaris tissue pulverization protocol. tissues were placed in an appropriate volume of qiazol and homogenized using lysing d matrix tubes. qiazol volume was scaled up according to the mass of tissue to be homogenized. total rna was next isolated and eluted into 80ul of 95 °c molecular biology grade water. rna quality was assessed using the agilent bioanalyzer using the rna pico kit. all rnas except for pancreas, 3/ hippocampus samples and uterus sample had rin values of or higher with an average rin value of. for samples that displayed rin values lower than the rna extraction protocol was repeated using additional qiazol/tissue homogenate. most samples achieved rin values greater than as shown in fig.  whole blood was placed into rnalater® for blood and all rnas were isolated using the mirneasy kit from qiagen following the manufacturer’s instructions. dorsal root ganglion rnas were provided by takeda pharmaceuticals. rnas were quantitated using a nanodrop and, if necessary, were concentrated by speed vac to obtain an appropriate concentration for library preparation prior to sequencing.fig. the tissues used for construction of the rat mirna body atlas and boxplots of the rin values associated with the tissues from male and female sprague dawley rats are illustrated above. rnas were extracted from rat tissues and were tested for rna quality using the agilent bioanalyzer. most rnas were of very high quality as judged by rin values. pancreas rnas had the lowest rin values and also had the fewest sequences that mapped to mature mirnas mirna seq library construction library construction was performed according to illumina’s truseq small rna sample prep protocol. library quality was assessed using the agilent bioanalyzer dna chip and quantified using ribo green. libraries were subjected to miseq analysis to ensure library quality. body atlas libraries were assigned to gels and lanes on the flow cells in a manner designed to mitigate batch effect. subsequent hiseq analysis was performed and 50 bp single end reads were generated on the illumina hiseq platform resulting in 4– million reads per sample. data analysis by institution read mapping and quantification data was analyzed by eli lilly, maastricht university and niehs using separate analysis methods. the output generated from each institution was then collated to create a venn diagram of the tissue enriched mirnas identified by each analysis method. read mapping and quantification for the rat and dog mirna body atlases was performed in the following manner by eli lilly: fastq files were processed to remove adaptor sequences and discard reads < 17 bp in length with cutadapt v.. with options –a tggaattctcgggtgccaagg–quality–base 33–q 20–match-read-wildcards –m all trimmed reads that contained an ‘n’ were discarded. identical sequences from the same sample were combined into a single sequence in the form expected by mirdeep. that is, with ids of the form “nnn_#_x#”, where nnn is a unique alphabetic code identifying the sample, the first ‘#’ is a unique integer, and ‘x#’ records the number of times this sequence occurred in the input for the sample. quantifier.pl from the mirdeep package v.. was used to generate mir alignment files against known mirs from mirbase with options–p–p < org > _hairpin.fa –m < org > _mature.fa–r trimmed_reads.fa–d. this was done separately for rat, mouse, human, and c. elegans known mirs. the.mrd files were then parsed with a custom perl script as follows. each isomir sequence in an alignment was associated with the corresponding mature mir identifier. if it aligned to a mirna precursor, but not with an expected mature mir sequence, it was identified as < mir > −pre to indicate this. usually these correspond with reverse strand mirs that have not yet been annotated as alternate mature forms. sometimes these appear to be microrna offset rnas, and sometimes they appear to represent incompletely processed sequence. if a given sequence was identified as aligning with more than one precursor, it was associated with all potential names as a composite name. that is, a sequence that aligned to both let-7a-5p and let-7f-5p was assigned to the composite mature mir let-7a-5p;let-7f-5p to indicate for subsequent analyses that this identification was ambiguous. all sequences that had not yet been identified against known rat mirs were then looked for in the analysis with respect to known mouse mirs, and so on against human and finally c. elegans mirs. this process allows for identification of conserved rodent mirs that have not yet been annotated in rat, and against conserved mammalian mirs that have not yet been annotated in rat and mouse. in parallel to the above identification of known mirs, we used mirdeep2’s novel mir identification process. additional details are described in additional file 1: supplementary methods. tissue specific and enriched mirna identification data analysis of the rat and dog sequencing data was performed by eli lilly as follows. tissue specific mirnas are defined as mirnas that are expressed at a raw count of greater than in out of rats in organ/ tissue only. for example, the brain may be the only organ with mir-x expression, but the hypothalamus may be the only tissue that has mir-x expressed. therefore, mir-x would be specific to the hypothalamus tissue. tissue enriched mirnas are defined as mirnas that are expressed at a raw count of greater than in out of rats in a few organs/tissues, typically 3– tissues. the definition of tissue enriched was purposely vague because a mirna may be enriched in several organs, but could still be a useful biomarker if toxicity is only present in of the organs the mirna is enriched in. sex specific mirnas are identified as mirnas that are expressed in tissue/organ at a raw count of greater than in 4/ males or females and in 0/ rats in the opposite sex. raw and normalized data were analyzed for tissue specific and tissue enriched mirnas. animal level data is obtained after preprocessing of the raw sequencing data. for each microrna, there are measurements from animals and potentially different tissues. these tissues were from organs. in order to find tissue specific mirnas, tissue level counts were obtained for each microrna by summing over all the animals. the tissue level counts were aggregated to organ level by selecting the maximum of the tissue counts to represent the organ count. for example, brain has different tissue types, the maximum of these tissue level counts will be used to represent the expression levels of brain. after these data processing steps, a count data matrix is generated, where each row corresponds to a mir and each column corresponds to an organ type. non-negative matrix factorization is a commonly used unsupervised learning algorithm to decompose a nonnegative matrix x into non-negative matrices w and h. in our case, x is the observed count matrix of dimension nx where n is the number of mirs. each column of w defines a mir group and each column of h represents the expression pattern of the mir group corresponding to the organ type. based on this decomposition, we can find a representative set of mirs that corresponds to the specific high expression pattern for the organ. to identify tissue specific mirs, it has to be organ specific. thus we start from each organ specific mir, and examine if it is also highly expressed in a particular type of tissue or a few tissues, which belong to the current organ. for example, a mirna may be specific to brain, but it may only be specific to the hippocampus. to determine the expression level of each mir, we fitted a two-component mixture of poisson distribution to the animal level tissue counts data. the component with larger intensity of the poisson distribution is naturally used to model high counts and the component with lower intensity is used to model background noise. the larger component of the poisson mixture model naturally corresponds to high expression level, whereas the smaller component of the poisson mixture model naturally corresponds to the low expression level. for an organ specific microrna to be tissue specific, it has to be expressed at a high level in tissue of this organ. after the mixture model is fitted and the intensities are estimated, we can determine which tissues the microrna is expressed in according to their estimated intensities and then declare tissue specific micrornas afterwards. additionally, the raw sequencing data was manually analyzed in each tissue for each mirna and statistical analysis of the manually identified tissue specific and enriched mirnas was then conducted to verify or dispute the tissue specificity or enrichment. read mapping and quantitation by maastricht university read mapping and quantification was performed in the following manner by maastricht university: mirdeep was used to predict putative novel mirnas. the raw fastq files from the rat samples were mapped to the rat genome. mirdeep output was then parsed to keep only the predicted mirna precursors with a score of or above. raw reads were trimmed off the 3′ adapter sequence using the java algorithm “adaptor remover” from the freely available uea srna workbench. the parameters were set to look for a perfect match with the first bases of the 3′ adapter to be considered a hit. the trimmed sequences with a size below or above 35 bp were discarded. during this trimming step, every identical sequence was pulled and each iteration was summed. the trimmed reads were then mapped to the list of rat precursor mirna obtained either from mirbase, or generated de novo from the mirdeep prediction using a fast short read mapping software patman, configured to allow no mismatches or gap. patman output was then parsed to 1) assign a unique name to each unique sequence consisting of the pre-mirna name, the pre-mirna arms and the mapping coordinates. 2) divide the total read count of each unique sequence by the number of assigned loci in the complete mapping. 3) create a mature output table by pulling the different isomir count for each mirna species. the parsed patman output were normalized with the trimmed mean of m-values method and filtered to remove all mirnas where the normalized read counts were fewer than reads in all the samples. total read counts were then computed for each individual mirna. tissue enriched mirna identification by maastricht university a mirna was called tissue specific or tissue enriched when the percentage of reads for a single tissue was respectively higher than or 50 % of the total read count. read mapping and quantitation by niehs in order to identify isomirs the raw count number of each given isomir was first converted to the percentage of expression compared to the mature count observed for its corresponding mirna. by comparing this percentage, we were able to report the mirna for which the most expressed isomirs differs between tissues. analysis of the data performed at niehs was executed in the following manner: raw fastq files with the mirna-seq reads were checked for quality using a tqe standard based on trimming the adapters, quality filtering at q ≤  and the elimination of short reads <  bases long. in some cases when minor adapter sequences were observed after initial process, additional rounds of processing were applied until all the criteria were met. in the end, reads longer than base pairs were also eliminated, resulting in quality filtered reads between 14–24 bp. reads passing the tqe filtering were then aligned to the rat mirbase database v. using bwa version.1-r104q with “-n 2” and default parameters otherwise. only perfect matches were further used in the quantification, which was performed by summarizing read count level measurements at each mature mirbase mirs. datapoints with zeros were imputed with. total counts for each mir were normalized by transcripts per million reads. prior to statistical analysis, a transformation from floats to integers was performed by ceiling the data. tissue enriched mirna identification by niehs to detect the tissue-enriched micrornas, we adopted the one-vs-rest strategy when applying statistic models. “differentially” expressed mirs in tissue against the combination of all other tissues were identified using a quasi-poisson approach with shrunken dispersion estimates. during the modeling process, the offset was computed on the 3rd quartile. significant demirs were detected at a nominal p-value < . rat and dog toxicity studies caerulein studies: studies were run at covance laboratories in greenfield, in and histopathology was performed by a board certified pathologist. due to the short half-life of caerulein, male sprague dawley rats per group were treated with intraperitoneal doses of vehicle or caerulein, 1 h apart, at and 50 μg/kg. blood was collected at and 48 h after dosing and examined for clinical chemistry, hematology and microrna changes with terminal serum taken at and 48 h. urine was collected by cystocentesis at and 48 h and examined for total protein and specific gravity. toxicokinetic analysis was performed at, and 48 h post dose. histopathology of adrenal glands, bone, brain stem, cerebellum, cerebrum, duodenum, heart, ileum, jejunum, kidneys, liver, lung, muscle, pancreas, spleen, stomach and thymus was examined at and 48 h post dose. for canine studies, beagle dogs per group, male and female obtained from covance research products, were given doses h apart of caerulein at or μg/kg. blood was collected at pre-dose, and h after the third dose and examined for clinical chemistry, hematology and mirna changes. urine was collected by cystocentesis and examined for changes in specific gravity, bun, creatinine and other urinary parameters at necropsy at and 72 h. histopathology on organs including adrenal glands, bone, brain stem, cerebellum, cerebrum, colon, epididymis, eyes, heart, ileum, jejunum, kidneys, liver, lung, muscle, sciatic nerve, ovaries, pancreas, prostate, spleen, stomach, testis, thymus, thyroid, uterus and vagina was performed at and 72 h post dose. animals were euthanized by barbiturate overdose or anesthesia and exsanguination at necropsy. qpcr of serum samples from caerulein studies for each sample, 100 μl of serum from all studies was subjected to total rna isolation, including mirnas, using the mirneasy tm kit from qiagen. upon the addition of qiazol 5 μl of 5 μm synthetic cel-mir-55-3p was spiked into each sample in the rat caerulein and dog caerulein studies. rnas were eluted in 80ul of 95 °c h2o and a water extraction control with spike-in was also included. reverse transcription was executed with the mircury lna microrna universal cdna synthesis kit. 10 μl of the extracted rnas were used in 20 μl rt reactions. a no template control and two -rt controls were also included. real-time pcr was conducted using exiqon primers and exilent sybr green master mix. samples were diluted 1: prior to qpcr and 5 μl of the diluted samples were used in 10 μl reactions. data was collected on a 7900ht real-time pcr instrument. statistical analysis of rat and dog caerulein qpcr data for each study, the fold change data were log transformed and fitted to a two-way anova model with repeated measures. because each animal has multiple measurements, we allow these measurements to be auto-correlated. once we obtained the parameter estimates, for the dog study, all treatment means and 95 % confidence intervals at each time point are compared against the pre-dose group, which have fc fixed at for the rat study, low dose and high dose groups are compared to untreated group at each time point. qpcr validation of tissue specific and enriched mirna tissue enriched and tissue specific mirnas were chosen for qpcr validation. the rnas from rats and tissues were normalized by mass and subjected to cdna synthesis using exiqon universal cdna synthesis kits and exilent sybr green mastermix. data was collected on a 7900ht or viia real-time pcr instrument. assays were designed to the most highly expressed mirna or isomir sequence present that mapped to the mature mirna. for example, many isomirs mapped to mir-xyz. for qpcr analysis we designed primers to the most highly expressed mirna sequence that mapped to mir-xyz as determined by deep sequencing. ct values of each mirna were compared to one another to determine tissue specificity or enrichment. a ct value was considered acceptable if the ct in a particular tissue was cts lower than the signal in the no rt or no template controls unless there was no signal in the no rt or ntc. if an assay did not display a signal that was not cts away from the no rt or ntc it was not included in the analysis. additionally, if the melt curves for an assay displayed more than one peak the assay was not included in the analysis. qpcr analysis of pancreas enriched mirnas total rna from all male and female rats used in the study were normalized by mass using 250 ng of rna per sample and subjected to cdna synthesis using exiqon universal cdna synthesis kits and exilent sybr green mastermix. data was collected on a 7900ht or viia real-time pcr instrument.  our first objective was to identify mirnas expressed in rat tissues. members of the hesi consortium were queried on which rat organs are most commonly injured during toxicology studies and therefore, would represent the most important organs to investigate for potential biomarkers of injury. this resulted in the list of organs detailed in fig.  dorsal root ganglion was added as an organ of investigative interest. the brain lacks serum and cerebrospinal fluid based biomarkers of injury, therefore the mirna content of the cerebrum, cerebellum, hippocampus and brainstem were interrogated. additionally, whole blood was included in the study in order to determine what the contribution of blood mirnas may be to highly vascularized tissues. total rna was isolated from the tissues of male and female sd rats and the quality was assessed using rin values to be sure that small rnas would be sequenced as opposed to degraded rnas. the rat microrna sequencing data, as analyzed by eli lilly, generated sequences that mapped to previously-identified mature micrornas in rat, mouse and human. this resulted in an increase of mirnas in the rat as compared to the mature rat mirnas in mirbase v while eli lilly’s data analysis serves as the main source of data for discussion in this article we also collaborated with hesi members to analyze the data in an attempt to maximize our ability to discover and confirm the tissue enrichment profile of mirnas.fig. workflow diagram. the small rna content of rat organs was interrogated by deep sequencing and the data was given to institutions for independent data analysis. reads were mapped to rat, mouse and human mirnas from mirbase v and v using mirdeep and tissue specific and tissue enriched mirnas were identified our second objective was to discover novel and confirm previously identified tissue specific and/or tissue enriched mirnas and/or isomirs that could be used as serum based biomarkers of organ injury in the rat. a variety of data analysis software packages and statistical methods exist to identify and predict mirnas and perform differential expression on mirna-seq data. since these methods have different strengths and weaknesses and the most suitable data analysis method is dependent on the output required, different methods were employed by scientists at eli lilly, niehs and maastricht university to analyze the sequencing data in order to maximize discovery of tissue enriched mirnas. the tissue enriched mirnas identified by each institution were included in a venn diagram in order to determine which tissue enriched mirnas were found in common between each data analysis method. tissue specific mirnas identified by lilly are reported in table  and were not included in the venn diagram. maastricht university identified tissue enriched mirnas, eli lilly identified tissue enriched mirnas and niehs identified tissue enriched mirnas. maastricht university and eli lilly had mirnas in common, maastricht university and niehs had mirnas in common, eli lilly and niehs had mirnas in common and all institutions had tissue enriched mirnas in common. differences in the results from each institution are due to the processing steps prior to alignment, the reference sequences the sequencing data were aligned to and the methods of identifying tissue specific and/or enriched mirnas. despite which approach was taken each analysis identified many of the same tissue enriched mirnas and the data has been made publically available for researchers in industry, academia or government to explore and potentially discover additional tissue specific and/or enriched mirnas which may serve as novel biomarkers for safety, pharmacodynamic or other effects.fig. venn diagram of data analysis methods. analysis of the mirna sequencing data was independently performed by eli lilly, maastricht university and the niehs using different data analysis methods. eli lilly identified tissue enriched mirnas, maastricht university identified tissue enriched mirnas and the niehs identified tissue enriched mirnas. the number of mirnas in common between maastricht university and eli lilly was the number of mirnas in common between maastricht university and the niehs was and the number of mirnas in common between eli lilly and the niehs was the number of mirnas in common between all analysis was 165 tissue specific mirnas were defined as mirnas there are expressed at a raw count of > in 8/ rats in tissue only. sequencing data was analyzed by non-negative matrix factorization which identified tissue specific mirnas. organs not listed in the table did not have tissue specific mirnas tissue specific and enriched mirnas identified by lilly are shown in and and the tissue enriched mirnas identified by niehs and maastricht university are shown in. the tissue specificity/enrichment of many other tissue enriched mirnas including muscle and heart enriched mirs-133a-3p, 499-5p, mir-1a-3p as well as pancreas enriched mirs-216a-5p, 217-5p and mir-375-3p were confirmed. the brain, testis and liver have the most tissue specific and enriched mirnas and interestingly, tissue enriched mirnas in the brain are often expressed in the intestine. the kidney has the fewest tissue specific and enriched mirnas of any organ examined, while skeletal muscle and heart share several of the same tissue enriched mirnas. while previously identified tissue specific and enriched mirnas were verified, previously unidentified tissue specific/enriched mirnas were discovered including, but not limited to, pancreas enriched mirs-148a-3p, 141-3p and 200c-3p. in order to verify that mirnas were, in fact, tissue specific or enriched we conducted qpcr verification of some of the tissue specific and/or enriched mirnas. as shown in fig.  of the mirnas tested correlated very well to the sequencing data with a correlation of −  to − , mirnas displayed a correlation of −  to −  and mirnas displayed poorer correlations from −  to − . mirs-196c-5p and 206-5p were not amplified by qpcr and, therefore, did not display tissue specificity or enrichment by qpcr.fig. heat map of tissue specific and tissue enriched mirnas. examples of tissue specific and enriched mirnas identified by this and previously published studies are displayed with tissues on the x-axis and mirnas listed on the y-axis. mirna expression from least to highest expression is indicated as a shift from blue to red with blue representing the least expression and dark red indicating maximal expression fig. correlation of qpcr analysis of tissue specific and/or enriched mirnas. tissue specific and/or enriched mirnas were chosen according to tissues that lack or require additional biomarkers of organ status and were measured in tissues from female rats by qpcr. the qpcr data was compared to the sequencing raw counts for each mirna. a perfect correlation is equal to −  because the ct value decreases as sequencing abundance increases. twenty seven of the thirty nine mirnas tested correlated very well to the sequencing data with a correlation of −  to − , mirnas displayed a correlation of −  to −  and mirnas displayed poorer correlations from −  to − . mirs-196c-5p and 206-5p were not amplified by qpcr and, therefore, did not display tissue specificity or enrichment by qpcr while each institution was able to identify tissue specific and enriched mirnas it is also imperative to understand whether the isomirs of the tissue specific and enriched mirnas are expressed in the same tissues as the parent mirna as isomir expression in other tissues could lead to inaccurate interpretation of serum mirna profiles with respect to toxicity. eli lilly investigated the isomir expression of pancreas tissue specific and enriched mirnas and the analysis revealed that isomirs generally mirror their parent mirna expression, but some isomirs are more tissue specific than others as shown for mir- in the intestines and mir-217-5p in the pancreas. for example, the mature mir-217-5p sequence is 5′tactgcatcaggaactgactgg-3′ and the isomir-217-5p sequence 5′-tactgcatcaggaactgactggac-3′ are both highly enriched in the pancreas, but the mir-217-5p isomir lacks expression in the stomach, intestine and ovary while other isomirs of mir-217-5p are expressed in the stomach, intestine and ovary. isomirs in other tissues such as isomirs of mir- in the intestine, mir-192-5p in the liver and mir-3473f-pre in the testis display similar isomir expression. it appears that interpretation of serum mirna toxicity biomarker data should not be compromised by expression of isomirs from tissues lacking parent mirna expression. rat and dog pancreas toxicity studies conserved pancreas tissue enriched mirnas were examined in the rat and dog models of caerulein pancreatic toxicity in order to begin to characterize mirnas as species translatable serum based biomarkers of pancreatic injury. rats were treated with vehicle, and 50 μg/kg of caerulein and the serum was examined for clinical chemistry and mirna changes at and 48 h with histopathologic assessment of organs at and 48 h. minimal and slight pancreatic necrosis was noted in vehicle treated rats at and 48 h which had no clinical chemistry or mirna correlates and were considered to be consistent with common background while minimal, slight, moderate and marked pancreatic necrosis was noted in the and 50 μg/kg treated rats that was interpreted to be dose related. clinical chemistry parameters including amylase and lipase, mir-122-5p, mir-133a-3p, 148a-3p, 208a-3p and pancreas mirnas conserved between rat and dog were examined for changes in the serum. mir-216a-5p, amylase and lipase were increased in the serum concurrently at 1 h, remained elevated until 8 h while mir-216a-5p remained elevated until 24 h. mir-217-5p displayed similar kinetics to mir-216a-5p except mir-217-5p generally had a larger dynamic range and remained elevated until 48 h. both mirs-216a-5p and 217-5p displayed much larger dynamic ranges than amylase or lipase and remained elevated longer. mir-216b-5p displayed serum increases that were kinetically similar to amylase and lipase, but had a dynamic range that was less than lipase and greater than amylase. mir-375-3p was increased from 4–24 h in the and 50 μg/kg groups and returned to approximately vehicle level by 48 h. the pancreas enriched mirnas conserved between rat and dog that were tested displayed increases in the serum similar to amylase, whereas mirs 320-3p, and were not increased.table pancreas enriched mirnas conserved between rat and dog tacagtactgtgataactgatc tissue enriched mirnas conserved between rat and dog were identified and are shown above. these mirnas were analyzed for changes in the serum in rat and dog caerulein toxicology studies. bases underlined indicate bases that are not conserved between the rat and dog isomirs with the highest sequencing counts fig. histopathologic assessment of the pancreases of rats treated with and 50ug/kg of caerulein. histopathologic assessment of organs was performed at and 48 h with dose, time and animal number indicated on the x-axis and the grade of injury indicated on the y-axis. rat pancreases displayed minimal and slight pancreatic necrosis in vehicle treated rats at and 48 h while dose responsive increases in the severity of necrosis was noted in the and 50ug/kg treatred rats at and 48 h fig. clinical chemistry and serum mirna analysis of caerulein treated rats. rats treated with and 50ug/kg of caerulein were examined for clinical chemistry and mirna changes. time is indicated on the x-axis and the log of the fold change, the mean as well as 95 % confidence interval are indicated on the y-axis. mirnas were incereased in the serum of rats with and ug/kg of caerulein coincident with increases in amylase and lipase with maximal increases occurring at 8–24 h for mirs-216a-5p, 217-5p and 375-3p. values are normalized to vehicle and to the cel-mir-55-3p spike in marshall beagle dogs were treated with vehicle, and 45 μg/kg of caerulein and the serum was examined for clinical chemistry and mirna changes at and 72 h with histopathologic assessment of organs at and 72 h. pancreatic necrosis was noted in dog at 24 h in the 15 μg/kg group and and 72 h in the 45 μg/kg group. amylase and lipase were not increased greater than two fold in the vehicle or 3 μg/kg treated dogs at any time point. amylase and lipase were increased in the 15 μg/kg treated dogs with the largest increases at 1 h and decreasing by 8–24 h while both were increased maximally at 4-8 h in the 45 μg/kg treated group. mir-216b-5p was increased in vehicle treated animal and mir-217-5p was increased in vehicle treated animals with no histopathologic or clinical chemistry correlates. in 3 μg/kg treated dogs mir-216b-5p and mir-141-3p were increased at some time points, but mir-216b-5p was not increased beyond the level of the vehicle treated dog and did not display consistent time dependent increases. mir-141-3p and mir-216a-5p were increased beyond vehicle treated dogs and did not display consistent time dependent increases. mir-375-3p did display increases as compared to vehicle reaching its maximum increase at 1 h and remaining elevated until 24 h. in 15 μg/kg treated dogs mirs-216a-5p, 216b-5p, mir-375-3p and 141-3p were increased in a time dependent manner with maximal increases at 4–8 h. mir-216a-5p displayed the greatest dynamic range reaching as high as fold above vehicle. in the 45 μg/kg treated dogs amylase, lipase and mirnas displayed dose and time dependent increases. mir-216a-5p was increased maximally at 4–8 h and remained elevated until 72 h while levels as high as fold above vehicle were observed in individual animals. mir-216b-5p reached levels as high as fold above vehicle in individual animals and remained elevated until 24 h. the remaining mirnas did increase and remained elevated generally as long, but not to the degree of mirs-216a-5p and 216b-5p. interestingly, liver enriched mir-122-5p was increased in a dose and time dependent manner while ast and alt were increased in the 45 μg/kg treated dogs while no histopatholgic correlate was observed. the rat caerulein study did not generate similar mir- increases. as demonstrated by these studies, mirnas can serve as species translatable serum based biomarkers of pancreatic injury.fig. histopathologic assessment of the pancreases of dogs treated with caerulein. marshall beagle dogs were treated with and ug/kg of caerulein and histopathologic assessment of organs was performed at and 72 h with dose, time and animal number indicated on the x-axis and the grade of injury indicated on the y-axis. pancreatic degeneration/necrosis was not noted in dogs treated with 3ug/kg of caerulein. pancreatic degeneration/necrosis was noted in dog at 24 h in the 15ug/kg group and displayed dose dependent increases in severity at and 72 h in the 45ug/kg group fig. clinical chemistry and serum mirna analysis of caerulein treated dogs. dogs treated with and ug/kg of caerulein were examined for clinical chemistry and serum mirna changes. time is indicated on the x-axis and the log of the fold change, the mean as well as 95 % confidence interval are indicated on the y-axis.. mirnas were increased in the serum of dogs trreated with 3ug/kg of caerulein prior to increases in amylase and lipase. mirnas were dose dependently increased in the serum of dogs treated with and ug/kg of caerulein and correlated to amylase and lipase increases while mir-216a-5p remained elevated until 24 h. values are normalized to vehicle and to the cel-mir-55-3p spike in discussion we have generated a rat mirna body atlas and identified an additional ~ potential mirnas in the rat as compared to the current number of rat mirnas in mirbase v we have also confirmed many previously identified tissue specific and enriched mirnas and have discovered many previously unknown rat tissue specific and enriched mirnas. rat and dog mirna atlases were compared in order to determine which mirnas are conserved and may represent species translatable biomarkers of organ injury. subsequently, the conserved pancreas enriched mirnas in rats and dogs were demonstrated in pancreas toxicology studies to be sensitive blood-based biomarkers of caerulein induced pancreatic injury in both species. the aims of this study were fold. first we sought to comprehensively catalogue the mirna and isomir content of the major organs of toxicologic interest in the rat in order to discover potential mirna biomarkers of organ injury. secondly, we wanted to identify tissue specific and tissue enriched mirnas that may serve as biomarkers of organ injury. thirdly, we wanted to test tissue specific and/or enriched mirnas that are conserved in the rat and dog in toxicology studies to determine whether conserved mirnas could be used as serum based biomarkers of organ injury in both species. while others have generated similar data using microarray or qpcr, next generation sequencing is advantageous in that it allows for the discovery of mirna sequences without the need for prior genome annotation. analysis of the sequencing data refined the tissue specific/enrichment of previously reported tissue specific mirnas and revealed many previously undiscovered tissue specific and enriched mirnas. identification of tissue specific and enriched mirs and isomirs independent data analysis, utilizing different techniques, was performed by members of hesi. approximately tissue enriched mirnas were found to be in common between the data analysis methods and additional tissue enriched mirnas were identified by each institution. in order to ensure that interpretation of serum mirna profiles is reflective of organ specific toxicity we sought to determine the tissue specificity and enrichment of the isomirs of tissue specific and enriched mirnas identified by eli lilly. eli lilly identified many sequences that mapped to mature mirnas which consisted of sequences that differed from the mature mirna sequence. these sequences, termed isomirs contained additional bases on the 5′ and 3′ ends as compared to the reference mirna sequences and were also evaluated for tissue specificity and enrichment. we determined that the isomirs of tissue specific and enriched mirnas, with some exceptions, did generally mirror the expression of the reference mirna sequence indicating that serum mirna profiles should accurately reflect organ toxicity. some isomirs did display increased tissue specificity/enrichment as compared to the parent mirna and these could serve as more specific markers of organ injury when faced with multiple organ toxicity. eli lilly’s definition of a tissue enriched mirna was purposefully vague because some mirnas may be expressed in multiple organs at a higher level than the remaining organs, but may still be reasonable candidate biomarkers of organ injury depending on the observed toxicity. it is possible that a mirna that is enriched in several tissues could be a useful marker of tissue injury if toxicity is restricted to a single organ that the mirna is expressed in. in the face of multiple organ toxicity this approach may not be reasonable unless additional biomarkers could improve interpretation of the injury. consequently, some mirnas may not appear to be particularly tissue enriched as shown in additional file 5: figure s rat and dog pancreas toxicity studies in order to characterize pancreas enriched mirnas conserved between rat and dog as potential biomarkers of pancreatic injury we conducted rat and dog caerulein studies in which animals were treated with caerulein and concurrent clinical chemistry and histopathology were assessed along with mirna changes in time course studies. amylase and lipase increases were noted from 1–8 h in rats in both and 50 μg/kg dose groups while pancreatic necrosis was noted at and 48 h. mir-375-3p has been reported to be enriched in islets and the mirna with the highest intra-islet expression and in our study was increased from 4–24 h in the and 50 μg/kg groups, returning to approximately vehicle level by 48 h. mir-216a-5p and mir-217-5p remained elevated in the serum of rats longer than amylase or lipase and had a much greater dynamic range which could be advantageous if detection of pancreatic injury is not able to be examined at earlier time points. dogs treated with caerulein displayed serum amylase and lipase increases in the 15 μg/kg group from 1–8 h while showing increases in the 50 μg/kg group from 1–24 h. mir-216a-5p levels increased prior to amylase and lipase and correlated well to amylase and lipase in the and 50 μg/kg groups, but remained elevated longer and to a much larger degree than amylase and lipase. mir-216b displayed similar kinetics to mir-216a-5p, but with a reduced dynamic range while mir-217-5p displayed increases similar to amylase. interestingly, mir-217-5p displayed very large increases for the longest amount of time in the rat and this result was not reflected in the dog. this difference could reflect species differences in the toxicity of the compound, differences in the relative amount of the mirna in dog pancreas or differences in the stability or clearance of the mirna in dog serum. mirs-141-3p and 200c-3p appeared to increase prior to amylase or lipase in the 3 μg/kg group and displayed increases similar to amylase in the and 45 μg/kg groups while mir-148a-3p displayed increases similar to amylase. mir-375-3p was increased prior to amylase and lipase at 3 μg/kg, correlated well with amylase and lipase in the and 45 μg/kg dose groups and had a much larger dynamic range in the dog study as compared to the rat study. the other conserved pancreas enriched micrornas tested were not correlated to clinical chemistry or histopathologic changes. these observations are important and demonstrate that mirnas may display different kinetics in different species. therefore, different mirnas may have more favorable characteristics as biomarkers in different species and should be thoroughly investigated in each. interestingly, liver enriched mir-122-5p was increased in a dose and time dependent manner while ast and alt were increased in the 45ug/kg treated dogs with no observed histopatholgic correlate. the rat caerulein study did not generate similar mir- increases. it is possible that, due to species differences in injury and physical proximity, pancreatic damage induced by caerulein in dogs caused some perturbation of the membranes of hepatocytes resulting in ast, alt and mir-122-5p increases in the serum while not producing similar results in the rat. additionally, individual rats and dogs displayed large differences in the amount of mirna increased in the serum. this may be due to individual animal variation in the pancreatic damage, differing stability of the mirnas or lack of mirna normalization similar to tubulin or gapdh used for mrna comparisons as all samples were normalized by serum volume. we tested several conserved pancreas enriched mirnas in rat and dog pancreas toxicity studies and determined that the mirnas perform as well or better than traditional markers of pancreas injury. based on these studies mir-216a-5p displayed the largest dynamic range and correlated well with amylase and lipase indicating that it is the best translatable mirna biomarker that correlated to pancreatic injury from rat to dog. it is possible that additional pancreas specific and/or enriched mirnas may supplement the interpretation of pancreatic injury or serve as better biomarkers than mir-216a-5p in pancreatic toxicity studies utilizing different pancreatic toxicants. the data generated in this study verified our rationale for generating a library of rat tissue/organ specific and enriched mirnas. these data can serve as a resource for future research by the scientific community for discovery and characterization of additional mirna biomarkers as well as a resource for more comprehensive isomir analysis. additional considerations while the work here has advanced the knowledge of the mirna and isomir content of rat organs it is important to note that the quantitative capacity of mirna sequencing may be inaccurate due to ligase bias as certain mirnas may be under or over-represented and potentially not detected. the sequencing data and subsequent analysis did identify many tissue enriched mirnas previously discovered through other methods, therefore it is reasonable to infer that ligase bias may not be dependent on the rna content of a particular sample, but may be dependent on the base composition of an individual mirna. ligase bias most certainly affected the quantitative capacity of mirna sequencing in this study as the quantitation of mirs-216a-5p, 375-3p and 148a-3p did differ between the sequencing data and qpcr. it is also possible that ligase bias could affect quantitation of isomirs and this could affect which mirna sequences are targeted for detection in the serum by qpcr or other technologies. qpcr assays in the caerulein studies conducted here were chosen based upon the sequences that aligned to the mature mirna with the highest expression level. this may have resulted in picking an isomir sequence that was expressed at a different level than the sequencing data indicated. if a mirna is not expressed at a high level in a tissue then it would be reasonable to conclude that little would be detectable in the serum and this could be an additional reason why some mirnas did not appear to perform as well as others in the toxicology studies conducted here. in addition, as methods for improving ligase bias improve, reanalysis of the mirna content of rat organs may need to be performed or analytical models may be developed that could account for ligase bias to gain a true understanding of mirna quantity in each organ and potentially detect additional mirnas. isomirs generally displayed more variation on the 3′ ends which can affect the selection of the most appropriate qpcr assay to use for detection, therefore selection of the sequence present in the greatest amount should result in optimal detection of injury if the mirna is released upon cellular injury, stable in the blood, present in sufficient amount and ligase bias is not a confounding factor. improvements in the quantitative capacity of mirna sequencing could also result in changing which mirna is defined as the mature mirna sequence in mirbase and other mirna cataloguing databases. an advantage of mirnas as biomarkers of organ injury is the potential for species translatability of the biomarker for organ toxicity using the same assay to detect the mirna. it is possible that expression of mirs and isomirs could change in the tissues of compound treated animals. therefore, it would be worthwhile to determine if the tissue specificity/enrichment and expression levels of mirna sequences change in response to toxic insult. ultimately, toxicology studies examining the kinetics of tissue specific and enriched mirnas utilizing a variety of compounds eliciting toxicity in a particular organ in different species will likely be needed to establish which micrornas are the best species translatable serum based biomarkers of organ injury. studies may also need to be conducted with a wide variety of toxicants to determine which mirnas are subject to the least amount of change in expression level in serum in order to determine if a serum mirna could be used for normalization. additionally, determining reference ranges for mirnas of interest will be a topic that will need to be investigated with additional toxicology studies in the future. this could assist researchers in determining what degree of change correlates to what degree of organ injury. clearly, mirnas will need to be rigorously evaluated in a manner similar to other biomarkers. while the major organs of toxicologic interest were examined, not all tissues in the rat or dog were examined. therefore, tissues such as the thyroid, parathyroid, lung, spleen, fat, bone, cartilage, skin, spinal cord, peripheral neurons and other tissues could also be examined in future studies for their mirna content. moreover, as shown in this study, different data analysis methods can generate different numbers of tissue enriched mirnas. these data may yield additional candidate biomarkers of organ injury if different data analysis methods are employed. as further characterization of serum mirna biomarkers is conducted, a qpcr based panel or a method utilizing multiplexing of mirnas representative of each organ with multiple mirnas per organ may be generated and used regularly in toxicology studies utilizing biomarkers as the major or a complimentary method for determining organ toxicity. this method could ultimately detect toxicity in organs with validated mirna biomarkers using a small volume of serum and may significantly reduce the cost of preclinical and clinical drug development.  using a next generation sequencing approach, we have generated a rat mirna body atlas. in that process we discovered an additional potential rat mirnas, verified previously identified tissue specific and enriched mirnas, identified many new tissue specific and enriched mirnas, and analyzed isomir expression of tissue specific and enriched mirnas. in addition, we provided proof of concept that pancreas enriched mirnas conserved between rat and dog could be used as serum based markers of pancreatic injury in rat and dog caerulein toxicology studies. in collaboration with takeda and hesi we have analyzed the dog mirna sequencing data and determined which tissue specific and enriched mirnas were conserved between the species. finally, we determined that mir-216a-5p appears to perform as well or better as a marker of pancreatic injury than amylase, lipase or other pancreas enriched mirnas in rat and dog caerulein toxicity studies. additional files additional file 1: supplementary methods. eli lilly mirdeep novel mirna identification additional file 2: sequences of mirnas and exiqon qpcr assay identification. additional file 3: table s list of mirnas in common in hesi analysis. tissue enriched mirnas mutually identified by independent data analysis techniques employed by eli lilly, niehs and maastricht university are listed including tissue in which the mirna is most abundantly expressed. sequences of mirnas and exiqon qpcr assay identification. additional file 4: table s tissue specific and enriched mirnas identified by lilly. tissue specific and tissue enriched mirnas identified by eli lilly are listed under the organ in which they are most abundantly expressed. tissue specific mirnas are in bold text and tissue enriched mirnas are in non-bold text. additional file 5: figure s tissue specific and enriched mirnas identified by lilly. tissue specific and enriched mirnas identified by eli lilly are displayed in a heat map with tissues listed on the x-axis and mirnas listed on the y-axis. expression levels are indicated by a range of colors with blue indicating low expression and red indicating high expression. p-values are included next to each mirna indicating the degree of tissue specificity and/or enrichment. some mirnas may be enriched in several tissues resulting in a high p-value. additional file 6: table s tissue enriched mirnas identified by niehs. tissue enriched mirnas identified by niehs are listed under the organ in which they are most abundantly expressed. additional file 7: table s tissue enriched mirnas identified by maastricht university. tissue enriched mirnas identified by maastricht university are listed under the organ in which they are most abundantly expressed. additional file 8: table s mirnas tested in qpcr validation. tissue enriched mirnas and their sequences are displayed. mirna abundance was measured through qpcr in order to correlate their tissue enrichment status to deep sequencing results. additional file 9: table s correlations of mirnas to sequencing data. correlation of tissue enriched mirna measurements using deep sequencing and qpcr are depicted with a correlation of − being highly correlated and − representing poor correlation. additional file 10: figure s isomirs of pancreas enriched mirnas. the isomirs of pancreas enriched mirnas and their sequences are displayed in a heat map with tissues listed on the x-axis and isomirs listed on the y-axis. expression levels are indicated by a range of colors with blue indicating low expression and red indicating high expression. isomirs of mir-217-5p display more expanded or restricted expression with respect to tissues. additional file 11: figure s isomirs of mir- in the intestine. the isomirs of select intestine enriched mirnas and their sequences are displayed in a heat map with tissues listed on the x-axis and isomirs listed on the y-axis. expression levels are indicated by a range of colors with blue indicating low expression and red indicating high expression. isomirs of mir- display more expanded or restricted expression with respect to tissues. additional file 12: figure s isomirs of mir-192-5p in the liver. the isomirs of select liver enriched mirnas and their sequences are displayed in a heat map with tissues listed on the x-axis and isomirs listed on the y-axis. expression levels are indicated by a range of colors with blue indicating low expression and red indicating high expression. isomirs of mir-192-5p display more expanded or restricted expression with respect to tissues. additional file 13: figure s mir-3473f-pre isomirs in the testis. the isomirs of select testis enriched mirnas and their sequences are displayed in a heat map with tissues listed on the x-axis and isomirs listed on the y-axis. expression levels are indicated by a range of colors with blue indicating low expression and red indicating high expression. isomirs of mir-3473f-pre display more expanded or restricted expression with respect to tissues. additional file 14: figure s mir-3473f-pre isomirs in the testis. the isomirs of select testis enriched mirnas and their sequences are displayed in a heat map with tissues listed on the x-axis and isomirs listed on the y-axis. expression levels are indicated by a range of colors with blue indicating low expression and red indicating high expression. isomirs of mir-3473f-pre display more expanded or restricted expression with respect to tissues. additional file 15: table s p-values of rat caerulein toxicity study. statistical analysis of mirna and clinical chemistry biomarkers measured in the rat caerulein toxicity study are presented in this table with their associated fold changes and p-values indicated at 1 h, 4 h, 8 h, 24 h and 48 h. additional file 16: table s p-values of dog caerulein toxicity study. statistical analysis of mirna biomarkers measured in the dog caerulein toxicity study are displayed in this table with their associated p-values, lsmean, lower confidence intervals, upper confidence intervals and t ratios indicated at 1 h, 4 h, 8 h, 24 h and 48 h for the each dose group. additional file 17: figure s representation of rat mirna sequencing data. rat mirna sequencing data is represented with the normalized reads per million on the y-axis and the tissues on the x-axis. each color within the figure represents an individual mirna while the height of the colored bar represents the percentage of reads that a particular mirna makes up within that tissue. mirs-216a-5p, 375-3p and 148a-3p are indicated in the yellow, green and pink bars in the pancreas tissue respectively. additional file 18: figure s qpcr of pancreas enriched mirnas. total rna from the pancreas tissue of each individual rat was normalized by mass and examined by qpcr for expression of mirs-216a-5p, 375-3p and 148a-3p. ct values are indicated on the y-axis and the rat pancreas samples are indicated on the x-axis. each mirna appears to be expressed at approximately the same ct value indicating they are present at roughly the same quantity in the pancreas. abbreviations alpalkaline phosphatase astaspartate aminotransferase altalanine aminotransferase bunblood urea nitrogen cnscentral nervous system co2carbon dioxoide creatinineurinary albumin kim-1kidney injury molecule-1 ctnicardiac troponin i ggtgamma-glutamyl transferase gldhglutamate dehydrogenase hesihealth and environmental sciences institute hmgb1high mobility group protein b1 mirmicrorna ngalneutrophil gelatinase associated lipocalin qpcrquantitative polymerase chain reaction riscrna induced silencing complex rnaribonucleic acid sdsprague dawley tmpd, tetramethyl-p-phenylenediamine this research was performed in collaboration with hesi genomics committee 
we generated sequences that mapped to mature mirnas in rat, mouse and human from mirbase and discovered an additional rat mirnas as compared to the current number of rat mirnas in mirbase version tissue specific and enriched mirnas were identified and a subset of these mirnas were validated by qpcr for tissue specificity or enrichment. the studies demonstrate that conserved tissue specific/enriched mirs-216a-5p, 375-3p, 148a-3p, 216b-5p and 141-3p are candidate biomarkers of pancreatic injury in the rat and dog. microrna tissue profiling in rodents has been published, but sample an insufficient number of organs of toxicologic interest using microarray or qpcr technologies for mirna detection. here we impart an improved rat microrna body atlas consisting of and tissues of toxicologic interest from male and female sprague dawley rats respectively, using illumina mirna sequencing. as an example of the power of this approach, we have conducted rat and dog pancreas toxicity studies and examined the levels of some tissue specific and enriched mirnas conserved between rat and dog in the serum of each species. they can be detected at low concentrations in circulation in animals and humans, generating interest in the utilization of mirnas as serum and/or plasma based biomarkers of tissue injury. several of the authors created a dog mirna body atlas and we collaborated to test mirnas conserved in rat and dog pancreas in caerulein toxicity studies utilizing both species. micrornas are stable in the serum and plasma, are detectable in a wide variety of body fluids, are conserved across veterinary species and humans and are expressed in a tissue specific manner. a rich data set is presented that more robustly defines the tissue specificity and enrichment profiles of previously published and undiscovered rat mirnas. a microrna body atlas for rat and dog was useful in identifying new candidate mirna biomarkers of organ toxicity in toxicologically relevant species. micrornas are ~19– nucleotide long rna molecules that fine tune gene expression through the inhibition of translation or degradation of the mrna through incorporation into the rna induced silencing complex.
 micrornas are ~19– nucleotide long rna molecules that fine tune gene expression through the inhibition of translation or degradation of the mrna through incorporation into the rna induced silencing complex. micrornas are stable in the serum and plasma, are detectable in a wide variety of body fluids, are conserved across veterinary species and humans and are expressed in a tissue specific manner. they can be detected at low concentrations in circulation in animals and humans, generating interest in the utilization of mirnas as serum and/or plasma based biomarkers of tissue injury. microrna tissue profiling in rodents has been published, but sample an insufficient number of organs of toxicologic interest using microarray or qpcr technologies for mirna detection. here we impart an improved rat microrna body atlas consisting of and tissues of toxicologic interest from male and female sprague dawley rats respectively, using illumina mirna sequencing. several of the authors created a dog mirna body atlas and we collaborated to test mirnas conserved in rat and dog pancreas in caerulein toxicity studies utilizing both species.  a rich data set is presented that more robustly defines the tissue specificity and enrichment profiles of previously published and undiscovered rat mirnas. we generated sequences that mapped to mature mirnas in rat, mouse and human from mirbase and discovered an additional rat mirnas as compared to the current number of rat mirnas in mirbase version tissue specific and enriched mirnas were identified and a subset of these mirnas were validated by qpcr for tissue specificity or enrichment. as an example of the power of this approach, we have conducted rat and dog pancreas toxicity studies and examined the levels of some tissue specific and enriched mirnas conserved between rat and dog in the serum of each species. the studies demonstrate that conserved tissue specific/enriched mirs-216a-5p, 375-3p, 148a-3p, 216b-5p and 141-3p are candidate biomarkers of pancreatic injury in the rat and dog.  a microrna body atlas for rat and dog was useful in identifying new candidate mirna biomarkers of organ toxicity in toxicologically relevant species. electronic supplementary material the online version of this article contains supplementary material, which is available to authorized users. keywords micrornabiomarkersdeep sequencingtissue specifictissue enrichedpancreatic toxicityissue-copyright-statement© the author 2016 