Prediction and characterization of Tomato leaf curl New Delhi virus (ToLCNDV) responsive novel microRNAs in Solanum lycopersicum

Profiling of Groundnut bud necrosis orthotospovirus-responsive microRNA and their targets in tomato based on deep sequencing

Tomato, an extensively cultivated vegetable crop produces miRNAs in response to infection with Groundnut bud necrosis orthotospovirus, a viral pathogen causing significant economic losses. High-throughput miRNA sequencing was performed on tomato leaves inoculated with GBNV and mock-inoculated leaves as controls. Analysis revealed 73 known miRNAs belonging to 24 miRNA families, with variable expression levels. Interestingly, 39 miRNAs were upregulated, and 34 were downregulated in response to GBNV infection. Stem-loop quantitative reverse transcription PCR validated the differential expression of selected miRNAs. Additionally, 30 miRNA encoded proteins were identified to be involved in disease resistance and susceptibility. The miRNA-target interactions were found to play significant roles in cellular and metabolic activities, as well as modulating signaling pathways during the plant-virus interaction. The findings shed light on the intricate regulatory network of miRNAs in tomato response to viral infection and may contribute to developing strategies for improving crop protection against viral diseases.

Transcriptional and epigenetic changes during tomato yellow leaf curl virus infection in tomato

BACKGROUND

Geminiviruses are DNA plant viruses that cause highly damaging diseases affecting crops worldwide. During the infection, geminiviruses hijack cellular processes, suppress plant defenses, and cause a massive reprogramming of the infected cells leading to major changes in the whole plant homeostasis. The advances in sequencing technologies allow the simultaneous analysis of multiple aspects of viral infection at a large scale, generating new insights into the molecular mechanisms underlying plant-virus interactions. However, an integrative study of the changes in the host transcriptome, small RNA profile and methylome during a geminivirus infection has not been performed yet. Using a time-scale approach, we aim to decipher the gene regulation in tomato in response to the infection with the geminivirus, tomato yellow leaf curl virus (TYLCV).

RESULTS

We showed that tomato undergoes substantial transcriptional and post-transcriptional changes upon TYLCV infection and identified the main altered regulatory pathways. Interestingly, although the principal plant defense-related processes, gene silencing and the immune response were induced, this cannot prevent the establishment of the infection. Moreover, we identified extra- and intracellular immune receptors as targets for the deregulated microRNAs (miRNAs) and established a network for those that also produced phased secondary small interfering RNAs (phasiRNAs). On the other hand, there were no significant genome-wide changes in tomato methylome at 14 days post infection, the time point at which the symptoms were general, and the amount of viral DNA had reached its maximum level, but we were able to identify differentially methylated regions that could be involved in the transcriptional regulation of some of the differentially expressed genes.

CONCLUSION

We have conducted a comprehensive and reliable study on the changes at transcriptional, post-transcriptional and epigenetic levels in tomato throughout TYLCV infection. The generated genomic information is substantial for understanding the genetic, molecular and physiological changes caused by TYLCV infection in tomato.

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.

The Sw5a gene confers resistance to ToLCNDV and triggers an HR response after direct AC4 effector recognition

Several attempts have been made to identify antiviral genes against Tomato leaf curl New Delhi virus (ToLCNDV) and related viruses. This has led to the recognition of Ty genes (Ty1-Ty6), which have been successful in developing virus-resistant crops to some extent. Owing to the regular appearance of resistance-breaking strains of these viruses, it is important to identify genes related to resistance. In the present study, we identified a ToLCNDV resistance (R) gene, SlSw5a, in a ToLCNDV-resistant tomato cultivar, H-88-78-1, which lacks the known Ty genes. The expression of SlSw5a is controlled by the transcription factor SlMyb33, which in turn is regulated by microRNA159 (sly-miR159). Virus-induced gene silencing of either SlSw5a or SlMyb33 severely increases the disease symptoms and viral titer in leaves of resistant cultivar. Moreover, in SlMyb33-silenced plants, the relative messenger RNA level of SlSw5a was reduced, suggesting SlSw5a is downstream of the sly-miR159-SlMyb33 module. We also demonstrate that SlSw5a interacts physically with ToLCNDV-AC4 (viral suppressor of RNA silencing) to trigger a hypersensitive response (HR) and generate reactive oxygen species at infection sites to limit the spread of the virus. The "RTSK" motif in the AC4 C terminus is important for the interaction, and its mutation completely abolishes the interaction with Sw5a and HR elicitation. Overall, our research reports an R gene against ToLCNDV and establishes a connection between the upstream miR159-Myb33 module and its downstream target Sw5a to activate HR in the tomato, resulting in geminivirus resistance.

Tuning Beforehand: A Foresight on RNA Interference (RNAi) and In Vitro-Derived dsRNAs to Enhance Crop Resilience to Biotic and Abiotic Stresses

Crop yield is severely affected by biotic and abiotic stresses. Plants adapt to these stresses mainly through gene expression reprogramming at the transcriptional and post-transcriptional levels. Recently, the exogenous application of double-stranded RNAs (dsRNAs) and RNA interference (RNAi) technology has emerged as a sustainable and publicly acceptable alternative to genetic transformation, hence, small RNAs (micro-RNAs and small interfering RNAs) have an important role in combating biotic and abiotic stresses in plants. RNAi limits the transcript level by either suppressing transcription (transcriptional gene silencing) or activating sequence-specific RNA degradation (post-transcriptional gene silencing). Using RNAi tools and their respective targets in abiotic stress responses in many crops is well documented. Many miRNAs families are reported in plant tolerance response or adaptation to drought, salinity, and temperature stresses. In biotic stress, the spray-induced gene silencing (SIGS) provides an intelligent method of using dsRNA as a trigger to silence target genes in pests and pathogens without producing side effects such as those caused by chemical pesticides. In this review, we focus on the potential of SIGS as the most recent application of RNAi in agriculture and point out the trends, challenges, and risks of production technologies. Additionally, we provide insights into the potential applications of exogenous RNAi against biotic stresses. We also review the current status of RNAi/miRNA tools and their respective targets on abiotic stress and the most common responsive miRNA families triggered by stress conditions in different crop species.

Integrated Analysis of microRNA and mRNA Transcriptome Reveals the Molecular Mechanism of Solanum lycopersicum Response to Bemisia tabaci and Tomato chlorosis virus

Tomato chlorosis virus (ToCV), is one of the most devastating cultivated tomato viruses, seriously threatened the growth of crops worldwide. As the vector of ToCV, the whitefly Bemisia tabaci Mediterranean (MED) is mainly responsible for the rapid spread of ToCV. The current understanding of tomato plant responses to this virus and B. tabaci is very limited. To understand the molecular mechanism of the interaction between tomato, ToCV and B. tabaci, we adopted a next-generation sequencing approach to decipher miRNAs and mRNAs that are differentially expressed under the infection of B. tabaci and ToCV in tomato plants. Our data revealed that 6199 mRNAs were significantly regulated, and the differentially expressed genes were most significantly associated with the plant-pathogen interaction, the MAPK signaling pathway, the glyoxylate, and the carbon fixation in photosynthetic organisms and photosynthesis related proteins. Concomitantly, 242 differentially expressed miRNAs were detected, including novel putative miRNAs. Sly-miR159, sly-miR9471b-3p, and sly-miR162 were the most expressed miRNAs in each sample compare to control group. Moreover, we compared the similarities and differences of gene expression in tomato plant caused by infection or co-infection of B. tabaci and ToCV. Taken together, the analysis reported in this article lays a solid foundation for further research on the interaction between tomato, ToCV and B. tabaci, and provide evidence for the identification of potential key genes that influences virus transmission in tomato plants.

MIR172d Is Required for Floral Organ Identity and Number in Tomato

MicroRNA172 (miR172) functions as a central regulator of flowering time and flower development by post-transcriptional repression of APETALA2-LIKE transcription factors. In the model crop Solanum lycopersicum (tomato), the miR172 family is still poorly annotated and information about the functions of specific members is lacking. Here, de-novo prediction of tomato miR172 coding loci identified seven genes (SlMIR172a-g), that code for four unique species of miR172 (sly-miR172). During reproductive development, sly-miR172s are differentially expressed, with sly-miR172c and sly-miR172d being the most abundant members in developing flowers, and are predicted to guide the cleavage of eight APETALA2-LIKE transcription factors. By CRISPR-Cas9 co-targeting of SlMIR172c and SlMIR172d we have generated a battery of loss-of-function and hypomorphic mutants (slmir172c-d^CR). The slmir172c-d^CR plants developed normal shoot but their flowers displayed graded floral organ abnormalities. Whereas slmir172c^CR loss-of-function caused only a slight greening of petals and stamens, hypomorphic and loss-of-function slmir172d^CR alleles were associated with the conversion of petals and stamens to sepaloids, which were produced in excess. Interestingly, the degrees of floral organ identity alteration and proliferation were directly correlated with the reduction in sly-miR172d activity. These results suggest that sly-miR172d regulates in a dose-dependent manner floral organ identity and number, likely by negatively regulating its APETALA2-like targets.

Differential miRNA profiles in South African cassava mosaic virus-infected cassava landraces reveal clues to susceptibility and tolerance to cassava mosaic disease

Specific miRNA families are involved in susceptibility or antiviral immunity in plants. Manihot esculenta Crantz (cassava) is a perennial plant that is an important food security crop in sub-Saharan Africa. Cassava is susceptible to several begomoviruses that cause cassava mosaic disease (CMD). In this study, we investigated the leaf miRNAome response in a tolerant (TME3) and susceptible (T200) cassava landrace challenged with South African cassava mosaic virus. RNAseq was performed on leaf samples at 12, 32 and 67 days post infection (dpi), representing early, symptomatic and late persistent stages of CMD infection. Significantly, distinct profiles of conserved miRNA family expression between the T200 and TME3 landraces at the three infection stages were observed. Notably at 12 days post SACMV infection, TME3 exhibited significant downregulation (log2fold<2.0) of 42 %, compared to 9% in T200, of the conserved miRNA families. This demonstrates an overall early response to SACMV in TME3 prior to symptom appearance not observed in T200, and expression of a large cohort of miRNA-regulated genes. Notably, at early infection, downregulation of mes-miR162 and 168 that target antiviral posttransriptional gene silencing (PTGS) regulators DCL1 and AGO1, respectively, was observed in TME3, and AGO1 and DCL1 expression was higher compared to T200 post infection. Early rapid responses prior to symptom development, including RNA silencing, may be key to establishing the tolerance/recovery phenotype exhibited by TME3 landrace later on at 67 dpi. At recovery, TME3 was hallmarked by a highly significant down-regulation of mes-miR167. MiR167 targets an auxin responsive factor which plays a role in auxin signaling and adaptive responses to stress, suggesting the importance of the auxin signaling in recovery of SACMV-induced symptoms. The gene targets of these miRNAs and their associated networks may provide clues to the molecular basis of CMD tolerance in perennial hosts such as cassava.

Virus and Viroid-Derived Small RNAs as Modulators of Host Gene Expression: Molecular Insights Into Pathogenesis

Virus-derived siRNAs (vsiRNAs) generated by the host RNA silencing mechanism are effectors of plant’s defense response and act by targeting the viral RNA and DNA in post-transcriptional gene silencing (PTGS) and transcriptional gene silencing (TGS) pathways, respectively. Contrarily, viral suppressors of RNA silencing (VSRs) compromise the host RNA silencing pathways and also cause disease-associated symptoms. In this backdrop, reports describing the modulation of plant gene(s) expression by vsiRNAs via sequence complementarity between viral small RNAs (sRNAs) and host mRNAs have emerged. In some cases, silencing of host mRNAs by vsiRNAs has been implicated to cause characteristic symptoms of the viral diseases. Similarly, viroid infection results in generation of sRNAs, originating from viroid genomic RNAs, that potentially target host mRNAs causing typical disease-associated symptoms. Pathogen-derived sRNAs have been demonstrated to have the propensity to target wide range of genes including host defense-related genes, genes involved in flowering and reproductive pathways. Recent evidence indicates that vsiRNAs inhibit host RNA silencing to promote viral infection by acting as decoy sRNAs. Nevertheless, it remains unclear if the silencing of host transcripts by viral genome-derived sRNAs are inadvertent effects due to fortuitous pairing between vsiRNA and host mRNA or the result of genuine counter-defense strategy employed by viruses to enhance its survival inside the plant cell. In this review, we analyze the instances of such cross reaction between pathogen-derived vsiRNAs and host mRNAs and discuss the molecular insights regarding the process of pathogenesis.

Novel targets for engineering Physostegia chlorotic mottle and tomato brown rugose fruit virus-resistant tomatoes: in silico prediction of tomato microRNA targets

Background

Physostegia chlorotic mottle virus (PhCMoV; genus: Alphanucleorhabdovirus, family: Rhabdoviridae) and tomato brown rugose fruit virus (ToBRFV; genus: Tobamovirus, family: Virgaviridae) are newly emerging plant viruses that have a dramatic effect on tomato production. Among various known virus-control strategies, RNAi-mediated defence has shown the potential to protect plants against various pathogens including viral infections. Micro(mi)RNAs play a major role in RNAi-mediated defence.

Methods

Using in silico analyses, we investigated the possibility of tomato-encoded miRNAs (TomiRNA) to target PhCMoV and ToBRFV genomes using five different algorithms, i.e., miRanda, RNAhybrid, RNA22, Tapirhybrid and psRNATarget.

Results

The results revealed that 14 loci on PhCMoV and 10 loci on ToBRFV can be targeted by the TomiRNAs based on the prediction of at least three algorithms. Interestingly, one TomiRNA, miR6026, can target open reading frames from both viruses, i.e., the phosphoprotein encoding gene of PhCMoV, and the two replicase components of ToBRFV. There are currently no commercially available PhCMoV- or ToBRFV-resistant tomato varieties, therefore the predicted data provide useful information for the development of PhCMoV- and ToBFRV-resistant tomato plants.

Potato Virus Y Infection Alters Small RNA Metabolism and Immune Response in Tomato

Potato virus Y (PVY) isolate PVY^C-to induces growth reduction and foliar symptoms in tomato, but new vegetation displays symptom recovery at a later stage. In order to investigate the role of micro(mi)RNA and secondary small(s)RNA-regulated mechanisms in tomato defenses against PVY, we performed sRNA sequencing from healthy and PVY^C-to infected tomato plants at 21 and 30 days post-inoculation (dpi). A total of 792 miRNA sequences were obtained, among which were 123 canonical miRNA sequences, many isomiR variants, and 30 novel miRNAs. MiRNAs were mostly overexpressed in infected vs. healthy plants, whereas only a few miRNAs were underexpressed. Increased accumulation of isomiRs was correlated with viral infection. Among miRNA targets, enriched functional categories included resistance (R) gene families, transcription and hormone factors, and RNA silencing genes. Several 22-nt miRNAs were shown to target R genes and trigger the production of 21-nt phased sRNAs (phasiRNAs). Next, 500 phasiRNA-generating loci were identified, and were shown to be mostly active in PVY-infected tissues and at 21 dpi. These data demonstrate that sRNA-regulated host responses, encompassing miRNA alteration, diversification within miRNA families, and phasiRNA accumulation, regulate R and disease-responsive genes. The dynamic regulation of miRNAs and secondary sRNAs over time suggests a functional role of sRNA-mediated defenses in the recovery phenotype.

Characterization and expression profiling of microRNAs in response to plant feeding in two host-plant strains of the lepidopteran pest Spodoptera frugiperda

Background

A change in the environment may impair development or survival of living organisms leading them to adapt to the change. The resulting adaptation trait may reverse, or become fixed in the population leading to evolution of species. Deciphering the molecular basis of adaptive traits can thus give evolutionary clues. In phytophagous insects, a change in host-plant range can lead to emergence of new species. Among them, Spodoptera frugiperda is a major agricultural lepidopteran pest consisting of two host-plant strains having diverged 3 MA, based on mitochondrial markers. In this paper, we address the role of microRNAs, important gene expression regulators, in response to host-plant change and in adaptive evolution.

Results

Using small RNA sequencing, we characterized miRNA repertoires of the corn (C) and rice (R) strains of S. frugiperda, expressed during larval development on two different host-plants, corn and rice, in the frame of reciprocal transplant experiments. We provide evidence for 76 and 68 known miRNAs in C and R strains and 139 and 171 novel miRNAs. Based on read counts analysis, 34 of the microRNAs were differentially expressed in the C strain larvae fed on rice as compared to the C strain larvae fed on corn. Twenty one were differentially expressed on rice compared to corn in R strain. Nine were differentially expressed in the R strain compared to C strain when reared on corn. A similar ratio of microRNAs was differentially expressed between strains on rice. We could validate experimentally by QPCR, variation in expression of the most differentially expressed candidates. We used bioinformatics methods to determine the target mRNAs of known microRNAs. Comparison with the mRNA expression profile during similar reciprocal transplant experiment revealed potential mRNA targets of these host-plant regulated miRNAs.

Conclusions

In the current study, we performed the first systematic analysis of miRNAs in Lepidopteran pests feeding on host-plants. We identified a set of the differentially expressed miRNAs that respond to the plant diet, or differ constitutively between the two host plant strains. Among the latter, the ones that are also deregulated in response to host-plant are molecular candidates underlying a complex adaptive trait.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5119-6) contains supplementary material, which is available to authorized users.

Suppression of nbe-miR166h-p5 attenuates leaf yellowing symptoms of potato virus X on Nicotiana benthamiana and reduces virus accumulation

MicroRNAs (miRNAs) play essential roles in plant development. There is increasing evidence that changed expression of miRNAs in virus-infected plants contributes to the development of viral symptoms. Here, we analysed the altered expression of miRNAs of Nicotiana benthamiana in response to Potato virus X (PVX) by Illumina Solexa sequencing. One of the 21 miRNAs significantly affected, nbe-miR166h-p5, was closely associated with viral symptoms. Using the Tobacco rattle virus-based miRNA suppression (VbMS) system, we found that the suppression of nbe-miR166h-p5 in plants caused leaves to turn dark green with increased chlorophyll. When PVX was inoculated on nbe-miR166h-p5-suppressed plants, the leaf yellowing symptom of PVX was largely attenuated with less reduction in chlorophyll content, and the accumulation of PVX was decreased. nbe-miR166h-p5 was also up-regulated in plants infected by Turnip mosaic virus (TuMV), and its suppression attenuated the leaf yellowing symptom of TuMV and decreased viral accumulation. Three potential targets of nbe-miR166h-p5 were identified. The results indicate the association of nbe-miR166h-p5 with symptoms of PVX and also with those of TuMV, providing useful information on the relationship between miRNA and viral infection.

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.

Identification and comparative analysis of microRNAs from tomato varieties showing contrasting response to ToLCV infections

Increasing incidence of viral infections in crop plants adversely affects their growth and yield. Tomato (Solanum lycopersicum) is considered to be a favorite host for viruses with over 50 species of begomoviruses naturally infecting this crop. Tomato leaf curl virus (ToLCV) is among the most widespread and devastating begomoviruses affecting tomato production. microRNAs (miRs) have been established as key regulators of gene expression and plant development. The miR pathways are disturbed during infection by viruses. Thus, comprehension of regulatory miR networks is crucial in understanding the effect of viral pathogenicity. To identify key miRs involved in ToLCV infection, a high throughput approach involving next generation sequencing was employed. Healthy and infected leaf tissues of two tomato varieties, differing in their susceptibility to ToLCV infection were analyzed. NGS data analysis followed by computational predictions, led to identification of 91 known miRs, 15 novel homologs and 53 novel miRs covering two different varieties of tomato, susceptible (Pusa Ruby) and tolerant (LA1777) to ToLCV infection. The cleaved targets of these miRs were identified using online available degradome libraries from leaf, flower and fruit of tomato and showed their involvement in various biological pathways through KEGG Orthology. With detailed comparative profiling of expression pattern of these miRs, we could associate the specific miRs with the resistant and infected genotypes. This study depicted that in depth analysis of miR expression patterns and their functions will help in identification of molecules that can be used for manipulation of gene expression to increase crop production and developing resistance against diseases.

Resistance to tomato leaf curl New Delhi virus in melon is controlled by a major QTL located in chromosome 11

Identification of three genomic regions and underlying candidate genes controlling the high level of resistance to ToLCNDV derived from a wild melon. SNP markers appropriated for MAS management of resistance. Tomato leaf curl New Delhi virus (ToLCNDV) is a bipartite begomovirus that severely affects melon crop (Cucumis melo) in the main production areas of Spain since 2012. In this work, we evaluated the degree of resistance of four accessions (two belonging to the subsp. agrestis var. momordica and two to the wild agrestis group) and their corresponding hybrids with a susceptible commercial melon belonging to the subsp. melo (Piel de Sapo, PS). The analysis using quantitative PCR (qPCR) allowed us to select one wild agrestis genotype (WM-7) with a high level of resistance and use it to construct segregating populations (F ₂ and backcrosses). These populations were phenotyped for symptom severity and virus content using qPCR, and genotyped with different sets of SNP markers. Phenotyping and genotyping results in the F ₂ and BC1s populations derived from the WM-7 × PS cross were used for QTL analysis. Three genomic regions controlling resistance to ToLCNDV were found, one major locus in chromosome 11 and two additional regions in chromosomes 12 and 2. The highest level of resistance (no or mild symptoms and very low viral titer) was obtained with the homozygous WM-7WM-7 genotype at the major QTL in chromosome 11, even with PSPS genotypes at the other two loci. The resistance derived from WM-7 is useful to develop new melon cultivars and the linked SNPs selected in this paper will be highly useful in marker-assisted breeding for ToLCNDV resistance in melon.

Identification of drought-responsive microRNAs in tomato using high-throughput sequencing

Drought is a major abiotic stress affecting crop productivity and quality. As a class of noncoding RNA, microRNA (miRNA) plays important roles in plant growth, development, and stress response. However, their response and roles in tomato drought stress is largely unknown. Here, by using high-throughput sequencing, we compared the miRNA profiles before and after drought treatment in two tomato genotypes: M82, a drought-sensitive cultivated tomato (Solanum lycopersicum), and IL2-5, a drought-tolerant introgression line derived from M82 and the tomato wild species S. pennellii (LA0716). A total of 108 conserved and 208 novel miRNAs were identified, among them, 32 and 68 were significantly changed in expression after stress. Further, 1936 putative target genes were predicted for those differentially-expressed miRNAs. Gene ontology and pathway analysis showed that many of the target genes were involved in stress resistance, such as genes in GO terms including response to stress, defense response, response to stimulus, phosphorylation, and signal transduction. Our results suggested that miRNAs play an essential role in the drought response of tomato. This work will help to further characterize specific miRNAs functioning in drought tolerance.

Tomato leaf curl New Delhi virus: a widespread bipartite begomovirus in the territory of monopartite begomoviruses

Tomato leaf curl New Delhi virus (ToLCNDV) is an exceptional Old World bipartite begomovirus. On the Indian subcontinent, a region in which monopartite DNA satellite-associated begomoviruses with mostly narrow geographical ranges predominate, it is widespread, with a geographical range also including the Far East, Middle East, North Africa and Europe. The success of ToLCNDV probably derives from its broad host range and highly flexible genomic configuration: its DNA-A component is capable of productively interacting with, and trans-replicating, diverse DNA-B components and betasatellites. An understanding of the capacity of ToLCNDV to infect a variety of hosts and spread across a broad and ecologically variable geographical range could illuminate the potential economic threats associated with similar begomoviral invasions. Towards this end, we used available ToLCNDV sequences to reconstruct the history of ToLCNDV spread.

TAXONOMY

Family Geminiviridae, Genus Begomovirus. ToLCNDV is a bipartite begomovirus. Following the revised begomovirus taxonomic criteria of 91% and 94% nucleotide identity for species and strain demarcation, respectively, ToLCNDV is a distinct species with two strains: ToLCNDV and ToLCNDV-Spain.

HOST RANGE

The primary cultivated host of ToLCNDV is tomato (Solanum lycopersicum), but the virus is also known to infect 43 other plant species from a range of families, including Cucurbitaceae, Euphorbiaceae, Solanaceae, Malvaceae and Fabaceae.

DISEASE SYMPTOMS

Typical symptoms of ToLCNDV infection in its various hosts include leaf curling, vein thickening, puckering, purpling/darkening of leaf margins, leaf area reduction, internode shortening and severe stunting.

Integrated miRNA and mRNA expression profiling reveals the response regulators of a susceptible tomato cultivar to early blight disease

Early blight, caused by the fungus Alternaria solani, is a devastating foliar disease of tomatoes, causes massive yield loss each year worldwide. Molecular basis of the compatible host–pathogen interaction was elusive. We adopted next generation sequencing approach to decipher miRNAs and mRNAs that are differentially expressed during Alternaria-stress in tomato. Some of the interesting findings were also validated by alternative techniques. Our analysis revealed 181 known-miRNAs, belonging to 121 miRNA families, of which 67 miRNAs showed at least 2-fold change in expression level with the majority being downregulated. Concomitantly, 5,450 mRNAs were significantly regulated in the same diseased tissues. Differentially expressed genes were most significantly associated with response to stimulus process, photosynthesis, biosynthesis of secondary metabolites, plant–pathogen interaction and plant hormone signal transduction pathways. GO term enrichment-based categorization of gene-functions further supported this observation, as terms related to pathogen perception, disease signal transduction, cellular metabolic processes including oxidoreductase and kinase activity were over represented. In addition, we have discovered 102 miRNA–mRNA pairs which were regulated antagonistically, and careful study of the targeted mRNAs depicted that multiple transcription factors, nucleotide-binding site leucine-rich repeats, receptor-like proteins and enzymes related to cellular ROS management were profoundly affected. These studies have identified key regulators of Alternaria-stress response in tomato and the subset of genes that are likely to be post-transcriptionally silenced during the infection.

Altered accumulation of osa-miR171b contributes to rice stripe virus infection by regulating disease symptoms

Viral infection affects the pattern of plant miRNA expression. It has been presumed that reduction of miR171 and several other miRNAs influences viral symptoms in plants. We here experimentally demonstrate the association of osa-miR171b with rice stripe virus (RSV) symptoms in rice. Inhibition of osa-miR171b caused stunting with reduced chlorophyll content in leaves similar to viral symptoms. Overexpression of osa-miR171b by an artificial miRNA extended vegetative growth and enhanced chlorophyll accumulation in leaves. Tillers were thicker, and panicles were longer with more spikelets in plants overexpressing osa-miR171b than in controls, but there were no differences in tiller numbers. Targets of osa-miR171b, OsSCL6-IIa, OsSCL6-IIb, and OsSCL6-IIc, were respectively up- and down-regulated in plants where osa-miR171b was inhibited or overexpressed. In plants overexpressing osa-miR171b, five positive regulators for heading development, Ehd1, Ehd2, Ehd3, Ehd4, and Hd3a were up-regulated, while the negative regulator Ghd7 was down-regulated. Plants overexpressing osa-miR171b were less susceptible to RSV and virus symptoms were attenuated. Taken together, the results reveal that a reduction of osa-miR171b in RSV-infected rice contributes to RSV symptoms, and provide more insight into the roles of osa-miR171b in rice.

Deep sequencing leads to the identification of eukaryotic translation initiation factor 5A as a key element in Rsv1-mediated lethal systemic hypersensitive response to Soybean mosaic virus infection in soybean

Rsv1, a single dominant resistance locus in soybean, confers extreme resistance to the majority of Soybean mosaic virus (SMV) strains, but is susceptible to the G7 strain. In Rsv1-genotype soybean, G7 infection provokes a lethal systemic hypersensitive response (LSHR), a delayed host defence response. The Rsv1-mediated LSHR signalling pathway remains largely unknown. In this study, we employed a genome-wide investigation to gain an insight into the molecular interplay between SMV G7 and Rsv1-genotype soybean. Small RNA (sRNA), degradome and transcriptome sequencing analyses were used to identify differentially expressed genes (DEGs) and microRNAs (DEMs) in response to G7 infection. A number of DEGs, DEMs and microRNA targets, and the interaction network of DEMs and their target mRNAs responsive to G7 infection, were identified. Knock-down of one of the identified DEGs, the eukaryotic translation initiation factor 5A (eIF5A), diminished the LSHR and enhanced viral accumulation, suggesting the essential role of eIF5A in the G7-induced, Rsv1-mediated LSHR signalling pathway. This work provides an in-depth genome-wide analysis of high-throughput sequencing data, and identifies multiple genes and microRNA signatures that are associated with the Rsv1-mediated LSHR.

Identification and validation of a virus-inducible ta-siRNA-generating TAS4 locus in tomato

Trans-acting small interfering RNAs (ta-siRNAs) are a class of endogenous small RNA, associated with posttranscriptional gene silencing. Their biogenesis requires an initial microRNA (miRNA)-mediated cleavage of precursor RNA. Around 20 different ta-siRNA-producing loci (TASs), whose sequences are conserved, are reported in plants. In tomato, two TAS gene families have been identified, which are found to target auxin response factor gene and bacterial spot disease resistance protein Bs4 gene. Using high-throughput computational and experimental approach, we identified a new locus-producing ta-siRNA in tomato. We have also identified the putative miRNA regulating the production of ta-siRNA from this locus. The ta-siRNAs generated from TAS4 were up-regulated upon infection with a DNA virus. The potential targets of ta-siRNAs were predicted to be variety of proteins including MYB transcription factors and cell cycle regulators for some of the ta-siRNAs produced.

The function of small RNAs in plant biotic stress response

Small RNAs (sRNAs) play essential roles in plants upon biotic stress. Plants utilize RNA silencing machinery to facilitate pathogen-associated molecular pattern-triggered immunity and effector-triggered immunity to defend against pathogen attack or to facilitate defense against insect herbivores. Pathogens, on the other hand, are also able to generate effectors and sRNAs to counter the host immune response. The arms race between plants and pathogens/insect herbivores has triggered the evolution of sRNAs, RNA silencing machinery and pathogen effectors. A great number of studies have been performed to investigate the roles of sRNAs in plant defense, bringing in the opportunity to utilize sRNAs in plant protection. Transgenic plants with pathogen-derived resistance ability or transgenerational defense have been generated, which show promising potential as solutions for pathogen/insect herbivore problems in the field. Here we summarize the recent progress on the function of sRNAs in response to biotic stress, mainly in plant-pathogen/insect herbivore interaction, and the application of sRNAs in disease and insect herbivore control.

Genome-Wide Investigation Using sRNA-Seq, Degradome-Seq and Transcriptome-Seq Reveals Regulatory Networks of microRNAs and Their Target Genes in Soybean during Soybean mosaic virus Infection

MicroRNAs (miRNAs) play key roles in a variety of cellular processes through regulation of their target gene expression. Accumulated experimental evidence has demonstrated that infections by viruses are associated with the altered expression profile of miRNAs and their mRNA targets in the host. However, the regulatory network of miRNA-mRNA interactions during viral infection remains largely unknown. In this study, we performed small RNA (sRNA)-seq, degradome-seq and as well as a genome-wide transcriptome analysis to profile the global gene and miRNA expression in soybean following infections by three different Soybean mosaic virus (SMV) isolates, L (G2 strain), LRB (G2 strain) and G7 (G7 strain). sRNA-seq analyses revealed a total of 253 soybean miRNAs with a two-fold or greater change in abundance compared with the mock-inoculated control. 125 transcripts were identified as the potential cleavage targets of 105 miRNAs and validated by degradome-seq analyses. Genome-wide transcriptome analysis showed that total 2679 genes are differentially expressed in response to SMV infection including 71 genes predicted as involved in defense response. Finally, complex miRNA-mRNA regulatory networks were derived using the RNAseq, small RNAseq and degradome data. This work represents a comprehensive, global approach to examining virus-host interactions. Genes responsive to SMV infection are identified as are their potential miRNA regulators. Additionally, regulatory changes of the miRNAs themselves are described and the regulatory relationships were supported with degradome data. Taken together these data provide new insights into molecular SMV-soybean interactions and offer candidate miRNAs and their targets for further elucidation of the SMV infection process.

Genome-wide identification of abiotic stress-regulated and novel microRNAs in mulberry leaf

As the most important food plant for sericultural industry, mulberry trees have to suffer from a wide range of abiotic and biotic stresses, such as drought and high salinity. MicroRNAs (miRNAs) have been proved to play important roles in abiotic stresses regulation in many plants. However, there are seldom reports on the miRNAs expression profiles upon abiotic challenges in mulberry. In this study, three small RNA libraries from mulberry leaf tissue with or without drought or salt treatment were constructed and deep sequenced. Total of 48 conserved miRNAs (including miRNA*) and 162 novel miRNAs were identified (processing precision value>0.1). A total of 270 and 1963 target genes were predicted for conserved miRNAs and novel miRNAs, respectively. 13 differentially expressed miRNAs were detected under drought or salt stresses by deep sequencing and qRT-PCR. 5' RLM-RACE validated Morus 013341 to be the target gene of miR-395a. Our results provided initial clue to further study molecular mechanism on abiotic stresses regulation in mulberry.