Identification and functional analysis of novel and conserved microRNAs in tomato

Identification of Tomato microRNAs in Late Response to Trichoderma atroviride

The tomato (Solanum lycopersicum) is an important crop worldwide and is considered a model plant to study stress responses. Small RNAs (sRNAs), 21-24 nucleotides in length, are recognized as a conserved mechanism for regulating gene expression in eukaryotes. Plant endogenous sRNAs, such as microRNA (miRNA), have been involved in disease resistance. High-throughput RNA sequencing was used to analyze the miRNA profile of the aerial part of 30-day-old tomato plants after the application of the fungus Trichoderma atroviride to the seeds at the transcriptional memory state. Compared to control plants, ten differentially expressed (DE) miRNAs were identified in those inoculated with Trichoderma, five upregulated and five downregulated, of which seven were known (miR166a, miR398-3p, miR408, miR5300, miR6024, miR6027-5p, and miR9471b-3p), and three were putatively novel (novel miR257, novel miR275, and novel miR1767). miRNA expression levels were assessed using real-time quantitative PCR analysis. A plant sRNA target analysis of the DE miRNAs predicted 945 potential target genes, most of them being downregulated (84%). The analysis of KEGG metabolic pathways showed that most of the targets harbored functions associated with plant-pathogen interaction, membrane trafficking, and protein kinases. Expression changes of tomato miRNAs caused by Trichoderma are linked to plant defense responses and appear to have long-lasting effects.

The applications of CRISPR/Cas-mediated microRNA and lncRNA editing in plant biology: shaping the future of plant non-coding RNA research

MAIN CONCLUSION

CRISPR/Cas technology has greatly facilitated plant non-coding RNA (ncRNA) biology research, establishing itself as a promising tool for ncRNA functional characterization and ncRNA-mediated plant improvement. Throughout the last decade, the promising genome editing tool clustered regularly interspaced short palindromic repeat (CRISPR)/CRISPR-associated proteins (Cas; CRISPR/Cas) has allowed unprecedented advances in the field of plant functional genomics and crop improvement. Even though CRISPR/Cas-mediated genome editing system has been widely used to elucidate the biological significance of a number of plant protein-coding genes, this technology has been barely applied in the functional analysis of those non-coding RNAs (ncRNAs) that modulate gene expression, such as microRNAs (miRNAs) and long non-coding RNAs (lncRNAs). Nevertheless, compelling findings indicate that CRISPR/Cas-based ncRNA editing has remarkable potential for deciphering the biological roles of ncRNAs in plants, as well as for plant breeding. For instance, it has been demonstrated that CRISPR/Cas tool could overcome the challenges associated with other approaches employed in functional genomic studies (e.g., incomplete knockdown and off-target activity). Thus, in this review article, we discuss the current status and progress of CRISPR/Cas-mediated ncRNA editing in plant science in order to provide novel prospects for further assessment and validation of the biological activities of plant ncRNAs and to enhance the development of ncRNA-centered protocols for crop improvement.

Sly-miR398 Participates in Cadmium Stress Acclimation by Regulating Antioxidant System and Cadmium Transport in Tomato (Solanum lycopersicum)

Cadmium (Cd) pollution is one of the major threats in agricultural production, and can cause oxidative damage and growth limitation in plants. MicroRNA398 (miR398) is involved in plant resistance to different stresses, and the post-transcriptional regulation of miR398 on CSDs plays a key role. Here, we report that miR398 was down-regulated in tomato in response to Cd stress. Simultaneously, CSD1 and SOD were up-regulated, with CSD2 unchanged, suggesting CSD1 is involved in miR398-induced regulation under Cd stress. In addition, the role of miR398 in Cd tolerance in tomato was evaluated using a transgenic line overexpressing MIR398 (miR398#OE) in which the down-expression of miR398 was disrupted. The results showed that Cd stress induced more significant growth inhibition, oxidative damage, and antioxidant enzymes disorder in miR398#OE than that in wild type (WT). Moreover, higher Cd concentration in the shoot and xylem sap, and net Cd influx rate, were observed in miR398#OE, which could be due to the increased Cd uptake genes (IRT1, IRT2, and NRAMP2) and decreased Cd compartmentalization gene HMA3. Overall, our results indicate that down-regulated miR398 plays a protective role in tomato against Cd stress by modulating the activity of antioxidant enzymes and Cd uptake and translocation.

miRNAs and lncRNAs in tomato: Roles in biotic and abiotic stress responses

Plants are continuously exposed to various biotic and abiotic stresses in the natural environment. To cope with these stresses, they have evolved a multitude of defenses mechanisms. With the rapid development of genome sequencing technologies, a large number of non-coding RNA (ncRNAs) have been identified in tomato, like microRNAs (miRNAs) and long non-coding RNAs (lncRNAs). Recently, more and more evidence indicates that many ncRNAs are involved in plant response to biotic and abiotic stresses in tomato. In this review, we summarize recent updates on the regulatory roles of ncRNAs in tomato abiotic/biotic responses, including abiotic (high temperature, drought, cold, salinization, etc.) and biotic (bacteria, fungi, viruses, insects, etc.) stresses. Understanding the molecular mechanisms mediated by ncRNAs in response to these stresses will help us to clarify the future directions for ncRNA research and resistance breeding in tomato.

Stu-miR827-Targeted StWRKY48 Transcription Factor Negatively Regulates Drought Tolerance of Potato by Increasing Leaf Stomatal Density

Stomata are specialized portals in plant leaves to modulate water loss from plants to the atmosphere by control of the transpiration, thereby determining the water-use efficiency and drought resistance of plants. Despite that the stomata developmental progression is well-understood at the molecular level, the experimental evidence that miRNA regulates stomata development is still lacking, and the underlying mechanism remains elusive. This study demonstrates the involvement of stu-miR827 in regulating the drought tolerance of potato due to its control over the leaf stomatal density. The expression analysis showed that stu-miR827 was obviously repressed by drought stresses and then rapidly increased after rewatering. Suppressing the expression of stu-miR827 transgenic potato lines showed an increase in stomatal density, correlating with a weaker drought resistance compared with wildtype potato lines. In addition, StWRKY48 was identified as the target gene of stu-miR827, and the expression of StWRKY48 was obviously induced by drought stresses and was greatly upregulated in stu-miR827 knockdown transgenic potato lines, suggesting its involvement in the drought stress response. Importantly, the expression of genes associated with stomata development, such as SDD (stomatal density and distribution) and TMM (too many mouths), was seriously suppressed in transgenic lines. Altogether, these observations demonstrated that suppression of stu-miR827 might lead to overexpression of StWRKY48, which may contribute to negatively regulating the drought adaptation of potato by increasing the stomatal density. The results may facilitate functional studies of miRNAs in the process of drought tolerance in plants.

MicroRNA398: A Master Regulator of Plant Development and Stress Responses

MicroRNAs (miRNAs) play crucial roles in plant development and stress responses, and a growing number of studies suggest that miRNAs are promising targets for crop improvement because they participate in the regulation of diverse, important agronomic traits. MicroRNA398 (miR398) is a conserved miRNA in plants and has been shown to control multiple stress responses and plant growth in a variety of species. There are many studies on the stress response and developmental regulation of miR398. To systematically understand its function, it is necessary to summarize the evolution and functional roles of miR398 and its target genes. In this review, we analyze the evolution of miR398 in plants and outline its involvement in abiotic and biotic stress responses, in growth and development and in model and non-model plants. We summarize recent functional analyses, highlighting the role of miR398 as a master regulator that coordinates growth and diverse responses to environmental factors. We also discuss the potential for fine-tuning miR398 to achieve the goal of simultaneously improving plant growth and stress tolerance.

Expression Patterns of miR398, miR167, and miR159 in the Interaction between Bread Wheat (Triticum aestivum L.) and Pathogenic Fusarium culmorum and Beneficial Trichoderma Fungi

Bread wheat (Triticum aestivum L.) is an agronomically significant cereal cultivated worldwide. Wheat breeding is limited by numerous abiotic and biotic stresses. One of the most deleterious factors is biotic stress provoked by the Fusarium culmorum fungus. This pathogen is a causative agent of Fusarium root rot and Fusarium head blight. Beneficial fungi Trichoderma atroviride and T. cremeum are strong antagonists of mycotoxigenic Fusarium spp. These fungi promote plant growth and enhance their tolerance of negative environmental conditions. The aim of the study was to determine and compare the spatial (in above- and underground organs) and temporal (early: 6 and 22 hpi; and late: 5 and 7 dpi reactions) expression profiles of three mature miRNAs (miR398, miR167, and miR159) in wheat plants inoculated with two strains of F. culmorum (KF846 and EW49). Moreover, the spatial expression patterns in wheat response between plants inoculated with beneficial T. atroviride (AN35) and T. cremeum (AN392) were assessed. Understanding the sophisticated role of miRNAs in wheat-fungal interactions may initiate a discussion concerning the use of this knowledge to protect wheat plants from the harmful effects of fungal pathogens. With the use of droplet digital PCR (ddPCR), the absolute quantification of the selected miRNAs in the tested material was carried out. The differential accumulation of miR398, miR167, and miR159 in the studied groups was observed. The abundance of all analyzed miRNAs in the roots demonstrated an increase in the early and reduction in late wheat response to F. culmorum inoculation, suggesting the role of these particles in the initial wheat reaction to the studied fungal pathogen. The diverse expression patterns of the studied miRNAs between Trichoderma-inoculated or F. culmorum-inoculated plants and control wheat, as well as between Trichoderma-inoculated and F. culmorum-inoculated plants, were noticed, indicating the need for further analysis.

Degradome Analysis of Tomato and Nicotiana benthamiana Plants Infected with Potato Spindle Tuber Viroid

Viroids are infectious non-coding RNAs that infect plants. During infection, viroid RNAs are targeted by Dicer-like proteins, generating viroid-derived small RNAs (vd-sRNAs) that can guide the sequence specific cleavage of cognate host mRNAs via an RNA silencing mechanism. To assess the involvement of these pathways in pathogenesis associated with nuclear-replicating viroids, high-throughput sequencing of sRNAs and degradome analysis were carried out on tomato and Nicotiana benthamiana plants infected by potato spindle tuber viroid (PSTVd). Both hosts develop similar stunting and leaf curling symptoms when infected by PSTVd, thus allowing comparative analyses. About one hundred tomato mRNAs potentially targeted for degradation by vd-sRNAs were initially identified. However, data from biological replicates and comparisons between mock and infected samples reduced the number of bona fide targets-i.e., those identified with high confidence in two infected biological replicates but not in the mock controls-to only eight mRNAs that encode proteins involved in development, transcription or defense. Somewhat surprisingly, results of RT-qPCR assays revealed that the accumulation of only four of these mRNAs was inhibited in the PSTVd-infected tomato. When these analyses were extended to mock inoculated and PSTVd-infected N. benthamiana plants, a completely different set of potential mRNA targets was identified. The failure to identify homologous mRNA(s) targeted by PSTVd-sRNA suggests that different pathways could be involved in the elicitation of similar symptoms in these two species. Moreover, no significant modifications in the accumulation of miRNAs and in the cleavage of their targeted mRNAs were detected in the infected tomato plants with respect to the mock controls. Taken together, these data suggest that stunting and leaf curling symptoms induced by PSTVd are elicited by a complex plant response involving multiple mechanisms, with RNA silencing being only one of the possible components.

Deciphering the role of microRNAs during Pi54 gene mediated Magnaporthe oryzae resistance response in rice

The broad-spectrum resistance gene Pi54 confers resistance to multiple isolates of Magnaporthe oryzae in rice. In order to decipher the molecular mechanism underlying the Pi54 mediated resistance in rice line Taipei309 ^Pi54 (carrying Pi54), miRNAome study was performed at 24 h post-inoculation (hpi) with M. oryzae. A total of 222 known miRNAs representing 101 miRNA families were found in this study. Of these, 29 and 24 miRNAs were respectively up- and down-regulated in the resistant Taipei309 ^Pi54 . Defence response (DR) genes, like, NBSGO35, and OsWAK129b, and genes related to transcription factors were up-regulated in Taipei309 ^Pi54 line. The vast array of miRNA candidates identified here are miR159c, miR167c, miR2100, miR2118o, miR2118l, miR319a, miR393, miR395l, miR397a, miR397b, miR398, miR439g, miR531b, miR812f, and miR815c, and they manifest their role in balancing the interplay between various DR genes during Pi54 mediated resistance. We also validated miRNA/target gene pairs involved in hormone signalling, and cross-talk among hormone pathways regulating the rice immunity. This study suggests that the Pi54 gene mediated blast resistance is influenced by several microRNAs through PTI and ETI components in the rice line Taipei309 ^Pi54 , leading to incompatible host-pathogen interaction.

Unique miRNAs and their targets in tomato leaf responding to combined drought and heat stress

BACKGROUND

Both drought and heat stress are serious global problems, leading to agricultural production loss. MicroRNAs (miRNAs) play important roles in plant species responding to individual drought and heat stress. However, the miRNAs and mRNAs in association with combined drought and heat in crops like tomato remains unclear.

RESULTS

We studied the crosstalk of miRNAs and their target genes in tomato plants grown under simultaneous drought and heat stress that frequently happen in field conditions. In total, 335 known miRNAs representing 55 miRNA families and 430 potential novel miRNAs were identified in Solanum lycopersicum L. using small RNA deep sequencing. Through expression analysis, miRNAs in association with drought, heat and the combination of these were investigated. In total, 61, 74 and 37 miRNAs were differentially regulated for combination (of both stresses) vs control, combination vs drought and combination vs heat, respectively. Target genes with different expression levels were found using degradome sequencing, which were mainly involved in transcription factor activity, sequence-specific DNA binding, transcription, regulation of transcription, nucleus, DNA binding etc. The quantitative real-time polymerase chain reaction (qRT-PCR) results confirmed the accuracy of sequencing.

CONCLUSIONS

Our study serves as valuable knowledge on how crop adapted to combined drought and heat stress by regulating miRNAs and mRNAs, which provide information for crop improvement to deal with future climate changes.

Identification and characterisation of microRNAs and their target genes in phosphate-starved Nicotiana benthamiana by small RNA deep sequencing and 5'RACE analysis

BACKGROUND

Phosphorus is an important macronutrient that is severely lacking in soils. In plants, specific microRNAs (miRNAs) essential for nutrient management and the regulation of stress responses are responsible for the control of many phosphate starvation responses. Further understanding of conserved and species-specific microRNA species has potential implications for the development of crops tolerant to soils with low phosphate.

RESULTS

This study identified and characterised phosphate starvation-responsive miRNAs in the native Australian tobacco Nicotiana benthamiana. Small RNA libraries were constructed and sequenced from phosphate-starved plant leaves, stems and roots. Twenty-four conserved miRNA families and 36 species-specific miRNAs were identified. The majority of highly phosphate starvation-responsive miRNAs were highly conserved, comprising of members from the miR399, miR827, and miR2111 families. In addition, two miRNA-star species were identified to be phosphate starvation-responsive. A total of seven miRNA targets were confirmed using RLM-5'RACE to be cleaved by five miRNA families, including two confirmed cleavage targets for Nbe-miR399 species, one for Nbe-miR2111, and two for Nbe-miR398. A number of N. benthamiana-specific features for conserved miRNAs were identified, including species-specific miRNA targets predicted or confirmed for miR399, miR827, and miR398.

CONCLUSIONS

Our results give an insight into the phosphate starvation-responsive miRNAs of Nicotiana benthamiana, and indicate that the phosphate starvation response pathways in N. benthamiana contain both highly conserved and species-specific components.

New insights into tomato microRNAs

Cultivated tomato, Solanum lycopersicum, is one of the most common fruits in the global food industry. Together with the wild tomato Solanum pennellii, it is widely used for developing better cultivars. MicroRNAs affect mRNA regulation, inhibiting its translation and/or promoting its degradation. Important proteins involved in these processes are ARGONAUTE and DICER. This study aimed to identify and characterize the genes involved in the miRNA processing pathway, miRNA molecules and target genes in both species. We validated the presence of pathway genes and miRNA in different NGS libraries and 6 miRNA families using quantitative RT-PCR. We identified 71 putative proteins in S. lycopersicum and 108 in S. pennellii likely involved in small RNAs processing. Of these, 29 and 32 participate in miRNA processing pathways, respectively. We identified 343 mature miRNAs, 226 pre-miRNAs in 87 families, including 192 miRNAs, which were not previously identified, belonging to 38 new families in S. lycopersicum. In S. pennellii, we found 388 mature miRNAs and 234 pre-miRNAs contained in 85 families. All miRNAs found in S. pennellii were unpublished, being identified for the first time in our study. Furthermore, we identified 2471 and 3462 different miRNA target in S. lycopersicum and S. pennellii, respectively.

Function identification of miR482b, a negative regulator during tomato resistance to Phytophthora infestans

Tomato is an important horticultural and economic crop cultivated worldwide. As Phytophthora infestans becomes a huge threat to tomato production, it is necessary to study the resistance mechanisms of tomato against P. infestans. Our previous research has found that miR482 might be involved in tomato-P. infestans interaction. In this study, miR482b precursor was cloned from Solanum pimpinellifolium "L3708" and miR482b was shown to decrease in abundance in tomato following P. infestans infection. Compared to wild-type tomato plants, tomato plants that overexpressed miR482b displayed more serious disease symptoms after P. infestans infection, with more necrotic cells, longer lesion diameters, and increased P. infestans abundance. Meanwhile, silencing of miR482b was performed by short tandem target mimic (STTM), resulting in enhancement of tomato resistance to P. infestans. Using miRNA and degradome data sets, NBS-LRR disease-resistance genes targeted by miR482b were validated. Negative correlation between the expression of miR482b and its target genes was found in all miR482b-overexpressing and -silencing tomato plants. Our results provide insight into tomato miR482b involved in the response to P. infestans infection, and demonstrate that miR482b-NBS-LRR is an important component in the network of tomato-P. infestans interaction.

Effective enhancement of resistance to Phytophthora infestans by overexpression of miR172a and b in Solanum lycopersicum

Overexpression of miR172a and b in tomato ( Solanum lycopersicum ) Zaofen No. 2 increased resistance to Phytophthora infestans infection by suppressing of an AP2/ERF transcription factor. The miR172 family has been shown to participate in the growth phase transition, flowering time control, abiotic and biotic stresses by regulating the expression of a small group of AP2/ERF transcription factors. In this study, the precursors of miR172a and b were cloned from tomato, Solanum pimpinellifolium L3708. We used the degradome sequencing to determine the cleavage site of miR172 to a member of the AP2/ERF transcription factor family (Solyc11g072600.1.1). qRT-PCR results showed that the expression of AP2/ERF was negatively correlated with the expression of miR172 in S. pimpinellifolium L3708 infected with Phytophthora infestans. Overexpression of miR172a and b in S. lycopersicum Zaofen No. 2 conferred greater resistance to P. infestans infection, as evidenced by decreased disease index, lesion sizes, and P. infestans abundance. The SOD and POD play important roles in scavenging late massive ROS in plant-pathogen interaction. Malonaldehyde (MDA) is widely recognized as an indicator of lipid peroxidation. Membrane damage in plants can be estimated by measuring leakage of electrolytes, which is evaluated by determining relative electrolyte leakage (REL). Less H₂O₂ and O₂^-, higher activities of POD and SOD, less MDA content and REL, and higher chlorophyll content and photosynthetic rate were also shown in transgenic plants after inoculation with P. infestans. Our results constitute the first step towards further investigations into the biological function and molecular mechanism of miR172-mediated silencing of AP2/ERF transcription factors in S. lycopersicum-P. infestans interaction and provide a candidate gene for breeding to enhance biotic stress-resistance in S. lycopersicum.

Cross-Kingdom Regulation of Putative miRNAs Derived from Happy Tree in Cancer Pathway: A Systems Biology Approach

MicroRNAs (miRNAs) are well-known key regulators of gene expression primarily at the post-transcriptional level. Plant-derived miRNAs may pass through the gastrointestinal tract, entering into the body fluid and regulate the expression of endogenous mRNAs. Camptotheca acuminata, a highly important medicinal plant known for its anti-cancer potential was selected to investigate cross-kingdom regulatory mechanism and involvement of miRNAs derived from this plant in cancer-associated pathways through in silico systems biology approach. In this study, total 33 highly stable putative novel miRNAs were predicted from the publically available 53,294 ESTs of C. acuminata, out of which 14 miRNAs were found to be regulating 152 target genes in human. Functional enrichment, gene-disease associations and network analysis of these target genes were carried out and the results revealed their association with prominent types of cancers like breast cancer, leukemia and lung cancer. Pathways like focal adhesion, regulation of lipolysis in adipocytes and mTOR signaling pathways were found significantly associated with the target genes. The regulatory network analysis showed the association of some important hub proteins like GSK3B, NUMB, PEG3, ITGA2 and DLG2 with cancer-associated pathways. Based on the analysis results, it can be suggested that the ingestion of the C. acuminata miRNAs may have a functional impact on tumorigenesis in a cross-kingdom way and may affect the physiological condition at genetic level. Thus, the predicted miRNAs seem to hold potentially significant role in cancer pathway regulation and therefore, may be further validated using in vivo experiments for a better insight into their mechanism of epigenetic action of miRNA.

Uncovering leaf rust responsive miRNAs in wheat (Triticum aestivum L.) using high-throughput sequencing and prediction of their targets through degradome analysis

Deep sequencing identified 497 conserved and 559 novel miRNAs in wheat, while degradome analysis revealed 701 targets genes. QRT-PCR demonstrated differential expression of miRNAs during stages of leaf rust progression. Bread wheat (Triticum aestivum L.) is an important cereal food crop feeding 30 % of the world population. Major threat to wheat production is the rust epidemics. This study was targeted towards identification and functional characterizations of micro(mi)RNAs and their target genes in wheat in response to leaf rust ingression. High-throughput sequencing was used for transcriptome-wide identification of miRNAs and their expression profiling in retort to leaf rust using mock and pathogen-inoculated resistant and susceptible near-isogenic wheat plants. A total of 1056 mature miRNAs were identified, of which 497 miRNAs were conserved and 559 miRNAs were novel. The pathogen-inoculated resistant plants manifested more miRNAs compared with the pathogen infected susceptible plants. The miRNA counts increased in susceptible isoline due to leaf rust, conversely, the counts decreased in the resistant isoline in response to pathogenesis illustrating precise spatial tuning of miRNAs during compatible and incompatible interaction. Stem-loop quantitative real-time PCR was used to profile 10 highly differentially expressed miRNAs obtained from high-throughput sequencing data. The spatio-temporal profiling validated the differential expression of miRNAs between the isolines as well as in retort to pathogen infection. Degradome analysis provided 701 predicted target genes associated with defense response, signal transduction, development, metabolism, and transcriptional regulation. The obtained results indicate that wheat isolines employ diverse arrays of miRNAs that modulate their target genes during compatible and incompatible interaction. Our findings contribute to increase knowledge on roles of microRNA in wheat-leaf rust interactions and could help in rust resistance breeding programs.

MiR1918 enhances tomato sensitivity to Phytophthora infestans infection

Late blight of tomato is caused by the oomycete pathogen Phytophthora infestans. In our previous work, we identified and characterized a miR1918 in P. infestans (pi-miR1918), and showed that its sequence is similar to the sequence of tomato miR1918 (sly-miR1918). In this study, we used Arabidopsis thaliana pre-miR159a as a backbone to synthesize pi-miR1918 via PCR and mutagenesis. The artificial pi-miR1918 was used to investigate the role of miR1918 in tomato-P. infestans interaction. Trangenic tomato plants that overexpressed the artificial pi-miR1918 displayed more serious disease symptoms than wild-type tomato plants after infection with P. infestans, as shown by increased number of necrotic cells, lesion sizes and number of sporangia per leaf. The target genes of pi-miR1918 and sly-miR1918 were also predicted for tomato and P. infestans, respectively. qPCR analysis of these targets also performed during tomato-P. infestans interaction. The expression of target gene, RING finger were negatively correlated with miR1918 in the all Lines of transgenic tomato plants. In addition, we used the 5′ RACE to determine the cleavage site of miR1918 to RING finger. These results suggested that miR1918 might be involved in the silencing of target genes, thereby enhancing the susceptibility of tomato to P. infestans infection.

Small RNA and degradome deep sequencing reveals drought-and tissue-specific micrornas and their important roles in drought-sensitive and drought-tolerant tomato genotypes

Drought stress has adverse impacts on plant production and productivity. MicroRNAs (miRNAs) are one class of noncoding RNAs regulating gene expression post-transcriptionally. In this study, we employed small RNA and degradome sequencing to systematically investigate the tissue-specific miRNAs responsible to drought stress, which are understudied in tomato. For this purpose, root and upground tissues of two different drought-responsive tomato genotypes (Lycopersicon esculentum as sensitive and L. esculentum var. cerasiforme as tolerant) were subjected to stress with 5% polyethylene glycol for 7 days. A total of 699 conserved miRNAs belonging to 578 families were determined and 688 miRNAs were significantly differentially expressed between different treatments, tissues and genotypes. Using degradome sequencing, 44 target genes were identified associated with 36 miRNA families. Drought-related miRNAs and their targets were enriched functionally by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses. Totally, 53 miRNAs targeted 23 key drought stress- and tissue development-related genes, including DRP (dehydration-responsive protein), GTs (glycosyltransferases), ERF (ethylene responsive factor), PSII (photosystem II) protein, HD-ZIP (homeodomain-leucine zipper), MYB and NAC-domain transcription factors. miR160, miR165, miR166, miR171, miR398, miR408, miR827, miR9472, miR9476 and miR9552 were the key miRNAs functioning in regulation of these genes and involving in tomato response to drought stress. Additionally, plant hormone signal transduction pathway genes were differentially regulated by miR169, miR172, miR393, miR5641, miR5658 and miR7997 in both tissues of both sensitive and tolerant genotypes. These results provide new insight into the regulatory role of miRNAs in drought response with plant hormone signal transduction and drought-tolerant tomato breeding.

Identification of Arbuscular Mycorrhiza (AM)-Responsive microRNAs in Tomato

A majority of land plants can form symbiosis with arbuscular mycorrhizal (AM) fungi. MicroRNAs (miRNAs) have been implicated to regulate this process in legumes, but their involvement in non-legume species is largely unknown. In this study, by performing deep sequencing of sRNA libraries in tomato roots and comparing with tomato genome, a total of 700 potential miRNAs were predicted, among them, 187 are known plant miRNAs that have been previously deposited in miRBase. Unlike the profiles in other plants such as rice and Arabidopsis, a large proportion of predicted tomato miRNAs was 24 nt in length. A similar pattern was observed in the potato genome but not in tobacco, indicating a Solanum genus-specific expansion of 24-nt miRNAs. About 40% identified tomato miRNAs showed significantly altered expressions upon Rhizophagus irregularis inoculation, suggesting the potential roles of these novel miRNAs in AM symbiosis. The differential expression of five known and six novel miRNAs were further validated using qPCR analysis. Interestingly, three up-regulated known tomato miRNAs belong to a known miR171 family, a member of which has been reported in Medicago truncatula to regulate AM symbiosis. Thus, the miR171 family likely regulates AM symbiosis conservatively across different plant lineages. More than 1000 genes targeted by potential AM-responsive miRNAs were provided and their roles in AM symbiosis are worth further exploring.

Identification of Novel and Conserved miRNAs in Leaves of In vitro Grown Citrus reticulata "Lugan" Plantlets by Solexa Sequencing

MicroRNAs (miRNAs) play essential roles in plant development, but the roles in the in vitro plant development are unknown. Leaves of ponkan plantlets derived from mature embryos at in vitro culture conditions were used to sequence small RNA fraction via Solexa sequencing, and the miRNAs expression was analyzed. The results showed that there were 3,065,625 unique sequences in ponkan, of which 0.79% were miRNAs. The RNA sequences with lengths of 18-25 nt derived from the library were analyzed, leading to the identification of 224 known miRNAs, of which the most abundant were miR157, miR156, and miR166. Three hundred and fifty-eight novel miRNA candidates were also identified, and the number of reads of ponkan novel miRNAs varied from 5 to 168,273. The expression of the most known miRNAs obtained was at low levels, which varied from 5 to 4,946,356. To better understand the role of miRNAs during the preservation of ponkan in vitro plantlet, the expression patterns of cre-miR156a/159b/160a/166a/167a/168a/171/398b were validated by quantitative real-time PCR (qPCR). The results showed that not only the development-associated miRNAs, e.g., cre-miR156/159/166/396, expressed highly at the early preservation period in the in vitro ponkan plantlet leaves but also the stress-related miRNAs, e.g., cre-miR171 and cre-miR398b, expressed highly at the same time. The expression levels of most tested miRNAs were found to decrease after 6 months and the amounts of these miRNAs were kept at low levels at 18 months. After analyzing the expression level of their targets during the reservation of the ponkan in vitro plantlet, development-associated cre-ARF6 and stress-related cre-CSD modules exhibited negative correlation with miR167 and miR398, respectively, indicating an involvement of the miRNAs in the in vitro development of ponkan and function in the conservation of ponkan germplasm.

Sp-miR396a-5p acts as a stress-responsive genes regulator by conferring tolerance to abiotic stresses and susceptibility to Phytophthora nicotianae infection in transgenic tobacco

KEY MESSAGE

Overexpression of Sp-miR396a-5p in tobacco increased tolerance to salt, drought, cold stress and susceptibility to Phytophthora nicotianae infection. MicroRNA396 (miR396) is one of the conserved microRNA families in plants, and it targeted growth-regulating factors (GRFs) family. The GRF transcription factors are associated with growth and stress responses. However, the molecular mechanisms of miR396 responding to environmental stresses are elusive. The purpose of this study was to explore the function of tomato miR396a-5p (Sp-miR396a-5p) in Solanaceae responses to abiotic and biotic stresses. We showed that Sp-miR396a-5p transcript levels were up-regulated under salt and drought stresses and down-regulated after Phytophthora infestans (P. infestans) infection. Consistently, overexpression of Sp-miR396a-5p in tobacco enhanced its tolerance to salt, drought and cold stresses. Additionally, the expression of Sp-miR396a-5p was found to be down-regulated under pathogen-related biotic stress. Tobacco plants overexpressing Sp-miR396a-5p showed increased susceptibility to Phytophthora nicotianae (P. nicotianae) infection. Physiological analysis indicated that Sp-miR396a-5p overexpression enhanced osmoregulation and decreased production of reactive oxygen species (ROS). Furthermore, four Sp-miR396a-5p target genes, NtGRF1, NtGRF3, NtGRF7 and NtGRF8, were down-regulated in these plants. Our results suggested that Sp-miR396a-5p plays critical roles in both abiotic stresses through targeting NtGRF7-regulated expression of osmotic stress-responsive genes and pathogen infection via the regulatory networks of NtGRF1 and NtGRF3.

The advance of tomato disease-related microRNAs

Tomato is a model plant for studying plant-pathogen interactions. As regulatory factors, microRNAs (miRNAs) have been widely identified and play crucial roles in tomato-pathogen interactions, including host defense and pathogen counter-defense. Here, the review summarizes the discoveries and highlights of miRNAs in tomato diseases. Roles of artificial miRNAs in disease resistance are further discussed. Hence, a better understanding of the contribution of miRNAs in tomato disease will shed light on strategies in enhancing tomato-pathogen resistance.

High-throughput sequencing reveals differential expression of miRNAs in tomato inoculated with Phytophthora infestans

The characterization and compare expression profiling of the miRNA transcriptome lay a solid foundation for unraveling the complex miRNA-mediated regulatory network in tomato resistance mechanisms against LB. MicroRNAs (miRNAs) are a class of small endogenous non-coding RNAs with 20-24 nt. They have been identified in many plants with their diverse regulatory roles in biotic stresses. The knowledge, that miRNAs regulate late blight (LB), caused by Phytophthora infestans, is rather limited. In this study, we used miRNA-Seq to investigate the miRNA expression difference between the tomatoes treated with and without P. infestans. A total of 42,714,516 raw reads were generated from two small RNA libraries by high-throughput sequencing. Finally, 207 known miRNAs and 67 new miRNAs were obtained. The differential expression profile of miRNAs in tomato was further analyzed with twofold change (P value ≤0.01). A total of 70 miRNAs were manifested to change significantly in samples treated with P. infestans, including 50 down-regulated miRNAs and 20 up-regulated miRNAs. Moreover, a total of 73 target genes were acquired for 28 differentially expressed miRNAs by psRNATarget analysis. By enrichment pathway analysis of target genes, plant-pathogen interaction was the most highly relevant pathway which played an important role in disease defense. In addition, 30 miRNAs were selected for qRT-PCR to validate their expression patterns. The expression patterns for targets of miR6027, miR5300, miR476b, miR159a, miR164a and miRn13 were selectively examined, and the results showed that there was a negative correlation on the expression patterns between miRNAs and their targets. The targets have previously been reported to be related with plant immune and involved in plant-pathogen interaction pathway in this study, suggesting these miRNAs might act as regulators in process of tomato resistance against P. infestans. These discoveries will provide us useful information to explain tomato resistance mechanisms against LB.