A genome-wide survey of microRNA truncation and 3' nucleotide addition events in larch (Larix leptolepis)

Diff isomiRs: Large-scale detection of differential isomiRs for understanding non-coding regulated stress omics in plants

Plants have an amazing ability to cope with wide variety of stresses by regulating the expression of genes and thus by altering the physiological status. In the past few years, canonical microRNA variants (isomiRs) have been shown to play pivotal roles by acting as regulators of the transcriptional machinery. In the present research, we present Diff isomiRs, a web-based exploratory repository of differential isomiRs across 16 sequenced plant species representing a total of 433 datasets across 21 different stresses and 158 experimental states. Diff isomiRs provides the high-throughput detection of differential isomiRs using mapping-based and model-based differential analysis revealing a total of 16,157 and 2,028 differential isomiRs, respectively. Easy-to-use and web-based exploration of differential isomiRs provides several features such as browsing of the differential isomiRs according to stress or species, as well as association of the differential isomiRs to targets and plant endogenous target mimics (PeTMs). Diff isomiRs also provides the relationship between the canonical miRNAs, isomiRs and the miRNA-target interactions. This is the first web-based large-scale repository for browsing differential isomiRs and will facilitate better understanding of the regulatory role of the isomiRs with respect to the canonical microRNAs. Diff isomiRs can be accessed at: www.mcr.org.in/diffisomirs.

Genome-wide identification and characterization of novel microRNAs in seed development of soybean

MicroRNAs (miRNAs) are important and ubiquitous regulators of gene expression in eukaryotes. However, the information about miRNAs population and their regulatory functions involving in soybean seed development remains incomplete. Base on the Dicer-like1-mediated cleavage signals during miRNA processing could be employed for novel miRNA discovery, a genome-wide search for miRNA candidates involved in seed development was carried out. As a result, 17 novel miRNAs, 14 isoforms of miRNA (isomiRs) and 31 previously validated miRNAs were discovered. These novel miRNAs and isomiRs represented tissue-specific expression and the isomiRs showed significantly higher abundance than that of their miRNA counterparts in different tissues. After target prediction and degradome sequencing data-based validation, 13 novel miRNA-target pairs were further identified. Besides, five targets of 22-nt iso-gma-miR393h were found to be triggered to produce secondary trans-acting siRNA (ta-siRNAs). Summarily, our results could expand the repertoire of miRNAs with potentially important functions in soybean.

Circulating miRNAs, isomiRs and small RNA clusters in human plasma and breast milk

Circulating small RNAs, including miRNAs but also isomiRs and other RNA species, have the potential to be used as non-invasive biomarkers for communicable and non-communicable diseases. This study aims to characterize and compare small RNA profiles in human biofluids. For this purpose, RNA was extracted from plasma and breast milk samples from 15 healthy postpartum mothers. Small RNA libraries were prepared with the NEBNext® small RNA library preparation kit and sequenced in an Illumina HiSeq2000 platform. miRNAs, isomiRs and clusters of small RNAs were annotated using seqBuster/seqCluster framework in 5 plasma and 10 milk samples that passed the initial quality control. The RNA yield was 81 ng/mL [standard deviation (SD): 41] and 3985 ng/mL (SD: 3767) for plasma and breast milk, respectively. Mean number of good quality reads was 4.04 million (M) (40.01% of the reads) in plasma and 12.5M (89.6%) in breast milk. One thousand one hundred eighty two miRNAs, 12,084 isomiRs and 1,053 small RNA clusters that included piwi-interfering RNAs (piRNAs), tRNAs, small nucleolar RNAs (snoRNA) and small nuclear RNAs (snRNAs) were detected. Samples grouped by biofluid, with 308 miRNAs, 1,790 isomiRs and 778 small RNA clusters differentially detected. In summary, plasma and milk showed a different small RNA profile. In both, miRNAs, piRNAs, tRNAs, snRNAs, and snoRNAs were identified, confirming the presence of non-miRNA species in plasma, and describing them for the first time in milk.

Plant IsomiR Atlas: Large Scale Detection, Profiling, and Target Repertoire of IsomiRs in Plants

microRNAs (miRNAs) play an important role as key regulators controlling the post-transcriptional events in plants across development, abiotic and biotic stress, tissue polarity and also in defining the evolutionary basis of the origin of the post-transcriptional machinery. Identifying patterns of regulated and co-regulated small RNAs, in particular miRNAs and their sequence variants with the availability of next generation sequencing approaches has widely demonstrated the role of miRNAs and their temporal regulation in maintaining plant development and their response to stress conditions. Although the role of canonical miRNAs has been widely explored and functional diversity is revealed, those works for isomiRs are still limited and urgent to be carried out across plants. This relative lack of information with respect to isomiRs might be attributed to the non-availability of large-scale detection of isomiRs across wide plant species. In the present research, we addressed this by developing Plant isomiR Atlas, which provides large-scale detection of isomiRs across 23 plant species utilizing 677 smallRNAs datasets and reveals a total of 98,374 templated and non-templated isomiRs from 6,167 precursors. Plant isomiR Atlas provides several visualization features such as species specific isomiRs, isomiRs and canonical miRNAs overlap, terminal modification classifications, target identification using psRNATarget and TargetFinder and also canonical miRNAs:target interactions. Plant isomiR Atlas will play a key role in understanding the regulatory nature of miRNAome and will accelerate to understand the functional role of isomiRs. Plant isomiR Atlas is available at www.mcr.org.in/isomir. One Sentence Summary  Plant isomiR Atlas will play a key role in understanding the regulatory nature of miRNAome and will accelerate the understanding and diversity of functional targets of plants isomiRs.

Transcriptome Analysis of mRNA and miRNA in Somatic Embryos of Larix leptolepis Subjected to Hydrogen Treatment

Hydrogen is a therapeutic antioxidant that has been used extensively in clinical trials. It also acts as a bioactive molecule that can alleviate abiotic stress in plants. However, the biological effects of hydrogen in somatic embryos and the underlying molecular basis remain largely unknown. In this study, the morphological and physiological influence of exogenous H₂ treatment during somatic embryogenesis was characterized in Larix leptolepis Gordon. The results showed that exposure to hydrogen increased the proportions of active pro-embryogenic cells and normal somatic embryos. We sequenced mRNA and microRNA (miRNA) libraries to identify global transcriptome changes at different time points during H₂ treatment of larch pro-embryogenic masses (PEMs). A total of 45,393 mRNAs and 315 miRNAs were obtained. Among them, 4253 genes and 96 miRNAs were differentially expressed in the hydrogen-treated libraries compared with the control. Further, a large number of the differentially expressed mRNAs and miRNAs were related to reactive oxygen species (ROS) homeostasis and cell cycle regulation. We also identified 4399 potential target genes for 285 of the miRNAs. The differential expression data and the mRNA-miRNA interaction network described here provide new insights into the molecular mechanisms that determine the performance of PEMs exposed to H₂ during somatic embryogenesis.

3' and 5' microRNA-end post-biogenesis modifications in plant transcriptomes: Evidences from small RNA next generation sequencing data analysis

The processing of miRNA from its precursors is a precisely regulated process and after biogenesis, the miRNAs are amenable to different kinds of modifications by the addition or deletion of nucleotides at the terminal ends. However, the mechanism and functions of such modifications are not well studied in plants. In this study, we have specifically analysed the terminal end non-templated miRNA modifications, using NGS data of rice, tomato and Arabidopsis small RNA transcriptomes from different tissues and physiological conditions. Our analysis reveals template independent terminal end modifications in the mature as well as passenger strands of the miRNA duplex. Interestingly, it is also observed that miRNA sequences terminating with a cytosine (C) at the 3' end undergo a higher percentage of 5' end modifications. The terminal end modifications did not correlate with the miRNA abundances and are independent of tissue types, physiological conditions and plant species. Our analysis indicates that the addition of nucleotides at miRNA ends is not influenced by the absence of RNA dependent RNA polymerase 6. Moreover the terminal end modified miRNAs are also observed amongst AGO1 bound small RNAs and have potential to alter target, indicating its important functional role in repression of gene expression.

Stress responsive miRNAs and isomiRs in cereals

Abiotic and biotic stress conditions are vital determinants in the production of cereals, the major caloric source in human nutrition. Small RNAs, miRNAs and isomiRs are central to post-transcriptional regulation of gene expression in a variety of cellular processes including development and stress responses. Several miRNAs have been identified using new technologies and have roles in stress responses in plants, including cereals. The overall knowledge about the cereal miRNA repertoire, as well as an understanding of complex miRNA mediated mechanisms of target regulation in response to stress conditions, is far from complete. Ongoing efforts that add to our understanding of complex miRNA machinery have implications in plant response to stress conditions. Additionally, sequence variants of miRNAs (isomiRNAs or isomiRs), regulation of their expression through dissection of upstream regulatory elements, the role of Processing-bodies (P-bodies) in miRNA exerted gene regulation and yet unveiled organellar plant miRNAs are newly emerging topics, which will contribute to the elucidation of the miRNA machinery and its role in cereal tolerance against abiotic and biotic stresses.

Evolutionary and expression analysis of miR-#-5p and miR-#-3p at the miRNAs/isomiRs levels

We mainly discussed miR-#-5p and miR-#-3p under three aspects: (1) primary evolutionary analysis of human miRNAs; (2) evolutionary analysis of miRNAs from different arms across the typical 10 vertebrates; (3) expression pattern analysis of miRNAs at the miRNA/isomiR levels using public small RNA sequencing datasets. We found that no bias can be detected between the numbers of 5p-miRNA and 3p-miRNA, while miRNAs from miR-#-5p and miR-#-3p show variable nucleotide compositions. IsomiR expression profiles from the two arms are always stable, but isomiR expressions in diseased samples are prone to show larger degree of dispersion. miR-#-5p and miR-#-3p have relative independent evolution/expression patterns and datasets of target mRNAs, which might also contribute to the phenomena of arm selection and/or arm switching. Simultaneously, miRNA/isomiR expression profiles may be regulated via arm selection and/or arm switching, and the dynamic miRNAome and isomiRome will adapt to functional and/or evolutionary pressures. A comprehensive analysis and further experimental study at the miRNA/isomiR levels are quite necessary for miRNA study.

Global microRNA modification in cotton (Gossypium hirsutum L.)

MicroRNAs (miRNAs) are small noncoding RNAs participating in versatile biological processes via post-transcriptionally gene regulation. However, how miRNAs are modified or degraded remains unknown, despite years of studies have unravelled much details of miRNA biogenesis and function. Here, we systematically investigated miRNA modification using six small RNA sequencing libraries generated from cotton seedling as well as cotton fibre at five developmental stages. Our results show that 1-2-nt truncation and addition on both 5' and 3' ends of miRNAs are the major modification forms. The 5' and 3' end miRNA modification was almost equal in the six development stages. Truncation was more common than addition on both 5' and 3' end. Structure analysis of the 5' and 3' ends of miRNAs and isomiRs shows that uridine is the preferential nucleotide at the first position of both 5' and 3' ends. According to analysis of nucleotides truncated and tailed from miRNAs, both miRNAs and isomiRs share a similar positional structure distribution at their 5' and 3' ends, respectively. Furthermore, opposite to previous reports, cytodine is more frequently truncated and tailed from the two ends of isomiRs, implying existence of a complex cytodine balance in isomiRs. Comparison of isomiR expression shows differential miRNA modification amongst the six developmental stages in terms of selective modification form, development-dependent modification and differential expression abundance. Our results globally uncovered miRNA modification features in cotton, which could contribute us to understanding miRNA's postmature modification and its regulatory function.

MicroRNA: a new target for improving plant tolerance to abiotic stress

MicroRNAs (miRNAs) are an extensive class of endogenous, small RNA molecules that sit at the heart of regulating gene expression in multiple developmental and signalling pathways. Recent studies have shown that abiotic stresses induce aberrant expression of many miRNAs, thus suggesting that miRNAs may be a new target for genetically improving plant tolerance to certain stresses. These studies have also shown that miRNAs respond to environmental stresses in a miRNA-, stress-, tissue-, and genotype-dependent manner. During abiotic stress, miRNAs function by regulating target genes within the miRNA-target gene network and by controlling signalling pathways and root development. Generally speaking, stress-induced miRNAs lead to down-regulation of negative regulators of stress tolerance whereas stress-inhibited miRNAs allow the accumulation and function of positive regulators. Currently, the majority of miRNA-based studies have focused on the identification of miRNAs that are responsive to different stress conditions and analysing their expression profile changes during these treatments. This has predominately been accomplished using deep sequencing technologies and other expression analyses, such as quantitative real-time PCR. In the future, more function and expression studies will be necessary in order to elucidate the common miRNA-mediated regulatory mechanisms that underlie tolerance to different abiotic stresses. The use of artificial miRNAs, as well as overexpression and knockout/down of both miRNAs and their targets, will be the best techniques for determining the specific roles of individual miRNAs in response to environmental stresses.

A challenge for miRNA: multiple isomiRs in miRNAomics

Accumulating evidence suggests that a single microRNA (miRNA) locus can generate a series of sequences during miRNA maturation process. These multiple sequences, called miRNA variants, or isomiRs, have different lengths and different 5' and 3' ends. Some of these isomiRs are detected as varied nucleotides and 3' additional non-template nucleotides. As physiological miRNA isoforms, they have drawn attention for possible regulatory biological roles. The present work mainly reviews miRNA/isomiR biogenesis, isomiR expression patterns, and functional and evolutionary implications, especially between isomiRs from homologous and clustered miRNA loci. The phenomenon of multiple isomiRs and their biological roles indicates that analysis performed at the miRNA and isomiR levels should be included in miRNA studies. This may enrich and complicate miRNA biogenesis and coding-non-coding RNA regulatory networks.