Exogenous dsRNA Induces RNA Interference of a Chalcone Synthase Gene in Arabidopsis thaliana

Recent investigations have shown the possibility of artificial induction of RNA interference (RNAi) via plant foliar treatments with naked double-stranded RNA (dsRNA) to silence essential genes in plant fungal pathogens or to target viral RNAs. Furthermore, several studies have documented the downregulation of plant endogenous genes via external application of naked gene-specific dsRNAs and siRNAs to the plant surfaces. However, there are limited studies on the dsRNA processing and gene silencing mechanisms after external dsRNA application. Such studies would assist in the development of innovative tools for crop improvement and plant functional studies. In this study, we used exogenous gene-specific dsRNA to downregulate the gene of chalcone synthase (CHS), the key enzyme in the flavonoid/anthocyanin biosynthesis pathway, in Arabidopsis. The nonspecific NPTII-dsRNA encoding the nonrelated neomycin phosphotransferase II bacterial gene was used to treat plants in order to verify that any observed effects and processing of AtCHS mRNA were sequence specific. Using high-throughput small RNA (sRNA) sequencing, we obtained six sRNA-seq libraries for plants treated with water, AtCHS-dsRNA, or NPTII-dsRNA. After plant foliar treatments, we detected the emergence of a large number of AtCHS- and NPTII-encoding sRNAs, while there were no such sRNAs after control water treatment. Thus, the exogenous AtCHS-dsRNAs were processed into siRNAs and induced RNAi-mediated AtCHS gene silencing. The analysis showed that gene-specific sRNAs mapped to the AtCHS and NPTII genes unevenly with peak read counts at particular positions, involving primarily the sense strand, and documented a gradual decrease in read counts from 17-nt to 30-nt sRNAs. Results of the present study highlight a significant potential of exogenous dsRNAs as a promising strategy to induce RNAi-based downregulation of plant gene targets for plant management and gene functional studies.     

Small RNAs beyond Model Organisms: Have We Only Scratched the Surface?

Small RNAs (sRNAs) are essential regulators in the adaptation of bacteria to environmental changes and act by binding targeted mRNAs through base complementarity. Approximately 550 distinct families of sRNAs have been identified since their initial characterization in the 1980s, accelerated by the emergence of RNA-sequencing. Small RNAs are found in a wide range of bacterial phyla, but they are more prominent in highly researched model organisms compared to the rest of the sequenced bacteria. Indeed, Escherichia coli and Salmonella enterica contain the highest number of sRNAs, with 98 and 118, respectively, with Enterobacteriaceae encoding 145 distinct sRNAs, while other bacteria families have only seven sRNAs on average. Although the past years brought major advances in research on sRNAs, we have perhaps only scratched the surface, even more so considering RNA annotations trail behind gene annotations. A distinctive trend can be observed for genes, whereby their number increases with genome size, but this is not observable for RNAs, although they would be expected to follow the same trend. In this perspective, we aimed at establishing a more accurate representation of the occurrence of sRNAs in bacteria, emphasizing the potential for novel sRNA discoveries.     

Analysis of Small RNAs in Streptococcus mutans under Acid Stress—A New Insight for Caries Research

Streptococcus mutans (S. mutans) is the major clinical pathogen responsible for dental caries. Its acid tolerance has been identified as a significant virulence factor for its survival and cariogenicity in acidic conditions. Small RNAs (sRNAs) are recognized as key regulators of virulence and stress adaptation. Here, we constructed three libraries of sRNAs with small size exposed to acidic conditions for the first time, followed by verification using qRT-PCR. The levels of two sRNAs and target genes predicted to be bioinformatically related to acid tolerance were further evaluated under different acid stress conditions (pH 7.5, 6.5, 5.5, and 4.5) at three time points (0.5, 1, and 2 h). Meanwhile, bacterial growth characteristics and vitality were assessed. We obtained 1879 sRNAs with read counts of at least 100. One hundred and ten sRNAs were perfectly mapped to reported msRNAs in S. mutans. Ten out of 18 sRNAs were validated by qRT-PCR. The survival of bacteria declined as the acid was increased from pH 7.5 to 4.5 at each time point. The bacteria can proliferate under each pH except pH 4.5 with time. The levels of sRNAs gradually decreased from pH 7.5 to 5.5, and slightly increased in pH 4.5; however, the expression levels of target mRNAs were up-regulated in acidic conditions than in pH 7.5. These results indicate that some sRNAs are specially induced at acid stress conditions, involving acid adaptation, and provide a new insight into exploring the complex acid tolerance for S. mutans.     

Glucosinolate and Trichome Defenses in a Natural Arabidopsis lyrata Population

Glucosinolates (GS) and trichomes contribute to plant resistance against insect herbivores in the model Arabidopsis thaliana. The functional and genetic characteristics of herbivore defense, however, can differ even between closely related species. In a quantitative genetic experiment with the out-crossing perennial Arabidopsis lyrata spp. petraea, we measured constitutive GS composition, trichome density, leaf thickness, and plant resistance in four different herbivore interactions. In a single population of A. lyrata, we found heritable variation for trichome density as well as GS amount and carbon side-chain elongation ratios associated with activity in methylthioalkylmalate synthase (MAM). Unexpectedly, heritabilities for indole GS in A. lyrata were high and less affected by differences in plant age and environment than aliphatic GS. We found significant heritability in plant resistance to the specialist Plutella xylostella and generalist Trichoplusia ni, but not to the specialists Pieris brassicae and Phyllotreta cruciferae. Analyses of phenotypic and genetic correlations between candidate defense traits and insect resistance suggested that A. lyrata resistance was conferred by a combination of indole GS amount and trichome density, and, to a lesser extent, aliphatic GS ratios and leaf thickness. Variation in the most abundant compound, the aliphatic 3-hydroxypropyl GS, had little impact on A. lyrata herbivore resistance. The contribution of defense traits to resistance depended on the experimental herbivory context, and resistances were weakly correlated. A diversified defense strategy is likely to be important for long-lived individuals of A. lyrata that are subject to attack by many different herbivores in nature.     

Function Analysis of the PR55/B Gene Related to Self-Incompatibility in Chinese Cabbage Using CRISPR/Cas9

Chinese cabbage, a major crop in Korea, shows self-incompatibility (SI). SI is controlled by the type 2A serine/threonine protein phosphatases (PP2As). The PP2A gene is controlled by regulatory subunits that comprise a 36 kDa catalyst C subunit, a 65 kDa regulatory A subunit, and a variety of regulatory B subunits (50–70 kDa). Among them, the PP2A 55 kDa B regulatory subunit (PR55/B) gene located in the A05 chromosome has 13 exons spanning 2.9 kb, and two homologous genes, Bra018924 and Bra014296, were found to be present on the A06 and A08 chromosome, respectively. In this study, we performed a functional analysis of the PR55/B gene using clustered regularly interspaced short palindromic repeats/CRISPR-associated system 9 (CRISPR/Cas9)-mediated gene mutagenesis. CRISPR/Cas9 technology can be used to easily introduce mutations in the target gene. Tentative gene-edited lines were generated by the Agrobacterium-mediated transfer and were selected by PCR and Southern hybridization analysis. Furthermore, pods were confirmed to be formed in flower pollination (FP) as well as bud pollination (BP) in some gene-edited lines. Seed fertility of gene-edited lines indicated that the PR55/B gene plays a key role in SI. Finally, self-compatible T-DNA-free T₂ gene-edited plants and edited sequences of target genes were secured. The self-compatible Chinese cabbage developed in this study is expected to contribute to Chinese cabbage breeding.     

The QseEF Two-Component System-GlmY Small RNA Regulatory Pathway Controls Swarming in Uropathogenic Proteus mirabilis

Bacterial sensing of environmental signals through the two-component system (TCS) plays a key role in modulating virulence. In the search for the host hormone-sensing TCS, we identified a conserved qseEGF locus following glmY, a small RNA (sRNA) gene in uropathogenic Proteus mirabilis. Genes of glmY-qseE-qseG-qseF constitute an operon, and QseF binding sites were found in the glmY promoter region. Deletion of glmY or qseF resulted in reduced swarming motility and swarming-related phenotypes relative to the wild-type and the respective complemented strains. The qseF mutant had decreased glmYqseEGF promoter activity. Both glmY and qseF mutants exhibited decreased flhDC promoter activity and mRNA level, while increased rcsB mRNA level was observed in both mutants. Prediction by TargetRNA2 revealed cheA as the target of GlmY. Then, construction of the translational fusions containing various lengths of cheA 5′UTR for reporter assay and site-directed mutagenesis were performed to investigate the cheA-GlmY interaction in cheA activation. Notably, loss of glmY reduced the cheA mRNA level, and urea could inhibit swarming in a QseF-dependent manner. Altogether, this is the first report elucidating the underlying mechanisms for modulation of swarming motility by a QseEF-regulated sRNA GlmY, involving expression of cheA, rcsB and flhDC in uropathogenic P. mirabilis.     

miRNA Profiling and Its Role in Multi-Omics Regulatory Networks Connected with Somaclonal Variation in Cucumber (Cucumis sativus L.)

The role of miRNAs in connection with the phenomenon of somaclonal variation, which occurs during plant in vitro culture, remains uncertain. This study aims to investigate the possible role of miRNAs in multi-omics regulatory pathways in cucumber somaclonal lines. For this purpose, we performed sRNA sequencing (sRNA-seq) from cucumber fruit samples identified 8, 10 and 44 miRNAs that are differentially expressed between somaclones (S1, S2, S3 lines) and the reference B10 line of Cucumis sativus. For miRNA identification, we use ShortStack software designed to filter miRNAs from sRNAs according to specific program criteria. The identification of predicted in-silico targets revealed 2,886 mRNAs encoded by 644 genes. The functional annotation of miRNA’s target genes and gene ontology classification revealed their association with metabolic processes, response to stress, multicellular organism development, biosynthetic process and catalytic activity. We checked with bioinformatic analyses for possible interactions at the level of target proteins, differentially expressed genes (DEGs) and genes affected by genomic polymorphisms. We assume that miRNAs can indirectly influence molecular networks and play a role in many different regulatory pathways, leading to somaclonal variation. This regulation is supposed to occur through the process of the target gene cleavage or translation inhibition, which in turn affects the proteome, as we have shown in the example of molecular networks. This is a new approach combining levels from DNA-seq through mRNA-seq, sRNA-seq and in silico PPI in the area of plants’ somaclonal variation.     

Quantification of Species-Preferential Micropylar Chemoattraction in Arabidopsis by Fluorescein Diacetate Staining of Pollen Tubes

Sexual reproduction between males and females of the same species is essential for species maintenance. Ovular micropylar guidance, the last step of pollen tube guidance in angiosperms, contributes to species-preferential reproduction. Previous studies using semi-in vivo attraction assays showed that species-preferential attractant peptides are secreted from the ovule through its micropyle. However, conventional semi-in vivo assays usually depend on transgenic pollen tubes expressing a fluorescent protein to determine whether the tubes are attracted to the ovule to precisely penetrate the micropyle. Here, we found that fluorescein diacetate (FDA) staining was suitable for evaluating the micropylar guidance rate of non-transgenic pollen tubes in semi-in vivo conditions. Micropylar guidance was quantified for ovules and pollen tubes of Arabidopsis thaliana and Arabidopsis lyrata by combining FDA staining with modified semi-in vivo assays. Our results using the simple staining method showed that the ovules of each species secrete species-preferential attractants, and that pollen tubes respond more strongly to attractants of their own species compared with those of closely related species. LURE-type CRP810 attractant peptides were shown to be responsible for micropylar attraction of A. thaliana in the semi-in vivo assay. The POLLEN-SPECIFIC RECEPTOR-LIKE KINASE 6 (PRK6) receptor for LURE1, as well as an unidentified receptor for other LURE-type attractants, are involved in the species-preferential response of these two Arabidopsis species.     

The Multiverse of Plant Small RNAs: How Can We Explore It?

Plant small RNAs (sRNAs) are a heterogeneous group of noncoding RNAs with a length of 20–24 nucleotides that are widely studied due to their importance as major regulators in various biological processes. sRNAs are divided into two main classes—microRNAs (miRNAs) and small interfering RNAs (siRNAs)—which differ in their biogenesis and functional pathways. Their identification and enrichment with new structural variants would not be possible without the use of various high-throughput sequencing (NGS) techniques, allowing for the detection of the total population of sRNAs in plants. Classifying sRNAs and predicting their functional role based on such high-performance datasets is a nontrivial bioinformatics task, as plants can generate millions of sRNAs from a variety of biosynthetic pathways. Over the years, many computing tools have been developed to meet this challenge. Here, we review more than 35 tools developed specifically for plant sRNAs over the past few years and explore some of their basic algorithms for performing tasks related to predicting, identifying, categorizing, and quantifying individual sRNAs in plant samples, as well as visualizing the results of these analyzes. We believe that this review will be practical for biologists who want to analyze their plant sRNA datasets but are overwhelmed by the number of tools available, thus answering the basic question of how to choose the right one for a particular study.     

smartPARE: An R Package for Efficient Identification of True mRNA Cleavage Sites

Degradome sequencing is commonly used to generate high-throughput information on mRNA cleavage sites mediated by small RNAs (sRNA). In our datasets of potato (Solanum tuberosum, St) and Phytophthora infestans (Pi), initial predictions generated high numbers of cleavage site predictions, which highlighted the need of improved analytic tools. Here, we present an R package based on a deep learning convolutional neural network (CNN) in a machine learning environment to optimize discrimination of false from true cleavage sites. When applying smartPARE to our datasets on potato during the infection process by the late blight pathogen, 7.3% of all cleavage windows represented true cleavages distributed on 214 sites in P. infestans and 444 sites in potato. The sRNA landscape of the two organisms is complex with uneven sRNA production and cleavage regions widespread in the two genomes. Multiple targets and several cases of complex regulatory cascades, particularly in potato, was revealed. We conclude that our new analytic approach is useful for anyone working on complex biological systems and with the interest of identifying cleavage sites particularly inferred by sRNA classes beyond miRNAs.     

Identification and Evolution of Functional Alleles of the Previously Described Pollen Specific Myrosinase Pseudogene AtTGG6 in Arabidopsis thaliana

Myrosinases are β-thioglucoside glucohydrolases and serve as defense mechanisms against insect pests and pathogens by producing toxic compounds. AtTGG6 in Arabidopsis thaliana was previously reported to be a myrosinase pseudogene but specifically expressed in pollen. However, we found that AlTGG6, an ortholog to AtTGG6 in A. lyrata (an outcrossing relative of A. thaliana) was functional, suggesting that functional AtTGG6 alleles may still exist in A. thaliana. AtTGG6 alleles in 29 A. thaliana ecotypes were cloned and sequenced. Results indicate that ten alleles were functional and encoded Myr II type myrosinase of 512 amino acids, and myrosinase activity was confirmed by overexpressing AtTGG6 in Pichia pastoris. However, the 19 other ecotypes had disabled alleles with highly polymorphic frame-shift mutations and diversified sequences. Thirteen frame-shift mutation types were identified, which occurred independently many times in the evolutionary history within a few thousand years. The functional allele was expressed specifically in pollen similar to the disabled alleles but at a higher expression level, suggesting its role in defense of pollen against insect pests such as pollen beetles. However, the defense function may have become less critical after A. thaliana evolved to self-fertilization, and thus resulted in loss of function in most ecotypes.     

Overlapping Functions of the Paralogous Proteins AtPAP2 and AtPAP9 in Arabidopsis thaliana

Arabidopsis thaliana purple acid phosphatase 2 (AtPAP2), which is anchored to the outer membranes of chloroplasts and mitochondria, affects carbon metabolism by modulating the import of some preproteins into chloroplasts and mitochondria. AtPAP9 bears a 72% amino acid sequence identity with AtPAP2, and both proteins carry a hydrophobic motif at their C-termini. Here, we show that AtPAP9 is a tail-anchored protein targeted to the outer membrane of chloroplasts. Yeast two-hybrid and bimolecular fluorescence complementation experiments demonstrated that both AtPAP9 and AtPAP2 bind to a small subunit of rubisco 1B (AtSSU1B) and a number of chloroplast proteins. Chloroplast import assays using [³⁵S]-labeled AtSSU1B showed that like AtPAP2, AtPAP9 also plays a role in AtSSU1B import into chloroplasts. Based on these data, we propose that AtPAP9 and AtPAP2 perform overlapping roles in modulating the import of specific proteins into chloroplasts. Most plant genomes contain only one PAP-like sequence encoding a protein with a hydrophobic motif at the C-terminus. The presence of both AtPAP2 and AtPAP9 in the Arabidopsis genome may have arisen from genome duplication in Brassicaceae. Unlike AtPAP2 overexpression lines, the AtPAP9 overexpression lines did not exhibit early-bolting or high-seed-yield phenotypes. Their differential growth phenotypes could be due to the inability of AtPAP9 to be targeted to mitochondria, as the overexpression of AtPAP2 on mitochondria enhances the capacity of mitochondria to consume reducing equivalents.     

Conserved Opposite Functions in Plant Resistance to Biotrophic and Necrotrophic Pathogens of the Immune Regulator SRFR1

Plant immunity is mediated in large part by specific interactions between a host resistance protein and a pathogen effector protein, named effector-triggered immunity (ETI). ETI needs to be tightly controlled both positively and negatively to enable normal plant growth because constitutively activated defense responses are detrimental to the host. In previous work, we reported that mutations in SUPPRESSOR OF rps4-RLD1 (SRFR1), identified in a suppressor screen, reactivated EDS1-dependent ETI to Pseudomonas syringae pv. tomato (Pto) DC3000. Besides, mutations in SRFR1 boosted defense responses to the generalist chewing insect Spodoptera exigua and the sugar beet cyst nematode Heterodera schachtii. Here, we show that mutations in SRFR1 enhance susceptibility to the fungal necrotrophs Fusarium oxysporum f. sp. lycopersici (FOL) and Botrytis cinerea in Arabidopsis. To translate knowledge obtained in AtSRFR1 research to crops, we generated SlSRFR1 alleles in tomato using a CRISPR/Cas9 system. Interestingly, slsrfr1 mutants increased expression of SA-pathway defense genes and enhanced resistance to Pto DC3000. In contrast, slsrfr1 mutants elevated susceptibility to FOL. Together, these data suggest that SRFR1 is functionally conserved in both Arabidopsis and tomato and functions antagonistically as a negative regulator to (hemi-) biotrophic pathogens and a positive regulator to necrotrophic pathogens.