RNA N6-Methyladenosine Affects Copper-Induced Oxidative Stress Response in Arabidopsis thaliana.

Recently, post-transcriptional regulation of mRNA mediated by N6-methyladenosine (m6A) has been found to have profound effects on transcriptome regulation during plant responses to various abiotic stresses. However, whether this RNA modification can affect an oxidative stress response in plants has not been studied. To assess the role of m6A modifications during copper-induced oxidative stress responses, m6A-IP-seq was performed in Arabidopsis seedlings exposed to high levels of copper sulfate. This analysis revealed large-scale shifts in this modification on the transcripts most relevant for oxidative stress. This altered epitranscriptomic mark is known to influence transcript abundance and translation; therefore we scrutinized these possibilities. We found an increased abundance of copper-enriched m6A-containing transcripts. Similarly, we also found increased ribosome occupancy of copper-enriched m6A-containing transcripts, specifically those encoding proteins involved with stress responses relevant to oxidative stressors. Furthermore, the significance of the m6A epitranscriptome on plant oxidative stress tolerance was uncovered by assessing germination and seedling development of the mta (N6-methyladenosine RNA methyltransferase A mutant complemented with ABI3:MTA) mutant exposed to high copper treatment. These analyses suggested hypersensitivity of the mta mutant compared to the wild-type plants in response to copper-induced oxidative stress. Overall, our findings suggest an important role for m6A in the oxidative stress response of Arabidopsis.

Understanding plant stress memory traits can provide a way for sustainable agriculture.

Being sessile, plants encounter various biotic and abiotic threats in their life cycle. To minimize the damages caused by such threats, plants have acquired sophisticated response mechanisms. One major such response includes memorizing the encountered stimuli in the form of a metabolite, hormone, protein, or epigenetic marks. All of these individually as well as together, facilitate effective transcriptional and post-transcriptional responses upon encountering the stress episode for a second time during the life cycle and in some instances even in the future generations. This review attempts to highlight the recent advances in the area of plant memory. A detailed understanding of plant memory has the potential to offer solutions for developing climate-resilient crops for sustainable agriculture.

High Daytime Temperature Responsive MicroRNA Profiles in Developing Grains of Rice Varieties with Contrasting Chalkiness.

High temperature impairs starch biosynthesis in developing rice grains and thereby increases chalkiness, affecting the grain quality. Genome encoded microRNAs (miRNAs) fine-tune target transcript abundances in a spatio-temporal specific manner, and this mode of gene regulation is critical for a myriad of developmental processes as well as stress responses. However, the role of miRNAs in maintaining rice grain quality/chalkiness during high daytime temperature (HDT) stress is relatively unknown. To uncover the role of miRNAs in this process, we used five contrasting rice genotypes (low chalky lines Cyp, Ben, and KB and high chalky lines LaGrue and NB) and compared the miRNA profiles in the R6 stage caryopsis samples from plants subjected to prolonged HDT (from the onset of fertilization through R6 stage of caryopsis development). Our small RNA analysis has identified approximately 744 miRNAs that can be grouped into 291 families. Of these, 186 miRNAs belonging to 103 families are differentially regulated under HDT. Only two miRNAs, Osa-miR444f and Osa-miR1866-5p, were upregulated in all genotypes, implying that the regulations greatly varied between the genotypes. Furthermore, not even a single miRNA was commonly up/down regulated specifically in the three tolerant genotypes. However, three miRNAs (Osa-miR1866-3p, Osa-miR5150-3p and canH-miR9774a,b-3p) were commonly upregulated and onemiRNA (Osa-miR393b-5p) was commonly downregulated specifically in the sensitive genotypes (LaGrue and NB). These observations suggest that few similarities exist within the low chalky or high chalky genotypes, possibly due to high genetic variation. Among the five genotypes used, Cypress and LaGrue are genetically closely related, but exhibit contrasting chalkiness under HDT, and thus, a comparison between them is most relevant. This comparison revealed a general tendency for Cypress to display miRNA regulations that could decrease chalkiness under HDT compared with LaGrue. This study suggests that miRNAs could play an important role in maintaining grain quality in HDT-stressed rice.

A high-resolution single-molecule sequencing-based Arabidopsis transcriptome using novel methods of Iso-seq analysis.

BACKGROUND: Accurate and comprehensive annotation of transcript sequences is essential for transcript quantification and differential gene and transcript expression analysis. Single-molecule long-read sequencing technologies provide improved integrity of transcript structures including alternative splicing, and transcription start and polyadenylation sites. However, accuracy is significantly affected by sequencing errors, mRNA degradation, or incomplete cDNA synthesis. RESULTS: We present a new and comprehensive Arabidopsis thaliana Reference Transcript Dataset 3 (AtRTD3). AtRTD3 contains over 169,000 transcripts-twice that of the best current Arabidopsis transcriptome and including over 1500 novel genes. Seventy-eight percent of transcripts are from Iso-seq with accurately defined splice junctions and transcription start and end sites. We develop novel methods to determine splice junctions and transcription start and end sites accurately. Mismatch profiles around splice junctions provide a powerful feature to distinguish correct splice junctions and remove false splice junctions. Stratified approaches identify high-confidence transcription start and end sites and remove fragmentary transcripts due to degradation. AtRTD3 is a major improvement over existing transcriptomes as demonstrated by analysis of an Arabidopsis cold response RNA-seq time-series. AtRTD3 provides higher resolution of transcript expression profiling and identifies cold-induced differential transcription start and polyadenylation site usage. CONCLUSIONS: AtRTD3 is the most comprehensive Arabidopsis transcriptome currently. It improves the precision of differential gene and transcript expression, differential alternative splicing, and transcription start/end site usage analysis from RNA-seq data. The novel methods for identifying accurate splice junctions and transcription start/end sites are widely applicable and will improve single-molecule sequencing analysis from any species.

mRNA N6 -methyladenosine is critical for cold tolerance in Arabidopsis.

Plants respond to low temperatures by altering the mRNA abundance of thousands of genes contributing to numerous physiological and metabolic processes that allow them to adapt. At the post-transcriptional level, these cold stress-responsive transcripts undergo alternative splicing, microRNA-mediated regulation and alternative polyadenylation, amongst others. Recently, m6 A, m5 C and other mRNA modifications that can affect the regulation and stability of RNA were discovered, thus revealing another layer of post-transcriptional regulation that plays an important role in modulating gene expression. The importance of m6 A in plant growth and development has been appreciated, although its significance under stress conditions is still underexplored. To assess the role of m6 A modifications during cold stress responses, methylated RNA immunoprecipitation sequencing was performed in Arabidopsis seedlings esposed to low temperature stress (4°C) for 24 h. This transcriptome-wide m6 A analysis revealed large-scale shifts in this modification in response to low temperature stress. Because m6 A is known to affect transcript stability/degradation and translation, we investigated these possibilities. Interestingly, we found that cold-enriched m6 A-containing transcripts demonstrated the largest increases in transcript abundance coupled with increased ribosome occupancy under cold stress. The significance of the m6 A epitranscriptome on plant cold tolerance was further assessed using the mta mutant in which the major m6 A methyltransferase gene was mutated. Compared to the wild-type, along with the differences in CBFs and COR gene expression levels, the mta mutant exhibited hypersensitivity to cold treatment as determined by primary root growth, biomass, and reactive oxygen species accumulation. Furthermore, and most importantly, both non-acclimated and cold-acclimated mta mutant demonstrated hypersensitivity to freezing tolerance. Taken together, these findings suggest a critical role for the epitranscriptome in cold tolerance of Arabidopsis.

Molecular and Biochemical Analysis of Duplicated Cytosolic CuZn Superoxide Dismutases of Rice and in silico Analysis in Plants.

Superoxide dismutases (SODs, EC 1.15.1.1) are ubiquitous antioxidant metalloenzymes important for oxidative stress tolerance and cellular redox environment. Multiple factors have contributed toward the origin and diversity of SOD isoforms among different organisms. In plants, the genome duplication events, responsible for the generation of multiple gene copies/gene families, have also contributed toward the SOD diversity. However, the importance of such molecular events on the characteristics of SODs has not been studied well. This study investigated the effects of divergence on important characteristics of two block-duplicated rice cytosolic CuZn SODs (OsCSD1, OsCSD4), along with in silico assessment of similar events in other plants. The analysis revealed heterogeneity in gene length, regulatory regions, untranslated regions (UTRs), and coding regions of two OsCSDs. An inconsistency in the database-predicted OsCSD1 gene structure was also identified and validated experimentally. Transcript analysis showed differences in the basal levels and stress responsiveness of OsCSD1 and OsCSD4, and indicated the presence of two transcription start sites in the OsCSD1. At the amino acid level, the two OsCSDs showed differences at 18 sites; however, both exist as a homodimer, displaying typical CuZn SOD characteristics, and enhancing the oxidative stress tolerance of Escherichia coli cells. However, OsCSD4 showed higher specific activity as well as stability. The comparison of the two OsCSDs with reported thermostable CSDs from other plants identified regions likely to be associated with stability, while the homology modeling and superposition highlighted structural differences. The two OsCSDs displayed heteromeric interaction capability and forms an enzymatically active heterodimer (OsCSD1:OsCSD4) on co-expression, which may have significance as both are cytosolic. In silico analysis of 74 plant genomes revealed the prevalence of block duplications for multiple CSD copies (mostly cytosolic). The divergence and clustering analysis of CSDs suggested the possibility of an ancestral duplication event in monocots. Conserved SOD features indicating retention of SOD function among CSD duplicates were evident in few monocots and dicots. In most other species, the CSD copies lacked critical features and may not harbor SOD function; however, other feature-associated functions or novel functions might be present. These aspects of divergent CSD copies encoding co-localized CSDs may have implications in plant SOD functions in the cytosol and other organelles.

Molecular basis of priming-induced acquired tolerance to multiple abiotic stresses in plants.

The growth, survival, and productivity of plants are constantly challenged by diverse abiotic stresses. When plants are exposed to stress for the first time, they can capture molecular information and store it as a form of memory, which enables them to competently and rapidly respond to subsequent stress(es). This process is referred to as a priming-induced or acquired stress response. In this review, we discuss how (i) the storage and retrieval of the information from stress memory modulates plant physiological, cellular, and molecular processes in response to subsequent stress(es), (ii) the intensity, recurrence, and duration of priming stimuli influences the outcomes of the stress response, and (iii) the varying responses at different plant developmental stages. We highlight current understanding of the distinct and common molecular processes manifested at the epigenetic, (post-)transcriptional, and post-translational levels mediated by stress-associated molecules and metabolites, including phytohormones. We conclude by emphasizing how unravelling the molecular circuitry underlying diverse priming-stimuli-induced stress responses could propel the use of priming as a management practice for crop plants. This practice, in combination with precision agriculture, could aid in increasing yield quantity and quality to meet the rapidly rising demand for food.

The neglected other half - role of the pistil in plant heat stress responses.

Heat stress coinciding with reproductive stage leads to a significant loss in reproductive organs viability, resulting in lower seed-set and crop productivity. Successful fertilization and seed formation are determined by the viability of male and female reproductive organs. The impact of heat stress on the male reproductive organ (pollen) is studied more often compared to the female reproductive organ (pistil). This is attributed to easier accessibility of the pollen coupled with the notion that the pistil's role in fertilization and seed-set under heat stress is negligible. However, depending on species and developmental stages, recent studies reveal varying degrees of sensitivity of the pistil to heat stress. Remarkably, in some cases, the vulnerability of the pistil is even greater than the pollen. This article summarizes the current knowledge of the impact of heat stress on three critical stages of pistil for successful seed-set, that is, female reproductive organ development (gametogenesis), pollen-pistil interactions including pollen capture on stigma and pollen tube growth in style, as well as fertilization and early embryogenesis. Further, future research directions are suggested to unravel molecular basis of heat stress tolerance in pistil, which is critical for sustaining crop yields under predicted warming scenarios.

Splice Variants of Superoxide Dismutases in Rice and Their Expression Profiles under Abiotic Stresses.

The superoxide dismutases (SODs) play vital roles in controlling cellular reactive oxygen species (ROS) that are generated both under optimal as well as stress conditions in plants. The rice genome harbors seven SOD genes (CSD1, CSD2, CSD3, CSD4, FSD1, FSD2, and MSD) that encode seven constitutive transcripts. Of these, five (CSD2, CSD3, CSD4, FSD1, and MSD) utilizes an alternative splicing (AS) strategy and generate seven additional splice variants (SVs) or mRNA variants, i.e., three for CSD3, and one each for CSD2, CSD4, FSD1, and MSD. The exon-intron organization of these SVs revealed variations in the number and length of exons and/or untranslated regions (UTRs). We determined the expression patterns of SVs along with their constitutive forms of SODs in rice seedlings exposed to salt, osmotic, cold, heavy metal (Cu+2) stresses, as well as copper-deprivation. The results revealed that all seven SVs were transcriptionally active in both roots and shoots. When compared to their corresponding constitutive transcripts, the profiles of five SVs were almost similar, while two specific SVs (CSD3-SV4 and MSD-SV2) differed significantly, and the differences were also apparent between shoots and roots suggesting that the specific SVs are likely to play important roles in a tissue-specific and stress-specific manner. Overall, the present study has provided a comprehensive analysis of the SVs of SODs and their responses to stress conditions in shoots and roots of rice seedlings.

Characterizing microRNAs and their targets in different organs of Camellia sinensis var. assamica.

To comprehensively annotate miRNAs and their targets in tea plant, Camellia sinensis, we sequenced small and messenger RNAs of 9 samples of Camellia sinensis var. assamica (YK-10), a diploid elite cultivar widely grown in southwest China. In order to identify targets of miRNAs, we sequenced two degradome sequencing profiles from leaves and roots of YK-10, respectively. By analyzing the small RNA-Seq profiles, we newly identified 137 conserved miRNAs and 23 species specific miRNAs in the genome of YK-10, which significantly improved the annotation of miRNAs in tea plant. Approximately 2000 differently expressed genes were identified when comparing RNA-Seq profiles of any two of the three organs selected in the study. Totally, more than 5000 targets of conserved miRNAs were identified in the two degradome profiles. Furthermore, our results suggest that a few miRNAs play roles in the biosynthesis pathways of theanine, caffeine and flavonoid. These results enhance our understanding of small RNA guided gene regulations in different organs of tea plant.

Epigenetics and epigenomics: underlying mechanisms, relevance, and implications in crop improvement.

Epigenetics is defined as changes in gene expression that are not associated with changes in DNA sequence but due to the result of methylation of DNA and post-translational modifications to the histones. These epigenetic modifications are known to regulate gene expression by bringing changes in the chromatin state, which underlies plant development and shapes phenotypic plasticity in responses to the environment and internal cues. This review articulates the role of histone modifications and DNA methylation in modulating biotic and abiotic stresses, as well as crop improvement. It also highlights the possibility of engineering epigenomes and epigenome-based predictive models for improving agronomic traits.

Molecular Characterization of a Date Palm Vascular Highway 1-Interacting Kinase (PdVIK) Under Abiotic Stresses.

The date palm (Khalas) is an extremophile plant that can adapt to various abiotic stresses including drought and salinity. Salinity tolerance is a complex trait controlled by numerous genes. Identification and functional characterization of salt-responsive genes from the date palm is fundamental to understand salinity tolerance at the molecular level in this plant species. In this study, a salt-inducible vascular highway 1-interacting kinase (PdVIK) that is a MAP kinase kinase kinase (MAPKKK) gene from the date palm, was functionally characterized using in vitro and in vivo strategies. PdVIK, one of the 597 kinases encoded by the date palm genome possesses an ankyrin repeat domain and a kinase domain. The recombinant PdVIK protein exhibited phosphotyrosine activity against myelin basic protein (MBP) substrate. Overexpression of PdVIK in yeast significantly improved its tolerance to salinity, LiCl, and oxidative stresses. Transgenic Arabidopsis seedlings overexpressing PdVIK displayed improved tolerance to salinity, osmotic, and oxidative stresses as assessed by root growth assay. The transgenic lines grown in the soil also displayed modulated salt response, compared to wild-type controls as evaluated by the overall plant growth and proline levels. Likewise, the transgenic lines exhibited drought tolerance by maintaining better relative water content (RWC) compared to non-transgenic control plants. Collectively, these results implicate the involvement of PdVIK in modulating the abiotic stress response of the date palm.

The Effect of Drought on Transcriptome and Hormonal Profiles in Barley Genotypes With Contrasting Drought Tolerance.

Like many cereal crops, barley is also negatively affected by drought stress. However, due to its simple genome as well as enhanced stress resilient nature compared to rice and wheat, barley has been considered as a model to decipher drought tolerance in cereals. In the present study, transcriptomic and hormonal profiles along with several biochemical features were compared between drought-tolerant (Otis) and drought-sensitive (Baronesse) barley genotypes subjected to drought to identify molecular and biochemical differences between the genotypes. The drought-induced decrease in the leaf relative water content, net photosynthesis, and biomass accumulation was relatively low in Otis compared to Baronesse. The hormonal profiles did not reveal significant differences for majority of the compounds other than the GA20 and the cis-zeatin-o-glucoside (c-ZOG), whose levels were greatly increased in Otis compared to Baronesse under drought. The major differences that emerged from the transcriptome analysis are; (1), the overall number of differentially expressed genes was relatively low in drought-tolerant Otis compared to drought-sensitive Baronesse; (2), a wax biosynthesis gene (CER1), and NAC transcription factors were specifically induced in Otis but not in Baronesse; (3), the degree of upregulation of betaine aldehyde dehydrogenase and a homeobox transcription factor (genes with proven roles in imparting drought tolerance), was greater in Otis compared to Baronesse; (4) the extent of downregulation of gene expression profiles for proteins of the reaction center photosystem II (PSII) (D1 and D2) was low in Otis compared to Baronesse; and, (5), alternative splicing (AS) was also found to differ between the genotypes under drought. Taken together, the overall transcriptional responses were low in drought-tolerant Otis but the genes that could confer drought tolerance were either specifically induced or greatly upregulated in the tolerant genotype and these differences could be important for drought tolerance in barley.

An improved method of constructing degradome library suitable for sequencing using Illumina platform.

BACKGROUND: Post-transcriptional gene regulation is one of the critical layers of overall gene expression programs and microRNAs (miRNAs) play an indispensable role in this process by guiding cleavage on the messenger RNA targets. The transcriptome-wide cleavages on the target transcripts can be identified by analyzing the degradome or PARE or GMUCT libraries. However, high-throughput sequencing of PARE or degradome libraries using Illumina platform, a widely used platform, is not so straightforward. Moreover, the currently used degradome or PARE methods utilize MmeI restriction site in the 5' RNA adapter and the resulting fragments are only 20-nt long, which often poses difficulty in distinguishing between the members of the same target gene family or distinguishing miRNA biogenesis intermediates from the primary miRNA transcripts belonging to the same miRNA family. Consequently, developing a method which can generate longer fragments from the PARE or degradome libraries which can also be sequenced easily using Illumina platform is ideal. RESULTS: In this protocol, 3' end of the 5'RNA adaptor of TruSeq small RNA library is modified by introducing EcoP15I recognition site. Correspondingly, the double-strand DNA (dsDNA) adaptor sequence is also modified to suit with the ends generated by the restriction enzyme EcoP15I. These modifications allow amplification of the degradome library by primer pairs used for small RNA library preparation, thus amenable for sequencing using Illumina platform, like small RNA library. CONCLUSIONS: Degradome library generated using this improved protocol can be sequenced easily using Illumina platform, and the resulting tag length is ~ 27-nt, which is longer than the MmeI generated fragment (20-nt) that can facilitate better accuracy in validating target transcripts belonging to the same gene family or distinguishing miRNA biogenesis intermediates of the same miRNA family. Furthermore, this improved method allows pooling and sequencing degradome libraries and small RNA libraries simultaneously using Illumina platform.

Water-deficit responsive microRNAs in the primary root growth zone of maize.

BACKGROUND: MicroRNA-mediated gene regulatory networks play a significant role in plant growth and development and environmental stress responses. RESULTS: We identified 79 microRNAs (miRNAs) and multiple miRNA variants (isomiRs) belonging to 26 miRNA families in the primary root growth zone of maize seedlings grown at one of three water potentials: well-watered (- 0.02 MPa), mild water deficit stress (- 0.3 MPa), and severe water deficit stress (- 1.6 MPa). The abundances of 3 miRNAs (mild stress) and 34 miRNAs representing 17 families (severe stress) were significantly different in water-deficit stressed relative to well-watered controls (FDR < 0.05 and validated by stem loop RT-qPCR). Degradome sequencing revealed 213 miRNA-regulated transcripts and trancriptome profiling revealed that the abundance of 77 (miRNA-regulated) were regulated by water-defecit stress. miR399e,i,j-3p was strongly regulated by water-defcit stress implicating the possibility of nutrient deficiency during stress. CONCLUSIONS: We have identified a number of maize miRNAs that respond to specific water deficits applied to the primary root growth zone. We have also identified transcripts that are targets for miRNA regulation in the root growth zone under water-deficit stress. The miR399e,i,j-3p that is known to regulate phosphate uptake in response to nutrient deficiencies responds to water-deficit stress, however, at the seedling stage the seed provides adequate nutrients for root growth thus miR399e,i,j-3p may play a separate role in water-deficit responses. A water-deficit regulated maize transcript, similar to known miR399 target mimics, was identified and we hypothesized that it is another regulatory player, moderating the role of miR399e,i,j-3p, in primary root growth zone water deficit responses.

Metabolomic analysis of date palm seedlings exposed to salinity and silicon treatments.

Silicon is known to promote plant growth as well as stress tolerance of plants. The current study was undertaken to assess the growth promoting effect of silicon on date palm seedling development as well as its ability to abate some of the negative effects of salinity. In this study, date palm seedlings were treated with silicon and sodium chloride salts, and the effect of these salts on some physiological parameters of the plants was determined. In addition, a global nontargeted metabolomics analysis was performed for the leaf and root tissues using liquid chromatography-mass spectrometry (LC-MS). The results showed that under non-stress conditions, silicon treatment enhanced the growth of the date palm seedlings, however, under salinity, silicon slightly mitigates the negative effects of salt stress on the date palm seedlings although it enhances the potassium accumulation under this condition. The global metabolomics analysis has identified a total of 1,101 significant differentially accumulated (p, q ≤ 0.05) metabolites in leaves and roots under silicon, salinity or their combination. A differential pairwise metabolic profile comparison revealed the accumulation of distinct metabolites in response to silicon and salinity treatments such as antioxidant compounds pyridoxine, cepharanthine, allithiamine, myristic acid and boldine; osmoregulators such as mucic acid; along with the accumulation of detoxification intermediates such as S-D-lactoylglutathione, beta-cyano-L-alanine and gamma-glutamyl-conjugates. In addition, histochemical analyses revealed that application of silicon significantly (p ≤ 0.05) enhanced the formation of the Casparian strip. Identification of the differentially accumulated metabolites could offer an insight into how silicon is able to promote growth and salinity tolerance in date palms.

Overexpression of a Metallothionein 2A Gene from Date Palm Confers Abiotic Stress Tolerance to Yeast and Arabidopsis thaliana.

Although the date palm tree is an extremophile with tolerance to drought and certain levels of salinity, the damage caused by extreme salt concentrations in the soil, has created a need to explore stress-responsive traits and decode their mechanisms. Metallothioneins (MTs) are low-molecular-weight cysteine-rich proteins that are known to play a role in decreasing oxidative damage during abiotic stress conditions. Our previous study identified date palm metallothionein 2A (PdMT2A) as a salt-responsive gene, which has been functionally characterized in yeast and Arabidopsis in this study. The recombinant PdMT2A protein produced in Escherichia coli showed high reactivity against the substrate 5'-dithiobis-2-nitrobenzoic acid (DTNB), implying that the protein has the property of scavenging reactive oxygen species (ROS). Heterologous overexpression of PdMT2A in yeast (Saccharomyces cerevisiae) conferred tolerance to drought, salinity and oxidative stresses. The PdMT2A gene was also overexpressed in Arabidopsis, to assess its stress protective function in planta. Compared to the wild-type control, the transgenic plants accumulated less Na+ and maintained a high K+/Na+ ratio, which could be attributed to the regulatory role of the transgene on transporters such as HKT, as demonstrated by qPCR assay. In addition, transgenic lines exhibited higher chlorophyll content, higher superoxide dismutase (SOD) activity and improved scavenging ability for reactive oxygen species (ROS), coupled with a better survival rate during salt stress conditions. Similarly, the transgenic plants also displayed better drought and oxidative stress tolerance. Collectively, both in vitro and in planta studies revealed a role for PdMT2A in salt, drought, and oxidative stress tolerance.

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.

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.