Mapping and expression studies of the mir17-92 cluster on pig Chromosome 11

We have identified the first porcine microRNA (miRNA) cluster (the mir17-92 cluster) and localized it to the q-arm of pig Chromosome 11. The miRNA cluster was found by sequence similarity search with human miRNA sequences against the pig genomic data generated within the Sino-Danish pig genome project. The resulting data contained three complete and two incomplete miRNA precursors of seven miRNAs from the human mir17-92 cluster. Because there is a 100% sequence identity between the four pig miRNAs and the corresponding human miRNAs, the sequences of three unavailable pig miRNAs were derived from the human data. The expression profiles of seven studied miRNAs were analyzed by hybridization to Northern blots containing five porcine tissues: cerebellum, cortex, hippocampus, kidney, and liver. In order to determine the localization of the mir17-92 cluster in the pig genome, we mapped it by PCR in the porcine somatic cell hybrid (SCH) panel and in the INRA-University of Minnesota (INRA-UMN) porcine radiation hybrid (IMpRH) panel. The PCR results enabled us to localize this cluster to the q-arm of pig Chromosome 11 and map it in relation to two microsatellites. Our study presents the first expression analyses of miRNAs in pig and adds information for further functional studies in this species.     

Suitable internal control microRNA genes for measuring miRNA abundance in pig milk during different lactation periods

Determination of an optimal set/number of internal control microRNA (miRNA) genes is a critical, but often undervalued, detail of quantitative gene expression analysis. No validated internal genes for miRNA quantitative PCR (q-PCR) in pig milk were available. We compared the expression stability of six porcine internal control miRNA genes in pig milk from different lactation periods (1 h, 3 days, 7 days, 14 days, 21 days, and 28 days postpartum), using an EvaGreen q-PCR approach. We found that using the three most stable internal control genes to calculate the normalization factor is sufficient for producing reliable q-PCR expression data. We also found that miRNAs are superior to ribosomal RNA (rRNA) and snRNA, which are commonly used as internal controls for normalizing miRNA q-PCR data. In terms of economic and experimental feasibility, we recommend the use of the three most stable internal control miRNA genes (miR-17, -107 and -103) for calculating the normalization factors for pig milk samples from different lactation periods. These results can be applied to future studies aimed at measuring miRNA abundance in porcine milk.     

Identification and Characterization of the miRNA Transcriptome of Ovis aries

The discovery and identification of Ovis aries (sheep) miRNAs will further promote the study of miRNA functions and gene regulatory mechanisms. To explore the microRNAome (miRNAome) of sheep in depth, samples were collected that included eight developmental stages: the longissimus dorsi muscles of Texel fetuses at 70, 85, 100, 120, and 135 days, and the longissimus dorsi muscles of Ujumqin fetuses at 70, 85, 100, 120, and 135 d, and lambs at 0 (birth), 35, and 70 d. These samples covered all of the representative periods of Ovis aries growth and development throughout gestation (about 150 d) and 70 d after birth. Texel and Ujumqin libraries were separately subjected to Solexa deep sequencing; 35,700,772 raw reads were obtained overall. We used ACGT101-miR v4.2 to analyze the sequence data. Following meticulous comparisons with mammalian mature miRNAs, precursor hairpins (pre-miRNAs), and the latest sheep genome, we substantially extended the Ovis aries miRNAome. The list of pre-miRNAs was extended to 2,319, expressing 2,914 mature miRNAs. Among those, 1,879 were genome mapped to unique miRNAs, representing 2,436 genome locations, and 1,754 pre-miRNAs were mapped to chromosomes. Furthermore, the Ovis aries miRNAome was processed using an elaborate bioinformatic analysis that examined multiple end sequence variation in miRNAs, precursors, chromosomal localizations, species-specific expressions, and conservative properties. Taken together, this study provides the most comprehensive and accurate exploration of the sheep miRNAome, and draws conclusions about numerous characteristics of Ovis aries miRNAs, including miRNAs and isomiRs.     

Global profiling of miRNAs and the hairpin precursors: insights into miRNA processing and novel miRNA discovery

MicroRNAs (miRNAs) constitute an important class of small regulatory RNAs that are derived from distinct hairpin precursors (pre-miRNAs). In contrast to mature miRNAs, which have been characterized in numerous genome-wide studies of different organisms, research on global profiling of pre-miRNAs is limited. Here, using massive parallel sequencing, we have performed global characterization of both mouse mature and precursor miRNAs. In total, 87 369 704 and 252 003 sequencing reads derived from 887 mature and 281 precursor miRNAs were obtained, respectively. Our analysis revealed new aspects of miRNA/pre-miRNA processing and modification, including eight Ago2-cleaved pre-miRNAs, eight new instances of miRNA editing and exclusively 5′ tailed mirtrons. Furthermore, based on the sequences of both mature and precursor miRNAs, we developed a miRNA discovery pipeline, miRGrep, which does not rely on the availability of genome reference sequences. In addition to 239 known mouse pre-miRNAs, miRGrep predicted 41 novel ones with high confidence. Similar as known ones, the mature miRNAs derived from most of these novel loci showed both reduced abundance following Dicer knockdown and the binding with Argonaute2. Evaluation on data sets obtained from Caenorhabditis elegans and Caenorhabditis sp.11 demonstrated that miRGrep could be widely used for miRNA discovery in metazoans, especially in those without genome reference sequences.     

Solexa sequencing identification of conserved and novel microRNAs in backfat of Large White and Chinese Meishan pigs

The domestic pig (Sus scrofa), an important species in animal production industry, is a right model for studying adipogenesis and fat deposition. In order to expand the repertoire of porcine miRNAs and further explore potential regulatory miRNAs which have influence on adipogenesis, high-throughput Solexa sequencing approach was adopted to identify miRNAs in backfat of Large White (lean type pig) and Meishan pigs (Chinese indigenous fatty pig). We identified 215 unique miRNAs comprising 75 known pre-miRNAs, of which 49 miRNA*s were first identified in our study, 73 miRNAs were overlapped in both libraries, and 140 were novelly predicted miRNAs, and 215 unique miRNAs were collectively corresponding to 235 independent genomic loci. Furthermore, we analyzed the sequence variations, seed edits and phylogenetic development of the miRNAs. 17 miRNAs were widely conserved from vertebrates to invertebrates, suggesting that these miRNAs may serve as potential evolutional biomarkers. 9 conserved miRNAs with significantly differential expressions were determined. The expression of miR-215, miR-135, miR-224 and miR-146b was higher in Large White pigs, opposite to the patterns shown by miR-1a, miR-133a, miR-122, miR-204 and miR-183. Almost all novel miRNAs could be considered pig-specific except ssc-miR-1343, miR-2320, miR-2326, miR-2411 and miR-2483 which had homologs in Bos taurus, among which ssc-miR-1343, miR-2320, miR-2411 and miR-2483 were validated in backfat tissue by stem-loop qPCR. Our results displayed a high level of concordance between the qPCR and Solexa sequencing method in 9 of 10 miRNAs comparisons except for miR-1a. Moreover, we found 2 miRNAs, miR-135 and miR-183, may exert impacts on porcine backfat development through WNT signaling pathway. In conclusion, our research develops porcine miRNAs and should be beneficial to study the adipogenesis and fat deposition of different pig breeds based on miRNAs.     

Identification of microRNA precursors with support vector machine and string kernel

MicroRNAs （miRNAs） are one family of short （21-23 nt） regulatory non-coding RNAs processed from long （70-110 nt） miRNA precursors （pre-miRNAs）. Identifying true and false precursors plays an important role in computational identification of miRNAs. Some numerical features have been extracted from precursor sequences and their secondary structures to suit some classification methods; however, they may lose some usefully discriminative information hidden in sequences and structures. In this study, pre-miRNA sequences and their secondary structures are directly used to construct an exponential kernel based on weighted Levenshtein distance between two sequences. This string kernel is then combined with support vector machine （SVM） for detecting true and false pre-miRNAs. Based on 331 training samples of true and false human pre-miRNAs, 2 key parameters in SVM are selected by 5-fold cross validation and grid search, and 5 realizations with different 5-fold partitions are executed. Among 16 independent test sets from 3 human, 8 animal, 2 plant, 1 virus, and 2 artificially false human pre-miRNAs, our method statistically outperforms the previous SVM-based technique on 11 sets, including 3 human, 7 animal, and 1 false human pre-miRNAs. In particular, premiRNAs with multiple loops that were usually excluded in the previous work are correctly identified in this study with an accuracy of 92.66%.     

Discovery of novel microRNAs in rat kidney using next generation sequencing and microarray validation

MicroRNAs (miRNAs) are small non-coding RNAs that regulate a variety of biological processes. The latest version of the miRBase database (Release 18) includes 1,157 mouse and 680 rat mature miRNAs. Only one new rat mature miRNA was added to the rat miRNA database from version 16 to version 18 of miRBase, suggesting that many rat miRNAs remain to be discovered. Given the importance of rat as a model organism, discovery of the completed set of rat miRNAs is necessary for understanding rat miRNA regulation. In this study, next generation sequencing (NGS), microarray analysis and bioinformatics technologies were applied to discover novel miRNAs in rat kidneys. MiRanalyzer was utilized to analyze the sequences of the small RNAs generated from NGS analysis of rat kidney samples. Hundreds of novel miRNA candidates were examined according to the mappings of their reads to the rat genome, presence of sequences that can form a miRNA hairpin structure around the mapped locations, Dicer cleavage patterns, and the levels of their expression determined by both NGS and microarray analyses. Nine novel rat hairpin precursor miRNAs (pre-miRNA) were discovered with high confidence. Five of the novel pre-miRNAs are also reported in other species while four of them are rat specific. In summary, 9 novel pre-miRNAs (14 novel mature miRNAs) were identified via combination of NGS, microarray and bioinformatics high-throughput technologies.     

Tissue-Specific Expression and Regulatory Networks of Pig MicroRNAome

Background

Despite the economic and medical importance of the pig, knowledge about its genome organization, gene expression regulation, and molecular mechanisms involved in physiological processes is far from that achieved for mouse and rat, the two most used model organisms in biomedical research. MicroRNAs (miRNAs) are a wide class of molecules that exert a recognized role in gene expression modulation, but only 280 miRNAs in pig have been characterized to date.

Results

We applied a novel computational approach to predict species-specific and conserved miRNAs in the pig genome, which were then subjected to experimental validation. We experimentally identified candidate miRNAs sequences grouped in high-confidence (424) and medium-confidence (353) miRNAs according to RNA-seq results. A group of miRNAs was also validated by PCR experiments. We established the subtle variability in expression of isomiRs and miRNA-miRNA star couples supporting a biological function for these molecules. Finally, miRNA and mRNA expression profiles produced from the same sample of 20 different tissue of the animal were combined, using a correlation threshold to filter miRNA-target predictions, to identify tissue-specific regulatory networks.

Conclusions

Our data represent a significant progress in the current understanding of miRNAome in pig. The identification of miRNAs, their target mRNAs, and the construction of regulatory circuits will provide new insights into the complex biological networks in several tissues of this important animal model.     

Epigenetic regulation of microRNAs in cancer: an integrated review of literature

MicroRNAs (miRNAs) belong to the heterogeneous class of non-coding RNAs (ncRNAs) that regulate the translation and degradation of target mRNAs, and control approximately 30% of human genes. MiRNA genes might be silenced in human tumors (oncomiRs) by aberrant hypermethylation of CpG islands that encompass or lie adjacent to miRNA genes and/or by histone modifications. We performed literature search for research articles describing epigenetically regulated miRNAs in cancer and identified 45 studies that were published between 2006 and 7/2010. The data from those papers are fragmented and methodologically heterogeneous and our work represents first systematic review towards to integration of diverse sets of information. We reviewed the methods used for detection of miRNA epigenetic regulation, which comprise bisulfite genomic sequencing PCR (BSP), bisulfite pyrosequencing, methylation specific PCR (MSP), combined bisulfite restriction analysis (COBRA), methylation sensitive single nucleotide primer extension (Ms-SNuPE), MassARRAY technique and some modifications of those methods. This integrative study revealed 122 miRNAs that were reported to be epigenetically regulated in 23 cancer types. Compared to protein coding genes, human oncomiRs showed an order of magnitude higher methylation frequency (11.6%; 122/1048 known miRNAs). Nearly half, (45%; 55/122) epigenetically regulated miRNAs were associated with different cancer types, but other 55% (67/122) miRNAs were present in only one cancer type and therefore representing cancer-specific biomarker potential. The data integration revealed miRNA epigenomic hot spots on the chromosomes 1q, 7q, 11q, 14q and 19q. CpG island analysis of corresponding miRNA precursors (pre-miRNAs) revealed that 20% (26/133) of epigenetically regulated miRNAs had a CpG island within the range of 5kb upstream, among them 14% (19/133) of miRNAs resided within the CpG island. Our integrative survey and analyses revealed candidate cancer-specific miRNA epigenetic signatures which provide the basis for new therapeutic strategies in cancer by targeting the epigenetic regulation of miRNAs.     

The miRNA aberrant expression dependence on DNA methylation in HeLa cells treated with mitomycin C

The dependence of expression of miRNAs and their precursors (pre-miRNAs) on the DNA methylation level in HeLa cells 8 days after mitomycin C treatment was studied. A massive parallel DNA sequencing method was applied to analyze miRNA expression. 5-Azacytidine (DNA methylation inhibitor) was added to the medium 6 days after mutagenic agent exposure. The results indicated that the change in expression for some mature miRNAs (39 of 61) was accompanied by the change in the expression of their pre-miRNAs, while there were no significant changes in the expression of pre-miRNA for other mature miRNAs (22 of 61). The aberrant expression was maintained by 8 of 61 mature miRNAs and 6 of 55 pre-miRNAs in the induced HeLa cells after 5-azacytidine treatment. In addition, the expression of more than 90% of miRNAs, which indicated a significant change in expression after mitomycin C treatment, does not depend or depends slightly on the DNA methylation level in HeLa cells without mitomycin C treatment. The results suggest that mitomycin C induces aberrant DNA methylation which affects maintenance of changes in the miRNA expression in cell generations after mutagen treatment.