An Improved microRNA Annotation of the Canine Genome

A review on microRNA detection and expression studies in dogs

MicroRNAs (miRNAs) are small non-coding RNAs that function by post-transcriptional regulation of gene expression. Their stability and abundance in tissue and body fluids makes them promising potential tools for both the diagnosis and prognosis of diseases and attractive therapeutic targets in humans and dogs. Studies of miRNA expression in normal and disease processes in dogs are scarce compared to studies published on miRNA expression in human disease. In this literature review, we identified 461 peer-reviewed papers from database searches using the terms "canine," "dog," "miRNA," and "microRNA"; we screened 244 for inclusion criteria and then included a total of 148 original research peer-reviewed publications relating to specific miRNA expression in canine samples. We found an overlap of miRNA expression changes between the four groups evaluated (normal processes, non-infectious and non-inflammatory conditions, infectious and/or inflammatory conditions, and neoplasia) in 39 miRNAs, 83 miRNAs in three of the four groups, 110 miRNAs in two of the three groups, where 158 miRNAs have only been reported in one of the groups. Additionally, the mechanism of action of these overlapping miRNAs varies depending on the disease process, elucidating a need for characterization of the mechanism of action of each miRNA in each disease process being evaluated. Herein we also draw attention to the lack of standardization of miRNA evaluation, consistency within a single evaluation method, and the need for standardized methods for a direct comparison.

Historical RNA expression profiles from the extinct Tasmanian tiger

Paleogenomics continues to yield valuable insights into the evolution, population dynamics, and ecology of our ancestors and other extinct species. However, DNA sequencing cannot reveal tissue-specific gene expression, cellular identity, or gene regulation, which are only attainable at the transcriptional level. Pioneering studies have shown that useful RNA can be extracted from ancient specimens preserved in permafrost and historical skins from extant canids, but no attempts have been made so far on extinct species. We extract, sequence, and analyze historical RNA from muscle and skin tissue of a ∼130-year-old Tasmanian tiger (Thylacinus cynocephalus) preserved in desiccation at room temperature in a museum collection. The transcriptional profiles closely resemble those of extant species, revealing specific anatomical features such as slow muscle fibers or blood infiltration. Metatranscriptomic analysis, RNA damage, tissue-specific RNA profiles, and expression hotspots genome-wide further confirm the thylacine origin of the sequences. RNA sequences are used to improve protein-coding and noncoding annotations, evidencing missing exonic loci and the location of ribosomal RNA genes while increasing the number of annotated thylacine microRNAs from 62 to 325. We discover a thylacine-specific microRNA isoform that could not have been confirmed without RNA evidence. Finally, we detect traces of RNA viruses, suggesting the possibility of profiling viral evolution. Our results represent the first successful attempt to obtain transcriptional profiles from an extinct animal species, providing thought-to-be-lost information on gene expression dynamics. These findings hold promising implications for the study of RNA molecules across the vast collections of natural history museums and from well-preserved permafrost remains.

MicroRNA profiling of cerebrospinal fluid from dogs with steroid responsive meningitis-arteritis and meningoencephalitis of unknown origin

Introduction

Non-infectious inflammatory diseases of the central nervous system in dogs, such as steroid responsive meningitis-arteritis (SRMA) and meningoencephalitis of unknown origin (MUO), represent a common clinical challenge that needs extensive and multimodal work-up to reach a presumptive diagnosis. Both diseases are presumably caused by dysregulations of the immune system, but further research is needed in order to understand the molecular mechanisms behind each disease and to optimize treatment.

Methods

By next-generation sequencing and subsequent quantitative real-time PCR (qPCR) verification, we designed a prospective case-control pilot study to analyze the small RNA profiles of cerebrospinal fluid from dogs suffering from MUO (N = 5), dogs suffering from SRMA (N = 8), and healthy dogs (N = 5) presented for elective euthanasia used as the Control group.

Results

Our results showed an overall enrichment in Y-RNA fragments across all samples, followed by microRNAs (miRNAs) and ribosomal RNAs as the major findings. Additional traces of short RNA reads mapped to long non-coding RNAs and protein-coding genes were also found. From the detected canine miRNAs, miR-21, miR-486, miR-148a, miR-99a, miR-191 and miR-92a were among the most abundant. Dogs with SRMA showed higher differences in miRNA abundance than dogs with MUO when compared to healthy dogs, and miR-142-3p was consistently detected as differentially upregulated in both diseases, although at a low concentration. Moreover, miR-405-5p and miR-503-5p showed different profiles between SRMA and MUO dogs. Subsequent qPCR analyses confirmed miR-142-5p, miR-191-5p and miR-92a-3p as significantly upregulated miRNAs in dogs with SRMA and/or MUO.

Discussion

Cerebrospinal fluid is a challenging biological material to use for profiling miRNAs due to the low content of circulating RNAs. Despite this, we could confirm several miRNAs being differentially abundant when comparing healthy dogs and dogs with MUO and SRMA, respectively. The results of this study indicate a potential role of miRNAs in the underlying molecular mechanisms of these diseases and establish the basis for further studies.

Urinary microRNAome in healthy cats and cats with pyelonephritis or other urological conditions

MicroRNAs (miRNAs) are short non-coding RNAs that regulate gene expression at the post-transcriptional level. miRNAs have been found in urine and have shown diagnostic potential in human nephropathies. Here, we aimed to characterize, for the first time, the feline urinary miRNAome and explore the use of urinary miRNA profiles as non-invasive biomarkers for feline pyelonephritis (PN). Thirty-eight cats were included in a prospective case-control study and classified in five groups: healthy Control cats (n = 11), cats with PN (n = 10), cats with subclinical bacteriuria or cystitis (SB/C, n = 5), cats with ureteral obstruction (n = 7) and cats with chronic kidney disease (n = 5). By small RNA sequencing we identified 212 miRNAs in cat urine, including annotated (n = 137) and putative novel (n = 75) miRNAs. The 15 most highly abundant urinary miRNAs accounted for nearly 71% of all detected miRNAs, most of which were previously identified in feline kidney. Ninety-nine differentially abundant (DA) miRNAs were identified when comparing Control cats to cats with urological conditions and 102 DA miRNAs when comparing PN to other urological conditions. Tissue clustering analysis revealed that the majority of urine samples clustered close to kidney, which confirm the likely cellular origin of the secreted urinary miRNAs. Relevant DA miRNAs were verified by quantitative real-time PCR (qPCR). Eighteen miRNAs discriminated Control cats from cats with a urological condition. Of those, seven miRNAs were DA by both RNAseq and qPCR methods between Control and PN cats (miR-125b-5p, miR-27a-3p, miR-21-5p, miR-27b-3p, miR-125a-5p, miR-17-5p and miR-23a-3p) or DA between Control and SB/C cats (miR-125b-5p). Six additional miRNAs (miR-30b-5p, miR-30c, miR-30e-5p, miR-27a-3p, miR-27b-39 and miR-222) relevant for discriminating PN from other urological conditions were identified by qPCR alone (n = 4) or by both methods (n = 2) (P<0.05). This panel of 13 miRNAs has potential as non-invasive urinary biomarkers for diagnostic of PN and other urological conditions in cats.

COMMD1 Exemplifies the Power of Inbred Dogs to Dissect Genetic Causes of Rare Copper-Related Disorders

Wilson's Disease is a rare autosomal recessive disorder in humans, often presenting with hepatic copper overload. Finding the genetic cause of a rare disease, especially if it is related to food constituents like the trace element copper, is a Herculean task. This review describes examples of how the unique population structure of in-bred dog strains led to the discovery of a novel gene and two modifier genes involved in inherited copper toxicosis. COMMD1, after the discovery in 2002, was shown to be a highly promiscuous protein involved in copper transport, protein trafficking/degradation, regulation of virus replication, and inflammation. Mutations in the ATP7A and ATP7B proteins in Labrador retrievers and Dobermann dogs resulted in a wide variation in hepatic copper levels in these breeds. To our knowledge, numerous dog breeds with inherited copper toxicosis of unknown genetic origin exist. Therefore, the possibility that men's best friend will provide new leads in rare copper storage diseases seems realistic.

A novel canine reference genome resolves genomic architecture and uncovers transcript complexity

We present GSD_1.0, a high-quality domestic dog reference genome with chromosome length scaffolds and contiguity increased 55-fold over CanFam3.1. Annotation with generated and existing long and short read RNA-seq, miRNA-seq and ATAC-seq, revealed that 32.1% of lifted over CanFam3.1 gaps harboured previously hidden functional elements, including promoters, genes and miRNAs in GSD_1.0. A catalogue of canine "dark" regions was made to facilitate mapping rescue. Alignment in these regions is difficult, but we demonstrate that they harbour trait-associated variation. Key genomic regions were completed, including the Dog Leucocyte Antigen (DLA), T Cell Receptor (TCR) and 366 COSMIC cancer genes. 10x linked-read sequencing of 27 dogs (19 breeds) uncovered 22.1 million SNPs, indels and larger structural variants. Subsequent intersection with protein coding genes showed that 1.4% of these could directly influence gene products, and so provide a source of normal or aberrant phenotypic modifications.

Whole genome analyses reveal significant convergence in obsessive-compulsive disorder between humans and dogs

Obsessive-compulsive disorder (OCD) represents a heterogeneous collection of diseases with diverse levels of phenotypic, genetic, and etiologic variability, making it difficult to identify the underlying genetic and biological mechanisms in humans. Domestic dogs exhibit several OCD-like behaviors. Using continuous circling as a representative phenotype for OCD, we screened two independent dog breeds, the Belgian Malinois and Kunming Dog and subsequently sequenced ten circling dogs and ten unaffected dogs for each breed. Using population differentiation analyses, we identified 11 candidate genes in the extreme tail of the differentiated regions between cases and controls. These genes overlap significantly with genes identified in a genome wide association study (GWAS) of human OCD, indicating strong convergence between humans and dogs. Through gene expressional analysis and functional exploration, we found that two candidate OCD risk genes, PPP2R2B and ADAMTSL3, affected the density and morphology of dendritic spines. Therefore, changes in dendritic spine may underlie some common biological and physiological pathways shared between humans and dogs. Our study revealed an unprecedented level of convergence in OCD shared between humans and dogs, and highlighted the importance of using domestic dogs as a model species for many human diseases including OCD.

Ancient microRNA profiles of a 14,300-year-old canid samples confirm taxonomic origin and give glimpses into tissue-specific gene regulation from the Pleistocene

DNA sequencing is the current key technology for historic or ancient biological samples and has led to many exciting discoveries in the field of paleogenomics. However, functional insights into tissue identity, cellular composition or gene regulation cannot be gained from DNA. Recent analyses have shown that, under favorable conditions, RNA can also be sequenced from ancient samples, enabling studies at the transcriptomic and regulatory level. Analyzing ancient RNA data from a Pleistocene canid, we find hundreds of intact microRNAs that are taxonomically informative, show tissue-specificity and have functionally predictive characteristics. With an extraordinary age of 14,300 years, these microRNA sequences are by far the oldest ever reported. The authenticity of the sequences is further supported by a) the presence of canid / Caniformia-specific sequences that never evolved outside of this clade, b) tissue-specific expression patterns (cartilage, liver and muscle) that resemble those of modern dogs and c) RNA damage patterns that are clearly distinct from those of fresh samples. By performing computational microRNA-target enrichment analyses on the ancient sequences, we predict microRNA functions consistent with their tissue pattern of expression. For instance, we find a liver-specific microRNA that regulates carbohydrate metabolism and starvation responses in canids. In summary, we show that straightforward paleotranscriptomic microRNA analyses can give functional glimpses into tissue identity, cellular composition and gene regulatory activity of ancient samples and biological processes that took place in the Pleistocene, thus holding great promise for deeper insights into gene regulation in extinct animals based on ancient RNA sequencing. .

The molecular effect of a polymorphic microRNA binding site of Wolfram syndrome 1 gene in dogs

Background

Although the molecular function of wolframin remains unclear, the lack of this protein is known to cause stress in the endoplasmic reticulum. Some variants in the Wolfram Syndrome 1 gene (WFS1) were associated with various neuropsychiatric disorders in humans, such as aggressiveness, impulsivity and anxiety.

Results

Here we present an in silico study predicting a single nucleotide polymorphism (rs852850348) in the canine WFS1 gene which was verified by direct sequencing and was genotyped by a PCR-based technique. We found that the rs852850348 polymorphism is located in a putative microRNA (cfa-miR-8834a and cfa-miR-1838) binding site. Therefore, the molecular effect of allelic variants was studied in a luciferase reporter system that allowed assessing gene expression. We demonstrated that the variant reduced the activity of the reporter protein expression in an allele-specific manner. Additionally, we performed a behavioral experiment and investigated the association with this locus to different performance in this test. Association was found between food possessivity and the studied WFS1 gene polymorphism in the Border collie breed.

Conclusions

Based on our findings, the rs852850348 locus might contribute to the genetic risk of possessivity behavior of dogs in at least one breed and might influence the regulation of wolframin expression.

Genetic Pathways of Aging and Their Relevance in the Dog as a Natural Model of Human Aging

Aging research has experienced a burst of scientific efforts in the last decades as the growing ratio of elderly people has begun to pose an increased burden on the healthcare and pension systems of developed countries. Although many breakthroughs have been reported in understanding the cellular mechanisms of aging, the intrinsic and extrinsic factors that contribute to senescence on higher biological levels are still barely understood. The dog, Canis familiaris, has already served as a valuable model of human physiology and disease. The possible role the dog could play in aging research is still an open question, although utilization of dogs may hold great promises as they naturally develop age-related cognitive decline, with behavioral and histological characteristics very similar to those of humans. In this regard, family dogs may possess unmatched potentials as models for investigations on the complex interactions between environmental, behavioral, and genetic factors that determine the course of aging. In this review, we summarize the known genetic pathways in aging and their relevance in dogs, putting emphasis on the yet barely described nature of certain aging pathways in canines. Reasons for highlighting the dog as a future aging and gerontology model are also discussed, ranging from its unique evolutionary path shared with humans, its social skills, and the fact that family dogs live together with their owners, and are being exposed to the same environmental effects.

Brain-enriched microRNAs circulating in plasma as novel biomarkers for Rett syndrome

Rett syndrome (RTT) is a severe neurodevelopmental disorder caused by mutations in the X-linked gene MECP2 (methyl-CpG-binding protein 2). Minimally invasive and accurate biomarkers of disease progression and treatment response could facilitate screening of therapeutic compounds in animal models, enrollment of better-defined participants into clinical trials, and treatment monitoring. In this study, we used a targeted approach based on analysis of brain-enriched microRNAs (miRNAs) circulating in plasma to identify miRNA biomarkers of RTT using Mecp2-mutant mice as a model system and human plasma samples. An “miRNA pair” approach, i.e. the ratio between two miRNAs, was used for data normalization. Specific miRNA pairs and their combinations (classifiers) analyzed in plasma differentiated wild-type from Mecp2 male and female mice with >90% accuracy. Individual miRNA pairs were more effective in distinguishing male (homozygous) animals than female (heterozygous) animals, suggesting that disease severity correlated with the levels of the miRNA biomarkers. In the human study, 30 RTT patients were compared with age-matched controls. The results of this study showed that miRNA classifiers were able to differentiate RTT patients from controls with 85–100% sensitivity. In addition, a comparison of various age groups demonstrated that the dynamics in levels of miRNAs appear to be associated with disease development (involvement of liver, muscle and lipid metabolism in the pathology). Importantly, certain miRNA biomarker pairs were common to both the animal models and human subjects, indicating the similarity between the underlying pathological processes. The data generated in this feasibility study suggest that circulating miRNAs have the potential to be developed as markers of RTT progression and treatment response. Larger clinical studies are needed to further evaluate the findings presented here.

Genome-wide microRNA screening in Nile tilapia reveals pervasive isomiRs' transcription, sex-biased arm switching and increasing complexity of expression throughout development

MicroRNAs (miRNAs) are key regulators of gene expression in multicellular organisms. The elucidation of miRNA function and evolution depends on the identification and characterization of miRNA repertoire of strategic organisms, as the fast-evolving cichlid fishes. Using RNA-seq and comparative genomics we carried out an in-depth report of miRNAs in Nile tilapia (Oreochromis niloticus), an emergent model organism to investigate evo-devo mechanisms. Five hundred known miRNAs and almost one hundred putative novel vertebrate miRNAs have been identified, many of which seem to be teleost-specific, cichlid-specific or tilapia-specific. Abundant miRNA isoforms (isomiRs) were identified with modifications in both 5p and 3p miRNA transcripts. Changes in arm usage (arm switching) of nine miRNAs were detected in early development, adult stage and even between male and female samples. We found an increasing complexity of miRNA expression during ontogenetic development, revealing a remarkable synchronism between the rate of new miRNAs recruitment and morphological changes. Overall, our results enlarge vertebrate miRNA collection and reveal a notable differential ratio of miRNA arms and isoforms influenced by sex and developmental life stage, providing a better picture of the evolutionary and spatiotemporal dynamics of miRNAs.

Identification of microRNAs associated with glioma diagnosis and prognosis

The sensitivity and specificity of microRNAs (miRNAs) for diagnosing glioma are controversial. We therefore performed a meta-analysis to systematically identify glioma-associated miRNAs. We initially screened five miRNA microarray datasets to evaluate the differential expression of miRNAs between glioma and normal tissues. We next compared the expression of the miRNAs in different organs and tissues to assess the sensitivity and specificity of the differentially expressed miRNAs in the diagnosis of glioma. Finally, pathway analysis was performed using GeneGO. We identified 27 candidate miRNAs associated with glioma initiation, progression, and patient prognosis. Sensitivity and specificity analysis indicated miR-15a, miR-16, miR-21, miR-23a, and miR-9 were up-regulated, while miR-124 was down-regulated in glioma. Ten signaling pathways showed the strongest association with glioma development and progression: the p53 pathway feedback loops 2, Interleukin signaling pathway, Toll receptor signaling pathway, Parkinson's disease, Notch signaling pathway, Cadherin signaling pathway, Apoptosis signaling pathway, VEGF signaling pathway, Alzheimer disease-amyloid secretase pathway, and the FGF signaling pathway. Our results indicate that the integration of miRNA, gene, and protein expression data can yield valuable biomarkers for glioma diagnosis and treatment. Indeed, six of the miRNAs identified in this study may be useful diagnostic and prognostic biomarkers in glioma.

Circulating brain-enriched microRNAs as novel biomarkers for detection and differentiation of neurodegenerative diseases

Background

Minimally invasive specific biomarkers of neurodegenerative diseases (NDs) would facilitate patient selection and disease progression monitoring. We describe the assessment of circulating brain-enriched microRNAs as potential biomarkers for Alzheimer’s disease (AD), frontotemporal dementia (FTD), Parkinson’s disease (PD), and amyotrophic lateral sclerosis (ALS).

Methods

In this case-control study, the plasma samples were collected from 250 research participants with a clinical diagnosis of AD, FTD, PD, and ALS, as well as from age- and sex-matched control subjects (n = 50 for each group), recruited from 2003 to 2015 at the University of Pennsylvania Health System, including the Alzheimer’s Disease Center, the Parkinson’s Disease and Movement Disorders Center, the Frontotemporal Degeneration Center, and the Amyotrophic Lateral Sclerosis Clinic. Each group was randomly divided into training and confirmation sets of equal size. To evaluate the potential of circulating microRNAs enriched in specific brain regions affected by NDs and present in synapses as biomarkers of NDs, the levels of 37 brain-enriched and inflammation-associated microRNAs in the plasma of all participants were measured using individual qRT-PCR. A “microRNA pair” approach was used for data normalization.

Results

MicroRNA pairs and their combinations (classifiers) capable of differentiating NDs from control and from each other were defined using independently and jointly analyzed training and confirmation datasets. AD, PD, FTD, and ALS are differentiated from control with accuracy of 0.89, 0.90, 0.88, and 0.83 (AUCs, 0.96, 0.96, 0.94, and 0.93), respectively; NDs are differentiated from each other with accuracy ranging from 0.77 (AUC, 0.87) for AD vs. FTD to 0.93 (AUC, 0.98) for AD vs. ALS. The data further indicate sex dependence of some microRNA markers. The average increase in accuracy in distinguishing ND from control for all and male/female groups is 0.06; the largest increase is for ALS, from 0.83 for all participants to 0.92/0.98 for male/female participants.

Conclusions

The work presented here suggests the possibility of developing microRNA-based diagnostics for detection and differentiation of NDs. Larger multicenter clinical studies are needed to further evaluate circulating brain-enriched microRNAs as biomarkers for NDs and to investigate their association with other ND biomarkers in clinical trial settings.

Electronic supplementary material

The online version of this article (doi:10.1186/s13195-017-0316-0) contains supplementary material, which is available to authorized users.

Characterization of a low expression haplotype in canine glutathione S-transferase (GSTT1) and its prevalence in golden retrievers

Glutathione S-transferase-theta (GSTT1) is a carcinogen detoxification enzyme, and low activity variants are associated with lymphoma in humans. We recently found a variant in the 3' untranslated region (UTR) of canine GSTT1, *101_102insT, which was predicted to change miRNA binding and was found in 5 of 17 golden retriever (GR) dogs with lymphoma but none of 14 healthy GRs. The aim of this study was to determine whether this variant led to decreased GSTT1 expression and was a discernible risk factor for lymphoma within the GR breed. On resequencing, *101_102insT appeared to be in complete linkage disequilibrium with 3 additional 3'UTR variants, leading to the inferred haplotype *3T>C; *101_102insT; *190C>A; *203T>C. In canine livers that were heterozygous for this variant haplotype, GSTT1 protein expression was significantly lower compared to the reference haplotype (densitometry .40 vs .64, P = .022), and GSTT1 transcript levels by qPCR were also significantly lower (fold difference .52, P = .012), without evidence of substantial allelic expression imbalance. The variant haplotype led to >50% decrease in expression in vitro (.31 ± .07 vs .64 ± .19; P = .019). We found no significant difference in minor allele frequencies between 71 GR dogs with lymphoma (MAF .162) and 33 healthy age-matched controls (MAF .136, P = .69). Our results indicate that the variant GSTT1 3'UTR haplotype containing *101_102insT reduces gene expression, which could lead to impaired carcinogen detoxification, but was not a detectable risk factor for lymphoma in GR dogs.

Restriction of Aerobic Metabolism by Acquired or Innate Arylsulfatase B Deficiency: A New Approach to the Warburg Effect

Aerobic respiration is required for optimal efficiency of metabolism in mammalian cells. Under circumstances when oxygen utilization is impaired, cells survive by anerobic metabolism. The malignant cell has cultivated the use of anerobic metabolism in an aerobic environment, the Warburg effect, but the explanation for this preference is not clear. This paper presents evidence that deficiency of the enzyme arylsulfatase B (ARSB; N-acetylgalactosamine 4-sulfatase), either innate or acquired, helps to explain the Warburg phenomenon. ARSB is the enzyme that removes 4-sulfate groups from the non-reducing end of chondroitin 4-sulfate and dermatan sulfate. Previous reports indicated reduced ARSB activity in malignancy and replication of the effects of hypoxia by decline in ARSB. Hypoxia reduced ARSB activity, since molecular oxygen is needed for post-translational modification of ARSB. In this report, studies were performed in human HepG2 cells and in hepatocytes from ARSB-deficient and normal C57BL/6J control mice. Decline of ARSB, in the presence of oxygen, profoundly reduced the oxygen consumption rate and increased the extracellular acidification rate, indicating preference for aerobic glycolysis. Specific study findings indicate that decline in ARSB activity enhanced aerobic glycolysis and impaired normal redox processes, consistent with a critical role of ARSB and sulfate reduction in mammalian metabolism.