Small RNAs, emerging regulators critical for the development of horticultural traits

Research advances of coloring mechanism regulated by MicroRNAs in plants

In plants, microRNAs (miRNAs) are a class of important small RNAs involved in their growth and development, and play a very significant role in regulating their tissue coloring. In this paper, the mechanisms on miRNA regulation of plant coloring are mainly reviewed from three aspects: macroscopic physiological and molecular foundations related to tissue coloring, miRNA biosynthesis and function, and specific analysis of miRNA regulation studies on leaf color, flower color, fruit color, and other tissue color formation in plants. Furthermore, we also systematically summarize the miRNA regulatory mechanisms identified on pigments biosynthesis and color formation in plants, and the regulatory mechanisms of these miRNAs mentioned on the existing researches can be divided into four main categories: directly targeting the related transcription factors, directly targeting the related structural genes, directly targeting the related long noncoding RNAs (LncRNAs) and miRNA-mediated production of trans-acting small interfering RNAs (ta-siRNAs). Together, these research results aim to provide a theoretical reference for the in-depth study of plant coloring mechanism and molecular breeding study of related plants in the future.

Exogenous Application of dsRNA-Inducing Silencing of the Fusarium oxysporum Tup1 Gene and Reducing Its Virulence

Fusarium oxysporum is a widespread soil-borne fungal pathogen that can infect various plants, causing wilt and root rot diseases. The root rot disease of Atractylodes macrocephala caused by F. oxysporum is among the most serious diseases associated with continuous cropping, significantly hindering its sustainable development. In this study, we aimed to investigate the effect of exogenous application of double-stranded RNA (dsRNA) on silencing the F. oxysporum Tup1 gene to reduce its virulence and to evaluate its potential application in controlling root rot disease in A. macrocephala. The Tup1 gene was amplified from the F. oxysporum genome, and different lengths of Tup1-dsRNA were designed and synthesized. The uptake of dsRNA by the fungus was verified using Tup1-dsRNA labeled with fluorescein, and in vitro dsRNA treatment experiments were conducted to assess its impact on the growth and virulence of F. oxysporum. Additionally, Tup1-dsRNA was applied to the roots of A. macrocephala to evaluate its effectiveness in controlling root rot disease. The experimental results showed that F. oxysporum could effectively uptake exogenously applied Tup1-dsRNA, significantly reducing Tup1 gene expression. All lengths of Tup1-dsRNA inhibited fungal growth and caused morphological changes in the fungal hyphae. Further plant experiments and Reverse Transcription Quantitative Polymerase Chain Reaction (RT-qPCR) analysis indicated that Tup1-dsRNA treatment significantly reduced the incidence of root rot disease in A. macrocephala, which was supported by the reduction in peroxidase (POD) and catalase (CAT) enzyme activities, malondialdehyde (MDA) content, and proline (Pro) levels in treated root tissues. This study demonstrated that exogenous dsRNA could reduce the virulence of F. oxysporum by silencing the Tup1 gene and effectively mitigate the root rot disease it causes in A. macrocephala. The successful application of Tup1-dsRNA provided strong evidence for the potential of RNA interference (RNAi) technology in plant disease control. Future research could further optimize the design and application of dsRNA to enhance its practical value in agriculture.

The Integrated Analysis of miRNome and Degradome Sequencing Reveals the Regulatory Mechanisms of Seed Development and Oil Biosynthesis in Pecan (Carya illinoinensis)

Pecan seed oil is a valuable source of essential fatty acids and various bioactive compounds; however, the functions of microRNAs and their targets in oil biosynthesis during seed development are still unknown. Here, we found that the oil content increased rapidly in the three early stages in three cultivars, and that oleic acid was the predominant fatty acid component in the mature pecan embryos. We identified, analyzed, and validated the expression levels of miRNAs related to seed development and oil biosynthesis, as well as their potential target genes, using small RNA sequencing data from three stages (120, 135, and 150 days after flowering). During the seed development process, 365 known and 321 novel miRNAs were discovered. In total, 91 known and 181 novel miRNAs were found to be differentially expressed, and 633 target genes were further investigated. The expression trend analysis revealed that the 91 known miRNAs were classified into eight groups, approximately two-thirds of which were up-regulated, whereas most novel miRNAs were down-regulated. The qRT-PCR and degradome sequencing data were used to identify five miRNA- target pairs. Overall, our study provides valuable insights into the molecular regulation of oil biosynthesis in pecan seeds.

The Biosynthesis Process of Small RNA and Its Pivotal Roles in Plant Development

In the realm of plant biology, small RNAs (sRNAs) are imperative in the orchestration of gene expression, playing pivotal roles across a spectrum of developmental sequences and responses to environmental stressors. The biosynthetic cascade of sRNAs is characterized by an elaborate network of enzymatic pathways that meticulously process double-stranded RNA (dsRNA) precursors into sRNA molecules, typically 20 to 30 nucleotides in length. These sRNAs, chiefly microRNAs (miRNAs) and small interfering RNAs (siRNAs), are integral in guiding the RNA-induced silencing complex (RISC) to selectively target messenger RNAs (mRNAs) for post-transcriptional modulation. This regulation is achieved either through the targeted cleavage or the suppression of translational efficiency of the mRNAs. In plant development, sRNAs are integral to the modulation of key pathways that govern growth patterns, organ differentiation, and developmental timing. The biogenesis of sRNA itself is a fine-tuned process, beginning with transcription and proceeding through a series of processing steps involving Dicer-like enzymes and RNA-binding proteins. Recent advances in the field have illuminated the complex processes underlying the generation and function of small RNAs (sRNAs), including the identification of new sRNA categories and the clarification of their involvement in the intercommunication among diverse regulatory pathways. This review endeavors to evaluate the contemporary comprehension of sRNA biosynthesis and to underscore the pivotal role these molecules play in directing the intricate performance of plant developmental processes.

Identification of miRNA858 long-loop precursors in seed plants

MicroRNAs (miRNAs) are a class of nonprotein-coding short transcripts that provide a layer of post-transcriptional regulation essential to many plant biological processes. MiR858, which targets the transcripts of MYB transcription factors, can affect a range of secondary metabolic processes. Although miR858 and its 187-nt precursor have been well studied in Arabidopsis (Arabidopsis thaliana), a systematic investigation of miR858 precursors and their functions across plant species is lacking due to a problem in identifying the transcripts that generate this subclass. By re-evaluating the transcript of miR858 and relaxing the length cut-off for identifying hairpins, we found in kiwifruit (Actinidia chinensis) that miR858 has long-loop hairpins (1,100 to 2,100 nt), whose intervening sequences between miRNA generating complementary sites were longer than all previously reported miRNA hairpins. Importantly, these precursors of miR858 containing long-loop hairpins (termed MIR858L) are widespread in seed plants including Arabidopsis, varying between 350 and 5,500 nt. Moreover, we showed that MIR858L has a greater impact on proanthocyanidin and flavonol levels in both Arabidopsis and kiwifruit. We suggest that an active MIR858L-MYB regulatory module appeared in the transition of early land plants to large upright flowering plants, making a key contribution to plant secondary metabolism.

MicroRNA482/2118 is lineage-specifically involved in gibberellin signalling via the regulation of GID1 expression by targeting noncoding PHAS genes and subsequently instigated phasiRNAs

MicroRNA482/2118 (miR482/2118) is a 22-nt miRNA superfamily, with conserved functions in disease resistance and plant development. It usually instigates the production of phased small interfering RNAs (phasiRNAs) from its targets to expand or reinforce its silencing effect. Using a new high-quality reference genome sequence and comprehensive small RNA profiling, we characterized a newly evolved regulatory pathway of miR482/2118 in litchi. In this pathway, miR482/2118 cleaved a novel noncoding trans-acting gene (LcTASL1) and triggered phasiRNAs to regulate the expression of gibberellin (GA) receptor gene GIBBERELLIN INSENSITIVE DWARF1 (GID1) in trans; another trans-acting gene LcTASL2, targeted by LcTASL1-derived phasiRNAs, produced phasiRNAs as well to target LcGID1 to reinforce the silencing effect of LcTASL1. We found this miR482/2118-TASL-GID1 pathway was likely involved in fruit development, especially the seed development in litchi. In vivo construction of the miR482a-TASL-GID1 pathway in Arabidopsis could lead to defects in flower and silique development, analogous to the phenotype of gid1 mutants. Finally, we found that a GA-responsive transcription factor, LcGAMYB33, could regulate LcMIR482/2118 as a feedback mechanism of the sRNA-silencing pathway. Our results deciphered a lineage-specifically evolved regulatory module of miR482/2118, demonstrating the high dynamics of miR482/2118 function in plants.

sRNAminer: A multifunctional toolkit for next-generation sequencing small RNA data mining in plants

Small RNAs (sRNAs), found extensively in plants, play an essential role in plant growth and development. Although various sRNA analysis tools have been developed for plants, the use of most of them depends on programming and command-line environments, which is a challenge for many wet-lab biologists. Furthermore, current sRNA analysis tools mostly focus on the analysis of certain type of sRNAs and are resource-intensive, normally demanding an immense amount of time and effort to learn the use of numerous tools or scripts and assemble them into a workable pipeline to get the final results. Here, we present sRNAminer, a powerful stand-alone toolkit with a user-friendly interface that integrates all common functions for the analysis of three major types of plant sRNAs: microRNAs (miRNAs), phased small interfering RNAs (phasiRNAs), and heterochromatic siRNAs (hc-siRNAs). We constructed a curated or "golden" set of MIRNA and PHAS loci, which was used to assess the performance of sRNAminer in comparison to other existing tools. The results showed that sRNAminer outperformed these tools in multiple aspects, highlighting its functionality. In addition, to enable an efficient evaluation of sRNA annotation results, we developed Integrative Genomics Viewer (IGV)-sRNA, a modified genome browser optimized from IGV and we incorporated it as a functional module in sRNAminer. IGV-sRNA can display a wealth of sRNA-specific features, enabling a more comprehensive understanding of sRNA data. sRNAminer and IGV-sRNA are both platform-independent software that can be run under all operating systems. They are now freely available at https://github.com/kli28/sRNAminer and https://gitee.com/CJchen/IGV-sRNA.

Identification and Analysis of the MIR399 Gene Family in Grapevine Reveal Their Potential Functions in Abiotic Stress

MiR399 plays an important role in plant growth and development. The objective of the present study was to elucidate the evolutionary characteristics of the MIR399 gene family in grapevine and investigate its role in stress response. To comprehensively investigate the functions of miR399 in grapevine, nine members of the Vvi-MIR399 family were identified based on the genome, using a miRBase database search, located on four chromosomes (Chr 2, Chr 10, Chr 15, and Chr 16). The lengths of the Vvi-miR399 precursor sequences ranged from 82 to 122 nt and they formed stable stem-loop structures, indicating that they could produce microRNAs (miRNAs). Furthermore, our results suggested that the 2 to 20 nt region of miR399 mature sequences were relatively conserved among family members. Phylogenetic analysis revealed that the Vvi-MIR399 members of dicots (Arabidopsis, tomato, and sweet orange) and monocots (rice and grapevine) could be divided into three clades, and most of the Vvi-MIR399s were closely related to sweet orange in dicots. Promoter analysis of Vvi-MIR399s showed that the majority of the predicted cis-elements were related to stress response. A total of 66.7% (6/9) of the Vvi-MIR399 promoters harbored drought, GA, and SA response elements, and 44.4% (4/9) of the Vvi-MIRR399 promoters also presented elements involved in ABA and MeJA response. The expression trend of Vvi-MIR399s was consistent in different tissues, with the lowest expression level in mature and young fruits and the highest expression level in stems and young leaves. However, nine Vvi-MIR399s and four target genes showed different expression patterns when exposed to low light, high light, heat, cold, drought, and salt stress. Interestingly, a putative target of Vvi-MIR399 targeted multiple genes; for example, seven Vvi-MIR399s simultaneously targeted VIT_213s0067g03280.1. Furthermore, overexpression of Vvi_MIR399e and Vvi_MIR399f in Arabidopsis enhanced tolerance to drought compared with wild-type (WT). In contrast, the survival rate of Vvi_MIR399d-overexpressed plants were zero after drought stress. In conclusion, Vvi-MIR399e and Vvi-MIR399f, which are related to drought tolerance in grapevine, provide candidate genes for future drought resistance breeding.

Plant biomarkers as early detection tools in stress management in food crops: a review

MAIN CONCLUSION

Plant Biomarkers are objective indicators of a plant's cellular state in response to abiotic and biotic stress factors. They can be explored in crop breeding and engineering to produce stress-tolerant crop species. Global food production safely and sustainably remains a top priority to feed the ever-growing human population, expected to reach 10 billion by 2050. However, abiotic and biotic stress factors negatively impact food production systems, causing between 70 and 100% reduction in crop yield. Understanding the plant stress responses is critical for developing novel crops that can adapt better to various adverse environmental conditions. Using plant biomarkers as measurable indicators of a plant's cellular response to external stimuli could serve as early warning signals to detect stresses before severe damage occurs. Plant biomarkers have received considerable attention in the last decade as pre-stress indicators for various economically important food crops. This review discusses some biomarkers associated with abiotic and biotic stress conditions and highlights their importance in developing stress-resilient crops. In addition, we highlighted some factors influencing the expression of biomarkers in crop plants under stress. The information presented in this review would educate plant researchers, breeders, and agronomists on the significance of plant biomarkers in stress biology research, which is essential for improving plant growth and yield toward sustainable food production.

Integrative transcriptomic analysis unveils lncRNA-miRNA-mRNA interplay in tomato plants responding to Ralstonia solanacearum

Ralstonia solanacearum, a bacterial plant pathogen, poses a significant threat to tomato (Solanum lycopersicum) production through destructive wilt disease. While noncoding RNA has emerged as a crucial regulator in plant disease, its specific involvement in tomato bacterial wilt remains limited. Here, we conducted a comprehensive analysis of the transcriptional landscape, encompassing both mRNAs and noncoding RNAs, in a tomato resistant line ('ZRS_7') and a susceptible line ('HTY_9') upon R. solanacearum inoculation using high-throughput RNA sequencing. Differential expression (DE) analysis revealed significant alterations in 7506 mRNAs, 997 lncRNAs, and 69 miRNAs between 'ZRS_7' and 'HTY_9' after pathogen exposure. Notably, 4548 mRNAs, 367 lncRNAs, and 26 miRNAs exhibited genotype-specific responses to R. solanacearum inoculation. GO and KEGG pathway analyses unveiled the potential involvement of noncoding RNAs in the response to bacterial wilt disease, targeting receptor-like kinases, cell wall-related genes, glutamate decarboxylases, and other key pathways. Furthermore, we constructed a comprehensive competing endogenous RNA (ceRNA) network incorporating 13 DE-miRNAs, 30 DE-lncRNAs, and 127 DEGs, providing insights into their potential contributions to the response against bacterial inoculation. Importantly, the characterization of possible endogenous target mimics (eTMs) of Sly-miR482e-3p via VIGS technology demonstrated the significant impact of eTM482e-3p-1 silencing on tomato's sensitivity to R. solanacearum. These findings support the existence of an eTM482e-3p-1-Sly-miR482e-3p-NBS-LRRs network in regulating tomato's response to the pathogen. Collectively, our findings shed light on the intricate interactions among lncRNAs, miRNAs, and mRNAs as underlying factors in conferring resistance to R. solanacearum in tomato.

The PcMYB44-mediated miR397-PcLACs module regulates defence-induced lignification in pear resistance to fungal disease

Diseases caused by Alternaria alternata and Botryosphaeria dothidea diminish pear yield and quality, and restrict the pear agricultural industry. Lignification is a conserved mechanism for plant resistance against pathogen invasion. The regulatory mechanisms underlying defence-induced lignification in pear in response to fungal pathogen infection remain unknown. In this study, analysis of lignification level and lignin content in pear revealed that A. alternata and B. dothidea induced lignification, and transcriptomics showed that lignin biosynthesis was affected. To explore whether laccases (LACs) mediated by miR397 regulate lignification in pear, we investigated the role of PcmiR397 in repressing the expression of PcLACs using 5'-RNA ligase-mediated-RACE and co-transformation in tobacco. Opposite expression patterns for PcmiR397 and PcLAC target genes were observed in pear in response to pathogens. Transient transformation in pear demonstrated that silencing PcmiR397 and overexpressing a single PcLAC enhanced resistance to pathogens via lignin synthesis. To further reveal the mechanism underpinning the PcMIR397 response of pear to pathogens, the PcMIR397 promoter was analysed, and pMIR397-1039 was found to be inhibited by pathogen infection. The transcription factor PcMYB44 was up-regulated, and it bound to the PcMIR397 promoter and inhibited transcription following pathogen infection. The results demonstrate the role of PcmiR397-PcLACs in broad-spectrum resistance to fungal disease, and the potential role of PcMYB44 involved in the miR397-PcLAC module in regulating defence-induced lignification. The findings provide valuable candidate gene resources and guidance for molecular breeding to improve resistance to fungal disease in pear.

Cotton RSG2 Mediates Plant Resistance against Verticillium dahliae by miR482b Regulation

Cotton Verticillium wilt, mainly caused by Verticillium dahliae, has a serious impact on the yield and quality of cotton fiber. Many microRNAs (miRNAs) have been identified to participate in plant resistance to V. dahliae infection, but the exploration of miRNA's function mechanism in plant defense is needed. Here, we demonstrate that the ghr-miR482b-GhRSG2 module mediates cotton plant resistance to V. dahliae infection. Based on the mRNA degradation data and GUS fusion experiments, ghr-miR482b directedly bonds to GhRSG2 mRNA to lead to its degradation. The knockdown and overexpression of ghr-miR482b through virus-induced gene silencing strategies enhanced (decreased by 0.39-fold in disease index compared with the control) and weakened (increased by 0.46-fold) the plant resistance to V. dahliae, respectively. In addition, silencing GhRSG2 significantly increased (increased by 0.93-fold in disease index) the plant sensitivity to V. dahliae compared with the control plants treated with empty vector. The expression levels of two SA-related disease genes, GhPR1 and GhPR2, significantly decreased in GhRSG2-silenced plants by 0.71 and 0.67 times, respectively, and in ghr-miR482b-overexpressed (OX) plants by 0.59 and 0.75 times, respectively, compared with the control, whereas the expression levels of GhPR1 and GhPR2 were significantly increased by 1.21 and 2.59 times, respectively, in ghr-miR482b knockdown (KD) plants. In sum, the ghr-miR482b-GhRSG2 module participates in the regulation of plant defense against V. dahliae by inducing the expression of PR1 and PR2 genes.

MicroRNA156ab regulates apple plant growth and drought tolerance by targeting transcription factor MsSPL13

Drought stress substantially reduces the productivity of apple plants and severely restricts the development of the apple industry. Malus sieversii, a wild apple with excellent drought resistance, is a valuable wild resource for rootstock improvement of cultivated apple (Malus domestica). miRNAs and their targets play essential roles in plant growth and stress responses, but their roles in drought stress responses in apple are unknown. Here, we demonstrate that microRNA156ab is upregulated in M. sieversii in response to drought stress. Overexpressing msi-miR156ab promoted auxin accumulation, maintained the growth of apple plants, and increased plant resistance to osmotic stress. Antioxidant enzyme activities and proline contents were also increased in miR156ab-OE transgenic apple lines, which improved drought resistance. The SQUAMOSA PROMOTER BINDING PROTEIN-LIKE (SPL) transcription factor MsSPL13 is the target of msi-miR156ab, as demonstrated by 5-RACE and dual luciferase assays. Heterologous expression of MsSPL13 decreased auxin contents and inhibited growth in Arabidopsis (Arabidopsis thaliana) under normal and stress conditions. The activities of antioxidant enzymes were also suppressed in MsSPL13-OE transgenic Arabidopsis, reducing drought resistance. We showed that MsSPL13 regulates the expression of the auxin-related genes MsYUCCA5, PIN-FORMED7 (MsPIN7), and Gretchen Hagen3-5 (MsGH3-5) by binding to the GTAC cis-elements in their promoters, thereby regulating auxin metabolism. Finally, we demonstrated that the miR156ab-SPL13 module is involved in mediating the difference in auxin metabolism and stress responses between the M. sieversii and M26 (M. domestica) rootstocks. Overall, these findings reveal that the miR156ab-SPL13 module enhances drought stress tolerance in apples by regulating auxin metabolism and antioxidant enzyme activities.

miRador: a fast and precise tool for the prediction of plant miRNAs

Plant microRNAs (miRNAs) are short, noncoding RNA molecules that restrict gene expression via posttranscriptional regulation and function in several essential pathways, including development, growth, and stress responses. Accurately identifying miRNAs in populations of small RNA sequencing libraries is a computationally intensive process that has resulted in the misidentification of inaccurately annotated miRNA sequences. In recent years, criteria for miRNA annotation have been refined with the aim to reduce these misannotations. Here, we describe miRador, a miRNA identification tool that utilizes the most up-to-date, community-established criteria for accurate identification of miRNAs in plants. We combined target prediction and Parallel Analysis of RNA Ends (PARE) data to assess the precision of the miRNAs identified by miRador. We compared miRador to other commonly used miRNA prediction tools and found that miRador is at least as precise as other prediction tools while being substantially faster than other tools. miRador should be broadly useful for the plant community to identify and annotate miRNAs in plant genomes.

The Tumorigenic Role of Circular RNA-MicroRNA Axis in Cancer

Circular RNAs (circRNAs) are a class of endogenous RNAs that control gene expression at the transcriptional and post-transcriptional levels. Recent studies have increasingly demonstrated that circRNAs act as novel diagnostic biomarkers and promising therapeutic targets for numerous cancer types by interacting with other non-coding RNAs such as microRNAs (miRNAs). The miRNAs are presented as crucial risk factors and regulatory elements in cancer by regulating the expression of their target genes. Some miRNAs are derived from transposable elements (MDTEs) that can transfer their location to another region of the genome. Genetic interactions between miRNAs and circular RNAs can form complex regulatory networks with various carcinogenic processes that play critical roles in tumorigenesis and cancer progression. This review focuses on the biological regulation of the correlative axis among circular RNAs, miRNAs, and their target genes in various cancer types and suggests the biological importance of MDTEs interacting with oncogenic or tumor-suppressive circRNAs in tumor progression.

Coordinating Diverse Functions of miRNA and lncRNA in Fleshy Fruit

Non-coding RNAs play vital roles in the diverse biological processes of plants, and they are becoming key topics in horticulture research. In particular, miRNAs and long non-coding RNAs (lncRNAs) are receiving increased attention in fruit crops. Recent studies in horticulture research provide both genetic and molecular evidence that miRNAs and lncRNAs regulate biological function and stress responses during fruit development. Here, we summarize multiple regulatory modules of miRNAs and lncRNAs and their biological roles in fruit sets and stress responses, which would guide the development of molecular breeding techniques on horticultural crops.

Integrated analysis of transcriptome and microRNAs associated with exogenous calcium-mediated enhancement of hypoxic tolerance in cucumber seedlings (Cucumis sativus L.)

Plants often suffer from hypoxic stress due to flooding caused by extreme weather. Hypoxia usually leads to restricted oxygen supply and alters metabolic patterns from aerobic to anaerobic. Cucumber roots are fragile and highly sensitive to damage from hypoxic stress. The purpose of this study was to investigate the regulatory mechanism of exogenous calcium alleviating hypoxic stress in cucumber through transcriptome and small RNAs analysis. Three treatments were performed in this paper, including untreated-control (CK), hypoxic stress (H), and hypoxic stress + exogenous calcium treatment (H + Ca²⁺). A large number of differentially expressed genes (DEGs) were identified, 1,463 DEGs between CK vs H, 3,399 DEGs between H vs H + Ca²⁺, and 5,072 DEGs between CK vs H + Ca²⁺, respectively. KEGG analysis of DEGs showed that exogenous calcium could activate hormone signaling pathways (ethylene, ABA, IAA and cytokinin), transcription factors (MYB, MYB-related, bHLH, bZIP, and WRKY), calcium signaling and glycolysis pathway to mitigating hypoxic stress in cucumber seedlings. Additionally, miRNA and their target genes were detected and predicted between treatments. The target genes of these miRNAs revealed that auxin, cellulose synthase, and mitochondrial ribosomal related genes (Csa2G315390, Csa6G141390, Csa4G053280, and Csa6G310480) probably play in the improvement of the hypoxic tolerance of cucumber seedlings through exogenous calcium application. In short, our data adds new information to the mechanism of exogenous calcium mitigation of hypoxic stress injury in cucumber seedlings at transcriptional and post-transcriptional levels.

The roles of microRNAs in horticultural plant disease resistance

The development of the horticultural industry is largely limited by disease and excessive pesticide application. MicroRNAs constitute a major portion of the transcriptomes of eukaryotes. Various microRNAs have been recognized as important regulators of the expression of genes involved in essential biological processes throughout the whole life cycle of plants. Recently, small RNA sequencing has been applied to study gene regulation in horticultural plants. In this review, we summarize the current understanding of the biogenesis and contributions of microRNAs in horticultural plant disease resistance. These microRNAs may potentially be used as genetic resources for improving disease resistance and for molecular breeding. The challenges in understanding horticultural plant microRNA biology and the possibilities to make better use of these horticultural plant gene resources in the future are discussed in this review.

Genome-Wide Identification and Characterization of Argonaute, Dicer-like and RNA-Dependent RNA Polymerase Gene Families and Their Expression Analyses in Fragaria spp

In the growth and development of plants, some non-coding small RNAs (sRNAs) not only mediate RNA interference at the post-transcriptional level, but also play an important regulatory role in chromatin modification at the transcriptional level. In these processes, the protein factors Argonaute (AGO), Dicer-like (DCL), and RNA-dependent RNA polymerase (RDR) play very important roles in the synthesis of sRNAs respectively. Though they have been identified in many plants, the information about these gene families in strawberry was poorly understood. In this study, using a genome-wide analysis and a phylogenetic approach, 13 AGO, six DCL, and nine RDR genes were identified in diploid strawberry Fragaria vesca. We also identified 33 AGO, 18 DCL, and 28 RDR genes in octoploid strawberry Fragaria × ananassa, studied the expression patterns of these genes in various tissues and developmental stages of strawberry, and researched the response of these genes to some hormones, finding that almost all genes respond to the five hormone stresses. This study is the first report of a genome-wide analysis of AGO, DCL, and RDR gene families in Fragaria spp., in which we provide basic genomic information and expression patterns for these genes. Additionally, this study provides a basis for further research on the functions of these genes and some evidence for the evolution between diploid and octoploid strawberries.

Gma-miR408 Enhances Soybean Cyst Nematode Susceptibility by Suppressing Reactive Oxygen Species Accumulation

Soybean cyst nematode (SCN, Heterodera glycine) is a serious damaging disease in soybean worldwide, thus resulting in severe yield losses. MicroRNA408 (miR408) is an ancient and highly conserved miRNA involved in regulating plant growth, development, biotic and abiotic stress response. Here, we analyzed the evolution of miR408 in plants and verified four miR408 members in Glycine max. In the current research, highly upregulated gma-miR408 expressing was detected during nematode migration and syncytium formation response to soybean cyst nematode infection. Overexpressing and silencing miR408 vectors were transformed to soybean to confirm its potential role in plant and nematode interaction. Significant variations were observed in the MAPK signaling pathway with low OXI1, PR1, and wounding of the overexpressing lines. Overexpressing miR408 could negatively regulate soybean resistance to SCN by suppressing reactive oxygen species accumulation. Conversely, silencing miR408 positively regulates soybean resistance to SCN. Overall, gma-miR408 enhances soybean cyst nematode susceptibility by suppressing reactive oxygen species accumulation.

Identification and characterization of miRNAs and PHAS loci related to the early development of the embryo and endosperm in Fragaria × ananassa

BACKGROUND

The strawberry fleshy fruit is actually enlarged receptacle tissue, and the successful development of the embryo and endosperm is essential for receptacle fruit set. MicroRNAs (miRNAs) and phased small interfering RNAs (phasiRNAs) play indispensable regulatory roles in plant growth and development. However, miRNAs and phasiRNAs participating in the regulation of strawberry embryo and endosperm development have yet to be explored.

RESULTS

Here, we performed genome-wide identification of miRNA and phasiRNA-producing loci (PHAS) in strawberry seeds with a focus on those involved in the development of the early embryo and endosperm. We found that embryos and endosperm have different levels of small RNAs. After bioinformatics analysis, the results showed that a total of 404 miRNAs (352 known and 52 novel) and 156 PHAS genes (81 21-nt and 75 24-nt genes) could be found in strawberry seed-related tissues, of which four and nine conserved miRNA families displayed conserved expression in the endosperm and embryo, respectively. Based on refined putative annotation of PHAS loci, some auxin signal-related genes, such as CM3, TAR2, AFB2, ASA1, NAC and TAS3, were found, which demonstrates that IAA biosynthesis is important for endosperm and embryo development during early fruit growth. Additionally, some auxin signal-related conserved (miR390-TAS3) and novel (miR156-ASA1) trigger-PHAS pairs were identified.

CONCLUSIONS

Taken together, these results expand our understanding of sRNAs in strawberry embryo and endosperm development and provide a genomic resource for early-stage fruit development.

MicroRNAs: emerging regulators in horticultural crops

Horticulture is one of the oldest agricultural practices with great popularity throughout the world. Horticultural crops include fruits, vegetables, ornamental plants, as well as medicinal and beverage plants. They are cultivated for food, specific nutrition, and medical use, or for aesthetic pleasure. MicroRNAs (miRNAs), which constitute a major class of endogenous small RNAs in plants, affect a multitude of developmental and physiological processes by imparting sequence specificity to gene regulation. Over the past decade, tens of thousands of miRNAs have been identified in more than 100 horticultural crops and their critical roles in regulating quality development of diverse horticultural crops have been demonstrated. Here, we review how miRNAs have emerged as important regulators and promising tools for horticultural crop improvement.

miR164g-MsNAC022 acts as a novel module mediating drought response by transcriptional regulation of reactive oxygen species scavenging systems in apple

Under drought stress, reactive oxygen species (ROS) overaccumulate as a secondary stress that impairs plant performance and thus severely reduces crop yields. The mitigation of ROS levels under drought stress is therefore crucial for drought tolerance. MicroRNAs (miRNAs) are critical regulators of plant development and stress responses. However, the complex molecular regulatory mechanism by which they function during drought stress, especially in drought-triggered ROS scavenging, is not fully understood. Here, we report a newly identified drought-responsive miRNA, miR164g, in the wild apple species Malus sieversii and elucidate its role in apple drought tolerance. Our results showed that expression of miR164g is significantly inhibited under drought stress and it can specifically cleave transcripts of the transcription factor MsNAC022 in M. sieversii. The heterologous accumulation of miR164g in Arabidopsis thaliana results in enhanced sensitivity to drought stress, while overexpression of MsNAC022 in Arabidopsis and the cultivated apple line 'GL-3' (Malus domestica Borkh.) lead to enhanced tolerance to drought stress by raising the ROS scavenging enzymes activity and related genes expression levels, particularly PEROXIDASE (MsPOD). Furthermore, we showed that expression of MsPOD is activated by MsNAC022 in transient assays. Interestingly, Part1 (P1) region is the key region for the positive regulation of MsPOD promoter by MsNAC022, and the different POD expression patterns in M. sieversii and M. domestica is attributed to the specific fragments inserted in P1 region of M. sieversii. Our findings reveal the function of the miR164g-MsNAC022 module in mediating the drought response of M. sieversii and lay a foundation for breeding drought-tolerant apple cultivars.

Genome-wide identification of small interfering RNAs from sRNA libraries constructed from soybean cyst nematode resistant and susceptible cultivars

Plant small-RNAs regulate various biological processes by manipulating the expression of target genes at the transcriptional and post-transcriptional levels. However, little is known about the response and the functional roles of sRNAs, particularly small-interfering RNAs (siRNAs), in the soybean-soybean cyst nematode interaction. In this study, siRNA data from 24 sRNA libraries constructed from SCN-infected and non-SCN-infected resistant and susceptible soybean roots were analysed in silico. A total of 26 novel siRNAs including 17 phasiRNAs and 9 nat-siRNAs, as well as two phasiRNAs that were differentially expressed (DE) in three comparisons, were identified. Then, using qRT-PCR, the expression of majority of siRNAs was found to be downregulated after SCN infection, and the expression patterns of DE siRNAs were confirmed. Further functional annotation analyses revealed that the target genes of these siRNA were highly related to disease resistance, which included the genes coding for the NB-ARC domain, leucine-rich repeats, and Hs1pro-1 homologous proteins. Overall, the present research identified novel siRNAs and annotated their target genes, thereby laying the foundation for deciphering the roles of siRNAs in the soybean-SCN interaction.

The evolutionary history of small RNAs in Solanaceae

The Solanaceae or "nightshade" family is an economically important group with remarkable diversity. To gain a better understanding of how the unique biology of the Solanaceae relates to the family's small RNA (sRNA) genomic landscape, we downloaded over 255 publicly available sRNA data sets that comprise over 2.6 billion reads of sequence data. We applied a suite of computational tools to predict and annotate two major sRNA classes: (1) microRNAs (miRNAs), typically 20- to 22-nucleotide (nt) RNAs generated from a hairpin precursor and functioning in gene silencing and (2) short interfering RNAs (siRNAs), including 24-nt heterochromatic siRNAs typically functioning to repress repetitive regions of the genome via RNA-directed DNA methylation, as well as secondary phased siRNAs and trans-acting siRNAs generated via miRNA-directed cleavage of a polymerase II-derived RNA precursor. Our analyses described thousands of sRNA loci, including poorly understood clusters of 22-nt siRNAs that accumulate during viral infection. The birth, death, expansion, and contraction of these sRNA loci are dynamic evolutionary processes that characterize the Solanaceae family. These analyses indicate that individuals within the same genus share similar sRNA landscapes, whereas comparisons between distinct genera within the Solanaceae reveal relatively few commonalities.

Role of phasiRNAs in plant-pathogen interactions: molecular perspectives and bioinformatics tools

The genome of an organism is regulated in concert with the organized action of various genetic regulators at different hierarchical levels. Small non-coding RNAs are one of these regulators, among which microRNAs (miRNAs), a distinguished sRNA group with decisive functions in the development, growth and stress-responsive activities of both plants as well as animals, are keenly explored over a good number of years. Recent studies in plants revealed that apart from the silencing activity exhibited by miRNAs on their targets, miRNAs of specific size and structural features can direct the phasing pattern of their target loci to form phased secondary small interfering RNAs (phasiRNAs). These trigger-miRNAs were identified to target both coding and long non-coding RNAs that act as potent phasiRNA precursors or PHAS loci. The phasiRNAs produced thereby exhibit a role in enhancing further downstream regulation either on their own precursors or on those transcripts that are distinct from their genetic source of origin. Hence, these tiny regulators can stimulate an elaborative cascade of interacting RNA networks via cis and trans-regulatory mechanisms. Our review focuses on the comprehensive understanding of phasiRNAs and their trigger miRNAs, by giving much emphasis on their role in the regulation of plant defense responses, together with a summary of the computational tools available for the prediction of the same.

CsbZIP2-miR9748-CsNPF4.4 Module Mediates High Temperature Tolerance of Cucumber Through Jasmonic Acid Pathway

High temperature stress seriously affects the growth of cucumber seedlings, and even leads to a decline in yield and quality. miRNAs have been shown to be involved in regulating the response to stress in plants, but little is known about its effects on cucumber high temperature stress tolerance. Here, we found that high temperature stress induced the expression of miR9748 in cucumber. Overexpression of cucumber miR9748 in Arabidopsis improved high temperature tolerance. Transcriptome analysis revealed that miR9748 might mediate high temperature tolerance through plant hormone signal pathway. 5' RNA ligase-mediated rapid amplification of cDNA ends (5' RLM-RACE) and transient transformation technology demonstrated that CsNPF4.4 was the target gene of miR9748. CsNPF4.4 overexpression plants decreased high temperature tolerance accompanied by reducing the content of jasmonic acid (JA), but alleviated by foliar application of methyl jasmonate, indicating that CsNPF4.4 negatively regulated high temperature stress tolerance through inhibition JA signal pathway. Furthermore, high temperature stress also increased the expression level of CsbZIP2. Yeast one-hybrid and dual-luciferase assays showed that CsbZIP2 directly bound to the promoter of MIR9748 to induce its expression. Taken together, our results indicated that CsbZIP2 directly regulated miR9748 expression to cleave CsNPF4.4 to mediate high temperature tolerance through JA pathway.

MicroRNA482/2118, a miRNA superfamily essential for both disease resistance and plant development

MicroRNAs (miRNAs) are a class of 21-24 nucleotides (nt) noncoding small RNAs ubiquitously distributed across the plant kingdom. miR482/2118, one of the conserved miRNA superfamilies originating from gymnosperms, has divergent main functions in core-angiosperms. It mainly regulates NUCLEOTIDE BINDING SITE-LEUCINE-RICH REPEAT (NBS-LRR) genes in eudicots, functioning as an essential component in plant disease resistance; in contrast, it predominantly targets numerous long noncoding RNAs (lncRNAs) in monocot grasses, which are vital for plant reproduction. Usually, miR482/2118 is 22-nt in length, which can trigger the production of phased small interfering RNAs (phasiRNAs) after directed cleavage. PhasiRNAs instigated from target genes of miR482/2118 enhance their roles in corresponding biological processes by cis-regulation on cognate genes and expands their function to other pathways via trans activity on different genes. This review summarizes the origin, biogenesis, conservation, and evolutionary characteristics of the miR482/2118 superfamily and delineates its diverse functions in disease resistance, plant development, stress responses, etc.

Transcriptome analysis of transgenic apple fruit overexpressing microRNA172 reveals candidate transcription factors regulating apple fruit development at early stages

Background

MicroRNA172 (miR172) has been proven to be critical for fruit growth, since elevated miR172 activity blocks the growth of apple (Malus x domestica Borkh.) fruit. However, it is not clear how overexpression of miR172 affects apple fruit developmental processes.

Methods

To answer this question, the present study, analyzed global transcriptional changes in miR172-overexpressing (miR172OX) and nongenetically modified wild-type (WT) apple fruit at two developmental stages and in different fruit tissues via RNA-seq. In addition, two cultivars, ‘Hanfu’ and ‘M9’, which have naturally fruit size variation, were included to identify miR172-dependent DEGs. qRT–PCRwas used to verify the reliability of our RNA-seq data.

Results

Overexpression of miR172 altered the expression levels of many cell proliferation- and cell expansion-related genes. Twenty-four libraries were generated, and 10,338 differentially expressed genes (DEGs) were detected between miR172OX and WT fruit tissues. ‘Hanfu’ and ‘M9’ are two common cultivars that bear fruit of different sizes (250 g and 75 g, respectively). Six libraries were generated, and 3,627 DEGs were detected between ‘Hanfu’ and ‘M9’. After merging the two datasets, 6,888 candidate miR172-specific DEGs were identified. The potential networks associated with fruit size triggered traits were defined among genes belonging to the families of hormone synthesis, signaling pathways, and transcription factors. Our comparative transcriptome analysis provides insights into transcriptome responses to miR172 overexpression in apple fruit and a valuable database for future studies to validate functional genes and elucidate the fruit developmental mechanisms in apple.

Identification and Characterization of Csa-miR395s Reveal Their Involvements in Fruit Expansion and Abiotic Stresses in Cucumber

The miR395 plays an indispensable role in biochemical processes by regulating their target genes. However, little is known about the roles of miR395 in cucumber fruit expansion and response to abiotic stresses. Here, 4 Csa-miR395s and 8 corresponding target genes were identified in the cucumber genome. Csa-miR395s were all located on the same chromosome (Chr 5). Csa-miR395a/b/c and Csa-miR395d were distributed in different branches without a closer genetic relationship. Massive cis-acting elements, including light, phytohormone, and stress response elements, were detected in the promoter regions of Csa-MIR395s, indicating that Csa-miR395s might be involved in complex regulatory networks to control cucumber growth and development and stress response. In addition, Csa-miR395a/b/c shared the same target genes, and Csa-miR395d had its specific target genes. Tissue-specific expression analysis showed that Csa-miR395a/b/c were all expressed in the leaf, root, ovary, and expanded fruit of cucumber and highly expressed in the expanded fruits compared to the ovary, while Csa2G215520 and Csa1G502860 (target genes of Csa-miR395a/b/c) presented a downregulated trend in the expanded fruit compared to the ovary. Meanwhile, the protein co-expression network revealed that these target genes had interactions in sulfur metabolism. These results suggested that Csa-miR395a/b/c targeting Csa2G215520 and Csa1G502860 might promote cucumber fruit expansion by affecting sulfur metabolism. Additionally, Quantitative Real-time PCR analysis validated that Csa-miR395s could be regulated by NaCl stress, and Csa-miR395a/b/c could respond to PEG stress, which further confirmed the reliability of cis-acting elements data. Taken together, our results could be helpful for further exploration of the functions of miR395s in cucumber fruit expansion and response to abiotic stresses.

Biogenesis, Trafficking, and Function of Small RNAs in Plants

Small RNAs (sRNAs) encoded by plant genomes have received widespread attention because they can affect multiple biological processes. Different sRNAs that are synthesized in plant cells can move throughout the plants, transport to plant pathogens via extracellular vesicles (EVs), and transfer to mammals via food. Small RNAs function at the target sites through DNA methylation, RNA interference, and translational repression. In this article, we reviewed the systematic processes of sRNA biogenesis, trafficking, and the underlying mechanisms of its functions.

Soybean miR159-GmMYB33 Regulatory Network Involved in Gibberellin-Modulated Resistance to Heterodera glycines

Soybean cyst nematode (SCN, Heterodera glycines) is an obligate sedentary biotroph that poses major threats to soybean production globally. Recently, multiple miRNAome studies revealed that miRNAs participate in complicated soybean-SCN interactions by regulating their target genes. However, the functional roles of miRNA and target genes regulatory network are still poorly understood. In present study, we firstly investigated the expression patterns of miR159 and targeted GmMYB33 genes. The results showed miR159-3p downregulation during SCN infection; conversely, GmMYB33 genes upregulated. Furthermore, miR159 overexpressing and silencing soybean hairy roots exhibited strong resistance and susceptibility to H. glycines, respectively. In particular, miR159-GAMYB genes are reported to be involve in GA signaling and metabolism. Therefore, we then investigated the effects of GA application on the expression of miR159-GAMYB module and the development of H. glycines. We found that GA directly controls the miR159-GAMYB module, and exogenous GA application enhanced endogenous biologically active GA1 and GA3, the abundance of miR159, lowered the expression of GmMYB33 genes and delayed the development of H. glycines. Moreover, SCN infection also results in endogenous GA content decreased in soybean roots. In summary, the soybean miR159-GmMYB33 module was directly involved in the GA-modulated soybean resistance to H. glycines.

Non-Coding RNAs in Response to Drought Stress

Drought stress causes changes in the morphological, physiological, biochemical and molecular characteristics of plants. The response to drought in different plants may vary from avoidance, tolerance and escape to recovery from stress. This response is genetically programmed and regulated in a very complex yet synchronized manner. The crucial genetic regulations mediated by non-coding RNAs (ncRNAs) have emerged as game-changers in modulating the plant responses to drought and other abiotic stresses. The ncRNAs interact with their targets to form potentially subtle regulatory networks that control multiple genes to determine the overall response of plants. Many long and small drought-responsive ncRNAs have been identified and characterized in different plant varieties. The miRNA-based research is better documented, while lncRNA and transposon-derived RNAs are relatively new, and their cellular role is beginning to be understood. In this review, we have compiled the information on the categorization of non-coding RNAs based on their biogenesis and function. We also discuss the available literature on the role of long and small non-coding RNAs in mitigating drought stress in plants.

Laccase Directed Lignification Is One of the Major Processes Associated With the Defense Response Against Pythium ultimum Infection in Apple Roots

Apple replant disease (ARD), incited by a pathogen complex including Pythium ultimum, causes stunted growth or death of newly planted trees at replant sites. Development and deployment of resistant or tolerant rootstocks offers a cost-effective, ecologically friendly, and durable approach for ARD management. Maximized exploitation of natural resistance requires integrated efforts to identify key regulatory mechanisms underlying resistance traits in apple. In this study, miRNA profiling and degradome sequencing identified major miRNA pathways and candidate genes using six apple rootstock genotypes with contrasting phenotypes to P. ultimum infection. The comprehensive RNA-seq dataset offered an expansive view of post-transcriptional regulation of apple root defense activation in response to infection from P. ultimum. Several pairs of miRNA families and their corresponding targets were identified for their roles in defense response in apple roots, including miR397-laccase, miR398-superoxide dismutase, miR10986-polyphenol oxidase, miR482-resistance genes, and miR160-auxin response factor. Of these families, the genotype-specific expression patterns of miR397 indicated its fundamental role in developing defense response patterns to P. ultimum infection. Combined with other identified copper proteins, the importance of cellular fortification, such as lignification of root tissues by the action of laccase, may critically contribute to genotype-specific resistance traits. Our findings suggest that quick and enhanced lignification of apple roots may significantly impede pathogen penetration and minimize the disruption of effective defense activation in roots of resistant genotypes. The identified target miRNA species and target genes consist of a valuable resource for subsequent functional analysis of their roles during interaction between apple roots and P. ultimum.

The role of SAXS and molecular simulations in 3D structure elucidation of a DNA aptamer against lung cancer

Aptamers are short, single-stranded DNA or RNA oligonucleotide molecules that function as synthetic analogs of antibodies and bind to a target molecule with high specificity. Aptamer affinity entirely depends on its tertiary structure and charge distribution. Therefore, length and structure optimization are essential for increasing aptamer specificity and affinity. Here, we present a general optimization procedure for finding the most populated atomistic structures of DNA aptamers. Based on the existed aptamer LC-18 for lung adenocarcinoma, a new truncated LC-18 (LC-18t) aptamer LC-18t was developed. A three-dimensional (3D) shape of LC-18t was reported based on small-angle X-ray scattering (SAXS) experiments and molecular modeling by fragment molecular orbital or molecular dynamic methods. Molecular simulations revealed an ensemble of possible aptamer conformations in solution that were in close agreement with measured SAXS data. The aptamer LC-18t had stronger binding to cancerous cells in lung tumor tissues and shared the binding site with the original larger aptamer. The suggested approach reveals 3D shapes of aptamers and helps in designing better affinity probes.

MicroRNAs in Woody Plants

MicroRNAs (miRNAs) are small (∼21-nucleotides) non-coding RNAs found in plant and animals. MiRNAs function as critical post-transcriptional regulators of gene expression by binding to complementary sequences in their target mRNAs, leading to mRNA destabilization and translational inhibition. Plant miRNAs have some distinct characteristics compared to their animal counterparts, including greater evolutionary conservation and unique miRNA processing methods. The lifecycle of a plant begins with embryogenesis and progresses through seed germination, vegetative growth, reproductive growth, flowering and fruiting, and finally senescence and death. MiRNAs participate in the transformation of plant growth and development and directly monitor progression of these processes and the expression of certain morphological characteristics by regulating transcription factor genes involved in cell growth and differentiation. In woody plants, a large and rapidly increasing number of miRNAs have been identified, but their biological functions are largely unknown. In this review, we summarize the progress of miRNA research in woody plants to date. In particular, we discuss the potential roles of these miRNAs in growth, development, and biotic and abiotic stresses responses in woody plants.

Identification and Expression Analysis of miR160 and Their Target Genes in Cucumber

miR160 plays a crucial role in various biological processes by regulating their target gene auxin response factor (ARF) in plants. However, little is known about miR160 and ARF in cucumber fruit expansion. Here, 4 Csa-MIR160 family members and 17 CsARFs were identified through a genome-wide search. Csa-miR160 showed a closer relationship with those in melon. Phylogenetic analysis revealed that CsARFs were divided into four classes and most of CsARFs presented a closer evolutionary relationship with those from tomato. Putative cis-elements analysis predicted that Csa-MIR160 and CsARFs were involved in light, phytohormone and stress response, which proved that they might take part in light, phytohormone and stress signaling pathway by the miR160-ARF module. In addition, CsARF5, CsARF11, CsARF13 and CsARF14 were predicted as the target genes of Csa-miR160. qRT-PCR revealed that Csa-miR160 and their target gene CsARFs were differentially expressed in differential cucumber tissues and developmental stages. Csa-miR160d was only expressed in the expanded cucumber fruit. CsARF5, CsARF11 and CsARF13 exhibited the lower expression in the expanded fruit than those in the ovary, while, CsARF14 showed the reverse trend. Our results suggested that Csa-miR160d might play a crucial role in cucumber fruit expansion by negatively targeting CsARF5, CsARF11 and CsARF13. This is the first genome-wide analysis of miR160 in cucumber. These findings provide useful information and resources for further studying the role of miR160 and ARF in cucumber fruit expansion.

Identification of miRNAs and Their Targets Involved in Flower and Fruit Development across Domesticated and Wild Capsicum Species

MicroRNAs (miRNAs) are regulators of the post-transcription stage of gene activity documented to play central roles in flower and fruit development in model plant species. However, little is known about their roles and differences in domesticated and wild Capsicum species. In this study, we used high-throughput sequencing to analyze the miRNA content at three developmental stages (flower, small fruit, and middle fruit) from two cultivated (C. baccatum and C. annuum) and two wild (C. chacoense and C. eximium) pepper species. This analysis revealed 22 known and 27 novel miRNAs differentially expressed across species and tissues. A number of stage- and species-specific miRNAs were identified, and Gene Ontology terms were assigned to 138 genes targeted by the miRNAs. Most Gene Ontology terms were for the categories "genetic information processing", "signaling and cellular processes", "amino acid metabolism", and "carbohydrate metabolism". Enriched KEGG analysis revealed the pathways amino acids, sugar and nucleotide metabolism, starch and sucrose metabolism, and fructose-mannose metabolism among the principal ones regulated by miRNAs during pepper fruit ripening. We predicted miRNA-target gene interactions regulating flowering time and fruit development, including miR156/157 with SPL genes, miR159 with GaMYB proteins, miR160 with ARF genes, miR172 with AP2-like transcription factors, and miR408 with CLAVATA1 gene across the different Capsicum species. In addition, novel miRNAs play an important role in regulating interactions potentially controlling plant pathogen defense and fruit quality via fructokinase, alpha-L-arabinofuranosidase, and aromatic and neutral amino acid transporter. Overall, the small RNA-sequencing results from this study represent valuable information that provides a solid foundation for uncovering the miRNA-mediated mechanisms of flower and fruit development between domesticated and wild Capsicum species.

sRNAanno-a database repository of uniformly annotated small RNAs in plants

Small RNAs (sRNAs) are essential regulatory molecules, and there are three major sRNA classes in plants: microRNAs (miRNAs), phased small interfering RNAs (phased siRNAs or phasiRNAs), and heterochromatic siRNAs (hc-siRNAs). Excluding miRNAs, the other two classes are not well annotated or available in public databases for most sequenced plant genomes. We performed a comprehensive sRNA annotation of 143 plant species that have fully sequenced genomes and next-generation sequencing sRNA data publicly available. The results are available via an online repository called sRNAanno ( www.plantsRNAs.org ). Compared with other public plant sRNA databases, we obtained was much more miRNA annotations, which are more complete and reliable because of the consistent and highly stringent criteria used in our miRNA annotations. sRNAanno also provides free access to genomic information for >22,721 PHAS loci and >22 million hc-siRNA loci annotated from these 143 plant species. Both miRNA and PHAS loci can be easily browsed to view their main features, and a collection of archetypal trans-acting siRNA 3 (TAS3) genes were annotated separately for quick access. To facilitate the ease of sRNA annotation, sRNAanno provides free service for sRNA annotations to the community. In summary, the sRNAanno database is a great resource to facilitate genomic and genetic research on plant small RNAs.

Relationship between Secondary Metabolism and miRNA for Important Flavor Compounds in Different Tissues of Tea Plant (Camellia sinensis) As Revealed by Genome-Wide miRNA Analysis

This study investigated the regulatory relationship between important flavor compounds and microRNA (miRNA) in nine different tissues of tea plant by analyzing the related metabolites, small RNAs (sRNAs), degradome, and coexpression network. A total of 272 differential expressed miRNAs (DEmiRNAs) were obtained, including 198 conserved miRNAs and 74 novel miRNAs. Meanwhile, the expression patterns of miR159-GAMYB, miR167-ARF, and miR396-GRF pairs were investigated by quantitative real-time polymerase chain reaction (qRT-PCR) and the target sites were verified by 5'RNA ligase-mediated RACE (5' RLM-RACE). Further coexpression analysis showed that the content of gallated catechins was significantly and negatively correlated with the expression of miR156, but positively correlated with the expression of miR166 and miR172. Additionally, the expression of miR169a, miR169l, and miR319h was shown to be positively correlated with the content of nongallated catechins and the experssion levels of ANRa, ANRb, and LARb. Moreover, important volatile compounds, such as linalool, geraniol, and 2-phenylethanol, were found to be highly positively correlated with the expression of miR171o, miRN71a, miRN71b, miRN71c, and miRN71d. Our data indicate that these miRNAs may play important roles in regulating the biosynthesis of flavor compounds in different tissues of tea plant.

Integrated analyses of the transcriptome and small RNA of the hemiparasitic plant Monochasma savatieri before and after establishment of parasite-host association

BACKGROUND

Monochasma savatieri is a medicinal root hemiparasitic herb that extracts water and nutrients from the host plant via a haustorium. M. savatieri exhibits an enhanced growth after the establishment of parasite-host associations, but little is known about the molecular mechanism responsible. In this study, endogenous hormones, RNA sequencing and small RNA sequencing analysis were performed on M. savatieri before and after establishment of parasite-host associations.

RESULTS

When grown with the host, decreased contents of jasmonic acid (JA) and indole-3-acetic acid (IAA) and increased abscisic acid (ABA) content were observed in M. savatieri with the established parasitic relationship. When grown with the host, 46,424 differentially expressed genes (DEGs) and 162 differentially expressed miRNAs (DEmiRs) were identified in the comparison between M. savatieri with the established parasitic relationship and without the established parasitic relationship. Analysis of Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) showed that these DEGs and targets of DEmiRs mostly participated in plant hormone signal transduction, starch and sucrose metabolism, carbohydrate metabolism, cell growth and death, and transport and catabolism. Furthermore, correlation analysis of mRNA and miRNA revealed that 10 miRNA-target pairs from novel_mir65, novel_mir40, novel_mir80, miR397-5p_1, novel_mir36, novel_mir25 and novel_mir17 may have important roles in regulating the parasitic development of M. savatieri.

CONCLUSIONS

Our study not only expands the understanding of enhanced growth in M. savatieri after the establishment of parasite-host associations, but also first provides abundant resources for future molecular and genetic studies in M. savatieri.

Illumina sequencing revealed roles of microRNAs in different aluminum tolerance of two citrus species

Self-germinated seedlings of Citrus sinensis and C. grandis were supplied with nutrient solution with 0 mM AlCl₃·6H₂O (control, -Al) or 1 mM AlCl₃·6H₂O (+Al) for 18 weeks. The DW (Dry weights) of leaf, stem, shoot and the whole plant of C. grandis were decreased and the ratio of root DW to shoot DW in C. grandis were increased by Al, whereas these parameters of C. sinensis were not changed by Al. Al treatment dramatically decreased the sulfur (S) content in C. grandis roots and the phosphorus (P) content in both C. sinensis and C. grandis roots. More Al was transported to shoots and leaves in C. grandis than in C. sinensis under Al treatment. Al treatment has more adverse effects on C. grandis than on C. sinensis, as revealed by the higher production of superoxide anion (O₂ ^·-), H₂O₂ and thiobarbituric acid reactive substace (TBARS) content in C. grandis roots. Via the Illumina sequencing technique, we successfully identified and quantified 12 and 16 differentially expressed miRNAs responding to Al stress in C. sinensis and C. grandis roots, respectively. The possible mechanism underlying different Al tolerance of C. sinensis and C. grandis were summarized as having following aspects: (a) enhancement of adventitious and lateral root development (miR160); (b) up-regulation of stress and signaling transduction related genes, such as SGT1, PLC and AAO (miR477, miR397 and miR398); (c) enhancement of citrate secretion (miR3627); (d) more flexible control of alternative glycolysis pathway and TCA cycle (miR3627 and miR482); (e) up-regulation of S-metabolism (miR172); (f) more flexible control of miRNA metabolism. For the first time, we showed that root development (miR160) and cell wall components (cas-miR5139, csi-miR12105) may play crucial roles in Al tolerance in citrus plants. In conclusion, our study provided a comprehensive profile of differentially expressed miRNAs in response to Al stress between two citrus plants differing in Al tolerance which further enriched our understanding of the molecular mechanism underlying Al tolerance in plants.

The Apple microR171i-SCARECROW-LIKE PROTEINS26.1 Module Enhances Drought Stress Tolerance by Integrating Ascorbic Acid Metabolism

Drought stress severely restricts crop yield and quality. Small noncoding RNAs play critical roles in plant growth, development, and stress responses by regulating target gene expression, but their roles in drought stress tolerance in apple (Malus domestica) are poorly understood. Here, we identified various small noncoding RNAs and their targets from the wild apple species Malus sieversii via high-throughput sequencing and degradome analysis. Forty known microRNAs (miRNAs) and eight new small noncoding RNAs were differentially expressed in response to 2 or 4 h of drought stress treatment. We experimentally verified the expression patterns of five selected miRNAs and their targets. We established that one miRNA, mdm-miR171i, specifically targeted and degraded SCARECROW-LIKE PROTEINS26 1 (MsSCL26 1) transcripts. Both knockout of mdm-miR171i and overexpression of MsSCL26 1 improved drought stress tolerance in the cultivated apple line 'GL-3' by regulating the expression of antioxidant enzyme genes, especially that of MONODEHYDROASCORBATE REDUCTASE, which functions in metabolism under drought stress. Transient expression analysis demonstrated that MsSCL26.1 activates MsMDHAR transcription by positively regulating the activity of the P1 region in its promoter. Therefore, the miR171i-SCL26 1 module enhances drought stress tolerance in apple by regulating antioxidant gene expression and ascorbic acid metabolism.

Drought-Induced Regulatory Cascades and Their Effects on the Nutritional Quality of Developing Potato Tubers

Competition for scarce water resources and the continued effects of global warming exacerbate current constraints on potato crop production. While plants’ response to drought in above-ground tissues has been well documented, the regulatory cascades and subsequent nutritive changes in developing tubers have been largely unexplored. Using the commercial Canadian cultivar “Vigor”, plants were subjected to a gradual drought treatment under high tunnels causing a 4 °C increase in the canopy temperature. Tubers were sampled for RNAseq and metabolite analysis. Approximately 2600 genes and 3898 transcripts were differentially expressed by at least 4-fold in drought-stressed potato tubers, with 75% and 69% being down-regulated, respectively. A further 229 small RNAs were implicated in gene regulation during drought. Expression of several small RNA clusters negatively correlated with expression of their six target patatin genes, suggesting involvement in the regulation of storage proteins during drought. The comparison of protein homologues between Solanum tuberosum L. and Arabidopsis thaliana L. indicated that down-regulated genes were associated with phenylpropanoid and carotenoid biosynthesis. As is indicative of reduced flow through the phenylpropanoid pathway, phenylalanine accumulated in drought-stressed tubers. This suggests that there may be nutritive implications to drought stress occurring during the potato tuber bulking phase in sensitive cultivars.

Skin Color in Apple Fruit (Malus × domestica): Genetic and Epigenetic Insights

Apple skin color is an important trait for organoleptic quality. In fact, it has a major influence on consumer choice. Skin color is, thus, one of the most important criteria taken into account by breeders. For apples, most novel varieties are so-called "mutants" or "sports" that have been identified in clonal populations. Indeed, many "sports" exist that show distinct phenotypic differences compared to the varieties from which they originated. These differences affect a limited number of traits of economic importance, including skin color. Until recently, the detailed genetic or epigenetic changes resulting in heritable phenotypic changes in sports was largely unknown. Recent technological advances and the availability of several high-quality apple genomes now provide the bases to understand the exact nature of the underlying molecular changes that are responsible for the observed phenotypic changes observed in sports. The present review investigates the molecular nature of sports affected in apple skin color giving arguments in favor of the genetic or epigenetic explanatory models.

Apple SERRATE negatively mediates drought resistance by regulating MdMYB88 and MdMYB124 and microRNA biogenesis

The function of serrate (SE) in miRNA biogenesis in Arabidopsis is well elucidated, whereas its role in plant drought resistance is largely unknown. In this study, we report that MdSE acts as a negative regulator of apple (Malus × domestica) drought resistance by regulating the expression levels of MdMYB88 and MdMYB124 and miRNAs, including mdm-miR156, mdm-miR166, mdm-miR172, mdm-miR319, and mdm-miR399. MdSE interacts with MdMYB88 and MdMYB124, two positive regulators of apple drought resistance. MdSE decreases the transcript and protein levels of MdMYB88 and MdMYB124, which directly regulate the expression of MdNCED3, a key enzyme in abscisic acid (ABA) biosynthesis. Furthermore, MdSE is enriched in the same region of the MdNECD3 promoter where MdMYB88/MdMYB124 binds. Consistently, MdSE RNAi transgenic plants are more sensitive to ABA-induced stomatal closure, whereas MdSE OE plants are less sensitive. In addition, under drought stress, MdSE is responsible for the biogenesis of mdm-miR399, a negative regulator of drought resistance, and negatively regulates miRNAs, including mdm-miR156, mdm-miR166, mdm-miR172, and mdm-miR319, which are positive regulators of drought resistance. Taken together, by revealing the negative role of MdSE, our results broaden our understanding of the apple drought response and provide a candidate gene for apple drought improvement through molecular breeding.

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.

Integrated annotations and analyses of small RNA-producing loci from 47 diverse plants

Plant endogenous small RNAs (sRNAs) are important regulators of gene expression. There are two broad categories of plant sRNAs: microRNAs (miRNAs) and endogenous short interfering RNAs (siRNAs). MicroRNA loci are relatively well-annotated but compose only a small minority of the total sRNA pool; siRNA locus annotations have lagged far behind. Here, we used a large data set of published and newly generated sRNA sequencing data (1333 sRNA-seq libraries containing more than 20 billion reads) and a uniform bioinformatic pipeline to produce comprehensive sRNA locus annotations of 47 diverse plants, yielding more than 2.7 million sRNA loci. The two most numerous classes of siRNA loci produced mainly 24- and 21-nucleotide (nt) siRNAs, respectively. Most often, 24-nt-dominated siRNA loci occurred in intergenic regions, especially at the 5′-flanking regions of protein-coding genes. In contrast, 21-nt-dominated siRNA loci were most often derived from double-stranded RNA precursors copied from spliced mRNAs. Genic 21-nt-dominated loci were especially common from disease resistance genes, including from a large number of monocots. Individual siRNA sequences of all types showed very little conservation across species, whereas mature miRNAs were more likely to be conserved. We developed a web server where our data and several search and analysis tools are freely accessible.

MicroRNA528, a hub regulator modulating ROS homeostasis via targeting of a diverse set of genes encoding copper-containing proteins in monocots

Plant microRNAs (miRNAs) regulate vital cellular processes, including responses to extreme temperatures with which reactive oxygen species (ROS) are often closely associated. In the present study, it was found that aberrant temperatures caused extensive changes in abundance to numerous miRNAs in banana fruit, especially the copper (Cu)-associated miRNAs. Among them, miR528 was significantly downregulated under cold stress and it was found to target genes encoding polyphenol oxidase (PPO), different from those identified in rice and maize. Expression of PPO genes was upregulated by > 100-fold in cold conditions, leading to ROS surge and subsequent peel browning of banana fruit. Extensive comparative genomic analyses revealed that the monocot-specific miR528 can potentially target a large collection of genes encoding Cu-containing proteins. Most of them are actively involved in cellular ROS metabolism, including not only ROS generating oxidases, but also ROS scavenging enzymes. It also was demonstrated that miR528 has evolved a distinct preference of target genes in different monocots, with its target site varying in position among/within gene families, implying a highly dynamic process of target gene diversification. Its broad capacity to target genes encoding Cu-containing protein implicates miR528 as a key regulator for modulating the cellular ROS homeostasis in monocots.