Post-transcriptional small RNA pathways in plants: mechanisms and regulations

Molecular Manipulation of the miR160/AUXIN RESPONSE FACTOR Expression Module Impacts Root Development in Arabidopsis thaliana

In Arabidopsis thaliana (Arabidopsis), microRNA160 (miR160) regulates the expression of AUXIN RESPONSE FACTOR10 (ARF10), ARF16 and ARF17 throughout development, including the development of the root system. We have previously shown that in addition to DOUBLE-STRANDED RNA BINDING1 (DRB1), DRB2 is also involved in controlling the rate of production of specific miRNA cohorts in the tissues where DRB2 is expressed in wild-type Arabidopsis plants. In this study, a miR160 overexpression transgene (MIR160B) and miR160-resistant transgene versions of ARF10 and ARF16 (mARF10 and mARF16) were introduced into wild-type Arabidopsis plants and the drb1 and drb2 single mutants to determine the degree of requirement of DRB2 to regulate the miR160 expression module as part of root development. Via this molecular modification approach, we show that in addition to DRB1, DRB2 is required to regulate the level of miR160 production from its precursor transcripts in Arabidopsis roots. Furthermore, we go on to correlate the altered abundance of miR160 or its ARF10, ARF16 and ARF17 target genes in the generated series of transformant lines with the enhanced development of the root system displayed by these plant lines. More specifically, promotion of primary root elongation likely stemmed from enhancement of miR160-directed ARF17 expression repression, while the promotion of lateral and adventitious root formation was the result of an elevated degree of miR160-directed regulation of ARF17 expression, and to a lesser degree, ARF10 and ARF16 expression. Taken together, the results presented in this study identify the requirement of the functional interplay between DRB1 and DRB2 to tightly control the rate of miR160 production, to in turn ensure the appropriate degree of miR160-directed ARF10, ARF16 and ARF17 gene expression regulation as part of normal root system development in Arabidopsis.

Comparative transcriptome analysis of resistant and susceptible watermelon genotypes reveals the role of RNAi, callose, proteinase, and cell wall in squash vein yellowing virus resistance

Watermelon (Citrullus lanatus) is the third largest fruit crop in the world in term of production. However, it is susceptible to several viruses. Watermelon vine decline (WVD), caused by whitefly-transmitted squash vein yellowing virus (SqVYV), is a disease that has caused over $60 million in losses in the US and continues to occur regularly in southeastern states. Understanding the molecular mechanisms underlying resistance to SqVYV is important for effective disease management. A time-course transcriptomic analysis was conducted on resistant (392291-VDR) and susceptible (Crimson Sweet) watermelon genotypes inoculated with SqVYV. Significantly higher levels of SqVYV were observed over time in the susceptible compared to the resistant genotype. The plasmodesmata callose binding protein (PDCB) gene, which is responsible for increased callose deposition in the plasmodesmata, was more highly expressed in the resistant genotype than in the susceptible genotype before and after inoculation, suggesting the inhibition of cell-to-cell movement of SqVYV. The potential role of the RNA interference (RNAi) pathway was observed in the resistant genotype based on differential expression of eukaryotic initiation factor (eIF), translin, DICER, ribosome inactivating proteins, RNA-dependent RNA polymerase (RDR), and Argonaute (AGO) genes after inoculation. The significant differential expression of hormone-related genes, including those involved in the ethylene, jasmonic acid, auxin, cytokinin, gibberellin, and salicylic acid signaling pathways, was observed, emphasizing their regulatory roles in the defense response. Genes regulating pectin metabolism, cellulose synthesis, cell growth and development, xenobiotic metabolism, and lignin biosynthesis were overexpressed in the susceptible genotype, suggesting that alterations in cell wall integrity and growth processes result in disease symptom development. These findings will be helpful for further functional studies and the development of SqVYV-resistant watermelon cultivars.

ZmDST44 Gene Is a Positive Regulator in Plant Drought Stress Tolerance

Improving drought tolerance in plants is essential for increasing crop yields under water-limited conditions. In this study, we investigated the functional role of the maize gene ZmDST44, which is targeted by the miRNA ZmmiR139. Our results indicate that ZmmiR139 regulates ZmDST44 by cleaving its mRNA, as confirmed by inverse expression patterns and 5'-RACE analysis. Overexpression of ZmDST44 in Arabidopsis, rice, and maize resulted in significant enhancements in drought tolerance. Transgenic plants exhibited reduced malondialdehyde (MDA) levels, increased proline accumulation, and upregulation of drought-responsive genes compared to wild-type plants. Transgenic Arabidopsis and rice showed improved drought resistance and higher post-drought recovery rates, and transgenic maize displayed lower sensitivity to drought stress. These findings suggest that ZmDST44 acts as a positive regulator of drought tolerance across different plant species and holds promise for developing drought-resistant crops through genetic engineering.

Citrus psorosis virus 24K protein inhibits the processing of miRNA precursors by interacting with components of the biogenesis machinery

Sweet orange (Citrus sinensis) is one of the most important fruit crops worldwide. Virus infections in this crop can interfere with cellular processes, causing dramatic economic losses. By performing RT-qPCR analyses, we demonstrated that citrus psorosis virus (CPsV)-infected orange plants exhibited higher levels of unprocessed microRNA (miRNA) precursors than healthy plants. This result correlated with the reported reduction of mature miRNAs species. The protein 24K, the CPsV suppressor of RNA silencing (VSR), interacts with miRNA precursors in vivo. Thus, this protein becomes a candidate responsible for the increased accumulation of unprocessed miRNAs. We analyzed 24K RNA-binding and protein-protein interaction domains and described patterns of its subcellular localization. We also showed that 24K colocalizes within nuclear D-bodies with the miRNA biogenesis proteins DICER-LIKE 1 (DCL1), HYPONASTIC LEAVES 1 (HYL1), and SERRATE (SE). According to the results of bimolecular fluorescence complementation and co-immunoprecipitation assays, the 24K protein interacts with HYL1 and SE. Thus, 24K may inhibit miRNA processing in CPsV-infected citrus plants by direct interaction with the miRNA processing complex. This work contributes to the understanding of how a virus can alter the regulatory mechanisms of the host, particularly miRNA biogenesis and function.IMPORTANCESweet oranges can suffer from disease symptoms induced by virus infections, thus resulting in drastic economic losses. In sweet orange plants, CPsV alters the accumulation of some precursors from the regulatory molecules called miRNAs. This alteration leads to a decreased level of mature miRNA species. This misregulation may be due to a direct association of one of the viral proteins (24K) with miRNA precursors. On the other hand, 24K may act with components of the cell miRNA processing machinery through a series of predicted RNA-binding and protein-protein interaction domains.

The plant siRNA landscape

Whereas micro (mi)RNAs are considered the clean, noble side of the small RNA world, small interfering (si)RNAs are often seen as a noisy set of molecules whose barbarian acronyms reflect a large diversity of often elusive origins and functions. Twenty-five years after their discovery in plants, however, new classes of siRNAs are still being identified, sometimes in discrete tissues or at particular developmental stages, making the plant siRNA world substantially more complex and subtle than originally anticipated. Focusing primarily on the model Arabidopsis, we review here the plant siRNA landscape, including transposable elements (TE)-derived siRNAs, a vast array of non-TE-derived endogenous siRNAs, as well as exogenous siRNAs produced in response to invading nucleic acids such as viruses or transgenes. We primarily emphasize the extraordinary sophistication and diversity of their biogenesis and, secondarily, the variety of their known or presumed functions, including via non-cell autonomous activities, in the sporophyte, gametophyte, and shortly after fertilization.

Identification of cell-type specific alternative transcripts in the multicellular alga Volvox carteri

BACKGROUND

Cell type specialization is a hallmark of complex multicellular organisms and is usually established through implementation of cell-type-specific gene expression programs. The multicellular green alga Volvox carteri has just two cell types, germ and soma, that have previously been shown to have very different transcriptome compositions which match their specialized roles. Here we interrogated another potential mechanism for differentiation in V. carteri, cell type specific alternative transcript isoforms (CTSAI).

METHODS

We used pre-existing predictions of alternative transcripts and de novo transcript assembly with HISAT2 and Ballgown software to compile a list of loci with two or more transcript isoforms, identified a small subset that were candidates for CTSAI, and manually curated this subset of genes to remove false positives. We experimentally verified three candidates using semi-quantitative RT-PCR to assess relative isoform abundance in each cell type.

RESULTS

Of the 1978 loci with two or more predicted transcript isoforms 67 of these also showed cell type isoform expression biases. After curation 15 strong candidates for CTSAI were identified, three of which were experimentally verified, and their predicted gene product functions were evaluated in light of potential cell type specific roles. A comparison of genes with predicted alternative splicing from Chlamydomonas reinhardtii, a unicellular relative of V. carteri, identified little overlap between ortholog pairs with alternative splicing in both species. Finally, we interrogated cell type expression patterns of 126 V. carteri predicted RNA binding protein (RBP) encoding genes and found 40 that showed either somatic or germ cell expression bias. These RBPs are potential mediators of CTSAI in V. carteri and suggest possible pre-adaptation for cell type specific RNA processing and a potential path for generating CTSAI in the early ancestors of metazoans and plants.

CONCLUSIONS

We predicted numerous instances of alternative transcript isoforms in Volvox, only a small subset of which showed cell type specific isoform expression bias. However, the validated examples of CTSAI supported existing hypotheses about cell type specialization in V. carteri, and also suggested new hypotheses about mechanisms of functional specialization for their gene products. Our data imply that CTSAI operates as a minor but important component of V. carteri cellular differentiation and could be used as a model for how alternative isoforms emerge and co-evolve with cell type specialization.

Genome-wide identification of DCL, AGO and RDR gene families and their associated functional regulatory element analyses in sunflower (Helianthus annuus)

RNA interference (RNAi) regulates a variety of eukaryotic gene expressions that are engaged in response to stress, growth, and the conservation of genomic stability during developmental phases. It is also intimately connected to the post-transcriptional gene silencing (PTGS) process and chromatin modification levels. The entire process of RNA interference (RNAi) pathway gene families mediates RNA silencing. The main factors of RNA silencing are the Dicer-Like (DCL), Argonaute (AGO), and RNA-dependent RNA polymerase (RDR) gene families. To the best of our knowledge, genome-wide identification of RNAi gene families like DCL, AGO, and RDR in sunflower (Helianthus annuus) has not yet been studied despite being discovered in some species. So, the goal of this study is to find the RNAi gene families like DCL, AGO, and RDR in sunflower based on bioinformatics approaches. Therefore, we accomplished an inclusive in silico investigation for genome-wide identification of RNAi pathway gene families DCL, AGO, and RDR through bioinformatics approaches such as (sequence homogeneity, phylogenetic relationship, gene structure, chromosomal localization, PPIs, GO, sub-cellular localization). In this study, we have identified five DCL (HaDCLs), fifteen AGO (HaAGOs), and ten RDR (HaRDRs) in the sunflower genome database corresponding to the RNAi genes of model plant Arabidopsis thaliana based on genome-wide analysis and a phylogenetic method. The analysis of the gene structure that contains exon-intron numbers, conserved domain, and motif composition analyses for all HaDCL, HaAGO, and HaRDR gene families indicated almost homogeneity among the same gene family. The protein-protein interaction (PPI) network analysis illustrated that there exists interconnection among identified three gene families. The analysis of the Gene Ontology (GO) enrichment showed that the detected genes directly contribute to the RNA gene-silencing and were involved in crucial pathways. It was observed that the cis-acting regulatory components connected to the identified genes were shown to be responsive to hormone, light, stress, and other functions. That was found in HaDCL, HaAGO, and HaRDR genes associated with the development and growth of plants. Finally, we are able to provide some essential information about the components of sunflower RNA silencing through our genome-wide comparison and integrated bioinformatics analysis, which open the door for further research into the functional mechanisms of the identified genes and their regulatory elements.

Genome-Wide Identification, Quantification, and Validation of Differentially Expressed miRNAs in Eggplant (Solanum melongena L.) Based on Their Response to Ralstonia solanacearum Infection

MicroRNAs (miRNAs), a type of short noncoding RNA molecule (21-23 nucleotides), mediate repressive gene regulation through RNA silencing at the posttranscriptional level and play an important role in the defense response to abiotic and biotic stresses. miRNAs of the plant system have been studied in model crops for their diverse regulatory role while less is known about their significance in other plants whose genome and transcriptome data are scarce in the database, including eggplant (Solanum melongena L.). In the present study, a next-generation sequencing platform was used for the sequencing of miRNA, and real-time quantitative PCR for miRNAs was used to validate the gene expression patterns of miRNAs in Solanum melongena plantlets infected with the bacterial wilt-causing pathogen Ralstonia solanacearum (R. solanacearum). Sequence analyses showed the presence of 375 miRNAs belonging to 29 conserved families. The miR414 is highly conserved miRNA across the plant system while miR5658 and miR5021 were found exclusively in Arabidopsis thaliana surprisingly, these miRNAs were found in eggplants too. The most abundant families were miR5658 and miR414. Ppt-miR414, hvu-miR444b, stu-miR8020, and sly miR5303 were upregulated in Pusa purple long (PPL) (susceptible) at 48 h postinfection, followed by a decline after 96 h postinfection. A similar trend was obtained in ath-miR414, stu-mir5303h, alymiR847-5p, far-miR1134, ath-miR5021, ath-miR5658, osa-miR2873c, lja-miR7530, stu-miR7997c, and gra-miR8741 but at very low levels after infection in the susceptible variety, indicating their negative role in the suppression of host immunity. On the other hand, osa-miR2873c was found to be slightly increased after 96 hpi from 48 hpi. Most of the miRNAs under study showed relatively lower expression in the resistant variety Arka Nidhi after infection than in the susceptible variety. These results shed light on a deeper regulatory role of miRNAs and their targets in regulation of the plant response to bacterial infection. The present experiment and their results suggested that the higher expression of miRNA leads to a decline in host mRNA and thus shows susceptibility.

Role of RNA silencing in plant-viroid interactions and in viroid pathogenesis

Viroids are small, single-stranded, non-protein coding and circular RNAs able to infect host plants in the absence of any helper virus. They may elicit symptoms in their hosts, but the underlying molecular pathways are only partially known. Here we address the role of post-transcriptional RNA silencing in plant-viroid-interplay, with major emphasis on the involvement of this sequence-specific RNA degradation mechanism in both plant antiviroid defence and viroid pathogenesis. This review is a tribute to the memory of Dr. Ricardo Flores, who largely contributed to elucidate this and other molecular mechanisms involved in plant-viroid interactions.

A beneficial fungal root endophyte triggers systemic RNA silencing and DNA methylation of a host reporter gene

A growing body of evidence suggests that RNA interference (RNAi) plays a pivotal role in the communication between plants and pathogenic fungi, where a bi-directional trans-kingdom RNAi is established to the advantage of either the host or the pathogen. Similar mechanisms acting during plant association with non-pathogenic symbiotic microorganisms have been elusive to this date. To determine whether root endophytes can induce systemic RNAi responses to their host plants, we designed an experimental reporter-based system consisting of the root-restricted, beneficial fungal endophyte, Fusarium solani strain K (FsK) and its host Nicotiana benthamiana. Since not all fungi encode the RNAi machinery, we first needed to validate that FsK does so, by identifying its core RNAi enzymes (2 Dicer-like genes, 2 Argonautes and 4 RNA-dependent RNA polymerases) and by showing its susceptibility to in vitro RNAi upon exogenous application of double stranded RNAs (dsRNAs). Upon establishing this, we transformed FsK with a hairpin RNA (hpRNA) construct designed to target a reporter gene in its host N. benthamiana. The hpRNA was processed by FsK RNAi machinery predominantly into 21-24-nt small RNAs that triggered RNA silencing but not DNA methylation in the fungal hyphae. Importantly, when the hpRNA-expressing FsK was used to inoculate N. benthamiana, systemic RNA silencing and DNA methylation of the host reporter gene was recorded. Our data suggest that RNAi signals can be translocated by root endophytes to their hosts and can modulate gene expression during mutualism, which may be translated to beneficial phenotypes.

Perspectives on plant virus diseases in a climate change scenario of elevated temperatures

Global food production is at risk from many abiotic and biotic stresses and can be affected by multiple stresses simultaneously. Virus diseases damage cultivated plants and decrease the marketable quality of produce. Importantly, the progression of virus diseases is strongly affected by changing climate conditions. Among climate-changing variables, temperature increase is viewed as an important factor that affects virus epidemics, which may in turn require more efficient disease management. In this review, we discuss the effect of elevated temperature on virus epidemics at both macro- and micro-climatic levels. This includes the temperature effects on virus spread both within and between host plants. Furthermore, we focus on the involvement of molecular mechanisms associated with temperature effects on plant defence to viruses in both susceptible and resistant plants. Considering various mechanisms proposed in different pathosystems, we also offer a view of the possible opportunities provided by RNA -based technologies for virus control at elevated temperatures. Recently, the potential of these technologies for topical field applications has been strengthened through a combination of genetically modified (GM)-free delivery nanoplatforms. This approach represents a promising and important climate-resilient substitute to conventional strategies for managing plant virus diseases under global warming scenarios. In this context, we discuss the knowledge gaps in the research of temperature effects on plant-virus interactions and limitations of RNA-based emerging technologies, which should be addressed in future studies.

Integrative analysis of transcriptome and miRNAome reveals molecular mechanisms regulating pericarp thickness in sweet corn during kernel development

Pericarp thickness affects the edible quality of sweet corn (Zea mays L. saccharata Sturt.). Therefore, breeding varieties with a thin pericarp is important for the quality breeding of sweet corn. However, the molecular mechanisms underlying the pericarp development remain largely unclear. We performed an integrative analysis of mRNA and miRNA sequencing to elucidate the genetic mechanism regulating pericarp thickness during kernel development (at 15 days, 19 days, and 23 days after pollination) of two sweet corn inbred lines with different pericarp thicknesses (M03, with a thinner pericarp and M08, with a thicker pericarp). A total of 2,443 and 1,409 differentially expressed genes (DEGs) were identified in M03 and M08, respectively. Our results indicate that phytohormone-mediated programmed cell death (PCD) may play a critical role in determining pericarp thickness in sweet corn. Auxin (AUX), gibberellin (GA), and brassinosteroid (BR) signal transduction may indirectly mediate PCD to regulate pericarp thickness in M03 (the thin pericarp variety). In contrast, abscisic acid (ABA), cytokinin (CK), and ethylene (ETH) signaling may be the key regulators of pericarp PCD in M08 (the thick pericarp variety). Furthermore, 110 differentially expressed microRNAs (DEMIs) and 478 differentially expressed target genes were identified. miRNA164-, miRNA167-, and miRNA156-mediated miRNA-mRNA pairs may participate in regulating pericarp thickness. The expression results of DEGs were validated by quantitative real-time PCR. These findings provide insights into the molecular mechanisms regulating pericarp thickness and propose the objective of breeding sweet corn varieties with a thin pericarp.

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

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

Integrative Analysis of the DICER-like (DCL) Genes From Peach (Prunus persica): A Critical Role in Response to Drought Stress

DICER-likes (DCLs) proteins are the core component for non-coding RNA (ncRNA) biogenesis, playing essential roles in some biological processes. The DCL family has been characterized in model plants, such as Arabidopsis, rice, and poplar. However, the evolutionary aspect and the expression mechanism under drought stress were scarce and have never been reported and characterized in one of the most important worldwide cultivated fruit trees, peach (Prunus persica). Eight DCLs genes in the Prunus persica genome were detected, in addition to 51 DCLs in the other seven Rosaceae genomes. The phylogenetic analysis with Arabidopsis thaliana and RTL1 gene as outgroups suggested that DCL members are divided into four clades: DCL1, DCL2, DCL3, and DCL4 with several gene gain/loss events of DCL gene copies through the evolutionary tract of the Rosacea family. The number of homologous DCL copies within each clade, along with the chromosomal location indicated gene duplication event of the DCL2 gene occurred once for the subfamily Amygdaloideae and twice for Pyrus communis and Prunus dulics and trice for the P. persica on Chromosome number 7 genes. Another duplication event was found for the DCL3 gene that occurred once for all the eight Rosaceae species with no match in A. thaliana. The DCL genetic similarity and activity was evaluated using BLASTp and previously published RNA-seq data among different tissues and over different time points of peach trees exposed to drought conditions. Finally, the expression pattern of PrupeDCLs in response to drought stress was identified, and two of these members, Prupe.7G047900 and Prupe.6G363600, were found as main candidate genes for response to drought stress. Our data presented here provide useful information for a better understanding of the molecular evolution of DCL genes in Rosaceae genomes, and the function of DCLs in P. persica.

A plant tethering system for the functional study of protein-RNA interactions in vivo

The sorting of RNA transcripts dictates their ultimate post-transcriptional fates, such as translation, decay or degradation by RNA interference (RNAi). This sorting of RNAs into distinct fates is mediated by their interaction with RNA-binding proteins. While hundreds of RNA binding proteins have been identified, which act to sort RNAs into different pathways is largely unknown. Particularly in plants, this is due to the lack of reliable protein-RNA artificial tethering tools necessary to determine the mechanism of protein action on an RNA in vivo. Here we generated a protein-RNA tethering system which functions on an endogenous Arabidopsis RNA that is tracked by the quantitative flowering time phenotype. Unlike other protein-RNA tethering systems that have been attempted in plants, our system circumvents the inadvertent triggering of RNAi. We successfully in vivo tethered a protein epitope, deadenylase protein and translation factor to the target RNA, which function to tag, decay and boost protein production, respectively. We demonstrated that our tethering system (1) is sufficient to engineer the downstream fate of an RNA, (2) enables the determination of any protein's function upon recruitment to an RNA, and (3) can be used to discover new interactions with RNA-binding proteins.

Multiple gene expression in plants using MIDAS-P, a versatile type II restriction-based modular expression vector

MIDAS-P is a plant expression vector with blue/white screening for iterative cloning of multiple, tandemly arranged transcription units (TUs). We have used the MIDAS-P system to investigate the expression of up to five genes encoding three anti-HIV proteins and the reporter gene DsRed in Nicotiana benthamiana plants. The anti-HIV cocktail was made up of a broadly neutralizing monoclonal antibody (VRC01), a lectin (Griffithsin), and a single-chain camelid nanobody (J3-VHH). Constructs containing different combinations of 3, 4, or 5 TUs encoding different components of the anti-HIV cocktail were assembled. Messenger RNA (mRNA) levels of the genes of interest decreased beyond two TUs. Coexpression of the RNA silencing suppressor P19 dramatically increased the overall mRNA and protein expression levels of each component. The position of individual TUs in 3 TU constructs did not affect mRNA or protein expression levels. However, their expression dropped to non-detectable levels in constructs with four or more TUs each containing the same promoter and terminator elements, with the exception of DsRed at the first or last position in 5 TU constructs. This drop was alleviated by co-expression of P19. In short, the MIDAS-P system is suitable for the simultaneous expression of multiple proteins in one construct.

Arabidopsis ABCG27 plays an essential role in flower and leaf development by modulating abscisic acid content

Abscisic acid (ABA) is a phytohormone that mediates stress responses and regulates plant development. Several ATP-binding cassette (ABC) transporters in the G subfamily of ABC (ABCG) proteins have been reported to transport ABA. We investigated whether there are any other ABCG proteins that mediate plant developmental processes regulated by ABA in Arabidopsis (Arabidopsis thaliana). The ABCG27 gene was upregulated in response to exogenous ABA treatment. The abcg27 knockout mutant exhibited two developmental defects: epinastic leaves and abnormally long pistils, which reduced fertility and silique length. ABCG27 expression was induced threefold when flower buds were exposed to exogenous ABA, and the promoter of ABCG27 had two ABA-responsive elements. ABA content in the pistil and true leaves were increased in the abcg27 knockout mutant. Detached abcg27 pistils exposed to exogenous ABA grew longer than those of the wild-type control. ABCG27 fused to GFP localized to the plasma membrane when expressed in Arabidopsis mesophyll protoplasts. A transcriptome analysis of the pistils and true leaves of the wild type and abcg27 knockout mutant revealed that the expression of organ development-related genes changed in the knockout mutant. In particular, the expression of trans-acting small interference (ta-si) RNA processing enzyme genes, which regulate flower and leaf development, was low in the knockout mutant. Together, these results suggest that ABCG27 most likely function as an ABA transporter at the plasma membrane, modulating ABA levels and thereby regulating the development of the pistils and leaves under normal, non-stressed conditions.

PtncRNAdb: plant transfer RNA-derived non-coding RNAs (tncRNAs) database

Specific endonucleolytic cleavage of tRNA molecules leads to the biogenesis of heterogeneously sized fragments called tRNA-derived non-coding RNAs (tncRNAs). The role of tncRNAs is well studied in human processes, and diseases including different types of cancers and other ailments. They are also generated under stress conditions in plants. Considering the potential role of tncRNAs in the plant system, we have developed a user-friendly, open-access web resource, PtncRNAdb (https://nipgr.ac.in/PtncRNAdb). PtncRNAdb consists of 4,809,503 tncRNA entries identified from ~ 2500 single-end small RNA-seq libraries from six plants, viz., Arabidopsis thaliana, Cicer arietinum, Zea mays, Oryza sativa, Medicago truncatula, and Solanum lycopersicum. It is provided with assorted options to search, browse, visualize, interpret, and download tncRNAs data. Users can perform query search using 'BLASTN' against PtncRNAdb entries. Highcharts have been included for better statistical PtncRNAdb data readability to the users. Additionally, PtncRNAdb includes 'DE tncRNAs' module for differentially expressed tncRNAs under various conditions. Their secondary structure, putative targets, interactive networks of target enrichment, and related publications are also incorporated for further interpretation of their biological functions. PtncRNAdb is an efficient, user-friendly, and exhaustive database, which will aid the ongoing research in plant tncRNAs as well as help in deciphering their role in gene regulation. We hope that it provides a promising platform for researchers to facilitate the understanding of tncRNAs, and their involvement in numerous pathways related to plant development and stress tolerance.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-022-03174-7.

Transcriptome Profiling Reveals a Petunia Transcription Factor, PhCOL4, Contributing to Antiviral RNA Silencing

RNA silencing is a common antiviral mechanism in eukaryotic organisms. However, the transcriptional regulatory mechanism that controls the RNA silencing process remains elusive. Here, we performed high-depth transcriptome analysis on petunia (Petunia hybrida) leaves infected with tobacco rattle virus (TRV) strain PPK20. A total of 7,402 differentially expressed genes (DEGs) were identified. Of them, some RNA silencing-related transcripts, such as RNA-dependent RNA polymerases (RDRs), Dicer-like RNase III enzymes (DCLs), and Argonautes (AGOs), were induced by viral attack. Furthermore, we performed TRV-based virus-induced gene silencing (VIGS) assay on 39 DEGs encoding putative transcription factors (TFs), using green fluorescent protein (GFP) and phytoene desaturase (PhPDS) as reporters. Results showed that the down-regulation of PhbHLH41, PhbHLH93, PhZPT4-3, PhCOL4, PhHSF-B3A, PhNAC90, and PhWRKY75 led to enhanced TRV accumulation and inhibited PhPDS-silenced photobleaching phenotype. In contrast, silencing of PhERF22 repressed virus accumulation and promoted photobleaching development. Thus, these TFs were identified as potential positive and negative regulators of antiviral RNA silencing, respectively. One positive regulator PhCOL4, belonging to the B-box zinc finger family, was selected for further functional characterization. Silencing and transient overexpression of PhCOL4 resulted in decreased and increased expression of several RNA silencing-related genes. DNA affinity purification sequencing analysis revealed that PhCOL4 targeted PhRDR6 and PhAGO4. Dual luciferase and yeast one-hybrid assays determined the binding of PhCOL4 to the PhRDR6 and PhAGO4 promoters. Our findings suggest that TRV-GFP-PhPDS-based VIGS could be helpful to identify transcriptional regulators of antiviral RNA silencing.

Effects of Elevated Temperature on the Susceptibility of Capsicum Plants to Capsicum Chlorosis Virus Infection

Capsicum, an important vegetable crop in Queensland, Australia, is vulnerable to both elevated temperatures and capsicum chlorosis virus (CaCV). Thus, it is imperative to understand the genetic responses of capsicum plants (Capsicum annuum) to CaCV under elevated temperature conditions. Here, we challenged susceptible plants (cv. Yolo Wonder) with CaCV and investigated the effects of elevated temperature on symptom expression, the accumulation of virus-derived short interfering RNA (vsiRNA) and viral RNA, and the expression of plant defense-associated genes. CaCV-inoculated plants initially showed more severe symptoms and higher viral concentrations at a higher temperature (HT, 35 °C) than at ambient temperature (AT, 25 °C). However, symptom recovery and reduced viral RNA accumulation were seen in the CaCV-infected plants grown at HT at later stages of infection. We also observed that HT enhanced the accumulation of vsiRNAs and that, concurrently, RNA interference (RNAi)-related genes, including Dicer-like2 (DCL2), DCL4, RNA-dependent RNA polymerase 1 (RdRp1), RdRp6, and Argonaute2 (AGO2), were upregulated early during infection. Moreover, continuous high levels of vsiRNAs were observed during later stages of CaCV infection at HT. Overall, our investigation suggests that HT facilitates CaCV replication during early infection stages. However, this appears to lead to an early onset of antiviral RNA silencing, resulting in a subsequent recovery from CaCV in systemic leaves.

Metabolic-based insecticide resistance mechanism and ecofriendly approaches for controlling of beet armyworm Spodoptera exigua: a review

The beet army worm, Spodoptera exigua, is a widely distributed polyphagous pest of economically important crops worldwide. The management of this pest insect continues to face many challenges. Despite synthetic chemicals posing a serious threat to the environment, these remain the conventional approach for controlling S. exigua in the field. An over-reliance on chemical control has not only led to selection for resistance to insecticides and to a reduction of natural enemies, but has also polluted various components of ecosystem. Given these increasing pressures on the ecosystem, there is a need to implement integrated pest management (IPM) approaches exploiting a wider range of tools (biotechnological approaches, microbial control, biological control, cultural control, and use of host plant resistance) for an alternative to chemical control. The IPM approach can not only reduce the hazard of chemical residues in the environment and associated health problems, but may also provide best strategies to control insect pests. This review synthesizes published information on insecticide resistance of S. exigua and explores alternative IPM approaches to control S. exigua.

Plant RNA-mediated gene regulatory network

Not all transcribed RNAs are protein-coding RNAs. Many of them are non-protein-coding RNAs in diverse eukaryotes. However, some of them seem to be non-functional and are resulted from spurious transcription. A lot of non-protein-coding transcripts have a significant function in the translation process. Gene expressions depend on complex networks of diverse gene regulatory pathways. Several non-protein-coding RNAs regulate gene expression in a sequence-specific system either at the transcriptional level or post-transcriptional level. They include a significant part of the gene expression regulatory network. RNA-mediated gene regulation machinery is evolutionarily ancient. They well-evolved during the evolutionary time and are becoming much more complex than had been expected. In this review, we are trying to summarizing the current knowledge in the field of RNA-mediated gene silencing.

Harnessing epigenetic variability for crop improvement: current status and future prospects

BACKGROUND

The epigenetic mechanisms play critical roles in a vast diversity of biological processes of plants, including development and response to environmental challenges. Particularly, DNA methylation is a stable epigenetic signature that supplements the genetics-based view of complex life phenomena. In crop breeding, the decrease in genetic diversity due to artificial selection of conventional breeding methods has been a long-standing concern. Therefore, the epigenetic diversity has been proposed as a new resource for future crop breeding, which will be hereinafter referred to as epibreeding.

DISCUSSION

The induction of methylome changes has been performed in plants by several methods including chemical drugs treatment and tissue culture. Target-specific epigenetic engineering has been also attempted by exogenous RNAi mediated by virus-induced gene silencing and grafting. Importantly, the new and innovative techniques including the CRISPR-Cas9 system have recently been adopted in epigenetic engineering of plant genomes, facilitating the efforts for epibreeding.

CONCLUSION

In this review, we introduce several examples of natural and induced epigenetic changes impacting on agronomic traits and discuss the methods for generating epigenomic diversity and site-specific epigenetic engineering.

Genome-Wide Identification of Wheat ZIP Gene Family and Functional Characterization of the TaZIP13-B in Plants

The ZIP (Zn-regulated, iron-regulated transporter-like protein) transporter plays an important role in regulating the uptake, transport, and accumulation of microelements in plants. Although some studies have identified ZIP genes in wheat, the significance of this family is not well understood, particularly its involvement under Fe and Zn stresses. In this study, we comprehensively characterized the wheat ZIP family at the genomic level and performed functional verification of three TaZIP genes by yeast complementary analysis and of TaZIP13-B by transgenic Arabidopsis. Totally, 58 TaZIP genes were identified based on the genome-wide search against the latest wheat reference (IWGSC_V1.1). They were then classified into three groups, based on phylogenetic analysis, and the members within the same group shared the similar exon-intron structures and conserved motif compositions. Expression pattern analysis revealed that the most of TaZIP genes were highly expressed in the roots, and nine TaZIP genes displayed high expression at grain filling stage. When exposed to ZnSO₄ and FeCl₃ solutions, the TaZIP genes showed differential expression patterns. Additionally, six ZIP genes responded to zinc-iron deficiency. A total of 57 miRNA-TaZIP interactions were constructed based on the target relationship, and three miRNAs were downregulated when exposed to the ZnSO₄ and FeCl₃ stresses. Yeast complementation analysis proved that TaZIP14-B, TaZIP13-B, and TaIRT2-A could transport Zn and Fe. Finally, overexpression of TaZIP13-B in Arabidopsis showed that the transgenic plants displayed better tolerance to Fe/Zn stresses and could enrich more metallic elements in their seeds than wild-type Arabidopsis. This study systematically analyzed the genomic organization, gene structure, expression profiles, regulatory network, and the biological function of the ZIP family in wheat, providing better understanding of the regulatory roles of TaZIPs and contributing to improve nutrient quality in wheat crops.

Influence of virus-host interactions on plant response to abiotic stress

Environmental factors play a significant role in controlling growth, development and defense responses of plants. Changes in the abiotic environment not only significantly alter the physiological and molecular pathways in plants, but also result in attracting the insect pests that carry a payload of viruses. Invasion of plants by viruses triggers the RNA silencing based defense mechanism in plants. In counter defense the viruses have gained the ability to suppress the host RNA silencing activities. A new paradigm has emerged, with the recognition that plant viruses also have the intrinsic capacity to modulate host plant response to environmental cues, in an attempt to favour their own survival. Thus, plant-virus interactions provide an excellent system to understand the signals in crosstalk between biotic (virus) and abiotic stresses. In this review, we have summarized the basal plant defense responses to pathogen invasion while emphasizing on the role of RNA silencing as a front line of defense response to virus infection. The emerging knowledge indicates overlap between RNA silencing with the innate immune responses during antiviral defense. The suppressors of RNA silencing serve as Avr proteins, which can be recognized by the host R proteins. The defense signals also function in concert with the phytohormones to influence plant responses to abiotic stresses. The current evidence on the role of virus induced host tolerance to abiotic stresses is also discussed.

Transfer RNA-derived non-coding RNAs (tncRNAs): Hidden regulation of plants' transcriptional regulatory circuits

The emergence of distinct classes of non-coding RNAs has led to better insights into the eukaryotic gene regulatory networks. Amongst them, the existence of transfer RNA (tRNA)-derived non-coding RNAs (tncRNAs) demands exploration in the plant kingdom. We have designed a methodology to uncover the entire perspective of tncRNAome in plants. Using this pipeline, we have identified diverse tncRNAs with a size ranging from 14 to 50 nucleotides (nt) by utilizing 2448 small RNA-seq samples from six angiosperms, and studied their various features, including length, codon-usage, cleavage pattern, and modified tRNA nucleosides. Codon-dependent generation of tncRNAs suggests that the tRNA cleavage is highly specific rather than random tRNA degradation. The nucleotide composition analysis of tncRNA cleavage positions indicates that they are generated through precise endoribonucleolytic cleavage machinery. Certain nucleoside modifications detected on tncRNAs were found to be conserved across the plants, and hence may influence tRNA cleavage, as well as tncRNA functions. Pathway enrichment analysis revealed that common tncRNA targets are majorly enriched during metabolic and developmental processes. Further distinct tissue-specific tncRNA clusters highlight their role in plant development. Significant number of tncRNAs differentially expressed under abiotic and biotic stresses highlights their potential role in stress resistance. In summary, this study has developed a platform that will help in the understanding of tncRNAs and their involvement in growth, development, and response to various stresses. The workflow, software package, and results are freely available at http://nipgr.ac.in/tncRNA.

A comprehensive overview of miRNA targeting drought stress resistance in plants

MicroRNAs (miRNAs) are essential nonprotein-coding genes. In a range of organisms, miRNAs has been reported to play an essential role in regulating gene expressions at post-transcriptional level. They participate in most of the stress responsive processes in plants. Drought is an ultimate abiotic stress that affects the crop production. Therefore understanding drought stress responses are essential to improve the production of agricultural crops. Throughout evolution, plants have developed their own defense systems to cope with the adversities of environmental stresses. Among defensive mechanisms include the regulations of gene expression by miRNAs. Drought stress regulates the expression of some of the functionally conserved miRNAs in different plants. The given properties of miRNAs provide an insight to genetic alterations and enhancing drought resistance in cereal crops. The current review gives a summary to regulatory mechanisms in plants as well as miRNAs response to drought stresses in cereal crops. Some possible approaches and guidelines for the exploitation of drought stress miRNA responses to improve cereal crops are also described.

Genome-wide identification of DCL, AGO and RDR gene families and their associated functional regulatory elements analyses in banana (Musa acuminata)

RNA silencing is mediated through RNA interference (RNAi) pathway gene families, i.e., Dicer-Like (DCL), Argonaute (AGO), and RNA-dependent RNA polymerase (RDR) and their cis-acting regulatory elements. The RNAi pathway is also directly connected with the post-transcriptional gene silencing (PTGS) mechanism, and the pathway controls eukaryotic gene regulation during growth, development, and stress response. Nevertheless, genome-wide identification of RNAi pathway gene families such as DCL, AGO, and RDR and their regulatory network analyses related to transcription factors have not been studied in many fruit crop species, including banana (Musa acuminata). In this study, we studied in silico genome-wide identification and characterization of DCL, AGO, and RDR genes in bananas thoroughly via integrated bioinformatics approaches. A genome-wide analysis identified 3 MaDCL, 13 MaAGO, and 5 MaRDR candidate genes based on multiple sequence alignment and phylogenetic tree related to the RNAi pathway in banana genomes. These genes correspond to the Arabidopsis thaliana RNAi silencing genes. The analysis of the conserved domain, motif, and gene structure (exon-intron numbers) for MaDCL, MaAGO, and MaRDR genes showed higher homogeneity within the same gene family. The Gene Ontology (GO) enrichment analysis exhibited that the identified RNAi genes could be involved in RNA silencing and associated metabolic pathways. A number of important transcription factors (TFs), e.g., ERF, Dof, C2H2, TCP, GATA and MIKC_MADS families, were identified by network and sub-network analyses between TFs and candidate RNAi gene families. Furthermore, the cis-acting regulatory elements related to light-responsive (LR), stress-responsive (SR), hormone-responsive (HR), and other activities (OT) functions were identified in candidate MaDCL, MaAGO, and MaRDR genes. These genome-wide analyses of these RNAi gene families provide valuable information related to RNA silencing, which would shed light on further characterization of RNAi genes, their regulatory elements, and functional roles, which might be helpful for banana improvement in the breeding program.

Identification of Genes and MicroRNAs Affecting Pre-harvest Sprouting in Rice (Oryza sativa L.) by Transcriptome and Small RNAome Analyses

Pre-harvest sprouting (PHS) is one of the primary problems associated with seed dormancy in rice (Oryza sativa L.). It causes yield loss and reduces grain quality under unpredictable humid conditions at the ripening stage, thus affecting the economic value of the rice crop. To resolve this issue, understanding the molecular mechanism underlying seed dormancy in rice is important. Recent studies have shown that seed dormancy is affected by a large number of genes associated with plant hormone regulation. However, understanding regarding the effect of heat stress on seed dormancy and plant hormones is limited. This study compared the transcriptome and small RNAome of the seed embryo and endosperm of two contrasting japonica rice accessions, PHS susceptible (with low seed dormancy) and PHS resistant (with high seed dormancy), at three different maturation stages. We found that 9,068 genes and 35 microRNAs (miRNAs) were differentially expressed in the embryo, whereas 360 genes were differentially expressed in the endosperm. Furthermore, we identified and verified the candidate genes associated with seed dormancy and heat stress-related responses in rice using quantitative real-time PCR. We newly discovered eight hormone-related genes, four heat shock protein-related genes, and two miRNAs potentially involved in PHS. These findings provide a strong foundation for understanding the dynamics of transcriptome and small RNAome of hormone- and heat stress-related genes, which affect PHS during seed maturation.

Epigenetic approaches to crop breeding: current status and perspectives

In order to tackle the cumulative adverse effects of global climate change, reduced farmland, and heightened needs of an ever-increasing world population, modern agriculture is in urgent search of solutions that can ensure world food security and sustainable development. Classical crop breeding is still a powerful method to obtain crops with valued agronomical traits, but its potential is gradually being compromised by the menacing decline of genetic variation. Resorting to the epigenome as a source of variation could serve as a promising alternative. Here, we discuss current status of epigenetics-mediated crop breeding (epibreeding), highlight its advances and limitations, outline currently available methodologies, and propose novel RNA-based strategies to modify the epigenome in a gene-specific and transgene-free manner.

Plant-derived xenomiRs and cancer: Cross-kingdom gene regulation

Exosomal microRNAs (miRNAs) critically regulate several major intracellular and metabolic activities, including cancer evolution. Currently, increasing evidence indicates that exosome harbor and transport these miRNAs from donor cells to neighboring and distantly related recipient cells, often in a cross-species manner. Several studies have reported that plant-based miRNAs can be absorbed into the serum of humans, where they hinder the expression of human disease-related genes. Moreover, few recent studies have demonstrated the role of these xenomiRs in cancer development and progression. However, the cross-kingdom gene regulation hypothesis remains highly debatable, and many follow up studies fail to reproduce the same. There are reports that show no effect of plant-derived miRNAs on mammalian cancers. The foremost cause of this controversy remains the lack of reproducibility of the results. Here, we reassess the latest developments in the field of cross-kingdom transference of miRNAs, emphasizing on the role of the diet-based xenomiRs on cancer progression.

A natural antisense transcript acts as a negative regulator for the maize drought stress response gene ZmNAC48

Although plant-specific NAC transcription factors play crucial roles in response to abiotic stress, few reports describe the regulation of NAC genes in maize (Zea mays) by the cis-natural antisense transcripts (cis-NATs). In this study, 521 NAC genes from Gramineae were classified, of which 51 NAC genes contained cis-NATs. ZmNAC48 and cis-NATZmNAC48 co-localized to the same cell nucleus, and both transcripts responded to drought stress. Arabidopsis plants overexpressing ZmNAC48 had improved drought tolerance, lower rate of water loss, enhanced stomatal closure, and higher rates of survival. Transient expression in both maize protoplasts and tobacco leaves indicated that cis-NATZmNAC48 reduced ZmNAC48 expression. Western blotting and ribosome profiling analyses confirmed that cis-NATZmNAC48 lacked protein coding potential. Furthermore, the cis-NAT-derived small-interfering RNAs (nat-siRNAs) generated from the overlapping regions of ZmNAC48 and cis-NATZmNAC48 were detected in maize and transgenic Arabidopsis. Cis-NATZmNAC48 overexpressing maize showed higher water loss rate, increased stomatal opening, and had more dead leaves. Expression of ZmNAC48 and nat-siRNA was decreased in these plants. Taken together, our study indicates that both ZmNAC48 and cis-NATZmNAC48 are involved in plant drought stress responses, and that the double-stranded RNA-dependent mechanism is involved in the interaction between cis-NATZmNAC48 and ZmNAC48. Additionally, cis-NATZmNAC48 may negatively regulate ZmNAC48 to affect stomatal closure of maize.

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.

Identification of Phosphorylated Amino Acids in Human TNRC6A C-Terminal Region and Their Effects on the Interaction with the CCR4-NOT Complex

Human GW182 family proteins have Argonaute (AGO)-binding domains in their N-terminal regions and silencing domains, which interact with RNA silencing-related proteins, in their C-terminal regions. Thus, they function as scaffold proteins between the AGO protein and RNA silencing-related proteins, such as carbon catabolite repressor4-negative on TATA (CCR4-NOT) or poly(A)-binding protein (PABP). Our mass spectrometry analysis and the phosphorylation data registered in PhosphoSitePlus, a post-translational modification database, suggested that the C-terminal region of a human GW182 family protein, TNRC6A, has at least four possible phosphorylation sites, which are located near the region interacting with the CCR4-NOT complex. Among them, two serine residues at amino acid positions 1332 and 1346 (S1332 and S1346) were certainly phosphorylated in human HeLa cells, but other two serine residues (S1616 and S1691) were not phosphorylated. Furthermore, it was revealed that the phosphorylation patterns of TNRC6A affect the interaction with the CCR4-NOT complex. When S1332 and S1346 were dephosphorylated, the interactions of TNRC6A with the CCR4-NOT complex were enhanced, and when S1616 and S1691 were phosphorylated, such interaction was suppressed. Thus, phosphorylation of TNRC6A was considered to regulate the interaction with RNA silencing-related factors that may affect RNA silencing activity.

R2R3-MYB transcription factors, StmiR858 and sucrose mediate potato flavonol biosynthesis

Flavonols and other phenylpropanoids protect plants from biotic and abiotic stress and are dietarily desirable because of their health-promoting properties. The ability to develop new potatoes (Solanum tuberosum) with optimal types and amounts of phenylpropanoids is limited by lack of knowledge about the regulatory mechanisms. Exogenous sucrose increased flavonols, whereas overexpression of the MYB StAN1 induced sucrolytic gene expression. Heterologous StAN1 protein bound promoter fragments from sucrolytic genes (SUSY1 and INV1). Two additional MYBs and one microRNA were identified that regulated potato flavonols. Overexpression analysis showed MYB12A and C increased amounts of flavonols and other phenylpropanoids. Endogenous flavonol amounts in light-exposed organs were much higher those in the dark. Expression levels of StMYB12A and C were high in flowers but low in tubers. Transient overexpression of miR858 altered potato flavonol metabolism. Endogenous StmiR858 expression was much lower in flowers than leaves and correlated with flavonol amounts in these organs. Collectively, these findings support the hypothesis that sucrose, MYBs, and miRNA control potato phenylpropanoid metabolism in a finely tuned manner that includes a feedback loop between sucrose and StAN1. These findings will aid in the development of potatoes with phenylpropanoid profiles optimized for crop performance and human health.

Genome-wide miRNA expression profiling in potato (Solanum tuberosum L.) reveals TOR-dependent post-transcriptional gene regulatory networks in diverse metabolic pathway

Target of rapamycin (TOR) operates as a hub of the signal transduction that integrates nutrient and energy signaling to promote cell proliferation and growth through mediating the transcriptional and post- transcriptional regulator networks in all eukaryotic species. MicroRNAs (miRNAs) are widespread classes of small, single-stranded, non-coding endogenous RNAs and are widely found in eukaryotes, which play a vital role in regulating gene expression by degrading targeted mRNAs or translational repression at post-transcriptional level. Recent studies found that there were necessarily close connections between miRNA and TOR pathways in mammals. However, there is little information about the interplay between the miRNA and TOR in plants. Thus, the aim of this study was to identify potential TOR-miRNA-mRNA regulatory networks in TOR signaling through global mRNA and microRNA expression profiling in potato. Based on the previous high-throughput transcriptome sequencing and filtering, a total of 2,899 genes were significantly differentially expressed in potato under TOR inhibitors treatment. Pathway analysis revealed that these genes were significantly enriched in multiple metabolic processes. Similarly, in the present study, suppression of TOR resulted in 41 miRNAs up-regulated and 45 down-regulated, revealing that TOR plays a crucial role in the regulation of miRNA regulatory network. Furthermore, integrated mRNA and miRNA expression profiling uncovered that these miRNAs participated in large-scale metabolic process in the TOR signal pathway in potato, such as regulation of autophagy and ubiquitination, and biosynthesis of secondary metabolites. Overall, the results shed new insight into TOR related post-transcriptional gene regulatory networks in potato and suggesting TOR-miRNA-targeting genes relevant networks as a potential genetic resource for potato improvement.

A Cotton Lignin Biosynthesis Gene, GhLAC4, Fine-Tuned by ghr-miR397 Modulates Plant Resistance Against Verticillium dahliae

Plant lignin is a component of the cell wall, and plays important roles in the transport potential of water and mineral nutrition and plant defence against biotic stresses. Therefore, it is necessary to identify lignin biosynthesis-related genes and dissect their functions and underlying mechanisms. Here, we characterised a cotton LAC, GhLAC4, which participates in lignin biosynthesis and plant resistance against Verticillium dahliae. According to degradome sequencing and GUS reporter analysis, ghr-miR397 was identified to directedly cleave the GhLAC4 transcript through base complementary. GhLAC4 knockdown and ghr-miR397 overexpression significantly reduced basal lignin content compared to the control, whereas ghr-miR397 silencing significantly increased basal lignin levels. Based on staining patterns and GC/MS analysis, GhLAC4 acted in G-lignin biosynthesis. Under V. dahliae infection, we found that G-lignin content in ghr-miR397-knockdowned plants significantly increased, compared to these plants under the mock treatment, while G-lignin contents in GhLAC4-silenced plants and ghr-miR397-overexpressed plants treated with pathogen were comparable with these plants treated with mock, indicating that GhLAC4 participates in defence-induced G-lignin biosynthesis in the cell wall. Knockdown of ghr-miR397 in plants inoculated with V. dahliae promoted lignin accumulation and increased plant resistance. The overexpression of ghr-miR397 and knockdown of GhLAC4 reduced lignin content and showed higher susceptibility of plants to the fungal infection compared to the control. The extract-free stems of ghr-miR397-knockdowned plants lost significantly less weight when treated with commercial cellulase and V. dahliae secretion compared to the control, while the stems of ghr-miR397-overexpressed and GhLAC4-silenced plants showed significantly higher loss of weight. These results suggest that lignin protects plant cell walls from degradation mediated by cellulase or fungal secretions. In summary, the ghr-miR397-GhLAC4 module regulates both basal lignin and defence-induced lignin biosynthesis and increases plant resistance against infection by V. dahliae.

Elevated Carbon Dioxide Levels Decreases Cucumber Mosaic Virus Accumulation in Correlation with Greater Accumulation of rgs-CaM, an Inhibitor of a Viral Suppressor of RNAi

Plant viruses cause a range of plant diseases symptoms that are often responsible for significant crop production losses and the severity and spread of the symptoms may be affected by climate change. While the increase in anthropogenic activities has caused a critical problem of increased CO2 levels in the atmosphere, these elevated CO2 levels have been reported to reduce virus disease severity in some plant species. In such instances, it is not clear if the plant defense mechanisms are being enhanced or virus-mediated mechanisms to overcome plant resistance are being defeated. Additionally, a few studies have been attempted in this area to determine if reduced disease is the norm or the exception under enhanced CO2 levels. In the present study, the effects of elevated CO2 levels (750 ppm vs. 390 ppm) on RNAi-mediated resistance of Nicotiana tabacum against the cucumber mosaic virus (CMV), and the activity of viral suppressor of RNAi (VSR) 2b protein of CMV were evaluated. On the one hand, our results showed that elevated CO2 decreased the transcription of dicer-like protein 2 (DCL2), DCL4, and argonaut 1 (AGO1) genes with functions related to RNAi-mediated resistance when infected by CMV, which is contradictory with the decreased CMV copy numbers under elevated CO2. On the other hand, we found that elevated CO2 increased the calcium concentration and expression of the calcium-binding protein rgs-CaM in tobacco plants when infected by CMV, which directly weakened the function of 2b protein, the VSR of CMV, and therefore decreased the infection efficiency of the virus and suppressed the severity of CMV in tobacco plants under elevated CO2. This study provides molecular insights into the ecological implications underlying the development of prevention strategies against plant virus infection in the context of climate change.

In silico identification and characterization of AGO, DCL and RDR gene families and their associated regulatory elements in sweet orange (Citrus sinensis L.)

RNA interference (RNAi) plays key roles in post-transcriptional and chromatin modification levels as well as regulates various eukaryotic gene expressions which are involved in stress responses, development and maintenance of genome integrity during developmental stages. The whole mechanism of RNAi pathway is directly involved with the gene-silencing process by the interaction of Dicer-Like (DCL), Argonaute (AGO) and RNA-dependent RNA polymerase (RDR) gene families and their regulatory elements. However, these RNAi gene families and their sub-cellular locations, functional pathways and regulatory components were not extensively investigated in the case of economically and nutritionally important fruit plant sweet orange (Citrus sinensis L.). Therefore, in silico characterization, gene diversity and regulatory factor analysis of RNA silencing genes in C. sinensis were conducted by using the integrated bioinformatics approaches. Genome-wide comparison analysis based on phylogenetic tree approach detected 4 CsDCL, 8 CsAGO and 4 CsRDR as RNAi candidate genes in C. sinensis corresponding to the RNAi genes of model plant Arabidopsis thaliana. The domain and motif composition and gene structure analyses for all three gene families exhibited almost homogeneity within the same group members. The Gene Ontology enrichment analysis clearly indicated that the predicted genes have direct involvement into the gene-silencing and other important pathways. The key regulatory transcription factors (TFs) MYB, Dof, ERF, NAC, MIKC_MADS, WRKY and bZIP were identified by their interaction network analysis with the predicted genes. The cis-acting regulatory elements associated with the predicted genes were detected as responsive to light, stress and hormone functions. Furthermore, the expressed sequence tag (EST) analysis showed that these RNAi candidate genes were highly expressed in fruit and leaves indicating their organ specific functions. Our genome-wide comparison and integrated bioinformatics analyses provided some necessary information about sweet orange RNA silencing components that would pave a ground for further investigation of functional mechanism of the predicted genes and their regulatory factors.

RNAi Mediated Hypoxia Stress Tolerance in Plants

Small RNAs regulate various biological process involved in genome stability, development, and adaptive responses to biotic or abiotic stresses. Small RNAs include microRNAs (miRNAs) and small interfering RNAs (siRNAs). MicroRNAs (miRNAs) are regulators of gene expression that affect the transcriptional and post-transcriptional regulation in plants and animals through RNA interference (RNAi). miRNAs are endogenous small RNAs that originate from the processing of non-coding primary miRNA transcripts folding into hairpin-like structures. The mature miRNAs are incorporated into the RNA-induced silencing complex (RISC) and drive the Argonaute (AGO) proteins towards their mRNA targets. siRNAs are generated from a double-stranded RNA (dsRNA) of cellular or exogenous origin. siRNAs are also involved in the adaptive response to biotic or abiotic stresses. The response of plants to hypoxia includes a genome-wide transcription reprogramming. However, little is known about the involvement of RNA signaling in gene regulation under low oxygen availability. Interestingly, miRNAs have been shown to play a role in the responses to hypoxia in animals, and recent evidence suggests that hypoxia modulates the expression of various miRNAs in plant systems. In this review, we describe recent discoveries on the impact of RNAi on plant responses to hypoxic stress in plants.

Characterization of Hibiscus Latent Fort Pierce Virus-Derived siRNAs in Infected Hibiscus rosa-sinensis in China

Although limited progress have been made about pathogen system of Hibiscus rosa-sinensis and Hibiscus latent Fort Pierce virus (HLFPV), interaction between plant host and pathogen remain largely unknown, which led to deficiency of effective measures to control disease of hibiscus plants caused by HLFPV. In this study, infection of HLFPV in Hibiscus rosa-sinensis was firstly confirmed for the first time by traditional electron microscopy, modern reverse transcription polymerase chain reaction and RNA-seq methods in China (HLFPV-Ch). Sequence properties analyzing suggested that the full-length sequences (6,465 nt) of HLFPV-Ch had a high sequence identity and a similar genomic structure with other tobamoviruses. It includes a 5'-terminal untranslated region (UTR), followed by four open reading frames encoding for a 128.5-kDa replicase, a 186.5-kDa polymerase, a 31-kDa movement protein, 17.6-kDa coat protein, and the last a 3'-terminal UTR. Furthermore, HLFPV-Ch-derived virus-derived siRNAs (vsiRNAs) ant its putative target genes, reported also for the first time, were identified and characterized from disease Hibiscus rosa-sinensis through sRNA-seq and Patmatch server to investigate the interaction in this pathogen systems. HLFPV-Ch-derived vsiRNAs demonstrated several general and specific characteristics. Gene Ontology classification revealed predicted target genes by vsiRNAs are involved in abroad range of cellular component, molecular function and biological processes. Taken together, for first time, our results certified the HLFPV infection in China and provide an insight into interaction between HLFPV and Hibiscus rosa-sinensis.

Interactions and links among the noncoding RNAs in plants under stresses

The complex interplay among sRNAs, lncRNAs and circRNAs has been implicated in plants under biotic and abiotic stresses. Here, we review current advances in our understanding of ncRNA interactions and links, which have considerable potential for improving the agronomic traits and the environmental adaptability of plants. Plants can respond to biotic or abiotic stresses. To cope with various conditions, numerous intricate molecular regulatory mechanisms have evolved in plants. Noncoding RNAs (ncRNAs) can be divided into small noncoding RNAs (sRNAs), long noncoding RNAs (lncRNAs) and circular RNAs (circRNAs). Emerging evidence has demonstrated that interplay among the ncRNAs acts as a novel layer in the regulatory mechanisms, which has attracted substantial interest. Links between sRNAs can affect plant immune responses and development in synergistic or antagonistic manners. Additionally, multiple interactions between lncRNAs and sRNAs are involved in crop breeding, disease resistance and high tolerance to environmental stresses. Here, we summarize current knowledge of the interactions and links among the ncRNAs in plant responses to stresses and the methods for identifying ncRNA interactions. Furthermore, challenges and prospects for further progress in elucidating ncRNA interactions and links are highlighted.

The effect of zebularine on the heat-activated retrotransposon ONSEN in Arabidopsis thaliana and Vigna angularis

The Ty1/copia-like retrotransposon ONSEN is conserved among Brassica species, as well as in beans, including adzuki bean (Vigna angularis (Willd.) Ohwi & Ohashi), which is one of the economically important crops in Japan. ONSEN has acquired a heat-responsive element that is recognized by plant heat stress defense factors, resulting in its transcription and the production of full-length extrachromosomal DNA under conditions with elevated temperatures. DNA methylation plays an important role in regulating the activation of this transposon in plants. Therefore, chemical inhibition of DNA methyltransferases has been utilized to study the effect of DNA methylation on transposon activation. To understand the effect of DNA methylation on ONSEN activation, Arabidopsis thaliana and adzuki bean seedlings were treated with zebularine, which is known to be an effective chemical demethylation agent. The results showed that ONSEN transcription levels were upregulated in zebularine-treated plants. Extrachromosomal DNA of ONSEN also accumulated in the treated plants.

NATpare: a pipeline for high-throughput prediction and functional analysis of nat-siRNAs

Natural antisense transcript-derived small interfering RNAs (nat-siRNAs) are a class of functional small RNA (sRNA) that have been found in both plant and animals kingdoms. In plants, these sRNAs have been shown to suppress the translation of messenger RNAs (mRNAs) by directing the RNA-induced silencing complex (RISC) to their sequence-specific mRNA target(s). Current computational tools for classification of nat-siRNAs are limited in number and can be computationally infeasible to use. In addition, current methods do not provide any indication of the function of the predicted nat-siRNAs. Here, we present a new software pipeline, called NATpare, for prediction and functional analysis of nat-siRNAs using sRNA and degradome sequencing data. Based on our benchmarking in multiple plant species, NATpare substantially reduces the time required to perform prediction with minimal resource requirements allowing for comprehensive analysis of nat-siRNAs in larger and more complex organisms for the first time. We then exemplify the use of NATpare by identifying tissue and stress specific nat-siRNAs in multiple Arabidopsis thaliana datasets.

Sound Vibration-Triggered Epigenetic Modulation Induces Plant Root Immunity Against Ralstonia solanacearum

Sound vibration (SV) is one of the several environmental stimuli that induce physiological changes in plants including changes in plant immunity. Immune activation is a complicated process involving epigenetic modifications, however, SV-induced epigenetic modifications remain unexplored. Here, we performed an integrative analysis comprising chromatin immunoprecipitation (ChIP) and microRNA sequencing (miRNA-seq) to understand the role of SV-mediated epigenetic modifications in immune activation in Arabidopsis thaliana against the root pathogen Ralstonia solanacearum. Plants exposed to SV (10 kHz) showed abundant H3K27me3 modification in the promoter regions of aliphatic glucosinolate biosynthesis and cytokinin signaling genes, leading to transcriptional changes that promote immunity. Additionally, 10 kHz SV down-regulated miR397b expression, thus activating three target LACCASE transcripts that mediate cell wall reinforcement via lignin accumulation. Taken together, SV triggers epigenetic modification of genes involved in secondary metabolite biosynthesis, defense hormone signaling, and pre-formed defense in A. thaliana, leading to the activation of plant immunity against R. solanacearum.

Genome-wide analysis of epigenetic and transcriptional changes associated with heterosis in pigeonpea

Hybrids are extensively used in agriculture to deliver an increase in yield, yet the molecular basis of heterosis is not well understood. Global DNA methylation analysis, transcriptome analysis and small RNA profiling were aimed to understand the epigenetic effect of the changes in gene expression level in the two hybrids and their parental lines. Increased DNA methylation was observed in both the hybrids as compared to their parents. This increased DNA methylation in hybrids showed that majority of the 24-nt siRNA clusters had higher expression in hybrids than the parents. Transcriptome analysis revealed that various phytohormones (auxin and salicylic acid) responsive hybrid-MPV DEGs were significantly altered in both the hybrids in comparison to MPV. DEGs associated with plant immunity and growth were overexpressed whereas DEGs associated with basal defence level were repressed. This antagonistic patterns of gene expression might contribute to the greater growth of the hybrids. It was also noticed that some common as well as unique changes in the regulatory pathways were associated with heterotic growth in both the hybrids. Approximately 70% and 67% of down-regulated hybrid-MPV DEGs were found to be differentially methylated in ICPH 2671 and ICPH 2740 hybrid, respectively. This reflected the association of epigenetic regulation in altered gene expressions. Our findings also revealed that miRNAs might play important roles in hybrid vigour in both the hybrids by regulating their target genes, especially in controlling plant growth and development, defence and stress response pathways. The above finding provides an insight into the molecular mechanism of pigeonpea heterosis.

Investigation of genes associated with petal variations between diploid and autotetraploid in Chinese cabbage (Brassica rapa L. ssp. pekinensis) by RNA-seq and sRNA-seq

Polyploidy promotes morphological, physiological, and reproductive diversity in plants. The imminent effect of chromosome doubling in plants is the enlargement of organs such as flowers and fruits, which increases the commercial value of crops. Flowering plays a vital role in the growth and development of angiosperms. Here, we prepared an isolated microspore culture of 'FT', a doubled haploid (DH) line of Chinese cabbage (Brassica rapa L. ssp. pekinensis), and obtained diploid and autotetraploid plants with the same genetic background. Compared with diploids, the autotetraploids were characterized by large floral organs, dark petals, delayed flowering, and reduced fertility. The indole-3-acetic acid (IAA) and jasmonic acid (JA) levels in autotetraploid petals were significantly higher and the abscisic acid (ABA) level was significantly lower than those in the diploid petals. The lutein level in autotetraploid petals was nearly two times higher than that in the diploid petals. A comparative transcriptome analysis revealed 14,412 differentially expressed genes (DEGs) between the diploids and autotetraploids, and they were enriched in 117 Gene Ontology terms and 110 Kyoto Encyclopedia of Genes and Genomes pathways. We detected 231 DEGs related to phytohormone signal transduction and 29 DEGs involved in carotenoid biosynthesis. An miRNA-target mRNA analysis showed that 32 DEGs regulated by 16 DEMs were associated with flowering timing (BraA03000336, BraA09004319, and BraA09000515), petal development (BraA05002408, BraA01004006, BraA09004069, and BraA04000966), flower opening (BraA07000350), and pollen development (BraA01000720, BraA09005727, and BraA01000253). This study provides information to help elucidate the molecular mechanisms underlying phenotypic variations induced by autopolyploidy in Chinese cabbage.