Evidence that miRNAs are different from other RNAs

Multi-Kernel Graph Attention Deep Autoencoder for MiRNA-Disease Association Prediction

Accumulating evidence indicates that microRNAs (miRNAs) can control and coordinate various biological processes. Consequently, abnormal expressions of miRNAs have been linked to various complex diseases. Recognizable proof of miRNA-disease associations (MDAs) will contribute to the diagnosis and treatment of human diseases. Nevertheless, traditional experimental verification of MDAs is laborious and limited to small-scale. Therefore, it is necessary to develop reliable and effective computational methods to predict novel MDAs. In this work, a multi-kernel graph attention deep autoencoder (MGADAE) method is proposed to predict potential MDAs. In detail, MGADAE first employs the multiple kernel learning (MKL) algorithm to construct an integrated miRNA similarity and disease similarity, providing more biological information for further feature learning. Second, MGADAE combines the known MDAs, disease similarity, and miRNA similarity into a heterogeneous network, then learns the representations of miRNAs and diseases through graph convolution operation. After that, an attention mechanism is introduced into MGADAE to integrate the representations from multiple graph convolutional network (GCN) layers. Lastly, the integrated representations of miRNAs and diseases are input into the bilinear decoder to obtain the final predicted association scores. Corresponding experiments prove that the proposed method outperforms existing advanced approaches in MDA prediction. Furthermore, case studies related to two human cancers provide further confirmation of the reliability of MGADAE in practice.

dsRID: in silico identification of dsRNA regions using long-read RNA-seq data

MOTIVATION

Double-stranded RNAs (dsRNAs) are potent triggers of innate immune responses upon recognition by cytosolic dsRNA sensor proteins. Identification of endogenous dsRNAs helps to better understand the dsRNAome and its relevance to innate immunity related to human diseases.

RESULTS

Here, we report dsRID (double-stranded RNA identifier), a machine-learning-based method to predict dsRNA regions in silico, leveraging the power of long-read RNA-sequencing (RNA-seq) and molecular traits of dsRNAs. Using models trained with PacBio long-read RNA-seq data derived from Alzheimer's disease (AD) brain, we show that our approach is highly accurate in predicting dsRNA regions in multiple datasets. Applied to an AD cohort sequenced by the ENCODE consortium, we characterize the global dsRNA profile with potentially distinct expression patterns between AD and controls. Together, we show that dsRID provides an effective approach to capture global dsRNA profiles using long-read RNA-seq data.

AVAILABILITY AND IMPLEMENTATION

Software implementation of dsRID, and genomic coordinates of regions predicted by dsRID in all samples are available at the GitHub repository: https://github.com/gxiaolab/dsRID.

Integration of miRNA dynamics and drought tolerant QTLs in rice reveals the role of miR2919 in drought stress response

To combat drought stress in rice, a major threat to global food security, three major quantitative trait loci for 'yield under drought stress' (qDTYs) were successfully exploited in the last decade. However, their molecular basis still remains unknown. To understand the role of secondary regulation by miRNA in drought stress response and their relation, if any, with the three qDTYs, the miRNA dynamics under drought stress was studied at booting stage in two drought tolerant (Sahbaghi Dhan and Vandana) and one drought sensitive (IR 20) cultivars. In total, 53 known and 40 novel differentially expressed (DE) miRNAs were identified. The primary drought responsive miRNAs were Osa-MIR2919, Osa-MIR3979, Osa-MIR159f, Osa-MIR156k, Osa-MIR528, Osa-MIR530, Osa-MIR2091, Osa-MIR531a, Osa-MIR531b as well as three novel ones. Sixty-one target genes that corresponded to 11 known and 4 novel DE miRNAs were found to be co-localized with the three qDTYs, out of the 1746 target genes identified. We could validate miRNA-mRNA expression under drought for nine known and three novel miRNAs in eight different rice genotypes showing varying degree of tolerance. From our study, Osa-MIR2919, Osa-MIR3979, Osa-MIR528, Osa-MIR2091-5p and Chr01_11911S14Astr and their target genes LOC_Os01g72000, LOC_Os01g66890, LOC_Os01g57990, LOC_Os01g56780, LOC_Os01g72834, LOC_Os01g61880 and LOC_Os01g72780 were identified as the most promising candidates for drought tolerance at booting stage. Of these, Osa-MIR2919 with 19 target genes in the qDTYs is being reported for the first time. It acts as a negative regulator of drought stress tolerance by modulating the cytokinin and brassinosteroid signalling pathway.

Turning off a few overexpressed genes in prostate cancer with microRNAs using a 7mer-seed match model

PURPOSE

This research aims to identify the miRNAs that could target the genes overexpressed in prostate cancer so that miRNA-based therapeutics could be developed.

METHODS

A 7mer-m8 model of microRNA targeting was utilized in order to analyse the relationship between microRNAs and overexpressed genes. The efficiency of miRNA binding was investigated using various parameters namely free energy (AMFE), GC and GC3 content, translation efficiency, cosine similarity metric, mRNA stability, free energy of RNA duplex, and base compositional difference. BLAST2GO software was used to elucidate the functional roles of the genes overexpressed in prostate cancer.

RESULTS

The current research reveals that the coding sequences of the genes were found targeted with multiple miRNAs. For instance, the HPN gene was targeted by the microRNA miR-4279 at two distinct sites i.e. 263-278 and 746-761 in the coding sequence. In the present study, it was observed that the target region of the genes exhibited a comparatively high GC and GC3 contents in comparison to the flanking regions. A low translational rate and weak relationship between RSCU and tRNA were obtained which may be due to the absence of optimal codons.

CONCLUSION

In this study, we have uncovered the human miRNAs that have potential for binding to the coding sequences of 14 most overexpressed genes in prostate cancer and thereby could silence those genes.

Potential role of the regulatory miR1119-MYC2 module in wheat (Triticum aestivum L.) drought tolerance

MicroRNA (miRNA)-target gene modules are essential components of plants' abiotic stress signalling pathways Little is known about the drought-responsive miRNA-target modules in wheat, but systems biology approaches have enabled the prediction of these regulatory modules and systematic study of their roles in responses to abiotic stresses. Using such an approach, we sought miRNA-target module(s) that may be differentially expressed under drought and non-stressed conditions by mining Expressed Sequence Tag (EST) libraries of wheat roots and identified a strong candidate (miR1119-MYC2). We then assessed molecular and physiochemical differences between two wheat genotypes with contrasting drought tolerance in a controlled drought experiment and assessed possible relationships between their tolerance and evaluated traits. We found that the miR1119-MYC2 module significantly responds to drought stress in wheat roots. It is differentially expressed between the contrasting wheat genotypes and under drought versus non-stressed conditions. We also found significant associations between the module's expression profiles and ABA hormone content, water relations, photosynthetic activities, H₂O₂ levels, plasma membrane damage, and antioxidant enzyme activities in wheat. Collectively, our results suggest that a regulatory module consisting of miR1119 and MYC2 may play an important role in wheat's drought tolerance.

Identification of Cumin (Cuminum cyminum) MicroRNAs through Deep Sequencing and Their Impact on Plant Secondary Metabolism

The pharmacological properties of plants lie in the content of secondary metabolites that are classified into different categories based on their biosynthesis, structures, and functions. MicroRNAs (miRNAs) are small non-coding RNA molecules that play crucial post-transcriptional regulatory roles in plants, including development and stress-response signaling; however, information about their involvement in secondary metabolism is still limited. Cumin is one of the most popular seeds from the plant Cuminum cyminum, with extensive applications in herbal medicine and cooking; nevertheless, no previous studies focus on the miRNA profile of cumin. In this study, the miRNA profile of C. cyminum and its association with the biosynthesis of secondary metabolites were determined using NGS technology. The sequencing data yielded 10,956,054 distinct reads with lengths ranging from 16 to 40 nt, of which 349 miRNAs were found to be conserved and 39 to be novel miRNAs. Moreover, this work identified 1959 potential target genes for C. cyminum miRNAs. It is interesting to note that several conserved and novel miRNAs have been found to specifically target important terpenoid backbone, flavonoid biosynthesis, and lipid/fatty acid pathways enzymes. We believe this investigation will aid in elucidating the implications of miRNAs in plant secondary metabolism.

Identification of Small RNAs Associated with Salt Stress in Chrysanthemums through High-Throughput Sequencing and Bioinformatics Analysis

The Chrysanthemum variety “Niu 9717” exhibits excellent characteristics as an ornamental plant and has good salt resistance. In this study, this plant was treated with 200 mM NaCl for 12 h followed by high-throughput sequencing of miRNA and degradome. Subsequently, the regulatory patterns of potential miRNAs and their target genes were searched to elucidate how Chrysanthemum miRNAs respond to salt. From the root and leaf samples, we identified a total of 201 known miRNAs belonging to 40 families; furthermore, we identified 79 new miRNAs, of which 18 were significantly differentially expressed (p < 0.05). The expressed miRNAs, which targeted a total of 144 mRNAs in the leaf and 215 mRNAs in the root, formed 144 and 226 miRNA–target pairs in roots and leaves, respectively. Combined with the miRNA expression profile, degradome and transcriptome data were then analyzed to understand the possible effects of the miRNA target genes and their pathways on salt stress. The identified genes were mostly located in pathways related to hormone signaling during plant growth and development. Overall, these findings suggest that conserved and novel miRNAs may improve salt tolerance through the regulation of hormone signal synthesis or expression of genes involved in hormone synthesis.

Dynamic Regulation of Grapevine's microRNAs in Response to Mycorrhizal Symbiosis and High Temperature

MicroRNAs (miRNAs) are non-coding small RNAs that play crucial roles in plant development and stress responses and can regulate plant interactions with beneficial soil microorganisms such as arbuscular mycorrhizal fungi (AMF). To determine if root inoculation with distinct AMF species affected miRNA expression in grapevines subjected to high temperatures, RNA-seq was conducted in leaves of grapevines inoculated with either Rhizoglomus irregulare or Funneliformis mosseae and exposed to a high-temperature treatment (HTT) of 40 °C for 4 h per day for one week. Our results showed that mycorrhizal inoculation resulted in a better plant physiological response to HTT. Amongst the 195 identified miRNAs, 83 were considered isomiRs, suggesting that isomiRs can be biologically functional in plants. The number of differentially expressed miRNAs between temperatures was higher in mycorrhizal (28) than in non-inoculated plants (17). Several miR396 family members, which target homeobox-leucine zipper proteins, were only upregulated by HTT in mycorrhizal plants. Predicted targets of HTT-induced miRNAs in mycorrhizal plants queried to STRING DB formed networks for Cox complex, and growth and stress-related transcription factors such as SQUAMOSA promoter-binding-like-proteins, homeobox-leucine zipper proteins and auxin receptors. A further cluster related to DNA polymerase was found in R. irregulare inoculated plants. The results presented herein provide new insights into miRNA regulation in mycorrhizal grapevines under heat stress and can be the basis for functional studies of plant-AMF-stress interactions.

A tissue-specific profile of miRNAs and their targets related to paeoniaflorin and monoterpenoids biosynthesis in Paeonia lactiflora Pall. by transcriptome, small RNAs and degradome sequencing

Paeonia lactiflora Pall. (Paeonia) has aroused many concerns due to its extensive medicinal value, in which monoterpene glucoside paeoniflorin and its derivatives are the active chemical components. However, little is known in the molecular mechanism of monoterpenoids biosynthesis, and the regulation network between small RNAs and mRNAs in monoterpenoids biosynthesis has not been investigated yet. Herein, we attempted to reveal the tissue-specific regulation network of miRNAs and their targets related to paeoniaflorin and monoterpenoids biosynthesis in Paeonia by combining mRNA and miRNA expression data with degradome analysis. In all, 289 miRNAs and 30177 unigenes were identified, of which nine miRNAs from seven miRNA families including miR396, miR393, miR835, miR1144, miR3638, miR5794 and miR9555 were verified as monoterpenoids biosynthesis-related miRNAs by degradome sequencing. Moreover, the co-expression network analysis showed that four monoterpenoids-regulating TFs, namely AP2, MYBC1, SPL12 and TCP2, were putatively regulated by five miRNAs including miR172, miR828, miR858, miR156 and miR319, respectively. The present study will improve our knowledge of the molecular mechanisms of the paeoniaflorin and monoterpenoids biosynthesis mediated by miRNA to a new level, and provide a valuable resource for further study on Paeonia.

Identification of Conserved and Novel MicroRNAs with their Targets in Garden Pea (Pisum Sativum L.) Leaves by High-Throughput Sequencing

MicroRNAs (miRNAs) are single-stranded, endogenous, non-coding RNAs of 20–24 nucleotides that play a significant role in post-transcriptional gene regulation. Various conserved and novel miRNAs have been characterized, especially from the plant species whose genomes were well-characterized; however, information on miRNA in economically important plants such as pea ( Pisum sativum L.) is limited. In this study, I have identified conserved and novel miRNA in garden pea plant leaves samples along with their targets by analyzing the next generation sequencing (NGS) data. The raw data obtained from NGS were processed and 1.38 million high-quality non-redundant reads were retained for analysis, this tremendous quantity of reads indicates a large and diverse small RNA population in pea leaves. After analyzing the deep sequencing data, 255 conserved and 11 novel miRNAs were identified in the garden pea leaves sample. Utilizing psRNATarget tool, the miRNA targets of conserved and novel miRNA were predicted. Further, the functional annotation of the miRNA targets were performed using blast2Go software and the target gene products were predicted. The miRNA target gene products along with GO_ID (Gene Ontology Identifier) were categorized into biological processes, cellular components, and molecular functions. The information obtained from this study will provide genomic resources that will help in understanding miRNA-mediated post-transcriptional gene regulation in garden peas.

Dynamic Analysis of microRNAs from Different Life Stages of Rhipicephalus microplus (Acari: Ixodidae) by High-Throughput Sequencing

MicroRNAs (miRNAs), which are small, noncoding RNA molecules, play an important regulatory role in gene expression at the posttranscriptional level. Relatively limited knowledge exists on miRNAs in Rhipicephalus microplus ticks in China; however, understanding the physiology of miRNA functions and expression at different developmental stages is important. In this study, three small RNA libraries were constructed for R. microplus eggs, larvae, and female adults; miRNAs were detected during these developmental stages by high-throughput sequencing, with 18,162,337, 8,090,736, and 11,807,326 clean reads, respectively. A total of 5132 known miRNAs and 31 novel miRNAs were identified. A total of 1736 differentially expressed miRNAs were significantly different at a p-value of <0.01; in female adults, 467 microRNAs were upregulated and 376 miRNAs downregulated compared to larval tick controls. Using larvae as controls, 218 upregulated and 203 downregulated miRNAs were detected in eggs; in eggs, 108 miRNAs were upregulated and 364 downregulated compared to female adults controls. To verify the reliability of the sequencing data, RT−qPCR was applied to compare expression levels of novel miRNAs. Some differentially expressed miRNAs are involved in developmental physiology, signal transduction, and cell-extracellular communications based on GO annotation and KEGG pathway analyses. Here, we provide a dynamic analysis of miRNAs in R. microplus and their potential targets, which has significance for understanding the biology of ticks and lays the foundation for improved understanding of miRNA functioning in the regulation of R. microplus development. These results can assist future miRNA studies in other tick species that have great significance for human and animal health.

Integrated analysis of small RNAs, transcriptome and degradome sequencing reveal the drought stress network in Agropyron mongolicum Keng

Agropyron mongolicum (A. mongolicum) is an excellent gramineous forage with extreme drought tolerance, which lives in arid and semiarid desert areas. However, the mechanism that underlies the response of microRNAs (miRNAs) and their targets in A. mongolicum to drought stress is not well understood. In this study, we analyzed the transcriptome, small RNAome (specifically the miRNAome) and degradome to generate a comprehensive resource that focused on identifying key regulatory miRNA-target circuits under drought stress. The most extended transcript in each collection is known as the UniGene, and a total of 41,792 UniGenes and 1,104 miRNAs were identified, and 99 differentially expressed miRNAs negatively regulated 1,474 differentially expressed target genes. Among them, eight miRNAs were unique to A. mongolicum, and there were 36 target genes. A weighted gene co-expression network analysis identified five hub genes. The miRNAs of five hub genes were screened with an integration analysis of the degradome and sRNAs, such as osa-miR444a-3p.2-MADS47, bdi-miR408-5p_1ss19TA-CCX1, tae-miR9774_L-2R-1_1ss11GT-carC, ata-miR169a-3p-PAO2, and bdi-miR528-p3_2ss15TG20CA-HOX24. The functional annotations revealed that they were involved in mediating the brassinosteroid signal pathway, transporting and exchanging sodium and potassium ions and regulating the oxidation-reduction process, hydrolase activity, plant response to water deprivation, abscisic acid (ABA) and the ABA-activated signaling pathway to regulate drought stress. Five hub genes were discovered, which could play central roles in the regulation of drought-responsive genes. These results show that the combined analysis of miRNA, the transcriptome and degradation group provides a useful platform to investigate the molecular mechanism of drought resistance in A. mongolicum and could provide new insights into the genetic engineering of Poaceae crops in the future.

CRISPR/Cas9-mediated fine-tuning of miRNA expression in tetraploid potato

MicroRNAs (miRNAs) are small noncoding RNAs, which modulate the abundance and spatiotemporal accumulation of target mRNAs at transcriptional and post-transcriptional levels and through that play important roles in several biological processes in plants. Here we show that in polyploid species, CRISPR/Cas9 system can be used for fine-tuning of miRNA expression, which can have broader range of applications compared to knock-out mutants. We established the complete pipeline for CRISPR-Cas9-mediated modulation of miRNA expression in potato. It consists of (1) design and assembly of dual sgRNA CRISPR/Cas9 constructs, (2) transient transfection of protoplasts following fast and efficient screening by high resolution melting analysis to select functional sgRNAs, and (3) stable transformation of potato explants with functional sgRNAs and selection of regenerated transgenic lines with desired mutations and desired miRNA abundance based on sequencing and RT-qPCR. We show that miRNA-editing using dual sgRNA approach results in different types of mutations among transgenic lines but also in different alleles of the same plant, which are target site-dependent. The most frequent were short deletions, but we also detected 1-nt insertions (T or G), deletions between two sgRNAs and larger deletions. miRNA abundance correlates with the frequency and type of introduced mutations, as more extensive mutations in more alleles result in lower miRNA abundance. Interestingly, some mutated loci can generate alternative miRNAs, now novel targets were however predicted for those. In all transgenic lines with Cas9 expression, we detected mutations, suggesting high efficiency of Cas9-editing. We confirmed the miRNA-editing efficiency of our optimised approach in two different potato genotypes and three different loci.

Identification of Long Noncoding RNAs Involved in Eyelid Pigmentation of Hereford Cattle

Several ocular pathologies in cattle, such as ocular squamous cell carcinoma and infectious keratoconjunctivitis, have been associated with low pigmentation of the eyelids. The main objective of this study was to analyze the transcriptome of eyelid skin in Hereford cattle using strand-specific RNA sequencing technology to characterize and identify long noncoding RNAs (lncRNAs). We compared the expression of lncRNAs between pigmented and unpigmented eyelids and analyzed the interaction of lncRNAs and putative target genes to reveal the genetic basis underlying eyelid pigmentation in cattle. We predicted 4,937 putative lncRNAs mapped to the bovine reference genome, enriching the catalog of lncRNAs in Bos taurus. We found 27 differentially expressed lncRNAs between pigmented and unpigmented eyelids, suggesting their involvement in eyelid pigmentation. In addition, we revealed potential links between some significant differentially expressed lncRNAs and target mRNAs involved in the immune response and pigmentation. Overall, this study expands the catalog of lncRNAs in cattle and contributes to a better understanding of the biology of eyelid pigmentation.

Identification of Known and Novel Arundo donax L. MicroRNAs and Their Targets Using High-Throughput Sequencing and Degradome Analysis

MicroRNAs (miRNAs) are a class of non-coding molecules involved in the regulation of a variety of biological processes. They have been identified and characterized in several plant species, but only limited data are available for Arundo donax L., one of the most promising bioenergy crops. Here we identified, for the first time, A. donax conserved and novel miRNAs together with their targets, through a combined analysis of high-throughput sequencing of small RNAs, transcriptome and degradome data. A total of 134 conserved miRNAs, belonging to 45 families, and 27 novel miRNA candidates were identified, along with the corresponding primary and precursor miRNA sequences. A total of 96 targets, 69 for known miRNAs and 27 for novel miRNA candidates, were also identified by degradome analysis and selected slice sites were validated by 5′-RACE. The identified set of conserved and novel candidate miRNAs, together with their targets, extends our knowledge about miRNAs in monocots and pave the way to further investigations on miRNAs-mediated regulatory processes in A. donax, Poaceae and other bioenergy crops.

miR390 family of Cymbidium goeringii is involved in the development of reproductive organs in transgenic Arabidopsis

BACKGROUND

miR390s is an ancient family with a high level of conservation among plant miRNAs. Through the auxin signal transduction pathway, miR390 participates in diverse biological processes of plant growth and development. As an important Chinese traditional orchid, Cymbidium goeringii has unique flower shape and elegant fragrance. But its development has been greatly restricted because of the low flower bud differentiation and the difficult reproduction. This study aims to provide guidance for the role of cgo-miR390 in reproductive organ development to enhance the ornamental and economic value of Cymbidium.

RESULTS

MIR390a, MIR390b and MIR390c of C. goeringii were cloned, and their length ranged from 130 to 150 nt. Each precursor sequence of cgo-miR390 contains 2 to 3 mature miRNAs. Three kinds of cgo-miR390s displayed distinct temporal and spatial expression patterns during floral development in C. goeringii. The overexpression of MIR390s alters morphology and function of stamens and pistils in Arabidopsis, such as enlargement of anther aspect ratio and separation of stylar and stigmas, which affects the development of fruits and seeds. In particular, the pollen amount decreased and the seed abortion rate increased in cgo-MIR390c-overexpressed plants.

CONCLUSIONS

cgo-miR390 family affected the development of reproductive organs in transgenic Arabidopsis. The study provides references for the genetic improvement for orchid with potentially great economic benefit.

Understanding the Molecular Mechanisms of Betel miRNAs on Human Health

BACKGROUND

Since ancient times, "betel leaf" (Piper betle) has been revered for its religious, cultural, and medicinal properties. Phytochemicals from the Piper betle are effective in a variety of conditions, including cancer. To date, however, no genomic study or evidence has been found to elucidate the regulatory mechanism that underpins its therapeutic properties. This is the first study of its kind to predict Piper betle miRNAs and also the first genomics source representation of Piper betle. According to previous research, miRNAs from the plants we eat can regulate gene expression. In line with this, our in-silico study revealed that Piper betle and human cross-kingdom control occurs.

METHOD

This study demonstrates the prediction and in-silico validation of Piper betle miRNAs from NGS-derived transcript sequences. The cross-kingdom regulation which can also be understood as inter-species RNA regulation was studied to identify human mRNA targets being controlled by Piper betle miRNAs. Functional annotation and gene-disease association of human targets were performed to understand the role of Piper betle miRNAs in human health and disease. The protein-protein interaction and expression study of targets was further carried out to decipher their role in cancer development.

RESULTS

Identified six Piper betle miRNAs belonging to miR156, miR164, miR172, and miR535 families were discovered to target 198 human mRNAs involved in various metabolic and disease processes. Angiogenesis and the cell surface signaling pathway were the most enriched gene ontology correlated with targets, both of which play a critical role in disease mechanisms, especially in the case of carcinoma. In an analysis of gene-disease interactions, 40 genes were found to be related to cancer. According to a protein-protein interaction, the CDK6 gene, which is thought to be a central regulator of cell cycle progression, was found as a hub protein, affecting the roles of CBFB, SAMD9, MDM4, AXIN2, and NOTCH2 onco genes. Further investigation revealed that pbe-miRNA164a can be used as a regulator to minimise disease severity in Acute Myeloid Leukemia, where CDK6 expression is highest compared to normal cells.

CONCLUSION

The predicted pbe-miRNA164a in this study can be a promising suppressor of CDK6 gene involved in tumour angiogenesis. In vivo validation of the pbe-miRNA164a mimic could pave the way for new opportunities to fight cancer and leverage the potential of Piper betle in the healthcare sector.

Binding sites of miRNA on the overexpressed genes of oral cancer using 7mer-seed match

The microRNAs having a length of ~ 19-22 nucleotides are the small, non-coding RNAs. The evolution of microRNAs in many disorders may hold the key to tackle complex challenges. Oral cancer belongs to the group of head and neck cancer. It occurs in the mouth region that appears as an ulcer. In this study, we collected information on the overexpressed genes of oral cancer. The coding sequences of the genes were derived from NCBI and the entire set of human microRNAs present in miRBASE 21 was retrieved. The human microRNAs that can target the overexpressed genes of oral cancer were determined with the aid of our in-house software. The interaction between microRNAs and the overexpressed genes was evaluated with 7mer-m8 model of microRNA targeting. The genes DKK1 and APLN paired with only one miRNA i.e., miR-447 and miR-6087, respectively. But the genes INHBA and MMP1 were found to be targeted by 2 miRNAs, while the genes FN1, FAP, TGFPI, COL4A1, COL4A2, and LOXL2 were found to be targeted by 16, 5, 9, 18, 29, and 11 miRNAs. Subsequently, several measures such as free energy, translation efficiency, and cosine similarity metric were used to estimate the binding process. It was found that the target region's stability was higher than the upstream and downstream zones. The overexpressed genes' GC contents were calculated, revealing that the codons in target miRNA region were overall GC rich as well as GC3 rich. Lastly, gene ontology was performed to better understand each gene's involvement in biological processes, molecular function, and cellular component. Our study showed the role of microRNAs in gene repression, which could possibly aid in the prognosis and diagnosis of oral cancer.

MicroRNA and Degradome Profiling Uncover Defense Response of Fraxinus velutina Torr. to Salt Stress

Soil salinization is a major environmental problem that seriously threatens the sustainable development of regional ecosystems and local economies. Fraxinus velutina Torr. is an excellent salt-tolerant tree species, which is widely planted in the saline-alkaline soils in China. A growing body of evidence shows that microRNAs (miRNAs) play important roles in the defense response of plants to salt stress; however, how miRNAs in F. velutina exert anti-salt stress remains unclear. We previously identified two contrasting F. velutina cuttings clones, salt-tolerant (R7) and salt-sensitive (S4) and found that R7 exhibits higher salt tolerance than S4. To identify salt-responsive miRNAs and their target genes, the leaves and roots of R7 and S4 exposed to salt stress were subjected to miRNA and degradome sequencing analysis. The results showed that compared with S4, R7 showed 89 and 138 differentially expressed miRNAs in leaves and roots, respectively. Specifically, in R7 leaves, miR164d, miR171b/c, miR396a, and miR160g targeting NAC1, SCL22, GRF1, and ARF18, respectively, were involved in salt tolerance. In R7 roots, miR396a, miR156a/b, miR8175, miR319a/d, and miR393a targeting TGA2.3, SBP14, GR-RBP, TCP2/4, and TIR1, respectively, participated in salt stress responses. Taken together, the findings presented here revealed the key regulatory network of miRNAs in R7 responding to salt stress, thereby providing new insights into improving salt tolerance of F. velutina through miRNA manipulation.

Identification of miRNA and their target genes in Cestrum nocturnum L. and Cestrum diurnum L. in stress responses

MicroRNAs (miRNAs) are small, highly conserved non-coding RNA molecules and products of primary miRNAs that regulate the target gene expression. Homology-based approaches were employed to identify miRNAs and their targets in Cestrum nocturnum L. and Cestrum diurnum L. A total of 32 and 12 miRNA candidates were identified in C. nocturnum and C. diurnum. These miRNAs belong to 26 and 10 miRNA families and regulate 1024 and 1007 target genes in C. nocturnum, and C. diurnum, respectively. The functional roles of these miRNAs have not been earlier elucidated in Cestrum. MiR815a, miR849, miR1089 and miR172 have a strong propensity to target genes controlling phytochrome-interacting factor 1 (PIF1), ubiquitin-specific protease 12 (UBP12), leucine-rich repeat (LRR) protein kinase and GAI, RGA, SCR (GRAS) family transcription factor in C. nocturnum. While miR5205a, miR1436 and miR530 regulate PATATIN-like protein 6 (PLP6), PHD finger transcription factor and myb domain protein 48 (MYB48) in C. diurnum. Overall, these miRNAs have regulatory responses in biotic and abiotic stresses in both plant species. Eight putative miRNAs and their target genes were selected for qRT-PCR validation. The validated results suggested the importance of miR815a, miR849, miR5205a, miR1089, miR172, miR1436, and miR530 in exerting control over stress responses in C. nocturnum and C. diurnum.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-022-01127-1.

Identification of microRNAs from Medicinal Plant Murraya koenigii by High-Throughput Sequencing and Their Functional Implications in Secondary Metabolite Biosynthesis

MicroRNAs (miRNAs) are small noncoding RNA molecules that play crucial post-transcriptional regulatory roles in plants, including development and stress-response signaling. However, information about their involvement in secondary metabolism is still limited. Murraya koenigii is a popular medicinal plant, better known as curry leaves, that possesses pharmaceutically active secondary metabolites. The present study utilized high-throughput sequencing technology to investigate the miRNA profile of M. koenigii and their association with secondary metabolite biosynthesis. A total of 343,505 unique reads with lengths ranging from 16 to 40 nt were obtained from the sequencing data, among which 142 miRNAs were identified as conserved and 7 as novel miRNAs. Moreover, 6078 corresponding potential target genes of M. koenigii miRNAs were recognized in this study. Interestingly, several conserved and novel miRNAs of M. koenigii were found to target key enzymes of the terpenoid backbone and the flavonoid biosynthesis pathways. Furthermore, to validate the sequencing results, the relative expression of eight randomly selected miRNAs was determined by qPCR. To the best of our knowledge, this is the first report of the M. koenigii miRNA profile that may provide useful information for further elucidation of the involvement of miRNAs in secondary metabolism. These findings might be crucial in the future to generate artificial-miRNA-based, genetically engineered M. koenigii plants for the overproduction of medicinally highly valuable secondary metabolites.

In Silico Investigation of Conserved miRNAs and Their Targets From the Expressed Sequence Tags in Neospora Caninum Genome

Neospora caninum is a protozoan parasite, the etiologic agent of Neosporosis—a common cause of abortion in cattle worldwide. Herd level prevalence of Neosporosis could be as high as 90%. However, there is no approved treatment and vaccines available for Neosporosis. MicroRNA (miRNA) based prophylaxis and therapeutics could be options for Neosporosis in cattle and other animals. The current study aimed to investigate the genome of Neospora caninum to identify and characterize the conserved miRNAs through Expressed Sequence Tags (ESTs) dependent homology search. A total of 1,041 mature miRNAs of reference organisms were employed against 336 non-redundant ESTs available in the genome of Neospora caninum. The study predicted one putative miRNA “nca-miR-9388-5p” of 19 nucleotides with MFEI value -1.51 kcal/mol and (A + U) content% 72.94% corresponding with its pre-miRNA. A comprehensive search for specific gene targets was performed and discovered 16 potential genes associated with different protozoal physiological functions. Significantly, the gene “Protein phosphatase” was found responsible for the virulence of Neospora caninum. The other genes were accounted for gene expression, vesicular transport, cell signaling, cell proliferation, DNA repair mechanism, and different developmental stages of the protozoon. Therefore, this study finding will provide pivotal information to future aspirants upon Bovine Neosporosis. It will also serve as the baseline information for further studies of the bioinformatics approach to identify other protozoal miRNAs.

Identification of novel microRNAs in Rous sarcoma Virus (RSV) and their target sites in tumor suppressor genes of chicken

A small non-coding, evolutionarily conserved regulatory RNA molecule known as microRNA (miRNA) regulates various cellular activities and pathways. MicroRNAs remain evolutionarily conserved in different species of same taxa. They are present in all organisms including viruses. Viral miRNAs are small, less conserved and less stable and have higher negative minimal folding free energy than miRNAs of different organisms. The size of viral precursor miRNA is approximately 60-119 nucleotides in length. The structure of the mature miRNA sequences is predicted by using higher negative MFE (ΔG) value. Rous sarcoma Virus (RSV), named after its inventor Peyton Rous, has been known for causing tumors in the chicken for which it is known as an oncogenic retrovirus. Using specific criteria we have predicted 5 potential miRNAs in RSV which targeted 8 tumor suppressor genes in Gallus gallus. This study aims to predict the potential miRNAs, secondary structures and their targets for better understanding of the regulatory network of Rous sarcoma virus miRNA in forming sarcoma.

Identification and profiling of conserved microRNAs in different developmental stages of crown imperial (Fritillaria imperialis L.) using high-throughput sequencing

BACKGROUND

Novel strategies for improvement of ornamental plants and their properties relay on miRNA control of differential plant gene expression modulation. Still, in response to the same abiotic stresses, some conserved miRNA families show different expression patterns in different plant species. In parallel, the use of deep sequencing technologies reveals new levels of complexity of regulatory networks in plants through identification of new miRNAs.

METHODS AND RESULTS

Fritillaria imperialis plants were collected from their natural habitats in Koohrang, Chaharmahal va Bakhtiari, Iran. Several tissues including stamen, pistil, petal, sepal, leaf, stem, bulb and fruit were collected during three developmental stages (stem elongation, flower development and seed head stages). Using RNAseq and qRT-PCR approach, this research revealed 21 conserved miRNAs, matching 15 miRNA families, in Fritilaria imperialis.

CONCLUSIONS

The expression of seven conserved miRNAs (Fim-miR156b, Fim-miR159, Fim-miR166a-5p, Fim-miR169d-5p, Fim-miR171c, Fim-miR393 and Fim-miR396e-3p) was further investigated in different tissues and three developmental stages, suggesting different roles for these miRNAs during growth and development of crown imperial. Gained knowledge from this research can open the door to find efficient ways to secure crown imperial survival, preservation and utilization and if proven useful may be applied in other plant species as well.

Genome-Wide Analysis of Alternative Splicing and Non-Coding RNAs Reveal Complicated Transcriptional Regulation in Cannabis sativa L

It is of significance to mine the structural genes related to the biosynthetic pathway of fatty acid (FA) and cellulose as well as explore the regulatory mechanism of alternative splicing (AS), microRNAs (miRNAs) and long non-coding RNAs (lncRNAs) in the biosynthesis of cannabinoids, FA and cellulose, which would enhance the knowledge of gene expression and regulation at post-transcriptional level in Cannabis sativa L. In this study, transcriptome, small RNA and degradome libraries of hemp 'Yunma No.1' were established, and comprehensive analysis was performed. As a result, a total of 154, 32 and 331 transcripts encoding key enzymes involved in the biosynthesis of cannabinoids, FA and cellulose were predicted, respectively, among which AS occurred in 368 transcripts. Moreover, 183 conserved miRNAs, 380 C. sativa-specific miRNAs and 7783 lncRNAs were predicted. Among them, 70 miRNAs and 17 lncRNAs potentially targeted 13 and 17 transcripts, respectively, encoding key enzymes or transporters involved in the biosynthesis of cannabinoids, cellulose or FA. Finally, the crosstalk between AS and miRNAs or lncRNAs involved in cannabinoids and cellulose was also predicted. In summary, all these results provided insights into the complicated network of gene expression and regulation in C. sativa.

Identification of putative microRNAs in the complete genome of Mycobacterium avium and their possible interaction with human transcripts

The grievous adversity regarding Mycobacterium avium is its ubiquitous nature. Isolation of the bacteria from drinking water, house dust, and soil, etc., is an alarming issue for the scientific community. The microRNAs are the molecular influencers of gene expression that act during the process of post transcription. A few reports claimed the existence of microRNAs or microRNA-like molecules in the prokaryotic species. Biogenesis of bacterial miRNAs requires their transport into the host cell. Subsequently, the host-encoded enzymes are exerted for the formation of bacterial mature miRNAs and their regulation. In our study, the screening of complete genome of Mycobacterium avium revealed six putative precursor microRNA sequences bearing typical secondary structures. The mature microRNAs were predicted in both arms of the secondary structures. A total of 12 possible mature microRNAs were identified in this study. The likely targets of the predicted mature miRNAs were searched in human 3' UTR. In the human transcriptome, 193 genes were possibly targeted by 12 mature miRNAs of Mycobacterium avium. The essential functionalities of the target genes included signal transduction, immune system, DNA binding, and response to stress.

Genome-wide prediction of cauliflower miRNAs and lncRNAs and their roles in post-transcriptional gene regulation

We have predicted miRNAs, their targets and lncRNAs from the genome of Brassica oleracea along with their functional annotation. Selected miRNAs and their targets are experimentally validated. Roles of these non-coding RNAs in post-transcriptional gene regulation are also deciphered. Cauliflower (Brassica oleracea var. Botrytis) is an important vegetable crop for its dietary and medicinal values with rich source of vitamins, dietary fibers, flavonoids and antioxidants. MicroRNAs (miRNAs) are small non-coding RNAs (ncRNAs), which regulate gene expression by inhibiting translation or by degrading messenger RNAs (mRNAs). On the other hand, long non-coding RNAs (lncRNAs) are responsible for the up regulation and the down regulation of transcription. Although the genome of cauliflower is reported, yet the roles of these ncRNAs in post-transcriptional gene regulation (PTGR) remain elusive. In this study, we have computationally predicted 355 miRNAs, of which 280 miRNAs are novel compared to miRBase 22.1. All the predicted miRNAs belong to 121 different families. We have also identified 934 targets of 125 miRNAs along with their functional annotation. These targets are further classified into biological processes, molecular functions and cellular components. Moreover, we have predicted 634 lncRNAs, of which 61 are targeted by 30 novel miRNAs. Randomly chosen 10 miRNAs and 10 lncRNAs are experimentally validated. Five miRNA targets including squamosa promoter-binding-like protein 9, homeobox-leucine zipper protein HDG12-like, NAC domain-containing protein 100, CUP-SHAPED COTYLEDON 1 and kinesin-like protein NACK2 of four miRNAs including bol-miR156a, bol-miR162a, bol-miR164d and bol-miR2673 are also experimentally validated. We have built network models of interactions between miRNAs and their target mRNAs, as well as between miRNAs and lncRNAs. Our findings enhance the knowledge of non-coding genome of cauliflower and their roles in PTGR, and might play important roles in improving agronomic traits of this economically important crop.

Identification and Characterization of Verticillium nonalfalfae-Responsive MicroRNAs in the Roots of Resistant and Susceptible Hop Cultivars

MicroRNAs are 21- to 24-nucleotide-long, non-coding RNA molecules that regulate gene expression at the post-transcriptional level. They can modulate various biological processes, including plant response and resistance to fungal pathogens. Hops are grown for use in the brewing industry and, recently, also for the pharmaceutical industry. Severe Verticillium wilt caused by the phytopathogenic fungus Verticillium nonalfalfae, is the main factor in yield loss in many crops, including hops (Humulus lupulus L.). In our study, we identified 56 known and 43 novel miRNAs and their expression patterns in the roots of susceptible and resistant hop cultivars after inoculation with V. nonalfalfae. In response to inoculation with V. nonalfalfae, we found five known and two novel miRNAs that are differentially expressed in the susceptible cultivar and six known miRNAs in the resistant cultivar. Differentially expressed miRNAs target 49 transcripts involved in protein localization and pigment synthesis in the susceptible cultivar, whereas they are involved in transcription factor regulation and hormone signalling in the resistant cultivar. The results of our study suggest that the susceptible and resistant hop cultivars respond differently to V. nonalfalfae inoculation at the miRNA level and that miRNAs may contribute to the successful defence of the resistant cultivar.

Evolution and Phylogeny of MicroRNAs - Protocols, Pitfalls, and Problems

MicroRNAs are important regulators in many eukaryotic lineages. Typical miRNAs have a length of about 22nt and are processed from precursors that form a characteristic hairpin structure. Once they appear in a genome, miRNAs are among the best-conserved elements in both animal and plant genomes. Functionally, they play an important role in particular in development. In contrast to protein-coding genes, miRNAs frequently emerge de novo. The genomes of animals and plants harbor hundreds of mutually unrelated families of homologous miRNAs that tend to be persistent throughout evolution. The evolution of their genomic miRNA complement closely correlates with important morphological innovation. In addition, miRNAs have been used as valuable characters in phylogenetic studies. An accurate and comprehensive annotation of miRNAs is required as a basis to understand their impact on phenotypic evolution. Since experimental data on miRNA expression are limited to relatively few species and are subject to unavoidable ascertainment biases, it is inevitable to complement miRNA sequencing by homology based annotation methods. This chapter reviews the state of the art workflows for homology based miRNA annotation, with an emphasis on their limitations and open problems.

Identification and function prediction of novel microRNAs in adenosine monophosphate activated protein kinase-activated Sertoli cells of immature boar

This study was carried out with the objective to identify function prediction of novel microRNAs (miRNAs) in immature boar Sertoli cells (SCs) treated with 5-aminoimidazole-4-carboxamide-1-β-D-ribofuranoside (AICAR), which is an agonist of adenosine monophosphate-activated protein kinase (AMPK) for regulating cellular energy homeostasis. Two small RNA libraries (control and AICAR treatment) prepared from immature boar SCs were constructed and sequenced by the Illumina small RNA deep sequencing. We identified 77 novel miRNAs and predicted 177 potential target genes for 26 differential novel miRNAs (four miRNAs up-regulation and 22 miRNAs down-regulation) in AICAR-treated SCs. Gene Ontology enrichment and Kyoto Encyclopedia of Genes and Genomes pathway suggested that target genes of differential novel miRNAs were implicated in many biological processes and metabolic pathways. Our findings provided useful information for the functional regulation of novel miRNAs and target mRNAs on AMPK-activated immature boar SCs.

Identifying potential human and medicinal plant microRNAs against SARS-CoV-2 3'UTR region: A computational genomics assessment

The coronavirus disease of 2019 (COVID-19) began as an outbreak and has taken a toll on human lives. The current pandemic requires scientific attention; hence we designed a systematic computational workflow to identify the cellular microRNAs (miRNAs) from human host possessing the capability to target and silence 3′UTR of SARS-CoV-2 genome. Based on this viewpoint, we extended our miRNA search to medicinal plants like Ocimum tenuiflorum, Zingiber officinale and Piper nigrum, which are well-known to possess antiviral properties, and are often consumed raw or as herbal decoctions. Such an approach, that makes use of miRNA of one species to interact and silence genes of another species including viruses is broadly categorized as cross-kingdom interactions. As a part of our genomics study on host-virus-plant interaction, we identified one unique 3′UTR conserved site ‘GGAAGAG’ amongst 5024 globally submitted SARS-CoV-2 complete genomes, which can be targeted by the human miRNA ‘hsa-miR-1236–3p’ and by Z. officinale miRNA ‘zof-miR2673b’. Additionally, we also predicted that the members of miR477 family commonly found in these three plant genomes possess an inherent potential to silence viral genome RNA and facilitate antiviral defense against SARS-CoV-2 infection. In conclusion, this study reveals a universal site in the SARS-CoV-2 genome that may be crucial for targeted therapeutics to cure COVID-19.

Identification and host response interaction study of SARS-CoV-2 encoded miRNA-like sequences: an in silico approach

COVID-19, a global pandemic caused by an RNA virus named SARS-CoV-2 has brought the world to a standstill in terms of infectivity, casualty, and commercial plummet. RNA viruses can encode microRNAs (miRNAs) capable of modulating host gene expression, and with that notion, we aimed to predict viral miRNA like sequences of MERS-CoV, SARS-CoV and SARS-CoV-2, analyze sequence reciprocity and investigate SARS-CoV-2 encoded potential miRNA-human genes interaction using bioinformatics tools. In this study, we retrieved 206 SARS-CoV-2 genomes, executed phylogenetic analysis, and the selected reference genome (MT434792.1) exhibited about 99% similarities among the retrieved genomes. We predicted 402, 137, and 85 putative miRNAs of MERS-CoV (NC_019843.3), SARS-CoV (NC_004718.3), and SARS-CoV-2 (MT434792.1) genome, respectively. Sequence similarity was analyzed among 624 miRNAs which revealed that the predicted miRNAs of SARS-CoV-2 share a cluster with the clad of miRNAs from MERS-CoV and SARS-CoV. Only SARS-CoV-2 derived 85 miRNAs were encountered for target prediction and 29 viral miRNAs seemed to target 119 human genes. Moreover, Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) analysis suggested the involvement of respective genes in various pathways and biological processes. Finally, we focused on eight putative miRNAs influencing 14 genes that are involved in the adaptive hypoxic response, neuroinvasion and hormonal regulation, and tumorigenic progression in patients with COVID-19. SARS-CoV-2 encoded miRNAs may cause misexpression of some critical regulators and facilitate viral neuroinvasion, altered hormonal axis, and tumorigenic events in the human host. However, these propositions need validation from future studies.

In Silico Identification and Functional Characterization of Conserved miRNAs in the Genome of Cryptosporidium parvum

Cryptosporidium parvum, a predominant causal agent of a fatal zoonotic protozoan diarrhoeal disease called cryptosporidiosis, bears a worldwide public health concern for childhood mortality and poses a key threat to the dairy and water industries. MicroRNAs (miRNAs), small but powerful posttranscriptional gene silencing RNA molecules, regulate a variety of molecular, biological, and cellular processes in animals and plants. As to the present date, there is a paucity of information regarding miRNAs of C. parvum; hence, this study was used to identify miRNAs in the organism using a comprehensible expressed sequence tag-based homology search approach consisting of a series of computational screening process from the identification of putative miRNA candidates to the functional annotation of the important gene targets in C. parvum. The results revealed a conserved miRNA that targeted 487 genes in the model organism (Drosophila melanogaster) and 85 genes in C. parvum, of which 11 genes had direct involvements in several crucial virulence factors such as environmental oocyst protection, excystation, locomotion, adhesion, invasion, stress protection, intracellular growth, and survival. Besides, 20 genes showed their association with various major pathways dedicated for the ribosomal biosynthesis, DNA repair, transportation, protein production, gene expression, cell cycle, cell proliferation, development, immune response, differentiation, and nutrient metabolism of the organism in the host. Thus, this study provides a strong evidence of great impact of identified miRNA on the biology, virulence, and pathogenesis of C. parvum. Furthermore, the study suggests that the detected miRNA could be a potential epigenomic tool for controlling the protozoon through silencing those virulent and pathway-related target genes.

Computational Identification of miRNAs and Temperature-Responsive lncRNAs From Mango (Mangifera indica L.)

Mango is a major tropical fruit in the world and is known as the king of fruits because of its flavor, aroma, taste, and nutritional values. Although various regulatory roles of microRNAs (miRNAs) and long non-coding RNAs (lncRNAs) have been investigated in many plants, there is yet an absence of such study in mango. This is the first study to provide information on non-coding RNAs (ncRNAs) of mango with the aims of identifying miRNAs and lncRNAs and discovering their potential functions by interaction prediction of the miRNAs, lncRNAs, and their target genes. In this analysis, about a hundred miRNAs and over 7,000 temperature-responsive lncRNAs were identified and the target genes of these ncRNAs were characterized. According to these results, the newly identified mango ncRNAs, like other plant ncRNAs, have a potential role in biological and metabolic pathways including plant growth and developmental process, pathogen defense mechanism, and stress-responsive process. Moreover, mango lncRNAs can target miRNAs to reduce the stability of lncRNAs and can function as molecular decoys or sponges of miRNAs. This paper would provide information about miRNAs and lncRNAs of mango and would help for further investigation of the specific functions of mango ncRNAs through wet lab experiments.

Identification and characterization of heat-responsive microRNAs at the booting stage in two rice varieties, 9311 and Nagina 22

MicroRNAs (miRNAs) are small, non-coding, regulatory RNAs that play important roles in abiotic stress responses in plants, but their regulatory roles in the adaptive response to heat stress at the booting stage in two rice varieties, 9311 and Nagina 22, remain largely unknown. In this study, 464 known miRNAs and 123 potential novel miRNAs were identified. Of these miRNAs, a total of 90 differentially expressed miRNAs were obtained with 9311 libraries as the control group, of which 54 were upregulated and 36 were downregulated. To gain insight into functional significance, 2773 potential target genes of these 90 differentially expressed miRNAs were predicted. GO enrichment analysis showed that the predicted target genes of differentially expressed miRNAs included NACs, LACs, CSD, and Hsp40. KEGG pathway analysis showed that the target genes of these differentially expressed miRNAs were significantly enriched in the plant hormone signal transduction pathway. The expression levels of 10 differentially expressed miRNAs and their target genes obtained by qRT-PCR were largely consistent with the sequencing results. This study lays a foundation for the elucidation of the miRNA-mediated regulatory mechanisms in rice at elevated temperatures.

Silencing lung cancer genes using miRNAs identified by 7mer-seed matching

Lung cancer (LC) is the main cause of cancer-associated deaths in both men and women globally with a very high mortality rate. The microRNAs (miRNAs) are a class of noncoding RNAs consisting of 18-25 nucleotides. They inhibit translation of protein through binding to complementary target mRNAs. The non-coding miRNAs are recognized as potent biomarkers for detection, development and treatment of malignancy. In this study, we screened a set of 12 genes over expressed in small cell lung cancer, non small cell lung cancer and the genes involved in both categories and their binding sites for human miRNAs as no work was reported yet. Screening of human miRNAs revealed that a few genes showed numerous miRNA binding sites. Free energy values of mRNA sequences revealed that they might acquire compact folded structure causing complexity for miRNAs to interact. GC content in the target site was relatively higher than that of their flanks. It was observed through analysis of cosine similarity metric and compAI parameters that the genes related to lung cancer were encoded with non optimal codons and thus might be translationally less efficient for producing polypeptides. Gene ontology analysis was carried out to understand the diverse functions of these 12 genes.

In silico identification of conserved miRNAs in the genome of fibre biogenesis crop Corchorus capsularis

Corchorus capsularis, commonly known as jute occupies the leading position in the production of natural fibre alongside lower environmental threat. Small noncoding ~21 to 24 nucleotides long microRNAs play significant roles in regulating the gene expression as well as different functions in cellular growth and development. Here, the study adopted a comprehensive in silico approach to identify and characterize the conserved miRNAs in the genome of C. capsularis including functional annotation of specific gene targets. Expressed Sequence Tags (ESTs) based homology search of 3350 known miRNAs of dicotyledons were allowed against 763 non-redundant ESTs of jute genome, resulted in the prediction of 5 potential miRNA candidates belonging five different miRNA families (miR1536, miR9567-3p, miR4391, miR11300, and miR8689). The putative miRNAs were composed of 18 nucleotides having a range of -0.49 to -1.56 MFEI values and 55%–61% of (A + U) content in their pre-miRNAs. A total of 1052 gene targets of putative miRNAs were identified and their functions were extensively analyzed. Most of the gene targets were involved in plant growth, cell cycle regulation, organelle synthesis, developmental process and environmental responses. Five gene targets, namely, NAC Domain Containing Protein, WRKY DNA binding protein, 3-dehydroquinate synthase, S-adenosyl-L-Met–dependent methyl transferase and Vascular-related NAC-Domain were found to be involved in the lignin biosynthesis, phenylpropanoid pathways and secondary wall formation. The present study might accelerate the more miRNA discovery, strengthening the complete understanding of miRNAs association in the cellular basis of lignin biosynthesis towards the production of high standard jute products.

Small RNA and degradome deep sequencing reveals important roles of microRNAs in cotton (Gossypium hirsutum L.) response to root-knot nematode Meloidogyne incognita infection

Investigation of cotton response to nematode infection will allow us to better understand the cotton immune defense mechanism and design a better biotechnological approach for efficiently managing pest nematodes in cotton. In this study, we firstly treated cotton by root knot nematode (RKN, Meloidogyne incognita) infections, then we employed the high throughput deep sequencing technology to sequence and genome-widely identify all miRNAs in cotton; finally, we analyzed the functions of these miRNAs in cotton response to RKN infections. A total of 266 miRNAs, including 193 known and 73 novel miRNAs, were identified by deep sequencing technology, which belong to 67 conserved and 66 novel miRNA families, respectively. A majority of identified miRNA families only contain one miRNA; however, miR482 family contains 14 members and some others contain 2-13 members. Certain miRNAs were specifically expressed in RKN-infected cotton roots and others were completely inhibited by RKN infection. A total of 50 miRNAs were differentially expressed after RKN infection, in which 28 miRNAs were up-regulated and 22 were inhibited by RKN treatment. Based on degradome sequencing, 87 gene targets were identified to be targeted by 57 miRNAs. These miRNA-targeted genes are involved in the interaction of cotton plants and nematode infection. Based on GO (gene ontology) and KEGG (Kyoto Encyclopedia of Genes and Genomes) analysis, 466 genes from all 636 miRNA targets were mapped to 6340 GO terms, 181 genes from 228 targets of differentially expressed miRNAs were mapped to 1588 GO terms. The GO terms were then categorized into the three main GO classes: biological processes, cellular components, and molecular functions. The targets of differentially expressed miRNAs were enriched in 43 GO terms, including 22 biological processes, 10 cellular components, and 11 molecular functions (p < 0.05). Many identified processes were associated with organism responses to the environmental stresses, including regulation of nematode larval development, response to nematode, and response to flooding. Our results will enhance the study and application of developing new cotton cultivars for nematode resistance.

Study of expressions of miRNAs in the spikelets based on their spatial location on panicle in rice cultivars provided insight into their influence on grain development

Development of rice cultivars bearing numerous spikelets by breeding approach to increase the yearly production of rice to approximately 800 million metric tons to feed the ever increasing population of the world accompanies poor grain filling in the inferior spikelets preventing achievement of the yield potential. As the initial stages of caryopses development are of much importance for grain filling, spatio-temporal expressions of the miRNAs were studied during these periods in the spikelets of a compact-panicle rice cultivar, Oryza sativa cv. Mahalaxmi, bearing numerous spikelets per panicle to understand the reason of poor grain filling at the level of the initial biochemical events. Differential expression of several known miRNAs between the superior and inferior spikelets suggested great difference in metabolism related to grain filling in the spikelets based on their spatial location on compact panicle. Expressions of five known and four novel miRNAs were validated by Northern. Their targets included the enzymes directly involved in starch biosynthesis like sucrose synthase, starch synthase and pullulanase, besides others. Spatio-temporal expression studies of these miRNAs in the spikelets of Mahalaxmi revealed a pattern of mostly a greater expression in the inferior spikelets compared with the superior ones concomitant with an inverse expression of the target genes, which was not observed in the lax-panicle cultivar Upahar. The study thus revealed that the grain filling in rice is greatly regulated by miRNAs, and these miRNAs or their target genes could be considered for biotechnological interventions for improving grain filling in the rice cultivars of interest.

Integrated miRNA-mRNA analysis reveals the roles of miRNAs in the replanting benefit of Achyranthes bidentata roots

The yield and quality of the medicinal plant Achyranthes bidentata can be increased when it is replanted into a field cultivated previously with the same crop, however, fundamental aspects of its biology (so-called "replanting benefit") still remain to be elucidated. miRNAs are sRNA molecules involved in the post-transcriptional regulation of gene expression in plant biological processes. Here, 267 conserved and 36 novel miRNAs were identified in A. bidentata roots. We compared the miRNA content of the roots (R1) from first-year planting with that of the roots (R2) of second-year replanting, and screened 21 differentially expressed (DE) miRNAs. Based on in silico functional analysis, integrated miRNA-mRNA datasets allowed the identification of 10 miRNA-target family modules, which might participate in the benefit. The expression profiles of the miRNA-target modules were potentially correlated with the presence of the replanting benefit. The indication was that the miRNA-responsive continuous monoculture could reprogram miRNA-mRNA expression patterns, which possibly promote the root growth and development, enhance its transport activity and strengthen its tolerance to various stresses, thereby improving A. bidentata productivity as observed in the replanting benefit. Our study provides basic data for further research on the molecular mechanisms of the benefit in A. bidentata.

MicroRNAs expression dynamics reveal post-transcriptional mechanisms regulating seed development in Phaseolus vulgaris L

The knowledge on post-transcriptional regulation mechanisms implicated in seed development (SD) is still limited, particularly in one of the most consumed grain legumes, Phaseolus vulgaris L. We explore for the first time the miRNA expression dynamics in P. vulgaris developing seeds. Seventy-two known and 39 new miRNAs were found expressed in P. vulgaris developing seeds. Most of the miRNAs identified were more abundant at 10 and 40 days after anthesis, suggesting that late embryogenesis/early filling and desiccation were SD stages in which miRNA action is more pronounced. Degradome analysis and target prediction identified targets for 77 expressed miRNAs. While several known miRNAs were predicted to target HD-ZIP, ARF, SPL, and NF-Y transcription factors families, most of the predicted targets for new miRNAs encode for functional proteins. MiRNAs-targets expression profiles evidenced that these miRNAs could tune distinct seed developmental stages. MiRNAs more accumulated at early SD stages were implicated in regulating the end of embryogenesis, postponing the seed maturation program, storage compound synthesis and allocation. MiRNAs more accumulated at late SD stages could be implicated in seed quiescence, desiccation tolerance, and longevity with still uncovered roles in germination. The miRNAs herein described represent novel P. vulgaris resources with potential application in future biotechnological approaches to modulate the expression of genes implicated in legume seed traits with impact in horticultural production systems.

Computational methods for annotation of plant regulatory non-coding RNAs using RNA-seq

Plant transcriptome encompasses numerous endogenous, regulatory non-coding RNAs (ncRNAs) that play a major biological role in regulating key physiological mechanisms. While studies have shown that ncRNAs are extremely diverse and ubiquitous, the functions of the vast majority of ncRNAs are still unknown. With ever-increasing ncRNAs under study, it is essential to identify, categorize and annotate these ncRNAs on a genome-wide scale. The use of high-throughput RNA sequencing (RNA-seq) technologies provides a broader picture of the non-coding component of transcriptome, enabling the comprehensive identification and annotation of all major ncRNAs across samples. However, the detection of known and emerging class of ncRNAs from RNA-seq data demands complex computational methods owing to their unique as well as similar characteristics. Here, we discuss major plant endogenous, regulatory ncRNAs in an RNA sample followed by computational strategies applied to discover each class of ncRNAs using RNA-seq. We also provide a collection of relevant software packages and databases to present a comprehensive bioinformatics toolbox for plant ncRNA researchers. We assume that the discussions in this review will provide a rationale for the discovery of all major categories of plant ncRNAs.

Regulatory protein genes and microRNAs in response to selenium stimuli in Pueraria lobata (Willd.) Ohwi

Regulatory protein genes and microRNAs (miRNAs) play important roles in response to abiotic and biotic stress, and the biosynthesis of secondary metabolites in plants. However, their responses to selenium (Se) stimuli have not been comprehensively studied in Pueraria lobata (Willd.) Ohwi, a selenocompound-rich medicinal and edible plant. In this study, we identified a total of 436/556/1161/624 transcription factors, 134/157/308/172 transcriptional regulators, and 341/456/250/518 protein kinases, which were co-expressed with at least one selenocompound-related structural gene/sulfate transporter or phosphate transporter/reactive oxygen species (ROS) scavenging structural gene/isoflavone-related structural gene, respectively. Then, we identified a total of 87 expressed miRNAs by Se disposure, in which 11 miRNAs, including miR171f-3p, miR390b-3P, miR-N111b, miR-N118, miR-N30, miR-N38-3P, miR-N61a, miR-N61b, miR-N80-3p, miR-N84-3P, and miR-N90.2-3P, were significantly upregulated. We also identified a total of 1172 target genes for the 87 expressed miRNAs. Gene Ontology enrichment analysis of these target genes showed that regulation of transcription, DNA-templated, integral component of membrane, nucleus, ATP binding, and plasma membrane are the top five subclassifications. Finally, we revealed that 5 miRNAs targeted 10 regulatory protein genes, which are highly correlated with at least one selenocompound-related structural gene or transporter gene; 5 miRNAs targeted 10 regulatory protein genes, which are highly correlated with at least one ROS scavenging structural gene; and 5 miRNAs targeted 9 regulatory protein genes, which are potentially involved in the isoflavone biosynthesis. Overall, the study provides us the comprehensive insight into the roles of regulatory proteins and miRNAs in response to Se stimuli in P. lobata.

Genome wide in-silico miRNA and target network prediction from stress responsive Horsegram (Macrotyloma uniflorum) accessions

Horsegram (Macrotyloma uniflorum (Lam.) Verdc.) is a drought hardy food and fodder legume of Indo-African continents with diverse germplasm sources demonstrating alternating mechanisms depicting contrasting adaptations to different climatic zones. Tissue specific expression of genes contributes substantially to location specific adaptations. Regulatory networks of such adaptive genes are elucidated for downstream translational research. MicroRNAs are small endogenous regulatory RNAs which alters the gene expression profiles at a particular time and type of tissue. Identification of such small regulatory RNAs in low moisture stress hardy crops can help in cross species transfer and validation confirming stress tolerance ability. This study outlined prediction of conserved miRNAs from transcriptome shotgun assembled sequences and EST sequences of horsegram. We could validate eight out of 15 of the identified miRNAs to demonstrate their role in deficit moisture stress tolerance mechanism of horsegram variety Paiyur1 with their target networks. The putative mumiRs were related to other food legumes indicating the presence of gene regulatory networks. Differential miRNA expression among drought specific tissues indicted the probable energy conservation mechanism. Targets were identified for functional characterization and regulatory network was constructed to find out the probable pathways of post-transcriptional regulation. The functional network revealed mechanism of biotic and abiotic stress tolerance, energy conservation and photoperiod responsiveness.

Computational Identification and Comparative Analysis of Conserved miRNAs and Their Putative Target Genes in the Juglans regia and J. microcarpa Genomes

MicroRNAs (miRNAs) are important factors for the post-transcriptional regulation of protein-coding genes in plants and animals. They are discovered either by sequencing small RNAs or computationally. We employed a sequence-homology-based computational approach to identify conserved miRNAs and their target genes in Persian (English) walnut, Juglans regia, and its North American wild relative, J. microcarpa. A total of 119 miRNA precursors (pre-miRNAs) were detected in the J. regia genome and 121 in the J. microcarpa genome and miRNA target genes were predicted and their functional annotations were performed in both genomes. In the J. regia genome, 325 different genes were targets; 87.08% were regulated by transcript cleavage and 12.92% by translation repression. In the J. microcarpa genome, 316 different genes were targets; 88.92% were regulated by transcript cleavage and 11.08% were regulated by translation repression. Totals of 1.3% and 2.0% of all resistance gene analogues (RGA) and 2.7% and 2.6% of all transcription factors (TFs) were regulated by miRNAs in the J. regia and J. microcarpa genomes, respectively. Juglans genomes evolved by a whole genome duplication (WGD) and consist of eight pairs of fractionated homoeologous chromosomes. Within each pair, the chromosome that has more genes with greater average transcription also harbors more pre-miRNAs and more target genes than its homoeologue. While only minor differences were detected in pre-miRNAs between the J. regia and J. microcarpa genomes, about one-third of the pre-miRNA loci were not conserved between homoeologous chromosome within each genome. Pre-miRNA and their corresponding target genes showed a tendency to be collocated within a subgenome.

Identification of candidate microRNAs from Ostreid herpesvirus-1 (OsHV-1) and their potential role in the infection of Pacific oysters (Crassostrea gigas)

Oyster production is an economic activity of great interest worldwide. Recently, oysters have been suffering significant mortalities from OsHV-1infection, which has resulted in substantial economic loses in several countries around the world. Understanding viral pathogenicity mechanisms is of central importance for the establishment of disease control measures. Thus, the present work aimed to identify and characterize miRNAs from OsHV-1 as well as to predict their target transcripts in the virus and the host. OsHV-1 genome was used for the in silico discovery of pre-miRNAs. Subsequently, viral and host target transcripts of the OsHV-1 miRNAs were predicted according to the base pairing interaction between mature miRNAs and mRNA 3' untranslated regions (UTRs). Six unique pre-miRNAs were found in different regions of the viral genome, ranging in length from 85 to 172 nucleotides. A complex network of self-regulation of viral gene expression mediated by the miRNAs was identified. These sequences also seem to have a broad ability to regulate the expression of host immune-related genes, especially those associated with pathogen recognition. Our results suggest that OsHV-1 encodes miRNAs with important functions in the infection process, inducing self-regulation of viral transcripts, as well as affecting the regulation of Pacific oyster transcripts related to immunity. Understanding the molecular basis of host-pathogen interactions can help mitigate the recurrent events of oyster mass mortalities by OsHV-1 observed worldwide.

Identification and the potential involvement of miRNAs in the regulation of artemisinin biosynthesis in A. annua

Micro RNAs (miRNAs) play crucial regulatory roles in multiple biological processes. Recently they have garnered the attention for their strong influence on the secondary metabolite production in plants. Their role in the regulation of artemisinin (ART) biosynthesis is, however, not fully elucidated. ART is a potent anti-malarial compound recommended by WHO for the treatment of drug-resistant malaria. It is produced by Artemisia annua (A. annua). The lower in planta content of ART necessitates a deep understanding of regulatory mechanisms involved in the biosynthesis of this metabolite. In this study, using modern high throughput small RNA-sequencing by Illumina Nextseq 500 platform for identification and stem-loop RT PCR for validation, miRNAs were identified in the leaf sample of A. annua plant. Here, we report a total of 121 miRNAs from A. annua that target several important genes and transcription factors involved in the biosynthesis of ART. This study revealed the presence of some important conserved miRNA families, miR396, miR319, miR399, miR858, miR5083 and miR6111 not identified so far in A. annua. The expression patterns and correlation between miRNAs and their corresponding targets at different developmental stages of the plant using real-time PCR indicate that they may influence ART accumulation. These findings thus, open new possibilities for the rational engineering of the secondary metabolite pathways in general and ART biosynthesis in particular.

Genome-wide identification of microRNAs involved in the regulation of fruit ripening and climacteric stages in melon (Cucumis melo)

Fruit ripening is influenced by multiple plant hormones and the regulation of genes. However, studies on posttranscriptional regulators (e.g., miRNAs) of fruit growth and ripening are limited. We used miRNA sequencing and degradome methods to identify miRNAs and their target genes in melon (Cucumis melo cv. Hetao melon). A total of 61 conserved miRNAs and 36 novel miRNAs were identified from fruit growth, ripening, climacteric, and postclimacteric developmental stage samples, of which 32 conserved miRNAs were differentially expressed between developmental stage samples. Sixty-two target genes of 43 conserved miRNAs and 1 novel miRNA were identified from degradome sequencing. To further investigate miRNA influencing fruit ripening, transgenic melon plants overexpressing pre-cme-miR393 (cme-miR393-OE) were generated and characterized. The results showed that fruit ripening was delayed in cme-miR393-OE transgenic lines compared to nontransgenic fruits. The target of cme-miR393 was also identified, and the expression of CmAFB2 was repressed in transgenic plants. These results provide evidence that miRNA regulates melon fruit ripening and provide potential targets to improve the horticultural traits of melon fruit.

Identification of Known and Novel MicroRNAs in Raspberry Organs Through High-Throughput Sequencing

MicroRNAs (miRNAs) are a class of small endogenous RNAs that play important regulatory roles in plants by negatively affecting gene expression. Studies on the identification of miRNAs and their functions in various plant species and organs have significantly contributed to plant development research. In the current study, we utilized high-throughput sequencing to detect the miRNAs in the root, stem, and leaf tissues of raspberry (Rubus idaeus). A total of more than 35 million small RNA reads ranging in size from 18 to 35 nucleotides were obtained, with 147 known miRNAs and 542 novel miRNAs identified among the three organs. Sequence verification and the relative expression profiles of the six known miRNAs were investigated by stem-loop quantitative real-time PCR. Furthermore, the potential target genes of the known and novel miRNAs were predicted and subjected to Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes pathway annotation. Enrichment analysis of the GO-associated biological processes and molecular functions revealed that these target genes were potentially involved in a wide range of metabolic pathways and developmental processes. Moreover, the miRNA target prediction revealed that most of the targets predicted as transcription factor-coding genes are involved in cellular and metabolic processes. This report is the first to identify miRNAs in raspberry. The detected miRNAs were analyzed by cluster analysis according to their expression, which revealed that these conservative miRNAs are necessary for plant functioning. The results add novel miRNAs to the raspberry transcriptome, providing a useful resource for the further elucidation of the functional roles of miRNAs in raspberry growth and development.

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.