A Comprehensive Prescription for Plant miRNA Identification

Elucidating circRNA-miRNA-mRNA competing endogenous regulatory RNA network during leaf rust pathogenesis in wheat (Triticum aestivum L.)

Advancements in bioinformatic tools and breakthroughs in high throughput RNA sequencing have unveiled the potential role of non-coding RNAs in influencing the overall expression of disease-responsive genes. Owing to the increasing need to develop resilient crop varieties against environmental constraints, our study explores the functional relationship of various non-coding RNAs in wheat during leaf rust pathogenesis. MicroRNAs (miRNAs) and circular RNAs (circRNAs) were retrieved from SAGE and RNA-Seq libraries, respectively, in the susceptible (HD2329) and resistant (HD2329 + Lr28) wheat Near-Isogenic Lines (NILs). Here we explored the previously published circRNAs for their differential expression and correlated the data with the differentially expressed miRNAs (DEMs) through various in silico methods to acquire the target miRNAs of circRNAs and the downstream target mRNAs of miRNAs. Finally, a competing endogenous RNA (ceRNAs) regulatory network was constructed and validated through RT-qPCR method. We have identified the ceRNA regulatory network of four differentially expressed circRNAs (DECs) and five DEMs to highlight their crucial roles in the robust enhancement of the temporal expression profiles of five defense responsive genes (mRNAs) in wheat NILs against leaf rust infection. The study confirms the synergistic expression of circRNAs and mRNAs with an antagonistic correlation with the expression profile of the corresponding miRNAs. The vital role of leaf rust-resistant gene Lr28 has also been highlighted for driving the efficiency of the circRNAs to upregulate target gene expression. Thus, understanding the circRNA-miRNA-target gene interaction during leaf rust pathogenesis can help to identify stress-specific regulatory biomarkers to enhance defense responses in wheat for improved resilience through multi-omics integration of transcriptomics, proteomics and metabolomics.

Plant miR6262 Modulates the Expression of Metabolic and Thermogenic Genes in Human Hepatocytes and Adipocytes

BACKGROUND

Edible plants have been linked to the mitigation of metabolic disturbances in liver and adipose tissue, including the decrease of lipogenesis and the enhancement of lipolysis and adipocyte browning. In this context, plant microRNAs could be key bioactive molecules underlying the cross-kingdom beneficial effects of plants. This study sought to explore the impact of plant-derived microRNAs on the modulation of adipocyte and hepatocyte genes involved in metabolism and thermogenesis.

METHODS

Plant miR6262 was selected as a candidate from miRBase for the predicted effect on the regulation of human metabolic genes. Functional validation was conducted after transfection with plant miRNA mimics in HepG2 hepatocytes exposed to free fatty acids to mimic liver steatosis and hMADs cells differentiated into brown-like adipocytes.

RESULTS

miR6262 decreases the expression of the predicted target RXRA in the fatty acids-treated hepatocytes and in brown-like adipocytes and affects the expression profile of critical genes involved in metabolism and thermogenesis, including PPARA, G6PC, SREBF1 (hepatocytes) and CIDEA, CPT1M and PLIN1 (adipocytes). Nevertheless, plant miR6262 mimic transfections did not decrease hepatocyte lipid accumulation or stimulate adipocyte browning.

CONCLUSIONS

these findings suggest that plant miR6262 could have a cross-kingdom regulation relevance through the modulation of human genes involved in lipid and glucose metabolism and thermogenesis in adipocytes and hepatocytes.

"Millet Models" for harnessing nuclear factor-Y transcription factors to engineer stress tolerance in plants: current knowledge and emerging paradigms

MAIN CONCLUSION

The main purpose of this review is to shed light on the role of millet models in imparting climate resilience and nutritional security and to give a concrete perspective on how NF-Y transcription factors can be harnessed for making cereals more stress tolerant. Agriculture faces significant challenges from climate change, bargaining, population, elevated food prices, and compromises with nutritional value. These factors have globally compelled scientists, breeders, and nutritionists to think of some options that can combat the food security crisis and malnutrition. To address these challenges, mainstreaming the climate-resilient and nutritionally unparalleled alternative crops like millet is a key strategy. The C4 photosynthetic pathway and adaptation to low-input marginal agricultural systems make millets a powerhouse of important gene and transcription factor families imparting tolerance to various kinds of biotic and abiotic stresses. Among these, the nuclear factor-Y (NF-Y) is one of the prominent transcription factor families that regulate diverse genes imparting stress tolerance. The primary purpose of this article is to shed light on the role of millet models in imparting climate resilience and nutritional security and to give a concrete perspective on how NF-Y transcription factors can be harnessed for making cereals more stress tolerant. Future cropping systems could be more resilient to climate change and nutritional quality if these practices were implemented.

The microRNAs in the antennae of Apolygus lucorum (Hemiptera: Miridae): Expression properties and targets prediction

MicroRNAs (miRNAs) regulate gene expression and contribute to numerous physiological processes. However, little is known about the functions of miRNAs in insect chemosensation. In this study, nine small RNA libraries were constructed and sequenced from the antennae of nymphs, adult males, and adult females of Apolygus lucorum. In total, 399 (275 known and 124 novel) miRNAs were identified. miR-7-5p_1 was the most abundant miRNA. Altogether, 69,708 target genes related to biogenesis, membrane, and binding activities were predicted. In particular, 15 miRNAs targeted 16 olfactory genes. Comparing the antennae of nymphs and adult males and females, 94 miRNAs were differentially expressed. Alternatively, a subset of differentially expressed miRNAs was verified by qPCR, supporting the reliability of the sequencing results. This study provides a global miRNA transcriptome for the antennae of A. lucorum and valuable information for further investigations of the functions of miRNAs in the regulation of chemosensation.

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.

Stress-induced miRNAs isolated from wheat have a unique therapeutic potential in ultraviolet-stressed human keratinocyte cells

Increasing evidence supports the existence of cross-kingdom gene regulation. However, the therapeutic potential of stress-specific plant miRNAs and their role in UV-related pathologies in human tissue remain largely unexplored. The aim of this study was to investigate the therapeutic potential and mechanisms of action of stress-induced miRNA cocktails (SI-WmiRs) from Einkorn wheat (Triticum monococcum L.) on human keratinocyte (HaCaT) cells exposed to a high dose of UV-B radiation. We used a biofactory approach and irradiated wheatgrass with UV-C for 240 min to obtain the specific SI-WmiRs that wheat produces to recover from UV stress. We followed the plant with molecular and biochemical analyses and extracted our SI-WmiRs at the most appropriate time (0 h and 6 h after UV-C application). Then, we applied the SI-WmiR cocktail to HaCaT cells exposed to high-dose of UV-B radiation. Our results show that UV-B radiation induced lipid peroxidation and DNA damage, as demonstrated by increased malondialdehyde (MDA) levels and changes in the RAPD band profile, respectively. UV stress also impaired IL6/JAK2/STAT3 signalling and activated the inflammatory mediators IL6 and TNF-α in HaCaT cells, leading to significant induction of apoptotic cell death. We found that SI-WmiR transfection prevents lipid peroxidation and oxidative stress-related DNA damage by increasing antioxidant (CuZn-SOD, Mn-SOD) and DNA repair (EXO1, SMUG1 and XRCC3) gene expression. In addition, SI-WmiRs regulated IL6/JAK2/STAT3 signalling by reducing JAK2 and STAT3 gene expression and phosphorylated protein levels compared to the control treatments. Moreover, SI-WmiRs inhibited pro-apoptotic BAX, Caspase 3 and Caspase 8 gene expression and protein levels to prevent apoptosis of UV-stressed HaCaT cells. Our results demonstrate that stress-induced wheat miRNAs produced using a biofactory approach have strong potential as a novel and effective alternative therapy for UV stress-related skin damage.

Computational screening of miRNAs and their targets in saffron (Crocus sativus L.) by transcriptome mining

A robust workflow for the identification of miRNAs and their targets in saffron was developed. MicroRNA-mediated gene regulation in saffron is potentially involved in several biological processes, including the biosynthesis of highly valuable apocarotenoids. Saffron (Crocus sativus L.) is the most expensive spice in the world and a major source of apocarotenoids. Even though miRNAs (20-24 nt non-coding small RNAs) are important regulators of gene expression at transcriptional and post-transcriptional levels, their role in saffron has not been thoroughly investigated. As a result, a workflow for computational identification of miRNAs and their targets can be useful to uncover the regulatory networks underlying biological processes in this valuable plant. The efficiency of several assembly tools such as Trans-ABySS, Trinity, Bridger, rnaSPAdes, and EvidentialGene was evaluated based on both reference-based and reference-free metrics using transcriptome data. A reliable workflow for computational identification of miRNAs and their targets in saffron was described. The EvidentialGene was found to be the most efficient de novo transcriptome assembler for saffron as a complex triploid model, followed by the Trinity. In total, 66 miRNAs from 19 different families that target 2880 genes, including several transcription factors involved in the flowering transition, were identified. Three of the identified targets were involved in the terpenoids backbone biosynthesis. CsCCD and CsUGT genes involved in the apocarotenoids biosynthetic pathway were targeted by csa-miR156g and csa-miR156b-3p, revealing a unique post-transcriptional regulation dynamic in saffron. The identified miRNAs and their targets add to our understanding of the many biological roles of miRNAs in saffron and shed new light on the control of the apocarotenoid biosynthetic pathway in this valuable plant.

Characterization of microRNAs from neem (Azadirachta indica) and their tissue-specific expression study in leaves and stem

Neem (Azadirachta indica) is a very popular traditional medicinal plant used since ancient times to treat numerous ailments. MicroRNAs (miRNAs) are highly conserved, non-coding, short RNA molecules that play important regulatory roles in plant development and metabolism. In this study, deploying a high stringent genome-wide computational-based approach and following a set of strict filtering norms a total of 44 potential conserved neem miRNAs belonging to 21 families and their corresponding 48 potential target transcripts were identified. Important targets include Squamosa promoter binding protein-like proteins, NAC, Scarecrow proteins, Auxin response factor, and F-box proteins. A biological network has also been developed to understand the miRNA-mediated gene regulation using the minimum free energy (MFE) values of the miRNA-target interaction. Moreover, six selected miRNAs were reported to be involved in secondary metabolism in other plant species (miR156a, miR156l, miR160, miR164, miR171, miR395) were validated by qPCR and their tissue-specific differential expression pattern was observed in leaves and stem. Except for ain-miR395, all the other miRNAs were found overexpressed in the stem as compared to leaves. To the best of our knowledge, this is the first report of neem miRNAs and we believe the finding of the present study will be useful for the functional genomic study of medicinal plants.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-021-02839-z.

Small RNA, Transcriptome and Degradome Analysis of the Transgenerational Heat Stress Response Network in Durum Wheat

Heat stress is a major limiting factor of grain yield and quality in crops. Abiotic stresses have a transgenerational impact and the mechanistic basis is associated with epigenetic regulation. The current study presents the first systematic analysis of the transgenerational effects of post-anthesis heat stress in tetraploid wheat. Leaf physiological traits, harvest components and grain quality traits were characterized under the impact of parental and progeny heat stress. The parental heat stress treatment had a positive influence on the offspring for traits including chlorophyll content, grain weight, grain number and grain total starch content. Integrated sequencing analysis of the small RNAome, mRNA transcriptome and degradome provided the first description of the molecular networks mediating heat stress adaptation under transgenerational influence. The expression profile of 1771 microRNAs (733 being novel) and 66,559 genes was provided, with differentially expressed microRNAs and genes characterized subject to the progeny treatment, parental treatment and tissue-type factors. Gene Ontology and KEGG pathway analysis of stress responsive microRNAs-mRNA modules provided further information on their functional roles in biological processes such as hormone homeostasis, signal transduction and protein stabilization. Our results provide new insights on the molecular basis of transgenerational heat stress adaptation, which can be used for improving thermo-tolerance in breeding.

MicroRNA annotation in plants: current status and challenges

Last two decades, the studies on microRNAs (miRNAs) and the numbers of annotated miRNAs in plants and animals have surged. Herein, we reviewed the current progress and challenges of miRNA annotation in plants. Via the comparison of plant and animal miRNAs, we pinpointed out the difficulties on plant miRNA annotation and proposed potential solutions. In terms of recalling the history of methods and criteria in plant miRNA annotation, we detailed how the major progresses made and evolved. By collecting and categorizing bioinformatics tools for plant miRNA annotation, we surveyed their advantages and disadvantages, especially for ones with the principle of mimicking the miRNA biogenesis pathway by parsing deeply sequenced small RNA (sRNA) libraries. In addition, we summarized all available databases hosting plant miRNAs, and posted the potential optimization solutions such as how to increase the signal-to-noise ratio (SNR) in these databases. Finally, we discussed the challenges and perspectives of plant miRNA annotations, and indicated the possibilities offered by an all-in-one tool and platform according to the integration of artificial intelligence.

Structural and Functional Characteristics of miRNAs in Five Strategic Millet Species and Their Utility in Drought Tolerance

Millets are the strategic food crops in arid and drought-prone ecologies. Millets, by virtue of nature, are very well-adapted to drought conditions and able to produce sustainable yield. Millets have important nutrients that can help prevent micro-nutrient malnutrition. As a result of the adverse effect of climate change and widespread malnutrition, millets have attained a strategic position to sustain food and nutritional security. Although millets can adapt well to the drought ecologies where other cereals fail completely, the yield level is very low under stress. There is a tremendous opportunity to increase the genetic potential of millet crops in dry lands when the genetics of the drought-tolerance mechanism is fully explained. MicroRNAs (miRNAs) are the class of small RNAs that control trait expression. They are part of the gene regulation but little studied in millets. In the present study, novel miRNAs and gene targets were identified from the genomic resources of pearl millet, sorghum, foxtail millet, finger millet, and proso millet through in silico approaches. A total of 1,002 miRNAs from 280 families regulating 23,158 targets were identified using different filtration criteria in five millet species. The unique as well as conserved structural features and functional characteristics of miRNA across millets were explained. About 84 miRNAs were conserved across millets in different species combinations, which explained the evolutionary relationship of the millets. Further, 215 miRNAs controlling 155 unique major drought-responsive genes, transcription factors, and protein families revealed the genetics of drought tolerance that are accumulated in the millet genomes. The miRNAs regulating the drought stress through specific targets or multiple targets showed through a network analysis. The identified genes regulated by miRNA genes could be useful in developing functional markers and used for yield improvement under drought in millets as well as in other crops.

Domestication of High-Copy Transposons Underlays the Wheat Small RNA Response to an Obligate Pathogen

Plant genomes have evolved several evolutionary mechanisms to tolerate and make use of transposable elements (TEs). Of these, transposon domestication into cis-regulatory and microRNA (miRNA) sequences is proposed to contribute to abiotic/biotic stress adaptation in plants. The wheat genome is derived at 85% from TEs, and contains thousands of miniature inverted-repeat transposable elements (MITEs), whose sequences are particularly prone for domestication into miRNA precursors. In this study, we investigate the contribution of TEs to the wheat small RNA immune response to the lineage-specific, obligate powdery mildew pathogen. We show that MITEs of the Mariner superfamily contribute the largest diversity of miRNAs to the wheat immune response. In particular, MITE precursors of miRNAs are wide-spread over the wheat genome, and highly conserved copies are found in the Lr34 and QPm.tut-4A mildew resistance loci. Our work suggests that transposon domestication is an important evolutionary force driving miRNA functional innovation in wheat immunity.

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, characterization, and expression analysis of cowpea (Vigna unguiculata [L.] Walp.) miRNAs in response to cowpea severe mosaic virus (CPSMV) challenge

Cowpea miRNAs and Argonaute genes showed differential expression patterns in response to CPSMV challenge Several biotic stresses affect cowpea production and yield. CPSMV stands out for causing severe negative impacts on cowpea. Plants have two main induced immune systems. In the basal system (PTI, PAMP-triggered immunity), plants recognize and respond to conserved molecular patterns associated with pathogens (PAMPs). The second type (ETI, Effector-triggered immunity) is induced after plant recognition of specific factors from pathogens. RNA silencing is another important defense mechanism in plants. Our research group has been using biochemical and proteomic approaches to learn which proteins and pathways are involved and could explain why some cowpea genotypes are resistant whereas others are susceptible to CPSMV. This current study was conducted to determine the role of cowpea miRNA in the interaction between a resistant cowpea genotype (BRS-Marataoã) and CPSMV. Previously identified and deposited plant microRNA sequences were used to find out all possible microRNAs in the cowpea genome. This search detected 617 mature microRNAs, which were distributed in 89 microRNA families. Next, 4 out of these 617 miRNAs and their possible target genes that encode the proteins Kat-p80, DEAD-Box, GST, and SPB9, all involved in the defense response of cowpea to CPSMV, had their expression compared between cowpea leaves uninoculated and inoculated with CPSMV. Additionally, the differential expression of genes that encode the Argonaute (AGO) proteins 1, 2, 4, 6, and 10 is reported. In summary, the studied miRNAs and AGO 2 and AGO4 associated genes showed differential expression patterns in response to CPSMV challenge, which indicate their role in cowpea defense.

Epstein Barr Virus Associated Lymphomas and Epithelia Cancers in Humans

Epstein Barr virus (EBV) is a cosmopolitan oncogenic virus, infecting about 90% of the world's population and it is associated to tumors originating from both epithelia and hematopoietic cells. Transmission of the virus is mainly through oral secretions; however, transmission through organ transplantation and blood transfusion has been reported. In order to evade immune recognition, EBV establishes latent infection in B lymphocytes where it expresses limited sets of proteins called EBV transcription programs (ETPs), including six nuclear antigens (EBNAs), three latent membrane proteins (LMP), and untranslated RNA called EBV encoded RNA (EBER), shown to efficiently transform B cells into lymphoblastic cells. These programs undergo different patterns of expression which determine the occurrence of distinct types of latency in the pathogenesis of a particular tumor. Hematopoietic cell derived tumors include but not limited to Burkitt's lymphoma, Hodgkin lymphoma, post-transplant lymphoproliferative disorders, and natural killer (NK)/T cell lymphoma. EBV undergoes lytic infection in epithelia cells for amplification of the viral particle for transmission where it expresses lytic stage genes. However, for reasons yet to be unveiled, EBV switches from the expression of lytic stage genes to the expression of ETPs in epithelia cells. The expression of the ETPs lead to the transformation of epithelia cells into permanently proliferating cells, resulting in epithelia cell derived malignancies such as nasopharyngeal cancer, gastric cancer, and breast cancer. In this review, we have summarized the current updates on EBV associated epithelial and B cell-derived malignancies, and the role of EBV latency gene products in the pathogenesis of the cancers, and have suggested areas for future studies when considering therapeutic measures.

Analysis of miRNAs in Two Wheat Cultivars Infected With Puccinia striiformis f. sp. tritici

MicroRNAs are small RNAs that regulate gene expression in eukaryotes. In this study, we analyzed the small RNA profiles of two cultivars that exhibit different reactions to stripe rust infection: one susceptible, the other partially resistant. Using small RNA libraries prepared from the two wheat cultivars infected with stripe rust fungus (Puccinia striiformis f. sp. tritici), we identified 182 previously known miRNAs, 91 variants of known miRNAs, and 163 candidate novel wheat miRNAs. Known miRNA loci were usually copied in all three wheat sub-genomes, whereas novel miRNA loci were often specific to a single sub-genome. DESeq2 analysis of differentially expressed microRNAs revealed 23 miRNAs that exhibit cultivar-specific differences. TA078/miR399b showed cultivar-specific differential regulation in response to infection. Using different target prediction algorithms, 145 miRNAs were predicted to target wheat genes, while 69 miRNAs were predicted to target fungal genes. We also confirmed reciprocal expression of TA078/miR399b and tae-miR9664 and their target genes in different treatments, providing evidence for miRNA-mediated regulation during infection. Both known and novel miRNAs were predicted to target fungal genes, suggesting trans-kingdom regulation of gene expression. Overall, this study contributes to the current repository of wheat miRNAs and provides novel information on the yet-uncharacterized roles for miRNAs in the wheat-stripe rust pathosystem.

Plant microRNAs: biogenesis, gene silencing, web-based analysis tools and their use as molecular markers

MicroRNAs (miRNAs) are tiny (20-24 nt bp) regulatory non-protein-coding RNA molecules that have been extensively characterized and found important for many physiological and developmental processes. The miss-expression of miRNAs leads to various defects in plants. MicroRNAs repress gene expression by directing mRNA degradation or translational arrest. Several proteins such as PP43A, HYL1, DCL, HST are indispensable role players in promoting miRNA biogenesis in plants. During miRNA biogenesis, lariat RNAs are produced as by-products of pre-mRNA splicing which have a negative role in regulation of miRNA homeostasis. By acting as a decoy and by sequestering to the dicing complex, lariat RNA can prevent the processing of miRNAs. A number of bioinformatic tools with different methodologies are available to identify and validate miRNAs and their targets. Many miRNAs have been reported in different crops for different traits; however, no reports are available on their use in plant breeding. Recently, researchers have developed trait specific miRNA-based molecular markers (miRNA-SSRs/SNP) for many quantitative traits in different plant species. In the future, these molecular markers can be used for plant breeding programs. In this review, a comprehensive up-to-date information is provided on the bioinformatic tools used for analysis of plant miRNAs and their targets, the number of miRNAs, their biogenesis, gene silencing mechanism and miRNA-based molecular markers.

Transcriptome-wide identification and characterization of microRNAs responsive to phosphate starvation in Populus tomentosa

miRNAs (microRNAs) are ~ 21-nt non-coding small RNAs (sRNAs) that play crucial regulatory roles in plant biotic and abiotic stress responses. Phosphorus (Pi) deficiency constrains plant growth and reduces yields worldwide. To identify tree miRNAs and evaluate their functions in the response to low Pi, we identified 261 known and 31 candidate novel miRNA families from three sRNA libraries constructed from Populus tomentosa subjected to sufficient or Pi deficiency condition or to restoration of a sufficient Pi level after Pi deficiency. Pi deficiency resulted in significant changes in the abundance of TPM (transcript per million) of 65 known and 3 novel miRNAs. Interestingly, four miRNAs responsive to low N-miR167, miR394, miR171, and miR857-were found to be involved in the response to low Pi. Thirty-five known and one novel miRNAs responded dynamically to Pi fluctuations, suggesting their involvement in the response to Pi deficiency. miRNA clusters comprising 36 miRNAs were identified in 10 chromosomes. Intriguingly, nine pairs of sense and antisense miRNAs transcribed from the same loci were detected in P. tomentosa, which is the first such report in woody plants. Moreover, target genes of the known miRNAs and novel miRNA candidates with significantly changed abundance were predicted, and their functions were annotated. Degradome sequencing supported the identified targets of miRNAs in P. tomentosa. These findings will enhance our understanding of universal and specific molecular regulatory mechanisms of trees under nutrition stress and may facilitate improvement of the Pi utilization efficiency of woody plants.

Chromosome-based survey sequencing reveals the genome organization of wild wheat progenitor Triticum dicoccoides

Wild emmer wheat (Triticum turgidum ssp. dicoccoides) is the progenitor of wheat. We performed chromosome-based survey sequencing of the 14 chromosomes, examining repetitive sequences, protein-coding genes, miRNA/target pairs and tRNA genes, as well as syntenic relationships with related grasses. We found considerable differences in the content and distribution of repetitive sequences between the A and B subgenomes. The gene contents of individual chromosomes varied widely, not necessarily correlating with chromosome size. We catalogued candidate agronomically important loci, along with new alleles and flanking sequences that can be used to design exome sequencing. Syntenic relationships and virtual gene orders revealed several small-scale evolutionary rearrangements, in addition to providing evidence for the 4AL-5AL-7BS translocation in wild emmer wheat. Chromosome-based sequence assemblies contained five novel miRNA families, among 59 families putatively encoded in the entire genome which provide insight into the domestication of wheat and an overview of the genome content and organization.

High-throughput sequencing of small RNAs revealed the diversified cold-responsive pathways during cold stress in the wild banana (Musa itinerans)

BACKGROUND

Cold stress is one of the most severe abiotic stresses affecting the banana production. Although some miRNAs have been identified, little is known about the role of miRNAs in response to cold stress in banana, and up to date, there is no report about the role of miRNAs in the response to cold stress in the plants of the cultivated or wild bananas.

RESULT

Here, a cold-resistant line wild banana (Musa itinerans) from China was used to profile the cold-responsive miRNAs by RNA-seq during cold stress. Totally, 265 known mature miRNAs and 41 novel miRNAs were obtained. Cluster analysis of differentially expressed (DE) miRNAs indicated that some miRNAs were specific for chilling or 0 °C treated responses, and most of them were reported to be cold-responsive; however, some were seldom reported to be cold-responsive in response to cold stress, e.g., miR395, miR408, miR172, suggesting that they maybe play key roles in response to cold stress. The GO and KEGG pathway enrichment analysis of DE miRNAs targets indicated that there existed diversified cold-responsive pathways, and miR172 was found likely to play a central coordinating role in response to cold stress, especially in the regulation of CK2 and the circadian rhythm. Finally, qPCR assays indicated the related targets were negatively regulated by the tested DE miRNAs during cold stress in the wild banana.

CONCLUSIONS

In this study, the profiling of miRNAs by RNA-seq in response to cold stress in the plants of the wild banana (Musa itinerans) was reported for the first time. The results showed that there existed diversified cold-responsive pathways, which provided insight into the roles of miRNAs during cold stress, and would be helpful for alleviating cold stress and cold-resistant breeding in bananas.

Deep sequencing and in silico analysis of small RNA library reveals novel miRNA from leaf Persicaria minor transcriptome

In current era, majority of microRNA (miRNA) are being discovered through computational approaches which are more confined towards model plants. Here, for the first time, we have described the identification and characterization of novel miRNA in a non-model plant, Persicaria minor (P. minor) using computational approach. Unannotated sequences from deep sequencing were analyzed based on previous well-established parameters. Around 24 putative novel miRNAs were identified from 6,417,780 reads of the unannotated sequence which represented 11 unique putative miRNA sequences. PsRobot target prediction tool was deployed to identify the target transcripts of putative novel miRNAs. Most of the predicted target transcripts (mRNAs) were known to be involved in plant development and stress responses. Gene ontology showed that majority of the putative novel miRNA targets involved in cellular component (69.07%), followed by molecular function (30.08%) and biological process (0.85%). Out of 11 unique putative miRNAs, 7 miRNAs were validated through semi-quantitative PCR. These novel miRNAs discoveries in P. minor may develop and update the current public miRNA database.

A large-scale multiomics analysis of wheat stem solidness and the wheat stem sawfly feeding response, and syntenic associations in barley, Brachypodium, and rice

The wheat stem sawfly (WSS), Cephus cinctus Norton (Hymenoptera: Cephidae), is an important pest of wheat and other cereals, threatening the quality and quantity of grain production. WSS larvae feed and develop inside the stem where they are protected from the external environment; therefore, pest management strategies primarily rely on host plant resistance. A major locus on the long arm of wheat chromosome 3B underlies most of the variation in stem solidness; however, the impact of stem solidness on WSS feeding has not been completely characterized. Here, we used a multiomics approach to examine the response to WSS in both solid- and semi-solid-stemmed wheat varieties. The combined transcriptomic, proteomic, and metabolomic data revealed that two important molecular pathways, phenylpropanoid and phosphate pentose, are involved in plant defense against WSS. We also detected a general downregulation of several key defense transcripts, including those encoding secondary metabolites such as DIMBOA, tricetin, and lignin, which suggested that the WSS larva might interfere with plant defense. We comparatively analyzed the stem solidness genomic region known to be associated with WSS tolerance in wild emmer, durum, and bread wheats, and described syntenic regions in the close relatives barley, Brachypodium, and rice. Additionally, microRNAs identified from the same genomic region revealed potential regulatory pathways associated with the WSS response. We propose a model outlining the molecular responses of the WSS–wheat interactions. These findings provide insight into the link between stem solidness and WSS feeding at the molecular level.

Plant Small Non-coding RNAs and Their Roles in Biotic Stresses

Non-coding RNAs (ncRNAs) have emerged as critical components of gene regulatory networks across a plethora of plant species. In particular, the 20-30 nucleotide small ncRNAs (sRNAs) play important roles in mediating both developmental processes and responses to biotic stresses. Based on variation in their biogenesis pathways, a number of different sRNA classes have been identified, and their specific functions have begun to be characterized. Here, we review the current knowledge of the biogenesis of the primary sRNA classes, microRNA (miRNA) and small nuclear RNA (snRNA), and their respective secondary classes, and discuss the roles of sRNAs in plant-pathogen interactions. sRNA mobility between species is also discussed along with potential applications of sRNA-plant-pathogen interactions in crop improvement technologies.

Assembly and Annotation of Transcriptome Provided Evidence of miRNA Mobility between Wheat and Wheat Stem Sawfly

Wheat Stem Sawfly (WSS), Cephus Cinctus Norton (Hymenoptera: Cephidae), is one of the most important pests, causing yield and economic losses in wheat and barley. The lack of information about molecular mechanisms of WSS for defeating plant's resistance prevents application of effective pest control strategies therefore, it is essential to identify the genes and their regulators behind WSS infestations. Long non-coding RNAs (lncRNAs) and microRNAs (miRNAs) are recognized with their regulatory functions on gene expression, tuning protein production by controlling transcriptional and post-transcriptional activities. A transcriptome-guided approach was followed in order to identify miRNAs, lncRNAs, and mRNA of WSS, and their interaction networks. A total of 1,893 were presented here as differentially expressed between larva and adult WSS insects. There were 11 miRNA families detected in WSS transcriptome. Together with the annotation of 1,251 novel mRNAs, the amount of genetic information available for WSS was expanded. The network between WSS miRNAs, lncRNAs, and mRNAs suggested miRNA-mediated regulatory roles of lncRNAs as competing endogenous RNAs. In the light of the previous evidence that small RNA molecules of a pathogen could suppress the immune response of host plant, we analyzed the putative interactions between larvae and wheat at the miRNA level. Overall, this study provides a profile of larva and adult WSS life stages in terms of coding and non-coding elements. These findings also emphasize the potential roles of wheat and larval miRNAs in wheat resistance to infestation and in the suppression of resistance which is critical for the development of effective pest control strategies.

RNA Sequencing and Co-expressed Long Non-coding RNA in Modern and Wild Wheats

There is an urgent need for the improvement of drought-tolerant bread and durum wheat. The huge and complex genome of bread wheat (BBAADD genome) stands as a vital obstruction for understanding the molecular mechanism underlying drought tolerance. However, tetraploid wheat (Triticum turgidum ssp., BBAA genome) is an ancestor of modern bread wheat and offers an important model for studying the drought response due to its less complex genome. Additionally, several wild relatives of tetraploid wheat have already shown a significant drought tolerance. We sequenced root transcriptome of three tetraploid wheat varieties with varying stress tolerance profiles, and built differential expression library of their transcripts under control and drought conditions. More than 5,000 differentially expressed transcripts were identified from each genotype. Functional characterization of transcripts specific to drought-tolerant genotype, revealed their association with osmolytes production and secondary metabolite pathways. Comparative analysis of differentially expressed genes and their non-coding RNA partners, long noncoding RNAs and microRNAs, provided valuable insight to gene expression regulation in response to drought stress. LncRNAs as well as coding transcripts share similar structural features in different tetraploid species; yet, lncRNAs slightly differ from coding transcripts. Several miRNA-lncRNA target pairs were detected as differentially expressed in drought stress. Overall, this study suggested an important pool of transcripts where their manipulations confer a better performance of wheat varieties under drought stress.

MicroRNAs in model and complex organisms

Non-coding RNAs such as microRNAs (miRNAs) are very tiny ribonucleotides having an essential role in gene regulation at both post-transcriptional and translational levels. They are very conserved and expressed in worms, flies, plants, and mammals in a sequence-specific manner. Furthermore, it is now possible to clone miRNAs using the new genome editing tool CRISPR/cas9, which shows benefit in control of untargeted effect. In this special issue, we tried to cover researches associated with functional roles of miRNAs accross model and complex organisms.

The 25-26 nt Small RNAs in Phytophthora parasitica Are Associated with Efficient Silencing of Homologous Endogenous Genes

Small RNAs (sRNAs) are important non-coding RNA regulators, playing key roles in developmental regulation, transposon suppression, environmental response, host-pathogen interaction and other diverse biological processes. However, their roles in oomycetes are poorly understood. Here, we performed sRNA sequencing and RNA sequencing of Phytophthora parasitica at stages of vegetative growth and infection of Arabidopsis roots to examine diversity and function of sRNAs in P. parasitica, a model hemibiotrophic oomycete plant pathogen. Our results indicate that there are two distinct types of sRNA-generating loci in P. parasitica genome, giving rise to clusters of 25-26 nt and 21 nt sRNAs, respectively, with no significant strand-biases. The 25-26 nt sRNA loci lie predominantly in gene-sparse and repeat-rich regions, and overlap with over 7000 endogenous gene loci. These overlapped genes are typically P. parasitica species-specific, with no homologies to the sister species P. infestans. They include approximately 40% RXLR effector genes, 50% CRN effector genes and some elicitor genes. The transcripts of most of these genes could not be detected at both the vegetative mycelium and infection stages as revealed by RNA sequencing, indicating that the 25-26 nt sRNAs are associated with efficient silencing of these genes. The 21 nt sRNA loci typically overlap with the exon regions of highly expressed genes, suggesting that the biogenesis of the 21 nt sRNAs may be dependent on the level of gene transcription and that these sRNAs do not mediate efficient silencing of homologous genes. Analyses of the published P. infestans sRNA and mRNA sequencing data consistently show that the 25-26 nt sRNAs, but not the 21 nt sRNAs, may mediate efficient gene silencing in Phytophthora.

Genome-wide analysis of microRNA targeting impacted by SNPs in cucumber genome

Background

microRNAs (miRNAs) are endogenous small RNAs that play important regulatory functions in plant development. Genetic variations in miRNAs sequences or their target-binding sites (microRNA-target interaction sites) can alter miRNA targets in animal and human. Whether these single nucleotide polymorphisms (SNPs) in plant are functional have not yet been determined.

Results

In this study, we constructed leaf, root, and stem-derived small libraries of cucumber (Cucumis sativus) line 9930 (cultivated China-group cucumber) and C. sativus var. hardwickii (wild India group cucumber). A total of 22 conserved miRNA families, nine less-conserved miRNA families, and 49 cucumber-specific miRNAs were identified in both line 9930 and hardwickii. We employed cucumber resequencing data to perform a genome-wide scan for SNPs in cucumber miRNA-target interaction sites, including miRNA mature sequences and miRNA-target binding sites. As a result, we identified a total of 19 SNPs in mature miRNA sequences and 113 SNPs in miRNA-target binding sites with the potential to affect miRNA-target interactions. Furthermore, we experimentally confirmed that these SNPs produced 14 9930-unique targets mRNAs and 15 hardwickii-unique targets mRNA for cucumber miRNAs. This is the first experimental validation of SNPs in miRNA-target interaction sites affecting miRNA-target binding in plants.

Conclusions

Our results indicate that SNPs can alter miRNA function and produce unique miRNA targets in cultivated and wild cucumbers. Therefore, miRNA-related SNPs may have played important in events that led to the agronomic differences between domestic and wild cucumber.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3665-y) contains supplementary material, which is available to authorized users.

Identification and characterization of durum wheat microRNAs in leaf and root tissues

MicroRNAs are a class of post-transcriptional regulators of plant developmental and physiological processes and responses to environmental stresses. Here, we present the study regarding the annotation and characterization of MIR genes conducted in durum wheat. We characterized the miRNAome of leaf and root tissues at tillering stage under two environmental conditions: irrigated with 100% (control) and 55% of evapotranspiration (early water stress). In total, 90 microRNAs were identified, of which 32 were classified as putative novel and species-specific miRNAs. In addition, seven microRNA homeologous groups were identified in each of the two genomes of the tetraploid durum wheat. Differential expression analysis highlighted a total of 45 microRNAs significantly differentially regulated in the pairwise comparisons leaf versus root. The miRNA families, miR530, miR395, miR393, miR5168, miR396 and miR166, miR171, miR319, and miR167, were the most expressed in leaves in comparison to roots. Putative microRNA targets were predicted for both five and three prime sequences derived from the stem-loop of the MIR gene. Gene ontology analysis showed significant overrepresented gene categories in microRNA targets belonging to transcription factors, phenylpropanoids, oxydases, and lipid binding-protein. This work represents one of the first genome wide characterization of MIR genes in durum wheat, identifying leaf and root tissue-specific microRNAs. This genomic identification of microRNAs together with the analysis of their expression profiles is a well-accepted starting point leading to a better comprehension of the role of MIR genes in the genus Triticum.

Characterization of miR061 and its target genes in grapevine responding to exogenous gibberellic acid

MicroRNAs (miRNAs), as an important growth regulator, are also involved in gibberellic acid (GA) signaling, revealing much relationship between miRNAs and GA in various plant responses. Grape is highly sensitive to GA₃, which plays a significant regulatory role in regulation of flower development, berry expansion, berry set, berry ripening, and seedlessness induction; further, it was found that grapevine miR061 (VvmiR061) is a GA₃ responsive miRNA. In this study, grapevine REV (VvREV) and HOX32 (VvHOX32), two target genes of VvmiR061, were predicted, verified, and cloned; homologous conservation was analyzed in various plants. The expression profiles of both VvmiR061 and its target genes (VvREV and VvHOX32) under GA₃ treatment were detected by qRT-PCR during grapevine flower and berry development. Results revealed that GA₃ treatment has upregulated the transcription of VvREV and VvHOX32, while it downregulated the expression of VvmiR061. The function of VvmiR061 in cleaving target genes VvREV and VvHOX32 was diminished by GA₃ treatment during flower developmental process. The results of this study exhibited the importance of VvmiR061 in regulating flower development and GA₃ signaling pathway and also contributed some to the knowledge of small RNA-mediated regulation in grape.