Premeiotic 24-nt phasiRNAs are present in the Zea genus and unique in biogenesis mechanism and molecular function.

Reproductive phasiRNAs (phased, small interfering RNAs) are broadly present in angiosperms and play crucial roles in sustaining male fertility. While the premeiotic 21-nt (nucleotides) phasiRNAs and meiotic 24-nt phasiRNA pathways have been extensively studied in maize (Zea mays) and rice (Oryza sativa), a third putative category of reproductive phasiRNAs-named premeiotic 24-nt phasiRNAs-have recently been reported in barley (Hordeum vulgare) and wheat (Triticum aestivum). To determine whether premeiotic 24-nt phasiRNAs are also present in maize and related species and begin to characterize their biogenesis and function, we performed a comparative transcriptome and degradome analysis of premeiotic and meiotic anthers from five maize inbred lines and three teosinte species/subspecies. Our data indicate that a substantial subset of the 24-nt phasiRNA loci in maize and teosinte are already highly expressed at the premeiotic phase. The premeiotic 24-nt phasiRNAs are similar to meiotic 24-nt phasiRNAs in genomic origin and dependence on DCL5 (Dicer-like 5) for biogenesis, however, premeiotic 24-nt phasiRNAs are unique in that they are likely i) not triggered by microRNAs, ii) not loaded by AGO18 proteins, and iii) not capable of mediating PHAS precursor cleavage. In addition, we also observed a group of premeiotic 24-nt phasiRNAs in rice using previously published data. Together, our results indicate that the premeiotic 24-nt phasiRNAs constitute a unique class of reproductive phasiRNAs and are present more broadly in the grass family (Poaceae) than previously known.

Coexpression network and trans-activation analyses of maize reproductive phasiRNA loci.

The anther-enriched phased, small interfering RNAs (phasiRNAs) play vital roles in sustaining male fertility in grass species. Their long non-coding precursors are synthesized by RNA polymerase II and are likely regulated by transcription factors (TFs). A few putative transcriptional regulators of the 21- or 24-nucleotide phasiRNA loci (referred to as 21- or 24-PHAS loci) have been identified in maize (Zea mays), but whether any of the individual TFs or TF combinations suffice to activate any PHAS locus is unclear. Here, we identified the temporal gene coexpression networks (modules) associated with maize anther development, including two modules highly enriched for the 21- or 24-PHAS loci. Comparisons of these coexpression modules and gene sets dysregulated in several reported male sterile TF mutants provided insights into TF timing with regard to phasiRNA biogenesis, including antagonistic roles for OUTER CELL LAYER4 and MALE STERILE23. Trans-activation assays in maize protoplasts of individual TFs using bulk-protoplast RNA-sequencing showed that two of the TFs coexpressed with 21-PHAS loci could activate several 21-nucleotide phasiRNA pathway genes but not transcription of 21-PHAS loci. Screens for combinatorial activities of these TFs and, separately, the recently reported putative transcriptional regulators of 24-PHAS loci using single-cell (protoplast) RNA-sequencing, did not detect reproducible activation of either 21-PHAS or 24-PHAS loci. Collectively, our results suggest that the endogenous transcriptional machineries and/or chromatin states in the anthers are necessary to activate reproductive PHAS loci.

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.

Dynamically expressed small RNAs, substantially driven by genomic structural variants, contribute to transcriptomic changes during tomato domestication.

Tomato has undergone extensive selections during domestication. Recent progress has shown that genomic structural variants (SVs) have contributed to gene expression dynamics during tomato domestication, resulting in changes of important traits. Here, we performed comprehensive analyses of small RNAs (sRNAs) from nine representative tomato accessions. We demonstrate that SVs substantially contribute to the dynamic expression of the three major classes of plant sRNAs: microRNAs (miRNAs), phased secondary short interfering RNAs (phasiRNAs), and 24-nucleotide heterochromatic siRNAs (hc-siRNAs). Changes in the abundance of phasiRNAs and 24-nucleotide hc-siRNAs likely contribute to the alteration of mRNA gene expression in cis during tomato domestication, particularly for genes associated with biotic and abiotic stress tolerance. We also observe that miRNA expression dynamics are associated with imprecise processing, alternative miRNA-miRNA* selections, and SVs. SVs mainly affect the expression of less-conserved miRNAs that do not have established regulatory functions or low abundant members in highly expressed miRNA families. Our data highlight different selection pressures on miRNAs compared to phasiRNAs and 24-nucleotide hc-siRNAs. Our findings provide insights into plant sRNA evolution as well as SV-based gene regulation during crop domestication. Furthermore, our dataset provides a rich resource for mining the sRNA regulatory network in tomato.

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.

CRISPR/Cas9-targeted mutagenesis of TaDCL4, TaDCL5 and TaRDR6 induces male sterility in common wheat.

Phased, small interfering RNAs (phasiRNAs) are important for plant anther development, especially for male sterility. PhasiRNA biogenesis is dependent on genes like RNA polymerase 6 (RDR6), DICER-LIKE 4 (DCL4), or DCL5 to produce 21- or 24 nucleotide (nt) double-strand small RNAs. Here, we generated mutants of DCL4, DCL5 and RDR6 using CRISPR/Cas9 system and studied their effects on plant reproductive development and phasiRNA production in wheat. We found that RDR6 mutation caused sever consequence throughout plant development starting from seed germination and the dcl4 mutants grew weaker with thorough male sterility, while dcl5 plants developed normally but exhibited male sterility. Correspondingly, DCL4 and DCL5, respectively, specified 21- and 24-nt phasiRNA biogenesis, while RDR6 contributed to both. Also, the three key genes evolved differently in wheat, with TaDCL5-A/B becoming non-functioning and TaRDR6-A being lost after polyploidization. Furthermore, we found that PHAS genes (phasiRNA precursors) identified via phasiRNAs diverged rapidly among sub-genomes of polyploid wheat. Despite no similarity being found among phasiRNAs of grasses, their targets were enriched for similar biological functions. In light of the important roles of phasiRNA pathways in gametophyte development, genetic dissection of the function of key genes may help generate male sterile lines suitable for hybrid wheat breeding.

Small noncoding RNA discovery and profiling with sRNAtools based on high-throughput sequencing.

Small noncoding RNAs (sRNA/sncRNAs) are generated from different genomic loci and play important roles in biological processes, such as cell proliferation and the regulation of gene expression. Next-generation sequencing (NGS) has provided an unprecedented opportunity to discover and quantify diverse kinds of sncRNA, such as tRFs (tRNA-derived small RNA fragments), phasiRNAs (phased, secondary, small-interfering RNAs), Piwi-interacting RNA (piRNAs) and plant-specific 24-nt short interfering RNAs (siRNAs). However, currently available web-based tools do not provide approaches to comprehensively analyze all of these diverse sncRNAs. This study presents a novel integrated platform, sRNAtools (https://bioinformatics.caf.ac.cn/sRNAtools), that can be used in conjunction with high-throughput sequencing to identify and functionally annotate sncRNAs, including profiling microRNAss, piRNAs, tRNAs, small nuclear RNAs, small nucleolar RNAs and rRNAs and discovering isomiRs, tRFs, phasiRNAs and plant-specific 24-nt siRNAs for up to 21 model organisms. Different modules, including single case, batch case, group case and target case, are developed to provide users with flexible ways of studying sncRNA. In addition, sRNAtools supports different ways of uploading small RNA sequencing data in a very interactive queue system, while local versions based on the program package/Docker/virtureBox are also available. We believe that sRNAtools will greatly benefit the scientific community as an integrated tool for studying sncRNAs.

Conserved and non-conserved triggers of 24-nucleotide reproductive phasiRNAs in eudicots.

Small RNAs play important roles in plant growth and development by modulating expression of genes and transposons. In many flowering plant species, male reproductive organs, the anthers, produce abundant phased small interfering RNAs (phasiRNAs). Two classes of reproductive phasiRNAs are generally known, mostly from monocots: (i) pre-meiotic 21-nucleotide (nt) phasiRNAs triggered by miR2118 and (ii) meiotic 24-nt phasiRNAs triggered by miR2275. Here, we describe conserved and non-conserved triggers of 24-nt phasiRNAs in several eudicots. We found that the abundant 24-nt phasiRNAs in the basal eudicot columbine (Aquilegia coerulea) are produced by the canonical trigger miR2275, as well as by other non-canonical triggers, miR482/2118 and miR14051. These triggering microRNAs (miRNAs) are localized in microspore mother cells and tapetal cells of meiotic and post-meiotic stage anthers. Furthermore, we identified a lineage-specific trigger (miR11308) of 24-nt phasiRNAs and an expanded number of 24-PHAS loci in wild strawberry (Fragaria vesca). We validated the presence of the miR2275-derived 24-nt phasiRNA pathway in rose (Rosa chinensis). Finally, we evaluated all eudicots that have been validated for the presence of 24-nt phasiRNAs as possible model systems in which to study the biogenesis and function of 24-nt phasiRNAs. We conclude that columbine (Aquilegia coerulea) would be a strong model because of its extensive number of 24-PHAS loci and its diversity of trigger miRNAs.

Heightened miR6024-NLR interactions facilitate necrotrophic pathogenesis in tomato.

KEY MESSAGE: miR6024 acts as a negative regulator of R genes, hence of Tomato plant immunity, and facilitates disease by the necrotrophic pathogen A. solani. Plant resistance genes or Nucleotide-binding leucine-rich repeat (NLR) genes, integral components of plant disease stress-signaling are targeted by variable groups of miRNAs. However, the significance of miRNA-mediated regulation of NLRs during a pathogen stress response, specifically for necrotrophic fungus, is poorly understood. A thorough examination of Tomato NLRs and miRNAs could map substantial interactions of which half the annotated NLRs were targets of Solanaceae-specific and conserved miRNAs, at the NB subdomain. The Solanaceae-specific miR6024 and its NLR targets analysed in different phytopathogenic stresses revealed differential and mutually antagonistic regulation. Interestingly, miR6024-targeted cleavage of a target NLR also triggered the generation of secondary phased siRNAs which could potentially amplify the defense signal. RNA-seq analysis of leaf tissues from miR6024 overexpressing Tomato plants evidenced a perturbation in the defense transcriptome with the transgenics showing unwarranted immune response-related genes' expression with or without infection with necrotrophic Alternaria solani, though no adverse effect could be observed in the growth and development of the transgenic plants. Transgenic plants exhibited constitutive downregulation of the target NLRs, aggravated disease phenotype with an enhanced lesion, greater ROS generation and hypersusceptibility to A. solani infection, thus establishing that miR6024 negatively impacts plant immune response during necrotrophic pathogenesis. Limited knowledge about the outcome of NLR-miRNA interaction during necrotrophic pathogenesis is a hindrance to the deployment of miRNAs in crop improvement programs. With the elucidation of the necrotrophic disease-synergistic role played by miR6024, it becomes a potent candidate for biotechnological manipulation for the rapid development of pathogen-tolerant solanaceous plants.

Temperature-sensitive male sterility in rice determined by the roles of AGO1d in reproductive phasiRNA biogenesis and function.

Phased secondary siRNAs (phasiRNAs) are broadly present in the reproductive tissues of flowering plants, with spatial-temporal specificity. However, the ARGONAUTE (AGO) proteins associated with phasiRNAs and their miRNA triggers remain elusive. Here, through histological and high-throughput sequencing analyses, we show that rice AGO1d, which is specifically expressed in anther wall cells before and during meiosis, associates with both miR2118 and miR2275 to mediate phasiRNA biogenesis. AGO1d preferentially binds to miR2118-triggered 21-nucleotide (nt) phasiRNAs with a 5'-terminal uridine, suggesting a dual role in phasiRNA biogenesis and function. Depletion of AGO1d causes a reduction of 21- and 24-nt phasiRNAs and temperature-sensitive male sterility. At lower temperatures, anthers of the ago1d mutant predominantly show excessive tapetal cells with little starch accumulation during pollen formation, possibly caused by the dysregulation of cell metabolism. These results uncover an essential role of AGO1d in rice anther development at lower temperatures and demonstrate coordinative roles of AGO proteins during reproductive phasiRNA biogenesis and function.

The red flesh of kiwifruit is differentially controlled by specific activation-repression systems.

Anthocyanins are visual cues for pollination and seed dispersal. Fruit containing anthocyanins also appeals to consumers due to its appearance and health benefits. In kiwifruit (Actinidia spp.) studies have identified at least two MYB activators of anthocyanin, but their functions in fruit and the mechanisms by which they act are not fully understood. Here, transcriptome and small RNA high-throughput sequencing were used to comprehensively identify contributors to anthocyanin accumulation in kiwifruit. Stable overexpression in vines showed that both 35S::MYB10 and MYB110 can upregulate anthocyanin biosynthesis in Actinidia chinensis fruit, and that MYB10 overexpression resulted in anthocyanin accumulation which was limited to the inner pericarp, suggesting that repressive mechanisms underlie anthocyanin biosynthesis in this species. Furthermore, motifs in the C-terminal region of MYB10/110 were shown to be responsible for the strength of activation of the anthocyanic response. Transient assays showed that both MYB10 and MYB110 were not directly cleaved by miRNAs, but that miR828 and its phased small RNA AcTAS4-D4(-) efficiently targeted MYB110. Other miRNAs were identified, which were differentially expressed between the inner and outer pericarp, and cleavage of SPL13, ARF16, SCL6 and F-box1, all of which are repressors of MYB10, was observed. We conclude that it is the differential expression and subsequent repression of MYB activators that is responsible for variation in anthocyanin accumulation in kiwifruit species.

24-nt phasiRNAs move from tapetal to meiotic cells in maize anthers.

In maize, 24-nt phased, secondary small interfering RNAs (phasiRNAs) are abundant in meiotic stage anthers, but their distribution and functions are not precisely known. Using laser capture microdissection, we analyzed tapetal cells, meiocytes and other somatic cells at several stages of anther development to establish the timing of 24-PHAS precursor transcripts and the 24-nt phasiRNA products. By integrating RNA and small RNA profiling plus single-molecule and small RNA FISH (smFISH or sRNA-FISH) spatial detection, we demonstrate that the tapetum is the primary site of 24-PHAS precursor and Dcl5 transcripts and the resulting 24-nt phasiRNAs. Interestingly, 24-nt phasiRNAs accumulate in all cell types, with the highest levels in meiocytes, followed by tapetum. Our data support the conclusion that 24-nt phasiRNAs are mobile from tapetum to meiocytes and to other somatic cells. We discuss possible roles for 24-nt phasiRNAs in anther cell types.

TarDB: an online database for plant miRNA targets and miRNA-triggered phased siRNAs.

BACKGROUND: In plants, microRNAs (miRNAs) are pivotal regulators of plant development and stress responses. Different computational tools and web servers have been developed for plant miRNA target prediction; however, in silico prediction normally contains false positive results. In addition, many plant miRNA target prediction servers lack information for miRNA-triggered phased small interfering RNAs (phasiRNAs). Creating a comprehensive and relatively high-confidence plant miRNA target database is much needed. RESULTS: Here, we report TarDB, an online database that collects three categories of relatively high-confidence plant miRNA targets: (i) cross-species conserved miRNA targets; (ii) degradome/PARE (Parallel Analysis of RNA Ends) sequencing supported miRNA targets; (iii) miRNA-triggered phasiRNA loci. TarDB provides a user-friendly interface that enables users to easily search, browse and retrieve miRNA targets and miRNA initiated phasiRNAs in a broad variety of plants. TarDB has a comprehensive collection of reliable plant miRNA targets containing previously unreported miRNA targets and miRNA-triggered phasiRNAs even in the well-studied model species. Most of these novel miRNA targets are relevant to lineage-specific or species-specific miRNAs. TarDB data is freely available at http://www.biosequencing.cn/TarDB . CONCLUSIONS: In summary, TarDB serves as a useful web resource for exploring relatively high-confidence miRNA targets and miRNA-triggered phasiRNAs in plants.

fIdentification of B. napus small RNAs responsive to infection by a necrotrophic pathogen.

BACKGROUND: Small RNAs are short non-coding RNAs that are key gene regulators controlling various biological processes in eukaryotes. Plants may regulate discrete sets of sRNAs in response to pathogen attack. Sclerotinia sclerotiorum is an economically important pathogen affecting hundreds of plant species, including the economically important oilseed B. napus. However, there are limited studies on how regulation of sRNAs occurs in the S. sclerotiorum and B. napus pathosystem. RESULTS: We identified different classes of sRNAs from B. napus using high throughput sequencing of replicated mock and infected samples at 24 h post-inoculation (HPI). Overall, 3999 sRNA loci were highly expressed, of which 730 were significantly upregulated during infection. These 730 up-regulated sRNAs targeted 64 genes, including disease resistance proteins and transcriptional regulators. A total of 73 conserved miRNA families were identified in our dataset. Degradome sequencing identified 2124 cleaved mRNA products from these miRNAs from combined mock and infected samples. Among these, 50 genes were specific to infection. Altogether, 20 conserved miRNAs were differentially expressed and 8 transcripts were cleaved by the differentially expressed miRNAs miR159, miR5139, and miR390, suggesting they may have a role in the S. sclerotiorum response. A miR1885-triggered disease resistance gene-derived secondary sRNA locus was also identified and verified with degradome sequencing. We also found further evidence for silencing of a plant immunity related ethylene response factor gene by a novel sRNA using 5'-RACE and RT-qPCR. CONCLUSIONS: The findings in this study expand the framework for understanding the molecular mechanisms of the S. sclerotiorum and B. napus pathosystem at the sRNA level.

Evolution and diversification of reproductive phased small interfering RNAs in Oryza species.

In grasses, two types of phased, small interfering RNAs (phasiRNAs) are expressed largely in young, developing anthers. They are 21 or 24 nucleotides (nt) in length and are triggered by miR2118 or miR2275, respectively. However, most of their functions and activities are not fully understood. We performed comparative genomic analysis of their source loci (PHAS) in five Oryza genomes and combined this with analysis of high-throughput sRNA and degradome datasets. In total, we identified 8216 21-PHAS and 626 24-PHAS loci. Local tandem and segmental duplications mainly contributed to the expansion and supercluster distribution of the 21-PHAS loci. Despite their relatively conserved genomic positions, PHAS sequences diverged rapidly, except for the miR2118/2275 target sites, which were under strong selection for conservation. We found that 21-nt phasiRNAs with a 5'-terminal uridine (U) demonstrated cis-cleavage at PHAS precursors, and these cis-acting sites were also variable among close species. miR2118 could trigger phasiRNA production from its own antisense transcript and the derived phasiRNAs might reversibly regulate miR2118 precursors. We hypothesised that successful initiation of phasiRNA biogenesis is conservatively maintained, while phasiRNA products diverged quickly and are not individually conserved. In particular, phasiRNA production is under the control of multiple reciprocal regulation mechanisms.

CHH DNA methylation increases at 24-PHAS loci depend on 24-nt phased small interfering RNAs in maize meiotic anthers.

Plant phased small interfering RNAs (phasiRNAs) contribute to robust male fertility; however, specific functions remain undefined. In maize (Zea mays), male sterile23 (ms23), necessary for both 24-nt phasiRNA precursor (24-PHAS) loci and Dicer-like5 (Dcl5) expression, and dcl5-1 mutants unable to slice PHAS transcripts lack nearly all 24-nt phasiRNAs. Based on sequence capture bisulfite-sequencing, we find that CHH DNA methylation of most 24-PHAS loci is increased in meiotic anthers of control plants but not in the ms23 and dcl5 mutants. Because dcl5-1 anthers express PHAS precursors, we conclude that the 24-nt phasiRNAs, rather than just activation of PHAS transcription, are required for targeting increased CHH methylation at these loci. Although PHAS precursors are processed into multiple 24-nt phasiRNA products, there is substantial differential product accumulation. Abundant 24-nt phasiRNA positions corresponded to high CHH methylation within individual loci, reinforcing the conclusion that 24-nt phasiRNAs contribute to increased CHH methylation in cis.

Small RNA profiling from meiotic and post-meiotic anthers reveals prospective miRNA-target modules for engineering male fertility in sorghum.

Understanding male gametophyte development is essential to augment hybrid production in sorghum. Although small RNAs are known to critically influence anther/pollen development, their roles in sorghum reproduction have not been deciphered yet. Here, we report small RNA profiling and high-confidence annotation of microRNAs (miRNAs) from meiotic and post-meiotic anthers in sorghum. We identified 262 miRNAs (82 known and 180 novel), out of which 58 (35 known and 23 novel) exhibited differential expression between two stages. Out of 35 differentially expressed known miRNAs, 13 are known to regulate anther/pollen development in other plant species. We also demonstrated conserved spatiotemporal patterns of 21- and 24-nt phasiRNAs and their respective triggers, miR2118 and miR2275, in sorghum anthers as evidenced in other monocots. miRNA target identification yielded 5622 modules, of which 46 modules comprising 16 known and 8 novel miRNA families with 38 target genes are prospective candidates for engineering male fertility in grasses.

Differential Response of Grapevine to Infection with 'Candidatus Phytoplasma solani' in Early and Late Growing Season through Complex Regulation of mRNA and Small RNA Transcriptomes.

Bois noir is the most widespread phytoplasma grapevine disease in Europe. It is associated with 'Candidatus Phytoplasma solani', but molecular interactions between the causal pathogen and its host plant are not well understood. In this work, we combined the analysis of high-throughput RNA-Seq and sRNA-Seq data with interaction network analysis for finding new cross-talks among pathways involved in infection of grapevine cv. Zweigelt with 'Ca. P. solani' in early and late growing seasons. While the early growing season was very dynamic at the transcriptional level in asymptomatic grapevines, the regulation at the level of small RNAs was more pronounced later in the season when symptoms developed in infected grapevines. Most differentially expressed small RNAs were associated with biotic stress. Our study also exposes the less-studied role of hormones in disease development and shows that hormonal balance was already perturbed before symptoms development in infected grapevines. Analysis at the level of communities of genes and mRNA-microRNA interaction networks revealed several new genes (e.g., expansins and cryptdin) that have not been associated with phytoplasma pathogenicity previously. These novel actors may present a new reference framework for research and diagnostics of phytoplasma diseases of grapevine.

Identification of tRFs and phasiRNAs in tomato (Solanum lycopersicum) and their responses to exogenous abscisic acid.

BACKGROUND: The non-coding small RNA tRFs (tRNA-derived fragments) and phasiRNAs (plant-specific) exert important roles in plant growth, development and stress resistances. However, whether the tRFs and phasiRNAs respond to the plant important stress hormone abscisic acid (ABA) remain enigma. RESULTS: Here, the RNA-sequencing was implemented to decipher the landscape of tRFs and phasiRNAs in tomato (Solanum lycopersicum) leaves and their responses when foliar spraying exogenous ABA after 24 h. In total, 733 tRFs and 137 phasiRNAs were detected. The tRFs were mainly derived from the tRNAAla transporting alanine, which tended to be cleaved at the 5'terminal guanine site and D loop uracil site to produce tRFAla with length of 20 nt. Most of phasiRNAs originated from NBS-LRR resistance genes. Expression analysis revealed that 156 tRFs and 68 phasiRNAs expressed differentially, respectively. Generally, exogenous ABA mainly inhibited the expression of tRFs and phasiRNAs. Furthermore, integrating analysis of target gene prediction and transcriptome data presented that ABA significantly downregulated the abundance of phsaiRNAs associated with biological and abiotic resistances. Correspondingly, their target genes such as AP2/ERF, WRKY and NBS-LRR, STK and RLK, were mainly up-regulated. CONCLUSIONS: Combined with the previous analysis of ABA-response miRNAs, it was speculated that ABA can improve the plant resistances to various stresses by regulating the expression and interaction of small RNAs (such as miRNAs, tRFs, phasiRNAs) and their target genes. This study enriches the plant tRFs and phasiRNAs, providing a vital basis for further investigating ABA response-tRFs and phasiRNAs and their functions in biotic and abiotic stresses.

Suppression of NB-LRR genes by miRNAs promotes nitrogen-fixing nodule development in Medicago truncatula.

Plant genomes contain two major classes of innate immune receptors to recognize different pathogens. The pattern recognition receptors perceive conserved pathogen-associated molecular patterns and the resistance genes with nucleotide-binding (NB) and leucine-rich repeat (LRR) domains recognize specific pathogen effectors. The precise regulation of resistance genes is important since the unregulated expression of NB-LRR genes can inhibit growth and may result in autoimmunity in the absence of pathogen infection. It was shown that a subset of miRNAs could target NB-LRR genes and act as an important regulator of plant immunity in the absence of pathogens. Plants not only interact with pathogens, but they can also establish symbiotic interactions with microbes. Nitrogen-fixing symbiotic interaction and nodule formation of legumes may also require the suppression of host defence to prevent immune responses. We found that upon symbiotic interactions, miRNAs repressing NB-LRR expression are upregulated in the developing nodules of Medicago truncatula. Furthermore, we show that the suppression of the activity of the NB-LRR genes targeted by these miRNAs is important during nodule development. Our results suggest that the downregulation of NB-LRR resistance genes in the developing nodule produces a suitable niche that facilitates bacterial colonization and the development of an N-fixing nodule.