Comprehensive analysis of small RNA-seq data reveals that combination of miRNA with its isomiRs increase the accuracy of target prediction in Arabidopsis thaliana

The intricacies of isomiRs: from classification to clinical relevance

MicroRNAs (miRNAs) and isoforms of their archetype, called isomiRs, regulate gene expression via complementary base-pair binding to messenger RNAs (mRNAs). The partially evolutionarily conserved isomiR sequence variations are differentially expressed among tissues, populations, and genders, and between healthy and diseased states. Aiming towards the clinical use of isomiRs as diagnostic biomarkers and for therapeutic purposes, several challenges need to be addressed, including (i) clarification of isomiR definition, (ii) improved annotation in databases with new standardization (such as the mirGFF3 format), and (iii) improved methods of isomiR detection, functional verification, and in silico analysis. In this review we discuss the respective challenges, and highlight the opportunities for clinical use of isomiRs, especially in the light of increasing amounts of next-generation sequencing (NGS) data.

One locus, several functional RNAs-emerging roles of the mechanisms responsible for the sequence variability of microRNAs

With the development of modern molecular genetics, the original "one gene-one enzyme" hypothesis has been outdated. For protein coding genes, the discovery of alternative splicing and RNA editing provided the biochemical background for the RNA repertoire of a single locus, which also serves as an important pillar for the enormous protein variability of the genomes. Non-protein coding RNA genes were also revealed to produce several RNA species with distinct functions. The loci of microRNAs (miRNAs), encoding for small endogenous regulatory RNAs, were also found to produce a population of small RNAs, rather than a single defined product. This review aims to present the mechanisms contributing to the astonishing variability of miRNAs revealed by the new sequencing technologies. One important source is the careful balance of arm selection, producing sequentially different 5p- or 3p-miRNAs from the same pre-miRNA, thereby broadening the number of regulated target RNAs and the phenotypic response. In addition, the formation of 5', 3' and polymorphic isomiRs, with variable end and internal sequences also leads to a higher number of targeted sequences, and increases the regulatory output. These miRNA maturation processes, together with other known mechanisms such as RNA editing, further increase the potential outcome of this small RNA pathway. By discussing the subtle mechanisms behind the sequence diversity of miRNAs, this review intends to reveal this engaging aspect of the inherited "RNA world", how it contributes to the almost infinite molecular variability among living organisms, and how this variability can be exploited to treat human diseases.

Stem and leaf rust-induced miRNAome in bread wheat near-isogenic lines and their comparative analysis

Wheat rusts remain a major threat to global wheat production and food security. The R-gene-mediated resistance has been employed as an efficient approach to develop rust-resistant varieties. However, evolution of new fungal races and infection strategies put forward the urgency of unravelling novel molecular players, including non-coding RNAs for plant response. This study identified microRNAs associated with Sr36 and Lr45 disease resistance genes in response to stem and leaf rust, respectively. Here, small RNA sequencing was performed on susceptible and resistant wheat near-isogenic lines inoculated with stem and leaf rust pathotypes. microRNA mining in stem rust-inoculated cultivars revealed a total of distinct 26 known and 7 novel miRNAs, and leaf rust libraries culminated with 22 known and 4 novel miRNAs. The comparative analysis between two disease sets provides a better understanding of altered miRNA profiles associated with respective R-genes and infections. Temporal differential expression pattern of miRNAs pinpoints their role during the progress of infection. Differential expression pattern of miRNAs among various treatments as well as time-course expression of miRNAs revealed stem and leaf rust-responsive miRNAs and their possible role in balancing disease resistance/susceptibility. Disclosure of guide strand, passenger strand and a variant of novel-Tae-miR02 from different subgenome origins might serve as a potential link between stem and leaf rust defence mechanisms downstream to respective R-genes. The outcome from the analysis of microRNA dynamics among two rust diseases and further characterization of identified microRNAs can contribute to significant novel insights on wheat-rust interactions and rust management. KEY POINTS: • Identification and comparative analysis of stem and leaf rust-responsive miRNAs. • Chromosomal location and functional prediction of miRNAs. • Time-course expression analysis of pathogen-responsive miRNAs.

Identification of Small RNAs During High Light Acclimation in Arabidopsis thaliana

The biological significance of non-coding RNAs (ncRNAs) has been firmly established to be important for the regulation of genes involved in stress acclimation. Light plays an important role for the growth of plants providing the energy for photosynthesis; however, excessive light conditions can also cause substantial defects. Small RNAs (sRNAs) are a class of non-coding RNAs that regulate transcript levels of protein-coding genes and mediate epigenetic silencing. Next generation sequencing facilitates the identification of small non-coding RNA classes such as miRNAs (microRNAs) and small-interfering RNAs (siRNAs), and long non-coding RNAs (lncRNAs), but changes in the ncRNA transcriptome in response to high light are poorly understood. We subjected Arabidopsis plants to high light conditions and performed a temporal in-depth study of the transcriptome data after 3 h, 6 h, and 2 days of high light treatment. We identified a large number of high light responsive miRNAs and sRNAs derived from NAT gene pairs, lncRNAs and TAS transcripts. We performed target predictions for differentially expressed miRNAs and correlated their expression levels through mRNA sequencing data. GO analysis of the targets revealed an overrepresentation of genes involved in transcriptional regulation. In A. thaliana, sRNA-mediated regulation of gene expression in response to high light treatment is mainly carried out by miRNAs and sRNAs derived from NAT gene pairs, and from lncRNAs. This study provides a deeper understanding of sRNA-dependent regulatory networks in high light acclimation.

A Network-Based Analysis Reveals the Mechanism Underlying Vitamin D in Suppressing Cytokine Storm and Virus in SARS-CoV-2 Infection

Background

SARS-CoV-2 causes ongoing pandemic coronavirus disease of 2019 (COVID-19), infects the cells of the lower respiratory tract that leads to a cytokine storm in a significant number of patients resulting in severe pneumonia, shortness of breathing, respiratory and organ failure. Extensive studies suggested the role of Vitamin D in suppressing cytokine storm in COVID-19 and reducing viral infection; however, the precise molecular mechanism is not clearly known. In this work, bioinformatics and systems biology approaches were used to understand SARS-CoV-2 induced cytokine pathways and the potential mechanism of Vitamin D in suppressing cytokine storm and enhancing antiviral response.

Results

This study used transcriptome data and identified 108 differentially expressed host genes (DEHGs) in SARS-CoV-2 infected normal human bronchial epithelial (NHBE) cells compared to control. Then, the DEHGs was integrated with the human protein-protein interaction data to generate a SARS-CoV-2 induced host gene regulatory network (SiHgrn). Analysis of SiHgrn identified a sub-network "Cluster 1" with the highest MCODE score, 31 up-regulated genes, and predominantly associated immune and inflammatory response. Interestingly, the iRegulone tool identified that "Cluster 1" is under the regulation of transcription factors STAT1, STAT2, STAT3, POU2F2, and NFkB1, collectively referred to as "host response signature network". Functional enrichment analysis with NDEx revealed that the "host response signature network" is predominantly associated with critical pathways, including "cytokines and inflammatory response", "non-genomic action of Vitamin D", "the human immune response to tuberculosis", and "lung fibrosis". Finally, in-depth analysis and literature mining revealed that Vitamin D binds with its receptor and could work through two different pathways: (i) it inhibits the expression of pro-inflammatory cytokines through blocking the TNF induced NFkB1 signaling pathway; and (ii) it initiates the expression of interferon-stimulating genes (ISGs) for antiviral defense program through activating the IFN-α induced Jak-STAT signaling pathway.

Conclusion

This comprehensive study identified the pathways associated with cytokine storm in SARS-CoV-2 infection. The proposed underlying mechanism of Vitamin D could be promising in suppressing the cytokine storm and inducing a robust antiviral response in severe COVID-19 patients. The finding in this study urgently needs further experimental validations for the suitability of Vitamin D in combination with IFN-α to control severe COVID-19.

pssRNAit: A Web Server for Designing Effective and Specific Plant siRNAs with Genome-Wide Off-Target Assessment

We report an advanced web server, the plant-specific small noncoding RNA interference tool pssRNAit, which can be used to design a pool of small interfering RNAs (siRNAs) for highly effective, specific, and nontoxic gene silencing in plants. In developing this tool, we integrated the transcript dataset of plants, several rules governing gene silencing, and a series of computational models of the biological mechanism of the RNA interference (RNAi) pathway. The designed pool of siRNAs can be used to construct a long double-strand RNA and expressed through virus-induced gene silencing (VIGS) or synthetic transacting siRNA vectors for gene silencing. We demonstrated the performance of pssRNAit by designing and expressing the VIGS constructs to silence Phytoene desaturase (PDS) or a ribosomal protein-encoding gene, RPL10 (QM), in Nicotiana benthamiana We analyzed the expression levels of predicted intended-target and off-target genes using reverse transcription quantitative PCR. We further conducted an RNA-sequencing-based transcriptome analysis to assess genome-wide off-target gene silencing triggered by the fragments that were designed by pssRNAit, targeting different homologous regions of the PDS gene. Our analyses confirmed the high accuracy of siRNA constructs designed using pssRNAit The pssRNAit server, freely available at https://plantgrn.noble.org/pssRNAit/, supports the design of highly effective and specific RNAi, VIGS, or synthetic transacting siRNA constructs for high-throughput functional genomics and trait improvement in >160 plant species.

Investigating the pathogenic SNPs in BLM helicase and their biological consequences by computational approach

The BLM helicase protein plays a vital role in DNA replication and the maintenance of genomic integrity. Variation in the BLM helicase gene resulted in defects in the DNA repair mechanism and was reported to be associated with Bloom syndrome (BS) and cancer. Despite extensive investigation of helicase proteins in humans, no attempt has previously been made to comprehensively analyse the single nucleotide polymorphism (SNPs) of the BLM gene. In this study, a comprehensive analysis of SNPs on the BLM gene was performed to identify, characterize and validate the pathogenic SNPs using computational approaches. We obtained SNP data from the dbSNP database version 150 and mapped these data to the genomic coordinates of the "NM_000057.3" transcript expressing BLM helicase (P54132). There were 607 SNPs mapped to missense, 29 SNPs mapped to nonsense, and 19 SNPs mapped to 3'-UTR regions. Initially, we used many consensus tools of SIFT, PROVEAN, Condel, and PolyPhen-2, which together increased the accuracy of prediction and identified 18 highly pathogenic non-synonymous SNPs (nsSNPs) out of 607 SNPs. Subsequently, these 18 high-confidence pathogenic nsSNPs were analysed for BLM protein stability, structure-function relationships and disease associations using various bioinformatics tools. These 18 mutants of the BLM protein along with the native protein were further investigated using molecular dynamics simulations to examine the structural consequences of the mutations, which might reveal their malfunction and contribution to disease. In addition, 28 SNPs were predicted as "stop gained" nonsense SNPs and one SNP was predicted as "start lost". Two SNPs in the 3'UTR were found to abolish miRNA binding and thus may enhance the expression of BLM. Interestingly, we found that BLM mRNA overexpression is associated with different types of cancers. Further investigation showed that the dysregulation of BLM is associated with poor overall survival (OS) for lung and gastric cancer patients and hence led to the conclusion that BLM has the potential to be used as an important prognostic marker for the detection of lung and gastric cancer.

A conserved sequence signature is essential for robust plant miRNA biogenesis

Micro (mi)RNAs are 20–22nt long non-coding RNA molecules involved in post-transcriptional silencing of targets having high base-pair complementarity. Plant miRNAs are processed from long Pol II-transcripts with specific stem-loop structures by Dicer-like (DCL) 1 protein. Although there were reports indicating how a specific region is selected for miRNA biogenesis, molecular details were unclear. Here, we show that the presence of specific GC-rich sequence signature within miRNA/miRNA* region is required for the precise miRNA biogenesis. The involvement of GC-rich signatures in precise processing and abundance of miRNAs was confirmed through detailed molecular and functional analysis. Consistent with the presence of the miRNA-specific GC signature, target RNAs of miRNAs also possess conserved complementary sequence signatures in their miRNA binding motifs. The selection of these GC signatures was dependent on an RNA binding protein partner of DCL1 named HYL1. Finally, we demonstrate a direct application of this discovery for enhancing the abundance and efficiency of artificial miRNAs that are popular in plant functional genomic studies.

Post-transcriptional regulation of several biological processes involved in latex production in Hevea brasiliensis

Background

Small RNAs modulate plant gene expression at both the transcriptional and post-transcriptional level, mostly through the induction of either targeted DNA methylation or transcript cleavage, respectively. Small RNA networks are involved in specific plant developmental processes, in signaling pathways triggered by various abiotic stresses and in interactions between the plant and viral and non-viral pathogens. They are also involved in silencing maintenance of transposable elements and endogenous viral elements. Alteration in small RNA production in response to various environmental stresses can affect all the above-mentioned processes. In rubber trees, changes observed in small RNA populations in response to trees affected by tapping panel dryness, in comparison to healthy ones, suggest a shift from a transcriptional to a post-transcriptional regulatory pathway. This is the first attempt to characterise small RNAs involved in post-transcriptional silencing and their target transcripts in Hevea.

Methods

Genes producing microRNAs (MIR genes) and loci producing trans-activated small interfering RNA (ta-siRNA) were identified in the clone PB 260 re-sequenced genome. Degradome libraries were constructed with a pool of total RNA from six different Hevea tissues in stressed and non-stressed plants. The analysis of cleaved RNA data, associated with genomics and transcriptomics data, led to the identification of transcripts that are affected by 20–22 nt small RNA-mediated post-transcriptional regulation. A detailed analysis was carried out on gene families related to latex production and in response to growth regulators.

Results

Compared to other tissues, latex cells had a higher proportion of transcript cleavage activity mediated by miRNAs and ta-siRNAs. Post-transcriptional regulation was also observed at each step of the natural rubber biosynthesis pathway. Among the genes involved in the miRNA biogenesis pathway, our analyses showed that all of them are expressed in latex. Using phylogenetic analyses, we show that both the Argonaute and Dicer-like gene families recently underwent expansion. Overall, our study underlines the fact that important biological pathways, including hormonal signalling and rubber biosynthesis, are subject to post-transcriptional silencing in laticifers.

Systematic Analysis of Alkaline/Neutral Invertase Genes Reveals the Involvement of Smi-miR399 in Regulation of SmNINV3 and SmNINV4 in Salvia miltiorrhiza

Alkaline/neutral invertases (NINVs), which irreversibly catalyze the hydrolysis of sucrose into fructose and glucose, play crucial roles in carbohydrate metabolism and plant development. Comprehensive insights into NINV genes are lacking in Salvia miltiorrhiza, a well-known traditional Chinese medicinal (TCM) plant with significant medicinal and economic value. Through genome-wide prediction, nine putative SmNINV genes, termed SmNINV1-SmNINV9, were identified. Integrated analysis of gene structures, sequence features, conserved domains, conserved motifs and phylogenetic trees revealed the conservation and divergence of SmNINVs. The identified SmNINVs were differentially expressed in roots, stems, leaves, flowers, and different root tissues. They also responded to drought, salicylic acid, yeast extract, and methyl jasmonate treatments. More importantly, computational prediction and experimental validation showed that SmNINV3 and SmNINV4 were targets of Smi-miR399, a conserved miRNA previously shown to affect Pi uptake and translocation through the cleavage of PHOSPHATE2 (PHO2). Consistently, analysis of 43 NINV genes and 26 miR399 sequences from Arabidopsis thaliana, Populus trichocarpa, Manihot esculenta, and Solanum lycopersicum showed that various AtNINV, PtNINV, MeNINV, and SlNINV genes were regulated by miR399. It indicates that the miR399-NINV module exists widely in plants. Furthermore, Smi-miR399 also cleaved SmPHO2 transcripts in S. miltiorrhiza, suggesting the complexity of NINVs, PHO2, and miR399 networks.

Integrated Network Analysis Reveals FOXM1 and MYBL2 as Key Regulators of Cell Proliferation in Non-small Cell Lung Cancer

Background: Loss of control on cell division is an important factor for the development of non-small cell lung cancer (NSCLC), however, its molecular mechanism and gene regulatory network are not clearly understood. This study utilized the systems bioinformatics approach to reveal the “driver-network” involve in tumorigenic processes in NSCLC.

Methods: A meta-analysis of gene expression data of NSCLC was integrated with protein-protein interaction (PPI) data to construct an NSCLC network. MCODE and iRegulone were used to identify the local clusters and its upstream transcription regulators involve in NSCLC. Pair-wise gene expression correlation was performed using GEPIA. The survival analysis was performed by the Kaplan-Meier plot.

Results: This study identified a local “driver-network” with highest MCODE score having 26 up-regulated genes involved in the process of cell proliferation in NSCLC. Interestingly, the “driver-network” is under the regulation of TFs FOXM1 and MYBL2 as well as miRNAs. Furthermore, the overexpression of member genes in “driver-network” and the TFs are associated with poor overall survival (OS) in NSCLC patients.

Conclusion: This study identified a local “driver-network” and its upstream regulators responsible for the cell proliferation in NSCLC, which could be promising biomarkers and therapeutic targets for NSCLC treatment.

Identifying High Confidence microRNAs in the Developing Seeds of Jatropha curcas

MicroRNAs (miRNAs) are endogenously short noncoding regulatory RNAs implicated in plant development and physiology. Nine small RNA (sRNA) libraries from three typical seed developmental stages (young, intermediate, and mature) were generated by deep sequencing to identify the miRNAs of J. curcas, a potential oilseed crop for the production of renewable oil. Strict criteria were adopted to identify 93 high confidence miRNAs including 48 conserved miRNAs and 45 novel miRNAs. Target genes of these miRNAs were involved in a broad range of physiological functions, including gene expression regulation, primary & secondary metabolism, growth & development, signal transduction, and stress response. About one third (29 out of 93) miRNAs showed significant changes in expression levels during the seed developmental process, indicating that the miRNAs might regulate its targets by their changes of transcription levels in seed development. However, most miRNAs were found differentially expressed in the late stage of seed development, suggesting that miRNAs play more important roles in the stage when seed accumulating organic matters and suffering dehydration stress. This study presents the first large scale identification of high confidence miRNAs in the developing seeds of J. curcas.

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.

Modulation of miRNA expression in natural populations of A. thaliana along a wide altitudinal gradient of Indian Himalayas

Plant populations growing along an altitudinal gradient are exposed to different environmental conditions. They are excellent resources to study regulatory mechanisms adopted by plants to respond to different environmental stresses. Regulation by miRNA is one of such strategies. Here, we report how different miRNAs are preferentially expressed in the three natural populations of A. thaliana originating from a wide altitudinal range. The expression level of miRNAs was mostly governed by temperature and radiation. Majority of the identified miRNAs expressed commonly in the three populations. However, 30 miRNAs expressed significantly at different level between the low and the high altitude populations. Most of these miRNAs regulate the genes associated with different developmental processes, abiotic stresses including UV, cold, secondary metabolites, etc. Further, the expression of miR397 and miR858 involved in lignin biosynthesis and regulation of secondary metabolites respectively, may be regulated by light intensity. A few miRNAs expressed at increasing level with the increase in the altitude of the site indicating environment driven tight regulation of these miRNAs. Further, several novel miRNAs and isomiR diversity specific to the Himalayas are reported which might have an adaptive advantage. To the best of our knowledge, this is the first report on miRNA expression from natural plant populations.

An Overview of miRNA and miRNA Target Analysis Tools

microRNA molecules have been shown to play various significant roles in many physiological and pathophysiological processes in living organisms. The tremendous interest in these molecules has led to the significant development and constant release of a number of computational tools useful for basic as well as advanced miRNA-related analyses. These approaches have various constantly evolving utilities, such as detection, target prediction, functional annotation, and many others. In this chapter, we provide an overview of several computational tools useful for broadly defined plant miRNA analysis.

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

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

Identification of the miRNAome of early mesoderm progenitor cells and cardiomyocytes derived from human pluripotent stem cells

MicroRNAs are small non-coding RNAs involved in post-transcriptional regulation of gene expression related to many cellular functions. We performed a small-RNAseq analysis of cardiac differentiation from pluripotent stem cells. Our analyses identified some new aspects about microRNA expression in this differentiation process. First, we described a dynamic expression profile of microRNAs where some of them are clustered according to their expression level. Second, we described the extensive network of isomiRs and ADAR modifications. Third, we identified the microRNAs families and clusters involved in the establishment of cardiac lineage and define the mirRNAome based on these groups. Finally, we were able to determine a more accurate miRNAome associated with cardiomyocytes by comparing the expressed microRNAs with other mature cells. MicroRNAs exert their effect in a complex and interconnected way, making necessary a global analysis to better understand their role. Our data expands the knowledge of microRNAs and their implications in cardiomyogenesis.

The miRNAome dynamics during developmental and metabolic reprogramming of tomato root infected with potato cyst nematode

Cyst-forming plant-parasitic nematodes are pests threatening many crops. By means of their secretions cyst nematodes induce the developmental and metabolic reprogramming of host cells that lead to the formation of a syncytium, which is the sole food source for growing nematodes. The in depth micro RNA (miRNA) dynamics in the syncytia induced by Globodera rostochiensis in tomato roots was studied. The miRNAomes were obtained from syncytia covering the early and intermediate developmental stages, and were the subject of differential expression analysis. The expression of 1235 miRNAs was monitored. The fold change (log₂FC) ranged from -7.36 to 8.38, indicating that this transcriptome fraction was very variable. Moreover, we showed that the DE (differentially expressed) miRNAs do not fully overlap between the selected time points, suggesting infection stage specific regulation by miRNA. The correctness of RNA-seq expression profiling was confirmed by qRT-PCR (quantitative Real Time Polymerase Chain Reaction) for seven miRNA species. Down- and up-regulated miRNA species, including their isomiRs, were further used to identify their potential targets. Among them there are a large number of transcription factors linked to different aspects of plant development belonging to gene families, such as APETALA2 (AP2), SQUAMOSA (MADS-box), MYB, GRAS, and AUXIN RESPONSE FACTOR (ARF). The substantial portion of potential target genes belong to the NB-LRR and RLK (RECEPTOR-LIKE KINASE) families, indicating the involvement of miRNA mediated regulation in defense responses. We also collected the evidence for target cleavage in the case of 29 miRNAs using one of three alternative methods: 5' RACE (5' Rapid Amplification of cDNA Ends), a search of tasiRNA within our datasets, and the meta-analysis of tomato degradomes in the GEO (Gene Expression Omnibus) database. Eight target transcripts showed a negative correlation with their respective miRNAs at two or three time points. These results indicate a large regulatory potential for miRNAs in tuning the development and defense responses.

High-throughput sequencing analysis revealed the regulation patterns of small RNAs on the development of A. comosus var. bracteatus leaves

Studies of the molecular mechanisms involved in the formation of the albino leaf cells are important for understanding the development of chimera leaves in Ananas comosus var. bracteatus. In this study, we identified a total of 163 novel miRNAs involved in the development of complete white (CWh) and complete green (CGr) leaves using high-throughput sequencing method. The potential miRNA target genes were predicted and annotated using the NR, Swiss-Prot, GO, COG, KEGG, KOG and Pfam databases. The main biological processes regulated by miRNAs were revealed. The miRNAs which potentially play important roles in the development of the leaves and the albino of the CWh leaf cells were selected and their expression patterns were analyzed. The expression levels of nine miRNAs and their potential target genes were studied using qRT-PCR. These results will help to elucidate the functional and regulatory roles of miRNAs in the formation of the albino cells and the development of the leaves of A. comosus var. bracteatus. These data may also be helpful as a resource for studies of small RNA in other leaf color chimeric plant species.

A Role for the F-Box Protein HAWAIIAN SKIRT in Plant microRNA Function

As regulators of gene expression in multicellular organisms, microRNAs (miRNAs) are crucial for growth and development. Although a plethora of factors involved in their biogenesis and action in Arabidopsis (Arabidopsis thaliana) has been described, these processes and their fine-tuning are not fully understood. Here, we used plants expressing an artificial miRNA target mimic (MIM) to screen for negative regulators of miR156. We identified a new mutant allele of the F-box gene HAWAIIAN SKIRT (HWS; At3G61590), hws-5, as a suppressor of the MIM156-induced developmental and molecular phenotypes. In hws plants, levels of some endogenous miRNAs are increased and their mRNA targets decreased. Plants constitutively expressing full-length HWS-but not a truncated version lacking the F-box domain-display morphological and molecular phenotypes resembling those of mutants defective in miRNA biogenesis and activity. In combination with such mutants, hws loses its delayed floral organ abscission ("skirt") phenotype, suggesting epistasis. Also, the hws transcriptome profile partially resembles those of well-known miRNA mutants hyl1-2, se-3, and ago1-27, pointing to a role in a common pathway. We thus propose HWS as a novel, F-box dependent factor involved in miRNA function.

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.

Naturally existing isoforms of miR-222 have distinct functions

Deep-sequencing reveals extensive variation in the sequence of endogenously expressed microRNAs (termed ‘isomiRs’) in human cell lines and tissues, especially in relation to the 3′ end. From the immunoprecipitation of the microRNA-binding protein Argonaute and the sequencing of associated small RNAs, we observe extensive 3′-isomiR variation, including for miR-222 where the majority of endogenously expressed miR-222 is extended by 1–5 nt compared to the canonical sequence. We demonstrate this 3′ heterogeneity has dramatic implications for the phenotype of miR-222 transfected cells, with longer isoforms promoting apoptosis in a size (but not 3′ sequence)-dependent manner. The transfection of longer miR-222 isomiRs did not induce an interferon response, but did downregulate the expression of many components of the pro-survival PI3K-AKT pathway including PIK3R3, a regulatory subunit whose knockdown phenocopied the expression of longer 222 isoforms in terms of apoptosis and the inhibition of other PI3K-AKT genes. As this work demonstrates the capacity for 3′ isomiRs to mediate differential functions, we contend more attention needs to be given to 3′ variance given the prevalence of this class of isomiR.

Two-tailed RT-qPCR: a novel method for highly accurate miRNA quantification

MicroRNAs are a class of small non-coding RNAs that serve as important regulators of gene expression at the posttranscriptional level. They are stable in body fluids and pose great potential to serve as biomarkers. Here, we present a highly specific, sensitive and cost-effective system to quantify miRNA expression based on two-step RT-qPCR with SYBR-green detection chemistry called Two-tailed RT-qPCR. It takes advantage of novel, target-specific primers for reverse transcription composed of two hemiprobes complementary to two different parts of the targeted miRNA, connected by a hairpin structure. The introduction of a second probe ensures high sensitivity and enables discrimination of highly homologous miRNAs irrespectively of the position of the mismatched nucleotide. Two-tailed RT-qPCR has a dynamic range of seven logs and a sensitivity sufficient to detect down to ten target miRNA molecules. It is capable to capture the full isomiR repertoire, leading to accurate representation of the complete miRNA content in a sample. The reverse transcription step can be multiplexed and the miRNA profiles measured with Two-tailed RT-qPCR show excellent correlation with the industry standard TaqMan miRNA assays (r² = 0.985). Moreover, Two-tailed RT-qPCR allows for rapid testing with a total analysis time of less than 2.5 hours.

Using RNA sequencing to identify putative competing endogenous RNAs (ceRNAs) potentially regulating fat metabolism in bovine liver

RNA sequencing has been extensively used to study specific gene expression patterns to discover potential key genes related to complex traits of interest in animals. Of note, a new regulatory mechanism builds a large-scale regulatory network among transcriptome, where lncRNAs act as competing endogenous RNAs (ceRNAs) to sponge miRNAs to regulate the expression of miRNA target genes post-transcriptionally. In this study, we sequenced the cDNA and sRNA libraries of nine liver samples from three Holstein cows during dry period, early lactation, and peak of lactation with HiSeq platform. As a result, we identified 665 genes, 57 miRNAs and 33 lncRNAs that displayed differential expression patterns across periods. Subsequently, a total of 41ceRNA pairs (lncRNA-mRNA) sharing 11 miRNAs were constructed including 30 differentially expressed genes. Importantly, 12 among them were presented in our large metabolic networks, and predicted to influence the lipid metabolism through insulin, PI3K-Akt, MAPK, AMPK, mTOR, and PPAR signaling pathways, thus, these genes were considered as the most promising candidates for milk fat formation. To our knowledge, this is first investigation to profile the ceRNA regulatory networks of liver transcriptome that could affect milk fat synthesis in bovine, providing a new view of the regulatory mechanism of RNAs.

Small RNA sequencing reveals a role for sugarcane miRNAs and their targets in response to Sporisorium scitamineum infection

BACKGROUND

Sugarcane smut caused by Sporisorium scitamineum leads to a significant reduction in cane yield and sucrose content. MicroRNAs (miRNAs) play an important role in regulating plant responses to biotic stress. The present study was the first to use two sugarcane genotypes, YA05-179 (smut-resistant) and ROC22 (smut-susceptible), to identify differentially expressed miRNAs in sugarcane challenged with S. scitamineum by using high-throughput sequencing.

RESULTS

The predicted target gene number corresponding to known differentially expressed miRNAs in YA05-179 was less than that in ROC22, however most of them were in common. Expression of differential miRNAs under S. scitamineum challenge was mostly downregulated, with similar trends in the two varieties. Gene ontology (GO) analysis showed that the target gene classification of known miRNAs was similar to that of the newly identified miRNAs. These were mainly associated with cellular processes and metabolic processes in the biological process category, as well as combination and catalytic activity in the molecular function category. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis revealed that these predicted target genes involved in a series of physiological and biochemical pathways or disease resistance-related physiological metabolism and signal transduction pathways, suggesting that the molecular interaction mechanism between sugarcane and S. scitamineum was a complex network system. These findings also showed certain predicted target genes of miR5671, miR5054, miR5783, miR5221, and miR6478 play roles in the mitogen-activated protein kinase (MAPK) signaling pathway, plant hormone signal transduction, and plant-pathogen interaction. Quantitative real-time PCR (qRT-PCR) analysis showed that majority of the known miRNAs and its predicted target genes followed a negatively regulated mode. Seven out of eight predicted target genes showed identical expression after 12 h treatment and reached the highest degree of matching at 48 h, indicating that the regulatory role of miRNAs on the target genes in sugarcane was maximized at 48 h after S. scitamineum challenge.

CONCLUSIONS

Taken together, our findings serve as evidence for the association of miRNA expression with the molecular mechanism underlying the pathogenesis of sugarcane smut, particularly on the significance of miRNA levels in relation to the cultivation of smut-resistant sugarcane varieties.