Distinct extremely abundant siRNAs associated with cosuppression in petunia

Functional evidence on the involvement of the MADS-box gene MdDAM4 in bud dormancy regulation in apple

Over the course of the year, temperate trees experience extremes in temperature and day length. In order to protect themselves from frost damage in winter, they enter a dormant state with no visible growth where all leaves are shed and buds are dormant. Also the young floral tissues need to withstand harsh winter conditions, as temperature fruit trees like apple develop their flower buds in the previous year of fruit development. So far, the genetic control of induction and release of dormancy is not fully understood. However, the transcription factor family of DORMANCY-Associated MADS-box (DAM) genes plays a major role in the control of winter dormancy. One of these genes is MdDAM4. This gene is expressed in the early phase of bud dormancy, but little is known about its function. Six transgenic apple lines were produced to study the function of MdDAM4 in apple. For plant transformation, the binary plasmid vector p9oN-35s-MdDAM4 was used that contains the coding sequence of MdDAM4 driven by the 35S promoter. Transgenicity of the lines was proven by PCR and southern hybridization. Based on siRNA sequencing and phenotypic observations, it was concluded that line M2024 overexpresses MdDAM4 whereas the gene is silenced in all other lines. Phenotyping of the transgenic lines provided evidence that the overexpression of MdDAM4 leads to an earlier induction and a later release of dormancy. Silencing this gene had exactly the opposite effects and thereby led to an increased duration of the vegetation period. Expression experiments revealed genes that were either potentially repressed or activated by MdDAM4. Among the potentially suppressed genes were several homologs of the cytokinin oxidase 5 (CKX5), five LOX homologs, and several expansins, which may indicate a link between MdDAM4 and the control of leaf senescence. Among the potentially activated genes is MdDAM1, which is in line with observed expression patterns during winter dormancy. MdDAM2, which shows little expression during endodormancy also appears to be activated by MdDAM4. Overall, this study provides experimental evidence with transgenic apple trees for MdDAM4 being an important regulator of the onset of bud dormancy in apple.

Contiguous identity between entire coding regions of transgenic and native genes rather than special regions is essential for a strong co-suppression

The discovery of co-suppression in plants has greatly boosted the study of gene silencing mechanisms, but its triggering mechanism has remained a mystery. In this study, we explored its possible trigger mechanism by using Fatty acid desaturase 2 (FAD2) and Fatty acid elongase 1 (FAE1) strong co-suppression systems. Analysis of small RNAs in FAD2 co-suppression lines showed that siRNAs distributed throughout the coding region of FAD2 with an accumulated peak. However, mutations of the peak siRNA-matched site and siRNA derived site had not alleviated the co-suppression of its transgenic lines. Synthetic FAD2 (AtFAD2sm), which has synonymous mutations in the entire coding region, failed to trigger any co-suppression. Furthermore, 5' and 3' portions of AtFAD2 and AtFAD2sm were swapped to form two hybrid genes, AtFAD2-3sm and AtFAD2-5sm. 80 % and 92 % of their transgenic lines exhibited co-suppression, respectively. Finally, FAE1s with different degrees of the continuous sequence identity compared with AtFAE1 were tested in their Arabidopsis transgenic lines, and the results showed the co-suppression frequency was reduced as their continuous sequence identity stepped down. This work suggests that contiguous identity between the entire coding regions of transgenic and native genes rather than a special region is essential for a strong co-suppression.

The unique dual targeting of AGO1 by two types of PRMT enzymes promotes phasiRNA loading in Arabidopsis thaliana

Arginine/R methylation (R-met) of proteins is a widespread post-translational modification (PTM), deposited by a family of protein arginine/R methyl transferase enzymes (PRMT). Regulations by R-met are involved in key biological processes deeply studied in metazoan. Among those, post-transcriptional gene silencing (PTGS) can be regulated by R-met in animals and in plants. It mainly contributes to safeguard processes as protection of genome integrity in germlines through the regulation of piRNA pathway in metazoan, or response to bacterial infection through the control of AGO2 in plants. So far, only PRMT5 has been identified as the AGO/PIWI R-met writer in higher eukaryotes. We uncovered that AGO1, the main PTGS effector regulating plant development, contains unique R-met features among the AGO/PIWI superfamily, and outstanding in eukaryotes. Indeed, AGO1 contains both symmetric (sDMA) and asymmetric (aDMA) R-dimethylations and is dually targeted by PRMT5 and by another type I PRMT in Arabidopsis thaliana. We showed also that loss of sDMA didn't compromise AtAGO1 subcellular trafficking in planta. Interestingly, we underscored that AtPRMT5 specifically promotes the loading of phasiRNA in AtAGO1. All our observations bring to consider this dual regulation of AtAGO1 in plant development and response to environment, and pinpoint the complexity of AGO1 post-translational regulation.

Frequency of cytosine methylation in the adjacent regions of soybean retrotransposon SORE-1 depends on chromosomal location

Mobilization of transposable elements (TEs) is suppressed by epigenetic mechanisms involving cytosine methylation. However, few studies have focused on clarifying relationships between epigenetic influences of TEs on the adjacent DNA regions and time after insertion of TEs into the genome and/or their chromosomal location. Here we addressed these issues using soybean retrotransposon SORE-1. We analyzed SORE-1, inserted in exon 1 of the GmphyA2 gene, one of the newest insertions in this family so far identified. Cytosine methylation was detected in this element but was barely present in the adjacent regions. These results were correlated, respectively, with the presence and absence of the production of short interfering RNAs. Cytosine methylation profiles of 74 SORE-1 elements in the Williams 82 reference genome indicated that methylation frequency in the adjacent regions of SORE-1 was profoundly higher in pericentromeric regions than in euchromatic chromosome arms and was only weakly correlated with the length of time after insertion into the genome. Notably, the higher level of methylation in the 5' adjacent regions of SORE-1 coincided with the presence of repetitive elements in pericentromeric regions. Together, these results suggest that epigenetic influence of SORE-1 on the adjacent regions is influenced by its location on the chromosome.

Phased secondary small interfering RNAs in Camellia sinensis var. assamica

Phased secondary small interfering RNAs (phasiRNAs) in plants play important roles in regulating genome stability, plant development and stress adaption. Camellia sinensis var. assamica has immense economic, medicinal and cultural significance. However, there are still no studies of phasiRNAs and their putative functions in this valuable plant. We identified 476 and 43 PHAS loci which generated 4290 twenty one nucleotide (nt) and 264 twenty four nt phasiRNAs, respectively. Moreover, the analysis of degradome revealed more than 35000 potential targets for these phasiRNAs. We identified several conserved 21 nt phasiRNA generation pathways in tea plant, including miR390 → TAS3, miR482/miR2118 → NB-LRR, miR393 → F-box, miR828 → MYB/TAS4, and miR7122 → PPR in this study. Furthermore, we found that some transposase and plant mobile domain genes could generate phasiRNAs. Our results show that phasiRNAs target genes in the same family in cis- or trans-manners, and different members of the same gene family may generate the same phasiRNAs. The phasiRNAs, generated by transposase and plant mobile domain genes, and their targets, suggest that phasiRNAs may be involved in the inhibition of transposable elements in tea plant. To summarize, these results provide a comprehensive view of phasiRNAs in Camellia sinensis var. assamica.

Indian Himalayan natural Arabidopsis thaliana accessions with abolished miR158 levels exhibit robust miR173-initiated trans-acting cascade silencing

Small RNAs (sRNAs) such as microRNAs (miRNAs) and small interfering RNAs (siRNAs) are short 20-24-nucleotide non-coding RNAs. They are key regulators of gene expression in plants and other organisms. Several 22-nucleotide miRNAs trigger biogenesis cascades of trans-acting secondary siRNAs, which are involved in various developmental and stress responses. Here we show that Himalayan Arabidopsis thaliana accessions having natural mutations in the miR158 locus exhibit robust cascade silencing of the pentatricopeptide repeat (PPR)-like locus. Furthermore, we show that these cascade sRNAs trigger tertiary silencing of a gene involved in transpiration and stomatal opening. The natural deletions or insertions in MIR158 led to improper processing of miR158 precursors, thereby blocking synthesis of mature miR158. Reduced miR158 levels led to increased levels of its target, a pseudo-PPR gene that is targeted by tasiRNAs generated by the miR173 cascade in other accessions. Using sRNA datasets derived from Indian Himalayan accessions, as well as overexpression and knockout lines of miR158, we show that absence of miR158 led to buildup of pseudo-PPR-derived tertiary sRNAs. These tertiary sRNAs mediated robust silencing of a gene involved in stomatal closure in Himalayan accessions lacking miR158 expression. We functionally validated the tertiary phasiRNA that targets NHX2, which encodes a Na⁺ -K⁺ /H⁺ antiporter protein, thereby regulating transpiration and stomatal conductance. Overall, we report the role of the miRNA-TAS-siRNA-pseudogene-tertiary phasiRNA-NHX2 pathway in plant adaptation.

The N-terminal extension of Arabidopsis ARGONAUTE 1 is essential for microRNA activities

microRNAs (miRNAs) regulate target gene expression through their ARGONAUTE (AGO) effector protein, mainly AGO1 in Arabidopsis thaliana. In addition to the highly conserved N, PAZ, MID and PIWI domains with known roles in RNA silencing, AGO1 contains a long, unstructured N-terminal extension (NTE) of little-known function. Here, we show that the NTE is indispensable for the functions of Arabidopsis AGO1, as a lack of the NTE leads to seedling lethality. Within the NTE, the region containing amino acids (a.a.) 91 to 189 is essential for rescuing an ago1 null mutant. Through global analyses of small RNAs, AGO1-associated small RNAs, and miRNA target gene expression, we show that the region containing a.a. 91-189 is required for the loading of miRNAs into AGO1. Moreover, we show that reduced nuclear partitioning of AGO1 did not affect its profiles of miRNA and ta-siRNA association. Furthermore, we show that the 1-to-90a.a. and 91-to-189a.a. regions of the NTE redundantly promote the activities of AGO1 in the biogenesis of trans-acting siRNAs. Together, we report novel roles of the NTE of Arabidopsis AGO1.

An efficient deep learning based predictor for identifying miRNA-triggered phasiRNA loci in plant

Phasic small interfering RNAs are plant secondary small interference RNAs that typically generated by the convergence of miRNAs and polyadenylated mRNAs. A growing number of studies have shown that miRNA-initiated phasiRNA plays crucial roles in regulating plant growth and stress responses. Experimental verification of miRNA-initiated phasiRNA loci may take considerable time, energy and labor. Therefore, computational methods capable of processing high throughput data have been proposed one by one. In this work, we proposed a predictor (DIGITAL) for identifying miRNA-initiated phasiRNAs in plant, which combined a multi-scale residual network with a bi-directional long-short term memory network. The negative dataset was constructed based on positive data, through replacing 60% of nucleotides randomly in each positive sample. Our predictor achieved the accuracy of 98.48% and 94.02% respectively on two independent test datasets with different sequence length. These independent testing results indicate the effectiveness of our model. Furthermore, DIGITAL is of robustness and generalization ability, and thus can be easily extended and applied for miRNA target recognition of other species. We provide the source code of DIGITAL, which is freely available at https://github.com/yuanyuanbu/DIGITAL.

Transcriptome-wide analyses of RNA m6A methylation in hexaploid wheat reveal its roles in mRNA translation regulation

N6-methyladenosine (m6A) is the most abundant RNA modification in eukaryotic messenger RNAs. m6A was discovered in wheat about 40 years ago; however, its potential roles in wheat remain unknown. In this study, we profiled m6As in spikelets transcriptome at the flowering stage of hexaploid wheat and found that m6As are evenly distributed across the A, B, and D subgenomes but their extents and locations vary across homeologous genes. m6As are enriched in homeologous genes with close expression levels and the m6A methylated genes are more conserved. The extent of m6A methylation is negatively correlated with mRNA expression levels and its presence on mRNAs has profound impacts on mRNA translation in a location-dependent manner. Specifically, m6As within coding sequences and 3'UTRs repress the translation of mRNAs while the m6As within 5'UTRs and start codons could promote it. The m6A-containing mRNAs are significantly enriched in processes and pathways of "translation" and "RNA transport," suggesting the potential role of m6As in regulating the translation of genes involved in translation regulation. Our data also show a stronger translation inhibition by small RNAs (miRNA and phasiRNA) than by m6A methylation, and no synergistical effect between the two was observed. We propose a secondary amplification machinery of translation regulation triggered by the changes in m6A methylation status. Taken together, our results suggest translation regulation as a key role played by m6As in hexaploid wheat.

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

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

Post-transcriptional gene silencing of the chalcone synthase gene CHS causes corolla lobe-specific whiting of Japanese gentian

Post-transcriptional gene silencing of the chalcone synthase gene CHS specifically suppresses anthocyanin biosynthesis in corolla lobes and is responsible for the formation of a stripe type bicolor in Japanese gentian. The flower of Japanese gentian is a bell-shaped corolla composed of lobes and plicae, which is painted uniformly blue. However, the gentian cultivar 'Hakuju' shows bicolor phenotype (blue-white stripe corolla), in which anthocyanin accumulation is suppressed only in corolla lobes. Expression analysis indicated that steady-state levels of chalcone synthase (CHS) transcripts were remarkably reduced in corolla lobes compared with plicae during petal pigmentation initiation. However, no significant difference in expression levels of other flavonoid biosynthetic structural and regulatory genes was detected in its lobes and plicae. On feeding naringenin in white lobes, anthocyanin accumulation was recovered. Northern blotting probed with CHS confirmed the abundant accumulation of small RNAs in corolla lobes. Likewise, small RNA-seq analysis indicated that short reads from its lobes were predominantly mapped onto the 2nd exon region of the CHS gene, whereas those from the plicae were scarcely mapped. Subsequent infection with the gentian ovary ringspot virus (GORV), which had an RNA-silencing activity, showed the recovery of partial pigmentation in lobes. Hence, these results strongly suggested that suppressing anthocyanin accumulation in the lobes of bicolored 'Hakuju' was attributed to the specific degradation of CHS mRNA in corolla lobes, which was through post-transcriptional gene silencing (PTGS). Herein, we revealed the molecular mechanism of strip bicolor formation in Japanese gentian, and showed that PTGS of CHS was also responsible for flower color pattern in a floricultural plant other than petunia and dahlia.

Post-transcriptional gene silencing of CYP76AD controls betalain biosynthesis in bracts of bougainvillea

Betalain is one of four major plant pigments and shares some features with anthocyanin; however, no plant has been found to biosynthesize both pigments. Previous studies have reported that anthocyanin biosynthesis in some plants is regulated by post-transcriptional gene-silencing (PTGS), but the importance of PTGS in betalain biosynthesis remains unclear. In this study, we report the occurrence of PTGS in betalain biosynthesis in bougainvillea (Bougainvillea peruviana) 'Thimma', which produces bracts of three different color on the same plant, namely pink, white, and pink-white. This resembles the unstable anthocyanin pigmentation phenotype that is associated with PTGS, and hence we anticipated the presence of PTGS in the betalain biosynthetic pathway. To test this, we analysed pigments, gene expression, small RNAs, and transient overexpression. Our results demonstrated that PTGS of BpCYP76AD1, a gene encoding one of the betalain biosynthesis enzymes, is responsible for the loss of betalain biosynthesis in 'Thimma'. Neither the genetic background nor DNA methylation in the BpCYP76AD1 sequence could explain the induction of PTGS, implying that another locus controls the unstable pigmentation. Our results indicate that naturally occurring PTGS contributes to the diversification of color patterns not only in anthocyanin biosynthesis but also in betalain biosynthesis.

Laccase Directed Lignification Is One of the Major Processes Associated With the Defense Response Against Pythium ultimum Infection in Apple Roots

Apple replant disease (ARD), incited by a pathogen complex including Pythium ultimum, causes stunted growth or death of newly planted trees at replant sites. Development and deployment of resistant or tolerant rootstocks offers a cost-effective, ecologically friendly, and durable approach for ARD management. Maximized exploitation of natural resistance requires integrated efforts to identify key regulatory mechanisms underlying resistance traits in apple. In this study, miRNA profiling and degradome sequencing identified major miRNA pathways and candidate genes using six apple rootstock genotypes with contrasting phenotypes to P. ultimum infection. The comprehensive RNA-seq dataset offered an expansive view of post-transcriptional regulation of apple root defense activation in response to infection from P. ultimum. Several pairs of miRNA families and their corresponding targets were identified for their roles in defense response in apple roots, including miR397-laccase, miR398-superoxide dismutase, miR10986-polyphenol oxidase, miR482-resistance genes, and miR160-auxin response factor. Of these families, the genotype-specific expression patterns of miR397 indicated its fundamental role in developing defense response patterns to P. ultimum infection. Combined with other identified copper proteins, the importance of cellular fortification, such as lignification of root tissues by the action of laccase, may critically contribute to genotype-specific resistance traits. Our findings suggest that quick and enhanced lignification of apple roots may significantly impede pathogen penetration and minimize the disruption of effective defense activation in roots of resistant genotypes. The identified target miRNA species and target genes consist of a valuable resource for subsequent functional analysis of their roles during interaction between apple roots and P. ultimum.

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

Background

In plants, microRNAs (miRNAs) are pivotal regulators of plant development and stress responses. Different computational tools and web servers have been developed for plant miRNA target prediction; however, in silico prediction normally contains false positive results. In addition, many plant miRNA target prediction servers lack information for miRNA-triggered phased small interfering RNAs (phasiRNAs). Creating a comprehensive and relatively high-confidence plant miRNA target database is much needed.

Results

Here, we report TarDB, an online database that collects three categories of relatively high-confidence plant miRNA targets: (i) cross-species conserved miRNA targets; (ii) degradome/PARE (Parallel Analysis of RNA Ends) sequencing supported miRNA targets; (iii) miRNA-triggered phasiRNA loci. TarDB provides a user-friendly interface that enables users to easily search, browse and retrieve miRNA targets and miRNA initiated phasiRNAs in a broad variety of plants. TarDB has a comprehensive collection of reliable plant miRNA targets containing previously unreported miRNA targets and miRNA-triggered phasiRNAs even in the well-studied model species. Most of these novel miRNA targets are relevant to lineage-specific or species-specific miRNAs. TarDB data is freely available at http://www.biosequencing.cn/TarDB.

Conclusions

In summary, TarDB serves as a useful web resource for exploring relatively high-confidence miRNA targets and miRNA-triggered phasiRNAs in plants.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07680-5.

Perspective: 50 years of plant chromosome biology

The past 50 years has been the greatest era of plant science discovery, and most of the discoveries have emerged from or been facilitated by our knowledge of plant chromosomes. At last we have descriptive and mechanistic outlines of the information in chromosomes that programs plant life. We had almost no such information 50 years ago when few had isolated DNA from any plant species. The important features of genes have been revealed through whole genome comparative genomics and testing of variants using transgenesis. Progress has been enabled by the development of technologies that had to be invented and then become widely available. Arabidopsis (Arabidopsis thaliana) and rice (Oryza sativa) have played extraordinary roles as model species. Unexpected evolutionary dramas were uncovered when learning that chromosomes have to manage constantly the vast numbers of potentially mutagenic families of transposons and other repeated sequences. The chromatin-based transcriptional and epigenetic mechanisms that co-evolved to manage the evolutionary drama as well as gene expression and 3-D nuclear architecture have been elucidated these past 20 years. This perspective traces some of the major developments with which I have become particularly familiar while seeking ways to improve crop plants. I draw some conclusions from this look-back over 50 years during which the scientific community has (i) exposed how chromosomes guard, readout, control, recombine, and transmit information that programs plant species, large and small, weed and crop, and (ii) modified the information in chromosomes for the purposes of genetic, physiological, and developmental analyses and plant improvement.

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

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

Widespread occurrence of microRNA-mediated target cleavage on membrane-bound polysomes

Background

Small RNAs (sRNAs) including microRNAs (miRNAs) and small interfering RNAs (siRNAs) serve as core players in gene silencing at transcriptional and post-transcriptional levels in plants, but their subcellular localization has not yet been well studied, thus limiting our mechanistic understanding of sRNA action.

Results

We investigate the cytoplasmic partitioning of sRNAs and their targets globally in maize (Zea mays, inbred line “B73”) and rice (Oryza sativa, cv. “Nipponbare”) by high-throughput sequencing of polysome-associated sRNAs and 3′ cleavage fragments, and find that both miRNAs and a subset of 21-nucleotide (nt)/22-nt siRNAs are enriched on membrane-bound polysomes (MBPs) relative to total polysomes (TPs) across different tissues. Most of the siRNAs are generated from transposable elements (TEs), and retrotransposons positively contributed to MBP overaccumulation of 22-nt TE-derived siRNAs (TE-siRNAs) as opposed to DNA transposons. Widespread occurrence of miRNA-mediated target cleavage is observed on MBPs, and a large proportion of these cleavage events are MBP-unique. Reproductive 21PHAS (21-nt phasiRNA-generating) and 24PHAS (24-nt phasiRNA-generating) precursors, which were commonly considered as noncoding RNAs, are bound by polysomes, and high-frequency cleavage of 21PHAS precursors by miR2118 and 24PHAS precursors by miR2275 is further detected on MBPs. Reproductive 21-nt phasiRNAs are enriched on MBPs as opposed to TPs, whereas 24-nt phasiRNAs are nearly completely devoid of polysome occupancy.

Conclusions

MBP overaccumulation is a conserved pattern for cytoplasmic partitioning of sRNAs, and endoplasmic reticulum (ER)-bound ribosomes function as an independent regulatory layer for miRNA-induced gene silencing and reproductive phasiRNA biosynthesis in maize and rice.

Role of phasiRNAs from two distinct phasing frames of GhMYB2 loci in cis- gene regulation in the cotton genome

BACKGROUND

Phased small interfering RNA (phasiRNA) is primarily derived from the 22-nt miRNA targeting loci. GhMYB2, a gene with potential roles in cotton fiber cell fate determination, is a target gene of miR828 and miR858 in the generation of phasiRNAs.

RESULTS

In the presented work, through the evaluation of phasing scores and phasiRNA distribution pattern, we found that phasiRNAs from GhMYB2 were derived from the 3' cleavage fragments of 22-nt miR828 and 21-nt miR858 respectively. These two miRNA targeting sites initiated two phasing frames on transcripts of one locus. By means of RNA ligase-mediated rapid amplification of cDNA ends (RLM-RACE), we further demonstrated that phasiRNAs derived from the two phasing frames played a role in cis-regulation of GhMYB2. The phasiRNAs derived from GhMYB2 were expressed in the somatic tissues, especially in anther and hypocotyl. We further employed our previous small RNA sequencing data as well as the degradome data of cotton fiber bearing ovules, anthers, hypocotyls and embryogenic calli tissues published in public databases, to validate the expression, phasing pattern and functions of phasiRNAs.

CONCLUSIONS

The presenting research provide insights of the molecular mechanism of phasiRNAs in regulation of GhMYB2 loci.

The MADS-Box Gene MdDAM1 Controls Growth Cessation and Bud Dormancy in Apple

Apple trees require a long exposure to chilling temperature during winter to acquire competency to flower and grow in the following spring. Climate change or adverse meteorological conditions can impair release of dormancy and delay bud break, hence jeopardizing fruit production and causing substantial economic losses. In order to characterize the molecular mechanisms controlling bud dormancy in apple we focused our work on the MADS-box transcription factor gene MdDAM1. We show that MdDAM1 silencing is required for the release of dormancy and bud break in spring. MdDAM1 transcript levels are drastically reduced in the low-chill varieties 'Anna' and 'Dorsett Golden' compared to 'Golden Delicious' corroborating its role as a key genetic factor controlling the release of bud dormancy in Malus species. The functional characterization of MdDAM1 using RNA silencing resulted in trees unable to cease growth in winter and that displayed an evergrowing, or evergreen, phenotype several years after transgenesis. These trees lost their capacity to enter in dormancy and produced leaves and shoots regardless of the season. A transcriptome study revealed that apple evergrowing lines are a genocopy of 'Golden Delicious' trees at the onset of the bud break with the significant gene repression of the related MADS-box gene MdDAM4 as a major feature. We provide the first functional evidence that MADS-box transcriptional factors are key regulators of bud dormancy in pome fruit trees and demonstrate that their silencing results in a defect of growth cessation in autumn. Our findings will help producing low-chill apple variants from the elite commercial cultivars that will withstand climate change.

Next-Generation Sequence Databases: RNA and Genomic Informatics Resources for Plants

We developed public web sites and resources for data access, display, and analysis of plant small RNAs. These web sites are interconnected with related data types. The current generation of these informatics tools was developed for Illumina data, evolving over more than 15 years of improvements. Our online databases have customized web interfaces to uniquely handle and display RNA-derived data from diverse plant species, ranging from Arabidopsis (Arabidopsis thaliana) to wheat (Triticum spp.), including many crop and model species. The web interface displays the abundance and genomic context of data for small RNAs, parallel analysis of RNA ends/degradome reads, RNA sequencing, and even chromatin immunoprecipitation sequencing data; it also provides information about potentially novel transcripts (antisense transcripts, alternative splice isoforms, and regulatory intergenic transcripts). Numerous options are included for downloading data as tables or via web services. Interpretation of these data is facilitated by the inclusion of extensive repeat or transposon data in our genome viewer. We have developed graphical and analytical tools, including a new viewer and a query page for the analysis of phased small RNAs; these are particularly useful for understanding the complex small RNA pathways of plants. These public databases are accessible at https://mpss.danforthcenter.org.

FIERY1 promotes microRNA accumulation by suppressing rRNA-derived small interfering RNAs in Arabidopsis

Plant microRNAs (miRNAs) associate with ARGONAUTE1 (AGO1) to direct post-transcriptional gene silencing and regulate numerous biological processes. Although AGO1 predominantly binds miRNAs in vivo, it also associates with endogenous small interfering RNAs (siRNAs). It is unclear whether the miRNA/siRNA balance affects miRNA activities. Here we report that FIERY1 (FRY1), which is involved in 5'-3' RNA degradation, regulates miRNA abundance and function by suppressing the biogenesis of ribosomal RNA-derived siRNAs (risiRNAs). In mutants of FRY1 and the nuclear 5'-3' exonuclease genes XRN2 and XRN3, we find that a large number of 21-nt risiRNAs are generated through an endogenous siRNA biogenesis pathway. The production of risiRNAs correlates with pre-rRNA processing defects in these mutants. We also show that these risiRNAs are loaded into AGO1, causing reduced loading of miRNAs. This study reveals a previously unknown link between rRNA processing and miRNA accumulation.

Small RNA discovery in the interaction between barley and the powdery mildew pathogen

Background

Plants encounter pathogenic and non-pathogenic microorganisms on a nearly constant basis. Small RNAs such as siRNAs and miRNAs/milRNAs influence pathogen virulence and host defense responses. We exploited the biotrophic interaction between the powdery mildew fungus, Blumeria graminis f. sp. hordei (Bgh), and its diploid host plant, barley (Hordeum vulgare) to explore fungal and plant sRNAs expressed during Bgh infection of barley leaf epidermal cells.

Results

RNA was isolated from four fast-neutron immune-signaling mutants and their progenitor over a time course representing key stages of Bgh infection, including appressorium formation, penetration of epidermal cells, and development of haustorial feeding structures. The Cereal Introduction (CI) 16151 progenitor carries the resistance allele Mla6, while Bgh isolate 5874 harbors the AVR_a6 avirulence effector, resulting in an incompatible interaction. Parallel Analysis of RNA Ends (PARE) was used to verify sRNAs with likely transcript targets in both barley and Bgh. Bgh sRNAs are predicted to regulate effectors, metabolic genes, and translation-related genes. Barley sRNAs are predicted to influence the accumulation of transcripts that encode auxin response factors, NAC transcription factors, homeodomain transcription factors, and several splicing factors. We also identified phasing small interfering RNAs (phasiRNAs) in barley that overlap transcripts that encode receptor-like kinases (RLKs) and nucleotide-binding, leucine-rich domain proteins (NLRs).

Conclusions

These data suggest that Bgh sRNAs regulate gene expression in metabolism, translation-related, and pathogen effectors. PARE-validated targets of predicted Bgh milRNAs include both EKA (effectors homologous to AVR_k1 and AVR_a10) and CSEP (candidate secreted effector protein) families. We also identified barley phasiRNAs and miRNAs in response to Bgh infection. These include phasiRNA loci that overlap with a significant proportion of receptor-like kinases, suggesting an additional sRNA control mechanism may be active in barley leaves as opposed to predominant R-gene phasiRNA overlap in many eudicots. In addition, we identified conserved miRNAs, novel miRNA candidates, and barley genome mapped sRNAs that have PARE validated transcript targets in barley. The miRNA target transcripts are enriched in transcription factors, signaling-related proteins, and photosynthesis-related proteins. Together these results suggest both barley and Bgh control metabolism and infection-related responses via the specific accumulation and targeting of genes via sRNAs.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5947-z) contains supplementary material, which is available to authorized users.

Analysis of microRNAs, phased small interfering RNAs and their potential targets in Rosarugosa Thunb

BACKGROUND

MicroRNAs (miRNAs) are small non-coding RNAs that play important roles by regulating other genes. Rosa rugosa Thunb. is an important ornamental and edible plant, yet there are only a few studies on the miRNAs and their functions in R. rugosa.

RESULTS

We sequenced 10 samll RNA profiles from the roots, petals, pollens, stamens, and leaves and 4 RNA-seq profiles in leaves and petals to analysis miRNA, phasiRNAs and mRNAs in R. rugosa. In addition, we acquired a degradome sequencing profile from leaf of R. rugosa to identify miRNA and phasiRNA targets using the SeqTar algorithm. We have identified 321 conserved miRNA homologs including primary transcripts for 25 conserved miRNAs, and 22 novel miRNAs. We identified 592 putative targets of the conserved miRNAs or tasiRNAs that showed significant accumulations of degradome reads. We found differential expression patterns of conserved miRNAs in five different tissues of R. rugosa. We identified three hundred and thirty nine 21 nucleotide (nt) PHAS loci, and forty nine 24 nt PHAS loci, respectively. Our results suggest that miR482 triggers generations of phasiRNAs by targeting nucleotide-binding, leucine-rich repeat (NB-LRR) disease resistance genes in R. rugosa. Our results also suggest that the deregulated genes in leaves and petals are significantly enriched in GO terms and KEGG pathways related to metabolic processes and photosynthesis.

CONCLUSIONS

These results significantly enhanced our knowledge of the miRNAs and phasiRNAs, as well as their potential functions, in R. rugosa.

The disease resistance protein SNC1 represses the biogenesis of microRNAs and phased siRNAs

Plants evolved an array of disease resistance genes (R genes) to fight pathogens. In the absence of pathogen infection, NBS-LRR genes, which comprise a major subfamily of R genes, are suppressed by a small RNA cascade involving microRNAs (miRNAs) that trigger the biogenesis of phased siRNAs (phasiRNAs) from R gene transcripts. However, whether or how R genes influence small RNA biogenesis is unknown. In this study, we isolate a mutant with global defects in the biogenesis of miRNAs and phasiRNAs in Arabidopsis thaliana and trace the defects to the over accumulation and nuclear localization of an R protein SNC1. We show that nuclear SNC1 represses the transcription of miRNA and phasiRNA loci, probably through the transcriptional corepressor TPR1. Intriguingly, nuclear SNC1 reduces the accumulation of phasiRNAs from three source R genes and concomitantly, the expression of a majority of the ~170R genes is up-regulated. Taken together, this study suggests an R gene-miRNA-phasiRNA regulatory module that amplifies plant immune responses.

A small RNA-based signaling cascade prevents the induction of plant resistance genes (R-genes) in the absence of pathogen challenge. Here Cai et al. show that nuclear accumulation of the R protein SNC1 can activate immunity by suppressing small RNA production and releasing R-gene repression.

Genome-wide identification and comprehensive analysis of microRNAs and phased small interfering RNAs in watermelon

BACKGROUND

MicroRNAs (miRNAs) are a class of endogenous small non-coding RNAs involved in the post-transcriptional gene regulation and play a critical role in plant growth, development and stress responses. Watermelon (Citrullus lanatus L.) is one of the important agricultural crops worldwide. However, the watermelon miRNAs and phasiRNAs and their functions are not well explored.

RESULTS

Here we carried out computational and experimental analysis of miRNAs and phased small interfering RNAs (phasiRNAs) in watermelon by analyzing 14 small RNA profiles from roots, leaves, androecium, petals, and fruits, and one published small RNA profile of mixed tissues. To identify the targets of miRNAs and phasiRNAs, we generated a degradome profile for watermelon leaf which is analyzed using the SeqTar algorithm. We identified 97 conserved pre-miRNAs, of which 58 have not been reported previously and 348 conserved mature miRNAs without precursors. We also found 9 novel pre-miRNAs encoding 18 mature miRNAs. One hundred and one 21 nucleotide (nt) PHAS loci, and two hundred and forty one 24 nt PHAS loci were also identified. We identified 127 conserved targets of the conserved miRNAs and TAS3-derived tasiRNAs by analyzing a degradome profile of watermelon leaf.

CONCLUSIONS

The presented results provide a comprehensive view of small regulatory RNAs and their targets in watermelon.

An intragenic approach to confer glyphosate resistance in chile (Capsicum annuum) by introducing an in vitro mutagenized chile EPSPS gene encoding for a glyphosate resistant EPSPS protein

Chile pepper (Capsicum annuum) is an important high valued crop worldwide, and when grown on a large scale has problems with weeds. One important herbicide used is glyphosate. Glyphosate inactivates the enzyme 5-enolpyruvylshikimate-3-phosphate synthase (EPSPS), a key enzyme in the synthesis of aromatic amino acids. A transgenic approach towards making glyphosate resistant plants, entails introducing copies of a gene encoding for glyphosate-resistant EPSPS enzyme into the plant. The main objective of our work was to use an intragenic approach to confer resistance to glyphosate in chile which would require using only chile genes for transformation including the selectable marker. Tobacco was used as the transgenic system to identify different gene constructs that would allow for the development of the intragenic system for chile, since chile transformation is inefficient. An EPSPS gene was isolated from chile and mutagenized to introduce substitutions that are known to make the encoded enzyme resistant to glyphosate. The promoter for EPSPS gene was isolated from chile and the mutagenized chile EPSPS cDNA was engineered behind both the CaMV35S promoter and the EPSPS promoter. The leaves from the transformants were checked for resistance to glyphosate using a cut leaf assay. In tobacco, though both gene constructs exhibited some degree of resistance to glyphosate, the construct with the CaMV35S promoter was more effective and as such chile was transformed with this gene construct. The chile transformants showed resistance to low concentrations of glyphosate. Furthermore, preliminary studies showed that the mutated EPSPS gene driven by the CaMV35S promoter could be used as a selectable marker for transformation. We have shown that an intragenic approach can be used to confer glyphosate-resistance in chile. However, we need a stronger chile promoter and a mutated chile gene that encodes for a more glyphosate resistant EPSPS protein.

Phased secondary small interfering RNAs in Panaxnotoginseng

Background

Recent results demonstrated that either non-coding or coding genes generate phased secondary small interfering RNAs (phasiRNAs) guided by specific miRNAs. Till now, there is no studies for phasiRNAs in Panax notoginseng (Burk.) F.H. Chen (P. notoginseng), an important traditional Chinese herbal medicinal plant species.

Methods

Here we performed a genome-wide discovery of phasiRNAs and its host PHAS loci in P. notoginseng by analyzing small RNA sequencing profiles. Degradome sequencing profile was used to identify the trigger miRNAs of these phasiRNAs and potential targets of phasiRNAs. We also used RLM 5’-RACE to validate some of the identified phasiRNA targets.

Results

After analyzing 24 small RNA sequencing profiles of P. notoginseng, 204 and 90 PHAS loci that encoded 21 and 24 nucleotide (nt) phasiRNAs, respectively, were identified. Furthermore, we found that phasiRNAs produced from some pentatricopeptide repeat-contain (PPR) genes target another layer of PPR genes as validated by both the degradome sequencing profile and RLM 5’-RACE analysis. We also found that miR171 with 21 nt triggers the generations of 21 nt phasiRNAs from its conserved targets.

Conclusions

We validated that some phasiRNAs generated from PPRs and TASL genes are functional by targeting other PPRs in trans. These results provide the first set of PHAS loci and phasiRNAs in P. notoginseng, and enhance our understanding of PHAS in plants.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-4331-0) contains supplementary material, which is available to authorized users.

Small RNAs, emerging regulators critical for the development of horticultural traits

Small RNAs (sRNAs) have been recently recognized as key genetic and epigenetic regulators in various organisms, ranging from the modification of DNA and histone methylations to the modulation of the abundance of coding or non-coding RNAs. In plants, major regulatory sRNAs are classified as respective microRNA (miRNA) and small interfering RNA (siRNA) species, with the former primarily engaging in posttranscriptional regulation while the latter in transcriptional one. Many of these characterized sRNAs are involved in regulation of diverse biological programs, processes, and pathways in response to developmental cues, environmental signals/stresses, pathogen infection, and pest attacks. Recently, sRNAs-mediated regulations have also been extensively investigated in horticultural plants, with many novel mechanisms unveiled, which display far more mechanistic complexity and unique regulatory features compared to those studied in model species. Here, we review the recent progress of sRNA research in horticultural plants, with emphasis on mechanistic aspects as well as their relevance to trait regulation. Given that major and pioneered sRNA research has been carried out in the model and other plants, we also discuss ongoing sRNA research on these plants. Because miRNAs and phased siRNAs (phasiRNAs) are the most studied sRNA regulators, this review focuses on their biogenesis, conservation, function, and targeted genes and traits as well as the mechanistic relation between them, aiming at providing readers comprehensive information instrumental for future sRNA research in horticulture crops.

miR3954 is a trigger of phasiRNAs that affects flowering time in citrus

In plant, a few 22-nt miRNAs direct cleavages of their targets and trigger the biogenesis of phased small interfering RNAs (phasiRNAs) in plant. In this study, we characterized a miRNA triggering phasiRNAs generation, miR3954, and explored its downstream target genes and potential function. Our results demonstrated that miR3954 showed specific expression in the flowers of citrus species, and it targeted a NAC transcription factor (Cs7 g22460) and two non-coding RNA transcripts (lncRNAs, Cs1 g09600 and Cs1 g09635). The production of phasiRNAs was detected from transcripts targeted by miR3954, and was further verified in both sequencing data and transient expression experiments. PhasiRNAs derived from the two lncRNAs targeted not only miR3954-targeted NAC gene but also additional NAC homologous genes. No homologous genes of these two lncRNAs were found in plants other than citrus species, implying that this miR3954-lncRNAs-phasiRNAs-NAC pathway is likely citrus-specific. Transgenic analysis indicated that the miR3954-overexpressing lines showed decreased transcripts of lncRNA, elevated abundance of phasiRNAs and reduced expression of NAC genes. Interestingly, the overexpression of miR3954 leads to early flowering in citrus plants. In summary, our results illustrated a model of the regulatory network of miR3954-lncRNA-phasiRNAs-NAC, which may be functionally involved in flowering in citrus.

Comprehensive transcriptome analyses reveal tomato plant responses to tobacco rattle virus-based gene silencing vectors

In plants, virus-induced gene silencing (VIGS) is a popular tool for functional genomic studies or rapidly assessing individual gene functions. However, molecular details regarding plant responses to viral vectors remain elusive, which may complicate experimental designs and data interpretation. To this end, we documented whole transcriptome changes of tomato elicited by the application of the most widely used tobacco rattle virus (TRV)-based vectors, using comprehensive genome-wide analyses. Our data illustrated multiple biological processes with functional implications, including (1) the enhanced activity of miR167 in guiding the cleavage of an auxin response factor; (2) reduced accumulation of phased secondary small interfering RNAs from two genomic loci; (3) altered expression of ~500 protein-coding transcripts; and (4) twenty long noncoding RNAs specifically responsive to TRV vectors. Importantly, we unraveled large-scale changes in mRNA alternative splicing patterns. These observations will facilitate future application of VIGS vectors for functional studies benefiting the plant research community and help deepen the understanding of plant-virus interactions.

Small RNA profiles from Panax notoginseng roots differing in sizes reveal correlation between miR156 abundances and root biomass levels

Plant genomes encode several classes of small regulatory RNAs (sRNAs) that play critical roles in both development and stress responses. Panax notoginseng (Burk.) F.H. Chen (P. notoginseng) is an important traditional Chinese herbal medicinal plant species for its haemostatic effects. Therefore, the root yield of P. notoginseng is a major economically important trait since the roots of P. notoginseng are the parts used to produce medicine. To identify sRNAs that are critical for the root biomass of P. notoginseng, we performed a comprehensive study of miRNA transcriptomes from P. notoginseng roots of different biomasses. We identified 675 conserved miRNAs, of which 180 pre-miRNAs are also identified, and three TAS3 loci in P. notoginseng. By using degradome sequencing, we identified 79 conserved miRNA:target or tasiRNA:target interactions, of which eight were further confirmed with the RLM 5'-RACE experiments. More importantly, our results revealed that a member of miR156 family and one of its SPL target genes have inverse expression levels, which is tightly correlated with greater root biomass contents. These results not only contributes to overall understanding of post-transcriptional gene regulation in roots of P. notoginseng but also could serve as markers for breeding P. notoginseng with greater root yield.

Comprehensive Transcriptome Analyses Reveal that Potato Spindle Tuber Viroid Triggers Genome-Wide Changes in Alternative Splicing, Inducible trans-Acting Activity of Phased Secondary Small Interfering RNAs, and Immune Responses

Many pathogens express noncoding RNAs (ncRNAs) during infection processes. In the most extreme case, pathogenic ncRNAs alone (such as viroids) can infect eukaryotic organisms, leading to diseases. While a few pathogenic ncRNAs have been implicated in regulating gene expression, the functions of most pathogenic ncRNAs in host-pathogen interactions remain unclear. Here, we employ potato spindle tuber viroid (PSTVd) infecting tomato as a system to dissect host interactions with pathogenic ncRNAs, using comprehensive transcriptome analyses. We uncover various new activities in regulating gene expression during PSTVd infection, such as genome-wide alteration in alternative splicing of host protein-coding genes, enhanced guided-cleavage activities of a host microRNA, and induction of the trans-acting function of phased secondary small interfering RNAs. Furthermore, we reveal that PSTVd infection massively activates genes involved in plant immune responses, mainly those in the calcium-dependent protein kinase and mitogen-activated protein kinase cascades, as well as prominent genes involved in hypersensitive responses, cell wall fortification, and hormone signaling. Intriguingly, our data support a notion that plant immune systems can respond to pathogenic ncRNAs, which has broad implications for providing new opportunities for understanding the complexity of immune systems in differentiating "self" and "nonself," as well as lay the foundation for resolving the long-standing question regarding the pathogenesis mechanisms of viroids and perhaps other infectious RNAs.IMPORTANCE Numerous pathogens, including viruses, express pathogenic noncoding transcripts during infection. In the most extreme case, pathogenic noncoding RNAs alone (i.e., viroids) can cause disease in plants. While some work has demonstrated that pathogenic noncoding RNAs interact with host factors for function, the biological significance of pathogenic noncoding RNAs in host-pathogen interactions remains largely unclear. Here, we apply comprehensive genome-wide analyses of plant-viroid interactions and discover several novel molecular activities underlying nuclear-replicating viroid infection processes in plants, including effects on the expression and function of host noncoding transcripts, as well as the alternative splicing of host protein-coding genes. Importantly, we show that plant immunity is activated upon infection of a nuclear-replicating viroid, which is a new concept that helps to understand viroid-based pathogenesis. Our finding has broad implications for understanding the complexity of host immune systems and the diverse functions of noncoding RNAs.

Mutations in Argonaute5 Illuminate Epistatic Interactions of the K1 and I Loci Leading to Saddle Seed Color Patterns in Glycine max

The soybean (Glycine max) seed coat has distinctive, genetically programmed patterns of pigmentation, and the recessive k1 mutation can epistatically overcome the dominant I and ii alleles, which inhibit seed color by producing small interfering RNAs (siRNAs) targeting chalcone synthase (CHS) mRNAs. Small RNA sequencing of dissected regions of immature seed coats demonstrated that CHS siRNA levels cause the patterns produced by the ii and ik alleles of the I locus, which restrict pigment to the hilum or saddle region of the seed coat, respectively. To identify the K1 locus, we compared RNA-seq data from dissected regions of two Clark isolines having similar saddle phenotypes mediated by CHS siRNAs but different genotypes (homozygous ik K1 versus homozygous ii k1). By examining differentially expressed genes, mapping information, and genome resequencing, we identified a 129-bp deletion in Glyma.11G190900 encoding Argonaute5 (AGO5), a member of the Argonaute family. Amplicon sequencing of several independent saddle pattern mutants from different genetic backgrounds revealed independent lesions affecting AGO5, thus establishing Glyma.11G190900 as the K1 locus. Nonfunctional AGO5 from k1 alleles leads to altered distributions of CHS siRNAs, thus explaining how the k1 mutation reverses the phenotype of the seed coat regions from yellow to pigmented, even in the presence of the normally dominant I or ii alleles.

Small RNAs Derived from the T-DNA of Agrobacterium rhizogenes in Hairy Roots of Phaseolus vulgaris

Agrobacterium rhizogenes is a pathogenic bacteria that causes hairy root disease by transferring bacterial DNA into the plant genome. It is an essential tool for industry and research due to its capacity to produce genetically modified roots and whole organisms. Here, we identified and characterized small RNAs generated from the transfer DNA (T-DNA) of A. rhizogenes in hairy roots of common bean (Phaseolus vulgaris). Distinct abundant A. rhizogenes T-DNA-derived small RNAs (ArT-sRNAs) belonging to several oncogenes were detected in hairy roots using high-throughput sequencing. The most abundant and diverse species of ArT-sRNAs were those of 21- and 22-nucleotides in length. Many T-DNA encoded genes constituted phasiRNA producing loci (PHAS loci). Interestingly, degradome analysis revealed that ArT-sRNAs potentially target genes of P. vulgaris. In addition, we detected low levels of ArT-sRNAs in the A. rhizogenes-induced calli generated at the wound site before hairy root emergence. These results suggest that RNA silencing targets several genes from T-DNA of A. rhizogenes in hairy roots of common bean. Therefore, the role of RNA silencing observed in this study has implications in our understanding and usage of this unique plant-bacteria interaction.

Biogenesis of phased siRNAs on membrane-bound polysomes in Arabidopsis

Small RNAs are central players in RNA silencing, yet their cytoplasmic compartmentalization and the effects it may have on their activities have not been studied at the genomic scale. Here we report that Arabidopsis microRNAs (miRNAs) and small interfering RNAs (siRNAs) are distinctly partitioned between the endoplasmic reticulum (ER) and cytosol. All miRNAs are associated with membrane-bound polysomes (MBPs) as opposed to polysomes in general. The MBP association is functionally linked to a deeply conserved and tightly regulated activity of miRNAs - production of phased siRNAs (phasiRNAs) from select target RNAs. The phasiRNA precursor RNAs, thought to be noncoding, are on MBPs and are occupied by ribosomes in a manner that supports miRNA-triggered phasiRNA production, suggesting that ribosomes on the rough ER impact siRNA biogenesis. This study reveals global patterns of cytoplasmic partitioning of small RNAs and expands the known functions of ribosomes and ER.

Identification of Arbuscular Mycorrhiza (AM)-Responsive microRNAs in Tomato

A majority of land plants can form symbiosis with arbuscular mycorrhizal (AM) fungi. MicroRNAs (miRNAs) have been implicated to regulate this process in legumes, but their involvement in non-legume species is largely unknown. In this study, by performing deep sequencing of sRNA libraries in tomato roots and comparing with tomato genome, a total of 700 potential miRNAs were predicted, among them, 187 are known plant miRNAs that have been previously deposited in miRBase. Unlike the profiles in other plants such as rice and Arabidopsis, a large proportion of predicted tomato miRNAs was 24 nt in length. A similar pattern was observed in the potato genome but not in tobacco, indicating a Solanum genus-specific expansion of 24-nt miRNAs. About 40% identified tomato miRNAs showed significantly altered expressions upon Rhizophagus irregularis inoculation, suggesting the potential roles of these novel miRNAs in AM symbiosis. The differential expression of five known and six novel miRNAs were further validated using qPCR analysis. Interestingly, three up-regulated known tomato miRNAs belong to a known miR171 family, a member of which has been reported in Medicago truncatula to regulate AM symbiosis. Thus, the miR171 family likely regulates AM symbiosis conservatively across different plant lineages. More than 1000 genes targeted by potential AM-responsive miRNAs were provided and their roles in AM symbiosis are worth further exploring.

Genome-Wide Mapping of Targets of Maize Histone Deacetylase HDA101 Reveals Its Function and Regulatory Mechanism during Seed Development

Histone deacetylases (HDACs) regulate histone acetylation levels by removing the acetyl group from lysine residues. The maize (Zea mays) HDACHDA101 influences several aspects of development, including kernel size; however, the molecular mechanism by which HDA101 affects kernel development remains unknown. In this study, we find that HDA101 regulates the expression of transfer cell-specific genes, suggesting that their misregulation may be associated with the defects in differentiation of endosperm transfer cells and smaller kernels observed in hda101 mutants. To investigate HDA101 function during the early stages of seed development, we performed genome-wide mapping of HDA101 binding sites. We observed that, like mammalian HDACs, HDA101 mainly targets highly and intermediately expressed genes. Although loss of HDA101 can induce histone hyperacetylation of its direct targets, this often does not involve variation in transcript levels. A small subset of inactive genes that must be negatively regulated during kernel development is also targeted by HDA101 and its loss leads to hyperacetylation and increased expression of these inactive genes. Finally, we report that HDA101 interacts with members of different chromatin remodeling complexes, such as NFC103/MSI1 and SNL1/SIN3-like protein corepressors. Taken together, our results reveal a complex genetic network regulated by HDA101 during seed development and provide insight into the different mechanisms of HDA101-mediated regulation of transcriptionally active and inactive genes.

Investigation of a miRNA-Induced Gene Silencing Technique in Petunia Reveals Alterations in miR173 Precursor Processing and the Accumulation of Secondary siRNAs from Endogenous Genes

MIGS (miRNA-induced gene silencing) is a straightforward and efficient gene silencing technique in Arabidopsis. It works by exploiting miR173 to trigger the production of phasiRNAs (phased small interfering RNAs). MIGS can be used in plant species other than Arabidopsis by co-expression of miR173 and target gene fragments fused to an upstream miR173 target site. However, the efficiency and technical mechanisms have not been thoroughly investigated in other plants. In this work, two vectors, pMIGS-chs and pMIGS-pds, were constructed and transformed into petunia plants. The transgenic plants showed CHS (chalcone synthase) and PDS (phytoene desaturase) gene-silencing phenotypes respectively, indicating that MIGS functions in petunia. MIGS-chs plants were used to investigate the mechanisms of this technique in petunia. Results of 5′- RACE showed that the miR173 target site was cleaved at the expected position and that endogenous CHS genes were cut at multiple positions. Small RNA deep sequencing analysis showed that the processing of Arabidopsis miR173 precursors in MIGS-chs transgenic petunia plants did not occur in exactly the same way as in Arabidopsis, suggesting differences in the machinery of miRNA processing between plant species. Small RNAs in-phase with the miR173 cleavage register were produced immediately downstream from the cleavage site and out-of-phase small RNAs were accumulated at relatively high levels from processing cycle 5 onwards. Secondary siRNAs were generated from multiple sites of endogenous CHS-A and CHS-J genes, indicating that miR173 cleavage induced siRNAs have the same ability to initiate siRNA transitivity as the siRNAs functioning in co-suppression and hpRNA silencing. On account of the simplicity of vector construction and the transitive amplification of signals from endogenous transcripts, MIGS is a good alternative gene silencing method for plants, especially for silencing a cluster of homologous genes with redundant functions.

Small RNAs in Plant Responses to Abiotic Stresses: Regulatory Roles and Study Methods

To survive under abiotic stresses in the environment, plants trigger a reprogramming of gene expression, by transcriptional regulation or translational regulation, to turn on protective mechanisms. The current focus of research on how plants cope with abiotic stresses has transitioned from transcriptomic analyses to small RNA investigations. In this review, we have summarized and evaluated the current methodologies used in the identification and validation of small RNAs and their targets, in the context of plant responses to abiotic stresses.

Post-transcriptional gene silencing triggered by sense transgenes involves uncapped antisense RNA and differs from silencing intentionally triggered by antisense transgenes

Although post-transcriptional gene silencing (PTGS) has been studied for more than a decade, there is still a gap in our understanding of how de novo silencing is initiated against genetic elements that are not supposed to produce double-stranded (ds)RNA. Given the pervasive transcription occurring throughout eukaryote genomes, we tested the hypothesis that unintended transcription could produce antisense (as)RNA molecules that participate to the initiation of PTGS triggered by sense transgenes (S-PTGS). Our results reveal a higher level of asRNA in Arabidopsis thaliana lines that spontaneously trigger S-PTGS than in lines that do not. However, PTGS triggered by antisense transgenes (AS-PTGS) differs from S-PTGS. In particular, a hypomorphic ago1 mutation that suppresses S-PTGS prevents the degradation of asRNA but not sense RNA during AS-PTGS, suggesting a different treatment of coding and non-coding RNA by AGO1, likely because of AGO1 association to polysomes. Moreover, the intended asRNA produced during AS-PTGS is capped whereas the asRNA produced during S-PTGS derives from 3′ maturation of a read-through transcript and is uncapped. Thus, we propose that uncapped asRNA corresponds to the aberrant RNA molecule that is converted to dsRNA by RNA-DEPENDENT RNA POLYMERASE 6 in siRNA-bodies to initiate S-PTGS, whereas capped asRNA must anneal with sense RNA to produce dsRNA that initiate AS-PTGS.

Analysis of the siRNA-Mediated Gene Silencing Process Targeting Three Homologous Genes Controlling Soybean Seed Oil Quality

In the past decade, RNA silencing has gained significant attention because of its success in genomic scale research and also in the genetic improvement of crop plants. However, little is known about the molecular basis of siRNA processing in association with its target transcript. To reveal this process for improving hpRNA-mediated gene silencing in crop plants, the soybean GmFAD3 gene family was chosen as a test model. We analyzed RNAi mutant soybean lines in which three members of the GmFAD3 gene family were silenced. The silencing levels of FAD3A, FAD3B and FAD3C were correlated with the degrees of sequence homology between the inverted repeat of hpRNA and the GmFAD3 transcripts in the RNAi lines. Strikingly, transgenes in two of the three RNAi lines were heavily methylated, leading to a dramatic reduction of hpRNA-derived siRNAs. Small RNAs corresponding to the loop portion of the hairpin transcript were detected while much lower levels of siRNAs were found outside of the target region. siRNAs generated from the 318-bp inverted repeat were found to be diced much more frequently at stem sequences close to the loop and associated with the inferred cleavage sites on the target transcripts, manifesting "hot spots". The top candidate hpRNA-derived siRNA share certain sequence features with mature miRNA. This is the first comprehensive and detailed study revealing the siRNA-mediated gene silencing mechanism in crop plants using gene family GmFAD3 as a test model.

Genome-wide identification of the Phaseolus vulgaris sRNAome using small RNA and degradome sequencing

BACKGROUND

MiRNAs and phasiRNAs are negative regulators of gene expression. These small RNAs have been extensively studied in plant model species but only 10 mature microRNAs are present in miRBase version 21, the most used miRNA database, and no phasiRNAs have been identified for the model legume Phaseolus vulgaris. Thanks to the recent availability of the first version of the common bean genome, degradome data and small RNA libraries, we are able to present here a catalog of the microRNAs and phasiRNAs for this organism and, particularly, we suggest new protagonists in the symbiotic nodulation events.

RESULTS

We identified a set of 185 mature miRNAs, including 121 previously unpublished sequences, encoded by 307 precursors and distributed in 98 families. Degradome data allowed us to identify a total of 181 targets for these miRNAs. We reveal two regulatory networks involving conserved miRNAs: those known to play crucial roles in the establishment of nodules, and novel miRNAs present only in common bean, suggesting a specific role for these sequences. In addition, we identified 125 loci that potentially produce phased small RNAs, with 47 of them having all the characteristics of being triggered by a total of 31 miRNAs, including 14 new miRNAs identified in this study.

CONCLUSIONS

We provide here a set of new small RNAs that contribute to the broader knowledge of the sRNAome of Phaseolus vulgaris. Thanks to the identification of the miRNA targets from degradome analysis and the construction of regulatory networks between the mature microRNAs, we present here the probable functional regulation associated with the sRNAome and, particularly, in N2-fixing symbiotic nodules.

A dynamic evolutionary and functional landscape of plant phased small interfering RNAs

Background

Secondary, phased small interfering RNAs (phasiRNAs) derived from protein-coding or noncoding loci (PHAS) are emerging as a new type of regulators of gene expression in plants. However, the evolution and function of these novel siRNAs in plant species remain largely unexplored.

Results

We systematically analyzed PHAS loci in 23 plant species covering major phylogenetic groups spanning alga, moss, gymnosperm, basal angiosperm, monocot, and dicot. We identified over 3,300 PHAS loci, among which ~1,600 were protein-coding genes. Most of these PHAS loci were novel and clade- or species-specific and showed distinct expression patterns in association with particular development stages, viral infection, or abiotic stresses. Unexpectedly, numerous PHAS loci produced phasiRNAs from introns or exon–intron junction regions. Our comprehensive analysis suggests that phasiRNAs predominantly regulate protein-coding genes from which they are derived and genes from the same families of the phasiRNA-deriving genes, in contrast to the dominant trans-regulatory mode of miRNAs. The stochastic occurrence of many PHAS loci in the plant kingdom suggests their young evolutionary origins.

Conclusions

Our study discovered an unprecedented diversity of protein-coding genes that produce phasiRNAs in a wide variety of plants, and set a kingdom-wide foundation for investigating the novel roles of phasiRNAs in shaping phenotype diversities of plants.

Electronic supplementary material

The online version of this article (doi:10.1186/s12915-015-0142-4) contains supplementary material, which is available to authorized users.

Enhanced formation of aromatic amino acids increases fragrance without affecting flower longevity or pigmentation in Petunia × hybrida

Purple Petunia × hybrida V26 plants accumulate fragrant benzenoid-phenylpropanoid molecules and anthocyanin pigments in their petals. These specialized metabolites are synthesized mainly from the aromatic amino acids phenylalanine. Here, we studied the profile of secondary metabolites of petunia plants, expressing a feedback-insensitive bacterial form of 3-deoxy-di-arabino-heptulosonate 7-phosphate synthase enzyme (AroG*) of the shikimate pathway, as a tool to stimulate the conversion of primary to secondary metabolism via the aromatic amino acids. We focused on specialized metabolites contributing to flower showy traits. The presence of AroG* protein led to increased aromatic amino acid levels in the leaves and high phenylalanine levels in the petals. In addition, the AroG* petals accumulated significantly higher levels of fragrant benzenoid-phenylpropanoid volatiles, without affecting the flowers' lifetime. In contrast, AroG* abundance had no effect on flavonoids and anthocyanins levels. The metabolic profile of all five AroG* lines was comparable, even though two lines produced the transgene in the leaves, but not in the petals. This implies that phenylalanine produced in leaves can be transported through the stem to the flowers and serve as a precursor for formation of fragrant metabolites. Dipping cut petunia stems in labelled phenylalanine solution resulted in production of labelled fragrant volatiles in the flowers. This study emphasizes further the potential of this metabolic engineering approach to stimulate the production of specialized metabolites and enhance the quality of various plant organs. Furthermore, transformation of vegetative tissues with AroG* is sufficient for induced production of specialized metabolites in organs such as the flowers.

Genome-wide analysis of leafbladeless1-regulated and phased small RNAs underscores the importance of the TAS3 ta-siRNA pathway to maize development

Maize leafbladeless1 (lbl1) encodes a key component in the trans-acting short-interfering RNA (ta-siRNA) biogenesis pathway. Correlated with a great diversity in ta-siRNAs and the targets they regulate, the phenotypes conditioned by mutants perturbing this small RNA pathway vary extensively across species. Mutations in lbl1 result in severe developmental defects, giving rise to plants with radial, abaxialized leaves. To investigate the basis for this phenotype, we compared the small RNA content between wild-type and lbl1 seedling apices. We show that LBL1 affects the accumulation of small RNAs in all major classes, and reveal unexpected crosstalk between ta-siRNA biogenesis and other small RNA pathways regulating transposons. Interestingly, in contrast to data from other plant species, we found no evidence for the existence of phased siRNAs generated via the one-hit model. Our analysis identified nine TAS loci, all belonging to the conserved TAS3 family. Information from RNA deep sequencing and PARE analyses identified the tasiR-ARFs as the major functional ta-siRNAs in the maize vegetative apex where they regulate expression of AUXIN RESPONSE FACTOR3 (ARF3) homologs. Plants expressing a tasiR-ARF insensitive arf3a transgene recapitulate the phenotype of lbl1, providing direct evidence that deregulation of ARF3 transcription factors underlies the developmental defects of maize ta-siRNA biogenesis mutants. The phenotypes of Arabidopsis and Medicago ta-siRNA mutants, while strikingly different, likewise result from misexpression of the tasiR-ARF target ARF3. Our data indicate that diversity in TAS pathways and their targets cannot fully account for the phenotypic differences conditioned by ta-siRNA biogenesis mutants across plant species. Instead, we propose that divergence in the gene networks downstream of the ARF3 transcription factors or the spatiotemporal pattern during leaf development in which these proteins act constitute key factors underlying the distinct contributions of the ta-siRNA pathway to development in maize, Arabidopsis, and possibly other plant species as well.

Mutations in maize leafbladeless1 (lbl1) that disrupt ta-siRNA biogenesis give rise to plants with thread-like leaves that have lost top/bottom polarity. We used genomic approaches to identify lbl1-dependent small RNAs and their targets to determine the basis for these polarity defects. This revealed substantial diversity in small RNA pathways across plant species and identified unexpected roles for LBL1 in the regulation of repetitive elements within the maize genome. We further show that only ta-siRNA loci belonging to the TAS3 family function in the maize vegetative apex. The TAS3-derived tasiR-ARFs are the main ta-siRNA active in the apex, and misregulation of their ARF3 targets emerges as the basis for the lbl1 leaf polarity defects. Supporting this, we show that plants expressing arf3a transcripts insensitive to tasiR-ARF-directed cleavage recapitulate the phenotypes observed in lbl1. The TAS3 ta-siRNA pathway, including the regulation of ARF3 genes, is conserved throughout land plant evolution, yet the phenotypes of plants defective for ta-siRNA biogenesis are strikingly different. Our data leads us to propose that divergence in the processes regulated by the ARF3 transcription factors or the spatiotemporal pattern during development in which these proteins act, underlies the diverse developmental contributions of this small RNA pathway across plants.

Genome-wide discovery and analysis of phased small interfering RNAs in Chinese sacred lotus

Phased small interfering RNA (phasiRNA) generating loci (briefly as PHAS) in plants are a novel class of genes that are normally regulated by microRNAs (miRNAs). Similar to miRNAs, phasiRNAs encoded by PHAS play important regulatory roles by targeting protein coding transcripts in plant species. We performed a genome-wide discovery of PHAS loci in Chinese sacred lotus and identified a total of 106 PHAS loci. Of these, 47 loci generate 21 nucleotide (nt) phasiRNAs and 59 loci generate 24 nt phasiRNAs, respectively. We have also identified a new putative TAS3 and a putative TAS4 loci in the lotus genome. Our results show that some of the nucleotide-binding, leucine-rich repeat (NB-LRR) disease resistance proteins and MYB transcription factors potentially generate phasiRNAs. Furthermore, our results suggest that some large subunit (LSU) rRNAs can derive putative phasiRNAs, which is potentially resulted from crosstalk between small RNA biogenesis pathways that are employed to process rRNAs and PHAS loci, respectively. Some of the identified phasiRNAs have putative trans-targets with less than 4 mismatches, suggesting that the identified PHAS are involved in many different pathways. Finally, the discovery of 24 nt PHAS in lotus suggests that there are 24 nt PHAS in dicots.

Distinctive profiles of small RNA couple inverted repeat-induced post-transcriptional gene silencing with endogenous RNA silencing pathways in Arabidopsis

The experimental induction of RNA silencing in plants often involves expression of transgenes encoding inverted repeat (IR) sequences to produce abundant dsRNAs that are processed into small RNAs (sRNAs). These sRNAs are key mediators of post-transcriptional gene silencing (PTGS) and determine its specificity. Despite its application in agriculture and broad utility in plant research, the mechanism of IR-PTGS is incompletely understood. We generated four sets of 60 Arabidopsis plants, each containing IR transgenes expressing different configurations of uidA and CHALCONE Synthase (At-CHS) gene fragments. Levels of PTGS were found to depend on the orientation and position of the fragment in the IR construct. Deep sequencing and mapping of sRNAs to corresponding transgene-derived and endogenous transcripts identified distinctive patterns of differential sRNA accumulation that revealed similarities among sRNAs associated with IR-PTGS and endogenous sRNAs linked to uncapped mRNA decay. Detailed analyses of poly-A cleavage products from At-CHS mRNA confirmed this hypothesis. We also found unexpected associations between sRNA accumulation and the presence of predicted open reading frames in the trigger sequence. In addition, strong IR-PTGS affected the prevalence of endogenous sRNAs, which has implications for the use of PTGS for experimental or applied purposes.

Roles of small RNAs in soybean defense against Phytophthora sojae infection

The genus Phytophthora consists of many notorious pathogens of crops and forestry trees. At present, battling Phytophthora diseases is challenging due to a lack of understanding of their pathogenesis. We investigated the role of small RNAs in regulating soybean defense in response to infection by Phytophthora sojae, the second most destructive pathogen of soybean. Small RNAs, including microRNAs (miRNAs) and small interfering RNAs (siRNAs), are universal regulators that repress target gene expression in eukaryotes. We identified known and novel small RNAs that differentially accumulated during P. sojae infection in soybean roots. Among them, miR393 and miR166 were induced by heat-inactivated P. sojae hyphae, indicating that they may be involved in soybean basal defense. Indeed, knocking down the level of mature miR393 led to enhanced susceptibility of soybean to P. sojae; furthermore, the expression of isoflavonoid biosynthetic genes was drastically reduced in miR393 knockdown roots. These data suggest that miR393 promotes soybean defense against P. sojae. In addition to miRNAs, P. sojae infection also resulted in increased accumulation of phased siRNAs (phasiRNAs) that are predominantly generated from canonical resistance genes encoding nucleotide binding-leucine rich repeat proteins and genes encoding pentatricopeptide repeat-containing proteins. This work identifies specific miRNAs and phasiRNAs that regulate defense-associated genes in soybean during Phytophthora infection.

Distinct and concurrent pathways of Pol II- and Pol IV-dependent siRNA biogenesis at a repetitive trans-silencer locus in Arabidopsis thaliana

Short interfering RNAs (siRNAs) homologous to transcriptional regulatory regions can induce RNA-directed DNA methylation (RdDM) and transcriptional gene silencing (TGS) of target genes. In our system, siRNAs are produced by transcribing an inverted DNA repeat (IR) of enhancer sequences, yielding a hairpin RNA that is processed by several Dicer activities into siRNAs of 21-24 nt. Primarily 24-nt siRNAs trigger RdDM of the target enhancer in trans and TGS of a downstream GFP reporter gene. We analyzed siRNA accumulation from two different structural forms of a trans-silencer locus in which tandem repeats are embedded in the enhancer IR and distinguished distinct RNA polymerase II (Pol II)- and Pol IV-dependent pathways of siRNA biogenesis. At the original silencer locus, Pol-II transcription of the IR from a 35S promoter produces a hairpin RNA that is diced into abundant siRNAs of 21-24 nt. A silencer variant lacking the 35S promoter revealed a normally masked Pol IV-dependent pathway that produces low levels of 24-nt siRNAs from the tandem repeats. Both pathways operate concurrently at the original silencer locus. siRNAs accrue only from specific regions of the enhancer and embedded tandem repeat. Analysis of these sequences and endogenous tandem repeats producing siRNAs revealed the preferential accumulation of siRNAs at GC-rich regions containing methylated CG dinucleotides. In addition to supporting a correlation between base composition, DNA methylation and siRNA accumulation, our results highlight the complexity of siRNA biogenesis at repetitive loci and show that Pol II and Pol IV use different promoters to transcribe the same template.

Undesired small RNAs originate from an artificial microRNA precursor in transgenic petunia (Petunia hybrida)

Although artificial microRNA (amiRNA) technology has been used frequently in gene silencing in plants, little research has been devoted to investigating the accuracy of amiRNA precursor processing. In this work, amiRNAchs1 (amiRchs1), based on the Arabidopsis miR319a precursor, was expressed in order to suppress the expression of CHS genes in petunia. The transgenic plants showed the CHS gene-silencing phenotype. A modified 5′ RACE technique was used to map small-RNA-directed cleavage sites and to detect processing intermediates of the amiRchs1 precursor. The results showed that the target CHS mRNAs were cut at the expected sites and that the amiRchs1 precursor was processed from loop to base. The accumulation of small RNAs in amiRchs1 transgenic petunia petals was analyzed using the deep-sequencing technique. The results showed that, alongside the accumulation of the desired artificial microRNAs, additional small RNAs that originated from other regions of the amiRNA precursor were also accumulated at high frequency. Some of these had previously been found to be accumulated at low frequency in the products of ath-miR319a precursor processing and some of them were accompanied by 3′-tailing variant. Potential targets of the undesired small RNAs were discovered in petunia and other Solanaceae plants. The findings draw attention to the potential occurrence of undesired target silencing induced by such additional small RNAs when amiRNA technology is used. No appreciable production of secondary small RNAs occurred, despite the fact that amiRchs1 was designed to have perfect complementarity to its CHS-J target. This confirmed that perfect pairing between an amiRNA and its targets is not the trigger for secondary small RNA production. In conjunction with the observation that amiRNAs with perfect complementarity to their target genes show high efficiency and specificity in gene silencing, this finding has an important bearing on future applications of amiRNAs in gene silencing in plants.