Identification of the highly accumulated microRNA*s in Arabidopsis (Arabidopsis thaliana) and rice (Oryza sativa)

Plant microRNAs (miRNAs) are crucial for the regulation of gene expression, which is involved in almost all the important biological processes. In the cytoplasm, the miRNA strand is selectively incorporated into a specific Argonaute (AGO)-associated gene silencing complex, while the miRNA* is degraded rapidly. Thus, most miRNA*s were thought to be biologically meaningless. Interestingly, several recent reports in both plants and animals have shaken this notion. Many miRNA*s were demonstrated to possess regulatory roles in gene expression. However, the low accumulation levels of most miRNA*s raise the question whether the activities of this small RNA (sRNA) species are widespread in plants. Here, by using publicly available sRNA high-throughput sequencing data, we found that the accumulation levels of several miRNA*s could be much higher than those of their miRNA partners in certain organs, mutants and/or AGO-associated silencing complexes of both Arabidopsis (Arabidopsis thaliana) and rice (Oryza sativa). Based on target prediction and degradome sequencing data-based validation, some of these highly accumulated miRNA*s were indicated to possess cleavage-based potential regulatory role on certain targets. Besides, some interesting biological interpretations were obtained based on the accumulation patterns of the miRNA*s, the annotations of the target genes, and literature mining. Taken together, the expanded list of the highly accumulated miRNA*s along with their potential target genes discovered in this study further strengthened the current notion that certain members of the miRNA* species are biologically relevant, which needs further inspection.     

Genetic framework for flowering-time regulation by ambient temperature-responsive miRNAs in Arabidopsis

Flowering is the primary trait affected by ambient temperature changes. Plant microRNAs (miRNAs) are small non-coding RNAs playing an important regulatory role in plant development. In this study, to elucidate the mechanism of flowering-time regulation by small RNAs, we identified six ambient temperature-responsive miRNAs (miR156, miR163, miR169, miR172, miR398 and miR399) in Arabidopsis via miRNA microarray and northern hybridization analyses. We also determined the expression profile of 120 unique miRNA loci in response to ambient temperature changes by miRNA northern hybridization analysis. The expression of the ambient temperature-responsive miRNAs and their target genes was largely anticorrelated at two different temperatures (16 and 23°C). Interestingly, a lesion in short vegetative phase (SVP), a key regulator within the thermosensory pathway, caused alteration in the expression of miR172 and a subset of its target genes, providing a link between a thermosensory pathway gene and miR172. The miR172-overexpressing plants showed a temperature-independent early flowering phenotype, suggesting that modulation of miR172 expression leads to temperature insensitivity. Taken together, our results suggest a genetic framework for flowering-time regulation by ambient temperature-responsive miRNAs under non-stress temperature conditions.     

The mRNA–miRNA–lncRNA Regulatory Network and Factors Associated with Prognosis Prediction of Hepatocellular Carcinoma

The aim of this study was to identify novel prognostic mRNA and microRNA (miRNA) biomarkers for hepatocellular carcinoma (HCC) using methods in systems biology. Differentially expressed mRNAs, miRNAs, and long non-coding RNAs (lncRNAs) were compared between HCC tumor tissues and normal liver tissues in The Cancer Genome Atlas (TCGA) database. Subsequently, a prognosis-associated mRNA co-expression network, an mRNA–miRNA regulatory network, and an mRNA–miRNA–lncRNA regulatory network were constructed to identify prognostic biomarkers for HCC through Cox survival analysis. Seven prognosis-associated mRNA co-expression modules were obtained by analyzing these differentially expressed mRNAs. An expression module including 120 mRNAs was significantly correlated with HCC patient survival. Combined with patient survival data, several mRNAs and miRNAs, including CHST4, SLC22A8, STC2, hsa-miR-326, and hsa-miR-21 were identified from the network to predict HCC patient prognosis. Clinical significance was investigated using tissue microarray analysis of samples from 258 patients with HCC. Functional annotation of hsa-miR-326 and hsa-miR-21-5p indicated specific associations with several cancer-related pathways. The present study provides a bioinformatics method for biomarker screening, leading to the identification of an integrated mRNA–miRNA–lncRNA regulatory network and their co-expression patterns in relation to predicting HCC patient survival.     

Identification of novel microRNA-like-coding sites on the long-stem microRNA precursors in Arabidopsis

Plant microRNA (miRNA) is a crucial regulator of gene expression. It has been reported that more than one miRNA/miRNA* duplex could be produced from a microRNA precursor (pre-miRNA). In this study, we performed a comprehensive search for the novel miRNA candidates on the pre-miRNAs of Arabidopsis. AGO1 enrichment, co-existence of the miRNA*-like coordinates, and unique genome-wide match sites were taken into consideration for candidate screening. As a result, 43 miRNA-like candidates derived from 25 pre-miRNAs were identified. Among these candidates, 31 strong candidates from 22 pre-miRNAs passed all the filtering steps. Interestingly, some of these miRNA-like candidates showed organ-specific expression patterns. After target prediction and degradome sequencing data-based validation, five miRNA candidate–target pairs (ath-miR863-5p.2–AT1G76550.1, ath-miR822.2–AT5G03552.1, ath-miR822.3–AT5G02350.1, sRNA4–AT1G66290.1 and sRNA6–AT1G66310.1) were identified, providing a basis for in-depth functional analysis of these miRNA candidates. These results could update the current understanding of the biogenesis and the action of the plant miRNAs.     

MicroRNA-target gene responses to root knot nematode (Meloidogyne incognita) infection in cotton (Gossypium hirsutum L.)

MicroRNAs (miRNAs) are a large class of small regulatory RNA molecules, however no study has been performed to elucidate the role of miRNAs in cotton (Gossypium hirsutum) response to the root knot nematode (RKN, Meloidogyne incognita) infection. We selected 28 miRNAs and 8 miRNA target genes to investigate the miRNA-target gene response to M. incognita infection. Our results show that RKN infection significantly affected the expression of several miRNAs and their targeted genes. After 10 days of RKN infection, expression fold changes on miRNA expressions ranged from down-regulated by 33% to upregulated by 406%; meanwhile the expression levels of miRNA target genes were 45.8% to 231%. Three miRNA-target pairs, miR159-MYB, miR319-TCP4 and miR167-ARF8, showed inverse expression patterns between gene targets and their corresponded miRNAs, suggesting miRNA-mediated gene regulation in cotton roots in response to RKN infection.     

Genome-Wide Identification of Reverse Complementary microRNA Genes in Plants

MicroRNAs (miRNAs) are ∼21-nucleotide small RNAs (sRNAs) with essential regulatory roles in plants. They are generated from stem-loop-structured precursors through two sequential Dicer-like 1 (DCL1)-mediated cleavages. To date, hundreds of plant miRNAs have been uncovered. However, the question, whether the sequences reverse complementary (RC) to the miRNA precursors could form hairpin-like structures and produce sRNA duplexes similar to the miRNA/miRNA* pairs has not been solved yet. Here, we interrogated this possibility in 16 plant species based on sRNA high-throughput sequencing data and secondary structure prediction. A total of 59 RC sequences with great potential to form stem-loop structures and generate miRNA/miRNA*-like duplexes were identified in ten plants, which were named as RC-miRNA precursors. Unlike the canonical miRNAs, only a few cleavage targets of the RC-miRNAs were identified in Arabidopsis (Arabidopsis thaliana) and rice (Oryza sativa), and none in Soybean (Glycine max) based on degradome data. Surprisingly, the genomic regions surrounding some of the RC-miRNA target recognition sites were observed to be specifically methylated in both Arabidopsis and rice. Taken together, we reported a new class of miRNAs, called RC-miRNAs, which were generated from the antisense strands of the miRNA precursors. Based on the results, we speculated that the mature RC-miRNAs might have subtle regulatory activity through target cleavages, but might possess short interfering RNA-like activity by guiding sequence-specific DNA methylation.     

Tissue and Process Specific microRNA–mRNA Co-Expression in Mammalian Development and Malignancy

An association between enrichment and depletion of microRNA (miRNA) binding sites, 3′ UTR length, and mRNA expression has been demonstrated in various developing tissues and tissues from different mature organs; but functional, context-dependent miRNA regulations have yet to be elucidated. Towards that goal, we examined miRNA–mRNA interactions by measuring miRNA and mRNA in the same tissue during development and also in malignant conditions. We identified significant miRNA-mediated biological process categories in developing mouse cerebellum and lung using non-targeted mRNA expression as the negative control. Although miRNAs in general suppress target mRNA messages, many predicted miRNA targets demonstrate a significantly higher level of co-expression than non-target genes in developing cerebellum. This phenomenon is tissue specific since it is not observed in developing lungs. Comparison of mouse cerebellar development and medulloblastoma demonstrates a shared miRNA–mRNA co-expression program for brain-specific neurologic processes such as synaptic transmission and exocytosis, in which miRNA target expression increases with the accumulation of multiple miRNAs in developing cerebellum and decreases with the loss of these miRNAs in brain tumors. These findings demonstrate the context-dependence of miRNA–mRNA co-expression.     

miRDeepFinder: a miRNA analysis tool for deep sequencing of plant small RNAs

miRDeepFinder is a software package developed to identify and functionally analyze plant microRNAs (miRNAs) and their targets from small RNA datasets obtained from deep sequencing. The functions available in miRDeepFinder include pre-processing of raw data, identifying conserved miRNAs, mining and classifying novel miRNAs, miRNA expression profiling, predicting miRNA targets, and gene pathway and gene network analysis involving miRNAs. The fundamental design of miRDeepFinder is based on miRNA biogenesis, miRNA-mediated gene regulation and target recognition, such as perfect or near perfect hairpin structures, different read abundances of miRNA and miRNA*, and targeting patterns of plant miRNAs. To test the accuracy and robustness of miRDeepFinder, we analyzed a small RNA deep sequencing dataset of Arabidopsis thaliana published in the GEO database of NCBI. Our test retrieved 128 of 131 (97.7%) known miRNAs that have a more than 3 read count in Arabidopsis. Because many known miRNAs are not associated with miRNA*s in small RNA datasets, miRDeepFinder was also designed to recover miRNA candidates without the presence of miRNA*. To mine as many miRNAs as possible, miRDeepFinder allows users to compare mature miRNAs and their miRNA*s with other small RNA datasets from the same species. Cleaveland software package was also incorporated into miRDeepFinder for miRNA target identification using degradome sequencing analysis. Using this new computational tool, we identified 13 novel miRNA candidates with miRNA*s from Arabidopsis and validated 12 of them experimentally. Interestingly, of the 12 verified novel miRNAs, a miRNA named AC1 spans the exons of two genes (UTG71C4 and UGT71C3). Both the mature AC1 miRNA and its miRNA* were also found in four other small RNA datasets. We also developed a tool, ??miRNA primer designer?? to design primers for any type of miRNAs. miRDeepFinder provides a powerful tool for analyzing small RNA datasets from all species, with or without the availability of genome information. miRDeepFinder and miRNA primer designer are freely available at http://www.leonxie.com/DeepFinder.php and at http://www.leonxie.com/miRNAprimerDesigner.php, respectively. A program (called RefFinder: http://www.leonxie.com/referencegene.php) was also developed for assessing the reliable reference genes for gene expression analysis, including miRNAs.     

Identification of Nitrogen Starvation-Responsive MicroRNAs in Arabidopsis thaliana

microRNAs (miRNAs) are a class of negative regulators that take part in many processes such as growth and development, stress responses, and metabolism in plants. Recently, miRNAs were shown to function in plant nutrient metabolism. Moreover, several miRNAs were identified in the response to nitrogen (N) deficiency. To investigate the functions of other miRNAs in N deficiency, deep sequencing technology was used to detect the expression of small RNAs under N-sufficient and -deficient conditions. The results showed that members from the same miRNA families displayed differential expression in response to N deficiency. Upon N starvation, the expression of miR169, miR171, miR395, miR397, miR398, miR399, miR408, miR827, and miR857 was repressed, whereas those of miR160, miR780, miR826, miR842, and miR846 were induced. miR826, a newly identified N-starvation-induced miRNA, was found to target the AOP2 gene. Among these N-starvation-responsive miRNAs, several were involved in cross-talk among responses to different nutrient (N, P, S, Cu) deficiencies. miR160, miR167, and miR171 could be responsible for the development of Arabidopsis root systems under N-starvation conditions. In addition, twenty novel miRNAs were identified and nine of them were significantly responsive to N-starvation. This study represents comprehensive expression profiling of N-starvation-responsive miRNAs and advances our understanding of the regulation of N homeostasis mediated by miRNAs.     

Computational identification of conserved microRNAs and their putative targets in the Hypericum perforatum L. flower transcriptome

MicroRNAs (miRNAs) have recently emerged as important regulators of gene expression in plants. Many miRNA families and their targets have been extensively studied in model species and major crops. We have characterized mature miRNAs along with their precursors and potential targets in Hypericum to generate a comprehensive list of conserved miRNA families and to investigate the regulatory role of selected miRNAs in biological processes that occur in the flower. St. John’s wort (Hypericum perforatum L., 2n = 4x = 32), a medicinal plant that produces pharmaceutically important metabolites with therapeutic activities, was chosen because it is regarded as an attractive model system for the study of apomixis. A computational in silico prediction of structure, in combination with an in vitro validation, allowed us to identify 7 pre-miRNAs, including miR156, miR166, miR390, miR394, miR396, and miR414. We demonstrated that H. perforatum flowers share highly conserved miRNAs and that these miRNAs potentially target dozens of genes with a wide range of molecular functions, including metabolism, response to stress, flower development, and plant reproduction. Our analysis paves the way toward identifying flower-specific miRNAs that may differentiate the sexual and apomictic reproductive pathways.