Analyses of RNA-Seq and sRNA-Seq data reveal a complex network of anti-viral defense in TCV-infected Arabidopsis thaliana

Impact of Exogenous dsRNA on miRNA Composition in Arabidopsis thaliana

The application of double-stranded RNAs (dsRNAs) to plant surfaces has emerged as a promising tool for manipulating gene expression in plants and pathogens, offering new opportunities for crop improvement. While research has shown the capability of exogenous dsRNAs to silence genes, the full spectrum of their impact, particularly on the intricate network of microRNAs (miRNAs), remains largely unexplored. Here, we show that the exogenous application of chalcone synthase (CHS)-encoding dsRNA to the rosette leaves of Arabidopsis thaliana induced extensive alterations in the miRNA profile, while non-specific bacterial neomycin phosphotransferase II (NPTII) dsRNA had a minimal effect. Two days after treatment, we detected 60 differentially expressed miRNAs among the 428 miRNAs found in the A. thaliana genome. A total of 59 miRNAs were significantly changed after AtCHS-dsRNA treatment compared with water and NPTII-dsRNA, and 1 miRNA was significantly changed after AtCHS-dsRNA and NPTII-dsRNA compared with the water control. A comprehensive functional enrichment analysis revealed 17 major GO categories enriched among the genes potentially targeted by the up- and downregulated miRNAs. These categories included processes such as aromatic compound biosynthesis (a pathway directly related to CHS activity), heterocycle biosynthesis, RNA metabolism and biosynthesis, DNA transcription, and plant development. Several predicted targets of upregulated and downregulated miRNAs, including APETALA2, SCL27, SOD1, GRF1, AGO2, PHB, and PHV, were verified by qRT-PCR. The analysis showed a negative correlation between the expression of miRNAs and the expression of their predicted targets. Thus, exogenous plant gene-specific dsRNAs induce substantial changes in the plant miRNA composition, ultimately affecting the expression of a wide range of genes. These findings have profound implications for our understanding of the effects of exogenously induced RNA interference, which can have broader effects beyond targeted mRNA degradation, affecting the expression of other genes through miRNA regulation.

Widespread changes to the translational landscape in a maize microRNA biogenesis mutant

MicroRNAs are short, non-coding RNAs that repress gene expression in both plants and animals and have diverse functions related to growth, development, and stress responses. The ribonuclease, DICER-LIKE1 (DCL1) is required for two steps in plant miRNA biogenesis: cleavage of the primary miRNAs (pri-miRNAs) to release a hairpin structure, called the precursor miRNA (pre-miRNA) and cleavage of the pre-miRNA to generate the miRNA/miRNA* duplex. The mature miRNA guides the RNA-induced silencing complex to target RNAs with complementary sequences, resulting in translational repression and/or RNA cleavage of target mRNAs. However, the relative contribution of translational repression versus mRNA degradation by miRNAs remains unknown at the genome-level in crops, especially in maize. The maize fuzzy tassel (fzt) mutant contains a hypomorphic mutation in DCL1 resulting in broad developmental defects. While most miRNAs are reduced in fzt, the levels of miRNA-targeted mRNAs are not dramatically increased, suggesting that translational regulation by miRNAs may be common. To gain insight into the repression mechanism of plant miRNAs, we combined ribosome profiling and RNA-sequencing to globally survey miRNA activities in maize. Our data indicate that translational repression contributes significantly to regulation of most miRNA targets and that approximately one-third of miRNA targets are regulated primarily at the translational level. Surprisingly, ribosomes appear altered in fzt mutants suggesting that DCL1 may also have a role in ribosome biogenesis. Thus, DICER-LIKE1 shapes the translational landscape in plants through both miRNA-dependent and miRNA-independent mechanisms.

Enhancing Pseudomonas syringae pv. Actinidiae sensitivity in kiwifruit by repressing the NBS-LRR genes through miRNA-215-3p and miRNA-29-3p identification

Kiwifruit bacterial canker, caused by Pseudomonas syringae pv. actinidiae (PSA), poses a grave threat to the global kiwifruit industry. In this study, we examined the role of microRNAs (miRNAs) in kiwifruit's response to PSA. Kiwifruit seedlings subjected to PSA treatment showed significant changes in both miRNA and gene expression compared to the control group. We identified 364 differentially expressed miRNAs (DEMs) and 7170 differentially expressed genes (DEGs). Further analysis revealed 180 miRNAs negatively regulating 641 mRNAs. Notably, two miRNAs from the miRNA482 family, miRNA-215-3p and miRNA-29-3p, were found to increase kiwifruit's sensitivity to PSA when overexpressed. These miRNAs were linked to the regulation of NBS-LRR target genes, shedding light on their role in kiwifruit's defence against PSA. This study offers insights into the miRNA482-NBS-LRR network as a crucial component in enhancing kiwifruit bioresistance to PSA infestation and provides promising candidate genes for further research.

RNAi activation with homologous and heterologous sequences that induce resistance against the begomovirus Pepper golden mosaic virus (PepGMV)

This study compared the transcriptional changes in Nicotiana benthamiana plants treated with homologous sequences derived from Pepper golden mosaic virus (PepGMV) and heterologous sequences that derived from another begomovirus, Tomato chino La Paz virus (ToChLPV) prior to infection by PepGMV. The results of microarray analyses identified upregulated genes associated with RNAi such as DCL2, DCL4, AGO3, AGO7, AGO10, NRPD2B (Pol IV), DRB3, CMT3, RDR6. The components that participate in different RNAi pathways were identified, including methylation induced by both constructs, as well as the code of these genes in Arabidopsis thaliana and its counterpart in N. benthamiana through different genome assembly. The expression of these genes was validated by quantitative reverse transcription polymerase chain reaction (RT-qPCR), where DCL3, DCL4, AGO1-1, AGO2, RDR6 and PPR1 showed increased expression during plant protection with the heterologous construct compared to those protected with the homologous construct. The results of this study confirmed the activation of the gene silencing mechanism at the transcriptional level with both constructs and established the possibility of their use as a protection system for both homologous and heterologous sequences.

HASTY modulates miRNA biogenesis by linking pri-miRNA transcription and processing

Post-transcriptional gene silencing mediated by microRNAs (miRNAs) modulates numerous developmental and stress response pathways. For the last two decades, HASTY (HST), the ortholog of human EXPORTIN 5, was considered to be a candidate protein that exports plant miRNAs from the nucleus to the cytoplasm. Here, we report that HST functions in the miRNA pathway independent of its cargo-exporting activity in Arabidopsis. We found that Arabidopsis mutants with impaired HST shuttling exhibit normal subcellular distribution of miRNAs. Interestingly, protein-protein interaction and microscopy assays showed that HST directly interacts with the microprocessor core component DCL1 through its N-terminal domain. Moreover, mass spectrometry analysis revealed that HST also interacts independently of its N-terminal domain with the mediator complex subunit MED37. Further experiments revealed that HST could act as a scaffold to facilitate the recruitment of DCL1 to genomic MIRNA loci by stabilizing the DCL1-MED37 complex, which in turn promotes the transcription and proper processing of primary miRNA transcripts (pri-miRNAs). Taken together, these results suggest that HST is likely associated with the formation of the miRNA biogenesis complex at MIRNA genes, promoting the transcription and processing of pri-miRNAs rather than the direct export of processed miRNAs from the nucleus.

Integrated mRNA and microRNA transcriptome analyses provide insights into paclobutrazol inhibition of lateral branching in herbaceous peony

Herbaceous peony (Paeonia lactiflora Pall.) is a new high-end cut flower, but a large number of lateral branches often appear in some excellent cultivars, which is inconvenient for cut flower production. In the present study, we analyzed the effects of paclobutrazol (PBZ) on the lateral branches of P. lactiflora and adopted a next-generation sequencing approach to identify miRNAs and mRNAs that were differentially expressed involved in the PBZ response. Our results indicate that PBZ may inhibit the production of lateral branches on P. lactiflora. There were 827 differentially expressed genes (DEGs) and 104 differentially expressed miRNAs (DEMs). Integrative analysis revealed 29 miRNA-mRNA interactions related to PBZ stress. Our results provided a wealth of genetic information and data on metabolic pathways for revealing the regulatory mechanism of PBZ inhibition of the development of lateral branches in P. lactiflora and provided a new possibility for reducing lateral branch formation in the production of herbaceous peony cut flowers.

Integrative expression network analysis of microRNA and gene isoforms in sacred lotus

Background

Gene expression is complex and regulated by multiple molecular mechanisms, such as miRNA-mediated gene inhibition and alternative-splicing of pre-mRNAs. However, the coordination of interaction between miRNAs with different splicing isoforms, and the change of splicing isoform in response to different cellular environments are largely unexplored in plants. In this study, we analyzed the miRNA and mRNA transcriptome from lotus (Nelumbo nucifera), an economically important flowering plant.

Results

Through RNA-seq analyses on miRNAs and their target genes (isoforms) among six lotus tissues, expression of most miRNAs seem to be negatively correlated with their targets and tend to be tissue-specific. Further, our results showed that preferential interactions between miRNAs and hub gene isoforms in one coexpression module which is highly correlated with leaf. Intriguingly, for many genes, their corresponding isoforms were assigned to different co-expressed modules, and they exhibited more divergent mRNA structures including presence and absence of miRNA binding sites, suggesting functional divergence for many isoforms is escalated by both structural and expression divergence. Further detailed functional enrichment analysis of miRNA targets revealed that miRNAs are involved in the regulation of lotus growth and development by regulating plant hormone-related pathway genes.

Conclusions

Taken together, our comprehensive analyses of miRNA and mRNA transcriptome elucidate the coordination of interaction between miRNAs and different splicing isoforms, and highlight the functional divergence of many transcript isoforms from the same locus in lotus.

iTRAQ-Based Proteomic Analysis of Watermelon Fruits in Response to Cucumber green mottle mosaic virus Infection

Cucumber green mottle mosaic virus (CGMMV) is an important viral pathogen on cucurbit plants worldwide, which can cause severe fruit decay symptoms on infected watermelon (usually called “watermelon blood flesh”). However, the molecular mechanism of this disease has not been well understood. In this study, we employed the isobaric tags for relative and absolute quantitation (iTRAQ) technique to analyze the proteomic profiles of watermelon fruits in response to CGMMV infection. A total of 595 differentially accumulated proteins (DAPs) were identified, of which 404 were upregulated and 191 were downregulated. Functional annotation analysis showed that these DAPs were mainly involved in photosynthesis, carbohydrate metabolism, secondary metabolite biosynthesis, plant–pathogen interaction, and protein synthesis and turnover. The accumulation levels of several proteins related to chlorophyll metabolism, pyruvate metabolism, TCA cycle, heat shock proteins, thioredoxins, ribosomal proteins, translation initiation factors, and elongation factors were strongly affected by CGMMV infection. Furthermore, a correlation analysis was performed between CGMMV-responsive proteome and transcriptome data of watermelon fruits obtained in our previous study, which could contribute to comprehensively elucidating the molecular mechanism of “watermelon blood flesh”. To confirm the iTRAQ-based proteome data, the corresponding transcripts of ten DAPs were validated by determining their abundance via quantitative reverse transcriptase-polymerase chain reaction (qRT-PCR). These results could provide a scientific basis for in-depth understanding of the pathogenic mechanisms underlying CGMMV-induced “watermelon blood flesh”, and lay the foundation for further functional exploration and verification of related genes and proteins.

Genome-Wide microRNA Profiling Using Oligonucleotide Microarray Reveals Regulatory Networks of microRNAs in Nicotiana benthamiana During Beet Necrotic Yellow Vein Virus Infection

Beet necrotic yellow vein virus (BNYVV) infections induce stunting and leaf curling, as well as root and floral developmental defects and leaf senescence in Nicotiana benthamiana. A microarray analysis with probes capable of detecting 1596 candidate microRNAs (miRNAs) was conducted to investigate differentially expressed miRNAs and their targets upon BNYVV infection of N. benthamiana plants. Eight species-specific miRNAs of N. benthamiana were identified. Comprehensive characterization of the N. benthamiana microRNA profile in response to the BNYVV infection revealed that 129 miRNAs were altered, including four species-specific miRNAs. The targets of the differentially expressed miRNAs were predicted accordingly. The expressions of miR164, 160, and 393 were up-regulated by BNYVV infection, and those of their target genes, NAC21/22, ARF17/18, and TIR, were down-regulated. GRF1, which is a target of miR396, was also down-regulated. Further genetic analysis of GRF1, by Tobacco rattle virus-induced gene silencing, assay confirmed the involvement of GRF1 in the symptom development during BNYVV infection. BNYVV infection also induced the up-regulation of miR168 and miR398. The miR398 was predicted to target umecyanin, and silencing of umecyanin could enhance plant resistance against viruses, suggesting the activation of primary defense response to BNYVV infection in N. benthamiana. These results provide a global profile of miRNA changes induced by BNYVV infection and enhance our understanding of the mechanisms underlying BNYVV pathogenesis.

A Regulatory Network for miR156-SPL Module in Arabidopsis thaliana

Vegetative phase changes in plants describes the transition between juvenile and adult phases of vegetative growth before flowering. It is one of the most fundamental mechanisms for plants to sense developmental signals, presenting a complex process involving many still-unknown determinants. Several studies in annual and perennial plants have identified the conservative roles of miR156 and its targets, SBP/SPL genes, in guiding the switch of plant growth from juvenile to adult phases. Here, we review recent progress in understanding the regulation of miR156 expression and how miR156-SPLs mediated plant age affect other processes in Arabidopsis. Powerful high-throughput sequencing techniques have provided rich data to systematically study the regulatory mechanisms of miR156 regulation network. From this data, we draw an expanded miR156-regulated network that links plant developmental transition and other fundamental biological processes, gaining novel and broad insight into the molecular mechanisms of plant-age-related processes in Arabidopsis.

Next generation sequencing and proteomics in plant virology: how is Colombia doing?

Crop production and trade are two of the most economically important activities in Colombia, and viral diseases cause a high negative impact to agricultural sector. Therefore, the detection, diagnosis, control, and management of viral diseases are crucial. Currently, Next-Generation Sequencing (NGS) and ‘Omic’ technologies constitute a right-hand tool for the discovery of novel viruses and for studying virus-plant interactions. This knowledge allows the development of new viral diagnostic methods and the discovery of key components of infectious processes, which could be used to generate plants resistant to viral infections. Globally, crop sciences are advancing in this direction. In this review, advancements in ‘omic’ technologies and their different applications in plant virology in Colombia are discussed. In addition, bioinformatics pipelines and resources for omics data analyses are presented. Due to their decreasing prices, NGS technologies are becoming an affordable and promising means to explore many phytopathologies affecting a wide variety of Colombian crops so as to improve their trade potential.

Antiviral ARGONAUTEs Against Turnip Crinkle Virus Revealed by Image-Based Trait Analysis

RNA-based silencing functions as an important antiviral immunity mechanism in plants. Plant viruses evolved to encode viral suppressors of RNA silencing (VSRs) that interfere with the function of key components in the silencing pathway. As effectors in the RNA silencing pathway, ARGONAUTE (AGO) proteins are targeted by some VSRs, such as that encoded by Turnip crinkle virus (TCV). A VSR-deficient TCV mutant was used to identify AGO proteins with antiviral activities during infection. A quantitative phenotyping protocol using an image-based color trait analysis pipeline on the PlantCV platform, with temporal red, green, and blue imaging and a computational segmentation algorithm, was used to measure plant disease after TCV inoculation. This process captured and analyzed growth and leaf color of Arabidopsis (Arabidopsis thaliana) plants in response to virus infection over time. By combining this quantitative phenotypic data with molecular assays to detect local and systemic virus accumulation, AGO2, AGO3, and AGO7 were shown to play antiviral roles during TCV infection. In leaves, AGO2 and AGO7 functioned as prominent nonadditive, anti-TCV effectors, whereas AGO3 played a minor role. Other AGOs were required to protect inflorescence tissues against TCV. Overall, these results indicate that distinct AGO proteins have specialized, modular roles in antiviral defense across different tissues, and demonstrate the effectiveness of image-based phenotyping to quantify disease progression.

MAC3A and MAC3B, Two Core Subunits of the MOS4-Associated Complex, Positively Influence miRNA Biogenesis

MAC3A and MAC3B are conserved U-box-containing proteins in eukaryotes. They are subunits of the MOS4-associated complex (MAC) that plays essential roles in plant immunity and development in Arabidopsis thaliana However, their functional mechanisms remain elusive. Here, we show that Arabidopsis MAC3A and MAC3B act redundantly in microRNA (miRNA) biogenesis. Lack of both MAC3A and MAC3B in the mac3b mac3b double mutant reduces the accumulation of miRNAs, causing elevated transcript levels of miRNA targets. mac3a mac3b also decreases the levels of primary miRNA transcripts (pri-miRNAs). However, MAC3A and MAC3B do not affect the promoter activity of genes encoding miRNAs (MIR genes), suggesting that they may not affect MIR transcription. This result, together with the fact that MAC3A associates with pri-miRNAs in vivo, indicates that MAC3A and MAC3B may stabilize pri-miRNAs. Furthermore, we find that MAC3A and MAC3B interact with the DCL1 complex that catalyzes miRNA maturation, promote DCL1 activity, and are required for the localization of HYL1, a component of the DCL1 complex. Besides MAC3A and MAC3B, two other MAC subunits, CDC5 and PRL1, also function in miRNA biogenesis. Based on these results, we propose that MAC functions as a complex to control miRNA levels through modulating pri-miRNA transcription, processing, and stability.

Unveiling Chloroplast RNA Editing Events Using Next Generation Small RNA Sequencing Data

Organellar RNA editing involves the modification of nucleotide sequences to maintain conserved protein functions, mainly by reverting non-neutral codon mutations. The loss of plastid editing events, resulting from mutations in RNA editing factors or through stress interference, leads to developmental, physiological and photosynthetic alterations. Recently, next generation sequencing technology has generated the massive discovery of sRNA sequences and expanded the number of sRNA data. Here, we present a method to screen chloroplast RNA editing using public sRNA libraries from Arabidopsis, soybean and rice. We mapped the sRNAs against the nuclear, mitochondrial and plastid genomes to confirm predicted cytosine to uracil (C-to-U) editing events and identify new editing sites in plastids. Among the predicted editing sites, 40.57, 34.78, and 25.31% were confirmed using sRNAs from Arabidopsis, soybean and rice, respectively. SNP analysis revealed 58.2, 43.9, and 37.5% new C-to-U changes in the respective species and identified known and new putative adenosine to inosine (A-to-I) RNA editing in tRNAs. The present method and data reveal the potential of sRNA as a reliable source to identify new and confirm known editing sites.

Evaluating the toxicity of silicon dioxide nanoparticles on neural stem cells using RNA-Seq

Neural stem cells are characterized by self-renewal and multipotency, and a capacity to regenerate in response to brain injury or neurodegenerative disease.