In silico detection and characteristics of novel microRNA genes in the Equus caballus genome using an integrated ab initio and comparative genomic approach

Effects of corn supplementation on serum and muscle microRNA profiles in horses

Laminitis associated with equine metabolic syndrome causes significant economic losses in the equine industry. Diets high in non-structural carbohydrates (NSC) have been linked to insulin resistance and laminitis in horses. Nutrigenomic studies analyzing the interaction of diets high in NSCs and gene expression regulating endogenous microRNAs (miRNA) are rare. This study's objectives were to determine whether miRNAs from dietary corn can be detected in equine serum and muscle and its impacts on endogenous miRNA. Twelve mares were blocked by age, body condition score, and weight and assigned to a control (mixed legume grass hay diet) and a mixed legume hay diet supplemented with corn. Muscle biopsies and serum were collected on Days 0 and 28. Transcript abundances were analyzed using qRT-PCR for three plant-specific and 277 endogenous equine miRNAs. Plant miRNAs were found in serum and skeletal muscle samples with a treatment effect (p < .05) with corn-specific miRNA being higher than control in serum after feeding. Endogenous miRNAs showed 12 different (p < .05) miRNAs in equine serum after corn supplementation, six (eca-mir16, -4863p, -4865p, -126-3p, -296, and -192) previously linked to obesity or metabolic disease. The results of our study indicate that dietary plant miRNAs can appear in circulation and tissues and may regulate endogenous genes.

Assessment and Distribution of Runs of Homozygosity in Horse Breeds Representing Different Utility Types

The present study reports runs of homozygosity (ROH) distribution in the genomes of six horse breeds (571 horses in total) representing three horse types (primitive, light, and draft horses) based on the 65k Equine BeadChip assay. Of major interest was the length, quantity, and frequency of ROH characteristics, as well as differences between horse breeds and types. Noticeable differences in the number, length and distribution of ROH between breeds were observed, as well as in genomic inbreeding coefficients. We also identified regions of the genome characterized by high ROH coverage, known as ROH islands, which may be signals of recent selection events. Eight to fourteen ROH islands were identified per breed, which spanned multiple genes. Many were involved in important horse breed characteristics, including WFIKNN2, CACNA1G, STXBP4, NOG, FAM184B, QDPR, LCORL, and the zinc finger protein family. Regions of the genome with zero ROH occurrences were also of major interest in specific populations. Depending on the breed, we detected between 2 to 57 no-ROH regions and identified 27 genes in these regions that were common for five breeds. These genes were involved in, e.g., muscle contractility (CACNA1A) and muscle development (miR-23, miR-24, miR-27). To sum up, the obtained results can be furthered analyzed in the topic of identification of markers unique for specific horse breed characteristics.

Comparative Analysis for Genetic Characterization in Korean Native Jeju Horse

Simple Summary

In modern times, horse breeds, mostly in horse racing, are the Thoroughbred varieties obtained by breeding three Godolphin Arabians with British mares in England. Especially in Jeju Island, Korea, Jeju horses have been introduced from Mongolia since the 13th century. They have contributed a lot to the agricultural community, but their population has been rapidly decreasing due to rapid agricultural industrialization. Therefore, we sympathize with Jeju horse-specific genetic variation and compare and analyze evolutionary correlations by utilizing Whole Genome Sequencing analysis to evaluate the genetic diversity of Jeju horses and preserve genetic information. We explored Jeju horse-specific genetic differences through a comparative analysis of large-capacity genomic data between the public database and a Thoroughbred variety. In order to adapt to the barren external environment, it is predicted that Jeju horses have experienced strong positive selection in the direction of accumulating many genetic variations, enough to cause functional differences in the eqCD1a6 gene to have an efficient immune function. In addition, we further validate the Jeju horse-specific single nucleotide polymorphisms in the aqCD1a6 gene by employing the digital PCR method, a diagnostic technique for genetic variations.

Abstract

The Jeju horse is a native Korean species that has been breeding on Jeju Island since the 13th century. Their shape has a distinct appearance from the representative species, Thoroughbred. Here, we performed a comparison of the Jeju horse and Thoroughbred horse for the identification of genome-wide structure variation by using the next-generation sequencing (NGS) technique. We generated an average of 95.59 Gb of the DNA sequence, resulting in an average of 33.74 X sequence coverage from five Jeju horses. In addition, reads obtained from WGRS data almost covered the horse reference genome (mapped reads 98.4%). Based on our results, we identified 1,244,064 single nucleotide polymorphisms (SNPs), 113,498 genomic insertions, and 114,751 deletions through bioinformatics analysis. Interestingly, the results of the WGRS comparison indicated that the eqCD1a6 gene contains signatures of positive natural selection in Jeju horses. The eqCD1a6 gene is known to be involved in immunity. The eqCD1a6 gene of Jeju horses commonly contained 296 variants (275 SNPs and 21 INDELs) that were compared with its counterpart of two Thoroughbred horses. In addition, we used LOAA, digital PCR, to confirm the possibility of developing a molecular marker for species identification using variant sites. As a result, it was possible to confirm the result of the molecular marker with high accuracy. Nevertheless, eqCD1a6 was shown to be functionally intact. Taken together, we have found significant genomic variation in these two different horse species.

Characterization of lncRNA-Perturbed TLR-Signaling Network Identifies Novel lncRNA Prognostic Biomarkers in Colorectal Cancer

Increasing evidence has suggested that long non-coding RNAs (lncRNAs) are critical regulators in the Toll-like receptors (TLR)-signaling network to modulate colorectal cancer (CRC) development and progression. However, the mechanism and clinical significance for lncRNAs regulating TLR signaling pathways in CRC remained largely unknown. In this study, we performed an integrative network analysis of transcriptomics by focusing on a lncRNA-perturbed TLR-signaling network, identifying 280 lncRNAs and 122 mRNAs. We found a profound phenomenon that abnormal expression of some lncRNAs can perturb the TLR-signaling network to contribute to CRC development and progression. Furthermore, we identified a novel TLR-related prognostic gene signature (TLRLncSig) composed of three lncRNAs (MCHR2, AC011472.4, and AC063944.1), and one mRNA (CDKN2B). Utilizing TLRLncSig could classify CRC patients of training set into two groups with significantly different overall survival. The prognostic value of the TLRLncSig was further validated in the other two independent CRC datasets with different platforms. Results of multivariate and stratification analysis indicated that the TLRLncSig is an independent prognostic factor, and our study underscores the clinical significance of TLR-related lncRNAs in CRC development and progression.

Overexpression of miR-20a promotes the progression of osteosarcoma by directly targeting QKI2

Osteosarcoma (OS) is the most common type of malignant primary bone neoplasm. Although the application of neoadjuvant chemotherapy has improved the 5-year survival rate of patients suffering from OS, prognosis remains poor. Therefore, it is important to elucidate the molecular mechanisms underlying the occurrence, progression and metastasis of OS. The RNA-binding protein Quaking (QKI) is a member of the STAR family of proteins, and can function as a tumor suppressor gene to suppress the occurrence and progression of a variety of tumors; however, the role of QKI in OS remains to be fully elucidated. In the present study, it was identified that the expression of QKI2 was downregulated in OS using western blot analysis. In addition, subsequent functional investigations, including MTT, Transwell invasion and migration assays, revealed that QKI2 inhibited the proliferation, invasion and migration of an OS cell line in vitro. By implementing a series of experimental techniques in molecular biology, including reverse transcription-quantitative polymerase chain reaction and a double fluorescence reporter assay, it was demonstrated that the expression of miR-20a was high and inhibited the expression of QKI2 in OS. In conclusion, it was revealed that aberrantly upregulated miR-20a inhibited the expression of QKI2 in OS by targeting QKI2 mRNA, subsequently promoting the proliferation, migration and invasion of OS cells.

Cloning and characterization of Bubaline mammary miRNAs: An in silico approach

MicroRNAs (miRNAs) are small non-coding RNAs ~ 19-25 nucleotides long that are involved in the regulation of gene expression. They negatively regulate the gene expression via inhibition or complete degradation of mRNAs by binding the complementary target sequences in 3' untranslated region. The present investigation was aimed at profiling of miRNAs expressed in the Bubaline mammary tissue at dry stage of lactation cycle. Small RNAs were isolated from freshly collected mammary tissues and T4 RNA ligase was used to ligate the enriched miRNAs with 3' and 5' linker sequences in two separate reactions. cDNA copies were synthesized from linkered small RNAs followed by the PCR amplification. The PCR products were resolved on 15% non-denaturing polyacrylamide gel electrophoresis by gelstar staining. The PCR products were cloned using pGEM®-T easy vector system and the desired clones (with linkered small RNA sequences) were confirmed using restriction digestion of plasmids with EcoRI. Out of 15 Bubaline small RNA sequences, eight sequences (Seq. ID I-VIII) matched the size range of miRNA molecules i.e., 18-26 nucleotides. The Bubaline small RNA sequences II and III showed partial alignment with various mammalian and non-mammalian miRNAs. The small RNA sequences obtained in the present study did not show any perfect match with already reported mRNA, rRNA or tRNA sequences in different databases. Hence, only the Bubaline small RNA sequences that showed partial homology with miRNAs were considered as putative Bubaline miRNAs. The present study established the basic repertoire of miRNAs expressed at dry stage of lactation in Bubaline mammary gland.

MiR-126-3p promotes the cell proliferation and inhibits the cell apoptosis by targeting TSC1 in the porcine granulosa cells

In mammalian ovaries, many studies demonstrated that the proliferation and apoptosis of granulosa cells are involved in folliculogenesis. Previous evidence suggests that miR-126-3p might get involved in the proliferation and apoptosis of granulosa cells, and tuberous sclerosis complex 1 (TSC1) gene was predicted as one target of miR-126-3p, and moreover, granulosa cell-specific TSC1 knockout stimulated folliculogenesis in mice. However, the molecular regulation of miR-126-3p on TSC1 and its effects on cell proliferation and apoptosis remain virtually unexplored in granulosa cells. Using porcine granulosa cells as a model, the luciferase report assay, mutation, deletion, Annexin-V/PI staining, and EdU assays were applied to investigate the molecular mechanism for miR-126-3p regulating the expression of TSC1 and their effects on the cell proliferation and apoptosis. We found that miR-126-3p showed a positive effect on cell proliferation and a negative effect on cell apoptosis in porcine granulosa cells, and knockdown of TSC1 significantly promoted cell proliferation and significantly inhibited cell apoptosis in porcine granulosa cells. Furthermore, miR-126-3p might target and repress the expressions of TSC1 at the post-transcriptional level, thereby promoting cell proliferation and inhibiting cell apoptosis of granulosa cells. These findings would provide of great insight in further exploring the molecular regulation of miR-126-3p and TSC1 on the functions of granulosa cells during the folliculogenesis in mammals.

Small RNA-Seq reveals novel miRNAs shaping the transcriptomic identity of rat brain structures

Small RNA-Seq of the rat central nervous system reveals known and novel miRNAs specifically regulated in brain structures and correlated with the expression of their predicted target genes, suggesting a critical role in the transcriptomic identity of brain structures.

In the central nervous system (CNS), miRNAs are involved in key functions, such as neurogenesis and synaptic plasticity. Moreover, they are essential to define specific transcriptomes in tissues and cells. However, few studies were performed to determine the miRNome of the different structures of the rat CNS, although a major model in neuroscience. Here, we determined by small RNA-Seq, the miRNome of the olfactory bulb, the hippocampus, the cortex, the striatum, and the spinal cord and showed the expression of 365 known miRNAs and 90 novel miRNAs. Differential expression analysis showed that several miRNAs were specifically enriched/depleted in these CNS structures. Transcriptome analysis by mRNA-Seq and correlation based on miRNA target predictions suggest that the specifically enriched/depleted miRNAs have a strong impact on the transcriptomic identity of the CNS structures. Altogether, these results suggest the critical role played by these enriched/depleted miRNAs, in particular the novel miRNAs, in the functional identities of CNS structures.

MicroRNAs in equine veterinary science

MicroRNAs are small noncoding RNAs that play a pivotal role in diverse cellular processes through post-transcriptional regulation of gene expression. The dysregulation of specific microRNAs is associated with disease development and progression. In this review, we summarise how microRNAs modulate gene expression, and explain microRNA nomenclature. We discuss the potential applications of microRNAs in equine disease diagnosis and treatment, in the context of the sum of current knowledge about microRNA expression in normal and diseased equine tissues.

Circulating microRNA profiles of Hendra virus infection in horses

Hendra virus (HeV) is an emerging zoonotic pathogen harbored by Australian mainland flying foxes. HeV infection can cause lethal disease in humans and horses, and to date all cases of human HeV disease have resulted from contact with infected horses. Currently, diagnosis of acute HeV infections in horses relies on the productive phase of infection when virus shedding may occur. An assay that identifies infected horses during the preclinical phase of infection would reduce the risk of zoonotic viral transmission during management of HeV outbreaks. Having previously shown that the host microRNA (miR)-146a is upregulated in the blood of HeV-infected horses days prior to the detection of viremia, we have profiled miRNAs at the transcriptome-wide level to comprehensively assess differences between infected and uninfected horses. Next-generation sequencing and the miRDeep2 algorithm identified 742 mature miRNA transcripts corresponding to 593 miRNAs in whole blood of six horses (three HeV-infected, three uninfected). Thirty seven miRNAs were differentially expressed in infected horses, two of which were validated by qRT-PCR. This study describes a methodology for the transcriptome-wide profiling of miRNAs in whole blood and supports the notion that measuring host miRNA expression levels may aid infectious disease diagnosis in the future.

Characterization of basal and lipopolysaccharide-induced microRNA expression in equine peripheral blood mononuclear cells using Next-Generation Sequencing

The innate immune response to lipopolysaccharide contributes substantially to the morbidity and mortality of gram-negative sepsis. Horses and humans share an exquisite sensitivity to lipopolysaccharide and thus the horse may provide valuable comparative insights into this aspect of the inflammatory response. MicroRNAs, small non-coding RNA molecules acting as post-transcriptional regulators of gene expression, have key roles in toll-like receptor signaling regulation but have not been studied in this context in horses. The central hypothesis of this study was that lipopolysaccharide induces differential microRNA expression in equine peripheral blood mononuclear cells in a manner comparable to humans. Illumina Next Generation Sequencing was used to characterize the basal microRNA transcriptome in isolated peripheral blood mononuclear cells from healthy adult horses, and to evaluate LPS-induced changes in microRNA expression in cells cultured for up to four hours. Selected expression changes were validated using quantitative reverse-transcriptase PCR. Only miR-155 was significantly upregulated by LPS, changing in parallel with supernatant tumor necrosis factor-α concentration. Eight additional microRNAs, including miR-146a and miR-146b, showed significant expression change with time in culture without a clear LPS effect. Target predictions indicated a number of potential immunity-associated targets for miR-155 in the horse, including SOCS1, TAB2 and elements of the PI3K signaling pathway, suggesting that it is likely to influence the acute inflammatory response to LPS. Gene alignment showed extensive conservation of the miR-155 precursor gene and associated promoter regions between horses and humans. The basal and LPS-stimulated microRNA expression pattern characterized here were similar to those described in human leukocytes. As well as providing a resource for further research into the roles of microRNAs in immune responses in horses, this will facilitate inter-species comparative study of the role of microRNAs in the inflammatory cascade during endotoxemia and sepsis.

Identification and Expression of Equine MER-Derived miRNAs

MicroRNAs (miRNAs) are single-stranded, small RNAs (21-23 nucleotides) that function in gene silencing and translational inhibition via the RNA interference mechanism. Most miRNAs originate from host genomic regions, such as intergenic regions, introns, exons, and transposable elements (TEs). Here, we focused on the palindromic structure of medium reiteration frequencies (MERs), which are similar to precursor miRNAs. Five MER consensus sequences (MER5A1, MER53, MER81, MER91C, and MER117) were matched with paralogous transcripts predicted to be precursor miRNAs in the horse genome (equCab2) and located in either intergenic regions or introns. The MER5A1, MER53, and MER91C sequences obtained from RepeatMasker were matched with the eca-miR-544b, eca-miR-1302, and eca-miR-652 precursor sequences derived from Ensembl transcript database, respectively. Each precursor form was anticipated to yield two mature forms, and we confirmed miRNA expression in six different tissues (cerebrum, cerebellum, lung, spleen, adrenal gland, and duodenum) of one thorough-bred horse. MER5A1-derived miRNAs generally showed significantly higher expression in the lung than in other tissues. MER91C-derived miRNA-5p also showed significantly higher expression in the duodenum than in other tissues (cerebellum, lung, spleen, and adrenal gland). The MER117-overlapped expressed sequence tag generated polycistronic miRNAs, which showed higher expression in the duodenum than other tissues. These data indicate that horse MER transposons encode miR-NAs that are expressed in several tissues and are thought to have biological functions.

Computational identification and characterization of novel microRNA in the mammary gland of dairy goat (Capra hircus)

Many studies have indicated that microRNAs (miRNAs) influence the development of the mammary gland by posttranscriptionally affecting their target genes. The objective of this research was to identify novel miRNAs in the mammary gland of dairy goats with a bioinformatics approach that was based on expressed sequence tag (EST) and genome survey sequence (GSS) analyses. We applied all known major mammals, miRNAs to search against the goat EST and GSS databases for the first time to identify new miRNAs. We, then, validated these newly predicted miRNAs with stem-loop reverse transcription followed by a SYBR Green polymerase chain reaction assay. Finally, 29 mature miRNAs were identified and verified, and of these, 14 were grouped into 13 families based on seed sequence identity and 85 potential target genes of newly verified miRNAs were subsequently predicted, most of which seemed to encode the proteins participating in regulation of metabolism, signal transduction, growth and development. The predicting accuracy of the new miRNAs was 70.37%, which confirmed that the methods used in this study were efficient and reliable. Detailed analyses of the sequence characteristics of the novel miRNAs of the goat mammary gland were performed. In conclusion, these results provide a reference for further identification of miRNAs in animals without a complete genome and thus improve the understanding of miRNAs in the caprine mammary gland.

Genome-wide identification and characterization of microRNA genes and their targets in large yellow croaker (Larimichthys crocea)

MicroRNAs (miRNAs or miRs) are a class of non-coding RNAs of 20-25 nucleotides (nt) in length, which regulates the expression of gene in eukaryotic organism. Studies has been confirmed that miRNA plays an important role in various biological and metabolic processes in both animals and plants. Predicting new miRNAs by computer based homology search analysis is an effective way to discover novel miRNAs. Though a large number of miRNAs have been reported in many fish species, reports of miRNAs in large yellow croaker (L. crocea) are limited especially via the computational-based approaches. In this paper, a method of comparative genomic approach by computational genomic homology based on the conservation of miRNA sequences and the stem-loop hairpin secondary structures of miRNAs was adopted. A total of 199 potential miRNAs were predicted representing 81 families. 12 of them were chose to be validated by real time RT-PCR, apart from miR-7132b-5p which was not detected. Results indicated that the prediction method that we used to identify the miRNAs was effective. Furthermore, 948 potential target genes were predicted. Gene ontology (GO) analysis revealed that 175, 287, and 486 target genes were involved in cellular components, biological processes and molecular functions, respectively. Overall, our findings provide a first computational identification and characterization of L. crocea miRNAs and their potential targets in functional analysis, and will be useful in laying the foundation for further characterization of their role in the regulation of diversity of physiological processes.

RNA sequencing as a powerful tool in searching for genes influencing health and performance traits of horses

RNA sequencing (RNA-seq) by next-generation technology is a powerful tool which creates new possibilities in whole-transcriptome analysis. In recent years, with the use of the RNA-seq method, several studies expanded transcriptional gene profiles to understand interactions between genotype and phenotype, supremely contributing to the field of equine biology. To date, in horses, massive parallel sequencing of cDNA has been successfully used to identify and quantify mRNA levels in several normal tissues, as well as to annotate genes. Moreover, the RNA-seq method has been applied to identify the genetic basis of several diseases or to investigate organism adaptation processes to the training conditions. The use of the RNA-seq approach has also confirmed that horses can be useful as a large animal model for human disease, especially in the field of immune response. The presented review summarizes the achievements of profiling gene expression in horses (Equus caballus).

The role of microRNAs in equine medicine: a review

The search for new markers of diseases in human as well as veterinary medicine is ongoing. Recently, microRNAs (miRNAs or miRs) have emerged as potential new biomarkers. MiRNAs are short sequences of RNA (∼22 nucleotides) that regulate gene expression via their target messenger RNA (mRNA). Circulating miRNAs in blood can be used as novel diagnostic markers for diseases due to their evolutionary conservation and stability. As a consequence of their systemic and manifold effects on the gene expression in various target organs, the concept that miRNAs could function as hormones has been suggested. This review summarizes the biogenesis, maturation, and stability of miRNAs and discusses their use as potential biomarkers in equine medicine. To date, over 700 equine miRNAs are identified with distinct subsets of miRNAs differentially expressed in a tissue-specific manner. A physiological involvement of various miRNAs in the regulation of cell survival, steroidogenesis, and differentiation during follicle selection and ovulation in the monovular equine ovary has been demonstrated. Furthermore, miRNAs might be used as novel diagnostic markers for myopathies such as polysaccharide storage myopathy and recurrent exertional rhabdomyolysis as well as osteochondrosis. Preliminary data indicate that miRNAs in blood might play important roles in equine glucose metabolism pathway. Of note, breed differences have been reported regarding the normal equine miRNA signature. For disease prevention, it is of utmost importance to identify disease-associated biomarkers which help detect diseases before symptoms appear. As such, circulating miRNAs represent promising novel diagnostic markers in equine medicine.

Characterisation of novel microRNAs in the Black flying fox (Pteropus alecto) by deep sequencing

Background

Bats are a major source of new and emerging viral diseases. Despite the fact that bats carry and shed highly pathogenic viruses including Ebola, Nipah and SARS, they rarely display clinical symptoms of infection. Host factors influencing viral replication are poorly understood in bats and are likely to include both pre- and post-transcriptional regulatory mechanisms. MicroRNAs are a major mechanism of post-transcriptional gene regulation, however very little is known about them in bats.

Results

This study describes 399 microRNAs identified by deep sequencing of small RNA isolated from tissues of the Black flying fox, Pteropus alecto, a confirmed natural reservoir of the human pathogens Hendra virus and Australian bat lyssavirus. Of the microRNAs identified, more than 100 are unique amongst vertebrates, including a subset containing mutations in critical seed regions. Clusters of rapidly-evolving microRNAs were identified, as well as microRNAs predicted to target genes involved in antiviral immunity, the DNA damage response, apoptosis and autophagy. Closer inspection of the predicted targets for several highly supported novel miRNA candidates suggests putative roles in host-virus interaction.

Conclusions

MicroRNAs are likely to play major roles in regulating virus-host interaction in bats, via dampening of inflammatory responses (limiting the effects of immunopathology), and directly limiting the extent of viral replication, either through restricting the availability of essential factors or by controlling apoptosis. Characterisation of the bat microRNA repertoire is an essential step towards understanding transcriptional regulation during viral infection, and will assist in the identification of mechanisms that enable bats to act as natural virus reservoirs. This in turn will facilitate the development of antiviral strategies for use in humans and other species.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-682) contains supplementary material, which is available to authorized users.

Structured RNAs and synteny regions in the pig genome

Background

Annotating mammalian genomes for noncoding RNAs (ncRNAs) is nontrivial since far from all ncRNAs are known and the computational models are resource demanding. Currently, the human genome holds the best mammalian ncRNA annotation, a result of numerous efforts by several groups. However, a more direct strategy is desired for the increasing number of sequenced mammalian genomes of which some, such as the pig, are relevant as disease models and production animals.

Results

We present a comprehensive annotation of structured RNAs in the pig genome. Combining sequence and structure similarity search as well as class specific methods, we obtained a conservative set with a total of 3,391 structured RNA loci of which 1,011 and 2,314, respectively, hold strong sequence and structure similarity to structured RNAs in existing databases. The RNA loci cover 139 cis-regulatory element loci, 58 lncRNA loci, 11 conflicts of annotation, and 3,183 ncRNA genes. The ncRNA genes comprise 359 miRNAs, 8 ribozymes, 185 rRNAs, 638 snoRNAs, 1,030 snRNAs, 810 tRNAs and 153 ncRNA genes not belonging to the here fore mentioned classes. When running the pipeline on a local shuffled version of the genome, we obtained no matches at the highest confidence level. Additional analysis of RNA-seq data from a pooled library from 10 different pig tissues added another 165 miRNA loci, yielding an overall annotation of 3,556 structured RNA loci. This annotation represents our best effort at making an automated annotation. To further enhance the reliability, 571 of the 3,556 structured RNAs were manually curated by methods depending on the RNA class while 1,581 were declared as pseudogenes. We further created a multiple alignment of pig against 20 representative vertebrates, from which RNAz predicted 83,859 de novo RNA loci with conserved RNA structures. 528 of the RNAz predictions overlapped with the homology based annotation or novel miRNAs. We further present a substantial synteny analysis which includes 1,004 lineage specific de novo RNA loci and 4 ncRNA loci in the known annotation specific for Laurasiatheria (pig, cow, dolphin, horse, cat, dog, hedgehog).

Conclusions

We have obtained one of the most comprehensive annotations for structured ncRNAs of a mammalian genome, which is likely to play central roles in both health modelling and production. The core annotation is available in Ensembl 70 and the complete annotation is available at http://rth.dk/resources/rnannotator/susscr102/version1.02.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-459) contains supplementary material, which is available to authorized users.

Virus versus host plant microRNAs: who determines the outcome of the interaction?

Considering the importance of microRNAs (miRNAs) in the regulation of essential processes in plant pathogen interactions, it is not surprising that, while plant miRNA sequences counteract viral attack via antiviral RNA silencing, viruses in turn have developed antihost defense mechanisms blocking these RNA silencing pathways and establish a counter-defense. In the current study, computational and stem-loop Reverse Transcription – Polymerase Chain Reaction (RT-PCR) approaches were employed to a) predict and validate virus encoded mature miRNAs (miRs) in 39 DNA-A sequences of the bipartite genomes of African cassava mosaic virus (ACMV) and East African cassava mosaic virus-Uganda (EACMV-UG) isolates, b) determine whether virus encoded miRs/miRs* generated from the 5′/3′ harpin arms have the capacity to bind to genomic sequences of the host plants Jatropha or cassava and c) investigate whether plant encoded miR/miR* sequences have the potential to bind to the viral genomes. Different viral pre-miRNA hairpin sequences and viral miR/miR* length variants occurring as isomiRs were predicted in both viruses. These miRNAs were located in three Open Reading Frames (ORFs) and in the Intergenic Region (IR). Moreover, various target genes for miRNAs from both viruses were predicted and annotated in the host plant genomes indicating that they are involved in biotic response, metabolic pathways and transcription factors. Plant miRs/miRs* from conserved and highly expressed families were identified, which were shown to have potential targets in the genome of both begomoviruses, representing potential plant miRNAs mediating antiviral defense. This is the first assessment of predicted viral miRs/miRs* of ACMV and EACMV-UG and host plant miRNAs, providing a reference point for miRNA identification in pathogens and their hosts. These findings will improve the understanding of host- pathogen interaction pathways and the function of viral miRNAs in Euphorbiaceous crop plants.

Identification and characterization of microRNAs in normal equine tissues by Next Generation Sequencing

The role of microRNAs (miRNAs) as a post-transcriptional gene regulator has been elucidated in a broad range of organisms including domestic animals. Characterization of miRNAs in normal tissues is an important step to investigate the functions of miRNAs in various physiological and pathological conditions. Using Illumina Next Generation Sequencing (NGS) technology, we identified a total of 292 known and 329 novel miRNAs in normal horse tissues including skeletal muscle, colon and liver. Distinct sets of miRNAs were differentially expressed in a tissue-specific manner. The miRNA genes were distributed across all the chromosomes except chromosomes 29 and 31 in the horse reference genome. In some chromosomes, multiple miRNAs were clustered and considered to be polycistronic transcript. A base composition analysis showed that equine miRNAs had a higher frequency of A+U than G+C. Furthermore, U tended to be more frequent at the 5′ end of miRNA sequences. This is the first experimental study that identifies and characterizes the global miRNA expression profile in normal horse tissues. The present study enriches the horse miRNA database and provides useful information for further research dissecting biological functions of miRNAs in horse.

Identification of novel microRNAs in primates by using the synteny information and small RNA deep sequencing data

Current technologies that are used for genome-wide microRNA (miRNA) prediction are mainly based on BLAST tool. They often produce a large number of false positives. Here, we describe an effective approach for identifying orthologous pre-miRNAs in several primates based on syntenic information. Some of them have been validated by small RNA high throughput sequencing data. This approach uses the synteny information and experimentally validated miRNAs of human, and incorporates currently available algorithms and tools to identify the pre-miRNAs in five other primates. First, we identified 929 potential pre-miRNAs in the marmoset in which miRNAs have not yet been reported. Then, we predicted the miRNAs in other primates, and we successfully re-identified most of the published miRNAs and found 721, 979, 650 and 639 new potential pre-miRNAs in chimpanzee, gorilla, orangutan and rhesus macaque, respectively. Furthermore, the miRNA transcriptome in the four primates have been re-analyzed and some novel predicted miRNAs have been supported by the small RNA sequencing data. Finally, we analyzed the potential functions of those validated miRNAs and explored the regulatory elements and transcription factors of some validated miRNA genes of interest. The results show that our approach can effectively identify novel miRNAs and some miRNAs that supported by small RNA sequencing data maybe play roles in the nervous system.

MicroRNAs in domestic livestock

microRNAs (miRNAs) are a class of small noncoding RNA that bind to complementary sequences in the untranslated regions of multiple target mRNAs resulting in posttranscriptional regulation of gene expression. The recent discovery and expression-profiling studies of miRNAs in domestic livestock have revealed both their tissue-specific and temporal expression pattern. In addition, breed-dependent expression patterns as well as single nucleotide polymorphisms in either the miRNA or in the target mRNA binding site have revealed associations with traits of economic importance and highlight the potential use of miRNAs in future genomic selection programs.

Stallion sperm transcriptome comprises functionally coherent coding and regulatory RNAs as revealed by microarray analysis and RNA-seq

Mature mammalian sperm contain a complex population of RNAs some of which might regulate spermatogenesis while others probably play a role in fertilization and early development. Due to this limited knowledge, the biological functions of sperm RNAs remain enigmatic. Here we report the first characterization of the global transcriptome of the sperm of fertile stallions. The findings improved understanding of the biological significance of sperm RNAs which in turn will allow the discovery of sperm-based biomarkers for stallion fertility. The stallion sperm transcriptome was interrogated by analyzing sperm and testes RNA on a 21,000-element equine whole-genome oligoarray and by RNA-seq. Microarray analysis revealed 6,761 transcripts in the sperm, of which 165 were sperm-enriched, and 155 were differentially expressed between the sperm and testes. Next, 70 million raw reads were generated by RNA-seq of which 50% could be aligned with the horse reference genome. A total of 19,257 sequence tags were mapped to all horse chromosomes and the mitochondrial genome. The highest density of mapped transcripts was in gene-rich ECA11, 12 and 13, and the lowest in gene-poor ECA9 and X; 7 gene transcripts originated from ECAY. Structural annotation aligned sperm transcripts with 4,504 known horse and/or human genes, rRNAs and 82 miRNAs, whereas 13,354 sequence tags remained anonymous. The data were aligned with selected equine gene models to identify additional exons and splice variants. Gene Ontology annotations showed that sperm transcripts were associated with molecular processes (chemoattractant-activated signal transduction, ion transport) and cellular components (membranes and vesicles) related to known sperm functions at fertilization, while some messenger and micro RNAs might be critical for early development. The findings suggest that the rich repertoire of coding and non-coding RNAs in stallion sperm is not a random remnant from spermatogenesis in testes but a selectively retained and functionally coherent collection of RNAs.

Comparative genomic analysis reveals evolutionary characteristics and patterns of microRNA clusters in vertebrates

MicroRNAs (miRNAs) are a class of small non-coding RNAs that can play important regulatory roles in many important biological processes. Although clustering patterns of miRNA clusters have been uncovered in animals, the origin and evolution of miRNA clusters in vertebrates are still poorly understood. Here, we performed comparative genomic analyses to construct 51 sets of orthologous miRNA clusters (SOMCs) across seven test vertebrate species, a collection of miRNA clusters from two or more species that are likely to have evolved from a common ancestral miRNA cluster, and used these to systematically examine the evolutionary characteristics and patterns of miRNA clusters in vertebrates. We found that miRNA clusters are continuously generated, and most of them tend to be conserved and maintained in vertebrate genomes, although some adaptive gains and losses of miRNA cluster have occurred during evolution. Furthermore, miRNA clusters appeared relatively early in the evolutionary history might suffer from more complicated adaptive gain-and-loss than those young miRNA clusters. Detailed analysis showed that genomic duplication events of ancestral miRNAs or miRNA clusters are likely to be major driving force and apparently contribute to origin and evolution of miRNA clusters. Comparison of conserved with lineage-specific miRNA clusters revealed that the contribution of duplication events for the formation of miRNA cluster appears to be more important for conserved miRNA clusters than lineage-specific. Our study provides novel insights for further exploring the origins and evolution of miRNA clusters in vertebrates at a genome scale.

Muscular microRNA expressions in healthy and myopathic horses suffering from polysaccharide storage myopathy or recurrent exertional rhabdomyolysis

REASONS FOR PERFORMING STUDY

MicroRNAs (miRNA) are small endogenous noncoding interfering RNA molecules (18-25 nucleotides) regarded as major regulators in eukaryotic gene expression. They play a role in developmental timing, cellular differentiation, signalling and apoptosis pathways. Because of the central function of miRNAs in the proliferation and differentiation of the myoblasts demonstrated in mouse and man, it is assumed that they could be present in equine muscles and their expression profile may be related to the muscle status.

OBJECTIVE

To identify miRNA candidates in the muscles of control and affected horses suffering from polysaccharide storage myopathy (PSSM) and recurrent exertional rhabdomyolysis (RER).

METHODS

Muscle biopsies were collected in the gluteus medius of horses allocated into 4 groups: French Trotters (3 control-TF vs. 3 RER-TF) and Norman Cob (5 control-Cob vs. 9 PSSM-Cob). Blood samples were collected for miRNA analysis. Total RNA were extracted and real time quantitative RT-QPCR analysis were conducted using 10 miRNA assays (mir-1-23-30-133-181-188-195-206-339-375).

RESULTS

All the miRNA candidates were significantly detected in the muscles and some in blood samples. Variance analysis revealed highly significant (P < 0.0001) effects of the miRNA type, breed and pathology on the miRNA expression. A specific miRNA profile was related to each myopathy: a higher expression of mir-1, 133, 23a, 30b, 195 and 339 in RER-TF vs. control-TF (P < 0.05); a higher expression of mir-195 in PSSM-Cob vs. control-Cob (P < 0.05). The miRNA profile was different between breeds for mir-181, 188 and 206 (P < 0.05). The mir-1, 133, 181, 195 and 206 were detected in blood of control-Cob and PSSM-Cob horses.

CONCLUSIONS

This first study about muscular miRNA profile in equine myopathies indicated that it is possible to discriminate pathological from control horses according to their miRNA profile. The RER miRNA profile was more specific and contrasted than the PSSM profile.

Genome-wide identification of novel microRNAs and their target genes in the human parasite Schistosoma mansoni

Mature microRNAs (miRNAs) are small, non-coding regulatory RNAs which can elicit post-transcriptional repression of mRNA levels of target genes. Here, we report the identification of 67 mature and 42 precursor miRNAs in the Schistosoma mansoni parasite. The evolutionarily conserved S. mansoni miRNAs consisted of 26 precursor miRNAs and 35 mature miRNAs, while we identified 16 precursor miRNAs and 32 mature miRNAs that displayed no conservation. These S. mansoni miRNAs are located on seven autosomal chromosomes and a sex (W) chromosome. miRNA expansion through gene duplication was suggested for at least two miRNA families miR-71 and mir-2. miRNA target finding analysis identified 389 predicted mRNA targets for the identified miRNAs and suggests that the sma-mir-71 may be involved in female sexual maturation. Given the important roles of miRNAs in animals, the identification and characterization of miRNAs in S. mansoni will facilitate novel approaches towards prevention and treatment of Schistosomiasis.

Origin and evolution of a placental-specific microRNA family in the human genome

Background

MicroRNAs (miRNAs) are a class of short regulatory RNAs encoded in the genome of DNA viruses, some single cell organisms, plants and animals. With the rapid development of technology, more and more miRNAs are being discovered. However, the origin and evolution of most miRNAs remain obscure. Here we report the origin and evolution dynamics of a human miRNA family.

Results

We have shown that all members of the miR-1302 family are derived from MER53 elements. Although the conservation scores of the MER53-derived pre-miRNA sequences are low, we have identified 36 potential paralogs of MER53-derived miR-1302 genes in the human genome and 58 potential orthologs of the human miR-1302 family in placental mammals. We suggest that in placental species, this miRNA family has evolved following the birth-and-death model of evolution. Three possible mechanisms that can mediate miRNA duplication in evolutionary history have been proposed: the transposition of the MER53 element, segmental duplications and Alu-mediated recombination. Finally, we have found that the target genes of miR-1302 are over-represented in transportation, localization, and system development processes and in the positive regulation of cellular processes. Many of them are predicted to function in binding and transcription regulation.

Conclusions

The members of miR-1302 family that are derived from MER53 elements are placental-specific miRNAs. They emerged at the early stage of the recent 180 million years since eutherian mammals diverged from marsupials. Under the birth-and-death model, the miR-1302 genes have experienced a complex expansion with some members evolving by segmental duplications and some by Alu-mediated recombination events.

A genome-wide SNP-association study confirms a sequence variant (g.66493737C>T) in the equine myostatin (MSTN) gene as the most powerful predictor of optimum racing distance for Thoroughbred racehorses

Background

Thoroughbred horses have been selected for traits contributing to speed and stamina for centuries. It is widely recognized that inherited variation in physical and physiological characteristics is responsible for variation in individual aptitude for race distance, and that muscle phenotypes in particular are important.

Results

A genome-wide SNP-association study for optimum racing distance was performed using the EquineSNP50 Bead Chip genotyping array in a cohort of n = 118 elite Thoroughbred racehorses divergent for race distance aptitude. In a cohort-based association test we evaluated genotypic variation at 40,977 SNPs between horses suited to short distance (≤ 8 f) and middle-long distance (> 8 f) races. The most significant SNP was located on chromosome 18: BIEC2-417495 ~690 kb from the gene encoding myostatin (MSTN) [P_unadj. = 6.96 × 10^-6]. Considering best race distance as a quantitative phenotype, a peak of association on chromosome 18 (chr18:65809482-67545806) comprising eight SNPs encompassing a 1.7 Mb region was observed. Again, similar to the cohort-based analysis, the most significant SNP was BIEC2-417495 (P_unadj. = 1.61 × 10^-9; P_Bonf. = 6.58 × 10^-5). In a candidate gene study we have previously reported a SNP (g.66493737C>T) in MSTN associated with best race distance in Thoroughbreds; however, its functional and genome-wide relevance were uncertain. Additional re-sequencing in the flanking regions of the MSTN gene revealed four novel 3' UTR SNPs and a 227 bp SINE insertion polymorphism in the 5' UTR promoter sequence. Linkage disequilibrium was highest between g.66493737C>T and BIEC2-417495 (r² = 0.86).

Conclusions

Comparative association tests consistently demonstrated the g.66493737C>T SNP as the superior variant in the prediction of distance aptitude in racehorses (g.66493737C>T, P = 1.02 × 10^-10; BIEC2-417495, P_unadj. = 1.61 × 10^-9). Functional investigations will be required to determine whether this polymorphism affects putative transcription-factor binding and gives rise to variation in gene and protein expression. Nonetheless, this study demonstrates that the g.66493737C>T SNP provides the most powerful genetic marker for prediction of race distance aptitude in Thoroughbreds.